Showing 1–50 of 87 results for author: Harsono, D

The Asymmetric Bipolar Fe II Jet and H2 Outflow of TMC1A Resolved with JWST's NIRSpec IFU

Authors: Korash Assani, Daniel Harsono, Jon Ramsey, Zhi-Yun Li, Per Bjerkeli, Klaus Pontoppidan, Łukasz Tychoniec, Hannah Calcutt, Lars Kristensen, Jes Jorgensen, Adele Plunkett, Martijn van Gelder, Logan Francis

Abstract: (abridged) Protostellar outflows exhibit large variations in their structure depending on the observed gas emission. This study analyzes the atomic jet and molecular outflow in the Class I protostar, TMC1A to characterize morphology and identify previously undetected spatial features with JWST's NIRSpec IFU. In addition to identifying a large number of Fe II and H2 lines, we have detected the bipo… ▽ More (abridged) Protostellar outflows exhibit large variations in their structure depending on the observed gas emission. This study analyzes the atomic jet and molecular outflow in the Class I protostar, TMC1A to characterize morphology and identify previously undetected spatial features with JWST's NIRSpec IFU. In addition to identifying a large number of Fe II and H2 lines, we have detected the bipolar Fe jet by revealing, for the first time, the presence of a red-shifted atomic jet. Similarly, the red-shifted component of the H2 slower wide-angle outflow is observed. Both Fe II and H2 red-shifted emission exhibit significantly lower flux densities compared to their blue-shifted counterparts. Additionally, we report the detection of a collimated high-velocity (100 km s-1), blue-shifted H2 outflow, suggesting the presence of a molecular jet in addition to the well-known wider angle low-velocity structure. The Fe II and H2 jets show multiple intensity peaks along the jet axis, which may be associated with ongoing or recent outburst events. In addition to the variation in their intensities, the H2 wide-angle outflow exhibits a "ring"-like structure. The blue-shifted H2 outflow also shows a left-right brightness asymmetry likely due to interactions with the surrounding ambient medium and molecular outflows. Using the Fe II line ratios, the extinction along the atomic jet is estimated to be between Av = 10-30 on the blue-shifted side, with a trend of decreasing extinction with distance from the protostar. A similar Av is found for the red-shifted side, supporting the argument for an intrinsic red-blue outflow lobe asymmetry rather than environmental effects such as extinction. This intrinsic difference revealed by the unprecedented sensitivity of JWST, suggests that younger outflows already exhibit the red-blue side asymmetry more commonly observed towards jets associated with Class II disks. △ Less

Submitted 28 April, 2024; originally announced April 2024.

Comments: 16 pages, 11 figures

Emergence of high-mass stars in complex fiber networks (EMERGE). I. Early ALMA Survey: observations and massive data reduction

Authors: A. Hacar, A. Socci, F. Bonanomi, D. Petry, M. Tafalla, D. Harsono, J. Forbrich, J. Alves, J. Grossschedl, J. R. Goicoechea, J. Pety, A. Burkert, G. X. Li

Abstract: (Abridged) Recent molecular surveys have revealed a rich gas organization of sonic-like fibers in all kind of environments prior to the formation of low- and high-mass stars. This paper introduces the EMERGE project aiming to investigate whether complex fiber arrangements could explain the origin of high-mass stars and clusters. We analyzed the EMERGE Early ALMA Survey including 7 star-forming reg… ▽ More (Abridged) Recent molecular surveys have revealed a rich gas organization of sonic-like fibers in all kind of environments prior to the formation of low- and high-mass stars. This paper introduces the EMERGE project aiming to investigate whether complex fiber arrangements could explain the origin of high-mass stars and clusters. We analyzed the EMERGE Early ALMA Survey including 7 star-forming regions in Orion (OMC-1/2/3/4 South, L1641N, NGC2023, and Flame Nebula) homogeneously surveyed in both molecular lines (N$_2$H$^+$ J=1-0, HNC J=1-0, plus HC3N J=10-9) and 3mm-continuum using a combination of interferometric ALMA mosaics and IRAM-30m single-dish (SD) maps. Based on our low-resolution (SD) observations, we describe the global properties of our sample covering a wide range of physical conditions including low-, intermediate, and high-mass star-forming regions in different evolutionary stages. Their comparison with ancillary YSO catalogs denotes N$_2$H$^+$ as the best proxy for the dense, star-forming gas in our targets showing a constant star formation efficiency and a fast time evolution of <1 Myr. While apparently clumpy and filamentary in our SD data, all targets show a much more complex fibrous substructure at the enhanced resolution of our ALMA+IRAM-30m maps. A large number of filamentary features at sub-parsec scales are clearly recognized in the high-density gas traced by N$_2$H$^+$ directly connected to the formation of individual protostars. This complex gas organization appears to extend further into the more diffuse gas traced by HNC. This paper presents the EMERGE Early ALMA survey including a first data release of continuum maps and spectral products for this project to be analysed in future papers of this series. A first look at these results illustrates the need of advanced data combination techniques to investigate the intrinsic multi-scale, gas structure of the ISM. △ Less

Submitted 12 March, 2024; originally announced March 2024.

Comments: 37 pages, 28 fgures, Accepted by A&A

Small and Large Dust Cavities in Disks around mid-M Stars in Taurus

Authors: Yangfan Shi, Feng Long, Gregory J. Herczeg, Daniel Harsono, Yao Liu, Paola Pinilla, Enrico Ragusa, Doug Johnstone, Xue-Ning Bai, Ilaria Pascucci, Carlo F. Manara, Gijs D. Mulders, Lucas A. Cieza

Abstract: High-angular resolution imaging by ALMA has revealed the near-universality and diversity of substructures in protoplanetary disks. However, disks around M-type pre-main-sequence stars are still poorly sampled, despite the prevalence of M-dwarfs in the galaxy. Here we present high-resolution (~50 mas, 8 au) ALMA Band 6 observations of six disks around mid-M stars in Taurus. We detect dust continuum… ▽ More High-angular resolution imaging by ALMA has revealed the near-universality and diversity of substructures in protoplanetary disks. However, disks around M-type pre-main-sequence stars are still poorly sampled, despite the prevalence of M-dwarfs in the galaxy. Here we present high-resolution (~50 mas, 8 au) ALMA Band 6 observations of six disks around mid-M stars in Taurus. We detect dust continuum emission in all six disks, 12CO in five disks, and 13CO line in two disks. The size ratios between gas and dust disks range from 1.6 to 5.1. The ratio of about 5 for 2M0436 and 2M0450 indicates efficient dust radial drift. Four disks show rings and cavities and two disks are smooth. The cavity sizes occupy a wide range: 60 au for 2M0412, and ~10 au for 2M0434, 2M0436 and 2M0508. Detailed visibility modeling indicates that small cavities of 1.7 and 5.7 au may hide in the two smooth disks 2M0450 and CIDA 12. We perform radiative transfer fitting of the infrared SEDs to constrain the cavity sizes, finding that micron-sized dust grains may have smaller cavities than millimeter grains. Planet-disk interactions are the preferred explanation to produce the large 60 au cavity, while other physics could be responsible for the three ~10 au cavities under current observations and theories. Currently, disks around mid-to-late M stars in Taurus show a higher detection frequency of cavities than earlier type stars, although a more complete sample is needed to evaluate any dependence of substructure on stellar mass. △ Less

Submitted 28 February, 2024; originally announced February 2024.

Comments: 25 pages, 6 figures. Accepted for publication in ApJ

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

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

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

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

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

Dual-Band Observations of the Asymmetric Ring around CIDA 9A: Dead or Alive?

Authors: Daniel Harsono, Feng Long, Paola Pinilla, Alessia A. Rota, Carlo F. Manara, Gregory J. Herczeg, Doug Johnstone, Giovanni Rosotti, Giuseppe Lodato, Francois Menard, Marco Tazzari, Yangfan Shi

Abstract: While the most exciting explanation of the observed dust asymmetries in protoplanetary disks is the presence of protoplanets, other mechanisms can also form the dust features. This paper presents dual-wavelength Atacama Large Millimeter/submillimeter Array (ALMA) observations of a large asymmetric dusty ring around the M-type star CIDA 9A. We detect a dust asymmetry in both 1.3 mm and 3.1 mm data.… ▽ More While the most exciting explanation of the observed dust asymmetries in protoplanetary disks is the presence of protoplanets, other mechanisms can also form the dust features. This paper presents dual-wavelength Atacama Large Millimeter/submillimeter Array (ALMA) observations of a large asymmetric dusty ring around the M-type star CIDA 9A. We detect a dust asymmetry in both 1.3 mm and 3.1 mm data. To characterize the asymmetric structure, a parametric model is used to fit the observed visibilities. We report a tentative azimuthal shift of the dust emission peaks between the observations at the two wavelengths. This shift is consistent with a dust trap caused by a vortex, which may be formed by an embedded protoplanet or other hydrodynamical instabilities, such as a dead zone. Deep high-spatial observations of dust and molecular gas are needed to constrain the mechanisms that formed the observed millimeter cavity and dust asymmetry in the protoplanetary disk around CIDA 9A. △ Less

Submitted 17 October, 2023; originally announced October 2023.

Comments: Accepted for publication in ApJ, 12 pages, 7 figures with appendix

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

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

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

Submitted 14 September, 2023; originally announced September 2023.

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

A high HDO/H$_{2}$O ratio in the Class I protostar L1551 IRS5

Authors: Audrey Andreu, Audrey Coutens, Fernando Cruz-Sáenz de Miera, Nicolas Houry, Jes K. Jørgensen, Ágnes Kóspál, Daniel Harsono

Abstract: Water is a very abundant molecule in star-forming regions. Its deuterium fractionation is an important tool for understanding its formation and evolution during the star and planet formation processes. While the HDO/H$_2$O ratio has been determined toward several Class 0 protostars and comets, the number of studies toward Class I protostars is limited. We aim to study the water deuteration toward… ▽ More Water is a very abundant molecule in star-forming regions. Its deuterium fractionation is an important tool for understanding its formation and evolution during the star and planet formation processes. While the HDO/H$_2$O ratio has been determined toward several Class 0 protostars and comets, the number of studies toward Class I protostars is limited. We aim to study the water deuteration toward the Class I binary protostar L1551 IRS5 and to investigate the effect of evolutionary stage and environment on variations in the water D/H ratio. Observations were made using the NOEMA interferometer. The HDO 3$_{1,2}$-2$_{2,1}$ transition at 225.9 GHz and the H$_2^{18}$O 3$_{1,3}$-2$_{2,0}$ transition at 203.4 GHz were covered with a spatial resolution of 0.5'' $\times$ 0.8'', while the HDO 4$_{2,2}$-4$_{2,3}$ transition at 143.7 GHz was observed with a resolution of 2.0'' $\times$ 2.5''. We used both LTE and non-LTE models. The three transitions are detected. The line profiles display two peaks, one at $\sim$6 km s$^{-1}$ and one at $\sim$9 km s$^{-1}$. We derive an HDO/H$_2$O ratio of (2.1 $\pm$ 0.8) $\times$ 10$^{-3}$ for the redshifted component and a lower limit of $>$ 0.3 $\times$ 10$^{-3}$ for the blueshifted component due to the blending with the redshifted CH$_3$OCH$_3$ emission. The HDO/H$_2$O in L1551 IRS5 is similar to the ratios in isolated Class 0 sources and to the Class I V883 Ori, while it is significantly higher than in the clustered Class 0 sources and the comets. This suggests that the chemistry of protostars in low source densities clouds share more similarities with the isolated sources than the protostars of very dense clusters. If Class 0 protostars with few sources around and isolated Class 0 objects are comparable in the HDO/H$_2$O ratio, it would mean that there is little water reprocessing from the Class 0 to Class I protostellar stage. △ Less

Submitted 4 September, 2023; originally announced September 2023.

Comments: Accepted in A&A, 10 pages, 8 figures

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

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)

JWST Peers into the Class I Protostar TMC1A: Atomic Jet and Spatially Resolved Dissociative Shock Region

Authors: Daniel Harsono, Per Bjerkeli, Jon Ramsey, Klaus Pontoppidan, Lars Kristensen, Jes Jørgensen, Hannah Calcutt, Zhi-Yun Li, Adele Plunkett

Abstract: Outflows and winds launched from young stars play a crucial role in the evolution of protostars and the early stages of planet formation. However, the specific details of the mechanism behind these phenomena, including how they affect the protoplanetary disk structure, are still debated. We present {\it JWST} NIRSpec Integral Field Unit (IFU) observations of atomic and H$_2$ lines from 1 -- 5.1… ▽ More Outflows and winds launched from young stars play a crucial role in the evolution of protostars and the early stages of planet formation. However, the specific details of the mechanism behind these phenomena, including how they affect the protoplanetary disk structure, are still debated. We present {\it JWST} NIRSpec Integral Field Unit (IFU) observations of atomic and H$_2$ lines from 1 -- 5.1 $μ$m toward the low-mass protostar TMC1A. For the first time, a collimated atomic jet is detected from TMC1A in the [Fe II] line at 1.644 $μ$m along with corresponding extended H$_2$ 2.12 $μ$m emission. Towards the protostar, we detected spectrally broad H I and He I emissions with velocities up to 300 km/s that can be explained by a combination of protostellar accretion and a wide-angle wind. The 2$μ$m continuum dust emission, H I, He I, and O I all show emission from the illuminated outflow cavity wall and scattered line emission. These observations demonstrate the potential of {\it JWST} to characterize and reveal new information about the hot inner regions of nearby protostars. In this case, a previously undetected atomic wind and ionized jet in a well-known outflow. △ Less

Submitted 14 June, 2023; originally announced June 2023.

Comments: Accepted for publication on ApJL, 10 pages, 6 figures including appendix

Triple spiral arms of a triple protostar system imaged in molecular lines

Authors: Jeong-Eun Lee, Tomoaki Matsumoto, Hyun-Jeong Kim, Seokho Lee, Daniel Harsono, Jaehan Bae, Neal J. Evans II, Shu-ichiro Inutsuka, Minho Choi, Ken'ichi Tatematsu, Jae-Joon Lee, Dan Jaffe

Abstract: Most stars form in multiple star systems. For a better understanding of their formation processes, it is important to resolve the individual protostellar components and the surrounding envelope and disk material at the earliest possible formation epoch because the formation history can be lost in a few orbital timescales. Here we present the ALMA observational results of a young multiple protostel… ▽ More Most stars form in multiple star systems. For a better understanding of their formation processes, it is important to resolve the individual protostellar components and the surrounding envelope and disk material at the earliest possible formation epoch because the formation history can be lost in a few orbital timescales. Here we present the ALMA observational results of a young multiple protostellar system, IRAS 04239+2436, where three well-developed large spiral arms were detected in the shocked SO emission. Along the most conspicuous arm, the accretion streamer was also detected in the SO$_2$ emission. The observational results are complemented by numerical magneto-hydrodynamic simulations, where those large arms only appear in magnetically weakened clouds. The numerical simulations also suggest that the large triple spiral arms are the result of gravitational interactions between compact triple protostars and the turbulent infalling envelope. △ Less

Submitted 10 June, 2023; originally announced June 2023.

Comments: Accepted for publication in ApJ

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

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

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

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

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

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

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

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

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

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

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

A Large Double-ring Disk around the Taurus M Dwarf J04124068+2438157

Authors: Feng Long, Bin B. Ren, Nicole L. Wallack, Daniel Harsono, Gregory J. Herczeg, Paola Pinilla, Dimitri Mawet, Michael C. Liu, Sean M. Andrews, Xue-Ning Bai, Sylvie Cabrit, Lucas A. Cieza, Doug Johnstone, Jarron M. Leisenring, Giuseppe Lodato, Yao Liu, Carlo F. Manara, Gijs D. Mulders, Enrico Ragusa, Steph Sallum, Yangfan Shi, Marco Tazzari, Taichi Uyama, Kevin Wagner, David J. Wilner , et al. (1 additional authors not shown)

Abstract: Planet formation imprints signatures on the physical structures of disks. In this paper, we present high-resolution ($\sim$50 mas, 8 au) Atacama Large Millimeter/submillimeter Array (ALMA) observations of 1.3 mm dust continuum and CO line emission toward the disk around the M3.5 star 2MASS J04124068+2438157. The dust disk consists only of two narrow rings at radial distances of 0.47 and 0.78 arcse… ▽ More Planet formation imprints signatures on the physical structures of disks. In this paper, we present high-resolution ($\sim$50 mas, 8 au) Atacama Large Millimeter/submillimeter Array (ALMA) observations of 1.3 mm dust continuum and CO line emission toward the disk around the M3.5 star 2MASS J04124068+2438157. The dust disk consists only of two narrow rings at radial distances of 0.47 and 0.78 arcsec ($\sim$70 and 116 au), with Gaussian $σ$ widths of 5.6 and 8.5 au, respectively. The width of the outer ring is smaller than the estimated pressure scale height by $\sim25\%$, suggesting dust trapping in a radial pressure bump. The dust disk size, set by the location of the outermost ring, is significantly larger (by $3σ$) than other disks with similar millimeter luminosity, which can be explained by an early formation of local pressure bump to stop radial drift of millimeter dust grains. After considering the disk's physical structure and accretion properties, we prefer planet--disk interaction over dead zone or photoevaporation models to explain the observed dust disk morphology. We carry out high-contrast imaging at $L'$ band using Keck/NIRC2 to search for potential young planets, but do not identify any source above $5σ$. Within the dust gap between the two rings, we reach a contrast level of $\sim$7 mag, constraining the possible planet below $\sim$2--4 $M_{\rm Jup}$. Analyses of the gap/ring properties suggest a $\sim$Saturn mass planet at $\sim$90 au is likely responsible for the formation of the outer ring, which can be potentially revealed with JWST. △ Less

Submitted 25 March, 2023; originally announced March 2023.

Comments: 15 pages, 5 figures. Accepted for publication in ApJ

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

Constraints on the non-thermal desorption of methanol in the cold core LDN 429-C

Authors: A. Taillard, V. WakelaM, P. Gratier, E. Dartois, M. Chabot, J. A. Noble, J. V. Keane, A. C. A. Boogert, D. Harsono

Abstract: Cold cores are an early step of star formation, characterized by densities > 10$^4$ cm$^{-3}$, low temperatures (< 15 K), and very low external UV radiation. We investigate the physico-chemical processes at play to tracing the origin of molecules that are predominantly formed via reactions on dust grain surfaces. We observed the cold core LDN 429-C with the NOEMA interferometer and the IRAM 30m si… ▽ More Cold cores are an early step of star formation, characterized by densities > 10$^4$ cm$^{-3}$, low temperatures (< 15 K), and very low external UV radiation. We investigate the physico-chemical processes at play to tracing the origin of molecules that are predominantly formed via reactions on dust grain surfaces. We observed the cold core LDN 429-C with the NOEMA interferometer and the IRAM 30m single dish telescope in order to obtain the gas-phase abundances of key species, including CO and CH$_3$OH. Comparing the observed gas phase of methanol to its solid phase previously observed with Spitzer allows us to put quantitative constraints on the efficiency of the non-thermal desorption of this species. With physical parameters determined from available Herschel data, we computed abundance maps of 11 detected molecules with a non-local thermal equilibrium radiative transfer model. These observations allowed us to probe the molecular abundances as a function of density and visual extinction, with the variation in temperature being restrained between 12 and 18 K. We then compared the observed abundances to the predictions of the Nautilus astrochemical model. We find that all molecules have lower abundances at high densities and visual extinctions with respect to lower density regions, except for methanol. Comparing these observations with a grid of chemical models based on the local physical conditions, we were able to reproduce these observations, allowing only the parameter time to vary. Comparing the observed gas-phase abundance of methanol with previous measurements of the methanol ice, we estimate a non-thermal desorption efficiency between 0.002% and 0.09%, increasing with density. The apparent increase in the desorption efficiency cannot be reproduced by our model unless the yield of cosmic-ray sputtering is altered due to the ice composition varying as a function of density. △ Less

Submitted 3 January, 2023; originally announced January 2023.

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)

The young embedded disk L1527 IRS: constraints on the water snowline and cosmic ray ionization rate from HCO+ observations

Authors: Merel L. R. van 't Hoff, Margot Leemker, John J. Tobin, Daniel Harsono, Jes K. Jørgensen, Edwin A. Bergin

Abstract: The water snowline in circumstellar disks is a crucial component in planet formation, but direct observational constraints on its location remain sparse due to the difficulty of observing water in both young embedded and mature protoplanetary disks. Chemical imaging provides an alternative route to locate the snowline, and HCO$^+$ isotopologues have been shown to be good tracers in protostellar en… ▽ More The water snowline in circumstellar disks is a crucial component in planet formation, but direct observational constraints on its location remain sparse due to the difficulty of observing water in both young embedded and mature protoplanetary disks. Chemical imaging provides an alternative route to locate the snowline, and HCO$^+$ isotopologues have been shown to be good tracers in protostellar envelopes and Herbig disks. Here we present $\sim$0.5$^{\prime\prime}$ resolution ($\sim$35 au radius) Atacama Large Millimeter/submillimeter Array (ALMA) observations of HCO$^+$ $J=4-3$ and H$^{13}$CO$^+$ $J=3-2$ toward the young (Class 0/I) disk L1527 IRS. Using a source-specific physical model with the midplane snowline at 3.4 au and a small chemical network, we are able to reproduce the HCO$^+$ and H$^{13}$CO$^+$ emission, but for HCO$^+$ only when the cosmic ray ionization rate is lowered to $10^{-18}$ s$^{-1}$. Even though the observations are not sensitive to the expected HCO$^+$ abundance drop across the snowline, the reduction in HCO$^+$ above the snow surface and the global temperature structure allow us to constrain a snowline location between 1.8 and 4.1 au. Deep observations are required to eliminate the envelope contribution to the emission and to derive more stringent constraints on the snowline location. Locating the snowline in young disks directly with observations of H$_2$O isotopologues may therefore still be an alternative option. With a direct snowline measurement, HCO$^+$ will be able to provide constraints on the ionization rate. △ Less

Submitted 18 April, 2022; originally announced April 2022.

Comments: Accepted for publication in ApJ, 15 pages, 6 figures and appendix

Disk Evolution Study Through Imaging of Nearby Young Stars (DESTINYS): A Panchromatic View of DO Tau's Complex Kilo-au Environment

Authors: Jane Huang, Christian Ginski, Myriam Benisty, Bin Ren, Alexander J. Bohn, Élodie Choquet, Karin I. Öberg, Álvaro Ribas, Jaehan Bae, Edwin A. Bergin, Til Birnstiel, Yann Boehler, Stefano Facchini, Daniel Harsono, Michiel Hogerheijde, Feng Long, Carlo F. Manara, François Ménard, Paola Pinilla, Christophe Pinte, Christian Rab, Jonathan P. Williams, Alice Zurlo

Abstract: While protoplanetary disks are often treated as isolated systems in planet formation models, observations increasingly suggest that vigorous interactions between Class II disks and their environments are not rare. DO Tau is a T Tauri star that has previously been hypothesized to have undergone a close encounter with the HV Tau system. As part of the DESTINYS ESO Large Programme, we present new VLT… ▽ More While protoplanetary disks are often treated as isolated systems in planet formation models, observations increasingly suggest that vigorous interactions between Class II disks and their environments are not rare. DO Tau is a T Tauri star that has previously been hypothesized to have undergone a close encounter with the HV Tau system. As part of the DESTINYS ESO Large Programme, we present new VLT/SPHERE polarimetric observations of DO Tau and combine them with archival HST scattered light images and ALMA observations of CO isotopologues and CS to map a network of complex structures. The SPHERE and ALMA observations show that the circumstellar disk is connected to arms extending out to several hundred au. HST and ALMA also reveal stream-like structures northeast of DO Tau, some of which are at least several thousand au long. These streams appear not to be gravitationally bound to DO Tau, and comparisons with previous Herschel far-IR observations suggest that the streams are part of a bridge-like structure connecting DO Tau and HV Tau. We also detect a fainter redshifted counterpart to a previously known blueshifted CO outflow. While some of DO Tau's complex structures could be attributed to a recent disk-disk encounter, they might be explained alternatively by interactions with remnant material from the star formation process. These panchromatic observations of DO Tau highlight the need to contextualize the evolution of Class II disks by examining processes occurring over a wide range of size scales. △ Less

Submitted 8 May, 2022; v1 submitted 4 April, 2022; originally announced April 2022.

Comments: 42 pages, 20 figures, accepted by ApJ, reduced data available at https://zenodo.org/record/6408903 (typos from v1 fixed)

Gas Disk Sizes from CO Line Observations: A Test of Angular Momentum Evolution

Authors: Feng Long, Sean M. Andrews, Giovanni Rosotti, Daniel Harsono, Paola Pinilla, David J. Wilner, Karin I Öberg, Richard Teague, Leon Trapman, Benoît Tabone

Abstract: The size of a disk encodes important information about its evolution. Combining new Submillimeter Array (SMA) observations with archival Atacama Large Millimeter Array (ALMA) data, we analyze mm continuum and CO emission line sizes for a sample of 44 protoplanetary disks around stars with masses of 0.15--2\,$M_{\odot}$ in several nearby star-forming regions. Sizes measured from $^{12}$CO line emis… ▽ More The size of a disk encodes important information about its evolution. Combining new Submillimeter Array (SMA) observations with archival Atacama Large Millimeter Array (ALMA) data, we analyze mm continuum and CO emission line sizes for a sample of 44 protoplanetary disks around stars with masses of 0.15--2\,$M_{\odot}$ in several nearby star-forming regions. Sizes measured from $^{12}$CO line emission span from 50 to 1000\,au. This range could be explained by viscous evolution models with different $α$ values (mostly of $10^{-4}-10^{-3}$) and/or a spread of initial conditions. The CO sizes for most disks are also consistent with MHD wind models that directly remove disk angular momentum, but very large initial disk sizes would be required to account for the very extended CO disks in the sample. As no CO size evolution is observed across stellar ages of 0.5--20\,Myr in this sample, determining the dominant mechanism of disk evolution will require a more complete sample for both younger and more evolved systems. We find that the CO emission is universally more extended than the continuum emission by an average factor of $2.9\pm1.2$. The ratio of the CO to continuum sizes does not show any trend with stellar mass, mm continuum luminosity, or the properties of substructures. The GO Tau disk has the most extended CO emission in this sample, with an extreme CO to continuum size ratio of 7.6. Seven additional disks in the sample show high size ratios ($\gtrsim4$) that we interpret as clear signs of substantial radial drift. △ Less

Submitted 30 March, 2022; originally announced March 2022.

Comments: Accepted for publication in ApJ

Planet Formation Theory in the Era of ALMA and Kepler: from Pebbles to Exoplanets

Authors: Joanna Drazkowska, Bertram Bitsch, Michiel Lambrechts, Gijs D. Mulders, Daniel Harsono, Allona Vazan, Beibei Liu, Chris W. Ormel, Katherine Kretke, Alessandro Morbidelli

Abstract: Our understanding of planet formation has been rapidly evolving in recent years. The classical planet formation theory, developed when the only known planetary system was our own Solar System, has been revised to account for the observed diversity of the exoplanetary systems. At the same time, the increasing observational capabilities of the young stars and their surrounding disks bring new constr… ▽ More Our understanding of planet formation has been rapidly evolving in recent years. The classical planet formation theory, developed when the only known planetary system was our own Solar System, has been revised to account for the observed diversity of the exoplanetary systems. At the same time, the increasing observational capabilities of the young stars and their surrounding disks bring new constraints on the planet formation process. In this chapter, we summarize the new information derived from the exoplanets population and the circumstellar disks observations. We present the new developments in planet formation theory, from dust evolution to the growth of planetary cores by accretion of planetesimals, pebbles, and gas. We review the state-of-the-art models for the formation of diverse planetary systems, including the population synthesis approach which is necessary to compare theoretical model outcomes to the exoplanet population. We emphasize that the planet formation process may not be spatially uniform in the disk and there are preferential locations for the formation of planetesimals and planets. Outside of these locations, a significant fraction of solids is not growing past the pebble-sizes. The reservoir of pebbles plays an important role in the growth of planetary cores in the pebble accretion process. The timescale of the emergence of massive planetary cores is an important aspect of the present models and it is likely that the cores within one disk form at different times. In addition, there is growing evidence that the first planetary cores start forming early, during the circumstellar disk buildup process. △ Less

Submitted 24 May, 2023; v1 submitted 18 March, 2022; originally announced March 2022.

Comments: 38 pages, 13 figures, review chapter accepted to Protostars and Planets VII, editors: Shu-ichiro Inutsuka, Yuri Aikawa, Takayuki Muto, Kengo Tomida, and Motohide Tamura, corrected label in Fig. 1

Tracing pebble drift and trapping using radial carbon depletion profiles in protoplanetary disks

Authors: J. A. Sturm, M. K. McClure, D. Harsono, S. Facchini, F. Long, M. Kama, E. A. Bergin, E. F. van Dishoeck

Abstract: The composition of planets may be largely determined by the chemical processing and accretion of icy pebbles in protoplanetary disks. Recent observations of protoplanetary disks hint at wide-spread depletion of gaseous carbon. The missing volatile carbon is likely frozen in CO and/or CO$_2$ ice on grains and locked into the disk through pebble trapping in pressure bumps or planetesimals. We presen… ▽ More The composition of planets may be largely determined by the chemical processing and accretion of icy pebbles in protoplanetary disks. Recent observations of protoplanetary disks hint at wide-spread depletion of gaseous carbon. The missing volatile carbon is likely frozen in CO and/or CO$_2$ ice on grains and locked into the disk through pebble trapping in pressure bumps or planetesimals. We present the results of the first successful ACA (Atacama Compact Array) [C I] $J$ = 1-0 mini-survey of seven protoplanetary disks. Using tailored azimuthally symmetric DALI (Dust And LInes) thermo-chemical disk models, supported by the [C I] $J$ = 1-0 and resolved CO isotopologue data, we determine the system-averaged elemental volatile carbon abundance in the outer disk of three sources. Six out of seven sources are detected in [C I] $J$ = 1-0 with ACA, four of which show a distinct disk component. Based on the modeling we find severe cold gaseous carbon depletion in the outer disk of DL Tau and moderate depletion in the outer disks of DR Tau and DO Tau. Combining the outer and inner disk carbon abundances, we demonstrate definitive evidence for radial drift in the disk of DL Tau, where the existence of multiple dust rings points to either short lived or leaky dust traps. We find dust locking in the compact and smooth disks of DO Tau and DR Tau, hinting at unresolved dust substructure. Comparing our results with stars of different ages and luminosities, we identify an observational evolutionary trend in gaseous carbon depletion that is consistent with dynamical models of CO depletion processes. Transport efficiency of solids in protoplanetary disks can significantly differ from what we expect based on the current resolved substructure in the continuum observations. This has important implications for our understanding of the impact of radial drift and pebble accretion on planetary compositions. △ Less

Submitted 11 January, 2022; originally announced January 2022.

Comments: 23 pages, 19 figures, accepted to A&A

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

Observational constraints on disc sizes in protoplanetary discs in multiple systems in the Taurus region. II. Gas disc sizes

Authors: A. A. Rota, C. F. Manara, A. Miotello, G. Lodato, S. Facchini, M. Koutoulaki, G. Herczeg, F. Long, M. Tazzari, S. Cabrit, D. Harsono, F. Menard, P. Pinilla, G. van der Plas, E. Ragusa, H. -W. Yen

Abstract: The formation of multiple stellar systems is a natural by-product of the star-formation process, and its impact on the properties of protoplanetary discs and on the formation of planets is still to be fully understood. To date, no detailed uniform study of the gas emission from a sample of protoplanetary discs around multiple stellar systems has been performed. Here we analyse new ALMA observation… ▽ More The formation of multiple stellar systems is a natural by-product of the star-formation process, and its impact on the properties of protoplanetary discs and on the formation of planets is still to be fully understood. To date, no detailed uniform study of the gas emission from a sample of protoplanetary discs around multiple stellar systems has been performed. Here we analyse new ALMA observations at a $\sim$21 au resolution of the molecular CO gas emission targeting discs in eight multiple stellar systems in the Taurus star-forming regions. $^{12}$CO gas emission is detected around all primaries and in seven companions. With these data, we estimate the inclination and the position angle for all primary discs and for five secondary or tertiary discs, and measure the gas disc radii of these objects with a cumulative flux technique on the spatially resolved zeroth moment images. When considering the radius including 95\% of the flux as a metric, the estimated gas disc size in multiple stellar systems is found to be on average $\sim 4.2$ times larger than the dust disc size. This ratio is higher than what was recently found in a population of more isolated and single systems. On the contrary, when considering the radius including 68\% of the flux, no difference between multiple and single discs is found in the distribution of ratios. This discrepancy is due to the sharp truncation of the outer dusty disc observed in multiple stellar systems. The measured gas disc sizes are consistent with tidal truncation models in multiple stellar systems assuming eccentricities of $\sim0.15$-$0.5$, as expected in typical binary systems. △ Less

Submitted 10 January, 2022; originally announced January 2022.

Comments: 25 pages, 29 figures, 9 tables. Accepted for publication in A&A

Journal ref: A&A 662, A121 (2022)

A cold accretion flow onto one component of a multiple protostellar system

Authors: Nadia M. Murillo, Ewine F. van Dishoeck, Alvaro Hacar, Daniel Harsono, Jes K. Jørgensen

Abstract: Context: Gas accretion flows transport material from the cloud core onto the protostar. In multiple protostellar systems, it is not clear if the delivery mechanism is preferential or evenly distributed among the components. Aims: Gas accretion flows within IRAS16293 is explored out to 6000 AU. Methods: ALMA Band 3 observations of low-$J$ transitions of HNC, cyanopolyynes (HC$_3$N, HC$_5$N), an… ▽ More Context: Gas accretion flows transport material from the cloud core onto the protostar. In multiple protostellar systems, it is not clear if the delivery mechanism is preferential or evenly distributed among the components. Aims: Gas accretion flows within IRAS16293 is explored out to 6000 AU. Methods: ALMA Band 3 observations of low-$J$ transitions of HNC, cyanopolyynes (HC$_3$N, HC$_5$N), and N$_2$H$^+$ are used to probe the cloud core structure at ~100 AU resolution. Additional Band 3 archival data provide low-$J$ HCN and SiO lines. These data are compared with the corresponding higher-$J$ lines from the PILS Band 7 data for excitation analysis. The HNC/HCN ratio is used as a temperature tracer. Results: The low-$J$ transitions of HC$_3$N, HC$_5$N, HNC and N$_2$H$^+$ trace extended and elongated structures from 6000 AU down to ~100 AU, without accompanying dust continuum emission. Two structures are identified: one traces a flow that is likely accreting toward the most luminous component of the system IRAS16293 A. Temperatures inferred from the HCN/HNC ratio suggest that the gas in this flow is cold, between 10 and 30 K. The other structure is part of an UV-irradiated cavity wall entrained by one of the outflows. The two outflows driven by IRAS16293 A present different molecular gas distributions. Conclusions: Accretion of cold gas is seen from 6000 AU scales onto IRAS16293 A, but not onto source B, indicates that cloud core material accretion is competitive due to feedback onto a dominant component in an embedded multiple protostellar system. The preferential delivery of material could explain the higher luminosity and multiplicity of source A compared to source B. The results of this work demonstrate that several different molecular species, and multiple transitions of each species, are needed to confirm and characterize accretion flows in protostellar cloud cores. △ Less

Submitted 7 November, 2021; originally announced November 2021.

Comments: 17 pages, 10 figures, accepted for publication in A&A

Journal ref: A&A 658, A53 (2022)

Imaging the water snowline around protostars with water and HCO$^+$ isotopologues

Authors: Merel L. R. van 't Hoff, Daniel Harsono, Martijn L. van Gelder, Tien-Hao Hsieh, John J. Tobin, Sigurd S. Jensen, Naomi Hirano, Jes K. Jørgensen, Edwin A. Bergin, Ewine F. van Dishoeck

Abstract: The water snowline location in protostellar envelopes provides crucial information about the thermal structure and the mass accretion process as it can inform about the occurrence of recent ($\lesssim$1,000 yr) accretion bursts. In addition, the ability to image water emission makes these sources excellent laboratories to test indirect snowline tracers such as H$^{13}$CO$^+$. We study the water sn… ▽ More The water snowline location in protostellar envelopes provides crucial information about the thermal structure and the mass accretion process as it can inform about the occurrence of recent ($\lesssim$1,000 yr) accretion bursts. In addition, the ability to image water emission makes these sources excellent laboratories to test indirect snowline tracers such as H$^{13}$CO$^+$. We study the water snowline in five protostellar envelopes in Perseus using a suite of molecular line observations taken with the Atacama Large Millimeter/submillimeter Array (ALMA) at $\sim$0.2$^{\prime\prime}-$0.7$^{\prime\prime}$ (60--210 au) resolution. B1-c provides a textbook example of compact H$_2^{18}$O ($3_{1,3}-2_{2,0}$) and HDO ($3_{1,2}-2_{2,1}$) emission surrounded by a ring of H$^{13}$CO$^+$ ($J=2-1$) and HC$^{18}$O$^+$ ($J=3-2$). Compact HDO surrounded by H$^{13}$CO$^+$ is also detected toward B1-bS. The optically thick main isotopologue HCO$^+$ is not suited to trace the snowline and HC$^{18}$O$^+$ is a better tracer than H$^{13}$CO$^+$ due to a lower contribution from the outer envelope. However, since a detailed analysis is needed to derive a snowline location from H$^{13}$CO$^+$ or HC$^{18}$O$^+$ emission, their true value as snowline tracer will lie in the application in sources where water cannot be readily detected. For protostellar envelopes, the most straightforward way to locate the water snowline is through observations of H$_2^{18}$O or HDO. Including all sub-arcsecond resolution water observations from the literature, we derive an average burst interval of $\sim$10,000 yr, but high-resolution water observations of a larger number of protostars is required to better constrain the burst frequency. △ Less

Submitted 15 October, 2021; originally announced October 2021.

Comments: Accepted for publication in ApJ

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)

The JCMT Transient Survey: Four Year Summary of Monitoring the Submillimeter Variability of Protostars

Authors: Yong-Hee Lee, Doug Johnstone, Jeong-Eun Lee, Gregory Herczeg, Steve Mairs, Carlos Contreras-Peña, Jennifer Hatchell, Tim Naylor, Graham S. Bell, Tyler L. Bourke, Colton Broughton, Logan Francis, Aashish Gupta, Daniel Harsono, Sheng-Yuan Liu, Geumsook Park, Spencer Plovie, Gerald H. Moriarty-Schieven, Aleks Scholz, Tanvi Sharma, Paula Stella Teixeira, Yao-Te Wang, Yuri Aikawa, Geoffrey C. Bower, Huei-Ru Vivien Chen , et al. (27 additional authors not shown)

Abstract: We present the four-year survey results of monthly submillimeter monitoring of eight nearby ($< 500 $pc) star-forming regions by the JCMT Transient Survey. We apply the Lomb-Scargle Periodogram technique to search for and characterize variability on 295 submillimeter peaks brighter than 0.14 Jy beam$^{-1}$, including 22 disk sources (Class II), 83 protostars (Class 0/I), and 190 starless sources.… ▽ More We present the four-year survey results of monthly submillimeter monitoring of eight nearby ($< 500 $pc) star-forming regions by the JCMT Transient Survey. We apply the Lomb-Scargle Periodogram technique to search for and characterize variability on 295 submillimeter peaks brighter than 0.14 Jy beam$^{-1}$, including 22 disk sources (Class II), 83 protostars (Class 0/I), and 190 starless sources. We uncover 18 secular variables, all of them protostars. No single-epoch burst or drop events and no inherently stochastic sources are observed. We classify the secular variables by their timescales into three groups: Periodic, Curved, and Linear. For the Curved and Periodic cases, the detectable fractional amplitude, with respect to mean peak brightness, is $\sim4$ % for sources brighter than $\sim$ 0.5 Jy beam$^{-1}$. Limiting our sample to only these bright sources, the observed variable fraction is 37 % (16 out of 43). Considering source evolution, we find a similar fraction of bright variables for both Class 0 and Class I. Using an empirically motivated conversion from submillimeter variability to variation in mass accretion rate, six sources (7 % of our full sample) are predicted to have years-long accretion events during which the excess mass accreted reaches more than 40 % above the total quiescently accreted mass: two previously known eruptive Class I sources, V1647 Ori and EC 53 (V371 Ser), and four Class 0 sources, HOPS 356, HOPS 373, HOPS 383, and West 40. Considering the full protostellar ensemble, the importance of episodic accretion on few years timescale is negligible, only a few percent of the assembled mass. However, given that this accretion is dominated by events of order the observing time-window, it remains uncertain as to whether the importance of episodic events will continue to rise with decades-long monitoring. △ Less

Submitted 22 July, 2021; originally announced July 2021.

Comments: Accepted for publication in the Astrophysical Journal

Which molecule traces what: chemical diagnostics of protostellar sources

Authors: Łukasz Tychoniec, Ewine F. van Dishoeck, Merel L. R. van 't Hoff, Martijn L. van Gelder, Benoît Tabone, Yuan Chen, Daniel Harsono, Charles L. H. Hull, Michiel R. Hogerheijde, Nadia M. Murillo, John J. Tobin

Abstract: The physical and chemical conditions in Class 0/I protostars are fundamental in unlocking the protostellar accretion process and its impact on planet formation. The aim is to determine which physical components are traced by different molecules at sub-arcsecond scales (100 - 400 au). We use a suite of Atacama Large Millimeter/submillimeter Array (ALMA) datasets in Band 6 (1 mm), Band 5 (1.8 mm) an… ▽ More The physical and chemical conditions in Class 0/I protostars are fundamental in unlocking the protostellar accretion process and its impact on planet formation. The aim is to determine which physical components are traced by different molecules at sub-arcsecond scales (100 - 400 au). We use a suite of Atacama Large Millimeter/submillimeter Array (ALMA) datasets in Band 6 (1 mm), Band 5 (1.8 mm) and Band 3 (3 mm) at spatial resolutions 0.5 - 3" for 16 protostellar sources. The protostellar envelope is well traced by C$^{18}$O, DCO$^+$ and N$_2$D$^+$, with the freeze-out of CO governing the chemistry at envelope scales. Molecular outflows are seen in classical shock tracers like SiO and SO, but ice-mantle products such as CH$_3$OH and HNCO released with the shock are also observed. The molecular jet is prominent not only in SiO and SO but also occasionally in H$_2$CO. The cavity walls show tracers of UV-irradiation such as C$_2$H c-C$_3$H$_2$ and CN. The hot inner envelope, apart from showing emission from complex organic molecules (COMs), also presents compact emission from small molecules like H$_2$S, SO, OCS and H$^{13}$CN, most likely related to ice sublimation and high-temperature chemistry. Sub-arcsecond millimeter-wave observations allow to identify those (simple) molecules that best trace each of the physical components of a protostellar system. COMs are found both in the hot inner envelope (high excitation lines) and in the outflows (lower-excitation lines) with comparable abundances. COMs can coexist with hydrocarbons in the same protostellar sources, but they trace different components. In the near future, mid-IR observations with JWST-MIRI will provide complementary information about the hottest gas and the ice mantle content, at unprecedented sensitivity and at resolutions comparable to ALMA for the same sources. △ Less

Submitted 13 July, 2021; v1 submitted 8 July, 2021; originally announced July 2021.

Comments: 38 pages, 33 figures, accepted for publication to A\&A

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

Resolved molecular line observations reveal an inherited molecular layer in the young disk around TMC1A

Authors: Daniel Harsono, Matthijs van der Wiel, Per Bjerkeli, Jon Ramsey, Hannah Calcutt, Lars Kristensen, Jes Jørgensen

Abstract: Physical processes that govern the star and planet formation sequence influence the chemical composition and evolution of protoplanetary disks. To understand the chemical composition of protoplanets, we need to constrain the composition and structure of the disks from whence they are formed. We aim to determine the molecular abundance structure of the young disk around the TMC1A protostar on au sc… ▽ More Physical processes that govern the star and planet formation sequence influence the chemical composition and evolution of protoplanetary disks. To understand the chemical composition of protoplanets, we need to constrain the composition and structure of the disks from whence they are formed. We aim to determine the molecular abundance structure of the young disk around the TMC1A protostar on au scales in order to understand its chemical structure and any possible implications for disk formation. We present spatially resolved Atacama Large Millimeter/submillimeter Array observations of CO, $HCO^{+}$, HCN, DCN, and SO line emission, as well as dust continuum emission, in the vicinity of TMC1A. Molecular column densities are estimated both under the assumption of optically thin emission from molecules in LTE as well as through more detailed non-LTE radiative transfer calculations. Resolved dust continuum emission from the disk is detected between 220 and 260 GHz. Rotational transitions from HCO$^{+}$, HCN, and SO are also detected from the inner 100 au region. From the derived $HCO^{+}$ abundance, we estimate the ionization fraction of the disk surface and find values that imply that the accretion process is not driven by the magneto-rotational instability. The molecular abundances averaged over the TMC1A disk are similar to its protostellar envelope and other, older Class II disks. We meanwhile find a discrepancy between the young disk's molecular abundances relative to Solar System objects. Abundance comparisons between the disk and its surrounding envelope for several molecular species reveal that the bulk of planet-forming material enters the disk unaltered. Differences in HCN and $H_2 O$ molecular abundances between the disk around TMC1A, Class II disks, and Solar System objects trace the chemical evolution during disk and planet formation. △ Less

Submitted 26 October, 2020; originally announced October 2020.

Comments: Accepted for publication on A&A; 15 pages, 13 figures

Journal ref: A&A 646, A72 (2021)

On the accuracy of the ALMA flux calibration in the time domain and across spectral windows

Authors: Logan Francis, Doug Johnstone, Gregory Herczeg, Todd R. Hunter, Daniel Harsono

Abstract: A diverse array of science goals require accurate flux calibration of observations with the Atacama Large Millimeter/Submillimeter array (ALMA), however, this goal remains challenging due to the stochastic time-variability of the ``grid'' quasars ALMA uses for calibration. In this work, we use 343.5 GHz (Band 7) ALMA Atacama Compact Array observations of four bright and stable young stellar object… ▽ More A diverse array of science goals require accurate flux calibration of observations with the Atacama Large Millimeter/Submillimeter array (ALMA), however, this goal remains challenging due to the stochastic time-variability of the ``grid'' quasars ALMA uses for calibration. In this work, we use 343.5 GHz (Band 7) ALMA Atacama Compact Array observations of four bright and stable young stellar objects over 7 epochs to independently assess the accuracy of the ALMA flux calibration and to refine the relative calibration across epochs. The use of these four extra calibrators allow us to achieve an unprecedented relative ALMA calibration accuracy of $\sim 3\%$. On the other hand, when the observatory calibrator catalog is not up-to-date, the Band 7 data calibrated by the ALMA pipeline may have a flux calibration poorer than the nominal 10\%, which can be exacerbated by weather-related phase decorrelation when self-calibration of the science target is either not possible or not attempted. We also uncover a relative flux calibration uncertainty between spectral windows of 0.8\%, implying that measuring spectral indices within a single ALMA band is likely highly uncertain. We thus recommend various methods for science goals requiring high flux accuracy and robust calibration, in particular, the observation of additional calibrators combined with a relative calibration strategy, and observation of solar system objects for high absolute accuracy. △ Less

Submitted 5 October, 2020; originally announced October 2020.

APEX-SEPIA660 Early Science: Gas at densities above $10^7$ cm$^{-3}$ towards OMC-1

Authors: A. Hacar, M. R. Hogerheijde, D. Harsono, S. Portegies Zwart, C. De Breuck, K. Torstensson, W. Boland, A. M. Baryshev, R. Hesper, J. Barkhof, J. Adema, M. E. Bekema, A. Koops, A. Khudchenko, R. Stark

Abstract: Context. The star formation rates and stellar densities found in young massive clusters suggest that these stellar systems originate from gas at densities n(H$_2$) $>10^7$ cm$^{-3}$. Until today, however, the physical characterization of this ultra high density material remains largely unconstrained in observations. Aims. We investigated the density properties of the star-forming gas in the OMC-1… ▽ More Context. The star formation rates and stellar densities found in young massive clusters suggest that these stellar systems originate from gas at densities n(H$_2$) $>10^7$ cm$^{-3}$. Until today, however, the physical characterization of this ultra high density material remains largely unconstrained in observations. Aims. We investigated the density properties of the star-forming gas in the OMC-1 region located in the vicinity of the Orion Nebula Cluster (ONC). Methods. We mapped the molecular emission at 652 GHz in OMC-1 as part of the APEX-SEPIA660 Early Science. Results. We detect bright and extended N$_2$H$^+$ (J=7-6) line emission along the entire OMC-1 region. Comparisons with previous ALMA data of the (J=1-0) transition and radiative transfer models indicate that the line intensities observed in this N$_2$H$^+$ (7-6) line are produced by large mass reservoirs of gas at densities n(H$_2$) $>10^7$ cm$^{-3}$. Conclusions. The first detection of this N$_2$H$^+$ (7-6) line at parsec-scales demonstrates the extreme density conditions of the star-forming gas in young massive clusters such as the ONC. Our results highlight the unique combination of sensitivity and mapping capabilities of the new SEPIA660 receiver for the study of the ISM properties at high frequencies. △ Less

Submitted 26 September, 2020; originally announced September 2020.

Comments: 8 pages, 5 Figures. Accepted for publication by A&A

Journal ref: A&A 644, A133 (2020)

Temperature structures of embedded disks: young disks in Taurus are warm

Authors: Merel L. R. van 't Hoff, Daniel Harsono, John J. Tobin, Arthur D. Bosman, Ewine F. van Dishoeck, Jes K. Jørgensen, Anna Miotello, Nadia M. Murillo, Catherine Walsh

Abstract: The chemical composition of gas and ice in disks around young stars set the bulk composition of planets. In contrast to protoplanetary disks (Class II), young disks that are still embedded in their natal envelope (Class 0 and I) are predicted to be too warm for CO to freeze out, as has been confirmed observationally for L1527 IRS. To establish whether young disks are generally warmer than their mo… ▽ More The chemical composition of gas and ice in disks around young stars set the bulk composition of planets. In contrast to protoplanetary disks (Class II), young disks that are still embedded in their natal envelope (Class 0 and I) are predicted to be too warm for CO to freeze out, as has been confirmed observationally for L1527 IRS. To establish whether young disks are generally warmer than their more evolved counterparts, we observed five young (Class 0/I and Class I) disks in Taurus with the Atacama Large Millimeter/submillimeter Array (ALMA), targeting C$^{17}$O $2-1$, H$_2$CO $3_{1,2}-2_{1,1}$, HDO $3_{1,2}-2_{2,1}$ and CH$_3$OH $5_K-4_K$ transitions at $0.48^{\prime\prime} \times 0.31^{\prime\prime}$ resolution. The different freeze-out temperatures of these species allow us to derive a global temperature structure. C$^{17}$O and H$_2$CO are detected in all disks, with no signs of CO freeze-out in the inner $\sim$100 au, and a CO abundance close to $\sim$10$^{-4}$. H$_2$CO emission originates in the surface layers of the two edge-on disks, as witnessed by the especially beautiful V-shaped emission pattern in IRAS~04302+2247. HDO and CH$_3$OH are not detected, with column density upper limits more than 100 times lower than for hot cores. Young disks are thus found to be warmer than more evolved protoplanetary disks around solar analogues, with no CO freeze-out (or only in the outermost part of $\gtrsim$100 au disks) or CO processing. However, they are not as warm as hot cores or disks around outbursting sources, and therefore do not have a large gas-phase reservoir of complex molecules. △ Less

Submitted 18 August, 2020; originally announced August 2020.

Comments: Accepted for publication in ApJ. 19 pages, 11 figures, 3 tables (+ appendix)

Dual-Wavelength ALMA Observations of Dust Rings in Protoplanetary Disks

Authors: Feng Long, Paola Pinilla, Gregory J. Herczeg, Sean M. Andrews, Daniel Harsono, Doug Johnstone, Enrico Ragusa, Ilaria Pascucci, David J. Wilner, Nathan Hendler, Jeff Jennings, Yao Liu, Giuseppe Lodato, Francois Menard, Gerrit van de Plas, Giovanni Dipierro

Abstract: We present new Atacama Large Millimeter/submillimeter Array (ALMA) observations for three protoplanetary disks in Taurus at 2.9\,mm and comparisons with previous 1.3\,mm data both at an angular resolution of $\sim0.''1$ (15\,au for the distance of Taurus). In the single-ring disk DS Tau, double-ring disk GO Tau, and multiple-ring disk DL Tau, the same rings are detected at both wavelengths, with r… ▽ More We present new Atacama Large Millimeter/submillimeter Array (ALMA) observations for three protoplanetary disks in Taurus at 2.9\,mm and comparisons with previous 1.3\,mm data both at an angular resolution of $\sim0.''1$ (15\,au for the distance of Taurus). In the single-ring disk DS Tau, double-ring disk GO Tau, and multiple-ring disk DL Tau, the same rings are detected at both wavelengths, with radial locations spanning from 50 to 120\,au. To quantify the dust emission morphology, the observed visibilities are modeled with a parametric prescription for the radial intensity profile. The disk outer radii, taken as 95\% of the total flux encircled in the model intensity profiles, are consistent at both wavelengths for the three disks. Dust evolution models show that dust trapping in local pressure maxima in the outer disk could explain the observed patterns. Dust rings are mostly unresolved. The marginally resolved ring in DS Tau shows a tentatively narrower ring at the longer wavelength, an observational feature expected from efficient dust trapping. The spectral index ($α_{\rm mm}$) increases outward and exhibits local minima that correspond to the peaks of dust rings, indicative of the changes in grain properties across the disks. The low optical depths ($τ\sim$0.1--0.2 at 2.9\,mm and 0.2--0.4 at 1.3\,mm) in the dust rings suggest that grains in the rings may have grown to millimeter sizes. The ubiquitous dust rings in protoplanetary disks modify the overall dynamics and evolution of dust grains, likely paving the way towards the new generation of planet formation. △ Less

Submitted 4 June, 2020; originally announced June 2020.

Comments: accepted for publication in ApJ

Spirals inside the millimeter cavity of transition disk SR 21

Authors: G. A. Muro-Arena, C. Ginski, C. Dominik, M. Benisty, P. Pinilla, A. J. Bohn, T. Moldenhauer, W. Kley, D. Harsono, T. Henning, R. G. van Holstein, M. Janson, M. Keppler, F. Ménard, L. M. Pérez, T. Stolker, M. Tazzari, M. Villenave, A. Zurlo, C. Petit, F. Rigal, O. Möller-Nilsson, M. Llored, T. Moulin, P. Rabou

Abstract: Hydrodynamical simulations of planet-disk interactions suggest that planets may be responsible for a number of the sub-structures frequently observed in disks in both scattered light and dust thermal emission. Despite the ubiquity of these features, direct evidence of planets embedded in disks and of the specific interaction features like spiral arms within planetary gaps still remain rare. In thi… ▽ More Hydrodynamical simulations of planet-disk interactions suggest that planets may be responsible for a number of the sub-structures frequently observed in disks in both scattered light and dust thermal emission. Despite the ubiquity of these features, direct evidence of planets embedded in disks and of the specific interaction features like spiral arms within planetary gaps still remain rare. In this study we discuss recent observational results in the context of hydrodynamical simulations in order to infer the properties of a putative embedded planet in the cavity of a transition disk. We imaged the transition disk SR 21 in H-band in scattered light with SPHERE/IRDIS and in thermal dust emission with ALMA band 3 (3mm) observations at a spatial resolution of 0.1". We combine these datasets with existing band 9 (430um) and band 7 (870um) ALMA continuum data. The Band 3 continuum data reveals a large cavity and a bright ring peaking at 53 au strongly suggestive of dust trapping.The ring shows a pronounced azimuthal asymmetry, with a bright region in the north-west that we interpret as a dust over-density. A similarly-asymmetric ring is revealed at the same location in polarized scattered light, in addition to a set of bright spirals inside the mm cavity and a fainter spiral bridging the gap to the outer ring. These features are consistent with a number of previous hydrodynamical models of planet-disk interactions, and suggest the presence of a ~1 MJup planet at 44 au and PA=11°. This makes SR21 the first disk showing spiral arms inside the mm cavity, as well as one for which the location of a putative planet can be precisely inferred. With the location of a possible planet being well-constrained by observations, it is an ideal candidate for follow-up observations to search for direct evidence of a planetary companion still embedded in its disk. △ Less

Submitted 18 March, 2020; originally announced March 2020.

Comments: 7 pages, 4 figures

Missing water in Class I protostellar disks

Authors: Daniel Harsono, Magnus Persson, Alyssa Ramos, Nadia Murillo, Luke Maud, Michiel Hogerheijde, Arthur Bosman, Lars Kristensen, Jes Jorgensen, Ted Bergin, Ruud Visser, Joe Mottram, Ewine van Dishoeck

Abstract: Water is a key volatile that provides insights into the initial stages of planet formation. The low water abundances inferred from water observations toward low-mass protostellar objects may point to a rapid locking of water as ice by large dust grains during star and planet formation. However, little is known about the water vapor abundance in newly formed planet-forming disks. We aim to determin… ▽ More Water is a key volatile that provides insights into the initial stages of planet formation. The low water abundances inferred from water observations toward low-mass protostellar objects may point to a rapid locking of water as ice by large dust grains during star and planet formation. However, little is known about the water vapor abundance in newly formed planet-forming disks. We aim to determine the water abundance in embedded Keplerian disks through spatially-resolved observations of H$_2^{18}$O lines to understand the evolution of water during star and planet formation. We present H$_2^{18}$O line observations with ALMA and NOEMA millimeter interferometers toward five young stellar objects. NOEMA observed the 3$_{1,3}$ - $2_{2,0}$ line (E$_{\rm up}$ = 203.7 K) while ALMA targeted the $4_{1,4}$ - $3_{2,1}$ line (E$_{\rm up}$ = 322.0 K). Water column densities are derived considering optically thin and thermalized emission. Our observations are sensitive to the emission from the known Keplerian disks around three out of the five Class I objects in the sample. No H$_2^{18}$O emission is detected toward any of our five Class I disks. We report upper limits to the integrated line intensities. The inferred water column densities in Class I disks are N < 10$^{15}$ cm$^{-2}$ on 100 au scales which include both disk and envelope. The upper limits imply a disk-averaged water abundance of $\lesssim 10^{-6}$ with respect to H$_2$ for Class I objects. After taking into account the physical structure of the disk, the upper limit to the water abundance averaged over the inner warm disk with $T>$ 100 K is between 10$^{-7}$ up to 10$^{-5}$. Water vapor is not abundant in warm protostellar envelopes around Class I protostars. Upper limits to the water vapor column densities in Class I disks are at least two orders magnitude lower than values found in Class 0 disk-like structures. △ Less

Submitted 26 February, 2020; originally announced February 2020.

Comments: 13 pages + 6 pages appendix, 7 figures, accepted for publication in aanda

Journal ref: A&A 636, A26 (2020)

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)

Observational constraints on dust disk sizes in tidally truncated protoplanetary disks in multiple systems in the Taurus region

Authors: C. F. Manara, M. Tazzari, F. Long, G. J. Herczeg, G. Lodato, A. A. Rota, P. Cazzoletti, G. van der Plas, P. Pinilla, G. Dipierro, S. Edwards, D. Harsono, D. Johnstone, Y. Liu, F. Menard, B. Nisini, E. Ragusa, Y. Boehler, S. Cabrit

Abstract: The impact of stellar multiplicity on the evolution of planet-forming disks is still the subject of debate. Here we present and analyze disk structures around ten multiple stellar systems that were included in an unbiased, high spatial resolution survey performed with ALMA of 32 protoplanetary disks in the Taurus star-forming region. At the unprecedented spatial resolution of ~0.12" we detect and… ▽ More The impact of stellar multiplicity on the evolution of planet-forming disks is still the subject of debate. Here we present and analyze disk structures around ten multiple stellar systems that were included in an unbiased, high spatial resolution survey performed with ALMA of 32 protoplanetary disks in the Taurus star-forming region. At the unprecedented spatial resolution of ~0.12" we detect and spatially resolve the disks around all primary stars, and those around eight secondary and one tertiary star. The dust radii of disks around multiple stellar systems are smaller than those around single stars in the same stellar mass range and in the same region. The disks in multiple stellar systems also show a steeper decay of the millimeter continuum emission at the outer radius than disks around single stars, suggestive of the impact of tidal truncation on the shape of the disks in multiple systems. However, the observed ratio between the dust disk radii and the observed separation of the stars in the multiple systems is consistent with analytic predictions of the effect of tidal truncation only if the eccentricities of the binaries are rather high (typically >0.5), or if the observed dust radii are a factor of two smaller than the gas radii, as is typical for isolated systems. Similar high-resolution studies targeting the gaseous emission from disks in multiple stellar systems are required to resolve this question. △ Less

Submitted 18 July, 2019; v1 submitted 8 July, 2019; originally announced July 2019.

Comments: Accepted on Astronomy & Astrophysics; 15 pages plus appendix, revised version after language editing

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

Compact Disks in a High Resolution ALMA Survey of Dust Structures in the Taurus Molecular Cloud

Authors: Feng Long, Gregory J. Herczeg, Daniel Harsono, Paola Pinilla, Marco Tazzari, Carlo F. Manara, Ilaria Pascucci, Sylvie Cabrit, Brunella Nisini, Doug Johnstone, Suzan Edwards, Colette Salyk, Francois Menard, Giuseppe Lodato, Yann Boehler, Gregory N. Mace, Yao Liu, Gijs D. Mulders, Nathanial Hendler, Enrico Ragusa, William J. Fischer, Andrea Banzatti, Elisabetta Rigliaco, Gerrit van der Plas, Giovanni Dipierro , et al. (2 additional authors not shown)

Abstract: We present a high-resolution ($\sim0.''12$, $\sim16$ au, mean sensitivity of $50~μ$Jy~beam$^{-1}$ at 225 GHz) snapshot survey of 32 protoplanetary disks around young stars with spectral type earlier than M3 in the Taurus star-forming region using Atacama Large Millimeter Array (ALMA). This sample includes most mid-infrared excess members that were not previously imaged at high spatial resolution,… ▽ More We present a high-resolution ($\sim0.''12$, $\sim16$ au, mean sensitivity of $50~μ$Jy~beam$^{-1}$ at 225 GHz) snapshot survey of 32 protoplanetary disks around young stars with spectral type earlier than M3 in the Taurus star-forming region using Atacama Large Millimeter Array (ALMA). This sample includes most mid-infrared excess members that were not previously imaged at high spatial resolution, excluding close binaries and highly extincted objects, thereby providing a more representative look at disk properties at 1--2 Myr. Our 1.3 mm continuum maps reveal 12 disks with prominent dust gaps and rings, 2 of which are around primary stars in wide binaries, and 20 disks with no resolved features at the observed resolution (hereafter smooth disks), 8 of which are around the primary star in wide binaries. The smooth disks were classified based on their lack of resolved substructures, but their most prominent property is that they are all compact with small effective emission radii ($R_{\rm eff,95\%} \lesssim 50$ au). In contrast, all disks with $R_{\rm eff,95\%}$ of at least 55 au in our sample show detectable substructures. Nevertheless, their inner emission cores (inside the resolved gaps) have similar peak brightness, power law profiles, and transition radii to the compact smooth disks, so the primary difference between these two categories is the lack of outer substructures in the latter. These compact disks may lose their outer disk through fast radial drift without dust trapping, or they might be born with small sizes. The compact dust disks, as well as the inner disk cores of extended ring disks, that look smooth at the current resolution will likely show small-scale or low-contrast substructures at higher resolution. The correlation between disk size and disk luminosity correlation demonstrates that some of the compact disks are optically thick at millimeter wavelengths. △ Less

Submitted 25 June, 2019; originally announced June 2019.

Comments: ApJ accepted, 25 pages, 11 figures

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

Authors: L. T. Maud, R. Cesaroni, M. S. N. Kumar, V. M. Rivilla, A. Ginsburg, P. D. Klaassen, D. Harsono, A. Sanchez-Monge, A. Ahmadi, V. Allen, M. T. Beltran, H. Beuther, R. Galvan-Madrid, C. Goddi, M. G. Hoare, M. R. Hogerheijde, K. G. Johnston, R. Kuiper, L. Moscadelli, T. Peters, L. Testi, F. F. S. van der Tak, W. J. de Wit

Abstract: We present the highest angular resolution (20x15mas - 44x33au) Atacama Large Millimeter/sub-millimeter Array (ALMA) observations currently possible of the proto-O-star G17.64+0.16 in Band 6. The Cycle 5 observations with baselines out to 16km probes scales <50au and reveal the rotating disc around G17.64+0.16, a massive forming O-type star. The disc has a ring-like enhancement in the dust emission… ▽ More We present the highest angular resolution (20x15mas - 44x33au) Atacama Large Millimeter/sub-millimeter Array (ALMA) observations currently possible of the proto-O-star G17.64+0.16 in Band 6. The Cycle 5 observations with baselines out to 16km probes scales <50au and reveal the rotating disc around G17.64+0.16, a massive forming O-type star. The disc has a ring-like enhancement in the dust emission, especially visible as arc structures to the north and south. The Keplerian kinematics are most prominently seen in the vibrationally excited water line, H2O (Eu=3461.9K). The mass of the central source found by modelling the Keplerian rotation is consistent with 45+/-10Mo. The H30alpha (231.9GHz) radio-recombination line and the SiO (5-4) molecular line were detected at up to the 10 sigma$ level. The estimated disc mass is 0.6-2.6Mo under the optically thin assumption. Analysis of the Toomre Q parameter, in the optically thin regime, indicates that the disc stability is highly dependent on temperature. The disc currently appears stable for temperatures >150K, this does not preclude that the substructures formed earlier through disc fragmentation. △ Less

Submitted 15 June, 2019; originally announced June 2019.

Comments: 9 Total, inc references, appendix, long author list, 4 figures (inc appendix)

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

An Inner Disk in the Large Gap of the Transition Disk SR 24S

Authors: Paola Pinilla, Myriam Benisty, Paolo Cazzoletti, Daniel Harsono, Laura M. Pérez, Marco Tazzari

Abstract: We report new Atacama Large Millimeter/sub-millimeter Array (ALMA) Band 3 observations at 2.75 mm of the TD around SR 24S with an angular resolution of $\sim$0.11''$\times$ 0.09'' and a peak signal-to-noise ratio of $\sim24$. We detect an inner disk and a mostly symmetric ring-like structure that peaks at $\sim$0.32'', that is $\sim$37 au at a distance of $\sim$114.4 pc. The full width at half max… ▽ More We report new Atacama Large Millimeter/sub-millimeter Array (ALMA) Band 3 observations at 2.75 mm of the TD around SR 24S with an angular resolution of $\sim$0.11''$\times$ 0.09'' and a peak signal-to-noise ratio of $\sim24$. We detect an inner disk and a mostly symmetric ring-like structure that peaks at $\sim$0.32'', that is $\sim$37 au at a distance of $\sim$114.4 pc. The full width at half maximum of this ring is $\sim$28 au. We analyze the observed structures by fitting the dust continuum visibilities using different models for the intensity profile, and compare with previous ALMA observations of the same disk at 0.45 mm and 1.30 mm. We qualitatively compare the results of these fits with theoretical predictions of different scenarios for the formation of a cavity or large gap. The comparison of the dust continuum structure between different ALMA bands indicates that photoevaporation and dead zone can be excluded as leading mechanisms for the cavity formation in SR 24S disk, leaving the planet scenario (single or multiple planets) as the most plausible mechanism. We compared the 2.75 mm emission with published (sub-)centimeter data and find that the inner disk is likely tracing dust thermal emission. This implies that any companion in the system should allow dust to move inwards throughout the gap and replenish the inner disk. In the case of one single planet, this puts strong constraints on the mass of the potential planet inside the cavity and the disk viscosity of about $\lesssim$5 $M_{\rm{Jup}}$ and $α\sim10^{-4}-10^{-3}$, respectively. △ Less

Submitted 25 April, 2019; originally announced April 2019.

Comments: Accepted to ApJ

The newborn planet population emerging from ring-like structures in discs

Authors: G. Lodato, G. Dipierro, E. Ragusa, F. Long, G. J. Herczeg, I. Pascucci, P. Pinilla, C. F. Manara, M. Tazzari, Y. Liu, G. D. Mulders, D. Harsono, Y. Boehler, F. Menard, D. Johnstone, C. Salyk, G. van der Plas, S. Cabrit, S. Edwards, W. J. Fischer, N. Hendler, B. Nisini, E. Rigliaco, H. Avenhaus, A. Banzatti , et al. (1 additional authors not shown)

Abstract: ALMA has observed a plethora of ring-like structures in planet forming discs at distances of 10-100 au from their host star. Although several mechanisms have been invoked to explain the origin of such rings, a common explanation is that they trace new-born planets. Under the planetary hypothesis, a natural question is how to reconcile the apparently high frequency of gap-carving planets at 10-100… ▽ More ALMA has observed a plethora of ring-like structures in planet forming discs at distances of 10-100 au from their host star. Although several mechanisms have been invoked to explain the origin of such rings, a common explanation is that they trace new-born planets. Under the planetary hypothesis, a natural question is how to reconcile the apparently high frequency of gap-carving planets at 10-100 au with the paucity of Jupiter mass planets observed around main sequence stars at those separations. Here, we provide an analysis of the new-born planet population emerging from observations of gaps in discs, under the assumption that the observed gaps are due to planets. We use a simple estimate of the planet mass based on the gap morphology, and apply it to a sample of gaps recently obtained by us in a survey of Taurus with ALMA. We also include additional data from recent published surveys, thus analysing the largest gap sample to date, for a total of 48 gaps. The properties of the purported planets occupy a distinctively different region of parameter space with respect to the known exo-planet population, currently not accessible through planet finding methods. Thus, no discrepancy in the mass and radius distribution of the two populations can be claimed at this stage. We show that the mass of the inferred planets conforms to the theoretically expected trend for the minimum planet mass needed to carve a dust gap. Finally, we estimate the separation and mass of the putative planets after accounting for migration and accretion, for a range of evolutionary times, finding a good match with the distribution of cold Jupiters. △ Less

Submitted 12 March, 2019; originally announced March 2019.

Comments: 9 pages, 5 figures, re-submitted to MNRAS after referee's comments

The Ring Structure in the MWC 480 Disk Revealed by ALMA

Authors: Yao Liu, Giovanni Dipierro, Enrico Ragusa, Giuseppe Lodato, Gregory J. Herczeg, Feng Long, Daniel Harsono, Yann Boehler, Francois Menard, Doug Johnstone, Ilaria Pascucci, Paola Pinilla, Colette Salyk, Gerrit van der Plas, Sylvie Cabrit, William J. Fischer, Nathan Hendler, Carlo F. Manara, Brunella Nisini, Elisabetta Rigliaco, Henning Avenhaus, Andrea Banzatti, Michael Gully-Santiago

Abstract: Gap-like structures in protoplanetary disks are likely related to planet formation processes. In this paper, we present and analyze high resolution (0.17*0.11 arcsec) 1.3 mm ALMA continuum observations of the protoplanetary disk around the Herbig Ae star MWC 480. Our observations for the first time show a gap centered at ~74au with a width of ~23au, surrounded by a bright ring centered at ~98au fr… ▽ More Gap-like structures in protoplanetary disks are likely related to planet formation processes. In this paper, we present and analyze high resolution (0.17*0.11 arcsec) 1.3 mm ALMA continuum observations of the protoplanetary disk around the Herbig Ae star MWC 480. Our observations for the first time show a gap centered at ~74au with a width of ~23au, surrounded by a bright ring centered at ~98au from the central star. Detailed radiative transfer modeling of both the ALMA image and the broadband spectral energy distribution is used to constrain the surface density profile and structural parameters of the disk. If the width of the gap corresponds to 4~8 times the Hill radius of a single forming planet, then the putative planet would have a mass of 0.4~3 M_Jup. We test this prediction by performing global three-dimensional smoothed particle hydrodynamic gas/dust simulations of disks hosting a migrating and accreting planet. We find that the dust emission across the disk is consistent with the presence of an embedded planet with a mass of ~2.3 M_Jup at an orbital radius of ~78au. Given the surface density of the best-fit radiative transfer model, the amount of depleted mass in the gap is higher than the mass of the putative planet, which satisfies the basic condition for the formation of such a planet. △ Less

Submitted 7 November, 2018; originally announced November 2018.

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

Journal ref: A&A 622, A75 (2019)

Gaps and Rings in an ALMA Survey of Disks in the Taurus Star-forming Region

Authors: Feng Long, Paola Pinilla, Gregory J. Herczeg, Daniel Harsono, Giovanni Dipierro, Ilaria Pascucci, Nathan Hendler, Marco Tazzari, Enrico Ragusa, Colette Salyk, Suzan Edwards, Giuseppe Lodato, Gerrit van de Plas, Doug Johnstone, Yao Liu, Yann Boehler, Sylvie Cabrit, Carlo F. Manara, Francois Menard, Gijs D. Mulders, Brunella Nisini, William J. Fischer, Elisabetta Rigliaco, Andrea Banzatti, Henning Avenhaus , et al. (1 additional authors not shown)

Abstract: Rings are the most frequently revealed substructure in ALMA dust observations of protoplanetary disks, but their origin is still hotly debated. In this paper, we identify dust substructures in 12 disks and measure their properties to investigate how they form. This subsample of disks is selected from a high-resolution ($\sim0.12''$) ALMA 1.33 mm survey of 32 disks in the Taurus star-forming region… ▽ More Rings are the most frequently revealed substructure in ALMA dust observations of protoplanetary disks, but their origin is still hotly debated. In this paper, we identify dust substructures in 12 disks and measure their properties to investigate how they form. This subsample of disks is selected from a high-resolution ($\sim0.12''$) ALMA 1.33 mm survey of 32 disks in the Taurus star-forming region, which was designed to cover a wide range of sub-mm brightness and to be unbiased to previously known substructures. While axisymmetric rings and gaps are common within our sample, spiral patterns and high contrast azimuthal asymmetries are not detected. Fits of disk models to the visibilities lead to estimates of the location and shape of gaps and rings, the flux in each disk component, and the size of the disk. The dust substructures occur across a wide range of stellar mass and disk brightness. Disks with multiple rings tend to be more massive and more extended. The correlation between gap locations and widths, the intensity contrast between rings and gaps, and the separations of rings and gaps could all be explained if most gaps are opened by low-mass planets (super-Earths and Neptunes) in the condition of low disk turbulence ($α=10^{-4}$). The gap locations are not well correlated with the expected locations of CO and N$_2$ ice lines, so condensation fronts are unlikely to be a universal mechanism to create gaps and rings, though they may play a role in some cases. △ Less

Submitted 14 October, 2018; originally announced October 2018.

Comments: 23 pages, 11 figures, ApJ accepted

Chasing discs around O-type (proto)stars - ALMA evidence for an SiO disc and disc wind from G17.64+0.16

Authors: L. T. Maud, R. Cesaroni, M. S. N. Kumar, F. F. S. van der Tak, V. Allen, M. G. Hoare, P. D. Klaassen, D. Harsono, M. R. Hogerheijde, Á. Sánchez-Monge, P. Schilke, A. Ahmadi, M. T. Beltrán, H. Beuther, T. Csengeri, S. Etoka, G. Fuller, R. Galván-Madrid, C. Goddi, Th. Henning, K. G. Johnston, R. Kuiper, S. Lumsden, L. Moscadelli, J. C. Mottram , et al. (6 additional authors not shown)

Abstract: We present high angular resolution 0.2 arcsec continuum and molecular emission line Atacama Large Millimeter/sub-millimeter Array (ALMA) observations of G17.64+0.16 in Band 6 (220GHz) taken as part of a campaign in search of circumstellar discs around (proto)-O-stars. At a resolution of 400au the main continuum core is essentially unresolved and isolated from other strong and compact emission peak… ▽ More We present high angular resolution 0.2 arcsec continuum and molecular emission line Atacama Large Millimeter/sub-millimeter Array (ALMA) observations of G17.64+0.16 in Band 6 (220GHz) taken as part of a campaign in search of circumstellar discs around (proto)-O-stars. At a resolution of 400au the main continuum core is essentially unresolved and isolated from other strong and compact emission peaks. At a resolution of 400au the main continuum core is essentially unresolved and isolated from other strong and compact emission peaks. We detect SiO (5-4) emission that is marginally resolved and elongated in a direction perpendicular to the large-scale outflow seen in the 13CO (2-1) line using the main ALMA array in conjunction with the Atacama Compact Array (ACA). Morphologically, the SiO appears to represent a disc-like structure. Using parametric models we show that the position-velocity profile of the SiO is consistent with the Keplerian rotation of a disc around an object between 10-30Mo in mass, only if there is also radial expansion from a separate structure. The radial motion component can be interpreted as a disc wind from the disc surface. Models with a central stellar object mass between 20 and 30Mo are the most consistent with the stellar luminosity (100000 Lo) and indicative of an O-type star. The H30a millimetre recombination line (231.9GHz) is also detected, but spatially unresolved, and is indicative of a very compact, hot, ionised region co-spatial with the dust continuum core. Accounting for all observables, we suggest that G17.64 is consistent with a O-type young stellar object in the final stages of protostellar assembly, driving a wind, but that has not yet developed into a compact HII region. The existance and detection of the disc in G17.64 is likely related to its isolated and possibly more evolved nature, traits which may underpin discs in similar sources. △ Less

Submitted 9 October, 2018; originally announced October 2018.

Comments: 23 pages, 13 fig (inc 3 Appendix figures)

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

Methanol and its relation to the water snowline in the disk around the young outbursting star V883 Ori

Authors: Merel L. R. van 't Hoff, John J. Tobin, Leon Trapman, Daniel Harsono, Patrick D. Sheehan, William J. Fischer, S. Thomas Megeath, Ewine F. van Dishoeck

Abstract: We report the detection of methanol in the disk around the young outbursting star V883 Ori with the Atacama Large Millimeter/submillimeter Array (ALMA). Four transitions are observed with upper level energies ranging between 115 and 459 K. The emission is spatially resolved with the 0.14" beam and follows the Keplerian rotation previously observed for C$^{18}$O. Using a rotational diagram analysis… ▽ More We report the detection of methanol in the disk around the young outbursting star V883 Ori with the Atacama Large Millimeter/submillimeter Array (ALMA). Four transitions are observed with upper level energies ranging between 115 and 459 K. The emission is spatially resolved with the 0.14" beam and follows the Keplerian rotation previously observed for C$^{18}$O. Using a rotational diagram analysis, we find a disk-averaged column density of $\sim10^{17}$ cm$^{-2}$ and a rotational temperature of $\sim90-100$ K, suggesting that the methanol has thermally desorbed from the dust grains. We derive outer radii between 120 and 140 AU for the different transitions, compared to the 360 AU outer radius for C$^{18}$O. Depending on the exact physical structure of the disk, the methanol emission could originate in the surface layers beyond the water snowline. Alternatively, the bulk of the methanol emission originates inside the water snowline, which can then be as far out as ~100 AU, instead of 42 AU as was previously inferred from the continuum opacity. In addition, these results show that outbursting young stars like V883 Ori are good sources to study the ice composition of planet forming material through thermally desorbed complex molecules, which have proven to be hard to observe in more evolved protoplanetary disks. △ Less

Submitted 24 August, 2018; originally announced August 2018.

Comments: 8 pages, 4 figures, 1 table. Accepted for publication in ApJL

Revised SED of the triple protostellar system VLA 1623-2417

Authors: Nadia M. Murillo, Daniel Harsono, Melissa McClure, Shih-Ping Lai, Michiel R. Hogerheijde

Abstract: VLA 1623$-$2417 is a triple protostellar system deeply embedded in Ophiuchus A. Sources A and B have a separation of 1.1", making their study difficult beyond the submillimeter regime. Lack of circumstellar gas emission suggested that VLA 1623$-$2417 B has a very cold envelope and is much younger than source A, generally considered the prototypical Class 0 source. We explore the consequences of ne… ▽ More VLA 1623$-$2417 is a triple protostellar system deeply embedded in Ophiuchus A. Sources A and B have a separation of 1.1", making their study difficult beyond the submillimeter regime. Lack of circumstellar gas emission suggested that VLA 1623$-$2417 B has a very cold envelope and is much younger than source A, generally considered the prototypical Class 0 source. We explore the consequences of new ALMA Band 9 data on the spectral energy distribution (SED) of VLA 1623$-$2417 and their inferred nature. Using dust continuum observations spanning from centimeter to near-infrared wavelengths, the SED of each component in VLA 1623$-$2417 is constructed and analysed. The ALMA Band 9 data presented here show that the SED of VLA 1623$-$2417 B does not peak at 850 $μ$m as previously expected, but instead presents the same shape as VLA 1623$-$2417 A at wavelengths shorter than 450 $μ$m. The results presented here indicate that the previous assumption that the flux in $Herschel$ and Spitzer observations is solely dominated by VLA 1623$-$2417 A is not valid, and instead, VLA 1623$-$2417 B most likely contributes a significant fraction of the flux at $λ~<$ 450 $μ$m. These results, however, do not explain the lack of circumstellar gas emission and puzzling nature of VLA 1623$-$2417 B. △ Less

Submitted 2 July, 2018; originally announced July 2018.

Comments: 6 pages, 2 figures, accepted to A&A letters to the editor

Journal ref: A&A 615, L14 (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

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

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

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

Submitted 31 May, 2018; originally announced May 2018.

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