Monte Carlo simulations of star clusters – III. A million-body star cluster

A revision of Stodółkiewicz's Monte Carlo code is used to simulate the evolution of million-body star clusters. The new method treats each superstar as a single star and follows the evolution and motion of all individual stellar objects. The evolution of N -body systems influenced by the tidal field of a parent galaxy and by stellar evolution is presented. All models consist of 1 000 000 stars. The process of energy generation is realized by means of appropriately modified versions of Spitzer's and Mikkola's formulae for the interaction cross-section between binaries and field stars and binaries themselves. The results presented are in good agreement with theoretical expectations and the results of other methods. During the evolution, the initial mass function (IMF) changes significantly. The local mass function around the half-mass radius closely resembles the actual global mass function. At the late stages of evolution, the mass of the evolved stars inside the core can be as high as 97 per cent of the total mass in this region. For the whole system, the evolved stars can compose up to 75 per cent of the total mass. The evolution of cluster anisotropy strongly depends on initial cluster concentration, IMF and the strength of the tidal field. The results presented are the first step in the direction of simulating the evolution of real globular clusters by means of the Monte Carlo method.     

Photometric study of the variable star population in the globular cluster NGC 6397

We present the results of a photometric survey for variable stars in the central region of the nearby globular cluster NGC 6397. Time-series photometry was obtained for 30 variable objects. The sample includes 12 new objects, of which six show periodic light curves and two are eclipsing binaries of unknown period. Six variables possess certain and three possess likely X-ray counterparts detected with the Chandra observatory. Among them, four are cataclysmic variables and one is a foreground eclipsing binary. The cataclysmic variable CV2 exhibited a likely dwarf nova type outburst in 2003 May. The cataclysmic variable CV3 was observed at  18.5 < V < 20.0  during five observing runs, but went into a low state in 2003 May when it reached   V > 22  . We have found that the light curve of the optical companion to the millisecond pulsar PSR J1740−5340 exhibits noticeable changes of its amplitude on a time-scale of a few months. A shallow eclipse with  Δ V = 0.03 mag  was detected in one of the cluster turn-off stars suggesting the presence of a large planet or brown dwarf in orbit.     

Monte Carlo simulations of star clusters – IV. Calibration of the Monte Carlo code and comparison with observations for the open cluster M67

We outline the steps needed in order to incorporate the evolution of single and binary stars into a particular Monte Carlo code for the dynamical evolution of a star cluster. We calibrate the results against N -body simulations, and present models for the evolution of the old open cluster M67 (which has been studied thoroughly in the literature with N -body techniques). The calibration is done by choosing appropriate free code parameters. We describe, in particular, the evolution of the binary, white dwarf and blue straggler populations, though not all channels for blue straggler formation are represented yet in our simulations. Calibrated Monte Carlo runs show good agreement with results of N -body simulations not only for global cluster parameters, but also for, for example, binary fraction, luminosity function and surface brightness. Comparison of Monte Carlo simulations with observational data for M67 shows that it is possible to get reasonably good agreement between them. Unfortunately, because of the large statistical fluctuations of the numerical data and uncertainties in the observational data the inferred conclusions about the cluster initial conditions are not firm.     

Monte Carlo simulations of star clusters — I. First Results

A revision of Stodoíkiewicz's Monte Carlo code is used to simulate evolution of star clusters. The new method treats each superstar as a single star and follows the evolution and motion of all individual stellar objects. The first calculations for isolated, equal-mass N -body systems with three-body energy generation according to Spitzer's formulae show good agreement with direct N -body calculations for N  = 2000, 4096 and 10 000 particles. The density, velocity, mass distributions, energy generation, number of binaries, etc., follow the N -body results. Only the number of escapers is slightly too high compared with N -body results, and there is no level-off anisotropy for advanced post-collapse evolution of Monte Carlo models as is seen in N -body simulations for N  ≤ 2000. For simulations with N  > 10 000 gravothermal oscillations are clearly visible. The calculations of N   2000, 4096, 10 000, 32 000 and 100 000 models take about 2, 6, 20, 130 and 2500 h, respectively. The Monte Carlo code is at least 10⁵ times faster than the N -body one for N  = 32 768 with special-purpose hardware. Thus it becomes possible to run several different models to improve statistical quality of the data and run individual models with N as large as 100 000. The Monte Carlo scheme can be regarded as a method which lies in the middle between direct N -body and Fokker–Planck models and combines most advantages of both methods.     

Monte Carlo simulations of star clusters – II. Tidally limited, multimass systems with stellar evolution

A revision of Stodółkiewicz's Monte Carlo code is used to simulate evolution of large star clusters. The new method treats each superstar as a single star and follows the evolution and motion of all individual stellar objects. A survey of the evolution of N -body systems influenced by the tidal field of a parent galaxy and by stellar evolution is presented. The process of energy generation is realized by means of appropriately modified versions of Spitzer's and Mikkola's formulae for the interaction cross-section between binaries and field stars and binaries themselves. The results presented are in good agreement with theoretical expectations and the results of other methods (Fokker–Planck, Monte Carlo and N -body). The initial rapid mass loss, resulting from stellar evolution of the most massive stars, causes expansion of the whole cluster and eventually leads to the disruption of less bound systems ( W ₀=3). Models with larger W ₀ survive this phase of evolution and then undergo core collapse and subsequent post-collapse expansion, like isolated models. The expansion phase is eventually reversed when tidal limitation becomes important. The results presented are the first major step in the direction of simulating evolution of real globular clusters by means of the Monte Carlo method.     

A stochastic Monte Carlo approach to model real star cluster evolution – II. Self-consistent models and primordial binaries

The new approach outlined in Paper I to follow the individual formation and evolution of binaries in an evolving, equal point-mass star cluster is extended for the self-consistent treatment of relaxation and close three- and four-body encounters for many binaries (typically a few per cent of the initial number of stars in the cluster mass). The distribution of single stars is treated as a conducting gas sphere with a standard anisotropic gaseous model. A Monte Carlo technique is used to model the motion of binaries, their formation and subsequent hardening by close encounters, and their relaxation (dynamical friction) with single stars and other binaries. The results are a further approach towards a realistic model of globular clusters with primordial binaries without using special hardware. We present, as our main result, the self-consistent evolution of a cluster consisting of 300 000 equal point-mass stars, plus 30 000 equal-mass binaries over several hundred half-mass relaxation times, well into the phase where most of the binaries have been dissolved and evacuated from the core. The cluster evolution is about three times slower than found by Gao et al. Other features are rather comparable. At every moment we are able to show the individual distribution of binaries in the cluster.     

A two-dimensional Kolmogorov–Smirnov test for crowded field source detection: ROSAT sources in NGC 6397

We present a two-dimensional version of the classical one-dimensional Kolmogorov–Smirnov (KS) test, extending an earlier idea due to Peacock and an implementation proposed by Fasano and Franceschini. The two-dimensional KS test is used to optimize the goodness of fit in an iterative source-detection scheme for astronomical images. The method is applied to a ROSAT /HRI X-ray image of the post-core-collapse globular cluster NGC 6397 to determine the most probable source distribution in the cluster core. Comparisons to other widely used source-detection methods, and to a Chandra image of the same field, show that our iteration scheme is superior in measuring statistics-limited sources in severely crowded fields.     

Radio observations of NGC 2808 and other globular clusters: constraints on intermediate-mass black holes

We present the results of a deep radio observation of the globular cluster NGC 2808. We show that there are no sources detected within the core of the cluster, placing constraints on both the pulsar population of the cluster and the mass of a possible intermediate-mass black hole in NGC 2808. We compare the results for this cluster with other constraints on intermediate-mass black holes derived from accretion measures. With the exception of G1 in M 31 which has previously shown radio emission, even with considerably more conservative assumptions, only the clusters with the poorest of observational constraints are consistent with falling on the   M _BH–σ  relation. This result is interpreted in terms of the fundamental differences between galaxies and globular clusters.     

The outer-halo globular cluster NGC 6229 – IV. Variable stars

We present time-series B , V photometry of NGC 6229, obtained with the purpose of providing the first extensive CCD variability study of this cluster. As a result, we were able to obtain periods for all NGC 6229 variables, with the exception of five stars from the candidate list of Borissova et al. located very close to the cluster centre. Two stars suspected to be variables by Carney, Fullton and Trammell are first-overtone RR Lyrae (RRc) stars, whereas seven of the 12 candidates of Borissova et al. are confirmed variables – three of them being fundamental RR Lyrae (RRab) pulsators, two first-overtone pulsators, one an eclipsing binary (most likely an Algol system) and one a bright star whose variability status could not be satisfactorily determined. A new image subtraction method (ISM) suggested by Alard has been employed which, together with the Welch–Stetson technique, has allowed us to discover twelve new RR Lyrae variables in the cluster, for which ephemerides are provided. Ten of these are RRabs, whereas the other two are RRcs. As originally suggested by Mayer, NGC 6229 is clearly an Oosterhoff type I globular cluster. We also confirm that V8 is a Population II Cepheid of the W Virginis type, and suspect that this is the case for V22 as well. The physical properties of the NGC 6229 RR Lyrae population are contrasted with those of M3 (NGC 5272) using several different methods, including a standard period-shift analysis. Possible differences between these two clusters are discussed.     

The globular cluster system of the young elliptical NGC 6702

We study the globular cluster (GC) system of the dust-lane elliptical galaxy NGC 6702, using B -, V - and I -band imaging observations carried out at the Keck telescope. This galaxy has a spectroscopic age of ≈2 Gyr suggesting recent star formation. We find strong evidence for a bimodal GC colour distribution, with the blue peak having a colour similar to that of the Galactic halo GCs. Assuming that the blue GCs are indeed old and metal-poor, we estimate an age of 2–5 Gyr and supersolar metallicity for the red GC subpopulation. Despite the large uncertainties, this is in reasonable agreement with the spectroscopic galaxy age. Additionally, we estimate a specific frequency of S _N =2.3±1.1 for NGC 6702. We predict that passive evolution of NGC 6702 will further increase its specific frequency to S _N ≈2.7 within 10 Gyr, in closer agreement to that of typical present-day ellipticals. We also discuss evidence that the merger/accretion event that took place a few Gyr ago involved a high gas fraction.     

Basic N-body modelling of the evolution of globular clusters – I. Time scaling

We consider the use of N -body simulations for studying the evolution of rich star clusters (i.e. globular clusters).The dynamical processes included in this study are restricted to gravitational (point-mass) interactions, the steady tidal field of a galaxy, and instantaneous mass loss resulting from stellar evolution. With evolution driven by these mechanisms, it is known that clusters fall roughly into two broad classes: those that dissipate promptly in the tidal field, as a result of mass loss; and those that survive long enough for their evolution to become dominated by two-body relaxation.
The time-scales of the processes we consider scale in different ways with the number of stars in the simulation, and the main aim of the paper is to suggest how the scaling of a simulation should be done so that the results are representative of the evolution of a 'real' cluster. We investigate three different ways of scaling time. One of these is appropriate to the first type of cluster, i.e. those that dissipate rapidly; similarly, a second scaling is appropriate only to the second (relaxation-dominated) type. We also develop a hybrid scaling, which is a satisfactory compromise for both types of cluster. Finally we present evidence that the widely used Fokker–Planck method produces models that are in good agreement with N -body models of those clusters that are relaxation-dominated, at least for N -body models with several thousand particles, but that the Fokker–Planck models evolve too fast for clusters that dissipate promptly.     

The evolution of binary fractions in globular clusters

We study the evolution of binary stars in globular clusters using a new Monte Carlo approach combining a population synthesis code ( startrack ) and a simple treatment of dynamical interactions in the dense cluster core using a new tool for computing three- and four-body interactions ( fewbody ). We find that the combination of stellar evolution and dynamical interactions (binary–single and binary–binary) leads to a rapid depletion of the binary population in the cluster core. The maximum binary fraction today in the core of a typical dense cluster such as 47 Tuc, assuming an initial binary fraction of 100 per cent, is only ∼ 5–10 per cent. We show that this is in good agreement with recent Hubble Space Telescope observations of close binaries in the core of 47 Tuc, provided that a realistic distribution of binary periods is used to interpret the results. Our findings also have important consequences for the dynamical modelling of globular clusters, suggesting that 'realistic models' should incorporate much larger initial binary fractions than has usually been the case in the past.