BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Hanno Rein

Hanno Rein contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
astro-ph.EPastro-ph.IMastro-ph.GAastro-ph.SRmath.DSphysics.comp-ph

Trust snapshot

Quick read

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

11 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
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

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