Open Access
Issue
A&A
Volume 615, July 2018
Article Number A100
Number of page(s) 7
Section Galactic structure, stellar clusters and populations
DOI https://doi.org/10.1051/0004-6361/201732568
Published online 20 July 2018

© ESO 2018

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

1 Introduction

In recent years, major progress has been made in our understanding of the Galactic bulge. Large photometric and/or spectroscopic data sets have allowed us to learn about the physical properties of the Galactic bulge components such as the three-dimensional structure, metallicity distribution, kinematics, and age, but their number and the corresponding mass budget is still under debate.

From spectroscopic studies of the bulge metallicity distribution function at different galactic longitudes and latitudes, Babusiaux et al. (2010), Hill et al. (2011), Uttenthaler et al. (2012), Rojas-Arriagada et al. (2014, 2017), Gonzalez et al. (2015), Schultheis et al. (2017), and Zoccali et al. (2017) found two main stellar components: a metal-rich component centred on [Fe/H] ~ 0.3 – 0.4 dex and a metal-poor component centred on [Fe/H] ~ –0.3 – –0.4 dex. They display distinct three-dimensional structure and kinematic behaviour. Metal-rich stars ([Fe/H] ≥ 0 dex) follow the X-shaped structure of a boxy peanut (McWilliam & Zoccali 2010; Nataf et al. 2010), show cylindrical rotation (Howard et al. 2008, 2009; Kunder et al. 2012; Ness et al. 2013b; Zoccali et al. 2014), and streaming motions (Rangwala et al. 2009; De Propris et al. 2011; Uttenthaler et al. 2012; Poleski et al. 2013; Babusiaux et al. 2014; Rojas-Arriagada et al. 2014), which argue for a secular formation scenario of the population. Although there is consensus that the metal-rich stars support the B/P structure and are formed via internal evolution of the thin disc (e.g. Combes & Sanders 1981; Raha et al. 1991; Athanassoula 2005; Debattista et al. 2006; Martinez-Valpuesta et al. 2006; Di Matteo et al. 2014), the origin of their metal-poor counterparts is controversial. Metal-poor stars ([Fe/H] < 0 dex) do not display the X-shaped structure. They are kinematically different from metal-rich stars. They generally rotate more slowly and are dynamically hotter. They were associated with a classical bulge population, formed via mergers or by a rapid dissipative collapse at early phases of Galaxy formation. However, the ARGOS spectroscopic survey (Freeman et al. 2013) data shed light on another possible origin of low-metallicity stars. Using the full ARGOS sample, Ness et al. (2013a) identified five different components, A, B, C, D, and E, whose metallicities peak at [Fe/H] 0.1 dex, –0.3 dex, –0.7 dex, –1.2 dex, and –1.7 dex, respectively, the fraction of stars in each component changing with Galactic latitude. The most significant components A ([Fe/H] ≥ 0), B (–0.5 ≤ [Fe/H] < 0), and C (–1.0 ≤ [Fe/H] < –0.5) are consistent with a common disc origin where component C was formed out of the thick disc (Ness et al. 2013b; Di Matteo et al. 2015; Portail et al. 2017). Di Matteo et al. (2014, 2015), and Di Matteo (2016) suggest that the B component is associated with the younger thick disc (Haywood et al. 2013). The more metal-poor components D and E represent less than about 6% of stars in the ARGOS survey and may be associated with the metal-weak thick disc or with either a possibly classical bulge or inner halo, respectively (Portail et al. 2017; Gerhard 2017). Recently, Bensby et al. (2017) analysed high-resolution spectra of microlensed dwarf and subgiant stars. A very wide metallicity distribution is observed, ranging from [Fe/H] ~–2 to 0.5 dex, with more than two significant peaks. Furthermore, the age distribution is also very wide. For [Fe/H] > 0 dex, the stars span all ages, from 1 Gyr to 12–13 Gyr. Below [Fe/H] ≤ –0.5 dex, most stars are 10 Gyr or older. Moreover, alpha-element abundance trends with metallicity of bulge stars with sub-solar metallicities, even though there are differences, show similarities with the local thick disc. The authors concluded that the observed age and abundance properties suggest a secular origin for the Galactic bulge, but they cannot rule out a small contributionof classical bulge or halo stars. On the other hand, studies based on red giants in the bulge also found abundance trends with metallicity to be similar to the nearby thick disc for sub-solar metallicities (e.g. Meléndez et al. 2008; Alves-Brito et al. 2010; Ryde et al. 2010, 2016; Gonzalez et al. 2011, 2015; Johnson et al. 2014; Jönsson et al. 2017; García Pérez et al. 2018).

The metallicity distribution observed from dwarfs and giants in the bulge displays less than about 5–6% of metal-poor stars with [Fe/H] < –1.0 dex; most of them have [Fe/H] values ≥ –1.0 dex. In what follows, stars with –1.0 ≤ [Fe/H] < 0 are called moderately metal-poor stars, and those with [Fe/H] ≥ 0 are called metal-rich stars. Stars with [Fe/H] < –1.0 dex are well traced by RRLyrae stars, whose origin is also controversial. Kunder et al. (2016) suggested that the observed spatial distribution and kinematics are consistent with a classical bulge origin, although they cannot rule out the possibility that they are the metal-poor tail of a more metal-rich halo-bulge population. Instead, Pérez-Villegas et al. (2017) showed that RRLyrae stars in the bulge might be the inner extension of the Galactic stellar halo.

We here investigate whether the origin of moderately metal-poor stars (–1.0 ≤ [Fe/H] < 0) is mainly consistent with a classical bulge origin or with a formation from a rather massive thick disc present in the inner Galaxy. Using N-body simulations, we have calculated the difference in the mean Galactocentric radial velocity (Δ VGC) between metal-rich and metal-poor stars in two models; one containing a small classical bulge, and the other consisting of only disc components. The kinematic signature given by each of the two models was obtained at different lines of sight. The predicted kinematic trends were compared with data from the ARGOS survey (Freeman et al. 2013) and from the APOGEE DR13 survey (Albareti et al. 2017; Majewski et al. 2017). The results show that the observed data are consistent with a thick-disc origin for moderately metal-poor stars, leaving little room for a classical bulge origin.

The paper is organised as follows: in Sect. 2, the models are briefly described. Section 3 shows the kinematic results of the simulations at different latitudes and longitudes. The comparison to observations is shown in Sect. 4. In Sects. 5 and 6, we discuss and summarise our results, respectively.

2 Simulations

Two different high-resolution simulations are considered in this paper. The first, called the “disc+bulge” here, consists of an isolated stellar disc and a classical bulge, with a bulge-to-disc ratio equal to B/D = 0.1, and contains no gas. This simulation has been extensively described and analysed in Di Matteo et al. (2014, 2015), and Gómez et al. (2016). The dark halo and the bulge are modelled as Plummer spheres (Binney & Tremaine 1987). The dark halo has a mass MH = 1.02 × 1011M and a characteristic radius rH = 10 kpc. The bulge has a mass MB = 9 × 109M and characteristic radius rB = 1.3 kpc. The stellar disc follows a Miyamoto-Nagai density profile (Binney & Tremaine 1987), with mass M* = 9 × 1010M and vertical and radial scale lengths given by h* = 0.5 kpc and a* = 4 kpc, respectively. The initial disc size is 13 kpc, and the Toomre parameter is set equal to Q = 1.8. The galaxy is represented by 30 720 000 particles redistributed among dark matter (10 240 000) and stars (20 480 000). This simulation reproduces the observed trends of the global kinematics and the global chemical characteristics of the Galactic bulge. However, it does not reproduce the chemo-kinematic relations satisfied by the individual bulge components. Di Matteo et al. (2015), assuming an initial radial metallicity profile in the disc similar to Martinez-Valpuesta & Gerhard (2013), showed that this model fails in reproducing the observed properties of stellar components in ARGOS data (see their Fig. 4).

That the Galactic thick disc seems to be as massive in the inner Galaxy as the Galactic thin disc (Fuhrmann et al. 2012; Haywood et al. 2013; Snaith et al. 2014) led Di Matteo et al. (2014) to suggest that the Milky Way bulge is the result of mapping the Galactic thin+thick disc into the boxy/peanut-shaped structure. Thus, we have considered a second simulation, called “thin+thick” here, which consists of three isolated stellar discs corresponding to a kinematically cold thin disc, to an intermediate disc with intermediate kinematics, and to a kinematically hot disc (hereafter called “thin”, “intermediate”, and “thick” discs). They have different scale heights and lengths, with masses and sizes in agreement with recent estimates for the Milky Way (Bensby et al. 2011, 2013; Bovy et al. 2012a,b, 2016; Haywood et al. 2013; Snaith et al. 2014; Hayden et al. 2015). The discs are modelled with Miyamoto-Nagai density distributions (Binney & Tremaine 1987) with masses M* = 2.55 × 1010M, 1.53  × 1010M, and 1.02 × 1010M; scale heights h* = 0.3 kpc, 0.6 kpc, and 0.9 kpc; and scale lengths a* = 4.7 kpc, 2.3 kpc, and 2.3 kpc for the thin-disc, the intermediate-disc, and the thick-disc components, respectively. The dark halo is modelled as a Plummer sphere (Binney & Tremaine 1987) with characteristic mass and radius MH = 1.61 × 1011M and rH =10 kpc. The modelled disc galaxy consists of 25 000 000 particles, redistributed among stars (20 000 000) and dark matter (5 000 000). Of the star particles, 10 000 000 belong to the thin disc, 6 000 000 to the intermediate disc, and 4 000 000 to the thick disc. The stellar disc is structured with scale heights and velocity dispersions at the solar vicinity similar to those observed for the thin disc, the young thick disc, and the old thick disc, respectively (Haywood et al. 2013). The intermediate and thick discs together represent 50% of the stellar mass, the remaining 50% being in the thin disc. All stellar components contribute to the stellar bar and to the boxy/peanut-shaped structure, but their fraction varies with the height above the plane. Thin-disc stars are more concentrated towards the galactic plane, the fraction of thick-disc stars increases with height above the plane, and that of intermediate-disc stars stays nearly constant (Di Matteo 2016), with proportions similar to those for populations A, B, and C by Ness et al. (2013a). At a given height above the plane, the boxy/peanut-shaped structure is more pronounced in the kinematically cold populations than in the hottest population (Di Matteo 2016; Fragkoudi et al. 2017a). Results on the kinematic characteristics predicted by this disc model are given in Di Matteo et al. (in prep.).

In an N-body simulation the length of the bar is difficult to control, and it is therefore necessary to rescale the model to match the length of the bar of the Milky Way. Both models were rescaled to match the Milky Way bar size and bulge velocities; see Di Matteo et al. (2015) for the disc+bulge simulation and Di Matteo et al. (in prep.) for the thin+thick simulation, the Sun being placed at 8 kpc from the Galactic centre. For consistency, we adopted a bar orientation relative to the Sun-Galactic centre of 20° as in Di Matteo et al. (2015) and Di Matteo (2016). As a consequence of the observed three-dimensional structure of the Milky Way box/peanut-bulge (e.g. Dwek et al. 1995; Cao et al. 2013; Wegg & Gerhard 2013; Bland-Hawthorn & Gerhard 2016), its extent in the y-direction is larger than in the x-direction, where x and y are measured on the Galactic plane from the Galactic centre, y in the direction Sun-Galactic centre, and x perpendicular to it. To avoid contamination from background and foreground stars, we selected stars inside |x| ≤2.5 kpc and |y| ≤ 3 kpc. The extent in the z-direction, the axis perpendicular to the Galactic plane, depends on the selected galactic latitude line of sight; it is of about 2 kpc at |b| = 10°.

The spatial structure of stellar populations with different chemical abundances in the Milky Way stellar disc is complex. Bovy et al. (2012a) showed that the Milky Way disc has a continuum of stellar mono-abundance populations, each having a simple spatial structure in both the vertical and the radial directions. Hayden et al. (2015) also showed that the stellar abundance distribution in the Galactic disc varies with cylindrical Galactocentric distance (R) and with height |z| (see their Fig. 4). In particular, stars with high-[Fe/H] values ([Fe/H] > 0 dex) are more confined to the mid-plane and to the inner disc (3 < R < 5 kpc). At least in the inner disc, the metal-poor mono-abundance populations have greater scale-heights, while the metal-rich mono-abundance populations have smaller scale-heights (Bovy et al. 2016). Because inner thin-disc stars are metal-rich and because we studied the mapping of the inner disc into the B/P bulge, we reasonably assumed that the disc stars in the simulation disc+bulge and the thin-disc stars in the simulation thin+thick may be considered as a proxy of metal-rich stars. On the other hand, bulge stars and intermediate-disc and thick-disc stars may be considered as a proxy of metal-poor stars. We calculated the difference in the mean Galactocentric radial velocity between metal-rich and metal-poor stars (Δ VGC = VGC (metal-rich) – VGC (metal-poor)) in both models at different latitudes (b = 0°, − 2°, − 4°, − 6°, − 8° and − 10°), and longitudes (l = 0°, ± 5°, ± 10° and ± 15°). The Δ VGC trends predicted by the two models are different (see next section), and when compared to observations, they may help to shed some light on the origin of moderately metal-poor stars.

3 Results

Figure 1 shows the rotation curve as deduced by radial velocity measurements in the model with a spheroid for the disc and the bulge components separately, as well as the difference in the mean Galactocentric radial velocity Δ VGC between disc and classical bulge stars. In Fig. 1, the size of the fields is Δl = Δb = 1° for b ≥ − 6° and Δ lb = 1°.5 at latitudes b < − 6°.

As shown in Di Matteo et al. (2014), the disc component displays approximately cylindrical rotation with a weak dependence of rotation speed on latitude; the velocity rotation curve is steeper towards lower latitudes. This effect has previously been noted by Howard et al. (2009) in BRAVA data. The difference in the velocity curves with latitude is similar to what is observed in ARGOS data (Ness et al. 2013a) and in GIBS data (Zoccali et al. 2014). Even though classical bulge stars acquire some angular momentumduring the bar formation and evolution (Saha et al. 2012; Saha & Gerhard 2013), their rotational velocities remain much lower than those of disc stars. This is clearly shown in Fig. 1 (right panel), where |Δ VGC| values are higher than about 50 km s−1 for |l| ≥   ~5°.

Figure 2 shows the rotational velocity for the three disc components in the thin+thick model. The modelled disc is madeof 20 000 000 particles, redistributed in the three components, compared to about 18 800 000 particles of the disc in the disc+bulge model. In order to obtain better number statistics, the adopted size of the fields was Δ l = Δb = 1°.5. We recall that thin-disc stars are concentrated towards the Galactic plane, and they constitute 40% of the total stellar density at |b| <  ~2° and contribute to less than 25% at |b|~ 8°−10° (Di Matteo 2016). The opposite trend is observed for thick-disc stars: their fraction increases with galactic latitude from about 20% at |b| < ~2° to more than 40% at |b|~ 8°−10° (Di Matteo 2016). Figure 2 shows that the hottest component (thick disc) exhibits cylindrical rotation with only weak latitude variations. The thin-disc component does not show cylindrical rotation for all latitudes, the rotation velocities being similar up to b ≥−2° and |l|≤ ~5°. The intermediate-disc component displays cylindrical rotation up to about b ≥ −4°, but the variation in rotational velocity as a function of the latitude is smaller than that corresponding to thin-disc stars. Our results are similar to those predicted by the M2M chemodynamical model of Portail et al. (2017) based on the ARGOS data and the APOGEE DR12 data: metal-rich stars (bin A, 0.5 ≥ [Fe/H] ≥ 0.0) and metal-poor stars (bin B,0.0 ≥ [Fe/H] ≥ –0.5) do not exhibit cylindrical rotation for all latitudes (see their Fig. 15).

Figure 3 shows the mean radial velocity differences Δ VGC between the modelled disc components. Intermediate-disc stars rotate slightly faster than thin-disc stars (left panel) for b < ~−4°, while for b > ~−4°, the trend is inverse, the maximum difference in |ΔVGC| being of about 10 km s−1. Portail et al. (2017) found a similar trend for stars in bins A and B and argued that it is likely the result of the complex orbital structure of the bar in the plane. Like intermediate-disc stars, thick-disc stars also display faster rotation than thin-disc stars (middle panel) only for b < ~−4°, but the maximum difference |ΔVGC| is of about 20 km s−1. Recently, Fragkoudi et al. (2017a), using N-body simulations with a thin disc and a thick disc, found that the rotational velocity of thick-disc stars in the outer parts of the bulge can be larger than those of thin-disc stars by about 20% for Milky Way-like orientation of the bar. This was interpreted as a consequence of the orbital structure of the thin-disc and thick-disc bars. The difference ΔVGC between intermediate-disc and thick-disc stars (right panel) shows that thick-disc stars rotate slightly faster than the intermediate-disc stars for almost all latitudes. Figure 1 (right panel) and Fig. 3 show that the trends of |Δ VGC| are different. In the disc+bulge model, at |l| ~ 5° |Δ VGC| is larger than about 40 km s−1, while in the thin+thick model, |Δ VGC| between the thin and the thick components is smaller than about 20 km s−1. At |l| >  ~5°, |Δ VGC| increases in both models, but only slightly in the thin+thick model compared to the disc+bulge model. As quoted before, we assume in what follows that disc stars in the disc+bulge simulation and thin-disc stars in the thin+thick simulation may be considered as a proxy of metal-rich stars. Similarly, bulge stars and intermediate-disc and thick-disc stars may be considered as a proxy of moderately metal-poor stars. The different trends observed in both simulations are compared in the next section with the kinematics of some of the available observational data of the Milky Way bulge.

thumbnail Fig. 1

Rotation curve of disc stars (left panel) and of classical bulge stars (middle panel) of the N-body model with a disc and a classical spheroid (B/D = 10%) (disc+bulge simulation). Δ VGC = VGC (disc) – VGC (classical bulge) between disc stars and classical bulge stars (right panel). Stars at |x| ≤ 2.5 kpc and |y| ≤ 3 kpc from the Galactic centre are selected. Different latitudes are shown for the modelled galaxy: b = 0°, − 2°, − 4°, − 6°, − 8°, and − 10°. The size of the fields is Δl = Δb = 1° for b ≥ − 6°. To obtain better number statistics for samples at latitudes b < −6°, Δ l = Δ b = 1. ° 5 have been adopted.

thumbnail Fig. 2

Rotation curve of thin-disc stars (left panel), intermediate-disc stars (middle panel) and thick-disc stars (right panel) of the N-body disc galaxy model (thin+thick simulation). Stars at |x| ≤ 2.5 kpc and |y| ≤ 3 kpc from the Galactic centre are selected. Different latitudes are shown for the modelled galaxy: b = 0°, − 2°, − 4°, − 6°, − 8°, and − 10°. The size of the fields is Δl = Δb = 1. °5.

thumbnail Fig. 3

ΔVGC values between stars in the different disc components in the thin+thick simulation. VGC (thin disc) – VGC(intermediate disc) (left panel); VGC(thin disc) – VGC(thick disc) (middle panel), and VGC(intermediate disc) - VGC(thick disc) (right panel). Stars at |x|≤ 2.5 kpc and |y| ≤ 3 kpc from the Galactic centre are selected. Different latitudes are shown: b = 0°, − 2°, − 4°, − 6°, − 8°, and − 10°. The size of the fields is Δl = Δb = 1. °5.

4 Comparison with observations

We compared the ΔVGC results of the two models with data from the ARGOS survey (Freeman et al. 2013) and from the APOGEE DR13 survey (Albareti et al. 2017; Majewski et al. 2017) along the lines of sight considered in this paper. Both surveys provide kinematic and metallicity measurements.

The ARGOS survey is a large medium-resolution spectroscopic survey that covers l between −15° and 15° and b between −5° and −10°. Of the metallicity components identified in the survey (A to E) (Ness et al. 2013a), A has a thin-disc origin, while C is likely to be the old-thick component that was in place before the bar formation (Di Matteo et al. 2014, 2015). We do not have individual ARGOS data but only average radial velocities for the A, B, and C components at different longitudes and latitudes selected within R ≤ 3.5 kpc (Ness et al. 2013b). The difference ΔVGC between the metal-rich component (A) and the metal-poor component (C) is shown in Fig. 4 (left panel) and is compared with ΔVGC values of the simulation with classical bulge (dotted lines, Fig. 1, right panel) and the values of the simulation with pure disc components (solid lines, Fig. 3, middle panel). ARGOS data are in good agreement with the thin+thick model, supporting the result that the metal-poor component C is associated with the thick disc (Ness et al. 2013a; Di Matteo et al. 2014, 2015; Di Matteo 2016) at the latitudes covered by the ARGOS survey. A similar conclusion has been obtained by Portail et al. (2017) for the C component outside the central kpc.

APOGEE data allow us to complement ARGOS data at latitudes near the Galactic plane. The APOGEE survey is a high-resolution spectroscopic survey in the near-infrared of Milky Way stellar populations. It provides both accurate radial velocities and reliable abundance measurements. Moreover, it allows us to observe dust-obscured regions of the Galactic bulge. We refer to Majewski et al. (2017) for a complete description of APOGEE, to Zasowski et al. (2013) for the target selection, and to García Pérez et al. (2016) for the ASPCAP: The APOGEE Stellar Parameter and Chemical Abundances Pipeline. To isolate the stars in the bulge region, we selected stars in a range of heliocentric distances between 4 and 12 kpc. Then we removed stars that have logg> 0.5 in order to reduce the metallicity bias against metal-rich stars in the inner Galaxy (Hayden et al. 2015; Ness et al. 2016). The distances were estimated using the Padova set of isochrones (Schultheis et al. 2017). For each star, the closest point on the isochrones was obtained from the stellar parameters [M/H], logg, and Teff. Finally, we kept stars with relative error in the distance smaller than 50%, although this cut does not change the results given that most of the stars have a relative error smaller than 40%. The distance selection of APOGEE stars is rather similar to the selection of ARGOS stars (distances between 4.5 and 11.5 kpc from the Sun, Ness et al. 2016), taking into account uncertainties in the star distance determinations in both surveys. Figure 4 (right panel) displays the results for four regions: one centred at l ~ 5° and b ~ −10° (110 stars), and the other three at low Galactic latitude (b ~−2°) and l ~ 5.2° (110 stars), 7.2° (91 stars) and 12.5° (137 stars), respectively. Each sample was divided into metal-rich stars ([M/H] ≥ 0) and metal-poor stars ([M/H] < 0), and [M/H] represents the global metallicity, which is closely correlated with the iron abundance. There is no systematic difference between [M/H] and [Fe/H] in the APOGEE data (Holtzman et al. 2015; García Pérez et al. 2016). Δ VGC was measured as the difference in the mean radial Galactocentric velocity of metal-rich stars minus metal-poor stars. Like the ARGOS data, the APOGEE data show a trend similar to the observed in the thin+thick model. This result suggests that most of the moderately metal-poor stars in the considered APOGEE samples are compatible with being predominantly thick-disc stars.

thumbnail Fig. 4

Comparison of ΔVGC values for the disc+bulge and the thin+thick simulations with ARGOS data: VGC (A component) – VGC (C component) (left panel) and APOGGE data VGC([M/H] >0) – VGC ([M/H] <0) (right panel). For details see the text.

5 Discussion

The kinematic trends obtained in Sect. 4 show that moderately metal-poor stars observed in the bulge region (–1.0 ≤ [Fe/H] < 0) in the ARGOS and in the APOGEE DR13 surveys are well reproduced with the co-spatial stellar discs model, which has a massive thick disc. It should be noted that Williams et al. (2016) found from the analysis of the kinematics of Gaia-ESO survey stars (Gilmore et al. 2012) “no discernible difference between the rotational signature” in the metal-rich stars ([M/H] ≥ 0) and the metal-poor ones ([M/H] < 0). On the other hand, from GIBS low spectral resolution data, Zoccali et al. (2017) showed that the difference in the Galactic radial velocity between metal-poor and metal-rich components seems marginal.

Elemental abundance trends for metal-poor stars are similar to those of local thick-disc stars. Studies based on red giants in the bulge found no evidence for a different behaviour of the alpha-elements versus iron trends in the bulge compared to the local thick disc for sub-solar metallicities (e.g. (Meléndez et al. 2008; Alves-Brito et al. 2010; Ryde et al. 2010; Gonzalez et al. 2011, 2015; Johnson et al. 2014; Jönsson et al. 2017). Bensby et al. (2017) observed in their analysis of micro-lensed dwarf and sub-giant stars that the alpha-elements (Mg, Si, Ca, and Ti) trends with [Fe/H] for metal-poor metallicities tend to follow the same trends as can be seen for local thick-disc stars and that the level of alpha-elements enhancement appears slightly higher in the bulge than in the thick disc. They argued that the observed offset would indicate that the star formation rate was slightly faster in the bulge than in the local thick disc, which is consistent with the much denser environment in the inner parts of the Galaxy. On the other hand, Rojas-Arriagada et al. (2017) obtained spectroscopic data of red clump stars from the fourth internal data release of the Gaia-ESO survey and found that for metal-poor stars, the [Mg/Fe] vs. [Fe/H] trend is comparable to the corresponding to the thick-disc sequence. Ryde et al. (2016) investigated the Galactic centre region at a projected Galactocentric distance of about 300 pc by observing [alpha/Fe] element trends of Mg and Si versus metallicity from high-resolution spectroscopic data for M giants. They found a wide range in metallicities from − 1.2 < [Fe/H] < 0.3. Their Fig. 4 shows that the obtained trends are rather similar, within the uncertainties, with the abundance trends based on micro-lensed dwarfs (Bensby et al. 2013) for the outer bulge.

Bensby et al. (2017) also analysed light odd-Z elements (Na and Al), iron-peak elements (Cr, Ni and Zn), and neutron-capture elements (Y and Ba). With the exception of the [Al/Fe] trend, which appeared to be placed at the upper envelope of the thick-disc trend, the other abundance-[Fe/H] trends display similarities between the thick disc and the metal-poor bulge. From high-resolution spectra, Van der Swaelmen et al. (2016) analysed heavy elements (Ba, La, Ce, Nd, and Eu) of bulge stars, most of which have [Fe/H] ≥ –1.0 dex. In particular, for metal-poor stars, they found that Ba and La are enhanced with respect to their thick-disc counterparts, but the observed [Eu/Fe] trend is in good agreement with the trend of Galactic thick disc derived by Bensby et al. (2005) and Reddy et al. (2006). Unlike the work of Bensby et al. (2017), which is based on homogeneous data sets analysed in a consistent manner, Van der Swaelmen et al. (2016) compared their determinations with that of the literature, which might explain the observed discrepancies.

Recently, a similar thin+thick model has allowed us to reproduce the metallicity gradient in the disc as derived by APOGEE DR13 data (Fragkoudi et al. 2017b). Moreover, the central concentration of metal-poor stars in the inner kpc of the Galaxy (e.g. Zoccali et al. 2017; Portail et al. 2017) and the metallicity distribution function of the bulge are also reproduced by the model (Fragkoudi et al. 2018).

To summarise, both kinematic and chemical results agree with the fact that moderately metal-poor stars are formed out of the thick-disc stars via secular evolution, without the necessity of a massive classical spheroid. However, this does not mean that the Milky Way bulge does not contain classical bulge stars, but that its mass contribution probably is a small percentage of the total bulge mass, as has been reported in several works (e.g. Shen et al. 2010; Kunder et al. 2012, 2016; Di Matteo et al. 2014; Debattista et al. 2017). Howes et al. (2015) found extremely metal-poor stars ([Fe/H] ≤ –2.3 dex) in the bulge. The observed stars follow tight orbits around the Galactic centre and have chemical compositions consistent with typical halo stars of the same metallicity, although they do not have the large carbon enhancements expected in halo stars. Whether most of the stars with [Fe/H] ≤ –1.0 dex are halo stars or classical bulge stars remains an open question and is beyond the scope of this paper.

6 Concluding remarks

The observed metallicity distribution function of Milky Way bulge stars show that most of them have [Fe/H] > –1 dex. We here explored the origin of moderately metal-poor stars (–1.0 ≤ [Fe/H] < 0) by comparing a kinematic signature given by the difference in the mean Galactocentric radial velocity between metal-rich and metal-poor stars in two different N-body models with ARGOS data and APOGEE DR13 data. One model consists of a disc and a small classical bulge (B/D = 0.1), and the other consists of a composite stellar disc (a kinematically cold thin-disc, an intermediate disc with intermediate kinematics, and a kinematically hot thick-disc) where the intermediate and thick discs together represent 50% of the stellar mass. As our models do not have chemical information, we reasonably assumed that the thin disc in both simulations may be considered as a proxy of metal-rich stars. Similarly, classical bulge stars and intermediate-disc and thick-disc stars may be considered as a proxy of metal-poor stars. The kinematic signature obtained at different latitudes and longitudes shows that the trends predicted by the models are different. We show that the observed kinematic trends from ARGOS data and APOGEE DR13 data are well reproduced with the co-spatial stellar discs model, which argues that the moderately metal-poor stars were formed from the thick-disc stars via secular evolution, leaving little room for a classical bulge origin. Our results give further evidence against the scenario that most of the metal-poor stars are classical bulge stars. If classical bulge stars exist, most of them probably have metallicities [Fe/H] < –1 dex, and their contribution to the mass of the bulge should be a small percentage of the total bulge mass.

Acknowledgments

This work has been supported by the ANR (Agence Nationale de la Recherche) through the MOD4Gaia project (ANR-15-CE31- 0007, P.I.: P. Di Matteo). FF is supported by a postdoctoral grant from the Centre National d’Etudes Spatiales (CNES).

References

  1. Albareti, F. D., et al. (SDSS Collaboration), 2017, ApJS, 233, 25 [NASA ADS] [CrossRef] [Google Scholar]
  2. Alves-Brito, A., Meléndez, J Asplund, M., Ramírez, I., & Young, D. 2010, A&A, 513, A35 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  3. Athanassoula, E. 2005, MNRAS, 358, 1477 [NASA ADS] [CrossRef] [Google Scholar]
  4. Babusiaux, C., Gómez, A., Hill, V., et al. 2010, A&A, 519, A77 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  5. Babusiaux, C., Katz, D., Hill, V., et al. 2014, A&A, 563, A15 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  6. Bensby, T., Feltzing, S., Lundstrom, I., & Ilyin, I., 2005, A&A, 433, 185 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  7. Bensby, T., Alves-Brito, A., Oey, M. S., Yong, D., & Meléndez J. 2011, ApJ, 735, L46 [NASA ADS] [CrossRef] [Google Scholar]
  8. Bensby, T., Yee, J. C., Feltzing, S., et al. 2013, A&A, 549, A147 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  9. Bensby, T., Feltzing, S., Gould, A., et al. 2017, A&A, 605, A89 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  10. Binney, J., & Tremaine, S. 1987, Nature, 326, 219 [Google Scholar]
  11. Bland-Hawthorn, J., & Gerhard, O. 2016, ARA&A, 54, 529 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  12. Bovy, J., Rix, H.-W., Liu, C., et al. 2012a, ApJ, 753, 148 [NASA ADS] [CrossRef] [Google Scholar]
  13. Bovy, J., Rix, H.-W., Hogg, D. W., et al. 2012b, ApJ, 755, 115 [NASA ADS] [CrossRef] [Google Scholar]
  14. Bovy, J., Rix, H.-W., Schlafly, E. F., et al. 2016, ApJ, 823, 30 [NASA ADS] [CrossRef] [Google Scholar]
  15. Cao, L., Mao, S., Nataf, D., Rattenbury, N. J., & Gould, A. 2013, MNRAS, 434, 595 [NASA ADS] [CrossRef] [Google Scholar]
  16. Combes, F., & Sanders, R. H. 1981, A&A, 96, 164 [NASA ADS] [Google Scholar]
  17. Debattista, V. P., Mayer, L., Carollo, C. M., et al. 2006, ApJ, 645, 209 [NASA ADS] [CrossRef] [Google Scholar]
  18. Debattista, V. P., Ness, M., Gonzalez, O. A., et al. 2017, MNRAS, 469, 1587 [NASA ADS] [CrossRef] [Google Scholar]
  19. De Propris, R., Rich, R. M., Kunder, A., et al. 2011, ApJ, 732, L36 [NASA ADS] [CrossRef] [Google Scholar]
  20. Di Matteo, P., Haywood, M., Gómez, A., et al. 2014, A&A, 567, A122 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  21. Di Matteo, P., Gómez, A., Haywood, M., et al. 2015, A&A, 577, A1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  22. Di Matteo, P. 2016, PASA, 33, e027 [NASA ADS] [CrossRef] [Google Scholar]
  23. Dwek, E., Arendt, R. G., Hauser, M. G., et al. 1995, ApJ, 445, 716 [CrossRef] [Google Scholar]
  24. Fragkoudi, F., Di Matteo, P., Haywood, M., et al. 2017a, A&A, 606, A47 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  25. Fragkoudi, F., Di Matteo, P., Haywood, M., et al. 2017b, A&A, 607, L4 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  26. Fragkoudi, F., Di Matteo, P., Haywood, M., et al. 2018, A&A, in press, DOI: 10.1051/0004-6361/201732509 [Google Scholar]
  27. Freeman, K., Ness, M., Wylie-de-Boer, E., et al. 2013, MNRAS, 428, 3660 [NASA ADS] [CrossRef] [Google Scholar]
  28. Fuhrmann, K., Chini, R., Hoffmeister, V. H., & Bernkopf, J. 2012, MNRAS, 420, 1423 [NASA ADS] [CrossRef] [Google Scholar]
  29. García Pérez, A. E., Allende Prieto, C., Holtzman,J. A., et al. 2016, AJ, 151, 144 [NASA ADS] [CrossRef] [Google Scholar]
  30. García Pérez, A. E., Ness, M., Robin, A. C., et al. 2018, ApJ, 852, 91 [NASA ADS] [CrossRef] [Google Scholar]
  31. Gerhard, O. 2017, IAU Symp. 334, in press, [arXiv:1710.01544] [Google Scholar]
  32. Gilmore, G., Randich, S., Asplund, M., et al. 2012, The Messenger, 147, 25 [NASA ADS] [Google Scholar]
  33. Gómez, A., Di Matteo, P., Stefanovitch, N., et al. 2016, A&A, 589, A122 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  34. Gonzalez, O. A., Rejkuba, M., Zoccali, M., et al. 2011, A&A, 530, A54 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  35. Gonzalez, O. A., Zoccali, M., Vasquez, S., et al. 2015, A&A, 584, A46 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  36. Hayden, M. R., Bovy, J., Holtzman, J. A., et al. 2015, ApJ, 808, 132 [Google Scholar]
  37. Haywood, M., Di Matteo, P., Lehnert, M. D., Katz, D., & Gómez A. 2013, A&A, 560, A109 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  38. Hill, V., Lecureur, A., Gómez, A., et al. 2011, A&A, 534, A80 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  39. Holtzman, J. A., Shetrone, M., Johnson, J. A., et al. 2015, AJ, 150, 148 [NASA ADS] [CrossRef] [Google Scholar]
  40. Howard, C. D., Rich, R. M., Reitzel, D. B., et al. 2008, ApJ, 688, 1060 [NASA ADS] [CrossRef] [Google Scholar]
  41. Howard, C. D., Rich, R. M., Clarkson, W., et al. 2009, ApJ, 702, L153 [NASA ADS] [CrossRef] [Google Scholar]
  42. Howes, L. M., Casey, A. R., Asplund, M., et al., 2015, Nature, 527, 484 [NASA ADS] [CrossRef] [Google Scholar]
  43. Johnson, C.I., Rich, R.M., Kobayashi, C., Kunder, A., & Koch, A. 2014, AJ, 148, 67 [NASA ADS] [CrossRef] [Google Scholar]
  44. Jönsson, H., Ryde, N., Schultheis, M., & Zoccali, M. 2017, A&A, 598, A101 [CrossRef] [EDP Sciences] [Google Scholar]
  45. Kunder, A., Koch, A., Rich R. M., et al. 2012, AJ, 143, 57 [NASA ADS] [CrossRef] [Google Scholar]
  46. Kunder, A., Rich, R.M., Koch, A., et al. 2016, ApJ, 821, L25 [NASA ADS] [CrossRef] [Google Scholar]
  47. Majewski, S. R., Schiavon, R. P., Frinchaboy, P. M., et al. 2017, AJ, 154, 94 [NASA ADS] [CrossRef] [Google Scholar]
  48. Martinez-Valpuesta, I., Shlosman, I., & Heller, C. 2006, ApJ, 637, 214 [NASA ADS] [CrossRef] [Google Scholar]
  49. Martinez-Valpuesta, I., & Gerhard, O. 2013, ApJ, 766, L3 [NASA ADS] [CrossRef] [Google Scholar]
  50. McWilliam, A., & Zoccali, M. 2010, ApJ, 724, 1491 [NASA ADS] [CrossRef] [Google Scholar]
  51. Meléndez, J., Asplund, M., Alves-Brito, A., et al. 2008, A&A, 484, L21 [NASA ADS] [CrossRef] [EDP Sciences] [MathSciNet] [Google Scholar]
  52. Nataf, D. M., Udalski, A., Gould, A., Fouqué, P., & Stanek, K. Z. 2010, ApJ, 721, L28 [NASA ADS] [CrossRef] [Google Scholar]
  53. Ness, M., Freeman, K., Athanassoula, E., et al. 2013a, MNRAS, 430, 836 [NASA ADS] [CrossRef] [Google Scholar]
  54. Ness, M., Freeman, K., Athanassoula, E., et al. 2013b, MNRAS, 432, 2092 [NASA ADS] [CrossRef] [Google Scholar]
  55. Ness, M., Zasowski, G., Johnson, J. A., et al. 2016, ApJ, 819, 2 [NASA ADS] [CrossRef] [Google Scholar]
  56. Pérez-Villegas, A, Portail, M., & Gerhard, O. 2017, MNRAS, 464, L80 [NASA ADS] [CrossRef] [Google Scholar]
  57. Portail, M., Wegg C., Gerhard, O., & Ness, M. 2017, MNRAS, 470, 1233 [NASA ADS] [CrossRef] [Google Scholar]
  58. Poleski, R., Udalski, A., Gould, A., et al. 2013, ApJ, 776, 76 [NASA ADS] [CrossRef] [Google Scholar]
  59. Raha, N., Sellwood, J. A., James, R. A., & Kahn, F. D. 1991, Nature, 352, 411 [NASA ADS] [CrossRef] [Google Scholar]
  60. Rangwala, N., Williams, T. B., & Stanek, K. Z. 2009, ApJ, 691, 1387 [NASA ADS] [CrossRef] [Google Scholar]
  61. Reddy, B. E., Lambert, D. L., & Allende Prieto C. 2006, MNRAS, 367, 1329 [NASA ADS] [CrossRef] [Google Scholar]
  62. Rojas-Arriagada, A., Recio-Blanco, A., Hill, V., et al. 2014, A&A, 569, A103 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  63. Rojas-Arriagada, A., Recio-Blanco, A., de Laverny, P., et al. 2017, A&A, 601, A140 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  64. Ryde, N., Gustafsson, B., Edvardsson, B., et al. 2010, A&A, 509, A20 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  65. Ryde, N., Schultheis, M., Grieco, V., et al. 2016, AJ, 151, 1 [NASA ADS] [CrossRef] [Google Scholar]
  66. Saha, K, Martinez-Valpuesta, I., & Gerhard, O. 2012, MNRAS, 421, 333 [NASA ADS] [Google Scholar]
  67. Saha, K., & Gerhard, O. 2013, MNRAS, 430, 2039 [NASA ADS] [CrossRef] [Google Scholar]
  68. Schultheis, M., Rojas-Arriagada, A., García Pérez, A. E., et al. 2017, A&A, 600, A14 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  69. Shen, J., Rich,R. M., Kormendy, J., et al. 2010, ApJ, 720, L72 [NASA ADS] [CrossRef] [Google Scholar]
  70. Snaith, O., Haywood, M., Di Matteo, P., et al. 2014, ApJ, 781, L31 [CrossRef] [Google Scholar]
  71. Uttenthaler, S., Schultheis, M., Nataf, D. M., et al. 2012, A&A, 546, A57 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  72. Van der Swaelmen, M., Barbuy, B., Hill, V., et al. 2016, A&A, 586, A1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  73. Wegg, C., & Gerhard, O. 2013, MNRAS, 435, 1874 [Google Scholar]
  74. Willams, A. A., Evans, N. W., Molloy, M., et al. 2016, ApJ, 824, L29 [NASA ADS] [CrossRef] [Google Scholar]
  75. Zasowski, G., Johnson, J. A., Frinchaboy,P. M., et al. 2013, AJ, 146, 81 [NASA ADS] [CrossRef] [Google Scholar]
  76. Zoccali, M., Gonzalez, O. A., Vasquez, S., et al. 2014, A&A, 562, A66 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  77. Zoccali, M., Vasquez, S., Gonzalez, O. A., et al. 2017, A&A, 599, A12 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]

All Figures

thumbnail Fig. 1

Rotation curve of disc stars (left panel) and of classical bulge stars (middle panel) of the N-body model with a disc and a classical spheroid (B/D = 10%) (disc+bulge simulation). Δ VGC = VGC (disc) – VGC (classical bulge) between disc stars and classical bulge stars (right panel). Stars at |x| ≤ 2.5 kpc and |y| ≤ 3 kpc from the Galactic centre are selected. Different latitudes are shown for the modelled galaxy: b = 0°, − 2°, − 4°, − 6°, − 8°, and − 10°. The size of the fields is Δl = Δb = 1° for b ≥ − 6°. To obtain better number statistics for samples at latitudes b < −6°, Δ l = Δ b = 1. ° 5 have been adopted.

In the text
thumbnail Fig. 2

Rotation curve of thin-disc stars (left panel), intermediate-disc stars (middle panel) and thick-disc stars (right panel) of the N-body disc galaxy model (thin+thick simulation). Stars at |x| ≤ 2.5 kpc and |y| ≤ 3 kpc from the Galactic centre are selected. Different latitudes are shown for the modelled galaxy: b = 0°, − 2°, − 4°, − 6°, − 8°, and − 10°. The size of the fields is Δl = Δb = 1. °5.

In the text
thumbnail Fig. 3

ΔVGC values between stars in the different disc components in the thin+thick simulation. VGC (thin disc) – VGC(intermediate disc) (left panel); VGC(thin disc) – VGC(thick disc) (middle panel), and VGC(intermediate disc) - VGC(thick disc) (right panel). Stars at |x|≤ 2.5 kpc and |y| ≤ 3 kpc from the Galactic centre are selected. Different latitudes are shown: b = 0°, − 2°, − 4°, − 6°, − 8°, and − 10°. The size of the fields is Δl = Δb = 1. °5.

In the text
thumbnail Fig. 4

Comparison of ΔVGC values for the disc+bulge and the thin+thick simulations with ARGOS data: VGC (A component) – VGC (C component) (left panel) and APOGGE data VGC([M/H] >0) – VGC ([M/H] <0) (right panel). For details see the text.

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.