Kinetic simulations of the cosmic ray pressure anisotropy instability: cosmic ray scattering rate in the saturated state

Xiaochen Sun,Xue-Ning Bai,Xihui Zhao
2024-09-13
Abstract:Cosmic ray (CR) feedback plays a vital role in shaping the formation and evolution of galaxies through their interaction with magnetohydrodynamic waves. In the CR self-confinement scenario, the waves are generated by the CR gyro-resonant instabilities via CR streaming or CR pressure anisotropy, and saturate by balancing wave damping. The resulting effective particle scattering rate by the waves, {\nu}eff, critically sets the coupling between the CRs and background gas, but the efficiency of CR feedback is yet poorly constrained. We employ 1D kinetic simulations under the Magnetohydrodynamic-Particle-In-Cell (MHD-PIC) framework with the adaptive {\delta}f method to quantify {\nu}eff for the saturated state of the CR pressure anisotropy instability (CRPAI) with ion-neutral friction. We drive CR pressure anisotropy by expanding/compressing box, mimicking background evolution of magnetic field strength, and the CR pressure anisotropy eventually reaches a quasi-steady state by balancing quasi-linear diffusion. At the saturated state, we measure {\nu}eff and the CR pressure anisotropy level, establishing a calibrated scaling relation with environmental parameters. The scaling relation is consistent with quasi-linear theory and can be incorporated to CR fluid models, in either the single-fluid or p-by-p treatments. Our results serve as a basis towards accurately calibrating the subgrid physics in macroscopic studies of CR feedback and transport.
High Energy Astrophysical Phenomena,Astrophysics of Galaxies
What problem does this paper attempt to address?
The problem that this paper attempts to solve is the effective scattering rate ($\nu_{\text{eff}}$) of cosmic rays and their coupling efficiency to the background gas after the wave growth and damping mechanisms are balanced by ion - neutral friction in the saturated state of cosmic - ray pressure - anisotropy instabilities (CRPAI). Specifically, the paper aims to quantify the level of cosmic - ray pressure anisotropy and the effective scattering rate in the saturated state of CRPAI through one - dimensional kinetic simulations (using the magnetohydrodynamic - particle - in - cell method, i.e., the MHD - PIC framework), and establish a calibration relationship between these parameters and environmental conditions. This helps to more accurately understand the micro - physical processes in the cosmic - ray feedback mechanism, thereby providing a more accurate sub - grid physical model for galactic simulations on the macro - scale. ### Main research objectives: 1. **Quantify the effective scattering rate of cosmic rays**: Measure the effective scattering rate $\nu_{\text{eff}}$ of cosmic rays by simulating the saturated state of CRPAI. 2. **Establish the calibration relationship of environmental parameters**: Determine the relationship between the level of cosmic - ray pressure anisotropy and environmental parameters (such as magnetic field strength, gas density, etc.). 3. **Verify theoretical predictions**: Compare the simulation results with the predictions of the quasi - linear theory (QLT) to calibrate and verify the theoretical model. ### Research background: - **Cosmic - ray feedback**: Cosmic rays play an important role in the formation and evolution of galaxies through their interaction with magnetohydrodynamic waves. - **Self - limiting mechanism**: In the CR self - limiting scenario, waves are generated by the gyroresonance instability caused by CR flow or CR pressure anisotropy and reach equilibrium through wave damping. - **Importance of scattering rate**: The effective particle scattering rate $\nu_{\text{eff}}$ determines the coupling degree between CR and the background gas, but this efficiency is currently unclear. ### Methods and techniques: - **MHD - PIC method**: Use the adaptive $\delta f$ method to reduce noise and simulate the saturated state of CRPAI. - **Expansion / compression box framework**: Drive CR pressure anisotropy by expanding or compressing the simulation box to simulate changes in the background magnetic field strength. - **Quasi - linear theory**: Predict the CR scattering rate and pressure anisotropy level based on QLT and compare them with the simulation results. ### Simulation settings: - **Initial conditions**: CRs are initialized with an isotropic $\kappa$ - distribution, and the background gas has a uniform density and magnetic field. - **Parameter selection**: Select reasonable physical parameters to ensure that the simulation results can be generalized to practical applications. - **Numerical method**: Use the ATHENA++ code for two - stage van Leer time integration, and the Roe solver and Boris pusher are used for the integration of MHD equations and CR particles respectively. ### Expected results: - **Quantitative measurement**: Provide quantitative measurements of the effective scattering rate and pressure anisotropy level of cosmic rays in the saturated state of CRPAI. - **Calibration relationship**: Establish a calibration relationship between these parameters and environmental conditions to provide a basis for macro - simulations. - **Theoretical verification**: Verify the predictions of the quasi - linear theory and improve the understanding of the cosmic - ray feedback mechanism. Through these studies, the paper aims to fill the current gap in the understanding of the micro - physical processes of CRPAI and provide a more accurate physical model for future galactic simulations.