Open Access
Issue
A&A
Volume 668, December 2022
Article Number A4
Number of page(s) 16
Section Galactic structure, stellar clusters and populations
DOI https://doi.org/10.1051/0004-6361/202243590
Published online 28 November 2022

© X. Fu 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 Gaia mission (see Gaia Collaboration 2018a, 2021, 2022) is revolutionising our knowledge of the Milky Way (MW) with its very precise and accurate astrometry and photometry of more than 1.8 billion stars. Many of these stars are found in stellar clusters, which are important components of the Galaxy. In particular, open clusters (OCs) could trace the formation history and chemical properties of the Galactic disc. They also provide very useful tests of stellar evolution models (see e.g. Semenova et al. 2020; Magrini et al. 2021, for two examples based on results obtained by the Gaia-ESO Survey). Characterising OCs is therefore a fundamental task. Such a process includes discovering them, separating cluster populations from underlying field interlopers, measuring radial velocities (RVs) and chemical abundances, and deriving distances and ages. All these tasks can be performed more effectively by combining the Gaia results with ground-based data (see e.g. Bragaglia 2018; Carrera et al. 2019; Zhong et al. 2020; Casali et al. 2020b; Alonso-Santiago et al. 2021; Spina et al. 2021).

The most commonly used catalogues of OC properties before Gaia were Dias et al. (2002, and its web updates) and Kharchenko et al. (2013), in which 2000 to 3000 objects were considered, respectively. In the Gaia era, mostly thanks to the precise astrometric information (parallax, ϖ′ and proper motion, PM) and Gaia’s full-sky coverage, many new results of OCs have been reported. For instance, membership of OC stars has been studied using Gaia data release 1 (DR1) and the Tycho-Gaia Astrometric Solution (TGAS) by Gaia Collaboration (2017), Cantat-Gaudin et al. (2018a), Yen et al. (2018), Randich et al. (2018). Using Gaia DR2, member stars in known OCs have been identified (e.g. Cantat-Gaudin et al. 2018b; Cantat-Gaudin & Anders 2020; Jackson et al. 2022) and new OCs have been discovered (e.g. Cantat-Gaudin et al. 2019; Castro-Ginard et al. 2018, 2019, 2020; Beccari et al. 2018; Ferreira et al. 2019; Liu & Pang 2019). The Early Data Release 3 of Gaia has already been used to detect new OCs (e.g. Castro-Ginard et al. 2022). In particular, Cantat-Gaudin et al. (2020) combine highly reliable cluster membership in OCs – known and identified by other works – and estimate the age, distance, and reddening for ∼2000 OCs; their data set will be used in the present paper. All these works, together with the revision of the OC census, also mean that many candidate clusters have not been confirmed, see for instance the discussions in Kos et al. (2018) and Cantat-Gaudin et al. (2018b).

The kinematical information of OCs is based on RV measurements of stellar spectra. While Gaia spectroscopic capabilities are limited (see e.g. Sartoretti et al. 2018; Katz et al. 2019), its instrument Radial Velocity Spectrometer (RVS, Cropper et al. 2018) collected data for several million bright stars. Matching the Gaia RVs to OC members (derived by Cantat-Gaudin et al. 2018b), Soubiran et al. (2018) could determine average RVs for nearly 900 clusters and derive their kinematics. The work has been extended by Tarricq et al. (2021), who also included data from ground-based surveys.

Gaia data have already been extensively adopted to clean the colour-magnitude diagram (CMD) of OCs, and to derive a more precise age (and distance, if not computed directly from the cluster ϖ), see for instance Randich et al. (2018), Yalyalieva et al. (2018), Dias et al. (2018), Choi et al. (2018) and the method described in Li & Shao (2022). An extensive derivation of cluster ages can be found in Gaia Collaboration (2018b), Bossini et al. (2019). In the latter paper, a Bayesian code was applied to the list of clusters in Cantat-Gaudin et al. (2018b) and they were able to obtain excellent age results for about 270 OCs. However, a common limitation of these works is the absence of metallicity information in the majority of Galactic clusters, which introduced degeneracies with reddening and age.

In fact, metallicity measured with high-resolution spectroscopy is available only for about 10% of the whole OC population. This low percentage of highly resolved OC metallicity is not due to a shortage of studies. To cite only a few, Magrini et al. (2018) combined Gaia-ESO Survey data with compilations from Netopil et al. (2016), while Donati et al. (2015), Reddy et al. (2016), Reddy & Lambert (2019), Casamiquela et al. (2017), Smiljanic et al. (2018), Bragaglia et al. (2018), Casali et al. (2020a) are based on private projects such as BOCCE, OCCASO, and SPA. Although Gaia will obtain the metallicity and some elemental abundances on a grand scale with RVS, to derive more detailed elemental abundances, ground-based surveys are fundamental, and especially in faint clusters. Examples are the high-resolution Gaia-ESO (Gilmore et al. 2012; Randich et al. 2022), APOGEE (Majewski et al. 2017), GALAH (De Silva et al. 2015), and the low-resolution LAMOST (Cui et al. 2012; Deng et al. 2012; Zhao et al. 2012) surveys. Future surveys such as WEAVE (Dalton et al. 2012) and 4MOST (de Jong et al. 2019) are also planning to observe the OC population.

The Gaia-ESO Survey targeted 62 OCs on purpose, observing from a few hundred to thousands of stars at an intermediate resolution, and roughly tens of members at high resolution in each of them (Bragaglia et al. 2022; Randich et al. 2022). About 20 more clusters from the ESO archive observations were also re-analysed homogeneously and included in the data release (see for instance Magrini et al. 2017; Bragaglia et al. 2022; Randich et al. 2022). The earlier data release of the main GALAH survey does not have OCs (see e.g. Buder et al. 2018), while their latest release covers 75 OCs (DR3, Buder et al. 2021), which were also analysed together with the APOGEE OCs to provide a homogeneous set (Spina et al. 2021). The APOGEE OC samples are mainly presented within the Open Cluster Abundances and Mapping (OCCAM) program (Donor et al. 2018, 2020), with a few to a few tens of stars in each OC. However, more OC stars have been serendipitously observed both by GALAH and by APOGEE, as found by Carrera et al. (2019). They cross-matched the OC member stars as defined by Cantat-Gaudin et al. (2018b) with the survey data releases and were able to retrieve RVs, metallicities, and chemical abundances for more than 100 OCs, many of them without previous determinations.

The same technique can be applied to the LAMOST survey, to extend the number of OCs with measured RV and metallicity, and investigate their chemical, kinematical, and dynamical properties on Galactic scales, with a catalogue of homogeneous analysis. Based on a previous data release (LAMOST DR5) and an earlier Gaia OC membership catalogue (Cantat-Gaudin et al. 2018b), Zhong et al. (2020) explored properties of 295 clusters and discussed their metallicity distributions in the MW. The latest LAMOST data release, DR81, includes 10 388 423 stellar spectra in total, which were observed between October 24th, 2011 and May 27th, 2020. This work is an updated and extended version of Zhong et al. (2020) on the LAMOST OC investigations based on Gaia.

This paper is organised as follows. In Sect. 2 we introduce the LAMOST and Gaia data adopted in this work, together with quality control methods. The catalogue results after the quality control (i.e. the RV and metallicity of the LAMOST OCs), are described in Sect. 3. In Sect. 4 we discuss the Galactic metallicity distribution obtained with our LAMOST OC catalogue, the dynamical properties of OCs, and the connection with the Galactic molecular clouds. Lastly, the main conclusions of this paper are summarised in Sect. 5.

2. Data and quality control

The Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST, Guo Shou Jing Telescope), located in Xinglong, China, is a quasi-meridian reflecting Schmidt telescope with an effective aperture of ∼4 m. With its 4000 fibres, it is one of the most efficient spectroscopic telescopes. In the low-resolution mode (R = 1800), its limiting magnitude is r = 19 mag.

In this work we adopt the LAMOST DR8 low-resolution catalogue with stellar parameters (namely, the LAMOST LRS A, F, G and K Star Catalog). In total, 6 478 063 spectra are published in the original LAMOST catalogue, including 100 468 A-type, 1 983 821 F-type, 3 249 746 G-type, and 1 144 028 K-type stars. All spectra in this catalogue have a criterion of g-band signal-to-noise ratio, S/N > 6 for dark night observations, or S/N > 15 for bright night observations. RV and stellar parameters (i.e. effective temperature Teff, surface gravity log(g), and iron abundance [Fe/H]) in this catalogue are determined with the official LAMOST Stellar Parameter pipeline (LASP, Wu et al. 2014), which uses ATLAS9 atmosphere models (Castelli & Kurucz 2003) and the Grevesse & Sauval (1998) Solar abundances.

To obtain RV and stellar parameters of OC member stars, we cross-matched the LAMOST and Cantat-Gaudin et al. (2020) catalogues, keeping only members with probabilities > 70%. For each star, we matched the above two catalogues with its Gaia source_id. The Cantat-Gaudin et al. (2020) catalogue already had the Gaia DR2 source_id, so the first step of our procedure was to match each LAMOST star with the Gaia DR2 data. We used the CDS X-match service in TOPCAT (Taylor 2005, 2017) to consider both the PM and the epoch of Gaia stars. Targets were identified with their RA and Dec. coordinates within 3.5 arcsec, because the fibre size of LAMOST is 3.3 arcsec (Zhao et al. 2012) and the dome seeing of LAMOST is sometimes slightly larger than the fibre scale (Luo et al. 2015). In the X-match procedure, all Gaia sources in the matching radius were considered because light from different sources cannot be resolved within the LAMOST fibre size. In some cases more than one Gaia sources were matched for a LAMOST observation. However, our determinations of the clusters’ velocity and [Fe/H] were not affected because a Monte Carlo sampling selection were applied in the procedure (see details in Sect. 2.1). After getting the Gaia DR2 source_id of each star in the LAMOST catalogue, we used the source_id as the identification to cross-match the LAMOST-Gaia table with the Gaia OC member star catalogue (Cantat-Gaudin et al. 2020). In total, 7570 stars from 386 clusters are in common. The spectra of these stars have a high S/N; the mean S/N in g-band is 109 and ∼78% stars have a g-band S/N > 30. Figure 1 shows the sky coverage of the LAMOST OC members (red dots). For comparison, all OC member stars of Cantat-Gaudin et al. (2020) are overlaid as grey dots. Since LAMOST operates from the northern hemisphere, the southernmost matched stars have a declination of ∼ − 7°. Among all the LAMOST OCs, 308 OCs have at least two matched stars, while the other 78 clusters have only one matched star in the table. We keep all these clusters for further analysis. In principle, even one star can represent the cluster because OCs are very homogeneous in their chemistry and RV.

thumbnail Fig. 1.

Sky coverage of the LAMOST open cluster (OC) members (red dots) and all the Cantat-Gaudin et al. (2020) OC member stars (grey dots).

To provide accurate chemical and kinematical information for these selected LAMOST OCs, we performed quality control on RV and stellar parameters, respectively.

2.1. Radial velocity quality control

The greatest advantage of LAMOST is the vast number of spectra, but the low resolution makes the stellar RV uncertainty not negligible compared to studies with high-resolution spectra (HRS). The median value of the RV uncertainty for our matched LAMOST OC member stars is 6.36 ± 3.20 km s−1, which is comparable to, or even larger than, the typical RV dispersion of a cluster. For instance, the RV dispersion of clusters NGC 2516, NGC 6705, and NGC 6633 are 1.0, 1.6, and 1.5 km s−1, respectively, as reported in the Gaia-ESO survey based on HRS (see e.g. Magrini et al. 2017).

To derive the average RV of clusters, the RV dispersion of each cluster, and a proper RV member quality control, we applied a Monte Carlo (MC) method, by considering each individual RV and its corresponding uncertainty from the LAMOST DR8 catalogue. For simplicity we assumed that all RV uncertainties follow Gaussian distributions and then we randomly sampled the RV of every member star in each cluster 5000 times. This allowed us to derive the mean and median of all the sampled RV values. Member stars with RV measurements within 2σ of the MC RV sampling are marked with a quality control flag of FLAG = 1, while those with > 2σ are marked as FLAG = 0. The mark of FLAG = 1 means we used the RV values of the star to determine the cluster RV.

Figure 2 shows an example of the RV quality control in cluster NGC 2548, where RV of all matched stars are plotted as Gaussian profiles. Dark-grey curves show selected RV members (FLAG = 1), while light-grey curves are stars we discarded (FLAG = 0) for the cluster RV determination. The thick blue curve is the Gaussian fit of all the stars with FLAG = 1 in this cluster. The thick blue vertical line represents the mean RV of the cluster (Vrad, mean). The two vertical dashed lines mark the 2-σ departure from Vrad, mean, which are adopted as thresholds of FLAG = 1 star selections. In the top left region of Fig. 2, we show the numbers of total matched stars in the cluster (Nstar), numbers of RV-selected (FLAG = 1) stars (Nstar, RV), the cluster mean RV (Vrad, mean), and the corresponding standard deviation (σVrad), which is adopted as the cluster RV uncertainty. In most clusters, Vrad, mean and the median RV of the cluster Vrad, med are very similar to each other, with a mean absolute difference of ∼0.85 km s−1. Indeed, 76% of clusters have a difference smaller than this mean value. In the rest part of this paper, we use Vrad, med to discuss the property of clusters’ RV Vrad.

thumbnail Fig. 2.

NGC 2548 as an example of the RV selection. The RVs of all the crossed-matched member stars in the cluster are plotted assuming a normal distribution. The light grey colour curves denote stars we discarded for further analysis (FLAG = 0), while the dark grey ones are the selected stars. The two vertical dashed lines mark the 2-σ selection criterion. The Monte Carlo model result of all the selected stars is illustrated with the blue Gaussian curve. Their mean RV value (Vrad, mean), the 1-σ RV dispersion (σVrad), as well as the number of member stars before (Nstar) and after the RV selection (Nstar, RV), are shown in the legend.

Stellar RV measurements of LAMOST are known to have a systematic offset compared to higher-resolution data, with a value2 of about 4–5 km s−1. In the ‘Survey of Surveys’ work by Tsantaki et al. (2022), where they compare RV measurements from different survey data, the median and mean difference between Gaia DR2 and LAMOST DR5 (RVGaia DR2 − RLAMOST) are 4.97 and 5.18 km s−1, respectively. In Fig. 3 we show the RV measurement difference (ΔRV = RVGaia DR2 − RLAMOST) distribution of all the FLAG = 1 OC members. The ΔRV has a median value of 4.66 km s−1, a mean value of 4.82 km s−1, and a standard deviation of 12.29 km s−1. All these values, based on cluster member stars after the RV quality control, are similar to the raw values reported by Tsantaki et al. (2022) before their RV correction.

thumbnail Fig. 3.

Histogram of the RV difference between Gaia DR3 measurements and LAMOST measurements for all the RV selected members. The median value of the difference ΔRV = RVGaia DR3 − RVLAMOST is 4.66 km s−1.

2.2. Stellar parameter quality control

Similar to the quality control of Vrad described in Sect. 2.1, it is also necessary to control the qualities of stellar parameters before further analysis. This control procedure is based not only on finding [Fe/H] outliers, but also checking member star evolution. In principle, the surface gravity log(g) and effective temperature Teff of member stars in the same cluster should follow the evolution of a simple stellar population in the Kiel diagram, where the stellar log(g) value is an index of the evolutionary phase. The [Fe/H] of these stars, on the other hand, should be almost a constant in different evolutionary phases3.

Therefore, such quality control could serve as useful examinations, or even calibrations, to check the accuracy of stellar parameter analysis. In these quality control processes for our LAMOST OC member stars, we find a systematic issue for the log(g) and [Fe/H] parameters of cool main sequence stars. Figure 4 presents cluster NGC 2632 as an example of such an issue. We find that main sequence stars with Teff ⪅ 5000 K do not follow the main sequence evolution trend. For reference, a theoretical isochrone from PARSEC v1.2S4 (Bressan et al. 2012; Chen et al. 2014, 2015; Tang et al. 2014) with a metallicity of Z = 0.013 and log(age) = 9.3 is plotted in the Kiel diagram (the upper panel of Fig. 4). Most of these stars have a log(g) value lower than the expected values (see stars in the grey part of the upper panel of Fig. 4). Their [Fe/H] values are also abnormally lower than those of other member stars (see the lower panel of Fig. 4).

thumbnail Fig. 4.

Kiel diagram (upper panel, colour-coded with [Fe/H]) and [Fe/H] as a function of Teff (lower panel, colour-coded with log(g)) of cluster NGC 2632. For reference, a PARSEC theoretical isochrone of Z = 0.013, log(age) = 9.3 is plotted in the Kiel diagram. The line-shaded region marks Teff < 4000 K, and the grey part marks 5000 K < Teff < 4000 K.

On the other hand, giant stars in this temperature regime seem unaffected. Member stars’ [Fe/H] values are almost constant for all giants, as well as for Teff ≳ 5000 K dwarfs. Only the dwarf stars show a declining trend at lower temperatures, which is seen in all of our clusters that have dwarf stars in this temperature range. To select OC member stars with a more secure [Fe/H] determination, we discarded dwarf stars cooler than 5000 K in our cluster [Fe/H] calculation.

It is always difficult to obtain stellar parameters for very hot and very cool stars. There are few metal lines that can be used in hot star spectra, and the continuum is relatively difficult to derive. For stars with cool temperature, molecular lines dominate their spectra and make the analysis challenging (see discussions in Jofré et al. 2019). Therefore, the second step of our stellar parameter quality control process was excluding [Fe/H] of member stars with Teff > 7500 K and Teff < 4000 K for the cluster metallicity determinations. The hot star criterion Teff = 7500 K was adopted from the LAMOST hot star stellar parameter work (Xiang et al. 2022).

After the two-step quality control process described above, we ended up with 355 OCs for further cluster [Fe/H] determination. The selected member stars are marked with FLAG=12 in our output catalogue, which means the parameters of these stars are good for both RV and [Fe/H] determination of the cluster. Among these clusters, 203 have at least three members with FLAG=12.

We then calculated the cluster [Fe/H] together with the corresponding uncertainty and scatter using a Monte Carlo method. The mean stellar [Fe/H] uncertainty of our OC member stars from the pipeline LASP is 0.07 dex. Assuming [Fe/H] uncertainty of member stars follow Gaussian distributions, we applied a random [Fe/H] sampling of 5000 times to each cluster, and obtained the median [Fe/H]med, mean [Fe/H]mean, and standard deviation σ[Fe/H] of all the sample values. The absolute difference between the cluster median [Fe/H]med and mean [Fe/H]mean values are very small, with a mean difference of ∼0.01 dex. About 70% of clusters have a difference smaller than this value. In the rest of this paper, we adopt the median value [Fe/H]med to discuss the cluster [Fe/H], and take the sample standard deviation σ[Fe/H] as the cluster [Fe/H] uncertainty.

Figure 5 shows the CMDs of member stars with Gaia DR2 photometry of two very typical OCs in our catalogue. They are Melotte 22 (Pleiades) in the upper two panels and NGC 2682 (M67) in the lower two panels. The left panel of each cluster is colour-coded with member stars’ [Fe/H], which means all the coloured ones are FLAG=12 stars. The right panel of each cluster is colour-coded with member stars’ RV, marking all the FLAG=1 stars in the cluster. The derived cluster Vrad and [Fe/H] values based on the FLAG=1 and FLAG=12 members, respectively, are shown in the figure.

thumbnail Fig. 5.

Colour-magnitude diagrams of Melotte 22 and NGC 2682 as examples of the clusters with [Fe/H] in the literature. The filled grey circles are high-quality member stars from Cantat-Gaudin et al. (2020), and open circles are all matched stars with LAMOST RV measurements. Right panel: stars good for the cluster Vrad determination (FLAG=1) are colour-coded with their RV values. Left panel: the colour-filled dots are stars good for [Fe/H] determination (FLAG=12), and the black dots are stars without [Fe/H] determination or with discarded [Fe/H]. The dashed curve in each sub-figure is the PARSEC isochrone with the Cantat-Gaudin et al. (2020) cluster parameters.

3. Results

3.1. Radial velocity

3.1.1. Clusters with Vrad in the literature

Catalogues of the Galactic OC RV have been compiled in many works (see e.g. Dias et al. 2002, 2021; Kharchenko et al. 2013; Conrad et al. 2014; Soubiran et al. 2018; Tarricq et al. 2021). We compared our cluster Vrad results to the two most recent and most complete large catalogues, from Tarricq et al. (2021) and Dias et al. (2021), together with the previous LAMOST OC catalogue from Zhong et al. (2020).

Among all of our 386 OCs, we found 308 in common with Tarricq et al. (2021), 185 in common with Dias et al. (2021), and 226 in common with Zhong et al. (2020). Figure 6 shows the comparison of the cluster Vrad between our results and these three catalogues. The distributions of the Vrad difference are also displayed as a histogram in each panel. Compared to the previous LAMOST OC results from Zhong et al. (2020), the median value of the Vrad difference is only 0.06 km s−1(see the right panel of Fig. 6), which indicates that the two catalogues are consistent in Vrad. The median Vrad difference between our results and those from Tarricq et al. (2021) and Dias et al. (2021) have a value around 6.2 km s−1 (see the left and middle panels of Fig. 6). These differences should, at least partially, be due to the aforementioned systematic offset of LAMOST RV.

thumbnail Fig. 6.

Comparison of our cluster Vrad to results in the literature. Left, middle, and right panels: comparison to 308, 185, and 226 OCs of Tarricq et al. (2021), Dias et al. (2021), and Zhong et al. (2020), respectively. The histograms of the cluster RV differences are shown as insets in each panel.

In total, 342 of the OCs in our catalogue already have a Vrad measurement in one of these three recent large catalogues. Their Vrad, σVrad, and number of stars used in the calculation, together with the literature values, are listed in Table 1.

Table 1.

Radial velocity of 342 OCs with known Vrad in the literature.

3.1.2. Clusters with newly obtained Vrad

The large amount of LAMOST spectra enabled us to obtain cluster Vrad values that were not reported in the literature before. We report here 44 clusters with newly obtained Vrad. Table 2 lists their information. These values, together with those listed in Table 1, will be employed later in the paper, for a kinematical analysis of the OC population.

Table 2.

Newly obtained cluster radial velocity of 44 OCs.

3.2. Metallicity

As groups of stars form together with the same composition, OCs are an ideal probe to trace the metallicity evolution of the Galactic disc. With OCs, in fact, not only we can improve the S/N by averaging over several member stars, but also we could deal with objects whose age can be determined with a lower uncertainty than that of field stars (with the only possible exception of stars for which a full asteroseismologic analysis is possible, see e.g. Rodrigues et al. 2017; Miglio et al. 2021).

In this section we report the metallicity of LAMOST OCs, which consists of 218 clusters with known metallicity collected from the literature, and 137 clusters with newly obtained spectroscopic metallicity. Here we use the iron abundance [Fe/H] as an index of the stellar metallicity.

3.2.1. Clusters with [Fe/H] in the literature

We collect cluster [Fe/H] measurements from HRS including results of large surveys such as Gaia-ESO (Spina et al. 2017; Bragaglia et al. 2022; Randich et al. 2022, and references therein), APOGEE (Donor et al. 2020; Carrera et al. 2019; Spina et al. 2021), GALAH (Carrera et al. 2019; Spina et al. 2021), large projects such as SPA (Origlia et al. 2019; Casali et al. 2020a; Zhang et al. 2021; Alonso-Santiago et al. 2021), OCCASO (Casamiquela et al. 2016, 2017, 2018), and high-quality (HQ) collections from Netopil et al. (2022). In general, the HRS results we use are similar to the HRS used in Zhang et al. (2021), but with more OCs from Spina et al. (2021) and the HQ results of Netopil et al. (2022). In total, we have 82 clusters in common with HRS [Fe/H] results. The left panel of Fig. 7 shows the comparison between our results based on LAMOST and the HRS values. The difference between the cluster [Fe/H] values of these two data sets is 0.01 ± 0.15 dex.

thumbnail Fig. 7.

Comparison of 82 clusters’ [Fe/H] to literature results with high-resolution spectra (HRS) in the left panel, and of 206 clusters to results from Zhong et al. (2020). The X symbols mark [Fe/H] values with less than three member stars, namely 38 clusters on the left (from HRS) and 75 clusters on the right (from Zhong et al. 2020). The open squares mark clusters with less than three FLAG=12 member stars from this work, which are 23 clusters on the left panel and 92 clusters on the right panel.

Clusters with less than three member stars for the [Fe/H] determination, both in our case (marked with open square of Fig. 7) and in the HRS works (marked with X), show a relatively large difference between the two sets of results. If we compare clusters with at least three [Fe/H] member stars in both data sets, their difference is 0.00 ± 0.11 dex, which shows a very good consistency. Therefore, in the following discussions on the Galactic metallicity evolution, we only consider clusters with at least three members and with FLAG=12.

We also compared our results to the previous LAMOST OC work by Zhong et al. (2020), which is based on an earlier data release and a different member star selection. There are 206 clusters in the Zhong et al. (2020) catalogue with [Fe/H] determination, and the right panel of Fig. 7 shows the [Fe/H] comparison. For all the matched clusters, the [Fe/H] value mean difference is 0.02 ± 0.23 dex. When only clusters with at least three [Fe/H] member stars are considered, their difference is 0.04 ± 0.06 dex.

Combining both literature results from HRS studies and the low-resolution spectra study of Zhong et al. (2020), 218 of our clusters have a reported metallicity. Table 3 lists their information.

Table 3.

218 clusters with [Fe/H] values in the literature.

3.2.2. Clusters with newly obtained [Fe/H]

Here we report [Fe/H] measurements of 137 OCs that do not have spectroscopic [Fe/H] in the literature. Among them, 63 clusters have at least three FLAG=12 member stars. Table 4 and Table 5 list our results for these clusters. Clusters with less than three FLAG=12 member stars are available in our final catalogue as well, but their metallicities should be used with caution. This constitutes an important addition to the number of OCs with metallicity determined on the basis of spectroscopy.

Table 4.

Newly obtained iron abundances for 63 OCs with at least three FLAG=12 member stars.

Table 5.

Newly obtained iron abundances for 74 OCs but with fewer than three FLAG=12 member stars.

4. Discussions

4.1. The Galactic metallicity distribution and dynamical properties of the LAMOST OCs

In order to study the metallicity distribution of the MW disc, one of the most common subjects investigated is the so-called radial metallicity gradient. The metallicity gradient is often displayed with [Fe/H] as a function of the Galactocentric radius RGC, indicating different levels of metal enrichment along the Galactic radius. Stellar OCs, being groups of stars in a simple stellar population, are ideal probes to trace the metal enrichment history of the MW.

With the homogeneous Vrad and [Fe/H] determination in our newly compiled LAMOST OC catalogue, we can investigate the evolution of the Galactic metallicity – even its gradient evolution – in the past 500 Myr. This 500-Myr range traces back to the time when the Sagittarius dwarf galaxy (Sgr dSph) had its last passage through the MW outer disc (see for instance the discussions in Xu et al. 2020). Sgr dSph is also the last relatively sizeable minor merger in the evolution history of our Galaxy known to date, which means the gravitational potential of the Galaxy can be considered as a constant in the past 500 Myr. A stable gravitational potential is further the guarantee of orbit calculations for these clusters.

To trace back the orbits of LAMOST OCs and their evolution, we used the publicly licensed code GALPOT5 (Dehnen & Binney 1998) and the MW gravitational potential from McMillan (2017). The gravitational potential takes into account the Galactic thick and thin stellar discs, a bulge component, a dark-matter halo, and a cold gas disc. We adopted the Solar motion of (U, V, W) = (11.1, 12.24, 7.25) km s−1 (Schönrich et al. 2010), which was insensitive to the metallicity gradient of the MW disc. The Galactic radius of the Sun is 8.2 kpc, the solar circular speed is 232.8 km s−1, both suggested by McMillan (2017). Proper motions and distances of clusters were adopted from Cantat-Gaudin et al. (2020).

Figure 8 shows the Galactic metallicity gradient traced by LAMOST OCs in six snapshots during the past 500 Myr, that is, from the present time (looking back time 0 Myr is in the uppermost left panel) back to the time of the last passage of Sgr dSph. It is very difficult for LAMOST to observe stars towards the Galactic centre, so almost all the LAMOST OCs are located with Galactocentric radii RGC ⪆ 8 kpc. As discussed earlier, only clusters with no less than three FLAG=12 stars are used in the following metallicity discussions. Here we also make the simplification that the cluster [Fe/H] does not change in its life time (i.e. the microscopic diffusion in the stellar surface layer is not considered). To investigate the metal enrichment in different epochs, we divided the LAMOST OCs into six age groups (colour-coded). A second-order polynomial fitting was applied to each age group to measure their metallicity trend:

[ F e / H ] = a × R GC 2 + b × R GC + c . $$ \begin{aligned} \mathtt [Fe/H] = {a \times R_{\rm GC}^2 + b \times R_{\rm GC} + c}. \end{aligned} $$(1)

thumbnail Fig. 8.

Clusters’ median [Fe/H] as a function of the Galactocentric radius RGC for all the LAMOST clusters with at least three FLAG=12 stars. Different panels show different look-back times according to the orbit calculations. Clusters are divided into six groups according to age from the Cantat-Gaudin et al. (2020) catalogue. The solid curves are the second-order polynomial MCMC fittings of each age group.

For each fitting, we used a Markov chain Monte Carlo (MCMC) method to consider the [Fe/H] uncertainty of each cluster. A number of 300 walkers and a 10 000-step MCMC were used in each fitting, based on the MCMC python package emcee (Foreman-Mackey et al. 2013). The parameters (a, b, c) to the second-order polynomial fits of the current time are listed in the RGC column of Table 6.

Table 6.

Second-order polynomial MCMC fitting parameters of the clusters’ [Fe/H] as a function of the Galactic radius (RGC), the Z component of the Galactic angular momentum (LZ), and the total Galactic energy (En) for six age groups.

In the literature, a linear fit has been widely adopted to describe the metallicity radial gradient (see e.g. Friel et al. 2002; Bragaglia & Tosi 2006; Carrera et al. 2019; Zhong et al. 2020; Zhang et al. 2021), but a break or a ‘knee’ point is needed at ∼12–14 kpc in order to properly take into account the outer disc clusters (see e.g. Reddy et al. 2016; Donor et al. 2020; Spina et al. 2022). Indeed, radial metallicity distribution studies in other galaxies suggest that breaks and changes of slopes are required in the linear fittings, and non-linear models often fit better than linear models (Scarano & Lépine 2013). For this reason we use the second-order polynomial model fittings to illustrate the metallicity trends. From the six snapshots of the [Fe/H] trend in Fig. 8, it is clear that, although generally speaking [Fe/H] is higher at smaller RGC, the metallicity gradients evolve among age groups, instead of being a constant curve. In other words, if we were observing OCs at different times in the past 500 Myr, we would have seen very different pictures of the metallicity gradient(s), covering multiple shapes of being steep, flat, with a turning point, and close to a straight line.

In order to interpret radial metallicity gradient variations for different age groups, we modelled the circular orbits of each cluster with GALPOT. We note that radial metallicity gradients are themselves trajectory projections of cluster orbits. We find that the clusters are not travelling in circular orbits. Figure 9 shows the eccentricity (ecc) of LAMOST OCs as a function of their Galactic apocentre distance. Although GALPOT does not give uncertainties of the orbit parameters, we calculated the possible orbit parameter changes according to the scatter of clusters’ Vrad.

thumbnail Fig. 9.

Clusters’ eccentricity as a function of their Galactic apocentric distance. The colour code is the log value of the cluster age in the left panel, and clusters’ [Fe/H] in the right panel. Only clusters with multiple members are displayed in the left panel, while those with at least three FLAG=12 stars are displayed in the right panel.

Most clusters have a non-circular orbit (ecc > 0). For instance, Berkeley 29 has at present the highest RGC = 20.44 kpc, being located on the outermost disc in our catalogue. However, its apocentric and pericentric distances are 20.54 kpc and 13.92 kpc, respectively, with an eccentricity of 0.19. The most eccentric cluster with at least three FLAG=12 member stars in our sample is Berkeley 32 (ecc = 0.31). It has a current RGC of 11.00 kpc, but it travels to 12.23 kpc and 6.42 kpc as its apocentre and pericentre, respectively.

Generally speaking, the youngest clusters are more likely to travel with a more circular orbit, while the oldest clusters have higher chances of having larger eccentricity (see the left panel of Fig. 9). The metal-poorer OCs, whether or not they have a circular orbit, can reach larger apocentric distances than the metal-richer clusters (see the right panel of Fig. 9). The non-circular orbit around the Galactic centre results in clusters that constantly and periodically alter their RGC between their pericentre and apocentre. The alteration period is shorter than one revolution time of the cluster, making the cluster travel back and forth several times in one revolution around the Galactic centre (see also discussions in Lépine et al. 2011). All these facts make the present-day Galactocentric radius RGC not an invariant in the Galactic metallicity distribution investigations.

Thanks to the astrometry results provided by Gaia and RV provided by spectroscopic facilities such as LAMOST, we can now use more invariant parameters such as the angular momenta and the dynamical energy to unfold the whole story of the Galactic [Fe/H] evolution in the Gaia era. In Fig. 10 we display the [Fe/H] distribution of the LAMOST OCs as a function of the Z component of their angular momentum (left panel) and the total dynamical energy (right panel). The corresponding uncertainties reflect the scatter of clusters’ Vrad. These two parameters take into account both the spatial positions and the motions, so they are constant for a cluster in a given gravitational potential. We calculated them with GALPOT following the same setups described before in the orbit integration part. The Z component of the angular momentum (LZ), namely the angular momentum component that describes the motion around the Galactic centre on the Galactic plane, is positive in the direction of the Galactic disc movement and has a larger value at the outer disc. The dynamical energy (Energy) describes the total energy of a cluster including kinetic and potential energy. For disc stars with normal motions, a location further from the Galactic centre usually leads to a higher value of the total energy.

thumbnail Fig. 10.

Clusters’ median [Fe/H] as a function of their Z component of the Galactic angular momentum (left panel) and the total Galactic energy (right panel) for all of the LAMOST clusters with no less than three FLAG=12 stars. The solid curves are the second-order polynomial MCMC fittings of each age group.

As shown in Fig. 8, OCs are divided into six groups based on their current age and only clusters with at least three FLAG=12 stars are included in the figure. The second-order polynomial MCMC fitting to each age group is calculated with the same settings as described earlier. The fitting parameters and their 16% and 84% distribution in the MCMC samples are listed in Table 6. The metallicity distribution along LZ or energy could be considered as the Galactic metallicity gradient form in the Gaia era.

From the [Fe/H] trend, we could find that the very old age groups, those with age > 4 Gyr (grey curves) and within an age range of [2 Gyr, 4 Gyr] (plum curves), both show a flat tail at their low-metallicity end. Such a flat metallicity tail has [Fe/H] ∼ −0.2 dex for the age > 4 Gyr group, starting from LZ ∼ 2.1 × 103 kpc km s−1 and energy ∼ − 1.48 × 105 km2 s−2, and [Fe/H] ∼ −0.4 dex for the [2 Gyr, 4 Gyr] group, starting from LZ ∼ 3.×103 kpc km s−1 and energy ∼ − 1.3 × 105 km2 s−2, respectively. The other four younger OC groups do not show such a flat tail feature.

The flat metallicity trend in the outer part of a Galactic disc is known in the literature, both in our own Galaxy (see e.g. Lépine et al. 2011; Andreuzzi et al. 2011; Spina et al. 2022) and in other galaxies (Scarano & Lépine 2013). However, here we show that the OC metallicity gradients appear to depend on their age, and the two flat tails display two different metallicity “plateaus”.

The nature of the flat metallicity tail is still an open question. Lépine et al. (2011) propose a hypothesis that clusters at the outer Galactic disc with a flat [Fe/H] trend may have their birth place at a smaller RGC. These clusters travel between their orbit pericentre and apocentre, and are observed in their present day position at a larger RGC. This hypothesis could explain the absence of young clusters at the outer Galactic disc and allow a simple trend of metallicity along the Galactic radius. Our metallicity trend with LZ and energy, however, considers both spatial position and motions, and thus should not be affected by orbit.

The existence of the flat metallicity tail on the LZ and energy plane, especially the two different flat metallicity tails, must come from an alternative mechanism. Another possible hypothesis proposed by Lépine et al. (2011) is the gas flow from the relative inner region (the corotation radius) of the Galaxy to the external regions, which brings gas with a slightly higher metallicity to the outer regions and flattens the metallicity gradients. The driver of the gas flow could be the gas interaction with the spiral potential perturbation (see detail discussions and simulations in Lépine et al. 2001). With the current set of OC data, we are not able to confirm or exclude this hypothesis. Detailed chemical abundances of the outer disc clusters may help to revive the initial gas composition of these clusters. Whether the gas flow hypothesis can also explain the other flat metallicity tail, that of the age > 4 Gyr clusters, requires further investigations.

Another key to understanding the metallicity evolution history of the flat metallicity tail is to investigate the origin of clusters in the Galactic outer disc, namely the clusters with high values of LZ and energy. Such clusters are visible on the right part of both panels in Fig. 10. We ask ourselves whether these clusters are native residents of our Galaxy, or if they have an extragalactic origin. For instance, the formation and origin of Berkeley 29 – which is the outermost known disc cluster with a RGC of ∼20 kpc – are discussed in the literature (Carraro & Bensby 2009). In their work, Carraro & Bensby (2009) find that the trailing tail of Sgr dSph passes close to the location of Berkeley 29, and their [Mg/Fe], [Ca/Fe] abundances are similar to each other at the same [Fe/H]. They suggest that Berkeley 29 was formed in Sgr dSph and was left on the Galactic disc during one of Sgr dSph’s passages through our Galaxy about 5 Gyr ago. In fact, the mean metallicity of dwarf satellite galaxies may be higher than that of the original outskirts stars of the MW disc. After the satellite galaxy was disrupted by the MW tidal field, the stripped stars distributed mainly on the outer region of MW due to their high kinetic energy and angular momentum. These processes will increase the overall metallicity of the outskirts and form a flat metallicity tail. In a follow-up paper (Chang, in prep.), we will discuss in detail the role of minor mergers in forming the flat metallicity tail.

One of the most popular and efficient methods to identify different Galactic structures is the integral of motion (IOM) space (see discussions in Helmi 2020, and references therein). The main idea of this method is using the total energy and the Z component of angular momentum as a long-term record. Stars or clusters with a similar origin tend to cluster in the energy-LZ space.

Figure 11 shows the IOM of LAMOST OCs. Most clusters lie on a compact region on the IOM space, which represents the dynamical property of the Galactic thin disc. However, a few outliers do exist, in which we highlight a cluster with a small value of both LZ and Energy in the inner part of the Galaxy (NGC 6791), clusters with large values of LZ and Energy in the outer part of the Galaxy (Berkeley 29, Berkeley 34, and Berkeley 19), and clusters showing a departure from the main thin disc property (Berkeley 32, UPK 39, Berkeley 17, and NGC 7128). These clusters must either have a different origin, or have experienced different evolution compared to the majority of normal OCs of the Galactic thin disc.

thumbnail Fig. 11.

Integral of motion (IOM) of LAMOST OCs with multiple members (Nflag = 1≥ 2). The [Fe/H] colour-code is applied to clusters with at least three FLAG=12 members. The empty symbols are clusters with less than three FLAG=12 members.

Similarly to the Galactic field stars, clusters inhabit and travel on both sides of the Galactic plane. The scale height of the Galactic disc is shaped by the disc stars and OCs. It increases quickly from the Solar RGC to the outer disc on both sides of the plane. This increase in Z distance departure from the Galactic plane is referred as the disc flare. The flare has already been traced with different types of stars, such as blue stragglers (see e.g. Thomas et al. 2019), red clump stars (see e.g. Wan et al. 2017), and red giants (see e.g. Wang et al. 2018).

Using simulations and LAMOST K giant stars, Xu et al. (2020) suggest that the last impact of Sgr dSph contributes to the flare, and they identify three different branches of the flare with rotation velocity Vϕ. According to their results, the boundary of the flare is constructed by stars with a similar Vϕ to the main part of the disc stars. The Vϕ value for the flare boundary stars is very different compared to halo stars at the same RGC and stars on the disc plane at larger RGC. In Fig. 12 we follow the analysis of Xu et al. (2020) to investigate the LAMOST OCs distribution on the Z − RGC plane. All clusters are colour-coded with their Galactic Vϕ.

thumbnail Fig. 12.

Present-day location of LAMOST OCs on the ZGC − RGC plane. Upper panel: all LAMOST OCs, while the middle and lower panels show clusters born after the last passage of Sgr dSph (age < 500 Myr), and clusters born before the passage (age > 500 Myr), respectively. Clusters are colour-coded with their azimuthal velocity Vϕ.

The upper panel shows all of the LAMOST OCs in our catalogue. It is apparent that the clusters’ Z distances increase towards the outer disc, indicating that OCs can also probe the disc flare. The flare clusters with similar Vϕ to clusters close to the Galactic plane show an asymmetric structure from the north side to the south side of the Galactic plane. Two sequences of clusters are apparent below the Galactic plane, similar to the south branch and the main branch identified by Xu et al. (2020), but less extended. Xu et al. (2020) simulated the disc evolution after the impact of Sgr dSph ∼500 Myr ago and conclude that the interaction between the Galactic disc and the dwarf galaxy contributed to the flare. To check if OCs are also under the influence of the Sgr dSph passage, we divided our LAMOST OCs into two categories: a young group with age < 500 Myr, representing clusters born after the last passage of Sgr dSph; and an old group with age > 500 Myr, representing clusters born before the passage. The two categories are shown in the middle and lower panels of Fig. 12, respectively. We expect clusters in the older category to travel in the Galaxy more like a solid body and behave similarly to field stars, while clusters in the young category should be less affected by the interaction directly because they were gas clouds, or even stars of a previous generation, when the impact took place. Indeed, clusters with age > 500 Myr shaped the OC flare boundary (see the lower panel of Fig. 12), while clusters with age < 500 Myr are more concentrated along the Galactic plane with a small Vϕ variation.

We do notice that according to radial migration simulations (see e.g. Minchev et al. 2012), secular evolution of the Galactic disc can increase the scale height in the outer disc and decrease the scale height in the inner region. Thus, it can also produce a flare structure. However, clusters move outwards by gaining angular momentum and azimuthal velocity Vϕ (Schönrich & Binney 2009) while the low Vϕ of the LAMOST OCs in the flare region are difficult to reconcile with radial migration. Considering their asymmetric distribution above and under the Galactic plane, we suggest that the disc perturbation introduced by the last impact of Sgr dSph contributed to the OC flare.

4.2. Connection with nearby molecular clouds

Young clusters are a link between the interstellar medium (ISM) and stellar evolution. Since most stars are born in stellar clusters, linking young clusters to their surrounding molecular clouds is of great help in understanding star formation on a Galactic scale.

Using Gaia DR2 distances and stellar optical and near-infrared photometry, Zucker et al. (2019) present a technique to determine distances of molecular clouds. Zucker et al. (2020) applied this method to ∼60 star-forming regions and molecular clouds within 2.5 kpc. To investigate the connections between these clouds and young clusters, we show in Fig. 13 the positions of the Zucker et al. (2020) clouds in the Cartesian coordinate X-Y plane (centred on the Sun) and over-plot the young clusters (age < 100 Myr) of the LAMOST OC sample. Many LAMOST young clusters show an overlap with the molecular clouds on the Galactic plane. With such a young age and a similar Galactic position, they are very likely associated with each other. Since it is easier to obtain metallicity from stars than from molecular clouds, a young cluster is a good probe to study the metallicity evolution of the surrounding molecular clouds. We colour-code clusters in Fig. 13 with [Fe/H] in the left panel and cluster age in the right panel. It is clear that even for clusters that share a similar position on the Galactic plane, their [Fe/H] and age differ. This indicates that there have been multiple star formation epochs during the past 100 Myr in the solar neighbourhood, with gas metallicity covering a Δ[Fe/H] ∼ 0.4 dex (i.e. about a 2.7 times difference in metallicity). Whether this difference is due to inhomogeneous mixing in the giant molecular cloud, or because of fast star formation and pollution in the past 100 Myr, requires further investigations.

thumbnail Fig. 13.

Positions of the age < 0.1 Gyr clusters on the Galactic plane, compared to the nearby molecular clouds from Zucker et al. (2020), Zucker (2019). The two dashed circles indicate a Galactic radius of 8.2 kpc and 10 kpc. The clusters are colour-coded with their median [Fe/H] value in the left panel, and with age in the right panel. Only clusters with more than three FLAG=12 stars are displayed.

To study the metallicity distribution along a sequence of nearby molecular clouds, we investigated OCs around the Radcliffe wave, which is a coherent gaseous wave-like structure in the solar neighbourhood reported by Alves et al. (2020) based on the molecular cloud distance catalogue (Zucker et al. 2019, 2020).

Following Alves et al. (2020), we selected 24 LAMOST young OCs along the Radcliffe wave. Table 7 lists the Galactic coordinates, Cartesian coordinate XYZ, age, and [Fe/H] of these clusters. They are located along the gas wave clouds in three-dimensional space. Figure 14 shows the spatial distribution of these clusters together with the Radcliffe wave clouds. In both the Cartesian X–Y (upper panel) and Y–Z (lower panel) frames, these clusters show a high spatial consistency with the wave clouds. The metallicities of these clusters, with the higher values in the more central part of the wave, cover a non-negligible range, from sub-solar to super-solar. The highest [Fe/H] is about 1.4 times the Solar metallicity, and the lowest [Fe/H] is about 0.6 times the Solar metallicity. The [Fe/H] distribution of these clusters could be considered as the metallicity distribution of the Radcliffe wave in the past star formation epochs.

thumbnail Fig. 14.

[Fe/H] distribution of clusters along the Radcliffe gas wave. Upper panel: a bird’s eye view and the lower panel is a side view along the Cartesian Y axis. Clusters are colour-coded with their [Fe/H] values.

Table 7.

Galactic coordinates, age, Cartesian coordinate XYZ (centred on the Sun), age, and [Fe/H] of the 24 very young LAMOST OCs along the Radcliffe wave.

5. Summary

Using high-quality OC membership based on Gaia data (Cantat-Gaudin et al. 2020) and spectroscopic results from LAMOST DR8, we obtained a LAMOST OC catalogue with 386 OCs. The RV and metallicity of these clusters were determined homogeneously. To our knowledge, this is the first RV determination for 44 clusters and the first spectroscopic [Fe/H] determination for 137 clusters. Among the clusters with newly obtained [Fe/H], 63 are based on at least three member stars with high-quality stellar parameter determinations.

The cluster parameter determinations were based on a Monte Carlo sampling method. Member stars used for cluster Vrad determinations are marked with flag=1, those used for both Vrad and [Fe/H] determinations are marked with FLAG=12. The final LAMOST OC parameter catalogue, together with a table of member stars with quality control flag and LAMOST stellar parameters, is available at CDS.

During the quality control process to select FLAG=12 member stars, we noticed a systematic issue of LAMOST cool main sequence stars. Both the surface gravity log(g) and iron abundance [Fe/H] were under-estimated for main sequence stars with Teff ⪅ 5000 K. The problem appears in all clusters that have member stars in this range. These stars were not considered in the cluster parameter determinations. We suggest taking this problem into account also in studies using LAMOST field stars.

Using LAMOST OCs as tracers, we further studied the Galactic metallicity distribution and the Galactic disc dynamical properties. We calculated the orbit of the clusters and traced their Galactic metallicity gradient evolution in the past 500 Myr. Since most of the clusters have a non-circular orbit around the Galactic centre, we suggest using the Z component of the angular momentum LZ or the total dynamical energy instead of the current Galactic radius RGC to describe the Galactic metallicity gradient. With these two forms of metallicity gradient, we find two flat metallicity trend tails for very old OCs > 4 Gyr and for OCs in the [2 Gyr, 4 Gyr] age range.

We also investigated the OC metallicity distribution in the IOM space. Most of the LAMOST OCs lie in a tight line following the Galactic disc dynamics, while some outlier clusters show different [Fe/H] compared to other clusters in similar IOM space.

LAMOST OCs can be used to trace the Galactic disc flare. The OC flare shows an asymmetric structure above and under the Galactic plane, similar to the field star flare in the literature using LAMOST K giants (Xu et al. 2020). Based on the morphology and Vϕ distribution of the flare, we suggest that the last impact of Sgr dSph contributed to the OC flare.

The very young LAMOST OCs (age < 100 Myr) were compared to nearby molecular clouds. These young clusters, which are spatially overlapping molecular clouds, can be used to probe star formation history and metallicity evolution of the molecular clouds. We find that even for young clusters with very similar Galactic position, their [Fe/H] and age could differ by a factor of three. This indicates possible inhomogeneous mixing in local ISM (see similar discussions in De Cia et al. 2021) or multiple star formations in OCs during short time scales. Using 24 very young clusters along the Radcliffe wave, we present the metallicity distributions of the gas wave. The nature of this metallicity distribution and its connection to the formation of the Radcliffe wave require further investigations.


2

For instance, the offset value is 3.78 km s−1 compared to SIMBAD literature results (Gao et al. 2015, LAMOST DR1), 4.54 km s−1 compared to APOGEE DR14 (Anguiano et al. 2018, LAMOST DR3), 4.9 km s−1 compared to GALAH DR2 (Zhong et al. 2020, LAMOST DR5), and 4.4 km s−1 compared to RAVE DR3 (Xiang et al. 2015, LAMOST LSP3 based on DR1 spectra).

3

This is not strictly true, as [Fe/H] has been shown to vary due to atomic diffusion. However, the variations are within about 0.1 dex, see for instance Bertelli Motta et al. (2018), Semenova et al. (2020) for two well studied OCs.

Acknowledgments

X.F. acknowledges the support of China Postdoctoral Science Foundation No. 2020M670023 and NSFC No. 12203100. X.F. and A.B. acknowledge funding from the Italian MIUR through PREMIALE 2016 MITiC. X.F. and H.Z. thanks the support of the National Key R&D Program of China No. 2019YFA0405500 and the National Natural Science Foundation of China (NSFC) under grant No.11973001, 12090040, and 12090044. X.F., Z.Y.Z., C.L., Y.C., and J.Z. acknowledge the science research grants from the China Manned Space Project with NO.CMS-CSST-2021-A08. J.Z. thanks the support of NSFC No. 12073060. Y.C. acknowledges the support of NSFC under grant No. 12003001. L.L. thanks the support of the UCAS Joint PHD Training Program. L.C. and J.Z. acknowledges the support of NSFC through grants 12090040and12090042. Z.Y.Z. acknowledges the support of NSFC under grants No. 12041305, 12173016, and the Program for Innovative Talents, Entrepreneur in Jiangsu. This work benefited from the International Space Science Institute (ISSI/ISSI-BJ) in Bern and Beijing, thanks to the funding of the team Chemical abundances in the ISM: the litmus test of stellar IMF variations in galaxies across cosmic time (P.I.D. Romano and Z.-Y. Zhang). Guoshoujing Telescope (the Large Sky Area Multi-Object Fiber Spectroscopic Telescope LAMOST) is a National Major Scientific Project built by the Chinese Academy of Sciences. Funding for the project has been provided by the National Development and Reform Commission. LAMOST is operated and managed by the National Astronomical Observatories, Chinese Academy of Sciences. This work has made use of data from the European Space Agency (ESA) mission Gaia (https://www.cosmos.esa.int/gaia), processed by the Gaia Data Processing and Analysis Consortium (DPAC, https://www.cosmos.esa.int/web/gaia/dpac/consortium). Funding for the DPAC has been provided by national institutions, in particular the institutions participating in the Gaia Multilateral Agreement. This research has made use of the TOPCAT catalogue handling and plotting tool (Taylor 2005, 2017); of the Simbad database and the VizieR catalogue access tool, CDS, Strasbourg, France (Ochsenbein et al. 2000); of NASA’s Astrophysics Data System, and is supported by High-performance Computing Platform of Peking University.

References

  1. Alonso-Santiago, J., Frasca, A., Catanzaro, G., et al. 2021, A&A, 656, A149 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  2. Alves, J., Zucker, C., Goodman, A. A., et al. 2020, Nature, 578, 237 [Google Scholar]
  3. Andreuzzi, G., Bragaglia, A., Tosi, M., & Marconi, G. 2011, MNRAS, 412, 1265 [NASA ADS] [Google Scholar]
  4. Anguiano, B., Majewski, S. R., Allende Prieto, C., et al. 2018, A&A, 620, A76 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  5. Beccari, G., Boffin, H. M. J., Jerabkova, T., et al. 2018, MNRAS, 481, L11 [Google Scholar]
  6. Bertelli Motta, C., Pasquali, A., Richer, J., et al. 2018, MNRAS, 478, 425 [NASA ADS] [CrossRef] [Google Scholar]
  7. Bossini, D., Vallenari, A., Bragaglia, A., et al. 2019, A&A, 623, A108 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  8. Bragaglia, A. 2018, in Astrometry and Astrophysics in the Gaia Sky, eds. A. Recio-Blanco, P. de Laverny, A. G. A. Brown, & T. Prusti, 330, 119 [NASA ADS] [Google Scholar]
  9. Bragaglia, A., & Tosi, M. 2006, AJ, 131, 1544 [Google Scholar]
  10. Bragaglia, A., Fu, X., Mucciarelli, A., Andreuzzi, G., & Donati, P. 2018, A&A, 619, A176 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  11. Bragaglia, A., Alfaro, E. J., Flaccomio, E., et al. 2022, A&A, 659, A200 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  12. Bressan, A., Marigo, P., Girardi, L., et al. 2012, MNRAS, 427, 127 [Google Scholar]
  13. Buder, S., Asplund, M., Duong, L., et al. 2018, MNRAS, 478, 4513 [Google Scholar]
  14. Buder, S., Sharma, S., Kos, J., et al. 2021, MNRAS, 506, 150 [NASA ADS] [CrossRef] [Google Scholar]
  15. Cantat-Gaudin, T., & Anders, F. 2020, A&A, 633, A99 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  16. Cantat-Gaudin, T., Vallenari, A., Sordo, R., et al. 2018a, A&A, 615, A49 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  17. Cantat-Gaudin, T., Jordi, C., Vallenari, A., et al. 2018b, A&A, 618, A93 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  18. Cantat-Gaudin, T., Krone-Martins, A., Sedaghat, N., et al. 2019, A&A, 624, A126 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  19. Cantat-Gaudin, T., Anders, F., Castro-Ginard, A., et al. 2020, A&A, 640, A1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  20. Carraro, G., & Bensby, T. 2009, MNRAS, 397, L106 [NASA ADS] [CrossRef] [Google Scholar]
  21. Carrera, R., Bragaglia, A., Cantat-Gaudin, T., et al. 2019, A&A, 623, A80 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  22. Casali, G., Magrini, L., Frasca, A., et al. 2020a, A&A, 643, A12 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  23. Casali, G., Spina, L., Magrini, L., et al. 2020b, A&A, 639, A127 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  24. Casamiquela, L., Carrera, R., Jordi, C., et al. 2016, MNRAS, 458, 3150 [NASA ADS] [CrossRef] [Google Scholar]
  25. Casamiquela, L., Carrera, R., Blanco-Cuaresma, S., et al. 2017, MNRAS, 470, 4363 [NASA ADS] [CrossRef] [Google Scholar]
  26. Casamiquela, L., Carrera, R., Balaguer-Núñez, L., et al. 2018, A&A, 610, A66 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  27. Castelli, F., & Kurucz, R. L. 2003, Modelling of Stellar Atmospheres, eds. N. Piskunov, W. W. Weiss, & D. F. Gray, 210, A20 [Google Scholar]
  28. Castro-Ginard, A., Jordi, C., Luri, X., et al. 2018, A&A, 618, A59 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  29. Castro-Ginard, A., Jordi, C., Luri, X., Cantat-Gaudin, T., & Balaguer-Núñez, L. 2019, A&A, 627, A35 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  30. Castro-Ginard, A., Jordi, C., Luri, X., et al. 2020, A&A, 635, A45 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  31. Castro-Ginard, A., Jordi, C., Luri, X., et al. 2022, A&A, 661, A118 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  32. Chen, Y., Girardi, L., Bressan, A., et al. 2014, MNRAS, 444, 2525 [Google Scholar]
  33. Chen, Y., Bressan, A., Girardi, L., et al. 2015, MNRAS, 452, 1068 [Google Scholar]
  34. Choi, J., Conroy, C., Ting, Y.-S., et al. 2018, ApJ, 863, 65 [NASA ADS] [CrossRef] [Google Scholar]
  35. Conrad, C., Scholz, R. D., Kharchenko, N. V., et al. 2014, A&A, 562, A54 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  36. Cropper, M., Katz, D., Sartoretti, P., et al. 2018, A&A, 616, A5 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  37. Cui, X.-Q., Zhao, Y.-H., Chu, Y.-Q., et al. 2012, Res. Astron. Astrophys., 12, 1197 [Google Scholar]
  38. Dalton, G., Trager, S. C., & Abrams, D. C. 2012, in Ground-based and Airborne Instrumentation for Astronomy IV, eds. I. S. McLean, S. K. Ramsay, & H. Takami, SPIE Conf. Ser., 8446, 84460P [NASA ADS] [Google Scholar]
  39. De Cia, A., Jenkins, E. B., Fox, A. J., et al. 2021, Nature, 597, 206 [NASA ADS] [CrossRef] [Google Scholar]
  40. de Jong, R. S., Agertz, O., Berbel, A. A., et al. 2019, The Messenger, 175, 3 [NASA ADS] [Google Scholar]
  41. De Silva, G. M., Freeman, K. C., Bland-Hawthorn, J., et al. 2015, MNRAS, 449, 2604 [NASA ADS] [CrossRef] [Google Scholar]
  42. Dehnen, W., & Binney, J. J. 1998, MNRAS, 298, 387 [Google Scholar]
  43. Deng, L.-C., Newberg, H. J., Liu, C., et al. 2012, Res. Astron. Astrophys., 12, 735 [Google Scholar]
  44. Dias, W. S., Alessi, B. S., Moitinho, A., & Lépine, J. R. D. 2002, A&A, 389, 871 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  45. Dias, W. S., Monteiro, H., Lépine, J. R. D., et al. 2018, MNRAS, 481, 3887 [NASA ADS] [CrossRef] [Google Scholar]
  46. Dias, W. S., Monteiro, H., Moitinho, A., et al. 2021, MNRAS, 504, 356 [NASA ADS] [CrossRef] [Google Scholar]
  47. Donati, P., Bragaglia, A., Carretta, E., et al. 2015, MNRAS, 453, 4185 [NASA ADS] [CrossRef] [Google Scholar]
  48. Donor, J., Frinchaboy, P. M., Cunha, K., et al. 2018, AJ, 156, 142 [NASA ADS] [CrossRef] [Google Scholar]
  49. Donor, J., Frinchaboy, P. M., Cunha, K., et al. 2020, AJ, 159, 199 [NASA ADS] [CrossRef] [Google Scholar]
  50. Ferreira, F. A., Santos, J. F. C., Corradi, W. J. B., Maia, F. F. S., & Angelo, M. S. 2019, MNRAS, 483, 5508 [Google Scholar]
  51. Foreman-Mackey, D., Hogg, D. W., Lang, D., & Goodman, J. 2013, PASP, 125, 306 [Google Scholar]
  52. Friel, E. D., Janes, K. A., Tavarez, M., et al. 2002, AJ, 124, 2693 [NASA ADS] [CrossRef] [Google Scholar]
  53. Gaia Collaboration (van Leeuwen, F., et al.) 2017, A&A, 601, A19 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  54. Gaia Collaboration (Babusiaux, C., et al.) 2018a, A&A, 616, A10 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  55. Gaia Collaboration (Brown, A. G. A., et al.) 2018b, A&A, 616, A1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  56. Gaia Collaboration (Brown, A. G. A., et al.) 2021, A&A, 649, A1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  57. Gaia Collaboration (Vallenari, A., et al.) 2022, A&A, in press, https://doi.org/10.1051/0004-6361/202243511 [Google Scholar]
  58. Gao, H., Zhang, H.-W., Xiang, M.-S., et al. 2015, Res. Astron. Astrophys., 15, 2204 [CrossRef] [Google Scholar]
  59. Gilmore, G., Randich, S., Asplund, M., et al. 2012, The Messenger, 147, 25 [NASA ADS] [Google Scholar]
  60. Grevesse, N., & Sauval, A. J. 1998, Space Sci. Rev., 85, 161 [Google Scholar]
  61. Helmi, A. 2020, ARA&A, 58, 205 [Google Scholar]
  62. Jackson, R. J., Jeffries, R. D., Wright, N. J., et al. 2022, MNRAS, 509, 1664 [Google Scholar]
  63. Jofré, P., Heiter, U., & Soubiran, C. 2019, ARA&A, 57, 571 [Google Scholar]
  64. Katz, D., Sartoretti, P., Cropper, M., et al. 2019, A&A, 622, A205 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  65. Kharchenko, N. V., Piskunov, A. E., Schilbach, E., Röser, S., & Scholz, R. D. 2013, A&A, 558, A53 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  66. Kos, J., de Silva, G., Buder, S., et al. 2018, MNRAS, 480, 5242 [Google Scholar]
  67. Lépine, J. R. D., Mishurov, Y. N., & Dedikov, S. Y. 2001, ApJ, 546, 234 [CrossRef] [Google Scholar]
  68. Lépine, J. R. D., Cruz, P., Scarano, S. J., et al. 2011, MNRAS, 417, 698 [CrossRef] [Google Scholar]
  69. Li, L., & Shao, Z. 2022, ApJ, 930, 44 [NASA ADS] [CrossRef] [Google Scholar]
  70. Liu, L., & Pang, X. 2019, ApJS, 245, 32 [CrossRef] [Google Scholar]
  71. Luo, A. L., Zhao, Y.-H., Zhao, G., et al. 2015, Res. Astron. Astrophys., 15, 1095 [Google Scholar]
  72. Magrini, L., Randich, S., Kordopatis, G., et al. 2017, A&A, 603, A2 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  73. Magrini, L., Spina, L., Randich, S., et al. 2018, A&A, 617, A106 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  74. Magrini, L., Lagarde, N., Charbonnel, C., et al. 2021, A&A, 651, A84 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  75. Majewski, S. R., Schiavon, R. P., Frinchaboy, P. M., et al. 2017, AJ, 154, 94 [Google Scholar]
  76. McMillan, P. J. 2017, MNRAS, 465, 76 [Google Scholar]
  77. Miglio, A., Chiappini, C., Mackereth, J. T., et al. 2021, A&A, 645, A85 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  78. Minchev, I., Famaey, B., Quillen, A. C., et al. 2012, A&A, 548, A127 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  79. Netopil, M., Paunzen, E., Heiter, U., & Soubiran, C. 2016, A&A, 585, A150 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  80. Netopil, M., Oralhan, İ. A., Çakmak, H., Michel, R., & Karataş, Y. 2022, MNRAS, 509, 421 [Google Scholar]
  81. Ochsenbein, F., Bauer, P., & Marcout, J. 2000, A&AS, 143, 23 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  82. Origlia, L., Dalessandro, E., Sanna, N., et al. 2019, A&A, 629, A117 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  83. Randich, S., Tognelli, E., Jackson, R., et al. 2018, A&A, 612, A99 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  84. Randich, S., Gilmore, G., Magrini, L., et al. 2022, A&A, 666, A121 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  85. Reddy, A. B. S., & Lambert, D. L. 2019, MNRAS, 485, 3623 [NASA ADS] [CrossRef] [Google Scholar]
  86. Reddy, A. B. S., Lambert, D. L., & Giridhar, S. 2016, MNRAS, 463, 4366 [NASA ADS] [CrossRef] [Google Scholar]
  87. Rodrigues, T. S., Bossini, D., Miglio, A., et al. 2017, MNRAS, 467, 1433 [NASA ADS] [Google Scholar]
  88. Sartoretti, P., Katz, D., Cropper, M., et al. 2018, A&A, 616, A6 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  89. Scarano, S., & Lépine, J. R. D. 2013, MNRAS, 428, 625 [NASA ADS] [CrossRef] [Google Scholar]
  90. Schönrich, R., & Binney, J. 2009, MNRAS, 396, 203 [Google Scholar]
  91. Schönrich, R., Binney, J., & Dehnen, W. 2010, MNRAS, 403, 1829 [Google Scholar]
  92. Semenova, E., Bergemann, M., Deal, M., et al. 2020, A&A, 643, A164 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  93. Smiljanic, R., Donati, P., Bragaglia, A., Lemasle, B., & Romano, D. 2018, A&A, 616, A112 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  94. Soubiran, C., Cantat-Gaudin, T., Romero-Gómez, M., et al. 2018, A&A, 619, A155 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  95. Spina, L., Randich, S., Magrini, L., et al. 2017, A&A, 601, A70 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  96. Spina, L., Ting, Y. S., De Silva, G. M., et al. 2021, MNRAS, 503, 3279 [CrossRef] [Google Scholar]
  97. Spina, L., Magrini, L., & Cunha, K. 2022, Universe, 8, 87 [NASA ADS] [CrossRef] [Google Scholar]
  98. Tang, J., Bressan, A., Rosenfield, P., et al. 2014, MNRAS, 445, 4287 [NASA ADS] [CrossRef] [Google Scholar]
  99. Tarricq, Y., Soubiran, C., Casamiquela, L., et al. 2021, A&A, 647, A19 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  100. Taylor, M. 2017, Informatics, 4, 18 [CrossRef] [Google Scholar]
  101. Taylor, M. B. 2005, ASP Conf. Ser., 347, 29 [Google Scholar]
  102. Thomas, G. F., Laporte, C. F. P., McConnachie, A. W., et al. 2019, MNRAS, 483, 3119 [Google Scholar]
  103. Tsantaki, M., Pancino, E., Marrese, P., et al. 2022, A&A, 659, A95 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  104. Wan, J.-C., Liu, C., & Deng, L.-C. 2017, Res. Astron. Astrophys., 17, 079 [CrossRef] [Google Scholar]
  105. Wang, H.-F., Liu, C., Xu, Y., Wan, J.-C., & Deng, L. 2018, MNRAS, 478, 3367 [NASA ADS] [CrossRef] [Google Scholar]
  106. Wu, Y., Du, B., Luo, A., Zhao, Y., & Yuan, H. 2014, in Statistical Challenges in 21st Century Cosmology, eds. A. Heavens, J. L. Starck, & A. Krone-Martins, 306, 340 [NASA ADS] [Google Scholar]
  107. Xiang, M. S., Liu, X. W., Yuan, H. B., et al. 2015, MNRAS, 448, 822 [NASA ADS] [CrossRef] [Google Scholar]
  108. Xiang, M.-S., Rix, H.-W., Ting, Y.-S., et al. 2022, A&A, 662, A66 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  109. Xu, Y., Liu, C., Tian, H., et al. 2020, ApJ, 905, 6 [NASA ADS] [CrossRef] [Google Scholar]
  110. Yalyalieva, L. N., Chemel, A. A., Glushkova, E. V., Dambis, A. K., & Klinichev, A. D. 2018, Astrophys. Bull., 73, 335 [NASA ADS] [CrossRef] [Google Scholar]
  111. Yen, S. X., Reffert, S., Schilbach, E., et al. 2018, A&A, 615, A12 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  112. Zhang, R., Lucatello, S., Bragaglia, A., et al. 2021, A&A, 654, A77 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  113. Zhao, G., Zhao, Y.-H., Chu, Y.-Q., Jing, Y.-P., & Deng, L.-C. 2012, Res. Astron. Astrophys., 12, 723 [Google Scholar]
  114. Zhong, J., Chen, L., Wu, D., et al. 2020, A&A, 640, A127 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  115. Zucker, C. 2019, Table A.1, Harvard Dataverse, V.2, https://doi.org/10.7910/DVN/07L7YZ [Google Scholar]
  116. Zucker, C., Speagle, J. S., Schlafly, E. F., et al. 2019, ApJ, 879, 125 [Google Scholar]
  117. Zucker, C., Speagle, J. S., Schlafly, E. F., et al. 2020, A&A, 633, A51 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]

All Tables

Table 1.

Radial velocity of 342 OCs with known Vrad in the literature.

Table 2.

Newly obtained cluster radial velocity of 44 OCs.

Table 3.

218 clusters with [Fe/H] values in the literature.

Table 4.

Newly obtained iron abundances for 63 OCs with at least three FLAG=12 member stars.

Table 5.

Newly obtained iron abundances for 74 OCs but with fewer than three FLAG=12 member stars.

Table 6.

Second-order polynomial MCMC fitting parameters of the clusters’ [Fe/H] as a function of the Galactic radius (RGC), the Z component of the Galactic angular momentum (LZ), and the total Galactic energy (En) for six age groups.

Table 7.

Galactic coordinates, age, Cartesian coordinate XYZ (centred on the Sun), age, and [Fe/H] of the 24 very young LAMOST OCs along the Radcliffe wave.

All Figures

thumbnail Fig. 1.

Sky coverage of the LAMOST open cluster (OC) members (red dots) and all the Cantat-Gaudin et al. (2020) OC member stars (grey dots).

In the text
thumbnail Fig. 2.

NGC 2548 as an example of the RV selection. The RVs of all the crossed-matched member stars in the cluster are plotted assuming a normal distribution. The light grey colour curves denote stars we discarded for further analysis (FLAG = 0), while the dark grey ones are the selected stars. The two vertical dashed lines mark the 2-σ selection criterion. The Monte Carlo model result of all the selected stars is illustrated with the blue Gaussian curve. Their mean RV value (Vrad, mean), the 1-σ RV dispersion (σVrad), as well as the number of member stars before (Nstar) and after the RV selection (Nstar, RV), are shown in the legend.

In the text
thumbnail Fig. 3.

Histogram of the RV difference between Gaia DR3 measurements and LAMOST measurements for all the RV selected members. The median value of the difference ΔRV = RVGaia DR3 − RVLAMOST is 4.66 km s−1.

In the text
thumbnail Fig. 4.

Kiel diagram (upper panel, colour-coded with [Fe/H]) and [Fe/H] as a function of Teff (lower panel, colour-coded with log(g)) of cluster NGC 2632. For reference, a PARSEC theoretical isochrone of Z = 0.013, log(age) = 9.3 is plotted in the Kiel diagram. The line-shaded region marks Teff < 4000 K, and the grey part marks 5000 K < Teff < 4000 K.

In the text
thumbnail Fig. 5.

Colour-magnitude diagrams of Melotte 22 and NGC 2682 as examples of the clusters with [Fe/H] in the literature. The filled grey circles are high-quality member stars from Cantat-Gaudin et al. (2020), and open circles are all matched stars with LAMOST RV measurements. Right panel: stars good for the cluster Vrad determination (FLAG=1) are colour-coded with their RV values. Left panel: the colour-filled dots are stars good for [Fe/H] determination (FLAG=12), and the black dots are stars without [Fe/H] determination or with discarded [Fe/H]. The dashed curve in each sub-figure is the PARSEC isochrone with the Cantat-Gaudin et al. (2020) cluster parameters.

In the text
thumbnail Fig. 6.

Comparison of our cluster Vrad to results in the literature. Left, middle, and right panels: comparison to 308, 185, and 226 OCs of Tarricq et al. (2021), Dias et al. (2021), and Zhong et al. (2020), respectively. The histograms of the cluster RV differences are shown as insets in each panel.

In the text
thumbnail Fig. 7.

Comparison of 82 clusters’ [Fe/H] to literature results with high-resolution spectra (HRS) in the left panel, and of 206 clusters to results from Zhong et al. (2020). The X symbols mark [Fe/H] values with less than three member stars, namely 38 clusters on the left (from HRS) and 75 clusters on the right (from Zhong et al. 2020). The open squares mark clusters with less than three FLAG=12 member stars from this work, which are 23 clusters on the left panel and 92 clusters on the right panel.

In the text
thumbnail Fig. 8.

Clusters’ median [Fe/H] as a function of the Galactocentric radius RGC for all the LAMOST clusters with at least three FLAG=12 stars. Different panels show different look-back times according to the orbit calculations. Clusters are divided into six groups according to age from the Cantat-Gaudin et al. (2020) catalogue. The solid curves are the second-order polynomial MCMC fittings of each age group.

In the text
thumbnail Fig. 9.

Clusters’ eccentricity as a function of their Galactic apocentric distance. The colour code is the log value of the cluster age in the left panel, and clusters’ [Fe/H] in the right panel. Only clusters with multiple members are displayed in the left panel, while those with at least three FLAG=12 stars are displayed in the right panel.

In the text
thumbnail Fig. 10.

Clusters’ median [Fe/H] as a function of their Z component of the Galactic angular momentum (left panel) and the total Galactic energy (right panel) for all of the LAMOST clusters with no less than three FLAG=12 stars. The solid curves are the second-order polynomial MCMC fittings of each age group.

In the text
thumbnail Fig. 11.

Integral of motion (IOM) of LAMOST OCs with multiple members (Nflag = 1≥ 2). The [Fe/H] colour-code is applied to clusters with at least three FLAG=12 members. The empty symbols are clusters with less than three FLAG=12 members.

In the text
thumbnail Fig. 12.

Present-day location of LAMOST OCs on the ZGC − RGC plane. Upper panel: all LAMOST OCs, while the middle and lower panels show clusters born after the last passage of Sgr dSph (age < 500 Myr), and clusters born before the passage (age > 500 Myr), respectively. Clusters are colour-coded with their azimuthal velocity Vϕ.

In the text
thumbnail Fig. 13.

Positions of the age < 0.1 Gyr clusters on the Galactic plane, compared to the nearby molecular clouds from Zucker et al. (2020), Zucker (2019). The two dashed circles indicate a Galactic radius of 8.2 kpc and 10 kpc. The clusters are colour-coded with their median [Fe/H] value in the left panel, and with age in the right panel. Only clusters with more than three FLAG=12 stars are displayed.

In the text
thumbnail Fig. 14.

[Fe/H] distribution of clusters along the Radcliffe gas wave. Upper panel: a bird’s eye view and the lower panel is a side view along the Cartesian Y axis. Clusters are colour-coded with their [Fe/H] values.

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.