M31 Globular Cluster System Research Articles

Prospects for Detecting Fast Radio Bursts in the Globular Clusters of Nearby Galaxies

The recent detection of a repeating fast radio burst (FRB) in an old globular cluster in M81 challenges traditional FRB formation mechanisms based on the magnetic activity of young neutron stars formed in core-collapse supernovae. Furthermore, the detection of this repeater in such a nearby galaxy implies a high local universe rate of similar events in globular clusters. Building off the properties inferred from the M81 FRB, we predict the number of FRB sources in nearby (d ≲ 20 Mpc) galaxies with large globular cluster systems known. Incorporating the uncertain burst energy distribution, we estimate the rate of bursts detectable in these galaxies by radio instruments such as FAST and MeerKat. Of all local galaxies, we find M87 is the best candidate for FRB detections. We predict that M87's globular cluster system contains FRB sources at present and that a dedicated radio survey (by either FAST or MeerKat) of hr has a 90% probability of detecting a globular cluster FRB in M87. The detection of even a handful of additional globular cluster FRBs would provide invaluable constraints on FRB mechanisms and population properties. Previous studies have demonstrated young neutron stars formed following the collapse of dynamically formed massive white dwarf binary mergers may provide the most natural mechanism for these bursts. We explore the white dwarf merger scenario using a suite of N-body cluster models, focusing in particular on such mergers in M87's clusters. We describe a number of outstanding features of this scenario that in principle may be testable with an ensemble of observed FRBs in nearby globular clusters.

The globular cluster system of the Auriga simulations

ABSTRACT We investigate whether the galaxy and star formation model used for the Auriga simulations can produce a realistic globular cluster (GC) population. We compare statistics of GC candidate star particles in the Auriga haloes with catalogues of the Milky Way (MW) and Andromeda (M31) GC populations. We find that the Auriga simulations do produce sufficient stellar mass for GC candidates at radii and metallicities that are typical for the MW GC system (GCS). We also find varying mass ratios of the simulated GC candidates relative to the observed mass in the MW and M31 GCSs for different bins of galactocentric radius metallicity (rgal–[Fe/H]). Overall, the Auriga simulations produce GC candidates with higher metallicities than the MW and M31 GCS and they are found at larger radii than observed. The Auriga simulations would require bound cluster formation efficiencies higher than 10 per cent for the metal-poor GC candidates, and those within the Solar radius should experience negligible destruction rates to be consistent with observations. GC candidates in the outer halo, on the other hand, should either have low formation efficiencies, or experience high mass-loss for the Auriga simulations to produce a GCS that is consistent with that of the MW or M31. Finally, the scatter in the metallicity as well as in the radial distribution between different Auriga runs is considerably smaller than the differences between that of the MW and M31 GCSs. The Auriga model is unlikely to give rise to a GCS that can be consistent with both galaxies.

The outer halo globular cluster system of M31 – III. Relationship to the stellar halo

We utilise the final catalogue from the Pan-Andromeda Archaeological Survey to investigate the links between the globular cluster system and field halo in M31 at projected radii $R_p=25-150$ kpc. In this region the cluster radial density profile exhibits a power-law decline with index $\Gamma=-2.37\pm0.17$, matching that for the stellar halo component with [Fe/H] $<-1.1$. Spatial density maps reveal a striking correspondence between the most luminous substructures in the metal-poor field halo and the positions of many globular clusters. By comparing the density of metal-poor halo stars local to each cluster with the azimuthal distribution at commensurate radius, we reject the possibility of no correlation between clusters and field overdensities with high confidence. We use our stellar density measurements and previous kinematic data to demonstrate that $\approx35-60\%$ of clusters exhibit properties consistent with having been accreted into the outskirts of M31 at late times with their parent dwarfs. Conversely, at least $\sim40\%$ of remote clusters show no evidence for a link with halo substructure. The radial density profile for this subgroup is featureless and closely mirrors that observed for the apparently smooth component of the metal-poor stellar halo. We speculate that these clusters are associated with the smooth halo; if so, their properties appear consistent with a scenario where the smooth halo was built up at early times via the destruction of primitive satellites. In this picture the entire M31 globular cluster system outside $R_p=25$ kpc comprises objects accumulated from external galaxies over a Hubble time of growth.

The kinematics of globular clusters systems in the outer halos of the Aquarius simulations

Stellar halos and globular cluster (GC) systems contain valuable information regarding the assembly history of their host galaxies. Motivated by the detection of a significant rotation signal in the outer halo GC system of M31, we investigate the likelihood of detecting such a rotation signal in projection, using cosmological simulations. To this end we select subsets of tagged particles in the halos of the Aquarius simulations to represent mock GC systems, and analyse their kinematics. We find that GC systems can exhibit a non-negligible rotation signal provided the associated stellar halo also has a net angular momentum. The ability to detect this rotation signal is highly dependent on the viewing perspective, and the probability of seeing a signal larger than that measured in M31 ranges from 10% to 90% for the different halos in the Aquarius suite. High values are found from a perspective such that the projected angular momentum of the GC system is within 40 deg of the rotation axis determined via the projected positions and line-of-sight velocities of the GCs. Furthermore, the true 3D angular momentum of the outer stellar halo is relatively well aligned, within 35 deg, with that of the mock GC systems. We argue that the net angular momentum in the mock GC systems arises naturally when the majority of the material is accreted from a preferred direction, namely along the dominant dark matter filament of the large-scale structure that the halos are embedded in. This, together with the favourable edge-on view of M31's disk suggests that it is not a coincidence that a large rotation signal has been measured for its outer halo GC system.

Stellar halos and the link to galaxy formation

AbstractI present a brief overview of how stellar halos may be used to constrain the process of galaxy formation. In particular, streams and substructure in stellar halos trace merger events but can also be used to determine the mass distribution of the host galaxy and hence put constraints on the nature of dark matter. Much of the focus of this contribution is on the Milky Way, but I also present an attempt to understand the kinematics of the globular cluster system of M31.

NONLINEAR COLOR-METALLICITY RELATIONS OF GLOBULAR CLUSTERS. V. NONLINEAR ABSORPTION-LINE INDEX VERSUS METALLICITY RELATIONS AND BIMODAL INDEX DISTRIBUTIONS OF M31 GLOBULAR CLUSTERS

Recent spectroscopy on the globular cluster (GC) system of M31 with unprecedented precision witnessed a clear bimodality in absorption-line index distributions of old GCs. Such division of extragalactic GCs, so far asserted mainly by photometric color bimodality, has been viewed as the presence of merely two distinct metallicity subgroups within individual galaxies and forms a critical backbone of various galaxy formation theories. Given that spectroscopy is a more detailed probe into stellar population than photometry, the discovery of index bimodality may point to the very existence of dual GC populations. However, here we show that the observed spectroscopic dichotomy of M31 GCs emerges due to the nonlinear nature of metallicity-to-index conversion and thus one does not necessarily have to invoke two separate GC subsystems. We take this as a close analogy to the recent view that metallicity-color nonlinearity is primarily responsible for observed GC color bimodality. We also demonstrate that the metallicity-sensitive magnesium line displays non-negligible metallicity-index nonlinearity and Balmer lines show rather strong nonlinearity. This gives rise to bimodal index distributions, which are routinely interpreted as bimodal metallicity distributions, not considering metallicity-index nonlinearity. Our findings give a new insight into the constitution of M31's GC system, which could change much of the current thought on the formation of GC systems and their host galaxies.

Metal abundance and kinematical properties of the M81 globular cluster system

We present metal abundance properties of 144 globular clusters associated with M81. These globulars represent the largest globular cluster sample in M81 till now. Our main results are: the distribution of metallicities is bimodal, with metallicity peaks at [Fe/H] ≈ −1.51 and −0.58, and the metal-poor globular clusters tend to be less spatially concentrated than the metal-rich ones; the metal-rich globular clusters in M81 do not demonstrate a centrally concentrated spatial distribution like the metal-rich ones in M31 do; like our Galaxy and M31, the globular clusters in M81 have a small radial metallicity gradient. These results are consistent with those obtained from a small sample of M81 globular clusters. In addition, this paper shows that there is evidence that a strong rotation of the M81 globular cluster system around the minor axis exists, and that rotation is present in the metal-rich globular cluster subsample, but the metal-poor globular cluster subsample shows no evidence of rotation. The most significant difference between the rotation of the metal-rich and metal-poor globular clusters occurs at intermediate projected galactocentric radii. Our results confirm the conclusion of Schroder et al. that M81's metal-rich globular clusters at intermediate projected radii are associated with a thick disk of M81.

CLASSICAL NOVAE IN ANDROMEDA: LIGHT CURVES FROM THE PALOMAR TRANSIENT FACTORY ANDGALEX

We present optical light curves of twenty-nine novae in M31 during the 2009 and 2010 observing seasons of the Palomar Transient Factory (PTF). The dynamic and rapid cadences in PTF monitoring of M31, from one day to even ten minutes, provide excellent temporal coverage of nova light curves, enabling us to record the photometric evolution of M31 novae in unprecedented detail. We also detect eight of these novae in the near ultraviolet (UV) band with the Galaxy Evolution Explorer (GALEX) satellite. Novae M31N2009-10b and 2010-11a show prominent UV emission peaking a few days prior to their optical maxima, possibly implying aspherical outbursts. Additionally, our blue-shifted spectrum of the recent outburst of PT And (M31N2010-12a) indicates that it is a recurrent nova in M31 and not a dwarf nova in the Milky Way as was previously assumed. Finally, we systematically searched for novae in all confirmed globular clusters of M31 and found only M31N 2010-10f associated with Bol 126. The specific nova rate in the M31 globular cluster system is thus about one per year which is not enhanced relative to the rate outside the globular cluster system.

WIDE-FIELD PRECISION KINEMATICS OF THE M87 GLOBULAR CLUSTER SYSTEM

We present the most extensive combined photometric and spectroscopic study to date of the enormous globular cluster (GC) system around M87, the central giant elliptical galaxy in the nearby Virgo cluster. Using observations from DEIMOS and LRIS at Keck, and Hectospec on the MMT, we derive new, precise radial velocities for 451 GCs around M87, with projected radii from ~ 5 to 185 kpc. We combine these measurements with literature data for a total sample of 737 objects, which we use for a re-examination of the kinematics of the GC system of M87. The velocities are analyzed in the context of archival wide-field photometry and a novel Hubble Space Telescope catalog of half-light radii, which includes sizes for 344 spectroscopically confirmed clusters. We use this unique catalog to identify 18 new candidate ultra-compact dwarfs, and to help clarify the relationship between these objects and true GCs. We find much lower values for the outer velocity dispersion and rotation of the GC system than in earlier papers, and also differ from previous work in seeing no evidence for a transition in the inner halo to a potential dominated by the Virgo cluster, nor for a truncation of the stellar halo. We find little kinematical evidence for an intergalactic GC population. Aided by the precision of the new velocity measurements, we see significant evidence for kinematical substructure over a wide range of radii, indicating that M87 is in active assembly. A simple, scale-free analysis finds less dark matter within ~85 kpc than in other recent work, reducing the tension between X-ray and optical results. In general, out to a projected radius of ~ 150 kpc, our data are consistent with the notion that M87 is not dynamically coupled to the Virgo cluster; the core of Virgo may be in the earliest stages of assembly.

STAR CLUSTERS IN M31. II. OLD CLUSTER METALLICITIES AND AGES FROM HECTOSPEC DATA

We present new high signal-to-noise spectroscopic data on the M31 globular cluster (GC) system, obtained with the Hectospec multifiber spectrograph on the 6.5 m MMT. More than 300 clusters have been observed at a resolution of 5 Å and with a median S/N of 75 per Å, providing velocities with a median uncertainty of 6 km s−1. The primary focus of this paper is the determination of mean cluster metallicities, ages, and reddenings. Metallicities were estimated using a calibration of Lick indices with [Fe/H] provided by Galactic GCs. These match well the metallicities of 24 M31 clusters determined from Hubble Space Telescope color–magnitude diagrams, the differences having an rms of 0.2 dex. The metallicity distribution is not generally bimodal, in strong distinction with the bimodal Galactic globular distribution. Rather, the M31 distribution shows a broad peak, centered at [Fe/H] = −1, possibly with minor peaks at [Fe/H] = −1.4, −0.7, and −0.2, suggesting that the cluster systems of M31 and the Milky Way had different formation histories. Ages for clusters with [Fe/H] > − 1 were determined using the automatic stellar population analysis program EZ_Ages. We find no evidence for massive clusters in M31 with intermediate ages, those between 2 and 6 Gyr. Moreover, we find that the mean ages of the old GCs are remarkably constant over about a decade in metallicity (−0.95≲ [Fe/H] ≲0.0).

DYNAMICAL EVOLUTION OF THE M87 GLOBULAR CLUSTER SYSTEM

We study the dynamical evolution of the M87 globular cluster (GC) system using the most advanced and realistic Fokker-Planck (FP) model.By comparing our FP models with both mass function (MF) and radial distribution (RD) of the observed GC system, we find the best-fit initial (at M87's age of 2-3 Gyr) MF and RD for three GC groups: all GCs, blue GCs, and red GCs. We estimate the initial total mass in GCs to be <TEX>$1.8^{+0.3}_{-0.2}{\times}10^{10}M_{\bigodot}$</TEX>, which is about 100 times larger than that of the Milky Way GC system. We also find that the fraction of the total mass currently in GCs is 34\%. When blue and red GCs are fitted separately, blue GCs initially have a larger total mass and a shallower radial distribution than red GCs. If one assumes that most of the significant major merger events of M87 have ended by the age of 2-3 Gyr, our finding that blue (metal-poor) GCs initially had a shallower radial distribution supports the major merger scenario for the origin of metallicity bimodality.

A systematic study of low-mass X-ray binaries in the M31 globular cluster system

We investigate low-mass X-ray binaries (LMXBs) in the M31 globular cluster (GC) system using data from the 2XMMi catalogue. These X-ray data are based on all publicly available XMM–Newton observations of the galaxy. This new survey provides the most complete and homogeneous X-ray survey of M31's GCs to date, covering >80 per cent of the confirmed old clusters in the galaxy. We associate 41 X-ray sources with confirmed old clusters in the recent M31 cluster catalogue of Peacock et al. Comparing these data with previous surveys of M31, it is found that three of these clusters are newly identified, including a bright transient source in the cluster B128. Four additional clusters, that are not detected in these 2XMMi data, have previously been associated with X-ray sources from Chandra or ROSAT observations. Including these clusters, we identify 45 clusters in M31 which are associated with X-ray emission. By considering the latest optical GC catalogues, we identify that three of the previously proposed X-ray clusters are likely to be background galaxies and two have stellar profiles. We consider the properties of LMXB hosting clusters and confirm significant trends between the presence of an LMXB and the metallicity, luminosity and stellar collision rate of a cluster. We consider the relationship between the luminosity and stellar collision rate of a cluster and note that LMXB hosting clusters have higher than average stellar collision rates for their mass. Our findings strongly suggest that the stellar collision rate is the dominant parameter related to the presence of LMXBs. This is consistent with the formation of LMXBs in GCs through dynamical interactions.

EVIDENCE FOR AN ACCRETION ORIGIN FOR THE OUTER HALO GLOBULAR CLUSTER SYSTEM OF M31

We use a sample of newly-discovered globular clusters from the Pan-Andromeda Archaeological Survey (PAndAS) in combination with previously-catalogued objects to map the spatial distribution of globular clusters in the M31 halo. At projected radii beyond ~30 kpc, where large coherent stellar streams are readily distinguished in the field, there is a striking correlation between these features and the positions of the globular clusters. Adopting a simple Monte Carlo approach, we test the significance of this association by computing the probability that it could be due to the chance alignment of globular clusters smoothly distributed in the M31 halo. We find the likelihood of this possibility is low, below 1%, and conclude that the observed spatial coincidence between globular clusters and multiple tidal debris streams in the outer halo of M31 reflects a genuine physical association. Our results imply that the majority of the remote globular cluster system of M31 has been assembled as a consequence of the accretion of cluster-bearing satellite galaxies. This constitutes the most direct evidence to date that the outer halo globular cluster populations in some galaxies are largely accreted.

SPECTROSCOPY OF M81 GLOBULAR CLUSTERS

We obtained spectra of 74 globular clusters in M81. These globular clusters had been identified as candidates in an HST ACS I-band survey. 68 of these 74 clusters lie within 7' of the M81 nucleus. 62 of these clusters are newly spectroscopically confirmed, more than doubling the number of confirmed M81 GCs from 46 to 108. We determined metallicities for our 74 observed clusters using an empirical calibration based on Milky Way globular clusters. We combined our results with 34 M81 globular cluster velocities and 33 metallicities from the literature and analyzed the kinematics and metallicity of the M81 globular cluster system. The mean of the total sample of 107 metallicities is -1.06 +/- 0.07, higher than either M31 or the Milky Way. We suspect this high mean metallicity is due to an overrepresentation of metal-rich clusters in our sample created by the spatial limits of the HST I-band survey. The metallicity distribution shows marginal evidence for bimodality, with metal-rich and metal-poor peaks approximately matching those of M31 and the Milky Way. The GC system as a whole, and the metal-poor GCs alone, show evidence of a radial metallicity gradient. The M81 globular cluster system as a whole shows strong evidence of rotation, with V_r(deprojected) = 108 +/- 22 km/s overall. This result is likely biased toward high rotational velocity due to overrepresentation of metal-rich, inner clusters. The rotation patterns among globular cluster subpopulations are roughly similar to those of the Milky Way: clusters at small projected radii and metal-rich clusters rotate strongly, while clusters at large projected radii and metal-poor clusters show weaker evidence of rotation.

Origin of rotational kinematics in the globular cluster system of M31: a new clue to the bulge formation

Abstract We propose that the rotational kinematics of the globular cluster system (GCS) in M31 can result from a past major merger event that could have formed its bulge component. We numerically investigate kinematical properties of globular clusters (GCs) in remnants of galaxy mergers between two discs with GCs in both their disc and halo components. We find that the GCS formed during major merging can show strongly rotational kinematics with the maximum rotational velocities of ∼140–170 km s−1 for a certain range of orbital parameters of merging. We also find that a rotating stellar bar, which can be morphologically identified as a boxy bulge if seen edge-on, can be formed in models for which the GCSs show strongly rotational kinematics. We thus suggest that the observed rotational kinematics of GCs with different metallicities in M31 can be closely associated with the ancient major merger event. We discuss whether the formation of the rotating bulge/bar in M31 can be due to the ancient merger.

STRUCTURAL PARAMETERS OF THE MESSIER 87 GLOBULAR CLUSTERS

We derive structural parameters for ~2000 globular clusters in the giant Virgo elliptical M87 using extremely deep Hubble Space Telescope images in F606W (V) and F814W (I) taken with the ACS/WFC. The cluster scale sizes (half-light radii r_h) and ellipticities are determined from PSF-convolved King-model profile fitting. We find that the r_h distribution closely resembles the inner Milky Way clusters, peaking at r_h~2.5 pc and with virtually no clusters more compact than r_h ~ 1 pc. The metal-poor clusters have on average an r_h 24% larger than the metal-rich ones. The cluster scale size shows a gradual and noticeable increase with galactocentric distance. Clusters are very slightly larger in the bluer waveband V a possible hint that we may be beginning to see the effects of mass segregation within the clusters. We also derived a color magnitude diagram for the M87 globular cluster system which show a striking bimodal distribution.

A new era for the globular cluster system in M31

AbstractThe globular cluster system in M31 is an ideal laboratory for studying the formation and evolution of M31 as well as the globular clusters themselves. There have been numerous surveys and studies of the globular clusters in M31. However, only recently has the entire body of M31 been searched for globular clusters using wide-field CCD images by our group. A new era for the M31 globular cluster system has begun with the advent of wide-field CCD surveys of M31. We have discovered more than 100 new globular clusters in M31. Our catalog currently includes more than 500 globular clusters confirmed either based on spectra or HST images, many more than in the Milky Way. We present the structure, kinematics and chemical abundance of the M31 globular cluster system based on this large sample, and the implications for the formation and evolution of M31.

The outer halo globular clusters of M31

Abstract We present Keck/HIRES spectra of three globular clusters in the outer halo of M31, at projected distances beyond ≈80 kpc from M31. The measured recession velocities for all three globular clusters confirm their association with the globular cluster system of M31. We find evidence for a declining velocity dispersion with radius for the globular cluster system. Their measured internal velocity dispersions, derived virial masses and mass-to-light ratios are consistent with those for the bulk of the M31 globular cluster system. We derive old ages and metallicities which indicate that all three belong to the metal-poor halo globular cluster subpopulation. We find indications that the radial gradient of the mean metallicity of the globular cluster system interior to 50 kpc flattens in the outer regions, however it is still more metal-poor than the corresponding field stars at the same (projected) radius.

COLOR BIMODALITY IN M87 GLOBULAR CLUSTERS

We present an analysis of a 50 orbit HST ACS observation of the M87 globular cluster system. We use the extraordinary depth of this dataset to test whether the colors and magnitudes show evidence for a mass-metallicity relation in globular cluster populations. We find only a weak or absent relation between the colors and magnitudes of the metal poor subpopulation of globular clusters. The weakness or absence of a color-magnitude relation is established over a wide range in luminosity from $M_V=-11$ to $M_V=-6$, encompassing most of the M87 globular clusters. The constancy of the colors of the metal-poor subpopulation seen in our 50 orbit observation is in contrast to suggestions from single orbit ACS data that the metal-poor globular clusters in M87 and several other galaxies show a "blue tilt." The formal best fit for the mass-metallicity relation for the metal-poor subpopulation in our much deeper data is $Z\propto M^{0.08\pm0.05}$. Our analysis of these data also shows a possible small "red tilt" in the metal-rich globular cluster subpopulation. While either of these small tilts may be real, they may also illustrate the limit to which mass-metallicity relations can be determined, even in such extraordinarily deep data. We specifically test for a wide range of systematic effects and find that while small tilts cannot be confirmed or rejected, the data place a strong upper limit to any tilt of $|0.20|\pm0.05$. This upper limit is much smaller than some earlier claims from single orbit data, and strongly limits self-enrichment within globular clusters. This mass-metallicity relation for globular clusters is also shallower than the relation for galaxies, suggesting that the formation mechanisms for these two types of objects are different.

The ACS Virgo Cluster Survey. XV. The Formation Efficiencies of Globular Clusters in Early‐Type Galaxies: The Effects of Mass and Environment

The fraction of stellar mass contained in globular clusters (GCs), also measured by number as the specific frequency, is a fundamental quantity that reflects both a galaxy's early star formation and its entire merging history. We present specific frequencies, luminosities, and mass fractions for the globular cluster systems of 100 early-type galaxies in the ACS Virgo Cluster Survey, the largest homogeneous catalog of its kind. We find the following: (1) GC mass fractions can be high in both giants and dwarfs but are universally low in galaxies with intermediate luminosities. (2) The behavior of specific frequency across galaxy mass is dominated by the blue GCs. (3) GC fractions of low-mass galaxies exhibit a dependence on environment. Nearly all dwarf galaxies with high GC fractions are within 1 Mpc of the cD galaxy M87, presenting the first strong evidence that GC formation in dwarfs is biased toward dense environments. (4) GC formation in central dwarfs is biased because their stars form earliest and most intensely. Comparisons to the Millennium Simulation show that central dwarfs have older stellar populations and form more stars at higher star formation rates (SFRs) and SFR surface densities. The SFR surface density in simulated dwarfs peaks before the total SFR, naturally producing GC populations that are older and more metal-poor than the field stars. (5) Dwarfs within ~40 kpc of the giant ellipticals M87 and M49 are red and have few or no GCs, suggesting that they have been tidally stripped and have contributed their GCs to the halos of their giant neighbors. The central dwarfs with high GC mass fractions are thus likely to be the survivors most similar to the protogalaxies that assembled the rich M87 globular cluster system.