Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userCosmic Ray Transport in the Ionized and Neutral ISM: MHD-PIC Simulations and Effective Fluid Treatments
103
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
Cosmic rays (CRs) have critical impacts in the multiphase interstellar medium (ISM), driving dynamical motions in low-density plasma and modifying the ionization state, temperature, and chemical composition of higher-density atomic and molecular gas. We present a study of CR propagation between the ionized ISM and a neutral cloud. Using one-dimensional magnetohydrodynamic particle-in-cell simulations which include ion-neutral drag to damp Alfv$acute{text{e}}$n waves in the cloud, we self-consistently evolve the kinetic physics of CRs and fluid dynamics of the multiphase gas. By introducing the cloud in our periodic domain, our simulations break translational symmetry and allow the emergence of spatial structure in the CR distribution function. A negative spatial gradient forms across the fully-ionized ISM region while a positive gradient forms across the neutral cloud. We connect our results with CR hydrodynamics formulations by computing the wave-particle scattering rates as predicted by quasilinear, fluid, and Fokker-Planck theory. For momenta where the mean free path is short relative to the box size, we find excellent agreement among all scattering rates. By exploring different cloud sizes and ion-neutral collision rates, we show that our results are robust. Our work provides a first-principles verification of CR hydrodynamics when particles stream down their pressure gradient, and opens a pathway toward comprehensive calibrations of transport coefficients from self-generated Alfv$acute{text{e}}$n wave scattering with CRs.
rate research
Read More
Recent in-situ and remote observations suggest that the transport regime associated with shock accelerated particles may be anomalous {i.e., the Mean Square Displacement (MSD) of such particles scales non-linearly with time}. We use self-consistent, hybrid PIC plasma simulations to simulate a quasi-parallel shock with parameters compatible with heliospheric shocks, and gain insights about the particle transport in such a system. For suprathermal particles interacting with the shock we compute the MSD separately in the upstream and downstream regions. Tracking suprathermal particles for sufficiently long times up and/or downstream of the shock poses problems in particle plasma simulations, such as statistically poor particle ensembles and trajectory fragments of variable length in time. Therefore, we introduce the use of time-averaged mean square displacement (TAMSD), which is based on single particle trajectories, as an additional technique to address the transport regime for the upstream and downstream regions. MSD and TAMSD are in agreement for the upstream energetic particle population, and both give a strong indication of superdiffusive transport, consistent with interplanetary shock observations. MSD and TAMSD are also in reasonable agreement downstream, where indications of anomalous transport are also found. TAMSD shows evidence of heterogeneity in the diffusion properties of the downstream particle population, ranging from subdiffusive behaviour of particles trapped in the strong magnetic field fluctuations generated at the shock, to superdiffusive behaviour of particles transmitted and moving away from the shock.
Particle acceleration in collisionless plasma systems is a central question in astroplasma and astroparticle physics. The structure of the acceleration regions, electron-ion energy equilibration, preacceleration of particles at shocks to permit further energization by diffusive shock acceleration, require knowledge of the distribution function of particles besides the structure and dynamic of electromagnetic fields, and hence a kinetic description is desirable. Particle-in-cell simulations offer an appropriate, if computationally expensive method of essentially conducting numerical experiments that explore kinetic phenomena in collisionless plasma. We review recent results of PIC simulations of astrophysical plasma systems, particle acceleration, and the instabilities that shape them.
We present the second data release for the HI-MaNGA programme of HI follow-up observations for the SDSS-IV MaNGA survey. This release contains measurements for 3669 unique galaxies, combining 2108 Green Bank Telescope observations with an updated crossmatch of the MaNGA sample with the ALFALFA survey. We combine these data with MaNGA spectroscopic measurements to examine relationships between HI-to-stellar mass ratio (M_HI/M_*) and average ISM/star formation properties probed by optical emission lines. M_HI/M_* is very weakly correlated with the equivalent width of Halpha, implying a loose connection between the instantaneous star formation rate and the HI reservoir, although the link between M_HI/M_* and star formation strengthens when averaged even over only moderate timescales (~30 Myrs). Galaxies with elevated HI depletion times have enhanced [OI]/Halpha and depressed Halpha surface brightness, consistent with more HI residing in a diffuse and/or shock heated phase which is less capable of condensing into molecular clouds. Of all optical lines, M_HI/M_* correlates most strongly with oxygen equivalent width, EW(O), which is likely a result of the existing correlation between M_HI/M_* and gas-phase metallicity. Residuals in the M_HI/M_*-EW(O) relation are again correlated with [OI]/Halpha and Halpha surface brightness, suggesting they are also driven by variations in the fraction of diffuse and/or shock-heated gas. We recover the strong anti-correlation between M_HI/M_* and gas-phase metallicity seen in previous studies. We also find a relationship between M_HI/M_* and [OI]/Halpha, suggesting that higher fractions of diffuse and/or shock-heated gas are more prevalent in gas-rich galaxies.
The transport of charged energetic particles in the presence of strong intermittent heliospheric turbulence is computationally analyzed based on known properties of the interplanetary magnetic field and solar wind plasma at 1 Astronomical Unit (AU). The turbulence is assumed to be static, composite, and quasi-three-dimensional with a varying energy distribution between a one-dimensional Alfvenic (slab) and a structured two-dimensional component. The spatial fluctuations of the turbulent magnetic field are modeled either as homogeneous with a Gaussian probability distribution function (PDF), or as intermittent on large and small scales with a q-Gaussian PDF. Simulations showed that energetic particle diffusion coefficients both parallel and perpendicular to the background magnetic field are significantly affected by intermittency in the turbulence. This effect is especially strong for parallel transport where for large-scale intermittency results show an extended phase of subdiffusive parallel transport during which cross-field transport diffusion dominates. The effects of intermittency are found to depend on particle rigidity and the fraction of slab energy in the turbulence, yielding a perpendicular to parallel mean free path ratio close to 1 for large-scale intermittency. Investigation of higher order transport moments (kurtosis) indicates that non-Gaussian statistical properties of the intermittent turbulent magnetic field are present in the parallel transport, especially for low rigidity particles at all times.
Astrophysical plasmas are turbulent and magnetized. The interaction between cosmic rays (CRs) and magnetohydrodynamic (MHD) turbulence is a fundamental astrophysical process. Based on the current understanding of MHD turbulence, we revisit the trapping of CRs by magnetic mirrors in the context of MHD turbulence. In compressible MHD turbulence, isotropic fast modes dominate both trapping and gyroresonant scattering of CRs. The presence of trapping significantly suppresses the pitch-angle scattering and the spatial diffusion of CRs along the magnetic field. The resulting parallel diffusion coefficient has a weaker dependence on CR energy at higher energies. In incompressible MHD turbulence, the trapping by pseudo-Alfv{e}n modes dominates over the gyroresonant scattering by anisotropic Alfv{e}n and pseudo-Alfv{e}n modes at all pitch angles and prevents CRs from diffusion.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
