This code allows the simulation of an Intermediate Mass Ratio Inspiral (IMRI), where a small compact object -- the secondary -- is around an orbit of a larger body -- the primary, usually a massive black hole (MBH).
The system loses energy and angular momentum due to dissipative forces and slowly inspirals. The dissipative forces include gravitation wave (GW) emission, dynamical friction with dark matter (DM), interaction with an accretion disk, or stellar diffusion in stellar cusps.
The main purpose of the code is to solve the differential equations that arise in the modeling of IMRIs. The code can track the semimajor-axis, eccentricity, periapse angle, and inclination angle of the Kepler orbit. This orbit changes over time due to different dissipative forces. The code is very modular in design such that different dissipative forces can be added easily and their effects studied in detail. Additionally, the code can explore the phase space flows of the forces, compute the GW signal and compare it to LISA sensitivity.
The code also includes a stochastic description of some of the interactions with the help of stochastic differential equations (SDEs). This can be used to calculate rates and population behavior.
For most of the functionality look at the dissertation here.
This code has been used in our publications 2112.09586 and 2211.05145. See this and this file respectively for plot generation.
The code includes inspiral models from
9402014 - gravitational wave emission
1408.3534, 2204.12508, 2305.17281, 1711.09706 - DM halos + dynamical friction + accretion
1908.10241, 2107.00741, 1807.07163 - Keplerian orbits and waveforms
2002.12811,2108.04154 - Halo Feedback
2207.10086, 10.1086/324713, 10.1093/mnras/stac1294, 2205.10382 - Accretion disk profile + interactions
1508.01390 - Stellar Diffusion
See examples and dissertation folder.
Clone the repository and run
pip install -e .
(the -e option allows you to continuously edit the files without recompiling, don't use if you don't need to edit the files)
If you would like to use the Stochastic module, download the torchsde fork and install it.
If you use this in your publications, please cite the two papers 2112.09586, and 2211.05145, and the ASCL entry.
See LICENSE File