Cosmic-ray Energy Density Research Articles

Suppressed Cosmic-Ray Energy Densities in Molecular Clouds from Streaming Instability-regulated Transport

Cosmic rays (CRs) are the primary driver of ionization in star-forming molecular clouds (MCs). Despite their potential impacts on gas dynamics and chemistry, no simulations of star cluster formation following the creation of individual stars have included explicit cosmic-ray transport (CRT) to date. We conduct the first numerical simulations following the collapse of a 2000M ⊙ MC and the subsequent star formation including CRT using the STAR FORmation in Gaseous Environments framework implemented in the GIZMO code. We show that when CRT is streaming-dominated, the CR energy in the cloud is strongly attenuated due to energy losses from the streaming instability. Consequently, in a Milky Way–like environment the median CR ionization rate in the cloud is low (ζ ≲ 2 × 10−19 s−1) during the main star-forming epoch of the calculation and the impact of CRs on the star formation in the cloud is limited. However, in high-CR environments, the CR distribution in the cloud is elevated (ζ ≲ 6 × 10−18), and the relatively higher CR pressure outside the cloud causes slightly earlier cloud collapse and increases the star formation efficiency by 50% to ∼13%. The initial mass function is similar in all cases except with possible variations in a high-CR environment. Further studies are needed to explain the range of ionization rates observed in MCs and explore star formation in extreme CR environments.

Examining a hadronic γ-ray scenario for the radiative shell and molecular clouds of the old GeV supernova remnant G298.6−0.0

Abstract Based on the 13.7 yr Fermi-LAT data, Yeung, Bamba, and Sano (2023, PASJ, 75, 384) claimed detection of two γ-ray sources (namely Src-NE and Src-NW) associated with the supernova remnant (SNR) G298.6−0.0, and interpreted it as an old GeV SNR interacting with molecular clouds (MCs). In this follow-up study, we refine the flux measurements below 2 GeV with Fermi-LAT event types of better angular reconstruction. Then, we report our γ-ray spectral fittings and cosmic-ray phenomenology in a hadronic scenario, considering both the shell and MC regions of SNR G298.6−0.0. We confirm that the γ-ray spectra of both Src-NE and Src-NW exhibit spectral breaks, at $1.50_{-0.50}^{+0.60}$ and $0.68_{-0.11}^{+0.32}\:$GeV, respectively. Src-NW has a harder broad-band photon index than Src-NE, suggesting an appreciable difference between the physical separations of their respective emission sites from SNR G298.6−0.0. The cosmic-ray spectrum responsible for Src-NE starts with a minimum energy $E_\mathrm{CR,min}=1.38_{-0.16}^{+0.47}\:$GeV, and has a proton index $\Gamma _\mathrm{CR}=2.57_{-0.21}^{+0.18}$ below the exponential cutoff energy $E_\mathrm{CR,max}=240_{-150}^{+240}\:$GeV. Accordingly, we argue that Src-NE is dominated by the SNR shell, while only a minor portion of lower-energy emission is contributed by the MCs interacting with the SNR. The cosmic-ray population for Src-NW starts at a higher energy such that the ECR, min ratio of Src-NW to Src-NE is ≳2. The high ECR, min, as well as the high cosmic-ray energy density required (∼26 eV cm−3), supports the interpretation that Src-NW is predominantly the γ-ray emission from the farther MCs being bombarded by protons that had earlier escaped from SNR G298.6−0.0. By comparing the high-energy features of G298.6−0.0 with those of analogical SNRs, especially SNR W 28 and SNR W 44, we further constrain the age of SNR G298.6−0.0 to be 10–30 kyr, under the assumption of a purely hadronic scenario.

Cosmic-Ray Acceleration of Galactic Outflows in Multiphase Gas

We investigate the dynamical interaction between cosmic rays (CRs) and the multiphase interstellar medium (ISM) using numerical magnetohydrodynamic (MHD) simulations with a two-moment CR solver and TIGRESS simulations of star-forming galactic disks. We previously studied the transport of CRs within TIGRESS outputs using a “postprocessing” approach, and we now assess the effects of the MHD backreaction to CR pressure. We confirm our previous conclusion that there are three quite different regimes of CR transport in multiphase ISM gas, while also finding that simulations with “live MHD” predict a smoother CR pressure distribution. The CR pressure near the midplane is comparable to other pressure components in the gas, but the scale height of CRs is far larger. Next, with a goal of understanding the role of CRs in driving galactic outflows, we conduct a set of controlled simulations of the extraplanar region above z = 500 pc, with imposed boundary conditions flowing from the midplane into this region. We explore a range of thermal and kinematic properties for the injected thermal gas, encompassing both hot, fast-moving outflows, and cooler, slower-moving outflows. The boundary conditions for CR energy density and flux are scaled from the supernova rate in the underlying TIGRESS model. Our simulations reveal that CRs efficiently accelerate extraplanar material if the latter is mostly warm/warm-hot gas, in which CRs stream at the Alfvén speed, and the effective sound speed increases as density decreases. In contrast, CRs have very little effect on fast, hot outflows where the Alfvén speed is small, even when the injected CR momentum flux exceeds the injected MHD momentum flux.

Steady states of the Parker instability

ABSTRACT We study the linear properties, non-linear saturation, and a steady, strongly non-linear state of the Parker instability in galaxies. We consider magnetic buoyancy and its consequences with and without cosmic rays. Cosmic rays are described using the fluid approximation with anisotropic, non-Fickian diffusion. To avoid unphysical constraints on the instability (such as boundary conditions often used to specify an unstable background state), non-ideal magnetohydrodynamic equations are solved for deviations from a background state representing an unstable magnetohydrostatic equilibrium. We consider isothermal gas and neglect rotation. The linear evolution of the instability is in broad agreement with earlier analytical and numerical models; but we show that most of the simplifying assumptions of the earlier work do not hold, such that they provide only a qualitative rather than quantitative picture. In its non-linear stage the instability has significantly altered the background state from its initial state. Vertical distributions of both magnetic field and cosmic rays are much wider, the gas layer is thinner, and the energy densities of both magnetic field and cosmic rays are much reduced. The spatial structure of the non-linear state differs from that of any linear modes. A transient gas outflow is driven by the weakly non-linear instability as it approaches saturation.

Steady states of the Parker instability: the effects of rotation

ABSTRACT We model the Parker instability in vertically stratified isothermal gas using non-ideal MHD three-dimensional simulations. Rotation, especially differential, more strongly and diversely affects the nonlinear state than the linear stage (where we confirm the most important conclusions of analytical models), and stronger than any linear analyses predict. Steady-state magnetic fields are stronger and cosmic ray energy density is higher than in comparable non-rotating systems. Transient gas outflows induced by the nonlinear instability persist longer, of order 2 Gyr, with rotation. Stratification combined with (differential) rotation drives helical flows, leading to mean-field dynamo. Consequently, the nonlinear state becomes oscillatory (while both the linear instability and the dynamo are non-oscillatory). The horizontal magnetic field near the mid-plane reverses its direction propagating to higher altitudes as the reversed field spreads buoyantly. The spatial pattern of the large-scale magnetic field may explain the alternating magnetic field directions in the halo of the edge-on galaxy NGC 4631. Our model is unique in producing a large-scale magnetic structure similar to such observation. Furthermore, our simulations show that the mean kinetic helicity of the magnetically driven flows has the sign opposite to that in the conventional non-magnetic flows. This has profound consequences for the nature of the dynamo action and large-scale magnetic field structure in the coronae of spiral galaxies that remain to be systematically explored and understood. We show that the energy density of cosmic rays and magnetic field strength are not correlated at scales of order 1 kiloparsec.

Galaxies at a Cosmic Ray Eddington Limit

Cosmic rays have been shown to be extremely important in the dynamics of diffuse gas in galaxies, helping to maintain hydrostatic equilibrium, and serving as a regulating force in star formation. In this paper, we address the influence of cosmic rays on galaxies by re-examining the theory of a cosmic ray Eddington limit, first proposed by Socrates et al. and elaborated upon by Crocker et al. and Huang & Davis. A cosmic ray Eddington limit represents a maximum cosmic ray energy density above which the interstellar gas cannot be in hydrostatic equilibrium, resulting in a wind. In this paper, we continue to explore the idea of a cosmic ray Eddington limit by introducing a general framework that accounts for the circumgalactic environment and applying it to five galaxies that we believe to be a good representative sample of the star-forming galaxy population, using different cosmic ray transport models to determine what gives each galaxy the best chance to reach this limit. We show that, while an Eddington limit for cosmic rays does exist, for our five galaxies, the limit either falls at star formation rates that are much larger or gas densities that are much lower than each galaxy’s measured values. This suggests that cosmic ray pressure is not the main factor limiting the luminosity of starburst galaxies.

Standard self-confinement and extrinsic turbulence models for cosmic ray transport are fundamentally incompatible with observations

ABSTRACT Models for cosmic ray (CR) dynamics fundamentally depend on the rate of CR scattering from magnetic fluctuations. In the ISM, for CRs with energies ∼MeV-TeV, these fluctuations are usually attributed either to ‘extrinsic turbulence’ (ET) – a cascade from larger scales – or ‘self-confinement’ (SC) – self-generated fluctuations from CR streaming. Using simple analytic arguments and detailed ‘live’ numerical CR transport calculations in galaxy simulations, we show that both of these, in standard form, cannot explain even basic qualitative features of observed CR spectra. For ET, any spectrum that obeys critical balance or features realistic anisotropy, or any spectrum that accounts for finite damping below the dissipation scale, predicts qualitatively incorrect spectral shapes and scalings of B/C and other species. Even if somehow one ignored both anisotropy and damping, observationally required scattering rates disagree with ET predictions by orders of magnitude. For SC, the dependence of driving on CR energy density means that it is nearly impossible to recover observed CR spectral shapes and scalings, and again there is an orders-of-magnitude normalization problem. But more severely, SC solutions with super-Alfvénic streaming are unstable. In live simulations, they revert to either arbitrarily rapid CR escape with zero secondary production, or to bottleneck solutions with far-too-strong CR confinement and secondary production. Resolving these fundamental issues without discarding basic plasma processes requires invoking different drivers for scattering fluctuations. These must act on a broad range of scales with a power spectrum obeying several specific (but plausible) constraints.

Magnetic fields on FIRE: Comparing B-fields in the multiphase ISM and CGM of simulated L* galaxies to observations

ABSTRACT 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.

The Coma Cluster at LOFAR Frequencies. II. The Halo, Relic, and a New Accretion Relic

We present LOw Frequency ARray observations of the Coma Cluster field at 144 MHz. The cluster hosts one of the most famous radio halos, a relic, and a low surface brightness bridge. We detect new features that allow us to make a step forward in the understanding of particle acceleration in clusters. The radio halo extends for more than 2 Mpc, which is the largest extent ever reported. To the northeast of the cluster, beyond the Coma virial radius, we discover an arc-like radio source that could trace particles accelerated by an accretion shock. To the west of the halo, coincident with a shock detected in the X-rays, we confirm the presence of a radio front, with different spectral properties with respect to the rest of the halo. We detect a radial steepening of the radio halo spectral index between 144 and 342 MHz, at ∼30′ from the cluster center, that may indicate a non-constant re-acceleration time throughout the volume. We also detect a mild steepening of the spectral index toward the cluster center. For the first time, a radial change in the slope of the radio–X-ray correlation is found, and we show that such a change could indicate an increasing fraction of cosmic-ray versus thermal energy density in the cluster outskirts. Finally, we investigate the origin of the emission between the relic and the source NGC 4789, and we argue that NGC 4789 could have crossed the shock originating the radio emission visible between its tail and the relic.

LHAASO and Galactic cosmic rays

Cosmic rays (CRs) are relativistic charged particles in the interstellar medium (ISM). They are mainly protons (hydrogen nuclei) with about 10% fraction of helium nuclei and smaller abundances of heavier elements. CRs are an important component in ISM. In our galaxy, the energy density of CRs is similar to that of the magnetic fields and interstellar radiation fields (ISRFs). CRs also determine the ionization rate and heating of gas in the dense core of molecular clouds, thus playing a leading role in the astro-chemistry processes therein and controling the star-forming processes.

Galactic center gamma-ray production by cosmic rays from stellar winds and Sgr A East

Context. The High Energy Stereoscopic System, the Major Atmospheric Gamma-ray Imaging Cherenkov Telescope, and the Very Energetic Radiation Imaging Telescope Array System have observed diffuse gamma-ray emission strongly correlated with the central molecular zone in the Galactic center. The most accepted scenario to generate this emission is via a hadronic interaction between cosmic rays (CRs) and ambient gas, where CRs are accelerated from a central and continuous source of 1 PeV protons. Aims. We explore the influence of the three-dimensional (3D) shape of the central molecular zone on the indirect observation of the CR energy density via gamma-ray detection. Methods. We simulated synthetic gamma-ray maps using a CR diffusion model with spherical injection, one isotropic diffusion coefficient, no advection, and mono-energetic particles of 1 PeV. Also, we used two different 3D gas distributions considering the observed gas column density, both with and without an inner cavity. Results. We find that when using a persistent CR source, a disk-like gas distribution is needed to reproduce the existing CR indirect observations. This is in agreement with the continuous gas distribution implied by some dynamical models and studies based on the comparison of emission and absorption molecular lines. However, it contradicts several models of the central molecular zone, which imply that this structure has a significant inner cavity. This tension can be reconciled by an additional, impulsive CR injection. Conclusions. If the central molecular zone has a cavity, a composite CR population, coming from the stellar winds of the Wolf-Rayet stars in the central 0.5 pc and the supernova Sgr A East, produces a good match to the observed gamma-ray morphology in the Galactic center.

The 1.28 GHz MeerKAT Galactic Center Mosaic

Abstract The inner ∼200 pc region of the Galaxy contains a 4 million M⊙ supermassive black hole (SMBH), significant quantities of molecular gas, and star formation and cosmic-ray energy densities that are roughly two orders of magnitude higher than the corresponding levels in the Galactic disk. At a distance of only 8.2 kpc, the region presents astronomers with a unique opportunity to study a diverse range of energetic astrophysical phenomena, from stellar objects in extreme environments, to the SMBH and star-formation-driven feedback processes that are known to influence the evolution of galaxies as a whole. We present a new survey of the Galactic center conducted with the South African MeerKAT radio telescope. Radio imaging offers a view that is unaffected by the large quantities of dust that obscure the region at other wavelengths, and a scene of striking complexity is revealed. We produce total-intensity and spectral-index mosaics of the region from 20 pointings (144 hr on-target in total), covering 6.5 square degrees with an angular resolution of 4″ at a central frequency of 1.28 GHz. Many new features are revealed for the first time due to a combination of MeerKAT’s high sensitivity, exceptional u, v-plane coverage, and geographical vantage point. We highlight some initial survey results, including new supernova remnant candidates, many new nonthermal filament complexes, and enhanced views of the Radio Arc bubble, Sagittarius A, and Sagittarius B regions. This project is a South African Radio Astronomy Observatory public legacy survey, and the image products are made available with this article.

Cosmic rays from massive star clusters: a close look at Westerlund 1

ABSTRACT We study the effect of cosmic ray (CR) acceleration in the massive compact star cluster Westerlund 1 in light of its recent detection in γ-rays. Recent observations reveal a 1/r radial distribution of the CR energy density. Here, we theoretically investigate whether or not this profile can help to distinguish between (1) continuous CR acceleration in the star cluster stellar wind-driven shocks and (2) discrete CR acceleration in multiple supernovae (SNe) shocks – which are often debated in the literature. Using idealized two-fluid simulations and exploring different acceleration sites and diffusion coefficients, we obtain the CR energy density profile and luminosity to find the best match for the γ-ray observations. We find that the inferred CR energy density profiles from observations of γ-ray luminosity and mass can be much different from the true radial profile. CR acceleration at either the cluster core region or the wind termination shock can explain the observations, if the diffusion coefficient is κcr ∼ 1027 cm2 s−1 and a fraction of ${\approx}10-20{{\ \rm per\ cent}}$ of the shock power/post-shock pressure is deposited into the CR component. We also study the possibility of discrete SNe explosions being responsible for CR acceleration and find that with an injection rate of 1 SN in every ∼0.03 Myr, one can explain the observed γ-ray profile. This multiple SN scenario is consistent with X-ray observations only if the thermal conductivity is close to the Spitzer value.

A GeV-TeV particle component and the barrier of cosmic-ray sea in the Central Molecular Zone

Cosmic rays are important probe of a number of fundamental physical problems such as the acceleration of high and very high energy particles in extreme astrophysical environments. The Galactic center is widely anticipated to be an important cosmic-ray source and the observations of some Imaging Atmospheric Cherenkov Telescopes did successfully reveal a component of TeV-PeV cosmic rays in the vicinity of the Galactic center. Here we report the identification of GeV-TeV cosmic rays in the central molecular zone with the γ-ray observations of the Fermi Large Area Telescope, whose spectrum and spatial gradient are consistent with that measured by the Imaging Atmospheric Cherenkov Telescopes but the corresponding cosmic-ray energy density is substantially lower than the so-called cosmic-ray sea component, suggesting the presence of a high energy particle accelerator at the Galactic center and the existence of a barrier that can effectively suppress the penetration of the particles from the cosmic-ray sea to the central molecular zone.

Cosmic rays and non-thermal emission in simulated galaxies – III. Probing cosmic-ray calorimetry with radio spectra and the FIR–radio correlation

ABSTRACT An extinction-free estimator of the star formation rate (SFR) of galaxies is critical for understanding the high-redshift universe. To this end, the nearly linear, tight correlation of far-infrared (FIR), and radio luminosity of star-forming galaxies is widely used. While the FIR is linked to massive star formation, which also generates shock-accelerated cosmic-ray (CR) electrons and radio synchrotron emission, a detailed understanding of the underlying physics is still lacking. Hence, we perform three-dimensional magnetohydrodynamical (MHD) simulations of isolated galaxies over a broad range of halo masses and SFRs using the moving-mesh code arepo, and evolve the CR proton energy density self-consistently. In post-processing, we calculate the steady-state spectra of primary, shock-accelerated and secondary CR electrons, which result from hadronic CR proton interactions with the interstellar medium. The resulting total radio luminosities correlate with the FIR luminosities as observed and are dominated by primary CR electrons if we account for anisotropic CR diffusion. The increasing contribution of secondary emission up to 30 per cent in starbursts is compensated by the larger bremsstrahlung and Coulomb losses. CR electrons are in the calorimetric limit and lose most of their energy through inverse Compton interactions with star light and cosmic microwave background (CMB) photons while less energy is converted into synchrotron emission. This implies steep steady-state synchrotron spectra in starbursts. Interestingly, we find that thermal free–free emission flattens the total radio spectra at high radio frequencies and reconciles calorimetric theory with observations while free–free absorption explains the observed low-frequency flattening towards the central regions of starbursts.

Probing the Sea of Cosmic Rays by Measuring Gamma-Ray Emission from Passive Giant Molecular Clouds with HAWC

Abstract The study of high-energy gamma rays from passive giant molecular clouds (GMCs) in our Galaxy is an indirect way to characterize and probe the paradigm of the “sea” of cosmic rays in distant parts of the Galaxy. By using data from the High Altitude Water Cerenkov (HAWC) Observatory, we measure the gamma-ray flux above 1 TeV of a set of these clouds to test the paradigm. We selected high galactic latitude clouds that are in HAWC’s field of view and that are within 1 kpc distance from the Sun. We find no significant excess emission in the cloud regions, nor when we perform a stacked log-likelihood analysis of GMCs. Using a Bayesian approach, we calculate 95% credible interval upper limits of the gamma-ray flux and estimate limits on the cosmic-ray energy density of these regions. These are the first limits to constrain gamma-ray emission in the multi-TeV energy range (>1 TeV) using passive high galactic latitude GMCs. Assuming that the main gamma-ray production mechanism is due to proton–proton interaction, the upper limits are consistent with a cosmic-ray flux and energy density similar to that measured at Earth.

The GADOT Galaxy Survey: Dense Gas and Feedback in Herschel-selected Starburst Galaxies at Redshifts 2 to 6

Abstract We report the detection of 23 OH+ 1 → 0 absorption, emission, or P-Cygni-shaped lines and CO(J = 9→8) emission lines in 18 Herschel-selected z = 2–6 starburst galaxies with the Atacama Large Millimeter/submillimeter Array and the NOrthern Extended Millimeter Array, taken as part of the Gas And Dust Over cosmic Time Galaxy Survey. We find that the CO(J = 9→8) luminosity is higher than expected based on the far-infrared luminosity when compared to nearby star-forming galaxies. Together with the strength of the OH+ emission components, this may suggest that shock excitation of warm, dense molecular gas is more prevalent in distant massive dusty starbursts than in nearby star-forming galaxies on average, perhaps due to an impact of galactic winds on the gas. OH+ absorption is found to be ubiquitous in massive high-redshift starbursts, and is detected toward 89% of the sample. The majority of the sample shows evidence for outflows or inflows based on the velocity shifts of the OH+ absorption/emission, with a comparable occurrence rate of both at the resolution of our observations. A small subsample appears to show outflow velocities in excess of their escape velocities. Thus, starburst-driven feedback appears to be important in the evolution of massive galaxies in their most active phases. We find a correlation between the OH+ absorption optical depth and the dust temperature, which may suggest that warmer starbursts are more compact and have higher cosmic-ray energy densities, leading to more efficient OH+ ion production. This is in agreement with a picture in which these high-redshift galaxies are “scaled-up” versions of the most intense nearby starbursts.

SILCC VI – Multiphase ISM structure, stellar clustering, and outflows with supernovae, stellar winds, ionizing radiation, and cosmic rays

ABSTRACT We present simulations of the multiphase interstellar medium (ISM) at solar neighbourhood conditions including thermal and non-thermal ISM processes, star cluster formation, and feedback from massive stars: stellar winds, hydrogen ionizing radiation computed with the novel treeray radiative transfer method, supernovae (SN), and the injection of cosmic rays (CR). N-body dynamics is computed with a 4th-order Hermite integrator. We systematically investigate the impact of stellar feedback on the self-gravitating ISM with magnetic fields, CR advection and diffusion, and non-equilibrium chemical evolution. SN-only feedback results in strongly clustered star formation with very high star cluster masses, a bi-modal distribution of the ambient SN densities, and low volume-filling factors (VFF) of warm gas, typically inconsistent with local conditions. Early radiative feedback prevents an initial starburst, reduces star cluster masses and outflow rates. Furthermore, star formation rate surface densities of $\Sigma _{\dot{M}_\star } = 1.4-5.9 \times 10^{-3}$$\mathrm{M}_\odot \, \mathrm{yr}^{-1}\, \mathrm{kpc}^{-2}$, VFFwarm = 60–80 per cent as well as thermal, kinetic, magnetic, and cosmic ray energy densities of the model including all feedback mechanisms agree well with observational constraints. On the short, 100 Myr, time-scales investigated here, CRs only have a moderate impact on star formation and the multiphase gas structure and result in cooler outflows, if present. Our models indicate that at low gas surface densities SN-only feedback only captures some characteristics of the star-forming ISM and outflows/inflows relevant for regulating star formation. Instead, star formation is regulated on star cluster scales by radiation and winds from massive stars in clusters, whose peak masses agree with solar neighbourhood estimates.

Cosmic rays across the star-forming galaxy sequence – I. Cosmic ray pressures and calorimetry

ABSTRACT In the Milky Way (MW), cosmic rays (CRs) are dynamically important in the interstellar medium (ISM), contribute to hydrostatic balance, and may help regulate star formation. However, we know far less about the importance of CRs in galaxies whose gas content or star formation rate (SFR) differ significantly from those of the MW. Here, we construct self-consistent models for hadronic CR transport, losses, and contribution to pressure balance as a function of galaxy properties, covering a broad range of parameters from dwarfs to extreme starbursts. While the CR energy density increases from ∼1 eV cm−3 to ∼1 keV cm−3 over the range from sub-MW dwarfs to bright starbursts, strong hadronic losses render CRs increasingly unimportant dynamically as the SFR surface density increases. In MW-like systems, CR pressure is typically comparable to turbulent gas and magnetic pressure at the galactic mid-plane, but the ratio of CR to gas pressure drops to ∼10−3 in dense starbursts. Galaxies also become increasingly CR calorimetric and gamma-ray bright in this limit. The degree of calorimetry at fixed galaxy properties is sensitive to the assumed model for CR transport, and in particular to the time CRs spend interacting with neutral ISM, where they undergo strong streaming losses. We also find that in some regimes of parameter space hydrostatic equilibrium discs cannot exist, and in Paper II of this series we use this result to derive a critical surface in the plane of star formation surface density and gas surface density beyond which CRs may drive large-scale galactic winds.

Photodissociation region diagnostics across galactic environments

ABSTRACT We present three-dimensional astrochemical simulations and synthetic observations of magnetized, turbulent, self-gravitating molecular clouds. We explore various galactic interstellar medium environments, including cosmic ray ionization rates in the range of ζCR = 10−17–$10^{-14}\, {\rm s}^{-1}$, far-UV intensities in the range of G0 = 1–103 and metallicities in the range of Z = 0.1–$2\, {\rm Z}_{\odot }$. The simulations also probe a range of densities and levels of turbulence, including cases where the gas has undergone recent compression due to cloud–cloud collisions. We examine: (i) the column densities of carbon species across the cycle of C ii, C i, and CO, along with O i, in relation to the H i-to-H2 transition; (ii) the velocity-integrated emission of [C ii] 158 μm, [13C ii] 158 μm, [C i] 609 μm and 370 μm, [O i] 63 μm and 146 μm, and of the first ten 12CO rotational transitions; (iii) the corresponding Spectral Line Energy Distributions; (iv) the usage of [C ii] and [O i] 63 μm to describe the dynamical state of the clouds; (v) the behaviour of the most commonly used ratios between transitions of CO and [C i]; and (vi) the conversion factors for using CO and C i as H2-gas tracers. We find that enhanced cosmic ray energy densities enhance all aforementioned line intensities. At low metallicities, the emission of [C ii] is well connected with the H2 column, making it a promising new H2 tracer in metal-poor environments. The conversion factors of XCO and XC i depend on metallicity and the cosmic ray ionization rate, but not on FUV intensity. In the era of ALMA, SOFIA, and the forthcoming CCAT-prime telescope, our results can be used to understand better the behaviour of systems in a wide range of galactic and extragalactic environments.