The protoplanetary system HD 169142 is one of the few cases where a potential candidate protoplanet has recently been detected by direct imaging in the near-infrared. To study the interaction between ...the protoplanet and the disk itself, observations of the gas and dust surface density structure are needed. This paper reports new ALMA observations of the dust continuum at 1.3 mm, 12CO, 13CO, and C18O J = 2−1 emission from the system HD 169142 (which is observed almost face-on) at an angular resolution of ~\hbox{$0\farcs3 \,\times\, 0\farcs2$}0 .̋ 3 × 0 .̋ 2 (~35 × 20 au). The dust continuum emission reveals a double-ring structure with an inner ring between \hbox{$0\farcs17{-}0\farcs28$}0 .̋ 17−0 .̋ 28 (~20−35 au) and an outer ring between \hbox{$0\farcs48{-}0\farcs64$}0 .̋ 48−0 .̋ 64 (~56−83 au). The size and position of the inner ring is in good agreement with previous polarimetric observations in the near-infrared and is consistent with dust trapping by a massive planet. No dust emission is detected inside the inner dust cavity (R ≲ 20 au) or within the dust gap (~35−56 au) down to the noise level. In contrast, the channel maps of the J = 2−1 line of the three CO isotopologs reveal gas inside the dust cavity and dust gap. The gaseous disk is also much larger than the compact dust emission; it extends to ~1\hbox{$\farcs5$}.̋ 5 (~180 au) in radius. This difference and the sharp drop of the continuum emission at large radii point to radial drift of large dust grains (>μm size). Using the thermo-chemical disk code dali, we modeled the continuum and the CO isotopolog emission to quantitatively measure the gas and dust surface densities. The resulting gas surface density is reduced by a factor of ~30−40 inward of the dust gap. The gas and dust distribution indicate that two giant planets shape the disk structure through dynamical clearing (dust cavity and gap) and dust trapping (double-ring dust distribution).
We present a new code for solving the molecular and atomic excitation and radiation transfer problem in a molecular gas and predicting emergent spectra. This code works in arbitrary three dimensional ...geometry using unstructured Delaunay latices for the transport of photons. Various physical models can be used as input, ranging from analytical descriptions over tabulated models to SPH simulations. To generate the Delaunay grid we sample the input model randomly, but weigh the sample probability with the molecular density and other parameters, and thereby we obtain an average grid point separation that scales with the local opacity. Our code does photon very efficiently so that the slow convergence of opaque models becomes traceable. When convergence between the level populations, the radiation field, and the point separation has been obtained, the grid is ray-traced to produced images that can readily be compared to observations. Because of the high dynamic range in scales that can be resolved using this type of grid, our code is particularly well suited for modeling of ALMA data. Our code can furthermore deal with overlapping lines of multiple molecular and atomic species.
Context.
The evolution of protoplanetary disks is dominated by the conservation of angular momentum, where the accretion of material onto the central star is fed by the 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 therefore allows us to distinguish between viscous stresses or disk winds as the main mechanism of disk evolution. Observationally, estimates of the size of the gaseous disk are based on the extent of CO submillimeter rotational emission, which is also affected by the changing physical and chemical conditions in the disk during the evolution.
Aims.
We study how the gas outer radius measured from the extent of the CO emission changes with time in a viscously expanding disk. We also 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.
Methods.
For a set of observationally informed initial conditions we calculated the viscously evolved density structure at several disk ages and used the thermochemical code
DALI
to compute synthetic emission maps, from which we measured gas outer radii in a similar fashion as observations.
Results.
The gas outer radii (
R
CO, 90%
) measured from our models match the expectations of a viscously spreading disk:
R
CO, 90%
increases with time and, for a given time,
R
CO, 90%
is larger for a disk with a higher viscosity
α
visc
. However, in the extreme case in which the disk mass is low (
M
disk
≤ 10
−4
M
⊙
) and
α
visc
is high (≥10
−2
),
R
CO, 90%
instead decreases with time as a result of CO photodissociation in the outer disk. For most disk ages,
R
CO, 90%
is up to ~12× larger than the characteristic size
R
c
of the disk, and
R
CO, 90%
/
R
c
is largest for the most massive disk. As a result of this difference, a simple conversion of
R
CO, 90%
to
α
visc
overestimates the true
α
visc
of the disk by up to an order of magnitude. Based on our models, we find that most observed gas outer radii in Lupus can be explained using viscously evolving disks that start out small (
R
c
(
t
= 0) ≃ 10 AU) and have a low viscosity (
α
visc
= 10
−4
−10
−3
).
Conclusions.
Current observations are consistent with viscous evolution, but expanding the sample of observed gas disk sizes to star-forming regions, both younger and older, would better constrain the importance of viscous spreading during disk evolution.
We investigate the apparent discrepancy between gas and dust outer radii derived from millimeter observations of protoplanetary disks. Using 230 and 345 GHz continuum and CO image data from the ...Submillimeter Array for four nearby disk systems (HD 163296, TW Hydrae, GM Aurigae, and MWC 480), we examine models of circumstellar disk structure and the effects of their treatment of the outer disk edge. We show that for these disks, models described by power laws in surface density and temperature that are truncated at an outer radius are incapable of reproducing both the gas and dust emission simultaneously: the outer radius derived from the dust continuum emission is always significantly smaller than the extent of the molecular gas disk traced by CO emission. However, a simple model motivated by similarity solutions of the time evolution of accretion disks that includes a tapered exponential edge in the surface density distribution (and the same number of free parameters) does much better at reproducing both the gas and dust emission. While this analysis does not rule out the disparate radii implied by the truncated power law models, a realistic alternative disk model, grounded in the physics of accretion, provides a consistent picture for the extent of both the gas and dust.
Context. The extent of the gas in protoplanetary discs is observed to be universally larger than the extent of the dust. This is often attributed to radial drift and grain growth of the millimetre ...grains, but line optical depth produces a similar observational signature. Aims. We investigate in which parts of the disc structure parameter space dust evolution and line optical depth are the dominant drivers of the observed gas and dust size difference. Methods. 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. Results. The relation between Rdust and dust evolution is non-monotonic and depends on the disc structure. The quantity Rgas is directly related to the radius where the CO column density drops below 1015 cm−2 and CO becomes photodissociated; Rgas is not affected by dust evolution but scales with the total CO content of the disc. While these cases are rare in current observations, Rgas/Rdust > 4 is a clear sign of dust evolution and radial drift in discs. For discs with a smaller Rgas/Rdust, identifying dust evolution from Rgas/Rdust requires modelling the disc structure including the total CO content. To minimize the uncertainties due to observational factors requires FWHMbeam < 1× the characteristic radius and a peak S/N > 10 on the 12CO emission moment zero map. For the dust outer radius to enclose most of the disc 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 contains most of the disc mass. Conclusions. To distinguish radial drift and grain growth from line optical depth effects based on size ratios requires discs to be observed at high enough angular resolution and the disc structure should to be modelled to account for the total CO content of the disc.
Context. Deuterium fractionation has been used to study the thermal history of prestellar environments. Their formation pathways trace different regions of the disk and may shed light into the ...physical structure of the disk, including locations of important features for planetary formation. Aims. We aim to constrain the radial extent of the main deuterated species; we are particularly interested in spatially characterizing the high and low temperature pathways for enhancing deuteration of these species. Methods. We observed the disk surrounding the Herbig Ae star HD 163296 using ALMA in Band 6 and obtained resolved spectral imaging data of DCO+ (J = 3 − 2), DCN (J = 3 − 2) and N2D+ (J = 3 − 2) with synthesized beam sizes of 0.̋53 × 0.̋42, 0.̋53 × 0.̋42, and 0.̋50 × 0.̋39, respectively. We adopted a physical model of the disk from the literature and use the 3D radiative transfer code LIME to estimate an excitation temperature profile for our detected lines. We modeled the radial emission profiles of DCO+, DCN, and N2D+, assuming their emission is optically thin, using a parametric model of their abundances and our excitation temperature estimates. Results. DCO+ can be described by a three-region model with constant-abundance rings centered at 70 AU, 150 AU, and 260 AU. The DCN radial profile peaks at about 60 AU and N2D+ is seen in a ring at 160 AU. Simple models of both molecules using constant abundances reproduce the data. Assuming reasonable average excitation temperatures for the whole disk, their disk-averaged column densities (and deuterium fractionation ratios) are 1.6–2.6×1012 cm-2 (0.04–0.07), 2.9–5.2×1012 cm-2 (~0.02), and 1.6–2.5×1011 cm-2 (0.34–0.45) for DCO+, DCN, and N2D+, respectively. Conclusions. Our simple best-fit models show a correlation between the radial location of the first two rings in DCO+ and the DCN and N2D+ abundance distributions that can be interpreted as the high and low temperature deuteration pathways regimes. The origin of the third DCO+ ring at 260 AU is unknown but may be due to a local decrease of ultraviolet opacity allowing the photodesorption of CO or due to thermal desorption of CO as a consequence of radial drift and settlement of dust grains. The derived Df values agree with previous estimates of 0.05 for DCO+/HCO+ and 0.02 for DCN/HCN in HD 163296, and 0.3−0.5 for N2D+/N2H+ in AS 209, a T Tauri disk. The high N2D+/N2H+ confirms N2D+ as a good candidate for tracing ionization in the cold outer disk.
Context.
Water is a key volatile that provides insight 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.
Aims.
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.
Methods.
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
up
∕
k
B
= 203.7 K) while ALMA targeted the 4
1,4
–
3
2,1
line (
E
up
∕
k
B
= 322.0 K). Water column densities were derived considering optically thin and thermalized emission. Our observations were sensitive to the emission from the known Keplerian disks around three out of the five Class I objects in the sample.
Results.
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
H
2
18
O
< 10
15
cm
−2
on 100 au scales, which include both the disk and envelope. The upper limits imply a disk-averaged water abundance of ≲10
−6
with respect to H
2
for Class I objects. After taking the physical structure of the disk into account, the upper limit to the water abundance averaged over the inner warm disk with
T
> 100 K is between ~10
−7
and 10
−5
.
Conclusions.
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 of magnitude lower than values found in Class 0 disk-like structures.
Context. Disks are observed around pre-main sequence stars, but how and when they form is still heavily debated. While disks around young stellar objects have been identified through thermal dust ...emission, spatially and spectrally resolved molecular line observations are needed to determine their nature. Only a handful of embedded rotationally supported disks have been identified to date. Aims. We identify and characterize rotationally supported disks near the end of the main accretion phase of low-mass protostars by comparing their gas and dust structures. Methods. Subarcsecond observations of dust and gas toward four Class I low-mass young stellar objects in Taurus are presented at significantly higher sensitivity than previous studies. The 13CO and C18O J = 2–1 transitions at 220 GHz were observed with the Plateau de Bure Interferometer at a spatial resolution of ≤0.8″ (56 AU radius at 140 pc) and analyzed using uv-space position velocity diagrams to determine the nature of their observed velocity gradient. Results. Rotationally supported disks (RSDs) are detected around 3 of the 4 Class I sources studied. The derived masses identify them as Stage I objects; i.e., their stellar mass is higher than their envelope and disk masses. The outer radii of the Keplerian disks toward our sample of Class I sources are ≤100 AU. The lack of on-source C18O emission for TMR1 puts an upper limit of 50 AU on its size. Flattened structures at radii >100 AU around these sources are dominated by infalling motion (υ ∝ r-1). A large-scale envelope model is required to estimate the basic parameters of the flattened structure from spatially resolved continuum data. Similarities and differences between the gas and dust disk are discussed. Combined with literature data, the sizes of the RSDs around Class I objects are best described with evolutionary models with an initial rotation of Ω = 10-14 Hz and slow sound speeds. Based on the comparison of gas and dust disk masses, little CO is frozen out within 100 AU in these disks. Conclusions. Rotationally supported disks with radii up to 100 AU are present around Class I embedded objects. Larger surveys of both Class 0 and I objects are needed to determine whether most disks form late or early in the embedded phase.
Context. The gas and dust in circumstellar disks provide the raw materials to form planets. The study of organic molecules and their building blocks in such disks offers insight into the origin of ...the prebiotic environment of terrestrial planets. Aims. We aim to determine the distribution of formaldehyde, H2CO, in the disk around HD 163296 to assess the contribution of gas- and solid-phase formation routes of this simple organic. Methods. Three formaldehyde lines were observed (H2CO 303–202, H2CO 322–221, and H2CO 321–220) in the protoplanetary disk around the Herbig Ae star HD 163296 with ALMA at ~0.5″ (60 AU) spatial resolution. Different parameterizations of the H2CO abundance were compared to the observed visibilities, using either a characteristic temperature, a characteristic radius or a radial power law index to describe the H2CO chemistry. Similar models were applied to ALMA Science Verification data of C18O. In each scenario, χ2 minimization on the visibilities was used to determine the best-fit model in each scenario. Results. H2CO 303–202 was readily detected via imaging, while the weaker H2CO 322–221 and H2CO 321–220 lines required matched filter analysis to detect. H2CO is present throughout most of the gaseous disk, extending out to ~550 AU. An apparent 50 AU inner radius of the H2CO emission is likely caused by an optically thick dust continuum. The H2CO radial intensity profile shows a peak at ~100 AU and a secondary bump at ~300 AU, suggesting increased production in the outer disk. In all modeling scenarios, fits to the H2CO data show an increased abundance in the outer disk. The overall best-fit H2CO model shows a factor of two enhancement beyond a radius of 270 ± 20 AU, with an inner abundance (relative to H2) of 2 − 5 × 10-12. The H2CO emitting region has a lower limit on the kinetic temperature of T> 20 K. The C18O modeling suggests an order of magnitude depletion of C18O in the outer disk and an abundance of 4 − 12 × 10-8 in the inner disk. Conclusions. There is a desorption front seen in the H2CO emission that roughly coincides with the outer edge of the 1.3 millimeter continuum. The increase in H2CO outer disk emission could be a result of hydrogenation of CO ices on dust grains that are then sublimated via thermal desorption or UV photodesorption. Alternatively, there could be more efficient gas-phase production of H2CO beyond ~300 AU if CO is photodisocciated in this region.
Context.
Many population studies have been performed over the past decade with the Atacama Large millimeter/submillimeter Array (ALMA) to understand the bulk properties of protoplanetary disks around ...young stars. The studied populations have mostly consisted of late spectral type (i.e., G, K & M) stars, with relatively few more massive Herbig stars (spectral types B, A & F). With
Gaia
-updated distances, now is a good time to use ALMA archival data for a Herbig disk population study and take an important step forward in our understanding of planet formation.
Aims.
The aim of this work is to determine the masses and sizes of all Herbig dust disks observed with ALMA to date in a volume-limited sample out to 450 pc. These masses and sizes are put in the context of the Lupus and Upper Sco T Tauri disk populations.
Methods.
ALMA Band 6 and Band 7 archival data of 36 Herbig stars are used, making this work 64% complete out to 225 pc, and 38% complete out to 450 pc also including Orion. Using stellar parameters and distances, the dust masses and sizes of the disks are determined via a curve-of-growth method. Survival analysis is used to obtain cumulative distributions of the dust masses and radii.
Results.
Herbig disks have a higher dust mass than the T Tauri disk populations of Lupus and Upper Sco by factors of ~3 and ~7 respectively. In addition, Herbig disks are often larger than the typical T Tauri disk. Although the masses and sizes of Herbig disks extend over a similar range to those of T Tauri disks, the distributions of masses and sizes of Herbig disks are significantly skewed toward higher values. Lastly, group I disks are more massive than group II disks. An insufficient number of group II disks are observed at sufficient angular resolution to determine whether or not they are also small in size compared to group I disks.
Conclusions.
Herbig disks are skewed towards more massive and larger dust disks compared to T Tauri disks. Based on this we speculate that these differences find their origin in an initial disk mass that scales with the stellar mass, and that subsequent disk evolution enlarges the observable differences, especially if (sub)millimeter continuum optical depth plays a role. Moreover, the larger disk masses and sizes of Herbig stars could be linked to the increasing prevalence of giant planets with host star mass.