Cosmic Ray Heating Research Articles

Fermi and eROSITA Bubbles as Persistent Structures of the Milky Way

The Fermi and eROSITA bubbles (FBs and eRBs), large diffuse structures in our Galaxy, can be the by-products of steady star formation activity. To simultaneously explain the star formation history of the Milky Way (MW) and the metallicity of ∼Z ⊙ at the Galactic disk, a steady Galactic wind driven by cosmic rays (CRs) is required. For tenuous gases with a density of ≲10−3 cm−3, CR heating dominates over radiative cooling, and the gas can maintain the virial temperature of ∼0.3 keV, ideal for escape from the Galactic system as the wind. A part of the wind falls back onto the disk like a Galactic fountain flow. We model the wind dynamics according to the Galactic evolution scenario and find that the scale height and surface brightness of the X-ray and the hadronic gamma-ray emissions from such fountain flow region can be consistent with the observed properties of the FBs and eRBs. This implies that the bubbles are persistent structures of the MW existing over (at least) the last ∼1 Gyr rather than evanescent structures formed by nontrivial, ∼10 Myr past Galactic center transient activities.

Effect of cosmic rays and ionizing radiation on observational ultraviolet plasma diagnostics in the circumgalactic medium

ABSTRACT The relevance of some galactic feedback mechanisms, in particular cosmic-ray (CR) feedback and the hydrogen ionizing radiation field, has been challenging to definitively describe in a galactic context, especially far outside the galaxy in the circumgalactic medium (CGM). Theoretical and observational uncertainties prevent conclusive interpretations of multiphase CGM properties derived from ultraviolet (UV) diagnostics. We conduct three-dimensional magnetohydrodynamic simulations of a section of a galactic disc with star formation and feedback, including radiative heating from stars, a UV background, and CR feedback. We utilize the temperature phases present in our simulations to generate Cloudy models to derive spatially and temporally varying synthetic UV diagnostics. We find that radiative effects without additional heating mechanisms are not able to produce synthetic diagnostics in the observed ranges. For low CR diffusivity $\kappa _{\rm {cr}}=10^{28} \rm {cm}^2 \rm {s}^{-1}$, CR streaming heating in the outflow helps our synthetic line ratios roughly match observed ranges by producing transitional temperature gas (T ∼ 105–106 K). High CR diffusivity $\kappa _{\rm {cr}}=10^{29} \rm {cm}^2 \rm {s}^{-1}$, with or without CR streaming heating, produced transitional temperature gas. The key parameter controlling the production of this gas phase remains unclear, as the different star formation history and outflow evolution itself influences these diagnostics. Our work demonstrates the use of UV plasma diagnostics to differentiate between galactic/circumgalactic feedback models.

ALMA Observations of Supernova Remnant N49 in the Large Magellanic Cloud. II. Non-LTE Analysis of Shock-heated Molecular Clouds

We present the first compelling evidence of shock-heated molecular clouds associated with the supernova remnant (SNR) N49 in the Large Magellanic Cloud (LMC). Using 12CO(J = 2–1, 3–2) and 13CO(J = 2–1) line emission data taken with the Atacama Large Millimeter/Submillimeter Array, we derived the H2 number density and kinetic temperature of eight 13CO-detected clouds using the large velocity gradient approximation at a resolution of 3.″5 (∼0.8 pc at the LMC distance). The physical properties of the clouds are divided into two categories: three of them near the shock front show the highest temperatures of ∼50 K with densities of ∼500–700 cm−3, while other clouds slightly distant from the SNR have moderate temperatures of ∼20 K with densities of ∼800–1300 cm−3. The former clouds were heated by supernova shocks, but the latter were dominantly affected by the cosmic-ray heating. These findings are consistent with the efficient production of X-ray recombining plasma in N49 due to thermal conduction between the cold clouds and hot plasma. We also find that the gas pressure is roughly constant except for the three shock-engulfed clouds inside or on the SNR shell, suggesting that almost no clouds have evaporated within the short SNR age of ∼4800 yr. This result is compatible with the shock-interaction model with dense and clumpy clouds inside a low-density wind bubble.

Signatures of cosmic ray heating in 21-cm observables

ABSTRACT Cosmic rays generated by supernovae carry away a significant portion of the lifetime energy emission of their parent star, making them a plausible mechanism for heating the early universe intergalactic medium (IGM). Following a review of the existing literature on cosmic ray heating, we develop a flexible model of this heating mechanism for use in 3D seminumerical 21-cm signal simulations and conduct the first investigations of the signatures it imprints on the 21-cm power spectrum and tomographic maps. We find that cosmic ray heating of the IGM is short-ranged, leading to heating clustered around star-forming sites, and a sharp contrast between heated regions of 21-cm emission and unheated regions of absorption. This contrast results in greater small-scale power for cosmic ray heated scenarios compared to what is found for X-ray heating, thus suggesting a way to test the nature of IGM heating with future 21-cm observations. Finally, we find an unexpectedly rich thermal history in models where cosmic rays can only escape efficiently from low-mass haloes, such as in scenarios where these energetic particles originate from population III star supernovae remnants. The interplay of heating and the Lyman–Werner feedback in these models can produce a local peak in the IGM kinetic temperature and, for a limited parameter range, a flattened absorption trough in the global 21-cm signal.

The impact of cosmic rays on thermal and hydrostatic stability in galactic haloes

ABSTRACT We investigate how cosmic rays (CRs) affect thermal and hydrostatic stability of circumgalactic (CGM) gas, in simulations with both CR streaming and diffusion. Local thermal instability can be suppressed by CR-driven entropy mode propagation, in accordance with previous analytic work. However, there is only a narrow parameter regime where this operates, before CRs overheat the background gas. As mass dropout from thermal instability causes the background density and hence plasma β ≡ Pg/PB to fall, the CGM becomes globally unstable. At the cool disc-to-hot−halo interface, a sharp drop in density boosts Alfven speeds and CR gradients, driving a transition from diffusive to streaming transport. CR forces and heating strengthen, while countervailing gravitational forces and radiative cooling weaken, resulting in a loss of both hydrostatic and thermal equilibrium. In lower β haloes, CR heating drives a hot, single-phase diffuse wind with velocities v ∝ (theat/tff)−1, which exceeds the escape velocity when theat/tff ≲ 0.4. In higher β haloes, where the Alfven Mach number is higher, CR forces drive multi-phase winds with cool, dense fountain flows and significant turbulence. These flows are CR dominated due to ‘trapping’ of CRs by weak transverse B-fields, and have the highest mass loading factors. Thus, local thermal instability can result in winds or fountain flows where either the heat or momentum input of CRs dominates.

Cosmic-ray driven galactic winds from the warm interstellar medium

ABSTRACT We study the properties of cosmic-ray (CR) driven galactic winds from the warm interstellar medium using idealized spherically symmetric time-dependent simulations. The key ingredients in the model are radiative cooling and CR-streaming-mediated heating of the gas. Cooling and CR heating balance near the base of the wind, but this equilibrium is thermally unstable, leading to a multiphase wind with large fluctuations in density and temperature. In most of our simulations, the heating eventually overwhelms cooling, leading to a rapid increase in temperature and a thermally driven wind; the exception to this is in galaxies with the shallowest potentials, which produce nearly isothermal $T \approx 10^4\,$ K winds driven by CR pressure. Many of the time-averaged wind solutions found here have a remarkable critical point structure, with two critical points. Scaled to real galaxies, we find mass outflow rates $\dot{M}$ somewhat larger than the observed star-formation rate in low-mass galaxies, and an approximately ‘energy-like’ scaling $\dot{M} \propto v_{\rm esc}^{-2}$. The winds accelerate slowly and reach asymptotic wind speeds of only ∼0.4vesc. The total wind power is $\sim 1~{{\ \rm per\ cent}}$ of the power from supernovae, suggesting inefficient preventive CR feedback for the physical conditions modelled here. We predict significant spatially extended emission and absorption lines from 104–105.5 K gas; this may correspond to extraplanar diffuse ionized gas seen in star-forming galaxies.

Resistive heating induced by streaming cosmic rays around a galaxy in the early Universe

ABSTRACTIt is expected that cosmic rays (CRs) escape from high-redshift galaxies at redshift z ∼ 10–20 because CRs are accelerated by supernova remnants of the first stars. Although ultraviolet and X-ray photons are widely considered the main source of heating of the intergalactic medium, CRs can also contribute to it. When the CRs propagate in the intergalactic medium, in addition to the heating process due to CR ionization, resistive heating occurs due to the electron return current induced by the streaming CRs. We evaluate the heating rate around a galaxy as a function of the distance from the galaxy. We find that the resistive heating induced by CRs dominates over the other heating processes in the vicinity of the galaxy ($r \lesssim 10^2 \, \mathrm{kpc}$) until the temperature reaches $T\sim 10^4 \, \mathrm{K}$. We also recalculate the strength of the magnetic field generated by streaming CRs under the presence of X-ray heating and show that achieved strength can be about one order of magnitude smaller when the X-ray heating is included. The presence of the ‘first’ CRs could be confirmed from the characteristic signature of CR heating imprinted on the 21-cm line map in future radio observations.

Impact of cosmic rays on the global 21-cm signal during cosmic dawn

ABSTRACT It is extremely important to understand the processes through which the thermal state of the inter-galactic medium (IGM) evolved in the early universe in order to study the evolution of HI 21-cm signal during cosmic dawn. Here, we consider the heating of the IGM due to cosmic ray (CR) protons generated by the supernovae from both early Pop III and Pop II stars. The low energy CR (cr protons from Pop III supernovae can escape from minihalos and heat the IGM via collision and ionization of hydrogen. Furthermore, high-energy protons generated in Pop II supernovae can escape the hosting halos and heat the IGM via magnetosonic Alfvén waves. We show that the heating due to these CR particles can significantly impact the IGM temperature and hence the global 21-cm signal at z ∼ 14−18. The depth, location, and duration of the 21-cm absorption profile are highly dependent on the efficiencies of CR heating. In particular, the Experiment to Detect the Epoch of Reionization Signature signal can be well fitted by the CR heating along with the Lyman-α coupling and the dark matter--baryon interaction that we consider to achieve a ‘colder IGM background’. Further, we argue that the properties of CRs and the nature of first generation of stars could be constrained by accurately measuring the global 21-cm absorption signal during the cosmic dawn.

Photochemistry and Heating/Cooling of the Multiphase Interstellar Medium with UV Radiative Transfer for Magnetohydrodynamic Simulations

We present an efficient heating/cooling method coupled with chemistry and UV radiative transfer that can be applied to numerical simulations of the interstellar medium (ISM). We follow the time-dependent evolution of hydrogen species (H2, H, H+), assume carbon/oxygen species (C, C+, CO, O, and O+) are in formation–destruction balance given the nonsteady hydrogen abundances, and include essential heating/cooling processes needed to capture the thermodynamics of all ISM phases. UV radiation from discrete point sources and the diffuse background is followed through adaptive ray tracing and a six-ray approximation, respectively, allowing for H2 self-shielding; cosmic-ray heating and ionization are also included. To validate our methods and demonstrate their application for a range of density, metallicity, and radiation fields, we conduct a series of tests, including the equilibrium curves of thermal pressure versus density, the chemical and thermal structure in photodissociation regions, H i-to-H2 transitions, and the expansion of H ii regions and radiative supernova remnants. Careful treatment of photochemistry and cosmic-ray ionization is essential for many aspects of ISM physics, including identifying the thermal pressure at which cold and warm neutral phases coexist. We caution that many current heating and cooling treatments used in galaxy formation simulations do not reproduce the correct thermal pressure and ionization fraction in the neutral ISM. Our new model is implemented in the MHD code Athena and incorporated in the TIGRESS simulation framework, for use in studying the star-forming ISM in a wide range of environments.

Tracing Interstellar Heating: An ALCHEMI Measurement of the HCN Isomers in NGC 253

We analyze HCN and HNC emission in the nearby starburst galaxy NGC 253 to investigate its effectiveness in tracing heating processes associated with star formation. This study uses multiple HCN and HNC rotational transitions observed using the Atacama Large Millimeter/submillimeter Array via the ALCHEMI Large Program. To understand the conditions and associated heating mechanisms within NGC 253's dense gas, we employ Bayesian nested sampling techniques applied to chemical and radiative transfer models, which are constrained using our HCN and HNC measurements. We find that the volume density and cosmic-ray ionization rate (CRIR) ζ are enhanced by about an order of magnitude in the galaxy’s central regions as compared to those further from the nucleus. In NGC 253's central giant molecular clouds (GMCs), where observed HCN/HNC abundance ratios are the lowest, n ∼ 105.5 cm−3 and ζ ∼ 10−12 s−1 (greater than 104 times the average Galactic rate). We find a positive correlation in the association of both density and CRIR with the number of star formation-related heating sources (supernova remnants, H ii regions, and super hot cores) located in each GMC, as well as a correlation between CRIRs and supernova rates. Additionally, we see an anticorrelation between the HCN/HNC ratio and CRIR, indicating that this ratio will be lower in regions where ζ is higher. Though previous studies suggested HCN and HNC may reveal strong mechanical heating processes in NGC 253's CMZ, we find cosmic-ray heating dominates the heating budget, and mechanical heating does not play a significant role in the HCN and HNC chemistry.

Towards a better understanding of ice mantle desorption by cosmic rays

ABSTRACT The standard model of cosmic ray heating-induced desorption of interstellar ices is based on a continuous representation of the sporadic desorption of ice mantle components from classical ($0.1\, \mu$m) dust grains. This has been re-evaluated and developed to include tracking the desorption through (extended) grain cooling profiles, consideration of grain size-dependencies and constraints to the efficiencies. A model was then constructed to study the true, sporadic, nature of the process with possible allowances from species co-desorption and whole mantle desorption from very small grains. The key results from the study are that the desorption rates are highly uncertain, but almost certainly significantly larger than have been previously determined. For typical interstellar grain size distributions it is found that the desorption is dominated by the contributions from the smallest grains. The sporadic desorption model shows that, if the interval between cosmic ray impacts is comparable to, or less than, the freeze-out time-scale, the continuous representation is inapplicable; chemical changes may occur on very long time-scales, resulting in strong gas phase chemical enrichments that have very non-linear dependences on the cosmic ray flux. The inclusion of even limited levels of species co-desorption and/or the contribution from very small grains further enhances the rates, especially for species such as H2O. In general, we find that cosmic ray heating is the dominant desorption mechanism in dark environments. These results may have important chemical implications for protostellar and protoplanetary environments.

The cosmic-ray staircase: the outcome of the cosmic-ray acoustic instability

ABSTRACT Recently, cosmic rays (CRs) have emerged as a leading candidate for driving galactic winds. Small-scale processes can dramatically affect global wind properties. We run two-moment simulations of CR streaming to study how sound waves are driven unstable by phase-shifted CR forces and CR heating. We verify linear theory growth rates. As the sound waves grow non-linear, they steepen into a quasi-periodic series of propagating shocks; the density jumps at shocks create CR bottlenecks. The depth of a propagating bottleneck depends on both the density jump and its velocity; ΔPc is smaller for rapidly moving bottlenecks. A series of bottlenecks creates a CR staircase structure, which can be understood from a convex hull construction. The system reaches a steady state between growth of new perturbations, and stair mergers. CRs are decoupled at plateaus, but exert intense forces and heating at stair jumps. The absence of CR heating at plateaus leads to cooling, strong gas pressure gradients and further shocks. If bottlenecks are stationary, they can drastically modify global flows; if their propagation times are comparable to dynamical times, their effects on global momentum and energy transfer are modest. The CR acoustic instability is likely relevant in thermal interfaces between cold and hot gas, as well as galactic winds. Similar to increased opacity in radiative flows, the build-up of CR pressure due to bottlenecks can significantly increase mass outflow rates, by up to an order of magnitude. It seeds unusual forms of thermal instability, and the shocks could have distinct observational signatures, on ∼kpc scales.

HI and H2 gas evolution over cosmic times: ColdSIM

We present first results of cold cosmic gas evolution obtained through a set of state-of-the-art numerical simulations (ColdSIM). We model time-dependent atomic and molecular non-equilibrium chemistry coupled to hydrodynamics, star formation, feedback effects, various UV backgrounds as suggested by the recent literature, HI and H2 self-shielding, H2 dust grain catalysis, photoelectric heating and cosmic-ray heating. By means of such nonequilibriumcalculations we are finally able to reproduce the latest HI and H2 observational data. Consistently with available determinations, neutral-gas mass density parameter results around Ωneutral ∼ 10−3 and increases from lower to higher redshift (z). The molecular-gas mass density parameter shows peak values of ΩH2 ∼ 10−4, while expected H2 fractions can be as high as 50% of the cold gas mass at z ∼ 4-8, in line with the latest measurements from high-z galaxies. These values agree with observations up to z ∼ 7 and both HI and H2 trends are well reproduced by our non-equilibrium H2-based star formation modelling. Corresponding H2 depletion times remain below the Hubble time and comparable to the dynamical time at all epochs. This implies that non-equilibrium molecular cooling is efficient at driving cold-gas collapse in a variety of environments and since the first half Gyr. Our findings suggest that, besides HI, non-equilibrium H2 analyses are key probes for assessing cold gas and the role of UV background radiation.

Atomic and molecular gas from the epoch of reionisation down to redshift 2

Context. Cosmic gas makes up about 90% of the baryonic matter in the Universe and H2 molecule is the most tightly linked to star formation. Aims. In this work we study cold neutral gas, its H2 component at different epochs, and corresponding depletion times. Methods. We perform state-of-the-art hydrodynamic simulations that include time-dependent atomic and molecular non-equilibrium chemistry coupled to star formation, feedback effects, different UV backgrounds presented in the recent literature and a number of additional processes occurring during structure formation (COLDSIM). We predict gas evolution and contrast the mass density parameters and gas depletion timescales. We also investigate their relation to cosmic expansion in light of the latest infrared and (sub)millimetre observations in the redshift range 2 ≲ z ≲ 7. Results. By performing updated non-equilibrium chemistry calculations we are able to broadly reproduce the latest HI and H2 observations. We find neutral-gas mass density parameters Ωneutral ≃ 10−3 and increasing from lower to higher redshift, in agreement with available HI data. Because of the typically low metallicities during the epoch of reionisation, time-dependent H2 formation is mainly led by the H− channel in self-shielded gas, while H2 grain catalysis becomes important in locally enriched sites at any redshift. Ultraviolet (UV) radiation provides free electrons and facilitates H2 build-up while heating cold metal-poor environments. Resulting H2 fractions can be as high as ∼50% of the cold gas mass at z ∼ 4–8, in line with the latest measurements from high-redshift galaxies. The H2 mass density parameter increases with time until a plateau of ΩH2 ≃ 10−4 is reached. Quantitatively, we find agreement between the derived ΩH2 values and the observations up to z ∼ 7 and both HI and H2 trends are better reproduced by our non-equilibrium H2-based star formation modelling. The predicted gas depletion timescales decrease at lower z in the whole time interval considered, with H2 depletion times remaining below the Hubble time and comparable to the dynamical time at all z. This implies that non-equilibrium molecular cooling is efficient at driving cold-gas collapse in a broad variety of environments and has done so since very early cosmic epochs. While the evolution of chemical species is clearly affected by the details of the UV background and gas self shielding, the assumptions on the adopted initial mass function, different parameterizations of H2 dust grain catalysis, photoelectric heating, and cosmic-ray heating can affect the results in a non-trivial way. In the Appendix, we show detailed analyses of individual processes, as well as simple numerical parameterizations and fits to account for them. Conclusions. Our findings suggest that, in addition to HI, non-equilibrium H2 observations are pivotal probes for assessing cold-gas cosmic abundances and the role of UV background radiation at different epochs.

The Launching of Cold Clouds by Galaxy Outflows. IV. Cosmic-Ray-driven Acceleration

Abstract We carry out a suite of simulations of the evolution of cosmic-ray (CR) driven, radiatively cooled cold clouds embedded in hot material, as found in galactic outflows. In such interactions, CRs stream toward the cloud at the Alfvén speed, which decreases dramatically at the cloud boundary, leading to a bottleneck in which pressure builds up in front of the cloud. At the same time, CRs stream along the sides of the cloud, forming a boundary layer where large filaments develop. Shear in this boundary layer is the primary mode of cloud destruction, which is relatively slow in all cases, but slowest in the cases with the lowest Alfvén speeds. Thus, the CR pressure in the bottleneck region has sufficient time to accelerate the cold clouds efficiently. Furthermore, radiative cooling has relatively little impact on these interactions. Our simulations are two-dimensional and limited by a simplified treatment of CR dynamics, the neglect of CR heating, and an idealized magnetic field geometry. Nevertheless, our results suggest that CRs, when acting as the primary source of momentum input, are capable of accelerating clouds to velocities comparable to those observed in galaxy outflows.

Chaotic cold accretion in giant elliptical galaxies heated by AGN cosmic rays

ABSTRACT Black hole feedback plays a central role in shaping the circumgalactic medium (CGM) of elliptical galaxies. We systematically study the impact of plasma physics on the evolution of ellipticals by performing three-dimensional non-ideal magnetohydrodynamic simulations of the interactions of active galactic nucleus (AGN) jets with the CGM including magnetic fields, and cosmic rays (CRs) and their transport processes. We find that the physics of feedback operating on large galactic scales depends very sensitively on plasma physics operating on small scales. Specifically, we demonstrate that (i) in the purely hydrodynamical case, the AGN jets initially maintain the atmospheres in global thermal balance. However, local thermal instability generically leads to the formation of massive cold discs in the vicinity of the central black hole in disagreement with observations; (ii) including weak magnetic fields prevents the formation of the discs because local B-field amplification in the precipitating cold gas leads to strong magnetic breaking, which quickly extracts angular momentum from the accreting clouds. The magnetic fields transform the cold clouds into narrow filaments that do not fall ballistically; (iii) when plasma composition in the AGN jets is dominated by CRs, and CR transport is neglected, the atmospheres exhibit cooling catastrophes due to inefficient heat transfer from the AGN to CGM despite Coulomb/hadronic CR losses being present; (iv) including CR streaming and heating restores agreement with the observations, i.e. cooling catastrophes are prevented and massive cold central discs do not form. The AGN power is reduced as its energy is utilized efficiently.

Cosmic rays or turbulence can suppress cooling flows (where thermal heating or momentum injection fail)

ABSTRACT The quenching ‘maintenance’ and ‘cooling flow’ problems are important from the Milky Way through massive cluster elliptical galaxies. Previous work has shown that some source of energy beyond that from stars and pure magnetohydrodynamic processes is required, perhaps from active galactic nuclei, but even the qualitative form of this energetic input remains uncertain. Different scenarios include thermal ‘heating’, direct wind or momentum injection, cosmic ray heating or pressure support, or turbulent ‘stirring’ of the intracluster medium (ICM). We investigate these in $10^{12}\!-\!10^{14}\, {\rm M}_{\odot }$ haloes using high-resolution non-cosmological simulations with the FIRE-2 (Feedback In Realistic Environments) stellar feedback model, including simplified toy energy injection models, where we arbitrarily vary the strength, injection scale, and physical form of the energy. We explore which scenarios can quench without violating observational constraints on energetics or ICM gas. We show that turbulent stirring in the central $\sim 100\,$ kpc, or cosmic ray injection, can both maintain a stable low-star formation rate halo for >Gyr time-scales with modest energy input, by providing a non-thermal pressure that stably lowers the core density and cooling rates. In both cases, associated thermal-heating processes are negligible. Turbulent stirring preserves cool-core features while mixing condensed core gas into the hotter halo and is by far the most energy efficient model. Pure thermal heating or nuclear isotropic momentum injection require vastly larger energy, are less efficient in lower mass haloes, easily overheat cores, and require fine tuning to avoid driving unphysical temperature gradients or gas expulsion from the halo centre.

Cosmic Snow Clouds: Self-gravitating Gas Spheres Manifesting Hydrogen Condensation

Abstract We present hydrostatic equilibrium models of spherical, self-gravitating clouds of helium and molecular hydrogen, focusing on the cold, high-density regime where solid- or liquid-hydrogen can form. The resulting structures have masses from 0.1 M ⊙ down to several ×10−8 M ⊙, and span a broad range of radii: 10−4 ≲ R(au) ≲ 107. Our models are fully convective, but all have a two-zone character with the majority of the mass in a small, condensate-free core, surrounded by a colder envelope where phase equilibrium obtains. Convection in the envelope is unusual in that it is driven by a mean-molecular-weight inversion, rather than by an entropy gradient. In fact, the entropy gradient is itself inverted, leading to the surprising result that envelope convection transports heat inward. In turn, that permits the outer layers to maintain steady-state temperatures below the cosmic microwave background. Among our hydrostatic equilibria we identify thermal equilibria appropriate to the Galaxy, in which radiative cooling from H2 is balanced by cosmic-ray heating. These equilibria are all thermally unstable, albeit with very long thermal timescales in some cases. The specific luminosities of all our models are very low, and they therefore describe a type of baryonic dark matter. Consequently such clouds are thermally fragile: when placed in a harsh radiation field, they will be unable to cool effectively and disruption will ensue as heat input drives a secular expansion. Disrupting clouds should leave trails of gas and H2 dust in their wake, which might make them easier to detect. Our models may be relevant to the cometary globules in the Helix Nebula and the G2 cloud orbiting Sgr A*.

Cosmic ray acceleration of cool clouds in the circumgalactic medium

ABSTRACT We investigate a mechanism for accelerating cool (104 K) clouds in the circumgalactic medium (CGM) with cosmic rays (CRs), possibly explaining some characteristics of observed high-velocity clouds (HVCs). Enforcing CRs to stream down their pressure gradient into a region of slow streaming speed results in significant build-up of CR pressure which can accelerate the CGM. We present the results of the first two-dimensional magnetohydrodynamic (MHD) simulations of such ‘CR bottlenecks,’ expanding on simpler simulations in 1D. Although much more investigation is required, we find two main results. First, radiative cooling in the interfaces of these clouds is sufficient to keep the cloud intact to CR wave heating. Secondly, cloud acceleration depends almost linearly with the injected CR flux at low values (comparable to that expected from a Milky Way-like star formation rate), but scales sublinearly at higher CR fluxes in 1D simulations. 2D simulations show hints of sublinear dependence at high CR fluxes but are consistent with pure linear dependence up to the CR fluxes tested. It may therefore be plausible to accelerate cool clouds in the CGM to speeds of hundreds of km s−1.

Starburst and post-starburst high-redshift protogalaxies

Quenching of star-formation has been identified in many starburst and post-starburst galaxies, indicating burst-like star-formation histories (SFH) in the primordial Universe. Galaxies undergoing violent episodes of star-formation are expected to be rich in high energy cosmic rays (CRs). We have investigated the role of these CRs in such environments, particularly how they could contribute to this burst-like SFH via quenching and feedback. These high energy particles interact with the baryon and radiation fields of their host via hadronic processes to produce secondary leptons. The secondary particles then also interact with ambient radiation fields to generate X-rays through inverse-Compton scattering. In addition, they can thermalise directly with the semi-ionised medium via Coulomb processes. Heating at a rate of ∼10−25 erg cm−3 s−1can be attained by Coulomb processes in a star-forming galaxy with one core-collapse SN event per decade, and this is sufficient to cause quenching of star-formation. At high-redshift, a substantial amount of CR secondary electron energy can be diverted into inverse-Compton X-ray emission. This yields an X-ray luminosity of above 1041 erg s−1by redshiftz = 7 which drives a further heating effect, operating over larger scales. This would be able to halt inflowing cold gas filaments, strangulating subsequent star-formation. We selected a sample of 16 starburst and post-starburst galaxies at 7 ≲ z ≲ 9 and determine the star-formation rates they could have sustained. We applied a model with CR injection, propagation and heating to calculate energy deposition rates in these 16 sources. Our calculations show that CR feedback cannot be neglected as it has the strength to suppress star-formation in these systems. We also show that their currently observed quiescence is consistent with the suffocation of cold inflows, probably by a combination of X-ray and CR heating.