An official website of the United States government
Abstract 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.
more »
« less
The stationary points of the hierarchical three-body problem
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 non-linear 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.
more »
« less
- Award ID(s):
- 1739160
- PAR ID:
- 10211967
- Date Published:
- Journal Name:
- Monthly Notices of the Royal Astronomical Society
- Volume:
- 499
- Issue:
- 2
- ISSN:
- 0035-8711
- Page Range / eLocation ID:
- 1682 to 1700
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Aims. We analyze the behavior of the argument of pericenter ω 2 of an outer particle in the elliptical restricted three-body problem, focusing on the ω 2 resonance or inverse Lidov-Kozai resonance. Methods. First, we calculated the contribution of the terms of quadrupole, octupole, and hexadecapolar order of the secular approximation of the potential to the outer particle’s ω 2 precession rate (d ω 2 ∕d τ ). Then, we derived analytical criteria that determine the vanishing of the ω 2 quadrupole precession rate (d ω 2 /d τ ) quad for different values of the inner perturber’s eccentricity e 1 . Finally, we used such analytical considerations and described the behavior of ω 2 of outer particles extracted from N-body simulations developed in a previous work. Results. Our analytical study indicates that the values of the inclination i 2 and the ascending node longitude Ω 2 associated with the outer particle that vanish (d ω 2 /d τ ) quad strongly depend on the eccentricity e 1 of the inner perturber. In fact, if e 1 < 0.25 (>0.40825), (d ω 2 /d τ ) quad is only vanished for particles whose Ω 2 circulates (librates). For e 1 between 0.25 and 0.40825, (d ω 2 /d τ ) quad can be vanished for any particle for a suitable selection of pairs (Ω 2 , i 2 ). Our analysis of the N-body simulations shows that the inverse Lidov-Kozai resonance is possible for small, moderate, and high values of e 1 . Moreover, such a resonance produces distinctive features in the evolution of a particle in the (Ω 2 , i 2 ) plane. In fact, if ω 2 librates and Ω 2 circulates, the extremes of i 2 at Ω 2 = 90° and 270° do not reach the same value, while if ω 2 and Ω 2 librate, the evolutionary trajectory of the particle in the (Ω 2 , i 2 ) plane shows evidence of an asymmetry with respect to i 2 = 90°. The evolution of ω 2 associated with the outer particles of the N-body simulations can be very well explained by the analytical criteria derived in our investigation.more » « less
-
Abstract Multiplanetary systems are prevalent in our Galaxy. The long-term stability of such systems may be disrupted if a distant inclined companion excites the eccentricity and inclination of the inner planets via the eccentric Kozai–Lidov mechanism. However, the star–planet and the planet–planet interactions can help stabilize the system. In this work, we extend the previous stability criterion that only considered the companion–planet and planet–planet interactions by also accounting for short-range forces or effects, specifically, relativistic precession induced by the host star. A general analytical stability criterion is developed for planetary systems with N inner planets and a relatively distant inclined perturber by comparing precession rates of relevant dynamical effects. Furthermore, we demonstrate as examples that in systems with two and three inner planets, the analytical criterion is consistent with numerical simulations using a combination of Gauss’s averaging method and direct N -body integration. Finally, the criterion is applied to observed systems, constraining the orbital parameter space of a possible undiscovered companion. This new stability criterion extends the parameter space in which an inclined companion of multiplanet systems can inhabit.more » « less
-
The most general bound binary black hole system has an eccentric orbit and precessing spins. The detection of such a system with significant eccentricity close to the merger would be a clear signature of dynamical formation. In order to study such systems, it is important to be able to evolve their spins and eccentricity from the larger separations at which the binary formed to the smaller separations at which it is detected, or vice versa. Knowledge of the precessional evolution of the binary’s orbital angular momentum can also be used to twist up aligned-spin eccentric waveform models to create a spin-precessing eccentric waveform model. In this paper, we present a new publicly available code to evolve eccentric, precessing binary black holes using orbit-averaged post-Newtonian (PN) equations from the literature. The spin-precession dynamics is 2PN accurate, i.e., with the leading spin-orbit and spin-spin corrections. The evolution of orbital parameters (orbital frequency, eccentricity, and periastron precession), which follow the quasi-Keplerian parametrization, is 3PN accurate in the point particle terms and includes the leading order spin-orbit and spin-spin effects. All the spin-spin terms include the quadrupole-monopole interaction. The eccentricity enhancement functions in the fluxes use the high-accuracy hyperasymptotic expansions from Loutrel and Yunes []. We discuss various features of the code and study the evolution of the orbital and spin-precession parameters of eccentric, precessing binary black holes. In particular, we study the dependence of the spin morphologies on eccentricity, where we find that the transition point from one spin morphology to another can depend nonmonotonically on eccentricity, and the fraction of binaries in a given morphology at a given point in the evolution of a population depends on the instantaneous eccentricity.more » « less
-
Abstract WASP-107 b seems to be a poster child of the long-suspected high-eccentricity migration scenario. It is on a 5.7 day, polar orbit. The planet is Jupiter-like in radius but Neptune-like in mass with exceptionally low density. WASP-107 c is on a 1100 day,e= 0.28 orbit with at least Saturn mass. Planet b may still have a residual eccentricity of 0.06 ± 0.04: the ongoing tidal dissipation leads to the observed internally heated atmosphere and hydrodynamic atmospheric erosion. We present a population synthesis study coupling octupole Lidov–Kozai oscillations with various short-range forces, while simultaneously accounting for the radius inflation and tidal disruption of the planet. We find that a high-eccentricity migration scenario can successfully explain nearly all observed system properties. Our simulations further suggest that the initial location of WASP-107 b at the onset of migration is likely within the snowline (<0.5 au). More distant initial orbits usually lead to tidal disruption or orbit crossing. WASP-107 b most likely lost no more than 20% of its mass during the high-eccentricity migration, i.e., it did not form as a Jupiter-mass object. More vigorous tidally induced mass loss leads to disruption of the planet during migration. We predict that the current-day mutual inclination between the planets b and c is substantial: at least 25°–55°, which may be tested with future Gaia astrometric observations. Knowing the current-day mutual inclination may further constrain the initial orbit of planet b. We suggest that the proposed high-eccentricity migration scenario of WASP-107 may be applicable to HAT-P-11, GJ-3470, HAT-P-18, and GJ-436, which have similar orbital architectures.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1093/mnras/staa2602
