V V Rahmani

C Klimenko

P Ledoux

R Petitjean

A V Srianand

D A Ivanchik

Varshalovich

S A Balashev
email: s.balashev@gmail.com

P Noterdaeme

H Rahmani

V V Klimenko

C Ledoux

P Petitjean

R Srianand

A V Ivanchik

D A Varshalovich

CO-dark molecular gas at high redshift: very large H 2 content and high pressure in a low metallicity damped Lyman-alpha system

Keywords: cosmology, observations -ISM, clouds -quasar, absorption lines

CO-dark molecular gas at high redshift: very large H-2 content and high pressure in a low-metallicity damped

Lyman alpha system

INTRODUCTION

The processes of star formation, galaxy evolution and AGN activities are directly connected with the evolution of different phases of the neutral interstellar medium (ISM). Recent kpcscale observations of nearby normal galaxies (non-starbursts) have shown that the surface star formation rate (Σ SFR ) in their inner parts correlates approximately linearly with the surface density of molecular hydrogen (Σ H 2 ), but not with that of atomic gas (see Bigiel et al. 2008;Leroy et al. 2013). In other words, the star formation in inner parts of nearby galaxies is closely associated with the molecular phase of the ISM. While such a tight correlation was found drogen is an essential ingredient for star formation or only a tracer of the cold and dense ISM (Krumholz et al. 2011;Krumholz 2012;Glover & Clark 2012a;Michalowski et al. 2015;Glover & Smith 2016). The basis for this question is that molecular hydrogen is not major coolant of the cold ISM even at low metallicities. However, the cooling of the gas through atomic emission lines depends on the metal abundance, to which the formation of H 2 is also linked, via the correlation between gas-phase and dust-phase metal abundances. In addition, H 2 is vital for other complex molecules to form, which can cool dense regions of the ISM and eventually initiate the collapse of the cloud, leading to the formation of stars.

It is not yet established whether the relation between star formation and H 2 surface/column density is valid in high-redshift galaxies, for which the metallicities are expected to be much lower on average (e.g. Rafelski et al. 2012). For example, recent works suggest that at low metallicities star formation can occur before the H i converts into H 2 (Krumholz et al. 2011;Hu et al. 2015). In addition, Σ SFR seems to remain well correlated with H i in dwarf irregular galaxies with metallicities ∼ 0.1 × Z (Roychowdhury et al. 2014). The analysis of damped Ly-α systems (DLAs) detected in the spectra of quasar and gamma-ray bursts (GRB) afterglows provides a unique method to probe the atomic and molecular hydrogen at high redshifts (e.g. Noterdaeme et al. 2015b). However, the relation with star-formation remains difficult to establish due to the paucity of direct detections of the associated galaxies, although global scaling relations (between e.g. SFR, metallicity and impact parameter) start to be revealed (Krogager et al. 2017). The measured SFR in the overall population of intervening DLAs appears to be very low (Fumagalli et al. 2015;Rafelski et al. 2016), while DLAs associated with GRB afterglows span a range of SFRs (Toy et al. 2016).

Finally, there is no unambiguous way to convert measurement of column densities towards one particular line of sight to the overall properties of the galaxy. In spite of that, it is known that H 2 is less abundant in DLAs at high redshifts (such DLAs typically show lower metallicity), i.e. in the diffuse ISM of high-redshift galaxies: only a small fraction of high-z DLAs (< 10 per cent; Noterdaeme et al. 2008;Balashev et al. 2014;Jorgenson et al. 2014;Noterdaeme et al. 2015a) show associated H 2 . This contrasts with observations in the local Universe, where the fraction of H i absorbers with associated H 2 is ∼ 90% for our Galaxy (Wakker 2005) and ∼ 70% for Magellanic Clouds (Tumlinson et al. 2002) with molecular fraction f > 10 -6 .

It is hard to attribute the difference in the H 2 detection rate between high redshifts and the local Universe to a pure observational effect (obscuration of the background quasar by dust associated with H 2 -bearing clouds), since it was found that in the majority of DLAs at high redshifts (z > 1.8), the dust extinction is negligible (e.g Murphy & Bernet 2016), although some highly reddened cases do exist (e.g. Srianand et al. 2008a;Noterdaeme et al. 2009;Ma et al. 2015;Krogager et al. 2016;Noterdaeme et al. 2017), in which cold gas is found. Therefore the low incidence rate of H 2 is most probably connected with the small covering factor of H 2bearing cold ISM at high redshifts.

The small incidence rate of H 2 absorption systems in high redshift DLAs also implies that only about 30 H 2 absorption systems at high redshifts have been identified and therefore the cold ISM remains least studied. Blind searches for H 2 absorption systems are quite ineffective because the detection of molecular absorption usually requires high-resolution (or intermediate in the case of very saturated H 2 absorption system) and high-S/N spectrum to overcome the effect of contamination by the Lyα forest. In addition, among the known high-z H 2 -bearing DLAs, only few have high H 2 column density, log N 19 1 , and high molecular fraction, f 0.1.

Studies of low-metallicity H 2 -bearing systems can be used to estimate the fraction of CO-dark molecular gas, which is not seen in CO emission (Wolfire et al. 2010). A recent work indicates that such gas is important for star formation, since it is more susceptible to turbulent compression and gravitational collapse than pure atomic gas (Glover & Smith 2016). It was measured that a large fraction of the molecular gas can be in CO-dark phase: observations of C ii in the Milky Way (MW) show that, on an average, the CO-dark gas can contain as much as ∼30% of the total H 2 content (Pineda et al. 2013). This fraction is expected to increase with decreasing metallicity (Glover & Clark 2012b). Indeed, at low metallicity the relative abundance of CO to H 2 is reduced by the reduced abundances of C and O together with a weaker dust extinction. Observationally, neither the fraction of CO-dark gas nor X CO , the conversion factor between CO emission and H 2 mass (see the recent review by Bolatto et al. 2013) is well constrained in galaxies both at high redshifts and with low metallicity (Shi et al. 2016, and references therein). Any measurements of these parameters are very valuable, since with the advent of state-of-the-art sub-mm facilities like ALMA, CO can now be detected in normal galaxies at high redshifts.

Several efficient techniques to search for molecular-rich absorption systems at high redshift have recently been proposed. These techniques use different target pre-selections in the Sloan Digital Sky Survey (SDSS; York 2000) spectroscopic data base. One technique based on the presence of strong neutral carbon lines in SDSS quasar spectra (see Ledoux et al. 2015) was proposed to select high-metallicity CO/H 2 -bearing systems. This led to the first studies of the translucent ISM at high redshifts (e.g. Srianand et al. 2008b;Noterdaeme et al. 2009Noterdaeme et al. , 2010)). It was also shown that the 2175Å bump can be detected in SDSS spectra (e.g. Wang et al. 2004;Srianand et al. 2008a;Jiang et al. 2011) and could potentially be used to preselect molecular-rich absorbers. Indeed the followup studies (e.g. Ma et al. 2015;Noterdaeme et al. 2017) showed that both C i and 2175Å bump selection techniques probe molecular bearing systems with high metallicities, up to super-solar. These are at the high end of the metallicity distribution of high-redshift DLAs and could therefore probe gas with chemical enrichment closer to what is seen in the local Universe than in typical high-z galaxies.

Another technique was proposed by Balashev et al. (2014) to directly search for H 2 absorption systems at z > 2.3 in SDSS quasar spectra. In short, the technique uses a template matching routine to automatically detect damped Lyman and Werner H 2 lines using the SDSS DLA catalogue (Noterdaeme et al. 2012c). A confidence is attributed to each candidate based on a resampling method and using control samples. While only about 3 per cent of SDSS spectra are suitable to search for H 2 , and only very high column density H 2 absorption systems (with log N 19.5) can be detected, it is worth mentioning that this technique does not depend on metal tracers, and should therefore probe high N(H 2 ) over a wide range of metallicities.

To date, we have detected ∼ 100 reliable H 2 candidates using SDSS DR12 and started follow-up observations with the Very Large Telescope (VLT). All of the eight candidates followed-up so far have been confirmed as strong H 2 systems.

In this paper, we present the study of the system with strongest H 2 absorption in our sample, at z abs = 2.786 towards the quasar SDSS J 084312.72+022117.28 (hereafter J 0843+0221) and postpone the study of the overall sample to a future work. H 2 lines towards J 0843+0221 are so saturated that they were initially identified as multiple DLA systems in the automatic DLA search by Noterdaeme et al. (2012c). We obtained high-resolution spectroscopic observations of the system using the Ultraviolet and Visual Echelle Spectrograph (UVES) on the VLT. Section 2 presents these observations and the data reduction. Section 3 presents the analysis of the system, including metal lines, associated molecules and dust extinction. We estimate the physical conditions in the gas in Section 4 and discuss our findings in Section 5. We conclude in Section 6.

OBSERVATIONS AND DATA REDUCTION

The observations of J 0843+0221 (z em =2.92) were carried out between 2014 January and March using UVES [START_REF] Dekker | Proc. SPIE, Optical and IR Telescope Instrumentation and Detectors[END_REF] under programme 092.A.0345 (PI: D.A. Varshalovich). The journal of the observations is presented in Table 1. We used two standard settings: 390+564 (with dichroic 1 and cross-dispersers CD#2 and CD#3) and 437+760 (with dichroic 2 and cross-dispersers CD#2 and CD#4) to cover most of the metal and molecular absorption lines associated with the DLA at z = 2.7865. A slit width of 1.0 arcsec and CCD readout with 2×2 binning were used for all the observations, resulting in the spectral resolution R ≈ 48 000.

We reduced the data using UVES Common Pipeline Library (CPL) data reduction pipeline release 6.5.1 2 with the optimal extraction method. To avoid extra rebinning of the pixels we used the final un-rebinned extracted spectrum of each order produced by CPL and binned it into a global grid of uniform pixel size of 2.5 km s -1 . We applied the CPL wavelength solutions to each order and merged the orders by implementing a weighted mean in the overlapping regions. All the exposures were shifted to the heliocentricvacuum frame correcting for the motion of the observatory along the quasar line of sight obtained at the exposure mid-point and using the air-to-vacuum relation from Edlén (1966). Then we combined all exposures using a weighted mean at each pixel. The standard deviation around the continuum level in our final combined spectrum is in agreement, within a few percent, with our error spectrum. The final spectrum covers the observed wavelength range of 328-947 nm with gaps of ∼ 7 and 11 nm centred at 564 and 760 nm, respectively.

ANALYSIS

We analysed absorption lines from the DLA using a code that simultaneously fits the absorption profiles of selected species using Voigt profile convolved with the instrumental response function. We used χ 2 likelihood function, which assumes the probability distribution function of uncertainties in the spectrum to be Gaussian. The shape of the likelihood function is defined by the Monte Carlo Markov chain approach with implementation of the affine-invariant ensemble sampler (Goodman & Weare 2010). Such a technique ensures that we confidently find the global maximum in the many parametric space and provides reliable estimates of statistical errors on the parameters. Furthermore, the correlation between different parameters can be easily traced by this method. We crosschecked our results for several species with an independent analy-2 http://www.eso.org/sci/facilities/paranal/instruments/uves/doc/ sis obtained from vpfit3 and found a good match between the two approaches. All the errors quoted here are given using 63.8% quantile interval, which formally corresponds to 1σ interval for normal distribution 4 . In each figure presenting a fit to absorption lines (except Figs 1 and Fig. 3), we grey out pixels of spectrum not used to constrain the fit and show the residuals (data minus model, divided by uncertainty).

Neutral hydrogen

We measured the H i column density by fitting the damped Ly-α absorption line with a Voigt profile. Higher Lyman series of the H i absorption lines are severely blended with the damped absorptions produced by associated H 2 . Because the Ly-α line is strongly damped and located near the Ly-α emission line of J 0843+0221 (δv ∼ 10000 km s -1 ), we simultaneously fitted the continuum (using first seven Chebyshev polynomials) and the DLA profile. We obtained log N(H i) ∼ 21.82 ± 0.11 at the best-fitted redshift, z DLA =2.78654 ± 0.00015 (see Fig. 1), confirming the value previously derived from SDSS data (log N(H i) = 21.8; Noterdaeme et al. 2014).

Neutral chlorine as a proxy to the velocity structure of H 2

The strong saturation of H 2 absorption lines prevents us from deriving the velocity structure of H 2 in its first five rotational levels. Several components can be distinguished in the J ≥ 5 levels but the velocity structure cannot be extracted unambiguously because of the relatively low S/N of the spectrum. Therefore we assume the velocity structure of the H 2 to be the same as that of Cl i, which is known to be a good tracer of H 2 -bearing gas (Jura 1974), both statistically in our Galaxy (Moomey et al. 2012) (Balashev et al. 2015), but also within a given absorbing cloud (Noterdaeme et al. 2017).

Absorption lines from two neutral chlorine transitions (Cl i λ1347 and 1088) are detected in two main components at z=2.786582 and 2.786459 (see Fig. 2). A third, weaker component can be distinguished in the red wing of Cl i λ1347, which is also seen in the line profiles of H 2 with high-J and low-ionization metals. In the following we will denote these components as A, B and C in order of increasing redshift. The Cl i λ1088 line is likely to be blended with unrelated absorption from the Ly-α forest since its apparent blue wing (v < -5 km s -1 ) has a higher optical depth than the corresponding region of Cl i λ1347, despite the oscillator strength of the former being smaller ( f =0.081) than that of the latter ( f =0.153; Schectman et al. 1993). Therefore we excluded Cl i λ1088 from the fit. The results of the fit are given in Table A1.

Molecular hydrogen

Molecular hydrogen is detected in the first 10 rotational levels of the ground (ν = 0) vibrational state. This is the first detection of such a large number of excited levels from lower vibrational state of H 2 in an intervening high-z DLA on QSO line of sight (there are two detections in DLA associated with GRB afterglows; Sheffer et al. 2009;Krühler et al. 2013).

In all detected Lyman (L) and Werner (W) bands of H 2 , the R(0), R(1) and P(1) lines are strongly damped and significantly blended with each other (Fig. 3). This affects the apparent quasar continuum. We therefore fitted simultaneously the damped absorption lines and the continuum using the first 10 Chebyshev polynomials. This is enough to take into account possible continuum variations due to the blaze function of our echelle spectrum, which is not flux calibrated. While we have no information on the column densities in the first two H 2 rotational levels for each component individually, the total H 2 column density is obtained with high accuracy thanks to the damping wings: we find log N(H 2 , J=0) = 20.71 ± 0.02 and log N(H 2 , J=1) = 21.06 ± 0.02 at the best-fitted redshift of z=2.786570(14). The best fit to the J = 0 and J = 1 H 2 absorption lines is shown in Fig. 3. We also checked that a fit using three components with velocity structure fixed to that of Cl i provides very similar total column densities of H 2 .

The measured excitation temperature T 01 = 123 +9 -8 K, is higher than what is usually obtained at high H 2 column densities (log N(H 2 ) > 19) both in the local Universe and at high redshifts. We will discuss this further in Sect. 4.2. The fit to J ≥ 2 rotational levels of H 2 is presented in Appendix A and its results are summarized in Table A1.

HD molecule

Absorption lines from deuterated molecular hydrogen (HD) are detected in components A and B in the first two rotational levels. This is only the second detection of HD, J=1 lines at high z (see Balashev et al. 2010, for the first detection towards Q 0812+3208). The Figure 3. Portion of our UVES spectrum of J 0843+0221 around the region of H 2 absorption lines from the DLA at z = 2.786. The red solid line shows the total fitted profile to the J = 0 and J = 1 absorption lines, while green and orange dashed lines represent the contribution of the J = 0 and J = 1 levels, respectively. This is the largest H 2 column density DLA detected in QSO spectra.

non-detection of HD in the component C is consistent with low H 2 column density in this component. Since HD components are severely blended with each other, we again used fixed redshifts of the components as determined from Cl i. For the typical physical conditions in the cold ISM, the excitation of the first rotational levels of HD from the ground level is expected to be dominated by collisions (see Abgrall et al. 1992;Balashev et al. 2009). Therefore we assumed that HD at J = 0 and J = 1 levels are co-spatial, and tied together their Doppler parameters. The fit to HD lines is shown in Fig. 4 and the parameters are given in Table A1. Unfortunately, only two lines from J = 0 rotational level and three lines from J = 1 level of HD can be used to constrain the fit. In addition, these lines are in the intermediate regime, leading to large uncertainties on the column densities. We obtain the total column density log N(HD) = 17.35 +0.15 -0.34 .

Neutral carbon

We detect C i lines corresponding to transitions from three finestructure levels of the 2s 2 2p 2 3P 0 C i ground state. As for HD, C i is only detected in components A and B. The fit to the C i lines is shown in Fig. 5 and the results are given in Table A1. The C ito-H 2 column density ratio is relatively low, (6.0 ± 0.4) × 10 -8 , compared to what has been seen so far in high-z H 2 -bearing systems (e.g. Srianand et al. 2005;Jorgenson et al. 2010). This can indicate an inefficient dust shielding of C i (i.e. higher UV flux and/or low dust content) and/or the low metallicity of the gas (see Sect. 3.7) that leads to a low abundance of carbon.

CO molecules

Carbon monoxide, which can only survive in the cold H 2 -bearing gas, has been detected at high redshift in only a few DLA systems. When measurable, these systems show high H 2 content (Srianand et al. 2008b;Noterdaeme et al. 2010Noterdaeme et al. , 2017)). Here, we do not detect CO absorption lines in spite of the extreme H 2 column density. Using several A-X CO bands as well as C-X(0 -0) and E-X(0 -0) CO bands we derive 3σ upper limits on the total CO column density assuming an excitation temperature 10.3 K, which is the cosmic microwave background (CMB) temperature at the absorption redshift z = 2.786 (see Noterdaeme et al. 2011). We derive log N(CO) < 13.1 and 13.0 for components A and B, respectively.

Metal content

Several metal species in different ionization stages (e.g. O i, Fe ii, Si iv) originating from the DLA are detected. While we concentrate here on the low-ionization metals, which provide the necessary information to derive the chemical abundances in the neutral gas, a comparison with highly ionized phase is also provided in Fig. A1, for completeness. The main components seen in the profiles of the low-ionization metal lines correspond to the three components identified in Cl i and H 2 . The redshifts of these components were left free during the fit and found to be in agreement (at 2σ level) with those obtained from the analysis of Cl i. For each component we tied the Doppler parameters of metal species using microturbulence assumption as b 2 = b 2 turb +2k B T/M, where M is the atomic mass of the species and k B is the Boltzmann constant. The turbulence Doppler parameter, b turb , and kinetic temperature of the gas, T , were left free during the fitting process. Our best fit is overlaid on the data in Fig. 6, and the corresponding parameters are given in Table A1. We note that a fit with independent b parameters for each metal species gives column densities that are in agreement within 1σ with the 'tied' fit, but have larger individual uncertainties. We provide here below some comments about individual species.

Zn ii

We estimate the column density of Zn ii from the Zn ii λ2026 and 2062 lines. Zn ii λ2026 is redshifted in a spectral region crowded with O 2 sky lines. Since the different exposures with 437 + 760 settings (that cover this spectral region) were taken at different times, the velocity of the observatory with respect to the line of sight is different, while O 2 lines remain stationary. As a result, we find that in three out of five exposures, the Zn ii λ2026 line is affected by sky

+2σ -2σ 0 1 Mg II 1240 A B C +2σ -2σ 0 1 Si II 1808 A B C +2σ -2σ 0 1 S II 1250 A B C +2σ -2σ 0 1 Ti II 1910 A B C +2σ -2σ 0 1 Ti II 1910 A B C +2σ -2σ 0 1 Normalized flux Cr II 2062 A B C +2σ -2σ 0 1 Cr II 2056 A B C +2σ -2σ 0 1 Mn II 1197 A B C +2σ -2σ 0 1 Fe II 2260 A B C +2σ -2σ -60 -30 0 30 60 v [km s -1 ] 0 1 Fe II 2249 A B C +2σ -2σ
Fe II 1611

A B C +2σ -2σ
Fe II 1142

A B C +2σ -2σ
Ni II 1751

A B C +2σ -2σ
Ni II 1709

A B C +2σ -2σ
Ni II 1741

A B C +2σ -2σ
Ni II 1703

A B C +2σ -2σ
Ni II 1454

A B C +2σ -2σ
Ni II 1370

A B C +2σ -2σ
Zn II 2062

A B C +2σ -2σ -60 -30 0 30 60 v [km s -1 ]
Zn II 2026 lines but not in the remaining two. We therefore used the combination of these two exposures to constrain the fit of Zn ii. Due to small velocity extent of the metal profile (∼ 50 km/s), there is no contamination of Zn ii λ2026 and 2062 by Mg i λ2026 and Cr ii λ2062, respectively.

Si ii

To determine the column densities of Si ii in the three main components we use the relatively weak Si ii λ1808 line. Other Si ii lines, Si ii λ1304, 1526 and 1260 lines are highly saturated and show at least four additional components. Using these lines (together with Si ii λ1808 line) we estimate Si ii column density in these additional components to be log N(Si ii) ∼ 15. This is significantly lower than what is measured in the main components log N(Si ii) ∼ 15.75. The situation is similar for other species such as Fe ii.

S ii

We estimate the column density of S ii using S ii λ1250 while two other transitions (S ii λ1253 and S ii λ1259) are partially blended with other lines. Interestingly, using the best-fitting parameters obtained from nearly saturated S ii λ1250 we found that the calculated profile of S ii λ1259, whose oscillator strength is three times larger than that of S ii λ1250, is deeper than the observed one. This is very likely the consequence of a partial coverage of the Ly-α emission line region of the quasar by the S ii absorbing region, since S ii λ1259 is redshifted on the top of the quasar Ly-α emission line (see Fig. 7). Partial coverage has already been detected in at least two H 2 -bearing intervening systems (e.g. Balashev et al. 2011;Klimenko et al. 2015). From the residual flux at the bottom of the S ii λ1259 line, we estimate that about ∼ 20 per cent of the quasar photons at the corresponding wavelengths are leaking through the S ii absorbing region of the DLA. Since Ly-α emission contributes to ∼ 80 per cent of emitted photons at that wavelength 5 , this means that about 25 per cent of the Ly-α photons arise from regions not covered by the DLA. Similar partial coverage (by H i) of Ly-α quasar emission has been found for proximate DLAs at the quasar redshift (see Finley et al. 2013), where further investigation show that the Ly-α emission can be very spatially extended (Fathivavsari et al. 2016;Borisova et al. 2016). Therefore such situation can be quite common.

P ii

Phosphorus is a particularly useful element for metal enrichment studies, since it is believed to be little dust depleted (De Cia et al. 2016). Unfortunately, the strongest absorption line of singly ionized phosphorus, P ii λ1152 is affected by blends. P ii λ1301 is partially blended and not saturated. Using a three-component fit with fixed redshift and Doppler parameter deduced from the overall metal fit, we obtained an upper limit on the total phosphorus column density of log N(P ii) < 13.72. The corresponding upper limit on the metallicity [P/H] < -1.28 is consistent with low depletion of the phosphorus in the DLA at z=2.7865.

Excited C ii* and Si ii*

The detection of excited Si ii* was recently reported in the two extremely strong DLAs (ESDLAs) at z = 2.21 towards J 1135-0010 (Kulkarni et al. 2012) and at z = 2.34 towards J 2140-0321 (Noterdaeme et al. 2015b). Neeleman et al. (2015) detected Si ii* in another seven DLAs, which also correspond to the high end of the H i column density distribution. The ratio N(Si ii * )/N(Si ii) can be used to estimate the number density (or electron density) since the excited state of Si ii is predominantly populated by collisions and depopulated by spontaneous transitions (Srianand & Petitjean 2000). Since the DLA studied here also has a very high column density (log N(H i) ≈ 21.8), we searched for Si ii* absorptions lines, which we tentatively detect in the strongest transition Si ii* λ1264 (see Fig. 8). Fixing the redshifts and Doppler parameters to that obtained for Si ii, our fitting procedure indicates log N(Si ii * ) = 11.32 +0.14 -0.21 and log N(Si ii * ) = 11.62 +0.08 -0.14 for components A and B, respectively. Because the detection is only tentative, we rather consider the 2σ upper limit on the total Si ii* column density (log N(Si ii * ) < 12) in the following. This translates to N(Si ii * )/N(Si ii) ≤ 1.7 × 10 -4 , i.e., similar to what is found towards J 2140-0321 but an order of magnitude lower than the ratio measured towards J 1135-0010.

We derive the column density of excited singly ionized carbon using the C ii* λ1335 line (see Fig. 8). Because this line is saturated, we used a three component fit with redshifts and Doppler parameters fixed to the values derived from the joint fit of other metal lines. We found that most of the C ii* column density arises from component B (see Table A1).

Extinction of the background quasar's light

We used the flux-calibrated SDSS spectrum of J 0843+0221 to measure the extinction due to the H 2 -bearing DLA at z = 2.7865. We followed a standard procedure similar to that presented in Srianand et al. (2008b): we fitted the observed spectrum of J 0843+0221 using the SDSS composite spectrum from SDSS Collaboration & Berk (2001), redhifted to z em = 2.92 and reddened using an extinction law applied at z abs = 2.7865. Since no 2175 Å bump is detected, we used the average Small Magellanic Cloud (SMC) extinction law (Gordon et al. 2003) and found a best-fitting value of A V = 0.07 (Fig. 9). In order to estimate the systematic uncertainty due to intrinsic variations of the spectral energy distribution of quasars, we used a control sample of quasar spectra from SDSS DR9 with no detected DLA and with emission redshifts in the range ∆z = 0.05 around that of J 0843+0221. We applied the same fitting procedure to each spectrum in the control sample and obtained the distribution of measured A V shown in the inlay of Fig. 9. The median and dispersion in the control sample are -0.02 and 0.10, respectively, which give the zero-point and systematic error for the measured A V in J 0843+0221. We thus obtain A V = 0.09 ± 0.10 towards J 0843+0221.

Such a low amount of extinction is incompatible with the value expected from extrapolating the canonical extinction-to-gas ratio seen in the MW: for the total hydrogen column density log N(H) = log N(H i) + log 2N(H 2 )) = 21.99 +0.08 -0.07 , we would expect A V ≈ 5 (Watson 2011). Using the SMC and Large Magellanic Cloud (LMC) extinction-to-gas ratio (De Cia et al. 2013) does not change significantly the expected value. The low A V measured here is however not surprising given the low metallicity of the DLA. Using the observed metallicity and depletion factor, we estimate the dust-to-gas ratio to be κ = [Zn/H] 1 -10 expected A V should be reduced by a factor of about 30 compared to the above extrapolations. We then expect A V = 5 × 0.03 ≈ 0.15, i.e. very close to the observed value.

PHYSICAL CONDITIONS

Metallicity and depletion pattern

For each element we determine the abundance relative to solar values from Asplund et al. (2009), without applying any ionization correction (i.e. we assume all metals to be in their main ionization state). This is a fairly reasonable assumption given the large H i column density in this DLA. Since we cannot determine the hydrogen column density in individual components, we de- Zn ii -1.52 +0.08 -0.10 rive gas-phase abundances using the total hydrogen column density N(H) = 21.99 +0.08 -0.07 and the total metal column densities, see Table 2. The observed gas-phase abundances of zinc and sulphur are usually considered as good indicators of the ISM metallicity since these are non-refractory elements for which dust depletion corrections are expected to be negligible. Their abundances, [Zn/H] = -1.52 +0.08 -0.10 and [S/H] = -1.54 +0.08 -0.08 , are indeed in very good agreement.

Other elements have lower observed gas-phase abundances. In particular, Fe ii, Cr ii, Ni ii and Ti ii are refractory elements that are known to severely deplete on to dust grains. In Fig. 10, we compare the observed depletion pattern in our DLA with typical values seen in the MW halo, the warm and the cold ISM (Welty et al. 1999). The observed depletion pattern is found to be intermediate between those observed for the warm and halo MW gas and to be much less depleted than the cold ISM gas in the MW. As noted by Noterdaeme et al. (2007), this is most likely a metallicity effect and does not mean that the DLA gas is either warm or in a halo. Indeed, the measured Fe ii depletion, [Fe/Zn] = -1.01 ± 0.13 is larger than typical values of ∼ -0.25 observed in DLAs at similar metallicities in the overall population of high-z DLAs (Rafelski et al. 2012;Quiret et al. 2016;De Cia et al. 2016) or in the sub-population of ESDLAs (Noterdaeme et al. 2014). Similarly, high depletion level, [Fe/P] = -1.21 +0.16 -0.17 , (for such a low-metallicity system) has recently been observed in the ESDLA (log N(H i) = 22.4 ± 0.1) towards J 2140-0321 (Noterdaeme et al. 2015b). The two systems . SDSS spectrum of J 0843+0221(black) compared to the unreddened SDSS composite (blue) and the composite spectrum reddened by our best-fitting value A V = 0.07. The cyan points indicate regions of the spectrum that were used to constrain the continuum fit. The inlay shows the distribution of measured A V for the control samples of spectra without having DLAs. This provides the systematic error of A V measurement. The red line indicates the extinction measured in the DLA at z=2.7865 towards J 0843+0221. In spite of the large H i and H 2 column densities the measured upper limit A V < 0.1 is a factor of ∼ 30 lower than the value expected from the extrapolation of relation between A V and hydrogen column density obtained in the SMC.

Mg Si P Ti Cr Mn Fe Ni -2 -1 0 [X/Zn] Figure 10
. Depletion pattern of different elements in the DLA at z=2.7865 (red dots with error bars). Blue dotted, magenta dash-dotted and black dashed lines show the typical depletion pattern inferred for MW halo, warm and cold ISM, respectively (Welty et al. 1999). The depletion pattern is observed to be mild compared to that of the cold ISM gas in the Galaxy.

actually share similar properties: the DLA towards J 2140-0321 also has low metallicity ([P/H]∼ -1.05) and high H 2 column density (log N(H 2 ) = 20.13 ± 0.07).

It is evident that the dust content scales not only with depletion, but also with the gas-phase metallicity. Fig. 11 shows the location of H 2 -bearing DLAs in the metallicity-N(H i) plane. There are 40 published z > 0 quasar H 2 absorption systems (Levshakov & Varshalovich 1985;Ge & Bechtold 1997;Ge et al. 2001;Ledoux et al. 2002;Levshakov et al. 2002;Reimers et al. 2003;Ledoux et al. 2003;Cui et al. 2005;Ledoux et al. 2006;Petitjean et al. 2006;Noterdaeme et al. 2007Noterdaeme et al. , 2008;;Srianand et al. 2008b;Jorgenson et al. 2009;Malec et al. 2009;Noterdaeme et al. 2010;Jorgenson et al. 2010;Srianand et al. 2010;Fynbo et al. 2011;Guimarães et al. 2012;Srianand et al. 2012;Crighton et al. 2013;Oliveira et al. 2014;Balashev et al. 2015;Noterdaeme et al. 2015b;Muzahid et al. 2015Muzahid et al. , 2016;;Krogager et al. 2016;Noterdaeme et al. 2017) and this paper. Also there are four H 2 -bearing systems detected in GRB afterglows (Prochaska et al. 2009;Krühler et al. 2013;Friis et al. 2015;[START_REF] Moorwood | [END_REF]. As observed by several authors (e.g. Petitjean et al. 2006), all H 2 detections are located in the upper part of the diagram, above a threshold value for overall metallicity ∼0.02 Z (i.e H 2 -bearing DLAs have statistically higher metallicities than overall population of DLAs). Among H 2 -bearing DLAs, the system towards J 0843+0221 has one of the lowest metallicity in spite of having the largest H 2 column density. This could be due to a dust selection effect against DLAs that would have similar hydrogen column densities but higher metallicities. Indeed, in our Galaxy, absorption systems with log N(H) ∼ 10 22 cm -2 are expected to have A V 3. Such an extinction would definitely preclude such systems from the SDSS. There are only two exceptional cases: the measurement in the DLA associated with GRB 080607 gives A V ≈ 3.2 mag (Prochaska et al. 2009) and the z = 0.89 lensing system towards PKS 1830-211 where very high column density of H 2 is expected given the plethora of other detected molecules [START_REF] Muller | [END_REF](Muller et al. , 2014)). In this latter peculiar case, the background blazar is completely extinguished in the optical bands and was selected from its millimetre emission and the foreground lensing galaxy (see Wiklind & Combes 1996).

As mentioned previously, chlorine is tightly connected with H 2 , with all chlorine is neutral in the H 2 -bearing gas while is singly ionized elsewhere. In addition, we expect little depletion of chlorine (Neufeld & Wolfire 2009) for the low metallicity and reddening (see Sect. 3.8) measured in this system. We can therefore use the abundance of neutral chlorine to constrain the metallicity in the H 2 -bearing gas assuming its molecular fraction, or, conversely, make assumptions on the metallicity in that phase to constrain the molecular fraction (Balashev et al. 2015). The measured log N(Cl i) = 13.63 +0.20 -0.05 translated to a molecular fraction in the H 2 phase of 0.8, under the assumption that the metallicity in that phase equals the overall metallicity. It is also likely that the metallicity is actually higher in the H 2 -bearing gas, meaning that the derived molecular fraction can be considered as a lower limit, i.e. the

Kinetic temperature

In Sect. 3.3, we found the excitation temperature of the first excited level of H 2 to be T 01 = 123 +9 -8 K. This is known to provide a good estimate for the kinetic temperature of the gas (e.g. Roy et al. 2006). In Fig. 12, we compare this value with measurements of T 01 in both distant DLAs and in the local Universe. The data for the local Universe were taken from Savage et al. (1977) for the MW, Gillmon et al. (2006) 2015)) but higher than what is measured locally. However, as can be seen in Fig. 12, T 01 clearly decreases as a function of N(H 2 ) so that the difference in the 'average' temperatures between the local and distant measurements is mostly due to a difference in the probed H 2 column densities. In the overlapping range of column densities (e.g. around log N(H 2 ) ∼ 18.5 -20.5), temperature measurements agree well between z = 0 and z > 0.

It is striking that the T 01 value measured in the DLA towards J 0843+0221 is significantly higher than T 01 in other H 2 strong systems. For log N(H 2 ) > 21, we would expect T 01 50 K when we observe more than twice this value. A possible explanation of this significant difference resides in the low metallicity of this system. The cooling through the metal line emission (mostly [C ii]157µm) can be much less efficient than at high metallicity and the equilibrium is shifted to higher temperatures. This possibility is supported by recent calculations of the metallicity dependence of temperature states in the neutral medium (e.g. Glover & Clark 2013). The measurement in the DLA towards J 0843+0221 is at least a factor of 2 higher than the value expected from extrapolating the trend of decreasing T 01 with increasing N(H 2 ).

However, there is no apparent trend in the measured T 01 as a function of metallicity (see colour code of the circles in Fig. 12). Since the system we are studying is unique both in terms of metallicity and N(H 2 ), it would be highly desirable to obtain additional measurements of molecular gas temperature in low-metallicity environments.

HD/2H 2 ratio

The relative abundance of HD to H 2 molecules is known to be sensitive to the local conditions in the gas (Balashev et al. 2010;Liszt 2015). Because of the low metallicity of our system, we can expect the astration of deuterium to be negligible (see Dvorkin et al. 2016). The intrinsic deuterium abundance should therefore be close to the primordial value, the latter being well constrained from the observations of atomic D i and H i lines in Lyman limit systems and (sub)-DLAs (e.g. Pettini & Cooke 2012;Noterdaeme et al. 2012b;Cooke et al. 2016;Balashev et al. 2016) or inferred from the analysis of the primordial anisotropy of CMB using the results of the big bang nucleosynthesis (Ade et al. 2014). We measure N(HD)/2N(H 2 )= 6.9 +2.9 -3.8 × 10 -5 . This value is consistent within errors with the measurements in Q 1232+082 and Q 0812+3208 and the primordial D/H ratio. Such a high HD/2H 2 ratio is expected to occur when all the D and H are in their molecular form. The abundance of neutral chlorine also indicates a very high molecular fraction (close to unity) in the H 2 -bearing gas. However, since the measured HD/2H 2 ratio is marginally higher than the intrinsic isotopic ratio (D/H), one can assume that the molecular fraction of hydrogen is slightly below unity.

Number density

The relative populations of the fine-structure levels of C i, C ii, Si ii and the rotational levels of HD can be used to estimate the number density in the medium. For the typical conditions of the cold neutral medium, these excited levels are mostly populated by collision with molecular and atomic hydrogen (depending on the actual molecular fraction in the medium) and depopulated by collisions and spontaneous transitions. In most cases, collisions with electrons can be neglected. For example, if the Si ii* is mostly populated by electron collisions, then the measured Si ii*/Si ii ratio would imply an ionization fraction > 0.01 which is very unlikely given the cold and molecular nature of the cloud. Similarly, radiative pumping by the UV background field is only important at low densities (Silva & Viegas 2002). A very high value of the UV field, more than 100 times the mean Galactic value (Draine & Bertoldi 1996) would be needed to explain the observed populations of the C i fine-structure levels. Therefore we include neither the electron collision nor the radiative pumping by UV photons in the analysis. This allows us to reduce the number of parameters in the calculations. In turn, we do take into account direct excitation by the CMB radiation. This does not add any additional free parameter since the redshift of the system (z ≈ 2.7865) is known and hence the value of the CMB field is fixed under standard cosmology. In any case, we found that radiative transitions by the CMB provide only a very small contribution to the excitation of the above species, even for C i, whose transition energies between fine-structure levels are close to the CMB temperature at the DLA redshift.

To calculate how the relative populations of the fine-structure levels of C i, C ii, Si ii and the rotational levels of HD depend on the physical conditions in the gas we used the code which is similar to the PopRatio (Silva & Viegas 2001, 2002) and RADEX (van der Tak et al. 2007) codes. We used data mostly provided along with PopRatio, which includes energy levels [START_REF] Moore | NSRDS-NBS[END_REF], transition probabilities (Galavis et al. 1997) and collisional coefficients (Schroder et al. 1991;Staemmler & Flower 1991;Abrahamsson et al. 2007) for C i and energy levels [START_REF] Moore | NSRDS-NBS[END_REF], transition probabilities (Galavis et al. 1997) and collisional coefficients (Wiesenfeld & Goldsmith 2014) for C ii and Si ii. The atomic data for HD were taken from Flower et al. (2000).

Neglecting collision with electrons, UV pumping and fixing the CMB field, we are left with three parameters to calculate the relative populations of excited levels: kinetic temperature, number density and molecular fraction. We then perform calculations using two extreme cases for the local molecular fraction: (i) f = 0.3, which corresponds to the overall observed molecular fraction and is hence a strict lower limit since it assumes all H i in the DLA to be associated with the H 2 -bearing phase and (ii) f = 1, where the only collision partner is H 2 . Note that the abundances of neutral chlorine and HD indicate that we should be closer to the fully molecularized case. We also assume the medium to be homogeneous, i.e. the calculation is performed using a one-zone model: the molecular fraction, number density and temperature are constants throughout the cloud.

To compare the calculated values with the observed populations we used a maximum likelihood function. Since for some species we do not have good observational constraints on the column densities for individual components, in the following we consider the total column densities instead. This should give a good approximation of average physical conditions in the gas. We note indeed that for C i, the population ratios in both components -on which we have good constraints -are very similar. While we have a measurement of C ii in its excited state (i.e. C ii*), we do not have a measurement of the ground state. From the zinc metallicity, and assuming solar intrinsic abundances, we estimate log N(C ii) = 16.93. We considered this value as an upper limit since carbon is likely more depleted in the diffuse ISM than zinc. Since Si ii* is only ten- . Constraints on the number density and kinetic temperature obtained from different excited species. Constraints from C i and HD are in orange and purple, respectively, while the excited ionized silicon and carbon provide respectively upper and lower limits on the density for a given temperature. The red region shows the joint estimate when also considering the excitation of H 2 . The top and bottom panels are for molecular fraction f =1.0 and f =0.3, respectively. The estimated number density in the DLA at z=2.7865 is higher than typical values obtained in the diffuse ISM both in the local Universe and at high redshifts.

tatively detected, we also treat Si ii* column density as an upper limit. The results of calculations are shown in Fig. 13.

The allowed regions in the temperature-density plane for the different species overlap for the two assumed molecular fractions, but in different regions of that plane. When adding the constraint on the kinetic temperature independently derived in Sect. 4.2, the number density is found to be log n = 2.58 +0.03 -0.06 and log n = 2.42 +0.06 -0.03

(where n is in units of cm -3 ) for the f = 1 and f = 0.3 cases, respectively. This is significantly higher than what is generally estimated in the neutral medium n ∼ 0.1 -10 cm -3 and also higher than what is usually derived in other molecular-rich DLA systems (e.g. Noterdaeme et al. 2017).

Thermal pressure

Using the estimates of the number density and kinetic temperature of the cloud, we derive the thermal pressure, P = nT = 4.7 +0.5 -0.6 × 10 4 cm -3 K. While the thermal pressure is not necessarily dominating the total pressure (that includes radiation, turbulence and magnetic field pressures), it is expected to correlate with other parameters of ISM, like turbulent pressure (Joung et al. 2009) as well as heating and star formation rates (Ostriker et al. 2010).

In our Galaxy typical thermal pressures are found to be in the range from 2 × 10 3 to 10 4 K cm -3 (Jenkins & Tripp 2011). Similar values were recently obtained in a sample of nearby galaxies (Herrera-Camus et al. 2017). However, Welty et al. (2016) measured slightly higher values in the Magellanic Clouds. This is consistent with the theoretical predictions that higher pressures are required for stable cold, neutral clouds in environments with low metallicity, low dust-to-gas ratios and enhanced radiation fields. Higher pressures are therefore also expected in DLAs at high redshifts. However, the measurements by Jorgenson et al. (2010) using C i and Neeleman et al. (2015) using C ii and Si ii give pressures consistent with MW values. We note that the measurements by Neeleman et al. (2015) are consistent with an increase of thermal pressure with increasing N(H), despite a large uncertainty on individual measurements.

We applied the same calculations as in the previous section to other known high-redshift C i-bearing systems to derive the pressure estimates shown as green points in Fig. 14. Note that we assumed the kinetic temperature for each system to equate the T 01 H 2 excitation temperature and that the cloud is fully molecularized. This a rather strong assumption that should in principle be tested for each individual cloud. Since the collisional coefficients for H 2 are slightly higher than for H i, the densities derived should be treated as upper limits and so the derived pressures. We see that most of the known C i-bearing systems have pressures P ∼ 10 4 K•cm -3 . Both the system studied here and the other H 2 absorption system in the ESDLA towards J 2140-0321 (Noterdaeme et al. 2015b) have higher pressure. This is consistent with the interpretation that ES-DLAs correspond to lines of sight that are statistically closer to the galaxy centre (Noterdaeme et al. 2014), where we expect the pressure to be higher.

UV flux

While the excitation of the high rotational levels of H 2 provides an indication of the strength of the UV field, it is not easy to obtain a quantitative estimate due to the complexity of the excitation. A main difficulty is that radiative pumping and self-shielding effects lead to inhomogeneous distribution of the population of different rotational levels of H 2 throughout the cloud (see e.g. Abgrall et al. 1992). High-J levels are dominant at the cloud surface, while low-J levels are mostly found in the shielded interior, where the physical conditions can be different (Balashev et al. 2009). This could explain the increasing Doppler parameter with increasing rotational level, as often observed for H 2 (Lacour et al. 2005;Noterdaeme et al. 2007;Balashev et al. 2009). The effects of geometry, which is likely more complex than the simple slab or spherical that is generally assumed, as well as the effects from turbulence can therefore play important roles.

Keeping this in mind, we consider the H 2 column densities in the different rotational levels (see Table A1) summed over the three components. In Fig. 15 we compare this 'average' excitation with that of other high-N(H 2 ) absorption systems: embedded in a typical UV field of the MW (Gry et al. 2002), a highly excited case towards the O9.5V runaway star HD 34078 (Boisse et al. 2005) and the ESDLA at z = 2.34 towards J 2140-0321 (Noterdaeme et al. 2015b). Both high-z DLAs show H 2 excitation higher than typical MW values, but lower than towards the OV star. Interestingly, the excitation towards the latter has been successfully modelled using two components, one little affected by the star and contributing to most of the J = 0 and 1 lines and another one with high density (n ∼ 10 4 cm -3 ) strongly illuminated by the star. While the situation is unlikely to be such extreme in the DLAs, the excitation diagram indicates an enhanced UV background and larger number densities compared to typical MW values. Indeed, photoionization models also indicate log n ∼ 2.5-3 towards J 2140-0321. For J 0843+0221 it can also be seen that the J = 2 and J = 3 rotational levels are consistent with the excitation temperature measured from J = 0 and J = 1, T 01 = 123 +9 -8 K, which suggests that the J = 2 and J = 3 levels are predominantly populated by collisions in the medium, as expected in strongly self-shielded regions (Abgrall et al. 1992).

An alternative way to constrain the UV field is the C ii* method (Wolfe et al. 2003). Assuming thermal equilibrium, the heating rate through photoelectric effect on dust grains equates the cooling rate dominated by the [C ii]158µm emission. From the column densities of C ii*, H i and H 2 we find that the cooling rate per hydrogen atom by C ii fine-structure transition is log l c = -26.0 +0.2 -0.2 , where l c is in units of erg s -1 (Pottasch et al. 1979). This is slightly above the values measured for typical DLAs with log N(H i) < 21.5 (Wolfe et al. 2008). Using the dust-to-gas ratio κ ∼ 0.03 and estimating the grain heating efficiency ∼ 0.03, we find that the UV field is a factor of 10-30 higher than the Habbing UV field. This is in qualitative agreement with the high-excitation diagram of H 2 . This is also in indirect agreement with the low CO/H 2 abundance since the photodestruction of CO scales with the UV flux.

On one hand, since our system has a very low dust content, the above method implies that a large UV flux is required to produce the estimated heating rate. This need not be the case if the heating is mainly due to cosmic rays (Dutta et al. 2014;Noterdaeme et al. 2017). In that case, balancing the cooling implies a cosmic ray ionization rate of about 3 × 10 -16 s -1 , i.e. 30 times higher than the canonical Galactic value 10 -17 s -1 . Note that usually higher values of cosmic ray ionization rate than canonical value are observed in the dense clouds (Dalgarno 2006). On the other hand, it is also 0 1000 2000 3000 4000 5000 expected that the cosmic ray ionization rate correlates with the UV flux, since both arise from star formation activity. Therefore, in the real situation, both dust grain heating and heating by cosmic rays are likely equally important. The factor of 30 derived for the UV flux and cosmic ionization rates alone are therefore upper limits and their actual values should be roughly a factor of 2 lower.

Energy, cm -1 -10 -9 -8 -7 -6 -5 -4 -3 -2 -1 0 log N i /g i -log N 0 /g 0 T 01 = 123K J0843+0221 HD 34078 diffuse MW J2140-0321

DISCUSSION

Properties of H 2 cloud

Using the derived physical conditions, we can discuss some properties of the H 2 -bearing cloud towards J 0843+0221 and compare them with measurements in the local Universe. In the following, we assume that the H 2 cloud is spherical and homogeneous. We again consider the two limiting cases for the molecular fraction, f = 0.3 and f = 1. The derived parameters of the cloud are given in Table 3. Given the hydrogen column density N H and volumic density, n, we find that the cloud size along the line of sight, l = N H /n, is between ∼ 1.4 and ∼ 12 pc (for f = 1 and f = 0.3, respectively). The upper limit is comparable with the characteristic sizes of giant molecular clouds (GMCs) observed in our Galaxy and in the Magellanic Clouds. The mass of the cloud can be estimated as M = 4/3πR 3 ρ, where R is the radius of the cloud, ρ = µn, where µ is the mean mass weight taken to be 2.15m H and 1.23m H (m H is the hydrogen atomic mass) for f = 1 and f = 0.3, respectively. Since the impact parameter between the line of sight and the cloud centre is unknown, we use the most probable radius of the cloud R = 3/4 • l. We then obtain masses in the range (1 -300) × 10 2 M , which corresponds to the lower end of the mass distribution of MW molecular clouds selected from 13 CO emission (Roman-Duval et al. 2010). The high allowed values reach the lower limit of ∼10 4 M , which usually is considered for GMCs (McKee & Ostriker 2007).

We found that the surface density, Σ = M/πR 2 (which in our

0 1 2 3 4 5 6 log M, M ⊙ 1 2 3 log Σ, M ⊙ pc -2 R = 0 .1 p c R = 1 p c R = 1 0 p c
Milky Way CO dark gas case is equal to µN H ), is constrained to be in the range ∼ 28 -112 M pc -2 (depending on the adopted molecular fraction).

13 CO HCO+ J 0843+0221 J 2140-0321 J 0000+0048 n = 1 0 1 c m -3 n = 1 0 2 c m -3 n = 1 0 3 c m -3 Figure 16.
In Fig. 16 we show the Σ-M dependence of molecular clouds. The allowed range for the DLA towards J 0843+0221 is located slightly below most of the HCO + -and 13 CO-selected clouds in the Milky-Way (Roman-Duval et al. 2010; Barnes et al. 2011). We note in passing that HCO + -selected clouds have higher Σ-values than 13 CO-selected ones for a given mass, due to HCO + tracing highdensity clumps in molecular clouds. The observed number density towards J 0843+0221 is comparable to the average value obtained in 13 CO-selected clouds.

While the upper limit of Σ = 112 M pc -2 ( f = 0.3 case) resembles the typical surface densities of MW clouds, the lower limit Σ = 28 M pc -2 ( f = 1 case) corresponds to the lower end of 13 CO-selected clouds in the MW. We note that in the case of the molecular cloud towards J 0000 + 0048 (Noterdaeme et al. 2017), the surface density remains lower than the MW clouds. In the case of the ESDLA towards J 2140 -0321 (log N(H i) = 22.4 Noterdaeme et al. 2015b), the upper range overlaps with MW points but this corresponds to the extreme case where all H i in the DLA is assumed to be associated with the H 2 phase, while modelling suggests that it is not (Noterdaeme et al. 2015b). Here, independently of the exact actual molecular fraction, the H 2 cloud appears to probe a low-metallicity analogue of 13 CO-selected molecular cloud.

Based on their data, Roman-Duval et al. (2010) found that most 13 CO clouds are gravitationally confined with virial parameter α vir < 1. To estimate the virial parameter, α vir = 5σ 2 R/(GM), one needs to know the 1D velocity dispersion in the medium, σ. Using σ ∼ 3 km s -1 from the Doppler parameter of the main Cl i component, we obtain α vir between ≈ 100 and 3 for f = 1 and f = 0.3, respectively. This means that the cloud is probably not self-gravitating, but pressure confined. This is especially true if we take into account the fact that b turb measured using low-ionization metals is significantly higher than that based on Cl i. The range of Doppler parameters suggests that turbulence in the cloud is supersonic with Mach number > 3.

Timescales

The low metallicity and dust-to-gas ratio in the DLA studied here as well as its marked difference with other H 2 -bearing DLAs call for a discussion of the times-cales of relevant chemical and thermodynamic processes. Indeed, the rates of these processes usually depend on the metal/dust content, and even for the solar-metallicity gas, the steady-state balance is sometimes questioned.

In the cold ISM, H 2 forms only on dust grains. The rate of such H 2 formation is relatively slow and therefore the dust content plays a crucial role. This can be estimated as t = 1/(Rn 2 D), where R is the formation rate coefficient, n is the number density in the cloud and D is the dust-to-gas ratio compared to the Galaxy one. We assumed that D is proportional to the metallicity measured in J 0843+0221 DLAs, i.e., we took D = 1/30. For the formation rate coefficient, we used the generally adopted value derived from Copernicus observations, R = 3 × 10 -17 cm 3 s -1 . Using our derived number density, we obtain formation time-scale to be ∼ 60 and ∼ 90 Myr for f = 1 and f = 0.3 cases, respectively (see Sect. 5.1). The dynamical time-scale can be estimated as t dyn = l/b turb , where l is the size of the cloud. Using the lower limit on the Doppler parameter b turb = 3 km s -1 from Cl i we obtained t dyn = 0.5 and 4.0 Myr (for f = 1 and f = 0.3, respectively). The free-fall time can be estimated as t ff = (3π/32/G/ρ) 1/2 and equals 1.8 and 2.7 Myr (for f = 1 and f = 0.3, respectively). Since the H 2 formation timescale is much larger than the dynamical and free fall time, this independently confirms that the H 2 -bearing cloud is likely pressure confined. We note that while a high dust content reduces the H 2 formation time-scale, the number density can have major importance since t ∝ n -2 . It has two direct consequences: first, in low-metallicity environments, H 2 may trace only the densest regions in the patchy medium of DLAs. High number densities will obviously severely reduce the covering fraction of the H 2 -bearing gas. Secondly, the formation time-scale of H 2 can be shortened by supersonic turbulences that can create pockets of compressed gas, where H 2 is quickly formed, and consequently redistributed to the adjacent medium.

In Sect. 4.2, we showed that the H 2 -bearing gas in the DLA at z=2.7865 has a kinetic temperature T ∼ 122K which is significantly higher than what is expected at such H 2 column density. We supposed that such high temperature was likely attributed to the low metallicity in the medium. The cooling time can be easily es-timated from the cooling rate, l c , derived from C ii*, as ∼ 2k B T/l c , where k B is the Boltzmann constant (the factor of 2 was chosen since the molecular fraction is not well constrained and the correct coefficients are 3/2 and 5/2 for monatomic and diatomic gas, respectively). This indicates a cooling time of about 0.1 Myr, which is at least an order shorter than the time-scales mentioned above. Hence the cloud is likely in thermal equilibrium and the high kinetic temperature is indeed due to the low metallicity and not to the formation of H 2 in a shock-compressed medium.

CO content, X CO and CO-dark gas

In spite of extremely high H 2 column density, we did not detect CO. This is not surprising given the facts that we measured a low associated extinction and possibly an enhanced UV flux: CO is easily photodestructed by UV photons in the diffuse medium and therefore at least moderate level of extinction A V 1 is required for this molecule to survive. Models indicate that the abundance of CO scales with the metallicity, but there are only scarce measurements at low metallicities (e.g. Shi et al. 2015). Our upper limit on N(CO) is therefore very valuable since it provides a stringent limit on X CO at Z ∼ 0.03Z for the first time.

Fig. 17 shows the dependence of the CO column density on that of H 2 measured in both the MW and high-redshift ISM. Values for high-z DLAs are from Srianand et al. (2008b); Noterdaeme et al. (2010Noterdaeme et al. ( , 2017) ) for quasar absorption systems and Prochaska et al. (2009) for the GRB-DLA6 . We also derived upper limits on the total CO column density in known strong H 2 absorption systems, using the covered A-X bands in the archival spectra and assuming that the population of CO rotational levels is dominated by the CMB excitation. We found N(CO) < 12.55 in the DLA towards Q 1232 + 082 with log N(H 2 ) = 19.57 ± 0.10 and N(CO) < 12.81 in the DLA towards Q 0812 + 3208 with log N(H 2 ) = 19.93 ± 0.05. A correlation between N(CO) and N(H 2 ) is naturally expected as we probe larger column densities of molecular gas but there is interesting deviation around the main trend. CO detections in DLAs tend to probe high CO/H 2 values. This could be due to selection effects since the measurements are close to the typical CO detection limits. In the case of detections, the high CO/H 2 ratio could result from a combination of small dust grains and high cosmic ray ionization rate (see Shaw et al. 2008;Rawlins et al. 2016;Noterdaeme et al. 2017). The case of the extreme H 2 absorber studied here is striking: our limit on N(CO) is 10 3 times lower than what the MW trend suggests or than the value measured in the GRB-DLA by Prochaska et al. (2009). This can be attributed to the much lower gas-phase abundances of C and O in the DLA towards J 0843+0221 together with the low dust content, which is crucial for CO since dust provides shielding from the UV field.

We calculated the lower limit on the X CO conversion factor for the optically thin regime (details can be found in Bolatto et al. 2013) to be X CO 2 -6 × 10 23 cm -2 /(km/s K), depending on the CO excitation temperature that must be in the range between the CMB temperature at z = 2.7865 and the kinetic temperature (from T 01 ). This value of X CO is 1000 -3000 times higher than the adopted MW value. Given this very high X CO value (or low CO/H 2 ratio), we can state that our high-z molecular cloud corresponds . The open and filled brown points correspond to direct (Sheffer et al. 2008;Burgh et al. 2009) and indirect (Federman et al. 1990) H 2 measurements in our Galaxy, respectively. The upper limit on N(CO) in the H 2 absorption system at z=2.7865 studied here is shown by a red circle and is a factor of 10 3 lower than values measured for similar H 2 column densities in the MW.

-1.5 -1.0 -0.5 0.0

[X/H] 10 0 10 1 10 2 10 3 10 4 X CO /X CO (MW), α CO /α CO (MW) Milky Way Leroy+11 Sandstrom+13 Israel+97 Shi+16 Shi+15 Elmegreen+13 J 0843+0221
Figure 18. Measurements of X CO and α CO conversion factors relative to adopted MW values (Bolatto et al. 2013) as a function of metallicity. The data for the local measurements of α CO are from Israel & P. 1997;Leroy et al. 2011;Sandstrom et al. 2013;Elmegreen et al. 2013;Shi et al. 2015Shi et al. , 2016. . to CO-dark gas. 7 Such a high X CO conversion factor is expected at low metallicity and enhanced UV flux in the medium, and any observational constraint is valuable. Fig. 18 summarizes the measurements of CO/H 2 conversion factors (X CO , α CO ) as a function of metallicity in different environments. While we are only able to put a lower limit on X CO that is similar to previous measurements at low metallicities, it is the first time a useful constraint is put at such 7 We note that, in the MW, the 'CO-dark gas' (or 'CO-faint gas', see Wolfire et al. 2010;Bolatto et al. 2013) is usually estimated to have Σ < 30 M pc -2 . This threshold value in the MW corresponds to log N(CO) ∼ 16.5 (see Fig. 17) a low metallicity (twice lower than any previous measurement in the MW).

Star formation rates

Following the formalism of Wolfe et al. (2008) and the estimates on the UV flux, we estimate the surface star formation rate, Σ SFR , in the DLA towards J 0843+0221 to be about 2 × 10 -2 M yr -1 kpc -2 , if stars are forming in situ. This value is comparable with the values of Σ SFR for the measured Σ H2 in local galaxies with kpc-resolution emission observations (Bigiel et al. 2011). This supports the idea that the H 2 -bearing system in J 0843+0221 probes a star-forming medium, similar to the GMCs in the local Universe.

Detecting more direct markers of star formation such as the nebular emission lines would be very valuable (see e.g. Krogager et al. 2017). Here, only the Ly-α line is covered by our UVES spectrum. This line can in principle be detected in the bottom of the damped Ly-α absorption that removes the flux of the quasar over a significant wavelength range (as in Noterdaeme et al. 2012a). However, we did not find any significant Ly-α emission in our spectrum. To estimate a limit on the Ly-α flux, we rebinned and coscaled the UVES spectrum with the flux-calibrated SDSS spectrum and obtained F Lyα < 0.8 × 10 -17 erg cm -2 s -1 (and Ly-α luminosity < 5.7 × 10 41 erg s -1 . Assuming no dust correction, this upper limit translates (see e.g. Rahmani et al. 2010) to SFR < 4 M yr -1 for in situ star formation activity.

CONCLUSION

We presented the analysis of the very high H 2 column density absorption system at z abs = 2.7865 identified in the spectrum of quasar J 0843+0221. We measured the total H 2 column density to be log N(H 2 ) = 21.21±0.02, which is a factor of 6 higher than what has been previously detected in intervening DLAs towards quasars. The relative population of J = 0, 1 levels of H 2 indicates a kinetic temperature in the H 2 -bearing medium of T = 123 +9 -8 K, which is at least twice higher than what is expected for such large H 2 column densities from extrapolating the trend seen in other measurements. This is probably due to the low metallicity of the medium. We find that metallicity in this ESDLA (log N(H i) ∼ 21.82 +0.11 -0.11 ) is one of the lowest, [Zn/H] ∼ -1.5, among H 2 -bearing DLAs. Although the measured metallicity is typical for DLAs at similar redshift, the depletion level in the DLA, [Fe/Zn] ∼ -1, is found to be significantly higher than in non H 2 -bearing DLAs at the same metallicity. This supports the idea that some amount of dust is required for the presence of H 2 . Even though the depletion and the total hydrogen column density are high, the dust content is still not enough to produce any significant associated extinction. We measure A V = 0.09 ± 0.10, i.e within the systematic uncertainties of the measurements. This is 30 times lower than what is expected using the average MW A V /N(H) ratio, but scales exactly with the measured metallicity.

The low metallicity implies a very low abundance of CO. We measure an upper limit on the column density of CO, log N(CO) < 13.35 that translates to a lower limit on X CO factor > 2 × 10 23 cm -2 / (km s -1 K). This is >10 3 times higher than the adopted MW value. While an increase of X CO with decreasing metallicity is expected, our system has the lowest metallicity in which X CO is constrained. The low CO and high H 2 content indicate that this system represents 'CO-dark/faint' gas, whose abundance in different galaxies is currently a matter of debate now (Grenier et al. 2005;Pineda et al. 2013;Smith et al. 2014).

The detection of fine-structure levels of C i, C ii, Si ii and the excitation of HD molecules allows us to constrain the number density in the cold gas to be ∼ 260 -380 cm -3 , depending on the adopted molecular fraction. Taking into account our estimate of the kinetic temperature, we found a relatively high thermal pressure in this system, ∼ 3 -5 × 10 4 K cm -3 , which is higher than typical pressures measured in the diffuse ISM at z = 0. Our measurement confirms the trend of increasing thermal pressure with increasing hydrogen column density. This is in line with the idea that high column density systems preferably probe gas at small impact parameters from the associated galaxy (Noterdaeme et al. 2014). The observed column density of C ii* implies a cooling rate by [C ii] 158 µm to be l c ∼ 10 -26 erg s -1 , which is slightly higher than what is usually observed in DLAs. This, together with the excitation diagram of H 2 rotational levels and the low CO abundance, indicates a strong UV field, roughly a few tens times higher than the Habing field.

The molecular fraction in the H 2 -bearing medium can be in range between 0.3 and 1. Assuming homogeneous density, we estimate the size of this cloud to be between ∼1 and 10 pc along the line of sight. The cloud is most probably pressure confined, based on the estimation of the virial parameter and characteristic timescales. Assuming spherical shape of the cloud, we derive the total gas mass between (1-300)×10 4 M , while the higher range reaches the commonly used threshold for GMC. The corresponding surface density is between ∼30 and 110 M pc -2 . This is very similar to the 13 CO-selected clouds studied in our Galaxy. This means that for the first time in quasar absorption line studies, we confidently probe a starforming molecular cloud similar to GMCs in the MW. A1 presents the best-fitting parameters to the absorption lines in the neutral (atomic and molecular) phases. We used the three components derived from Cl i with fixed redshifts (given in Table A1) to fit the J ≥ 2 lines of H 2 and obtain column densities in each component after locally adjusting the continuum level using spline interpolation. each component we tied the Doppler parameters (b) of the different rotational levels. The results of the analysis are presented in Table A1 and the profiles of fitted H 2 lines are shown in Fig. A2. Since it has been shown by different authors that the Doppler parameters of H 2 lines tend to increase with the rotational level (Lacour et al. 2005;Noterdaeme et al. 2007;Balashev et al. 2009), we also performed a fit with independent b-values for each rotational level. We indeed found the increase of Doppler parameter with J for the two stronger components of H 2 absorption system. Nevertheless, the derived column densities remain consistent with the fit using tied Doppler parameters, although the uncertainties become much larger due to this additional freedom.

The excitation diagrams of H 2 for the components A, B and C are shown in Fig. A3. The J = 0, 1 points represent the total column densities in the system and so are strict upper limits for the individual components. Focusing on J ≥ 2 levels, it appears that the components A and B have similar diagrams, with slightly larger column densities for the component B but slightly higher excitation temperature for the component A. The component C is the weakest of the three, but apparently has much higher excitation, as suggested by the shallow slope between J ≥ 2 levels. Interestingly, this component can therefore contribute significantly to the total high-J column densities, but little to the overall H 2 column density, which mainly arises from the components A and B. . Excitation diagrams of H 2 rotational levels for A (green circles), B (orange circles) and C (violet circles) components of H 2 absorption system at z=2.786. The black points correspond to the J = 0 and 1 levels, where we are only able to measure total column densities.

Figure 1 .

 1 Figure 1. A portion of our UVES spectrum around the DLA at z = 2.786. The solid red line represents the fitted Ly-α profile with associated uncertainty as shaded region. The red dashed line shows the reconstructed continuum, which was fitted together with the H i line profile. The measured H i column density is log N(H i)=21.82 ± 0.11, i.e. this DLA corresponds to a sub-class of DLA known as extremely strong DLAs (ESDLAs; Noterdaeme et al. 2014).

Figure 2 .

 2 Figure 2. Fit to the Cl i absorption lines at z=2.7865. In this and following figures presenting fits to absorption lines, pixels used (resp. not used) to constrain the fit are represented in black (resp. grey). The synthetic profile is shown in red and the contribution from each component is shown in green, orange and purple for components A, B and C, respectively. The small filled circles at the top of each panel show the residuals and the ±2 σ error levels are shown by green lines. We used component structure measured in Cl i lines to fit H 2 and HD lines.

Figure 4 .

 4 Figure 4. Fit to the HD lines in the J = 0, 1 rotational levels. Lines are as in Fig. 2. This is the second detection of lines from J=1 rotational levels of HD. Due to saturation of HD lines and low S/N the measured HD column density has large uncertainties.

Figure 5 .

 5 Figure 5. Fit to the C i absorption lines. Lines are as in Fig. 2.

Figure 6 .

 6 Figure 6. Fit to the low-ionization metal lines. Lines are as in Fig. 2.

Figure 7 .

 7 Figure 7. Portion of the UVES spectrum around the positions of S ii absorption lines. The red lines show the calculated S ii profiles based on the best-fitting parameters of S ii λ1250. The dashed violet line indicates the local QSO continuum used for S ii line fitting. The horizontal blue dashed line indicates the observed bottom of the S iiλ1259 line, which is located at ∼20% of the total flux in that region. Inner panel: portion of the UVES spectrum around Lyα and N v emission lines. Green dashed line indicates the estimate of the intrinsic QSO continuum. The violet box defines the region shown in the outer panel.

Figure 8 .

 8 Figure 8. Fit to C ii* and Si ii* lines, associated with DLAs at z = 2.7865.

  Figure9. SDSS spectrum of J 0843+0221(black) compared to the unreddened SDSS composite (blue) and the composite spectrum reddened by our best-fitting value A V = 0.07. The cyan points indicate regions of the spectrum that were used to constrain the continuum fit. The inlay shows the distribution of measured A V for the control samples of spectra without having DLAs. This provides the systematic error of A V measurement. The red line indicates the extinction measured in the DLA at z=2.7865 towards J 0843+0221. In spite of the large H i and H 2 column densities the measured upper limit A V < 0.1 is a factor of ∼ 30 lower than the value expected from the extrapolation of relation between A V and hydrogen column density obtained in the SMC.

Figure 11 .

 11 Figure 11. Metallicity versus H i column density in DLAs. The green and blue circles indicate the H 2 -bearing DLAs at high redshifts detected towards QSO and in GRB afterglows, respectively. The size of each circle depends on H 2 column density. The grey crosses correspond to non-H 2bearing DLAs. The dashed black lines show constant value of A V calculated in the assumption to fixed scalefactor between metallicity and dust abundance corresponds to the MW measurements. The shaded region reflects systematic bias due to obscuration.

  for high-latitude MW clouds and Welty et al. (2012) for the SMC and LMC. The excitation temperature in the DLA at z=2.7865 is similar to what is typically measured in high-z and low-z H 2 -bearing DLAs ( T 01 ∼ 150 K; Srianand et al. (2005); Muzahid et al. (

Figure 12 .

 12 Figure 12. Measurements of the H 2 excitation temperature, T 01 at different N(H 2 ). The grey points represent measurements in the Local Group (MW, SMC and LMC). Filled circles represent measurements at z > 0, where the colour corresponds to the average metallicity of each associated DLA. The measurement in the DLA towards J 0843+0221 is at least a factor of 2 higher than the value expected from extrapolating the trend of decreasing T 01 with increasing N(H 2 ).

  Figure13. Constraints on the number density and kinetic temperature obtained from different excited species. Constraints from C i and HD are in orange and purple, respectively, while the excited ionized silicon and carbon provide respectively upper and lower limits on the density for a given temperature. The red region shows the joint estimate when also considering the excitation of H 2 . The top and bottom panels are for molecular fraction f =1.0 and f =0.3, respectively. The estimated number density in the DLA at z=2.7865 is higher than typical values obtained in the diffuse ISM both in the local Universe and at high redshifts.

Figure 14 .

 14 Figure14. Estimated thermal pressure versus total hydrogen column density. Circles represent values derived from the C i fine-structure levels together with the H 2 excitation temperature, T 01 . The size of each circle corresponds to the H 2 column density. Squares correspond to the average values found in the Local Group (MW, SMC and LMC). Measurement in the DLA towards J 0843+0221 is consistent with an increase of the thermal pressures with increasing log N(H).

Figure 15 .

 15 Figure 15. Comparison of H 2 rotational excitation in different H 2 absorption systems. The black dashed line correspond to the temperature T 01 = 123 K measured using J=0,1 levels of H 2 in J 0843+0221. The excitation of H 2 suggests the UV flux is enhanced in the H 2 -bearing medium of the DLA at z=2.7865.

  Figure17. Dependence of N(CO) to N(H 2 ) column densities for highredshift DLAs detected in quasar and GRB afterglow (green and blue points, respectively). The open and filled brown points correspond to direct(Sheffer et al. 2008;Burgh et al. 2009) and indirect(Federman et al. 1990) H 2 measurements in our Galaxy, respectively. The upper limit on N(CO) in the H 2 absorption system at z=2.7865 studied here is shown by a red circle and is a factor of 10 3 lower than values measured for similar H 2 column densities in the MW.

Fig

  Fig.A1presents the ionization structure of the DLA at z = 2.7865, which is found to be quite typical of what is seen in other DLAs. TableA1presents the best-fitting parameters to the absorption lines in the neutral (atomic and molecular) phases.We used the three components derived from Cl i with fixed redshifts (given in TableA1) to fit the J ≥ 2 lines of H 2 and obtain column densities

a

  The redshifts of the components A,B and C were fixed during the fit and correspond to the best-fitting values from fit to Cl i lines (see the text). b To fit H i and H 2 J=0, 1 lines we used one-component model, due to high saturation of the absorption lines. c To fit low-ionization metal absorption lines we tied b parameters using microturbulence Doppler parameter, b turb , and kinetic temperature, T , using standard relation (see text). d The redshifts of the components A,B and C were fixed during the fit and correspond to the best-fitting values from fit to low-ionization metal lines (see the text).

Figure A2 .

 A2 Figure A2. Fit to the H 2 absorption lines. Lines of transitions from rotational levels J = 2 to 9 of the ν = 0 vibrational level. Lines are as in Fig.2.

Table 1 .

 1 and at high redshifts Journal of the observations † For this exposure calibration files were not attached

	No.	Date	Starting time (UT)	Setting	Exposure (s)	Airmass	Seeing (arcsec)
	1	07.01.2014 †	07:19:05	437+760	4800	1.162→1.377	1.00→0.78
	2	28.01.2014	03:06:44	437+760	4800	1.265→1.131	0.83→1.10
	3	24.02.2014	01:07:58	390+564	4800	1.306→1.141	0.97→0.92
	4	26.02.2014	01:34:21	437+760	4800	1.211→1.122	0.81→1.00
	5	26.02.2014	02:58:31	437+760	4800	1.121→1.188	0.98→1.24
	6	27.02.2014	01:32:14	390+564	4800	1.207→1.122	0.99→0.96
	7	27.02.2014	02:55:32	390+564	4800	1.121→1.190	0.94→0.74
	8	27.02.2014	04:16:46	390+564	4800	1.192→1.466	0.69→0.75
	9	02.03.2014	01:02:50	437+760	4800	1.249→1.127	0.69→0.84
	Total				43200

Table 2 .

 2 Summary of the overall gas-phase abundances

	Species (X)	[X/H]
	Cl i	-1.86 +0.21 -0.10
	Mg ii	-1.65 +0.09 -0.10
	Si ii	-1.73 +0.08 -0.10
	S ii	-1.54 +0.08 -0.08
	Ti ii	-1.96 +0.10 -0.11
	Cr ii	-2.40 +0.09 -0.10
	Mn ii	-1.75 +0.10 -0.10
	Fe ii	-2.53 +0.08 -0.10
	Ni ii	-2.52 +0.09 -0.10

  Estimates of surface density of cold gas, Σ, versus mass of the cloud M. The blue and green points correspond to measurements of molecular clouds in MW using 13 CO (Roman-Duval et al. 2010) and HCO + (Barnes et al. 2011). The lines connect the values for extreme cases (only H 2 or all H i associated with H 2 ) for three strong H 2 -bearing DLAs. The cloud associated with DLA towards J 0843+0221 is most likely similar to 13 CO-selected cloud in the MW.

Table 3 .

 3 Properties of the H 2 -bearing cloud at z=2.7865. All values given above are only indicative and do not take into account the uncertainties, which in some cases can be large.

	Parameter	f = 1	f = 0.3
	n, 10 2 cm -3	3.8	2.6
	P, 10 4 K cm -3	4.7	3.2
	l, pc	1.4	12
	Σ, M pc -2	28	112
	M, M	94	2.9×10 4
	α vir	115	3
	t H 2 , Myr	63	91
	t ff , Myr	1.8	2.7
	t dyn , Myr	0.5	4.0
	t cool , Myr	0.1	0.1

  Comparison of the line profiles of different ionization species to probe the ionization structure of the DLA at z=2.7865. We keep only lines of interest, not used in the fit. The orange vertical line corresponds to redshift z=2.786574 (component B). The ionization energy of the species is given in each panel.

		1
		0	Cl I 1347				12.968 eV
		1
		0	O I 1302				13.618 eV
		1
		0	Fe II 2382				16.199 eV
		1
		0	Si II 1304				16.346 eV
		1
		0	Al II 1670				18.829 eV
	Normalized flux	0 1	C II 1334				24.384 eV
		1
		0	Al III 1854				28.448 eV
		1
		0	Si IV 1402				45.142 eV
		1
		0	Si IV 1393				45.142 eV
		1
		0	C IV 1550				64.494 eV
		1
		0	C IV 1548				64.494 eV
		-240	-160	-80	0 v [km/s]	80	160	240
	Figure A1. APPENDIX A: IONIZATION STRUCTURE, FITTING
	PARAMETERS AND DETAILS ON H 2 J≥2 ROTATIONAL
	LEVELS

Table A1 .

 A1 Fitting results. Each section of Table delimited horizontal lines contains species that were fitted individually.

		Component A		Component B		Component C		total
	species	z	log N	b, km s -1	z	log N	b, km s -1	z	log N	b, km s -1	log N
	Cl i	2.786459(33) 13.21 +0.18 -0.11	7.9 +4.3 -1.3	2.786582(9) 13.35 +0.28 -0.16	2.5 +1.7 -1.2	2.786720(55) 12.62 +0.26 -0.30	3.9 +8.5 -2.1	13.63 +0.20 -0.05
	H i			One component fit with z = 2.786540(150)			21.82 ± 0.11
	H 2 , J=0			One component fit with z = 2.786570(14)			20.71 ± 0.02
	H 2 , J=1					"					21.06 ± 0.02
	H 2 , total										21.21 ± 0.02
	H, total										21.99 +0.08 -0.07
	H 2 , J=2	2.786459 a 19.47 +0.03 -0.05	8.8 +0.1 -0.3	2.786582 a 19.10 +0.04 -0.19	6.1 +0.2 -0.2	2.786720 a 16.00 +0.30 -1.70	9.4 +0.9 -0.4	19.59 +0.02 -0.02
	H 2 , J=3	"	17.67 +0.09 -0.20	"	"	18.71 +0.02 -0.03	"	"	13.29 +0.77 -0.14	"	18.74 +0.02 -0.03
	H 2 , J=4	"	16.44 +0.08 -0.05	"	"	17.14 +0.15 -0.13	"	"	14.80 +0.08 -0.05	"	17.23 +0.13 -0.13
	H 2 , J=5	"	15.89 +0.05 -0.05	"	"	15.85 +0.10 -0.05	"	"	15.30 +0.06 -0.04	"	16.23 +0.05 -0.04
	H 2 , J=6	"	14.76 +0.02 -0.03	"	"	13.80 +0.15 -0.38	"	"	14.45 +0.04 -0.05	"	14.95 +0.02 -0.02
	H 2 , J=7	"	14.57 +0.05 -0.02	"	"	13.85 +0.18 -0.41	"	"	14.54 +0.04 -0.05	"	14.89 +0.02 -0.02
	H 2 , J=8	"	14.04 +0.07 -0.07	"	"	13.28 +0.41 -0.49	"	"	13.88 +0.09 -0.18	"	14.28 +0.06 -0.07
	H 2 , J=9	"	13.85 +0.12 -0.09	"	"	13.70 +0.08 -0.25	"	"	13.38 +0.22 -0.20	"	14.14 +0.06 -0.07
	H 2 , J=10	"	< 13.3	"	"	< 13.4	"	"	< 12.9	"	< 13.6
	H 2 , J=0, ν = 1	"	< 13.7	"	"	< 13.4	"	"	< 13.4	"	< 13.9
	H 2 , J=1, ν = 1	"	< 13.7	"	"	< 13.6	"	"	< 13.5	"	< 13.9
	H 2 , J=2, ν = 1	"	< 13.7	"	"	< 13.5	"	"	< 13.7	"	< 14.0
	HD, J=0	2.786459 a 15.15 +0.27 -0.12	6.4 +0.5 -0.6	2.786582 a 17.34 +0.13 -0.38	1.5 +0.6 -0.4	2.786582 a	< 14.0		17.34 +0.13 -0.37
	HD, J=1	"	14.99 +0.03 -0.05	"	"	15.81 +0.75 -0.66	"	"	< 14.1		15.87 +0.72 -0.49
	HD, total										17.35 +0.15 -0.34
	CO	2.786459 a	< 13.1		2.786582 a	< 13.0					< 13.4
	C i	2.786435(9) 13.09 +0.05 -0.04	4.7 +1.2 -0.7	2.786578(4) 13.33 +0.05 -0.03	3.6 +0.4 -0.6				13.53 +0.04 -0.02
	C i *	"	13.13 +0.08 -0.05	"	"	13.44 +0.02 -0.05	"				13.61 +0.02 -0.02
	C i * *	'	12.90 +0.10 -0.10	"	"	13.14 +0.04 -0.04	"				13.34 +0.04 -0.03
	b c turb , km s -1			4.2 +1.1 -1.1			9.1 +0.6 -0.5			6.8 +2.8 -1.8
	T c , 10 4 K			0.8 +1.2 -0.6			1.1 +0.7 -0.7			1.2 +0.7 -0.7
	Mg ii	2.786377(6) 14.67 +0.48 -1.55		2.786555(4) 15.91 +0.03 -0.05		2.786775(3) 12.72 +1.91 -0.48		15.94 +0.06 -0.05
	Si ii	"	14.94 +0.07 -0.13		"	15.62 +0.04 -0.06		"	14.90 +0.14 -0.14		15.77 +0.04 -0.05
	S ii	"	13.67 +0.34 -0.13		"	15.51 +0.04 -0.02		"	14.59 +0.23 -0.12		15.57 +0.05 -0.02
	Ti ii	"	12.03 +0.18 -0.00		"	12.84 +0.07 -0.07		"	12.24 +0.24 -0.78		12.98 +0.07 -0.07
	Cr ii	"	12.52 +0.12 -0.16		"	13.08 +0.05 -0.05		"	12.24 +0.16 -0.94		13.23 +0.05 -0.06
	Mn ii	"	12.32 +0.26 -0.21		"	13.60 +0.07 -0.06		"	12.64 +0.21 -0.48		13.67 +0.07 -0.06
	Fe ii	"	14.27 +0.07 -0.12		"	14.80 +0.05 -0.04		"	14.04 +0.16 -0.32		14.96 +0.04 -0.04
	Ni ii	"	12.91 +0.11 -0.08		"	13.50 +0.04 -0.06		"	12.97 +0.15 -0.15		13.69 +0.05 -0.05
	Zn ii	"	12.12 +0.10 -0.22		"	12.91 +0.04 -0.05		"	12.12 +0.06 -0.31		13.03 +0.03 -0.06
	P ii										< 13.72
	Si ii*	2.786377 d 11.32 +0.14 -0.21		2.786555 d 11.62 +0.08 -0.14					11.80 +0.07 -0.11
	C ii*	2.786377 d	< 13.6		2.786555 d 15.55 +0.16 -0.18

http://www.ast.cam.ac.uk/∼rfc/vpfit.html

In many cases, the derived likelihood function is actually asymmetric, resulting in asymmetric error bars.MNRAS 000, 1-20 (2017)

MNRAS 000, 1-20 (2017)

The rest arising from the compact continuum emission close to the central black hole (hence very likely completely covered as also indicated by the saturated S ii λ1253 line, see Fig.7).

A few additional CO measurements are available(Noterdaeme et al. 2009(Noterdaeme et al. , 2011;;Ma et al. 2015) for quasar DLAs but their lower redshifts prohibit direct H 2 measurements from ground-based telescopes.MNRAS 000, 1-20 (2017)

Acknowledgments:

Based on observations collected at the European Organisation for Astronomical Research in the Southern Hemishpere under ESO programme 092.A.0345(A). We are very grateful to the anonymous referee for the detailed and careful reading our manuscript and many useful comments. This work was supported by the Russian Science Foundation grant 14-12-00955. SB thanks the French Government, Campus France and the Institut Lagrange de Paris for support during his visits at IAP. SB and RS thank IAP for hospitality and support during the time part of this work was done. PN, PPJ and RS gratefully acknowledge support from Indo-French Centre for the Promotion of Advanced Research under project No. 5504-B. PN acknowledges support from the 'Programme National de Cosmologie and Galaxies' (PNCG) funded by CNRS/INSU-IN2P3-INP, CEA and CNES, France.