Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userPlanetesimal Formation at the Boundary Between Steady Super/Sub-Keplerian Flow Created by Inhomogeneous Growth of Magnetorotational Instability
100
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
We have studied formation of planetesimals at a radial pressure bump in a protoplanetary disk created by radially inhomogeneous magnetorotational instability (MRI), through three-dimensional resistive MHD simulations including dust particles. In our previous papers, we showed that the inhomogeneous MRI developing in non-uniform structure of magnetic field or magnetic resistivity can transform the local gas flow in the disk to a quasi-steady state with local rigid rotation that is no more unstable against the MRI. Since the outer part of the rigid rotation is super-Keplerian flow, a quasi-static pressure bump is created and dust concentration is expected there. In this paper, we perform simulations of the same systems, adding dust particles that suffer gas drag and modulate gas flow via the back-reaction of the gas drag (dust drag). We use O(10^7) super-particles, each of which represents many dust particles with sizes of centimeter to meter. We have found that the dust drag suppresses turbulent motion to decrease the velocity dispersion of the dust particles while it broadens the dust concentrated regions to limit peaky dust concentration, compared with the simulation without the dust drag. We found that reduction in the velocity dispersion) is dominated over the suppression in particle concentration. For meter-size particles with the friction time ~1/Omega, where Omega is Keplerian frequency, the gravitational instability of the dust particles that may lead to planetesimal formation is expected. For such a situation, we further introduced the self-gravity of dust particles to the simulation to demonstrate that several gravitationally bound clumps are actually formed. Through analytical arguments, we found that the planetesimal formation from meter-sized dust particles can be possible at ~5AU, if dust spatial density is a few times larger than that in the minimum mass solar nebula.
rate research
Read More
Recent years have seen growing interest in the streaming instability as a candidate mechanism to produce planetesimals. However, these investigations have been limited to small-scale simulations. We now present the results of a global protoplanetary disk evolution model that incorporates planetesimal formation by the streaming instability, along with viscous accretion, photoevaporation by EUV, FUV, and X-ray photons, dust evolution, the water ice line, and stratified turbulence. Our simulations produce massive (60-130 $M_oplus$) planetesimal belts beyond 100 au and up to $sim 20 M_oplus$ of planetesimals in the middle regions (3-100 au). Our most comprehensive model forms 8 $M_oplus$ of planetesimals inside 3 au, where they can give rise to terrestrial planets. The planetesimal mass formed in the inner disk depends critically on the timing of the formation of an inner cavity in the disk by high-energy photons. Our results show that the combination of photoevaporation and the streaming instability are efficient at converting the solid component of protoplanetary disks into planetesimals. Our model, however, does not form enough early planetesimals in the inner and middle regions of the disk to give rise to giant planets and super-Earths with gaseous envelopes. Additional processes such as particle pileups and mass loss driven by MHD winds may be needed to drive the formation of early planetesimal generations in the planet forming regions of protoplanetary disks.
A critical step toward the emergence of planets in a protoplanetary disk consists in accretion of planetesimals, bodies 1-1000 km in size, from smaller disk constituents. This process is poorly understood partly because we lack good observational constraints on the complex physical processes that contribute to planetesimal formation. In the outer solar system, the best place to look for clues is the Kuiper belt, where icy planetesimals survived to this day. Here we report evidence that Kuiper belt planetesimals formed by the streaming instability, a process in which aerodynamically concentrated clumps of pebbles gravitationally collapse into 100-km-class bodies. Gravitational collapse was previously suggested to explain the ubiquity of equal-size binaries in the Kuiper belt. We analyze new hydrodynamical simulations of the streaming instability to determine the model expectations for the spatial orientation of binary orbits. The predicted broad inclination distribution with 80% of prograde binary orbits matches the observations of trans-Neptunian binaries. The formation models which imply predominantly retrograde binary orbits can be ruled out. Given its applicability over a broad range of protoplanetary disk conditions, it is expected that the streaming instability seeded planetesimal formation also elsewhere in the solar system, and beyond.
The magnetorotational (MRI) dynamo has long been considered one of the possible drivers of turbulent angular momentum transport in astrophysical accretion disks. However, various numerical results suggest that this dynamo may be difficult to excite in the astrophysically relevant regime of magnetic Prandtl number (Pm) significantly smaller than unity, for reasons currently not well understood. The aim of this article is to present the first results of an ongoing numerical investigation of the role of both linear and nonlinear dissipative effects in this problem. Combining a parametric exploration and an energy analysis of incompressible nonlinear MRI dynamo cycles representative of the transitional dynamics in large aspect ratio shearing boxes, we find that turbulent magnetic diffusion makes the excitation and sustainment of this dynamo at moderate magnetic Reynolds number (Rm) increasingly difficult for decreasing Pm. This results in an increase in the critical Rm of the dynamo for increasing kinematic Reynolds number (Re), in agreement with earlier numerical results. Given its very generic nature, we argue that turbulent magnetic diffusion could be an important determinant of MRI dynamo excitation in disks, and may also limit the efficiency of angular momentum transport by MRI turbulence in low Pm regimes.
Debris discs are commonly swathed in gas which can be observed in UV, in fine structure lines in FIR, and in resolved maps of CO emission. Carbon and oxygen are overabundant in such gas, but it is severely depleted in hydrogen. As a consequence, its ionisation fraction is remarkably high, suggesting magnetohydrodynamic (MHD) processes may be important. In particular, the gas may be subject to the magnetorotational instability (MRI), and indeed recent modelling of $beta$ Pictoris requires an anomalous viscosity to explain the gass observed radial structure. In this paper we explore the possibility that the MRI is active in debris-disc gas and responsible for the observed mass transport. We find that non-ideal MHD and dust-gas interactions play a subdominant role, and that linear instability is viable at certain radii. However, owing to low gas densities, the outer parts of the disc could be stabilised by a weak ambient magnetic field, though it is difficult to constrain such a field. Even if the MRI is stabilised by too strong a field, a magnetocentrifugal wind may be launched in its place and this could lead to equivalent (non-turbulent) transport. Numerical simulations of the vertically stratified MRI in conditions appropriate to the debris disc gas should be able to determine the nature of the characteristic behaviour at different radii, and decide on the importance of the MRI (and MHD more generally) on the evolution of these discs.
The collapse of dust particle clouds directly to km-sized planetesimals is a promising way to explain the formation of planetesimals, asteroids and comets. In the past, this collapse has been studied in stratified shearing box simulations with super-solar dust-to-gas ratio epsilon, allowing for streaming instability (SI) and gravitational collapse. This paper studies the non-stratified SI under dust-to-gas ratios from epsilon=0.1 up to epsilon=1000 without self-gravity. The study covers domain sizes of L=0.1 H, 0.01 H and 0.001 H, in terms of gas disk scale height H, using the PencilCode. They are performed in radial-azimuthal (2-d) and radial-vertical (2.5-d) extent. The used particles of St=0.01 and 0.1 mark the upper end of the expected dust growth. SI-activity is found up to very high dust-to-gas ratios, providing fluctuations in the local dust-to-gas ratios and turbulent particle diffusion delta. We find an SI-like instability that operates in r-varphi even when vertical modes are suppressed. This new azimuthal streaming instability (aSI) shows similar properties and appearance as the SI. Both, SI and aSI, show diffusivity at epsilon=100 only to be two orders of magnitude lower than at epsilon=1, suggesting a delta ~ epsilon^{-1} relation that is shallow around epsilon = 1. The (a)SI ability to concentrate particles is found to be uncorrelated with its strength in particle turbulence. Finally, we performed a resolution study to test our findings of the aSI. This paper stresses out the importance of properly resolving the (a)SI at high dust-to-gas ratios and planetesimal collapse simulations, leading else wise to potentially incomplete results.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
