H2 Column Density Research Articles

Revisiting the Velocity Dispersion–Size Relation in Molecular Cloud Structures

Structures in molecular ISM are observed to follow a power-law relation between the velocity dispersion and spatial size, known as Larson’s first relation, which is often attributed to the turbulent nature of molecular ISM and imprints the dynamics of molecular cloud structures. Using the 13CO (J = 1–0) data from the Milky Way Imaging Scroll Painting survey, we built a sample with 360 structures having relatively accurate distances obtained from either the reddened background stars with Gaia parallaxes or associated maser parallaxes, spanning from 0.4 to ∼15 kpc. Using this sample and about 0.3 million pixels, we analyzed the correlations between velocity dispersion, surface/column density, and spatial scales. Our structure-wise results show power-law indices smaller than 0.5 in both the σ v –R eff and σ v –R eff · Σ relations. In the pixel-wise results, the σvpix is statistically scaling with the beam physical size (R s ≡ ΘD/2) in form of σvpix∝Rs0.43±0.03 . Meanwhile, σvpix in the inner Galaxy is statistically larger than the outer side. We also analyzed correlations between σvpix and the H2 column density N(H2), finding that σvpix stops increasing with N(H2) after ≳1022 cm−2. The structures with and without high-column-density (>1022 cm−2) pixels show different σvpix∝N(H2)ξ relations, where the mean (std) ξ values are 0.38 (0.14) and 0.62 (0.27), respectively.

The Spatial Correlation between CN-line- and Dust-continuum-emitting Regions in High-mass Star-forming Clouds

Measuring the strength of a three-dimensional (3D) magnetic field vector is challenging as it is not easy to recognize whether its line-of-sight (LOS) and plane-of-sky (POS) components are obtained from the same region. CN (N = 1–0) emission has been used to get the LOS component of a magnetic field (B LOS) from its Zeeman splitting lines, while dust continuum emission has been used to get the POS component of a magnetic field (B POS). We use the CN (N = 1–0) data observed with the Taeduk Radio Astronomy Observatory 14 m telescope and the dust continuum data from the Herschel archive toward six high-mass star-forming regions in order to test whether CN line and dust continuum emission can trace a similar region and thus can be used for inferring 3D magnetic field strength. Our comparison between CN and H2 column densities for all targets indicates that CN line emission tends to be strong toward bright continuum regions. The positions of peak CN column densities are particularly well correlated with those of peak H2 column densities, at least over the H2 column density of 8.0 × 1022 cm−2 within one or two telescope beam sizes in all targets, implying that CN-line- and dust-continuum-emitting regions are likely spatially coincident. This enabled us to make the reliable measurement of the 3D magnetic field strengths of five targets by taking a vector sum of their B LOS and B POS, helping to decide the magnetical criticality of the targets as supercritical or transcritical.

Investigating the Star-forming Sites in the Outer Galactic Arm

We aim to investigate the global star formation scenario in star-forming sites AFGL 5157, [FSR2007] 0807 (hereafter FSR0807), [HKS2019] E70 (hereafter E70), [KPS2012] MWSC 0620 (hereafter KPS0620), and IRAS 05331+3115 in the outer Galactic arm. The distribution of young stellar objects in these sites coincides with a higher extinction and H2 column density, which agrees with the notion that star formation occurs inside the dense molecular cloud cores. We have found two molecular structures at different velocities in this direction; one contains AFGL 5157 and FSR0807, and the other contains E70, [KPS2012] MWSC 0620, and IRAS 05331+3115. All these clusters in our target region are in different evolutionary stages and might form stars through different mechanisms. The E70 cluster seems to be the oldest in our sample; AFGL 5157 and FSR0807 formed later, and KPS0620 and IRAS 05331+3115 are the youngest sites. AFGL 5157 and FSR0807 are physically connected and have cold filamentary structures and dense hub regions. Additionally, the near-infrared photometric analysis shows signatures of massive star formation in these sites. KPS0620 also seems to have cold filamentary structures with the central hub but lacks signatures of massive stars. Our analysis suggests molecular gas flow and the hub filamentary star formation scenario in these regions. IRAS 05331+3115 is a single clump of molecular gas favoring low-mass star formation. Our study suggests that the selected area is a menagerie of star-forming sites where the formation of the stars happens through different processes.

Temperature stratification in a molecular shock: Analysis of the emission of H2 pure rotational lines in IC443G

Context. Supernovae remnants (SNRs) represent a major source of feedback from stars on the interstellar medium of galaxies. During the latest stage of supernova explosions (which lasts 10–100 kyr), shock waves produced by the initial blast modify the chemistry of gas and dust, inject kinetic energy in the surroundings, and may alter star formation characteristics. Simultaneously, γ-ray emission is generated by the interaction between the ambient medium and cosmic rays, in particular those locally accelerated in the early stages of the explosion. Aims. We aim to estimate the total molecular mass, local density, and total column density of H2 and the temperature structure in a shocked clump interacting with the supernova remnant IC443 located in a region where cosmic rays interact with the interstellar medium. Measuring the mass of the dense and neutral component of the medium is a prerequisite to understanding the chemistry, energetics, and GeV to TeV γ-ray emission. Methods. Assuming that the emission of H2 pure rotational lines is produced by a collection of gas layers with variable temperature, we compared Spitzer/IRS emission maps for the ν = 0–0 S(0) to S(7) lines with a thermal admixture model. Our description is based on a power-law distribution of thermalized components with temperatures varying between Tmin = 25 K and Tmax = 1500 K. Results. Our thermal admixture model allows the level populations of H2 to be described by a power-law distribution dN = ΛT−ΓdT, with Γ ~ 2.2−4.7. We measured a total mass MH2 = 220−80+110 M⊙ across the Spitzer/IRS field of observations. Conclusions. Our analysis shows that an estimate of the cold molecular gas temperature is paramount to accurately constraining the H2 mass, although the mass remains affected by significant uncertainties due to the assumptions on the gas temperature distribution.

The CO-dark molecular gas in the cold H I arc

The CO-dark molecular gas (DMG), which refers to the molecular gas not traced by CO emission, is crucial for the evolution of the interstellar medium (ISM). While the gas properties of DMG have been widely explored in the Solar neighborhood, whether or not they are similar in the outer disk regions of the Milky Way is still not well understood. In this Letter, we confirm the existence of DMG toward a cold H I arc structure at 13 kpc away from the Galactic center with both OH emission and H I narrow self-absorption (HINSA). This is the first detection of HINSA in the outer disk region, in which the HINSA fraction (NHINSA/NH2 = 0.022 ± 0.011) is an order of magnitude higher than the average value observed in nearby evolved dark clouds, but is consistent with that of the early evolutionary stage of dark clouds. The inferred H2 column density from both extinction and OH emission (NH2 ≈ 1020 cm−2) is an order of magnitude higher than previously estimated. Although the ISM environmental parameters are expected to be different between the outer Galactic disk regions and the Solar neighborhood, we find that the visual extinction (AV = 0.19 ± 0.03 mag), H2-gas density (nH2 = 91 ± 46 cm−3), and molecular fraction (58% ± 28%) of the DMG are rather similar to those of nearby diffuse molecular clouds. The existence of DMG associated with the expanding H I supershell supports a scenario where the expansion of supershells may trigger the formation of molecular clouds within a crossing timescale of the shock wave (∼106 yr).

Variation in XCO Factor in N55 Region

The XCO factor is defined as XCO=N(H2)/W12CO. It is useful for estimating cloud mass. However, there is only limited research on how the XCO factor varies within a single cloud. Employing 12CO(J=1-0) and 13CO(J=1-0) spectral data, we computed an XCO factor of 3.6 ×1020cm−2 (K km s−1)−1 for luminous gas of the N55 region. Our analysis revealed a V-shaped correlation between the XCO factor and H2 column densities, while the relationship with excitation temperature exhibited obscurity. This suggests that the CO-to-H2 conversion is not consistent on small scale (∼1 pc). Additionally, we found that star formation activity has little influence on the variability in the XCO factor.

FAUST

Context. The origin of the chemical diversity observed around low-mass protostars probably resides in the earliest history of these systems. Aims. We aim to investigate the impact of protostellar feedback on the chemistry and grain growth in the circumstellar medium of multiple stellar systems. Methods. In the context of the ALMA Large Program FAUST, we present high-resolution (50 au) observations of CH3OH, H2CO, and SiO and continuum emission at 1.3 mm and 3 mm towards the Corona Australis star cluster. Results. Methanol emission reveals an arc-like structure at ∼1800 au from the protostellar system IRS7B along the direction perpendicular to the major axis of the disc. The arc is located at the edge of two elongated continuum structures that define a cone emerging from IRS7B. The region inside the cone is probed by H2CO, while the eastern wall of the arc shows bright emission in SiO, a typical shock tracer. Taking into account the association with a previously detected radio jet imaged with JVLA at 6 cm, the molecular arc reveals for the first time a bow shock driven by IRS7B and a two-sided dust cavity opened by the mass-loss process. For each cavity wall, we derive an average H2 column density of ∼7 × 1021 cm−2, a mass of ∼9 × 10−3 M⊙, and a lower limit on the dust spectral index of 1.4. Conclusions. These observations provide the first evidence of a shock and a conical dust cavity opened by the jet driven by IRS7B, with important implications for the chemical enrichment and grain growth in the envelope of Solar System analogues.

Ultraviolet H2luminescence in molecular clouds induced by cosmic rays

Context. Galactic cosmic rays (CRs) play a crucial role in ionisation, dissociation, and excitation processes within dense cloud regions where UV radiation is absorbed by dust grains and gas species. CRs regulate the abundance of ions and radicals, leading to the formation of more and more complex molecular species, and determine the charge distribution on dust grains. A quantitative analysis of these effects is essential for understanding the dynamical and chemical evolution of star-forming regions.Aims. The CR-induced photon flux has a significant impact on the evolution of the dense molecular medium in its gas and dust components. This study evaluates the flux of UV photons generated by CRs to calculate the photon-induced dissociation and ionisation rates of a vast number of atomic and molecular species, as well as the integrated UV photon flux.Methods. To achieve these goals, we took advantage of recent developments in the determination of the spectra of secondary electrons, in the calculation of state-resolved excitation cross sections of H2by electron impact, and of photodissociation and photoionisation cross sections.Results. We calculated the H2level population of each rovibrational level of theX, B, C, B′,D, B″,D′, andastates. We then computed the UV photon spectrum of H2in its line and continuum components between 72 and 700 nm, with unprecedented accuracy, as a function of the CR spectrum incident on a molecular cloud, the H2column density, the isomeric H2composition, and the dust properties. The resulting photodissociation and photoionisation rates are, on average, lower than previous determinations by a factor of about 2, with deviations of up to a factor of 5 for the photodissociation of species such as AlH, C2H2, C2H3, C3H3, LiH, N2, NaCl, NaH, O2+, S2, SiH, l-C4, and l-C5H. A special focus is given to the photoionisation rates of H2, HF, and N2, as well as to the photodissociation of H2, which we find to be orders of magnitude higher than previous estimates. We give parameterisations for both the photorates and the integrated UV photon flux as a function of the CR ionisation rate, which implicitly depends on the H2column density, as well as the dust properties.

Absorption of Millimeter-band CO and CN in the Early Universe: Molecular Clouds in the Radio Galaxy B2 0902+34 at Redshift 3.4

Using the Karl G. Jansky Very Large Array, we have detected absorption lines due to carbon monoxide, CO(J = 0 → 1), and the cyano radical, CN(N = 0 → 1), associated with radio galaxy B2 0902+34 at redshift z = 3.4. The detection of millimeter-band absorption observed 1.5 Gyr after the Big Bang facilitates studying molecular clouds down to gas masses inaccessible to emission-line observations. The CO absorption in B2 0902+34 has a peak optical depth of τ ≥ 8.6% and consists of two components, one of which has the same redshift as previously detected 21 cm absorption of neutral hydrogen (H i) gas. Each CO component traces an integrated H2 column density of NH2 ≳ 3 × 1020 cm−2. CN absorption is detected for both CO components, as well as for a blueshifted component not detected in CO, with CO/CN line ratios ranging from ≲0.4 to 2.4. We discuss the scenario that the absorption components originate from collections of small and dense molecular clouds that are embedded in a region with more diffuse gas and high turbulence, possibly within the influence of the central active galactic nucleus or a starburst region. The degree of reddening in B2 0902+34, with rest-frame color B − K ∼ 4.2, is lower than the very red colors (B − K > 6) found among other known redshifted CO absorption systems at z < 1. Nevertheless, when including the many nondetections from the literature, a potential correlation between the absorption-line strength and B − K color is evident, giving weight to the argument that the red colors of CO absorbers are due to a high dust content.

ALMA High-resolution Spectral Survey of Thioformaldehyde (H2CS) toward Massive Protoclusters

Investigating the temperature and density structures of gas in massive protoclusters is crucial for understanding the chemical properties therein. In this study, we present observations of the continuum and thioformaldehyde (H2CS) lines at 345 GHz of 11 massive protoclusters using the Atacama Large Millimeter/submillimeter Array telescope. High spatial resolution and sensitivity observations have detected 145 continuum cores from the 11 sources. H2CS line transitions are observed in 72 out of 145 cores, including line-rich cores, warm cores, and cold cores. The H2 column densities of the 72 cores are estimated from the continuum emission, which are larger than the density threshold value for star formation, suggesting that H2CS can be widely distributed in star-forming cores with different physical environments. The rotation temperature and column density of H2CS are derived using the XCLASS software. The results show that the H2CS abundances increase as temperature rises and higher gas temperatures are usually associated with higher H2CS column densities. The abundances of H2CS are positively correlated with its column density, suggesting that the H2CS abundances are enhanced from cold cores, warm cores, and line-rich cores in star-forming regions.

Formaldehyde observations of the Perseus Molecular Cloud

ABSTRACT Large-scale observations of the Perseus Molecular Cloud (MC) with Nanshan 26-m telescope are presented using the 6 cm ortho-H2CO (110–111) transition. As a probe of dense gas at low temperatures, the H2CO absorption extends over the main parts of the Perseus MC. A comparison of the H2CO absorption line with the 12CO and 13CO (J = 1–0) emissions shows that the H2CO and CO are similarly distributed over the Perseus but that H2CO correlates better with 13CO. Comparison with the Herschel-derived H2 column density shows that both the 13CO and H2CO column densities vary linearly with the H2 column density. The main parameters of H2CO absorption line data show a log-normal distribution, which suggests that the strong non-thermal line-broadening results from large-scale supersonic turbulence related to the star formation. Formaldehyde absorption serves well as a tracer of star formation activity and also the H2CO-to-13CO and H2CO-to-H2 abundances systematically trace the star formation activity in the six subregions of Perseus MC. The H2CO abundances anticorrelate with the number of prestellar and protostellar cores and the IRIS 12 μm flux in the six subregions and reveal the star formation history in the Perseus MC.

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

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

Characterizing the line emission from molecular clouds

Aims. We aim to characterize and compare the molecular-line emission of three clouds whose star-formation rates span one order of magnitude: California, Perseus, and Orion A. Methods. We used stratified random sampling to select positions representing the different column density regimes of each cloud and observed them with the IRAM 30 m telescope. We covered the 3 mm wavelength band and focused our analysis on CO, HCN, CS, HCO+, HNC, and N2H+. Results. We find that the line intensities depend most strongly on the H2 column density, with which they are tightly correlated. A secondary effect, especially visible in Orion A, is a dependence of the line intensities on the gas temperature. We explored a method that corrects for temperature variations and show that, when it is applied, the emission from the three clouds behaves very similarly. CO intensities vary weakly with column density, while the intensity of traditional dense-gas tracers such as HCN, CS, and HCO+ varies almost linearly with column density. N2H+ differs from all other species in that it traces only cold dense gas. The intensity of the rare HCN and CS isotopologs reveals additional temperature-dependent abundance variations. Overall, the clouds have similar chemical compositions that, as the depth increases, are sequentially dominated by photodissociation, gas-phase reactions, molecular freeze-out, and stellar feedback in the densest parts of Orion A. Our observations also allowed us to calculate line luminosities for each cloud, and a comparison with literature values shows good agreement. We used our HCN(1–0) data to explore the behavior of the HCN conversion factor, finding that it is dominated by the emission from the outermost cloud layers. It also depends strongly on the gas kinetic temperature. Finally, we show that the HCN/CO ratio provides a gas volume density estimate, and that its correlation with the column density resembles that found in extragalactic observations.

Protonated hydrogen cyanide as a tracer of pristine molecular gas

Context. Protonated hydrogen cyanide, HCNH+, plays a fundamental role in astrochemistry because it is an intermediary in gas-phase ion-neutral reactions within cold molecular clouds. However, the impact of the environment on the chemistry of HCNH+ remains poorly understood. Aims. We aim to study HCNH+, HCN, and HNC, as well as two other chemically related ions, HCO+ and N2H+, in different star formation regions in order to investigate how the environment influences the chemistry of HCNH+. Methods. With the IRAM 30 m and APEX 12 m telescopes, we carried out HCNH+, H13CN, HN13C, H13CO+, and N2H+ imaging observations toward two dark clouds, the Serpens filament and Serpens South, both of which harbor sites of star formation that include protostellar objects and regions that are quiescent. Results. We report the first robust distribution of HCNH+ in the Serpens filament and in Serpens South. Our data suggest that HCNH+ is abundant in cold and quiescent regions but is deficient in active star-forming regions. The observed HCNH+ fractional abundances relative to H2 range from 3.1 × 10−11 in protostellar cores to 5.9 × 10−10 in prestellar cores, and the HCNH+ abundance generally decreases with increasing H2 column density, which suggests that HCNH+ coevolves with cloud cores. Our observations and modeling results suggest that the abundance of HCNH+ in cold molecular clouds is strongly dependent on the H2 number density. The decrease in the abundance of HCNH+ is caused by the fact that its main precursors (e.g., HCN and HNC) undergo freeze-out as the number density of H2 increases. However, current chemical models cannot explain other observed trends, such as the fact that the abundance of HCNH+ shows an anticorrelation with that of HCN and HNC but a positive correlation with that of N2H+ in the southern part of Serpens South’s northern clump. This indicates that additional chemical pathways have to be invoked for the formation of HCNH+ via molecules such as N2 in regions in which HCN and HNC freeze out. Conclusions. Both the fact that HCNH+ is most abundant in molecular cores prior to gravitational collapse and the fact that low-J HCNH+ transitions have very low H2 critical densities make this molecular ion an excellent probe of pristine molecular gas.

A global view on star formation: The GLOSTAR Galactic plane survey

Context. Cygnus X is one of the closest and most active high-mass star-forming regions in our Galaxy, making it one of the best laboratories for studying massive star formation. Aims. We aim to investigate the properties of molecular gas structures on different linear scales with the 4.8 GHz formaldehyde (H2CO) absorption line in Cygnus X. Methods. As part of the GLOSTAR Galactic plane survey, we performed large-scale (7º×3º) simultaneous H2CO (11,0–11,1) spectral line and radio continuum imaging observations toward Cygnus X at λ ~6 cm with the Karl G. Jansky Very Large Array and the Effelsberg 100 m radio telescope. We used auxiliary HI, 13CO (1–0), dust continuum, and dust polarization data for our analysis. Results. Our Effelsberg observations reveal widespread H2CO (11,0–11,1) absorption with a spatial extent of ≳50 pc in Cygnus X for the first time. On large scales of 4.4 pc, the relative orientation between the local velocity gradient and the magnetic field tends to be more parallel at H2 column densities of ≳1.8×1022 cm−2. On the smaller scale of 0.17 pc, our VLA+Effelsberg combined data reveal H2CO (11,0–11,1) absorption only towards three bright HII regions. Our observations demonstrate that H2CO (11,0–11,1) is optically thin in general. The kinematic analysis supports the assertion that molecular clouds generally exhibit supersonic motions on scales of 0.17−4.4 pc. We show a non-negligible contribution of the cosmic microwave background radiation to the extended absorption features in Cygnus X. Our observations suggest that H2CO (11,0–11,1) can trace molecular gas with H2 column densities of ≳5 × 1021 cm−2 (i.e., AV ≳ 5). The ortho-H2CO fractional abundance with respect to H2 has a mean value of 7.0 × 10−10. A comparison of the velocity dispersions on different linear scales suggests that the velocity dispersions of the dominant −3 km s−1 velocity component in the prominent DR21 region are nearly identical on scales of 0.17−4.4 pc, which deviates from the expected behavior of classic turbulence.

Predicting reliable H2 column density maps from molecular line data using machine learning

ABSTRACT The total mass estimate of molecular clouds suffers from the uncertainty in the H2-CO conversion factor, the so-called XCO factor, which is used to convert the 12CO (1–0) integrated intensity to the H2 column density. We demonstrate the machine learning’s ability to predict the H2 column density from the 12CO, 13CO, and C18O (1–0) data set of four star-forming molecular clouds: Orion A, Orion B, Aquila, and M17. When the training is performed on a subset of each cloud, the overall distribution of the predicted column density is consistent with that of the Herschel column density. The total column density predicted and observed is consistent within 10 per cent, suggesting that the machine learning prediction provides a reasonable total mass estimate of each cloud. However, the distribution of the column density for values &amp;gt;∼2 × 1022 cm−2, which corresponds to the dense gas, could not be predicted well. This indicates that molecular line observations tracing the dense gas are required for the training. We also found a significant difference between the predicted and observed column density when we created the model after training the data on different clouds. This highlights the presence of different XCO factors between the clouds, and further training in various clouds is required to correct for these variations. We also demonstrated that this method could predict the column density towards the area not observed by Herschel if the molecular line and column density maps are available for the small portion, and the molecular line data are available for the larger areas.

An ALMA Glimpse of Dense Molecular Filaments Associated with High-mass Protostellar Systems in the Large Magellanic Cloud

Recent millimeter/submillimeter facilities have revealed the physical properties of filamentary molecular clouds in relation to high-mass star formation. A uniform survey of the nearest, face-on star-forming galaxy, the Large Magellanic Cloud (LMC), complements the Galactic knowledge. We present ALMA survey data with a spatial resolution of ∼0.1 pc in the 0.87 mm continuum and HCO+ (4–3) emission toward 30 protostellar objects with luminosities of 104–105.5 L ⊙ in the LMC. The spatial distributions of the HCO+ (4–3) line and thermal dust emission are well correlated, indicating that the line effectively traces dense, filamentary gas with an H2 volume density of ≳105 cm−3 and a line mass of ∼103–104 M ⊙ pc−1. Furthermore, we obtain an increase in the velocity line widths of filamentary clouds, which follows a power-law dependence on their H2 column densities with an exponent of ∼0.5. This trend is consistent with observations toward filamentary clouds in nearby star-forming regions within ≲1 kpc from us and suggests enhanced internal turbulence within the filaments due to surrounding gas accretion. Among the 30 sources, we find that 14 are associated with hub-filamentary structures, and these complex structures predominantly appear in protostellar luminosities exceeding ∼5 × 104 L ⊙. The hub-filament systems tend to appear in the latest stages of their natal cloud evolution, often linked to prominent H ii regions and numerous stellar clusters. Our preliminary statistics suggest that the massive filaments accompanied by hub-type complex features may be a necessary intermediate product in forming extremely luminous high-mass stellar systems capable of ultimately dispersing the parent cloud.

Characterizing three-dimensional magnetic field, turbulence, and self-gravity in the star-forming region L1688

ABSTRACT Interaction of three-dimensional magnetic fields, turbulence, and self-gravity in the molecular cloud is crucial in understanding star formation but has not been addressed so far. In this work, we target the low-mass star-forming region L1688 and use the spectral emissions of 12CO, 13CO, C18O, and H i, as well as polarized dust emissions. To obtain the 3D direction of the magnetic field, we employ the novel polarization fraction analysis. In combining with the plane-of-the-sky (POS) magnetic field strength derived from the Davis–Chandrasekhar–Fermi (DCF) method and the new differential measure analysis (DMA) technique, we present the first measurement of L1688’s three-dimensional magnetic field, including its orientation and strength. We find that L1688’s magnetic field has two statistically different inclination angles. The low-intensity tail has an inclination angle ≈55° on average, while that of the central dense clump is ≈30°. We find the global mean value of total magnetic field strength is Btot ≈ $135 \,\mathrm{\mu }{\rm G}$ from DCF and Btot ≈ $75 \,\mathrm{\mu }{\rm G}$ from DMA. We use the velocity gradient technique (VGT) to separate the magnetic fields’ POS orientation associated with L1688 and its foreground/background. The magnetic fields’ orientations are statistically coherent. The probability density function of H2 column density and VGT reveal that L1688 is potentially undergoing gravitational contraction at large scale ≈1.0 pc and gravitational collapse at small scale ≈0.2 pc. The gravitational contraction mainly along the magnetic field resulting in an approximate power-law relation $B_{\rm tot}\propto n_{\rm H}^{1/2}$ when volume density nH is less than approximately 6.0 × 103 cm−3.

Gas, dust, and the CO-to-molecular gas conversion factor in low-metallicity starbursts

The factor relating CO emission to molecular hydrogen column density, XCO, is still subject to uncertainty, in particular at low metallicity. In this paper, to quantify XCO at two different spatial resolutions, we exploited a dust-based method together with ALMA 12-m and ACA data and H I maps of three nearby metal-poor starbursts, NGC 625, NGC 1705, and NGC 5253. Dust opacity at 250 pc resolution was derived based on dust temperatures estimated by fitting two-temperature modified blackbodies to Herschel PACS data. By using the HI maps, we were then able to estimate dust-to-gas ratios in the regions dominated by atomic gas, and, throughout the galaxy, to infer total gas column densities and H2 column densities as the difference with HI. Finally, from the ACA CO(1–0) maps, we derived XCO. We used a similar technique with 40 pc ALMA 12-m data for the three galaxies, but instead derived dust attenuation at 40 pc resolution from reddening maps based on VLT/MUSE data. At 250 pc resolution, we find XCO ∼ 1022 − 1023 cm−2/K km s−1, 5–1000 times the Milky Way value, with much larger values than would be expected from a simple metallicity dependence. Instead, at 40 pc resolution, XCO again shows large variation, but is roughly consistent with a power-law metallicity dependence, given the Z ∼ 1/3 Z⊙ metal abundances of our targets. The large scatter in both estimations could imply additional parameter dependence, which we have investigated by comparing XCO with the observed velocity-integrated brightness temperatures, ICO, as predicted by recent simulations. Indeed, larger XCO is significantly correlated with smaller ICO, but with slightly different slopes and normalizations than predicted by theory. Such behavior can be attributed to the increasing fraction of CO-faint (or dark) H2 gas with lower spatial resolution (larger beams). This confirms the idea the XCO is multivariate, depending not only on metallicity but also on the CO brightness temperature and beam size. Future work is needed to consolidate these empirical results by sampling galaxies with different metal abundances observed at varying spatial resolutions.

ATLASGAL: 3 mm class I methanol masers in high-mass star formation regions

Context. Class I methanol masers are known to be associated with shocked outflow regions around massive protostars, indicating a possible link between the maser properties and those of their host clumps. Aims. The main goals of this study are (1) to search for new class I methanol masers, (2) to statistically study the relationship between class I masers and shock tracers, (3) to compare the properties between class I masers and their host clumps, also as a function of their evolutionary stage, and (4) to constrain the physical conditions that excite multiple class I masers simultaneously. Methods. We analysed the 3 mm wavelength spectral line survey of 408 clumps identified by the APEX Telescope Large Area Survey of the Galaxy (ATLASGAL), which were observed with the IRAM 30-meter telescope, focusing on the class I methanol masers with frequencies near 84, 95, and 104.3 GHz. Results. We detect narrow maser-like features towards 54, 100, and 4 sources in the maser lines near 84, 95, and 104.3 GHz, respectively. Among them, 50 masers at 84 GHz, 29 masers at 95 GHz, and 4 rare masers at 104.3 GHz are new discoveries. The new detections increase the number of known 104.3 GHz masers from five to nine. The 95 GHz class I methanol maser is generally stronger than the 84 GHz maser counterpart. We find nine sources showing class I methanol masers, but no SiO emission, indicating that class I methanol masers might be the only signpost of protostellar outflow activity in extremely embedded objects at the earliest evolutionary stage. Class I methanol masers that are associated with sources that show SiO line wings are more numerous and stronger than those without such wings. The total integrated intensity of class I methanol masers is well correlated with the integrated intensity and velocity coverage of the SiO (2−1) emission. The properties of class I methanol masers are positively correlated with the bolometric luminosity, clump mass, and peak H2 column density of their associated clumps, but are uncorrelated with the luminosity-to-mass ratio, dust temperature, and mean H2 volume density. Conclusions. We suggest that the properties of class I masers are related to shocks traced by SiO. Based on our observations, we conclude that class I methanol masers at 84 and 95 GHz can trace a similar evolutionary stage to the H2O maser, and appear prior to 6.7 and 12.2 GHz methanol and OH masers. Despite their small number, the 104.3 GHz class I masers appear to trace a shorter and more evolved stage compared to the other class I masers.