Cluster Half-mass Radius Research Articles

Determining the dynamical age of the LMC globular cluster NGC 1835 using the ‘dynamical clock’

In the context of the study of the size–age relationship observed in star clusters in the Large Magellanic Cloud (LMC) and the investigation of its origin, we present the determination of the structural parameters and the dynamical age of the massive cluster NGC 1835. We used the powerful combination of optical and near-ultraviolet images acquired with the WFC3 on board the HST to construct the star density profile from resolved star counts, determining the values of the core, half-mass, and tidal radii through comparison with the King model family. The same data also allowed us to evaluate the dynamical age of the cluster by using the ‘dynamical clock’. This is an empirical method that quantifies the level of the central segregation of blue stragglers stars (BSSs) within the cluster half-mass radius by means of the Arh+ parameter, which is defined as the area enclosed between the cumulative radial distribution of BSSs and that of a reference (lighter) population. The results confirm that NGC 1835 is a very compact cluster with a core radius of only 0.84 pc. The estimated value of Arh+ (0.30 ± 0.04) is the largest measured so far in the LMC clusters, providing evidence of a highly dynamically evolved stellar system. NGC 1835 fits nicely into the correlation between Arh+ and the central relaxation time and in the anti-correlation between Arh+ and the core radius defined by the Galactic and Magellanic Cloud clusters investigated to date.

Empirical Measurement of the Dynamical Ages of Three Globular Clusters and Some Considerations on the Use of the Dynamical Clock * *Based on observations collected at the Hubble Space Telescope, under proposals GO12517 (PI: Ferraro), GO13410 (PI: Pallanca), GO15232 (PI: Ferraro).

We have used the “dynamical clock” to measure the level of dynamical evolution reached by three Galactic globular clusters (namely, NGC 3201, NGC 6316, and NGC 6440). This is an empirical method that quantifies the level of central segregation of blue straggler stars (BSSs) within the cluster half-mass radius by means of the parameter, defined as the area enclosed between the cumulative radial distribution of BSSs and that of a lighter population. The total sample with homogeneous determinations of currently includes 59 clusters: 52 old GCs in the Milky Way (including the three investigated here), five old clusters in the Large Magellanic Cloud, and two young systems in the Small Magellanic Cloud. The three objects studied here nicely nest into the correlation between and the central relaxation time defined by the previous sample, thus proving and consolidating the use of the dynamical clock as an excellent tracer of the stage of dynamical evolution of a star cluster in different galactic environments. Finally, we discuss the advantages of using the dynamical clock as an indicator of the dynamical ages of star clusters, compared to the present-day central relaxation time.

The growth of H ii regions around massive stars: the role of metallicity and dust

ABSTRACT Gas metallicity (Z) and the related dust-to-gas ratio (fd) can influence the growth of H ii regions via metal line cooling and ultraviolet (UV) absorption. We model these effects in star-forming regions containing massive stars. We compute stellar feedback from photoionization and radiation pressure (RP) using Monte Carlo radiative transfer coupled with hydrodynamics, including stellar and diffuse radiation fields. We follow a $10^{5}\, \mathrm{M}_{\odot }$ turbulent cloud with Z/Z⊙ = 2, 1, 0.5, and 0.1, and $f_\textrm{d} = 0.01\, Z/Z_{\odot }$ with a cluster-sink particle method for star formation. The models evolve for at least 1.5 Myr under feedback. Lower Z results in higher temperatures and therefore larger H ii regions. For Z ≥ Z⊙, RP (Prad) can dominate locally over the gas pressure (Pgas) in the inner half-parsec around sink particles. Globally, the ratio of Prad/Pgas is around 1 (2 Z⊙), 0.3 (Z⊙), 0.1 (0.5 Z⊙), and 0.03 (0.1 Z⊙). In the solar model, excluding RP results in an ionized volume several times smaller than the fiducial model with both mechanisms. Excluding RP and UV attenuation by dust results in a larger ionized volume than the fiducial case. That is, UV absorption hinders growth more than RP helps it. The radial expansion velocity of ionized gas reaches +15 km s−1 outwards, while neutral gas has inward velocities for most of the runtime, except for 0.1 Z⊙ that exceeds +4 km s−1. Z and fd do not significantly alter the star formation efficiency, rate, or cluster half-mass radius, with the exception of 0.1 Z⊙ due to the earlier expulsion of neutral gas.

Hunting for Runaways from the Orion Nebula Cluster

Abstract We use Gaia DR2 to hunt for runaway stars from the Orion Nebula Cluster (ONC). We search a region extending 45° around the ONC and out to 1 kpc to find sources that have overlapped in angular position with the cluster in the last ∼10 Myr. We find ∼17,000 runaway/walkaway candidates that satisfy this 2D traceback condition. Most of these are expected to be contaminants, e.g., caused by Galactic streaming motions of stars at different distances. We thus examine six further tests to help identify real runaways, namely: (1) possessing young stellar object (YSO) colors and magnitudes based on Gaia optical photometry; (2) having IR excess consistent with YSOs based on 2MASS and Wide-field Infrared Survey Explorer photometry; (3) having a high degree of optical variability; (4) having closest approach distances well-constrained to within the cluster half-mass radius; (5) having ejection directions that avoid the main Galactic streaming contamination zone; and (6) having a required radial velocity (RV) for 3D overlap of reasonable magnitude (or, for the 7% of candidates with measured RVs, satisfying 3D traceback). Thirteen sources, not previously noted as Orion members, pass all these tests, while another twelve are similarly promising, except they are in the main Galactic streaming contamination zone. Among these 25 ejection candidates, ten with measured RVs pass the most restrictive 3D traceback condition. We present full lists of runaway/walkaway candidates, estimate the high-velocity population ejected from the ONC, and discuss its implications for cluster formation theories via comparison with numerical simulations.

The star cluster survivability after gas expulsion is independent of the impact of the Galactic tidal field

We study the impact of the tidal field on the survivability of star clusters following instantaneous gas expulsion. Our model clusters are formed with a centrally-peaked star-formation efficiency profile as a result of star-formation taking place with a constant efficiency per free-fall time. We define the impact of the tidal field as the ratio of the cluster half-mass radius to its Jacobi radius immediately after gas expulsion, $\lambda = r_{h}/R_{J}$. We vary $\lambda$ by varying either the Galactocentric distance, or the size (hence volume density) of star clusters. We propose a new method to measure the violent relaxation duration, in which we compare the total mass-loss rate of star clusters with their stellar evolutionary mass-loss rate. That way, we can robustly estimate the bound mass fraction of our model clusters at the end of violent relaxation. The duration of violent relaxation correlates linearly with the Jacobi radius, when considering identical clusters at different Galactocentric distances. In contrast, it is nearly constant for the solar neighbourhood clusters, slightly decreasing with $\lambda$. The violent relaxation does not last longer than 50 Myr in our simulations. Identical model clusters placed at different Galactocentric distances have the same final bound fraction, despite experiencing different impacts of the tidal field. The solar neighbourhood clusters with different densities experience only limited variations of their final bound fraction. In general, we conclude that the cluster survivability after instantaneous gas expulsion, as measured by their bound mass fraction at the end of violent relaxation, $F_{bound}$, is independent of the impact of the tidal field, $\lambda$.

Neutron stars and millisecond pulsars in star clusters: implications for the diffuse γ-radiation from the Galactic Centre

Globular clusters (GCs) are the ideal environment for the formation of neutron stars (NSs) and millisecond pulsars (MSPs). NSs origin and evolution provide a useful information on stellar dynamics and evolution in star clusters, and are among the most interesting astrophysical objects, being precursors of several high-energy phenomena such as gravitational waves and gamma-ray bursts. Due to a large velocity kick that they receive at birth, most of the NSs escape the local field, affecting the evolution and dynamics of their parent cluster. In this paper, we study the origin and dynamical evolution of NSs within GCs with different initial masses, metallicities and primordial binary fractions. We find that the radial profile of NSs is shaped by the BH content of the cluster, which partially quenches the NS segregation until most of the BHs are ejected from the system. Independently on the cluster mass and initial configuration, the NSs map the average stellar population, as their average radial distance is $\approx 60-80\%$ of the cluster half-mass radius. Finally, by assuming a recycling fraction of $f_\mathrm{rec}=0.1$ and an average MSP gamma-ray emission of $L_\gamma=2\times 10^{33}$ erg s$^{-1}$, we show that the typical gamma-ray emission from our GCs agrees with observations and supports the MSP origin of the gamma-ray excess signal observed by the Fermi-LAT telescope in the Galactic Centre.

The signatures of the parental cluster on field planetary systems

Due to the high stellar densities in young clusters, planetary systems formed in these environments are likely to have experienced perturbations from encounters with other stars. We carry out direct $N$-body simulations of multi-planet systems in star clusters to study the combined effects of stellar encounters and internal planetary dynamics. These planetary systems eventually become part of the Galactic field population the parental cluster dissolves, which is where most presently-known exoplanets are observed. We show that perturbations induced by stellar encounters lead to distinct signatures in the field planetary systems, most prominently, the excited orbital inclinations and eccentricities. Planetary systems that form within the cluster's half-mass radius are more prone to such perturbations. The orbital elements are most strongly excited in the outermost orbit, but the effect propagates to the entire planetary system through secular evolution. Planet ejections may occur long after a stellar encounter. The surviving planets in these reduced systems tend to have, on average, higher inclinations and larger eccentricities compared to systems that were perturbed less strongly. As soon as the parental star cluster dissolves, external perturbations stop affecting the escaped planetary systems, and further evolution proceeds on a relaxation time scale. The outer regions of these ejected planetary systems tend to relax so slowly that their state carries the memory of their last strong encounter in the star cluster. Regardless of the stellar density, we observe a robust anticorrelation between multiplicity and mean inclination/eccentricity. We speculate that the "Kepler dichotomy" observed in field planetary systems is a natural consequence of their early evolution in the parental cluster.

A Chandra X-Ray Census of the Interacting Binaries in Old Open Clusters—Collinder 261

Abstract We present the first X-ray study of Collinder 261 (Cr 261), which at an age of 7 Gyr is one of the oldest open clusters known in the Galaxy. Our observation with the Chandra X-Ray Observatory is aimed at uncovering the close interacting binaries in Cr 261, and reaches a limiting X-ray luminosity of (0.3–7 keV) for stars in the cluster. We detect 107 sources within the cluster half-mass radius r h , and we estimate that among the sources with , ∼26 are associated with the cluster. We identify a mix of active binaries and candidate active binaries, candidate cataclysmic variables, and stars that have “straggled” from the main locus of Cr 261 in the color–magnitude diagram. Based on a deep optical source catalog of the field, we estimate that Cr 261 has an approximate mass of 6500 M ⊙, roughly the same as the old open cluster NGC 6791. The X-ray emissivity of Cr 261 is similar to that of other old open clusters, supporting the trend that they are more luminous in X-rays per unit mass than old populations of higher (globular clusters) and lower (the local neighborhood) stellar density. This implies that the dynamical destruction of binaries in the densest environments is not solely responsible for the observed differences in X-ray emissivity.

On the mass of the Galactic star cluster NGC 4337

Only a small number of galactic open clusters survives for longer than few hundred million years. Longer lifetimes are routinely explained in term of larger initial masses, particularly quiet orbits, and off-plane birth-places. We derive in this work the actual mass of NGC 4337, one of the few open clusters in the Milky Way inner disk that managed to survive for about 1.5 Gyr. We derive its mass in two different ways. First, we exploit an unpublished photometric data set in the UBVI passbands to estimate -using star counts- the cluster luminosity profile, and luminosity and mass function, and hence its actual mass both from the luminosity profile and from the mass function.This data-set is also used to infer crucial cluster parameters, as the cluster half-mass radius and distance. Second, we make use of a large survey of cluster star radial velocities to derive dynamical estimates for the cluster mass. Under the assumption of virial equilibrium and neglecting the external gravitational field leads to values for the mass significantly larger than those obtained by mean of observed density distribution or with the mass function but still marginally compatible with the inferred values of the invisible mass in form of both low mass stars or remnants of high mass stars in the cluster. Finally, we derive the cluster initial mass by computing the mass loss experienced by the cluster during its lifetime, and adopting the various estimates of the actual mass.

On the link between energy equipartition and radial variation in the stellar mass function of star clusters

We make use of $N$-body simulations to determine the relationship between two observable parameters that are used to quantify mass segregation and energy equipartition in star clusters. Mass segregation can be quantified by measuring how the slope of a cluster's stellar mass function $\alpha$ changes with clustercentric distance r, and then calculating $\delta_\alpha = \frac{d \alpha(r)}{d ln(r/r_m)}$ where $r_m$ is the cluster's half-mass radius. The degree of energy equipartition in a cluster is quantified by $\eta$, which is a measure of how stellar velocity dispersion $\sigma$ depends on stellar mass m via $\sigma(m) \propto m^{-\eta}$. Through a suite of $N$-body star cluster simulations with a range of initial sizes, binary fractions, orbits, black hole retention fractions, and initial mass functions, we present the co-evolution of $\delta_\alpha$ and $\eta$. We find that measurements of the global $\eta$ are strongly affected by the radial dependence of $\sigma$ and mean stellar mass and the relationship between $\eta$ and $\delta_\alpha$ depends mainly on the cluster's initial conditions and the tidal field. Within $r_m$, where these effects are minimized, we find that $\eta$ and $\delta_\alpha$ initially share a linear relationship. However, once the degree of mass segregation increases such that the radial dependence of $\sigma$ and mean stellar mass become a factor within $r_m$, or the cluster undergoes core collapse, the relationship breaks down. We propose a method for determining $\eta$ within $r_m$ from an observational measurement of $\delta_\alpha$. In cases where $\eta$ and $\delta_\alpha$ can be measured independently, this new method offers a way of measuring the cluster's dynamical state.

The dynamical evolution of accreted star clusters in the Milky Way

We perform $N$-body simulations of star clusters in time-dependant galactic potentials. Since the Milky Way was built-up through mergers with dwarf galaxies, its globular cluster population is made up of clusters formed both during the initial collapse of the Galaxy and in dwarf galaxies that were later accreted. Throughout a dwarf-Milky Way merger, dwarf galaxy clusters are subject to a changing galactic potential. Building on our previous work, we investigate how this changing galactic potential affects the evolution of a cluster's half mass radius. In particular, we simulate clusters on circular orbits around a dwarf galaxy that either falls into the Milky Way or evaporates as it orbits the Milky Way. We find that the dynamical evolution of a star cluster is determined by whichever galaxy has the strongest tidal field at the position of the cluster. Thus, clusters entering the Milky Way undergo changes in size as the Milky Way tidal field becomes stronger and that of the dwarf diminishes. We find that ultimately accreted clusters quickly become the same size as a cluster born in the Milky Way on the same orbit. Assuming their initial sizes are similar, clusters born in the Galaxy and those that are accreted cannot be separated based on their current size alone.

The size of star clusters accreted by the Milky Way

We perform N-body simulations of a cluster that forms in a dwarf galaxy and is then accreted by the Milky Way to investigate how a cluster's structure is affected by a galaxy merger. We find that the cluster's half-mass radius will respond quickly to this change in potential. When the cluster is placed on an orbit in the Milky Way with a stronger tidal field the cluster experiences a sharp decrease in size in response to increased tidal forces. Conversely, when placed on an orbit with a weaker tidal field, the cluster expands since tidal forces decrease and no longer limit the expansion due to internal effects. In all cases, we find that the cluster's half-mass radius will eventually be indistinguishable from a cluster that has always lived in the Milky Way on that orbit. These adjustments occur within 1–2 half-mass relaxation times of the cluster in the dwarf galaxy. We also find this effect to be qualitatively independent of the time that the cluster is taken from the dwarf galaxy. In contrast to the half-mass radius, we show the core radius of the cluster is not affected by the potential the cluster lives in. Our work suggests that structural properties of accreted clusters are not distinct from clusters born in the Milky Way. Other cluster properties, such as metallicity and horizontal branch morphology, may be the only way to identify accreted star clusters in the Milky Way.

Dynamical evolution and spatial mixing of multiple population globular clusters

Numerous spectroscopic and photometric observational studies have provided strong evidence for the widespread presence of multiple stellar populations in globular clusters. In this paper, we study the long-term dynamical evolution of multiple population clusters, focusing on the evolution of the spatial distributions of the first- (FG) and second-generation (SG) stars. In previous studies, we have suggested that SG stars formed from the ejecta of FG AGB stars are expected initially to be concentrated in the cluster inner regions. Here, by means of N-body simulations, we explore the time-scales and the dynamics of the spatial mixing of the FG and the SG populations and their dependence on the SG initial concentration. Our simulations show that, as the evolution proceeds, the radial profile of the SG/FG number ratio, NSG/NFG, is characterized by three regions: (1) a flat inner part; (2) a declining part in which FG stars are increasingly dominant and (3) an outer region where the NSG/NFG profile flattens again (the NSG/NFG profile may rise slightly again in the outermost cluster regions). Until mixing is complete and the NSG/NFG profile is flat over the entire cluster, the radial variation of NSG/NFG implies that the fraction of SG stars determined by observations covering a limited range of radial distances is not, in general, equal to the SG global fraction, (NSG/NFG)glob. The distance at which NSG/NFG equals (NSG/NFG)glob is approximately between 1 and 2 cluster half-mass radii. The time-scale for complete mixing depends on the SG initial concentration, but in all cases complete mixing is expected only for clusters in advanced evolutionary phases, having lost at least 60–70 per cent of their mass due to two-body relaxation (in addition to the early FG loss due to the cluster expansion triggered by SNII ejecta and gas expulsion).The results of our simulations suggest that in many Galactic globular clusters the SG should still be more spatially concentrated than the FG.

On the mass—radius relation of hot stellar systems

Abstract Most globular clusters have half-mass radii of a few pc with no apparent correlation with their masses. This is different from elliptical galaxies, for which the Faber–Jackson relation suggests a strong positive correlation between mass and radius. Objects that are somewhat in between globular clusters and low-mass galaxies, such as ultracompact dwarf galaxies, have a mass–radius relation consistent with the extension of the relation for bright ellipticals. Here we show that at an age of 10 Gyr a break in the mass–radius relation at ∼106 M⊙ is established because objects below this mass, i.e. globular clusters, have undergone expansion driven by stellar evolution and hard binaries. From numerical simulations we find that the combined energy production of these two effects in the core comes into balance with the flux of energy that is conducted across the half-mass radius by relaxation. An important property of this ‘balanced’ evolution is that the cluster half-mass radius is independent of its initial value and is a function of the number of bound stars and the age only. It is therefore not possible to infer the initial mass–radius relation of globular clusters, and we can only conclude that the present day properties are consistent with the hypothesis that all hot stellar systems formed with the same mass–radius relation and that globular clusters have moved away from this relation because of a Hubble time of stellar and dynamical evolution.

THE DYNAMICAL STATE OF THE GLOBULAR CLUSTER M10 (NGC 6254)

Studying the radial variation of the stellar mass function in globular clusters (GCs) has proved a valuable tool to explore the collisional dynamics leading to mass segregation and core collapse. In order to study the radial dependence of the luminosity and mass function of M 10, we used ACS/HST deep high resolution archival images, reaching out to approximately the cluster's half-mass radius (rhm), combined with deep WFPC2 images that extend our radial coverage to more than 2 rhm. From our photometry, we derived a radial mass segregation profile and a global mass function that we compared with those of simulated clusters containing different energy sources (namely hard binaries and/or an IMBH) able to halt core collapse and to quench mass segregation. A set of direct N-body simulations of GCs, with and without an IMBH of mass 1% of the total cluster mass, comprising different initial mass functions (IMFs) and primordial binary fractions, was used to predict the observed mass segregation profile and mass function. The mass segregation profile of M 10 is not compatible with cluster models without either an IMBH or primordial binaries, as a source of energy appears to be moderately quenching mass segregation in the cluster. Unfortunately, the present observational uncertainty on the binary fraction in M10 does not allow us to confirm the presence of an IMBH in the cluster, since an IMBH, a dynamically non-negligible binary fraction (~ 5%), or both can equally well explain the radial dependence of the cluster mass function.

ASPITZER SPACE TELESCOPEATLAS OF ω CENTAURI: THE STELLAR POPULATION, MASS LOSS, AND THE INTRACLUSTER MEDIUM

We present a Spitzer Space Telescope imaging survey of the most massive Galactic globular cluster, ω Centauri, and investigate stellar mass loss at low metallicity and the intracluster medium (ICM). The survey covers approximately 3.2× the cluster half-mass radius at 3.6, 4.5, 5.8, 8, and 24 μm, resulting in a catalog of over 40,000 point sources in the cluster. Approximately 140 cluster members ranging 1.5 dex in metallicity show a red excess at 24 μm, indicative of circumstellar dust. If all of the dusty sources experience mass loss, the cumulative rate of loss is estimated at 2.9-4.2 ×10–7 M ☉ yr–1, 63%-66% of which is supplied by three asymptotic giant branch stars at the tip of the red giant branch (RGB). There is little evidence for strong mass loss lower on the RGB. If this material had remained in the cluster center, its dust component (1 × 10–4 M ☉) would be detectable in our 24 and 70 μm images. While no dust cloud located at the center of ω Cen is apparent, we do see four regions of very faint, diffuse emission beyond two half-mass radii at 24 μm. It is unclear whether these dust clouds are foreground emission or are associated with ω Cen. In the latter case, these clouds may be the ICM in the process of escaping from the cluster.

Star Cluster Life-times: Dependence on Mass, Radius and Environment

AbstractThe dissolution time (tdis) of clusters in a tidal field does not scale with the “classical” expression for the relaxation time. First, the scaling with N, and hence cluster mass, is shallower due to the finite escape time of stars. Secondly, the cluster half-mass radius is of little importance. This is due to a balance between the relative tidal field strength and internal relaxation, which have an opposite effect on tdis, but of similar magnitude. When external perturbations, such as encounters with giant molecular clouds (GMC) are important, tdis for an individual cluster depends strongly on radius. The mean dissolution time for a population of clusters, however, scales in the same way with mass as for the tidal field, due to the weak dependence of radius on mass. The environmental parameters that determine tdis are the tidal field strength and the density of molecular gas. We compare the empirically derived tdis of clusters in six galaxies to theoretical predictions and argue that encounters with GMCs are the dominant destruction mechanism. Finally, we discuss a number of pitfalls in the derivations of tdis from observations, such as incompleteness, with the cluster system of the SMC as particular example.

Why is the mass function of NGC 6218 flat?

We have used the FORS-1 camera on the VLT to study the main sequence (MS) of the globular cluster NGC 6218 in the V and R bands. The observations cover an area of around the cluster centre and probe the stellar population out to the cluster's half-mass radius (). The colour–magnitude diagram (CMD) that we derive in this way reveals a narrow and well defined MS extending down to the detection limit at , or about 6 magnitudes below the turn-off, corresponding to stars of ~0.25 . The luminosity function (LF) obtained with these data shows a marked radial gradient, in that the ratio of lower- and higher-mass stars increases monotonically with radius. The mass function (MF) measured at the half-mass radius, and as such representative of the cluster's global properties, is surprisingly flat. Over the range , the number of stars per unit mass follows a power-law distribution of the type , where, for comparison, Salpeter's IMF would be . We expect that such a flat MF does not represent the cluster's IMF but is the result of severe tidal stripping of the stars from the cluster due to its interaction with the Galaxy's gravitational field. Our results cannot be reconciled with the predictions of recent theoretical models that imply a relatively insignificant loss of stars from NGC 6218 as measured by its expected very long time to disruption. They are more consistent with the orbital parameters based on the Hipparcos reference system that imply a much higher degree of interaction of this cluster with the Galaxy than assumed by those models. Our results indicate that, if the orbit of a cluster is known, the slope of its MF could be useful in discriminating between the various models of the Galactic potential.

The Dynamical State and Blue Straggler Population of the Globular Cluster NGC 6266 (M62)

We have used a proper combination of multiband high-resolution Hubble Space Telescope WFPC2 and wide-field ground-based observations to image the Galactic globular cluster NGC 6266 (M62). The extensive photometric data set allows us to determine the center of gravity and to construct the most extended radial profile ever published for this cluster including, for the first time, detailed star counts in the very inner region. The star density profile is well reproduced by a standard King model with an extended core (∼19'') and a modest value of the concentration parameter (c = 1.5), indicating that the cluster has not yet experienced core collapse. The millisecond pulsar population (whose members are all in binary systems) and the X-ray-emitting population (more than 50 sources within the cluster half-mass radius) suggest that NGC 6266 is in a dynamical phase particularly active in generating binaries through dynamical encounters. UV observations of the central region have been used to probe the population of blue straggler stars, whose origin might be also affected by dynamical interactions. The comparison with other globular clusters observed with a similar strategy shows that the blue straggler content in NGC 6266 is relatively low, suggesting that the formation channel that produces binary systems hosting neutron stars or white dwarfs is not effective in significantly increasing the blue straggler population. Moreover, an anticorrelation between millisecond pulsar content and blue straggler specific frequency in globular clusters seems to be emerging with increasing evidence.

X‐Ray Sources and Their Optical Counterparts in the Globular Cluster M4

We report on the Chandra X-Ray Observatory ACIS-S3 imaging observation of the Galactic globular cluster M4 (NGC 6121). We detect 12 X-ray sources inside the core and 19 more within the cluster half-mass radius. The limiting luminosity of this observation is LX ≈ 1029 ergs s-1 for sources associated with the cluster, the deepest X-ray observation of a globular cluster to date. We identify six X-ray sources with known objects and use ROSAT observations to show that the brightest X-ray source is variable. Archival data from the Hubble Space Telescope allow us to identify optical counterparts to 16 X-ray sources. Based on the X-ray and optical properties of the identifications and the information from the literature, we classify two (possibly three) sources as cataclysmic variables, one X-ray source as a millisecond pulsar, and 12 sources as chromospherically active binaries. Comparison of M4 with 47 Tuc and NGC 6397 suggests a scaling of the number of active binaries in these clusters with the cluster (core) mass.