Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userOn the evolution of vortices in massive protoplanetary discs
81
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
It is expected that a pressure bump can be formed at the inner edge of a dead-zone, and where vortices can develop through the Rossby Wave Instability (RWI). It has been suggested that self-gravity can significantly affect the evolution of such vortices. We present the results of 2D hydrodynamical simulations of the evolution of vortices forming at a pressure bump in self-gravitating discs with Toomre parameter in the range $4-30$. We consider isothermal plus non-isothermal disc models that employ either the classical $beta$ prescription or a more realistic treatment for cooling. The main aim is to investigate whether the condensating effect of self-gravity can stabilize vortices in sufficiently massive discs. We confirm that in isothermal disc models with ${cal Q} gtrsim 15$, vortex decay occurs due to the vortex self-gravitational torque. For discs with $3lesssim {cal Q} lesssim 7$, the vortex develops gravitational instabilities within its core and undergoes gravitational collapse, whereas more massive discs give rise to the formation of global eccentric modes. In non-isothermal discs with $beta$ cooling, the vortex maintains a turbulent core prior to undergoing gravitational collapse for $beta lesssim 0.1$, whereas it decays if $beta ge 1$. In models that incorpore both self-gravity and a better treatment for cooling, however, a stable vortex is formed with aspect ratio $chi sim 3-4$. Our results indicate that self-gravity significantly impacts the evolution of vortices forming in protoplanetary discs, although the thermodynamical structure of the vortex is equally important for determining its long-term dynamics.
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.
In order to circumvent the loss of solid material through radial drift towards the central star, the trapping of dust inside persistent vortices in protoplanetary discs has often been suggested as a process that can eventually lead to planetesimal formation. Although a few special cases have been discussed, exhaustive studies of possible quasi-steady configurations available for dust-laden vortices and their stability have yet to be undertaken, thus their viability or otherwise as locations for the gravitational instability to take hold and seed planet formation is unclear. In this paper we generalise and extend the well known Kida solution to obtain a series of steady state solutions with varying vorticity and dust density distributions in their cores, in the limit of perfectly coupled dust and gas. We then present a local stability analysis of these configurations, considering perturbations localised on streamlines. Typical parametric instabilities found have growthrates of $~0.05Omega_P$, where $Omega_P$ is the angular velocity at the centre of the vortex. Models with density excess can exhibit many narrow parametric instability bands while those with a concentrated vorticity source display internal shear which significantly affects their stability. However, the existence of these parametric instabilities may not necessarily prevent the possibility of dust accumulation in vortices.
We perform a population synthesis of protoplanetary discs including infall with a total of $50,000$ simulations using a 1D vertically integrated viscous evolution code, studying a large parameter space in final stellar mass. Initial conditions and infall locations are chosen based on the results from a radiation-hydrodynamic population synthesis of circumstellar discs. We also consider a different infall prescription based on a magnetohydrodynamic (MHD) collapse simulation in order to assess the influence of magnetic fields on disc formation. The duration of the infall phase is chosen to produce a stellar mass distribution in agreement with the observationally determined stellar initial mass function. We find that protoplanetary discs are very massive early in their lives. When averaged over the entire stellar population, the discs have masses of $sim 0.3$ and $0.1,mathrm{M_odot}$ for systems based on hydrodynamic or MHD initial conditions, respectively. In systems with final stellar mass $sim 1,mathrm{M_odot}$, we find disc masses of $sim 0.7,mathrm{M_odot}$ for the `hydro case and $sim 0.2,mathrm{M_odot}$ for the `MHD case at the end of the infall phase. Furthermore, the inferred total disc lifetimes are long, $approx 5-7,mathrm{Myr}$ on average, despite our choice of a high value of $10^{-2}$ for the background viscosity $alpha$-parameter. In addition, fragmentation is common in systems that are simulated using hydrodynamic cloud collapse, with more fragments of larger mass formed in more massive systems. In contrast, if disc formation is limited by magnetic fields, fragmentation is suppressed entirely.
Gravitational coupling between protoplanetary discs and planets embedded in them leads to the emergence of spiral density waves, which evolve into shocks as they propagate through the disc. We explore the performance of a semi-analytical framework for describing the nonlinear evolution of the global planet-driven density waves, focusing on the low planet mass regime (below the so-called thermal mass). We show that this framework accurately captures the (quasi-)self-similar evolution of the wave properties expressed in terms of properly rescaled variables, provided that certain theoretical inputs are calibrated using numerical simulations (an approximate, first principles calculation of the wave evolution based on the inviscid Burgers equation is in qualitative agreement with simulations but overpredicts wave damping at the quantitative level). We provide fitting formulae for such inputs, in particular, the strength and global shape of the planet-driven shock accounting for nonlinear effects. We use this nonlinear framework to theoretically compute vortensity production in the disc by the global spiral shock and numerically verify the accuracy of this calculation. Our results can be used for interpreting observations of spiral features in discs, kinematic signatures of embedded planets in CO line emission (kinks), and for understanding the emergence of planet-driven vortices in protoplanetary discs.
Protoplanetary disc systems observed at radio wavelengths often show excess emission above that expected from a simple extrapolation of thermal dust emission observed at short millimetre wavelengths. Monitoring the emission at radio wavelengths can be used to help disentangle the physical mechanisms responsible for this excess, including free-free emission from a wind or jet, and chromospheric emission associated with stellar activity. We present new results from a radio monitoring survey conducted with Australia Telescope Compact Array over the course of several years with observation intervals spanning days, months and years, where the flux variability of 11 T Tauri stars in the Chamaeleon and Lupus star forming regions was measured at 7 and 15 mm and 3 and 6 cm. Results show that for most sources are variable to some degree at 7 mm, indicating the presence of emission mechanisms other than thermal dust in some sources. Additionally, evidence of grain growth to cm-sized pebbles was found for some sources that also have signs of variable flux at 7 mm. We conclude that multiple processes contributing to the emission are common in T Tauri stars at 7 mm and beyond, and that a detection at a single epoch at radio wavelengths should not be used to determine all processes contributing to the emission.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
