First ALMA maps of HCO, an important precursor of complex organic molecules, towards IRAS 16293-2422

The formyl radical HCO has been proposed as the basic precursor of many complex organic molecules such as methanol (CH3OH) and glycolaldehyde (CH2OHCHO). Using ALMA, we have mapped, for the first time at high angular resolution (∼1 arcsec, ∼140 au), HCO towards the solar-type protostellar binary IRAS 16293–2422, where numerous complex organic molecules have been previously detected. We also detected several lines of the chemically related species H2CO, CH3OH, and CH2OHCHO. The observations revealed compact HCO emission arising from the two protostars. The line profiles also show redshifted absorption produced by foreground material of the circumbinary envelope that is infalling towards the protostars. Additionally, IRAM 30 m single-dish data revealed a more extended HCO component arising from the common circumbinary envelope. The comparison between the observed molecular abundances and our chemical model suggests that whereas the extended HCO from the envelope can be formed via gas-phase reactions during the cold collapse of the natal core, the HCO in the hot corinos surrounding the protostars is predominantly formed by the hydrogenation of CO on the surface of dust grains and subsequent thermal desorption during the protostellar phase. The derived abundance of HCO in the dust grains is high enough to produce efficiently more complex species such as H2CO, CH3OH, and CH2OHCHO by surface chemistry. We found that the main formation route of CH2OHCHO is the reaction between HCO and CH2OH

The GUAPOS project: III. Characterization of the O- and N-bearing complex organic molecules content and search for chemical differentiation

Context. The G31.41+0.31 Unbiased ALMA sPectral Observational Survey (GUAPOS) project targets the hot molecular core (HMC) G31.41+0.31 (G31) to reveal the complex chemistry of one of the most chemically rich high-mass star-forming regions outside the Galactic center (GC). Aims. In the third paper of the project we present a study of nine O-bearing (CH3OH, 13CH3OH, CH318OH, CH3CHO, CH3OCH3, CH3COCH3, C2H5OH, aGg′-(CH2OH)2, and gGg′-(CH2OH)2) and six N-bearing (CH3CN,13CH3CN, CH313CN, C2H3CN, C2H5CN, and C2H513CN) complex organic molecules toward G31. The aim of this work is to characterize the abundances in G31 and to compare them with the values estimated in other sources. Moreover, we searched for a possible chemical segregation between O-bearing and N-bearing species in G31, which hosts four compact sources as seen with higher angular resolution data. In the discussion we also include the three isomers of C2H4O2 and the O- and N-bearing molecular species NH2CHO, CH3NCO, CH3C(O)CH2, and CH3NHCHO, which were analyzed in previous GUAPOS papers. Methods. Observations were carried out with the interferometer ALMA and cover the entire Band 3 from 84 to 116 GHz (∼32 GHz bandwidth) with an angular resolution of 1.2″ × 1.2″ (∼4400 au × 4400 au) and a spectral resolution of ∼0.488 MHz (∼1.3- 1.7 km s-1). The transitions of the 14 molecular species were analyzed with the tool SLIM of MADCUBA to determine the physical parameters of the emitting gas. Moreover, we analyzed the morphology of the emission of the molecular species. Results. The values of abundances with respect to H2 in G31 range from 10-6 to 10-10 for the different species. We compared the abundances with respect to methanol of O-bearing, N-bearing, and O- and N-bearing COMs in G31 with 27 other sources, including other hot molecular cores inside and outside the GC, hot corinos, shocked regions, envelopes around young stellar objects, and quiescent molecular clouds, and with chemical models. Conclusions. From the comparison with other sources there is not a unique template for the abundances in hot molecular cores, pointing toward the importance of the thermal history for the chemistry of the various sources. The abundances derived from the chemical models are in good agreement, within a factor of 10, with those of G31. From the analysis of the maps we derived the peak positions of all the molecular species toward G31. Different species peak at slightly different positions, and this, together with the different central velocities of the lines obtained from the spectral fitting, point to chemical differentiation of selected O-bearing species

Towards the prediction of molecular parameters from astronomical emission lines using Neural Networks

Molecular astronomy is a field that is blooming in the era of large observatories such as the Atacama Large Millimeter/Submillimeter Array (ALMA). With modern, sensitive, and high spectral resolution radio telescopes like ALMA and the Square Kilometer Array, the size of the data cubes is rapidly escalating, generating a need for powerful automatic analysis tools. This work introduces MolPred, a pilot study to perform predictions of molecular parameters such as excitation temperature (Tex) and column density (log(N)) from input spectra by the use of neural networks. We used as test cases the spectra of CO, HCO+, SiO and CH3CN between 80 and 400 GHz. Training spectra were generated with MADCUBA, a state-of-the-art spectral analysis tool. Our algorithm was designed to allow the generation of predictions for multiple molecules in parallel. Using neural networks, we can predict the column density and excitation temperature of these molecules with a mean absolute error of 8.5% for CO, 4.1% for HCO+, 1.5% for SiO and 1.6% for CH3CN. The prediction accuracy depends on the noise level, line saturation, and number of transitions. We performed predictions upon real ALMA data. The values predicted by our neural network for this real data differ by 13% from the MADCUBA values on average. Current limitations of our tool include not considering linewidth, source size, multiple velocity components, and line blending

The distribution and origin of C₂H in NGC 253 from ALCHEMI

Context. Observations of chemical species can provide insights into the physical conditions of the emitting gas however it is important to understand how their abundances and excitation vary within different heating environments. C2H is a molecule typically found in PDR regions of our own Galaxy but there is evidence to suggest it also traces other regions undergoing energetic processing in extragalactic environments. / Aims. As part of the ALCHEMI ALMA large program, we map the emission of C2H in the central molecular zone of the nearby starburst galaxy NGC 253 at 1.6″ (28 pc) resolution and characterize it to understand its chemical origins. / Methods. We used spectral modeling of the N = 1−0 through N = 4−3 rotational transitions of C2H to derive the C2H column densities towards the dense clouds in NGC 253. We then use chemical modeling, including photodissociation region (PDR), dense cloud, and shock models to investigate the chemical processes and physical conditions that are producing the molecular emission. / Results. We find high C2H column densities of ∼1015 cm−2 detected towards the dense regions of NGC 253. We further find that these column densities cannot be reproduced if it is assumed that the emission arises from the PDR regions at the edge of the clouds. Instead, we find that the C2H abundance remains high even in the high visual extinction interior of these clouds and that this is most likely caused by a high cosmic-ray ionization rate

Starburst Energy Feedback Seen through HCO+/HOC+ Emission in NGC 253 from ALCHEMI

Molecular abundances are sensitive to the UV photon flux and cosmic-ray ionization rate. In starburst environments, the effects of high-energy photons and particles are expected to be stronger. We examine these astrochemical signatures through multiple transitions of HCO+ and its metastable isomer HOC+ in the center of the starburst galaxy NGC 253 using data from the Atacama Large Millimeter/submillimeter Array large program ALMA Comprehensive High-resolution Extragalactic Molecular inventory. The distribution of the HOC+(1−0) integrated intensity shows its association with "superbubbles," cavities created either by supernovae or expanding H ii regions. The observed HCO+/HOC+ abundance ratios are ∼10–150, and the fractional abundance of HOC+ relative to H2 is ∼1.5 × 10−11–6 × 10−10, which implies that the HOC+ abundance in the center of NGC 253 is significantly higher than in quiescent spiral arm dark clouds in the Galaxy and the Galactic center clouds. Comparison with chemical models implies either an interstellar radiation field of G0 ≳ 103 if the maximum visual extinction is ≳5, or a cosmic-ray ionization rate of ζ ≳ 10−14 s−1 (3–4 orders of magnitude higher than that within clouds in the Galactic spiral arms) to reproduce the observed results. From the difference in formation routes of HOC+, we propose that a low-excitation line of HOC+ traces cosmic-ray dominated regions, while high-excitation lines trace photodissociation regions. Our results suggest that the interstellar medium in the center of NGC 253 is significantly affected by energy input from UV photons and cosmic rays, sources of energy feedback

Reconstructing the shock history in the CMZ of NGC 253 with ALCHEMI

Context: HNCO and SiO are well-known shock tracers and have been observed in nearby galaxies, including the nearby (D = 3.5 Mpc) starburst galaxy NGC 253. The simultaneous detection of these two species in regions where the star-formation rate is high may be used to study the shock history of the gas. // Aims: We perform a multi-line molecular study of NGC 253 using the shock tracers SiO and HNCO and aim to characterize its gas properties. We also explore the possibility of reconstructing the shock history in the central molecular zone (CMZ) of the galaxy. // Methods: Six SiO transitions and eleven HNCO transitions were imaged at high resolution 1.″6 (28 pc) with the Atacama Large Millimeter/submillimeter Array (ALMA) as part of the ALCHEMI Large Programme. Both non local thermaldynamic equilibrium (non-LTE) radiative transfer analysis and chemical modeling were performed in order to characterize the gas properties and investigate the chemical origin of the emission. // Results: The nonLTE radiative transfer analysis coupled with Bayesian inference shows clear evidence that the gas traced by SiO has different densities and temperatures than that traced by HNCO, with an indication that shocks are needed to produce both species. Chemical modeling further confirms such a scenario and suggests that fast and slow shocks are responsible for SiO and HNCO production, respectively, in most GMCs. We are also able to infer the physical characteristics of the shocks traced by SiO and HNCO for each GMC. // Conclusions: Radiative transfer and chemical analysis of the SiO and HNCO in the CMZ of NGC 253 reveal a complex picture whereby most of the GMCs are subjected to shocks. We speculate on the possible shock scenarios responsible for the observed emission and provide potential history and timescales for each shock scenario. Observations of higher spatial resolution for these two species are required in order to quantitatively differentiate between the possible scenarios

The feedback of an HC HII region on its parental molecular core: The case of core A1 in the star-forming region G24.78+0.08

Context: G24.78+0.08 is a well known high-mass star-forming region, where several molecular cores harboring OB young stellar objects are found inside a clump of size ≈1 pc. This article focuses on the most prominent of these cores, A1, where an intense hypercompact (HC) HII region has been discovered by previous observations. Aims: Our aim is to determine the physical conditions and the kinematics of core A1, and study the interaction of the HII region with the parental molecular core. Methods: We combine ALMA 1.4 mm high-angular resolution (≈0.′′2) observations of continuum and line emission with multi-epoch Very Long Baseline Interferometry data of water 22 GHz and methanol 6.7 GHz masers. These observations allow us to study the gas kinematics on linear scales from 10 to 104 au, and to accurately map the physical conditions of the gas over core A1. Results: The 1.4 mm continuum is dominated by free-free emission from the intense HC HII region (size ≈1000 au) observed to the North of core A1 (region A1N). Analyzing the H30α line, we reveal a fast bipolar flow in the ionized gas, covering a range of LSR velocities (VLSR) of ≈60 km s−1. The amplitude of the VLSR gradient, 22 km s−1 mpc−1, is one of the highest so far observed towards HC HII regions. Water and methanol masers are distributed around the HC HII region in A1N, and the maser three-dimensional (3D) velocities clearly indicate that the ionized gas is expanding at high speed (≥200 km s−1) into the surrounding molecular gas. The temperature distribution (in the range 100–400 K) over core A1, traced with molecular (CH3OCHO, 13CH3CN, 13CH3OH, and CH3CH2CN) transitions with level energy in the range 30 K ≤ Eu/k ≤ 300 K, reflects the distribution of shocks produced by the fast-expansion of the ionized gas of the HII region. The high-energy (550 K ≤ Eu/k ≤ 800 K) transitions of vibrationally excited CH3CN are likely radiatively pumped, and their rotational temperature can significantly differ from the kinetic temperature of the gas. Over core A1, the VLSR maps from both the 1.4 mm molecular lines and the 6.7 GHz methanol masers consistently show a VLSR gradient (amplitude ≈0.3 km s−1 mpc−1) directed approximately S–N. Rather than gravitationally supported rotation of a massive toroid, we interpret this velocity gradient as a relatively slow expansion of core A1

Kinematics of Galactic Centre clouds shaped by shear-seeded solenoidal turbulence

The Central Molecular Zone (CMZ; the central ∼500 pc of the Galaxy) is a kinematically unusual environment relative to the Galactic disc, with high velocity dispersions and a steep size-linewidth relation of the molecular clouds. In addition, the CMZ region has a significantly lower star formation rate (SFR) than expected by its large amount of dense gas. An important factor in explaining the low SFR is the turbulent state of the star-forming gas, which seems to be dominated by rotational modes. However, the turbulence driving mechanism remains unclear. In this work, we investigate how the Galactic gravitational potential affects the turbulence in CMZ clouds. We focus on the CMZ cloud G0.253+0.016 (‘the Brick’), which is very quiescent and unlikely to be kinematically dominated by stellar feedback. We demonstrate that several kinematic properties of the Brick arise naturally in a cloud-scale hydrodynamics simulation that takes into account the Galactic gravitational potential. These properties include the line-of-sight velocity distribution, the steepened size-linewidth relation, and the predominantly solenoidal nature of the turbulence. Within the simulation, these properties result from the Galactic shear in combination with the cloud’s gravitational collapse. This is a strong indication that the Galactic gravitational potential plays a crucial role in shaping the CMZ gas kinematics, and is a major contributor to suppressing the SFR by inducing predominantly solenoidal turbulent modes

Substructures in the Keplerian disc around the O-type (proto-)star G17.64+0.16

Contains fulltext : 205611.pdf (publisher's version ) (Open Access

Accelerating infall and rotational spin-up in the hot molecular core G31.41+0.31

As part of our effort to search for circumstellar disks around high-mass stellar objects, we observed the well-known core G31.41 +0.31 with ALMA at 1.4 mm with an angular resolution of ~0.′′22 (~1700 au). The dust continuum emission has been resolved into two cores namely Main and NE. The Main core, which has the stronger emission and is the more chemically rich, has a diameter of ~5300 au, and is associated with two free-free continuum sources. The Main core looks featureless and homogeneous in dust continuum emission and does not present any hint of fragmentation. Each transition of CH₃CN and CH₃OCHO, both ground and vibrationally excited, as well as those of CH₃CN isotopologues, shows a clear velocity gradient along the NE–SW direction, with velocity linearly increasing with distance from the center, consistent with solid-body rotation. However, when comparing the velocity field of transitions with different upper level energies, the rotation velocity increases with increasing energy of the transition, which suggests that the rotation speeds up toward the center. Spectral lines towardtoward the dust continuum peak show an inverse P-Cygni profile that supports the existence of infall in the core. The infall velocity increases with the energy of the transition suggesting that the infall is accelerating toward the center of the core, consistent with gravitational collapse. Despite the monolithic appearance of the Main core, the presence of red-shifted absorption, the existence of two embedded free-free sources at the center, and the rotational spin-up are consistent with an unstable core undergoing fragmentation with infall and differential rotation due to conservation of angular momentum. Therefore, the most likely explanation for the monolithic morphology is that the large opacity of the dust emission prevents the detection of any inhomogeneity in the core