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.     

Mass transfer from a giant star to a main-sequence companion and its contribution to long-orbital-period blue stragglers

Binary population synthesis shows that mass transfer from a giant star to a main-sequence (MS) companion may account for some observed long-orbital-period blue stragglers. However, little attention is paid to this blue straggler formation scenario as dynamical instability often happens when the mass donor is a giant star. In this paper, we have studied the critical mass ratio, q _c, for dynamically stable mass transfer from a giant star to a MS companion using detailed evolution calculations. The results show that a more evolved star is generally less stable for Roche lobe overflow. Meanwhile,   q _c  almost linearly increases with the amount of the mass and angular momentum lost during mass transfer, but has little dependance on stellar wind. To conveniently use the result, we give a fit of q _c as a function of the stellar radius at the onset of Roche lobe overflow and of the mass-transfer efficiency during the Roche lobe overflow.
To examine the formation of blue stragglers from mass transfer between giants and MS stars, we have performed Monte Carlo simulations with various q _c. The simulations show that some binaries with the mass donor on the first giant branch may contribute to blue stragglers with q _c obtained in this paper but will not from previous q _c. Meanwhile, from our q _c, blue stragglers from the mass transfer between an asymptotic giant branch star and a MS companion may be more numerous and have a wider range of orbital periods than those from the other q _c.     

Ratios of star cluster core and half-mass radii: a cautionary note on intermediate-mass black holes in star clusters

There is currently much interest in the possible presence of intermediate-mass black holes (IMBHs) in the cores of globular clusters (GCs). Based on theoretical arguments and simulation results it has previously been suggested that a large core radius – or particularly a large ratio of the core radius to half-mass radius – is a promising indicator for finding such a black hole (BH) in a star cluster. In this study N -body models of 100 000 stars with and without primordial binaries are used to investigate the long-term structural evolution of star clusters. Importantly, the simulation data are analysed using the same processes by which structural parameters are extracted from observed star clusters. This gives a ratio of the core and half-mass (or half-light) radii that are directly comparable to the Galactic GC sample. As a result, it is shown that the ratios observed for the bulk of this sample can be explained without the need for an IMBH. Furthermore, it is possible that clusters with large core to half-light radius ratios harbour a BH binary (comprising stellar mass BHs) rather than a single massive BH. This work does not rule out the existence of IMBHs in the cores of at least some star clusters.     

The effect of gas loss on the formation of bound stellar clusters

The effect of gas ejection on the structure and binding energy of newly formed stellar clusters is investigated. The star formation efficiency (SFE), necessary for forming a gravitationally bound stellar cluster, is determined.
Two sets of numerical N -body simulations are presented. As a first simplified approach we treat the residual gas as an external potential. The gas expulsion is approximated by reducing the gas mass to zero on a given time-scale, which is treated as a free parameter. In a second set of simulations we use smoothed particle hydrodynamics (SPH) to follow the dynamics of the outflowing residual gas self-consistently. We investigate cases where gas outflow is induced by an outwards propagating shock front and where the whole gas cloud is heated homogeneously, leading to ejection.
If the stars are in virial equilibrium with the gaseous environment initially, bound clusters only form in regions where the local SFE is larger than 50 per cent or where the gas expulsion time-scale is long compared with the dynamical time-scale. A small initial velocity dispersion of the stars leads to a compaction of the cluster during the expulsion phase and reduces the SFE needed to form bound clusters to less than 10 per cent.     

A robust statistical estimation of the basic parameters of single stellar populations – I. Method

The colour–magnitude diagrams of resolved single stellar populations, such as open and globular clusters, have provided the best natural laboratories to test stellar evolution theory. Whilst a variety of techniques have been used to infer the basic properties of these simple populations, systematic uncertainties arise from the purely geometrical degeneracy produced by the similar shape of isochrones of different ages and metallicities. Here we present an objective and robust statistical technique which lifts this degeneracy to a great extent through the use of a key observable: the number of stars along the isochrone. Through extensive Monte Carlo simulations we show that, for instance, we can infer the four main parameters (age, metallicity, distance and reddening) in an objective way, along with robust confidence intervals and their full covariance matrix. We show that systematic uncertainties due to field contamination, unresolved binaries, initial or present-day stellar mass function are either negligible or well under control. This technique provides, for the first time, a proper way to infer with unprecedented accuracy the fundamental properties of simple stellar populations, in an easy-to-implement algorithm.     

Brown dwarf populations in open clusters

We present the results of multiple simulations of open clusters, modelling the dynamics of a population of brown dwarf members. We consider the effects of a large range of primordial binary populations, including the possibilities of having brown dwarf members contained within a binary system. We also examine the effects of various cluster diameters and masses. Our examination of a population of wide binary systems containing brown dwarfs, reveals evidence for exchange reactions whereby the brown dwarf is ejected from the system and replaced by a heavier main-sequence star. We find that there exists the possibility of hiding a large fraction of the brown dwarfs contained within the primordial binary population. We conclude that it is probable that the majority of brown dwarfs are contained within primordial binary systems which then hides a large proportion of them from detection.     

Binary formation in stellar clusters

We consider how the tidal potential of a stellar cluster or a dense molecular cloud affects the fragmentation of gravitationally unstable molecular cloud cores. We find that molecular cloud cores which would collapse to form a single star in the absence of tidal shear, can be forced to fragment if they are subjected to tides. This may enhance the frequency of binaries in star-forming regions such as Ophiuchus and the frequency of binaries with separations ≲100 au in the Orion Trapezium Cluster. We also find that clouds which collapse to form binary systems in the absence of a tidal potential will form bound binary systems if exposed to weak tidal shear. However, if the tidal shear is sufficiently strong, even though the cloud still collapses to form two fragments, the fragments are pulled apart while they are forming by the tidal shear and two single stars are formed. This sets an upper limit for the separation of binaries that form near dense molecular clouds or in stellar clusters.     

High-resolution simulations of stellar collisions between equal-mass main-sequence stars in globular clusters

We performed high-resolution simulations of two stellar collisions relevant for stars in globular clusters. We considered one head-on collision and one off-axis collision between two 0.6-M_⊙ main-sequence stars. We show that a resolution of about 100 000 particles is sufficient for most studies of the structure and evolution of blue stragglers. We demonstrate conclusively that collision products between main-sequence stars in globular clusters do not have surface convection zones larger than 0.004 M_⊙ after the collision, nor do they develop convection zones during the 'pre-main-sequence' thermal relaxation phase of their post-collision evolution. Therefore, any mechanism which requires a surface convection zone (i.e. chemical mixing or angular momentum loss via a magnetic wind) cannot operate in these stars. We show that no disc of material surrounding the collision product is produced in off-axis collisions. The lack of both a convection zone and a disc proves a continuing problem for the angular momentum evolution of blue stragglers in globular clusters.     

A search for variable stars in the globular cluster M3

We describe here the results of a photometric time-sequence survey of the globular cluster M3 (NGC 5272), in a search for contact and detached eclipsing binary stars. We have discovered only one likely eclipsing binary and one SX Phe-type star in spite of monitoring 4077 stars with V  < 20.0 and observing 25 blue straggler stars (BSS). The newly identified SX Phe star, V237, shows a light curve with a variable amplitude. Variable V238 shows variability either with a period of 0.49 d or with a period of 0.25 d. On the cluster colour–magnitude diagram, the variable occupies a position a few hundredths of a magnitude to the blue of the base of the red giant branch. V238 is a likely descendant of a binary blue straggler.   As a side result we obtained high-quality data for 42 of the previously known RR Lyrae variables, including 33 of Bailey type ab, seven type c and two double-mode pulsators. We used equations that relate the physical properties of RRc stars to their pulsation periods and Fourier parameters in order to derive masses, luminosities, temperatures and helium parameters for five of the RRc stars. We also tested equations that relate the [Fe/H], M _V and temperature of RRab stars to pulsation period and Fourier parameters. We derived [Fe/H]= −1.42 in good agreement with spectroscopic determinations.     

Multisite campaign on the open cluster M67 – I. Observations and photometric reductions

We report on an ambitious multisite campaign aimed at detecting stellar variability, particularly solar-like oscillations, in the red giant stars in the open cluster M67 (NGC 2682). During the six-week observing run, which comprised 164 telescope nights, we used nine 0.6-m to 2.1-m class telescopes located around the world to obtain uninterrupted time series photometry. We outline here the data acquisition and reduction, with emphasis on the optimization of the signal-to-noise ratio of the low-amplitude (50–500 μmag) solar-like oscillations. This includes a new and efficient method for obtaining the linearity profile of the CCD response at ultrahigh precision (∼10 parts per million). The noise in the final time series is 0.50 mmag per minute integration for the best site, while the noise in the Fourier spectrum of all sites combined is 20 μmag. In addition to the red giant stars, this data set proves to be very valuable for studying high-amplitude variable stars such as eclipsing binaries, W UMa systems and δ Scuti stars.     

The formation of a bound star cluster: from the Orion nebula cluster to the Pleiades

Direct N -body calculations are presented of the formation of Galactic clusters using GasEx , which is a variant of the code Nbody6 . The calculations focus on the possible evolution of the Orion nebula cluster (ONC) by assuming that the embedded OB stars explosively drove out 2/3 of its mass in the form of gas about 0.4 Myr ago. A bound cluster forms readily and survives for 150 Myr despite additional mass loss from the large number of massive stars, and the Galactic tidal field. This is the very first time that cluster formation is obtained under such realistic conditions. The cluster contains about 1/3 of the initial 10⁴ stars, and resembles the Pleiades cluster to a remarkable degree, implying that an ONC-like cluster may have been a precursor of the Pleiades. This scenario predicts the present expansion velocity of the ONC, which will be measurable by upcoming astrometric space missions. These missions should also detect the original Pleiades members as an associated expanding young Galactic-field subpopulation. The results arrived at here suggest that Galactic clusters form as the nuclei of expanding OB associations.
The results have wide implications, also for the formation of globular clusters and the Galactic-field and halo stellar populations. In view of this, the distribution of binary orbital periods and the mass function within and outside the model ONC and Pleiades is quantified, finding consistency with observational constraints. Advanced mass segregation is evident in one of the ONC models. The calculations show that the primordial binary population of both clusters could have been much the same as is observed in the Taurus–Auriga star-forming region. The computations also demonstrate that the binary proportion of brown dwarfs is depleted significantly for all periods, whereas massive stars attain a high binary fraction.     

Mass segregation in young stellar clusters

We investigate the evolutionary effect of dynamical mass segregation in young stellar clusters. Dynamical mass segregation acts on a time-scale of order the relaxation time of a cluster. Although some degree of mass segregation occurs earlier, the position of massive stars in rich young clusters generally reflects the cluster's initial conditions. In particular, the positions of the massive stars in the Trapezium cluster in Orion cannot be due to dynamical mass segregation, but indicate that they formed in, or near, the centre of the cluster. Implications of this for cluster formation and for the formation of high-mass stars are discussed.     

Modelling the components of binaries in the Hyades: the dependence of the mixing-length parameter on stellar mass

We present our findings based on a detailed analysis of the binaries of the Hyades, in which the masses of the components are well known. We fit the models of the components of a binary system to observations so as to give the observed total V and B − V of that system and the observed slope of the main sequence in the corresponding parts. According to our findings, there is a very definite relationship between the mixing-length parameter and the stellar mass. The fitting formula for this relationship can be given as  α= 9.19( M /M_⊙− 0.74)^0.053− 6.65  , which is valid for stellar masses greater than  0.77 M_⊙  . While no strict information is gathered for the chemical composition of the cluster, as a result of degeneracy in the colour–magnitude diagram, by adopting   Z = 0.033  and using models for the components of 70 Tau and θ² Tau we find the hydrogen abundance to be   X = 0.676  and the age to be 670 Myr. If we assume that   Z = 0.024  , then   X = 0.718  and the age is 720 Myr. Our findings concerning the mixing-length parameter are valid for both sets of the solution. For both components of the active binary system V818 Tau, the differences between radii of the models with   Z = 0.024  and the observed radii are only about 4 per cent. More generally, the effective temperatures of the models of low-mass stars in the binary systems studied are in good agreement with those determined by spectroscopic methods.