Make binary stellar evolution dynamical again!
Language: python script for binary evolution with N-body dynamics
version: 1.e-6
A collaborative effort to construct a binary population synthesis code in which the orbital evolution is calculated using a general-purpose N-body code and the stellar evolution by a general purpose stellar evolution code.
The objective is to eventually, be able to resplace the binary population synthesis modeules in N-body codes by more general-purpose (but slower) prescription for the binary evolution in which the dynamics is taken care of by the N-body code.
Initial discussions started on 27 January 2020 between: Nicola Giacobbo, Iorio Giuliano, Simon Portegies Zwart, Steven Rieder, Alessandro Trani, and Long Wang
Current content are 3 scripts:
-dynbin_simple.py: Simplest form of integrating a binary without any mass loss.
-dynbin_massloss.py: Simplest for of binary evolution in which the binary orbit is integrated with a 4th-order Hermite scheme. Stellar mass loss (currently constant with time) is directly incorporated into the N-body code via a channel.
-dynbin_massloss_bridge.py: Binary evolution in which the binary orbit is integrated with a 4th-order Hermite scheme and the stellar mass loss (currently constant with time) is incorporated using the drift-operator in the classic bridge.
Licence: MIT