Binary Evolution and Neutron Stars in Globular Clusters

Abstract

We investigate the dynamical formation and evolution of binaries containing neutron stars in dense globular clusters. Our numerical simulations combine a simple Monte Carlo prescription for stellar dynamics, a sophisticated binary population synthesis code, and a small-N-body integrator for computing 3-body and 4-body interactions. Our results suggest that there is no ``retention problem,'' i.e., that, under standard assumptions, globular clusters can retain enough neutron stars to produce the observed numbers of millisecond pulsars. We also identify the dominant evolutionary and dynamical channels through which globular clusters produce their two main types of binary millisecond pulsars