Damped Lyman Alpha Research Articles

On the Nature of the Radio Calibrator and Gamma-Ray Emitting NLS1 Galaxy 3C 286 and Its Multiwavelength Variability

The quasar 3C 286, a well-known calibrator source in radio astronomy, was found to exhibit exceptional multiwavelength properties. Its rich and complex optical emission-line spectrum revealed its narrow-line Seyfert 1 (NLS1) nature. Given its strong radio emission, this makes 3C 286 one of the radio-loudest NLS1 galaxies known to date. 3C 286 is also one of very few known compact steep-spectrum (CSS) sources detected in the gamma-ray regime. Observations in the X-ray regime, rarely carried out so far, revealed evidence for variability, raising the question whether it is driven by the accretion disk or jet. 3C 286 is also well known for its damped Lyman alpha system from an intervening absorber at z = 0.692, triggering a search for the corresponding X-ray absorption along the line-of-sight. Here, we present new observations in the radio, X-ray, optical, and UV bands. The nature of the X-ray variability is addressed. Spectral evidence suggests that it is primarily driven by the accretion disk (not the jet), and the X-ray spectrum is well fit by a powerlaw plus soft excess model. The radio flux density and polarization remain constant at the Effelsberg telescope resolution, reconfirming the use of 3C 286 as radio calibrator. The amount of reddening/absorption along the line-of-sight intrinsic to 3C 286 is rigorously assessed. None is found, validating the derivation of a high Eddington ratio (L/LEdd ∼ 1) and of the very high radio-loudness index of 3C 286. Based on the first deep Chandra image of 3C 286, tentative evidence for hard X-ray emission from the SW radio lobe is reported. A large variety of models for the gamma-ray emission of 3C 286 are briefly discussed.

Neutral carbon in diffuse interstellar medium: abundance matching with H2 for damped Lyman alpha systems at high redshifts

ABSTRACT We present a study of C i/H2 relative abundance in the diffuse cold neutral medium (CNM). Using the chemical and thermal balance model, we calculate the dependence of C i/H2 on the main parameters of the medium: hydrogen number density, metallicity, strength of the UV field, and cosmic ray ionization rate (CRIR). We show that the observed relative C i and H2 column densities in damped Lyman alpha systems (DLAs) at high redshifts can be reproduced within our model assuming the typically expected conditions in the diffuse CNM. Using additional observed information on metallicity, H i column density, and excitation of C i fine-structure levels, as well as temperature, we estimated for a wide range metallicities in the CNM at high redshifts that CRIRs are in the range from ∼10−16 to a $\rm few \times 10^{-15}\, \rm s^{-1}$, hydrogen number densities are in the range ∼10−103 cm−3, and the UV field is in the range from 10−2 to a $\rm few \times 10^2$ of the Mathis field. We argue that because the observed quantities used in this work are quite homogeneous and much less affected by radiative transfer effects (in comparison with, for example, the dissociation of HD and UV pumping of H2 rotational levels), our estimates are quite robust against the assumption of the exact geometrical model of the cloud and local sources of the UV field.

Frequency–Redshift Relation of the Cosmic Microwave Background

We point out that a modified temperature–redshift relation (T-z relation) of the cosmic microwave background (CMB) cannot be deduced by any observational method that appeals to an a priori thermalisation to the CMB temperature T of the excited states in a probe environment of independently determined redshift z. For example, this applies to quasar-light absorption by a damped Lyman-alpha system due to atomic as well as ionic fine-splitting transitions or molecular rotational bands. Similarly, the thermal Sunyaev-Zel’dovich (thSZ) effect cannot be used to extract the CMB’s T-z relation. This is because the relative line strengths between ground and excited states in the former and the CMB spectral distortion in the latter case both depend, apart from environment-specific normalisations, solely on the dimensionless spectral variable x=hνkBT. Since the literature on extractions of the CMB’s T-z relation always assumes (i) ν(z)=(1+z)ν(z=0), where ν(z=0) is the observed frequency in the heliocentric rest frame, the finding (ii) T(z)=(1+z)T(z=0) just confirms the expected blackbody nature of the interacting CMB at z>0. In contrast to the emission of isolated, directed radiation, whose frequency–redshift relation (ν-z relation) is subject to (i), a non-conventional ν-z relation ν(z)=f(z)ν(z=0) of pure, isotropic blackbody radiation, subject to adiabatically slow cosmic expansion, necessarily has to follow that of the T-z relation T(z)=f(z)T(z=0) and vice versa. In general, the function f(z) is determined by the energy conservation of the CMB fluid in a Friedmann–Lemaitre–Robertson–Walker universe. If the pure CMB is subject to an SU(2) rather than a U(1) gauge principle, then f(z)=1/41/3(1+z) for z≫1, and f(z) is non-linear for z∼1.

A successful search for intervening 21 cm H i absorption in galaxies at 0.4 < z <1.0 with the Australian square kilometre array pathfinder (ASKAP)

ABSTRACT We have used the Australian Square Kilometre Array Pathfinder (ASKAP) radio telescope to search for intervening 21 cm neutral hydrogen (H i) absorption along the line of sight to 53 bright radio continuum sources. Our observations are sensitive to H i column densities typical of Damped Lyman Alpha absorbers (DLAs) in cool gas with an H i spin temperature below about 300–500 K. The six-dish Boolardy Engineering Test Array (BETA) and twelve-antenna Early Science array (ASKAP-12) covered a frequency range corresponding to redshift 0.4 < z < 1.0 and 0.37 < z < 0.77, respectively, for the H i line. Fifty of the 53 radio sources observed have reliable optical redshifts, giving a total redshift path Δz = 21.37. This was a spectroscopically untargeted survey, with no prior assumptions about the location of the lines in redshift space. Four intervening H i lines were detected, two of them new. In each case, the estimated H i column density lies above the DLA limit for H i spin temperatures above 50–80 K, and we estimate a DLA number density at redshift z ∼ 0.6 of $n(z)=0.19^{+0.15 }_{ -0.09}$. This value lies somewhat above the general trend of n(z) with redshift seen in optical DLA studies. Although the current sample is small, it represents an important proof of concept for the much larger 21 cm First Large Absorption Survey in H i (FLASH) project to be carried out with the full 36-antenna ASKAP telescope, probing a total redshift path $\Delta z\sim \, 50,000$.

The Sloan Digital Sky Survey Quasar Catalog: Sixteenth Data Release

We present the final Sloan Digital Sky Survey IV (SDSS-IV) quasar catalog from Data Release 16 of the extended Baryon Oscillation Spectroscopic Survey (eBOSS). This catalog comprises the largest selection of spectroscopically confirmed quasars to date. The full catalog includes two sub-catalogs: a "superset" of all SDSS-IV/eBOSS objects targeted as quasars containing 1,440,615 observations and a quasar-only catalog containing 750,414 quasars, including 225,082 new quasars appearing in an SDSS data release for the first time, as well as known quasars from SDSS-I/II/III. We present automated identification and redshift information for these quasars alongside data from visual inspections for 320,161 spectra. The quasar-only catalog is estimated to be 99.8% complete with 0.3% to 1.3% contamination. Automated and visual inspection redshifts are supplemented by redshifts derived via principal component analysis and emission lines. We include emission line redshifts for H$\alpha$, H$\beta$, Mg II, C III], C IV, and Ly$\alpha$. Identification and key characteristics generated by automated algorithms are presented for 99,856 Broad Absorption Line quasars and 35,686 Damped Lyman Alpha quasars. In addition to SDSS photometric data, we also present multi-wavelength data for quasars from GALEX, UKIDSS, WISE, FIRST, ROSAT/2RXS, XMM-Newton, and Gaia. Calibrated digital optical spectra for these quasars can be obtained from the SDSS Science Archive Server.

Deep Hubble Space Telescope Imaging on the Extended Lyα Emission of a QSO at z = 2.19 with a Damped Lyman Alpha System as a Natural Coronagraph

Abstract Recent observations suggest that proximate damped Lyα (PDLA) systems can be used to study the host galaxies of quasi-stellar objects (QSOs), because the PDLAs can block the bright point-spread function from central QSOs. Using six orbits of narrowband imaging with Hubble Space Telescope (HST)/WFC3, we present the first high-resolution narrowband image of the Lyα emission in the PDLA trough of the QSO SDSS J115432.67–021537.9. We detect one major component and one minor component in the narrowband imaging. Combining the /WFC3 imaging with deep Magellan/MagE spectra, we measure that the Lyα flux = , which is among the luminous (≈2.7 ) Lyα emitters at 2.19. The -based star formation rate is ∼7 . These observational results favor that the star formation from the host galaxy could be the main mechanism powering the Lyα emission. This new method sheds light on the study of the kinematic structure and the spatial distribution of the extended Lyα emitting regions around the QSO host.

Multiwavelength observations of GRB 140629A

Aims.We investigate the long gamma-ray burst (GRB) 140629A through multiwavelength observations to derive the properties of the dominant jet and its host galaxy.Methods.The afterglow and host galaxy observations were taken in the optical (Swift/UVOT and various facilities worldwide), infrared (Spitzer), and X-rays (Swift/XRT) between 40 s and 3 yr after the burst trigger.Results.Polarisation observations by the MASTER telescope indicate that this burst is weakly polarised. The optical spectrum contains absorption features, from which we confirm the redshift of the GRB as originating atz = 2.276 ± 0.001. We performed spectral fitting of the X-rays to optical afterglow data and find there is no strong spectral evolution. We determine the hydrogen column densityNHto be 7.2 × 1021cm−2along the line of sight. The afterglow in this burst can be explained by a blast wave jet with a long-lasting central engine expanding into a uniform medium in the slow cooling regime. At the end of energy injection, a normal decay phase is observed in both the optical and X-ray bands. An achromatic jet break is also found in the afterglow light curves ∼0.4 d after trigger. We fit the multiwavelength data simultaneously with a model based on a numerical simulation and find that the observations can be explained by a narrow uniform jet in a dense environment with an opening angle of 6.7° viewed 3.8° off-axis, which released a total energy of 1.4 × 1054erg. Using the redshift and opening angle, we find GRB 140629A follows both the Ghirlanda and Amati relations. From the peak time of the light curve, identified as the onset of the forward shock (181s after trigger), the initial Lorentz factor (Γ0) is constrained in the range 82–118. Fitting the host galaxy photometry, we find the host to be a low mass, star-forming galaxy with a star formation rate of log (SFR) 1.1+0.9−0.4M⊙yr−1. We obtain a value of the neutral hydrogen density by fitting the optical spectrum, logNHI = 21.0 ± 0.3, classifying this host as a damped Lyman-alpha. High ionisation lines (N V, Si IV) are also detected in the spectrum.

A metal-line strength indicator for damped Lyman alpha (DLA) systems at low signal-to-noise

The Baryon Oscillation Spectroscopic Survey of SDSS-III has provided an unprecedentedly large sample of Damped \lya systems (DLAs), the largest repositories of neutral hydrogen in the Universe. This DLA sample has been used to determine the DLA bias factor from their cross-correlation with the \lya forest absorption in \cite{FontRibera2012,Perez2018}, showing that DLAs are associated with relatively massive halos. However, the low resolution and signal-to-noise of BOSS spectra do not allow precise measurements of the DLA metal lines. We define a metal strength parameter, $S$, based on combining equivalent widths of 17 metal lines, that can be measured with an optimal signal-to-noise ratio for individual DLAs in BOSS spectra, allowing for the classification of the DLA population into subgroups of different $S$. We present the distribution of this DLA metal strength and the dependence of its mean value on $N_{\rm HI}$ and redshift. We search for systematic effects and variations in the catalogue purity by examining the dependence of the $S$ distribution on the spectral signal-to-noise and the estimated error on $S$. A catalogue of DLAs with measured equivalent widths for the selected 17 metal lines and the value of $S$ is made publicly available, which will be used to measure the dependence of the DLA bias factor on the $S$ parameter. The relation of the metal strength on the gas metal abundances and velocity dispersion can be constrained by studying the stacked metal absorption spectra of DLAs as a function of $S$, allowing for future determinations of the dependence of the bias factor on the metallicity and velocity dispersion of DLAs.

Multicomponent H2 in DLA at zabs = 2.05: physical conditions through observations and numerical models★

We perform detailed spectroscopic analysis and numerical modelling of an H2-bearing damped Lyman-alpha absorber (DLA) at zabs = 2.05 towards the quasar FBQS J2340-0053. Metal absorption features arise from fourteen components spread over $\Delta v_{90}$ = 114 km s$^{-1}$, seven of which harbour H2. Column densities of atomic and molecular species are derived through Voigt profile analysis of their absorption lines. We measure total N(H I), N(H2) and N(HD) to be 20.35+/-0.05, 17.99+/-0.05 and 14.28+/-0.08 (log cm$^{-2}$) respectively. H2 is detected in the lowest six rotational levels of the ground vibrational state. The DLA has metallicity, Z = 0.3 Z$_\sun$ ([S/H] = -0.52+/-0.06) and dust-to-gas ratio, $\kappa$ = 0.34+/-0.07. Numerical models of the H2 components are constrained individually to understand the physical structure of the DLA. We conclude that the DLA is subjected to the metagalactic background radiation and cosmic ray ionization rate of $\sim$ 10$^{-15.37}$ s$^{-1}$. Dust grains in this DLA are smaller than grains in the Galactic interstellar medium. The inner molecular regions of the H2 components have density, temperature and gas pressure in the range 30-120 cm$^{-3}$, 140-360 K and 7,000-23,000 cm$^{-3}$ K respectively. Micro-turbulent pressure is a significant constituent of the total pressure, and can play an important role in these innermost regions. Our H2 component models enable us to constrain component-wise N(H I), and elemental abundances of sulphur, silicon, iron and carbon. We deduce the line-of-sight thickness of the H2-bearing parts of the DLA to be 7.2 pc.

The cosmological bias factor of damped Lyman alpha systems: dependence on metal line strength

We measure the cosmological bias factor of DLAs from their cross-correlation with the Ly$\alpha$ forest absorption, as a function of the DLA metal strength, defined from an average of equivalent widths of the strongest detectable low-ionization metal lines. A clear increase of the bias factor with metal strength is detected, as expected from a relation of metallicity and velocity dispersion with host halo mass. The relation is stronger after the metal strength is corrected for the HI column density, to make it more related to metallicity instead of metal column density. After correcting for the effects of measurement errors of the metal strength parameter, we find that the bias factor of DLAs with the weakest metal lines is close to unity, consistent with an origin in dwarf galaxies with host halo masses $\sim10^{10} M_{\odot}$, whereas the most metal rich DLAs have a bias factor as large as $b_{\rm DLA} \sim 3$, indicative of massive galaxies or galaxy groups in host halos with masses $\sim 10^{12} M_{\odot}$. Our result confirms the physical origin of the relation of bias factors measured from cross-correlation studies to the host halos of the absorbers.

Multiple Absorption Systems in the Lines of Sight to Quadruply Lensed Quasar H1413+1143 As a Probe of the Circumgalactic medium around Dwarf Galaxy

AbstractWe present the first measurement of differences in MgII absorption strength in multiple intervening absorbers, which are also identified as (sub-)Damped Lyman alpha absorption systems, in the four spectra of the quadruply lensed quasar H1413+1143, often referred to the “Cloverleaf”, from highly spatial resolution and high signal-to-noise spectroscopy with an optical multi-mode spectrograph, the Kyoto tridimentional spectrograph II on board the Subaru telescope. The detection of significant MgII absorptions in multiple components in the spatially-resolved spectra suggests that chemical enrichment differs at least on scale of about 10 kpc within the separation of sightlines. For, a DLA system at redshift zabs = 1.66, the rest equivalent widths of MgII absorption lines change by factors up to 6, which is similar to those of HI absorption lines. This suggests that (inhomogeneous) cold absorbers which give rise to strong HI/MgII absorptions dwell on a scale within 10 kpc in the circumgalactic medium (CGM).

Oxygen and zinc abundances in 417 Galactic bulge red giants

Context. Oxygen and zinc in the Galactic bulge are key elements for the understanding of the bulge chemical evolution. Oxygen-to-iron abundance ratios provide a most robust indicator of the star formation rate and chemical evolution of the bulge. Zinc is enhanced in metal-poor stars, behaving as an α-element, and its production may require nucleosynthesis in hypernovae. Most of the neutral gas at high redshift is in damped Lyman-alpha systems (DLAs), where Zn is also observed to behave as an α-element. Aims. The aim of this work is the derivation of the α-element oxygen, together with nitrogen, and the iron-peak element zinc abundances in 417 bulge giants, from moderate resolution (R ~ 22 000) FLAMES-GIRAFFE spectra. For stars in common with a set of UVES spectra with higher resolution (R ~ 45 000), the data are intercompared. The results are compared with literature data and chemodynamical models. Methods. We studied the spectra obtained for a large sample of red giant stars, chosen to be one magnitude above the horizontal branch, using FLAMES-GIRAFFE on the Very Large Telescope. We computed the O abundances using the forbidden [OI] 6300.3 Å and Zn abundances using the Zn I 6362.34 Å lines. Stellar parameters for these stars were established in a previous work from our group. Results. We present oxygen abundances for 358 stars, nitrogen abundances for 403 stars and zinc abundances were derived for 333 stars. Having oxygen abundances for this large sample adds information in particular at the moderate metallicities of −1.6 < [Fe/H] < −0.8. Zn behaves as an α-element, very similarly to O, Si, and Ca. It shows the same trend as a function of metallicity as the α-elements, i.e., a turnover around [Fe/H] ~ − 0.6, and then decreasing with increasing metallicity. The results are compared with chemodynamical evolution models of O and Zn enrichment for a classical bulge. DLAs also show an enhanced zinc-to-iron ratio, suggesting they may be enriched by hypernovae.

Correlations in the three-dimensional Lyman-alpha forest contaminated by high column density absorbers

Correlations measured in three dimensions in the Lyman-alpha forest are contaminated by the presence of the damping wings of high column density (HCD) absorbing systems of neutral hydrogen (HI; having column densities $N(\mathrm{HI}) > 1.6\times10^{17}\,\mathrm{atoms}\,\mathrm{cm}^{-2}$), which extend significantly beyond the redshift-space location of the absorber. We measure this effect as a function of the column density of the HCD absorbers and redshift by measuring 3D flux power spectra in cosmological hydrodynamical simulations from the Illustris project. Survey pipelines exclude regions containing the largest damping wings. We find that, even after this procedure, there is a scale-dependent correction to the 3D Lyman-alpha forest flux power spectrum from residual contamination. We model this residual using a simple physical model of the HCD absorbers as linearly biased tracers of the matter density distribution, convolved with their Voigt profiles and integrated over the column density distribution function. We recommend the use of this model over existing models used in data analysis, which approximate the damping wings as top-hats and so miss shape information in the extended wings. The simple 'linear Voigt model' is statistically consistent with our simulation results for a mock residual contamination up to small scales ($|k| < 1\,h\,\mathrm{Mpc}^{-1}$). It does not account for the effect of the highest column density absorbers on the smallest scales (e.g., $|k| > 0.4\,h\,\mathrm{Mpc}^{-1}$ for small damped Lyman-alpha absorbers; HCD absorbers with $N(\mathrm{HI}) \sim 10^{21}\,\mathrm{atoms}\,\mathrm{cm}^{-2}$). However, these systems are in any case preferentially removed from survey data. Our model is appropriate for an accurate analysis of the baryon acoustic oscillations feature. It is additionally essential for reconstructing the full shape of the 3D flux power spectrum.

Baseline metal enrichment from Population III star formation in cosmological volume simulations

We utilize the hydrodynamic and N-body code GIZMO coupled with our newly developed sub-grid Population III (Pop III) Legacy model, designed specifically for cosmological volume simulations, to study the baseline metal enrichment from Pop III star formation at z > 7. In this idealized numerical experiment, we only consider Pop III star formation. We find that our model Pop III star formation rate density (SFRD), which peaks at ∼ 10− 3 M⊙ yr− 1 Mpc− 1 near z ∼ 10, agrees well with previous numerical studies and is consistent with the observed estimates for Pop II SFRDs. The mean Pop III metallicity rises smoothly from z = 25 to 7, but does not reach the critical metallicity value, Zcrit = 10−4 Z⊙, required for the Pop III to Pop II transition in star formation mode until z ≃ 7. This suggests that, while individual haloes can suppress in situ Pop III star formation, the external enrichment is insufficient to globally terminate Pop III star formation. The maximum enrichment from Pop III star formation in star-forming dark matter haloes is Z ∼ 10−2 Z⊙, whereas the minimum found in externally enriched haloes is Z ≳ 10−7 Z⊙. Finally, mock observations of our simulated IGM enriched with Pop III metals produce equivalent widths similar to observations of an extremely metal-poor damped Lyman alpha system at z = 7.04, which is thought to be enriched by Pop III star formation only.

Study of the 2H(p,γ)3He reaction in the BBN energy range at LUNA

Using Big Bang Nucleosynthesis with the recent cosmological parameters obtained by the Planck collaboration, a primordial deuterium abundance value D/H = (2.65 ± 0.07) × 10−5 is obtained. This one is a little bit in tension with astronomical observations on metal- poor damped Lyman alpha systems where D/H = (2.53 ± 0.04) × 105. In order to reduce the BBN calculation uncertainty, a measurement of the 2H(p,γ)3He cross section in the energy range 10-300 keV with a 3% accuracy is thus desirable. Thanks to the low background of the underground Gran Sasso Laboratories, and to the experience accumulated in more than twenty years of scientific activity, LUNA (Laboratory for Underground Nuclear Astrophysics) planned to measure the 2H(p,γ)3He fusion cross section at the BBN energy range in 2015-2016. A feasibility test of the measurement has been recently performed at LUNA. In this paper, the results obtained will be shown. Possible cosmological outcomes from the future LUNA data will be also discussed.

An HST/COS legacy survey of high-velocity ultraviolet absorption in the Milky Way’s circumgalactic medium and the Local Group

To characterize the absorption properties of this circumgalactic medium (CGM) and its relation to the LG we present the so-far largest survey of metal absorption in Galactic high-velocity clouds (HVCs) using archival ultraviolet (UV) spectra of extragalactic background sources. The UV data are obtained with the Cosmic Origins Spectrograph (COS) onboard the Hubble Space Telescope (HST) and are supplemented by 21 cm radio observations of neutral hydrogen. Along 270 sightlines we measure metal absorption in the lines of SiII, SiIII, CII, and CIV and associated HI 21 cm emission in HVCs in the velocity range |v_LSR|=100-500 km s^-1. With this unprecedented large HVC sample we were able to improve the statistics on HVC covering fractions, ionization conditions, small-scale structure, CGM mass, and inflow rate. For the first time, we determine robustly the angular two point correlation function of the high-velocity absorbers, systematically analyze antipodal sightlines on the celestial sphere, and compare the absorption characteristics with that of Damped Lyman alpha absorbers (DLAs) and constrained cosmological simulations of the LG. Our study demonstrates that the Milky Way CGM contains sufficient gaseous material to maintain the Galactic star-formation rate at its current level. We show that the CGM is composed of discrete gaseous structures that exhibit a large-scale kinematics together with small-scale variations in physical conditions. The Magellanic Stream clearly dominates both the cross section and mass flow of high-velocity gas in the Milky Way's CGM. The possible presence of high-velocity LG gas underlines the important role of the local cosmological environment in the large-scale gas-circulation processes in and around the Milky Way (abridged).

Neutral hydrogen in the post-reionization universe

AbstractThe evolution of neutral hydrogen (HI) across redshifts is a powerful probe of cosmology, large scale structure in the universe and the intergalactic medium. Using a data-driven halo model to describe the distribution of HI in the post-reionization universe (z ∼ 5 to 0), we obtain the best-fitting parameters from a rich sample of observational data: low redshift 21-cm emission line studies, intermediate redshift intensity mapping experiments, and higher redshift Damped Lyman Alpha (DLA) observations. Our model describes the abundance and clustering of neutral hydrogen across redshifts 0 - 5, and is useful for investigating different aspects of galaxy evolution and for comparison with hydrodynamical simulations. The framework can be applied for forecasting future observations with neutral hydrogen, and extended to the case of intensity mapping with molecular and other line transitions at intermediate redshifts.

The SDSS-DR12 large-scale cross-correlation of damped Lyman alpha systems with the Lyman alpha forest

We present a measurement of the DLA mean bias from the cross-correlation of DLA and the Ly$\alpha$ forest, updating earlier results of Font-Ribera et al. 2012 with the final BOSS Data Release and an improved method to address continuum fitting corrections. Our cross-correlation is well fitted by linear theory with the standard $\Lambda CDM$ model, with a DLA bias of $b_{\rm DLA} = 1.99\pm 0.11$; a more conservative analysis, which removes DLA in the Ly$\beta$ forest and uses only the cross-correlation at $r> 10{\rm h^{-1}\,Mpc}$, yields $b_{\rm DLA} = 2.00\pm 0.19$. This assumes the cosmological model from \cite{Planck2015} and the Ly$\alpha$ forest bias factors of Bautista et al. 2017, and includes only statistical errors obtained from bootstrap analysis. The main systematic errors arise from possible impurities and selection effects in the DLA catalogue, and from uncertainties in the determination of the Ly$\alpha$ forest bias factors and a correction for effects of high column density absorbers. We find no dependence of the DLA bias on column density or redshift. The measured bias value corresponds to a host halo mass $\sim 4\cdot10^{11} {\rm M_{\odot}}$ if all DLA were hosted in halos of a similar mass. In a realistic model where host halos over a broad mass range have a DLA cross section $\Sigma(M_h) \propto M_h^{\alpha}$ down to $M_h > M_{\rm min} =10^{8.5} {\rm M_{\odot}}$, we find that $\alpha > 1$ is required to have $b_{\rm DLA}> 1.7$, implying a steeper relation or higher value of $M_{\rm min}$ than is generally predicted in numerical simulations of galaxy formation.

The HighAVQuasar Survey: Az = 2.027 metal-rich damped Lyman-αabsorber towards a red quasar atz = 3.21

To fully exploit the potential of quasars as probes of cosmic chemical evolution and the internal gas dynamics of galaxies it is important to understand the selection effects behind the quasar samples and in particular if the selection criteria exclude foreground galaxies with certain properties (most importantly a high dust content). Here we present spectroscopic follow-up from the 10.4-m GTC telescope of a dust-reddened quasar, eHAQ0111+0641, from the extended High A_V Quasar (HAQ) survey. We find that the z=3.21 quasar has a foreground Damped Lyman-alpha Absorber (DLA) at z=2.027 along the line of sight. The DLA has very strong metal lines due to a moderately high metallicity (with an inferred lower limit of 25% of the solar metallicity), but a very large gas column density along the line-of-sight in its host galaxy. This discovery is further evidence that there is a dust bias affecting the census of metals, caused by the combined effect of dust obscuration and reddening, in existing samples of z>2 DLAs. The case of eHAQ0111+0641 illustrates that dust bias is not only caused by dust obscuration, but also dust reddening.

The statistical properties of neutral gas at z &lt; 1.65 from UV measurements of Damped Lyman Alpha systems

We derive the statistical properties of neutral gas at redshifts 0.11<z<1.65 from UV measurements of quasar Lyman alpha absorption lines corresponding to 369 MgII systems with $W^{\lambda2796}_{0} \ge 0.3$ \AA. In addition to the 41 damped Lyman alpha (DLA) systems presented in Rao et al. (2006), the current DLA sample includes 29 newly discovered DLAs. Of these, 26 were found in our HST ACS prism survey for DLAs (Turnshek et al. 2015) and three were found in a GALEX archival search. In addition, an HST COS Cycle 19 survey yielded no DLAs that could be used for this study. Formally, this DLA sample includes 70 systems with $N_{\rm HI}\ge 2\times 10^{20}$ atoms cm$^{-2}$. We find that the incidence of DLAs, or the product of their gas cross section and their comoving number density, can be described by $n_{\rm DLA}(z) = (0.027 \pm 0.007) (1+z)^{(1.682 \pm 0.200)}$ over the redshift range 0<z<5. The cosmic mass density of neutral gas can be described by $\Omega_{\rm DLA}(z) = (4.77 \pm 1.60)\times10^{-4} (1 + z)^{(0.64\pm 0.27)}$. The low-redshift column density distribution function is well-fitted by a power law of the form $f(N) \sim N^\beta$ with $\beta = -1.46 \pm 0.20$. It is consistent with the high-redshift as well as z=0 estimates at the high column density end but, lies between them at the low column density end. We discuss possible $N_{\rm HI}$ and metallicity bias in MgII-selected DLA samples and show that such biases do not exist in the current data at z<1.65. Thus, at least at z<1.65, DLAs found through MgII selection statistically represent the true population of DLAs. However, we caution that studies of DLA metallicities should take into the account the relative incidence of DLAs with respect to $W^{\lambda2796}_{0}$ (or gas velocity spread) in order to correctly measure the mean neutral-gas cosmic metallicity of the universe.