No Arabic abstract
We present new determinations of disk surface density, independent of an assumed dust opacity, for a sample of 7 bright, diverse protoplanetary disks using measurements of disk dust lines. We develop a robust method for determining the location of dust lines by modeling disk interferometric visibilities at multiple wavelengths. The disks in our sample have newly derived masses that are 9-27% of their host stellar mass, substantially larger than the minimum mass solar nebula. All are stable to gravitational collapse except for one which approaches the limit of Toomre-Q stability. Our mass estimates are 2-15 times larger than estimates from integrated optically thin dust emission. We derive depleted dust-to-gas ratios with typical values of ~$10^{-3}$ in the outer disk. Using coagulation models we derive dust surface density profiles that are consistent with millimeter dust observations. In these models, the disks formed with an initial dust mass that is a factor of ~10 greater than is presently observed. Of the three disks in our sample with resolved CO line emission, the masses of HD 163296, AS 209, and TW Hya are roughly 3, 115, and 40 times more massive than estimates from CO respectively. This range indicates that CO depletion is not uniform across different disks and that dust is a more robust tracer of total disk mass. Our method of determining surface density using dust lines is robust even if particles form as aggregates and is useful even in the presence of dust substructure caused by pressure traps. The low Toomre-Q values observed in this sample indicate that at least some disks do not accrete efficiently.
We present a novel method for determining the surface density of protoplanetary disks through consideration of disk dust lines which indicate the observed disk radial scale at different observational wavelengths. This method relies on the assumption that the processes of particle growth and drift control the radial scale of the disk at late stages of disk evolution such that the lifetime of the disk is equal to both the drift timescale and growth timescale of the maximum particle size at a given dust line. We provide an initial proof of concept of our model through an application to the disk TW Hya and are able to estimate the disk dust-to-gas ratio, CO abundance, and accretion rate in addition to the total disk surface density. We find that our derived surface density profile and dust-to-gas ratio are consistent with the lower limits found through measurements of HD gas. The CO ice line also depends on surface density through grain adsorption rates and drift and we find that our theoretical CO ice line estimates have clear observational analogues. We further apply our model to a large parameter space of theoretical disks and find three observational diagnostics that may be used to test its validity. First we predict that the dust lines of disks other than TW Hya will be consistent with the normalized CO surface density profile shape for those disks. Second, surface density profiles that we derive from disk ice lines should match those derived from disk dust lines. Finally, we predict that disk dust and ice lines will scale oppositely, as a function of surface density, across a large sample of disks.
Dust evolution in protoplanetary disks from small dust grains to pebbles is key to the planet formation process. The gas in protoplanetary disks should influence the vertical distribution of small dust grains ($sim$1 $mu m$) in the disk.Utilizing archival near-infrared polarized light and millimeter observations, we can measure the scale height and the flare parameter $beta$ of the small dust grain scattering surface and $^{12}$CO gas emission surface for three protoplanetary disks IM Lup, HD 163296, and HD 97048 (CU Cha). For two systems, IM Lup and HD 163296, the $^{12}$CO gas and small dust grains at small radii from the star have similar heights but at larger radii ($>$100 au) the dust grain scattering surface height is lower than the $^{12}$CO gas emission surface height. In the case of HD 97048, the small dust grain scattering surface has similar heights to the $^{12}$CO gas emission surface at all radii. We ran a protoplanetary disk radiative transfer model of a generic protoplanetary disk with TORUS and showed that there is no difference between the observed scattering surface and $^{12}$CO emission surface. We also performed analytical modeling of the system and found that gas-to-dust ratios larger than 100 could explain the observed difference in IM Lup and HD 163296. This is the first direct comparison of observations of gas and small dust grain heights distribution in protoplanetary disks. Future observations of gas emission and near-infrared scattered light instruments are needed to look for similar trends in other protoplanetary disks.
Mid-infrared imaging traces the sub-micron and micron sized dust grains in protoplanetary disks and it offers constraints on the geometrical properties of the disks and potential companions, particularly if those companions have circumplanetary disks. We use the VISIR instrument and its upgrade NEAR on the VLT to take new mid-infrared images of five (pre-)transition disks and one circumstellar disk with proposed planets and obtain the deepest resolved mid-infrared observations to date in order to put new constraints on the sizes of the emitting regions of the disks and the presence of possible companions. We derotate and stack the data to find the disk properties. Where available we compare the data to ProDiMo (Protoplanetary Disk Model) radiation thermo-chemical models to achieve a deeper understanding of the underlying physical processes within the disks. We apply the circularised PSF subtraction method to find upper limits on the fluxes of possible companions and model companions with circumplanetary disks. We resolve three of the six disks and calculate position angles, inclinations and (upper limits to) sizes of emission regions in the disks, improving upper limits on two of the unresolved disks. In all cases the majority of the mid-IR emission comes from small inner disks or the hot inner rims of outer disks. We refine the existing ProDiMo HD 100546 model SED fit in the mid-IR by increasing the PAH abundance relative to the ISM, adopting coronene as the representative PAH, and increase the outer cavity radius to 22.3 AU. We produce flux estimates for putative planetary-mass companions and circumplanetary disks, ruling out the presence of planetary-mass companions with $L > 0.0028 L_{odot}$ for $a > 180$ AU in the HD 100546 system. Upper limits of 0.5 mJy-30 mJy are obtained at 8 $mu$m-12 $mu$m for potential companions in the different disks.
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 $Delta_{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.
The material that formed the present-day Solar System originated in feeding zones in the inner Solar Nebula located at distances within ~20 AU from the Sun, known as the planet-forming zone. Meteoritic and cometary material contain abundant evidence for the presence of a rich and active chemistry in the planet-forming zone during the gas-rich phase of Solar System formation. It is a natural conjecture that analogs can be found amoung the zoo of protoplanetary disks around nearby young stars. The study of the chemistry and dynamics of planet formation requires: 1) tracers of dense gas at 100-1000 K and 2) imaging capabilities of such tracers with 5-100 (0.5-20 AU) milli-arcsec resolution, corresponding to the planet-forming zone at the distance of the closest star-forming regions. Recognizing that the rich infrared (2-200 micron) molecular spectrum recently discovered to be common in protoplanetary disks represents such a tracer, we present a new general raytracing code, RADLite, that is optimized for producing infrared line spectra and images from axisymmetric structures. RADLite can consistently deal with a wide range of velocity gradients, such as those typical for the inner regions of protoplanetary disks. The code is intended as a backend for chemical and excitation codes, and can rapidly produce spectra of thousands of lines for grids of models for comparison with observations. Such radiative transfer tools will be crucial for constraining both the structure and chemistry of planet-forming regions, including data from current infrared imaging spectrometers and extending to the Atacama Large Millimeter Array and the next generation of Extremely Large Telescopes, the James Webb Space Telescope and beyond.