Showing 1–50 of 125 results for author: Kristensen, L E

Linking ice and gas in the Coronet cluster in Corona Australis

Authors: G. Perotti, J. K. Jørgensen, W. R. M. Rocha, A. Plunkett, E. Artur de la Villarmois, L. E. Kristensen, M. Sewiło, P. Bjerkeli, H. J. Fraser, S. B. Charnley

Abstract: During the journey from the cloud to the disc, the chemical composition of the protostellar envelope material can be either preserved or processed to varying degrees depending on the surrounding physical environment. This works aims to constrain the interplay of solid (ice) and gaseous methanol (CH$_3$OH) in the outer regions of protostellar envelopes located in the Coronet cluster in Corona Austr… ▽ More During the journey from the cloud to the disc, the chemical composition of the protostellar envelope material can be either preserved or processed to varying degrees depending on the surrounding physical environment. This works aims to constrain the interplay of solid (ice) and gaseous methanol (CH$_3$OH) in the outer regions of protostellar envelopes located in the Coronet cluster in Corona Australis (CrA), and assess the importance of irradiation by the Herbig Ae/Be star R CrA. CH$_3$OH is a prime test-case as it predominantly forms as a consequence of the solid-gas interplay (hydrogenation of condensed CO molecules onto the grain surfaces) and it plays an important role in future complex molecular processing. We present 1.3 mm Submillimeter Array (SMA) and Atacama Pathfinder Experiment (APEX) observations towards the envelopes of four low-mass protostars in the Coronet. Eighteen molecular transitions of seven species are identified. We calculate CH$_3$OH gas-to-ice ratios in this strongly irradiated cluster and compare them with ratios determined towards protostars located in less irradiated regions such as the Serpens SVS 4 cluster in Serpens Main and the Barnard 35A cloud in the $λ$ Orionis region. The CH$_3$OH gas-to-ice ratios in the Coronet vary by one order of magnitude (from 1.2$\times$10$^{-4}$ to 3.1$\times$10$^{-3}$) which is similar to less irradiated regions as found in previous studies. We find that the CH$_3$OH gas-to-ice ratios estimated in these three regions are remarkably similar despite the different UV radiation field intensities and formation histories. This result suggests that the overall CH$_3$OH chemistry in the outer regions of low-mass envelopes is relatively independent of variations in the physical conditions and hence that it is set during the prestellar stage. △ Less

Submitted 25 August, 2023; originally announced August 2023.

Comments: 20 pages, 17 figures, accepted for publication in Astronomy & Astrophysics

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

Protostellar Interferometric Line Survey of the Cygnus-X region (PILS-Cygnus) -- The role of the external environment in setting the chemistry of protostars

Authors: S. J. van der Walt, L. E. Kristensen, H. Calcutt, J. K. Jørgensen, R. T. Garrod

Abstract: (Abridged) Molecular lines are commonly detected towards protostellar sources. However, to get a better understanding of the chemistry of these sources we need unbiased molecular surveys over a wide frequency range for as many sources as possible to shed light on the origin of this chemistry, particularly any influence from the external environment. We present results from the PILS-Cygnus survey o… ▽ More (Abridged) Molecular lines are commonly detected towards protostellar sources. However, to get a better understanding of the chemistry of these sources we need unbiased molecular surveys over a wide frequency range for as many sources as possible to shed light on the origin of this chemistry, particularly any influence from the external environment. We present results from the PILS-Cygnus survey of ten intermediate- to high-mass protostellar sources in the nearby Cygnus-X complex, through high angular resolution interferometric observations over a wide frequency range. Using the Submillimeter Array (SMA), a spectral line survey of ten sources was performed in the frequency range 329-361 GHz, with an angular resolution of $\sim$1\farcs5, ($\sim$2000 AU, source distance of 1.3 kpc). Spectral modelling was performed to identify molecular emission and determine column densities and excitation temperatures for each source. We detect CH$_3$OH towards nine of the ten sources, CH$_3$OCH$_3$ and CH$_3$OCHO towards three sources, and CH$_3$CN towards four sources. Towards five sources the chemistry is spatially differentiated (different species peak at different positions and are offset from the peak continuum emission). The chemical properties of each source do not correlate with their position in the Cygnus-X complex, nor do the distance or direction to the nearest OB associations. However, the five sources located in the DR21 filament do appear to show less line emission compared to the five sources outside the filament. This work shows how important wide frequency coverage observations are combined with high angular resolution observations for studying the protostellar environment. Based on the ten sources observed here, the external environment appears to only play a minor role in setting the chemical environment on these small scales ($<$ 2000 AU). △ Less

Submitted 2 August, 2023; originally announced August 2023.

Comments: 19 pages, 77 pages appendices

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

Shock excitation of H$_2$ in the James Webb Space Telescope era

Authors: L. E. Kristensen, B. Godard, P. Guillard, A. Gusdorf, G. Pineau des Forets

Abstract: (Abridged) H2 is the most abundant molecule in the Universe. Thanks to its widely spaced energy levels, it predominantly lights up in warm gas, T > 100 K, such as shocked regions, and it is one of the key targets of JWST observations. These include shocks from protostellar outflows, all the way up to starburst galaxies and AGN. Shock models are able to simulate H2 emission. We aim to explore H2 ex… ▽ More (Abridged) H2 is the most abundant molecule in the Universe. Thanks to its widely spaced energy levels, it predominantly lights up in warm gas, T > 100 K, such as shocked regions, and it is one of the key targets of JWST observations. These include shocks from protostellar outflows, all the way up to starburst galaxies and AGN. Shock models are able to simulate H2 emission. We aim to explore H2 excitation using such models, and to test over which parameter space distinct signatures are produced in H2 emission. We present simulated H2 emission using the Paris-Durham shock code over an extensive grid of 14,000 plane-parallel stationary shock models, a large subset of which are exposed to an external UV radiation field. The grid samples 6 input parameters: preshock density, shock velocity, transverse magnetic field strength, UV radiation field strength, cosmic-ray-ionization rate, and PAH abundance. Physical quantities, such as temperature, density, and width, have been extracted along with H2 integrated line intensities. The strength of the transverse magnetic field, set by the scaling factor, b, plays a key role in the excitation of H2. At low values of b (<~ 0.3, J-type shocks), H2 excitation is dominated by vibrationally excited lines; at higher values (b >~ 1, C-type shocks), rotational lines dominate the spectrum for shocks with an external radiation field comparable to (or lower than) the solar neighborhood. Shocks with b >= 1 can be spatially resolved with JWST for nearby objects. When the input kinetic energy flux increases, the excitation and integrated intensity of H2 increases similarly. An external UV field mainly serves to increase the excitation, particularly for shocks where the input radiation energy is comparable to the input kinetic energy flux. These results provide an overview of the energetic reprocessing of input kinetic energy flux and the resulting H2 line emission. △ Less

Submitted 9 July, 2023; originally announced July 2023.

Comments: Published in A&A

Journal ref: A&A (2023), 675, A86

Episodic infall towards a compact disk in B335?

Authors: Per Bjerkeli, Jon P. Ramsey, Daniel Harsono, Adele Plunkett, Zhi-Yun Li, Matthijs H. D., van der Wiel, Hannah Calcutt, Jes K. Jørgensen, Lars E. Kristensen

Abstract: Previous observations of B335 have presented evidence of ongoing infall in various molecular lines, e.g., HCO$^+$, HCN, CO. There have been no confirmed observations of a rotationally supported disk on scales greater than ~12~au. The presence of an outflow in B335 suggests that also a disk should be present or in formation. To constrain the earliest stages of protostellar evolution and disk format… ▽ More Previous observations of B335 have presented evidence of ongoing infall in various molecular lines, e.g., HCO$^+$, HCN, CO. There have been no confirmed observations of a rotationally supported disk on scales greater than ~12~au. The presence of an outflow in B335 suggests that also a disk should be present or in formation. To constrain the earliest stages of protostellar evolution and disk formation, we aim to map the region where gas falls inwards and observationally constrain its kinematics. Furthermore, we aim to put strong limits on the size and orientation of any disk-like structure in B335. We use high angular resolution $^{13}$CO data from ALMA, and combine it with shorter-baseline archival data to produce a high-fidelity image of the infall in B335. We also revisit the imaging of high-angular resolution Band 6 continuum data to study the dust distribution in the immediate vicinity of B335. Continuum emission shows an elliptical structure (10 by 7 au) with a position angle 5 degrees east of north, consistent with the expectation for a forming disk in B335. A map of the infall velocity (as estimated from the $^{13}$CO emission), shows evidence of asymmetric infall, predominantly from the north and south. Close to the protostar, infall velocities appear to exceed free-fall velocities. 3D radiative transfer models, where the infall velocity is allowed to vary within the infall region, can explain the observed kinematics. The data suggests that a disk has started to form in B335 and that gas is falling towards that disk. However, kinematically-resolved line data towards the disk itself is needed to confirm the presence of a rotationally supported disk around this young protostar. The measured high infall velocities are not easily reconcilable with a magnetic braking scenario and suggest that there is a pressure gradient that allows the infall velocity to vary in the region. △ Less

Submitted 26 June, 2023; originally announced June 2023.

Comments: 14 pages, 11 figures

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

Investigating the Impact of Metallicity on Star Formation in the Outer Galaxy. I. VLT/KMOS Survey of Young Stellar Objects in Canis Major

Authors: Dominika Itrich, Agata Karska, Marta Sewiło, Lars E. Kristensen, Gregory J. Herczeg, Suzanne Ramsay, William J. Fischer, Benoît Tabone, Will R. M. Rocha, Maciej Koprowski, Ngân Lê, Beata Deka-Szymankiewicz

Abstract: The effects of metallicity on the evolution of protoplanetary disks may be studied in the outer Galaxy where the metallicity is lower than in the solar neighbourhood. We present the VLT/KMOS integral field spectroscopy in the near-infrared of $\sim$120 candidate young stellar objects (YSOs) in the CMa-$\ell$224 star-forming region located at a Galactocentric distance of 9.1 kpc. We characterise th… ▽ More The effects of metallicity on the evolution of protoplanetary disks may be studied in the outer Galaxy where the metallicity is lower than in the solar neighbourhood. We present the VLT/KMOS integral field spectroscopy in the near-infrared of $\sim$120 candidate young stellar objects (YSOs) in the CMa-$\ell$224 star-forming region located at a Galactocentric distance of 9.1 kpc. We characterise the YSO accretion luminosities and accretion rates using the hydrogen Br$γ$ emission and find the median accretion luminosity of $\log{(L_{\rm acc})} = -0.82^{+0.80}_{-0.82} L_\odot$. Based on the measured accretion luminosities, we investigate the hypothesis of star formation history in the CMa-$\ell$224. Their median values suggest that Cluster C, where most of YSO candidates have been identified, might be the most evolved part of the region. The accretion luminosities are similar to those observed toward low-mass YSOs in the Perseus and Orion molecular clouds, and do not reveal the impact of lower metallicity. Similar studies in other outer Galaxy clouds covering a wide range of metallicities are critical to gain a complete picture of star formation in the Galaxy. △ Less

Submitted 8 June, 2023; originally announced June 2023.

Comments: Accepted for publication in APJS, 51 pages, 37 figures, 6 tables

Far-infrared line emission from the outer Galaxy cluster Gy 3-7 with SOFIA/FIFI-LS: Physical conditions and UV fields

Authors: N. Le, A. Karska, M. Figueira, M. Sewiło, A. Mirocha, Ch. Fischer, M. Kaźmierczak-Barthel, R. Klein, M. Gawroński, M. Koprowski, K. Kowalczyk, W. J. Fischer, K. M. Menten, F. Wyrowski, C. König, L. E. Kristensen

Abstract: (abridged) Far-infrared (FIR) line emission provides key information about the gas cooling and heating due to shocks and UV radiation associated with the early stages of star formation. Gas cooling via FIR lines might, however, depend on metallicity. We aim to quantify the FIR line emission and determine the spatial distribution of the CO rotational temperature, ultraviolet (UV) radiation field, a… ▽ More (abridged) Far-infrared (FIR) line emission provides key information about the gas cooling and heating due to shocks and UV radiation associated with the early stages of star formation. Gas cooling via FIR lines might, however, depend on metallicity. We aim to quantify the FIR line emission and determine the spatial distribution of the CO rotational temperature, ultraviolet (UV) radiation field, and H2 number density toward the embedded cluster Gy 3-7 in the CMa-l224 star-forming region, whose metallicity is expected to be intermediate between that of the LMC and the Solar neighborhood. By comparing the total luminosities of CO and [O I] toward Gy 3-7 with values found for low- and high-mass protostars extending over a broad range of metallicities, we also aim to identify the possible effects of metallicity on the FIR line cooling within our Galaxy. We studied SOFIA/FIFI-LS spectra of Gy 3-7 covering several FIR lines. The spatial extent of CO high-J (J>14) emission resembles that of the elongated 160 um continuum emission detected with Herschel. The CO transitions from J=14-13 to J=16-15 are detected throughout the cluster and show a median rotational temperature of 170+/-30 K on Boltzmann diagrams. Comparisons to other protostars observed with Herschel show a good agreement with intermediate-mass sources in the inner Galaxy. Assuming an origin of the [O I] and high-J CO emission in UV-irradiated C-shocks, we obtained pre-shock H2 number densities of 10^4-5 cm-3 and UV radiation field strengths of 0.1-10 Habing fields. Far-IR line observations reveal ongoing star formation in Gy 3-7, dominated by intermediate-mass Class 0/I young stellar objects. The ratio of molecular-to-atomic far-IR line emission shows a decreasing trend with bolometric luminosities of the protostars. However, it does not indicate that the low-metallicity has an impact on the line cooling in Gy 3-7. △ Less

Submitted 17 April, 2023; originally announced April 2023.

Comments: 26 pages, 23 figures

Journal ref: A&A 674, A64 (2023)

Star-formation-rate estimates from water emission

Authors: K. M. Dutkowska, L. E. Kristensen

Abstract: (Abridged) The star-formation rate (SFR) quantitatively describes the star-formation process in galaxies. Current ways to calibrate this rate do not usually employ observational methods accounting for the low-mass end of stellar populations as their signatures are too weak. Accessing the bulk of protostellar activity within galactic star-forming regions can be achieved by tracing signposts of ongo… ▽ More (Abridged) The star-formation rate (SFR) quantitatively describes the star-formation process in galaxies. Current ways to calibrate this rate do not usually employ observational methods accounting for the low-mass end of stellar populations as their signatures are too weak. Accessing the bulk of protostellar activity within galactic star-forming regions can be achieved by tracing signposts of ongoing star formation. One such signpost is molecular outflows, which are bright in molecular emission. We propose to utilize the protostellar outflow emission as a tracer of the SFR. In this work, we introduce a novel version of the galaxy-in-a-box model, which can be used to relate molecular emission from star formation in galaxies with the SFR. We measured the predicted para-H2O emission at 988 GHz and corresponding SFRs for galaxies with LFIR = $10^8$ - $10^{11}$ L$_\odot$ in a distance-independent manner, and compared them with expectations from observations. We evaluated the derived results by varying the star formation efficiency, the free-fall time scaling factor, and the initial mass function. For the chosen H2O transition, relying on the current Galactic observations and star formation properties, we are underestimating the total galactic emission, while overestimating the SFRs, particularly for more starburst-like configurations. The current version of the galaxy-in-a-box model accounts for a limited number of processes and configurations, that is, it focuses on ongoing star formation in massive young clusters in a spiral galaxy. Therefore, the inferred results, which underestimate the emission and overestimate the SFR, are not surprising: known sources of emission are not included in the model. To improve the results, the next version of the model needs to include a more detailed treatment of the entire galactic ecosystem and other processes that would contribute to the emission. △ Less

Submitted 13 April, 2023; originally announced April 2023.

Comments: Accepted for publication in A&A. 11 pages, 6 figures

Journal ref: A&A 674, A95 (2023)

Simulation of CH$_3$OH ice UV photolysis under laboratory conditions

Authors: W. R. M. Rocha, P. Woitke, S. Pilling, W. -F. Thi, J. K. Jørgensen, L. E. Kristensen, G. Perotti, I. Kamp

Abstract: Methanol is the most complex molecule securely identified in interstellar ices and is a key chemical species for understanding chemical complexity in astrophysical environments. Important aspects of the methanol ice photochemistry are still unclear such as the branching ratios and photo-dissociation cross-sections at different temperatures and irradiation fluxes. This work aims at a quantitative a… ▽ More Methanol is the most complex molecule securely identified in interstellar ices and is a key chemical species for understanding chemical complexity in astrophysical environments. Important aspects of the methanol ice photochemistry are still unclear such as the branching ratios and photo-dissociation cross-sections at different temperatures and irradiation fluxes. This work aims at a quantitative agreement between laboratory experiments and astrochemical modelling of the CH3OH ice UV photolysis. This work also allows us to better understand which processes govern the methanol ice photochemistry present in laboratory experiments. We use ProDiMo to simulate the conditions of laboratory measurements. The simulations start with simple chemistry consisting only of methanol ice and helium to mimic the residual gas in the experimental chamber. A surface chemical network enlarged by photo-dissociation reactions is used to study the chemical reactions within the ice. Additionally, different surface chemistry parameters (surface competition, tunnelling, thermal diffusion and reactive desorption) are adopted to check those that reproduce the experimental results. The chemical models with ProDiMo can reproduce the methanol ice destruction via UV photodissociation at temperatures of 20, 30, 50 and 70 K as observed in the experiments. We note that the results are sensitive to different branching ratios after photolysis and to the mechanisms of reactive desorption. In the simulations of a molecular cloud at 20 K, we observed an increase in the methanol gas abundance of one order of magnitude, with a similar decrease in the solid-phase abundance. Comprehensive astrochemical models provide new insights into laboratory experiments as the quantitative understanding of the processes that govern the reactions within the ice. Ultimately, those insights can help to better interpret astronomical observations. △ Less

Submitted 8 February, 2023; originally announced February 2023.

Comments: Accepted for publication in A&A on 08-Feb-2023

Journal ref: A&A 673, A70 (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

Physical properties of accretion shocks toward the Class I protostellar system Oph-IRS 44

Authors: E. Artur de la Villarmois, V. V. Guzmán, J. K. Jørgensen, L. E. Kristensen, E. A. Bergin, D. Harsono, N. Sakai, E. F. van Dishoeck, S. Yamamoto

Abstract: (Abridged) Physical processes such as accretion shocks are thought to be common in the protostellar phase, where the envelope component is still present, and they can release molecules from the dust to the gas phase, altering the original chemical composition of the disk. Consequently, the study of accretion shocks is essential for a better understanding of the physical processes at disk scales an… ▽ More (Abridged) Physical processes such as accretion shocks are thought to be common in the protostellar phase, where the envelope component is still present, and they can release molecules from the dust to the gas phase, altering the original chemical composition of the disk. Consequently, the study of accretion shocks is essential for a better understanding of the physical processes at disk scales and their chemical output. The purpose of this work is to assess the characteristics of accretion shocks traced by sulfur-related species. We present ALMA high angular resolution observations (0.1") of the Class I protostar Oph-IRS 44. The continuum emission at 0.87 mm is observed, together with sulfur-related species such as SO, SO$_{2}$, and $^{34}$SO$_{2}$. Six lines of SO$_{2}$, two lines of $^{34}$SO$_{2}$, and one line of SO are detected toward IRS 44. The emission of all the detected lines peaks at ~0.1" (~14 au) from the continuum peak and we find infalling-rotating motions inside 30 au. However, only redshifted emission is seen between 50 and 30 au. Colder and more quiescent material is seen toward an offset region located at a distance of ~400 au from the protostar, and we do not find evidence of a Keplerian profile in these data. Accretion shocks are the most plausible explanation for the high temperatures, high densities, and velocities found for the SO$_{2}$ emission. When material enters the disk--envelope system, it generates accretion shocks that increase the dust temperature and desorb SO$_{2}$ molecules from dust grains. High-energy SO$_{2}$ transitions (~200 K) seem to be the best tracers of accretion shocks that can be followed up by future higher angular resolution ALMA observations and compared to other species to assess their importance in releasing molecules from the dust to the gas phase. △ Less

Submitted 6 September, 2022; originally announced September 2022.

Comments: 16 pages, 16 figures

Journal ref: A&A 667, A20 (2022)

Water emission tracing active star formation from the Milky Way to high-$z$ galaxies

Authors: K. M. Dutkowska, L. E. Kristensen

Abstract: (Abridged) The question of how most stars in the Universe form remains open. While star formation predominantly occurs in young massive clusters, the current framework focuses on isolated star formation. One way to access the bulk of protostellar activity within star-forming clusters is to trace signposts of active star formation with emission from molecular outflows. These outflows are bright in… ▽ More (Abridged) The question of how most stars in the Universe form remains open. While star formation predominantly occurs in young massive clusters, the current framework focuses on isolated star formation. One way to access the bulk of protostellar activity within star-forming clusters is to trace signposts of active star formation with emission from molecular outflows. These outflows are bright in water emission, providing a direct observational link between nearby and distant galaxies. We propose to utilize the knowledge of local star formation as seen with molecular tracers to explore the nature of star formation in the Universe. We present a large-scale statistical galactic model of emission from galactic active star-forming regions. Our model is built on observations of well-resolved nearby clusters. By simulating emission from molecular outflows, which is known to scale with mass, we create a proxy that can be used to predict the emission from clustered star formation at galactic scales. We evaluated the impact of the most important global-star formation parameters (i.e., initial stellar mass function (IMF), molecular cloud mass distribution, star formation efficiency (SFE), and free-fall time efficiency) on simulation results. We observe that for emission from the para-H2O 202 - 111 line, the IMF and molecular cloud mass distribution have a negligible impact on the emission, both locally and globally, whereas the opposite holds for the SFE and free-fall time efficiency. Moreover, this water transition proves to be a low-contrast tracer of star formation. The fine-tuning of the model and adaptation to morphologies of distant galaxies should result in realistic predictions of observed molecular emission and make the galaxy-in-a-box model a tool to analyze and better understand star formation throughout cosmological times. △ Less

Submitted 3 June, 2022; originally announced June 2022.

Comments: Accepted for publication in A&A. 16 pages, 13 figures

Journal ref: A&A 667, A135 (2022)

The Detection of Deuterated Water in the Large Magellanic Cloud with ALMA

Authors: Marta Sewiło, Agata Karska, Lars E. Kristensen, Steven B. Charnley, C. -H. Rosie Chen, Joana M. Oliveira, Martin Cordiner, Jennifer Wiseman, Álvaro Sánchez-Monge, Jacco Th. van Loon, Remy Indebetouw, Peter Schilke, Emmanuel Garcia-Berrios

Abstract: We report the first detection of deuterated water (HDO) toward an extragalactic hot core. The HDO 2$_{11}$-2$_{12}$ line has been detected toward hot cores N105-2A and 2B in the N105 star-forming region in the low-metallicity Large Magellanic Cloud (LMC) dwarf galaxy with the Atacama Large Millimeter/submillimeter Array (ALMA). We have compared the HDO line luminosity ($L_{\rm HDO}$) measured towa… ▽ More We report the first detection of deuterated water (HDO) toward an extragalactic hot core. The HDO 2$_{11}$-2$_{12}$ line has been detected toward hot cores N105-2A and 2B in the N105 star-forming region in the low-metallicity Large Magellanic Cloud (LMC) dwarf galaxy with the Atacama Large Millimeter/submillimeter Array (ALMA). We have compared the HDO line luminosity ($L_{\rm HDO}$) measured toward the LMC hot cores to those observed toward a sample of seventeen Galactic hot cores covering three orders of magnitude in $L_{\rm HDO}$, four orders of magnitude in bolometric luminosity ($L_{\rm bol}$), and a wide range of Galactocentric distances (thus metallicities). The observed values of $L_{\rm HDO}$ for the LMC hot cores fit very well into the $L_{\rm HDO}$ trends with $L_{\rm bol}$ and metallicity observed toward the Galactic hot cores. We have found that $L_{\rm HDO}$ seems to be largely dependent on the source luminosity, but metallicity also plays a role. We provide a rough estimate of the H$_2$O column density and abundance ranges toward the LMC hot cores by assuming that HDO/H$_2$O toward the LMC hot cores is the same as that observed in the Milky Way; the estimated ranges are systematically lower than Galactic values. The spatial distribution and velocity structure of the HDO emission in N105-2A is consistent with HDO being the product of the low-temperature dust grain chemistry. Our results are in agreement with the astrochemical model predictions that HDO is abundant regardless of the extragalactic environment and should be detectable with ALMA in external galaxies. △ Less

Submitted 9 May, 2022; originally announced May 2022.

Comments: 21 pages, 2 tables, 9 figures (including appendices); Accepted for publication in the Astrophysical Journal

A census of young stellar objects in two line-of-sight star forming regions toward IRAS 22147+5948 in the outer Galaxy

Authors: Agata Karska, Maciej Koprowski, Aleksandra Solarz, Ryszard Szczerba, Marta Sewiło, Natasza Siódmiak, Davide Elia, Marcin Gawroński, Konrad Grzesiak, Bosco H. K. Yung, William J. Fischer, Lars E. Kristensen

Abstract: (abridged) Star formation in the outer Galaxy, i.e., outside of the Solar circle, has been lightly studied in part due to low CO brightness of molecular clouds linked with the negative metallicity gradient. Recent infrared surveys provide an overview of dust emission in large sections of the Galaxy, but suffer from cloud confusion and poor spatial resolution at far-infrared wavelengths. We aim to… ▽ More (abridged) Star formation in the outer Galaxy, i.e., outside of the Solar circle, has been lightly studied in part due to low CO brightness of molecular clouds linked with the negative metallicity gradient. Recent infrared surveys provide an overview of dust emission in large sections of the Galaxy, but suffer from cloud confusion and poor spatial resolution at far-infrared wavelengths. We aim to develop a methodology to identify and classify young stellar objects (YSOs) in star-forming regions in the outer Galaxy, and use it to solve a long-standing confusion with the distance and evolutionary status of IRAS 22147+5948. We use Support Vector Machine learning algorithm to complement standard color-color and color-magnitude diagrams in search for YSOs in the IRAS 22147 region using publicly available data from the `Spitzer Mapping of the Outer Galaxy' survey. The agglomerative hierarchical clustering algorithm is used to identify clusters, along with the Robitaille et al. (2017) code to calculate physical properties of individual YSOs. We identify 13 Class I and 13 Class II YSO candidates using the color-color diagrams, and additional 2 and 21 sources, respectively, using the machine learning techniques. Spectral energy distributions of 23 sources are modelled with a star and a passive disk, corresponding to Class II objects. Models of 3 sources include envelopes typical for Class I objects. The objects are grouped in 2 clusters located at the distance of ~2.2 kpc, and 5 clusters at ~5.6 kpc. The spatial extent of CO, radio continuum, and dust emission confirms the origin of YSOs in two distinct star-forming regions along a similar line-of-sight. The outer Galaxy might serve as a unique laboratory of star formation across environments on condition that complementary methods and ancillary data are used to properly account for cloud confusion and distance uncertainties. △ Less

Submitted 11 May, 2022; v1 submitted 25 April, 2022; originally announced April 2022.

Comments: 29 pages, 17 figures, 9 tables

Journal ref: A&A 663, A133 (2022)

Connecting Galactic and extragalactic outflows: From the Cygnus-X cluster to active galaxies

Authors: I. M. Skretas, L. E. Kristensen

Abstract: Molecular outflows are often detected originating from both protostellar and extragalactic sources. Studies of low-mass, isolated high-mass, and extragalactic sources reveal scaling relations connecting the force carried by an outflow and the properties of the source that drives it. The aim of this work is to examine the effects of clustered star formation on the protostellar outflows and their sc… ▽ More Molecular outflows are often detected originating from both protostellar and extragalactic sources. Studies of low-mass, isolated high-mass, and extragalactic sources reveal scaling relations connecting the force carried by an outflow and the properties of the source that drives it. The aim of this work is to examine the effects of clustered star formation on the protostellar outflows and their scaling relations and to explore the possibility that outflows varying greatly scale and energetics are consistent with being launched by the same physical processes. High-angular resolution CO J = 3-2 observations were used of ten high-mass protostars in the Cygnus-X molecular cloud, obtained at the SMA as part of the PILS-Cygnus survey. From these data, the outflow force was measured. In addition, an extended sample of protostellar and extragalactic outflow-force measurements was assembled from existing literature to allow for a direct comparison of the scaling relations of the two types of outflows. Molecular outflows were detected originating from all ten sources of the PILS-Cygnus survey, and their outflow forces are found in close agreement with measurements from the literature. The comparison of the protostellar and extragalactic sources reveals, with 95\% confidence, that Class 0 protostars and extragalactic sources follow the same outflow force--bolometric luminosity correlation. The close agreement between the Cygnus-X sources and sources of similar envelope mass and bolometric luminosity suggests that clustered star formation has no significant effect on protostellar outflows. We find a strong indication that protostellar and extragalactic outflows are consistent with having a similar launch mechanism. The existence of such a mechanism would enable the development of a single universal outflow launch model, although more observations are required in order to verify this connection. △ Less

Submitted 9 February, 2022; originally announced February 2022.

Comments: 37 pages, 39 figures, 12 tables (including appendices); Accepted for publication in Astronomy and Astrophysics

Journal ref: A&A 660, A39 (2022)

ALMA Observations of Molecular Complexity in the Large Magellanic Cloud: The N105 Star-Forming Region

Authors: Marta Sewiło, Martin Cordiner, Steven B. Charnley, Joana M. Oliveira, Emmanuel Garcia Berrios, Peter Schilke, Jacob L. Ward, Jennifer Wiseman, Remy Indebetouw, Kazuki Tokuda, Jacco Th. van Loon, Álvaro Sánchez-Monge, Veronica Allen, C. -H. Rosie Chen, Roya Hamedani Golshan, Agata Karska, Lars E. Kristensen, Stan E. Kurtz, Toshikazu Onishi, Sarolta Zahorecz

Abstract: The Large Magellanic Cloud (LMC) is the nearest laboratory for detailed studies on the formation and survival of complex organic molecules (COMs), including biologically important ones, in low-metallicity environments--typical for earlier cosmological epochs. We report the results of 1.2 mm continuum and molecular line observations of three fields in the star-forming region N105 with the Atacama L… ▽ More The Large Magellanic Cloud (LMC) is the nearest laboratory for detailed studies on the formation and survival of complex organic molecules (COMs), including biologically important ones, in low-metallicity environments--typical for earlier cosmological epochs. We report the results of 1.2 mm continuum and molecular line observations of three fields in the star-forming region N105 with the Atacama Large Millimeter/submillimeter Array (ALMA). N105 lies at the western edge of the LMC bar with on-going star formation traced by H$_2$O, OH, and CH$_3$OH masers, ultracompact H II regions, and young stellar objects. Based on the spectral line modeling, we estimated rotational temperatures, column densities, and fractional molecular abundances for twelve 1.2 mm continuum sources. We identified sources with a range of chemical make-ups, including two bona fide hot cores and four hot core candidates. The CH$_3$OH emission is widespread and associated with all the continuum sources. COMs CH$_3$CN and CH$_3$OCH$_3$ are detected toward two hot cores in N105 together with smaller molecules typically found in Galactic hot cores (e.g., SO$_2$, SO, and HNCO) with the molecular abundances roughly scaling with metallicity. We report a tentative detection of the astrobiologically relevant formamide molecule (NH$_2$CHO) toward one of the hot cores; if confirmed, this would be the first detection of NH$_2$CHO in an extragalactic sub-solar metallicity environment. We suggest that metallicity inhomogeneities resulting from the tidal interactions between the LMC and the Small Magellanic Cloud (SMC) might have led to the observed large variations in COM abundances in LMC hot cores. △ Less

Submitted 24 January, 2022; originally announced January 2022.

Comments: 75 pages, 48 figures, 7 tables (including appendices); Accepted for publication in the Astrophysical Journal

Atomic Shocks in the Outflow of L1551 IRS 5 Identified with SOFIA-upGREAT Observations of [OI]

Authors: Yao-Lun Yang, Neal J. Evans II, Agata Karska, Lars E. Kristensen, Rebeca Aladro, Jon P. Ramsey, Joel D. Green, Jeong-Eun Lee

Abstract: We present velocity resolved SOFIA/upGREAT observations of [OI] and [CII] lines toward a Class I protostar, L1551 IRS 5, and its outflows. The SOFIA observations detect [OI] emission toward only the protostar and [CII] emission toward the protostar and the red-shifted outflow. The [OI] emission has a width of $\sim$100 km s$^{-1}$ only in the blue-shifted velocity, suggesting an origin in shocked… ▽ More We present velocity resolved SOFIA/upGREAT observations of [OI] and [CII] lines toward a Class I protostar, L1551 IRS 5, and its outflows. The SOFIA observations detect [OI] emission toward only the protostar and [CII] emission toward the protostar and the red-shifted outflow. The [OI] emission has a width of $\sim$100 km s$^{-1}$ only in the blue-shifted velocity, suggesting an origin in shocked gas. The [CII] lines are narrow, consistent with an origin in a photodissociation region. Differential dust extinction from the envelope due to the inclination of the outflows is the most likely cause of the missing red-shifted [OI] emission. Fitting the [OI] line profile with two Gaussian components, we find one component at the source velocity with a width of $\sim$20 km s$^{-1}$ and another extremely broad component at -30 km s$^{-1}$ with a width of 87.5 km s$^{-1}$, the latter of which has not been seen in L1551 IRS 5. The kinematics of these two components resemble cavity shocks in molecular outflows and spot shocks in jets. Radiative transfer calculations of the [OI], high-J CO, and H$_2$O lines in the cavity shocks indicate that [OI] dominates the oxygen budget, making up more than 70% of the total gaseous oxygen abundance and suggesting [O]/[H] of $\sim$1.5$\times$10$^{-4}$. Attributing the extremely broad [OI] component to atomic winds, we estimate the intrinsic mass loss rate of (1.3$\pm$0.8)$\times$ 10$^{-6}$ M$_{\odot}$ yr$^{-1}$. The intrinsic mass loss rates derived from low-J CO, [OI], and HI are similar, supporting the model of momentum-conserving outflows, where the atomic wind carries most momentum and drives the molecular outflows. △ Less

Submitted 22 November, 2021; originally announced November 2021.

Comments: 19 pages, 14 figures, accepted for publication in The Astrophysical Journal

Linking High- and Low-Mass Star Formation: Observation-Based Continuum Modelling and Physical Conditions

Authors: R. L. Pitts, L. E. Kristensen, J. K. Jørgensen, S. J. van der Walt

Abstract: Astronomers have yet to establish whether high-mass protostars form from high-mass prestellar cores, similar to their lower-mass counterparts, or from lower-mass fragments at the heart of a pre-protostellar cluster undergoing large-scale collapse. Part of the uncertainty is due to a shortage of envelope structure data on protostars of a few tens of solar masses, where we expect to see a transition… ▽ More Astronomers have yet to establish whether high-mass protostars form from high-mass prestellar cores, similar to their lower-mass counterparts, or from lower-mass fragments at the heart of a pre-protostellar cluster undergoing large-scale collapse. Part of the uncertainty is due to a shortage of envelope structure data on protostars of a few tens of solar masses, where we expect to see a transition from intermediate-mass star formation to the high-mass process. We sought to derive the masses, luminosities, and envelope density profiles for eight sources in Cygnus-X, whose mass estimates in the literature placed them in the sampling gap. Combining these sources with similarly evolved sources in the literature enabled us to perform a meta-analysis of protostellar envelope parameters over six decades in source luminosity. We performed spectral energy distribution (SED) fitting on archival broadband photometric continuum data from 1.2 to 850 $μ$m, to derive bolometric luminosities for our eight sources plus initial mass and radius estimates for modelling density and temperature profiles with the radiative transfer package Transphere. The envelope masses, densities at 1000 AU, outer envelope radii, and density power law indices as functions of bolometric luminosity all follow established trends in the literature spanning six decades in luminosity. Most of our sources occupy an intermediate to moderately high range of masses and luminosities, which helps to more firmly establish the continuity between low- and high-mass star formation mechanisms. Our density power law indices are consistent with observed values in literature, which show no discernible trends with luminosity. Finally, we show that the trends in all of the envelope parameters for high-mass protostars are statistically indistinguishable from trends in the same variables for low- and intermediate-mass protostars. △ Less

Submitted 21 October, 2021; originally announced October 2021.

Comments: 17 pages, 6 figures (+10 pages of Appendices); accepted to A&A

Journal ref: A&A 657, A70 (2022)

Signatures of UV radiation in low-mass protostars I. Origin of HCN and CN emission in the Serpens Main region

Authors: A. Mirocha, A. Karska, M. Gronowski, L. E. Kristensen, Ł. Tychoniec, D. Harsono, M. Figueira, M. Gładkowski, M. Żółtowski

Abstract: Context: Ultraviolet radiation (UV) influences the physics and chemistry of star-forming regions, but its properties and significance in the immediate surroundings of low-mass protostars are still poorly understood. Aims: We aim to extend the use of the CN/HCN ratio, already established for high-mass protostars, to the low-mass regime to trace and characterize the UV field around low-mass protosta… ▽ More Context: Ultraviolet radiation (UV) influences the physics and chemistry of star-forming regions, but its properties and significance in the immediate surroundings of low-mass protostars are still poorly understood. Aims: We aim to extend the use of the CN/HCN ratio, already established for high-mass protostars, to the low-mass regime to trace and characterize the UV field around low-mass protostars on $\sim 0.6\times0.6$ pc scales. Methods: We present $5'\times5'$ maps of the Serpens Main Cloud encompassing 10 protostars observed with the EMIR receiver at the IRAM 30 m telescope in CN 1-0, HCN 1-0, CS 3-2, and some of their isotopologues. The radiative-transfer code RADEX and the chemical model Nahoon are used to determine column densities of molecules, gas temperature and density, and the UV field strength, $G_\mathrm{0}$. Results: The spatial distribution of HCN and CS are well-correlated with CO 6-5 emission that traces outflows. The CN emission is extended from the central protostars to their immediate surroundings also tracing outflows, likely as a product of HCN photodissociation. The ratio of CN to HCN total column densities ranges from $\sim$1 to 12 corresponding to G$_0$ $\approx$ $10^{1}-10^{3}$ for gas densities and temperatures typical for outflows of low-mass protostars. Conclusions: UV radiation associated with protostars and their outflows is indirectly identified in a significant part of the Serpens Main low-mass star-forming region. Its strength is consistent with the values obtained from the OH and H$_2$O ratios observed with Herschel and compared with models of UV-illuminated shocks. From a chemical viewpoint, the CN to HCN ratio is an excellent tracer of UV fields around low- and intermediate-mass star-forming regions. △ Less

Submitted 14 September, 2021; originally announced September 2021.

Comments: 32 pages, 25 figures, accepted by A\&A

Journal ref: A&A 656, A146 (2021)

Protostellar Interferometric Line Survey of the Cygnus X region (PILS-Cygnus) -- First results: observations of CygX-N30

Authors: S. J. van der Walt, L. E. Kristensen, J. K. Jørgensen, H. Calcutt, S. Manigand, M. el Akel, R. T. Garrod, K. Qiu

Abstract: (Abridged) Complex organic molecules (COMs) are commonly detected in and near star-forming regions. However, the dominant process in the release of these COMs from the icy grains - where they predominately form - to the gas phase is still an open question. We investigate the origin of COM emission in a protostellar source, CygX-N30, through high-angular-resolution interferometric observations over… ▽ More (Abridged) Complex organic molecules (COMs) are commonly detected in and near star-forming regions. However, the dominant process in the release of these COMs from the icy grains - where they predominately form - to the gas phase is still an open question. We investigate the origin of COM emission in a protostellar source, CygX-N30, through high-angular-resolution interferometric observations over a continuous broad frequency range. We used 32 GHz Submillimeter Array observations with continuous frequency coverage from 329 to 361 GHz at an angular resolution of ~1" to do a line survey and obtain a chemical inventory of the source. The line emission was used to determine column densities and excitation temperatures for the COMs. We mapped out the intensity distribution of the different species and identified approximately 400 lines that can be attributed to 29 different molecular species and their isotopologues. We find that the molecular peak emission is along a linear gradient, coinciding with the axis of red- and blueshifted H2CO and CS emission. Chemical differentiation is detected along this gradient, with the O-bearing molecular species peaking towards one component of the system and the N- and S-bearing species peaking towards the other. The inferred column densities and excitation temperatures are compared to other sources where COMs are abundant. The origin of the observed COM emission is probably a combination of the young stellar sources along with accretion of infalling material onto a disc-like structure surrounding a young protostar. The low D/H ratio observed (<0.1%) likely reflects a pre-stellar phase where COMs formed on the ices at warm temperatures (~ 30 K), with inefficient deuterium fractionation. The observations and results presented here demonstrate the importance of good frequency coverage and high angular resolution when disentangling the origin of COM emission. △ Less

Submitted 8 September, 2021; originally announced September 2021.

Comments: 61 pages, 60 figures. To be published in Astronomy & Astrophysics

Journal ref: A&A 655, A86 (2021)

Fitting infrared ice spectra with genetic modelling algorithms. Presenting the ENIIGMA fitting tool

Authors: Will R. M. Rocha, Giulia Perotti, Lars E. Kristensen, Jes K. Jørgensen

Abstract: Context. A variety of laboratory ice spectra simulating different chemical environments, ice morphology as well as thermal and energetic processing are demanded to provide an accurate interpretation of the infrared spectra of protostars. To answer which combination of laboratory data best fit the observations, an automated statistically-based computational approach becomes necessary. Aims. To intr… ▽ More Context. A variety of laboratory ice spectra simulating different chemical environments, ice morphology as well as thermal and energetic processing are demanded to provide an accurate interpretation of the infrared spectra of protostars. To answer which combination of laboratory data best fit the observations, an automated statistically-based computational approach becomes necessary. Aims. To introduce a new approach, based on evolutionary algorithms, to search for molecules in ice mantles via spectral decomposition of infrared observational data with laboratory ice spectra. Methods. A publicly available and open-source fitting tool, called ENIIGMA (dEcompositioN of Infrared Ice features using Genetic Modelling Algorithms), is introduced. The tool has dedicated Python functions to carry out continuum determination of the protostellar spectra, silicate extraction, spectral decomposition and statistical analysis to calculate confidence intervals and quantify degeneracy. As an assessment of the code, several tests were conducted with known ice samples and constructed mixtures. A complete analysis of the Elias 29 spectrum was performed as well. Results. The ENIIGMA fitting tool can identify the correct ice samples and their fractions in all checks with known samples tested in this paper. Concerning the Elias 29 spectrum, the broad spectral range between 2.5-20 $μ$m was successfully decomposed after continuum determination and silicate extraction. This analysis allowed the identification of different molecules in the ice mantle, including a tentative detection of CH$_3$CH$_2$OH. Conclusions. The ENIIGMA is a toolbox for spectroscopy analysis of infrared spectra that is well-timed with the launch of the James Webb Space Telescope. Additionally, it allows for exploring the different chemical environments and irradiation fields in order to correctly interpret astronomical observations. △ Less

Submitted 18 July, 2021; originally announced July 2021.

Comments: 23 pages, 19 figures, 3 tables. Accepted for publication in A&A

Journal ref: A&A 654, A158 (2021)

Linking ice and gas in the Lambda Orionis Barnard 35A cloud

Authors: G. Perotti, J. K. Jørgensen, H. J. Fraser, A. N. Suutarinen, L. E. Kristensen, W. R. M. Rocha, P. Bjerkeli, K. M. Pontoppidan

Abstract: Dust grains play an important role in the synthesis of molecules in the interstellar medium, from the simplest species to complex organic molecules. How some of these solid-state molecules are converted into gas-phase species is still a matter of debate. Our aim is to directly compare ice and gas abundances of methanol (CH$_3$OH) and CO, and to investigate the relationship between ice and gas in l… ▽ More Dust grains play an important role in the synthesis of molecules in the interstellar medium, from the simplest species to complex organic molecules. How some of these solid-state molecules are converted into gas-phase species is still a matter of debate. Our aim is to directly compare ice and gas abundances of methanol (CH$_3$OH) and CO, and to investigate the relationship between ice and gas in low-mass protostellar envelopes. We present Submillimeter Array and Atacama Pathfinder EXperiment observations of gas-phase CH$_3$OH and CO towards the multiple protostellar system IRAS05417+0907 located in the B35A cloud. We use archival AKARI ice data toward the same target to calculate CH$_3$OH and CO gas-to-ice ratios. The CO isotopologues emissions are extended, whereas the CH$_3$OH emission is compact and traces the giant outflow emanating from IRAS05417+0907. A discrepancy between submillimeter dust emission and H$_2$O ice column density is found for B35A$-$4 and B35A$-$5, similar to what has previously been reported. B35A$-$2 and B35A$-$3 are located where the submillimeter dust emission peaks and show H$_2$O column densities lower than for B35A$-$4. The difference between the submillimeter continuum emission and the infrared H$_2$O ice observations suggests that the distributions of dust and H$_2$O ice differ around the young stellar objects in this dense cloud. The reason for this may be that the sources are located in different environments resolved by the interferometric observations: B35A$-$2, B35A$-$3 and in particular B35A$-$5 are situated in a shocked region plausibly affected by sputtering and heating impacting the submillimeter dust emission pattern, while B35A$-$4 is situated in a more quiescent part of the cloud. Gas and ice maps are essential to connect small-scale variations in the ice composition with large-scale astrophysical phenomena probed by gas observations. △ Less

Submitted 3 May, 2021; originally announced May 2021.

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

Journal ref: A&A 650, A168 (2021)

ALMA observations of doubly deuterated water: Inheritance of water from the prestellar environment

Authors: S. S. Jensen, J. K. Jørgensen, L. E. Kristensen, A. Coutens, E. F. van Dishoeck, K. Furuya, D. Harsono, M. V. Persson

Abstract: Establishing the origin of the water D/H ratio in the Solar System is central to our understanding of the chemical trail of water during the star and planet formation process. Recent modeling suggests that comparisons of the D$_2$O/HDO and HDO/H$_2$O ratios are a powerful way to trace the chemical evolution of water and, in particular, determine whether the D/H ratio is inherited from the molecula… ▽ More Establishing the origin of the water D/H ratio in the Solar System is central to our understanding of the chemical trail of water during the star and planet formation process. Recent modeling suggests that comparisons of the D$_2$O/HDO and HDO/H$_2$O ratios are a powerful way to trace the chemical evolution of water and, in particular, determine whether the D/H ratio is inherited from the molecular cloud or established locally. We seek to determine the D$_2$O column density and derive the D$_2$O/HDO ratios in the warm region toward the low-mass Class 0 sources B335 and L483. The results are compared with astrochemical models and previous observations to determine their implications for the chemical evolution of water. We present ALMA observations of the D$_2$O transition at 316.8 GHz toward B335 and L483 at $<$0.5" ($<$ 100 au) resolution, probing the inner warm envelope gas. The column densities of D$_2$O, HDO, and H$_{2}^{18}$O are determined by synthetic spectrum modeling and direct Gaussian fitting, under the assumption of a single excitation temperature and similar spatial extent for the three water isotopologs. D$_2$O is detected toward both sources in the inner warm envelope. The derived D$_2$O/HDO ratios is $(1.0\pm0.2)\times10^{-2}$ for L483 and $(1.4\pm0.1)\times10^{-2}$ for B335. The high D$_2$O/HDO ratios are a strong indication of chemical inheritance of water from the prestellar phase down to the inner warm envelope. This implies that the local cloud conditions in the prestellar phase, such as temperatures and timescales, determine the water chemistry at later stages and could provide a source of chemical differentiation in young systems. In addition, the observed D$_2$O/H$_2$O ratios support an observed dichotomy in the deuterium fractionation of water toward isolated and clustered protostars, namely, a higher D/H ratio toward isolated sources △ Less

Submitted 20 June, 2021; v1 submitted 27 April, 2021; originally announced April 2021.

Comments: Accepted for publication in A&A. Revision fixes typo in title

Journal ref: A&A 650, A172 (2021)

Water in star-forming regions (WISH): Physics and chemistry from clouds to disks as probed by Herschel spectroscopy

Authors: E. F. van Dishoeck, L. E. Kristensen, the WISH team

Abstract: (abridged) Data and results from the WISH key program are summarized, designed to provide a legacy data set to address its physics and chemistry. WISH targeted ~80 sources along the two axes of luminosity and evolutionary stage: from low- to high-mass protostars and from pre-stellar cores to protoplanetary disks. Lines of H2O, HDO, OH, CO and [O I] were observed with the HIFI and PACS instruments,… ▽ More (abridged) Data and results from the WISH key program are summarized, designed to provide a legacy data set to address its physics and chemistry. WISH targeted ~80 sources along the two axes of luminosity and evolutionary stage: from low- to high-mass protostars and from pre-stellar cores to protoplanetary disks. Lines of H2O, HDO, OH, CO and [O I] were observed with the HIFI and PACS instruments, complemented by molecules that probe UV, X-ray or grain chemistry. Most of the far-infrared water emission from protostars is found to be compact, originating from warm outflowing and shocked gas at high density and temperature in at least two physical components. This gas is not probed by low-J CO lines, only by J>14. Water is a significant, but not dominant, coolant. Its abundance is universally low, of order H2O/H2=2E-6, pointing to shock and outflow cavity models that include UV radiation at 100-1000 times the ISRF. In cold quiescent pre-stellar cores and envelopes, the water abundance structure is accurately probed through velocity-resolved line profiles, confirming basic chemistry networks. The gaseous HDO/H2O ratio of 0.025, much higher than that of bulk ice, is representative of the outer photodesorbed ice layers and cold chemistry. Water abundances in the inner hot cores are high, but with variations from 5E-6 to 2E-4. Combined analyses of water gas and ice show that up to 50% of the oxygen budget may be missing, with possible explanations discussed. Water vapor emission from disks is weak, indicating that water ice is locked up in larger pebbles early on and that these pebbles have settled and drifted inward by the Class II stage. Quantitatively, many oceans of water ice are available. Extragalactic low-J H2O emission is mostly compact and collisionally excited. Prospects for future mid- to far-infrared missions are given. △ Less

Submitted 3 February, 2021; originally announced February 2021.

Comments: 58 pages, 39 figures, accepted for publication in A&A

Journal ref: A&A 648, A24 (2021)

Time-Variable Radio Recombination Line Emission in W49A

Authors: C. G. De Pree, D. J. Wilner, L. E. Kristensen, R. Galván-Madrid, W. M. Goss, R. S. Klessen, M. -M. Mac Low, T. Peters, A. Robinson, S. Sloman, M. Rao

Abstract: We present new Jansky Very Large Array (VLA) images of the central region of the W49A star-forming region at 3.6~cm and at 7~mm at resolutions of 0\farcs15 (1650 au) and 0\farcs04 (440 au), respectively. The 3.6~cm data reveal new morphological detail in the ultracompact \ion{H}{2} region population, as well as several previously unknown and unresolved sources. In particular, source A shows elonga… ▽ More We present new Jansky Very Large Array (VLA) images of the central region of the W49A star-forming region at 3.6~cm and at 7~mm at resolutions of 0\farcs15 (1650 au) and 0\farcs04 (440 au), respectively. The 3.6~cm data reveal new morphological detail in the ultracompact \ion{H}{2} region population, as well as several previously unknown and unresolved sources. In particular, source A shows elongated, edge-brightened, bipolar lobes, indicative of a collimated outflow, and source E is resolved into three spherical components. We also present VLA observations of radio recombination lines at 3.6~cm and 7~mm, and IRAM Northern Extended Millimeter Array (NOEMA) observations at 1.2~mm. Three of the smallest ultracompact \ion{H}{2} regions (sources A, B2 and G2) all show broad kinematic linewidths, with $Δ$V$_{FWHM}\gtrsim$40~km~s$^{-1}$. A multi-line analysis indicates that broad linewidths remain after correcting for pressure broadening effects, suggesting the presence of supersonic flows. Substantial changes in linewidth over the 21 year time baseline at both 3.6 cm and 7 mm are found for source G2. At 3.6 cm, the linewidth of G2 changed from 31.7$\pm$1.8 km s$^{-1}$ to 55.6$\pm$2.7 km s$^{-1}$, an increase of $+$23.9$\pm$3.4 km s$^{-1}$. The G2 source was previously reported to have shown a 3.6~cm continuum flux density decrease of 40\% between 1994 and 2015. This source sits near the center of a very young bipolar outflow whose variability may have produced these changes. △ Less

Submitted 9 October, 2020; originally announced October 2020.

Comments: 18 pages, 7 figures, Accepted for publication in The Astronomical Journal

Linking ice and gas in the Serpens low-mass star-forming region

Authors: G. Perotti, W. R. M. Rocha, J. K. Jørgensen, L. E. Kristensen, H. J. Fraser, K. M. Pontoppidan

Abstract: The interaction between dust, ice, and gas during the formation of stars produces complex organic molecules. While observations indicate that several species are formed on ice-covered dust grains and are released into the gas phase, the exact chemical interplay between solid and gas phases and their relative importance remain unclear. Our goal is to study the interplay in regions of low-mass star… ▽ More The interaction between dust, ice, and gas during the formation of stars produces complex organic molecules. While observations indicate that several species are formed on ice-covered dust grains and are released into the gas phase, the exact chemical interplay between solid and gas phases and their relative importance remain unclear. Our goal is to study the interplay in regions of low-mass star formation through ice- and gas-mapping and by directly measuring gas-to-ice ratios. This provides constraints on the routes that lead to the chemical complexity that is observed in both phases. We present observations of gas-phase methanol (CH$_3$OH) and carbon monoxide at 1.3 mm towards ten low-mass young protostars in the Serpens SVS4 cluster from the SubMillimeter Array and the Atacama Pathfinder EXperiment telescope. We used archival data from the Very Large Telescope to derive abundances of ice H$_2$O, CO, and CH$_3$OH towards the same region. Finally, we constructed gas-ice maps of SVS4 and directly measured CO and CH$_3$OH gas-to-ice ratios. The CH$_3$OH gas-to-ice ratio agrees with values that were previously reported for embedded Class 0/I low-mass protostars. The CO gas-maps trace an extended gaseous component that is not sensitive to the effect of freeze-out. We find that there is no straightforward correlation between CO and CH$_3$OH gas with their ice counterparts in the cluster. This is likely related to the complex morphology of SVS4: the Class 0 protostar SMM4 and its envelope lie in the vicinity, and the outflow associated with SMM4 intersects the cluster. This study serves as a pathfinder for future observations with ALMA and the James Webb Space Telescope that will provide high-sensitivity gas-ice maps of molecules more complex than methanol. Such comparative maps will be essential to constrain the chemical routes that regulate the chemical complexity in star-forming regions. △ Less

Submitted 6 August, 2020; originally announced August 2020.

Comments: 25 pages, 20 figures, accepted for publication by A&A

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

Steady Wind-Blown Cavities within Infalling Rotating Envelopes: Application to the Broad Velocity Component in Young Protostars

Authors: Lichen Liang, Doug Johnstone, Sylvie Cabrit, Lars E. Kristensen

Abstract: Wind-driven outflows are observed around a broad range of accreting objects throughout the Universe, ranging from forming low-mass stars to super-massive black holes. We study the interaction between a central isotropic wind and an infalling, rotating, envelope, determining the steady-state cavity shape formed at their interface under the assumption of weak mixing. The shape of the resulting wind-… ▽ More Wind-driven outflows are observed around a broad range of accreting objects throughout the Universe, ranging from forming low-mass stars to super-massive black holes. We study the interaction between a central isotropic wind and an infalling, rotating, envelope, determining the steady-state cavity shape formed at their interface under the assumption of weak mixing. The shape of the resulting wind-blown cavity is elongated and self-similar, with a physical size determined by the ratio between wind ram pressure and envelope thermal pressure. We compute the growth of a warm turbulent mixing-layer between the shocked wind and the deflected envelope, and calculate the resultant broad line profile, under the assumption of a linear (Couette-type) velocity profile across the layer. We then test our model against the warm broad velocity component observed in CO $J$=16--15 by Herschel/HIFI in the protostar Serpens-Main SMM1. Given independent observational constraints on the temperature and density of the dust envelope around SMM1, we find an excellent match to all its observed properties (line profile, momentum, temperature) and to the SMM1 outflow cavity width for a physically reasonable set of parameters: a ratio of wind to infall mass-flux $\simeq 4\%$, a wind speed $v_{\rm w} \simeq 30$ km/s, an interstellar abundance of CO and H$_2$, and a turbulent entrainment efficiency consistent with laboratory experiments. The inferred ratio of ejection to disk accretion rate, $\simeq 6-20\%$, is in agreement with current disk wind theories. Thus, the model provides a new framework to reconcile the modest outflow cavity widths in protostars with the large observed flow velocities. Being self-similar, it is applicable over a broader range of astrophysical contexts as well. △ Less

Submitted 27 July, 2020; originally announced July 2020.

Comments: 31 pages, 21 figures, accepted to ApJ for publication (comments are welcome)

Feedback of molecular outflows from protostars in NGC 1333, revealed by Herschel and Spitzer spectro-imaging observations

Authors: Odysseas Dionatos, Lars. E. Kristensen, Mario Tafalla, Manuel Güdel, Magnus Persson

Abstract: Large scale spectral maps of star forming regions enable the comparative study of the gas excitation around an ensemble of sources at a common frame of reference, providing direct insights in the multitude of processes involved. In this paper we employ spectral-line maps to decipher the excitation, the kinematical and dynamical processes in NGC 1333 as revealed by a number of different emission li… ▽ More Large scale spectral maps of star forming regions enable the comparative study of the gas excitation around an ensemble of sources at a common frame of reference, providing direct insights in the multitude of processes involved. In this paper we employ spectral-line maps to decipher the excitation, the kinematical and dynamical processes in NGC 1333 as revealed by a number of different emission lines, aiming to set a reference for the applicability of tracers in constraining diverse physical processes. We reconstruct line maps for H$_2$ , CO, H$_2$O and C$^+$ using data obtained with the Spitzer-IRS and Herschel HIFI-SPIRE. We compare the morphological features of the maps and derive the gas excitation for regions of interest employing LTE and non-LTE methods. We also calculate the kinematical and dynamical properties for each outflow tracer consistently for all outflows in NGC 1333. We finally measure the water abundance in outflows with respect to carbon monoxide and molecular hydrogen. CO and H$_2$ are highly excited around B-stars and at lower levels trace protostellar outflows. H$_2$O emission is dominated by a moderately fast component associated with outflows. Intermediate J CO lines appear brightest at the locations traced by a narrow H$_2$O component, indicating that beyond the dominating collisional processes, a secondary, radiative excitation component can also be active. The morphology, kinematics, excitation and abundance variations of water are consistent with its excitation and partial dissociation in shocks. Water abundance ranges between 5 x 10$^{-7}$ and 10$^{-5}$, with the lower values being more representative. Water is brightest and most abundant around IRAS 4A which is consistent with the latter hosting a hot corino source. Finally, the outflow mass flux is found highest for CO and decreases by one and two orders of magnitude for H$_2$ and H$_2$O, respectively. △ Less

Submitted 1 June, 2020; originally announced June 2020.

Comments: 19 pages, 14 figures, 5 tables, A&A accepted

Journal ref: A&A 641, A36 (2020)

Mass Assembly of Stellar Systems and their Evolution with the SMA (MASSES) -- Full Data Release

Authors: Ian W. Stephens, Tyler L. Bourke, Michael M. Dunham, Philip C. Myers, Riwaj Pokhrel, John J. Tobin, Héctor G. Arce, Sarah I. Sadavoy, Eduard I. Vorobyov, Jaime E. Pineda, Stella S. R. Offner, Katherine I. Lee, Lars E. Kristensen, Jes K. Jørgensen, Mark A. Gurwell, Alyssa A. Goodman

Abstract: We present and release the full dataset for the Mass Assembly of Stellar Systems and their Evolution with the SMA (MASSES) survey. This survey used the Submillimeter Array (SMA) to image the 74 known protostars within the Perseus molecular cloud. The SMA was used in two array configurations to capture outflows for scales $>$30$^{\prime\prime}$ ($>$9000 au) and to probe scales down to $\sim$1… ▽ More We present and release the full dataset for the Mass Assembly of Stellar Systems and their Evolution with the SMA (MASSES) survey. This survey used the Submillimeter Array (SMA) to image the 74 known protostars within the Perseus molecular cloud. The SMA was used in two array configurations to capture outflows for scales $>$30$^{\prime\prime}$ ($>$9000 au) and to probe scales down to $\sim$1$^{\prime\prime}$ ($\sim$300 au). The protostars were observed with the 1.3 mm and 850 $μ$m receivers simultaneously to detect continuum at both wavelengths and molecular line emission from CO(2-1), $^{13}$CO(2-1), C$^{18}$O(2-1), N$_2$D$^+$(3-2), CO(3-2), HCO$^+$(4-3), and H$^{13}$CO$^+$(4-3). Some of the observations also used the SMA's recently upgraded correlator, SWARM, whose broader bandwidth allowed for several more spectral lines to be observed (e.g., SO, H$_2$CO, DCO$^+$, DCN, CS, CN). Of the main continuum and spectral tracers observed, 84% of the images and cubes had emission detected. The median C$^{18}$O(2-1) linewidth is $\sim$1.0 km s$^{-1}$, which is slightly higher than those measured with single-dish telescopes at scales of 3000-20000 au. Of the 74 targets, six are suggested to be first hydrostatic core candidates, and we suggest that L1451-mm is the best candidate. We question a previous continuum detection toward L1448 IRS2E. In the SVS13 system, SVS13A certainly appears to be the most evolved source, while SVS13C appears to be hotter and more evolved than SVS13B. The MASSES survey is the largest publicly available interferometric continuum and spectral line protostellar survey to date, and is largely unbiased as it only targets protostars in Perseus. All visibility ($uv$) data and imaged data are publicly available at https://dataverse.harvard.edu/dataverse/full_MASSES/. △ Less

Submitted 19 November, 2019; originally announced November 2019.

Comments: Accepted to ApJS

Chemical and kinematic structure of extremely high-velocity molecular jets in the Serpens Main star-forming region

Authors: Łukasz Tychoniec, Charles L. H. Hull, Lars E. Kristensen, John J. Tobin, Valentin J. M. Le Gouellec, Ewine F. van Dishoeck

Abstract: The fastest molecular component to the protostellar outflows -- extremely high-velocity (EHV) molecular jets -- are still puzzling since they are seen only rarely. The first aim is to analyze the interaction between the EHV jet and the slow outflow by comparing their outflow force content. The second aim is to analyze the chemical composition of the different outflow velocity components and to rev… ▽ More The fastest molecular component to the protostellar outflows -- extremely high-velocity (EHV) molecular jets -- are still puzzling since they are seen only rarely. The first aim is to analyze the interaction between the EHV jet and the slow outflow by comparing their outflow force content. The second aim is to analyze the chemical composition of the different outflow velocity components and to reveal the spatial location of molecules. ALMA 3 mm and 1.3 mm observations of five outflow sources at 130 -- 260 au resolution in the Serpens Main cloud are presented. Observations of CO, SiO, H$_2$CO and HCN reveal the kinematic and chemical structure of those flows. Three velocity components are distinguished: the slow and the fast wing, and the EHV jet. Out of five sources, three have the EHV component. Comparison of outflow forces reveals that only the EHV jet in the youngest source Ser-emb 8 (N) has enough momentum to power the slow outflow. The SiO abundance is generally enhanced with velocity, while HCN is present in the slow and the fast wing, but disappears in the EHV jet. For Ser-emb 8 (N), HCN and SiO show a bow-shock shaped structure surrounding one of the EHV peaks suggesting sideways ejection creating secondary shocks upon interaction with the surroundings. Also, the SiO abundance in the EHV gas decreases with distance from this protostar, whereas that in the fast wing increases. H$_2$CO is mostly associated with low-velocity gas but also appears surprisingly in one of the bullets in the Ser-emb~8~(N) EHV jet. The high detection rate suggests that the presence of the EHV jet may be more common than previously expected. The origin and temporal evolution of the abundances of SiO, HCN and H$_2$CO through high-temperature chemistry are discussed. The data are consistent with a low C/O ratio in the EHV gas versus high C/O ratio in the fast and slow wings. △ Less

Submitted 17 October, 2019; originally announced October 2019.

Comments: 25 pages, 20 figures, Accepted for publication in the Astronomy & Astrophysics

Journal ref: A&A 632, A101 (2019)

ALMA observations of water deuteration: A physical diagnostic of the formation of protostars

Authors: S. S. Jensen, J. K. Jørgensen, L. E. Kristensen, K. Furuya, A. Coutens, E. F. van Dishoeck, D. Harsono, M. V. Persson

Abstract: The deuterium fractionation of water can serve as a tracer for the chemical and physical evolution of water during star formation and can constrain the origin of water in Solar System bodies. We determine the HDO/H$_2$O ratio in the inner warm gas toward three low-mass Class 0 protostars selected to be in isolated cores, i.e., not associated with any cloud complexes. Previous sources for which the… ▽ More The deuterium fractionation of water can serve as a tracer for the chemical and physical evolution of water during star formation and can constrain the origin of water in Solar System bodies. We determine the HDO/H$_2$O ratio in the inner warm gas toward three low-mass Class 0 protostars selected to be in isolated cores, i.e., not associated with any cloud complexes. Previous sources for which the HDO/H$_2$O ratio have been established were all part of larger star-forming complexes. Targeting these isolated protostars allows comparison of the water chemistry in isolated and clustered regions to determine the influence of local cloud environment. We present ALMA observations of the HDO $3_{1,2}$-$2_{2,1}$ and $2_{1,1}$-$2_{1,2}$ transitions at 225.897 GHz and 241.562 GHz along with the H$_2^{18}$O $3_{1,3}$-$2_{2,0}$ transition at 203.407 GHz. The high angular resolution (0\farcs3-1\farcs3) allow the study of the inner warm envelope gas. Model-independent estimates for the HDO/H$_2$O ratios are obtained and compared with previous determinations in the warm gas toward low-mass protostars. We detect the targeted water transitions toward the three sources with S/N > 5. We determine the HDO/H$_2$O ratio toward L483, B335 and BHR71-IRS1 to be ($2.2\pm0.4$)$\times 10^{-3}$, ($1.7\pm0.3$)$\times 10^{-3}$, and ($1.8\pm0.4$)$\times 10^{-3}$, respectively, assuming $T_\mathrm{ex} = 124$ K. The degree of water deuteration of these isolated protostars are a factor of 2-4 higher relative to Class 0 protostars that are members of known nearby clustered star-forming regions. The results indicate that the water deuterium fractionation is influenced by the local cloud environment. This effect can be explained by variations in either collapse timescales or temperatures, which depends on local cloud dynamics and could provide a new method to decipher the history of young stars. △ Less

Submitted 23 September, 2019; originally announced September 2019.

Comments: 13 pages, 6 figures. Accepted for publication in A&A (18/09/2019)

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

Kinematics around the B335 protostar down to au scales

Authors: Per Bjerkeli, Jon P. Ramsey, Daniel Harsono, Hannah Calcutt, Lars E. Kristensen, Matthijs H. D. van der Wiel, Jes K. Jørgensen, Sébastien Muller, Magnus V. Persson

Abstract: Context. The relationship between outflow launching and formation of accretion disks around young stellar objects is still not entirely understood, which is why spectrally and spatially resolved observations are needed. Recently, the Atacama Large Millimetre/sub-millimetre Array (ALMA) has carried out long-baseline observations towards a handful of sources, revealing connections between outflows a… ▽ More Context. The relationship between outflow launching and formation of accretion disks around young stellar objects is still not entirely understood, which is why spectrally and spatially resolved observations are needed. Recently, the Atacama Large Millimetre/sub-millimetre Array (ALMA) has carried out long-baseline observations towards a handful of sources, revealing connections between outflows and the inner regions of disks. Aims. Here we aim to determine the small-scale kinematic and morphological properties of the outflow from the isolated protostar B335 for which no Keplerian disk has, so far, been observed on scales down to 10 au. Methods. We use ALMA in its longest-baseline configuration to observe emission from CO isotopologs, SiO, SO$_2$ and CH$_3$OH. The proximity of B335 provides a resolution of ~3 au (0.03''). We also combine our long-baseline data with archival data to produce a high-fidelity image covering scales up to 700 au (7''). Results. $^{12}$CO has a X-shaped morphology with arms ~50 au in width that we associate with the walls of an outflow cavity, similar to what is observed on larger scales. Long-baseline continuum emission is confined to <7 au of the protostar, while short-baseline continuum emission follows the $^{12}$CO outflow and cavity walls. Methanol is detected within ~30 au of the protostar. SiO is also detected in the vicinity of the protostar, but extended along the outflow. Conclusions. The $^{12}$CO outflow shows no clear signs of rotation at distances $\gtrsim$30 au from the protostar. SiO traces the protostellar jet on small scales, but without obvious rotation. CH$_3$OH and SO$_2$ trace a region <16 au in diameter, centred on the continuum peak, which is clearly rotating. Using episodic, high-velocity, $^{12}$CO features, we estimate the launching radius of the outflow to be <0.1 au and dynamical timescales on the order of a few years. △ Less

Submitted 11 September, 2019; originally announced September 2019.

Comments: Abstract very slightly abridged. 15 pages, 13 figures, accepted for publication in A&A

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

Characterizing magnetic field morphologies in three Serpens protostellar cores with ALMA

Authors: Valentin J. M. Le Gouellec, Charles L. H. Hull, Anaëlle J. Maury, Josep M. Girart, Łukasz Tychoniec, Lars E. Kristensen, Zhi-Yun Li, Fabien Louvet, Paulo C. Cortes, Ramprasad Rao

Abstract: With the aim of characterizing the dynamical processes involved in the formation of young protostars, we present high angular resolution ALMA dust polarization observations of the Class 0 protostellar cores Serpens SMM1, Emb 8(N), and Emb 8. With spatial resolutions ranging from 150 to 40 au at 870 μm, we find unexpectedly high values of the polarization fraction along the outflow cavity walls in… ▽ More With the aim of characterizing the dynamical processes involved in the formation of young protostars, we present high angular resolution ALMA dust polarization observations of the Class 0 protostellar cores Serpens SMM1, Emb 8(N), and Emb 8. With spatial resolutions ranging from 150 to 40 au at 870 μm, we find unexpectedly high values of the polarization fraction along the outflow cavity walls in Serpens Emb8(N). We use 3 mm and 1 mm molecular tracers to investigate outflow and dense gas properties and their correlation with the polarization. These observations allow us to investigate the physical processes involved in the Radiative Alignment Torques (RATs) acting on dust grains along the outflow cavity walls, which experience irradiation from accretion processes and outflow shocks. The inner core of SMM1-a presents a polarization pattern with a poloidal magnetic field at the bases of the two lobes of the bipolar outflow. To the south of SMM1-a we see two polarized filaments, one of which seems to trace the redshifted outflow cavity wall. The other may be an accretion streamer of material infalling onto the central protostar. We propose that the polarized emission we see at millimeter wavelengths along the irradiated cavity walls can be reconciled with the expectations of RAT theory if the aligned grains present at '<' 500 au scales in Class 0 envelopes have grown larger than the 0.1 μm size of ISM dust grains. Our observations allow us to constrain the star-forming sources magnetic field morphologies within the central cores, along the outflow cavity walls, and in possible accretion streamers. △ Less

Submitted 30 August, 2019; originally announced September 2019.

The Bridge: a transient phenomenon of forming stellar multiples

Authors: Michael Kuffmeier, Hannah Calcutt, Lars E. Kristensen

Abstract: Observations with modern instruments such as Herschel reveal that stars form clustered inside filamentary arms of ~1 pc length embedded in Giant Molecular Clouds. On smaller scales (~1000 au), observations of, e.g., IRAS 16293--2422 show signs of filamentary `bridge' structures connecting young protostars in their birth environment. We investigate the formation of a deeply embedded triple protoste… ▽ More Observations with modern instruments such as Herschel reveal that stars form clustered inside filamentary arms of ~1 pc length embedded in Giant Molecular Clouds. On smaller scales (~1000 au), observations of, e.g., IRAS 16293--2422 show signs of filamentary `bridge' structures connecting young protostars in their birth environment. We investigate the formation of a deeply embedded triple protostellar system, and characterize its connection to the filamentary structure present on GMC scales. Using the magnetohydrodynamical code RAMSES, we carry out zoom-in simulations of low-mass star formation starting from Giant-Molecular-Cloud-scales. Colliding flows of gas in the filamentary arms induce the formation of two protostellar companions at distances of ~1000 au from the primary. After their birth, the stellar companions quickly (~10 kyr) approach and orbit the primary on eccentric orbits with separations of ~100 au. The colliding flows induce transient structures lasting for up to a few 10 kyr connecting two forming protostellar objects that are kinematically quiescent along the line-of-sight. Colliding flows compress the gas and trigger the formation of stellar companions via turbulent fragmentation. Our results suggest that protostellar companions form with a wide separation of ~1000 au. Smaller separations (~100 au or less) are a consequence of subsequent migration and capturing. Associated with the formation phase of the companion, the turbulent environment induces arc- and bridge-like structures. These bridges can become kinematically quiescent, when the velocity components of the colliding flows eliminate each other. However, the gas in bridges still contributes to stellar accretion later. Our results demonstrate: bridge-like structures are a transient phenomenon of stellar multiple formation. △ Less

Submitted 3 July, 2019; originally announced July 2019.

Comments: 12 pages, 9 figures, accepted for publication in A&A

Journal ref: A&A 628, A112 (2019)

Revealing the chemical structure of the Class I disc Oph-IRS 67

Authors: E. Artur de la Villarmois, L. E. Kristensen, J. K. Jørgensen

Abstract: (Abridged) The purpose of this paper is to explore the structure of a line-rich Class I protobinary source, Oph-IRS 67, and analyse the differences and similarities with Class 0 and Class II sources. We present a systematic molecular line study of IRS 67 with the Submillimeter Array (SMA) on 1 - 2" (150 - 300 AU) scales. The wide instantaneous band-width of the SMA observations (~30 GHz) provide d… ▽ More (Abridged) The purpose of this paper is to explore the structure of a line-rich Class I protobinary source, Oph-IRS 67, and analyse the differences and similarities with Class 0 and Class II sources. We present a systematic molecular line study of IRS 67 with the Submillimeter Array (SMA) on 1 - 2" (150 - 300 AU) scales. The wide instantaneous band-width of the SMA observations (~30 GHz) provide detections of a range of molecular transitions that trace different physics, such as CO isotopologues, sulphur-bearing species, deuterated species, and carbon-chain molecules. We see significant differences between different groups of species. For example, the CO isotopologues and sulphur-bearing species show a rotational profile and are tracing the larger-scale circumbinary disc structure, while CN, DCN, and carbon-chain molecules peak at the southern edge of the disc at blue-shifted velocities. In addition, the cold gas tracer DCO+ is seen beyond the extent of the circumbinary disc. The detected molecular transitions can be grouped into three main components: cold regions far from the system, the circumbinary disc, and a UV-irradiated region likely associated with the surface layers of the disc that are reached by the UV radiation from the sources. The different components are consistent with the temperature structure derived from the ratio of two H2CO transitions, that is, warm temperatures are seen towards the outflow direction, lukewarm temperatures are associated with the UV-radiated region, and cold temperatures are related with the circumbinary disc structure. The chemistry towards IRS 67 shares similarities with both Class 0 and Class II sources, possibly due to the high gas column density and the strong UV radiation arising from the binary system. IRS 67 is, therefore, highlighting the intermediate chemistry between deeply embedded sources and T-Tauri discs. △ Less

Submitted 3 June, 2019; originally announced June 2019.

Comments: 15 pages, 12 Figures, 5 Tables

Journal ref: A&A 627, A37 (2019)

Physical and chemical fingerprint of protostellar disc formation

Authors: E. Artur de la Villarmois, J. K. Jørgensen, L. E. Kristensen, E. A. Bergin, D. Harsono, N. Sakai, E. F. van Dishoeck, S. Yamamoto

Abstract: (Abridged) The purpose of this paper is to explore and compare the physical and chemical structure of Class I low-mass protostellar sources on protoplanetary disc scales. We present a study of the dust and gas emission towards a representative sample of 12 Class I protostars from the Ophiuchus molecular cloud with the Atacama Large Millimeter/submillimeter Array (ALMA). The continuum at 0.87 mm an… ▽ More (Abridged) The purpose of this paper is to explore and compare the physical and chemical structure of Class I low-mass protostellar sources on protoplanetary disc scales. We present a study of the dust and gas emission towards a representative sample of 12 Class I protostars from the Ophiuchus molecular cloud with the Atacama Large Millimeter/submillimeter Array (ALMA). The continuum at 0.87 mm and molecular transitions from C17O, C34S, H13CO+, CH3OH, SO2 , and C2H were observed at high angular resolution (0.4", ~60 au diameter) towards each source. Disc and stellar masses are estimated from the continuum flux and position-velocity diagrams, and six of the sources show disc-like structures. Towards the more luminous sources, compact emission and large line widths are seen for transitions of SO2 that probe warm gas (Eu ~200 K). In contrast, C17O emission is detected towards the least evolved and less luminous systems. No emission of CH3OH is detected towards any of the continuum peaks, indicating an absence of warm CH3OH gas towards these sources. A power-law relation is seen between the stellar mass and the bolometric luminosity, corresponding to a mass accretion rate of (2.4 +/- 0.6) x 10^-7 Msun/year for the Class I sources. This mass accretion rate is lower than the expected value if the accretion is constant in time and rather points to a scenario of accretion occurring in bursts. The differentiation between C17O and SO2 suggests that they trace different physical components: C17O traces the densest and colder regions of the disc-envelope system, while SO2 may be associated with regions of higher temperature, such as accretion shocks. The lack of warm CH3OH emission suggests that there is no hot-core-like region around any of the sources and that the CH3OH column density averaged over the disc is low. △ Less

Submitted 30 April, 2019; originally announced April 2019.

Comments: 20 pages, 16 figures, 8 tables

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

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: 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: 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)

Protostellar half-life: new methodology and estimates

Authors: L. E. Kristensen, M. M. Dunham

Abstract: (Abridged) Protostellar systems evolve from prestellar cores, through the deeply embedded stage and then disk-dominated stage, before they end up on the main sequence. Knowing how much time a system spends in each stage is crucial for understanding how stars and associated planetary systems form, because a key constraint is the time available to form such systems. Equally important is understandin… ▽ More (Abridged) Protostellar systems evolve from prestellar cores, through the deeply embedded stage and then disk-dominated stage, before they end up on the main sequence. Knowing how much time a system spends in each stage is crucial for understanding how stars and associated planetary systems form, because a key constraint is the time available to form such systems. Equally important is understanding what the spread in these time scales is. The most commonly used method for inferring protostellar ages is to assume the lifetime of one evolutionary stage, and then scale this to the relative number of protostars in the other stages, i.e., assuming steady state. This method does not account for the underlying age distribution and apparent stochasticity of star formation, nor that relative populations are not in steady state. To overcome this, we propose a new scheme where the lifetime of each protostellar stage follows a distribution based on the formalism of sequential nuclear decay. The main assumptions are: Class 0 sources follow a straight path to Class III sources, the age distribution follows a binomial distribution, and the star-formation rate is constant. The results are that the half-life of Class 0, Class I, and Flat sources are (2.4+/-0.2)%, (4.4+/-0.3)%, and (4.3+/-0.4)% of the Class II half-life, respectively, which translates to 47+/-4, 88+/-7, and 87+/-8 kyr, respectively, for a Class II half-life of 2 Myr for protostars in the Gould Belt clouds with more than 100 protostars. The mean age of these clouds is 1.2+/-0.1 Myr, and the star formation rate is (8.3+/-0.5)x10^-4 Msun/yr. The critical parameters in arriving at these numbers are the assumed half-life of the Class II stage, and the assumption that the star-formation rate and half-lives are constant. This method presents a first step in moving from steady-state to non-steady-state solutions of protostellar populations. △ Less

Submitted 30 July, 2018; originally announced July 2018.

Comments: Accepted for publication in A&A

Journal ref: A&A 618, A158 (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)

Evidence for the start of planet formation in a young circumstellar disk

Authors: Daniel Harsono, Per Bjerkeli, Matthijs H. D. van der Wiel, Jon P. Ramsey, Luke T. Maud, Lars E. Kristensen, Jes K. Jørgensen

Abstract: The growth of dust grains in protoplanetary disks is a necessary first step towards planet formation. This growth has been inferred via observations of thermal dust emission towards mature protoplanetary systems (age >2 million years) with masses that are, on average, similar to Neptune3. In contrast, the majority of confirmed exoplanets are heavier than Neptune. Given that young protoplanetary di… ▽ More The growth of dust grains in protoplanetary disks is a necessary first step towards planet formation. This growth has been inferred via observations of thermal dust emission towards mature protoplanetary systems (age >2 million years) with masses that are, on average, similar to Neptune3. In contrast, the majority of confirmed exoplanets are heavier than Neptune. Given that young protoplanetary disks are more massive than their mature counterparts, this suggests that planet formation starts early, but evidence for grain growth that is spatially and temporally coincident with a massive reservoir in young disks remains scarce. Here, we report observations on a lack of emission of carbon monoxide isotopologues within the inner ~15 au of a very young (age ~100,000 years) disk around the Solar-type protostar TMC1A. By using the absence of spatially resolved molecular line emission to infer the gas and dust content of the disk, we conclude that shielding by millimeter-size grains is responsible for the lack of emission. This suggests that grain growth and millimeter-size dust grains can be spatially and temporally coincident with a mass reservoir sufficient for giant planet formation. Hence, planet formation starts during the earliest, embedded phases in the life of young stars. △ Less

Submitted 25 June, 2018; originally announced June 2018.

Comments: Accepted for publication in Nature Astronomy, 3 figures, 3 extended figures

Mass Assembly of Stellar Systems and their Evolution with the SMA - 1.3 mm Subcompact Data Release

Authors: Ian W. Stephens, Michael M. Dunham, Philip C. Myers, Riwaj Pokhrel, Tyler L. Bourke, Eduard I. Vorobyov, John J. Tobin, Sarah I. Sadavoy, Jaime E. Pineda, Stella S. R. Offner, Katherine I. Lee, Lars E. Kristensen, Jes K. Jørgensen, Alyssa A. Goodman, Héctor G. Arce, Mark Gurwell

Abstract: We present the Mass Assembly of Stellar Systems and their Evolution with the SMA (MASSES) survey, which uses the Submillimeter Array (SMA) interferometer to map the continuum and molecular lines for all 74 known Class 0/I protostellar systems in the Perseus molecular cloud. The primary goal of the survey is to observe an unbiased sample of young protostars in a single molecular cloud so that we ca… ▽ More We present the Mass Assembly of Stellar Systems and their Evolution with the SMA (MASSES) survey, which uses the Submillimeter Array (SMA) interferometer to map the continuum and molecular lines for all 74 known Class 0/I protostellar systems in the Perseus molecular cloud. The primary goal of the survey is to observe an unbiased sample of young protostars in a single molecular cloud so that we can characterize the evolution of protostars. This paper releases the MASSES 1.3 mm data from the subcompact configuration ($\sim$4$^{\prime\prime}$ or $\sim$1000 au resolution), which is the SMA's most compact array configuration. We release both $uv$ visibility data and imaged data for the spectral lines CO(2-1),$^{13}$CO(2-1), C$^{18}$O(2-1), and N$_2$D$^+$(3-2), as well as for the 1.3 mm continuum. We identify the tracers that are detected toward each source. We also show example images of continuum and CO(2-1) outflows, analyze C$^{18}$O(2-1) spectra, and present data from the SVS 13 star-forming region. The calculated envelope masses from the continuum show a decreasing trend with bolometric temperature (a proxy for age). Typical C$^{18}$O(2-1) linewidths are 1.45 km s$^{-1}$, which is higher than the C$^{18}$O linewidths detected toward Perseus filaments and cores. We find that N$_2$D$^+$(3-2) is significantly more likely to be detected toward younger protostars. We show that the protostars in SVS 13 are contained within filamentary structures as traced by C$^{18}$O(2-1) and N$_2$D$^+$(3-2). We also present the locations of SVS 13A's high velocity (absolute line-of-sight velocities $>$150 km s$^{-1}$) red and blue outflow components. Data can be downloaded from https://dataverse.harvard.edu/dataverse/MASSES . △ Less

Submitted 19 June, 2018; originally announced June 2018.

Comments: Accepted ApJS

CO in Protostars (COPS): $Herschel$-SPIRE Spectroscopy of Embedded Protostars

Authors: Yao-Lun Yang, Joel D. Green, Neal J. Evans II, Jeong-Eun Lee, Jes K Jorgensen, Lars E. Kristensen, Joseph C. Mottram, Gregory Herczeg, Agata Karska, Odysseas Dionatos, Edwin A. Bergin, Jeroen Bouwman, Ewine F. van Dishoeck, Tim A. van Kempen, Rebecca L. Larson, Umut A. Yildiz

Abstract: We present full spectral scans from 200-670$μ$m of 26 Class 0+I protostellar sources, obtained with $Herschel$-SPIRE, as part of the "COPS-SPIRE" Open Time program, complementary to the DIGIT and WISH Key programs. Based on our nearly continuous, line-free spectra from 200-670 $μ$m, the calculated bolometric luminosities ($L_{\rm bol}$) increase by 50% on average, and the bolometric temperatures (… ▽ More We present full spectral scans from 200-670$μ$m of 26 Class 0+I protostellar sources, obtained with $Herschel$-SPIRE, as part of the "COPS-SPIRE" Open Time program, complementary to the DIGIT and WISH Key programs. Based on our nearly continuous, line-free spectra from 200-670 $μ$m, the calculated bolometric luminosities ($L_{\rm bol}$) increase by 50% on average, and the bolometric temperatures ($T_{\rm bol}$) decrease by 10% on average, in comparison with the measurements without Herschel. Fifteen protostars have the same Class using $T_{\rm bol}$ and $L_{\rm bol}$/$L_{\rm submm}$. We identify rotational transitions of CO lines from J=4-3 to J=13-12, along with emission lines of $^{13}$CO, HCO$^+$, H$_{2}$O, and [CI]. The ratios of $^{12}$CO to $^{13}$CO indicate that $^{12}$CO emission remains optically thick for $J_{\rm up}$ < 13. We fit up to four components of temperature from the rotational diagram with flexible break points to separate the components. The distribution of rotational temperatures shows a primary population around 100 K with a secondary population at $\sim$350 K. We quantify the correlations of each line pair found in our dataset, and find the strength of correlation of CO lines decreases as the difference between $J$-level between two CO lines increases. The multiple origins of CO emission previously revealed by velocity-resolved profiles are consistent with this smooth distribution if each physical component contributes to a wide range of CO lines with significant overlap in the CO ladder. We investigate the spatial extent of CO emission and find that the morphology is more centrally peaked and less bipolar at high-$J$ lines. We find the CO emission observed with SPIRE related to outflows, which consists two components, the entrained gas and shocked gas, as revealed by our rotational diagram analysis as well as the studies with velocity-resolved CO emission. △ Less

Submitted 11 May, 2018; v1 submitted 2 May, 2018; originally announced May 2018.

Comments: 50 pages, 18 figures, accepted to ApJS. Revised for Table 6 and Figure 4

The ALMA-PILS survey: Complex nitriles towards IRAS 16293--2422

Authors: H. Calcutt, J. K. Jørgensen, H. S. P. Müller, L. E. Kristensen, A. Coutens, T. L. Bourke, R. T. Garrod, M. V. Persson, M. H. D. van der Wiel, E. F. van Dishoeck, S. F. Wampfler

Abstract: Complex organic molecules are readily detected in the inner regions of the gaseous envelopes of forming protostars. In particular, molecules that contain nitrogen are interesting due to the role nitrogen plays in the development of life and the compact scales such molecules have been found to trace around forming protostars. The goal of this work is to determine the inventory of one family of nitr… ▽ More Complex organic molecules are readily detected in the inner regions of the gaseous envelopes of forming protostars. In particular, molecules that contain nitrogen are interesting due to the role nitrogen plays in the development of life and the compact scales such molecules have been found to trace around forming protostars. The goal of this work is to determine the inventory of one family of nitrogen-bearing organic molecules, complex nitriles (molecules with a $-$CN functional group) towards two hot corino sources in the low-mass protostellar binary IRAS 16293$-$2422. This work explores the abundance differences between the two sources, the isotopic ratios, and the spatial extent derived from molecules containing the nitrile functional group. Using data from the Protostellar Interferometric Line Survey (PILS) obtained with ALMA we determine abundances and excitation temperatures for the detected nitriles. We also present a new method for determining the spatial structure of sources with high line density and large velocity gradients $-$ Velocity-corrected INtegrated emission (VINE) maps. We detect methyl cyanide (CH$_3$CN) as well as 5 of its isotopologues, including the detection of CHD$_2$CN which is the first detection in the ISM. We also detect ethyl cyanide (C$_2$H$_5$CN), vinyl cyanide (C$_2$H$_3$CN), and cyanoacetylene (HC$_3$N). We find that abundances are similar between IRAS 16293A and IRAS 16293B on small scales except for vinyl cyanide which is only detected towards the latter source. This suggests an important difference between the sources either in their evolutionary stage or warm-up timescales. We also detect a spatially double-peaked emission for the first time in molecular emission in the A source, suggesting that this source is showing structure related to a rotating toroid of material. △ Less

Submitted 24 April, 2018; originally announced April 2018.

Comments: A&A accepted 24th April 2018

Journal ref: A&A 616, A90 (2018)

Chemistry of a newly detected circumbinary disk in Ophiuchus

Authors: E. Artur de la Villarmois, L. E. Kristensen, J. K. Jørgensen, E. A. Bergin, C. Brinch, S. Frimann, D. Harsono, N. Sakai, S. Yamamoto

Abstract: (Abridged) Astronomers recently started discovering exoplanets around binary systems. Therefore, understanding the formation and evolution of circumbinary disks is crucial for a complete scenario of planet formation. The aim of this paper is to present the detection of a circumbinary disk around Oph-IRS67 and analyse its structure. We present high-angular-resolution (0.4", 60 AU) observations of C… ▽ More (Abridged) Astronomers recently started discovering exoplanets around binary systems. Therefore, understanding the formation and evolution of circumbinary disks is crucial for a complete scenario of planet formation. The aim of this paper is to present the detection of a circumbinary disk around Oph-IRS67 and analyse its structure. We present high-angular-resolution (0.4", 60 AU) observations of C17O, H13CO+ , C34S, SO2, C2H and c-C3H2 molecular transitions with ALMA at 0.8 mm. The spectrally and spatially resolved maps reveal the kinematics of the circumbinary disk as well as its chemistry. Molecular abundances are estimated using RADEX. The continuum emission reveals the presence of a circumbinary disk around the two sources. This disk has a diameter of ~620 AU and is well traced by C17O and H13CO+ emission. C2H and c-C3H2 trace a higher-density region which is spatially offset from the sources (~430 AU). Finally, SO2 shows compact emission around one of the sources, Oph-IRS67 B. The molecular transitions which trace the circumbinary disk are consistent with a Keplerian profile on disk scales (< 200 AU) and an infalling profile for envelope scales (> 200 AU). The Keplerian fit leads to a mass of 2.2 Msun. Inferred CO abundances w.r.t. H2 are comparable to the canonical ISM value of 2.7e-4. This study proves the first detection of the circumbinary disk associated with Oph-IRS67. The disk is chemically differentiated from the nearby high-density region. The lack of methanol emission suggests the extended disk dominates the mass budget in the inner- most regions of the protostellar envelope, generating a flat density profile where less material is exposed to high temperatures. Thus, complex organic molecules would be associated with lower column densities. Finally, Oph-IRS67 is a promising candidate for the detection of both circumstellar disks with higher-angular-resolution observations. △ Less

Submitted 26 February, 2018; originally announced February 2018.

Comments: 19 pages, 14 figures, 6 tables

Journal ref: A&A 614, A26 (2018)

The Herschel-PACS legacy of low-mass protostars: Properties of warm and hot gas and its origin in far-UV illuminated shocks

Authors: Agata Karska, Michael J. Kaufman, Lars E. Kristensen, Ewine F. van Dishoeck, Gregory J. Herczeg, Joseph C. Mottram, Łukasz Tychoniec, Johan E. Lindberg, Neal J. Evans II, Joel D. Green, Yao-Lun Yang, Antoine Gusdorf, Dominika Itrich, Natasza Siódmiak

Abstract: Recent observations from Herschel allow the identification of important mechanisms responsible for the heating of gas surrounding low-mass protostars and its subsequent cooling in the far-infrared (FIR). Shocks are routinely invoked to reproduce some properties of the far-IR spectra, but standard models fail to reproduce the emission from key molecules, e.g. H$_2$O. Here, we present the Herschel-P… ▽ More Recent observations from Herschel allow the identification of important mechanisms responsible for the heating of gas surrounding low-mass protostars and its subsequent cooling in the far-infrared (FIR). Shocks are routinely invoked to reproduce some properties of the far-IR spectra, but standard models fail to reproduce the emission from key molecules, e.g. H$_2$O. Here, we present the Herschel-PACS far-IR spectroscopy of 90 embedded low-mass protostars (Class 0/I). The Herschel-PACS spectral maps covering $\sim55-210$ $μ$m with a field-of-view of $\sim$50'' are used to quantify the gas excitation conditions and spatial extent using rotational transitions of H$_{2}$O, high-$J$ CO, and OH, as well as [O I] and [C II]. We confirm that a warm ($\sim$300 K) CO reservoir is ubiquitous and that a hotter component ($760\pm170$ K) is frequently detected around protostars. The line emission is extended beyond $\sim$1000 AU spatial scales in 40/90 objects, typically in molecular tracers in Class 0 and atomic tracers in Class I objects. High-velocity emission ($\gtrsim90$ km s$^{-1}$) is detected in only 10 sources in the [O I] line, suggesting that the bulk of [O I] arises from gas that is moving slower than typical jets. Line flux ratios show an excellent agreement with models of $C$-shocks illuminated by UV photons for pre-shock densities of $\sim$$10^5$ cm$^{-3}$ and UV fields 0.1-10 times the interstellar value. The far-IR molecular and atomic lines are a unique diagnostic of feedback from UV emission and shocks in envelopes of deeply embedded protostars. △ Less

Submitted 9 February, 2018; originally announced February 2018.

Comments: Accepted by ApJS, 55 pages, 170 figures

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)

Shocks and UV radiation around low-mass protostars: the Herschel-PACS legacy

Authors: Agata Karska, Michael J. Kaufman, Lars E. Kristensen, Ewine F. van Dishoeck

Abstract: Far-infrared spectroscopy reveals gas cooling and its underlying heating due to physical processes taking place in the surroundings of protostars. These processes are reflected in both the chemistry and excitation of abundant molecular species. Here, we present the Herschel-PACS far-IR spectroscopy of 90 embedded low-mass protostars from the WISH (van Dishoeck et al. 2011), DIGIT (Green et al. 201… ▽ More Far-infrared spectroscopy reveals gas cooling and its underlying heating due to physical processes taking place in the surroundings of protostars. These processes are reflected in both the chemistry and excitation of abundant molecular species. Here, we present the Herschel-PACS far-IR spectroscopy of 90 embedded low-mass protostars from the WISH (van Dishoeck et al. 2011), DIGIT (Green et al. 2013), and WILL surveys (Mottram et al. 2017). The $5\times5$ spectra covering the $\sim50''\times50''$ field-of-view include rotational transitions of CO, H$_2$O, and OH lines, as well as fine-structure [O I] and [C II] in the $\sim$50-200 $μ$m range. The CO rotational temperatures (for $J_\mathrm{u}\geq14)$ are typically $\sim$300 K, with some sources showing additional components with temperatures as high as $\sim$1000 K. The H$_2$O / CO and H$_2$O / OH flux ratios are low compared to stationary shock models, suggesting that UV photons may dissociate some H$_2$O and decrease its abundance. Comparison to C shock models illuminated by UV photons shows good agreement between the line emission and the models for pre-shock densities of $10^5$ cm$^{-3}$ and UV fields 0.1-10 times the interstellar value. The far-infrared molecular and atomic lines are the unique diagnostic of shocks and UV fields in deeply-embedded sources. △ Less

Submitted 8 November, 2017; originally announced November 2017.

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

Alignment Between Protostellar Outflows and Filamentary Structure

Authors: Ian W. Stephens, Michael M. Dunham, Philip C. Myers, Riwaj Pokhrel, Sarah I. Sadavoy, Eduard I. Vorobyov, John J. Tobin, Jaime E. Pineda, Stella S. R. Offner, Katherine I. Lee, Lars E. Kristensen, Jes K. Jørgensen, Alyssa A. Goodman, Tyler L. Bourke, Héctor G. Arce, Adele L. Plunkett

Abstract: We present new Submillimeter Array (SMA) observations of CO(2-1) outflows toward young, embedded protostars in the Perseus molecular cloud as part of the Mass Assembly of Stellar Systems and their Evolution with the SMA (MASSES) survey. For 57 Perseus protostars, we characterize the orientation of the outflow angles and compare them with the orientation of the local filaments as derived from… ▽ More We present new Submillimeter Array (SMA) observations of CO(2-1) outflows toward young, embedded protostars in the Perseus molecular cloud as part of the Mass Assembly of Stellar Systems and their Evolution with the SMA (MASSES) survey. For 57 Perseus protostars, we characterize the orientation of the outflow angles and compare them with the orientation of the local filaments as derived from $Herschel$ observations. We find that the relative angles between outflows and filaments are inconsistent with purely parallel or purely perpendicular distributions. Instead, the observed distribution of outflow-filament angles are more consistent with either randomly aligned angles or a mix of projected parallel and perpendicular angles. A mix of parallel and perpendicular angles requires perpendicular alignment to be more common by a factor of $\sim$3. Our results show that the observed distributions probably hold regardless of the protostar's multiplicity, age, or the host core's opacity. These observations indicate that the angular momentum axis of a protostar may be independent of the large-scale structure. We discuss the significance of independent protostellar rotation axes in the general picture of filament-based star formation. △ Less

Submitted 25 July, 2017; originally announced July 2017.

Comments: Accepted to ApJ

ALMA observations of dust polarization and molecular line emission from the Class 0 protostellar source Serpens SMM1

Authors: Charles L. H. Hull, Josep M. Girart, Łukasz Tychoniec, Ramprasad Rao, Paulo C. Cortés, Riwaj Pokhrel, Qizhou Zhang, Martin Houde, Michael M. Dunham, Lars E. Kristensen, Shih-Ping Lai, Zhi-Yun Li, Richard L. Plambeck

Abstract: We present high angular resolution dust polarization and molecular line observations carried out with the Atacama Large Millimeter/submillimeter Array (ALMA) toward the Class 0 protostar Serpens SMM1. By complementing these observations with new polarization observations from the Submillimeter Array (SMA) and archival data from the Combined Array for Research in Millimeter-wave Astronomy (CARMA) a… ▽ More We present high angular resolution dust polarization and molecular line observations carried out with the Atacama Large Millimeter/submillimeter Array (ALMA) toward the Class 0 protostar Serpens SMM1. By complementing these observations with new polarization observations from the Submillimeter Array (SMA) and archival data from the Combined Array for Research in Millimeter-wave Astronomy (CARMA) and the James Clerk Maxwell Telescopes (JCMT), we can compare the magnetic field orientations at different spatial scales. We find major changes in the magnetic field orientation between large (~0.1 pc) scales -- where the magnetic field is oriented E-W, perpendicular to the major axis of the dusty filament where SMM1 is embedded -- and the intermediate and small scales probed by CARMA (~1000 AU resolution), the SMA (~350 AU resolution), and ALMA (~140 AU resolution). The ALMA maps reveal that the redshifted lobe of the bipolar outflow is shaping the magnetic field in SMM1 on the southeast side of the source; however, on the northwestern side and elsewhere in the source, low velocity shocks may be causing the observed chaotic magnetic field pattern. High-spatial-resolution continuum and spectral-line observations also reveal a tight (~130 AU) protobinary system in SMM1-b, the eastern component of which is launching an extremely high-velocity, one-sided jet visible in both CO(2-1) and SiO(5-4); however, that jet does not appear to be shaping the magnetic field. These observations show that with the sensitivity and resolution of ALMA, we can now begin to understand the role that feedback (e.g., from protostellar outflows) plays in shaping the magnetic field in very young, star-forming sources like SMM1. △ Less

Submitted 12 July, 2017; originally announced July 2017.

Comments: 15 pages, 6 figures, 4 tables, 1 appendix. Accepted for publication in the Astrophysical Journal. Materials accessible in the online version of the (open-access) ApJ article include the FITS files used to make the ALMA image in Figure 1(d), and a full, machine-readable version of Table 3