Search

Dense cold molecular cores/filaments are surrounded by an envelope of translucent gas. Some of the low-J emission lines of CO and HCO$^+$ isotopologues are more sensitive to the conditions either in the translucent environment or in the dense cold one. We propose a cloud model composed of three homogeneous slabs of gas along each LoS, representing an envelope and a shielded inner layer. IRAM-30m data from the ORION-B large program toward the Horsehead nebula are used to demonstrate the method's capability. We use the non-LTE radiative transfer code RADEX to model the line profiles from the kinetic temperature $T_{kin}$, the volume density $n_{H_2}$, kinematics and chemical properties of the different layers. We then use a maximum likelihood estimator to simultaneously fit the lines of the CO and HCO$^+$ isotopologues. We constrain column density ratios to limit the variance on the estimates. This simple heterogeneous model provides good fits of the fitted lines over a large part of the cloud. The decomposition of the intensity into three layers allows to discuss the distribution of the estimated physical/chemical properties along the LoS. About 80$\%$ the CO integrated intensity comes from the envelope, while $\sim55\%$ of that of the (1-0) and (2-1) lines of C$^{18}$O comes from the inner layer. The $N(^{13}CO)/N(C^{18}O)$ in the envelope increases with decreasing $A_v$, and reaches $25$ in the pillar outskirts. The envelope $T_{kin}$ varies from 25 to 40 K, that of the inner layer drops to $\sim 15$ K in the western dense core. The inner layer $n_{H_2}$ is $\sim 3\times10^4\,\text{cm}^{-3}$ toward the filament and it increases by a factor $10$ toward dense cores. The proposed method correctly retrieves the physical/chemical properties of the Horsehead nebula and offers promising prospects for less supervised model fits of wider-field datasets.

Observations of ionic, atomic, or molecular lines are performed to improve our understanding of the interstellar medium (ISM). However, the potential of a line to constrain the physical conditions of the ISM is difficult to assess quantitatively, because of the complexity of the ISM physics. The situation is even more complex when trying to assess which combinations of lines are the most useful. Therefore, observation campaigns usually try to observe as many lines as possible for as much time as possible. We search for a quantitative statistical criterion to evaluate the constraining power of a (or combination of) tracer(s) with respect to physical conditions in order to improve our understanding of the statistical relationships between ISM tracers and physical conditions and helps observers to motivate their observation proposals. The best tracers are obtained by comparing the mutual information between a physical parameter and different sets of lines. We apply this method to simulations of radio molecular lines emitted by a photodissociation region similar to the Horsehead Nebula that would be observed at the IRAM 30m telescope. We search for the best lines to constrain the visual extinction $A_v^{tot}$ or the far UV illumination $G_0$. The most informative lines change with the physical regime (e.g., cloud extinction). Short integration time of the CO isotopologue $J=1-0$ lines already yields much information on the total column density most regimes. The best set of lines to constrain the visual extinction does not necessarily combine the most informative individual lines. Precise constraints on $G_0$ are more difficult to achieve with molecular lines. They require spectral lines emitted at the cloud surface (e.g., [CII] and [CI] lines). This approach allows one to better explore the knowledge provided by ISM codes, and to guide future observation campaigns.

The study of the chemical composition of the interstellar medium (ISM) requires a strong synergy between laboratory astrophysics, modeling, and observations. In particular, astrochemical models have been developed for decades now and include an increasing number of processes studied in the laboratory or theoretically. These models follow the chemistry both in the gas phase and at the surface of interstellar grains. Since 2012, we have provided complete gas-phase chemical networks for astrochemical codes that can be used to model various environments of the ISM. Our aim is to introduce the new up-to-date astrochemical network kida.uva.2024 together with the ice chemical network and the fortran code to compute time dependent compositions of the gas, the ice surface, and the ice mantles under physical conditions relevant for the ISM. The gas-phase chemical reactions, as well as associated rate coefficients, included in kida.uva.2024 were carefully selected from the KIDA online database and represent the most recent values. The model predictions for cold core conditions and for when considering only gas-phase processes were computed as a function of time and compared to the predictions obtained with the previous version, kida.uva.2014. In addition, key chemical reactions were identified. The model predictions, including both gas and surface processes, were compared to the molecular abundances as observed in the cold core TMC1-CP. Many gas-phase reactions were revised or added to produce kida.uva.2024. The new model predictions are different by several orders of magnitude for some species. The agreement of this new model with observations in TMC-1 (CP) is, however, similar to the one obtained with the previous network., Comment: Accepted for publication in A&A

Resolved molecular line observations are essential for gaining insight into the physical and chemical structure of protoplanetary disks, particularly in cold, dense regions where planets form and acquire their chemical compositions. However, tracing these regions is challenging because most molecules freeze onto grain surfaces and are not observable in the gas phase. We investigated cold molecular chemistry in the triple stellar T Tauri disk GG Tau A, which harbours a massive gas and dust ring and an outer disk, using ALMA Band 7 observations. We present high angular resolution maps of N2H+ and DCO+ emission, with upper limits reported for H2D+, 13CS, and SO2. The radial intensity profile of N2H+ shows most emission near the ring outer edge, while DCO+ exhibits double peaks, one near the ring inner edge and the other in the outer disk. With complementary observations of lower-lying transitions, we constrained the molecular surface densities and rotation temperatures. We compared the derived quantities with model predictions across different cosmic ray ionization (CRI) rates, carbon-to-oxygen (C/O) ratios, and stellar UV fluxes. Cold molecular chemistry, affecting N2H+, DCO+, and H2D+ abundances, is most sensitive to CRI rates, while stellar UV flux and C/O ratios have minimal impact on these three ions. Our best model requires a low cosmic ray ionization rate of 1e-18 s-1. However, it fails to match the low temperatures derived from N2H+ and DCO+, 12 to 16 K, which are much lower than the CO freezing temperature., Comment: Accepted for publication in the Astrophysical Journal; 26 Pages (15 figures and 8 tables)

Robust radiative transfer techniques are requisite for efficiently extracting the physical and chemical information from molecular rotational lines.We study several hypotheses that enable robust estimations of the column densities and physical conditions when fitting one or two transitions per molecular species. We study the extent to which simplifying assumptions aimed at reducing the complexity of the problem introduce estimation biases and how to detect them.We focus on the CO and HCO+ isotopologues and analyze maps of a 50 square arcminutes field. We used the RADEX escape probability model to solve the statistical equilibrium equations and compute the emerging line profiles, assuming that all species coexist. Depending on the considered set of species, we also fixed the abundance ratio between some species and explored different values. We proposed a maximum likelihood estimator to infer the physical conditions and considered the effect of both the thermal noise and calibration uncertainty. We analyzed any potential biases induced by model misspecifications by comparing the results on the actual data for several sets of species and confirmed with Monte Carlo simulations. The variance of the estimations and the efficiency of the estimator were studied based on the Cram{\'e}r-Rao lower bound.Column densities can be estimated with 30% accuracy, while the best estimations of the volume density are found to be within a factor of two. Under the chosen model framework, the peak 12CO(1--0) is useful for constraining the kinetic temperature. The thermal pressure is better and more robustly estimated than the volume density and kinetic temperature separately. Analyzing CO and HCO+ isotopologues and fitting the full line profile are recommended practices with respect to detecting possible biases.Combining a non-local thermodynamic equilibrium model with a rigorous analysis of the accuracy allows us to obtain an efficient estimator and identify where the model is misspecified. We note that other combinations of molecular lines could be studied in the future., Comment: Astronomy and Astrophysics - A\&A, In press

We present 5 deg^2 (~250 pc^2) HCN, HNC, HCO+, and CO J=1-0 maps of the Orion B GMC, complemented with existing wide-field [CI] 492 GHz maps, as well as new pointed observations of rotationally excited HCN, HNC, H13CN, and HN13C lines. We detect anomalous HCN J=1-0 hyperfine structure line emission almost everywhere in the cloud. About 70% of the total HCN J=1-0 luminosity arises from gas at A_V < 8 mag. The HCN/CO J=1-0 line intensity ratio shows a bimodal behavior with an inflection point at A_V < 3 mag typical of translucent gas and UV-illuminated cloud edges. We find that most of the HCN J=1-0 emission arises from extended gas with n(H2) ~< 10^4 cm^-3, even lower density gas if the ionization fraction is > 10^-5 and electron excitation dominates. This result explains the low-A_V branch of the HCN/CO J=1-0 intensity ratio distribution. Indeed, the highest HCN/CO ratios (~0.1) at A_V < 3 mag correspond to regions of high [CI] 492 GHz/CO J=1-0 intensity ratios (>1) characteristic of low-density PDRs. Enhanced FUV radiation favors the formation and excitation of HCN on large scales, not only in dense star-forming clumps. The low surface brightness HCN and HCO+ J=1-0 emission scale with I_FIR (a proxy of the stellar FUV radiation field) in a similar way. Together with CO J=1-0, these lines respond to increasing I_FIR up to G0~20. On the other hand, the bright HCN J=1-0 emission from dense gas in star-forming clumps weakly responds to I_FIR once the FUV radiation field becomes too intense (G0>1500). The different power law scalings (produced by different chemistries, densities, and line excitation regimes) in a single but spatially resolved GMC resemble the variety of Kennicutt-Schmidt law indexes found in galaxy averages. As a corollary for extragalactic studies, we conclude that high HCN/CO J=1-0 line intensity ratios do not always imply the presence of dense gas., Comment: accepted for publication in A&A. 24 pages, 18 figures, plus Appendix. Abridged Abstract

The interpretation of observations of atomic and molecular tracers in the galactic and extragalactic interstellar medium (ISM) requires comparisons with state-of-the-art astrophysical models to infer some physical conditions. Usually, ISM models are too time-consuming for such inference procedures, as they call for numerous model evaluations. As a result, they are often replaced by an interpolation of a grid of precomputed models. We propose a new general method to derive faster, lighter, and more accurate approximations of the model from a grid of precomputed models. These emulators are defined with artificial neural networks (ANNs) designed and trained to address the specificities inherent in ISM models. Indeed, such models often predict many observables (e.g., line intensities) from just a few input physical parameters and can yield outliers due to numerical instabilities or physical bistabilities. We propose applying five strategies to address these characteristics: 1) an outlier removal procedure; 2) a clustering method that yields homogeneous subsets of lines that are simpler to predict with different ANNs; 3) a dimension reduction technique that enables to adequately size the network architecture; 4) the physical inputs are augmented with a polynomial transform to ease the learning of nonlinearities; and 5) a dense architecture to ease the learning of simple relations. We compare the proposed ANNs with standard classes of interpolation methods to emulate the Meudon PDR code, a representative ISM numerical model. Combinations of the proposed strategies outperform all interpolation methods by a factor of 2 on the average error, reaching 4.5% on the Meudon PDR code. These networks are also 1000 times faster than accurate interpolation methods and require ten to forty times less memory. This work will enable efficient inferences on wide-field multiline observations of the ISM.

Context. The availability of large bandwidth receivers for millimeter radio telescopes allows the acquisition of position-position-frequency data cubes over a wide field of view and a broad frequency coverage. These cubes contain much information on the physical, chemical, and kinematical properties of the emitting gas. However, their large size coupled with inhomogenous signal-to-noise ratio (SNR) are major challenges for consistent analysis and interpretation.Aims. We search for a denoising method of the low SNR regions of the studied data cubes that would allow to recover the low SNR emission without distorting the signals with high SNR.Methods. We perform an in-depth data analysis of the 13 CO and C 17 O (1 -- 0) data cubes obtained as part of the ORION-B large program performed at the IRAM 30m telescope. We analyse the statistical properties of the noise and the evolution of the correlation of the signal in a given frequency channel with that of the adjacent channels. This allows us to propose significant improvements of typical autoassociative neural networks, often used to denoise hyperspectral Earth remote sensing data. Applying this method to the 13 CO (1 -- 0) cube, we compare the denoised data with those derived with the multiple Gaussian fitting algorithm ROHSA, considered as the state of the art procedure for data line cubes.Results. The nature of astronomical spectral data cubes is distinct from that of the hyperspectral data usually studied in the Earth remote sensing literature because the observed intensities become statistically independent beyond a short channel separation. This lack of redundancy in data has led us to adapt the method, notably by taking into account the sparsity of the signal along the spectral axis. The application of the proposed algorithm leads to an increase of the SNR in voxels with weak signal, while preserving the spectral shape of the data in high SNR voxels.Conclusions. The proposed algorithm that combines a detailed analysis of the noise statistics with an innovative autoencoder architecture is a promising path to denoise radio-astronomy line data cubes. In the future, exploring whether a better use of the spatial correlations of the noise may further improve the denoising performances seems a promising avenue. In addition

(Abridged) Massive stars can determine the evolution of molecular clouds with their strong ultraviolet (UV) radiation fields. Moreover, UV radiation is relevant in setting the thermal gas pressure in star-forming clouds, whose influence can extend from the rims of molecular clouds to entire star-forming galaxies. Probing the fundamental structure of nearby molecular clouds is therefore crucial to understand how massive stars shape their surrounding medium and how fast molecular clouds are destroyed, specifically at their UV-illuminated edges, where models predict an intermediate zone of neutral atomic gas between the molecular cloud and the surrounding ionized gas whose size is directly related to the exposed physical conditions. We present the highest angular resolution (~$0.5$", corresponding to $207$ au) and velocity-resolved images of the molecular gas emission in the Horsehead nebula, using CO J=3-2 and HCO$^+$ J=4-3 observations with ALMA. We find that CO and HCO$^+$ are present at the edge of the cloud, very close to the ionization (H$^+$/H) and dissociation fronts (H/H$_2$), suggesting a very thin layer of neutral atomic gas (<$650$ au) and a small amount of CO-dark gas ($A_V=0.006-0.26$ mag) for stellar UV illumination conditions typical of molecular clouds in the Milky Way. The new ALMA observations reveal a web of molecular gas filaments with an estimated thermal gas pressure of $P_{\mathrm{th}} = (2.3 - 4.0) \times 10^6$ K cm$^{-3}$, and the presence of a steep density gradient at the cloud edge that can be well explained by stationary isobaric PDR models with pressures consistent with our estimations. However, in the HII region and PDR interface, we find $P_{\mathrm{th,PDR}} > P_{\mathrm{th,HII}}$, suggesting the gas is slightly compressed. Therefore, dynamical effects cannot be completely ruled out and even higher angular observations will be needed to unveil their role., Comment: 15 pages, 1 table, 9 figures; Accepted for publication in A&A

Ice is ubiquitous in the interstellar medium. We model the formation of the main constituents of interstellar ices, including H2O, CO2 , CO, and CH3 OH. We strive to understand what physical or chemical parameters influence the final composition of the ice and how they benchmark to what has already been observed, with the aim of applying these models to the preparation and analysis of JWST observations. We used the Nautilus gas-grain model, which computes the gas and ice composition as a function of time for a set of physical conditions, starting from an initial gas phase composition. All important processes (gas-phase reactions, gas-grain interactions, and grain surface processes) are included and solved with the rate equation approximation. We first ran an astrochemical code for fixed conditions of temperature and density mapped in the cold core L429-C to benchmark the chemistry. One key parameter was revealed to be the dust temperature. When the dust temperature is higher than 12 K, CO2 will form efficiently at the expense of H2O, while at temperatures below 12 K, it will not form. Whatever hypothesis we assumed for the chemistry (within realistic conditions), the static simulations failed to reproduce the observed trends of interstellar ices in our target core. In a second step, we simulated the chemical evolution of parcels of gas undergoing different physical and chemical situations throughout the molecular cloud evolution and starting a few 1e7 yr prior to the core formation (dynamical simulations). Our dynamical simulations satisfactorily reproduce the main trends already observed for interstellar ices. Moreover, we predict that the apparent constant ratio of CO2/H2O observed to date is probably not true for regions of low AV , and that the history of the evolution of clouds plays an essential role, even prior to their formation., Comment: Accepted for publication in A&A

Cold cores are an early step of star formation, characterized by densities > 10$^4$ cm$^{-3}$, low temperatures (< 15 K), and very low external UV radiation. We investigate the physico-chemical processes at play to tracing the origin of molecules that are predominantly formed via reactions on dust grain surfaces. We observed the cold core LDN 429-C with the NOEMA interferometer and the IRAM 30m single dish telescope in order to obtain the gas-phase abundances of key species, including CO and CH$_3$OH. Comparing the observed gas phase of methanol to its solid phase previously observed with Spitzer allows us to put quantitative constraints on the efficiency of the non-thermal desorption of this species. With physical parameters determined from available Herschel data, we computed abundance maps of 11 detected molecules with a non-local thermal equilibrium radiative transfer model. These observations allowed us to probe the molecular abundances as a function of density and visual extinction, with the variation in temperature being restrained between 12 and 18 K. We then compared the observed abundances to the predictions of the Nautilus astrochemical model. We find that all molecules have lower abundances at high densities and visual extinctions with respect to lower density regions, except for methanol. Comparing these observations with a grid of chemical models based on the local physical conditions, we were able to reproduce these observations, allowing only the parameter time to vary. Comparing the observed gas-phase abundance of methanol with previous measurements of the methanol ice, we estimate a non-thermal desorption efficiency between 0.002% and 0.09%, increasing with density. The apparent increase in the desorption efficiency cannot be reproduced by our model unless the yield of cosmic-ray sputtering is altered due to the ice composition varying as a function of density.

Understanding the initial properties of star-forming material and how they affect the star formation process is key. From an observational point of view, the feedback from young high-mass stars on future star formation properties is still poorly constrained. In the framework of the IRAM 30m ORION-B large program, we obtained observations of the translucent and moderately dense gas, which we used to analyze the kinematics over a field of 5 deg^2 around the filamentary structures. We used the ROHSA algorithm to decompose and de-noise the C18O(1-0) and 13CO(1-0) signals by taking the spatial coherence of the emission into account. We produced gas column density and mean velocity maps to estimate the relative orientation of their spatial gradients. We identified three cloud velocity layers at different systemic velocities and extracted the filaments in each velocity layer. The filaments are preferentially located in regions of low centroid velocity gradients. By comparing the relative orientation between the column density and velocity gradients of each layer from the ORION-B observations and synthetic observations from 3D kinematic toy models, we distinguish two types of behavior in the dynamics around filaments: (i) radial flows perpendicular to the filament axis that can be either inflows (increasing the filament mass) or outflows and (ii) longitudinal flows along the filament axis. The former case is seen in the Orion B data, while the latter is not identified. We have also identified asymmetrical flow patterns, usually associated with filaments located at the edge of an HII region. This is the first observational study to highlight feedback from HII regions on filament formation and, thus, on star formation in the Orion B cloud. This simple statistical method can be used for any molecular cloud to obtain coherent information on the kinematics., Comment: 45 pages. Abridged abstract. Accepted by Astronomy and Astrophysics

The observational counterparts of theoretically predicted first hydrostatic cores (FHSC) have been searched for in the interstellar medium for nearly two decades now. Distinguishing them from other types of more evolved but still embedded objects remains a challenge because these objects have a short lifetime, are small, and embedded in a dense cocoon. One possible lead to finding them is the characterization of the outflows that are launched by these objects. We observed the L1451-mm FHSC candidate with the NOEMA interferometer in order to study the emission of several molecules. A nonlocal thermodynamic equilibrium analysis of the CH$_3$OH detected lines was performed to retrieve the physical conditions of the emitting region around the central source, together with the CH$_3$OH, SiO, CS, and H$_2$CO column densities. Of the targeted molecules, we detected lines of c-C$_3$H$_2$, CH$_3$OH, CS, C$^{34}$S, SO, DCN, DCO$^+$, H$_2$CO, HC$_3$N, HDCO, and SiO. One of the methanol lines appears to be a maser line. The detection of this class I maser and the SiO line in L1451-mm support the presence of a low-velocity and compact outflow. The excitation conditions of the thermal lines of methanol are also compatible with shocks (H$_2$ density of $\sim 3\times 10^6$~cm$^{-3}$ and a temperature higher than 40~K). Although these low-velocity outflows are theoretically predicted by some models of FHSC, these models also predict the shock temperature to be below 20~K, that is, not evaporating methanol. In addition, the predicted velocities would not erode the grains and release silicon in the gas phase. We therefore conclude that these new observations favor the hypothesis that L1451-mm would be at a very early protostellar stage, launching an outflow nearly on the plane of the sky with a higher velocity than is observed., Comment: Accepted for publication in A&A

The accretion rates needed to fuel the central black hole in a galaxy can be achieved via viscous torques in thick disks and rings, which can be resolved by millimetre interferometry within the inner ~20pc of the active galaxy NGC1068 at comparable scales and sensitivity to single dish observations of the Circumnuclear Disk (CND) in the Galactic Center. To interpret observations of these regions and determine the physical properties of their gas distribution, we present a modelling effort that includes (i) a simple dynamical simulations involving partially inelastic collisions between disk gas clouds, (ii) an analytical model of a turbulent clumpy gas disk calibrated by the dynamical model and observations, (iii) local turbulent and cosmic ray gas heating and cooling via H2O, H2, and CO emission, and (iv) determination of the molecular abundances. We also consider photodissociation regions (PDR) where gas is directly illuminated by the central engine. We compare the resulting model datacubes of the CO, HCN, HCO+, and CS brightness temperatures to available observations. In both cases the kinematics can be explained by one or two clouds colliding with a pre-existing ring, in a prograde sense for the CND and retrograde for NGC1068. And, with only dense disk clouds, the line fluxes can be reproduced to within a factor of about two. To avoid self-absorption of the intercloud medium, turbulent heating at the largest scales, comparable to the disk height, has to be decreased by a factor of 50-200. Our models indicate that turbulent mechanical energy input is the dominant gas heating mechanism within the thick gas disks. In N1068, while the bulk of the AGN X-ray radiation is absorbed in a layer of Compton-thick gas inside the dust sublimation radius, the optical/UV radiation may enhance the molecular line emission from photodissociation regions by ~50% at the inner edge of the gas ring., Comment: 36 pages, 32 figures, accepted for publication in A&A

The ISM is a turbulent, multi-phase, and multi-scale medium following scaling relations. Analytical models of galactic gaseous disks need to take into account the multi-scale and multi-phase nature of the interstellar medium. They can be described as clumpy star-forming accretion disks in vertical hydrostatic equilibrium, with the mid-plane pressure balancing the gravity of the gaseous and stellar disk. ISM turbulence is taken into account by applying Galactic scaling relations to the cold atomic and molecular gas phases. Turbulence is maintained through energy injection by supernovae. With the determination of the gas mass fraction at a given spatial scale, the equilibrium gas temperature between turbulent heating and line cooling, the molecular abundances, and the molecular line emission can be calculated. The resulting model radial profiles of IR, H{\sc i}, CO, HCN, and HCO$^+$ emission are compared to THINGS, HERACLES, EMPIRE, SINGS, and GALEX observations of 17 local spiral galaxies. The Toomre parameter, which measures the stability against star formation (cloud collapse), exceeds unity in the inner disk of a significant number of galaxies. In two galaxies it also exceeds unity in the outer disk. Therefore, in spiral galaxies $Q_{\rm tot}=1$ is not ubiquitous. The model gas velocity dispersion is consistent with the observed H{\sc i} velocity dispersion where available. Within our model HCN and HCO$^+$ is already detectable in relatively low-density gas ($\sim 1000$~cm$^{-3}$). CO and HCN conversion factors and molecular gas depletion time were derived. Both conversion factors are consistent with values found in the literature. Whereas in the massive galaxies the viscous timescale greatly exceeds the star formation timescale, the viscous timescale is smaller than the star formation timescale within $\rm{R}~\sim~2~\rm{R}_{\rm d}$, the disk scale length, in the low-mass galaxies., Comment: 22 pages, 10 figures (without appendix) 44 pages, 27 figures