Open Access
Issue
A&A
Volume 665, September 2022
Article Number A46
Number of page(s) 13
Section Galactic structure, stellar clusters and populations
DOI https://doi.org/10.1051/0004-6361/202243609
Published online 07 September 2022

© T. Matsuno et al. 2022

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

This article is published in open access under the Subscribe-to-Open model. Subscribe to A&A to support open access publication.

1. Introduction

The Helmi streams were discovered by Helmi et al. (1999) as a kinematic substructure among halo stars, which is an expected signature of past galaxy accretion. They were thus one of the first pieces of evidence for the hierarchical formation of the Milky Way expected from the standard cosmology. The streams were discovered as an overdensity in the plane of angular momenta in a sample of < 100 low-metallicity stars. Later studies also confirmed their existence (Chiba & Beers 2000; Kepley et al. 2007; Myeong et al. 2018; Koppelman et al. 2018). Their progenitor was likely a dwarf galaxy with a stellar mass of ∼108M and a mean metallicity of [Fe/H] ∼ − 1.5 that was accreted to the Milky Way around 5 − 8 Gyr ago (Kepley et al. 2007; Koppelman et al. 2019).

Even though the streams were discovered more than 20 years ago, their chemical properties have not received much attention until recently. Roederer et al. (2010) studied the chemical abundances of 12 stars in the stream, concluding that the stars have similar abundance ratios as the rest of the halo stars. In recent years, Limberg et al. (2021), Gull et al. (2021), Aguado et al. (2021), and Nissen et al. (2021) conducted chemical analysis for the Helmi stream members. While Gull et al. (2021) reach the same conclusion as Roederer et al. (2010), the other studies hint at possible chemical peculiarities of the streams. For example, the Sr abundance seems to be much lower than typical for halo stars at [Fe/H] ≲  − 2 (Aguado et al. 2021); two stars in a binary system in the Helmi streams, namely G112-43 and G112-44, have enhanced ratios of [Mn/Fe], [Ni/Fe], [Cu/Fe], ande [Zn/Fe] compared to other accreted stars in the Milky Way (Nissen et al. 2021); and the [α/Fe] ratio seems to be decreasing at [Fe/H] ≳ − 2 (Limberg et al. 2021; Aguado et al. 2021).

Many of the previous comparisons made use of heterogeneous literature compilation. They were therefore limited by systematic uncertainties, which can be as large as 0.4 dex in [Fe/H] and 0.2 dex in [X/Fe]. This can easily obscure or falsely create abundance peculiarities, as we will see later in this paper. Another complication arises because the literature sample was often a mixture of accreted and in situ stars. As Nissen & Schuster (2010, hereafter, NS10) showed from their homogeneously analyzed samples, the Milky Way halo contains two major populations, accreted and in situ stars, which differ in abundance ratios. Without understanding the origin of the comparison sample, it remains unclear if the chemical peculiarities of the Helmi stream stars, if any, indicate simply that they are accreted or that the progenitor had a unique evolution. Nissen et al. (2021) is the only exception in the sense that they compared Helmi stream members with accreted stars using homogeneously derived chemical abundance. However, larger samples could aid our understanding of the interesting abundance pattern reported in the two stars in a binary system.

Throughout this study, we aim to precisely characterize the chemical abundances of the Helmi stream stars and to compare them with homogeneous abundances of other components in the Milky Way halo, including Gaia-Enceladus, in situ stars (NS10; Reggiani et al. 2017, hereafter R17), and Sequoia (Matsuno et al. 2022, hereafter, Paper I). The precise chemical abundance ratios tell us about the progenitor’s properties, but the chemical distinctness of kinematic substructures, if found, would also enable us to “chemically tag” stars originating from the same progenitor. In Paper I, we show that a combination of a line-by-line abundance analysis for a carefully selected sample and homogenization of literature abundances using standard stars is a powerful way to precisely characterize the chemical properties of kinematic substructures.

We studied the chemical abundance of the Helmi stream stars in the same way as in Paper I. Namely, we analyzed high signal-to-noise ratio spectra for 11 stars in the Helmi streams, measured their chemical abundance through a line-by-line analysis with careful treatment of uncertainty, and compared them with other halo stars in NS10 and R17 after homogenizing all the abundances. We describe the observation, data reduction, targets, and stars in the literature in Sect. 2. In Sect. 3, after briefly explaining our approach for the chemical abundance measurement and homogenization, we present chemical abundance ratios of the Helmi stream stars. After providing discussions in Sect. 4, we summarize our findings in Sect. 5.

2. Data

2.1. Observation

We obtained high-resolution spectra for ten stars with the High Dispersion Spectrograph (Noguchi et al. 2002) on the Subaru telescope1. The selection and properties of the targets are described in the following subsection. The observations were conducted on August 9-10, 2020, with the standard setup of the HDS (StdYd), which yields a wavelength coverage from ∼ 4000 Å to ∼ 6800 Å. We used the image slicer #2 (Tajitsu et al. 2012), which provides a resolving power of R ∼ 80 000. The data reduction was performed using an IRAF script2, hdsql3, including CCD linearity correction, scattered light subtraction, aperture extraction, flat-fielding, wavelength calibration, and heliocentric velocity correction. Information about the spectra is summarized in Table 1.

Table 1.

Summary of the data.

We also searched archives for high signal-to-noise ratio and high-resolution spectra for stars that satisfy our selection criteria described in the following subsection. We find that LP894-3 satisfies our selection and has spectra in the Subaru telescope archive. We, therefore, analyzed the archival spectrum of this star. The spectrum used is the same as that used by Ishigaki et al. (2012).

2.2. Targets

Targets were selected based on the orbital angular momenta of the stars following Koppelman et al. (2019). Namely, the selection box was defined with the component of the angular momentum along the z-axis of the Milky Way (Lz) and Lperp, which is given by L x 2 + L y 2 $ \sqrt{L_x^2+L_{\mathit{y}}^2} $. We also restricted our sample to those around the turn-off region using Gaia DR2 photometry. The exact selection criteria were 1000 < Lz/kpc km s−1 < 1500 and 1750 < Lperp/kpc km s−1 < 2600 (the kinematic selection, corresponding to the box A of Koppelman et al. 2019), and 0.5 < GBP − GRP < 0.95 and 3.2 < Gabs < 5.8 (the color-magnitude selection), where GBP, GRP, and Gabs are GaiaBP and RP magnitude, and the absolute magnitude in the GaiaG band, respectively.

Initially, for the target selection, we used Gaia DR2 astrometry, photometry, and radial velocity (Gaia Collaboration 2018), and Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST) survey DR5 radial velocity. We assumed U = −11.1 km s−1, W = 7.3 km s−1 (Schönrich et al. 2010), R0 = 8.21 kpc (McMillan 2017), and z0 = 20.8 pc (Bennett & Bovy 2019) for the solar velocity and position. Adopting the proper motion measurement for Sgr A* (Reid & Brunthaler 2004), we assumed the sun is moving at 245.3 km s−1 in the direction of the Galactic rotation. When calculating the orbital energies of the stars, we adopted the Milky Way potential of McMillan (2017).

Information about the kinematic and photometric properties of the targets is summarized in Tables 2 and 3, and visualized in Figs. 1 and 2. We note that we updated the kinematic and photometric information with Gaia EDR3 photometry and astrometry (Gaia Collaboration 2021). The radial velocities were also updated to the values we measured from our high-resolution spectra by comparing the observed wavelengths of neutral Fe lines with those from laboratory experiments.

thumbnail Fig. 1.

Distribution of stars in kinematic spaces, namely E − Lz, Lperp − Lz, and vy − vz. The insets in the left two panels show better distribution of the stars that are part of the Helmi streams. The Helmi stream selection from Koppelman et al. (2019) is shown as a dashed box in the inset of the middle panel. Stars are hatched according to the subgroups found by Dodd et al. (2022). Three stars (2447_5952, 4998_5552, and LP894-3) are not associated with either of the subgroups. An open symbol is assigned only to 2447_5952 since it has much larger energy than the other Helmi stream stars. We used filled symbols for stars from NS10 and R17 if they satisfy the kinematic selection for Gaia-Enceladus (see Paper I).

thumbnail Fig. 2.

Location of the stars in the color-magnitude diagram. The symbols follow those of Fig. 1. The four isochrones are from the PARSEC tracks with the age of 12 Gyr and [Fe/H]= − 2.0,   − 1.5,   − 1.0,  and  − 0.5.

Table 2.

Property of the targets.

Table 3.

Kinematics of the targets.

Table 3 also includes information about subgroups among the Helmi stream stars. We consider here two ways of subdividing the Helmi stream stars. One is by the sign of vz; Kepley et al. (2007) and Koppelman et al. (2019) used the asymmetry in the number of stars in each subgroup to estimate the accretion time of the progenitor of the Helmi streams. Another subdivision follows that by Dodd et al. (2022) in the space of Lz and Lperp. They show that these two groups of stars with slightly different kinematics reflect the effect of an orbital resonance in the Milky Way potential. In either definition, the two subgroups are considered parts of the same substructure, the Helmi streams. In all figures, symbols are hatched according to the subdivision by Dodd et al. (2022).

2.3. Literature

We compared the chemical abundance of the Helmi stream stars with those in literature. Since our approach described in Sect. 3 is identical to that used in Paper I, we can naturally compare our abundance in the present study with Gaia-Enceladus and in situ stars from NS10 and R17 and with Sequoia stars from Paper I. Since abundances from these studies were homogenized using standard stars (see Paper I), precise comparisons are possible. We refer to Paper I for more details.

Among the stars studied by NS10 and R17, we find two stars that satisfy our kinematic selection criteria for the Helmi streams. These two stars are G112-43 and G112-44 and were studied in detail by Nissen et al. (2021) as part of the Helmi streams. Although they are considered to form a binary system, their motions within the system do not affect the estimated orbital properties within the Milky Way since the separation between the two stars is large (> 2000 au).

3. Chemical abundances

3.1. Abundance analysis

We provide here a brief description of our analysis, and readers can refer to Paper I for a complete and detailed explanation. The abundance measurements were based on equivalent width measurements, where we assumed a Voigt profile for the line shape. We applied spectral synthesis for Si, Mn, Zn, and Y. Hyperfine structure splitting was considered for Na, Mn, and Ba. Abundances were estimated based on the assumption of local thermodynamic equilibrium (LTE) using MOOG (Sneden 1973) and based on one-dimensional plane-parallel MARCS atmosphere models (Gustafsson et al. 2008). We applied non-LTE corrections of Lind et al. (2011) to the obtained Na abundance through the INSPECT database45.

The effective temperature (Teff) was determined by minimizing the correlation between the excitation potentials and line-by-line abundances of neutral iron lines. The microturbulent velocity (vt) was determined similarly by minimizing the correlation between reduced equivalent widths and abundances. The surface gravity (log g) was determined from temperature, luminosity, and stellar mass, where the luminosity was obtained by applying a bolometric correction of Casagrande & VandenBerg (2018) to the absolute magnitude in the GaiaG band, and mass was derived from isochrone fitting in the color-magnitude diagram. This method of log g determination has an advantage over the spectroscopic method, requiring the ionization equilibrium between Fe I and Fe II, in that it does not strongly depend on the assumed Teff. For the stellar parameter determination, we used the simple mean of iron abundances derived from individual Fe II lines because Fe II abundances are less affected by the non-LTE effect and temperature uncertainty than Fe I abundances. We took correlated uncertainties into consideration for the estimates of uncertainties. The derived parameters are presented in Table 4.

Table 4.

Stellar parameters.

The elemental abundances are the weighted average of the abundances derived from the individual lines. The weights were determined following the prescription by Ji et al. (2020). The weights on individual lines as well as line-by-line abundances and their sensitivities to stellar parameters are summarized in Table 5. When considering [X/Fe], we used the Fe abundance from the same ionization stage as the species X.

Table 5.

Linelist and line-by-line abundance.

All of the analysis was conducted relative to HD59392, for which we adopted Teff = 6012 K, log g = 3.954, vt = 1.4 km s−1, and [Fe/H]sp = −1.6 (NS10; Paper I). The elemental abundance of this star was taken from NS10 and Nissen & Schuster (2011). In Paper I, we validated our results using a few standard stars and confirm that our abundances are on the same scale as NS10, Nissen & Schuster (2011) and R17. We note that one of the standard stars was G112-43 in the Helmi streams, for which Nissen & Schuster (2011) and Nissen et al. (2021) reported high Zn abundance. Since we do not find any systematic offset in the derived abundance for this star when comparing with NS10 and Nissen & Schuster (2011), our analysis also confirms its high Zn abundance.

In case no lines were detected for a species, we provided 1-σ and 3-σ upper limits on the abundance, which were estimated from the expected uncertainty on equivalent widths from the equation in Cayrel (1988). We used the Si I line at 6237 Å and the Y II line at 4884 Å to place upper limits. The derived abundances and upper limits are presented in Table A.1.

As an additional test, we investigated the effect of slightly different atmospheric structures due to different α-element abundances in model atmospheres, since the sample includes stars showing different [α/Fe] compared to typical halo stars. We used two MARCS model atmospheres with Teff = 6000 K, log g = 4.0, vt = 1.5 km s−1, and [Fe/H]= − 1.5 differing in [α/Fe], one of which has the “standard composition” at this metallicity with [α/Fe]= + 0.4 and the other has a solar [α/Fe]. Assuming the chemical composition of HD59392, we first computed equivalent widths of absorption lines using the standard α-enhanced model atmosphere. We then derived abundances with the other model from the computed equivalent widths. The maximum difference between the assumed and derived abundances is found in strong lines to be ∼0.02 dex. Since the variation in [α/Fe] among our sample is ∼0.2 dex and since every line is affected in the same direction, the effect of different [α/Fe] in model atmospheres is at most 0.01 dex when discussing [X/H] and even smaller when discussing [X/Fe].

3.2. Results

Figure 3 shows [Na/Fe] and [Mg/Fe] ratios of the Helmi stream stars together with the halo stars from NS10 and R17 and the Sequoia stars from Paper I, all of which are on the same abundance scale. We note that we indicate Gaia-Enceladus stars in the literature with filled symbols in figures. At −1.8 ≲ [Fe/H] ≲  − 1.6, the Helmi stream stars show similar Na and Mg abundances to Sequoia stars, having clearly lower [Na/Fe] and [Mg/Fe] values compared to other halo stars including Gaia-Enceladus stars. The Helmi stream stars with [Fe/H] ≳ − 1.5 also have lower values in these abundance ratios than in situ halo stars. Although they are less separated from Gaia-Enceladus stars at this high metallicity, they still tend to show lower [Mg/Fe] values compared to Gaia-Enceladus stars. The most metal-poor star 6615_9776 also seems to show a lower abundance of Na and Mg compared to other halo stars. We, however, note that this star is in the metallicity range, [Fe/H]<  − 2.1, where R17 measured the abundance of stars using a standard star that is different from what they used for stars with [Fe/H] > − 2.1. Hence, there is no guarantee that the abundance of 6615_9976 and R17 stars with similar metallicity are on the same scale. Readers can refer to Paper I for the effect of different standard stars (see the results of abundance comparison for HIP28104 in Sect. 3 of Paper I).

thumbnail Fig. 3.

[Na/Fe] and [Mg/Fe] abundance ratios of the Helmi stream stars (red symbols). Also plotted are halo stars from NS10 and R17, with Gaia-Enceladus stars shown with filled symbols. G112-43 and G112-44 from NS10 are shown with red stars since their kinematics are consistent with the Helmi streams. We also show kinematically selected Sequoia stars from Paper I.

Figure 4 shows Si, Ca, and Ti abundances, which are also known to be different between stars with different origins (NS10; Paper I). The behavior of the Helmi stream stars in Ca and Ti is similar to that in Na and Mg; the stars show clearly lower abundance at low metallicity than most halo stars, while they become more indistinguishable from Gaia-Enceladus stars at higher metallicity. For Si, the picture is unclear because the Si lines become too weak at [Fe/H] ≲  − 1.5.

thumbnail Fig. 4.

α-element abundances of the Helmi stream stars. The symbols follow those of Fig. 3. When a star does not have any detectable lines for an element, we indicate the abundances corresponding to the 1-σ and 3-σ upper limits on the equivalent widths, respectively shown with the location of the symbols and the upper end of the error bars.

The observed trends in the α-elements confirm the results of Limberg et al. (2021) and Aguado et al. (2021) in the sense that the Helmi stream stars have lower α-element abundance compared to general Milky Way stars. We added here another important result; the α-element abundances of the Helmi stream stars are even lower than another accreted population, the Gaia-Enceladus stars, especially at [Fe/H] ≲  − 1.5. Moreover, contrary to the conclusion of Limberg et al. (2021) and Aguado et al. (2021) that the [α/Fe] ratios decrease with metallicity, the evolution in α-elements (Mg, Ca, and Ti) with metallicity observed in the present study is rather flat from [Fe/H] ∼ − 2.1 to [Fe/H] ∼ − 1.2. The low α abundance of the Helmi stream stars is also in contrast with Roederer et al. (2010), who concluded that the Helmi stream stars have a similar abundance as other halo stars. We consider that our precise and homogeneous abundance comparisons enable us to detect the abundance difference between the Helmi stream stars and Gaia-Enceladus stars and precisely depict the chemical evolutionary track followed by the Helmi streams. We further discuss the importance of the homogeneous abundance analysis in Sect. 4.2.

The abundances of elements near the iron-peak, namely Cr, Mn, Ni, and Zn, are shown in Fig. 5. It has been reported that Cr, Ni and Zn abundances are different between accreted stars and in situ stars (Nissen & Schuster 2011). We, therefore, expect that these abundances are different among the various kinematic substructures with different α-element abundances. We, however, do not see apparent abundance differences between the Helmi stream stars and the other halo stars in Cr and Ni at our measurement uncertainty. On the other hand, Zn appears depleted in the Helmi streams at low metallicity, while there is a large variation at high metallicity. The Helmi stream stars tend to show a low value of [Mn/Fe] at [Fe/H] ≲  − 1.5 while they seem to be consistent with typical Gaia-Enceladus stars at higher metallicity.

thumbnail Fig. 5.

Abundances of Cr, Mn, Ni, and Zn. Symbols follow Fig. 3.

The high Zn abundance of G112-43 and G112-44 was reported in Nissen & Schuster (2011) and confirmed in Nissen et al. (2021). It is also supported for G112-43 by our analysis in Paper I as described in Sect. 3.1. The other Helmi stream stars in the present study, on the other hand, do not share a similar Zn enhancement. We attempt to interpret the Zn abundance evolution of the Helmi streams in Sect. 4.1.

Abundances of neutron-capture elements are shown in Fig. 6. The Helmi stream stars are clearly deficient in Y at low metallicity. This behavior of Y resembles that of Sr reported by Aguado et al. (2021), which is expected since both are light neutron-capture elements and since these two elements usually show similar trends with [Fe/H]. It is hard to conclude if the Ba abundance evolution of the Helmi streams is distinct from that of Gaia-Enceladus or other halo stars because of large scatter. However, we remark that the Helmi stream stars tend to be on the lower side in [Ba/Fe] at given [Fe/H] compared to other halo stars.

thumbnail Fig. 6.

Abundances of Y and Ba. Symbols follow Fig. 3.

As described in Sect. 2.2, symbols in Figs. 36 are hatched according to the subgroups discussed by Dodd et al. (2022). Both subgroups show similar abundance peculiarities compared to other halo stars (lower abundances of α-elements and Y). This result supports that both subgroups can be regarded as a part of the Helmi streams. We made a similar inspection between stars with positive and negative vz and reached the same conclusion.

We note that the only star with an exceptionally high value of En (2447_5952) has a slightly different elemental abundance than the other Helmi stream stars. The different behavior of this star is especially prominent in Ca, Ti and neutron-capture elements (Figs. 4 and 6). This might indicate that the Helmi streams do not extend toward high En, which is consistent with the recent identification of Helmi stream members based on clustering analysis of halo stars (Lövdal et al. 2022; Ruiz-Lara et al. 2022).

4. Discussion

4.1. Properties of the Helmi stream progenitor

The chemical properties of the Helmi streams can be summarized as follows. The abundances of elements usually produced by explosions of massive stars (e.g., Na and α-elements) are generally lower than most of the other halo stars, including Gaia-Enceladus stars. While the difference to Gaia-Enceladus stars becomes less clear at [Fe/H] ≳ − 1.5, the Helmi stream stars still tend to be on the lower side of the distribution in [Mg/Fe], [Ca/Fe], and [Ti/Fe]. Zn and Y also show similar trends as these elements, suggesting they originate in massive stars. This chemical abundance pattern is similar to that seen for Sequoia stars (Paper I).

The low abundance of α-elements and Na are usually attributed to Fe enrichments from type Ia supernovae (SNe Ia). Interestingly, SNe Ia seem to have started to operate in the Helmi stream progenitor already at very low metallicity ([Fe/H]<  − 2.2) since even the lowest-metallicity star in the sample (6615_9776) has low abundances of α-elements and Na. Assuming that every system has the same abundance ratios among Na, α-elements, and Fe before the onset of SNe Ia, we can tentatively provide an upper limit for the “knee” metallicity as [Fe/H]knee ≲ −2.4 from the ∼0.2 dex lower [Mg/Fe] of 6615_9776 than the halo stars from R17.

Another interesting feature in α-elements is that there is no apparent decreasing trend in [α/Fe] with metallicity. Moreover, a few stars at [Fe/H] ∼ − 1.2 seem to have higher [α/Fe] than the stars at [Fe/H] ∼ − 1.8, although it remains to be seen if they are actually the Helmi stream members rather than contaminants from heated thick disk or other accreted populations. If this almost flat or slightly increasing trend is confirmed with a larger sample of homogeneous abundances, it would resemble the observed trend in the Large Magellanic Could (LMC) by Nidever et al. (2020). These authors suggest that such a chemical evolution trend cannot be reproduced in models with a constant star formation efficiency. They suggest that such a trend needs a starburst at a late phase. Although the progenitor of the Helmi streams might have experienced a similar starburst at a late phase, the old age of the majority of the Helmi stream stars (Koppelman et al. 2019; Ruiz-Lara et al. 2022) suggests that the transition from quiescent to bursty star formation occurred on a shorter timescale than that in LMC. A future chemodynamical study with precise stellar age as well as chemical evolution modeling would be necessary to confirm this scenario.

The Y abundance of the Helmi stream stars is particularly low at [Fe/H] ≲  − 1.8. Together with the results of Aguado et al. (2021) on Sr abundance, this suggests that the Helmi stream stars have very low abundances of light neutron-capture elements. Such a low light neutron-capture element abundance is not common among the Milky Way stars but is seen in low-mass dwarf galaxies such as Draco and Ursa Minor dwarf spheroids and most of the ultra-faint dwarf galaxies (e.g., Frebel & Norris 2015).

Several nucleosynthesis processes have been proposed as production sites of these light neutron-capture elements, including r-process in neutron star mergers (Wanajo et al. 2014; Watson et al. 2019), in magneto-rotational supernovae (Winteler et al. 2012), or in collapsars (Siegel et al. 2019), s-process in low- to intermediate-mass stars (Karakas & Lattanzio 2014), weak s-process in rapidly rotating massive stars (Frischknecht et al. 2012; Choplin et al. 2018), and weak r-process in electron-capture supernovae (Wanajo et al. 2011). It is not yet clear which process is the dominant source of the elements in the early Universe (see discussions by, e.g., Côté et al. 2019; Prantzos et al. 2018; Kobayashi et al. 2020). Since we here discuss the low-metallicity end of the sample, the production of neutron-capture elements would not be dominated by low- to intermediate-mass stars (e.g., de los Reyes et al. 2022). One possible explanation for the low light neutron-capture element abundances of the Helmi streams is that, as a result of the low stellar mass of the galaxy, the progenitor did not experience rare r-process nucleosynthesis events, such as neutron star mergers, electron capture supernovae, and magneto-rotational supernovae. In this case, a small amount of light neutron-capture elements could be produced by rapidly rotating massive stars (Hirai et al. 2019; Tarumi et al. 2021). Another explanation is that the progenitor dwarf galaxy had a small number of rotating massive stars. The small number of rotating massive stars might be a result of the top-light initial mass function in dwarf galaxies (Weidner & Kroupa 2005), or different distribution of initial rotation velocity of stars. The observational indication by Gull et al. (2021) that metal-poor stars of the Helmi streams show r-process abundance pattern in neutron-capture elements heavier than Ba might favor the second possibility. However, it is necessary to investigate the abundance pattern of light neutron-capture elements in order to understand the cause of the low light neutron-capture element abundance of the Helmi streams. A larger sample of low-metallicity Helmi stream stars with neutron-capture element abundances would also be welcomed. They would enable us to constrain the property of the nucleosynthesis processes, such as their event rates, by studying how neutron-capture elements were enriched as a function of metallicity (e.g., Tsujimoto et al. 2017).

The [Zn/Fe] values observed at high metallicity are also noteworthy. While the binary pair G112-43 and G112-44 have high Zn abundance (Nissen et al. 2021), other Helmi stream stars do not share such a high abundance, resulting in a large star-to-star variation in [Zn/Fe]. Although Nissen et al. (2021) noted that the binary pair is also enhanced in Mn, Ni, and Cu, the measurement uncertainties in the present study are not high enough to see if there are significant scatters in Mn and Ni abundances among the Helmi stream stars. A significant dispersion in [Zn/Fe] at high metallicity is also reported in the Sculptor dwarf spheroidal galaxy (Skúladóttir et al. 2017). While Hirai et al. (2018) provided a theoretical calculation of Zn enrichments in dwarf galaxies, assuming electron capture supernovae as one of the sources of Zn, they find it challenging to produce a large scatter at high metallicity. Based on iron-group element abundances, Nissen et al. (2021) suggest pure helium detonation type Ia supernovae could be a promising nucleosynthesis event producing the unique abundance pattern of the binary pair. It remains to be seen whether chemical evolution models including this type of SNe Ia can explain the Zn abundance variation within a galaxy. A larger sample of stars with precise multi-element abundances would be helpful to establish, for example, if the binary pair is a chemical outlier or a tail of [Zn/Fe] distribution in the system and if the Zn abundance correlates with abundances of other elements.

4.2. Helmi streams in the literature

In this section, we highlight the importance of homogeneous chemical abundance through a comparison of our results with those presented in the literature, by Roederer et al. (2010), by Limberg et al. (2021), by Gull et al. (2021), and by Aguado et al. (2021)6. All but Limberg et al. (2021) conducted new observations of the Helmi stream stars, while Limberg et al. (2021) used the data from GALAH DR3. All, including Limberg et al. (2021), compared the chemical abundance of the Helmi stream stars with a literature compilation without homogenizing the chemical abundances.

Figure 7 compares the abundances of the Helmi stream stars from those literature studies and from our present work, with those of stars in the Milky Way from various literature sources (SAGA database; Suda et al. 2008, 2011; Yamada et al. 2013). To confirm the membership to the Helmi streams, we recomputed angular momenta of the stars in Roederer et al. (2010), Gull et al. (2021), and Aguado et al. (2021) consistently and reselected the Helmi stream members. As there is no report of radial velocity in Aguado et al. (2021), we adopted radial velocity measurements by SDSS or LAMOST for their stars. For the stars in common between Roederer et al. (2010) and Gull et al. (2021), we prioritized the measurements by the latter study over the former. We also included in the figure abundances that we obtained for three stars studied by Aguado et al. (2021) by reanalyzing their UVES spectra. For this reanalysis, we redetermined stellar parameters following the same procedure as described in Paper I and estimated the abundances using these parameters.

thumbnail Fig. 7.

Chemical abundance of the Helmi stream stars from literature, namely Roederer et al. (R10; 2010), Aguado et al. (A21; 2021), and Gull et al. (G21; 2021). We also include Helmi stream stars in the present study and in NS10. The background gray points are from the SAGA database (Suda et al. 2008, 2011; Yamada et al. 2013). The UVES spectra of three stars from Aguado et al. (2021) were reanalysed and abundances from the reanalysis are included in the figure. Abundances from our reanalysis and those from Aguado et al. (2021) are connected with a solid line for each star.

We first emphasize the importance of homogeneity in abundance by using the three stars from Aguado et al. (2021) that we reanalyzed. As can be seen in Fig. 7, stars from Aguado et al. (2021) seem to have higher α-elements abundances compared to the Helmi stream stars in the present study if we adopt the reported abundances. On the other hand, our reanalysis shows that the three stars actually have similar abundances as the sample in the present study, although the uncertainties are larger due to different spectral quality. To investigate the effect of different stellar parameters, we fixed Teff and log g to the values adopted by Aguado et al. (2021), determined vt and metallicity, and then measured abundances instead of redetermining the four stellar parameters. The results do not change significantly compared to what we obtained from the full reanalysis. The source of difference thus could be due to different atomic data, line selection, and/or spectral synthesis software.

The different abundance from Aguado et al. (2021) and our reanalysis is not too surprising since it has been known that systematic uncertainty can be significant when abundances from different studies are compared without homogenization. For example, G112-43 is a well studied star and its abundance is available in a number of literature (Ryan et al. 2001; Charbonnel & Primas 2005; Zhang et al. 2009; Ishigaki et al. 2010, 2012, 2013; Nissen & Schuster 2010, 2011; Yan et al. 2016)7. The reported metallicity ([Fe/H]I) and [Mg/Fe] of this star respectively range from −1.52 (Ishigaki et al. 2010) to −1.14 (Zhang et al. 2009), and from 0.16 (Zhang et al. 2009) to 0.37 (Ishigaki et al. 2010). These ranges are significantly larger than the typical uncertainty reported. We note that we derive [Fe/H]I = −1.254 and [Mg/Fe] = 0.179 in Paper I and NS10 (adopted for figures) derive [Fe/H]I = −1.25 and [Mg/Fe] = 0.21.

Figure 7 demonstrates that our sample has the smallest uncertainty, which was made possible thanks to the high quality of the spectra. However, even with the high precision, the chemical distinctness of the Helmi stream stars is less clear than what we see in Figs. 36. This is because the comparison sample from the SAGA database is not necessarily on the same abundance scale as ours, because its distribution is broadened due to systematic uncertainty, and because it can contain poor quality data.

In conclusion, systematic uncertainties have hampered a clear chemical characterization of the Helmi stream stars in previous studies. Thanks to our approach of abundance analysis and homogenization of abundances from literature (see, Paper I), we can detect clear but small abundance difference between the Helmi stream stars and other halo stars that include Gaia-Enceladus stars.

5. Summary

Through a differential abundance analysis, we show that the Helmi stream stars clearly depict lower [X/Fe] in Na, Mg, Ca, Ti, Zn, and Y, compared to the other halo stars, including stars accreted from Gaia-Enceladus. While the distinction from Gaia-Enceladus stars is clearer at [Fe/H] ≲  − 1.5, a part of the Helmi stream stars start to overlap with Gaia-Enceladus stars in elemental abundance ratios at high metallicities.

This is the first time that we see the chemical distinctness of the Helmi stream stars over a wide metallicity range. This result was made possible thanks to the precise abundances derived in this study and to the homogenization of abundances across NS10, R17 and the present study. As in Paper I, the present investigation of the Helmi streams highlights the importance of homogenized abundances for chemical identification of individual building blocks of the Milky Way.

The observed low α-element abundances of the Helmi streams are likely due to the Fe enrichments from SNe Ia at low metallicity. There seems to be a contribution of SNe Ia to the chemical evolution of the progenitor already at very low metallicity ([Fe/H] ≲  − 2.4), indicating that the star formation proceeded with low efficiency in early times. Star formation at later times might have been bursty since [α/Fe] remains almost constant or slightly increases toward high metallicities according to Nidever et al. (2020); such a feature may require a late starburst.

The origin of the extremely low Y abundances at low metallicity remains unclear. Our results, as well as those of Aguado et al. (2021) on Sr, strongly suggest that the Helmi streams have very low abundances of light neutron-capture elements at low metallicity. It would be interesting to obtain detailed abundance patterns over light neutron-capture elements in order to constrain the nucleosynthesis processes that produce this feature.

Further observations would also be able to aid understanding of the Zn abundance variation at high metallicity. The binary pair G112-43 and G112-44 are the only stars with clear Zn enhancements, and all the other Helmi stream stars at [Fe/H] ≳ − 1.5 have Zn abundances comparable to or lower than Gaia-Enceladus stars. It would be interesting to investigate the detailed [Zn/Fe] distribution with a lager number of stars.


1

Gaia EDR3 1909569058536197760 was also observed but not analyzed since it is a double-lined spectroscopic binary.

2

IRAF is distributed by the National Optical Astronomy Observatory, which is operated by the Association of Universities for Research in Astronomy (AURA) under a cooperative agreement with the National Science Foundation.

4

Data obtained from the INSPECT database, version 1.0 (www.inspect-stars.com).

5

The Na abundances in the comparison literature are in LTE (NS10) and in non-LTE (R17). This inconsistency does not affect our conclusions since we put all the abundances onto the same scale and since NS10 used Na lines that require small non-LTE corrections (≲0.1 dex).

6

We consider five stars observed with high-resolution spectrographs (HORuS and UVES) for the Aguado et al. (2021) sample.

7

The data are collected through the SAGA database.

Acknowledgments

We thank the anonymous referee for his/her constructive comments, which helped us with improving the clarity of the discussion. We thank David Aguado for sharing their UVES spectra with us. This research has been supported by a Spinoza Grant from the Dutch Research Council (NWO). W.A. and M.N.I. were supported by JSPS KAKENHI Grant Number 21H04499. MNI is supported by JSPS KAKENHI Grant Number 20H05855. Z.Y. acknowledges support from the French National Research Agency (ANR) funded project “Pristine” (ANR-18-CE31-0017) and the European Research Council (ERC) under the European Unions Horizon 2020 research and innovation programme (grant agreement No. 834148). This research is based in part on data collected at Subaru Telescope, which is operated by the National Astronomical Observatory of Japan. We are honored and grateful for the opportunity of observing the Universe from Maunakea, which has the cultural, historical and natural significance in Hawaii. Part of the data were retrieved from the JVO portal (http://jvo.nao.ac.jp/portal/) operated by ADC/NAOJ.

References

  1. Aguado, D. S., Myeong, G. C., Belokurov, V., et al. 2021, MNRAS, 500, 889 [NASA ADS] [Google Scholar]
  2. Bennett, M., & Bovy, J. 2019, MNRAS, 482, 1417 [Google Scholar]
  3. Casagrande, L., & VandenBerg, D. A. 2018, MNRAS, 479, L102 [NASA ADS] [CrossRef] [Google Scholar]
  4. Cayrel, R. 1988, in The Impact of Very High S/N Spectroscopy on Stellar Physics, eds. G. Cayrel de Strobel, & M. Spite, 132, 345 [Google Scholar]
  5. Charbonnel, C., & Primas, F. 2005, A&A, 442, 961 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  6. Chiba, M., & Beers, T. C. 2000, AJ, 119, 2843 [NASA ADS] [CrossRef] [Google Scholar]
  7. Choplin, A., Hirschi, R., Meynet, G., et al. 2018, A&A, 618, A133 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  8. Côté, B., Eichler, M., Arcones, A., et al. 2019, ApJ, 875, 106 [Google Scholar]
  9. de los Reyes, M. A. C., Kirby, E. N., Ji, A. P., & Nuñez, E. H. 2022, ApJ, 925, 66 [NASA ADS] [CrossRef] [Google Scholar]
  10. Dodd, E., Helmi, A., & Koppelman, H. H. 2022, A&A, 659, A61 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  11. Frebel, A., & Norris, J. E. 2015, ARA&A, 53, 631 [NASA ADS] [CrossRef] [Google Scholar]
  12. Frischknecht, U., Hirschi, R., & Thielemann, F. K. 2012, A&A, 538, L2 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  13. Gaia Collaboration (Brown, A. G. A., et al.) 2018, A&A, 616, A1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  14. Gaia Collaboration (Brown, A. G. A., et al.) 2021, A&A, 649, A1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  15. Green, G. M., Schlafly, E., Zucker, C., Speagle, J. S., & Finkbeiner, D. 2019, ApJ, 887, 93 [NASA ADS] [CrossRef] [Google Scholar]
  16. Gull, M., Frebel, A., Hinojosa, K., et al. 2021, ApJ, 912, 52 [NASA ADS] [CrossRef] [Google Scholar]
  17. Gustafsson, B., Edvardsson, B., Eriksson, K., et al. 2008, A&A, 486, 951 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  18. Helmi, A., White, S. D. M., de Zeeuw, P. T., & Zhao, H. 1999, Nature, 402, 53 [Google Scholar]
  19. Hirai, Y., Saitoh, T. R., Ishimaru, Y., & Wanajo, S. 2018, ApJ, 855, 63 [NASA ADS] [CrossRef] [Google Scholar]
  20. Hirai, Y., Wanajo, S., & Saitoh, T. R. 2019, ApJ, 885, 33 [NASA ADS] [CrossRef] [Google Scholar]
  21. Ishigaki, M., Chiba, M., & Aoki, W. 2010, PASJ, 62, 143 [NASA ADS] [Google Scholar]
  22. Ishigaki, M. N., Chiba, M., & Aoki, W. 2012, ApJ, 753, 64 [NASA ADS] [CrossRef] [Google Scholar]
  23. Ishigaki, M. N., Aoki, W., & Chiba, M. 2013, ApJ, 771, 67 [NASA ADS] [CrossRef] [Google Scholar]
  24. Ji, A. P., Li, T. S., Hansen, T. T., et al. 2020, AJ, 160, 181 [NASA ADS] [CrossRef] [Google Scholar]
  25. Karakas, A. I., & Lattanzio, J. C. 2014, PASA, 31, e030 [NASA ADS] [CrossRef] [Google Scholar]
  26. Kepley, A. A., Morrison, H. L., Helmi, A., et al. 2007, AJ, 134, 1579 [NASA ADS] [CrossRef] [Google Scholar]
  27. Kobayashi, C., Karakas, A. I., & Lugaro, M. 2020, ApJ, 900, 179 [Google Scholar]
  28. Koppelman, H., Helmi, A., & Veljanoski, J. 2018, ApJ, 860, L11 [NASA ADS] [CrossRef] [Google Scholar]
  29. Koppelman, H. H., Helmi, A., Massari, D., Roelenga, S., & Bastian, U. 2019, A&A, 625, A5 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  30. Limberg, G., Santucci, R. M., Rossi, S., et al. 2021, ApJ, 913, L28 [NASA ADS] [CrossRef] [Google Scholar]
  31. Lind, K., Asplund, M., Barklem, P. S., & Belyaev, A. K. 2011, A&A, 528, A103 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  32. Lövdal, S. S., Ruiz-Lara, T., Koppelman, H. H., et al. 2022, A&A, in press, https://doi.org/10.1051/0004-6361/202243060 [Google Scholar]
  33. Matsuno, T., Koppelman, H. H., Helmi, A., et al. 2022, A&A, 661, A103 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  34. McMillan, P. J. 2017, MNRAS, 465, 76 [Google Scholar]
  35. Myeong, G. C., Evans, N. W., Belokurov, V., Amorisco, N. C., & Koposov, S. E. 2018, MNRAS, 475, 1537 [NASA ADS] [CrossRef] [Google Scholar]
  36. Nidever, D. L., Hasselquist, S., Hayes, C. R., et al. 2020, ApJ, 895, 88 [Google Scholar]
  37. Nissen, P. E., & Schuster, W. J. 2010, A&A, 511, L10 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  38. Nissen, P. E., & Schuster, W. J. 2011, A&A, 530, A15 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  39. Nissen, P. E., Silva-Cabrera, J. S., & Schuster, W. J. 2021, A&A, 651, A57 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  40. Noguchi, K., Aoki, W., Kawanomoto, S., et al. 2002, PASJ, 54, 855 [NASA ADS] [CrossRef] [Google Scholar]
  41. Prantzos, N., Abia, C., Limongi, M., Chieffi, A., & Cristallo, S. 2018, MNRAS, 476, 3432 [Google Scholar]
  42. Reggiani, H., Meléndez, J., Kobayashi, C., Karakas, A., & Placco, V. 2017, A&A, 608, A46 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  43. Reid, M. J., & Brunthaler, A. 2004, ApJ, 616, 872 [Google Scholar]
  44. Roederer, I. U., Sneden, C., Thompson, I. B., Preston, G. W., & Shectman, S. A. 2010, ApJ, 711, 573 [NASA ADS] [CrossRef] [Google Scholar]
  45. Ruiz-Lara, T., Matsuno, T., Sofie Lövdal, S., et al. 2022, A&A, in press, https://doi.org/10.1051/0004-6361/202243061 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  46. Ryan, S. G., Kajino, T., Beers, T. C., et al. 2001, ApJ, 549, 55 [NASA ADS] [CrossRef] [Google Scholar]
  47. Schlegel, D. J., Finkbeiner, D. P., & Davis, M. 1998, ApJ, 500, 525 [Google Scholar]
  48. Schönrich, R., Binney, J., & Dehnen, W. 2010, MNRAS, 403, 1829 [Google Scholar]
  49. Siegel, D. M., Barnes, J., & Metzger, B. D. 2019, Nature, 569, 241 [Google Scholar]
  50. Skúladóttir, Á., Tolstoy, E., Salvadori, S., Hill, V., & Pettini, M. 2017, A&A, 606, A71 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  51. Sneden, C. 1973, ApJ, 184, 839 [Google Scholar]
  52. Suda, T., Katsuta, Y., Yamada, S., et al. 2008, PASJ, 60, 1159 [NASA ADS] [Google Scholar]
  53. Suda, T., Yamada, S., Katsuta, Y., et al. 2011, MNRAS, 412, 843 [NASA ADS] [Google Scholar]
  54. Tajitsu, A., Aoki, W., & Yamamuro, T. 2012, PASJ, 64, 77 [NASA ADS] [CrossRef] [Google Scholar]
  55. Tarumi, Y., Suda, T., van de Voort, F., et al. 2021, MNRAS, 505, 3755 [NASA ADS] [CrossRef] [Google Scholar]
  56. Tsujimoto, T., Matsuno, T., Aoki, W., Ishigaki, M. N., & Shigeyama, T. 2017, ApJ, 850, L12 [Google Scholar]
  57. Wanajo, S., Janka, H.-T., & Müller, B. 2011, ApJ, 726, L15 [NASA ADS] [CrossRef] [Google Scholar]
  58. Wanajo, S., Sekiguchi, Y., Nishimura, N., et al. 2014, ApJ, 789, L39 [Google Scholar]
  59. Watson, D., Hansen, C. J., Selsing, J., et al. 2019, Nature, 574, 497 [NASA ADS] [CrossRef] [Google Scholar]
  60. Weidner, C., & Kroupa, P. 2005, ApJ, 625, 754 [NASA ADS] [CrossRef] [Google Scholar]
  61. Winteler, C., Käppeli, R., Perego, A., et al. 2012, ApJ, 750, L22 [Google Scholar]
  62. Yamada, S., Suda, T., Komiya, Y., Aoki, W., & Fujimoto, M. Y. 2013, MNRAS, 436, 1362 [Google Scholar]
  63. Yan, H. L., Shi, J. R., Nissen, P. E., & Zhao, G. 2016, A&A, 585, A102 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  64. Zhang, L., Ishigaki, M., Aoki, W., Zhao, G., & Chiba, M. 2009, ApJ, 706, 1095 [Google Scholar]

Appendix A: Additional table

Table A.1.

Abundances of the Helmi stream stars

All Tables

Table 1.

Summary of the data.

Table 2.

Property of the targets.

Table 3.

Kinematics of the targets.

Table 4.

Stellar parameters.

Table 5.

Linelist and line-by-line abundance.

Table A.1.

Abundances of the Helmi stream stars

All Figures

thumbnail Fig. 1.

Distribution of stars in kinematic spaces, namely E − Lz, Lperp − Lz, and vy − vz. The insets in the left two panels show better distribution of the stars that are part of the Helmi streams. The Helmi stream selection from Koppelman et al. (2019) is shown as a dashed box in the inset of the middle panel. Stars are hatched according to the subgroups found by Dodd et al. (2022). Three stars (2447_5952, 4998_5552, and LP894-3) are not associated with either of the subgroups. An open symbol is assigned only to 2447_5952 since it has much larger energy than the other Helmi stream stars. We used filled symbols for stars from NS10 and R17 if they satisfy the kinematic selection for Gaia-Enceladus (see Paper I).

In the text
thumbnail Fig. 2.

Location of the stars in the color-magnitude diagram. The symbols follow those of Fig. 1. The four isochrones are from the PARSEC tracks with the age of 12 Gyr and [Fe/H]= − 2.0,   − 1.5,   − 1.0,  and  − 0.5.

In the text
thumbnail Fig. 3.

[Na/Fe] and [Mg/Fe] abundance ratios of the Helmi stream stars (red symbols). Also plotted are halo stars from NS10 and R17, with Gaia-Enceladus stars shown with filled symbols. G112-43 and G112-44 from NS10 are shown with red stars since their kinematics are consistent with the Helmi streams. We also show kinematically selected Sequoia stars from Paper I.

In the text
thumbnail Fig. 4.

α-element abundances of the Helmi stream stars. The symbols follow those of Fig. 3. When a star does not have any detectable lines for an element, we indicate the abundances corresponding to the 1-σ and 3-σ upper limits on the equivalent widths, respectively shown with the location of the symbols and the upper end of the error bars.

In the text
thumbnail Fig. 5.

Abundances of Cr, Mn, Ni, and Zn. Symbols follow Fig. 3.

In the text
thumbnail Fig. 6.

Abundances of Y and Ba. Symbols follow Fig. 3.

In the text
thumbnail Fig. 7.

Chemical abundance of the Helmi stream stars from literature, namely Roederer et al. (R10; 2010), Aguado et al. (A21; 2021), and Gull et al. (G21; 2021). We also include Helmi stream stars in the present study and in NS10. The background gray points are from the SAGA database (Suda et al. 2008, 2011; Yamada et al. 2013). The UVES spectra of three stars from Aguado et al. (2021) were reanalysed and abundances from the reanalysis are included in the figure. Abundances from our reanalysis and those from Aguado et al. (2021) are connected with a solid line for each star.

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.