Related papers: Nii-body: Bayesian Inference of Multiplanet Dynami…
Transiting exoplanets in multi-planet systems exhibit non-Keplerian orbits as a result of the gravitational influence from companions which can cause the times and durations of transits to vary. The amplitude and periodicity of the transit…
We construct an advanced model for interacting multiple stellar systems in which we compute all trajectories with a numerical N-body integrator, namely the Bulirsch--Stoer from the SWIFT package. We can then derive various observables:…
Given that secular perturbations in a binary system not only excite high orbital eccentricities but also alter the planetary orbital inclination, the classical Keplerian orbital model is no longer applicable for orbital retrieval. The…
Here we present an open source Python-based Bayesian orbit retrieval code (Nii) that implements an automatic parallel tempering Markov chain Monte Carlo (APT-MCMC) strategy. Nii provides a module to simulate the observations of a…
While Keplerian orbits account for the majority of the astrometric motion of directly-imaged planets, perturbations due to N-body interactions allow us to directly constrain exoplanet masses in multiplanet systems. This has the potential to…
In the 20+ years of Doppler observations of stars, scientists have uncovered a diverse population of extrasolar multi-planet systems. A common technique for characterizing the orbital elements of these planets is Markov chain Monte Carlo…
We describe work in progress by a collaboration of astronomers and statisticians developing a suite of Bayesian data analysis tools for extrasolar planet (exoplanet) detection, planetary orbit estimation, and adaptive scheduling of…
We describe a technique for solving for the orbital elements of multiple planets from radial velocity (RV) and/or astrometric data taken with 1 m/s and microarcsecond precision, appropriate for efforts to detect Earth-massed planets in…
Multiple stellar systems are common especially among O and B stars. In order to accurately describe their dynamics, interactions among components must be accounted for. In this work, we describe the new dynamical model in Phoebe, which…
In a Keplerian system, a large number of bodies orbit a central mass. Accretion disks, protoplanetary disks, asteroid belts, and planetary rings are examples. Simulations of these systems require algorithms that are computationally…
We show that short-term perturbations among massive planets in multiple planet systems can result in radial velocity variations of the central star which differ substantially from velocity variations derived assuming the planets are…
When fitting N-body models to astronomical data - including transit times, radial velocity, and astrometric positions at observed times - the derivatives of the model outputs with respect to the initial conditions can help with model…
Cosmological N-body simulations are done on massively parallel computers. This necessitates the use of simple time integrators, and, additionally, of mesh-grid approximations of the potentials. Recently, Adamek et al. (2015);…
$N$-body integrations are used to model a wide range of astrophysical dynamics, but they suffer from errors which make their orbits diverge exponentially in time from the correct orbits. Over long time-scales, their reliability needs to be…
Entangled quantum many-body systems can be used as sensors that enable the estimation of parameters with a precision larger than that achievable with ensembles of individual quantum detectors. Typically, the parameter estimation strategy…
We present TES, a new n-body integration code for the accurate and rapid propagation of planetary systems in the presence of close encounters. TES builds upon the classic Encke method and integrates only the perturbations to Keplerian…
If an orbit is fitted from combined RV and astrometric data, the orbit should be physically consistent with both data sets. The Keplerian orbit of a planet is a highly nonlinear function of seven parameters. The astrometric orbit problem…
The $N$-body problem is of historical significance because it was the first implementation of the Newtonian dynamical laws for the description of our Solar System. Motivated by this, the project's goal is to revisit this problem for small…
Exoplanet transits contain substantial information about the architecture of a system. By fitting transit light curves we can extract dynamical parameters and place constraints on the properties of the planets and their host star. Having a…
Static structure models, which map mass-radius constraints to bulk planet composition, are frequently used to categorise exoplanets due to their computational efficiency and the high-level insight they offer into planetary properties.…