Cosmic Ray Transport Research Articles

Gas and Dust Temperature in Prestellar Cores Revisited: New Limits on Cosmic-Ray Ionization Rate

Abstract We develop a self-consistent model for the equilibrium gas temperature and size-dependent dust temperature in cold, dense, prestellar cores, assuming an arbitrary power-law size distribution of dust grains. Compact analytical expressions applicable to a broad range of physical parameters are derived and compared with predictions of the commonly used standard model. It is suggested that combining the theoretical results with observations should allow us to constrain the degree of dust evolution and the cosmic-ray ionization rate in dense cores, and to help with discriminating between different regimes of cosmic-ray transport in molecular clouds. In particular, assuming a canonical MRN distribution of grain sizes, our theory demonstrates that the gas-temperature measurements in the prestellar core L1544 are consistent with an ionization rate as high as ∼10−16 s−1, an order of magnitude higher than previously thought.

An Improved Treatment of Neutral Sheet Drift in the Inner Heliosphere

Abstract Drifts due to the curvature and gradients in the heliospheric magnetic field, as well as along the heliospheric current sheet, have long been known to play a significant role in the transport of galactic cosmic rays. Recently, there has been greater interest in the role these drifts play in the transport of solar energetic particles. This study proposes an approach to modeling particle drift velocities in particle transport codes that, while being relatively simple to implement and computationally inexpensive, also models drift effects accurately across a broader range of energies than previous approaches.

Correlation of Long-term Cosmic-Ray Modulation with Solar Activity Parameters

Abstract In this work, we analyze the long-term cosmic-ray modulation observed by the Hermanus neutron monitor, which is the detector with the longest cosmic-ray record, from 1957 July. For our study we use the force-field approximation to the cosmic-ray transport equation, and the newest results on the mean free paths from the scattering theory. We compare the modulation parameter (ϕ) with different rigidity (P) dependences: P, P 2, and P 2/3. We correlate them with solar and interplanetary parameters. We found that (1) these rigidity dependences properly describe the modulation, (2) long-term cosmic-ray variations are better correlated with the magnitude of the heliospheric magnetic field (HMF) than the sunspot number, solar wind speed, and tilt angle of the HMF, and (3) the theoretical dependence of the parallel mean free path on the magnetic field variance is in agreement with the modulation parameter and therefore with the neutron monitor record. We also found that the force-field approximation is not able to take into account the effects of three-dimensional particle transport, showing a poor correlation with the perpendicular mean free path.

Role of cosmic-ray streaming and turbulent damping in driving galactic winds

ABSTRACT Large-scale galactic winds driven by stellar feedback are one phenomenon that influences the dynamical and chemical evolution of a galaxy, redistributing material throughout the circumgalatic medium. Non-thermal feedback from galactic cosmic rays (CRs) – high-energy charged particles accelerated in supernovae and young stars – can impact the efficiency of wind driving. The streaming instability limits the speed at which they can escape. However, in the presence of turbulence, the streaming instability is subject to suppression that depends on the magnetization of turbulence given by its Alfvén Mach number. While previous simulations that relied on a simplified model of CR transport have shown that super-Alfvénic streaming of CRs enhances galactic winds, in this paper we take into account a realistic model of streaming suppression. We perform three-dimensional magnetohydrodynamic simulations of a section of a galactic disc and find that turbulent damping dependent on local magnetization of turbulent interstellar medium (ISM) leads to more spatially extended gas and CR distributions compared to the earlier streaming calculations, and that scale heights of these distributions increase for stronger turbulence. Our results indicate that the star formation rate increases with the level of turbulence in the ISM. We also find that the instantaneous wind mass loading is sensitive to local streaming physics with the mass loading dropping significantly as the strength of turbulence increases.

On the Growth and Saturation of the Gyroresonant Streaming Instabilities

Abstract The self-regulation of cosmic-ray (CR) transport in the interstellar and intracluster media has long been viewed through the lenses of linear and quasi-linear kinetic plasma physics. Such theories are believed to capture the essence of CR behavior in the presence of self-generated turbulence but cannot describe potentially critical details arising from the nonlinearities of the problem. We utilize the particle-in-cell numerical method to study the time-dependent nonlinear behavior of the gyroresonant streaming instabilities, self-consistently following the combined evolution of particle distributions and self-generated wave spectra in one-dimensional periodic simulations. We demonstrate that the early growth of instability conforms to the predictions from linear physics, but that the late-time behavior can vary depending on the properties of the initial CR distribution. We emphasize that the nonlinear stages of instability depend strongly on the initial anisotropy of CRs—highly anisotropic CR distributions do not efficiently reduce to Alfvénic drift velocities, owing to reduced production of left-handed resonant modes. We derive estimates for the wave amplitudes at saturation and the timescales for nonlinear relaxation of the CR distribution and then demonstrate the applicability of these estimates to our simulations. Bulk flows of the background plasma due to the presence of resonant waves are observed in our simulations, confirming the microphysical basis of CR-driven winds.

Spectrum of Electrons Accelerated in Supernova Remnants.

Using a semianalytic model of nonlinear diffusive shock acceleration, we model the total spectrum of cosmic ray (CR) electrons accelerated by supernova remnants (SNRs). Because electrons experience synchrotron losses in the amplified magnetic fields characteristic of SNRs, they exhibit substantially steeper spectra than protons. In particular, we find that the difference between the electron and proton spectral index (power law slope) ranges from 0.1 to 0.4. Our findings must be reckoned with theories of Galactic CR transport, which often assume that electrons and protons are injected with the same slope, and may especially have implications for the observed "positron excess".

Breaks in interstellar spectra of positrons and electrons derived from time-dependent AMS data

Until fairly recently, it was widely accepted that local cosmic ray spectra were largely featureless power laws, containing limited information on their acceleration and transport. This viewpoint is currently being revised in the light of evidence for a variety of spectral breaks in the fluxes of cosmic ray nuclei. Here, we focus on cosmic ray electrons and positrons which at the highest energies must be of local origin due to strong radiative losses. We consider a pure diffusion model for their Galactic transport and determine its free parameters by fitting data in a wide energy range: measurements of the interstellar spectrum by Voyager at MeV energies, radio synchrotron data (sensitive to GeV electrons and positrons) and local observations by AMS up to ~ 1 TeV. For the first time, we also model the time-dependent fluxes of cosmic ray electrons and positrons at GeV energies recently presented by AMS, treating solar modulation in a simple extension of the widely used force-field approximation. We are able to reproduce all the available measurements to date. Our model of the interstellar spectrum of cosmic ray electrons and positrons requires the presence of a number of spectral breaks, both in the source spectra and the diffusion coefficients. While we remain agnostic as to the origin of these spectral breaks, their presence will inform future models of the microphysics of cosmic ray acceleration and transport.

Warped diffusive radio halo around the quiescent spiral edge-on galaxy NGC 4565

Context. Cosmic rays play a pivotal role in launching galactic winds, particularly in quiescently star-forming galaxies where the hot gas alone is not sufficient to drive a wind. Except for the Milky Way, not much is known about the transport of cosmic rays in galaxies. Aims. In this Letter, we present low-frequency observations of the nearby edge-on spiral galaxy NGC 4565 using the LOw-Frequency ARray (LOFAR). With our deep 144 MHz observations, we obtain a clean estimate of the emission originating from old cosmic-ray electrons (CRe), which is almost free from contamination by thermal emission. Methods. We measured vertical profiles of the non-thermal radio continuum emission that we fitted with Gaussian and exponential functions. The different profile shapes correspond to 1D cosmic-ray transport models of pure diffusion and advection, respectively. Results. We detect a warp in the radio continuum that is reminiscent of the previously known H I warp. Because the warp is not seen at GHz-frequencies in the radio continuum, its minimum age must be about 100 Myr. The warp also explains the slight flaring of the thick radio disc that can otherwise be well described by a Gaussian profile with an FWHM of 65 arcsec (3.7 kpc). Conclusions. The diffusive radio halo together with the extra-planar X-ray emission may be remnants of enhanced star-forming activity in the past where the galaxy had a galactic wind, as GHz-observations indicate only a weak outflow in the last 40 Myr. NGC 4565 could be in transition from an outflow- to an inflow-dominated phase.

Effects of re-acceleration and source grammage on secondary cosmic rays spectra

ABSTRACT The ratio between secondary and primary cosmic ray (CR) particles is the main source of information about CR propagation in the Galaxy. Primary CRs are thought to be accelerated mainly in supernova remnant shocks and then released in the interstellar medium. Here, they produce secondary particles by occasional collisions with interstellar matter. As a result, the ratio between the fluxes of secondary and primary particles carries information about the amount of matter CRs have encountered during their journey from their sources to the Earth. Recent measurements by AMS-02 revealed an unexpected behaviour of two main secondary-to-primary ratios, the Boron-to-Carbon ratio and the antiproton-to-proton ratio. In this work, we discuss how such anomalies may reflect the action of two phenomena that are usually overlooked, namely the fact that some fraction of secondary particles can be produced within the acceleration region, and the non-negligible probability that secondary particles encounter an accelerator (and are re-accelerated) during propagation. Both effects must be taken into account in order to correctly extract information about CR transport from secondary-to-primary ratios.

Cosmic-ray transport from AMS-02 boron to carbon ratio data: Benchmark models and interpretation

This article aims at establishing new benchmark scenarios for Galactic cosmic-ray propagation in the GV-TV rigidity range, based on fits to the AMS-02 B/C data with the USINE v3.5 propagation code. We employ a new fitting procedure, cautiously taking into account data systematic error correlations in different rigidity bins and considering Solar modulation potential and leading nuclear cross-section as nuisance parameters. We delineate specific low, intermediate, and high-rigidity ranges that can be related to both features in the data and peculiar microphysics mechanisms resulting in spectral breaks. We single out a scenario which yields excellent fits to the data and includes all the presumably relevant complexity, the BIG model. This model has two limiting regimes: (i) the SLIM model, a minimal diffusion-only setup, and (ii) the QUAINT model, a convection-reacceleration model where transport is tuned by non-relativistic effects. All models lead to robust predictions in the high-energy regime ($\gtrsim10$GV), i.e. independent of the propagation scenario: at $1\sigma$, the diffusion slope $\delta$ is $[0.43-0.53]$, whereas $K_{10}$, the diffusion coefficient at 10GV, is $[0.26-0.36]$kpc$^2$Myr$^{-1}$; we confirm the robustness of the high-energy break, with a typical value $\Delta_h\sim 0.2$. We also find a hint for a similar (reversed) feature at low rigidity around the B/C peak ($\sim 4$GV) which might be related to some effective damping scale in the magnetic turbulence.

Numerical modeling of cosmic-ray transport in the heliosphere and interpretation of the proton-to-helium ratio in Solar Cycle 24

Thanks to space-borne experiments of cosmic-ray (CR) detection, such as the AMS and PAMELA missions in low-Earth orbit, or the Voyager-1 spacecraft in the interstellar space, a large collection of multi-channel and time-resolved CR data has become available. Recently, the AMS experiment has released new precision data, on the proton and helium fluxes in CRs, measured on monthly basis during its first six years of mission. The AMS data reveal a remarkable long-term behavior in the temporal evolution of the proton-to-helium ratio at rigidity R≡p/Z≲,3 GV. As we have argued in a recent work, such a behavior may reflect the transport properties of low-rigidity CRs in the inteplanetary space. In particular, it can be caused by mass/charge dependence of the CR diffusion coefficient. In this paper, we present our developments in the numerical modeling of CR transport in the Milky Way and in the heliosphere. Within our model, and with the help of approximated analytical solutions, we describe in details the relations between the properties of CR diffusion and the time-dependent evolution of the proton-to-helium ratio.

Constraints on cosmic-ray transport in galaxy clusters from radio and γ-ray observations

ABSTRACT The nature of cosmic rays (CRs) and their transport in galaxy clusters is probed by several observations. Radio observations reveal synchrotron radiation of cosmic-ray electrons (CRe) spiralling around cluster magnetic fields. γ-ray observations reveal hadronic reactions of cosmic-ray protons (CRp) with gas nuclei that produce pions. No such cluster-wide γ-ray signal has been measured, putting an upper limit on the density of CRp in clusters. But the presence of CRe implies some source of CRp, and consequently there must be some CRp-loss mechanism. We quantify the observational constraints on this mechanism assuming that losses are dominated by CR transport, ultimately deriving lower limits on this transport. Using the Coma cluster as an example, we find that bulk outward speeds of 10–100 km s−1 are sufficient to reduce γ-radiation below current upper limits. These speeds are sub-Alfvénic and are consistent with a self-confinement model for CR transport if the magnetic field is coherent on large scales. If the transport is diffusive, we require minimum diffusion coefficients of 1031–1032 cm2 s−1. This is consistent with CRs free streaming at the speed of light along a field tangled on length-scales of a few kpc. We find that a model of the Coma cluster with a tangled field and the self-confinement picture together can be consistent with observations if the relative acceleration efficiency of CR protons is less than 15 times more than that of electrons of the same energy. This value is 3–6 times lower than the same quantity for Galactic cosmic rays.

Fickian and non-Fickian diffusion of cosmic rays

Abstract Fluid approximations to cosmic ray (CR) transport are often preferred to kinetic descriptions in studies of the dynamics of the interstellar medium (ISM) of galaxies, because they allow simpler analytical and numerical treatments. Magnetohydrodynamic simulations of the ISM usually incorporate CR dynamics as an advection–diffusion equation for CR energy density, with anisotropic, magnetic-field-aligned diffusion with the diffusive flux assumed to obey Fick’s law. We compare test particle and fluid simulations of CRs in a random magnetic field. We demonstrate that a non-Fickian prescription of CR diffusion, which corresponds to the telegraph equation for the CR energy density, can be easily calibrated to match the test particle simulations with great accuracy. In particular, we consider a random magnetic field in the fluid simulation that has a lower spatial resolution than that used in the particle simulation to demonstrate that an appropriate choice of the diffusion tensor can account effectively for the unresolved (subgrid) scales of the magnetic field. We show that the characteristic time that appears in the telegraph equation can be physically interpreted as the time required for the particles to reach a diffusive regime and we stress that the Fickian description of the CR fluid is unable to describe complex boundary or initial conditions for the CR energy flux.

13th International Conference on Numerical Modeling of Space Plasma Flows (ASTRONUM-2018)

PREFACEThe Center for Space Plasma and Aeronomic Research (CSPAR) at the University of Alabama in Huntsville (UAH) and Maison de la Simulation, a joint laboratory at the French Alternative Energies and Atomic Energy Commission (CEA), National Center for Scientific Research (CNRS), University of Paris-Sud, and University of Versailles, organized the 13th International Conference on Numerical Modeling of Space Plasma Flows (ASTRONUM-2018) on June 25—29, 2018 in Panama City Beach, Florida, USA.The Program Committee consisted of Tahar Amari (CNRS Ecole Polytechnique, France), Edouard Audit (CEA/CNRS Maison de la Simulation, Gif-sur-Yvette, France, co-chair), Amitava Bhattacharjee (Princeton University, USA), Phillip Colella (Lawrence Berkeley National Laboratory, USA), Anthony Mezzacappa (University of Tennessee, Knoxville, USA), Ewald Müller (Max-Planck-Institute for Astrophysics, Garching, Germany), Nikolai Pogorelov (University of Alabama in Huntsville/CSPAR, USA, chair), Kazunari Shibata (Kyoto University, Japan), James Stone (Princeton University, USA), Jon Linker (Predictive Science, Inc., USA), and Gary P. Zank (University of Alabama in Huntsville, USA).The conference attracted 67 scientists from 14 countries representing different branches of the plasma simulation community. The distinctive feature of this conference is a combination of diverse research topics, all of which are essential for performing high-resolution, continuum mechanics and particle, simulations of physical phenomena in space physics and astrophysics. Among such topics were software packages for modeling and analyzing plasma flows; advanced numerical methods for space and astrophysical flows; large-scale fluid-based, kinetic, and hybrid simulations; turbulence and cosmic ray transport; and magnetohydrodynamics. The applications discussed included cosmology and galaxy formation, supernova explosions, physics of the Sun-heliosphere-magnetosphere environments, theinterstellar medium and star formation, stellar physics, experimental plasma physics, astrophysical accretion, numerical methods for ideal and non-ideal, relativistic and nonrelativistic MHD, etc. The proceedings volume is structured so that it covers all of these topics.

Magnetohydrodynamic Particle-in-cell Simulations of the Cosmic-Ray Streaming Instability: Linear Growth and Quasi-linear Evolution

Abstract The gyro-resonant cosmic-ray (CR) streaming instability is believed to play a crucial role in CR transport, leading to the growth of Alfvén waves at small scales that scatter CRs, and impacts the interaction of CRs with the interstellar medium (ISM) on large scales. However, extreme scale separation (λ ≪ pc), low CR number density (n CR/n ISM ∼ 10−9), and weak CR anisotropy (∼v A/c) pose strong challenges for proper numerical studies of this instability on a microphysical level. Employing the recently developed magnetohydrodynamic particle-in-cell method, which has unique advantages to alleviate these issues, we conduct 1D simulations that quantitatively demonstrate the growth and saturation of the instability in the parameter regime consistent with realistic CR streaming in the large-scale ISM. Our implementation of the δf method dramatically reduces Poisson noise and enables us to accurately capture wave growth over a broad spectrum equally shared between left- and right-handed Alfvén modes. We are also able to accurately follow the quasi-linear diffusion of CRs subsequent to wave growth, which is achieved by employing phase randomization across periodic boundaries. Full isotropization of the CRs in the wave frame requires the pitch angles of most CRs to efficiently cross 90° and can be captured in simulations with relatively high wave amplitude and/or spatial resolution. We attribute this crossing to nonlinear wave–particle interaction (rather than mirror reflection) by investigating individual CR trajectories. We anticipate that our methodology will open up opportunities for future investigations that incorporate additional physics.

Ultrahigh energy cosmic rays and neutrinos from light nuclei composition

The baryonic mass composition of ultrahigh energy ($\gtrsim 10^{18}$ eV) cosmic rays (UHECRs) at injection accompanied by their interactions on universal photon backgrounds during propagation directly governs the UHECR flux on the Earth. Secondary neutrinos and photons produced in these interactions serve as crucial astrophysical messengers of UHECR sources. A modeling of the latest data obtained by the Pierre Auger Observatory (PAO) suggests a mixed element composition of UHECRs with the sub-ankle spectrum being explained by a different class of sources than the super-ankle region ($> 10^{18.7}$ eV). In this work, we obtain two kinds of fit to the UHECR spectrum -- one with a single population of sources comprising of $^1$H and $^2$He, over an energy range commencing at $\approx 10^{18}$ eV -- another for a mixed composition of representative nuclei $^1$H, $^4$He, $^{14}$N and $^{28}$Si at injection, for which a fit is obtained from above $\approx 10^{18.7}$ eV. In both cases, we consider the source emissivity evolution to be a simple power-law in redshift. We test the credibility of H+He composition by varying the source properties over a wide range of values and compare the results to that obtained for H+He+N+Si composition, using the Monte Carlo simulation tool CRPropa 3. The secondary electrons and photons are propagated using the cosmic ray transport code DINT. We place limits on the source spectral index, source evolution index and cutoff rigidity of the source population in each case by fitting the UHECR spectrum. Cosmogenic neutrino fluxes can further constrain the abundance fraction and maximum source redshift in case of light nuclei injection model.

Cosmic ray transport and radiative processes in nuclei of starburst galaxies

AbstractThe high rate of star formation and supernova explosions of starburst galaxies make them interesting sources of high-energy radiation. Depending on the level of turbulence present in their interstellar medium, the bulk of cosmic rays produced inside starburst galaxies may lose most of their energy before escaping, thereby making these sources behave as calorimeters, at least up to some maximum energy. Contrary to previous studies, here we investigate in detail the conditions under which cosmic ray confinement may be effective for electrons and nuclei and we study the implications of cosmic ray confinement in terms of multifrequency emission from starburst nuclei and production of high-energy neutrinos. The general predictions are then specialized to three cases of active starbursts, namely, M82, NGC 253, and Arp220. Both primary and secondary electrons, as well as electron–positron pairs produced by gamma-ray absorption inside starburst galaxies are taken into account. Electrons and positrons produced as secondary products of hadronic interactions are found to be responsible for most of the emission of leptonic origin. In particular, synchrotron emission of very high energy secondary electrons produces an extended emission of hard X-rays that represents a very interesting signature of hadronic process in starburst galaxies, potentially accessible to current and future observations in the X-ray band. A careful understanding of both the production and absorption of gamma-rays in starburst galaxies is instrumental to the assessment of the role of these astrophysical sources as sources of high-energy astrophysical neutrinos.

The HelMod model in the works for inner and outer heliosphere: From AMS to Voyager probes observations

HelMod is a Monte Carlo code developed to describe the transport of Galactic Cosmic Rays (GCRs) through the heliosphere from the interstellar space to the Earth. In the current HelMod version 4 the modulation process, based on Parker’s equation, is applied to the propagation of GCRs in the inner and outer heliosphere, i.e., including the heliosheath. HelMod was proved to reproduce protons, nuclei and electrons cosmic rays spectra observed during solar cycles 23–24 by several detectors, for instance, PAMELA, BESS and AMS-02. In particular, the unprecedented accuracy of AMS-02 observations allowed one a better tuning of the description regarding the solar modulation mechanisms implemented in HelMod. In addition, HelMod demonstrated to be capable of reproducing the fluxes observed by the Voyager probes in the inner and outer regions of heliosphere up to its border.

CHANG-ES

Context. The observation of total and linearly polarized synchrotron radiation of spiral galaxies in the radio continuum reveals the distribution and structure of their magnetic fields. By observing these, information about the proposed dynamo processes that preserve the large-scale magnetic fields in spiral galaxies can be gained. Additionally, by analyzing the synchrotron intensity, the transport processes of cosmic rays into the halo of edge-on spiral galaxies can be investigated. Aims. We analyze the magnetic field geometry and the transport processes of the cosmic rays of the edge-on spiral starburst galaxy NGC 4666 from CHANG-ES radio data in two frequencies; 6 GHz (C-band) and 1.5 GHz (L-band). Supplementary X-ray data are used to investigate the hot gas in NGC 4666. Methods. We determine the radio scale heights of total power emission at both frequencies for this galaxy. We show the magnetic field orientations derived from the polarization data. Using rotation measure (RM) synthesis we further study the behavior of the RM values along the disk in C-band to investigate the large-scale magnetic-field pattern. We use the revised equipartition formula to calculate a map of the magnetic field strength. Furthermore, we model the processes of cosmic-ray transport into the halo with the 1D SPINNAKER model. Results. The extended radio halo of NGC 4666 is box-shaped and is probably produced by the previously observed supernova-driven superwind. This is supported by our finding of an advective cosmic-ray transport such as that expected for a galactic wind. The scale-height analysis revealed an asymmetric halo above and below the disk as well as between the two sides of the major axis. A central point source as well as a bubble structure is seen in the radio data for the first time. Our X-ray data show a box-shaped hot halo around NGC 4666 and furthermore confirm the AGN nature of the central source. NGC 4666 has a large-scale X-shaped magnetic field in the halo, as has been observed in other edge-on galaxies. The analysis furthermore revealed that the disk of NGC 4666 shows hints of field reversals along its radius, which is the first detection of this phenomenon in an external galaxy.

The Implications of Simple Estimates of the 2D Outerscale Based on Measurements of Magnetic Islands for the Modulation of Galactic Cosmic-Ray Electrons

Abstract The behavior of the 2D turbulence power spectrum at the lowest wavenumbers has a significant effect on the perpendicular diffusion coefficients of charged particles in the heliosphere derived from various scattering theories, and subsequently used to model the transport of cosmic rays (CRs) and solar energetic particles. In this regard, the lengthscale at which the energy-containing range begins, as opposed to that at which the inertial range commences, is of particular interest. This 2D outerscale has, however, never before been directly observed. Recently, direct measurements of magnetic islands in the solar wind have been reported by various authors. Assuming that these may provide an estimate of the 2D ultrascale, the direct calculation of the 2D outerscale becomes possible, should an observationally motivated form for the 2D turbulence power spectrum be employed. This study presents the results of such a calculation and provides comparisons of these with previous estimates of the 2D outerscale. Furthermore, the sensitivity of galactic CR electron intensities, calculated using a 3D ab initio CR modulation model, is demonstrated, and conclusions are drawn therefrom.