Issue |
A&A
Volume 642, October 2020
|
|
---|---|---|
Article Number | L19 | |
Number of page(s) | 8 | |
Section | Letters to the Editor | |
DOI | https://doi.org/10.1051/0004-6361/202039233 | |
Published online | 22 October 2020 |
Letter to the Editor
VVVX-Gaia discovery of a low luminosity globular cluster in the Milky Way disk
1
Departamento de Ciencias Físicas, Facultad de Ciencias Exactas, Universidad Andrés Bello, Fernández Concha 700, Las Condes, Santiago, Chile
e-mail: elisaritagarro1@gmail.com
2
Vatican Observatory, Vatican City State 00120, Italy
3
Centro de Astronomía (CITEVA), Universidad de Antofagasta, Av. Angamos 601, Antofagasta, Chile
4
Millennium Institute of Astrophysics, Santiago, Chile
5
Departamento de Astronomía, Av. Cisternas 1200 Norte, La Serena, Chile
6
Universidade de São Paulo, IAG, Rua do Matão 1226, Cidade Universitária, São Paulo 05508-900, Brazil
7
Observatorio Astronómico, Universidad Nacional de Córdoba, Laprida 854, Córdoba, Argentina
8
Gemini Observatory, Northern Operations Center, 670 A’ohoku Place, Hilo, HI 96720, USA
9
Instituto de Alta Investigación, Universidad de Tarapacá, Casilla 7D, Arica, Chile
10
European Southern Observatory, Karl-Schwarszchild-Strasse 2, 85748 Garching bei Muenchen, Germany
11
Department of Astronomy, University of Hertfordshire, Hertfordshire, UK
12
Mount Saint Vincent University, Halifax, Nova Scotia, Canada
13
Saint Mary’s University, Halifax, Nova Scotia, Canada
14
Instituto de Astronomía, Universidad Católica del Norte, Av. Angamos 0610, Antofagasta, Chile
15
Departamento de Física, Universidade Federal de Santa Catarina, Trinidade 88040-900 Florianópolis, SC, Brazil
16
Institute of Astronomy, University of Cambridge, Madingley Road, Cambridge CB3 0HA, UK
17
Instituto de Astrofísica de Canarias, Vía Láctea s/n, 38200 La Laguna, Tenerife, Spain
18
INAF-Osservatorio Astronomico di Capodimonte, Salita Moiariello 16, 80131 Naples, Italy
19
Universidad de Atacama, Copiapó, Chile
Received:
21
August
2020
Accepted:
29
September
2020
Context. Milky Way globular clusters (GCs) are difficult to identify at low Galactic latitudes because of high differential extinction and heavy star crowding. The new deep near-infrared (IR) images and photometry from the VISTA Variables in the Via Láctea Extended Survey (VVVX) allow us to chart previously unexplored regions.
Aims. Our long term aim is to complete the census of Milky Way GCs. The immediate goals are to estimate the astrophysical parameters for the newly discovered GC candidates, measuring their reddenings, extinctions, distances, total luminosities, proper motions, sizes, metallicities, and ages.
Methods. We used the near-IR VVVX survey database, in combination with the optical photometry and proper motions (PMs) from Gaia Data Release 2 (DR2), and with the Two Micron All Sky Survey photometry to search for and characterise new GCs within the southern Galactic plane (|b| < 5°).
Results. We report the detection of a heretofore unknown Galactic GC at RA = 14:09:00.0; Dec = −65:37:12 (J2000) corresponding to l = 310.828 deg; and b = −3.944 deg in galactic coordinates. We calculate a reddening of E(J − Ks) = (0.3 ± 0.03) mag and an extinction of AKs = (0.15 ± 0.01) mag for this new GC. Its distance modulus and corresponding distance were measured as (m − M) = (15.93 ± 0.03) mag and D = (15.5 ± 1.0) kpc, respectively. We also estimate the metallicity and age by comparison with known GCs and by fitting PARSEC and Dartmouth isochrones, finding [Fe/H] = ( − 0.70 ± 0.2) dex and t = (11.0 ± 1.0) Gyr. The mean GC PMs from Gaia DR2 are μα* = ( − 4.68 ± 0.47) mas yr−1 and μδ = ( − 1.34 ± 0.45) mas yr−1. The total luminosity of our cluster is estimated to be MKs = ( − 7.76 ± 0.5) mag. The core and tidal radii from the radial density profile are rc ∼ 2.1′ (4.6 pc) and rt = 6.5′ (14.6 pc) at the cluster distance.
Conclusions. We have found a new low luminosity, old and metal-rich GC, situated in the far side of the Galactic disk at RG = 11.2 kpc from the Galactic centre and at z = 1.0 kpc below the plane. Interestingly, the location, metallicity, and age of this GC are coincident with the Monoceros ring structure.
Key words: Galaxy: disk / Galaxy: stellar content / globular clusters: general / infrared: stars / surveys
© ESO 2020
1. Introduction
While most Galactic globular clusters (GCs) have been found using optical wavelength observations, large infrared surveys have made it possible to discover many more of them in recent decades. The main surveys are the Two Micron All Sky Survey (2MASS) observed in both hemispheres with two 2 m telescopes, one located at the Cerro Tololo Inter-American Observatory and the other at Kitt Peak National Observatory (Skrutskie et al. 2006; Cutri et al. 2003), the Galactic Legacy Infrared Midplane Survey Extraordinaire (GLIMPSE) based on observations with the NASA/ESA Spitzer Space Telescope (Benjamin et al. 2005), the UKIDSS Galactic Plane Survey (GPS) based on observations at the UKIRT 4 m telescope in Hawaii (Lucas et al. 2008), the NASA Wide Infrared Satellite Experiment (WISE; Wright et al. 2010), and the VISTA Variables in the Via Láctea Survey based on observations at the VISTA 4 m telescope (VVV; Minniti et al. 2010; Saito et al. 2012). This last survey, in particular, was expanded, becoming the VISTA Variables in the Via Láctea eXtended Survey (VVVX) with observations ranging from 2017 to 2021, also using the Visual-Infrared Camera (VIRCAM) at the VISTA telescope at ESO Paranal Observatory. The VVV has already enabled the discovery and characterisation of numerous new clusters (e.g. Minniti et al. 2011, 2017; Borissova et al. 2011; Moni Bidin et al. 2011; Ivanov et al. 2017; Gran et al. 2019; Camargo & Minniti 2019), and its extension VVVX is delivering on the promise to find even more in the extended area (e.g. Borissova et al. 2018).
However, the search for new GCs is still an intricate undertaking, since many considerations about the formation and evolution of GCs have to be included in this type of study. Clusters are not likely to contribute to the stellar halo much (e.g. Martell et al. 2011; Koch et al. 2019; Reina-Campos et al. 2020), but it is known that the clusters are disrupted, especially those with an orbit going through the bulge and disk (e.g. Baumgardt & Makino 2003; Kruijssen et al. 2011). Hence in the dense disk and bulge areas, there may be many small clusters that have lost part of their mass or had fewer stars from the beginning given the cluster initial mass function that predicts many more low mass clusters that are forming than are observed today among old GCs. Therefore, the detection of some GCs may be complicated due to their low luminosities and also the presence of dust and field stars that probably cover these objects.
2. Discovery of Garro01
The search and physical characterisation of Galactic GCs is one of the scientific goals of our VVVX survey, which already yields some results (Borissova et al. 2018, and in prep.; Barbá et al. 2019; Minniti et al., in prep.; Garro et al., in prep.; Obasi et al., in prep.). Here, we report on the discovery and confirmation of VVVX-GC-140900-653712 (hereafter Garro01 for short), a distant new GC deeply embedded in the Milky Way (MW) disk. Following the same procedure described by Minniti et al. (2017), this cluster was found in the VVVX database as a clear overdensity of red giant stars above the background. This GC lies in a complex low-latitude region, with variable extinction, high stellar density, containing evidence for star formation towards the Galactic plane, and also close (at ∼1°) in projection to the Circinus Seyfert 2 galaxy that is very extended on the sky (Freeman et al. 1977). From the first visual analysis, for example, Garro01 looks like a low luminosity GC compared with our previous finding of VVV−GC−05, which is a metal-poor [Fe/H] = −1.3 dex GC located in the same region of the sky at l = 330°, b = −1.8°, and D = 7.5 kpc (Minniti et al. 2017; Contreras Ramos et al. 2018). We also observe Garro01 as a clear excess of stars above the background in the optical, where the Gaia Data Release 2 (DR2; Gaia Collaboration 2018a) source density map of the region (Fig. A.1) clearly shows a nearly circular region. We have visually selected it, similar to a circle with r ∼ 2.5′ centred at equatorial coordinates RA = 14:09:00.0 and Dec = −65:37:12 (J2000), corresponding to Galactic coordinates l = 310.8278°, b = −3.9442°.
3. Physical characterisation of the new GC Garro01
We used a combination of near-infrared (IR) and optical data, obtained with the VVVX, the 2MASS, and Gaia DR2 surveys in order to measure the main parameters for Garro01. While we downloaded the available data from the respective repositories for 2MASS and Gaia DR2, we used a preliminary version of the point spread function (PSF) photometry for VVVX, which we are developing for its whole footprint based on Alonso-García et al. (2018) VVV PSF photometry. This VVVX PSF photometry was already used by Barbá et al. (2019) for the discovery and characterisation of FSR-1758, a giant metal-poor retrograde GC (Villanova et al. 2019), and by Minniti et al. (2018) for the discovery of extra tidal RR Lyrae stars in the metal-poor bulge GC NGC 6266 (M 62). For our analysis, we put the VVVX PSF photometry of our region of interest in the 2MASS JHKs photometric system.
We built a clean, decontaminated catalogue of highly-probable cluster members, benefiting from the precise astrometry and proper motions (PMs) from Gaia DR2, but also matching 2MASS+Gaia and VVVX+Gaia catalogues as a means to include both brighter and fainter sources (since stars with Ks < 11 mag are saturated in VVVX photometry). We first discarded all nearby stars with a parallax > 0.5 mas. Then, we inspected the vector PM (VPM) diagram (Fig. A.2), which shows a sharp peak on top of the broad stellar background distribution, which we interpret as the cluster mean PM. These mean PMs, as measured by Gaia DR2, are as follows: μα* = ( − 4.68 ± 0.47) mas yr−1 and μδ = ( − 1.35 ± 0.45) mas yr−1. Therefore, as probable cluster members, we selected all stars within 1 mas yr−1 from the mean cluster PMs. Figure A.2 shows the VPM diagrams for sources with Ks < 15 mag matched in the 2MASS+Gaia catalogues and sources with Ks > 13 mag matched in the VVVX+Gaia catalogues, respectively. In the same figure, the stars located outside of the 1 mas yr−1 circle selected on the VPM diagram, but within the selected radius of 2.5′ of the cluster centre, are shown in yellow. These stars are predominantly coming from the surrounding field that includes disk foreground and background stars.
The Gaia+2MASS colour-magnitude diagrams (CMDs) are shown in Fig. A.3, highlighting the tight red giant branch (RGB) of the PM-selected cluster members. The fact that the cluster RGB is preferentially fainter and redder than the field stars indicates that the field stars mostly belong to a foreground, less reddened population in the Galactic plane. This rules out a window in the dust distribution; as such, windows were also detected similar to stellar overdensities (e.g. Minniti 2018; Saito et al. 2020). For comparison, the fiducial RGB from the GC 47 Tuc (taken from Gaia Collaboration 2018b; Cohen et al. 2015) is overplotted in red in GaiaG versus (BP − RP) and in VVVX Ks versus (J − Ks), respectively; this was applied for all other colour indexes (moving it to ΔG = 3.8 mag and Δ(BP − RP) = 0.39 mag), showing excellent agreement with the RGB of Garro01. We bring attention to the overdensity of stars around G ∼ 17.70 in the Gaia CMD, which is well fitted with the slanted line corresponding to the location of the red clump (RC) core He burning stars in 47 Tuc.
Figure A.4 shows the luminosity functions of G and Ks bands for a 3′ radius field centred on the cluster for all the red giant stars with J − Ks > 0.7 mag (top panel), as well as for the PM-selected red giants (bottom panel). These luminosity functions clearly exhibit the peaks due to the cluster RC giants at GRC = (17.70 ± 0.05) mag and KsRC = (14.48 ± 0.05) mag, respectively.
The resulting PM-cleaned CMDs from the Gaia+2MASS+VVVX data are displayed in Fig. A.5. These CMDs show that the RGB is narrow and well defined, and the cluster RC is clearly seen. It is also important to note the presence of blue star residual contamination. These stars are too bright to be blue stragglers belonging to Garro01, and we argue that this contamination arises from a fraction of Galactic foreground field stars with similar PMs as the GC.
We estimate the reddening and extinction in this field following the maps of Ruiz-Dern et al. (2018) in the near-IR, those of Schlafly & Finkbeiner (2011) in the optical, and by taking advantage of the mean magnitude of RC giants. In detail, we assume the RC stars calibration of Ruiz-Dern et al. (2018), where the absolute magnitude in the Ks-band is MKs = ( − 1.601 ± 0.009) mag and the relative colour is (J − Ks)0 = (0.66 ± 0.02) mag, obtaining E(J − Ks) = (0.30 ± 0.03) mag and AKs = (0.15 ± 0.01) mag, respectively. In particular, we note that within ∼10′ from the GC centre, the reddening seems to be fairly uniform, with variations of ΔE(J − Ks) < 0.05 mag.
The GC distance modulus can then be measured by adopting AKs/E(J − Ks) = 0.5 mag (Minniti et al. 2018), which yields (m − M)0 = (15.93 ± 0.03) mag for Garro01, and therefore the heliocentric distance is D = (15.5 ± 1.0) kpc, placing this GC at RG = 11.2 kpc from the Galactic centre (Bland-Hawthorn & Gerhard 2016), assuming R0 = 8.3 kpc (Dékány et al. 2013). We would like to point out that we also compared our results with other K-band calibration methods in order to achieve more robust parameters. We used RC magnitude theoretical calibrations from Salaris & Girardi (2002) and Alves et al. (2002), finding a perfect agreement ((m − M)0 = 15.87 mag, D = 15.0 kpc, and (m − M)0 = 15.93 mag, D = 15.3 kpc, respectively). Simultaneously, we derived the parameters using Gaia DR2 photometry. Adopting the extinction AKs = (0.15 ± 0.01) mag, we obtain the equivalent value in the G-band of AG = (1.2 ± 0.2) mag (employing the V-band extinction value AV = 1.4 mag from Schlafly & Finkbeiner 2011) and so the reddening of E(BP − RP) = (0.39 ± 0.06) mag. These values were used to measure the distance modulus of (m − M)0 = (16.0 ± 0.2) mag and the heliocentric distance D = (15.9 ± 1.0) kpc, which is in excellent agreement with the 2MASS+VVVX photometry. The distance measurement of D = 15.5 kpc yields a scale of 1′ = 2.25 pc, placing this cluster at a height of about 1.0 kpc below the Galactic plane.
Furthermore, we carried out the radial density profile of Garro01 in order to determine its size. First of all, we checked the central position, because an offset in its position may alter the radial density profile. Indeed, we found that the new centre is RA = 212.23 deg and Dec = −65.628 deg (shifted to ΔRA = 2.1′ and ΔDec = 0.48′ from the initial coordinates). Using the resulting centre, we divided our sample into six radial bins (out to a radius of 1.5′), and we calculated the area in each bin. After that, we derived the density inside the circular annuli as the total number of stars over the area. Finally, we overplotted the King profile (King 1962) that well-reproduces density points; obtaining a core radius of rc = (2.1 ± 1.5) arcmin is equivalent to a physical size of rc = (4.6 ± 3.1) pc and a tidal radius of arcmin, corresponding to
pc, which is consistent with the typical galactic GC sizes as listed in the 2010 Harris (1996) compilation. We note that the Poisson errors are very large due to a poor statistics, and the background level was supposed to be zero since the catalogue is assumed to be clean from contaminants.
The integrated GC absolute magnitude is estimated by coadding the RGB stars from the PM decontaminated diagrams and employing the GC size of 6.5 arcmin. We find MKs = ( − 7.76 ± 0.5) mag, which is equivalent to MV = ( − 5.26 ± 1.0) mag for typical GC mean colours of V − K = (2.5 ± 0.94) mag. Certainly, this represents a lower limit since sub-giant branch and main-sequence stars are below the observational limit. Indeed, benefiting from a strong resemblance 47 Tuc, we also estimated the total luminosity for Garro01. Briefly, we calculated the integrated absolute magnitude, coadding the RGB 47 Tuc stars and scaling to the value MV = −9.42 mag from the 2010 Harris (1996) catalogue. This results in an absolute magnitude of MV = −5.62 mag for Garro01, placing this GC on the low luminosity tail of the MW GCLF, ∼2 mag fainter than the peak of the MW GCLF (MV = ( − 7.4 ± 0.2) mag from Harris 1991; Ashman & Zepf 1998).
Finally, the cluster metallicity was derived following two different methods. Firstly, we compared Garro01 with the fiducial GC 47 Tuc (age t = 11.8 Gyr, [Fe/H] = −0.72 dex, [α/Fe]= + 0.4 dex, and D = 4.5 kpc from Brogaard et al. 2017), making reddening adjustments in order for evolutionary sequences of both GCs to coincide. Secondly, we used the fitting-isochrone method, preferring both PARSEC (Marigo et al. 2017) and Dartmouth (Dotter et al. 2008) isochrones. In order to have a good fit of the stellar isochrones in the near-IR and optical CMDs (Fig. A.5), we adopted the resulting values of reddening, extinction, and distance modulus calculated from the VVVX and Gaia photometries. Both methods provide a metallicity of [Fe/H] = ( − 0.70 ± 0.1) dex and α-enhanced between 0 and +0.4 dex. However, the cluster appears to be slightly younger than 47 Tuc because we estimate an age of t = (11 ± 1) Gyr (from stellar isochrones). Briefly, the best fitting age and metallicity was obtained by iteratively fitting isochrones by comparing our data with isochrones generated with different ages and metallicities and selecting the best by-eye fit. First, we fixed the age and varied the metallicity. Then we fixed the metallicity and searched for the best fitting age. In this way, we also deduced an estimation of the metallicity and age errors, until the fitting-isochrones did not reproduce all evolutionary sequences in both optical and near-IR CMDs. It is important to note that this is a rough estimation of the age since the main-sequence turn-off is below to the magnitude lower limit.
In addition to the comparison with the GC 47 Tuc, it is appropriate to make some comparisons with other known GCs. For this reason, we considered Lynga 7 ([Fe/H] = −1.01 dex, and D = 8.0 kpc, Bica et al. 2016), and NGC 5927 ([Fe/H] = −0.32 dex and D = 8.2 kpc, Pérez-Villegas et al. 2020), which are located in the same region of the Galactic plane, and also the 9 Gyr old open cluster NGC 6791 ([Fe/H] = +0.3 dex, D = 4.1 kpc, Boesgaard et al. 2009). These comparisons do not match as well as 47 Tuc, which exhibits obvious differences in the colour of the RGB, in particular, and also the position of the RC.
We also searched for more evidence of an old GC population, such as the RR Lyrae stars, finding two RR Lyrae candidates from the Gaia DR2 catalogue (Clementini et al. 2019) within 10′ of this GC. These are Gaia DR2 5851208911444645504, located 8.5 arcmin away from the GC centre at RA = 14:10:16.06 and Dec = −65:40:28.6 (J2000), with a mean magnitude of G = 16.94 mag, an amplitude of A = 0.73 mag, and a period of P = 0.570059 days; and Gaia DR2 5851209426894334592, located 8.6 arcmin away at RA = 14:10:23.49 and Dec = −65:36:56.0 (J2000), with a mean magnitude of G = 16.33 mag, an amplitude of A = 0.94 mag, and a period of P = 0.595121 days. Both RR Lyrae are compatible with fundamental mode pulsators (RRab type). Their distances were measured using the period-luminosity relation from Clementini et al. (2019). We believe that only Gaia DR2 5851208911444645504 may be a real cluster member, with a heliocentric distance of DRRL = 15.9 kpc; whereas for Gaia DR2 5851209426894334592, we calculated a much shorter distance of DRRL = 12.4 kpc. However, for both stars, their Gaia DR2 PMs of μα* = ( − 9.685 ± 0.082) mas yr−1 and μδ = ( − 0.603 ± 0.102) mas yr−1 as well as μα* = ( − 5.866 ± 0.085) mas yr−1 and μδ = ( − 4.190 ± 0.098) mas yr−1, respectively, are > 4σ different from the mean GC PM. This result calls into question the cluster membership for both variable stars, therefore, until the RR Lyrae membership is better established with additional data, we prefer to be cautious and adopt the distance measured using the RC giants, D = 15.5 kpc. The physical parameters of the new GC Garro01 are summarised in Table A.1, along with their respective uncertainties (1σ).
4. Discussion
The discovery of this new GC demonstrates that the Galactic GC census is incomplete, so many more objects similar to this may still be found at low Galactic latitudes, hidden in very crowded, heavily reddened regions (e.g. Ivanov et al. 2000a,b, 2005; Borissova et al. 2007; Kurtev et al. 2008; Longmore et al. 2011; Valcheva et al. 2015). The known low luminosity GCs are very rare and located at high Galactic latitudes (Van Den Bergh 2003), and the low latitude ones would be particularly difficult to find.
In addition, we have to consider that the history of GCs is complex, as well as that of our Galaxy. It is known that the MW has suffered from past merging events (e.g. Ibata et al. 1994; Helmi et al. 1999, 2018; Belokurov et al. 2018; Myeong et al. 2018). Further, photometric detections of stellar streams and substructures, especially in the Galactic halo (Bell et al. 2008), are consistent with predictions from cosmological simulations that link stellar streams to accreted dwarf satellite galaxies (Bullock & Johnston 2005). As a consequence, many GCs may be associated with these merging events (e.g. Massari et al. 2019; Vasiliev 2019).
A known structure close to the position of Garro01 is the Monoceros ring (MRi, which is also known as the Galactic anti-centre stellar structure, GASS). It was found by Newberg et al. (2002) and could represent a fitting example of a past merging event. However, the nature of MRi is still a topic of discussion. Martin et al. (2005), Peñarrubia et al. (2005), Morganson et al. (2016), and Guglielmo et al. (2017) modelled the remnant of a tidally-disrupted satellite galaxy, the Canis Major dwarf galaxy (Ibata et al. 2003; Yanny et al. 2003; Crane et al. 2003; Rocha-Pinto et al. 2003; Frinchaboy et al. 2006; Grillmair 2006; Conn et al. 2008; Sollima et al. 2011). These simulations reveal that the MRi feature should be present in the field studied here, as shown, for example, by Conn et al. (2008, see their Figs. 12 and 13). Also, spectroscopic studies have revealed that this feature is moderately metal-poor, with a mean [Fe/H] = −0.8 dex (Li et al. 2012), which is in agreement with the result presented here of [Fe/H] = ( − 0.7 ± 0.1) dex. However, an alternative scenario (e.g. Momany et al. 2006; Hammersley & López-Corredoira 2011; Kalberla et al. 2014; Sheffield et al. 2018) provides support for the idea that the MRi is merely the Galactic warp and flare, composed of stars from the MW disk kicked out to their current location due to interactions between a satellite galaxy and the disk.
These results are controversial, revealing the complexity of the regions under study. This fact makes it all the more valuable to have star clusters that may be associated with the MRi structure. Indeed, Frinchaboy et al. (2006) searched for open and GCs associated with the MRi structure, concluding that Berkeley 29, Saurer 1, Tombaugh 2, Arp-Madore 2, Palomar 1, NGC 2808, and NGC 5286 were possibly associated, while BH176 is not. We note that the classification of BH 176 as a GC has been recently questioned by Sharina et al. (2014) and Vásquez et al. (2018), and it might be a unique case. On the other hand, other studies (Carraro et al. 2007; Moitinho et al. 2006) concluded that old open clusters are not consistent with the models for the MRi, which is probably just related to the Galactic warp and flare. Using more recent photometric and spectroscopic observations, Frinchaboy et al. (2014) argue that Berkeley 29 ([Fe/H] = −0.44 dex, D = 13.4 kpc, and RG = 21.1 kpc) and Saurer 1 ([Fe/H] = −0.38 dex, D = 13.1 kpc, and RG = 20.2 kpc) are associated with the MRi, while again BH176 is not.
Our main goal is to understand the real nature of Garro01 and derive all parameters that describe it. For this reason, we searched for a possible association to other GCs in that region, including Lynga 7, NGC 5927, 47 Tuc, BH176, and others. As mentioned in Sect. 2, we found very close similarity with 47 Tuc, so we conclude that Garro01 is a new Galactic GC. Also, given the metallicity ([Fe/H] = −0.7 dex) and the location (D = 15.5 kpc, RG = 11.2 kpc) of Garro01 in close proximity of the MRi structure, it is appealing to consider its possible association with the MRi.
The field star contamination, in the same field where Garro01 is located, overwhelms the GC stars, but the vast majority of the field stars along this line of sight seem to be brighter and less reddened and, therefore, are located in the foreground (Fig. A.3). They also have PMs whose mean values are different from those of the new Garro01 (Fig. A.2). Nonetheless, it is crucial to obtain radial velocities (RVs) for this cluster in order to constrain its orbital properties and to definitely confirm its association with the MRi structure. Also, we note that, according to the maps of Momany et al. (2006), there should be no warp in the specific direction of the field studied here at a Galactic longitude of l = 311 deg.
Unfortunately, reliable Gaia DR2 RVs are not available for this cluster. A thorough search revealed that there is an excess of stars with RV > 200 km s−1 in this field, although all of these sources are too bright to be cluster members. Although the lack of RV measurements prevent us from calculating the orbit, the low Galactic height (z ∼ − 1.1 kpc) and the very small vertical proper motion directed towards the Galactic plane (μb = +0.223 mas yr−1) suggest a possible association with the old Galactic disk. The combination of old age, a high metallicity, and a disk orbit closely resembles 47 Tuc and three other GCs (namely M 107, NGC 6362, and E3), whose polar paths cross in two very well-defined points on the sky (de la Fuente Marcos et al. 2015). Nevertheless, we checked that the polar path of Garro01 has no similarity with theirs, and it does not pass near these intersections. Therefore, the association of this new GC with the possible Sagittarius stream proposed by de la Fuente Marcos et al. (2015) should therefore be excluded.
5. Conclusions
We report the discovery of the GC Garro01 in the VVVX near-IR images. This is a low luminosity GC, located at a low latitude in the MW disk, with large reddening and high foreground stellar density. We complement our VVVX data set with 2MASS near-IR photometry as well as Gaia optical photometry and PMs.
The optical and near-IR photometry allowed us to construct CMDs in order to estimate principal GC parameters, such as age, metallicity, distance, absorption, and reddening. Most field stars are less reddened and therefore located in the foreground disk.
The optical and near-IR CMDs show that this new GC is more metal-poor and older than the canonical old metal-rich open cluster NGC 6791 (D = 4.1 kpc, t = 9 Gyr, and [Fe/H] = 0.3 dex, Boesgaard et al. 2009). On the other hand, we find a close similarity with the physical parameters of 47 Tuc (D = 4.5 kpc, t = 11.8 Gyr, and [Fe/H] = −0.7 dex, Brogaard et al. 2017), arguing for an old GC at D = 15.5 kpc. However, this cluster is intrinsically 4.16 mag fainter than 47 Tuc, implying that its total mass is about M ≈ 104 M⊙ (if 47 Tuc has M = 7 × 105 M⊙, Marks & Kroupa 2010), on the low mass end of the Galactic GCs (Baumgardt et al. 2019). We also note the presence of two candidate RR Lyrae stars within 10 arcmin from the cluster centre. One of them has a distance consistent with the measured RC distance for the cluster. Its PMs, however, suggest that both RR Lyrae are not cluster members.
The Galactic GC census is incomplete at low latitudes (Ivanov et al. 2005; Minniti et al. 2017). Furthermore, the discovery of Garro01 suggests that more low luminosity GCs may still be uncovered at low Galactic latitudes.
Another interesting implication is that the location, distance, and metallicity of this GC match those of the MRi structure, a potential accretion event recently identified in the Galactic plane. The possibility of a physical association with this structure must be confirmed with follow-up spectroscopy to measure radial velocities and detailed chemical abundances.
Acknowledgments
ERG acknowledges support from an UNAB PhD scholarship. J.A.-G. acknowledges support from Fondecyt Regular 1201490 and from ANID, Millennium Science Initiative ICN12_009, awarded to the Millennium Institute of Astrophysics (MAS). DM and MG are supported by Proyecto FONDECYT Regular No. 1170121. DM is also supported by the BASAL Center for Astrophysics and Associated Technologies (CATA) through grant AFB 170002. We gratefully acknowledge the use of data from the ESO Public Survey program IDs 179.B-2002 and 198.B-2004 taken with the VISTA telescope and data products from the Cambridge Astronomical Survey Unit. PWL acknowledges support by STFC Consolidated Grant ST/R00905/1. SRA acknowledges support from the FONDECYT Iniciacion project 11171025, the FONDECYT Regular project 1201490, and the CONICYT + PAI “Concurso Nacional Insercion de Capital Humano Avanzado en la Academia 2017” project PAI 79170089. R.K.S. acknowledges support from CNPq/Brazil through project 305902/2019-9.
References
- Alonso-García, J., Saito, R. K., Hempel, M., et al. 2018, A&A, 619, A4 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Alves, D. R., Rejkuba, M., Minniti, D., & Cook, K. H. 2002, ApJ, 573, L51 [NASA ADS] [CrossRef] [Google Scholar]
- Ashman, K. M., & Zepf, S. E. 1998, Cambridge Astrophysics Series (Cambridge: Cambridge University Press), 30 [Google Scholar]
- Barbá, R. H., Minniti, D., Geisler, D., et al. 2019, ApJ, 870, L24 [NASA ADS] [CrossRef] [Google Scholar]
- Baumgardt, H., & Makino, J. 2003, MNRAS, 340, 227 [NASA ADS] [CrossRef] [Google Scholar]
- Baumgardt, H., Hilker, M., Sollima, A., & Bellini, A. 2019, MNRAS, 482, 5138 [Google Scholar]
- Bell, E. F., Zucker, D. B., Belokurov, V., et al. 2008, ApJ, 680, 295 [NASA ADS] [CrossRef] [Google Scholar]
- Belokurov, V., Erkal, D., Evans, N. W., Koposov, S. E., & Deason, A. J. 2018, MNRAS, 478, 611 [NASA ADS] [CrossRef] [Google Scholar]
- Benjamin, R. A., Churchwell, E., Babler, B. L., et al. 2005, ApJ, 630, L149 [NASA ADS] [CrossRef] [Google Scholar]
- Bica, E., Ortolani, S., & Barbuy, B. 2016, PASA, 33, e028 [NASA ADS] [CrossRef] [Google Scholar]
- Bland-Hawthorn, J., & Gerhard, O. 2016, ARA&A, 54, 529 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Boesgaard, A. M., Jensen, E. E. C., & Deliyannis, C. P. 2009, AJ, 137, 4949 [NASA ADS] [CrossRef] [Google Scholar]
- Borissova, J., Ivanov, V. D., Stephens, A. W., et al. 2007, A&A, 474, 121 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Borissova, J., Bonatto, C., Kurtev, R., et al. 2011, A&A, 532, A131 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Borissova, J., Ivanov, V. D., Lucas, P. W., et al. 2018, MNRAS, 481, 3902 [CrossRef] [Google Scholar]
- Brogaard, K., VandenBerg, D. A., Bedin, L. R., et al. 2017, MNRAS, 468, 645 [NASA ADS] [CrossRef] [Google Scholar]
- Bullock, J. S., & Johnston, K. V. 2005, ApJ, 635, 931 [NASA ADS] [CrossRef] [Google Scholar]
- Camargo, D., & Minniti, D. 2019, MNRAS, 484, L90 [CrossRef] [Google Scholar]
- Carraro, G., Moitinho, A., Zoccali, M., Vázquez, R. A., & Baume, G. 2007, AJ, 133, 1058 [NASA ADS] [CrossRef] [Google Scholar]
- Clementini, G., Ripepi, V., Molinaro, R., et al. 2019, A&A, 622, A60 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Cohen, R. E., Hempel, M., Mauro, F., et al. 2015, AJ, 150, 176 [CrossRef] [Google Scholar]
- Conn, B. C., Lane, R. R., Lewis, G. F., et al. 2008, MNRAS, 390, 1388 [NASA ADS] [Google Scholar]
- Contreras Ramos, R., Minniti, D., Gran, F., et al. 2018, ApJ, 863, 79 [NASA ADS] [CrossRef] [Google Scholar]
- Crane, J. D., Majewski, S. R., Rocha-Pinto, H. J., et al. 2003, ApJ, 594, L119 [NASA ADS] [CrossRef] [Google Scholar]
- Cutri, R. M., Skrutskie, M. F., van Dyk, S., et al. 2003, VizieR Online Data Catalog: II/246 [Google Scholar]
- de la Fuente Marcos, R., de la Fuente Marcos, C., Moni Bidin, C., Ortolani, S., & Carraro, G. 2015, A&A, 581, A13 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Dékány, I., Minniti, D., Catelan, M., et al. 2013, ApJ, 776, L19 [NASA ADS] [CrossRef] [Google Scholar]
- Dotter, A., Chaboyer, B., Jevremović, D., et al. 2008, ApJS, 178, 89 [NASA ADS] [CrossRef] [Google Scholar]
- Freeman, K. C., Karlsson, B., Lynga, G., et al. 1977, A&A, 55, 445 [NASA ADS] [Google Scholar]
- Frinchaboy, P. M., Muñoz, R. R., Phelps, R. L., Majewski, S. R., & Kunkel, W. E. 2006, AJ, 131, 922 [NASA ADS] [CrossRef] [Google Scholar]
- Frinchaboy, P. M., Thompson, B. A., O’Connell, J., et al. 2014, Am. Astron. Soc. Meet. Abstr., 223, 442.12 [Google Scholar]
- Gaia Collaboration (Brown, A. G. A., et al.) 2018a, A&A, 616, A1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Gaia Collaboration (Babusiaux, C., et al.) 2018b, A&A, 616, A10 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Gran, F., Zoccali, M., Contreras Ramos, R., et al. 2019, A&A, 628, A45 [CrossRef] [EDP Sciences] [Google Scholar]
- Grillmair, C. J. 2006, ApJ, 651, L29 [NASA ADS] [CrossRef] [Google Scholar]
- Guglielmo, M., Lane, R. R., Conn, B. C., et al. 2017, MNRAS, 474, 4584 [CrossRef] [Google Scholar]
- Hammersley, P. L., & López-Corredoira, M. 2011, A&A, 527, A6 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Harris, W. E. 1991, ARA&A, 29, 543 [NASA ADS] [CrossRef] [Google Scholar]
- Harris, W. E. 1996, AJ, 112, 1487 [NASA ADS] [CrossRef] [Google Scholar]
- Helmi, A., White, S. D. M., de Zeeuw, P. T., & Zhao, H. 1999, Nature, 402, 53 [NASA ADS] [CrossRef] [Google Scholar]
- Helmi, A., Babusiaux, C., Koppelman, H. H., et al. 2018, Nature, 563, 85 [NASA ADS] [CrossRef] [Google Scholar]
- Ibata, R. A., Gilmore, G., & Irwin, M. J. 1994, Nature, 370, 194 [NASA ADS] [CrossRef] [Google Scholar]
- Ibata, R. A., Irwin, M. J., Lewis, G. F., Ferguson, A. M. N., & Tanvir, N. 2003, MNRAS, 340, L21 [NASA ADS] [CrossRef] [Google Scholar]
- Ivanov, V. D., Borissova, J., & Vanzi, L. 2000a, A&A, 362, L1 [NASA ADS] [Google Scholar]
- Ivanov, V. D., Borissova, J., Alonso-Herrero, A., & Russeva, T. 2000b, AJ, 119, 2274 [NASA ADS] [CrossRef] [Google Scholar]
- Ivanov, V. D., Kurtev, R., & Borissova, J. 2005, A&A, 442, 195 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Ivanov, V. D., Piatti, A. E., Beamín, J.-C., et al. 2017, A&A, 600, A112 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Kalberla, P. M. W., Kerp, J., Dedes, L., & Haud, U. 2014, ApJ, 794, 90 [NASA ADS] [CrossRef] [Google Scholar]
- King, I. 1962, AJ, 67, 471 [NASA ADS] [CrossRef] [Google Scholar]
- Koch, A., Grebel, E. K., & Martell, S. L. 2019, A&A, 625, A75 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Kruijssen, J. M. D., Pelupessy, F. I., Lamers, H. J. G. L. M., Portegies Zwart, S. F., & Icke, V. 2011, MNRAS, 414, 1339 [NASA ADS] [CrossRef] [Google Scholar]
- Kurtev, R., Ivanov, V. D., Borissova, J., & Ortolani, S. 2008, A&A, 489, 583 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Li, J., Newberg, H. J., Carlin, J. L., et al. 2012, ApJ, 757, 151 [NASA ADS] [CrossRef] [Google Scholar]
- Longmore, A. J., Kurtev, R., Lucas, P. W., et al. 2011, MNRAS, 416, 465 [NASA ADS] [Google Scholar]
- Lucas, P. W., Hoare, M. G., Longmore, A., et al. 2008, MNRAS, 391, 136 [NASA ADS] [CrossRef] [Google Scholar]
- Marigo, P., Girardi, L., Bressan, A., et al. 2017, ApJ, 835, 77 [NASA ADS] [CrossRef] [Google Scholar]
- Marks, M., & Kroupa, P. 2010, MNRAS, 406, 2000 [NASA ADS] [Google Scholar]
- Martell, S. L., Smolinski, J. P., Beers, T. C., & Grebel, E. K. 2011, A&A, 534, A136 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Martin, N. F., Ibata, R. A., Conn, B. C., et al. 2005, MNRAS, 362, 906 [NASA ADS] [CrossRef] [Google Scholar]
- Massari, D., Koppelman, H. H., & Helmi, A. 2019, A&A, 630, L4 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Minniti, D. 2018, in The Vatican Observatory, Castel Gandolfo: 80th Anniversary Celebration, eds. G. Gionti, & J. B. Kikwaya Eluo, 51, 63 [Google Scholar]
- Minniti, D., Lucas, P. W., Emerson, J. P., et al. 2010, New Astron., 15, 433 [Google Scholar]
- Minniti, D., Hempel, M., Toledo, I., et al. 2011, A&A, 527, A81 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Minniti, D., Palma, T., Dékány, I., et al. 2017, ApJ, 838, L14 [NASA ADS] [CrossRef] [Google Scholar]
- Minniti, D., Saito, R. K., Gonzalez, O. A., et al. 2018, A&A, 616, A26 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Moitinho, A., Vázquez, R. A., Carraro, G., et al. 2006, MNRAS, 368, L77 [NASA ADS] [CrossRef] [Google Scholar]
- Momany, Y., Zaggia, S., Gilmore, G., et al. 2006, A&A, 451, 515 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Moni Bidin, C., Mauro, F., Geisler, D., et al. 2011, A&A, 535, A33 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Morganson, E., Conn, B., Rix, H.-W., et al. 2016, ApJ, 825, 140 [CrossRef] [Google Scholar]
- Myeong, G. C., Evans, N. W., Belokurov, V., Sanders, J. L., & Koposov, S. E. 2018, ApJ, 856, L26 [CrossRef] [Google Scholar]
- Newberg, H. J., Yanny, B., Rockosi, C., et al. 2002, ApJ, 569, 245 [NASA ADS] [CrossRef] [Google Scholar]
- Peñarrubia, J., Martínez-Delgado, D., Rix, H. W., et al. 2005, ApJ, 626, 128 [NASA ADS] [CrossRef] [Google Scholar]
- Pérez-Villegas, A., Barbuy, B., Kerber, L. R. O., et al. 2020, MNRAS, 491, 3251 [Google Scholar]
- Reina-Campos, M., Hughes, M. E., Kruijssen, J. M. D., et al. 2020, MNRAS, 493, 3422 [NASA ADS] [CrossRef] [Google Scholar]
- Rocha-Pinto, H. J., Majewski, S. R., Skrutskie, M. F., & Crane, J. D. 2003, ApJ, 594, L115 [NASA ADS] [CrossRef] [Google Scholar]
- Ruiz-Dern, L., Babusiaux, C., Arenou, F., Turon, C., & Lallement, R. 2018, A&A, 609, A116 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Saito, R. K., Hempel, M., Minniti, D., et al. 2012, A&A, 537, A107 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Saito, R. K., Minniti, D., Benjamin, R. A., et al. 2020, MNRAS, 494, L32 [CrossRef] [Google Scholar]
- Salaris, M., & Girardi, L. 2002, MNRAS, 337, 332 [NASA ADS] [CrossRef] [Google Scholar]
- Schlafly, E. F., & Finkbeiner, D. P. 2011, ApJ, 737, 103 [NASA ADS] [CrossRef] [Google Scholar]
- Sharina, M. E., Donzelli, C. J., Davoust, E., Shimansky, V. V., & Charbonnel, C. 2014, A&A, 570, A48 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Sheffield, A. A., Price-Whelan, A. M., Tzanidakis, A., et al. 2018, ApJ, 854, 47 [CrossRef] [Google Scholar]
- Skrutskie, M. F., Cutri, R. M., Stiening, R., et al. 2006, AJ, 131, 1163 [Google Scholar]
- Sollima, A., Martínez-Delgado, D., Valls-Gabaud, D., & Peñarrubia, J. 2011, ApJ, 726, 47 [NASA ADS] [CrossRef] [Google Scholar]
- Valcheva, A. T., Ovcharov, E. P., Lalova, A. D., et al. 2015, MNRAS, 446, 730 [CrossRef] [Google Scholar]
- Van Den Bergh, S. 2003, ApJ, 590, 797 [NASA ADS] [CrossRef] [Google Scholar]
- Vasiliev, E. 2019, MNRAS, 484, 2832 [NASA ADS] [CrossRef] [Google Scholar]
- Vásquez, S., Saviane, I., Held, E. V., et al. 2018, A&A, 619, A13 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Villanova, S., Monaco, L., Geisler, D., et al. 2019, ApJ, 882, 174 [NASA ADS] [CrossRef] [Google Scholar]
- Wright, E. L., Eisenhardt, P. R. M., Mainzer, A. K., et al. 2010, AJ, 140, 1868 [Google Scholar]
- Yanny, B., Newberg, H. J., Grebel, E. K., et al. 2003, ApJ, 588, 824 [NASA ADS] [CrossRef] [Google Scholar]
Appendix A: Figures and table
![]() |
Fig. A.1. Left panel: Gaia stellar density map for a r = 2 degree field in Galactic coordinates, indicating the new GC with an arrow, along with the size of the HI emission of the Circinus galaxy (Freeman et al. 1977). Middle and right panel: density maps of the zoomed 42′×42′ region around Garro01 using optical Gaia data and VVVX near-IR data, respectively, where the circle indicates the approximate cluster size. It is important to note that the redder areas are representative of overdensities, while the yellower areas are lower densities. |
![]() |
Fig. A.2. Vector PM diagrams for the bright sources (Ks < 15 mag) matched in the 2MASS+Gaia catalogues (right panel), and the fainter sources (Ks > 13 mag) matched in the VVVX+Gaia catalogues (left panel). The black circle indicates the cluster selection, while blue crosses represent the position of the two RR Lyrae found within 10′ of the cluster centre. The colour bars indicate a higher (towards the black colour) and lower (towards the yellow colour) concentration. |
![]() |
Fig. A.3. Optical and near-IR CMDs for the GC Garro01 field from 2MASS (black points) and Gaia or VVVX (yellow points), showing the PM-selected cluster members. We note that the well-defined cluster RGB located lower and to the red side of the diagrams, indicating that the fiducial RGB from GC 47 Tuc from Cohen et al. (2015) and Gaia Collaboration (2018b) is overplotted in the top panels, is appropriately shifted to the cluster distance and reddening; this is a distant and reddened cluster. |
![]() |
Fig. A.4. Top panel: luminosity functions for GC Garro01, showing the VVVX Ks-band photometry in red on the left, and the Gaia DR2 G-band photometry in blue on the right. The arrows point to the location of the RC. Bottom panel: same diagram showing only the counts of PM-selected sources. |
![]() |
Fig. A.5. PM decontaminated CMDs for the GC Garro01, showing the VVVX near-IR CMD (top left), the Gaia DR2 optical CMD (top right), and the Gaia-VVVX optical-IR CMDs (bottom panels). Red points are VVVX-Gaia matched sources, while grey points are 2MASS-Gaia matched sources. The black dotted line shows the fit to a PARSEC isochrone with a metallicity of [Fe/H] = −0.70 dex and an age of t = 11.0 Gyr. We note that the well defined cluster RGB and RC are present in all of the CMDs, as is also the case for the residual field contamination on the blue side arising from a fraction of Galactic foreground field stars with similar PMs as the GC. |
Final GC physical parameters.
All Tables
All Figures
![]() |
Fig. A.1. Left panel: Gaia stellar density map for a r = 2 degree field in Galactic coordinates, indicating the new GC with an arrow, along with the size of the HI emission of the Circinus galaxy (Freeman et al. 1977). Middle and right panel: density maps of the zoomed 42′×42′ region around Garro01 using optical Gaia data and VVVX near-IR data, respectively, where the circle indicates the approximate cluster size. It is important to note that the redder areas are representative of overdensities, while the yellower areas are lower densities. |
In the text |
![]() |
Fig. A.2. Vector PM diagrams for the bright sources (Ks < 15 mag) matched in the 2MASS+Gaia catalogues (right panel), and the fainter sources (Ks > 13 mag) matched in the VVVX+Gaia catalogues (left panel). The black circle indicates the cluster selection, while blue crosses represent the position of the two RR Lyrae found within 10′ of the cluster centre. The colour bars indicate a higher (towards the black colour) and lower (towards the yellow colour) concentration. |
In the text |
![]() |
Fig. A.3. Optical and near-IR CMDs for the GC Garro01 field from 2MASS (black points) and Gaia or VVVX (yellow points), showing the PM-selected cluster members. We note that the well-defined cluster RGB located lower and to the red side of the diagrams, indicating that the fiducial RGB from GC 47 Tuc from Cohen et al. (2015) and Gaia Collaboration (2018b) is overplotted in the top panels, is appropriately shifted to the cluster distance and reddening; this is a distant and reddened cluster. |
In the text |
![]() |
Fig. A.4. Top panel: luminosity functions for GC Garro01, showing the VVVX Ks-band photometry in red on the left, and the Gaia DR2 G-band photometry in blue on the right. The arrows point to the location of the RC. Bottom panel: same diagram showing only the counts of PM-selected sources. |
In the text |
![]() |
Fig. A.5. PM decontaminated CMDs for the GC Garro01, showing the VVVX near-IR CMD (top left), the Gaia DR2 optical CMD (top right), and the Gaia-VVVX optical-IR CMDs (bottom panels). Red points are VVVX-Gaia matched sources, while grey points are 2MASS-Gaia matched sources. The black dotted line shows the fit to a PARSEC isochrone with a metallicity of [Fe/H] = −0.70 dex and an age of t = 11.0 Gyr. We note that the well defined cluster RGB and RC are present in all of the CMDs, as is also the case for the residual field contamination on the blue side arising from a fraction of Galactic foreground field stars with similar PMs as the GC. |
In the text |
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.