Binary Population Synthesis with BSE

Read the paper for details: Millisecond Pulsars from Accretion Induced Collapse naturally explain the Galactic Center Gamma-ray Excess

Supplementary reading:

• Nature Astronomy publication
• Nature Astronomy blog post

Getting started

Using from scratch:

• Install Python 3
• Visit the AMUSE docs and install all the dependancies
• Install AMUSE Framework and BSE with pip
  

  pip install amuse-framework
  pip install amuse-bse
  

• Clone this repository

  git clone https://github.com/gautam-404/Binary-Evolution.git
  

Running the code

• Run

  cd Binary-Evolution
  python main.py
  

  In case the python program throws an MPI exception (like MPI_ERR_SPAWN: could not spawn processes), explicitly use Open-MPI instead of MPICH (refer this doc) and run echo 'export OMPI_MCA_rmaps_base_oversubscribe=yes' >> ~/.zshrc. For Bash shell, replace ~/.zshrc with ~/.bashrc. Restart your terminal.

The python code will create an output folder with evolution histories of all binary systems at the path "~/OutputFiles". See /src/BPS_eval.py which contains necessary functions that produce the data to see the total gamma-ray luminosity produced by the population as a function of time.

Contact for more info.