# Neutral atomic-carbon quasar absorption-line systems at $z> 1.5$. Sample selection, H i content, reddening, and 2175 Å extinction feature

C. Ledoux, P. Noterdaeme, P. Petitjean, R. Srianand

► To cite this version:

C. Ledoux, P. Noterdaeme, P. Petitjean, R. Srianand. Neutral atomic-carbon quasar absorption-line systems at $z> 1.5$. Sample selection, H i content, reddening, and 2175 Å extinction feature. Astronomy and Astrophysics - A&A, 2015, 580, 10.1051/0004-6361/201424122. insu-03644896

HAL Id: insu-03644896

https://insu.hal.science/insu-03644896

Submitted on 20 Apr 2022

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Neutral atomic-carbon quasar absorption-line systems at $z > 1.5$
## Sample selection, H1 content, reddening, and 2175 Å extinction feature*

C. Ledoux<fnref n="1" />, P. Noterdaeme<fnref n="2" />, P. Petitjean<fnref n="2" />, and R. Srianand<fnref n="3" />

Received 2 May 2014 / Accepted 20 April 2015

## ABSTRACT

We present the results of a search for cold gas at high redshift along quasar lines of sight carried out without any a priori assumption on the neutral atomic-hydrogen content of the absorption-line systems. To do this, we systematically looked for neutral-carbon (C I) λλ1560, 1656 transition lines in 41 696 low-resolution quasar spectra ($1.5 < z_{em} < 4.46$) from the SDSS-II - Data Release 7 - database. C I absorption lines should indeed probe the shielded gas in the neutral interstellar medium of galaxies more efficiently than traditional tracers such as neutral atomic-hydrogen (H I) damped Lyman-$\alpha$ (DLA) and/or Mg II systems. We built up a sample of 66 C I absorbers with redshifts in the range $1.5 < z < 3.1$ and rest-frame equivalent widths $0.1 < W_r(\lambda\text{1560}) < 1.7$ Å. The completeness limit of our survey is $W_{r,\text{lim}}(\lambda\text{1560}) \approx 0.4$ Å. C I systems stronger than that are more than one hundred times rarer than DLAs at $z_{abs} = 2.5$. The number of C I systems per unit redshift is found to increase significantly below $z = 2$. We suggest that these C I absorbers are closely related to the process of star formation and the production of dust in galaxies and that their cosmic evolution is driven by the interplay between dust shielding and the evolution of the ultra-violet background at ~10 eV. We derive the neutral atomic-hydrogen content of the C I systems observable from the southern hemisphere from VLT/UVES spectroscopy and find that a majority of them are sub-DLAs with $N(H\text{I}) \sim 10^{20}$ atoms cm$^{-2}$. The dust content of these absorbers is yet significant as seen from the redder optical colours of the corresponding background quasars and their reddened spectral energy distributions, with $E(B-V)$ values up to ~0.3. The overall $N(H\text{I})$ distribution of C I systems is, however, relatively flat. As a consequence, among the C I systems classifying as DLAs, there is a probable excess of strong DLAs with $\log N(H\text{I}) > 21$ (atoms cm$^{-2}$) compared to systematic DLA surveys. Whilst the dust content of such systems is significant, their dust-to-gas ratio must still be limited. Indeed, strong DLAs having large amounts of shielded gas and dust producing stronger reddening and extinction of the background quasars, if they exist, should have been missed in the current magnitude-limited quasar sample. We study the empirical relations between $W_r(C\text{I})$, $N(H\text{I})$, $E(B-V)$, and the strength of the 2175 Å extinction feature. The latter is detected in about 30% of the C I absorbers. We show that, at equal amount of reddening, the 2175 Å feature is weak compared to Galactic lines of sight. This is probably the consequence of current or past star formation in the vicinity of the C I systems. We also find that the strongest C I systems tend to have the largest amounts of dust and that the metallicity of the gas and its molecular fraction is likely to be high in a large number of cases. The C I-absorber sample presented here hence provides ideal targets for detailed studies of the dust composition and molecular species at high redshift.

**Key words.** cosmology: observations – quasars: absorption lines – galaxies: ISM – dust, extinction

## 1. Introduction

Understanding the mechanisms of star formation at high redshifts is central to our knowledge of how galaxies formed and subsequently evolved chemically. This is especially true at $z \sim 2$ when the cosmic star-formation activity was highest (Madau & Dickinson 2014). Stars form out of cold gas, metals, and dust in molecular clouds (e.g. Snow & McCall 2006) in the interstellar medium (ISM) of galaxies. In turn, the radiative and mechanical feedback from stars has a strong impact on the physical state of the ISM. Studying the ISM at high redshifts and, in

particular, deriving the physical properties of the diffuse molecular phase in galaxies are therefore crucial for understanding how stars formed in the early Universe.

The best way to derive physical properties accurately is to detect the tracers of the cold gas in absorption (see Muller et al. 2014). The neutral, shielded, and possibly cold gas clouds at high redshifts can be searched for in the radio domain by targeting the neutral atomic-hydrogen (H I) 21-cm absorption line (e.g., Gupta et al. 2009). However, systematic – blind – surveys have to await the increased sensitivity of new facilities, such as MeerKAT/SKA and ASKAP (Booth et al. 2009; Duffy et al. 2012). On the other hand, this gas can be efficiently traced in the optical wave bands by detecting the redshifted H I damped Lyman-$\alpha$ (Prochaska et al. 2005; Noterdaeme et al. 2009b, 2012) and/or strong Mg II (e.g., Quider et al. 2011; Budzynski & Hewett 2011) lines imprinted in the spectra of bright enough background sources, such as quasars (also

<footnote n="1">¹ European Southern Observatory, Alonso de Córdova 3107, Casilla 19001, Vitacura, Santiago 19, Chile
e-mail: cledoux@eso.org</footnote>

<footnote n="2">² Institut d'Astrophysique de Paris, CNRS and UPMC Paris 6, UMR 7095, 98bis boulevard Arago, 75014 Paris, France</footnote>

<footnote n="3">³ Inter-University Centre for Astronomy and Astrophysics, Post Bag 4, Ganeshkhind, 411007 Pune, India</footnote>

<footnote>* Based on data from the Sloan Digital Sky Survey (SDSS) and dedicated follow-up observations carried out at the European Southern Observatory (ESO) under programmes 082.A-0544 and 083.A-0454 (P.I. Ledoux) using the Ultra-violet and Visual Echelle Spectrograph (UVES) installed at the Nasmyth-B focus of the Very Large Telescope (VLT), Unit-2 (Kueyen), on Cerro Paranal, Chile.</footnote>
referred to as quasi-stellar objects – QSOs) or the rapidly fad-
ing γ-ray burst (GRB) afterglows (for the latter, see Fynbo et al.
2009, and references therein).

Damped Lyman-α systems (hereafter DLAs) observed in
QSO spectra have column densities of N(HI) ≥ 2 ×
10²⁰ atoms cm⁻² and are known to contain most of the neu-
tral gas in the Universe in the redshift range 0 < z < 5 (see
Wolfe et al. 2005, for a review). It has been shown, however,
that DLAs typically probe warm (T ≥ 3000 K) and diffuse
(n_H < 1 cm⁻³) neutral gas (e.g., Petitjean et al. 2000; Srianand
et al. 2012). The metallicity of DLAs is generically low, i.e.,
on average about 1/30th of solar (Pettini 2006; Rafelski et al.
2012), and their dust-to-gas ratio is typically less than one-
tenth of what is observed in the Galactic ISM (e.g., Vladilo
et al. 2008). This probably explains the low detection rates of
molecular hydrogen (H₂) in DLAs where only about 10% of the
QSO lines of sight intercept H₂-bearing gas down to a limit of
N(H₂) ~ 10¹⁴ molecules cm⁻² (e.g., Noterdaeme et al. 2008a;
Balashev et al. 2014; for searches for H₂ in DLAs originating in
the host galaxies of GRBs, see Ledoux et al. 2009; Krühler et al.
2013).

Based on the observed correlation between metallicity and
dust depletion in DLAs (Ledoux et al. 2003), DLAs with high
metallicity are expected to contain more dust and therefore larger
amount of H₂ (see Petitjean et al. 2006). However, even in DLAs
with the highest metallicities, typical dust signatures, such as
reddening of the background QSOs, the 2175 Å extinction fea-
ture (hereafter also called ultraviolet [UV] bump) or diffuse in-
terstellar bands, are not apparent. In the rare cases with H₂ de-
tections, the inferred molecular fractions are low and typical
of what is seen in Galactic diffuse atomic gas with f(H₂) ≡
2N(H₂)/[2N(H₂) + N(HI)] ≤ 0.01 and often much lower than
this (see Noterdaeme et al. 2008a). The primary reason for this
is that the cold and dusty phases are missed probably because of
their reduced cross-sections relative to the more pervasive warm
neutral ISM (Zwaan & Prochaska 2006). Direct evidence for the
relatively small physical sizes (≤0.15 pc) of H₂-detected clouds
in DLAs has recently come from the observation of partial cover-
age of the QSO broad-line emitting region (Balashev et al. 2011;
Klimenko et al. 2015).

The clouds where H₂ is detected in DLAs are found to have
kinetic temperatures in the range T ~ 70–200 K and particle
densities n<sub>H</sub> ~ 1–100 cm<sup>-3</sup> (e.g., Srianand et al. 2005). When
detected, H<sub>2</sub> is usually coincident with neutral atomic carbon
(CI; see also Ge & Bechtold 1999). This is because the ion-
ization potential of neutral carbon (11.26 eV) is similar to the
average energy of the Lyman-Werner photons that dissociate
H<sub>2</sub>. Therefore, shielding of UV photons is essential for these
species to remain at detectable levels. Carbon monoxide (CO)
has long escaped detection even in DLAs with detected H<sub>2</sub>,
down to N(CO) ~ 10<sup>12</sup> molecules cm<sup>-2</sup> (see, e.g., Petitjean et al.
2002). This is not surprising since with a dissociation energy of
11.09 eV, CO needs to be even more shielded than H<sub>2</sub> and CI to
be detected. After UV absorption bands from CO were detected
for the first time at high redshift in a sub-DLA towards the QSO
SDSS J143912.05+111740.6 (Srianand et al. 2008b), it became
clear that the best places to detect CO in absorption at high red-
shift are the systems with strong CI absorption. Following this
strategy allowed us to detect carbon monoxide in five additional
systems (Noterdaeme et al. 2009a, 2010a, 2011).

In Galactic translucent interstellar clouds, CO starts to be
produced in significant amounts when neutral atomic carbon be-
comes the dominant carbon species and a large portion of hy-
drogen turns molecular (Snow & McCall 2006). The strength

of the C I absorption is expected to be such that it could be de-
tected even in a low-resolution spectrum. We therefore embarked
on a systematic search for C I absorption in QSO spectra from
the SDSS-II – Data Release 7 (hereafter DR 7) – database. In
this paper, we present the results of this search and the basic
properties of the detected C I absorbers. Here we equally refer to
QSO absorption-line systems detected through C I absorption as
“C I systems” or “C I absorbers”.

In Sect. 2, we describe our selection and identification of C I absorbers. We discuss the properties of the sample in terms of intervening C I-absorber number per unit redshift, proximate systems, and C I rest-frame equivalent widths in Sects. 3.3, 3.4, and 3.1, respectively. We then assess the impact of the C I absorbers on their respective background QSOs both from the observed QSO optical colours (Sect. 4.1) and the reddening these systems induce on the QSO spectral energy distributions (Sect. 4.2). In Sect. 5, we present the H I column-density distribution of the C I systems from VLT/UVES spectroscopy. In Sect. 6, we investigate empirical relations between neutral atomic-carbon and neutral atomic-hydrogen contents, QSO reddening, and the strength of possible 2175 Å extinction features (whose measurements are described in Sect. 4.3). We summarize our findings and conclude in Sect. 7.

Throughout this paper, we assume a standard Λ cold dark-matter cosmology with H₀ = 70 km s⁻¹ Mpc⁻¹, ΩΛ = 0.7 and ΩM = 0.3.

## 2. C I absorption-line selection and identification

We systematically searched for C I absorption lines in high-
redshift QSO spectra from the Sloan Digital Sky Survey (York
et al. 2000) – DR 7 (Abazajian et al. 2009) – quasar catalogue
(Schneider et al. 2010). This survey imposed an i-band magni-
tude cut of 19.1 for QSO candidates whose colours indicate a
probable redshift lower than ~3. The spectra cover the wave-
length range 3800–9200 Å at a resolving power R ~ 2000.

We implemented a dedicated IDL procedure to detect and
identify absorption-line features in SDSS QSO spectra auto-
matically. Since SDSS spectra are log-lambda-binned, the pix-
els have a constant velocity size (≈69 km s⁻¹). This makes it
straightforward to cross-correlate the spectra with an emission-
or absorption-line template. We used the method introduced by
Noterdaeme et al. (2010b) to search for [O III] λλ4959, 5007
emission and Mg II absorption lines. We first normalized the
spectra iteratively using Savitzky-Golay filtering. This consists
of smoothing the spectra by convolving them with a Savitsky-
Golay kernel that preserves the sharp QSO emission-line peaks
but ignores narrow features, such as metal absorption lines, bad
CCD pixels, and sky emission-line residuals. Deviant pixels and
their neighbours are then masked out, and the resulting data
is convolved again in the same way, and so on and so forth.
This procedure has the most advantages because no assump-
tion is required about the functional form of the QSO contin-
uum (i.e., power law or other), and in addition it is extremely
fast computationally. We then cross-correlated the normalized
spectra with a synthetic profile of C I λλ1560, 1656 absorption
lines. We looked for the positive correlation signal, together with
peak absorptions detected at more than 2σ and 2.5σ, and dif-
fering by less than a factor of three. The simultaneous detec-
tions of the Si II λ1526 and Al II λ1670 absorption lines were re-
quired to support the identifications of the two features as C I,
hence minimize the probability of coincidence. Spurious detec-
tions (~50%), most of them close to the detection limit, were
<figure><img src="image_2.png" /></figure>

**Fig. 1.** SDSS spectrum of the $z_{em} = 1.94$ QSO J0815+2640. The C I $\lambda\lambda$1560,1656 absorption lines detected at $z = 1.681$ are indicated, along with the presence of several low-ionization metal lines at the same redshift. Among the latter, the simultaneous detections of Si II $\lambda$1526 and Al II $\lambda$1670 lines were required by our detection algorithm to minimize the probability of coincidence. In this system, relatively weak Fe II $\lambda$1608 absorption is also observed.

identified visually and removed from the sample. In total, we find 66 systems, one of which is shown in Fig. 1.

The search for C I lines was limited to the regions of the spectra redwards of the QSO Lyman-$\alpha$ emission line to avoid the spurious coincidences that are frequent in the Lyman-$\alpha$ forest. The wavelength range above 7200 Å was also not considered in the search to avoid regions of the spectra that are heavily affected by residuals from the sky emission-line subtraction. We requested that the search window encompasses the wavelengths of the Si II $\lambda$1526 and Al II $\lambda$1670 absorption lines of the putative C I systems so that the validity of a system does not rely solely on the detection of two transitions (see above). For a given line of sight, the redshift lower bound ($z_{min}$) of the C I search is therefore the highest value between $z = 3820/1526 - 1 \approx 1.50$ and $(1 + z_{em}) \times 1215/1526 - 1 \approx 0.8 \times z_{em} - 0.2$, where $z_{em}$ is the QSO emission redshift. The redshift upper bound ($z_{max}$) of the search is the lowest value between $z_{em} + 0.1$ (to not exclude a priori proximate systems with infalling velocities of up to +5000 km s$^{-1}$) and $z = 7200/1656 - 1 \approx 3.35$. In order to avoid too many false positives at a low signal-to-noise ratio (S/N), we requested the median S/N per pixel to be higher than four for a given spectrum to be actually scanned. This resulted in a sample of 41 696 QSOs with $1.5 < z_{em} < 4.46$ whose spectra were searched for intervening or proximate C I absorbers. We did not initially reject broad absorption-line (BAL) quasars because our procedure follows the QSO continuum locally and can detect narrow absorption lines embedded in broad and not-fully-saturated troughs. Regions of deep absorption are avoided de facto when we study the number of intervening C I absorbers per unit redshift in Sect. 3.3 because they have low S/N per pixel.

Table 1 summarises our C I sample, which we refer to in the following as the overall sample. QSO names with J2000 coordinates are given for each absorber. No line of sight is found to feature more than one system<fnref n="1" />. The SDSS plate and fibre numbers, as well as the MJD, are also provided in the table as useful cross-references. The QSO emission redshifts derived by the SDSS team are indicated, together with the absorption redshifts and rest-frame equivalent widths of the C I lines. The last were

<figure><img src="image_3.png" /></figure>

**Fig. 2.** Rest-frame equivalent widths of C I absorption lines measured in SDSS spectra and associated $1\sigma$ uncertainties. The upper dashed line shows the 1:1 relation expected in case of heavily saturated profiles. The lower dashed line shows the expectation for absorption lines on the linear part of the curve of growth.

carefully determined by us for each individual system. In the next column of the table, we specify the average S/Ns per pixel in the regions the two C I lines are detected. At $z_{abs} > 2.2$, the Lyman-$\alpha$ line of the systems is also covered by the SDSS spectra. We therefore provide in the table a determination of the total neutral atomic-hydrogen column density of these systems following the method developed by Noterdaeme et al. (2009b). The reliability of the latter is confirmed in a number of cases with follow-up high-resolution VLT/UVES spectroscopy (see the last column of Table 1, and Sect. 5). This column also has $N(H\alpha)$ measurements from UVES for several low-redshift systems for which Lyman-$\alpha$ absorption is not covered in the SDSS spectrum.

## 3. Sample properties

It must be noted that in this survey C I systems are found without any presumptions on the presence of neutral atomic hydrogen; i.e., the C I systems found in this work do not necessarily have to be DLAs. Moreover, DLA absorbers can be observed in SDSS spectra only when their redshifts are higher than 2.2, while C I lines can be identified down to $z_{abs} \approx 1.50$. As stated before, the reality of the identified C I systems in our sample was checked by visual inspection. Therefore, we believe the C I detections are secure. We discuss the completeness of our survey in Sect. 3.2.

### 3.1. Line equivalent widths

Because we rely on them in the analysis, we here seek to verify the robustness and accuracy of the equivalent-width measurements of C I-absorption lines performed in SDSS spectra. For this purpose, we plot in Fig. 2 the measured rest-frame equivalent widths of the C I $\lambda\lambda$1560, 1656 lines versus each other. It appears that except in one case, the strengths of the two lines

<footnote n="1">¹ There is, however, a second C I system towards SDSS J234023.67-005327.1, which with $z_{abs} = 1.36$, falls below the redshift cut-off of our survey. This system happens to be detected in 21-cm absorption so is also related to cold gas (see Gupta et al. 2009; Kanekar et al. 2010).</footnote>
<table><caption>**Table 1.** C I systems identified in SDSS-DR7 QSO spectra.</caption><thead><tr><th>QSO name</th><th>Plate</th><th>MJD</th><th>Fibre</th><th>z<sub>em</sub></th><th>z<sub>abs</sub></th><th>W<sub>r</sub>(λ1560)<sup>a</sup><br>[Å]</th><th>W<sub>r</sub>(λ1656)<sup>a</sup><br>[Å]</th><th>S/N<sup>b</sup></th><th>log N(HI)<sub>SDSS</sub><sup>c</sup><br>[atoms cm<sup>-2</sup>]</th><th>log N(HI)<sub>UVES</sub><sup>d</sup><br>[atoms cm<sup>-2</sup>]</th></tr></thead><tbody><tr><td>J021606.13-002104.9</td><td>0405</td><td>51816</td><td>114</td><td>2.22</td><td>1.737</td><td>0.33 ± 0.10</td><td>0.77 ± 0.08</td><td>14</td><td>...</td><td>...</td></tr><tr><td>J030049.24-072137.8</td><td>0458</td><td>51929</td><td>437</td><td>2.11</td><td>1.536</td><td>0.79 ± 0.15</td><td>1.17 ± 0.21</td><td>5</td><td>...</td><td>...</td></tr><tr><td>J080801.74+330009.2</td><td>0825</td><td>52289</td><td>166</td><td>1.90</td><td>1.888</td><td>0.24 ± 0.05</td><td>0.39 ± 0.05</td><td>18</td><td>...</td><td>...</td></tr><tr><td>J081116.06+083837.7</td><td>2570</td><td>54081</td><td>566</td><td>2.07</td><td>1.906</td><td>0.23 ± 0.06</td><td>0.21 ± 0.06</td><td>15</td><td>...</td><td>20.10 ± 0.10</td></tr><tr><td>J081540.60+264021.6</td><td>1266</td><td>52709</td><td>355</td><td>1.94</td><td>1.681</td><td>1.70 ± 0.21</td><td>1.33 ± 0.16</td><td>5</td><td>...</td><td>...</td></tr><tr><td>J082003.40+155932.9</td><td>2272</td><td>53713</td><td>378</td><td>1.95</td><td>1.547</td><td>0.72 ± 0.15</td><td>0.77 ± 0.11</td><td>7</td><td>...</td><td>...</td></tr><tr><td>J085206.65+193548.4</td><td>2281</td><td>53711</td><td>452</td><td>1.91</td><td>1.788</td><td>0.34 ± 0.06</td><td>0.86 ± 0.07</td><td>13</td><td>...</td><td>20.00 ± 0.20</td></tr><tr><td>J085437.60+031734.8</td><td>0565</td><td>52225</td><td>480</td><td>2.24</td><td>1.567</td><td>0.40 ± 0.09</td><td>0.43 ± 0.07</td><td>12</td><td>...</td><td>...</td></tr><tr><td>J085726.79+185524.3</td><td>2281</td><td>53711</td><td>630</td><td>1.89</td><td>1.730</td><td>0.20 ± 0.02</td><td>0.43 ± 0.04</td><td>33</td><td>...</td><td>19.60 ± 0.20</td></tr><tr><td>J090558.75+553039.2</td><td>0450</td><td>51908</td><td>470</td><td>1.57</td><td>1.579</td><td>0.13 ± 0.05</td><td>0.15 ± 0.05</td><td>16</td><td>...</td><td>...</td></tr><tr><td>J090942.56+532126.5</td><td>0553</td><td>51999</td><td>425</td><td>2.07</td><td>1.628</td><td>0.47 ± 0.08</td><td>0.70 ± 0.09</td><td>12</td><td>...</td><td>...</td></tr><tr><td>J091516.27+071750.9</td><td>1194</td><td>52703</td><td>506</td><td>2.38</td><td>1.681</td><td>0.50 ± 0.10</td><td>0.95 ± 0.13</td><td>11</td><td>...</td><td>...</td></tr><tr><td>J091721.37+015448.1</td><td>0473</td><td>51929</td><td>426</td><td>2.18</td><td>2.107</td><td>0.67 ± 0.16</td><td>0.93 ± 0.16</td><td>5</td><td>...</td><td>...</td></tr><tr><td>J092759.79+154321.8</td><td>2579</td><td>54068</td><td>364</td><td>1.81</td><td>1.731</td><td>0.61 ± 0.10</td><td>0.69 ± 0.09</td><td>9</td><td>...</td><td>21.10 ± 0.20</td></tr><tr><td>J094529.37+552525.7</td><td>0556</td><td>51991</td><td>045</td><td>2.24</td><td>1.867</td><td>0.74 ± 0.12</td><td>0.91 ± 0.13</td><td>7</td><td>...</td><td>...</td></tr><tr><td>J095114.97+545736.5</td><td>0769</td><td>52282</td><td>458</td><td>1.80</td><td>1.613</td><td>0.77 ± 0.07</td><td>0.78 ± 0.06</td><td>14</td><td>...</td><td>...</td></tr><tr><td>J101544.34+310617.2</td><td>1953</td><td>53358</td><td>604</td><td>1.56</td><td>1.596</td><td>0.20 ± 0.09</td><td>0.34 ± 0.09</td><td>9</td><td>...</td><td>...</td></tr><tr><td>J104705.76+205734.5</td><td>2478</td><td>54097</td><td>558</td><td>2.01</td><td>1.775</td><td>0.59 ± 0.14</td><td>0.96 ± 0.13</td><td>7</td><td>...</td><td>...</td></tr><tr><td>J105436.96+542343.7</td><td>0907</td><td>52373</td><td>129</td><td>1.50</td><td>1.513</td><td>0.49 ± 0.08</td><td>0.87 ± 0.08</td><td>10</td><td>...</td><td>...</td></tr><tr><td>J105746.42+662250.5</td><td>0490</td><td>51929</td><td>179</td><td>2.34</td><td>2.088</td><td>0.30 ± 0.06</td><td>0.48 ± 0.06</td><td>13</td><td>...</td><td>...</td></tr><tr><td>J111756.53+143716.9</td><td>1753</td><td>53383</td><td>353</td><td>2.09</td><td>2.001</td><td>0.67 ± 0.09</td><td>0.65 ± 0.08</td><td>10</td><td>...</td><td>19.80 ± 0.15</td></tr><tr><td>J112250.95+143732.5</td><td>1753</td><td>53383</td><td>477</td><td>2.05</td><td>1.554</td><td>0.17 ± 0.04</td><td>0.19 ± 0.03</td><td>28</td><td>...</td><td>...</td></tr><tr><td>J112913.29-023740.9</td><td>0327</td><td>52294</td><td>308</td><td>1.86</td><td>1.623</td><td>0.60 ± 0.10</td><td>0.99 ± 0.12</td><td>8</td><td>...</td><td>...</td></tr><tr><td>J113341.29-005740.1</td><td>0282</td><td>51658</td><td>215</td><td>1.68</td><td>1.706</td><td>1.36 ± 0.25</td><td>1.24 ± 0.24</td><td>4</td><td>...</td><td>...</td></tr><tr><td>J114141.84+444206.1</td><td>1367</td><td>53083</td><td>499</td><td>1.96</td><td>1.903</td><td>0.48 ± 0.11</td><td>0.76 ± 0.13</td><td>7</td><td>...</td><td>...</td></tr><tr><td>J115153.87+150945.0</td><td>1762</td><td>53415</td><td>495</td><td>3.05</td><td>2.400</td><td>0.36 ± 0.11</td><td>0.37 ± 0.11</td><td>7</td><td>20.15 ± 0.33<sup>e</sup></td><td>20.15 ± 0.10</td></tr><tr><td>J115705.51+615521.7</td><td>0777</td><td>52320</td><td>107</td><td>2.51</td><td>2.460</td><td>1.22 ± 0.14</td><td>1.30 ± 0.14</td><td>7</td><td>21.80 ± 0.20<sup>f</sup></td><td>...</td></tr><tr><td>J120935.79+671715.7</td><td>0493</td><td>51957</td><td>178</td><td>2.03</td><td>1.843</td><td>0.47 ± 0.11</td><td>0.48 ± 0.11</td><td>8</td><td>...</td><td>...</td></tr><tr><td>J122825.67+303038.6</td><td>2234</td><td>53823</td><td>430</td><td>3.91</td><td>3.098</td><td>0.17 ± 0.07</td><td>0.26 ± 0.07</td><td>11</td><td>21.05 ± 0.21<sup>e</sup></td><td>...</td></tr><tr><td>J123714.61+064759.6</td><td>1628</td><td>53474</td><td>193</td><td>2.78</td><td>2.691</td><td>0.37 ± 0.06</td><td>0.77 ± 0.08</td><td>14</td><td>20.15 ± 0.28<sup>e</sup></td><td>20.00 ± 0.15<sup>g</sup></td></tr><tr><td>J124708.43+500320.8</td><td>1278</td><td>52735</td><td>039</td><td>2.27</td><td>2.135</td><td>1.05 ± 0.10</td><td>1.42 ± 0.11</td><td>9</td><td>...</td><td>...</td></tr><tr><td>J124829.71+284858.1</td><td>2239</td><td>53726</td><td>219</td><td>1.54</td><td>1.513</td><td>0.32 ± 0.07</td><td>0.65 ± 0.07</td><td>12</td><td>...</td><td>...</td></tr><tr><td>J124841.45+302433.0</td><td>2239</td><td>53726</td><td>469</td><td>2.06</td><td>1.691</td><td>0.57 ± 0.06</td><td>0.87 ± 0.06</td><td>15</td><td>...</td><td>...</td></tr><tr><td>J125552.60+223424.4</td><td>2649</td><td>54212</td><td>507</td><td>1.82</td><td>1.526</td><td>0.25 ± 0.11</td><td>1.18 ± 0.28</td><td>8</td><td>...</td><td>...</td></tr><tr><td>J130225.28+211158.6</td><td>2650</td><td>54505</td><td>296</td><td>1.76</td><td>1.656</td><td>0.56 ± 0.12</td><td>0.49 ± 0.10</td><td>8</td><td>...</td><td>...</td></tr><tr><td>J130628.87+281550.8</td><td>2009</td><td>53904</td><td>288</td><td>2.10</td><td>2.012</td><td>0.28 ± 0.10</td><td>0.56 ± 0.10</td><td>8</td><td>...</td><td>...</td></tr><tr><td>J130828.43+584000.6</td><td>0958</td><td>52410</td><td>224</td><td>3.09</td><td>2.473</td><td>0.39 ± 0.07</td><td>0.61 ± 0.14</td><td>11</td><td>20.00 ± 0.15</td><td>...</td></tr><tr><td>J131129.11+222552.6</td><td>2651</td><td>54507</td><td>398</td><td>3.14</td><td>3.093</td><td>0.33 ± 0.05</td><td>0.51 ± 0.08</td><td>16</td><td>20.50 ± 0.10</td><td>20.62 ± 0.10</td></tr><tr><td>J131400.57+054319.5</td><td>0850</td><td>52338</td><td>620</td><td>1.89</td><td>1.583</td><td>0.23 ± 0.10</td><td>0.70 ± 0.10</td><td>8</td><td>...</td><td>...</td></tr><tr><td>J134122.51+185214.0</td><td>2642</td><td>54232</td><td>494</td><td>2.00</td><td>1.544</td><td>0.16 ± 0.03</td><td>0.15 ± 0.04</td><td>27</td><td>...</td><td>...</td></tr><tr><td>J134601.10+064408.4</td><td>1803</td><td>54152</td><td>074</td><td>2.09</td><td>1.512</td><td>0.79 ± 0.20</td><td>1.33 ± 0.28</td><td>6</td><td>...</td><td>...</td></tr><tr><td>J135122.00+461828.5</td><td>1466</td><td>53083</td><td>416</td><td>1.81</td><td>1.606</td><td>0.99 ± 0.15</td><td>1.58 ± 0.17</td><td>10</td><td>...</td><td>...</td></tr><tr><td>J141550.47+300146.9</td><td>2129</td><td>54252</td><td>288</td><td>2.08</td><td>1.676</td><td>1.35 ± 0.10</td><td>1.61 ± 0.12</td><td>13</td><td>...</td><td>...</td></tr><tr><td>J141606.79+180403.2</td><td>2759</td><td>54534</td><td>359</td><td>2.13</td><td>1.622</td><td>0.47 ± 0.11</td><td>0.80 ± 0.13</td><td>8</td><td>...</td><td>...</td></tr><tr><td>J143243.93+330746.7</td><td>1841</td><td>53491</td><td>513</td><td>2.09</td><td>2.058</td><td>0.28 ± 0.05</td><td>0.44 ± 0.05</td><td>17</td><td>...</td><td>...</td></tr><tr><td>J143657.87+291100.6</td><td>2138</td><td>53757</td><td>259</td><td>1.77</td><td>1.596</td><td>0.35 ± 0.11</td><td>0.56 ± 0.11</td><td>8</td><td>...</td><td>...</td></tr><tr><td>J143912.05+111740.6</td><td>1711</td><td>53535</td><td>374</td><td>2.58</td><td>2.418</td><td>0.50 ± 0.08</td><td>0.68 ± 0.09</td><td>10</td><td>20.27 ± 0.24<sup>e</sup></td><td>20.10 ± 0.10<sup>h</sup></td></tr><tr><td>J144929.27+333811.0</td><td>1646</td><td>53498</td><td>065</td><td>2.17</td><td>2.021</td><td>0.47 ± 0.12</td><td>0.62 ± 0.12</td><td>7</td><td>...</td><td>...</td></tr><tr><td>J145432.54+343523.9</td><td>1384</td><td>53121</td><td>269</td><td>1.61</td><td>1.580</td><td>0.59 ± 0.16</td><td>0.59 ± 0.13</td><td>6</td><td>...</td><td>...</td></tr><tr><td>J145953.25+012944.2</td><td>0538</td><td>52029</td><td>039</td><td>1.66</td><td>1.623</td><td>0.36 ± 0.17</td><td>0.62 ± 0.16</td><td>5</td><td>...</td><td>...</td></tr><tr><td>J150738.73+415530.6</td><td>1291</td><td>52735</td><td>505</td><td>1.79</td><td>1.674</td><td>0.50 ± 0.11</td><td>0.85 ± 0.16</td><td>7</td><td>...</td><td>...</td></tr><tr><td>J152209.12+083020.0</td><td>1721</td><td>53857</td><td>192</td><td>1.93</td><td>1.627</td><td>0.32 ± 0.10</td><td>0.71 ± 0.11</td><td>8</td><td>...</td><td>...</td></tr><tr><td>J160320.76+170117.8</td><td>2200</td><td>53875</td><td>318</td><td>1.99</td><td>1.890</td><td>0.16 ± 0.07</td><td>0.43 ± 0.08</td><td>12</td><td>...</td><td>19.20 ± 0.15</td></tr><tr><td>J160457.51+220300.5</td><td>2205</td><td>53793</td><td>328</td><td>1.98</td><td>1.641</td><td>1.34 ± 0.12</td><td>1.60 ± 0.13</td><td>7</td><td>...</td><td>...</td></tr><tr><td>J161526.65+264813.8</td><td>1576</td><td>53496</td><td>394</td><td>2.18</td><td>2.118</td><td>0.28 ± 0.06</td><td>0.56 ± 0.12</td><td>15</td><td>...</td><td>20.55 ± 0.15</td></tr><tr><td>J162321.43+135532.4</td><td>2202</td><td>53566</td><td>132</td><td>1.75</td><td>1.751</td><td>0.39 ± 0.09</td><td>0.60 ± 0.10</td><td>9</td><td>...</td><td>19.95 ± 0.15</td></tr><tr><td>J164610.20+232923.0</td><td>1424</td><td>52912</td><td>301</td><td>2.06</td><td>1.998</td><td>0.36 ± 0.05</td><td>0.58 ± 0.08</td><td>15</td><td>...</td><td>19.68 ± 0.15</td></tr><tr><td>J170542.92+354340.4</td><td>0974</td><td>52427</td><td>561</td><td>2.01</td><td>2.038</td><td>0.70 ± 0.08</td><td>1.40 ± 0.10</td><td>10</td><td>...</td><td>...</td></tr><tr><td>J212329.47-005053.0</td><td>0987</td><td>52523</td><td>103</td><td>2.26</td><td>2.060</td><td>0.11 ± 0.02</td><td>0.11 ± 0.02</td><td>41</td><td>...</td><td>19.25 ± 0.20<sup>i</sup></td></tr></tbody></table>

Notes. (a) Rest-frame equivalent width of the absorption detected at the position of the C I lines (measurements drawn from SDSS spectra).

(b) Average S/N per pixel around the C I lines. (c) Neutral atomic-hydrogen column density measured from Lyα absorption in SDSS spectra.

(d) N(HI) value derived from follow-up UVES spectroscopy. (e) Noterdaeme et al. (2009b). (f) Wang et al. (2012). (g) Noterdaeme et al. (2010a).

(h) Noterdaeme et al. (2008b). (i) Milutinovic et al. (2010). (j) Ledoux et al. (2006).
<table><caption>**Table 1.** continued.</caption><thead><tr><th>QSO name</th><th>Plate</th><th>MJD</th><th>Fibre</th><th>z<sub>em</sub></th><th>z<sub>abs</sub></th><th>W<sub>r</sub>(λ1560)<sup>a</sup><br>[Å]</th><th>W<sub>r</sub>(λ1656)<sup>a</sup><br>[Å]</th><th>S/N<sup>b</sup></th><th>log N(H I)<sub>SDSS</sub><sup>c</sup><br>[atoms cm<sup>-2</sup>]</th><th>log N(H I)<sub>UVES</sub><sup>d</sup><br>[atoms cm<sup>-2</sup>]</th></tr></thead><tbody><tr><td>J222910.16+141402.2</td><td>0738</td><td>52 521</td><td>382</td><td>2.11</td><td>1.586</td><td>0.33 ± 0.13</td><td>0.47 ± 0.13</td><td>6</td><td>...</td><td>...</td></tr><tr><td>J225719.04–100104.7</td><td>0724</td><td>52 254</td><td>150</td><td>2.08</td><td>1.836</td><td>0.64 ± 0.05</td><td>0.89 ± 0.05</td><td>17</td><td>...</td><td>20.38 ± 0.07</td></tr><tr><td>J233133.05–090246.6</td><td>0646</td><td>52 523</td><td>619</td><td>2.44</td><td>1.734</td><td>0.41 ± 0.04</td><td>0.35 ± 0.04</td><td>20</td><td>...</td><td>...</td></tr><tr><td>J233156.49–090802.0</td><td>0646</td><td>52 523</td><td>616</td><td>2.66</td><td>2.143</td><td>1.34 ± 0.13</td><td>1.17 ± 0.11</td><td>8</td><td>...</td><td>21.15 ± 0.15</td></tr><tr><td>J233633.81–105841.5</td><td>0647</td><td>52 553</td><td>201</td><td>2.04</td><td>1.829</td><td>0.22 ± 0.05</td><td>0.40 ± 0.07</td><td>15</td><td>...</td><td>...</td></tr><tr><td>J234023.67–005327.1</td><td>0385</td><td>51 877</td><td>204</td><td>2.09</td><td>2.054</td><td>0.11 ± 0.02</td><td>0.17 ± 0.02</td><td>37</td><td>...</td><td>20.33 ± 0.03</td></tr><tr><td>J235057.87–005210.0</td><td>0386</td><td>51 788</td><td>137</td><td>3.03</td><td>2.426</td><td>0.62 ± 0.08</td><td>0.63 ± 0.08</td><td>11</td><td>20.39 ± 0.27<sup>e</sup></td><td>20.50 ± 0.10<sup>j</sup></td></tr></tbody></table>

are in the expected range. This gives confidence in the derived values and their associated uncertainties. In the case of the outlier seen in the lower part of the plot (i.e., at $z_{abs} = 1.526$ towards SDSS J125552.60+223424.4), an unrelated blend to the λ1656 line is a probable reason for the observed deviation.

All the systems are located within about $2\sigma$ of the boundaries defined by the optically thin regime, on one hand, and the relation expected for heavily saturated profiles, on the other. We note that because of their large equivalent widths ($W_r \ge 0.4$ Å), most of the absorbers, especially those with equivalent-width ratios consistent with the optically thin regime, are probably made of numerous velocity components.

Here, we assumed that the C I ground state is solely responsible for the absorption lines while in reality the absorption from the two fine-structure energy levels of the neutral-carbon ground state ($^3P_1$ and $^3P_2$) could contribute mildly to the measured equivalent widths. However, this will affect the equivalent widths of both of the λ1560 and λ1656 lines in the same way so that any departure from the assumed relations due to this blending will be small.

## 3.2. Completeness

Before discussing the number of C I absorbers per unit redshift ($n_{C1}$; see following section), we first need to estimate the completeness of the sample. Given the resolving power $R$ of the SDSS spectra, the C I λ1560 line rest-frame equivalent width limit is given by

$$W_{r,\text{lim}}(\lambda1560) \approx n \times \frac{1560}{R} \times S/N^{-1} \quad (1)$$

where $n = 2$ is the number of standard deviations above which the peak absorption must be detected and $S/N > 4$ is the limit on the S/N per pixel at the corresponding line position. The FWHM of the lines is sampled by two velocity pixels with constant values. Our survey should therefore be complete down to $W_{r,\text{lim}}(\lambda1560) \approx 0.4$ Å.

We checked the exact level of completeness of our survey at this equivalent-width limit by implementing the following procedure. For this purpose, we used the same data set that we used for calculating $n_{C1}$ in Sect. 3.3, i.e., the same quasar sample, the same [$z_{min}, z_{max}$] values, and the same mean $S/N > 4$ limit. We then randomly selected 1000 QSO spectra and introduced an artificial C I system of rest-frame equivalent width $W_r(\lambda1560)$ at different positions in the spectra where the local S/N at both C I lines is greater than four. The distribution of the equivalent width ratio of the λ1560 and λ1656 lines is assumed to be a normal distribution with a dispersion corresponding to what is seen in Fig. 2. However, neither the equivalent width ratio nor the exact number of artificial systems used in the simulation has any significant impact on the completeness we infer.

<figure><img src="image_4.png" /><figcaption>**Fig. 3.** Completeness of our C I search as a function of C I λ1560 rest-frame equivalent width. Each point is based on a total of 40 000 artificial C I systems introduced in 1000 randomly selected SDSS QSO spectra. The vertical line shows our chosen completeness limit of $W_r(\lambda1560) = 0.4$ Å.</figcaption></figure>

We implemented about 40 000 C I systems that we sought to recover by using the same automatic procedure as described in Sect. 2. We varied $W_r$ over the range 0.1–1.0 Å and defined the completeness as the ratio of the number of recovered systems to the total number of systems introduced in the spectra. The results are displayed in Fig. 3. It can be seen that the completeness is greater than 80% for $W_r(\lambda1560) \ge 0.4$ Å.

## 3.3. Number of absorbers per unit redshift

We calculated the sensitivity function, $g(z)$, of our survey, i.e., the number of lines of sight that probe a given redshift $z$ and that have $S/N > 4$ at the expected positions of both C I lines. This function is shown in Fig. 4. It combines the [$z_{min}, z_{max}$] pairs, previously defined in Sect. 2, for all the lines of sight. We further excluded the regions with velocities relative to the QSO emission redshifts lower than 5000 km s$^{-1}$, which could, in principle, be influenced by the quasar (see Sect. 3.4). The uncertainties on $z_{em}$ are of the order of 500 km s$^{-1}$ (see, e.g., Pâris et al. 2012) and are therefore small enough not to affect the statistics. The total statistical absorption path length probed by the QSO sample over $z = 1.50–3.35$ is $\Delta z \approx 13 000$ with an average redshift
<figure><img src="image_5.png" /><figcaption>**Fig. 4.** Statistical sensitivity function of our C I survey, *g(z)* (solid curve; left-hand axis). The dips seen at *z* ~ 2.35 (resp. 2.6) are caused by C I λ1656 (resp. C I λ1560) falling over bad pixels at the junction between the two SDSS CCD chips. Dips at higher redshifts (*z* ~ 2.8, 3.05) are due to telluric lines (e.g., [O I] λ6300). The histogram (right-hand axis) shows the C I-absorber number counts in different redshift bins for intervening systems with rest-frame equivalent width above our completeness limit, i.e., *W*<sub>r</sub>(λ1560) ≥ 0.4 Å.</figcaption></figure>

<table><caption>**Table 2.** Intervening C I-absorber number counts.</caption><thead><tr><th>C I samples</th><th>W<sub>r</sub>(λ1560)<br>[Å]</th><th>z<sub>min</sub></th><th>z<sub>max</sub></th><th>Δz</th><th># of C I</th><th>n<sub>C I</sub><sup>a</sup><br>(×10<sup>-3</sup>)</th></tr></thead><tbody><tr><td rowspan="2">All systems</td><td rowspan="2">≥0.4</td><td>1.5</td><td>1.9</td><td>6268</td><td>22</td><td>3.5 ± 0.8</td></tr><tr><td>1.9</td><td>3.35</td><td>6710</td><td>8</td><td>1.2 ± 0.4</td></tr><tr><td rowspan="2">Strongest</td><td rowspan="2">≥0.64</td><td>1.5</td><td>1.9</td><td>6268</td><td>10</td><td>1.6 ± 0.5</td></tr><tr><td>1.9</td><td>3.35</td><td>6710</td><td>4</td><td>0.6 ± 0.3</td></tr><tr><td rowspan="2">Weaker</td><td rowspan="2">≥0.4, &lt;0.64</td><td>1.5</td><td>1.9</td><td>6268</td><td>12</td><td>1.9 ± 0.6</td></tr><tr><td>1.9</td><td>3.35</td><td>6710</td><td>4</td><td>0.6 ± 0.3</td></tr></tbody></table>

**Notes.** (a) Prior to correcting for incompleteness (i.e., ~18% at $W_r(\lambda1560) = 0.4$ Å).

$\langle z \rangle = 1.9$. From Fig. 4, it is apparent that the sensitivity of the survey is an increasing function towards lower redshifts. One therefore expects a larger number of absorbers to be found at $z < 2$. This is what is observed in practice as indicated by the redshift histogram of the detected intervening C I absorbers overplotted in the same figure. On the other hand, the small number of C I systems found at $z_{abs} > 2.2$ (i.e., eight systems out of a total of 66 systems or, equivalently, 12% of the sample) is striking.

To investigate this further, we calculated the number of intervening C I absorbers per unit redshift, $n_{C I}$, in two redshift bins of roughly equal total absorption path length (with boundary redshift $z = 1.9$). Here, we only considered the systems with rest-frame equivalent widths above the completeness limit of the survey, i.e., $W_r(\lambda1560) \ge 0.4$ Å. The results are summarized in Table 2 and shown in Fig. 5. In the table, we also separated the strongest from the weaker absorbers (around a median rest-frame equivalent width of 0.64 Å) but there is no obvious difference between the redshift evolution of these two groups.

<figure><img src="image_6.png" /><figcaption>**Fig. 5.** Evolution of the number of intervening C I absorbers per unit redshift, $n_{C I}$, for systems with rest-frame equivalent widths $W_r(\lambda1560) \ge 0.4$ Å. The measurements (solid) are corrected for incompleteness at the latter limit. For clarity purposes, the uncorrected data (dashed) are displayed with a slight offset in abscissa. The curve shows the expected redshift behaviour of a non-evolving population.</figcaption></figure>

The $n_{C I} \sim 1.4 \times 10^{-3}$ we measure in the higher redshift bin ($1.9 < z_{abs} < 3.35$), taking the effect of incompleteness estimated in Sect. 3.2 into account, implies that C I systems with $W_r(\lambda1560) \ge 0.4$ Å are more than one hundred times rarer than DLAs at $z_{abs} = 2.5$ (see Noterdaeme et al. 2012). An evolution of $n_{C I}$, with nearly thrice as many systems below $z = 1.9$ than above that is also observed. Compared to the redshift behaviour of a non-evolving population, this is significant at the $4.3\sigma$ level. Such an evolution is interesting and should be studied further since it depends on the balance between dust shielding and the UV radiation field. This may imply a strong evolution of the shielding of 10 eV photons by dust between $z = 2.5$ and $z = 1.5$.

### 3.4. Proximate systems

There are 14 C I systems with velocities relative to the QSO emission redshifts that are lower than 5000 km s<sup>-1</sup>. These could be associated with the QSO host galaxy or nearby environment. Six systems even have absorption redshifts greater than the corresponding QSO emission redshifts (by up to ~4000 km s<sup>-1</sup>), which is difficult to explain by high peculiar velocities in intervening systems. Imposing the same data-quality cuts and minimum equivalent widths as in the previous section, we find that the incidence of C I absorbers at small velocity differences from the quasars is consistent with that of intervening systems. However, associated errors are large due to small number statistics.

Because of the clustering of galaxies around the massive QSO host galaxies, an excess of proximate C I systems could be expected. However, C I is a fragile species that can easily be photo-ionized by the intense UV radiation emitted by the QSO engine. Interestingly, the lack of a significant excess of proximate systems was also observed by Prochaska et al. (2008)
when considering DLAs. In this case, the abundance of prox-
imate DLAs is only a factor of two larger than that of the
overall DLA population. A similarly low over-abundance fac-
tor was observed by Finley et al. (2013) for strong DLAs with
log N(HI) > 21.3 (atoms cm⁻²). This is much less than expected
based on clustering arguments alone.

In this work, we do not observe that the properties of proximate C I systems are different from those of intervening C I systems. This is true for redshift, C I equivalent-width, H I content, reddening, and UV bump-strength distributions. Nevertheless, because of the possibly different origins of these absorbers and the possible requirement of strong dust shielding from the nearby QSOs, we distinguish in the following proximate C I systems from the rest of the population and comment, whenever possible, on proximate C I systems of interest.

4. Evidence of dust

4.1. QSO optical colours

To assess the impact of the C I absorbers on their background QSOs and check for the existence of dust in these systems, we first consider the observed colours of these QSOs and compare them with the colours of the overall QSO population used as a control sample.

In Fig. 6, we show the distributions of (g − r), (r − i), and (r − z) colours for the 41 696 QSOs whose spectra were searched for C I absorption. In the upper panels of this figure, it is ap-
parent that the lines of sight with detected C I absorption do not
distribute in the same way as the other lines of sight. They are
all displaced together towards redder optical colours compared
to the average loci of the QSO redshift sequences. The effect
is most easily seen in the lower panels of Fig. 6, which com-
pare the colour histograms of the two QSO populations (i.e., the
C I-detected QSO sample and the overall QSO sample). The two-
sided Kolmogorov-Smirnov test probability that the two samples
are drawn from the same parent distribution is as low as ≤10⁻¹⁰.
The typical colour excess is ~0.15 mag, i.e., about five times
more than the mean (r − z) colour excess of 0.03 mag derived
by Vladilo et al. (2008) in zₐbs ≈ 2.8 DLAs from SDSS DR 5. A
similar result for DLAs was found by Khare et al. (2012) based
on (g − i) colours of SDSS-DR 7 QSOs. This is clear evidence
for the presence of dust among C I absorbers.

4.2. Reddening

Motivated by the unequivocal signature of dust in the form of a
colour excess of the background QSOs with detected C I absorp-
tion, we now aim at constraining the properties and the nature of
dust in these systems.

For each of the 66 QSOs with foreground C I absorbers,
we derived the QSO reddening, E(B − V), following the same
approach as used in Srianand et al. (2008a) and Noterdaeme
et al. (2009a, 2010a), for example. First, we corrected the
QSO spectra for Galactic reddening using the extinction maps
from Schlegel et al. (1998). We then fitted the spectra with the
SDSS QSO composite spectrum from Vanden Berk et al. (2001)
shifted to that QSO emission redshift and reddened with ei-
ther a Small Magellanic Cloud (SMC), Large Magellanic Cloud
(LMC), LMC2 super shell, or Milky Way (MW) extinction law
(Gordon et al. 2003) at the C I-absorber redshift. Our procedure
is illustrated in the left-hand panel of Fig. 7. The fit with the low-
est χ² value indicates the most representative extinction law for

a given absorber. The latter is specified in column “Best fit” of
Table 3, and the corresponding E(B − V) value is given in the
preceding column.

For each QSO line of sight exhibiting C I absorption, we de-
fined a control sample (hereafter denoted as “C.S.”) made of
SDSS-DR 7 QSOs from the searched sample having an emission
redshift within ±0.05 and a z-band magnitude within ±0.1 mag
from those of the QSO under consideration. In some instances,
this resulted in a sample of less than 30 QSOs in which case
we increased the above maximum magnitude difference by steps
of 0.01 until the number of QSOs in the control sample reached
(or exceeded) 30. We then applied the exact same fitting proce-
dure as described in the previous paragraph to each QSO spectra
from the control sample. Table 3 lists the number of QSOs and
the median reddening and standard deviation of the distribution
of E(B − V) values in each control sample (see the upper right
panel of Fig. 7 for an illustration). The values given in Table 3
correspond to the most representative extinction law previously
determined for that particular C I-detected QSO spectrum.

In the left-hand panel of Fig. 8, we show the histogram of reddening for the sample of 66 QSO lines of sight with detected C I absorbers compared to the cumulative control sample (calculated as the sum of the normalized distributions of individual control samples). As in Sect. 4.1, an offset between the two samples is apparent. The mean reddening induced by C I systems is 0.065 mag. A tail in the histogram of the C I-detected lines of sight is observed with E(B − V) values up to ~0.3 mag.

4.3. The 2175 Å extinction feature

A number of C I-detected QSO spectra are best matched by an extinction law exhibiting the absorption feature at rest-frame wavelength 2175 Å. To measure the strength of this UV bump (denoted A<sub>bump</sub>), we followed a prescription similar to the one used by Jiang et al. (2010) where the observed QSO spectrum was fitted with the SDSS QSO composite spectrum reddened via a parametrized pseudo-extinction law made of a smooth component and a Drude component. However, here we fixed the wavelength and width of the bump to the Galactic values determined by Fitzpatrick & Massa (2007). Both of these quantities indeed show little variation from line of sight to line of sight through the Galaxy and the Magellanic Clouds. This then limits the number of free parameters and prevents the fit from diverging towards very wide and shallow solutions that could be non-physical. Indeed, imperfect matching of the observed QSO continuum by the smooth component is expected owing to intrinsic QSO-shape variations (e.g., Pitman et al. 2000).

The fitting process is illustrated in the left-hand panel of
Fig. 7. The shaded area represents the measure of the bump
strength. This is the difference between the above best-fit func-
tion and the same function but considering only its smooth com-
ponent (i.e., with the Drude component set to zero). The A<sub>bump</sub>
values are listed for each absorber in Table 3. As previously
done for determining reddening (see Sect. 4.2), we also defined
a QSO control sample whose measured A<sub>bump</sub> distribution is
shown in the lower right-hand panel of Fig. 7 (i.e., for the given
QSO emission redshift and z-band magnitude). Table 3 gives the
median and standard deviation of this distribution for each con-
trol sample.

The histogram of bump strengths in the C I-absorber sample
is displayed in the right-hand panel of Fig. 8. One can see from
this figure that more than a quarter of the C I systems feature ab-
sorption at 2175 Å. This strengthens the result from the previous
section that significant reddening of the background QSOs by
<figure><img src="image_7.png" /><figcaption>**Fig. 6.** *Upper panels:* colours of the 41 696 SDSS QSOs whose spectra were searched for C I absorption versus QSO emission redshift. *(g − r)*, *(r − i)* and *(r − z)* colours are displayed in the *left-*, *middle-*, and *right-hand panels*, respectively. The thick blue line shows the median QSO colour as a function of redshift. Thin blue lines show one standard deviation away from that median. The red circles indicate the 66 lines of sight passing through C I-detected gas. *Lower panels:* distribution of colour excesses (defined as, e.g.: *(g − r) − ⟨g − r⟩*<sub>QSO</sub> for Δ*(g − r)*)) for the 66 QSOs whose spectra were found to exhibit C I absorption (thick red histogram) compared to the whole QSO sample (thin-hashed histogram, scaled down by a factor of 632 to have the same area). In each panel (from left to right: Δ*(g − r)*, Δ*(r − i)* and Δ*(r − z)*), vertical dotted red lines indicate the median colour excesses of the sample of QSOs with detected C I absorption. The median values for the whole QSO sample are shown by vertical dashed lines.</figcaption></figure>

dust is induced by some of the C I absorbers. We come back to
this and quantify the effect in Sect. 6.

In the following, we use the median *E*(B − *V*) and *A*<sub>bump</sub> val-
ues of the control samples, i.e., ⟨*E*(B − *V*)⟩<sub>C.S.</sub> and ⟨*A*<sub>bump</sub>⟩<sub>C.S.</sub>, to
define the exact colour excess and bump strength towards a given

C I-detected QSO line of sight: E(B − V) = E(B − V)measured −
⟨E(B − V)⟩C.S., and likewise for A<sub>bump</sub>. These zero-point correc-
tions are usually almost negligible (see Table 3). In addition, the
standard deviations of E(B − V) and A<sub>bump</sub> values in each control
sample provide an estimate of the uncertainty due to intrinsic
<figure><img src="image_8.png" /></figure>

**Fig. 7.** Left-hand panel: illustration of reddening measurement based on QSO SED fitting. The blue curve shows our best fit to the QSO continuum using the QSO template spectrum from Vanden Berk et al. (2001) and the best-fitting LMC extinction law at the redshift of the C I absorber. The non-extinguished QSO template is displayed in grey. The additional presence of a 2175 Å absorption feature is shown by the light-red shaded area (i.e., A<sub>bump</sub>). The broken orange line at zero ordinate indicates the wavelength regions devoid of QSO emission lines, strong absorption lines and sky residuals, that were used during the χ<sup>2</sup> minimization. Upper right-hand panel: distribution of reddening measured in the QSO control sample assuming either a SMC (red) or the best-fitting extinction law (in this case, LMC; blue). The reddening values of the line of sight under consideration are indicated by downward arrows. Lower right-hand panel: same as above but for the distribution of 2175 Å bump strengths.

<figure><img src="image_9.png" /></figure>

**Fig. 8.** Left-hand panel: histogram of reddening for the sample of 66 QSOs whose spectra were found to exhibit C I absorption (multi-colour). The colour/hash coding for the different C I-detected lines of sight relates to the best-fitting extinction law (red: SMC; blue: LMC; green: LMC2; purple: MW). The distribution of reddening from the QSO control sample (C.S.), i.e., the sum of the normalized distributions of individual QSO control samples, is displayed with dashed lines. Right-hand panel: same as in the left-hand panel but for the histogram of 2175 Å bump strengths. The distribution from the QSO control sample is well represented by a Gaussian profile (solid curve).
<figure><img src="image_10.png" /></figure>

**Fig. 9.** Neutral atomic-hydrogen column-density distribution of C I absorbers (red-filled histogram; left-hand axis). Proximate systems are displayed in salmon. The N(HI) distribution of intervening DLAs from SDSS DR7 (Noterdaeme et al. 2009b, blue-hashed histogram) is over-plotted using a different scaling (right-hand axis) so that the areas of the two histograms above log N(HI) = 20.3 are equal. The curve is a fit of the distribution in the sub-DLA regime from Prochaska et al. (2014). The inset shows the cumulative N(HI) distributions of the intervening C I and DLA samples starting at log N(HI) = 20.3 (atoms cm⁻²) and the two-sided Kolmogorov-Smirnov test probability that the two distributions come from the same parent population.

QSO-shape variations (see Pitman et al. 2000), hence the sig-
nificance of the reddening induced by each C I absorber and the
significance of associated 2175 Å absorption, respectively.

5. HI content

As part of a spectroscopic campaign, which we will describe in
a companion paper, we followed up the C I absorbers from the
overall sample that are observable from the southern hemisphere
using VLT/UVES. We present in Fig. 9 the H I column-density
distribution of this C I absorber subsample (referred to in the fol-
lowing as the H I subsample) and compare it with the distribution
of N(HI) from systematic DLA and/or sub-DLA surveys.

We secured H I column-density measurements for most
of the systems in the overall sample that have a declina-
tion of δ < +28 deg: i.e., 14 out of 16 systems at redshifts
z<sub>abs</sub> > 1.8 (the two exceptions being the lines of sight towards
SDSS J091721.37+015448.1 and J233633.81–105841.5) and
four out of eight systems at z<sub>abs</sub> ≈ 1.75 (see Table 1). While
H I column densities derived from UVES spectroscopy are usu-
ally more accurate, the last two columns of Table 1 show that
they confirm those derived directly from SDSS spectra as seen
from the five systems at z<sub>abs</sub> > 2.2 where this measurement
could be done from both datasets. For this reason, we comple-
ment our UVES measurements with the values we derived us-
ing SDSS spectra for the three systems at z<sub>abs</sub> > 2.2, for which
high-resolution spectroscopic data are not available because the
background QSOs are too far north for the VLT to observe them.
The H I subsample thus comprises a total of 21 systems.

<figure><img src="image_11.png" /></figure>

**Fig. 10.** C I λ1560 rest-frame equivalent width versus N(HI), the neutral atomic-hydrogen column density of C I absorbers. Intervening (resp. proximate) systems are displayed as filled (resp. empty) black circles. Data points are circled in red as they correspond to the subsample of C I systems with measured N(HI), i.e., the H I subsample. The dashed (resp. dashed-dotted) curve materializes the limit above which systems are expected to have a metallicity larger than one-tenth of solar (resp. in excess of solar; see text). The completeness limit of the survey is indicated by the horizontal line.

In Fig. 9, we compare the observed H I column-density distribution of C I-selected absorbers with that of H I-selected DLAs (from SDSS DR7 as well; Noterdaeme et al. 2009b). In this figure, we also show the expected number of sub-DLAs using the fitted distribution function from Prochaska et al. (2014). We find that a large number of C I absorbers have neutral atomic-hydrogen column densities slightly below the conventional DLA limit (N(HI) = 2 × 10²⁰ atoms cm⁻²) and are therefore classified as strong sub-DLAs. However, the number of C I absorbers among sub-DLAs is much less than among DLAs, indicating that efficient shielding is much more difficult to obtain below the DLA limit. Though rare, the existence of C I absorbers with low neutral atomic-hydrogen column densities supports the presence of dust in these systems. The dust-to-gas ratio in these systems has to be high enough so that the absorption of UV photons by dust allows C I to be present in large amounts.

No C I system is found with log N(HI) < 19 (atoms cm⁻²). This is a regime where shielding of UV photons becomes extremely difficult even in the presence of dust. However, it is possible that such systems are missed in our search. Indeed, as seen in Fig. 10, there is a tendency for neutral atomic-hydrogen column density to increase with C I equivalent width. The gradually decreasing completeness fraction of the survey below W<sub>r</sub>(λ1560) ≈ 0.4 Å (see Sect. 2) would therefore preclude low N(HI) systems from appearing in our sample.

Above the DLA limit, where the incompleteness of our sur-
vey is less of a problem, it appears that C I-selected absorbers do
not follow the statistics of H I-selected DLAs. Although there
are only nine C I systems with log N(HI) > 20.3 (atoms cm⁻²),
it is apparent that the overall N(HI) distribution of C I systems
is relatively flat. A two-sided Kolmogorov-Smirnov test applied
<table><caption>**Table 3.** SED fitting of QSO spectra with detected C1 absorbers.</caption><thead><tr><th rowspan="2">QSO name</th><th rowspan="2">z<sub>em</sub></th><th rowspan="2">z<sub>abs</sub></th><th rowspan="2">E(B-V) [mag]</th><th rowspan="2">Best fit</th><th rowspan="2">A<sub>bump</sub></th><th colspan="2">C.S.<sup>a</sup>: E(B-V) [mag]</th><th colspan="2">A<sub>bump</sub></th></tr><tr><th>N<sub>QSO</sub></th><th>Median</th><th>σ</th><th>Median</th><th>σ</th></tr></thead><tbody><tr><td>J021606.13-002104.9</td><td>2.22</td><td>1.737</td><td>0.039</td><td>SMC</td><td>0.097</td><td>173</td><td>-0.008</td><td>0.053</td><td>-0.009</td><td>0.123</td></tr><tr><td>J030049.24-072137.8</td><td>2.11</td><td>1.536</td><td>0.075</td><td>SMC</td><td>0.160</td><td>448</td><td>-0.004</td><td>0.048</td><td>0.002</td><td>0.136</td></tr><tr><td>J080801.74+330009.2</td><td>1.90</td><td>1.888</td><td>0.029</td><td>SMC</td><td>-0.032</td><td>151</td><td>0.001</td><td>0.044</td><td>-0.028</td><td>0.112</td></tr><tr><td>J081116.06+083837.7</td><td>2.07</td><td>1.906</td><td>0.081</td><td>SMC</td><td>-0.027</td><td>247</td><td>-0.003</td><td>0.049</td><td>-0.031</td><td>0.127</td></tr><tr><td>J081540.60+264021.6</td><td>1.94</td><td>1.681</td><td>0.109</td><td>LMC</td><td>0.364</td><td>428</td><td>-0.012</td><td>0.043</td><td>-0.027</td><td>0.142</td></tr><tr><td>J082003.40+155932.9</td><td>1.95</td><td>1.547</td><td>0.203</td><td>LMC2</td><td>0.420</td><td>373</td><td>-0.004</td><td>0.069</td><td>-0.012</td><td>0.142</td></tr><tr><td>J085206.65+193548.4</td><td>1.91</td><td>1.788</td><td>0.050</td><td>SMC</td><td>0.032</td><td>580</td><td>-0.003</td><td>0.044</td><td>-0.020</td><td>0.112</td></tr><tr><td>J085437.60+031734.8</td><td>2.24</td><td>1.567</td><td>0.070</td><td>SMC</td><td>0.263</td><td>192</td><td>-0.003</td><td>0.054</td><td>0.050</td><td>0.106</td></tr><tr><td>J085726.79+185524.3</td><td>1.89</td><td>1.730</td><td>0.016</td><td>SMC</td><td>0.063</td><td>81</td><td>-0.004</td><td>0.053</td><td>-0.064</td><td>0.107</td></tr><tr><td>J090558.75+553039.2</td><td>1.57</td><td>1.579</td><td>-0.019</td><td>SMC</td><td>0.070</td><td>692</td><td>-0.002</td><td>0.041</td><td>0.011</td><td>0.080</td></tr><tr><td>J090942.56+532126.5</td><td>2.07</td><td>1.628</td><td>0.003</td><td>SMC</td><td>0.003</td><td>527</td><td>-0.005</td><td>0.044</td><td>-0.019</td><td>0.140</td></tr><tr><td>J091516.27+071750.9</td><td>2.38</td><td>1.681</td><td>-0.010</td><td>SMC</td><td>0.117</td><td>40</td><td>-0.017</td><td>0.077</td><td>0.019</td><td>0.206</td></tr><tr><td>J091721.37+015448.1</td><td>2.18</td><td>2.107</td><td>0.121</td><td>SMC</td><td>0.227</td><td>86</td><td>-0.014</td><td>0.031</td><td>0.006</td><td>0.114</td></tr><tr><td>J092759.79+154321.8</td><td>1.81</td><td>1.731</td><td>0.217</td><td>LMC2</td><td>0.331</td><td>357</td><td>-0.008</td><td>0.059</td><td>-0.018</td><td>0.092</td></tr><tr><td>J094529.37+552525.7</td><td>2.24</td><td>1.867</td><td>0.120</td><td>SMC</td><td>0.015</td><td>314</td><td>-0.006</td><td>0.043</td><td>-0.016</td><td>0.130</td></tr><tr><td>J095114.97+545736.5</td><td>1.80</td><td>1.613</td><td>-0.004</td><td>SMC</td><td>-0.027</td><td>397</td><td>-0.004</td><td>0.042</td><td>-0.015</td><td>0.109</td></tr><tr><td>J101544.34+310617.2</td><td>1.56</td><td>1.596</td><td>0.016</td><td>SMC</td><td>-0.006</td><td>685</td><td>0.000</td><td>0.043</td><td>0.012</td><td>0.085</td></tr><tr><td>J104705.76+205734.5</td><td>2.01</td><td>1.775</td><td>0.159</td><td>LMC2</td><td>0.277</td><td>282</td><td>-0.021</td><td>0.055</td><td>-0.036</td><td>0.133</td></tr><tr><td>J105436.96+542343.7</td><td>1.50</td><td>1.513</td><td>0.014</td><td>SMC</td><td>-0.008</td><td>319</td><td>0.001</td><td>0.037</td><td>0.007</td><td>0.072</td></tr><tr><td>J105746.42+662250.5</td><td>2.34</td><td>2.088</td><td>0.040</td><td>SMC</td><td>0.112</td><td>277</td><td>-0.010</td><td>0.037</td><td>-0.046</td><td>0.120</td></tr><tr><td>J111756.53+143716.9</td><td>2.09</td><td>2.001</td><td>0.049</td><td>SMC</td><td>0.085</td><td>560</td><td>-0.005</td><td>0.036</td><td>-0.030</td><td>0.097</td></tr><tr><td>J112250.95+143732.5</td><td>2.05</td><td>1.554</td><td>-0.028</td><td>SMC</td><td>0.032</td><td>106</td><td>0.003</td><td>0.047</td><td>-0.026</td><td>0.138</td></tr><tr><td>J112913.29-023740.9</td><td>1.86</td><td>1.623</td><td>0.020</td><td>SMC</td><td>-0.151</td><td>603</td><td>-0.003</td><td>0.045</td><td>-0.026</td><td>0.124</td></tr><tr><td>J113341.29-005740.1</td><td>1.68</td><td>1.706</td><td>0.266</td><td>SMC</td><td>0.062</td><td>386</td><td>-0.002</td><td>0.042</td><td>0.002</td><td>0.094</td></tr><tr><td>J114141.84+444206.1</td><td>1.96</td><td>1.903</td><td>0.149</td><td>LMC2</td><td>0.368</td><td>398</td><td>-0.023</td><td>0.060</td><td>-0.028</td><td>0.109</td></tr><tr><td>J115153.87+150945.0</td><td>3.05</td><td>2.400</td><td>-0.007</td><td>SMC</td><td>0.110</td><td>92</td><td>-0.007</td><td>0.033</td><td>-0.105</td><td>0.176</td></tr><tr><td>J115705.51+615521.7</td><td>2.51</td><td>2.460</td><td>0.215</td><td>SMC</td><td>0.370</td><td>148</td><td>-0.014</td><td>0.033</td><td>-0.087</td><td>0.101</td></tr><tr><td>J120935.79+671715.7</td><td>2.03</td><td>1.843</td><td>0.094</td><td>MW</td><td>0.394</td><td>300</td><td>-0.013</td><td>0.028</td><td>-0.039</td><td>0.119</td></tr><tr><td>J122825.67+303038.6</td><td>3.91</td><td>3.098</td><td>-0.034</td><td>SMC</td><td>0.067</td><td>43</td><td>-0.008</td><td>0.061</td><td>0.412</td><td>0.471</td></tr><tr><td>J123714.61+064759.6</td><td>2.78</td><td>2.691</td><td>0.143</td><td>LMC2</td><td>0.431</td><td>35</td><td>-0.026</td><td>0.074</td><td>-0.168</td><td>0.156</td></tr><tr><td>J124708.43+500320.8</td><td>2.27</td><td>2.135</td><td>0.075</td><td>LMC</td><td>0.319</td><td>357</td><td>-0.012</td><td>0.029</td><td>-0.037</td><td>0.124</td></tr><tr><td>J124829.71+284858.1</td><td>1.54</td><td>1.513</td><td>-0.035</td><td>SMC</td><td>-0.049</td><td>1035</td><td>-0.002</td><td>0.041</td><td>0.022</td><td>0.092</td></tr><tr><td>J124841.45+302433.0</td><td>2.06</td><td>1.691</td><td>-0.025</td><td>SMC</td><td>-0.019</td><td>583</td><td>-0.004</td><td>0.042</td><td>-0.024</td><td>0.136</td></tr><tr><td>J125552.60+223424.4</td><td>1.82</td><td>1.526</td><td>0.036</td><td>SMC</td><td>0.141</td><td>898</td><td>-0.004</td><td>0.041</td><td>-0.015</td><td>0.134</td></tr><tr><td>J130225.28+211158.6</td><td>1.76</td><td>1.656</td><td>-0.004</td><td>SMC</td><td>-0.139</td><td>452</td><td>-0.012</td><td>0.036</td><td>-0.002</td><td>0.114</td></tr><tr><td>J130628.87+281550.8</td><td>2.10</td><td>2.012</td><td>0.078</td><td>SMC</td><td>0.015</td><td>465</td><td>-0.008</td><td>0.034</td><td>-0.030</td><td>0.099</td></tr><tr><td>J130828.43+584000.6</td><td>3.09</td><td>2.473</td><td>0.020</td><td>SMC</td><td>-0.158</td><td>54</td><td>-0.010</td><td>0.037</td><td>-0.144</td><td>0.211</td></tr><tr><td>J131129.11+222552.6</td><td>3.14</td><td>3.093</td><td>0.036</td><td>SMC</td><td>0.067</td><td>82</td><td>-0.011</td><td>0.036</td><td>-0.054</td><td>0.345</td></tr><tr><td>J131400.57+054319.5</td><td>1.89</td><td>1.583</td><td>0.046</td><td>SMC</td><td>0.270</td><td>340</td><td>-0.006</td><td>0.044</td><td>-0.013</td><td>0.120</td></tr><tr><td>J134122.51+185214.0</td><td>2.00</td><td>1.544</td><td>-0.016</td><td>SMC</td><td>0.000</td><td>32</td><td>0.003</td><td>0.053</td><td>-0.001</td><td>0.185</td></tr><tr><td>J134601.10+064408.4</td><td>2.09</td><td>1.512</td><td>0.015</td><td>SMC</td><td>0.034</td><td>118</td><td>-0.015</td><td>0.040</td><td>0.036</td><td>0.130</td></tr><tr><td>J135122.00+461828.5</td><td>1.81</td><td>1.606</td><td>0.142</td><td>LMC2</td><td>0.408</td><td>680</td><td>-0.005</td><td>0.067</td><td>-0.027</td><td>0.123</td></tr><tr><td>J141550.47+300146.9</td><td>2.08</td><td>1.676</td><td>0.074</td><td>SMC</td><td>0.226</td><td>83</td><td>-0.001</td><td>0.045</td><td>-0.033</td><td>0.245</td></tr><tr><td>J141606.79+180403.2</td><td>2.13</td><td>1.622</td><td>0.056</td><td>SMC</td><td>-0.027</td><td>108</td><td>-0.016</td><td>0.035</td><td>0.018</td><td>0.150</td></tr><tr><td>J143243.93+330746.7</td><td>2.09</td><td>2.058</td><td>0.010</td><td>SMC</td><td>-0.005</td><td>352</td><td>-0.003</td><td>0.041</td><td>-0.039</td><td>0.095</td></tr><tr><td>J143657.87+291100.6</td><td>1.77</td><td>1.596</td><td>0.147</td><td>LMC</td><td>0.514</td><td>812</td><td>-0.003</td><td>0.049</td><td>-0.014</td><td>0.111</td></tr><tr><td>J143912.05+111740.6</td><td>2.58</td><td>2.418</td><td>0.172</td><td>LMC2</td><td>0.271</td><td>52</td><td>-0.033</td><td>0.145</td><td>-0.080</td><td>0.098</td></tr><tr><td>J144929.27+333811.0</td><td>2.17</td><td>2.021</td><td>0.018</td><td>SMC</td><td>-0.157</td><td>89</td><td>-0.009</td><td>0.029</td><td>0.025</td><td>0.137</td></tr><tr><td>J145432.54+343523.9</td><td>1.61</td><td>1.580</td><td>0.293</td><td>LMC2</td><td>0.631</td><td>346</td><td>-0.002</td><td>0.058</td><td>-0.002</td><td>0.089</td></tr><tr><td>J145953.25+012944.2</td><td>1.66</td><td>1.623</td><td>0.157</td><td>LMC2</td><td>0.331</td><td>842</td><td>-0.006</td><td>0.057</td><td>0.011</td><td>0.089</td></tr><tr><td>J150738.73+415530.6</td><td>1.79</td><td>1.674</td><td>0.056</td><td>LMC</td><td>0.251</td><td>806</td><td>-0.005</td><td>0.040</td><td>-0.003</td><td>0.111</td></tr><tr><td>J152209.12+083020.0</td><td>1.93</td><td>1.627</td><td>0.017</td><td>SMC</td><td>-0.057</td><td>169</td><td>-0.016</td><td>0.035</td><td>-0.031</td><td>0.134</td></tr><tr><td>J160320.76+170117.8</td><td>1.99</td><td>1.890</td><td>-0.010</td><td>SMC</td><td>0.030</td><td>406</td><td>-0.004</td><td>0.043</td><td>-0.033</td><td>0.108</td></tr><tr><td>J160457.51+220300.5</td><td>1.98</td><td>1.641</td><td>0.263</td><td>LMC2</td><td>0.735</td><td>225</td><td>-0.007</td><td>0.075</td><td>-0.041</td><td>0.151</td></tr><tr><td>J161526.65+264813.8</td><td>2.18</td><td>2.118</td><td>0.035</td><td>SMC</td><td>0.120</td><td>324</td><td>-0.005</td><td>0.042</td><td>-0.037</td><td>0.095</td></tr><tr><td>J162321.43+135532.4</td><td>1.75</td><td>1.751</td><td>0.031</td><td>SMC</td><td>0.034</td><td>816</td><td>-0.005</td><td>0.036</td><td>0.006</td><td>0.093</td></tr><tr><td>J164610.20+232923.0</td><td>2.06</td><td>1.998</td><td>-0.027</td><td>SMC</td><td>-0.031</td><td>284</td><td>-0.004</td><td>0.043</td><td>-0.037</td><td>0.092</td></tr></tbody></table>

Notes. (a) C.S. = control sample.
<table><caption>**Table 3.** continued.</caption><thead><tr><th rowspan="2">QSO name</th><th rowspan="2">z<sub>em</sub></th><th rowspan="2">z<sub>abs</sub></th><th rowspan="2">E(B-V)<br>[mag]</th><th rowspan="2">Best fit</th><th rowspan="2">A<sub>bump</sub></th><th colspan="2">C.S.<sup>a</sup>:<br>N<sub>QSO</sub></th><th colspan="2">E(B-V)<br>[mag]</th><th colspan="2">A<sub>bump</sub></th></tr><tr><th>Median</th><th>σ</th><th>Median</th><th>σ</th></tr></thead><tbody><tr><td>J170542.92+354340.4</td><td>2.01</td><td>2.038</td><td>0.146</td><td>LMC2</td><td>0.364</td><td>210</td><td>-0.009</td><td>0.075</td><td>-0.030</td><td>0.077</td></tr><tr><td>J212329.47-005053.0</td><td>2.26</td><td>2.060</td><td>-0.011</td><td>SMC</td><td>0.000</td><td>32</td><td>0.002</td><td>0.041</td><td>-0.023</td><td>0.123</td></tr><tr><td>J222910.16+141402.2</td><td>2.11</td><td>1.586</td><td>-0.024</td><td>SMC</td><td>0.026</td><td>538</td><td>-0.004</td><td>0.041</td><td>0.006</td><td>0.136</td></tr><tr><td>J225719.04-100104.7</td><td>2.08</td><td>1.836</td><td>0.045</td><td>SMC</td><td>0.177</td><td>192</td><td>-0.007</td><td>0.039</td><td>-0.040</td><td>0.129</td></tr><tr><td>J233133.05-090246.6</td><td>2.44</td><td>1.734</td><td>0.101</td><td>LMC</td><td>0.393</td><td>38</td><td>0.002</td><td>0.089</td><td>0.084</td><td>0.215</td></tr><tr><td>J233156.49-090802.0</td><td>2.66</td><td>2.143</td><td>0.136</td><td>SMC</td><td>0.234</td><td>46</td><td>-0.008</td><td>0.078</td><td>-0.015</td><td>0.104</td></tr><tr><td>J233633.81-105841.5</td><td>2.04</td><td>1.829</td><td>0.004</td><td>SMC</td><td>-0.103</td><td>484</td><td>-0.003</td><td>0.045</td><td>-0.026</td><td>0.123</td></tr><tr><td>J234023.67-005327.1</td><td>2.09</td><td>2.054</td><td>0.057</td><td>SMC</td><td>-0.023</td><td>34</td><td>-0.001</td><td>0.047</td><td>-0.071</td><td>0.113</td></tr><tr><td>J235057.87-005210.0</td><td>3.03</td><td>2.426</td><td>0.003</td><td>SMC</td><td>0.001</td><td>42</td><td>-0.013</td><td>0.067</td><td>-0.098</td><td>0.169</td></tr></tbody></table>

to all absorbers with $log N(HI) > 20.3$ (atoms cm$^{-2}$) gives a probability of only 17% that the two distributions come from the same parent population (see inset of Fig. 9). This could be explained by a larger number of velocity components in higher HI column-density gas, thereby increasing the probability of detecting C I. Moreover, large amounts of shielded gas are probably the consequence of the line of sight passing through the absorbing galaxy at a low impact parameter, in which case we can expect the N(HI) distribution to be flatter than for the overall DLA population.

The strongest DLA found among the C I absorbers in the HI subsample has $N(HI) = 10^{21.8}$ atoms cm$^{-2}$. It is, however, located at $z_{abs} \approx z_{em}$. Even ignoring proximate systems, it is still surprising that three intervening DLAs with $N(HI) \ge 10^{21}$ atoms cm$^{-2}$ are present in such a small absorber sample. From DLA statistics alone, the probability of randomly selecting three DLAs that are that strong in a sample of six DLAs is only 6%. There is therefore a probable excess of strong DLAs amongst C I absorbers. While the dust content of these systems is significant (see Sects. 4.1 and 4.2), their dust-to-gas ratio must be limited. Indeed, dust reddening and extinction of the background QSOs will inevitably reduce the incidence of strong and dusty DLAs in magnitude-limited QSO samples. This implies that the actual proportion of strong DLAs amongst C I systems in general is likely to be even higher than what we found.

## 6. Empirical relations in the sample

Based on the results presented in the previous sections, we now study the empirical relations between the different quantities measured in this work: neutral atomic-carbon and neutral atomic-hydrogen contents, the reddening that C I-selected absorbers induce on their background QSOs, and the strength of possible 2175 Å extinction features. Because the C I λ1560 transition line is weaker, so it exhibits less saturation than C I λ1656, we adopted the equivalent width of the former as a proxy for the amount of neutral atomic carbon in the systems.

In Fig. 10, we plot the C I λ1560 rest-frame equivalent width versus $log N(HI)$ for the C I absorbers from the HI subsample. Both quantities appear to be weakly correlated. A Kendall rank-correlation test indicates that the significance of the correlation is only $1.8\sigma$. There is therefore a tendency for strong DLAs to have larger C I equivalent widths, but at the same time, for a given HI column density, the C I content can vary substantially from one system to another. High values of $W_r(\lambda1560)$ are observed in DLAs but also in sub-DLAs. The amount of shielded and probably cold gas could actually be large in some of these

sub-DLAs. From the optically thin approximation applied to the C I λ1560 absorption line and assuming the ionization equilibrium relation, $N(CI)/N(CII) \sim 0.01$, which is valid for the cold neutral medium (CNM; see, e.g., Jenkins & Tripp 2011), a lower limit on the gas metallicity can be derived:

$$[X/H] \ge 18.35 + \log \left( \frac{W_r(\lambda1560)}{0.01 \times N(HI)} \right). \quad (2)$$

For $W_r(\lambda1560) = 0.4$ Å and $log N(HI) = 20$ (atoms cm$^{-2}$), the metallicity should be close to solar. More generally, the curves in Fig. 10 were calculated using the above equation and assuming metallicities of one-tenth of solar and solar. Within measurement uncertainties, most of the C I systems lie in between these two curves. If the medium probed by the line of sight is a mixture of cold and warm gas, the metallicity of the systems will be even higher. However, if part of the hydrogen is in molecular form, the metallicity will be lower. For the whole C I-absorber sample, Eq. (2) implies a metallicity distribution ranging between $[X/H] = -1.4$ and metallicities in excess of solar with a median value of $[X/H] \approx -0.5$. This means that the metallicities of C I absorbers would on average be at least ten times greater than those of typical DLAs (for the latter, see, e.g., Rafelski et al. 2012). This should be confirmed by accurate measurements of metal column densities.

In Fig. 11, we display the relation between $E(B-V)$ and $log N(HI)$ among the C I absorbers from the HI subsample. Here again, the data points are highly scattered. Most of the systems are associated with low, albeit consistently non-zero, QSO reddening. Since most of the $N(HI)$ values are relatively low, the measured amounts of reddening, with median $E(B-V) \sim 0.045$, are actually remarkable. This departs from what is observed in the overall DLA population where the reddening is usually negligible (see, e.g., Vladilo et al. 2008; Khare et al. 2012). Vladilo et al. (2008) have shown that DLAs at $z_{abs} \approx 2.8$ typically induce a reddening $E(B-V) \sim 5 \times 10^{-3}$ mag. Apart from a few outliers, most of the C I systems have reddening properties that are consistent with those of the Galactic ISM. This is represented in Fig. 11. In the Galaxy, the reddening induced along a line of sight is indeed directly proportional to the neutral atomic-hydrogen column density with $E(B-V)/N(H) = 1.63 \times 10^{-22}$ mag atoms$^{-1}$ cm$^2$ (Gudennavar et al. 2012). Only two C I systems with large $N(HI)$ are more consistent with what is seen in typical DLAs and/or along SMC lines of sight, where the above ratio is lower than in the Galaxy. Two other C I systems with low $N(HI)$ might also deviate from the Galactic relation being consistent with a ten-times higher ratio. We caution, however, that the uncertainties on the reddening measurements
<figure><img src="image_12.png" /><figcaption>**Fig. 11.** Reddening of the background QSOs with detected C1 absorption versus N(HI) of the C1 absorbers. Symbol conventions are the same as in Fig. 10. The solid curve shows the relation observed in the local ISM where $E(B-V)/N(H) = 1.63 \times 10^{-22} \text{ mag atoms}^{-1} \text{ cm}^2$ (Gudennavar et al. 2012). The dotted curve illustrates a reddening per hydrogen atom ten times higher than that. The dashed curves correspond to the observations of the Magellanic Clouds (Gordon et al. 2003). The long-dashed curve shows the relation derived for typical DLAs at $z_{\text{abs}} \approx 2.8$ (e.g., Vladilo et al. 2008; Khare et al. 2012).</figcaption></figure>

are fairly large. If real, this would imply in the latter systems the existence of a grain chemistry more evolved than in the Galaxy with a higher ratio of big grains to very small grains (e.g., Pei 1992). This is opposite to the trend observed in the Magellanic Clouds. In such systems, strong 2175 Å absorption is expected. Interestingly, this is what is observed in practice because these two C1 systems exhibit two of the strongest three $A_{\text{bump}}$ values of the HI subsample.

To investigate the characteristics of the C1 absorbers further, we look in Fig. 12 in more detail at the properties of dust in these systems. In the left-hand panel of this figure, $E(B-V)$ and $W_r(\lambda_{1560})$ are found to be correlated with each other at the $4.4\sigma$ significance level. This is noteworthy because different degrees of saturation of the C1$\lambda_{1560}$ line are expected to produce scatter in this relation. This implies that the neutral-carbon content of the C1 systems is intimately related to the reddening induced along the line of sight or, equivalently, that the amounts of shielded gas and dust are tightly inter-connected. We also note that two of the highest three $E(B-V)$ values in this plot correspond to systems located at $z_{\text{abs}} \approx z_{\text{em}}$. This may, however, be the result of the small sample size because the reddening induced by the other proximate systems in our sample varies substantially from one system to another.

The relation between $E(B-V)$ and the UV bump strength, $A_{\text{bump}}$, which we previously determined independently from $E(B-V)$ (see Sect. 4.3), is shown in the right-hand panel of Fig. 12. It can be seen that both quantities are tightly correlated ($6.0\sigma$). A linear least-squares fit (linear correlation coefficient of $r = 0.77$), taking errors in both parameters into account, gives $E(B-V) \approx 0.43 \times A_{\text{bump}}$. The 2175 Å extinction feature is detected at more than $2\sigma$ (95% confidence level) in about 30%

of the C1 systems. In such cases, we find $A_{\text{bump}} \sim 0.4$ and $E(B-V) \sim 0.2$ mag or, equivalently, $A_V \sim 0.6$ mag. These values are comparable to what Budzynski & Hewett (2011) have found when targeting the strongest Mg II systems from SDSS DR6 [1 < $W_r(\lambda_{2796})$ < 5 Å] where the 2175 Å absorption is detected on a statistical basis alone (see also Jiang et al. 2011, for candidate 2175 Å absorption in similar systems). Interestingly, our measured UV bump strengths are also comparable to what has been observed along GRB lines of sight at similarly low levels of extinction, e.g., towards GRB 080605 (see Zafar et al. 2012).

In the local Universe, 2175 Å absorption can be observed with reddening values as low as ~0.2 mag (see, e.g., Fitzpatrick & Massa 2007). Even at such low levels of reddening, the UV bump is significantly stronger along Galactic lines of sight, i.e., by up to a factor of ten, than in the present C1 absorber sample and/or through GRB host galaxies. This discrepancy may be explained if these high-redshift systems probe regions of the ISM affected by a star formation that is more vigorous than in the Galaxy. A similar argument was proposed by Gordon et al. (2003) with the aim of explaining the variety of LMC and SMC extinction curves. In fact, most of the C1 absorbers at the high end of the reddening tail in our sample are best-fit using the extinction law of the LMC2 super-shell near the 30 Dor star-forming region<fnref n="2" /> (see left panel of Fig. 8). This means that the far-UV rise of the extinction curve is enhanced and the carriers of the 2175 Å absorption are depleted compared to Galactic lines of sight. This is probably the consequence of a high UV flux and/or the mechanical feedback from stars (e.g., Fischera & Dopita 2011) in the vicinity of the C1 systems. In contrast, the lack of a UV bump in typical DLAs (e.g., Khare et al. 2012) is probably intrinsic to their low dust and metal contents since the lines of sight are likely to pass at high impact parameters from the absorbing galaxy.

The main outlier in the right-hand panel of Fig. 12, which exhibits high reddening but no UV bump, is a proximate system. This is consistent with the above picture where the enhanced UV radiation field from the QSO and/or star-forming regions within the QSO host galaxy are expected to deplete the carriers of the 2175 Å absorption.

## 7. Conclusions

In this work, we presented a new population of QSO absorbers that was selected directly from the properties of the shielded gas, namely the strongest C1 absorbers, detected in low-resolution QSO spectra from the SDSS-II DR7 database. These C1 absorbers, with $W_r(\lambda_{1560}) \ge 0.4$ Å, are more than one hundred times rarer than DLAs at $z_{\text{abs}} = 2.5$. Their number per unit redshift is increasing significantly below $z_{\text{abs}} = 2$, probably coupled to an increase in the star-formation efficiency at these redshifts. Gupta et al. (2012) report a similarly high detection rate of 21-cm absorbers towards even lower redshifts among strong Mg II systems, which they argue must be related to the evolution of the CNM filling factor in the latter systems.

The HI column-density distribution of C1-selected absorbers is flatter than that of HI-selected absorbers. While sub-DLAs have a much larger cross-section than DLAs, this can be understood as the shielding of the gas being more difficult at

<footnote n="2">² Super-giant shells with sizes approaching 1 kpc form the largest structures seen in the ISM of galaxies where large amounts of kinetic energy are contributed by multiple supernovae explosions and energetic stellar winds.</footnote>
<figure><img src="image_13.png" /><figcaption>**Fig. 12.** Left-hand panel: reddening of the background QSOs with detected C I absorption (overall sample) versus C I λ1560 rest-frame equivalent width. Symbol conventions are the same as in Fig. 10. The data points circled in red correspond to the H I subsample. Right-hand panel: same as in the left-hand panel with the exception that reddening is plotted against the 2175 Å bump strength.</figcaption></figure>

low H I column densities and there are probably fewer clouds along the line of sight. Cold and dusty gas as traced by C I absorbers is also more likely to be found at low impact parameters from the absorbing galaxies where a flatter N(H I) distribution is expected. Indeed, despite a likely bias against strong DLAs with large amounts of dust, we find there is a probable excess of strong DLAs with log N(H I) > 21 (atoms cm⁻²) among C I systems compared to systematic DLA searches. This is reminiscent of the N(H I) distribution of DLAs within GRB host galaxies, which is skewed towards extremely strong DLAs (see Fig. 10 in Fynbo et al. 2009).

The reddening and therefore the presence of dust along the QSO lines of sight with detected C I absorption is directly related to the amount of shielded gas but depends weakly on the total H I column density. The latter can indeed vary by more than a factor of ten for the same C I rest-frame equivalent width. This is probably the consequence of the shielded gas being clumpy, while H I absorption samples simultaneously warm diffuse neutral clouds and cold, high-metallicity dusty pockets of gas. The presence of dust inducing significant reddening of the background QSOs and/or 2175 Å extinction features is ubiquitous in about 30% of the C I absorbers. Several systems like these have been found before (see, e.g., Srianand et al. 2008a; Wang et al. 2012). Here, we find that the UV bump is weak compared to Galactic lines of sight exhibiting the same amount of reddening. We interpret this as the consequence of star formation in the vicinity of the systems.

It is likely that the metal and molecular contents of C I absorbers are high and actually higher than those of most DLAs studied until now. High-resolution spectroscopic follow-up observations of the present sample therefore opens up the door to systematic searches for carbon monoxide (CO; see Noterdaeme et al. 2011) and molecules like CN and CH, as well as for diffuse interstellar bands at high redshift. Such a spectroscopic campaign will be presented in a companion paper. The typical reddening induced by C I absorbers along with the relation between reddening and shielded-gas column density imply that the

extinction could be high in some DLAs with C I absorption. If strong dusty DLAs exist, they probably have been missed in the current magnitude-limited QSO samples (see also Boissé et al. 1998; Vladilo et al. 2008). Some of the QSO lines of sight identified here, as well as those that may be found by extending the present survey to even larger databases<fnref n="3" />, will result in exceedingly long integration times on high-resolution spectrographs installed on 8–10 m class telescopes. These will, however, be targets of choice for the coming generation of extremely large telescopes.

*Acknowledgements.* P.N. acknowledges support from the ESO Chile visiting scientist programme. R.S. and P.P.J. gratefully acknowledge support from the Indo-French Centre for the Promotion of Advanced Research (Centre Franco-Indien pour la Promotion de la Recherche Avancée) under contract No. 4304-2. The authors of this paper also acknowledge the tremendous effort put forth by the Sloan Digital Sky Survey team to produce and release the SDSS survey. Funding for SDSS and SDSS-II has been provided by the Alfred P. Sloan Foundation, the Participating Institutions, the National Science Foundation, the US Department of Energy, the National Aeronautics and Space Administration, the Japanese Monbukagakusho, the Max Planck Society, and the Higher Education Funding Council for England. The SDSS Web Site is http://www.sdss.org/. The SDSS is managed by the Astrophysical Research Consortium for the Participating Institutions. The Participating Institutions are the American Museum of Natural History, Astrophysical Institute Potsdam, University of Basel, University of Cambridge, Case Western Reserve University, University of Chicago, Drexel University, Fermilab, the Institute for Advanced Study, the Japan Participation Group, Johns Hopkins University, the Joint Institute for Nuclear Astrophysics, the Kavli Institute for Particle Astrophysics and Cosmology, the Korean Scientist Group, the Chinese Academy of Sciences (LAMOST), Los Alamos National Laboratory, the Max-Planck-Institute for Astronomy (MPIA), the Max-Planck-Institute for Astrophysics (MPA), New Mexico State University, Ohio State University, University of Pittsburgh, University of Portsmouth, Princeton University, the United States Naval Observatory, the University of Washington.

<footnote n="3">³ From the Baryon Oscillation Spectroscopic Survey, which is part of SDSS-III, relatively few additional C I systems are expected since the bulk of the new QSOs is at $z_{em} \sim 3$, which provides shorter C I-absorption path length.</footnote>
References

Abazajian, K. N., Adelman-McCarthy, J. K., Agüeros, M. A., et al. 2009, ApJS,
182, 543

Balashev, S. A., Petitjean, P., Ivanchik, A. V., et al. 2011, MNRAS, 418, 357

Balashev, S. A., Klimenko, V. V., Ivanchik, A. V., et al. 2014, MNRAS, 440, 225

Boissé, P., Le Brun, V., Bergeron, J., & Deharveng, J.-M. 1998, A&A, 333, 841

Booth, R. S., de Blok, W. J. G., Jonas, J. L., & Fanaroff, B. 2009, ArXiv e-prints
[arXiv:0910.2935]

Budzynski, J. M., & Hewett, P. C. 2011, MNRAS, 416, 1871

Duffy, A. R., Meyer, M. J., Staveley-Smith, L., et al. 2012, MNRAS, 426, 3385

Finley, H., Petitjean, P., Pâris, I., et al. 2013, A&A, 558, A111

Fischera, J., & Dopita, M. 2011, A&A, 533, A117

Fitzpatrick, E. L., & Massa, D. 2007, ApJ, 663, 320

Fynbo, J. P. U., Jakobsson, P., Prochaska, J. X., et al. 2009, ApJS, 185, 526

Ge, J., & Bechtold, J. 1999, in Highly Redshifted Radio Lines, eds. C. L. Carilli, S. J. E. Radford, K. M. Menten, & G. I. Langston, ASP Conf. Ser., 156, 121

Gordon, K. D., Clayton, G. C., Misselt, K. A., Landolt, A. U., & Wolff, M. J. 2003, ApJ, 594, 279

Gudennavar, S. B., Bubbly, S. G., Preethi, K., & Murthy, J. 2012, ApJS, 199, 8

Gupta, N., Srianand, R., Petitjean, P., Noterdaeme, P., & Saikia, D. J. 2009, MNRAS, 398, 201

Gupta, N., Srianand, R., Petitjean, P., et al. 2012, A&A, 544, A21

Jenkins, E. B., & Tripp, T. M. 2011, ApJ, 734, 65

Jiang, P., Ge, J., Prochaska, J. X., et al. 2010, ApJ, 720, 328

Jiang, P., Ge, J., Zhou, H., Wang, J., & Wang, T. 2011, ApJ, 732, 110

Kanekar, N., Prochaska, J. X., Ellison, S. L., & Chengalur, J. N. 2010, ApJ, 712, L148

Khare, P., vanden Berk, D., York, D. G., Lundgren, B., & Kulkarni, V. P. 2012, MNRAS, 419, 1028

Klimenko, V. V., Balashev, S. A., Ivanchik, A. V., et al. 2015, MNRAS, 448, 280

Krühler, T., Ledoux, C., Fynbo, J. P. U., et al. 2013, A&A, 557, A18

Ledoux, C., Petitjean, P., & Srianand, R. 2003, MNRAS, 346, 209

Ledoux, C., Petitjean, P., Fynbo, J. P. U., Møller, P., & Srianand, R. 2006, A&A, 457, 71

Ledoux, C., Vreeswijk, P. M., Smette, A., et al. 2009, A&A, 506, 661

Madau, P., & Dickinson, M. 2014, ARA&A, 52, 415

Milutinovic, N., Ellison, S. L., Prochaska, J. X., & Tumlinson, J. 2010, MNRAS, 408, 2071

Muller, S., Combes, F., Guélin, M., et al. 2014, A&A, 566, A112

Noterdaeme, P., Ledoux, C., Petitjean, P., & Srianand, R. 2008a, A&A, 481, 327

Noterdaeme, P., Petitjean, P., Ledoux, C., Srianand, R., & Ivanchik, A. 2008b, A&A, 491, 397

Noterdaeme, P., Ledoux, C., Srianand, R., Petitjean, P., & López, S. 2009a, A&A, 503, 765

Noterdaeme, P., Petitjean, P., Ledoux, C., & Srianand, R. 2009b, A&A, 505, 1087

Noterdaeme, P., Petitjean, P., Ledoux, C., et al. 2010a, A&A, 523, A80

Noterdaeme, P., Srianand, R., & Mohan, V. 2010b, MNRAS, 403, 906

Noterdaeme, P., Petitjean, P., Srianand, R., Ledoux, C., & López, S. 2011, A&A, 526, L7

Noterdaeme, P., Petitjean, P., Carithers, W. C., et al. 2012, A&A, 547, L1

Pâris, I., Petitjean, P., Aubourg, É., et al. 2012, A&A, 548, A66

Pei, Y. C. 1992, ApJ, 395, 130

Petitjean, P., Ledoux, C., Noterdaeme, P., & Srianand, R. 2006, A&A, 456, L9

Petitjean, P., Srianand, R., & Ledoux, C. 2000, A&A, 364, L26

Petitjean, P., Srianand, R., & Ledoux, C. 2002, MNRAS, 332, 383

Pettini, M. 2006, in The Fabulous Destiny of Galaxies: Bridging Past and Present, eds. V. Le Brun, A. Mazure, S. Arnouts, & D. Burgarella, 319

Pitman, K. M., Clayton, G. C., & Gordon, K. D. 2000, PASP, 112, 537

Prochaska, J. X., Herbert-Fort, S., & Wolfe, A. M. 2005, ApJ, 635, 123

Prochaska, J. X., Hennawi, J. F., & Herbert-Fort, S. 2008, ApJ, 675, 1002

Prochaska, J. X., Madau, P., O'Meara, J. M., & Fumagalli, M. 2014, MNRAS, 438, 476

Quider, A. M., Nestor, D. B., Turnshek, D. A., et al. 2011, AJ, 141, 137

Rafelski, M., Wolfe, A. M., Prochaska, J. X., Neeleman, M., & Mendez, A. J. 2012, ApJ, 755, 89

Schlegel, D. J., Finkbeiner, D. P., & Davis, M. 1998, ApJ, 500, 525

Schneider, D. P., Richards, G. T., Hall, P. B., et al. 2010, AJ, 139, 2360

Snow, T. P., & McCall, B. J. 2006, ARA&A, 44, 367

Srianand, R., Petitjean, P., Ledoux, C., Ferland, G., & Shaw, G. 2005, MNRAS, 362, 549

Srianand, R., Gupta, N., Petitjean, P., Noterdaeme, P., & Saikia, D. J. 2008a, MNRAS, 391, L69

Srianand, R., Noterdaeme, P., Ledoux, C., & Petitjean, P. 2008b, A&A, 482, L39

Srianand, R., Gupta, N., Petitjean, P., et al. 2012, MNRAS, 421, 651

Vanden Berk, D. E., Richards, G. T., Bauer, A., et al. 2001, AJ, 122, 549

Vladilo, G., Prochaska, J. X., & Wolfe, A. M. 2008, A&A, 478, 701

Wang, J.-G., Zhou, H.-Y., Ge, J., et al. 2012, ApJ, 760, 42

Wolfe, A. M., Gawiser, E., & Prochaska, J. X. 2005, ARA&A, 43, 861

York, D. G., Adelman, J., Anderson, Jr., J. E., et al. 2000, AJ, 120, 1579

Zafar, T., Watson, D., Elíasdóttir, Á., et al. 2012, ApJ, 753, 82

Zwaan, M. A., & Prochaska, J. X. 2006, ApJ, 643, 675