A unified N-body and statistical treatment of stellar dynamics. III - Early postcollapse evolution of globular clusters

The period immediately following the core collapse phase in the evolution of a globular cluster is studied using a hybrid N-body/Fokker-Planck stellar dynamical code. Several core oscillations of the type predicted in earlier work are seen. The oscillations are driven by the formation, hardening, and ejection of binaries by three-body processes, and appear to decay on a timescale of about 10 to the 7th yr, for the choice of 'typical' cluster parameters made here. There is no evidence that they are gravothermal in nature. The mechanisms responsible for the decay are discussed in some detail. The distribution of hard binaries produced by the oscillations is compared with theoretical expectations and the longer term evolution of the system is considered.