Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userOn the Dynamics of the Inclination Instability
176
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
Axisymmetric disks of eccentric Kepler orbits are vulnerable to an instability which causes orbits to exponentially grow in inclination, decrease in eccentricity, and cluster in their angle of pericenter. Geometrically, the disk expands to a cone shape which is asymmetric about the mid-plane. In this paper, we describe how secular gravitational torques between individual orbits drive this inclination instability. We derive growth timescales for a simple two-orbit model using a Gauss $N$-ring code, and generalize our result to larger $N$ systems with $N$-body simulations. We find that two-body relaxation slows the growth of the instability at low $N$ and that angular phase coverage of orbits in the disk is important at higher $N$. As $N to infty$, the e-folding timescale converges to that expected from secular theory.
rate research
Read More
The interaction between a planet located in the inner region of a disc and the warped outer region is studied. We consider the stage of evolution after the planet has cleared-out a gap, so that the planetary orbit evolves only under the gravitational potential from the disc. We develop a secular analysis and compute the evolution of the orbital elements by solving Lagranges equations valid to second order in the eccentricity. We also perform numerical simulations with the full disc potential. In general, the interaction between the disc and the planet leads to the precession of the orbit. The orbital plane therefore becomes tilted relative to the discs inner parts, with no change in the eccentricity. When the inclination approaches 90 degrees, there is an instability and the eccentricity increases. In this case, both the inclination and the eccentricity develop large variations, with the orbit becoming retrograde. As the eccentricity reaches high values, we would expect tidal capture on a short orbit of the planet by the star to occur. This instability happens when the disc is severely warped, or if there is a significant amount of mass in a ring inclined by at least 45 degrees relative to the initial orbital plane. The inclination of the orbit does not depend on the semimajor axis nor on the planets mass. However, for a significant inclination to be generated on a timescale of at most a few Myr, the planet should be beyond the snow line. The process described here would therefore produce two distinct populations of inclined planets: one with objects beyond the snow line with at most moderate eccentricities, and another with objects on short circularized orbits.
Disks of bodies orbiting a much more massive central object are extremely common in astrophysics. When the orbits comprising such disks are eccentric, we show they are susceptible to a new dynamical instability. Gravitational forces between bodies in the disk drive exponential growth of their orbital inclinations and clustering in their angles of pericenter, expanding an initially thin disk into a conical shape by giving each orbit an identical tilt with respect to the disk plane. This new instability dynamically produces the unusual distribution of orbits observed for minor planets beyond Neptune, suggesting that the instability has shaped the outer Solar System. It also implies a large initial disk mass (1-10 Earth masses) of scattered bodies at hundreds of AU; we predict increasing numbers of detections of minor planets clustered in their angles of pericenter with high inclinations.
The Streaming Instability (SI) is a mechanism to concentrate solids in protoplanetary disks. Nonlinear particle clumping from the SI can trigger gravitational collapse into planetesimals. To better understand the numerical robustness of the SI, we perform a suite of vertically-stratified 3D simulations with fixed physical parameters known to produce strong clumping. We vary the numerical implementation, namely the computational domain size and the vertical boundary conditions (vBCs), comparing newly-implemented outflow vBCs to the previously-used periodic and reflecting vBCs. We find strong particle clumping by the SI is mostly independent of the vBCs. However, peak particle densities are higher in larger simulation domains due to a larger particle mass reservoir. We report SI-triggered zonal flows, i.e., azimuthally-banded radial variations of gas pressure. These structures have low amplitudes, insufficient to halt particle radial drift, confirming that particle trapping in gas pressure maxima is not the mechanism of the SI. We find that outflow vBCs produce artificially large gas outflow rates at vertical boundaries. However, the outflow vBCs reduce artificial reflections at vertical boundaries, allowing more particle sedimentation, and showing less temporal variation and better convergence with box size. The radial spacing of dense particle filaments is $sim0.15$ gas scale heights ($H$) for all vBCs, which sets the feeding zone for planetesimal growth in self-gravitating simulations. Our results validate the use of the outflow vBCs in SI simulations, even with vertical boundaries close ($leq 0.4H$) to the disk midplane. Overall, our study demonstrates the numerical robustness of nonlinear particle clumping by the SI.
Several recent studies have suggested that circumstellar disks in young stellar binaries may be driven into misalignement with their host stars due to secular gravitational interactions between the star, disk and the binary companion. The disk in such systems is twisted/warped due to the gravitational torques from the oblate central star and the external companion. We calculate the disk warp profile, taking into account of bending wave propagation and viscosity in the disk. We show that for typical protostellar disk parameters, the disk warp is small, thereby justifying the flat-disk approximation adopted in previous theoretical studies. However, the viscous dissipation associated with the small disk warp/twist tends to drive the disk toward alignment with the binary or the central star. We calculate the relevant timescales for the alignment. We find the alignment is effective for sufficiently cold disks with strong external torques, especially for systems with rapidly rotating stars, but is ineffective for the majority of star-disk-binary systems. Viscous warp driven alignment may be necessary to account for the observed spin-orbit alignment in multi-planet systems if these systems are accompanied by an inclined binary companion.
The streaming instability is a leading candidate mechanism to explain the formation of planetesimals. Yet, the role of this instability in the driving of turbulence in protoplanetary disks, given its fundamental nature as a linear hydrodynamical instability, has so far not been investigated in detail. We study the turbulence that is induced by the streaming instability as well as its interaction with the vertical shear instability. For this purpose, we employ the FLASH Code to conduct two-dimensional axisymmetric global disk simulations spanning radii from $1$ au to $100$ au, including the mutual drag between gas and dust as well as the radial and vertical stellar gravity. If the streaming instability and the vertical shear instability start their growth at the same time, we find the turbulence in the dust mid-plane layer to be primarily driven by the streaming instability. It gives rise to vertical gas motions with a Mach number of up to ${sim}10^{-2}$. The dust scale height is set in a self-regulatory manner to about $1%$ of the gas scale height. In contrast, if the vertical shear instability is allowed to saturate before the dust is introduced into our simulations, then it continues to be the main source of the turbulence in the dust layer. The vertical shear instability induces turbulence with a Mach number of ${sim}10^{-1}$ and thus impedes dust sedimentation. Nonetheless, we find the vertical shear instability and the streaming instability in combination to lead to radial dust concentration in long-lived accumulations which are significantly denser than those formed by the streaming instability alone. Thus, the vertical shear instability may promote planetesimal formation by creating weak overdensities that act as seeds for the streaming instability.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
