Abstract

ABSTRACT We derive a consistent set of moment equations for cosmic ray (CR)-magnetohydrodynamics, assuming a gyrotropic distribution function (DF). Unlike previous efforts, we derive a closure, akin to the M1 closure in radiation hydrodynamics (RHD), that is valid in both the nearly isotropic DF and/or strong-scattering regimes, and the arbitrarily anisotropic DF or free-streaming regimes, as well as allowing for anisotropic scattering and transport/magnetic field structure. We present the appropriate two-moment closure and equations for various choices of evolved variables, including the CR phase space DF f, number density n, total energy e, kinetic energy ϵ, and their fluxes or higher moments, and the appropriate coupling terms to the gas. We show that this naturally includes and generalizes a variety of terms including convection/fluid motion, anisotropic CR pressure, streaming, diffusion, gyro-resonant/streaming losses, and re-acceleration. We discuss how this extends previous treatments of CR transport including diffusion and moment methods and popular forms of the Fokker–Planck equation, as well as how this differs from the analogous M1-RHD equations. We also present two different methods for incorporating a reduced speed of light (RSOL) to reduce time-step limitations: In both, we carefully address where the RSOL (versus true c) must appear for the correct behaviour to be recovered in all interesting limits, and show how current implementations of CRs with an RSOL neglect some additional terms.

Highlights

  • Cosmic rays (CRs) could play a potentially crucial role in the interstellar and circum-galactic medium, star and galaxy formation, and our understanding of high-energy astro-particle and plasma physics

  • These applications have generally relied on moment-based approaches, where one begins by assuming that the CR distribution function (DF) f is gyrotropic, averages over the micro-scale Lorentz forces and scattering processes, considers moments of the distribution function in terms of the remaining momentum direction, the pitch angle μ

  • Just like with radiation hydrodynamics (RHD), it is important that the reduced speed of light (RSOL) appear only in the dynamical equations for the CRs, not in the terms that couple to the gas that are written in terms of physical quantities

Read more

Summary

INTRODUCTION

Cosmic rays (CRs) could play a potentially crucial role in the interstellar and circum-galactic medium, star and galaxy formation, and our understanding of high-energy astro-particle and plasma physics. In planet/star/galaxy formation models the resolution scales are vastly larger than CR gyro-radii for CRs with energies TeV (which contain most of the energy/pressure, and dominate the interactions with the non-relativistic matter) As such, these applications have generally relied on moment-based approaches, where one begins by assuming that the CR distribution function (DF) f is gyrotropic (symmetric around the magnetic-field direction), averages over the micro-scale Lorentz forces and scattering processes, considers moments of the distribution function in terms of the remaining momentum direction, the pitch angle μ. These applications have generally relied on moment-based approaches, where one begins by assuming that the CR distribution function (DF) f is gyrotropic (symmetric around the magnetic-field direction), averages over the micro-scale Lorentz forces and scattering processes, considers moments of the distribution function in terms of the remaining momentum direction, the pitch angle μ The simplest of these – “zeroth moment methods” – correspond to pure diffusion models. We define the pitch-angle-averaging operations, pitch-angle moments of f , and DF-weighted pitch angle moments: X μ

DERIVATION OF THE CR TRANSPORT MOMENTS EQUATIONS
Scattering Terms
Focused Transport Equation to Leading Order
The Close-to-Isotropic-DF Case
The Maximally-Anisotropic-DF Case
Co-Moving Expressions to Leading Order
Spectrally-Integrated Expressions
DF Equation in Finite-Volume Form
Equations for the Mean Evolution of a CR “Group”
Propagation With Bent Fields in A Simple Geometry
Comparison to the M1 RHD Closure
Behavior of the CR Closures
Propagation in a Non-Trivial Field Geometry
Summary
Relation to Previous CR Moments Formulations
Relation to the Isotropic FP Equation
Relation to the M1 RHD Equations
Out of Equilibrium Behaviors and Timescales
Which Speed of Light Enters the Closure Relation?
Rigidity-Dependent RSOL
SUMMARY
Full Text
Published version (Free)

Talk to us

Join us for a 30 min session where you can share your feedback and ask us any queries you have

Schedule a call