A direct N-body integrator for modelling the chaotic, tidal dynamics of multi-body extrasolar systems: TIDYMESS
arXiv:2209.03955 · doi:10.1093/mnras/stad1133
Abstract
Tidal dissipation plays an important role in the dynamical evolution of moons, planets, stars and compact remnants. The interesting complexity originates from the interplay between the internal structure and external tidal forcing. Recent and upcoming observing missions of exoplanets and stars in the Galaxy help to provide constraints on the physics of tidal dissipation. It is timely to develop new N-body codes, which allow for experimentation with various tidal models and numerical implementations. We present the open-source N-body code TIDYMESS, which stands for ``TIdal DYnamics of Multi-body ExtraSolar Systems''. This code implements a creep deformation law for the bodies, parametrized by their fluid Love numbers and fluid relaxation times. Due to tidal and centrifugal deformations, we approximate the general shape of a body to be an ellipsoid. We calculate the associated gravitational field to quadruple order, from which we derive the gravitational accelerations and torques. The equations of motion for the orbits, spins and deformations are integrated directly using a fourth-order integration method based on a symplectic composition. We implement a novel integration method for the deformations, which allows for a time step solely dependent on the orbits, and not on the spin periods or fluid relaxation times. This feature greatly speeds up the calculations, while also improving the consistency when comparing different tidal regimes. We demonstrate the capabilities and performance of TIDYMESS, particularly in the niche regime of parameter space where orbits are chaotic and tides become non-linear.
17 pages, 6 figures
References in corpus (18)
- On the tidal evolution of Hot Jupiters on inclined orbits
- The triple evolution dynamical instability: Stellar collisions in the field and the formation of exotic binaries
- Tidal evolution of hierarchical and inclined systems
- Deformation and tidal evolution of close-in planets and satellites using a Maxwell viscoelastic rheology
- Stellar triples on the edge; Comprehensive overview of the evolution of destabilised triples leading to stellar and binary exotica
- Traditional formation scenarios fail to explain 4:3 mean motion resonances
- The fate of planetesimal discs in young open clusters: implications for 1I/'Oumuamua, the Kuiper belt, the Oort cloud and more
- Complete spin and orbital evolution of close-in bodies using a Maxwell viscoelastic rheology
- Gargantuan chaotic gravitational three-body systems and their irreversibility to the Planck length
- Stellar collisions in flattened and rotating Pop. III star clusters
- Chaos in self-gravitating many-body systems: Lyapunov time dependence of and the influence of general relativity
- Towards time symmetric N-body integration
- The Origin of Chaos in the Orbit of Comet 1P/Halley
- Tidal excitation of the obliquity of Earth-like planets in the habitable zone of M-dwarf stars
- GRIT: a package for structure-preserving simulations of gravitationally interacting rigid-bodies
- Tidal spin down rates of homogeneous triaxial viscoelastic bodies
- Modeling gravitational few-body problems with TSUNAMI and OKINAMI
- Reversible time-step adaptation for the integration of few-body systems