ABSTRACT We present the first high-resolution sub-millimeter survey of both dust and gas for a large population of protoplanetary disks. Characterizing fundamental properties of protoplanetary disks ...on a statistical level is critical to understanding how disks evolve into the diverse exoplanet population. We use the Atacama Large Millimeter/Submillimeter Array (ALMA) to survey 89 protoplanetary disks around stars with in the young (1-3 Myr), nearby (150-200 pc) Lupus complex. Our observations cover the 890 m continuum and the 13CO and C18O 3-2 lines. We use the sub-millimeter continuum to constrain to a few Martian masses (0.2-0.4 M⊕) and the CO isotopologue lines to constrain to roughly a Jupiter mass (assuming an interstellar medium (ISM)-like abundance). Of 89 sources, we detect 62 in continuum, 36 in 13CO, and 11 in C18O at significance. Stacking individually undetected sources limits their average dust mass to Lunar masses (0.03 M⊕), indicating rapid evolution once disk clearing begins. We find a positive correlation between and M*, and present the first evidence for a positive correlation between and M*, which may explain the dependence of giant planet frequency on host star mass. The mean dust mass in Lupus is 3× higher than in Upper Sco, while the dust mass distributions in Lupus and Taurus are statistically indistinguishable. Most detected disks have and gas-to-dust ratios , assuming an ISM-like abundance; unless CO is very depleted, the inferred gas depletion indicates that planet formation is well underway by a few Myr and may explain the unexpected prevalence of super-Earths in the exoplanet population.
Context. The formation of planets strongly depends on the total amount as well as on the spatial distribution of solids in protoplanetary disks. Thanks to the improvements in resolution and ...sensitivity provided by ALMA, measurements of the surface density of mm-sized grains are now possible on large samples of disks. Such measurements provide statistical constraints that can be used to inform our understanding of the initial conditions of planet formation. Aims. We aim to analyze spatially resolved observations of 36 protoplanetary disks in the Lupus star forming complex from our ALMA survey at 890 μm, aiming to determine physical properties such as the dust surface density, the disk mass and size, and to provide a constraint on the temperature profile. Methods. We fit the observations directly in the uv-plane using a two-layer disk model that computes the 890 μm emission by solving the energy balance at each disk radius. Results. For 22 out of 36 protoplanetary disks we derive robust estimates of their physical properties. The sample covers stellar masses between ~0.1 and ~ 2 M⊙, and we find no trend in the relationship between the average disk temperatures and the stellar parameters. We find, instead, a correlation between the integrated sub-mm flux (a proxy for the disk mass) and the exponential cut-off radii (a proxy of the disk size) of the Lupus disks. Comparing these results with observations at similar angular resolution of Taurus-Auriga and Ophiuchus disks found in literature and scaling them to the same distance, we observe that the Lupus disks are generally fainter and larger at a high level of statistical significance. Considering the 1–2 Myr age difference between these regions, it is possible to tentatively explain the offset in the disk mass-size relation with viscous spreading, however with the current measurements other mechanisms cannot be ruled out.
Aims. We attempt to explain grain growth to mm sized particles and their retention in the outer regions of protoplanetary disks, as observed at sub-mm and mm wavelengths, by investigating whether ...strong inhomogeneities in the gas density profiles can decelerate excessive radial drift and help the dust particles to grow. Methods. We use coagulation/fragmentation and disk-structure models, to simulate the evolution of dust in a bumpy surface density profile, which we mimic with a sinusoidal disturbance. For different values of the amplitude and length scale of the bumps, we investigate the ability of this model to produce and retain large particles on million-year timescales. In addition, we compare the pressure inhomogeneities considered in this work with the pressure profiles that come from magnetorotational instability. Using the Common Astronomy Software Applications ALMA simulator, we study whether there are observational signatures of these pressure inhomogeneities that can be seen with ALMA. Results. We present the conditions required to trap dust particles and the corresponding calculations predicting the spectral slope in the mm-wavelength range, to compare with current observations. Finally, we present simulated images using different antenna configurations of ALMA at different frequencies, to show that the ring structures will be detectable at the distances of either the Taurus Auriga or Ophiucus star-forming regions.
Context. Transitional disks are thought to be a late evolutionary stage of protoplanetary disks whose inner regions have been depleted of dust. The mechanism responsible for this depletion is still ...under debate. To constrain the various models it is mandatory to have a good understanding of the properties of the gas content in the inner part of the disk. Aims. Using X-Shooter broad band − UV to near-infrared − medium-resolution spectroscopy, we derive the stellar, accretion, and wind properties of a sample of 22 transitional disks. The analysis of these properties allows us to place strong constraints on the gas content in a region very close to the star (≲0.2 AU) that is not accessible with any other observational technique. Methods. We fitted the spectra with a self-consistent procedure to simultaneously derive spectral type, extinction, and accretion properties of the targets. From the continuum excess at near-infrared wavelength we distinguished whether our targets have dust free inner holes. By analyzing forbidden emission lines, we derived the wind properties of the targets. We then compared our findings with results for classical T Tauri stars. Results. The accretion rates and wind properties of 80% of the transitional disks in our sample, which is strongly biased toward stongly accreting objects, are comparable to those of classical T Tauri stars. Thus, there are (at least) some transitional disks with accretion properties compatible with those of classical T Tauri stars, irrespective of the size of the dust inner hole. Only in two cases are the mass accretion rates much lower, while the wind properties remain similar. We detected no strong trend of the mass accretion rates with the size of the dust-depleted cavity or with the presence of a dusty optically thick disk very close to the star. These results suggest that, close to the central star, there is a gas-rich inner disk with a density similar to that of classical T Tauri star disks. Conclusions. The sample analyzed here suggests that, at least for some objects, the process responsible of the inner disk clearing allows for a transfer of gas from the outer disk to the inner region. This should proceed at a rate that does not depend on the physical mechanisms that produces the gap seen in the dust emission and results in a gas density in the inner disk similar to that of unperturbed disks around stars of similar mass.
We present new 3 mm ATCA data of two Class I young stellar objects (YSOs) in the Ophiucus star forming region: Elias29 and WL12. For our analysis we compare them with archival 1.1 mm SMA data. In the ...(u,v) plane the two sources present a similar behavior: a nearly constant non-zero emission at long baselines, which suggests the presence of an unresolved component and an increase of the fluxes at short baselines, related to the presence of an extended envelope. Our data analysis leads to unusually low values of the spectral index α1.1−3 mm, which may indicate that mm-sized dust grains have already formed both in the envelopes and in the disk-like structures at such early stages. To explore the possible scenarios for the interpretation of the sources we perform a radiative transfer modeling using a Monte Carlo code, in order to take into account possible deviations from the Rayleigh-Jeans and optically thin regimes. Comparison between the model outputs and the observations indicates that dust grains may form aggregates up to millimeter size already in the inner regions of the envelopes of Class I YSOs. Moreover, we conclude that the embedded disk-like structures in our two Class I YSOs are probably very compact, in particular in the case of WL12, with outer radii down to tens of AU.
The mass accretion rate, M sub(acc), is a key quantity for the understanding of the physical processes governing the evolution of accretion discs around young low-mass (Mlow *< or =2.0 M sub(middot ...in circle)) stars and substellar objects (YSOs). We present here the results of a study of the stellar and accretion properties of the (almost) complete sample of class II and transitional YSOs in the Lupus I, II, III and IV clouds, based on spectroscopic data acquired with the VLT/X-shooter spectrograph. Our study combines the dataset from our previous work with new observations of 55 additional objects. We have investigated 92 YSO candidates in total, 11 of which have been definitely identified with giant stars unrelated to Lupus. The stellar and accretion properties of the 81 bona fide YSOs, which represent more than 90% of the whole class II and transition disc YSO population in the aforementioned Lupus clouds, have been homogeneously and self-consistently derived, allowing for an unbiased study of accretion and its relationship with stellar parameters. The accretion luminosity, L sub(acc), increases with the stellar luminosity, Llow *, with an overall slope of ~1.6, similar but with a smaller scatter than in previous studies. There is a significant lack of strong accretors below Llow *approximate 0.1L sub(middot in circle), where L sub(acc) is always lower than 0.01Llow *. We argue that the L sub(acc)-Llow * slope is not due to observational biases, but is a true property of the Lupus YSOs. The logM sub(acc)- logMlow * correlation shows a statistically significant evidence of a break, with a steeper relation for Mlow *< or =0.2M sub(middot in circle) and a flatter slope for higher masses. The bimodality of the M sub(acc)- Mlow * relation is confirmed with four different evolutionary models used to derive the stellar mass. The bimodal behaviour of the observed relationship supports the importance of modelling self-gravity in the early evolution of the more massive discs, but other processes, such as photo-evaporation and planet formation during the YSO's lifetime, may also lead to disc dispersal on different timescales depending on the stellar mass. The sample studied here more than doubles the number of YSOs with homogeneously and simultaneously determined L sub(acc) and luminosity, L sub(line), of many permitted emission lines. Hence, we also refined the empirical relationships between L sub(acc) and L sub(line) on a more solid statistical basis.
Context. Studies of the formation and evolution of young stars and their disks rely on knowledge of the stellar parameters of the young stars. The derivation of these parameters is commonly based on ...comparison with photospheric template spectra. Furthermore, chromospheric emission in young active stars impacts the measurement of mass accretion rates, a key quantity for studying disk evolution. Aims. Here we derive stellar properties of low-mass (M⋆≲ 2 M⊙) pre-main sequence stars without disks, which represent ideal photospheric templates for studies of young stars. We also use these spectra to constrain the impact of chromospheric emission on the measurements of mass accretion rates. The spectra are reduced, flux-calibrated, and corrected for telluric absorption, and are made available to the community. Methods. We derive the spectral type for our targets by analyzing the photospheric molecular features present in their VLT/X-shooter spectra by means of spectral indices and comparison of the relative strength of photospheric absorption features. We also measure effective temperature, gravity, projected rotational velocity, and radial velocity from our spectra by fitting them with synthetic spectra with the ROTFIT tool. The targets have negligible extinction (AV< 0.5 mag) and spectral type from G5 to K6, and from M6.5 to M8. They thus complement the library of photospheric templates presented in our previous publication. We perform synthetic photometry on the spectra to derive the typical colors of young stars in different filters. We measure the luminosity of the emission lines present in the spectra and estimate the noise due to chromospheric emission in the measurements of accretion luminosity in accreting stars. Results. We provide a calibration of the photospheric colors of young pre-main sequence stars as a function of their spectral type in a set of standard broad-band optical and near-infrared filters. The logarithm of the noise on the accretion luminosity normalized to the stellar luminosity is roughly constant and equal to ~−2.3 for targets with masses larger than 1 solar mass, and decreases with decreasing temperatures for lower-mass stars. For stars with masses of ~ 1.5 M⊙ and ages of ~ 1−5 Myr, the chromospheric noise converts to a limit of measurable mass accretion rates of ~ 3 × 10-10M⊙/yr. The limit on the mass accretion rate set by the chromospheric noise is of the order of the lowest measured values of mass accretion rates in Class II objects.
We present the most sensitive 3 mm-survey to date of protoplanetary disks carried out in the Taurus-Auriga star forming region (average rms of about 0.3 mJy), using the IRAM PdBI. With our high ...detection rate of 17/19 we provide the first detections at wavelengths longer than about 1 mm for 12 sources. This enables us to statistically study the mm SED slopes and dust properties of faint disks and compare them to brighter disks using a uniform analysis method. With these new data and literature measurements at sub-millimeter and millimeter wavelengths, we analyze the dust properties of a sample of 21 isolated disks around T Tauri stars in the Taurus-Auriga star forming region. Together with the information about the disks' spatial extension from sub/mm-mm interferometric studies, we derive from the observed sub-mm/mm spectral energy distribution constraints on the dust opacity law at these wavelengths, using two-layer flared disk models and a self-consistent dust model that properly takes into account the variation of the dust opacity with grain growth. We find evidence for the presence of dust particles in the disk midplane that have grown to sizes as large as at least 1 millimeter in all the disks of our sample, confirming what was previously observed on smaller brighter objects. This indicates that the dust coagulation from ISM dust to mm-sized grains is a very fast process in protoplanetary disks, which appears to occur before a young stellar object enters the class II evolutionary stage. The amount of these large grains in the disk outer regions is also found to be stationary throughout the whole class II evolutionary stage, which indicates that mechanisms slowing down the dust inward migration play an important role in the Taurus-Auriga protoplanetary disks. Another result is that the spectral index between 1 and 3 mm for the six faintest disks in our sample is on average smaller than for the brighter disks, indicating either that these fainter yet unmapped disks are spatially much less extended than the brighter spatially resolved disks, or that fainter disks have typically larger dust grains in their outer regions. Considering that these fainter disks are more representative of the bulk of the disk population than the brighter ones, this may have important consequences for the theories of planetesimal formation and disk formation and evolution. Finally, we investigate the relations between the derived dust properties, namely dust mass and grain growth, and the properties of the central star, like its mass, age and mass accretion rate.
This paper discusses the properties of the inner puffed-up rim that forms in circumstellar disks when dust evaporates. We argue that the rim shape is controlled by a fundamental property of ...circumstellar disks, namely their very large vertical density gradient, through the dependence of grain evaporation temperature on gas density. As a result, the bright side of the rim is curved, rather than vertical, as expected when a constant evaporation temperature is assumed. We have computed a number of rim models that take into account this effect in a self-consistent way. The results show that the curved rim (as the vertical rim) emits most of its radiation in the near and mid-IR, and provides a simple explanation for the observed values of the near-IR excess (the “3 μm bump” of Herbig Ae stars). Contrary to the vertical rim, for curved rims the near-IR excess does not depend much on the inclination, being maximum for face-on objects. We then computed synthetic images of the curved rim seen under different inclinations; face-on rims are seen as bright, centrally symmetric rings on the sky; increasing the inclination, the rim takes an elliptical shape, with one side brighter than the other.
ABSTRACT
We analyse spatially resolved ALMA observations at 0.9, 1.3, and 3.1 mm for the 26 brightest protoplanetary discs in the Lupus star-forming region. We characterize the discs multiwavelength ...brightness profiles by fitting the interferometric visibilities in a homogeneous way, obtaining effective disc sizes at the three wavelengths, spectral index profiles, and optical depth estimates. We report three fundamental discoveries: first, the millimetre continuum size–luminosity relation already observed at 0.9 mm is also present at 1.3 mm with an identical slope, and at 3.1 mm with a steeper slope, confirming that emission at longer wavelengths becomes increasingly optically thin. Second, when observed at 3.1 mm the discs appear to be only 9 per cent smaller than when observed at 0.9 mm, in tension with models of dust evolution that predict a starker difference. Third, by forward modelling the sample of measurements with a simple parametric disc model, we find that the presence of large grains ($a_\mathrm{max}\gt 1\,$ mm) throughout the discs is the most favoured explanation for all discs as it reproduces simultaneously their spectral indices, optical depth, luminosity, and radial extent in the 0.9–1.3 mm wavelength range. We also find that the observations can be alternatively interpreted with the discs being dominated by optically thick, unresolved, substructures made of mm-sized grains with a high scattering albedo.