An exosolar planetary system N-body simulator II

We present a general N-body exasolar system simulator in anticipation of upcoming searches for exoplanets and even exomoons by next generation telescopes such as James Webb Space Telescope. For habitable zones, traditionally defined by temperature, we here address the essential problem of dynamical stability of planetary orbits. Illustrative examples are presented on P-type orbits in stellar binary systems, that should be fairly common as in Kepler 16b. Specific attention is paid to reduced orbital lifetimes of exoplanets in the habitable zone by the stellar binary, that is propoesed by Maurice van Putten (2017). Especially, we focused on a classic work of complex three-body problem that is well known by Dvorak(1986). We charge his elliptic restricted three-body problem to extend unrestricted three-body problem to look into dynamical motions in view of circumbinary planet, furthermore, we suggest that opposite angular orientation of the planet is relative to the stability of orbits. In here, counter-rotation case is relatively more faster than co-rotation case for being stable. As a result, we find that various initial conditions and thresholds to approach dynamical stability and unstability with unexpectable isolated islands over enormous parameter space. Even, superkeplerian effect of binary is important to habitability of the exoplanet and we can verify that superfaster binary doesn't effect on th planet and increases survivality of planet around the binary.