skip to main content


Title: A Dynamical Systems Approach to the Theory of Circumbinary Orbits in the Circular Restricted Problem
Abstract

To better understand the orbital dynamics of exoplanets around close binary stars, i.e., circumbinary planets (CBPs), we applied techniques from dynamical systems theory to a physically motivated set of solutions in the Circular Restricted Three-Body Problem (CR3BP). We applied Floquet theory to characterize the linear dynamical behavior—static, oscillatory, or exponential—surrounding planar circumbinary periodic trajectories (limit cycles). We computed prograde and retrograde limit cycles and analyzed their geometries, stability bifurcations, and dynamical structures. Orbit and stability calculations are exact computations in the CR3BP and reproducible through the open-source Python packagepyraa. The periodic trajectories (doi.org/10.5281/zenodo.7532982) produce a set of noncrossing, dynamically cool circumbinary orbits conducive to planetesimal growth. For mass ratiosμ∈ [0.01, 0.50], we found recurring features in the prograde families. These features include (1) an innermost near-circular trajectory, inside which solutions have resonant geometries, (2) an innermost stable trajectory (ac≈ 1.61 − 1.85abin) characterized by a tangent bifurcating limit cycle, and (3) a region of dynamical instability (a≈ 2.1abin; Δa≈ 0.1abin), the exclusion zone, bounded by a pair of critically stable trajectories bifurcating limit cycles. The exterior boundary of the exclusion zone is consistent with prior determinations ofacaround a circular binary. We validate our analytic results withN-body simulations and apply them to the Pluto–Charon system. The absence of detected CBPs in the inner stable region, between the prograde exclusion zone andac, suggests that the exclusion zone may inhibit the inward migration of CBPs.

 
more » « less
PAR ID:
10399921
Author(s) / Creator(s):
;
Publisher / Repository:
DOI PREFIX: 10.3847
Date Published:
Journal Name:
The Astronomical Journal
Volume:
165
Issue:
4
ISSN:
0004-6256
Format(s):
Medium: X Size: Article No. 140
Size(s):
Article No. 140
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    We present fully relativistic predictions for the electromagnetic emission produced by accretion disks surrounding spinning and nonspinning supermassive binary black holes on the verge of merging. We use the codeBothrosto post-process data from 3D general relativistic magnetohydrodynamic simulations via ray-tracing calculations. These simulations model the dynamics of a circumbinary disk and the mini-disks that form around two equal-mass black holes orbiting each other at an initial separation of 20 gravitational radii, and evolve the system for more than 10 orbits in the inspiral regime. We model the emission as the sum of thermal blackbody radiation emitted by an optically thick accretion disk and a power-law spectrum extending to hard X-rays emitted by a hot optically thin corona. We generate time-dependent spectra, images, and light curves at various frequencies to investigate intrinsic periodic signals in the emission, as well as the effects of the black hole spin. We find that prograde black hole spin makes mini-disks brighter since the smaller innermost stable circular orbit angular momentum demands more dissipation before matter plunges to the horizon. However, compared to mini-disks in larger separation binaries with spinning black holes, our mini-disks are less luminous: unlike those systems, their mass accretion rate is lower than in the circumbinary disk, and they radiate with lower efficiency because their inflow times are shorter. Compared to a single black hole system matched in mass and accretion rate, these binaries have spectra noticeably weaker and softer in the UV. Finally, we discuss the implications of our findings for the potential observability of these systems.

     
    more » « less
  2. Abstract

    We analyze accretion-rate time series for equal-mass binaries in coplanar gaseous disks spanning a continuous range of orbital eccentricities up to 0.8 for both prograde and retrograde systems. The dominant variability timescales match those of previous investigations; the binary orbital period is dominant for prograde binaries withe≳ 0.1, with a 5 × longer “lump” period taking over fore≲ 0.1. This lump period fades and drops from 5 × to 4.5 × the binary period aseapproaches 0.1, where it vanishes. For retrograde orbits, the binary orbital period dominates ate≲ 0.55 and is accompanied by a 2 × longer timescale periodicity at higher eccentricities. The shape of the accretion-rate time series varies with binary eccentricity. For prograde systems, the orientation of an eccentric disk causes periodic trading of accretion between the binary components in a ratio that we report as a function of binary eccentricity. We present a publicly available tool,binlite, that can rapidly (≲0.01 s) generate templates for the accretion-rate time series onto either binary component for choice of binary eccentricity below 0.8. As an example use case, we build lightcurve models where the accretion rate through the circumbinary disk and onto each binary component sets contributions to the emitted specific flux. We combine these rest-frame, accretion-variability lightcurves with observer-dependent Doppler boosting and binary self-lensing. This allows a flexible approach to generating lightcurves over a wide range of binary and observer parameter space. We envisionbinliteas the access point to a living database that will be updated with state-of-the-art hydrodynamical calculations as they advance.

     
    more » « less
  3. Abstract

    Many studies have recently documented the orbital response of eccentric binaries accreting from thin circumbinary disks, characterizing the change in the binary semimajor axis and eccentricity. We extend these calculations to include the precession of the binary’s longitude of periapse induced by the circumbinary disk, and we characterize this precession continuously with binary eccentricityebfor equal mass components. This disk-induced apsidal precession is prograde with a weak dependence on the binary eccentricity wheneb≲ 0.4 and decreases approximately linearly foreb≳ 0.4; yet at allebbinary precession is faster than the rates of change to the semimajor axis and eccentricity by an order of magnitude. We estimate that such precession effects are likely most important for subparsec separated binaries with masses ≲107M, like LISA precursors. We find that accreting, equal-mass LISA binaries withM< 106M(and the most massiveM∼ 107Mbinaries out toz∼ 3) may acquire a detectable phase offset due to the disk-induced precession. Moreover, disk-induced precession can compete with general relativistic precession in a vacuum, making it important for observer-dependent electromagnetic searches for accreting massive binaries—like Doppler boost and binary self-lensing models—after potentially only a few orbital periods.

     
    more » « less
  4. 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
  5. Abstract

    Hydrodynamical interactions between binaries and circumbinary disks (CBDs) play an important role in a variety of astrophysical systems, from young stellar binaries to supermassive black hole binaries. Previous simulations of CBDs have mostly employed locally isothermal equations of state. We carry out 2D viscous hydrodynamic simulations of CBDs around equal-mass, circular binaries, treating the gas thermodynamics by thermal relaxation toward equilibrium temperature (the constant-βcooling ansatz, whereβis the cooling time in units of the local Keplerian time). As an initial study, we use the grid-based codeAthena++on a polar grid, covering an extended disk outside the binary co-orbital region. We find that with a longer cooling time, the accretion variability is gradually suppressed, and the morphology of the CBD becomes more symmetric. The disk also shows evidence of hysteresis behavior depending on the initial conditions. Gas cooling also affects the rate of angular momentum transfer between the binary and the CBD, where given our adopted disk thickness and viscosity (H/r∼ 0.1 andα∼ 0.1), the binary orbit expands while undergoing accretion for mostβvalues between 0 and 4.0 except over a narrow range of intermediateβvalues. The validity of using a polar grid excising the central domain is also discussed.

     
    more » « less