On the deficiency of Argon in DLA systems

The ESO UVES Advanced Data Products Quasar Sample - IV. On the deficiency of Argon in DLA systems


In this work, we study argon abundances in the interstellar medium of high-redshift galaxies (  ) detected as Damped Ly absorbers (DLA) in the spectra of background quasars. We use high-resolution quasar spectra obtained from the ESO-UVES advanced data products (EUADP) database. We present 3 new measurements and 5 upper limits of Ar i. We further compiled DLAs/sub-DLA data from the literature with measurements available of argon and -capture elements (S or Si), making up a total of 37 systems, i.e. the largest DLA argon sample investigated so far. We confirm that argon is generally deficient in DLAs, with a mean value [Ar/]  dex (standard error of the mean). The [Ar/] ratios show a weak, positive trend with increasing H i column density and increasing absorption redshift, and a weak, negative trend with dust-free metallicity, [S/H]. Detailed analysis of the abundance ratios indicates that Ar i ionisation, rather than dust depletion or nucleosynthetic evolution, is responsible for the argon deficiency. Altogether, the observational evidence is consistent with a scenario of argon ionisation dominated by quasar metagalactic radiation modulated by local H i self-shielding inside the DLA host galaxies. Our measurements and limits of argon abundances suggest that the cosmic reionisation of He ii is completed above , but more measurements at   3.5 are required to probe the final stages of this process of cosmic reionisation.

Galaxies: formation – galaxies: evolution – galaxies: abundances – galaxies: ISM – quasars: absorption lines – intergalactic medium

1 Introduction

High-resolution quasar spectra show a large number of absorption lines produced in the diffuse gas intercepted by the line of sight, from redshift up to the quasar’s redshift. Most quasar absorption lines originate in the intergalactic medium (IGM) and contain information on the physical properties of a large volume of the Universe over a large fraction of look-back time. A small fraction of quasar absorption systems, those with the highest values of neutral hydrogen column density, N(H i) atoms cm, originate in the diffuse gas of intervening galaxies. These absorbers, called Damped Ly absorbers (DLAs) (Wolfe et al., 1986, 2005), contain information on the physical state and chemical composition of the interstellar medium (ISM) in high-redshift galaxies. By studying DLAs we can cast light on the astrophysical processes that take place in their host galaxies and, in particular, on the interaction between the metagalactic radiation propagating through the IGM and the interstellar gas at high redshift. In this work, we focus our attention on a particular type of absorption lines in DLAs, the absorption of neutral argon, which offers a way to probe this interaction.

Owing to the high value of its first ionisation potential ( eV), argon is expected to be mostly neutral in interstellar H i regions, which are opaque to photons with energies just above the H i ionisation threshold ( eV). Nevertheless, the ionisation fraction of argon is very sensitive to high energy ionising photons that are able to leak through the neutral gas. This sensitivity is due to the high ratio of photoionisation to recombination rates, typically one order of magnitude higher for Ar i than for H i (Sofia & Jenkins, 1998). The fraction of Ar i is expected to vary according to the energy distribution of the radiation field. When the radiation field is hard, Ar i is predicted to be deficient relative to other “low ions” typical of H i regions. On the other hand, normal (i.e. solar) abundances are expected when the ionising spectrum is soft. In our own Galaxy, neutral argon is found to be deficient in low-density interstellar regions (Sofia & Jenkins, 1998; Jenkins et al., 2000; Lehner et al., 2003; Jenkins, 2013). However, owing to the saturation of Ar i lines, Galactic studies can only probe regions with relatively low H i column density, typically associated to warm interstellar gas. On the other hand, thanks to the low metallicity of DLA galaxies (e.g. Wolfe et al., 2005), Ar i lines can be unsaturated in DLAs in spite of the high column densities of these absorbers, N(H i) atoms cm (Wolfe et al., 1986). Therefore, Ar i lines offer a way to probe the ionisation state of the neutral gas in high-redshift galaxies detected in absorption. In principle, by measuring the ionisation state of the gas at different redshifts, we can track the evolution of the comoving density of cosmic sources and variations of the transparency of the IGM resulting from the processes of cosmic reionisation.

According to our current understanding of the reionisation history of the Universe (Fan et al., 2006; McQuinn et al., 2009; Jarosik et al., 2011; Haardt & Madau, 2012), hydrogen was reionised by the ultraviolet radiation of the first stars and quasars, starting at ; with the increasing overlap of IGM H ii regions, hydrogen reionisation was completed around , leaving the IGM completely transparent to photons with  eV. During the early stages of H i reionisation, the ionising radiation field was sufficiently hard to ionise He i ( eV), but not He ii ( eV). At a later stage, when the density of hard UV-emitting quasars became sufficiently high, also He ii was reionised. This process possibly started around (Bolton et al., 2012), leading to the complete overlap of IGM He iii regions around (Furlanetto & Oh, 2008; McQuinn et al., 2009). With the completion of He ii reionisation the IGM became completely transparent to photons with  eV.

Probing the processes of cosmic reionisation requires a variety of experimental tools. Observations of Ar i in DLAs provide a rare possibility to probe the effects of internal and external radiation fields acting inside high-redshift galaxies. Since the stellar radiation originated inside the galaxy is softer than that of the quasar background, one can in principle discriminate the ionisation sources by taking advantage of the sensitivity of Ar i to high energy photons. This type of study helps to understand the last stages of cosmic reionisation that take place in a redshift interval where DLAs can be easily observed.

Neutral argon can be measured from the resonant transitions at 1048 Å and 1066 Å. In DLAs these lines are observable from ground in the visible and near ultraviolet when the absorption redshift is   2. Argon lines are difficult to detect because they can be weak in the DLAs of lower metallicity and because they fall in the Ly forest of the spectrum. The first measurement of Ar i abundance in a DLA system, namely the absorber at =3.39 towards Q 0000-2620, did not show evidence for a deficiency of argon (Molaro et al., 2001). The first systematic study, based on a sample of 10 measurements and 5 limits of Ar i, indicated that the absorber at =3.39 towards Q 0000-2620 is exceptional, since argon is deficient in most DLAs (Vladilo et al. 2003; hereafter Paper I). An evolutionary trend between argon deficiency and redshift was found in the redshift interval covering the epoch of He ii reionisation. However, the limited size of the sample made difficult to derive firm conclusions on this and other trends of argon abundances.

In the present work we reassess the current status of argon abundances in DLAs after the accumulation of many years of observations obtained with the European Southern Observatory (ESO) Ultraviolet-Visual Echelle Spectrograph (UVES). This spectrograph, optimized for the near ultraviolet (Dekker et al., 2000) and fed by the Unit 2 of the Very Large Telescope (VLT), is an ideal instrument for this type of research. To enlarge the sample of argon measurements, we use the Advanced Data Products (EUADP) sample (Zafar et al., 2013b, a), together with an updated compilation of literature data. The new sample and the analysis of the data are presented in §2 and notes on new discoveries are provided in §3. In §4 we present the results, showing that argon is generally deficient in DLAs §5 we discuss the nature of argon deficiency and its possible causes and in §6 we discuss the nature of the ionising sources and the evidence of evolution. The conclusions of our work are provided in §7.

2 Data sample

2.1 The EUADP data

This work is based on high-resolution () quasar spectra obtained from the ESO-UVES instrument and collected in the EUADP sample. Details on building and processing of the EUADP sample are presented in Zafar et al. (2013b). This EUADP sample of 250 quasar spectra has been used to study DLAs/sub-DLAs (Zafar et al., 2013a). Indeed, 197 DLAs and sub-DLAs have been reported in these quasar spectra, thus providing the ideal parent sample for the study of rare elements in quasar absorbers. In particular, the EUADP sample allows for the observations of Ar i at redshifts 1.86 4.74. In this work, we search for such Ar i lines in 140 quasar absorbers. Ar i is not detected in more than one-third of the cases. This is due to the lack of wavelength coverage in the relevant region, the presence of a Lyman Limit System, or low signal-to-noise ratio in the region of interest. Moreover, 9 measurements, 1 upper and 1 lower limit of Ar i have been already reported in the literature for the EUADP DLAs/sub-DLAs (see Table 1). Of the remaining systems, more than half were rejected due to strong blending of Ar i with the lines from the Ly forest. We find 3 new detections of Ar i never reported before. In addition to these measurements, we report 5 new upper limits of Ar i in this EUADP DLAs/sub-DLAs sample.

2.2 Fitting method

The ionic column densities (Ar i, S ii, Si ii, etc) were derived using the -minimization routine FITLYMAN within the MIDAS environment (Fontana & Ballester, 1995). Laboratory wavelengths and oscillator strengths were taken from Morton (2003). The global fit returns the best fit parameters for the central wavelength, column density and Doppler turbulent broadening, as well as 1 errors on each of these quantities.

Given that the two lines of Ar i transitions are located in the Ly forest, great care was taken to assess the contamination of the lines by interlopers. This possible blending would tend to overestimate the column density of Ar i. A comparison with radial velocity profiles of low ionisation species (S ii, Si ii, and/or Fe ii) has been made to uncover the transitions which are free from contamination. These low ionisation species were chosen as reference absorption lines because they lie outside the Ly forest and are thus free from possible blending. The agreement with the radial velocity profiles of the Ar i lines is generally satisfactory. Upper limits were derived in the cases with evidences for contamination of Ar i transitions by Ly forest interlopers. This is done by fixing the and parameters of Ar i lines to those of the other detected low ionisation lines and fitting for column density. The details of each individual system are given in §3.

2.3 Literature data

In addition to the measurements presented here, we made an investigation in the literature for quasar absorbers with Ar i abundance determination. We find 17 measurements, 11 upper and 1 lower limit of Ar i. Incidentally, some of the objects are part of the EUADP DLA/sub-DLA sample. The combined sample makes a total of 37 DLAs/sub-DLAs (20 measurements together with 16 upper and 1 lower limit). The targets, the basic data for the DLAs/sub-DLAs investigated, and Ar i total column densities are listed in Table 1 (see also Table 2).

Sample Measurements Upper limits Lower limits
Literature 17 11 1
Table 2: Summary of the number of systems with detected Ar i in the two samples.

3 Notes on individual objects

In the following, we briefly describe each system for which we report here for the first time the Ar i abundance. Additional information on the UVES spectra and the determination of H i column densities for these absorbers can be found in Zafar et al. (2013b, a).

J 1113-1533, =3.2655: this DLA has an H i column density of log N(H i) . The four-component structure of this system is derived from the fit to the S ii  1250, 1259 absorption lines. The argon abundance is derived from unsaturated, non-blended Ar i  1066 line (see Fig. 1). The redshift of the components, Doppler parameters and column densities are listed in Table 3. The Si ii lines are either saturated or blended. We derive a lower limit from Si ii  1304, 1526 transitions: log N(Si ii) by using the redshift and -value of the components obtained from the S ii transitions.

Comp. Ion log N
km s cm
1 3.2653 ii
Ar i
2 3.2655 ii
Ar i
3 3.2657 ii
Ar i
4 3.2658 ii
Ar i
Table 3: Results of voigt profile decomposition of the =3.2655 DLA towards J 1113-1533.
Figure 1: Voigt-profile fits (red overlay) of Ar i and S ii for the DLA at =3.2655 (zero velocity) in J 1113-1533. Normalised quasar spectrum and error spectrum are shown in black and blue dotted lines, respectively. The blue overlay represents the profile of the transition which is not fit and obtained from the fit of the other transitions in the system. The locations of the voigt profile components

are indicated by the red tick marks.

J 1356-1101, =2.9669: this DLA has an H i column density of log N(H i) . The four-component structure of the system is derived from the fit to the Si ii  1808 absorption line. The argon abundance is derived from unsaturated, non-blended Ar i  1066 line (see Fig. 2). The fit parameters are described in Table 4. Note that Noterdaeme et al. (2008) derived log N(Si ii), consistent with our results within the errors.

Comp. Ion log N
km s cm
1 2.9666 Si ii
Ar i
2 2.9668 Si ii
Ar i
3 2.9669 Si ii
Ar i
4 2.9670 Si ii
Ar i
Table 4: Results of voigt profile decomposition of the =2.9669 DLA towards J 1356-1101.
Figure 2: Voigt-profile fits (red overlay) of Ar i and Si ii for the DLA at =2.9669 (zero velocity) in J 1356-1101.

B 2332-094, =3.0572: this sub-DLA has an H i column density of log N(H i) . The metal profile is spread over several components. However, unsaturated, non-blended Ar i  1048 line is seen in only three of the strongest components (see Fig. 3). This three-component structure is used to fit the S ii  1259, Fe ii  1608, and Si ii  1808 absorption lines. The parameters for the fit are summarised in Table 5. Note that Petitjean et al. (2008); Prochaska et al. (2003) derived log N(S ii), log N(Fe ii), and log N(Si ii) by fitting all the components in this system. Logically, our determinations of the total column densities are slightly lower than the one reported by other authors.

Comp. Ion log N
km s cm
1 3.0572 ii
Si ii
Fe ii
Ar i
2 3.0573 ii
Si ii
Fe ii
Ar i
3 3.0574 ii
Si ii
Fe ii
Ar i
Table 5: Results of voigt profile decomposition of the =3.0572 DLA towards B 2332-094.
Figure 3: Voigt-profile fits (red overlay) of Ar i, S ii, Fe ii, and Si ii (see labels) for the DLA at =3.0572 (zero velocity) in B 2332-094.

4 Argon abundances in DLA systems

In this section we discuss the elemental abundances of argon obtained from the measurements presented in Table 1. We follow the common convention expressing the number abundance ratios of two elements X and Y according to the relation , where is the ratio of the column densities measured in the gas phase and is a reference abundance based on the solar composition. Here we adopt the set of solar photospheric abundances of Asplund et al. (2009).

Unfortunately, it is difficult to obtain accurate solar abundances for “noble gas” elements such as argon. The value that we adopt here, , is an average of results obtained from a variety of techniques (Asplund et al., 2009). In Paper I, we adopted the solar abundances of Grevesse & Sauval (1998) for most elements and the argon solar value proposed by (Sofia & Jenkins, 1998), , which is marginally consistent with the value proposed by Asplund et al. (2009). Taking into account these differences, the [Ar/S] and [Ar/Si] ratios that we calculate now are higher by +0.04 dex and +0.07 dex, respectively, with respect to those calculated in Paper I. The uncertainty of the reference value affects the level of argon abundance that we derive, but not the correlation trends that we present below.

Figure 4: The [Ar/S] and [Ar/Si] versus S-based metallicities in DLA systems. Solid squares and arrows correspond to measurements and limits, respectively. The open squares represent proximate DLAs. A small chemical evolution could be seen in both [Ar/S] and [Ar/Si].

4.1 Argon deficiency in DLA systems

In our sample, the absolute abundance of argon (i.e., the abundance relative to hydrogen) spans the range [Ar/H] , with a mean value  dex (standard error of the mean; 20 measurements). Broadly speaking, this result is in line with the fact that DLA systems are metal poor. However, the mean value [Ar/H] is lower than the mean metallicity of DLA systems, [M/H]  dex (Wolfe et al., 2005), indicating that argon is underabundant with respect to other metals. To evidentiate this underabundance, we use relative abundances, i.e., the abundances of argon relative to other products of stellar nucleosynthesis. By using relative abundances, we can offset differences in the level of metallicity attained by individual DLA galaxies, and we can focus our attention on effects other than metallicity, such as dust and ionisation effects. In practice, we use sulphur and silicon as reference elements to measure the relative abundances. The existence of a relatively large number of S ii and Si ii column density measurements in DLA systems is one reason for choosing these two elements. In addition, silicon and sulphur have nucleosynthetic history similar to that of argon and negligible (or low) dust depletion (see §5).

From the analysis of our sample of measurements, rather than limits, we obtain the following results. The relative abundance with respect to sulphur spans the range [Ar/S] , with mean value  dex (standard error of the mean; 18 systems). Using silicon as a reference, we find values [Ar/Si] , with a mean value  dex (17 systems).

The upper limits are consistent with these values and in most cases are stringent enough to confirm the general deficiency of argon. Therefore, the results obtained from the updated Ar i sample confirm the general deficiency of argon in DLA systems found in Paper I.

The standard deviation of the argon abundance ratios is  dex and  dex for [Ar/S] and [Ar/Si], respectively. The spread of these abundance ratios is somewhat higher than that of other relative abundances in DLA systems. For instance, we find a scatter of 0.14 dex and 0.17 dex for [Zn/S] and [Si/Fe], respectively, for the current sample of DLA systems with column densities of S ii, Zn ii, Si ii and Fe ii obtained from high-resolution spectra.

Proximate DLA systems

A few absorption systems of our sample are proximate DLAs, i.e. their relative velocity separation from the background quasar is  km s (e.g., Møller et al., 1998; Ellison et al., 2010; Zafar et al., 2011). Specifically, these are the absorber at towards B 0528-2505, at towards J 13373152, and that at towards J 1604+3951. These systems will be treated separately in the subsequent analysis since they may not be representative of the typical population of DLA systems (Ellison et al., 2010). With the exclusion of these three cases, the mean values and standard errors of the mean are  dex (15 systems) and  dex (14 systems), confirming that argon is significantly deficient in classic DLA systems.

4.2 Trends of argon abundances

The relatively large sample of the present study allows us to reassess the presence of trends between the relative abundances of argon and other quantities that we can measure in DLAs. In particular, we review the evidence for trends with metallicity, hydrogen column density and absorption redshift. In Tables 6 and 7 we provide the results of a correlation analysis performed for the complete sample and the sample without proximate DLAs, respectively. In these Tables we give the Pearson’s correlation coefficient, and the Spearman’s rank correlation coefficient, , obtained from the datasets of [Ar/S] and [Ar/Si] measurements. We also give the probability of null hypothesis, and obtained for the Pearson’s and Spearman’s tests, respectively. The -values roughly indicate the probability of an uncorrelated population producing datasets that have a Pearson or Spearman correlation at least as extreme as the one computed from the measured datasets.

Trend with metallicity

The absorption lines of the S ii triplet at 1250 Å are usually detected in the quasar spectra where Ar i is detected. For this reason, and because sulphur is not affected by dust depletion (see §5.1), we use the absolute abundance of sulphur, [S/H], as a metallicity indicator in our study.

In Fig. 4 we plot our measurements and limits of [Ar/S] and [Ar/Si] versus [S/H]. Visual inspection shows that argon becomes more deficient as the metallicity increases. At the lowest values of metallicity of our sample, the relative abundances are approximately solar, whereas at the highest metallicities they are deficient by at least  dex with respect to the solar value.

Some evidence for an anti-correlation between [Ar/S] and metallicity is present in both samples (Tables 6 and 7), with correlation coefficients and null hypothesis probabilities . A similar trend is found for [Ar/Si], albeit with a weaker correlation, (; /). The Spearman tests yield similar results. The existence of a trend between argon abundances and metallicity in DLAs is reported here for the first time.

Trend with hydrogen column density

In Fig. 5 we plot our measurements and limits of [Ar/S] and [Ar/Si] versus log N(H i). Visual inspection shows that argon has a tendency to become more abundant as the H i column density increases, albeit with a large spread. At the lowest H i column densities of our sample the relative abundances are deficient by more than  dex , whereas at the highest column densities they may attain approximately solar values.

The statistical analysis of [Ar/S] versus log N(H i) for the complete sample indicate some evidence for a positive correlation (; ). The evidence is stronger for the sample without proximate DLAs (; ). A positive correlation is also present for [Ar/Si], but less robust: we find () for the complete sample, and () for the sample without proximate DLAs. The Spearman tests yield similar results, with slightly larger null-hypothesis probabilities for the [Ar/Si] versus log N(H i) trends.

In the sample of Paper I there was no statistical evidence for the existence of a correlation with N(H i), even though the DLA with highest N(H i) did show the highest argon abundance.

Trend with absorption redshift

In Fig. 6 we plot our measurements and limits of [Ar/S] and [Ar/Si] versus . Visual inspection shows that argon has a tendency to become more abundant as the redshift increases, but with a large spread. The statistical analysis of [Ar/S] versus  shows a very weak evidence for a positive correlation (; ) for both samples (with and without proximate DLAs). A similar, positive trend is found for the [Ar/Si] ratios, with a slightly higher correlation: we find () for the complete sample and () for the sample without proximate DLAs. The Spearman tests yield slightly lower correlation coefficients and larger probabilities of null-hypothesis. The weakness of the correlations with absorption redshift contrasts with the high confidence level found in Paper I for the same type of analysis. The larger size of the present sample reveals a more complex situation, with a modest trend and a large scatter.

[Ar/S] [S/H] 18 0.5% 0.3%
[Ar/Si] [S/H] 15 6.3% 5.6%
[Ar/S] log N(H i) 18 1.8% 1.8%
[Ar/Si] log N(H i) 17 14.7% 21.2%
[Ar/S] 18 40.9% +0.13 60.7%
[Ar/Si] 17 16.4% +0.21 41.4%
Table 6: Correlation analysis of argon abundances for the complete sample of measurements in DLAs.
[Ar/S] [S/H] 15 1.7% 1.1%
[Ar/Si] [S/H] 12 14.2% 10.6%
[Ar/S] log N(H i) 15 0.5% 0.5%
[Ar/Si] log N(H i) 14 6.3% 10.4%
[Ar/S] 15 40.1% +0.13 64.8%
[Ar/Si] 14 12.5% +0.22 45.4%
Table 7: Correlation analysis of argon abundances for the sample without proximate DLAs.

5 Nature of the argon deficiency

The abundance of argon, as any other elemental abundance measured in DLAs, is influenced by three types of effects at work in the diffuse gas of the host galaxy: dust depletion, nucleosynthesis and ionisation processes. Disentangling the relative contribution of these effects is crucial for a correct interpretation of the measured abundances.

5.1 Dust Depletion

Figure 5: The [Ar/S] and [Ar/Si] ratio in DLA systems plotted versus H i column density. Solid squares and arrows correspond to measurements and limits, respectively. The open squares illustrate proximate DLAs. Curves correspond to the predictions of idealized CLOUDY photoionisation models calculated for solar abundance ratios. Dashed curve: model predictions obtained for a layer of gas at atoms cm embedded in the HM extragalactic radiation field (e.g. Haardt & Madau, 2012); dotted curve: predictions obtained by doubling the intensity of the radiation field in the same model. See §5.3.

Dust depletion is due to the incorporation of part of the atoms of the ISM into solid phase (dust grains). Elements classified as volatiles and refractories according to condensation calculations (Lodders, 2003) are characterized by low and high levels of depletion, respectively (Savage & Sembach, 1996). Dust depletion in DLAs is known to be moderate, similar to that found in low-density components of the Galactic ISM (Pettini et al., 1994, 1997; Ledoux et al., 2002; Vladilo et al., 2011). For this reason, depletions of well-known volatile elements, such as zinc or sulphur, or even moderately refractory elements, such as silicon, have negligible or modest effects on DLA abundance measurements. However, depletions of refractory elements, such as iron, do affect the abundance measurements and will not be considered in our work, which is focussed on ionisation effects. Argon, as a “noble gas”, has little tendency to take part to chemical reactions, including those that lead to the condensation of dust grains. Detailed calculations of condensation temperatures indicate, indeed, that argon is an extremely volatile element (Lodders, 2003). In principle, argon atoms could form the stable hydride ArH through the reaction Ar+H ArH+H (Duley, 1980; Sofia & Jenkins, 1998). However, this reaction requires molecular hydrogen, but the fraction of molecular hydrogen, i, is generally small or negligible in DLAs (see e.g., Srianand et al., 2012). The absorber at towards B 1232+0815 is the only system of our sample with (Noterdaeme et al., 2008). This systems shows a modest argon depletion (see Table 1), suggesting that the molecular route for argon depletion may require a high fraction of molecular gas. Therefore, we do not expect the ArH production to play a role in the other molecular systems of our sample, which have at most (see Table 1).

To test whether argon is depleted or not from an experimental point of view, in Fig. 7 we plot the ratio [Ar/S] versus [Zn/Fe]. The ratio [Zn/Fe] is a well known indicator of dust depletion in the ISM (e.g., Savage & Sembach, 1996) commonly employed in studies of DLA systems (e.g., Pettini et al. 1994, 1997; Ledoux et al. 2002; Ellison & Lopez 2009; Vladilo et al. 2011; see however Rafelski et al. 2012). If argon is depleted into dust, we would expect a trend of decreasing [Ar/S] with increasing [Zn/Fe], assuming sulphur undepleted. Fig. 7 shows instead a large scatter without a clear trend. A linear regression analysis of [Ar/S] versus [Zn/Fe] for the sample without proximate DLA systems yields a Pearson’s correlation coefficient (; 8 systems). Therefore, even if the statistics are too small to draw firm conclusions, there is no evidence of trend.

Based on the theoretical considerations discussed above, and given the lack of a trend with [Zn/Fe], in the rest of the discussion we assume that the dust depletion of argon is negligible.

5.2 Nucleosynthesis

The abundance ratios measured in DLAs tend to change in the course of the evolution of the host galaxies as a result of the enrichment of the diffuse gas by heavy elements produced by stellar nucleosynthesis and ejected by supernova explosions. Studies of Galactic chemical evolution indicate that the relative abundances of elements with a similar nucleosynthetic origin, such as -capture elements, undergo a modest evolution as the global level of metallicity increases. Therefore, relative abundance changes due to ionisation effects for the -elements S, Si and Ar considered here are expected to be small.

Unfortunately, to our knowledge, there are no measurements of argon abundances in the samples of metal-poor stars used to trace the effects of Galactic chemical evolution. Argon abundances can be measured in early-type stars (e.g., Lanz et al., 2008), but these stars are short-lived and therefore not suitable for tracing the chemical evolution of the Galaxy. Observational studies of argon chemical evolution are based on measurements of emission lines from H ii regions rather than on stellar measurements. By comparing measurements obtained in extragalactic H ii regions of different metallicity, Henry & Worthey (1999) find that abundance ratios Ne/O, S/O, and Ar/O appear to be universally constant and independent of metallicity. Similar conclusions are obtained from a study of O, Ar and S abundances in metal-poor H ii regions Izotov & Thuan (1999). From the analysis of the data published in Table 1 of that paper, and adopting the solar reference abundances of Asplund et al. (2009), we find a mean value  dex, without evidence of evolution of [Ar/S] with metallicity [S/H]. Based on these observational results, we conclude that nucleosynthetic effects are unable to explain the general under-abundance of [Ar/] that we observe in DLAs.

Figure 6: The [Ar/S] and [Ar/Si] ratios in DLA systems plotted versus redshift (bottom axis) and loockback time (top axis), for km s Mpc, , and . The open squares correspond to proximate DLAs.

5.3 Ionisation effects

The argon abundance ratios are affected by photoionisation processes induced by the radiation field that pervades the diffuse gas of DLA galaxies. The ionisation ratios measured in DLA systems can be modelled with ionisation equilibrium codes, such as CLOUDY (Ferland et al., 1998), assuming plane parallel geometry with photons incident on one side of the cloud. We used CLOUDY (version 13.03; Ferland et al. 2013) to perform two sets of calculations: one with the Haardt & Madau (HM) extragalactic background (Table HM05, see e.g. Haardt & Madau 2012) at , and one with the radiation field of a star with  K (Kurucz, 1991), representative of stellar sources internal to the DLA host galaxy. In all calculations we adopted metal abundances 0.1 the solar reference values (Asplund et al., 2009). We included interstellar grains for heating and cooling processes with a dust-to-gas ratio 0.1 of the Galactic value, as well as cosmic ray heating typical of Galactic interstellar gas. We adopted a total hydrogen density atoms cm, a value representative of Galactic warm gas. The calculations were stopped at an assigned value of H i column density. The highest values of H i column density are difficult to model with CLOUDY since the deepest layers tend to become cooler than 1000 K, the default value of the temperature threshold which is present in the code. To increase the possibility of tracking high density layers, we lowered the temperature threshold to 100 K. To constrain the model, we required the predicted hydrogen ionisation to be and the Al iii/Al ii ratio to match the observed data. To this end, we searched for measurements of Al iii and Al ii column densities for the DLAs of our sample, checking that the absorption lines of both ions shared the same velocity profiles as the Ar i lines. Unfortunately, only a few DLAs in our sample have Al iii and Al ii measurements3. Based on these data, we adopted log Al iii/Al ii  dex as a reference value. The fractions of Ar i, Si ii and S ii predicted by the ionisation models were used to calculate the [Ar/S] and [Ar/Si] ratios.

In the case of the extragalactic radiation field, the model was tuned by scaling the intensity of the HM field. Remarkably, the model is quite successful in matching the experimental data even taking the HM field at face value. At log N(H i)=20.75 dex the model yields , log Al iii/Al ii= dex, [Ar/S]  dex and [Ar/Si]= dex. The abundance of neutral argon tends to increase with increasing H i column density. The predictions of [Ar/S] and [Ar/Si] versus log N(H i) are shown with a dashed line in Fig. 5. One can see that the model yields values of [Ar/S] and [Ar/Si] in the range of the observations and, in addition, predicts the existence of a trend between argon deficiency and H i column density consistent with the observed one (Fig. 5). The scatter of the experimental [Ar/S] and [Ar/Si] ratios at any value of N(H i) can be modelled, to some extent, invoking local variations of the extragalactic radiation field and of the hydrogen density of the cloud. As an example, we plot in Fig. 5 the results obtained by doubling the intensity of the HM radiation field (dotted line). Similar results are obtained by decreasing by a factor of two. Variations of the cloud geometry, that cannot be accounted for by CLOUDY, may certainly introduce additional scatter.

In the case of the stellar radiation field, the model was tuned by varying the ionisation parameter4. However, with the stellar field CLOUDY is not able to finish the calculations because the gas temperature falls below the minimum allowed value before the total N(H i) is attained. To bypass this problem, Vladilo et al. (2001) considered a two-region model composed of a mildly ionised interface and an inner neutral region that only hosts dominant ionisation species. With this type of model, the CLOUDY calculations can be performed in the interface, thanks to its low H i column density, and the ionisation ratios can be calculated by averaging the column densities over the two regions (see Vladilo et al. 2001 for more details). No significant argon underabundance was found with this model tuning the stellar field to match the measured Al iii/Al ii ratios (Vladilo et al., 2003). In the present work we have explored the predictions of the stellar field model as a function of depth in the cloud in the range of H i column densities for which the CLOUDY calculations are not interrupted. We find that the ionisation of Ar is negligible in the layers where the Al iii/Al ii ratio attains values comparable to the observed ones.

In conclusion, ionisation processes offer a key for understanding the origin of the argon deficiency, at variance with dust depletion and nucleosynthetic evolution effects. The general level of argon deficiency and the trend with H i column density can be explained by photoionisation models with radiation field dominated by extragalactic background. The fact that the stellar field is not able to produce a significant underabundance of argon must be related to the paucity of stellar photons with energies of several tens of eV. In this energy range5, where the spectrum of the extragalactic background does not decline, argon ionisation is particularly efficient because the ratio of photoionisation-to-recombination rates of Ar i is high (see, e.g., Sofia & Jenkins, 1998, Fig. 3).

6 Sources of argon ionisation

Based on the above discussion, we now assume that the measured Ar i deficiency is due to ionisation and we investigate the nature of the ionising sources.

6.1 External sources

The comparison between observational data and model predictions in Fig. 5 suggests that the metagalactic UV/X-ray background impinging on the DLA host galaxies is the most natural source of ionisation. The trend of Ar abundances with redshift (§4.2) provides an observational test for probing the hypothesis of external ionisation.

Interpretation of the trend with redshift

In the redshift interval typical of most DLA observations () the hydrogen ionisation rate in the IGM has a weak tendency to decrease with redshift, albeit with large uncertainties (see e.g. left panel of Fig. 8 in Haardt & Madau, 2012). The CLOUDY model with HM extragalactic radiation (§5.3) predicts the [Ar/S] and [Ar/Si] ratios to increase by 0.09 dex and 0.13 dex, respectively, when the redshift increases from to (case log N(H i)  dex). The measurements shown in Fig. 6 and the correlation analysis in Tables 6 and 7 are consistent with a modest rise of argon abundances with redshift, in line with the model predictions. The large scatter observed at any given redshift can be naturally explained in terms of attenuation of the external field by H i absorption inside the cloud. This possibility is supported by the trend between argon deficiency and H i column density (Fig. 5), which indicates that local absorption does affect the level of ionisation. A possible example of local attenuation is provided by the two DLAs6 with small [Ar/S] deficiency and low redshift that are visible in the left panel of Fig. 6. The high H i column density of these two absorbers may provide a strong H i self-shielding, attenuating argon ionisation in a redshift interval where the external field is expected to be relatively strong. If we exclude from our sample all DLAs with log N(H i) , the Pearson’s correlation coefficient of [Ar/S] versus increases from =+0.21 () to =+0.34 (). These considerations suggest that a weak trend with redshift may be present, but it is probably masked by local H i self-shielding.

Figure 7: The [Zn/Fe] ratios against [Ar/S] abundance ratios in DLA systems. The open squares represent proximate DLAs.

Probing the end of cosmic He ii reionisation

In the upper range of the DLA redshift interval the cosmic reionisation of He ii is expected to be completed (Furlanetto & Oh, 2008; McQuinn et al., 2009). When cosmic He iii regions become completely overlapped, the photons with  eV should travel undisturbed through the IGM, boosting the efficiency of argon ionisation processes. The measurements versus redshift shown in Fig. 6 do not show a sudden rise of ionisation level that could be associated with the completion of He ii reionisation. Revealing this effect is quite difficult given the low statistics of DLA systems at high redshift. Our new data sample includes a total of 8 Ar i measurements at (5 more than in Paper I), but only one of these lies at . Therefore, the redshift coverage of our sample is still marginal if one wants to follow the evolution of He ii reionisation, which apparently starts around (Bolton et al., 2012) and is completed around (Dall’Aglio et al., 2008; Syphers et al., 2009; McQuinn et al., 2009). At our statistics is relatively large and suggests that hard photons are able to ionise argon, with the exception of the most self-shielded cases mentioned above. This result is consistent with a scenario in which the IGM becomes transparent to hard photons at some stage above . However, at high redshift, say , we have only two absorbers to probe earlier stages, and they give contrasting results. The DLA at towards Q 0000-2620 (Molaro et al., 2001) is characterized by nearly solar argon abundances, whereas the DLA at towards J 1443+2724 (Noterdaeme et al., 2008) shows evidence of argon deficiency. All together, our results are consistent with a scenario in which the He ii reionisation is completed above , but a larger statistics of data at is required to probe the expected rise of IGM photons with  eV.

6.2 Testing the hypothesis of external ionisation

The comparison of the argon abundances in DLA systems with similar measurements in proximate DLA systems and in the Galactic ISM provide ways of testing the hypothesis that the ionising field is external.

Clues from proximate DLA systems

If external ionisation dominates, we expect to find a strong argon deficiency in the proximate DLAs of our sample due to their vicinity to a quasar source of external hard photons. In Fig. 5 we plot with different symbols the argon abundances versus N(H i) for the proximate DLAs (open squares) and the classic DLAs of our sample (filled squares). One can see that the three proximate absorbers have a stronger argon deficiency compared to the classic DLA systems of similar N(H i). This suggests that, for a given level of H i self-shielding, the resulting argon ionisation is stronger in the proximate DLAs due to the vicinity of the quasar. In Fig. 6 one can see that the same three proximate absorbers also have a stronger argon deficiency at a given redshift, compared to classic DLA systems. This suggests that, for a given level of metagalactic background at a specific redshift, their ionisation is stronger due to excess radiation of the nearby quasar. The analysis of the upper limits of [Ar/], rather than the measurements, is consistent with these results. In fact, the Ar i upper limit for the proximate DLA system at towards J 1131+6044 (Ellison et al., 2010) yields the strongest underabundance of any other absorber of our sample, with [Ar/Si]  dex.

The results of these tests are in line with the hypothesis that argon ionisation is very sensitive to external radiation with a quasar-like spectrum. Clearly, it is necessary to extend this test to a larger sample in order to obtain a robust conclusion.

Figure 8: The [S/H] ratios versus H i column density in DLA systems. The open squares represent proximate DLAs.

Clues from Galactic interstellar studies

A recent study of argon abundance in the Galactic ISM (Jenkins, 2013) indicates an average argon deficiency [Ar i/O i]  dex, with the argon solar abundance7 adopted here (Asplund et al., 2009). The mean [Ar/] that we measure in DLA systems is  dex lower than that measured in the ISM. The stronger underabundance in DLA systems is remarkable, considering that the H i column densities of the Galactic lines of sight with Ar i measurements are significanty lower8 than the H i column densities in DLA systems. Therefore, H i shielding of ionising radiation should be less effective in the Galactic regions investigated by Jenkins (2013) than in the DLAs of our sample. As a result, for a given flux of ionising photons, argon should be more ionised in the ISM than in DLAs. Since this is not the case, the flux of hard photons must be more intense in DLA systems than in the Galactic neutral gas. The comoving specific emissivity of IGM ionising photons at the typical redshift of our sample () is at least one order of magnitude higher than at present time (see e.g. left panel of Fig. 15 in Haardt & Madau, 2012). The higher background at high redshift offers a natural explanation for the higher level of argon ionisation found in DLAs. This explanation is consistent with the hypothesis that the ionising source of argon in DLAs is external.

6.3 Internal sources of ionisation

Finding sources of Ar i ionisation internal to the DLA host galaxies faces the difficulty that the spectrum of ionising photons must include a significant hard component in order to explain the general level of argon deficiency (§5.3). Stellar sources, characterized by relatively soft spectra, may contribute only to some extent to the observed level of argon ionisation. In principle, bubbles of hot interstellar gas () created by supernova explosions could emit hard photons in the extreme UV/soft X. In fact, ionising sources of this type have been invoked by Jenkins (2013) to explain the argon deficiency observed in warm neutral regions of the Galactic ISM. In the case of the Galactic ISM, internal sources must be invoked given the weakness of the quasar background at . In our case this requirement is less compelling since the observed Ar deficiency can be reproduced assuming that the gas is ionised by the extragalactic radiation field at the redshift of the absorbers (§5.3).

Part of the DLA systems of our sample has a measurement of the C ii column density, which can be used to estimate the [C ii] 158 m cooling rate of the gas (Wolfe et al., 2008), N(C ii)/N(H i) , where N(C ii) is the column density of the excited state in the term of C, is the Einstein coefficient for spontaneous photon decay to the ground state ( s), and is the energy of this transition ( K). According to Wolfe et al. (2008), DLAs show a bimodal distribution of values, with two peaks at and , representative of “low-cool” DLAs and “high-cool” DLAs, respectively. In Fig. 9 we plot the [Ar/S] and [Ar/Si] ratios versus for our sample. There are only 7 systems where C ii has been measured and most of the Ar i measurements (5 out 7) belong to the “low” group. Even if the statistics at our disposal is quite small, one can see that in both DLA populations there are systems with both high and low Ar depletion: no trend is found between between Ar depletion and . Dutta et al. (2014), by means of CLOUDY models of metal-poor DLAs with measured C ii, and assuming a simple relation between SFR and CR ionisation rate, argued that in the “low-cool” population the radiation field is similar to that of the Milky Way, while in “high cool” DLAs the radiation field is more than 4 times higher. The lack of a trend between Ar depletion and observed in our sample suggests that in situ SFR should not play a relevant role in determining the different degrees of Ar photoionisation.

Figure 9: The [Ar/] ratios in DLA systems plotted versus cooling rate 6.3). The squares and triangles correspond to the [Ar/S] and [Ar/Si] ratios, respectively. The open squares correspond to proximate DLAs.

6.4 Interpretation of the trend with metallicity

As the metallicity increases, cooling processes become more efficient, so that the gas should become less ionised and neutral argon more abundant. In line with this theoretical expectation, our CLOUDY model with HM extragalactic radiation (§5.3) predicts the [Ar/S] and [Ar/Si] ratios to increase by 0.04 dex and 0.06 dex, respectively, when the metallicity increases from 0.01 to 1.0 the solar reference level (case , log N(H i) dex). Even if the effect is almost undetectable, it is clear that the model does not predict the observed decrease of [Ar/] with metallicity (Fig. 4, Tables 6-7). This suggests that the idealized model of a gas layer directly exposed to the extragalactic background needs to be refined. The neutral gas layer intercepted by the line of sight will be, in general, surrounded by gas of the host galaxy that may absorb part of the extragalactic radiation. The trend with metallicity could be due to the different degree of transparency of the gas that surrounds the observed line of sight. The trend would be indirectly induced by the links that are known to exist between star formation efficiency and metallicity. The DLA galaxies of our sample span a relatively small interval of lookback time (top axis of Fig. 6), so that a higher metallicity implies a higher star formation efficiency. In turn, a higher star formation efficiency implies a stronger removal of neutral gas (astration) and a stronger production of ionised gas (supernova remnants). Since ionised gas is transparent to hard photons while neutral gas is not, DLAs in high-metallicity environments will be more directly exposed to the extragalactic radiation than DLAs in low metallicity environments. This scenario would yield a stronger Ar i ionisation with increasing metallicity, in line with the observed trend. The trend may also be reinforced by local production of ionising photons from the hot gas heated by supernovae, as suggested by Jenkins (2013) for the Galactic ISM.

Bias effects

An additional possibility to explain the trend with metallicity is the well-known anti-correlation between metallicity and H i column density (Boissé et al., 1998; Vladilo & Péroux, 2005). As a result of this trend, DLA systems with high metallicity would experience more argon ionisation because they would have, on average, low H i self-shielding. In Fig. 8 one can see that also in our sample the highest values of metallicity [S/H] are found at low H i column density, whereas systems of high H i column density typically have low metallicities. Clearly, this trend can contribute to the trend of argon deficiency with metallicity that we observe. The existence of this bias makes less compelling, but does not rule out, the interpretation based on the transparency of the DLA galaxies to ionising photons.

7 Summary and conclusions

We have used the EUADP quasar database (Zafar et al., 2013b, a) to search for Ar i lines at the absorption redshift of DLAs/sub-DLAs identified in the quasar spectra. We report 3 new measurements and 5 upper limits of Ar i column density from the EUADP data. We combine these data with the literature high-resolution measurements published to date. This combined sample comprises of 37 systems, i.e. the largest argon sample collected so far. We have used this sample to investigate the relative abundance of argon in DLAs. The main results that we obtain can be summarized as follows.

  • Argon has a mean absolute abundance  dex (standard error of the mean), lower than the mean metallicity of DLAs, [M/H]  dex. When compared with elements of similar nucleosynthetic origin, such as the -capture elements silicon and sulphur, we find a mean deficiency  dex and  dex. The scatter of these [Ar/] ratios is somewhat higher than the scatter of other abundance ratios measured in DLAs, such as the [Zn/S] and [Si/Fe] ratios.

  • Both [Ar/S] and [Ar/Si] ratios show evidence for a positive correlation with H i column density, N(H i), and absorption redshift, , and for a negative correlation with the dust-free metallicity, [S/H]. These trends are characterized by modest correlation coefficients, the trend with redshift being particularly weak (Tables 6-7).

  • Detailed analysis of the measured abundance ratios indicates that the argon deficiencies are due to ionisation processes, rather than dust depletion or nucleosynthesis. The positive correlation between [Ar/] and log N(H i) is in line with the expectation that argon ionisation becomes less important with increasing self-shielding by neutral hydrogen. The level of argon deficiency is broadly in agreement with the one predicted by simple models of DLA ionisation by extragalactic field at the redshift of the absorbers.

  • The behaviour of the Ar/ ratios with redshift and N(H i) can be interpreted as the combination of two factors: a weakly evolving ionising background, with origin external to the DLA host galaxies, and an attenuation of the external radiation field according to H i self-shielding of individual absorbers. The combination of these two factors increases the scatter of each one of the two correlations with N(H i) and .

  • The argon abundances of the few proximate DLAs present in our sample lie at the lower boundary of the abundances measured in classical DLAs at the same H i column density or same absorption redshift. This result is consistent with the proposed scenario of argon ionisation by external quasar photons, with enhanced ionisation of proximate DLAs due to the vicinity of the background quasar.

  • The average [Ar/] abundance in DLAs is  dex lower than that found in warm interstellar regions of the Milky Way Jenkins (2013). The fact that the argon deficiency is stronger in DLAs than in the Galactic ISM is consistent with the strong rise of the metagalactic ionising background known to take place from to , the typical redshift of our sample.

  • The negative correlation between [Ar/] ratios and metallicity that we find (Fig. 4) cannot be explained with a simple model of gas layer directly exposed to the extragalactic background. We argue that DLAs in high-metallicity environments are likely to be more directly exposed to the extragalactic radiation than DLAs in low metallicity environments. Local production of ionising photons in hot gas might enhance argon ionisation in metal-rich environments.

  • Measurements of neutral argon in DLAs provide an independent test of the evolution of ionising sources in the Universe. Our measurements and upper limits of argon abundances suggest that the He ii reionisation takes place at redshift higher than . The possibility of setting firmer constraints on the redshift evolution of the He ii reionisation with our method requires a larger database of argon measurements at .


We would like to thank the ESO staff for making the UVES Advanced Data Products available to the community. We are thankful to the anonymous referee for his/her constructive comments. CP would like to thanks the BINGO! (‘history of Baryons: INtergalactic medium/Galaxies cO-evolution’) project by the Agence Nationale de la Recherche (ANR) under the allocation ANR-08-BLAN-0316-01. We are thankful P. Bonifacio for helpful discussions. G.V., M.C., and V.D. acknowledge support from the PRIN INAF “The 1 Billion Year Old Universe: Probing Primordial Galaxies and the Intergalactic Medium at the Edge of Reionisation”. P. M. acknowledges the international team #272 lead by C. M. Coppola “EUROPA - Early Universe: Research on Plasma Astrochemistry” at ISSI (International Space Science Institute) in Bern.

QSO log N(H i) log N(Ar i) log N(S ii) log N(Si ii) log N(Fe ii) Ref. [Ar/S] [Ar/Si] [S/H] [Zn/Fe] log
cm cm cm cm cm
Q 0000-2620 4.111 3.3901 1
QXO 0001 3.0000 2
LBQS 0010-0012 2.145 2.0250 3
B 0027-1836 2.550 2.4020 0.02 4
B 0100+1300 2.686 2.3090 5
J 0140-0839 3.716 3.6940 6
J 0142+0023 3.370 3.3477 6
Q 0201+365 3.610 2.4630 2
B 0336-017 3.197 3.0620 2, 7
B 0347-383 3.222 3.0250 7, 8, 9
B 0450-1310B 2.250 2.0667 10
B 0528-2505 2.756 2.8112 9, 11, 12
B 0642-5038 3.090 2.6590 3, 9
HS 0741+4741 3.221 3.0170 2, 7
B 0841+129 2.495 2.3745 10
B 0841+129 2.495 2.4760 10
Q 0930+2858 3.425 3.2350 2
B 1036-2257 3.103 2.7778 3
J 1113-1533 3.370 3.2665 3
J 1131+6044 2.921 2.8754 6
J 1211+0422 2.541 2.3766 13
B 1232+0815 2.570 2.3377 9, 14
J 1337+3152 3.174 3.1735 15
J 1340+1106 2.914 2.7958 16
BRI 1346-03 3.992 3.7360 7
J 1356-1101 3.006 2.9669 3
B 1409+0930 2.838 2.6681 3
Q 1425+6039 3.192 2.8268 2
J 1443+2724 4.420 4.2240 7, 9, 17
J 1604+3951 3.130 3.1670 6
GB 1759+7539 3.050 2.6260 18
HE 2243-6031 3.010 2.3310 19
J 2321+1421 2.538 2.5731 6
B 2332-094 3.330 3.0572 3
B 2342+3417 3.010 2.9092 2, 3
B 2343+125 2.763 2.4310 4, 9
Q 2344+1228 2.515 2.5380 7

References: 1: Molaro et al. (2001); 2: Prochaska et al. (2002a); 3: This work; 4: Noterdaeme et al. (2007); 5: Dessauges-Zavadsky et al. (2004); 6: Ellison et al. (2010); 7: Prochaska et al. (2001); 8: Levshakov et al. (2002); 9: Noterdaeme et al. (2008); 10: Dessauges-Zavadsky et al. (2006); 11: Centurión et al. (2003); 12: Vladilo et al. (2003); 13: Lehner et al. (2008); 14: Balashev et al. (2011); 15 : Srianand et al. (2010); 16: Cooke et al. (2011), Cooke (2014, priv. comm.); 17: Ledoux et al. (2006); 18: Prochaska et al. (2002b); 19: Lopez et al. (2002)
Sum of voigt profile components seen in Ar i. The total N(H i) is scaled down to match the fraction of S ii or Si ii measured in the Ar i voigt profile components.
Emission redshift, , of the quasar is not reported by Prochaska et al. (2002a).
Cases with measurements of Al iii/Al ii ratio: B 0450-1310 (log Al iii/Al ii  dex); HS 0741+4741 (log Al iii/Al ii  dex); J 1340+1106 (log Al iii/Al ii  dex).

Table 1: Metal column densities and abundance ratios of DLA/sub-DLAs seen in Ar i. The new estimates presented in this work are shown in bold-face. The proximate DLA cases are shown in italics.


  1. pagerange: The ESO UVES Advanced Data Products Quasar Sample - IV. On the deficiency of Argon in DLA systems1
  2. pubyear: 2013
  3. This happens because the Al ii line 1676Å is often saturated in DLAs. The few measurements of the Al iii/Al ii ratio available only in our sample are indicated with the footnote “c” in Table 1.
  4. The ionisation parameter is the dimensionless ratio , where is the total surface flux of ionising photons (cm s) and the H particle density (cm).
  5. Photons in this energy range can penetrate the H i layers because the H i photoionisation cross section scales as and therefore neutral layers tend to become transparent to ionising photons when  eV.
  6. The DLAs at towards Q 0027-1836 and at towards Q 0100-1300, with log N(H i) and  dex, respectively.
  7. The ISM argon deficiency quoted by Jenkins (2013) is [Ar i/O i]  dex; this author adopts solar reference values and .
  8. We estimate log N(H i) for the Galactic lines of sight from the equivalent widths of the O i lines published by Jenkins (2013), assuming solar abundances. The Ar i absorption lines become saturated in a gas of solar composition when log N(H i) . Therefore it is practically impossible to compare Ar i measurements in ISM and DLA sightlines with same H i column density.


  1. Asplund, M., Grevesse, N., Sauval, A. J., & Scott, P. 2009, ARAA, 47, 481
  2. Balashev, S. A., Petitjean, P., Ivanchik, A. V., et al. 2011, MNRAS, 418, 357
  3. Boissé, P., Le Brun, V., Bergeron, J., & Deharveng, J.-M. 1998, A&A, 333, 841
  4. Bolton, J. S., Becker, G. D., Raskutti, S., et al. 2012, MNRAS, 419, 2880
  5. Centurión, M., Molaro, P., Vladilo, G., et al. 2003, A&A, 403, 55
  6. Cooke, R., Pettini, M., Steidel, C. C., Rudie, G. C., & Nissen, P. E. 2011, MNRAS, 417, 1534
  7. Dall’Aglio, A., Wisotzki, L., & Worseck, G. 2008, A&A, 491, 465
  8. Dekker, H., D’Odorico, S., Kaufer, A., Delabre, B., & Kotzlowski, H. 2000, in SPIE Conference Series, Vol. 4008, Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, ed. M. Iye & A. F. Moorwood, 534–545
  9. Dessauges-Zavadsky, M., Calura, F., Prochaska, J. X., D’Odorico, S., & Matteucci, F. 2004, A&A, 416, 79
  10. Dessauges-Zavadsky, M., Prochaska, J. X., D’Odorico, S., Calura, F., & Matteucci, F. 2006, A&A, 445, 93
  11. Duley, W. W. 1980, MNRAS, 190, 683
  12. Dutta, R., Srianand, R., Rahmani, H., et al. 2014, MNRAS, 440, 307
  13. Ellison, S. L. & Lopez, S. 2009, MNRAS, 397, 467
  14. Ellison, S. L., Prochaska, J. X., Hennawi, J., et al. 2010, MNRAS, 406, 1435
  15. Fan, X., Strauss, M. A., Richards, G. T., et al. 2006, AJ, 131, 1203
  16. Ferland, G. J., Korista, K. T., Verner, D. A., et al. 1998, PASP, 110, 761
  17. Ferland, G. J., Porter, R. L., van Hoof, P. A. M., et al. 2013, RMxAA, 49, 137
  18. Fontana, A. & Ballester, P. 1995, The Messenger, 80, 37
  19. Furlanetto, S. R. & Oh, S. P. 2008, ApJ, 681, 1
  20. Grevesse, N. & Sauval, A. J. 1998, SSR, 85, 161
  21. Haardt, F. & Madau, P. 2012, ApJ, 746, 125
  22. Henry, R. B. C. & Worthey, G. 1999, PASP, 111, 919
  23. Izotov, Y. I. & Thuan, T. X. 1999, ApJ, 511, 639
  24. Jarosik, N., Bennett, C. L., Dunkley, J., et al. 2011, ApJS, 192, 14
  25. Jenkins, E. B. 2013, ApJ, 764, 25
  26. Jenkins, E. B., Oegerle, W. R., Gry, C., et al. 2000, ApJL, 538, L81
  27. Kurucz, R. L. 1991, in Precision Photometry: Astrophysics of the Galaxy, ed. A. G. D. Philip, A. R. Upgren, & K. A. Janes, 27
  28. Lanz, T., Cunha, K., Holtzman, J., & Hubeny, I. 2008, ApJ, 678, 1342
  29. Ledoux, C., Bergeron, J., & Petitjean, P. 2002, A&A, 385, 802
  30. Ledoux, C., Petitjean, P., & Srianand, R. 2006, ApJL, 640, L25
  31. Lehner, N., Howk, J. C., Prochaska, J. X., & Wolfe, A. M. 2008, MNRAS, 390, 2
  32. Lehner, N., Jenkins, E. B., Gry, C., et al. 2003, ApJ, 595, 858
  33. Levshakov, S. A., Dessauges-Zavadsky, M., D’Odorico, S., & Molaro, P. 2002, ApJ, 565, 696
  34. Lodders, K. 2003, ApJ, 591, 1220
  35. Lopez, S., Reimers, D., D’Odorico, S., & Prochaska, J. X. 2002, A&A, 385, 778
  36. McQuinn, M., Lidz, A., Zaldarriaga, M., et al. 2009, ApJ, 694, 842
  37. Møller, P., Warren, S. J., & Fynbo, J. U. 1998, A&A, 330, 19
  38. Molaro, P., Levshakov, S. A., D’Odorico, S., Bonifacio, P., & Centurión, M. 2001, ApJ, 549, 90
  39. Morton, D. C. 2003, ApJS, 149, 205
  40. Noterdaeme, P., Ledoux, C., Petitjean, P., et al. 2007, A&A, 474, 393
  41. Noterdaeme, P., Ledoux, C., Petitjean, P., & Srianand, R. 2008, A&A, 481, 327
  42. Petitjean, P., Ledoux, C., & Srianand, R. 2008, A&A, 480, 349
  43. Pettini, M., Smith, L. J., Hunstead, R. W., & King, D. L. 1994, ApJ, 426, 79
  44. Pettini, M., Smith, L. J., King, D. L., & Hunstead, R. W. 1997, ApJ, 486, 665
  45. Prochaska, J. X., Gawiser, E., Wolfe, A. M., Cooke, J., & Gelino, D. 2003, ApJS, 147, 227
  46. Prochaska, J. X., Henry, R. B. C., O’Meara, J. M., et al. 2002a, PASP, 114, 933
  47. Prochaska, J. X., Howk, J. C., O’Meara, J. M., et al. 2002b, ApJ, 571, 693
  48. Prochaska, J. X., Wolfe, A. M., Tytler, D., et al. 2001, ApJS, 137, 21
  49. Rafelski, M., Wolfe, A. M., Prochaska, J. X., Neeleman, M., & Mendez, A. J. 2012, ApJ, 755, 89
  50. Savage, B. D. & Sembach, K. R. 1996, ARAA, 34, 279
  51. Sofia, U. J. & Jenkins, E. B. 1998, ApJ, 499, 951
  52. Srianand, R., Gupta, N., Petitjean, P., Noterdaeme, P., & Ledoux, C. 2010, MNRAS, 405, 1888
  53. Srianand, R., Gupta, N., Petitjean, P., et al. 2012, MNRAS, 421, 651
  54. Syphers, D., Anderson, S. F., Zheng, W., et al. 2009, ApJS, 185, 20
  55. Vladilo, G., Abate, C., Yin, J., Cescutti, G., & Matteucci, F. 2011, A&A, 530, A33
  56. Vladilo, G., Centurión, M., Bonifacio, P., & Howk, J. C. 2001, ApJ, 557, 1007
  57. Vladilo, G., Centurión, M., D’Odorico, V., & Péroux, C. 2003, A&A, 402, 487
  58. Vladilo, G. & Péroux, C. 2005, A&A, 444, 461
  59. Wolfe, A. M., Gawiser, E., & Prochaska, J. X. 2005, ARAA, 43, 861
  60. Wolfe, A. M., Prochaska, J. X., Jorgenson, R. A., & Rafelski, M. 2008, ApJ, 681, 881
  61. Wolfe, A. M., Turnshek, D. A., Smith, H. E., & Cohen, R. D. 1986, ApJS, 61, 249
  62. Zafar, T., Møller, P., Ledoux, C., et al. 2011, A&A, 532, A51
  63. Zafar, T., Péroux, C., Popping, A., et al. 2013a, A&A, 556, A141
  64. Zafar, T., Popping, A., & Péroux, C. 2013b, A&A, 556, A140
Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
Add comment
Loading ...
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test description