Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userThermal evolution of protoplanetary disks: from $beta$-cooling to decoupled gas and dust temperatures
87
0
0.0
(
0
)
Authors
Eduard I. Vorobyov
Ask ChatGPT about the research
No Arabic abstract
Aims: We explore the long-term evolution of young protoplanetary disks with different approaches to computing the thermal structure determined by various cooling and heating processes in the disk and its surroundings. Methods: Numerical hydrodynamics simulations in the thin-disk limit were complemented with three thermal evolution schemes: a simplified $beta$-cooling approach with and without irradiation, in which the rate of disk cooling is proportional to the local dynamical time, a fiducial model with equal dust and gas temperatures calculated taking viscous heating, irradiation, and radiative cooling into account, and also a more sophisticated approach allowing decoupled dust and gas temperatures. Results: We found that the gas temperature may significantly exceed that of dust in the outer regions of young disks thanks to additional compressional heating caused by the infalling envelope material in the early stages of disk evolution and slow collisional exchange of energy between gas and dust in low-density disk regions. The outer envelope however shows an inverse trend with the gas temperatures dropping below that of dust. The global disk evolution is only weakly sensitive to temperature decoupling. Nevertheless, separate dust and gas temperatures may affect the chemical composition, dust evolution, and disk mass estimates. Constant-$beta$ models without stellar and background irradiation fail to reproduce the disk evolution with more sophisticated thermal schemes because of intrinsically variable nature of the $beta$-parameter. Constant-$beta$ models with irradiation can better match the dynamical and thermal evolution, but the agreement is still incomplete. Conclusions: Models allowing separate dust and gas temperatures are needed when emphasis is placed on the chemical or dust evolution in protoplanetary disks, particularly in sub-solar metallicity environments.
rate research
Read 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.
We present new Atacama Large Millimeter/submillimeter Array (ALMA) observations for three protoplanetary disks in Taurus at 2.9,mm and comparisons with previous 1.3,mm data both at an angular resolution of $sim0.1$ (15,au for the distance of Taurus). In the single-ring disk DS Tau, double-ring disk GO Tau, and multiple-ring disk DL Tau, the same rings are detected at both wavelengths, with radial locations spanning from 50 to 120,au. To quantify the dust emission morphology, the observed visibilities are modeled with a parametric prescription for the radial intensity profile. The disk outer radii, taken as 95% of the total flux encircled in the model intensity profiles, are consistent at both wavelengths for the three disks. Dust evolution models show that dust trapping in local pressure maxima in the outer disk could explain the observed patterns. Dust rings are mostly unresolved. The marginally resolved ring in DS Tau shows a tentatively narrower ring at the longer wavelength, an observational feature expected from efficient dust trapping. The spectral index ($alpha_{rm mm}$) increases outward and exhibits local minima that correspond to the peaks of dust rings, indicative of the changes in grain properties across the disks. The low optical depths ($tausim$0.1--0.2 at 2.9,mm and 0.2--0.4 at 1.3,mm) in the dust rings suggest that grains in the rings may have grown to millimeter sizes. The ubiquitous dust rings in protoplanetary disks modify the overall dynamics and evolution of dust grains, likely paving the way towards the new generation of planet formation.
We perform a comprehensive demographic study of the CO extent relative to dust of the disk population in the Lupus clouds, in order to find indications of dust evolution and possible correlations with other properties. We increase up to 42 the number of disks of the region with measured CO and dust sizes ($R_{mathrm{CO}}$, $R_{mathrm{dust}}$) from observations with the Atacama Large Millimeter/submillimeter Array (ALMA). The sizes are obtained from modeling the ${^{12}}$CO $J = 2-1$ line emission and continuum emission at $sim 0.89$ mm with an empirical function (Nuker profile or Gaussian function). The CO emission is more extended than the dust continuum, with a $R_{68%}^{mathrm{CO}}$/$R_{68%}^{mathrm{dust}}$ median value of 2.5, for the entire population and for a sub-sample with high completeness. 6 disks, around $15%$ of the Lupus disk population have a size ratio above 4. Based on thermo-chemical modeling, this value can only be explained if the disk has undergone grain growth and radial drift. These disks do not have unusual properties in terms of stellar mass ($M_{star}$), disk mass ($M_{mathrm{disk}}$), CO and dust sizes ($R_{mathrm{CO}}$, $R_{mathrm{dust}}$), and mass accretion. We search for correlations between the size ratio and $M_{star}$, $M_{mathrm{disk}}$, $R_{mathrm{CO}}$ and $R_{mathrm{dust}}$: only a weak monotonic anti-correlation with the $R_{mathrm{dust}}$ is found. The lack of strong correlations is remarkable and suggests that the bulk of the population may be in a similar evolutionary stage, independent of the stellar and disk properties. These results should be further investigated, since the optical depth difference between CO and dust continuum may play a role in the inferred size ratios. Lastly, the CO emission for the majority of the disks is consistent with optically thick emission and an average CO temperature of around 30 K.
Recent high angular resolution observations of protoplanetary disks at different wavelengths have revealed several kinds of structures, including multiple bright and dark rings. Embedded planets are the most used explanation for such structures, but there are alternative models capable of shaping the dust in rings as it has been observed. We assume a disk around a Herbig star and investigate the effect that ice lines have on the dust evolution, following the growth, fragmentation, and dynamics of multiple dust size particles, covering from 1 $mu$m to 2 m sized objects. We use simplified prescriptions of the fragmentation velocity threshold, which is assumed to change radially at the location of one, two, or three ice lines. We assume changes at the radial location of main volatiles, specifically H$_2$O, CO$_2$, and NH$_3$. Radiative transfer calculations are done using the resulting dust density distributions in order to compare with current multiwavelength observations. We find that the structures in the dust density profiles and radial intensities at different wavelengths strongly depend on the disk viscosity. A clear gap of emission can be formed between ice lines and be surrounded by ring-like structures, in particular between the H$_2$O and CO$_2$ (or CO). The gaps are expected to be shallower and narrower at millimeter emission than at near-infrared, opposite to model predictions of particle trapping. In our models, the total gas surface density is not expected to show strong variations, in contrast to other gap-forming scenarios such as embedded giant planets or radial variations of the disk viscosity.
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. We analyze spatially resolved observations of 36 protoplanetary disks in the Lupus star forming complex from our ALMA survey at 890 micron, 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. We fit the observations directly in the uv-plane using a two-layer disk model that computes the 890 micron emission by solving the energy balance at each disk radius. 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 Solar masses, and we find no trend 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/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/disk size relation with viscous spreading, however with the current measurements other mechanisms cannot be ruled out.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
