No Arabic abstract
We investigate the evolution of grains composed of an ice shell surrounding an olivine core as they pass through a spiral shock in a protoplanetary disk. We use published three-dimensional radiation-hydrodynamics simulations of massive self-gravitating protoplanetary disks to extract the thermodynamics of spiral shocks in the region between 10 and 20 AU from the central star. As the density wave passes, it heats the grains, causing them to lose their ice shell and resulting in a lowering of the grain opacity. In addition, since grains of different sizes will have slightly different temperatures, there will be a migration of ice from the hotter grains to the cooler ones. The rate of migration depends on the temperature of the background gas, so a grain distribution that is effectively stable for low temperatures, can undergo an irreversible change in opacity if the gas is temporarily heated to above $sim 150$,K. We find that the opacity can drop more, and at a significantly faster rate throughout the spiral shocks relative to the prediction of standard dust grains models adopted in hydrodynamical calculations of protoplanetary disks. This would lead to faster gas cooling within spiral arms. We discuss the implications of our results on the susceptibility of disk fragmentation into sub-stellar objects at distances of a few tens of astronomical units.
We investigate how the detectability of signatures of self-gravity in a protoplanetary disc depends on its temporal evolution. We run a one-dimensional model for secular timescales to follow the disc mass as a function of time. We then combine this with three-dimensional global hydrodynamics simulations that employ a hybrid radiative transfer method to approximate realistic heating and cooling. We simulate ALMA continuum observations of these systems, and find that structures induced by the gravitational instability (GI) are readily detectable when $q=M_mathrm{disc}/M_*gtrsim 0.25$ and $R_mathrm{outer}lesssim 100$ au. The high accretion rate generated by gravito-turbulence in such a massive disc drains its mass to below the detection threshold in $sim10^4$ years, or approximately 1 % of the typical disc lifetime. Therefore, discs with spiral arms detected in ALMA dust observations, if generated by self-gravity, must either be still receiving infall to maintain a high $q$ value, or have just emerged from their natal envelope. Detection of substructure in systems with lower $q$ is possible, but would require a specialist integration with the most extended configuration over several days. This disfavours the possibility of GI-caused spiral structure in systems with $q<0.25$ being detected in relatively short integration times, such as those found in the DSHARP ALMA survey (Andrews et al. 2018; Huang et al. 2018). We find no temporal dependence of detectability on dynamical timescales
In this dissertation, I describe theoretical and numerical studies that address the three-dimensional behavior of spiral shocks in protoplanetary disks and the controversial topic of gas giant formation by disk instability. For this work, I discuss characteristics of gravitational instabilities (GIs) in bursting and asymptotic phase disks; outline a theory for the three-dimensional structure of spiral shocks, called shock bores, for isothermal and adiabatic gases; consider convection as a source of cooling for protoplanetary disks; investigate the effects of opacity on disk cooling; use multiple analyses to test for disk stability against fragmentation; test the sensitivity of GI behavior to radiation boundary conditions; measure shock strengths and frequencies in GI-bursting disks; evaluate temperature fluctuations in unstable disks; and investigate whether spiral shocks can form chondrules when GIs activate. The numerical methods developed for these studies are discussed, including a radiation transport routine that explicitly couples the low and high optical depth regimes and a routine that models ortho and parahydrogen. Finally, I explore the hypothesis that chondrule formation and the FU Ori phenomenon are driven by GI activation in dead zones.
Recent observations of protoplanetary disks, as well as simulations of planet-disk interaction, have suggested that a single planet may excite multiple spiral arms in the disk, in contrast to the previous expectations based on linear theory (predicting a one-armed density wave). We re-assess the origin of multiple arms in the framework of linear theory, by solving for the global two-dimensional response of a non-barotropic disk to an orbiting planet. We show that the formation of a secondary arm in the inner disk, at about half of the orbital radius of the planet, is a robust prediction of linear theory. This arm becomes stronger than the primary spiral at several tenths of the orbital radius of the planet. Several additional, weaker spiral arms may also form in the inner disk. On the contrary, a secondary spiral arm is unlikely to form in the outer disk. Our linear calculations, fully accounting for the global behavior of both the phases and amplitudes of perturbations, generally support the recently proposed WKB phase argument for the secondary arm origin (as caused by the intricacy of constructive interference of azimuthal harmonics of the perturbation at different radii). We provide analytical arguments showing that the process of a single spiral wake splitting up into multiple arms is a generic linear outcome of wave propagation in differentially rotating disks. It is not unique to planet-driven waves and occurs also in linear calculations of spiral wakes freely propagating with no external torques. These results are relevant for understanding formation of multiple rings and gaps in protoplanetary disks.
The long-term evolution of a circumstellar disk starting from its formation and ending in the T Tauri phase was simulated numerically with the purpose of studying the evolution of dust in the disk with distinct values of viscous alpha-parameter and dust fragmentation velocity v_frag. We solved numerical hydrodynamics equations in the thin-disk limit, which are modified to include a dust component consisting of two parts: sub-micron-sized dust and grown dust with a maximum radius a_r. The former is strictly coupled to the gas, while the latter interacts with the gas via friction. The conversion of small to grown dust, dust growth, and dust self-gravity are also considered. We found that the process of dust growth known for the older protoplanetary phase also holds for the embedded phase of disk evolution. The dust growth efficiency depends on the radial distance from the star - a_r is largest in the inner disk and gradually declines with radial distance. In the inner disk, a_r is limited by the dust fragmentation barrier. The process of small-to-grown dust conversion is very fast once the disk is formed. The total mass of grown dust in the disk (beyond 1 AU) reaches tens or even hundreds of Earth masses already in the embedded phase of star formation and even a greater amount of grown dust drifts in the inner, unresolved 1 AU of the disk. Dust does not usually grow to radii greater than a few cm. A notable exception are models with alpha <= 10^{-3}, in which case a zone with reduced mass transport develops in the inner disk and dust can grow to meter-sized boulders in the inner 10 AU. Grown dust drifts inward and accumulates in the inner disk regions. This effect is most pronounced in the alpha <= 10^{-3} models where several hundreds of Earth masses can be accumulated in a narrow region of several AU from the star by the end of embedded phase. (abridged).
Recent ALMA surveys of protoplanetary disks have shown that for most disks the extent of the gas emission is greater than the extent of the thermal emission of the millimeter-sized dust. Both line optical depth and the combined effect of radially dependent grain growth and radial drift may contribute to this observed effect. For a sample of 10 disks from the Lupus survey we investigate how well dust-based models without radial dust evolution reproduce the observed 12CO outer radius, and determine whether radial dust evolution is required to match the observed gas-dust size difference. We used the thermochemical code DALI to obtain 12CO synthetic emission maps and measure gas and dust outer radii (Rco, Rmm) using the same methods as applied to the observations, which were compared to observations on a source-by-source basis. For 5 disks we find that the observed gas-dust size difference is larger than the gas-dust size difference due to optical depth, indicating that we need both dust evolution and optical depth effects to explain the observed gas-dust size difference. For the other 5 disks the observed gas-dust size difference can be explained using only line optical depth effects. We also identify 6 disks not included in our initial sample but part of a survey of the same star-forming region that show significant 12CO emission beyond 4 x Rmm. These disks, for which no Rco is available, likely have gas-dust size differences greater than 4 and are difficult to explain without substantial dust evolution. Our results suggest that radial drift and grain growth are common features among both bright and fain disks. The effects of radial drift and grain growth can be observed in disks where the dust and gas radii are significantly different, while more detailed models and deeper observations are needed to see this effect in disks with smaller differences.