Issue |
A&A
Volume 580, August 2015
|
|
---|---|---|
Article Number | A85 | |
Number of page(s) | 15 | |
Section | Galactic structure, stellar clusters and populations | |
DOI | https://doi.org/10.1051/0004-6361/201526305 | |
Published online | 06 August 2015 |
The Gaia-ESO Survey: Insights into the inner-disc evolution from open clusters⋆
1 INAF–Osservatorio Astrofisico di Arcetri, Largo E. Fermi, 5, 50125 Firenze, Italy
e-mail: laura@arcetri.astro.it
2 INAF–Osservatorio Astronomico di Bologna, via Ranzani 1, 40127 Bologna, Italy
3 Dipartimento di Fisica e Astronomia, Universitá di Bologna, via Ranzani 1, 40127 Bologna, Italy
4 Instituto de Astrofísica e Ciências do Espaço, Universidade do Porto, CAUP, Rua das Estrelas, 4150-762 Porto, Portugal
5 Department for Astrophysics, Nicolaus Copernicus Astronomical Center, ul. Rabiańska 8, 87-100 Toruń, Poland
6 Observatoire de Genève, Université de Genève, 1290 Versoix, Switzerland
7 Institute of Theoretical Physics and Astronomy, Vilnius University, Gostauto 12, 01108 Vilnius, Lithuania
8 Astronomy Department, Indiana University, 727 East 3rd Street, Bloomington, IN 47405, USA
9 Department of Physics and Kavli Institute for Astrophysics and Space Research, Massachusetts Institute of Technology, 77 Massachusetts Avenue, Cambridge, MA 02139, USA
10 Dipartimento di Fisica e Astronomia, Universitá di Padova, vicolo Osservatorio 3, 35122 Padova, Italy
11 INAF–Osservatorio Astronomico di Padova, vicolo Osservatorio 5, 35122 Padova, Italy
12 ASI Science Data Center, via del Politecnico SNC, 00133 Roma, Italy
13 Departamento de Astronomía, Casilla 160-C, Universidad de Concepción, Concepción, Chile
14 Centro de Estudios de Física del Cosmos de Aragón (CEFCA), Plaza San Juan 1, 44001 Teruel, Spain
15 Institute of Astronomy, University of Cambridge, Madingley Road, Cambridge CB3 0HA, UK
16 GEPI, Observatoire de Paris, CNRS, Université Paris Diderot, 5 place Jules Janssen, 92190 Meudon, France
17 Dept. of Astronomy and Theoretical physics, Lund university, Box 43, 22100 Lund, Sweden
18 INAF–Osservatorio Astronomico di Palermo, Piazza del Parlamento 1, 90134, Palermo, Italy
19 Department of Physics and Astronomy, Uppsala University, Box 516, 75120 Uppsala, Sweden
20 Laboratoire Lagrange (UMR 7293), Université de Nice Sophia Antipolis, CNRS, Observatoire de la Côte d’Azur, CS 34229, 06304 Nice Cedex 4, France
21 European Southern Observatory, Alonso de Cordova 3107 Vitacura, Santiago de Chile, Chile
22 Instituto de Astrofísica de Andalucía-CSIC, Apdo. 3004, 18080 Granada, Spain
23 Astrophysics Research Institute, Liverpool John Moores University, 146 Brownlow Hill, Liverpool L3 5RF, UK
24 Departamento de Ciencias Físicas, Universidad Andrés Bello, República 220, 837-0134 Santiago, Chile
25 Departamento de Física e Astronomia, Faculdade de Ciências, Universidade do Porto, Rua do Campo Alegre, 4169-007 Porto, Portugal
Received: 13 April 2015
Accepted: 13 May 2015
Context. The inner disc, which links the thin disc with the bulge, has been somewhat neglected in the past because of the intrinsic difficulties in its study, among which crowding and high extinction. Open clusters located in the inner disc are among thebest tracers of its chemistry at different ages and distances.
Aims. We analyse the chemical patterns of four open clusters located within 7 kpc of the Galactic centre and of field stars to infer the properties of the inner disc with the Gaia-ESO survey idr2/3 data release.
Methods. We derive the parameters of the newly observed cluster, Berkeley 81, finding an age of about 1 Gyr and a Galactocentric distance of ~5.4 kpc. We construct the chemical patterns of clusters and we compare them with those of field stars in the solar neighbourhood and in the inner-disc samples.
Results. Comparing the three populations we observe that inner-disc clusters and field stars are both, on average, enhanced in [O/Fe], [Mg/Fe], and [Si/Fe]. Using the idr2/3 results of M67, we estimate the non-local thermodynamic equilibrium (NLTE) effect on the abundances of Mg and Si in giant stars. After empirically correcting for NLTE effects, we note that NGC 6705 and Be 81 still have a high [α/Fe].
Conclusions. The location of the four open clusters and of the field population reveals that the evolution of the metallicity [Fe/H] and of [α/Fe] can be explained within the framework of a simple chemical evolution model: both [Fe/H] and [α/Fe] of Trumpler 20 and of NGC 4815 are in agreement with expectations from a simple chemical evolution model. On the other hand, NGC 6705, and to a lesser degree Berkeley 81, have higher [α/Fe] than expected for their ages, location in the disc, and metallicity. These differences might originate from local enrichment processes as explained in the inhomogeneous evolution framework.
Key words: Galaxy: abundances / open clusters and associations: general / open clusters and associations: individual: Trumpler 20 / Galaxy: disk / open clusters and associations: individual: NGC 6705 / globular clusters: individual: NGC 4815
© ESO, 2015
1. Introduction
Open clusters are among the best tracers of the chemical evolution of the Galactic thin disc (see e.g. Friel 1995; Bragaglia & Tosi 2006) from the very distant outskirts (e.g. Sestito et al. 2008; Yong et al. 2012) to the inner regions close to the Galactic bulge (see e.g. Magrini et al. 2010, 2014). In particular, the clusters in the so-called inner disc, i.e. the part of the Galactic thin disc located at Galactocentric distances ≤8 kpc1, have been explored relatively little: indeed, few clusters are known and have been spectroscopically observed (e.g. Carretta et al. 2005, 2007; Sestito et al. 2007; Magrini et al. 2010, 2014; Carraro et al. 2014a,b,c) both because of observational limits, i.e. high reddening and crowding towards the Galactic centre, and because of the high mortality of clusters in regions where the density of stars is higher (Freeman 1970; van der Kruit 2002). However, this region is of great importance because it constitutes a connection between the properties of the bulge and of the thin/thick disc.
The Gaia-ESO Survey (Gilmore et al. 2012; Randich et al. 2013) is providing high-resolution spectroscopy with FLAMES at VLT and homogeneous stellar parameters and abundances for a large sample of Galactic open clusters. The sample includes clusters located throughout the disc, at different distances from the Sun, and with ages ranging from their initial phases (with the exception of the embedded ones), i.e. a few Myr after their formation, up to several Gyr (Randich et al., in prep.). The selection of targets to be observed within clusters is as uniform as possible: Gaia-ESO observes unbiased samples of cluster candidates with GIRAFFE, i.e. their selection is based only on the location of their colour–magnitude diagrams (CMDs) and no other information, such as radial velocity, is used to ensure their membership; more secure cluster members are observed with UVES (see Bragaglia et al., in preparation for more details).
Observations for several clusters were completed in the first 18 months of the Gaia-ESO Survey. Analysis is now completed and the products have been released to the consortium as internal data release 2, hereafter idr2/3. The sample includes four very interesting intermediate-age/old clusters in the inner disc: NGC 6705, Trumpler 20, NGC 4815, and Berkeley 81. The first three have been discussed in a series of previous paper using Gaia-ESO idr1 recommended values (e.g. Magrini et al. 2014; Tautvaišienė et al. 2015). In the present paper, we consider the updated idr2/3 analysis, with important differences in the determination of the final parameters with respect to idr1 (see Smiljanic et al. 2014, for details); we also add to the sample Berkeley 81 (hereafter Be 81), which is the innermost cluster in the idr2/3 sample. We aim at comparing its abundance pattern with both field stars located in the solar neighbourhood and toward the Galactic centre, and with the other three clusters located in the inner disc.
Be 81 is an intermediate-age cluster located at approximately 6 kpc from the Galactic centre (l = 34.505, b = −2.068). It has been recently studied by Donati et al. (2014a), who, comparing observed and synthetic CMDs, found an age close to 1 Gyr, a large reddening E(B−V) = 0.9, a distance modulus of 12.4, and a metallicity close to solar. Previous literature data on Be 81 are quite scarce: BVI photometry of a portion of Be 81 was obtained by Sagar & Griffiths (1998), from which they derived an age of ~1 Gyr and a reddening E(B−V) = 1. The only determination of the metallicity of this cluster is by Warren & Cole (2009) through calcium triplet spectroscopy. They obtained a rather uncertain sub-solar metallicity ([Fe/H] = −0.11 ± 0.15 dex), somewhat unusual for its internal location in the Galactic disc. This is likely related to the limiting magnitude of Warren & Cole (2009), which did not allow them to reach stars as faint as the cluster giant branch and so they likely were not able to identify any of the cluster member stars.
The paper is structured as follows. In Sect. 2 we describe the analysis of Be 81 in idr2/3. In Sect. 3 we derive its fundamental parameters, whereas in Sect. 4 we compare the chemical patterns of field and cluster populations. In Sect. 5 we discuss the chemical patterns, in particular the α-enhancement, in the context of our current knowledge of the disc’s evolution. In Sect. 6 we give our summary and conclusions.
Fig. 1 Upper left panel: in grey the Be 81 field of view of the photometric data in Donati et al. (2014a), obtained with LBC at LBT; in black the stars observed by the Gaia-ESO survey, all inside a radius of 6 arcmin. Lower left panel: distribution of the radial velocities vs. distance from the cluster centre. The red line is the cluster systemic radial velocity. Right panel: CMD of the inner 6 arcmin, with stars observed indicated by larger symbols (orange squares HR09B; red crosses HR15N; blue circles UVES). |
2. Sample and analysis of Be 81
In Be 81, 14 stars were observed with UVES using the setup U580 (R = 47 000, λλ 4800–5800 Å); 266 stars were instead targeted with GIRAFFE using HR15N (R = 17 000, λλ 6470–6790 Å) and HR09B (R = 25 900, λλ 5143–5356 Å) (16 stars observed with both setups, 80 with HR09B only, and 170 HR15N with only). In Fig. 1 (right panel) we show the CMD of the inner 6 arcmin of the cluster. The UVES targets are candidate red clump giant stars, while the candidate main sequence stars were observed with GIRAFFE. In the upper left panel, the field of view of the photometric data in Donati et al. (2014a) obtained with the Large Binocular Camera at the Large Binocular Telescope (LBC at LBT) is presented. Finally in the lower left panel, we show the distribution of the radial velocities vs. distance from the cluster centre.
Parameters and their uncertanities for the UVES targets in Be 81.
Parameters and their uncertainties for the GIRAFFE RV members Be 81 from WG 10.
The spectrum analysis of the Gaia-ESO data is organised in several working groups (WGs) that provide astrophysical parameters and element abundances. In particular, the analysis of the UVES and GIRAFFE spectra discussed in the present paper was performed by three WGs: one dedicated to FGK-type stars observed with UVES (WG11; see Smiljanic et al. 2014 for details), one for the same type of stars targeted with GIRAFFE (WG10; Recio-Blanco et al., in prep.), and another one devoted to warmer stars observed with both UVES and GIRAFFE (WG13; Blomme et al., in prep.).
The atmospheric parameters obtained by the different WGs were then combined by a WG dedicated to the homogenisation (WG 15) that provides a final set of recommended parameters (Hourihane et al., in prep.). More specifically, the homogenisation process compares the results of the different WG for common targets and calibrators and, based on this comparison, derives and applies offsets to obtain a common scale for all the Gaia-ESO Survey results. We mention that the idr2/3 WG 15 homogenisation process has not affected the recommended values of effective temperature (Teff), surface gravity (log g), and microturbolent velocity (ξ), while it has derived and applied offsets in metallicity [Fe/H] and in radial velocity (RV). Furthermore, the homogenisation flow of idr2/3 has not involved the homogenisation of the elemental abundances through different WG. Therefore, the elemental abundances used in the present work are those provided as recommended values by WG11 (Smiljanic et al. 2014).
All the targets observed in Be 81 field have recommended atmospheric parameters and radial velocities in the final table of WG 15. The UVES targets also have element abundance determinations from WG11. In Table 1 we report the information about the UVES stars: their identification names (Cname)2 and id3, equatorial coordinates, radial velocities, stellar parameters with associated errors, RV from Hayes & Friel (2014) with associated errors, and membership flag (described in Sect. 2.2). The idr2/3 RVs are in good agreement with those of Hayes & Friel (2014). In Table 2 we show the information about confirmed members (based on their RV) observed with GIRAFFE and analysed by WG10 (see Sect. 2.2 for the membership probability definition). Their stellar parameters are obtained with spectral indices (Damiani et al. 2014) and thus the microturbulent velocity ξ is not determined. Stellar parameters of WG 13 are not used here because of their large uncertainties, and therefore they are not presented in Table 2.
2.1. The field comparison sample
The solar neighbourhood.
For the MW field stars used as comparison in the present work, we considered the analysis of WG11. Specifically, the UVES solar neighbourhood sample aims at obtaining an unbiased sample of ~5000 G-stars within ~2 kpc of the Sun in the turn-off phase (see Gilmore et al. in prep. for a description of the samples). The final purpose is to quantify the local elemental abundance distribution functions in detail. For this paper we used the whole sample of 861 Milky Way solar neighbourhood stars with recommended parameters.
The bulge and inner disc.
UVES stars observed in parallel with the GIRAFFE targets dedicated to the study of the Galactic Bulge are expected to be evolved stars belonging to both bulge and inner-disc populations. During the first 18 months of the Gaia-ESO Survey, 80 stars of the bulge/inner-disc were observed and are included in the present discussion.
2.2. Membership
The determination of the RVs for GIRAFFE and UVES spectra, respectively, can be found in Koposov et al. (in prep.) and (Sacco et al. 2014). We have used the recommended radial velocity of the UVES and GIRAFFE sample to define the cluster systemic velocity.
The resulting RV distribution, including both GIRAFFE and UVES targets, is shown in Fig. 2. A Gaussian fit of the peak of the distribution identifying the cluster signature with respect to the field stars gives a central velocity of 45.8 km s-1 and a full width half maximum (FWHM) of 3.1 km s-1. To avoid a large contamination by field stars, we define a star as RV cluster member if it has a RV within 1σ of the average cluster radial velocity. By doing so we identify 56 RV members, 14 observed with UVES and 42 targeted with GIRAFFE (5 observed with both HR09B and HR15N, 11 with HR09B, 9 with HR15N and analysed by WG13, and 18 others observed with HR15N and analysed by WG10). In summary, 20% of the observed targets (almost 100% for UVES and 16% for GIRAFFE) have a RV consistent with membership. For UVES, the contamination by field stars has been reduced by observing high probability member stars whose radial velocity were obtained by Hayes & Friel (2014) to assist Gaia-ESO observations. For GIRAFFE, the low fraction of recovered cluster members is due to the location of the cluster towards the Galactic centre where a large contamination by field stars is expected.
Fig. 2 Distribution of radial velocities of the stars observed in the field of Be 81. The solid line is a Gaussian fit of the distribution of the member stars. |
In Fig. 3 we show [Fe/H] vs. atmospheric parameters of the UVES stars. We note that stars observed with UVES are selected to be red giant stars in the clump and so they are expected to have similar atmospheric parameters. The average error (as given in the recommended value table) is plotted on the left side of each panel. The 1σ and 2σ areas around the median value4 are also marked. We note that the star with the lowest metallicity (19015978-0028183) is beyond 1σ from the median Teff and 2σ from the median log g, likely being a non-member with radial velocity consistent with the cluster systemic velocity. This is not surprising given the expected high contamination of field stars. The most metal-rich star (19013651-0027021) is also somewhat deviant and might be a non-member star as well. In the next sections we discuss the chemical homogeneity of cluster member stars and how it can help to disentangle member vs. non-member stars.
In Fig. 4 we plot the histogram of the Gaia-ESO recommended metallicities, [Fe/H], of the RV member stars observed with UVES and GIRAFFE. The two distributions are compatible: the median of the UVES stars is [Fe/H] = 0.23 ± 0.08 dex, whereas for the GIRAFFE stars is 0.16 ± 0.44 dex. The UVES metallicities are much more concentrated in the bin centred at 0.25 dex, while the GIRAFFE metallicities are more dispersed owing to the lower resolution and smaller spectral range. We note that the errors in Table 2 are based on a single determination of stellar parameters, while the errors in Table 1 are computed with the MAD of typically ten determinations. Thus the errors on stellar parameters and metallicities of the GIRAFFE stars are likely underestimated.
In Fig. 5 we present the HR diagram of the candidate member stars based on RV with a PARSEC isochrone (Bressan et al. 2012) suitable for the age and metallicity of Be 81 (see next section). We show the UVES and GIRAFFE stars with different colours. The most metal-poor UVES star is located in the upper part of the plot, quite far away from the location of the other clump stars. The most metal-rich star is located just below the main locus of the clump stars. The GIRAFFE targets were selected among main sequence (MS) stars. Several members follow, within the errors, the PARSEC isochrone, while other stars have very uncertain parameters and are not in agreement with the isochrone.
Fig. 3 [Fe/H] vs. atmospheric parameters of UVES stars in Be 81. The dotted lines mark the 1σ region, while the dashed lines indicate the 2σ area. The typical error bar on each star are shown on the left side of each panel. |
3. Fundamental parameters of Be 81
3.1. Chemical properties of Be 81
The high-resolution and high-quality UVES spectra allow WG11 to determine the detailed composition of the stars of Be 81. In Table 3 we report the elemental abundances of the UVES stars. These are graphically presented in Fig. 6, where we show [El/Fe] of each individual UVES star. Their median values, computed excluding the metal-poor non-member star are shown in Table 4. The abundances are scaled consistently to the solar abundances obtained in the analysis of WG11. The adopted iron abundances are those listed in the [Fe/H] column of idr2/3 table of recommended parameters5. In Table 5 we show the atmospheric parameters and elemental abundances obtained for the selected solar spectrum. For elements without a WG11 determination, we scaled with the solar abundance of Grevesse et al. (2007).
Interestingly, the abundance distribution confirms the compactness of the cluster in most abundance ratios, whereas the most metal-poor star has a different chemical pattern. This is demonstrated in Fig. 6, which shows the dispersions around the median values (dashed lines) computed with the Robust sigma IDL routine and indicated with σ in Table 4, as well as the typical error bar (black crosses in the left part of each panel and indicated with ⟨ δ ⟩) computed as the average of the individual stellar errors. The error on the parameters and abundances of each star are discussed in Smiljanic et al. (2014), where they have been estimated with the method-to-method dispersion. For most elements these two quantities are comparable indicating a high degree of chemical homogeneity of the cluster, as already noticed for NGC 6705, Tr 20, and NGC 4815 with idr1 data by Magrini et al. (2014). Be 81 is slightly enhanced in some of the α-elements and in one of the neutron capture elements, having [O/Fe], [Mg/Fe], [Si/Fe], and [Eu/Fe] > 0.1 dex, while it has [Ti/Fe] and [Ca/Fe] close to zero. In this context, the most metal-poor star is clearly distinguishable not only because of its lower [Fe/H], but also for differing more than 1σ in several elemental abundances, such as [O/Fe], [Na/Fe], [Mg/Fe], [Al/Fe], [Si/Fe], [Ti/Fe], [Cr/Fe], and [Eu/Fe]. Taken all together these chemical signatures clearly indicate that 19015978-0028183 is not a member star, while the most metal-rich star (19013651-0027021) is consistent within 1σ in all panels with the median values of the cluster and thus can be considered a bona fide member star.
Fig. 4 Distribution of [Fe/H] of the stars observed in the field of Be 81. The cyan rectangles show the GIRAFFE metallicities of RV member stars, while the blue ones indicate the UVES metallicities (including star 19015978-0028183). |
Fig. 5 HR diagram with the PARSEC isochrone for an age of 0.98 Gyr and a metallicity Z = 0.025 (red curve) and the recommended Gaia-ESO parameters of UVES stars (green squares) and of GIRAFFE stars (blue circles). The stars that are members on the basis of their radial velocity but with [Fe/H] outside 1σ from the cluster average metallicity are shown with empty symbols. |
3.2. Age, distance, and reddening of Be 81
Thanks to accurate metallicity measurements and membership based on RV determinations, we are able to derive improved age, distance, and average reddening values for Be 81 by using the classical approach of isochrone fitting. As was done for the open cluster Trumpler 20 (see Donati et al. 2014b), we chose four different sets of isochrones to take into account how the results change when adopting different models, and we quantified the effect of different input physics on the derived parameters: PARSEC (Bressan et al. 2012), BASTI (Pietrinferni et al. 2004), Dartmouth (Dotter et al. 2008), and Victoria-Regina (VandenBerg et al. 2006). We adopted the B, V, and I photometric catalogue described in Donati et al. (2014a). The accurate spectroscopic analysis of the cluster metallicity from the UVES spectra clearly indicates a super-solar chemical mixture for Be 81, in contrast to previous findings based only on Ca-triplet spectra. Using this information, we are able to select a more appropriate set of isochrones and thus to obtain a more accurate determination of the cluster evolutionary status, reddening, and distance modulus.
The best-fitting isochrone is chosen by visual examination as the one that can describe at the same time the main age-sensitive evolutionary phases: the luminosity and colour of the main sequence turn-off (MSTO), red hook (RH), and red clump (RC), defined in Donati et al. (2014a), in two different photometric colours (B−V and V−I, adopting RV = 3.1, AB/AV = 1.29719, and AI/AV = 0.60329 taken from Cardelli et al. 1989 and O’Donnell 1994). We started by using the metallicity resulting from the UVES spectroscopic estimate, i.e. [Fe/H] = + 0.23 dex (see Sect. 2.2), and we converted it to Z taking into account the different solar abundances of the four sets of isochrones. The difference in luminosity between the MSTO and RC is adopted as the most reliable constraint on the cluster age. In fact, while the colour difference between the same two indicators is also related to age, it is used here as a secondary age indicator because colour properties are more affected by theoretical uncertainties, like colour transformations and the super-adiabatic convection, while luminosity constraints are more reliable. The errors on the estimated cluster parameters are mainly due to the uncertainties in the definition of the age indicators because of the heavy contamination by field interlopers, while the photometric error does not have a significant impact on the error budget (see Donati et al. 2014a). We make use of the RV to select the most probable cluster members among the many stars observed in order to highlight the most probable cluster sequences in the CMD and to have more stringent constraints on the best-fitting models. Moreover, we tweak the age determination looking for the best-fitting isochrone in the theoretical plane Teff, log g by using the RC star parameters of the UVES targets (see Table 1).
Abundances in the 12+log (X/H) form of the UVES stars in Be 81 from idr2/3 recommended values.
Median elemental abundances with their dispersions and average errors.
idr2/3 solar parameters and abundances.
Fig. 6 Abundance ratios of stars in Be 81. The dotted lines indicate the median+ σ and median−σ regions in both [El/Fe] and [Fe/H]. The typical error bar on each star are shown on the left side of each panel. |
For the PARSEC set (see Fig. 7, left panel), for which Z⊙ = 0.015, we adopted Z = 0.025. The luminosity and colour of the RC are nicley reproduced, as are those of the RH and upper MS. In Fig. 5 we show, as an example, the isochrone fit in the theoretical plane which, together with the photometric data, was used to determine the best age, i.e. 1.0 ± 0.1 Gyr. We note that in the case of the V,V−I CMD the RC and MS colours of the model are slightly bluer than in the observations, and this is also evident in the best fit with the three other isochrones. As mentioned above, this effect may be due to uncertainties in the temperature-colour transformations or to the differences in chemical mixtures between models and observations or, finally, to small systematics in the B and I photometry that in principle might have a larger impact on colours, of the order of some hundredths of magnitude. Considering, however, that metallicity is fixed by spectroscopy and age by the difference in luminosity between RC and TO, we actually used V,B−V and V,V−I photometry to adjust the distance modulus and reddening, and deemed this fit quite satisfying and the small mismatch of the RC colour acceptable. In the case of the BASTI isochrones (Z = 0.03), the best fit is obtained for a younger age of 0.75 ± 0.1 Gyr. The Dartmouth isochrones and the Victoria-Regina isochrones do not include the evolved phases after the RGB for the age of Be 81 and thus the cluster’s parameters for these two models are more loosely constrained. We proceeded by choosing the best-fitting solution as the one that best matches the MS and MSTO morphologies (see Fig. 7, right panel). We used models with [Fe/H] = +0.23 for both sets of isochrones. The results obtained with the four sets are given in Table 6, where we indicate age, distance modulus, reddening, distance from the Sun and the Galactocentric distance from the Galactic plane, and mass of the stars at the MSTO. The cluster’s centre coordinates are taken from Donati et al. (2014a): RA = 19:01:42.8 and Dec = −0:27:07.7.
Since from spectroscopic analysis it turned out that Be 81 stars show an α-enhanced chemical pattern ([α/Fe] ~ +0.2, see Sect. 3.1), we also looked at the model predictions when an α-enhanced mixture is taken into account. Slightly younger ages are obtained. The results are reported in Table 6.
Fig. 7 CMD obtained for stars inside 3′ using the photometry from Donati et al. (2014a) and best-fit isochrone for different evolutionary models (from left to right: PARSEC, BASTI, DARTMOUTH, and Victoria-Regina). GES target non-members are plotted as black crosses, while members are highlighted with orange (GIRAFFE) and blue (UVES) points. See Table 6 for the adopted parameters for the isochrone fitting. |
Results, errors, and estimated systematic uncertainties using different evolutionary models with [Fe/H] ~ 0.2.
Cluster parameters.
With respect to literature values (see Donati et al. 2014a; Sagar & Griffiths 1998), we found an acceptable agreement with the determination of the age, distance, and reddening of the cluster.
4. A comparison with inner disc clusters and field stars
The large variety of available chemical elements belonging to different nucleosynthesis channels allows us to compare, in a fully consistent way, the detailed chemical patterns of open clusters and of field stars observed and analysed homogeneously by the Gaia-ESO Survey. In addition to Be 81, three other old/intermediate-age clusters are available in idr2/3, namely Tr 20, NGC 6705, and NGC 4815. Their stellar parameters, abundances, and radial velocities were also provided in the first data release, idr1, and analysed in several papers: Donati et al. (2014b) derived membership, age, distance, and reddening of Tr 20; Friel et al. (2014) derived the parameters of NGC 4815 and studied its chemical composition; Cantat-Gaudin et al. (2014) studied the properties of NGC 6705 and searched for possible inhomogeneities in its chemical composition; finally Magrini et al. (2014) searched for unique signatures in the chemical patterns of the three clusters as possible indications of their birthplace. Their CNO abundances have been discussed by Tautvaišienė et al. (2015) using idr2/3 results. Their main parameters are summarised in Table 7.
Here we use the idr2/3 stellar parameters and abundances of NGC 4815, NGC 6705, and Tr 20 to compare their chemical patterns with Be 81 and with field stars. Although shown in Fig. 6, in the following we do not discuss Al and Na in terms of Galactic chemical evolution. They are indeed potentially interesting elements for comparison with the field samples and deduction of the chemical evolution of the inner disc. However, the abundances derived in idr2/3 are all computed in local thermodynamic equilibrium (LTE), and those elements are affected by relevant non-LTE (NLTE) effects. All stars observed with UVES in Be 81 are in the same evolutionary phase (RC), which means that the internal dispersion in abundance is not affected by this problem. However, the absolute value of Na and Al abundances and the comparison with stars of different parameters could be unreliable without these corrections. Furthermore, both Na and Al are subject to extra-mixing processes in stars of the mass involved here (see e.g. Charbonnel & Lagarde 2010; Karakas & Lattanzio 2014). These elements will be discussed in a forthcoming paper (Smiljanic et al., in prep.) and are not be used here.
As is well known, the ratio of α-elements to iron, [α/Fe], is commonly used to trace the star-formation timescale in a system because it is sensitive to the ratio of SNII (massive stars) to SNIa (intermediate mass binary systems with mass transfer) that have occurred in the lifetime of the system. However, the so-called α-elements, which are generally associated with the explosion of massive stars, do not have exactly the same place of production, and some of them (e.g. Si, Ca, and Ti) also have a significant contribution from SNIa. In addition, there are also differences in the way these elements are produced in massive stars: O and Mg are produced during the hydrostatic He burning, and their yields are not expected to be significantly affected by the SNII explosion conditions, while Si, Ca, and Ti are mostly produced during the SNII explosion. This distinction is also seen in the observations (e.g. Fulbright et al. 2007), where Si, Ca, and Ti usually track one another, but O and Mg often show different trends with [Fe/H]. The “pure” α-elements (O and Mg) are produced almost entirely by massive stars, and thus in shorter timescales. In the bulge, the observed Mg overabundances and the general α-enhancement agree with the predictions of the chemical evolution models of Matteucci & Brocato (1990) and Grieco et al. (2012) and they are explained with a rapid formation timescale for the bulge. A summary of the place of production and of the contribution of SNII and SNIa is shown in Table 8. The table is computed using the information given in Woosley et al. (2002) and in Iwamoto et al. (1999, model W7) and summarises the origin of the main isotope of each element.
Fig. 9 Abundance ratios of the idr2/3 old and intermediate-age open clusters located in the inner disc as a function of the atomic number. In black are the abundance ratios of Be 81, in yellow NGC 6705, in green NGC 4815, and in red Tr 20. The dashed line marks the solar abundance, and the shaded area the 1σ regions with σ = 0.1 dex. |
Stellar nucleosynthesis.
Fig. 10 Abundance ratios vs. atomic numbers of the idr2/3 stars in the solar neighbourhood and in the inner disc/bulge (selected in the [Fe/H] range from −0.10 to +0.33 dex, i.e. including the metallicity ranges of the four clusters) compared with the chemical pattern of NGC 6705, Trumpler 20, NGC 4815, and Be 81. The continuous lines mark the solar abundances, and the shaded area the 1σ region with σ = 0.1 dex. |
In Fig. 9 we show the median [El/Fe] abundance ratios vs. the atomic number in the four old/intermediate-age clusters located within 7 kpc of the Galactic centre, whereas in Fig. 10 we show the chemical pattern of the four clusters compared to the median abundance ratios of stars with −0.10 < [Fe/H] <+ 0.33 dex in the solar neighbourhood [382 stars] and in the inner-disc/bulge [42 stars] samples, i.e. embracing the maximum metallicity range of the four clusters within 1σ. We have excluded all stars with errors on [Fe/H] larger than 0.15 dex.
In both figures, the available elements include the α-elements (O, Mg, Si, Ca, Ti), the refractory elements (Sc, V), the iron-peak elements (Cr, Ni), and the neutron capture elements (Y, Eu). From Fig. 9, we can make some general considerations on the similarities to the chemical patterns of the inner clusters.
Considering the errors that affect the abundance ratios, we define an element as enhanced with respect to solar if [El/Fe] ≥ 0.1 dex. In addition, we consider that the abundance of an element is globally enhanced in the inner cluster sample if at least three of the four clusters have [El/Fe] ≥ 0.1 dex for that element.
NGC 6705, NGC 4815, and Be 81 have enhanced abundances of the α-elements O, Mg, and Si, with NGC 6705 having the highest abundance ratios of [Mg/Fe] and [Si/Fe], while [O/Fe] does not exceed a 0.1 dex enhancement. The group of elements from Ca to Ni shows solar abundance and they have similar values in the four clusters. The elements Y and Eu are very similar in all clusters, close to solar value for [Y/Fe] and slightly super-solar for [Eu/Fe]. From Fig. 10, we find that the chemical pattern of Be 81 is nearly identical to that of inner-disc/bulge stars: both populations show similar α-enhancements for the lower atomic number α elements (O, Mg, and Si), while Ca, Ti, and iron-peak elements are almost solar. For comparison, the solar neighbourhood stars in the same metallicity range are plotted: they all have [Mg/Fe], [Si/Fe], and [O/Fe] close to zero.
Fig. 11 Abundance ratios [Si/H] and [Mg/H] vs. Teff of member stars in M 67 from idr2/3 analysis. The continuous lines mark the ±1σ from the average value, whereas the dotted lines indicate the ±2σ region. |
5. Inner clusters and disc evolution
For the first time, thanks to the Gaia-ESO Survey, we have at our disposal a significant number of open clusters together with field stars, homogeneously analysed and for which elements from different nucleosynthesis channels are available. They allow us to analyse the typical chemical patterns in different regions of the disc of our Galaxy. Clusters are unique tools because their distances and ages can be derived with great accuracy, allowing reliable placement in space and time in our Galaxy. The family of clusters here described belongs to an important but poorly investigated region of our Galaxy: the inner disc.
There are two main results:
-
1.
As noted in Magrini et al. (2014),there are real differences in the abundance patterns of the clusterseven if they are located at similar distances from the Galacticcentre. One very noticeable difference is the higher α-enhancement of NGC 6705, which might suggest an inner birthplace for this cluster, or, as we see in next sections, possibly an episode of local enrichment.
-
2.
Despite the small but real differences, a common behaviour of the family of old/intermediate-age clusters analysed so far can be seen; they are almost solar in most of the iron-peak and heavy α-elements, while they are enhanced in the light α-elements O, Mg, and Si. We discuss in the next sections to what extent this enhancement is real or due to NLTE effects.
5.1. Is the α-enhancement seen in the inner clusters real?
Three of the clusters in our inner disc sample have [O/Fe], [Si/Fe], and [Mg/Fe] ≥ 0.1 dex. The [Si/Fe] ratio is very high in NGC 6705, while it is around 0.1 dex for the remaining two clusters. First, it is necessary to investigate if the O, Mg, and Si enhancements are real or related to other effects, such as NLTE. All stars analysed in the four clusters are indeed giant stars, and thus NLTE might affect the abundance of some elements. If the NLTE effect is large, the comparison between the giant stars in clusters and the Sun or the Milky Way field solar-neighbourhood dwarf stars can be misleading. We clarify that our definition of NLTE is totally empirical: in the following we indicate as a “NLTE effect” the sum of all un-modelled physics plus potential model biases.
Oxygen. In Be 81, and in the other clusters and field stars, oxygen abundance was determined from the forbidden [O I] line at 6300.31 Å, as described in detail by Tautvaišienė et al. (2015). Lines of [O I] are considered to be very good indicators of oxygen abundances because they are not only insensitive to NLTE effects, but also give similar oxygen abundance results with 3D and 1D model atmospheres (cf. Asplund et al. 2004 and Pereira et al. 2009). Thus, the comparison of giants and dwarfs should not be affected by any such bias.
Magnesium. Mg is less affected than the other α-elements by model dependent differences in the yields of SN II with progenitors of different masses. Therefore, it constitutes one of the best tracers of the evolutionary timescales of the Galaxy. However, the NLTE effects in Mg are not easily quantifiable. From the theoretical point of view, Merle et al. (2011) computed the Mg I abundance corrections for several lines in the optical interval as functions of atmospheric parameters for cool evolved stars. Their NLTE corrections are computed as a function of stellar parameters and metallicity. At solar metallicity for a star with Teff = 5000 K and log g = 2 the NLTE corrections of different Mg I lines vary from −0.03 to 0.18 dex. From the observational point of view, several works in the literature have analysed dwarf and giant stars in the same cluster to investigate if there is any difference in the abundance ratios related to the type of the star analysed. A study of dwarf and giant stars in the Hyades open cluster (Schuler et al. 2009) showed that Mg abundances are overabundant by 0.2–0.5 dex in the giants relative to the dwarfs. Such large enhancements are not predicted by stellar models, and they explained them as likely being due to NLTE effects, significantly larger than those available in the literature. A similar result was obtained by Yong et al. (2005) comparing abundances in dwarf and giant stars in M 67. However, Pasquini et al. (2004) and Blanco-Cuaresma et al. (2015) found contrary results for the open cluster IC 4651 ([Fe/H] = 0.11).
Silicon. The few studies of its NLTE behaviour indicate strong effects only in low metallicity stars. Recently, Bergemann et al. (2013) found that the effect of NLTE is a reduction of the number densities of atoms in a minority ionisation stage, thus, for a given abundance, the line equivalent width is smaller than that given by LTE. This effect is larger in metal-poor stars. At solar metallicity, the opposite behaviour is encountered, showing negative NLTE effects such that the LTE abundance is higher than the NLTE one.
In order to better investigate the question for Si and Mg, we test the nature of the α-enhancement in our results using the results for one of the Gaia-ESO calibrator open clusters (see Pancino et al. in prep.), M 67. The data for M 67 are, indeed, perfectly suited because they contain stars belonging to different evolutionary phases that have been analysed with the same tools and methods as the inner-disc clusters and field stars of idr2/3.
In Fig. 11, we show [Mg/Fe] and [Si/Fe] vs. Teff in member stars of M 67 analysed in idr2/3. The ±1σ and ±2σ regions from the average are shown. The sample of stars with Teff> 5700 K includes TO, MS, and SG stars. The stars with Teff< 5700 K are instead giant stars with log g ranging from 2.97 to 3.46. We note that for both elements a slight overabundance in giants is present with respect to TO, MS, and SG stars. The average [Mg/H] for stars with Teff< 5700 K is 0.04 ± 0.02 and for [Si/H] is 0.08 ± 0.01, while the stars with Teff> 5700 K have average values of −0.09 ± 0.04 and −0.03 ± 0.03 from [Mg/H] and [Si/H], respectively. We can use these differences to roughly estimate the level of NLTE effects: 0.13 dex and 0.11 dex for [Mg/H] and [Si/H], respectively. If we consider the abundances with respect to iron, the NLTE effects are 0.12 dex and 0.09 dex for [Mg/Fe] and [Si/Fe], respectively.
It is worth mentioning that M 67 has been used by Önehag et al. (2014) to study the effect of selective atomic diffusion of chemical elements in stars in different evolutionary phases. More specifically, they investigated the detailed chemical composition of 14 stars located on the main sequence, at the turn-off point, and on the early sub-giant branch and found that the heavy-element abundances are reduced in the TO and main sequence (typically by ≤0.05 dex) when compared to the abundances of the sub-giants (even lower for ratios of metals). We note that the effect of diffusion in our data is circumscribed in the ± 1σ area.
We finally mention that Blanco-Cuaresma et al. (2015) have studied the differences in the abundance ratios in cluster giant and dwarf member stars. They show the chemical differences between dwarfs and giants in several clusters, including IC 4651, M 67, NGC 2447, and NGC 3680. They find that, on average, [Mg/Fe] and [Si/Fe] are enhanced in giant stars by ~0.06 dex and ~0.10 dex with respect to dwarfs. In particular, the effect they measured for M 67 is 0.07 and 0.10 dex in [Mg/Fe] and [Si/Fe], respectively, thus comparable with what we find with the Gaia-ESO recommended values.
In conclusion, no NLTE has to be applied for oxygen, and the three clusters for which this element has been measured show a [O/Fe] ≥ 0.1 dex. The enhancement in [Mg/Fe] and [Si/Fe] observed in the inner-disc/bulge population and in Be 81, Tr 20, and NGC 4815 with respect to the solar neighbourhood dwarf population can be related to different evolutionary phases of the observed stars (see Figs. 9 and 10). On the other hand, [Mg/Fe] and [Si/Fe] are remarkably higher in NGC 6705 than in the other inner-disc clusters and in the field populations. We can exclude that this enhancement is due solely to NLTE effects. The abundance ratios of NGC 6705 clearly indicate an α-enhanced population whose origin needs to be investigated.
5.2. Comparing the inner clusters with high-α metal-rich population
Once shown that the enhancements of [O/Fe], [Si/Fe], and [Mg/Fe] in NGC 6705 are genuine, we can try to compare its chemical pattern with those of a population that presents similar chemical characteristics, namely the high-α metal-rich stars (HAMR) with [Fe/H] = 0.03 ± 0.12.
The HAMR stars–whose α-enhancement, computed with Mg, Si, and Ti, is in the range 0–0.2 dex– are rare objects in the solar neighbourhood, and their exact nature is still under debate. Adibekyan et al. (2011) identified these stars and suggested that they may form a distinct family, based on a “gap” in metallicity between HAMR and classical thick disc stars. Other authors also confirmed the existence of the α-enhanced metal-rich stars but with different conclusions regarding their origin (e.g. Gazzano et al. 2013; Nissen et al. 2014; Bensby et al. 2014; Adibekyan et al. 2014). A comparison of the main physical properties of HAMR stars with those from the simulations of Roškar et al. (2012) suggests that the HAMR stars might have their origin in the central part of the Galaxy, within about 2 kpc (Adibekyan et al. 2013). Looking at Figs. 9 and 10, one would claim that stars in Be 81 and in the inner-disc clusters are similar to the α-enhanced metal-rich stars discovered in the solar neighbourhood. They have indeed solar or super-solar [Fe/H] (see Magrini et al. 2014, and Table 4 for Be 81) and they are, in general, enhanced in the light α-elements, such as O, Mg, and Si. At a first glance, the chemical patterns of Be 81 shown in Fig. 9 nicely resemble the metal-rich α-enhanced population identified in the solar neighbourhood. We note that HAMR stars studied by Adibekyan et al. (2011) have Teff> 4900 K and only five stars have log g< 3.8 dex, therefore their sample is essentially composed of dwarf stars. If confirmed, the existence of HAMR stars in an open cluster at 5–7 kpc from the Galactic centre, might put important constraints on the stellar migration and on the birthplace of HAMR stars.
However, we should consider that there are important differences between HAMR stars and the inner disc clusters. First, there is a large difference in age between the two populations. The age of HAMR stars in the solar neighbourhood is not precisely determined. However, they are on average believed to be quite old. Adibekyan et al. (2011) showed that the HAMR and thick disc family stars have almost the same age; they are on average older than thin disc stars by about 3 Gyr. At the same time, they found that HAMR stars are quite different from the thick disc stars; they are more metal-rich and with different chemical patterns. They also show properties in common with metal-rich bulge stars (e.g. Fulbright et al. 2007; Bensby et al. 2014; Ness et al. 2013). On the other hand, the age and distance of NGC 6705 are well known by isochrone fitting (see Table 7), and it is much younger.
Fig. 12 [Fe/H] vs. RGC (upper panel) and RGC vs. [α/Fe] (lower panel) in the chemical evolution model of Magrini et al. (2009). The three curves are the curve at present time (blue), ~5 Gyr ago (red) and ~10 Gyr ago (black). The green shaded area indicates where and when HAMR stars might be born, the black filled circles are the four open clusters and the yellow squares are the median values of the inner disc field stars in the metallicity range from −0.1 to +0.33 dex, whose [α/Fe] are computed removing the NLTE effects. Finally, the red triangles are the median values of the solar neighbourhood stars in the metallicity range from −0.1 to +0.33 dex. |
Second, in a classical chemical evolution view of the Milky Way disc (e.g. Tosi 1988; Matteucci & Francois 1989; Chiappini et al. 1997, 2000; Portinari et al. 1998; Boissier & Prantzos 1999) we do not expect relatively young stars born in the inner disc to show an α-enhancement. It is instructive to use a recent and simple chemical evolution model that does not include the effect of radial migration to study the most likely birthplace of the HAMR stars. In Fig. 12 we show a plot of [Fe/H] and [α/Fe] vs. RGC in the chemical evolution model of Magrini et al. (2009). The three curves show the prediction for [Fe/H] and [α/Fe] as a function of RGC at present, 5 Gyr ago, and 10 Gyr ago. The model predicts that the formation of a family of α-enhanced metal-rich stars is very limited in space and in time. They formed in the inner part of the disc and during the first epochs of disc formation (black curve), then the contribution of SNIa decreased the [α/Fe] ratio to solar values. Thus, following stellar populations have higher [Fe/H] but lower [α/Fe] than the oldest ones. Similar results can be obtained even with more sophisticated models, those by e.g. Minchev et al. (2013), where the combination of high [α/Fe] and young ages is not predicted. For instance, in Fig. 2 of Minchev et al. (2013), stars with [O/Fe] = 0.2 are older than 7 Gyr, whatever their birth location. Stars with slightly smaller [O/Fe] can be younger on the condition that they are born in the very outer disc, but then they would be metal poor.
We also plot [Fe/H] and [α/Fe] for the four clusters and for the solar neighbourhood and inner-disc samples in the same metallicity range of the clusters. We compute [α/Fe] considering [O/Fe], [Mg/Fe], and [Si/Fe] the latter two corrected for the NLTE effects. For the field populations we assumed a distance of 8 ± 1 kpc for the solar neighbourhood stars, and of 5 ± 2 kpc for the inner-disc sample. The uncertainties on the x-axes reflect their interval of distances.
The four clusters are in good agreement with the model prediction in the [Fe/H] vs. RGC plane; however, in the [α/Fe] vs. RGC plane NGC 6705 is displaced with respect to what is expected for its age. Stars with the age of ~0.3 Gyr should have solar or sub-solar [α/Fe] given their age and place of birth (blue curve of Fig. 12) and thus their α-enhancement is totally unexpected. The location of Be 81 is also a bit surprising. As we discussed above, this cluster does not have a high [α/Fe] per se. However, if we consider its [α/Fe] in conjunction with its [Fe/H], we find it a bit higher than expected. On the other hand, [α/Fe] in Tr 20 and NGC 4815 are consistent within the errors with the model curves (in between the red and blue curves in Fig. 12).
Therefore, the stars in NGC 4815 and Tr 20 have a [α/Fe] that is fully consistent with their age and location in the disc as indicated by the predictions of the chemical evolution model of Magrini et al. (2009). Be 81 has a slightly higher [α/Fe] than expected for its [Fe/H] and RGC. However, the most striking case is that of NGC 6705, which shows a genuine enhancement of [α/Fe] unlikely related to the HAMR stellar population in the solar neighbourhood since they were born in a completely different epoch. Owing to the young age of NGC 6705, its α-enhancement cannot be simply explained by disc chemical evolution, and the hypothesis of Magrini et al. (2014) of radial migration from the inner disc is also unlikely given the young age of this cluster.
Very recently, several studies revealed the presence of relatively young α-enhanced stars. Chiappini et al. (2015) report the discovery of a group of apparently young CoRoT red-giant stars exhibiting enhanced [α/Fe] abundance ratios. These stars are more numerous in their inner-disc sample. Martig et al. (2014) has analysed a sample of giants with seismic parameters from Kepler and also found a small percentage of stars with [α/Fe] > 0.2 dex and ages below 4 Gyr. As in the case of our inner disc clusters, the existence of this kind of star cannot be explained in the framework of standard chemical evolution models of the Milky Way. Their presence might indicate that the chemical-enrichment history of the Galactic disc is more complex. Chiappini et al. (2015) analyse some possible interpretations of the nature of these stars: i) they might have been formed close to the end of the Galactic bar, near co-rotation where gas stays inert for longer times and where the mass return from older inner-disc stellar generations is expected to be highest with the net effect of in situ gas dilution; ii) they might originate from a recent gas-accretion event. However, there is a substantial difference between the stars discussed in Chiappini et al. (2015) and the inner-disc open clusters: Chiappini’s stars tend to be more metal poor than our clusters, with [Fe/H] ranging from ~−0.6 to ~0.1 dex. Might the inner-disc open cluster be a younger tail of the “young” metal-poor α-enhanced population? Is the high [α/Fe] of both populations related to their place of birth, near to co-rotation?
In the next section, we try to estimate another possibility: the effect of a very local enrichment in the molecular cloud from which NGC 6705 was born to explain its chemical pattern.
5.3. The local enrichment hypothesis: the effect of a supernova explosion
Can a single SNII explosion explain the chemical pattern of NGC 6705 with its high [α/Fe]? Using Gaia-ESO data, Bergemann et al. (2014) found a large scatter in [Fe/H] and [Mg/Fe] of stars at each age towards the inner disc. Malinie et al. (1993) explained the observed scatter at each Galactocentric radius and age as being due to the fact that the disc contains chemical inhomogeneities that can avoid remixing long enough. This allows the activation of the star formation process in both enriched and unenriched portions, with typical timescales of 108–109 yr. Gilmore et al. (1990) suggested that the scatter could be a consequence of self-enrichment in giant molecular cloud complexes. Their lifetime is shorter than the timescale of Galactic evolution, but much longer than the lifetime of the massive stars. Therefore, self-enrichment could occur inside a giant molecular cloud, and after the disruption of the cloud the metallicity inhomogeneity survives. It spreads only through a Galactic annulus, but it can be easily detected if the cluster formed from the molecular cloud is still observable.
Cluster enrichment by SN II using yields of WW95.
Starting from the assumption that the interstellar medium from which NGC 6705 was born had a mass comparable to its present time mass (~3500 M⊙, Cantat-Gaudin et al. 2014), we compute the enrichment due to the explosion of a massive star of 15, 18, or 25 M⊙ (yields from Woosley & Weaver 1995) in the hypothesis of a uniform and instantaneous recycling of the medium. The enrichment has been computed assuming solar initial composition (Grevesse et al. 2007).
In Table 9, we present the results of this calculation: the explosion of a massive star with a mass 15–18 M⊙ might be able to explain the local enrichment of NGC 6705. The yields of Woosley & Weaver (1995) can also explain the higher enrichment in [Si/Fe] with respect to the other elements. The hypothesis of local enrichment is also supported by Roy & Kunth (1995) where characteristic timescales for homogenisation are short, 106–108 yr on scales of 100–1000 pc. The number of stars with a mass in the range 15–18 M⊙ computed with the initial mass function of Cantat-Gaudin et al. (2014) in the whole area of the cluster is about two. Given the present total mass of NGC 6705, the expected number of stars with 15 <M< 18 M⊙ is thus consistent with an early enrichment by one or more type II supernovae with such massive progenitors. However, the hypothesis that the SN was a cluster member implies at least three assumptions, namely i) a prolonged star formation (at least 107 Myr: timescales of SNII with progenitors of 15 <M< 18 M⊙), where massive stars have formed first and have enriched the medium from which the lower mass stars have subsequently formed; ii) that the explosion of the SN has not dissipated the gas from which the low mass stars have formed; and iii) that no low mass stars have formed from the pristine (not yet enriched) gas, or that those stars have escaped the clusters. Similar problems are also faced in globular cluster formation and in their self-enrichment.
6. Summary and conclusions
We analyse the Gaia-ESO results (idr2/3) for four inner disc clusters, including the old cluster Be 81 which has not been studied in the Gaia-ESO survey so far. Using the improvements obtained in the metallicity measurement and the membership with RV determinations, we derive age, distance, and average reddening of Be 81 by using the classical approach of isochrone fitting. The best-fit isochrone gives an age of ~1 Gyr, a reddening E(B−V) =0.83 and a Galactocentric distance of 5.45 kpc. We compare the chemical pattern of Be 81 and of the other three old open clusters of idr2/3 with the solar-neighbourhood and inner-disc/bulge stars. We find that the inner-disc/bulge giant stars and three open clusters show an enhancement in the light α-elements (O, Si, Mg). We use the calibrator cluster M 67 to infer how much of this enhancement is related to NLTE effects, finding that the enhancement of [Mg/Fe] and [Si/Fe] in giants with respect to dwarfs is ~0.10 dex. Applying these empirical NLTE corrections to the chemical patterns of open clusters, we find that Tr 20, NGC 4815, and Be 81 do not show any remarkable α-enhancement. On the other hand, NGC 6705 still shows an uncommon pattern, with a significant α-enhancement. We compare the cluster [α/Fe] with those of the HAMR stars, finding that they belong to very different epochs of the lifetime of the Galaxy. We also relate them to the recently discovered population of young (≤4 Gyr) and α-enhanced stars Martig et al. 2014; Chiappini et al. 2015 mainly located towards the Galactic centre. The inner-disc clusters, with ages between ~0.3 to ~1.5 Gyr, and the young α-enhanced stars might belong to an evolutionary sequence of the same population, born in the inner part of the disc. The common chemical features of some of the clusters with the young α-enhanced stars might indicate peculiar conditions for the star formation near to the co-rotation. In our discussion, we also tentatively explain the very high values of [O/Fe], [Mg/Fe], and [Si/Fe] in NGC 6705 with an episode of local enrichment due to a type II supernova in the mass range 15–18 M⊙, as, contemplated by e.g. Malinie et al. (1993) to explain the observed inhomogeneities at each Galactocentric radius.
In the whole paper the adopted RGC, ⊙ is 8 kpc (see Malkin 2013).
id of the original photometric catalogue of Donati et al. (2014a), also available at the CDS.
Acknowledgments
The Authors thank an anonymous referee whose comments and suggestions improved the quality of the paper. The results presented here benefited from discussions in three Gaia-ESO workshops supported by the ESF (European Science Foundation) through the GREAT (Gaia Research for European Astronomy Training) Research Network Program (Science meetings 3855, 4127 and 4415). This work was partially supported by the Gaia Research for European Astronomy Training (GREAT-ITN) Marie Curie network, funded through the European Union Seventh Framework Programme [FP7/2007-2013] under grant agreement n. 264895. 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 the Local Group Galaxies” (prot. 2010LY5N2T). V.A. acknowledges the support from the Fundação para a Ciência e a Tecnologia, FCT (Portugal) in the form of the fellowship SFRH/BPD/70574/2010. S.V. gratefully acknowledges the support provided by Fondecyt. reg. 1130721. T.B. was supported by the project grant “The New Milky” from the Knut and Alice Wallenberg foundation. U.H. acknowledges support from the Swedish National Space Board (SNSB/Rymdstyrelsen). G.T. acknowledges the support from the Research Council of Lithuania (grant No. MIP-082/2015). This research has made use of the SIMBAD database, operated at CDS, Strasbourg, France.
References
- Adibekyan, V. Z., Santos, N. C., Sousa, S. G., & Israelian, G. 2011, A&A, 535, L11 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Adibekyan, V. Z., Figueira, P., Santos, N. C., et al. 2013, A&A, 554, A44 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Adibekyan, V. Z., González Hernández, J. I., Delgado Mena, E., et al. 2014, A&A, 564, L15 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Asplund, M., Grevesse, N., Sauval, A. J., Allen de Prieto, C., & Kiselman, D. 2004, A&A, 417, 751 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Bensby, T., Feltzing, S., & Oey, M. S. 2014, A&A, 562, A71 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Bergemann, M., Kudritzki, R.-P., Würl, M., et al. 2013, ApJ, 764, 115 [NASA ADS] [CrossRef] [Google Scholar]
- Bergemann, M., Ruchti, G. R., Serenelli, A., et al. 2014, A&A, 565, A89 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Blanco-Cuaresma, S., Soubiran, C., Heiter, U., et al. 2015, A&A, 577, A47 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Boissier, S., & Prantzos, N. 1999, MNRAS, 307, 857 [NASA ADS] [CrossRef] [Google Scholar]
- Bragaglia, A., & Tosi, M. 2006, AJ, 131, 1544 [NASA ADS] [CrossRef] [Google Scholar]
- Bressan, A., Marigo, P., Girardi, L., et al. 2012, MNRAS, 427, 12739 [NASA ADS] [CrossRef] [Google Scholar]
- Cantat-Gaudin, T., Vallenari, A., Zaggia, S., et al. 2014, A&A, 569, A17 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Cardelli, J. A., Clayton, G. C., & Mathis, J. S. 1989, ApJ, 345, 245 [NASA ADS] [CrossRef] [Google Scholar]
- Carraro, G., Monaco, L., & Villanova, S. 2014a, A&A, 568, A86 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Carraro, G., Giorgi, E. E., Costa, E., & Vázquez, R. A. 2014b, MNRAS, 441, L36 [NASA ADS] [CrossRef] [Google Scholar]
- Carraro, G., Villanova, S., Monaco, L., et al. 2014c, A&A, 562, A39 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Carretta, E., Bragaglia, A., Gratton, R. G., & Tosi, M. 2005, A&A, 441, 131 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Carretta, E., Bragaglia, A., & Gratton, R. G. 2007, A&A, 473, 129 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Charbonnel, C., & Lagarde, N. 2010, A&A, 522, A10 [NASA ADS] [CrossRef] [EDP Sciences] [MathSciNet] [PubMed] [Google Scholar]
- Chiappini, C., Matteucci, F., & Gratton, R. 1997, ApJ, 477, 765 [NASA ADS] [CrossRef] [Google Scholar]
- Chiappini, C., Matteucci, F., & Padoan, P. 2000, ApJ, 528, 711 [NASA ADS] [CrossRef] [Google Scholar]
- Chiappini, C., Anders, F., Rodrigues, T. S., et al. 2015, A&A, 576, L12 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Damiani, F., Prisinzano, L., Micela, G., et al. 2014, A&A, 566, A50 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Donati, P., Beccari, G., Bragaglia, A., Cignoni, M., & Tosi, M. 2014a, MNRAS, 437, 1241 [NASA ADS] [CrossRef] [Google Scholar]
- Donati, P., Cantat Gaudin, T., Bragaglia, A., et al. 2014b, A&A, 561, A94 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Dotter, A., Chaboyer, B., Jevremović, D., et al. 2008, ApJS, 178, 89 [NASA ADS] [CrossRef] [Google Scholar]
- Freeman, K. C. 1970, ApJ, 160, 811 [NASA ADS] [CrossRef] [Google Scholar]
- Friel, E. D. 1995, ARA&A, 33, 381 [NASA ADS] [CrossRef] [Google Scholar]
- Friel, E. D., Donati, P., Bragaglia, A., et al. 2014, A&A, 563, A117 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Fulbright, J. P., McWilliam, A., & Rich, R. M. 2007, ApJ, 661, 1152 [NASA ADS] [CrossRef] [Google Scholar]
- Gazzano, J.-C., Kordopatis, G., Deleuil, M., et al. 2013, A&A, 550, A125 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Gilmore, G., King, I. R., van der Kruit, P. C., & Buser, R. 1990, Science, 250, 703 [NASA ADS] [Google Scholar]
- Gilmore, G., Randich, S., Asplund, M., et al. 2012, The Messenger, 147, 25 [NASA ADS] [Google Scholar]
- Grevesse, N., Asplund, M., & Sauval, A. J. 2007, Space Sci. Rev., 130, 105 [Google Scholar]
- Grieco, V., Matteucci, F., Pipino, A., & Cescutti, G. 2012, A&A, 548, A60 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Hayes, C. R., & Friel, E. D. 2014, AJ, 147, 69 [NASA ADS] [CrossRef] [Google Scholar]
- Karakas, A. I., & Lattanzio, J. C. 2014, PASA, 31, e030 [NASA ADS] [CrossRef] [Google Scholar]
- Iwamoto, K., Brachwitz, F., Nomoto, K., et al. 1999, ApJS, 125, 439 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
- Magrini, L., Sestito, P., Randich, S., & Galli, D. 2009, A&A, 494, 95 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Magrini, L., Randich, S., Zoccali, M., et al. 2010, A&A, 523, A11 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Magrini, L., Randich, S., Romano, D., et al. 2014, A&A, 563, A44 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Malinie, G., Hartmann, D. H., Clayton, D. D., & Mathews, G. J. 1993, ApJ, 413, 633 [NASA ADS] [CrossRef] [Google Scholar]
- Martig, M., Rix, H.-W., Silva Aguirre, V., et al. 2014, MNRAS, 451, 2230 [Google Scholar]
- Matteucci, F., & Brocato, E. 1990, ApJ, 365, 539 [NASA ADS] [CrossRef] [Google Scholar]
- Matteucci, F., & Francois, P. 1989, MNRAS, 239, 885 [NASA ADS] [CrossRef] [Google Scholar]
- Malkin, Z. 2013, IAU Symp., 289, 406 [Google Scholar]
- Merle, T., Thévenin, F., Pichon, B., & Bigot, L. 2011, MNRAS, 418, 863 [NASA ADS] [CrossRef] [Google Scholar]
- Minchev, I., Chiappini, C., & Martig, M. 2013, A&A, 558, A9 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Ness, M., Freeman, K., Athanassoula, E., et al. 2013, MNRAS, 432, 2092 [NASA ADS] [CrossRef] [Google Scholar]
- Nissen, P. E., Chen, Y. Q., Carigi, L., Schuster, W. J., & Zhao, G. 2014, A&A, 568, A25 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- O’Donnell, J. E. 1994, ApJ, 422, 158 [NASA ADS] [CrossRef] [Google Scholar]
- Önehag, A., Gustafsson, B., & Korn, A. 2014, A&A, 562, A102 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Pasquini, L., Randich, S., Zoccali, M., et al. 2004, A&A, 424, 951 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Pereira, T. M. D., Asplund, M., & Kiselman, D. 2009, A&A, 508, 1403 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Pietrinferni, A., Cassisi, S., Salaris, M., & Castelli, F. 2004, ApJ, 612, 168 [NASA ADS] [CrossRef] [Google Scholar]
- Randich, S., Gilmore, G., & Gaia-ESO Consortium 2013, The Messenger, 154, 47 [NASA ADS] [Google Scholar]
- Portinari, L., Chiosi, C., & Bressan, A. 1998, A&A, 334, 505 [NASA ADS] [Google Scholar]
- Roškar, R., Debattista, V. P., Quinn, T. R., & Wadsley, J. 2012, MNRAS, 426, 2089 [NASA ADS] [CrossRef] [Google Scholar]
- Roy, J.-R., & Kunth, D. 1995, A&A, 294, 432 [NASA ADS] [Google Scholar]
- Sestito, P., Randich, S., & Bragaglia, A. 2007, A&A, 465, 185 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Sestito, P., Bragaglia, A., Randich, S., et al. 2008, A&A, 488, 943 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Sacco, G. G., Morbidelli, L., Franciosini, E., et al. 2014, A&A, 565, A113 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Sagar, R., & Griffiths, W. K. 1998, MNRAS, 299, 1 [NASA ADS] [CrossRef] [Google Scholar]
- Schuler, S. C., King, J. R., & The, L.-S. 2009, ApJ, 701, 837 [NASA ADS] [CrossRef] [Google Scholar]
- Smiljanic, R., Korn, A. J., Bergemann, M., et al. 2014, A&A, 570, A122 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Tautvaišienė, G., Drazdauskas, A., Mikolaitis, Š., et al. 2015, A&A, 573, A55 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Tosi, M. 1988, A&A, 197, 33 [NASA ADS] [Google Scholar]
- VandenBerg, D. A., Bergbusch, P. A., & Dowler, P. D. 2006, ApJS, 162, 375 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
- van der Kruit, P. C. 2002, in The Dynamics, Structure and History of Galaxies: A Workshop in Honour of Professor Ken Freeman, 273, 7 [NASA ADS] [Google Scholar]
- Warren, S. R., & Cole, A. A. 2009, MNRAS, 393, 272 [NASA ADS] [CrossRef] [Google Scholar]
- Woosley, S. E., & Weaver, T. A. 1995, ApJS, 101, 181 [NASA ADS] [CrossRef] [Google Scholar]
- Woosley, S. E., Heger, A., & Weaver, T. A. 2002, Rev. Mod. Phys., 74, 1015 [NASA ADS] [CrossRef] [Google Scholar]
- Yong, D., Carney, B. W., & Teixera de Almeida, M. L. 2005, AJ, 130, 597 [NASA ADS] [CrossRef] [Google Scholar]
- Yong, D., Carney, B. W., & Friel, E. D. 2012, AJ, 144, 95 [NASA ADS] [CrossRef] [Google Scholar]
All Tables
Abundances in the 12+log (X/H) form of the UVES stars in Be 81 from idr2/3 recommended values.
Results, errors, and estimated systematic uncertainties using different evolutionary models with [Fe/H] ~ 0.2.
All Figures
Fig. 1 Upper left panel: in grey the Be 81 field of view of the photometric data in Donati et al. (2014a), obtained with LBC at LBT; in black the stars observed by the Gaia-ESO survey, all inside a radius of 6 arcmin. Lower left panel: distribution of the radial velocities vs. distance from the cluster centre. The red line is the cluster systemic radial velocity. Right panel: CMD of the inner 6 arcmin, with stars observed indicated by larger symbols (orange squares HR09B; red crosses HR15N; blue circles UVES). |
|
In the text |
Fig. 2 Distribution of radial velocities of the stars observed in the field of Be 81. The solid line is a Gaussian fit of the distribution of the member stars. |
|
In the text |
Fig. 3 [Fe/H] vs. atmospheric parameters of UVES stars in Be 81. The dotted lines mark the 1σ region, while the dashed lines indicate the 2σ area. The typical error bar on each star are shown on the left side of each panel. |
|
In the text |
Fig. 4 Distribution of [Fe/H] of the stars observed in the field of Be 81. The cyan rectangles show the GIRAFFE metallicities of RV member stars, while the blue ones indicate the UVES metallicities (including star 19015978-0028183). |
|
In the text |
Fig. 5 HR diagram with the PARSEC isochrone for an age of 0.98 Gyr and a metallicity Z = 0.025 (red curve) and the recommended Gaia-ESO parameters of UVES stars (green squares) and of GIRAFFE stars (blue circles). The stars that are members on the basis of their radial velocity but with [Fe/H] outside 1σ from the cluster average metallicity are shown with empty symbols. |
|
In the text |
Fig. 6 Abundance ratios of stars in Be 81. The dotted lines indicate the median+ σ and median−σ regions in both [El/Fe] and [Fe/H]. The typical error bar on each star are shown on the left side of each panel. |
|
In the text |
Fig. 7 CMD obtained for stars inside 3′ using the photometry from Donati et al. (2014a) and best-fit isochrone for different evolutionary models (from left to right: PARSEC, BASTI, DARTMOUTH, and Victoria-Regina). GES target non-members are plotted as black crosses, while members are highlighted with orange (GIRAFFE) and blue (UVES) points. See Table 6 for the adopted parameters for the isochrone fitting. |
|
In the text |
Fig. 8 Same as Fig. 7 but in the V,V−I CMD. |
|
In the text |
Fig. 9 Abundance ratios of the idr2/3 old and intermediate-age open clusters located in the inner disc as a function of the atomic number. In black are the abundance ratios of Be 81, in yellow NGC 6705, in green NGC 4815, and in red Tr 20. The dashed line marks the solar abundance, and the shaded area the 1σ regions with σ = 0.1 dex. |
|
In the text |
Fig. 10 Abundance ratios vs. atomic numbers of the idr2/3 stars in the solar neighbourhood and in the inner disc/bulge (selected in the [Fe/H] range from −0.10 to +0.33 dex, i.e. including the metallicity ranges of the four clusters) compared with the chemical pattern of NGC 6705, Trumpler 20, NGC 4815, and Be 81. The continuous lines mark the solar abundances, and the shaded area the 1σ region with σ = 0.1 dex. |
|
In the text |
Fig. 11 Abundance ratios [Si/H] and [Mg/H] vs. Teff of member stars in M 67 from idr2/3 analysis. The continuous lines mark the ±1σ from the average value, whereas the dotted lines indicate the ±2σ region. |
|
In the text |
Fig. 12 [Fe/H] vs. RGC (upper panel) and RGC vs. [α/Fe] (lower panel) in the chemical evolution model of Magrini et al. (2009). The three curves are the curve at present time (blue), ~5 Gyr ago (red) and ~10 Gyr ago (black). The green shaded area indicates where and when HAMR stars might be born, the black filled circles are the four open clusters and the yellow squares are the median values of the inner disc field stars in the metallicity range from −0.1 to +0.33 dex, whose [α/Fe] are computed removing the NLTE effects. Finally, the red triangles are the median values of the solar neighbourhood stars in the metallicity range from −0.1 to +0.33 dex. |
|
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.