Modeling gravitational few-body problems with TSUNAMI and OKINAMI

TL;DR

This paper introduces two numerical codes, TSUNAMI and OKINAMI, for modeling gravitational few-body systems, addressing recent astrophysical challenges like exoplanets and gravitational waves.

Contribution

The paper presents novel computational tools for simulating three- and few-body gravitational systems with advanced features and applications.

Findings

TSUNAMI incorporates regularization, tidal forces, and post-Newtonian corrections.

OKINAMI efficiently models stable hierarchical triples using secular, double-averaged methods.

Applications include gravitational-wave astronomy, planetary science, and escape dynamics.

Abstract

In recent years, an increasing amount of attention is being paid to the gravitational few-body problem and its applications to astrophysical scenarios. Among the main reasons for this renewed interest there is large number of newly discovered exoplanets and the detection of gravitational waves. Here, we present two numerical codes to model three- and few-body systems, called TSUNAMI and OKINAMI. The TSUNAMI code is a direct few-body code with algorithmic regularization, tidal forces and post-Newtonian corrections. OKINAMI is a secular, double-averaged code for stable hierarchical triples. We describe the main methods implemented in our codes, and review our recent results and applications to gravitational-wave astronomy, planetary science and statistical escape theories.