BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Daniel Tamayo

Daniel Tamayo appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-ph.EPastro-ph.IMastro-ph.SRmath.CAmath.DSmath.NAnlin.CDphysics.comp-ph

Catalog footprint

What is connected

17works
8topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

17 published item(s)

preprint2026arXiv

Democratic heliocentric coordinates underestimate the rate of instabilities in long-term integrations of the Solar System

Wisdom-Holman (WH) integrators are symplectic operator-splitting methods widely used for long-term N-body simulations of planetary systems. Most implementations use either Jacobi coordinates or democratic heliocentric coordinates (DHC) for the Hamiltonian splitting, resulting in slightly different algorithms. In this paper we report results from numerical experiments, which show that integrations of the Solar System using DHC coordinates with typical timesteps of a few days suppress instabilities of the planet Mercury. We further show that this is due to an eccentricity dependent artificial numerical precession introduced by the DHC splitting. While the DHC splitting converges to the correct results at shorter timesteps of ~0.6 days, we argue that Jacobi coordinates remain reliable to significantly longer timesteps when orbits become moderately eccentric, and are thus a better choice when the innermost planet can reach high eccentricities.

0 citations0 reviewsNo code linkedOpen access
preprint2026arXiv

Self-consistent Dynamical and Chaotic Tides in the REBOUNDx framework

At high eccentricities, tidal forcing excites vibrational modes within orbiting bodies known as dynamical tides. In this paper, we implement the coupled evolution of these modes with the body's orbit in the \texttt{REBOUNDx} framework, an extension to the popular $N$-body integrator \texttt{REBOUND}. We provide a variety of test cases relevant to exoplanet dynamics and demonstrate overall agreement with prior studies of dynamical tides in the secular regime. Our implementation is readily applied to various high-eccentricity scenarios and allows for fast and accurate $N$-body investigations of astrophysical systems for which dynamical tides are relevant.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Stability Constrained Characterization of the 23 Myr-old V1298 Tau System: Do Young Planets Form in Mean Motion Resonance Chains?

A leading theoretical expectation for the final stages of planet formation is that disk migration should naturally drive orbits into chains of mean motion resonances (MMRs). In order to explain the dearth of MMR chains observed at Gyr ages ($<1\%$), this picture requires such configurations to destabilize and scramble period ratios following disk dispersal. Strikingly, the only two known stars with three or more planets younger than $\lesssim 100$ Myrs, HR 8799 and V1298 Tau, have been suggested to be in such MMR chains, given the orbits' near-integer period ratios. We incorporate recent transit and radial velocity observations of the V1298 Tau system, and investigate constraints on the system's orbital architecture imposed by requiring dynamical stability on timescales much shorter than the system's age. We show that the recent radial-velocity mass measurement of V1298 Tau $b$ places it within a factor of two of the instability limit, and that this allows us to set significantly lower limits on the eccentricity ($e_b \leq 0.17$ at $99.7\%$ confidence). Additionally, we rule out a resonant chain configuration for V1298 Tau at $\gtrsim 99\%$ confidence. Thus, if the $\sim 23$ Myr-old V1298 Tau system did form as a resonant chain, it must have undergone instability and rearrangement shortly after disk dispersal. We expect that similar stability constrained characterization of future young multi-planet systems will be valuable in informing planet formation models.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Stellar Evolution and Tidal Dissipation in REBOUNDx

To study the post-main sequence evolution of the Solar system and exoplanetary systems more accurately and efficiently, we introduce two new features to REBOUNDx, an extended library for the N-body integrator REBOUND. The first is a convenient parameter interpolator for coupling different physics and integrators using numerical splitting schemes. The second implements a constant time lag model for tides without evolving spins. We demonstrate various uses of these features using stellar evolution data from MESA (Modules for Experiments in Stellar Astrophysics) as an example. The results of our tests agree with several studies in the literature on post-main sequence orbital evolution, and our convergence and performance studies respectively demonstrate our implementations' accuracy and efficiency. These additional effects are publicly available as of REBOUNDx's latest release.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Direct Imaging of Irregular Satellite Disks in Scattered Light

Direct imaging surveys have found that long-period super-Jupiters are rare. By contrast, recent modeling of the widespread gaps in protoplanetary disks revealed by ALMA suggests an abundant population of smaller Neptune to Jupiter-mass planets at large separations. The thermal emission from such lower-mass planets is negligible at optical and near-infrared wavelengths, leaving only their weak signals in reflected light. Planets do not scatter enough light at these large orbital distances, but there is a natural way to enhance their reflecting area. Each of the four giant planets in our solar system hosts swarms of dozens of irregular satellites, gravitationally captured planetesimals that fill their host planets' spheres of gravitational influence. What we see of them today are the leftovers of an intense collisional evolution. At early times, they would have generated bright circumplanetary debris disks. We investigate the properties and detectability of such irregular satellite disks (ISDs) following models for their collisional evolution from Kennedy and Wyatt 2011. We find that ISD brightnesses peak in the $10-100$ AU range probed by ALMA, and can render planets detectable over a wide range of parameters with upcoming high-contrast instrumentation. We argue that future instruments with wide fields of view could simultaneously characterize the atmospheres of known close-in planets, and reveal the population of long-period $\sim$ Neptune-Jupiter mass exoplanets inaccessible to other detection methods. This provides a complementary and compelling science case that would elucidate the early lives of planetary systems.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Fundamental limits from chaos on instability time predictions in compact planetary systems

Instabilities in compact planetary systems are generically driven by chaotic dynamics. This implies that an instability time measured through direct N-body integration is not exact, but rather represents a single draw from a distribution of equally valid chaotic trajectories. In order to characterize the "errors" on reported instability times from direct N-body integrations, we investigate the shape and parameters of the instability time distributions (ITDs) for ensembles of shadow trajectories that are initially perturbed from one another near machine precision. We find that in the limit where instability times are long compared to the Lyapunov (chaotic) timescale, ITDs approach remarkably similar lognormal distributions with standard deviations ~0.43 $\pm$ 0.16 dex, despite the instability times varying across our sample from $10^4-10^8$ orbits. We find excellent agreement between these predictions, derived from ~450 closely packed configurations of three planets, and a much wider validation set of ~10,000 integrations, as well as on ~20,000 previously published integrations of tightly packed five-planet systems, and a seven-planet resonant chain based on TRAPPIST-1, despite their instability timescales extending beyond our analyzed timescale. We also test the boundary of applicability of our results on dynamically excited versions of our Solar System. These distributions define the fundamental limit imposed by chaos on the predictability of instability times in such planetary systems. It provides a quantitative estimate of the intrinsic error on an N-body instability time imprinted by chaos, approximately a factor of 3 in either direction.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Predicting the long-term stability of compact multiplanet systems

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Radial profiles of the Phoebe ring, a vast debris disk around Saturn

We present observations at optical wavelengths with the Cassini Spacecraft's Imaging Science System of the Phoebe ring, a vast debris disk around Saturn that seems to be collisionally generated by its irregular satellites. The analysis reveals a radial profile from 80-260 Saturn radii ($R_S$) that changes behavior interior to $\approx 110 R_S$, which we attribute to either the moon Iapetus sweeping up small particles, or to orbital instabilities that cause the ring to flare up vertically. Our study yields an integrated I/F at 0.635 $μ$m along Saturn's shadow in the Phoebe ring's midplane from 80-250 $R_S$ of $2.7^{+0.9}_{-0.3} \times 10^{-9}$. We develop an analytical model for the size-dependent secular dynamics of retrograde Phoebe ring grains, and compare this model to the observations. This analysis implies that 1) the "Phoebe" ring is partially sourced by debris from irregular satellites beyond Phoebe's orbit and 2) the scattered light signal is dominated by small grains ($\lesssim 20μ$m in size). If we assume that the Phoebe ring is generated through steady-state micrometeoroid bombardment, this implies a power-law size distribution with index $> 4$, which is unusually steep among solar system rings. This suggests either a steep size distribution of ejecta when material is initially released, or a subsequent process that preferentially breaks up large grains.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Second-order variational equations for N-body simulations

First-order variational equations are widely used in N-body simulations to study how nearby trajectories diverge from one another. These allow for efficient and reliable determinations of chaos indicators such as the Maximal Lyapunov characteristic Exponent (MLE) and the Mean Exponential Growth factor of Nearby Orbits (MEGNO). In this paper we lay out the theoretical framework to extend the idea of variational equations to higher order. We explicitly derive the differential equations that govern the evolution of second-order variations in the N-body problem. Going to second order opens the door to new applications, including optimization algorithms that require the first and second derivatives of the solution, like the classical Newton's method. Typically, these methods have faster convergence rates than derivative-free methods. Derivatives are also required for Riemann manifold Langevin and Hamiltonian Monte Carlo methods which provide significantly shorter correlation times than standard methods. Such improved optimization methods can be applied to anything from radial-velocity/transit-timing-variation fitting to spacecraft trajectory optimization to asteroid deflection. We provide an implementation of first and second-order variational equations for the publicly available REBOUND integrator package. Our implementation allows the simultaneous integration of any number of first and second-order variational equations with the high-accuracy IAS15 integrator. We also provide routines to generate consistent and accurate initial conditions without the need for finite differencing.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Tatooine's Future: The Eccentric Response of Kepler's Circumbinary Planets to Common-Envelope Evolution of their Host Stars

Inspired by the recent Kepler discoveries of circumbinary planets orbiting nine close binary stars, we explore the fate of the former as the latter evolve off the main sequence. We combine binary star evolution models with dynamical simulations to study the orbital evolution of these planets as their hosts undergo common-envelope stages, losing in the process a tremendous amount of mass on dynamical timescales. Five of the systems experience at least one Roche-lobe overflow and common-envelope stages (Kepler-1647 experiences three), and the binary stars either shrink to very short orbits or coalesce; two systems trigger a double-degenerate supernova explosion. Kepler's circumbinary planets predominantly remain gravitationally bound at the end of the common-envelope phase, migrate to larger orbits, and may gain significant eccentricity; their orbital expansion can be more than an order of magnitude and can occur over the course of a single planetary orbit. The orbits these planets can reach are qualitatively consistent with those of the currently known post-common-envelope, eclipse-time variations circumbinary candidates. Our results also show that circumbinary planets can experience both modes of orbital expansion (adiabatic and non-adiabatic) if their host binaries undergo more than one common-envelope stage; multiplanet circumbinary systems like Kepler-47 can experience both modes during the same common-envelope stage. Additionally, unlike Mercury orbiting the Sun, a circumbinary planet with the same semi-major axis can survive the common envelope evolution of a close binary star with a total mass of 1 MSun.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Could Jupiter or Saturn Have Ejected a Fifth Giant Planet?

Models of the dynamical evolution of the early solar system following the dispersal of the gaseous protoplanetary disk have been widely successful in reconstructing the current orbital configuration of the giant planets. Statistically, some of the most successful dynamical evolution simulations have initially included a hypothetical fifth giant planet, of ice giant mass, which gets ejected by a gas giant during the early solar system's proposed instability phase. We investigate the likelihood of an ice giant ejection event by either Jupiter or Saturn through constraints imposed by the current orbits of their wide-separation regular satellites Callisto and Iapetus respectively. We show that planetary encounters that are sufficient to eject an ice giant, often provide excessive perturbations to the orbits of Callisto and Iapetus making it difficult to reconcile a planet ejection event with the current orbit of either satellite. Quantitatively, we compute the likelihood of reconciling a regular Jovian satellite orbit with the current orbit of Callisto following an ice giant ejection by Jupiter of ~ 42% and conclude that such a large likelihood supports the hypothesis of a fifth giant planet's existence. A similar calculation for Iapetus reveals that it is much more difficult for Saturn to have ejected an ice giant and reconcile a Kronian satellite orbit with that of Iapetus (likelihood ~ 1%), although uncertainties regarding the formation of Iapetus, on its unusual orbit, complicates the interpretation of this result.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Dynamical Stability of Imaged Planetary Systems in Formation: Application to HL Tau

A recent ALMA image revealed several concentric gaps in the protoplanetary disk surrounding the young star HL Tau. We consider the hypothesis that these gaps are carved by planets, and present a general framework for understanding the dynamical stability of such systems over typical disk lifetimes, providing estimates for the maximum planetary masses. We collect these easily evaluated constraints into a workflow that can help guide the design and interpretation of new observational campaigns and numerical simulations of gap opening in such systems. We argue that the locations of resonances should be significantly shifted in massive disks like HL Tau, and that theoretical uncertainties in the exact offset, together with observational errors, imply a large uncertainty in the dynamical state and stability in such disks. This presents an important barrier to using systems like HL Tau as a proxy for the initial conditions following planet formation. An important observational avenue to breaking this degeneracy is to search for eccentric gaps, which could implicate resonantly interacting planets. Unfortunately, massive disks like HL Tau should induce swift pericenter precession that would smear out any such eccentric features of planetary origin. This motivates pushing toward more typical, less massive disks. For a nominal non-resonant model of the HL Tau system with five planets, we find a maximum mass for the outer three bodies of approximately 2 Neptune masses. In a resonant configuration, these planets can reach at least the mass of Saturn. The inner two planets' masses are unconstrained by dynamical stability arguments.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

WHFast: A fast and unbiased implementation of a symplectic Wisdom-Holman integrator for long term gravitational simulations

We present WHFast, a fast and accurate implementation of a Wisdom-Holman symplectic integrator for long-term orbit integrations of planetary systems. WHFast is significantly faster and conserves energy better than all other Wisdom-Holman integrators tested. We achieve this by significantly improving the Kepler-solver and ensuring numerical stability of coordinate transformations to and from Jacobi coordinates. These refinements allow us to remove the linear secular trend in the energy error that is present in other implementations. For small enough timesteps we achieve Brouwer's law, i.e. the energy error is dominated by an unbiased random walk due to floating-point round-off errors. We implement symplectic correctors up to order eleven that significantly reduce the energy error. We also implement a symplectic tangent map for the variational equations. This allows us to efficiently calculate two widely used chaos indicators the Lyapunov characteristic number (LCN) and the Mean Exponential Growth factor of Nearby Orbits (MEGNO). WHFast is freely available as a flexible C package, as a shared library, and as an easy-to-use python module.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

First observations of the Phoebe ring in optical light

The Phoebe ring, Saturn's largest and faintest ring, lies far beyond the planet's well-known main rings. It is primarily sourced by collisions with Saturn's largest irregular satellite Phoebe, perhaps through stochastic macroscopic collisions, or through more steady micrometeoroid bombardment. The ring was discovered with the Spitzer Space Telescope at 24 $μ$m and has a normal optical depth of $\sim 2 \times 10^{-8}$ (Verbiscer et al. 2009). We report the first observations of sunlight scattered off the Phoebe ring using the Cassini spacecraft's ISS camera at optical wavelengths. We find that material between $\approx 130-210$ Saturnian radii ($R_S$) from the planet produces an I/F of $1.7 \pm 0.1 \times 10^{-11}$ per $R_S$ of the line-of-sight distance through the disk. Combining our measurements with the Spitzer infrared data, we can place constraints on the ring-particles' light-scattering properties. Depending on the particles' assumed phase function, the derived single-scattering albedo can match either photometric models of Phoebe's dark regolith or brighter sub-surface material excavated by macroscopic impacts on Phoebe.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Chaotic dust dynamics and implications for the hemispherical color asymmetries of the Uranian satellites

Dust grains generated by the Uranian irregular satellites will undergo chaotic large-amplitude eccentricity oscillations under the simultaneous action of radiation forces and the highly misaligned quadrupole potentials of the oblate planet and distant Sun. From a suite of orbital histories, we estimate collision proba- bilities of dust particles with the regular satellites and argue that this process may explain the observed hemispherical color asymmetries of the outermost four regular satellites of Uranus.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Consequences of an Eccentric Orbit for Fomalhaut b

Fomalhaut b is currently the least massive, directly imaged exoplanet candidate. New observation epochs reveal this object to be moving on a highly eccentric orbit (Kalas et al. 2013), which sets important new constraints. I consider scenarios where Fomalhaut b is the only object interacting with the debris disk, and ones involving an additional unseen planet. I also investigate the possibility that Fomalhaut b is merely a transient dust cloud in light of the revised eccentric orbit. I argue that the scenario best able to match the observational constraints is a super-Earth Fomalhaut b surrounded by a vast cloud of dust that is generated by a population of irregular satellites, with an undetected $\sim$ Saturn-mass planet orbiting interior to the disk and driving the secular dynamics. Testable predictions are summarized that could differentiate between this scenario and other possibilities.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Finding the trigger to Iapetus' odd global albedo pattern: Dynamics of dust from Saturn's irregular satellites

The leading face of Saturn's moon Iapetus, Cassini Regio, has an albedo only one tenth that on its trailing side. The origin of this enigmatic dichotomy has been debated for over forty years, but with new data, a clearer picture is emerging. Motivated by Cassini radar and imaging observations, we investigate Soter's model of dark exogenous dust striking an originally brighter Iapetus by modeling the dynamics of the dark dust from the ring of the exterior retrograde satellite Phoebe under the relevant perturbations. In particular, we study the particles' probabilities of striking Iapetus, as well as their expected spatial distribution on the Iapetian surface. We find that, of the long-lived particles (greater than about 5 microns), most particle sizes (greater than about 10 microns) are virtually certain to strike Iapetus, and their calculated distribution on the surface matches up well with Cassini Regio's extent in its longitudinal span. The satellite's polar regions are observed to be bright, presumably because ice is deposited there. Thus, in the latitudinal direction we estimate polar dust deposition rates to help constrain models of thermal migration invoked to explain the bright poles (Spencer & Denk 2010). We also analyze dust originating from other irregular outer moons, determining that a significant fraction of that material will eventually coat Iapetus--perhaps explaining why the spectrum of Iapetus' dark material differs somewhat from that of Phoebe. Finally we track the dust particles that do not strike Iapetus, and find that most land on Titan, with a smaller fraction hitting Hyperion. As has been previously conjectured, such exogenous dust, coupled with Hyperion's chaotic rotation, could produce Hyperion's roughly isotropic, moderate-albedo surface.

0 citations0 reviewsNo code linkedOpen access