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

Hanno Rein

Hanno Rein 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.SRastro-ph.GAphysics.comp-phmath.DSmath.NAmath.CAnlin.CDphysics.data-anphysics.flu-dynphysics.soc-ph

Catalog footprint

What is connected

34works
12topics
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

34 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
preprint2024arXiv

A new timestep criterion for N-body simulations

We derive a new criterion for estimating characteristic dynamical timescales in N-body simulations. The criterion uses the second, third, and fourth derivatives of particle positions: acceleration, jerk, and snap. It can be used for choosing timesteps in integrators with adaptive step size control. For any two-body problem the criterion is guaranteed to determine the orbital period and pericenter timescale regardless of eccentricity. We discuss why our criterion is the simplest derivative-based expression for choosing adaptive timesteps with the above properties and show its superior performance over existing criteria in numerical tests. Because our criterion uses lower order derivatives, it is less susceptible to rounding errors caused by finite floating point precision. This significantly decreases the volume of phase space where an adaptive integrator fails or gets stuck due to unphysical timestep estimates. For example, our new criterion can accurately estimate timesteps for orbits around a 50m sized Solar System object located at 40AU from the coordinate origin when using double floating point precision. Previous methods where limited to objects larger than 10km. We implement our new criterion in the high order IAS15 integrator which is part of the freely available N-body package REBOUND.

0 citations0 reviewsNo code linkedOpen access
preprint2023arXiv

Thermal hysteresis and front propagation in dense planetary rings

Saturn's rings are composed of icy grains, most in the mm to m size ranges, undergoing several collisions per orbit. Their collective behaviour generates a remarkable array of structure over many orders of magnitude, much of it not well understood. On the other hand, the collisional properties and parameters of individual ring particles are poorly constrained; usually N-body simulations and kinetic theory employ hard-sphere models with a coefficient of restitution $ε$ that is constant or a decreasing function of impact speed. Due to plastic deformation of surface regolith, however, it is likely that $ε$ will be more complicated, at the very least a non-monotonic function. We undertake N-body simulations with the REBOUND code with non-monotonic $ε$ laws to approximate surfaces that are friable but not sticking. Our simulations reveal that such ring models can support two thermally stable steady states for the same (dynamical) optical depth: a cold and a warm state. If the ring breaks up into radial bands of one or the other state, we find that warmer states tend to migrate into the colder states via a coherent travelling front. We also find stationary `viscous' fronts, which connect states of different optical depth, but the same angular momentum flux. We discuss these preliminary results and speculate on their implications for structure formation in Saturn's B and C-rings, especially with respect to structures that appear in Cassini images but not in occultations.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

An Implementation of Stochastic Forces for the N-body code REBOUND

We describe the implementation of a new module which can be used to simulate physical systems in which the motion of particles is affected by stochastic forces. Such forces are expected to be present in turbulent circumstellar disks or remnant planetesimal disks. Our implementation offers a convenient way to generate correlated noise with a user-specified amplitude and auto-correlation time for each particle. The module has minimal memory requirements and is freely available within the REBOUNDx additional effects library.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

On the long-term stability of the Solar System in the presence of weak perturbations from stellar flybys

The architecture and evolution of planetary systems are shaped in part by stellar flybys. Within this context, we look at stellar encounters which are too weak to immediately destabilize a planetary system but are nevertheless strong enough to measurably perturb the system's dynamical state. We estimate the strength of such perturbations on secularly evolving systems using a simple analytic model and confirm those estimates with direct N-body simulations. We then run long-term integrations and show that even small perturbations from stellar flybys can influence the stability of planetary systems over their lifetime. We find that small perturbations to the outer planets' orbits are transferred between planets, increasing the likelihood that the inner planetary system will destabilize. Specifically, our results for the Solar System show that relative perturbations to Neptune's semi-major axis of order 0.1% are strong enough to increase the probability of destabilizing the Solar System within 5 Gyrs by one order of magnitude.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The vertical structure of debris disks and the impact of gas

The vertical structure of debris disks provides clues about their dynamical evolution and the collision rate of the unseen planetesimals. Thanks to the ever-increasing angular resolution of contemporary instruments and facilities, we are beginning to constrain the scale height of a handful of debris disks, either at near-infrared or millimeter wavelengths. Nonetheless, this is often done for individual targets only. We present here the geometric modeling of eight disks close to edge-on, all observed with the same instrument (SPHERE) and using the same mode (dual-beam polarimetric imaging). Motivated by the presence of CO gas in two out of the eight disks, we then investigate the impact that gas can have on the scale height by performing N-body simulations including gas drag and collisions. We show that gas can quickly alter the dynamics of particles (both in the radial and vertical directions), otherwise governed by gravity and radiation pressure. We find that, in the presence of gas, particles smaller than a few tens of microns can efficiently settle toward the midplane at the same time as they migrate outward beyond the birth ring. For second generation gas ($M_\mathrm{gas} \leq 0.1$ $M_\oplus$), the vertical settling should be best observed in scattered light images compared to observations at millimeter wavelengths. But if the gas has a primordial origin ($M_\mathrm{gas} \geq 1$ $M_\oplus$), the disk will appear very flat both at near-infrared and sub-mm wavelengths. Finally, far beyond the birth ring, our results suggest that the surface brightness profile can be as shallow as $\sim -2.25$.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

BEBOP III. Observations and an independent mass measurement of Kepler-16 (AB) b -- the first circumbinary planet detected with radial velocities

The radial velocity method is amongst the most robust and most established means of detecting exoplanets. Yet, it has so far failed to detect circumbinary planets despite their relatively high occurrence rates. Here, we report velocimetric measurements of Kepler-16A, obtained with the SOPHIE spectrograph, at the Observatoire de Haute-Provence's 193cm telescope, collected during the BEBOP survey for circumbinary planets. Our measurements mark the first radial velocity detection of a circumbinary planet, independently determining the mass of Kepler-16~(AB)~b to be $0.313 \pm 0.039\,{\rm M}_{\rm Jup}$, a value in agreement with eclipse timing variations. Our observations demonstrate the capability to achieve photon-noise precision and accuracy on single-lined binaries, with our final precision reaching $\rm 1.5~m\,s^{-1}$ on the binary and planetary signals. Our analysis paves the way for more circumbinary planet detections using radial velocities which will increase the relatively small sample of currently known systems to statistically relevant numbers, using a method that also provides weaker detection biases. Our data also contain a long-term radial velocity signal, which we associate with the magnetic cycle of the primary star.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Embedded operator splitting methods for perturbed systems

It is common in classical mechanics to encounter systems whose Hamiltonian $H$ is the sum of an often exactly integrable Hamiltonian $H_0$ and a small perturbation $εH_1$ with $ε\ll1$. Such near-integrability can be exploited to construct particularly accurate operator splitting methods to solve the equations of motion of $H$. However, in many cases, for example in problems related to planetary motion, it is computationally expensive to obtain the exact solution to $H_0$. In this paper we present a new family of embedded operator splitting (EOS) methods which do not use the exact solution to $H_0$, but rather approximate it with yet another, embedded operator splitting method. Our new methods have all the desirable properties of classical methods which solve $H_0$ directly. But in addition they are very easy to implement and in some cases faster. When applied to the problem of planetary motion, our EOS methods have error scalings identical to that of the often used Wisdom-Holman method but do not require a Kepler solver, nor any coordinate transformations, or the allocation of memory. The only two problem specific functions that need to be implemented are the straight-forward kick and drift steps typically used in the standard second order leap-frog method.

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

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

The Astropy Problem

The Astropy Project (http://astropy.org) is, in its own words, "a community effort to develop a single core package for Astronomy in Python and foster interoperability between Python astronomy packages." For five years this project has been managed, written, and operated as a grassroots, self-organized, almost entirely volunteer effort while the software is used by the majority of the astronomical community. Despite this, the project has always been and remains to this day effectively unfunded. Further, contributors receive little or no formal recognition for creating and supporting what is now critical software. This paper explores the problem in detail, outlines possible solutions to correct this, and presents a few suggestions on how to address the sustainability of general purpose astronomical software.

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

High Order Harmonics in Light Curves of Kepler Planets

The Kepler mission was launched in 2009 and has discovered thousands of planet candidates. In a recent paper, Esteves et al. (2013) found a periodic signal in the light curves of KOI-13 and HAT-P-7, with a frequency triple the orbital frequency of a transiting planet. We found similar harmonics in many systems with a high occurrence rate. At this time, the origins of the signal are not entirely certain. We look carefully at the possibility of errors being introduced through our data processing routines but conclude that the signal is real. The harmonics on multiples of the orbital frequency are a result of non-sinusoidal periodic signals. We speculate on their origin and generally caution that these harmonics could lead to wrong estimates of planet albedos, beaming mass estimates, and ellipsoidal variations.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Reanalysis of radial velocity data from the resonant planetary system HD128311

The multi-planetary system HD128311 hosts at least two planets. Its dynamical formation history has been studied extensively in the literature. We reanalyse the latest radial velocity data for this system with the affine-invariant Markov chain Monte Carlo sampler EMCEE. Using the high order integrator IAS15, we perform a fully dynamical fit, allowing the planets to interact during the sampling process. A stability analysis using the MEGNO indicator reveals that the system is located in a stable island of the parameter space. In contrast to a previous study, we find that the system is locked in a 2:1 mean motion resonance. The resonant angle $φ_1$ is librating with a libration amplitude of approximately 37°. The existence of mean motion resonances has important implication for planet formation theories. Our results confirm predictions of models involving planet migration and stochastic forces.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Tides Alone Cannot Explain Kepler Planets Close to 2:1 MMR

A number of Kepler planet pairs lie just wide of first-order mean motion resonances (MMRs). Tides have been frequently proposed to explain these pileups, but it is still an ongoing discussion. We contribute to this discussion by calculating an optimistic theoretical estimate on the minimum initial eccentricity required by Kepler planets to explain the current observed spacing, and compliment these calculations with N-body simulations. In particular, we investigate 27 Kepler systems having planets within 6% of the 2:1 MMR, and find that the initial eccentricities required to explain the observed spacings are unreasonable from simple dynamical arguments. Furthermore, our numerical simulations reveal resonant tugging, an effect which conspires against the migration of resonant planets away from the 2:1 MMR, requiring even higher initial eccentricities in order to explain the current Kepler distribution. Overall, we find that tides alone cannot explain planets close to 2:1 MMR, and additional mechanisms are required to explain these systems.

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

IAS15: A fast, adaptive, high-order integrator for gravitational dynamics, accurate to machine precision over a billion orbits

We present IAS15, a 15th-order integrator to simulate gravitational dynamics. The integrator is based on a Gauß-Radau quadrature and can handle conservative as well as non-conservative forces. We develop a step-size control that can automatically choose an optimal timestep. The algorithm can handle close encounters and high-eccentricity orbits. The systematic errors are kept well below machine precision and long-term orbit integrations over $10^9$ orbits show that IAS15 is optimal in the sense that it follows Brouwer's law, i.e. the energy error behaves like a random walk. Our tests show that IAS15 is superior to a mixed-variable symplectic integrator (MVS) and other popular integrators, including high-order ones, in both speed and accuracy. In fact, IAS15 preserves the symplecticity of Hamiltonian systems better than the commonly-used nominally symplectic integrators to which we compared it. We provide an open-source implementation of IAS15. The package comes with several easy-to-extend examples involving resonant planetary systems, Kozai-Lidov cycles, close encounters, radiation pressure, quadrupole moment, and generic damping functions that can, among other things, be used to simulate planet-disc interactions. Other non-conservative forces can be added easily.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Some inconvenient truths about biosignatures involving two chemical species on Earth-like exoplanets

The detection of strong thermochemical disequilibrium in the atmosphere of an extrasolar planet is thought to be a potential biosignature. In this article we present a new kind of false positive that can mimic a disequilibrium or any other biosignature that involves two chemical species. We consider a scenario where the exoplanet hosts a moon that has its own atmosphere and neither of the atmospheres is in chemical disequilibrium. Our results show that the integrated spectrum of the planet and the moon closely resembles that of a single object in strong chemical disequilibrium. We derive a firm limit on the maximum spectral resolution that can be obtained for both directly-imaged and transiting planets. The spectral resolution of even idealized space-based spectrographs that might be achievable in the next several decades is in general insufficient to break the degeneracy. Both chemical species can only be definitively confirmed in the same object if absorption features of both chemicals can be unambiguously identified and their combined depth exceeds 100%.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Large-scale N-body simulations of the viscous overstability in Saturn's rings

We present results from large-scale particle simulations of the viscous overstability in Saturn's rings. The overstability generates a variety of structure on scales covering a few hundred metres to several kilometres, including axisymmetric wavetrains and larger-scale modulations. Such patterns have been observed in Saturn's rings by the Cassini spacecraft. Our simulations model the collisional evolution of particles in a co-rotating patch of the disk. These are the largest N-body simulations of the viscous overstability yet performed. The radial box size is five orders of magnitude larger than a typical particle radius, and so describes a 20-50 km radial portion of the rings. Its evolution is tracked for more than 10,000 orbits. In agreement with hydrodynamics, our N-body simulations reveal that the viscous overstability exhibits a rich set of dynamics characterised by nonlinear travelling waves with wavelengths of a few hundred meters. In addition, wave defects, such as sources and shocks, punctuate this bed of waves and break them up into large-scale divisions of radial width ~5 km. We find that the wavelength of the travelling waves is positively correlated with the mean optical depth. In order to assess the role of the numerical boundary conditions and also background ring structure, we include simulations of broad spreading rings and simulations with a gradient in the background surface density. Overall, our numerical results and approach provide a tool with which to interpret Cassini occultation observations of microstructure in Saturn's rings. We present an example of such a synthetic occultation observation and discuss what features to expect. We also make the entire source code freely available.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

SMACK: A New Algorithm for Modeling Collisions and Dynamics of Planetesimals in Debris Disks

We present the Superparticle Model/Algorithm for Collisions in Kuiper belts and debris disks (SMACK), a new method for simultaneously modeling, in 3-D, the collisional and dynamical evolution of planetesimals in a debris disk with planets. SMACK can simulate azimuthal asymmetries and how these asymmetries evolve over time. We show that SMACK is stable to numerical viscosity and numerical heating over 10^7 yr, and that it can reproduce analytic models of disk evolution. We use SMACK to model the evolution of a debris ring containing a planet on an eccentric orbit. Differential precession creates a spiral structure as the ring evolves, but collisions subsequently break up the spiral, leaving a narrower eccentric ring.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The formation of systems with closely spaced low-mass planets and the application to Kepler-36

The Kepler-36 system consists of two planets that are spaced unusually close together, near the 7:6 mean motion resonance. While it is known that mean motion resonances can easily form by convergent migration, Kepler-36 is an extreme case due to the close spacing and the relatively high planet masses of 4 and 8 times that of the Earth. In this paper, we investigate whether such a system can be obtained by interactions with the protoplanetary disc. These discs are thought to be turbulent and exhibit density fluctuations which might originate from the magneto-rotational instability. We adopt a realistic description for stochastic forces due to these density fluctuations and perform both long term hydrodynamical and N-body simulations. Our results show that planets in the Kepler-36 mass range can be naturally assembled into a closely spaced planetary system for a wide range of migration parameters in a turbulent disc similar to the minimum mass solar nebula. The final orbits of our formation scenarios tend to be Lagrange stable, even though large parts of the parameter space are chaotic and unstable.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A proposal for community driven and decentralized astronomical databases and the Open Exoplanet Catalogue

I present a new kind of astronomical database based on small text files and a distributed version control system. This encourages the community to work collaboratively. It creates a decentralized, completely open and democratic way of managing small to medium sized heterogeneous astronomical databases and catalogues. The use of the XML file format allows an easy to parse and read, yet dynamic and extendable database structure. The Open Exoplanet Catalogue is based on these principles and presented as an example. It is a catalogue of all discovered extra-solar planets. It is the only catalogue that can correctly represent the orbital structure of planets in arbitrary binary, triple and quadruple star systems, as well as orphan planets.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Period ratios in multi-planetary systems discovered by Kepler are consistent with planet migration

The Kepler planet candidates are an interesting testbed for planet formation scenarios. We present results from N-body simulations of multi-planetary systems that resemble those observed by Kepler. We add both smooth (Type I/II) and stochastic migration forces. The observed period ratio distribution is inconsistent with either of those two scenarios on its own. However, applying both stochastic and smooth migration forces to the planets simultaneously results in a period ratio distribution that is similar to the observed one. This is a natural scenario if planets form in a turbulent proto-planetary disk where these forces are always present. We show how the observed period ratio and eccentricity distribution can constrain the relative strength of these forces, a parameter which has been notoriously hard to predict for decades. We make the source code of our simulations and the initial conditions freely available to enable the community to expand this study and include effect other than planetary migration.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Planet-disk interaction in highly inclined systems

We study the interaction of a proto-planetary disk and a planet on a highly inclined orbit in the linear regime. The evolution of the planet is dominated by dynamical friction for planet masses above several Earth-masses. Smaller planets are dominated by aerodynamic drag, especially for very high inclinations and retrograde orbits. The time-scales associated with migration and inclination damping are calculated. For certain values of the inclination, the inclination damping time-scale is longer than the migration time-scale and the disk lifetime. This result shows that highly inclined planets can not (re-)align with the proto-planetary disk. We discuss the dependence of numerical simulations on the gravitational softening parameter. We find only a logarithmic dependence, making global three dimensional simulations of this process computationally feasible. A large fraction of Hot Jupiters is on highly inclined orbits with respect to the rotation axis of the star. On the other hand small-mass planetary systems discovered by the Kepler mission have low mutual inclinations. This shows that there are two distinct formation mechanisms at work. The process that creates inclined Hot Jupiters does not operate on small mass planets because the damping timescales are so long that these systems would still be inclined today.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Stochastic flights of propellers

Kilometer-sized moonlets in Saturn's A ring create S-shaped wakes called "propellers" in surrounding material. The Cassini spacecraft has tracked the motions of propellers for several years and finds that they deviate from Keplerian orbits having constant semimajor axes. The inferred orbital migration is known to switch sign. We show using a statistical test that the time series of orbital longitudes of the propeller Blériot is consistent with that of a time-integrated Gaussian random walk. That is, Blériot's observed migration pattern is consistent with being stochastic. We further show, using a combination of analytic estimates and collisional N-body simulations, that stochastic migration of the right magnitude to explain the Cassini observations can be driven by encounters with ring particles 10-20 m in radius. That the local ring mass is concentrated in decameter-sized particles is supported on independent grounds by occultation analyses.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The gravitational instability of a stream of co-orbital particles

We describe the dynamics of a stream of equally spaced macroscopic particles in orbit around a central body (e.g. a planet or star). A co-orbital configuration of small bodies may be subject to gravitational instability, which takes the system to a spreading, disordered and collisional state. We detail the linear instability's mathematical and physical features using the shearing sheet model and subsequently track its nonlinear evolution with local N-body simulations. This model provides a convenient tool with which to understand the gravitational and collisional dynamics of narrow belts, such as Saturn's F-ring and the streams of material wrenched from tidally disrupted bodies. In particular, we study the tendency of these systems to form long-lived particle aggregates. Finally, we uncover an unexpected connection between the linear dynamics of the gravitational instability and the magnetorotational instability.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Traditional formation scenarios fail to explain 4:3 mean motion resonances

At least two multi-planetary systems in a 4:3 mean motion resonance have been found by radial velocity surveys. These planets are gas giants and the systems are only stable when protected by a resonance. Additionally the Kepler mission has detected at least 4 strong candidate planetary systems with a period ratio close to 4:3. This paper investigates traditional dynamical scenarios for the formation of these systems. We systematically study migration scenarios with both N-body and hydro-dynamic simulations. We investigate scenarios involving the in-situ formation of two planets in resonance. We look at the results from finely tuned planet-planet scattering simulations with gas disk damping. Finally, we investigate a formation scenario involving isolation-mass embryos. Although the combined planet-planet scattering and damping scenario seems promising, none of the above scenarios is successful in forming enough systems in 4:3 resonance with planetary masses similar to the observed ones. This is a negative result but it has important implications for planet formation. Previous studies were successful in forming 2:1 and 3:2 resonances. This is generally believed to be evidence of planet migration. We highlight the main differences between those studies and our failure in forming a 4:3 resonance. We also speculate on more exotic and complicated ideas. These results will guide future investigators toward exploring the above scenarios and alternative mechanisms in a more general framework.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

REBOUND: An open-source multi-purpose N-body code for collisional dynamics

REBOUND is a new multi-purpose N-body code which is freely available under an open-source license. It was designed for collisional dynamics such as planetary rings but can also solve the classical N-body problem. It is highly modular and can be customized easily to work on a wide variety of different problems in astrophysics and beyond. REBOUND comes with three symplectic integrators: leap-frog, the symplectic epicycle integrator (SEI) and a Wisdom-Holman mapping (WH). It supports open, periodic and shearing-sheet boundary conditions. REBOUND can use a Barnes-Hut tree to calculate both self-gravity and collisions. These modules are fully parallelized with MPI as well as OpenMP. The former makes use of a static domain decomposition and a distributed essential tree. Two new collision detection modules based on a plane-sweep algorithm are also implemented. The performance of the plane-sweep algorithm is superior to a tree code for simulations in which one dimension is much longer than the other two and in simulations which are quasi-two dimensional with less than one million particles. In this work, we discuss the different algorithms implemented in REBOUND, the philosophy behind the code's structure as well as implementation specific details of the different modules. We present results of accuracy and scaling tests which show that the code can run efficiently on both desktop machines and large computing clusters.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Symplectic integrators in the shearing sheet

The shearing sheet is a model dynamical system that is used to study the small-scale dynamics of astrophysical disks. Numerical simulations of particle trajectories in the shearing sheet usually employ the leapfrog integrator, but this integrator performs poorly because of velocity-dependent (Coriolis) forces. We describe two new integrators for this purpose; both are symplectic, time-reversible and second-order accurate, and can easily be generalized to higher orders. Moreover, both integrators are exact when there are no small-scale forces such as mutual gravitational forces between disk particles. In numerical experiments these integrators have errors that are often several orders of magnitude smaller than competing methods. The first of our new integrators (SEI) is well-suited for disks in which the typical inter-particle separation is large compared to the particles' Hill radii (e.g., planetary rings), and the second (SEKI) is designed for disks in which the particles are on bound orbits or the separation is smaller than the Hill radius (e.g., irregular satellites of the giant planets).

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Stochastic orbital migration of small bodies in Saturn's rings

Many small moonlets, creating propeller structures, have been found in Saturn's rings by the Cassini spacecraft. We study the dynamical evolution of such 20-50m sized bodies which are embedded in Saturn's rings. We estimate the importance of various interaction processes with the ring particles on the moonlet's eccentricity and semi-major axis analytically. For low ring surface densities, the main effects on the evolution of the eccentricity and the semi-major axis are found to be due to collisions and the gravitational interaction with particles in the vicinity of the moonlet. For large surface densities, the gravitational interaction with self-gravitating wakes becomes important. We also perform realistic three dimensional, collisional N-body simulations with up to a quarter of a million particles. A new set of pseudo shear periodic boundary conditions is used which reduces the computational costs by an order of magnitude compared to previous studies. Our analytic estimates are confirmed to within a factor of two. On short timescales the evolution is always dominated by stochastic effects caused by collisions and gravitational interaction with self-gravitating ring particles. These result in a random walk of the moonlet's semi-major axis. The eccentricity of the moonlet quickly reaches an equilibrium value due to collisional damping. The average change in semi-major axis of the moonlet after 100 orbital periods is 10-100m. This translates to an offset in the azimuthal direction of several hundred kilometres. We expect that such a shift is easily observable.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The effects of stochastic forces on the evolution of planetary systems and Saturn's rings

The increasing number of extra-solar planets opens a new opportunity for studies of the formation of planetary systems. Resonant systems are of particular interest because their dynamical configuration provides constraints on the unobservable formation and migration phase. In this thesis, formation scenarios for the planetary systems HD128311 and HD45364 are presented. N-body simulations of two planets and two dimensional hydrodynamical simulations of proto-planetary discs are used to model the migration phase and the capture into resonance. The results indicate that the proto-planetary disc has a larger surface density than previously thought. Proto-planets are exposed to stochastic forces, generated by density fluctuations in the disc. A generic model of both a single planet, and two planets in resonance, being stochastically forced is presented. The system GJ876, for example, is stable for reasonable strengths of the stochastic forces, but systems with lighter planets can get disrupted. Even if they are not disrupted completely, stochastic forces create characteristic, observable libration patterns. Turbulence plays also an important role during the early phases of the planet formation process. Due to the large separation of scales, the gravitational collapse in the core accretion model is very hard to model numerically. A scaled method is presented, that allows for the correct treatment of self-gravity in a marginally collisional system by taking into account the relevant small scale processes. Interestingly, this system is dynamically very similar to Saturn's rings. The stochastic migration of small bodies in Saturn's rings is also studied. Analytic predictions of the interactions of a moonlet/propeller with ring particles are compared to collisional N-body simulations with up to a million particles. The random walk is fast enough to be directly observable by the Cassini spacecraft.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The Validity of the Super-Particle Approximation during Planetesimal Formation

The formation mechanism of planetesimals in protoplanetary discs is hotly debated. Currently, the favoured model involves the accumulation of meter-sized objects within a turbulent disc, followed by a phase of gravitational instability. At best one can simulate a few million particles numerically as opposed to the several trillion meter-sized particles expected in a real protoplanetary disc. Therefore, single particles are often used as super-particles to represent a distribution of many smaller particles. It is assumed that small scale phenomena do not play a role and particle collisions are not modeled. The super-particle approximation can only be valid in a collisionless or strongly collisional system, however, in many recent numerical simulations this is not the case. In this work we present new results from numerical simulations of planetesimal formation via gravitational instability. A scaled system is studied that does not require the use of super-particles. We find that the scaled particles can be used to model the initial phases of clumping if the properties of the scaled particles are chosen such that all important timescales in the system are equivalent to what is expected in a real protoplanetary disc. Constraints are given for the number of particles needed in order to achieve numerical convergence. We compare this new method to the standard super-particle approach. We find that the super-particle approach produces unreliable results that depend on artifacts such as the gravitational softening in both the requirement for gravitational collapse and the resulting clump statistics. Our results show that short range interactions (collisions) have to be modelled properly.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Mean-Motion Resonances of High Order in Extrasolar Planetary Systems

Many multi-planet systems have been discovered in recent years. Some of them are in mean-motion resonances (MMR). Planet formation theory was successful in explaining the formation of 2:1, 3:1 and other low resonances as a result of convergent migration. However, higher order resonances require high initial orbital eccentricities in order to be formed by this process and these are in general unexpected in a dissipative disk. We present a way of generating large initial eccentricities using additional planets. This procedure allows us to form high order MMRs and predict new planets using a genetic N-body code.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

The Dynamical Origin of the Multi-Planetary System HD45364

The recently discovered planetary system HD45364 which consists of a Jupiter and Saturn mass planet is very likely in a 3:2 mean motion resonance. The standard scenario to form planetary commensurabilities is convergent migration of two planets embedded in a protoplanetary disc. When the planets are initially separated by a period ratio larger than two, convergent migration will most likely lead to a very stable 2:1 resonance for moderate migration rates. To avoid this fate, formation of the planets close enough to prevent this resonance may be proposed. However, such a simultaneous formation of the planets within a small annulus, seems to be very unlikely. Rapid type III migration of the outer planet crossing the 2:1 resonance is one possible way around this problem. In this paper, we investigate this idea in detail. We present an estimate for the required convergent migration rate and confirm this with N-body and hydrodynamical simulations. If the dynamical history of the planetary system had a phase of rapid inward migration that forms a resonant configuration, we predict that the orbital parameters of the two planets are always very similar and hence should show evidence of that. We use the orbital parameters from our simulation to calculate a radial velocity curve and compare it to observations. Our model can explain the observational data as good as the previously reported fit. The eccentricities of both planets are considerably smaller and the libration pattern is different. Within a few years, it will be possible to observe the planet-planet interaction directly and thus distinguish between these different dynamical states.

0 citations0 reviewsNo code linkedOpen access