Towards an N-body Model for the Globular Cluster M4

Abstract

This paper describes an N-body model for the dynamical evolution of the nearby globular cluster M4. The initial conditions, with N = 484710 particles, were generated from a published study of this cluster with a Monte Carlo code. With the Monte Carlo code, these initial conditions led, after 12 Gyr of dynamical and stellar evolution, to a model which resembles M4 in terms of its surface brightness and velocity dispersion profiles, and its local luminosity function. Though the N-body model reported here is marred by some errors, its evolution can be compared with that of the published Monte Carlo model, with a result from the synthetic evolution code EMACSS, and with M4 itself.