VVV CL001: Likely the Most Metal-Poor Surviving Globular Cluster in the Inner Galaxy

We present the first high-resolution abundance analysis of the globular cluster VVV~CL001, which resides in a region dominated by high interstellar reddening towards the Galactic Bulge. Using \textit{H}-band spectra acquired by the Apache Point Observatory Galactic Evolution Experiment (APOGEE), we identified two potential members of the cluster, and estimate from their Fe I lines that the cluster has an average metallicity of [Fe/H] = $-2.45$ with an uncertainty due to systematics of 0.24 dex. We find that the light-(N), $\alpha$-(O, Mg, Si), and Odd-Z (Al) elemental abundances of the stars in VVV~CL001 follow the same trend as other Galactic metal-poor globular clusters. This makes VVV~CL001 possibly the most metal-poor globular cluster identified so far within the Sun's galactocentric distance and likely one of the most metal-deficient clusters in the Galaxy after ESO280-SC06. Applying statistical isochrone fitting, we derive self-consistent age, distance, and reddening values, yielding an estimated age of $11.9^{+3.12}_{-4.05}$ Gyr at a distance of $8.22^{+1.84}_{-1.93}$ kpc, revealing that VVV~CL001 is also an old GC in the inner Galaxy. The Galactic orbit of VVV~CL001 indicates that this cluster lies on a halo-like orbit that appears to be highly eccentric. Both chemistry and dynamics support the hypothesis that VVV~CL001 could be an ancient fossil relic left behind by a massive merger event during the early evolution of the Galaxy, likely associated with either the Sequoia or the \textit{Gaia}-Enceladus-Sausage structures.


INTRODUCTION
The inner 1 Milky Way (MW) contain a significant population of globular clusters (GCs) with a wide range Corresponding author: José G. Fernández-Trincado jose.fernandez@uda.cl of metallicities (−2.37 < [Fe/H] 0; Harris 1996Harris , 2010, most of which are still poorly explored due to the large foreground extinction and high field-star densities that complicate the analysis, especially at low Galactic latitudes. These limitations have been mitigated with wide field near-infrared imaging such as the Vista Variables in the Via Lactea survey (VVV), and its extension the VVVX survey (Minniti et al. 2010;Smith et al. 2018), which has expanded the family of Galactic GCs to more than 300 candidates by the inclusion of objects in the inner Galaxy (see e.g. Minniti et al. 2017b,a;Palma et al. 2019;Garro et al. 2020), including those in the bulge (Geisler et al. in prep.).
In this context, the high-resolution (R > 22, 0000) capabilities of the near-IR multi-fiber spectrographs of the Apache Point Observatory Galactic Evolution Experiment (APOGEE-2; Majewski et al. 2017) allow measurement of new parameters (radial velocity, metallicity, and detailed chemical abundances for many species) with high precision for a large number of Galactic GCs in a homogeneous way (see e.g. Mészáros et al. 2015;Schiavon et al. 2017;Masseron et al. 2019;Fernández-Trincado et al. 2019dMészáros et al. 2020).
Very metal-poor GCs are preferentially associated with the oldest components of galaxies, and are often used as cosmic clocks to track the enrichment history of their host galaxy (Geisler et al. 1995). Thus, beyond the intrinsic scientific value of identifying such ancient objects, the measurements of their chemical composition can therefore provide insight into the build-up of the chemical elements in the earliest epoch of the Milky Way.
To date, the most metal-poor MW GC known is ESO280-SC06 (located ∼15 kpc from the Galactic center), with an estimated metallicity of [Fe/H]= −2.48 +0.06 −0.11 (Simpson 2018), and associated to the Gaia-Enceladus-Sausage (GES; Massari et al. 2019). Here, we report another possible extreme case, VVV CL001, which becomes possibly the most metal-poor GC known inside the Sun's Galactocentric distance, and likely the most metal-poor GC in the Galaxy. It has survived the strong tidal field of the inner MW during its interaction, and its existence implies that GCs below that of the empirical metallicity floor (e.g. Wan et al. 2020;Larsen et al. 2020) have been formed, but only a few exceptional cases have survived during Galactic evolution.
In this Letter, we make use of APOGEE-2 spectra to provide the first spectroscopic study of VVV CL001.

OBSERVATIONAL DATA
We use an interim release data product of the APOGEE-2 survey (Majewski et al. 2017) part of the Sloan Digital Sky Survey IV (SDSS-IV; Blanton et al. 2017) that includes data taken after the release of DR16 (Ahumada et al. 2020), to investigate, for the first time, the chemical composition of VVV CL001.
The (300-fiber) APOGEE instruments are highresolution (R ∼ 22, 500), near-infrared (NIR) spectrographs (Wilson et al. 2019) observing all the components of the MW (halo, disc, and bulge/bar) from the Northern Hemisphere on the 2.5m telescope at Apache Point Observatory (APO, APOGEE-2N; Gunn et al. 2006) and from the Southern Hemisphere on the Irénée du Pont 2.5m telescope (Bowen & Vaughan 1973) at Las Campanas Observatory (LCO, APOGEE-2S). Each instrument records most of the H -band (1.51µm -1.69µm). See Nidever et al. (2015), García Pérez et al. (2016), and Holtzman et al. (2018 for details regarding the data reduction process and standard stellar parameter estimates. As of February 2020, the dual APOGEE-2 instruments have observed some ∼600,000 sources across the MW, targeting these stars according to methods described in Zasowski et al. (2013) and Zasowski et al. (2017), with updates to the targeting plan described in Beaton et al. and Santana et al. (in prep.).
The present study focuses on VVV CL001, a GC discovered by the VVV survey (see e.g., Minniti et al. 2011). This GC lies in the direction of the Galactic bulge, and is strongly dominated by large foreground extinction, E(B−V) 2.0 mag , hampering the observations of this object in the optical bands.
VVV CL001 stars fall on the same APOGEE-2 plugplates as those associated with the very nearby GC UKS 1, as shown in panel (a) of Figure 1 (see also Fernández-Trincado et al. 2020e, for details regarding to UKS 1), however the high radial velocities (RVs) of potential VVV CL001 members allow us to cleanly distinguish the UKS 1 sources from VVV CL001 stars. Thus, two high-confidence VVV CL001 members, based on proper motions (see panel (b) of Figure 1) from Gaia Early Data Release 3 (Gaia EDR3: Gaia Collaboration et al. 2020), location in the color-magnitude diagrams (CMDs), and RVs from APOGEE-2 were identified 0.1 from the cluster centre. Panels (c) and (d) of Figure 1 show that both our derived [Fe/H] and RV of the VVV CL001 stars are clearly distinct as compared to the foreground and background stars (hereafter field stars). The [Fe/H] and RVs are at least 1.5 dex and ∼ 150 km s −1 offset from the field stars, respectively. Both of these offsets are very large and imply that VVV CL001 is a truly extreme GC. The two potential cluster members from APOGEE-2 are red giant branch (RGB) stars close to the tip of the giant branch as shown in panel (e) of Figure 1. In order to have a self-consistent method for age derivation via statistical isochrone fitting, we use the SIRIUS code (Souza et al. 2020), and the most probable cluster members in the VVV catalogue located inside 1.5 from the cluster center and that have proper motions compatible with that of VVV CL001, as well as those sources with RVs information (see Section 4), which have been marked as triangle (for Baumgardt's data set) and square (for APOGEE-2 data set) open symbols in panel (e) of Figure 1. Due to the quality of our data, we applied some assumptions to obtain an age distribution: a uniform prior in age between 1 and 15 Gyr combined with a slow drop above the age of the universe (13.7 Gyr; Planck Collaboration et al. 2016); the metallicity was variated around the value determined with high-resolution spectroscopy in the present work; the isochrone is limited to log g < 4.5 representing the RGB region. We dereddened and extinction-corrected the VVV+2MASS Ks s and J − K s colors with the bulge-specific reddening maps from Gonzalez et al. (2011Gonzalez et al. ( , 2012 assuming the reddening law of Cardelli et al. (1989). We also noticed that in a ∼2 × 2 area the differential reddening across the field do not affect the CMD of VVV CL001, with a negligible variation of 0.03 mag in K s . Finally, we adopted the Dart-  provides a reasonable solution both in the over-plotted isochrone (panel e) and the posterior distributions of the corner plot (panel f). To represent the distributions, we adopt the median as the most probable value and the uncertainties calculated from the 16 th and 84 th percentiles. We found an age of 11.9 +3.12 −4.05 Gyr and a probable distance of ∼ 8.22 +1.84 −1.93 kpc. Also, we want to stress that without the adopted assumptions, the internal error of the age determination could increase and give not a clear age distribution. Consequently, our probable solutions, within 1 − σ, fit well the central part of the CMD, providing confidence that the age estimate is a reasonable determination for VVV CL001.

ELEMENTAL ABUNDANCES
As the ASPCAP/APOGEE-2 pipeline (García Pérez et al. 2016) does not provide [Fe/H] and [X/Fe] determinations for VVV CL001 stars, we followed the same technique as described in Fernández-Trincado et al. (2019b,a,c,d, 2020a, and carried out a consistent chemical-abundance analysis for the two VVV CL001 stars with the BACCHUS code (Masseron et al. 2016). The spectra of our sample, in general, have a signal-tonoise (S/N) that is appropriate for elemental-abundance determinations; see Table 1. The atmospheric parameters (T eff , log g, and ξ t ), metallicity ([Fe/H]), and elemental-abundance ratios ([X/Fe]) were derived with the BACCHUS code. Figure 2 shows the good quality of the APOGEE-2 spectrum compared to the best-fit synthetic spectra for selected atomic and molecule lines of the two stars in VVV CL001, from which the chemical-abundance ratios were determined. For each spectrum we were able to identify reliable lines for six chemical species (N, O, Mg, Al, Si, and Fe). The same figure also reveals the scarcely detectable Fe I lines. The resulting elemental-abundance ratios are listed in Table 1.
The resulting chemical abundances are also displayed in Figure 3, and compared to four GCs at similar metallicity taken from the APOGEE-2 GC sample of Mészáros et al. (2020). We find that the two program stars exhibit an iron abundance ratio of −2.47 to −2.44, suggesting that VVV CL001 has a mean metallicity [Fe/H] = −2.45, with an uncertainty due to systematics of 0.24 dex, which makes this cluster possibly the most metal-poor GC identified so far within the Sun's Galactocentric distance, and with an extreme metallicity close to the apparent "floor" in the empirical metallicity distribution function for GCs in the MW and Local Universe (e.g., Geisler et al. 1995;Simpson 2018;Kruijssen 2019;Larsen et al. 2020;Wan et al. 2020). Note that VVV CL001 has the lowest Fe abundance amongst this sample, which includes the lowest metallicity GCs observed by APOGEE. With our limited sample, we do not find evidence for an intrinsic Fe-abundance spread.
Regarding the α-elements (O, Mg, and Si ), VVV CL001 displays a modest α-element enhancement, a clear signature of the fast enrichment provided by supernovae (SNe) II events, and compatible with other Galactic metal-poor GCs at similar metallicity (see Figure 3). Furthermore, no Mg-Al anti-correlation is evident in Figure 3. Elemental abundance of stars in VVV CL001 (black symbols) and four comparison metal-poor GCs. Each violin representation show the univariate kernel density estimate of the abundance ratios of each cluster as determined from Mészáros et al. (2020). Each GC has been slightly offset horizontally to distinguish them. our small sample; the two VVV CL001 stars exhibit an aluminum deficit, which places them within the definition of first-generation stars according to the criteria developed by Mészáros et al. (2020). However, one star in our sample displays a very high enrichment in nitrogen, and is considered to likely belong to some of the families of second-generation 2 stars with low-aluminum enrichment (see e.g. Mészáros et al. 2020), indicating the possible evidence for multiple stellar populations (MPs) in VVV CL001. However, we caution the reader that the error bars in either N or Al abundance (or both) are not properly understood or inconclusive, therefore more cluster members need to be followed-up of in order to confirm or discard the existence of MPs in this VVV CL001.
We made use of the state-of-art MW model -GravPot16 to predict the orbital path of VVV CL001 in a steady-state gravitational Galactic model that includes a "boxy/peanut" bar structure.
For the orbit computations, we adopt the same Galactic model configuration, solar position and velocity vector as described in Fernández-Trincado et al. (2020c), except for the angular velocity of the bar, for which we employed the recommended value of 41 km s −1 kpc −1 (Sanders et al. 2019).
The most likely orbital parameters and their uncertainties are estimated using a simple Monte Carlo approach. An ensemble of ten thousand orbits were calculated under variations of the observational parameters according to their estimated errors (assumed as 1−σ variation), where the errors were assumed to follow a Gaussian distribution.
To compute the orbits, we adopt a mean RV of −325.95 ± 6.6 km s −1 computed from the combined APOGEE-2 RVs of the two program stars and Hol-  Note-The errors were determined in the same manner, as described in Fernández-Trincado et al. (2020a) by variating the atmospheric parameters one at a time by the typical, albeit conservative, values of ∆T eff = ±100 K, ∆ log g = ±100 cgs, and ∆ξt = ±0.05 km s −1 . Thus, the reported uncertainties are defined as ger Baumgardt's RV compilation 3 of 34 stars, which lie in the RV range between −341.13 to −310.22 km s −1 . From those stars, we select objects having a Renormalized Unit Weight Error (RUWE) below 1.4, as extracted from Gaia EDR3, which allows us to discard sources with problematic astrometric solutions (see e.g. Lindegren et al. 2018). This reduces the 36 stars with RV information to 27 stars with both reliable proper-motions and RV information, which were considered to compute the nominal proper motion of VVV CL001 after applying a 3−σ clipping to the data, e.g., only 25 out of these 27 sources lies inside 3−σ of the nominal proper motion of the cluster, as shown in panel (b) of Figure  1 . From this procedure, the nominal proper motion of VVV CL001 is (µ α cos(δ), µ δ ) = (−3.41, −1.97) mas yr −1 , with an assumed uncertainty of 0.5 mas yr −1 . The heliocentric distance (d ) of VVV CL001 remains uncertain; for this reason we assume three possible estimate (5.5, 8.0, and 10.5 kpc) close to the best isochrone fit, as shown in panel (f) of Figure 1. The main orbital elements and the ensemble of orbits of VVV CL001 are displayed in panels (a) to (c) and (d) of Figure 4. The orbital parameters reveal that VVV CL001 lies on a radial and highly eccentric (>0.8) 3 https://people.smp.uq.edu.au/HolgerBaumgardt halo-like orbit with rather small excursions above the Galactic plane (Z max < 3 kpc), and pericentric (r peri ) distances below 1 kpc. We also find that VVV CL001 exhibits a retrograde and prograde sense at the same time when assuming a close heliocentric distance, while it is exclusively retrograde at and beyond the Bulge, similar to other GC GCs in the inner Galaxy (see Pérez-Villegas et al. 2020). Therefore, a more robust heliocentric distance estimation will better constraint these dynamics scenarios.
Panel (d) of Figure 4 shows that the orbital energy configuration of VVV CL001 is comparable to that of Galactic GCs associated with the major accretion events such as Sequoia (Seq.) and GES (see e.g., Myeong et al. 2018;Massari et al. 2019), indicating that VVV CL001 could be the fossil relic of one of these accreted dwarf galaxies. Baumgardt & Hilker (2018) performed N -body simulations of star clusters, and found that they could reproduce the surface-density profile of VVV CL001, finding a present-day mass of 9×10 4 M .

MASS
With the available RV data, we match the line-ofsight dispersion profiles to the N -body simulations of VVV CL001, as shown in panel (e) of Figure 4, and thus determine the most likely mass of the cluster from kinematics constraint. We adopted three radial bins (with bin centers of 6 , 15 , and 30 ), and chosen to ensure that at least ten stars were in each bin; resulting in the three blue dots shown in panel (e) of 4. With the new data in panel (e) of Figure 4, we find σ 0 ∼ 6.6 km s −1 . This yields a present-day estimated mass of ∼2.1×10 5 M , which suggests that VVV CL001 is two times more massive than previously thought, and approximately four times more massive than ESO280-SC06 (the most metalpoor GC known in the MW).

CONCLUDING REMARKS
We have performed the first near-IR high-resolution spectral analysis of two likely members of the GC VVV CL001, a cluster obscured by the heavy extinction and high field-star density in the direction of the Galactic Bulge. Based on high S/N APOGEE spectra, we measure a mean [Fe/H] metallicity of −2.45 ± 0.1, which makes VVV CL001 possibly the most metal-poor GC known, only slightly higher than ESO280-SC06 at −2.48.
VVV CL001 is very close in projection to UKS 1, which motivated Minniti et al. (2011) to suggest the possibility that they could be gravitationally bound. However, the RV of VVV CL001 is too large, and the orbits are very different, which allows us to rule out the binarycluster scenario.
We find that the α-element abundances of VVV CL001 are typical of the lowest metallicity GCs known. Spectra for more members are required in order to confirm the presence of MPs in this cluster.
For the derivation of the age and reddening, we employed the new code SIRIUS (Souza et al. 2020). As shown in panel (e) of Figure 1, we derived a median age of ∼ 11.9 +3.12 −4.05 Gyr, indicating that the very metal-poor GC VVV CL001 is among the oldest and most massive (∼2.1×10 5 M ) MW clusters.
A dynamical analysis of VVV CL001 reveals that this object has a radial and highly eccentric halo-like orbit confined inside the Sun's galactocentric distance. Both its metallicity and orbit favor the interpretation of VVV CL001 being a GC that belongs to an early accretion event in the MW corresponding to the either the Seq or GES dwarf galaxies. Finally, multi-band photometry in the near-IR will be useful to identify other stellar tracers, such as RR Lyrae stars (if any) toward VVV CL001, which will significantly help to constrain the distance, age, and origin of the cluster.
We thank the anonymous referee for helpful comments that greatly improved the paper. We warmly thank Holger Baumgardt for providing his published numerical N -body modeling of the line-of-sight velocity dispersion of VVV CL001. J.G.F-T is supported by