Issue |
A&A
Volume 616, August 2018
|
|
---|---|---|
Article Number | A142 | |
Number of page(s) | 8 | |
Section | Galactic structure, stellar clusters and populations | |
DOI | https://doi.org/10.1051/0004-6361/201833255 | |
Published online | 05 September 2018 |
Abundance of zinc in the red giants of Galactic globular cluster 47 Tucanae
1
Astronomical Observatory of Vilnius University, Saul˙etekio al. 3, Vilnius 10257, Lithuania
e-mail algimantas.cerniauskas@tfai.vu.lt
2
GEPI, Observatoire de Paris, Université PSL, CNRS, 5 Place Jules Janssen, 92190 Meudon, France
3
Zentrum für Astronomie der Universität Heidelberg, Landessternwarte, Königstuhl 12, 69117 Heidelberg, Germany
4
Leibniz-Institut für Astrophysik Potsdam, An der Sternwarte 16, 14482 Potsdam, Germany
Received:
18
April
2018
Accepted:
24
May
2018
Aims. We investigate possible relations between the abundances of zinc and the light elements sodium, magnesium, and potassium in the atmospheres of red giant branch (RGB) stars of the Galactic globular cluster 47 Tuc and study connections between the chemical composition and dynamical properties of the cluster RGB stars.
Methods. The abundance of zinc was determined in 27 RGB stars of 47 Tuc using 1D local thermal equilibrium (LTE) synthetic line profile fitting to the high-resolution 2dF/HERMES spectra obtained with the Anglo-Australian Telescope (AAT). Synthetic spectra used in the fitting procedure were computed with the SYNTHE code and 1D ATLAS9 stellar model atmospheres.
Results. The average 1D LTE zinc-to-iron abundance ratio and its RMS variations due to star-to-star abundance spread determined in the sample of 27 RGB stars is 〈[Zn/Fe]〉1D LTE = 0.11 ± 0.09. We did not detect any statistically significant relations between the abundances of zinc and those of light elements. Neither did we find any significant correlation or anticorrelation between the zinc abundance in individual stars and their projected distance from the cluster center. Finally, no statistically significant relation between the absolute radial velocities of individual stars and the abundance of zinc in their atmospheres was detected. The obtained average [Zn/Fe]1DLTE ratio agrees well with those determined in this cluster in earlier studies and nearly coincides with that of Galactic field stars at this metallicity. All these results suggest that nucleosynthesis of zinc and light elements proceeded in separate, unrelated pathways in 47 Tuc.
Key words: globular clusters: individual: NGC 104 / stars: late-type / stars: atmospheres / stars: abundances / techniques: spectroscopic / convection
© ESO 2018
1. Introduction
Zinc has traditionally been considered as a transition element between the Fe-peak and light s-process species. It is thought that zinc is produced through a variety of channels: explosive nucleosynthesis in type II supernovae (SNe II) and type Ia SNe (SNe Ia), hypernovae, the main and the weak components of the s-process (see, e.g., Sneden et al. 1991; Matteucci et al. 1993; Mishenina et al. 2002; Kobayashi et al. 2006).
Observational and theoretical studies of the Galactic evolution of zinc performed during the past decade have confirmed that zinc is synthesized through multiple channels in stars that belong to different Galactic populations. These investigations have shown that the [Zn/Fe] ratio reaches up to ≈+0.5 in halo stars with [Fe/H] < −3.0 (Cayrel et al. 2004). In the thin and thick disks, however, its value decreases with increasing metallicity and becomes sub-solar at [Fe/H] ≈ 0.0, with significant scatter in [Zn/Fe] ratios at this metallicity in giant stars that most likely belong to the thin disk (e.g., Duffau et al. 2017). Nissen & Schuster (2011) have found that solar neighborhood stars with low and high-element abundances were characterized with correspondingly low and high [Zn/Fe] ratios. More recently, Duffau et al. (2017) have shown that stars located at galactocentric distances <7 kpc seem be to be depleted in zinc, down to [Zn/Fe] ≈ −0.3, which agrees with the values found in the Galactic bulge stars (Barbuy et al. 2015; da Silveira et al. 2018). These and other observations suggest a related origin of zinc and α-elements (e.g., da Silveira et al. 2018). However, [Zn/Fe] trends at low metallicity point to a production of zinc in high-energy core-collapse supernovae, that is, hypernovae (Kobayashi et al. 2006). The significant scatter observed in [Zn/Fe] ratios at solar metallicity is somewhat more difficult to explain and may require a dilution with essentially zinc-free type SN Ia ejecta (see Duffau et al. 2017, for adiscussion).
That zinc may be produced through several channels makes it a potentially interesting tracer in the context of Galactic globular cluster (GGCs) evolution. One of the scenarios that was suggested to explain relations between the abundances of light elements observed in GGCs (such as the Na–O, Mg–Al anticorrelation; see, e.g., Carretta et al. 2009) assumes that atmospheres of the second-generation stars could have been polluted with elements synthesized in asymptotic giant branch (AGB) stars (Ventura et al. 2001). Since zinc may be produced during the s-process nucleosynthesis that takes place in intermediatemass AGB stars (3–6 M⊙, e.g., Karakas et al. 2009), relations between its abundance and those of light elements (e.g., sodium) may point to AGB stars as possible polluters. Therefore, a study of the zinc abundance in the atmospheres of GGC stars characterized by different sodium abundances may reveal additional clues about the possible polluters of second-generation stars in the GGCs. As far as individual GGCs are concerned, one potentially interesting target is 47 Tuc, which not only exhibits various correlations and anticorrelations between the abundances of light elements, such as Li, O, Na, Mg, and Al (Carretta et al. 2009; D’Orazietal. 2010; Dobrovolskas et al. 2014), but also shows indications for connections between their abundances in the atmospheres of individual cluster stars and kinematical properties of these stars (see, e.g., Kučinskas et al. 2014).
Only a few determinations of zinc abundances in 47 Tuc are available, and they are typically based on relatively small samples of stars. For example, Thygesen et al. (2014) determined the zinc-to-iron ratio in 13 red giant branch (RGB) stars of this cluster, with a median value of [Zn/Fe]1D LTE = 0.26 ± 0.13 and no evidence for an intrinsic zinc abundance spread. More recently, based on the study of 19 RGB stars in 47 Tuc, Duffau et al. (2017) obtained a mean value of 〈[Zn/Fe]〉1D NLTE =0.17 ± 0.10 (1D NLTE-LTE abundance corrections for zinc are small, <0.1 dex, see Sect. 2.2; thus, we directly compared local thermal equilibrium (LTE) and NLTE estimates here). This relatively limited information about the zinc abundance in 47 Tuc is rather unfortunate, especially because at the metallicity of 47 Tuc, [Fe/H]1D LTE = − 0.76 (Carretta et al. 2009), Zn I lines located at 472.2 nm and 481.0 nm are quite strong in the spectra of its RGB stars (7–10 pm) and are only weakly blended in their wings, which makes cluster RGB stars good targets for investigating the zinc abundance.
The main goal of this study was therefore to determine the zinc abundance in a larger sample of RGB stars in 47 Tuc and to search for possible relations between the zinc abundance and those of light elements, primarily sodium and magnesium (but also potassium, see Sect. 3 for a discussion). In addition, we also aimed to study possible connections between the kinematical properties of RGB stars and the zinc abundance in their atmospheres.
The paper is structured as follows: spectroscopic data and the determination of the atmospheric parameters of the target stars are presented in Sect. 2.1, while the method we used to determine the zinc abundance and its uncertainties is described in Sect. 2.2. The influence of convection on the Zn I line formation in red giant atmospheres is investigated in Sect. 2.3. The determined zinc abundances, as well as their possible relations with those of light elements and kinematical properties of RGB stars, are discussed in Sect. 3. Finally, the main results obtained in this work are summarized in Sect. 4.
2. Methodology
The abundance of zinc was determined using 1D hydrostatic ATLAS9 model atmospheres and a 1D LTE abundance analysis methodology. A brief description of all steps involved in the abundance analysis is provided below.
2.1. Spectroscopic data and atmospheric parameters of the sample stars
We used the same sample of RGB stars of 47 Tuc as in Černiauskas et al. (2017; hereafter Paper I) as well as the same set of archival high-resolution 2dF/HERMES spectra obtained with the Anglo-Australian Telescope (AAT)1. The spectra were acquired during the science verification phase of the GALAH survey (De Silva et al. 2015) in the wavelength range of 471.5–490.0 nm (HERMES/BLUE) using a spectral resolution of R ~ 28 000 and exposure time of 1200 s. The typical signal-to-noise ratio (per pixel) in this wavelength range was S /N ≈ 50. Further details regarding the reduction and continuum normalization of the observed spectra are provided in Paper I.
The final sample of RGB stars for which zinc abundances were determined contained 27 stars (see Fig. 1). Their effective temperatures and surface gravities were taken from Paper I: the former were determined using photometric data from Bergbusch & Stetson (2009) and Teff − (V − I) calibration of Ramírez & Meléndez (2005), while the latter were estimated using the classical relation between the surface gravity, mass, effective temperature, and luminosity.
![]() |
Fig. 1. Upper part of the color-magnitude diagram of 47 Tuc constructed using VI photometry from Bergbusch & Stetson (2009). RGB stars used in this study are marked as filled green circles. |
2.2. Determination of the 1D LTE zinc abundance in the atmospheres of RGB stars in 47 Tuc
Zinc abundances in individual RGB stars were determined using two Zn I lines located at 472.21 nm and 481.05 nm. Excitation potentials and oscillator strengths of these spectral lines were taken from the VALD-3 database (Piskunov et al. 1995; Kupka et al. 2011). Broadening constants were computed using the classical prescription from Castelli (2005). All atomic line parameters are provided in Table 1.
Atomic parameters of Zn I lines that were used in this study.
The spectra of all RGB stars used in this study were carefully inspected for blends and/or possible contamination of Zn I lines with telluric lines. We verified that Zn I lines were not affected by telluric lines, therefore we did not correct the spectra for telluric absorption. Since Zn I lines were of different quality in the spectra of different stars, the lines were grouped into three classes according to their quality, which was determined by visual inspection:
-
A-class: strong or moderately strong lines with well-resolved line profiles;
-
B-class: lines that were moderately blended or insufficiently resolved in the line wings;
-
C-class: weak, poorly resolved, or significantly blended lines.
Zinc abundances were determined using synthetic 1D LTE spectra that were computed using the SYNTHE package and 1D hydrostatic ATLAS9 model atmospheres (Kurucz 1993), with both SYNTHE and ATLAS9 packages used in the form of the Linux port by Sbordone et al. (2004) and Sbordone (2005). Several typical examples of synthetic spectra fits are provided in Fig. 2. As in Paper I, a fixed value of microturbulence velocity, ξ = 1.5 km s−1, was used to compute all synthetic spectra. The macroturbulence velocity was varied as a free fitting parameter to obtain the best match to the observed Zn I lines, with its value changing in the range from 1 to 6 km s−1 in different stars of our sample. We verified that there is no dependence of the determined zinc abundances on the effective temperature (Fig. 3).
![]() |
Fig. 2. Typical fits of synthetic spectra (red solid lines) to the observed Zn I line profiles (marked by vertical ticks) in the 2dF/HERMES spectra (filled gray circles) of the target RGB stars S1800 (top row; Teff = 4510 K, log g = 1.90) and S167 (bottom row; Teff = 4390 K, log g = 1.60). We also indicate zinc abundances A(Zn) determined from each spectral line, together with their errors (Sect. 2.2). The quality class of the Zn I line is indicated at the top of each panel (see Sect. 2.2 for the line class definition). |
![]() |
Fig. 3. [Zn/Fe] abundance ratios determined in the sample of 27 RGB stars in 47 Tuc and plotted versus the effective temperature of individual stars. |
The solar zinc abundance used in our analysis was . As described in Appendix A, this value was obtained from the same spectral lines that were used in the analysis of zinc abundance in 47 Tuc (Table 1). The obtained value agrees well with the previous determinations of the solar zinc abundance, for example,
obtained using the same Zn I lines by Takeda et al. (2005).
As discussed in Paper I, we did not determine the iron abundance in individual RGB stars because the number of available Fe I lines was small and the quality of the spectra did not allow obtaining reliable estimates of the iron abundance. Instead, the fixed value of [Fe/H]1D LTE = −0.76 taken from Carretta et al. (2009) was used for all stars in this study. To check the validity of this assumption, we determined the iron abundance for three stars with the best-quality spectra in our sample, S167, S1490, and S1563 (Teff = 4385, 4560, 4695 K and log g = 1.60, 2.05, 2.30, respectively). In each case, we could identify six Fe I lines suitable for the abundance analysis. The iron-to-hydrogen ratios determined in these stars were [Fe/H]1D LTE = −0.77, −0.73, and –0.76, respectively, with the mean value, 〈[Fe/H]1D LTE〉 = −0.75, nearly identical to [Fe/H]1D LTE = −0.76 used in this work. As explained above, however, determination of iron abundances for all stars in our sample was not feasible.
We did not correct the obtained zinc abundances for NLTE effects. To our knowledge, the only available source of published 1D NLTE–LTE abundance corrections for the two Zn I lines used in our study is Takeda et al. (2005). However, the authors do not provide corrections at the effective temperatures and gravities bracketing those of our sample stars because of numerical problems in their NLTE computations at these atmospheric parameters. Nevertheless, by extrapolating their available data, we estimated that abundance corrections for our sample stars will be on the order of Δ1D NLTE−LTE ≈ −0.05 dex and are therefore of little significance in the context of present study.
Uncertainties in the determined zinc abundances were evaluated using the same methodology as in Paper I, which is described in Appendix B. Our basic assumption was that these uncertainties were governed by the errors in the determined stellar atmospheric parameters and uncertainties in the line profile fitting. We stress that we did not take into account errors due to uncertainties in the atomic line parameters. Therefore, the obtained error estimates (Table A.2) only provide a lower limit for the uncertainties in the zinc abundances determined from individual Zn I lines because they do not account for other possible sources of systematic errors. Moreover, the individual errors are considered to be uncorrelated.
2.3. 3D–1D LTE zinc abundance corrections
The influence of 3D hydrodynamical effects on the formation of Zn I lines has so far been investigated only for dwarf and subgiant stars (Nissen et al. 2004). The authors have found that the role of these effects was very small, with the resulting 3D–1D LTE abundance corrections significantly below 0.1 dex.
To investigate the role of 3D effects in the formation of Zn I lines in the atmospheres of RGB stars in 47 Tuc, we used 3D hydrodynamical CO5BOLD (Freytag et al. 2012) and 1D hydrostatic LHD (Caffau et al. 2008) model atmospheres. For this, we used models with atmospheric parameters similar to those of the median object in our RGB star sample, Teff ≈ 4490 K, log g = 2.0, and [M/H] = −1.0. This choice is justified because the range covered by the atmospheric parameters of the sample stars is small ( Teff ≈ 300 K, Δ log g ≈ 0.3), which means that the obtained corrections should be applicable to all RGB stars in our sample. In addition to the atmospheric parameters, the CO5BOLD and LHD models share identical chemical composition, opacities, and equation of state (see Paper I for details). Spectral line synthesis computations were carried out with the Linfor3D spectral synthesis package2.
The 3D–1D LTE abundance corrections, Δ3D−1D LTE, were computed for two values of the line equivalent width, W (corresponding to “weak” and “strong” spectral lines) that bracketed the range measured in the observed spectra of the sample RGB stars. Equivalent widths used for the weakest lines were 8.5 pm and 8 pm, while for the strongest lines we used 10.5 pm and 10 pm for Zn 472.21 and Zn 481.05, respectively. The microturbulence velocity in the 3D model atmosphere was determined by applying Method 13 described in Steffen et al. (2013) and was subsequently used in the spectral line synthesis with the LHD model atmospheres (see Paper I).
The obtained 3D–1D LTE abundance corrections are provided in Table 2. For both lines, they do not exceed 0.05 dex and show little dependence on the line strength, which allows us to conclude that the influence of convection on the formation of Zn I lines in the atmospheres of RGB stars is minor. The value of the 3D–1D LTE abundance correction for Zn I line at 481.05 nm is indeed very similar to that obtained using the same CO5BOLD and LHD models by Duffau et al. (2017).
3D–1D abundance corrections, Δ3D−1D LTE, computed for different strengths of Zn I lines used in this work (see text for details).
We note that the computed abundance corrections were not used to obtain the 3D-corrected abundances by adding the 3D–1D LTE corrections to the determined 1D LTE abundances of zinc. The main reason for this was that the obtained Δ3D−1D LTE corrections were small. Therefore, if applied, they result in a very small and nearly uniform shift of abundances determined in all RGB stars, with no effect on the intrinsic abundance spread of zinc and/or various possible relations in the [Zn/Fe]–[Na/Fe] and [Zn/Fe]–[Mg/Fe] planes.
3. Results and discussion
The mean value of the zinc-to-iron abundance ratio we obtained in the sample of 27 RGB stars is 〈[Zn/Fe]〉1D LTE = 0.11 ± 0.09, where the number after the ± sign is the RMS abundance variation due to the star-to-star scatter. Since this variation is not insignificant, the question arises whether it might be a result of intrinsic scatter in the zinc abundance. To answer this question, we applied the maximum-likelihood (ML) technique to evaluate the mean zinc-to-iron abundance ratio, 〈[Zn/Fe]〉, as well as its intrinsic spread, , in our sample of RGB stars. For this, we followed the procedure used in Paper I, which was based on the prescription of Mucciarelli et al. (2012, 2015). The average zinc-to-iron ratio obtained in the ML test is 〈[Zn/Fe]〉 = 0.11 and its uncertainty is ±0.03. We note that the latter value is the error of the mean and therefore is smaller than the RMS variation due to the star-to-star abundance scatter. The determined intrinsic abundance scatter is
, thus there is no intrinsic variation in the zinc abundance in the RGB stars of 47 Tuc. We stress, however, that this non-detection may be at least partly due to significant uncertainties of the zinc abundance measurements in individual stars, which in turn were caused by the relatively low S/N in the observed spectra.
To investigate possible relations between the abundance of zinc and those of other light elements, we used the abundances of Na, Mg, and K determined in our sample RGB stars in Paper I (all abundances of light elements were determined using the 1D NLTE methodology). It is well established that the sodium abundance shows a significant intrinsic star-to-star variation and is anticorrelated with that of oxygen in the vast majority of GGCs studied so far. Therefore, the relation between abundances of Na and Zn would indicate that the two elements were synthesized by the same polluters. In some GGCs, the Mg– Al correlation is observed and points to the Mg–Al cycle, which occurs at higher temperatures than the synthesis of Na. Thus, in a similar way, the relation between Zn and Mg abundances may indicate that the nucleosynthesis of the two elements is connected. The case with potassium is less clear. Recently, Mucciarelli et al. (2017) claimed a detection of a K–Na correlation and K–O anticorrelation in a sample of 144 RGB stars in 47 Tuc. However, our own study of the K abundance in RGB and TO stars in 47 Tuc (32 and 75 objects, respectively) does not support this claim (Paper I and Černiauskas et al. 2018; hereafter Paper II; we note that our data may suggest a weak K–Na anticorrelation in TO stars, but its statistical significance is low). It is worthwhile noting that, if confirmed, relations between the abundance of K and other light elements would be very difficult to explain from the theoretical point of view because K is synthesized at temperatures ≈2 × 108 K, at which all sodium should be destroyed (e.g., Prantzos et al. 2017). Despite this, we wished to verify whether a relation between the abundance of zinc and potassium in the sample of our RGB stars existed, which in case of a positive detection would add another piece of information that would need to be explained by the chemical evolution models of the GGCs. The yields from AGBs depend on the stellar mass, thus in order to reproduce general trends of the light chemical elements (O, Na, Mg, and Al) seen in different globular clusters, the mass range of AGB stars should be in the range of 5–9 M⊙ (D’Ercole et al. 2010; Ventura et al. 2012). Zinc may be synthesized in slightly less massive AGB stars, 3–6 M⊙ (Karakas et al. 2009), thus there may be a slight overlap in stellar masses in which both zinc and light elements could synthesized in AGB stars. In such a situation, relations between the abundances of zinc and those of light elements would be expected.
The determined [Zn/Fe] ratios are plotted versus various light element element-to-iron ratios in Fig. 4. Similar to Paper I, we used Student’s t-test to verify the validity of the null hypothesis, meaning that there is no correlation in the different panels of Fig. 4. For this, using each x–y dataset shown in the three panels of Fig. 4, we computed the two-tailed probability, p, that the t-value in the given dataset could attain the value we obtained when there is no correlation in the given x–y plane. In all panels of Fig. 4, Pearson’s correlation coefficients were computed by taking errors on both x and y axes into account (they are marked as p in Fig. 4). In addition, we also carried out non-parametric Spearman’s and Kendall’s τ rank-correlation tests, in this case without accounting for errors in the determined abundances. Finally, we also computed p-values using Pearson’s correlation coefficients obtained without taking abundance errors into account (marked as pno_errors in Fig. 4). These tests were also performed using the data shown in Figs. 5–7 (see discussion below). The p-values determined in all tests are provided in the corresponding panels of Figs. 4–7.
![]() |
Fig. 4. [Zn/Fe] abundance ratios determine in the sample of 27 RGB stars in 47 Tuc and plotted versus [Na/Fe], [Mg/Fe], and [K/Fe] ratios obtained in the same stars in Paper I. The color on the left and right sides of the symbols denote quality class of the spectral lines used to determined abundance ratios on the y and x axes, respectively (see Sect. 2.2). We also list the values of two-tailed probabilities, p, computed using Pearson’s parametric correlation coefficients (with and without abundance errors, p and pno_errors, respectively), as well as Spearman’s (pSpearman) and Kendall’s (pKendall) non-parametric rank-order correlation coefficients (see Sect. 3). |
![]() |
Fig. 5. [Zn/Fe]1D LTE abundance ratios determined in the sample of 27 RGB stars in 47 Tuc and plotted versus the projected radial distance of individual stars (small filled symbols). The symbol colors are identical to those shown in Fig. 4. Large open circles are averages obtained in non-overlapping Δr/rh = 1 bins. |
The results of these tests do not show any indications of statistically significant relations between the abundances of zinc and the other light elements shown in Fig. 4: in all cases, the p-values are ≥0.64. Similarly, there is no statistically significant relation between the determined [Zn/Fe] ratios and the normalized distance from the cluster center, r/rh (Fig. 5): all p-values are 0.60 (here, r is the projected distance from the cluster center and rh is the half-light radius of 47 Tuc taken from Trager et al. 1993, rh = 174″). Neither do we find a statistically significant relation between the [Zn/Fe] ratios and absolute radial velocities of RGB stars, |Δυr| ≡ |υrad − 〈υrad〉clust|, where υrad is radial velocity of the individual star and 〈υrad〉clust = −18.6km s−1 is the mean radial velocity of the sample (Fig. 6). In this test, the absolute radial velocities were taken from Paper I, while the obtained p-values are all ≥0.49. Finally, we tested whether there may be a relation between the radial velocity dispersion of individual stars (computed in non-overlapping 0.1 dex wide bins of [Zn/Fe] ratios) and their average [Zn/Fe] values (Fig. 7). Again, we find no statistically significant relation, with p-values for all tests ≥0.33.
![]() |
Fig. 6. Absolute radial velocities determined for the sample of 27 RGB stars in 47 Tuc and plotted versus [Zn/Fe]1D LTE ratios of individual stars (small filled symbols). The symbol colors are identical to those shown in Fig. 4. Large open circles are averages obtained in non-overlapping Δ[Zn/Fe]1D LTE = 0.1 bins. |
![]() |
Fig. 7. Radial velocity dispersion of RGB stars in 47 Tuc plotted versus the [Zn/Fe]1D LTE ratio (both quantities were computed in nonoverlapping Δ[Zn/Fe]1D LTE = 0.1 bins). The black solid line shows a linear fit to the data. |
We also performed all statistical tests mentioned above using the three sub-samples of zinc abundances sorted according to the quality class of Zn I lines from which the abundances were determined. For none of the three subsamples (A–C, marked with different symbols in Figs. 4–6) did we find statistically significant relations in any of the planes shown in Figs. 4–6.
Our results therefore indicate that there is no relation between the nucleosynthesis of zinc and that of light elements in 47 Tuc. The sample-averaged zinc-to-iron abundance ratio determined in this study, 〈[Zn/Fe]〉1D LTE = 0.11 ± 0.09, agrees to within less than one sigma to those obtained by Thygesen et al. (2014), [Zn/Fe]1D LTE = 0.26 ± 0.13 (13 RGB stars), and Duffau et al. (2017), 〈[Zn/Fe]〉1D LTE = 0.17 ± 0.10 (19 RGB stars). Our value is nearly identical to the average [Zn/Fe] ratio typical to Galactic field stars at the metallicity of 47 Tuc (e.g., Mishenina et al. 2002; Barbuy et al. 2015; Duffau etal. 2017; da Silveira et al. 2018). It also agrees well with the [Zn/Fe] ratio determined at this galactocentric distance for the Galactic red giants of similar metallicity (Duffau et al. 2017; Ernandes et al. 2018). All this evidence suggests that the production of zinc in 47 Tuc most likely followed the same pathways as in Galactic field stars of the same metallicity, and it most likely occurred through α-element nucleosynthesis. We find no evidence that zinc could have been synthesized in 47 Tuc during the s-process nucleosynthesis.
4. Summary and conclusions
We determined the abundance of zinc in the atmospheres of 27 RGB stars in the Galactic globular cluster 47 Tuc. The abundances were obtained using archival 2dF/HERMES spectra (471.5–490.0 nm, R = 28 000, S/N ≈ 50) that were obtained with the Anglo-Australian Telescope. Spectroscopic data were analyzed using 1D ATLAS9 model atmospheres and the 1D LTE abundance analysis methodology. The 1D LTE spectral line synthesis was performed with the SYNTHE package. The obtained sample-averaged zinc-to-iron abundance ratio is 〈[Zn/Fe]〉1D LTE = 0.11 ± 0.09, here the value following the ± sign is RMS abundance variation due to star-to-star scatter. We also used 3D hydrodynamical CO5 BOLD and 1D hydrostatic LHD model atmospheres to compute the 3D–1D LTE abundance corrections for Zn I lines. The obtained 3D–1D LTE abundance corrections are always positive and are in the range 0.04–0.05 dex, indicating that the influence of convection on the formation of Zn I lines in the atmospheres of RGB stars in 47 Tuc is minor. Applying these corrections would lead to a 3D-corrected estimate of the average zinc-to-iron abundance ratio of 〈[Zn/Fe]〉1D LTE = 0.16 ± 0.09.
The averaged 〈[Zn/Fe]〉1D LTE = 0.11 ± 0.09 ratio derived in this study agrees well with the mean [Zn/Fe] ratios in 47 Tuc determined in the previous studies that were based on smaller samples of stars, [Zn/Fe]1D LTE = 0.26 ± 0.13 (Thygesen et al. 2014, 13 RGB stars) and 〈[Zn/Fe]〉1D LTE = 0.17 ± 0.10, Duffau et al. (2017, 19 RGB stars). Our results show no statistically significant relations between the [Zn/Fe]1D LTE ratio and abundances of light elements (Na, Mg, K), as well as between the [Zn/Fe]1D LTE ratio and kinematic properties of the RGB stars. The sampleaveraged [Zn/Fe]1D LTE ratio coincides with the mean [Zn/Fe] ratio obtained in Galactic field stars at the metallicity of 47 Tuc. All these facts indicate that the nucleosynthesis of zinc and that of light elements, such as Na and Mg, has proceeded separately in 47 Tuc. The obtained average [Zn/Fe]1D LTE ratio suggests that in 47 Tuc zinc is slightly enhanced at this metallicity and most likely originated during the α-element nucleosynthesis.
The raw spectra were downloaded from the AAT data archive http://site.aao.gov.au/arc-bin/wdb/aat_database/observation_log/make
Acknowledgments
Based on data acquired through the Australian Astronomical Observatory, under program 2013B/13. This work was supported by grants from the Research Council of Lithuania (MIP-089/2015, TAP LZ 06/2013). H.G.L. acknowledgesfinancialsupport by the Sonderforschungsbereich SFB 881 “The Milky Way System” (subprojects A4) of the German Research Foundation (DFG).
References
- Andreasen, D.T., Sousa, S. G., Delgado Mena, E. et al. 2016, A&A, 585, A143 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Asplund, M., Grevesse, N., Sauval, A. J., et al. 2009, ARA&A, 47, 481 [NASA ADS] [CrossRef] [Google Scholar]
- Barbuy, B., Friaça, A. C. S., da Silveira, C. R., et al. 2015, A&A, 580, A40 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Bergbusch, P. A., & Stetson, P. B. 2009, AJ, 138, 1455 [NASA ADS] [CrossRef] [Google Scholar]
- Caffau, E., Ludwig, H.-G., Steffen, M., et al. 2008, A&A, 488, 1031 [CrossRef] [EDP Sciences] [Google Scholar]
- Carretta, E., Bragaglia, A., Gratton, R., et al. 2009, A&A, 505, 117 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Castelli, F. 2005, Mem. Soc. Astron. It., 8, 44 [Google Scholar]
- Cayrel, R., Depagne, E., Spite, M., et al. 2004, A&A, 416, 1117 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Černiauskas, A., Kucčinskas, A., Klevas, J., et al. 2017, A&A, 604, A35 (Paper I) [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Černiauskas, A., Kucčinskas, A., Klevas, J., et al. 2018, A&A, 615, A173 (Paper II) [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Cordero, M. J., Pilachowski, C. A., Johnson, C. I., et al. 2014, ApJ, 780, 94 [NASA ADS] [CrossRef] [Google Scholar]
- da Silveira, C. R., Barbuy, B., Friaça, A. C. S., et al. 2018, A&A, 614, A149 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- D’Ercole, A., D’Antona, F., Ventura, P., Vesperini, E., & McMillan, S. L. W. 2010, MNRAS, 407, 854 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
- De Silva, G. M., Freeman, K. C., Bland-Hawthorn, J., et al. 2015, MNRAS, 449, 2604 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
- Dobrovolskas, V., Kucčinskas, A., Bonifacio, P., et al. 2014, A&A, 565, A121 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- D’Orazi, V., Lucatello, S., Gratton, R., Bragaglia, A., Carretta, E., et al. 2010, ApJ, 713L, 1 [Google Scholar]
- Duffau, S., Caffau, E., Sbordone, L., et al. 2017, A&A, 604, A128 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Ernandes, H., Barbuy, B., Alves-Brito, A., et al. 2018, A&A, 616, A18 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Freytag, B., Steffen, M., Ludwig, H.-G., et al. 2012, J. Comput. Phys., 231, 919 [Google Scholar]
- Karakas, A. I., van Raai, M. A., Lugaro, M., Sterling, N. C., & Dinerstein, H. L. 2009, ApJ, 690, 1130 [NASA ADS] [CrossRef] [Google Scholar]
- Kobayashi, C., Umeda, H., Nomoto, K., et al. 2006, ApJ, 653, 1145 [NASA ADS] [CrossRef] [Google Scholar]
- Kucčinskas, A., Dobrovolskas, V., & Bonifacio, P. 2014, A&A, 568, L4 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Kupka, F., Dubernet, M.-L., & VAMDC Collaboration 2011, Balt. Astron., 20, 503 [Google Scholar]
- Kurucz, R. L. 1993, ATLAS9 Stellar Atmosphere Programs and 2 km s−1 Grid, CD-ROM No. 13 (Cambridge, Mass.) [Google Scholar]
- Kurucz, R. L. 2006 ArXiv e-prints [arXiv:astro-ph/0605029] [Google Scholar]
- Lodders, K., Palme, H., & Gail, H.-P. 2009, in Abundances of the Elements in the Solar System, ed. J. E. Trümper, Landolt Börnstein, New Series (Berlin, Heidelberg, New York: Springer), VI/4B, 560 [Google Scholar]
- Matteucci, F., Raiteri, C. M., Busso, M., et al. 1993, A&A, 272, 421 [NASA ADS] [Google Scholar]
- Mishenina, T. V., Kovtyukh, V. V., Soubiran, C., Travaglio, C., & Busso, M. 2002, A&A, 396, 189 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Mucciarelli, A., Bellazzini, M., Ibata, R., et al. 2012, MNRAS, 426, 2889 [NASA ADS] [CrossRef] [Google Scholar]
- Mucciarelli, A., Bellazzini, M., Merle, T., et al. 2015, ApJ, 801, 68 [NASA ADS] [CrossRef] [Google Scholar]
- Mucciarelli, A., Merle, T., & Bellazzini, M. 2017, A&A, 600, A104 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Nissen, P. E., & Schuster, W. J. 2011, A&A, 530, A15 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Nissen, P. E., Chen, Y. Q., Asplund, M., & Pettini, M. 2004, A&A, 415, 993 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Piskunov, N. E., Kupka, F., Ryabchikova, T. A., Weiss, W. W., & Jeffery, C. S. 1995, A&AS, 112, 525 [NASA ADS] [Google Scholar]
- Prantzos, N., Charbonnel, C., & Iliadis, C. 2017, A&A, 608, A28 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Ramírez, I., & Meléndez, J. 2005, ApJ, 626, 465 [NASA ADS] [CrossRef] [Google Scholar]
- Sbordone, L. 2005, Mem. Soc. Astron. It., 8, 61 [Google Scholar]
- Sbordone, L., Bonifacio, P., Castelli, F., & Kurucz, R. L. 2004, Mem. Soc. Astron. It., 5, 93 [Google Scholar]
- Sneden, C., Gratton, R. G., & Crocker, D. A. 1991, A&A, 246, 354 [NASA ADS] [Google Scholar]
- Steffen, M., Caffau, E., & Ludwig, H.-G. 2013, Mem. Soc. Astron. It. Suppl., 24, 37 [Google Scholar]
- Takeda, Y., Hashimoto, O., Taguchi, H., et al. 2005, PASJ, 57, 751 [NASA ADS] [Google Scholar]
- Thygesen, A. O., Sbordone, L., Andrievsky, S., et al. 2014, A&A, 572, A108 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Trager, S.C., Djorgovski, S.G., & King, I.R. 1993, ASP Conf. Ser., 50, 347 [Google Scholar]
- Ventura, P., D’Antona, F., Mazzitelli, I., & Gratton, R. 2001, ApJ, 550, L65 [NASA ADS] [CrossRef] [Google Scholar]
- Ventura, P., D’Antona, F., Di Criscienzo, M., et al. 2012, ApJ, 761, L30 [NASA ADS] [CrossRef] [Google Scholar]
Appendix A: Zinc abundance in the Sun
We determined the solar zinc abundance using the same Zn I lines and their atomic parameters (Table 1) as used by us in the analysis of the zinc abundance in RGB stars in 47 Tuc. Synthetic Zn I line profiles were computed using the SYNTHE spectral synthesis package and 1D hydrostatic ATLAS9 model atmospheres (see Sect. 2.2). Atomic line parameters were taken from the VALD-3 atomic database (Piskunov et al. 1995; Kupka et al. 2011).
To determine the solar zinc abundance, we used the rereduced Kitt Peak Solar Flux atlas of Kurucz (2006). The solar spectrum covers a wavelength range of 300–1000 nm, with R = 523 000, and S/N ≈ 3500 in the wavelength range around Zn I lines. The solar 1D LTE zinc abundance was determined using the spectral lines located at 472.2 nm and 481.0 nm, with the microturbulence velocity set to 1.0 km s−1. Each line was fit independently with the synthetic spectrum (Fig. A.1) that was computed using an ATLAS9 model atmosphere of the Sun (Teff = 5777 K and log g = 4.43, as recommended by Andreasen et al. 2015). The solar abundances of zinc obtained from the 472.2 nm and 481.0 nm Zn I lines were and
, respectively, with a mean value of
(the value after the “±” sign is the RMS variation). This agrees well with those determined in the earlier studies, for example,
obtained using the same Zn I lines by Takeda et al. (2005),
derived by Lodders et al. (2009), and
recommended by Asplund et al. (2009).
![]() |
Fig. A.1. Zn I lines (filled gray circles) observed at 472.2 nm and 481.0 nm in the solar spectrum. The synthetic solar spectrum is shown as a green solid line. Numbers above the spectrum correspond to the atomic numbers of the line-producing chemical element. |
We further estimated the uncertainty of the solar zinc abundance determination using the same procedure and uncertainties in the atmospheric parameters of the Sun and spectral line fitting procedure as described in Paper II. Errors from the individual uncertainties are provided in Table A.1. The final value of the solar zinc abundance obtained from the two Zn I lines is .
Total error in the abundance of zinc determined from the Zn I lines in the atmosphere of the Sun, σ(A)tot.
Errors in the abundance of zinc determined from the two Zn I lines used in this study and different spectral line quality (A–C, see Sect. 3).
Appendix B: Uncertainties in zinc abundances determined in the RGB stars in 47 Tuc
The uncertainty in the determined zinc abundances was estimated following the same procedure as in Paper I. For this, we assumed that the error in the determined abundances occurs due to an inaccurate determination of the atmospheric parameters (Teff, log g, ξt) and the spectral line profile fitting, with the additional assumption that the latter depends on the goodness of line profile fit and the placement of continuum level. These uncertainties were estimated in the following way:
-
Errors in the determination of atmospheric parameters: Following Paper I, we assumed that the error in the determination of the effective temperature of the RGB stars was ±65 K. The errors in the determined zinc abundance resulting from this uncertainty, σ(Teff), are provided in Table A.2, Col. 4, for both Zn I lines. For the uncertainty in surface gravity, we adopted a conservative (and, in our view, more realistic) value ±0.1 dex; the resulting errors in the zinc abundances, σ(log g), are listed in Table A.2, Col. 5. Since we did not derive individual values of the microturbulence velocity for individual RGB stars, we adopted an uncertainty in the microturbulence velocity of ±0.15 km s−1. This value is based on the RMS variation of microturbulence velocity determined in RGB stars in 47 Tuc by Carretta et al. (2009, 58 objects) and Cordero et al. (2014, 81 objects), with stars occupying the same range in the effective temperature as those used in our sample. The obtained abundance uncertainties that occur due to errors in microturbulence velocity, σ(ξt), are listed in Table A.2, Col. 6.
-
Errors in the spectral line profile fitting: The error in continuum determination was estimated in the same way as in Paper I, that is, by measuring the dispersion at the continuum level (inverse signal-to-noise ratio) in the spectral regions deemed to be free of spectral lines (see Appendix B in Paper I for details). The errors in the line profile fitting were computed individually for each line quality class (see Sect. 2.2) using RMS values as a measure of differences between the observed and best-fit synthetic line profiles. These RMS values were used to compute the change in the determined abundance of zinc. The final errors due to continuum determination and spectral line profile fitting are shown in Table A.2, Cols. 7 and 8, respectively.
The errors in the determined zinc abundances resulting from various uncertainty sources are listed in Table A.2. The total uncertainties in zinc abundance (obtained by adding individual components in quadratures, Col. 9 in Table A.2) are shown in Figs. 3–5. These errors were also used in the ML analysis aimed to determine the average value of the zinc abundance and its intrinsic spread in the sample of RGB stars (see Sect. 3).
We stress that in this procedure, errors resulting from uncertainties in the atomic line parameters were ignored. Therefore, the total abundance uncertainties listed in Table A.2 are only lower limit estimates because they do not account for the various systematic uncertainties that are unavoidable in the abundance analysis procedure.
All Tables
3D–1D abundance corrections, Δ3D−1D LTE, computed for different strengths of Zn I lines used in this work (see text for details).
Total error in the abundance of zinc determined from the Zn I lines in the atmosphere of the Sun, σ(A)tot.
Errors in the abundance of zinc determined from the two Zn I lines used in this study and different spectral line quality (A–C, see Sect. 3).
All Figures
![]() |
Fig. 1. Upper part of the color-magnitude diagram of 47 Tuc constructed using VI photometry from Bergbusch & Stetson (2009). RGB stars used in this study are marked as filled green circles. |
In the text |
![]() |
Fig. 2. Typical fits of synthetic spectra (red solid lines) to the observed Zn I line profiles (marked by vertical ticks) in the 2dF/HERMES spectra (filled gray circles) of the target RGB stars S1800 (top row; Teff = 4510 K, log g = 1.90) and S167 (bottom row; Teff = 4390 K, log g = 1.60). We also indicate zinc abundances A(Zn) determined from each spectral line, together with their errors (Sect. 2.2). The quality class of the Zn I line is indicated at the top of each panel (see Sect. 2.2 for the line class definition). |
In the text |
![]() |
Fig. 3. [Zn/Fe] abundance ratios determined in the sample of 27 RGB stars in 47 Tuc and plotted versus the effective temperature of individual stars. |
In the text |
![]() |
Fig. 4. [Zn/Fe] abundance ratios determine in the sample of 27 RGB stars in 47 Tuc and plotted versus [Na/Fe], [Mg/Fe], and [K/Fe] ratios obtained in the same stars in Paper I. The color on the left and right sides of the symbols denote quality class of the spectral lines used to determined abundance ratios on the y and x axes, respectively (see Sect. 2.2). We also list the values of two-tailed probabilities, p, computed using Pearson’s parametric correlation coefficients (with and without abundance errors, p and pno_errors, respectively), as well as Spearman’s (pSpearman) and Kendall’s (pKendall) non-parametric rank-order correlation coefficients (see Sect. 3). |
In the text |
![]() |
Fig. 5. [Zn/Fe]1D LTE abundance ratios determined in the sample of 27 RGB stars in 47 Tuc and plotted versus the projected radial distance of individual stars (small filled symbols). The symbol colors are identical to those shown in Fig. 4. Large open circles are averages obtained in non-overlapping Δr/rh = 1 bins. |
In the text |
![]() |
Fig. 6. Absolute radial velocities determined for the sample of 27 RGB stars in 47 Tuc and plotted versus [Zn/Fe]1D LTE ratios of individual stars (small filled symbols). The symbol colors are identical to those shown in Fig. 4. Large open circles are averages obtained in non-overlapping Δ[Zn/Fe]1D LTE = 0.1 bins. |
In the text |
![]() |
Fig. 7. Radial velocity dispersion of RGB stars in 47 Tuc plotted versus the [Zn/Fe]1D LTE ratio (both quantities were computed in nonoverlapping Δ[Zn/Fe]1D LTE = 0.1 bins). The black solid line shows a linear fit to the data. |
In the text |
![]() |
Fig. A.1. Zn I lines (filled gray circles) observed at 472.2 nm and 481.0 nm in the solar spectrum. The synthetic solar spectrum is shown as a green solid line. Numbers above the spectrum correspond to the atomic numbers of the line-producing chemical element. |
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.