Efficient cold outflows driven by cosmic rays in high-redshift galaxies and their global effects on the IGM

The origin of the most energetic cosmic ray particles is one of the most important open problems in astrophysics. Despite a big experimental effort done in the past years, the sources of these very energetic particles remain unidentified. Therefore, their distribution on the Universe and even their space density are still unknown. It has been shown that different spatial configurations of the sources lead to different energy spectra and composition profiles (in the case of sources injecting heavy nuclei) at Earth. These ensemble fluctuations are more important at the highest energies, because only nearby sources, which are necessarily few, can contribute to the flux observed at Earth. This is due to the interaction of the cosmic rays with the low energy photons of the radiation field, present in the intergalactic medium, during propagation. It is believed that the intergalactic medium is permeated by a turbulent magnetic field. Although at present it is still unknown, there are several constraints for its intensity and coherence length obtained from different observational techniques. Charged cosmic rays are affected by the intergalactic magnetic field because of the bending of their trajectories during propagation through the intergalactic medium. In this work, the influence of the intergalactic magnetic field on the ensemble fluctuations is studied. Sources injecting only protons and only iron nuclei are considered. The ensemble fluctuations are studied for different values of the density of sources compatible with the constraints recently obtained from cosmic ray data. Also, the possible detection of the ensemble fluctuations in the context of the future JEM-EUSO mission is discussed.

The role of cosmic-ray-modified contact discontinuities and pressure balance structures in two-fluid cosmic-ray hydrodynamics in one Cartesian space dimension are investigated by means of analytic and numerical solution examples, as well as by weakly nonlinear asymptotics. The fundamental wave modes of the two-fluid cosmic-ray hydrodynamic equations in the long-wavelength limit consist of the backward and forward propagating cosmic-ray-modified sound waves, with sound speed dependent on both the cosmic-ray and thermal gas pressures; the contact discontinuity; and a pressure balance mode in which the sum ofthe cosmic ray and thermal gas pressure perturbations is zero. The pressure balance mode, like the contact discontinuity is advected with the background flow. The interaction of the pressure balance mode with the contact discontinuity is investigated by means of the method of multiple scales. The thermal gas and cosmic-ray pressure perturbations satisfy a linear diffusion equation, and entropy perturbations arising from nonisentropic initial conditions for the thermal gas are frozen into the fluid. The contact discontinuity and pressure balance eigenmodes both admit nonzero perturbations in the thermal gas, whereas the cosmic-ray-modified sound waves are isentropic. The total entropy perturbation is shared between the contact discontinuity and pressure balance eigenmodes, and examples are given in which there is a transfer of entropy between the two modes. In particular, N-wave type density disturbances are obtained which arise as a result of the entropy transfer between the two modes. A weakly nonlinear geometric optics perturbation expansion is used to study the long timescale evolution of the short-wavelength entropy wave and the thermal gas sound waves in a slowly varying, large-scale background flow. The weakly nonlinear geometric optics expansion is also used to generalize previous studies of squeezing instability for short-wavelength sound waves in the two fluid model, by including a weakly nonlinear wave steepening term that leads to shock formation, as well as the effect of long time and space dependence of the background flow. Implications of cosmic-ray-modified pressure balance structures and contact discontinuities in models of the interaction of traveling interplanetary shocks and compression and rarefraction waves with the solar wind termination shock are briefly discussed.

We explore the redshift evolution of the radio luminosity function (RLF) of star-forming galaxies using galform, a semi-analytic model of galaxy formation and a dynamo model of the magnetic field evolving in a galaxy. Assuming energy equipartition between the magnetic field and cosmic rays, we derive the synchrotron luminosity of each sample galaxy. In a model where the turbulent speed is correlated with the star formation rate, the RLF is in fair agreement with observations in the redshift range 0 ≤ z ≤ 2. At larger redshifts, the structure of galaxies, their interstellar matter, and turbulence appear to be rather different from those at z ≲ 2, so that the turbulence and magnetic field models applicable at low redshifts become inadequate. The strong redshift evolution of the RLF at 0 ≤ z ≤ 2 can be attributed to an increased number, at high redshift, of galaxies with large disc volumes and strong magnetic fields. On the other hand, in models where the turbulent speed is a constant or an explicit function of z, the observed redshift evolution of the RLF is poorly captured. The evolution of the interstellar turbulence and outflow parameters appear to be major (but not the only) drivers of the RLF changes. We find that both the small- and large-scale magnetic fields contribute to the RLF but the small-scale field dominates at high redshifts. Polarization observations will therefore be important to distinguish these two components and understand better the evolution of galaxies and their non-thermal constituents.

The physics of magnetic fields (B) and cosmic rays (CRs) have recently been included in simulations of galaxy formation. However, significant uncertainties remain in how these components affect galaxy evolution. To understand their common observational tracers, we analyse the magnetic fields in a set of high-resolution, magnetohydrodynamic, cosmological simulations of Milky-Way-like galaxies from the FIRE-2 project. We compare mock observables of magnetic field tracers for simulations with and without CRs to observations of Zeeman splitting and rotation/dispersion measures. We find reasonable agreement between simulations and observations in both the neutral and the ionized interstellar medium (ISM). We find that the simulated galaxies with CRs show weaker ISM |B| fields on average compared to their magnetic-field-only counterparts. This is a manifestation of the effects of CRs in the diffuse, low density inner circumgalactic medium (CGM). We find that equipartition between magnetic and cosmic ray energy densities may be valid at large (> 1 kpc) scales for typical ISM densities of Milky-Way-like galaxies, but not in their haloes. Within the ISM, the magnetic fields in our simulated galaxies follow a power-law scaling with gas density. The scaling extends down to neutral hydrogen number densities < 300 cm−3, in contrast to observationally derived models, but consistent with the observational measurements. Finally, we generate synthetic rotation measure (RM) profiles for projections of the simulated galaxies and compare to observational constraints in the CGM. While consistent with upper limits, improved data are needed to detect the predicted CGM RMs at 10–200 kpc and better constrain theoretical predictions.

This paper investigates the distribution and implications of cosmic-ray electrons within the intergalactic medium (IGM). Utilizing a synthesis model of the extragalactic background, we evolve the spectrum of Compton-included cosmic rays. The energy density distribution of cosmic-ray electrons peaks at redshift z ≈ 2, and peaks in the ∼MeV range. The fractional contribution of cosmic-ray pressure to the general IGM pressure progressively increases toward lower redshift. At mean density, the ratio of cosmic-ray electron to thermal pressure in the IGM PCRe/Pth is 0.3% at z = 2, rising to 1.0% at z = 1, and 1.8% at z = 0.1 (considering only the cosmic rays produced locally by Compton scattering). We compute the linear Landau damping rate of plasma oscillations in the IGM caused by the ∼MeV cosmic-ray electrons, and find it to be of order ∼10−6 s−1 for wavenumbers 1.2 ≲ ck/ω p ≲ 5 at z = 2 and mean density (where ω p is the plasma frequency). This strongly affects the fate of TeV e + e − pair beams produced by blazars, which are potentially unstable to oblique instabilities involving plasma oscillations with wavenumber ck/ωp≈secθ (θ being the angle between the beam and wavevector). Linear Landau damping is at least thousands of times faster than either pair beam instability growth or collisional effects; it thus turns off the pair beam instability except for modes with very small θ (ck/ω p → 1, where linear Landau damping is kinematically suppressed). This leaves open the question of whether the pair beam instability is turned off entirely, or can still proceed via the small-θ modes.

Because cosmic rays are charged particles scrambled by magnetic fields, combining direct measurements with other observations is crucial to understanding their origin and propagation. As energetic particles traverse matter and electromagnetic fields, they leave marks in the form of neutral interaction products. Among those, γ rays trace interactions of nuclei that inelastically collide with interstellar gas, as well as of leptons that undergo Bremsstrahlung and inverse-Compton scattering. Data collected by the Fermi large area telescope (LAT) are therefore telling us the story of cosmic rays along their journey from sources through their home galaxies. Supernova remnants emerge as a notable γ-ray source population, and older remnants interacting with interstellar matter finally show strong evidence of the presence of accelerated nuclei. Yet the maximum energy attained by shock accelerators is poorly constrained by observations. Cygnus X, a massive star-forming region established by the LAT as housing cosmic-ray sources, provides a test case to study the impact of wind-driven turbulence on the early propagation. Interstellar emission resulting from the large-scale propagation of cosmic rays in the Milky Way is revealed in unprecedented detail that challenges some of the simple assumptions used for the modeling. Moreover, the cosmic-ray induced γ-ray luminosities of galaxies-scale quasi-linearly with their massive-star formation rates: the overall normalization of that relation below the calorimetric limit suggests that for most systems, a substantial fraction of energy in cosmic rays escapes into the intergalactic medium. The nuclear production models and the distribution of target gas and radiation fields, not determined precisely enough yet, are key to exploiting the full potential of γ-ray data. Nevertheless, data being collected by Fermi and complementary multiwavelength/multimessenger observations are bringing us ever closer to solving the cosmic-ray mystery.

Galactic magnetic fields and cosmic rays are tightly bound up together. The fact that the energy density of cosmic rays in the local region of the Galaxy is approximately the same order of magnitude as the energy density of the magnetic field has led many authors to postulate that this is a universal phenomenon, i.e., that there is equipartition. Gamma ray observations of the Magellanic Clouds lead us to doubt the presence of such equipartition there, however. Doubt is raised, therefore, about the application of equipartition, and a fixed electron to proton ratio, by radio astronomers to observations of synchrotron radiation from galaxies in general.A consequence of the lack of equipartition may be that in some galaxies the cosmic ray intensity is increasing (or decreasing). The positive link between the strength of galactic magnetic fields and the star formation rate implies that dynamo action may not be the only source of galactic magnetic fields.

Using data from the Voyager and Pioneer spacecraft at distances up to 60 AU from the Sun to help define the effects of solar modulation, we have made a new estimate of the local interstellar cosmic-ray spectra and the associated energy density and ionization rate. The energy density, including the low-energy electron contribution, is ~1.80 eV cm-3, equivalent energetically to that in a 7 μG magnetic field. The ionization rate from this same cosmic-ray population, including a significant contribution from low-energy electrons, is (3-4) × 10-17 s-1 (H atom). This is most likely a minimum value and is not too different from estimates of the typical ionization rate required to heat the interstellar medium and maintain the diffuse ionized gas range from (3-30) × 10-17 s-1 (H atom) as a function of locale in the Galaxy. This suggests that cosmic rays could be a significant contributor to interstellar heating. The cosmic-ray ionization rate will be more or less uniform throughout the Galactic disk and also perpendicular to the disk according to the scale height of cosmic rays which may be ~1 kpc, thus providing a source of ionization at high Z values. The recent observation of low-energy particle acceleration at the heliospheric termination shock, which produces local ionization rates comparable to those due to Galactic cosmic rays, raises the interesting possibility that "hot spots," where the ionization rate due to locally accelerated particles may exceed 10-16 s-1 (H atom), may occur near massive stars with more powerful stellar winds.

Motivated by recent experimental proposals to search for extragalactic cosmic rays (including anti-matter from distant galaxies), we study particle propagation through the intergalactic medium (IGM). We first use estimates of the magnetic field strength between galaxies to constrain the mean free path for diffusion of particles through the IGM. We then develop a simple analytic model to describe the diffusion of cosmic rays. Given the current age of galaxies, our results indicate that, in reasonable models, a completely negligible number of particles can enter our Galaxy from distances greater than $\sim 100$ Mpc for relatively low energies ($E$ $< 10^6$ GeV/n). We also find that particle destruction in galaxies along the diffusion path produces an exponential suppression of the possible flux of extragalactic cosmic rays. Finally, we use gamma ray constraints to argue that the distance to any hypothetical domains of anti-matter must be roughly comparable to the horizon scale.

The presence of heavy particles in cosmic rays makes it plausible that cosmic radiation is generated by the acceleration of charged particles in extended electromagnetic fields. The energy density of cosmic rays ($3\ifmmode\times\else\texttimes\fi{}{10}^{\ensuremath{-}10}\frac{M{c}^{2}}{{\mathrm{cm}}^{3}}$) would lead to very great total amounts of energy if one assumes that cosmic radiation extends throughout intergalactic or throughout interstellar space. This would require very efficient methods of cosmic-ray production. It is therefore of interest to investigate the possibility that cosmic rays are generated in the vicinity of the sun and are confined to the solar system by extended magnetic fields. The strength of these fields is estimated to be ${10}^{\ensuremath{-}5}$ gauss. According to this picture, cosmic rays circulate in the neighborhood of the planetary system for thousands of years, during which time the radiation becomes isotropic. The electron component of the radiation is eliminated by collisions with solar light quanta.

Understanding the physics of galaxy formation is an outstanding problem in modern astrophysics. Recent cosmological simulations have demonstrated that feedback by star formation, supernovae and active galactic nuclei appears to be critical in obtaining realistic disk galaxies and to slow down star formation to the small observed rates. However the particular physical processes underlying these feedback processes still remain elusive. In particular, many of these simulations neglected magnetic fields and relativistic particle populations (so-called cosmic rays). Those are known to provide a pressure support comparable to the thermal gas in our Galaxy and couple dynamically and thermally to the gas, which seriously questions their neglect. Early attempts to include magnetic fields and cosmic rays put the plasma kinetics of wave-particle interactions into the limelight for explaining the power of galactic winds and for regulating star formation in and cosmic accretion onto galaxies. After introducing the underlying physical concepts, I will present our recent efforts to model kinetic plasma physics in macroscopic hydrodynamical models of galaxy formation. In particular, I will explain how cosmic rays interact with and propagate through the magnetized plasma in the interstellar and circumgalactic media. I will elucidate the role of plasma instabilities in regulating the momentum transfer from cosmic rays to the background plasma and how we can observationally test these theoretical considerations using new high-sensitivity MeerKAT observations. I will then demonstrate that cosmic rays play a decisive role in the formation and evolution of spiral galaxies by providing feedback that regulates star formation and drives gas out in galactic winds. Comparing cosmic ray spectra of electrons and protons to observational data and studying the correlation of the far-infrared emission with the gamma-ray and radio emission from galaxies enables us to test the cosmic ray feedback and dynamo models for the growth of galactic magnetic fields. This argues that a complete understanding of galaxy formation necessarily includes plasma astrophysics.

The Epoch of Reionization (EoR) depends on the complex astrophysics governing the birth and evolution of the first galaxies and structures in the intergalactic medium. EoR models rely on cosmic microwave background (CMB) observations, and in particular the large-scale E-mode polarization power spectra (EE PS), to help constrain their highly uncertain parameters. However, rather than directly forward-modelling the EE PS, most EoR models are constrained using a summary statistic – the Thompson scattering optical depth, τe. Compressing CMB observations to τe requires adopting a basis set for the EoR history. The common choice is the unphysical, redshift-symmetric hyperbolic tangent (tanh) function, which differs in shape from physical EoR models based on hierarchical structure formation. Combining public EoR and CMB codes, 21cmfast and class, here we quantify how inference using the τe summary statistic impacts the resulting constraints on galaxy properties and EoR histories. Using the last Planck 2018 data release, we show that the marginalized constraints on the EoR history are more sensitive to the choice of the basis set (tanh versus physical model) than to the CMB likelihood statistic (τe versus PS). For example, EoR histories implied by the growth of structure show a small tail of partial reionization extending to higher redshifts. However, biases in inference using τe are negligible for the Planck 2018 data. Using EoR constraints from high-redshift observations including the quasar dark fraction, galaxy UV luminosity functions, and CMB EE PS, our physical model recovers $\tau _\mathrm{ e} = 0.0569_{-0.0066}^{+0.0081}$.

The diffusion of cosmic rays (CRs) in turbulent magnetic fields is fundamental to understanding various astrophysical processes. We explore the CR diffusion in the magnetic fluctuations amplified by the nonlinear turbulent dynamo in the absence of a strong mean magnetic field. Using test particle simulations, we identify three distinct CR diffusion regimes: mirroring, wandering, and magnetic moment scattering (MMS). With highly inhomogeneous distribution of the dynamo-amplified magnetic fields, we find that the diffusion of CRs is also spatially inhomogeneous. Our results reveal that lower-energy CRs preferentially undergo the mirror and wandering diffusion in the strong-field regions, and the MMS diffusion in the weak-field regions. The former two diffusion mechanisms play a more important role toward lower CR energies, resulting in a relatively weak energy dependence of the overall CR mean free path (MFP). In contrast, higher-energy CRs predominantly undergo the MMS diffusion, for which the incomplete particle gyration, i.e., the limit case of mirroring, in strong fields has a more significant effect than the scattering by small-scale field tangling/reversal. Compared with lower-energy CRs, they are more poorly confined in space and their MFPs have a stronger energy dependence. We stress the fundamental role of magnetic field inhomogeneity of nonlinear turbulent dynamo in causing the different diffusion behavior of CRs compared to that in sub-Alfvénic magnetohydrodynamic turbulence.

The circumgalactic medium (CGM) plays a critical role in galaxy evolution, influencing gas flows, feedback processes, and galactic dynamics. Observations have shown a substantial cold gas reservoir in the CGM, but the mechanisms driving its formation and evolution remain unclear. Cosmic rays (CRs), as a source of non-thermal pressure, are increasingly recognised as key regulators of cold gas dynamics. This study explores how CRs affect cold clouds that condense from the hot CGM through thermal instability (TI). Using three-dimensional CR magnetohydrodynamics simulations with the moving-mesh code Arepo, we assessed the impact of various CR transport models on cold gas evolution. Under purely advective CR transport, CR pressure significantly suppressed the collapse of thermally unstable regions, altering the CGM's structure. In contrast, our realistic CR transport models revealed that CRs escape collapsing regions via anisotropic streaming and diffusion along magnetic fields, reducing their ability to prevent collapse and diminishing their impact on the thermal structure of the cold CGM. The ratio of the CR escape timescale to the cloud collapse timescale emerged as a critical factor in determining the influence of CRs on TI. The CRs remained confined within cold clouds when effective CR diffusion was slow, thereby maximising their pressure support and inhibiting collapse. The fast and effective CR diffusion realised in our two-moment CR-magnetohydrodynamics model facilitated rapid CR escape, diminishing their stabilising effect. This realistic CR transport model shows a wide dynamic range of the effective CR diffusion coefficient; its CR-energy-weighted median ranges from 1029 to 1030 cm2 s−1 for thermally to CR-dominated atmospheres, respectively. In addition to these CR transport-related effects, we demonstrated that a high numerical resolution is crucial, as it is necessary to avoid spuriously large clouds formed in low-resolution simulations, which would result in overly long CR escape times and artificially amplified CR pressure support.

We study the change in cosmic-ray pressure, the change in cosmic-ray density, and the level of cosmic-ray induced heating via Alfven-wave damping when cosmic rays move from a hot ionized plasma to a cool cloud embedded in that plasma. The general analysis method outlined here can apply to diffuse clouds in either the ionized interstellar medium or in galactic winds. We introduce a general-purpose model of cosmic-ray diffusion building upon the hydrodynamic approximation for cosmic rays (from McKenzie & Voelk and Breitschwerdt and collaborators). Our improved method self-consistently derives the cosmic-ray flux and diffusivity under the assumption that the streaming instability is the dominant mechanism for setting the cosmic-ray flux and diffusion. We find that, as expected, cosmic rays do not couple to gas within cool clouds (cosmic rays exert no forces inside of cool clouds), that the cosmic-ray density does not increase within clouds (it may slightly decrease in general, and decrease by an order of magnitude in some cases), and that cosmic-ray heating (via Alfven-wave damping and not collisional effects as for ~10 MeV cosmic rays) is only important under the conditions of relatively strong (10 micro-Gauss) magnetic fields or high cosmic-ray pressure (~10^{-11} ergs cm^{-3}).