No Arabic abstract
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.
We investigate the impact of cosmic rays (CRs) on the circumgalactic medium (CGM) in FIRE-2 simulations, for ultra-faint dwarf through Milky Way (MW)-mass halos hosting star-forming (SF) galaxies. Our CR treatment includes injection by supernovae, anisotropic streaming and diffusion along magnetic field lines, collisional and streaming losses, with constant parallel diffusivity $kappasim3times10^{29},mathrm{cm^2 s^{-1}}$ chosen to match $gamma$-ray observations. With this, CRs become more important at larger halo masses and lower redshifts, and dominate the pressure in the CGM in MW-mass halos at $zlesssim 1-2$. The gas in these CR-dominated halos differs significantly from runs without CRs: the gas is primarily cool (a few $sim10^{4},$K), and the cool phase is volume-filling and has a thermal pressure below that needed for virial or local thermal pressure balance. Ionization of the low and mid ions in this diffuse cool gas is dominated by photo-ionization, with O VI columns $gtrsim 10^{14.5},mathrm{cm^{-2}}$ at distances $gtrsim 150,mathrm{kpc}$. CR and thermal gas pressure are locally anti-correlated, maintaining total pressure balance, and the CGM gas density profile is determined by the balance of CR pressure gradients and gravity. Neglecting CRs, the same halos are primarily warm/hot ($Tgtrsim 10^{5},$K) with thermal pressure balancing gravity, collisional ionization dominates, O VI columns are lower and Ne VIII higher, and the cool phase is confined to dense filaments in local thermal pressure equilibrium with the hot phase.
We use adaptive mesh refinement cosmological simulations to study the spatial distribution and covering fraction of OVI absorption in the circumgalactic medium (CGM) as a function of projected virial radius and azimuthal angle. We compare these simulations to an observed sample of 53 galaxies from the Multiphase Galaxy Halos Survey. Using Mockspec, an absorption line analysis pipeline, we generate synthetic quasar absorption line observations of the simulated CGM. To best emulate observations, we studied the averaged properties of 15,000 mock samples each of 53 sightlines having a distribution of $D/R_{vir}$ and sightline orientation statistically consistent with the observations. We find that the OVI covering fraction obtained for the simulated galaxies agrees well with the observed value for the inner halo ($D/R_{vir} leq 0.375$) and is within $1.1sigma$ in the outer halo ($D/R_{vir} > 0.75$), but is underproduced within $0.375 < D/R_{vir} leq 0.75$. The observed bimodal distribution of OVI covering fraction with azimuthal angle, showing higher frequency of absorption along the projected major and minor axes of galaxies, is not reproduced in the simulations. Further analysis reveals the spatial-kinematic distribution of OVI absorbing gas is dominated by outflows in the inner halo mixed with a inflowing gas that originates from further out in the halo. Though the CGM of the individual simulated galaxies exhibit spatial structure, the flat azimuthal distribution occurs because the individual simulated galaxies do not develop a CGM structure that is universal from galaxy to galaxy.
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.
We analyse the properties of circumgalactic gas around simulated galaxies in the redshift range z >= 3, utilising a new sample of cosmological zoom simulations. These simulations are intended to be representative of the observed samples of Lyman-alpha emitters recently obtained with the MUSE instrument (halo masses ~10^10-10^11 solar masses). We show that supernova feedback has a significant impact on both the inflowing and outflowing circumgalactic medium by driving outflows, reducing diffuse inflow rates, and by increasing the neutral fraction of inflowing gas. By temporally stacking simulation outputs we find that significant net mass exchange occurs between inflowing and outflowing phases: none of the phases are mass-conserving. In particular, we find that the mass in neutral outflowing hydrogen declines exponentially with radius as gas flows outwards from the halo centre. This is likely caused by a combination of both fountain-like cycling processes and gradual photo/collisional ionization of outflowing gas. Our simulations do not predict the presence of fast-moving neutral outflows in the CGM. Neutral outflows instead move with modest radial velocities (~ 50 kms^-1), and the majority of the kinetic energy is associated with tangential rather than radial motion.
We present an analysis of the role of feedback in shaping the neutral hydrogen (HI) content of simulated disc galaxies. For our analysis, we have used two realisations of two separate Milky Way-like (~L*) discs - one employing a conservative feedback scheme (MUGS), the other significantly more energetic (MaGICC). To quantify the impact of these schemes, we generate zeroth moment (surface density) maps of the inferred HI distribution; construct power spectra associated with the underlying structure of the simulated cold ISM, in addition to their radial surface density and velocity dispersion profiles. Our results are compared with a parallel, self-consistent, analysis of empirical data from THINGS (The HI Nearby Galaxy Survey). Single power-law fits (P~k^gamma) to the power spectra of the stronger-feedback (MaGICC) runs (over spatial scales corresponding to 0.5 kpc to 20 kpc) result in slopes consistent with those seen in the THINGS sample (gamma = -2.5). The weaker-feedback (MUGS) runs exhibit shallower power law slopes (gamma = -1.2). The power spectra of the MaGICC simulations are more consistent though with a two-component fit, with a flatter distribution of power on larger scales (i.e., gamma = -1.4 for scales in excess of 2 kpc) and a steeper slope on scales below 1 kpc (gamma = -5), qualitatively consistent with empirical claims, as well as our earlier work on dwarf discs. The radial HI surface density profiles of the MaGICC discs show a clear exponential behaviour, while those of the MUGS suite are essentially flat; both behaviours are encountered in nature, although the THINGS sample is more consistent with our stronger (MaGICC) feedback runs.