celmech: A Python package for celestial mechanics

Abstract

We present celmech, an open-source Python package designed to facilitate a wide variety of celestial mechanics calculations. The package allows users to formulate and integrate equations of motion incorporating user-specified terms from the classical disturbing function expansion of the interaction potential between pairs of planets. The code can be applied, for example, to isolate the contribution of particular resonances to a system's dynamical evolution and develop simple analytical models with the minimum number of terms required to capture a particular dynamical phenomenon. Equations and expressions can be easily manipulated by leveraging the extensive symbolic mathematics capabilities of the sympy Python package. The celmech package is designed to interface seamlessly with the popular $N$-body code REBOUND to facilitate comparisons between calculation results and direct $N$-body integrations. The code is extensively documented and numerous example Jupyter notebooks illustrating its use are available online.

submitted to AAS Journals, code and example notebooks available at http://github.com/shadden/celmech ; documentation available at http://celmech.readthedocs.io

References in corpus (38)

• SciPy 1.0--Fundamental Algorithms for Scientific Computing in Python
• REBOUND: An open-source multi-purpose N-body code for collisional dynamics
• batman: BAsic Transit Model cAlculatioN in Python
• 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
• EXOFAST: A fast exoplanetary fitting suite in IDL
• RadVel: The Radial Velocity Modeling Toolkit
• Extracting Planet Mass and Eccentricity From TTV data
• REBOUNDx: A Library for Adding Conservative and Dissipative Forces to Otherwise Symplectic N-body Integrations
• exoplanet: Gradient-based probabilistic inference for exoplanet data & other astronomical time series
• First order resonance overlap and the stability of close two planet systems
• TTVFast: An efficient and accurate code for transit timing inversion problems
• Three Body Resonance Overlap in Closely Spaced Multiple Planet Systems
• Analytical Treatment of Planetary Resonances
• Obliquity-Driven Sculpting of Exoplanetary Systems
• A Criterion for the Onset of Chaos in Systems of Two Eccentric Planets
• Transit timing to first order in eccentricity
• A general model of resonance capture in planetary systems: First and second order resonances
• The path to instability in compact multi-planetary systems
• Numerical and Analytical Modelling of Transit Time Variations
• Effects of Unseen Additional Planetary Perturbers on Compact Extrasolar Planetary Systems
• The Effect of Conjunctions on the Transit Timing Variations of Exoplanets
• Analytical model of multi-planetary resonant chains and constraints on migration scenarios
• An Integrable Model for the Dynamics of Planetary Mean Motion Resonances
• Secular spin-axis dynamics of exoplanets
• Transit timing variations for planets near eccentricity-type mean motion resonances
• A study of symplectic integrators for planetary system problems: error analysis and comparisons
• Stability of resonant configurations during the migration of planets and constraints on disk-planet interactions
• The origin of chaos in the Solar System through computer algebra
• Prospects for TTV Detection and Dynamical Constraints with TESS
• Transit Probabilities in Secularly Evolving Planetary Systems
• The Criterion for Chaos in Three-Planet Systems
• Non-Trivial Oblique Spin Equilibria of Super-Earths in Multi-planetary Systems
• Resonant Laplace-Lagrange theory for extrasolar systems in mean-motion resonance
• A differentiable N-body code for transit timing and dynamical modeling. I. Algorithm and derivatives
• Modeling Radial Velocity Data of Resonant Planets to Infer Migration Histories
• Long-term instability of the inner Solar System: numerical experiments
• A Criterion for the Onset of Chaos in Compact, Eccentric Multiplanet Systems

Cited by in corpus (5)

• Resonant and Ultra-short-period Planet Systems are at Opposite Ends of the Exoplanet Age Distribution
• The instability mechanism of compact multiplanet systems
• Higher-Order Mean-Motion Resonances Can Form in Type-I Disk Migration
• Secular Dynamics of Compact Three-Planet Systems
• General Relativity Can Prevent a Runaway Greenhouse on Potentially Habitable Planets Orbiting White Dwarfs