Observations and theoretical work suggest that globular clusters may be born with initially very large binary fractions. We present first results from our newly modified Monte-Carlo cluster evolution code, which treats binary interactions exactly via direct N-body integration. It is shown that binary scattering interactions generate significantly less energy than predicted by the recipes that have been used in the past to model them in approximate cluster evolution methods. The new result that the cores of globular clusters in the long-lived binary-burning phase are smaller than previously predicted weakens the agreement with observations, thus implying that more than simply stellar dynamics is at work in shaping the globular clusters we observe today.

We study the dynamical evolution of globular clusters using our 2D Monte Carlo code with the inclusion of primordial binary interactions for equal-mass stars. We use approximate analytical cross sections for energy generation from binary-binary and binary-single interactions. After a brief period of slight contraction or expansion of the core over the first few relaxation times, all clusters enter a much longer phase of stable "binary burning" lasting many tens of relaxation times. The structural parameters of our models during this phase match well those of most observed globular clusters. At the end of this phase, clusters that have survived tidal disruption undergo deep core collapse, followed by gravothermal oscillations. Our results clearly show that the presence of even a small fraction of binaries in a cluster is sufficient to support the core against collapse significantly beyond the normal core collapse time predicted without the presence of binaries. For tidally truncated systems, collapse is easily delayed sufficiently that the cluster will undergo complete tidal disruption before core collapse. As a first step toward the eventual goal of computing all interactions exactly using dynamical three- and four-body integration, we have incorporated an exact treatment of binary-single interactions in our code. We show that results using analytical cross sections are in good agreement with those using exact three-body integration, even for small binary fractions where binary-single interactions are energetically most important.

We report the results of analysis of observations of the Vela Pulsar by PCA on RXTE. Our data consists of two parts. The first part contains observations at 1, 4, and 9 days after the glitch in 1996 and has 27000 sec. total exposure time. The second part of observations were performed three months after this glitch and have a total exposure time of 93000 sec. We found pulsations in both sets. The observed spectrum is a power-law with no apparent change in flux or count rate. The theoretical expectations of increase in flux due to internal heating after a glitch are smaller than the uncertainty of the observations.