Search

Deuterium fractionation is a well-established evolutionary tracer in low-mass star formation, but its applicability to the high-mass regime remains an open question. The abundances and ratios of deuterated species have often been proposed as reliable evolutionary indicators for different stages of the high-mass star formation. We investigate the role of N$_2$H$^+$ and key deuterated molecules as tracers of the different stages of the high-mass star formation, and test whether their abundance ratios can serve as reliable evolutionary indicators. We conducted APEX observations of o-H$_2$D$^+$ (1$_{10}$-1$_{11}$), N$_2$H$^+$ (4-3), and N$_2$d$^+$ (3-2) in 40 high-mass clumps at different evolutionary stages, selected from the ATLASGAL survey. Molecular column densities ($N$) and abundances ($X$), were derived through spectral line modelling, both under local thermodynamic equilibrium (LTE) and non-LTE conditions. The $N$(o-H$_2$D$^+$) show the smallest deviation from LTE results when derived under non-LTE assumptions. In contrast, N$_2$D$^+$ shows the largest discrepancy between the $N$ derived from LTE and non-LTE. In all the cases discussed, we found that $X$(o-H$_2$D$^+$) decreases more significantly with time than in the case of $X$(N$_2$D$^+$); whereas $X$(N$_2$H$^+$) increases slightly. Therefore, the validity of the recently proposed $X$(o-H$_2$D$^+$)/$X$(N$_2$D$^+$) ratio as a reliable evolutionary indicator was not observed for this sample. While the deuteration fraction derived from N$_2$D$^+$ and N$_2$H$^+$ clearly decreases with clump evolution, the interpretation of this trend is complex, given the different distribution of the two tracers. Our results suggest that a careful consideration of the observational biases and beam-dilution effects are crucial for an accurate interpretation of the evolution of the deuteration process during the high-mass star formation process., Comment: 13 pages, 5 figures, 3 tables. Accepted for publication in A&A

ALMA-IMF observed 15 massive protoclusters capturing multiple spectral lines and the continuum emission. We focus on the G351.77 protocluster ($\sim$ 2500 M$_{\odot}$, estimated from single-dish continuum observations) located at 2 kpc. We trace the dense gas emission and kinematics with N$_2$H$^+$ (1-0) at $\sim$ 4 kau resolution. We estimate an N$_2$H$^+$ relative abundance $\sim (1.7 \pm 0.5) \times 10^{-10}$. We decompose the N$_2$H$^+$ emission into up to two velocity components, highlighting the kinematic complexity in the dense gas. By examining the position-velocity (PV) diagrams on small scales, we observe clear inflow signatures (V-shapes) associated with dense cores. The most prominent V-shape has a mass inflow rate of $\sim 9.8 \times 10^{-4}$ M$_{\odot}$ yr$^{-1}$ and a short timescale of $\sim$ 15.6 kyr. We also observe V-shapes without associated cores. This suggests both that cores or centers of accretion exist below the 1.3 mm detection limit, and that the V-shapes may be viable tracers of very early accretion and star formation on $\sim$ 4 kau scales. The large-scale PV diagram shows that the protocluster is separated into 2 principal velocity structures. Combined with smaller scale DCN and H$_2$CO emission, we propose a scenario of larger scale slow contraction with rotation in the center based on simple toy models. This scenario leads the suggestion of outside-in evolution of the protocluster as it collapses. The gas depletion times implied by the V-shapes are short ($\sim$ 0.3 Myr), requiring either very fast cluster formation, and/or continuous mass feeding of the protocluster. The latter is possible via the Mother Filament G351.77 is forming out of. The similarities in the properties of G351.77 and the recently published work in G353.41 indicate that many of the physical conditions inferred via the ALMA-IMF N$_2$H$^+$ observations may be generic to protoclusters., Comment: Submitted in A&A, 28 pages, 31 figures, 4 interactive figures, 7 tables

The increasing number of observations towards different environments in the Milky Way, as well as theoretical and experimental works, are improving our knowledge of the astrochemical processes in the interstellar medium (ISM). In this chapter we report some of the main projects to study the chemical complexity and isotopic ratios across the Galaxy. High-sensitivity spectral surveys covering broad bandwidths towards Galactic Center molecular clouds (e.g. G+0.693-0.027) and star-forming regions (e.g. the hot core G31.41+0.31) are revealing very rich astrochemical reservoirs, which include molecules of prebiotic interest. At the same time, isotopic ratios (e.g. $^{12}$C/$^{13}$C and $^{14}$N/$^{15}$N) can give important information on the Galactic chemical evolution, as well as on chemical local processes due to the physical conditions of the molecular clouds. We also highlight the role of cosmic rays as a key agent affecting the interstellar chemistry described above., Comment: 9 pages, 1 figures - Chapter II of the National Congress of (proto-)Planetary Astrochemistry (CNAP II) Proceedings - Accepted for publication on MemSAIt, in print

The first generation of ELT instruments includes an optical-infrared high-resolution spectrograph, indicated as ELT-HIRES and recently christened ANDES (ArmazoNes high Dispersion Echelle Spectrograph). ANDES consists of three fibre-fed spectrographs ([U]BV, RIZ, YJH) providing a spectral resolution of $\sim$100,000 with a minimum simultaneous wavelength coverage of 0.4-1.8 $\mu$m with the goal of extending it to 0.35-2.4 $\mu$m with the addition of a U arm to the BV spectrograph and a separate K band spectrograph. It operates both in seeing- and diffraction-limited conditions and the fibre feeding allows several, interchangeable observing modes including a single conjugated adaptive optics module and a small diffraction-limited integral field unit in the NIR. Modularity and fibre-feeding allow ANDES to be placed partly on the ELT Nasmyth platform and partly in the Coud\'e room. ANDES has a wide range of groundbreaking science cases spanning nearly all areas of research in astrophysics and even fundamental physics. Among the top science cases, there are the detection of biosignatures from exoplanet atmospheres, finding the fingerprints of the first generation of stars, tests on the stability of Nature's fundamental couplings, and the direct detection of the cosmic acceleration. The ANDES project is carried forward by a large international consortium, composed of 35 Institutes from 13 countries, forming a team of almost 300 scientists and engineers which include the majority of the scientific and technical expertise in the field that can be found in ESO member states., Comment: SPIE astronomical telescope and instrumentation 2024, in press

We aim to determine the $^{14}$N/$^{15}$N and $^{12}$C/$^{13}$C ratios for HCN in six starless and prestellar cores and compare the results between the direct method using radiative transfer modeling and the indirect double isotope method assuming a fixed $^{12}$C/$^{13}$C ratio. We present IRAM 30m observations of the HCN 1-0, HCN 3-2, HC15N 1-0 and H13CN 1-0 transitions toward six embedded cores. The ${}^{14}$N/${}^{15}$N ratio was derived using both the indirect double isotope method and directly through non-local thermodynamic equilibrium (NLTE) 1D radiative transfer modeling of the HCN emission. The latter also provides the ${}^{12}$C/${}^{13}$C ratio, which we compared to the local interstellar value. The derived ${}^{14}$N/${}^{15}$N ratios using the indirect method are generally in the range of 300-550. This result could suggest an evolutionary trend in the nitrogen fractionation of HCN between starless cores and later stages of the star formation process. However, the direct method reveals lower fractionation ratios of around $\sim$250, mainly resulting from a lower ${}^{12}$C/${}^{13}$C ratio in the range $\sim$20-40, as compared to the local interstellar medium value of 68. This study reveals a significant difference between the nitrogen fractionation ratio in HCN derived using direct and indirect methods. This can influence the interpretation of the chemical evolution and reveal the pitfalls of the indirect double isotope method for fractionation studies. However, the direct method is challenging, as it requires well-constrained source models to produce accurate results. No trend in the nitrogen fractionation of HCN between earlier and later stages of the star formation process is evident when the results of the direct method are considered., Comment: Accepted for publication in A&A

Cosmic rays (CRs) heavily impact the chemistry and physics of cold and dense star-forming regions. However, characterising their ionisation rate is still challenging from an observational point of view. In the past, a few analytical formulas have been proposed to infer the cosmic-ray ionization rate $\zeta_2$ from molecular line observations. These have been derived from the chemical kinetics of the involved species, but they have not been validated using synthetic data processed with a standard observative pipeline. We aim to bridge this gap. We perform the radiative transfer on a set of three-dimensional magneto-hydrodynamical simulations of prestellar cores, exploring different initial $\zeta_2$, evolutionary stages, types of radiative transfer (e.g. assuming local-thermodynamic-equilibrium conditions), and telescope responses. We then compute the column densities of the involved tracers to determine $\zeta_2$, using, in particular, the equation proposed by Bovino et. al (2020) and by Caselli et al. (1998) both used nowadays. Our results confirm that the method of Bovino et al. (2020) accurately retrieves the actual $\zeta_2$ within a factor of $2-3$, in the physical conditions explored in our tests. Since we also explore a non-local thermodynamic equilibrium radiative transfer, this work indirectly offers insights into the excitation temperatures of common transitions at moderate volume densities ($n\approx 10^5 \, \rm cm^{-3}$). We have also performed a few tests using the formula proposed by Caselli et al. (1998), which overestimates the actual $\zeta_2$ by at least two orders of magnitudes. We also consider a new derivation of this method, which, however, still leads to large overestimates., Comment: Accepted for publication in A&A (15th Feb.)

Aims. We aim to improve a previous model for the prediction of column densities and deuterium fractions of non- and singly deuterated methanol. Thereby, we try to identify crucial chemical and physical parameters, for which the study of deuteration could provide valuable additional constraints. Methods. We employed a gas-grain chemical code to devise a model that is in agreement with the observed column density and deuterium fraction profiles of the innermost region of the pre-stellar core L1544. For that purpose, we developed a new treatment of reactive desorption, deriving an individual reactive desorption efficiency for every product species in a chemical reaction, that depends on the reaction enthalpy and type of underlying surface. Furthermore, we explored several options to promote the diffusion of hydrogen and deuterium atoms over the surface of interstellar dust grains, in order to increase methanol formation. Results. Our fiducial model employs diffusion by quantum tunneling of hydrogen and deuterium atoms, resulting in CH$_3$OH and CH$_2$DOH column densities that are approximately an order of magnitude lower than the observed values, which improves the results compared to the previous model by a factor 10. The $N$(CH$_2$DOH)/$N$(CH$_3$OH) ratio is reproduced within a factor of 1.2 for the centre and 1.8 for the position of the methanol peak. Given the large uncertainties that chemical models typically have, we consider our predictions to be in agreement with the observations. In general, we conclude that a diffusion process with a high diffusion rate needs to be employed to obtain methanol column densities that are in accordance with the observed values. Also, we find that the introduction of abstraction reactions into the methanol formation scheme suppresses deuteration, when used in combination with a high diffusion rate.

Studies of dense core morphologies and their orientations with respect to gas flows and the local magnetic field have been limited to only a small sample of cores with spectroscopic data. Leveraging the Green Bank Ammonia Survey alongside existing sub-millimeter continuum observations and Planck dust polarization, we produce a cross-matched catalogue of 399 dense cores with estimates of core morphology, size, mass, specific angular momentum, and magnetic field orientation. Of the 399 cores, 329 exhibit 2D $\mathrm{v}_\mathrm{LSR}$ maps that are well fit with a linear gradient, consistent with rotation projected on the sky. We find a best-fit specific angular momentum and core size relationship of $J/M \propto R^{1.82 \pm 0.10}$, suggesting that core velocity gradients originate from a combination of solid body rotation and turbulent motions. Most cores have no preferred orientation between the axis of core elongation, velocity gradient direction, and the ambient magnetic field orientation, favouring a triaxial and weakly magnetized origin. We find, however, strong evidence for a preferred anti-alignment between the core elongation axis and magnetic field for protostellar cores, revealing a change in orientation from starless and prestellar populations that may result from gravitational contraction in a magnetically-regulated (but not dominant) environment. We also find marginal evidence for anti-alignment between the core velocity gradient and magnetic field orientation in the L1228 and L1251 regions of Cepheus, suggesting a preferred orientation with respect to magnetic fields may be more prevalent in regions with locally ordered fields., Comment: 33 pages, 28 figures, accepted to MNRAS

In the centre of pre-stellar cores, deuterium fractionation is enhanced due to the low temperatures and high densities. Therefore, the chemistry of deuterated molecules can be used to study the earliest stages of star formation. We analyse the deuterium fractionation of simple molecules, comparing the level of deuteration in the envelopes of the pre-stellar core L1544 in Taurus and the protostellar core HH211 in Perseus. We used single-dish observations of CCH, HCN, HNC, HCO$^+$, and their $^{13}$C-, $^{18}$O- and D-bearing isotopologues, detected with the Onsala 20m telescope. We derived the column densities and the deuterium fractions of the molecules. Additionally, we used radiative transfer simulations and results from chemical modelling to reproduce the observed molecular lines. We used new collisional rate coefficients for HNC, HN$^{13}$C, DNC, and DCN that consider the hyperfine structure of these molecules. We find high levels of deuteration for CCH (10%) in both sources, consistent with other carbon chains, and moderate levels for HCN (5-7%) and HNC (8%). The deuterium fraction of HCO$^+$ is enhanced towards HH211, most likely caused by isotope-selective photodissociation of C$^{18}$O. Similar levels of deuteration show that the process is likely equally efficient towards both cores, suggesting that the protostellar envelope still retains the chemical composition of the original pre-stellar core. The fact that the two cores are embedded in different molecular clouds also suggests that environmental conditions do not have a significant effect on the deuteration within dense cores. Radiative transfer modelling shows that it is necessary to include the outer layers of the cores to consider the effects of extended structures. Besides HCO$^+$ observations, HCN observations towards L1544 also require the presence of an outer diffuse layer where the molecules are relatively abundant., Comment: 27 pages, 17 figures, accepted for publication in A&A

In spite of the advent of extremely large telescopes in the UV/optical/NIR range, the current generation of 8-10m facilities is likely to remain competitive at ground-UV wavelengths for the foreseeable future. The Cassegrain U-Band Efficient Spectrograph (CUBES) has been designed to provide high-efficiency (>40%) observations in the near UV (305-400 nm requirement, 300-420 nm goal) at a spectral resolving power of R>20,000, although a lower-resolution, sky-limited mode of R ~ 7,000 is also planned. CUBES will offer new possibilities in many fields of astrophysics, providing access to key lines of stellar spectra: a tremendous diversity of iron-peak and heavy elements, lighter elements (in particular Beryllium) and light-element molecules (CO, CN, OH), as well as Balmer lines and the Balmer jump (particularly important for young stellar objects). The UV range is also critical in extragalactic studies: the circumgalactic medium of distant galaxies, the contribution of different types of sources to the cosmic UV background, the measurement of H2 and primordial Deuterium in a regime of relatively transparent intergalactic medium, and follow-up of explosive transients. The CUBES project completed a Phase A conceptual design in June 2021 and has now entered the Phase B dedicated to detailed design and construction. First science operations are planned for 2028. In this paper, we briefly describe the CUBES project development and goals, the main science cases, the instrument design and the project organization and management., Comment: Proceedings for the HACK100 conference, Trieste, June 2022. arXiv admin note: substantial text overlap with arXiv:2208.01672

Prestellar cores, the birthplace of Sun-like stars, form from the fragmentation of the filamentary structure that composes molecular clouds, from which they must inherit at least partially the kinematics. Furthermore, when they are on the verge of gravitational collapse, they show signs of subsonic infall motions. How extended these motions are, which depends on how the collapse occurs, remains largely unknown. We want to investigate the kinematics of the envelope that surrounds the prototypical prestellar core L1544, studying the cloud-core connection. To our aims, we observed the $\rm HCO^+$(1-0) transition in a large map. \hcop is expected to be abundant in the envelope, making it an ideal probe of the large-scale kinematics in the source. We modelled the spectrum at the dust peak by means of a non local-thermodynamical-equilibrium radiative transfer. In order to reproduce the spectrum at the dust peak, a large ($\sim 1\, \rm pc$) envelope is needed, with low density (tens of $\rm cm^{-3}$ at most) and contraction motions, with an inward velocity of $\approx 0.05\,\rm km \, s^{-1}$. We fitted the data cube using the Hill5 model, which implements a simple model {for the optical depth and excitation temperature profiles along the line-of-sight,} in order to obtain a map of the infall velocity. This shows that the infall motions are extended, with typical values in the range $0.1-0.2\,\rm km \, s^{-1}$. Our results suggest that the contraction motions extend in the diffuse envelope surrounding the core, which is consistent with recent magnetic field measurements in the source, which showed that the envelope is magnetically supercritical., Comment: Accepted for publication on ApJ, 24 Oct. 2022

We aim to quantify the effect of chemistry on the infall velocity in the prestellar core L1544. Previous observational studies have found evidence for double-peaked line profiles for the rotational transitions of several molecules, which cannot be accounted for with the models presently available for the physical structure of the source, without ad hoc up-scaling of the infall velocity. We ran one-dimensional hydrodynamical simulations of the collapse of a core with L1544-like properties (in terms of mass and outer radius), using a state-of-the-art chemical model with a very large chemical network combined with an extensive description of molecular line cooling, determined via radiative transfer simulations, with the aim of determining whether these expansions of the simulation setup (as compared to previous models) can lead to a higher infall velocity. After running a series of simulations where the simulation was sequentially simplified, we found that the infall velocity is almost independent of the size of the chemical network or the approach to line cooling. We conclude that chemical evolution does not have a large impact on the infall velocity, and that the higher infall velocities that are implied by observations may be the result of the core being more dynamically evolved than what is now thought, or alternatively the average density in the simulated core is too low. However, chemistry does have a large influence on the lifetime of the core, which varies by about a factor of two across the simulations and grows longer when the chemical network is simplified. Therefore, although the model is subject to several sources of uncertainties, the present results clearly indicate that the use of a small chemical network leads to an incorrect estimate of the core lifetime, which is naturally a critical parameter for the development of chemical complexity in the precollapse phase., Comment: 12 pages incl. Appendix; 8 figures; accepted to A&A; replacement to reflect the change of the title of the paper post-acceptance

Spectral line profiles of several molecules observed towards the pre-stellar core L1544 appear double-peaked. For abundant molecular species this line morphology has been linked to self-absorption. However, the physical process behind the double-peaked morphology for less abundant species is still under debate. In order to understand the cause behind the double-peaked spectra of optically thin transitions and their link to the physical structure of pre-stellar cores, we present high-sensitivity and high-spectral resolution HC$^{17}$O$^+$ $J =$1-0 observations towards the dust peak in L1544. We observed the HC$^{17}$O$^+$ (1-0) spectrum with the Institut de Radioastronomie Millim\'etrique (IRAM) 30m telescope. By using new state-of-the-art collisional rate coefficients, a physical model for the core and the fractional abundance profile of HC$^{17}$O$^+$, the hyperfine structure of this molecular ion is modelled for the first time with the radiative transfer code LOC applied to the predicted chemical structure of a contracting pre-stellar core. We applied the same analysis to the chemically related C$^{17}$O molecule. The observed HC$^{17}$O$^+$(1-0) and C$^{17}$O(1-0) lines have been successfully reproduced with a non-local thermal equilibrium (LTE) radiative transfer model applied to chemical model predictions for a contracting pre-stellar core. An upscaled velocity profile (by 30%) is needed to reproduce the HC$^{17}$O$^+$(1-0) observations. The double peaks observed in the HC$^{17}$O$^+$(1-0) hyperfine components are due to the contraction motions at densities close to the critical density of the transition ($\sim$10$^{5}$ cm$^{-3}$) and to the fact that the HCO$^{+}$ fractional abundance decreases toward the centre.

High-mass star formation theories make distinct predictions on the properties of the prestellar seeds of high-mass stars. Observations of the early stages of high-mass star formation can provide crucial constraints, but they are challenging and scarce. We investigate the properties of the prestellar core population embedded in the high-mass clump AGAL014.492-00.139, and we study the kinematics at the clump and the clump-to-core scales. We have analysed an extensive dataset acquired with the ALMA interferometer. Applying a dendrogram analysis to the Band o-$\rm H_2D^+$ data, we identified 22 cores. We have fitted their average spectra in local-thermodinamic-equilibrium conditions, and we analysed their continuum emission at $0.8 \, \rm mm$. The cores have transonic to mildly supersonic turbulence levels and appear mostly low-mass, with $M_\mathrm{core}< 30 \, \rm M_\odot$. Furthermore, we have analysed Band 3 observations of the $\rm N_2H^+$ (1-0) transition, which traces the large scale gas kinematics. Using a friend-of-friend algorithm, we identify four main velocity coherent structures, all of which are associated with prestellar and protostellar cores. One of them presents a filament-like structure, and our observations could be consistent with mass accretion towards one of the protostars. In this case, we estimate a mass accretion rate of $ \dot{M}_\mathrm{acc}\approx 2 \times 10^{-4} \rm \, M_\odot \, yr^{-1}$. Our results support a clump-fed accretion scenario in the targeted source. The cores in prestellar stage are essentially low-mass, and they appear subvirial and gravitationally bound, unless further support is available for instance due to magnetic fields., Comment: Accepted for publication in ApJ

In the era of Extremely Large Telescopes, the current generation of 8-10m facilities are likely to remain competitive at ground-UV wavelengths for the foreseeable future. The Cassegrain U-Band Efficient Spectrograph (CUBES) has been designed to provide high-efficiency (>40%) observations in the near UV (305-400 nm requirement, 300-420 nm goal) at a spectral resolving power of R>20,000 (with a lower-resolution, sky-limited mode of R ~ 7,000). With the design focusing on maximizing the instrument throughput (ensuring a Signal to Noise Ratio (SNR) ~20 per high-resolution element at 313 nm for U ~18.5 mag objects in 1h of observations), it will offer new possibilities in many fields of astrophysics, providing access to key lines of stellar spectra: a tremendous diversity of iron-peak and heavy elements, lighter elements (in particular Beryllium) and light-element molecules (CO, CN, OH), as well as Balmer lines and the Balmer jump (particularly important for young stellar objects). The UV range is also critical in extragalactic studies: the circumgalactic medium of distant galaxies, the contribution of different types of sources to the cosmic UV background, the measurement of H2 and primordial Deuterium in a regime of relatively transparent intergalactic medium, and follow-up of explosive transients. The CUBES project completed a Phase A conceptual design in June 2021 and has now entered the detailed design and construction phase. First science operations are planned for 2028., Comment: SPIE proceedings, SPIE Astronomical Telescopes + Instrumentation 2022, Montr\'eal, Canada; 20 pages, 13 figures, 2 tables

The Beam Expander Testing X-ray facility (BEaTriX) is a unique X-ray apparatus now operated at the Istituto Nazionale di Astrofisica (INAF), Osservatorio Astronomico di Brera (OAB), in Merate, Italy. It has been specifically designed to measure the point spread function (PSF) and the effective area (EA) of the X-ray mirror modules (MMs) of the Advanced Telescope for High-ENergy Astrophysics (ATHENA), based on silicon pore optics (SPO) technology, for verification before integration into the mirror assembly. To this end, BEaTriX generates a broad, uniform, monochromatic, and collimated X-ray beam at 4.51 keV. [...] In BEaTriX, a micro-focus X-ray source with a titanium anode is placed in the focus of a paraboloidal mirror, which generates a parallel beam. A crystal monochromator selects the 4.51 keV line, which is expanded to the final size by a crystal asymmetrically cut with respect to the crystalline planes. [...] After characterization, the BEaTriX beam has the nominal dimensions of 170 mm x 60 mm, with a vertical divergence of 1.65 arcsec and a horizontal divergence varying between 2.7 and 3.45 arcsec, depending on the monochromator setting: either high collimation or high intensity. The flux per area unit varies from 10 to 50 photons/s/cm2 from one configuration to the other. The BEaTriX beam performance was tested using an SPO MM, whose entrance pupil was fully illuminated by the expanded beam, and its focus was directly imaged onto the camera. The first light test returned a PSF and an EA in full agreement with expectations. As of today, the 4.51 keV beamline of BEaTriX is operational and can characterize modular X-ray optics, measuring their PSF and EA with a typical exposure of 30 minutes. [...] We expect BEaTriX to be a crucial facility for the functional test of modular X-ray optics, such as the SPO MMs for ATHENA., Comment: Accepted on 29 June 2022 for publication in Astronomy & Astrophysics, Sect. 13, Astronomical instrumentation of Astronomy and Astrophysics. Version after language correction

Context: In the centre of pre-stellar cores, the deuterium fractionation is enhanced due to the cold temperatures and high densities. Therefore, the chemistry of deuterated molecules can be used to probe the evolution and the kinematics in the earliest stages of star formation. Aims: We analyse emission maps of cyclopropenylidene, c-C$_3$H$_2$, to study the distribution of the deuteration throughout the prototypical pre-stellar core L1544. Methods: We use single-dish observations of c-C$_3$H$_2$, c-H$^{13}$CC$_2$H, c-C$_3$HD, and c-C$_3$D$_2$ towards the pre-stellar core L1544, performed at the IRAM 30m telescope. We derive the column density and deuterium fraction maps, and compare these observations with non-LTE radiative transfer simulations. Results: The highest deuterium fractions are found close to the dust peak at the centre of L1544, where the increased abundance of H$_2$D$^+$ ions drives the deuteration process. The peak values are N(c-C$_3$HD)/N(c-C$_3$H$_2)=0.17\pm0.01$, N(c-C$_3$D$_2$)/N(c-C$_3$H$_2)=0.025\pm0.003$ and N(c-C$_3$D$_2$)/N(c-C$_3$HD$)=0.16\pm0.03$, which is consistent with previous single point observations. The distributions of c-C$_3$HD and c-C$_3$D$_2$ indicate that the deuterated forms of c-C$_3$H$_2$ in fact trace the dust peak and not the c-C$_3$H$_2$ peak. Conclusions: The N(c-C$_3$D$_2$)/N(c-C$_3$HD) map confirms that the process of deuteration is more efficient towards the centre of the core and demonstrates that carbon-chain molecules are still present at high densities. This is likely caused by an increased abundance of He$^+$ ions destroying CO, which increases the amount of carbon atoms in the gas phase., Comment: 13 pages, 9 figures

Most stars in the Galaxy, including the Sun, were born in high-mass star-forming regions. It is hence important to study the chemical processes in these regions to better understand the chemical heritage of both the Solar System and most stellar systems in the Galaxy. The molecular ion HCNH+ is thought to be a crucial species in ion-neutral astrochemical reactions, but so far it has been detected only in a handful of star-forming regions, and hence its chemistry is poorly known. We have observed with the IRAM-30m Telescope 26 high-mass star-forming cores in different evolutionary stages in the J=3-2 rotational transition of HCNH+. We report the detection of HCNH+ in 16 out of 26 targets. This represents the largest sample of sources detected in this molecular ion so far. The fractional abundances of HCNH+, [HCNH+], w.r.t. H2, are in the range 0.9 - 14 X $10^{-11}$, and the highest values are found towards cold starless cores. The abundance ratios [HCNH+]/[HCN] and [HCNH+]/[HCO+] are both < 0.01 for all objects except for four starless cores, for which they are well above this threshold. These sources have the lowest gas temperature in the sample. We run two chemical models, a "cold" one and a "warm" one, which attempt to match as much as possible the average physical properties of the cold(er) starless cores and of the warm(er) targets. The reactions occurring in the latter case are investigated in this work for the first time. Our predictions indicate that in the warm model HCNH+ is mainly produced by reactions with HCN and HCO+, while in the cold one the main progenitor species of HCNH+ are HCN+ and HNC+. The results indicate that the chemistry of HCNH+ is different in cold/early and warm/evolved cores, and the abundance ratios [HCNH+]/[HCN] and [HCNH+]/[HCO+] is a useful astrochemical tool to discriminate between different evolutionary phases in the process of star formation., Comment: 13 pages, 3 tables, 12 figures (+3 figures in appendix). Accepted for publication in Astronomy & Astrophysics

Context. The different theoretical models concerning the formation of high-mass stars make distinct predictions regarding their progenitors, i.e. the high-mass prestellar cores. However, so far no conclusive observation of such objects has been made. Aims. We aim to study the very early stages of high-mass star formation in two infrared-dark, massive clumps, to identify the core population that they harbour. Methods. We obtained ALMA observations of continuum emission at 0.8mm and of the ortho-$\rm H_2D^+$ transition at 372GHz towards the two clumps. We use the SCIMES algorithm to identify cores in the position-position-velocity space, finding 16 cores. We model their observed spectra in the LTE approximation, deriving the centroid velocity, linewidth, and column density maps. We also study the correlation between the continuum and molecular data, which in general do not present the same structure. Results. We report for the first time the detection of ortho-$\rm H_2D^+$ in high-mass star-forming regions performed with an interferometer. The molecular emission shows narrow and subsonic lines, suggesting that locally the temperature of the gas is less than 10K. From the continuum emission we estimate the cores' total masses, and compare them with the respective virial masses. We also compute the volume density values, which are found to be higher than $10^{6}\, \rm cm^{-3}$. Conclusions. Our data confirm that ortho-$\rm H_2D^+$ is an ideal tracer of cold and dense gas. Interestingly, almost all the $\rm H_2D^+$-identified cores are less massive than 13M_sun , with the exception of one core in AG354. Furthermore, most of them are subvirial and larger than their Jeans masses. These results are difficult to explain in the context of the turbulent accretion models, which predict massive and virialised prestellar cores., Comment: Accepted for publications in A&A

SOXS (Son Of X-Shooter) is a medium resolution (~4500) wide-band (0.35 - 2.0 {\mu}m) spectrograph which passed the Final Design Review in 2018. The instrument is planned to be installed at the NTT in La Silla and it is mainly composed by five different optomechanical subsystems (Common Path, NIR spectrograph, UV-VIS spectrograph, Camera, and Calibration) and other mechanical subsystems (Interface flange, Platform, cable corotator, and cooling). It is currently in the procurement and integration phase. In this paper we present the post-FDR modifications in the mechanical design due to the various iterations with the manufacturers and the actual procurement status. The last part describes the strategy used to keep under control the mechanical interfaces between the subsystems.

Aims. We report on ESPRESSO high-resolution transmission spectroscopic observations of two primary transits of the highly-irradiated, ultra-hot Jupiter-size planet WASP-76b. We investigate the presence of several key atomic and molecular features of interest that may reveal the atmospheric properties of the planet. Methods. We extracted two transmission spectra of WASP-76b with R approx 140,000 using a procedure that allowed us to process the full ESPRESSO wavelength range (3800-7880 A) simultaneously. We observed that at a high signal-to-noise ratio, the continuum of ESPRESSO spectra shows wiggles that are likely caused by an interference pattern outside the spectrograph. To search for the planetary features, we visually analysed the extracted transmission spectra and cross-correlated the observations against theoretical spectra of different atomic and molecular species. Results. The following atomic features are detected: Li I, Na I, Mg I, Ca II, Mn I, K I, and Fe I. All are detected with a confidence level between 9.2 sigma (Na I) and 2.8 sigma (Mg I). We did not detect the following species: Ti I, Cr I, Ni I, TiO, VO, and ZrO. We impose the following 1 sigma upper limits on their detectability: 60, 77, 122, 6, 8, and 8 ppm, respectively. Conclusions. We report the detection of Li I on WASP-76b for the first time. In addition, we found the presence of Na I and Fe I as previously reported in the literature. We show that the procedure employed in this work can detect features down to the level of ~ 0.1 % in the transmission spectrum and ~ 10 ppm by means of a cross-correlation method. We discuss the presence of neutral and singly ionised features in the atmosphere of WASP-76b., Comment: 20 pages, 19 figures, accepted for publication in Astronomy and Astrophysics

Context. The nitrogen isotopic ratio is considered an important diagnostic tool of the star formation process, and $N_2H^+$ is particularly important because it is directly linked to molecular nitrogen $N_2$. However, theoretical models still lack to provide an exhaustive explanation for the observed $^{14}N/^{15}N$ values. Aims. Recent theoretical works suggest that the $^{14}N/^{15}N$ behaviour is dominated by two competing reactions that destroy $ N_2H^+$: dissociative recombination and reaction with CO. When CO is depleted from the gas phase, if $N_2H^+$ recombination rate is lower with respect to the $N^{15}NH^+$ one, the rarer isotopologue is destroyed faster. This implies that the $N_2H^+$ isotopic ratio in protostars should be lower than the one in prestellar cores, and consistent with the elemental value of ~440. We aim to test this hypothesis, producing the first sample of $N_2H^+ / N^{15}NH^+$ measurements in low mass protostars. Methods. We observe the $N_2H^+$ and $N^{15}NH^+$ lowest rotational transition towards six young stellar objects in Perseus and Taurus molecular clouds. We model the spectra with a custom python code using a constant $T_{ex}$ approach to fit the observations. We discuss in appendix the validity of this hypothesis. The derived column densities are used to compute the nitrogen isotopic ratio. Results. Our analysis yields an average of $\rm ^{14}N/^{15}N|_{pro} = 420 \pm 15$ in the protostellar sample. This is consistent with the protosolar value of 440, and significantly lower than the average value previously obtained in a sample of prestellar objects. Conclusions. Our results are in agreement with the hypothesis that, when CO is depleted from the gas-phase, dissociative recombinations with free electrons destroy $N^{15}NH^+$ faster than $N_2H^+$, leading to high isotopic ratios in prestellar cores, where CO is frozen on dust grains., Comment: Accepted on A&A on 09 Oct 2020

ESPRESSO is the new high-resolution spectrograph of ESO's Very-Large Telescope (VLT). It was designed for ultra-high radial-velocity precision and extreme spectral fidelity with the aim of performing exoplanet research and fundamental astrophysical experiments with unprecedented precision and accuracy. It is able to observe with any of the four Unit Telescopes (UT) of the VLT at a spectral resolving power of 140,000 or 190,000 over the 378.2 to 788.7 nm wavelength range, or with all UTs together, turning the VLT into a 16-m diameter equivalent telescope in terms of collecting area, while still providing a resolving power of 70,000. We provide a general description of the ESPRESSO instrument, report on the actual on-sky performance, and present our Guaranteed-Time Observation (GTO) program with its first results. ESPRESSO was installed on the Paranal Observatory in fall 2017. Commissioning (on-sky testing) was conducted between December 2017 and September 2018. The instrument saw its official start of operations on October 1st, 2018, but improvements to the instrument and re-commissioning runs were conducted until July 2019. The measured overall optical throughput of ESPRESSO at 550 nm and a seeing of 0.65 arcsec exceeds the 10% mark under nominal astro-climatic conditions. We demonstrate a radial-velocity precision of better than 25 cm/s during one night and 50 cm/s over several months. These values being limited by photon noise and stellar jitter show that the performanceis compatible with an instrumental precision of 10 cm/s. No difference has been measured across the UTs neither in throughput nor RV precision. The combination of the large collecting telescope area with the efficiency and the exquisite spectral fidelity of ESPRESSO opens a new parameter space in RV measurements, the study of planetary atmospheres, fundamental constants, stellar characterisation and many other fields., Comment: 26 pages, 28 figures

Deuterium fractionation processes in the interstellar medium (ISM) have been shown to be highly efficient in the family of nitrogen hydrides. To date, observations were limited to ammonia (NH$_2$D, NHD$_2$, ND$_3$) and imidogen radical (ND) isotopologues. We want to explore the high frequency windows offered by the \emph{Herschel Space Observatory} to search for deuterated forms of amidogen radical NH$_2$ and to compare the observations against the predictions of our comprehensive gas-grain chemical model. Making use of the new molecular spectroscopy data recently obtained at high frequencies for NHD and ND$_2$, both isotopologues have been searched for in the spectral survey towards the class 0 IRAS 16293-2422, a source in which NH$_3$, NH and their deuterated variants have been previously detected. We used the observations carried out with HIFI (Heterodyne Instrument for the Far Infrared) in the framework of the key program "Chemical Herschel surveys of star forming regions" (CHESS). We report the first detection of interstellar NHD and ND$_2$. Both species are observed in absorption against the continuum of the protostar. From the analysis of their hyperfine structure, accurate excitation temperature and column density values have been determined. The latter were combined with the column density of the parent species NH$_2$ to derive the deuterium fractionation in amidogen. The amidogen D/H ratio measured in the low-mass protostar IRAS 16293-2422 is comparable to the one derived for the related species imidogen and much higher than that observed for ammonia. Additional observations of these species will give more insights into the mechanism of ammonia formation and deuteration in the ISM. We finally indicate the current possibilities to further explore these species at submillimeter wavelengths., Comment: 11 pages, 5 figures, 7 tables. Accepted for publication in A&A

The bright star $\pi$ Men was chosen as the first target for a radial velocity follow-up to test the performance of ESPRESSO, the new high-resolution spectrograph at the ESO's Very-Large Telescope (VLT). The star hosts a multi-planet system (a transiting 4 M$_\oplus$ planet at $\sim$0.07 au, and a sub-stellar companion on a $\sim$2100-day eccentric orbit) which is particularly appealing for a precise multi-technique characterization. With the new ESPRESSO observations, that cover a time span of 200 days, we aim to improve the precision and accuracy of the planet parameters and search for additional low-mass companions. We also take advantage of new photometric transits of $\pi$ Men c observed by TESS over a time span that overlaps with that of the ESPRESSO follow-up campaign. We analyse the enlarged spectroscopic and photometric datasets and compare the results to those in the literature. We further characterize the system by means of absolute astrometry with Hipparcos and Gaia. We used the spectra of ESPRESSO for an independent determination of the stellar fundamental parameters. We present a precise characterization of the planetary system around $\pi$ Men. The ESPRESSO radial velocities alone (with typical uncertainty of 10 cm/s) allow for a precise retrieval of the Doppler signal induced by $\pi$ Men c. The residuals show an RMS of 1.2 m/s, and we can exclude companions with a minimum mass less than $\sim$2 M$_\oplus$ within the orbit of $\pi$ Men c). We improve the ephemeris of $\pi$ Men c using 18 additional TESS transits, and in combination with the astrometric measurements, we determine the inclination of the orbital plane of $\pi$ Men b with high precision ($i_{b}=45.8^{+1.4}_{-1.1}$ deg). This leads to the precise measurement of its absolute mass $m_{b}=14.1^{+0.5}_{-0.4}$ M$_{Jup}$, and shows that the planetary orbital planes are highly misaligned., Comment: Accepted for publication on A&A

We characterized the transiting planetary system orbiting the G2V star K2-38 using the new-generation echelle spectrograph ESPRESSO. We carried out a photometric analysis of the available K2 photometric light curve of this star to measure the radius of its two known planets. Using 43 ESPRESSO high-precision radial velocity measurements taken over the course of 8 months along with the 14 previously published HIRES RV measurements, we modeled the orbits of the two planets through a MCMC analysis, significantly improving their mass measurements. Using ESPRESSO spectra, we derived the stellar parameters, $T_{\rm eff}$=5731$\pm$66, $\log g$=4.38$\pm$0.11~dex, and $[Fe/H]$=0.26$\pm$0.05~dex, and thus the mass and radius of K2-38, $M_{\star}$=1.03 $^{+0.04}_{-0.02}$~M$_{\oplus}$ and $R_{\star}$=1.06 $^{+0.09}_{-0.06}$~R$_{\oplus}$. We determine new values for the planetary properties of both planets. We characterize K2-38b as a super-Earth with $R_{\rm P}$=1.54$\pm$0.14~R$_{\rm \oplus}$ and $M_{\rm p}$=7.3$^{+1.1}_{-1.0}$~M$_{\oplus}$, and K2-38c as a sub-Neptune with $R_{\rm P}$=2.29$\pm$0.26~R$_{\rm \oplus}$ and $M_{\rm p}$=8.3$^{+1.3}_{-1.3}$~M$_{\oplus}$. We derived a mean density of $\rho_{\rm p}$=11.0$^{+4.1}_{-2.8}$~g cm$^{-3}$ for K2-38b and $\rho_{\rm p}$=3.8$^{+1.8}_{-1.1}$~g~cm$^{-3}$ for K2-38c, confirming K2-38b as one of the densest planets known to date. The best description for the composition of K2-38b comes from an iron-rich Mercury-like model, while K2-38c is better described by a rocky model with a H2 envelope. The maximum collision stripping boundary shows how giant impacts could be the cause for the high density of K2-38b. The irradiation received by each planet places them on opposite sides of the radius valley. We find evidence of a long-period signal in the radial velocity time-series whose origin could be linked to a 0.25-3~M$_{\rm J}$ planet or stellar activity., Comment: 16 pages, 14 figures

We aim to confirm the presence of Proxima b using independent measurements obtained with the new ESPRESSO spectrograph, and refine the planetary parameters taking advantage of its improved precision. We analysed 63 spectroscopic ESPRESSO observations of Proxima taken during 2019. We obtained radial velocity measurements with a typical radial velocity photon noise of 26 cm/s. We ran a joint MCMC analysis on the time series of the radial velocity and full-width half maximum of the cross-correlation function to model the planetary and stellar signals present in the data, applying Gaussian process regression to deal with stellar activity. We confirm the presence of Proxima b independently in the ESPRESSO data. The ESPRESSO data on its own shows Proxima b at a period of 11.218 $\pm$ 0.029 days, with a minimum mass of 1.29 $\pm$ 0.13 Me. In the combined dataset we measure a period of 11.18427 $\pm$ 0.00070 days with a minimum mass of 1.173 $\pm$ 0.086 Me. We find no evidence of stellar activity as a potential cause for the 11.2 days signal. We find some evidence for the presence of a second short-period signal, at 5.15 days with a semi-amplitude of merely 40 cm/s. If caused by a planetary companion, it would correspond to a minimum mass of 0.29 $\pm$ 0.08 Me. We find that the FWHM of the CCF can be used as a proxy for the brightness changes and that its gradient with time can be used to successfully detrend the radial velocity data from part of the influence of stellar activity. The activity-induced radial velocity signal in the ESPRESSO data shows a trend in amplitude towards redder wavelengths. Velocities measured using the red end of the spectrograph are less affected by activity, suggesting that the stellar activity is spot-dominated. The data collected excludes the presence of extra companions with masses above 0.6 Me at periods shorter than 50 days., Comment: 25 pages, 26 figures

We carried out a parameter-space exploration of the ammonia abundance in the pre-stellar core L1544, where it has been observed to increase toward the center of the core with no signs of freeze-out onto grain surfaces. We considered static and dynamical physical models coupled with elaborate chemical and radiative transfer calculations, and explored the effects of varying model parameters on the (ortho+para) ammonia abundance profile. None of our models are able to reproduce the inward-increasing tendency in the observed profile; ammonia depletion always occurs in the center of the core. In particular, our study shows that including the chemical desorption process, where exothermic association reactions on the grain surface can result in the immediate desorption of the product molecule, leads to ammonia abundances that are over an order of magnitude above the observed level in the innermost 15000 au of the core - at least when one employs a constant efficiency for the chemical desorption process irrespective of the ice composition. Our results seemingly constrain the chemical desorption efficiency of ammonia on water ice to below 1%. It is increasingly evident that time-dependent effects must be considered so that the results of chemical models can be reconciled with observations., Comment: Accepted to MNRAS

L1521E is a dense starless core in Taurus that was found to have relatively low molecular depletion by earlier studies, thus suggesting a recent formation. We aim to characterize the chemical structure of L1521E and compare it to the more evolved L1544 pre-stellar core. We have obtained $\sim$2.5$\times$2.5 arcminute maps toward L1521E using the IRAM-30m telescope in transitions of various species. We derived abundances for the species and compared them to those obtained toward L1544. We estimated CO depletion factors. Similarly to L1544, $c$-C$_3$H$_2$ and CH$_3$OH peak at different positions. Most species peak toward the $c$-C$_3$H$_2$ peak. The CO depletion factor derived toward the $Herschel$ dust peak is 4.3$\pm$1.6, which is about a factor of three lower than that toward L1544. The abundances of sulfur-bearing molecules are higher toward L1521E than toward L1544 by factors of $\sim$2-20. The abundance of methanol is similar toward the two cores. The higher abundances of sulfur-bearing species toward L1521E than toward L1544 suggest that significant sulfur depletion takes place during the dynamical evolution of dense cores, from the starless to pre-stellar stage. The CO depletion factor measured toward L1521E suggests that CO is more depleted than previously found. Similar CH$_3$OH abundances between L1521E and L1544 hint that methanol is forming at specific physical conditions in Taurus, characterized by densities of a few $\times$10$^4$ cm$^{-3}$ and $N$(H$_2$)$\gtrsim$10$^{22}$ cm$^{-2}$, when CO starts to catastrophically freeze-out, while water can still be significantly photodissociated, so that the surfaces of dust grains become rich in solid CO and CH$_3$OH, as already found toward L1544. Methanol can thus provide selective crucial information about the transition region between dense cores and the surrounding parent cloud., Comment: Accepted for publication in A&A, abstract abridged

Understanding the early stages of star formation is a research field of ongoing development, both theoretically and observationally. In this context, molecular data have been continuously providing observational constraints on the gas dynamics at different excitation conditions and depths in the sources. We have investigated the Barnard 59 core, the only active site of star formation in the Pipe Nebula, to achieve a comprehensive view of the kinematic properties of the source. These information were derived by simultaneously fitting ammonia inversion transition lines (1,1) and (2,2). Our analysis unveils the imprint of protostellar feedback, such as increasing line widths, temperature and turbulent motions in our molecular data. Combined with complementary observations of dust thermal emission, we estimate that the core is gravitationally bound following a virial analysis. If the core is not contracting, another source of internal pressure, most likely the magnetic field, is supporting it against gravitational collapse and limits its star formation efficiency., Comment: 18 pages, 18 figures