Free Access
Issue
A&A
Volume 611, March 2018
Article Number A51
Number of page(s) 17
Section Interstellar and circumstellar matter
DOI https://doi.org/10.1051/0004-6361/201730797
Published online 26 March 2018

© ESO 2018

1 Introduction

Theoretical works on heating and cooling in the interstellar medium (ISM) predict the existence of two thermodynamically stable phases in the neutral atomic gas (Field et al. 1969; McKee & Ostriker 1977): the warm neutral medium (WNM) and the cold neutral medium (CNM). The volume density and kinetic temperature in the WNM are in the range 0.03–1.3 cm−3, and 4100–8800 K, respectively, while in the CNM their ranges are 5–120 cm−3, and 40–200 K, respectively (Wolfire et al. 2003). Given the densities of WNM and CNM, collisions between electrons, ions, and H atoms are able to thermalize the 21 cm transition in the CNM, but not in the WNM. So the spin temperature is expected to be close to the kinetic one in the CNM and less than the kinetic one in the WNM (Deguchi & Watson 1985; Liszt 2001). The 21 cm emission line of the atomic hydrogen traces the whole neutral atomic gas from WNM to CNM, but, without independentmeasurements of the 21 cm line in absorption toward distant radio sources to determine the spin temperature, one cannot retrieve the exact column density of hydrogen. Applying a uniform spin temperature, TS , over a cloudcomplex or along sight lines provides an estimate of the proportion of CNM and WNM, and only an average correction of the column densities in the dense CNM. The amplitude of column density correction with respect to the optically thin case increases with decreasing spin temperature and increasing AV . It reaches up to 35% at AV ~ 3 (Liszt 2014).

The molecular phase is mainly composed of molecular hydrogen and is commonly traced by the J = 1 → 0 line emission of 12 CO molecules at115 GHz, hereafter simply referred to as CO lines. In their study of CO photo-dissociation and chemistry, van Dishoeck & Black (1988) predicted that CO emission may not trace all the molecular gas in translucent clouds. In the diffuse envelopes near the H I–H2 transition, H2 molecules survive more efficiently than CO against dissociation by UV radiation (Wolfire et al. 2010). Furthermore, the brightness of the 12 CO line depends on the abundance and the level of collisional excitation of CO molecules, which are both low in the diffuse H2 envelopes. Large quantities ofH2 can therefore be CO-dark below the sensitivity level of the current large-scale CO surveys (of order 1 K km s−1).

Chemical species precursors to CO formation, such as C II, CH, and OH, are used to trace the gas at the atomic-to-molecular transition. Magnani et al. (2003) found that the CO line intensity in the MBM16 translucent cloud weakly correlates with the stellar reddening by dust, E(BV ), contrary to the CH line intensity at 3335 MHz. Similarly, [C II] line intensity at 158 microns and OH line intensity at 18 cm trace a larger fraction of molecular gas than CO line intensity in diffuse environments (Velusamy et al. 2010; Barriault et al. 2010). However, [C II] line emission also arises from the widely spread regions of atomic and ionized gas.

Blitz et al. (1990) have compiled a catalogue of clouds detected by dust emission in the far infrared (FIR) at 100 μm with IRAS, but not seen in 12CO observations. These clouds at low extinction (<0.25 mag) were proposed to be made of the diffuse molecular gas studied here. Röhser et al. (2014) have studied intermediate-velocity clouds (IVCs) with and without a far-infrared (FIR) excess compared to the H I column densities. They suggest that dim and bright IVCs in the FIR probe different stages of the transition from atomic to molecular phase.

Several studies have used the additional information provided by tracers of the total gas column density and have compared them to H I and 12CO line intensities in order to map the dark gas in the Galaxy at different scales. This has been done using γ-ray emission (Grenier et al. 2005), thermal emission from large dust grains (Planck Collaboration XIX 2011), and stellar reddening by dust (Paradis et al. 2012). In these studies, the H I spin temperature is not fully constrained for each line of sight, so the additional gas may include a fraction of the CNM in addition to diffuse molecular gas. The suggestion that optically thick H I fully or largely accounts for the dark gas (Fukui et al. 2015) is challenged by H I absorption, dust extinction, and γ-ray observations that we detail in the discussion part of the paper. The additional gas is thus likely to include both optically thick H I and CO-dark H2 at the transition between the atomic and molecular phases. In the absence of extensive emission tracers for this gas, we refer to this transition as the dark neutral medium (DNM).

Simulations indicate that the CO line brightness rapidly saturates at large visual extinctions, near 4 mag (Shetty et al. 2011). Similarly, observations show that CO-line-absorption measurements at high optical depth reveal features not seen in emission (Liszt & Pety 2012). In Remy et al. (2017) we detected additional gas seen with both dust and γ-ray emissions towards dense regions where the 12CO line emission saturates. We refer to this additional, CO-saturated H2 component as COsat.

In this paper we use the DNM and COsat maps derived jointly from γ-ray emission and dust optical depth at 353 GHz in addition to the H I and CO emission data in order to investigate the H I–H2 transition from diffuse to dense molecular clouds. We provide constraints on the range in N H tot $N_{\textrm{H}}^{\textrm{tot}}$ and AV required for H2 formation. We follow the evolution of the contribution of CO-dark H2 to the totalmolecular gas inside clouds of different types and test the range of H2 column densities where 12 CO line emission saturates. In order to compare clouds of different types, we have selected a set of local clouds spanning between 140 and 420 pc in distance and with masses between 600 and 34 400 solar masses. The sample includes the local anti-centre clouds of Cetus, Taurus, Auriga, Perseus, California analysed in Remy et al. (2017) and the Chamaeleon cloud analysed in Planck Collaboration Int. XXVIII (2015). These two analyses are based on the same gas tracers and use the same method.

The paper is organised as follows. Details on the method and on the cloud sample are given in the following section. The results and discussions are presented in Sect. 3. We focus on the transitions between the different gas phases in Sect. 3.1. The evolution of the relative contributions of the CO-dark, CO-bright, and CO-saturated H2 to the molecular phase across a cloud is discussed in Sect. 3.2. Cloud-to-cloud variations of the average mass fractions are discussed in Sect. 3.3. The main conclusions and possible follow-up studies are listed in the last section.

2 Analyses and data

The results presented here are based on the same dust and γ-ray analysis procedures as applied to two broad regions in the anticentre (Remy et al. 2017) and in the Chamaeleon directions (Planck Collaboration Int. XXVIII 2015) in the sky. The analysis method and multi-wavelength data used to build the γ-ray and dust models are detailed in these papers. We summarise their main features below.

To trace the total gas we have used the dust optical depth, τ353 , inferred at 353 GHz from the spectral energy distribution of the thermal emission of the large grains, which has been recorded between 353 and 3000 GHz by Planck and IRAS (Planck Collaboration XI 2014). We have also used six years of Pass 8 photon data from the Fermi Large Area Telescope (LAT) between 0.4 and 100 GeV. Technical details on the photon selection, instrument response functions, energy bands, and the ancillary γ-ray data (local interstellar spectrum, point sources, inverse Compton and isotropic components) that we have used are given in Remy et al. (2017).

In order to display the results as a function of the visual extinction AV in both regions, we have used the all-sky AJ map constructed by Juvela & Montillaud (2016) from the 2MASS extinction data, with the NICEST method at 12.0 arcmin resolution. AJ values are converted into AV with a factor 3.55 according to the extinction law of Cardelli et al. (1989). Other AV datasets are available, but they do not cover either the anticentre or Chamaeleon regions using the same extraction method. We have not modelled AV as a total gas tracer as we have done with γ rays and τ353 because of significant discrepancies in the AV values obtained with different methods in the diffuse parts of the clouds (AV ≲ 1 mag) that are the focus of this study (Rowles & Froebrich 2009; Green et al. 2015; Juvela & Montillaud 2016).

2.1 Summary of γ-ray and dust models

Galactic cosmic rays (CRs) interact with gas and low-energy radiation fields to produce γ rays. At the energies relevant for the LAT observations, the particle diffusion lengths in the ISM exceed the cloud dimensions, and there is no spectral indication ofvariations in CR flux inside the clouds studied (Planck Collaboration Int. XXVIII 2015; Remy et al. 2017). Therefore, the interstellar γ radiation can be modelled, to first order, as a linear combination of gas column densities in the various phases and the different regions seen along the lines of sight. The model also includes a contribution from the large-scale Galactic inverse-Compton (IC) emission, anisotropic intensity to account for the instrumental and extragalactic backgrounds, and point sources (see Eq. (7) of Remy et al. 2017).

The large dust grains responsible for the thermal emission seen in the far IR and at sub-mm wavelengths are supposed to be well mixed with the interstellar gas. For a uniform mass emission coefficient of the grains, the dust optical depth should approximatively scale with the total gas column densities. Therefore, we have also modelled the optical depth map as a linear combinationof gas column densities in the various phases and in the different cloud complexes (see Eq. (6) of Remy et al. 2017).

The γ-ray intensity and the dust optical depth have been jointly modelled as a linear combination of H I-bright, CO-bright, and ionized gas components (plus the non-gaseous ancillary components of the γ-ray model mentioned above). The coefficients of the linear combination are used to scale the dust emission and γ-ray emission of gaseous origin into gas column density, and to independently probe the different gas phases and the different clouds separated in position-velocity. We have used 21-cm H I and 2.6-mm CO emission lines to trace the atomic and CO-bright molecular hydrogen, respectively,and we have used 70-GHz free-free emission to trace the ionized gas in the anticentre region. The contribution of ionized gas in the Chamaeleon region is negligible (Planck Collaboration Int. XXVIII 2015). The only way the spatial distributions of the γ-ray intensity and of the dust optical depth can correlate is via the underlying gas structure (Grenier et al. 2005; Planck Collaboration Int. XXVIII 2015). So the joint information from dust emission and γ rays is used to reveal the gas not seen, or poorly traced, by H I, free-free, and 12CO emissions, namely the opaque H I and diffuse H2 present in the DNM, as well as the dense H2 to be added when the 12CO line emission saturates (COsat). The two templates of additional gas used in the γ-ray model are inferred from the dust model. In the anticentre analysis we have separated the DNM and COsat components in directions inside and outside of the 7 K km s−1 contour in WCO intensity (see Sect. 2.3 for details). Such a separation had not been done for the analysis of the Chamaleon region (Planck Collaboration Int. XXVIII 2015), so for the present study we have updated the analysis of the Chamaeleon region in order to separate the DNM and COsat components in terms of direction using the same cut.

The residual maps (data minus model) presented in Planck Collaboration Int. XXVIII (2015) and Remy et al. (2017) show that the best-fit models adequately describe both the Planck-IRAS observations in dust optical depth and the LAT observations in γ-ray intensity. Moreover, jackknife tests have shown that the fitted parameters of the γ-ray and dust models apply statistically well to the whole region, both in the anticentre and Chamaeleon directions.

The residuals obtained between γ-ray data and the best fits are consistent with noise at all angular scales and they show that the linear model provides an excellent fit to the γ-ray data in the overall energy band, as well as in the separate energy bands. Significant positive residuals remain at small angular scales in the dust fit. They are likely due to the limitations in angular resolution and in sensitivity of the γ-ray DNM template used in the dust model. Small-scale clumps in the residual structure can also reflect localised variations in dust properties per gas nucleon that are not accounted for in the linear models. These effects are discussed in Remy et al. (2017).

2.2 H I and CO data and cloud selection

H I and CO emission in the anticentre and Chamaeleon regions have been mapped by different surveys. Their relative calibrations have been verified. The Chamaeleon and anticentre analyses use the CO data from the Nanten and the Center for Astrophysics surveys (Mizuno et al. 2001; Dame et al. 2001; Dame & Thaddeus 2004), respectively. After correction, the integrated CO intensities, WCO, recorded in the Chamaeleon region by the two surveys are close to a one-to-one correlation (Planck Collaboration Int. XXVIII 2015). The H I survey GASS II (Kalberla et al. 2010) covers the Chamaeleon region. The GALFA and EBHIS surveys (Peek et al. 2011; Winkel et al. 2016) encompass the anticentre region. We have verified the good correlation between EHBIS and GALFAcolumn densities in the anticentre region (Remy et al. 2017). According to Kalberla & Haud (2015) and Winkel et al. (2016), the brightness temperature of GASS II is on average 4% higher than that of EBHIS. This discrepancy does not affect the column density profiles we present through single clouds (each covered by a single survey), but it induces a small shift between the gas fractions measured in the Chamaeleon region relative to the anticentre ones. Given the uncertainties in the γ-ray and dust models of the order of 5 to10%, a 4% shift does not significantly bias our results. We can therefore compare the H I column densities and CO integrated intensities measured inboth regions and merge the two samples of clouds to study cloud-to-cloud variations.

We have used H I and CO lines to kinematically separate cloud complexes along the lines of sight in both regions and to separate thenearby clouds from the Galactic backgrounds. Details of the pseudo-Voigt line decompositions and selection of longitude, latitude, and velocity boundaries of the complexes are given in Planck Collaboration Int. XXVIII (2015) and Remy et al. (2017). The analyses resulted in the separation of six different complexes in the anticentre region (Cetus, South Taurus, North Taurus, Main Taurus Perseus, and California; see Fig. 1 of Remy et al. 2017) and two complexes in the Chamaeleon region (Chamaeleon and IVA; see Fig. 2 of Planck Collaboration Int. XXVIII 2015). The complexes separated in position-velocity form coherent entities in H I and 12CO to which we can associate specific distances compiled from the literature (Planck Collaboration Int. XXVIII 2015; Remy et al. 2017).

We have selected 15 clouds or sub-regions within the broader complexes in order to study phase transitions in different entities. The contours and names of the sub-regions are given in Fig. 1. When referring to those clouds, we use the same names and colour scheme throughout the paper. The set of sub-regions in the Chamaeleon region follow that studied in Planck Collaboration Int. XXVIII (2015). The boundaries of the sub-regions have been chosen to avoid directions where the different complexes overlap along the lines of sight, so that we can study the relative contributions of the different gas phases to the total column density. This is particularly important for the DNM contributions which cannot be kinematically separated along the lines of sight. Inside a given cloud area, we associate all the DNM gas to the cloud. The morphological coherence between the H I, DNM, and CO structures within the area strengthens the DNM association. The areas around the CO clouds are wide enough to encompass the extent of their DNM phase. One cannot rule out the presence along the line-of-sight of an isolated translucent (DNM) cloud without CO counterpart in addition to the DNM envelopes surrounding the visible CO clouds. However, recent simulations (Valdivia et al. 2016) have shown that the different gas phases are well intermixed, so we expect most of the DNM gas to be spatially associated with the H I and CO phases present in the area.

We do not attempt to leave the parameters of the dust and γ-ray models free in the 15 smaller clouds because of the current γ-ray statistics, the separation in velocity of the different complexes, and the broad extension of the local H I structures. We obtain excellent γ-ray residuals with only six local complexes, so there is no obvious need for further subdivision. One should also note that the 15 sub-regions do not cover the whole analysis regions because we have excluded directions where clouds overlap in velocity/distance. It is not possible to model the 15 clouds without also modelling the remaining regions.

The 15 sub-regions are presented in Table 1, together with their masses inferred from the H I and CO data, their star-formation rates (SFR), and their XCO factors. All quoted masses are directly derived from the column density maps, taking into account the helium contribution. The H I masses have been calculated for optically thin conditions in the Chamaeleon region and for a spin temperature of 400 K in the anticentre region. Those conditions best match the interstellar γ-ray maps. We have applied a different CO-to-H2 conversion factor, X CO = N H 2 / W CO $X_{\textrm{CO}}= N_{\textrm{H}_2}/W_{\textrm{CO}}$, for each of their parent complexes (as we do not leave the model parameters free in the 15 clouds). The SFR have been derived from the mass of dense molecular gas seen in CO emission at AV > 8 mag according to the empirical relation expected from the “microphysics” of prestellar core formation within filaments (André et al. 2014; Könyves et al. 2015; Shimajiri et al. 2017): SFR = (4.5 ± 2.5) × 10−8 M yr M H 2 A V >8 $M^{A_{\textrm{V}}>8}_{{\textrm{H}_2}} $). This relation is in good agreement with measurements in Local clouds (Lada et al. 2010; Shimajiri et al. 2017). Table 1 shows that the dataset includes clouds ranging from low-column-density translucent clouds to star-forming molecular clouds belonging to the solar neighbourhood, with distances ranging between 140 and 420 pc and with a large span of masses between 200 and 18 300 solar masses in the CO-bright phase.

thumbnail Fig. 1

Integrated CO intensity, WCO, overlaid with the boundaries of the selected substructures in the anticentre (left) and Chamaeleon (right) regions.

Table 1

Distance to the clouds, mass in their H I and CO phases, SFR, and XCO factors.

2.3 Estimation of the gas column density

Several studies have shown that the dust opacity, τ353NH, rises as the gas, whether atomic or molecular, becomes denser. This opacity increase is likely due to a change in emission properties of the dust related to a chemical or structural evolution of the grains (Martin et al. 2012; Köhler et al. 2015). In the dense molecular filaments of the Taurus and Chamaeleon clouds, the opacity can rise by a factor of 2 to 4 above the value found in the diffuse ISM (Stepnik et al. 2003; Planck Collaboration XXV 2011; Ysard et al. 2013; Planck Collaboration Int. XXVIII 2015; Remy et al. 2017). Such a rise severely limits the use of the dust thermal emission to trace the total gas in dense media. The limit of the linear regime is estimated to be N H tot <2× 10 21 $N_{\textrm{H}}^{\textrm{tot}} <2\times10^{21}$ cm−2 in the Chamaeleon clouds and N H tot <3× 10 21 $N_{\textrm{H}}^{\textrm{tot}} <3\times10^{21}$ cm−2 in the anticentre clouds.

In all the nearby clouds studied so far and across all the gas phases we find no evidence of spectral evolution of the γ-ray emissivity of the gas relative to the local average. The latter is referred to as the local interstellar spectrum (LIS). This consistency is verified in the denser molecular parts (COsat). At the multi-GeV to TeV particle energies relevant for the LAT observations, the same CR flux pervades the clouds, from the diffuse atomic outskirts to the dense parts seen in 12CO line emission. Therefore the same CR flux pervades through the bulk of a cloud volume and mass. The γ-ray emission efficiently traces the total column density of gas, contrary to the dust optical depth which is affected by the evolution of the grain emission properties. For this reason we preferentially estimate N H tot $N_{\textrm{H}}^{\textrm{tot}}$ using γ rays.

The total column density is the sum of the H I, DNM, CO, and COsat contributions scaled according to the best-fit parameter of the γ-ray model (details for each components are given in the following paragraphs). The total molecular column density is the sum of the molecular DNM, CO, and COsat column densities derived from the γ-ray fit: N H 2 tot = N H 2 DNM + N H 2 CO + N H 2 COsat $N_{{\textrm{H}_2}}^{\textrm{tot}}=N_{{\textrm{H}_2}}^{\textrm{DNM}}+N_{{\textrm{H}_2}}^{\textrm{CO}}+N_{{\textrm{H}_2}}^{\textrm{COsat}}$. As the DNM gas can be composed of both opaque H I and diffuse H2, we have assumed a composition varying between 50 and 100% of molecular hydrogen (see discussion at the beginning of Sect. 3.2). By default, the results are given for a half molecular DNM composition, except if specified otherwise.

For the calculation of the H I column densities, NHI, we did not attempt to set a different H I spin temperature for each complex, but we have used a uniform average spin temperature. In the following the calculations are performed at the best-fit spin temperature of 400 K in the anticentre clouds, and for an optically thin emission in the Chamaeleon clouds (Planck Collaboration Int. XXVIII 2015; Remy et al. 2017).

The column density map of DNM component map derived from the dust distribution and τ353 analysis has been scaled in mass with the γ-ray emission under the assumption of a uniform CR flux permeating the H I and DNM phases. This assumption is substantiated by the uniformity of the CR spectrum found in the H I, DNM, and CO phases. By construction, the DNM maps spatially exclude the denser CO regions (WCO > 7 K km s−1) that are attributed to the COsat component. This separation that we have adopted between the DNM and COsat components is only a directional boundary. The 7 K km s−1 cut in WCO intensity is chosen deeply enough into CO clouds to ensure that the CO-bright H2 column density dominates the NH in the other phases. This choice is specific to the present set of clouds to ensure a uniform cut in all six local cloud complexes. Because of the complex three-dimensional (3D) structure of the clouds, the DNM can be detected along lines of sight intercepting CO emission in excess of 1 K km s−1, but the fraction of DNM gas decreases from the cloud outskirts to its dense, CO-intense parts. The choice of 7 K km s−1 ensures that the additional gas detected by the dust and γ rays towards the dense molecular regions is predominantly, in column density, in the COsat phase. We show in Sect. 3 that the DNM gas contributions towards these directions do not exceed 10% of the COsat column density.

The quantity of gas in the molecular phase is usually estimated via the 12CO line emission and scaled into hydrogen column density via the N H 2 / W CO $N_{\textrm{H}_2}/W_{\textrm{CO}}$ conversion factor, XCO. Thanks to the good penetration of CRs up to the CO-bright phase, the H I and CO related parameters of our γ-ray model allow the derivation ofmean XCO values for each cloud complex independently. In a previous paper (Remy et al. 2017), we demonstrated that these average XCO ratios significantly decrease from the diffuse clouds, subject to heavier photo-dissociation, to the more compact, dense, and well-shielded CO clouds. According to theory (Bell et al. 2006; Glover & Mac Low 2011), the differences in mean XCO values reflect intrinsic differences in CO abundance and excitation conditions, and therefore in XCO spatial gradients across the various clouds. This is why we have used the individual XCO values of Table 1 in our models. We have not attempted to model spatial gradients in XCO across the clouds since there is no consensus on the theoretical prescription for how XCO should scale with WCO (Bell et al. 2006; Liszt et al. 2010; Glover & Mac Low 2011; Shetty et al. 2011; Liszt & Pety 2012; Bertram et al. 2016).

3 Results and discussion

The results of the γ-ray and dust models obtained from the updated analysis of the Chameleon region compare very well with previously published studies. The XCO factor of (0.65 ± 0.02) × 1020 cm−2 K−1 km−1 s derived from the γ-ray fit has not changed. The new dust opacities found in the H I, DNM, and CO phases are τ 353 /N ¯ H HI =(15.0±0.2)× 10 27 $\overline{\tau_{353}/N}_{\textrm{H}}^{\textrm{HI}} = (15.0 \pm 0.2) \times 10^{-27}$ cm2 H−1, τ 353 /N ¯ H DNM =(17.9±0.5)× 10 27 $\overline{\tau_{353}/N}_{\textrm{H}}^{\textrm{DNM}} = (17.9 \pm 0.5) \times 10^{-27}$ cm2 H−1, and τ 353 /N ¯ H CO =(32.9±1.0)× 10 27 $\overline{\tau_{353}/N}_{\textrm{H}}^{\textrm{CO}} = (32.9 \pm 1.0) \times 10^{-27}$ cm2 H−1, respectively. Compared to values from the previously cited analyses the − 8%, − 4%, and + 3% changes, respectively, are not significant.

The hydrogen column-density maps obtained in the H I, DNM, CO, and COsat phases seen in the anticentre and Chamaeleon regions are shown in Figs. 2 and 3, respectively. The CO map is based on the WCO data and on the XCO ratios derived in γ rays for each individual cloud complex. In both regions, the column densities span from 1020 cm−2 in the atomic gas to a few times 1022 cm−2 in the molecular cores. We find comparable ranges of column densities in each phase of neutral gas. The means of the column density distributions in the anticentre clouds are larger than in the smaller Chamaeleon clouds, typically by a factor of two in the H I and DNM phases, and by a factor of three in the CO and COsat phases.

Figures 2 and 3 illustrate that both analyses exhibit significant quantities of gas not linearly traced by NHI column densities and WCO intensities, but jointly revealed by the dust mixed with gas and by the CRs spreading through it. The separation of the additional gas into DNM and COsat components according to the WCO intensity in each direction (for WCO > 7 K km s−1) highlights the structural differences between the spatially extended, diffuse DNM gas gathering at the H I and CO interface, and the more compact, denser filaments and clumps of molecular gas in regions of large WCO intensities. Figure 4 shows that the additional gas detected in the COsat component is likely due to large CO opacities because of the good spatial correspondence between the distributions of the COsat hydrogen column densities and of the intensities of the optically thinner 13CO J = 1 → 0 lines observed in the main Taurus and Perseus clouds (Narayanan et al. 2008; Ridge et al. 2006). Further exploiting the quantitative relation between the observed 13CO intensities, the saturating 12CO intensities, and the COsat column densities is beyond the scope of this paper. It requires careful modelling of the 12CO and 13CO radiative transfer and of the non-linear evolution of the dust emissivities that are used with γ rays to trace the additional gas (see the discussions of the evolution of τ353NH opacities at large N H tot $N_{\textrm{H}}^{\textrm{tot}}$ in Remy et al. 2017). We note that clouds with equivalent WCO intensities may have different quantities of additional H2 in the COsat component. For example, we observe comparable WCO intensities in the filaments of Musca and ChaEastI, and in the clumps of ChaI and ChaEastII, but we find hardly any COsat gas in ChaEastIand ChaEastII. Musca and ChaI, conversely, apparently require additional column densities as large as 1022 cm−2.

The use of a 13CO template covering the whole region would add a morphological constraint on the fit to trace the dense molecular gas not seen in 12CO emission. This would greatly help to separate COsat from DNM along the lines of sight. Moreover it would particularly improve the fit of the dust model in the dense H2 regions where the flaring of the τ353NH ratios or, alternatively, the loss of precise E(BV ) measurements, can significantly bias the estimation of the gas column density from the dust tracers. In the absence of an independent template for COsat, such as 13CO emission, our analysis cannot separate the additional COsat gas from the foreground and background DNM present along the lines of sight toward the dense regions of large WCO intensities. Therefore, the column densities of our COsat map include a small contribution from the DNM gas. The interpolation of the DNM H2 column densities across the COsat regions indicates that the COsat column density could be lowered by about 5–10% after subtraction of the DNM contributions.

DNM structures spatially extend between those of H I and CO clouds. The comparison between different clouds in the maps shows environmental or evolutionary differences in the DNM content of clouds. Thisis illustrated by differences among the high-latitude translucent molecular clouds. MBM18 (l = 189. °1 b = −36°), MBM16 (l = 170. °6 b = −37. °3), MBM12 (l = 159. °4 b = −34. °3), and the chain of CO clouds along Cetus (at b < −40°) exhibit a rich DNM component contrary to MBM8 (l = 151. °7 b = −38. °7) and MBM6 (l = 145. °1 b = −37. °3). The proportion of optically thick H I and CO-dark H2 in the DNM composition may also vary from cloud to cloud. An interesting structure is the DNM complex located at 170° < l < 190° and − 40° < b < −20°. This is one of the largest DNM structures identifiable in the early all-sky DNM maps (Grenier et al. 2005; Planck Collaboration XIX 2011; Paradis et al. 2012). It extends spatially well beyond the edges of CO emission with column densities 10 20 < N H DNM < 10 21 $10^{20}<N_{\textrm{H}}^{\textrm{DNM}}<10^{21}$ cm−2. Its morphology is comparable to that of the H I gas in South Taurus. Its structure contrasts with the DNM found close to the CO edges of the MBM clouds and the California cloud where the column densities can reach larger values near 2 × 1021 cm−2. We note that the fraction of the shell in California cloud missing in CO at l = 159. °5 and b = −10° is filled by DNM gas.

thumbnail Fig. 2

Hydrogen column density maps of the HI, DNM (derived from the τ353 and γ-ray fits), CO-bright, and COsat (derived from the τ353 and γ-ray fits) components in the Chamaeleon region. The XCO and τ353NH scale factors used in their derivation come from the γ-ray analysis for optically thin H I.

thumbnail Fig. 3

Hydrogen column density maps of the HI, DNM (derived from the τ353 and γ-ray fits), CO-bright, and COsat (derived from the τ353 and γ-ray fits) in the anticentre region. The XCO and τ353NH scale factors used in their derivation come from the γ-ray analysis with an H I spin temperature of 400 K.

thumbnail Fig. 4

Hydrogen column density map of the COsat component from the τ353 analysis overlaid with contours of 13CO line intensities at 2 and 4 K km s−1 in the Taurus and Perseus clouds (Narayanan et al. 2008; Ridge et al. 2006).

3.1 Phase transitions

In order to study the transition between the different neutral gas phases we have derived the ratios of the hydrogen column densities in the H I, DNM, CO-bright, and COsat components over the total column density. The ratios do not depend on cloud distances. We have plotted the ratios in the four neutral gas phases as a function of the total column density N H tot $N_{\textrm{H}}^{\textrm{tot}}$ and of the visual extinction AV in Figs. 5 and 6, respectively. Because of systematic uncertainties in the current measurements of low extinctions, we do not attempt to study the phase transitions below 0.3 mag in AV .

We find that in all fifteen clouds, the H I fraction in the total column density decreases for column densities greater than (0.6–1) × 1021 cm−2, in agreement with the theoretical limit of 1021 cm−2 necessary to shield H2 against UV dissociation for clouds of solar metallicity according to the models of Krumholz et al. (2009) and Liszt (2014). This limit is also confirmed by FUSE observations of H2 (Gillmon et al. 2006), by OH observations (Barriault et al. 2010), and by dust observations in Perseus (Lee et al. 2012), which measure NHI column densities leveling off at (0.3–0.5), (0.4–0.5), (0.8–1.4) × 1021 cm−2, respectively.We note that, even though the H I mass of the individual clouds spans a tenfold range, the H I fractional decrease occurs over a small range in N H tot $N_{\textrm{H}}^{\textrm{tot}}$, only a factor of four difference from cloud to cloud. We also note that the slope of the decline is approximately the same in the different clouds. The column density fractions in the H I often remain large (about 30%) towards the CO-bright phase because of the turbulent 3D structure of the clouds (Valdivia et al. 2016) and because of the atomic envelopes surrounding the molecular clouds.

The CO-dark and CO-bright regimes in clouds are expected to be offset in depth because H2 is better shelf-shielded against UV dissociation than CO at low extinction. The combination of dust screening and gas shielding is required at large extinction to protect CO molecules (Visser et al. 2009). The profiles of Fig. 6 show transitions to CO-bright emission that tend to be deeper into the clouds than for the H I–H2 transition. At NH < 1021 cm−2 (AV < 1 mag), the decrease in H I column-density fraction is larger than the increase in CO fraction. In this range the DNM fractions exceed the CO ones.

The DNM fraction in the total column density often peaks in the (1–3) × 1021 cm−2 range. Compared to the higher degree of organisation of the H I and CO-bright transitions, the DNM fractions show a large diversity of profiles that may be due to the large spatial overlap between the HI-bright, DNM, and CO-bright phases when integrated along the lines of sight.

We note that the DNM column densities are not very sensitive to the sensitivity threshold in CO observations for two reasons. Firstly, optically thick H I contributes with diffuse CO-dark H2 to the DNM column density. The H I proportion is unknown, but likely not zero. Secondly, a large fraction of the DNM map would have easily been detected in 12CO emission if the DNM gas were mostly H2 with a CO abundance and level of excitation equivalent to the conditions prevailing in the CO-bright parts. In order to check this assertion, we have applied an XCO ratio of 1020 cm−2 K−1 km−1 s, which characterises the diffuse CO clouds in our sample, to translate the DNM column densities to equivalent WCO intensities. We have then calculated the fraction of the total DNM solid angle where the equivalent WCO intensity exceeds 1 K km s−1 level that would have been easily detected. We find a fraction of 58% in solid angle, which gathers 74% of the DNM mass. In those directions the lack of CO emission is not due to the sensitivity limit of the CO survey, but to the genuine under-abundance and/or under-excitation of CO molecules.

The CO sensitivity threshold of the present observations are respectively 0.4 and 1.5 K km s−1 for the anticentre and Chamaeleon regions (Dame et al. 2001; Planck Collaboration Int. XXVIII 2015). At that sensitivity level, the onset of CO intensities occurs for total gas column densities ranging between 0.6 and 2.5 × 1021 cm−2, in close correspondence with the transitional drop in H I. We note again a small (factor-of-four) variation in the CO-bright transition from cloud to cloud.

The CO-bright transition occurs at total N H 2 $N_{\textrm{H}_2}$ column densities of (1.5–4) × 1020 cm−2. Simulations of photo-dissociation regions (PDR) by Levrier et al. (2012) show that the column density of CO molecules rises steeply after the molecular transition, at hydrogen volume densities of about 100 cm−3 and at N H 2 =(3 $N_{\textrm{H}_2} = (3$–6) × 1020 cm−2 depending on the density profile in the PDR region. A uniform density profile shifts the transition to larger N H 2 $N_{\textrm{H}_2}$. The value of 3 × 1020 cm−2 corresponding to the non-uniform case agrees reasonably well with our results, as well as with CO and H2 observations by the HST and FUSE, which indicate a transition at N H 2 =2.5× 10 20 $N_{\textrm{H}_2} = 2.5 \times 10^{20}$ cm−2 (Sheffer et al. 2008).

The saturation of the CO emission lines corresponds to a fractional increase of the column densities in the COsat component, but they contribute only 10% to 30% more gas. The saturation occurs in the main molecular clouds beyond (1.5 − 2.5) × 1021 cm−2 in N H tot $N_{\textrm{H}}^{\textrm{tot}}$.

thumbnail Fig. 5

Evolution with N H tot $N_{\textrm{H}}^{\textrm{tot}}$ of the fractions of the total hydrogen column density in the H I (top left), DNM (top right), CO (bottom left), and COsat (bottom right) components. The number of points is the same for all the curves. The points are equally spaced between the minimum and maximum of each profile. The colours refer to the same clouds as in Fig. 1.

thumbnail Fig. 6

Evolution with AV of the fractions of the total hydrogen column density in the H I (top left), DNM (top right), CO (bottom left), and COsat (bottom right) components. The colours refer to the same clouds as in Fig. 1.

3.2 Evolution of the CO-dark H2 fraction in column density within the clouds

This section focuses on the transitions inside the molecular phase from the CO-dark to CO-bright and then CO-saturated regimes in each of the clouds in our sample. We distinguish the diffuse CO-dark H2 present in the DNM from the dense H2 implied by the COsat component to the bright CO cores in order to isolate the fraction of dark molecular gas lying at the diffuse atomic-molecular interface.

We first discuss the transition from CO-dark to CO-bright H2 by measuring the change in the fractions of molecular DNM in the total H2 column density, N H 2 $N_{\textrm{H}_2}$. We define the CO-dark H2 fraction in column density as: \beginalign* f COdk H 2 = N H 2 DNM N H 2 tot = N H 2 DNM N H 2 DNM + N H 2 CO + N H 2 COsat nn\endalign* \begin{align&#x002A;} f_{{\textrm{COdk \, H}_2}} = \frac{N_{{\textrm{H}_2}}^{\textrm{DNM}}}{N_{{\textrm{H}_2}}^{\textrm{tot}}}=\frac{N_{{\textrm{H}_2}}^{\textrm{DNM}}}{N_{{\textrm{H}_2}}^{\textrm{DNM}}+N_{{\textrm{H}_2}}^{\textrm{CO}}+N_{{\textrm{H}_2}}^{\textrm{COsat}}}\cdot \end{align&#x002A;}(1)

Instead of the generic “dark gas fraction” (fDG) term found in the literature for column density or mass ratios and which often implicitly assumes that the whole DNM is molecular, we use the more explicit designation of “CO-dark H2 fraction” to restrict the fraction to the molecular part of the DNM (excluding the opaque H I) and to exclude issues with saturation of the WCO intensities (COsat components).

As the DNM lies at the interface between the atomic and molecular phases, it is expected to be composed of optically thick H I and diffuse H2 (not seen in CO at the level of sensitivity used here), both mixed with a proportion varying as the gas becomes denser and according to the ambient conditions (e.g. H2 volume density and kinetic temperature, intensity of the interstellar radiation field, ISRF, dust content, flux of ionizing low-energy CRs, Visser et al. 2009). The suggestion that optically thick H I fully or largely accounts for the additional gas in the DNM (Fukui et al. 2015) is challenged by several observations. Fukui et al. (2015) estimated a spin temperature of TS = 20−40 K and an optical depth τmax > 0.5 for 85% of their sky coverage (|b| > 15°). Stanimirović et al. (2014) measured H I absorption against 26 radio continuum sources toward the Perseus cloud and found that 54% of directions have τmax > 0.5 and only 15% of lines of sight have a spin temperature lower than 40 K. Their observation suggests that H I absorption traces mostly the central cloud regions where CO is bright, and traces the CO-dark envelopes only to a smaller degree, where DNM gas is more abundant (see Figs. 2 and 3). Corrections of the H I column densities are too small to fully explain the changes observed in the NHIE(BV ) ratios at E(BV ) > 0.1, so larger column density must be attributed to the onset of H2 formation (Liszt 2014). Analysis of the dust thermal emission as a function of Galactic latitude and of H I spin temperaturesuggests that the DNM is less than 50% atomic (Planck Collaboration XIX 2011). Moreover, doubling the H I densities in the local ISM because of ubiquitously large H I opacities would yield a twice lower γ-ray emissivity per gas nucleon, therefore a twice lower CR flux inferred from γ rays in the local ISM. This is clearly at variance with the direct CR measurements performed in the solar system at energies of 0.1 to 1 TeV for which solar modulation is negligible (Grenier et al. 2015). Large atomic fractions in the DNM probably being the exception for our set of clouds, we have performed the calculations with a 50% or 100% molecular composition of the DNM.

Figure 7 shows the evolution profile of the CO-dark H2 fraction in column density, for a half-molecular DNM composition, as a function of the total gas column density, N H tot $N_{\textrm{H}}^{\textrm{tot}}$, and of visual extinction, AV, in each of the clouds. The CO-dark diffuse molecular hydrogen exceeds 70% of the total H2 column density at N H tot <6× 10 20 $N_{\textrm{H}}^{\textrm{tot}} < 6\times 10^{20}$ cm−2, except in the tauM2 region which is centred on the extensively bright part of the Taurus CO cloud. Velusamy et al. (2010) observed 53 transition clouds with H I and 12CO emissions and no detectable 13CO emission, thus containing negligible COsat H2 . In about 30% of them, the [C II] line emission at 1.9 THz indicates CO-dark H2 fractions in column density exceeding 60%.

The fraction starts to decrease in the 0.4–0.9 mag range in visual extinction or in the (0.3–1) × 1021 cm−2 range in total gas column density. It then declines with rather comparable slopes in the different clouds. We note that CO emission in the very compact molecular filament of ChaEastI is detected at significantly lower column densities than in the other clouds, possibly because of larger H2 volume densities, but one cannot relate filamentary CO structures to low DNM abundances since the other compact filament in our sample, Musca, retains large DNM column densities up to 2 mag of extinction.

By integrating the gas over the entire anticentre or Chamaeleon regions, we find the average profiles displayed in Fig. 8 for a fully molecular composition of the DNM to allow comparisons with Xu et al. (2016). Systematic uncertainties in the current measurements of low extinctions limit our study of the CO-dark H2 fraction to AV > 0.3 mag. The rms dispersion reflects the differences in individual profiles of Fig. 7, due on the one hand to a modest difference, within a factor of 2 to 3, in the amount of extinction required to overcome CO destruction in each cloud and, on the other hand, to the difference in DNM abundance prior to the transition from CO-dark to CO-bright H2 . On average, the CO-dark H2 dominates the column densities up to 1.5 mag in visual extinction.

The individual and average f COdk H 2 $f_{{\textrm{COdk \, H}_2}}$ profiles do not compare well with the model prediction of Wolfire et al. (2010) that CO emission reaches an optical depth of 1 at AV ≳ 0.7 mag, well beyond the 0.1 mag required for the formation of H2 . Most of the profiles shown in Figs. 6 and 7 indicate a transition to the CO-bright regime deeper into the clouds, at AV between 1 and 2.5 mag, even though Wolfire et al. (2010) have modelled more massive clouds, with a total column density N H tot $N_{\textrm{H}}^{\textrm{tot}}$ of 1.6 × 1022 cm−2, that should screen CO molecules more effectively than the more translucent clouds in our sample. The onset of CO emission is more consistent with the approximate range near 1.5–2.5 mag where the dust optical depth and dust extinction are seen to recover a linear correlation with H I and CO intensities (Planck Collaboration XIX 2011; Paradis et al. 2012). We note that the cloud closest to the model prediction is the very compact filament of ChaEastI, suggesting that the volume density gradient through the cloud plays an important role in addition to the integral screening through the gas columns.

We parametrize the decrease of the CO-dark H2 fraction as afunction of AV with a half-Gaussian profile defined as: f COdk H 2 ={ f max if A V A V 0 f max ×exp[ ( A V A V 0 σ A V ) 2 ] if A V > A V 0 . $ f_{{\textrm{COdk \, H}_2}}= \left\lbrace \begin{array}{ccc} f_{\textrm{max}} & \mbox{if} & \textit{A}_{\textrm{V}} \leqslant \textit{A}_{{\textrm{V}_0}} \\ f_{\textrm{max}} \times \textrm{exp}\left[ - \left( \frac{\textit{A}_{\textrm{V}}-\textit{A}_{{\textrm{V}_0}}}{\sigma_{\textit{A}_{\textrm{V}}}} \right) ^2\right] & \mbox{if} & \textit{A}_{\textrm{V}} > \textit{A}_{{\textrm{V}_0}}.\\ \end{array}\right .$(2)

The results of the fits for the two analysis regions are given by the red curves in Fig. 8. We find fmax = 0.91, A V 0 =0.55 $A_{{\textrm{V}_0}}=0.55$ mag, σ A V =1.27 $\sigma_{A_{\textrm{V}}}=1.27$ mag in the anticentre region and fmax = 0.98, A V 0 =0.43 $A_{{\textrm{V}_0}}=0.43$ mag, σ A V =0.97 $\sigma_{A_{\textrm{V}}}=0.97$ mag in the Chameleon region. These parameters are derived for a fully molecular DNM, but they change only by a few percent when using a half molecular DNM composition.

Xu et al. (2016) have studied the PDR border of the main Taurus CO cloud in an approximately 0. ° 5-wide square, centred on l = 172. °8 and b = −13. °2, at a high angular resolution of 0.05°. This edge lies within the tauM2 sub-region of our sample. They observed OH, 12CO, and 13CO line emissions and inferred the N H 2 $N_{\textrm{H}_2}$ column densities from the AV extinction observations of Pineda et al. (2010), assuming that the hydrogen is predominately molecular and using the relation N H 2 / A V =9.4× 10 20 $N_{\textrm{H}_2}/A_{\textrm{V}} = 9.4 \times 10^{20}$ cm−2 mag−1. Their CO-dark H2 fraction is defined as the fraction of molecular gas not traced by 12CO or 13 CO line emission, so it compares with our measurements of f COdk H 2 $f_{{\textrm{COdk \, H}_2}}$ in the case of a fully molecular DNM composition. They modelled the variation of the CO-dark H2 fraction as a function of extinction with a Gaussian profile. Given the very different origins of both datasets, the close-up measurement of Xu et al. (2016) agrees remarkably well with the mean evolution of the CO-dark H2 fraction we find in the whole anticentre region. The steeper decline of their CO-dark H2 fraction compared to our mean trend remains compatible with the cloud-to-cloud dispersion in the sample. The discrepancy with the tauM2 profile (Fig. 7), however, calls for further investigation in order to understand whether it is due to methodological differences (angular resolution, total gas tracers versus radiative modelling, etc.) or to genuine spatial changes in the PDR profiles across tauM2.

Based on observations of [C II] lines at 1.9 THz with Herschel, Langer et al. (2014) estimated the column densities of CO-dark H2 in a large sample of clouds in the Galactic disc and measured the CO-dark H2 fraction in the total gas column density. As their estimates of the associated NHI column densities are very low, their CO-dark H2 fractions in the total gasare close to the f COdk H 2 $f_{{\textrm{COdk \, H}_2}}$ fraction we have defined in the molecular phase. They found a very broad range of CO-dark H2 fractions across the whole 1021−22 cm−2 interval. The spread largely encompasses the profiles shown in Fig. 7 for N H tot > 10 21 $N_{\textrm{H}}^{\textrm{tot}} > 10^{21}$ cm−2. We note, however, that the transitions we observe in the nearby clouds tend to occur at lower N H tot $N_{\textrm{H}}^{\textrm{tot}}$ column densities than the bulk of the Galactic sample probed with [C II] lines. They also occur outside the minimum and maximum fraction bounds predicted by Visser et al. (2009) and shown in Fig. 17 of Langer et al. (2014). Our observations show marked fraction declines at two to three times lower N H tot $N_{\textrm{H}}^{\textrm{tot}}$ column densities than the minimum model prediction, even though the nearby clouds probed in our study have the same solar metallicity and total extinction range (AV < 4 mag) as the gas slabs modelled by Visser et al. (2009). Radiative modelling of those slabs is required to compare the model profiles with our data and with the recorded WCO intensities in order to disentangle the chemical or radiative origin of this discrepancy.

Figure 9 illustrates the relationship between the CO-dark H2 fraction and the structure of the molecular clouds. Generally f COdk H 2 $f_{{\textrm{COdk \, H}_2}}$ progressively decreases from the DNM outskirts to the CO-bright edges, but local fluctuations imply a more complex pattern that is difficult to interpret in the 2D projection. Extended or elongated DNM-rich structures are more susceptible to photo-dissociation along the short axes of the cloud. The transition from the atomic to molecular phase is expected to occur at different column densities according to the strength of the ambient ISRF, outside and inside the cloud, according to the 3D spatial distribution of the gas and shadowing effects, and according to the 3D structure in density which impacts the ambient H2 and CO formation rates. Theoretical predictions by Krumholz et al. (2009) indicate that the fraction of molecular gas in a galaxy depends primarily on its column density and is, to a good approximation, independent of the strength of the interstellar radiation field. According to the spherical PDR model of Wolfire et al. (2010), an increased ISRF intensity implies a modest diminution in the thickness of the diffuse H2 envelope surrounding the CO cloud, but the ratio of H2 and CO depths remains constant, so the CO-dark H2 fraction remains stable against changes in ISRF intensity. Yet, this stability depends on the choice of volume density profile inside the cloud (nH ∝ 1∕r). For this profile they find only 30% changes in the CO-dark H2 fraction for UV radiation fields ranging between 3 and 30 times larger than the local value of Draine (1978).

We do not have a direct observational measure of the strength of the ISRF, but the specific power radiated by large dust grains per gas nucleon provides an indirect estimate of their heating rate by stellar radiation if the grains have reached thermal equilibrium. We used the radiance map of Planck Collaboration XI (2014) and our distributions of N H tot $N_{\textrm{H}}^{\textrm{tot}}$ column densities to map the specific power of the grains across the anticentre and Chamaeleon regions. We find only a 20% dispersion around the mean value of the specific power in the 0.4 < AV < 1.5 mag range that iscritical for the H I–H2 transition. We find no evidence of a spatial correlation between the distributions of the grain specific powers and those of the CO-dark H2 fractions shown in Fig. 9. As the strength of the ISRF potentially plays a marginal role in the relative contributions of CO-dark and CO-bright H2, the density structure of a cloud in 3D and the internal level of turbulence may be the key parameters to explain the small spread (by a factor of 4) in N H tot $N_{\textrm{H}}^{\textrm{tot}}$ and AV thresholds required to overcome CO destruction and sustain a population of CO molecules dense enough to be detected by the current CO surveys. Simulation of H2 formation within dynamically evolving turbulent molecular clouds shows that H2 is rapidly formed in dense clumps and dispersed in low-density and warm media due to the complex structure of molecular clouds (Valdivia et al. 2016).

We now turn to the transition from CO-bright to CO-saturated H2 that occurs in the cores of CO clouds when the optical thickness to 12CO emission is large enough. Figure 5 indicates that the CO-bright part of the total hydrogen column density saturates at AV ~ 2 mag whereas the COsat fraction still increases. We follow the variation of the COsat fractions in the total H2 column-density in Fig. 8. They reach 40 % at AV > 4 mag and 50% at AV > 6 mag. As mentioned in Sect. 3, our analysis cannot separate the additional COsat gas from the foreground and background DNM present along the lines of sight towards the dense regions of large WCO intensities, that is, large AV extinctions. The interpolation of the DNM H2 column densities across the COsat regions indicates that the COsat H2 fractions could be lowered by about 5–10% after subtraction of the DNM contributions.

thumbnail Fig. 7

Evolution of the fraction f COdk H 2 $f_{{\textrm{COdk \, H}_2}}$ of molecular DNM in the total H2 column-density as a function of the total gas column density N H tot $N_{\textrm{H}}^{\textrm{tot}}$ (top) and visual extinction AV (bottom). The DNM gas is assumed to be 50% molecular. The colours refer to the same clouds as in Fig. 1.

thumbnail Fig. 8

Evolution with visual extinction AV of the fraction f COdk H 2 $f_{{\textrm{COdk \, H}_2}}$ of molecular DNM gas (grey) or additional COsat H2 (green) in the total H2 column-density. The DNM is assumed to be 100% molecular. The shaded areas and error bars respectively give the standard deviation in the sample and the standard error of the mean. The black line gives the Gaussian profile of Xu et al. (2016) and the red line shows the best fit by a half-Gaussian to the f COdk H 2 $f_{{\textrm{COdk \, H}_2}}$ fraction.

thumbnail Fig. 9

Fraction of the molecular DNM in the total H2 column-density ( f COdk H 2 $f_{{\textrm{COdk \, H}_2}}$) derived fromthe dust (left) and γ-ray (right) analyses for the anticentre (top) and Chamaeleon (bottom) regions. The DNM is assumed to be 50% molecular. The white contours outline the shapes of the CO clouds at the 7 K km s−1 level chosen to separate DNM and COsat components.

3.3 Average CO-dark H2 fractions in mass and properties of the clouds

This section focuses on the variations of the average CO-dark H2 fractions in mass from cloud to cloud. We integrated the gas masses over the extent of each cloud in each of the H I, DNM, CO, and COsat phases. The masses are directly derived from the column density maps. All quoted masses include the Helium contribution. We used the column density maps derived from the γ-ray analyses, which are more robust against dust evolution (see details of their derivation in Sect. 2.3). The distance used for each cloud is the same as that of its respective parent complex given in Table 1. The mass ratios are not affected by distance uncertainties. Table 2 lists the relative contributions of the different phases to the total gas mass of clouds. The number distribution for each gas phase is shown in Fig. 10. We have also calculated the CO-dark H2 fraction in mass, integrated over each cloud, as: F COdk H 2 = M H 2 DNM M H 2 tot = M H 2 DNM M H 2 DNM + M H 2 CO + M H 2 COsat $ F_{{\textrm{COdk \, H}_2}}= \frac{M_{{\textrm{H}_2}}^{\textrm{DNM}}}{M_{{\textrm{H}_2}}^{\textrm{tot}}}=\frac{M_{{\textrm{H}_2}}^{\textrm{DNM}}}{M_{{\textrm{H}_2}}^{\textrm{DNM}}+M_{{\textrm{H}_2}}^{\textrm{CO}}+M_{{\textrm{H}_2}}^{\textrm{COsat}}}\cdot $(3) Table 2 provides the results for 50% and 100% molecular DNM composition.

Figure 10 highlights that atomic hydrogen is the dominant form of gas by mass in all the sampledclouds. The DNM is the second most important contributor, with mass fractions often close to 20%, in close agreement with the average of 25% found in the sample of 53 clouds exhibiting CO-dark H2 excesses in [C II] line emission (Velusamy et al. 2010). Bright CO clouds such as the star-forming tauM2 region or the compact ChaEastI filament gather only ~10% of their mass at the DNM interface. On the contrary, the fraction increases up to 27% in faint CO clouds such as tauS1 and tauS2. The values in Table 2 show that the relative amount of molecular gas present in the CO-dark and CO-bright regions varies greatly from cloud to cloud. Generally, the bright CO clouds, with a mean WCO intensity above 4 K km s−1 or a maximum intensity above 20 K km s−1, have more mass in the CO-bright region than in the DNM envelopes, contrary to the CO-faint and rather diffuse clouds. This is reflected by the average CO-dark H2 fractions varying from 10% to 50% in the brighter molecular clouds (tauM, Aur, Cal, Cha) and between 50% and 90% in the diffuse clouds (Cet, tauS, tauN).

Chen et al. (2015) studied a wider region of the Galactic anticentre that includes the clouds studied here and the Orion clouds, at an angular resolution of 0. °1. They modelled the stellar extinction data as a linear combination of H I column densities from the GALFA survey and CO intensities from the Planck type 3 map, corrected for 13CO contamination by multiplying the map by 0.86 (Planck Collaboration XIII 2014). Chen et al. (2015) used the differences between the data and model to derive the column density of DNM gas and then calculate the CO-dark H2 mass fraction. They found this mass fraction to be 24% in Taurus , 31% in Orion, and 47% in Perseus. For the high-latitude clouds TauE2 and TauE3, which approximately correspond to the tauS2 and tauS1 clouds of our sample, they found values of 80% and 77%. For a 50–100% molecular composition of the DNM, we find CO-dark H2 mass fractions of 11–19% in Taurus(tauM2), 0.82–0.90% in tauS2, and 0.79–0.88% in tauS1, in close agreement with the results of Chen et al. (2015).

We further discuss the relative contribution of the diffuse CO-dark H2 to the total molecular mass according to the cloud properties and location. Figure 11 shows the evolution of the CO-dark H2 mass fraction, F COdk H 2 $F_{{\textrm{COdk \, H}_2}}$, as a function of:

  • the surface fraction in solid angle of dense regions with large WCO intensity inside a cloud, SF dense = S( W CO>5 Kkms 1 ) S( W CO>1 Kkms 1 ) $\textrm{SF}_{\textrm{dense}}=\frac{\textit{S}(W_{\textrm{CO}>5 {\textrm{K km s}^{-1}}})}{\textit{S}(W_{\textrm{CO}>1 {\textrm{K km s}^{-1}}})}$; this parameter reflects the compactness of the molecular gas;

  • the maximum CO intensity, W CO max $W_{\textrm{CO}}^{\textrm{max}}$, recordedin the cloud;

  • the average visual extinction, A V ¯ $\overline{A_{\textrm{V}}}$, of the cloud;

  • the H2 mass found in the CO-bright phase (excluding the COsat component), M H 2 CO $M_{{\textrm{H}_2}}^{\textrm{CO}}$;

  • the cloud height above the Galactic plane.

We find a significant trend for F COdk H 2 $F_{{\textrm{COdk \, H}_2}}$ to decrease from 0.8 to 0.1 with increasing SFdense and W CO max $W_{\textrm{CO}}^{\textrm{max}}$. A negative slope is detected with a confidence level larger than 18σ for these diagnostics, indicating that the fraction tends to decrease from diffuse to more compact clouds. Figure 11 shows that the diffuse clouds, rich in CO-dark H2, also tend to lie at larger heights above the Galactic plane than the brighter and more compact CO clouds. There is a possible bias due to velocity crowding and an increased difficulty in separating the different complexes and gas phases in the dust and γ-ray models as the latitude decreases. We note, however, that the DNM-rich cloud of Musca lies close to the plane. Conversely, the medium-latitude clouds of TauM2, ChaII-III, and ChaEastI, well separated from the Galactic disc background, are DNM poor.

Using the thermal emission of dust grains over the whole sky at |b| > 5° and masking out CO regions with WCO intensities larger than 1 K km s−1, Planck Collaboration XIX (2011) have found an average F COdk H 2 $F_{{\textrm{COdk \, H}_2}}$ mass fraction of 0.54 close to the value of 0.62 obtained at |b| > 10° from dust visual extinction (Paradis et al. 2012). These high averages indeed reflect the local predominance of molecular clouds with moderate masses and diffuse structures in the sample of nearby clouds seen at medium latitudes. Mizuno et al. (2016) studied several nearby molecular clouds at intermediate latitudes, namely MBM 53, 54, and 55 and the Pegasus loop. They used the dust optical depth and the dust radiance, derived from Planck and IRAS data, in combination with γ-ray observations from Fermi LAT in order to monitor the variations of the dust emission properties as a function of the dust temperature, and to estimate the total gas column density. They found that the gas mass in the DNM not traced by H I and CO surveys is up to five times larger than the molecular gas mass traced by CO emission. Therefore, according to the definition in Eq. (3), their CO-dark H2 fraction reaches up to 75–83% for a 50–100% molecular DNM composition. Although their method differs from ours, this result agrees with the fraction we found in similar translucent clouds of the anticentre region (Cet, tauN, tauS in Table 2).

While it is clear in Fig. 11 that obscured clouds with A V ¯ $\overline{A_{\textrm{V}}}$ in excess of 1.2 mag have low CO-dark H2 mass fractions, below 20%, the large scatter at low A V ¯ $\overline{A_{\textrm{V}}}$ indicates that this observable cannot be used to decide whether the CO-dark H2 dominates over the CO-bright component in the molecular census. For this purpose, the peak CO intensity in a cloud may turn out to be more useful. Adding other clouds to the sample will confirm if clouds with W CO max $W_{\textrm{CO}}^{\textrm{max}}$ less than approximately 20 K km s−1 have more molecular hydrogen in the DNM than in the CO-bright parts. The model of Wolfire et al. (2010) predicts that, for a given cloud, F COdk H 2 $F_{{\textrm{COdk \, H}_2}}$ is a strong decreasing function of the mean extinction of the cloud, A V ¯ $\overline{A_{\textrm{V}}}$. Yet, their modelled fractions remain very large, near 75%, at A V ¯ =2 $\overline{A_{\textrm{V}}} = 2$ mag, at variance with the low fractions we find already above ~1.2 mag.

Except for the outlier CO filament of ChaEastI, compact and filamentary, which has particularly low mass, we find a significant decrease of the F COdk H 2 $F_{{\textrm{COdk \, H}_2}}$ fraction with the logarithm of the M H 2 CO $M_{{\textrm{H}_2}}^{\textrm{CO}}$ mass. The negative slope is preferred relative to a constant at a confidence level exceeding 19σ. Figure 12 further shows that the mass present in the DNM envelopes approximately scales with the square root of the molecular mass seen in CO. We find a correlation coefficient of 0.86 between the two sets of masses on a logarithmic scale. A power-law fit yields: M H DNM =(62±7) M H 2 CO 0.51±0.02 , $ M_{\textrm{H}}^{\textrm{DNM}} = (62 \pm 7) {M_{\textrm{H_2}}^{\textrm{CO}}}^{\,0.51 \pm 0.02},$(4)

if we include all the clouds, and M H DNM =(166±19) M H 2 CO 0.38±0.02 , $ M_{\textrm{H}}^{\textrm{DNM}} = (166 \pm 19) {M_{\textrm{H_2}}^{\textrm{CO}}}^{\,0.38 \pm 0.02} ,$(5)

if we exclude the outlier filament of ChaEastI.

The slopes are detected at the 13σ and 24σ confidence levels, respectively. They compare well with the DNM mass evolution, M H DNM M H 2 CO 0.4 $M_{\textrm{H}}^{\textrm{DNM}} \propto {M_{{\textrm{H}_2}}^{\textrm{CO}}}^{\,0.4}$, found at latitudes |b| > 5° in broader and more massive molecular complexes than the clouds studied here (Grenier et al. 2005). The former relation had been derived with a coarser derivation of the dust optical depth (from IRAS and DIRBE) and lower-resolution data in γ rays (from EGRET) and in H I lines (LAB survey, Kalberla et al. 2005). The data points in Fig. 12 strengthen the reliability of a power-law relation between the DNM mass and the CO-bright H2 component. They also extend its application to lower CO masses. According to Larson’s laws, gravitationally bound molecular clouds have masses that approximately scale with the square of the cloud linear size (Larson 1981; Lombardi et al. 2010). For bound clouds, the scaling relation of Eq. (4) implies that the mass present in the DNM envelopes should scale with the linear size of the molecular clouds. This can be observationally tested by using the kinematical information of the H I and CO lines when firm distances to these clouds can be inferred from the stellar reddening information from Gaia. This relation between masses in the DNM and CO phases can also serve to extrapolate the amount of DNM gas beyond the local ISM and in external galaxies. Its confirmation is of prime importance, but it requires further tests in giant molecular clouds of larger CO masses than in Fig. 12, in compact clouds of small size, and in clouds of different metallicity.

The PDR model of Levrier et al. (2012) yields a mass fraction of 32% in a cloudlet with only 9.5 M in CO and illuminated by a standard UV field. The model of Wolfire et al. (2010) yields a narrow range of fractions between 25% and 33% for giant molecular clouds with masses of (0.1–3) × 106 M much largerthan those probed in our sample and with a UV radiation field 3 to 30 times larger than the local value of Draine (1978). The theoretical studies therefore point to rather stable CO-dark H2 fractions that should little depend on cloud mass or radiation field. Our results challenge this view by showing a wide range of mass fractions that relate to cloud properties such as the CO diffuseness (SF dense) and the CO brightness (peak WCO intensity or integral of WCO intensities that gives the M H 2 CO $M_{{\textrm{H}_2}}^{\textrm{CO}}$ mass).

Cloud-dependent variations call for care in the interpretation of averages taken over large ensembles of clouds with different evolutionaryor structural properties. For instance, the difference noted between the local clouds seen towards the inner and outer Galaxy by Paradis et al. (2012) does not relate to the Galactocentric distance, but to different combinations of clouds in different states. This could explain the discrepancy noted by Smith et al. (2014) between an “apparent decrease” of F COdk H 2 $F_{{\textrm{COdk \, H}_2}}$ with Galactocentric distance in Paradis et al. (2012) and an increase with this distance in the much broader sample of clouds studied with [C II] lines across the Milky Way (Pineda et al. 2013). Moreover, the present trend for larger F COdk H 2 $F_{{\textrm{COdk \, H}_2}}$ fractions in clouds lying higher above the Galactic plane, if confirmed, should also raise concerns for radial studies across the Galactic disc in constant latitude bands. The rise observed in F COdk H 2 $F_{{\textrm{COdk \, H}_2}}$ from about 10% at 4 kpc up to 60–80% near 10 kpc is not biased by such height variations because it is based on [C II] line measurements taken along the Galactic plane, at b = 0° (Pineda et al. 2013).

In their simulation of the Milky Way, Smith et al. (2014) use a uniform IRSF, metallicity, and mean gas surface density, so they cannot reproduce a global variation of F COdk H 2 $F_{{\textrm{COdk \, H}_2}}$ with Galactocentric radius. In conditions akin to those in the local ISM, they obtain an average fraction of 42% that is almost insensitive to the WCO intensity threshold used to separate CO-bright and CO-dark gas. They have, however, found notable differences in F COdk H 2 $F_{{\textrm{COdk \, H}_2}}$ between the arm and inter-arm regions. Generally the fraction is anti-correlated with the total molecular N H 2 $N_{\textrm{H}_2}$ column density. The molecular gas is predominantly CO-bright in the spiral arms ( 0.1< F COdk H 2 <0.5 $0.1 < F_{{\textrm{COdk \, H}_2}} <0.5$), but the CO-dark H2 mass fraction exceeds 80% in the inter-arm regions where the clouds are sheared out by the differential rotation of the disc following their passage through the spiral arm. The stretched cloud structures become more permeable to dissociating radiation. Interestingly, we find in our sample a wide range of CO-dark H2 mass fractions responding to properties such as SFdense, M H 2 CO $M_{{\textrm{H}_2}}^{\textrm{CO}}$, and Galactic height, that can reflect the transitory/sheared versus gravitationally well-bound state of a cloud.

Figure 13 shows the distribution of WCO intensities as a function of N H 2 $N_{\textrm{H}_2}$ column densities observed in the anticentre and Chamaeleon regions. The data points above the grey line in this figure correspond to directions where XCO < 1020 cm−2 K−1 km−1s (as in the Chamaeleon, main Taurus, California, and Perseus clouds). The data points gathering along the horizontal axis, at null or very low WCO up to large N H 2 $N_{\textrm{H}_2}$ values, correspond to the DNM gas. In the simulation of Smith et al. (2014) with local-ISM characteristics, CO molecules predominantly form in spiral-arm clouds, which emit mostly at N H 2 > 10 21 $N_{\textrm{H}_2} >10^{21}$ cm−2 with WCO > 10 K km s−1 as shown bythe blue and red contours in Fig. 13. The comparison highlights that the simulation produces underluminous CO clouds at low column densities compared to the observations. This was also noted by Levrier et al. (2012). The CO deficit in their simulation reached a factor of ten in column density in the peripheral regions exposed to the UV radiation, where N H 2 2× 10 20 $N_{\textrm{H}_2} \lesssim 2\times 10^{20}$ cm−2. So, there seems to be a general problem with our understanding of chemistry at low gas density in the photo-dominated regions (PDR), either in thetreatment of the UV attenuation or in the absence of warm chemistry driven by intermittent energy deposition in regions where the interstellar turbulence is dissipated (Godard et al. 2014). This under-prediction of CO emission biases the simulated values of CO-dark H2 fractions upward in column density and in mass. Additionally simulations do not predict bright CO emission at larger N H 2 $N_{\textrm{H}_2}$ column densities of a few 1021 cm−2, as observed in nearby molecular clouds. This explains why simulations obtain twice larger XCO factors thanthe values measured in several nearby clouds. We discussed this discrepancy in Remy et al. (2017). These specifics of Smith et al. (2014) simulation explain why the simulated CO emission gathers primarily in giant molecular clouds in the spiral arms and why the simulations yield a sharp contrast in CO-dark H2 abundance inside and outside of the arms.

Table 2

Fractions of the total mass in the different gas phases for the 15 clouds and their CO-dark H2 mass fractions.

thumbnail Fig. 10

Number distributions of the HI, DNM, CO, and COsat fractions of the total mass found in the sample of clouds for each phase.

thumbnail Fig. 11

Average fraction of the molecular DNM in the total H2 mass ( F COdk H 2 $F_{{\textrm{COdk \, H}_2}}$) as a functionof: the surface fraction of dense gas SFdense in the cloud (top left), the maximum CO integrated intensity W CO max $W_{\textrm{CO}}^{\textrm{max}}$ recorded in the cloud (top right), the average visual extinction A V ¯ $\overline{A_{\textrm{V}}}$ over the cloud (middle left), the logarithm of the H2 mass in the CO phase M H 2 CO $M_{{\textrm{H}_2}}^{\textrm{CO}}$ (middle right), and the height above the Galactic plane (bottom left). The DNM is assumed to be 50% molecular. The colours refer to the same clouds as in Fig. 1.

thumbnail Fig. 12

Gas mass in the DNM (both atomic and molecular) as a function of the H2 mass in the CO-bright phase. The dashed and dotted lines correspond to the best power-law fit respectively including and excluding the compact ChaEastI filament (black point). The colours refer to the same clouds as in Fig. 1.

thumbnail Fig. 13

Number distribution of the integrated 12CO intensity, WCO, as a functionof the total H2 column-density derived from the γ-ray fits for the anticentre and Chamaeleon regions. In each direction N H 2 tot $N_{{\textrm{H}_2}}^{\textrm{tot}}$ is the sum of N H 2 DNM $N_{{\textrm{H}_2}}^{\textrm{DNM}}$, N H 2 CO $N_{{\textrm{H}_2}}^{\textrm{CO}}$, and N H 2 COsat $N_{{\textrm{H}_2}}^{\textrm{COsat}}$ contributions. The DNM is assumed to be 50% molecular. We use the XCO factors derived from γ rays to estimate N H 2 CO $N_{{\textrm{H}_2}}^{\textrm{CO}}$. The grey line is an arbitrary reference corresponding to XCO = 1020 cm−2 K−1 km−1 s. The blue and red contours delimit the least and most populated parts of the distribution found in a simulation of the Milky Way by Smith et al. (2014).

4 Conclusions

We have analysed the gas, dust, and CR content of several nearby clouds in the anticentre and Chamaeleon regions. As described in Planck Collaboration Int. XXVIII (2015) and Remy et al. (2017), in order to trace the total gas column density, we modelled the γ-ray emission and dust optical depth at 353 GHz as linear combinations of several gaseous components associated with different phases (ionized, H I, DNM, CO, and COsat media). We selected a subset of 15 clouds with no overlap in direction to avoid confusion between the separate phases of the various clouds. Thereby, we have been able to follow the relative contributions of the different gas phases to the total mass and column density, and to study phase transitions in the clouds. The main results are summarised as follows:

  • On the gas content: H I is the dominant form of gas by mass in our sample; it contributes 50% to 80% of the total mass. The capability of CO lines to trace the total H2 gas varies widely from cloud to cloud. The DNM envelopes, composed of dense optically thick H I and diffuse CO-dark H2 at the transition between the atomic and molecular phases, contribute about 20% of the cloud masses in our sample. The DNM mass often exceeds the H2 mass seen in CO in diffuse clouds (see below) and it approximately scales with the square root of the CO-bright mass. The CO-saturated H2 gas, which is seen in addition to the WCO-inferred component in dense regions of large optical thickness to 12CO lines, generallycontributes less than 10% to the total cloud mass. It accounts for 40% of the molecular column densities at AV > 3 mag.

  • On phase transitions: the transition from the H I to H2 composition and the onset of visible CO emission both occur in narrow intervals of total gas column density in the selected clouds, at (0.6–1) × 1021 cm−2 and (0.6–2.5) × 1021 cm−2, respectively. We observed a systematic decrease in the H I fraction in the total column density above NH ≈ 1021 cm−2 (AV ≈ 0.4 mag). The CO intensities rise in the 1–2.5 mag range. The 12CO emission significantly saturates at AV > 3 mag.

  • On the CO-dark H2 fractions: the f COdk H 2 $f_{{\textrm{COdk \, H}_2}}$ fraction in column density is observed to decrease from the edge to the core of an individual cloud or a broader complex. The f COdk H 2 $f_{{\textrm{COdk \, H}_2}}$ decrease with increasing AV follows a half-Gaussian profile peaking at 0.4–0.6 mag and with a standard deviation close to 1 mag. The dust specific power reflecting the strength of the ISRF varies only marginally across the DNM envelopes and it does not correlate with spatial changes in f COdk H 2 $f_{{\textrm{COdk \, H}_2}}$. In the solar neighbourhood, differences in the F COdk H 2 $F_{{\textrm{COdk \, H}_2}}$ mass fraction primarily stem from the cloud structure: F COdk H 2 $F_{{\textrm{COdk \, H}_2}}$ significantly decreases from diffuse to dense molecular clouds. The DNM molecular mass often exceeds the H2 mass seen in CO for clouds with peak W CO max $W_{\textrm{CO}}^{\textrm{max}}$ intensities lower than about 15–20 K km s−1, with surface fractions of bright CO intensities SFdense less than about 0.3, or with average visual extinctions A V ¯ $\overline{A_{\textrm{V}}}$ lower than about 1 mag. We note a trend for larger F COdk H 2 $F_{{\textrm{COdk \, H}_2}}$ fractions in clouds located at higher latitudes above the Galactic plane.

In future studies, it will be interesting to compare maps of the total gas derived jointly from γ rays and dust optical depth, or dust extinction, to the gas traced by other species. In diffuse environments, CO precursors such as CH (Gerin et al. 2010), OH (Allen et al. 2015), and HCO (Liszt et al. 2014) can probe the CO dark H2 in the DNM with less confusion than [C II] lines. In dense molecular environments, CO isotopologues such as 13CO and C18 O, and rarer molecules routinely probe the CO-saturated molecular gas at large volume densities. Observations of these tracers are, however, still sparse and their angular extent is often inadequate for comparisons with gas maps inferred from γ-ray and dust studies. Dedicated campaigns towards individual nearby clouds are required to fully explore the relative merits and limitations of the different tracers and to constrain the DNM composition.

Additionally, extending our sample to other clouds in the local ISM, with other masses and various heights above the Galactic plane, then at larger distance in the Local Arm and in between the Local and Perseus arms, is necessary to better characterise the trends in F COdk H 2 $F_{{\textrm{COdk \, H}_2}}$ discussed here. The results would greatly improve our view of the relative distributions of CO-dark and CO-bright H2 gas across the Galaxy.

Acknowledgements

The Fermi-LAT Collaboration acknowledges generous ongoing support from a number of agencies and institutes that have supported both the development and the operation of the LAT as well as scientific data analysis. These include the National Aeronautics and Space Administration and the Department of Energy in the United States, the Commissariat à l’Énergie Atomique and the Centre National de la Recherche Scientifique/Institut National de Physique Nucléaire et de Physique des Particules in France, the Agenzia Spaziale Italiana and the Istituto Nazionale di Fisica Nucleare in Italy, the Ministry of Education, Culture, Sports, Science and Technology (MEXT), High Energy Accelerator Research Organization (KEK) and Japan Aerospace Exploration Agency (JAXA) in Japan, and the K. A. Wallenberg Foundation, the Swedish Research Council and the Swedish National Space Board in Sweden. Additional support for science analysis during the operations phase is gratefully acknowledged from the Istituto Nazionale di Astrofisica in Italy and the Centre National d’Études Spatiales in France. The authors acknowledge the support of the Agence Nationale de la Recherche (ANR) under award number STILISM ANR-12- BS05-0016.

References

  1. Allen, R. J., Hogg, D. E., & Engelke, P. D. 2015, AJ, 149, 123 [NASA ADS] [CrossRef] [Google Scholar]
  2. André, P., Di Francesco, J., Ward-Thompson, D., et al. 2014, Protostars and Planets VI (Tucson: University of Arizona Press), 27 [Google Scholar]
  3. Barriault, L., Joncas, G., Lockman, F. J., & Martin, P. G. 2010, MNRAS, 407, 2645 [NASA ADS] [CrossRef] [Google Scholar]
  4. Bell, T. A., Roueff, E., Viti, S., & Williams, D. A. 2006, MNRAS, 371, 1865 [NASA ADS] [CrossRef] [Google Scholar]
  5. Bertram, E., Glover, S. C. O., Clark, P. C., Ragan, S. E., & Klessen, R. S. 2016, MNRAS, 455, 3763 [NASA ADS] [CrossRef] [Google Scholar]
  6. Blitz, L., Bazell, D., & Desert, F. X. 1990, ApJ, 352, L13 [NASA ADS] [CrossRef] [Google Scholar]
  7. Cardelli, J. A., Clayton, G. C., & Mathis, J. S. 1989, ApJ, 345, 245 [NASA ADS] [CrossRef] [Google Scholar]
  8. Chen, B.-Q., Liu, X.-W., Yuan, H.-B., Huang, Y., & Xiang, M.-S. 2015, MNRAS, 448, 2187 [NASA ADS] [CrossRef] [Google Scholar]
  9. Dame, T. M., & Thaddeus, P., 2004, in Milky Way Surveys: The Structure and Evolution of our Galaxy, eds. D. Clemens, R. Shah, & T. Brainerd, ASP Conf. Ser., 317, 66 [NASA ADS] [Google Scholar]
  10. Dame, T. M., Hartmann, D., & Thaddeus, P. 2001, ApJ, 547, 792 [NASA ADS] [CrossRef] [Google Scholar]
  11. Deguchi, S., & Watson, W. D. 1985, ApJ, 290, 578 [NASA ADS] [CrossRef] [Google Scholar]
  12. Draine, B. T. 1978, ApJS, 36, 595 [NASA ADS] [CrossRef] [Google Scholar]
  13. Field, G. B., Goldsmith, D. W., & Habing, H. J. 1969, ApJ, 155, L149 [NASA ADS] [CrossRef] [Google Scholar]
  14. Fukui, Y., Torii, K., Onishi, T., et al. 2015, ApJ, 798, 6 [NASA ADS] [CrossRef] [Google Scholar]
  15. Gerin, M., de Luca, M., Goicoechea, J. R., et al. 2010, A&A, 521, L16 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  16. Gillmon, K., Shull, J. M., Tumlinson, J., & Danforth, C. 2006, ApJ, 636, 891 [NASA ADS] [CrossRef] [Google Scholar]
  17. Glover, S. C. O., & Mac Low, M.-M. 2011, MNRAS, 412, 337 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
  18. Godard, B., Falgarone, E., & Pineau des Forêts, G. 2014, A&A, 570, A27 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  19. Green, G. M., Schlafly, E. F., Finkbeiner, D. P., et al. 2015, ApJ, 810, 25 [NASA ADS] [CrossRef] [Google Scholar]
  20. Grenier, I. A., Casandjian, J.-M., & Terrier, R. 2005, Science, 307, 1292 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
  21. Grenier, I. A., Black, J. H., & Strong, A. W. 2015, ARA&A, 53, 199 [NASA ADS] [CrossRef] [Google Scholar]
  22. Juvela, M., & Montillaud, J. 2016, A&A, 585, A38 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  23. Kalberla, P. M. W., & Haud, U. 2015, A&A, 578, A78 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  24. Kalberla, P. M. W., Burton, W. B., Hartmann, D., et al. 2005, A&A, 440, 775 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  25. Kalberla, P. M. W., McClure-Griffiths, N. M., Pisano, D. J., et al. 2010, A&A, 521, A17 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  26. Köhler, M., Ysard, N., & Jones, A. P. 2015, A&A, 579, A15 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  27. Könyves, V., André, P., Men’shchikov, A., et al. 2015, A&A, 584, A91 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  28. Krumholz, M. R., McKee, C. F., & Tumlinson, J. 2009, ApJ, 693, 216 [NASA ADS] [CrossRef] [Google Scholar]
  29. Lada, C. J., Lombardi, M., & Alves, J. F. 2010, ApJ, 724, 687 [NASA ADS] [CrossRef] [Google Scholar]
  30. Langer, W. D., Velusamy, T., Pineda, J. L., Willacy, K., & Goldsmith, P. F. 2014, A&A, 561, A122 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  31. Larson, R. B. 1981, MNRAS, 194, 809 [NASA ADS] [CrossRef] [Google Scholar]
  32. Lee, M.-Y., Stanimirović, S., Douglas, K. A., et al. 2012, ApJ, 748, 75 [NASA ADS] [CrossRef] [Google Scholar]
  33. Levrier, F., Le Petit, F., Hennebelle, P., et al. 2012, A&A, 544, A22 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  34. Liszt, H. 2001, A&A, 371, 698 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  35. Liszt, H. 2014, ApJ, 783, 17 [NASA ADS] [CrossRef] [Google Scholar]
  36. Liszt, H. S., Pety, J., & Lucas, R. 2010, A&A, 518, A45 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  37. Liszt, H. S., & Pety, J. 2012, A&A, 541, A58 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  38. Liszt, H. S., Pety, J., Gerin, M., & Lucas, R. 2014, A&A, 564, A64 [Google Scholar]
  39. Lombardi, M., Alves, J., & Lada, C. J. 2010, A&A, 519, L7 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  40. Magnani, L., Chastain, R. J., Kim, H. C., et al. 2003, ApJ, 586, 1111 [NASA ADS] [CrossRef] [Google Scholar]
  41. Martin, P. G., Roy, A., Bontemps, S., et al. 2012, ApJ, 751, 28 [NASA ADS] [CrossRef] [Google Scholar]
  42. McKee, C. F., & Ostriker, J. P. 1977, ApJ, 218, 148 [NASA ADS] [CrossRef] [Google Scholar]
  43. Mizuno, A., Yamaguchi, R., Tachihara, K., et al. 2001, PASJ, 53, 1071 [NASA ADS] [CrossRef] [Google Scholar]
  44. Mizuno, T., Abdollahi, S., Fukui, Y., et al. 2016, ApJ, 833, 278 [NASA ADS] [CrossRef] [Google Scholar]
  45. Narayanan, G., Heyer, M. H., Brunt, C., et al. 2008, ApJS, 177, 341 [NASA ADS] [CrossRef] [Google Scholar]
  46. Paradis, D., Dobashi, K., Shimoikura, T., et al. 2012, A&A, 543, A103 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  47. Peek, J. E. G., Heiles, C., Douglas, K. A., et al. 2011, ApJS, 194, 20 [NASA ADS] [CrossRef] [Google Scholar]
  48. Pineda, J. L., Goldsmith, P. F., Chapman, N., et al. 2010, ApJ, 721, 686 [NASA ADS] [CrossRef] [Google Scholar]
  49. Pineda, J. L., Langer, W. D., Velusamy, T., & Goldsmith, P. F. 2013, A&A, 554, A103 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  50. Planck Collaboration XIX 2011, A&A, 536, A19 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  51. Planck Collaboration XXV 2011, A&A, 536, A25 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  52. Planck Collaboration XI 2014, A&A, 571, A11 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  53. Planck Collaboration XIII 2014, A&A, 571, A13 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  54. Planck Collaboration Int XXVIII 2015, A&A, 582, A31 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  55. Remy, Q., Grenier, I. A., Marshall, D. J., & Casandjian, J. M. 2017, A&A, 601, A78 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  56. Ridge, N. A., Di Francesco, J., Kirk, H., et al. 2006, AJ, 131, 2921 [NASA ADS] [CrossRef] [Google Scholar]
  57. Röhser, T., Kerp, J., Winkel, B., Boulanger, F., & Lagache, G. 2014, A&A, 564, A71 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  58. Rowles, J., & Froebrich, D. 2009, MNRAS, 395, 1640 [NASA ADS] [CrossRef] [Google Scholar]
  59. Sheffer, Y., Rogers, M., Federman, S. R., et al. 2008, ApJ, 687, 1075 [NASA ADS] [CrossRef] [Google Scholar]
  60. Shetty, R., Glover, S. C., Dullemond, C. P., & Klessen, R. S. 2011, MNRAS, 412, 1686 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
  61. Shimajiri, Y., André, P., Braine, J., et al. 2017, A&A, 604, A74 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  62. Smith, R. J., Glover, S. C. O., Clark, P. C., Klessen, R. S., & Springel, V. 2014, MNRAS, 441, 1628 [NASA ADS] [CrossRef] [Google Scholar]
  63. Stanimirović,S., Murray, C. E., Lee, M.-Y., Heiles, C., & Miller, J. 2014, ApJ, 793, 132 [NASA ADS] [CrossRef] [Google Scholar]
  64. Stepnik, B., Abergel, A., Bernard, J.-P., et al. 2003, A&A, 398, 551 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  65. Valdivia, V., Hennebelle, P., Gérin, M., & Lesaffre, P. 2016, A&A, 587, A76 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  66. van Dishoeck, E. F., & Black, J. H. 1988, ApJ, 334, 771 [NASA ADS] [CrossRef] [Google Scholar]
  67. Velusamy, T., Langer, W. D., Pineda, J. L., et al. 2010, A&A, 521, L18 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  68. Visser, R., van Dishoeck, E. F., & Black, J. H. 2009, A&A, 503, 323 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  69. Winkel, B., Kerp, J., Flöer, L., et al. 2016, A&A, 585, A41 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  70. Wolfire, M. G., McKee, C. F., Hollenbach, D., & Tielens, A. G. G. M. 2003, ApJ, 587, 278 [NASA ADS] [CrossRef] [Google Scholar]
  71. Wolfire, M. G., Hollenbach, D., & McKee, C. F. 2010, ApJ, 716, 1191 [NASA ADS] [CrossRef] [Google Scholar]
  72. Xu, D., Li, D., Yue, N., & Goldsmith, P. F. 2016, ApJ, 819, 22 [NASA ADS] [CrossRef] [Google Scholar]
  73. Ysard, N., Abergel, A., Ristorcelli, I., et al. 2013, A&A, 559, A133 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]

All Tables

Table 1

Distance to the clouds, mass in their H I and CO phases, SFR, and XCO factors.

Table 2

Fractions of the total mass in the different gas phases for the 15 clouds and their CO-dark H2 mass fractions.

All Figures

thumbnail Fig. 1

Integrated CO intensity, WCO, overlaid with the boundaries of the selected substructures in the anticentre (left) and Chamaeleon (right) regions.

In the text
thumbnail Fig. 2

Hydrogen column density maps of the HI, DNM (derived from the τ353 and γ-ray fits), CO-bright, and COsat (derived from the τ353 and γ-ray fits) components in the Chamaeleon region. The XCO and τ353NH scale factors used in their derivation come from the γ-ray analysis for optically thin H I.

In the text
thumbnail Fig. 3

Hydrogen column density maps of the HI, DNM (derived from the τ353 and γ-ray fits), CO-bright, and COsat (derived from the τ353 and γ-ray fits) in the anticentre region. The XCO and τ353NH scale factors used in their derivation come from the γ-ray analysis with an H I spin temperature of 400 K.

In the text
thumbnail Fig. 4

Hydrogen column density map of the COsat component from the τ353 analysis overlaid with contours of 13CO line intensities at 2 and 4 K km s−1 in the Taurus and Perseus clouds (Narayanan et al. 2008; Ridge et al. 2006).

In the text
thumbnail Fig. 5

Evolution with N H tot $N_{\textrm{H}}^{\textrm{tot}}$ of the fractions of the total hydrogen column density in the H I (top left), DNM (top right), CO (bottom left), and COsat (bottom right) components. The number of points is the same for all the curves. The points are equally spaced between the minimum and maximum of each profile. The colours refer to the same clouds as in Fig. 1.

In the text
thumbnail Fig. 6

Evolution with AV of the fractions of the total hydrogen column density in the H I (top left), DNM (top right), CO (bottom left), and COsat (bottom right) components. The colours refer to the same clouds as in Fig. 1.

In the text
thumbnail Fig. 7

Evolution of the fraction f COdk H 2 $f_{{\textrm{COdk \, H}_2}}$ of molecular DNM in the total H2 column-density as a function of the total gas column density N H tot $N_{\textrm{H}}^{\textrm{tot}}$ (top) and visual extinction AV (bottom). The DNM gas is assumed to be 50% molecular. The colours refer to the same clouds as in Fig. 1.

In the text
thumbnail Fig. 8

Evolution with visual extinction AV of the fraction f COdk H 2 $f_{{\textrm{COdk \, H}_2}}$ of molecular DNM gas (grey) or additional COsat H2 (green) in the total H2 column-density. The DNM is assumed to be 100% molecular. The shaded areas and error bars respectively give the standard deviation in the sample and the standard error of the mean. The black line gives the Gaussian profile of Xu et al. (2016) and the red line shows the best fit by a half-Gaussian to the f COdk H 2 $f_{{\textrm{COdk \, H}_2}}$ fraction.

In the text
thumbnail Fig. 9

Fraction of the molecular DNM in the total H2 column-density ( f COdk H 2 $f_{{\textrm{COdk \, H}_2}}$) derived fromthe dust (left) and γ-ray (right) analyses for the anticentre (top) and Chamaeleon (bottom) regions. The DNM is assumed to be 50% molecular. The white contours outline the shapes of the CO clouds at the 7 K km s−1 level chosen to separate DNM and COsat components.

In the text
thumbnail Fig. 10

Number distributions of the HI, DNM, CO, and COsat fractions of the total mass found in the sample of clouds for each phase.

In the text
thumbnail Fig. 11

Average fraction of the molecular DNM in the total H2 mass ( F COdk H 2 $F_{{\textrm{COdk \, H}_2}}$) as a functionof: the surface fraction of dense gas SFdense in the cloud (top left), the maximum CO integrated intensity W CO max $W_{\textrm{CO}}^{\textrm{max}}$ recorded in the cloud (top right), the average visual extinction A V ¯ $\overline{A_{\textrm{V}}}$ over the cloud (middle left), the logarithm of the H2 mass in the CO phase M H 2 CO $M_{{\textrm{H}_2}}^{\textrm{CO}}$ (middle right), and the height above the Galactic plane (bottom left). The DNM is assumed to be 50% molecular. The colours refer to the same clouds as in Fig. 1.

In the text
thumbnail Fig. 12

Gas mass in the DNM (both atomic and molecular) as a function of the H2 mass in the CO-bright phase. The dashed and dotted lines correspond to the best power-law fit respectively including and excluding the compact ChaEastI filament (black point). The colours refer to the same clouds as in Fig. 1.

In the text
thumbnail Fig. 13

Number distribution of the integrated 12CO intensity, WCO, as a functionof the total H2 column-density derived from the γ-ray fits for the anticentre and Chamaeleon regions. In each direction N H 2 tot $N_{{\textrm{H}_2}}^{\textrm{tot}}$ is the sum of N H 2 DNM $N_{{\textrm{H}_2}}^{\textrm{DNM}}$, N H 2 CO $N_{{\textrm{H}_2}}^{\textrm{CO}}$, and N H 2 COsat $N_{{\textrm{H}_2}}^{\textrm{COsat}}$ contributions. The DNM is assumed to be 50% molecular. We use the XCO factors derived from γ rays to estimate N H 2 CO $N_{{\textrm{H}_2}}^{\textrm{CO}}$. The grey line is an arbitrary reference corresponding to XCO = 1020 cm−2 K−1 km−1 s. The blue and red contours delimit the least and most populated parts of the distribution found in a simulation of the Milky Way by Smith et al. (2014).

In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.