Massive Globular Clusters Research Articles

The giants that were born swiftly – implications of the top-heavy stellar initial mass function on the birth conditions of globular clusters

ABSTRACT Recent results suggest that the initial mass function (IMF) of globular clusters (GCs) is metallicity and density dependent. Here, it is studied how this variation affects the initial masses and the numbers of core-collapse supernovae (CCSNe) required to reproduce the observed iron spreads in GCs. The IMFs of all of the investigated GCs were top-heavy implying larger initial masses compared to previous results computed assuming an invariant canonical IMF. This leads to more CCSNe being required to explain the observed iron abundance spreads. The results imply that the more massive GCs formed at smaller Galactocentric radii, possibly suggesting in situ formation of the population II halo. The time until star formation ended within a proto-GC is computed to be 3.5–4 Myr, being slightly shorter than the 4 Myr obtained using the canonical IMF. Therefore, the impact of the IMF on the time for which IMF lasts is small.

J-PLUS: Detecting and studying extragalactic globular clusters

Context. Extragalactic globular clusters (GCs) are key objects in studies of galactic histories. The advent of wide-field surveys, such as the Javalambre Photometric Local Universe Survey (J-PLUS), offers new possibilities for the study of these systems. Aims. We performed the first study of GCs in J-PLUS to recover information on the history of NGC 1023, taking advantage of wide-field images and 12 filters. Methods. We developed the semiautomatic pipeline GCFinder for detecting GC candidates in J-PLUS images, which can also be adapted to similar surveys. We studied the stellar population properties of a sub-sample of GC candidates using spectral energy distribution (SED) fitting. Results. We found 523 GC candidates in NGC 1023, about 300 of which are new. We identified subpopulations of GC candidates, where age and metallicity distributions have multiple peaks. By comparing our results with the simulations, we report a possible broad age-metallicity relation, supporting the notion that NGC 1023 has experienced accretion events in the past. With a dominating age peak at 1010 yr, we report a correlation between masses and ages that suggests that massive GC candidates are more likely to survive the turbulent history of the host galaxy. Modeling the light of NGC 1023, we find two spiral-like arms and detect a displacement of the galaxy’s photometric center with respect to the outer isophotes and center of GC distribution (~700pc and ~1600pc, respectively), which could be the result of ongoing interactions between NGC 1023 and NGC 1023A. Conclusions. By studying the GC system of NGC 1023 with J-PLUS, we showcase the power of multi-band surveys for these kinds of studies and we find evidence to support the complex accretion history of the host galaxy.

Detection of a 105 M⦿ Black Hole in M31's Most Massive Globular Cluster

Detection of a 105 M⦿ Black Hole in M31's Most Massive Globular Cluster

EROSITA study of the globular cluster 47 Tucanae

Aims. We present the results of the analysis of five observations of the globular clutser 47 Tucanae (47 Tuc) with the extended Roentgen Survey with an Imaging Telescope Array (eROSITA) on board the Spektrum-Roentgen-Gamma (Spektr-RG, SRG). We study the X-ray population in the field of one of the most massive globular clusters in our Milky Way. We focused on the classification of point-like sources in the field of 47 Tuc. The unresolved dense core of 47 Tuc (1.7 radius) and also sources that show extended emission are excluded from this study. Methods. We applied different methods of X-ray spectral and timing analysis together with multi-wavelength studies to classify the X-rays sources in the field of 47 Tuc. Results. We detected 888 point-like sources in the energy range of 0.2–5.0 keV. We identified 126 background active galactic nuclei and 25 foreground stars. One of the foreground stars is classified as a variable M dwarf. We also classified 14 X-ray sources as members of 47 Tuc, including one symbiotic star, two quiescent low-mass X-ray binaries, and four cataclysmic variables. There are also five X-ray sources that can either be a cataclysmic variable or a contact binary, and also one X-ray source can be an active binary (type RS CVn). We identified one X-ray binary that belongs to the Small Magellanic Cloud. Moreover, we calculated the X-ray luminosity function of 47 Tuc. No significant population that seems to belong to the globular cluster is observed in the energy range of 0.5–2.0 keV using eROSITA observations.

Testing feedback from star clusters in simulations of the Milky Way formation

ABSTRACT We present a suite of galaxy formation simulations that directly model star cluster formation and disruption. Starting from a model previously developed by our group, here we introduce several improvements to the prescriptions for cluster formation and feedback, then test these updates using a large suite of cosmological simulations of Milky Way mass galaxies. We perform a differential analysis with the goal of understanding how each of the updates affects star cluster populations. Two key parameters are the momentum boost of supernova feedback fboost and star formation efficiency per free-fall time ϵff. We find that fboost has a strong influence on the galactic star formation rate, with higher values leading to less star formation. The efficiency ϵff does not have a significant impact on the global star formation rate, but dramatically changes cluster properties, with increasing ϵff leading to a higher maximum cluster mass, shorter age spread of stars within clusters, and higher integrated star formation efficiencies. We also explore the redshift evolution of the observable cluster mass function, finding that most massive clusters have formed at high redshift z > 4. Extrapolation of cluster disruption to z = 0 produces good agreement with both the Galactic globular cluster mass function and age–metallicity relation. Our results emphasize the importance of using small-scale properties of galaxies to calibrate subgrid models of star cluster formation and feedback.

Revisiting the relation between the number of globular clusters and galaxy mass for low-mass galaxies

ABSTRACT Using a new method to estimate total galaxy mass (MT) and two samples of low-luminosity galaxies containing measurements of the number of globular clusters (GCs) per galaxy (NGC), we revisit the NGC–MT relation using a total of 203 galaxies, 157 of which have MT ≤ 1010 M⊙. We find that the relation is nearly linear, NGC ∝ MT0.92 ± 0.08 down to at least MT ∼ 108.75 M⊙. Because the relationship extends to galaxies that average less than one GC per galaxy and to a mass range in which mergers are relatively rare, the relationship cannot be solely an emergent property of hierarchical galaxy formation. The character of the radial GC distribution in low-mass galaxies, and the lack of mergers at these galaxy masses, also appears to challenge models in which the GCs form in central, dissipatively concentrated high-density, high-pressure regions and are then scattered to large radius. The slight difference between the fitted power-law exponent and a value of one leaves room for a shallow MT-dependent variation in the mean mass per GC that would allow the relation between total mass in GCs and MT to be linear.

A Population of Luminous Globular Clusters and Stripped Nuclei with Elevated Mass to Light Ratios around NGC 5128* * This paper includes data gathered with the 6.5 m Magellan Telescope at Las Campanas Observatory, Chile.

The dense central regions of tidally disrupted galaxies can survive as ultracompact dwarfs (UCDs) that hide among the luminous globular clusters (GCs) in the halo of massive galaxies. An exciting confirmation of this model is the detection of overmassive black holes in the centers of some UCDs, which also lead to elevated dynamical mass-to-light ratios (M/L dyn). Here we present new high-resolution spectroscopic observations of 321 luminous GC candidates in the massive galaxy NGC 5128/Centaurus A. Using these data we confirm 27 new luminous GCs, and measure velocity dispersions for 57 luminous GCs (with g-band luminosities between 2.5 × 105 and 2.5 × 107 L ⊙), of which 48 are new measurements. Combining these data with size measurements from Gaia, we determine the M/L dyn for all 57 luminous GCs. We see a clear bimodality in the M/L dyn distribution, with a population of normal GCs with mean M/L dyn = 1.51 ± 0.31, and a second population of ∼20 GCs with elevated mean M/L dyn = 2.68 ± 0.22. We show that black holes with masses ∼4%–18% of the luminous GCs can explain the elevated mass-to-light ratios. Hence, it is plausible that the NGC 5128 sources with elevated M/L dyn are mostly stripped galaxy nuclei that contain massive central black holes, though future high spatial resolution observations are necessary to confirm this hypothesis for individual sources. We also present a detailed discussion of an extreme outlier, VHH81-01, one of the largest and most massive GC in NGC 5128, making it an exceptionally strong candidate to be a tidally stripped nucleus.

A possible solution to the Milky Way’s binary-deficient retrograde stellar population

Context. The fraction of field binaries on retrograde orbits about the Milky Way is significantly lower compared to its prograde counterpart. Chemical and dynamical evidence suggests that the retrograde stellar population originates from ω Centauri, which is either the most massive globular cluster (GC) of the Milky Way or the putative core of a former dwarf galaxy. Aims. Star formation conditions required to produce the retrograde binary population are constrained assuming that the retrograde stellar population originates from ω Centauri’s progenitor. Methods. We match the observed low binary fraction with dynamical population synthesis models, including a universal initial binary population and dynamical processing in star clusters, making use of the publicly available binary population synthesis tool BIPOS1. Results. It is found that either the GC progenitor of ω Cen must have formed with a stellar density of ≈108 M⊙ pc−3 or the ω Centauri dwarf galaxy’s progenitor star cluster population must have formed in an extreme starburst with a star formation rate exceeding 1000 M⊙ yr−1 and probably a top-heavy embedded-cluster mass function with suppressed low-mass cluster formation. The separation and mass-ratio distribution for retrograde field binaries are predicted for comparison with future observations. Conclusions. A viable solution for the deficiency of binaries on retrograde orbits is presented, and star formation conditions for ω Centauri as well as orbital parameter distributions for the Milky Way’s retrograde binary population are predicted. The dwarf galaxy origin for ω Centauri is tentatively preferred within the present context.

Survey of Multiple Populations in Globular Clusters among Very-low-mass Stars

Abstract Recent work has shown that near-infrared (NIR) Hubble Space Telescope (HST) photometry allows us to disentangle multiple populations (MPs) among M dwarfs of globular clusters (GCs) and to investigate this phenomenon in very-low-mass (VLM) stars. Here, we present the color–magnitude diagrams of nine GCs and the open cluster NGC 6791 in the F110W and F160W bands of HST, showing that the main sequences (MSs) below the knee are either broadened or split, thus providing evidence of MPs among VLM stars. In contrast, the MS of NGC 6791 is consistent with a single population. The color distribution of M dwarfs dramatically changes between different GCs, and the color width correlates with the cluster mass. We conclude that the MP ubiquity, variety, and dependence on GC mass are properties common to VLM and more-massive stars. We combined UV, optical, and NIR observations of NGC 2808 and NGC 6121 (M4) to identify MPs along with a wide range of stellar masses (∼0.2–0.8  ⊙ ), from the MS turnoff to the VLM regime, and measured, for the first time, their mass functions (MFs). We find that the fraction of MPs does not depend on the stellar mass and that their MFs have similar slopes. These findings indicate that the properties of MPs do not depend on stellar mass. In a scenario where the second generations formed in higher-density environments than the first generations, the possibility that the MPs formed with the same initial MF would suggest that it does not depend on the environment.

The Mg–K Anticorrelation in ω Centauri

Abstract We present [K/Fe] abundance ratios for a sample of 450 stars in ω Centauri, using high-resolution spectra acquired with the multiobject spectrograph FLAMES@VLT. Abundances for Fe, Na, and Mg were also derived. We detected intrinsic K variations in the analyzed stars. Moreover, [K/Fe] shows a significant correlation with [Na/Fe] and an anticorrelation with [Mg/Fe]. The presence of a clear-cut Mg–K anticorrelation makes ω Centauri the third stellar system, after NGC 2419 and NGC 2808, hosting a subpopulation of stars with [Mg/Fe] < 0.0 dex, K-enriched in the case of ω Centauri by ∼0.3 dex with respect to Mg-rich stars ([Mg/Fe] > 0.0 dex). The correlation/anticorrelation between K and other light elements involved in chemical anomalies supports the idea that the spread in [K/Fe] can be associated with the same self-enrichment process typical of globular clusters. We suggest that significant variations in K abundances perhaps can be found in the most massive and/or metal-poor globular clusters as a manifestation of an extreme self-enrichment process. Theoretical models face problems explaining K production in globular clusters. Indeed, models where asymptotic giant branch stars are responsible for the Mg–K anticorrelation only qualitatively agree with the observations. Finally, we discovered a peculiar star with an extraordinary K overabundance ([K/Fe] = +1.60 dex) with respect to the other stars with similar [Mg/Fe]. We suggest that this K-rich star could be formed from the pure ejecta of AGB stars before dilution with pristine material.

Metallicity of the globular cluster NGC 6388 based on high-resolution spectra of more than 160 giant stars

NGC 6388 is one of the most massive Galactic globular clusters (GC) and it is an old, metal-rich Galactic bulge cluster. By exploiting previous spectroscopic observations, we were able to bypass the uncertainties in membership related to the contamination from strong field stars. We present the abundance analysis of 12 new giant stars with UVES spectra and 150 giants with GIRAFFE spectra acquired at the ESO-VLT. We derived radial velocities, atmospheric parameters, and iron abundances for all the stars. When combined with the previous data, we obtained a grand total of 185 stars homogeneously analysed in NGC 6388 from high-resolution spectroscopy. The average radial velocity of the 185 stars is 81.2 ± 0.7, rms 9.4 km s−1. We obtained an average metallicity [Fe/H] = −0.480 dex, rms = 0.045 dex (35 stars), and [Fe/H] = −0.488 dex, rms = 0.040 dex (150 stars) from the UVES and GIRAFFE samples, respectively. Comparing these values to the internal errors in abundance, we excluded the presence of a significant intrinsic metallicity spread within the cluster. Since about a third of giants in NGC 6388 is claimed to belong to the ‘anomalous red giants’ in the HST pseudo-colour map defining the so-called type-II GCs, we conclude that either enhanced metallicity is not a necessary requisite to explain this classification (as also suggested by the null iron spread for NGC 362) or NGC 6388 is not a type-II globular cluster.

The spatial distribution of globular clusters in dwarf spheroidal galaxies and the timing problem

ABSTRACT The dynamical friction time-scale of massive globular clusters (GCs) in the inner regions of cuspy dark haloes in dwarf spheroidal (dSph) galaxies can be much shorter than the Hubble time. This implies that a small fraction of the GCs is expected to be caught close to the centre of these galaxies. We compare the radial distribution of GCs predicted in simple Monte Carlo models with that of a sample of 38 spectroscopically confirmed GCs plus 17 GC candidates associated mainly with low-luminosity dSph galaxies. If dark matter haloes follow an NFW profile, the observed number of off-centre GCs at projected distances less than one-half the galaxy effective radius is significantly higher than the models predict. This timing problem can be viewed as a fine-tuning of the starting GC distances. As a result of the short sinking time-scale for GCs in the central regions, the radial distribution of GCs is expected to evolve significantly during the next 1–2 Gyr. However, dark matter haloes with cores of size comparable to the galaxy effective radii can lead to a slow orbital inspiral of GCs in the central regions of these galaxies, providing a simple solution to the timing problem. We also examine any indication of mass segregation in the summed distribution of our sample of GCs.

The physics governing the upper truncation mass of the globular cluster mass function

ABSTRACT The mass function of globular cluster (GC) populations is a fundamental observable that encodes the physical conditions under which these massive stellar clusters formed and evolved. The high-mass end of star cluster mass functions are commonly described using a Schechter function, with an exponential truncation mass Mc, *. For the GC mass functions in the Virgo galaxy cluster, this truncation mass increases with galaxy mass (M*). In this paper, we fit Schechter mass functions to the GCs in the most massive galaxy group ($M_{\mathrm{200}} = 5.14 \times 10^{13} \, {\rm M}_{\odot }$) in the E-MOSAICS simulations. The fiducial cluster formation model in E-MOSAICS reproduces the observed trend of Mc, * with M* for the Virgo cluster. We therefore examine the origin of the relation by fitting Mc, * as a function of galaxy mass, with and without accounting for mass loss by two-body relaxation, tidal shocks and/or dynamical friction. In the absence of these mass-loss mechanisms, the Mc, *-M* relation is flat above $M_* \gt 10^{10}\, {\rm M}_{\odot }$. It is therefore the disruption of high-mass GCs in galaxies with $M_{*}\sim 10^{10} \, {\rm M}_{\odot }$ that lowers the Mc, * in these galaxies. High-mass GCs are able to survive in more massive galaxies, since there are more mergers to facilitate their redistribution to less-dense environments. The Mc, * − M* relation is therefore a consequence of both the formation conditions of massive star clusters and their environmentally dependent disruption mechanisms.

Detection of a 100,000 M ⊙ black hole in M31's Most Massive Globular Cluster: A Tidally Stripped Nucleus

Abstract We investigate the presence of a central black hole (BH) in B023-G078, M31's most massive globular cluster. We present high-resolution, adaptive-optics assisted, integral-field spectroscopic kinematics from Gemini/NIFS that show a strong rotation (∼20 km s−1) and a velocity dispersion rise toward the center (37 km s−1). We combine the kinematic data with a mass model based on a two-component fit to HST ACS/HRC data of the cluster to estimate the mass of a putative BH. Our dynamical modeling suggests a >3σ detection of a BH component of 9.1 − 2.8 + 2.6 × 10 4 M ⊙ (1σ uncertainties). The inferred stellar mass of the cluster is 6.22 − 0.05 + 0.03 × 10 6 M ⊙ , consistent with previous estimates, thus the BH makes up 1.5% of its mass. We examine whether the observed kinematics are caused by a collection of stellar mass BHs by modeling an extended dark mass as a Plummer profile. The upper limit on the size scale of the extended mass is 0.56 pc (95% confidence), which does not rule out an extended mass. There is compelling evidence that B023-G078 is the tidally stripped nucleus of a galaxy with a stellar mass >109 M ⊙, including its high-mass, two-component luminosity profile, color, metallicity gradient, and spread in metallicity. Given the emerging evidence that the central BH occupation fraction of >109 M ⊙ galaxies is high, the most plausible interpretation of the kinematic data is that B023-G078 hosts a central BH. This makes it the strongest BH detection in a lower-mass (<107 M ⊙) stripped nucleus, and one of the few dynamically detected intermediate-mass BHs.

Rebounding Cores to Build Star Cluster Multiple Populations

Abstract We present a novel approach to the riddle of star cluster multiple populations. Stars form from molecular cores. But not all cores form stars. Following their initial compression, such “failed” cores re-expand, rather than collapsing. We propose that their formation and subsequent dispersal regulate the gas density of cluster-forming clumps and, therefore, their core and star formation rates. Clumps for which failed cores are the dominant core type experience star formation histories with peaks and troughs (i.e., discrete star formation episodes). In contrast, too few failed cores results in smoothly decreasing star formation rates. We identify three main parameters shaping the star formation history of a clump: the star and core formation efficiencies per free-fall time, and the timescale on which failed cores return to the clump gas. The clump mass acts as a scaling factor. We use our model to constrain the density and mass of the Orion Nebula Cluster progenitor clump, and to caution that the star formation histories of starburst clusters may contain close-by peaks concealed by stellar age uncertainties. Our model generates a great variety of star formation histories. Intriguingly, the chromosome maps and O–Na anticorrelations of old globular clusters also present diverse morphologies. This prompts us to discuss our model in the context of globular cluster multiple stellar populations. More massive globular clusters exhibit stronger multiple stellar population patterns, which our model can explain if the formation of the polluting stars requires a given stellar mass threshold.

Spectroscopy and Photometry of the Least Massive Type II Globular Clusters: NGC 1261 and NGC 6934*

Abstract Recent work has revealed two classes of globular clusters (GCs), dubbed Type I and Type II. Type II GCs are characterized by both a blue and a red red giant branch composed of stars with different metallicities, often coupled with distinct abundances in the slow neutron-capture elements (s-elements). Here we continue the chemical tagging of Type II GCs by adding the two least massive clusters of this class, NGC 1261 and NGC 6934. Based on both spectroscopy and photometry, we find red stars in NGC 1261 to be slightly enhanced in [Fe/H] by ∼0.1 dex and confirm that red stars of NGC 6934 are enhanced in iron by ∼0.2 dex. Neither NGC 1261 nor NGC 6934 show internal variations in the s-elements, which suggests a GC mass threshold for the occurrence of s-process enrichment. We found a significant correlation between the additional Fe locked in the red stars of Type II GCs and the present-day mass of the cluster. Nevertheless, most Type II GCs retained a small fraction of Fe produced by SNe II, lower than the 2%; NGC 6273, M54, and ω Centauri are remarkable exceptions. In the Appendix, we infer for the first time chemical abundances of lanthanum, assumed as representative of the s-elements, in M54, the GC located in the nucleus of the Sagittarius dwarf galaxy. Red-sequence stars are marginally enhanced in [La/Fe] by 0.10 ± 0.06 dex, in contrast with the large [La/Fe] spread of most Type II GCs. We suggest that different processes are responsible for the enrichment in iron and s-elements in Type II GCs.

The impact of primordial binary on the dynamical evolution of intermediate massive star clusters

ABSTRACT Observations found that star clusters contain a large fraction of binaries. Tight binaries are an important heating source that influences the long-term dynamical evolution of star clusters. However, due to the limitation of N-body tool, previous theoretical modelling for globular clusters (GCs) by using direct N-body simulations has not investigated how a large fraction of primordial binaries affect their long-term evolution. In this work, by using the high-performance N-body code, petar, we carry out star-by-star models for intermediate massive GCs (N = 100 000) with the primordial binary fraction varying from 0 to 1. We find that when a stellar-mass black hole (BH) subsystem exists, the structural evolution of GCs (core and half-mass radii) only depends on the properties of massive primordial binaries, because they affect the number of BH binaries (BBHs), which dominate the binary heating process. Low-mass binaries including double white dwarf binaries (BWDs) have almost no influence on the dynamics. Meanwhile, only gravitational wave mergers from BBHs are strongly affected by dynamical interactions, while low-mass mergers from BWDs show no difference in the isolated environment (field) and in GCs. Low-mass binaries become important only after most BHs escape and the core collapse of light stars occurs. Our result suggests that for N-body modelling of GCs with a BH subsystem dominating binary heating, it is not necessary to include low-mass binaries. These binaries can be studied separately by using standalone binary stellar evolution codes. This way can significantly reduce the computing cost.

Correction to: UV photometry of spotted stars in the horizontal branch of the globular cluster NGC 2808 using AstroSat

A recent study of hot (20,000 to 30,000 K) extreme horizontal branch (EHB) stars in globular clusters (GCs) has led to the discovery of their variability. It is suggested that this variability is driven by the projected rotation of magnetic spots on the stellar surfaces and is expected to have higher amplitudes at shorter wavelengths. Here, we present the analysis of such hot stars in the massive GC NGC 2808 using the Ultraviolet Imaging Telescope (UVIT), aboard AstroSat. We use the UVIT data in combination with the Hubble Space Telescope UV Globular Cluster Survey (HUGS) data for the central region (within $\sim$ $2.7'$ $\times$ $2.7'$) and ground-based optical photometry for the outer parts of the cluster. We generate the far-UV (FUV) - optical colour-magnitude diagrams (CMDs) and in these we find a population of hot EHB stars fainter than the zero-age horizontal branch (ZAHB) model. A comparison of our FUV magnitudes of the already reported variable EHB stars (vEHBs) shows that the longest period vEHBs are the faintest, along with a tentative correlation between rotation period and UV magnitude of spotted stars. In order to firmly establish any correlation, further study is essential.

Evolution of globular-cluster systems of ultra-diffuse galaxies due to dynamical friction in MOND gravity

Context. Dynamical friction can be used to distinguish Newtonian gravity and modified Newtonian dynamics (MOND) because it works differently in these frameworks. This concept, however, has yet to be explored very much with MOND. Previous simulations showed weaker dynamical friction during major mergers for MOND than for Newtonian gravity with dark matter. Analytic arguments suggest the opposite for minor mergers. In this work, we verify the analytic predictions for MOND by high-resolution N-body simulations of globular clusters (GCs) moving in isolated ultra-diffuse galaxies (UDGs). Aims. We test the MOND analog of the Chandrasekhar formula for the dynamical friction proposed by Sánchez-Salcedo on a single GC. We also explore whether MOND allows GC systems of isolated UDGs to survive without sinking into nuclear star clusters. Methods. The simulations are run using the adaptive-mesh-refinement code Phantom of Ramses. The mass resolution is 20 M⊙ and the spatial resolution 50 pc. The GCs are modeled as point masses. Results. Simulations including a single GC reveal that, as long as the apocenter of the GC is over about 0.5 effective radii, the Sánchez-Salcedo formula works excellently, with an effective Coulomb logarithm increasing with orbital circularity. Once the GC reaches the central kiloparsec, its sinking virtually stops, likely because of the core stalling mechanism. In simulations with multiple GCs, many of them sink toward the center, but the core stalling effect seems to prevent them from forming a nuclear star cluster. The GC system ends up with a lower velocity dispersion than the stars of the galaxy. By scaling the simulations, we extend these results to most UDG parameters, as long as these UDGs are not external-field dominated. We verify analytically that approximating the GCs by point masses has little effect if the GCs have the usual properties, but for massive GCs such as those observed in the NGC 1052-DF2 galaxy, further simulations with resolved GCs are desirable.

CAPOS: The bulge Cluster APOgee Survey

We present results from a study of 15 red giant members of the intermediate-metallicity globular cluster (GC) FSR 1758 using high-resolution, near-infrared spectra collected with the Apache Point Observatory Galactic Evolution Experiment II survey (APOGEE-2) that were obtained as part of CAPOS (the bulge Cluster APOgee Survey). Since its very recent discovery as a massive GC in the bulge region, evoking the name Sequoia, this has been an intriguing object with a highly debated origin, and initially led to the suggestion of a purported progenitor dwarf galaxy of the same name. In this work, we use new spectroscopic and astrometric data to provide additional clues as to the nature of FSR 1758. Our study confirms the GC nature of FSR 1758, and as such we report the existence of the characteristic N-C anticorrelation and Al-N correlation for the first time. We thereby reveal the existence of the multiple-population phenomenon, similar to that observed in virtually all GCs. Furthermore, the presence of a population with strongly enriched aluminum makes it unlikely that FSR 1758 is the remnant nucleus of a dwarf galaxy because Al-enhanced stars are uncommon in dwarf galaxies. We find that FSR 1758 is slightly more metal rich than previously reported in the literature; this source has a mean metallicity [Fe/H] between −1.43 to −1.36, depending on the adopted atmospheric parameters and a scatter within observational error, again pointing to its GC nature. Overall, the α-enrichment (≳ + 0.3 dex), Fe-peak (Fe, Ni), light (C, N), and odd-Z (Al) elements follow the trend of intermediate-metallicity GCs. Isochrone fitting in the Gaia bands yields an estimated age of ∼11.6 Gyr. We used the exquisite kinematic data, including our CAPOS radial velocities and Gaia eDR3 proper motions, to constrain the N-body density profile of FSR 1758, and found that it is as massive (∼2.9 ± 0.6 × 105 M⊙) as NGC 6752. We confirm a retrograde and eccentric orbit for FSR 1758. A new examination of its dynamical properties with the GravPot16 model favors an association with the Gaia-Enceladus-Sausage accretion event. Thus, paradoxically, the cluster that gave rise to the name of the Sequoia dwarf galaxy does not appear to belong to this specific merging event.