Showing 1–20 of 20 results for author: Trapman, L

High-contrast JWST-MIRI spectroscopy of planet-forming disks for the JDISC Survey

Authors: Klaus M. Pontoppidan, Colette Salyk, Andrea Banzatti, Ke Zhang, Ilaria Pascucci, Karin I. Oberg, Feng Long, Carlos Munoz-Romero, John Carr, Joan Najita, Geoffrey A. Blake, Nicole Arulanantham, Sean Andrews, Nicholas P. Ballering, Edwin Bergin, Jenny Calahan, Douglas Cobb, Maria Jose Colmenares, Annie Dickson-Vandervelde, Anna Dignan, Joel Green, Phoebe Heretz, Greg Herczeg, Anusha Kalyaan, Sebastian Krijt , et al. (4 additional authors not shown)

Abstract: The JWST Disk Infrared Spectral Chemistry Survey (JDISCS) aims to understand the evolution of the chemistry of inner protoplanetary disks using the Mid-InfraRed Instrument (MIRI) on the James Webb Space Telescope (JWST). With a growing sample of >30 disks, the survey implements a custom method to calibrate the MIRI Medium Resolution Spectrometer (MRS) to contrasts of better than 1:300 across its 4… ▽ More The JWST Disk Infrared Spectral Chemistry Survey (JDISCS) aims to understand the evolution of the chemistry of inner protoplanetary disks using the Mid-InfraRed Instrument (MIRI) on the James Webb Space Telescope (JWST). With a growing sample of >30 disks, the survey implements a custom method to calibrate the MIRI Medium Resolution Spectrometer (MRS) to contrasts of better than 1:300 across its 4.9-28 micron spectral range. This is achieved using observations of Themis-family asteroids as precise empirical reference sources. High spectral contrast enables precise retrievals of physical parameters, searches for rare molecular species and isotopologues, and constraints on the inventories of carbon- and nitrogen-bearing species. JDISCS also offers significant improvements to the MRS wavelength and resolving power calibration. We describe the JDISCS calibrated data and demonstrate its quality using observations of the disk around the solar-mass young star FZ Tau. The FZ Tau MIRI spectrum is dominated by strong emission from warm water vapor. We show that the water and CO line emission originates from the disk surface and traces a range of gas temperatures of ~500-1500 K. We retrieve parameters for the observed CO and H2O lines, and show that they are consistent with a radial distribution represented by two temperature components. A high water abundance of n(H2O)~10^-4 fills the disk surface at least out to the 350 K isotherm at 1.5 au. We search the FZ Tau environs for extended emission detecting a large (radius of ~300 au) ring of emission from H2 gas surrounding FZ Tau, and discuss its origin. △ Less

Submitted 16 January, 2024; v1 submitted 28 November, 2023; originally announced November 2023.

Comments: Accepted for publication in the Astrophysical Journal

How large is a disk -- what do protoplanetary disk gas sizes really mean?

Authors: Leon Trapman, Giovanni Rosotti, Ke Zhang, Benoit Tabone

Abstract: It remains unclear what mechanism is driving the evolution of protoplanetary disks. Direct detection of the main candidates, either turbulence driven by magnetorotational instability or magnetohydrodynamical disk winds, has proven difficult, leaving the time evolution of the disk size as one of the most promising observables able to differentiate between these two mechanisms. But to do so successf… ▽ More It remains unclear what mechanism is driving the evolution of protoplanetary disks. Direct detection of the main candidates, either turbulence driven by magnetorotational instability or magnetohydrodynamical disk winds, has proven difficult, leaving the time evolution of the disk size as one of the most promising observables able to differentiate between these two mechanisms. But to do so successfully, we need to understand what the observed gas disk size actually traces. We studied the relation between $R_{\rm CO,\ 90\%}$, the radius that encloses 90% of the $^{12}$CO flux, and $R_c$, the radius that encodes the physical disk size, in order to provide simple prescriptions for conversions between these two sizes. For an extensive grid of thermochemical models we calculate $R_{\rm CO,\ 90\%}$ from synthetic observations and relate properties measured at this radius, such as the gas column density, to bulk disk properties, such as $R_c$ and the disk mass $M_{\rm disk}$. We found an empirical correlation between the gas column density at $R_{\rm CO,\ 90\%}$ and disk mass: $N_{\rm gas}(R_{\rm CO,\ 90\%}) \approx 3.73\times10^{21}(M_{\rm disk}/\mathrm{M}_{\odot})^{0.34}\ \mathrm{cm}^{-2}$. Using this correlation we derive an analytical prescription of $R_{\rm CO,\ 90\%}$ that only depends on $R_c$ and $M_{\rm disk}$. We derive $R_c$ for disks in Lupus, Upper Sco, Taurus and DSHARP, finding that disks in the older Upper Sco region are significantly smaller ($\langle R_c \rangle$ = 4.8 au) than disks in the younger Lupus and Taurus regions ($\langle R_c \rangle$ = 19.8 and 20.9 au, respectively). This temporal decrease in $R_c$ goes against predictions of both viscous and wind-driven evolution, but could be a sign of significant external photoevaporation having truncated disks in Upper Sco. △ Less

Submitted 14 July, 2023; originally announced July 2023.

Comments: 23 pages, 12 figures, accepted in ApJ

Substructures in Compact Disks of the Taurus Star-forming Region

Authors: Shangjia Zhang, Matt Kalscheur, Feng Long, Ke Zhang, Deryl E. Long, Edwin A. Bergin, Zhaohuan Zhu, Leon Trapman

Abstract: Observations of substructure in protoplanetary disks have largely been limited to the brightest and largest disks, excluding the abundant population of compact disks which are likely sites of planet formation. Here, we reanalyze ~0.1'', 1.33 mm ALMA continuum observations of 12 compact protoplanetary disks in the Taurus star-forming region. By fitting visibilities directly, we identify substructur… ▽ More Observations of substructure in protoplanetary disks have largely been limited to the brightest and largest disks, excluding the abundant population of compact disks which are likely sites of planet formation. Here, we reanalyze ~0.1'', 1.33 mm ALMA continuum observations of 12 compact protoplanetary disks in the Taurus star-forming region. By fitting visibilities directly, we identify substructures in 6 of the 12 compact disks. We then compare the substructures identified in the full Taurus sample of 24 disks in single star systems and the ALMA DSHARP survey, differentiating between compact (R_eff,90% < 50 au) and extended (R_eff,90% > 50 au) disk sources. We find that substructures are detected at nearly all radii in both small and large disks. Tentatively, we find fewer wide gaps in intermediate-sized disks with R_eff,90% between 30 and 90 au. We perform a series of planet-disk interaction simulations to constrain the sensitivity of our visibility-fitting approach. Under an assumption of planet-disk interaction, we use the gap widths and common disk parameters to calculate potential planet masses within the Taurus sample. We find that the young planet occurrence rate peaks near Neptune masses, similar to the DSHARP sample. For 0.01 $M_J/M_\odot$ $\lesssim$ $M_p/M_*$ $\lesssim$ 0.1 $M_J/M_\odot$, the rate is 17.4$\pm$8.3%; for 0.1 $M_J/M_\odot$ $\lesssim$ $M_p/M_*$ $\lesssim$ 1 $M_J/M_\odot$, it is 27.8$\pm$8.3%. Both of them are consistent with microlensing surveys. For gas giants more massive than 5 $M_J$, the occurrence rate is 4.2$\pm$4.2%, consistent with direct imaging surveys. △ Less

Submitted 9 May, 2023; v1 submitted 5 May, 2023; originally announced May 2023.

Comments: 30 pages, 15 figures, 6 tables, accepted to be published in The Astrophysical Journal

An analytic solution to measure the gas size in protoplanetary discs in the viscous self-similar scenario

Authors: Claudia Toci, Giuseppe Lodato, Francesco Gerardo Livio, Giovanni Rosotti, Leon Trapman

Abstract: In order to understand which mechanism is responsible for accretion in protoplanetary discs, a robust knowledge of the observed disc radius using gas tracers such as $^{12}$CO and other CO isotopologues is pivotal. Indeed, the two main theories proposed, viscous accretion and wind-driven accretion, predict different time evolution for the disc radii. In this Letter, we present an analytical soluti… ▽ More In order to understand which mechanism is responsible for accretion in protoplanetary discs, a robust knowledge of the observed disc radius using gas tracers such as $^{12}$CO and other CO isotopologues is pivotal. Indeed, the two main theories proposed, viscous accretion and wind-driven accretion, predict different time evolution for the disc radii. In this Letter, we present an analytical solution for the evolution of the disc radii in viscously evolving protoplanetary discs using $^{12}$CO as a tracer, under the assumption that the $^{12}$CO radius is the radius where the surface density of the disc is equal to the threshold for CO photo-dissociation. We discuss the properties of the solution and the limits of its applicability as a simple numerical prescription to evaluate the observed disc radii of populations of discs. Our results suggest that, in addition to photo-dissociation, also freeze out plays an important role in setting the disc size. We find an effective reduction of the CO abundance by about two orders of magnitude at the location of CO photo-dissociation, which however should not be interpreted as the bulk abundance of CO in the disc. The use of our analytical solution will allow to compute disc sizes for large quantities of models without using expensive computational resources such as radiative transfer calculations. △ Less

Submitted 10 November, 2022; v1 submitted 26 October, 2022; originally announced October 2022.

Comments: 5 pages, 3 figures

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

A novel way of measuring the gas disk mass of protoplanetary disks using N2H+ and C18O

Authors: L. Trapman, K. Zhang, M. R. L. van 't Hoff, M. R. Hogerheijde, E. A. Bergin

Abstract: Measuring the gas mass of protoplanetary disks, the reservoir available for giant planet formation, has proven to be difficult. We currently lack a far-infrared observatory capable of observing HD, and the most common gas mass tracer, CO, suffers from a poorly constrained CO-to-H$_2$ ratio. Expanding on previous work, we investigate if N2H+, a chemical tracer of CO poor gas, can be used to observa… ▽ More Measuring the gas mass of protoplanetary disks, the reservoir available for giant planet formation, has proven to be difficult. We currently lack a far-infrared observatory capable of observing HD, and the most common gas mass tracer, CO, suffers from a poorly constrained CO-to-H$_2$ ratio. Expanding on previous work, we investigate if N2H+, a chemical tracer of CO poor gas, can be used to observationally measure the CO-to-H$_2$ ratio and correct CO-based gas masses. Using disk structures obtained from the literature, we set up thermochemical models for three disks, TW Hya, DM Tau and GM Aur, to examine how well the CO-to-H$_2$ ratio and gas mass can be measured from N2H+ and C18O line fluxes. Furthermore, we compare these gas masses to independently gas masses measured from archival HD observations. The N2H+ (3-2)/C18O (2-1) line ratio scales with the disk CO-to-H$_2$ ratio. Using these two lines, we measure 4.6e-3 Msun < Mdisk < 1.1e-1 Msun for TW Hya, 1.5e-2 Msun < Mdisk < 9.6e-2 Msun for GM Aur and 3.1e-2 Msun < Mdisk < 9.6e-2 Msun for DM Tau. These gas masses agree with values obtained from HD within their respective uncertainties. The uncertainty on the N2H+ + C18O gas mass can be reduced by observationally constraining the cosmic ray ionization rate in disks. These results demonstrate the potential of using the combination of N2H+ and C18O to measure gas masses of protoplanetary disks. △ Less

Submitted 24 January, 2022; originally announced January 2022.

Comments: 10 pages, 5 figures, accepted for publication in ApJL

Effect of MHD wind-driven disk evolution on the observed sizes of protoplanetary disks

Authors: Leon Trapman, Benoit Tabone, Giovanni Rosotti, Ke Zhang

Abstract: It is still unclear whether the evolution of protoplanetary disks, a key ingredient in the theory of planet formation, is driven by viscous turbulence or magnetic disk winds. As viscously evolving disks expand outward over time, the evolution of disk sizes is a discriminant test for studying disk evolution. However, it is unclear how the observed disk size changes over time if disk evolution is dr… ▽ More It is still unclear whether the evolution of protoplanetary disks, a key ingredient in the theory of planet formation, is driven by viscous turbulence or magnetic disk winds. As viscously evolving disks expand outward over time, the evolution of disk sizes is a discriminant test for studying disk evolution. However, it is unclear how the observed disk size changes over time if disk evolution is driven by magnetic disk winds. Combining the thermochemical code DALI with the analytical wind-driven disk evolution model presented in Tabone et al. (2021a), we study the time evolution of the observed gas outer radius as measured from CO rotational emission ($R_{\rm CO, 90\%}$). The evolution of $R_{\rm CO, 90\%}$ is driven by the evolution of the disk mass, as the physical radius stays constant over time. For a constant $α_{\rm DW}$, an extension of the $α-$Shakura-Sunyaev parameter to wind-driven accretion, $R_{\rm CO, 90\%}$ decreases linearly with time. Its initial size is set by the disk mass and the characteristic radius $R_c$, but only $R_c$ affects the evolution of $R_{\rm CO, 90\%}$, with a larger $R_c$ resulting in a steeper decrease of $R_{\rm CO, 90\%}$. For a time-dependent $α_{\rm DW}$ $R_{\rm CO, 90\%}$ stays approximately constant during most of the disk lifetime until $R_{\rm CO, 90\%}$ rapidly shrinks as the disk dissipates. The constant $α_{\rm DW}$-models are able to reproduce the observed gas disk sizes in the $\sim1-3$ Lupus and $\sim5-11$ Myr old Upper Sco star-forming regions. However, they likely overpredict the gas disk size of younger $(\lessapprox0.7\ \mathrm{Myr})$ disks. △ Less

Submitted 1 December, 2021; originally announced December 2021.

Comments: 22 pages, 16 figures, accepted in ApJ

Molecules with ALMA at Planet-forming Scales (MAPS) V: CO gas distributions

Authors: Ke Zhang, Alice S. Booth, Charles J. Law, Arthur D. Bosman, Kamber R. Schwarz, Edwin A. Bergin, Karin I. Öberg, Sean M. Andrews, Viviana V. Guzmán, Catherine Walsh, Chunhua Qi, Merel L. R. van 't Hoff, Feng Long, David J. Wilner, Jane Huang, Ian Czekala, John D. Ilee, Gianni Cataldi, Jennifer B. Bergner, Yuri Aikawa, Richard Teague, Jaehan Bae, Ryan A. Loomis, Jenny K. Calahan, Felipe Alarcón , et al. (10 additional authors not shown)

Abstract: Here we present high resolution (15-24 au) observations of CO isotopologue lines from the Molecules with ALMA on Planet-forming Scales (MAPS) ALMA Large Program. Our analysis employs $^{13}$CO and C$^{18}$O ($J$=2-1), (1-0), and C$^{17}$O (1-0) line observations of five protoplanetary disks. We retrieve CO gas density distributions, using three independent methods: (1) a thermo-chemical modeling f… ▽ More Here we present high resolution (15-24 au) observations of CO isotopologue lines from the Molecules with ALMA on Planet-forming Scales (MAPS) ALMA Large Program. Our analysis employs $^{13}$CO and C$^{18}$O ($J$=2-1), (1-0), and C$^{17}$O (1-0) line observations of five protoplanetary disks. We retrieve CO gas density distributions, using three independent methods: (1) a thermo-chemical modeling framework based on the CO data, the broadband spectral energy distribution, and the mm-continuum emission; (2) an empirical temperature distribution based on optically thick CO lines; and (3) a direct fit to the C$^{17}$O hyperfine lines. Results from these methods generally show excellent agreement. The CO gas column density profiles of the five disks show significant variations in the absolute value and the radial shape. Assuming a gas-to-dust mass ratio of 100, all five disks have a global CO-to-H$_2$ abundance of 10-100 times lower than the ISM ratio. The CO gas distributions between 150-400 au match well with models of viscous disks, supporting the long-standing theory. CO gas gaps appear to be correlated with continuum gap locations, but some deep continuum gaps do not have corresponding CO gaps. The relative depths of CO and dust gaps are generally consistent with predictions of planet-disk interactions, but some CO gaps are 5-10 times shallower than predictions based on dust gaps. This paper is part of the MAPS special issue of the Astrophysical Journal Supplement. △ Less

Submitted 23 September, 2021; v1 submitted 13 September, 2021; originally announced September 2021.

Comments: This paper is part of the MAPS special issue of the Astrophysical Journal Supplement. 36 pages, 21 figures, accepted for publication in ApJS

On the secular evolution of the ratio between gas and dust radii in protoplanetary discs

Authors: Claudia Toci, Giovanni Rosotti, Giuseppe Lodato, Leonardo Testi, Leon Trapman

Abstract: A key problem in protoplanetary disc evolution is understanding the efficiency of dust radial drift. This process makes the observed dust disc sizes shrink on relatively short timescales, implying that discs started much larger than what we see now. In this paper we use an independent constraint, the gas radius (as probed by CO rotational emission), to test disc evolution models. In particular, we… ▽ More A key problem in protoplanetary disc evolution is understanding the efficiency of dust radial drift. This process makes the observed dust disc sizes shrink on relatively short timescales, implying that discs started much larger than what we see now. In this paper we use an independent constraint, the gas radius (as probed by CO rotational emission), to test disc evolution models. In particular, we consider the ratio between the dust and gas radius, $R_{\rm CO}/R_{\rm dust}$. We model the time evolution of protoplanetary discs under the influence of viscous evolution, grain growth, and radial drift. Then, using the radiative transfer code RADMC with approximate chemistry, we compute the dust and gas radii of the models and investigate how $R_{\rm CO}/R_{\rm dust}$ evolves. Our main finding is that, for a broad range of values of disc mass, initial radius, and viscosity, $R_{\rm CO}/R_{\rm dust}$ becomes large (>5) after only a short time (<1 Myr) due to radial drift. This is at odds with measurements in young star forming regions such as Lupus, which find much smaller values, implying that dust radial drift is too efficient in these models. Substructures, commonly invoked to stop radial drift in large, bright discs, must then be present, although currently unresolved, in most discs. △ Less

Submitted 21 July, 2021; originally announced July 2021.

Comments: 17 pages, 14 figures

CO isotopolog line fluxes of viscously evolving disks: cold CO conversion insufficient to explain observed low fluxes

Authors: Leon Trapman, Arthur D. Bosman, Giovanni Rosotti, Michiel R. Hogerheijde, Ewine F. van Dishoeck

Abstract: Protoplanetary disks are thought to evolve viscously, where the disk mass - the reservoir available for planet formation - decreases over time as material is accreted onto the central star. Observations show a correlation between dust mass and the stellar accretion rate, as expected from viscous theory. However, the gas mass inferred from 13CO and C18O line fluxes, which should be a more direct me… ▽ More Protoplanetary disks are thought to evolve viscously, where the disk mass - the reservoir available for planet formation - decreases over time as material is accreted onto the central star. Observations show a correlation between dust mass and the stellar accretion rate, as expected from viscous theory. However, the gas mass inferred from 13CO and C18O line fluxes, which should be a more direct measure, shows no such correlation. Using thermochemical DALI models, we investigate how 13CO and C18O J=3-2 line fluxes change over time in a viscously evolving disk. We also investigate if the chemical conversion of CO through grain-surface chemistry combined with viscous evolution can explain the observations of disks in Lupus. The 13CO and C18O 3-2 line fluxes increase over time due to their optically thick emitting regions growing in size as the disk expands viscously. The C18O 3-2 emission is optically thin throughout the disk for only a subset of our models (Mdisk (t = 1 Myr) < 1e-3 Msun). For these disks the integrated C18O flux decreases with time, similar to the disk mass. The C18O 3-2 fluxes for the bulk of the disks in Lupus (with Mdust < 5e-5 Msun) can be reproduced to within a factor of ~2 with viscously evolving disks in which CO is converted into other species through grain-surface chemistry driven by a cosmic-ray ionization rate zeta_cr ~ 5e-17 - 1e-16 s^-1. However, explaining the stacked C18O upper limits requires a lower average abundance than our models can produce and they cannot explain the observed 13CO fluxes, which, for most disks, are more than an order of magnitude fainter than what our models predict. Reconciling the 13CO fluxes of viscously evolving disks with the observations requires either a combination of efficient vertical mixing and a high zeta_cr or low mass disks (Mdust < 3e-5 Msun) being much thinner and/or smaller than their more massive counterparts. △ Less

Submitted 9 March, 2021; originally announced March 2021.

Comments: 21 pages, 14 figures, accepted in A&A

Journal ref: A&A 649, A95 (2021)

Chemically tracing the water snowline in protoplanetary disks with HCO$^+$

Authors: M. Leemker, M. L. R. van 't Hoff, L. Trapman, M. L. van Gelder, M. R. Hogerheijde, D. Ruíz-Rodríguez, E. F. van Dishoeck

Abstract: [Abridged] Planet formation is expected to be enhanced around snowlines in protoplanetary disks, in particular around the water snowline. However, the close proximity of the water snowline to the host star and water in the Earth's atmosphere makes a direct detection of the water snowline in protoplanetary disks challenging. Following earlier work on protostellar envelopes, the aim of this research… ▽ More [Abridged] Planet formation is expected to be enhanced around snowlines in protoplanetary disks, in particular around the water snowline. However, the close proximity of the water snowline to the host star and water in the Earth's atmosphere makes a direct detection of the water snowline in protoplanetary disks challenging. Following earlier work on protostellar envelopes, the aim of this research is to investigate the validity of HCO$^+$ and H$^{13}$CO$^+$, as tracers of the water snowline in protoplanetary disks, as HCO$^+$ is destroyed by gas-phase water. Two small chemical networks are used to predict the HCO$^+$ abundance in a typical Herbig Ae disk. Subsequently, the corresponding emission profiles are modelled for H$^{13}$CO$^+$ and HCO$^+$ $J=2-1$, which provides the best balance between brightness and optical depth effects of the continuum emission. The HCO$^+$ abundance jumps by two orders of magnitude just outside the water snowline at 4.5 AU. We find that the emission of H$^{13}$CO$^+$ and HCO$^+$ is ring-shaped due to three effects: destruction of HCO$^+$ by gas-phase water, continuum optical depth, and molecular excitation effects. The presence of gas-phase water causes an additional drop of only $\sim$13% and 24% in the center of the disk, for H$^{13}$CO$^+$ and HCO$^+$, respectively. For the much more luminous outbursting source V883Ori, our models predict that the effect of dust and excitation are not limiting if the snowline is located outside $\sim$40 AU. Our analysis of ALMA observations of HCO$^+$ $J=3-2$ is consistent with the water snowline located around 100 AU. The HCO$^+$ abundance drops steeply around the water snowline, but dust and excitation can conceal the drop in HCO$^+$ emission due to the water snowline. Therefore, locating the water snowline with HCO$^+$ in Herbig disks is very difficult, but it is possible for outbursting sources like V883Ori. △ Less

Submitted 24 November, 2020; originally announced November 2020.

Comments: Accepted for publication in Astronomy and Astrophysics. 13 pages, 11 figures + appendices

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

Observed sizes of planet-forming disks trace viscous evolution

Authors: L. Trapman, G. Rosotti, A. D. Bosman, M. R. Hogerheijde, E. F. van Dishoeck

Abstract: The evolution of protoplanetary disks is dominated by the conservation of angular momentum, where the accretion of material onto the central star is driven by viscous expansion of the outer disk or by disk winds extracting angular momentum without changing the disk size. Studying the time evolution of disk sizes allows us therefore to distinguish between viscous stresses or disk winds as the main… ▽ More The evolution of protoplanetary disks is dominated by the conservation of angular momentum, where the accretion of material onto the central star is driven by viscous expansion of the outer disk or by disk winds extracting angular momentum without changing the disk size. Studying the time evolution of disk sizes allows us therefore to distinguish between viscous stresses or disk winds as the main mechanism of disk evolution. Observationally, estimates of the disk gaseous outer radius are based on the extent of the CO rotational emission, which, during the evolution, is also affected by the changing physical and chemical conditions in the disk. We use physical-chemical DALI models to study how the extent of the CO emission changes with time in a viscously expanding disk and investigate to what degree this observable gas outer radius is a suitable tracer of viscous spreading and whether current observations are consistent with viscous evolution. We find that the gas outer radius (R_co) measured from our models matches the expectations of a viscously spreading disk: R_co increases with time and for a given time R_co is larger for a disk with a higher viscosity alpha_visc. However, in the extreme case where the disk mass is low (less than 10^-4 Msun) and alpha_visc is high (larger than 10^-2), R_co will instead decrease with time as a result of CO photodissociation in the outer disk. For most disk ages R_co is up to 12x larger than the characteristic size R_c of the disk, and R_co/R_c is largest for the most massive disk. As a result of this difference, a simple conversion of R_co to alpha_visc will overestimate the true alpha_visc of the disk by up to an order of magnitude. We find that most observed gas outer radii in Lupus can be explained using a viscously evolving disk that starts out small (R_c = 10 AU) and has a low viscosity (alpha_visc = 10^-4 - 10^-3). △ Less

Submitted 22 May, 2020; originally announced May 2020.

Comments: 19 pages, 17 figures, accepted in A&A

Journal ref: A&A 640, A5 (2020)

Constraining the radial drift of millimeter-sized grains in the protoplanetary disks in Lupus

Authors: L. Trapman, M. Ansdell, M. R. Hogerheijde, S. Facchini, C. F. Manara, A. Miotello, J. P. Williams, S. Bruderer

Abstract: Recent ALMA surveys of protoplanetary disks have shown that for most disks the extent of the gas emission is greater than the extent of the thermal emission of the millimeter-sized dust. Both line optical depth and the combined effect of radially dependent grain growth and radial drift may contribute to this observed effect. For a sample of 10 disks from the Lupus survey we investigate how well du… ▽ More Recent ALMA surveys of protoplanetary disks have shown that for most disks the extent of the gas emission is greater than the extent of the thermal emission of the millimeter-sized dust. Both line optical depth and the combined effect of radially dependent grain growth and radial drift may contribute to this observed effect. For a sample of 10 disks from the Lupus survey we investigate how well dust-based models without radial dust evolution reproduce the observed 12CO outer radius, and determine whether radial dust evolution is required to match the observed gas-dust size difference. We used the thermochemical code DALI to obtain 12CO synthetic emission maps and measure gas and dust outer radii (Rco, Rmm) using the same methods as applied to the observations, which were compared to observations on a source-by-source basis. For 5 disks we find that the observed gas-dust size difference is larger than the gas-dust size difference due to optical depth, indicating that we need both dust evolution and optical depth effects to explain the observed gas-dust size difference. For the other 5 disks the observed gas-dust size difference can be explained using only line optical depth effects. We also identify 6 disks not included in our initial sample but part of a survey of the same star-forming region that show significant 12CO emission beyond 4 x Rmm. These disks, for which no Rco is available, likely have gas-dust size differences greater than 4 and are difficult to explain without substantial dust evolution. Our results suggest that radial drift and grain growth are common features among both bright and fain disks. The effects of radial drift and grain growth can be observed in disks where the dust and gas radii are significantly different, while more detailed models and deeper observations are needed to see this effect in disks with smaller differences. △ Less

Submitted 15 April, 2020; originally announced April 2020.

Comments: 17 pages, 11 figures, accepted in A&A

Journal ref: A&A 638, A38 (2020)

Mass constraints for 15 protoplanetary disks from HD 1-0

Authors: M. Kama, L. Trapman, D. Fedele, S. Bruderer, M. R. Hogerheijde, A. Miotello, E. F. van Dishoeck, C. Clarke, E. A. Bergin

Abstract: Hydrogen deuteride (HD) rotational line emission can provide reliable protoplanetary disk gas mass measurements, but it is difficult to observe and detections have been limited to three T-Tauri disks. No new data have been available since the \emph{Herschel} Space Observatory mission ended in 2013. We set out to obtain new disk gas mass constraints by analysing upper limits on HD 1-0 emission in \… ▽ More Hydrogen deuteride (HD) rotational line emission can provide reliable protoplanetary disk gas mass measurements, but it is difficult to observe and detections have been limited to three T-Tauri disks. No new data have been available since the \emph{Herschel} Space Observatory mission ended in 2013. We set out to obtain new disk gas mass constraints by analysing upper limits on HD 1-0 emission in \emph{Herschel}/PACS archival data from the DIGIT key programme. With a focus on the Herbig Ae/Be disks, whose stars are more luminous than T Tauris, we determine upper limits for HD in data previosly analysed for its line detections. Their significance is studied with a grid of models run with the DALI physical-chemical code, customised to include deuterium chemistry. Nearly all the disks are constrained to $M_{\rm gas}\leq0.1\,$M$_{\odot}$, ruling out global gravitational instability. A strong constraint is obtained for the HD 163296 disk mass, $M_{\rm gas} \leq 0.067\,$M$_{\odot}$, implying $Δ_{\rm g/d}\leq100$. This HD-based mass limit is towards the low end of CO-based mass estimates for the disk, highlighting the large uncertainty in using only CO and suggesting that gas-phase CO depletion in HD 163296 is at most a factor of a few. The $M_{\rm gas}$ limits for HD 163296 and HD 100546, both bright disks with massive candidate protoplanetary systems, suggest disk-to-planet mass conversion efficiencies of $M_{\rm p}/(M_{\rm gas} + M_{\rm p})\approx10$ to $40\,$% for present-day values. Near-future observations with SOFIA/HIRMES will be able to detect HD in the brightest Herbig~Ae/Be disks within $150\,$pc with $\approx10\,$h integration time. △ Less

Submitted 26 December, 2019; originally announced December 2019.

Comments: Accepted for publication in A&A

Journal ref: A&A 634, A88 (2020)

Gas vs dust sizes of protoplanetary disks: effects of dust evolution

Authors: L. Trapman, S. Facchini, M. R. Hogerheijde, E. F. van Dishoeck, S. Bruderer

Abstract: The extent of the gas in protoplanetary disks is observed to be universally larger than the extent of the dust. This is often attributed to radial drift and grain growth of the mm grains, but line optical depth produces a similar observational signature. We investigate in what parts of the disk structure parameter space dust evolution and line optical depth are the dominant drivers of the observed… ▽ More The extent of the gas in protoplanetary disks is observed to be universally larger than the extent of the dust. This is often attributed to radial drift and grain growth of the mm grains, but line optical depth produces a similar observational signature. We investigate in what parts of the disk structure parameter space dust evolution and line optical depth are the dominant drivers of the observed gas and dust size difference. Using the thermochemical model DALI with dust evolution included we ran a grid of models aimed at reproducing the observed gas and dust size dichotomy. The relation between R_dust and dust evolution is non-monotonic and depends on the disk structure. R_gas is directly related to the radius where the CO column density drops below 10^15 cm^-2 and CO becomes photodissociated. R_gas is not affected by dust evolution but scales with the total CO content of the disk. R_gas/R_dust > 4 is a clear sign for dust evolution and radial drift in disks, but these cases are rare in current observations. For disks with a smaller R_gas/R_dust, identifying dust evolution from R_gas/R_dust requires modelling the disk structure including the total CO content. To minimize the uncertainties due to observational factors requires FWHM_beam < 1x the characteristic radius and a peak SNR > 10 on the 12CO emission moment zero map. For the dust outer radius to enclose most of the disk mass, it should be defined using a high fraction (90-95%) of the total flux. For the gas, any radius enclosing > 60% of the 12CO flux will contain most of the disk mass. To distinguish radial drift and grain growth from line optical depth effects based on size ratios requires disks to be observed at high enough angular resolution and the disk structure should to be modelled to account for the total CO content of the disk. △ Less

Submitted 14 March, 2019; originally announced March 2019.

Comments: 18 pages, 27 figures, accepted in A&A

The ALMA Lupus protoplanetary disk survey: evidence for compact gas disks and molecular rings from CN

Authors: S. E. van Terwisga, E. F. van Dishoeck, P. Cazzoletti, S. Facchini, L. Trapman, J. P. Williams, C. F. Manara, A. Miotello, N. van der Marel, M. Ansdell, M. R. Hogerheijde, M. Tazzari, L. Testi

Abstract: The cyanide radical CN is abundant in protoplanetary disks, with line fluxes often comparable to those of $^{13}$CO. It is known to be sensitive to UV irradiation of the upper disk atmosphere, with models predicting ring-shaped emission. We seek to characterize the CN emission from 94 Class II disks in the Lupus star forming region, compare it to observations in other regions, and interpret our ob… ▽ More The cyanide radical CN is abundant in protoplanetary disks, with line fluxes often comparable to those of $^{13}$CO. It is known to be sensitive to UV irradiation of the upper disk atmosphere, with models predicting ring-shaped emission. We seek to characterize the CN emission from 94 Class II disks in the Lupus star forming region, compare it to observations in other regions, and interpret our observations with a grid of models. The CN emission morphology is discussed for two primordial disks, Sz 71 and Sz 98, and modeled in more detail. ALMA observed CN N=3-2 in Lupus disks down to sensitivities better than previous surveys. Models constructed with the physico-chemical code DALI are used to study the integrated fluxes of the disks and resolved emission of CN in disks without (dust) substructures. We find that CN N=3-2 is bright, and detected in 38% of sources, but its disk-integrated flux is not strongly correlated to either $^{13}$CO or continuum flux. Compared to pre-ALMA single-dish surveys, no significant difference in the CN flux distributions in Lupus and Taurus-Auriga is found, although $ρ$ Ophiuchus disks may be fainter on average. We find ring-shaped CN emission with peak radii of ~50 AU in two resolved disks. We conclude that a large fraction of sources is faint in CN; only exponential gas surface density cutoffs at $R_{c} \leq 15$ AU can reconcile observations with models. This is the first observational evidence of such a compact gas disk population in Lupus. Absolute intensities and the emission morphology of CN are reproduced by DALI models without the need for any continuum substructure; they are unrelated to the CO snowline location. These observations and the successful modeling of these rings provide a new probe of the structure and conditions in disks, and particularly their incident UV radiation field, if disk size is determined by the data. △ Less

Submitted 7 November, 2018; originally announced November 2018.

Comments: 22 pages, 14 figures, A&A accepted

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

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

V1094 Sco: a rare giant multi-ringed disk around a T Tauri star

Authors: S. E. van Terwisga, E. F. van Dishoeck, M. Ansdell, N. van der Marel, L. Testi, J. P. Williams, S. Facchini, M. Tazzari, M. R. Hogerheijde, L. Trapman, C. F. Manara, A. Miotello, L. T. Maud, D. Harsono

Abstract: A wide variety of ring-like dust structures has been detected in protoplanetary disks, but their origin and frequency are still unclear. We characterize the structure of an extended, multi-ringed disk discovered serendipitously in the ALMA Lupus disk survey and put it in the context of the Lupus disk population. ALMA observations in Band 6 at 234 GHz and Band 7 at 328 GHz at 0.3" resolution toward… ▽ More A wide variety of ring-like dust structures has been detected in protoplanetary disks, but their origin and frequency are still unclear. We characterize the structure of an extended, multi-ringed disk discovered serendipitously in the ALMA Lupus disk survey and put it in the context of the Lupus disk population. ALMA observations in Band 6 at 234 GHz and Band 7 at 328 GHz at 0.3" resolution toward the K6 star V1094 Sco in Lupus III are presented, and its disk structure is analyzed. The spectral index $α_{mm}$ is determined in the inner 150 AU of the disk. The ALMA continuum data show a very extended disk with two gap/ring pairs. The gaps are located at 100 AU and 170 AU, the bright rings at 130 AU and 220 AU. Continuum emission is detected out to a 300 AU distance, similar to IM Lup but a factor of 5 larger than typically found for Lupus disks at this sensitivity and resolution. The bright central region of the disk (within 35 AU) is possibly optically thick at 1 mm wavelengths, and has a brightness temperature of only 13 K. The spectral index increases between the inner disk and the first ring, at the location of the first gap. Due to the low temperature of the disk midplane, snow lines can be excluded as the drivers behind the ring and gap formation in this disk. Disks the size of V1094 Sco are rare, and only 2.1+-1.5% of disks in Lupus show continuum emission beyond 200 AU. Possible connections between the large primordial disk population, transition disks, and exoplanets are discussed. △ Less

Submitted 8 May, 2018; originally announced May 2018.

Comments: 12 pages, 7 figures. Accepted for publication in Astronomy & Astrophysics

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

ALMA Survey of Lupus Protoplanetary Disks II: Gas Disk Radii

Authors: M. Ansdell, J. P. Williams, L. Trapman, S. E. van Terwisga, S. Facchini, C. F. Manara, N. van der Marel, A. Miotello, M. Tazzari, M. Hogerheijde, G. Guidi, L. Testi, E. F. van Dishoeck

Abstract: We present ALMA Band 6 observations of a complete sample of protoplanetary disks in the young (1-3 Myr) Lupus star-forming region, covering the 1.33 mm continuum and the 12CO, 13CO, and C18O J=2-1 lines. The spatial resolution is 0.25 arcsec with a medium 3-sigma continuum sensitivity of 0.30 mJy, corresponding to M_dust ~ 0.2 M_earth. We apply "Keplerian masking" to enhance the signal-to-noise ra… ▽ More We present ALMA Band 6 observations of a complete sample of protoplanetary disks in the young (1-3 Myr) Lupus star-forming region, covering the 1.33 mm continuum and the 12CO, 13CO, and C18O J=2-1 lines. The spatial resolution is 0.25 arcsec with a medium 3-sigma continuum sensitivity of 0.30 mJy, corresponding to M_dust ~ 0.2 M_earth. We apply "Keplerian masking" to enhance the signal-to-noise ratios of our 12CO zero-moment maps, enabling measurements of gas disk radii for 22 Lupus disks; we find that gas disks are universally larger than mm dust disks by a factor of two on average, likely due to a combination of the optically thick gas emission as well as the growth and inward drift of the dust. Using the gas disk radii, we calculate the dimensionless viscosity parameter, alpha_visc, finding a broad distribution and no correlations with other disk or stellar parameters, suggesting that viscous processes have not yet established quasi-steady states in Lupus disks. By combining our 1.33 mm continuum fluxes with our previous 890 micron continuum observations, we also calculate the mm spectral index, alpha_mm, for 70 Lupus disks; we find an anti-correlation between alpha_mm and mm flux for low-mass disks (M_dust < 5), followed by a flattening as disks approach alpha_mm = 2, which could indicate faster grain growth in higher-mass disks, but may also reflect their larger optically thick components. In sum, this work demonstrates the continuous stream of new insights into disk evolution and planet formation that can be gleaned from unbiased ALMA disk surveys. △ Less

Submitted 3 May, 2018; v1 submitted 15 March, 2018; originally announced March 2018.

Comments: 22 pages, 10 figures; accepted to ApJ

HD far infrared emission as a measure of protoplanetary disk mass

Authors: Leon Trapman, Anna Miotello, Mihkel Kama, Ewine F. van Dishoeck, Simon Bruderer

Abstract: Protoplanetary disks around young stars are the sites of planet formation. While the dust mass can be estimated using standard methods, determining the gas mass - and thus the amount of material available to form giant planets - has proven to be very difficult. Hydrogen deuteride (HD) is a promising alternative to the commonly-used gas mass tracer, CO. We aim to examine the robustness of HD as tra… ▽ More Protoplanetary disks around young stars are the sites of planet formation. While the dust mass can be estimated using standard methods, determining the gas mass - and thus the amount of material available to form giant planets - has proven to be very difficult. Hydrogen deuteride (HD) is a promising alternative to the commonly-used gas mass tracer, CO. We aim to examine the robustness of HD as tracer of the disk gas mass, specifically the effect of gas mass on the HD FIR emission and its sensitivity to the vertical structure. Deuterium chemistry reactions relevant for HD were implemented in the thermochemical code DALI and models were run for a range of disk masses and vertical structures. The HD J=1-0 line intensity depends directly on the gas mass through a sublinear power law relation with a slope of ~0.8. Assuming no prior knowledge about the vertical structure of a disk and using only the HD 1-0 flux, gas masses can be estimated to within a factor of 2 for low mass disks (M$_{\rm disk} < 10^{-3}$ M$_\odot$). For more massive disks, this uncertainty increases to more than an order of magnitude. Adding the HD 2-1 line or independent information about the vertical structure can reduce this uncertainty to a factor of ~3 for all disk masses. For TW Hya, using the radial and vertical structure from Kama et al. 2016b the observations constrain the gas mass to $6\cdot10^{-3}$ M$_\odot$ < M$_{\rm disk} < 9\cdot10^{-3}$ M$_\odot$. Future observations require a 5$σ$ sensitivity of $1.8\cdot10^{-20}$ W m$^{-2}$ ($2.5\cdot10^{-20}$ W m$^{-2}$) and a spectral resolving power R > 300 (1000) to detect HD 1-0 (HD 2-1) for all disk masses above $10^{-5}$ M$_\odot$ with a line-to-continuum ratio > 0.01. These results show that HD can be used as an independent gas mass tracer with a relatively low uncertainty and should be considered as an important science goal for future FIR missions. △ Less

Submitted 22 May, 2017; originally announced May 2017.

Comments: 18 pages, 13 figures, accepted in A&A

Journal ref: A&A 605, A69 (2017)