Cover Image

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

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 × 105 stars (with 2.0 × 104 in primordial hard binaries) performing high-resolution direct N-body (nbody6++gpu)...

Full description

Saved in:
Bibliographic Details
Main Authors: Kamlah, Albrecht (Author) , Leveque, Agostino (Author) , Spurzem, Rainer (Author) , Arca Sedda, Manuel (Author) , Askar, Aabbas (Author) , Banerjee, Sambaran (Author) , Berczik, Peter (Author) , Giersz, Mirek (Author) , Hurley, Jarrod (Author) , Belloni, Diogo (Author) , Kühmichel, Lars (Author) , Wang, Long (Author)
Format: Article (Journal)
Language:English
Published: 2022
In: Monthly notices of the Royal Astronomical Society
Year: 2022, Volume: 511, Issue: 3, Pages: 4060-4089
ISSN:1365-2966
DOI:10.1093/mnras/stab3748
Online Access:Verlag, lizenzpflichtig, Volltext: https://doi.org/10.1093/mnras/stab3748
Get full text
Author Notes:A.W.H Kamlah, A. Leveque, R. Spurzem, M. Arca Sedda, A. Askar, S. Banerjee, P. Berczik, M. Giersz, J. Hurley, D. Belloni, L. Kühmichel and L. Wang
Description
Summary: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 × 105 stars (with 2.0 × 104 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.
Item Description:Published: 27 December 2021
Gesehen am 11.04.2022
Physical Description:Online Resource
ISSN:1365-2966
DOI:10.1093/mnras/stab3748

Similar Items