No Arabic abstract
The variety of observed protoplanetary disks in polarimetric light motivates a taxonomical study to constrain their evolution and establish the current framework of this type of observations. We classified 58 disks with available polarimetric observations into six major categories (Ring, Spiral, Giant, Rim, Faint, and Small disks) based on their appearance in scattered light. We re-calculated the stellar and disk properties from the newly available GAIA DR2 and related these properties with the disk categories. More than a half of our sample shows disk sub-structures. For the remaining sources, the absence of detected features is due to their faintness, to their small size, or to the disk geometry. Faint disks are typically found around young stars and typically host no cavity. There is a possible dichotomy in the near-IR excess of sources with spiral-disks (high) and ring-disks (low). Like spirals, shadows are associated with a high near-IR excess. If we account for the pre-main sequence evolutionary timescale of stars with different mass, spiral arms are likely associated to old disks. We also found a loose, shallow declining trend for the disk dust mass with time. Protoplanetary disks may form sub-structures like rings very early in their evolution but their detectability in scattered light is limited to relatively old sources (more than 5 Myr) where the recurrently detected disk cavities allow to illuminate the outer disk. The shallow decrease of disk mass with time might be due to a selection effect, where disks observed thus far in scattered light are typically massive, bright transition disks with longer lifetime than most disks. Our study points toward spirals and shadows being generated by planets of fraction-to-few Jupiter masses that leave their (observed) imprint on both the inner disk near the star and the outer disk cavity.
Context: Vertically hydrostatic protoplanetary disk models are based on the assumption that the main heating source, stellar irradiation, does not vary much with time. However, it is known that accreting young stars are variable sources of radiation. This is particularly evident for outbursting sources such as EX Lupi and FU Orionis stars. Aim: We investigate how such outbursts affect the vertical structure of the outer regions of the protoplanetary disk, in particular their appearance in scattered light at optical and near-infrared wavelengths. Methods: We employ the 3D FARGOCA radiation-hydrodynamics code, in polar coordinates, to compute the time-dependent behavior of the axisymmetric disk structure. The outbursting inner disk region is not included explicitly. Instead, its luminosity is added to the stellar luminosity and is thus included in the irradiation of the outer disk regions. For time snapshots of interest we insert the density structure into the RADMC-3D radiative transfer code and compute the appearance of the disk at optical/near-infrared wavelengths. Results: We find that, depending on the amplitude of the outbursts, the vertical structure of the disk can become highly dynamic, featuring circular surface waves of considerable amplitude. These hills and valleys on the disks surface show up in the scattered light images as bright and dark concentric rings. Initially these rings are small and act as standing waves, but they subsequently lead to outward propagating waves, like the waves produced by a stone thrown into a pond. These waves continue long after the actual outburst has died out. Conclusions: We propose that some of the multi-ringed structures seen in optical/infrared images of several protoplanetary disks may have their origin in outbursts that occurred decades or centuries ago.
Three-dimensional hydrodynamic numerical simulations have demonstrated that the structure of a protoplanetary disc may be strongly affected by a planet orbiting in a plane that is misaligned to the disc. When the planet is able to open a gap, the disc is separated into an inner, precessing disc and an outer disc with a warp. In this work, we compute infrared scattered light images to investigate the observational consequences of such an arrangement. We find that an inner disc misaligned by a less than a degree to the outer disc is indeed able to cast a shadow at larger radii. In our simulations a planet of around 6 Jupiter masses inclined by around 2 degrees is enough to warp the disc and cast a shadow with a depth of more than 10% of the average flux at that radius. We also demonstrate that warp in the outer disc can cause a variation in the azimuthal brightness profile at large radii. Importantly, this latter effect is a function of the distance from the star and is most prominent in the outer disc. We apply our model to the TW Hya system, where a misaligned, precessing inner disc has been invoked to explain an recently observed shadow in the outer disc. Consideration of the observational constraints suggest that an inner disc precessing due to a misaligned planet is an unlikely explanation for the features found in TW Hya.
High-contrast scattered light observations have revealed the surface morphology of several dozens of protoplanetary disks at optical and near-infrared wavelengths. Inclined disks offer the opportunity to measure part of the phase function of the dust grains that reside in the disk surface which is essential for our understanding of protoplanetary dust properties and the early stages of planet formation. We aim to construct a method which takes into account how the flaring shape of the scattering surface of an (optically thick) protoplanetary disk projects onto the image plane of the observer. This allows us to map physical quantities (scattering radius and scattering angle) onto scattered light images and retrieve stellar irradiation corrected (r^2-scaled) images and dust phase functions. We apply the method on archival polarized intensity images of the protoplanetary disk around HD 100546 that were obtained with VLT/SPHERE in R-band and VLT/NACO in H- and Ks-band. The brightest side of the r^2-scaled R-band polarized intensity image of HD 100546 changes from the far to the near side of the disk when a flaring instead of a geometrically flat disk surface is used for the r^2-scaling. The decrease in polarized surface brightness in the scattering angle range of ~40-70 deg is likely a result of the dust phase function and degree of polarization which peak in different scattering angle regimes. The derived phase functions show part of a forward scattering peak which indicates that large, aggregate dust grains dominate the scattering opacity in the disk surface. Projection effects of a protoplanetary disk surface need to be taken into account to correctly interpret scattered light images. Applying the correct scaling for the correction of stellar irradiation is crucial for the interpretation of the images and the derivation of the dust properties in the disk surface layer.
Spatially resolved observations of protoplanetary discs are revealing that their inner regions can be warped or broken from the outer disc. A few mechanisms are known to lead to such 3D structures; among them, the interaction with a stellar companion. We perform a 3D SPH simulation of a circumbinary disc misaligned by $60^circ$ with respect to the binary orbital plane. The inner disc breaks from the outer regions, precessing as a rigid body, and leading to a complex evolution. As the inner disc precesses, the misalignment angle between the inner and outer discs varies by more than $100^circ$. Different snapshots of the evolution are post-processed with a radiative transfer code, in order to produce observational diagnostics of the process. Even though the simulation was produced for the specific case of a circumbinary disc, most of the observational predictions hold for any disc hosting a precessing inner rim. Synthetic scattered light observations show strong azimuthal asymmetries, where the pattern depends strongly on the misalignment angle between inner and outer disc. The asymmetric illumination of the outer disc leads to azimuthal variations of the temperature structure, in particular in the upper layers, where the cooling time is short. These variations are reflected in asymmetric surface brightness maps of optically thick lines, as CO $J$=3-2. The kinematical information obtained from the gas lines is unique in determining the disc structure. The combination of scattered light images and (sub-)mm lines can distinguish between radial inflow and misaligned inner disc scenarios.
We perform calculations of our one-dimensional, two-zone disk model to study the long-term evolution of the circumstellar disk. In particular, we adopt published photoevaporation prescriptions and examine whether the photoevaporative loss alone, coupled with a range of initial angular momenta of the protostellar cloud, can explain the observed decline of the frequency of optically-thick dusty disks with increasing age. In the parameter space we explore, disks have accreting and/or non-accreting transitional phases lasting of $lesssim20 %$ of their lifetime, which is in reasonable agreement with observed statistics. Assuming that photoevaporation controls disk clearing, we find that initial angular momentum distribution of clouds needs to be weighted in favor of slowly rotating protostellar cloud cores. Again, assuming inner disk dispersal by photoevaporation, we conjecture that this skewed angular momentum distribution is a result of fragmentation into binary or multiple stellar systems in rapidly-rotating cores. Accreting and non-accreting transitional disks show different evolutionary paths on the $dot{M}-R_{rm wall}$ plane, which possibly explains the different observed properties between the two populations. However, we further find that scaling the photoevaporation rates downward by a factor of 10 makes it difficult to clear the disks on the observed timescales, showing that the precise value of the photoevaporative loss is crucial to setting the clearing times. While our results apply only to pure photoevaporative loss (plus disk accretion), there may be implications for models in which planets clear disks preferentially at radii of order 10 AU.