Showing 1–34 of 34 results for author: Hadden, S

Action-Angle Variables for Axisymmetric Potentials via Birkhoff Normalization

Authors: Sam Hadden

Abstract: We describe a method for calculating action-angle variables in axisymmetric galactic potentials using Birkhoff normalization, a technique from Hamiltonian perturbation theory. An advantageous feature of this method is that it yields explicit series expressions for both the forward and inverse transformations between the action-angle variables and position-velocity data. It also provides explicit e… ▽ More We describe a method for calculating action-angle variables in axisymmetric galactic potentials using Birkhoff normalization, a technique from Hamiltonian perturbation theory. An advantageous feature of this method is that it yields explicit series expressions for both the forward and inverse transformations between the action-angle variables and position-velocity data. It also provides explicit expressions for the Hamiltonian and dynamical frequencies as functions of the action variables. We test this method by examining orbits in a Miyamoto-Nagai model potential and compare it to the popular Stäckel approximation method. When vertical actions are not too large, the Birkhoff normalization method achieves fractional errors smaller than a part in $10^{3}$ and outperforms the Stäckel approximation. We also show that the range over which Birkhoff normalization provides accurate results can be extended by constructing Padé approximants from the perturbative series expressions developed with the method. Numerical routines in Python for carrying out the Birkhoff normalization procedure are made available. △ Less

Submitted 25 April, 2024; originally announced April 2024.

Comments: Submitted to AAS journals. Code available at https://github.com/shadden/AA_variables_via_Birkhoff_Normalization

The instability mechanism of compact multiplanet systems

Authors: Caleb Lammers, Sam Hadden, Norman Murray

Abstract: To improve our understanding of orbital instabilities in compact planetary systems, we compare suites of $N$-body simulations against numerical integrations of simplified dynamical models. We show that, surprisingly, dynamical models that account for small sets of resonant interactions between the planets can accurately recover $N$-body instability times. This points towards a simple physical pict… ▽ More To improve our understanding of orbital instabilities in compact planetary systems, we compare suites of $N$-body simulations against numerical integrations of simplified dynamical models. We show that, surprisingly, dynamical models that account for small sets of resonant interactions between the planets can accurately recover $N$-body instability times. This points towards a simple physical picture in which a handful of three-body resonances (3BRs), generated by interactions between nearby two-body mean motion resonances (MMRs), overlap and drive chaotic diffusion, leading to instability. Motivated by this, we show that instability times are well-described by a power law relating instability time to planet separations, measured in units of fractional semi-major axis difference divided by the planet-to-star mass ratio to the $1/4$ power, rather than the frequently-adopted $1/3$ power implied by measuring separations in units of mutual Hill radii. For idealized systems, the parameters of this power-law relationship depend only on the ratio of the planets' orbital eccentricities to the orbit-crossing value, and we report an empirical fit to enable quick instability time predictions. This relationship predicts that observed systems comprised of three or more sub-Neptune mass planets must be spaced with period ratios $\mathcal{P} \gtrsim 1.35$ and that tightly-spaced systems ($\mathcal{P} \lesssim 1.5$) must possess very low eccentricities ($e \lesssim 0.05$) to be stable for more than $10^9$ orbits. △ Less

Submitted 26 March, 2024; originally announced March 2024.

Comments: 20 pages, 13 figures. Submitted to AAS Journals

Masses, Revised Radii, and a Third Planet Candidate in the "Inverted" Planetary System Around TOI-1266

Authors: Ryan Cloutier, Michael Greklek-McKeon, Serena Wurmser, Collin Cherubim, Erik Gillis, Andrew Vanderburg, Sam Hadden, Charles Cadieux, Étienne Artigau, Shreyas Vissapragada, Annelies Mortier, Mercedes López-Morales, David W. Latham, Heather Knutson, Raphaëlle D. Haywood, Enric Pallé, René Doyon, Neil Cook, Gloria Andreuzzi, Massimo Cecconi, Rosario Cosentino, Adriano Ghedina, Avet Harutyunyan, Matteo Pinamonti, Manu Stalport , et al. (18 additional authors not shown)

Abstract: Is the population of close-in planets orbiting M dwarfs sculpted by thermally driven escape or is it a direct outcome of the planet formation process? A number of recent empirical results strongly suggest the latter. However, the unique architecture of the TOI-1266 system presents a challenge to models of planet formation and atmospheric escape given its seemingly "inverted" architecture of a larg… ▽ More Is the population of close-in planets orbiting M dwarfs sculpted by thermally driven escape or is it a direct outcome of the planet formation process? A number of recent empirical results strongly suggest the latter. However, the unique architecture of the TOI-1266 system presents a challenge to models of planet formation and atmospheric escape given its seemingly "inverted" architecture of a large sub-Neptune ($P_b=10.9$ days, $R_{p,b}=2.62\pm 0.11\, \mathrm{R}_{\oplus}$) orbiting interior to that of the system's smaller planet ($P_c=18.8$ days, $R_{p,c}=2.13\pm 0.12\, \mathrm{R}_{\oplus}$). Here we present revised planetary radii based on new TESS and diffuser-assisted ground-based transit observations, and characterize both planetary masses using a set of 145 radial velocity measurements from HARPS-N ($M_{p,b}=4.23\pm 0.69\, \mathrm{M}_{\oplus}, M_{p,c}=2.88\pm 0.80\, \mathrm{M}_{\oplus}$). Our analysis also reveals a third planet candidate ($P_d=32.3$ days, $M_{p,d}\sin{i} = 4.59^{+0.96}_{-0.94}\, \mathrm{M}_{\oplus}$), which if real, would form a chain of near 5:3 period ratios, although the system is likely not in a mean motion resonance. Our results indicate that TOI-1266 b and c are among the lowest density sub-Neptunes around M dwarfs and likely exhibit distinct bulk compositions of a gas-enveloped terrestrial ($X_{\mathrm{env},b}=5.5\pm 0.7$%) and a water-rich world (WMF$_c=59\pm 14$%), which is supported by hydrodynamic escape models. If distinct bulk compositions are confirmed through atmospheric characterization, the system's unique architecture would represent an interesting test case of inside-out sub-Neptune formation at pebble traps. △ Less

Submitted 6 November, 2023; v1 submitted 20 October, 2023; originally announced October 2023.

Comments: Accepted for publication in MNRAS. 21 pages. Our spectroscopic time series are included in the arXiv source files as table6.csv

Scattered Disk Dynamics: The Mapping Approach

Authors: Sam Hadden, Scott Tremaine

Abstract: We derive, and discuss the properties of, a symplectic map for the dynamics of bodies on nearly parabolic orbits. The orbits are perturbed by a planet on a circular, coplanar orbit interior to the pericenter of the parabolic orbit. The map shows excellent agreement with direct numerical integrations and elucidates how the dynamics depends on perturber mass and pericenter distance. We also use the… ▽ More We derive, and discuss the properties of, a symplectic map for the dynamics of bodies on nearly parabolic orbits. The orbits are perturbed by a planet on a circular, coplanar orbit interior to the pericenter of the parabolic orbit. The map shows excellent agreement with direct numerical integrations and elucidates how the dynamics depends on perturber mass and pericenter distance. We also use the map to explore the onset of chaos, statistical descriptions of chaotic transport, and sticking in mean-motion resonances. We discuss implications of our mapping model for the dynamical evolution of the solar system's scattered disk and other highly eccentric trans-Neptunian objects. △ Less

Submitted 1 September, 2023; originally announced September 2023.

Comments: submitted to MNRAS

Mercury's chaotic secular evolution as a subdiffusive process

Authors: Dorian S. Abbot, Robert J. Webber, David M. Hernandez, Sam Hadden, Jonathan Weare

Abstract: Mercury's orbit can destabilize, generally resulting in a collision with either Venus or the Sun. Chaotic evolution can cause g1 to decrease to the approximately constant value of g5 and create a resonance. Previous work has approximated the variation in g1 as stochastic diffusion, which leads to a phenomological model that can reproduce the Mercury instability statistics of secular and N-body mod… ▽ More Mercury's orbit can destabilize, generally resulting in a collision with either Venus or the Sun. Chaotic evolution can cause g1 to decrease to the approximately constant value of g5 and create a resonance. Previous work has approximated the variation in g1 as stochastic diffusion, which leads to a phenomological model that can reproduce the Mercury instability statistics of secular and N-body models on timescales longer than 10 Gyr. Here we show that the diffusive model underpredicts the Mercury instability probability by a factor of 3-10,000 on timescales less than 5 Gyr, the remaining lifespan of the Solar System. This is because g1 exhibits larger variations on short timescales than the diffusive model would suggest. To better model the variations on short timescales, we build a new subdiffusive phenomological model for g1. Subdiffusion is similar to diffusion but exhibits larger displacements on short timescales and smaller displacements on long timescales. We choose model parameters based on the behavior of the g1 trajectories in the N-body simulations, leading to a tuned model that can reproduce Mercury instability statistics from 1-40 Gyr. This work motivates fundamental questions in Solar System dynamics: Why does subdiffusion better approximate the variation in g1 than standard diffusion? Why is there an upper bound on g1, but not a lower bound that would prevent it from reaching g5? △ Less

Submitted 12 April, 2024; v1 submitted 20 June, 2023; originally announced June 2023.

Comments: accepted at ApJ

Intra-system uniformity: a natural outcome of dynamical sculpting

Authors: Caleb Lammers, Sam Hadden, Norman Murray

Abstract: There is evidence that exoplanet systems display intra-system uniformity in mass, radius, and orbital spacing (like ``peas in a pod'') when compared with the system-to-system variations of planetary systems. This has been interpreted as the outcome of the early stages of planet formation, indicative of a picture in which planets form at characteristic mass scales with uniform separations. In this… ▽ More There is evidence that exoplanet systems display intra-system uniformity in mass, radius, and orbital spacing (like ``peas in a pod'') when compared with the system-to-system variations of planetary systems. This has been interpreted as the outcome of the early stages of planet formation, indicative of a picture in which planets form at characteristic mass scales with uniform separations. In this paper, we argue instead that intra-system uniformity in planet sizes and orbital spacings likely arose from the dynamical sculpting of initially-overly-packed planetary systems (in other words, the giant impact phase). With a suite of $N$-body simulations, we demonstrate that systems with random initial masses and compact planet spacings naturally develop intra-system uniformity, in quantitative agreement with observations, due to collisions between planets. Our results suggest that the pre-giant impact planet mass distribution is fairly wide and provide evidence for the prevalence of dynamical sculpting in shaping the observed population of exoplanets. △ Less

Submitted 17 July, 2023; v1 submitted 5 April, 2023; originally announced April 2023.

Comments: 6 pages, 5 figures. Accepted for publication in MNRAS Letters

Self-Consistent Spin, Tidal and Dynamical Equations of Motion in the REBOUNDx Framework

Authors: Tiger Lu, Hanno Rein, Daniel Tamayo, Sam Hadden, Rosemary Mardling, Sarah C. Millholland, Gregory Laughlin

Abstract: We have introduced self-consistent spin, tidal and dynamical equations of motion into REBOUNDx, a library of additional effects for the popular N-body integrator REBOUND. The equations of motion used are derived from the constant time lag approximation to the equilibrium tide model of tidal friction. These effects will allow the study of a variety of systems where the full dynamical picture cannot… ▽ More We have introduced self-consistent spin, tidal and dynamical equations of motion into REBOUNDx, a library of additional effects for the popular N-body integrator REBOUND. The equations of motion used are derived from the constant time lag approximation to the equilibrium tide model of tidal friction. These effects will allow the study of a variety of systems where the full dynamical picture cannot be encapsulated by point particle dynamics. We provide several test cases and benchmark the code's performance against analytic predictions. The open-source code is available in the most recent release of REBOUNDx. △ Less

Submitted 28 February, 2023; originally announced March 2023.

Comments: Accepted to ApJ. 15 pages, 6 figures

Simple physics and integrators accurately reproduce Mercury instability statistics

Authors: Dorian S. Abbot, David M. Hernandez, Sam Hadden, Robert J. Webber, Georgios P. Afentakis, Jonathan Weare

Abstract: The long-term stability of the Solar System is an issue of significant scientific and philosophical interest. The mechanism leading to instability is Mercury's eccentricity being pumped up so high that Mercury either collides with Venus or is scattered into the Sun. Previously, only three five-billion-year $N$-body ensembles of the Solar System with thousands of simulations have been run to assess… ▽ More The long-term stability of the Solar System is an issue of significant scientific and philosophical interest. The mechanism leading to instability is Mercury's eccentricity being pumped up so high that Mercury either collides with Venus or is scattered into the Sun. Previously, only three five-billion-year $N$-body ensembles of the Solar System with thousands of simulations have been run to assess long-term stability. We generate two additional ensembles, each with 2750 members, and make them publicly available at \texttt{https://archive.org/details/@dorianabbot}. We find that accurate Mercury instability statistics can be obtained by (1) including only the Sun and the 8 planets, (2) using a simple Wisdom-Holman scheme without correctors, (3) using a basic representation of general relativity, and (4) using a time step of 3.16 days. By combining our Solar System ensembles with previous ensembles we form a 9,601-member ensemble of ensembles. In this ensemble of ensembles, the logarithm of the frequency of a Mercury instability event increases linearly with time between 1.3 and 5 Gyr, suggesting that a single mechanism is responsible for Mercury instabilities in this time range and that this mechanism becomes more active as time progresses. Our work provides a robust estimate of Mercury instability statistics over the next five billion years, outlines methodologies that may be useful for exoplanet system investigations, and provides two large ensembles of publicly available Solar System integrations that can serve as testbeds for theoretical ideas as well as training sets for artificial intelligence schemes. △ Less

Submitted 21 February, 2023; v1 submitted 30 December, 2022; originally announced December 2022.

Comments: accepted at ApJ

Journal ref: Ap.J. 944 (2023) 190

celmech: A Python package for celestial mechanics

Authors: Sam Hadden, Daniel Tamayo

Abstract: We present celmech, an open-source Python package designed to facilitate a wide variety of celestial mechanics calculations. The package allows users to formulate and integrate equations of motion incorporating user-specified terms from the classical disturbing function expansion of the interaction potential between pairs of planets. The code can be applied, for example, to isolate the contributio… ▽ More We present celmech, an open-source Python package designed to facilitate a wide variety of celestial mechanics calculations. The package allows users to formulate and integrate equations of motion incorporating user-specified terms from the classical disturbing function expansion of the interaction potential between pairs of planets. The code can be applied, for example, to isolate the contribution of particular resonances to a system's dynamical evolution and develop simple analytical models with the minimum number of terms required to capture a particular dynamical phenomenon. Equations and expressions can be easily manipulated by leveraging the extensive symbolic mathematics capabilities of the sympy Python package. The celmech package is designed to interface seamlessly with the popular $N$-body code REBOUND to facilitate comparisons between calculation results and direct $N$-body integrations. The code is extensively documented and numerous example Jupyter notebooks illustrating its use are available online. △ Less

Submitted 20 May, 2022; originally announced May 2022.

Comments: submitted to AAS Journals, code and example notebooks available at http://github.com/shadden/celmech ; documentation available at http://celmech.readthedocs.io

The Effects of Disk Induced Apsidal Precession on Planets Captured into Mean Motion Resonance

Authors: Zachary Murray, Sam Hadden, Matthew J. Holman

Abstract: The process of migration into resonance capture has been well studied for planetary systems where the gravitational potential is generated exclusively by the star and planets. However, massive protoplanetary disks add a significant perturbation to these models. In this paper we consider two limiting cases of disk-induced precession on migrating planets and find that small amounts of precession sig… ▽ More The process of migration into resonance capture has been well studied for planetary systems where the gravitational potential is generated exclusively by the star and planets. However, massive protoplanetary disks add a significant perturbation to these models. In this paper we consider two limiting cases of disk-induced precession on migrating planets and find that small amounts of precession significantly affect the equilibrium reached by migrating planets. We investigate these effects with a combination of semi-analytic models of the resonance and numerical integrations. We also consider the case of the disk's dispersal, which can excite significant libration amplitude and can cause ejection from resonance for large enough precession rates. Both of these effects have implications for interpreting the known exoplanet population and may prove to be important considerations as the population of well-characterized exoplanet systems continues to grow. △ Less

Submitted 25 April, 2022; originally announced April 2022.

Stepsize errors in the $N$-body problem: discerning Mercury's true possible long-term orbits

Authors: David M. Hernandez, Richard E. Zeebe, Sam Hadden

Abstract: Numerical integrations of the Solar System have been carried out for decades. Their results have been used, for example, to determine whether the Solar System is chaotic, whether Mercury's orbit is stable, or to help discern Earth's climate history. We argue that all of the past studies we consider in this work are affected by numerical chaos to different degrees, affecting the possible orbits and… ▽ More Numerical integrations of the Solar System have been carried out for decades. Their results have been used, for example, to determine whether the Solar System is chaotic, whether Mercury's orbit is stable, or to help discern Earth's climate history. We argue that all of the past studies we consider in this work are affected by numerical chaos to different degrees, affecting the possible orbits and instability probability of Mercury, sometimes significantly. We show how to eliminate the effects of numerical chaos by resolving Mercury's pericentre passage. We also show that several higher order symplectic maps do not exhibit significant differences in resolving pericentre passage of Mercury (at fixed time step), making their advantages suspect for calculating long-term orbits. Resolving pericentre passage affects a wide array of orbital numerical studies, like exoplanet studies, studies of the galactic centre, and other $N$-body problems. △ Less

Submitted 13 December, 2021; v1 submitted 16 November, 2021; originally announced November 2021.

Comments: 6 pages, 6 figures. Matches accepted MNRAS version

The Criterion for Chaos in Three-Planet Systems

Authors: Jeremy Rath, Sam Hadden, Yoram Lithwick

Abstract: We establish the criterion for chaos in three-planet systems, for systems similar to those discovered by the Kepler spacecraft. Our main results are as follows: (i) The simplest criterion, which is based on overlapping mean motion resonances MMR's), only agrees with numerical simulations at a very crude level. (ii) Much greater accuracy is attained by considering neighboring MMR's that do not over… ▽ More We establish the criterion for chaos in three-planet systems, for systems similar to those discovered by the Kepler spacecraft. Our main results are as follows: (i) The simplest criterion, which is based on overlapping mean motion resonances MMR's), only agrees with numerical simulations at a very crude level. (ii) Much greater accuracy is attained by considering neighboring MMR's that do not overlap. We work out the width of the chaotic zones around each of the neighbors, and also provide simple approximate expressions for the widths. (iii) Even greater accuracy is provided by the overlap of three-body resonances (3BR's), which accounts for fine-grained structure seen in maps from N-body simulations, and also predicts the Lyapunov times. Previous studies conflict on whether overlap of MMR's or of 3BR's drive interplanetary chaos. We show that both do, and in fact they are merely different ways of looking at the same effect. (iv) We compare both criteria with high-resolution maps of chaos from N-body simulations, and show that they agree at a high level of detail. △ Less

Submitted 6 October, 2021; originally announced October 2021.

Rare Event Sampling Improves Mercury Instability Statistics

Authors: Dorian S. Abbot, Robert J. Webber, Sam Hadden, Darryl Seligman, Jonathan Weare

Abstract: Due to the chaotic nature of planetary dynamics, there is a non-zero probability that Mercury's orbit will become unstable in the future. Previous efforts have estimated the probability of this happening between 3 and 5 billion years in the future using a large number of direct numerical simulations with an N-body code, but were not able to obtain accurate estimates before 3 billion years in the f… ▽ More Due to the chaotic nature of planetary dynamics, there is a non-zero probability that Mercury's orbit will become unstable in the future. Previous efforts have estimated the probability of this happening between 3 and 5 billion years in the future using a large number of direct numerical simulations with an N-body code, but were not able to obtain accurate estimates before 3 billion years in the future because Mercury instability events are too rare. In this paper we use a new rare event sampling technique, Quantile Diffusion Monte Carlo (QDMC), to estimate that the probability of a Mercury instability event in the next 2 billion years is approximately $10^{-4}$ in the REBOUND N-body code. We show that QDMC provides unbiased probability estimates at a computational cost of up to 100 times less than direct numerical simulation. QDMC is easy to implement and could be applied to many problems in planetary dynamics in which it is necessary to estimate the probability of a rare event. △ Less

Submitted 27 December, 2021; v1 submitted 16 June, 2021; originally announced June 2021.

Comments: Abbot et al 2021 ApJ 923 236

How Close are Compact Multi-Planet Systems to the Stability Limit?

Authors: Samuel W. Yee, Daniel Tamayo, Samuel Hadden, Joshua N. Winn

Abstract: Transit surveys have revealed a significant population of compact multi-planet systems, containing several sub-Neptune-mass planets on close-in, tightly-packed orbits. These systems are thought to have formed through a final phase of giant impacts, which would tend to leave systems close to the edge of stability. Here, we assess this hypothesis, comparing observed eccentricities in systems exhibit… ▽ More Transit surveys have revealed a significant population of compact multi-planet systems, containing several sub-Neptune-mass planets on close-in, tightly-packed orbits. These systems are thought to have formed through a final phase of giant impacts, which would tend to leave systems close to the edge of stability. Here, we assess this hypothesis, comparing observed eccentricities in systems exhibiting transit-timing variations (TTVs), with the maximum eccentricities compatible with long-term stability. We use the machine-learning classifier SPOCK (Tamayo et al. 2020) to rapidly classify the stability of numerous initial configurations and hence determine these stability limits. While previous studies have argued that multi-planet systems are often maximally packed, in the sense that they could not host any additional planets, we find that the existing planets in these systems have measured eccentricities below the limits allowed by stability by a factor of 2--10. We compare these results against predictions from the giant impact theory of planet formation, derived from both $N$-body integrations and theoretical expectations that in the absence of dissipation, the orbits of such planets should be distributed uniformly throughout the phase space volume allowed by stability. We find that the observed systems have systematically lower eccentricities than this scenario predicts, with a median eccentricity about 4 times lower than predicted. These findings suggest that if such systems formed through giant impacts, then some dissipation must occur to damp their eccentricities. This may take place during formation, perhaps through interactions with the natal gas disk or a leftover population of planetesimals, or over longer timescales through the coupling of tidal and secular processes. △ Less

Submitted 13 May, 2021; originally announced May 2021.

Comments: 14 pages, 7 figures, 1 machine-readable table. Accepted to AJ

Significant improvement in planetary system simulations from statistical averaging

Authors: David M. Hernandez, Eric Agol, Matthew J. Holman, Sam Hadden

Abstract: Symplectic integrators are widely used for the study of planetary dynamics and other $N$-body problems. In a study of the outer Solar system, we demonstrate that individual symplectic integrations can yield biased errors in the semi-major axes and possibly other orbital elements. The bias is resolved by studying an ensemble of initial conditions of the outer Solar system. Such statistical sampling… ▽ More Symplectic integrators are widely used for the study of planetary dynamics and other $N$-body problems. In a study of the outer Solar system, we demonstrate that individual symplectic integrations can yield biased errors in the semi-major axes and possibly other orbital elements. The bias is resolved by studying an ensemble of initial conditions of the outer Solar system. Such statistical sampling could significantly improve measurement of planetary system properties like their secular frequencies. We also compared the distributions of action-like variables between high and low accuracy integrations; traditional statistical metrics are unable to distinguish the distribution functions. △ Less

Submitted 5 April, 2021; originally announced April 2021.

Comments: 4 pages, 1 figure. Accepted by RNAAS

Journal ref: 2021 Res. Notes AAS 5 77

A Bayesian neural network predicts the dissolution of compact planetary systems

Authors: Miles Cranmer, Daniel Tamayo, Hanno Rein, Peter Battaglia, Samuel Hadden, Philip J. Armitage, Shirley Ho, David N. Spergel

Abstract: Despite over three hundred years of effort, no solutions exist for predicting when a general planetary configuration will become unstable. We introduce a deep learning architecture to push forward this problem for compact systems. While current machine learning algorithms in this area rely on scientist-derived instability metrics, our new technique learns its own metrics from scratch, enabled by a… ▽ More Despite over three hundred years of effort, no solutions exist for predicting when a general planetary configuration will become unstable. We introduce a deep learning architecture to push forward this problem for compact systems. While current machine learning algorithms in this area rely on scientist-derived instability metrics, our new technique learns its own metrics from scratch, enabled by a novel internal structure inspired from dynamics theory. Our Bayesian neural network model can accurately predict not only if, but also when a compact planetary system with three or more planets will go unstable. Our model, trained directly from short N-body time series of raw orbital elements, is more than two orders of magnitude more accurate at predicting instability times than analytical estimators, while also reducing the bias of existing machine learning algorithms by nearly a factor of three. Despite being trained on compact resonant and near-resonant three-planet configurations, the model demonstrates robust generalization to both non-resonant and higher multiplicity configurations, in the latter case outperforming models fit to that specific set of integrations. The model computes instability estimates up to five orders of magnitude faster than a numerical integrator, and unlike previous efforts provides confidence intervals on its predictions. Our inference model is publicly available in the SPOCK package, with training code open-sourced. △ Less

Submitted 11 January, 2021; originally announced January 2021.

Comments: 8 content pages, 7 appendix and references. 8 figures. Source code at: https://github.com/MilesCranmer/bnn_chaos_model inference code at https://github.com/dtamayo/spock

Mildly-Hierarchical triple dynamics and applications to the outer solar system

Authors: Hareesh Gautham Bhaskar, Gongjie Li, Sam Hadden, Matthew J. Payne, Matthew J. Holman

Abstract: Three-body interactions are ubiquitous in astrophysics. For instance, Kozai-Lidov oscillations in hierarchical triple systems have been studied extensively and applied to a wide range of astrophysical systems. However, mildly-hierarchical triples also play an important role, but they are less explored. In this work we consider the secular dynamics of a test particle in a mildly-hierarchical config… ▽ More Three-body interactions are ubiquitous in astrophysics. For instance, Kozai-Lidov oscillations in hierarchical triple systems have been studied extensively and applied to a wide range of astrophysical systems. However, mildly-hierarchical triples also play an important role, but they are less explored. In this work we consider the secular dynamics of a test particle in a mildly-hierarchical configuration. We find the limit within which the secular approximation is reliable, present resonances and chaotic regions using surface of sections, and characterize regions of phase space that allow large eccentricity and inclination variations. Finally, we apply the secular results to the outer solar system. We focus on the distribution of extreme trans-neptunian objects (eTNOs) under the perturbation of a possible outer planet (Planet-9), and find that in addition to a low inclination Planet-9, a polar or a counter-orbiting one could also produce pericenter clustering of eTNOs, while the polar one leads to a wider spread of eTNO inclinations. △ Less

Submitted 30 November, 2020; v1 submitted 10 August, 2020; originally announced August 2020.

Comments: Accepted for publication in AJ

Predicting the long-term stability of compact multiplanet systems

Authors: Daniel Tamayo, Miles Cranmer, Samuel Hadden, Hanno Rein, Peter Battaglia, Alysa Obertas, Philip J. Armitage, Shirley Ho, David Spergel, Christian Gilbertson, Naireen Hussain, Ari Silburt, Daniel Jontof-Hutter, Kristen Menou

Abstract: We combine analytical understanding of resonant dynamics in two-planet systems with machine learning techniques to train a model capable of robustly classifying stability in compact multi-planet systems over long timescales of $10^9$ orbits. Our Stability of Planetary Orbital Configurations Klassifier (SPOCK) predicts stability using physically motivated summary statistics measured in integrations… ▽ More We combine analytical understanding of resonant dynamics in two-planet systems with machine learning techniques to train a model capable of robustly classifying stability in compact multi-planet systems over long timescales of $10^9$ orbits. Our Stability of Planetary Orbital Configurations Klassifier (SPOCK) predicts stability using physically motivated summary statistics measured in integrations of the first $10^4$ orbits, thus achieving speed-ups of up to $10^5$ over full simulations. This computationally opens up the stability constrained characterization of multi-planet systems. Our model, trained on $\approx 100,000$ three-planet systems sampled at discrete resonances, generalizes both to a sample spanning a continuous period-ratio range, as well as to a large five-planet sample with qualitatively different configurations to our training dataset. Our approach significantly outperforms previous methods based on systems' angular momentum deficit, chaos indicators, and parametrized fits to numerical integrations. We use SPOCK to constrain the free eccentricities between the inner and outer pairs of planets in the Kepler-431 system of three approximately Earth-sized planets to both be below 0.05. Our stability analysis provides significantly stronger eccentricity constraints than currently achievable through either radial velocity or transit duration measurements for small planets, and within a factor of a few of systems that exhibit transit timing variations (TTVs). Given that current exoplanet detection strategies now rarely allow for strong TTV constraints (Hadden et al., 2019), SPOCK enables a powerful complementary method for precisely characterizing compact multi-planet systems. We publicly release SPOCK for community use. △ Less

Submitted 14 July, 2020; v1 submitted 13 July, 2020; originally announced July 2020.

Comments: Published week of July 13th in Proceedings of the National Academy of Sciences: https://www.pnas.org/cgi/doi/10.1073/pnas.2001258117. Check out simple usage of package (and regenerate paper figures) at: https://github.com/dtamayo/spock

Modeling Radial Velocity Data of Resonant Planets to Infer Migration Histories

Authors: Sam Hadden, Matthew J. Payne

Abstract: A number of giant planet pairs discovered by the radial velocity method with period ratios $\lesssim 2$ may reside in mean motion resonances. Convergent orbital migration and resonant capture at the time of formation would naturally explain the present-day resonant orbital configurations of these systems. Planets that experience smooth migration and eccentricity damping forces due to a proto-plane… ▽ More A number of giant planet pairs discovered by the radial velocity method with period ratios $\lesssim 2$ may reside in mean motion resonances. Convergent orbital migration and resonant capture at the time of formation would naturally explain the present-day resonant orbital configurations of these systems. Planets that experience smooth migration and eccentricity damping forces due to a proto-planetary disk should not only capture into mean motion resonances but also end up in a specific dynamical configuration within the resonance, sometimes referred to as apsidal corotation resonance (ACR). Here we develop a method for testing the hypothesis that a planet pair resides in an ACR by directly fitting radial velocity data. The ACR hypothesis strongly restricts the number of free parameters describing the radial velocity signal and we compare fits using this highly restricted model to fits using a more conventional two-planet RV model by using nested sampling simulations. We apply our method to HD 45364 and HD 33844, two systems hosting giant planet pairs in 3:2 and 5:3 resonances, respectively. We demonstrate that the observations of both systems support an ACR configuration and we use the results of our ACR model fits to constrain possible migration histories of these systems. △ Less

Submitted 28 July, 2020; v1 submitted 6 May, 2020; originally announced May 2020.

Comments: Revised in response to referee report and accepted to AJ; Code for RV fitting available online at https://github.com/shadden/ResonantPlanetPairsRVModeling

Are long-term $N$-body simulations reliable?

Authors: David M. Hernandez, Sam Hadden, Junichiro Makino

Abstract: $N$-body integrations are used to model a wide range of astrophysical dynamics, but they suffer from errors which make their orbits diverge exponentially in time from the correct orbits. Over long time-scales, their reliability needs to be established. We address this reliability by running a three-body planetary system over about $200$ e-folding times. Using nearby initial conditions, we can cons… ▽ More $N$-body integrations are used to model a wide range of astrophysical dynamics, but they suffer from errors which make their orbits diverge exponentially in time from the correct orbits. Over long time-scales, their reliability needs to be established. We address this reliability by running a three-body planetary system over about $200$ e-folding times. Using nearby initial conditions, we can construct statistics of the long-term phase-space structure and compare to rough estimates of resonant widths of the system. We compared statistics for a wide range of numerical methods, including a Runge--Kutta method, Wisdom--Holman method, symplectic corrector methods, and a method by Laskar and Robutel. "Improving" an integrator did not increase the phase space accuracy, but simply increasing the number of initial conditions did. In fact, the statistics of a higher order symplectic corrector method were inconsistent with the other methods in one test. △ Less

Submitted 6 February, 2020; v1 submitted 18 October, 2019; originally announced October 2019.

Comments: 14 pages, 13 figures, 4 tables. Matches accepted MNRAS version

An Integrable Model for the Dynamics of Planetary Mean Motion Resonances

Authors: Sam Hadden

Abstract: I consider the dynamics of mean motion resonances between pairs of co-planar planets and derive a new integrable Hamiltonian model for planets' resonant motion. The new model generalizes previously-derived integrable Hamiltonians for first-order resonances to treat higher-order resonances by exploiting a surprising near-symmetry of the full, non-integrable Hamiltonians of higher-order resonances.… ▽ More I consider the dynamics of mean motion resonances between pairs of co-planar planets and derive a new integrable Hamiltonian model for planets' resonant motion. The new model generalizes previously-derived integrable Hamiltonians for first-order resonances to treat higher-order resonances by exploiting a surprising near-symmetry of the full, non-integrable Hamiltonians of higher-order resonances. Whereas past works have frequently relied on truncated disturbing function expansions to derive integrable approximations to resonant motion, I show that no such expansion is necessary, thus enabling the new model to accurately capture the dynamics of both first- and higher-order resonances for eccentricities up to orbit-crossing. I demonstrate that predictions of the new integrable model agree well with numerical integrations of resonant planet pairs. Finally, I explore the secular evolution of resonant planets' eccentricities. I show that the secular dynamics are governed by conservation of an AMD-like quantity. I also demonstrate that secular frequencies depend on planets' resonant libration amplitude and this generally gives rise to a secular resonance inside the mean motion resonance at large libration amplitudes. Outside of the secular resonance the long-term dynamics are characterized small adiabatic modulations of the resonant motion while inside the secular resonance planets can experience large variations of the resonant trajectory over secular timescales. The integrable model derived in this work can serve as a framework for analyzing the dynamics of planetary MMRs in a wide variety of contexts. △ Less

Submitted 11 September, 2019; originally announced September 2019.

Comments: submitted to AAS Journals; code available at github.com/shadden/IntegrableModelforResonances; comments welcome

The L 98-59 System: Three Transiting, Terrestrial-Sized Planets Orbiting a Nearby M-dwarf

Authors: Veselin B. Kostov, Joshua E. Schlieder, Thomas Barclay, Elisa V. Quintana, Knicole D. Colon, Jonathan Brande, Karen A. Collins, Adina D. Feinstein, Samuel Hadden, Stephen R. Kane, Laura Kreidberg, Ethan Kruse, Christopher Lam, Elisabeth Matthews, Benjamin T. Montet, Francisco J. Pozuelos, Keivan G. Stassun, Jennifer G. Winters, George Ricker, Roland Vanderspek, David Latham, Sara Seager, Joshua Winn, Jon M. Jenkins, Dennis Afanasev , et al. (90 additional authors not shown)

Abstract: We report the Transiting Exoplanet Survey Satellite (TESS) discovery of three terrestrial-sized planets transiting L 98-59 (TOI-175, TIC 307210830) -- a bright M dwarf at a distance of 10.6 pc. Using the Gaia-measured distance and broad-band photometry we find that the host star is an M3 dwarf. Combined with the TESS transits from three sectors, the corresponding stellar parameters yield planet ra… ▽ More We report the Transiting Exoplanet Survey Satellite (TESS) discovery of three terrestrial-sized planets transiting L 98-59 (TOI-175, TIC 307210830) -- a bright M dwarf at a distance of 10.6 pc. Using the Gaia-measured distance and broad-band photometry we find that the host star is an M3 dwarf. Combined with the TESS transits from three sectors, the corresponding stellar parameters yield planet radii ranging from 0.8REarth to 1.6REarth. All three planets have short orbital periods, ranging from 2.25 to 7.45 days with the outer pair just wide of a 2:1 period resonance. Diagnostic tests produced by the TESS Data Validation Report and the vetting package DAVE rule out common false positive sources. These analyses, along with dedicated follow-up and the multiplicity of the system, lend confidence that the observed signals are caused by planets transiting L 98-59 and are not associated with other sources in the field. The L 98-59 system is interesting for a number of reasons: the host star is bright (V = 11.7 mag, K = 7.1 mag) and the planets are prime targets for further follow-up observations including precision radial-velocity mass measurements and future transit spectroscopy with the James Webb Space Telescope; the near resonant configuration makes the system a laboratory to study planetary system dynamical evolution; and three planets of relatively similar size in the same system present an opportunity to study terrestrial planets where other variables (age, metallicity, etc.) can be held constant. L 98-59 will be observed in 4 more TESS sectors, which will provide a wealth of information on the three currently known planets and have the potential to reveal additional planets in the system. △ Less

Submitted 28 May, 2019; v1 submitted 19 March, 2019; originally announced March 2019.

Comments: 27 pages, 22 figures, AJ accepted

Near-resonance in a system of sub-Neptunes from TESS

Authors: Samuel N. Quinn, Juliette C. Becker, Joseph E. Rodriguez, Sam Hadden, Chelsea X. Huang, Timothy D. Morton, Fred Adams, David Armstrong, Jason D. Eastman, Jonathan Horner, Stephen R. Kane, Jack J. Lissauer, Joseph D. Twicken, Andrew Vanderburg, Rob Wittenmyer, George R. Ricker, Roland K. Vanderspek, David W. Latham, Sara Seager, Joshua N. Winn, Jon M. Jenkins, Eric Agol, Khalid Barkaoui, Charles A. Beichman, François Bouchy , et al. (56 additional authors not shown)

Abstract: We report the Transiting Exoplanet Survey Satellite ($TESS$) detection of a multi-planet system orbiting the $V=10.9$ K0 dwarf TOI 125. We find evidence for up to five planets, with varying confidence. Three high signal-to-noise transit signals correspond to sub-Neptune-sized planets ($2.76$, $2.79$, and $2.94\ R_{\oplus}$), and we statistically validate the planetary nature of the two inner plane… ▽ More We report the Transiting Exoplanet Survey Satellite ($TESS$) detection of a multi-planet system orbiting the $V=10.9$ K0 dwarf TOI 125. We find evidence for up to five planets, with varying confidence. Three high signal-to-noise transit signals correspond to sub-Neptune-sized planets ($2.76$, $2.79$, and $2.94\ R_{\oplus}$), and we statistically validate the planetary nature of the two inner planets ($P_b = 4.65$ days, $P_c = 9.15$ days). With only two transits observed, we report the outer object ($P_{.03} = 19.98$ days) as a high signal-to-noise ratio planet candidate. We also detect a candidate transiting super-Earth ($1.4\ R_{\oplus}$) with an orbital period of only $12.7$ hours and a candidate Neptune-sized planet ($4.2\ R_{\oplus}$) with a period of $13.28$ days, both at low signal-to-noise. This system is amenable to mass determination via radial velocities and transit timing variations, and provides an opportunity to study planets of similar size while controlling for age and environment. The ratio of orbital periods between TOI 125 b and c ($P_c/P_b = 1.97$) is slightly smaller than an exact 2:1 commensurability and is atypical of multiple planet systems from $Kepler$, which show a preference for period ratios just $wide$ of first-order period ratios. A dynamical analysis refines the allowed parameter space through stability arguments and suggests that, despite the nearly commensurate periods, the system is unlikely to be in resonance. △ Less

Submitted 25 January, 2019; originally announced January 2019.

Comments: Submitted to AAS Journals. 20 pages, 10 figures, 5 tables

Prospects for TTV Detection and Dynamical Constraints with TESS

Authors: Sam Hadden, Thomas Barclay, Matthew J. Payne, Matthew J. Holman

Abstract: We consider the potential for the Transiting Exoplanet Survey Satellite (TESS) to detect transit timing variations (TTVs) during both its nominal and extended mission phases. Building on previous estimates of the overall yield of planetary systems from the TESS mission, we predict that during its nominal two-year mission, TESS will observe measurable TTVs in $\sim30$ systems, from which… ▽ More We consider the potential for the Transiting Exoplanet Survey Satellite (TESS) to detect transit timing variations (TTVs) during both its nominal and extended mission phases. Building on previous estimates of the overall yield of planetary systems from the TESS mission, we predict that during its nominal two-year mission, TESS will observe measurable TTVs in $\sim30$ systems, from which $\mathcal{O}(10)$ planet will get precise mass measurements from TTVs alone, $\sim 5$ planets will have significant constraints placed on their masses from TTVs, and over a dozen systems will be singly-transiting TTV systems. We consider a number of different extended mission scenarios, and predict that in a typical scenario, an extended mission will allow TESS to increase the number of systems with measurable TTVs to a total of $\sim90$, from which $\sim 15$ planets will have precise mass measurements, another $\sim 15$ will have significant constraints placed on their masses, and $\sim60$ will be singly-transiting TTV systems. We also describe how follow-up transit observations of multi-planet systems discovered by the TESS mission can be optimally planned to maximize TTV mass and eccentricity constraints. △ Less

Submitted 3 September, 2019; v1 submitted 5 November, 2018; originally announced November 2018.

Comments: Revised in response to referee report; includes updated TTV yield projections; accepted to AJ; code available at github.com/shadden/TTV2Fast2Furious

Prospects for Refining Kepler TTV Masses using TESS Observations

Authors: Max Goldberg, Sam Hadden, Matthew J. Payne, Matthew J. Holman

Abstract: In this paper we investigate systems previously identified to exhibit transit timing variations (TTVs) in Kepler data, with the goal of predicting the expected improvements to the mass and eccentricity constraints that will arise from combining Kepler data with future data from the TESS mission. We advocate for the use of the Kullback-Leibler (KL) divergence as a means to quantify improvements in… ▽ More In this paper we investigate systems previously identified to exhibit transit timing variations (TTVs) in Kepler data, with the goal of predicting the expected improvements to the mass and eccentricity constraints that will arise from combining Kepler data with future data from the TESS mission. We advocate for the use of the Kullback-Leibler (KL) divergence as a means to quantify improvements in the measured constraints. Compared to the original Kepler data, the TESS data will have a lower signal-to-noise ratio, rendering some of the planetary transits undetectable, and lowering the accuracy with which the transit mid-time can be estimated. Despite these difficulties, out of the 55 systems (containing 143 planets) investigated, we predict that the collection of short-cadence data by TESS will be of significant value (i.e. it will improve the mass uncertainty such that the KL divergence is > 0.1) for approximately 6 - 14 planets during the nominal mission, with the range primarily driven by the uncertain precision with which transit mid-times will be recovered from TESS data. In an extended mission this would increase to a total of approximately 12 - 25 planets. △ Less

Submitted 5 October, 2018; originally announced October 2018.

Comments: 11 pages, 8 figures, submitted to AAS journals

A Compact Multi-Planet System With A Significantly Misaligned Ultra Short Period Planet

Authors: Joseph E. Rodriguez, Juliette C. Becker, Jason D. Eastman, Sam Hadden, Andrew Vanderburg, Tali Khain, Samuel N. Quinn, Andrew Mayo, Courtney D. Dressing, Joshua E. Schlieder, David R. Ciardi, David W. Latham, Saul Rappaport, Fred C. Adams, Perry Berlind, Allyson Bieryla, Michael L. Calkins, Gilbert A. Esquerdo, Martti H. Kristiansen, Mark Omohundro, Hans Martin Schwengeler, Keivan G. Stassun, Ivan Terentev

Abstract: We report the discovery of a compact multi-planet system orbiting the relatively nearby (78pc) and bright ($K=8.9$) K-star, K2-266 (EPIC248435473). We identify up to six possible planets orbiting K2-266 with estimated periods of P$_b$ = 0.66, P$_{.02}$ = 6.1, P$_c$ = 7.8, P$_d$ = 14.7, P$_e$ = 19.5, and P$_{.06}$ = 56.7 days and radii of R$_P$ = 3.3 R$_{\oplus}$, 0.646 R$_{\oplus}$, 0.705 R… ▽ More We report the discovery of a compact multi-planet system orbiting the relatively nearby (78pc) and bright ($K=8.9$) K-star, K2-266 (EPIC248435473). We identify up to six possible planets orbiting K2-266 with estimated periods of P$_b$ = 0.66, P$_{.02}$ = 6.1, P$_c$ = 7.8, P$_d$ = 14.7, P$_e$ = 19.5, and P$_{.06}$ = 56.7 days and radii of R$_P$ = 3.3 R$_{\oplus}$, 0.646 R$_{\oplus}$, 0.705 R$_{\oplus}$, 2.93 R$_{\oplus}$, 2.73 R$_{\oplus}$, and 0.90 R$_{\oplus}$, respectively. We are able to confirm the planetary nature of two of these planets (d & e) from analyzing their transit timing variations ($m_d= 8.9_{-3.8}^{+5.7} M_\oplus$ and $m_e=14.3_{-5.0}^{+6.4} M_\oplus$), confidently validate the planetary nature of two other planets (b & c), and classify the last two as planetary candidates (K2-266.02 & .06). From a simultaneous fit of all 6 possible planets, we find that K2-266 b's orbit has an inclination of 75.32$^{\circ}$ while the other five planets have inclinations of 87-90$^{\circ}$. This observed mutual misalignment may indicate that K2-266 b formed differently from the other planets in the system. The brightness of the host star and the relatively large size of the sub-Neptune sized planets d and e make them well-suited for atmospheric characterization efforts with facilities like the Hubble Space Telescope and upcoming James Webb Space Telescope. We also identify an 8.5-day transiting planet candidate orbiting EPIC248435395, a co-moving companion to K2-266. △ Less

Submitted 29 October, 2018; v1 submitted 21 June, 2018; originally announced June 2018.

Comments: 18 pages, 12 figures, 7 tables, Accepted for Publication in the Astronomical Journal

The Secular Dynamics of TNOs and Planet Nine Interactions

Authors: Gongjie Li, Samuel Hadden, Matthew Payne, Matthew J. Holman

Abstract: The existence of Planet Nine has been suggested to explain the pericenter clustering of extreme trans-Neptunian objects (TNOs). However, the underlying dynamics involving Planet Nine, test particles and Neptune is rich, and it remains unclear which dynamical processes lead to the alignment and how they depends on the properties of Planet Nine. Here, we investigate the secular interactions between… ▽ More The existence of Planet Nine has been suggested to explain the pericenter clustering of extreme trans-Neptunian objects (TNOs). However, the underlying dynamics involving Planet Nine, test particles and Neptune is rich, and it remains unclear which dynamical processes lead to the alignment and how they depends on the properties of Planet Nine. Here, we investigate the secular interactions between an eccentric outer perturber and TNOs starting in a near-coplanar configuration. We find that a large number of TNOs could survive outside of mean motion resonances at 4Gyr, which differs from previous results obtained in the exact coplanar case with Neptune being treated as a quadrupole potential. In addition, secular dynamics leads to the orbital clustering seen in N-body simulations. We find that a near coplanar Planet Nine can flip TNO orbital planes, and when this happens, the geometrical longitudes of pericenter of the TNOs librate around $180^\circ$ during the flip. Orbital precession caused by the inner giant planets can suppress the flips while keeping the longitude of pericenter librating when $30\lesssim r_p \lesssim 80$ AU $\&$ $a\gtrsim 250$ AU. This results in the alignment of the pericenter of the low inclination TNOs ($i \lesssim 40^\circ$). We find the anti-aligned population and the flipped orbits could be produced by an eccentric ($e_9 \gtrsim 0.4$) outer planet of $\sim 10M_{\oplus}$ in a wide $a_9 \gtrsim 400 \sim 800$ AU orbit. Future surveys on the high inclination TNOs will help further constrain the properties of possible outer planets. △ Less

Submitted 11 October, 2018; v1 submitted 18 June, 2018; originally announced June 2018.

Comments: 21 pages, 26 figures, accepted by AJ

A Criterion for the Onset of Chaos in Systems of Two Eccentric Planets

Authors: Sam Hadden, Yoram Lithwick

Abstract: We derive a criterion for the onset of chaos in systems consisting of two massive, eccentric, coplanar planets. Given the planets' masses and separation, the criterion predicts the critical eccentricity above which chaos is triggered. Chaos occurs where mean motion resonances overlap, as in Wisdom (1980)'s pioneering work. But whereas Wisdom considered only nearly circular planets, and hence exami… ▽ More We derive a criterion for the onset of chaos in systems consisting of two massive, eccentric, coplanar planets. Given the planets' masses and separation, the criterion predicts the critical eccentricity above which chaos is triggered. Chaos occurs where mean motion resonances overlap, as in Wisdom (1980)'s pioneering work. But whereas Wisdom considered only nearly circular planets, and hence examined only first order resonances, we extend his results to arbitrarily eccentric planets (up to crossing orbits) by examining resonances of all orders. We thereby arrive at a simple expression for the critical eccentricity. We do this first for a test particle in the presence of a planet, and then generalize to the case of two massive planets, based on a new approximation to the Hamiltonian (Hadden, in prep). We then confirm our results with detailed numerical simulations. Finally, we explore the extent to which chaotic two-planet systems eventually result in planetary collisions. △ Less

Submitted 9 August, 2018; v1 submitted 22 March, 2018; originally announced March 2018.

Comments: 17 pages, 13 figures. Accepted for publication in AJ. Revised in response to community and referee comments

Chaotic Dynamics of Trans-Neptunian Objects Perturbed by Planet Nine

Authors: Sam Hadden, Gongjie Li, Matthew J. Payne, Matthew J. Holman

Abstract: Observations of clustering among the orbits of the most distant trans-Neptunian objects (TNOs) has inspired interest in the possibility of an undiscovered ninth planet lurking in the outskirts of the solar system. Numerical simulations by a number of authors have demonstrated that, with appropriate choices of planet mass and orbit, such a planet can maintain clustering in the orbital elements of t… ▽ More Observations of clustering among the orbits of the most distant trans-Neptunian objects (TNOs) has inspired interest in the possibility of an undiscovered ninth planet lurking in the outskirts of the solar system. Numerical simulations by a number of authors have demonstrated that, with appropriate choices of planet mass and orbit, such a planet can maintain clustering in the orbital elements of the population of distant TNOs, similar to the observed sample. However, many aspects of the rich underlying dynamical processes induced by such a distant eccentric perturber have not been fully explored. We report the results of our investigation of the dynamics of coplanar test-particles that interact with a massive body on an circular orbit (Neptune) and a massive body on a more distant, highly eccentric orbit (the putative Planet Nine). We find that a detailed examination of our idealized simulations affords tremendous insight into the rich test-particle dynamics that are possible. In particular, we find that chaos and resonance overlap plays an important role in particles' dynamical evolution. We develop a simple mapping model that allows us to understand, in detail, the web of overlapped mean-motion resonances explored by chaotically evolving particles. We also demonstrate that gravitational interactions with Neptune can have profound effects on the orbital evolution of particles. Our results serve as a starting point for a better understanding of the dynamical behavior observed in more complicated simulations that can be used to constrain the mass and orbit of Planet 9. △ Less

Submitted 1 May, 2018; v1 submitted 18 December, 2017; originally announced December 2017.

Comments: Revised to address referree comments; accepted to AJ

The Gold Standard: Accurate Stellar and Planetary Parameters for Eight Kepler M dwarf Systems Enabled by Parallaxes

Authors: Andrew W. Mann, Trent Dupuy, Philip S. Muirhead, Marshall C. Johnson, Michael C. Liu, Megan Ansdell, Paul A. Dalba, Jonathan J. Swift, Sam Hadden

Abstract: We report parallaxes and proper motions from the Hawaii Infrared Parallax Program for eight nearby M dwarf stars with transiting exoplanets discovered by Kepler. We combine our directly measured distances with mass-luminosity and radius-luminosity relationships to significantly improve constraints on the host stars' properties. Our astrometry enables the identification of wide stellar companions t… ▽ More We report parallaxes and proper motions from the Hawaii Infrared Parallax Program for eight nearby M dwarf stars with transiting exoplanets discovered by Kepler. We combine our directly measured distances with mass-luminosity and radius-luminosity relationships to significantly improve constraints on the host stars' properties. Our astrometry enables the identification of wide stellar companions to the planet hosts. Within our limited sample, all the multi-transiting planet hosts (three of three) appear to be single stars, while nearly all (four of five) of the systems with a single detected planet have wide stellar companions. By applying strict priors on average stellar density from our updated radius and mass in our transit fitting analysis, we measure the eccentricity probability distributions for each transiting planet. Planets in single-star systems tend to have smaller eccentricities than those in binaries, although this difference is not significant in our small sample. In the case of Kepler-42bcd, where the eccentricities are known to be $\simeq$0, we demonstrate that such systems can serve as powerful tests of M dwarf evolutionary models by working in $L_*-ρ_*$ space. The transit-fit density for Kepler-42bcd is inconsistent with model predictions at 2.1$σ$ (22%), but matches more empirical estimates at 0.2$σ$ (2%), consistent with earlier results showing model radii of M dwarfs are under-inflated. Gaia will provide high precision parallaxes for the entire Kepler M dwarf sample, and TESS will identify more planets transiting nearby, late-type stars, enabling significant improvements in our understanding of the eccentricity distribution of small planets and the parameters of late-type dwarfs. △ Less

Submitted 21 May, 2017; v1 submitted 3 May, 2017; originally announced May 2017.

Comments: Accepted to AJ. Updated with corrections from proofs (May 22, 2017)

Kepler Planet Masses and Eccentricities from TTV Analysis

Authors: Sam Hadden, Yoram Lithwick

Abstract: We conduct a uniform analysis of the transit timing variations (TTVs) of 145 planets from 55 Kepler multiplanet systems to infer planet masses and eccentricities. Eighty of these planets do not have previously reported mass and eccentricity measurements. We employ two complementary methods to fit TTVs: Markov chain Monte Carlo simulations based on N-body integration and an analytic fitting approac… ▽ More We conduct a uniform analysis of the transit timing variations (TTVs) of 145 planets from 55 Kepler multiplanet systems to infer planet masses and eccentricities. Eighty of these planets do not have previously reported mass and eccentricity measurements. We employ two complementary methods to fit TTVs: Markov chain Monte Carlo simulations based on N-body integration and an analytic fitting approach. Mass measurements of 49 planets, including 12 without previously reported masses, meet our criterion for classification as robust. Using mass and radius measurements, we infer the masses of planets' gaseous envelopes for both our TTV sample as well as transiting planets with radial velocity observations. Insight from analytic TTV formulae allows us to partially circumvent degeneracies inherent to inferring eccentricities from TTV observations. We find that planet eccentricities are generally small, typically a few percent, but in many instances are non-zero. △ Less

Submitted 10 November, 2016; originally announced November 2016.

Comments: 62 pages, 64 figures, submitted to AAS Journals

Numerical and Analytical Modelling of Transit Time Variations

Authors: Sam Hadden, Yoram Lithwick

Abstract: We develop and apply methods to extract planet masses and eccentricities from observed transit time variations (TTVs). First, we derive simple analytic expressions for the TTV that include the effects of both first- and second-order resonances. Second, we use N-body Markov chain Monte Carlo (MCMC) simulations, as well as the analytic formulae, to measure the masses and eccentricities of ten planet… ▽ More We develop and apply methods to extract planet masses and eccentricities from observed transit time variations (TTVs). First, we derive simple analytic expressions for the TTV that include the effects of both first- and second-order resonances. Second, we use N-body Markov chain Monte Carlo (MCMC) simulations, as well as the analytic formulae, to measure the masses and eccentricities of ten planets discovered by Kepler that have not previously been analyzed. Most of the ten planets have low densities. Using the analytic expressions to partially circumvent degeneracies, we measure small eccentricities of a few percent or less. △ Less

Submitted 8 October, 2015; originally announced October 2015.

Comments: 19 pages, 18 figures, submitted to ApJ

Densities and Eccentricities of 139 Kepler Planets from Transit Time Variations

Authors: Sam Hadden, Yoram Lithwick

Abstract: We extract densities and eccentricities of 139 sub-Jovian planets by analyzing transit time variations (TTVs) obtained by the Kepler mission through Quarter 12. We partially circumvent the degeneracies that plague TTV inversion with the help of an analytical formula for the TTV. From the observed TTV phases, we find that most of these planets have eccentricities of order a few percent. More precis… ▽ More We extract densities and eccentricities of 139 sub-Jovian planets by analyzing transit time variations (TTVs) obtained by the Kepler mission through Quarter 12. We partially circumvent the degeneracies that plague TTV inversion with the help of an analytical formula for the TTV. From the observed TTV phases, we find that most of these planets have eccentricities of order a few percent. More precisely, the r.m.s. eccentricity is 0.018^{+0.005}_{-0.004}, and planets smaller than 2.5R_\earth are around twice as eccentric as those bigger than 2.5R_\earth. We also find a best-fit density-radius relationship ρ~3 g/cm^3 \times (R/3R_\earth)^{-2.3} for the 56 planets that likely have small eccentricity and hence small statistical correction to their masses. Many planets larger than 2.5R_\earth are less dense than water, implying that their radii are largely set by a massive hydrogen atmosphere. △ Less

Submitted 8 April, 2014; v1 submitted 29 October, 2013; originally announced October 2013.

Comments: 8 pages, 7 figures, 2 tables; accepted to ApJ; Improved figures and data analysis; Tables 1 and 2 available in ancillary data and from https://sites.google.com/a/u.northwestern.edu/shadden/ttvs

Relativistic magnetohydrodynamics in one dimension

Authors: Maxim Lyutikov, Samuel Hadden

Abstract: We derive a number of solution for one-dimensional dynamics of relativistic magnetized plasma that can be used as benchmark estimates in relativistic hydrodynamic and magnetohydrodynamic numerical codes. First, we analyze the properties of simple waves of fast modes propagating orthogonally to the magnetic field in relativistically hot plasma. The magnetic and kinetic pressures obey different eq… ▽ More We derive a number of solution for one-dimensional dynamics of relativistic magnetized plasma that can be used as benchmark estimates in relativistic hydrodynamic and magnetohydrodynamic numerical codes. First, we analyze the properties of simple waves of fast modes propagating orthogonally to the magnetic field in relativistically hot plasma. The magnetic and kinetic pressures obey different equations of state, so that the system behaves as a mixture of gases with different polytropic indices. We find the self-similar solutions for the expansion of hot strongly magnetized plasma into vacuum. Second, we derive linear hodograph and Darboux equations for the relativistic Khalatnikov potential, which describe arbitrary one-dimensional isentropic relativistic motion of cold magnetized plasma and find their general and particular solutions. The obtained hodograph and Darboux equations are very powerful: system of highly non-linear, relativistic, time dependent equations describing arbitrary (not necessarily self-similar) dynamics of highly magnetized plasma reduces to a single linear differential equation. △ Less

Submitted 1 December, 2011; originally announced December 2011.

Comments: accepted by Phys. Rev. E