skip to main content

Title: Self-consistent Ring Model in Protoplanetary Disks: Temperature Dips and Substructure Formation
Abstract Rings and gaps are ubiquitous in protoplanetary disks. Larger dust grains will concentrate in gaseous rings more compactly due to stronger aerodynamic drag. However, the effects of dust concentration on the ring’s thermal structure have not been explored. Using MCRT simulations, we self-consistently construct ring models by iterating the ring’s thermal structure, hydrostatic equilibrium, and dust concentration. We set up rings with two dust populations having different settling and radial concentration due to their different sizes. We find two mechanisms that can lead to temperature dips around the ring. When the disk is optically thick, the temperature drops outside the ring, which is the shadowing effect found in previous studies adopting a single-dust population in the disk. When the disk is optically thin, a second mechanism due to excess cooling of big grains is found. Big grains cool more efficiently, which leads to a moderate temperature dip within the ring where big dust resides. This dip is close to the center of the ring. Such a temperature dip within the ring can lead to particle pileup outside the ring and feedback to the dust distribution and thermal structure. We couple the MCRT calculations with a 1D dust evolution model more » and show that the ring evolves to a different shape and may even separate to several rings. Overall, dust concentration within rings has moderate effects on the disk’s thermal structure, and a self-consistent model is crucial not only for protoplanetary disk observations but also for planetesimal and planet formation studies. « less
; ; ;
Award ID(s):
1753168 1716259
Publication Date:
Journal Name:
The Astrophysical Journal
Page Range or eLocation-ID:
Sponsoring Org:
National Science Foundation
More Like this
  1. Knowledge of the midplane temperature of protoplanetary disks is one of the key ingredients in theories of dust growth and planet formation. However, direct measurement of this quantity is complicated, and often depends on the fitting of complex models to data. In this paper we demonstrate a method to directly measure the midplane gas temperature from an optically thick molecular line if the disk is moderately inclined. The only model assumption that enters is that the line is very optically thick, specifically in the midplane region where we wish to measure the temperature. Freeze-out of the molecule onto dust grains could thwart this. However, in regions that are expected to be warm enough to avoid freeze-out, this method should work. We apply the method to the CO 2–1 line channel maps of the disk around HD 163296. We find that the midplane temperature between 100 and 400 au drops only mildly from 25 K down to 18 K. While we see no direct evidence of the midplane being optically thin due to strong CO depletion by freeze-out, we cannot rule it out either. The fact that the inferred temperatures are close to the expected CO freeze-out temperature could be anmore »indication of this. Incidently, for the disk around HD 163296 we also find dynamic evidence for a rather abrupt outer edge of the disk, suggestive of outside-in photoevaporation or truncation by an unseen companion.« less
  2. Abstract High spatial resolution CO observations of midinclination (≈30°–75°) protoplanetary disks offer an opportunity to study the vertical distribution of CO emission and temperature. The asymmetry of line emission relative to the disk major axis allows for a direct mapping of the emission height above the midplane, and for optically thick, spatially resolved emission in LTE, the intensity is a measure of the local gas temperature. Our analysis of Atacama Large Millimeter/submillimeter Array archival data yields CO emission surfaces, dynamically constrained stellar host masses, and disk atmosphere gas temperatures for the disks around the following: HD 142666, MY Lup, V4046 Sgr, HD 100546, GW Lup, WaOph 6, DoAr 25, Sz 91, CI Tau, and DM Tau. These sources span a wide range in stellar masses (0.50–2.10 M ⊙ ), ages (∼0.3–23 Myr), and CO gas radial emission extents (≈200–1000 au). This sample nearly triples the number of disks with mapped emission surfaces and confirms the wide diversity in line emitting heights ( z / r ≈ 0.1 to ≳0.5) hinted at in previous studies. We compute the radial and vertical CO gas temperature distributions for each disk. A few disks show local temperature dips or enhancements, some of which correspondmore »to dust substructures or the proposed locations of embedded planets. Several emission surfaces also show vertical substructures, which all align with rings and gaps in the millimeter dust. Combining our sample with literature sources, we find that CO line emitting heights weakly decline with stellar mass and gas temperature, which, despite large scatter, is consistent with simple scaling relations. We also observe a correlation between CO emission height and disk size, which is due to the flared structure of disks. Overall, CO emission surfaces trace ≈2–5× gas pressure scale heights (H g ) and could potentially be calibrated as empirical tracers of H g .« less
  3. We present high-resolution millimeter continuum ALMA observations of the disks around the T Tauri stars LkCa 15 and 2MASS J16100501-2132318 (hereafter, J1610). These transition disks host dust-depleted inner regions, which have possibly been carved by massive planets, and they are of prime interest to the study of the imprints of planet-disk interactions. While at moderate angular resolution, they appear as a broad ring surrounding a cavity, the continuum emission resolves into multiple rings at a resolution of ~60 × 40 mas (~7.5 au for LkCa 15, ~6 au for J1610) and ~7 μ Jy beam −1 rms at 1.3 mm. In addition to a broad extended component, LkCa 15 and J1610 host three and two narrow rings, respectively, with two bright rings in LkCa 15 being radially resolved. LkCa 15 possibly hosts another faint ring close to the outer edge of the mm emission. The rings look marginally optically thick, with peak optical depths of ~0.5 (neglecting scattering), in agreement with high angular resolution observations of full disks. We performed hydrodynamical simulations with an embedded, sub-Jovian-mass planet and show that the observed multi-ringed substructure can be qualitatively explained as the outcome of the planet-disk interaction. We note, however, that themore »choice of the disk cooling timescale alone can significantly impact the resulting gas and dust distributions around the planet, leading to different numbers of rings and gaps and different spacings between them. We propose that the massive outer disk regions of transition disks are favorable places for planetesimals, and possibly second-generation planet formation of objects with a lower mass than the planets carving the inner cavity (typically few M Jup ), and that the annular substructures observed in LkCa 15 and J1610 may be indicative of planetary core formation within dust-rich pressure traps. Current observations are compatible with other mechanisms contributing to the origin of the observed substructures, in particular with regard to narrow rings generated (or facilitated) at the edge of the CO and N 2 snowlines.« less
  4. Vortices are one of the most promising mechanisms to locally concentrate millimeter dust grains and allow the formation of planetesimals through gravitational collapse. The outer disk around the binary system HD 142527 is known for its large horseshoe structure with azimuthal contrasts of ~3–5 in the gas surface density and of ~50 in the dust. Using 13 CO and C 18 O J = 3–2 transition lines, we detect kinematic deviations to the Keplerian rotation, which are consistent with the presence of a large vortex around the dust crescent, as well as a few spirals in the outer regions of the disk. Comparisons with a vortex model suggest velocity deviations up to 350 m s −1 after deprojection compared to the background Keplerian rotation, as well as an extension of ±40 au radially and ~200° azimuthally, yielding an azimuthal-to-radial aspect ratio of ~5. Another alternative for explaining the vortex-like signal implies artificial velocity deviations generated by beam smearing in association with variations of the gas velocity due to gas pressure gradients at the inner and outer edges of the circumbinary disk. The two scenarios are currently difficult to differentiate and, for this purpose, would probably require the use of multiplemore »lines at a higher spatial resolution. The beam smearing effect, due to the finite spatial resolution of the observations and gradients in the line emission, should be common in observations of protoplanetary disks and may lead to misinterpretations of the gas velocity, in particular around ring-like structures.« less
  5. Context. T Tauri stars are low-mass young stars whose disks provide the setting for planet formation. Despite this, their structure is poorly understood. We present new infrared interferometric observations of the SU Aurigae circumstellar environment that offer resolution that is three times higher and a better baseline position angle coverage than previous observations. Aims. We aim to investigate the characteristics of the circumstellar material around SU Aur, constrain the disk geometry, composition and inner dust rim structure. Methods. The CHARA array offers unique opportunities for long baseline observations, with baselines up to 331 m. Using the CLIMB three-telescope combiner in the K -band allows us to measure visibilities as well as closure phase. We undertook image reconstruction for model-independent analysis, and fitted geometric models such as Gaussian and ring distributions. Additionally, the fitting of radiative transfer models constrain the physical parameters of the disk. For the first time, a dusty disk wind is introduced to the radiative transfer code TORUS to model protoplanetary disks. Our implementation is motivated by theoretical models of dusty disk winds, where magnetic field lines drive dust above the disk plane close to the sublimation zone. Results. Image reconstruction reveals an inclined disk with slight asymmetrymore »along its minor-axis, likely due to inclination effects obscuring the inner disk rim through absorption of incident star light on the near-side and thermal re-emission and scattering of the far-side. Geometric modelling of a skewed ring finds the inner rim at 0.17 ± 0.02 au with an inclination of 50.9 ± 1.0° and minor axis position angle 60.8 ± 1.2°. Radiative transfer modelling shows a flared disk with an inner radius at 0.18 au which implies a grain size of 0.4 μ m assuming astronomical silicates and a scale height of 15.0 at 100 au. Among the tested radiative transfer models, only the dusty disk wind successfully accounts for the K -band excess by introducing dust above the mid-plane.« less