Kinematic data rebuild the nuclear star cluster as the most metal-rich region of the Galaxy

Context. The Galactic centre (GC) is located at only 8kpc from Earth and constitutes a unique template with which to understand Galactic nuclei. Nevertheless, the high crowding and extinction towards the GC hamper the study of its main stellar components, the nuclear stellar disc (NSD) and the nuclear star cluster (NSC). Aims. Recent work has suggested that the NSD and the NSC can be distinguished along the line of sight towards the NSC via the di ﬀ erent extinction of their stars. This motivated us to analyse the proper motion, radial velocity, and the metallicity distributions of the di ﬀ erent extinction groups. Methods. We use photometric, kinematic, and metallicity data to distinguish between probable NSD and NSC stars in a region centred on the NSC. Results. We detected two di ﬀ erent extinction groups of stars and obtained that they have signiﬁcantly di ﬀ erent proper motion distributions, in agreement with the expected kinematics for the NSD and the NSC. We also derived radial velocity maps that appear to be di ﬀ erent for the two components, and found di ﬀ erent metallicities for each of them, with the largest one measured for the most extinguished group of stars. We obtained that the metallicity distribution of each extinction group is best ﬁtted by a bimodal distribution, indicating the presence of two metallicity components for each of them (a broad one slightly below solar metallicity, and a more metal-rich, narrower one that is largest for the high-extinction group of stars). Conclusions. We conclude that both extinction groups are distinct GC components with di ﬀ erent kinematics and metallicity, and correspond to the NSD and the NSC. Therefore, it is possible to distinguish them via their di ﬀ erent extinction. The high mean metal-licity, [M / H] ∼ 0 . 3dex, obtained for the metal-rich stars of the NSC supports the proposition that this component is arguably the most metal-rich region of the Galaxy.


Introduction
The centre of the Milky Way is the closest galaxy nucleus and the only one where we can resolve individual stars down to milliparsec scales. Besides the supermassive black hole, Sagittarius A*, two main stellar structures outline the Galactic centre (GC): (1) the nuclear star cluster (NSC), a massive stellar cluster (∼ 2.5 × 10 7 M , e.g. Launhardt et al. 2002;Schödel et al. 2014;Feldmeier et al. 2014) placed at the heart of the Galaxy with an effective radius of ∼ 5 pc (e.g. Graham & Spitler 2009;Schödel 2011;Feldmeier-Krause et al. 2017;Gallego-Cano et al. 2020), and (2) the nuclear stellar disc (NSD), a much larger stellar structure (with a radius ∼ 200 pc, e.g. Launhardt et al. 2002;Nishiyama et al. 2013;Gallego-Cano et al. 2020;Sormani et al. 2020Sormani et al. , 2022 surrounding the NSC, and partially overlapping with the dense gas from the central molecular zone (e.g. Henshaw et al. 2022).
In spite of being placed at only 8 kpc from Earth (e.g. Gravity Collaboration et al. 2018;Do et al. 2019), the study of the GC is hampered by the high extinction and the extreme source crowding (e.g. Nishiyama et al. 2008;Schödel et al. 2010;Nogueras-Lara et al. 2018a, 2020b, 2021b. Therefore, disentangling the GC components and analysing their main properties is a formidable challenge. Recent studies point towards a different stellar population and formation scenario for the NSC and the NSD (Nogueras-Lara et al. 2020a;Schödel et al. 2020;Nogueras-Lara et al. 2021c). In this way, despite being domi-nated by old stars (more than 80 % of their stellar population is older than 8 Gyr), the NSC contains a stellar population of intermediate age (up to 15 % of the stellar mass was formed ∼ 3 Gyr ago; Schödel et al. 2020) that is not present in the NSD. Analogously, a significant mass fraction of the NSD stars (∼ 5 %) formed about 1 Gyr ago, corresponding to a time during which the NSC experienced no significant star formation activity (Pfuhl et al. 2011;Nogueras-Lara et al. 2020a;Schödel et al. 2020). On the other hand, the metallicity of the NSC stars seems to be higher than that of the NSD (Schultheis et al. 2021;Feldmeier-Krause 2022).
Recently, Nogueras-Lara et al. (2021c) showed that NSD and NSC stars are subject to significantly different reddening. Here, we follow up on their work by including kinematic and spectroscopic data. We find that the kinematics and metallicity distributions are different for each of the extinction groups, and that they are in agreement with the expected results for the NSD and the NSC. coordinates 17 h 45 m 38 s , -29 • 00 12 ). The J data belong to the GALACTICNUCLEUS survey (GNS, Nogueras-Lara et al. 2018a, a high-angular resolution (∼ 0.2 ) JHK s catalogue specially designed to observe the GC. The K s -band data were obtained in 2013 and form part of a pilot study for the GNS survey (Nogueras-Lara et al. 2018a) that was obtained under excellent observing conditions (K s seeing ∼ 0.4 ), allowing photometry ∼ 1 mag deeper than for the corresponding GNS field (see Table in Nogueras-Lara et al. 2018a). To correct potential saturation problems affecting bright stars in K s (e.g. Nogueras-Lara et al. 2019), we used the SIRIUS/IRSF survey (Nagayama et al. 2003;Nishiyama et al. 2006) to replace the K s photometry of stars with K s < 11.5 mag (for further details, see Nogueras-Lara et al. 2020a).

Proper motions
We used a publicly available proper motion catalogue of the GC (Shahzamanian et al. 2019(Shahzamanian et al. , 2022, that overlaps with the region analysed in this paper. This catalogue was specifically designed to study the GC and constitutes an unprecedented kinematic data set for the NSD region. The proper motions were computed using two data sets (the GNS H-band data, Nogueras-Lara et al. 2018a, and the HST Paschen-α survey, Wang et al. 2010;Dong et al. 2011), with a timeline of ∼ 7−8 years between them. The high angular resolution of the used data sets (∼ 0.2 Shahzamanian et al. 2022) allows the catalogue to supersede previous surveys covering the same area (e.g. the VIRAC survey Smith et al. 2018) -that are limited to a narrow magnitude range in the analysed region because of saturation and seeing-limited resolution (Smith et al. 2018;Shahzamanian et al. 2022).
Given the complex procedure of matching the data from the two different epochs, their different tiling patterns, and the different size of the detectors, the proper motion catalogue does not have a homogeneous coverage across its surveyed field (see Fig. 2 in Shahzamanian et al. 2022). Nevertheless, we verified that 80 % of the region analysed in this study is properly covered as shown in Fig. 1.
The computed proper motions are relative and not absolute, meaning that they were computed assuming an average zero motion between the stellar positions from the used reference catalogue to calculate the proper motions. This only means that there is an offset between proper motions considering an absolute relative frame and the ones in the catalogue. Shahzamanian et al. (2022) carried out a detailed comparison of their catalogue with absolute calibrated work (Libralato et al. 2021), and concluded that the proper motions perfectly agree once converted to the reference frame. In any case, in this paper, we do not pursue the calculation of absolute proper motions, but aim to analyse the proper motion distribution of the NSC and the NSD. Given that we are using the same pointings for the NSC and the NSD (instead of comparing different GC regions), the use of relative proper motions does not pose any problem.

Metallicity and radial velocity
We used metallicity and radial velocity KMOS (Sharples et al. 2013) medium-resolution spectra data for around 1 000 stars in the NSC region (see Fig. 1) obtained by Feldmeier-Krause et al. (2017. These latter authors derived the radial velocities using the IDL routine pPXF (Cappellari & Emsellem 2004) using the high-resolution spectra of Wallace & Hinkle (1996) as a template. They then applied a full spectral fitting method with PHOENIX models (Husser et al. 2013) to derive the metallicities and other stellar parameters. The absolute metallicities were calibrated using existing empirical spectra for [M/H] < 0.3 dex. This might produce an overestimation of metallicities [M/H] > 0.3 dex. To check this, Schultheis et al. (2021) derived the metallicities using the data from Feldmeier-Krause et al. (2020), applying an alternative method using CO and NaI index, and empirical spectra calibrated for [M/H] < 0.6 dex (for further details of the method, see Fritz et al. 2021). Schultheis et al. (2021) concluded that there is no systematic bias for stars with [M/H] < 0.5 dex, whereas stars with [M/H] > 0.5 dex show significantly higher metallicities when applying the methodology of Feldmeier-Krause et al. (2017). On the other hand, this metallicity cut at [M/H] = 0.5 dex also agrees with the most metal-rich stars detected in the NSC when using high-resolution spectroscopy (e.g. Do et al. 2015;Rich et al. 2017;Thorsbro et al. 2020).

Disentangling the stellar populations
Recent work by Nogueras-Lara et al. (2021c) used HK s photometry to show that the stellar populations from the NSC and the NSD can be identified along the line of sight towards the NSC via their significantly different extinction (A K s NS D ∼ 1.7 mag, and A K s NS C ∼ 2.3 mag, Nogueras-Lara et al. 2021c). In the present work, we use J data instead of H to improve the identification of the NSD and the NSC, given that the extinction is higher for shorter wavelengths (e.g. Nogueras-Lara et al. 2021b), and thus the use of J-band data increases the difference between both stellar populations. Figure 2 shows the colour-magnitude diagram (CMD) K s versus J − K s , where we visually detect two red clump (stars in their helium burning core sequence; Girardi 2016) bumps with different extinction (see Fig. 16 of Nogueras-Lara et al. 2018a). These bumps correspond to the NSD (lowextinction group) and the NSC (high-extinction group), according to previous work (Nogueras-Lara et al. 2021c). To analyse the proper motion distribution of each of the extinction groups, we cross-correlated the HAWK-I photometric data and the proper motion catalogue (Shahzamanian et al. 2022). Figure 3 shows a CMD K s versus J − K s of the stars with available proper motions in the field. We found about 7 000 common stars and checked that more than 90 % of them have proper motion uncertainties below 0.7 mas/yr for both components of the proper motions (µ l and µ b , where l and b refer to Galactic longitude and latitude, respectively), as indicated in the right panels of Fig. 3.
We defined two boxes in the CMD (Fig. 3, left panel) according to the colour cut at J − K s ∼ 5 mag, which divides the two red clump bumps with different extinction (Fig. 2). The colour cut applied to distinguish between each extinction group accounts for the predominantly old stellar population ( 80 % of the stellar mass) that dominates the NSD and the NSC (Nogueras-Lara et al. 2020a; Schödel et al. 2020;Nogueras-Lara et al. 2021c). In this way, to reduce the confusion between both extinction groups, the shape of the boxes defined in Fig. 3 is in agreement with the slope of an old stellar isochrone (for further details, see Fig. 4 in Nogueras-Lara et al. 2018b). Moreover, the blue cut of the first box (J − K s ∼ 4 mag) was chosen to remove the foreground stellar population, which mainly belongs to the Galactic disc (and to some extent to the Galactic bulge/bar), and presents a significantly lower extinction than GC stars (Nogueras-Lara et al. 2021a). We also defined a brightness cut for each box at K s = 14.2 mag, according to the detection limit of the proper motion catalogue. Figure 4 shows the comparison between the proper motion distribution of each of the extinction groups (Fig. 3). We observed a clearly different distribution of the proper motions parallel to the Galactic plane (µ l ) between both groups. In this way, the stars belonging to the low extinction group present a larger µ l in comparison to those from the high-extinction group. On the other hand, the distribution of the proper motion component perpen- dicular to the Galactic plane (µ b ) appears to be similar for both groups.

Proper motion distribution
We further analysed the proper motion distribution using the SCIKIT-LEARN Python function GaussianMixture (GMM, Pedregosa et al. 2011) to obtain the probability density function that gives rise to the underlying distributions based on Gaussian models. First of all, we applied the Akaike information criterion (AIC, Akaike 1974) to obtain the number of Gaussian models that best describe the probability density function of the data. Following the results of Shahzamanian et al. (2022) for the NSD, we tried up to three Gaussian models for the proper motion component parallel to the Galactic plane (µ l ), and two models for the perpendicular component (µ b ). We obtained that the µ l distribution of the low-extinction group of stars is best described by the combination of two Gaussian models, whereas a three-Gaussian model is favoured for the high-extinction group. The µ b distribution of both groups is similar and best fitted by a single-Gaussian model. Figure 5 shows the obtained results. To estimate the parameters and the associated uncertainties of each Gaussian model, we resorted to a Jackknife approach repeating the GMM modelling 1 000 times, randomly dropping 20 % of the used stars for each iteration. Table 1 shows the obtained results.

Discussion
We find that the µ l distribution of each of the extinction groups is significantly different. Therefore, their stellar populations have different kinematics, reinforcing the conclusion that they correspond to different GC components that we identify as the NSD and the NSC. We explain the presence of two Gaussian models in the low-extinction group because of the observation of the stellar population from the NSD that is in front of the NSC. In this way, the Gaussian component with µ l = 2.25 ± 0.10 mas/yr corresponds to the NSD stellar population moving eastwards, and agrees within the uncertainties with the value obtained by Shahzamanian et al. (2022) when analysing the NSD kinematics (µ l = 1.99 ± 0.13 mas/yr). On the other hand, the secondary Gaussian (µ l = 0.23 ± 0.21 mas/yr) is probably due to stars belonging to the Galactic bulge/bar that present a broader distri- bution, in agreement with previous work (Clarkson et al. 2008;Kunder et al. 2012;Soto et al. 2014;Shahzamanian et al. 2022). A third Gaussian component corresponding to the stellar population moving westwards (e.g. Shahzamanian et al. 2022) is not visible due to the presence of the NSC, due to its high stellar density and also the increase in interstellar reddening towards the farther edge of the NSD.
We interpret the three Gaussian components detected in the µ l distribution of the high-extinction stellar group as the result of rotation of the NSC, which produces two main peaks due to stars moving eastwards (µ l = 1.80 ± 0.06 mas/yr) and westwards (µ l = −2.50 ± 0.05 mas/yr). The third peak centred on µ l = −0.36 ± 0.04 mas/yr is probably due to the differential rotation of the NSC. In this way, stars rotating slower than the eastwards and westwards population producing the previously described peaks are likely giving rise to the central population. This is compatible with the presence of a radial dependence of the rotation velocity, which when measured was found to be smaller (∼ 0.5 mas/yr) for the innermost parsec of the NSC (Trippe et al. 2008;Schödel et al. 2009). Moreover, some contamination from the NSD and the Galactic bulge/bar might also contribute to the observed distribution. However, due to the higher density of stars present in the NSC in comparison with the NSD, in combination with the significantly larger extinction due to the given colour cut, we expect that the contamination from the Galactic bulge/bar will be less important for this stellar population than for the previous case of the NSD.
The fact that the third peak is not centred on zero is probably because the proper motions in the used catalogue were computed as relative ones assuming that the mean proper motion of   all the stars in a given field is zero (Shahzamanian et al. 2022).
In this way, due to the presence of the NSC in the studied region, it is not possible to observe stars from the far side of the NSD, and thus the reference frame for µ l is dominated by a majority of stars moving eastwards; the positive component of µ l is therefore underestimated, whereas the negative one is overestimated. This scenario is consistent with the larger absolute value of µ l obtained for the stars moving westwards in comparison with the stars moving eastwards (µ l = −2.50 ± 0.05 mas/yr, and 1.80 ± 0.06 mas/yr, respectively). More specifically, assuming that the µ l component of the proper motions is shifted µ l = +0.36 ± 0.04 mas/yr, we obtain that the components due to the rotation of the NSC are equal in absolute value within the uncertainties, and that the third Gaussian peak is centred on zero.
Analysing the µ l distribution of the NSC for different colour cuts, we verified that the stellar population moving westwards is more extinguished. This is compatible with the lower number of stars belonging to this Gaussian component obtained in the GMM modelling (Table 1), which is probably due to the extinction by dust within the NSC (Chatzopoulos et al. 2015). Therefore, we obtained that the direction of rotation is the same for both the NSC and the NSD, in agreement with previous work (e.g. Feldmeier et al. 2014;Schönrich et al. 2015;Shahzamanian et al. 2022). Moreover, the values obtained for the eastwardmoving component of the NSD and the NSC indicate a faster rotation of the NSD compared to the NSC, as expected (e.g. Fig. 18 in Sormani et al. 2022).
On the other hand, both extinction groups show similar distributions of µ b (except for a somewhat larger dispersion for stars from the high-extinction group), which can be best reproduced by a single-Gaussian model approximately centred on zero. Although the presence of some contamination from the Galactic bulge/bar in the low-extinction group (NSD) might have led to a two-Gaussian distribution (Shahzamanian et al. 2022), this was probably not observed because of the low number of stars from the Galactic bulge/bar in the sample. However, we measured a broadening of the detected single-Gaussian component (σ µb = 2.39 ± 0.05 mas/yr) in comparison with the expected one for the NSD according to previous work (σ µb = 1.499 ± 0.0002 mas/yr, Shahzamanian et al. 2022), which is probably due to the influence of the expected Galactic bulge/bar stars. This is in agreement with the results of Sormani et al. (2022), who estimated the contamination of stars from the Galactic bulge/bar for different fields distributed across the NSD. These authors obtained that the contribution from Galactic bulge/bar stars is less important for their fields closest to the NSC ( 20 % of the stars with H − K s > 1.3 mag belong to the Galactic bulge/bar stars; see Table 2 in Sormani et al. 2022). Although the regions analysed in Sormani et al. (2022) avoid the NSC, which is precisely the target field of this work, we expect that the relative contribution of Galactic bulge/bar will be lower than in their innermost fields, because of the presence of a significantly higher number of stars belonging to the NSC. On the other hand, in this work we use the near-infrared bands J and K s , which allow us to remove the foreground stellar population (from the Galactic disc and also from the bulge/bar) in a more efficient way than the H and K s bands used by Sormani et al. (2022), because of the longer wavelength base line.

Radial velocities
We analysed the distribution of radial velocities of spectroscopically characterised stars from Feldmeier-Krause et al. (2017. We searched for stars found in our photometric catalogue that also have known metallicities (about 1 000 common stars), and classified them following the boxes defined in Fig. 2. We obtained that 10 % of the stars with known metallicities belong to the low-extinction group. We used a photometric criterion instead of a proper motion analysis because of the low number of common detections between the proper motion catalogue and the stars with known metallicities (see the gaps of the proper motion catalogue in Fig. 1).
We built a radial velocity map for the stars belonging to each of the extinction groups. We defined a pixel size of ∼ 50 , and ∼ 25 for the radial velocity maps corresponding to the NSD and the NSC, given the significantly lower number of stars detected for the NSD (around 10 % of the total number of stars). We computed the radial velocity for a given pixel calculating the mean value of the radial velocities from the stars located within the pixel, imposing a 3-sigma clipping algorithm to remove outliers. We only computed a value for a given pixel if at least four stars were present within it. Figure 6 shows the obtained maps. The radial velocity distribution appears to be different between the NSD and the NSC. In this way, the NSD map is dominated by positive values. This might be due to the NSD geometry not yet being well known, the low number of stars used to produce this map (∼ 60 stars), and/or the possible contamination from the Galactic bulge/bar that can bias the results. On the other hand, the NSC map clearly shows a rotation pattern that is also compatible with the differential rotation suggested in Sect. 3.3 to explain the three Gaussian components observed for the µ l distribution.
Article number, page 5 of 8 A&A proofs: manuscript no. 44411corr . The x-and y-axes indicate the distance in arcseconds from Sagittarius A* (indicated by a black star). Cross-shaped pixels indicate that there are not enough stars to compute a radial velocity value. The pixel size is larger in the upper panel because of the low number of NSD stars in the region (∼ 60 stars).

Metallicity
We studied the metallicity distribution of the stars in the target region using the previously computed list of common stars between our photometric catalogue and stars with known metallicities (Feldmeier-Krause et al. 2017. Figure 7 represents the metallicity distribution of the stars corresponding to each of the extinction groups. Our results show different metallicity distributions for each extinction group, with the high-extinction group being more metal rich than the low-extinction group. This indicates the presence of two components with different metallicity and agrees with the results obtained for the NSD and the NSC (e.g. Feldmeier-Krause et al. 2017Fritz et al. 2021;Schultheis et al. 2021;Feldmeier-Krause 2022). The NSC metallicity obtained after removing the contribution from the NSD is somewhat higher than that obtained without considering the presence of the NSD (as done in previous work; e.g. Feldmeier-Krause et al. 2017). This is because NSD stars are less metal rich on average and bias the NSC sample when they are not removed.
We further study the metallicity distribution of each component restricting the analysis to stars with metallicity [M/H] < 0.5 dex. In this way, we removed stars whose metallicity might have been overestimated because of the calibration used in Feldmeier-Krause et al. 2017. We applied the GMM method explained in Sect. 3.2 to each of the extinction groups. We used the AIC to distinguish between models considering one and two Gaussians. We obtained that the metallicity distribution of both extinction groups is best represented by a two-Gaussian model (Fig. 8). To estimate the mean value and uncertainty of each component, we generated 1 000 Monte Carlo samples for the metallicity distribution, assuming Gaussian uncertainties for the metallicity of each star and randomly varying its value considering a Gaussian distribution. Table 2  Normalised #stars mean = 0.12 dex median = 0.15 dex = 0.35 dex mean = 0.29 dex median = 0.34 dex = 0.41 dex Fig. 7. Metallicity distribution of the stars belonging to the NSD (in blue) and the NSC (in orange). The mean, median, and standard deviation of each distribution are shown. The systematic uncertainty of the metallicity measurements is ∼ 0.25 dex according to Feldmeier-Krause et al. (2017. determined averaging over the results of the 1 000 Monte Carlo samples, applying a 3-sigma outlier-resistant criterion to remove outliers. For both extinction groups, the metallicity distributions present a broad component centred around [M/H] ∼ -0.20 dex, and a more metal-rich narrower component whose metallicity is higher for the NSC stellar population. Our results agree with the recent work by Schultheis et al. (2021), where two stellar components with different metallicities were found for both the NSC and the NSD, probably indicating differences in the formation of these two components. Part of the data used in Schultheis et al. (2021) to analyse the NSC corresponds to the same data set that we used in our analysis (Feldmeier-Krause et al. 2020). Nevertheless, Schultheis et al. (2021) assumed all the stars in the sample to be part of the NSC. Here, we verified that 10 % of the stars that we used belong to the NSD and have a lower mean metallicity in comparison to the stars from the NSC. Therefore, considering them as NSC stars might slightly bias the results towards lower metallicity values for the NSC. On the other hand, we are able to disentangle the stellar populations from the NSD and the NSC using the same data set and line of sight for both components, reinforcing the results obtained in Schultheis et al. (2021) when comparing different data sets and lines of sight.
High-resolution spectroscopic studies have previously been carried out for some stars in the NSC. In particular, Ryde et al. (2016) found a metal-poor ([Fe/H] ∼ −1.0 dex) red giant star at a projected distance of 1.5 pc from Sagittarius A*, and determined that it probably belongs to one of the nuclear components (NSD or NSC), confirming the presence of metal-poor stars in the GC. On the other hand, Do et al. (2018) targeted two stars found to be very metal rich in previous spectral-template fitting  work (Do et al. 2015). These authors confirmed that at least one of them has an unusually high metallicity [M/H] > 0.6 dex, confirming the presence of these kinds of stars in the GC. Moreover, Rich et al. (2017) determined the metallicity of 17 stars in the NSC to span the range of −0.5 < [Fe/H] < +0.5 dex, with these values also being compatible with the range obtained in this work. Finally, Thorsbro et al. (2020) found that the most metal-rich stars in their NSC sample reach [Fe/H] = +0.5 dex. These values end up in higher metallicities when considering the overall metallicity, [M/H], and are in agreement with the most metal-rich stars that we use from Feldmeier-Krause et al. (2017. We would like to stress that the use of JK s photometry and proper motions can be very helpful for reanalysing previous work where some of the stars that were considered to belong to the NSC might be part of the NSD.

Conclusion
In this study, we analysed the kinematics and metallicity of the stars belonging to two extinction groups identified along the line of sight towards the NSC (Nogueras-Lara et al. 2018a, 2021c). We detected two kinematically distinct components associated to the each of the extinction groups, which we confirmed as the NSD and the NSC. In this way, our results show the potential of proper motions to disentangle the stellar population belonging to each of the GC structures. We analysed their proper motion components parallel (µ l ) and perpendicular (µ b ) to the Galactic plane, and find that the µ l distributions of the NSD and the NSC are best fitted by two-and three-Gaussian models, respectively. In this way, we explained the NSD proper motion distribution as the combination of the stellar population from the closest edge of the NSD (stars rotating eastwards), and some contamination from Galactic bulge/bar stars that cannot be easily removed from the sample because of their large extinction. We concluded that the presence of the NSC impedes the detection of stars from the far side of the NSD that are rotating westwards (Shahzamanian et al. 2022). The µ l distribution of the NSC shows three Gaussian components, which we explain as a consequence of the rotation of the NSC. These components correspond to stars moving eastwards (positive µ l ) and westwards (negative µ l ), and stars moving with relatively slower velocities from the innermost regions of the NSC. We obtained relatively low values for the rotation of the NSC in comparison with the NSD, as expected (e.g. Sormani et al. 2022). The µ b distributions of both the NSD and the NSC seem to be similar except for a somewhat higher dispersion for the NSC values. We created radial velocity maps using spectroscopically characterised stars from each of the extinction groups and found that the velocity distributions seem to be different between the NSD and the NSC. On the other hand, we observed a velocity pattern for the NSC compatible with differential rotation.
We also analysed stars with known metallicities from each of the extinction groups. We found that they follow different distributions, with the group corresponding to the NSD being less metal rich, in agreement with previous studies (e.g. Feldmeier-Krause et al. 2017Fritz et al. 2021;Schultheis et al. 2021;Feldmeier-Krause 2022). We obtained that both components can be best described by a two-Gaussian model with a less metalrich, wider component, and a predominant metal-rich narrower component that might have a different origin, in agreement with Schultheis et al. (2021). Moreover, we measure a mean value of the NSC metal-rich stellar population of [M/H] ∼ 0.3 dex, which arguably makes the NSC the most metal-rich region of the Galaxy.
Our results confirm that the NSC and the NSD can be distinguished along the line of sight towards the NSC via their different extinction and agree with previous work on the NSD and the NSC suggesting different formation scenarios and stellar populations (e.g. Nogueras-Lara et al. 2020a;Schödel et al. 2020;Schultheis et al. 2021;Feldmeier-Krause 2022;Nogueras-Lara et al. 2021c).