Search

High-velocity outflows are ubiquitous in compact, massive ( M _* ∼ 10 ^11 M _⊙ ), z ∼ 0.5 galaxies with extreme star formation surface densities (Σ _SFR ∼ 2000 M _⊙ yr ^−1 kpc ^−2 ). We have previously detected and characterized these outflows using Mg ii absorption lines. To probe their full extent, we present Keck/KCWI integral field spectroscopy of the [O ii ] and Mg ii emission nebulae surrounding all of the 12 galaxies in this study. We find that [O ii ] is more effective than Mg ii in tracing low surface brightness, extended emission in these galaxies. The [O ii ] nebulae are spatially extended beyond the stars, with radial extent R _90 between 10 and 40 kpc. The nebulae exhibit nongravitational motions, indicating galactic outflows with maximum blueshifted velocities ranging from −335 to −1920 km s ^−1 . The outflow kinematics correlate with the bursty star formation histories of these galaxies. Galaxies with the most recent bursts of star formation (within the last

We analyze the O vi content and kinematics for 126 H i -selected absorbers at 0.14 ≲ z ≲ 0.73 for which the metallicities of their cool photoionized phase have been determined. We separate the absorbers into 100 strong Ly α forest systems (SLFSs with 15 ≲ log N (H i ) < 16.2) and 26 partial Lyman Limit systems (pLLSs with 16.2 ≤ log N (H i ) ≤ 17.2). The sample is drawn from the COS CGM Compendium (CCC) and has O vi coverage in signal-to-noise ratio ≥ 8 Hubble Space Telescope/COS G130M/G160M QSO spectra, yielding a 2 σ completeness level of $\mathrm{log}N$ (O vi ) ≥ 13.6. The O vi detection rates differ substantially between low-metallicity (LM; [X/H] ≤ −1.4) and high-metallicity (HM; [X/H] > −1.4) SLFSs, with 20% and 60% detection rates, respectively. The O vi detection frequency for the HM and LM pLLSs is, however, similar at ∼60%. The SLFSs and pLLSs without detected O vi are consistent with the absorbing gas being in a single phase, while those with O vi trace multiphase gas. We show that the O vi velocity widths and column densities have different distributions in LM and HM gas. We find a strong correlation between O vi column density and metallicity. The strongest ( $\mathrm{log}N$ (O vi ) ≳ 14) and broadest O vi absorbers are nearly always associated with HM absorbers, while weaker O vi absorbers are found in both LM and HM absorbers. From comparisons with galaxy-selected and blind O vi surveys, we conclude absorbers with $\mathrm{log}N$ (O vi ) ≳ 14 most likely arise in the CGM of star-forming galaxies. Absorbers with weak O vi likely trace the extended CGM or intergalactic medium, while those without O vi likely originate in the intergalactic medium.

We report on the analysis of a multiphase Lyman limit system (LLS) at z = 0.39047 identified toward the background quasar FBQS J0209–0438. The O VI doublet lines associated with this absorber have a different profile from the low-ionization metals and H i . Ly α has a very broad H i ( b ≈ 150 km s ^−1 ) component well-aligned with one of the O VI components. The Doppler b -parameters for the broad H i and O VI indicate gas with T = (0.8 − 2.0) × 10 ^6 K and a total hydrogen column density that is an order of magnitude larger than the cooler phase of gas responsible for the LLS. Observations by the Very Large Telescope MUSE show two moderately star-forming galaxies within ρ ≲ 105 kpc and ∣Δ v ∣ ≲ 130 km s ^−1 of the absorber, one of them a dwarf galaxy ( M _* ≈ 10 ^6 M _⊙ ) overlapping with the quasar point-spread function, and the other a larger galaxy ( R _1/2 ≈ 4 kpc) with M _* ≈ 3 × 10 ^10 M _⊙ and M _h ≈ 7 × 10 ^11 M _⊙ , and the dwarf galaxy within its virial radius. Although the absorption is aligned with the extended major axis of the larger galaxy, the line-of-sight velocity of the absorbing gas is inconsistent with corotating accretion. The metallicity inferred for the LLS is lower than the gas phase [O/H] of the two galaxies. The mixture of cool and warm/hot gas phases for the absorbing gas and its proximity and orientation to the galaxy pair points to the LLS being a high-velocity gas in the combined halo environment of both galaxies.

The recent detection of large column density absorption lines from highly ionized gas in a few directions through the circumgalactic medium (CGM) of the Milky Way (MW) has been puzzling. The inferred temperature from these absorption lines far exceeds the virial temperature of the MW, and the column densities are also too large to be easily explained. In this paper, we propose a novel idea to explain these observations and claim that they may not have originated from the CGM, but from a totally different type of source, namely, stellar ejecta from supernovae (SNe) above the Galactic disk that happen to lie in the line of sight to the background quasars. About ∼20% of massive OB stars (progenitors of core-collapse supernovae) are known to be runaway stars that have high ejection velocities near the Galactic plane and can end up exploding as SNe above the Galactic disk. We show that the associated reverse shock in the supernova remnant in the early nonradiative phase can heat the ejecta to temperatures of ≳10 ^7 K and can naturally explain the observed high column density of ions in the observed “supervirial” phase along with α -enriched supersolar abundance that is typical of core-collapse supernovae. However, SNe from runaway stars has a covering fraction of ≲0.7% and thus can only explain the observations along limited sightlines.

The presence of the ≈10 ^6 K gas in the circumgalactic medium of the Milky Way (MW) has been well established. However, the location and the origin of the newly discovered hot gas at “supervirial (SV)” temperatures of ≈10 ^7 K have been puzzling. This hot gas has been detected in both absorption and emission; here, we focus on the emitting gas only. We show that both the “virial” and the SV temperature gas, as observed in emission, occupy disk-like extraplanar regions, in addition to the diffuse virial temperature gas filling the halo of the MW. We perform idealized hydrodynamical simulations to show that the ≈10 ^7 K emitting gas is likely to be produced by stellar feedback in and around the Galactic disk. We further show that the emitting gas at both SV and virial temperatures in the extraplanar regions is metal enriched and is not in hydrostatic equilibrium with the halo but is continuously evolving.

Feedback-driven massive outflows play a crucial role in galaxy evolution by regulating star formation and influencing the dynamics of surrounding media. Extracting outflow properties from spectral lines is a notoriously difficult process for a number of reasons, including the possibility that a substantial fraction of the outflow is carried by dense gas in a very narrow range in velocity. This gas can hide in spectra with insufficient resolution. Empirically motivated analysis based on the apparent optical depth method, commonly used in the literature, neglects the contribution of this gas, and may therefore underestimate the true gas column density. More complex semianalytical line transfer (e.g., SALT) models, on the other hand, allow for the presence of this gas by modeling the radial density and velocity of the outflows as power laws. Here we compare the two approaches to quantify the uncertainties in the inferences of outflow properties based on 1D “down-the-barrel” spectra, using the UV spectra of the CLASSY galaxy sample. We find that empirical modeling may significantly underestimate the column densities relative to SALT analysis, particularly in the optically thick regime. We use simulations to show that the main reason for this discrepancy is the presence of a large amount of dense material at low velocities, which can be hidden by the finite spectral resolution of the data. The SALT models in turn could overestimate the column densities if the assumed power laws of the density profiles are not a property of actual outflows.

We analyze measurements of the thermal Sunyaev–Zeldovich (tSZ) effect arising in the circumgalactic medium (CGM) of L * galaxies, reported by J. N. Bregman et al. (B+22) and S. Das et al. (D+23). In our analysis, we use the Y. Faerman et al. CGM models, a new power-law model (PLM), and the TNG100 simulation. For a given M _vir , our PLM has four parameters: the fraction, f _hCGM , of the halo baryon mass in hot CGM gas, the ratio, ϕ _T , of the actual gas temperature at the virial radius to the virial temperature, and the power-law indices, a _P _,th and a _n for the thermal electron pressure and the hydrogen nucleon density. The B+22 Compton- y profile implies steep electron pressure slopes ( a _P _,th ≃ 2). For isothermal conditions, the temperature is at least 1.1 × 10 ^6 K, with a hot CGM gas mass of up to 3.5 × 10 ^11 M _⊙ for a virial mass of 2.75 × 10 ^12 M _⊙ . However, if isothermal, the gas must be expanding out of the halos. An isentropic equation of state is favored for which hydrostatic equilibrium is possible. The B+22 and D+23 results are consistent with each other and with recent (0.5–2 keV) CGM X-ray observations of Milky Way mass systems. For M _vir ≃ 3 × 10 ^12 M _⊙ , the scaled Compton pressure integrals, $E{(z)}^{-2/3}{Y}_{500}/{M}_{\mathrm{vir},12}^{5/3}$ , lie in the narrow range, 2.5 × 10 ^−4 –5.0 × 10 ^−4 kpc ^2 , for all three sets of observations. TNG100 underpredicts the tSZ parameters by factors ∼0.5 dex for the L * galaxies, suggesting that the feedback strengths and CGM gas losses are overestimated in the simulated halos at these mass scales.

We present initial results associating galaxies in the Multi-Unit Spectroscopic Explorer (MUSE) Ultra Deep Field (MUDF) with gas seen in absorption along the line of sight to two bright quasars in this field to explore the dependence of metals in the circumgalactic medium (CGM) on galaxy properties. The MUDF includes ∼140 hr of Very Large Telescope (VLT)/MUSE data and 90 orbits of Hubble Space Telescope (HST)/G141M grism observations alongside VLT/Ultraviolet and Visual Echelle Spectrograph spectroscopy of the two quasars and several bands of HST imaging. We compare the metal absorption around galaxies in this field as a function of impact parameter, azimuthal angle, and galaxy metallicity across redshifts 0.5 < z < 3.2. Due to the depth of our data and a large field of view, our analysis extends to low stellar masses (600 kpc). We find a correlation between the absorber equivalent width and the number of nearby galaxies, but do not detect a significant anticorrelation with the impact parameter. Our full sample does not show any significant change in absorber incidence as a function of azimuthal angle. However, we do find a bimodality in the azimuthal angle distribution of absorption at small impact parameters (

Circumgalactic dust grains trace the circulation of mass and metals between star-forming regions and gaseous galactic halos, giving insight into feedback and tidal stripping processes. We perform a search for ultraviolet (UV) reflection nebulae produced by extraplanar dust around 551 nearby ( D < 100 Mpc), edge-on disk galaxies using archival near-ultraviolet and far-ultraviolet images from the Galaxy Evolution Explorer (GALEX), accounting for the point-spread function (FWHM = 4″–5″). We detect extraplanar emission ubiquitously in stacks of galaxies binned by morphology and star formation rate, with scale heights of h _h = 1–2.3 kpc and ≈10% of the total (reddened) flux in the galaxy found beyond the B -band isophotal level of μ _B = 25 mag arcsec ^−2 . This emission is detected in 7% of the individual galaxies, and an additional one-third have at least 5% of their total flux found beyond μ _B = 25 mag arcsec ^−2 in a disk component. The extraplanar luminosities and colors are consistent with reflection nebulae rather than stellar halos and indicate that, on average, disk galaxies have an extraplanar dust mass of 5%–15% of that in their interstellar medium. This suggests that recycled material composes at least a third of the inner circumgalactic medium ( R < 10 kpc) in ∼ L * galaxies.

Mapping stars and gas in nearby galaxies is fundamental for understanding their growth and the impact of their environment. This issue is addressed by comparing the stellar “edges” of galaxies D _stellar , defined as the outermost diameter where in situ star formation significantly drops, with the gaseous distribution parameterized by the neutral atomic hydrogen diameter measured at 1 M _⊙ pc ^−2 , D _HI . By sampling a broad H i mass range 10 ^5 M _⊙ < M _HI < 10 ^11 M _⊙ , we find several dwarf galaxies with M _HI < 10 ^9 M _⊙ from the field and Fornax Cluster that are distinguished by D _stellar ≫ D _HI . For the cluster dwarfs, the average H i surface density near D _stellar is ∼0.3 M _⊙ pc ^−2 , reflecting the impact of quenching and outside-in gas removal from ram pressure and tidal interactions. In comparison, D _stellar / D _HI ranges between 0.5 and 2 in dwarf field galaxies, consistent with the expectations from stellar feedback. Only more massive disk galaxies in the field can thus be characterized by the common assumption that D _stellar ≲ D _HI . We discover a break in the D _stellar – M _⋆ relation at m _break ∼ 4 × 10 ^8 M _⊙ that potentially differentiates the low-mass regime, where the influence of stellar feedback and environmental processes more prominently regulates the sizes of nearby galaxies. Our results highlight the importance of combining deep optical and H i imaging for understanding galaxy evolution.

We explore the influence of radio-mode feedback on the properties of the cool circumgalactic medium (CGM). To this end, we assemble a statistical sample of approximately 30,000 radio galaxies with background quasars by combining optical spectroscopic measurements of luminous red galaxies and quasars from the year 1 data set of the Dark Energy Spectroscopic Instrument and radio sources from the LOw-Frequency ARray Two-metre Sky Survey (LoTSS) DR2 catalog and the Very Large Array Sky Survey (VLASS) quick-look catalog. Galaxies with similar optical properties but with no radio counterparts in LoTSS and VLASS are selected as the control group. We measure the cool CGM properties of radio galaxies and their control samples traced by Mg ii absorption lines, including covering fraction, rest equivalent width, and gas kinematics. Our results show no significant difference in the properties of gas around radio galaxies and their control sample, indicating that the operating radio-mode feedback of massive galaxies does not produce detectable effects on the properties of the cool CGM. Finally, we show that the CGM of radio galaxies contains a nonnegligible amount of cool gas with approximately 10 ^10 M _⊙ . This abundance can place a stringent constraint on the radio-mode feedback models.

We have performed wide-field, ultra-low-surface-brightness H α emission-line mapping around NGC 1068 with the newly commissioned Circumgalactic H α Spectrograph. NGC 1068 is notable for its active galactic nucleus, which globally ionizes gas in the disk and halo. Line-emitting diffuse ionized gas is distributed throughout the galactic disk and large-scale ionized filaments are found well beyond the disk, aligned with the cone angle of the central jet. We report the discovery of a new ribbon of ionized gas around NGC 1068 beyond even the known outer filamentary structure, located 20 kpc from the galaxy. The H α surface brightness of this ribbon is on the order of the bright telluric lines, ranging from 4 to 16 R , with fainter regions on the order of the sky background continuum. Unlike previous extended emission, the ribbon is not as well aligned with the current axis of the central jet. It is not associated with any galactic structure or known tidal features in the halo of NGC 1068, though it may originate from a larger distribution of unmapped neutral atomic or molecular gas in the halo. The morphology of the ribbon emission in H α is correlated with extended UV emission around NGC 1068. H α -to-UV flux ratios in the ribbon are comparable to extended emission-line ratios in the halos of NGC 5128, NGC 253, and M82. The H α excess in the ribbon gas suggests ionization by slow shocks or a mixture of in situ star formation and photoionization and collisional ionization processes.

Using stacked emission-line flux measurements of cool circumgalactic gas (CGM) in lower-mass galaxies (10 ^9.0 ≤ M _* / M _⊙ ≤ 10 ^10.2 ), we measure the dependence of the emission characteristics on orientation relative to the disk plane as a function of radius and compare to what we found previously for massive ( M _* > 10 ^10.4 M _⊙ ) early-type galaxies. Although the line ratios (the lower [N ii ]/H α and [O iii ]/H β ) suggest an overall softer ionizing source than in the more massive galaxies, consistent with previous findings, we find the same ionization hardening signature (a higher [N ii ]/H α ratio in the inner polar region) along the polar direction at small radii that we found for the more massive galaxies. The line ratio in the inner polar bin is distinct from that measured for the inner planar bin with 99.99%, confidence and with >99.9% confidence we conclude that it lies outside the star formation regime. The effective hardening of the ionization of the CGM along the polar axis, at small radii, could indicate either relic effects of active galactic nucleus activity or shock ionization. In either case, this signature appears to be ubiquitous across the stellar mass range we are able to explore with our spectral stacking technique and currently available archival data.

Observations of metallic doublet emission lines, particularly Mg ii λ λ 2796, 2803, provide crucial information for understanding galaxies and their circumgalactic medium. This study explores the effects of resonant scattering on the Mg ii doublet lines and the stellar continuum in spherical and cylindrical geometries. Our findings show that under certain circumstances, resonance scattering can cause an increase in the doublet flux ratio and the escaping flux of the lines beyond what is expected in optically thin spherical media. As expected, the doublet ratio is consistently lower than the intrinsic ratio when the scattering medium is spherically symmetric and dusty. However, if the scattering medium has a disk shape, such as face-on disk galaxies, and is viewed face-on, the doublet ratio is predicted to be higher than 2. It is also shown that doublet ratios as low as those observed in compact star-forming galaxies cannot be explained solely by pure dust attenuation of intrinsic Mg ii emission lines in spherical models unless dust opacity deviates markedly from that expected based on the dust-to-Mg ^+ gas ratio of our Galaxy. The importance of the continuum-pumped emission lines and expanding media is discussed to understand observational aspects, including doublet flux ratios, which can be lower than 1.5 or higher than 2, as well as symmetric or asymmetric line profiles. It is also discussed that the diffuse warm neutral medium may be an important source of Mg ii emission. These results provide insight into the complexity of the shape and orientation of distant, spatially unresolved galaxies.

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

The circumgalactic medium (CGM) of star-forming dwarf galaxies plays a key role in regulating the galactic baryonic cycle. We investigate how susceptible the CGM of dwarf satellite galaxies is to ram pressure stripping in Milky Way–like environments. In a suite of hydrodynamical wind tunnel simulations, we model an intermediate-mass dwarf satellite galaxy ( M _* = 10 ^7.2 M _⊙ ) with a multiphase interstellar medium (ISM; M _ISM = 10 ^7.9 M _⊙ ) and CGM ( M _CGM,vir = 10 ^8.5 M _⊙ ) along two first-infall orbits to more than 500 Myr past pericenter of a Milky Way–like host. The spatial resolution is ∼79 pc in the star-forming ISM and 316−632 pc in the CGM. Our simulations show that the dwarf satellite CGM removal is fast and effective: more than 95% of the CGM mass is ram pressure stripped within a few hundred megayears, even under a weak ram pressure orbit where the ISM stripping is negligible. The conditions for CGM survival are consistent with the analytical halo gas stripping predictions in McCarthy et al. We also find that including the satellite CGM does not effectively shield its galaxy, and therefore the ISM stripping rate is unaffected. Our results imply that a dwarf galaxy CGM is unlikely to be detected in satellite galaxies; and that the star formation of gaseous dwarf satellites is likely devoid of replenishment from a CGM.

The origin of our Galaxy’s high-velocity clouds (HVCs) remains a mystery after many decades of effort. In this paper, we use the TNG50 simulation of the IllustrisTNG project to identify cool, dense clouds that match observations of Galactic H i HVCs. We track these clouds back in time to determine their origin. For a TNG50 Milky Way−like galaxy, we find that only 17% of HVCs can be tracked directly to the disk and 21% to material stripped out of satellites. The majority of HVCs (62%) arise from warm and hot circumgalactic gas that cools through thermal instability. They then obtain their anomalous velocities through interactions with the turbulent circumgalactic medium. At TNG50 resolution, we do not see evidence for HVCs forming out of very low metallicity intergalactic material. Instead, low-metallicity HVCs are most likely associated with satellites. These results suggest that Galactic HVCs are highly heterogeneous in their origin and can provide insight into the physical processes that shape the circumgalactic medium, such as disk outflows, satellite accretion, and thermal instabilities.

We present a suite of high-resolution numerical simulations to study the evolution and survival of dust in hot galactic winds. We implement a novel dust framework in the Cholla hydrodynamics code and use wind tunnel simulations of cool, dusty clouds to understand how thermal sputtering affects the dust content of galactic winds. Our simulations illustrate how various regimes of cloud evolution impact dust survival, dependent on cloud size, wind properties, and dust grain size. We find that significant amounts of dust can survive in winds in all scenarios, even without shielding from the cool phase of outflows. We present an analytic framework that explains this result, along with an analysis of the impact of cloud evolution on the total fraction of dust survival. Using these results, we estimate that 60% of 0.1 μ m dust that enters a starburst-driven wind could survive to populate both the hot and cool phases of the halo, based on a simulated distribution of cloud properties. We also investigate how these conclusions depend on grain size, exploring grains from 0.1 μ m to 10 Å. Under most circumstances, grains smaller than 0.01 μ m cannot withstand hot-phase exposure, suggesting that the small grains observed in the circumgalactic medium (CGM) are either formed in situ due to the shattering of larger grains, or must be carried there in the cool phase of outflows. Finally, we show that the dust-to-gas ratio of clouds declines as a function of distance from the galaxy due to cloud–wind mixing and condensation. These results provide an explanation for the vast amounts of dust observed in the CGMs of galaxies and beyond.

The dispersion measure of fast radio bursts (FRBs), arising from the interactions with free electrons along the propagation path, constitutes a unique probe of the cosmic baryon distribution. Their constraining power is further enhanced in combination with observations of the foreground large-scale structure and intervening galaxies. In this work, we present the first constraints on the partition of the cosmic baryons between the intergalactic medium (IGM) and circumgalactic medium (CGM), inferred from the FLIMFLAM spectroscopic survey. In its first data release, the FLIMFLAM survey targeted galaxies in the foreground of eight localized FRBs. Using Bayesian techniques, we reconstruct the underlying ∼Mpc-scale matter density field that is traced by the IGM gas. Simultaneously, deeper spectroscopy of intervening foreground galaxies (at impact parameters b _⊥ ≲ r _200 ) and the FRB host galaxies constrains the contribution from the CGM. Applying Bayesian parameter inference to our data and assuming a fiducial set of priors, we infer the IGM cosmic baryon fraction to be ${f}_{\mathrm{igm}}={0.59}_{-0.10}^{+0.11}$ and a CGM gas fraction of ${f}_{\mathrm{gas}}={0.55}_{-0.29}^{+0.26}$ for 10 ^10 M _⊙ ≲ M _halo ≲ 10 ^13 M _⊙ halos. The mean FRB host dispersion measure (rest-frame) in our sample is $\langle {\mathrm{DM}}_{\mathrm{host}}\rangle ={90}_{-19}^{+29}\,\mathrm{pc}\ \ {\mathrm{cm}}^{-3}$ , of which $\langle {\mathrm{DM}}_{\mathrm{host}}^{\mathrm{unk}}\rangle ={69}_{-19}^{+28}\,\mathrm{pc}\ \ {\mathrm{cm}}^{-3}$ arises from the host galaxy interstellar medium (ISM) and/or the FRB progenitor environment. While our current f _igm and f _gas uncertainties are too broad to constrain most galactic feedback models, this result marks the first measurement of the IGM and CGM baryon fractions, as well as the first systematic separation of the FRB host dispersion measure into two components: arising from the halo and from the inner ISM/FRB engine.

The properties of warm-hot gas around ∼ L _* galaxies can be studied with absorption lines from highly ionized metals. We predict Ne viii column densities from cosmological zoom-in simulations of halos with masses in ∼10 ^12 and ∼10 ^13 M _☉ from the Feedback in Realistic Environments (FIRE) project. Ne viii traces the volume-filling, virial-temperature gas in ∼10 ^12 M _☉ halos. In ∼10 ^13 M _☉ halos the Ne viii gas is clumpier, and biased toward the cooler part of the warm-hot phase. We compare the simulations to observations from the COS Absorption Survey of Baryon Harbors (or CASBaH) and COS Ultraviolet Baryon Survey (or CUBS). We show that when inferring halo masses from stellar masses to compare simulated and observed halos, it is important to account for the scatter in the stellar-mass–halo-mass relation, especially at M _⋆ ≳ 10 ^10.5 M _☉ . Median Ne viii columns in the fiducial FIRE-2 model are about as high as observed upper limits allow, while the simulations analyzed do not reproduce the highest observed columns. This suggests that the median Ne viii profiles predicted by the simulations are consistent with observations, but that the simulations may underpredict the scatter. We find similar agreement with analytical models that assume a product of the halo gas fraction and metallicity (relative to solar) ∼0.1, indicating that observations are consistent with plausible circumgalactic medium temperatures, metallicities, and gas masses. Variants of the FIRE simulations with a modified supernova feedback model and/or active galactic nuclei feedback included (as well as some other cosmological simulations from the literature) more systematically underpredict Ne viii columns. The circumgalactic Ne viii observations therefore provide valuable constraints on simulations that otherwise predict realistic galaxy properties.

We present CloudFlex , an open-source tool for predicting absorption-line signatures of cool gas in galaxy halos with small-scale structure. Motivated by analyses of ∼10 ^4 K material in hydrodynamical simulations of turbulent, multiphase media, we model cool gas structures as complexes of cloudlets sampled from a power-law distribution of mass $\propto {m}_{\mathrm{cl}}^{-\alpha }$ with velocities drawn from a turbulent velocity field. The user may specify α , the lower limit of the cloudlet mass distribution ( ${m}_{\mathrm{cl},\min }$ ), and several other parameters that set the mass, size, and velocity distribution of the complex. This permits investigation of the relation between these parameters and absorption-line observables. As a proof-of-concept, we calculate the Mg ii λ 2796 absorption induced by the cloudlets in background quasi-stellar object (QSO) spectra. We demonstrate that, at fixed metallicity, the covering fraction of sight lines with equivalent widths W _2796 < 0.3 Å increases significantly with decreasing ${m}_{\mathrm{cl},\min }$ , cloudlet number density ( n _cl ), and complex size. We then use this framework to predict the halo-scale W _2796 distribution around ∼ L ^* galaxies. We show that the observed incidences of W _2796 > 0.3 Å sight lines with impact parameters 10 kpc < R _⊥ < 50 kpc in projected QSO–galaxy studies are consistent with our model over much of parameter space. However, they are underpredicted by models with ${m}_{\mathrm{cl},\min }\geqslant 100{M}_{\odot }$ and n _cl ≥ 0.03 cm ^−3 , in keeping with a picture in which the inner cool circumgalactic medium (CGM) is dominated by numerous low-mass cloudlets ( m _cl ≲ 100 M _⊙ ) with a volume filling factor ≲1%. When used to model absorption-line data sets built from multi-sight line and/or spatially extended background probes, CloudFlex enables detailed constraints on the size and velocity distributions of structures comprising the photoionized CGM.

We present Keck Cosmic Web Imager integral field observations of extended Ly α emission in the circumgalactic medium of 27 typical star-forming galaxies at z ∼ 2, drawn from the Multi-Object Spectrometer for Infra-Red Exploration (MOSFIRE) Deep Evolution Field (MOSDEF) survey. Using composite spectra in two bins of star formation rate (SFR), star formation rate surface density (Σ _SFR ), and other galactic properties, we measure spatial variations in the Ly α profile across three regions in the Ly α halo. We find single-peaked, redshifted profiles are ubiquitous within a central 7 kpc radius region. Further out in the halo (7–14 and 14–21 kpc), the Ly α profile of the resonantly scattered emission exhibits more diversity, either transitioning to a double-peaked profile or remaining single peaked across the halo. We find a shorter scale length of the Ly α halo surface brightness profile for composite halos with faster winds. The composites have a similar average inclination, suggesting those with faster winds clear channels in the interstellar medium (ISM), reducing the fraction of Ly α photons resonantly scattered to large radii. A uniform expanding shell radiative transfer model reproduces the shape but not the normalization of the observed double-peaked Ly α profiles. Models that adopt a more realistic, clumpy ISM are likely needed to reproduce both the shape and normalization of the Ly α profiles.

Jetted active galactic nuclei (AGN) hosting extended photoionized nebulae provide us with a unique view of the timescales associated with AGN activity. Here, we present a new green bean galaxy at z = 0.304458 ± 0.000007 with large-scale jet-induced radio emission. The spectral energy distributions of the radio components show steep spectral indices ( α = −0.85 to −0.92 for the extended regions, and α = −1.02 for the faint radio core), and spectral age modeling of the extended radio emission indicates that the lobes are >6 Myr old. It is unclear whether the jet is active, or is a remnant with an off-time of 2–3 Myr. Several detached clouds lie around the host galaxy up to 37.8 kpc away from the nucleus, and their ionization profile indicates a decline (∼2 dex) in the AGN ionizing photon production over the past ∼0.15 Myr. Furthermore, we measure a blueshift for one of the clouds that is spatially coincident with the path of the radio jet. The cloud is likely illuminated by the photoionizing AGN, and potentially underwent an interaction with the relativistic jet. Our multiwavelength analysis suggests that RGB1 was in a phase of jet production prior to the radiatively efficient accretion phase traced by the detached cloud emission. It is unclear whether RGB1 transitioned into a low-excitation radio galaxy or an inactive galaxy over the past ∼0.15 Myr, or whether the extended radio and optical emission trace distinct accretion phases that occurred in succession.

We present an analysis of Hubble Space Telescope COS/G160M observations of C IV in the inner circumgalactic medium (CGM) of a novel sample of eight z ∼ 0, L ≈ L ^⋆ galaxies, paired with UV-bright QSOs at impact parameters ( R _proj ) between 25 and 130 kpc. The galaxies in this stellar-mass-controlled sample (log _10 M _⋆ / M _⊙ ∼ 10.2–10.9 M _⊙ ) host supermassive black holes (SMBHs) with dynamically measured masses spanning log _10 M _BH / M _⊙ ∼ 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 C IV column density measurements ( N _C IV ; average log _10 N _C IV,CH = 13.94 ± 0.09 cm ^−2 ) are largely consistent with existing measurements from other surveys of N _C IV in the CGM (average log _10 N _C IV,Lit = 13.90 ± 0.08 cm ^−2 ), but do not show obvious variation as a function of the SMBH mass. By 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 _BH / M _⊙ > 7.0, where the CGM C IV content shows a clear dependence on galaxy sSFR but not M _BH . Our results do not indicate an obvious causal link between CGM C IV and the mass of the galaxy’s SMBH; however, through comparisons to the EAGLE, Romulus25, and IllustrisTNG simulations, we find that our sample is likely too small to constrain such causality.

We present a sample of 305 QSO candidates having ∣ b ∣ < 30°, the majority with GALEX magnitudes near-UV < 18.75. To generate this sample, we apply UV–IR color selection criteria to photometric data from the Ultraviolet Galactic Plane Survey as part of GALEX-CAUSE, the Million Quasars Catalog, Gaia DR2, and Pan-STARRS DR1. 165 of these 305 candidate UV-bright active galactic nuclei (AGN; 54%) have published spectroscopic redshifts from 45 different surveys, confirming them as AGN. We further obtained low-dispersion, optical, long-slit spectra with the Apache Point Observatory 3.5 m, MDM 2.4 m, and MDM 1.3 m telescopes for 84 of the candidates, and confirm 86% ( N = 72) as AGN, generally with z < 0.6. Of these 72 confirmed AGN, 25 are newly discovered low-latitude QSOs without any previous spectroscopy. These sources fill a gap in the Galactic latitude coverage of the available samples of known UV-bright QSO background probes. Along with a description of the confirmed QSO properties, we provide the fully reduced, flux- and wavelength-calibrated spectra of 72 low-latitude QSOs through the Mikulski Archive for Space Telescopes. Future Hubble Space Telescope/Cosmic Origins Spectrograph spectroscopy of these low-Galactic-latitude QSOs has the potential to transform our view of the Milky Way and Local Group circumgalactic medium.

NGC 777 provides an example of a phenomenon observed in some group-central ellipticals, in which the temperature profile shows a central peak, despite the short central cooling time of the intragroup medium. We use deep Chandra X-ray observations of the galaxy, supported by uGMRT 400 MHz radio imaging, to investigate the origin of this hot core. We confirm the centrally peaked temperature profile and find that the entropy and cooling time both monotonically decline to low values (2.62 ${}_{-0.18}^{+0.19}$ keV cm ^2 and ${71.3}_{-13.1}^{+12.8}$ Myr, respectively) in the central ∼700 pc. Faint diffuse radio emission surrounds the nuclear point source, with no clear jets or lobes but extending to ∼10 kpc on the northwest–southeast axis. This alignment and extent agree well with a previously identified filamentary H α + [N ii ] nebula. While cavities are not firmly detected, we see X-ray surface brightness decrements on the same axis at 10–20 kpc radii, which are consistent with the intragroup medium having been pushed aside by expanding radio lobes. Any such outburst must have occurred long enough ago for lobe emission to have faded below detectability. Cavities on this scale would be capable of balancing radiative cooling for at least ∼240 Myr. We consider possible causes of the centrally peaked temperature profile, including gravitational heating of gas as the halo relaxes after a period of active galactic nucleus jet activity, and heating by particles leaking from the remnant relativistic plasma of the old radio jets.

We present an extension of the massively parallel, GPU native, astrophysical hydrodynamics code Cholla to magnetohydrodynamics (MHD). Cholla solves the ideal MHD equations in their Eulerian form on a static Cartesian mesh utilizing the Van Leer + constrained transport integrator, the HLLD Riemann solver, and reconstruction methods at second and third order. Cholla’s MHD module can perform ≈260 million cell updates per GPU-second on an NVIDIA A100 while using the HLLD Riemann solver and second order reconstruction. The inherently parallel nature of GPUs combined with increased memory in new hardware allows Cholla’s MHD module to perform simulations with resolutions ∼500 ^3 cells on a single high-end GPU (e.g., an NVIDIA A100 with 80 GB of memory). We employ GPU direct Message Passing Interface to attain excellent weak scaling on the exascale supercomputer Frontier, while using 74,088 GPUs and simulating a total grid size of over 7.2 trillion cells. A suite of test problems highlights the accuracy of Cholla’s MHD module and demonstrates that zero magnetic divergence in solutions is maintained to round off error. We also present new testing and CI tools using GoogleTest, GitHub Actions, and Jenkins that have made development more robust and accurate and ensure reliability in the future.

We present statistics of z ∼ 6–9 galaxy outflows indicated by spatially extended gas emission and broad lines. With a total of 61 spectroscopically confirmed galaxies at z ∼ 6–9 in the JWST CEERS, GLASS, and ERO data, we find four galaxies with [O iii ] + H β ionized gas emission significantly extended beyond the kiloparsec-scale stellar components on the basis of the emission line images constructed by the subtraction of NIRCam broadband (line on/off-band) images. By comparison with low- z galaxies, the fraction of galaxies with spatially extended gas, 4/18, at z ∼ 6–9 is an order of magnitude higher than that at z ∼ 0–1, which can be explained by events triggered by frequent major mergers at high redshift. We also investigate medium- and high-resolution NIRSpec spectra of 30 galaxies at z ∼ 6–9, and identify five galaxies with broad (140−800 km s ^−1 ) lines in the [O iii ] forbidden line emission, suggestive of galaxy outflows. One galaxy at z = 6.38 shows both spatially extended gas emission and the broad lines, while none of the galaxies with spatially extended gas emission or broad lines present a clear signature of active galactic nuclei (AGN) in either the line diagnostics or type 1 AGN line broadening (>1000 km s ^−1 ), which suggests that outflows are mainly driven by stellar feedback. The existence of galaxies with/without spatially extended gas emission or broad lines may be explained by different viewing angles toward outflows or the fact that these are galaxies in the early, late, or post phases of galaxy outflows at high redshift, where relatively large fractions of such galaxies indicate longer-duration and/or more frequent outflows at early cosmic epochs.

The Astro2020 Decadal Survey has identified the mapping of the circumgalactic medium (CGM; the gaseous plasma around galaxies) as a key objective. We explore the prospects for characterizing the CGM in and around nearby galaxy halos with a future large-grasp X-ray microcalorimeter. We create realistic mock observations from hydrodynamical simulations (EAGLE, IllustrisTNG, and Simba) that demonstrate a wide range of potential measurements, which will address the open questions in galaxy formation and evolution. By including all background and foreground components in our mock observations, we show why it is impossible to perform these measurements with current instruments, such as X-ray CCDs, and why only microcalorimeters will allow us to distinguish the faint CGM emission from the bright Milky Way (MW) foreground emission lines. We find that individual halos of MW mass can, on average and depending on star formation rate, be traced out to large radii, around R _500 , and for larger galaxies even out to R _200 , using prominent emission lines, such as O vii , or O viii . Furthermore, we show that emission-line ratios for individual halos can reveal the radial temperature structure. Substructure measurements show that it will be possible to relate azimuthal variations to the feedback mode of the galaxy. We demonstrate the ability to construct temperature, velocity, and abundance ratio maps from spectral fitting for individual galaxy halos, which reveal rotation features, active galactic nucleus outbursts, and enrichment.

The GAMA J0913−0107 system is a rare conjunction of a submillimeter galaxy (SMG) at z ≈ 2.7 and two background QSOs with projected separations

We examined archival Far Ultraviolet Spectroscopic Explorer data to search for far-ultraviolet emission lines in the starburst galaxy M82. The observations were made in an outflow region that extends beyond the galactic disk. We found the O vi λ λ 1032, 1038 emission lines from the galaxy’s southern outflow region. The O vi lines suggest that the outflowing warm-hot gas is undergoing radiative cooling. We measured a radial velocity of ∼420 km s ^−1 from the O vi lines, which is faster than the velocity seen in H α observations. The O vi λ 1038 emission line seems to be blended with the C ii λ 1037 emission line, which has a radial velocity of ∼300 km s ^−1 , similar to what is observed in H α observations. The outflow medium of M82 appears to be composed of gas in multiple phases with varying temperatures and kinematics. Future spectroscopic observations in high energy regimes covering a wider spatial area are necessary to understand better the properties of the warm-hot gas medium in the outflow.

We present a statistical study of the properties of diffuse H i in 10 nearby galaxies, comparing the H i detected by the single-dish telescope FAST (FEASTS program) and the interferometer Very Large Array (THINGS program), respectively. The THINGS observation missed H i with a median of 23% due to the short-spacing problem of interferometry and limited sensitivity. We extract the diffuse H i by subtracting the dense H i , which is obtained from the THINGS data with a uniform flux-density threshold, from the total H i detected by FAST. Among the sample, the median diffuse-H i fraction is 34%, and more diffuse H i is found in galaxies exhibiting more prominent tidal-interaction signatures. The diffuse H i we detected seems to be distributed in disk-like layers within a typical thickness of 1 kpc, different from the more halo-like diffuse H i detected around NGC 4631 in a previous study. Most of the diffuse H i is cospatial with the dense H i and has a typical column density of 10 ^17.7 –10 ^20.1 cm ^−2 . The diffuse and dense H i exhibit a similar rotational motion, but the former lags by a median of 25% in at least the inner disks, and its velocity dispersions are typically twice as high. Based on a simplified estimation of circumgalactic medium properties and assuming pressure equilibrium, the volume density of diffuse H i appears to be constant within each individual galaxy, implying its role as a cooling interface. Comparing with existing models, these results are consistent with a possible link between tidal interactions, the formation of diffuse H i , and gas accretion.

This paper presents a newly established sample of 103 unique galaxies or galaxy groups at 0.4 ≲ z ≲ 0.7 from the Cosmic Ultraviolet Baryon Survey (CUBS) for studying the warm-hot circumgalactic medium (CGM) probed by both O vi and Ne viii absorption. The galaxies and associated neighbors are identified at 40 km s ^−1 for galaxies of $\mathrm{log}{M}_{\mathrm{star}}/{M}_{\odot }\gt 9$ within the virial radius, suggesting a more dynamic warm-hot halo around these galaxies. Finally, the warm-hot CGM probed by O vi and Ne viii is suggested to be the dominant phase in sub- L * galaxies with $\mathrm{log}{M}_{\mathrm{star}}/{M}_{\odot }\approx 9\mbox{--}10$ based on their high ionization fractions in the CGM.

We present Selective Amplification of Ultra Noisy Astronomical Signal ( SAUNAS ), a pipeline designed for detecting diffuse X-ray emission in the data obtained with the Advanced CCD Imaging Spectrometer (ACIS) of the Chandra X-ray Observatory. SAUNAS queries the available observations in the Chandra archive and performs photometric calibration, point-spread function modeling and deconvolution, point-source removal, adaptive smoothing, and background correction. This pipeline builds on existing and well-tested software including CIAO , VorBin , and LIRA . We characterize the performance of SAUNAS through several quality performance tests and demonstrate the broad applications and capabilities of SAUNAS using two galaxies already known to show X-ray-emitting structures. SAUNAS successfully detects the 30 kpc X-ray superwind of NGC 3079 using Chandra/ACIS data sets, matching the spatial distribution detected with more sensitive XMM-Newton observations. The analysis performed by SAUNAS reveals an extended low surface brightness source in the field of UGC 5101 in the 0.3–1.0 keV and 1.0–2.0 keV bands. This source is potentially a background galaxy cluster or a hot gas plume associated with UGC 5101. SAUNAS demonstrates its ability to recover previously undetected structures in archival data, expanding exploration into the low surface brightness X-ray Universe with Chandra/ACIS.

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 smoothed particle hydrodynamics 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 log( M _* / M _⊙ ) = 9.5–11.5. We measure the fraction of metals remaining in the interstellar medium (ISM) and CGM of each galaxy and calculate the expected mass of each SMBH based on the M _BH – σ relation (Kormendy & Ho 2013). The deviation of each SMBH from its expected mass, Δ M _BH , is compared to the potential of its host via σ . We find that SMBHs with accreted mass above M _BH – σ are more effective at removing metals from the ISM than undermassive SMBHs in star-forming galaxies. Overall, overmassive SMBHs suppress the total star formation of their host galaxies and more effectively move metals from the ISM into the CGM. However, we see little to no evacuation of gas from the CGM out of their halos, in contrast with other simulations. Finally, we predict that C iv column densities in the CGM of L _* galaxies are unlikely to depend on host galaxy SMBH mass. Our results show that the scatter in the low-mass end of the M _BH – σ relation may indicate how effective an SMBH is in the local redistribution of mass in its host galaxy.

The extraplanar diffuse ionized gas (eDIG) represents ionized gases traced by optical/UV lines beyond the stellar extent of galaxies. We herein introduce a novel multislit narrow-band spectroscopy method to conduct spatially resolved spectroscopy of the eDIG around a sample of nearby edge-on disk galaxies (eDIG-CHANGES). In this paper, we introduce the project design and major scientific goals, as well as a pilot study of NGC 3556 (M108). The eDIG is detected to a vertical extent of a few kiloparsecs above the disk, comparable to the X-ray and radio images. We do not see significant vertical variation of the [N ii ]/H α line ratio. A rough examination of the pressure balance between different circumgalactic medium phases indicates the magnetic field is in a rough pressure balance with the X-ray emitting hot gas and may play an important role in the global motion of both the eDIG and the hot gas in the lower halo. At the location of an Hubble Space Telescope/Cosmic Origins Spectrograph observed UV bright background active galactic nucleus ∼29 kpc from the center of NGC 3556, the magnetic pressure is much lower than that of the hot gas and the ionized gas traced by UV absorption lines, although the extrapolation of the pressure profiles may cause some biases in this comparison. By comparing the position–velocity diagrams of the optical and CO lines, we also find the dynamics of the two gas phases are consistent with each other, with no evidence of a global inflow/outflow and a maximum rotation velocity of ∼150 km s ^−1 .

Most diffuse baryons, including the circumgalactic medium (CGM) surrounding galaxies and the intergalactic medium (IGM) in the cosmic web, remain unmeasured and unconstrained. Fast radio bursts (FRBs) offer an unparalleled method to measure the electron dispersion measures (DMs) of ionized baryons. Their distribution can resolve the missing baryon problem and constrain the history of feedback theorized to impart significant energy to the CGM and IGM. We analyze the Cosmology and Astrophysics with Machine Learning Simulations using three suites, IllustrisTNG, SIMBA, and Astrid, each varying six parameters (two cosmological and four astrophysical feedback), for a total of 183 distinct simulation models. We find significantly different predictions between the fiducial models of the suites owing to their different implementations of feedback. SIMBA exhibits the strongest feedback, leading to the smoothest distribution of baryons and reducing the sight-line-to-sight-line variance in DMs between z = 0 and 1. Astrid has the weakest feedback and the largest variance. We calculate FRB CGM measurements as a function of galaxy impact parameter, with SIMBA showing the weakest DMs due to aggressive active galactic nucleus (AGN) feedback and Astrid the strongest. Within each suite, the largest differences are due to varying AGN feedback. IllustrisTNG shows the most sensitivity to supernova feedback, but this is due to the change in the AGN feedback strengths, demonstrating that black holes, not stars, are most capable of redistributing baryons in the IGM and CGM. We compare our statistics directly to recent observations, paving the way for the use of FRBs to constrain the physics of galaxy formation and evolution.

We present an all-southern sky survey for Mg ii λλ 2796, 2803 doublet absorbers in 951 z < 4 active galactic nucleus/quasar spectra from the Australian Dark Energy Survey. The spectral resolution ranges from R = 1400 to 1700 over the wavelengths 3700–8800 Å. The survey has a 5 σ detection completeness of 50% and above for rest-frame equivalent widths W _r (2796) ≥ 0.3 Å. We studied 656 Mg ii absorption systems over the redshift range 0.33 ≤ z ≤ 2.19 with equivalent widths 0.3 ≤ W _r (2796) ≤ 3.45 Å. The equivalent width distribution is well fit by an exponential function with W _* = 0.76 ± 0.04 Å, and the redshift path density exhibits very little evolution. Overall, our findings are consistent with the large, predominantly northern sky, surveys of Mg ii absorbers. We developed and implemented a Monte Carlo model informed by a high-resolution Mg ii survey for determining the Mg ii mass density, Ω _Mg _II . We found Ω _Mg _II ∼ 5 × 10 ^−7 with no evidence of evolution over a ∼7 Gyr time span following Cosmic Noon. Incorporating measurements covering 2.0 ≤ z ≤ 6.4 from the literature, we extended our insights into Mg ii mass density evolution from the end of reionization well past the Cosmic Noon epoch. The presented Monte Carlo model has potential for advancing our knowledge of the evolution of mass densities of metal ions common to quasar absorption line studies, as it exploits the efficiency of large low-resolution surveys while requiring only small samples from expensive high-resolution surveys.

We report the discovery of large, ionized, [O ii ]-emitting circumgalactic nebulae around the majority of 30 UV-luminous quasars at z = 0.4–1.4 observed with deep, wide-field integral field spectroscopy with the Multi-Unit Spectroscopy Explorer (MUSE) by the Cosmic Ultraviolet Baryon Survey and MUSE Quasar Blind Emitters Survey. Among the 30 quasars, seven (23%) exhibit [O ii ]-emitting nebulae with major axis sizes greater than 100 kpc, 20 greater than 50 kpc (67%), and 27 (90%) greater than 20 kpc. Such large, optically emitting nebulae indicate that cool, dense, and metal-enriched circumgalactic gas is common in the halos of luminous quasars at intermediate redshift. Several of the largest nebulae exhibit morphologies that suggest interaction-related origins. We detect no correlation between the sizes and cosmological-dimming-corrected surface brightnesses of the nebulae and quasar redshift, luminosity, black hole mass, or radio-loudness, but find a tentative correlation between the nebulae and rest-frame [O ii ] equivalent width in the quasar spectra. This potential trend suggests a relationship between interstellar medium content and gas reservoirs on CGM scales. The [O ii ]-emitting nebulae around the z ≈ 1 quasars are smaller and less common than Ly α nebulae around z ≈ 3 quasars. These smaller sizes can be explained if the outer regions of the Ly α halos arise from scattering in more neutral gas, by evolution in the cool circumgalactic medium content of quasar-host halos, by lower-than-expected metallicities on ≳50 kpc scales around z ≈ 1 quasars, or by changes in quasar episodic lifetimes between z = 3 and 1.

Turbulent radiative mixing layers play an important role in many astrophysical contexts where cool (≲10 ^4 K) clouds interact with hot flows (e.g., galactic winds, high-velocity clouds, infalling satellites in halos and clusters). The fate of these clouds (as well as many of their observable properties) is dictated by the competition between turbulence and radiative cooling; however, turbulence in these multiphase flows remains poorly understood. We have investigated the emergent turbulence arising in the interaction between clouds and supersonic winds in hydrodynamic enzo-e simulations. In order to obtain robust results, we employed multiple metrics to characterize the turbulent velocity, v _turb . We find four primary results when cooling is sufficient for cloud survival. First, v _turb manifests clear temperature dependence. Initially, v _turb roughly matches the scaling of sound speed on temperature. In gas hotter than the temperature where cooling peaks, this dependence weakens with time until v _turb is constant. Second, the relative velocity between the cloud and wind initially drives rapid growth of v _turb . As it drops (from entrainment), v _turb starts to decay before it stabilizes at roughly half its maximum. At late times, cooling flows appear to support turbulence. Third, the magnitude of v _turb scales with the ratio between the hot phase sound-crossing time and the minimum cooling time. Finally, we find tentative evidence for a length scale associated with resolving turbulence. Underresolving this scale may cause violent shattering and affect the cloud’s large-scale morphological properties.

We present a new technique to identify associations of H i emission in the Magellanic Stream (MS) and ultraviolet (UV) absorbers from 92 QSO sight lines near the MS. We quantify the level of associations of individual H i elements to the main H i body of the Stream using Wasserstein distance-based models, and derive characteristic spatial and kinematic distances of the H i emission in the MS. With the emission-based model, we further develop a comparison metric, which identifies the dominant associations of individual UV absorbers with respect to the MS and nearby galaxies. For ionized gas associated with the MS probed by C ii , C iv , Si ii , Si iii , Si iv , we find that the ion column densities are generally ∼0.5 dex higher than those that are not associated, and that the gas is more ionized toward the tail of the MS as indicated by the spatial trend of the C ii /C iv ratios. For nearby galaxies, we identify potential new absorbers associated with the circumgalactic medium of M33 and NGC 300, and affirm the associations of absorbers with IC 1613 and the Wolf–Lundmark–Mellote galaxy. For M31, we find the previously identified gradient in column densities as a function of the impact parameter, and that absorbers with higher column densities beyond M31's virial radius are more likely to be associated with the MS. Our analysis of absorbers associated with the Magellanic Clouds reveals the presence of continuous and blended diffuse ionized gas between the Stream and the Clouds. Our technique can be applied to future applications of identifying associations within physically complex gaseous structures.

We have imaged the entirety of eight (plus one partial) Milky Way (MW)–like satellite systems, a total of 42 (45) satellites, from the Satellites Around Galactic Analogs II catalog in both H α and H i with the Canada–France–Hawaii Telescope and the Jansky Very Large Array. In these eight systems we have identified four cases where a satellite appears to be currently undergoing ram pressure stripping (RPS) as its H i gas collides with the circumgalactic medium (CGM) of its host. We also see a clear suppression of gas fraction ( M _HI / M _* ) with decreasing (projected) satellite–host separation—to our knowledge, the first time this has been observed in a sample of MW-like systems. Comparisons to the Auriga, A Project Of Simulating The Local Environment, and TNG50 cosmological zoom-in simulations show consistent global behavior, but they systematically underpredict gas fractions across all satellites by roughly 0.5 dex. Using a simplistic RPS model, we estimate the average peak CGM density that satellites in these systems have encountered to be $\mathrm{log}{\rho }_{\mathrm{cgm}}/{\rm{g}}\,{\mathrm{cm}}^{-3}\approx -27.3$ . Furthermore, we see tentative evidence that these satellites are following a specific star formation rate to gas fraction relation that is distinct from field galaxies. Finally, we detect one new gas-rich satellite in the UGC 903 system with an optical size and surface brightness meeting the standard criteria to be considered an ultra-diffuse galaxy.

The XMM-Newton Line Emission Analysis Program (X-LEAP) is designed to study diffuse X-ray emissions from the Milky Way (MW) hot gas, as well as emissions from the foreground solar wind charge exchange (SWCX). This paper presents an all-sky survey of spectral feature intensities corresponding to the O vii , O viii , and iron L-shell (Fe-L) emissions. These intensities are derived from 5418 selected XMM-Newton observations with long exposure times and minimal contamination from point or extended sources. For 90% of the measured intensities, the values are within ≈2–18 photons cm ^−2 s ^−1 sr ^−1 (line unit (LU)), ≈0–8 LU, and ≈0–9 LU, respectively. We report long-term variations in O vii and O viii intensities over 22 yr, closely correlating with the solar cycle and attributed to SWCX emissions. These variations contribute ∼30% and ∼20% to the observed intensities on average and peak at ≈4 and ≈1 LU during solar maxima. We also find evidence of short-term and spatial variations in SWCX, indicating the need for a more refined SWCX model in future studies. In addition, we present SWCX- and absorption-corrected all-sky maps for a better view of the MW hot gas emission. These maps show a gradual decrease in oxygen intensity moving away from the Galactic center and a concentration of Fe-L intensity in the Galactic bubbles and disk.

This paper presents the multiscale temperature structures in the Milky Way (MW) hot gas, as part of the XMM-Newton Line Emission Analysis Program, surveying the O vii , O viii , and Fe-L band emission features in the XMM-Newton archive. In particular, we define two temperature tracers, I _OVIII / I _OVII (O87) and I _FeL /( I _OVII + I _OVIII ) (FeO). These two ratios cannot be explained simultaneously using single-temperature collisional ionization models, which indicates the need for multitemperature structures in hot gas. In addition, we show three large-scale features in the hot gas: the eROSITA bubbles around the Galactic center (GC), the disk, and the halo. In the eROSITA bubbles, the observed line ratios can be explained by a log-normal temperature distribution with a median of $\mathrm{log}T/{\rm{K}}\approx 6.4$ and a scatter of σ _T ≈ 0.2 dex. Beyond the bubbles, the line ratio dependence on the Galactic latitude suggests higher temperatures around the midplane of the MW disk. The scale height of the temperature variation is estimated to be ≈2 kpc assuming an average distance of 5 kpc for the hot gas. The halo component is characterized by the dependence on the distance to the GC, showing a temperature decline from $\mathrm{log}\,T/{\rm{K}}\,\approx 6.3$ to 5.8. Furthermore, we extract the autocorrelation and cross-correlation functions to investigate the small-scale structures. O87 and FeO ratios show a consistent autocorrelation scale of ≈5° (i.e., ≈400 pc at 5 kpc), which is consistent with the expected physical sizes of X-ray bubbles associated with star-forming regions or supernova remnants. Finally, we examine the cross-correlation between the hot and UV-detected warm gas and show an intriguing anticorrelation.

We present the results from an analysis of multiwavelength archival data on the multiphase outflow in the starburst galaxy NGC 1808. We report the detection at 70 and 100 μ m of dust filaments that extend up to ∼13 kpc from the galactic midplane and trace an edge-brightened biconical structure along the minor axis of the galaxy. The inner filaments are roughly cospatial with previously identified optical dust filaments, extraplanar polycyclic aromatic hydrocarbon emission, and neutral and ionized gaseous outflows. The 70/160 μ m flux ratio, a proxy for the dust temperature, is elevated along the edges of the cones, indicating that the dusty medium has been driven out of the central regions of these cones and is possibly shock-heated by a large-scale outflow. We establish lower limits on the extraplanar dust mass and mean height above the stellar disk of $\mathrm{log}({M}_{d}/{M}_{\odot })\,=\,6.48$ and ∣ z ∣ ∼ 5 kpc. The energy requirement of (5.1–9.6) × 10 ^56 erg needed to lift the dusty material, assuming a Milky Way–like dust-to-gas ratio, can be supplied by the current starburst, with a measured star formation rate of 3.5–5.4 M _⊙ yr ^−1 over a timescale of (4–26) ξ ^−1 Myr, where ξ is the efficiency of energy transfer. We conclude that a starburst-driven outflow is the most likely mechanism by which the dust features were formed.

We present Atacama Large Millimeter/submillimeter Array (ALMA) deep spectroscopy for a lensed galaxy at z _spec = 8.496 with $\mathrm{log}({M}_{\mathrm{star}}/{M}_{\odot })\sim 7.8$ whose optical nebular lines and stellar continuum are detected by JWST/NIRSpec and NIRCam Early Release Observations in the field of SMACS J0723.3–7327. Our ALMA spectrum shows [O iii ] 88 μ m and [C ii ] 158 μ m line detections at 4.0 σ and 4.5 σ , respectively. The redshift and position of the [O iii ] line coincide with those of the JWST source, while the [C ii ] line is blueshifted by 90 km s ^−1 with a spatial offset of 0.″5 (≈0.5 kpc in the source plane) from the centroid of the JWST source. The NIRCam F444W image, including [O iii ] λ 5007 and H β line emission, spatially extends beyond the stellar components by a factor of >8. This indicates that the z = 8.5 galaxy has already experienced strong outflows as traced by extended [O iii ] λ 5007 and offset [C ii ] emission, which would promote ionizing photon escape and facilitate reionization. With careful slit-loss corrections and the removal of emission spatially outside the galaxy, we evaluate the [O iii ] 88 μ m/ λ 5007 line ratio, and derive the electron density n _e by photoionization modeling to be ${220}_{-130}^{+230}$ cm ^−3 , which is comparable with those of z ∼ 2–3 galaxies. We estimate an [O iii ] 88 μ m/[C ii ] 158 μ m line ratio in the galaxy of >4, as high as those of known z ∼ 6–9 galaxies. This high [O iii ] 88 μ m/[C ii ] 158 μ m line ratio is generally explained by the high n _e as well as the low metallicity ( ${Z}_{\mathrm{gas}}/{Z}_{\odot }={0.04}_{-0.02}^{+0.02}$ ), high ionization parameter ( $\mathrm{log}U\gt -2.27$ ), and low carbon-to-oxygen abundance ratio (log(C/O) = [−0.52: −0.24]) obtained from the JWST/NIRSpec data; further [C ii ] follow-up observations will constrain the covering fraction of photodissociation regions.

Testing the standard cosmological model (ΛCDM) at small scales is challenging. Galaxies that inhabit low-mass dark matter halos provide an ideal test bed for dark matter models by linking observational properties of galaxies at small scales (low mass, low velocity) to low-mass dark matter halos. However, the observed kinematics of these galaxies do not align with the kinematics of the dark matter halos predicted to host them, obscuring our understanding of the low-mass end of the galaxy–halo connection. We use deep H i observations of low-mass galaxies at high spectral resolution in combination with cosmological simulations of dwarf galaxies to better understand the connection between dwarf galaxy kinematics and low-mass halos. Specifically, we use H i line widths to directly compare to the maximum velocities in a dark matter halo and find that each deeper measurement approaches the expected one-to-one relationship between the observed kinematics and the predicted kinematics in ΛCDM. We also measure baryonic masses and place these on the baryonic Tully–Fisher relation (BTFR). Again, our deepest measurements approach the theoretical predictions for the low-mass end of this relation, a significant improvement on similar measurements based on line widths measured at 50% and 20% of the peak. Our data also hint at the rollover in the BTFR predicted by hydrodynamical simulations of ΛCDM for low-mass galaxies.

We report new observations of the cool diffuse gas around 29, 2.3 < z < 6.3 galaxies using deep JWST/NIRCam slitless grism spectroscopy around the sight line to the quasar J0100+2802. The galaxies span a stellar mass range of $7.1\leqslant {\rm{log}}{M}_{\ast }/{M}_{\odot }\leqslant 10.7$ , and star formation rates (SFRs) of $-0.1\lt \mathrm{log}$ SFR/ M _⊙ yr ^−1 < 2.3. We find galaxies for seven Mg ii absorption systems within 300 kpc of the quasar sight line. The Mg ii radial absorption profile falls off sharply with radius, with most of the absorption extending out to 2–3 R _200 of the host galaxies. Six out of seven Mg ii absorption systems are detected around galaxies with ${\rm{log}}{M}_{\ast }/{M}_{\odot }\gt $ 9. The Mg ii absorption kinematics are shifted from the systemic redshifts of host galaxies with a median absolute velocity ≈ 135 km s ^−1 and standard deviation ≈ 85 km s ^−1 . The high kinematic offset and large radial separation ( R > 1.3 R _200 ), suggest that five out of the seven Mg ii absorption systems are not gravitationally bound to their host galaxy. In contrast, most of the cool circumgalactic medium at z < 1 is gravitationally bound. The high incidence of unbound Mg ii gas in this work suggests that toward the end of reionization, galaxy halos are in a state of remarkable disequilibrium, and are highly efficient in enriching the intergalactic medium. The two strongest Mg ii absorption systems are detected at z ∼ 4.22 and 4.5, the former associated with a merging galaxy system and the latter associated with three kinematically close galaxies. Both of these galaxies reside in local galaxy overdensities, indicating the presence of cool Mg ii absorption in two “protogroups” at z > 4.

Turbulent motions in the circumgalactic medium play a critical role in regulating the evolution of galaxies, yet their detailed characterization remains elusive. Using two-dimensional velocity maps constructed from spatially extended [O ii ] and [O iii ] emission, Chen et al. measured the velocity structure functions (VSFs) of four quasar nebulae at z ≈ 0.5–1.1. One of these exhibits a spectacular Kolmogorov relation. Here, we carry out an ensemble study using an expanded sample incorporating four new nebulae from three additional quasi-stellar object (QSO) fields. The VSFs measured for all eight nebulae are best explained by subsonic turbulence revealed by the line-emitting gas, which in turn strongly suggests that the cool gas ( T ∼ 10 ^4 K) is dynamically coupled to the hot ambient medium. Previous work demonstrates that the largest nebulae in our sample reside in group environments with clear signs of tidal interactions, suggesting that environmental effects are vital in seeding and enhancing the turbulence within the gaseous halos, ultimately promoting the formation of the extended nebulae. No discernible differences are observed in the VSF properties between radio-loud and radio-quiet QSO fields. We estimate the turbulent heating rate per unit volume, Q _turb , in the QSO nebulae to be ∼10 ^−26 –10 ^−22 erg cm ^−3 s ^−1 for the cool phase and ∼10 ^−28 –10 ^−25 erg cm ^−3 s ^−1 for the hot phase. This range aligns with measurements in the intracluster medium and star-forming molecular clouds but is ∼10 ^3 times higher than the Q _turb observed inside cool gas clumps on scales ≲1 kpc using absorption-line techniques. We discuss the prospect of bridging the gap between emission and absorption studies by pushing the emission-based VSF measurements to below ≈10 kpc.

We analyze the circumgalactic medium (CGM) for eight commonly-used cosmological codes in the AGORA collaboration. The codes are calibrated to use identical initial conditions, cosmology, heating and cooling, and star formation thresholds, but each evolves with its own unique code architecture and stellar feedback implementation. Here, we analyze the results of these simulations in terms of the structure, composition, and phase dynamics of the CGM. We show properties such as metal distribution, ionization levels, and kinematics are effective tracers of the effects of the different code feedback and implementation methods, and as such they can be highly divergent between simulations. This is merely a fiducial set of models, against which we will in the future compare multiple feedback recipes for each code. Nevertheless, we find that the large parameter space these simulations establish can help disentangle the different variables that affect observable quantities in the CGM, e.g., showing that abundances for ions with higher ionization energy are more strongly determined by the simulation’s metallicity, while abundances for ions with lower ionization energy are more strongly determined by the gas density and temperature.