Light‐element abundance variations in globular clusters

Star‐to‐star variations in abundances of the light elements carbon, nitrogen, oxygen, and sodium have been observed in stars of all evolutionary phases in all Galactic globular clusters that have been thoroughly studied. The data available for studying this phenomenon, and the hypotheses as to its origin, have both co‐evolved with observing technology; once high‐resolution spectra were available even for main‐sequence stars in globular clusters, scenarios involving multiple closely spaced stellar generations enriched by feedback from moderate‐ and high‐mass stars began to gain traction in the literature. This paper briefly reviews the observational history of globular cluster abundance inhomogeneities, discusses the presently favored models of their origin, and considers several aspects of this problem that require further study (© 2011 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Neutron star retention and millisecond pulsar production in globular clusters

We investigate the conditions by which neutron star retention in globular clusters is favoured. We find that neutron stars formed in massive binaries are far more likely to be retained. Such binaries are likely to then evolve into contact before encountering other stars, possibly producing a single neutron star after a common envelope phase. A large fraction of the single neutron stars in globular clusters are then likely to exchange into binaries containing moderate-mass main-sequence stars, replacing the lower-mass components of the original systems. These binaries will become intermediate-mass X-ray binaries (IMXBs), once the moderate-mass star evolves off the main sequence, as mass is transferred on to the neutron star, possibly spinning it up in the process. Such systems may be responsible for the population of millisecond pulsars (MSPs) that has been observed in globular clusters. Additionally, the period of mass-transfer (and thus X-ray visibility) in the vast majority of such systems will have occurred 5–10 Gyr ago, thus explaining the observed relative paucity of X-ray binaries today, given the MSP population.     

Stellar kinematics in the centers of globular clusters

The high central stellar densities in globular clusters provide a unique environment to study the fundamental dynamical process of two‐body relaxation. This process is the main driver of the dynamical evolution in the center of a globular cluster and has a profound effect on the structure of the cluster and on its stellar environment. We have obtained stellar absorption line spectra with STIS to measure the radial velocities of individual stars in the crowded center of the globular cluster M15. These data increase the number of stars with known radial velocities within the central arcsec by a factor of about three and significantly improve the constraints on the mass distribution in M15. The data provide the most detailed look of the central structure of any globular cluster and show that there is a compact dark central mass component. Similar studies using ground based facilities can be efficiently performed by employing Integral Field Units. We have started a project to better constrain the central mass density in the globular cluster M3 using the GMOS‐IFU on Gemini North. The data will also allow us to better understand the central rotation which is neither explained nor predicted by any globular cluster model. (© 2004 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

The period–luminosity relation for type II Cepheids in globular clusters

We report the result of our near-infrared observations ( JHK _s) for type II Cepheids (including possible RV Tau stars) in galactic globular clusters. We detected variations of 46 variables in 26 clusters (10 new discoveries in seven clusters) and present their light curves. Their periods range from 1.2 d to over 80 d. They show a well-defined period–luminosity relation at each wavelength. Two type II Cepheids in NGC 6441 also obey the relation if we assume the horizontal branch stars in NGC 6441 are as bright as those in metal-poor globular clusters in spite of the high metallicity of the cluster. This result supports the high luminosity which has been suggested for the RR Lyr variables in this cluster. The period–luminosity relation can be reproduced using the pulsation equation     assuming that all the stars have the same mass. Cluster RR Lyr variables were found to lie on an extrapolation of the period–luminosity relation. These results provide important constraints on the parameters of the variable stars.
Using Two Micron All-Sky Survey (2MASS) data, we show that the type II Cepheids in the Large Magellanic Cloud (LMC) fit our period–luminosity relation within the expected scatter at the shorter periods. However, at long periods (   P > 40  d, i.e. in the RV Tau star range) the LMC field variables are brighter by about one magnitude than those of similar periods in galactic globular clusters. The long-period cluster stars also differ from both these LMC stars and galactic field RV Tau stars in a colour–colour diagram. The reasons for these differences are discussed.     

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.     

球状星团成员星的空间相关性分析

以球状星团NGC (New General Catalogue) 104、NGC 5139、NGC 6121为实验样区, 选取了视差等10个恒星参数, 通过引入地学中的空间分析理论和相应的分析框架为定量描述球状星团成员星的空间分布特征提出了一种基于地学的研究范式. 通过计算全局和局部莫兰(Moran)指数得到球状星团成员星各恒星参数的空间分布特征. 研究结果表明: 球状星团NGC 104、NGC 5139、NGC 6121成员星的各恒星参数在总体上呈现出空间正相关特性, 表现出空间集聚特征, 但不同恒星参数之间存在差异; 局部空间分布也呈现聚集特征, 而不同的成员星呈现出不同的空间分布特性和趋势. 总体而言, 用地学空间相关分析系统地定量化描述球状星团成员星空间分布特征, 能够为球状星团的研究提供新的思路.     

Galactic globular clusters as a test for very-low-mass stellar models

We make use of the 'Next Generation' model atmospheres of Allard et al. and Hauschildt, Allard & Baron to compute theoretical models for low- and very-low-mass stars for selected metallicities in the range Z =0.0002 to 0.002. On this basis, we present theoretical predictions covering the sequence of H-burning stars as observed in Galactic globulars from the faint end of the main sequence up to, and beyond, the cluster turn-off. The role played by the new model atmospheres is discussed, showing that present models appear in excellent agreement with models by Baraffe et al. as computed on a quite similar physical basis. One finds that the theoretical mass–luminosity relations based on this updated set of models are in good agreement with the empirical data provided by Henry & McCarthy. Comparison with HST observation discloses that the location on the colour–magnitude diagram of the lower main sequence in Galactic globular clusters appears again in good agreement with the predicted sensitive dependence of these sequences on the cluster metallicity.     

Effects of external tidal field on the evolution of the outer regions of multi-mass star clusters

We present N -body simulations (including an initial mass function) of globular clusters in the Galaxy in order to study effects of the tidal field systematically on the properties of the outer parts of globular clusters. Using nbody6 , which correctly takes into account the two-body relaxation, we investigate the development of tidal tails of globular clusters in the Galactic tidal field. For simplicity, we have employed only the spherical components (bulge and halo) of the Galaxy, and ignored the effects of stellar evolution which could have been important in the very early phase of the cluster evolution. The total number of stars in our simulations is about 20 000, which is much smaller than the realistic number of stars. All simulations had been done for several orbital periods in order to understand the development of the tidal tails. In our scaled-down models, the relaxation time is sufficiently short to show the mass segregation effect, but we did not go far enough to see the core collapse, and the fraction of stars lost from the cluster at the end of the simulations is only ∼10 per cent. The radial distribution of extra-tidal stars can be described by a power law with a slope around −3 in surface density. The directions of tidal tails are determined by the orbits and locations of the clusters. We find that the length of tidal tails increases towards the apogalacticon and decreases towards the perigalacticon. This is an anti-correlation with the strength of the tidal field, caused by the fact that the time-scale for the stars to respond to the potential is similar to the orbital time-scale of the cluster. The escape of stars in the tidal tails towards the pericentre could be another reason for the decrease of the length of tidal tails. We find that the rotational angular velocity of tidally induced clusters shows quite different behaviour from that of initially rotating clusters.     

Globular clusters

Southern hemisphere astronomers are well placed to observe globular clusters in the Galaxy, and significant work has been done on them in South Africa. This brief review covers research on colour-magnitude diagrams and spectroscopic analyses of cluster stars as well as on miscellaneous stars representing late evolutionary stages.     

On the existence of differences in luminosity between horizontal branch stars in globular clusters and in the field

The discrepancy between a long distance-scale derived from Hipparcos -based distances to globular clusters via main-sequence fitting to local subdwarfs, and a short distance-scale derived from the absolute magnitude of field RR Lyraes via statistical parallaxes and the Baade–Wesselink method could be accounted for whether an intrinsic difference of about ∼0.1–0.2 mag was found to exist between horizontal branch (HB) stars populating the sparse general field and the dense globular clusters. In this paper we discuss the possible existence of such a systematic difference comparing the period-shifts observed for field and cluster RR Lyraes. Various approaches based on different parameters and data sets for both cluster and field variables were used in order to establish the size of such a hypothetical difference, if any. We find that on the whole very small not significant differences exist between the period–metallicity distributions of field and cluster RR Lyraes, thus confirming with a more quantitative approach, the qualitative conclusions by Catelan . This observational evidence translates into a very small difference between the horizontal branch luminosity of field and cluster stars, unless RR Lyraes in globular clusters are about 0.06 M_⊙ more massive than field RR Lyraes at same metallicity, which is to be proven.     

Expected planets in globular clusters

We argue that all transient searches for planets in globular clusters have a very low detection probability. Planets of low-metallicity stars typically do not reside at small orbital separations. The dependence of planetary system properties on metallicity is clearly seen when the quantity   I _e ≡ M _p[ a (1 − e )]²  is considered;   M _p, a   and e are the planet mass, semimajor axis and eccentricity, respectively. In high-metallicity systems, there is a concentration of systems at high and low values of I _e , with a low-populated gap near   I _e ∼ 0.3 M _J au²  , where M _J is Jupiter's mass. In low-metallicity systems, the concentration is only at the higher range of I _e , with a tail to low values of I _e . Therefore, it is still possible that planets exist around main-sequence stars in globular clusters, although at small numbers because of the low metallicity, and at orbital periods of ≳10 d. We discuss the implications of our conclusions on the role that companions can play in the evolution of their parent stars in globular clusters, for example, influencing the distribution of horizontal branch stars on the Hertzsprung–Russell diagram of some globular clusters, and in forming low-mass white dwarfs.     

Planets in 47 Tuc

We consider the survivability of planetary systems in the globular cluster 47 Tucanae. We compute the cross-sections for the breakup of planetary systems via encounters with single stars and binaries. We also compute the cross-sections to leave planets on eccentric orbits. We find that wider planetary systems ( d ≳0.3 au) are likely to be broken up in the central regions of 47 Tucanae (within the half-mass radius of the cluster). However, tighter systems and those in less-dense regions may survive. Tight systems will certainly survive in less-dense clusters where subsequent surveys should be conducted.     

Constraints on the initial conditions of globular clusters

N -body simulations are made with a variety of initial conditions, in particular clumpy and flattened distributions, to attempt to constrain the possible initial conditions of globular clusters, using the observations that young LMC globular clusters appear relaxed after only 20 to 40 Myr. It is found that violent relaxation is able to erase most of the initial substructure in only ≈ 6 crossing times. However, initially very clumpy distributions (≲ 100 clumps) form clusters that are too concentrated to resemble real globular clusters. Such clusters also often have large clumps in long-lasting (≳ 30 crossing times) orbits which do not appear in observed cluster profiles. It is also found that even modest amounts of initial flattening produce clusters that are too elliptical to resemble real globular clusters. In such a scenario, cloud–cloud collisions and similar energetic processes would be unlikely to produce sufficiently spherical globular clusters. It is suggested that globular clusters form from roughly spherical initial conditions with star formation occurring either smoothly or in many small clumps.     

The connection between globular cluster systems and the host galaxies

A large number of early-type galaxies are now known to possess blue and red subpopulations of globular clusters. We have compiled a data base of 28 such galaxies exhibiting bimodal globular cluster colour distributions. After converting to a common V – I colour system, we investigate correlations between the mean colour of the blue and red subpopulations with galaxy velocity dispersion. We support previous claims that the mean colours of the blue globular clusters are unrelated to their host galaxy. They must have formed rather independently of the galaxy potential they now inhabit. The mean blue colour is similar to that for halo globular clusters in our Galaxy and M31. The red globular clusters, on the other hand, reveal a strong correlation with galaxy velocity dispersion. Furthermore, in well-studied galaxies the red subpopulation has similar, and possibly identical, colours to the galaxy halo stars. Our results indicate an intimate link between the red globular clusters and the host galaxy; they share a common formation history. A natural explanation for these trends would be the formation of the red globular clusters during galaxy collapse.     

Galaxy evolution: Inhomogeneous halo collapse

In an effort to study the halo globular clusters we perform a simulation of inhomogeneous halo collapse using a SPH/N-body code. We include dark matter particles, (cooling) gas, and stars. Supernova feedback is included as a heat source for the gas in the local environment of star forming regions. We find this model cannot reproduce the age distribution and spatial distribution of halo globular clusters. This revised version was published online in August 2006 with corrections to the Cover Date.     

Self-enrichment by asymptotic giant branch stars in globular clusters: comparison between intermediate and high metallicities

We present theoretical evolutionary sequences of intermediate-mass stars  ( M = 3 − 6.5 M_⊙)  with metallicity   Z = 0.004  . Our goal is to test whether the self-enrichment scenario by massive asymptotic giant branch stars may work for the high-metallicity globular clusters, after previous works by the same group showed that the theoretical yields by this class of objects can reproduce the observed trends among the abundances of some elements, namely the O–Al and O–Na anticorrelations, at intermediate metallicities, i.e.  [Fe/H]=−1.3  . We find that the increase in the metallicity favours only a modest decrease of the luminosity and the temperature at the bottom of the envelope for the same core mass, and also the efficiency of the third dredge-up is scarcely altered. On the contrary, differences are found in the yields, due to the different impact that processes with the same efficiency have on the overall abundance of envelopes with different metallicities. We expect the same qualitative patterns as in the intermediate-metallicity case, but the slopes of some of the relationships among the abundances of some elements are different. We compare the sodium–oxygen anticorrelation for clusters of intermediate metallicity ( Z ≈ 10⁻³) and clusters of metallicity large as in these new models. Although the observational data are still too scarce, the models are consistent with the observed trends, provided that only stars of   M ≳ 5 M_⊙  contribute to self-enrichment.     

The environment of formation as a second parameter for globular cluster classification

We perform an evolutionary multivariate analysis of a sample of 54 Galactic globular clusters with high-quality colour–magnitude diagrams and well-determined ages. The four parameters adopted for the analysis are: metallicity, age, maximum temperature on the horizontal branch and absolute V magnitude. Our cladistic analysis breaks the sample into three novel groups. An a posteriori kinematical analysis puts groups 1 and 2 in the halo, and group 3 in the thick disc. The halo and disc clusters separately follow a luminosity–metallicity relation of much weaker slope than galaxies. This property is used to propose a new criterion for distinguishing halo and disc clusters. A comparison of the distinct properties of the two halo groups with those of Galactic halo field stars indicates that the clusters of group 1 originated in the inner halo, while those of group 2 formed in the outer halo of the Galaxy. The inner halo clusters were presumably initially the most massive one, which allowed the formation of more strongly helium-enriched second generation stars, thus explaining the presence of Cepheids and of very hot horizontal-branch stars exclusively in this group. We thus conclude that the 'second parameter' is linked to the environment in which globular clusters form, the inner halo favouring the formation of the most massive clusters which subsequently become more strongly self-enriched than their counterparts of the galactic outer halo and disc.     

Radial velocities of five globular clusters obtained with AAOmega

Using the recently commissioned multi‐object spectrograph AAOmega on the 3.9m AAT we have obtained mediumresolution near‐infrared spectra for 10 500 stars in and around five southern globular clusters. The targets were 47 Tuc, M12, M30, M55 and NGC 288. We have measured radial velocities to ± 1 kms‐ 1 with the cross correlation method and estimated metallicity, effective temperature, surface gravity and rotational velocity for each star by fitting synthetic model spectra. An analysis of the velocity maps and velocity dispersion of member stars revealed systemic rotation in four of the target clusters. (© 2007 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Pulsational MV versus [Fe/H] relation(s) for globular cluster RR Lyrae variables

We use the results from recent computations of updated non-linear convective pulsating models to constrain the distance modulus of Galactic globular clusters through the observed periods of first-overtone (RR _c ) pulsators. The resulting relation between the mean absolute magnitude of RR Lyrae stars 〈 M _V (RR)〉 and the heavy element content [Fe/H] appears well in the range of several previous empirical calibrations, but with a non-linear dependence on [Fe/H] so that the slope of the relation increases when moving towards larger metallicities. On this ground, our results suggest that metal-poor ([Fe/H]<−1.5) and metal-rich ([Fe/H]>−1.5) variables follow two different linear 〈 M _V (RR)〉−[Fe/H] relations. Application to RR Lyrae stars in the metal-poor globular clusters of the Large Magellanic Cloud (LMC) provides an LMC distance modulus of the order of 18.6 mag, thus supporting the 'long' distance scale. The comparison with recent predictions based on updated stellar evolution theory is briefly presented and discussed.     

Monte Carlo simulations of star clusters – VI. The globular cluster NGC 6397

We describe Monte Carlo models for the dynamical evolution of the nearby globular cluster NGC 6397. The code includes treatments of two-body relaxation, most kinds of three- and four-body interactions involving primordial binaries and those formed dynamically, the Galactic tide and the internal evolution of both single and binary stars. We arrive at a set of initial parameters for the cluster which, after 12 Gyr of evolution, gives a model with a fairly satisfactory match to the surface brightness profile, the velocity dispersion profile and the luminosity function in two fields. We describe in particular those aspects of the evolution which distinguish this cluster from M4, which has a roughly similar mass and Galactocentric distance, but a qualitatively different surface brightness profile. Within the limitations of our modelling, we conclude that the most plausible explanation for the difference is fluctuations: both clusters are post-collapse objects, but sometimes have resolvable cores and sometimes not.