Modelling compact binary populations in globular clusters

SNIC 2019/6-61


SNIC Small Compute

Principal Investigator:

Alexey Bobrick


Lunds universitet

Start Date:


End Date:


Primary Classification:

10305: Astronomy, Astrophysics and Cosmology




I am developing code for high-performance modelling of compact binary populations in globular clusters. Globular clusters are dense groups of stars, containing millions of stars, and reaching stellar densities of about 100 000 stars per cubic light-year in the centre. Due to such high stellar densities, the single and binary stars in the centre frequently encounter each other, merge and participate in other complex interactions (e.g. exchanges, break-ups and others). We have developed a Monte Carlo code to produce populations of compact binaries (made of white dwarfs and neutron stars) in such environments. The code simulates several physical processes: stellar evolution, binary evolution, few-body interactions, globular cluster evolution, triple stellar evolution and others, and allows to model the full Milky Way globular cluster system in about 30 minutes on a 20 node machine. High efficiency of the code enables us to explore the parameter space in great detail, as well as to obtain excellent statical uncertainties for each cluster. I would like to request a small allocation on Lunarc to be able to run tests and verifications of the code efficiently, in preparation for production runs on a larger allocation we have in our group.