Researcher profile

Sam Hadden

Sam Hadden contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
7works
0followers
4topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

7 published item(s)

preprint2022arXiv

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

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.

preprint2021arXiv

Rare Event Sampling Improves Mercury Instability Statistics

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.

preprint2021arXiv

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

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.

preprint2021arXiv

The Criterion for Chaos in Three-Planet Systems

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.

preprint2020arXiv

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

$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.

preprint2020arXiv

Modeling Radial Velocity Data of Resonant Planets to Infer Migration Histories

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.

preprint2019arXiv

An Integrable Model for the Dynamics of Planetary Mean Motion Resonances

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.