Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userThe effects of disc warping on the inclination of planetary orbits
151
0
0.0
(
0
)
Authors
Caroline Terquem
Ask ChatGPT about the research
No Arabic abstract
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.
rate research
Read More
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.
Many exoplanets are discovered in binary star systems in internal or in circumbinary orbits. Whether the planet can be habitable or not depends on the possibility to maintain liquid water on its surface, and therefore on the luminosity of its host stars and on the dynamical properties of the planetary orbit. The trajectory of a planet in a double star system can be determined, approximating stars and planets with point masses, by solving numerically the equations of motion of the classical three-body system. In this study, we analyze a large data set of planetary orbits, made up with high precision long integration at varying: the mass of the planet, its distance from the primary star, the mass ratio for the two stars in the binary system, and the eccentricity of the star motion. To simulate the gravitational dynamics, we use a 15th order integration scheme (IAS15, available within the REBOUND framework), that provides an optimal solution for long-term integration. In our data analysis, we evaluate if an orbit is stable or not and also provide the statistics of different types of instability: collisions with the primary or secondary star and planets ejected away from the binary star system. Concerning the stability, we find a significant number of orbits that are only marginally stable, according to the classification introduced by Musielak et al. 2005. For planets of negligible mass, we estimate the critical semi-major axis $a_c$ as a function of the mass ratio and the eccentricity of the binary, in agreement with the results of Holman and Wiegert 1999. However, we find that for very massive planets (Super-Jupiters) the critical semi-major axis decrease in some cases by a few percent, compared to cases in which the mass of the planet is negligible.
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.
Recent observations of several protoplanetary discs have found evidence of departures from flat, circular motion in the inner regions of the disc. One possible explanation for these observations is a disc warp, which could be induced by a planet on a misaligned orbit. We present three-dimensional numerical simulations of the tidal interaction between a protoplanetary disc and a misaligned planet. For low planet masses we show that our simulations accurately model the evolution of inclined planet orbit (up to moderate inclinations). For a planet massive enough to carve a gap, the disc is separated into two components and the gas interior and exterior to the planet orbit evolve separately, forming an inner and outer disc. Due to the inclination of the planet, a warp develops across the planet orbit such that there is a relative tilt and twist between these discs. We show that when other parameters are held constant, the relative inclination that develops between the inner and outer disc depends on the outer radius of the total disc modelled. For a given disc mass, our results suggest that the observational relevance of the warp depends more strongly on the mass of the planet rather than the inclination of the orbit.
Numerical simulations of planets embedded in protoplanetary gaseous discs are a precious tool for studying the planetary migration ; however, some approximations have to be made. Most often, the selfgravity of the gas is neglected. In that case, it is not clear in the literature how the material inside the Roche lobe of the planet should be taken into account. Here, we want to address this issue by studying the influence of various methods so far used by different authors on the migration rate. We performed high-resolution numerical simulations of giant planets embedded in discs. We compared the migration rates with and without gas selfgravity, testing various ways of taking the circum-planetary disc (CPD) into account. Different methods lead to significantly different migration rates. Adding the mass of the CPD to the perturbing mass of the planet accelerates the migration. Excluding a part of the Hill sphere is a very touchy parameter that may lead to an artificial suppression of the type III, runaway migration. In fact, the CPD is smaller than the Hill sphere. We recommend excluding no more than a 0.6 Hill radius and using a smooth filter. Alternatively, the CPD can be given the acceleration felt by the planet from the rest of the protoplanetary disc. The gas inside the Roche lobe of the planet should be very carefully taken into account in numerical simulations without any selfgravity of the gas. The entire Hill sphere should not be excluded. The method used should be explicitly given. However, no method is equivalent to computing the full selfgravity of the gas.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
