Search

Metals in the diffuse, ionized gas at the boundary between the Milky Way's interstellar medium (ISM) and circumgalactic medium (CGM), known as the disk-halo interface (DHI), are valuable tracers of the feedback processes that drive the Galactic fountain. However, metallicity measurements in this region are challenging due to obscuration by the Milky Way ISM and uncertain ionization corrections that affect the total hydrogen column density. In this work, we constrain the ionization corrections to neutral hydrogen column densities using precisely measured electron column densities from the dispersion measure of pulsars that lie in the same globular clusters as UV-bright targets with high-resolution absorption spectroscopy. We address the blending of absorption lines with the ISM by jointly fitting Voigt profiles to all absorption components. We present our metallicity estimates for the DHI of the Milky Way based on detailed photoionization modeling to the absorption from ionized metal lines and ionization-corrected total hydrogen columns. Generally, the gas clouds show a large scatter in metallicity, ranging between $0.04-3.2\ Z_{\odot}$, implying that the DHI consists of a mixture of gaseous structures having multiple origins. We estimate the inflow and outflow timescales of the DHI ionized clouds to be $6 - 35$ Myr. We report the detection of an infalling cloud with super-solar metallicity that suggests a Galactic fountain mechanism, whereas at least one low-metallicity outflowing cloud ($Z < 0.1\ Z_{\odot}$) poses a challenge for Galactic fountain and feedback models., Comment: 28 pages, 15 figures, accepted for publication in ApJ

We investigate the co-evolution of the angular momentum of Milky Way-like galaxies, their circumgalactic gas, and their dark matter halos using zoom-in simulations from the Figuring Out Gas & Galaxies in Enzo (FOGGIE) suite. We examine how the magnitude and orientation of the angular momentum varies over time within the halo and between the components of mass. From z~2 to today, and in general across the simulated halos, the specific angular momenta of the central galaxies and the cool gas in their circumgalactic media (T < 10^5 K) increase together. Over that same period, the specific angular momenta of the hot (>10^6 K) and dark components of the halo change minimally. By z~1, the central galaxies have generally lost association with the angular momentum of their full dark matter halo -- both in magnitude and orientation. We find a wide distribution of angular momentum orientations in the halo, varying by up to 180 degrees over small (~tens of kpc) scales and between the different components of mass. The net angular momenta of the galaxies, their circumgalactic gas, and their dark matter halos are generally misaligned with one another at all cosmic times. The present-day orientation of the central galaxies are established at late times (after z=1), after the rates of cosmic accretion and mergers decline and the disks are able to settle and stabilize their orientation., Comment: 35 pages, 14 figures, submitted to ApJ

We investigate the prospects for detecting and constraining density and temperature inhomogeneities in the circumgalactic medium (CGM) using absorption measurements of metal ions. Distributions in the gas thermal properties could arise from turbulence, gas cooling from the hot phase, and mixing between the cool and hot phases. Focusing on these physically motivated models, we parameterize each with a single parameter for simplicity and provide empirical and theoretical estimates for reasonable parameter values. We then construct the probability distribution functions for each of these scenarios, calculate the effective ion fractions, and fit our models to the COS-Halos absorption measurements to infer the gas densities and metallicities. We find that the models we consider (i) produce similarly good fits to the observations with or without distributions in the gas thermal properties, and (ii) result in detectable changes in the column densities only at the boundaries of reasonable parameter values. We show that He II self-shielding can have a larger effect on the ion fractions than density and temperature fluctuations. As a result, uncertainties in cloud geometry and their spatial distribution, affecting the details of radiation transfer, may obscure the effect of inhomogeneities., Comment: Submitted to ApJ, 11 figures and one table

We present an analysis of \textit{HST}/COS/G160M observations of CIV in the inner circumgalactic medium (CGM) of a novel sample of eight z$\sim$0, L$\approx$L$^{\star}$ galaxies, paired with UV-bright QSOs at impact parameters ($R_\mathrm{proj}$) between 25-130 kpc. The galaxies in this stellar-mass-controlled sample (log$_{10}$M$_{\star}$/M$_{\odot}$ $\sim$ 10.2-10.9 M$_{\odot}$) host super-massive black holes (SMBHs) with dynamically-measured masses spanning log$_{10}$M$_\mathrm{BH}$/M$_{\odot}$ $\sim$ 6.8-8.4; this allows us to compare our results with models of galaxy formation where the integrated feedback history from the SMBH alters the CGM over long timescales. We find that the \ion{C}{IV} column density measurements (N$_{\rm C IV}$) (average log$_{10}$N$_{\rm C IV, CH}$ = 13.94$\pm$0.09 cm$^{-2}$) are largely consistent with existing measurements from other surveys of N$_{\rm C IV}$ in the CGM (average log$_{10}$N$_{\rm C IV, Lit}$ = 13.90$\pm$0.08 cm$^{-2}$), but do not show obvious variation as a function of the SMBH mass. In contrast, specific star-formation rate (sSFR) is highly correlated with the ionized content of the CGM. We find a large spread in sSFR for galaxies with log$_{10}$M$_\mathrm{BH}$/M$_{\odot}$ $>$ 7.0, where the CGM \ion{C}{IV} content shows clear dependence on galaxy sSFR but not M$_\mathrm{BH}$. Our results do not indicate an obvious causal link between CGM CIV and the mass of the galaxy's SMBH; however through comparisons to the EAGLE, Romulus25, $\&$ IllustrisTNG simulations, we find that our sample is likely too small to constrain such causality., Comment: Accepted to ApJ, 28 pages, 7 figures

The cold ($\sim 10^{4}\,{\rm K}$) component of the circumgalactic medium (CGM) accounts for a significant fraction of all galactic baryons. However, using current galaxy-scale simulations to determine the origin and evolution of cold CGM gas poses a significant challenge, since it is computationally infeasible to directly simulate a galactic halo alongside the sub-pc scales that are crucial for understanding the interactions between cold CGM gas and the surrounding ''hot'' medium. In this work, we introduce a new approach: the Cold Gas Subgrid Model (CGSM), which models unresolved cold gas as a second fluid in addition to the standard ''normal'' gas fluid. The CGSM tracks the total mass density and bulk momentum of unresolved cold gas, deriving the properties of its unresolved cloudlets from the resolved gas phase. The interactions between the subgrid cold fluid and the resolved fluid are modeled by prescriptions from high-resolution simulations of ''cloud crushing'' and thermal instability. Through a series of idealized tests, we demonstrate the CGSM's ability to overcome the resolution limitations of traditional hydrodynamics simulations, successfully capturing the correct cold gas mass, its spatial distribution, and the timescales for cloud destruction and growth. We discuss the implications of using this model in cosmological simulations to more accurately represent the microphysics that govern the galactic baryon cycle., Comment: Submitted to MNRAS

The interaction between supermassive black hole (SMBH) feedback and the circumgalactic medium (CGM) continues to be an open question in galaxy evolution. In our study, we use SPH simulations to explore the impact of SMBH feedback on galactic metal retention and the motion of metals and gas into and through the CGM of L$_{*}$ galaxies. We examine 140 galaxies from the 25 Mpc cosmological volume, Romulus25, with stellar masses between 3 $\times$ 10$^{9}$ - 3 $\times$ 10$^{11}$ M$_{\odot}$. We measure the fraction of metals remaining in the ISM and CGM of each galaxy, and calculate the expected mass of its SMBH based on the $M-\sigma$ relation. The deviation of each SMBH from its expected mass, $\Delta M_{BH}$ is compared to the potential of its host via $\sigma$. We find that SMBHs with accreted mass above the empirical $M-\sigma$ relation are about 15\% more effective at removing metals from the ISM than under-massive SMBHs in star forming galaxies. Over-massive SMBHs suppress the overall star formation of their host galaxies and more effectively move metals from the ISM into the CGM. However, we see little evidence for the evacuation of gas from their halos, in contrast with other simulations. Finally, we predict that C IV column densities in the CGM of L$_{*}$ galaxies may depend on host galaxy SMBH mass. Our results show that the scatter in the low mass end of $M-\sigma$ relation may indicate how effective a SMBH is at the local redistribution of mass in its host galaxy., Comment: 16 pages, 14 figures, Submitted to ApJ

We present an analytic model for the cool, $T \approx 10^4$ K, circumgalactic medium (CGM), describing the gas distribution, thermal and ionization state. Our model assumes (total) pressure equilibrium with the ambient warm/hot CGM, photoionization by the metagalactic radiation field, and allows for non-thermal pressure support, parametrized by the ratio of thermal pressures, $\eta = P_{\rm hot,th}/P_{\rm cool,th}$. We apply the model to the COS-Halos data set and find that a nominal model with $\eta = 3$, gas distribution out to $r \approx 0.6 R_{\rm vir}$, and $M_{\rm cool} = 3 \times 10^9~{\rm M_{\odot}}$, corresponding to a volume filling fraction of $f_{\rm V,cool} \approx 1\%$, reproduces the mean measured column densities of HI and low/intermediate metal ions (CII, CIII, SiII, SiIII, MgII). Variation of $\pm 0.5$ dex in the non-thermal pressure or gas mass encompasses $\sim 2/3$ of the scatter between objects. Our nominal model underproduces the measured CIV and SiIV columns, and we show these can be reproduced with (i) a cool phase with $M_{\rm cool} \approx 10^{10}~{\rm M_{\odot}}$ and $\eta \approx 5$, or (ii) an additional component at intermediate temperatures, of cooling or mixing gas, with $M \approx 1.5 \times 10^{10}~{\rm M_{\odot}}$ and occupying $\sim 1/2$ of the total CGM volume. For cool gas with $f_{\rm V,cool} \approx 1\%$ we provide an upper limit on the cloud sizes, $R_{\rm cl} \lesssim 0.5$ kpc. Our results suggest that for the average galaxy CGM, the mass and non-thermal support in the cool phase are lower than estimated in previous works, and extreme scenarios for galactic feedback and non-thermal support may not be necessary. We estimate the rates of cool gas depletion and replenishment, and find accretion onto the galaxy can be entirely offset by condensation, outflows, and IGM accretion, allowing $\dot{M}_{\rm cool}\sim0$ over long timescales., Comment: Submitted to ApJ, 17 pages + appendices, comments welcome

Dwarf galaxies are found to have lost most of their metals via feedback processes; however, there still lacks consistent assessment on the retention rate of metals in their circumgalactic medium (CGM). Here we investigate the metal content in the CGM of 45 isolated dwarf galaxies with $M_*=10^{6.5-9.5}~M_\odot$ ($M_{\rm 200m}=10^{10.0-11.5}~M_\odot$) using {\it HST}/COS. While H I (Ly$\alpha$) is ubiquitously detected ($89\%$) within the CGM, we find low detection rates ($\approx5\%-22\%$) in C II, C IV, Si II, Si III, and Si IV, largely consistent with literature values. Assuming these ions form in the cool ($T\approx10^4$ K) CGM with photoionization equilibrium, the observed H I and metal column density profiles can be best explained by an empirical model with low gas density and high volume filling factor. For a typical galaxy with $M_{\rm 200m}=10^{10.9}~M_\odot$ (median of the sample), our model predicts a cool gas mass of $M_{\rm CGM,cool}\sim10^{8.4}~M_\odot$, corresponding to $\sim2\%$ of the galaxy's baryonic budget. Assuming a metallicity of $0.3Z_\odot$, we estimate that the dwarf galaxy's cool CGM likely harbors $\sim10\%$ of the metals ever produced, with the rest either in more ionized states in the CGM or transported to the intergalactic medium. We further examine the EAGLE simulation and show that H I and low ions may arise from a dense cool medium, while C IV arises from a diffuse warmer medium. Our work provides the community with a uniform dataset on dwarf galaxies' CGM that combines our recent observations, additional archival data and literature compilation, which can be used to test various theoretical models of dwarf galaxies., Comment: Finalized version. Accepted for publication in ApJ

We combine datasets from the CGM$^{2}$ and CASBaH surveys to model a transition point, $R_{\rm cross}$, between circumgalactic and intergalactic media (CGM and IGM, respectively). In total, our data consist of 7244 galaxies at z < 0.5 with precisely measured spectroscopic redshifts, all having impact parameters of 0.01 - 20 comoving Mpc from 28 QSO sightlines with high-resolution UV spectra that cover H I Ly$\alpha$. Our best-fitting model is an exclusionary two-component model that combines a 3D absorber-galaxy cross correlation function with a simple Gaussian profile at inner radii to represent the CGM. By design, this model gives rise to a determination of $R_{\rm cross}$ as a function of galaxy stellar mass, which can be interpreted as the boundary between the CGM and IGM. For galaxies with $10^8 \leq M_{\star}/M_{\odot} \leq 10^{10.5}$, we find that $R_{\rm cross}(M_{\star}) \approx 2 \pm 0.6 R_{\rm vir}$. Additionally, we find excellent agreement between $R_{\rm cross}(M_{\star})$ and the theoretically-determined splashback radius for galaxies in this mass range. Overall, our results favor models of galaxy evolution at z < 0.5 that distribute $T \approx 10^{4}$K gas to distances beyond the virial radius., Comment: Submitted to ApJ, 17 pages, 8 figures

This study addresses how the incidence rate of strong O VI absorbers in a galaxy's circumgalactic medium (CGM) depends on galaxy mass and, independently, on the amount of star formation in the galaxy. We use HST/COS absorption spectroscopy of quasars to measure O VI absorption within 400 projected kpc and 300 km s$^{-1}$ of 52 $M_{*}\sim 10^{10}$ $M_\odot$ galaxies. The galaxies have redshifts $0.12

We use medium-resolution Keck/Echellette Spectrograph and Imager spectroscopy of bright quasars to study cool gas traced by CaII 3934,3969 and NaI 5891,5897 absorption in the interstellar/circumgalactic media of 21 foreground star-forming galaxies at redshifts 0.03 < z < 0.20 with stellar masses 7.4 < log M_*/M_sun < 10.6. The quasar-galaxy pairs were drawn from a unique sample of Sloan Digital Sky Survey quasar spectra with intervening nebular emission, and thus have exceptionally close impact parameters (R_perp < 13 kpc). The strength of this line emission implies that the galaxies' star formation rates (SFRs) span a broad range, with several lying well above the star-forming sequence. We use Voigt profile modeling to derive column densities and component velocities for each absorber, finding that column densities N(CaII) > 10^12.5 cm^-2 (N(NaI) > 10^12.0 cm^-2) occur with an incidence f_C(CaII) = 0.63^+0.10_-0.11 (f_C(NaI) = 0.57^+0.10_-0.11). We find no evidence for a dependence of f_C or the rest-frame equivalent widths W_r(CaII K) or W_r(NaI 5891) on R_perp or M_*. Instead, W_r(CaII K) is correlated with local SFR at >3sigma significance, suggesting that CaII traces star formation-driven outflows. While most of the absorbers have velocities within +/-50 km/s of the host redshift, their velocity widths (characterized by Delta v_90) are universally 30-177 km/s larger than that implied by tilted-ring modeling of the velocities of interstellar material. These kinematics must trace galactic fountain flows and demonstrate that they persist at R_perp > 5 kpc. Finally, we assess the relationship between dust reddening and W_r(CaII K) (W_r(NaI 5891)), finding that 33% (24%) of the absorbers are inconsistent with the best-fit Milky Way E(B-V)-W_r relations at >3sigma significance., Comment: 38 pages, 16 figures, 4 tables. Accepted to ApJ

The circumgalactic medium (CGM) is often assumed to exist in or near hydrostatic equilibrium with the regulation of accretion and the effects of feedback treated as perturbations to a stable balance between gravity and thermal pressure. We investigate global hydrostatic equilibrium in the CGM using four highly-resolved $L^*$ galaxies from the Figuring Out Gas & Galaxies In Enzo (FOGGIE) project. The FOGGIE simulations were specifically targeted at fine spatial and mass resolution in the CGM ($\Delta x \lesssim 1$ kpc $h^{-1}$ and $M \simeq 200M_\odot$). We develop a new analysis framework that calculates the forces provided by thermal pressure gradients, turbulent pressure gradients, ram pressure gradients of large-scale radial bulk flows, centrifugal rotation, and gravity acting on the gas in the CGM. Thermal and turbulent pressure gradients vary strongly on scales of $\lesssim5$ kpc throughout the CGM. Thermal pressure gradients provide the main supporting force only beyond $\sim 0.25R_{200}$, or $\sim50$ kpc at $z=0$. Within $\sim0.25R_{200}$, turbulent pressure gradients and rotational support provide stronger forces than thermal pressure. More generally, we find that global equilibrium models are neither appropriate nor predictive for the small scales probed by absorption line observations of the CGM. Local conditions generally cannot be derived by assuming a global equilibrium, but an emergent global equilibrium balancing radially inward and outward forces is obtained when averaging over the non-equilibrium local conditions on large scales in space and time. Approximate hydrostatic equilibrium holds only at large distances from galaxies even when averaging out small-scale variations., Comment: 36 pages, 24 figures. Accepted to ApJ. See https://foggie.science/cgm_forces.html for animated Figures 6 and 7

The bimodal absorption system imaging campaign (BASIC) aims to characterize the galaxy environments of a sample of 36 HI-selected partial Lyman limit systems (pLLSs) and Lyman limit systems (LLSs) in 23 QSO fields at $z \lesssim 1$. These pLLSs/LLSs provide a unique sample of absorbers with unbiased and well-constrained metallicities, allowing us to explore the origins of metal-rich and low-metallicity circumgalactic medium (CGM) at $z<1$. Here we present Keck/KCWI and VLT/MUSE observations of 11 of these QSO fields (19 pLLSs) that we combine with HST/ACS imaging to identify and characterize the absorber-associated galaxies. We find 23 unique absorber-associated galaxies, with an average of one associated galaxy per absorber. For seven absorbers, all with $<10\%$ solar metallicities, we find no associated galaxies with $\log M_\star \gtrsim 9.0$ within $\rho/R_{vir}$ and $|\Delta v|/v_{esc} \le$ 1.5 with respect to the absorber. We do not find any strong correlations between the metallicities or HI column densities of the gas and most of the galaxy properties, except for the stellar mass of the galaxies: the low-metallicity ([X/H] $\le -1.4$) systems have a probability of $0.39^{+0.16}_{-0.15}$ for having a host galaxy with $\log M_\star \ge 9.0$ within $\rho/R_{vir} \le 1.5$, while the higher metallicity absorbers have a probability of $0.78^{+0.10}_{-0.13}$. This implies metal-enriched pLLSs/LLSs at $z<1$ are typically associated with the CGM of galaxies with $\log M_\star > 9.0$, whereas low-metallicity pLLSs/LLSs are found in more diverse locations, with one population arising in the CGM of galaxies and another more broadly distributed in overdense regions of the universe. Using absorbers not associated with galaxies, we estimate the unweighted geometric mean metallicity of the intergalactic medium to be [X/H] $\lesssim -2.1$ at $z<1$, which is lower than previously estimated., Comment: Key figures: 13 and 18. Modified abstract. Accepted for publication in ApJ. Updated after referee comments, expanded subsection 4.4

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., Comment: 26 pages including 19 figures and appendix. Accepted by ApJ. Updated version includes additional figures, tests, and expanded discussion. Supplemental figures at irynabutsky.me/data

The classical definition of the virial temperature of a galaxy halo excludes a fundamental contribution to the energy partition of the halo: the kinetic energy of non-thermal gas motions. Using simulations of low-redshift, $\sim L^*$ galaxies from the FOGGIE project (Figuring Out Gas & Galaxies In Enzo) that are optimized to resolve low-density gas, we show that the kinetic energy of non-thermal motions is roughly equal to the energy of thermal motions. The simulated FOGGIE halos have $\sim 2\times$ lower bulk temperatures than expected from a classical virial equilibrium, owing to significant non-thermal kinetic energy that is formally excluded from the definition of $T_\mathrm{vir}$. We derive a modified virial temperature explicitly including non-thermal gas motions that provides a more accurate description of gas temperatures for simulated halos in virial equilibrium. Strong bursts of stellar feedback drive the simulated FOGGIE halos out of virial equilibrium, but the halo gas cannot be accurately described by the standard virial temperature even when in virial equilibrium. Compared to the standard virial temperature, the cooler modified virial temperature implies other effects on halo gas: (i) the thermal gas pressure is lower, (ii) radiative cooling is more efficient, (iii) O VI absorbing gas that traces the virial temperature may be prevalent in halos of a higher mass than expected, (iv) gas mass estimates from X-ray surface brightness profiles may be incorrect, and (v) turbulent motions make an important contribution to the energy balance of a galaxy halo., Comment: 30 pages, 14 figures, accepted to ApJ

Quasar absorption systems encode a wealth of information about the abundances, ionization structure, and physical conditions in intergalactic and circumgalactic media. Simple (often single-phase) photoionization models are frequently used to decode such data. Using five discrete absorbers from the COS Absorption Survey of Baryon Harbors (CASBaH) that exhibit a wide range of detected ions (e.g., Mg II, S II--S VI, O II--O VI, Ne VIII), we show several examples where single-phase ionization models cannot reproduce the full set of measured column densities. To explore models that can self-consistently explain the measurements and kinematic alignment of disparate ions, we develop a Bayesian multiphase ionization modeling framework that characterizes discrete phases by their unique physical conditions and also investigates variations in the shape of the UV flux field, metallicity, and relative abundances. Our models require at least two (but favor three) distinct ionization phases ranging from $T \approx 10^{4}$ K photoionized gas to warm-hot phases at $T \lesssim 10^{5.8}$ K. For some ions, an apparently single absorption "component" includes contributions from more than one phase, and up to 30% of the H I is not from the lowest ionization phase. If we assume that all of the phases are photoionized, we cannot find solutions in thermal pressure equilibrium. By introducing hotter, collisionally ionized phases, however, we can achieve balanced pressures. The best models indicate moderate metallicities, often with sub-solar N/$\alpha$, and, in two cases, ionizing flux fields that are softer and brighter than the fiducial Haardt & Madau UV background model., Comment: Accepted for publication in MNRAS

Using 10 sightlines observed with the Hubble Space Telescope/Cosmic Origins Spectrograph, we study the circumgalactic medium (CGM) and outflows of IC1613, which is a low-mass ($M_*\sim10^8~M_\odot$), dwarf irregular galaxy on the outskirts of the Local Group. Among the sightlines, 4 are pointed towards UV-bright stars in IC1613, and the other 6 sightlines are background QSOs at impact parameters from 6 kpc ($<0.1R_{200}$) to 61 kpc ($0.6R_{200}$). We detect a number of Si II, Si III, Si IV, C II, and C IV absorbers, most of which have velocities less than the escape velocity of IC1613 and thus are gravitationally bound. The line strengths of these ion absorbers are consistent with the CGM absorbers detected in dwarf galaxies at low redshifts. Assuming that Si II, Si III, and Si IV comprise nearly 100% of the total silicon, we find 3% ($\sim$8$\times$10$^3~{\rm M_\odot}$), 2% ($\sim$7$\times$10$^3~{\rm M_\odot}$), and 32--42% [$\sim$(1.0--1.3)$\times$10$^5~{\rm M_\odot}$] of the silicon mass in the stars, interstellar medium, and within $0.6R_{200}$ of the CGM of IC1613. We also estimate the metal outflow rate to be ${\rm \dot{M}_{out, Z}\geq1.1\times10^{-5}~M_\odot~yr^{-1}}$ and the instantaneous metal mass loading factor to be $\eta_{\rm Z}\geq0.004$, which are in broad agreement with available observation and simulation values. This work is the first time a dwarf galaxy of such low mass is probed by a number of both QSO and stellar sightlines, and it shows that the CGM of low-mass gas-rich galaxies can be a large reservoir enriched with metals from past and ongoing outflows., Comment: Resubmitted to ApJ with the referee's comments addressed. Minor changes were made to the original draft

From our position embedded within the Milky Way's interstellar medium (ISM), we have limited ability to detect gas at low relative velocities in the extended Galactic halo because those spectral lines are blended with much stronger signals from dense foreground gas. As a result, the content of the Milky Way's circumgalactic medium (CGM) is poorly constrained at $|v_{\rm LSR}|$ $\lesssim$ 150 km s$^{-1}$. To overcome this complication, the QuaStar Survey applies a spectral differencing technique using paired quasar-star sightlines to measure the obscured content of the Milky Way's CGM for the first time. We present measurements of the CIV doublet ($\lambda\lambda$ 1548 \r{A}, 1550 \r{A}), a rest-frame UV metal line transition detected in HST/COS spectra of 30 halo-star/quasar pairs evenly distributed across the sky at Galactic latitudes $|b|>30^\circ$. The 30 halo stars have well-constrained distances (d$\approx$5-14 kpc), and are paired with quasars separated by $<$ 2.8$^\circ$. We argue that the difference in absorption between the quasar and stellar sightlines originates primarily in the Milky Way's extended CGM beyond $\sim$10 kpc. For the Milky Way's extended, low velocity CGM ($|v|<$150 km/s), we place an upper limit on the mean CIV column density of $\rm \Delta logN_{LVCGM} < 13.39$ and find a covering fraction of $f_{\rm CIV,LVCGM} (\rm logN>13.65)=$ 20% [6/30], a value significantly lower than the covering fraction for star-forming galaxies at low redshift. Our results suggest either that the bulk of Milky Way's CIV-traced CGM lies at low Galactic latitudes, or that the Milky Way's CGM is lacking in warm, ionized material compared to low-redshift ($z < 0.1$) star-forming galaxy halos., Comment: 15 pages, 6 figures, 1 table. Accepted to ApJ

Using the N-body+Smoothed particle hydrodynamics code, ChaNGa, we identify two merger-driven processes\textemdash disk disruption and supermassive black hole (SMBH) feedback\textemdash which work together to quench L$^*$ galaxies for over 7 Gyr. Specifically, we examine the cessation of star formation in a simulated Milky Way (MW) analog, driven by an interaction with two minor satellites. Both interactions occur within $\sim$100 Myr of each other, and the satellites both have masses 5 to 20 times smaller than that of their MW-like host galaxy. Using the genetic modification process of \cite{Roth2016}, we generate a set of four zoom-in, MW-mass galaxies all of which exhibit unique star formation histories due to small changes to their assembly histories. In two of these four cases, the galaxy is quenched by $z = 1$. Because these are controlled modifications, we are able to isolate the effects of two closely-spaced minor merger events, the relative timing of which determines whether the MW-mass main galaxy quenches. This one-two punch works to: 1. fuel the primary halo's supermassive black hole (SMBH) at its peak accretion rate; and 2. disrupt the cold, gaseous disk of the host galaxy. The end result is that feedback from the SMBH thoroughly and abruptly ends the galaxy's star formation by $z\approx1$. We search for and find a similar quenching event in {\sc Romulus25}, a hydrodynamical $(25\,\mathrm{Mpc})^3$ volume simulation, demonstrating that the mechanism is common enough to occur even in a small sample of MW-mass quenched galaxies at $z=0$., Comment: 12 pages, 9 figures, accepted to ApJ

We present initial results from the \textit{COS and Gemini Mapping the Circumgalactic Medium} (\mbox{CGMCGM} $\equiv$ CGM$^{2}$) survey. The CGM$^{2}$ survey consists of 1689 galaxies, all with high-quality Gemini GMOS spectra, within 1 Mpc of twenty-two $z \lesssim 1$ quasars, all with S/N$\sim$10 {\emph{HST/COS}} G130M$+$G160M spectra. For 572 of these galaxies having stellar masses $10^{7} M_{\odot} < M_{\star} < 10^{11} M_{\odot}$ and $z \lesssim 0.5$, we show that the \ion{H}{1} covering fraction above a threshold of \NHI$>10^{14} $cm$^{-2}$ is $\gtrsim 0.5$ within 1.5 virial radii ($R_{\rm vir} \sim R_{200m}$). We examine the \ion{H}{1} kinematics and find that the majority of absorption lies within $\pm$ 250 km s$^{-1}$ of the galaxy systemic velocity. We examine \ion{H}{1} covering fractions over a range of impact parameters to infer a characteristic size of the CGM, $R^{14}_{\rm CGM}$, as a function of galaxy mass. $R^{14}_{\rm CGM}$ is the impact parameter at which the probability of observing an absorber with \NHI $>$ 10$^{14}$ cm$^{-2}$ is $>$ 50\%. In this framework, the radial extent of the CGM of $M_{\star} > 10^{9.9} M_{\odot}$ galaxies is $R^{14}_{\rm CGM} = 346^{+57}_{-53}$ kpc or $R^{14}_{\rm CGM} \simeq 1.2R_{\rm vir}$. Intermediate-mass galaxies with $10^{9.2} < M_{\star}/M_{\odot} < 10^{9.9}$ have an extent of $R^{14}_{\rm CGM} = 353^{+64}_{-50}$ kpc or $R^{14}_{\rm CGM} \simeq 2.4R_{\rm vir}$. Low-mass galaxies, $M_{\star} < 10^{9.2} M_{\odot}$, show a smaller physical scale $R^{14}_{\rm CGM} = 177_{-65}^{+70}$ kpc and extend to $R^{14}_{\rm CGM} \simeq 1.6R_{\rm vir}$. Our analysis suggests that using $R_{\rm vir}$ as a proxy for the characteristic radius of the CGM likely underestimates its extent., Comment: 30 pages, 14 figures. Submitted to ApJ, Comments welcome

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., Comment: Accepted for publication in ApJ; 28 pages, 15 figures. Simulation movies can be found at https://www.irynabutsky.me/movies

Project AMIGA (Absorption Maps In the Gas of Andromeda) is a large ultraviolet Hubble Space Telescope program, which has assembled a sample of 43 QSOs that pierce the circumgalactic medium (CGM) of Andromeda (M31) from R=25 to 569 kpc (25 of them probing gas from 25 kpc to about the virial radius-Rvir = 300 kpc-of M31). Our large sample provides an unparalleled look at the physical conditions and distribution of metals in the CGM of a single galaxy using ions that probe a wide range of gas phases (Si II, Si III, Si IV, C II, C IV, and O VI, the latter being from the Far Ultraviolet Spectroscopic Explorer). We find that Si III and O VI have near unity covering factor maintained all the way out to 1.2Rvir and 1.9Rvir, respectively. We show that Si III is the dominant ion over Si II and Si IV at any R. While we do not find that the properties of the CGM of M31 depend strongly on the azimuth, we show that they change remarkably around 0.3-0.5Rvir, conveying that the inner regions of the CGM of M31 are more dynamic and have more complicated multi-phase gas-structures than at R>0.5Rvir. We estimate the metal mass of the CGM within Rvir as probed by Si II, Si III, and Si IV is 2x10^7 Msun and by O VI is >8x10^7 Msun, while the baryon mass of the 10^4-10^5.5 K gas is ~4x10^10 (Z/0.3 Zsun)^(-1) Msun within Rvir. We show that different zoom-in cosmological simulations of L* galaxies better reproduce the column density profile of O VI with R than Si III or the other studied ions. We find that observations of the M31 CGM and zoom-in simulations of L* galaxies have both lower ions showing higher column density dispersion and dependence on R than higher ions, indicating that the higher ionization structures are larger and/or more broadly distributed., Comment: Submitted to the Astrophysical Journal. Comments welcome

We present new calculations of the mass inflow and outflow rates around the Milky Way, derived from a catalog of ultraviolet metal-line high velocity clouds (HVCs). These calculations are conducted by transforming the HVC velocities into the Galactic Standard of Rest (GSR) reference frame, identifying inflowing (v_GSR < 0 km/s) and outflowing (v_GSR > 0 km/s) populations, and using observational constraints on the distance, metallicity, dust content, covering fractions, and total hydrogen column density of each population. After removing HVCs associated with the Magellanic Stream and the Fermi Bubbles, we find inflow and outflow rates in cool (T~10^4 K) ionized gas of dM_in/dt >~ 0.53+/-0.17 (d/12 kpc) (Z/0.2 Z_sun)^-1 M_sun/yr and dM_out/dt >~ 0.16+/-0.06 (d/12 kpc) (Z/0.5 Z_sun)^-1 M_sun/yr. The excess of inflowing over outflowing gas suggests that the Milky Way is currently in an inflow-dominated phase, but the presence of substantial mass flux in both directions supports a Galactic fountain model, in which gas is constantly recycled between the disk and the halo. We also find that the metal flux in both directions (in and out) is indistinguishable. By comparing the outflow rate to the Galactic star formation rate, we present the first estimate of the mass loading factor (etc_HVC) of the disk-wide Milky Way wind, finding eta_HVC >~ 0.10+/-0.06 (d/12 kpc) (Z/0.5 Z_sun)^-1. Including the contributions from low- and intermediate-velocity clouds and from hot gas would increase these inflow and outflow estimates., Comment: Accepted for publication in ApJ, 10 pages, 1 figure, 3 tables

We describe the survey for galaxies in the fields surrounding 9 sightlines to far-UV bright, z~1 quasars that define the COS Absorption Survey of Baryon Harbors (CASBaH) program. The photometry and spectroscopy that comprise the dataset come from a mixture of public surveys (SDSS, DECaLS) and our dedicated efforts on private facilities (Keck, MMT, LBT). We report the redshifts and stellar masses for 5902 galaxies within ~10 comoving-Mpc (cMpc) of the sightlines with a median of z=0.28 and M_* ~ 10^(10.1) Msun. This dataset, publicly available as the CASBaH specDB, forms the basis of several recent and ongoing CASBaH analyses. Here, we perform a clustering analysis of the galaxy sample with itself (auto-correlation) and against the set of OVI absorption systems (cross-correlation) discovered in the CASBaH quasar spectra with column densities N(O^+5) >= 10^(13.5)/cm^2. For each, we describe the measured clustering signal with a power-law correlation function xi(r) = (r/r_0)^(-gamma) and find that (r_0,gamma) = (5.48 +/- 0.07 h_100^-1 Mpc, 1.33 +/- 0.04) for the auto-correlation and (6.00 +/- 1 h^-1 Mpc, 1.25 +/- 0.18) for galaxy-OVI cross-correlation. We further estimate a bias factor of b_gg = 1.3 +/- 0.1 from the galaxy-galaxy auto-correlation indicating the galaxies are hosted by halos with mass M_halo ~ 10^(12.1 +/- 0.05) Msun. Finally, we estimate an OVI-galaxy bias factor b_OVI = 1.0 +/- 0.1 from the cross-correlation which is consistent with OVI absorbers being hosted by dark matter halos with typical mass M_halo ~ 10^(11) Msun. Future works with upcoming datasets (e.g., CGM^2) will improve upon these results and will assess whether any of the detected OVI arises in the intergalactic medium., Comment: 25 pages, 27 figures. Published in ApJS. All of the galaxy spectra are provided as a specDB file (https://github.com/specdb/specdb)

We present a novel absorption line survey using 54 blue horizontal branch stars (BHBs) in the Milky Way halo as background sources for detecting gas flows at the disk-halo interface. Distance measurements to high-latitude ($b$ > 60{\deg}) background stars at 3.1-13.4 kpc, combined with unprecedented spatial sampling and spectral resolution, allow us to examine the 3-dimensional spatial distribution and kinematics of gas flows near the disk. We detect absorption signatures of extraplanar CaII and NaI in Keck HIRES spectra and find that their column densities exhibit no trend with distance to the background sources, indicating that these clouds lie within 3.1 kpc of the disk. We calculate covering fractions of $f_{CaII}$ = 63%, $f_{NaI}$ = 26%, and $f_{HI}$ = 52%, consistent with a picture of the CGM that includes multi-phase clouds containing small clumps of cool gas within hotter, more diffuse gas. Our measurements constrain the scale of any substructure within these cool clouds to <0.5 kpc. CaII and NaI absorption features exhibit an intermediate-velocity (IV) component inflowing at velocities of -75 km/s < v < -25 km/s relative to the local standard of rest, consistent with previously-studied HI structures in this region. We report the new detection of an inflow velocity gradient ${\Delta}v_z$ ~ 6-9 km/s/kpc across the Galactic plane. These findings place constraints on the physical and kinematic properties of CGM gas flows through the disk-halo interface, and support a galactic fountain model in which cold gas rains back onto the disk., Comment: 18 pages, 15 figures, 2 tables, and a 10-page appendix with 9 additional figures. Submitted to ApJ

We address the spatial scale, ionization structure, mass and metal content of gas at the Milky Way disk-halo interface detected as absorption in the foreground of seven closely-spaced, high-latitude halo blue horizontal branch stars (BHBs) with heights z = 3 - 14 kpc. We detect transitions that trace multiple ionization states (e.g. CaII, FeII, SiIV, CIV) with column densities that remain constant with height from the disk, indicating that the gas most likely lies within z < 3.4 kpc. The intermediate ionization state gas traced by CIV and SiIV is strongly correlated over the full range of transverse separations probed by our sightlines, indicating large, coherent structures greater than 1 kpc in size. The low ionization state material traced by CaII and FeII does not exhibit a correlation with either N$_{\rm HI}$ or transverse separation, implying cloudlets or clumpiness on scales less than 10 pc. We find that the observed ratio log(N_SiIV/ N_CIV), with a median value of -0.69+/-0.04, is sensitive to the total carbon content of the ionized gas under the assumption of either photoionization or collisional ionization. The only self-consistent solution for photoionized gas requires that Si be depleted onto dust by 0.35 dex relative to the solar Si/C ratio, similar to the level of Si depletion in DLAs and in the Milky Way ISM. The allowed range of values for the areal mass infall rate of warm, ionized gas at the disk-halo interface is 0.0003 < dM_gas / dtdA [M_sun kpc^-2 yr^-] < 0.006. Our data support a physical scenario in which the Milky Way is fed by complex, multiphase processes at its disk-halo interface that involve kpc-scale ionized envelopes or streams containing pc-scale, cool clumps., Comment: 12 figures, 4 tables, Accepted to ApJ. Revision features areal mass accretion rate calculation

Davies et al. (2019) established that for L^* galaxies the fraction of baryons in the circumgalactic medium (CGM) is inversely correlated with the mass of their central supermassive black holes (BHs) in the EAGLE hydrodynamic simulation. The interpretation is that, over time, a more massive BH has provided more energy to transport baryons beyond the virial radius, which additionally reduces gas accretion and star formation. We continue this research by focusing on the relationship between the 1) BH masses, 2) physical and observational properties of the CGM, and 3) galaxy colours for Milky Way-mass systems. The ratio of the cumulative BH feedback energy over the gaseous halo binding energy is a strong predictor of the CGM gas content, with BHs injecting >~10x the binding energy resulting in gas-poor haloes. Observable tracers of the CGM, including CIV, OVI, and HI absorption line measurements, are found to be effective tracers of the total z~0 CGM halo mass. We use high-cadence simulation outputs to demonstrate that BH feedback pushes baryons beyond the virial radius within 100 Myr timescales, but that CGM metal tracers take longer (0.5-2.5 Gyr) to respond. Secular evolution of galaxies results in blue, star-forming or red, passive populations depending on the cumulative feedback from BHs. The reddest quartile of galaxies with M_*=10^{10.2-10.7} M_solar (median u-r = 2.28) has a CGM mass that is 2.5x lower than the bluest quartile (u-r=1.59). We propose strategies for observing the predicted lower CGM column densities and covering fractions around galaxies hosting more massive BHs using the Cosmic Origins Spectrograph on Hubble., Comment: 15 pages, 10 figures, 1 table, submitted to MNRAS

Quasar absorption-line studies in the ultraviolet (UV) can uniquely probe the nature of the multiphase cool-warm (10^4 < T < 10^6 K) gas in and around galaxy clusters, promising to provide unprecedented insights into 1) interactions between the circumgalactic medium (CGM) associated with infalling galaxies and the hot (T > 10^6 K) X-ray emitting intracluster medium (ICM), 2) the stripping of metal-rich gas from the CGM, and 3) a multiphase structure of the ICM with a wide range of temperatures and metallicities. In this work, we present results from a high-resolution simulation of a ~10^14 solar mass galaxy cluster to study the physical properties and observable signatures of this cool-warm gas in galaxy clusters. We show that the ICM becomes increasingly multiphased at large radii, with the cool-warm gas becoming dominant in cluster outskirts. The diffuse cool-warm gas also exhibits a wider range of metallicity than the hot X-ray emitting gas. We make predictions for the covering fractions of key absorption-line tracers, both in the ICM and in the CGM of cluster galaxies, typically observed with the Cosmic Origins Spectrograph aboard the Hubble Space Telescope (HST). We further extract synthetic spectra to demonstrate the feasibility of detecting and characterizing the thermal, kinematic, and chemical composition of the cool-warm gas using H I, O VI, and C IV lines, and we predict an enhanced population of broad Ly-alpha absorbers tracing the warm gas. Lastly, we discuss future prospects of probing the multiphase structure of the ICM beyond HST., Comment: 15 pages, 12 figures

Galaxies evolve under the influence of gas flows between their interstellar medium and their surrounding gaseous halos known as the circumgalactic medium (CGM). The CGM is a major reservoir of galactic baryons and metals, and plays a key role in the long cycles of accretion, feedback, and recycling of gas that drive star formation. In order to fully understand the physical processes at work within galaxies, it is therefore essential to have a firm understanding of the composition, structure, kinematics, thermodynamics, and evolution of the CGM. In this white paper we outline connections between the CGM and galactic star formation histories, internal kinematics, chemical evolution, quenching, satellite evolution, dark matter halo occupation, and the reionization of the larger-scale intergalactic medium in light of the advances that will be made on these topics in the 2020s. We argue that, in the next decade, fundamental progress on all of these major issues depends critically on improved empirical characterization and theoretical understanding of the CGM. In particular, we discuss how future advances in spatially-resolved CGM observations at high spectral resolution, broader characterization of the CGM across galaxy mass and redshift, and expected breakthroughs in cosmological hydrodynamic simulations will help resolve these major problems in galaxy evolution., Comment: Astro2020 Decadal Science White Paper

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., Comment: 19 pages, 13 figures, submitted to ApJ. Simulation movies available at http://chummels.org/tempest

As galaxies evolve, they must enrich and exchange gas with the surrounding medium, but the timing of these processes and how much gas is involved remain poorly understood. In this work, we leverage metals as tracers of past gas flows to constrain the history of metal ejection and redistribution in M31. This roughly $L*$ galaxy is a unique case where spatially resolved measurements of the gas-phase and stellar metallicity, dust extinction, and neutral ISM gas content are all available, enabling a census of the present metal mass. We combine spatially resolved star formation histories from the Panchromatic Hubble Andromeda Treasury survey with a metal production model to calculate the history of metal production in M31. We find that $1.8\times10^9 \, M_\odot$ of metals, or 62\% of the metal mass formed within $r < 19 \,\mathrm{kpc}$, is missing from the disk in our fiducial model, implying that the M31 disk has experienced significant gaseous outflows over its lifetime. Under a conservative range of model assumptions, we find that between 3\% and 88\% of metals have been lost ($1.9\times10^7 - 6.4\times10^9 \, M_\odot$), so metals are missing even when all model parameters are chosen to favor metal retention. We show that the missing metal mass could be harbored in M31's CGM if the majority of the metals reside in a hot gas phase. Finally, we find that some metal mass produced in the past 1.5 Gyr in the central $\sim5 \,\mathrm{kpc}$ has likely been redistributed to larger radii within the disk., Comment: In press at ApJ; updated to match published version. Figure 7 shows key result

We examine the effects of SMBH feedback on the CGM using a cosmological hydrodynamic simulation \citep[{\sc Romulus25};][]{Tremmel2017} and a set of four zoom-in `genetically modified' Milky Way-mass galaxies sampling different evolutionary paths. By tracing the distribution of metals in the circumgalactic medium (CGM), we show that \ion{O}{6} is a sensitive indicator of supermassive black hole (SMBH) feedback. First, we calculate the column densities of \ion{O}{6} in simulated Milky Way-mass galaxies and compare them with observations from the COS-Halos Survey. Our simulations show column densities of \ion{O}{6} in the CGM consistent with those of COS-Halos star forming and quenched galaxies. These results contrast with those from previous simulation studies which typically underproduce CGM column densities of \ion{O}{6}. We determine that a galaxy's star formation history and assembly record have little effect on the amount of \ion{O}{6} in its CGM. Instead, column densities of \ion{O}{6} are closely tied to galaxy halo mass and BH growth history. The set of zoom-in, genetically modified Milky Way-mass galaxies indicates that the SMBH drives highly metal-enriched material out into its host galaxy's halo which in turn elevates the column densities of \ion{O}{6} in the CGM., Comment: 15 pages, 12 figures, Accepted to ApJ; Note: Figures have transparency that may not render in browser

We survey the highly ionized circumgalactic media (CGM) of 29 blindly selected galaxies at 0.49 < z_(gal) < 1.44 based on high-S/N ultraviolet spectra of z > 1 QSOs and the galaxy database from the COS Absorption Survey of Baryon Harbors (CASBaH). We detect the Ne VIII doublet in nine of the galaxies, and for gas with N(Ne VIII) > 10^13.3 cm^-2 (> 10^13.5 cm^-2), we derive a Ne VIII covering fraction f_c = 75 +15/-25% (44 +22/-20%) within impact parameter (rho) < 200 kpc of M_* = 10^(9.5-11.5) Msol galaxies and f_c = 70 +16/-22% (f_c = 42 +20/-17%) within rho < 1.5 virial radii. We estimate the mass in Ne VIII-traced gas to be M_gas(Ne VIII) > 10^9.5 Msol (Z/Zsol)^-1, or 6-20% of the expected baryonic mass if the Ne VIII absorbers have solar metallicity. Ionizing Ne VII to Ne VIII requires 207 eV, and photons with this energy are scarce in the CGM. However, for the median halo mass and redshift of our sample, the virial temperature is close to the peak temperature for the Ne VIII ion, and the Ne VIII-bearing gas is plausibly collisionally ionized near this temperature. Moreover, we find that photoionized Ne VIII requires cool and low-density clouds that would be highly underpressured (by approximately two orders of magnitude) relative to the putative, ambient virialized medium, complicating scenarios where such clouds could survive. Thus, more complex (e.g., non-equilibrium) models may be required; this first statistical sample of Ne VIII absorber/galaxy systems will provide stringent constraints for future CGM studies., Comment: Published in ApJL, Volume 877, Issue 2, Article L20

We analyze the equivalent widths of HI Ly-$\alpha$ ($W_{Ly\alpha}$) absorption from the inner (R < 160 kpc) circumgalactic medium (CGM) of 85 galaxies at $z \sim 0$ with stellar masses $M*$ ranging $\rm{8 \leq log M* / M_{\odot} \leq 11.6}$. Across three orders of magnitude in stellar mass, the CGM of present-day galaxies exhibits a very high covering fraction of cool hydrogen gas ($f_C = 87\pm 4$\%) indicating that the CGM is ubiquitous in modern, isolated galaxies. When HI Ly-$\alpha$ is detected, its equivalent width declines with increasing radius regardless of the galaxy mass, but the scatter in this trend correlates closely with $M*$. Using the radial and stellar mass correlations, we construct a planar surface describing the cool CGM of modern galaxies: $\log W^{\rm{s}}_{HI 1215} \; = \; (0.34 \pm 0.02) -( 0.0026 \pm 0.0005)\times (R) + (0.286 \pm 0.002) \times \log (M*/M_{\odot})$. The RMS scatter around this bivariate relation is $\sim$0.2 dex. We interpret the explicit correlation between $W_{Ly\alpha}$ and $M*$ to arise from the underlying dark matter halo mass ($M_{halo}$), thereby suggesting a CGM fundamental plane between $W_{Ly\alpha}$, $R$ and $M_{halo}$. This correlation can be used to estimate the underlying dark matter halo mass from observations of saturated HI Ly-$\alpha$ in the CGM of a modern galaxy., Comment: 8 pages, 4 figures, Accepted for publication in ApJ

Radiation in the extreme ultraviolet (EUV) and soft X-ray holds clues to the location of the missing baryons, the energetics in stellar feedback processes, and the cosmic enrichment history. Additionally, EUV and soft X-ray photons help determine the ionization state of most intergalactic and circumgalactic metals, shaping the rate at which cosmic gas cools. Unfortunately, this band is extremely difficult to probe observationally due to absorption from the Galaxy. In this paper, we model the contributions of various sources to the cosmic EUV and soft X-ray backgrounds. We bracket the contribution from (1) quasars, (2) X-ray binaries, (3) hot interstellar gas, (4) circumgalactic gas, (5) virialized gas, and (6) supersoft sources, developing models that extrapolate into these bands using both empirical and theoretical inputs. While quasars are traditionally assumed to dominate these backgrounds, we discuss the substantial uncertainty in their contribution. Furthermore, we find that hot intrahalo gases likely emit an O(1) fraction of this radiation at low redshifts, and that interstellar and circumgalactic emission potentially contribute tens of percent to these backgrounds at all redshifts. We estimate that uncertainties in the angular-averaged background intensity impact the ionization corrections for common circumgalactic and intergalactic metal absorption lines by ~0.3-1 dex, and we show that local emissions are comparable to the cosmic background only at r_prox = 10-100 kpc from Milky Way-like galaxies., Comment: 13 pages, 9 figures, accepted to ApJ. New calculation added in section 2.3

The gas surrounding galaxies outside their disks or interstellar medium and inside their virial radii is known as the circumgalactic medium (CGM). In recent years this component of galaxies has assumed an important role in our understanding of galaxy evolution owing to rapid advances in observational access to this diffuse, nearly invisible material. Observations and simulations of this component of galaxies suggest that it is a multiphase medium characterized by rich dynamics and complex ionization states. The CGM is a source for a galaxy's star-forming fuel, the venue for galactic feedback and recycling, and perhaps the key regulator of the galactic gas supply. We review our evolving knowledge of the CGM with emphasis on its mass, dynamical state, and coevolution with galaxies. Observations from all redshifts and from across the electromagnetic spectrum indicate that CGM gas has a key role in galaxy evolution. We summarize the state of this field and pose unanswered questions for future research., Comment: 46 pages, 11 figures. Authors' draft version; see http://www.annualreviews.org/doi/10.1146/annurev-astro-091916-055240 for a well-edited final version and fun supplementary figures

We explore the circumgalactic metal content traced by commonly observed low ion absorbers, including C II, Si II, Si III, Si IV, and Mg II. We use a set of cosmological hydrodynamical zoom simulations run with the EAGLE model and including a non-equilibrium ionization and cooling module that follows 136 ions. The simulations of z~0.2 L* (M_200=10^11.7-10^12.3 Msol) haloes hosting star-forming galaxies and group-sized (M_200=10^12.7-10^13.3 Msol) haloes hosting mainly passive galaxies reproduce key trends observed by the COS-Halos survey-- low ion column densities show 1) little dependence on galaxy specific star formation rate, 2) a patchy covering fraction indicative of 10^4 K clouds with a small volume filling factor, and 3) a declining covering fraction as impact parameter increases from 20-160 kpc. Simulated Si II, Si III, Si IV, C II, and C III column densities show good agreement with observations, while Mg II is under-predicted. Low ions trace a significant metal reservoir, ~10^8 Msol, residing primarily at 10-100 kpc from star-forming and passive central galaxies. These clouds tend to flow inwards and most will accrete onto the central galaxy within the next several Gyr, while a small fraction are entrained in strong outflows. A two-phase structure describes the inner CGM (<0.5 R_200) with low-ion metal clouds surrounded by a hot, ambient medium. This cool phase is separate from the O VI observed by COS-Halos, which arises from the outer CGM (>0.5 R_200) tracing virial temperature gas around L* galaxies. Physical parameters derived from standard photo-ionization modelling of observed column densities (e.g. aligned Si II/Si III absorbers) are validated against our simulations. Our simulations therefore support previous ionization models indicating that cloud covering factors decline while densities and pressures show little variation with increasing impact parameter., Comment: 26 pages, 19 figures, submitted to MNRAS

Observations reveal massive amounts of OVI around star-forming $L_*$ galaxies, with covering fractions of near unity extending to the host halo's virial radius. This OVI absorption is typically kinematically centered upon photoionized gas, with line widths that are suprathermal and kinematically offset from the galaxy. We discuss various scenarios and whether they could result in the observed phenomenology (cooling gas flows, boundary layers, shocks, virialized gas, photoionized clouds in thermal equilibrium). If predominantly collisionally ionized, as we argue is most probable, the OVI observations require that the circumgalactic medium (CGM) of $L_*$ galaxies holds nearly all the associated baryons within a virial radius ($\sim 10^{11}M_\odot$) and hosts massive flows of cooling gas with $\approx30[nT/30{\rm~cm^{-3}K}]~M_\odot~$yr$^{-1}$, which must be largely prevented from accreting onto the host galaxy. Cooling and feedback energetics considerations require $10 <\langle nT\rangle<100{\rm~cm^{-3}K}$ for the warm and hot halo gases. We argue that virialized gas, boundary layers, hot winds, and shocks are unlikely to directly account for the bulk of the OVI. Furthermore, we show that there is a robust constraint on the number density of many of the photoionized $\sim10^4$K absorption systems that yields upper bounds in the range $n<(0.1-3)\times10^{-3}(Z/0.3)$cm$^{-3}$, where $Z$ is the metallicity, suggestive that the dominant pressure in some photoionized clouds is nonthermal. This constraint, which requires minimal ionization modeling, is in accord with the low densities inferred from more complex photoionization modeling. The large amount of cooling gas that is inferred could re-form these clouds in a fraction of the halo dynamical time, as some arguments require, and it requires much of the feedback energy available from supernovae and stellar winds to be dissipated in the CGM., Comment: 16 pages, matches version accepted to ApJ

We analyze new far-ultraviolet spectra of 13 quasars from the z~0.2 COS-Halos survey that cover the HI Lyman limit of 14 circumgalactic medium (CGM) systems. These data yield precise estimates or more constraining limits than previous COS-Halos measurements on the HI column densities NHI. We then apply a Monte-Carlo Markov Chain approach on 32 systems from COS-Halos to estimate the metallicity of the cool (T~10^4K) CGM gas that gives rise to low-ionization state metal lines, under the assumption of photoionization equilibrium with the extragalactic UV background. The principle results are: (1) the CGM of field L* galaxies exhibits a declining HI surface density with impact parameter Rperp (at >99.5%$ confidence), (2) the transmission of ionizing radiation through CGM gas alone is 70+/-7%; (3) the metallicity distribution function of the cool CGM is unimodal with a median of 1/3 Z_Sun and a 95% interval from ~1/50 Z_Sun to over 3x solar. The incidence of metal poor (<1/100 Z_Sun) gas is low, implying any such gas discovered along quasar sightlines is typically unrelated to L* galaxies; (4) we find an unexpected increase in gas metallicity with declining NHI (at >99.9% confidence) and, therefore, also with increasing Rperp. The high metallicity at large radii implies early enrichment; (5) A non-parametric estimate of the cool CGM gas mass is M_CGM_cool = 9.2 +/- 4.3 10^10 Msun, which together with new mass estimates for the hot CGM may resolve the galactic missing baryons problem. Future analyses of halo gas should focus on the underlying astrophysics governing the CGM, rather than processes that simply expel the medium from the halo., Comment: 26 pages with figures. Accepted to ApJ; see https://github.com/pyigm/pyigm for all published COS-Halos measurements

The total contribution of diffuse halo gas to the galaxy baryon budget strongly depends on its dominant ionization state. In this paper, we address the physical conditions in the highly ionized circumgalactic medium (CGM) traced by OVI absorption lines observed in COS-Halos spectra. We analyze the observed ionic column densities, absorption-line widths and relative velocities, along with the ratios of NV/OVI for 39 fitted Voigt profile components of OVI. We compare these quantities with the predictions given by a wide range of ionization models. Photoionization models that include only extragalactic UV background radiation are ruled out; conservatively, the upper limits to NV/OVI and measurements of N$_{\rm OVI}$ imply unphysically large path lengths > 100 kpc. Furthermore, very broad OVI absorption (b > 40 km s$^{-1}$) is a defining characteristic of the CGM of star-forming L* galaxies. We highlight two possible origins for the bulk of the observed OVI: (1) highly structured gas clouds photoionized primarily by local high energy sources or (2) gas radiatively cooling on large scales behind a supersonic wind. Approximately 20% of circumgalactic OVI does not align with any low-ionization state gas within $\pm$50 km s$^{-1}$ and is found only in halos with M$_{\rm halo}$ < 10$^{12}$ M$_{\odot}$. We suggest that this type of unmatched OVI absorption traces the hot corona itself at a characteristic temperature of 10$^{5.5}$ K. We discuss the implications of these very distinct physical origins for the dynamical state, gas cooling rates, and total baryonic content of L* gaseous halos., Comment: 25 pages, 13 figures. Accepted for publication in ApJ October 6, 2016

We develop a new method to constrain the physical conditions in the cool (~10^4 K) circumgalactic medium (CGM) from measurements of ionic column densities, by assuming that the cool CGM spans a large range of gas densities and that small high-density clouds are hierarchically embedded in large low-density clouds. The new method combines the information available from different sightlines during the photoionization modeling, thus yielding tighter constraints on CGM properties compared to traditional methods which model each sightline individually. Applying this new technique to the COS-Halos survey of low-redshift ~L* galaxies, we find that we can reproduce all observed ion columns in all 44 galaxies in the sample, from the low-ions to OVI, with a single universal density structure for the cool CGM. The gas densities span the range 50 < \rho/\rho_mean < 5x10^5 (\rho_mean is the cosmic mean), while the physical size of individual clouds scales as ~\rho^-1, from ~35 kpc of the low density OVI clouds to ~6 pc of the highest density low-ion clouds. The deduced cloud sizes are too small for this density structure to be driven by self-gravity, thus its physical origin is unclear. The implied cool CGM mass within the virial radius is 1.3x10^10 M_sun (~1% of the halo mass), distributed rather uniformly over the four decades in density. The mean cool gas density profile scales as R^-1.0, where R is the distance from the galaxy center. We construct a 3D model of the cool CGM based on our results, which we argue provides a benchmark for the CGM structure in hydrodynamic simulations. Our results can be tested by measuring the coherence scales of different ions., Comment: Revised version submitted to ApJ following a positive referee report. Main change from previous version is the addition of a discussion sub-section on the origin of OVI. 23 pages, 11 figures

Using Hubble Space Telescope Cosmic Origins Spectrograph observations of 89 QSO sightlines through the Sloan Digital Sky Survey footprint, we study the relationships between C IV absorption systems and the properties of nearby galaxies as well as large-scale environment. To maintain sensitivity to very faint galaxies, we restrict our sample to 0.0015 < z < 0.015, which defines a complete galaxy survey to L > 0.01 L* or stellar mass log M_* > 8 Msun. We report two principal findings. First, for galaxies with impact parameter rho < 1 rvir, C IV detection strongly depends on the luminosity/stellar mass of the nearby galaxy. C IV is preferentially associated with galaxies with log M_* > 9.5 Msun; lower mass galaxies rarely exhibit significant C IV absorption (covering fraction f = 9 +12-6% for 11 galaxies with log M_* < 9.5 Msun). Second, C IV detection within the log M_* > 9.5 Msun population depends on environment. Using a fixed-aperture environmental density metric for galaxies with rho < 160 kpc at z < 0.055, we find that 57+/-12% (8/14) of galaxies in low-density regions (regions with fewer than seven L > 0.15 L* galaxies within 1.5 Mpc) have affiliated C IV absorption; however, none (0/7) of the galaxies in denser regions show C IV. Similarly, the C IV detection rate is lower for galaxies residing in groups with dark-matter halo masses of log Mhalo > 12.5 Msun. In contrast to C IV, H I is pervasive in the CGM without regard to mass or environment. These results indicate that C IV absorbers with log N(C IV) > 13.5 cm^-2 trace the halos of log M_* > 9.5 Msun galaxies but also reflect larger scale environmental conditions., Comment: 26 pages, 13 figures. ApJ, in press

We analyze the mass, temperature, metal enrichment, and OVI abundance of the circumgalactic medium (CGM) around $z\sim 0.2$ galaxies of mass $10^9 M_\odot

To investigate the evolution of metal-enriched gas over recent cosmic epochs as well as to characterize the diffuse, ionized, metal-enriched circumgalactic medium (CGM), we have conducted a blind survey for C IV absorption systems in 89 QSO sightlines observed with the Hubble Space Telescope (HST) Cosmic Origins Spectrograph (COS). We have identified 42 absorbers at z < 0.16, comprising the largest uniform blind sample size to date in this redshift range. Our measurements indicate an increasing C IV absorber number density per comoving path length (dN/dX = 7.5 +/- 1.1) and modestly increasing mass density relative to the critical density of the Universe (Omega(C IV) = 10.0 +/- 1.5 x 10^-8 ) from z ~ 1.5 to the present epoch, consistent with predictions from cosmological hydrodynamical simulations. Furthermore, the data support a functional form for the column density distribution function that deviates from a single power-law, also consistent with independent theoretical predictions. As the data also probe heavy element ions in addition to C IV at the same redshifts, we identify, measure, and search for correlations between column densities of these species where components appear aligned in velocity. Among these ion-ion correlations, we find evidence for tight correlations between C II and Si II, C II and Si III, and C IV and Si IV, suggesting that these pairs of species arise in similar ionization conditions. However, the evidence for correlations decreases as the difference in ionization potential increases. Finally, when controlling for observational bias, we find only marginal evidence for a correlation (86.8% likelihood) between the Doppler line width b(C IV) and column density N(C IV)., Comment: 32 pages, 18 figures, 4 tables. Resubmitted to ApJ after referee review

Modern analyses of structure formation predict a universe tangled in a 'cosmic web' of dark matter and diffuse baryons. These theories further predict that at low-z, a significant fraction of the baryons will be shock-heated to $T \sim 10^{5}-10^{7}$K yielding a warm-hot intergalactic medium (WHIM), but whose actual existence has eluded a firm observational confirmation. We present a novel experiment to detect the WHIM, by targeting the putative filaments connecting galaxy clusters. We use HST/COS to observe a remarkable QSO sightline that passes within $\Delta d = 3$ Mpc from the 7 inter-cluster axes connecting 7 independent cluster-pairs at redshifts $0.1 \le z \le 0.5$. We find tentative excesses of total HI, narrow HI (NLA; Doppler parameters $b<50$ km/s), broad HI (BLA; $b \ge 50$ km/s) and OVI absorption lines within rest-frame velocities of $\Delta v \lesssim 1000$ km/s from the cluster-pairs redshifts, corresponding to $\sim 2$, $\sim 1.7$, $\sim 6$ and $\sim 4$ times their field expectations, respectively. Although the excess of OVI likely comes from gas close to individual galaxies, we conclude that most of the excesses of NLAs and BLAs are truly intergalactic. We find that the covering fractions, $f_c$, of BLAs close to cluster-pairs are $\sim 4-7$ times higher than the random expectation (at the $\sim 2 \sigma$ c.l.), whereas the $f_c$ of NLAs and OVI are not significantly enhanced. We argue that a larger relative excess of BLAs compared to those of NLAs close to cluster-pairs may be a signature of the WHIM in inter-cluster filaments. By extending the present analysis to tens of sightlines our experiment offers a promising route to detect the WHIM., Comment: MNRAS in press; 21 pages, 10 figures, 5 tables and 15-page appendix with additional 5 figures and 4 tables

We analyze the low-redshift (z~0.2) circumgalactic medium by comparing absorption-line data from the COS-Halos Survey to absorption around a matched galaxy sample from two cosmological hydrodynamic simulations. The models include different prescriptions for galactic outflows, namely hybrid energy/momentum driven wind (ezw), and constant winds (cw). We extract for comparison direct observables including equivalent widths, covering factors, ion ratios, and kinematics. Both wind models are generally in good agreement with these observations for HI and certain low ionization metal lines, but show poorer agreement with higher ionization metal lines including SiIII and OVI that are well-observed by COS-Halos. These discrepancies suggest that both wind models predict too much cool, metal-enriched gas and not enough hot gas, and/or that the metals are not sufficiently well-mixed. This may reflect our model assumption of ejecting outflows as cool and unmixing gas. Our ezw simulation includes a heuristic prescription to quench massive galaxies by super-heating its ISM gas, which we show yields sufficient low ionisation absorption to be broadly consistent with observations, but also substantial OVI absorption that is inconsistent with data, suggesting that gas around quenched galaxies in the real Universe does not cool. At impact parameters of <50 kpc, recycling winds dominate the absorption of low ions and even HI, while OVI almost always arises from metals ejected longer than 1 Gyr ago. The similarity between the wind models is surprising, since we show that they differ substantially in their predicted amount and phase distribution of halo gas. We show that this similarity owes mainly to our comparison here at fixed stellar mass rather than at fixed halo mass in our previous works, which suggests that CGM properties are more closely tied to the stellar mass of galaxies rather than halo mass., Comment: 21 pages, 11 figures, submitted to MNRAS

KA1858+4850 is a narrow-line Seyfert 1 galaxy at redshift 0.078 and is among the brightest active galaxies monitored by the Kepler mission. We have carried out a reverberation mapping campaign designed to measure the broad-line region size and estimate the mass of the black hole in this galaxy. We obtained 74 epochs of spectroscopic data using the Kast Spectrograph at the Lick 3-m telescope from February to November of 2012, and obtained complementary V-band images from five other ground-based telescopes. We measured the H-beta light curve lag with respect to the V-band continuum light curve using both cross-correlation techniques (CCF) and continuum light curve variability modeling with the JAVELIN method, and found rest-frame lags of lag_CCF = 13.53 (+2.03, -2.32) days and lag_JAVELIN = 13.15 (+1.08, -1.00) days. The H-beta root-mean-square line profile has a width of sigma_line = 770 +/- 49 km/s. Combining these two results and assuming a virial scale factor of f = 5.13, we obtained a virial estimate of M_BH = 8.06 (+1.59, -1.72) x 10^6 M_sun for the mass of the central black hole and an Eddington ratio of L/L_Edd ~ 0.2. We also obtained consistent but slightly shorter emission-line lags with respect to the Kepler light curve. Thanks to the Kepler mission, the light curve of KA1858+4850 has among the highest cadences and signal-to-noise ratios ever measured for an active galactic nucleus; thus, our black hole mass measurement will serve as a reference point for relations between black hole mass and continuum variability characteristics in active galactic nuclei.