Damped Lyα Research Articles

Improved Astrophysical and Computational Oscillator Strengths for Ultraviolet P ii Lines

Abundance measurements for the volatile element phosphorus are important for measuring the metallicity in interstellar and circumgalactic gas, where their accuracies are limited by uncertainties in the oscillator strengths. We report updated oscillator strength values for two resonant transitions of the dominant ion P II, the transitions at 961.041 and 963.801 Å, which have historically shown large uncertainties. Using a combination of observational measurements and highly accurate quasi-relativistic Hartree–Fock theoretical calculations, we present an updated oscillator strength of f = 0.147 ± 0.021 for the poorly constrained P II resonant transition at 961.401 Å, which arises from the ground electronic state 3s 23p 2 3P0 to the excited level 3s 23p3d 3D 1o . This result utilizes archival optical spectra obtained with the Very Large Telescope for the quasar PKS 0528–250, which has a damped Lyα absorber at z = 2.811. We calculate a theoretical f-value = 0.153 for 961.401 Å consistent with our empirically derived value, and calculate a theoretical f-value = 1.79 for 963.801 Å. We also present theoretical oscillator strengths for the P II resonant transitions at 972.779, 1124.945, 1152.818, 1301.874, and 1532.533 Å, as well as for multiple P II fine-structure and excited-level transitions. The updated f-value for the P II 961 Å transition will be useful in future studies of P abundances, especially in sight lines where the 1152 Å line is saturated.

The rich galactic environment of a H2-absorption-selected quasar

We present the first Very Large Telescope (VLT) Multi Unit Spectroscopic Explorer (MUSE) observations of a quasar featuring a proximate molecular absorption system, SDSS J125917.31+030922.5. The proximate damped Lyα absorption acts as a natural coronagraph, removing the quasar emission over ∼40 Å in wavelength, and allows us to detect extended Lyα emission without the necessity of subtracting the quasar emission. This natural coronagraph permits the investigation of the quasar environment down to its inner regions (r < 20 kpc), where galaxy interactions or feedback processes should have the most noticeable effects. Our observations reveal a dense environment, with a highly asymmetric Lyα emission within 2″ (∼15 kpc), possibly shaped by a companion galaxy, and a southern extension of the nebulae to about 50 kpc, with rotation-like kinematic signature. The width of the Lyα emission is broadest closer to the quasar, indicating perturbed kinematics as expected if interactions and significant gas flows are present. The foreground absorbing system itself is redshifted by ≈400 km/s relative to the background quasar, and therefore is likely arising from gas moving toward the quasar. Finally two additional Lyα emitters are detected with > 10σ significance at 96 and 223 kpc from the quasar, making this field overdense relative to other similar observations of quasars at z ∼ 3. Our results support the hypothesis that quasars with proximate neutral and molecular absorption trace rich environments where galaxy interactions are at play and motivates further studies of H2-selected quasars to shed light on feeding and feedback processes.

Detecting Novel Damped Lyα Absorbers in the Intergalactic Medium

Damped Lyα (DLA) absorbers are a class of quasar absorption line systems that can be used to trace hydrogen evolution and study gas, composition, kinematics, and stellar formation in galaxies or the intergalactic medium in the early universe. In this study, the number of DLA absorbers found with a deep learning algorithm is augmented using metal Mg ii absorbers as signposts. The artificial training data was generated by inserting synthetic absorption lines into spectra. Using previously identified DLAs as a test set, the model returned a 94.6% accuracy score and 1197 new DLA absorbers, which nearly doubles the previously known number of samples within the search catalog. Additionally, the column density and Doppler broadening are measured on the newly found DLA absorbers through Voigt profile fitting.

Circumnuclear Gas of a Quasar in Absorption and Emission Lines

We present a detailed analysis of the unusual absorption- and emission-line spectrum of the quasar SDSS J074850.39+442439.0 (hereafter J0748+4424) at a redshift of z = 2.18. The archival SDSS optical spectrum is abundant in narrow absorption lines (NALs) originating from mostly singly ionized iron, nickel, and silicon ions at both ground and excited levels. With the aid of the photoionization simulations, we find that these NALs can be reasonably well reproduced by a homogeneous gas slab with a neutral hydrogen column density of 1021.6–1021.9 cm−2 and a number density of 106.0–106.6 cm−3 illuminated by the quasar central engine with an ionization parameter of 10−2.8–10−2.1. We infer the gas is located at a distance of ∼30–130 pc from the black hole. This circumnuclear gas can reveal itself in emission lines of a width of ∼1000 km s−1. Such intermediate-width emission lines (IELs) are indeed dominant in the observed Lyα emission, as the Lyα broad emission line (BEL) is mostly eaten up by the expected damped Lyα absorption (DLA). The IELs, though rather weak, are also detected in Hα and [O iii] emission lines. The IEL intensities largely agree with the predicted values by the simplified photoionization models using the parameters derived from the absorption lines and a covering fraction of 10%. The intrinsic DLA in J0748+4424, while itself contains abundant information on the circumnuclear environment, may serve as a “natural coronagraph” at the line of sight, exposing the IELs that are otherwise overwhelmed by the BELs.

Impact of Self-shielding Minihalos on the Lyα Forest at High Redshift

Dense gas in minihalos with masses of 106−108 M ⊙ can shield themselves from reionization for ∼100 Myr after being exposed to the UV background. These self-shielded systems, often unresolved in cosmological simulations, can introduce strong absorption in quasar spectra. This paper is the first systematic study on the impact of these systems on the Lyα forest. We first derive the H i column density profile of photoevaporating minihalos by conducting 1D radiation–hydrodynamics simulations. We utilize these results to estimate the Lyα opacity from minihalos in a large-scale simulation that cannot resolve self-shielding. When the ionization rate of the background radiation is 0.03 × 10−12 s−1, as expected near the end of reionization at z ∼ 5.5, we find that the incidence rate of damped Lyα absorbers increases by a factor of ∼2−4 compared to at z = 4.5. The Lyα flux is, on average, suppressed by ∼3% of its mean due to minihalos. The absorption features enhance the 1D power spectrum up to ∼5% at k ∼ 0.1 h Mpc−1 (or 10−3 km−1 s), which is comparable to the enhancement caused by inhomogeneous reionization. The flux is particularly suppressed in the vicinity of large halos along the line-of-sight direction at separations of up to 10 h −1 Mpc at r ⊥ ≲ 2 h −1 Mpc. However, these effects become much smaller for higher ionizing rates (≳0.3 × 10−12 s−1) expected in the post-reionization Universe. Our findings highlight the need to consider minihalo absorption when interpreting the Lyα forest at z ≳ 5.5. Moreover, the sensitivity of these quantities to the ionizing background intensity can be exploited to constrain the intensity itself.

No Top-heavy Stellar Initial Mass Function Needed: The Ionizing Radiation of GS9422 Can Be Powered by a Mixture of an Active Galactic Nucleus and Stars

JWST is producing high-quality rest-frame optical and UV spectra of faint galaxies at z > 4 for the first time, challenging models of galaxy and stellar populations. One galaxy recently observed at z = 5.943, GS9422, has nebular line and UV continuum emission that appears to require a high ionizing photon production efficiency. This has been explained with an exotic stellar initial mass function (IMF; Cameron et al. 2023a), 10–30x more top-heavy than a Salpeter IMF. Here we suggest an alternate explanation to this exotic IMF. We use a new flexible neural net emulator for CLOUDY, Cue, to infer the shape of the ionizing spectrum directly from the observed emission line fluxes. By describing the ionizing spectrum with a piecewise power law, Cue is agnostic to the source of the ionizing photons. Cue finds that the ionizing radiation from GS9422 can be approximated by a double power law characterized by QHeIIQH=−1.5 , which can be interpreted as a combination of young, metal-poor stars and a low-luminosity active galactic nucleus with F ν ∝ λ 2 in a 65%/35% ratio. This suggests a significantly lower nebular continuum contribution to the observed UV flux (24%) than a top-heavy IMF (≳80%), and hence, necessitates a damped Lyα absorber to explain the continuum turnover bluewards of ∼1400 Å. While current data cannot rule out either scenario, given the immense impact the proposed top-heavy IMF would have on models of galaxy formation, it is important to propose viable alternative explanations and to further investigate the nature of peculiar high-z nebular emitters.

Nebular dominated galaxies: insights into the stellar initial mass function at high redshift

ABSTRACT We identify a low-metallicity ($12+\log ({\rm O}/{\rm H})=7.59$) Ly $\alpha$-emitting galaxy at $z=5.943$ with evidence of a strong Balmer jump, arising from nebular continuum. While Balmer jumps are sometimes observed in low-redshift star-forming galaxies, this galaxy also exhibits a steep turnover in the UV continuum. Such turnovers are typically attributed to absorption by a damped Ly $\alpha$ system (DLA); however, the shape of the turnover and the high observed Ly $\alpha$ escape fraction ($f_{\rm esc,Ly\alpha }~\sim 27~{{\ \rm per\ cent}}$) is also consistent with strong nebular two-photon continuum emission. Modelling the UV turnover with a DLA requires extreme column densities ($N_{\rm HI}\,\,\gt\,\, 10^{23}$ cm$^{-2}$), and simultaneously explaining the high $f_{\rm esc,Ly\alpha }$ requires a fine-tuned geometry. In contrast, modelling the spectrum as primarily nebular provides a good fit to both the continuum and emission lines, motivating scenarios in which (a) we are observing only nebular emission or (b) the ionizing source is powering extreme nebular emission that outshines the stellar emission. The nebular-only scenario could arise if the ionizing source has ‘turned off’ more recently than the recombination time-scale ($\sim$1000 yr), hence we may be catching the object at a very specific time. Alternatively, hot stars with $T_{\rm eff}\gtrsim 10^5$ K (e.g. Wolf–Rayet or low-metallicity massive stars) produce enough ionizing photons such that the two-photon emission becomes visible. While several stellar SEDs from the literature fit the observed spectrum well, the hot-star scenario requires that the number of $\gtrsim 50~{\rm M}_\odot$ stars relative to $\sim 5\!-\!50~{\rm M}_\odot$ stars is significantly higher than predicted by typical stellar initial mass functions (IMFs). The identification of more galaxies with similar spectra may provide evidence for a top-heavy IMF at high redshift.

METAL-Z: Measuring Dust Depletion in Low-metallicity Dwarf Galaxies

The cycling of metals between interstellar gas and dust is a critical aspect of the baryon cycle of galaxies, yet our understanding of this process is limited. This study focuses on understanding dust depletion effects in the low-metallicity regime (<20% Z ⊙) typical of cosmic noon. Using medium-resolution UV spectroscopy from the Cosmic Origins Spectrograph on board the Hubble Space Telescope, gas-phase abundances and depletions of iron and sulfur were derived toward 18 sight lines in local dwarf galaxies IC 1613 and Sextans A. The results show that the depletion of Fe and S is consistent with that found in the Milky Way (MW), LMC, and SMC. The depletion level of Fe increases with gas column density, indicating dust growth in the interstellar medium. The level of Fe depletion decreases with decreasing metallicity, resulting in the fraction of iron in gas ranging from 3% in the MW to 9% in IC 1613 and ∼19% in Sextans A. The dust-to-gas and dust-to-metal ratios (D/G, D/M) for these dwarf galaxies were estimated based on the MW relations between the depletion of Fe and other elements. The study finds that D/G decreases only slightly sublinearly with metallicity, with D/M decreasing from 0.41 ± 0.05 in the MW to 0.11 ± 0.11 at 0.10 Z ⊙ (at log N(H) = 21 cm−2). The trend of D/G versus metallicity using depletion in local systems is similar to that inferred in Damped Lyα systems from abundance ratios but lies higher than the trend inferred from far-IR measurements in nearby galaxies.

Polarimetry of the Lyα envelope of the radio-quiet quasar SDSS J124020.91+145535.6

The radio quiet quasar SDSS J1240+1455 lies at a redshift of z = 3.11, is surrounded by a Lyα blob (LAB), and is absorbed by a proximate damped Lyα system. In order to better define the morphology of the blob and determine its emission mechanism, we gathered deep narrow-band images isolating the Lyα line of this object in linearly polarized light. We provide a deep intensity image of the blob, showing a filamentary structure extending up to 16″ (or 122 physical kpc) in diameter. No significant polarization signal could be extracted from the data, but 95% probability upper limits were defined through simulations. They vary between ∼3% in the central 0.75″ disk (after subtraction of the unpolarized quasar continuum) and ∼10% in the 3.8 − 5.5″ annulus. The low polarization suggests that the Lyα photons are emitted mostly in situ, by recombination and de-excitation in a gas largely ionized by the quasar ultraviolet light, rather than by a central source and scattered subsequently by neutral hydrogen gas. This blob shows no detectable polarization signal, contrary to LAB1, a brighter and more extended blob that is not related to the nearby active galactic nucleus (AGN) in any obvious way, and where a significant polarization signal of about 18% was detected.

A Survey of Lyα Emission around Damped Lyα Absorbers at z ≈ 2 with the Keck Cosmic Web Imager

We present Keck Cosmic Web Imager Lyα integral field spectroscopy of the fields surrounding 14 damped Lyα absorbers (DLAs) at z ≈ 2. Of these 14 DLAs, nine have high metallicities ([M/H] > − 0.3), and four of those nine feature a CO-emitting galaxy at an impact parameter ≲30 kpc. Our search reaches median Lyα line flux sensitivities of ∼2 × 10−17 erg s−1 cm−2 over apertures of ∼6 kpc and out to impact parameters of ∼50 kpc. We recover the Lyα flux of three known Lyα-emitting H i-selected galaxies in our sample. In addition, we find two Lyα emitters at impact parameters of ≈50–70 kpc from the high-metallicity DLA at z ≈ 1.96 toward QSO B0551-366. This field also contains a massive CO-emitting galaxy at an impact parameter of ≈15 kpc. Apart from the field with QSO B0551-366, we do not detect significant Lyα emission in any of the remaining eight high-metallicity DLA fields. Considering the depth of our observations and our ability to recover previously known Lyα emitters, we conclude that H i-selected galaxies associated with high-metallicity DLAs at z ≈ 2 are dusty and therefore might feature low Lyα escape fractions. Our results indicate that complementary approaches—using Lyα, CO, Hα, and [C ii] 158 μm emission—are necessary to identify the wide range of galaxy types associated with z ≈ 2 DLAs.

HST FUV Spectroscopy of Super Star Cluster A in the Green Pea Analog Mrk 71: Revealing the Presence of Very Massive Stars

Mrk 71 is a low-metallicity (Z = 0.16 Z ☉) starburst region in the local dwarf galaxy NGC 2366, hosting two super star clusters (SSCs A and B), and it is recognized as a Green Pea (GP) analog with SSC A responsible for the GP properties. We present STIS and FOS far-ultraviolet (FUV) spectra of the embedded SSC Mrk 71-A obtained with the Hubble Space Telescope. The STIS FUV spectrum shows the characteristic features of very massive stars (VMS; masses >100 M ⊙) and we derive an age of 1 ± 1 Myr by comparison with the Charlot & Bruzual suite of spectral population synthesis models with upper mass limits of 300 and 600 M ⊙. We compare the STIS spectrum with all known SSC spectra exhibiting VMS signatures: NGC 5253-5, R136a, NGC 3125-A1, and the z = 2.37 Sunburst cluster. We find that the cluster mass-loss rates and wind velocities, as characterized by the C iv P Cygni profiles and the He ii emission line strengths, are very similar over Z = 0.16–0.4 Z☉. This agrees with predictions that the optically thick winds of VMS will be enhanced near the Eddington limit and show little metallicity dependence. We find very strong damped Lyα absorption with N(H i) =1022.2 cm−2 associated with Mrk 71-A. We discuss the natal environment of this young SSC in terms of radiatively driven winds, catastrophic cooling, and recent models where the cluster is surrounded by highly pressurized clouds with large neutral columns.

DESI z ≳ 5 Quasar Survey. I. A First Sample of 400 New Quasars at z ∼ 4.7–6.6

We report the first results of a high-redshift (z ≳ 5) quasar survey using the Dark Energy Spectroscopic Instrument (DESI). As a DESI secondary target program, this survey is designed to carry out a systematic search and investigation of quasars at 4.8 < z < 6.8. The target selection is based on the DESI Legacy Imaging Surveys (the Legacy Surveys) DR9 photometry, combined with the Pan-STARRS1 data and J-band photometry from public surveys. A first quasar sample has been constructed from the DESI Survey Validation 3 (SV3) and first-year observations until 2022 May. This sample includes more than 400 new quasars at redshift 4.7 ≤ z < 6.6, down to 21.5 magnitude (AB) in the z band, discovered from 35% of the entire target sample. Remarkably, there are 220 new quasars identified at z ≥ 5, more than one-third of existing quasars previously published at this redshift. The observations so far result in an average success rate of 23% at z > 4.7. The current spectral data set has already allowed analysis of interesting individual objects (e.g., quasars with damped Lyα absorbers and broad absorption line features), and statistical analysis will follow the survey’s completion. A set of science projects will be carried out leveraging this program, including quasar luminosity function, quasar clustering, intergalactic medium, quasar spectral properties, intervening absorbers, and properties of early supermassive black holes. Additionally, a sample of 38 new quasars at z ∼ 3.8–5.7 discovered from a pilot survey in the DESI SV1 is also published in this paper.

Quasar Factor Analysis—An Unsupervised and Probabilistic Quasar Continuum Prediction Algorithm with Latent Factor Analysis

Since their first discovery, quasars have been essential probes of the distant Universe. However, due to our limited knowledge of its nature, predicting the intrinsic quasar continua has bottlenecked their usage. Existing methods of quasar continuum recovery often rely on a limited number of high-quality quasar spectra, which might not capture the full diversity of the quasar population. In this study, we propose an unsupervised probabilistic model, quasar factor analysis (QFA), which combines factor analysis with physical priors of the intergalactic medium to overcome these limitations. QFA captures the posterior distribution of quasar continua through generatively modeling quasar spectra. We demonstrate that QFA can achieve the state-of-the-art performance, ∼2% relative error, for continuum prediction in the Lyα forest region compared to previous methods. We further fit 90,678 2 < z < 3.5, signal-to-noise ratio >2 quasar spectra from Sloan Digital Sky Survey Data Release 16 and found that for ∼30% quasar spectra where the continua were ill-determined with previous methods, QFA yields visually more plausible continua. QFA also attains ≲1% error in the 1D Lyα power spectrum measurements at z ∼ 3 and ∼4% in z ∼ 2.4. In addition, QFA determines latent factors representing more physical motivation than principal component analysis. We investigate the evolution of the latent factors and report no significant redshift or luminosity dependency except for the Baldwin effect. The generative nature of QFA also enables outlier detection robustly; we showed that QFA is effective in selecting outlying quasar spectra, including damped Lyα systems and potential Type II quasar spectra.

CLASSY VII Lyα Profiles: The Structure and Kinematics of Neutral Gas and Implications for LyC Escape in Reionization-era Analogs* * Based on observations made with the NASA/ESA Hubble Space Telescope, obtained from the Data Archive at the Space Telescope Science Institute, which is operated by the Association of Universities for Research in Astronomy, Inc., under NASA contract NAS 5-26555.

Lyα line profiles are a powerful probe of interstellar medium (ISM) structure, outflow speed, and Lyman-continuum escape fraction. In this paper, we present the Lyα line profiles of the Cosmic Origins Spectrograph (COS) Legacy Archive Spectroscopic SurveY, a sample rich in spectroscopic analogs of reionization-era galaxies. A large fraction of the spectra show a complex profile, consisting of a double-peaked Lyα emission profile in the bottom of a damped, Lyα absorption trough. Such profiles reveal an inhomogeneous ISM. We successfully fit the damped Lyα absorption and the Lyα emission profiles separately, but with complementary covering factors, a surprising result because this approach requires no Lyα exchange between high-N H i and low-N H i paths. The combined distribution of column densities is qualitatively similar to the bimodal distributions observed in numerical simulations. We find an inverse relation between Lyα peak separation and the [O iii]/[O ii] flux ratio, confirming that the covering fraction of Lyman-continuum-thin sightlines increases as the Lyα peak separation decreases. We combine measurements of Lyα peak separation and Lyα red peak asymmetry in a diagnostic diagram, which identifies six Lyman-continuum leakers in the COS Legacy Archive Spectrocopy SurveY (CLASSY) sample. We find a strong correlation between the Lyα trough velocity and the outflow velocity measured from interstellar absorption lines. We argue that greater vignetting of the blueshifted Lyα peak, relative to the redshifted peak, is the source of the well-known discrepancy between shell-model parameters and directly measured outflow properties. The CLASSY sample illustrates how scattering of Lyα photons outside the spectroscopic aperture reshapes Lyα profiles because the distances to these compact starbursts span a large range.

On the Metallicities and Kinematics of the Circumgalactic Media of Damped Lyα Systems at z ∼ 2.5* * The data presented herein were obtained at the W. M. Keck Observatory, which is operated as a scientific partnership among the California Institute of Technology, the University of California and the National Aeronautics and Space Administration. The Observatory was made possible by

We use medium- and high-resolution spectroscopy of close pairs of quasars to analyze the circumgalactic medium (CGM) surrounding 32 damped Lyα absorption systems (DLAs). The primary quasar sightline in each pair probes an intervening DLA in the redshift range 1.6 < z abs < 3.5, such that the secondary sightline probes absorption from Lyα and a large suite of metal-line transitions (including O i, C ii, C iv, Si ii, and Si iv) in the DLA host galaxy’s CGM at transverse distances 24 kpc ≤ R ⊥ ≤ 284 kpc. Analysis of Lyα in the CGM sightlines shows an anticorrelation between R ⊥ and H i column density (N HI) with 99.8% confidence, similar to that observed around luminous galaxies. The incidences of C ii and Si ii with N > 1013 cm−2 within 100 kpc of DLAs are larger by 2σ than those measured in the CGM of Lyman break galaxies (C f (N CII ) > 0.89 and ). Metallicity constraints derived from ionic ratios for nine CGM systems with negligible ionization corrections and N HI > 1018.5 cm−2 show a significant degree of scatter (with metallicities/limits across the range ), suggesting inhomogeneity in the metal distribution in these environments. Velocity widths of C iv λ1548 and low-ionization metal species in the DLA versus CGM sightlines are strongly (>2σ) correlated, suggesting that they trace the potential well of the host halo over R ⊥ ≲ 300 kpc scales. At the same time, velocity centroids for C iv λ1548 differ in DLA versus CGM sightlines by >100 km s−1 for ∼50% of velocity components, but few components have velocities that would exceed the escape velocity assuming dark matter host halos of ≥1012 M ⊙.

Metal-enriched Neutral Gas Reservoir around a Strongly Lensed Low-mass Galaxy at z = 4 Identified by JWST/NIRISS and VLT/MUSE

Direct observations of low-mass, low-metallicity galaxies at z ≳ 4 provide an indispensable opportunity for detailed inspection of the ionization radiation, gas flow, and metal enrichment in sources similar to those that reionized the universe. Combining the James Webb Space Telescope (JWST), Very Large Telescope/MUSE, and Atacama Large Millimeter/submillimeter Array, we present detailed observations of a strongly lensed, low-mass (≈107.6 M ⊙) galaxy at z = 3.98 (also see Vanzella et al.). We identify strong narrow nebular emission, including C iv λ λ1548, 1550, He ii λ1640, O iii] λ λ1661, 1666, [Ne iii] λ3868, [O ii] λ3727, and the Balmer series of hydrogen from this galaxy, indicating a metal-poor H ii region (≲0.12 Z ⊙) powered by massive stars. Further, we detect a metal-enriched damped Lyα system (DLA) associated with the galaxy with the H i column density of N H I ≈ 1021.8 cm−2. The metallicity of the associated DLA may reach the supersolar metallicity (≳Z ⊙). Moreover, thanks to JWST and gravitational lensing, we present the resolved UV slope (β) map at the spatial resolution of ≈100 pc at z = 4, with steep UV slopes reaching β ≈ −2.5 around three star-forming clumps. Combining with low-redshift analogs, our observations suggest that low-mass, low-metallicity galaxies, which dominate reionization, could be surrounded by a high covering fraction of the metal-enriched, neutral-gaseous clouds. This implies that the metal enrichment of low-mass galaxies is highly efficient, and further supports that in low-mass galaxies, only a small fraction of ionizing radiation can escape through the interstellar or circumgalactic channels with low-column-density neutral gas.

KODIAQ-Z: Metals and Baryons in the Cool Intergalactic and Circumgalactic Gas at 2.2 ≲ z ≲ 3.6

We present the KODIAQ-Z survey aimed to characterize the cool, photoionized gas at 2.2 ≲ z ≲ 3.6 in 202 H i-selected absorbers with 14.6 ≤ < 20 that probe the interface between galaxies and the intergalactic medium (IGM). We find that gas with at 2.2 ≲ z ≲ 3.6 can be metal-rich (−1.6 ≲ [X/H] ≲ − 0.2) as seen in damped Lyα absorbers (DLAs); it can also be very metal-poor ([X/H] < − 2.4) or even pristine ([X/H] < − 3.8), which is not observed in DLAs but is common in the IGM. For absorbers, the frequency of pristine absorbers is about 1%–10%, while for absorbers it is 10%–20%, similar to the diffuse IGM. Supersolar gas is extremely rare (<1%) at these redshifts. The factor of several thousand spread from the lowest to highest metallicities and large metallicity variations (a factor of a few to >100) between absorbers separated by less than Δv < 500 km s−1 imply that the metals are poorly mixed in gas. We show that these photoionized absorbers contribute to about 14% of the cosmic baryons and 45% of the cosmic metals at 2.2 ≲ z ≲ 3.6. We find that the mean metallicity increases with N H i , consistent with what is found in z < 1 gas. The metallicity of gas in this column density regime has increased by a factor ∼8 from 2.2 ≲ z ≲ 3.6 to z < 1, but the contribution of the absorbers to the total metal budget of the universe at z < 1 is a quarter of that at 2.2 ≲ z ≲ 3.6. We show that FOGGIE cosmological zoom-in simulations have a similar evolution of [X/H] with N H i , which is not observed in lower-resolution simulations. In these simulations, very metal-poor absorbers with [X/H] < − 2.4 at z ∼ 2–3 are tracers of inflows, while higher-metallicity absorbers are a mixture of inflows and outflows.

METAL: The Metal Evolution, Transport, and Abundance in the Large Magellanic Cloud Hubble Program. IV. Calibration of Dust Depletions versus Abundance Ratios in the Milky Way and Magellanic Clouds and Application to Damped Lyα Systems

The evolution of the metal content of the universe can be tracked through rest-frame UV spectroscopy of damped Lyα systems (DLAs). Gas-phase abundances in DLAs must be corrected for dust depletion effects, which can be accomplished by calibrating the relation between abundance ratios such as [Zn/Fe] and depletions (the fraction of metals in gas, as opposed to dust). Using samples of gas-phase abundances and depletions in the Milky Way (MW), LMC, and SMC, we demonstrate that the relation between [Zn/Fe] and other abundance ratios does not change significantly between these local galaxies and DLAs, indicating that [Zn/Fe] should trace depletions of heavy elements in those systems. The availability of photospheric abundances in young massive stars, a proxy for the total (gas+dust) metallicity of neutral gas, in the MW, LMC, and SMC allows us to calibrate the relation between [Zn/Fe] and depletions in these nearby galaxies. We apply the local calibrations of depletions to DLAs. We find that the fraction of metals in dust, the dust-to-gas ratio, and total abundances are 2–5 times lower than inferred from previous depletion calibrations based on MW measurements and a different formalism. However, the trend of dust abundance versus metallicity remains only slightly sublinear for all existing depletion calibrations, contrary to what is inferred from far-IR (FIR), 21 cm, and CO emission in nearby galaxies and predicted by chemical evolution models. Observational constraints on the FIR dust opacity and depletions at metallicities lower than 20% solar will be needed to resolve this tension.

The ALMA REBELS Survey: The Cosmic H i Gas Mass Density in Galaxies at z ≈ 7

The neutral atomic gas content of individual galaxies at large cosmological distances has until recently been difficult to measure due to the weakness of the hyperfine H i 21 cm transition. Here we estimate the H i gas mass of a sample of main-sequence star-forming galaxies at z ∼ 6.5–7.8 surveyed for [C ii] 158 μm emission as part of the Reionization Era Bright Emission Line Survey (REBELS), using a recent calibration of the [C ii]-to-H i conversion factor. We find that the H i gas mass excess in galaxies increases as a function of redshift, with an average of M Hi /M ⋆ ≈ 10, corresponding to H i gas mass fractions of f Hi = M Hi /(M ⋆ + M Hi ) = 90%, at z ≈ 7. Based on the [C ii] 158 μm luminosity function (LF) derived from the same sample of galaxies, we further place constraints on the cosmic H i gas mass density in galaxies (ρ Hi ) at this redshift, which we measure to be . This estimate is substantially lower by a factor of ≈10 than that inferred from an extrapolation of damped Lyα absorber (DLA) measurements and largely depends on the exact [C ii] LF adopted. However, we find this decrease in ρ Hi to be consistent with recent simulations and argue that this apparent discrepancy is likely a consequence of the DLA sight lines predominantly probing the substantial fraction of H i gas in high-z galactic halos, whereas [C ii] traces the H i in the ISM associated with star formation. We make predictions for this buildup of neutral gas in galaxies as a function of redshift, showing that at z ≳ 5, only ≈10% of the cosmic H i gas content is confined in galaxies and associated with the star-forming ISM.

Emergence of a new H i 21-cm absorption component at z ∼ 1.1726 towards the γ-ray blazar PKS 2355-106

ABSTRACT We report the emergence of a new H i 21-cm absorption at zabs = 1.172 635 in the damped Lyα absorber (DLA) towards the γ-ray blazar PKS 2355-106 (zem∼1.639) using science verification observations (2020 June) from the MeerKAT Absorption Line Survey (MALS). Since 2006, this DLA is known to show a narrow H i 21-cm absorption at zabs = 1.173019 coinciding with a distinct metal absorption-line component. We do not detect significant H i 21-cm optical depth variations from this known H i component. A high-resolution optical spectrum (2010 August) shows a distinct Mg i absorption at the redshift of the new H i 21-cm absorber. However, this component is not evident in the profiles of singly ionized species. We measure the metallicity ([Zn/H] = −(0.77 ± 0.11) and [Si/H]= −(0.96 ± 0.11)) and depletion ([Fe/Zn] = −(0.63 ± 0.16)) for the full system. Using the apparent column density profiles of Si ii, Fe ii, and Mg i, we show that the depletion and the N(Mg i)/N(Si ii) column density ratio systematically vary across the velocity range. The region with high depletion tends to have a slightly larger N(Mg i)/N(Si ii) ratio. The two H i 21-cm absorbers belong to this velocity range. The emergence of zabs = 1.172 635 can be understood if there is a large optical depth gradient over a length-scale of ∼0.35 pc. However, the gas producing the zabs = 1.173 019 component must be nearly uniform over the same scale. Systematic uncertainties introduced by the absorption-line variability has to be accounted for in experiments measuring the variations of fundamental constants and cosmic acceleration even when the radio emission is apparently compact as in PKS 2355-106.