Showing 1–50 of 54 results for author: Quarles, B

More Realistic Planetesimal Masses Alter Kuiper Belt Formation Models and Add Stochasticity

Authors: Nathan A. Kaib, Alex Parsells, Simon Grimm, Billy Quarles, Matthew S. Clement

Abstract: We perform simulations here that include the gravitational effects of the primordial planetesimal belt consisting of ~10^5 massive bodies. In our simulations, Neptune unlocks from resonance with the other giant planets and begins to migrate outward due to interactions with planetesimals before a planetary orbital instability is triggered, and afterward, residual Neptunian migration completes the f… ▽ More We perform simulations here that include the gravitational effects of the primordial planetesimal belt consisting of ~10^5 massive bodies. In our simulations, Neptune unlocks from resonance with the other giant planets and begins to migrate outward due to interactions with planetesimals before a planetary orbital instability is triggered, and afterward, residual Neptunian migration completes the formation of the modern Kuiper belt. Our present work exhibits a number of notable differences from prior work. First, Neptune's planetary resonance unlocking requires the Neptunian 3:2 mean motion resonance to sweep much of the primordial disk interior to 30 au prior to the giant planet instability. The pre-instability population of planetesimals is consequently lower in semimajor axis, eccentricity, and inclination, and this effect persists after the instability. Second, direct scattering between Pluto-mass bodies and other small bodies removes material from Neptunian resonances more efficiently than resonant dropout resulting from small changes in Neptune's semimajor axis during scattering between Pluto-mass bodies and Neptune. Thus, the primordial population of Pluto-mass bodies may be as few as ~200 objects. Finally, our simulation end states display a wide variety of orbital distributions, and clear relationships between final bulk Kuiper belt properties and Neptune's migration or initial planetesimal properties largely elude us. In particular, we find that the rapid, stochastic planetary orbital evolution occurring during the giant planet instability can significantly alter final Kuiper belt properties such as its inclination dispersion and the prominence of resonant populations. This complicates using modern Kuiper belt properties to confidently constrain early solar system events and conditions, including planetary orbital migration and the primordial Kuiper belt's characteristics. △ Less

Submitted 18 March, 2024; originally announced March 2024.

Comments: Accepted to Icarus; 20 pages, 8 figures, 4 tables

A Reply to: Large Exomoons unlikely around Kepler-1625 b and Kepler-1708 b

Authors: David Kipping, Alex Teachey, Daniel A. Yahalomi, Ben Cassese, Billy Quarles, Steve Bryson, Brad Hansen, Judit Szulágyi, Chris Burke, Kevin Hardegree-Ullman

Abstract: Recently, Heller & Hippke argued that the exomoon candidates Kepler-1625 b-i and Kepler-1708 b-i were allegedly 'refuted'. In this Matters Arising, we address these claims. For Kepler-1625 b, we show that their Hubble light curve is identical to that previously published by the same lead author, in which the moon-like dip was recovered. Indeed, our fits of their data again recover the moon-like di… ▽ More Recently, Heller & Hippke argued that the exomoon candidates Kepler-1625 b-i and Kepler-1708 b-i were allegedly 'refuted'. In this Matters Arising, we address these claims. For Kepler-1625 b, we show that their Hubble light curve is identical to that previously published by the same lead author, in which the moon-like dip was recovered. Indeed, our fits of their data again recover the moon-like dip with improved residuals than that obtained by Heller & Hippke. Their fits therefore appear to have somehow missed this deeper likelihood maximum, as well producing apparently unconverged posteriors. Consequently, their best-fitting moon is the same radius as the planet, Kepler-1625 b; a radically different signal from that which was originally claimed. The authors then inject this solution into the Kepler data and remark, as a point of concern, how retrievals obtain much higher significances than originally reported. However, this issue stems from the injection of a fundamentally different signal. We demonstrate that their Hubble light curve exhibits ~20% higher noise and discards 11% of the useful data, which compromises its ability to recover the subtle signal of Kepler-1625 b-i. For Kepler-1708 b-i it was claimed that the exomoon model's Bayes factor is highly sensitive to detrending choices, yielding reduced evidence with a biweight filter versus the original claim. We use their own i) detrended light curve and ii) biweight filter code to investigate these claims. For both, we recover the original moon signal, to even higher confidence than before. The discrepancy is explained by comparing to their quoted fit metrics, where we again demonstrate that the Heller & Hippke regression definitively missed the deeper likelihood maximum corresponding to Kepler-1708 b-i. We conclude that both candidates remain viable but certainly demand further observations. △ Less

Submitted 18 January, 2024; originally announced January 2024.

Comments: Under consideration by Nature Astronomy as Matters Arising

Stability of coorbital planets around binaries

Authors: Stefan Adelbert, Anna B. T. Penzlin, Christoph M. Schäfer, Wilhelm Kley, Billy Quarles, Rafael Sfair

Abstract: In previous hydrodynamical simulations, we found a mechanism for nearly circular binary stars, like Kepler-413, to trap two planets in a stable 1:1 resonance. Therefore, the stability of coorbital configurations becomes a relevant question for planet formation around binary stars. Here, we investigate the coorbital planet stability using a Kepler-413 analogue as example and then expanding the para… ▽ More In previous hydrodynamical simulations, we found a mechanism for nearly circular binary stars, like Kepler-413, to trap two planets in a stable 1:1 resonance. Therefore, the stability of coorbital configurations becomes a relevant question for planet formation around binary stars. Here, we investigate the coorbital planet stability using a Kepler-413 analogue as example and then expanding the parameters to study general n-body stability of planet pairs in eccentric horseshoe orbits around binaries. The stability is tested by evolving the planet orbits for $10^5$ binary periods with varying initial semi-major axes and planet eccentricities. The unstable region of a single circumbinary planet is used as a comparison to the investigated coorbital configurations in this work. We confirm previous findings on the stability of single planets and find a first order linear relation between orbit eccentricity and pericentre to identify stable orbits for various binary configurations. Such a linear relation is also found for the stability of 1:1 resonant planets around binaries. Stable orbits for eccentric horseshoe configurations exist with a pericentre closer than seven binary separations and, in the case of Kepler-413, the pericentre of the first stable orbit can be approximated by $r_{c,peri} = (2.88 e_p + 2.46) a_{bin}$. △ Less

Submitted 11 October, 2023; originally announced October 2023.

Comments: 9 pages, 10 figures, accepted for publication by A&A

Journal ref: A&A 680, A29 (2023)

Mass Derivation of planets K2-21b and K2-21c from Transit Timing Variations

Authors: Maryame El Moutamid, Kevin B. Stevenson, Billy Quarles, Nikole K. Lewis, Erik Petigura Daniel Fabrycky, Jacob L. Bean, Diana Dragomir, Kristin S. Sotzenvand Michael W. Werner

Abstract: While various indirect methods are used to detect exoplanets, one of the most effective and accurate methods is the transit method, which measures the brightness of a given star for periodic dips when an exoplanet is passing in front of the parent star. For systems with multiple transiting planets, the gravitational perturbations between planets affect their transit times. The difference in transi… ▽ More While various indirect methods are used to detect exoplanets, one of the most effective and accurate methods is the transit method, which measures the brightness of a given star for periodic dips when an exoplanet is passing in front of the parent star. For systems with multiple transiting planets, the gravitational perturbations between planets affect their transit times. The difference in transit times allows a measurement of the planet masses and orbital eccentricities. These parameters help speculating on the formation, evolution and stability of the system. Using Transit Timing Variations (TTVs), we measure the masses and eccentricities of two planets orbiting K2-21, a relatively bright K7 dwarf star. These two planets exhibit measurable TTVs, have orbital periods of about 9.32 days and 15.50 days, respectively, and a period ratio of about 1.66, which is relatively near to the 5:3 mean motion resonance. We report that the inner and outer planets in the K2-21 system have properties consistent with the presence of a hydrogen and helium dominated atmospheres, as we estimate their masses to be 1.59^{+0.52}_{-0.44} M_E and 3.88^{+1.22}_{-1.07} M_E and densities of 0.22^{+0.05}_{-0.04} rho_E and 0.34^{+0.08}_{-0.06} rho_E, respectively (M_E and rho_E are the mass and density of Earth, respectively). Our results show that the inner planet is less dense than the outer planet; one more counter-intuitive exoplanetary system such as Kepler-105, LTT 1445, TOI-175 and Kepler-279 systems. △ Less

Submitted 18 May, 2023; originally announced May 2023.

Statistical Analysis of the Dearth of Super-eccentric Jupiters in the Kepler Sample

Authors: Jonathan M. Jackson, Rebekah I. Dawson, Billy Quarles, Jiayin Dong

Abstract: Hot Jupiters may have formed in situ, or been delivered to their observed short periods through one of two categories of migration mechanisms: disk migration or high-eccentricity migration. If hot Jupiters were delivered by high-eccentricity migration, we would expect to observe some "super-eccentric" Jupiters in the process of migrating. We update a prediction for the number of super-eccentric Ju… ▽ More Hot Jupiters may have formed in situ, or been delivered to their observed short periods through one of two categories of migration mechanisms: disk migration or high-eccentricity migration. If hot Jupiters were delivered by high-eccentricity migration, we would expect to observe some "super-eccentric" Jupiters in the process of migrating. We update a prediction for the number of super-eccentric Jupiters we would expect to observe in the Kepler sample if all hot Jupiters migrated through high-eccentricity migration and estimate the true number observed by Kepler. We find that the observations fail to match the prediction from high-eccentricity migration with 94.3% confidence and show that high-eccentricity migration can account for at most ~62% of the hot Jupiters discovered by Kepler. △ Less

Submitted 14 December, 2022; originally announced December 2022.

Comments: 12 pages, 5 figures, 2 tables, submitted to AAS journals and revised following comments from referee

Moon-packing around an Earth-mass Planet

Authors: Suman Satyal, Billy Quarles, Marialis Rosario-Franco

Abstract: All 4 giant planets in the Solar System host systems of multiple moons, whereas the terrestrial planets only host up to 2 moons. The Earth can capture small asteroids as temporary satellites, which begs the question as to how many moons could stably orbit the Earth, or an Earth-mass exoplanet. We perform a series of N-body simulations of closely-spaced equal mass moons in nested orbits around an E… ▽ More All 4 giant planets in the Solar System host systems of multiple moons, whereas the terrestrial planets only host up to 2 moons. The Earth can capture small asteroids as temporary satellites, which begs the question as to how many moons could stably orbit the Earth, or an Earth-mass exoplanet. We perform a series of N-body simulations of closely-spaced equal mass moons in nested orbits around an Earth-mass planet orbiting a Sun-like star. The innermost moon begins near the host planets Roche radius, and the system is packed until the outermost moon begins near the stability limit for single moons. The initial spacing of the moons follows an iterative scheme commonly used for studies of compact planetary systems around single stars. For 3-moons system, we generate MEGNO maps to calculate periodic and chaotic regions and to identify the destabilizing MMRs. Our calculations show that the maximum number of moons depends on the assumed masses of the satellites (Ceres-, Pluto-, and Luna-mass) that could maintain stable orbits in a tightly-packed environment. Through our N-body simulations, we find stable configurations for up to 7 $\pm$ 1 Ceres-mass, 4 $\pm$ 1 Pluto-mass, and 3 $\pm$ 1 Luna-mass moons. However, outward tidal migration will likely play a substantial role in the number of moons on stable orbits over the 10 Gyr stellar lifetime of a Sun-like star. △ Less

Submitted 6 August, 2022; originally announced August 2022.

Comments: Published at MNRAS

Habitable planet formation around low-mass stars: Rapid accretion, rapid debris removal and the essential contribution of external giants

Authors: Matthew S. Clement, Elisa V. Quintana, Billy L. Quarles

Abstract: In recent years a paradigm shift has occurred in exoplanet science, wherein low-mass stars are increasingly viewed as a foundational pillar of the search for potentially habitable worlds in the solar neighborhood. However, the formation processes of this rapidly accumulating sample of planet systems are still poorly understood. Moreover, it is unclear whether tenuous primordial atmospheres around… ▽ More In recent years a paradigm shift has occurred in exoplanet science, wherein low-mass stars are increasingly viewed as a foundational pillar of the search for potentially habitable worlds in the solar neighborhood. However, the formation processes of this rapidly accumulating sample of planet systems are still poorly understood. Moreover, it is unclear whether tenuous primordial atmospheres around these Earth-analogs could have survived the intense epoch of heightened stellar activity that is typical for low-mass stars. We present new simulations of in-situ planet formation across the M-dwarf mass spectrum, and derive leftover debris populations of small bodies that might source delayed volatile delivery. We then follow the evolution of this debris with high-resolution models of real systems of habitable zone planets around low-mass stars such as TRAPPIST-1, Proxima Centauri and TOI-700. While debris in the radial vicinity of the habitable zone planets is removed rapidly, thus making delayed volatile delivery highly unlikely, we find that material ubiquitously scattered into an exo-asteroid belt region during the planet formation process represents a potentially lucrative reservoir of icy small bodies. Thus, the presence of external ~Neptune-Saturn mass planets capable of dynamically perturbing these asteroids would be a sign that habitable zone worlds around low-mass stars might have avoided complete desiccation. However, we also find that such giant planets significantly limit the efficiency of asteroidal implantation during the planet formation process. In the coming decade, long-baseline radial velocity studies and Roman Space Telescope microlensing observations will undoubtedly further constrain this process. △ Less

Submitted 11 February, 2022; originally announced February 2022.

Comments: 19 pages, 9 figures, 3 tables, accepted for publication in ApJ

An Exomoon Survey of 70 Cool Giant Exoplanets and the New Candidate Kepler-1708 b-i

Authors: David Kipping, Steve Bryson, Chris Burke, Jessie Christiansen, Kevin Hardegree-Ullman, Billy Quarles, Brad Hansen, Judit Szulágyi, Alex Teachey

Abstract: Exomoons represent a crucial missing puzzle piece in our efforts to understand extrasolar planetary systems. To address this deficiency, we here describe an exomoon survey of 70 cool, giant transiting exoplanet candidates found by Kepler. We identify only one which exhibits a moon-like signal that passes a battery of vetting tests: Kepler-1708 b. We show that Kepler-1708 b is a statistically valid… ▽ More Exomoons represent a crucial missing puzzle piece in our efforts to understand extrasolar planetary systems. To address this deficiency, we here describe an exomoon survey of 70 cool, giant transiting exoplanet candidates found by Kepler. We identify only one which exhibits a moon-like signal that passes a battery of vetting tests: Kepler-1708 b. We show that Kepler-1708 b is a statistically validated Jupiter-sized planet orbiting a Sun-like quiescent star at ~1.6AU. The signal of the exomoon candidate, Kepler-1708 b-i, is a 4.8-sigma effect and is persistent across different instrumental detrending methods, with a 1% false-positive probability via injection-recovery. Kepler-1708 b-i is ~2.6 Earth radii and is located in an approximately coplanar orbit at ~12 planetary radii from its ~1.6AU Jupiter-sized host. Future observations will be necessary to validate or reject the candidate. △ Less

Submitted 12 January, 2022; originally announced January 2022.

Comments: 115 pages, 24 figures, 4 tables. Published in Nature Astronomy. Posteriors available at https://doi.org/10.5061/dryad.18931zcz9

Updated Studies on Exomoons in the HD 23079 System

Authors: Oshina Jagtap, Billy Quarles, Manfred Cuntz

Abstract: We re-evaluate the outer edge of orbital stability for possible exomoons orbiting the radial velocity planet discovered in the HD 23079 system. In this system, a solar-type star hosts a Jupiter-mass planet in a nearly circular orbit in the outer stellar habitable zone. The outer stability limit of exomoons is deduced using $N$-body and tidal migration simulations considering a large range of initi… ▽ More We re-evaluate the outer edge of orbital stability for possible exomoons orbiting the radial velocity planet discovered in the HD 23079 system. In this system, a solar-type star hosts a Jupiter-mass planet in a nearly circular orbit in the outer stellar habitable zone. The outer stability limit of exomoons is deduced using $N$-body and tidal migration simulations considering a large range of initial conditions, encompassing both prograde and retrograde orbits. In particular, we extend previous works by evaluating many values in the satellite mean anomaly to identify and exclude regions of quasi-stability. Future observations of this system can make use of our results through a scale factor relative to the currently measured minimum mass. Using a constant time lag tidal model (Hut 1981), we find that plausible tidal interactions within the system are insufficient to induce significant outward migration toward the theoretical stability limit. While current technologies are incapable of detecting exomoons in this system, we comment on the detectability of putative moons through Doppler monitoring within direct imaging observations in view of future research capacities. △ Less

Submitted 19 October, 2021; originally announced October 2021.

Comments: 10 pages, 5 figures, 2 tables; accepted for publication in PASA

Exploring tidal obliquity variations with SMERCURY-T

Authors: Steven M. Kreyche, Jason W. Barnes, Billy L. Quarles, John E. Chambers

Abstract: We introduce our new code, SMERCURY-T, which is based on existing codes SMERCURY (Lissauer et al. 2012) and Mercury-T (Bolmont et al. 2015). The result is a mixed-variable symplectic N-body integrator that can compute the orbital and spin evolution of a planet within a multi-planet system under the influence of tidal spin torques from its star. We validate our implementation by comparing our exper… ▽ More We introduce our new code, SMERCURY-T, which is based on existing codes SMERCURY (Lissauer et al. 2012) and Mercury-T (Bolmont et al. 2015). The result is a mixed-variable symplectic N-body integrator that can compute the orbital and spin evolution of a planet within a multi-planet system under the influence of tidal spin torques from its star. We validate our implementation by comparing our experimental results to that of a secular model. As we demonstrate in a series of experiments, SMERCURY-T allows for the study of secular spin-orbit resonance crossings and captures for planets within complex multi-planet systems. These processes can drive a planet's spin state to evolve along vastly different pathways on its road toward tidal equilibrium, as tidal spin torques dampen the planet's spin rate and evolve its obliquity. Additionally, we show the results of a scenario that exemplifies the crossing of a chaotic region that exists as the overlap of two spin-orbit resonances. The test planet experiences violent and chaotic swings in its obliquity until its eventual escape from resonance as it tidally evolves. All of these processes are and have been important over the obliquity evolution of many bodies within the Solar System and beyond, and have implications for planetary climate and habitability. SMERCURY-T is a powerful and versatile tool that allows for further study of these phenomena. △ Less

Submitted 7 September, 2021; originally announced September 2021.

Comments: 21 pages, 9 figures, 3 tables; accepted by The Planetary Science Journal

Milankovitch Cycles for a Circumstellar Earth-analog within $α$ Centauri-like Binaries

Authors: Billy Quarles, Gongjie Li, Jack J. Lissauer

Abstract: An Earth-analog orbiting within the habitable zone of $α$ Centauri B was shown to undergo large variations in its obliquity, or axial tilt, which affects the planetary climate by altering the radiative flux for a given latitude. We examine the potential implications of these obliquity variations for climate through Milankovitch cycles using an energy balance model with ice sheets. Similar to previ… ▽ More An Earth-analog orbiting within the habitable zone of $α$ Centauri B was shown to undergo large variations in its obliquity, or axial tilt, which affects the planetary climate by altering the radiative flux for a given latitude. We examine the potential implications of these obliquity variations for climate through Milankovitch cycles using an energy balance model with ice sheets. Similar to previous studies, the largest amplitude obliquity variations from spin-orbit resonances induce snowball states within the habitable zone, while moderate variations can allow for persistent ice caps or an ice belt. Particular outcomes for the global ice distribution can depend on the planetary orbit, obliquity, spin precession, binary orbit, and which star the Earth-analog orbits. An Earth-analog with an inclined orbits relative to the binary orbit can periodically transition through several global ice distribution states and risk runaway glaciation when periods of ice caps and an ice belt overlap. When determining the potential habitability for planets in stellar binaries, more care must be taken due to the orbital and spin dynamics. △ Less

Submitted 20 November, 2021; v1 submitted 28 August, 2021; originally announced August 2021.

Comments: 25 pages, 16 figures, 1 table; fixed missing reference; accepted for publication in MNRAS; public repository at https://github.com/saturnaxis/Ice-ages-in-AlphaCen

Journal ref: Monthly Notices of the Royal Astronomical Society, 2021, stab3179

TIC 172900988: A Transiting Circumbinary Planet Detected in One Sector of TESS Data

Authors: Veselin B. Kostov, Brian P. Powell, Jerome A. Orosz, William F. Welsh, William Cochran, Karen A. Collins, Michael Endl, Coel Hellier, David W. Latham, Phillip MacQueen, Joshua Pepper, Billy Quarles, Lalitha Sairam, Guillermo Torres, Robert F. Wilson, Serge Bergeron, Pat Boyce, Allyson Bieryla, Robert Buchheim, Caleb Ben Christiansen, David R. Ciardi, Kevin I. Collins, Dennis M. Conti, Scott Dixon, Pere Guerra , et al. (64 additional authors not shown)

Abstract: We report the first discovery of a transiting circumbinary planet detected from a single sector of TESS data. During Sector 21, the planet TIC 172900988b transited the primary star and then 5 days later it transited the secondary star. The binary is itself eclipsing, with a period of P = 19.7 days and an eccentricity of e = 0.45. Archival data from ASAS-SN, Evryscope, KELT, and SuperWASP reveal a… ▽ More We report the first discovery of a transiting circumbinary planet detected from a single sector of TESS data. During Sector 21, the planet TIC 172900988b transited the primary star and then 5 days later it transited the secondary star. The binary is itself eclipsing, with a period of P = 19.7 days and an eccentricity of e = 0.45. Archival data from ASAS-SN, Evryscope, KELT, and SuperWASP reveal a prominent apsidal motion of the binary orbit, caused by the dynamical interactions between the binary and the planet. A comprehensive photodynamical analysis of the TESS, archival and follow-up data yields stellar masses and radii of M1 = 1.2384 +/- 0.0007 MSun and R1 = 1.3827 +/- 0.0016 RSun for the primary and M2 = 1.2019 +/- 0.0007 MSun and R2 = 1.3124 +/- 0.0012 RSun for the secondary. The radius of the planet is R3 = 11.25 +/- 0.44 REarth (1.004 +/- 0.039 RJup). The planet's mass and orbital properties are not uniquely determined - there are six solutions with nearly equal likelihood. Specifically, we find that the planet's mass is in the range of 824 < M3 < 981 MEarth (2.65 < M3 < 3.09 MJup), its orbital period could be 188.8, 190.4, 194.0, 199.0, 200.4, or 204.1 days, and the eccentricity is between 0.02 and 0.09. At a V = 10.141 mag, the system is accessible for high-resolution spectroscopic observations, e.g. Rossiter-McLaughlin effect and transit spectroscopy. △ Less

Submitted 27 August, 2021; v1 submitted 18 May, 2021; originally announced May 2021.

Comments: 57 pages, 30 figures, 25 tables; Accepted AJ

Exomoons in Systems with a Strong Perturber: Applications to $α$ Cen AB

Authors: Billy Quarles, Siegfried Eggl, Marialis Rosario-Franco, Gongjie Li

Abstract: The presence of a stellar companion can place constraints on occurrence and orbital evolution of satellites orbiting exoplanets, i.e., exomoons. In this work we revise earlier orbital stability limits for retrograde orbits in the case of a three body system consisting of star-planet-satellite. The latter reads $a_{\rm sat}^{\rm crit} \approx 0.668(1-1.236e_{\rm p})$ for $e_p \leq 0.8$ in units of… ▽ More The presence of a stellar companion can place constraints on occurrence and orbital evolution of satellites orbiting exoplanets, i.e., exomoons. In this work we revise earlier orbital stability limits for retrograde orbits in the case of a three body system consisting of star-planet-satellite. The latter reads $a_{\rm sat}^{\rm crit} \approx 0.668(1-1.236e_{\rm p})$ for $e_p \leq 0.8$ in units of the Hill Radius and represents the lower critical orbit as a function of the planetary eccentricity $e_{\rm p}$. A similar formula is determined for exomoons hosted by planets in binary star systems, where $e_{\rm p}$ is replaced with the components of free and forced eccentricity from secular orbit evolution theory. By exploring the dynamics of putative exomoons in $α$ Centauri AB we find that the outer stability limit can be much less than half the Hill Radius due to oscillations in the planetary orbital eccentricity caused by the gravitational interaction with the binary star. We show, furthermore, how the resulting truncation of the outer stability limit can affect the outward tidal migration and potential observability of exomoons through transit timing variations (TTVs). Typical TTV (RMS) amplitudes induced by exomoons in binary systems are $\lesssim$10 min and appear more likely for planets orbiting the less massive stellar component. A GitHub repository (saturnaxis/exomoon-in-binaries) is available to reproduce figures. △ Less

Submitted 23 May, 2021; v1 submitted 30 April, 2021; originally announced May 2021.

Comments: 16 pages, 11 figures, 3 tables; accepted in Astronomical Journal; for the GitHub repository, see https://github.com/saturnaxis/exomoon-in-binaries

Journal ref: AJ 162 58 (2021)

Warm Jupiters in TESS Full-Frame Images: A Catalog and Observed Eccentricity Distribution for Year 1

Authors: Jiayin Dong, Chelsea X. Huang, Rebekah I. Dawson, Daniel Foreman-Mackey, Karen A. Collins, Samuel N. Quinn, Jack J. Lissauer, Thomas G. Beatty, Billy Quarles, Lizhou Sha, Avi Shporer, Zhao Guo, Stephen R. Kane, Lyu Abe, Khalid Barkaoui, Zouhair Benkhaldoun, Rafael A. Brahm, Francois Bouchy, Theron W. Carmichael, Kevin I. Collins, Dennis M. Conti, Nicolas Crouzet, Georgina Dransfield, Phil Evans, Tianjun Gan , et al. (35 additional authors not shown)

Abstract: Warm Jupiters -- defined here as planets larger than 6 Earth radii with orbital periods of 8--200 days -- are a key missing piece in our understanding of how planetary systems form and evolve. It is currently debated whether Warm Jupiters form in situ, undergo disk or high eccentricity tidal migration, or have a mixture of origin channels. These different classes of origin channels lead to differe… ▽ More Warm Jupiters -- defined here as planets larger than 6 Earth radii with orbital periods of 8--200 days -- are a key missing piece in our understanding of how planetary systems form and evolve. It is currently debated whether Warm Jupiters form in situ, undergo disk or high eccentricity tidal migration, or have a mixture of origin channels. These different classes of origin channels lead to different expectations for Warm Jupiters' properties, which are currently difficult to evaluate due to the small sample size. We take advantage of the \TESS survey and systematically search for Warm Jupiter candidates around main-sequence host stars brighter than the \TESS-band magnitude of 12 in the Full-Frame Images in Year 1 of the \TESS Prime Mission data. We introduce a catalog of 55 Warm Jupiter candidates, including 19 candidates that were not originally released as \TESS Objects of Interest (TOIs) by the \TESS team. We fit their \TESS light curves, characterize their eccentricities and transit-timing variations (TTVs), and prioritize a list for ground-based follow-up and \TESS Extended Mission observations. Using hierarchical Bayesian modeling, we find the preliminary eccentricity distributions of our Warm-Jupiter-candidate catalog using a Beta distribution, a Rayleigh distribution, and a two-component Gaussian distribution as the functional forms of the eccentricity distribution. Additional follow-up observations will be required to clean the sample of false positives for a full statistical study, derive the orbital solutions to break the eccentricity degeneracy, and provide mass measurements. △ Less

Submitted 5 April, 2021; originally announced April 2021.

Comments: 30 pages, 8 figures, 6 tables. submitted to ApJS, revised in response to referee report

Application of Orbital Stability and Tidal Migration Constraints for Exomoon Candidates

Authors: Billy Quarles, Gongjie Li, Marialis Rosario-Franco

Abstract: Satellites of extrasolar planets, or exomoons, are on the frontier of detectability using current technologies and theoretical constraints should be considered in their search. In this Letter, we apply theoretical constraints of orbital stability and tidal migration to the six candidate KOI systems proposed by Fox & Wigert (2020) to identify whether these systems can potentially host exomoons. The… ▽ More Satellites of extrasolar planets, or exomoons, are on the frontier of detectability using current technologies and theoretical constraints should be considered in their search. In this Letter, we apply theoretical constraints of orbital stability and tidal migration to the six candidate KOI systems proposed by Fox & Wigert (2020) to identify whether these systems can potentially host exomoons. The host planets orbit close to their respective stars and the orbital stability extent of exomoons is limited to only $\sim$40% of the host planet's Hill radius ($\sim$20 R$_{\rm p}$). Using plausible tidal parameters from the solar system, we find that four out of six systems would either tidally disrupt their exomoons or lose them to outward migration within the system lifetimes. The remaining two systems (KOI 268.01 and KOI 1888.01) could host exomoons that are within 25 R$_{\rm p}$ and less than $\sim$3% of the host planet's mass. However, a recent independent transit timing analysis by Kipping (2020) found that these systems fail rigorous statistical tests to validate them as candidates. Overall, we find the presence of exomoons in these systems that are large enough for TTV signatures to be unlikely given the combined constraints of observational modeling, tidal migration, and orbital stability. Software to reproduce our results is available in the GitHub repository: Multiversario/satcand. △ Less

Submitted 30 September, 2020; originally announced September 2020.

Comments: 6 pages, 4 figures, 1 table; Accepted for publication in ApJL; for the GitHub repository, see https://github.com/Multiversario/satcand

Journal ref: ApJL 902 L20 (2020)

Effects of Flux Variation on the Surface Temperatures of Earth-like Circumbinary Planets

Authors: Srisurya Karthik Yadavalli, Billy Quarles, Gongjie Li, Nader Haghighipour

Abstract: The Kepler Space telescope has uncovered around thirteen circumbinary planets (CBPs) that orbit a pair of stars and experience two sources of stellar flux. We characterize the top-of-atmosphere flux and surface temperature evolution in relation to the orbital short-term dynamics between the central binary star and an Earth-like CBP analog. We compare the differential evolution of an Earth-like CBP… ▽ More The Kepler Space telescope has uncovered around thirteen circumbinary planets (CBPs) that orbit a pair of stars and experience two sources of stellar flux. We characterize the top-of-atmosphere flux and surface temperature evolution in relation to the orbital short-term dynamics between the central binary star and an Earth-like CBP analog. We compare the differential evolution of an Earth-like CBP's flux and surface temperature with that of an equivalent single-star (ESS) system to uncover the degree by which the potential habitability of the planet could vary. For a Sun-like primary, we find that the flux variation over a single planetary orbit is greatest when the dynamical mass ratio is ~0.3 for a G-K spectral binary. Using a latitudinal energy balance model, we show that the ice-albedo feedback plays a substantial role in Earth-like CBP habitability due to the interplay between flux redistribution (via obliquity) and changes in the total flux (via binary gyration). We examine the differential evolution of flux and surface temperature for Earth-like analogs of the habitable zone CBPs (4 Kepler and 1 hypothetical system) and find that these analogs are typically warmer than their ESS counterparts. △ Less

Submitted 21 September, 2020; v1 submitted 11 August, 2020; originally announced August 2020.

Comments: 12 pages, 14 figures, 5 tables; revised relative to first review at MNRAS, updated table 1 with Kopparapu et al. (2014) values

Multiple Transits during a Single Conjunction: Identifying Transiting Circumbinary Planetary Candidates from TESS

Authors: Veselin B. Kostov, William F. Welsh, Nader Haghighipour, Eric Agol, Daniel C. Fabrycky, Billy Quarles, Gongjie Li, Sean M. Mills, Laurance R. Doyle, Tsevi Mazeh, Jerome A. Orosz, David Martin, Brian Powell

Abstract: We present results of a study on identifying circumbinary planet candidates that produce multiple transits during one conjunction with eclipsing binary systems. The occurrence of these transits enables us to estimate the candidates' orbital periods, which is crucial as the periods of the currently known transiting circumbinary planets are significantly longer than the typical observational baselin… ▽ More We present results of a study on identifying circumbinary planet candidates that produce multiple transits during one conjunction with eclipsing binary systems. The occurrence of these transits enables us to estimate the candidates' orbital periods, which is crucial as the periods of the currently known transiting circumbinary planets are significantly longer than the typical observational baseline of TESS. Combined with the derived radii, it also provides valuable information needed for follow-up observations and subsequent confirmation of a large number of circumbinary planet candidates from TESS. Motivated by the discovery of the 1108-day circumbinary planet Kepler-1647, we show the application of this technique to four of Kepler's circumbinary planets that produce such transits. Our results indicate that in systems where the circumbinary planet is on a low-eccentricity orbit, the estimated planetary orbital period is within <10-20% of the true value. This estimate is derived from photometric observations spanning less than 5% of the planet's period, demonstrating the strong capability of the technique. Capitalizing on the current and future eclipsing binaries monitored by NASA's TESS mission, we estimate that hundreds of circumbinary planets candidates producing multiple transits during one conjunction will be detected in the TESS data. Such a large sample will enable statistical understanding of the population of planets orbiting binary stars and shed new light on their formation and evolution. △ Less

Submitted 6 August, 2020; originally announced August 2020.

Comments: 16 pages, 6 figures, 2 tables, AJ accepted

Journal ref: AJ 160 174 (2020)

Evolution of $α$ Centauri B's protoplanetary disc

Authors: Rebecca G. Martin, Jack J. Lissauer, Billy Quarles

Abstract: With hydrodynamical simulations we examine the evolution of a protoplanetary disc around $α$ Centauri B including the effect of the eccentric orbit binary companion $α$ Centauri A. The initially circular orbit disc undergoes two types of eccentricity growth. First, the eccentricity oscillates on the orbital period of the binary, $P_{\rm orb}$, due to the eccentricity of the binary orbit. Secondly,… ▽ More With hydrodynamical simulations we examine the evolution of a protoplanetary disc around $α$ Centauri B including the effect of the eccentric orbit binary companion $α$ Centauri A. The initially circular orbit disc undergoes two types of eccentricity growth. First, the eccentricity oscillates on the orbital period of the binary, $P_{\rm orb}$, due to the eccentricity of the binary orbit. Secondly, for a sufficiently small disc aspect ratio, the disc undergoes global forced eccentricity oscillations on a time-scale of around $20\,P_{\rm orb}$. These oscillations damp out through viscous dissipation leaving a quasi-steady eccentricity profile for the disc that oscillates only on the binary orbital period. The time-averaged global eccentricity is in the range 0.05-0.1, with no precession in the steady state. The periastrons of the gas particles are aligned to one another. The higher the disc viscosity, the higher the disc eccentricity. With $N$-body simulations we examine the evolution of a disc of planetesimals that forms with the orbital properties of the quasi-steady protoplanetary disc. We find that the average magnitude of the eccentricity of particles increases and their periastrons become misaligned to each other once they decouple from the gas disc. The low planetesimal collision velocity required for planet formation suggests that for planet formation to have occurred in a disc of planetesimals formed from a protoplanetary disc around $α$ Centauri B, said disc's viscosity must be have been small and planet formation must have occurred at orbital radii smaller than about $2.5\,\rm au$. Planet formation may be easier with the presence of gas. △ Less

Submitted 9 June, 2020; originally announced June 2020.

Comments: Accepted for publication in MNRAS

Orbital Stability of Exomoons and Submoons with Applications to Kepler 1625b-I

Authors: Marialis Rosario-Franco, Billy Quarles, Zdzislaw Musielak, Manfred Cuntz

Abstract: An intriguing question in the context of dynamics arises: Could a moon possess a moon itself? Such a configuration does not exist in the Solar System, although this may be possible in theory. Kollmeier et al. (2019) determined the critical size of a satellite necessary to host a long-lived sub-satellite, or submoon. However, the orbital constraints for these submoons to exist are still undetermine… ▽ More An intriguing question in the context of dynamics arises: Could a moon possess a moon itself? Such a configuration does not exist in the Solar System, although this may be possible in theory. Kollmeier et al. (2019) determined the critical size of a satellite necessary to host a long-lived sub-satellite, or submoon. However, the orbital constraints for these submoons to exist are still undetermined. Domingos et al. (2006) indicated that moons are stable out to a fraction of the host planet Hill radius $R_{H,p}$, which in turn depends on the eccentricity of its host's orbit. Motivated by this, we simulate a system of exomoons and submoons for $10^5$ planetary orbits, while considering many initial orbital phases to obtain the critical semimajor axis in terms of $R_{H,p}$ or the hosts satellite's Hill radius $R_{H,sat}$, respectively. We find that, assuming circular coplanar orbits, the stability limit for exomoons is 0.40 $R_{H,p}$ and for a submoon is 0.33 $R_{H,sat}$. Additionally, we discuss the observational feasibility of detecting these sub-satellites through photometric, radial velocity, or direct imaging observations using the Neptunes-sized exomoon candidates Kepler 1625b-I (Teachey et al. 2018) and identify how stability can shape the identification of future candidates. △ Less

Submitted 13 May, 2020; originally announced May 2020.

Comments: Accepted to AJ. 15 pages, 5 figures, 2 tables

TOI-1338: TESS' First Transiting Circumbinary Planet

Authors: Veselin B. Kostov, Jerome A. Orosz, Adina D. Feinstein, William F. Welsh, Wolf Cukier, Nader Haghighipour, Billy Quarles, David V. Martin, Benjamin T. Montet, Guillermo Torres, Amaury H. M. J. Triaud, Thomas Barclay, Patricia Boyd, Cesar Briceno, Andrew Collier Cameron, Alexandre C. M. Correia, Emily A. Gilbert, Samuel Gill, Michael Gillon, Jacob Haqq-Misra, Coel Hellier, Courtney Dressing, Daniel C. Fabrycky, Gabor Furesz, Jon Jenkins , et al. (43 additional authors not shown)

Abstract: We report the detection of the first circumbinary planet found by TESS. The target, a known eclipsing binary, was observed in sectors 1 through 12 at 30-minute cadence and in sectors 4 through 12 at two-minute cadence. It consists of two stars with masses of 1.1 MSun and 0.3 MSun on a slightly eccentric (0.16), 14.6-day orbit, producing prominent primary eclipses and shallow secondary eclipses. Th… ▽ More We report the detection of the first circumbinary planet found by TESS. The target, a known eclipsing binary, was observed in sectors 1 through 12 at 30-minute cadence and in sectors 4 through 12 at two-minute cadence. It consists of two stars with masses of 1.1 MSun and 0.3 MSun on a slightly eccentric (0.16), 14.6-day orbit, producing prominent primary eclipses and shallow secondary eclipses. The planet has a radius of ~6.9 REarth and was observed to make three transits across the primary star of roughly equal depths (~0.2%) but different durations -- a common signature of transiting circumbinary planets. Its orbit is nearly circular (e ~ 0.09) with an orbital period of 95.2 days. The orbital planes of the binary and the planet are aligned to within ~1 degree. To obtain a complete solution for the system, we combined the TESS photometry with existing ground-based radial-velocity observations in a numerical photometric-dynamical model. The system demonstrates the discovery potential of TESS for circumbinary planets, and provides further understanding of the formation and evolution of planets orbiting close binary stars. △ Less

Submitted 16 April, 2020; originally announced April 2020.

Comments: 35 pages, 21 figures, 6 tables

Retrograde-rotating exoplanets experience obliquity excitations in an eccentricity-enabled resonance

Authors: Steven M. Kreyche, Jason W. Barnes, Billy L. Quarles, Jack J. Lissauer, John E. Chambers, Matthew M. Hedman

Abstract: Previous studies have shown that planets that rotate retrograde (backwards with respect to their orbital motion) generally experience less severe obliquity variations than those that rotate prograde (the same direction as their orbital motion). Here we examine retrograde-rotating planets on eccentric orbits and find a previously unknown secular spin-orbit resonance that can drive significant obliq… ▽ More Previous studies have shown that planets that rotate retrograde (backwards with respect to their orbital motion) generally experience less severe obliquity variations than those that rotate prograde (the same direction as their orbital motion). Here we examine retrograde-rotating planets on eccentric orbits and find a previously unknown secular spin-orbit resonance that can drive significant obliquity variations. This resonance occurs when the frequency of the planet's rotation axis precession becomes commensurate with an orbital eigenfrequency of the planetary system. The planet's eccentricity enables a participating orbital frequency through an interaction in which the apsidal precession of the planet's orbit causes a cyclic nutation of the planet's orbital angular momentum vector. The resulting orbital frequency follows the relationship $f = 2 \dot{\varpi} - \dotΩ$, where $\dot{\varpi}$ and $\dotΩ$ are the rates of the planet's changing longitude of periapsis and ascending node, respectively. We test this mechanism by simulating cases of a simple Earth-Jupiter system, and confirm the predicted resonance. Over the course of 100 Myr, the test Earths with rotation axis precession rates near the predicted resonant frequency experienced pronounced obliquity variations of order $10^\circ$-$30^\circ$. These variations can be significant, and suggest that while retrograde rotation is a stabilizing influence most of the time, retrograde rotators can experience large obliquity variations if they are on eccentric orbits and enter this spin-orbit resonance. △ Less

Submitted 7 April, 2020; v1 submitted 30 March, 2020; originally announced March 2020.

Comments: Accepted to The Planetary Science Journal. Figures 2 and 4b have associated video files hosted at https://github.com/SMKreyche/Kreyche-et-al.-2020-Videos

Kepler-1661 b: A Neptune-sized Kepler Transiting Circumbinary Planet around a Grazing Eclipsing Binary

Authors: Quentin J Socia, William F Welsh, Jerome A Orosz, William D Cochran, Michael Endl, Billy Quarles, Donald R Short, Guillermo Torres, Gur Windmiller, Mitchell Yenawine

Abstract: We report the discovery of a Neptune-size (R_p = 3.87 +/- 0.06 R_Earth) transiting circumbinary planet, Kepler-1661 b, found in the Kepler photometry. The planet has a period of ~175 days and its orbit precesses with a period of only 35 years. The precession causes the alignment of the orbital planes to vary, and the planet is in a transiting configuration only ~7% of the time as seen from Earth.… ▽ More We report the discovery of a Neptune-size (R_p = 3.87 +/- 0.06 R_Earth) transiting circumbinary planet, Kepler-1661 b, found in the Kepler photometry. The planet has a period of ~175 days and its orbit precesses with a period of only 35 years. The precession causes the alignment of the orbital planes to vary, and the planet is in a transiting configuration only ~7% of the time as seen from Earth. As with several other Kepler circumbinary planets, Kepler-1661 b orbits close to the stability radius, and is near the (hot) edge of habitable zone. The planet orbits a single-lined, grazing eclipsing binary, containing a 0.84 M_Sun and 0.26 M_Sun pair of stars in a mildly eccentric (e=0.11), 28.2-day orbit. The system is fairly young, with an estimated age of ~1-3 Gyrs, and exhibits significant starspot modulations. The grazing-eclipse configuration means the system is very sensitive to changes in the binary inclination, which manifests itself as a change in the eclipse depth. The starspots contaminate the eclipse photometry, but not in the usual way of inducing spurious eclipse timing variations. Rather, the starspots alter the normalization of the light curve, and hence the eclipse depths. This can lead to spurious eclipse depth variations, which are then incorrectly ascribed to binary orbital precession. △ Less

Submitted 9 January, 2020; originally announced January 2020.

Comments: 37 pages, 15 figures

Orbital Stability of Circumstellar Planets in Binary Systems

Authors: Billy Quarles, Gongjie Li, Veselin Kostov, Nader Haghighipour

Abstract: Planets that orbit only one of the stars in stellar binary systems (i.e., circumstellar) are dynamically constrained to a limited range of orbital parameters and thus understanding conditions on their stability is of great importance in exoplanet searches. We perform $\sim$700 million N-body simulations to identify how stability regions depend on properties of the binary, as well as, the starting… ▽ More Planets that orbit only one of the stars in stellar binary systems (i.e., circumstellar) are dynamically constrained to a limited range of orbital parameters and thus understanding conditions on their stability is of great importance in exoplanet searches. We perform $\sim$700 million N-body simulations to identify how stability regions depend on properties of the binary, as well as, the starting planetary inclination and mean longitude relative to the binary orbit. Moreover, we provide grid interpolation maps and lookup tables for the community to use our results. Through Monte-Carlo methods we determine that planets with a semimajor axis $a_p$ $\lesssim$8\% of the binary semimajor axis $a_{bin}$ will likely be stable given the known distribution of binary star parameters. This estimate varies in the Lidov-Kozai regime or for retrograde orbits to 4\% or 10\% of $a_{bin}$, respectively. Our method to quickly determine the circumstellar stability limit is important for interpreting observations of binaries using direct imaging with JWST, photometry with TESS, or even astrometry with Gaia. △ Less

Submitted 23 December, 2019; originally announced December 2019.

Comments: 28 pages, 10 figures, 6 tables; accepted for publication in the Astronomical Journal. Numerical tools and data are available in a Github (see https://github.com/saturnaxis/ThreeBody_Stability) and Zenodo (see http://doi.org/10.5281/zenodo.3579202) repository, respectively

Journal ref: AJ 159 3 (2020)

Obliquity Evolution of Circumstellar Planets in Sun-like Stellar Binaries

Authors: Billy Quarles, Gongjie Li, Jack J. Lissauer

Abstract: Changes in planetary obliquity, or axial tilt, influence the climates on Earth-like planets. In the solar system, the Earth's obliquity is stabilized due to interactions with our moon and the resulting {small amplitude variations ($\sim$2.4\degree)} are beneficial for advanced life. Most Sun-like stars have at least one stellar companion and the habitability of circumstellar exoplanets is shaped b… ▽ More Changes in planetary obliquity, or axial tilt, influence the climates on Earth-like planets. In the solar system, the Earth's obliquity is stabilized due to interactions with our moon and the resulting {small amplitude variations ($\sim$2.4\degree)} are beneficial for advanced life. Most Sun-like stars have at least one stellar companion and the habitability of circumstellar exoplanets is shaped by their stellar companion. We show that a stellar companion can dramatically change whether {Earth-like obliquity stability is} possible through planetary orbital precession relative to the binary orbit or resonant pumping of the obliquity through spin-orbit interactions. We present a new formalism for the planetary spin precession that accounts for orbital misalignments between the planet and binary. Using numerical modeling in $α$ Centauri AB we show: a stark contrast between the planetary obliquity variations depending on the host star, planetary neighbors limit the possible spin states for {Earth-like obliquity stability}, and the presence of a moon can destabilize the obliquity, defying our Earth-based expectations. An Earth-like rotator orbiting the primary star would experience {small} obliquity variations for 87\%, 74\%, or 54\% of Solar type binaries, depending on the mass of the primary (0.8, 1.0, or 1.2 M$_\odot$, respectively). Thus, Earth-like planets likely experience much larger obliquity variations, with more extreme climates, unless they are in specific states, such as orbiting nearly planar with the binary and rotating retrograde (backwards) like Venus. △ Less

Submitted 30 April, 2021; v1 submitted 19 November, 2019; originally announced November 2019.

Comments: 19 pages, 14 figures, 2 tables; Published in ApJ; corrected typo in Eq. 4

Journal ref: ApJ 886 56 (2019)

Searching for Planets Orbiting Alpha Centauri A with the James Webb Space Telescope

Authors: Charles Beichman, Marie Ygouf, Jorge Llop Sayson, Dimitri Mawet, Yuk Yung, Elodie Choquet, Pierre Kervella, Anthony Boccaletti, Ruslan Belikov, Jack J. Lissauer, Billy Quarles, Pierre-Olivier Lagage, Daniel Dicken, Renyu Hu, Bertrand Mennesson, Mike Ressler, Eugene Serabyn, John Krist, Eduardo Bendek, Jarron Leisenring, Laurent Pueyo

Abstract: Alpha Centauri A is the closest solar-type star to the Sun and offers an excellent opportunity to detect the thermal emission of a mature planet heated by its host star. The MIRI coronagraph on JWST can search the 1-3 AU (1"-2") region around alpha Cen A which is predicted to be stable within the alpha Cen AB system. We demonstrate that with reasonable performance of the telescope and instrument,… ▽ More Alpha Centauri A is the closest solar-type star to the Sun and offers an excellent opportunity to detect the thermal emission of a mature planet heated by its host star. The MIRI coronagraph on JWST can search the 1-3 AU (1"-2") region around alpha Cen A which is predicted to be stable within the alpha Cen AB system. We demonstrate that with reasonable performance of the telescope and instrument, a 20 hr program combining on-target and reference star observations at 15.5 um could detect thermal emission from planets as small as ~5 RE. Multiple visits every 3-6 months would increase the geometrical completeness, provide astrometric confirmation of detected sources, and push the radius limit down to ~3 RE. An exozodiacal cloud only a few times brighter than our own should also be detectable, although a sufficiently bright cloud might obscure any planet present in the system. While current precision radial velocity (PRV) observations set a limit of 50-100 ME at 1-3 AU for planets orbiting alpha Cen A, there is a broad range of exoplanet radii up to 10 RE consistent with these mass limits. A carefully planned observing sequence along with state-of-the-art post-processing analysis could reject the light from alpha Cen A at the level of ~10^-5 at 1"-2" and minimize the influence of alpha Cen B located 7-8" away in the 2022-2023 timeframe. These space-based observations would complement on-going imaging experiments at shorter wavelengths as well as PRV and astrometric experiments to detect planets dynamically. Planetary demographics suggest that the likelihood of directly imaging a planet whose mass and orbit are consistent with present PRV limits is small, ~5%, and possibly lower if the presence of a binary companion further reduces occurrence rates. However, at a distance of just 1.34 pc, alpha Cen A is our closest sibling star and certainly merits close scrutiny. △ Less

Submitted 21 October, 2019; originally announced October 2019.

Comments: PASP in press

Could there be an undetected inner planet near the stability limit in Kepler-1647?

Authors: Ziqian Hong, Billy Quarles, Gongjie Li, Jerome A. Orosz

Abstract: Kepler-1647b is the most recently discovered planet that transits two stars, i.e., a circumbinary planet (CBP). Due to its large orbital separation, Kepler-1647b stands out from the rest of the Kepler CBPs, which mostly reside on much tighter orbits near the stability limit. The large separation of Kepler-1647b challenges inward disk migration as a dominant formation pathway, suggested by the othe… ▽ More Kepler-1647b is the most recently discovered planet that transits two stars, i.e., a circumbinary planet (CBP). Due to its large orbital separation, Kepler-1647b stands out from the rest of the Kepler CBPs, which mostly reside on much tighter orbits near the stability limit. The large separation of Kepler-1647b challenges inward disk migration as a dominant formation pathway, suggested by the other Kepler CBPs. In this paper, we consider the possibility of an undetected planet near the stability limit by examining observational consequences of such a planet. We calculate the transit probability of the putative planet, transit timing variations (TTVs) of the known planet, and eclipsing timing variations (ETVs) of the host binary caused by the putative planet. We find the presence of a $\gtrsim$30M$_{\oplus}$ inner planet to be highly unlikely near the stability limit. In addition, we provide future TTV observation windows, which will further constrain possible undetected planets with lower masses. △ Less

Submitted 17 May, 2019; v1 submitted 16 April, 2019; originally announced April 2019.

Comments: 11 pages, 6 figures. Accepted for publication in AJ

Discovery of a Third Transiting Planet in the Kepler-47 Circumbinary System

Authors: Jerome A. Orosz, William F. Welsh, Nader Haghighipour, Billy Quarles, Donald R. Short, Sean M. Mills, Suman Satyal, Guillermo Torres, Eric Agol, Daniel C. Fabrycky, Daniel Jontof-Hutter, Gur Windmiller, Tobias W. A. Müller, Tobias C. Hinse, William D. Cochran, Michael Endl, Eric B. Ford, Tsevi Mazeh, Jack J. Lissauer

Abstract: Of the nine confirmed transiting circumbinary planet systems, only Kepler-47 is known to contain more than one planet. Kepler-47 b (the "inner planet") has an orbital period of 49.5 days and a radius of about $3\,R_{\oplus}$. Kepler-47 c (the "outer planet") has an orbital period of 303.2 days and a radius of about $4.7\,R_{\oplus}$. Here we report the discovery of a third planet, Kepler-47 d (the… ▽ More Of the nine confirmed transiting circumbinary planet systems, only Kepler-47 is known to contain more than one planet. Kepler-47 b (the "inner planet") has an orbital period of 49.5 days and a radius of about $3\,R_{\oplus}$. Kepler-47 c (the "outer planet") has an orbital period of 303.2 days and a radius of about $4.7\,R_{\oplus}$. Here we report the discovery of a third planet, Kepler-47 d (the "middle planet"), which has an orbital period of 187.4 days and a radius of about $7\,R_{\oplus}$. The presence of the middle planet allows us to place much better constraints on the masses of all three planets, where the $1σ$ ranges are less than $26\,M_{\oplus}$, between $7-43\,M_{\oplus}$, and between $2-5\,M_{\oplus}$ for the inner, middle, and outer planets, respectively. The middle and outer planets have low bulk densities, with $ρ_{\rm middle} < 0.68$ g cm$^{-3}$ and $ρ_{\rm outer} < 0.26$ g cm$^{-3}$ at the $1σ$ level. The two outer planets are "tightly packed," assuming the nominal masses, meaning no other planet could stably orbit between them. All of the orbits have low eccentricities and are nearly coplanar, disfavoring violent scattering scenarios and suggesting gentle migration in the protoplanetary disk. △ Less

Submitted 15 April, 2019; originally announced April 2019.

Comments: 68 pages, 30 figures

Journal ref: AJ, 157, 174 (2019)

Instabilities in the Early Solar System due to a Self-gravitating Disk

Authors: Billy Quarles, Nathan Kaib

Abstract: Modern studies of the early solar system routinely invoke the possibility of an orbital instability among the giant planets triggered by gravitational interactions between the planets and a massive exterior disk of planetesimals. Previous works have suggested that this instability can be substantially delayed ($\sim$100s Myr) after the formation of the giant planets. Bodies in the disk are typical… ▽ More Modern studies of the early solar system routinely invoke the possibility of an orbital instability among the giant planets triggered by gravitational interactions between the planets and a massive exterior disk of planetesimals. Previous works have suggested that this instability can be substantially delayed ($\sim$100s Myr) after the formation of the giant planets. Bodies in the disk are typically treated in a semi-active manner, wherein their gravitational force on the planets is included, but interactions between the planetesimals are ignored. We perform $N$-body numerical simulations using GENGA, which makes use of GPUs to allow for the inclusion of all gravitational interactions between bodies. Although our simulated Kuiper belt particles are more massive than the probable masses of real primordial Kuiper belt objects, our simulations indicate that the self-stirring of the primordial Kuiper belt is very important to the dynamics of the giant planet instability. We find that interactions between planetesimals dynamically heat the disk and typically prevent the outer solar system instability from being delayed by more than a few tens of million years after giant planet formation. Longer delays occur in a small fraction of systems that have at least 3.5 AU gaps between the planets and planetesimal disk. Our final planetary configurations match the solar system at a rate consistent with other previous works in most regards. Pre-instability heating of the disk typically yields final Jovian eccentricities comparable to the modern solar system value, which has been a difficult constraint to match in past works. △ Less

Submitted 20 December, 2018; originally announced December 2018.

Comments: 16 pages, 17 figures, 3 tables; accepted for publication in AJ

Stability Limits of Circumbinary Planets: Is There a Pile-up in the Kepler CBPs?

Authors: Billy Quarles, Suman Satyal, Veselin Kostov, Nathan Kaib, Nader Haghighipour

Abstract: The stability limit for circumbinary planets (CBPs) is not well defined and can depend on initial parameters defining either the planetary orbit or the inner binary orbit. We expand on the work of Holman & Wiegert (1999, AJ 117, 621) to develop numerical tools for quick, easy, and accurate determination of the stability limit. The results of our simulations, as well as our numerical tools, are ava… ▽ More The stability limit for circumbinary planets (CBPs) is not well defined and can depend on initial parameters defining either the planetary orbit or the inner binary orbit. We expand on the work of Holman & Wiegert (1999, AJ 117, 621) to develop numerical tools for quick, easy, and accurate determination of the stability limit. The results of our simulations, as well as our numerical tools, are available to the community through $\texttt{Zenodo}$ and $\texttt{GitHub}$, respectively. We employ a grid interpolation method based on $\sim$150 million full N-body simulations of initially circular, coplanar systems and compare to the 9 known Kepler CBP systems. Using a formalism from planet packing studies, we find that 55% of the Kepler CBP systems allow for an additional equal-mass planet to potentially exist on an interior orbit relative to the observed planet. Therefore, we do $\textit{not}$ find strong evidence for a pile-up in the Kepler CBP systems and more detections are needed to adequately characterize the formation mechanisms for the CBP population. Observations from the Transiting Exoplanet Survey Satellite are expected to substantially increase the number of detections using the unique geometry of CBP systems, where multiple transits can occur during a single conjunction. △ Less

Submitted 24 February, 2018; originally announced February 2018.

Comments: 18 pages, 8 figures, 4 tables; Accepted in Astrophysical Journal

Journal ref: B. Quarles et al 2018 ApJ 856 150

The Habitable Zone of Kepler-16: Impact of Binarity and Climate Models

Authors: S. Y. Moorman, B. L. Quarles, Zh. Wang, M. Cuntz

Abstract: We continue to investigate the binary system Kepler-16, consisting of a K-type main-sequence star, a red dwarf, and a circumbinary Saturnian planet. As part of our study, we describe the system's habitable zone based on different climate models. We also report on stability investigations for possible Earth-mass Trojans while expanding a previous study by B. L. Quarles and collaborators given in 20… ▽ More We continue to investigate the binary system Kepler-16, consisting of a K-type main-sequence star, a red dwarf, and a circumbinary Saturnian planet. As part of our study, we describe the system's habitable zone based on different climate models. We also report on stability investigations for possible Earth-mass Trojans while expanding a previous study by B. L. Quarles and collaborators given in 2012. For the climate models we carefully consider the relevance of the system's parameters. Furthermore, we pursue new stability simulations for the Earth-mass objects starting along the orbit of Kepler-16b. The eccentricity distribution as obtained prefers values close to circular, whereas the inclination distribution remains flat. The stable solutions are distributed near the co-orbital Lagrangian points, thus enhancing the plausibility that Earth-mass Trojans might be able to exist in the Kepler-16(AB) system. △ Less

Submitted 21 February, 2018; v1 submitted 21 January, 2018; originally announced February 2018.

Comments: Accepted by: International Journal of Astrobiology; 34 pages, 11 figures, 8 tables

Journal ref: International Journal of Astrobiology 18 (2019) 79-89

Long-Term Stability of Tightly Packed Multi-Planet Systems in Prograde, Coplanar, Circumstellar Orbits within the $α$ Centauri AB System

Authors: Billy Quarles, Jack J. Lissauer

Abstract: We perform long-term simulations, up to ten billion years, of closely-spaced configurations of 2 -- 6 planets, each as massive as the Earth, traveling on nested orbits about either stellar component in $α$ Centauri AB. The innermost planet initially orbits at either the inner edge of its star's empirical habitable zone (HZ) or the inner edge of its star's conservative HZ. Although individual plane… ▽ More We perform long-term simulations, up to ten billion years, of closely-spaced configurations of 2 -- 6 planets, each as massive as the Earth, traveling on nested orbits about either stellar component in $α$ Centauri AB. The innermost planet initially orbits at either the inner edge of its star's empirical habitable zone (HZ) or the inner edge of its star's conservative HZ. Although individual planets on low inclination, low eccentricity, orbits can survive throughout the habitable zones of both stars, perturbations from the companion star require that the minimum spacing of planets in multi-planet systems within the habitable zones of each star must be significantly larger than the spacing of similar multi-planet systems orbiting single stars in order to be long-lived. The binary companion induces a forced eccentricity upon the orbits of planets in orbit around either star. Planets on appropriately-phased circumstellar orbits with initial eccentricities equal to their forced eccentricities can survive on more closely spaced orbits than those with initially circular orbits, although the required spacing remains higher than for planets orbiting single stars. A total of up to nine planets on nested prograde orbits can survive for the current age of the system within the empirical HZs of the two stars, with five of these orbiting $α$ Centauri B and four orbiting $α$ Centauri A. △ Less

Submitted 18 January, 2018; originally announced January 2018.

Comments: 24 pages, 15 figures, 6 tables; accepted for publication in Astronomical Journal

Journal ref: 2018 AJ 155 130

Long-Term Stability of Planets in the $α$ Centauri System, II: Forced Eccentricities

Authors: Billy Quarles, Jack J. Lissauer, Nathan Kaib

Abstract: We extend our study of the extent of the regions within the $α$ Centauri AB star system where small planets are able to orbit for billion-year timescales (Quarles & Lissauer 2016, AJ 151, 111) to investigate the effects of minimizing the forced eccentricity of initial trajectories. We find that initially prograde, circumstellar orbits require a piecewise quadratic function to accurately approximat… ▽ More We extend our study of the extent of the regions within the $α$ Centauri AB star system where small planets are able to orbit for billion-year timescales (Quarles & Lissauer 2016, AJ 151, 111) to investigate the effects of minimizing the forced eccentricity of initial trajectories. We find that initially prograde, circumstellar orbits require a piecewise quadratic function to accurately approximate forced eccentricity as a function of semimajor axis, but retrograde orbits can be modeled using a linear function. Circumbinary orbits in the $α$ Centauri AB system are less affected by the forced eccentricity. {Planets on circumstellar orbits that begin with eccentricity vectors near their forced values are generally stable, up to $\sim$$10^9$ yr, out to a larger semimajor axis than are planets beginning on circular orbits. The amount by which the region of stability expands is much larger for retrograde orbits than it is for prograde orbits. The location of the stability boundary for two planet systems on prograde, circular orbits is much more sensitive to the initial eccentricity state than it is for analogous single planet systems. △ Less

Submitted 18 January, 2018; originally announced January 2018.

Comments: 11 pages, 10 figures, 3 tables; Accepted in Astronomical Journal

Journal ref: B. Quarles, J. J. Lissauer, and N. Kaib 2018 AJ 155 64

Obliquity Evolution of the Potentially Habitable Exoplanet Kepler-62f

Authors: Billy Quarles, Jason W. Barnes, Jack J. Lissauer, John Chambers

Abstract: Variations in the axial tilt, or obliquity, of terrestrial planets can affect their climates and therefore their habitability. Kepler-62f is a 1.4 R$_\oplus$ planet orbiting within the habitable zone of its K2 dwarf host star (Borucki et al. 2013). We perform N-body simulations that monitor the evolution of obliquity of Kepler-62f for 10 million year timescales to explore the effects on model assu… ▽ More Variations in the axial tilt, or obliquity, of terrestrial planets can affect their climates and therefore their habitability. Kepler-62f is a 1.4 R$_\oplus$ planet orbiting within the habitable zone of its K2 dwarf host star (Borucki et al. 2013). We perform N-body simulations that monitor the evolution of obliquity of Kepler-62f for 10 million year timescales to explore the effects on model assumptions, such as the masses of the Kepler-62 planets and the possibility of outer bodies. Significant obliquity variation occurs when the rotational precession frequency overlaps with one or more of the secular orbital frequencies, but most variations are limited to $\lesssim$10$^\circ$. Moderate variations ($\sim$10$^\circ - 20^\circ$) can occur over a broader range of initial obliquities when the relative nodal longitude ($ΔΩ$) overlaps with the frequency and phase of a given secular mode. However, we find that adding outer gas giants on long period orbits ($\gtrsim$ 1000 days) can produce large ($\sim$60$^\circ$) variations in obliquity if Kepler-62f has a very rapid (4 hr) rotation period. The possibility of giant planets on long period orbits impacts the climate and habitability of Kepler-62f through variations in the latitudinal surface flux, where the timescale for large variation can occur on million year timescales. △ Less

Submitted 23 July, 2019; v1 submitted 22 October, 2017; originally announced October 2017.

Comments: 15 pages, 17 figures, 5 tables; Revised relative to Referee Reports. Submitted to Astrobiology

Journal ref: Astrobiology, 20, 73 (2020)

Plausible Compositions of the Seven TRAPPIST-1 Planets Using Long-term Dynamical Simulations

Authors: Billy Quarles, Elisa V. Quintana, Eric D. Lopez, Joshua E. Schlieder, Thomas Barclay

Abstract: TRAPPIST-1 is a nearby ultra-cool dwarf that is host to a remarkable planetary system consisting of seven transiting planets. The orbital properties and radii of the planets have been well-constrained, and recently the masses of the inner six planets have been measured with additional ground and space-based photometric observations. Large uncertainties in these mass measurements have prevented a r… ▽ More TRAPPIST-1 is a nearby ultra-cool dwarf that is host to a remarkable planetary system consisting of seven transiting planets. The orbital properties and radii of the planets have been well-constrained, and recently the masses of the inner six planets have been measured with additional ground and space-based photometric observations. Large uncertainties in these mass measurements have prevented a robust analysis of the planetary compositions. Here we perform many thousands of N-body dynamical simulations with planet properties perturbed from the observed values and identify those that are stable for millions of years. This allows us to identify self-consistent orbital solutions that can be used in future studies. From our range of dynamical masses, we find that most of the planets are consistent with an Earth-like composition, where TRAPPIST-1f is likely to have a volatile-rich envelope. △ Less

Submitted 9 June, 2017; v1 submitted 7 April, 2017; originally announced April 2017.

Comments: 4 pages, 2 figure, 3 tables, accepted to ApJL; Changes made due to first review including modification and addition of figures

Long-Term Stability of Planets in the $α$ Centauri System

Authors: Billy Quarles, Jack J. Lissauer

Abstract: We evaluate the extent of the regions within the $α$ Centauri AB star system where small planets are able to orbit for billion-year timescales, and we calculate the positions on the sky plane where planets on stable orbits about either stellar component may appear. We confirm the qualitative results of Wiegert and Holman (AJ 113, 1445, 1997) regarding the approximate size of the regions of stable… ▽ More We evaluate the extent of the regions within the $α$ Centauri AB star system where small planets are able to orbit for billion-year timescales, and we calculate the positions on the sky plane where planets on stable orbits about either stellar component may appear. We confirm the qualitative results of Wiegert and Holman (AJ 113, 1445, 1997) regarding the approximate size of the regions of stable orbits, which are larger for retrograde orbits relative to the binary than for prograde orbits. Additionally, we find that mean motion resonances with the binary orbit leave an imprint on the limits of orbital stability, and the effects of the Lidov-Kozai mechanism are also readily apparent. △ Less

Submitted 17 April, 2016; originally announced April 2016.

Comments: 6 pages, 11 figures, 2 tables, Published in AJ

Journal ref: 2016 The Astronomical Journal 151 111

Obliquity Variability of a Potentially Habitable Early Venus

Authors: Jason W. Barnes, Billy Quarles, Jack J. Lissauer, John Chambers, Matthew M. Hedman

Abstract: Venus currently rotates slowly, with its spin controlled by solid-body and atmospheric thermal tides. However, conditions may have been far different 4 billion years ago, when the Sun was fainter and most of the carbon within Venus could have been in solid form, implying a low-mass atmosphere. We investigate how the obliquity would have varied for a hypothetical rapidly rotating Early Venus. The o… ▽ More Venus currently rotates slowly, with its spin controlled by solid-body and atmospheric thermal tides. However, conditions may have been far different 4 billion years ago, when the Sun was fainter and most of the carbon within Venus could have been in solid form, implying a low-mass atmosphere. We investigate how the obliquity would have varied for a hypothetical rapidly rotating Early Venus. The obliquity variation structure of an ensemble of hypothetical Early Venuses is simpler than that Earth would have if it lacked its large Moon (Lissauer et al., 2012), having just one primary chaotic regime at high prograde obliquities. We note an unexpected long-term variability of up to $\pm7^\circ$ for retrograde Venuses. Low-obliquity Venuses show very low total obliquity variability over billion-year timescales -- comparable to that of the real Moon-influenced Earth. △ Less

Submitted 29 January, 2016; originally announced February 2016.

Comments: in revision at _Astrobiology_

Kepler Eclipsing Binary Stars. VII. The Catalog of Eclipsing Binaries Found in the Entire Kepler Data-Set

Authors: Brian Kirk, Kyle Conroy, Andrej Prša, Michael Abdul-Masih, Angela Kochoska, Gal Matijevič, Kelly Hambleton, Thomas Barclay, Steven Bloemen, Tabetha Boyajian, Laurance R. Doyle, B. J. Fulton, Abe Johannes Hoekstra, Kian Jek, Stephen R. Kane, Veselin Kostov, David Latham, Tsevi Mazeh, Jerome A. Orosz, Joshua Pepper, Billy Quarles, Darin Ragozzine, Avi Shporer, John Southworth, Keivan Stassun , et al. (25 additional authors not shown)

Abstract: The primary Kepler Mission provided nearly continuous monitoring of ~200,000 objects with unprecedented photometric precision. We present the final catalog of eclipsing binary systems within the 105 square degree Kepler field of view. This release incorporates the full extent of the data from the primary mission (Q0-Q17 Data Release). As a result, new systems have been added, additional false posi… ▽ More The primary Kepler Mission provided nearly continuous monitoring of ~200,000 objects with unprecedented photometric precision. We present the final catalog of eclipsing binary systems within the 105 square degree Kepler field of view. This release incorporates the full extent of the data from the primary mission (Q0-Q17 Data Release). As a result, new systems have been added, additional false positives have been removed, ephemerides and principal parameters have been recomputed, classifications have been revised to rely on analytical models, and eclipse timing variations have been computed for each system. We identify several classes of systems including those that exhibit tertiary eclipse events, systems that show clear evidence of additional bodies, heartbeat systems, systems with changing eclipse depths, and systems exhibiting only one eclipse event over the duration of the mission. We have updated the period and galactic latitude distribution diagrams and included a catalog completeness evaluation. The total number of identified eclipsing and ellipsoidal binary systems in the Kepler field of view has increased to 2878, 1.3% of all observed Kepler targets. An online version of this catalog with downloadable content and visualization tools is maintained at http://keplerEBs.villanova.edu. △ Less

Submitted 29 December, 2015; originally announced December 2015.

Comments: 52 pages, 14 figures, aastex

Planetary Candidates Observed by Kepler. VII. The First Fully Uniform Catalog Based on The Entire 48 Month Dataset (Q1-Q17 DR24)

Authors: Jeffrey L. Coughlin, F. Mullally, Susan E. Thompson, Jason F. Rowe, Christopher J. Burke, David W. Latham, Natalie M. Batalha, Aviv Ofir, Billy L. Quarles, Christopher E. Henze, Angie Wolfgang, Douglas A. Caldwell, Stephen T. Bryson, Avi Shporer, Joseph Catanzarite, Rachel Akeson, Thomas Barclay, William J. Borucki, Tabetha S. Boyajian, Jennifer R. Campbell, Jessie L. Christiansen, Forrest R. Girouard, Michael R. Haas, Steve B. Howell, Daniel Huber , et al. (10 additional authors not shown)

Abstract: We present the seventh Kepler planet candidate catalog, which is the first to be based on the entire, uniformly processed, 48 month Kepler dataset. This is the first fully automated catalog, employing robotic vetting procedures to uniformly evaluate every periodic signal detected by the Q1-Q17 Data Release 24 (DR24) Kepler pipeline. While we prioritize uniform vetting over the absolute correctness… ▽ More We present the seventh Kepler planet candidate catalog, which is the first to be based on the entire, uniformly processed, 48 month Kepler dataset. This is the first fully automated catalog, employing robotic vetting procedures to uniformly evaluate every periodic signal detected by the Q1-Q17 Data Release 24 (DR24) Kepler pipeline. While we prioritize uniform vetting over the absolute correctness of individual objects, we find that our robotic vetting is overall comparable to, and in most cases is superior to, the human vetting procedures employed by past catalogs. This catalog is the first to utilize artificial transit injection to evaluate the performance of our vetting procedures and quantify potential biases, which are essential for accurate computation of planetary occurrence rates. With respect to the cumulative Kepler Object of Interest (KOI) catalog, we designate 1,478 new KOIs, of which 402 are dispositioned as planet candidates (PCs). Also, 237 KOIs dispositioned as false positives (FPs) in previous Kepler catalogs have their disposition changed to PC and 118 PCs have their disposition changed to FP. This brings the total number of known KOIs to 8,826 and PCs to 4,696. We compare the Q1-Q17 DR24 KOI catalog to previous KOI catalogs, as well as ancillary Kepler catalogs, finding good agreement between them. We highlight new PCs that are both potentially rocky and potentially in the habitable zone of their host stars, many of which orbit solar-type stars. This work represents significant progress in accurately determining the fraction of Earth-size planets in the habitable zone of Sun-like stars. The full catalog is publicly available at the NASA Exoplanet Archive. △ Less

Submitted 18 February, 2016; v1 submitted 18 December, 2015; originally announced December 2015.

Comments: Accepted to the Astrophysical Journal Supplement Series. 30 pages, 9 figures, 7 tables. We make the DR24 robovetter decision code publicly available at http://github.com/JeffLCoughlin/robovetter, with input and output examples provided using the same data as contained in the full paper's tables

Kepler-1647b: the largest and longest-period Kepler transiting circumbinary planet

Authors: Veselin B. Kostov, Jerome A. Orosz, William F. Welsh, Laurance R. Doyle, Daniel C. Fabrycky, Nader Haghighipour, Billy Quarles, Donald R. Short, William D. Cochran, Michael Endl, Eric B. Ford, Joao Gregorio, Tobias C. Hinse, Howard Isaacson, Jon M. Jenkins, Eric L. N. Jensen, Stephen Kane, Ilya Kull, David W. Latham, Jack J. Lissauer, Geoffrey W. Marcy, Tsevi Mazeh, Tobias W. A. Muller, Joshua Pepper, Samuel N. Quinn , et al. (6 additional authors not shown)

Abstract: We report the discovery of a new Kepler transiting circumbinary planet (CBP). This latest addition to the still-small family of CBPs defies the current trend of known short-period planets orbiting near the stability limit of binary stars. Unlike the previous discoveries, the planet revolving around the eclipsing binary system Kepler-1647 has a very long orbital period (~1100 days) and was at conju… ▽ More We report the discovery of a new Kepler transiting circumbinary planet (CBP). This latest addition to the still-small family of CBPs defies the current trend of known short-period planets orbiting near the stability limit of binary stars. Unlike the previous discoveries, the planet revolving around the eclipsing binary system Kepler-1647 has a very long orbital period (~1100 days) and was at conjunction only twice during the Kepler mission lifetime. Due to the singular configuration of the system, Kepler-1647b is not only the longest-period transiting CBP at the time of writing, but also one of the longest-period transiting planets. With a radius of 1.06+/-0.01 RJup it is also the largest CBP to date. The planet produced three transits in the light-curve of Kepler-1647 (one of them during an eclipse, creating a syzygy) and measurably perturbed the times of the stellar eclipses, allowing us to measure its mass to be 1.52+/-0.65 MJup. The planet revolves around an 11-day period eclipsing binary consisting of two Solar-mass stars on a slightly inclined, mildly eccentric (e_bin = 0.16), spin-synchronized orbit. Despite having an orbital period three times longer than Earth's, Kepler-1647b is in the conservative habitable zone of the binary star throughout its orbit. △ Less

Submitted 19 May, 2016; v1 submitted 1 December, 2015; originally announced December 2015.

Comments: 61 pages, 21 figures

The three-body problem

Authors: Z. E. Musielak, B. Quarles

Abstract: The three-body problem, which describes three masses interacting through Newtonian gravity without any restrictions imposed on the initial positions and velocities of these masses, has attracted the attention of many scientists for more than 300 years. In this paper, we present a review of the three-body problem in the context of both historical and modern developments. We describe the general and… ▽ More The three-body problem, which describes three masses interacting through Newtonian gravity without any restrictions imposed on the initial positions and velocities of these masses, has attracted the attention of many scientists for more than 300 years. In this paper, we present a review of the three-body problem in the context of both historical and modern developments. We describe the general and restricted (circular and elliptic) three-body problems, different analytical and numerical methods of finding solutions, methods for performing stability analysis, search for periodic orbits and resonances, and application of the results to some interesting astronomical and space dynamical settings. We also provide a brief presentation of the general and restricted relativistic three-body problems, and discuss their astronomical applications. △ Less

Submitted 10 August, 2015; originally announced August 2015.

Comments: 49 pages, 10 figures, Published in Reports on Progress in Physics

Journal ref: Reports on Progress in Physics, Volume 77, Issue 6, article id. 065901 (2014)

Planetary Candidates Observed by Kepler VI: Planet Sample from Q1-Q16 (47 Months)

Authors: F. Mullally, Jeffrey L. Coughlin, Susan E. Thompson, Jason Rowe, Christopher Burke, David W. Latham, Natalie M. Batalha, Stephen T. Bryson, Jessie Christiansen, Christopher E. Henze, Aviv Ofir, Billy Quarles, Avi Shporer, Vincent Van Eylen, Christa Van Laerhoven, Yash Shah, Angie Wolfgang, W. J. Chaplin, Ji-Wei Xie, Rachel Akeson, Vic Argabright, Eric Bachtell, Thomas Barclay William J. Borucki, Douglas A. Caldwell, Jennifer R. Campbell , et al. (34 additional authors not shown)

Abstract: \We present the sixth catalog of Kepler candidate planets based on nearly 4 years of high precision photometry. This catalog builds on the legacy of previous catalogs released by the Kepler project and includes 1493 new Kepler Objects of Interest (KOIs) of which 554 are planet candidates, and 131 of these candidates have best fit radii <1.5 R_earth. This brings the total number of KOIs and planet… ▽ More \We present the sixth catalog of Kepler candidate planets based on nearly 4 years of high precision photometry. This catalog builds on the legacy of previous catalogs released by the Kepler project and includes 1493 new Kepler Objects of Interest (KOIs) of which 554 are planet candidates, and 131 of these candidates have best fit radii <1.5 R_earth. This brings the total number of KOIs and planet candidates to 7305 and 4173 respectively. We suspect that many of these new candidates at the low signal-to-noise limit may be false alarms created by instrumental noise, and discuss our efforts to identify such objects. We re-evaluate all previously published KOIs with orbital periods of >50 days to provide a consistently vetted sample that can be used to improve planet occurrence rate calculations. We discuss the performance of our planet detection algorithms, and the consistency of our vetting products. The full catalog is publicly available at the NASA Exoplanet Archive. △ Less

Submitted 6 February, 2015; originally announced February 2015.

Comments: 18 pages, to be published in the Astrophysical Journal Supplement Series

Planetary Candidates Observed by Kepler V: Planet Sample from Q1-Q12 (36 Months)

Authors: Jason F. Rowe, Jeffrey L. Coughlin, Victoria Antoci, Thomas Barclay, Natalie M. Batalha, William J. Borucki, Christopher J. Burke, Steven T. Bryson, Douglas A. Caldwell, Jennifer R. Campbell, Joseph H. Catanzarite, Jessie L. Christiansen, William Cochran, Ronald L. Gilliland, Forrest R. Girouard, Michael R. Haas, Krzysztof G. Helminiak, Christopher E. Henze, Kelsey L. Hoffman, Steve B. Howell, Daniel Huber, Roger C. Hunter, Hannah Jang-Condell, Jon M. Jenkins, Todd C. Klaus , et al. (21 additional authors not shown)

Abstract: The Kepler mission discovered 2842 exoplanet candidates with 2 years of data. We provide updates to the Kepler planet candidate sample based upon 3 years (Q1-Q12) of data. Through a series of tests to exclude false-positives, primarily caused by eclipsing binary stars and instrumental systematics, 855 additional planetary candidates have been discovered, bringing the total number known to 3697. We… ▽ More The Kepler mission discovered 2842 exoplanet candidates with 2 years of data. We provide updates to the Kepler planet candidate sample based upon 3 years (Q1-Q12) of data. Through a series of tests to exclude false-positives, primarily caused by eclipsing binary stars and instrumental systematics, 855 additional planetary candidates have been discovered, bringing the total number known to 3697. We provide revised transit parameters and accompanying posterior distributions based on a Markov Chain Monte Carlo algorithm for the cumulative catalogue of Kepler Objects of Interest. There are now 130 candidates in the cumulative catalogue that receive less than twice the flux the Earth receives and more than 1100 have a radius less than 1.5 Rearth. There are now a dozen candidates meeting both criteria, roughly doubling the number of candidate Earth analogs. A majority of planetary candidates have a high probability of being bonafide planets, however, there are populations of likely false-positives. We discuss and suggest additional cuts that can be easily applied to the catalogue to produce a set of planetary candidates with good fidelity. The full catalogue is publicly available at the NASA Exoplanet Archive. △ Less

Submitted 29 January, 2015; v1 submitted 28 January, 2015; originally announced January 2015.

Comments: Accepted for publication, ApJS

Dynamical Evolution of the Earth-Moon Progenitors - Whence Theia?

Authors: Billy Quarles, Jack J. Lissauer

Abstract: We present integrations of a model Solar System with five terrestrial planets (beginning ~30-50 Myr after the formation of primitive Solar System bodies) in order to determine the preferred regions of parameter space leading to a giant impact that resulted in the formation of the Moon. Our results indicate which choices of semimajor axes and eccentricities for Theia (the proto-Moon) at this epoch… ▽ More We present integrations of a model Solar System with five terrestrial planets (beginning ~30-50 Myr after the formation of primitive Solar System bodies) in order to determine the preferred regions of parameter space leading to a giant impact that resulted in the formation of the Moon. Our results indicate which choices of semimajor axes and eccentricities for Theia (the proto-Moon) at this epoch can produce a late Giant Impact, assuming that Mercury, Venus, and Mars are near the current orbits. We find that the likely semimajor axis of Theia, at the epoch when our simulations begin, depends on the assumed mass ratio of Earth-Moon progenitors (8/1, 4/1, or 1/1). The low eccentricities of the terrestrial planets are most commonly produced when the progenitors have similar semimajor axes at the epoch when our integrations commence. Additionally, we show that mean motion resonances among the terrestrial planets and perturbations from the giant planets can affect the dynamical evolution of the system leading to a late Giant Impact. △ Less

Submitted 27 October, 2014; originally announced October 2014.

Comments: 17 pages, Accepted for publication in Icarus

KIC 9632895 - The 10th Kepler Transiting Circumbinary Planet

Authors: William F. Welsh, Jerome A. Orosz, Donald R. Short, Nader Haghighipour, Lars A. Buchhave, Laurance R. Doyle, Daniel C. Fabrycky, Tobias Cornelius Hinse, Stephen Kane, Veselin Kostov, Tsevi Mazeh, Sean M. Mills, Tobias W. A. Mueller, Billy Quarles, Samuel N. Quinn, Darin Ragozzine, Avi Shporer, Jason H. Steffen, Lev Tal-Or, Guillermo Torres, Gur Windmiller, William J. Borucki

Abstract: We present the discovery of KIC 9632895b, a 6.2 Earth-radius planet in a low-eccentricity, 240.5-day orbit about an eclipsing binary. The binary itself consists of a 0.93 and 0.194 solar mass pair of stars with an orbital period of 27.3 days. The plane of the planet's orbit is rapidly precessing, and its inclination only becomes sufficiently aligned with the primary star in the latter portion of t… ▽ More We present the discovery of KIC 9632895b, a 6.2 Earth-radius planet in a low-eccentricity, 240.5-day orbit about an eclipsing binary. The binary itself consists of a 0.93 and 0.194 solar mass pair of stars with an orbital period of 27.3 days. The plane of the planet's orbit is rapidly precessing, and its inclination only becomes sufficiently aligned with the primary star in the latter portion of the Kepler data. Thus three transits are present in the latter half of the light curve, but none of the three conjunctions that occurred during the first half of the light curve produced transits. The precession period is ~103 years, and during that cycle, transits are visible only ~8% of the time. This has the important implication that for every system like KIC 9632895 that we detect, there are ~12 circumbinary systems that exist but are not currently exhibiting transits. The planet's mass is too small to noticeably perturb the binary, consequently its mass is not measurable with these data; but our photodynamical model places a 1-sigma upper limit of 16 Earth masses. With a period 8.8 times that of the binary, the planet is well outside the dynamical instability zone. It does, however, lie within the habitable zone of the binary, and making it the third of ten Kepler circumbinary planets to do so. △ Less

Submitted 4 September, 2014; originally announced September 2014.

Comments: 13 figures; 5 tables; submitted to ApJ

Kepler Eclipsing Binary Stars. V. Identification of 31 Eclipsing Binaries in the K2 Engineering Data-set

Authors: Kyle E. Conroy, Andrej Prša, Keivan G. Stassun, Steven Bloemen, Mahmoud Parvizi, Billy Quarles, Tabetha Boyajian, Thomas Barclay, Avi Shporer, David W. Latham, Michael Abdul-Masih

Abstract: Over 2500 eclipsing binaries were identified and characterized from the ultra-precise photometric data provided by the Kepler space telescope. Kepler is now beginning its second mission, K2, which is proving to again provide ultra-precise photometry for a large sample of eclipsing binary stars. In the 1951 light curves covering 12 days in the K2 engineering data-set, we have identified and determi… ▽ More Over 2500 eclipsing binaries were identified and characterized from the ultra-precise photometric data provided by the Kepler space telescope. Kepler is now beginning its second mission, K2, which is proving to again provide ultra-precise photometry for a large sample of eclipsing binary stars. In the 1951 light curves covering 12 days in the K2 engineering data-set, we have identified and determined the ephemerides for 31 eclipsing binaries that demonstrate the capabilities for eclipsing binary science in the upcoming campaigns in K2. Of those, 20 are new discoveries. We describe both manual and automated approaches to harvesting the complete set of eclipsing binaries in the K2 data, provide identifications and details for the full set of eclipsing binaries present in the engineering data-set, and discuss the prospects for application of eclipsing binary searches in the K2 mission. △ Less

Submitted 14 July, 2014; originally announced July 2014.

Comments: 12 pages, 2 figures, submitted to PASP

An Earth-sized Planet in the Habitable Zone of a Cool Star

Authors: Elisa V. Quintana, Thomas Barclay, Sean N. Raymond, Jason F. Rowe, Emeline Bolmont, Douglas A. Caldwell, Steve B. Howell, Stephen R. Kane, Daniel Huber, Justin R. Crepp, Jack J. Lissauer, David R. Ciardi, Jeffrey L. Coughlin, Mark E. Everett, Christopher E. Henze, Elliott Horch, Howard Isaacson, Eric B. Ford, Fred C. Adams, Martin Still, Roger C. Hunter, Billy Quarles, Franck Selsis

Abstract: The quest for Earth-like planets represents a major focus of current exoplanet research. While planets that are Earth-sized and smaller have been detected, these planets reside in orbits that are too close to their host star to allow liquid water on their surface. We present the detection of Kepler-186f, a 1.11+\-0.14 Earth radius planet that is the outermost of five planets - all roughly Earth-si… ▽ More The quest for Earth-like planets represents a major focus of current exoplanet research. While planets that are Earth-sized and smaller have been detected, these planets reside in orbits that are too close to their host star to allow liquid water on their surface. We present the detection of Kepler-186f, a 1.11+\-0.14 Earth radius planet that is the outermost of five planets - all roughly Earth-sized - that transit a 0.47+\-0.05 Rsun star. The intensity and spectrum of the star's radiation places Kepler-186f in the stellar habitable zone, implying that if Kepler-186f has an Earth-like atmosphere and H2O at its surface, then some of this H2O is likely to be in liquid form. △ Less

Submitted 22 April, 2014; originally announced April 2014.

Comments: Main text and supplemental material combined. This is the author's version of the work. It is posted here by permission of the AAAS for personal use, not for redistribution. The definitive version was published in Science on 18 April 2014, Vol. 344, #6181. Published version available from http://www.sciencemag.org/content/344/6181/277

Chaotic dynamics of the planet in HD 196885 AB

Authors: Suman Satyal, Tobias Cornelius Hinse, Billy Quarles, Joaquin Noyola

Abstract: Depending on the planetary orbit around the host star(s), a planet could orbit either one or both stars in a binary system as S-type or P-type, respectively. We have analysed the dynamics of the S-type planetary system in HD 196885 AB with an emphasis on a planet with a higher orbital inclination relative to the binary plane. The mean exponential growth factor of nearby orbits (MEGNO) maps are use… ▽ More Depending on the planetary orbit around the host star(s), a planet could orbit either one or both stars in a binary system as S-type or P-type, respectively. We have analysed the dynamics of the S-type planetary system in HD 196885 AB with an emphasis on a planet with a higher orbital inclination relative to the binary plane. The mean exponential growth factor of nearby orbits (MEGNO) maps are used as an indicator to determine regions of periodicity and chaos for the various choices of the planet's semimajor axis, eccentricity and inclination with respect to the previously determined observational uncertainties. We have quantitatively mapped out the chaotic and quasi-periodic regions of the system's phase space which indicate a likely regime of the planet's inclination. In addition, we inspect the resonant angle to determine whether alternation between libration and circulation occurs as a consequence of Kozai oscillations, a probable mechanism that can drive the planetary orbit to a very large inclination. Also, we demonstrate the possible higher mass limit of the planet and improve upon the current dynamical model based on our analysis. △ Less

Submitted 18 June, 2014; v1 submitted 6 January, 2014; originally announced January 2014.

Comments: 10 pages, 9 figures (Accepted for publication at MNRAS)

On the Possibility of Habitable Moons in the System of HD 23079: Results from Orbital Stability Studies

Authors: M. Cuntz, B. Quarles, J. Eberle, A. Shukayr

Abstract: The aim of our study is to investigate the possibility of habitable moons orbiting the giant planet HD 23079b, a Jupiter-mass planet, which follows a low-eccentricity orbit in the outer region of HD 23079's habitable zone. We show that HD 23079b is able to host habitable moons in prograde and retrograde orbits, as expected, noting that the outer stability limit for retrograde orbits is increased b… ▽ More The aim of our study is to investigate the possibility of habitable moons orbiting the giant planet HD 23079b, a Jupiter-mass planet, which follows a low-eccentricity orbit in the outer region of HD 23079's habitable zone. We show that HD 23079b is able to host habitable moons in prograde and retrograde orbits, as expected, noting that the outer stability limit for retrograde orbits is increased by nearly 90% compared to that of prograde orbits, a result consistent with previous generalized studies. For the targeted parameter space it was found that the outer stability limit for habitable moons varies between 0.05236 and 0.06955 AU (prograde orbits) and between 0.1023 and 0.1190 AU (retrograde orbits) depending on the orbital parameters of the Jupiter-type planet if a minimum mass is assumed. These intervals correspond to 0.306 and 0.345 (prograde orbits) and 0.583 and 0.611 (retrograde orbits) of the planet's Hill radius. Larger stability limits are obtained if an increased value for the planetary mass m_p is considered; they are consistent with the theoretically deduced relationship of m_p^{1/3}. Finally, we compare our results to the statistical formulae of Domingos et al. (2006) [MNRAS 373, 1227], indicating both concurrence and limitations. △ Less

Submitted 8 April, 2013; v1 submitted 3 April, 2013; originally announced April 2013.

Comments: 15 pages, 8 figures, 6 tables; Publ. Astron. Soc. Australia (PASA); in press. arXiv admin note: text overlap with arXiv:1104.3092

Application of chaos indicators in the study of dynamics of S-type extrasolar planets in stellar binaries

Authors: S. Satyal, B. Quarles, T. C. Hinse

Abstract: The orbits of two individual planets in two known binary star systems, γCephei and HD 196885 are numerically integrated using various numerical techniques to assess the chaotic or quasi-periodic nature of the dynamical system considered. The Hill stability (HS) function which measures the orbital perturbation of a planet around the primary star due to the secondary star is calculated for each syst… ▽ More The orbits of two individual planets in two known binary star systems, γCephei and HD 196885 are numerically integrated using various numerical techniques to assess the chaotic or quasi-periodic nature of the dynamical system considered. The Hill stability (HS) function which measures the orbital perturbation of a planet around the primary star due to the secondary star is calculated for each system. The maximum Lyapunov exponent (MLE) time series are generated to measure the divergence/convergence rate of stable manifolds, which are used to differentiate between chaotic and non-chaotic orbits. Then, we calculate dynamical Mean Exponential Growth factor of Nearby Orbits (MEGNO) maps from solving the variational equations along with the equations of motion. These maps allow us to accurately differentiate between stable and unstable dynamical systems. The results obtained from the analysis of HS, MLE, and MEGNO maps are analysed for their dynamical variations and resemblance. The HS test for the planets shows stability and quasi-periodicity for at least ten million years. The MLE and the MEGNO maps have also indicated the local quasi-periodicity and global stability in a relatively short integration period. The orbital stability of the systems is tested using each indicator for various values of planet inclinations (i_{pl} \le 25^\circ) and binary eccentricities. The reliability of HS criterion is also discussed based on its stability results compared with the MLE and MEGNO maps. △ Less

Submitted 17 May, 2013; v1 submitted 16 November, 2012; originally announced November 2012.

Comments: 12 pages, 12 figures, 2 tables; accepted for publication in MNRAS

Study of Resonances for the Restricted 3-Body Problem

Authors: Billy Quarles, Zdzislaw Musielak, Manfred Cuntz

Abstract: Our aim is to identify and classify mean-motion resonances (MMRs) for the coplanar circular restricted three-body problem (CR3BP) for mass ratios between 0.10 and 0.50. Our methods include the maximum Lyapunov exponent, which is used as an indicator for the location of the resonances, the Fast Fourier Transform (FFT) used for determining what kind of resonances are present, and the inspection of t… ▽ More Our aim is to identify and classify mean-motion resonances (MMRs) for the coplanar circular restricted three-body problem (CR3BP) for mass ratios between 0.10 and 0.50. Our methods include the maximum Lyapunov exponent, which is used as an indicator for the location of the resonances, the Fast Fourier Transform (FFT) used for determining what kind of resonances are present, and the inspection of the orbital elements to classify the periodicity. We show that the 2:1 resonance occurs the most frequently. Among other resonances, the 3:1 resonance is the second most common, and furthermore both 3:2 and 5:3 resonances occur more often than the 4:1 resonance. Moreover, the resonances in the coplanar CR3BP are classified based on the behaviour of the orbits. We show that orbital stability is ensured for high values of resonance (i.e., high ratios) where only a single resonance is present. The resonances attained are consistent with the previously established resonances for the Solar System, i.e., specifically, in regards to the asteroid belt. Previous work employed digital filtering and Lyapunov characteristic exponents to determine stochasticity of the eccentricity, which is found to be consistent with our usage of Lyapunov exponents as an alternate approach based on varying the mass ratio instead of the eccentricity. Our results are expected to be of principal interest to future studies, including augmentations to observed or proposed resonances, of extra-solar planets in binary stellar systems. △ Less

Submitted 4 May, 2012; originally announced May 2012.

Comments: 10 pages, 9 figures, 5 tables; accepted by Astronomical Notes (AN)