Nauyaca: a new tool to determine planetary masses and orbital elements through transit timing analysis
arXiv:2110.03830 · doi:10.3847/1538-3881/ac2744
Abstract
Transit Timing Variations (TTVs) is currently the most successful method to determine dynamical masses and orbital elements for Earth-sized transiting planets. Precise mass determination is fundamental to restrict planetary densities and thus infer planetary compositions. In this work, we present Nauyaca, a Python package dedicated to find planetary masses and orbital elements through the fitting of observed mid-transit times from a N-body approach. The fitting strategy consists in performing a sequence of minimization algorithms (optimizers) that are used to identify high probability regions in the parameter space. These results from optimizers are used for initialization of a Markov chain Monte Carlo (MCMC) method, using an adaptive Parallel-Tempering algorithm. A set of runs are performed in order to obtain posterior distributions of planetary masses and orbital elements. In order to test the tool, we created a mock catalog of synthetic planetary systems with different number of planets where all of them transit. We calculate their mid-transit times to give them as an input to Nauyaca, testing statistically its efficiency in recovering the planetary parameters from the catalog. For the recovered planets, we find typical dispersions around the real values of 1-14 M for masses, between 10-110 seconds for periods and between 0.01-0.03 for eccentricities. We also investigate the effects of the signal-to-noise and number of transits in the correct determination of the planetary parameters. Finally, we suggest choices of the parameters that govern the tool, for the usage with real planets, according to the complexity of the problem and computational facilities.
arXiv commentAccepted for publication in AJ, 23 pages, 10 figures, 7 tables
References in corpus (41)
- SciPy 1.0--Fundamental Algorithms for Scientific Computing in Python
- emcee: The MCMC Hammer
- Accurate masses and radii of normal stars: modern results and applications
- Parallel Tempering: Theory, Applications, and New Perspectives
- REBOUND: An open-source multi-purpose N-body code for collisional dynamics
- DHP Framework: Digital Health Passports Using Blockchain -- Use case on international tourism during the COVID-19 pandemic
- On detecting terrestrial planets with timing of giant planet transits
- Architecture of Kepler's Multi-transiting Systems: II. New investigations with twice as many candidates
- WHFast: A fast and unbiased implementation of a symplectic Wisdom-Holman integrator for long term gravitational simulations
- Extracting Planet Mass and Eccentricity From TTV data
- Data analysis recipes: Using Markov Chain Monte Carlo
- Dynamic temperature selection for parallel-tempering in Markov chain Monte Carlo simulations
- WASP-47: A Hot Jupiter System with Two Additional Planets Discovered by K2
- Kepler Mission Stellar and Instrument Noise Properties
- Architecture of Planetary Systems Based on Kepler Data: Number of Planets and Coplanarity
- Kepler Planet Masses and Eccentricities from TTV Analysis
- Very Low-Density Planets around Kepler-51 Revealed with Transit Timing Variations and an Anomaly Similar to a Planet-Planet Eclipse Event
- The Detection and Characterization of a Nontransiting Planet by Transit Timing Variations
- KOI-142, the King of Transit Variations, is a Pair of Planets near the 2:1 Resonance
- TTVFast: An efficient and accurate code for transit timing inversion problems
- Secure TTV Mass Measurements: Ten Kepler Exoplanets between 3 and 8 Earth Masses with Diverse Densities and Incident Fluxes
- The Planetary Mass-Radius Relation and its Dependence on Orbital Period as Measured by Transit Timing Variations and Radial Velocities
- Transit timing to first order in eccentricity
- TRADES: A new software to derive orbital parameters from observed transit times and radial velocities. Revisiting Kepler-11 and Kepler-9
- Numerical and Analytical Modelling of Transit Time Variations
- Sensitivity bias in the mass-radius distribution from Transit Timing Variations and Radial Velocity measurements
- The Effect of Conjunctions on the Transit Timing Variations of Exoplanets
- Systemic: A Testbed for Characterizing the Detection of Extrasolar Planets. II. Numerical approaches to the Transit Timing Inverse Problem
- Fast Inversion Method for Determination of Planetary Parameters from Transit Timing Variations
- Eccentric Companions to Kepler-448b and Kepler-693b: Clues to the Formation of Warm Jupiters
- Physical properties of the WASP-67 planetary system from multi-colour photometry
- Quantifying the challenges of detecting unseen planetary companions with transit timing variations
- Testing a Fast Dynamical Indicator: The MEGNO
- A Simplified Photodynamical Model for Planetary Mass Determination in Low-Eccentricity Multi-Transiting Systems
- Modal Decomposition of TTV - Inferring Planet Masses and Eccentricities
- Masses of Kepler-46b, c from Transit Timing Variations
- The Efficiency of Geometric Samplers for Exoplanet Transit Timing Variation Models
- A super-Earth and a mini-Neptune around Kepler-59
- Masses of the Kepler-419 Planets from Transit Timing Variations Analysis
- How to Find a Planet from Transit Variations
- An alternative stable solution for the Kepler-419 system, obtained with the use of a genetic algorithm