English

celmech: A Python package for celestial mechanics

Earth and Planetary Astrophysics 2022-10-26 v1 Instrumentation and Methods for Astrophysics

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 NN-body code REBOUND to facilitate comparisons between calculation results and direct NN-body integrations. The code is extensively documented and numerous example Jupyter notebooks illustrating its use are available online.

Keywords

Cite

@article{arxiv.2205.10385,
  title  = {celmech: A Python package for celestial mechanics},
  author = {Sam Hadden and Daniel Tamayo},
  journal= {arXiv preprint arXiv:2205.10385},
  year   = {2022}
}

Comments

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

R2 v1 2026-06-24T11:23:52.333Z