Open Access
Issue
A&A
Volume 653, September 2021
Article Number A33
Number of page(s) 15
Section Galactic structure, stellar clusters and populations
DOI https://doi.org/10.1051/0004-6361/202040073
Published online 03 September 2021

© F. Nogueras-Lara et al. 2021

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.

Open Access funding provided by Max Planck Society.

1. Introduction

The position of the Sun in the Galactic disc allows us to perform a detailed analysis of the stellar population in its close vicinity (e.g., Bland-Hawthorn & Gerhard 2016; Ruiz-Lara et al. 2020). Nevertheless, this position significantly hampers the study of the spiral arm structure of the Milky Way and its properties, given the interstellar dust that characterises the low-Galactic-latitude lines of sight (e.g., Chen et al. 2017). However, the characterisation of the spiral arms of the Milky Way is of fundamental interest because it is needed in order to understand our Galaxy in a wider context of galactic morphology, dynamics, and evolution.

Much effort during recent decades has gone into properly characterising the main parameters (number, shape, inter-arm separation) of the spiral arm structure (e.g., Bland-Hawthorn & Gerhard 2016, and references therein). Nevertheless, the whole structure remains unclear (e.g., Momany et al. 2006; Hou & Han 2014; Khoperskov et al. 2020). According to recent studies based on trigonometric parallaxes and proper motions of molecular masers associated with young high-mass stars, the Milky Way appears as a four-arm spiral galaxy with some extra segments and spurs (Reid et al. 2019).

In particular, the analysis of the line of sight towards the Galactic centre (GC) is hampered by the extreme cumulative extinction (with an average value of AV ≳ 30 mag, corresponding to AKs ≳ 2.5 mag, see e.g., Nishiyama et al. 2008; Schödel et al. 2010; Fritz et al. 2011; Nogueras-Lara et al. 2018a, 2020b) due to the dust and gas present in the Galactic disc. Therefore, the analysis of its stellar population is restricted to infrared wavelengths, which suffer least from interstellar reddening. The precise measurements of the Gaia survey (e.g., Gaia Collaboration 2018) within the Galactic disc are limited to distances ≲3 kpc from the Sun, therefore most studies of the Milky Way disc structure used radio observations of masers and/or star forming regions; their distances are measured by means of trigonometric parallaxes (e.g., Reid et al. 2016, 2019; Bland-Hawthorn & Gerhard 2016; Wu et al. 2019). Moreover, infrared surveys such as the COBE/DIRBE survey and the infrared Spitzer/GLIMPSE survey have been used to trace the spiral arm structure (e.g., Drimmel & Spergel 2001; Churchwell et al. 2009).

The analysis of red clump (RC) stars (red giant stars in their helium-core-burning sequence, e.g., Girardi 2016) using near infrared (NIR) photometry is key to analysing the structure and stellar population towards the innermost regions of the galaxy. In this way, NIR data from the 2MASS (Skrutskie et al. 2006) and UKIDSS (Lucas et al. 2008) surveys have been used to trace the spiral structure and to outline the shape of the Galactic bulge and bar (e.g., Cabrera-Lavers et al. 2008; Francis & Anderson 2012; Robin et al. 2012). More recently, the VVV survey (Minniti et al. 2010; Saito et al. 2012) allowed improvement of the study towards the obscured regions of the Galactic bar (e.g., Gonzalez et al. 2011; Minniti et al. 2014), and also permitted the spiral arm structure to be traced beyond the Galactic bulge (e.g., Gonzalez et al. 2018; Saito et al. 2020).

With the present study, we aim to characterise the spiral arms along the line of sight towards the GC, computing their distance and the extinction using NIR photometry from the GALACTICNUCLEUS survey (Nogueras-Lara et al. 2018a, 2019a). It is specially designed to characterise the structure and the stellar population of the nuclear bulge of the Milky Way with unprecedented detail. The high angular resolution of the data and the fact that it is several magnitudes deeper than any existing NIR catalogue for the GC (e.g., Nogueras-Lara et al. 2018a, 2019a) make this survey the best of its kind and the only one that reveals the details of the structure of the Galactic disc towards this extremely crowded and highly extinguished line of sight. We also used the Gaia survey to cross-check our results for the stars belonging to the closest spiral arms.

2. Data

For the analysis presented in this paper, we used the GALACTICNUCLEUS catalogue (Nogueras-Lara et al. 2018a, 2019a). This is a NIR JHKs survey carried out with the HAWK-I instrument (Kissler-Patig et al. 2008) at the ESO VLT unit telescope 4. It is a single-epoch catalogue consisting of 49 pointings in all three bands (JHKs), and includes observations of more than three million sources in the nuclear stellar disc (NSD), the innermost Galactic bulge, and the transition region between the bulge and the NSD. It uses the speckle holography technique (Schödel et al. 2013) to achieve high angular resolution of ∼0.2″, reaching 5 σ detections at J ∼ 22 mag, H ∼ 21 mag, and Ks ∼ 21 mag. The survey supersedes all previous photometric surveys for the GC region by several magnitudes. The uncertainties on the photometry are below 0.05 mag at J ∼ 21 mag, H ∼ 19 mag, and Ks ∼ 18 mag. The zero point (ZP) systematic uncertainty is below 0.04 mag in all three bands.

As we aim to analyse the stellar population between the Earth and the GC, we used five regions from GALACTICNUCLEUS that correspond to different lines of sight given their different Galactic latitudes. In particular, we used the central region, the transition regions east (TE) and west (TW), and the inner bulge regions north (IBN) and south (IBS). Figure 1 shows a scheme of the studied area. Given the relatively low number of stars in the TE, TW, IBS, and IBN regions, we combined the transition regions and the inner bulge ones to end up with central, transition, and inner bulge regions.

thumbnail Fig. 1.

False colour image from the GALACTICNUCLEUS survey (using JHKs bands) showing the regions used in this study. The red arrow marks the position of the supermassive black hole, Sagittarius A* (Sgr A*). The names of the regions are also indicated in the figure: central, transition east (TE), transition west (TW), inner bulge south (IBS), and inner bulge north (IBN). The distances between the regions correspond to real distances, as indicated by the scaling arrows.

The stellar population belonging to the GC is highly reddened (e.g., Schödel et al. 2010; Nogueras-Lara et al. 2018a, 2019a), which allows us to clearly distinguish it from the foreground population belonging to the stellar disc using a simple colour cut (Nogueras-Lara et al. 2020a). Figure 2 shows the colour-magnitude diagram (CMD) Ks versus J − Ks of the central region of the GALACTICNUCLEUS survey. This region allows us to easily differentiate between the stellar disc and the GC because it is the most extinguished region and the one with the largest number of stars in the GALACTICNUCLEUS catalogue. The most prominent feature (around Ks ∼ 15, J − Ks ∼ 5 mag in the upper panel) corresponds to RC stars belonging to the GC (see Fig. 14 in Nogueras-Lara et al. 2019a). Namely, the GC stellar population comprises stars from the nuclear bulge (NB, e.g., Launhardt et al. 2002; Nogueras-Lara et al. 2020a) and the inner Galactic bulge (e.g., Nogueras-Lara et al. 2018b), the stellar population of this latter being less reddened than the one pertaining to the NB (see Fig. 2 for details). Stars located at redder colours correspond to the Galactic disc. In this way, and as an initial assumption that is later corroborated by our analysis, the red dashed line in the figure indicates the approximate colour division between the foreground population and GC stars.

thumbnail Fig. 2.

Colour-magnitude diagram Ks vs. J − Ks from the central region of the GALACTICNUCLEUS survey. The red dashed line indicates the colour cut between the stellar population belonging to the GC: Galactic bulge (GB) + Nuclear bulge (NB), and the Galactic disc (GD). The contours outline the population belonging to the GD.

thumbnail Fig. 3.

Model fitting performance for a given MC sample (left panel). Middle panel: best fit for the MC sample (). Right panel: residuals after subtracting the model from the data. The colour bar indicates the number of stars per bin using a linear scale.

3. Colour-magnitude diagram fitting

Previous work using the GALACTICNUCLEUS data suggested the presence of spiral arms towards the GC by visual inspection of the CMDs (e.g., Figs. 14 and 14 from Nogueras-Lara et al. 2018a, 2019a, respectively). A number of studies based on different spiral arm tracers detect four spiral arms along the line of sight towards the GC (e.g., Hou & Han 2014; Vallée 2016; Reid et al. 2019). In this way, we analysed our data, constructing a simple synthetic model considering the four spiral arms scenario as a fiducial model.

3.1. Models

We created synthetic models for each of the spiral arms assuming the star formation history (SFH) derived by Ruiz-Lara et al. (2020) for the thin disc in a 2-kpc-radius bubble around the Sun. In this way, we combined stellar populations of 13, 6, 2, 1, and 0.1 Gyr and scaled them in agreement with the SFH shown in Fig. 4 of Ruiz-Lara et al. (2020). The metallicity of the models was selected taking into account the age–metallicity relation for the solar neighbourhood, assuming the mean values for different age distributions from Fig. 18 in Feltzing et al. (2001). In particular, we used [M/H]= − 0.07 for all the models except for the oldest one ([M/H]= − 0.28). This metallicity distribution is also in agreement with previous studies of the Galactic disc (Bergemann et al. 2014; Ruiz-Lara et al. 2020).

We used PARSEC models (release v1.2S + COLIBRI S_35, Bressan et al. 2012; Chen et al. 2014, 2015; Tang et al. 2014; Marigo et al. 2017; Pastorelli et al. 2019) and created the synthetic CMDs using the web tool CMD 3.31. We used the CMD H versus J − H for the model fitting because in this way we avoid the Ks band that is more affected by saturation and crowding than the H band. Additionally, using the CMD H versus J − H we can more clearly separate the different components along the line of sight because their mean colour differences due to extinction are significantly larger at J − H than the scatter due to measurement uncertainties (Nogueras-Lara et al. 2019a). The model was generated as follows:

The main parameters of the model are the distance to each spiral arm and its associated extinction (d and AH, respectively). The distance only affects the y-axis of the CMD (the colours are not distance dependent). On the other hand, to compute the extinction in the J band (AJ) given AH, we used the relation AH/AJ = 1.87 ± 0.03 obtained by Nogueras-Lara et al. (2020b).

– We simulated the photometric uncertainties of the data assuming a Gaussian distribution with a standard deviation of 0.05 mag for each of the bands. This value considers the statistical uncertainties and also the possible ZP variations (∼0.04 mag) between the different pointings that were combined to get the final photometric catalogue (Nogueras-Lara et al. 2019a).

– We simulated a background of stars, whose relative contribution is a free parameter of the fit. This background follows an exponential distribution to account for the Galactic disc component. We used a scale length of 2600 pc (Bland-Hawthorn & Gerhard 2016) and placed stars at distances between 100 and 5500 pc from the Earth. The interstellar extinction was assumed to linearly increase from AH = 0 (corresponding to 0 pc) to AH = 2.025 mag. This value was computed using the GC extinction of A1.61 μm = 3.40 mag (at ∼8 kpc) obtained in Nogueras-Lara et al. (2020b), and assuming that it increases linearly with the distance from the Earth to the GC. Given the complexity of the extinction towards the innermost regions of the Galaxy (e.g., Nataf et al. 2016; Alonso-García et al. 2017; Nogueras-Lara et al. 2020b), we used the three-dimensional extinction maps obtained by Schultheis et al. (2014) to assess the assumption of a linear increase in extinction. These authors obtained a steep rise in AKs for different lines of sight near the GC, with a flattening occurring at ∼4 − 6 kpc (see their Fig. 8). The steep behaviour is compatible with the linear increase that we assumed, and it is probably due to the presence of gas and dust in the Galactic disc. On the other hand, the flattening is likely produced by the lower dust and gas content associated to the Galactic bulge in comparison to the Galactic disc (e.g., Englmaier & Gerhard 1999). The distance ∼4 − 6 kpc was obtained using the extinction curve obtained by Nishiyama et al. (2009) (AJ/AKs ∼ 3.02). Using the more recent extinction curve by Nogueras-Lara et al. (2020b) (AJ/AKs ∼ 3.44) would increase the estimated distance where the flattening occurs by ∼1 kpc, matching our assumption of an upper limit for the contribution of background stars at 5500 kpc from the Earth.

– The numbers of stars in each spiral arm were defined as a ratio to the number in the first arm.

– We considered that the observed area increases as the square of the distance, implying that the number of stars is proportional to the square radius.

– We assumed that all the spiral arms have the same width (∼400 pc from mid-arm to dust lane, Vallée 2014). To simulate this, we used a Gaussian with a standard deviation of 400 pc to randomly generate the width effect for each spiral arm.

– The extinction to each spiral arm was a free parameter and depends on each of the models considered.

– After having generated the models, we applied a Gaussian uncertainty distribution to the simulated J and H photometry. This accounts for the possible differential extinction towards the features and also for some possible additional measurement uncertainties. For this, we used Gaussians with a standard deviation of 0.05 mag.

3.2. Fitting procedure

For the model fit, we restricted the data to stars that belong to the foreground population as indicated in Fig. 2. We only used the central region, because the lower extinction towards bulge and transition regions makes it difficult to separate Galactic disc stars from GC stars and because the lower number of stars towards these regions would lead to lower quality fits. Interstellar extinction towards the GC is significantly higher at low latitudes, which facilitates the separation of disc and GC stars (Schödel et al. 2014; Nogueras-Lara et al. 2018b).

We only considered stars with J − H ∈ [0.3, 2.3] and H ∈ [10, 18.3] mag, where the completeness in H band is more than 80% according to Nogueras-Lara et al. (2020b). To estimate the uncertainty of the fits and to significantly decrease the necessary computing time, we implemented a Monte Carlo (MC) simulation creating 500 independent data sets, randomly sampling 5000 of the accepted stars for each of them.

To compare the model with the data, we binned the selected region of the CMD. We used the Python function numpy.histogram_bin_edges (Harris et al. 2020) to compute the bin widths for each axis. We ended up with 0.09 and 0.18 mag for the J − H and H axes, respectively.

We generated a grid of models to fit the MC samples, varying the free parameters (see Table 1) and identifying the best ones via χ2 minimisation. The distance parameters were selected to homogeneously sample the parameter space in steps of 350 pc. The extinction values also cover the space homogeneously in steps of 0.15 mag up to an upper limit of 2.3 mag.

Table 1.

Parameters to create the models.

We used the Poisson maximum likelihood parameter χ2 = 2Σimi − ni + niln(ni/mi), where mi is the number of the stars of the model, ni is the number of stars in the observations, and the subindex i indicates the bin of the CMD (e.g., Mighell 1999; Dolphin 2002; Pfuhl et al. 2011). Our grid approach allowed us to: (1) not have to choose any starting parameters for spiral arm distances and extinctions, (2) explore the whole parameter space avoiding false minima that might appear when using non-linear least-squares approaches, and (3) reduce the computing time to make the computation feasible.

3.3. Properties of the spiral arms

We computed the reduced χ2 (), considering the fixed parameters for each of the MC samples. The minimum values from the MC runs show a Gaussian distribution. We obtained a mean , where the uncertainty corresponds to the standard deviation of the values.

Table 2 shows the best solutions for the parameters fixed averaging over the results obtained for the 500 MC realisations. The uncertainties were obtained considering the maximum between the standard deviation of the parameter distribution for each of the 500 MC samples and the half step of the grid (175 pc and 0.08 mag, for the distance and the extinction, respectively).

Table 2.

Model fitting results.

We addressed the following sources of systematic errors:

– We analysed the CMD Ks vs. H − Ks. We used the synthetic population described in Sect. 3.1 and kept the same distance parameters in steps of 350 pc to sample the whole parameter space. We defined a grid of AKs in agreement with the AH grid used for the analysis of the CMD H versus J − H. In this way, we created the models assuming four extinction values for each spiral arm varying from 0.075 to 1.2 mag in steps of 0.075 mag. We also considered the lower Ks-band completeness (mainly due to crowding) in comparison to the H and J bands used previously (80% completeness at Ks ∼ 16.3 mag, see Nogueras-Lara et al. 2020b). For this, we applied a completeness correction selecting a reference level with a completeness of 50% (Ks = 18.3 mag) and randomly removing stars from levels of completeness above in steps of 1%, as described in Sect. 3.1.2 of Nogueras-Lara et al. (2020b). We selected stars fulfilling Ks ∈ [11.5, 18.3] and H − Ks ∈ [0, 1.05] to account for the significant saturation in Ks for stars brighter than 11.5 mag (Nogueras-Lara et al. 2019a) and also to minimise the contamination from the inner Galactic bulge (see Fig. 2). We then sampled the data to create the 500 MC realisations in order to apply the model fitting. Table 2 shows the results. We concluded that there is no significant variation in the distance parameters within the uncertainties with respect to the results obtained for the J − H case.

– We also analysed the systematic errors due to the chosen bin width. For this, we repeated the fit varying the bin size of the CMD (±25%). The results are presented in Table A.1. We did not observe any significant difference within the uncertainties.

– We tested the influence of the ZP variability on the results. We repeated the fit adding and subtracting the ZP uncertainty (0.04 mag for J and H bands) independently for each band. The results are presented in Table A.2 and agree well with our results.

– To study the influence of the selected models, we repeated the analysis using BaSTI models (Pietrinferni et al. 2004, 2006) to generate the synthetic population described in Sect. 3.1. Because the synthetic population generated using the BaSTI web tool2 only considers stars with masses larger than 0.5 M, we restricted the lower limit of the H cut to 17 mag instead of 18.3 mag to avoid gaps when analysing the closest spiral arms. The results are presented in Table A.3. We did not observe any significant difference within the uncertainties.

– We also tested the possible systematic errors introduced by the chosen SFH. For this, we used a BaSTI (Pietrinferni et al. 2004, 2006) model that considers a stellar population similar to that expected for the local disc according to Rocha-Pinto et al. (2000) as it is implemented in the BaSTI web tool. The selection of the stars was also limited to H = 17 mag given the limitations of the BaSTI models. The results are presented in Table A.3. We did not observe any significant difference within the uncertainties.

– Finally, to check the method, we generated three simulated data sets considering the same SFH but more realistic situations. We varied the relative number of stars between arms, the width of the spiral arms, and the relative contribution of the stellar background. We applied the method described previously, now using 100 MC samples. The results are presented in Table A.4. We find that the extinction of all the features is always determined within the uncertainties in all the cases. On the other hand, we measured some deviation in the results for the distances of the first simulated spiral arm. The determination of the distance is more dependent on the precise shape of the spiral arm features given that it only affects the y-axis (the colour, x-axis, does not depend on distance). Thus, it is more influenced by the scatter of the data, or by a less defined spiral arm feature. In this sense, the determination of the distance for the first spiral arm might be more affected because its number of stars is lower in comparison with the other spiral arms. Nevertheless, the measured deviations are always within 1.5 σ. In the case of real data, we believe that the distance to the first spiral arm is properly determined given the agreement with Gaia (Sect. 4) data and also with previous work. It might be possible that the scatter of the data is overestimated when creating the simulations. We also measured some deviation when computing the influence of the background population. Nevertheless this is expected given our simple approach, which is designed to simply identify the extinction and distance to the spiral arms.

3.4. The role of extinction

The extreme extinction towards the GC, in particular towards the NB (e.g., Nogueras-Lara et al. 2020a), allows us to clearly distinguish the foreground population and trace the spiral arms as shown above. Nevertheless, some contamination from the bulge is expected for the fourth spiral arm, in particular when using the Ks band. We overcame this difficulty when analysing the J and H bands, where the effect of extinction is larger (e.g., Nishiyama et al. 2009), making it easier to distinguish between components with different reddening. Comparing the results with the ones obtained when applying the technique to the CMD Ks versus H − Ks, allows us to assess the consistency between the distance and extinction results. We ended up with a larger distance for the fourth spiral arm when using the CMD Ks versus H − Ks, although this result remains compatible within the uncertainties with the results from the CMD H versus J − H (Table 2). Therefore, we conclude that the obtained parameters are not significantly affected by this contamination.

A similar analysis is not possible at higher latitudes given the lower extinction towards the Galactic bulge at these latitudes. The inner bulge fields analysed in Sect. 4 correspond to an average extinction of AKs ∼ 1.3 mag (averaging over the values obtained by Nogueras-Lara et al. 2018b, for two inner-bulge regions located at ∼0.6° and ∼0.4° to the Galactic north of the Milky Way centre), which is only about half of the total extinction towards the central field (latitude ∼0°). Moreover, the larger size of the central field and the lower density of stars given the projection effects for higher latitudes also affect the observed number of stars from the spiral arms in the transition and the inner bulge regions.

Therefore, the stars in the third and fourth (and to a certain degree even the second) spiral arm from the Earth overlap with the GC stars in the CMDs for the bulge fields. Figure 4 illustrates this effect showing the CMDs H versus J − H and Ks versus H − Ks for the studied regions. We over-plotted PARSEC isochrones of 1 Gyr and [M/H]= − 0.07 according to the parameters of distance and extinction obtained for the J − H and H − Ks data of the central region (Table 2). We observed good agreement between the position of the isochrones and the over-densities in the CMDs, within the uncertainties. Given the lower extinction, the features corresponding to the GC are less extinguished for the transition and the inner bulge regions than for the central one. In this sense, the middle and lower panels of Fig. 4 show how the RC features from the inner bulge or the GC stellar population are contaminated by the presence of stars from the main sequence belonging to the fourth spiral arm. This effect is visible in the CMDs H versus J − H, where the RC feature shows an excess of stars belonging to the isochrone associated to that spiral arm. The zoomed-in region in Fig. 5, corresponding to the inner bulge regions, depicts this effect. In this particular case, the different extinction from the IBN and IBS (Nogueras-Lara et al. 2018b) magnifies the overlap.

thumbnail Fig. 4.

CMDs H vs. J − H and Ks vs. H − Ks for the central, transition, and inner bulge regions of the GALACTICNUCLEUS survey, from the upper to the lower panels, respectively. The colour code refers to stellar densities, using a power-stretch scale (different in each case to stress the spiral arm features). Cyan lines correspond to a PARSEC isochrone of 1 Gyr and [M/H = −0.07] according to the extinction and distance parameters obtained for the spiral arm. The blue and green error bars indicate the 1 σ uncertainties of the position of the isochrones and the photometric ZP, respectively. The zoomed-in region corresponds to the red clump features (e.g., Girardi 2016) from the inner Galactic bulge (‘1’ and ‘2’ associated to different extinctions from the IBS and IBS regions) and the overlap with the fourth spiral arm.

3.5. Discussion

A detailed two-dimensional fit of the experimental CMDs is very complex given the high number of free variables to be considered. Namely, the presence of stars with different ages and/or metallicities, the different contributions from each spiral arm depending on its corresponding distance and/or width, the differential extinction, and the possible contamination from the GC stellar population for the innermost arms contribute to this increase in complexity.

Nevertheless, we did not aim to characterise the stellar population and/or metallicities of the analysed stars, but to compute the distance to the spiral arms and their extinction given the average features in the CMDs. We demonstrate that our method is good enough for our purposes, addressing potential sources of systematic errors and trying different stellar populations and models. Our results agree well with the distances from the recent work by Reid et al. (2019). From Fig. 1 of this latter paper we obtained d1 ∼ 1.4, d2 ∼ 2.7, d3 ∼ 3.7, and d4 ∼ 4.7 kpc, which are fully compatible with our results.

The obtained structures would correspond to the Sagittarius–Carina, the Scutum–Centaurus–OSC, the Norma–Outer, and the 3 kpc arms. As reference values for the distance towards the spiral arms, we used the results obtained for the CMD H versus J − H because of the lower contamination from the GC and also because the shape of the isochrones is more sensitive to changes in the y-axis as can be seen in Fig. 4.

4. Analysis of Gaia sources

To cross-check the obtained distances to the spiral arms, we cross-correlated the regions from the GALACTICNUCLEUS survey with the Gaia DR2 catalogue (Gaia Collaboration 2018). The mean offset between the two catalogues is ΔRA = 0.05 ± 0.02 arcsec and ΔDec = 0.04 ± 0.02 arcsec. We defined a maximum distance of ∼0.15″ to cross-correlate the catalogues. Moreover, we only accepted stars whose distances was measured with more than 3 σ significance, and removed all the stars with distances larger than 10 000 pc. In this way, we avoid erroneous identifications and stars with spurious astrometric solutions (Gaia Collaboration 2018).

The high extinction towards the GC and the optical photometry from Gaia only allowed us to obtain common stars whose colours are restricted to J − H ≲ 1 mag for faint stars (H ≲ 16 mag) and J − H ≲ 1.5 mag for some bright stars around H ≲ 12 mag. To obtain the distances to the detected common stars, we used the Gaia parallaxes taking into account that, in general, they cannot be directly inverted to get the distance to a given source (e.g., Luri et al. 2018; Bailer-Jones et al. 2018). We applied the near-linear trend of distance bias (δp = −0.054) derived by Schönrich et al. (2019). To remove any additional problem related to the distance to individual stars, we computed average distances building histograms of the underlying distributions corresponding to each of the analysed regions (Fig. 5). We ended up with ∼5000, ∼800, and ∼400 stars for the central, transition, and inner bulge regions, respectively.

thumbnail Fig. 5.

Histograms of the distances to the common stars between Gaia and GALACTICNUCLEUS detected in the central, transition, and inner bulge regions (upper, middle, and lower panel, respectively). The solid green lines indicate the two-Gaussian fits, whereas the red and cyan dashed lines show each of the two Gaussian components. The mean values and their errors for each component are indicated for each panel in units of kpc.

Given the high extinction along the line of sight towards the GC, we can only detect stars located at ≲3 kpc from Earth using Gaia data for this line of sight, restricting the analysis to the closest spiral arms. We detected a bimodal distribution that can be fitted well with a two-Gaussian model. We computed the average distance to each of the features as the mean value of each of the Gaussians. The uncertainties were estimated considering the error of the mean. The results and the associated uncertainties are over-plotted in Fig. 5. We find that the two-Gaussian distribution fits the data better than a single Gaussian when selecting stars with distances of < 5000 pc and using the SCIKIT-LEARN Python function GaussianMixture (GMM Pedregosa et al. 2011) to compute the Bayesian information criterion (Schwarz 1978) and the Akaike information criterion (Akaike 1974). We find that the two-Gaussian distribution is preferred. We did not try any more complex models because of the low completeness of the data beyond 3000 pc, and so no more components are expected (nor visually identified) in the distributions.

Our results agree well with the presence of two spiral arms, particularly in the case of the central and the transition regions, where the number of stars is larger than for the inner bulge, and where the distance values obtained are somewhat smaller. The secondary Gaussian feature is significantly smaller than the first one, as expected. This is because only a small fraction of the stars belonging to the second spiral arm were detected in the Gaia survey given the extinction and the larger distance. We computed a mean value for the distance to each of the features averaging over the obtained results and ended up with d1_Gaia = 1.5 ± 0.1 kpc and d2_Gaia = 2.8 ± 0.2 kpc. The uncertainties refer to the standard deviation of the measurements towards the three different regions. The obtained results agree within the uncertainties with the distances obtained for the first and the second spiral arms in the previous section and in the literature (see Sect. 1).

5. Analysis of the colour–colour diagram

Figure 6 shows the colour–colour diagram (CCD) J − H versus H − Ks of the central region of the GALACTICNUCLEUS survey. The highest density region in the CCD (H − Ks > 1.3 mag and J − H > 2.65 mag) corresponds to the GC. These are stars from the NB and the inner Galactic bulge. Excluding these, there are two main over-density features (J − H ∼ 1 mag and J − H ∼ 1.8 mag) that can be identified in the right panel of Fig. 6. According to the analyses carried out in the previous sections, these overdensities correspond to the spiral arm structure along the line of sight from the Earth to the GC. We built histograms of the J − H and H − Ks distributions to further analyse these features (right panels in Fig. 6):

thumbnail Fig. 6.

Left panel: colour–colour diagram J − H vs. H − Ks corresponding to the central region of the GALACTICNUCLEUS survey. The black arrows indicate the over-density features that correspond to the spiral arms detected in Sect. 3. The blue dashed lines mark the stellar population belonging to the GC, where the density is significantly different from the rest. Right panel: zoom into the foreground population from the spiral arms and histograms corresponding to the underlying distributions. The red solid line indicates the result of a four-Gaussian fit (green dashed lines for each individual Gaussian). The black arrows show the mean location of the spiral arms. The red arrow indicates the overlap between the first and the second spiral arms in one single feature in the histogram J − H.

J − H: Given the overlap between the isochrones corresponding to the first and the second spiral arms that we observed in the CMDs H versus J − H (left panels Fig. 4), it is not possible to distinguish the four spiral arms in the histogram corresponding to the J − H distribution. We observed a prominent narrow peak around J − H ∼ 1 mag and a secondary extended peak around J − H ∼ 1.75 mag. According to the synthetic CMD fitting, the first peak is probably the result of the overlap between the first and the second spiral arms. On the other hand, the extended peak might be the result of the overlap between the third and the fourth spiral arms.

H − Ks: There are two clear features that probably correspond to the overlap between the first and second, and third and fourth spiral arms. Given that the overlap between the isochrones of the spiral arms is not as significant as for J − H (see Fig. 4), we tried to fit this histogram using a four-Gaussian model to check whether it is compatible with the structure obtained previously. This is a simplistic approach that only attempts to detect over-densities associated to the spiral arms in the CCD.

We used the stars detected in all three bands (JHKs) within the same magnitude ranges used in Sect. 3 for the analysis of the CMDs H versus J − H and Ks versus H − Ks. We fitted the model using the Python SCIPY routine curve-fit and the initial values according to the expected positions of the spiral arm features. We only imposed that the four Gaussians cannot take negative values for the fit. We obtained that this simple model agrees well with the results from previous sections. We obtained mean values of (H − Ks)1 ∼ 0.2, (H − Ks)2 ∼ 0.4, (H − Ks)3 ∼ 0.8, and (H − Ks)4 ∼ 0.9 mag, where the subindex indicates the number of each spiral arm. The feature corresponding to the fourth spiral arm is smaller than the others. Nevertheless, this is expected given the colour cuts applied to limit possible confusion with the stars belonging to the inner Galactic bulge, whose extinction is close to that from this spiral arm (Nogueras-Lara et al. 2018b).

Therefore, the analysis of the CCD allows us to distinguish between GC and spiral arm stellar populations and agrees with the structure derived previously.

6. Extinction

According to the previous analysis, we find that the CMD Ks versus H − Ks is the best choice to select the reference stars with which to analyse the extinction corresponding to each feature. This is because the stellar population belonging to each spiral arm can be disentangled more easily there than in the CMD using the J band, where the isochrones from different spiral arms overlap (see Fig. 4).

We only used the central region to analyse the extinction for two reasons: (1) The average extinction of the GC stellar population is maximum for this region, allowing us to detect the innermost spiral arms. (2) The number of available reference stars is largest given the covered area, which implies that more stars are available to analyse the extinction.

To select the reference stars belonging to each of the spiral arms, we generated a synthetic stellar population according to the best fit obtained for the CMD Ks versus H − Ks (see Sect. 3.3). The synthetic model indicates the regions where the probability of finding stars from a given spiral arm is maximum. In this way, we defined a selection box to identify them in the real data (left panel of Fig. 7). Figure 7 shows the chosen stars with which we analyse the extinction.

thumbnail Fig. 7.

Colour–magnitude diagrams Ks vs. H − Ks. Left panel: best model-fit solution obtained for the CMD Ks vs. H − Ks (see Table 2). Different colours indicate the stars contributing from each of the spiral arms. In particular, purple, green, orange, red, and black dots indicate stars from the first, second, third, and fourth spiral arms, and background sources respectively. Black solid parallelograms mark the selection of reference (mainly main-sequence) stars to analyse the extinction corresponding to each of the spiral arms. Right panel: real data considering completeness correction as indicated in Sect. 3.3. The red dashed line indicates the saturation limit of the Ks photometry (Ks = 11.5 mag). The red solid parallelograms depict the reference main sequence stars according to the best model fit as shown in the left panel. The isochrones corresponding to a stellar population of 1 Gyr and [M/H]= − 0.07 are over-plotted in red (see Sect. 3.3). The blue dashed lines indicate the stars added to the reference main sequence stars to create a Ks luminosity function for the second and third features (see Sect. 7).

6.1. Extinction curve

We computed the ratios between extinctions AJ/AH and AH/AKs using the selected stars belonging to the spiral arms, and then used these ratios to derive the extinction curve. To select the reference stars in the J band avoiding the overlap of the spiral arm features (see Fig. 7), we used the reference stars in the CMD Ks versus H − Ks that were also detected in J. To calculate the extinctions Ai (where i indicates the photometric band), we computed the intrinsic colours (J − H and H − Ks) of each of the reference stars interpolating from a PARSEC isochrone (Fig. 7) corresponding to the best-fit parameters obtained in Sect. 3.

Given the completeness due to crowding of the Ks-band data, we applied a completeness correction as explained in Sect. 3.3. As this approach is based on the removal of stars above a reference completeness limit of 50%, we generated 100 MC samples, randomly removing stars to correct for completeness (see Sect. 3.1.2 of Nogueras-Lara et al. (2020b) for further details). We computed the ratios between extinctions following the equation:

(1)

where i and j are the photometric bands, the subindex ‘0’ indicates intrinsic colour and extj is the mean extinction to a spiral arm as indicated in Table 2. We averaged over all the reference stars for a given spiral arm and computed the uncertainties quadratically propagating the uncertainties corresponding to the magnitudes involved in Eq. (1). In particular, we considered the ZP systematic uncertainty of 0.04 mag for each band in the calculation, and the uncertainty of the mean extinction to a given spiral arm.

We obtained the final values and uncertainties as the mean of the individual results of each of the 100 MC iterations. Table 3 summarises the results. We find that the extinction ratios AJ/AH and AH/AKs do not vary for the different spiral arms considered, within the uncertainties. We used a weighted mean to compute mean values of the ratios taking into account the larger uncertainties obtained for the first spiral arm, where the lower number of stars and the possible contamination by stars close to the Sun might influence the results. We ended up with AJ/AH = 1.89 ± 0.11 and AH/AKs = 1.86 ± 0.11. These values are in agreement with the recent results obtained for the GC (AJ/AH = 1.87 ± 0.03 and AH/AKs = 1.84 ± 0.03) by Nogueras-Lara et al. (2019b, 2020b). Therefore, our results point towards a more complex extinction curve in the NIR than the simple power-law approach that is generally accepted (e.g., Nishiyama et al. 2006). This more complex behaviour of the extinction curve might explain the different values obtained in the literature (e.g., Nishiyama et al. 2009; Stead & Hoare 2009; Fritz et al. 2011; Alonso-García et al. 2017), as discussed in Sect. 7 in Nogueras-Lara et al. (2020b).

Table 3.

Extinction curve using spiral arm reference stars.

These values are in agreement with the results obtained for the GC (AJ/AH = 1.87 ± 0.03 and AH/AKs = 1.84 ± 0.03) by Nogueras-Lara et al. (2020b). The constant extinction ratios for stars belonging to different structures and distances indicate that the extinction curve in the NIR bands JHKs towards the GC does not vary significantly with distance. Therefore, the extinction law obtained for the GC (e.g., Nogueras-Lara et al. 2019b, 2020b) is also applicable when studying the foreground population corresponding to the spiral arm structure from the Earth to the GC.

6.2. Extinction maps

Once we had determined the spiral arm structure, we built Ks-extinction maps corresponding to each of the spiral arms in order to further analyse how the extinction varies along the line of sight towards the GC. We used all the available stars within the selection box without applying the completeness correction. This is because detecting less stars due to completeness problems will not influence the values of the extinction. On the contrary, using all the stars allows us to increase the number of reference stars to improve the quality of the extinction maps.

6.2.1. Methodology

We built the extinction maps using stars within the reference boxes shown in Fig. 7. We converted the observed stellar colour H − Ks into extinction following the equation:

(2)

where H and Ks are the observed magnitudes, AH/AKs is the relation between extinctions obtained in Nogueras-Lara et al. (2020b) and is equal to 1.84 ± 0.03, and (H − Ks)0 refers to the intrinsic colour, computed for each star interpolating from a PARSEC isochrone as explained in Sect. 6.1. Equation (2) was adapted from Eq. (5) in Nogueras-Lara et al. (2018a) to avoid using the effective wavelengths, given that they vary for different stellar types and extinctions (e.g., Nogueras-Lara et al. 2020b).

To create the extinction maps, we applied the approach described in Nogueras-Lara et al. (2018a, 2020a), but using the reference stars previously specified instead of RC stars. Namely, we defined a pixel size of 20″ (larger than in the Nogueras-Lara et al. 2018a, 2020a given the significantly lower number of reference stars, and the significantly lower influence of the differential extinction for the spiral arm than for the GC stellar population), and computed the extinction for each pixel using at least five reference stars within a radius of 30″ from the centre of each pixel. We also applied an inverse distance weight method (p = 0.25) to consider the different distances of each reference star to the centre of the pixel (see Sect. 7 of Nogueras-Lara et al. 2018a for further details). We computed the statistical uncertainties considering the variation of the extinction for a given pixel using a jackknife algorithm. Moreover, we estimated the systematic uncertainty taking into account the uncertainties of AH/AKs and the ZP systematic uncertainty (0.04 mag for both H and Ks bands Nogueras-Lara et al. 2019a).

Figure 8 and Table 4 show the results. We obtained that the extinction maps for the different spiral arms are quite homogeneous and do not present any significant extinction variations. The largest variations are measured for the first spiral arm, where the lower number of stars and the possible influence of stars in the vicinity of the Sun (introducing large differential distances between the stars within this map) might contaminate our results. Moreover, the lower value of the mean extinction for the first layer makes any relative difference more visible due to the denominator of Eq. (2).

thumbnail Fig. 8.

Left column: extinction maps corresponding to each spiral arm. Right column: associated uncertainties obtained using a jackknife algorithm (see main text for further details). White pixels indicate regions where the number of reference stars is not enough to compute an extinction value. The white rectangle in the upper part of the fourth spiral arm corresponds to a region of poor quality data in H band (see Table A.1. in Nogueras-Lara et al. 2019a). Given the increase in extinction for further spiral arms, the colour scale is different for each of the extinction maps. The extinction range is always the same, 0.35 mag, in order to make small extinction variations within the same extinction layer visible.

Table 4.

Extinction values from the extinction maps.

To check the homogeneity of the derived extinction layers, Fig. 9 shows histograms of the extinction computed for the spiral arm reference stars using the corresponding extinction maps. We obtained that a Gaussian model properly fits the data for all of the spiral arms. The standard deviations of the distributions (over-plotted in Fig. 9) are similar in all the cases, but are somewhat larger for the first spiral arm as expected given the reasons explained above. We concluded that the extinction is homogeneous within each spiral arm feature and that there is no overlap between them. Therefore, they appear to be independent extinction layers. We also compared them with an extinction map obtained for a region of the NB (Nogueras-Lara et al. 2020a). In this case, the Gaussian fit is not adequate given the inhomogeneity and the larger extinction variations within the line of sight. On the other hand, we compared the mean values obtained in Table 4 with the results from the analysis of the CMD Ks versus H − Ks (Table 2) and checked that they agree within the uncertainties.

thumbnail Fig. 9.

Histograms of the AKs values per pixel in the extinction maps. The purple, green, orange, and red histograms correspond to the spiral arms. The grey histogram was obtained using the extinction map in Nogueras-Lara et al. (2020a), corresponding to the central region of the NSD. Given the smaller area of the NSD extinction map, the spiral arm histograms were properly scaled. Filled coloured contours show Gaussian fits, whose mean and standard deviation are indicated in the figure.

6.2.2. Inter-arm extinction

We computed the differential extinction per kiloparsec (AKs/d) for each of the spiral arms. We used the results from Table 4 and the distances to each of the spiral arms obtained in Sect. 3.5. We find 0.07 ± 0.01, 0.28 ± 0.07, 0.43 ± 0.12, and 0.45 ± 0.24 mag kpc−1 for each of the spiral arms, from the closest to the furthest, respectively. The uncertainties were computed quadratically propagating the uncertainties of the AKs and the distances. As the main goal is to compare the differential extinction per kiloparsec for the spiral arms, we did not consider the systematic uncertainty of AKs because it would influence all the measurements in the same direction (adding them would increase the uncertainties of AKs/d up to 0.04, 0.11, 0.14, and 0.29 mag kpc−1). Figure 10 shows the relation between the extinction and the distance to a given spiral arm. The red dashed lines join the experimental points and correspond to linear fits between consecutive spiral arms, whose slope is the differential extinction per kiloparsec computed previously.

thumbnail Fig. 10.

Extinction vs. distance. The blue points correspond to the results obtained for the distance (see Table 2) and extinction AKs (Table 4) obtained for the spiral arms features. The dashed red line shows a linear fit of the data points, two by two. The green error bar indicates the systematic uncertainty of the extinction, which does not affect the relative variation between arms because all the arms are affected in the same way.

We find an increase in the differential extinction per kiloparsec along the line of sight towards the GC. We also find that the increase in ratio (ri, j = (AKs/d)i/(AKs/d)j) between consecutive arms (i,j) is r1, 2 = 4.00 ± 0.07, r2, 3 = 1.54 ± 0.24, and r3, 4 = 1.04 ± 0.58. This might indicate an increase in material between arms that would be in agreement with the exponential thin disc of the Milky Way (e.g., Bland-Hawthorn & Gerhard 2016). Other possibilities, such as different grain composition of the interstellar dust, and/or different spiral arm width for the observed line of sight, might also be possible. A number of recent studies aim at analysing the dust distribution along the Galactic plane (e.g., Rezaei et al. 2018; Green et al. 2019; Lallement et al. 2019; Hottier et al. 2020). These are normally limited by a 3 kpc sphere imposed by the limitations of the second release of the Gaia survey. Nevertheless, we can compare their findings with our results for the closest spiral arms. Although the dust distribution seems to be very patchy and does not obviously follow continuous spiral arm footprints, it is possible to see an increase in extinction in the line of sight towards the GC (see Fig. 7 in Hottier et al. 2020) and an increase in dust density (Fig. 4 in Rezaei et al. 2018), which might be in agreement with our results. Deeper photometry and distance measurements are needed for a more sensitive analysis of the innermost regions of the Galactic disc. Our analysis and the GALACTICNUCLEUS survey might be useful for the next steps in tracing and measuring the dust distribution and extinction along the line of sight towards the GC.

7. Stellar population analysis

7.1. Luminosity function

We created de-reddened Ks-luminosity functions (KLFs) and fitted them with a linear combination of theoretical luminosity functions (Parsec models) to determine the stellar population and the SFH of the detected spiral arm features. The analysis of the KLFs is based on the relative weight of the main features appearing on the luminosity functions, such as the RC, the red giant branch bump, or the asymptotic giant branch bump (e.g., Nogueras-Lara et al. 2020a; Schödel et al. 2020). Given that the uncertainties of the luminosity functions are proportionally related to the detected number of stars (∼#stars1/2), a sufficient number of stars is needed to apply the proposed methodology. Therefore, we excluded the first spiral arm from the analysis. Moreover, some contamination from the Galactic bulge might affect the fourth spiral arm given the proximity in extinction and the use of the CMD Ks versus H − Ks, to select the stars belonging to each spiral arm. Hence, we restricted the analysis to the second and the third spiral arms.

First of all, to correct possible saturation problems and also include stars with Ks < 11.5 mag, we used the SIRIUS IRSF catalogue (e.g., Nagayama et al. 2003; Nishiyama et al. 2006) to replace the Ks photometry of stars with Ks < 11.5 mag. We accounted for possible deviations of the ZP, correcting the photometry from SIRIUS with respect to non-saturated common bright stars between both catalogues. We also included the photometry of bright stars that were not detected in Ks band in the GALACTICNUCLEUS catalogue because of saturation. We selected stars within the boxes used to create the extinction maps and also included the stars belonging to the ascending giant branch following the isochrones, as indicated by the blue dashed lines in Fig. 7. We established a lower limit of Ks = 9 mag to account for the saturation of the SIRIUS survey (Matsunaga et al. 2009). In this way, we are able to include more red giant stars whose main features in the KLF, in particular the RC feature (e.g., Girardi 2016), are very useful for disentangling the different stellar populations (e.g., Nogueras-Lara et al. 2018b, 2020a).

We corrected the reddening, applying the extinction maps derived previously for the corresponding spiral arms (Sect. 6.2). The KLFs were computed selecting the bin-width using the Python function numpy.histogram (Harris et al. 2020). We used the completeness solution due to crowding obtained in Nogueras-Lara et al. (2020b) for the Ks band. We limited the correction to magnitudes where the completeness is > 70%.

To fit the KLFs, we used a similar approach as in Nogueras-Lara et al. (2020a), Schödel et al. (2020). We generated 1000 MC samples from the original KLFs, assuming Gaussian uncertainties and fit them with a linear combination of 14 theoretical PARSEC luminosity functions3. We used the following models (assuming metallicities in agreement with Feltzing et al. 2001 for the different ages, as in previous sections, and a Kroupa initial mass function) to sample the age space: 12, 9, 7, 5, 3, 2, 1, 0.7, 0.5, 0.3, 0.1, 0.06, 0.03, and 0.01 Gyr. We considered a free parameter in the fit for a Gaussian smoothing to take into account different distances to the stars and/or possible photometric uncertainties. We also fitted the distance modulus allowing it to vary in a 3 σ range from the values measured in Table 2 for the four spiral arms J − H case.

The best fit for each of the MC realisations was obtained via χ2 minimisation. To avoid the degeneracies between similar ages (that is more important for old stellar populations), we created larger bins for the ages adding the contribution from several models. Figure 11 shows the selected bins and the results obtained for the second and third spiral arms when averaging over the 1000 MC realisations, removing 3 σ outliers. The uncertainties were computed as the standard deviation of the underlying distributions. Figure 12 shows the best fit obtained for the original KLFs, in agreement with the results of the MC simulations.

thumbnail Fig. 11.

SFHs derived for the second and the third spiral arms, as derived from the MC simulations. The error bars show the 1 σ uncertainties. The absence of any bar and the error bar indicate that there is no contribution for a given age bin.

thumbnail Fig. 12.

Best-fit solutions obtained for the original (no MC samples) of the KLFs for the second and third spiral arms (left and right panels, respectively). The legend indicates the colour coding corresponding to the selected age bins. The is indicated for each of the cases.

7.2. Systematic uncertainty

To check our results, we explored potential sources of systematic uncertainty. In particular, we varied the bin width of the original KLFs considering a variation of ±25%. We repeated the MC simulations without obtaining any significant variation within the uncertainties. We also analysed the influence of the models in our results. For this, we used the updated BaSTI models4 (Pietrinferni et al. 2004, 2013; Bedin et al. 2005; Cordier et al. 2007; Hidalgo et al. 2018) with solar metallicity, a Salpeter initial mass function, and slightly different ages (13, 11, 9, 7, 5, 3, 1, 0.8, 0.5, 0.25, 0.15, 0.1, 0.08, and 0.03 Gyr). We obtained similar results within the uncertainties for all the stellar ages, except for the second spiral arm, where the bin for ages between 0.6 and 1.5 Gyr appears to have a somewhat larger (∼5 ± 3%) contribution. This is not surprising given the lower Ks range used to derive the SFH. In particular, the KLF does not cover the asymptotic giant branch bump and the RC feature is not fully included, which might introduce some bias when deriving the SFH. On the other hand, these features are present in the KLF belonging to the third spiral arm (∼9.5 and ∼11.5 mag, respectively).

7.3. Results

We found relatively similar stellar populations for both spiral arms. We identified two main episodes of star formation: (1) the majority of the stellar mass (∼60 − 70%) was formed more than 6 Gyr ago. (2) A second important burst appears for ages between 1.5 and 4 Gyr, where ∼20 − 30% of the stellar mass originated. We also measured some significant star formation at more recent ages for the case of the third spiral arm, accounting for ∼10% of the stellar mass.

Our results agree well with previous stellar population studies on the Milky Way disc. In this sense, Rocha-Pinto et al. (2000) detected three enhanced episodes of star formation of the Milky Way disc around 0–1 Gyr, 2–5 Gyr, and 7–9 Gyr ago that agree well with our findings, in particular with the SFH for the third spiral arm, where we also detected clear signatures of recent star formation. On the other hand, the more recent work by Ruiz-Lara et al. (2020) analysed the stellar population within a radius of ∼2 kpc from the Sun and found three significant star formation events around 5.7, 1.9, and 1.0 Gyr ago. The first two agree with the SFH derived in this work, whereas the third one (which is lower in their analysis) appears to be shifted towards younger ages in our results, and/or not to be very important for the case of the second spiral arm. Overall, our results point towards a similar SFH in the solar neighbourhood to that in the spiral arms towards the GC that we analyse here. This also indicates that the models used to fit the spiral arms in Sect. 3 are adequate.

8. Conclusions

In the present study, we computed the distance to the spiral arms along the line of sight towards the GC and analysed the extinction in the NIR using the data from the GALACTICNUCLEUS survey (Nogueras-Lara et al. 2018a, 2019a). We investigated the spiral arm structure by fitting a synthetic model to the CMD H versus J − H. We applied a χ2 minimisation technique and computed the distances and extinction to the spiral-arm features, obtaining: d1 = 1.6 ± 0.2 kpc, d2 = 2.6 ± 0.2 kpc, d3 = 3.9 ± 0.3 kpc, and d4 = 4.5 ± 0.2 kpc, and AH1 = 0.35 ± 0.08 mag, AH2 = 0.77 ± 0.08 mag, AH3 = 1.68 ± 0.08 mag, and AH4 = 2.30 ± 0.08 mag. The sub-indices refer to the corresponding spiral arms.

We checked that our results are robust and addressed potential sources of systematic uncertainty, repeating the calculation using different ages for the synthetic models, bin sizes, photometric bands, SFHs, and simulated models. We also considered the possible effect of the ZP systematic error. Moreover, we analysed the CCD distribution of the population in front of the GC corresponding to the central region of the GALACTICNUCLEUS survey and find that it is compatible with the four-spiral-arms structure. Finally, we cross-correlated the GALACTICNUCLEUS data with the Gaia DR2 catalogue (Gaia Collaboration 2018) and estimated distances corresponding to the two nearest spiral arms. These are consistent with our results.

We calculated the ratios between the extinction values AH and AKs, selecting stars belonging to each of the spiral-arm features and using the mean AKs extinctions derived with the synthetic fitting of the CMD. We conclude that the extinction curve does not significantly change with distance and the mean absolute extinction. Combining the values for the different spiral arms, we obtain AJ/AH = 1.91 ± 0.11 and AH/AKs = 1.86 ± 0.11. These values are compatible, within the uncertainties, with the extinction ratios of AJ/AH = 1.87 ± 0.03 AH/AKs = 1.84 ± 0.03 derived for the GC (Nogueras-Lara et al. 2020b), pointing towards a constant extinction law in the studied wavelength regime.

We also created AKs-extinction maps using the features in the CMD Ks versus H − Ks for the central region of the GALACTICNUCLEUS survey. We find that the extinction maps are homogeneous within the uncertainties and do not present any significant variations in the studied line of sight, in contrast to the extinction maps corresponding to the innermost regions of the Galaxy (e.g., Nogueras-Lara et al. 2018a, 2020a). Moreover, we find that the ratio of extinction to distance increases towards the GC, indicating a higher dust density for the inner spiral arms, and/or different dust composition or spiral arm width. We also find that ∼50% of the cumulative extinction along the line of sight towards the GC occurs within a distance of ∼5 kpc from Earth.

Finally, we studied the stellar populations and the SFH of the second and third spiral arms, analysing their KLFs and fitting them with a linear combination of theoretical models of different ages (Parsec KLFs). We ended up with similar stellar populations for both spiral arms, detecting two main episodes of star formation around > 6 Gyr (∼60 − 70% of the stellar mass) and 1.5 − 4 Gyr (∼20 − 30% of the stellar mass), and also more recent star formation (∼10%) for ages < 1 Gyr for the third spiral arm. Our results agree well with those of previous studies of the Galactic disc stellar population (e.g., Rocha-Pinto et al. 2000) and also suggest a similar SFH with respect to the stars within a 2 kpc radius around the Sun (Ruiz-Lara et al. 2020).


Acknowledgments

This work is based on observations made with ESO Telescopes at the La Silla Paranal Observatory under programme ID 195.B-0283. We thank the staff of ESO for their great efforts and helpfulness. F.N.-L. and N.N. gratefully acknowledge support by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) – Project-ID 138713538 – SFB 881 (“The Milky Way System”, subproject B8). R.S. acknowledges financial support from the State Agency for Research of the Spanish MCIU through the “Center of Excellence Severo Ochoa” award for the Instituto de Astrofísica de Andalucía (SEV-2017-0709). R. S. acknowledges financial support from national project PGC2018-095049-B-C21 (MCIU/AEI/FEDER, UE).

References

  1. Akaike, H. 1974, IEEE Trans. Autom. Control, 19, 716 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
  2. Alonso-García, J., Minniti, D., Catelan, M., et al. 2017, ApJ, 849, L13 [Google Scholar]
  3. Bailer-Jones, C. A. L., Rybizki, J., Fouesneau, M., Mantelet, G., & Andrae, R. 2018, AJ, 156, 58 [NASA ADS] [CrossRef] [Google Scholar]
  4. Bedin, L. R., Cassisi, S., Castelli, F., et al. 2005, MNRAS, 357, 1038 [NASA ADS] [CrossRef] [Google Scholar]
  5. Bergemann, M., Ruchti, G. R., Serenelli, A., et al. 2014, A&A, 565, A89 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  6. Bland-Hawthorn, J., & Gerhard, O. 2016, ARA&A, 54, 529 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  7. Bressan, A., Marigo, P., Girardi, L., et al. 2012, MNRAS, 427, 127 [Google Scholar]
  8. Cabrera-Lavers, A., González-Fernández, C., Garzón, F., Hammersley, P. L., & López-Corredoira, M. 2008, A&A, 491, 781 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  9. Chen, Y., Girardi, L., Bressan, A., et al. 2014, MNRAS, 444, 2525 [NASA ADS] [CrossRef] [Google Scholar]
  10. Chen, Y., Bressan, A., Girardi, L., et al. 2015, MNRAS, 452, 1068 [NASA ADS] [CrossRef] [Google Scholar]
  11. Chen, B. Q., Liu, X. W., Yuan, H. B., et al. 2017, MNRAS, 464, 2545 [NASA ADS] [CrossRef] [Google Scholar]
  12. Churchwell, E., Babler, B. L., Meade, M. R., et al. 2009, PASP, 121, 213 [NASA ADS] [CrossRef] [Google Scholar]
  13. Cordier, D., Pietrinferni, A., Cassisi, S., & Salaris, M. 2007, AJ, 133, 468 [NASA ADS] [CrossRef] [Google Scholar]
  14. Dolphin, A. E. 2002, MNRAS, 332, 91 [NASA ADS] [CrossRef] [Google Scholar]
  15. Drimmel, R., & Spergel, D. N. 2001, ApJ, 556, 181 [NASA ADS] [CrossRef] [Google Scholar]
  16. Englmaier, P., & Gerhard, O. 1999, MNRAS, 304, 512 [NASA ADS] [CrossRef] [Google Scholar]
  17. Feltzing, S., Holmberg, J., & Hurley, J. R. 2001, A&A, 377, 911 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  18. Francis, C., & Anderson, E. 2012, MNRAS, 422, 1283 [NASA ADS] [CrossRef] [Google Scholar]
  19. Fritz, T. K., Gillessen, S., Dodds-Eden, K., et al. 2011, ApJ, 737, 73 [NASA ADS] [CrossRef] [Google Scholar]
  20. Gaia Collaboration (Brown, A. G. A., et al.) 2018, A&A, 616, A1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  21. Girardi, L. 2016, ARA&A, 54, 95 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  22. Gonzalez, O. A., Rejkuba, M., Minniti, D., et al. 2011, A&A, 534, L14 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  23. Gonzalez, O. A., Minniti, D., Valenti, E., et al. 2018, MNRAS, 481, L130 [NASA ADS] [CrossRef] [Google Scholar]
  24. Green, G. M., Schlafly, E., Zucker, C., Speagle, J. S., & Finkbeiner, D. 2019, ApJ, 887, 93 [NASA ADS] [CrossRef] [Google Scholar]
  25. Harris, C. R., Millman, K. J., van der Walt, S. J., et al. 2020, Nature, 585, 357 [CrossRef] [PubMed] [Google Scholar]
  26. Hidalgo, S. L., Pietrinferni, A., Cassisi, S., et al. 2018, ApJ, 856, 125 [Google Scholar]
  27. Hottier, C., Babusiaux, C., & Arenou, F. 2020, A&A, 641, A79 [CrossRef] [EDP Sciences] [Google Scholar]
  28. Hou, L. G., & Han, J. L. 2014, A&A, 569, A125 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  29. Khoperskov, S., Gerhard, O., Di Matteo, P., et al. 2020, A&A, 634, L8 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  30. Kissler-Patig, M., Pirard, J.-F., Casali, M., et al. 2008, A&A, 491, 941 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  31. Lallement, R., Babusiaux, C., Vergely, J. L., et al. 2019, A&A, 625, A135 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  32. Launhardt, R., Zylka, R., & Mezger, P. G. 2002, A&A, 384, 112 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  33. Lucas, P. W., Hoare, M. G., Longmore, A., et al. 2008, MNRAS, 391, 136 [NASA ADS] [CrossRef] [Google Scholar]
  34. Luri, X., Brown, A. G. A., Sarro, L. M., et al. 2018, A&A, 616, A9 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  35. Marigo, P., Girardi, L., Bressan, A., et al. 2017, ApJ, 835, 77 [NASA ADS] [CrossRef] [Google Scholar]
  36. Matsunaga, N., Kawadu, T., Nishiyama, S., et al. 2009, MNRAS, 399, 1709 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
  37. Mighell, K. J. 1999, ApJ, 518, 380 [Google Scholar]
  38. Minniti, D., Lucas, P. W., Emerson, J. P., et al. 2010, New Astron., 15, 433 [Google Scholar]
  39. Minniti, D., Saito, R. K., Gonzalez, O. A., et al. 2014, A&A, 571, A91 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  40. Momany, Y., Zaggia, S., Gilmore, G., et al. 2006, A&A, 451, 515 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  41. Nagayama, T., Nagashima, C., Nakajima, Y., et al. 2003, in Instrument Design and Performance for Optical/Infrared Ground-based Telescopes, eds. M. Iye, A. F. M. Moorwood, et al., Proc. SPIE, 4841, 459 [Google Scholar]
  42. Nataf, D. M., Gonzalez, O. A., Casagrande, L., et al. 2016, MNRAS, 456, 2692 [NASA ADS] [CrossRef] [Google Scholar]
  43. Nishiyama, S., Nagata, T., Kusakabe, N., et al. 2006, ApJ, 638, 839 [Google Scholar]
  44. Nishiyama, S., Nagata, T., Tamura, M., et al. 2008, ApJ, 680, 1174 [NASA ADS] [CrossRef] [Google Scholar]
  45. Nishiyama, S., Tamura, M., Hatano, H., et al. 2009, ApJ, 696, 1407 [NASA ADS] [CrossRef] [Google Scholar]
  46. Nogueras-Lara, F., Gallego-Calvente, A. T., Dong, H., et al. 2018a, A&A, 610, A83 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  47. Nogueras-Lara, F., Schödel, R., Dong, H., et al. 2018b, A&A, 620, A83 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  48. Nogueras-Lara, F., Schödel, R., Gallego-Calvente, A. T., et al. 2019a, A&A, 631, A20 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  49. Nogueras-Lara, F., Schödel, R., Najarro, F., et al. 2019b, A&A, 630, L3 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  50. Nogueras-Lara, F., Schödel, R., Gallego-Calvente, A. T., et al. 2020a, Nat. Astron., 4, 377 [CrossRef] [Google Scholar]
  51. Nogueras-Lara, F., Schödel, R., Neumayer, N., et al. 2020b, A&A, 641, A141 [EDP Sciences] [Google Scholar]
  52. Pastorelli, G., Marigo, P., Girardi, L., et al. 2019, MNRAS, 485, 5666 [Google Scholar]
  53. Pedregosa, F., Varoquaux, G., Gramfort, A., et al. 2011, J. Mach. Learn. Res., 12, 2825 [Google Scholar]
  54. Pfuhl, O., Fritz, T. K., Zilka, M., et al. 2011, ApJ, 741, 108 [NASA ADS] [CrossRef] [Google Scholar]
  55. Pietrinferni, A., Cassisi, S., Salaris, M., & Castelli, F. 2004, ApJ, 612, 168 [NASA ADS] [CrossRef] [Google Scholar]
  56. Pietrinferni, A., Cassisi, S., Salaris, M., & Castelli, F. 2006, ApJ, 642, 797 [NASA ADS] [CrossRef] [Google Scholar]
  57. Pietrinferni, A., Cassisi, S., Salaris, M., & Hidalgo, S. 2013, A&A, 558, A46 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  58. Reid, M. J., Dame, T. M., Menten, K. M., & Brunthaler, A. 2016, ApJ, 823, 77 [NASA ADS] [CrossRef] [Google Scholar]
  59. Reid, M. J., Menten, K. M., Brunthaler, A., et al. 2019, ApJ, 885, 131 [Google Scholar]
  60. Rezaei, K. h. S., Bailer-Jones, C. A. L., Schlafly, E. F., & Fouesneau, M. 2018, A&A, 616, A44 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  61. Robin, A. C., Marshall, D. J., Schultheis, M., & Reylé, C. 2012, A&A, 538, A106 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  62. Rocha-Pinto, H. J., Scalo, J., Maciel, W. J., & Flynn, C. 2000, A&A, 358, 869 [NASA ADS] [Google Scholar]
  63. Ruiz-Lara, T., Gallart, C., Bernard, E. J., & Cassisi, S. 2020, Nat. Astron., 142, 166 [Google Scholar]
  64. Saito, R. K., Hempel, M., Minniti, D., et al. 2012, A&A, 537, A107 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  65. Saito, R. K., Minniti, D., Benjamin, R. A., et al. 2020, MNRAS, 494, L32 [CrossRef] [Google Scholar]
  66. Schödel, R., Najarro, F., Muzic, K., & Eckart, A. 2010, A&A, 511, A18 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  67. Schödel, R., Yelda, S., Ghez, A., et al. 2013, MNRAS, 429, 1367 [NASA ADS] [CrossRef] [Google Scholar]
  68. Schödel, R., Feldmeier, A., Kunneriath, D., et al. 2014, A&A, 566, A47 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  69. Schödel, R., Nogueras-Lara, F., Gallego-Cano, E., et al. 2020, A&A, 641, A102 [CrossRef] [EDP Sciences] [Google Scholar]
  70. Schönrich, R., McMillan, P., & Eyer, L. 2019, MNRAS, 487, 3568 [Google Scholar]
  71. Schultheis, M., Zasowski, G., Allende Prieto, C., et al. 2014, AJ, 148, 24 [NASA ADS] [CrossRef] [Google Scholar]
  72. Schwarz, G. 1978, Ann. Stat., 6, 461 [Google Scholar]
  73. Skrutskie, M. F., Cutri, R. M., Stiening, R., et al. 2006, AJ, 131, 1163 [Google Scholar]
  74. Stead, J. J., & Hoare, M. G. 2009, MNRAS, 400, 731 [NASA ADS] [CrossRef] [Google Scholar]
  75. Tang, J., Bressan, A., Rosenfield, P., et al. 2014, MNRAS, 445, 4287 [NASA ADS] [CrossRef] [Google Scholar]
  76. Vallée, J. P. 2014, AJ, 148, 5 [NASA ADS] [CrossRef] [Google Scholar]
  77. Vallée, J. P. 2016, AJ, 151, 55 [NASA ADS] [CrossRef] [Google Scholar]
  78. Wu, Y. W., Reid, M. J., Sakai, N., et al. 2019, ApJ, 874, 94 [NASA ADS] [CrossRef] [Google Scholar]

Appendix A: Systematic errors of the synthetic CMD model fitting

Table A.1.

Influence of bin width on the four-spiral-arm model using the CMD H vs. J − H.

Table A.2.

Influence of the ZP systematic errors on the four-spiral-arm model using the CMD H vs. J − H.

Table A.3.

Influence of the synthetic stellar populations on the four spiral arm model using the CMD H vs. J − H.

Table A.4.

Simulations to test the reliability of the model fitting.

All Tables

Table 1.

Parameters to create the models.

Table 2.

Model fitting results.

Table 3.

Extinction curve using spiral arm reference stars.

Table 4.

Extinction values from the extinction maps.

Table A.1.

Influence of bin width on the four-spiral-arm model using the CMD H vs. J − H.

Table A.2.

Influence of the ZP systematic errors on the four-spiral-arm model using the CMD H vs. J − H.

Table A.3.

Influence of the synthetic stellar populations on the four spiral arm model using the CMD H vs. J − H.

Table A.4.

Simulations to test the reliability of the model fitting.

All Figures

thumbnail Fig. 1.

False colour image from the GALACTICNUCLEUS survey (using JHKs bands) showing the regions used in this study. The red arrow marks the position of the supermassive black hole, Sagittarius A* (Sgr A*). The names of the regions are also indicated in the figure: central, transition east (TE), transition west (TW), inner bulge south (IBS), and inner bulge north (IBN). The distances between the regions correspond to real distances, as indicated by the scaling arrows.

In the text
thumbnail Fig. 2.

Colour-magnitude diagram Ks vs. J − Ks from the central region of the GALACTICNUCLEUS survey. The red dashed line indicates the colour cut between the stellar population belonging to the GC: Galactic bulge (GB) + Nuclear bulge (NB), and the Galactic disc (GD). The contours outline the population belonging to the GD.

In the text
thumbnail Fig. 3.

Model fitting performance for a given MC sample (left panel). Middle panel: best fit for the MC sample (). Right panel: residuals after subtracting the model from the data. The colour bar indicates the number of stars per bin using a linear scale.

In the text
thumbnail Fig. 4.

CMDs H vs. J − H and Ks vs. H − Ks for the central, transition, and inner bulge regions of the GALACTICNUCLEUS survey, from the upper to the lower panels, respectively. The colour code refers to stellar densities, using a power-stretch scale (different in each case to stress the spiral arm features). Cyan lines correspond to a PARSEC isochrone of 1 Gyr and [M/H = −0.07] according to the extinction and distance parameters obtained for the spiral arm. The blue and green error bars indicate the 1 σ uncertainties of the position of the isochrones and the photometric ZP, respectively. The zoomed-in region corresponds to the red clump features (e.g., Girardi 2016) from the inner Galactic bulge (‘1’ and ‘2’ associated to different extinctions from the IBS and IBS regions) and the overlap with the fourth spiral arm.

In the text
thumbnail Fig. 5.

Histograms of the distances to the common stars between Gaia and GALACTICNUCLEUS detected in the central, transition, and inner bulge regions (upper, middle, and lower panel, respectively). The solid green lines indicate the two-Gaussian fits, whereas the red and cyan dashed lines show each of the two Gaussian components. The mean values and their errors for each component are indicated for each panel in units of kpc.

In the text
thumbnail Fig. 6.

Left panel: colour–colour diagram J − H vs. H − Ks corresponding to the central region of the GALACTICNUCLEUS survey. The black arrows indicate the over-density features that correspond to the spiral arms detected in Sect. 3. The blue dashed lines mark the stellar population belonging to the GC, where the density is significantly different from the rest. Right panel: zoom into the foreground population from the spiral arms and histograms corresponding to the underlying distributions. The red solid line indicates the result of a four-Gaussian fit (green dashed lines for each individual Gaussian). The black arrows show the mean location of the spiral arms. The red arrow indicates the overlap between the first and the second spiral arms in one single feature in the histogram J − H.

In the text
thumbnail Fig. 7.

Colour–magnitude diagrams Ks vs. H − Ks. Left panel: best model-fit solution obtained for the CMD Ks vs. H − Ks (see Table 2). Different colours indicate the stars contributing from each of the spiral arms. In particular, purple, green, orange, red, and black dots indicate stars from the first, second, third, and fourth spiral arms, and background sources respectively. Black solid parallelograms mark the selection of reference (mainly main-sequence) stars to analyse the extinction corresponding to each of the spiral arms. Right panel: real data considering completeness correction as indicated in Sect. 3.3. The red dashed line indicates the saturation limit of the Ks photometry (Ks = 11.5 mag). The red solid parallelograms depict the reference main sequence stars according to the best model fit as shown in the left panel. The isochrones corresponding to a stellar population of 1 Gyr and [M/H]= − 0.07 are over-plotted in red (see Sect. 3.3). The blue dashed lines indicate the stars added to the reference main sequence stars to create a Ks luminosity function for the second and third features (see Sect. 7).

In the text
thumbnail Fig. 8.

Left column: extinction maps corresponding to each spiral arm. Right column: associated uncertainties obtained using a jackknife algorithm (see main text for further details). White pixels indicate regions where the number of reference stars is not enough to compute an extinction value. The white rectangle in the upper part of the fourth spiral arm corresponds to a region of poor quality data in H band (see Table A.1. in Nogueras-Lara et al. 2019a). Given the increase in extinction for further spiral arms, the colour scale is different for each of the extinction maps. The extinction range is always the same, 0.35 mag, in order to make small extinction variations within the same extinction layer visible.

In the text
thumbnail Fig. 9.

Histograms of the AKs values per pixel in the extinction maps. The purple, green, orange, and red histograms correspond to the spiral arms. The grey histogram was obtained using the extinction map in Nogueras-Lara et al. (2020a), corresponding to the central region of the NSD. Given the smaller area of the NSD extinction map, the spiral arm histograms were properly scaled. Filled coloured contours show Gaussian fits, whose mean and standard deviation are indicated in the figure.

In the text
thumbnail Fig. 10.

Extinction vs. distance. The blue points correspond to the results obtained for the distance (see Table 2) and extinction AKs (Table 4) obtained for the spiral arms features. The dashed red line shows a linear fit of the data points, two by two. The green error bar indicates the systematic uncertainty of the extinction, which does not affect the relative variation between arms because all the arms are affected in the same way.

In the text
thumbnail Fig. 11.

SFHs derived for the second and the third spiral arms, as derived from the MC simulations. The error bars show the 1 σ uncertainties. The absence of any bar and the error bar indicate that there is no contribution for a given age bin.

In the text
thumbnail Fig. 12.

Best-fit solutions obtained for the original (no MC samples) of the KLFs for the second and third spiral arms (left and right panels, respectively). The legend indicates the colour coding corresponding to the selected age bins. The is indicated for each of the cases.

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.