Showing 1–50 of 526 results for author: Caselli, P

Probing the Physics of Star-Formation (ProPStar) III. No evidence for dissipation of turbulence down to 20 mpc (4 000 au) scale

Authors: Jaime E. Pineda, Juan D. Soler, Stella Offner, Eric W. Koch, Dominique M. Segura-Cox, Roberto Neri, Michael Kuffmeier, Alexei V. Ivlev, Maria Teresa Valdivia-Mena, Olli Sipilä, Maria Jose Maureira, Paola Caselli, Nichol Cunningham, Anika Schmiedeke, Caroline Gieser, Michael Chen, Silvia Spezzano

Abstract: Context. Turbulence is a key component of molecular cloud structure. It is usually described by a cascade of energy down to the dissipation scale. The power spectrum for subsonic incompressible turbulence is $k^{-5/3}$, while for supersonic turbulence it is $k^{-2}$. Aims. We aim to determine the power spectrum in an actively star-forming molecular cloud, from parsec scales down to the expected ma… ▽ More Context. Turbulence is a key component of molecular cloud structure. It is usually described by a cascade of energy down to the dissipation scale. The power spectrum for subsonic incompressible turbulence is $k^{-5/3}$, while for supersonic turbulence it is $k^{-2}$. Aims. We aim to determine the power spectrum in an actively star-forming molecular cloud, from parsec scales down to the expected magnetohydrodynamic (MHD) wave cutoff (dissipation scale). Methods. We analyze observations of the nearby NGC 1333 star-forming region in three different tracers to cover the different scales from $\sim$10 pc down to 20 mpc. The largest scales are covered with the low density gas tracer $^{13}$CO (1-0) obtained with single dish, the intermediate scales are covered with single-dish observations of the C$^{18}$O (3-2) line, while the smallest scales are covered in H$^{13}$CO$^+$ (1-0) and HNC (1-0) with a combination of NOEMA interferometer and IRAM 30m single dish observations. The complementarity of these observations enables us to generate a combined power spectrum covering more than two orders of magnitude in spatial scale. Results. We derive the power spectrum in an active star-forming region spanning more than 2 decades of spatial scales. The power spectrum of the intensity maps shows a single power-law behavior, with an exponent of 2.9$\pm$0.1 and no evidence of dissipation. Moreover, there is evidence for the power-spectrum of the ions to have more power at smaller scales than the neutrals, which is opposite from theoretical expectations. Conclusions. We show new possibilities of studying the dissipation of energy at small scales in star-forming regions provided by interferometric observations. △ Less

Submitted 31 August, 2024; originally announced September 2024.

Comments: 6 pages (+2 appendix), 5 figures. Accepted to A&A Letters. GitHub repo with code available at https://github.com/jpinedaf/NGC1333_NOEMA_turbulence

Impact of ice growth on the physical and chemical properties of dense cloud cores

Authors: O. Sipilä, P. Caselli, M. Juvela

Abstract: We investigated the effect of time-dependent ice growth on dust grains on the opacity and hence on the dust temperature in a collapsing molecular cloud core, with the aim of quantifying the effect of the dust temperature variations on ice abundances as well as the evolution of the collapse. We employed a one-dimensional collapse model that self-consistently and time-dependently combines hydrodynam… ▽ More We investigated the effect of time-dependent ice growth on dust grains on the opacity and hence on the dust temperature in a collapsing molecular cloud core, with the aim of quantifying the effect of the dust temperature variations on ice abundances as well as the evolution of the collapse. We employed a one-dimensional collapse model that self-consistently and time-dependently combines hydrodynamics with chemical and radiative transfer simulations. The dust opacity was updated on-the-fly based on the ice growth as a function of location in the core. The results of the fully dynamical model were compared against simulations assuming fixed ice thickness. We found that the ice thickness increases fast and reaches a saturation value of approximately 90 monolayers in the central core (volume density $\sim$$10^4\,\rm cm^{-3}$), and several tens of monolayers at a volume density of $\sim$$10^3\,\rm cm^{-3}$, after only a few $10^5\,\rm yr$ of evolution. The results thus exclude the adoption of thin ($\sim$10 monolayer) ices in molecular cloud simulations except at very short timescales. The differences in abundances and dust temperature between the fully dynamic simulation and those with fixed dust opacity are small; abundances change between the solutions generally within a factor of two. The assumptions on the dust opacity do have an effect on the collapse dynamics through the influence of the photoelectric effect on the gas temperature, and the simulations take a different time to reach a common central density. In conclusion, carrying out chemical simulations using a dust temperature corresponding to a fixed opacity seems to be a good approximation. Still, although at least in the present case its effect on the overall results is limited - as long as the grains are monodisperse - ice growth should be considered to obtain the most accurate representation of the collapse dynamics. △ Less

Submitted 30 August, 2024; originally announced August 2024.

Comments: Accepted for publication in A&A; abstract abridged to meet arXiv requirements

FAUST XIX. D$_2$CO in the outflow cavities of NGC\,1333 IRAS\,4A: recovering the physical structure of its original prestellar core

Authors: Layal Chahine, Cecilia Ceccarelli, Marta De Simone, Claire J. Chandler, Claudio Codella, Linda Podio, Ana López-Sepulcre, Brian Svoboda, Giovanni Sabatini, Nami Sakai, Laurent Loinard, Charlotte Vastel, Nadia Balucani, Albert Rimola, Piero Ugliengo, Yuri Aikawa, Eleonora Bianchi, Mathilde Bouvier, Paola Caselli, Steven Charnley, Nicolás Cuello, Tomoyuki Hanawa, Doug Johnstone, Maria José Maureira, Francois Ménard , et al. (3 additional authors not shown)

Abstract: Molecular deuteration is a powerful diagnostic tool for probing the physical conditions and chemical processes in astrophysical environments. In this work, we focus on formaldehyde deuteration in the protobinary system NGC\,1333 IRAS\,4A, located in the Perseus molecular cloud. Using high-resolution ($\sim$\,100\,au) ALMA observations, we investigate the [D$_2$CO]/[HDCO] ratio along the cavity wal… ▽ More Molecular deuteration is a powerful diagnostic tool for probing the physical conditions and chemical processes in astrophysical environments. In this work, we focus on formaldehyde deuteration in the protobinary system NGC\,1333 IRAS\,4A, located in the Perseus molecular cloud. Using high-resolution ($\sim$\,100\,au) ALMA observations, we investigate the [D$_2$CO]/[HDCO] ratio along the cavity walls of the outflows emanating from IRAS\,4A1. Our analysis reveals a consistent decrease in the deuteration ratio (from $\sim$\,60-20\% to $\sim$\,10\%) with increasing distance from the protostar (from $\sim$\,2000\,au to $\sim$\,4000\,au). Given the large measured [D$_2$CO]/[HDCO], both HDCO and D$_2$CO are likely injected by the shocks along the cavity walls into the gas-phase from the dust mantles, formed in the previous prestellar phase. We propose that the observed [D$_2$CO]/[HDCO] decrease is due to the density profile of the prestellar core from which NGC\,1333 IRAS\,4A was born. When considering the chemical processes at the base of formaldehyde deuteration, the IRAS\,4A's prestellar precursor had a predominantly flat density profile within 3000\,au and a decrease of density beyond this radius. △ Less

Submitted 28 August, 2024; originally announced August 2024.

Assessing realistic binding energies of some essential interstellar radicals with amorphous solid water: A fully quantum chemical approach

Authors: Milan Sil, Arghyadeb Roy, Prasanta Gorai, Naoki Nakatani, Takashi Shimonishi, Kenji Furuya, Natalia Inostroza-Pino, Paola Caselli, Ankan Das

Abstract: In the absence of laboratory data, state-of-the-art quantum chemical approaches can provide estimates on the binding energy (BE) of interstellar species with grains. Without BE values, contemporary astrochemical models are compelled to utilize wild guesses, often delivering misleading information. Here, we employ a fully quantum chemical approach to estimate the BE of 7 diatomic radicals - CH, NH,… ▽ More In the absence of laboratory data, state-of-the-art quantum chemical approaches can provide estimates on the binding energy (BE) of interstellar species with grains. Without BE values, contemporary astrochemical models are compelled to utilize wild guesses, often delivering misleading information. Here, we employ a fully quantum chemical approach to estimate the BE of 7 diatomic radicals - CH, NH, OH, SH, CN, NS, and NO, playing a crucial role in shaping the interstellar chemical composition, using a suitable amorphous solid water model as a substrate since water is the principal constituent of interstellar ice in dense and shielded regions. While the BEs are compatible with physisorption, the binding of CH in some sites shows chemisorption, in which a chemical bond to an oxygen atom of a water molecule is formed. While no structural change can be observed for the CN radical, it is believed that the formation of a hemibonded system between the outer layer of the water cluster and the radical is the reason for the unusually large BE in one of the binding sites considered in our study. A significantly lower BE for NO, consistent with recent calculations, is obtained, which is useful in explaining the recently observed HONO/NH2OH and HONO/HNO ratio in the low-mass hot corino, IRAS 16293-2422B, with chemical models. △ Less

Submitted 26 August, 2024; originally announced August 2024.

Re-evaluation of the cosmic-ray ionization rate in diffuse clouds

Authors: M. Obolentseva, A. V. Ivlev, K. Silsbee, D. A. Neufeld, P. Caselli, G. Edenhofer, N. Indriolo, T. G. Bisbas, D. Lomeli

Abstract: All current estimates of the cosmic-ray (CR) ionization rate rely on assessments of the gas density along the probed sight lines. Until now, these have been based on observations of different tracers, with C$_2$ being the most widely used in diffuse molecular clouds for this purpose. However, three-dimensional dust extinction maps have recently reached sufficient accuracy as to give an independent… ▽ More All current estimates of the cosmic-ray (CR) ionization rate rely on assessments of the gas density along the probed sight lines. Until now, these have been based on observations of different tracers, with C$_2$ being the most widely used in diffuse molecular clouds for this purpose. However, three-dimensional dust extinction maps have recently reached sufficient accuracy as to give an independent measurement of the gas density on parsec scales. In addition, they allow us to identify the gas clumps along each sight line, thus localizing the regions where CR ionization is probed. We re-evaluate H$_3^+$ observations, which are often considered as the most reliable method to measure the H$_2$ ionization rate $ζ_{\rm H_2}$ in diffuse clouds. The peak density values derived from the extinction maps for 12 analyzed sight lines turn out to be, on average, an order of magnitude lower than the previous estimates, and agree with the values obtained from revised analysis of C$_2$ data. We use the extinction maps in combination with the 3D-PDR code to self-consistently compute the H$_3^+$ and H$_2$ abundances in the identified clumps for different values of $ζ_{\rm H_2}$. For each sight line, we obtain the optimum value by comparing the simulation results with observations. We show that $ζ_{\rm H_2}$ is systematically reduced with respect to the earlier estimates by a factor of $\approx 9$ on average, to $\approx6\times10^{-17}$ s$^{-1}$, primarily as a result of the density reduction. We emphasize that these results have profound consequences for all available measurements of the ionization rate. △ Less

Submitted 21 August, 2024; originally announced August 2024.

Comments: Accepted for publication in ApJ

The densities in diffuse and translucent molecular clouds: estimates from observations of C$_2$ and from 3-dimensional extinction maps

Authors: David A. Neufeld, Daniel E. Welty, Alexei V. Ivlev, Paola Caselli, Gordian Edenhofer, Nick Indriolo, Marta Obolentseva, Kedron Silsbee, Paule Sonnentrucker, Mark G. Wolfire

Abstract: Newly-computed collisional rate coefficients for the excitation of C$_2$ in collisions with H$_2$, presented recently by Najar and Kalugina (2020), are significantly larger than the values adopted previously in models for the excitation of the C$_2$ molecule, a widely used probe of the interstellar gas density. With these new rate coefficients, we have modeled the C$_2$ rotational distributions in… ▽ More Newly-computed collisional rate coefficients for the excitation of C$_2$ in collisions with H$_2$, presented recently by Najar and Kalugina (2020), are significantly larger than the values adopted previously in models for the excitation of the C$_2$ molecule, a widely used probe of the interstellar gas density. With these new rate coefficients, we have modeled the C$_2$ rotational distributions inferred from visible and ultraviolet absorption observations of electronic transitions of C$_2$ towards a collection of 46 nearby background sources. The inferred gas densities in the foreground interstellar clouds responsible for the observed C$_2$ absorption are a factor 4 to 7 smaller than those inferred previously, a direct reflection of the larger collisional rate coefficients computed by Najar and Kalugina (2020). These lower density estimates are generally in good agreement with the peak densities inferred from 3D extinction maps for the relevant sightlines. In cases where H$_3^+$ absorption has also been observed and used to estimate the cosmic-ray ionization rate (CRIR), our estimates of the latter will also decrease accordingly because the H$_3^+$ abundance is a function of the ratio of the CRIR to the gas density. △ Less

Submitted 1 September, 2024; v1 submitted 20 August, 2024; originally announced August 2024.

Comments: 46 pages, 11 figures, 1 figure set with 23 figures, accepted for publication in the Astrophysical Journal. Version 2: An error is corrected in Table 2, in which values listed in two of the columns were jumbled in Version 1. There are no other changes

Advancing spectroscopic understanding of HOCS$^+$: Laboratory investigations and astronomical implications

Authors: Valerio Lattanzi, Miguel Sanz-Novo, Víctor M. Rivilla, Mitsunori Araki, Hayley A Bunn, Jesús Martín-Pintado, Izaskun Jiménez-Serra, Paola Caselli

Abstract: Sulphur-bearing species play crucial roles in interstellar chemistry, yet their precise characterisation remains challenging. Here, we present laboratory experiments aimed at extending the high-resolution spectroscopy of protonated carbonyl sulphide (HOCS$^+$), a recently detected molecular ion in space. Using a frequency-modulated free-space absorption spectrometer, we detected rotational transit… ▽ More Sulphur-bearing species play crucial roles in interstellar chemistry, yet their precise characterisation remains challenging. Here, we present laboratory experiments aimed at extending the high-resolution spectroscopy of protonated carbonyl sulphide (HOCS$^+$), a recently detected molecular ion in space. Using a frequency-modulated free-space absorption spectrometer, we detected rotational transitions of HOCS$^+$ in an extended negative glow discharge with a mixture of H$_2$ and OCS, extending the high-resolution rotational characterisation of the cation well into the millimetre wave region (200-370 GHz). Comparisons with prior measurements and quantum chemical calculations revealed an overall agreement in the spectroscopic parameters. With the new spectroscopic dataset in hand, we re-investigated the observations of HOCS$^+$ towards G+0.693-0.027, which were initially based solely on K$_a$ = 0 lines contaminated by HNC$^{34}$S. This re-investigation enabled the detection of weak K$_a$ = 0 transitions, free from HNC$^{34}$S contamination. Our high-resolution spectroscopic characterisation also provides valuable insights for future millimetre and submillimetre astronomical observations of these species in different interstellar environments. In particular, the new high-resolution catalogue will facilitate the search for this cation in cold dark clouds, where very narrow line widths are typically observed. △ Less

Submitted 7 August, 2024; originally announced August 2024.

Comments: 9 pages, 3 figures, and 3 tables. Accepted for publication in Astronomy & Astrophysics

The GRAVITY young stellar object survey XIV : Investigating the magnetospheric accretion-ejection processes in S CrA N

Authors: GRAVITY Collaboration, H. Nowacki, K. Perraut, L. Labadie, J. Bouvier, C. Dougados, M. Benisty, J. A. Wojtczak, A. Soulain, E. Alecian, W. Brandner, A. Caratti o Garatti, R. Garcia Lopez, V. Ganci, J. Sánchez-Bermúdez, J. -P. Berger, G. Bourdarot, P. Caselli, Y. Clénet, R. Davies, A. Drescher, A. Eckart, F. Eisenhauer, M. Fabricius, H. Feuchtgruber , et al. (31 additional authors not shown)

Abstract: The dust- and gas-rich protoplanetary disks around young stellar systems play a key role in star and planet formation. While considerable progress has recently been made in probing these disks on large scales of a few tens of astronomical units (au), the central au needs to be more investigated. We aim at unveiling the physical processes at play in the innermost regions of the strongly accreting T… ▽ More The dust- and gas-rich protoplanetary disks around young stellar systems play a key role in star and planet formation. While considerable progress has recently been made in probing these disks on large scales of a few tens of astronomical units (au), the central au needs to be more investigated. We aim at unveiling the physical processes at play in the innermost regions of the strongly accreting T Tauri Star S CrA N by means of near-infrared interferometric observations. The K-band continuum emission is well reproduced with an azimuthally-modulated dusty ring. As the star alone cannot explain the size of this sublimation front, we propose that magnetospheric accretion is an important dust-heating mechanism leading to this continuum emission. The differential analysis of the Hydrogen Br$γ$ line is in agreement with radiative transfer models combining magnetospheric accretion and disk winds. Our observations support an origin of the Br$γ$ line from a combination of (variable) accretion-ejection processes in the inner disk region. △ Less

Submitted 5 August, 2024; originally announced August 2024.

FAUST. XVIII. Evidence for annular substructure in a very young Class 0 disk

Authors: M. J. Maureira, J. E. Pineda, H. B. Liu, L. Testi, D. Segura-Cox, C. Chandler, D. Johnstone, P. Caselli, G. Sabatini, Y. Aikawa, E. Bianchi, C. Codella, N. Cuello, D. Fedele, R. Friesen, L. Loinard, L. Podio, C. Ceccarelli, N. Sakai, S. Yamamoto

Abstract: When the planet formation process begins in the disks surrounding young stars is still an open question. Annular substructures such as rings and gaps in disks are intertwined with planet formation, and thus their presence or absence is commonly used to investigate the onset of this process. Current observations show a limited number of disks surrounding protostars exhibiting annular substructures,… ▽ More When the planet formation process begins in the disks surrounding young stars is still an open question. Annular substructures such as rings and gaps in disks are intertwined with planet formation, and thus their presence or absence is commonly used to investigate the onset of this process. Current observations show a limited number of disks surrounding protostars exhibiting annular substructures, all of them in the Class I stage. The lack of observed features in most of these sources may indicate a late emergence of substructures, but it could also be an artifact of these disks being optically thick. To mitigate the problem of optical depth, we investigate substructures within a very young Class 0 disk characterized by a low inclination using observations at longer wavelengths. We use 3 mm ALMA observations tracing dust emission at a resolution of 7 au to search for evidence of annular substructures in the disk around the deeply embedded Class 0 protostar Oph A SM1. The observations reveal a nearly face-on disk (i$\sim$16$^{\circ}$) extending up to 40 au. The radial intensity profile shows a clear deviation from a smooth profile near 30 au, which we interpret as the presence of either a gap at 28 au or a ring at 34 au with Gaussian widths of $σ=1.4^{+2.3}_{-1.2}$ au and $σ=3.9^{+2.0}_{-1.9}$ au, respectively. The 3 mm emission at the location of the possible gap or ring is determined to be optically thin, precluding the possibility that this feature in the intensity profile is due to the emission being optically thick. Annular substructures resembling those in the more evolved Class I and II disks could indeed be present in the Class 0 stage, earlier than previous observations suggested. Similar observations of embedded disks in which the high optical depth problem can be mitigated are clearly needed to better constrain the onset of substructures in the embedded stages. △ Less

Submitted 29 July, 2024; originally announced July 2024.

Comments: 12 pages, 7 figures, 2 tables. Accepted for publication in A&A Letters

Journal ref: A&A 689, L5 (2024)

Role of NH3 Binding Energy in the Early Evolution of Protostellar Cores

Authors: S. Kakkenpara Suresh, O. Sipila, P. Caselli, F. Dulieu

Abstract: NH$_{3}$(ammonia) plays a critical role in the chemistry of star and planet formation, yet uncertainties in its binding energy (BE) values complicate accurate estimates of its abundances. Recent research suggests a multi-binding energy approach, challenging the previous single-value notion. In this work, we use different values of NH$_{3}$ binding energy to examine its effects on the NH$_{3}$ abun… ▽ More NH$_{3}$(ammonia) plays a critical role in the chemistry of star and planet formation, yet uncertainties in its binding energy (BE) values complicate accurate estimates of its abundances. Recent research suggests a multi-binding energy approach, challenging the previous single-value notion. In this work, we use different values of NH$_{3}$ binding energy to examine its effects on the NH$_{3}$ abundances and, consequently, in the early evolution of protostellar cores. Using a gas-grain chemical network, we systematically vary the values of NH$_{3}$ binding energies in a model Class 0 protostellar core and study the effects of these binding energies on the NH$_{3}$ abundances. Our simulations indicate that abundance profiles of NH$_{3}$ are highly sensitive to the binding energy used, particularly in the warmer inner regions of the core. Higher binding energies lead to lower gas-phase NH$_{3}$ abundances, while lower values of binding energy have the opposite effect. Furthermore, this BE-dependent abundance variation of NH$_{3}$ significantly affects the formation pathways and abundances of key species such as HNC, HCN, and CN. Our tests also reveal that the size variation of the emitting region due to binding energy becomes discernible only with beam sizes of 10 arcsec or less. These findings underscore the importance of considering a range of binding energies in astrochemical models and highlight the need for higher resolution observations to better understand the subtleties of molecular cloud chemistry and star formation processes. △ Less

Submitted 25 July, 2024; originally announced July 2024.

FAUST XVII: Super deuteration in the planet forming system IRS 63 where the streamer strikes the disk

Authors: L. Podio, C. Ceccarelli, C. Codella, G. Sabatini, D. Segura-Cox, N. Balucani, A. Rimola, P. Ugliengo, C. J. Chandler, N. Sakai, B. Svoboda, J. Pineda, M. De Simone, E. Bianchi, P. Caselli, A. Isella, Y. Aikawa, M. Bouvier, E. Caux, L. Chahine, S. B. Charnley, N. Cuello, F. Dulieu, L. Evans, D. Fedele , et al. (33 additional authors not shown)

Abstract: Recent observations suggest that planets formation starts early, in protostellar disks of $\le10^5$ yrs, which are characterized by strong interactions with the environment, e.g., through accretion streamers and molecular outflows. To investigate the impact of such phenomena on disk physical and chemical properties it is key to understand what chemistry planets inherit from their natal environment… ▽ More Recent observations suggest that planets formation starts early, in protostellar disks of $\le10^5$ yrs, which are characterized by strong interactions with the environment, e.g., through accretion streamers and molecular outflows. To investigate the impact of such phenomena on disk physical and chemical properties it is key to understand what chemistry planets inherit from their natal environment. In the context of the ALMA Large Program Fifty AU STudy of the chemistry in the disk/envelope system of Solar-like protostars (FAUST), we present observations on scales from ~1500 au to ~60 au of H$_2$CO, HDCO, and D$_2$CO towards the young planet-forming disk IRS~63. H$_2$CO probes the gas in the disk as well as in a large scale streamer (~1500 au) impacting onto the South-East (SE) disk side. We detect for the first time deuterated formaldehyde, HDCO and D$_2$CO, in a planet-forming disk, and HDCO in the streamer that is feeding it. This allows us to estimate the deuterium fractionation of H$_2$CO in the disk: [HDCO]/[H$_2$CO]$\sim0.1-0.3$ and [D$_2$CO]/[H$_2$CO]$\sim0.1$. Interestingly, while HDCO follows the H$_2$CO distribution in the disk and in the streamer, the distribution of D$_2$CO is highly asymmetric, with a peak of the emission (and [D]/[H] ratio) in the SE disk side, where the streamer crashes onto the disk. In addition, D$_2$CO is detected in two spots along the blue- and red-shifted outflow. This suggests that: (i) in the disk, HDCO formation is dominated by gas-phase reactions similarly to H$_2$CO, while (ii) D$_2$CO was mainly formed on the grain mantles during the prestellar phase and/or in the disk itself, and is at present released in the gas-phase in the shocks driven by the streamer and the outflow. These findings testify on the key role of streamers in the build-up of the disk both concerning the final mass available for planet formation and its chemical composition. △ Less

Submitted 5 July, 2024; originally announced July 2024.

Comments: 12 pages, 10 figures, accepted for publication on A&A

Journal ref: A&A 688, L22 (2024)

Unveiling the role of magnetic fields in an accreting filament onto a young protocluster

Authors: Farideh S. Tabatabaei, Elena Redaelli, Daniele Galli, Paola Caselli, Gabriel A. P. Franco, Ana Duarte-Cabral, Marco Padovani

Abstract: In order to develop a more comprehensive picture of star formation, it is essential to understand the physical relationship between dense cores and the filaments embedding them. There is evidence that magnetic fields play a crucial role in this context. We aim to understand how magnetic fields influence the properties and kinematics of an isolated filament located east of the Barnard 59 clump, bel… ▽ More In order to develop a more comprehensive picture of star formation, it is essential to understand the physical relationship between dense cores and the filaments embedding them. There is evidence that magnetic fields play a crucial role in this context. We aim to understand how magnetic fields influence the properties and kinematics of an isolated filament located east of the Barnard 59 clump, belonging to the Pipe Nebula. We use near infrared polarization observations to determine the magnetic field configuration, and we apply the Davis Chandrasekhar Fermi method to infer the magnetic field strength in the plane of the sky. Furthermore, we use complementary data from the James Clerk Maxwell Submillimetre Telescope (JCMT) of C18O and 13CO J=3-2 transition to determine the filament's kinematics. Finally, we model the radial density profile of the filament with polytropic cylindrical models. Our results indicate that the filament is stable to radial collapse and is radially supported by agents other than thermal pressure. In addition, based on previous observations of emission lines on this source, we suggest that gas is flowing toward the hub, while C18O (3-2) non-thermal motions indicate that the cloud is in a quiescent state. △ Less

Submitted 26 June, 2024; originally announced June 2024.

Comments: 12 pages, accepted at A&A journal

Journal ref: A&A 688, A98 (2024)

PRODIGE -- Planet-forming disks in Taurus with NOEMA

Authors: R. Franceschi, Th. Henning, G. V. Smirnov-Pinchukov, D. A. Semenov, K. Schwarz, A. Dutrey, E. Chapillon, U. Gorti, S. Guilloteau, V. Piétu, S. van Terwisga, L. Bouscasse, P. Caselli, G. Gieser, T. -H. Hsieh, A. Lopez-Sepulcre, D. M. Segura-Cox, J. E. Pineda, M. J. Maureira, M. T. Valdivia-Mena

Abstract: We aim to constrain the gas density and temperature distributions as well as gas masses in several T Tauri protoplanetary disks located in Taurus. We use the 12CO, 13CO, and C18O (2-1) isotopologue emission observed at 0.9 with the IRAM NOrthern Extended Millimeter Array (NOEMA) as part of the MPG-IRAM Observatory Program PRODIGE (PROtostars and DIsks: Global Evolution PIs: P. Caselli & Th. Hennin… ▽ More We aim to constrain the gas density and temperature distributions as well as gas masses in several T Tauri protoplanetary disks located in Taurus. We use the 12CO, 13CO, and C18O (2-1) isotopologue emission observed at 0.9 with the IRAM NOrthern Extended Millimeter Array (NOEMA) as part of the MPG-IRAM Observatory Program PRODIGE (PROtostars and DIsks: Global Evolution PIs: P. Caselli & Th. Henning). Our sample consists of Class II disks with no evidence of strong radial substructures. We use thesedata to constrain the thermal and chemical structure of these disks through theoretical models for gas emission. To fit the combined optically thick and thin CO line data in Fourier space, we developed the DiskCheF code, which includes the parameterized disk physical structure, machine-learning (ML) accelerated chemistry, and the RADMC-3D line radiative transfer module. A key novelty of DiskCheF is the fast and feasible ML-based chemistry trained on the extended grid of the disk physical-chemical models precomputed with the ANDES2 code. This ML approach allows complex chemical kinetics models to be included in a time-consuming disk fitting without the need to run a chemical code. We present a novel approach to incorporate chemistry into disk modeling without the need to explicitly calculate a chemical network every time. Using this new disk modeling tool, we successfully fit the 12CO, 13CO, and C18O (2-1) data from the CI, CY, DL, DM, DN, and IQ Tau disks. The combination of optically thin and optically thick CO lines allows us to simultaneously constrain the disk temperature and mass distribution, and derive the CO-based gas masses. These values are in reasonable agreement with the disk dust masses rescaled by a factor of 100 as well as with other indirect gas measurements. △ Less

Submitted 24 June, 2024; originally announced June 2024.

Comments: 25 pages, 19 figures

Overtone Transition $2ν_1$ of $\text{HCO}^+$ and $\text{HOC}^+$: Origin, Radiative Lifetime, Collisional Quenching

Authors: Miguel Jiménez-Redondo, Liliia Uvarova, Petr Dohnal, Miroslava Kassayová, Paola Caselli, Pavol Jusko

Abstract: We present spectra of the first overtone vibration transition of $\text{C-H}$/$\text{O-H}$ stretch ($2ν_1$) in $\text{HCO}^+$ and $\text{HOC}^+$, recorded using a laser induced reaction action scheme inside a cryogenic 22 pole radio frequency trap. Band origins have been located at 6078.68411(19) and 6360.17630(26) $\text{cm}^{-1}$, respectively. We introduce a technique based on mass selective ej… ▽ More We present spectra of the first overtone vibration transition of $\text{C-H}$/$\text{O-H}$ stretch ($2ν_1$) in $\text{HCO}^+$ and $\text{HOC}^+$, recorded using a laser induced reaction action scheme inside a cryogenic 22 pole radio frequency trap. Band origins have been located at 6078.68411(19) and 6360.17630(26) $\text{cm}^{-1}$, respectively. We introduce a technique based on mass selective ejection from the ion trap for recording background free action spectra. Varying the number density of the neutral action scheme reactant ($\text{CO}_2$ and Ar, respectively) and collisional partner reactant inside the ion trap, permitted us to estimate the radiative lifetime of the state to be 1.53(34) and 1.22(34) ms, respectively, and the collisional quenching rates of $\text{HCO}^+$ ($2ν_1$) with He, H$_2$, and N$_2$. △ Less

Submitted 13 June, 2024; originally announced June 2024.

Comments: To appear in ChemPhysChem

High contrast at short separation with VLTI/GRAVITY: Bringing Gaia companions to light

Authors: N. Pourré, T. O. Winterhalder, J. -B. Le Bouquin, S. Lacour, A. Bidot, M. Nowak, A. -L. Maire, D. Mouillet, C. Babusiaux, J. Woillez, R. Abuter, A. Amorim, R. Asensio-Torres, W. O. Balmer, M. Benisty, J. -P. Berger, H. Beust, S. Blunt, A. Boccaletti, M. Bonnefoy, H. Bonnet, M. S. Bordoni, G. Bourdarot, W. Brandner, F. Cantalloube , et al. (151 additional authors not shown)

Abstract: Since 2019, GRAVITY has provided direct observations of giant planets and brown dwarfs at separations of down to 95 mas from the host star. Some of these observations have provided the first direct confirmation of companions previously detected by indirect techniques (astrometry and radial velocities). We want to improve the observing strategy and data reduction in order to lower the inner working… ▽ More Since 2019, GRAVITY has provided direct observations of giant planets and brown dwarfs at separations of down to 95 mas from the host star. Some of these observations have provided the first direct confirmation of companions previously detected by indirect techniques (astrometry and radial velocities). We want to improve the observing strategy and data reduction in order to lower the inner working angle of GRAVITY in dual-field on-axis mode. We also want to determine the current limitations of the instrument when observing faint companions with separations in the 30-150 mas range. To improve the inner working angle, we propose a fiber off-pointing strategy during the observations to maximize the ratio of companion-light-to-star-light coupling in the science fiber. We also tested a lower-order model for speckles to decouple the companion light from the star light. We then evaluated the detection limits of GRAVITY using planet injection and retrieval in representative archival data. We compare our results to theoretical expectations. We validate our observing and data-reduction strategy with on-sky observations; first in the context of brown dwarf follow-up on the auxiliary telescopes with HD 984 B, and second with the first confirmation of a substellar candidate around the star Gaia DR3 2728129004119806464. With synthetic companion injection, we demonstrate that the instrument can detect companions down to a contrast of $8\times 10^{-4}$ ($Δ\mathrm{K}= 7.7$ mag) at a separation of 35 mas, and a contrast of $3\times 10^{-5}$ ($Δ\mathrm{K}= 11$ mag) at 100 mas from a bright primary (K<6.5), for 30 min exposure time. With its inner working angle and astrometric precision, GRAVITY has a unique reach in direct observation parameter space. This study demonstrates the promising synergies between GRAVITY and Gaia for the confirmation and characterization of substellar companions. △ Less

Submitted 6 June, 2024; originally announced June 2024.

Comments: 16 pages, 14 figures. Submitted to A&A

A low cosmic-ray ionisation rate in the prestellar core Ophiuchus/H-MM1. Mapping of the molecular ions ortho-H2D+, N2H+, and DCO+

Authors: Jorma Harju, Charlotte Vastel, Olli Sipilae, Elena Redaelli, Paola Caselli, Jaime E. Pineda, Arnaud Belloche, Friedrich Wyrowski

Abstract: (abridged) We have mapped the prestellar core H-MM1 in Ophiuchus in rotational lines of ortho-H2D+ (oH2D+), N2H+, and DCO+ at the wavelength 0.8 mm with the Large APEX sub-Millimeter Array (LAsMA) multibeam receiver of the Atacama Pathfinder EXperiment (APEX) telescope. We also ran a series of chemistry models to predict the abundance distributions of the observed molecules, and to estimate the ef… ▽ More (abridged) We have mapped the prestellar core H-MM1 in Ophiuchus in rotational lines of ortho-H2D+ (oH2D+), N2H+, and DCO+ at the wavelength 0.8 mm with the Large APEX sub-Millimeter Array (LAsMA) multibeam receiver of the Atacama Pathfinder EXperiment (APEX) telescope. We also ran a series of chemistry models to predict the abundance distributions of the observed molecules, and to estimate the effect of the cosmic-ray ionisation rate on their abundances. The three line maps show different distributions. The oH2D+ map is extended and outlines the general structure of the core, while N2H+ mainly shows the density maxima, and the DCO+ emission peaks are shifted towards one edge of the core where a region of enhanced desorption has been found previously. According to the chemical simulation, the fractional oH2D+ abundance remains relatively high in the centre of the core, and its column density correlates strongly with the cosmic-ray ionisation rate. Simulated line maps constrain the cosmic-ray ionisation rate per hydrogen molecule to be low, between 5e-18/s and 1e-17/s in the H-MM1 core. This estimate agrees with the gas temperature measured in the core. Modelling line emission of oH2D+ provides a straightforward method of determining the cosmic-ray ionisation rate in dense clouds, where the primary ion, H3+, is not observable. △ Less

Submitted 30 May, 2024; originally announced May 2024.

Comments: accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 688, A117 (2024)

Multiple chemical tracers finally unveil the intricate NGC\,1333 IRAS\,4A outflow system. FAUST XVI

Authors: Layal Chahine, Cecilia Ceccarelli, Marta De Simone, Claire J. Chandler, Claudio Codella, Linda Podio, Ana López-Sepulcre, Nami Sakai, Laurent Loinard, Mathilde Bouvier, Paola Caselli, Charlotte Vastel, Eleonora Bianchi, Nicolás Cuello, Francesco Fontani, Doug Johnstone, Giovanni Sabatini, Tomoyuki Hanawa, Ziwei E. Zhang, Yuri Aikawa, Gemma Busquet, Emmanuel Caux, Aurore Durán, Eric Herbst, François Ménard , et al. (32 additional authors not shown)

Abstract: The exploration of outflows in protobinary systems presents a challenging yet crucial endeavour, offering valuable insights into the dynamic interplay between protostars and their evolution. In this study, we examine the morphology and dynamics of jets and outflows within the IRAS\,4A protobinary system. This analysis is based on ALMA observations of SiO(5--4), H$_2$CO(3$_{0,3}$--2$_{0,3}$), and H… ▽ More The exploration of outflows in protobinary systems presents a challenging yet crucial endeavour, offering valuable insights into the dynamic interplay between protostars and their evolution. In this study, we examine the morphology and dynamics of jets and outflows within the IRAS\,4A protobinary system. This analysis is based on ALMA observations of SiO(5--4), H$_2$CO(3$_{0,3}$--2$_{0,3}$), and HDCO(4$_{1,4}$--3$_{1,3}$) with a spatial resolution of $\sim$150\,au. Leveraging an astrochemical approach involving the use of diverse tracers beyond traditional ones has enabled the identification of novel features and a comprehensive understanding of the broader outflow dynamics. Our analysis reveals the presence of two jets in the redshifted emission, emanating from IRAS\,4A1 and IRAS\,4A2, respectively. Furthermore, we identify four distinct outflows in the region for the first time, with each protostar, 4A1 and 4A2, contributing to two of them. We characterise the morphology and orientation of each outflow, challenging previous suggestions of bends in their trajectories. The outflow cavities of IRAS\,4A1 exhibit extensions of 10$''$ and 13$''$ with position angles (PA) of 0$^{\circ}$ and -12$^{\circ}$, respectively, while those of IRAS\,4A2 are more extended, spanning 18$''$ and 25$''$ with PAs of 29$^{\circ}$ and 26$^{\circ}$. We propose that the misalignment of the cavities is due to a jet precession in each protostar, a notion supported by the observation that the more extended cavities of the same source exhibit lower velocities, indicating they may stem from older ejection events. △ Less

Submitted 21 May, 2024; originally announced May 2024.

Characterization of carbon dioxide on Ganymede and Europa supported by experiments: Effects of temperature, porosity, and mixing with water

Authors: L. Schiltz, B. Escribano, G. M. Muñoz Caro, S. Cazaux, C. del Burgo Olivares, H. Carrascosa, I. Boszhuizen, C. González Díaz, Y. -J. Chen, B. M. Giuliano, P. Caselli

Abstract: The surfaces of icy moons are primarily composed of water ice that can be mixed with other compounds, such as carbon dioxide. The carbon dioxide (CO$_2)$ stretching fundamental band observed on Europa and Ganymede appears to be a combination of several bands that are shifting location from one moon to another. We investigate the cause of the observed shift in the CO$_2$ stretching absorption band… ▽ More The surfaces of icy moons are primarily composed of water ice that can be mixed with other compounds, such as carbon dioxide. The carbon dioxide (CO$_2)$ stretching fundamental band observed on Europa and Ganymede appears to be a combination of several bands that are shifting location from one moon to another. We investigate the cause of the observed shift in the CO$_2$ stretching absorption band experimentally. We also explore the spectral behaviour of CO$_2$ ice by varying the temperature and concentration.} %H$_2$O:CO$_2$ deposition ratios. We analyzed pure CO$_2$ ice and ice mixtures deposited at 10 K under ultra-high vacuum conditions using Fourier-transform infrared (FTIR) spectroscopy and temperature programmed desorption (TPD) experiments. Laboratory ice spectra were compared to JWST observation of Europa's and Ganymede's leading hemispheres. The simulated IR spectra were calculated using density functional theory (DFT) methods, exploring the effect of porosity in CO$_2$ ice. Pure CO$_2$ and CO$_2$-water ice show distinct spectral changes and desorption behaviours at different temperatures, revealing intricate CO$_2$ and H$_2$O interactions. The number of discernible peaks increases from two in pure CO$_2$ to three in CO$_2$-water mixtures. △ Less

Submitted 17 May, 2024; originally announced May 2024.

Comments: Accepted for publication in Astronomy & Astrophysics

FAUST XV. A disk wind mapped by CH$_3$OH and SiO in the inner 300 au of the NGC 1333 IRAS 4A2 protostar

Authors: M. De Simone, L. Podio, L. Chahine, C. Codella, C. J. Chandler, C. Ceccarelli, A. Lopez-Sepulcre, L. Loinard, B. Svoboda, N. Sakai, D. Johnstone, F. Menard, Y. Aikawa, M. Bouvier, G. Sabatini, A. Miotello, C. Vastel, N. Cuello, E. Bianchi, P. Caselli, E. Caux, T. Hanawa, E. Herbst, D. Segura-Cox, Z. Zhang , et al. (1 additional authors not shown)

Abstract: Context. Understanding the connection between outflows, winds, accretion and disks in the inner protostellar regions is crucial for comprehending star and planet formation process. Aims. We aim to we explore the inner 300 au of the protostar IRAS 4A2 as part of the ALMA FAUST Large Program. Methods. We analysed the kinematical structures of SiO and CH$_3$OH emission with 50 au resolution. Results.… ▽ More Context. Understanding the connection between outflows, winds, accretion and disks in the inner protostellar regions is crucial for comprehending star and planet formation process. Aims. We aim to we explore the inner 300 au of the protostar IRAS 4A2 as part of the ALMA FAUST Large Program. Methods. We analysed the kinematical structures of SiO and CH$_3$OH emission with 50 au resolution. Results. The emission arises from three zones: i) a very compact and unresolved region ($<$50 au) dominated by the ice sublimation zone, at $\pm$1.5 km s$^{-1}$ with respect to vsys, traced by methanol; ii) an intermediate region (between 50 au and 150 au) traced by both SiO and CH$_3$OH, between 2 and 6 km s$^{-1}$ with respect to vsys, with an inverted velocity gradient (with respect to the large scale emission), whose origin is not clear; iii) an extended region ($>$150 au) traced by SiO, above 7 km s$^{-1}$ with respect to vsys, and dominated by the outflow. In the intermediate region we estimated a CH$_3$OH/SiO abundance ratio of about 120-400 and a SiO/H$_2$ abundance of 10$^{-8}$. We explored various possibilities to explain the origin of this region such as, rotating disk/inner envelope, jet on the plane of the sky/precessing, wide angle disk wind. Conclusions. We propose that CH$_3$OH and SiO in the inner 100 au probe the base of a wide-angle disk wind. The material accelerated in the wind crosses the plane of the sky, giving rise to the observed inverted velocity gradient, and sputtering the grain mantles and cores releasing CH$_3$OH and SiO. This is the first detection of a disk wind candidate in SiO, and the second ever in CH$_3$OH. △ Less

Submitted 30 April, 2024; originally announced April 2024.

Comments: Accepted in A&A Letters

Journal ref: A&A 686, L13 (2024)

Quasi-equilibrium chemical evolution in starless cores

Authors: Jonathan Rawlings, Eric Keto, Paola Caselli

Abstract: The chemistry of H2O, CO and other small molecular species in an isolated pre-stellar core, L1544, has been assessed in the context of a comprehensive gas-grain chemical model, coupled to an empirically constrained physical/dynamical model. Our main findings are (i) that the chemical network remains in near equilibrium as the core evolves towards star formation and the molecular abundances change… ▽ More The chemistry of H2O, CO and other small molecular species in an isolated pre-stellar core, L1544, has been assessed in the context of a comprehensive gas-grain chemical model, coupled to an empirically constrained physical/dynamical model. Our main findings are (i) that the chemical network remains in near equilibrium as the core evolves towards star formation and the molecular abundances change in response to the evolving physical conditions. The gas-phase abundances at any time can be calculated accurately with equilibrium chemistry, and the concept of chemical clocks is meaningless in molecular clouds with similar conditions and dynamical time scales, and (ii) A comparison of the results of complex and simple chemical networks indicates that the abundances of the dominant oxygen and carbon species, H2O, CO, C, and C+ are reasonably approximated by simple networks. In chemical equilibrium, the time-dependent differential terms vanish and a simple network reduces to a few algebraic equations. This allows rapid calculation of the abundances most responsible for spectral line radiative cooling in molecular clouds with long dynamical time scales. The dust ice mantles are highly structured and the ice layers retain a memory of the gas-phase abundances at the time of their deposition. A complex (gas-phase and gas-grain) chemical structure therefore exists, with cosmic-ray induced processes dominating in the inner regions. The inferred H2O abundance profiles for L1544 require that the outer parts of the core and also any medium exterior to the core are essentially transparent to the interstellar radiation field. △ Less

Submitted 24 April, 2024; originally announced April 2024.

Comments: 19 pages, 6 figures

A broad linewidth, compact, millimeter-bright molecular emission line source near the Galactic Center

Authors: Adam Ginsburg, John Bally, Ashley T. Barnes, Cara Battersby, Nazar Budaiev, Natalie O. Butterfield, Paola Caselli, Laura Colzi, Katarzyna M. Dutkowska, Pablo García, Savannah Gramze, Jonathan D. Henshaw, Yue Hu, Desmond Jeff, Izaskun Jiménez-Serra, Jens Kauffmann, Ralf S. Klessen, Emily M. Levesque, Steven N. Longmore, Xing Lu, Elisabeth A. C. Mills, Mark R. Morris, Francisco Nogueras-Lara, Tomoharu Oka, Jaime E. Pineda , et al. (15 additional authors not shown)

Abstract: A compact source, G0.02467-0.0727, was detected in ALMA \threemm observations in continuum and very broad line emission. The continuum emission has a spectral index $α\approx3.3$, suggesting that the emission is from dust. The line emission is detected in several transitions of CS, SO, and SO$_2$ and exhibits a line width FWHM $\approx160$ \kms. The line profile appears Gaussian. The emission is w… ▽ More A compact source, G0.02467-0.0727, was detected in ALMA \threemm observations in continuum and very broad line emission. The continuum emission has a spectral index $α\approx3.3$, suggesting that the emission is from dust. The line emission is detected in several transitions of CS, SO, and SO$_2$ and exhibits a line width FWHM $\approx160$ \kms. The line profile appears Gaussian. The emission is weakly spatially resolved, coming from an area on the sky $\lesssim1"$ in diameter ($\lesssim10^4$ AU at the distance of the Galactic Center; GC). The centroid velocity is $v_{LSR}\approx40$-$50$ \kms, which is consistent with a location in the Galactic Center. With multiple SO lines detected, and assuming local thermodynamic equilibrium (LTE) conditions, $T_\mathrm{LTE} = 13$ K, which is colder than seen in typical GC clouds, though we cannot rule out low-density, subthermally excited, warmer gas. Despite the high velocity dispersion, no emission is observed from SiO, suggesting that there are no strong ($\gtrsim10~\mathrm{km~s}^{-1}$) shocks in the molecular gas. There are no detections at other wavelengths, including X-ray, infrared, and radio. We consider several explanations for the Millimeter Ultra-Broad Line Object (MUBLO), including protostellar outflow, explosive outflow, collapsing cloud, evolved star, stellar merger, high-velocity compact cloud, intermediate mass black hole, and background galaxy. Most of these conceptual models are either inconsistent with the data or do not fully explain it. The MUBLO is, at present, an observationally unique object. △ Less

Submitted 1 May, 2024; v1 submitted 11 April, 2024; originally announced April 2024.

Comments: Accepted to ApJL

Four-of-a-kind? Comprehensive atmospheric characterisation of the HR 8799 planets with VLTI/GRAVITY

Authors: E. Nasedkin, P. Mollière, S. Lacour, M. Nowak, L. Kreidberg, T. Stolker, J. J. Wang, W. O. Balmer, J. Kammerer, J. Shangguan, R. Abuter, A. Amorim, R. Asensio-Torres, M. Benisty, J. -P. Berger, H. Beust, S. Blunt, A. Boccaletti, M. Bonnefoy, H. Bonnet, M. S. Bordoni, G. Bourdarot, W. Brandner, F. Cantalloube, P. Caselli , et al. (73 additional authors not shown)

Abstract: With four companions at separations from 16 to 71 au, HR 8799 is a unique target for direct imaging, presenting an opportunity for the comparative study of exoplanets with a shared formation history. Combining new VLTI/GRAVITY observations obtained within the ExoGRAVITY program with archival data, we perform a systematic atmospheric characterisation of all four planets. We explore different levels… ▽ More With four companions at separations from 16 to 71 au, HR 8799 is a unique target for direct imaging, presenting an opportunity for the comparative study of exoplanets with a shared formation history. Combining new VLTI/GRAVITY observations obtained within the ExoGRAVITY program with archival data, we perform a systematic atmospheric characterisation of all four planets. We explore different levels of model flexibility to understand the temperature structure, chemistry and clouds of each planet using both petitRADTRANS atmospheric retrievals and fits to self-consistent radiative-convective equilibrium models. Using Bayesian Model Averaging to combine multiple retrievals, we find that the HR 8799 planets are highly enriched in metals, with [M/H] $\gtrsim$1, and have stellar to super-stellar C/O ratios. The C/O ratio increases with increasing separation from $0.55^{+0.12}_{-0.10}$ for d to $0.78^{+0.03}_{-0.04}$ for b, with the exception of the innermost planet which has a C/O ratio of $0.87\pm0.03$. By retrieving a quench pressure and using a disequilibrium chemistry model we derive vertical mixing strengths compatible with predictions for high-metallicity, self-luminous atmospheres. Bayesian evidence comparisons strongly favour the presence of HCN in HR 8799 c and e, as well as CH$_{4}$ in HR 8799 c, with detections at $>5σ$ confidence. All of the planets are cloudy, with no evidence for patchiness. The clouds of c, d and e are best fit by silicate clouds lying above a deep iron cloud layer, while the clouds of the cooler HR 8799 b are more likely composed of Na$_{2}$S. With well defined atmospheric properties, future exploration of this system is well positioned to unveil further detail in these planets, extending our understanding of the composition, structure, and formation history of these siblings. △ Less

Submitted 17 July, 2024; v1 submitted 4 April, 2024; originally announced April 2024.

Comments: 45 pages, 25 figures, accepted by Astronomy & Astrophysics

Journal ref: A&A 687, A298 (2024)

Comprehensive laboratory constraints on thermal desorption of interstellar ice analogues

Authors: F. Kruczkiewicz, F. Dulieu, A. V. Ivlev, P. Caselli, B. M. Giuliano, C. Ceccarelli, P. Theulé

Abstract: To explain grain growth and destruction in warm media, ice mantle formation and sublimation in cold media, and gas line emission spectroscopy, astrochemical models must mimic the gas--solid abundance ratio. Ice-sublimation mechanisms determine the position of snow lines and the nature of gas emitted by and locked inside planetary bodies in star-forming regions. To interpret observations from the i… ▽ More To explain grain growth and destruction in warm media, ice mantle formation and sublimation in cold media, and gas line emission spectroscopy, astrochemical models must mimic the gas--solid abundance ratio. Ice-sublimation mechanisms determine the position of snow lines and the nature of gas emitted by and locked inside planetary bodies in star-forming regions. To interpret observations from the interplanetary and extragalactic interstellar mediums, gas phase abundances must be modelled correctly. This study presents comprehensive thermal desorption data for interstellar ice analogues, aiming to refine astrochemical models by generating a set of benchmarks to evaluate both the kinetics and thermodynamics in astrochemical models. Our experiments focused on temperature-programmed desorption of pure and mixed ices, including Ar, CO, CO2, NH3, CH3OH, and H2O, under ultrahigh vacuum (1 x 10^-10 hPa) and low temperatures (10 K). Each experiment includes the experimental parameters, ice desorption kinetics for pure species, and the desorption yield (gas--solid ratio) for ice mixtures. From the desorption yields, we find common trends in the trapping of molecules when their abundance is compared to water: compact amorphous water ices are capable of trapping up to 20 % of volatiles (Ar, CO, and CO2), ~ 3 % of CH3OH, and ~ 5% NH3 in relation to the water content within the ice matrix; ammonium formate is not trapped in the water ice films, and compact amorphous water ice formed in situ has similar trapping capabilities to a compact amorphous water ice deposited using molecular beams. Our results highlight the limited trapping capacity of compact amorphous water ice for gases, crucial for understanding the formation of interstellar complex organic molecules. △ Less

Submitted 3 April, 2024; originally announced April 2024.

Comments: Accepted for publication (A&A); 21 pages, 12 figures

Probing the Physics of Star-Formation (ProPStar): II. The first systematic search for streamers toward protostars

Authors: María Teresa Valdivia-Mena, Jaime E. Pineda, Paola Caselli, Dominique M. Segura-Cox, Anika Schmiedeke, Silvia Spezzano, Stella Offner, Alexei V. Ivlev, Michael Küffmeier, Nichol Cunningham, Roberto Neri, María José Maureira

Abstract: The detection of narrow channels of accretion toward protostellar disks, known as streamers, have increased in number in the last few years. However, it is unclear if streamers are a common feature around protostars that were previously missed, or if they are a rare phenomenon. Our goals are to obtain the incidence of streamers toward a region of clustered star formation and to trace the origins o… ▽ More The detection of narrow channels of accretion toward protostellar disks, known as streamers, have increased in number in the last few years. However, it is unclear if streamers are a common feature around protostars that were previously missed, or if they are a rare phenomenon. Our goals are to obtain the incidence of streamers toward a region of clustered star formation and to trace the origins of their gas, to determine if they originate within the filamentary structure of molecular clouds or from beyond. We used combined observations of the nearby NGC 1333 star-forming region, carried out with the NOEMA interferometer and the IRAM 30m single dish. Our observations cover the area between the IRAS 4 and SVS 13 systems. We traced the chemically fresh gas within NGC 1333 with HC3N molecular gas emission and the structure of the fibers in this region with N2H+ emission. We fit multiple velocity components in both maps and used clustering algorithms to recover velocity-coherent structures. We find streamer candidates toward 7 out of 16 young stellar objects within our field of view. This represents an incidence of approximately 40\% of young stellar objects with streamer candidates when looking at a clustered star forming region. The incidence increases to about 60\% when we considered only embedded protostars. All streamers are found in HC3N emission. Given the different velocities between HC3N and N2H+ emission, and the fact that, by construction, N2H+ traces the fiber structure, we suggest that the gas that forms the streamers comes from outside the fibers. This implies that streamers can connect cloud material that falls to the filaments with protostellar disk scales. △ Less

Submitted 2 April, 2024; originally announced April 2024.

Comments: 22 pages, 15 figures, accepted for publication in Astronomy and Astrophysics. GitHub repository with code available at https://github.com/tere-valdivia/propstar_gasflow

Journal ref: A&A 687, A71 (2024)

FAUST XIII. Dusty cavity and molecular shock driven by IRS7B in the Corona Australis cluster

Authors: G. Sabatini, L. Podio, C. Codella, Y. Watanabe, M. De Simone, E. Bianchi, C. Ceccarelli, C. J. Chandler, N. Sakai, B. Svoboda, L. Testi, Y. Aikawa, N. Balucani, M. Bouvier, P. Caselli, E. Caux, L. Chahine, S. Charnley, N. Cuello, F. Dulieu, L. Evans, D. Fedele, S. Feng, F. Fontani, T. Hama , et al. (32 additional authors not shown)

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

Submitted 2 April, 2024; v1 submitted 26 March, 2024; originally announced March 2024.

Comments: 12 pages, 8 figures, 3 tables. Accepted Letter in Astronomy & Astrophysics

PRODIGE -- Envelope to Disk with NOEMA III. The origin of complex organic molecule emission in SVS13A

Authors: T. -H. Hsieh, J. E. Pineda, D. M. Segura-Cox, P. Caselli, M. T. Valdivia-Mena, C. Gieser, M. J. Maureira, A. Lopez-Sepulcre, L. Bouscasse, R. Neri, Th. Möller, A. Dutrey, A. Fuente, D. Semenov, E. Chapillon, N. Cunningham, Th. Henning, V. Pietu, I. Jimenez-Serra, S. Marino, C. Ceccarelli

Abstract: Complex Organic Molecules (COMs) have been found toward low-mass protostars but the origins of the COM emission are still unclear. It can be associated with, for example, hot corinos, outflows, and/or accretion shock/disk atmosphere. We have conducted NOEMA observations toward SVS13A from the PROtostars & DIsks: Global Evolution (PRODIGE) program. Our previous \ce{DCN} observations reveal a possib… ▽ More Complex Organic Molecules (COMs) have been found toward low-mass protostars but the origins of the COM emission are still unclear. It can be associated with, for example, hot corinos, outflows, and/or accretion shock/disk atmosphere. We have conducted NOEMA observations toward SVS13A from the PROtostars & DIsks: Global Evolution (PRODIGE) program. Our previous \ce{DCN} observations reveal a possible infalling streamer, which may affect the chemistry of the central protobinary by inducing accretion outbursts and/or shocked gas. Here, we further analyze six O-bearing COMs: CH3OH, aGg'-(CH2OH)2, C2H5OH, CH2(OH)CHO, CH3CHO, and CH3OCHO. Although the COM emission is not spatially resolved, we constrain the source sizes to $\lesssim0.3-0.4$ arcsec (90$-$120 au) by conducting uv-domain Gaussian fitting. Interestingly, the high-spectral resolution data reveal complex line profiles with multiple peaks showing differences between these six O-bearing COMs. The LTE fitting unveils differences in excitation temperatures and emitting areas among these COMs. We further conduct multiple-velocity-component LTE fitting to decompose the line emission into different kinematic components. Up to 6 velocity components are found from the LTE modeling. The temperature, column density, and source size of these components from each COM are obtained. We find a variety in excitation temperatures ($100-500$ K) and source sizes (D$\sim10-70$ au) from these kinematic components from different COMs. The emission of each COM can trace several components and different COMs most likely trace different regions. Given this complex structure, we suggest that the central region is inhomogeneous and unlikely to be heated by only protostellar radiation. We conclude that accretion shocks induced by the large-scale infalling streamer likely exist and contribute to the complexity of the COM emission. △ Less

Submitted 25 March, 2024; originally announced March 2024.

Comments: 29 pages, 18 figures, accepted to A&A

Journal ref: A&A 686, A289 (2024)

Millimeter and submillimeter spectroscopy of the deuterated molecular ion SD+

Authors: Mitsunori Araki, Valerio Lattanzi, Christian P. Endres, Paula Caselli

Abstract: Seven rotational and fine-structure transitions of the deuterated molecular ion SD+ in the X 3S- ground electronic state have been measured in the 271-863 GHz region in the laboratory. This ion has been produced by DC-glow discharge using a mixture of D2S and argon in a free space cell in a temperature range of -140 to -160C. The rotational, centrifugal distortion, spin-spin interaction, and hyper… ▽ More Seven rotational and fine-structure transitions of the deuterated molecular ion SD+ in the X 3S- ground electronic state have been measured in the 271-863 GHz region in the laboratory. This ion has been produced by DC-glow discharge using a mixture of D2S and argon in a free space cell in a temperature range of -140 to -160C. The rotational, centrifugal distortion, spin-spin interaction, and hyperfine constants have been determined; the standard deviation of the residuals in the fitting is 109 kHz. The set of obtained spectroscopic parameters provides a list of accurate sub-millimeter rest frequencies of SD+ for astronomical detection. We have investigated lines of SD+ toward the quasar PKS 1830-211 using the ALMA archive, as the z = 0.89 molecular absorber exists in front of this quasar. A data set covering the 297 GHz region includes the N_J = 2_3-1_2 transition at 561 GHz due to redshift, providing an upper limit of the column density Ntot = 3 x 10^12 cm-2 for SD+. △ Less

Submitted 22 March, 2024; originally announced March 2024.

Combining Gaia and GRAVITY: Characterising five new Directly Detected Substellar Companions

Authors: T. O. Winterhalder, S. Lacour, A. Mérand, A. -L. Maire, J. Kammerer, T. Stolker, N. Pourré, C. Babusiaux, R. Abuter, A. Amorim, R. Asensio-Torres, W. O. Balmer, M. Benisty, J. -P. Berger, H. Beust, S. Blunt, A. Boccaletti, M. Bonnefoy, H. Bonnet, M. S. Bordoni, G. Bourdarot, W. Brandner, F. Cantalloube, P. Caselli, B. Charnay , et al. (74 additional authors not shown)

Abstract: Precise mass constraints are vital for the characterisation of brown dwarfs and exoplanets. Here we present how the combination of data obtained by Gaia and GRAVITY can help enlarge the sample of substellar companions with measured dynamical masses. We show how the Non-Single-Star (NSS) two-body orbit catalogue contained in Gaia DR3 can be used to inform high-angular-resolution follow-up observati… ▽ More Precise mass constraints are vital for the characterisation of brown dwarfs and exoplanets. Here we present how the combination of data obtained by Gaia and GRAVITY can help enlarge the sample of substellar companions with measured dynamical masses. We show how the Non-Single-Star (NSS) two-body orbit catalogue contained in Gaia DR3 can be used to inform high-angular-resolution follow-up observations with GRAVITY. Applying the method presented in this work to eight Gaia candidate systems, we detect all eight predicted companions, seven of which were previously unknown and five are of a substellar nature. Among the sample is Gaia DR3 2728129004119806464 B, which - detected at an angular separation of (34.01 $\pm$ 0.15) mas from the host - is the closest substellar companion ever imaged. This translates to a semi-major axis of (0.938 $\pm$ 0.023) AU. WT 766 B, detected at a greater angular separation, was confirmed to be on an orbit exhibiting an even smaller semi-major axis of (0.676 $\pm$ 0.008) AU. The GRAVITY data were then used to break the host-companion mass degeneracy inherent to the Gaia NSS orbit solutions as well as to constrain the orbital solutions of the respective target systems. Knowledge of the companion masses enabled us to further characterise them in terms of their ages, effective temperatures, and radii via the application of evolutionary models. The inferred ages exhibit a distinct bias towards values younger than what is to be expected based on the literature. The results serve as an independent validation of the orbital solutions published in the NSS two-body orbit catalogue and show that the combination of astrometric survey missions and high-angular-resolution direct imaging holds great promise for efficiently increasing the sample of directly imaged companions in the future, especially in the light of Gaia's upcoming DR4 and the advent of GRAVITY+. △ Less

Submitted 24 June, 2024; v1 submitted 19 March, 2024; originally announced March 2024.

Comments: Accepted for publication in Astronomy and Astrophysics

Journal ref: A&A 688, A44 (2024)

FAUST XI: Enhancement of the complex organic material in the shocked matter surrounding the [BHB2007] 11 protobinary system

Authors: C. Vastel, T. Sakai, C. Ceccarelli, I. Jiménez-Serra, F. Alves, N. Balucani, E. Bianchi, M. Bouvier, P. Caselli, C. J. Chandler, S. Charnley, C. Codella, M. De Simone, F. Dulieu, L. Evans, F. Fontani, B. Lefloch, L. Loinard, F. Menard, L. Podio, G. Sabatini, N. Sakai, S. Yamamoto

Abstract: iCOMs are species commonly found in the interstellar medium. They are believed to be crucial seed species for the build-up of chemical complexity in star forming regions as well as our own Solar System. Thus, understanding how their abundances evolve during the star formation process and whether it enriches the emerging planetary system is of paramount importance. We use data from the ALMA Large P… ▽ More iCOMs are species commonly found in the interstellar medium. They are believed to be crucial seed species for the build-up of chemical complexity in star forming regions as well as our own Solar System. Thus, understanding how their abundances evolve during the star formation process and whether it enriches the emerging planetary system is of paramount importance. We use data from the ALMA Large Program FAUST to study the compact line emission towards the [BHB2007] 11 proto-binary system (sources A and B), where a complex structure of filaments connecting the two sources with a larger circumbinary disk has previously been detected. More than 45 CH3OCHO lines are clearly detected, as well as 8 CH3OCH3 transitions , 1 H2CCO transition and 4 t-HCOOH transitions. We compute the abundance ratios with respect to CH3OH for CH3OCHO, CH3OCH3, H2CCO, t-HCOOH (as well as an upper limit for CH3CHO) through a radiative transfer analysis. We also report the upper limits on the column densities of nitrogen bearing iCOMs, N(C2H5CN) and N(C2H3CN). The emission from the detected iCOMs and their precursors is compact and encompasses both protostars, which are separated by only 0.2" (~ 28 au). The integrated intensities tend to align with the Southern filament, revealed by the high spatial resolution observations of the dust emission at 1.3 mm. A PV and 2D analysis are performed on the strongest and uncontaminated CH3OCH3 transition and show three different spatial and velocity regions, two of them being close to 11B (Southern filament) and the third one near 11A. All our observations suggest that the detected methanol, as well as the other iCOMs, are generated by the shocked gas from the incoming filaments streaming towards [BHB2007] 11A and 11B, respectively, making this source one of the few where chemical enrichment of the gas caused by the streaming material is observed. △ Less

Submitted 12 March, 2024; originally announced March 2024.

Comments: Accepted in A&A

Fractionation in young cores: Direct determinations of nitrogen and carbon fractionation in HCN

Authors: S. S. Jensen, S. Spezzano, P. Caselli, O. Sipilä, E. Redaelli, K. Giers, J. Ferrer Asensio

Abstract: 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… ▽ More 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. △ Less

Submitted 7 March, 2024; originally announced March 2024.

Comments: Accepted for publication in A&A

The Reservoir of the Per-emb-2 Streamer

Authors: Kotomi Taniguchi, Jaime E Pineda, Paola Caselli, Tomomi Shimoikura, Rachel K. Friesen, Dominique M. Segura-Cox, Anika Schmiedeke

Abstract: Streamers bring gas from outer regions to protostellar systems and could change the chemical composition around protostars and protoplanetary disks. We have carried out mapping observations of carbon-chain species (HC$_3$N, HC$_5$N, CCH, and CCS) in the 3mm and 7mm bands toward the streamer flowing to the Class 0 young stellar object (YSO) Per-emb-2 with the Nobeyama 45m radio telescope. A region… ▽ More Streamers bring gas from outer regions to protostellar systems and could change the chemical composition around protostars and protoplanetary disks. We have carried out mapping observations of carbon-chain species (HC$_3$N, HC$_5$N, CCH, and CCS) in the 3mm and 7mm bands toward the streamer flowing to the Class 0 young stellar object (YSO) Per-emb-2 with the Nobeyama 45m radio telescope. A region with a diameter of $\sim0.04$ pc is located north with a distance of $\sim 20,500$ au from the YSO. The streamer connects to this north region which is the origin of the streamer. The reservoir has high density and low temperature ($n_{\rm {H}_2} \approx 1.9 \times 10^4$ cm$^{-3}$, $T_{\rm {kin}} = 10$ K), which are similar to those of early stage starless cores. By comparisons with the observed abundance ratios of CCS/HC$_3$N to the chemical simulations, the reservoir and streamer are found to be chemically young. The total mass available for the streamer is derived to be $24-34$ M$_{\odot}$. If all of the gas in the reservoir will accrete onto the Per-emb-2 protostellar system, the lifetime of the streamer has been estimated at ($1.1 - 3.2$)$\times10^{5}$ yr, suggesting that the mass accretion via the streamer would continue until the end of the Class I stage. △ Less

Submitted 7 March, 2024; v1 submitted 29 February, 2024; originally announced February 2024.

Comments: Accepted for The Astrophysical Journal (ApJ), 18 pages, 10 figures, 6 tables

Probing the physics of star formation (ProPStar): I. First resolved maps of the electron fraction and cosmic-ray ionization rate in NGC 1333

Authors: Jaime E. Pineda, Olli Sipilä, Dominique M. Segura-Cox, Maria Teresa Valdivia-Mena, Roberto Neri, Michael Kuffmeier, Alexei V. Ivlev, Stella S. R. Offner, Maria Jose Maureira, Paola Caselli, Silvia Spezzano, Nichol Cunningham, Anika Schmiedeke, Mike Chen

Abstract: Electron fraction and cosmic-ray ionization rates (CRIR) in star-forming regions are important quantities in astrochemical modeling and are critical to the degree of coupling between neutrals, ions, and electrons, which regulates the dynamics of the magnetic field. However, these are difficult quantities to estimate. We aim to derive the electron fraction and CRIR maps of an active star-forming re… ▽ More Electron fraction and cosmic-ray ionization rates (CRIR) in star-forming regions are important quantities in astrochemical modeling and are critical to the degree of coupling between neutrals, ions, and electrons, which regulates the dynamics of the magnetic field. However, these are difficult quantities to estimate. We aim to derive the electron fraction and CRIR maps of an active star-forming region. We combined observations of the nearby NGC 1333 star-forming region carried out with the NOEMA interferometer and IRAM 30-m single dish to generate high spatial dynamic range maps of different molecular transitions. We used the DCO$^+$ and H$^{13}$CO$^+$ ratio (in addition to complementary data) to estimate the electron fraction and produce cosmic-ray ionization rate maps. We derived the first large-area electron fraction and CRIR resolved maps in a star-forming region, with typical values of $10^{-6.5}$ and $10^{-16.5}$ s$^{-1}$, respectively. The maps present clear evidence of enhanced values around embedded young stellar objects (YSOs). This provides strong evidence for locally accelerated cosmic rays. We also found a strong enhancement toward the northwest region in the map that might be related either to an interaction with a bubble or to locally generated cosmic rays by YSOs. We used the typical electron fraction and derived a MHD turbulence dissipation scale of 0.054 pc, which could be tested with future observations. We found a higher cosmic-ray ionization rate compared to the canonical value for $N({\rm H_2})=10^{21}-10^{23}$ cm$^{-2}$ of $10^{-17}$ s$^{-1}$ in the region, and it is likely generated by the accreting YSOs. The high value of the electron fraction suggests that new disks will form from gas in the ideal-MHD limit. This indicates that local enhancements of $ζ({\rm H_2})$, due to YSOs, should be taken into account in the analysis of clustered star formation. △ Less

Submitted 25 February, 2024; originally announced February 2024.

Comments: 13 pages, 13 Figures. Accepted to A&A. Abridged abstract. GitHub repo with code available at https://github.com/jpinedaf/NGC1333_NOEMA_Ions

Journal ref: A&A 686, A162 (2024)

PRODIGE -- Planet-forming disks in Taurus with NOEMA. I. Overview and first results for 12CO, 13CO, and C18O

Authors: D. Semenov, Th. Henning, S. Guilloteau, G. Smirnov-Pinchukov, A. Dutrey, E. Chapillon, V. Pietu, R. Franceschi, K. Schwarz, S. van Terwisga, L. Bouscasse, P. Caselli, C. Ceccarelli, N. Cunningham, A. Fuente, C. Gieser, T. -H. Hsieh, A. Lopez-Sepulcre, D. M. Segura-Cox, J. E. Pineda, M. J. Maureira, Th. Moeller, M. Tafalla, M. T. Valdivia-Mena

Abstract: We are performing a line survey of 8 planet-forming Class II disks in Taurus with the IRAM NOrthern Extended Millimeter Array (NOEMA), as a part of the MPG-IRAM Observatory Program PRODIGE (PROtostars and DIsks: Global Evolution; PIs: P. Caselli and Th. Henning). Compact and extended disks around T Tauri stars CI, CY, DG, DL, DM, DN, IQ Tau, and UZ Tau E are observed in ~80 lines from >20 C-, O,-… ▽ More We are performing a line survey of 8 planet-forming Class II disks in Taurus with the IRAM NOrthern Extended Millimeter Array (NOEMA), as a part of the MPG-IRAM Observatory Program PRODIGE (PROtostars and DIsks: Global Evolution; PIs: P. Caselli and Th. Henning). Compact and extended disks around T Tauri stars CI, CY, DG, DL, DM, DN, IQ Tau, and UZ Tau E are observed in ~80 lines from >20 C-, O,- N-, and S-bearing species. The observations in four spectral settings at 210-280 GHz with $1σ$ rms sensitivity of ~ 8-12 mJy/beam at 0.9" and 0.3 km/s resolution will be completed in 2024. The uv-visibilities are fitted with the DiskFit model to obtain key stellar and disk properties. In this paper, the combined $^{12}$CO, $^{13}$CO and C$^{18}$O $J = 2-1$ data are presented. We find that the CO fluxes and disk masses inferred from dust continuum tentatively correlate with the CO emission sizes. We constrain dynamical stellar masses, geometries, temperatures, the CO column densities and gas masses for each disk. The best-fit temperatures at 100 au are ~ 17-37 K, and decrease radially with the power-law exponent q ~ 0.05-0.76. The inferred CO column densities decrease radially with the power-law exponent p ~ 0.2-3.1. The gas masses estimated from $^{13}$CO (2-1) are ~ $0.001-0.2 M_\textrm{Sun}$. The best-fit CO column densities point to severe CO freeze-out in the disks. The DL Tau disk is an outlier, and has either stronger CO depletion or lower gas mass than the rest of the sample. The CO isotopologue ratios are roughly consistent with the observed values in disks and the low-mass star-forming regions. △ Less

Submitted 27 February, 2024; v1 submitted 22 February, 2024; originally announced February 2024.

Comments: 47 pages, 10 figures, 6 tables, accepted for publication by A&A on 22/02/2024

Journal ref: A&A 685, A126 (2024)

Testing analytical methods to derive the cosmic-ray ionisation rate in cold regions via synthetic observations

Authors: E. Redaelli, S. Bovino, A. Lupi, T. Grassi, D. Gaete-Espinoza, G. Sabatini, P. Caselli

Abstract: 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 $ζ_2$ from molecular line observations. These have been derived from the chemical kinetics of the i… ▽ More 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 $ζ_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 $ζ_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 $ζ_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 $ζ_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 $ζ_2$ by at least two orders of magnitudes. We also consider a new derivation of this method, which, however, still leads to large overestimates. △ Less

Submitted 16 February, 2024; originally announced February 2024.

Comments: Accepted for publication in A&A (15th Feb.)

FAUST XII. Accretion streamers and jets in the VLA 1623--2417 protocluster

Authors: C. Codella, L. Podio, M. De Simone, C. Ceccarelli, S. Ohashi, C. J. Chandler, N. Sakai, J. E. Pineda, D. M. Segura-Cox, E. Bianchi, N. Cuello, A. López-Sepulcre, D. Fedele, P. Caselli, S. Charnley, D. Johnstone, Z. E. Zhang, M. J. Maureira, Y. Zhang, G. Sabatini, B. Svoboda, I. Jiménez-Serra, L. Loinard, S. Mercimek, N. Murillo , et al. (1 additional authors not shown)

Abstract: The ALMA interferometer has played a key role in revealing a new component of the Sun-like star forming process: the molecular streamers, i.e. structures up to thousands of au long funneling material non-axisymmetrically to disks. In the context of the FAUST ALMA LP, the archetypical VLA1623-2417 protostellar cluster has been imaged at 1.3 mm in the SO(5$_6$--4$_5$), SO(6$_6$--5$_5$), and SiO(5--4… ▽ More The ALMA interferometer has played a key role in revealing a new component of the Sun-like star forming process: the molecular streamers, i.e. structures up to thousands of au long funneling material non-axisymmetrically to disks. In the context of the FAUST ALMA LP, the archetypical VLA1623-2417 protostellar cluster has been imaged at 1.3 mm in the SO(5$_6$--4$_5$), SO(6$_6$--5$_5$), and SiO(5--4) line emission at the spatial resolution of 50 au. We detect extended SO emission, peaking towards the A and B protostars. Emission blue-shifted down to 6.6 km s$^{-1}$ reveals for the first time a long ($\sim$ 2000 au) accelerating streamer plausibly feeding the VLA1623 B protostar. Using SO, we derive for the first time an estimate of the excitation temperature of an accreting streamer: 33$\pm$9 K. The SO column density is $\sim$ 10$^{14}$ cm$^{-2}$, and the SO/H$_2$ abundance ratio is $\sim$ 10$^{-8}$. The total mass of the streamer is 3 $\times$ 10$^{-3}$ $Msun$, while its accretion rate is 3--5 $\times$ 10$^{-7}$ Msun yr$^{-1}$. This is close to the mass accretion rate of VLA1623 B, in the 0.6--3 $\times$ 10$^{-7}$ Msun yr$^{-1}$ range, showing the importance of the streamer in contributing to the mass of protostellar disks. The highest blue- and red-shifted SO velocities behave as the SiO(5--4) emission, the latter species detected for the first time in VLA1623-2417: the emission is compact (100-200 au), and associated only with the B protostar. The SO excitation temperature is $\sim$ 100 K, supporting the occurrence of shocks associated with the jet, traced by SiO. △ Less

Submitted 15 February, 2024; originally announced February 2024.

Comments: accepted by MNRAS

The GRAVITY young stellar object survey XIII. Tracing the time-variable asymmetric disk structure in the inner AU of the Herbig star HD98922

Authors: GRAVITY Collaboration, V. Ganci, L. Labadie, K. Perraut, A. Wojtczak, J. Kaufhold, M. Benisty, E. Alecian, G. Bourdarot, W. Brandner, A. Caratti o Garatti, C. Dougados, R. Garcia Lopez, J. Sanchez-Bermudez, A. Soulain, A. Amorim, J. -P. Berger, P. Caselli, Y. Clénet, A. Drescher, A. Eckart, F. Eisenhauer, M. Fabricius, H. Feuchtgruber, P. Garcia , et al. (30 additional authors not shown)

Abstract: Temporal variability in the photometric and spectroscopic properties of protoplanetary disks is common in YSO. However, evidence pointing toward changes in their morphology over short timescales has only been found for a few sources, mainly due to a lack of high cadence observations at mas resolution. We combine GRAVITY multi-epoch observations of HD98922 at mas resolution with PIONIER archival da… ▽ More Temporal variability in the photometric and spectroscopic properties of protoplanetary disks is common in YSO. However, evidence pointing toward changes in their morphology over short timescales has only been found for a few sources, mainly due to a lack of high cadence observations at mas resolution. We combine GRAVITY multi-epoch observations of HD98922 at mas resolution with PIONIER archival data covering a total time span of 11 years. We interpret the interferometric visibilities and spectral energy distribution with geometrical models and through radiative transfer techniques. We investigated high-spectral-resolution quantities to obtain information on the properties of the HI BrG-line-emitting region. The observations are best fitted by a model of a crescent-like asymmetric dust feature located at 1 au and accounting for 70% of the NIR emission. The feature has an almost constant magnitude and orbits the central star with a possible sub-Keplerian period of 12 months, although a 9 month period is another, albeit less probable, solution. The radiative transfer models show that the emission originates from a small amount of carbon-rich (25%) silicates, or quantum-heated particles located in a low-density region. Among different possible scenarios, we favor hydrodynamical instabilities in the inner disk that can create a large vortex. The high spectral resolution differential phases in the BrG-line show that the hot-gas component is offset from the star and in some cases is located between the star and the crescent feature. The scale of the emission does not favor magnetospheric accretion as a driving mechanism. The scenario of an asymmetric disk wind or a massive accreting substellar or planetary companion is discussed. With this unique observational data set for HD98922, we reveal morphological variability in the innermost 2 au of its disk region. △ Less

Submitted 31 January, 2024; originally announced January 2024.

Comments: 45 pages, 20 figures, accepted by and to be published in Astronomy & Astrophysics (A&A)

Massive clumps in W43-main: Structure formation in an extensively shocked molecular cloud

Authors: Yuxin Lin, Friedrich Wyrowski, Hauyu Baobab Liu, Yan Gong, Olli Sipilä, Andrés F. Izquierdo, Timea Csengeri, Adam Ginsburg, Guang-Xing Li, Silvia Spezzano, Jaime E. Pineda, Silvia Leurini, Paola Caselli, Karl M. Menten

Abstract: W43-main is a massive molecular complex located at the interaction of the Scutum arm and the Galactic bar undergoing starburst activities. We aim to investigate the gas dynamics, in particular, the prevailing shock signatures from the cloud to clump scale and assess the impact of shocks on the formation of dense gas and early-stage cores. We have carried out NOEMA and IRAM-30m observations at 3 mm… ▽ More W43-main is a massive molecular complex located at the interaction of the Scutum arm and the Galactic bar undergoing starburst activities. We aim to investigate the gas dynamics, in particular, the prevailing shock signatures from the cloud to clump scale and assess the impact of shocks on the formation of dense gas and early-stage cores. We have carried out NOEMA and IRAM-30m observations at 3 mm with an angular resolution of $\sim$0.1 pc towards five massive clumps in W43 main. We use CH$_{3}$CCH and H$_{2}$CS lines to trace the extended gas temperature and CH$_{3}$OH lines to probe the volume density of the dense gas ($\gtrsim$10$^{5}$ cm$^{-3}$). The emission of SiO (2-1) is extensive across the region ($\sim$4 pc) and is mostly contained within a low-velocity regime, hinting at a large-scale origin of the shocks. The position-velocity maps of multiple tracers show systematic spatio-kinematic offsets supporting the cloud-cloud collision/merging scenario. We identify an additional extended velocity component in CCH emission, which coincides with one of the velocity components of the larger scale $^{13}$CO (2-1) emission, likely representing an outer, less dense gas layer in the cloud merging process. We find that the V-shaped, asymmetric SiO wings are tightly correlated with localised gas density enhancements, which is direct evidence of dense gas formation and accumulation in shocks. We resolve two categories of NH$_{2}$D cores: ones exhibiting only subsonic to transonic velocity dispersion, and the others with an additional supersonic velocity dispersion. The centroid velocities of the latter cores are correlated with the shock front seen by SiO. The kinematics of the $\sim$0.1 pc NH$_{2}$D cores are heavily imprinted by shock activities, and may represent a population of early-stage cores forming around the shock interface. △ Less

Submitted 30 January, 2024; originally announced January 2024.

Comments: 35 pages, 34 figures (including appendices); Accepted in A&A

Accurate ab initio spectroscopic studies of promising interstellar ethanolamine iminic precursors

Authors: D. Alberton, N. Inostroza-Pino, Ryan C. Fortenberry, V. Lattanzi, C. Endres, J. Fuentealba Zamponi, P. Caselli

Abstract: Context. The detection of NH2CH2CH2OH (ethanolamine) in molecular cloud G+0.693-0.027 adds an additional player to the prebiotic molecules discovered so far in the interstellar medium (ISM). As this molecule might be formed through condensed-phase hydrogenation steps, detecting one or more of the molecules involved might help to elucidate the chemical pathway leading to its production. Aims. The c… ▽ More Context. The detection of NH2CH2CH2OH (ethanolamine) in molecular cloud G+0.693-0.027 adds an additional player to the prebiotic molecules discovered so far in the interstellar medium (ISM). As this molecule might be formed through condensed-phase hydrogenation steps, detecting one or more of the molecules involved might help to elucidate the chemical pathway leading to its production. Aims. The chemical path involves the formation of four chemical species. In this work, we study the energies of the isomers involved, indicate the best candidates for detection purposes, and provide the distortion constants of the most energetically favoured isomers undetected so far. Methods. We used highly accurate CCSD(T)-F12/cc-pCVTZ-F12 computations to predict the lowest energy isomers as well as their spectroscopic constants, taking corrections for core electron correlation and scalar relativity into account. Results. We studied 14 isomers. We find that the lowest energy isomer proposed in previous studies is not the actual minimum. We provide a set of rotational and distortion constants of the two new most stable isomers together with their fundamental vibrational frequencies in order to guide the search for these important astrochemical precursors of prebiotic molecules in the ISM. △ Less

Submitted 25 January, 2024; originally announced January 2024.

Polarized Light from Massive Protoclusters (POLIMAP). I. Dissecting the role of magnetic fields in the massive infrared dark cloud G28.37+0.07

Authors: C-Y Law, Jonathan C. Tan, Raphael Skalidis, Larry Morgan, Duo Xu, Felipe de Oliveira Alves, Ashley T. Barnes, Natalie Butterfield, Paola Caselli, Giuliana Cosentino, Francesco Fontani, Jonathan D. Henshaw, Izaskun Jimenez-Serra, Wanggi Lim

Abstract: Magnetic fields may play a crucial role in setting the initial conditions of massive star and star cluster formation. To investigate this, we report SOFIA-HAWC+ $214\:μ$m observations of polarized thermal dust emission and high-resolution GBT-Argus C$^{18}$O(1-0) observations toward the massive Infrared Dark Cloud (IRDC) G28.37+0.07. Considering the local dispersion of $B$-field orientations, we p… ▽ More Magnetic fields may play a crucial role in setting the initial conditions of massive star and star cluster formation. To investigate this, we report SOFIA-HAWC+ $214\:μ$m observations of polarized thermal dust emission and high-resolution GBT-Argus C$^{18}$O(1-0) observations toward the massive Infrared Dark Cloud (IRDC) G28.37+0.07. Considering the local dispersion of $B$-field orientations, we produce a map of $B$-field strength of the IRDC, which exhibits values between $\sim0.03 - 1\:$mG based on a refined Davis-Chandrasekhar-Fermi (r-DCF) method proposed by Skalidis \& Tassis. Comparing to a map of inferred density, the IRDC exhibits a $B-n$ relation with a power law index of $0.51\pm0.02$, which is consistent with a scenario of magnetically-regulated anisotropic collapse. Consideration of the mass-to-flux ratio map indicates that magnetic fields are dynamically important in most regions of the IRDC. A virial analysis of a sample of massive, dense cores in the IRDC, including evaluation of magnetic and kinetic internal and surface terms, indicates consistency with virial equilibrium, sub-Alfvénic conditions and a dominant role for $B-$fields in regulating collapse. A clear alignment of magnetic field morphology with direction of steepest column density gradient is also detected. However, there is no preferred orientation of protostellar outflow directions with the $B-$field. Overall, these results indicate that magnetic fields play a crucial role in regulating massive star and star cluster formation and so need to be accounted for in theoretical models of these processes. △ Less

Submitted 21 January, 2024; originally announced January 2024.

Comments: Submitted to ApJ, comments welcome

The GRAVITY young stellar object survey XII. The hot gas disk component in Herbig Ae/Be stars

Authors: GRAVITY Collaboration, R. Garcia Lopez, A. Natta, R. Fedriani, A. Caratti o Garatti, J. Sanchez-Bermudez, K. Perraut, C. Dougados, Y. -I. Bouarour, J. Bouvier, W. Brandner, P. Garcia, M. Koutoulaki, L. Labadie, H. Linz, E. Al'ecian, M. Benisty, J. -P. Berger, G. Bourdarot, P. Caselli, Y. Clenet, P. T. de Zeeuw, R. Davies, A. Eckart, F. Eisenhauer , et al. (24 additional authors not shown)

Abstract: The region of protoplanetary disks closest to a star (within 1-2\,au) is shaped by a number of different processes, from accretion of the disk material onto the central star to ejection in the form of winds and jets. Optical and near-IR emission lines are potentially good tracers of inner disk processes if very high spatial and/or spectral resolution are achieved. In this paper, we exploit the cap… ▽ More The region of protoplanetary disks closest to a star (within 1-2\,au) is shaped by a number of different processes, from accretion of the disk material onto the central star to ejection in the form of winds and jets. Optical and near-IR emission lines are potentially good tracers of inner disk processes if very high spatial and/or spectral resolution are achieved. In this paper, we exploit the capabilities of the VLTI-GRAVITY near-IR interferometer to determine the location and kinematics of the hydrogen emission line Bracket gamma. We present VLTI-GRAVITY observations of the Bracket gamma line for a sample of 26 stars of intermediate mass (HAEBE), the largest sample so far analysed with near-IR interferometry. The Bracket gamma line was detected in 17 objects. The emission is very compact (in most cases only marginally resolved), with a size of 10-30R* (1-5 mas). About half of the total flux comes from even smaller regions, which are unresolved in our data. For eight objects, it was possible to determine the position angle (PA) of the line-emitting region, which is generally in agreement with that of the inner-dusty disk emitting the K-band continuum. The position-velocity pattern of the Bracket gamma line-emitting region of the sampled objects is roughly consistent with Keplerian rotation. The exception is HD~45677, which shows more extended emission and more complex kinematics. The most likely scenario for the Bracket gamma origin is that the emission comes from an MHD wind launched very close to the central star, in a region well within the dust sublimation radius. An origin in the bound gas layer at the disk surface cannot be ruled out, while accreting matter provides only a minor fraction of the total flux. These results show the potential of near-IR spectro-interferometry to study line emission in young stellar objects. △ Less

Submitted 15 January, 2024; originally announced January 2024.

Comments: A&A accepted

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

Authors: I. V. Petrashkevich, A. F. Punanova, P. Caselli, O. Sipilä, J. E. Pineda, R. K. Friesen, M. G. Korotaeva, A. I. Vasyunin

Abstract: In this work, we study deuterium fractionation in four starless cores in the low-mass star-forming region L1688 in the Ophiuchus molecular cloud. We study how the deuterium fraction ($R_D$) changes with environment, compare deuteration of ions and neutrals, core centre and its envelope, and attempt to reproduce the observed results with a gas-grain chemical model. We chose high and low gas density… ▽ More In this work, we study deuterium fractionation in four starless cores in the low-mass star-forming region L1688 in the Ophiuchus molecular cloud. We study how the deuterium fraction ($R_D$) changes with environment, compare deuteration of ions and neutrals, core centre and its envelope, and attempt to reproduce the observed results with a gas-grain chemical model. We chose high and low gas density tracers to study both core centre and the envelope. With the IRAM 30m antenna, we mapped N$_2$H$^+$(1-0), N$_2$D$^+$(1-0), H$^{13}$CO$^+$ (1-0) and (2-1), DCO$^+$(2-1), and $p$-NH$_2$D(1$_{11}$-1$_{01}$) towards the chosen cores. The missing $p$-NH$_3$ and N$_2$H$^+$(1-0) data were taken from the literature. To measure the molecular hydrogen column density, dust and gas temperature within the cores, we used the Herschel/SPIRE dust continuum emission data, the GAS survey data (ammonia), and the COMPLETE survey data to estimate the upper limit on CO depletion. We present the deuterium fraction maps for three species towards four starless cores. Deuterium fraction of the core envelopes traced by DCO$^+$/H$^{13}$CO$^+$ is one order of magnitude lower ($\sim$0.08) than that of the core central parts traced by the nitrogen-bearing species ($\sim$0.5). Deuterium fraction increases with the gas density as indicated by high deuterium fraction of high gas density tracers and low deuterium fraction of lower gas density tracers and by the decrease of $R_D$ with core radii, consistent with the predictions of the chemical model. Our model results show a good agreement with observations for $R_D$(N$_2$D$^+$/N$_2$H$^+$) and R$_D$(DCO$^+$/HCO$^+$) and underestimate the $R_D$(NH$_2$D/NH$_3$). △ Less

Submitted 10 January, 2024; v1 submitted 9 January, 2024; originally announced January 2024.

Comments: 30 pages with 22 figures, accepted for publication in MNRAS

MSC Class: 70F15 (Primary) ACM Class: J.2

Infall of material onto the filaments in Barnard 5

Authors: Spandan Choudhury, Jaime E. Pineda, Paola Caselli, Michael Chun-Yuan Chen, Stella S. R. Offner, Maria Teresa Valdivia-Mena

Abstract: Aims. We aim to study the structure and kinematics of the two filaments inside the subsonic core Barnard 5 in Perseus using high-resolution ($\approx$ 2400 au) NH3 data and a multi-component fit analysis. Methods. We used observations of NH3 (1,1) and (2,2) inversion transitions using the Very Large Array (VLA) and the Green Bank Telescope (GBT). We smoothed the data to a beam of 8'' to reliably… ▽ More Aims. We aim to study the structure and kinematics of the two filaments inside the subsonic core Barnard 5 in Perseus using high-resolution ($\approx$ 2400 au) NH3 data and a multi-component fit analysis. Methods. We used observations of NH3 (1,1) and (2,2) inversion transitions using the Very Large Array (VLA) and the Green Bank Telescope (GBT). We smoothed the data to a beam of 8'' to reliably fit multiple velocity components towards the two filamentary structures identified in B5. Results. Along with the core and cloud components, which dominate the flux in the line of sight, we detected two components towards the two filaments showing signs of infall. We also detected two additional components that can possibly trace new material falling into the subsonic core of B5. Conclusions. Following comparison with previous simulations of filament formation scenarios in planar geometry, we conclude that either the formation of the B5 filaments is likely to be rather cylindrically symmetrical or the filaments are magnetically supported. We also estimate infall rates of $1.6\times10^{-4}\,M_\odot\,yr^{-1}$ and $1.8\times10^{-4}\,M_\odot\,yr^{-1}$ (upper limits) for the material being accreted onto the two filaments. At these rates, the filament masses can change significantly during the core lifetime. We also estimate an upper limit of $3.5\times10^{-5}\,M_\odot\,yr^{-1}$ for the rate of possible infall onto the core itself. Accretion of new material onto cores indicates the need for a significant update to current core evolution models, where cores are assumed to evolve in isolation. △ Less

Submitted 28 December, 2023; originally announced December 2023.

Comments: 15 pages, 13 figures, accepted for publication on Astronomy and Astrophysics. The scripts used for the analysis presented in this paper can be found at https://github.com/SpandanCh/Barnard5_filaments

VLTI/GRAVITY Provides Evidence the Young, Substellar Companion HD 136164 Ab formed like a "Failed Star"

Authors: William O. Balmer, L. Pueyo, S. Lacour, J. J. Wang, T. Stolker, J. Kammerer, N. Pourré, M. Nowak, E. Rickman, S. Blunt, A. Sivaramakrishnan, D. Sing, K. Wagner, G. -D. Marleau, A. -M. Lagrange, R. Abuter, A. Amorim, R. Asensio-Torres, J. -P. Berger, H. Beust, A. Boccaletti, A. Bohn, M. Bonnefoy, H. Bonnet, M. S. Bordoni , et al. (71 additional authors not shown)

Abstract: Young, low-mass Brown Dwarfs orbiting early-type stars, with low mass ratios ($q\lesssim0.01$), appear intrinsically rare and present a formation dilemma: could a handful of these objects be the highest mass outcomes of ``planetary" formation channels (bottom up within a protoplanetary disk), or are they more representative of the lowest mass ``failed binaries" (formed via disk fragmentation, or c… ▽ More Young, low-mass Brown Dwarfs orbiting early-type stars, with low mass ratios ($q\lesssim0.01$), appear intrinsically rare and present a formation dilemma: could a handful of these objects be the highest mass outcomes of ``planetary" formation channels (bottom up within a protoplanetary disk), or are they more representative of the lowest mass ``failed binaries" (formed via disk fragmentation, or core fragmentation)? Additionally, their orbits can yield model-independent dynamical masses, and when paired with wide wavelength coverage and accurate system age estimates, can constrain evolutionary models in a regime where the models have a wide dispersion depending on initial conditions. We present new interferometric observations of the $16\,\mathrm{Myr}$ substellar companion HD~136164~Ab (HIP~75056~Ab) with VLTI/GRAVITY and an updated orbit fit including proper motion measurements from the Hipparcos-Gaia Catalogue of Accelerations. We estimate a dynamical mass of $35\pm10\,\mathrm{M_J}$ ($q\sim0.02$), making HD~136164~Ab the youngest substellar companion with a dynamical mass estimate. The new mass and newly constrained orbital eccentricity ($e=0.44\pm0.03$) and separation ($22.5\pm1\,\mathrm{au}$) could indicate that the companion formed via the low-mass tail of the Initial Mass Function. Our atmospheric fit to the \texttt{SPHINX} M-dwarf model grid suggests a sub-solar C/O ratio of $0.45$, and $3\times$ solar metallicity, which could indicate formation in the circumstellar disk via disk fragmentation. Either way, the revised mass estimate likely excludes ``bottom-up" formation via core accretion in the circumstellar disk. HD~136164~Ab joins a select group of young substellar objects with dynamical mass estimates; epoch astrometry from future \textit{Gaia} data releases will constrain the dynamical mass of this crucial object further. △ Less

Submitted 13 December, 2023; originally announced December 2023.

Comments: Accepted for publication in The Astronomical Journal. 9 figures, 3 tables

Experimental study of the binding energy of NH3 on different types of ice and its impact on the snow line of NH3 and H2O

Authors: S. Kakkenpara Suresh, F. Dulieu, J. Vitorino, P. Caselli

Abstract: N-bearing molecules (like N2H+ or NH3) are excellent tracers of high-density, low-temperature regions like dense cloud cores and could shed light into snowlines in protoplanetary disks and the chemical evolution of comets. However, uncertainties exist about the grain surface chemistry of these molecules -- which could play an important role in their formation and evolution. This study explores exp… ▽ More N-bearing molecules (like N2H+ or NH3) are excellent tracers of high-density, low-temperature regions like dense cloud cores and could shed light into snowlines in protoplanetary disks and the chemical evolution of comets. However, uncertainties exist about the grain surface chemistry of these molecules -- which could play an important role in their formation and evolution. This study explores experimentally the behaviour of NH$_3$ on surfaces mimicking grains under interstellar conditions alongside other major interstellar ice components (ie. H$_2$O, CO, CO$_2$). We performed co-deposition experiments using the Ultra High Vacuum (UHV) setup VENUS (VErs des NoUvelles Syntheses) of NH$_3$ along with other adsorbates (here, H$_2$O, $^{13}$CO and CO$_2$) and performed Temperature Programmed Desorption (TPD) and Temperature Programmed-During Exposure Desorption (TP-DED) experiments. We obtained binding Energy (BE) distribution of NH$_3$ on Crystalline Ice(CI) and compact-Amorphous Solid Water (c-ASW) by analyses of the TPD profiles of NH3 on the substrates. We observe a significant delay in the desorption and a decrease in the desorption rate of NH$_3$ when H$_2$O is introduced into the co-deposited mixture of NH$_3$-$^{13}$Co or NH$_3$-CO$_2$, absent without H$_2$O. Secondly, H$_2$O traps nearly 5-9 per cent of the co-deposited NH3, released during water's amorphous-to-crystalline phase change. Thirdly, for CI, we obtained a BE distribution between 3780K-4080K, and c-ASW between 3780K-5280K -- using a pre-exponential factor A = 1.94$\times 10^{15}$/s. We conclude that NH$_3$ behaviour is significantly influenced by the presence of H$_2$O due to the formation of hydrogen bonds, in line with quantum calculations. This interaction preserves NH$_3$ on grain surfaces to higher temperatures making it available to the central protostar in protoplanetary disks. It also explains why NH$_3$ freeze out in pre-stellar cores is efficient. △ Less

Submitted 30 November, 2023; originally announced November 2023.

A Collection of German Science Interests in the Next Generation Very Large Array

Authors: M. Kadler, D. A. Riechers, J. Agarwal, A. -K. Baczko, H. Beuther, F. Bigiel, T. Birnstiel, B. Boccardi, D. J. Bomans, L. Boogaard, T. T. Braun, S. Britzen, M. Brüggen, A. Brunthaler, P. Caselli, D. Elsässer, S. von Fellenberg, M. Flock, C. M. Fromm, L. Fuhrmann, P. Hartogh, M. Hoeft, R. P. Keenan, Y. Kovalev, K. Kreckel , et al. (66 additional authors not shown)

Abstract: The Next Generation Very Large Array (ngVLA) is a planned radio interferometer providing unprecedented sensitivity at wavelengths between 21 cm and 3 mm. Its 263 antenna element array will be spatially distributed across North America to enable both superb low surface brightness recovery and sub-milliarcsecond angular resolution imaging. The project was developed by the international astronomy com… ▽ More The Next Generation Very Large Array (ngVLA) is a planned radio interferometer providing unprecedented sensitivity at wavelengths between 21 cm and 3 mm. Its 263 antenna element array will be spatially distributed across North America to enable both superb low surface brightness recovery and sub-milliarcsecond angular resolution imaging. The project was developed by the international astronomy community under the lead of the National Radio Astronomy Observatory (NRAO), and is anticipated to be built between 2027 and 2037. Two workshops have been held in 2022 and 2023 with the goal to discuss and consolidate the scientific interests in the ngVLA within the German astronomical community. This community paper constitutes a collection of 48 science ideas which the German community aims to pursue with the ngVLA in the 2030s. This is not a complete list and the ideas are not developed at the level of a "Science Book", such that the present document is mainly meant provide a basis for further discussion within the community. As such, additional contributions are welcome, and will be considered for inclusion in future revisions. △ Less

Submitted 18 June, 2024; v1 submitted 16 November, 2023; originally announced November 2023.

Comments: Version 2.0 (status June 18, 2024): 169 pages, comments and future contributions welcome [v2.0: 7 new science cases added, some minor revisions to other chapters]

Rotational spectroscopic characterisation of the [D2,C,S] system: an update from the laboratory and theory

Authors: Natalia Inostroza-Pino, Valerio Lattanzi, C. Zachary Palmer, Ryan C. Fortenberry, Diego Mardones, Paola Caselli, Oko E. Godwin, Timothy J. Lee

Abstract: The synergy between high-resolution rotational spectroscopy and quantum-chemical calculations is essential for exploring future detection of molecules, especially when spectroscopy parameters are not available yet. By using highly correlated ab initio quartic force fields (QFFs) from explicitly correlated coupled cluster theory, a complete set of rotational constants and centrifugal distortion con… ▽ More The synergy between high-resolution rotational spectroscopy and quantum-chemical calculations is essential for exploring future detection of molecules, especially when spectroscopy parameters are not available yet. By using highly correlated ab initio quartic force fields (QFFs) from explicitly correlated coupled cluster theory, a complete set of rotational constants and centrifugal distortion constants for D$_2$CS and cis/trans-DCSD isomers have been produced. Comparing our new ab initio results for D$_2$CS with new rotational spectroscopy laboratory data for the same species, the accuracy of the computed B and C rotational constants is within 0.1% while the A constant is only slightly higher. Additionally, quantum chemical vibrational frequencies are also provided, and these spectral reference data and new experimental rotational lines will provide additional references for potential observation of these deuterated sulfur species with either ground-based radio telescopes or space-based infrared observatories. △ Less

Submitted 15 November, 2023; originally announced November 2023.

Comments: 15 pages, 2 figures and 1 appendix table; published in Molecular Physics (Tim Lee Memorial Issue)

Nuclear spin ratios of deuterated ammonia in prestellar cores. LAsMA observations of H-MM1 and Oph D

Authors: Jorma Harju, Jaime E. Pineda, Olli Sipilä, Paola Caselli, Arnaud Belloche, Friedrich Wyrowski, Wiebke Riedel, Elena Redaelli, Anton I. Vasyunin

Abstract: We determine the ortho/para ratios of NH2D and NHD2 in two dense, starless cores, where their formation is supposed to be dominated by gas-phase reactions, which, in turn, is predicted to result in deviations from the statistical spin ratios. The Large APEX sub-Millimeter Array (LAsMA) multibeam receiver of the Atacama Pathfinder EXperiment (APEX) telescope was used to observe the prestellar cores… ▽ More We determine the ortho/para ratios of NH2D and NHD2 in two dense, starless cores, where their formation is supposed to be dominated by gas-phase reactions, which, in turn, is predicted to result in deviations from the statistical spin ratios. The Large APEX sub-Millimeter Array (LAsMA) multibeam receiver of the Atacama Pathfinder EXperiment (APEX) telescope was used to observe the prestellar cores H-MM1 and Oph D in Ophiuchus in the ground-state lines of ortho and para NH2D and NHD2. The fractional abundances of these molecules were derived employing 3D radiative transfer modelling, using different assumptions about the abundance profiles as functions of density. We also ran gas-grain chemistry models with different scenarios concerning proton or deuteron exchanges and chemical desorption from grains to find out if one of these models can reproduce the observed spin ratios. The observationally deduced ortho/para ratios of NH2D and NHD2 are in both cores within 10% of their statistical values 3 and 2, respectively, and taking 3-sigma limits, deviations from these of about 20% are allowed. Of the chemistry models tested here, the model that assumes proton hop (as opposed to full scrambling) in reactions contributing to ammonia formation, and a constant efficiency of chemical desorption, comes nearest to the observed abundances and spin ratios. The nuclear spin ratios derived here are in contrast with spin-state chemistry models that assume full scrambling in proton donation and hydrogen abstraction reactions leading to deuterated ammonia. The efficiency of chemical desorption influences strongly the predicted abundances of NH3, NH2D, and NHD2, but has a lesser effect on their ortho/para ratios. For these the proton exchange scenario in the gas is decisive. We suggest that this is because of rapid re-processing of ammonia and related cations by gas-phase ion-molecule reactions. △ Less

Submitted 14 November, 2023; originally announced November 2023.

Comments: accepted for publication in Astronomy & Astrophysics

Exploring the dust grain size and polarization mechanism in the hot and massive Class 0 disk IRAS 16293-2422 B

Authors: Joaquin Zamponi, María José Maureira, Hauyu Baobab Liu, Bo Zhao, Dominique Segura-Cox, Chia-Lin Ko, Paola Caselli

Abstract: Multiwavelength dust continuum and polarization observations arising from self-scattering have been used to investigate grain sizes in young disks. However, the polarization by self-scattering is low in face-on optically thick disks and puts some of the size constraints from polarization on hold, particularly for the younger and more massive disks. The 1.3 mm emission detected toward the hot (… ▽ More Multiwavelength dust continuum and polarization observations arising from self-scattering have been used to investigate grain sizes in young disks. However, the polarization by self-scattering is low in face-on optically thick disks and puts some of the size constraints from polarization on hold, particularly for the younger and more massive disks. The 1.3 mm emission detected toward the hot ($\gtrsim$400 K) Class 0 disk IRAS 16293-2422 B has been attributed to self-scattering, predicting grain sizes between 200-2000 $μ$m. We investigate the effects of grain size in the resultant flux and polarization fractions from self-scattering using a hot and massive Class 0 disk model and compare with observations. We compared new and archival high-resolution observations between 1.3 and 18 mm to a set of synthetic models. We have developed a new public tool to automate this process called Synthesizer. This is an easy-to-use program to generate synthetic observations from numerical simulations. Optical depths are in the range of 130 to 2 from 1.3 to 18 mm, respectively. Predictions from significant grain growth populations, including millimetric grains are comparable to the observations at all wavelengths. The polarization fraction produced by self-scattering reaches a maximum of $\sim$0.1% at 1.3 mm for a maximum grain size of 100 $μ$m, being an order of magnitude lower than that observed with ALMA. From the comparison of Stokes I fluxes, we conclude that significant grain growth could be present in the young Class 0 disk IRAS 16293 B, particularly in the inner hot region ($<10$ au, $T>$ 300 K) where refractory organics evaporate. The polarization produced by self-scattering in our model is not high enough to explain the observations at 1.3 and 7 mm, and effects like dichroic extinction or polarization reversal of elongated aligned grains remain other possible but untested scenarios. △ Less

Submitted 4 November, 2023; originally announced November 2023.

Comments: Accepted for publication in A&A

The role of turbulence in high-mass star formation: Subsonic and transonic turbulence are ubiquitously found at early stages

Authors: Chao Wang, Ke Wang, Feng-Wei Xu, Patricio Sanhueza, Hauyu Baobab Liu, Qizhou Zhang, Xing Lu, F. Fontani, Paola Caselli, Gemma Busquet, Jonathan C. Tan, Di Li, J. M. Jackson, Thushara Pillai, Paul T. P. Ho, Andrés E. Guzmán, Nannan Yue

Abstract: Context. Traditionally, supersonic turbulence is considered to be one of the most likely mechanisms to slow down the gravitational collapse in dense clumps, thereby enabling the formation of massive stars. However, several recent studies have raised differing points of view based on observations carried out with sufficiently high spatial and spectral resolution. These studies call for a re-evaluat… ▽ More Context. Traditionally, supersonic turbulence is considered to be one of the most likely mechanisms to slow down the gravitational collapse in dense clumps, thereby enabling the formation of massive stars. However, several recent studies have raised differing points of view based on observations carried out with sufficiently high spatial and spectral resolution. These studies call for a re-evaluation of the role turbulence plays in massive star-forming regions. Aims. Our aim is to study the gas properties, especially the turbulence, in a sample of massive star-forming regions with sufficient spatial and spectral resolution, which can both resolve the core fragmentation and the thermal line width. Methods. We observed NH3 metastable lines with the Very Large Array (VLA) to assess the intrinsic turbulence. Results. Analysis of the turbulence distribution histogram for 32 identified NH3 cores reveals the presence of three distinct components. Furthermore, our results suggest that (1) sub- and transonic turbulence is a prevalent (21 of 32) feature of massive star-forming regions and those cold regions are at early evolutionary stage. This investigation indicates that turbulence alone is insufficient to provide the necessary internal pressure required for massive star formation, necessitating further exploration of alternative candidates; and (2) studies of seven multi-core systems indicate that the cores within each system mainly share similar gas properties and masses. However, two of the systems are characterized by the presence of exceptionally cold and dense cores that are situated at the spatial center of each system. Our findings support the hub-filament model as an explanation for this observed distribution △ Less

Submitted 7 February, 2024; v1 submitted 27 October, 2023; originally announced October 2023.

Comments: 34 pages, 15 figures, 4 tables. Accepted for publication on A&A

The evolution of sulphur-bearing molecules in high-mass star-forming cores

Authors: F. Fontani, E. Roueff, L. Colzi, P. Caselli

Abstract: To understand the chemistry of sulphur (S) in the interstellar medium, models need to be tested by observations of S-bearing molecules in different physical conditions. We analyse observations obtained with the IRAM 30m telescope towards 15 well-known cores classified in the three main evolutionary stages of the high-mass star-formation process: high-mass starless cores,high-mass protostellar obje… ▽ More To understand the chemistry of sulphur (S) in the interstellar medium, models need to be tested by observations of S-bearing molecules in different physical conditions. We analyse observations obtained with the IRAM 30m telescope towards 15 well-known cores classified in the three main evolutionary stages of the high-mass star-formation process: high-mass starless cores,high-mass protostellar objects, and ultracompact HII regions. We detected rotational lines of SO, SO+, NS, C34S, 13CS, SO2, CCS, H2S, HCS+, OCS, H2CS, and CCCS. We also analyse for the first time lines of the NO molecule to complement the analysis. From a local thermodynamic equilibrium approach, we have derived column densities of each species and excitation temperatures. Based on a statistical analysis on the line widths and the excitation temperatures, we find that: NS, C34S, 13CS, CCS, and HCS+ trace cold, quiescent, and likely extended material; OCS, and SO2 trace warmer, more turbulent, and likely denser and more compact material; SO and perhaps SO+ trace both quiescent and turbulent material depending on the target. The abundances of SO, SO2, and H2S show the strongest positive correlations with the kinetic temperature, believed to be an evolutionary indicator. Moreover, the sum of all molecular abundances show an enhancement of gaseous S from the less evolved to the more evolved stages. These trends could be due to the increasing amount of S sputtered from dust grains owing to the increasing protostellar activity with evolution. The average abundances in each evolutionary group increase especially in the oxygen-bearing molecules, perhaps due to the increasing abundance of atomic oxygen with evolution owing to photodissociation of water in gas phase. Our observational work represents a test-bed for theoretical studies aimed at modelling the chemistry of sulphur during the evolution of high-mass star-forming cores △ Less

Submitted 16 October, 2023; originally announced October 2023.

Comments: 16 pages. 3 appendices. Accepted for publication in A&A