Open Access
Issue
A&A
Volume 690, October 2024
Article Number A313
Number of page(s) 11
Section Galactic structure, stellar clusters and populations
DOI https://doi.org/10.1051/0004-6361/202450946
Published online 17 October 2024

© The Authors 2024

Licence Creative CommonsOpen Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

This article is published in open access under the Subscribe to Open model. Subscribe to A&A to support open access publication.

1 Introduction

The heart of the Milky Way stands out as our closest galaxy centre, and it is the only one where we can study stars down to milli-parsec scales (e.g. Dong et al. 2011; Nogueras-Lara et al. 2019; Do et al. 2019; Gravity Collaboration 2020). Therefore, it constitutes a unique template for analysing the stellar population of a galaxy nucleus, which can provide a better understanding of galactic nuclei and their role in galaxy formation and evolution. Nevertheless, high extinction and extreme source crowding hamper analysis of the Galactic centre stellar population and primarily limit it to the near-infrared regime (e.g. Nishiyama et al. 2006, 2008; Schultheis et al. 2009; Fritz et al. 2011; Alonso-García et al. 2017; Thorsbro et al. 2020; Nogueras-Lara et al. 2021b; Sanders et al. 2022b).

The Galactic centre is roughly outlined by the nuclear stellar disc (NSD), a flat and likely axisymmetric stellar structure with a scale length of ~ 100 pc and a scale height of ~40 pc (e.g. Gallego-Cano et al. 2020; Nogueras-Lara 2022a; Sormani et al. 2022). It contains a total stellar mass of ~ 109 M (Launhardt et al. 2002; Sormani et al. 2022) and is characterised by a predominantly old stellar population and a bursty star-formation history (e.g. Nogueras-Lara et al. 2020; Sanders et al. 2022a; Schödel et al. 2023), which makes it different from the surrounding and much larger Galactic bar/bulge (e.g. Zoccali et al. 2003; Valenti et al. 2013; Wegg & Gerhard 2013; Wegg et al. 2015; Portail et al. 2017; Surot et al. 2019). Additionally, the NSD contains ~106 M of young stars that formed in situ in the last 30Myr (e.g. Matsunaga et al. 2013; Nogueras-Lara et al. 2020, 2022; Nogueras-Lara 2024), and it has contributed to up to 10% of the star-forming activity observed throughout the entire Milky Way over the past ~100 Myr (e.g. Mezger et al. 1996; Mauerhan et al. 2010; Matsunaga et al. 2011; Crocker et al. 2011; Nogueras-Lara et al. 2020).

The NSD is a rotating structure (Lindqvist et al. 1992; Schönrich et al. 2015; Fritz et al. 2021; Shahzamanian et al. 2022) that appears to be kinematically distinct from the Galactic bar/bulge (Schönrich et al. 2015; Schultheis et al. 2021). However, Zoccali et al. (2024) recently reported that they were unable to definitively identify the NSD as a distinct kinematic feature when comparing it to nearby control fields in the Galactic bar/bulge, contradicting previous work. Their analysis used proper motions derived from the VISTA Variables in the Vía Láctea survey (Minniti et al. 2010) and APOGEE DR17 data (Majewski et al. 2017; Abdurro’uf et al. 2022).

The NSD stellar population exhibits a broad metallicity distribution that has been previously fitted using two Gaussian models accounting for metal-rich and metal-poor stars, respectively (Schultheis et al. 2021; Nogueras-Lara 2022b; Nogueras-Lara et al. 2023a). The metal-rich component (mean metallicity [M/H] = 0.12, Schultheis et al. 2021) constitutes approximately 70–80% of the NSD’s stars and shows a lower line-of-sight velocity dispersion in comparison with the metal-poor component (mean metallicity [M/H] = −0.22, Schultheis et al. 2021). This suggests that metal-rich stars have a different origin and probably formed in situ (Schultheis et al. 2021) from the dense accumulation of gas and dust that comprises the central molecular zone (CMZ; e.g. Henshaw et al. 2023). This is further supported by indications that stars in the NSD and the gas in the CMZ appear to be co-rotating (Schönrich et al. 2015; Schultheis et al. 2021).

In this paper, we aim to study the proper motion distribution of NSD stars with different metallicities to assess previous analysis using line-of-sight velocities and to confirm whether they have different kinematics and origin. Additionally, we utilise the recent orbital analysis of NSD stars by Nieuwmunster et al. (2024) to determine whether stars with different metallici-ties have different orbits. Analysing stars with different kinematics and orbits also helps shed light on the potential kinematic difference between the NSD and the Galactic bar/bulge, and it is crucial to understanding their relation as well as the NSD formation mechanism.

thumbnail Fig. 1

Scheme of the observed region overplotted on a Spitzer/IRAC image at 8 µm. (Stolovy et al. 2006). This wavelength highlights the bright and dusty clouds that outline the NSD. a) KMOS fields in which the colour scale indicates the star counts. b) Number of stars per field with computed orbits (Nieuwmunster et al. 2024). The two northernmost fields do not contain any stars due to the selection criteria used to create the orbital catalogue (|l| < 1.5° and |b| < 0.3°; for further details, see Nieuwmunster et al. 2024). The dotted contours show the 25, 50, and 75% levels of the NSD self-consistent model obtained by Sormani et al. (2022).

2 Data

2.1 Metallicity, photometry, and proper motions

We used the spectroscopic survey of the NSD (Fritz et al. 2021) obtained with the KMOS instrument at the VLT (Sharples et al. 2013). The chosen sample includes the primary targets of the survey contained within the central region of the NSD shown in Fig. 1 (2303 stars). The metallicities were derived with an uncertainty of ~0.3 dex and a calibration uncertainty of ~0.1 dex (Fritz et al. 2021).

We also used H and Ks photometry from the SIRIUS/IRSF survey of the Galactic centre (Nagayama et al. 2003; Nishiyama et al. 2006, 2013), which is available for all the targets of the KMOS spectroscopic survey of the NSD (see Table E1 in Fritz et al. 2021). The mean zero-point uncertainties are ~0.03 mag in both bands, and the 10 σ limiting magnitudes are H = 16.6 mag and Ks = 15.6 mag.

For the kinematical analysis, we used proper motions from a preliminary version of the second version of VIRAC (VIRAC2, Smith et al. 2024). This data set was also used by Sormani et al. (2022) and Nieuwmunster et al. (2024). It contains absolute proper motions of stars in the spectroscopic survey, using Gaia as an astrometric reference. Sormani et al. (2022) describes the data set in more detail.

We reduced the contamination from foreground stars (mainly from the spiral arms along the line of sight towards the Galactic centre and up to a certain extent to the Galactic bar/bulge) by using their significantly different extinction (e.g. Nogueras-Lara et al. 2021a). We excluded them in all fields by applying a colour cut, (HKs) > max(1.3, −0.0233 Ks+1.63), as done in Schultheis et al. (2021) and Nieuwmunster et al. (2024). Additionally, we removed stars with high proper motion uncertainties (> 1 mas/yr). This can lead to an effective magnitude cut at K ~ 10 mag (Fig. 2) because the uncertainties are high for too bright stars (see Sect. 3 in Sormani et al. 2022). We ended up with a total of 1497 stars in the target region for our subsequent analysis.

2.2 Orbital catalogue

We utilised the orbital catalogue obtained by Nieuwmunster et al. (2024), which was derived from the previously described KMOS spectroscopic survey of the NSD (Fritz et al. 2021) and the preliminary proper motions from VIRAC2 (Smith et al. 2024). The orbital catalogue also implemented the previously mentioned colour cut to mitigate potential contamination from foreground sources and applied a strict proper motion uncertainty cut, l,b < 0.6 mas/yr, resulting in a sample of 1130 stars. Nieuwmunster et al. (2024) used Ks luminosity functions to estimate the completeness of the data to be ~ 15% at Ks ~ 15 mag. Thus, the stellar sample is limited, and the fraction of star orbits may not fully represent the overall NSD stellar population.

The orbits of the stars were integrated in a non-axisymmetric potential, combining the potentials corresponding to the Galactic bar/bulge density (Launhardt et al. 2002), the NSD (Sormani et al. 2022), and the nuclear star cluster (Chatzopoulos et al. 2015). They also assumed a distance distribution considering a mean Galactic centre distance of 8.2 kpc and an NSD scale length of ~100 pc (e.g. Gallego-Cano et al. 2020; Sormani et al. 2022; Nogueras-Lara 2022a). The final catalogue comprises several orbital families (i.e. chaotic, z-tube, x-tube, banana, fish, saucer, pretzel, 5:4, and 5:6 orbits), which are present in the NSD region. Examples of each of the orbital families are shown in Appendix A of Nieuwmunster et al. (2024).

Figure 2 shows a colour-magnitude diagram of Ks versus HKs of the observed fields. The parallelogram-shaped area is due to the selection criteria for the target selection of the KMOS survey. The targets were chosen to have an extinction-free apparent magnitude in the range K0 ∈ [7.0,9.5] mag (for further details, see Fritz et al. 2021).

3 Metallicity distribution and proper motions

To investigate the proper motion distribution of stars with different metallicities, we first analysed the metallicity distribution of the target stars by repeating the analysis in Schultheis et al. (2021). Namely, we applied a Gaussian mixture model (GMM; Pedregosa et al. 2011) to obtain the number of Gaussians that best reproduce the metallicity distribution. We distinguished between single-Gaussian and two-Gaussian models by applying the Bayesian information criterion (Schwarz 1978) and the Akaike information criterion (Akaike 1974). We found that the data is best represented by a two-Gaussian model. Figure 3 shows the obtained results. To estimate the mean metallicity and the associated uncertainty for each component, we resorted to Monte Carlo (MC) simulations and generated 1000 MC samples by randomly varying the metallicity of each star, assuming Gaussian uncertainties. We performed the GMM analysis on each of the 1000 samples, assuming a two-Gaussian distribution, and computed the final values along with their associated uncertainties as the mean and standard deviation. The results, shown in Table 1 are in perfect agreement with those reported by Schultheis et al. (2021).

We also employed the GMM approach to calculate the likelihood of each star being associated with each of the Gaussian distributions. To analyse the proper motion distribution of stars belonging to each metallicity component, we assumed that a star belongs to either the metal-rich or the metal-poor component if its likelihood of being associated with that component is larger than 70%. Given that the stellar metallicity of the KMOS catalogue was calibrated using empirical spectra with [M/H] < 0.6 dex, we excluded stars with higher metallicities from the subsequent analysis. Figure 3 shows the final selection of stars. We ended up with 1249 metal-rich and 106 metal-poor stars.

We built histograms for the proper motion component parallel (µl) and perpendicular (µb) to the Galactic plane. To estimate the mean and the standard deviation of each distribution and their corresponding uncertainties, we applied a 3σ clipping algorithm to remove potential outliers, and we used a bootstrap resampling method with 5000 iterations. Table 1 and Fig. 4 show the obtained results. To assess potential sources of systematic uncertainties, we repeated the calculation, considering stars with uncertainties below 1.2 mas/yr and 0.8 mas/yr as well as by adjusting the cut to remove the foreground population (i.e. HKs ~ 1.1 mag). We also repeated the calculation using all stars within the colour cut HKs ≳ 1.3 mag and with proper motion uncertainties below 1 mas/yr, assuming that stars with [M/H] < −0.25 dex are metal poor, while those with −0.25 dex< [M/H] < 0.6 dex are metal-rich. Finally, we explored the effect of the effective cut at Ks ~ 10 mag when considering stars with proper motion uncertainties below 1 mas/yr (see Fig. 2). For this, we defined a new selection cut parallel to the original one (following the reddening vector, Fritz et al. 2011) in order to exclude the stars with Ks > 10 mag, keeping only the ones below the black dashed line in Fig. 2. In all cases, the results were consistent within the uncertainties.

We concluded that although the mean proper motions are similar for both stellar populations, there is a significant difference in the velocity dispersion of the µb component. Namely, the metal-poor component exhibits a µb standard deviation that is approximately a factor of two larger than that of the metal-rich group. Regarding the µl component, we found that the velocity dispersion is slightly larger for the metal-poor stellar population, although there is not a clearly significant difference. This is not surprising because the proper motion component parallel to the Galactic plane is much more affected by the extinction along the line of sight and because of the fact that the NSD is a rotating structure. Therefore, the presence of stars moving eastwards and westwards depends on the completeness of the data and will be different for the different regions analysed (e.g. Schönrich et al. 2015; Martínez-Arranz et al. 2022; Nogueras-Lara 2022a and Figs. A1–A4 in Sormani et al. 2022, which show the µl and µb decomposition between the Galactic bar and the NSD).

We also explored the proper motion orientation, α = atan(µb/µl), of metal-rich and metal-poor stars. We subtracted the proper motion velocity of Sgr A* (µl = −6.396 ± 0.071 mas/yr and µb = −0.239 ± 0.045 mas/yr, see Gordon et al. 2023) to obtain Galactic proper motions relative to the Galactic centre. To ensure α ∈ [−180,180] degrees, we adjusted the output based on the orientation of the proper motions and adopted the origin pointing east. Figure 5 shows the obtained results. We computed the mean orientation values, their standard deviations, and the associated uncertainties by applying a bootstrap resampling method with 5000 iterations. We obtained αrich=5±1°,σαrich=30±1°${\sigma _{{\alpha _{rich}}}} = 30 \pm 1^\circ $, and αpoor = 3 + 5°, σαpoor=48±5°${\sigma _{{\alpha _{poor}}}} = 48 \pm 5^\circ $, for the metal-rich and the metal-poor component, respectively. Although the orientation values for both populations are compatible, the metal-poor stellar population shows a significantly larger standard deviation.

thumbnail Fig. 2

Colour-magnitude diagram Ks versus HKs. The grey dots represent all the stars in the KMOS catalogue in the analysed region after excluding the foreground population. Stars with proper motions below 1 mas/yr and stars with orbits are overplotted with different colours and symbols, as indicated in the legend. The black dashed line corresponds to the cut applied in Sects. 3, 4.2 and 5 to analyse the effect of missing bright stars at Ks ~ 10 mag when applying the proper motion quality cut.

Table 1

Characterisation of the two-Gaussian decomposition of the NSD metallicity distribution.

thumbnail Fig. 3

Two-Gaussian decomposition of the NSD metallicity distribution. The orange line shows the result of the GMM analysis, whereas the blue dashed lines depict each of the Gaussian components. Red and blue bars correspond to stars with a more than 70% probability of belonging to each of the components.

thumbnail Fig. 4

Proper motion distribution of metal-rich and metal-poor stars in the NSD following the GMM analysis (see Sect. 3).

thumbnail Fig. 5

Proper motion orientation of metal-rich and metal-poor stars. The uncertainties correspond to Poisson errors (i.e. #stars/bin$\sqrt {{\rm{\# stars/bin}}} $) scaled to the normalised number of stars in each bin.

thumbnail Fig. 6

Relative fraction of stars with z-tube and chaotic/box orbits for the metal-rich and metal-poor stellar populations defined in Sect. 3. The label “other” corresponds to stars following banana and x-tube orbits (see Nieuwmunster et al. (2024) for further details). The uncertainty bars for the metal-rich population are smaller to the red circles and are therefore not visible in the figure.

4 Metallicity distribution and orbits

4.1 Orbits for each metallicity group

We used the orbital catalogue obtained by Nieuwmunster et al. (2024) to explore whether stars with different metallicities and kinematics undergo different orbits. The authors used a statistical approach to obtain the orbits, and it may not be valid on an individual star basis. However, this does not pose any problem because our objective is to understand the mean behaviour of stars with different metallicities. Nieuwmunster et al. (2024) found that z-tube orbits (including all orbital families with a circulation around the z-axis, such as z-tube, fish, saucer, pretzel, 5:4, and 5:6 orbits) dominate the orbital families, comprising approximately 70% of the stars in the catalogue, and they are directly associated with the NSD. Conversely, chaotic/box orbits constitute about 25% of the sample and are likely related to Galactic bar orbits.

We analysed the fraction of stars with these orbits belonging to the previously defined metal-rich and metal-poor stellar populations. We obtained that the orbits are different for the different metallicity components, as shown in Fig. 6. The uncertainties were estimated based on an assumed membership probability 60% and 80% for each of the metallicity groups, while the mean values correspond to the 70% probability. We concluded that z-tube orbits dominate the metal-rich component (~70% of the metal-rich stars belong to this orbital family) and are less frequent in the metal-poor stellar population, in which the chaotic/box orbits dominate and correspond to ~50% of the stars with defined orbits in the metal-poor component.

thumbnail Fig. 7

Relative fraction of stars with z-tube and chaotic orbits. The contours show the 50 and 75% levels of the NSD self-consistent model obtained by Sormani et al. (2022).

4.2 Metallicity for z-tube and chaotic/box orbits

We also explored the metallicity distribution of all the stars in the catalogue of orbits with z-tube and chaotic/box orbits. To begin, we analysed the spatial distribution of z-tube and chaotic/box orbits. Figure 7 shows the obtained results. We found that the relative fraction of stars undergoing chaotic orbits is more common towards the edge of the NSD, as was also discussed by Nieuwmunster et al. (2024), supporting their relation with potential contaminants from the bar. In contrast, stars with z-tube orbits and, in general, stars without chaotic orbits, are more centrally concentrated and seem to better trace the NSD.

We computed the metallicity distribution for stars with z-tube and chaotic/box orbits and observed that they are different, with z-tube stars on average being more metal rich (Fig. 8). We also applied the previously described GMM approach (see Sect. 3) to compare a single- and a two-Gaussian model. In all cases, we found that a two-Gaussian model is preferred. Figure 8 shows the results for each of the orbital families. We also computed the parameters associated with the metal-rich and metal-poor stellar populations by resorting to MC simulations with 1000 iterations, as explained in Sect. 3. The results are presented in Table 2. We examined the effect of potential systematic uncertainties related to the effective cut at Ks ~ 10 mag resulting from the exclusion of stars with proper motion uncertainties above 0.6 mas/yr in the orbital catalogue (Nieuwmunster et al. 2024). For this, we repeated the analysis, applying the cut shown in Fig. 2 as previously explained in Sect. 3. We did not find any significant difference within the uncertainties.

We observed that stars exhibiting chaotic/box orbits, regardless of their metallicity, are generally more metal poor compared to those following z-tube orbits (Fig. 8). The mean metallicity values obtained for both metallicity components with chaotic/box orbits align well with the results from Schultheis et al. (2021) for the inner regions of the Galactic bar/bulge (0.04 and −0.55 dex for the metal-rich and metal-poor components, respectively) based on analysis of APOGEE data (Ahumada et al. 2020; Rojas-Arriagada et al. 2020). Similarly, the metallicity values obtained for stars undergoing z-tube orbits closely match our results for the NSD (see Table 1).

thumbnail Fig. 8

Metallicity distribution of all stars in the orbital catalogue (Nieuwmunster et al. 2024) following z-tube orbits (including all orbital families with a circulation around the z-axis) and chaotic/box orbits. The solid lines show the result of the GMM analysis, whereas the dashed lines correspond to each Gaussian model.

Table 2

Characterisation of the two-Gaussian metallicity decomposition of stars with different orbits.

5 Extinction of stars with different metallicities along the line of sight

The NSD formed ≳8 Gyr ago and appears to be dominated by an old stellar population (>80% of the stellar mass is older than ≳7Gyr, e.g. Nogueras-Lara et al. 2020, 2022; Sanders et al. 2022a; Schödel et al. 2023). Assuming that metal-rich and metal-poor stars have a similar age allowed us to use the HKs colour as a proxy for the extinction towards these stars (e.g. Nogueras-Lara et al. 2018, 2021b). Moreover, there is a statistical correlation between distance along the line of sight and extinction in the NSD (e.g. Nogueras-Lara 2022a; Nogueras-Lara et al. 2023b), which makes it possible to estimate whether the spatial distribution is similar for both stellar populations by analysing their extinction. We would therefore expect a significant difference in extinction if one of the components is located in a specific region along the line of sight and thus the stars would not be homogeneously distributed.

We built a colour-magnitude diagram Ks versus HKs to analyse the colour distribution (Fig. 9) of the stars belonging to each stellar population whose proper motions were previously analysed. We computed the mean colour and its associated uncertainty for the metal-rich and the metal-poor stars by applying a bootstrap resampling method with 5000 iterations. Table 1 shows the obtained results. We found that the metal-poor stars show a somewhat smaller HKs colour, which might indicate that they are more abundant close to the NSD edge. Nevertheless, the colours for both stellar populations (and their standard deviations) are similar within the uncertainties, so there is no clear difference in the spatial distribution along the line of sight. We detected metal-poor and metal-rich stars with both large and small extinction, as indicated by the HKs standard deviations, showing that both stellar populations are present all over the NSD. We assessed potential systematic uncertainties as previously explained in Sect. 3. Namely, we varied the proper motion uncertainty cut to 0.8 mas/yr and 1.2 mas/yr, we modified the foreground cut to HKs ~ 1.1 mag, and we applied a more restrictive selection cut, as indicated in Fig. 2. The results were always consistent within the uncertainties, with the only difference being a somewhat smaller extinction for both stellar populations when assuming the more restrictive selection criterion. In any case, this was expected because the new cut removes the most reddened stars.

thumbnail Fig. 9

Colour-magnitude diagram of the metal-rich and metal-poor stars used in Sect. 3. The black arrow shows the direction of the red-dening vector.

6 Galactic longitude analysis of stars with different metallicities

We studied the metallicity distribution of KMOS stars at different longitudes from the Galactic centre. In this way, we investigated whether the two-Gaussian metallicity profile is homogeneously distributed across the field and whether metal-rich and metal-poor stars have different spatial distributions.

To avoid mixing stellar populations at significantly different latitudes, whose metallicity might also vary vertically, we restricted our analysis to stars within a latitude of ±20pc from the centre of the NSD (Sgr A*), as shown in Fig. 10. For our analysis, we defined five regions with different longitudes, as colour-coded in the upper panel of Fig. 10.

We applied the GMM approach described in Sect. 3 to check whether the metallicity distribution of each of the defined regions is best represented by one or two Gaussian models. We found that a two-Gaussian distribution is preferred for all the cases. To estimate the mean metallicity and the associated uncertainty for each component and region, we followed the previously described method and employed MC simulations, building 1000 samples, randomly varying the metallicity of each star, and assuming Gaussian uncertainties. Figure A.1 and Table A.1 show the obtained results. We measured a metallicity gradient with increasing metallicities towards the centre of the NSD that is present in both the metal-rich and metal-poor components. Figure 10 shows the detected gradient.

To check whether the metallicity gradients are real and to estimate their values, we employed MC simulations with 1000 iterations. In each iteration, we varied the metallicity of each point, assuming Gaussian uncertainties. We assumed a Galactic centre distance of 8.25 kpc (Gravity Collaboration 2020) and computed the gradients for the negative and positive longitudes, separately for the metal-rich and metal-poor populations. We estimated the final gradient as the mean value of the gradients obtained for the positive and negative longitudes. Our results indicate real gradients with increasing metallicity towards the centre of the NSD for both stellar populations. We ended up with values of (3.4 ± 1.7) ⋅ 10−4dex/pc for the metal-rich population and (1.5 ± 0.5) ⋅ 10−3 dex/pc for the metal-poor population. These gradients are significant at approximately the 2σ and 3σ levels, respectively. However, the metal-poor trend might be influenced by the implicit assumption of the model and the higher metallicity uncertainties of the metal-poor stars (see Fig. 10 in Fritz et al. 2021). In any case, our result aligns with previous work finding a global metallicity gradient (without distinguishing between metal-rich and metal-poor stars) towards the innermost regions of the NSD (Feldmeier-Krause 2022).

We also found that the relative weight of the components does not change spatially within the uncertainties. The relative weight of the metal-rich component is slightly larger than the one obtained when combining all the stars without defining different regions at different latitudes. This is probably due to the metallicity gradient that shifted the mean metallicity of the metal-rich component towards lower metallicities when combining all the stars. Therefore, the metal-rich component encompasses more stars, and its relative weight is somewhat higher.

We also repeated the previous analysis on proper motions and HKs colour for the metal-rich and metal-poor stellar populations in each region. We kept stars whose likelihood of belonging to either group is larger than 70%, restricted the sample to targets whose proper motion uncertainty is below 1 mas/yr, and limited the selection to stars with [M/H] < 0.6 dex, as previously explained. Figure A.2 and Table A.1 show the results. In spite of the relatively low number of stars for the metal-poor population (~ 15–25 stars for each spatial group), we found good agreement with the results obtained when analysing all the stars in Sect. 3. Although there is not a significant difference between the mean proper motions for each stellar population (only in some cases we found a difference larger than 1 σ), we detected that the velocity dispersion of the µb component of the metal-poor stellar population is ≳2 σ larger than that of the metal-rich stellar population in nearly all cases, except at l = 0°, where it is ~1σ larger. This supports the metal-poor component being kinematically hotter than the metal-rich one regardless of the spatial distribution. The bottom panel of Fig. 10 shows the obtained results.

Finally, the colour analysis is also consistent with the previous results (Sect. 5). Thus, the difference between HKs for the metal-rich and the metal-poor stellar population is always ≲ 1 σ.

A similar analysis along the Galactic latitude, as presented in this section, is not feasible due to the low number of stars observed at higher NSD latitudes (see Fig. 1). This limitation is particularly true for metal-poor stars in the innermost fields.

thumbnail Fig. 10

Line-of-sight analysis of the NSD metallicity and kinematics (the corresponding values are included in Table. A.1). Top panel: Target KMOS stars divided into five colour-coded regions at different longitudes. The grey star indicates the position of Sgr A*. Middle panels: mean metallicity obtained for the metal-rich (top) and the metal-poor (bottom) components of the two-Gaussian distribution. We plot metal-rich and metal-poor stars separately to highlight the detected gradients more effectively. Bottom panel: velocity dispersion of µb for the metal-rich and metal-poor components.

7 Discussion

7.1 Correlation between metallicity, velocity dispersion, and orbits

We distinguished between metal-rich and metal-poor stars by applying a GMM approach and obtained that a two-Gaussian model best reproduces the data, as found in previous work (e.g. Schultheis et al. 2021; Nogueras-Lara 2022b; Nogueras-Lara et al. 2023a). In principle, the presence of these two Gaussian components is just a mathematical description of the underlying metallicity distribution and is not necessarily related with the origin of metal-poor and metal-rich stars in the NSD. Actually, the metallicity distribution might be also explained by chemical evolution models (e.g. Grieco et al. 2015; Friske & Schönrich 2023), and thus the presence of metal-rich and metal-poor stars would be a consequence of the natural star-formation process in the NSD. However, the obtained correlation between metal-poor stars, high velocity dispersion, and chaotic/box orbits indicates a different origin for a significant fraction of metal-poor stars in comparison to the metal-rich ones. Our results suggest that the metal-poor stars are related to the innermost regions of the Galactic bar/bulge, which overlaps with the NSD. In this way, the obtained fraction of stars following chaotic orbits in Fig. 7 aligns with the estimated ~20–30% (decreasing for larger NSD radii) of contaminants from the Galactic bar/bulge reported in Sormani et al. (2022). Moreover, the metallicity distribution obtained for stars following chaotic orbits (likely associated with the Galactic bar/bulge) matches previous results for the innermost regions of the Galactic bar/bulge (e.g. Schultheis et al. 2021). Additionally, the higher velocity dispersion found for the metal-poor stellar population is in agreement with previous µb velocity dispersion measurements ~2.5–3 mas/yr obtained for the Galactic bar/bulge in an inner Galactic region far away from any NSD contamination (l = 0° and −4.5° < b < −2.5°, e.g. Kozłowski et al. 2006; Rattenbury et al. 2007; Sanders et al. 2019).

Another possibility is that a fraction of the metal-poor stars might be remnants of dispersed, young gas-rich stellar clusters within a 1.5 kpc radius of the Galactic centre, formed at z>4. These clusters have been suggested as precursors to the old metal-poor stars currently observed at the Galactic centre (e.g. van Donkelaaret al. 2024). The accretion ofstellarclusters, combined with in situ star formation, has been proposed as one of the formation channels of the nuclear star cluster at the centre of the Galaxy (e.g. Tremaine et al. 1975; Capuzzo-Dolcetta 1993; Antonini et al. 2012; Antonini 2013; Perets & Mastrobuono-Battisti 2014; Gnedin et al. 2014; Arca-Sedda et al. 2015, 2020; Tsatsi et al. 2017). Thus, van Donkelaar et al. (2024) used hydro-dynamical simulations to show that the infalling clusters can also be accreted by the NSD contributing to its growth.

We also found that the proper motion orientation of both metallicity components shows a different standard deviation, with the distribution of metal-poor stars being more scattered in comparison to the metal-rich component. This is compatible with the more chaotic motion that is expected under the influence of the Galactic bar/bulge. On the contrary, metal-rich stars are kinematically cooler and show orbits tracing the NSD (e.g. Schultheis et al. 2021; Nieuwmunster et al. 2024).

The similar extinction observed for both the metal-poor and metal-rich components suggests a widespread distribution of both types of stars across the region. Hence, there is not a specific region along the line of sight in which either metal-rich or metal-poor stars are concentrated. This observation is also supported by the spatial analysis in Sect. 5, given that we found metal-rich and metal-poor stars distributed all over the field and with similar relative weights within the uncertainties.

7.2 The NSD as a distinct kinematic structure

Recently, Zoccali et al. (2024) analysed the proper motion distribution of the NSD compared to control fields in the innermost region of the Galactic bar/bulge. Their findings did not reveal any clear, distinct kinematical signature of the NSD. This result challenges previous studies that identified the NSD as kinematically different from its surroundings based on line-of-sight velocities (Schönrich et al. 2015) and velocity dispersions (Schultheis et al. 2021).

In this paper we limited our analysis to the NSD without any comparison with the surrounding region. Nevertheless, we found the presence of two kinematically distinct components with different metallicities and orbits. We detected that the metal-poor component is dominated by chaotic/box orbits that are likely associated with the Galactic bar and have a metallicity distribution also in agreement with previous work on the inner bar/bulge (carried out in fields without contamination from the NSD, e.g. Schultheis et al. 2021). On the other hand, the dominant metal-rich component is kinematically cooler, given its lower velocity dispersion, and undergoes z-tube orbits, better tracing the NSD (e.g. Nieuwmunster et al. 2024). The fact that we observed these two different stellar populations with different orbits and kinematics - with one more closely associated with the Galactic bar/bulge - suggests that the NSD is a different kinematic structure, and this agrees with previous work (Schönrich et al. 2015; Schultheis et al. 2021), showing that it is possible to kinematically distinguish the NSD. Additionally, having a metal-rich and kinematically cool NSD different from the surrounding Galactic bar/bulge also agrees with the results of previous works on external Milky Way-like galaxies, where these structures have been identified via their different kinematics and metallicities (e.g. Gadotti et al. 2019; Bittner et al. 2020; Gadotti et al. 2020).

We would like to emphasise that Zoccali et al. (2024) do not exclude the possibility of a kinematically different NSD in comparison to the surrounding Galactic bar/bulge. They only state that the region corresponding to the NSD is dominated by stars moving eastwards, which is compatible with the CMZ hiding the counter-rotating stars, but this does not mean that the NSD is kinematically different. Nevertheless, analysing the proper motion component parallel to the Galactic plane is very challenging due to the extreme extinction that dominates this region. An in-detail analysis of velocity dispersions, especially the proper motion component perpendicular to the Galactic plane, including fields from the innermost Galactic bar/bulge, is likely to confirm that the NSD is kinematically distinct.

7.3 The metallicity gradient and the NSD formation

Stellar population analyses along and across the line of sight towards the NSD have suggested the presence of radial age and metallicity gradients towards the centre of the NSD (e.g. Nogueras-Lara et al. 2020, 2022; Nogueras-Lara et al. 2023b,a; Feldmeier-Krause 2022; Nogueras-Lara 2024). In this way, the innermost regions of the NSD contain a predominantly old (≳7 Gyr) and metal-rich stellar population, while a significant intermediate-age stellar population (~40% of the stellar mass with ages between 2-7 Gyr) is found towards the outer edge of the NSD, where the metallicity is lower. In contrast, Sanders et al. (2024) analysed the period-age relation of Mira variables in the NSD and found only a weak indication of inside-out formation. They attributed this to potential dynamical mixing and/or the smoothing effect of their spline model, which might hide a significant population of Mira variables with ages around 5–6 Gyr (periods ~400 days). These variables could potentially trace the gradient, but they have large uncertainties.

We observed that the metal-rich and metal-poor components of the NSD are uniformly distributed across the field, and both exhibit a metallicity gradient with a higher metallicity towards the centre of the NSD. These metallicity gradients correspond to an increase of (3.4 ± 1.7) ≳ 10−4 dex/pc and (1.5 ± 0.5) ⋅ 10−3 dex/pc for the metal-rich and metal-poor stellar populations, respectively. To compare our results with the metallicity gradient obtained by Feldmeier-Krause (2022), which does not distinguish between metal-rich and metal-poor stars, we estimated the metallicity increase per parsec using the metallicity values at the edge and centre of a region located ±19 pc in latitude from the centre of the NSD, as specified in their Fig. 11. We obtained a metallicity gradient of approximately 5 × 10−4 dex/pc, which roughly corresponds to the mean of the gradients we found for the metal-rich and metal-poor stars.

The detected gradients align with the currently accepted inside-out formation scenario of the NSD through bar-driven mechanisms (e.g. Bittner et al. 2020; Nogueras-Lara et al. 2023b). Similar gradients have been found in Milky Way-like galaxies analysed in the TIMER survey (e.g. Gadotti et al. 2019; Bittner et al. 2020) and can be attributed to the growth of the galaxy bar that funnels gas from different galactic radii towards the centre as it evolves (the more internal, the more metal rich). Furthermore, in evolved systems, it is anticipated that metallic-ity will rise towards the centre due to the enhanced difficulty for stellar feedback processes to expel gas as it becomes more tightly bound within the galaxy’s potential well. This is compounded by the potential rapid increase in metallicity in star-forming systems.

Interestingly, the metallicity gradients that we found when analysing different NSD lines of sight align well with the previous results from Nogueras-Lara et al. (2023a), who also detected metallicity gradients for metal-rich and metal-poor stars along the line of sight towards the nuclear star cluster (which also contains stars from the NSD that dominate the stellar population in front of the nuclear star cluster). These similarities may suggest a common formation scenario for both the nuclear star cluster and the disc, as discussed by Nogueras-Lara et al. (2023a).

Additionally, the metal-rich stellar population follows a rotation pattern similar to that of the gas in the CMZ (Schultheis et al. 2021). Thus, it is likely that the metal-rich stars formed in situ from gas funnelled towards the Galactic centre via bar-driven processes, contributing to the inside-out growth of the NSD. This also justifies why these stars undergo z-tube orbits and have a lower velocity dispersion in comparison to the metal-poor component. Conversely, the metallicity gradient observed for metal-poor stars, which are kinematically hotter and follow Galactic bar-related orbits, may be caused by the combination of Galactic bar interlopers and/or stars accreted from clusters within a radius of 1.5 kpc at z > 4 (van Donkelaar et al. 2024).

8 Conclusion

In this paper, we have investigated the proper motion and orbit distribution of stars with different metallicities in the NSD. We obtained that metal-rich stars show a lower velocity dispersion in comparison to metal-poor ones and that they are dominated by stars following z-tube orbits that are representative of the NSD (Nieuwmunster et al. 2024). Conversely, metal-poor stars are kinematically hotter, and a significant fraction of them undergo chaotic/bar orbits that are likely related to the Galactic bar. We also found that the proper motion orientation of metal-rich stars shows a sharper profile parallel to the Galactic plane in comparison with stars that are metal poor.

We used the statistical correlation between distance and extinction in the NSD (Nogueras-Lara 2022a) to explore whether stars with different metallicities show a special location along the line of sight. We concluded that the average near-infrared colour of both stellar populations is similar within the uncertainties, and there is no significant difference in extinction. Therefore, we did not detect any preferred location for stars with different metallicities, and they seem to be present all over the observed regions.

Finally, we analysed the metallicity distribution of stars at different longitudes and found that metal-rich and metal-poor stars are present throughout the NSD. We also detected the presence of metallicity gradients for both stellar populations. These gradients have increasing metallicities towards the central regions of the NSD. Our results support an inside-out formation of the nuclear stellar disc, as inferred from studies of external Milky Way-like galaxies (e.g. Gadotti et al. 2019; Bittner et al. 2020). We interpret the gradient in the metal-rich population as being a consequence of in situ formation from gas funnelled towards the Galactic centre through bar-driven mechanisms (e.g. Sormani & Barnes 2019). This justifies the dominant z-tube orbits and the lower velocity dispersion and kinematics that align with the rotation velocity of gas from the CMZ (Schultheis et al. 2021).

Acknowledgements

This work is based on observations made with ESO Telescopes at the La Silla Paranal Observatory under program ID 0101.B-0354. F.N.-L. gratefully acknowledges the sponsorship provided by the European Southern Observatory through a research fellowship. F.F. is supported by a UKRI Future Leaders Fellowship (grant no. MR/X033740/1). M.C.S. acknowledges financial support from the European Research Council under the ERC Starting Grant “GalFlow” (grant 101116226). B.T. acknowledges the financial support from the Wenner-Gren Foundation (WGF2022-0041).

Appendix A Parameters of the spatial analysis of the two-Gaussian metallicity distribution

thumbnail Fig. A.1

Two-Gaussian decomposition of the NSD metallicity distribution for different longitudes. The orange line shows the result of the GMM analysis, whereas the blue dashed lines depict each of the Gaussian components.

Table A.1

Characterisation of the two-Gaussian decomposition of the NSD metallicity distribution at different longitudes.

thumbnail Fig. A.2

Proper motion distribution of metal-rich and metal-poor stars in the NSD for different longitudes following the GMM analysis in Sect. 6.

References

  1. Abdurro’uf, Accetta, K., Aerts, C., et al. 2022, ApJS, 259, 35 [NASA ADS] [CrossRef] [Google Scholar]
  2. Ahumada, R., Allende Prieto, C., Almeida, A., et al. 2020, ApJS, 249, 3 [NASA ADS] [CrossRef] [Google Scholar]
  3. Akaike, H. 1974, Automat. Contrl. IEEE Trans., 19, 716 [NASA ADS] [CrossRef] [Google Scholar]
  4. Alonso-García, J., Minniti, D., Catelan, M., et al. 2017, ApJ, 849, L13 [Google Scholar]
  5. Antonini, F. 2013, ApJ, 763, 62 [Google Scholar]
  6. Antonini, F., Capuzzo-Dolcetta, R., Mastrobuono-Battisti, A., & Merritt, D. 2012, ApJ, 750, 111 [NASA ADS] [CrossRef] [Google Scholar]
  7. Arca-Sedda, M., Capuzzo-Dolcetta, R., Antonini, F., & Seth, A. 2015, ApJ, 806, 220 [NASA ADS] [CrossRef] [Google Scholar]
  8. Arca Sedda, M., Gualandris, A., Do, T., et al. 2020, ApJ, 901, L29 [Google Scholar]
  9. Bittner, A., Sánchez-Blázquez, P., Gadotti, D. A., et al. 2020, A&A, 643, A65 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  10. Capuzzo-Dolcetta, R. 1993, ApJ, 415, 616 [NASA ADS] [CrossRef] [Google Scholar]
  11. Chatzopoulos, S., Fritz, T. K., Gerhard, O., et al. 2015, MNRAS, 447, 948 [Google Scholar]
  12. Crocker, R. M., Jones, D. I., Aharonian, F., et al. 2011, MNRAS, 413, 763 [NASA ADS] [CrossRef] [Google Scholar]
  13. Do, T., Hees, A., Ghez, A., et al. 2019, Science, 365, 664 [Google Scholar]
  14. Dong, H., Wang, Q. D., Cotera, A., et al. 2011, MNRAS, 417, 114 [Google Scholar]
  15. Feldmeier-Krause, A. 2022, MNRAS, 513, 5920 [Google Scholar]
  16. Friske, J. K. S., & Schönrich, R. 2023, arXiv e-prints [arXiv:2310.01899] [Google Scholar]
  17. Fritz, T. K., Gillessen, S., Dodds-Eden, K., et al. 2011, ApJ, 737, 73 [Google Scholar]
  18. Fritz, T. K., Patrick, L. R., Feldmeier-Krause, A., et al. 2021, A&A, 649, A83 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  19. Gadotti, D. A., Sánchez-Blázquez, P., Falcón-Barroso, J., et al. 2019, MNRAS, 482, 506 [Google Scholar]
  20. Gadotti, D. A., Bittner, A., Falcón-Barroso, J., et al. 2020, A&A, 643, A14 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  21. Gallego-Cano, E., Schödel, R., Nogueras-Lara, F., et al. 2020, A&A, 634, A71 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  22. Gnedin, O. Y., Ostriker, J. P., & Tremaine, S. 2014, ApJ, 785, 71 [Google Scholar]
  23. Gordon, D., de Witt, A., & Jacobs, C. S. 2023, AJ, 165, 49 [NASA ADS] [CrossRef] [Google Scholar]
  24. Gravity Collaboration (Abuter, R., et al.) 2020, A&A, 636, L5 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  25. Grieco, V., Matteucci, F., Ryde, N., Schultheis, M., & Uttenthaler, S. 2015, MNRAS, 450, 2094 [NASA ADS] [CrossRef] [Google Scholar]
  26. Henshaw, J. D., Barnes, A. T., Battersby, C., et al. 2023, ASP Conf. Ser., 534, 83 [NASA ADS] [Google Scholar]
  27. Kozłowski, S., Wozniak, P. R., Mao, S., et al. 2006, MNRAS, 370, 435 [CrossRef] [Google Scholar]
  28. Launhardt, R., Zylka, R., & Mezger, P. G. 2002, A&A, 384, 112 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  29. Lindqvist, M., Habing, H. J., & Winnberg, A. 1992, A&A, 259, 118 [NASA ADS] [Google Scholar]
  30. Majewski, S. R., Schiavon, R. P., Frinchaboy, P. M., et al. 2017, AJ, 154, 94 [NASA ADS] [CrossRef] [Google Scholar]
  31. Martínez-Arranz, Á., Schödel, R., Nogueras-Lara, F., & Shahzamanian, B. 2022, A&A, 660, L3 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  32. Matsunaga, N., Kawadu, T., Nishiyama, S., et al. 2011, Nature, 477, 188 [Google Scholar]
  33. Matsunaga, N., Feast, M. W., Kawadu, T., et al. 2013, MNRAS, 429, 385 [CrossRef] [Google Scholar]
  34. Mauerhan, J. C., Cotera, A., Dong, H., et al. 2010, ApJ, 725, 188 [Google Scholar]
  35. Mezger, P. G., Duschl, W. J., & Zylka, R. 1996, A&A Rev., 7, 289 [NASA ADS] [CrossRef] [Google Scholar]
  36. Minniti, D., Lucas, P. W., Emerson, J. P., et al. 2010, New A, 15, 433 [Google Scholar]
  37. Nagayama, T., Nagashima, C., Nakajima, Y., et al. 2003, Proc. SPIE, 4841, 459 [NASA ADS] [CrossRef] [Google Scholar]
  38. Nieuwmunster, N., Schultheis, M., Sormani, M., et al. 2024, A&A, 685, A93 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  39. Nishiyama, S., Nagata, T., Kusakabe, N., et al. 2006, ApJ, 638, 839 [Google Scholar]
  40. Nishiyama, S., Nagata, T., Tamura, M., et al. 2008, ApJ, 680, 1174 [Google Scholar]
  41. Nishiyama, S., Yasui, K., Nagata, T., et al. 2013, ApJ, 769, L28 [Google Scholar]
  42. Nogueras-Lara, F. 2022a, A&A, 668, L8 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  43. Nogueras-Lara, F. 2022b, A&A, 666, A72 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  44. Nogueras-Lara, F. 2024, A&A, 681, L21 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  45. Nogueras-Lara, F., Gallego-Calvente, A. T., Dong, H., et al. 2018, A&A, 610, A83 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  46. Nogueras-Lara, F., Schödel, R., Gallego-Calvente, A. T., et al. 2019, A&A, 631, A20 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  47. Nogueras-Lara, F., Schödel, R., Gallego-Calvente, A. T., et al. 2020, Nat. Astron., 4, 377 [Google Scholar]
  48. Nogueras-Lara, F., Schödel, R., & Neumayer, N. 2021a, A&A, 653, A33 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  49. Nogueras-La?a, F., Schödel, R., & Neumayer, N. 2021b, A&A, 653, A133 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  50. Nogueras-La?a, F., Schödel, R., & Neumayer, N. 2022, Nat. Astron., 6, 1178 [NASA ADS] [CrossRef] [Google Scholar]
  51. Nogueras-Lara, F., Feldmeier-Krause, A., Schödel, R., et al. 2023a, A&A, 680, A75 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  52. Nogueras-Lara, F., Schultheis, M., Najarro, F., et al. 2023b, A&A, 671, L10 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  53. Pedregosa, F., Varoquaux, G., Gramfort, A., et al. 2011, J. Mach. Learn. Res., 12, 2825 [Google Scholar]
  54. Perets, H. B., & Mastrobuono-Battisti, A. 2014, ApJ, 784, L44 [NASA ADS] [CrossRef] [Google Scholar]
  55. Portail, M., Gerhard, O., Wegg, C., & Ness, M. 2017, MNRAS, 465, 1621 [NASA ADS] [CrossRef] [Google Scholar]
  56. Rattenbury, N. J., Mao, S., Debattista, V. P., et al. 2007, MNRAS, 378, 1165 [NASA ADS] [CrossRef] [Google Scholar]
  57. Rojas-Arriagada, A., Zasowski, G., Schultheis, M., et al. 2020, MNRAS, 499, 1037 [Google Scholar]
  58. Sanders, J. L., Smith, L., Evans, N. W., & Lucas, P. 2019, MNRAS, 487, 5188 [CrossRef] [Google Scholar]
  59. Sanders, J. L., Matsunaga, N., Kawata, D., et al. 2022a, MNRAS, 517, 257 [CrossRef] [Google Scholar]
  60. Sanders, J. L., Smith, L., González-Fernández, C., Lucas, P., & Minniti, D. 2022b, MNRAS, 514, 2407 [CrossRef] [Google Scholar]
  61. Sanders, J. L., Kawata, D., Matsunaga, N., et al. 2024, MNRAS, [arXiv:2311.00035] [Google Scholar]
  62. Schödel, R., Nogueras-Lara, F., Hosek, M., et al. 2023, A&A, 672, L8 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  63. Schönrich, R., Aumer, M., & Sale, S. E. 2015, ApJ, 812, L21 [Google Scholar]
  64. Schultheis, M., Sellgren, K., Ramírez, S., et al. 2009, A&A, 495, 157 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  65. Schultheis, M., Fritz, T. K., Nandakumar, G., et al. 2021, A&A, 650, A191 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  66. Schwarz, G. 1978, Annal. Stat., 6, 461 [Google Scholar]
  67. Shahzamanian, B., Schödel, R., Nogueras-Lara, F., et al. 2022, A&A, 662, A11 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  68. Sharples, R., Bender, R., Agudo Berbel, A., et al. 2013, The Messenger, 151, 21 [NASA ADS] [Google Scholar]
  69. Smith, L. C., Lucas, P. W., Koposov, S. E., et al. 2024, MNRAS, submitted [Google Scholar]
  70. Sormani, M. C., & Barnes, A. T. 2019, MNRAS, 484, 1213 [NASA ADS] [CrossRef] [Google Scholar]
  71. Sormani, M. C., Sanders, J. L., Fritz, T. K., et al. 2022, MNRAS, 512, 1857 [CrossRef] [Google Scholar]
  72. Stolovy, S., Ramirez, S., Arendt, R. G., et al. 2006, J. Phys. Conf. Ser., 54, 176 [CrossRef] [Google Scholar]
  73. Surot, F., Valenti, E., Hidalgo, S. L., et al. 2019, A&A, 623, A168 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  74. Thorsbro, B., Ryde, N., Rich, R. M., et al. 2020, ApJ, 894, 26 [NASA ADS] [CrossRef] [Google Scholar]
  75. Tremaine, S. D., Ostriker, J. P., & Spitzer, Jr., L. 1975, ApJ, 196, 407 [NASA ADS] [CrossRef] [Google Scholar]
  76. Tsatsi, A., Mastrobuono-Battisti, A., van de Ven, G., et al. 2017, MNRAS, 464, 3720 [Google Scholar]
  77. Valenti, E., Zoccali, M., Renzini, A., et al. 2013, A&A, 559, A98 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  78. van Donkelaar, F., Mayer, L., Capelo, P. R., et al. 2024, MNRAS, 529, 4104 [NASA ADS] [CrossRef] [Google Scholar]
  79. Wegg, C., & Gerhard, O. 2013, MNRAS, 435, 1874 [Google Scholar]
  80. Wegg, C., Gerhard, O., & Portail, M. 2015, MNRAS, 450, 4050 [NASA ADS] [CrossRef] [Google Scholar]
  81. Zoccali, M., Renzini, A., Ortolani, S., et al. 2003, A&A, 399, 931 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  82. Zoccali, M., Rojas-Arriagada, A., Valenti, E., et al. 2024, A&A, 684, A214 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]

All Tables

Table 1

Characterisation of the two-Gaussian decomposition of the NSD metallicity distribution.

Table 2

Characterisation of the two-Gaussian metallicity decomposition of stars with different orbits.

Table A.1

Characterisation of the two-Gaussian decomposition of the NSD metallicity distribution at different longitudes.

All Figures

thumbnail Fig. 1

Scheme of the observed region overplotted on a Spitzer/IRAC image at 8 µm. (Stolovy et al. 2006). This wavelength highlights the bright and dusty clouds that outline the NSD. a) KMOS fields in which the colour scale indicates the star counts. b) Number of stars per field with computed orbits (Nieuwmunster et al. 2024). The two northernmost fields do not contain any stars due to the selection criteria used to create the orbital catalogue (|l| < 1.5° and |b| < 0.3°; for further details, see Nieuwmunster et al. 2024). The dotted contours show the 25, 50, and 75% levels of the NSD self-consistent model obtained by Sormani et al. (2022).

In the text
thumbnail Fig. 2

Colour-magnitude diagram Ks versus HKs. The grey dots represent all the stars in the KMOS catalogue in the analysed region after excluding the foreground population. Stars with proper motions below 1 mas/yr and stars with orbits are overplotted with different colours and symbols, as indicated in the legend. The black dashed line corresponds to the cut applied in Sects. 3, 4.2 and 5 to analyse the effect of missing bright stars at Ks ~ 10 mag when applying the proper motion quality cut.

In the text
thumbnail Fig. 3

Two-Gaussian decomposition of the NSD metallicity distribution. The orange line shows the result of the GMM analysis, whereas the blue dashed lines depict each of the Gaussian components. Red and blue bars correspond to stars with a more than 70% probability of belonging to each of the components.

In the text
thumbnail Fig. 4

Proper motion distribution of metal-rich and metal-poor stars in the NSD following the GMM analysis (see Sect. 3).

In the text
thumbnail Fig. 5

Proper motion orientation of metal-rich and metal-poor stars. The uncertainties correspond to Poisson errors (i.e. #stars/bin$\sqrt {{\rm{\# stars/bin}}} $) scaled to the normalised number of stars in each bin.

In the text
thumbnail Fig. 6

Relative fraction of stars with z-tube and chaotic/box orbits for the metal-rich and metal-poor stellar populations defined in Sect. 3. The label “other” corresponds to stars following banana and x-tube orbits (see Nieuwmunster et al. (2024) for further details). The uncertainty bars for the metal-rich population are smaller to the red circles and are therefore not visible in the figure.

In the text
thumbnail Fig. 7

Relative fraction of stars with z-tube and chaotic orbits. The contours show the 50 and 75% levels of the NSD self-consistent model obtained by Sormani et al. (2022).

In the text
thumbnail Fig. 8

Metallicity distribution of all stars in the orbital catalogue (Nieuwmunster et al. 2024) following z-tube orbits (including all orbital families with a circulation around the z-axis) and chaotic/box orbits. The solid lines show the result of the GMM analysis, whereas the dashed lines correspond to each Gaussian model.

In the text
thumbnail Fig. 9

Colour-magnitude diagram of the metal-rich and metal-poor stars used in Sect. 3. The black arrow shows the direction of the red-dening vector.

In the text
thumbnail Fig. 10

Line-of-sight analysis of the NSD metallicity and kinematics (the corresponding values are included in Table. A.1). Top panel: Target KMOS stars divided into five colour-coded regions at different longitudes. The grey star indicates the position of Sgr A*. Middle panels: mean metallicity obtained for the metal-rich (top) and the metal-poor (bottom) components of the two-Gaussian distribution. We plot metal-rich and metal-poor stars separately to highlight the detected gradients more effectively. Bottom panel: velocity dispersion of µb for the metal-rich and metal-poor components.

In the text
thumbnail Fig. A.1

Two-Gaussian decomposition of the NSD metallicity distribution for different longitudes. The orange line shows the result of the GMM analysis, whereas the blue dashed lines depict each of the Gaussian components.

In the text
thumbnail Fig. A.2

Proper motion distribution of metal-rich and metal-poor stars in the NSD for different longitudes following the GMM analysis in Sect. 6.

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.