Free Access
Issue
A&A
Volume 591, July 2016
Article Number A37
Number of page(s) 8
Section Galactic structure, stellar clusters and populations
DOI https://doi.org/10.1051/0004-6361/201527654
Published online 07 June 2016

© ESO, 2016

1. Introduction

Ever since the seminal work of Janes (1979), open clusters have been used to trace the metallicity distribution in the Galactic disk. In that work, Janes found that the distribution of metallicity ([Fe/H]) with Galactocentric distance (RGC) was consistent with a gradient with metallicity decreasing with increasing distance from the Galactic center at a rate of 0.05 dex kpc-1. Such a distribution has been shown in chemical evolution models to represent the inside-out growth of the Galactic disk (e.g., Chiappini et al. 1997). The term “gradient” has come to be used as short hand to refer to the (azimuthally averaged) change of [Fe/H] as a function of RGC and has generally been described as a linear function, as in Janes (1979). However, it has been pointed out several times over the decades that the metallicity distribution with RGC shown by open clusters and other stellar populations can be described by other functions, as well (e.g., Twarog et al. 1997; Andrievsky et al. 2002). Indeed, as the number of clusters studied by high-resolution spectroscopy has grown over the past decade, it now seems clear that the form of the radial dependence on mean metallicity in the disk changes with Galactocentric radius. While a linearly decreasing function may well describe the change in [Fe/H] as a function of RGC through the solar neighborhood out to RGC ~ 10−13 kpc, the clusters in the outer disk have roughly constant [Fe/H] ~ −0.3 and remarkably small dispersion (e.g., Heiter et al. 2014).

Much work has been done in recent years to better characterize the gradient in the outer disk, and the nature of the transition between the inner and outer disk (e.g., Carraro et al. 2007; Yong et al. 2012; Bragaglia et al. 2008; Jacobson et al. 2011; Cantat-Gaudin et al. 2016). It can be said, however, that the gradient in the inner disk has received comparatively less attention. This is not surprising given that the line of sight toward the Galactic center is much more crowded, confused, and extincted. Nevertheless, this is an important part of the Galaxy where the bulge/bar meets both the thin and thick disks.

There is some evidence that the nature of the gradient changes inside the solar circle. Some of the first evidence came from Cepheids: Andrievsky et al. (2002) found that the gradient inside RGC ~ 7−8 kpc was steeper than that outside the solar circle. Recent work by Genovali et al. (2013, 2014) has reinforced the increase in [Fe/H] seen in the inner disk to super solar metallicities of ~0.4 dex. The Cepheid sample of Martin et al. (2015), which extends to as close as 3 kpc from the Galactic center, indicates that 0.4 dex is the upper bound of the metallicity range in the inner disk, and the gradient plateaus at this level inside of RGC = 5 kpc. In contrast to the evidence from Cepheids, a study of inner disk red clump stars has found no evidence for a gradient at all: red clump stars spanning RGC ~ 4.5−6.5 kpc along a line of sight towards the Galactic center have a mean [Fe/H] = + 0.17 with no sign of a change with RGC (Hill et al. 2012).

Regarding the open cluster population, the number of confirmed clusters inside RGC ~ 7 kpc that we could use to probe the inner disk gradient is small relative to that in the solar neighborhood and beyond. Many works have noted the dearth of old (>1 Gyr) clusters in the inner region (see, e.g., Bonatto et al. 2006) and while new surveys in the infrared are revealing promising candidates for future study (e.g., Froebrich et al. 2007; Borissova et al. 2014), many of these candidate clusters are embedded and highly obscured. Relatively few have fundamental parameters determined, and of those even fewer are accessible to traditional optical techniques for abundance analysis. However, there have been efforts to study the known open clusters in the inner disk specifically to determine whether the older population of open clusters shows the same steepening of the gradient that is seen in the younger Cepheids. With their sample of three clusters all inside RGC ~ 7 kpc, along with a selection of clusters from the literature, Magrini et al. (2010) did indeed find evidence that the gradient rises more sharply towards the galactic center. While this work greatly increased the number of inner disk open clusters subject to detailed chemical abundance study, the inner disk cluster sample that has been well studied is still small relative to that in the solar neighborhood and beyond.

In the last few years, a number of spectroscopic surveys have been planned and begun, aimed at investigating and constraining in great detail the chemo-dynamical evolution of the Milky Way: for example, the extentions to SDSS (York et al. 2000) SEGUE (Yanny et al. 2009) and APOGEE (Majewski et al. 2015); RAVE (Steinmetz et al. 2006); the Gaia-ESO Survey (Gilmore et al. 2012; Randich et al. 2013); GALAH (De Silva et al. 2015). Several multi-object spectrographs have been designed and/or built that will also carry out their own spectroscopic surveys: LAMOST (Zhao et al. 2006; Cui et al. 2012); WEAVE (Dalton et al. 2012; Balcells 2014); MOONS (Cirasuolo et al. 2011); 4MOST (de Jong et al. 2014). In particular, the SDSS III APOGEE Survey (Ahn et al. 2014; Majewski et al. 2015) has done much to increase our picture of the metallicity distribution of the disk, based on near-infrared high-resolution spectroscopy of red giant stars. While the APOGEE Survey includes many open clusters, their sample predominantly has objects located at RGC ≳ 7 kpc (Frinchaboy et al. 2013; Cunha et al. 2016), and in some cases only a few members per cluster were observed. Instead, insight into the inner disk metallicity distribution comes from a study of APOGEE field stars, which found that the gradient of ~−0.09 dex kpc-1 in the Solar neighborhood transitions around RGC = 6 kpc, becoming shallower interior to that value (Hayden et al. 2014, 2015).

The Gaia-ESO Survey (GES) is a large-scale, high-resolution spectroscopic survey being carried out with the FLAMES instrument on the VLT. In addition to probing field stars in all major stellar components of the Milky Way, the survey includes a large number of open clusters spanning a wide range in age and location in the Galactic disk (Gilmore et al. 2012, Randich et al. 2013). Stars of all evolutionary states are observed in GES, from pre-main sequence stars to red giants, and elements from all groups of the Periodic Table, from Li to the neutron-capture species, are measurable in the FLAMES spectra. Regarding open clusters, with the complementary capabilities of the GIRAFFE and UVES spectographs, GES stellar samples can range from 102 to 103 members per cluster (e.g., Donati et al. 2014b; Cantat-Gaudin et al. 2014; Frasca et al. 2015; Prisinzano et al. 2016) allowing not only for robust average measures of cluster chemistry and dynamics, but also their variations as a function of stellar evolutionary state (Smiljanic et al. 2016).

Though the survey is far from complete, the first internal data releases from GES (data releases iDR1, iDR2/3, iDR4) span roughly 2.5 out of five years of observations and provide an opportunity to explore the form of the abundance gradient in the inner disk from high-resolution spectroscopy of both open clusters and field stars. Complementary studies of the thin and thick disk gradients and age-metallicity relationships using GES iDR1 data for field stars have appeared in Bergemann et al. (2014), Recio-Blanco et al. (2014) and Mikolaitis et al. (2014). Bergemann et al. (2014) used UVES data of field FGK stars to explore the age-metallicity relationships and gradients in the thin and thick disk populations. To obtain ages, they limited their sample to unevolved stars and so probed only the solar neighborhood at distances of 7 <RGC< 9 kpc. Recio-Blanco et al. (2014) used the lower resolution GIRAFFE spectra of FGK stars in the first data release to study the kinematic and spatial gradients of the thin and thick disk populations over the range 5 <RGC< 12 kpc. Mikolaitis et al. (2014) used a subset of the Recio-Blanco et al. (2014) sample and the distances calculated therein that spanned a range of 4 <RGC< 12 kpc and found a rather shallow metallicity gradient for the thin disk (~−0.03 to −0.04 dex kpc-1).

Table 1

Cluster parameters.

While current sample sizes of field stars in spectroscopic surveys such as GES and APOGEE reveal much about the state of our Galaxy, open clusters within these surveys are complementary tracers of the chemistry of the thin disk because their ages and distances can be precisely determined, they span a wide range in age and location within the Galaxy, and their chemical abundances can be robustly measured from large samples of stellar members. In this work, we utilize the iDR4 GES results obtained from UVES observations of stars in 12 intermediate-aged (>100 Myr, <3 Gyr) open clusters. These 12 are from a total of 26 open clusters included in iDR4. These results are placed in the context of those of other inner disk populations to explore the nature of the abundance gradient in this relatively understudied region of the Galaxy. This paper is organized as follows: Sect. 2 presents the GES observations used in this work; Sect. 3 presents the discussion of the inner disk metallicity gradient and Sect. 4 compares it to chemical evolution models. A summary and conclusions are available in Sect. 5.

2. GES data

The Gaia-ESO survey began on 31 December 2011. The survey makes use of the FLAMES instrument with both GIRAFFE and the fiber link to UVES. Here, we used only the UVES spectra obtained with the 580 nm setup, which have spectral resolution Rλ/ Δλ = 47 000, signal-to-noise (S/N) ≳ 50 and span 4760–6840 Å. UVES data were processed with the GES pipeline as described in Sacco et al. (2014). Stellar parameters and abundances for FGK-type stars observed with UVES were determined by up to 14 independent groups, and homogenized in a process described in Smiljanic et al. (2014). Results from working groups specializing in different stellar types were then combined to result in a set of recommended parameters utilized by the GES collaboration (Hourihane et al., in prep.). These recommended parameters, available in the “recommendedastroanalysis” table within GESiDR4Final, were used in this work.

Analysis has now been completed and the data products for all survey observations taken before July 2014 have been released to the consortium as GESviDR4Final. This data release included some updated analysis methods and homogenization across multiple working groups. As a result, this release included newly determined stellar parameters and abundances for the previously analyzed objects, so the parameters presented here differ slightly from those in earlier GES publications.

Twelve intermediate-aged clusters inside the solar circle (for this work we adopted RGC = 8.0 kpc for the Sun) have been targeted with UVES observations and included in GESviDR4Final (Table 1). Of these, results from earlier internal data releases were published for the clusters Trumpler 20, NGC 4815, NGC 6705 (M 11) and Berkeley 81. References for these works are given in Table 1. For these objects, we adopted the ages and distances determined from those published papers. For five of the remaining eight clusters we used values from the literature as given in Table 1. Three of these clusters are nearby and have distances that rest on Hipparcos measurements (NGC 2516, NGC 3532, and NGC 6633). Two others (NGC 6005 and Pismis 18) have parameters that rely on traditional methods of fitting isochrones, but with models that are consistent with the metallicities derived here; we explore the impact of errors on their distances below. For the remaining three clusters, we estimated new cluster ages, distances, and E(BV) values taking into account the new GES metallicities, which are significantly higher than previously assumed metallicities (Be 44 and NGC 6802), or for which we already had analysis underway (Trumpler 23, Overbeek et al., in prep.). Analysis for these clusters was intended only to give preliminary values that primarily take into account the revised metallicity and membership information from GES; a more complete treatment and redetermination of cluster parameters will be presented in future publications (e.g., Tang et al., in prep. for NGC 6802). For these clusters, we used radial velocities measured from both UVES and GIRAFFE observations in the cluster field to determine estimates of cluster systemic velocities. First, a mean radial velocity was calculated from all stars in the cluster field. Probable radial velocity members were chosen as those stars within approximately 2 sigma about the mean velocity. Using this indication of membership in conjuction with published photometry (Janes & Hoq 2011, for Be 44 and NGC 6802 and Carraro et al. 2006, for Trumpler 23), we fit PARSEC isochrones (Bressan et al. 2012) at the derived GES metallicities (see below), to obtain estimates of cluster reddening, distance, and age which are summarized in Table 1.

Table 2

Parameters for cluster members.

For clusters available in earlier GES releases, we determined membership independently of those previously published. Not only have the recommended parameters and abundances varied slightly with each data release as a consequence of the calibration process, but in the case of some clusters, additional fields were observed and include more cluster members than were available in iDR1 and iDR2/3. Mean cluster metallicities were based on an evaluation of individual star cluster membership using radial velocities and abundances from UVES spectra in GESviDR4Final. We used the radial velocities in combination with the derived [Fe/H] values from UVES observations as a guide for determining membership; cluster members clump tightly in a diagram of radial velocity versus [Fe/H], aiding the elimination of potential non-member stars near the cluster velocity. The majority of [Fe/H] values came from evolved stars in the red clump or along the giant branch, but for the closest clusters NGC 2516, 3532 and 6633, include determinations from dwarf stars as well. The resulting mean cluster abundances and standard deviations about the mean and radial velocities, based on these likely members, are given in Table 1. The GES ID, coordinates, and the GES recommended parameters for all stars considered cluster members and used in forming the mean metallicities are given in Table 2.

The location of these inner disk GES clusters relative to the plane of the Milky Way disk is shown in the top panel of Fig. 1. All clusters are within ~150 pc of the plane, and appear to be thin disk objects. The bottom panel of Fig. 1 shows the range of cluster age as a function of RGC. Clusters at similar RGC can vary in age by as much as 1 Gyr.

thumbnail Fig. 1

Top panel: locations of the GES open clusters (red triangles) relative to the Galactic plane (z = 0, solid line), as a function of their Galactocentric radii. All clusters lie within 150 pc of the mid-plane, and thus are members of the thin disk. Bottom panel: cluster ages (in Gyr) plotted as a function of RGC.

thumbnail Fig. 2

Metallicity gradient as shown by the GES open clusters (red triangles) in comparison to literature studies. The open cluster literature compilation of Netopil et al. (2016) is shown as open gray squares, while clusters observed by the APOGEE survey (Frinchaboy et al. 2013) are filled gray squares. The solid blue line shows the metallicity gradient as determined by APOGEE field giant stars within 250 pc of the Galactic mid-plane (Hayden et al. 2014). The red dashed line is least-squares fit to the GES sample, with the value of the slope given in the caption. [Fe/H] = 0 is indicated by the solid black line, to guide the eye.

3. Inner disk metallicity gradient

With nine open clusters inside RGC = 7 kpc, the GES sample is a factor of three times larger than previous homogeneous samples in this part of the Galactic disk (Magrini et al. 2010) and allows us to see more clearly the trend of metallicity inside the solar radius. As shown in Fig. 2, the sample of 12 GES clusters (in red triangles) shows a steady decrease in [Fe/H] with increasing Galactocentric distance over the range of RGC ~ 5.5−8 kpc. The vertical error bars represent the standard deviations about the mean abundance determined from those stars judged to be members, and given in Table 1. As discussed in Magrini et al. (2014) and Magrini et al. (2015), in light of the highly uniform abundance scale for the GES determinations and the internal chemical homogeneity of the clusters, the differences in metallicity seen among the GES clusters are reliable and significant. The current sample reinforces this fact; as can be seen from the stellar parameters given in Table 2, typical uncertainties in [Fe/H] for individual stars are 0.10 dex. The rms deviations about the cluster means, on the other hand, are much smaller, from 0.02 to 0.06 dex, indicating that there is no evidence of an internal spread in abundance within each cluster. There is a clear decrease in mean abundance from the innermost point sampled, at RGC = 5.5 kpc, out to the clusters close to the solar neighborhood. A linear regression analysis on the GES sample shows a metallicity gradient with a slope of −0.10 ± 0.02 dex kpc-1.

In addition to this clear trend, however, there is significant variation in [Fe/H] around RGC = 7 kpc. The 5 clusters at RGC = 7 kpc cover a range of 0.2 dex, from Be 44 at [Fe/H] = + 0.17 ± 0.04 to NGC 4815, with [Fe/H] = −0.03 ± 0.06, well in excess of the scatter about the mean cluster metallicities. Since stars in these two clusters have similar atmospheric parameters and were analyzed homogeneously, this dispersion in [Fe/H] is likely real. Of the five clusters around RGC = 7 kpc, two have l~ 54°, and the other three cluster around l~ 305° (Table 1). However these two groups do not also cluster together in metallicity: the [Fe/H] dispersion at each Galactic longitude is 0.07 and 0.13 dex, respectively.

As mentioned previously, this homogeneous GES cluster sample explores the inner part of the disk with considerably larger statistics than in previous studies. More clusters will be available in future GES releases, sampling also the outer parts of the disk. At this stage, in order to study the behavior of the metallicity distribution across a wider range of Galactocentric distances, we also consider the large compilation of high-resolution open cluster metallicities in Netopil et al. (2016). This compilation is an update to that of Heiter et al. (2014), who determined weighted average [Fe/H] values for clusters based on rigorous selection criteria over a limited range of atmospheric parameters of literature measures. Though the numbers of open clusters observed homogeneously as part of large spectroscopic surveys (e.g., in GES; APOGEE – Frinchaboy et al. 2013; LAMOST – Zhang et al. 2015) is growing, the Heiter et al. (2014) sample and its update is currently the closest one can come to a uniform sample formed from inhomogeneous literature measurements. Eight clusters in the current work are also in Netopil et al. (2016); we use the [Fe/H], age and RGC values given in Table 1 for the clusters in common.

Figure 2 shows the metallicity distribution with RGC for both the GES inner disk clusters (red triangles), and the Netopil et al. (2016) sample (open gray squares). We note that though clusters beyond RGC = 14 kpc have been studied, we limit our discussion to objects inside this limit, as the gradient appears to change around this Galactocentric radius (e.g., Yong et al. 2005, 2012; Carraro et al. 2007; Pancino et al. 2010; Andreuzzi et al. 2011; Donati et al. 2015; Netopil et al. 2016; see also Twarog et al. 1997). For comparison, we also plot open clusters observed as part of the APOGEE Survey (Frinchaboy et al. 2013; filled gray squares)1.

The linear gradient traced by the GES clusters, −0.10 ± 0.02 dex kpc-1 is shown as the dashed red line. The gradient described by the APOGEE cluster sample is just as steep, if not steeper; Frinchaboy et al. (2013) reported −0.20 ± 0.08 dex kpc-1 for clusters inside RGC = 10 kpc 2. However, these cluster samples overlap in a very limited region of RGC ~ 7−8 kpc, and given the small sample sizes these slopes may not fully reflect the metallicity distribution in the inner disk. Fitting their larger literature sample Netopil et al. (2016) find a linear metallicity gradient of −0.085 ± 0.017 dex kpc-1 for the 64 clusters inside 9 kpc, a slope quite consistent with the relationship traced by the GES clusters. They note the strong influence of the three clusters with much higher metallicity at RGC ~ 7 kpc (NGC 6253, 6583, and 6791) and omitting those find a more modest slope of −0.061 ± 0.015 dex kpc-1.

We note that, while there is overlap in the distribution of metallicities of the various samples, the GES abundances appear on the lower boundary of the distribution of Netopil et al. (2016) and APOGEE values at the same RGC. The slightly lower [Fe/H] values for the GES clusters around RGC ~ 8 kpc may indicate a ~0.1 dex zero-point offset between the GES and Netopil et al. (2016) metallicity scales. It is also possible that the small numbers of the GES clusters in the solar neighborhood shown here do not fully reflect its intrinsic abundance distribution, and a larger sample of clusters would include objects of solar metallicity. For the eight clusters in common between the GES and the Netopil et al. (2016) study, the difference in [Fe/H] is only 0.04 ± 0.05 (s.d.), suggesting that the systematic differences are small. Given these caveats, the metallicity distribution of the GES clusters is consistent with that of the Netopil et al. (2016) sample.

Also shown in Fig. 2 is the gradient determined from 2 × 104 giant stars in the APOGEE Survey (Hayden et al. 2014). Specifically, we plot the gradient determined from a subsample of APOGEE stars with low [α/Fe] measures, confined to within 250 pc of the Galactic mid-plane (i.e., likely thin disk giant stars). As mentioned in the Introduction, Hayden et al. (2014) found that the metallicity gradient of 0.09 dex kpc-1 broke and became shallower inside RGC = 6 kpc. Only one GES object lies inside (though on the edge of) this break, and its [Fe/H] value is consistent with both slopes shown by the APOGEE data.

Hayden et al. (2014) also traced the vertical abundance gradients in the disk, finding that they may be steeper in the inner Galaxy, with values of 0.4 dex kpc-1 for the low [α/Fe] stars inside RGC = 7 kpc. Figure 1 shows that the GES clusters reach at most 150 pc from the Galactic plane, so the impact of this gradient, if it applied to the cluster population, would be at most ~0.05 dex. Correcting for this small effect, though, would serve to steepen the gradient by a small amount, to 0.12 dex kpc-1 since the two innermost clusters are also among those most distant from the Galactic plane.

While the importance of homogeneous [Fe/H] measures to the determination of the metallicity gradient cannot be overstated, distances must also be homogeneous. Given that the GES clusters and the Netopil et al. (2016) sample have different distance scales (though both adopt RGC = 8.0 kpc for the Sun), we calculate the GES cluster metallicity gradient again adopting Netopil et al. (2016)RGC values for objects in common. This is shown in Fig. 3. The gradient steepens slightly to −0.11 ± 0.02 dex kpc-1, as the Netopil et al. (2016)RGC values are generally closer to the Sun than our values. We noted earlier that two of the clusters in the GES sample have larger than typical uncertainties in distance; these have little impact on the derived linear gradient. Changing their distances by the quoted errors changes the gradient by at most 0.01. We conclude that, in general, the metallicity gradients shown by the GES sample, the Netopil et al. (2016) sample, and the APOGEE field star sample are consistent with one another.

thumbnail Fig. 3

Same as in Fig. 2, but with distances determined by Netopil et al. (2016) for GES clusters in their compilation. Note the magnitude of the GES cluster gradient has changed by 0.01 dex kpc-1.

The motivation of this work was to determine if the GES inner disk objects showed an indication of the metallicity gradient steepening inside the solar circle as previous studies have indicated (e.g., Genovali et al. 2014; Magrini et al. 2010). Given the consistency of the slopes of the different samples in Fig. 2 and Fig. 3, there is no indication that the open clusters’ metallicity gradient changes within the range of RGC ~ 6−8 kpc, and in fact the APOGEE field star data indicate a shallower slope in the innermost region studied.

4. Comparison to model predictions

Open clusters have long been used to help constrain chemical evolution models of the Milky Way disk, through both the shape of their radial gradients and their evolution with time (e.g., Tosi 1982; Chiappini et al. 2001; Cescutti et al. 2007; Magrini et al. 2009). To date, the observational evidence for the time evolution of the gradient has been mixed. Generally, open cluster studies have found relatively little evidence of a change of the gradient with time, though there are some indications that the younger clusters have a slightly shallower distribution than the older clusters (a flattening of the gradient with time: Friel et al. 2002; Magrini et al. 2009; Carrera & Pancino 2011; Andreuzzi et al. 2011). However, studies of other objects or combinations of different stellar populations have indicated the opposite. Stanghellini & Haywood (2010), for example, found evidence that the gradient has grown steeper with time based on an analysis of disk planetary nebulae.

thumbnail Fig. 4

GES open clusters (black-lined triangles) and Netopil et al. (2016) clusters (filled squares), color-coded by age. The colored solid lines represent the metallicity gradients for different-aged stars in the chemical evolution model of Minchev et al. (2014). The gradient shown by the clusters is much steeper than that of the model.

The chemical evolution model of Minchev et al. (2013, 2014) allows for the movement of star particles from their birthplace due to radial migration effects. (We take their work as an example; other dynamical models of Milky Way-like galaxies that trace stellar chemistry do exist, e.g., Roškar et al. 2008; Schönrich & Binney 2009; Loebman et al. 2011; Bird et al. 2012; Kubryk et al. 2013.) Figure 4 shows the metallicity gradients predicted by their model for different-aged objects within 250 pc of the Galactic plane (taken from Fig. 10 and Table 1 of Minchev et al. 2014). We also plot the GES and Netopil et al. (2016) open cluster samples color-coded with the same age ranges as the model gradients. The majority of open clusters are <2 Gyr old (magenta symbols), and show a steeper gradient than the model at all age ranges. We note also that the open clusters, almost regardless of age, do not align with the model predictions, either in the value of the metallicity or the dependance on age. It has been shown that open cluster metallicity does not correlate with cluster age as would arise from a general increase in the metallicity of the disk with time (e.g., Friel et al. 2002) and this general impression is reinforced in Fig. 4; almost all the younger clusters, with ages <2 Gyr, have lower metallicities than predicted in the model. Conversely, the oldest of the clusters lie at higher metallicities than predicted.

Netopil et al. (2016), investigating possible age-metallicity relationships among their homogenized literature sample, found evidence, in fact, of an increase in metallicity with age. Comparing metallicities of clusters with ages less than 0.5 Gyr to those with ages between 1 and 2.5 Gyr, corrected for the radial gradient, they found that the younger group had metallicities ~0.07 dex lower than the older group. This held even when they only considered clusters within 7 ≤ RGC ≤ 9 kpc. We can carry out the same exercise with the GES clusters, and find an identical result. Residuals from the linear gradient show a clear correlation with cluster age (Fig. 5); the younger clusters are at systematically lower metallicities, while the older ones are at systematically higher metallicities. Using the same age groupings adopted by Netopil et al. (2016), we find, as they did, that clusters with ages less than 0.5 Gyr show a mean residual of −0.021 ± 0.026 while clusters older than 1.0 Gyr have a mean residual of + 0.041 ± 0.037. That this trend is evident in our homogeneously analyzed (and smaller) cluster sample that extends to RGC ~ 5.5 kpc is worth emphasizing.

As Netopil et al. (2016) point out, this increase in metallicity with age, contrary to what would be expected from simple expectations of chemical evolution, can be explained by the effect of radial migration, as more metal-rich and older objects from the inner disk move outward. An analogous result has been seen in a comparison of metallicity gradients traced by planetary nebulae and H II regions in M 31: the older nebulae exhibit a flatter gradient and higher oxygen abundances than the younger H II regions (Magrini et al. 2016). In their comparison to the chemical evolution model of Minchev et al. (2013), Netopil et al. (2016) found that the trend of increasing metallicity with age was qualitatively consistent with the effect of radial migration, but was stronger than predicted by the model. The same holds true for our cluster sample.

thumbnail Fig. 5

Residuals in [Fe/H] from the linear gradient as a function of cluster age.

However, in the larger context, the exact magnitude and effect of radial migration on the chemical abundance trends seen in the Milky Way have yet to be firmly established. Observationally, evidence appears to be growing that radial migration effects cannot be ignored. The lack of an age-metallicity relation for stars in the solar neighborhood can be naturally explained by migration processes, for example (e.g., Haywood 2006). Also, the shapes of metallicity distribution functions in different radial bins of APOGEE red giant stars led Hayden et al. (2015) to conclude that “migration is of global importance in the evolution of the disk” (p. 12). Furthermore, some radial migration models can reproduce observations remarkably well (e.g., Loebman et al. 2016). However, radial migration models still vary from one to another on the assumption of initial conditions, as well as how the bar and spiral arms system are treated. Still more assumptions are made when a chemical evolution model is applied to a dynamical one (Minchev et al. 2013, are clear on the sensitivity of their results to assumptions made in their model). Adding chemistry to a dynamical model introduces additional uncertainties that complicate comparison of models to one another and to observations.

As a result, different models can produce conflicting predictions. For example, Minchev et al. (2013) found that the Sun most likely formed near somewhere between RGC = 5.5−7 kpc, while Martínez-Barbosa et al. (2015) found the probability that the Sun formed inside RGC = 8 kpc to be no higher than 30%. Furthermore, they found that when radial migration effects were strong, the Sun was more likely to migrate from the outer disk, and this only occurred in specific cases.

In an analysis of FGK stars near the Sun, Haywood et al. (2013) argue for the lack of a detectable influence of radial migration, at least in the solar neighborhood. They argue that the ellipticity of stellar orbits, plus the Sun’s nearness to the transition to the metal-poor outer disk (RGC ~ 10 kpc), are sufficient to explain the presence of both metal-poor and metal-rich stars in the solar neighborhood. They furthermore argue that the Milky Way disk likely formed outside-in, rather than inside-out, based on the different age distributions of the metal-poor and metal-rich stars near the Sun.

Finally, it is still uncertain whether the effect of radial migration on open clusters is the same as that for field stars. Migration has been proposed to explain the presence of clusters with [Fe/H] ~ + 0.4 dex near RGC ~ 7 kpc, such as NGC 6791; however, a study by Jílková et al. (2012) concluded it was very unlikely. On the other hand, Fujii & Baba (2012) used N-body simulations of clusters in a Galactic disk and found a timescale of ~100 Myr for the radial migration of open clusters. In that time period, they demonstrated a cluster could move ~1.5 kpc from its birth location, with minimal mass loss. However, their simulation did not include the effects of the bar, and so it is not clear whether their general conclusions apply to clusters inside the Solar circle, as we consider here. Martinez-Medina et al. (2016) demonstrated that spiral arms can lift open clusters >1 kpc above the Galactic mid-plane, which can explain the tendency of older (>1 Gyr) clusters to have large z-distances, but they do not explore variations of cluster RGC in their model.

In conclusion, to fully gauge the impact of radial migration on the formation and evolution of the disk, the full range of observations of different stellar populations at different locations in the disk (correcting for observational biases) must be compared with multiple chemo-dynamical models that are built upon different assumptions and have different uncertainties. This should be possible within the next decade, after the many ongoing and planned spectroscopic surveys of the Galaxy are complete.

5. Summary

We have presented an analysis of the inner disk metallicity gradient based on metallicities of 12 intermediate-aged clusters observed as part of the Gaia-ESO Survey. Our sample triples the number of clusters inside RGC = 7 kpc compared to previous homogeneous literature studies, and also spans a wider range in age (~1.5 Gyr) than other studies of objects in this part of the disk (e.g., Cepheids in Andrievsky et al. 2002).

The metallicity gradient of our sample, which spans 5.5 <RGC< 8 kpc, is −0.10 ± 0.02 dex kpc-1, consistent within uncertainties to the gradient measured by APOGEE red giant stars (Hayden et al. 2014). As our sample stars all lie within 150 pc of the Galactic plane, the impact of a vertical metallicity gradient on our results is expected to be less than 0.05 dex, and would only change the radial gradient by 0.01 dex kpc-1. The metallicity gradient shown by our inner disk sample is also consistent with that shown by the larger cluster sample of Netopil et al. (2016) in the range of 6 <RGC< 14 kpc. Therefore, the GES clusters do not support previous claims in the literature that the metallicity gradient steepens inside the solar circle.

We have also found that the GES clusters exhibit a trend of increasing metallicity with cluster age: after correcting for the effect of the gradient, clusters older than 1 Gyr are 0.06 dex more metal-rich than clusters aged <0.5 Gyr. Two clusters that exemplify this effect, Berkeley 44 (age 1.6 Gyr, [Fe/H] = + 0.17 ± 0.04) and NGC 4815 (age 570 Myr, [Fe/H] = −0.03 ± 0.06), both reside at RGC = 6.9 kpc, indicating that the dispersion in [Fe/H] with RGC is ~0.2 dex.

A comparison of the inner disk metallicity gradient to the predictions of the chemical evolution model of Minchev et al. (2014) found relatively poor agreement (Fig. 4). The age-metallicity relation for our sample is consistent with the effects of radial migration predicted by Minchev et al. (2014), as also discussed in Netopil et al. (2016), but the slope and temporal evolution of the gradient predicted by the model is at odds with the values shown by open clusters in general. We look forward to a comparison of cluster and model results when the Gaia-ESO Survey is complete.


1

We show only those APOGEE clusters in which more than one star was observed.

2

Frinchaboy et al. (2013) adopted RGC = 8.5 kpc for the Sun, as opposed to 8.0 kpc in this work.

Acknowledgments

The authors thank the referee for comments that improved the presentation of this research. This work was partly supported by the European Union FP7 programme through ERC grant number 320360 and by the Leverhulme Trust through grant RPG-2012-541. We acknowledge the support from INAF and Ministero dell’ Istruzione, dell’ Università e della Ricerca (MIUR) in the form of the grant “Premiale VLT 2012” and “The Chemical and Dynamical Evolution of the Milky Way and Local Group Galaxies” (prot. 2010LY5N2T). The results presented here benefit from discussions held during the Gaia-ESO workshops and conferences supported by the ESF (European Science Foundation) through the GREAT Research Network Programme. F.J.E. acknowledges financial support from the ARCHES project (7th Framework of the European Union, n 313146). S.V. gratefully acknowledges the support provided by Fondecyt reg. 1130721. U.H. acknowledges support from the Swedish National Space Board (SNSB). D.G. gratefully acknowledges support from the Chilean BASAL Centro de Excelencia en Astrofísica y Tecnologías Afines (CATA) grant PFB-06/2007.

References

  1. Ahn, C. P., Alexandroff, R., Allen de Prieto, C., et al. 2014, ApJS, 211, 17 [NASA ADS] [CrossRef] [Google Scholar]
  2. Andreuzzi, G., Bragaglia, A., Tosi, M., & Marconi, G. 2011, MNRAS, 412, 1265 [NASA ADS] [Google Scholar]
  3. Andrievsky, S. M., Bersier, D., Kovtyukh, V. V., et al. 2002, A&A, 384, 140 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  4. Balcells, M. 2014, EAS Pub. Ser., 67, 227 [CrossRef] [EDP Sciences] [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. Bird, J. C., Kazantzidis, S., & Weinberg, D. H. 2012, MNRAS, 420, 913 [NASA ADS] [CrossRef] [Google Scholar]
  7. Bonatto, C., Kerber, L. O., Bica, E., & Santiago, B. X. 2006, A&A, 446, 121 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  8. Borissova, J., Chené, A.-N., Ramírez Alegría, S., et al. 2014, A&A, 569, A24 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  9. Bragaglia, A., Sestito, P., Villanova, S., et al. 2008, A&A, 480, 79 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  10. Bressan, A., Marigo, P., Girardi, L., et al. 2012, MNRAS, 427, 127 [NASA ADS] [CrossRef] [Google Scholar]
  11. Cantat-Gaudin, T., Vallenari, A., Zaggia, S., et al. 2014, A&A, 569, A17 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  12. Cantat-Gaudin, T., Donati, P., Vallenari, A., et al. 2016, A&A, 588, A120 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  13. Carraro, G., Janes, K. A., Costa, E., & Méndez, R. A. 2006, MNRAS, 368, 1078 [NASA ADS] [CrossRef] [Google Scholar]
  14. Carraro, G., Geisler, D., Villanova, S., Frinchaboy, P. M., & Majewski, S. R. 2007, A&A, 476, 217 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  15. Carrera, R., & Pancino, E. 2011, A&A, 535, A30 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  16. Cescutti, G., Matteucci, F., François, P., & Chiappini, C. 2007, A&A, 462, 943 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  17. Chiappini, C., Matteucci, F., & Gratton, R. 1997, ApJ, 477, 765 [NASA ADS] [CrossRef] [Google Scholar]
  18. Chiappini, C., Matteucci, F., & Romano, D. 2001, ApJ, 554, 1044 [NASA ADS] [CrossRef] [Google Scholar]
  19. Cirasuolo, M., Afonso, J., Bender, R., et al. 2011, The Messenger, 145, 11 [Google Scholar]
  20. Clem, J. L., Landolt, A. U., Hoard, D. W., & Wachter, S. 2011, AJ, 141, 115 [NASA ADS] [CrossRef] [Google Scholar]
  21. Cui, X.-Q., Zhao, Y.-H., Chu, Y.-Q., et al. 2012, Res. Astron. Astrophys., 12, 1197 [NASA ADS] [CrossRef] [Google Scholar]
  22. Cunha, K., Frinchaboy, P. M., Souto, D., et al. 2016, Astron. Nachr., in press [arXiv:1601.03099] [Google Scholar]
  23. Dalton, G., Trager, S. C., Abrams, D. C., et al. 2012, in Ground-based and Airborne Instrumentation for Astronomy IV, Proc. SPIE, 8446, 84460P [CrossRef] [Google Scholar]
  24. de Jong, R. S., Barden, S., Bellido-Tirado, O., et al. 2014, in Ground-based and Airborne Instrumentation for Astronomy V, Proc. SPIE, 9147, 91470M [Google Scholar]
  25. De Silva, G. M., Freeman, K. C., Bland-Hawthorn, J., et al. 2015, MNRAS, 449, 2604 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
  26. Donati, P., Beccari, G., Bragaglia, A., Cignoni, M., & Tosi, M. 2014a, MNRAS, 437, 1241 [NASA ADS] [CrossRef] [Google Scholar]
  27. Donati, P., Cantat Gaudin, T., Bragaglia, A., et al. 2014b, A&A, 561, A94 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  28. Donati, P., Bragaglia, A., Carretta, E., et al. 2015, MNRAS, 453, 4185 [NASA ADS] [CrossRef] [Google Scholar]
  29. Frasca, A., Biazzo, K., Lanzafame, A. C., et al. 2015, A&A, 575, A4 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  30. Friel, E. D., Janes, K. A., Tavarez, M., et al. 2002, AJ, 124, 2693 [NASA ADS] [CrossRef] [Google Scholar]
  31. Friel, E. D., Donati, P., Bragaglia, A., et al. 2014, A&A, 563, A117 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  32. Frinchaboy, P. M., Thompson, B., Jackson, K. M., et al. 2013, ApJ, 777, L1 [NASA ADS] [CrossRef] [Google Scholar]
  33. Froebrich, D., Scholz, A., & Raftery, C. L. 2007, MNRAS, 374, 399 [Google Scholar]
  34. Fujii, M. S., & Baba, J. 2012, MNRAS, 427, L16 [NASA ADS] [Google Scholar]
  35. Genovali, K., Lemasle, B., Bono, G., et al. 2013, A&A, 554, A132 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  36. Genovali, K., Lemasle, B., Bono, G., et al. 2014, A&A, 566, A37 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  37. Gilmore, G., Randich, S., Asplund, M., et al. 2012, The Messenger, 147, 25 [NASA ADS] [Google Scholar]
  38. Hayden, M. R., Holtzman, J. A., Bovy, J., et al. 2014, AJ, 147, 116 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
  39. Hayden, M. R., Bovy, J., Holtzman, J. A., et al. 2015, ApJ, 808, 132 [NASA ADS] [CrossRef] [Google Scholar]
  40. Haywood, M. 2006, MNRAS, 371, 1760 [NASA ADS] [CrossRef] [Google Scholar]
  41. 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]
  42. Heiter, U., Soubiran, C., Netopil, M., & Paunzen, E. 2014, A&A, 561, A93 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  43. Hill, V., Babusiaux, C., Gómez, A., et al. 2012, in EpJ Web Conf., 19, 6001 [Google Scholar]
  44. Jacobson, H. R., Pilachowski, C. A., & Friel, E. D. 2011, AJ, 142, 59 [NASA ADS] [CrossRef] [Google Scholar]
  45. Janes, K. A. 1979, ApJS, 39, 135 [NASA ADS] [CrossRef] [Google Scholar]
  46. Janes, K. A., & Hoq, S. 2011, AJ, 141, 92 [NASA ADS] [CrossRef] [Google Scholar]
  47. Jeffries, R. D., Totten, E. J., Harmer, S., & Deliyannis, C. P. 2002, MNRAS, 336, 1109 [NASA ADS] [CrossRef] [Google Scholar]
  48. Jílková, L., Carraro, G., Jungwiert, B., & Minchev, I. 2012, A&A, 541, A64 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  49. Kubryk, M., Prantzos, N., & Athanassoula, E. 2013, MNRAS, 436, 1479 [NASA ADS] [CrossRef] [Google Scholar]
  50. Loebman, S. R., Roškar, R., Debattista, V. P., et al. 2011, ApJ, 737, 8 [NASA ADS] [CrossRef] [Google Scholar]
  51. Loebman, S. R., Debattista, V. P., Nidever, D. L., et al. 2016, ApJ, 818, L6 [NASA ADS] [CrossRef] [Google Scholar]
  52. Magrini, L., Sestito, P., Randich, S., & Galli, D. 2009, A&A, 494, 95 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  53. Magrini, L., Randich, S., Zoccali, M., et al. 2010, A&A, 523, A11 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  54. Magrini, L., Randich, S., Romano, D., et al. 2014, A&A, 563, A44 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  55. Magrini, L., Randich, S., Donati, P., et al. 2015, A&A, 580, A85 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  56. Magrini, L., Coccato, L., Stanghellini, L., Casasola, V., & Galli, D. 2016, A&A, 588, A91 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  57. Majewski, S. R., Schiavon, R. P., Frinchaboy, P. M., et al. 2015, AJ, submitted [arXiv:1509.05420] [Google Scholar]
  58. Martin, R. P., Andrievsky, S. M., Kovtyukh, V. V., et al. 2015, MNRAS, 449, 4071 [NASA ADS] [CrossRef] [Google Scholar]
  59. Martínez-Barbosa, C. A., Brown, A. G. A., & Portegies Zwart, S. 2015, MNRAS, 446, 823 [NASA ADS] [CrossRef] [Google Scholar]
  60. Martinez-Medina, L. A., Pichardo, B., Moreno, E., Peimbert, A., & Velazquez, H. 2016, ApJ, 817, L3 [NASA ADS] [CrossRef] [Google Scholar]
  61. Mikolaitis, Š., Hill, V., Recio-Blanco, A., et al. 2014, A&A, 572, A33 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  62. Minchev, I., Chiappini, C., & Martig, M. 2013, A&A, 558, A9 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  63. Minchev, I., Chiappini, C., & Martig, M. 2014, A&A, 572, A92 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  64. Netopil, M., Paunzen, E., Heiter, U., & Soubiran, C. 2016, A&A, 585, A150 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  65. Pancino, E., Carrera, R., Rossetti, E., & Gallart, C. 2010, A&A, 511, A56 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  66. Piatti, A. E., Clariá, J. J., Bica, E., Geisler, D., & Minniti, D. 1998, AJ, 116, 801 [NASA ADS] [CrossRef] [Google Scholar]
  67. Prisinzano, L., Damiani, F., Micela, G., et al. 2016, A&A, 589, A70 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  68. Randich, S., Gilmore, G., & Gaia-ESO Consortium. 2013, The Messenger, 154, 47 [NASA ADS] [Google Scholar]
  69. Recio-Blanco, A., de Laverny, P., Kordopatis, G., et al. 2014, A&A, 567, A5 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  70. Roškar, R., Debattista, V. P., Quinn, T. R., Stinson, G. S., & Wadsley, J. 2008, ApJ, 684, L79 [NASA ADS] [CrossRef] [Google Scholar]
  71. Sacco, G. G., Morbidelli, L., Franciosini, E., et al. 2014, A&A, 565, A113 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  72. Schönrich, R., & Binney, J. 2009, MNRAS, 396, 203 [NASA ADS] [CrossRef] [Google Scholar]
  73. Smiljanic, R., Korn, A. J., Bergemann, M., et al. 2014, A&A, 570, A122 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  74. Smiljanic, R., Romano, D., Bragaglia, A., et al. 2016, A&A, 589, A115 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  75. Stanghellini, L., & Haywood, M. 2010, ApJ, 714, 1096 [NASA ADS] [CrossRef] [Google Scholar]
  76. Steinmetz, M., Zwitter, T., Siebert, A., et al. 2006, AJ, 132, 1645 [NASA ADS] [CrossRef] [Google Scholar]
  77. Sung, H., Bessell, M. S., Lee, B.-W., & Lee, S.-G. 2002, AJ, 123, 290 [NASA ADS] [CrossRef] [Google Scholar]
  78. Tosi, M. 1982, ApJ, 254, 699 [NASA ADS] [CrossRef] [Google Scholar]
  79. Twarog, B. A., Ashman, K. M., & Anthony-Twarog, B. J. 1997, AJ, 114, 2556 [NASA ADS] [CrossRef] [Google Scholar]
  80. Yanny, B., Rockosi, C., Newberg, H. J., et al. 2009, AJ, 137, 4377 [NASA ADS] [CrossRef] [Google Scholar]
  81. Yong, D., Carney, B. W., & Teixera de Almeida, M. L. 2005, AJ, 130, 597 [NASA ADS] [CrossRef] [Google Scholar]
  82. Yong, D., Carney, B. W., & Friel, E. D. 2012, AJ, 144, 95 [NASA ADS] [CrossRef] [Google Scholar]
  83. York, D. G., Adelman, J., Anderson, Jr., J. E., et al. 2000, AJ, 120, 1579 [Google Scholar]
  84. Zhang, B., Chen, X.-Y., Liu, C., et al. 2015, Res. Astron. Astrophys., 15, 1197 [NASA ADS] [CrossRef] [Google Scholar]
  85. Zhao, G., Chen, Y.-Q., Shi, J.-R., et al. 2006, Chin. J. Astron. Astrophys., 6, 265 [Google Scholar]

All Tables

Table 1

Cluster parameters.

Table 2

Parameters for cluster members.

All Figures

thumbnail Fig. 1

Top panel: locations of the GES open clusters (red triangles) relative to the Galactic plane (z = 0, solid line), as a function of their Galactocentric radii. All clusters lie within 150 pc of the mid-plane, and thus are members of the thin disk. Bottom panel: cluster ages (in Gyr) plotted as a function of RGC.

In the text
thumbnail Fig. 2

Metallicity gradient as shown by the GES open clusters (red triangles) in comparison to literature studies. The open cluster literature compilation of Netopil et al. (2016) is shown as open gray squares, while clusters observed by the APOGEE survey (Frinchaboy et al. 2013) are filled gray squares. The solid blue line shows the metallicity gradient as determined by APOGEE field giant stars within 250 pc of the Galactic mid-plane (Hayden et al. 2014). The red dashed line is least-squares fit to the GES sample, with the value of the slope given in the caption. [Fe/H] = 0 is indicated by the solid black line, to guide the eye.

In the text
thumbnail Fig. 3

Same as in Fig. 2, but with distances determined by Netopil et al. (2016) for GES clusters in their compilation. Note the magnitude of the GES cluster gradient has changed by 0.01 dex kpc-1.

In the text
thumbnail Fig. 4

GES open clusters (black-lined triangles) and Netopil et al. (2016) clusters (filled squares), color-coded by age. The colored solid lines represent the metallicity gradients for different-aged stars in the chemical evolution model of Minchev et al. (2014). The gradient shown by the clusters is much steeper than that of the model.

In the text
thumbnail Fig. 5

Residuals in [Fe/H] from the linear gradient as a function of cluster age.

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.