No Arabic abstract
Large reservoirs of cold (~ 10^4 K) gas exist out to and beyond the virial radius in the circumgalactic medium (CGM) of all types of galaxies. Photoionization modeling suggests that cold CGM gas has significantly lower densities than expected by theoretical predictions based on thermal pressure equilibrium with hot CGM gas. In this work, we investigate the impact of cosmic ray physics on the formation of cold gas via thermal instability. We use idealized three-dimensional magnetohydrodynamic simulations to follow the evolution of thermally unstable gas in a gravitationally stratified medium. We find that cosmic ray pressure lowers the density and increases the size of cold gas clouds formed through thermal instability. We develop a simple model for how the cold cloud sizes and the relative densities of cold and hot gas depend on cosmic ray pressure. Cosmic ray pressure can help counteract gravity to keep cold gas in the CGM for longer, thereby increasing the predicted cold mass fraction and decreasing the predicted cold gas inflow rates. Efficient cosmic ray transport, by streaming or diffusion, redistributes cosmic ray pressure from the cold gas to the background medium, resulting in cold gas properties that are in-between those predicted by simulations with inefficient transport and simulations without cosmic rays. We show that cosmic rays can significantly reduce galactic accretion rates and resolve the tension between theoretical models and observational constraints on the properties of cold CGM gas.
We use hydrodynamical simulations of two Milky Way-mass galaxies to demonstrate the impact of cosmic-ray pressure on the kinematics of cool and warm circumgalactic gas. Consistent with previous studies, we find that cosmic-ray pressure can dominate over thermal pressure in the inner 50 kpc of the circumgalactic medium (CGM), creating an overall cooler CGM than that of similar galaxy simulations run without cosmic rays. We generate synthetic sightlines of the simulated galaxies CGM and use Voigt profile fitting methods to extract ion column densities, Doppler-b parameters, and velocity centroids of individual absorbers. We directly compare these synthetic spectral line fits with HST/COS CGM absorption-line data analyses, which tend to show that metallic species with a wide range of ionization potential energies are often kinematically aligned. Compared to the Milky-Way simulation run without cosmic rays, the presence of cosmic-ray pressure in the inner CGM creates narrower OVI absorption features and broader SiIII absorption features, a quality which is more consistent with observational data. Additionally, because the cool gas is buoyant due to nonthermal cosmic-ray pressure support, the velocity centroids of both cool and warm gas tend to align in the simulated Milky Way with feedback from cosmic rays. Our study demonstrates that detailed, direct comparisons between simulations and observations, focused on gas kinematics, have the potential to reveal the dominant physical mechanisms that shape the CGM.
Traditional cosmological hydrodynamics simulations fail to spatially resolve the circumgalatic medium (CGM), the reservoir of tenuous gas surrounding a galaxy and extending to its virial radius. We introduce the technique of Enhanced Halo Resolution (EHR), enabling more realistic physical modeling of the simulated CGM by consistently forcing gas refinement to smaller scales throughout the virial halo of a simulated galaxy. We investigate the effects of EHR in the Tempest simulations, a suite of Enzo-based cosmological zoom simulations following the evolution of an L* galaxy, resolving spatial scales of 500 comoving pc out to 100 comoving kpc in galactocentric radius. Among its many effects, EHR (1) changes the thermal balance of the CGM, increasing its cool gas content and decreasing its warm/hot gas content; (2) preserves cool gas structures for longer periods; and (3) enables these cool clouds to exist at progressively smaller size scales. Observationally, this results in a boost in low ions like H I and a drop in high ions like O VI throughout the CGM. These effects of EHR do not converge in the Tempest simulations, but extrapolating these trends suggests that the CGM in reality is a mist consisting of ubiquitous, small, long-lived, cool clouds suspended in a hot medium at the virial temperature of the halo. Additionally, we explore the physical mechanisms to explain why EHR produces the above effects, proposing that it works both by (1) better sampling the distribution of CGM phases enabling runaway cooling in the denser, cooler tail of the phase distribution; and (2) preventing cool gas clouds from artificially mixing with the ambient hot halo and evaporating. Evidence is found for both EHR mechanisms occurring in the Tempest simulations.
Cold ($Tsim 10^{4} mathrm{K}$) gas is very commonly found in both galactic and cluster halos. There is no clear consensus on its origin. Such gas could be uplifted from the central galaxy by galactic or AGN winds. Alternatively, it could form in situ by thermal instability. Fragmentation into a multi-phase medium has previously been shown in hydrodynamic simulations to take place once $t_mathrm{cool}/t_mathrm{ff}$, the ratio of the cooling time to the free-fall time, falls below a threshold value. Here, we use 3D plane-parallel MHD simulations to investigate the influence of magnetic fields. We find that because magnetic tension suppresses buoyant oscillations of condensing gas, it destabilizes all scales below $l_mathrm{A}^mathrm{cool} sim v_mathrm{A} t_mathrm{cool}$, enhancing thermal instability. This effect is surprisingly independent of magnetic field orientation or cooling curve shape, and sets in even at very low magnetic field strengths. Magnetic fields critically modify both the amplitude and morphology of thermal instability, with $delta rho/rho propto beta^{-1/2}$, where $beta$ is the ratio of thermal to magnetic pressure. In galactic halos, magnetic fields can render gas throughout the entire halo thermally unstable, and may be an attractive explanation for the ubiquity of cold gas, even in the halos of passive, quenched galaxies.
Thermal instability (TI) can strongly affect the structure and dynamics of the interstellar medium (ISM) in the Milky Way and other disk galaxies. Thermal conduction plays an important role in the TI by stabilizing small scales and limiting the size of the smallest condensates. In the magnetized ISM, however, heat is conducted anisotropically (primarily along magnetic field lines). We investigate the effects of anisotropic thermal conduction on the nonlinear regime of the TI by performing two-dimensional magnetohydrodynamic simulations. We present models with magnetic fields of different initial geometries and strengths, and compare them to hydrodynamic models with isotropic conduction. We find anisotropic conduction does not significantly alter the overall density and temperature statistics in the saturated state of the TI. However, it can strongly affect the shapes and sizes of cold clouds formed by the TI. For example, for uniform initial fields long filaments of cold gas are produced that are reminiscent of some observed HI clouds. For initially tangled fields, such filaments are not produced. We also show that anisotropic conduction suppresses turbulence generated by evaporative flows from the surfaces of cold blobs, which may have implications for mechanisms for driving turbulence in the ISM.
The majority of galactic baryons reside outside of the galactic disk in the diffuse gas known as the circumgalactic medium (CGM). While state-of-the art simulations excel at reproducing galactic disk properties, many struggle to drive strong galactic winds or to match the observed ionization structure of the CGM using only thermal supernova feedback. To remedy this, recent studies have invoked non-thermal cosmic ray (CR) stellar feedback prescriptions. However, numerical schemes of CR transport are still poorly constrained. We explore how the choice of CR transport affects the multiphase structure of the simulated CGM. We implement anisotropic CR physics in the astrophysical simulation code, {sc Enzo} and simulate a suite of isolated disk galaxies with varying prescriptions for CR transport: isotropic diffusion, anisotropic diffusion, and streaming. We find that all three transport mechanisms result in strong, metal-rich outflows but differ in the temperature and ionization structure of their CGM. Isotropic diffusion results in a spatially uniform, warm CGM that underpredicts the column densities of low-ions. Anisotropic diffusion develops a reservoir of cool gas that extends further from the galactic center, but disperses rapidly with distance. CR streaming projects cool gas out to radii of 200 kpc, supporting a truly multiphase medium. In addition, we find that streaming is less sensitive to changes in constant parameter values like the CR injection fraction, transport velocity, and resolution than diffusion. We conclude that CR streaming is a more robust implementation of CR transport and motivate the need for detailed parameter studies of CR transport.