We investigate the dynamics of a nonzero mass, circular orbit planet around an eccentric orbit binary for various values of the binary eccentricity, binary mass fraction, planet mass, and planet semi-major axis by means of numerical simulations. Previous studies investigated the secular dynamics mainly by approximate analytic methods. In the stationary inclination state, the planet and binary precess together with no change in relative tilt. For both prograde and retrograde planetary orbits, we explore the conditions for planetary orbital libration versus circulation and the conditions for stationary inclination. As was predicted by analytic models, for sufficiently high initial inclination, a prograde planet's orbit librates about the stationary tilted state. For a fixed binary eccentricity, the stationary angle is a monotonically decreasing function of the ratio of the planet-to-binary angular momentum j. The larger j, the stronger the evolutionary changes in the binary eccentricity and inclination. We also calculate the critical tilt angle that separates the circulating from the librating orbits for both prograde and retrograde planet orbits. The properties of the librating orbits and stationary angles are quite different for prograde versus retrograde orbits. The results of the numerical simulations are in very good quantitative agreement with the analytic models. Our results have implications for circumbinary planet formation and evolution.

We study the orbital stability of a non-zero mass, close-in circular orbit planet around an eccentric orbit binary for various initial values of the binary eccentricity, binary mass fraction, planet mass, planet semimajor axis, and planet inclination by means of numerical simulations that cover 5 × 104 binary orbits. For small binary eccentricity, the stable orbits that extend closest to the binary (most stable orbits) are nearly retrograde and circulating. For high binary eccentricity, the most stable orbits are highly inclined and librate near the so-called generalized polar orbit which is a stationary orbit that is fixed in the frame of the binary orbit. For more extreme mass ratio binaries, there is a greater variation in the size of the stability region (defined by initial orbital radius and inclination) with planet mass and initial inclination, especially for low binary eccentricity. For low binary eccentricity, inclined planet orbits may be unstable even at large orbital radii (separation ${\gt}5 \, a_{\rm b}$). The escape time for an unstable planet is generally shorter around an equal mass binary compared with an unequal mass binary. Our results have implications for circumbinary planet formation and evolution and will be helpful for understanding future circumbinary planet observations.

We investigate the formation mechanism for the observed nearly polar aligned (perpendicular to the binary orbital plane) debris ring around the eccentric orbit binary 99 Herculis. An initially inclined nonpolar debris ring or disc will not remain flat and will not evolve to a polar configuration, due to the effects of differential nodal precession that alter its flat structure. However, a gas disc with embedded well coupled solids around the eccentric binary may evolve to a polar configuration as a result of pressure forces that maintain the disc flatness and as a result of viscous dissipation that allows the disc to increase its tilt. Once the gas disc disperses, the debris disc is in a polar aligned state in which there is little precession. We use three-dimensional hydrodynamical simulations, linear theory, and particle dynamics to study the evolution of a misaligned circumbinary gas disc and explore the effects of the initial disc tilt, mass, and size. We find that for a wide range of parameter space, the polar alignment timescale is shorter than the lifetime of the gas disc. Using the observed level of alignment of 3 • from polar, we place an upper limit on the mass of the gas disc of about 0.014 M ⊙ at the time of gas dispersal. We conclude that the polar debris disc around 99 Her can be explained as the result of an initially moderately inclined gas disc with embedded solids. Such a disc may provide an environment for the formation of polar planets.

We investigate the orbital dynamics of circumbinary planetary systems with two planets around a circular or eccentric orbit binary. The orbits of the two planet are initially circular and coplanar to each other, but misaligned with respect to the binary orbital plane. The binary-planet and planet-planet interactions result in complex planet tilt oscillations. We use analytic models and numerical simulations to explore the effects of various values of the planet semi-major axes, binary eccentricity, and initial inclination. Around a circular orbit binary, secular tilt oscillations are driven by planet-planet interactions and are periodic. In that case, planets undergo mutual libration if close together and circulation if far apart with an abrupt transition at a critical separation. Around an eccentric orbit binary, secular tilt oscillations are driven by both planet-planet interactions and binary-planet interactions. These oscillations generally display more than one frequency and are generally not periodic. The transition from mutual planet libration to circulation is not sharp and there is a range of separations for which the planets are on orbits that are sometimes mutually librating and sometimes circulating. In addition, at certain separations, there are resonances for which tilt oscillations are complicated but periodic. For planets that are highly misaligned with respect to an eccentric orbit binary, there are stationary (non-oscillating) tilt configurations that are generalisations of polar configurations for the single planet case. Tilt oscillations of highly inclined planets occur for initial tilts that depart from the stationary configuration.

With n-body simulations we investigate the stability of tilted circumbinary planetary systems consisting of two nonzero mass planets. The planets are initially in circular orbits that are coplanar to each other, as would be expected if they form in a flat but tilted circumbinary gas disc and decouple from the disc within a time difference that is much less than the disc nodal precession period. We constrain the parameters of stable multiple planet circumbinary systems. Both planet-planet and planet-binary interactions can cause complex planet tilt oscillations which can destabilise the orbits of one or both planets. The system is considerably more unstable than the effects of these individual interactions would suggest, due to the interplay between these two interactions. The stability of the system is sensitive to the binary eccentricity, the orbital tilt and the semi-major axes of the two circumbinary planets. With an inner planet semi-major axis of 5 ab, where ab is semi-major axis of the binary, the system is generally stable if the outer planet is located at ≳ 8 ab, beyond the 2:1 mean motion resonance with the inner planet. For larger inner planet semi-major axis the system is less stable because the von-Zeipel–Kozai–Lidov mechanism plays a significant role, particularly for low binary-eccentricity cases. For the unstable cases, the most likely outcome is that one planet is ejected and the other remains bound on a highly eccentric orbit. Therefore we suggest that this instability is an efficient mechanism for producing free-floating planets.