The s-process enriched star HD 55496: origin from a globular cluster or from the tidal disruption of a dwarf galaxy?
C. B. Pereira, N. A. Drake, F. Roig

TL;DR
This study analyzes the chemical composition and orbital dynamics of HD 55496, a metal-poor s-process enriched star, suggesting it originated either from a globular cluster or from the tidal disruption of a dwarf galaxy, with evidence leaning towards the latter.
Contribution
It provides a detailed chemical and dynamical analysis of HD 55496, proposing a novel origin scenario involving dwarf galaxy disruption based on its peculiar abundance pattern and orbital properties.
Findings
HD 55496 has [Fe/H] = -1.55 and is s-process enriched.
It exhibits a Na-O anti-correlation and is aluminum rich.
Dynamical analysis suggests a low probability of globular cluster origin, supporting a dwarf galaxy origin.
Abstract
We present a new abundance analysis of HD 55496, previously known as a metal-poor barium star. We found that HD 55496 has a metallicity [Fe/H] and is s-process enriched. We find that HD 55496 presents four chemical peculiarities: (i) a Na-O abundance anti-correlation; (ii) it is aluminum rich; (iii) it is carbon poor for a s-process enriched star and (iv) the heavy 2 s-process peak elements, such as Ba, La, Ce, and Nd, present smaller abundances than the lighter s-process elements, such as Sr, Y and Zr, which is not usually observed among the chemically-peculiar binary stars at this metallicity. The heavy-element abundance pattern suggests that the main source of the neutrons is the Ne(,n)Mg reaction. Taken all these abundance evidence together into consideration, this strongly suggests that HD 55496 is a "second generation of globular cluster star"…
| Parameter | Value | Ref. |
| (K) | 4 70070 | 1 |
| 4 850 | 2 | |
| 4 800 | 3 | |
| 4 750 | 4 | |
| (dex) | 1.90.2 | 1 |
| 2.05 | 2 | |
| 2.8 | 3 | |
| 2.5 | 4 | |
| [Fe/H] (dex) | 1.550.11 | 1 |
| 1.49 | 2 | |
| 1.55 | 3 | |
| (km s-1) | 1.50.3 | 1 |
| 1.5 | 2 | |
| 2.7 | 3 | |
| 1.5 | 4 |
| Species | ⊙ | n(#) | [X/H] | [X/Fe] | ||
|---|---|---|---|---|---|---|
| Li | 3.31 | syn | 0.00 | 0.15 | — | — |
| C (CH) | 8.52 | syn | 6.62 | 0.05 | 1.90 | 0.35 |
| N (CN) | 7.92 | syn | 8.17 | 0.05 | 0.25 | 1.80 |
| O | 8.83 | syn | 7.43 | 0.10 | 1.40 | 0.15 |
| Na i | 6.33 | 8 | 5.59 | 0.13 | 0.74 | 0.81 |
| Mg i | 7.58 | 5 | 6.73 | 0.14 | 0.85 | 0.70 |
| Al i | 6.47 | 6 | 5.66 | 0.12 | 0.81 | 0.74 |
| Si i | 7.55 | 6 | 6.66 | 0.09 | 0.89 | 0.66 |
| Ca i | 6.36 | 16 | 5.39 | 0.12 | 0.97 | 0.58 |
| Ti i | 5.02 | 24 | 3.72 | 0.07 | 1.30 | 0.25 |
| Fe i | 7.50 | 70 | 5.95 | 0.11 | 1.55 | — |
| Fe ii | 7.50 | 10 | 5.92 | 0.09 | 1.58 | — |
| Cr i | 5.67 | 10 | 4.06 | 0.11 | 1.61 | 0.06 |
| Ni i | 6.25 | 35 | 4.70 | 0.09 | 1.55 | 0.00 |
| Sr i | 2.97 | 1 | 2.04 | 0.07 | 0.93 | 0.62 |
| Y ii | 2.24 | 7 | 1.53 | 0.13 | 0.71 | 0.84 |
| Zr i | 2.60 | 6 | 2.20 | 0.06 | 0.40 | 1.15 |
| Zr ii | 2.60 | 4 | 2.28 | 0.06 | 0.32 | 1.23 |
| Mo i | 1.92 | 1 | 1.47 | 0.10 | 0.45 | 1.10 |
| Ba ii | 2.13 | syn | 1.20 | 0.06 | 0.93 | 0.62 |
| La ii | 1.17 | syn | 0.07 | 0.08 | 1.24 | 0.31 |
| Ce ii | 1.58 | 10 | 0.54 | 0.10 | 1.04 | 0.51 |
| Nd ii | 1.50 | 21 | 0.34 | 0.10 | 1.16 | 0.39 |
| Sm ii | 1.01 | 6 | 0.35 | 0.06 | 1.36 | 0.19 |
| Pb i | 1.95 | syn | 1.75 | 0.05 | 0.20 | 1.35 |
| Species | [Fe/H] | |||||||
|---|---|---|---|---|---|---|---|---|
| 70 K | 0.2 | 0.3 km s-1 | 0.1 | 3 mÅ | ||||
| Li | 0.15 | 0.10/0.10 | 0.05/0.05 | 0.00/0.00 | 0.00/0.00 | — | 0.19/0.19 | 0.15 |
| Fe i | 0.01 | 0.09/0.09 | 0.01/0.01 | 0.11/0.13 | 0.01/0.01 | 0.07/0.07 | 0.16/0.17 | 0.11 |
| Fe ii | 0.03 | 0.02/0.04 | 0.09/0.07 | 0.03/0.05 | 0.02/0.01 | 0.09/0.09 | 0.10/0.12 | 0.09 |
| Na i | 0.05 | 0.06/0.06 | 0.02/0.02 | 0.05/0.05 | 0.01/0.01 | 0.05/0.06 | 0.09/0.11 | 0.13 |
| Mg i | 0.06 | 0.07/0.06 | 0.04/0.04 | 0.07/0.07 | 0.00/0.01 | 0.05/0.04 | 0.11/0.12 | 0.14 |
| Al i | 0.05 | 0.04/0.05 | 0.01/0.01 | 0.02/0.02 | 0.01/0.01 | 0.06/0.07 | 0.09/0.10 | 0.12 |
| Si i | 0.04 | 0.02/0.01 | 0.02/0.01 | 0.01/0.02 | 0.00/0.00 | 0.09/0.10 | 0.10/0.11 | 0.09 |
| Ca i | 0.03 | 0.07/0.07 | 0.03/0.02 | 0.12/0.12 | 0.02/0.01 | 0.06/0.07 | 0.16/0.16 | 0.12 |
| Ti i | 0.01 | 0.12/0.13 | 0.02/0.02 | 0.06/0.07 | 0.01/0.02 | 0.07/0.07 | 0.15/0.17 | 0.07 |
| Cr i | 0.03 | 0.10/0.10 | 0.02/0.01 | 0.08/0.10 | 0.01/0.02 | 0.08/0.09 | 0.13/0.17 | 0.11 |
| Ni i | 0.01 | 0.08/0.09 | 0.02/0.00 | 0.05/0.06 | 0.01/0.00 | 0.07/0.08 | 0.12/0.13 | 0.09 |
| Sr ii | 0.07 | 0.12/0.13 | 0.03/0.03 | 0.17/0.24 | 0.02/0.02 | 0.08/0.08 | 0.23/0.30 | 0.07 |
| Y ii | 0.05 | 0.02/0.01 | 0.06/0.06 | 0.17/0.22 | 0.02/0.01 | 0.08/0.07 | 0.21/0.24 | 0.13 |
| Zr i | 0.02 | 0.13/0.13 | 0.01/0.02 | 0.01/0.02 | 0.00/0.02 | 0.10/0.11 | 0.16/0.17 | 0.06 |
| Zr ii | 0.03 | 0.01/0.00 | 0.07/0.06 | 0.13/0.19 | 0.02/0.01 | 0.08/0.07 | 0.17/0.21 | 0.06 |
| Mo i | 0.10 | 0.11/0.12 | 0.02/0.01 | 0.01/0.01 | 0.01/0.01 | 0.13/0.18 | 0.20/0.24 | 0.10 |
| Ba ii | 0.03 | 0.01/0.00 | 0.10/0.10 | 0.20/0.30 | 0.00/0.00 | — | 0.23/0.31 | 0.06 |
| La ii | 0.04 | 0.05/0.05 | 0.10/0.10 | 0.00/0.00 | 0.00/0.00 | — | 0.12/0.12 | 0.08 |
| Ce ii | 0.03 | 0.02/0.02 | 0.08/0.07 | 0.05/0.09 | 0.02/0.02 | 0.08/0.09 | 0.13/0.15 | 0.10 |
| Nd ii | 0.02 | 0.02/0.03 | 0.07/0.08 | 0.04/0.05 | 0.02/0.03 | 0.08/0.09 | 0.12/0.14 | 0.10 |
| Sm ii | 0.02 | 0.03/0.02 | 0.08/0.07 | 0.02/0.04 | 0.03/0.02 | 0.09/0.09 | 0.13/0.13 | 0.06 |
| Pb i | 0.05 | 0.10/0.20 | 0.00/0.00 | 0.10/0.10 | 0.00/0.00 | — | 0.11/0.22 | 0.05 |
| Species | ||||||||
|---|---|---|---|---|---|---|---|---|
| 70 K | 0.2 | 0.3 km s-1 | 0.20 | 0.20 | 0.20 | |||
| C | 0.05 | 0.10/0.10 | 0.00/0.05 | 0.00/0.00 | — | 0.00 | 0.15 | 0.19/0.19 |
| N | 0.02 | 0.10/0.10 | 0.05/0.05 | 0.00/0.00 | 0.10 | — | 0.20 | 0.29/0.29 |
| O | 0.10(?) | 0.05/0.05 | 0.10/0.10 | 0.00/0.05 | 0.00 | 0.00 | — | 0.15/0.16 |
| RV (km s-1) | Ref. |
|---|---|
| 315.630.23 | 1 |
| 316.210.16 | 2 |
| 315.280.80 | 3 |
| 322 | 4 |
| 315.30317.72 | 5 |
| 31710 | 6 |
| 317 | 7 |
| Name | [Fe/H] | |||
|---|---|---|---|---|
| (%) | (%) | km s-1 | ||
| NGC 5139 | 6.0 | 16.2 | -1.445 | |
| NGC 6402 | 4.0 | 14.8 | -1.145 | |
| NGC 6749 | 4.6 | 16.3 | -1.600 | |
| Terzan 4 | 3.4 | 15.3 | -1.045 |
| Element | (Å) | (eV) | Wλ(mÅ) | |
|---|---|---|---|---|
| Fe i | 5133.69 | 4.18 | 0.20 | 89 |
| 5150.84 | 0.99 | 3.00 | 114 | |
| 5151.91 | 1.01 | 3.32 | 106 | |
| 5159.06 | 4.28 | 0.65 | 30 | |
| 5162.27 | 4.18 | 0.07 | 82 | |
| 5194.94 | 1.56 | 2.09 | 123 | |
| 5198.71 | 2.22 | 2.14 | 87 | |
| 5232.94 | 2.94 | 0.08 | 139 | |
| 5242.49 | 3.63 | 0.97 | 58 | |
| 5250.21 | 0.12 | 4.92 | 85 | |
| 5281.79 | 3.04 | 0.83 | 101 | |
| 5302.31 | 3.28 | 0.74 | 90 | |
| 5307.36 | 1.61 | 2.97 | 84 | |
| 5322.04 | 2.28 | 2.84 | 43 | |
| 5339.93 | 3.27 | 0.68 | 95 | |
| 5341.02 | 1.61 | 1.95 | 135 | |
| 5364.87 | 4.45 | 0.23 | 71 | |
| 5367.47 | 4.42 | 0.43 | 75 | |
| 5369.96 | 4.37 | 0.54 | 82 | |
| 5373.71 | 4.47 | 0.71 | 22 | |
| 5393.17 | 3.24 | 0.72 | 91 | |
| 5400.50 | 4.37 | 0.10 | 59 | |
| 5410.91 | 4.47 | 0.40 | 71 | |
| 5445.04 | 4.39 | 0.04 | 62 | |
| 5487.75 | 4.32 | 0.65 | 43 | |
| 5506.78 | 0.99 | 2.80 | 130 | |
| 5522.45 | 4.21 | 1.40 | 14 | |
| 5554.90 | 4.55 | 0.38 | 41 | |
| 5560.21 | 4.43 | 1.04 | 15 | |
| 5563.60 | 4.19 | 0.84 | 46 | |
| 5567.39 | 2.61 | 2.56 | 38 | |
| 5569.62 | 3.42 | 0.49 | 95 | |
| 5572.84 | 3.40 | 0.28 | 104 | |
| 5576.09 | 3.43 | 0.85 | 76 | |
| 5624.02 | 4.39 | 1.33 | 12 | |
| 5686.53 | 4.55 | 0.45 | 26 | |
| 5691.50 | 4.30 | 1.37 | 10 | |
| 5717.83 | 4.28 | 0.97 | 25 | |
| 5731.76 | 4.26 | 1.15 | 22 | |
| 5762.99 | 4.21 | 0.41 | 58 | |
| 5791.02 | 3.21 | 2.27 | 36 | |
| 5883.82 | 3.96 | 1.21 | 34 | |
| 5916.25 | 2.45 | 2.99 | 37 | |
| 6024.06 | 4.55 | 0.06 | 53 | |
| 6027.05 | 4.08 | 1.09 | 29 | |
| 6056.01 | 4.73 | 0.40 | 25 | |
| 6065.48 | 2.61 | 1.53 | 106 | |
| 6079.01 | 4.65 | 0.97 | 25 | |
| 6136.61 | 2.45 | 1.40 | 115 | |
| 6137.69 | 2.59 | 1.40 | 107 | |
| 6151.62 | 2.18 | 3.29 | 31 | |
| 6157.73 | 4.08 | 1.11 | 36 | |
| 6170.51 | 4.79 | 0.38 | 28 | |
| 6173.34 | 2.22 | 2.88 | 59 | |
| 6191.56 | 2.43 | 1.42 | 114 | |
| 6200.31 | 2.60 | 2.44 | 51 | |
| 6213.43 | 2.22 | 2.48 | 74 | |
| 6230.72 | 2.56 | 1.28 | 121 | |
| 6252.56 | 2.40 | 1.72 | 107 | |
| 6265.13 | 2.18 | 2.55 | 85 | |
| 6322.69 | 2.59 | 2.43 | 59 | |
| 6393.60 | 2.43 | 1.43 | 120 | |
| 6411.65 | 3.65 | 0.66 | 82 | |
| 6419.95 | 4.73 | 0.09 | 41 | |
| 6421.35 | 2.28 | 2.01 | 100 | |
| 6430.85 | 2.18 | 2.01 | 108 | |
| 6592.91 | 2.72 | 1.47 | 93 | |
| 6593.87 | 2.44 | 2.42 | 68 |
| (Å) | Element | (eV) | Ref | Wλ(mÅ) | |
| 4664.81 | 1Na i | 2.10 | 1.55 | T94 | 28 |
| 4668.60 | 1Na i | 2.10 | 1.30 | T94 | 29 |
| 4982.81 | 1Na i | 2.10 | 0.95 | T94 | 54 |
| 5682.65 | 1Na i | 2.10 | 0.70 | T94 | 77 |
| 5688.22 | 1Na i | 2.10 | 0.37 | T94 | 97 |
| 6154.22 | 1Na i | 2.10 | 1.51 | T94 | 23 |
| 6160.75 | 1Na i | 2.10 | 1.21 | T94 | 40 |
| 8183.26 | 1Na i | 2.10 | 0.22 | T94 | 168 |
| 4702.99 | 1Mg i | 4.35 | 0.52 | A2004 | 146 |
| 4730.04 | 1Mg i | 4.34 | 2.39 | R03 | 35 |
| 5528.40 | 1Mg i | 4.35 | 0.49 | A2004 | 157 |
| 5711.10 | 1Mg i | 4.34 | 1.68 | R99 | 78 |
| 8736.04 | 1Mg i | 5.94 | 0.34 | WSM1 | 82 |
| 6696.03 | 1Al i | 3.14 | 1.48 | MR94 | 18 |
| 6698.67 | 1Al i | 3.14 | 1.63 | R03 | 17 |
| 7835.32 | 1Al i | 4.02 | 0.58 | R03 | 28 |
| 7836.13 | 1Al i | 4.02 | 0.40 | R03 | 25 |
| 8772.88 | 1Al i | 4.02 | 0.25 | R03 | 41 |
| 8773.91 | 1Al i | 4.02 | 0.07 | R03 | 59 |
| 5793.08 | 1Si i | 4.93 | 2.06 | R03 | 46 |
| 6125.03 | 1Si i | 5.61 | 1.54 | E93 | 15 |
| 6131.58 | 1Si i | 5.62 | 1.69 | C2003 | 10 |
| 6145.08 | 1Si i | 5.62 | 1.43 | E93 | 43 |
| 7800.00 | 1Si i | 6.18 | 0.72 | E93 | 27 |
| 8742.45 | 1Si i | 5.87 | 0.51 | E93 | 65 |
| 5581.80 | 1Ca i | 2.52 | 0.67 | C2003 | 76 |
| 5601.29 | 1Ca i | 2.52 | 0.52 | C2003 | 77 |
| 5857.46 | 1Ca i | 2.93 | 0.11 | C2003 | 99 |
| 6102.73 | 1Ca i | 1.88 | 0.79 | D2002 | 113 |
| 6122.23 | 1Ca i | 1.89 | 0.32 | D2002 | 142 |
| 6166.44 | 1Ca i | 2.52 | 1.14 | R03 | 48 |
| 6169.04 | 1Ca i | 2.52 | 0.80 | R03 | 74 |
| 6169.56 | 1Ca i | 2.53 | 0.48 | DS91 | 88 |
| 6439.08 | 1Ca i | 2.52 | 0.47 | D2002 | 137 |
| 6449.82 | 1Ca i | 2.52 | 0.50 | C2003 | 85 |
| 6455.60 | 1Ca i | 2.51 | 1.29 | R03 | 40 |
| 6464.68 | 1Ca i | 2.52 | 2.41 | C2003 | 13 |
| 6471.66 | 1Ca i | 2.51 | 0.69 | S86 | 78 |
| 6493.79 | 1Ca i | 2.52 | 0.11 | DS91 | 113 |
| 6499.65 | 1Ca i | 2.52 | 0.81 | C2003 | 72 |
| 6717.69 | 1Ca i | 2.71 | 0.52 | C2003 | 84 |
Peer Reviews
No public reviews on file for this paper yet. If you reviewed it on a platform where reviews are public (OpenReview, ICLR, NeurIPS, ICML), you can paste yours below so the community can read it here.
Videos
No videos yet. Explain this paper in a talk, walkthrough, or lecture? Add one.
The s-process enriched star HD
55496: origin from a globular cluster or from the tidal disruption of a dwarf galaxy? ††thanks: Based on the observations made with the 2.2m telescope at the European Southern Observatory (La Silla, Chile) under the agreement between Observatório Nacional (Brazil) and European Southern Observatory (ESO)
C.B. Pereira1, N. A. Drake1,2,3, F. Roig1
1Observatório Nacional/MCTIC, Rua Gen. José Cristino, 77, 20921-400, Rio de Janeiro, Brazil
2 Laboratory of Observational Astrophysics, Saint Petersburg State University, Universitetski pr. 28, Petrodvoretz 198504, Saint Petersburg,
Russia
3 Laboratório Nacional de Astrofísica/MCTIC, Rua Estados Unidos, 154, 37500-000, Itajubá, Brazil E-mail:[email protected]:[email protected]:[email protected]
(Accepted xxx. Received xxx; in original form xxxx)
Abstract
We present a new abundance analysis of HD 55496, previously known as a metal-poor barium star. We found that HD 55496 has a metallicity [Fe/H] = and is s-process enriched. We find that HD 55496 presents four chemical peculiarities : (i) a Na-O abundance anti-correlation; (ii) it is aluminum rich; (iii) it is carbon poor for a s-process enriched star and (iv) the heavy 2nd s-process peak elements, such as Ba, La, Ce, and Nd, present smaller abundances than the lighter s-process elements, such as Sr, Y and Zr, which is not usually observed among the chemically-peculiar binary stars at this metallicity. The heavy-element abundance pattern suggests that the main source of the neutrons is the 22Ne(,n)25Mg reaction. Taken all these abundance evidence together into consideration, this strongly suggests that HD 55496 is a “second generation of globular cluster star” formed from gas already strongly enriched in s-process elements and now is a field halo object. Our dynamical analysis, however, indicates that the past encounter probabilities with the known globular clusters is very small (). This evidence, together with the fact of having a retrograde motion, points to a halo intruder possibly originated from the tidal disruption of a dwarf galaxy.
keywords:
nuclear reactions, nucleosynthesis — stars: abundances — stars: individual: HD 55496 stars: chemically peculiar — stars: evolution — stars: fundamental parameters
††pubyear: 2017††pagerange: The s-process enriched star HD 55496: origin from a globular cluster or from the tidal disruption of a dwarf galaxy? ††thanks: Based on the observations made with the 2.2m telescope at the European Southern Observatory (La Silla, Chile) under the agreement between Observatório Nacional (Brazil) and European Southern Observatory (ESO)–The s-process enriched star HD 55496: origin from a globular cluster or from the tidal disruption of a dwarf galaxy? ††thanks: Based on the observations made with the 2.2m telescope at the European Southern Observatory (La Silla, Chile) under the agreement between Observatório Nacional (Brazil) and European Southern Observatory (ESO)
1 Introduction
It is well known that “second generation of globular cluster stars” present an anti-correlation of the abundances among the light elements, such as a depletion of carbon, oxygen and magnesium accompanied by an overabundance of other light elements, such as sodium, aluminum and nitrogen (e.g. Shetrone 1996, Carretta et al. 2010, just to name a few of the many works dedicated to investigate the chemical abundances in globular clusters). Such abundance pattern or the anti-correlations between these abundances have not been observed in the past among the stars of the halo of low metallicity. Thus, any star found in the field halo in the Galaxy presenting such abundance pattern could be considered as a second generation candidate star of a globular cluster that has escaped from the cluster (Fernández-Trincado et al. 2016, Carretta et al. 2010). Finding these kind of stars within field halo of the Galaxy is not an easy task and few of them have been found in the last years, thanks to large spectroscopic surveys (Ramírez et al. 2012; Schiavon et al. 2017; Fernández-Trincado, 2016, 2017, Martell et al. 2016), and to photometric plus dynamical analysis searching for tidal debris from Cen (Majewski et al. 2012).
In this work we present a new abundance and dynamical analysis of the star HD 55496. This star, first noticed as a barium star by MacConnel et al. (1972) and later classified as a "metal-poor barium star" by Luck & Bond (1991), was observed during the high-resolution spectroscopic survey dedicated to investigate and to determine the heavy-abundance abundance pattern of a sample of barium stars. Barium stars represent among the chemically peculiar stars the largest sample, and therefore they are useful targets to constrain nucleosynthesis process in the asymptotic giant branch (AGB) stars (de Castro et al. 2016, Cseh et al. 2018). Being a barium star, HD 55496 should also be a binary star, but analysis of radial velocity data collected over 14 years were inconclusive (Jorissen et al. 2015).
As we shall see, HD 55496 is an atypical object to be classified as a barium star. Not only due to the inconclusive results of the radial velocity variation, but also due to its peculiar chemical abundance of the elements created by the s-process and the high sodium and aluminum abundances. This chemical peculiarity, turns HD 55496 to be an attractive object for an investigation whether it can be a “second generation of globular cluster star” which is now a field halo object.
This work is organized as follows: In Section 2 we describe the observations. In Section 3 we describe the determination of the atmospheric parameters and the chemical abundances. In Section 4, we discuss the results obtained for the luminosity, for the radial velocity, for the chemical abundances and for the kinematical analysis. Finally, our conclusions are presented in Section 5.
2 Observations
The high-resolution spectrum of HD 55496 was obtained with the FEROS (Fibre-fed Extended Range Optical Spectrograph) echelle spectrograph (Kaufer et al. 1999) at the 2.2 m ESO telescope in La Silla (Chile), during the night of April 3, 2007. The exposure time for this observation was 900 sec. Technical details about the FEROS spectrograph is reported in Santrich et al. (2013).
Since HD 55496 has a high sodium and aluminum abundance it is important to compare the spectrum of HD 55496 with another star whose atmospheric parameters are similar to that of HD 55496 and not either sodium nor aluminum enriched. In this case, the differences in the strengths of absorption lines of sodium and aluminum between the two stars will be only due to a high abundance of sodium and aluminum in HD 55496. Figures 1, 2 and 3 show the spectra of HD 55496 in the spectral regions of the interest of the present work, that is around the Na i and Al i absorption lines in comparison with the spectrum of CD-27∘8864. CD-27∘8864, is a metal-poor star with similar effective temperature ( = 4 650 K), surface gravity ( = 1.6), iron abundance ([Fe/H] = 1.54), and microturbulent velocity 1.6 km s*-1* analyzed in Pereira et al. (2019).
3 Analysis
3.1 Atmospheric Parameters
For the determination of the atmospheric parameters it is first necessary to measure the equivalent widths of the absorption lines of Fe i and Fe ii. The absorption lines selected for this task were selected from Lambert et al. (1996) and that was the same source of lines used in the work of Moriz et al. (2017) in the study of the CD-50∘776 recently identified as a new CEMP-s. 111CEMP stands for Carbon-Enhanced Metal-Poor. The suffix ′′s′′ means that the star is s-process enriched. Table 1 shows the values of the equivalent widths of the Fe i and Fe ii absorption lines with their respective excitation potentials ( in eV) and the values of the transtions. The equivalent widths were determined using the splot task in IRAF by adjusting Gaussian profiles to the observed profiles.
In addition, the determination of the stellar atmospheric parameters, such as effective temperature (), surface gravity (), microturbulence (), and metallicity ([Fe/H]) (we use the notation [X/H]=), were also done in the same way as in Roriz et al. (2017). In brief, it consists in using the local thermodynamic equilibrium (hereafter LTE) model atmospheres of Kurucz (1993) and the spectral analysis code moog, version 2013 (Sneden 1973). Table 2 shows the final adopted atmospheric parameters with their respective errors. Previous atmospheric parameters determination for HD 55496 were obtained by Karinkuzhi & Goswami (2015) and earlier by Luck & Bond (1991). Our results are similar to those of Karinkuzhi & Goswami (2015).
The errors in temperature and microturbulent velocity were estimated considering, respectively, the uncertainty in the value of the slope of the relationships between the abundance of iron, given by the Fe i lines, and the excitation potential and the abundance of iron and the reduced equivalent width. For gravity, the uncertainty was estimated considering that the difference between the Fe i and Fe ii abundances differed by 1 of the standard deviation of the mean value of Fe i abundances. We found for the temperature, surface gravity and for the microturbulent velocity uncertainties of and km s*-1*, respectively.
3.2 Abundance analysis
The abundances of chemical elements Na, Mg, Al, Si, Ca, Ti, Cr, Ni, Sr, Y, Zr, Mo, Ce, Nd and Sm were determined by measuring the equivalent widths of their absorption lines considering local-thermodynamic-equilibrium (LTE) model-atmosphere. We used the line-synthesis code moog, (Sneden, 1973) for the calculations and Table 3 shows the atomic lines used to obtain the abundances of the elements and Table 4 provides the results, the number of lines employed for each species, , and the standard deviations. Our abundances were determined differentially with respect to the solar abundances of Grevesse & Sauval (1998).
The abundances of carbon, nitrogen, oxygen and lithium were determined by means of the spectral synthesis technique. For carbon, the abundance was determined using the CH lines of the system at ÅḞor the determination of the abundance of nitrogen and the 12C/13C isotopic ratio, we used the 12CN lines of the (2, 0) band of the CN red system in the 7994 – 8020 Å wavelength range. For the abundance of oxygen, we used the forbidden line at 6300.3 Å where we adopted = 9.72 from Allende Prieto et al. (2001). The lithium abundance was derived using the Li i 6708 Å resonance doublet.
The oscillator strengths and wavelengths of the CH and CN molecular lines and the Li i doublet were taken from Drake & Pereira (2008). For the other elements which were abundances were determined using spectral synthesis, such as barium, lanthanum, and lead, the line lists were constructed using the databases of Kurucz222http://kurucz.harvard.edu/ and VALD.333http://vald.astro.uu.se/ Figures 4, 5 and 6 show the observed and synthetic spectra for the spectral regions where the abundances of carbon, nitrogen and lead were obtained.
The abundances of barium, lanthanum and lead were also determined by means of the spectral synthesis technique. The determination of barium abundance was obtained using the Ba ii lines at , and . Hyperfine and isotope splitting were taken from McWilliam (1998) and the lead abundance was derived from the Pb i line at Å where the isotopic shifts and hyperfine splitting were taken from van Eck et al. (2003). For lanthanum, the abundance was determined using the La ii lines at , , and where the isotopic shifts and hyperfine splitting were taken from Roederer & Thompson (2015) and Karinkuzhi et al. (2018).
3.3 Abundance uncertainties
Table 5 shows the uncertainties in the abundances of the chemical elements due to errors in temperature, surface gravity, microturbulence and metallicity determined previously in Section 3.1. Abundance errors were estimated by changing these parameters individually according to the uncertainties shown in Table 2. The uncertainties in the measurements of equivalent widths were estimated using the expression given by Cayrel (1988). For the Feros spectrum having a nominal resolution of 48 000 and a S/N around 100, the estimated error is of the order of 3 mÅ,. These error estimates were then taken into account in the measurements of the equivalent widths of the absorption lines of the chemical elements, provinding other values for the elemental abundances. The eighth column of Table 5 shows the final uncertainty due to uncertainties in atmospheric parameters (temperature, surface gravity, metallicity and microturbulence velocity), equivalent widths and the dispersion of the abundances. The total uncertainty was calculated as the root quadratic sum of each of these parameters considering that each individual uncertainty is independent.
For the elements analyzed through spectral synthesis technique we used the same methodology that was used to calculate the uncertainties of the abundances obtained based on measurements of the equivalent widths. We vary the temperature, surface gravity, metallicity and microturbulence velocity and then compute independently the variation of the abundances obtained by each value. In addition, for carbon, nitrogen and oxygen, we consider that their respective abundances are interdependent, that is, the uncertainty in the determination of abundance of oxygen affects the abundance of carbon and vice-versa. The uncertainties in carbon abundance also affect the abundance of nitrogen since the CN molecule was used for the determination of nitrogen abundance. For the light elements, carbon, nitrogen and oxygen, the variations of the abundance due to changes in effective temperature, surface gravity and micro-turbulent velocities are given in Table 6. Abundances based on molecular lines are little affected by microturbulence velocity variation because they are weak lines. As for surface gravity, the molecules can be affected since the molecular equilibrium is sensitive to the gas pressure. However, considering variations of only 0.2 dex for the surface gravity, the error introduced in abundance is not significant. In fact similar uncertainties were found by Tautvaiienė et al. (2000).
Table 5 shows the well known behavior that the neutral elements are more sensitive to the temperature variations, while singly-ionized elements are more sensitive to the surface gravity variations (Gray 1992). For the elements whose abundance is based on stronger lines, such as strontium, yttrium and calcium the error introduced by the microturbulence is important. This is because the equivalent widths of the moderately strong lines which lie in the saturated part of the curve of growth are larger than those synthetic lines calculated considering only thermal and damping broadening. Once microturbulent velocity is introduced in the line computation, it delays the saturation of an absorption line in the curve of growth and also becomes the major source of uncertainty in the abundance analysis based in these strong lines. For stronger absorption lines with saturated cores, the wings grow increasing the line strength and become responsible for the most of their equivalent widths. Since the equivalent width is also proportional to the Doppler broadening which is related to microturbulent velocity, one should expect that microturbulence also affects the line strength.
4 Results and Discussion
4.1 Luminosity
The luminosity of HD 55496 can be estimated based on the recently Gaia DR2 distance determination which is 49411 pc (Bailer-Jones 2018). Considering the relationship by Alonso et al. (1999) between bolometric correction (BC) and effective temperature, BC results to be 0.38. Assuming and (Bessel, 1998), we obtain log ( = 2.14 - 2.38 or = 140 - 240 L*⊙*. This range in luminosity is due to uncertain interstellar absorption because the relatively low galactic latitude ( = 5.96∘) of HD 55496. In fact, using the calibrations between , galactic coordinates and distances given by Chen et al. (1998), we obtain = 0.16 or 0.61.
The luminosity above obtained is not enough for a star develop helium-burning through the thermal pulses during AGB phase and then becoming a star rich in elements formed by slow neutron capture. Theoretical calculations show that for a star to develop the first thermal pulse the required luminosity should be around 1 800 (Lattanzio, 1986) or 1 400 (Vassiliadis & Wood, 1993).
We can also estimate the distance through the derived effective temperature and surface gravity. The relation between the distance of a star to the Sun, , and its temperature, gravity, mass, magnitude, bolometric correction (), and interstellar absorption (), is given by:
[TABLE]
Using the values of the temperature and surface gravity given in Table 2 and the values of the -magnitude and the bolometric correction given above and assuming a mass of for HD 55496 (the low metallicity suggests an old and hence lower mass nature) equation (1) can be written as
[TABLE]
Considering the three possible values for the interstellar absorption above mentioned 0.16 and 0.61 the obtained distance is, respectively, 465 and 378 parsecs. The first value, obtained for AV = 0.16 is in agreement with the recent distance determination of the Gaia DR2 release.
4.2 Radial Velocity
Table 7 shows all known measurements of radial velocity of HD 55496 available in the literature and determined in this work. We determined the radial velocity measuring the Doppler shifts of selected absorption lines. It is not clear whether the radial velocity of HD 55496 presents variations due to orbital motion. A similar conclusion has already been raised by Jorissen et al. (2005). In Figure 7 we show the radial velocity measurements obtained by Jorissen et al. (2005) where we also include our value of 315.630.23 km s*-1*. It seems that there is no significant variation of the radial velocity. Systematic radial velocity monitoring is necessary to confirm the possible binary nature of HD 55496. If HD 55496 would indeed be a binary star, the observed s-process overabundance could be due to a mass transfer happened in the past, since HD 55496 is not luminous enough to be an AGB star.
4.3 Elemental Abundances
In following Subsections we discuss the abundance pattern of HD 55496 comparing it with previous studies done for some stars in the halo and globular cluster stars, and also with the other chemically peculiar stars where heavy-element overabundances have already been reported in the literature. Figure 8 shows the abundance pattern of HD 55496 for all the elements analyzed in this work.
4.3.1 Carbon, Nitrogen, Lithium and
the isotopic ratio
As shown in Table 4, the carbon abundance is low compared to a s-process enriched star of similar metallicity such as CD-621346, with [Fe/H] = 1.590.08 and a [C/Fe] ratio of 0.86 (Pereira et al. 2012). The [N/Fe] ratio is higher than the field giant stars of similar metallicity (Gratton et al. 2000). HD 55496 has a low ratio (12C/13C = 6.01.0). The high nitrogen abundance and the low strongly suggest that significant mixing has occurred in the atmosphere of HD 55496. In fact, the [N/Fe] ratio versus [Fe/H] of dwarf stars shows no trend in the metallicity range [Fe/H] , and [N/Fe] is (Clegg et al. 1981; Tomkin & Lambert 1984; Carbon et al. 1987). As a star becomes a giant, nuclear processed material by the CN cycle which creates more 14N and more 13C relative to 12C, due to the deepening of its convective envelope, is brought from the interior to the outer layers of the star changing the surface composition.
The excess of CN is the difference between the total CN abundance and the expected primordial value. In non-evolved stars, the abundance of nitrogen scales with the iron abundance that is for metalicities between 0.3 and 2.0, therefore the [N/Fe] ratio is basically zero (Wheeler et al. 1989). For carbon we used the trend between the [C/H] versus [Fe/H] from Masseron et al. (2006). In HD 55496, the (CN) abundance sum shows an enhancement of 1.1 dex compared to stars with similar metallicity, which is seen in Figure 9 notwithstanding the low carbon abundance. The value of (CN) means that the star did not experience the triple-alpha process which would result in the carbon enhanced (see Table 4). Therefore, HD 55496 is not at the AGB phase, a conclusion already obtained in Section 4.1. In addition, the C/O ratio, which is 0.15, gives another evidence that HD 55496 is not carbon enriched.
Nitrogen can be enhanced through two mixing episodes, cool bottom processing (CBP) (Wasserburg et al. 1995)444An extra mixing mechanism that takes material from the cool base of the convective zone to the deeper layers of the stellar interior where hydrogen burning occurs. or hot bottom burning (HBB) (Sackmann & Boothroyd 1992)555Hot Bottom Burning takes place in the stars more massive than 4.0 M*⊙* where the outer convection zone reaches temperatures high enough for proton-capture nucleosynthesis occurring at the base of the outer envelope favoring the conversion of carbon to nitrogen through the CN-cycle. Lithium can also be enhanced thanks to these two mixing process, but as seen in Table 4, lithium is not enhanced. The absence of a high lithium abundance rules out CBP, since Sackmann & Boothroyd (1999) showed that CBP can account for the high lithium abundances. Therefore we should consider HBB as the source of the abundance peculiarities. In Section 4.3.2 we will also provide further support for this conclusion based on the analysis of the abundances of the elements sodium and aluminum.
4.3.2 Sodium, Aluminum, - and iron peak elements
In this Section we comment on the abundances of sodium, aluminum elements, alpha-elements and iron-peak elements in comparison with field halo stars and globular cluster stars. Inspecting Table 4, the chemical abundances of HD 55496 present an interesting aspect concerning the abundance of Na and Al. The abundance of these elements is higher compared with field stars of similar metallicity (Norris et al. 2001; Carretta et al. 2010).
In field halo stars with metallicities down to 1.0, the [Na/Fe] ratio is not higher than 0.5 (Timmes et al. 1995). However there are some exceptions, a few stars having [Na/Fe] ratios between 0.6 and 0.7 (Roederer et al. 2014). It is therefore possible that they may represent a population evaporated from globular clusters. The abundance of aluminum ([Al/Fe]) in extremely metal-poor stars is never higher than 0.5 (Fullbright 2000; Carretta et al. 2002). Contrary to field stars, both ratios in stars of globular clusters can reach ratios higher than 1.0 (Shetrone 1996; Carretta et al. 2006, just to mention a few references). The [Na/Fe] and [Al/Fe] ratios in HD 54496, with values of 0.83 and 0.76, respectively, are similar to ratios usually observed in globular cluster stars of the "second generation". In addition, the [O/Fe] ratio in HD 55496 is lower than in field stars of similar metallicity. Figure 10 shows the position of HD 55496 in the diagram [Na/Fe] and [Al/Fe] ratios versus the [O/Fe] ratio for some globular cluster stars. Based on the position on these two figures, the most likely explanation for the origin of HD 55496, is that it is an escaped globular cluster star.
As far as the other alpha-elements is concerned, the [Si/Fe], [Ca/Fe], and [Ti/Fe] ratios follow the same trend of halo stars with similar metallicity. The abundances of chromium and nickel follow, as expected, the iron abundance trend. Finally, HD 55496 shows a [Mg/Fe] ratio of +0.72, which is a value previously found in other globular clusters (Gratton et al. 2012; Carretta et al. 2013; Gratton et al. 2014; Gratton et al. 2015).
4.3.3 The s-process elements
One of the most interesting aspects of the abundance pattern of HD 55496 for a star to be considered as a “escaped globular cluster star” is the abundance of the s-process elements. From Table 4, it is evident that the light s-process elements (Sr, Y, Zr, and Mo) show a very high abundance, with a mean of [ls/Fe] = 0.960.26. The abundances of the heavy s-process elements, given by Ba, La, Ce and Nd, provide a lower [hs/Fe] ratio, which is 0.460.14, therefore the [hs/ls] ratio is 0.500.30. A similar [hs/ls] ratio of 0.35 was obtained by Karinkuzhi & Goswami (2015) for HD 55496. In addition, the [s/Fe] ratio ( ’s’ represents the mean abundance of the s-process elements Sr, Y, Zr, Mo, Ba, La, Ce and Pb) is 0.800.28. This explains why HD 55496 was previously classified as a barium star by several authors (MacConnell et al. 1972; Bond 1974; Catchpole et al. 1977; Lu et al. 1979; Sneden 1983; Luck & Bond 1991; Karinkuzhi & Goswami 2015).
According to model predictions by Busso et al (2001), for metallicities between 0.5 and 2.0 for a 1.5 AGB star and for metallicities between 0.5 and 1.0 for a 3.0 AGB star the [hs/ls] ratio has a maximum value around 1.2, which means that heavy s-process elements are the dominant products of neutron capture in AGB stars at those metallicities lower than solar. In AGB stars at near solar metallicities, the light s-process elements are the dominant products of neutron capture. Therefore the index [hs/ls] is an important measure of the neutron capture efficiency and has been widely used in the AGB nucleosynthesis models. This anti-correlation between the [hs/ls] versus metallicity is due to the operation of the 13C()16O reaction, as this neutron source is anti-correlated with the metallicity (Wallerstein 1997).
Negative [hs/ls] indexes are also predicted by models of Busso et al. (2001) for stars with [Fe/H] 1.0 however, none of the chemically peculiar binary stars where the overabundances of the elements created by the s-process are due to a mass-transfer during the AGB phase of the primary star, have negative [hs/ls] index at the same metallicity of HD 55496. This is shown in Figure 11 where the values of the [hs/ls] index for barium stars, CH stars666Chemically peculiar binary systems enhanced in carbon and s-process elements. CH stars belong to the halo population. and CEMP-s stars are displayed. HD 55496 falls in a different position compared to the other chemically peculiar binary stars, which is below the main trend of the anti-correlation between the [hs/ls] and [Fe/H]. This may raise the possibility that other enrichment processes, rather than the mass transfer hypothesis, were responsible for the overabundance in light s-process elements observed in HD55496. Another evidence that can probably rule out the hypothesis of mass-transfer is due to the inconclusive results of the behavior of the radial velocity (Section 4.2), not permitting to confirm or reject a binary nature of HD 55496. In addition, as seen in Section 4.1, we may disregard self-enrichment of HD 55496 in s-process elements because of its low luminosity to be considered an AGB star. Therefore, we suggest as the most probable explanation that the observed overabundance of the light s-process elements was due to the s-process enriched medium where HD 55496 was born and where AGB stars have played a dominant role in the chemical evolution.
In order to verify this possibility we should search among the AGB models those that predict negative [hs/ls] index for low metallicity. Models of massive AGB stars between 5.0 and 8.0 M*⊙* predict negative [hs/ls] ratios considering the reaction 22Ne(,n)25Mg as the most likely neutron source for the origin of the s-process elements (Karakas & Lattanzio 2014; Cristallo et al. 2015). In Figure 12 we show the position of HD 55496 in the diagrams [hs/ls] and [Pb/hs] versus iron abundance together with AGB model predictions given by Cristallo et al. (2015) for AGB stars between 4.0 and 6.0 M*⊙* at [Fe/H] 1.0. It is clear that based on the position of HD 55496 in both diagrams, its abundance pattern can be reproduced by an AGB star with 5.0 - 6.0 M*⊙. Therefore, if indeed HD 55496 is a “second generation” escaped globular cluster star, as suggested by the overabundances of sodium and aluminum, it is also possible that intermediate-mass AGB stars between 3.0 and 6.0 M⊙* were responsible for the production of the elements created by s-process during the time of formation of HD 55496 in a globular cluster.
4.4 Kinematics and Dynamical Origin
In order to investigate whether HD 55496 is a bound or unbound object to the Galaxy and also to investigate the globular cluster origin of HD 55496 we should first not only compare its velocity in the Galactocentric Reference Frame system () with the escape Galactic velocity but also to perform a dynamical analysis of its orbital evolution.
Knowning the distance, the proper motion and the radial velocity of HD 55496 we are able to calculate the space velocity components (,,) using the algorithm by Johnson & Soderblom (1987). We also consider the peculiar solar motion of (8.5, 13.4, 6.5) km s*-1*, as derived by Coşkunoğlu et al. (2011) and the proper motion of HD 55496 taken from Høg et al (2000). A Galactocentric solar distance of 8.5 kpc and a local standard of rest (LSR) rotation velocity relative to the GRF of 220 km s*-1* was adopted in our calculations. The results for ,, are respectively 130.0 km s*-1*, 55.3km s*-1* and 4.4km s*-1*. Considering the uncertainty of 11 pc in the distance (Section 4.1), we obtain a modulus of the star velocity in the GRF of = 1410.5 km s*-1*. This velocity is much lower than the escape Galactic velocity of km s*-1*, which corresponds to the value provided by the Galactic gravitational potential of Allen & Santillan (1991) at an heliocentric distance of 494 pc. Therefore HD 55496 appears to be a bound object to the Galaxy, however presents , which indicate that is a halo intruder object.
We performed a numerical simulation of the orbit of HD 55496 using a Bulirsh-Stoer integrator that solves the equations of motion of the star in the potential of Allen & Santillan (1991). The star shows a very elongated trajectory, with an eccentricity of . It reaches the priapsis at 0.87 kpc with a velocity of 477 km s*-1*, and the apoapsis at 10.45 kpc with a velocity of 40 km s*-1*. In order to check a possible globular cluster origin of this star, we applied the same procedure as in Pereira et al. (2017). The orbit of the star was simulated over 12 Gyr into the past, together with the orbits of 142 globular clusters from the catalog of Kharchenko et al. (2013). During the simulation, we tracked the mutual distances between the star and each cluster, and computed the encounter probabilities within a given distance . The simulation was repeated 5000 times considering different orbital parameters for the star and the clusters within their 1 uncertainties. Table 8 summarizes the results for values of corresponding to 1.5 and 5.0 cluster tidal radii. The low encounter probabilities indicate that the origin of HD 55496 is unlikely to be related to any globular cluster. The few cases displaying the largest probabilities are also associated to very large mutual encounter velocities, indicating that even in such cases the star probably did not detach from the cluster by tidal forces. From the dynamical point of view, HD 55496 has similarities with TYC 5619-109-1 (Pereira et al. 2017), and may represent another example of a possible leftover from a tidally disrupted dwarf galaxy.
5 Conclusions
Based on high-resolution optical spectroscopic data, we presented a new chemical analysis of HD 55496, a star previously known as barium star. We determined abundances of light elements, Na, Al, alpha-elements, iron-peak elements, and s-process elements. We showed that HD 55496 is a metal-poor star characterized by an enhancement of nitrogen, sodium, aluminum and s-process elements including lead. The oxygen abundance is low when compared to field metal poor stars of similar metallicity.
Our results show that HD 55496 present a Na–O anti-correlation usually observed among globular cluster stars. The star is also aluminum rich. HD 55496 presents a ratio [Al/Mg] = 0.04 which could be classified as a value between to a primordial and intermediate globular cluster populations according to Carretta et al. (2012). But since the [Al/Fe] ratio is 0.76, which is a “typical” value for an intermediate population, HD 55496 owed, probably, its abundance pattern due to a different polluters (Carretta et al. 2012).
As was classified as a barium star, HD 55496 is s-process enriched but its heavy-element abundance pattern strongly differs from those of the low-metallicity and chemically-peculiar binary stars, in sense that the abundance of Sr, Y and Zr is higher than the abundance of Ba, La, Ce and Nd. In fact, the value of the [hs/ls] ratio is 0.50, never observed for a chemically-peculiar binary star at the metallicity of HD 55496. This suggests that the main source of the neutrons is not the 13C()16O reaction but the 22Ne(,n)25Mg reaction. If this is the case, the origin of the heavy-element abundance pattern was provided by stars with masses higher than 4.0 M*⊙*, or by the intermediate mass stars which have polluted HD 55496, just like these stars were responsible for the abundances seen among the “second generation of the globular cluster stars”. In the earlier paper (Pereira et al. 2017) we showed that TYC 5619-109-1 presented a chemistry that may have been formed in a globular cluster. In addition, it was also showed that TYC 5619-109-1 is carbon poor and s-process enriched. However, HD 55496 differs from TYC 5619-109-1 in sense that the HD 55496 seems to owe its heavy-element abundance pattern due a pollution by massive stars while for TYC 5619-109-1 the contamination may have happened by a mass transfer or by a pollution of a gas already strongly enriched in s-process elements. Both HD 55496 and TYC 5619-109-1 are up to now the only two known s-process enriched globular cluster escaped stars. Finally, the kinematical and dynamical analysis shows that HD 55496 may have originated from a globular cluster however having a retrograde motion, we can not rule out the possibility that may have originated from a tidally disrupted dwarf galaxy. The low encounter probabilities with a globular clusters leaves open the origin of HD 55496.
6 Acknowledgements
We thank Verne Smith, Simone Daflon and João Victor S. Silva for fruitful discussions. NAD acknowledges Russian Foundation for Basic Research (RFBR) according to the research projects 18-02-00554 and 18-52-06004. Finally we thank the referee for the suggestions and comments that improved the quality of the paper.
The reference list from the paper itself. Each links out to its DOI / PubMed record.
- 1[1] Allen, C. & Santillan, A., 1991, R Mx AA, 22, 255
- 2[2] Allende Prieto, C., Lambert, D.L. & Asplund, M., 2001, Ap J, 556, L 63
- 3[3] Alonso, A., Arribas, S. & Martínez-Roger, C. 1999, A&AS, 140, 261
- 4[4] Antipova, L.I., Boyarchuk, A.A., Pakhomov, Yu. V. & Panchuk, V.E., 2004, A Rep, 48, 597
- 5[5] Aoki, W., Norris, J.E., Ryan, S.G., Beers, T.C.; Christlieb, N., 2004, Ap J, 608, 917
- 6[6] Bailer-Jones, C.A.L., 2018, A&A, 609, 8
- 7[7] Barbuy, B., Spite, M., Spite, F., Hill, V., Cayrel, R., 2005, A&A, 429, 1031
- 8[8] Beers, T.C. & Sommer-Larsen, J., 1995, Ap JS, 96, 175
