N-body simulations of the open cluster M67 -- I: global cluster properties
Abstract
The most recent Gaia data releases have seen the construction of many open cluster catalogues containing thousands of open clusters. Simulations of these clusters are needed in order to interpret observations, both existing and upcoming. M67 is an ideal open cluster to calibrate these simulations due to its advanced age and large, diverse binary populations. Moreover it is one of the most extensively observed open clusters with observations in several wavelengths. We aim to reproduce the global cluster properties of M67 with modern N-body simulations, as well as comparing its observed properties in the largest homogeneous open cluster catalogue with other available observations. We carry out several NBODY6++GPU simulations with varying initial conditions, including metallicity, half-mass radius, and the distribution of orbital binary parameters. Our models broadly reproduce the observed global properties of M67, and we are able to recreate the cluster's surface brightness profile, radial mass profile, and colour-magnitude diagram. Our favoured model indicates an initial half-mass radius of around 1.6 pc, although its total observable cluster mass may be too high. Further simulations are planned for a subsequent study, with a focus on recreating the present-day binary populations of M67.
Related
- Identifying heating processes in simulations with an entropy-based scheme: A single jet episode in a galaxy cluster
- Constraining Effective Viscosity in the Intracluster Medium via the Thermal Sunyaev-Zeldovich Effect -- Predictions from the SLOW Constrained Coma Cluster Simulations
- Black Hole-Galaxy Correlations in Cluster Zoomed-in Simulations: GIZMO-SIMBA and TNG-Cluster
- Unresolved triples in the open cluster NGC 2437 (M 46) confirmed by the KMOS VVVX-GalCen Spectroscopic Survey
- Physics-based phenomenological modeling of binary black hole hierarchical formation 1: Synthetic universes from globular cluster simulations for GWTC
BibTeX
@misc{andersonbaldwin2026nbody,
title = {{$N$-body simulations of the open cluster M67 -- I: global cluster properties}},
author = {Jasmine S. M. Anderson-Baldwin and Jarrod R. Hurley},
howpublished = {arXiv},
year = {2026},
month = oct,
doi = {10.48550/arxiv.2610.22946},
eprint = {2610.22946},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2610.22946},
note = {CC BY 4.0},
}
RIS
TY - UNPB
TI - N-body simulations of the open cluster M67 -- I: global cluster properties
AU - Jasmine S. M. Anderson-Baldwin
AU - Jarrod R. Hurley
PY - 2026
DA - 2026/10/16
PB - arXiv
DO - 10.48550/arxiv.2610.22946
UR - https://arxiv.org/abs/2610.22946
AB - The most recent Gaia data releases have seen the construction of many open cluster catalogues containing thousands of open clusters. Simulations of these clusters are needed in order to interpret observations, both existing and upcoming. M67 is an ideal open cluster to calibrate these simulations due to its advanced age and large, diverse binary populations. Moreover it is one of the most extensively observed open clusters with observations in several wavelengths. We aim to reproduce the global cluster properties of M67 with modern N-body simulations, as well as comparing its observed properties in the largest homogeneous open cluster catalogue with other available observations. We carry out several NBODY6++GPU simulations with varying initial conditions, including metallicity, half-mass radius, and the distribution of orbital binary parameters. Our models broadly reproduce the observed global properties of M67, and we are able to recreate the cluster's surface brightness profile, radial mass profile, and colour-magnitude diagram. Our favoured model indicates an initial half-mass radius of around 1.6 pc, although its total observable cluster mass may be too high. Further simulations are planned for a subsequent study, with a focus on recreating the present-day binary populations of M67.
ER -