Issue |
A&A
Volume 638, June 2020
|
|
---|---|---|
Article Number | A154 | |
Number of page(s) | 11 | |
Section | Galactic structure, stellar clusters and populations | |
DOI | https://doi.org/10.1051/0004-6361/202037923 | |
Published online | 29 June 2020 |
The HR 1614 moving group is not a dissolving cluster
1
Lund Observatory, Department of Astronomy and Theoretical Physics, Box 43, 221 00 Lund, Sweden
e-mail: iryna@astro.lu.se, tbensby@astro.lu.se
2
Research School of Astronomy & Astrophysics, Mount Stromlo Observatory, The Australian National University, Canberra, ACT 2611, Australia
Received:
11
March
2020
Accepted:
8
May
2020
Context. The HR 1614 is an overdensity in velocity space and has for a long time been known as an old (∼2 Gyr) and metal-rich ([Fe/H] ≈ +0.2) nearby moving group that has a dissolving open cluster origin. The existence of such old and metal-rich groups in the solar vicinity is quite unexpected since the vast majority of nearby moving groups are known to be young.
Aims. In the light of new and significantly larger data sets than ever before (astrometric, photometric, and spectroscopic), we aim to re-investigate the properties and origin of the HR 1614 moving group. If the HR 1614 overdensity is a dissolving cluster, its stars should represent a single-age and single-elemental abundance population.
Methods. To identify and characterise the HR 1614 moving group we use astrometric data from Gaia DR2; distances, extinction, and reddening corrections from the StarHorse code; elemental abundances from the GALAH and APOGEE spectroscopic surveys; and photometric metallicities from the SkyMapper survey. Bayesian ages were estimated for the SkyMapper stars. Since the Hercules stream is the closest kinematical structure to the HR 1614 moving group in velocity space and as its origin is believed to be well-understood, we use the Hercules stream for comparison purposes. Stars that are likely to be members of the two groups were selected based on their space velocities.
Results. The HR 1614 moving group is located mainly at negative U velocities, does not form an arch of constant energy in the U − V space, and is tilted in V. We find that the HR 1614 overdensity is not chemically homogeneous, but that its stars exist at a wide range of metallicities, ages, and elemental abundance ratios. They are essentially similar to what is observed in the Galactic thin and thick discs, a younger population (around 3 Gyr) that is metal-rich (−0.2 ≤ [Fe/H] ≤ 0.4) and alpha-poor. These findings are very similar to what is seen for the Hercules stream, which is believed to have a dynamical origin and consists of regular stars from the Galactic discs.
Conclusions. The HR 1614 overdensity has a wide spread in metallicity, [Mg/Fe], and age distributions resembling the general properties of the Galactic disc. It should therefore not be considered a dissolving open cluster, or an accreted population. Based on the kinematic and chemical properties of the HR 1614 overdensity we suggest that it has a complex origin that could be explained by combining several different mechanisms such as resonances with the Galactic bar and spiral structure, phase mixing of dissolving spiral structure, and phase mixing due to an external perturbation.
Key words: stars: kinematics and dynamics / Galaxy: kinematics and dynamics / Galaxy: formation / Galaxy: evolution
© ESO 2020
1. Introduction
The Milky Way has a complex structure (e.g. Bland-Hawthorn & Gerhard 2016). By studying how stars move in the Galaxy we can trace the formation history of the Milky Way (e.g. Freeman & Bland-Hawthorn 2002). Stars that share a common motion are usually called kinematic structures, and studying their origin is one of the ways of obtaining more information about the formation and evolution of the Galactic disc. The analysis of data from the HIPPARCOS mission (Perryman et al. 1997) revealed a rich structure of the local velocity distribution (e.g. Dehnen 2000; Arifyanto & Fuchs 2006). Later, an even more complex picture was discovered by the Gaia mission (Gaia Collaboration 2016). The local velocity field is composed of dozens of kinematic structures that together form arches and ridges in different velocity projections (e.g. Antoja et al. 2018; Gaia Collaboration 2018c; Ramos et al. 2018; Khanna et al. 2019; Kushniruk & Bensby 2019). The origin of the non-smooth local velocity distribution is directly linked to major formation processes of the Milky Way. For example, kinematic structures can have different origins, including resonances with the Galactic bar (e.g. Dehnen 2000; Monari et al. 2018), resonances with the spiral arms (e.g. Quillen et al. 2018a), phase mixing as a result of a dynamical interaction with a merging dwarf galaxy (e.g. Minchev et al. 2009; Antoja et al. 2018; Laporte et al. 2019), phase mixing as a result of transient spiral arms and a perturbation with a dwarf satellite (e.g. Khanna et al. 2019), and phase mixing due to transient spiral structure (e.g. Hunt et al. 2019).
Another explanation for some kinematic structures is that open clusters dissolve with time and form moving groups. This idea was first introduced by Olin Eggen (Eggen 1965). Since stars in a moving group originate from the same cluster, they share similar chemical composition, ages, and motion. Two well-known moving groups are the Pleiades and the Hyades, which are visible to the naked eye. More contended examples of moving groups are the Arcturus and the HR 1614 moving groups that were discovered in Eggen (1971, 1978). The groups were later confirmed with HIPPARCOS and Gaia data (e.g. Dehnen 1998; Feltzing & Holmberg 2000; Ramos et al. 2018; Kushniruk & Bensby 2019). For a long time the Arcturus structure was considered a dissolving open cluster, later an accreted stellar population. However, Kushniruk & Bensby (2019) used the Gaia Data Release 2 (DR2) data and found no signatures of it being a dissolving cluster or an accreted population. Instead, the group was likely caused by phase mixing induced by a merger, as was proposed by Minchev et al. (2009).
A discovery of a group of stars with mean radial velocity component U ≃ 0 km s−1 and mean rotational velocity component V ≃ −60 km s−1 and higher than solar metallicity was first reported in Eggen (1978). The group was named after the star HR 16141 which is one of its member stars. It was proposed that the overdensity is a dissolving old open cluster. Later Eggen (1992) estimated that the HR 1614 moving group is about 5 Gyr old. The origin of the group was re-investigated in Feltzing & Holmberg (2000) using HIPPARCOS data. Candidate member stars of the group were carefully selected by examining Hertzsprung–Russell diagrams (HR diagram) for different slices of the U − V space velocity distribution. A slice with younger and more metal-rich stars was found. From dynamical simulations of a disrupting open cluster they could explain the tilt of the group in the U − V space. They also found the group to be about 2 Gyr old, and that its metallicity is [Fe/H] ≃ 0.19. Later, De Silva et al. (2007) performed a spectroscopic study of stars that were assigned as members of the group by Feltzing & Holmberg (2000) and confirmed that HR 1614 is a 2 Gyr metal-rich group with [Fe/H] ≥ 0.25. It was also found to be chemically homogeneous with a scatter of only about 0.01 dex in various chemical elements. So, from the literature we know that HR 1614 is a metal-rich, 2 Gyr moving group that is scattered around the Sun.
When searching for the Arcturus stream Kushniruk & Bensby (2019) were able to clearly identify the HR 1614 moving group in the velocity distribution. The question is whether HR 1614 can still be considered a dissolved open cluster if the much larger and more precise astrometric sample from Gaia DR2 is used to identify and characterise its properties in more detail. In this paper we aim to do just that, to analyse kinematic, chemical, and photometric properties of the 1614 moving group in order to constrain its origin using data from Gaia DR2 and spectroscopic surveys such as APOGEE and GALAH.
The paper is structured in the following way. The data sets used are described in Sect. 2. The selection of candidate member stars of the HR 1614 moving group is described in Sect. 3, where we also analyse the HR diagrams, and the metallicity and age distributions of the HR 1614 moving group stars. In Sect. 4 we discuss the possible origins of HR 1614, and finally our findings are summarised in Sect. 5.
2. Data
In this work we use proper motions, sky positions, and radial velocities from the Gaia DR2 (Gaia Collaboration 2018b) catalogue, Bayesian extinction corrections and distances from the StarHorse code (Anders et al. 2019). We cross-match the StarHorse and Gaia DR2 catalogues by Gaia IDs where the radial velocity is not NULL, where the internal StarHorse quality flag SH_OUTFLAG is set to ‘00000’, and where the Gaia quality flag SH_GAIAFLAG is set to ‘000’. These are recommended quality flags that are explained in detail in Anders et al. (2019). In addition, we also cut out stars with parallax_over_error > 10 and visibility_periods_used > 8 as suggested by Gaia Collaboration (2018a). The final data query is as follows:
SELECT s.*, g.* FROM gdr2.gaia_source AS g, gdr2_contrib.starhorse AS s WHERE g.source_id = s.source_id AND g.radial_velocity IS NOT NULL AND s.SH_OUTFLAG LIKE `00000' AND s.SH_GAIAFLAG LIKE `000' AND g.parallax_over_error > 10 AND g.visibility_periods_used > 8
The query listed above gives us a sample of 4 790 725 stars. We use the galpy2 package (Bovy 2015) to calculate: 1) space velocities U, V, and W3 corrected for peculiar motion of the Sun with (U⊙, V⊙, W⊙) = (11.10, 12.24, 7.25) km s−1 (Schönrich et al. 2010) and 2) Galactocentric cylindrical coordinates R, ϕ, and Z4 with R⊙ = 8.34 kpc (Reid et al. 2014), ϕ⊙ = 0°, and Z⊙ = 14 pc (Binney et al. 1997).
To characterise the groups chemically, we use elemental abundances from both the APOGEE (Majewski et al. 2017) and GALAH (De Silva et al. 2015) spectroscopic surveys. APOGEE provides stellar atmospheric parameters and elemental abundances derived from near-IR (H-band), high-resolution (R ∼ 23 000) spectra. We use the abundances from APOGEE DR16 (Ahumada et al. 2019), which includes stars in both the northern and southern hemispheres. After applying the same cuts for quality assurance as in Feuillet et al. (2019) and cross-matching our sample with APOGEE by Gaia IDs, we have a sample of 126 690 red giants. GALAH is also a high-resolution (R ∼ 28 000) spectroscopic survey operating at optical wavelengths in the southern hemisphere that provides stellar atmospheric parameters and elemental abundances. We use GALAH DR2 Buder et al. (2018), setting the flag_cannon and flag_x_fe flags to zero as recommended. Cross-matching GALAH with our set of stars results in a sample of 231 725 stars. We also use data from the SkyMapper (Casagrande et al. 2019) survey, which provided photometric metallicities for 907 893 stars after cross-matching it with our sample of Gaia IDs. To navigate on HR diagrams presented in this paper, we use PARSEC 1.2S isochrones5 re-derived for the Gaia DR2 photometric system (Maíz Apellániz & Weiler 2018).
3. Analysis
3.1. The old HR 1614 moving group
The nature of the HR 1614 overdensity was previously investigated by Feltzing & Holmberg (2000). They used HIPPARCOS data and divided the U − V space into seven smaller boxes. Individual HR diagrams for each box revealed one box with a population that was younger and more metal-rich than the Sun, which was connected to the HR 1614 moving group.
In this study we repeat the same procedure, but with our much larger stellar sample. Figure 1 shows the U − V distribution of 581 190 stars close to the Sun within a region defined by R⊙ ± 0.2 kpc and ϕ⊙ ± 1.5°. This data is sliced into seven boxes, as was done in Feltzing & Holmberg (2000). The boxes are shown with dashed lines in Fig. 1. Box 5 is where the group is supposed to be located according to Feltzing & Holmberg (2000). With Gaia DR2/StarHorse we clearly see that box 5 contains two velocity overdensities, whereas in Feltzing & Holmberg (2000) two separate overdensities were not observed (probably due to a significantly smaller stellar sample). The inclined overdensity in box 5 in Fig. 1 is what we link to the HR 1614 moving group (marked by a black ellipse). The horizontal overdensity in box 5 is known as the Hercules stream (marked by a red ellipse). The identification of the groups is based on results of the wavelet transform performed in Kushniruk & Bensby (2019).
Fig. 1. U − V distribution of 581 190 stars inside a region defined in R⊙ ± 0.2 kpc and ϕ⊙ ± 1.5° (see Sect. 3). The color scale is proportional to the number of stars, as indicated in the bar at the top. Stars in the HR 1614 moving group from Feltzing & Holmberg (2000) are shown as blue crosses and their best-selected candidates as black open circles. We identify stars within the red and black ellipses as members of the Hercules stream and the HR 1614 moving group, respectively (see discussion in Sect. 3). Dashed lines and numbers show the same boxes as in Feltzing & Holmberg (2000). |
The Hercules stream has been studied in many works: it is likely caused by resonances with the bar (e.g. Dehnen 2000; Bensby et al. 2007; Wegg et al. 2015), and thus, will be used as a benchmark group in this study. Due to its resonant origin, the Hercules stream mainly consists of stars from the Galactic disk. Fig. 1 show stars that are potential members of the HR 1614 group taken from Table 1 in Feltzing & Holmberg (2000). These are the same stars as the ones shown as filled black circles in Table 1 in Feltzing & Holmberg (2000). Interestingly, although Feltzing & Holmberg (2000) did not see HR 1614 as a separate overdensity in their box 5, the circled crosses that show stars from their paper and our black ellipse are tilted at similar angles.
Figure 2 shows individual HR diagrams for boxes 1–7 using extinction corrected colours and absolute magnitudes from the StarHorse catalogue. Two sets of isochrones are shown: one with an age of 2 Gyr and metallicity [M/H] = +0.2 and one with an age of 8 Gyr and [M/H] = +0.0. It is clear that none of the HR diagrams can be fitted with only one isochrone. This probably means that the stars in all boxes represent an underlying population of stars with a wide range in metallicity and/or age. Similarly to what was seen in Feltzing & Holmberg (2000), boxes 2 and 5 appear to contain younger stars compared to the other boxes. These are also the boxes with the largest number of stars. As was shown in Fig. 1, in our larger data set we clearly see that box 5 consists of two different groups, Hercules and HR 1614, that were not observed as separate overdensities with the HIPPARCOS data (Feltzing & Holmberg 2000).
Fig. 2. Hertzsprung–Russel diagrams for stars inside the seven boxes as shown in Fig. 1. The dashed line is an isochrone for [M/H] = 0 and 8 Gyr, the solid line shows an isochrone for [M/H] = + 0.2 and 8 Gyr. The Colour bar shows the number of stars in each bin. |
3.2. The new HR 1614 moving group
A more robust way to select candidate member stars in the HR 1614 moving group is to apply the results of the wavelet transform analysis presented in Kushniruk & Bensby (2019). They analysed the Gaia DR2 radial velocity sample in the U − V space by individually exploring smaller sub-samples at different R and ϕ. Their analysis is based on the discrete wavelet transform, an algorithm that decomposes data into a set of wavelet coefficients. These coefficients contain information about the location and significance of velocity overdensities present in the data. More details on the wavelet transform can be found in Starck & Murtagh (2002), and about the procedure in Kushniruk et al. (2017) and Kushniruk & Bensby (2019).
The locations of the sub-samples (regions 00, 01a, 11a, 21a, and 31a) from Kushniruk & Bensby (2019) in Cartesian Galactic X and Y coordinate system are shown in the top left plot of Fig. 3. The other plots in Fig. 3 show wavelet transform maps for the corresponding regions. The HR 1614 moving group is observed as four overdensities with mean U = −18 km s−1 and V = −65 km s−1 in region 00. The group is also visible in regions 11a and 31a at the same location and is slightly shifted downwards in V in region 01a and shifted upwards in region 21a. The shift occurs due to the location of the regions at different Galacticentric radii.
Fig. 3. Top left: location of five regions (00, 01a, 11a, 21a, and 31a) in the X − Y space explored in this work. Other plots: wavelet transform maps that show location of the kinematic overdensities in the U − V space. Their centres are shown with blue crosses. Lines connect eponymous overdensities into structures. The names of the structures are listed in the legend. The colour bar shows the normalised wavelet coefficients in each bin. |
The U − V distribution of stars in Fig. 1 is shown for region 00. A black ellipse drawn around the HR 1614 moving group is centred at U = −18 km s−1 and V = −65 km s−1, and a red ellipse around the Hercules stream is centred at U = −37 km s−1 and V = −49 km s−1. The positions of the two groups were taken from Kushniruk & Bensby (2019), and the height and width of the ellipse were arbitrarily selected to fit the U − V distribution.
Figure 4 shows HR diagrams for the stars inside the black and red ellipses that we link to the HR 1614 moving group and the Hercules stream. The red circles show stars from Feltzing & Holmberg (2000) that we also could identify in our StarHorse sample. Most of these stars follow an [M/H] = 0.2 and 2 Gyr isochrone, the parameters of HR 1614 that were found by Feltzing & Holmberg (2000) and De Silva et al. (2007). At the same time it is clearly seen that stars from ellipses around Hercules and HR 1614 are also composed of multiple stellar populations and cannot be fitted with only one isochrone. From now on we refer to stars inside the black and red ellipses in Fig. 1 as the HR 1614 moving group and the Hercules stream, respectively.
Fig. 4. Hertzsprung–Russel diagrams for 12 654 stars from the HR 1614 moving group (left) and 36 656 stars from Hercules (right) selected from the black and red ellipses shown in Fig. 1. An isochrone with [M/H] = + 0.2 and 2 Gyr is shown as a black solid line, and one with [M/H] = 0 and 8 Gyr is shown as a dashed line. Red dots in the top left plot show all the stars from Table 1 in Feltzing & Holmberg (2000). The colour bar shows the number of stars in each bin for both plots. |
3.2.1. HR diagrams
Figures 5 and 6 show HR diagrams of the HR 1614 moving group and the Hercules stream, respectively, for six metallicity bins in the range −0.6 < [Fe/H] < 0.6. The top panels in each figure show the APOGEE and GALAH data (red and blue dots, respectively), and the bottom panels show stars with SkyMapper data (green dots). We plot isochrones with a metallicity corresponding to the mean metallicity of the stars in each bin. The HR diagrams show that the HR 1614 moving group and the Hercules stream cover a wide range of metallicities and ages. An age–metallicity gradient appears to be present in both groups. Metal-poor stars tend to be older and metal-rich stars seem to be younger, based on a visual inspection of the HR diagrams. It is worth noting that stars in the most metal-poor and metal-rich bins do not perfectly match the isochrones. Therefore, in Fig. 7 we compare metallicities for stars in common between the GALAH and SkyMapper, and the APOGEE and SkyMapper surveys. The difference between SkyMapper and GALAH or APOGEE metallicities increases for metal-rich SkyMapper targets. As the SkyMapper metallicities are estimated based on photometric parameters and calibrated on data from the GALAH survey. The SkyMapper metallicities are overestimated for targets more metal-rich than ≃0.5 dex, where the GALAH survey does not provide any stars (see Fig. 11 in Casagrande et al. (2019) for more details).
Fig. 5. Hertzsprung–Russel diagrams for the stars in the HR 1614 moving group selected from the cross-matched samples between StarHorse and APOGEE (red), and StarHorse and GALAH (blue), and StarHorse and SkyMapper for different metallicity bins. Black lines show isochrones for a mean metallicity of each bin and cover ages 0.2 Gyr and from 2 to 12 Gyr with the steps of 2 Gyr. |
Fig. 7. Difference between GALAH and SkyMapper metallicities as a function of SkyMapper metallicity (left). Difference between APOGEE and SkyMapper metallicities as a function of SkyMapper metallicity (right). Both plots are colour-coded by [Mg/Fe] values taken from the GALAH and APOGEE samples, respectively, and are shown for stars in the HR 1614 moving group and the Hercules stream. Blue crosses show typical (mean) error. |
Overall, the HR diagrams presented in Figs. 4–6 show that the HR 1614 moving group is a more complex structure than previously thought, and that it does not consist of single-age and single-elemental abundance population of stars, similarly to the Hercules stream. This finding contradicts the hypothesis that the HR 1614 overdensity has a moving group origin, in which case its stars would show homogeneity in age and elemental abundance.
3.2.2. Abundance trends and metallicity distribution
Figure 8 shows the [Mg/Fe] versus [Fe/H] diagrams for stars from the GALAH and APOGEE catalogues that are located in region 00. Candidate members of the HR 1614 moving group (left panels) and the Hercules stream (right panels) are shown as red dots. The GALAH data (lower panels) generally show larger uncertainties than the APOGEE data (upper panels), and therefore we only show a single error bar representing the mean uncertainty for the GALAH stars. Background distributions as well as stars from the groups are composed of two sequences, low- and high-alpha stellar populations, that can be recognised as the chemically defined Galactic thin and thick discs (e.g. Bensby et al. 2011).
Fig. 8. [Mg/Fe]–[Fe/H] diagrams for HR 1614 (left) and Hercules (right) selected from the APOGEE (top) and GALAH (bottom) surveys are shown as red dots. The corresponding uncertainties are shown as error bars for APOGEE distributions and a typical (mean) error is shown for GALAH stars. Background black distributions show stars from APOGEE and GALAH in region 00. |
The metallicity distributions for the HR 1614 moving group and the Hercules stream with GALAH, APOGEE, and SkyMapper data are shown in Fig. 9. The distributions for the HR 1614 moving group look more spiky due to lower number of stars in the group compared to the Hercules stream. Again, both groups cover a wide range of metallicities and are likely composed of stars that come from different stellar populations. We do not observe the overdensity around [Fe/H] = 0.2 in the metallicity distribution of the HR 1614 moving group. These findings question the existence of a single metal-rich population of stars reported in Feltzing & Holmberg (2000) and De Silva et al. (2007) as well as the dissolving cluster origin of the HR 1614 overdensity.
Fig. 9. Metallicity distributions for HR 1614 (left) and Hercules (right) from the GALAH (blue), APOGEE (red), and SkyMapper (green) surveys. Both regular histograms (upper three panels) and normalised cumulative histograms (bottom panel) are shown. |
3.2.3. Age distribution
Figure 10 shows the age–metallicity probability distributions for the HR 1614 moving group and the Hercules streams. To determine these distributions, we first calculated the two-dimensional probability distributions in age and metallicity, 𝒢i(τ, [Fe/H]), for each star i following the method by Howes et al. (2019). In short, the 𝒢 function of a star is calculated by fitting a grid of isochrones (PARSEC) to the observed G magnitude, GBP − GRP colour, metallicity, and distance. Then the age–metallicity distribution of the sample as a whole, ϕ(τ, [Fe/H]), is estimated by maximising the likelihood
Fig. 10. Age–metallicity probability distributions for the HR 1614 moving group (left) and Hercules stream (right) for stars from the SkyMapper survey. The upper panels show the normalised age distributions after summing over the metallicity dimension. |
using an inversion algorithm. The inversion is subject to regularisation which means there is a single free parameter governing the smoothness of ϕ. Figure 10 shows the results for one choice of this parameter, but we have tested a wide range of values and find that the results are similar enough to not affect our conclusions. This method will be described in detail in an upcoming publication (Sahlholdt & Lindegren, in prep.).
Both groups are composed of at least two stellar populations. The first population is metal-poor and old with the centre at [Fe/H] ≈ −0.2 and an age of about 8 Gyr. The second population is metal-rich and young with the centre at [Fe/H] ≈ +0.1 and an age of about 3 Gyr. These two populations can most likely be associated with the Galactic thin and thick discs. There is also a third clump at a very young age, less than about 1 Gyr. These are stars with ages that are located at the edge of the grid and are artefacts of the age estimation procedure.
From this analysis it is clear that the HR 1614 overdensity is not a 2 Gyr old stellar population, but rather is composed of a mix of stars from the Galactic disc. The presence of an age bi-modality in the HR 1614 overdensity also contradicts the idea of its dissolving cluster origin.
4. Origin of the HR 1614 overdensity
4.1. Dissolving open cluster origin
Stars that originate from a dissolving open cluster were formed from the same molecular cloud, and thus, retain similar chemical compositions and ages (Freeman & Bland-Hawthorn 2002). As was proposed in Eggen (1971, 1992, 1998), for a long time HR 1614 was considered a classical moving group. The analysis of HIPPARCOS stars accompanied by photometric and spectroscopic data by Feltzing & Holmberg (2000) and De Silva et al. (2007) showed that HR 1614 was a unique stellar population with a metallicity [Fe/H] ≈ 0.2, and age of about 2 Gyr with a very small scatter in abundances in other elements. However, we do not see any evidence that the HR 1614 overdensity is a mono-age and mono-abundance stellar population in our study. The analysis of HR diagrams, metallicity, and age distributions presented in Sect. 3 clearly shows that HR 1614 is a mix of two stellar populations that resemble properties of the Galactic thin and thick discs. A small scatter in abundances observed in De Silva et al. (2007) is, possibly, a result of the stellar sample analysed in their work. It is mainly composed of stars from Feltzing & Holmberg (2000), who selected metal-rich stars from their Box 5 as possible members of the HR 1614 moving group. This selection effect can potentially explain the small scatter in abundances reported in De Silva et al. (2007).
In this work we selected 12 654 targets in the HR 1614 moving group from our data sample. This is a significantly larger number of stars-members of the group than ever analysed before. Pre-Gaia works studied stars that are located roughly within 100 pc around the Sun and found that the HR 1614 overdensity is a single age and abundance population (e.g. Feltzing & Holmberg 2000; De Silva et al. 2007). Since our stellar sample covers a larger volume around the Sun, we check if young stars with higher than solar metallicities are located closer to the Sun in Fig. 11. We do not observe any peculiarities in distributions shown for SkyMapper stars in Fig. 11. Both distributions look similar to younger stars in the metal-rich part of the diagram and older stars in the metal-poor part. It also shows that most of the stars are located within a 500 pc radius around the Sun. There is no unique, young, and nearby stellar population with higher-than-solar metallicity neither in the HR 1614 nor in the Hercules velocity overdensities. This means that HR 1614 should not be considered a dissolving open cluster any longer.
Fig. 11. Distance relative to the Sun vs. metallicity diagram colour-cored by age for SkyMapper stars in the HR 1614 moving group (left) and the Hercules stream (right). |
4.2. Accretion origin
Stars that originate from a dwarf galaxy or satellite and were accreted to the Milky Way retain unique chemical composition, ages, and velocities that are usually different from the background Milky Way stars (e.g. Tolstoy et al. 2009; Ruchti et al. 2015; Helmi 2020). Recent examples of such populations are Gaia-Sequoia (e.g. Myeong et al. 2019), Gaia-Sausage/Gaia-Enceladus (e.g. Belokurov et al. 2018; Helmi et al. 2018). Based on the metallicity, velocity, and age distributions discussed in Sect. 3, it is clear that the HR 1614 and Hercules overdensities are not accreted stellar populations. The dynamical and chemical properties of the groups are similar to the vast majority of the Milky Way disc stars. It is unlikely that HR 1614 is an accreted stellar population.
4.3. Resonant origin
Stars that were influenced by resonances with the bar and/or spiral arms can be observed as kinematic structures that consist of mixed stellar populations. Hercules is a group that is widely discussed in the context of resonances. As was first proposed in Dehnen (2000), the Hercules stream is possibly caused by the outer Lindblad resonance (OLR) of the Galactic bar. Later it was proposed that the Hercules stream is likely a result of a corotation resonance (CR) with the bar (e.g. Pérez-Villegas et al. 2017; Binney 2020). Recent results from Hunt et al. (2019) show that it is possible to reproduce the Hercules stream in simulations by combining bar and multiple spiral structures. In their Fig. 2 Hunt et al. (2019) call the bottom-most structure a part of the Hercules stream. If we convert V velocity of the HR 1614 group to the vertical component of the angular momentum Lz, assuming that the Sun is located at 8.34 kpc from the Galactic centre, we get ∼1340 kpc km s−1. This is similar to the value found in Hunt et al. (2019). Hunt et al. (2019) and we use the same values that we do for the peculiar motion of the Sun. It is exactly the same position in the distribution where the HR 1614 overdensity is supposed to be located, and thus, should be called the HR 1614 moving group instead.
Many studies show that the Hercules stream conserves a vertical component of angular momentum with Galactocentric radii which is a signature of a resonant origin (e.g. Ramos et al. 2018; Hunt et al. 2019; Kushniruk & Bensby 2019). In Figs. 1 and 3 we clearly see that the Hercules stream is not tilted in the U − V space. It covers a wide range of negative U velocities, but its V velocity is −50 km s−1 and is almost constant. This means that the vertical component of angular momentum Lz is conserved (taking into account that V is proportional to Lz). Unlike the Hercules stream, the HR 1614 overdensity is tilted in the U − V space. This tilt was explained in Feltzing & Holmberg (2000) in dynamical simulations of a dissolving open cluster. At the same time Feltzing & Holmberg (2000) found that the HR 1614 overdensity is homogeneous in age and metallicity, which we do not observe in our work. The HR 1614 overdensity is also present in other four regions shown in Fig. 3. It almost retain Lz with R. Taking into account that the HR 1614 overdensity is mainly located at negative U and its proximity to the Hercules stream, the resonant origin of the HR 1614 overdensity should not be excluded.
4.4. Phase mixing origin
The concept of phase mixing implies that a phase-space (position and momentum variables) full of overdensities will evolve with time to a stationary state (Tremaine 1999). A moving groups is one example of a phase mixing process. In the recent literature, phase mixing is mainly discussed in the context of internal or external perturbation mechanisms that can induce the process on large scales (e.g. Laporte et al. 2019; Hunt et al. 2019). This discussion is especially pertinent after the discovery of a phase spiral by Antoja et al. (2018), which is strong evidence of an ongoing phase mixing in the Galactic disc.
Symmetric arches of constant energy around U = 0 km s−1 with 20 km s−1 separation in V were predicted in the model that was first proposed by Minchev et al. (2009), where a perturbation with a dwarf galaxy, that took place roughly 2 Gyr ago, leads to a phase wrapping. This external perturbation event could have contributed to a velocity distribution and can be observed as overdensities aligned across lines of constant energy. Minchev et al. (2009) discuss the HR 1614 overdensity in the context of a dissolving cluster and try to fit it into their ringing theory. They state that the HR 1614 moving group is elongated in U and curved in V in their Figure 3, which makes is consistent with a wrapping origin. The fact that the moving group is not symmetric and is elongated towards negative U is explained by Minchev et al. (2009) as being due to the proximity to the Hercules stream, which is most likely caused by the Galactic bar. In summary, Minchev et al. (2009) suggest that the HR 1614 overdensity is a dissolving open cluster that formed earlier than an external perturbation event took place and was distorted by phase wrapping, forming almost an arch of constant energy, and by resonances with the bar, causing elongation towards negative U.
Since the arches predicted by Minchev et al. (2009) are symmetrical and the observed velocity distribution has non-asymmetries, Quillen et al. (2018b) proposed another explanation, where asymmetric arches are formed due to stars crossing spiral arms. At the same time their theory cannot explain vertical velocity distribution, and thus external perturbation is not excluded. Hunt et al. (2019) explored various combinations of spirals and bar resonances, and stars crossing dissolving transient spirals and found that there are many models that can successfully explain the main velocity overdensities. Some simulations reproduce a velocity overdensity at the place where the HR 1614 moving group is supposed to be. Khanna et al. (2019) showed that phase mixing of disrupting spiral arms and phase mixing due to external perturbation can generate arches and ridges, but they do not talk specifically about the HR 1614 moving group in their work.
In our work we clearly see that the HR 1614 overdensity is not a dissolving open cluster due to chemical and age spread in the group of stars. We also observe that the HR 1614 overdensity is located mainly at negative U. At the same time the group is tilted in V, but does not form a complete arch of constant energy. This points towards a complex origin of the HR 1614 overdensity. The stars could be influenced by transiting spirals and bar that would shift it towards negative U. Phase mixing due to disrupting spiral structure or ringing due to external perturbation could have contributed to the tilt in V.
5. Conclusions
In this paper we revised the origin of the HR 1614 moving group. We analysed a combination of Gaia and StarHorse catalogues and investigated kinematic and photometric properties of stars in a small volume near to the Sun. To study chemical properties of the group we accompanied kinematic data with elemental abundances from the APOGEE and GALAH surveys and photometric metallicities SkyMapper survey. For SkyMapper stars Bayesian ages were calculated for stellar populations. Combining these surveys allowed us to analyse a significantly larger data set with higher precision in the astrometric, photometric, and spectroscopic measurements compared to previous studies. Using the methods developed in Kushniruk & Bensby (2019) we have improved the selection of potential members. These two advances allowed us to better explore the possibility that HR 1614 is or is not a dissolving open cluster. We also compared the properties of the HR 1614 moving group with the Hercules stream which is the nearest velocity overdensity to the HR 1614 moving group. The Hercules stream is most likely of resonant origin, making it a valuable comparison to determine whether HR 1614 is a single stellar population. The main results of the paper are the following:
-
The HR 1614 overdensity is clearly present in the U − V velocity distribution at U ≃ −20, V ≃ −60 km s−1. This location is consistent with the results from previous studies (e.g. Eggen 1998; Feltzing & Holmberg 2000; Ramos et al. 2018; Kushniruk & Bensby 2019).
-
The HR diagrams, metallicity, and age distributions show that HR 1614 consists of two stellar populations that have properties similar to the thin and thick discs.
-
The HR 1614 overdensity is mainly present at negative U velocities, which could point towards a resonant origin similarly to the Hercules stream. At the same time the HR 1614 overdensity is slightly tilted in V, which means that the group has a variation in its vertical angular momentum distribution.
-
The HR 1614 overdensity does not form a complete arch of constant energy in the U − V space, which would be an indication of a phase mixing origin.
Based on the above, it is clear that the HR 1614 overdensity is neither a dissolving open cluster nor an accreted stellar population. We conclude that the HR 1614 overdensity has a complex origin which is a combination of various dynamical mechanisms. As in the case of the Hercules stream, the HR 1614 overdensity could have formed due to different types of resonances such as the CR or OLR with the Galactic bar. This together with phase mixing due to disrupting spiral arms could explain the shape and location of the group in the U − V space. Phase mixing due to external perturbation with a dwarf galaxy or a satellite that happened 2 Gyr ago, as proposed in Minchev et al. (2009), is also possible. The mixed populations that we observe in the HR 1614 overdensity and the Hercules stream are older than a potential perturbation event. Also in Kushniruk & Bensby (2019) we observed kinematic structures with roughly 20 km s−1 separation in V, especially at lower V, which is consistent with a ringing event proposed in Minchev et al. (2009).
In order to disentangle how many different dynamical mechanisms have actually contributed to the formation of velocity structures like the HR 1614 overdensity, further studies are required. Numerical simulations combined with detailed investigation of elemental abundances from spectroscopic surveys like WEAVE (Dalton et al. 2014), 4MOST (de Jong et al. 2019), and the Gaia-ESO survey (Gilmore et al. 2012), and kinematics from upcoming Gaia data releases will provide more information about the origin of kinematic structures.
Available at https://github.com/jobovy/galpy
Available at http://stev.oapd.inaf.it/cgi-bin/cmd
Acknowledgments
T.B. was funded by grant No. 2018-04857 from the Swedish Research Council. L.C. acknowledges support from the Australian Research Council Future Fellowship FT160100402. C.L.S., D.F., and S.F. were supported by the project grant The New Milky Way from the Knut and Alice Wallenberg foundation and by the grant 2016-03412 from the Swedish Research Council.
References
- Ahumada, R., Allende Prieto, C., Almeida, A., et al. 2019, ArXiv e-prints [arXiv:1912.02905] [Google Scholar]
- Anders, F., Khalatyan, A., Chiappini, C., et al. 2019, A&A, 628, A94 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Antoja, T., Helmi, A., Romero-Gómez, M., et al. 2018, Nature, 561, 360 [NASA ADS] [CrossRef] [Google Scholar]
- Arifyanto, M. I., & Fuchs, B. 2006, A&A, 449, 533 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Belokurov, V., Erkal, D., Evans, N. W., Koposov, S. E., & Deason, A. J. 2018, MNRAS, 478, 611 [NASA ADS] [CrossRef] [Google Scholar]
- Bensby, T., Oey, M. S., Feltzing, S., & Gustafsson, B. 2007, ApJ, 655, L89 [NASA ADS] [CrossRef] [Google Scholar]
- Bensby, T., Alves-Brito, A., Oey, M. S., Yong, D., & Meléndez, J. 2011, ApJ, 735, L46 [NASA ADS] [CrossRef] [Google Scholar]
- Binney, J. 2020, MNRAS, 495, 895 [CrossRef] [Google Scholar]
- Binney, J., Gerhard, O., & Spergel, D. 1997, MNRAS, 288, 365 [NASA ADS] [CrossRef] [Google Scholar]
- Bland-Hawthorn, J., & Gerhard, O. 2016, Ann. Rev. Astron. Astrophys., 54, 529 [Google Scholar]
- Bovy, J. 2015, ApJS, 216, 29 [NASA ADS] [CrossRef] [Google Scholar]
- Buder, S., Asplund, M., Duong, L., et al. 2018, MNRAS, 478, 4513 [NASA ADS] [CrossRef] [Google Scholar]
- Casagrande, L., Wolf, C., Mackey, A. D., et al. 2019, MNRAS, 482, 2770 [NASA ADS] [Google Scholar]
- Dalton, G., Trager, S., Abrams, D. C., et al. 2014, in SPIE Conf. Ser., 9147, 91470L [Google Scholar]
- de Jong, R. S., Agertz, O., Berbel, A. A., et al. 2019, Messenger, 175, 3 [Google Scholar]
- De Silva, G. M., Freeman, K. C., Bland-Hawthorn, J., et al. 2015, MNRAS, 449, 2604 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
- De Silva, G. M., Freeman, K. C., Bland-Hawthorn, J., Asplund, M., & Bessell, M. S. 2007, AJ, 133, 694 [NASA ADS] [CrossRef] [Google Scholar]
- Dehnen, W. 1998, AJ, 115, 2384 [NASA ADS] [CrossRef] [Google Scholar]
- Dehnen, W. 2000, AJ, 119, 800 [NASA ADS] [CrossRef] [Google Scholar]
- Eggen, O. J. 1965, Moving Groups of Stars (University of Chicago Press), 111 [Google Scholar]
- Eggen, O. J. 1971, PASP, 83, 762 [NASA ADS] [CrossRef] [Google Scholar]
- Eggen, O. J. 1978, ApJ, 222, 203 [NASA ADS] [CrossRef] [Google Scholar]
- Eggen, O. J. 1992, AJ, 104, 1493 [NASA ADS] [CrossRef] [Google Scholar]
- Eggen, O. J. 1998, AJ, 115, 2397 [NASA ADS] [CrossRef] [Google Scholar]
- Feltzing, S., & Holmberg, J. 2000, A&A, 357, 153 [NASA ADS] [Google Scholar]
- Feuillet, D. K., Frankel, N., Lind, K., et al. 2019, MNRAS, 489, 1742 [NASA ADS] [CrossRef] [Google Scholar]
- Freeman, K., & Bland-Hawthorn, J. 2002, ARA&A, 40, 487 [NASA ADS] [CrossRef] [Google Scholar]
- Gaia Collaboration (Prusti, T., et al.) 2016, A&A, 595, A1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Gaia Collaboration (Katz, D., et al.) 2018a, A&A, 616, A11 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Gaia Collaboration (Brown, A. G. A., et al.) 2018b, A&A, 616, A1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Gaia Collaboration (Babusiaux, C., et al.) 2018c, A&A, 616, A10 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Gilmore, G., Randich, S., Asplund, M., et al. 2012, Messenger, 147, 25 [Google Scholar]
- Helmi, A. 2020, ARA&A, submitted [Google Scholar]
- Helmi, A., Babusiaux, C., Koppelman, H. H., et al. 2018, Nature, 563, 85 [NASA ADS] [CrossRef] [Google Scholar]
- Howes, L. M., Lindegren, L., Feltzing, S., Church, R. P., & Bensby, T. 2019, A&A, 622, A27 [CrossRef] [EDP Sciences] [Google Scholar]
- Hunt, J. A. S., Bub, M. W., Bovy, J., et al. 2019, MNRAS, 490, 1026 [NASA ADS] [CrossRef] [Google Scholar]
- Khanna, S., Sharma, S., Tepper-Garcia, T., et al. 2019, MNRAS, 489, 4962 [CrossRef] [Google Scholar]
- Kushniruk, I., & Bensby, T. 2019, A&A, 631, A47 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Kushniruk, I., Schirmer, T., & Bensby, T. 2017, A&A, 608, A73 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Laporte, C. F. P., Minchev, I., Johnston, K. V., & Gómez, F. A. 2019, MNRAS, 485, 3134 [Google Scholar]
- Maíz Apellániz, J., & Weiler, M. 2018, A&A, 619, A180 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Majewski, S. R., Schiavon, R. P., Frinchaboy, P. M., et al. 2017, AJ, 154, 94 [NASA ADS] [CrossRef] [Google Scholar]
- Minchev, I., Quillen, A. C., Williams, M., et al. 2009, MNRAS, 396, L56 [NASA ADS] [CrossRef] [Google Scholar]
- Monari, G., Famaey, B., Minchev, I., et al. 2018, Res. Notes Am. Astron. Soc., 2, 32 [NASA ADS] [CrossRef] [Google Scholar]
- Myeong, G. C., Vasiliev, E., Iorio, G., Evans, N. W., & Belokurov, V. 2019, MNRAS, 488, 1235 [NASA ADS] [CrossRef] [Google Scholar]
- Pérez-Villegas, A., Portail, M., Wegg, C., & Gerhard, O. 2017, ApJ, 840, L2 [NASA ADS] [CrossRef] [Google Scholar]
- Perryman, M. A. C., Lindegren, L., Kovalevsky, J., et al. 1997, A&A, 500, 501 [NASA ADS] [Google Scholar]
- Quillen, A. C., Carrillo, I., Anders, F., et al. 2018a, MNRAS, 480, 3132 [NASA ADS] [CrossRef] [Google Scholar]
- Quillen, A. C., De Silva, G., Sharma, S., et al. 2018b, MNRAS, 478, 228 [NASA ADS] [CrossRef] [Google Scholar]
- Ramos, P., Antoja, T., & Figueras, F. 2018, A&A, 619, A72 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Reid, M. J., Menten, K. M., Brunthaler, A., et al. 2014, ApJ, 783, 130 [Google Scholar]
- Ruchti, G. R., Read, J. I., Feltzing, S., et al. 2015, MNRAS, 450, 2874 [NASA ADS] [CrossRef] [Google Scholar]
- Schönrich, R., Binney, J., & Dehnen, W. 2010, MNRAS, 403, 1829 [NASA ADS] [CrossRef] [Google Scholar]
- Starck, J.-L., & Murtagh, F. 2002, Astronomical Image and Data Analysis, (Berlin; New York: Springer), oCLC: 679368657 [Google Scholar]
- Tolstoy, E., Hill, V., & Tosi, M. 2009, ARA&A, 47, 371 [NASA ADS] [CrossRef] [Google Scholar]
- Tremaine, S. 1999, MNRAS, 307, 877 [NASA ADS] [CrossRef] [Google Scholar]
- Wegg, C., Gerhard, O., & Portail, M. 2015, MNRAS, 450, 4050 [NASA ADS] [CrossRef] [Google Scholar]
All Figures
Fig. 1. U − V distribution of 581 190 stars inside a region defined in R⊙ ± 0.2 kpc and ϕ⊙ ± 1.5° (see Sect. 3). The color scale is proportional to the number of stars, as indicated in the bar at the top. Stars in the HR 1614 moving group from Feltzing & Holmberg (2000) are shown as blue crosses and their best-selected candidates as black open circles. We identify stars within the red and black ellipses as members of the Hercules stream and the HR 1614 moving group, respectively (see discussion in Sect. 3). Dashed lines and numbers show the same boxes as in Feltzing & Holmberg (2000). |
|
In the text |
Fig. 2. Hertzsprung–Russel diagrams for stars inside the seven boxes as shown in Fig. 1. The dashed line is an isochrone for [M/H] = 0 and 8 Gyr, the solid line shows an isochrone for [M/H] = + 0.2 and 8 Gyr. The Colour bar shows the number of stars in each bin. |
|
In the text |
Fig. 3. Top left: location of five regions (00, 01a, 11a, 21a, and 31a) in the X − Y space explored in this work. Other plots: wavelet transform maps that show location of the kinematic overdensities in the U − V space. Their centres are shown with blue crosses. Lines connect eponymous overdensities into structures. The names of the structures are listed in the legend. The colour bar shows the normalised wavelet coefficients in each bin. |
|
In the text |
Fig. 4. Hertzsprung–Russel diagrams for 12 654 stars from the HR 1614 moving group (left) and 36 656 stars from Hercules (right) selected from the black and red ellipses shown in Fig. 1. An isochrone with [M/H] = + 0.2 and 2 Gyr is shown as a black solid line, and one with [M/H] = 0 and 8 Gyr is shown as a dashed line. Red dots in the top left plot show all the stars from Table 1 in Feltzing & Holmberg (2000). The colour bar shows the number of stars in each bin for both plots. |
|
In the text |
Fig. 5. Hertzsprung–Russel diagrams for the stars in the HR 1614 moving group selected from the cross-matched samples between StarHorse and APOGEE (red), and StarHorse and GALAH (blue), and StarHorse and SkyMapper for different metallicity bins. Black lines show isochrones for a mean metallicity of each bin and cover ages 0.2 Gyr and from 2 to 12 Gyr with the steps of 2 Gyr. |
|
In the text |
Fig. 6. Same as Fig. 5, but for the Hercules stream. |
|
In the text |
Fig. 7. Difference between GALAH and SkyMapper metallicities as a function of SkyMapper metallicity (left). Difference between APOGEE and SkyMapper metallicities as a function of SkyMapper metallicity (right). Both plots are colour-coded by [Mg/Fe] values taken from the GALAH and APOGEE samples, respectively, and are shown for stars in the HR 1614 moving group and the Hercules stream. Blue crosses show typical (mean) error. |
|
In the text |
Fig. 8. [Mg/Fe]–[Fe/H] diagrams for HR 1614 (left) and Hercules (right) selected from the APOGEE (top) and GALAH (bottom) surveys are shown as red dots. The corresponding uncertainties are shown as error bars for APOGEE distributions and a typical (mean) error is shown for GALAH stars. Background black distributions show stars from APOGEE and GALAH in region 00. |
|
In the text |
Fig. 9. Metallicity distributions for HR 1614 (left) and Hercules (right) from the GALAH (blue), APOGEE (red), and SkyMapper (green) surveys. Both regular histograms (upper three panels) and normalised cumulative histograms (bottom panel) are shown. |
|
In the text |
Fig. 10. Age–metallicity probability distributions for the HR 1614 moving group (left) and Hercules stream (right) for stars from the SkyMapper survey. The upper panels show the normalised age distributions after summing over the metallicity dimension. |
|
In the text |
Fig. 11. Distance relative to the Sun vs. metallicity diagram colour-cored by age for SkyMapper stars in the HR 1614 moving group (left) and the Hercules stream (right). |
|
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.