Molecular Column Density Research Articles

CHEMOUT: CHEMical complexity in star-forming regions of the OUTer Galaxy. IV. ALMA observations of organic species at a galactocentric radius of ~ 23 kpc

Single-dish observations suggest that the abundances of organic species in star-forming regions of the outer Galaxy, which are characterised by sub-solar metallicities, are comparable to those found in the local Galaxy. To understand this counter-intuitive result and avoid a misleading interpretation due to beam dilution effects at these large distances, spatially resolved molecular emission maps are needed to correctly link the measured abundances and local physical properties. We observed several organic molecules with the Atacama Large Millimeter Array towards WB89-671, the source with the largest galactocentric distance (23.4 kpc) of the project CHEMical complexity in star-forming regions of the OUTer Galaxy (CHEMOUT) at a resolution of $ 15000$ au. We compared the observed molecular abundances with chemical model predictions. We detected emission of C$_4$H HCO HCS$^+$, CS, HN13C, and SO. The emission morphology is complex, extended, and different in each tracer. In particular, the most intense emission in and arises from two millimeter-continuum infrared-bright cores. The most intense and SO emission predominantly arises from the part of the filament that lacks continuum sources. The narrow line widths across the filament indicate quiescent gas in spite of the two embedded protostars. The derived molecular column densities are comparable with those in local star-forming regions, and they suggest an anti-correlation between hydrocarbons, ions, HCO, and on the one hand, and and SO on the other. The static chemical models that match the observed column densities best favour low-energy conditions that are expected at large galactocentric radii, but they also favour carbon elemental abundances that exceed those derived by extrapolating the C/H galactocentric gradient at 23 kpc by three times. This would indicate a flatter C/H trend at large galactocentric radii, which is in line with a flat abundance of organics. However, to properly reproduce the chemical composition of each region, models should include dynamical evolution.

Molecular Oxygen Abundance in Galactic Massive Star Formation Regions Based on SWAS Observations

Molecular oxygen abundance is a key parameter in understanding the chemical network of the interstellar medium. We estimate the molecular oxygen column density and abundance for a sample of Galactic massive star formation regions based on observations from the Submillimiter Wave Astronomy Satellite (SWAS) survey. We obtained an averaged O2 spectrum based on this sample using the (SWAS) survey data (O2, 487.249 GHz, N = 3–1, J = 3–2). No emission or absorption feature is seen around the supposed central velocity with a total integration time of t total = 8.67 × 103 hr and an rms noise per channel of 1.45 mK. Assuming a kinetic temperature T kin = 30 K, we derive the 3σ upper limit of the O2 column density to be 3.3 × 1015 cm−2, close to the lowest values reported in Galactic massive star formation regions in previous studies. The corresponding O2 abundance upper limit is 6.7 × 10−8, lower than all previous results based on SWAS observations and is close to the lowest reported value in massive star formation regions. On a galactic scale, our statistical results confirm a generally low O2 abundance for Galactic massive star formation regions. This abundance is also lower than results reported in extragalactic sources.

Modifications of astrophysical ices induced by cosmic rays

Aims. Astrophysical ices on dust grain mantles in the interstellar medium (ISM) and dense circumstellar envelopes (CSEs) are continuously exposed to galactic cosmic rays (GCRs). In a laboratory setting, we studied the physical and chemical modifications of ice layers induced by energetic heavy ions as GCR analogues. The ice layers used have a molecular composition similar to that of icy grain mantles. Methods. Mixtures of H2O:CO:CH3OH molecules (percentages 73:24:3, 68:30:3, and 58:38:3) were condensed on a substrate at 15 K and irradiated with 40 MeV 58Ni11+ ion beams. Irradiation-induced modifications were followed using the mid-infrared absorption spectroscopy technique. Results. We observed the evolution of infrared bands of CO2, HCO, HCOOH, CH4, H2CO,H2O2, and more complex synthesised molecules. From the molecular column densities, cross-sections and sputtering yields were determined and compared to published results of water and carbon monoxide. Analysis of the chemical modifications reveals that the precursors are easily destroyed when they are in a molecular mixture, while others are desorbed. Conclusions. The main radiolitic modifications induced by GCR irradiations are molecular decomposition and sputtering. Extrapolation to astrophysical radiation conditions shows a strong dependence on the intensity of the GCR distributions at low energies, which allows the analysis of the ice evolution at timescales comparable to those of the ISM and CSE.

ALMA Reveals Spatially Resolved Properties of Molecular Gas in the Host Galaxy of FRB 20191001A at z = 0.2340

We report the detection of the CO(2–1) emission line with a spatial resolution of 0.″9 (3.5 kpc) from the host galaxy of the fast radio burst (FRB), FRB 20191001A at z = 0.2340, using the Atacama Large Millimeter/submillimeter Array. This is the first detection of spatially resolved CO emission from the host galaxy of an FRB at a cosmological distance. The inferred molecular gas mass of the host galaxy is (2.3 ± 0.4) × 1010 M ⊙, indicating that it is gas-rich, as evidenced by the measured molecular gas fraction μ gas = 0.50 ± 0.22. This molecular gas mass and the star formation rate of the host, SFR = (8.06 ± 2.42) M ⊙ yr−1, differ from those observed in the other FRB host galaxies with the average M gas = 9.6 × 108 M ⊙ and SFR = 0.90M ⊙ yr−1. This lends further credibility to the hypothesis that FRBs may originate from single or multiple progenitors across a diverse range of galaxy environments. Based on the observed velocity field modeling, we find that the molecular gas disk is dominated by an ordered circular rotation, despite the fact that the host galaxy has a gas-rich companion galaxy with a projected separation of ∼25 kpc. The formation of the FRB’s progenitor might not have been triggered by this interaction. We derive the 3σ upper limit of the molecular gas column density at the FRB detection site to be <2.1 × 1021 cm−2 with a 3σ upper limit.

Impact of nonconvergence and various approximations of the partition function on the molecular column densities in the interstellar medium (Corrigendum)

Impact of nonconvergence and various approximations of the partition function on the molecular column densities in the interstellar medium <i>(Corrigendum)</i>

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.

Discovery of Thionylimide, HNSO, in Space: The first N-, S-, and O-bearing Interstellar Molecule

We present the first detection in space of thionylimide (HNSO) toward the Galactic center molecular cloud G + 0.693-0.027, thanks to the superb sensitivity of an ultradeep molecular line survey carried out with the Yebes 40 m and IRAM 30 m telescopes. This molecule is the first species detected in the interstellar medium containing, simultaneously, N, S, and O. We have identified numerous K a = 0, 1, and 2 transitions belonging to HNSO covering from J up = 2 to J up =10, including several completely unblended features. We derive a molecular column density of N = (8 ± 1)×1013 cm−2, yielding a fractional abundance relative to H2 of ∼6 × 10−10, which is about ∼37 and ∼4.8 times less abundant than SO and SO2, respectively. Although there are still many unknowns in the interstellar chemistry of NSO-bearing molecules, we propose that HNSO is likely formed through the reaction of the NSO radical and atomic H on the surface of icy grains, with alternative routes also deserving exploration. Finally, HNSO appears as a promising link between N, S, and O interstellar chemistry, and its discovery paves the route to the detection of a new family of molecules in space.

Cold molecules in H I 21 cm absorbers across redshifts ∼0.1–4

Absorption lines at high-redshift in front of quasars are quite rare in the millimeter (mm) domain. Only five associated and five intervening systems have been reported in the literature. Nevertheless, these discoveries provide very useful information that is complementary to emission lines, allowing, for instance, to distinguish between inflows and outflows. These lines are also good candidates for studying the variations of the fundamental constants of physics. Here we report the findings of our search for CO and other molecules in emission and absorption in front of a sample of 30 targets, comprising 16 associated and 14 intervening H I 21 cm absorbers. The observations were made with the IRAM-30 m telescope simultaneously at 3 mm and 2 mm, exploring several lines of the CO ladder and HCO+, depending on the redshift. We detected eight targets in emission, of which five are new. The derived molecular gas masses range from 109 to 7 × 1011 M⊙ and the highest redshift detection (z = 3.387) corresponds to a relatively average-metallicity damped Lyman-α absorber for this redshift. We also report four new detections in absorption. Two of the associated CO absorption line detections at high-redshift (z = 1.211 and 1.275) result from high-spatial-resolution follow-up observations with NOEMA. The disparity between the mm molecular and H I 21 cm absorption lines for these and another intervening system detected in HNC at z = 1.275 is attributable to radio and mm sight lines tracing different media. We compare the atomic and molecular column densities of 14 known high-redshift (z &gt; 0.1) molecular absorption line systems. The associated H I absorption lines are broad and exhibit multiple components, and the molecular absorption generally corresponds to the broader and weaker 21 cm absorption component. This indicates two distinct phases: one near galaxy centers with a larger CO-to-H I abundance ratio, and another with lower molecular abundance in the outer regions of the galaxy. In comparison, intervening absorption profiles correspond primarily to H I-dominated gas structure in galaxy outskirts, except for gas at low impact parameters in gravitationally lensed systems. The comparison of interferometric and single-dish observations presented here shows that the detection of absorption requires sufficient spatial resolution to overcome the dilution by emission and will be an important criterion for mm follow-up of 21 cm absorbers from ongoing large-scale surveys.

MINDS. JWST/MIRI Reveals a Dynamic Gas-rich Inner Disk inside the Cavity of SY Cha

SY Cha is a T Tauri star surrounded by a protoplanetary disk with a large cavity seen in the millimeter continuum but has the spectral energy distribution of a full disk. Here we report the first results from JWST/Mid-InfraRed Instrument (MIRI) Medium Resolution Spectrometer (MRS) observations taken as part of the MIRI mid-INfrared Disk Survey (MINDS) GTO Program. The much improved resolution and sensitivity of MIRI-MRS compared to Spitzer enables a robust analysis of the previously detected H2O, CO, HCN, and CO2 emission as well as a marginal detection of C2H2. We also report the first robust detection of mid-infrared OH and rovibrational CO emission in this source. The derived molecular column densities reveal the inner disk of SY Cha to be rich in both oxygen- and carbon-bearing molecules. This is in contrast to PDS 70, another protoplanetary disk with a large cavity observed with JWST, which displays much weaker line emission. In the SY Cha disk, the continuum, and potentially the line, flux varies substantially between the new JWST observations and archival Spitzer observations, indicative of a highly dynamic inner disk.

Deuterium fractionation in cold dense cores in the low-mass star-forming region L1688

ABSTRACT In this work, we study deuterium fractionation in four starless cores in the low-mass star-forming region L1688 in the Ophiuchus molecular cloud. We study how the deuterium fraction (RD) changes with environment, compare deuteration of ions and neutrals, core centre and its envelope, and attempt to reproduce the observed results with a gas–grain chemical model. We chose high and low gas density tracers to study both core centre and the envelope. With the IRAM 30 m antenna, we mapped N2H+(1–0), N2D+(1–0), H13CO+ (1–0) and (2–1), DCO+(2–1), and p-NH2D(111–101) towards the chosen cores. The missing p-NH3 and N2H+(1–0) data were taken from the literature. To measure the molecular hydrogen column density, dust and gas temperature within the cores, we used the Herschel/SPIRE dust continuum emission data, the Green Bank Ammonia Survey data (NH3), and the COMPLETE survey data to estimate the upper limit on CO depletion. We present the deuterium fraction maps for three species towards four starless cores. Deuterium fraction of the core envelopes traced by DCO+/H13CO+ is one order of magnitude lower (∼0.08) than that of the core central parts traced by the nitrogen-bearing species (∼0.5). Deuterium fraction increases with the gas density as indicated by high deuterium fraction of high gas density tracers and low deuterium fraction of lower gas density tracers and by the decrease of RD with core radii, consistent with the predictions of the chemical model. Our model results show a good agreement with observations for RD(N2D+/N2H+) and RD(DCO+/HCO+) and underestimate the RD(NH2D/NH3).

HCNH+ abundance in cold dense clouds based on the first hyperfine resolved rate coefficients

The protonated form of hydrogen cyanide, HCNH+, holds significant importance in astrochemistry, serving as an intermediate species in ion-neutral reactions occurring in the cold molecular clouds. Although it plays a crucial role in the chemistry of HCN and HNC, the excitation rate coefficients of this molecular cation by the dominant interstellar colliders have not been thoroughly investigated, leading to limitations in the radiative transfer models used to derive its abundance. In this work, we present the first hyperfine-resolved excitation rate coefficients for HCNH+ induced by collisions with both He and H2 at low temperatures, addressing a crucial requirement for precise modeling of HCNH+ abundance in typical cold dense molecular clouds. Using non-local thermodynamic equilibrium (non-LTE) radiative transfer calculations, we reproduced the 1 → 0 and 2 → 1 observational spectra of HCNH+ fairly well and derived updated molecular column densities. For the TMC-1 molecular cloud, the new HCNH+ abundance is twice as large as suggested by previous LTE modeling, whereas the column density of this molecular cation is improved only by 10% in the case of the L483 proto-star. The factor of two in the case of TMC-1 most likely arises from an error in the early analysis of observational spectra rather than an effect of the LTE assumption, given that the HCNH+ lines are predominantly thermalized at densities higher than 2 × 104 cm−3. For multiline studies of clouds of moderate densities, we strongly recommend using the collisional rate coefficients reported in this work.

A survey of molecular line emission towards Herbig Be star V645 Cyg

ABSTRACT We present a survey of molecular line emission towards the molecular cloud surrounding Herbig Be star V645 Cyg. The survey was performed with the 20-m Onsala space telescope at 3 and 4 mm. We detected emission lines of 33 molecules and their isotopologues from diatomic molecules to four COMs up to seven atoms. Using detected lines, we estimated molecular column densities and abundances relative to molecular hydrogen in local thermodynamic equilibrium (LTE) approximation for all molecules except for methanol, for which we obtained physical parameters using a non-LTE model. Moreover, in the basement of the non-LTE model of methanol line emission, we consider that there is a weak maser effect in the additional spectral component of 51–40 E methanol line at 84.521 GHz. We compared the molecular abundances with values found in several astrochemical templates: molecular clouds, hot cores, and photodissociation regions, and found that signatures of these different types can be found towards V645 Cyg. We also obtained maps of the cloud in several molecular emission lines. The peaks of CO and CH3OH emission are shifted from the direction of the star, but the CS, HCO+, HNC, HCN, and N2H+ emission peaks are observed directly towards the star. Exploring the gas kinematics around V645 Cyg, we found that velocity structure in the ambient molecular cloud on the scale ≈1.6–2.0 pc is not the same as within ≈0.5 pc found previously by other authors.

ALMA-IMF. VIII. Combination of Interferometric Continuum Images with Single-dish Surveys and Structural Analysis of Six Protoclusters

We present the combination of ALMA-IMF and single-dish continuum images from the MUSTANG-2 Galactic Plane Survey (MGPS90) at 3 mm and the Bolocam Galactic Plane Survey (BGPS) at 1 mm. Six and 10 out of the 15 ALMA-IMF fields are combined with MGPS90 and BGPS, respectively. The combination is made via the feathering technique. We used the dendrogram algorithm throughout the combined images, and performed further analysis in the six fields with the combination in both bands (G012.80, W43-MM1, W43-MM2, W43-MM3, W51-E, W51-IRS2). In these fields, we calculated spectral index maps and used them to separate regions dominated by dust or free–free emission, and then performed further structural analysis. We report the basic physical parameters of the dust-dominated (column densities, masses) and ionized (emission measures, hydrogen ionization photon rates) structures. We also searched for multiscale relations in the dust-dominated structures across the analyzed fields, finding that the fraction of mass in dendrogram leaves (which we label leaf mass efficiency (LME)) as a function of molecular gas column density follows a similar trend: a rapid, exponential-like growth, with maximum values approaching 100% in most cases. The observed behavior of the LME with the gas column is tentatively interpreted as an indicator of large star formation activity within the ALMA-IMF protoclusters. W51-E and G012.80 stand out as cases with comparatively large and reduced potential for further star formation, respectively.

Modeling the chemical evolution and kinetics of pure H2O Ices under various types of radiation employing the PROCODA code

Water is one of the most abundant molecules in space, especially in cold environments, where it is the main constituting of astrophysical ices. The space ionizing radiation affects these ices and induces chemical changes, including desorption to gas-phase, which increase the complexity of the interstellar medium. In this work, we employed the PROCODA code to investigate the behavior of several pure water ices under different type of ionizing radiation such as UV, X-rays, electrons and cosmic rays analogues. Here, we employ molecular column densities from laboratory and solved a set of coupled chemical reactions to calculated effective reaction rates (ERCs) and characterize the chemical equilibrium of water ices under high radiation fluences. Briefly, we monitored the evolution of nine species (including the observed ones H2O, H2O2, and O3, and the predicted ones H, O, H2, OH, O2, and HO2). A discussion on the branching ratio for the considered reactions with the type of ionizing radiation is provided. Among the results, we observed that approximately 63% of the modeled molecules quantified at chemical equilibrium were non-observed species in the X-rays experiment, highlighting the importance of this work in providing insights into the processes that occur on the surface of icy interstellar grains exposed to cosmic radiation, including the formation and destruction of water ice. Accurate modeling of these processes can lead to a better understanding of the chemical evolution of interstellar and circumstellar environments, as well as offer insight into the formation and composition of celestial objects such as comets.

The Perseus ALMA Chemical Survey (PEACHES)

Context. Sulfur chemistry is poorly understood in the process of low-mass star and planet formation, where the main carriers of sulfur in both the gas and the dust phase are still unknown. Furthermore, the chemical evolution of sulfur-bearing species is not fully understood given that simple S-bearing molecules, such as SO and SO2, are commonly seen in embedded Class 0/I sources but hardly detected in more evolved Class II disks. Despite the fact that simple S-bearing molecules are usually detected toward embedded sources, large surveys of S-bearing molecules with high angular resolution and sensitive observations are currently lacking. Aims. The goal of this work is to present an unbiased survey of simple sulfur-bearing species in protostars and provide new statistics on detection rates, emitting regions, and molecular column densities. In addition, we investigate the role of S-bearing molecules in accretion processes and the connection between (non-)detection of complex organic molecules (COMs) and S-related species. Methods. We present the observations of sulfur-bearing species (CS, SO,34SO, and SO2) that are part of the Perseus ALMA Chemical Survey (PEACHES). We analyzed a total of 50 Class 0/I sources with observations that have an average angular resolution of about 0″.6 (∼180 au) in ALMA band 6. Results. Class 0 sources show detection rates of 97% for CS, 86% for SO, 31% for 34SO, and 44% for SO2, while Class I sources present detection rates of 71% for CS, 57% for SO, 36% for 34SO, and 43% for SO2. When 34SO is detected, the SO/34SO ratio is lower than the canonical value of 22, suggesting optically thick emission, and the lowest values are found for those sources that are rich in COMs. When SO2 is detected, those sources that show CS and SO emission parallel to the outflow direction are usually very rich in COMs, while for sources where the CS and SO emission is perpendicular to the outflow direction, only a few or no COMs are detected. When CH3OH and SO2 are detected, the comparison between CH3OH and SO2 abundances shows a positive trend and CH3OH is between 10 and 100 times more abundant than SO2. The SO2 abundances toward the PEACHES sample are, on average, two orders of magnitude lower than values from the Ophiuchus star-forming region and comparable with sources in Taurus. Conclusions. The SO/34SO ratio seems to be a good tracer of the inner high-density envelope and it could be used in the future to infer the presence of multiple COMs. The detection of multiple COMs seems to be related to the presence of collimated outflows (seen in CS and SO emission), where a high column density of warm material is expected close to the protostar, and SO2 emission seems to trace the warm gas in those sources where CH3OH is also detected. The difference in SO2 abundances between different star-forming regions might indicate that the sulfur depletion in the gas-phase could depend on the external UV radiation toward the molecular cloud. Finally, the SO2 emission detected in different evolutionary stages seems to arise from different physical mechanisms: high column density of warm material in Class 0 sources, shocks in Class I/II, and exposure to UV radiation from the protostar in more evolved Class II disks.

Variations of the HCO+, HCN, HNC, N2H+, and NH3 deuterium fractionation in high-mass star-forming regions

ABSTRACT We use spectra and maps of the J = 1 − 0 and J = 2 − 1 DCO+, DCN, DNC, $\rm N_2D^+$ lines, and 111−101 ortho- and para-NH2D lines, obtained with the Institut de Radioastronomie Millimétrique (IRAM)-30 m telescope, as well as observations of their hydrogenated isotopologues to study deuteration processes in five high-mass star-forming regions. The temperature was estimated from CH 3CCH lines, also observed with the IRAM-30 m telescope, and from NH 3 lines, observed with the 100 m radio telescope in Effelsberg, as well as using the integrated intensity ratios of the J = 1 − 0 H13CN and HN13C lines and their main isotopologues. Applying a non-local thermodynamic equilibrium radiative transfer model with radex, the gas density and the molecular column densities were estimated. D/H ratios are 0.001–0.05 for DCO+, 0.001–0.02 for DCN, 0.001–0.05 for DNC, and 0.02–0.4 for NH2D. The D/H ratios decrease with increasing temperature in the range of 20–40 K and slightly vary at densities $n(\rm H_2) \sim 10^4\!-\!10^6$ cm−3. The deuterium fraction of $\rm N_2H^{+}$ is 0.008–0.1 at temperatures in the range of 20–25 K and at a density of ∼105 cm−3. We also estimate relative abundances and find ∼10−11–10−9 for DCO+ and DNC, ∼10−11–10−10 for $\rm N_2D^+$, and ∼10−10–10−8 for NH2D. The relative abundances of these species decrease with increasing temperature. However, the DCN/H2 ratio is almost constant (∼10−10). The observational results agree with the predictions of chemical models (although in some cases there are significant differences).

Probing infrared dark clouds with class I methanol masers and thermal molecular emission using the Onsala 20 meter telescope

ABSTRACT Infrared dark clouds (IRDCs) represent the earliest stage of high-mass star formation and host molecular cores at different states of activity from quiet state without any IR-signature of star formation to active state with IR-signatures and pronounced maser activity. Many IRDCs show indications of shocked gas associated with protostar outflows and can contain class I methanol masers (cIMMs). The aim of this study is to probe different types of cores in IRDCs with cIMMs and thermal molecular emission and to investigate the relationship between presence of cIMMs and physical conditions in IRDCs. For a sample of 37 molecular cores, using the 20-m Onsala radio telescope, we collected molecular line data at 44, 85, and 97 GHz for more than 15 species including CH3OH, CH3CCH, and CS. Kinetic temperature of the gas and molecular column densities were obtained. Methanol emission at 44 GHz was detected in 29 sources, with 4 sources are being new discoveries. None of the cores in quiescent state show emission at 44 GHz. Our results testify that cIMM emission is a reliable marker of advanced state of molecular cores. The higher detection rate for intermediate and IR-quiet sources suggests that cIMMs most readily trace the early stages of star formation characterized by moderate IR-signatures. We found that masers with higher flux densities tend to be associated with emission in the CH3CCH lines with higher integrated intensities. Sources undetected at 44 GHz have 4 times lower integrated intensities of CH3CCH and exhibit poorer molecular spectra than the most sources with cIMMs.

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.

Distributions of the Density and Kinetic Temperature of the Molecular Gas in the Central Region of NGC 613 Using Hierarchical Bayesian Inference

We present position–position–velocity (PPV) cubes of the physical and chemical properties of the molecular medium in the central 1.2 kpc region of the active galaxy NGC 613 at a PPV resolution of 0.″8 × 0.″8 × 10 km s−1 (0.″8 = ∼68 pc). We used eight molecular lines (13CO(1–0), C18O(1–0), HCN(1–0), HCO+(1–0), CS(2–1), HCN(4–3), HCO+(4–3), and CS(7–6)) obtained with the Atacama Large Millimeter/submillimeter Array. Non-LTE calculation with hierarchical Bayesian inference was used to construct PPV cubes of the gas kinetic temperature (T kin), molecular hydrogen volume density (), column densities (), and fractional abundances of four molecules (12C18O, HCN, HCO+, and CS). The derived , , and T kin ranged from 103.21−3.85 cm−3, 1020.8−22.1 cm−2, and 102.33−2.64 K, respectively. Our first application of the non-LTE method with the hierarchical Bayesian inference to external galaxies yielded compatible results compared with the previous studies of this galaxy, demonstrating the efficacy of this method for application to other galaxies. We examined the correlation between gas surface density (converted from ) and the star formation rate ΣSFR obtained from the 110 GHz continuum flux map and found two distinct sequences in the –ΣSFR diagram; the southwestern subregion of the star-forming ring exhibited a ∼0.5 dex higher star formation efficiency (SFE; ΣSFR/) than the eastern subregion. However, they exhibited no systematic difference in , which is often argued as a driver of SFE variation. We suggest that the deficiency of molecular gas in the southwestern subregion, where no significant gas supply is evident along the offset ridges in the bar, is responsible for the elevated SFE.

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.