Showing 1–50 of 51 results for author: Drozdovskaya, M

Rotational spectroscopy of CH$_3$OD with a reanalysis of CH$_3$OD toward IRAS 16293$-$2422

Authors: V. V. Ilyushin, H. S. P. Müller, M. N. Drozdovskaya, J. K. Jørgensen, S. Bauerecker, C. Maul, R. Porohovoi, E. A. Alekseev, O. Dorovskaya, O. Zakharenko, F. Lewen, S. Schlemmer, R. M. Lees

Abstract: We have started a measurement campaign of numerous methanol isotopologs in low-lying torsional states in order to provide extensive line lists for radio astronomical observations from an adequate spectroscopic model and to investigate how the intricate vibration-torsion-rotation interactions manifest themselves in the spectra of different isotopic species. After CD$_3$OH and CD$_3$OD, we turn our… ▽ More We have started a measurement campaign of numerous methanol isotopologs in low-lying torsional states in order to provide extensive line lists for radio astronomical observations from an adequate spectroscopic model and to investigate how the intricate vibration-torsion-rotation interactions manifest themselves in the spectra of different isotopic species. After CD$_3$OH and CD$_3$OD, we turn our focus to CH$_3$OD, which is an important species for studying deuteration in prestellar cores and envelopes that enshroud protostars. Notably, deuteration is frequently viewed as a diagnostic tool for star formation. The measurements used in this study were obtained in two spectroscopic laboratories and cover large fractions of the 34 GHz--1.35 THz range. As done in previous studies, we employed a torsion-rotation Hamiltonian model for our analysis that is based on the rho-axis method. The resulting model describes the ground and first excited torsional states of CH$_3$OD well up to quantum numbers $J \leqslant 51$ and $K_a \leqslant 18$. We derived a line list for radio astronomical observations from this model that is accurate up to at least 1.35~THz and should be sufficient for all types of radio astronomical searches for this methanol isotopolog in these two lowest torsional states. This line list was applied to a reinvestigation of CH$_3$OD in data from the Protostellar Interferometric Line Survey of IRAS 16293--2422 obtained with the Atacama Large Millimeter/submillimeter Array. The new accurately determined value for the column density of CH$_3$OD implies that the deuteration in methanol differs in its two functional groups by a factor of $\sim$7.5. △ Less

Submitted 12 June, 2024; originally announced June 2024.

Comments: 12 pages including figures, tables, references and appendix. Astron. Astrophys., in press

JWST MIRI MRS Images Disk Winds, Water, and CO in an Edge-On Protoplanetary Disk

Authors: Nicole Arulanantham, M. K. McClure, Klaus Pontoppidan, Tracy L. Beck, J. A. Sturm, D. Harsono, A. C. A. Boogert, M. Cordiner, E. Dartois, M. N. Drozdovskaya, C. Espaillat, G. J. Melnick, J. A. Noble, M. E. Palumbo, Y. J. Pendleton, H. Terada, E. F. van Dishoeck

Abstract: We present JWST MIRI MRS observations of the edge-on protoplanetary disk around the young sub-solar mass star Tau 042021, acquired as part of the Cycle 1 GO program "Mapping Inclined Disk Astrochemical Signatures (MIDAS)." These data resolve the mid-IR spatial distributions of H$_2$, revealing X-shaped emission extending to ~200 au above the disk midplane with a semi-opening angle of $35 \pm 5$ de… ▽ More We present JWST MIRI MRS observations of the edge-on protoplanetary disk around the young sub-solar mass star Tau 042021, acquired as part of the Cycle 1 GO program "Mapping Inclined Disk Astrochemical Signatures (MIDAS)." These data resolve the mid-IR spatial distributions of H$_2$, revealing X-shaped emission extending to ~200 au above the disk midplane with a semi-opening angle of $35 \pm 5$ degrees. We do not velocity-resolve the gas in the spectral images, but the measured semi-opening angle of the H$_2$ is consistent with an MHD wind origin. A collimated, bipolar jet is seen in forbidden emission lines from [Ne II], [Ne III], [Ni II], [Fe II], [Ar II], and [S III]. Extended H$_2$O and CO emission lines are also detected, reaching diameters between ~90 and 190 au, respectively. Hot molecular emission is not expected at such radii, and we interpret its extended spatial distribution as scattering of inner disk molecular emission by dust grains in the outer disk surface. H I recombination lines, characteristic of inner disk accretion shocks, are similarly extended, and are likely also scattered light from the innermost star-disk interface. Finally, we detect extended PAH emission at 11.3 microns co-spatial with the scattered light continuum, making this the first low-mass T Tauri star around which extended PAHs have been confirmed, to our knowledge. MIRI MRS line images of edge-on disks provide an unprecedented window into the outflow, accretion, and scattering processes within protoplanetary disks, allowing us to constrain the disk lifetimes and accretion and mass loss mechanisms. △ Less

Submitted 20 March, 2024; v1 submitted 19 February, 2024; originally announced February 2024.

Comments: Accepted for publication in ApJL on March 13th, 2024

A deep search for large complex organic species toward IRAS16293-2422 B at 3 mm with ALMA

Authors: P. Nazari, J. S. Y. Cheung, J. Ferrer Asensio, N. M. Murillo, E. F. van Dishoeck, J. K. Jørgensen, T. L. Bourke, K. -J. Chuang, M. N. Drozdovskaya, G. Fedoseev, R. T. Garrod, S. Ioppolo, H. Linnartz, B. A. McGuire, H. S. P. Müller, D. Qasim, S. F. Wampfler

Abstract: Complex organic molecules (COMs) have been detected ubiquitously in protostellar systems. However, at shorter wavelengths (~0.8mm) it is more difficult to detect larger molecules than at longer wavelengths (~3mm) because of the increase of millimeter dust opacity, line confusion, and unfavorable partition function. We aim to search for large molecules (>8 atoms) in the ALMA Band 3 spectrum of IRAS… ▽ More Complex organic molecules (COMs) have been detected ubiquitously in protostellar systems. However, at shorter wavelengths (~0.8mm) it is more difficult to detect larger molecules than at longer wavelengths (~3mm) because of the increase of millimeter dust opacity, line confusion, and unfavorable partition function. We aim to search for large molecules (>8 atoms) in the ALMA Band 3 spectrum of IRAS 16293-2422 B. We search for more than 70 molecules and identify as many lines as possible in the spectrum. The spectral settings were set to specifically target three-carbon species such as propanol and glycerol. We identify lines of 31 molecules including many oxygen-bearing COMs such as CH3OH and c-C2H4O and a few nitrogen- and sulfur-bearing ones such as HOCH2CN and CH3SH. The largest detected molecules are gGg-(CH2OH)2 and CH3COCH3. We do not detect glycerol or propanol but provide upper limits for them which are in line with previous laboratory and observational studies. The line density in Band 3 is only ~2.5 times lower in frequency space than in Band 7. From the detected lines in Band 3 at a $\gtrsim 6σ$ level, ~25-30% of them could not be identified indicating the need for more laboratory data of rotational spectra. We find similar column densities and column density ratios of COMs (within a factor ~2) between Band 3 and Band 7. The effect of dust optical depth for IRAS 16293-2422 B at an off-source location on column densities and column density ratios is minimal. Moreover, for warm protostars, long wavelength spectra are not only crowded, but also take longer integration times to reach the same sensitivity limit. The 3mm search has not yet resulted in detection of larger and more complex molecules in warm sources. A full deep ALMA Band 2-3 (i.e., 3-4 mm) survey is needed to assess whether low frequency data have the potential to reveal more complex molecules in warm sources. △ Less

Submitted 9 January, 2024; originally announced January 2024.

Comments: Accepted for publication in A&A

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

Spatially resolving the volatile sulfur abundance in the HD 100546 protoplanetary disk

Authors: Luke Keyte, Mihkel Kama, Ko-Ju Chuang, L. Ilsedore Cleeves, Maria N. Drozdovskaya, Kenji Furuya, Jonathan Rawlings, Oliver Shorttle

Abstract: Volatile elements play a crucial role in the formation of planetary systems. Their abundance and distribution in protoplanetary disks provide vital insights into the connection between formation processes and the atmospheric composition of individual planets. Sulfur, being one of the most abundant elements in planet-forming environments, is of great significance, and now observable in exoplanets w… ▽ More Volatile elements play a crucial role in the formation of planetary systems. Their abundance and distribution in protoplanetary disks provide vital insights into the connection between formation processes and the atmospheric composition of individual planets. Sulfur, being one of the most abundant elements in planet-forming environments, is of great significance, and now observable in exoplanets with JWST. However, planetary formation models currently lack vital knowledge regarding sulfur chemistry in protoplanetary disks. Developing a deeper understanding of the major volatile sulfur carriers in disks is essential to building models that can meaningfully predict planetary atmospheric composition, and reconstruct planetary formation pathways. In this work, we combine archival observations with new data from ALMA and APEX, covering a range of sulfur-bearing species/isotopologs. We interpret this data using the DALI thermo-chemical code, for which our model is highly refined and disk-specific. We find that volatile sulfur is heavily depleted from the cosmic value by a factor of 1000, with a disk-averaged abundance of S/H = 1e-8. We show that the gas-phase sulfur abundance varies radially by 3 orders of magnitude, with the highest abundances inside the inner dust ring and coincident with the outer dust ring at 150 to 230 au. Extracting chemical abundances from our models, we find OCS, H2CS, and CS to be the dominant molecular carriers in the gas phase. We also infer the presence of a substantial OCS ice reservoir. We relate our results to the potential atmospheric composition of planets in HD 100546, and the wider exoplanet population △ Less

Submitted 21 December, 2023; originally announced December 2023.

Comments: Accepted in MNRAS (20/12/2023)

A JWST inventory of protoplanetary disk ices: The edge-on protoplanetary disk HH 48 NE, seen with the Ice Age ERS program

Authors: J. A. Sturm, M. K. McClure, T. L. Beck, D. Harsono, J. B. Bergner, E. Dartois, A. C. A. Boogert, J. E. Chiar, M. A. Cordiner, M. N. Drozdovskaya, S. Ioppolo, C. J. Law, H. Linnartz, D. C. Lis, G. J. Melnick, B. A. McGuire, J. A. Noble, K. I. Öberg, M. E. Palumbo, Y. J. Pendleton, G. Perotti, K. M. Pontoppidan, D. Qasim, W. R. M. Rocha, H. Terada , et al. (2 additional authors not shown)

Abstract: Ices are the main carriers of volatiles in protoplanetary disks and are crucial to our understanding of the chemistry that ultimately sets the organic composition of planets. The ERS program Ice Age on the JWST follows the ice evolution through all stages of star and planet formation. JWST/NIRSpec observations of the edge-on Class II protoplanetary disk HH~48~NE reveal spatially resolved absorptio… ▽ More Ices are the main carriers of volatiles in protoplanetary disks and are crucial to our understanding of the chemistry that ultimately sets the organic composition of planets. The ERS program Ice Age on the JWST follows the ice evolution through all stages of star and planet formation. JWST/NIRSpec observations of the edge-on Class II protoplanetary disk HH~48~NE reveal spatially resolved absorption features of the major ice components H$_2$O, CO$_2$, CO, and multiple weaker signatures from less abundant ices NH$_3$, OCN$^-$, and OCS. Isotopologue $^{13}$CO$_2$ ice has been detected for the first time in a protoplanetary disk. Since multiple complex light paths contribute to the observed flux, the ice absorption features are filled in by ice-free scattered light. The $^{12}$CO$_2$/$^{13}$CO$_2$ ratio of 14 implies that the $^{12}$CO$_2$ feature is saturated, without the flux approaching 0, indicative of a very high CO$_2$ column density on the line of sight, and a corresponding abundance with respect to hydrogen that is higher than ISM values by a factor of at least a few. Observations of rare isotopologues are crucial, as we show that the $^{13}$CO$_2$ observation allows us to determine the column density of CO$_2$ to be at an order of magnitude higher than the lower limit directly inferred from the observed optical depth. Radial variations in ice abundance, e.g., snowlines, are significantly modified since all observed photons have passed through the full radial extent of the disk. CO ice is observed at perplexing heights in the disk, extending to the top of the CO-emitting gas layer. We argue that the most likely interpretation is that we observe some CO ice at high temperatures, trapped in less volatile ices like H$_2$O and CO$_2$. Future radiative transfer models will be required to constrain the implications on our current understanding of disk physics and chemistry. △ Less

Submitted 14 September, 2023; originally announced September 2023.

Comments: 16 pages, 8 figures, accepted for publication in A&A

Evolution of Chemistry in the envelope of Hot Corinos (ECHOS). I. Extremely young sulphur chemistry in the isolated Class 0 object B335

Authors: G. Esplugues, M. Rodríguez-Baras, D. San Andrés, D. Navarro-Almaida, A. Fuente, P. Rivière-Marichalar, Á. Sánchez-Monge, M. N. Drozdovskaya, S. Spezzano, P. Caselli

Abstract: Within the project Evolution of Chemistry in the envelope of HOt corinoS (ECHOS), we present a study of sulphur chemistry in the envelope of the Class 0 source B335 through observations in the spectral range 7, 3, and 2 mm. We have modelled observations assuming LTE and LVG approximation. We have also used the code Nautilus to study the time evolution of sulphur species. We have detected 20 sulphu… ▽ More Within the project Evolution of Chemistry in the envelope of HOt corinoS (ECHOS), we present a study of sulphur chemistry in the envelope of the Class 0 source B335 through observations in the spectral range 7, 3, and 2 mm. We have modelled observations assuming LTE and LVG approximation. We have also used the code Nautilus to study the time evolution of sulphur species. We have detected 20 sulphur species with a total gas-phase S abundance similar to that found in the envelopes of other Class 0 objects, but with significant differences in the abundances between sulphur carbon chains and sulphur molecules containing oxygen and nitrogen. Our results highlight the nature of B335 as a source especially rich in sulphur carbon chains unlike other Class 0 sources. The low presence or absence of some molecules, such as SO and SO+, suggests a chemistry not particularly influenced by shocks. We, however, detect a large presence of HCS+ that, together with the low rotational temperatures obtained for all the S species (<15 K), reveals the moderate or low density of the envelope of B335. We also find that observations are better reproduced by models with a sulphur depletion factor of 10 with respect to the sulphur cosmic elemental abundance. The comparison between our model and observational results for B335 reveals an age of 10$^4$$<$t$<$10$^5$ yr, which highlights the particularly early evolutionary stage of this source. B335 presents a different chemistry compared to other young protostars that have formed in dense molecular clouds, which could be the result of accretion of surrounding material from the diffuse cloud onto the protostellar envelope of B335. In addition, the analysis of the SO2/C2S, SO/CS, and HCS+/CS ratios within a sample of prestellar cores and Class 0 objects show that they could be used as good chemical evolutionary indicators of the prestellar to protostellar transition. △ Less

Submitted 4 September, 2023; originally announced September 2023.

Journal ref: A&A 678, A199 (2023)

CoCCoA: Complex Chemistry in hot Cores with ALMA. Selected oxygen-bearing species

Authors: Y. Chen, M. L. van Gelder, P. Nazari, C. L. Brogan, E. F. van Dishoeck, H. Linnartz, J. K. Jørgensen, T. R. Hunter, O. H. Wilkins, G. A. Blake, P. Caselli, K. -J. Chuang, C. Codella, I. Cooke, M. N. Drozdovskaya, R. T. Garrod, S. Ioppolo, M. Jin, B. M. Kulterer, N. F. W. Ligterink, A. Lipnicky, R. Loomis, M. G. Rachid, S. Spezzano, B. A. McGuire

Abstract: Complex organic molecules (COMs) have been observed to be abundant in the gas phase toward protostars. Deep line surveys have been carried out only for a limited number of well-known high-mass star forming regions using the Atacama Large Millimeter/submillimeter Array (ALMA), which has unprecedented resolution and sensitivity. Statistical studies on oxygen-bearing COMs (O-COMs) in high-mass protos… ▽ More Complex organic molecules (COMs) have been observed to be abundant in the gas phase toward protostars. Deep line surveys have been carried out only for a limited number of well-known high-mass star forming regions using the Atacama Large Millimeter/submillimeter Array (ALMA), which has unprecedented resolution and sensitivity. Statistical studies on oxygen-bearing COMs (O-COMs) in high-mass protostars using ALMA are still lacking. With the recent CoCCoA survey, we are able to determine the column density ratios of six O-COMs with respect to methanol (CH$_3$OH) in a sample of 14 high-mass protostellar sources to investigate their origin through ice and/or gas-phase chemistry. The selected species are: acetaldehyde (CH$_3$CHO), ethanol (C$_2$H$_5$OH), dimethyl ether (DME, CH$_3$OCH$_3$), methyl formate (MF, CH$_3$OCHO), glycolaldehyde (GA, CH$_2$OHCHO), and ethylene glycol (EG, (CH$_2$OH)$_2$). DME and MF have the highest and most constant ratios within one order of magnitude, while the other four species have lower ratios and exhibit larger scatter by 1-2 orders of magnitude. We compare the O-COM ratios of high-mass CoCCoA sources with those of 5 low-mass protostars available from the literature, along with the results from experiments and simulations. We find that the O-COM ratios with respect to methanol are on the same level in both the high- and low-mass samples, which suggests that these species are mainly formed in similar environments during star formation, probably in ice mantles on dust grains during early pre-stellar stages. Current simulations and experiments can reproduce most observational trends with a few exceptions, and hypotheses exist to explain the differences between observations and simulations/experiments, such as the involvement of gas-phase chemistry and different emitting areas of molecules. △ Less

Submitted 4 August, 2023; originally announced August 2023.

Comments: Accepted for publication in Astronomy & Astrophysics. 17 pages and 7 figures in the main body, 15 pages and 8 figures in the appendix

Journal ref: A&A 678, A137 (2023)

Low NH$_{3}$/H$_{2}$O ratio in comet C/2020 F3 (NEOWISE) at 0.7 au from the Sun

Authors: Maria N. Drozdovskaya, Dominique Bockelée-Morvan, Jacques Crovisier, Brett A. McGuire, Nicolas Biver, Steven B. Charnley, Martin A. Cordiner, Stefanie N. Milam, Cyrielle Opitom, Anthony J. Remijan

Abstract: A lower-than-solar elemental nitrogen content has been demonstrated for several comets, including 1P/Halley and 67P/C-G with independent in situ measurements of volatile and refractory budgets. The recently discovered semi-refractory ammonium salts in 67P/C-G are thought to be the missing nitrogen reservoir in comets. The thermal desorption of ammonium salts from cometary dust particles leads to t… ▽ More A lower-than-solar elemental nitrogen content has been demonstrated for several comets, including 1P/Halley and 67P/C-G with independent in situ measurements of volatile and refractory budgets. The recently discovered semi-refractory ammonium salts in 67P/C-G are thought to be the missing nitrogen reservoir in comets. The thermal desorption of ammonium salts from cometary dust particles leads to their decomposition into ammonia and a corresponding acid. The NH$_{3}$/H$_{2}$O ratio is expected to increase with decreasing heliocentric distance with evidence for this in near-infrared observations. NH$_{3}$ has been claimed to be more extended than expected for a nuclear source. Here, the aim is to constrain the NH$_{3}$/H$_{2}$O ratio in comet C/2020 F3 (NEOWISE) during its July 2020 passage. OH emission from comet C/2020 F3 (NEOWISE) was monitored for 2 months with NRT and observed from GBT on 24 July and 11 August 2020. Contemporaneously with the 24 July 2020 OH observations, the NH$_{3}$ hyperfine lines were targeted with GBT. The concurrent GBT and NRT observations allowed the OH quenching radius to be determined at $\left(5.96\pm0.10\right)\times10^{4}$ km on 24 July 2020, which is important for accurately deriving $Q(\text{OH})$. C/2020 F3 (NEOWISE) was a highly active comet with $Q(\text{H}_{2}\text{O}) \approx 2\times10^{30}$ molec s$^{-1}$ one day before perihelion. The $3σ$ upper limit for $Q_{\text{NH}_{3}}/Q_{\text{H}_{2}\text{O}}$ is $<0.29\%$ at $0.7$ au from the Sun. The obtained NH$_{3}$/H$_{2}$O ratio is a factor of a few lower than measurements for other comets at such heliocentric distances. The abundance of NH$_{3}$ may vary strongly with time depending on the amount of water-poor dust in the coma. Lifted dust can be heated, fragmented, and super-heated; whereby, ammonium salts, if present, can rapidly thermally disintegrate and modify the NH$_{3}$/H$_{2}$O ratio. △ Less

Submitted 21 July, 2023; originally announced July 2023.

Comments: Accepted for publication in A&A; 18 pages, 8 figures, 6 tables

Journal ref: A&A 677, A157 (2023)

The edge-on protoplanetary disk HH 48 NE II. Modeling ices and silicates

Authors: J. A. Sturm, M. K. McClure, J. B. Bergner, D. Harsono, E. Dartois, M. N. Drozdovskaya, S. Ioppolo, K. I. Öberg, C. J. Law, M. E. Palumbo, Y. J. Pendleton, W. R. M. Rocha, H. Terada, R. G. Urso

Abstract: The abundance and distribution of ice in protoplanetary disks (PPD) is critical to understand the linkage between the composition of circumstellar matter and the composition of exoplanets. Edge-on PPDs are a useful tool to constrain such ice composition and its location in the disk, as ice spectral signatures can be observed in absorption against the continuum emission arising from the warmer cent… ▽ More The abundance and distribution of ice in protoplanetary disks (PPD) is critical to understand the linkage between the composition of circumstellar matter and the composition of exoplanets. Edge-on PPDs are a useful tool to constrain such ice composition and its location in the disk, as ice spectral signatures can be observed in absorption against the continuum emission arising from the warmer central disk regions. The aim of this work is to model ice absorption features in PPDs and determine how well the abundance of the main ice species across the disk can be determined within the uncertainty of the physical parameter space. The edge-on PPD around HH 48 NE, a target of the JWST ERS program IceAge, is used as a reference system. We use RADMC-3D to raytrace the mid-infrared continuum. Using a constant parameterized ice abundance, ice opacities are added to the dust opacity in regions wherever the disk is cold enough for the main carbon, oxygen and nitrogen carriers to freeze out. The global abundance of the main ice carriers in HH 48 NE can be determined within a factor of 3, when taking the uncertainty of the physical parameters into account. Ice features in PPDs can be saturated at an optical depth <1, due to local saturation. Spatially observed ice optical depths cannot be directly related to column densities due to radiative transfer effects. Vertical snowlines will not be a clear transition due to the radially increasing height of the snowsurface, but their location may be constrained from observations using radiative transfer modeling. Radial snowlines are not really accesible. Not only the ice abundance, but also inclination, settling, grain size distribution and disk mass have strong impact on the observed ice absorption features in disks. Relative changes in ice abundance can be inferred from observations only if the source structure is well constrained △ Less

Submitted 8 May, 2023; v1 submitted 3 May, 2023; originally announced May 2023.

Comments: 18 pages, 12 figures, accepted for publication in Astronomy & Astrophysics

The edge-on protoplanetary disk HH 48 NE I. Modeling the geometry and stellar parameters

Authors: J. A. Sturm, M. K. McClure, C. J. Law, D. Harsono, J. B. Bergner, E. Dartois, M. N. Drozdovskaya, S. Ioppolo, K. I. Öberg, M. E. Palumbo, Y. J. Pendleton, W. R. M. Rocha, H. Terada, R. G. Urso

Abstract: Context. Observations of edge-on disks are an important tool for constraining general protoplanetary disk properties that cannot be determined in any other way. However, most radiative transfer models cannot simultaneously reproduce the spectral energy distributions (SEDs) and resolved scattered light and submillimeter observations of these systems, due to the differences in geometry and dust prop… ▽ More Context. Observations of edge-on disks are an important tool for constraining general protoplanetary disk properties that cannot be determined in any other way. However, most radiative transfer models cannot simultaneously reproduce the spectral energy distributions (SEDs) and resolved scattered light and submillimeter observations of these systems, due to the differences in geometry and dust properties at different wavelengths. Aims. We simultaneously constrain the geometry of the edge-on protoplanetary disk HH 48 NE and the characteristics of the host star. HH 48 NE is part of the JWST early release science program Ice Age. This work serves as a stepping stone towards a better understanding of the disk physical structure and icy chemistry in this particular source. This kind of modeling lays the groundwork for studying other edge-on sources to be observed with the JWST. Methods. We fit a parameterized dust model to HH 48 NE by coupling the radiative transfer code RADMC-3D and an MCMC framework. The dust structure was fitted independently to a compiled SED, a scattered light image at 0.8 $μ$m and an ALMA dust continuum observation at 890 $μ$m. Results. We find that 90% of the dust mass in HH 48 NE is settled to the disk midplane, less than in average disks, and that the atmospheric layers of the disk contain exclusively large grains (0.3-10 $μ$m). The exclusion of small grains in the upper atmosphere likely has important consequences for the chemistry due to the deep penetration of high-energy photons. The addition of a relatively large cavity (ca. 50 au in radius) is necessary to explain the strong mid-infrared emission, and to fit the scattered light and continuum observations simultaneously. △ Less

Submitted 8 May, 2023; v1 submitted 3 May, 2023; originally announced May 2023.

Comments: 16 pages, 8 figures, accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 677, A17 (2023)

Azimuthal C/O Variations in a Planet-Forming Disk

Authors: Luke Keyte, Mihkel Kama, Alice S. Booth, Edwin A. Bergin, L. Ilsedore Cleeves, Ewine F. van Dishoeck, Maria N. Drozdovskaya, Kenji Furuya, Jonathan Rawlings, Oliver Shorttle, Catherine Walsh

Abstract: The elemental carbon-to-oxygen ratio (C/O) in the atmosphere of a giant planet is a promising diagnostic of that planet's formation history in a protoplanetary disk. Alongside efforts in the exoplanet community to measure C/O in planetary atmospheres, observational and theoretical studies of disks are increasingly focused on understanding how the gas-phase C/O varies both with radial location and… ▽ More The elemental carbon-to-oxygen ratio (C/O) in the atmosphere of a giant planet is a promising diagnostic of that planet's formation history in a protoplanetary disk. Alongside efforts in the exoplanet community to measure C/O in planetary atmospheres, observational and theoretical studies of disks are increasingly focused on understanding how the gas-phase C/O varies both with radial location and between disks. This is mostly tied to the icelines of major volatile carriers such as CO and H2O. Using ALMA observations of CS and SO, we have unearthed evidence for an entirely novel type of C/O variation in the protoplanetary disk around HD 100546: an azimuthal variation from a typical, oxygen-dominated ratio (C/O=0.5) to a carbon-dominated ratio (C/O>1.0). We show that the spatial distribution and peculiar line kinematics of both CS and SO molecules can be well-explained by azimuthal variations in the C/O ratio. We propose a shadowing mechanism that could lead to such a chemical dichotomy. Our results imply that tracing the formation history of giant exoplanets using their atmospheric C/O ratios will need to take into account time-dependent azimuthal C/O variations in a planet's accretion zone. △ Less

Submitted 15 March, 2023; originally announced March 2023.

Comments: Accepted in Nature Astronomy

ALMA ACA study of the H$_2$S/OCS ratio in low-mass protostars

Authors: Tanya Kushwahaa, Maria N. Drozdovskaya, Łukasz Tychoniec, Benoît Tabone

Abstract: The identification of the main sulfur reservoir on its way from the diffuse interstellar medium to the cold dense star-forming cores and eventually to protostars is a long-standing problem. Despite sulfur's astrochemical relevance, the abundance of S-bearing molecules in dense cores and regions around protostars is still insufficiently constrained. The goal of this investigation is to derive the g… ▽ More The identification of the main sulfur reservoir on its way from the diffuse interstellar medium to the cold dense star-forming cores and eventually to protostars is a long-standing problem. Despite sulfur's astrochemical relevance, the abundance of S-bearing molecules in dense cores and regions around protostars is still insufficiently constrained. The goal of this investigation is to derive the gas-phase H$_2$S/OCS ratio for several low-mass protostars, which could provide crucial information about the physical and chemical conditions in the birth cloud of Sun-like stars. Using ALMA ACA Band 6 observations, H$_2$S, OCS, and their isotopologs are searched for in 10 Class 0/I protostars with different source properties such as age, mass, and environmental conditions. An LTE model is used to fit synthetic spectra to the detected lines and to derive the column densities based solely on optically thin lines. The H$_2$S and OCS column densities span four orders of magnitude across the sample. The H$_2$S/OCS ratio is found to be in the range from 0.2 to above 9.7. IRAS 16293-2422 A and Ser-SMM3 have the lowest ratio, while BHR71-IRS1 has the highest. Only the H$_2$S/OCS ratio of BHR71-IRS1 agress within uncertainties with the ratio in comet 67P/C$-$G. The determined gas-phase H$_2$S/OCS ratios can be below the upper limits on the solid-state ratios by as much as an order of magnitude. The H$_2$S/OCS ratio depends significantly on the environment of the birth cloud, such as UV-irradiation and heating received prior to the formation of a protostar. The highly isolated birth environment of BHR71-IRS1 is hypothesized to be the reason for its high gaseous H$_2$S/OCS ratio due to lower rates of photoreactions and more efficient hydrogenation reactions under such dark, cold conditions. The gaseous inventory of S-bearing molecules in BHR71-IRS1 appears to be most similar to that of interstellar ices. △ Less

Submitted 18 February, 2023; originally announced February 2023.

Comments: Accepted for publication in A&A; 32 pages, 24 figures, 18 tables

Journal ref: A&A 672, A122 (2023)

An Ice Age JWST inventory of dense molecular cloud ices

Authors: M. K. McClure, W. R. M. Rocha, K. M. Pontoppidan, N. Crouzet, L. E. U. Chu, E. Dartois, T. Lamberts, J. A. Noble, Y. J. Pendleton, G. Perotti, D. Qasim, M. G. Rachid, Z. L. Smith, Fengwu Sun, Tracy L Beck, A. C. A. Boogert, W. A. Brown, P. Caselli, S. B. Charnley, Herma M. Cuppen, H. Dickinson, M. N. Drozdovskaya, E. Egami, J. Erkal, H. Fraser , et al. (17 additional authors not shown)

Abstract: Icy grain mantles are the main reservoir of the volatile elements that link chemical processes in dark, interstellar clouds with the formation of planets and composition of their atmospheres. The initial ice composition is set in the cold, dense parts of molecular clouds, prior to the onset of star formation. With the exquisite sensitivity of JWST, this critical stage of ice evolution is now acces… ▽ More Icy grain mantles are the main reservoir of the volatile elements that link chemical processes in dark, interstellar clouds with the formation of planets and composition of their atmospheres. The initial ice composition is set in the cold, dense parts of molecular clouds, prior to the onset of star formation. With the exquisite sensitivity of JWST, this critical stage of ice evolution is now accessible for detailed study. Here we show the first results of the Early Release Science program "Ice Age" that reveal the rich composition of these dense cloud ices. Weak ices, including, $^{13}$CO$_2$, OCN$^-$, $^{13}$CO, OCS, and COMs functional groups are now detected along two pre-stellar lines of sight. The $^{12}$CO$_2$ ice profile indicates modest growth of the icy grains. Column densities of the major and minor ice species indicate that ices contribute between 2 and 19% of the bulk budgets of the key C, O, N, and S elements. Our results suggest that the formation of simple and complex molecules could begin early in a water-ice rich environment. △ Less

Submitted 22 January, 2023; originally announced January 2023.

Comments: To appear in Nature Astronomy on January 23rd, 2023. 33 pages, 16 figures, 3 tables; includes extended and supplemental data sections. Part of the JWST Ice Age Early Release Science program's science enabling products. Enhanced spectra downloadable on Zenodo at the following DOI: 10.5281/zenodo.7501239

Interstellar Heritage and the Birth Environment of the Solar System

Authors: Edwin A. Bergin, Conel Alexander, Maria Drozdovskaya, Matthieu Gounelle, Susanne Pfalzner

Abstract: In this chapter, we explore the origins of cometary material and discuss the clues cometary composition provides in the context of the origin of our solar system. The review focuses on both cometary refractory and volatile materials, which jointly provide crucial information about the processes that shaped the solar system into what it is today. Both areas have significantly advanced over the past… ▽ More In this chapter, we explore the origins of cometary material and discuss the clues cometary composition provides in the context of the origin of our solar system. The review focuses on both cometary refractory and volatile materials, which jointly provide crucial information about the processes that shaped the solar system into what it is today. Both areas have significantly advanced over the past decade. We also view comets more broadly and discuss compositions considering laboratory studies of cometary materials, including interplanetary dust particles and meteoritic material that are potential cometary samples, along with meteorites, and in situ/remote studies of cometary comae. In our review, we focus on key areas from elemental/molecular compositions, isotopic ratios, carbonaceous and silicate refractories, short-lived radionuclides, and solar system dynamics that can be used as probes of the solar birth environment. We synthesize this data that points towards the birth of our solar system in a clustered star-forming environment. △ Less

Submitted 12 January, 2023; originally announced January 2023.

Comments: 31 pages, 7 figures, 2 tables; Chapter in press for the book Comets III, edited by K. Meech and M. Combi, University of Arizona Press

The mid-infrared molecular inventory towards Orion IRc2

Authors: Sarah Nickerson, Naseem Rangwala, Sean W. J. Colgan, Curtis DeWitt, Jose S. Monzon, Xinchuan Huang, Kinsuk Acharyya, Maria N. Drozdovskaya, Ryan C. Fortenberry, Eric Herbst, Timothy J. Lee

Abstract: We present the first high spectral resolution mid-infrared survey in the Orion BN/KL region, covering 7.2 to 28.3 micron. With SOFIA/EXES we target the enigmatic source Orion IRc2. While this is in the most prolifically studied massive star-forming region, longer wavelengths and molecular emission lines dominated previous spectral surveys. The mid-infrared observations in this work access differen… ▽ More We present the first high spectral resolution mid-infrared survey in the Orion BN/KL region, covering 7.2 to 28.3 micron. With SOFIA/EXES we target the enigmatic source Orion IRc2. While this is in the most prolifically studied massive star-forming region, longer wavelengths and molecular emission lines dominated previous spectral surveys. The mid-infrared observations in this work access different components and molecular species in unprecedented detail. We unambiguously identify two new kinematic components, both chemically rich with multiple molecular absorption lines. The "blue clump" has vLSR = -7.1 \pm 0.7 km/s and the "red clump" 1.4 \pm 0.5 km/s. While the blue and red clumps have similar temperatures and line widths, molecular species in the blue clump have higher column densities. They are both likely linked to pure rotational H2 emission also covered by this survey. This work provides evidence for the scenario that the blue and red clumps are distinct components unrelated to the classic components in the Orion BN/KL region. Comparison to spectroscopic surveys towards other infrared targets in the region show that the blue clump is clearly extended. We analyze, compare, and present in depth findings on the physical conditions of C2H2, 13CCH2, CH4, CS, H2O, HCN, H13CN, HNC, NH3, and SO2 absorption lines and an H2 emission line associated with the blue and red clumps. We also provide limited analysis of H2O and SiO molecular emission lines towards Orion IRc2 and the atomic forbidden transitions [FeII], [SI], [SIII], and [NeII]. △ Less

Submitted 28 November, 2022; originally announced November 2022.

Comments: Accepted to ApJ; 44 pages, 14 figures, 13 tables

Fevering Interstellar Ices Have More CH3OD

Authors: Beatrice M. Kulterer, Maria N. Drozdovskaya, Stefano Antonellini, Catherine Walsh, Tom J. Millar

Abstract: Mono-deuterated methanol is thought to form during the prestellar core stage of star formation. Observed variations in the CH2DOH/CH3OD ratio suggest that its formation is strongly dependent on the surrounding cloud conditions. Thus, it is a potential tracer of the physical conditions before the onset of star formation. A single-point physical model representative of a typical prestellar core is c… ▽ More Mono-deuterated methanol is thought to form during the prestellar core stage of star formation. Observed variations in the CH2DOH/CH3OD ratio suggest that its formation is strongly dependent on the surrounding cloud conditions. Thus, it is a potential tracer of the physical conditions before the onset of star formation. A single-point physical model representative of a typical prestellar core is coupled to chemical models to investigate potential formation pathways towards deuterated methanol at the prestellar stage. Simple addition reactions of H and D are not able to reproduce observed abundances. The implementation of an experimentally verified abstraction scheme leads to the efficient formation of methyl-deuterated methanol, but lacks sufficient formation of hydroxy-deuterated methanol. CH3OD is most likely formed at a later evolutionarymstage, potentially from H-D exchange reactions in warm ices between HDO (and D2O) and CH3OH. The CH2DOH/CH3OD ratio is not an appropriate tracer of the physical conditions during the prestellar stage, but might be better suited as a tracer of ice heating. △ Less

Submitted 3 April, 2022; originally announced April 2022.

Comments: Accepted for publication in ACS Earth and Space Chemistry; 56 pages, 7 figures

Successive deuteration in low-mass star-forming regions: the case of D$_{2}$-methanol (CHD$_{2}$OH) in IRAS 16293-2422

Authors: Maria N. Drozdovskaya, Laurent H. Coudert, Laurent Margulès, Audrey Coutens, Jes K. Jørgensen, Sébastien Manigand

Abstract: Accurate quantification of the column density of di-deuterated methanol is a key missing puzzle piece in the otherwise thoroughly constrained family of D-bearing methanol in the deeply embedded low-mass protostellar system and astrochemical template source IRAS16293-2422. A spectroscopic dataset for astrophysical purposes is built for CHD$_{2}$OH and made publicly available to facilitate accurate… ▽ More Accurate quantification of the column density of di-deuterated methanol is a key missing puzzle piece in the otherwise thoroughly constrained family of D-bearing methanol in the deeply embedded low-mass protostellar system and astrochemical template source IRAS16293-2422. A spectroscopic dataset for astrophysical purposes is built for CHD$_{2}$OH and made publicly available to facilitate accurate characterization of this species in astrochemical surveys. The newly computed line list and partition function are used to search for CHD$_{2}$OH towards IRAS16293-2422 A and B in data from ALMA-PILS. Only non-blended, optically thin lines of CHD$_{2}$OH are used for the synthetic spectral fitting. The constructed spectroscopic database contains line frequencies and strengths for 7417 transitions in the 0 to 500 GHz frequency range. ALMA-PILS observations in the 329-363 GHz range are used to identify 105 unique, non-blended, optically thin line frequencies of CHD$_{2}$OH for synthetic spectral fitting. The derived excitation temperatures and column densities yield high D/H ratios of CHD$_{2}$OH in IRAS 16293-2422 A and B of 7.5$\pm$1.1% and 7.7$\pm$1.2%, respectively. Deuteration in IRAS 16293-2422 is not higher than in other low-mass star-forming regions. Di-deuterated molecules consistently have higher D/H ratios than their monodeuterated counterparts in all low-mass protostars, which may be a natural consequence of H-D substitution reactions as seen in laboratory experiments. The Solar System's natal cloud, as traced by comet 67P/Churyumov-Gerasimenko, may have had a lower initial abundance of D, been warmer than the cloud of IRAS16293-2422, or been partially reprocessed. In combination with accurate spectroscopy, a careful spectral analysis, and a consideration of the underlying assumptions, successive deuteration is a robust window on the physicochemical provenance of star-forming systems. △ Less

Submitted 18 January, 2022; originally announced January 2022.

Comments: Accepted for publication in A&A; 24 pages, 18 figures, 3 tables; the full line list is made available as an electronic Table at the Centre de Données astronomiques de Strasbourg (CDS) with association to this A&A article

Journal ref: A&A 659, A69 (2022)

Molecule dependent oxygen isotopic ratios in the coma of comet 67P/Churyumov-Gerasimenko

Authors: K. Altwegg, H. Balsiger, M. Combi, J. De Keyser, M. N. Drozdovskaya, S. A. Fuselier, T. I. Gombosi, N. Hänni, M. Rubin, M. Schuhmann, I. Schroeder, S. Wampfler

Abstract: The ratios of the three stable oxygen isotopes 16O, 17O and 18O on Earth and, as far as we know in the solar system, show variations on the order of a few percent at most, with a few outliers in meteorites. However, in the interstellar medium there are some highly fractionated oxygen isotopic ratios in some specific molecules. The goal of this work is to investigate the oxygen isotopic ratios in d… ▽ More The ratios of the three stable oxygen isotopes 16O, 17O and 18O on Earth and, as far as we know in the solar system, show variations on the order of a few percent at most, with a few outliers in meteorites. However, in the interstellar medium there are some highly fractionated oxygen isotopic ratios in some specific molecules. The goal of this work is to investigate the oxygen isotopic ratios in different volatile molecules found in the coma of comet 67P/Churyumov-Gerasimenko and compare them with findings from interstellar clouds in order to assess commonalities and differences. To accomplish this goal, we analyzed data from the ROSINA instrument on Rosetta during its mission around the comet. 16O/18O ratios could be determined for O2, methanol, formaldehyde, carbonyl sulfide and sulfur monoxide/dioxide. For O2 the 16O/17O ratio is also available. Some ratios are strongly enriched in the heavy isotopes, especially for sulfur bearing molecules and formaldehyde, whereas for methanol the ratios are compatible with the ones in the solar system. O2 falls in-between, but its oxygen isotopic ratios clearly differ from water, which likely rules out an origin of O2 from water, be it by radiolysis, dismutation during sublimation or the Eley-Rideal process from water ions hitting the nucleus as postulated in the literature. △ Less

Submitted 10 December, 2020; originally announced December 2020.

Journal ref: 2020, Monthly Notices of the Royal Astronomical Society, 498(4), 5855-5862

Prestellar grain-surface origins of deuterated methanol in comet 67P/Churyumov-Gerasimenko

Authors: Maria N. Drozdovskaya, Isaac R. H. G. Schroeder, Martin Rubin, Kathrin Altwegg, Ewine F. van Dishoeck, Beatrice M. Kulterer, Johan De Keyser, Stephen A. Fuselier, Michael Combi

Abstract: Deuterated methanol is one of the most robust windows astrochemists have on the individual chemical reactions forming deuterium-bearing molecules and the physicochemical history of the regions where they reside. The first-time detection of mono- and di-deuterated methanol in a cometary coma is presented for comet 67P/Churyumov-Gerasimenko using Rosetta-ROSINA data. D-methanol (CH3OD and CH2DOH com… ▽ More Deuterated methanol is one of the most robust windows astrochemists have on the individual chemical reactions forming deuterium-bearing molecules and the physicochemical history of the regions where they reside. The first-time detection of mono- and di-deuterated methanol in a cometary coma is presented for comet 67P/Churyumov-Gerasimenko using Rosetta-ROSINA data. D-methanol (CH3OD and CH2DOH combined) and D2-methanol (CH2DOD and CHD2OH combined) have an abundance of 5.5+/-0.46 and 0.00069+/-0.00014 per cent relative to normal methanol. The data span a methanol deuteration fraction (D/H ratio) in the 0.71-6.6 per cent range, accounting for statistical corrections for the location of D in the molecule and including statistical error propagation in the ROSINA measurements. It is argued that cometary CH2DOH forms from CO hydrogenation to CH3OH and subsequent H-D substitution reactions in CH3-R. CHD2OH is likely produced from deuterated formaldehyde. Meanwhile, CH3OD and CH2DOD, could form via H-D exchange reactions in OH-R in the presence of deuterated water ice. Methanol formation and deuteration is argued to occur at the same epoch as D2O formation from HDO, with formation of mono-deuterated water, hydrogen sulfide, and ammonia occurring prior to that. The cometary D-methanol/methanol ratio is demonstrated to agree most closely with that in prestellar cores and low-mass protostellar regions. The results suggest that cometary methanol stems from the innate cold (10-20 K) prestellar core that birthed our Solar System. Cometary volatiles individually reflect the evolutionary phases of star formation from cloud to core to protostar. △ Less

Submitted 23 October, 2020; originally announced October 2020.

Comments: Accepted for publication in MNRAS; 29 pages, 8 figures, 4 tables

The First Mid-Infrared Detection of HNC in the Interstellar Medium: Probing the Extreme Environment Towards the Orion Hot Core

Authors: Sarah Nickerson, Naseem Rangwala, Sean Colgan, Curtis DeWitt, Xinchuan Huang, Kinsuk Acharyya, Maria Drozdovskaya, Ryan C. Fortenberry, Eric Herbst, Timothy J. Lee

Abstract: We present the first mid-infrared (MIR) detections of HNC and H13CN in the interstellar medium, and numerous, resolved HCN rovibrational transitions. Our observations span 12.8 to 22.9 micron towards the hot core Orion IRc2, obtained with the Echelon-Cross-Echelle Spectrograph aboard the Stratospheric Observatory for Infrared Astronomy (SOFIA). Exceptional, ~5 km/s, resolution distinguishes indivi… ▽ More We present the first mid-infrared (MIR) detections of HNC and H13CN in the interstellar medium, and numerous, resolved HCN rovibrational transitions. Our observations span 12.8 to 22.9 micron towards the hot core Orion IRc2, obtained with the Echelon-Cross-Echelle Spectrograph aboard the Stratospheric Observatory for Infrared Astronomy (SOFIA). Exceptional, ~5 km/s, resolution distinguishes individual rovibrational transitions of the HNC and HCN P, Q, and R branches; and the H13CN R branch. This allows direct measurement of the species' excitation temperatures, column densities, and relative abundances. HNC and H13CN exhibit a local standard rest velocity of -7 km/s that may be associated with an outflow from nearby radio source I and an excitation temperature of about 100 K. We resolve two velocity components for HCN, the primary component also being at -7 km/s with temperature 165 K. The hottest component, which had never before been observed, is at 1 km/s with temperature 309 K. This is the closest component to the hot core's centre measured to date. The derived 12C/13C=13 is below expectation for Orion's Galactocentric distance, but the derived HCN/HNC=72 is expected for this extreme environment. Compared to previous sub-mm and mm observations, our SOFIA line survey of this region shows that the resolved MIR molecular transitions are probing a distinct physical component and isolating the chemistry closest to the hot core. △ Less

Submitted 28 August, 2020; originally announced August 2020.

Comments: 20 pages, 6 figures, submitted to ApJ, comments welcome

Physicochemical models: source-tailored or generic?

Authors: Beatrice M. Kulterer, Maria N. Drozdovskaya, Audrey Coutens, Sébastien Manigand, Gwendoline Stéphan

Abstract: Physicochemical models can be powerful tools to trace the chemical evolution of a protostellar system and allow to constrain its physical conditions at formation. The aim of this work is to assess whether source-tailored modelling is needed to explain the observed molecular abundances around young, low-mass protostars or if, and to what extent, generic models can improve our understanding of the c… ▽ More Physicochemical models can be powerful tools to trace the chemical evolution of a protostellar system and allow to constrain its physical conditions at formation. The aim of this work is to assess whether source-tailored modelling is needed to explain the observed molecular abundances around young, low-mass protostars or if, and to what extent, generic models can improve our understanding of the chemistry in the earliest stages of star formation. The physical conditions and the abundances of simple, most abundant molecules based on three models are compared. After establishing the discrepancies between the calculated chemical output, the calculations are redone with the same chemical model for all three sets of physical input parameters. With the differences arising from the chemical models eliminated, the output is compared based on the influence of the physical model. Results suggest that the impact of the chemical model is small compared to the influence of the physical conditions, with considered timescales having the most drastic effect. Source-tailored models may be simpler by design; however, likely do not sufficiently constrain the physical and chemical parameters within the global picture of star-forming regions. Generic models with more comprehensive physics may not provide the optimal match to observations of a particular protostellar system, but allow a source to be studied in perspective of other star-forming regions. △ Less

Submitted 11 August, 2020; originally announced August 2020.

The ALMA-PILS survey: First detection of the unsaturated 3-carbon molecules Propenal (C$_2$H$_3$CHO) and Propylene (C$_3$H$_6$) towards IRAS 16293$-$2422 B

Authors: S. Manigand, A. Coutens, J. -C. Loison, V. Wakelam, H. Calcutt, H. S. P. Müller, J. K. Jørgensen, V. Taquet, S. F. Wampfler, T. L. Bourke, B. M. Kulterer, E. F. van Dishoeck, M. N. Drozdovskaya, N. F. W. Ligterink

Abstract: Complex organic molecules with three carbon atoms are found in the earliest stages of star formation. In particular, propenal (C$_2$H$_3$CHO) is a species of interest due to its implication in the formation of more complex species and even biotic molecules. This study aims to search for the presence of C$_2$H$_3$CHO and other three-carbon species such as propylene (C$_3$H$_6$) in the hot corino re… ▽ More Complex organic molecules with three carbon atoms are found in the earliest stages of star formation. In particular, propenal (C$_2$H$_3$CHO) is a species of interest due to its implication in the formation of more complex species and even biotic molecules. This study aims to search for the presence of C$_2$H$_3$CHO and other three-carbon species such as propylene (C$_3$H$_6$) in the hot corino region of the low-mass protostellar binary IRAS 16293--2422 to understand their formation pathways. We use ALMA observations in Band 6 and 7 from various surveys to search for the presence of C$_3$H$_6$ and C$_2$H$_3$CHO towards the protostar IRAS 16293--2422 B (IRAS 16293B). We report the detection of both C$_3$H$_6$ and C$_2$H$_3$CHO towards IRAS 16293B, however, no unblended lines were found towards the other component of the binary system, IRAS 16293A. We derive column density upper limits for C$_3$H$_8$, HCCCHO, n-C$_3$H$_7$OH, i-C$_3$H$_7$OH, C$_3$O, and cis-HC(O)CHO towards IRAS 16293B. We then use a three-phase chemical model to simulate the formation of these species in a typical prestellar environment followed by its hydrodynamical collapse until the birth of the central protostar. Different formation paths, such as successive hydrogenation and radical-radical additions on grain surfaces, are tested and compared to the observational results. The simulations reproduce the abundances within one order of magnitude from those observed towards IRAS 16293B, with the best agreement found for a rate of $10^{-12}$ cm$^3$ s$^{-1}$ for the gas-phase reaction C$_3$ + O $\rightarrow$ C$_2$ + CO. Successive hydrogenations of C$_3$, HC(O)CHO, and CH$_3$OCHO on grain surfaces are a major and crucial formation route of complex organics molecules, whereas both successive hydrogenation pathways and radical-radical addition reactions contribute to the formation of C$_2$H$_5$CHO. △ Less

Submitted 8 July, 2020; originally announced July 2020.

Comments: Accepted for publication in A&A. 24 pages, 17 figures

Journal ref: A&A 645, A53 (2021)

The ALMA-PILS survey: Inventory of complex organic molecules towards IRAS 16293-2422 A

Authors: S. Manigand, J. K. Jørgensen, H. Calcutt, H. S. P. Müller, N. F. W. Ligterink, A. Coutens, M. N. Drozdovskaya, E. F. van Dishoeck, S. F. Wampfler

Abstract: Complex organic molecules (COM) are detected in many sources in the warm inner regions of envelopes surrounding deeply embedded protostars. Exactly how these COM form remains an open question. This study aims to constrain the formation of complex organic molecules through comparisons of their abundances towards the Class 0 protostellar binary IRAS 16293-2422 (IRAS16293). We utilised observations f… ▽ More Complex organic molecules (COM) are detected in many sources in the warm inner regions of envelopes surrounding deeply embedded protostars. Exactly how these COM form remains an open question. This study aims to constrain the formation of complex organic molecules through comparisons of their abundances towards the Class 0 protostellar binary IRAS 16293-2422 (IRAS16293). We utilised observations from the ALMA Protostellar Interferometric Line Survey of IRAS16293. The species identification and the rotational temperature and column density estimation were derived by fitting the extracted spectra towards IRAS16293 A and IRAS16293 B with synthetic spectra. The majority of the work in this paper pertains to the analysis of IRAS16293 A for a comparison with the results from the other binary component, which have already been published. We detect 15 different COM, as well as 16 isotopologues towards the most luminous companion protostar IRAS16293 A. Tentative detections of an additional 11 isotopologues are reported. We also searched for and report on the first detections of CH3OCH2OH and t-C2H5OCH3 towards IRAS16293 B and the follow-up detection of CH2DCHO and CH3CDO. Twenty-four lines of CHD2OH are also identified. The comparison between the two protostars of the binary system shows significant differences in abundance for some of the species, which are partially correlated to their spatial distribution. The spatial distribution is consistent with the sublimation temperature of the species; those with higher expected sublimation temperatures are located in the most compact region of the hot corino towards IRAS16293 A. This spatial differentiation is not resolved in IRAS16293 B and will require observations at a higher angular resolution. In parallel, the list of identified CHD2OH lines shows the need of accurate spectroscopic data including their line strength. △ Less

Submitted 17 January, 2020; originally announced January 2020.

Comments: Accepted for publication in A&A. 31 pages, 17 figures

Journal ref: A&A 635, A48 (2020)

ALMA and ROSINA detections of phosphorus-bearing molecules: the interstellar thread between star-forming regions and comets

Authors: V. M. Rivilla, M. N. Drozdovskaya, K. Altwegg, P. Caselli, M. T. Beltrán, F. Fontani, F. F. S. van der Tak, R. Cesaroni, A. Vasyunin, M. Rubin, F. Lique, S. Marinakis, L. Testi, the ROSINA team

Abstract: To understand how Phosphorus-bearing molecules are formed in star-forming regions, we have analysed ALMA observations of PN and PO towards the massive star-forming region AFGL 5142, combined with a new analysis of the data of the comet 67P/Churyumov-Gerasimenko taken with the ROSINA instrument onboard Rosetta. The ALMA maps show that the emission of PN and PO arises from several spots associated w… ▽ More To understand how Phosphorus-bearing molecules are formed in star-forming regions, we have analysed ALMA observations of PN and PO towards the massive star-forming region AFGL 5142, combined with a new analysis of the data of the comet 67P/Churyumov-Gerasimenko taken with the ROSINA instrument onboard Rosetta. The ALMA maps show that the emission of PN and PO arises from several spots associated with low-velocity gas with narrow linewidths in the cavity walls of a bipolar outflow. PO is more abundant than PN in most of the spots, with the PO/PN ratio increasing as a function of the distance to the protostar. Our data favor a formation scenario in which shocks sputter phosphorus from the surface of dust grains, and gas-phase photochemistry induced by UV photons from the protostar allows efficient formation of the two species in the cavity walls. Our analysis of the ROSINA data has revealed that PO is the main carrier of P in the comet, with PO/PN>10. Since comets may have delivered a significant amount of prebiotic material to the early Earth, this finding suggests that PO could contribute significantly to the phosphorus reservoir during the dawn of our planet. There is evidence that PO was already in the cometary ices prior to the birth of the Sun, so the chemical budget of the comet might be inherited from the natal environment of the Solar System, which is thought to be a stellar cluster including also massive stars. △ Less

Submitted 26 November, 2019; originally announced November 2019.

Comments: Accepted in MNRAS (20 pages, 12 figures)

The ALMA-PILS survey: propyne (CH$_3$CCH) in IRAS 16293$-$2422

Authors: H. Calcutt, E. R. Willis, J. K. Jørgensen, P. Bjerkeli, N. F. W. Ligterink, A. Coutens, H. S. P. Müller, R. T. Garrod, S. F. Wampfler, M. N. Drozdovskaya

Abstract: Context. Propyne (CH$_3$CCH) has been detected in a variety of environments, from Galactic star-forming regions to extragalactic sources. Such molecules are excellent tracers of the physical conditions in star-forming regions. Aims. This study explores the emission of CH$_3$CCH in the low-mass protostellar binary, IRAS 16293$-$2422, examining the spatial scales traced by this molecule, as well a… ▽ More Context. Propyne (CH$_3$CCH) has been detected in a variety of environments, from Galactic star-forming regions to extragalactic sources. Such molecules are excellent tracers of the physical conditions in star-forming regions. Aims. This study explores the emission of CH$_3$CCH in the low-mass protostellar binary, IRAS 16293$-$2422, examining the spatial scales traced by this molecule, as well as its formation and destruction pathways. Methods. ALMA observations from the Protostellar Interferometric Line Survey (PILS) are used to determine the abundances and excitation temperatures of CH$_3$CCH towards both protostars, exploring spatial scales from 70 to 2400 au. The three-phase chemical kinetics model MAGICKAL is also used, to explore the chemical reactions of this molecule. Results. CH$_3$CCH is detected towards both IRAS 16293A and IRAS 16293B and is found to trace the hot corino component around each source in the PILS dataset. Eighteen transitions above 3$σ$ are detected, enabling robust excitation temperatures and column densities to be determined in each source. In IRAS 16293A, an excitation temperature of 90 K and a column density of 7.8$\times$10$^{15}$ cm$^{-2}$ best fits the spectra. In IRAS 16293B, an excitation temperature of 100 K and 6.8$\times$10$^{15}$ cm$^{-2}$ best fits the spectra. The chemical modelling finds that in order to reproduce the observed abundances, both gas-phase and grain-surface reactions are needed. Conclusions. CH$_3$CCH is a molecule whose brightness and abundance in many different regions can be utilised to provide a benchmark of molecular variation with the physical properties of star-forming regions. It is essential when making such comparisons, that the abundances are determined with a good understanding of the spatial scale of the emitting region, to ensure that accurate abundances are derived. △ Less

Submitted 29 September, 2019; originally announced September 2019.

Comments: Accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 631, A137 (2019)

Ingredients for Solar-like Systems: protostar IRAS 16293-2422 B versus comet 67P/Churyumov-Gerasimenko

Authors: Maria N. Drozdovskaya, Ewine F. van Dishoeck, Martin Rubin, Jes K. Jørgensen, Kathrin Altwegg

Abstract: Our modern day Solar System has $4.6\times10^9$ yrs of evolution behind it with just a few relics of its birth conditions remaining. Comets are thought to be some of the most pristine tracers of the initial ingredients that were combined to produce the Earth and the other planets. Other low-mass protostars may be analogous to our proto-Sun and hence, could be used to study the building blocks nece… ▽ More Our modern day Solar System has $4.6\times10^9$ yrs of evolution behind it with just a few relics of its birth conditions remaining. Comets are thought to be some of the most pristine tracers of the initial ingredients that were combined to produce the Earth and the other planets. Other low-mass protostars may be analogous to our proto-Sun and hence, could be used to study the building blocks necessary to form Solar-like systems. This study tests this idea on the basis of new high sensitivity, high spatial resolution ALMA data on the protoplanetary disc-scales (~70 au) of IRAS 16293-2422 and the bulk composition of comet 67P/Churyumov-Gerasimenko, as determined for the first time with the unique in situ monitoring carried out by Rosetta. The comparative analysis of the observations from the Protostellar Interferometric Line Survey (PILS) and the measurements made with Rosetta Orbiter Spectrometer for Ion and Neutral Analysis (ROSINA) shows that the relative abundances of CHO-, N-, and S-bearing molecules correlate, with some scatter, between protostellar and cometary data. A tentative correlation is seen for the first time for P- and Cl-bearing compounds. The results imply that the volatile composition of cometesimals and planetesimals is partially inherited from the pre- and protostellar phases of evolution. △ Less

Submitted 29 August, 2019; originally announced August 2019.

Comments: Accepted for publication in MNRAS; 36 pages, 15 figures, 4 tables

Volatile species in comet 67P/Churyumov-Gerasimenko -- investigating the link from the ISM to the terrestrial planets

Authors: Martin Rubin, David V. Bekaert, Michael W. Broadley, Maria N. Drozdovskaya, Susanne F. Wampfler

Abstract: Comets contain abundant amounts of organic and inorganic species. Many of the volatile molecules in comets have also been observed in the interstellar medium and some of them even with similar relative abundances, indicating formation under similar conditions or even sharing a common chemical pathway. There is a growing amount of evidence that suggests comets inherit and preserve substantial fract… ▽ More Comets contain abundant amounts of organic and inorganic species. Many of the volatile molecules in comets have also been observed in the interstellar medium and some of them even with similar relative abundances, indicating formation under similar conditions or even sharing a common chemical pathway. There is a growing amount of evidence that suggests comets inherit and preserve substantial fractions of materials inherited from previous evolutionary phases, potentially indicating that commonplace processes occurred throughout comet-forming regions. Through impacts, part of this material has also been transported to the inner planetary system, including the terrestrial planets. While comets have been ruled out as a major contributor to terrestrial ocean water, substantial delivery of volatile species to the Earth's atmosphere, and as a consequence also organic molecules to its biomass, appears more likely. Comets contain many species of pre-biotic relevance and molecules that are related to biological processes on Earth, and have hence been proposed as potential indicators for the presence of biological processes in the search of extraterrestrial life. While the delivery of cometary material to Earth may have played a crucial role in the emergence of life, the presence of such alleged biosignature molecules in the abiotical environment of comets complicates the detection of life elsewhere in the universe. △ Less

Submitted 8 August, 2019; originally announced August 2019.

Comments: This document is the unedited Author's version of a Submitted Work that was subsequently accepted for publication in ACS Earth and Space Chemistry, copyright American Chemical Society after peer review

AMBITION -- Comet Nucleus Cryogenic Sample Return (White paper for ESA's Voyage 2050 programme)

Authors: D. Bockelée-Morvan, G. Filacchione, K. Altwegg, E. Bianchi, M. Bizzarro, J. Blum, L. Bonal, F. Capaccioni, C. Codella, M. Choukroun, H. Cottin, B. Davidsson, M. C. De Sanctis, M. Drozdovskaya, C. Engrand, M. Galand, C. Güttler, P. Henri, A. Herique, S. Ivanoski, R. Kokotanekova, A. -C. Levasseur-Regourd, K. E. Miller, A. Rotundi, M. Schönbächler , et al. (5 additional authors not shown)

Abstract: This white paper proposes that AMBITION, a Comet Nucleus Sample Return mission, be a cornerstone of ESA's Voyage 2050 programme. We summarise some of the most important questions still open in cometary science after the successes of the Rosetta mission, many of which require sample analysis using techniques that are only possible in laboratories on Earth. We then summarise measurements, instrument… ▽ More This white paper proposes that AMBITION, a Comet Nucleus Sample Return mission, be a cornerstone of ESA's Voyage 2050 programme. We summarise some of the most important questions still open in cometary science after the successes of the Rosetta mission, many of which require sample analysis using techniques that are only possible in laboratories on Earth. We then summarise measurements, instrumentation and mission scenarios that can address these questions, with a recommendation that ESA select an ambitious cryogenic sample return mission. Rendezvous missions to Main Belt comets and Centaurs are compelling cases for M-class missions, expanding our knowledge by exploring new classes of comets. AMBITION would engage a wide community, drawing expertise from a vast range of disciplines within planetary science and astrophysics. With AMBITION, Europe will continue its leadership in the exploration of the most primitive Solar System bodies. △ Less

Submitted 25 July, 2019; originally announced July 2019.

Comments: 8 figures, 26 pages

Elemental and molecular abundances in comet 67P/Churyumov-Gerasimenko

Authors: Martin Rubin, Kathrin Altwegg, Hans Balsiger, Jean-Jacques Berthelier, Michael R. Combi, Johan De Keyser, Maria Drozdovskaya, Björn Fiethe, Stephen A. Fuselier, Sébastien Gasc, Tamas I. Gombosi, Nora Hänni, Kenneth C. Hansen, Urs Mall, Henri Rème, Isaac R. H. G. Schroeder, Markus Schuhmann, Thierry Sémon, Jack H. Waite, Susanne F. Wampfler, Peter Wurz

Abstract: Comets are considered to be some of the most pristine and unprocessed solar system objects accessible to in-situ exploration. Investigating their molecular and elemental composition takes us on a journey back to the early period of our solar system and possibly even further. In this work, we deduce the bulk abundances of the major volatile species in comet 67P/Churyumov-Gerasimenko, the target of… ▽ More Comets are considered to be some of the most pristine and unprocessed solar system objects accessible to in-situ exploration. Investigating their molecular and elemental composition takes us on a journey back to the early period of our solar system and possibly even further. In this work, we deduce the bulk abundances of the major volatile species in comet 67P/Churyumov-Gerasimenko, the target of the European Space Agency's Rosetta mission. The basis are measurements obtained with the ROSINA instrument suite on board the Rosetta orbiter during a suitable period of high outgassing near perihelion. The results are combined with both gas and dust composition measurements published in the literature. This provides an integrated inventory of the major elements present in the nucleus of 67P/Churyumov-Gerasimenko. Similar to comet 1P/Halley, which was visited by ESA's Giotto spacecraft in 1986, comet 67P/Churyumov-Gerasimenko also shows near-solar abundances of oxygen and carbon, whereas hydrogen and nitrogen are depleted compared to solar. Still, the degree of devolatilization is lower than that of inner solar system objects, including meteorites and the Earth. This supports the idea that comets are among the most pristine objects in our solar system. △ Less

Submitted 12 August, 2019; v1 submitted 25 July, 2019; originally announced July 2019.

Comments: Accepted for publication in MNRAS

The ALMA-PILS survey: Gas dynamics in IRAS 16293$-$2422 and the connection between its two protostars

Authors: M. H. D. van der Wiel, S. K. Jacobsen, J. K. Jørgensen, T. L. Bourke, L. E. Kristensen, P. Bjerkeli, N. M. Murillo, H. Calcutt, H. S. P. Müller, A. Coutens, M. N. Drozdovskaya, C. Favre, S. F. Wampfler

Abstract: [Abridged] The majority of stars form in binary or higher order systems. The Class 0 protostellar system IRAS16293-2422 contains two protostars, 'A' and 'B', separated by ~600 au and embedded in a single, 10^4 au scale envelope. Their relative evolutionary stages have been debated. We aim to study the relation and interplay between the two protostars A and B at spatial scales of 60 to ~1000 au. We… ▽ More [Abridged] The majority of stars form in binary or higher order systems. The Class 0 protostellar system IRAS16293-2422 contains two protostars, 'A' and 'B', separated by ~600 au and embedded in a single, 10^4 au scale envelope. Their relative evolutionary stages have been debated. We aim to study the relation and interplay between the two protostars A and B at spatial scales of 60 to ~1000 au. We selected molecular gas line transitions of CO, H2CO, HCN, CS, SiO, and CCH from the ALMA-PILS spectral imaging survey (329-363 GHz) and used them as tracers of kinematics, density, and temperature in the IRAS16293-2422 system. The angular resolution of the PILS data set allows us to study these quantities at a resolution of 0.5 arcsec (60 au [..]). Line-of-sight velocity maps of both optically thick and optically thin molecular lines reveal: (i) new manifestations of previously known outflows emanating from protostar A; (ii) a kinematically quiescent bridge of dust and gas spanning between the two protostars, with an inferred density between 4 10^4 and 3 10^7 cm^-3; and (iii) a separate, straight filament seemingly connected to protostar B seen only in CCH, with a flat kinematic signature. Signs of various outflows, all emanating from source A, are evidence of high-density and warmer gas; none of them coincide spatially and kinematically with the bridge. We hypothesize that the bridge arc is a remnant of filamentary substructure in the protostellar envelope material from which protostellar sources A and B have formed. One particular morphological structure appears to be due to outflowing gas impacting the quiescent bridge material. The continuing lack of clear outflow signatures unambiguously associated to protostar B and the vertically extended shape derived for its disk-like structure lead us to conclude that source B may be in an earlier evolutionary stage than source A. △ Less

Submitted 29 March, 2019; originally announced March 2019.

Comments: Accepted for publication in A&A. 27 pages, 24 figures

Journal ref: A&A 626, A93 (2019)

The ALMA-PILS survey: First detection of nitrous acid (HONO) in the interstellar medium

Authors: A. Coutens, N. F. W. Ligterink, J. -C. Loison, V. Wakelam, H. Calcutt, M. N. Drozdovskaya, J. K. Jørgensen, H. S. P. Müller, E. F. van Dishoeck, S. F. Wampfler

Abstract: Nitrogen oxides are thought to play a significant role as a nitrogen reservoir and to potentially participate in the formation of more complex species. Until now, only NO, N$_2$O and HNO have been detected in the interstellar medium. We report the first interstellar detection of nitrous acid (HONO). Twelve lines were identified towards component B of the low-mass protostellar binary IRAS~16293--24… ▽ More Nitrogen oxides are thought to play a significant role as a nitrogen reservoir and to potentially participate in the formation of more complex species. Until now, only NO, N$_2$O and HNO have been detected in the interstellar medium. We report the first interstellar detection of nitrous acid (HONO). Twelve lines were identified towards component B of the low-mass protostellar binary IRAS~16293--2422 with the Atacama Large Millimeter/submillimeter Array, at the position where NO and N$_2$O have previously been seen. A local thermodynamic equilibrium model was used to derive the column density ($\sim$ 9 $\times$ 10$^{14}$ cm$^{-2}$ in a 0.5'' beam) and excitation temperature ($\sim$ 100 K) of this molecule. HNO, NO$_2$, NO$^+$, and HNO$_3$ were also searched for in the data, but not detected. We simulated the HONO formation using an updated version of the chemical code Nautilus and compared the results with the observations. The chemical model is able to reproduce satisfactorily the HONO, N$_2$O, and NO$_2$ abundances, but not the NO, HNO, and NH$_2$OH abundances. This could be due to some thermal desorption mechanisms being destructive and therefore limiting the amount of HNO and NH$_2$OH present in the gas phase. Other options are UV photodestruction of these species in ices or missing reactions potentially relevant at protostellar temperatures. △ Less

Submitted 8 March, 2019; originally announced March 2019.

Comments: Accepted in A&A Letters

Journal ref: A&A 623, L13 (2019)

Deuterated methyl mercaptan (CH$_3$SD): Laboratory rotational spectroscopy and search toward IRAS 16293-2422 B

Authors: Olena Zakharenko, Frank Lewen, Vadim V. Ilyushin, Maria N. Drozdovskaya, Jes K. Jørgensen, Stephan Schlemmer, Holger S. P. Müller

Abstract: Methyl mercaptan (also known as methanethiol), CH3SH, has been found in the warm and dense parts of high -- as well as low -- mass star-forming regions. The aim of the present study is to obtain accurate spectroscopic parameters of the S-deuterated methyl mercaptan CH$_3$SD to facilitate astronomical observations by radio telescope arrays at (sub)millimeter wavelengths. We have measured the rotati… ▽ More Methyl mercaptan (also known as methanethiol), CH3SH, has been found in the warm and dense parts of high -- as well as low -- mass star-forming regions. The aim of the present study is to obtain accurate spectroscopic parameters of the S-deuterated methyl mercaptan CH$_3$SD to facilitate astronomical observations by radio telescope arrays at (sub)millimeter wavelengths. We have measured the rotational spectrum associated with the large-amplitude internal rotation of the methyl group of methyl mercaptan using an isotopically enriched sample in the 150-510 GHz frequency range using the Köln millimeter wave spectrometer. The analysis of the spectra has been performed up to the second excited torsional state. We present modeling results of these data with the RAM36 program. CH$_3$SD was searched for, but not detected, in data from the Atacama Large Millimeter/submillimeter Array (ALMA) Protostellar Interferometric Line Survey (PILS) of the deeply embedded protostar IRAS 16293-2422. The derived upper limit corresponds to a degree of deuteration of at most ~18%. △ Less

Submitted 12 March, 2019; v1 submitted 6 December, 2018; originally announced December 2018.

Comments: Astron. Astrophys., 8 pages; article appeared; CDMS links updated

Journal ref: A&A 621, A114 (2019)

The ALMA-PILS survey: The first detection of doubly-deuterated methyl formate (CHD2OCHO) in the ISM

Authors: S. Manigand, H. Calcutt, J. K. Jørgensen, V. Taquet, H. S. P. Müller, A. Coutens, S. F. Wampfler, N. F. W. Ligterink, M. N. Drozdovskaya, L. E. Kristensen, M. H. D. van der Wiel, T. L. Bourke

Abstract: Studies of deuterated isotopologues of complex organic molecules can provide important constraints on their origin in regions of star formation. In particular, the abundances of deuterated species are very sensitive to the physical conditions in the environment where they form. Due to the low temperatures in regions of star formation, these isotopologues are enhanced to significant levels, making… ▽ More Studies of deuterated isotopologues of complex organic molecules can provide important constraints on their origin in regions of star formation. In particular, the abundances of deuterated species are very sensitive to the physical conditions in the environment where they form. Due to the low temperatures in regions of star formation, these isotopologues are enhanced to significant levels, making detections of multiply-deuterated species possible. However, for complex organic species, only the multiply-deuterated variants of methanol and methyl cyanide have been reported so far. The aim of this paper is to initiate the characterisation of multiply-deuterated variants of complex organic species with the first detection of doubly-deuterated methyl formate, CHD2OCHO. We use ALMA observations from the Protostellar Interferometric Line Survey (PILS) of the protostellar binary IRAS 16293-2422, in the spectral range of 329.1 GHz to 362.9 GHz. We report the first detection of doubly-deuterated methyl formate CHD2OCHO in the ISM. The D/H ratio of CHD2OCHO is found to be 2-3 times higher than the D/H ratio of CH2DOCHO for both sources, similar to the results for formaldehyde from the same dataset. The observations are compared to a gas-grain chemical network coupled to a dynamical physical model, tracing the evolution of a molecular cloud until the end of the Class 0 protostellar stage. The overall D/H ratio enhancements found in the observations are of the same order of magnitude as the predictions from the model for the early stages of Class 0 protostars. However, the higher D/H ratio of CHD2OCHO compared to the D/H ratio of CH2DOCHO is still not predicted by the model. This suggests that a mechanism is enhancing the D/H ratio of singly- and doubly-deuterated methyl formate that is not in the model, e.g. mechanisms for H-D substitutions. △ Less

Submitted 22 November, 2018; originally announced November 2018.

Comments: Accepted for publication in A&A

Journal ref: A&A 623, A69 (2019)

The ALMA-PILS survey: Isotopic composition of oxygen-containing complex organic molecules toward IRAS 16293-2422B

Authors: J. K. Jørgensen, H. S. P. Müller, H. Calcutt, A. Coutens, M. N. Drozdovskaya, K. I. Öberg, M. V. Persson, V. Taquet, E. F. van Dishoeck, S. F. Wampfler

Abstract: This paper presents a systematic survey of the deuterated and 13C isotopologues of a variety of oxygen-bearing complex organic molecules on Solar System scales toward the protostar IRAS 16293-2422B. We use the data from an unbiased molecular line survey between 329 and 363 GHz from the Atacama Large Millimeter/submillimeter Array (ALMA). The observations probe scales of 60 AU where most of the org… ▽ More This paper presents a systematic survey of the deuterated and 13C isotopologues of a variety of oxygen-bearing complex organic molecules on Solar System scales toward the protostar IRAS 16293-2422B. We use the data from an unbiased molecular line survey between 329 and 363 GHz from the Atacama Large Millimeter/submillimeter Array (ALMA). The observations probe scales of 60 AU where most of the organic molecules have sublimated off dust grains and are present in the gas-phase. The complex organic molecules can be divided into two groups with one group, the simpler species, showing a D/H ratio of approximately 2% and the other, the more complex species, D/H ratios of 4-8%. This division may reflect the formation time of each species in the ices before or during warm-up/infall of material through the protostellar envelope. No significant differences are seen in the deuteration of different functional groups for individual species, possibly a result of the short time-scale for infall through the innermost warm regions where exchange reactions between different species may be taking place. The species show differences in excitation temperatures between 125 K and 300 K. This likely reflects the binding energies/sublimation temperatures of the individual species, in good agreement to what has previously been found for high-mass sources. For dimethyl ether, the 12C/13C ratio is found to be lower by up to a factor of 2 compared to typical ISM values similar to what has previously been inferred for glycolaldehyde. The results point to the importance of ice surface chemistry for the formation of these complex organic molecules at different stages in the evolution of embedded protostars and demonstrate the use of accurate isotope measurements for understanding the history of individual species. △ Less

Submitted 27 August, 2018; originally announced August 2018.

Comments: Accepted for publication in A&A. Abstract abridged

Journal ref: A&A 620, A170 (2018)

The ALMA-PILS survey: Stringent limits on small amines and nitrogen-oxides towards IRAS 16293-2422B

Authors: N. F. W. Ligterink, H. Calcutt, A. Coutens, L. E. Kristensen, T. L. Bourke, M. N. Drozdovskaya, H. S. P. Müller, S. F. Wampfler, M. H. D. van der Wiel, E. F. van Dishoeck, J. K. Jørgensen

Abstract: Hydroxylamine (NH$_{2}$OH) and methylamine (CH$_{3}$NH$_{2}$) have both been suggested as precursors to the formation of amino acids and are therefore of interest to prebiotic chemistry. Their presence in interstellar space and formation mechanisms, however, are not well established. We aim to detect both amines and their potential precursor molecules NO, N$_{2}$O and CH$_{2}$NH towards the low-ma… ▽ More Hydroxylamine (NH$_{2}$OH) and methylamine (CH$_{3}$NH$_{2}$) have both been suggested as precursors to the formation of amino acids and are therefore of interest to prebiotic chemistry. Their presence in interstellar space and formation mechanisms, however, are not well established. We aim to detect both amines and their potential precursor molecules NO, N$_{2}$O and CH$_{2}$NH towards the low-mass protostellar binary IRAS 16293--2422, in order to investigate their presence and constrain their interstellar formation mechanisms around a young Sun-like protostar. ALMA observations from the unbiased, high angular resolution and sensitivity Protostellar Interferometric Line Survey (PILS) are used. Spectral transitions of the molecules under investigation are searched for with the CASSIS line analysis software. CH$_2$NH and N$_{2}$O are detected for the first time towards a low-mass source, the latter molecule through confirmation with the single-dish TIMASSS survey. NO is also detected. CH$_{3}$NH$_{2}$ and NH$_{2}$OH are not detected and stringent upper limit column densities are determined. The non-detection of CH$_{3}$NH$_{2}$ and NH$_{2}$OH limits the importance of formation routes to amino acids involving these species. The detection of CH$_{2}$NH makes amino acid formation routes starting from this molecule plausible. The low abundances of CH$_2$NH and CH$_{3}$NH$_{2}$ compared to Sgr B2 indicate that different physical conditions influence their formation in low- and high-mass sources. △ Less

Submitted 1 August, 2018; originally announced August 2018.

Comments: Accepted for publication in A&A

Journal ref: A&A 619, A28 (2018)

The ALMA-PILS survey: first detection of methyl isocyanide (CH$_3$NC) in a solar-type protostar

Authors: H. Calcutt, M. R. Fiechter, E. R. Willis, H. S. P. Müller, R. T. Garrod, J. K. Jørgensen, S. F. Wampfler, T. L. Bourke, A. Coutens, M. N. Drozdovskaya, N. F. W. Ligterink, L. E. Kristensen

Abstract: Methyl isocyanide (CH$_3$NC) is the isocyanide with the largest number of atoms confirmed in the interstellar medium (ISM), but it is not an abundant molecule, having only been detected towards a handful of objects. Conversely, its isomer, methyl cyanide (CH$_3$CN), is one of the most abundant complex organic molecules detected in the ISM, with detections in a variety of low- and high-mass sources… ▽ More Methyl isocyanide (CH$_3$NC) is the isocyanide with the largest number of atoms confirmed in the interstellar medium (ISM), but it is not an abundant molecule, having only been detected towards a handful of objects. Conversely, its isomer, methyl cyanide (CH$_3$CN), is one of the most abundant complex organic molecules detected in the ISM, with detections in a variety of low- and high-mass sources. We use ALMA observations from the Protostellar Interferometric Line Survey (PILS) to search for methyl isocyanide and compare its abundance with that of its isomer methyl cyanide. We use a new line catalogue from the Cologne Database for Molecular Spectroscopy (CDMS) to identify methyl isocyanide lines. We also model the chemistry with an updated version of the three-phase chemical kinetics model {\em MAGICKAL}, presenting the first chemical modelling of methyl isocyanide to date. We detect methyl isocyanide for the first time in a solar-type protostar, IRAS 16293$-$2422 B, and present upper limits for its companion protostar, IRAS 16293$-$2422 A. Methyl isocyanide is found to be at least 20 times more abundant in source B compared to source A, with a CH$_3$CN/CH$_3$NC abundance ratio of 200 in IRAS 16293--2422 B and >5517 in IRAS 16293$-$2422 A. We also present the results of a chemical model of methyl isocyanide chemistry in both sources, and discuss the implications on methyl isocyanide formation mechanisms and the relative evolutionary stages of both sources. △ Less

Submitted 8 July, 2018; originally announced July 2018.

Comments: Accepted to A&A on 5th July 2018

Journal ref: A&A 617, A95 (2018)

Linking interstellar and cometary O$_2$: a deep search for $^{16}$O$^{18}$O in the solar-type protostar IRAS 16293--2422

Authors: Vianney Taquet, Ewine F. van Dishoeck, Matthew Swayne, Daniel Harsono, Jes K. Jørgensen, Luke Maud, Niels F. W. Ligterink, Holger S. P. Müller, Claudio Codella, Kathrin Altwegg, Andre Bieler, Audrey Coutens, Maria N. Drozdovskaya, Kenji Furuya, Magnus V. Persson, Merel L. R. van 't Hoff, Catherine Walsh, Suzanne F. Wampfler

Abstract: Recent measurements carried out at comet 67P/C-G with the ${\it Rosetta}$ probe revealed that molecular oxygen, O$_2$, is the fourth most abundant molecule in comets. Models show that O$_2$ is likely of primordial nature, coming from the interstellar cloud from which our Solar System was formed. However, gaseous O$_2$ is an elusive molecule in the interstellar medium with only one detection toward… ▽ More Recent measurements carried out at comet 67P/C-G with the ${\it Rosetta}$ probe revealed that molecular oxygen, O$_2$, is the fourth most abundant molecule in comets. Models show that O$_2$ is likely of primordial nature, coming from the interstellar cloud from which our Solar System was formed. However, gaseous O$_2$ is an elusive molecule in the interstellar medium with only one detection towards quiescent molecular clouds, in the $ρ$ Oph A core. We perform a deep search for molecular oxygen, through the $2_1 - 0_1$ rotational transition at 234 GHz of its $^{16}$O$^{18}$O isotopologue, towards the warm compact gas surrounding the nearby Class 0 protostar IRAS 16293--2422 B with the ALMA interferometer. The targeted $^{16}$O$^{18}$O transition is surrounded by two brighter transitions at $\pm 1$ km s$^{-1}$ relative to the expected $^{16}$O$^{18}$O transition frequency. After subtraction of these two transitions, residual emission at a 3$σ$ level remains, but with a velocity offset of $0.3 - 0.5$ km s$^{-1}$ relative to the source velocity, rendering the detection "tentative". We derive the O$_2$ column density for two excitation temperatures $T_{\rm ex}$ of 125 and 300 K, as indicated by other molecules, in order to compare the O$_2$ abundance between IRAS16293 and comet 67P/C-G. Assuming that $^{16}$O$^{18}$O is not detected and using methanol CH$_3$OH as a reference species, we obtain a [O$_2$]/[CH$_3$OH] abundance ratio lower than $2-5$, depending on the assumed $T_{\rm ex}$, a three to four times lower abundance than the [O$_2$]/[CH$_3$OH] ratio of $5-15$ found in comet 67P/C-G. Such a low O$_2$ abundance could be explained by the lower temperature of the dense cloud precursor of IRAS16293 with respect to the one at the origin of our Solar System that prevented an efficient formation of O$_2$ in interstellar ices. △ Less

Submitted 31 May, 2018; originally announced May 2018.

Comments: 12 pages, 3 figures, accepted for publication in Astronomy & Astrophysics

Tracing the cold and warm physico-chemical structure of deeply embedded protostars: IRAS 16293-2422 versus VLA 1623-2417

Authors: Nadia M. Murillo, Ewine F. van Dishoeck, Matthijs H. D. van der Wiel, Jes K. Jørgensen, Maria N. Drozdovskaya, Hannah Calcutt, Daniel Harsono

Abstract: Much attention has been placed on the dust distribution in protostellar envelopes, but there are still many unanswered questions regarding the structure of the gas. We aim to start identifying the factors that determine the chemical structure of protostellar regions, by studying and comparing low-mass embedded systems in key molecular tracers. The cold and warm chemical structures of two embedded… ▽ More Much attention has been placed on the dust distribution in protostellar envelopes, but there are still many unanswered questions regarding the structure of the gas. We aim to start identifying the factors that determine the chemical structure of protostellar regions, by studying and comparing low-mass embedded systems in key molecular tracers. The cold and warm chemical structures of two embedded Class 0 systems, IRAS16293 and VLA1623 are characterized through interferometric observations. DCO+, N2H+ and N2D+ are used to trace the spatial distribution and physics of the cold regions of the envelope, while c-C3H2 and C2H from models of the chemistry are expected to trace the warm (UV-irradiated) regions. Both sources show a number of striking similarities and differences. DCO+ consistently traces the cold material at the disk-envelope interface, where gas and dust temperatures are lowered due to disk shadowing. N2H+ and N2D+, also tracing cold gas, show low abundances towards VLA1623, but for IRAS16293, the distribution of N2D+ is consistent with the same chemical models that reproduce DCO+. c-C3H2 and C2H show different spatial distributions for the two systems. For IRAS16293, c-C3H2 traces the outflow cavity wall, while C2H is found in the envelope material but not the outflow cavity wall. In contrast, toward VLA1623 both molecules trace the outflow cavity wall. Finally, hot core molecules are abundantly observed toward IRAS16293 but not toward VLA1623. We identify temperature as one of the key factors in determining the chemical structure of protostars as seen in gaseous molecules. More luminous protostars, such as IRAS16293, will have chemical complexity out to larger distances than colder protostars, such as VLA1623. Additionally, disks in the embedded phase have a crucial role in controlling both the gas and dust temperature of the envelope, and consequently the chemical structure. △ Less

Submitted 14 May, 2018; originally announced May 2018.

Comments: 24 pages, 15 figures, accepted to A&A

Journal ref: A&A 617, A120 (2018)

The ALMA-PILS survey: The sulphur connection between protostars and comets: IRAS 16293-2422 B and 67P/Churyumov-Gerasimenko

Authors: Maria N. Drozdovskaya, Ewine F. van Dishoeck, Jes K. Jørgensen, Ursina Calmonte, Matthijs H. D. van der Wiel, Audrey Coutens, Hannah Calcutt, Holger S. P. Müller, Per Bjerkeli, Magnus V. Persson, Susanne F. Wampfler, Kathrin Altwegg

Abstract: The evolutionary past of our Solar System can be pieced together by comparing analogous low-mass protostars with remnants of our Protosolar Nebula - comets. Sulphur-bearing molecules may be unique tracers of the joint evolution of the volatile and refractory components. ALMA Band 7 data from the large unbiased Protostellar Interferometric Line Survey (PILS) are used to search for S-bearing molecul… ▽ More The evolutionary past of our Solar System can be pieced together by comparing analogous low-mass protostars with remnants of our Protosolar Nebula - comets. Sulphur-bearing molecules may be unique tracers of the joint evolution of the volatile and refractory components. ALMA Band 7 data from the large unbiased Protostellar Interferometric Line Survey (PILS) are used to search for S-bearing molecules in the outer disc-like structure, 60 au from IRAS 16293-2422 B, and are compared with data on 67P/C-G stemming from the ROSINA instrument aboard Rosetta. Species such as SO$_{2}$, SO, OCS, CS, H$_{2}$CS, H$_{2}$S and CH$_{3}$SH are detected via at least one of their isotopologues towards IRAS 16293-2422 B. The search reveals a first-time detection of OC$^{33}$S towards this source and a tentative first-time detection of C$^{36}$S towards a low-mass protostar. The data show that IRAS 16293-2422 B contains much more OCS than H$_{2}$S in comparison to 67P/C-G; meanwhile, the SO/SO$_{2}$ ratio is in close agreement between the two targets. IRAS 16293-2422 B has a CH$_{3}$SH/H$_{2}$CS ratio in range of that of our Solar System (differences by a factor of 0.7-5.3). It is suggested that the levels of UV radiation during the initial collapse of the systems may have varied and have potentially been higher for IRAS 16293-2422 B due to its binary nature; thereby, converting more H$_{2}$S into OCS. It remains to be conclusively tested if this also promotes the formation of S-bearing complex organics. Elevated UV levels of IRAS 16293-2422 B and a warmer birth cloud of our Solar System may jointly explain the variations between the two low-mass systems. △ Less

Submitted 21 February, 2018; v1 submitted 8 February, 2018; originally announced February 2018.

Comments: Accepted for publication in MNRAS; 47 pages, 43 figures, 6 tables

Journal ref: Mon. Not. R. Astron. Soc. 476 (2018) 4949-4964

The chemical connection between 67P/C-G and IRAS 16293-2422

Authors: Maria Nikolayevna Drozdovskaya, Ewine F. van Dishoeck, Martin Rubin, Jes Kristian Jørgensen, Kathrin Altwegg

Abstract: The chemical evolution of a star- and planet-forming system begins in the prestellar phase and proceeds across the subsequent evolutionary phases. The chemical trail from cores to protoplanetary disks to planetary embryos can be studied by comparing distant young protostars and comets in our Solar System. One particularly chemically rich system that is thought to be analogous to our own is the low… ▽ More The chemical evolution of a star- and planet-forming system begins in the prestellar phase and proceeds across the subsequent evolutionary phases. The chemical trail from cores to protoplanetary disks to planetary embryos can be studied by comparing distant young protostars and comets in our Solar System. One particularly chemically rich system that is thought to be analogous to our own is the low-mass IRAS 16293-2422. ALMA-PILS observations have made the study of chemistry on the disk scales (< 100 AU) of this system possible. Under the assumption that comets are pristine tracers of the outer parts of the innate protosolar disk, it is possible to compare the composition of our infant Solar System to that of IRAS 16293-2422. The Rosetta mission has yielded a wealth of unique in situ measurements on comet 67P/C-G, making it the best probe to date. Herein, the initial comparisons in terms of the chemical composition and isotopic ratios are summarized. Much work is still to be carried out in the future as the analysis of both of these data sets is still ongoing. △ Less

Submitted 8 February, 2018; originally announced February 2018.

Comments: To appear in "Astrochemistry VII -- Through the Cosmos from Galaxies to Planets", proceedings of the IAU Symposium No. 332, 2017, Puerto Varas, Chile. M. Cunningham, T. Millar and Y. Aikawa, eds. (6 pages)

First detection of cyanamide (NH$_2$CN) towards solar-type protostars

Authors: A. Coutens, E. R. Willis, R. T. Garrod, H. S. P. Müller, T. L. Bourke, H. Calcutt, M. N. Drozdovskaya, J. K. Jørgensen, N. F. W. Ligterink, M. V. Persson, G. Stéphan, M. H. D. van der Wiel, E. F. van Dishoeck, S. F. Wampfler

Abstract: Searches for the prebiotically-relevant cyanamide (NH$_2$CN) towards solar-type protostars have not been reported in the literature. We here present the first detection of this species in the warm gas surrounding two solar-type protostars, using data from the Atacama Large Millimeter/Submillimeter Array Protostellar Interferometric Line Survey (PILS) of IRAS 16293-2422 B and observations from the… ▽ More Searches for the prebiotically-relevant cyanamide (NH$_2$CN) towards solar-type protostars have not been reported in the literature. We here present the first detection of this species in the warm gas surrounding two solar-type protostars, using data from the Atacama Large Millimeter/Submillimeter Array Protostellar Interferometric Line Survey (PILS) of IRAS 16293-2422 B and observations from the IRAM Plateau de Bure Interferometer of NGC1333 IRAS2A. We furthermore detect the deuterated and $^{13}$C isotopologues of NH$_2$CN towards IRAS 16293-2422 B. This is the first detection of NHDCN in the interstellar medium. Based on a local thermodynamic equilibrium analysis, we find that the deuteration of cyanamide ($\sim$ 1.7%) is similar to that of formamide (NH$_2$CHO), which may suggest that these two molecules share NH$_2$ as a common precursor. The NH$_2$CN/NH$_2$CHO abundance ratio is about 0.2 for IRAS 16293-2422 B and 0.02 for IRAS2A, which is comparable to the range of values found for Sgr B2. We explored the possible formation of NH$_2$CN on grains through the NH$_2$ + CN reaction using the chemical model MAGICKAL. Grain-surface chemistry appears capable of reproducing the gas-phase abundance of NH$_2$CN with the correct choice of physical parameters. △ Less

Submitted 27 December, 2017; originally announced December 2017.

Comments: Accepted in A&A

Journal ref: A&A 612, A107 (2018)

The ALMA-PILS survey: 3D modeling of the envelope, disks and dust filament of IRAS 16293-2422

Authors: S. K. Jacobsen, J. K. Jørgensen, M. H. D. van der Wiel, H. Calcutt, T. L. Bourke, C. Brinch, A. Coutens, M. N. Drozdovskaya, L. E. Kristensen, H. S. P. Müller, S. F. Wampfler

Abstract: Context. The Class 0 protostellar binary IRAS 16293-2422 is an interesting target for (sub)millimeter observations due to, both, the rich chemistry toward the two main components of the binary and its complex morphology. Its proximity to Earth allows the study of its physical and chemical structure on solar system scales using high angular resolution observations. Such data reveal a complex morpho… ▽ More Context. The Class 0 protostellar binary IRAS 16293-2422 is an interesting target for (sub)millimeter observations due to, both, the rich chemistry toward the two main components of the binary and its complex morphology. Its proximity to Earth allows the study of its physical and chemical structure on solar system scales using high angular resolution observations. Such data reveal a complex morphology that cannot be accounted for in traditional, spherical 1D models of the envelope. Aims. The purpose of this paper is to study the environment of the two components of the binary through 3D radiative transfer modeling and to compare with data from the Atacama Large Millimeter/submillimeter Array. Such comparisons can be used to constrain the protoplanetary disk structures, the luminosities of the two components of the binary and the chemistry of simple species. Methods. We present 13CO, C17O and C18O J=3-2 observations from the ALMA Protostellar Interferometric Line Survey (PILS), together with a qualitative study of the dust and gas density distribution of IRAS 16293-2422. A 3D dust and gas model including disks and a dust filament between the two protostars is constructed which qualitatively reproduces the dust continuum and gas line emission. Results and conclusions. Radiative transfer modeling of source A and B, with the density solution of an infalling, rotating collapse or a protoplanetary disk model, can match the constraints for the disk-like emission around source A and B from the observed dust continuum and CO isotopologue gas emission. If a protoplanetary disk model is used around source B, it has to have an unusually high scale-height in order to reach the dust continuum peak emission value, while fulfilling the other observational constraints. Our 3D model requires source A to be much more luminous than source B; LA ~ 18 $L_\odot$ and LB ~ 3 $L_\odot$. △ Less

Submitted 21 December, 2017; v1 submitted 19 December, 2017; originally announced December 2017.

Comments: Accepted for publication in Astronomy and Astrophysics, 20 pages, 12 figures

Journal ref: A&A 612, A72 (2018)

Outflows, infall and evolution of a sample of embedded low-mass protostars. The William Herschel Line Legacy (WILL) survey

Authors: J. C. Mottram, E. F. van Dishoeck, L. E. Kristensen, A. Karska, I. San José-García, S. Khanna, G. J. Herczeg, Ph. Andr, S. Bontemps, S. Cabrit, M. T. Carney, M. N. Drozdovskaya, M. M. Dunham, N. J. Evans, D. Fedele, J. D. Green, D. Harsono, D. Johnstone, J. K. Jørgensen, V. Könyves, B. Nisini, M. V. Persson, M. Tafalla, R. Visser, U. A. Yıldız

Abstract: [Abridged] We present spectroscopic observations in H$_{2}$O, CO and related species with \textit{Herschel} HIFI and PACS, as well as ground-based follow-up with the JCMT and APEX in CO, HCO$^{+}$ and isotopologues, of a sample of 49 nearby ($d<$500\,pc) candidate protostars. These data are used to study the outflow and envelope properties of these sources. We also compile their continuum SEDs in… ▽ More [Abridged] We present spectroscopic observations in H$_{2}$O, CO and related species with \textit{Herschel} HIFI and PACS, as well as ground-based follow-up with the JCMT and APEX in CO, HCO$^{+}$ and isotopologues, of a sample of 49 nearby ($d<$500\,pc) candidate protostars. These data are used to study the outflow and envelope properties of these sources. We also compile their continuum SEDs in order to constrain their physical properties. Water emission is dominated by shocks associated with the outflow, rather than the cooler, slower entrained outflowing gas probed by ground-based CO observations. These shocks become less energetic as sources evolve from Class 0 to Class I. The fraction of mass in the outflow relative to the total envelope (i.e. $M_{\mathrm{out}}/M_{\mathrm{env}}$) remains broadly constant between Class 0 and I. The median value ($\sim$1$\%$) is consistent with a core to star formation efficiency on the order of 50$\%$ and an outflow duty cycle on the order of 5$\%$. Entrainment efficiency, as probed by $F_{\mathrm{CO}}/\dot{M}_{\mathrm{acc}}$, is also invariant with source properties and evolutionary stage. The median value (6.3\kms{}) suggests an entrainment efficiency of between 30 and 60$\%$ if the wind is launched at $\sim$1AU. $L$[O\,{\sc i}] is strongly correlated with $L_{\mathrm{bol}}$ but not with $M_{\mathrm{env}}$, while low-$J$ CO is more closely correlated with the latter than the former. This suggests that [O\,{\sc i}] traces the present-day accretion activity while CO traces time-averaged accretion over the dynamical timescale of the outflow. $L$[O\,{\sc i}] does not vary from Class 0 to Class I, unlike CO and H$_{2}$O. This is likely due to the ratio of atomic to molecular gas in the wind increasing as the source evolves, balancing out the decrease in mass accretion rate. Infall signatures are detected in HCO$^{+}$ and H$_{2}$O in a few sources. △ Less

Submitted 19 January, 2017; v1 submitted 17 January, 2017; originally announced January 2017.

Comments: Accepted to A&A, version after language editor corrections. 47 pages, 18 figures, 15 tables. The figures have been converted to pdf due to file-size considerations, which may lead to some degradation

The ALMA-PILS survey: First detections of ethylene oxide, acetone and propanal toward the low-mass protostar IRAS 16293-2422

Authors: J. M. Lykke, A. Coutens, J. K. Jørgensen, M. H. D. van der Wiel, R. T. Garrod, H. S. P. Müller, P. Bjerkeli, T. L. Bourke, H. Calcutt, M. N. Drozdovskaya, C. Favre, E. C. Fayolle, S. K. Jacobsen, K. I. Öberg, M. V. Persson, E. F. van Dishoeck, S. F. Wampfler

Abstract: One of the open questions in astrochemistry is how complex organic and prebiotic molecules are formed. Aims. Our aim is to start the process of compiling an inventory of oxygen-bearing complex organic molecules toward the solar-type Class 0 protostellar binary IRAS16293-2422 from an unbiased spectral survey with ALMA (PILS). Here we focus on the new detections of ethylene oxide (c-C$_2$H$_4$O), ac… ▽ More One of the open questions in astrochemistry is how complex organic and prebiotic molecules are formed. Aims. Our aim is to start the process of compiling an inventory of oxygen-bearing complex organic molecules toward the solar-type Class 0 protostellar binary IRAS16293-2422 from an unbiased spectral survey with ALMA (PILS). Here we focus on the new detections of ethylene oxide (c-C$_2$H$_4$O), acetone (CH$_3$COCH$_3$), and propanal (C$_2$H$_5$CHO). Methods. With ALMA, we surveyed the spectral range from 329 to 363 GHz at 0.5$"$ (60 AU diameter) resolution. Using a simple model for the molecular emission in LTE, the excitation temperatures and column densities of each species were constrained. Results. We successfully detect propanal (44 lines), ethylene oxide (20 lines) and acetone (186 lines) toward one component of the protostellar binary, IRAS16293B. The high resolution maps demonstrate that the emission for all investigated species originates from the compact central region close to the protostar. This, along with a derived common excitation temperature of $\sim$ 125 K, is consistent with a coexistence of these molecules in the same gas. Conclusions. The observations mark the first detections of acetone, propanal and ethylene oxide toward a low-mass protostar. The relative abundance ratios of the two sets of isomers (CH$_3$COCH$_3$/C$_2$H$_5$CHO $\sim$ 8 and CH$_3$CHO/c-C$_2$H$_4$O $\sim$ 12) are comparable to previous observations toward high-mass protostars. The majority of observed abundance ratios from these results as well as those measured toward high-mass protostars are up to an order of magnitude above the predictions from chemical models. This may reflect either missing reactions or uncertain rates in the chemical networks. The physical conditions, such as temperatures or densities, used in the models, may not be applicable to solar-type protostars either. △ Less

Submitted 22 November, 2016; originally announced November 2016.

Comments: Accepted in A&A

Journal ref: Astron. Astrophys. 597 (2017), Art. No. A53

Water delivery from cores to disks: deuteration as a probe of the prestellar inheritance of H2O

Authors: K. Furuya, M. N. Drozdovskaya, R. Visser, E. F. van Dishoeck, C. Walsh, D. Harsono, U. Hincelin, V. Taquet

Abstract: We investigate the delivery of regular and deuterated forms of water from prestellar cores to circumstellar disks. We adopt a semi-analytical axisymmetric two-dimensional collapsing core model with post-processing gas-ice astrochemical simulations, in which a layered ice structure is considered. The physical and chemical evolutions are followed until the end of the main accretion phase. When mass… ▽ More We investigate the delivery of regular and deuterated forms of water from prestellar cores to circumstellar disks. We adopt a semi-analytical axisymmetric two-dimensional collapsing core model with post-processing gas-ice astrochemical simulations, in which a layered ice structure is considered. The physical and chemical evolutions are followed until the end of the main accretion phase. When mass averaged over the whole disk, a forming disk has a similar H2O abundance and HDO/H2O abundance ratio as their precollapse values (within a factor of 2), regardless of time in our models. Consistent with previous studies, our models suggest that interstellar water ice is delivered to forming disks without significant alteration. On the other hand, the local vertically averaged H2O ice abundance and HDO/H2O ice ratio can differ more, by up to a factor of several, depending on time and distance from a central star. Key parameters for the local variations are the fluence of stellar UV photons en route into the disk and the ice layered structure, the latter of which is mostly established in the prestellar stages. We also find that even if interstellar water ice is destroyed by stellar UV and (partly) reformed prior to disk entry, the HDO/H2O ratio in reformed water ice is similar to the original value. This finding indicates that some caution is needed in discussions on the prestellar inheritance of H2O based on comparisons between the observationally derived HDO/H2O ratio in clouds/cores and that in disks/comets. Alternatively, we propose that the ratio of D2O/HDO to HDO/H2O better probes the prestellar inheritance of H2O. It is also found that icy organics are more enriched in deuterium than water ice in forming disks. The differential deuterium fractionation in water and organics is inherited from the prestellar stages. △ Less

Submitted 24 October, 2016; originally announced October 2016.

Comments: 20 pages, 15 figures, accepted by A&A

Journal ref: A&A 599, A40 (2017)

The ALMA Protostellar Interferometric Line Survey (PILS): First results from an unbiased submillimeter wavelength line survey of the Class 0 protostellar binary IRAS 16293-2422 with ALMA

Authors: J. K. Jørgensen, M. H. D. van der Wiel, A. Coutens, J. M. Lykke, H. S. P. Müller, E. F. van Dishoeck, H. Calcutt, P. Bjerkeli, T. L. Bourke, M. N. Drozdovskaya, C. Favre, E. C. Fayolle, R. T. Garrod, S. K. Jacobsen, K. I. Öberg, M. V. Persson, S. F. Wampfler

Abstract: The inner regions of the envelopes surrounding young protostars are characterised by a complex chemistry, with prebiotic molecules present on the scales where protoplanetary disks eventually may form. This paper introduces a systematic survey, "Protostellar Interferometric Line Survey (PILS)" of the Class 0 protostellar binary IRAS 16293-2422 using the Atacama Large Millimeter/submillimeter Array… ▽ More The inner regions of the envelopes surrounding young protostars are characterised by a complex chemistry, with prebiotic molecules present on the scales where protoplanetary disks eventually may form. This paper introduces a systematic survey, "Protostellar Interferometric Line Survey (PILS)" of the Class 0 protostellar binary IRAS 16293-2422 using the Atacama Large Millimeter/submillimeter Array (ALMA). The survey covers the full frequency range from 329 to 363 GHz (0.8 mm) with additional targeted observations at 3.0 and 1.3 mm. More than 10,000 features are detected toward one component in the protostellar binary. Glycolaldehyde, its isomers, methyl formate and acetic acid, and its reduced alcohol, ethylene glycol, are clearly detected. For ethylene glycol both lowest state conformers, aGg' and gGg', are detected, the latter for the first time in the ISM. The abundance of glycolaldehyde is comparable to or slightly larger than that of ethylene glycol. In comparison to the Galactic Center, these two species are over-abundant relative to methanol, possibly an indication of formation at low temperatures in CO-rich ices. Both 13C and deuterated isotopologues of glycolaldehyde are detected, also for the first time ever in the ISM. For the deuterated species, a D/H ratio of approximately 5% is found with no differences between the deuteration in the different functional groups of glycolaldehyde. Measurements of the 13C-species lead to a 12C:13C ratio of approximately 30, lower than the typical ISM value. This low ratio may reflect an enhancement of 13CO in the ice due to either ion-molecule reactions in the gas before freeze-out or differences in the temperatures where 12CO and 13CO ices sublimate. The results reinforce the importance of low-temperature grain surface chemistry for the formation of prebiotic molecules seen here in the gas after sublimation of the entire ice mantle. △ Less

Submitted 29 July, 2016; originally announced July 2016.

Comments: Accepted for publication in A&A. Abstract abridged

Journal ref: A&A 595, A117 (2016)

Cometary ices in forming protoplanetary disc midplanes

Authors: Maria N. Drozdovskaya, Catherine Walsh, Ewine F. van Dishoeck, Kenji Furuya, Ulysse Marboeuf, Amaury Thiabaud, Daniel Harsono, Ruud Visser

Abstract: Low-mass protostars are the extrasolar analogues of the natal Solar System. Sophisticated physicochemical models are used to simulate the formation of two protoplanetary discs from the initial prestellar phase, one dominated by viscous spreading and the other by pure infall. The results show that the volatile prestellar fingerprint is modified by the chemistry en route into the disc. This holds re… ▽ More Low-mass protostars are the extrasolar analogues of the natal Solar System. Sophisticated physicochemical models are used to simulate the formation of two protoplanetary discs from the initial prestellar phase, one dominated by viscous spreading and the other by pure infall. The results show that the volatile prestellar fingerprint is modified by the chemistry en route into the disc. This holds relatively independent of initial abundances and chemical parameters: physical conditions are more important. The amount of CO2 increases via the grain-surface reaction of OH with CO, which is enhanced by photodissociation of H2O ice. Complex organic molecules are produced during transport through the envelope at the expense of CH3OH ice. Their abundances can be comparable to that of methanol ice (few % of water ice) at large disc radii (R > 30 AU). Current Class II disc models may be underestimating the complex organic content. Planet population synthesis models may underestimate the amount of CO2 and overestimate CH3OH ices in planetesimals by disregarding chemical processing between the cloud and disc phases. The overall C/O and C/N ratios differ between the gas and solid phases. The two ice ratios show little variation beyond the inner 10 AU and both are nearly solar in the case of pure infall, but both are sub-solar when viscous spreading dominates. Chemistry in the protostellar envelope en route to the protoplanetary disc sets the initial volatile and prebiotically-significant content of icy planetesimals and cometary bodies. Comets are thus potentially reflecting the provenances of the midplane ices in the Solar Nebula. △ Less

Submitted 26 July, 2016; originally announced July 2016.

Comments: Accepted for publication in MNRAS; 19 pages, 7 figures, 6 tables

The ALMA-PILS survey: First detections of deuterated formamide and deuterated isocyanic acid in the interstellar medium

Authors: A. Coutens, J. K. Jørgensen, M. H. D. van der Wiel, H. S. P. Müller, J. M. Lykke, P. Bjerkeli, T. L. Bourke, H. Calcutt, M. N. Drozdovskaya, C. Favre, E. C. Fayolle, R. T. Garrod, S. K. Jacobsen, N. F. W. Ligterink, K. I. Öberg, M. V. Persson, E. F. van Dishoeck, S. F. Wampfler

Abstract: Formamide (NH$_2$CHO) has previously been detected in several star-forming regions and is thought to be a precursor for different prebiotic molecules. Its formation mechanism is still debated, however. Observations of formamide, related species and their isopotologues may provide useful clues to the chemical pathways leading to their formation. The Protostellar Interferometric Line Survey (PILS) r… ▽ More Formamide (NH$_2$CHO) has previously been detected in several star-forming regions and is thought to be a precursor for different prebiotic molecules. Its formation mechanism is still debated, however. Observations of formamide, related species and their isopotologues may provide useful clues to the chemical pathways leading to their formation. The Protostellar Interferometric Line Survey (PILS) represents an unbiased high angular resolution and sensitivity spectral survey of the low-mass protostellar binary IRAS 16293-2422 with the Atacama Large Millimeter/submillimeter Array (ALMA). We detect for the first time the three singly deuterated forms of NH$_2$CHO (NH$_2$CDO, cis- and trans-NHDCHO) as well as DNCO towards the component B of this binary source. The images reveal that the different isotopologues are all present in the same region. Based on observations of the $^{13}$C isotopologues of formamide and a standard $^{12}$C/$^{13}$C ratio, the deuterium fractionation is found to be similar for the three different forms with a value of about 2%. The DNCO/HNCO ratio is also comparable to the D/H ratio of formamide ($\sim$1%). These results are in agreement with the hypothesis that NH$_2$CHO and HNCO are chemically related through grain surface formation. △ Less

Submitted 12 May, 2016; v1 submitted 9 May, 2016; originally announced May 2016.

Comments: Accepted in A&A Letters

Journal ref: A&A 590, L6 (2016)

The complex chemistry of outflow cavity walls exposed: the case of low-mass protostars

Authors: Maria N. Drozdovskaya, Catherine Walsh, Ruud Visser, Daniel Harsono, Ewine F. van Dishoeck

Abstract: Complex organic molecules are ubiquitous companions of young low-mass protostars. Recent observations suggest that their emission stems, not only from the traditional hot corino, but also from offset positions. In this work, 2D physicochemical modelling of an envelope-cavity system is carried out. Wavelength-dependent radiative transfer calculations are performed and a comprehensive gas-grain chem… ▽ More Complex organic molecules are ubiquitous companions of young low-mass protostars. Recent observations suggest that their emission stems, not only from the traditional hot corino, but also from offset positions. In this work, 2D physicochemical modelling of an envelope-cavity system is carried out. Wavelength-dependent radiative transfer calculations are performed and a comprehensive gas-grain chemical network is used to simulate the physical and chemical structure. The morphology of the system delineates three distinct regions: the cavity wall layer with time-dependent and species-variant enhancements; a torus rich in complex organic ices, but not reflected in gas-phase abundances; and the remaining outer envelope abundant in simpler solid and gaseous molecules. Strongly irradiated regions, such as the cavity wall layer, are subject to frequent photodissociation in the solid phase. Subsequent recombination of the photoproducts leads to frequent reactive desorption, causing gas-phase enhancements of several orders of magnitude. This mechanism remains to be quantified with laboratory experiments. Direct photodesorption is found to be relatively inefficient. If radicals are not produced directly in the icy mantle, the formation of complex organics is impeded. For efficiency, a sufficient number of FUV photons needs to penetrate the envelope; and elevated cool dust temperatures need to enable grain-surface radical mobility. As a result, a high stellar luminosity and a sufficiently wide cavity favor chemical complexity. Furthermore within this paradigm, complex organics are demonstrated to have unique lifetimes and be grouped into early (formaldehyde, ketene, methanol, formic acid, methyl formate, acetic acid, glycolaldehyde) and late (acetaldehyde, dimethyl ether, ethanol) species. △ Less

Submitted 21 May, 2015; originally announced May 2015.

Comments: Accepted for publication in MNRAS; 24 pages, 20 figures, 3 tables

Shockingly low water abundances in Herschel / PACS observations of low-mass protostars in Perseus

Authors: A. Karska, L. E. Kristensen, E. F. van Dishoeck, M. N. Drozdovskaya, J. C. Mottram, G. J. Herczeg, S. Bruderer, S. Cabrit, N. J. Evans II, D. Fedele, A. Gusdorf, J. K. Jorgensen, M. J. Kaufman, G. J. Melnick, D. A. Neufeld, B. Nisini, G. Santangelo, M. Tafalla, S. F. Wampfler

Abstract: Protostars interact with their surroundings through jets and winds impacting on the envelope and creating shocks, but the nature of these shocks is still poorly understood. Our aim is to survey far-infrared molecular line emission from a uniform and significant sample of deeply-embedded low-mass young stellar objects in order to characterize shocks and the possible role of ultraviolet radiation in… ▽ More Protostars interact with their surroundings through jets and winds impacting on the envelope and creating shocks, but the nature of these shocks is still poorly understood. Our aim is to survey far-infrared molecular line emission from a uniform and significant sample of deeply-embedded low-mass young stellar objects in order to characterize shocks and the possible role of ultraviolet radiation in the immediate protostellar environment. Herschel/PACS spectral maps of 22 objects in the Perseus molecular cloud were obtained as part of the `William Herschel Line Legacy' survey. Line emission from H$_\mathrm{2}$O, CO, and OH is tested against shock models from the literature. Observed line ratios are remarkably similar and do not show variations with source physical parameters. Observations show good agreement with the shock models when line ratios of the same species are compared. Ratios of various H$_\mathrm{2}$O lines provide a particularly good diagnostic of pre-shock gas densities, $n_\mathrm{H}\sim10^{5}$ cm$^{-3}$, in agreement with typical densities obtained from observations of the post-shock gas. The corresponding shock velocities, obtained from comparison with CO line ratios, are above 20 km\,s$^{-1}$. However, the observations consistently show one-to-two orders of magnitude lower H$_\mathrm{2}$O-to-CO and H$_\mathrm{2}$O-to-OH line ratios than predicted by the existing shock models. The overestimated model H$_\mathrm{2}$O fluxes are most likely caused by an overabundance of H$_\mathrm{2}$O in the models since the excitation is well-reproduced. Illumination of the shocked material by ultraviolet photons produced either in the star-disk system or, more locally, in the shock, would decrease the H$_\mathrm{2}$O abundances and reconcile the models with observations. Detections of hot H$_\mathrm{2}$O and strong OH lines support this scenario. △ Less

Submitted 16 September, 2014; originally announced September 2014.

Comments: 28 pages, 12 figures, accepted to Astronomy & Astrophysics

Journal ref: A&A 572, A9 (2014)