N-body Simulations with Live Stellar Evolution

TL;DR

This paper introduces a new N-body simulation code that integrates live stellar evolution, demonstrating its accuracy and flexibility compared to existing methods, and outlining future enhancements for modeling stars and binaries.

Contribution

The paper presents a novel combined N-body and stellar evolution code, integrating NBODY6 and STARS, enabling more accurate and flexible simulations of star clusters.

Findings

Results closely match those from NBODY6 and SSE when physics are aligned.

The new code allows easy modification of stellar physics, such as convective overshooting.

Successful simulation of 10,000 star clusters demonstrating the code's capabilities.

Abstract

An N-body code containing live stellar evolution through combination of the software packages NBODY6 and STARS is presented. Operational details of the two codes are outlined and the changes that have been made to combine them discussed. We have computed the evolution of clusters of 10 000 stars using the combined code and we compare the results with those obtained using NBODY6 and the synthetic stellar evolution code SSE. We find that, providing the physics package within STARS is set up correctly to match the parameters of the models used to construct SSE, the results are very similar. This provides a good indication that the new code is working well. We also demonstrate how this physics can be changed simply in the new code with convective overshooting as an example. Similar changes in SSE would require considerable reworking of the model fits. We conclude by outlining proposed future development of the code to include more complete models of single stars and binary star systems.