Researcher profile

Daniel Tamayo

Daniel Tamayo contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
7works
0followers
3topics
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)

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.

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.

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&#39; near-integer period ratios. We incorporate recent transit and radial velocity observations of the V1298 Tau system, and investigate constraints on the system&#39;s orbital architecture imposed by requiring dynamical stability on timescales much shorter than the system&#39;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.

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&#39; accuracy and efficiency. These additional effects are publicly available as of REBOUNDx&#39;s latest release.

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&#39; 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.

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 &#34;errors&#34; 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.

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&#39; 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.