A test particle orbit around an eccentric binary has two stationary states in which there is no nodal precession: coplanar and polar. Nodal precession of a misaligned test particle orbit centres on one of these stationary states. A lowmass circumbinary disc undergoes the same precession and moves towards one of these states through dissipation within the disc. For a massive particle orbit, the stationary polar alignment occurs at an inclination less than 90°, which is the progradepolar stationary inclination. A sufficiently high angular momentum particle has an additional higher inclination stationary state, the retrogradepolar stationary inclination. Misaligned particle orbits close to the retrogradepolar stationary inclination are not nested like the orbits close to the other stationary points. We investigate the evolution of a gas disc that begins close to the retrogradepolar stationary inclination. With hydrodynamical disc simulations, we find that the disc moves through the unnested crescent shape precession orbits and eventually moves towards the progradepolar stationary inclination, thus increasing the parameter space over which circumbinary discs move towards polar alignment. If protoplanetary discs form with an isotropic orientation relative to the binary orbit, then polar discs may be more common than coplanar discs around eccentric binaries, even formore »
Mutually misaligned circumbinary planets may form in a warped or broken gas disk or from later planet–planet interactions. With numerical simulations and analytic estimates we explore the dynamics of two circumbinary planets with a large mutual inclination. A coplanar inner planet causes prograde apsidal precession of the binary and the stationary inclination for the outer planet is higher for larger outer planet orbital radius. In this case a coplanar outer planet always remains coplanar. On the other hand, a polar inner planet causes retrograde apsidal precession of the binary orbit and the stationary inclination is smaller for larger outer planet orbital radius. For a range of outer planet semimajor axes, an initially coplanar orbit is librating meaning that the outer planet undergoes large tilt oscillations. Circumbinary planets that are highly inclined to the binary are difficult to detect—it is unlikely for a planet to have an inclination below the transit detection limit in the presence of a polar inner planet. These results suggest that there could be a population of circumbinary planets that are undergoing large tilt oscillations.
 Publication Date:
 NSFPAR ID:
 10399857
 Journal Name:
 The Astrophysical Journal Letters
 Volume:
 945
 Issue:
 1
 Page Range or eLocationID:
 Article No. L11
 ISSN:
 20418205
 Publisher:
 DOI PREFIX: 10.3847
 Sponsoring Org:
 National Science Foundation
More Like this

ABSTRACT 
Abstract Systems with ultrashortperiod (USP) planets tend to possess larger mutual inclinations compared to those with planets located farther from their host stars. This could be explained due to precession caused by stellar oblateness at early times when the host star was rapidly spinning. However, stellar oblateness reduces over time due to the decrease in the stellar rotation rate, and this may further shape the planetary mutual inclinations. In this work, we investigate in detail how the final mutual inclination varies under the effect of a decreasing J 2 . We find that different initial parameters (e.g., the magnitude of J 2 and planetary inclinations) will contribute to different final mutual inclinations, providing a constraint on the formation mechanisms of USP planets. In general, if the inner planets start in the same plane as the stellar equator (or coplanar while misaligned with the stellar spin axis), the mutual inclination decreases (or increases then decreases) over time due to the decay of the J 2 moment. This is because the inner orbit typically possesses less orbital angular momentum than the outer ones. However, if the outer planet is initially aligned with the stellar spin while the inner one is misaligned, themore »

ABSTRACT We study the stationary points of the hierarchical three body problem in the planetary limit (m1, m2 ≪ m0) at both the quadrupole and octupole orders. We demonstrate that the extension to octupole order preserves the principal stationary points of the quadrupole solution in the limit of small outer eccentricity e2 but that new families of stable fixed points occur in both prograde and retrograde cases. The most important new equilibria are those that branch off from the quadrupolar solutions and extend to large e2. The apsidal alignment of these families is a function of mass and inner planet eccentricity, and is determined by the relative directions of precession of ω1 and ω2 at the quadrupole level. These new equilibria are also the most resilient to the destabilizing effects of relativistic precession. We find additional equilibria that enable libration of the inner planet argument of pericentre in the limit of radial orbits and recover the nonlinear analogue of the Laplace–Lagrange solutions in the coplanar limit. Finally, we show that the chaotic diffusion and orbital flips identified with the eccentric Kozai–Lidov mechanism and its variants can be understood in terms of the stationary points discussed here.

ABSTRACT The asteroid belt is characterized by an extreme low total mass of material on dynamically excited orbits. The Nice model explains many peculiar qualities of the Solar system, including the belt’s excited state, by invoking an orbital instability between the outer planets. However, previous studies of the Nice model’s effect on the belt’s structure struggle to reproduce the innermost asteroids’ orbital inclination distribution. Here, we show how the final phase of giant planet migration sculpts the asteroid belt, in particular its inclination distribution. As interactions with leftover planetesimals cause Saturn to move away from Jupiter, its rate of orbital precession slows as the two planets’ mutual interactions weaken. When the planets approach their modern separation, where Jupiter completes just short of five orbits for every two of Saturn’s, Jupiter’s eccentric forcing on Saturn strengthens. We use numerical simulations to show that the absence of asteroids with orbits that precess between 24 and 28 arcsec yr−1 is related to the inclination problem. As Saturn’s precession speeds back up, highinclination asteroids are excited on to planet crossing orbits and removed from the inner main belt. Through this process, the asteroid belt’s orbital structure is reshaped, leading to markedly improved simulation outcomes.

ABSTRACT We describe the first gridbased simulations of the polar alignment of a circumbinary disc. We simulate the evolution of an inclined disc around an eccentric binary using the gridbased code athena++ . The use of a gridbased numerical code allows us to explore lower disc viscosities than have been examined in previous studies. We find that the disc aligns to a polar orientation when the α viscosity is high, while discs with lower viscosity nodally precess with little alignment over 1000 binary orbital periods. The timescales for polar alignment and disc precession are compared as a function of disc viscosity, and are found to be in agreement with previous studies. At very low disc viscosities (e.g. α = 10−5), anticyclonic vortices are observed along the inner edge of the disc. These vortices can persist for thousands of binary orbits, creating azimuthally localized overdensities and multiple pairs of spiral arms. The vortex is formed at ∼3–4 times the binary semimajor axis, close to the inner edge of the disc, and orbits at roughly the local Keplerian speed. The presence of a vortex in the disc may play an important role in the evolution of circumbinary systems, such as driving episodic accretionmore »