No Arabic abstract
We study the structure of passively heated disks around T Tauri and Herbig Ae stars, and present a vectorized Monte Carlo dust radiative transfer model of protoplanetary disks. The vectorization provides a speed up factor of 100 when compared to a scalar version of the code. Disks are composed of either fluffy carbon and silicate grains of various sizes or dust of the diffuse ISM. The IR emission and the midplane temperature derived by the MC method differ from models where the radiative transfer is solved in slab geometry of small ring segments. In the MC treatment, dusty halos above the disks are considered. Halos lead to an enhanced IR emission and warmer midplane temperature than do pure disks. Under the assumption of hydrostatic equilibrium we find that the disk in the inner rim puffs up, followed by a shadowed region. The shadow reduces the temperature of the midplane and decreases the height of the extinction layer of the disk. It can be seen as a gap in the disk unless the surface is again exposed to direct stellar radiation. There the disk puffs up a second time, a third time and so forth. Therefore several gaps and ring-like structures are present in the disk surface and appear in emission images. They result from shadows in the disks and are present without the need to postulate the existence of any companion or planet. As compared to Herbig Ae stars, such gaps and ring-like structures are more pronounced in regions of terrestrial planets around T Tauri stars.
Planets form in young circumstellar disks called protoplanetary disks. However, it is still difficult to catch planet formation in-situ. Nevertheless, from recent ALMA/SPHERE data, encouraging evidence of the direct and indirect presence of embedded planets has been identified in disks around young stars: co-moving point sources, gravitational perturbations, rings, cavities, and emission dips or shadows cast on disks. The interpretation of these observations needs a robust physical framework to deduce the complex disk geometry. In particular, protoplanetary disk models usually assume the gas pressure scale-height given by the ratio of the sound speed over the azimuthal velocity $H/r = c_{srm }/v_{rm k}$. By doing so, textit{radiative} pressure fields are often ignored, which could lead to a misinterpretation of the real vertical structure of such disks. We follow the evolution of a gaseous disk with an embedded Jupiter mass planet through hydrodynamical simulations, computing the disk scale-height including radiative pressure, which was derived from a generalization of the stellar atmosphere theory. We focus on the vertical impact of the radiative pressure in the vicinity of circumplanetary disks, where temperatures can reach $gtrsim 1000$ K for an accreting planet, and radiative forces can overcome gravitational forces from the planet. The radiation-pressure effects create a vertical optically thick column of gas and dust at the proto-planet location, casting a shadow in scattered light. This mechanism could explain the peculiar illumination patterns observed in some disks around young stars such as HD 169142 where a moving shadow has been detected, or the extremely high aspect-ratio $H/r sim 0.2$ observed in systems like AB Aur and CT Cha.
We study the origin of tail-like structures recently detected around the disk of SU Aurigae and several FU~Orionis-type stars. Dynamic protostellar disks featuring ejections of gaseous clumps and quiescent protoplanetary disks experiencing a close encounter with an intruder star were modeled using the numerical hydrodynamics code FEOSAD. Both the gas and dust dynamics were taken into account, including dust growth and mutual friction between the gas and dust components. Only plane-of-the-disk encounters were considered. Ejected clumps produce a unique type of tail that is characterized by a bow-shock shape. Such tails originate from the supersonic motion of ejected clumps through the dense envelope that often surrounds young gravitationally unstable protostellar disks. The ejected clumps either sit at the head of the tail-like structure or disperse if their mass is insufficient to withstand the head wind of the envelope. On the other hand, close encounters with quiescent protoplanetary disks produce three types of the tail-like structure; we define these as pre-collisional, post-collisional, and spiral tails. These tails can in principle be distinguished from one another by particular features of the gas and dust flow in and around them. We find that the brown-dwarf-mass intruders do not capture circumintruder disks during the encounter, while the subsolar-mass intruders can acquire appreciable circumintruder disks with elevated dust-to-gas ratios, which can ease their observational detection. However, this is true only for prograde collisions; the retrograde intruders fail to collect appreciable amounts of gas or dust from the disk of the target. The predicted mass of dust in the model tail-like structures is higher than what was inferred for similar structures in SU~Aur, FU~Ori, and Z~CMa, making their observational detection feasible. Abridged.
High resolution ALMA observations of protoplanetary disks have revealed that many, if not all primordial disks consist of ring-like dust structures. The origin of these dust rings remains unclear, but a common explanation is the presence of planetary companions that have cleared gaps along their orbit and trapped the dust at the gap edge. A signature of this scenario is a decrease of gas density inside these gaps. In recent work, Isella et al. 2016 derived drops in gas density consistent with Saturn-mass planets inside the gaps in the HD163296 disk through spatially resolved CO isotopologue observations. However, as CO abundance and temperature depends on a large range of factors, the interpretation of CO emission is non-trivial. We use the physical-chemical code DALI to show that the gas temperature increases inside dust density gaps, implying that any gaps in the gas, if present, would have to be much deeper, consistent with planet masses higher than a Jupiter mass. Furthermore, we show that a model with increased grain growth at certain radii, as expected at a snowline, can reproduce the dust rings in HD163296 equally well without the need for companions. This scenario can explain both younger and older disks with observed gaps, as gaps have been seen in systems as young <1 Myr. While the origin of the rings in HD163296 remains unclear, these modeling results demonstrate that care has to be taken when interpreting CO emission in protoplanetary disk observations.
Recent millimeter and infrared observations have shown that gap and ring-like structures are common in both dust thermal emission and scattered-light of protoplanetary disks. We investigate the impact of the so-called Thermal Wave Instability (TWI) on the millimeter and infrared scattered-light images of disks. We perform 1+1D simulations of the TWI and confirm that the TWI operates when the disk is optically thick enough for stellar light, i.e., small-grain-to-gas mass ratio of $gtrsim0.0001$. The mid-plane temperature varies as the waves propagate and hence gap and ring structures can be seen in both millimeter and infrared emission. The millimeter substructures can be observed even if the disk is fully optically thick since it is induced by the temperature variation, while density-induced substructures would disappear in the optically thick regime. The fractional separation between TWI-induced ring and gap is $Delta r/r sim$ 0.2-0.4 at $sim$ 10-50 au, which is comparable to those found by ALMA. Due to the temperature variation, snow lines of volatile species move radially and multiple snow lines are observed even for a single species. The wave propagation velocity is as fast as $sim$ 0.6 ${rm au~yr^{-1}}$, which can be potentially detected with a multi-epoch observation with a time separation of a few years.
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.