Populations In Globular Clusters Research Articles

NGC 6388 reloaded: some like it hot, but not too much. New constraints on the first-generation polluters

Multiple stellar populations in globular clusters (GCs) are defined and recognized by their chemical signature. Second-generation stars show the effects of nucleosynthesis in the more massive stars of the earliest component that formed in the first star formation burst. High-temperature H-burning produces the whole pattern of (anti-) correlations in proton-capture elements that are widely found in GCs. However, it is still debated where this burning occurred. Here we introduce new powerful diagnostic plots to detect evidence (if any) of products from proton-capture reactions that occur at very high temperatures. To test these detectors of high-temperature H-burning plots, we show that stringent constraints can be placed on the temperature range of the first-generation polluters that contributed to shaping the chemistry of multiple stellar population in the massive bulge GC NGC 6388. Using the largest sample to date (185 stars) of giants with detailed abundance ratios in a single GC (except ω Cen), we may infer that the central temperature of part of the polluters must have been comprised between ∼100 and ∼150 MK if we consider hydrostatic H-burning in the core of massive stars. A much more narrow range (110−120 MK) is inferred if the polluters can be identified in massive asymptotic giant branch stars.

H-band discovery of additional second-generation stars in the Galactic bulge globular cluster NGC 6522 as observed by APOGEE and Gaia

We present an elemental abundance analysis of high-resolution spectra for five giant stars spatially located within the innermost regions of the bulge globular cluster NGC 6522 and derive Fe, Mg, Al, C, N, O, Si, and Ce abundances based on H-band spectra taken with the multi-object APOGEE-north spectrograph from the SDSS-IV Apache Point Observatory Galactic Evolution Experiment (APOGEE) survey. Of the five cluster candidates, two previously unremarked stars are confirmed to have second-generation (SG) abundance patterns, with the basic pattern of depletion in C and Mg simultaneous with enrichment in N and Al as seen in other SG globular cluster populations at similar metallicity. In agreement with the most recent optical studies, the NGC 6522 stars analyzed exhibit (when available) only mild overabundances of the s-process element Ce, contradicting the idea that NGC 6522 stars are formed from gas enriched by spinstars and indicating that other stellar sources such as massive AGB stars could be the primary polluters of intra-cluster medium. The peculiar abundance signatures of SG stars have been observed in our data, confirming the presence of multiple generations of stars in NGC 6522.

X-Ray Binary Luminosity Function Scaling Relations for Local Galaxies Based on Subgalactic Modeling

Abstract We present new Chandra constraints on the X-ray luminosity functions (XLFs) of X-ray binary (XRB) populations, as well as their scaling relations, for a sample of 38 nearby galaxies (D = 3.4–29 Mpc). Our galaxy sample is drawn primarily from the Spitzer Infrared Nearby Galaxies Survey (SINGS) and contains a wealth of Chandra (5.8 Ms total) and multiwavelength data, allowing for star formation rates (SFRs) and stellar masses (M ⋆) to be measured on subgalactic scales. We divided the 2478 X-ray-detected sources into 21 subsamples in bins of specific SFR (sSFR ≡ SFR/M ⋆) and constructed XLFs. To model the XLF dependence on sSFR, we fit a global XLF model, containing contributions from high-mass XRBs (HMXBs), low-mass XRBs (LMXBs), and background sources from the cosmic X-ray background that respectively scale with SFR, M ⋆, and sky area. We find an HMXB XLF that is more complex in shape than previously reported and an LMXB XLF that likely varies with sSFR, potentially due to an age dependence. When applying our global model to XLF data for each individual galaxy, we discover a few galaxy XLFs that significantly deviate from our model beyond statistical scatter. Most notably, relatively low-metallicity galaxies have an excess of HMXBs above ≈1038 erg s−1, and elliptical galaxies that have relatively rich populations of globular clusters (GCs) show excesses of LMXBs compared to the global model. Additional modeling of how the XRB XLF depends on stellar age, metallicity, and GC specific frequency is required to sufficiently characterize the XLFs of galaxies.

Globular clusters as probes of dark matter cusp-core transformations

Abstract Bursty star formation in dwarf galaxies can slowly transform a steep dark matter cusp into a constant density core. We explore the possibility that globular clusters (GCs) retain a dynamical memory of this transformation. To test this, we use the nbody6df code to simulate the dynamical evolution of GCs, including stellar evolution, orbiting in static and time-varying potentials for a Hubble time. We find that GCs orbiting within a cored dark matter halo, or within a halo that has undergone a cusp-core transformation, grow to a size that is substantially larger (Reff > 10 pc) than those in a static cusped dark matter halo. They also produce much less tidal debris. We find that the cleanest signal of an historic cusp-core transformation is the presence of large GCs with tidal debris. However, the effect is small and will be challenging to observe in real galaxies. Finally, we qualitatively compare our simulated GCs with the observed GC populations in the Fornax, NGC 6822, IKN, and Sagittarius dwarf galaxies. We find that the GCs in these dwarf galaxies are systematically larger (〈Reff〉 ≃ 7.8 pc), and have substantially more scatter in their sizes than in situ metal-rich GCs in the Milky Way and young massive star clusters forming in M83 (〈Reff〉 ≃ 2.5 pc). We show that the size, scatter, and survival of GCs in dwarf galaxies are all consistent with them having evolved in a constant density core, or a potential that has undergone a cusp-core transformation, but not in a dark matter cusp.

Kinematical evolution of multiple stellar populations in star clusters

ABSTRACT We present the results of a suite of N-body simulations aimed at understanding the fundamental aspects of the long-term evolution of the internal kinematics of multiple stellar populations in globular clusters. Our models enable us to study the cooperative effects of internal, relaxation-driven processes and external, tidally induced perturbations on the structural and kinematic properties of multiple-population globular clusters. To analyse the dynamical behaviour of the multiple stellar populations in a variety of spin-orbit coupling conditions, we have considered three reference cases in which the tidally perturbed star cluster rotates along an axis oriented in different directions with respect to the orbital angular momentum vector. We focus specifically on the characterization of the evolution of the degree of differential rotation and anisotropy in the velocity space, and we quantify the process of spatial and kinematic mixing of the two populations. In light of recent and forthcoming explorations of the internal kinematics of this class of stellar systems by means of line-of sight and astrometric measurements, we also investigate the implications of projection effects and spatial distribution of the stars adopted as tracers. The kinematic and structural richness emerging from our models further emphasizes the need and the importance of observational studies aimed at building a complete kinematical picture of the multiple population phenomenon.

The Relationship between Globular Cluster Mass, Metallicity, and Light-element Abundance Variations

We investigate aluminum abundance variations in the stellar populations of globular clusters using both literature measurements of sodium and aluminum and APOGEE measurements of nitrogen and aluminum abundances. For the latter, we show that the Payne is the most suitable of the five available abundance pipelines for our purposes. Our combined sample of 42 globular clusters spans approximately 2 dex in [Fe/H] and 1.5 dex in $\log{M_{GC}/M_{\odot}}$. We find no fewer than five globular clusters with significant internal variations in nitrogen and/or sodium with little-to-no corresponding variation in aluminum, and that the minimum present-day cluster mass for aluminum enrichment in metal-rich systems is $\log{M_{GC}/M_{\odot}} \approx 4.50 + 2.17(\rm{[Fe/H]}+1.30)$. We demonstrate that the slopes of the [Al/Fe] vs [Na/Fe] and [Al/Fe] vs [N/Fe] relations for stars without field-like abundances are approximately log-linearly dependent on both the metallicity and the stellar mass of the globular clusters. In contrast, the relationship between [Na/Fe] and [N/Fe] shows no evidence of such dependencies. This suggests that there were (at least) two classes of non-supernovae chemical polluters that were common in the early universe, and that their relative contributions within globular clusters somehow scaled with the metallicity and mass of globular clusters. The first of these classes is predominantly responsible for the CNO and NeNa abundance variations, and likewise the second for the MgAl abundance variations. 47 Tuc and M4 are particularly striking examples of this dichotomy. As an auxiliary finding, we argue that abundance variations among Terzan 5 stars are consistent with it being a normal globular cluster.

The evolution of the UV luminosity function of globular clusters in the E-MOSAICS simulations

We present the evolution of the rest-frame ultraviolet (UV) properties of the globular cluster (GC) populations and their host galaxies formed in the E-MOSAICS suite of cosmological hydrodynamical simulations. We compute the luminosities of all clusters associated with 25 simulated Milky Way-mass galaxies, discussed in previous works, in the rest-frame UV and optical bands by combining instantaneous cluster properties (age, mass, metallicity) with simple stellar population models, from redshifts $z=0$ to 10. Due to the rapid fading of young stellar populations in the UV, most of the simulated galaxies do not host GCs bright enough to be individually identified in deep Hubble Space Telescope (HST) observations, even in highly magnified systems. The median age of the most UV-luminous GCs is $<10$ Myr (assuming no extinction), increasing to $\gtrsim 100$ Myr for red optical filters. We estimate that these GCs typically only contribute a few per cent of the total UV luminosity of their host galaxies at any epoch. We predict that the number density of UV-bright proto-GCs (or cluster clumps) will peak between redshifts $z=1-3$. In the main progenitors of Milky Way-mass galaxies, $10$-$20$ per cent of the galaxies at redshifts $1 \lesssim z \lesssim 3$ have clusters brighter than $M_{\rm UV} < -15$, and less than $10$ per cent at other epochs. The brightest cluster in the galaxy sample at $z>2$ is typically $M_{\rm UV} \sim -16$, consistent with the luminosities of observed compact, high-redshift sources.

On the Orbital Decay of Globular Clusters in NGC 1052-DF2: Testing a Baryon-only Mass Model

Abstract The dark matter content of the ultra-diffuse galaxy NGC 1052-DF2, as inferred from globular cluster (GC) and stellar kinematics, carries a considerable amount of uncertainty, with current constraints also allowing for the complete absence of dark matter. We test the viability of such a scenario by examining whether in a “baryon-only” mass model the observed GC population experiences rapid orbital decay due to dynamical friction. Using a suite of 50 multi-GC N-body simulations that match observational constraints on both the stellar component of NGC 1052-DF2 and its GC population but differ in the initial line-of-sight positions and the tangential velocities of the GCs, we show that there is a substantial amount of realization-to-realization variance in the evolution of the GCs. Nevertheless, over ∼10 Gyr, some of the GCs experience significant orbital evolution. Others evolve less. A combination of reduced dynamical friction in the galaxy core and GC–GC scattering keeps the GCs afloat, preventing them from sinking all the way to the galaxy center. While the current phase-space coordinates of the GCs are not unlikely for a baryon-only mass model, the GC system does evolve over time. Therefore, if NGC 1052-DF2 has no dark matter, some of its GCs must have formed farther out, and the GC system must have been somewhat more extended in the past. The presence of a low-mass cuspy halo, while allowed by the kinematics, seems improbable, as significantly shorter inspiral timescales in the central region would quickly lead to the formation of a nuclear star cluster.

Thirteen facts that you need to know on multiple populations in globular clusters

AbstractI review the methods, mostly developed in the last decade, that are commonly used to identify and characterize multiple populations (MPs) in Globular Clusters (GCs) based on photometry. I summarize the results from the recent surveys of MPs with the Hubble Space Telescope (HST) and ground-based facilities and provide a list of the main properties of MPs as inferred from these studies.

A unified construction of stellar evolution and chemical evolution models for the multiple populations in globular clusters

AbstractRecent investigations of multiple stellar populations in globular clusters (GCs) suggest that the horizontal-branch (HB) morphology and mean period of type ab RR Lyrae variables are mostly sensitive to helium abundance, while the star formation timescale has the greatest effect on our chemical evolution model constructed to reproduce the Na-O anti-correlation of GCs. Therefore, by combining the results from synthetic HB model with those from chemical evolution model, we could put better constraints on star formation history and chemical evolution in GCs with multiple populations. From such efforts made for four GCs, M4, M5, M15, and M80, we find that consistent results can be obtained from these two independent models.

Spectroscopic studies of stellar populations in globular clusters and field stars: Implications for globular cluster and Milky Way halo formation

AbstractWe review spectroscopic results concerning multiple stellar populations in globular clusters. The cluster initial mass is the most important parameter determining the fraction of second generation stars. The threshold for the onset of the multiple population phenomenon is 1–3×105 M⊙. Nucleosynthesis is influenced by metallicity: Na/O and Mg/Al anti-correlations are more extended in metal-poor than in metal-rich clusters. Massive clusters are more complex systems than the smaller ones, with several populations characterized by different chemical compositions. The high Li abundance observed in the intermediate second generation stars strongly favours intermediate mass AGB stars as polluters for this class of stars; however, it is well possible that the polluters of extreme second generation stars, that often do not have measurable Li, may be fast rotating massive stars or super-massive stars. The mass budget factor should be a function of the cluster mass, and needs to be large only in massive clusters.

Multiple populations in globular clusters: New insights from chemical evolution and horizontal-branch models

AbstractIn order to investigate the origin of multiple populations in globular clusters (GCs), we have constructed new chemical evolution models for proto-GCs where the supernova blast waves undergo blowout without expelling the ambient gas. Chemical enrichments in our models are then dictated by the winds of massive stars together with the asymptotic-giant-branch stars ejecta. We find that the observed Na-O anti-correlation can be reproduced when multiple episodes of starburst and enrichment are allowed to continue in proto-GCs. The “mass budget problem” is mostly resolved by our models without ad-hoc assumptions on star formation efficiency, initial mass function, and significant loss of first-generation stars. Interestingly, ages and chemical abundances predicted by this chemical evolution model are in good agreements with those independently obtained from our stellar evolution model for the horizontal-branch. We also discuss observational evidence for the GC-like multiple populations in the Milky Way bulge.

Kinematics of multiple stellar populations in Globular Clusters with Gaia

AbstractThe internal dynamics of multiple stellar populations in Globular Clusters (GCs) provides unique constraints on the physical processes responsible for their formation. Specifically, the present-day kinematics of cluster stars, such as rotation and velocity dispersion, seems to be related to the initial configuration of the system. In recent work (Milone et al. 2018), we analyzed for the first time the kinematics of the different stellar populations in NGC 0104 (47 Tucanae) over a large field of view, exploiting the Gaia Data Release 2 proper motions combined with multi-band ground-based photometry. In this paper, based on the work by Cordoni et al. (2019), we extend this analysis to six GCs, namely NGC 0288, NGC 5904 (M 5), NGC 6121 (M 4), NGC 6752, NGC 6838 (M 71) and further explore NGC 0104. Among the analyzed clusters only NGC 0104 and NGC 5904 show significant rotation on the plane of the sky. Interestingly, multiple stellar populations in NGC 5904 exhibit different rotation curves.

Supermassive stars as the origin of the multiple populations in globular clusters

AbstractGlobular clusters (GCs) display anomalous light-elements abundances (HeCNONaMgAl), resembling the yields of hot-hydrogen burning, but there is no consensus yet on the origin of these ubiquitous multiple populations. We present a model in which a super-massive star (SMS, ≳103 M⊙) forms via stellar collisions during GC formation and pollutes the intra-cluster medium. The growth of the SMS finds a balance with the wind mass loss rate, such that the SMS can produce a significant fraction of the total GC mass in processed material, thereby overcoming the so-called mass-budget problem that plagues other models. Because of continuous rejuvenation, the SMS acts as a ‘conveyer-belt’ of hot-hydrogen burning yields with (relatively) low He abundances, in agreement with empirical constraints. Additionally, the amount of processed material per unit of GC mass correlates with GC mass, addressing the specific mass budget problem. We discuss uncertainties and tests of this new self-enrichment scenario.

Internal kinematics of multiple stellar populations in globular clusters

AbstractSpectroscopy and photometry have revealed existence, complexity and properties of the multiple stellar populations (mPOPs) hosted in Galactic globular clusters. However, the conundrum of the formation and evolution of mPOPs is far from being completely exploited: the available pieces of information seem not enough to shed light on these topics. Astrometry, and in particular high-precision proper motions, can provide us the sought-after answers about how mPOPs formed and have evolved in these ancient stellar systems. In the following, I present a brief overview of the observational results on the internal kinematics of the mPOPs in some GCs thanks to Hubble Space Telescope high-precision proper motions.

Unexpected kinematics of multiple populations in globular clusters

AbstractPresent-day structural and kinematical properties of multiple populations (MPs) can provide useful information about the physical mechanisms driving the formation and early evolution of globular clusters (GCs). As part of a large project aimed at characterizing the kinematics of MPs, here we present a detailed multi-epoch analysis of the low-mass GC NGC6362. We find that MPs in this system show significant differences in their line-of-sight velocity dispersion profiles. This result is totally unexpected given the dynamical age and fraction of mass lost by NGC6362. We also find that the binary fraction is remarkably larger in the first (FP) than in the second population (SP). We show that such a difference can efficiently inflate the velocity dispersion of FP at intermediate/large cluster-centric distances with respect to SP. Indeed, our results nicely match the predictions of state-of-the art N-body simulations of the co-evolution of MPs in GCs including the effect of binaries.

Dynamical effects on the stellar mass function of multiple stellar populations in globular clusters

AbstractWe present a brief summary of the results of a study of the effects of dynamical evolution on the stellar mass function of multiple-population globular clusters. Theoretical studies have predicted that the process of multiple-population cluster formation results in a system in which second-generation (2G) stars are initially more centrally concentrated than first-generation (1G) stars. In the study presented here, we have explored the implications of the initial differences between the 2G and 1G structural properties for the evolution of the local (measured at different distances from a cluster center) and global mass function. We have studied both systems in which 1G and 2G stars start with the same initial mass function (IMF) and systems in which 1G and 2G stars have different IMFs. Finally we have explored the evolution of the spatial mixing and found that the multiscale nature of the clusters studied leads to a dependence of the mixing rate on the stellar mass.

Kinematical evolution of Globular Clusters

AbstractWe present several results of the study of the evolution of globular clusters’ internal kinematics, as driven by two-body relaxation and the interplay between internal angular momentum and the external Galactic tidal field. Via a large suite of N-body simulations, we explored the three-dimensional velocity space of tidally perturbed clusters, by characterizing their degree of velocity dispersion anisotropy and their rotational properties. These studies have shown that a cluster’s kinematical properties contain distinct imprints of the cluster’s initial structural properties, dynamical history, and tidal environment. Building on this fundamental understanding, we then studied the dynamics of multiple stellar populations in globular clusters, with attention to the largely unexplored role of angular momentum.

Mass loss of different stellar populations in globular clusters

AbstractOnce the age and metallicity are fixed, the colour distribution of horizontal branch stars in a globular cluster depends on few parameters: the helium abundance of the population and the mass lost during the pre-HB stages. These parameters are usually derived from the HB itself, hence they are degenerate. Breaking this degeneracy and understanding their role is a tricky and challenging problem that no study has solved yet. Combining the information obtained from the chromosome maps and the analysis of multi-band photometry with state of the art stellar evolution models, we can obtain a solid estimate of Y for the various stellar populations in a GC. We will then have, for the first time, the possibility to break the parameters’ degeneracy on the HB, understand the role of the mass loss, and lay the foundation to build another piece of the multiple populations mosaic.

The Role of “black hole burning” in the evolution of dense star clusters

AbstractAs self-gravitating systems, dense star clusters exhibit a natural diffusion of energy from their innermost to outermost regions, leading to a slow and steady contraction of the core until it ultimately collapses under gravity. However, in spite of the natural tendency toward “core collapse,” the globular clusters (GCs) in the Milky Way exhibit a well-observed bimodal distribution in core radii separating the core-collapsed and non-core-collapsed clusters. This suggests an internal energy source is at work, delaying the onset of core collapse in many clusters. Over the past decade, a large amount of work has suggested that stellar black holes (BHs) play a dynamically-significant role in clusters throughout their entire lifetimes. Here we review our latest understanding of BH populations in GCs and demonstrate that, through their dynamical interaction with their host cluster, BHs can naturally explain the distinction between core-collapsed and non-core-collapsed clusters through a process we call “black hole burning.”