4.7 Article

Preparing the next gravitational million-body simulations: evolution of single and binary stars in nbody6++gpu, mocca, and mcluster

Journal

MONTHLY NOTICES OF THE ROYAL ASTRONOMICAL SOCIETY
Volume 511, Issue 3, Pages 4060-4089

Publisher

OXFORD UNIV PRESS
DOI: 10.1093/mnras/stab3748

Keywords

methods: numerical; software: development; software: documentation; binaries: general; stars: general; globular clusters: general

Funding

  1. Gauss Centre for Supercomputing e.V.
  2. computing and network department of NAOC
  3. DFG Priority Program 'Exploring the Diversity of Extrasolar Planets' [SP 345/20-1, SP 345/22-1]
  4. Volkswagen Foundation [90411, 97778]
  5. Polish National Science Center (NCN) [UMO-2016/23/B/ST9/02732]
  6. NRF of Ukraine 'Leading and Young Scientists Research Support' - 'Astrophysical Relativistic Galactic Objects (ARGO): life cycle of active nucleus' [2020.02/0346]
  7. National Academy of Sciences of Ukraine under the Main Astronomical Observatory GPU computing cluster project [13.2021.MM]
  8. Science Committee of the Ministry of Education and Science of the Republic of Kazakhstan [AP08856149]
  9. Ministry of Education and Science of Ukraine [M86-22.11.2021]
  10. Chinese Academy of Sciences (CAS) through the Silk Road Project at NAOC
  11. Swedish Research Council [2017-04217]
  12. Alexander von Humboldt Stiftung
  13. Deutsche Forschungsgemeinschaft (DFG
  14. German Research Foundation) [BA 4281/6-1]
  15. ESO/Gobierno de Chile
  16. SAo Paulo Research Foundation (FAPESP) [2017/14289-3]
  17. JSPS
  18. Australian Research Council Centre of Excellence for Gravitational Wave Discovery (OzGrav) [CE170100004]
  19. CAS
  20. Royal Physiographic Society of Lund
  21. Walter Gyllenberg Foundation

Ask authors/readers for more resources

This paper presents the implementation of updated stellar evolution recipes in several codes and tests them through numerical simulations of star clusters. The results show differences in density, remnant masses, and binary fractions between the different models. The mocca models produce more black holes and helium white dwarfs, while nbody6++gpu models have a larger amount of white dwarf-white dwarf binaries.
We present the implementation of updated stellar evolution recipes in the codes nbody6++gpu, mocca, and mcluster. We test them through numerical simulations of star clusters containing 1.1 x 10(5) stars (with 2.0 x 10(4) in primordial hard binaries) performing high-resolution direct N-body (nbody6++gpu) and Monte Carlo (mocca) simulations to an age of 10 Gyr. We compare models implementing either delayed or core-collapse supernovae mechanisms, a different mass ratio distribution for binaries, and white dwarf (WD) natal kicks enabled/disabled. Compared to nbody6++gpu, the mocca models appear to be denser, with a larger scatter in the remnant masses, and a lower binary fraction on average. The mocca models produce more black holes (BHs) and helium WDs, while nbody6++gpu models are characterized by a much larger amount of WD-WD binaries. The remnant kick velocity and escape speed distributions are similar for the BHs and neutron stars (NSs), and some NSs formed via electron-capture supernovae, accretion-induced collapse, or merger-induced collapse escape the cluster in all simulations. The escape speed distributions for the WDs, on the other hand, are very dissimilar. We categorize the stellar evolution recipes available in nbody6++gpu into four levels: the one implemented in previous nbody6++gpu and mocca versions (level A), state-of-the-art prescriptions (level B), some in a testing phase (level C), and those that will be added in future versions of our codes.

Authors

I am an author on this paper
Click your name to claim this paper and add it to your profile.

Reviews

Primary Rating

4.7
Not enough ratings

Secondary Ratings

Novelty
-
Significance
-
Scientific rigor
-
Rate this paper

Recommended

No Data Available
No Data Available