Issue |
A&A
Volume 692, December 2024
|
|
---|---|---|
Article Number | A212 | |
Number of page(s) | 17 | |
Section | Galactic structure, stellar clusters and populations | |
DOI | https://doi.org/10.1051/0004-6361/202450726 | |
Published online | 13 December 2024 |
When LAMOST meets Gaia DR3
Exploring the metallicity of open clusters
1
Yunnan Observatories, Chinese Academy of Sciences,
Kunming
650216,
PR China
2
School of Physical Science and Technology, Xinjiang University,
Urumqi
830046,
PR China
3
INAF-Osservatorio Astronomico di Padova,
vicolo Osservatorio 5,
35122
Padova,
Italy
4
Dipartimento di Fisica e Astronomia “Galileo Galilei”, Universita di Padova,
Vicolo Osservatorio 3,
35122
Padova,
Italy
5
Department of Physics, Stellenbosch University,
Matieland 7602,
Stellenbosch,
South Africa
6
National Institute for Theoretical and Computational Sciences (NITheCS),
Stellenbosch,
South Africa
7
Astrophysics Department, American Museum of Natural History,
Central Park West,
New York,
NY
10024,
USA
8
Institute for Frontiers in Astronomy and Astrophysics, Beijing Normal University,
Beijing
102206,
PR China
9
School of Physics and Astronomy, Beijing Normal University,
Beijing
100875,
PR China
10
Institute for Advanced Studies, Technische Universität München,
Lichtenbergstraße 2 a,
85748
Garching bei München,
Germany
11
Purple Mountain Observatory, Chinese Academy of Sciences,
Nanjing
210023,
China
★ Corresponding author; ruyuan.zhang@ynao.ac.cn, yl4331@columbia.edu, guosufen@ynao.ac.cn
Received:
15
May
2024
Accepted:
9
October
2024
Context. Open clusters (OCs) are excellent probes as their age and abundance can be tightly constrained, allowing us to explore the distribution of metallicity and composition across the disk of the Milky Way. By conducting a comprehensive analysis of the metallicity of OCs, we can obtain valuable information about the history of their chemical enrichment. Moreover, by observing stars in different regions of the Milky Way, we can identify significant spatial structures in their chemical composition and abundance. This enables us to understand stellar birth radii through chemical tagging. Nevertheless, it remains challenging to infer the original positions of OCs using current data alone.
Aims. The aim of this study is to investigate the distribution of metallicity in the solar neighborhood using a large dataset from Gaia DR3 combined with LAMOST spectra. With accurate ages and metallicity measurements, we can determine birth radii for the stars and attempt to understand their migration pattern.
Methods. We chose a total of 1131 OCs within 3 kpc of the Sun from the Gaia DR3 and LAMOST DR8 low-resolution spectral database (R=1800). We used an artificial neural network to correct the LAMOST data by incorporating high-resolution spectral data from GALAH DR3 (R=28 000). The average metallicity of the OCs was determined based on the reliable [Fe/H] values for their members. We then examined the distribution of metallicity across different regions within the Galaxy and inferred birth radii of the OCs from their age and metallicity.
Results. The correction method presented here can partially eliminate the systematic offset for LAMOST data. We discuss the metal- licity trend as a function of Galactocentric distance and the guiding radii. We also compare these observational results with those from chemo-dynamic simulations. Values derived from observational metallicity data are slightly lower than predicted values when the uncertainties are not considered. However, the metallicity gradients are consistent with previous calculations. Finally, we investigated the birthplace of OCs and find hints that the majority of OCs near the Sun have migrated from the outer Galactic disk.
Key words: stars: abundances / stars: evolution / galaxies: star clusters: general
© The Authors 2024
Open 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
Galactic archaeology seeks to unravel the past of the Milky Way and to construct a timeline of its evolution. By analyzing observational data for stars, we aim to understand the Galaxy’s environment during its formation and subsequent evolution. The kinematics, age, and chemical composition of the stars allow us to derive tight constraints on their birthplace. Stars inherit their chemical composition from the interstellar medium (ISM) where they form, including elements produced through various processes such as supernovae, stellar winds, and neutron star mergers. This enables the identification of stars with similar chemical patterns as originating from the same molecular cloud at specific moments in time, a technique known as chemical tagging (see e.g., Freeman & Bland-Hawthorn 2002). While molecular clouds within the same galaxy may exhibit unique compositions, accurately tracing stars to their original ISM requires precise chemical analysis (De Silva et al. 2006). The method of weak chemical labeling proposed by Vanessa Hill (see e.g., Gilmore et al. 2022), combined with age and dynamic information (e.g., Wang et al. 2020b,c and references therein), offers a promising approach to exploring the original locations of stars, overcoming challenges in analyzing large samples.
Determining the precise ages of stars is vital for constructing a timeline of the formation of the Milky Way. However, age measurement, especially for older stars, poses challenges and entails greater uncertainties. In addition to traditional isochrone fitting, age determination relies on combining photometric data across multiple bands with spectrum analysis using sky survey data (Sanders & Das 2018; Queiroz et al. 2018; Kordopatis et al. 2023; Wang et al. 2024). Furthermore, a new method based on aster-oseismology provides a more reliable approach to determining stellar age (Chaplin & Miglio 2013; Chaplin et al. 2014, 2020; Anders et al. 2023). Another promising method is the use of “chemical clocks”, which rely on tracking the changes in element abundances within stars as a result of nuclear reactions occurring in their interiors. These reactions can either increase or decrease the abundances of certain elements as the star ages (Hayden et al. 2022; Viscasillas Vazquez et al. 2022; Moya et al. 2022). The integration of stellar age with chemical abundance could be used as a tool for establishing a comprehensive event calendar.
Open clusters (OCs) are stellar groups held together by gravity and formed from a single molecular cloud. They are predominantly young – typically less than 1 Gyr (Fujii & Portegies Zwart 2016) – with the oldest existing for around 8 Gyr before dissipation (Salaris et al. 2004; Kharchenko et al. 2013). The age, distance, metallicity, and composition of OCs are relatively easy to determine compared to those of field stars because OCs consist of stars with similar ages and chemical compositions. The spatial concentration of cluster members also facilitates accurate distance measurements using methods such as main sequence fitting and isochrone fitting. Furthermore, the stars in an OC formed within the same ISM. This simplifies the process of determining their metallicity and chemical composition through spectroscopic analysis. Mostly found in the Milky Way disk, OCs share metallicities of approximately the solar value, and are distributed across the galactic thin disk. Over time, OCs move away from their birth location via radial migration (Sellwood & Binney 2002). Studies reveal a correlation between the ages of OCs and their spatial distribution, with older clusters situated farther from the Galactic plane, influenced by vertical heating mechanisms driven by Galactic structures such as warps and minor mergers (see e.g., van den Bergh & McClure 1980; Friel 1995; Cantat-Gaudin et al. 2020; Martig et al. 2014; Mackereth et al. 2019; Ting & Rix 2019; Sharma et al. 2021). These dynamic processes, involving radial migration and heating, can elevate OCs to higher vertical positions, showcasing the intricate interplay between objects and the Galactic structure (Sellwood & Binney 2002; Bird et al. 2012; Minchev & Famaey 2010; Minchev et al. 2011, 2013).
Open clusters offer a means to study the chemical, structural, and dynamic evolution of the Galactic disk across various positions in the Galaxy, and the [Fe/H] trend is crucial for constraining chemical evolution models and for the construction of a timeline of Galactic events. Various studies (Reddy et al. 2016; Genovali et al. 2014; Casamiquela et al. 2019; Donor et al. 2020; Cantat-Gaudin et al. 2020; Netopil et al. 2022; Myers et al. 2022) have analyzed high-quality data for OC samples and stars in order to use these objects as tracers to detect the behavior and gradient of chemical composition across galactic space. The metallicity gradient derived from OCs typically has a value of around −0.060 dex kpc−1, but flattens toward the outer disk (Friel et al. 2010; Donor et al. 2020; Myers et al. 2022). Myers et al. (2022) compared the metallicity trend observed in APOGEE DR17 with results predicted from chemo-dynamical models, discovering that observational results exhibit a systematic bias for older OCs when compared to predictions. These authors suggested that this bias may be attributed to the migration of OCs.
The work of Minchev & Famaey (2010) and Minchev et al. (2018) focuses on mechanisms of disk migration. Many studies use large survey samples to further explore galactic disk migration. Chen & Zhao (2020) and Zhang et al. (2021b) selected samples of OCs based on the High-resolution spectra (HRS) data provided by Netopil et al. (2016), and found evidence of migration. These samples consist of three types of migration modes, and these authors conducted their analysis of migration based on the distribution of metallicity. Viscasillas Vázquez et al. (2023) investigated the interaction between field stars and OCs and galactic disk migration using Gaia-DR3 data.
Currently, a few thousand OCs have been identified using photometric data from large Sky surveys, and less than 10% of these have been the target of high-resolution spectroscopic observations prior to recent advancements. There are several major surveys – such as Gaia-ESO, APOGEE, and GALAH –1 that have collected information about OCs using high-resolution spectral data, although the sample sizes are limited. Gaia Data Release 3 (DR3) has provided information on approximately 5.6 million stars through its Radial Velocity Spectrometer (RVS) (Cropper et al. 2018; Katz et al. 2023). Using the General Stellar Parametriser spectroscopy module (GSP-Spec), stellar atmospheric parameters and composition abundance were determined (Recio-Blanco et al. 2023). Compared to Gaia DR2, parallax and proper motion accuracies have improved by 30% and 100%, respectively (Marton et al. 2023), and the Gaia DR3 target magnitude range spans from stars brighter than 13.6 in G band to faint stars of around 25 in G band (Riello et al. 2021; Gaia Collaboration 2023)
The Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST)2 is located in Xinglong, China. It has a diameter of 4 meters and consists of two spectral modules. LAMOST is the only instrument capable of obtaining data for the outer-disk OCs. Thanks to its powerful capability, LAMOST has observed millions of stars, covering a wide range of disk positions (Deng et al. 2012). This dataset includes the latest Gaia DR3 photometric data and features an extensive number of OC samples with metallicity measurements, significantly expanding the sample size for outer-disk OCs.
This paper is organized as follows: Section 2 presents the sample selection and the reference data, providing a detailed comparison of the radial velocity (RV) of the samples and the correction of metallicity and α-elements using an ANN. In Sect. 3, we discuss the differences between our results and those gathered from the literature. Section 4 examines the distribution of metallicity across various Galactic spatial regions and the migration patterns of OCs. Finally, we summarize our findings in Sect. 5.
2 Data and samples
Gaia DR3 data can be used to improve the identification of OC membership and can provide more accurate distances for distant clusters. In their catalog, Hunt & Reffert (2023) reported over 7000 clusters and stellar associations, including 2378 newly identified clusters and 134 globular clusters spanning the entire sky. Additionally, Cavallo et al. (2024) developed an ANN model using a new feature-extraction technique called the QuadTree algorithm (Schiappacasse-Ulloa et al. 2023). This model was used to re-estimate the parameters of OCs based on Hunt’s catalog. However, the determination of metallicity in this study relied on photometric data, and the metallicities obtained are hence of modest accuracy, as shown in Fig. 7 of Cavallo et al. (2024). Although Gaia has expanded the number of OCs with available chemical abundance data, such information is limited to objects within a relatively bright magnitude limit within which Gaia RVS can effectively collect information. In their work, Viscasillas Vázquez et al. (2023) used the Gaia DR3 samples with chemical abundance data from GSP-Spec. These authors selected clusters and field stars that covered the Galactocentric range of 7.5–9 kpc and investigated the radial migration of objects. However, the availability of samples in the outer range of the disk is limited.
LAMOST is a highly effective telescope, and is able to observe approximately 4000 targets per exposure. The low-resolution module has a resolution of 1800 and covers the wavelength range of 3700–9000 Å. The medium-resolution module has a resolution of 7500. In those modules, LAMOST uses two cameras: the blue camera captures spectra from 4950 to 5350 Å, while the red camera captures spectra from 6300 to 6800 Å. LAMOST has already released its DR8 dataset for the international community, which includes low- and medium-resolution spectra for over 5 million objects (Liu et al. 2020). The LAMOST stellar parameter pipeline (LASP) analyzes LAM-OST’s low-resolution spectra and determines key atmospheric parameters. Spectral analysis involves extracting 1D spectra post-data reduction, including corrections and combining sub-exposures. It also involves detecting initial parameters using correlation function interpolation (CFI, Du et al. 2012) and deriving final parameters with ULySS; (Koleva et al. 2009; Wu et al. 2011a,b) by minimizing χ2,3. The details of the spectral template used to measure atmospheric parameters are summarized in the work of Bai et al. (2021).
Additionally, LAMOST has entered the phase of medium-resolution observations, leading to an increase of approximately 2 million stellar samples. However, at present, there are no publicly accessible pipelines specifically developed for the analysis of LAMOST medium-resolution spectra (Luo et al. 2015; Ren et al. 2021), and therefore we do not analyze these here. In the present study, a total of 1131 OCs were selected, which includes a total of 8640 stars. Of these, more than 1020 OCs have less than three stars, while the remaining OCs typically contain five stars.
Recently, Cavallo et al. (2024) undertook a new analysis of stellar associations made up of more than ten stars using the catalog compiled by Hunt & Reffert (2023). The authors focused on approximately 6300 OCs, and used an ANN to determine their fundamental parameters. These parameters were derived from Gaia and 2MASS data, including age, extinction, distance, and metallicity. We selected the samples for the present study by matching the positions of objects within 3 arcsec (Cui et al. 2012) between the LAMOST catalog and the catalog of Cavallo et al. (2024). Furthermore, we cross-checked the citation ID number from LAMOST with the Gaia DR3 source ID to ensure consistency. Of the identified stars, over 8000 were found to belong to approximately 1000 OCs listed in Cavallo’s catalog.
2.1 Radial velocity
The RVs employed in this section are derived from the low-resolution spectra (LRS) obtained from the LAMOST Data Release 8 catalog. LASP uses cross-correlation to identify spectral classes and determine redshifts or RVs. To tackle the challenges of detection, a low-order polynomial is used to align spectra. Currently, SDSS templates are being used; however, LAMOST templates will be updated after the first data are released (Luo et al. 2015). These velocities are influenced by the S/N and are determined through the optimal fit of χ2. However, it is important to note that RV measurements derived from LRS are affected by significant uncertainties, which warrant attention. There are noticeable differences between RVs obtained from LRS and those derived from high-resolution spectra (Fu et al. 2022). As the members of the OCs are identified from the literature, the process of identifying membership is not included in the present study. The comparison results are taken directly from the LASP DR8 catalog, where the uncertainty associated with the LAMOST RVs peaks around 5 km/s. We noted the radial-velocity scatter in a single OC and used the LASP outcome, considering different signal-to-noise ratios (S/Ns). The RV dispersion values obtained depend on S/N and spectral type, as indicated by Bai et al. (2021). The corresponding standard deviation ranges for FGK stars are approximately 7.67–13.18 km/s, 8.60–6.59 km/s, 5.30–6.19 km/s, and 3.93–5.17 km/s for S/N values of 5–10, 10–20, 20–40, and >40, respectively. The RV of an OC becomes more uniform as the S/N increases. The accuracy of RV measurements for individual stars relies on their S/N, with higher S/N measurements being considered more reliable. Stars with significant RV measurement discrepancies are flagged as potential binaries and will be excluded from further selection unless validated as single stars. Ultimately, the RV of an OC is calculated by averaging the RVs of its individual members.
The RVs used for comparison are derived from high-resolution spectra obtained through various projects, including the Gaia-ESO Survey (GES), Open Clusters Chemical Abundances from Spanish Observatories (OCCASO), GALactic Archaeology with HERMES (GALAH), Stellar Populations and Abundances (SPA), and One Star to Tag Them All (OSTTA) datasets (Casali et al. 2019; Casamiquela et al. 2016; Spina et al. 2021; Frasca et al. 2019; Casali et al. 2020; D’Orazi et al. 2020; Zhang et al. 2021a; Carrera et al. 2022). Additionally, we include the sample studied by Donor et al. (2020) using APOGEE DR16 in our RV comparison. In this study, to ensure data quality and an adequate number of samples for comparison, we imposed a minimum S/N of 30 in the final data cut. Table 1 presents the radial velocities for the currently considered open clusters. Figure 1 displays the RV distribution per OC and the offset between the HRS result and our work. On average, there is a difference of approximately 5.59 km s−1, with a scatter of 14.01 km s−1. Tsantaki et al. (2022) assessed the differences in RV measurements between LAMOST and other surveys (APOGEE, GALAH, RAVE), finding that the RV bias between LAMOST and high-resolution spectra is between 4 and 5 km s−1 for spectra with a S/N of greater than 30. This value is slightly lower than but comparable to our findings.
We find a greater scatter between the RVs of the OCs studied here compared with other studies; this is verified again in this section. For example, the difference in RV between FSR 0951 and Berkeley 32 exceeds 50 km/s. The present study only includes two stars from FSR 0951, which increases the uncertainty in the RV measurements. Berkeley 23, an old OC, has six members, but only two stars have RV measurements that can be compared with those from the literature. Most of the large offsets in Fig. 1 correspond to clusters that contain only one or two member stars, such as Collinder 107, Alessi 62, Berkeley 19, Berkeley 31, Berkeley 17, SFR 0667, ASCC 123, and UPK 422.
This highlights the importance of considering the number of stars within OCs when making comparisons.
Gaia DR3 provides the stellar RV for more than 300 000 stars, which include objects found in 595 of the OCs in our sample. Figure 2 shows the RVs of Gaia DR3 compared with the RVs for the sample of OCs studied in the present paper. The systematic offset is around 6.075 km s−1, but the significant dispersion in the RV residuals between the two databases may be attributed to RV measurement errors and unverified binaries.
Radial velocities for the OCs studied here compared to those from the literature.
![]() |
Fig. 1 Difference in RVs between 101 cross-matched targets from previous works, including those by Casali et al. (2019), Casamiquela et al. (2016), Spina et al. (2021), Frasca et al. (2019), Casali et al. (2020), D’Orazi et al. (2020), Zhang et al. (2021a), and Carrera et al. (2022), and the current study is shown. The gray lines indicate the mean offset value of 5.598 km/s with a standard deviation of 14.013 km/s for this comparison. The x-axis represents the RV of the present study, incorporating the error for the overlapping sample. |
![]() |
Fig. 2 Difference distribution in RVs for 595 cross-matching targets between Gaia DR3 and the present study. The mean value of the offset is 6.075 with a standard deviation of 11.496 km s−1 as indicated by the gray lines in the plot. The x-axis represents the LAMOST RV and the error of the overlapped sample. The y-axis includes the RV uncertainties from Gaia DR3. |
2.2 Metallicity determination
This section focuses on controlling the quality of the data. The analysis presented here was carried out using the LAMOST official pipeline LASP, which has now published eight data releases. The LAMOST General catalog contains the atmospheric parameters (effective temperature, surface gravity, and metallicity) of F, G, and K stars. Currently, this catalog has published spectral data for over 6 million stars. To ensure high quality, as mentioned above, only the final results with a spectral S/N of greater than 30 were selected. However, the metallicity values show a systematic bias, especially for low-temperature dwarfs with an effective temperature of below 5000 K and surface gravity higher than 4.3 dex (Fu et al. 2022; Soubiran et al. 2022). Similar properties are derived for these stars when using the data-driven DD-Payne method (Xiang et al. 2019), which relies on observed spectra for training. The spectral analysis of cool stars presents larger uncertainties. This is because the analysis becomes more challenging as the stars get cooler, mainly due to blending (Casali et al. 2020), which leads to underestimation of metallicity in the LAMOST data. This could be inherited from the determination template mentioned in Sect. 2. To tackle this issue, we employed a calibration technique that involves adjusting the metallicity abundance. This was done using an empirical correction relationship derived from wide binary stars. Furthermore, the metallicity of cooler stars was calibrated by referencing their hotter companions (Niu et al. 2023). However, based on the current data, the average standard deviation of [Fe/H] for the chosen OCs is approximately 0.16 dex. Additionally, over 10% of the clusters exhibit a metallicity dispersion of greater than 0.25 dex, surpassing the desired chemical abundance precision of 0.05 dex (Freeman & Bland-Hawthorn 2002). In this section, we aim to correct the stellar metallicity using a similar method to that used by Spina et al. (2021). We use an ANN to build a model that corrects the iron abundance for individual stars.
In the remaining samples, Soubiran et al. (2022) evaluated the accuracy of the measurements of [Fe/H] as provided by multiple surveys, including APOGEE,GALAH, RAVE, LAMOST, SEGUE, and Gaia-ESO. These authors used samples of OCs to investigate internal uncertainties. They found that the residuals of [Fe/H] changed by between −0.3 and −0.15 dex within a temperature range from 4000 K to 6500 K. In this section, we correct the data from LAMOST using the data from the GALAH survey. The GALAH survey aims to investigate the origin of the Milky Way, including its dynamic processes and chemical composition. The survey provides 30 different chemical abundance measurements for around 600 000 stars (Buder et al. 2021).
For our training samples, we selected around 30 000 targets that were observed in both the LAMOST and GALAH surveys. When choosing these samples, we avoided stars that are too cool and for which the spectra are dominated by molecular lines (Jofré et al. 2019), and those that are too hot and lack metal lines for analysis. Our selection criteria align with those mentioned by Spina et al. (2021): the effective temperature must fall within the range of 4000 K and 7000 K, the surface gravity must be >0 dex, the metallicity must be between −1 and 0.5 dex, and the microturbulence should be smaller than 2.5 km/s. Additional criteria include | Teff,LAMOST-Teff, GALAH |< 150 K,| log gLAMOST-log gGALAH |< 0.3 dex, and |[Fe/H]LAMOST-[Fe/H]GALAH |< 0.1 dex.
We used our ANN to achieve the correction. An ANN is a computational model that takes inspiration from biological neural networks. The general structure of an ANN is made up of an input layer, several hidden layers, and an output layer. An ANN is trained on a data set to learn a mapping between the input and the output. In the training process, the hyperparameters (weight w and bias b) of the ANN are optimized by minimizing a loss function. For the correction here, we trained an ANN to learn a mapping between the LAMOST data and the GALAH data. To this end, we fed the ANN model with the LAMOST metallicity and [Fe/H] uncertainties, expecting to obtain the corresponding calibration values. We built and trained the ANN model using the code CoLFI4, which is a framework for likelihood-free inference with neural density estimators (Wang et al. 2020a, 2022, 2023b). The ANN model adopted here contains a hidden layer with 128 neurons (see Fig. 3). We used the rectified linear unit (ReLU, Nair & Hinton 2010) as the nonlinear function, and the least absolute deviation as the loss function.
The training set, validation set, and test set for the ANN model were divided using the traditional fraction of 6:2:2. The loss functions of the training and validation sets are depicted in Fig. 4, indicating the absence of overfitting. Next, we used the well-trained ANN model to correct the data from LAMOST. The corrected data are presented in Fig. 5, which shows a greater concentration around the relational function. The correction is particularly effective for stars with significant differences in [Fe/H] between the two databases.
2.3 Open cluster metallicity
The accuracy of OC parameters can be enhanced by examining the star that shows the greatest scatter. Furthermore, the uncertainties can be reduced by averaging the results obtained from multiple stars. We aim to estimate the metallicity of OCs using the method described in Sect. 2.2. For OCs with the highest precision, the range of [Fe/H] dispersion within OCs is reported to be between 0.04 dex and 0.08 dex (Jofré et al. 2019). However, it is important to note that the metallicity dispersion tends to increase with age. Therefore, caution should be exercised when using individual metallicities that exhibit a significant deviation from the mean value of OCs. To ensure data accuracy, we applied a data-cleaning process for each OC, excluding members with [Fe/H] values exceeding three standard deviations of the OC’s mean metallicity. The final metallicity for each OC was then determined using the averaged values of remaining stars post-cleaning
In order to validate the applicability of the ANN model described in Sect. 2.2 to the current sample, it is important to have a training sample that covers a similar range as the parameters that need to be corrected. The input table contains the stellar metallicity values ranging from −1 to 0.5 dex. Figure 6 displays the overlapping portion of the histogram between the model’s training data and the data that require correction. The plot demonstrates that the training set and the metallicity values to be corrected are comparable.
We verified the coverage of the input samples in the training set before eventually using the correction prediction. The final metallicity for each OC is checked by removing membership outliers, and even though they are in the multi-evolution phase, their final metallicity shows less dispersion. Figure 7 shows the member distribution of OCs on a color-magnitude diagrams (CMD) for six typical OCs in this study, with the color code indicating metallicity. The isochrone used is from the PARSEC V2.05 isochrone (Nguyen et al. 2022), and the OC parameters are based on Cavallo et al. (2024) as shown in Fig. 7. The members of the selected OCs span an age range from 130 Myr to 5.3 Gyr, and they are located in different evolution phases with a metallicity scatter of less than 0.2 dex. However, uncertainties and systematic errors in age measurements can cause a slight offset between targets and isochrones. We examined whether or not there are variations in the fit differences when different metallicity bins are used, as shown in Fig. 8 for Collinder 69 as an example. It was concluded that any potential systematic bias between the observational data and isochrones can be ignored.
![]() |
Fig. 3 Structure of the ANN used in this work. The number of neurons in the hidden layer is 128. |
![]() |
Fig. 4 Loss function of the training and validation sets. |
![]() |
Fig. 5 Comparison of Metallicity Measurements from LAMOST DR8 and GALAH. The right panel shows the iron relation for 30 000 targets. The y-axis represents the metallicity from LAMOST dr8, while the x-axis represents the corresponding value from GALAH. The left plot displays the metallicity residuals of LAMOST and GALAH, both before and after correction. The blue dots represent the original data selected from LAMOST DR8, while the orange dot represents the LAMOST result after correction. |
![]() |
Fig. 6 Histogram of the training set and the metallicity that needed to be corrected for the LASP result. The orange section represents the training data that were used to build the ANN model, while the blue section represents the metallicity that needed to be corrected for the selected samples. |
3 Comparison
3.1 Metallicity comparison with the literature for individual stars
This section presents the iron abundances of the OCs as determined using high-resolution spectra mainly from three datasets: GALAH, APOGEE DR17 (Abdurro’uf et al. 2022), and PASTEL Soubiran et al. (2016, 2022). The last of these three is a collection of atmospheric parameters from high-resolution spectra (R > 25 000) with S/N > 50. The PASTEL catalog comprises high-resolution data collected since the 1980s (Cayrel de Strobel et al. 1980, 1981, 1985, 1992, 1997, 2001) and is continually updated with recent high-resolution observations.
In this section, there are 68 stars in common with the PASTEL table, and the comparison results are shown in Fig. 9. The average difference in metallicity is 0.000 dex, with a scatter of 0.093 dex for corrected [Fe/H] and LASP metallicity. The standard deviation for corrected [Fe/H] and the LASP metallicity is 0.096 dex. The difference in metallicity between PASTEL and other high-resolution surveys depends on the stellar iron ratio. However, for open clusters, the average deviation for both PASTEL and other surveys is approximately 0.04 dex. According to Soubiran et al. (2022), the discrepancy and scatter are around −0.01 ±0.05 dex for APOGEE, 0.00±0.07 dex for GALAH, and 0.02±0.05 dex for Gaia-ESO. The dispersion we find is larger than what is reported in the literature. However, due to the limited number of samples that overlap, there are no noticeable patterns or biases between the corrected data and the original results.
Additionally, a total of 592 cross-match samples were chosen between this study and APOGEE DR17 (Abdurro’uf et al. 2022), and the comparison results are listed in Fig. 10. The average difference for this comparison is approximately −0.053 dex with a standard deviation of 0.140 dex for the original result, and −0.020 dex with an associated error of 0.160 dex for the corrected result. There is no systematic bias affecting our comparison of the metallicities provided by APOGEE DR17 and LAMOST LRS. However, there is an offset for stars with effective temperatures outside of a certain range (i.e., <5000 K and 6500 K >) and for dwarfs with log g > 4.5 dex in the original results. This bias has been improved to some degree after correction.
Furthermore, 453 stars were selected that are in common with the sample of Spina et al. (2021). In Fig. 11, the [Fe/H] residual to effective temperature is shown. The blue point represents the difference between the literature value and the original value from the LAMOST catalog, while the orange point represents the difference between Spina’s result and the [Fe/H] derived from LAMOST after correction. The blue point indicates a systematically lower offset compared to the GALAH result, especially for targets with temperatures cooler than 5000 K and with surface gravity of around 4.5 dex. The differences have been calibrated after the correction procedure described in Sect. 2.2, as shown by the orange dots in Fig. 11. Both systematic offsets are close to 0 dex before and after correction. However, for stars with Teff < 5000 K and log g ranging from 4.0 dex to 5.0 dex, the bias decreases in Fig. 11. Furthermore, the comparison indicates no systematic lower bias for hotter stars (above 6500 K).
Additionally, we examined the iron abundance of the final members, and stars with [Fe/H] values that fall outside 3σ of the standard deviation for OCs are excluded. Figure 12 displays the distribution of the [Fe/H] membership for six OCs (Collinder 69, Melotte 111, Melotte 20, Melotte 22, NGC 2682, and NGC 6819) that were studied in Spina et al. (2021) and selected as samples for internal comparisons within the OC. The blue, orange, and green points represent the iron abundance values from Spina et al. (2021), the corrected values, and the original values, respectively. For all selected OCs samples, the [Fe/H] dispersion is decreased after correction.
![]() |
Fig. 7 Kiel diagram displaying the member distribution for six selected OCs, color-coded according to metallicity. The x-axis corresponds to the effective temperature, while the y-axis corresponds to the surface gravity. The gray dashed track represents the PARSEC isochrone, which is based on the OC parameters from Cavallo et al. (2024). |
![]() |
Fig. 8 Kiel diagram and isochrone for Collinder 69, divided into six age bins. |
![]() |
Fig. 9 [Fe/H] residual for 68 common stars between this work and that from the PASTEL catalog. |
![]() |
Fig. 10 Residuals between APOGEE DR17 and LAMOST LRS calculated based on the effective temperature and surface gravity. The blue dots represent the difference between HRS and the original LRS result, while the orange curve shows the difference between HRS and LRS after correction. The gray-shaded region includes samples with an effective temperature of below 5000 K or above 6500 K and dwarfs with a surface gravity of around 4.5 dex. |
3.2 Comparison of OCs in terms of metallicity
We looked for differences between the metallicities derived in the present study and those derived by other authors from the HRS of the OCs. We include OC samples observed by APOGEE (Donor et al. 2020), Gaia-ESO (Casali et al. 2020), and GALAH DR3(Spina et al. 2021). Additionally, we include data from the SPA project (Frasca et al. 2019; Casali et al. 2020; D’Orazi et al. 2020; Zhang et al. 2021a), OCCASO (Casamiquela et al. 2016), OSTTA (Casamiquela et al. 2016), and the high-quality HRS results from Netopil et al. (2022).
A total of 88 samples have been selected for this comparison. Figure 13 shows the linear relationships between the metallicities derived in the present study and those from other studies based on HRS. The bottom panel of Fig. 13 displays the residual metallicity of LAMOST’s iron abundance. According to the literature, the average difference in this comparison is approximately −0.023 dex, with a dispersion of 0.114 dex. The red circles represent OCs with LAMOST membership of fewer than three stars, while symbols with green crosses represent HRS samples with fewer than three stars within a single OC. This suggests that the accuracy of average metallicity changes with the number of stars for each OC.
3.3 Metallicity of OCs as derived with Monte Carlo sampling
In this section, we use the Monte Carlo (MC) sampling method to perform parameter quality control and determine the metallicity of OCs for data comparison. The objective is to prevent any outliers in [Fe/H] within a single OC. This is because there is a slight dispersion (within 0.1 dex) in the metallicity of OC memberships, even in different evolution phases (Bertelli Motta et al. 2018; Semenova et al. 2020). The method primarily focuses on the accuracy of the stellar parameters, membership selection, evolution stage checking, and data cuts, which are discussed in the previous section. The basic parameter assumptions are based on Fu et al. (2022). The uncertainties in [Fe/H] follow a Gaussian distribution, with a mean uncertainty of 0.07 dex, which is the typical error for the LASP pipeline. For each cluster, we randomly selected 5000 samples to determine the mean [Fe/H]mean and median [Fe/H]med values, finding errors of around 0.12 dex and 0.15 dex, respectively. The mean value is used as the OC metallicity, and the standard deviation is the final uncertainty. Figure 14 shows the correlation between OC metallicity based on ANN and MC. Generally, there is good agreement between the OC metallicity obtained from the two methods, with a systematic bias of approximately −0.005 dex and a standard deviation of 0.064 dex. However, the [Fe/H] derived from ANN is used for further discussion for each OC.
![]() |
Fig. 11 Residuals for metallicities between HRS (Spina et al. 2021) and LRS with respect to effective temperature and surface gravity. The blue dots are the difference between HRS and the original result from LRS, and the orange circles show the difference between HRS and LRS after correction. The gray shadow highlights stars with Teff < 5000 K and log g > 4.5 dex. |
![]() |
Fig. 12 Membership metallicity distribution vs. Teff and log g for six example OCs analyzed in the present study. The blue dots are the stars studied by Spina et al. (2021). The green points are from LASP DR8 and the orange points are the metallicity for the present study after correction. |
![]() |
Fig. 13 Comparison of [Fe/H] Distribution between LAMOST and HRS Data. Left panel: distribution of [Fe/H] for both our results and results from the literature based on HRS. The metallicities from these two datasets follow the function [Fe/H]HRS = 0·729[Fe/H]LAMOST − 0.01. The dashed line represents a 1:1 linear relationship. Right panel: residual metallicities between the two datasets, with an average difference of −0.0023±0.114 dex. The red circles represent clusters with fewer than three members in the LAMOST data, while the green crosses represent HRS OCs with fewer than three members. |
![]() |
Fig. 14 Linear regression relationship between OC metallicities determined with an ANN model (x-axis) and an MCMC method (y-axis). |
4 Discussion
Most of the targets observed by LAMOST are located in the anti-Galactic center direction, around a galactocentric distance of 8 kpc, ranging from areas near the Sun to farther regions, but with a lack of samples between the Galactic center and the Sun. These samples are ideal for studying the chemical properties and migration patterns of the solar neighborhood. In a study by Schönrich & Binney (2009), the authors discussed the migration of stars from the inner disk to the vicinity of the Sun in the Galactic disk. Galactic disk. Their model suggests that this migration is linked to the formation of the thin and thick disks. This section of the discussion focuses on the chemical properties as they vary with galactocentric distance (RGC) and guiding distance (R𝑔uis). We will discuss a method using a semi-empirical function to estimate the birthplaces of OCs later, particularly for samples located in the solar neighborhood.
4.1 Guiding radii of LAMOST OCs
We used Gaia DR3 and the reference catalog to gather photometric information about OCs, which includes some basic details. However, instead of using the current galactocentric radius, we can obtain more accurate metallicity gradients by using the guiding center radius of a cluster, which is not affected by orbital blurring (Netopil et al. 2022; Spina et al. 2021; Zhang et al. 2021b). The guiding distance is calculated as R𝑔uidin𝑔= Lz/Vc(R) where Lz represents the angular momentum in the z-axis in the Galactic coordinate system, and Vc(R) represents the circular velocity at the cylindrical radius R. The Cavallo catalog does not provide the guiding distance of the OCs. However, in a study by Myers et al. (2022), the authors calculated the guiding distance using the Milky Way model mentioned by Price-Whelan et al. (2021). Additionally, Spina et al. (2021) discuss the distribution of [Fe/H]−R𝑔uidin𝑔 and calculate R𝑔uidin𝑔 using the GALPY package6. In this section, we present the dynamical properties of our sample of OCs as calculated using the GALPY package (Bovy 2015). The inputs for the calculation of R𝑔uidin𝑔 include observational parameters such as coordinates, distance, proper motions, and RVs. The work by Spina et al. (2021) discusses the angular momenta JR, LZ, and JZ, as well as the guiding radius R𝑔uid with respect to the MWpotential2014 Bovy (2015). The information for the current samples of OCs, including their metallicities and spatial distributions, can be found in Table A1.
4.2 Metallicity distribution across Galactic space
This section examines the metallicity properties across various regions of the Galaxy’s spatial distribution. We discuss the relationship between the metallicity distribution [Fe/H]−RGC with Galactocentric distance and the spatial distribution of metallicity abundance with cluster guiding radius [Fe/H]−R𝑔uid. The model and the observational data indicate that the distribution of metallicity shows a decreasing trend with increasing distance from the Galactic center. In this section, we discuss the characteristics of metallicity as it changes with galactocentric distance. The age and galactocentric distance of the OC samples in this section are based on the catalog provided by Cavallo et al. (2024), which contains information for over 1000 OCs younger than 1 Gyr, with the oldest OC being 9 Gyr old. These OCs cover a wide range of galactocentric distances (RGC) from 4 kpc to 16 kpc, with more than 100 OCs showing a RGC of greater than 11 kpc.
For clarity, the error associated with the OC metallicity is not plotted in Figure 15 and the subsequent figures, which shows the metallicity distribution versus RGC for all age and distance ranges compared with a similar distribution based on the literature. Despite the presence of over 1000 OCs in this study, only a limited number of these clusters are located within 8 kpc. This scarcity contributes to discrepancies between our analysis results and those reported in the literature. In this section, the Bayesian linear regression between metallicity and galactocentric distance with PYMC3 package Salvatier et al. (2016), adopted a similar estimate to Spina et al. (2021) and Cavallo et al. (2024). A simple linear regression relationship is used in the present fitting: Yi = α * Xi + β. Here, Yi and Xi are RGC and the OC metallicity following a normal distribution:
(1)
and the values are RGC and [Fe/H] for each OC, which have the maximum possibility in the normal distribution. The is the calculated standard deviation based on the RGC estimated using the ANN model (Cavallo et al. 2024), and the dispersion of [Fe/H] is derived in the metallicity calculation for each OC. The priors for α and β follow the distributions N(−0.068 dex/kpc, 0.1 dex/kpc) and N(0.5 dex, 1 dex), respectively. The prior free parameter ϵ follows the positive part of the half-Cauchy distribution, where γ is 1, based on the assumptions outlined in Spina et al. (2021). This study utilizes 10 000 samples in the simulation, employing the methods mentioned above. The final parameters for the linear regression are listed in Table 2, and we employed the same fitting method as reported in the literature. We show the gradient of [Fe/H]−RGC from other work in Table 3; the trend is around −0.060 dex kpc−1 (Jacobson et al. 2016; Donor et al. 2020; Spina et al. 2021; Zhang et al. 2021a; Myers et al. 2022). The metallicity gradient is relatively flat in the present calculation, with a value of −0.045 ± 0.006 dex kpc−1; however, the trend from these latter works is based on HRS, and the measurement we obtain in the present study is less precise. The trend increases to −0.054 ± 0.004 dex kpc−1 in the studies that use the same fitting method as that employed here. Moreover, the metallicity trend derived from Friel et al. (2002), Reddy et al. (2016), Carrera et al. (2019), Casamiquela et al. (2019), and Cavallo et al. (2024) is comparable with this measurement considering the uncertainties. We adopted the same Bayesian fitting method to fit the linear regression parameters of [Fe/H] and OC guiding distance. The parameters of the linear relationship are listed in Table. 2.
Figures 15 and 16 show the metallicity distribution as a function of either distance from galactic center or guiding distance for the present samples, respectively. In Fig. 15, the orange dots represent the OCs from the present work, while the blue points represent OCs selected from the literature where HRS were used to derive metallicity. The light shadow and dark shadow are the 95% and 68% confidence intervals for the final Bayesian regression. In this calculation, only 111 OCs contain more than three stars, which may lead to uncertainties in metallicity determination. We used the MCMC method mentioned in Spina et al. (2021) to determine that the metallicity gradient flattens as the Galactocentric distance reaches 9.129 kpc. The metallicity gradient in the inner disk is −0.041 ± 0.078 dex kpc−1. However, the processes of churning and blurring can lead to significant dispersion in the samples, as noted by Spina et al. (2021). In contrast, the outer disk shows a metallicity trend of −0.029 ± 0.013 dex kpc−1. Donor et al. (2020) mention a “Keen” at around 12 kpc, but the limited number of OCs within 7 kpc may account for the observed difference. All of the samples mentioned here are situated within 1 kpc from the Galactic midplane. The metallicity trend, which can be observed in Fig. 16, follows the gradient [Fe/H]−Rցuide, corresponding to RGC. The value of this trend is −0.037 ± 0.003 dex kpc−1, suggesting that the gradient is slightly flatter than the [Fe/H]−RGC relationship, but the two generally align. We observed a change in the trend of [Fe/H]−Rցuide at a guiding distance of approximately 9.518 kpc. At this distance, the gradient will shift from −0.062 ± 0.041 to −0.025 ± 0.011. To evaluate the alignment between our fitted results and the observed data, we conducted a Kolmogorov-Smirnov (KS) test, considering scenarios with and without metallicity uncertainties. We note that the maximum discrepancy between the cumulative distribution functions (CDFs) of the fitted and observed data is approximately 29%. Additionally, the inclusion or exclusion of [Fe/H]-related uncertainties appears to have little impact on the results. We conducted a similar comparison in Spina et al. (2021) for the posterior distributions of fitting parameters shown in Fig. 17. However, as opposed to these latter authors, who reported an identical posterior distribution for α and β, we find a clear difference in the posterior distribution.
Linear regression parameters and the 95% confidence intervals. The left panel shows the fitting results for the present sample and the left panel shows the linear regression parameters for selected samples from previous works.
![]() |
Fig. 15 Metallicity distribution as a function of galactocentric distance for OCs of all ages. The orange dots represent the sample from this work, while the blue points represent OCs selected from the literature where HRS was used to determine metallicity. The linear regression models for the present samples are shown in dark orange and orange, representing the 68% and 95% confidence intervals, respectively. Similarly, the confidence intervals for the literature samples are shown in dark blue and blue, representing the 68% and 95% confidence intervals, respectively. |
![]() |
Fig. 16 Metallicity distribution as a function of guiding distance for our sample of OCs. The dark and light orange areas are the 68% and 95% confidence intervals for our sample, respectively. |
Slope of the metallicity gradient from selected literature papers.
4.3 Metallicity distribution compared with a chemo-dynamical model
The chemical properties we presently observe are the outcome of intricate evolutionary processes, disk asymmetries, and non-equilibrium processes, as mentioned in studies of the Gaia dataset, such as Wang et al. (2023c) and Frankel et al. (2018). For instance, according to Sellwood & Binney (2002), migration within the disk is a crucial process that needs to be taken into account when modeling galactic evolution. In their work, Minchev et al. (2013, 2014a,b) considered the impact of migration and kinematic heating on galactic disk modeling, but did not address the issues related to self-consistent simulation.
In this section, we discuss the predictions based on chemodynamical models (Minchev et al. 2014a, e.g. the MCM model) and compare these to observational findings. Minchev et al. (2014a) provide the galactic metallicity distribution based on the galactocentric distance split into 17 age bins. MCM takes into account the distance from the disk midplane and migration, combining chemical evolution and dynamics.
Figures B1 and B2 show observational results compared with predictions split by age and color-coded according to the number of OCs in each cluster. The chemo-dynamic track is separated into |z| < 0.3 kpc and 0.3 kpc < |z| < 0.8 kpc (Minchev et al. 2014a). In the comparison plots, the red line represents the predicted track with |z| < 0.3 kpc, while the blue line represents the predicted track for 0.3 kpc < |z| < 0.8 kpc. According to the model, the midplane distance of OCs for these samples is limited to 0.8 kpc. The cross symbols represent OCs with fewer than three members. As the plot shows, most of the OCs of the compared sample are younger than 0.6 Gyr. The metallicity distribution is systematically lower for the youngest age bins (0.0 < age < 0.3 Gyr and 0.3 < age < 0.6 Gyr), and OCs with fewer than three members show a greater dispersion compared to the chemo-dynamic prediction. The observational data are consistent with the predictions for ages of greater than 1.2 Gyr, as mentioned in previous work (Zhang et al. 2021a).
![]() |
Fig. 17 Posterior distributions of parameters α and β obtained using Bayesian linear regression analyses of [Fe/H]−RGC and [Fe/H]−Rցuide. |
4.4 Birth radii of open clusters
Sellwood & Binney (2002) first suggested that transient spiral arms can permanently change the angular momentum of stars without heating their orbits, causing them to migrate away from their birth location over time. Moreover, due to the exponential density profile of the Galaxy, a population of field stars will preferably move from the inner disk to the outer disk. Reddy et al. (2016) suggest that OCs follow a similar trend to field stars, in that those stars currently in the outer disk were born within 12 kpc and have migrated. This process can flatten the birth metallicity gradient over time (Minchev et al. 2013); as a result, radial migration has to be taken into account when trying to understand the birth environment of OCs (Minchev et al. 2018; Ratcliffe et al. 2023). The present work adopted an empirical method from Lu et al. (2022a) to infer the birth radii of their OC sample. This method is built upon the previous method developed by Minchev et al. (2018) and the assumption that the ISM metallicity gradient (which is a function of lookback time) has a linear relationship with the metallicity range in mono-age populations.
For the present study, we used age and metallicity measurements to estimate the birth radii for the OCs of our sample following the method described in Lu et al. (2022a). The birth radii, Rb(age, [Fe/H]), are calculated using
(2)
where [Fe/H] is the observed metallicity, and [Fe/H] (0, τ) and ∆[Fe/H] (τ) represent the metallicity at the Galactic center and the metallicity gradient at a given lookback time, τ, respectively. ∆[Fe/H] (τ) and [Fe/H] (0,τ) are taken from Table A1 in Lu et al. (2022a). Figure 18 shows the migration distance for the OCs in this study as a function of age. Surprisingly, according to this analysis, most OCs have migrated from the outer disk to their current guiding radius. However, the oldest OCs have migrated from the inner Galaxy. This trend, if real, could suggest that the youngest clusters most likely cannot survive the tidal disruption of the inner Galaxy. However, as the relation used to infer birth radii is only calibrated with field stars of > 1.5 Gyr old (Xiang & Rix 2022), inferring birth radii for these young clusters can be biased. As a result, we do not draw any conclusions on this topic here, and propose that further testing is necessary, with particular attention to any possible bias that might affect the derived birth radii of young field stars.
5 Summary
This study uses low-resolution spectra obtained from the LAMOST to re-evaluate the metallicity abundance of over 1131 open clusters (OCs) identified in the Gaia DR3 catalog (Cavallo et al. 2024). The selected OCs in this study span an age range from a few million years to 9.5 Gyr, with the majority of OCs being younger than 600 Myr. These samples are primarily located in the solar neighborhood, covering a galactocentric distance ranging from 4.68 kpc to 16.7 kpc. The key findings of this study are summarized as follows:
We used the LASP pipeline in the LAMOST LRS DR8 to determine the average RV for a set of selected OCs. The mean difference between the RVs obtained from the LASP measurement and those derived from the HRS measurement is 5.598 ± 14.031 km/s. Additionally, there is an offset of 6.075 ± 11.496 km/s between our results and the determination from Gaia DR3;
The stars with a temperature of below 5000 K and the dwarfs with a surface gravity of around 4.5 dex consistently exhibit lower [Fe/H] values compared to the original LAMOST DR8 catalog. To determine the metallicity of the OCs, we used an ANN that incorporates the correction method introduced by Niu et al. (2023). This approach helped us re-identify the metallicity of the stars and minimize bias in the process. Consequently, the discrepancy in metallicity between our work and the HRS result is compatible after corrections are made;
We examined the distribution of metallicity across the spatial extent of the galaxy. The metallicity gradient, determined through Bayesian linear regression, is calculated as −0.0045 * RGC + 0.304 km dex−1. This gradient is relatively flat. Notably, this trend is consistent with the results obtained from literature samples and the previous calculation presented in Table 3;
We compared the observed trend of metallicity with the chemo-dynamical predictions (Minchev et al. 2013, 2014a,b) based on age. The results from LAMOST indicate a lower metallicity compared to the predictions when the observational uncertainties are taken into account. Moreover, the metallicities of observed OCs are not taken into account. Furthermore, the metallicities of observed OCs in the youngest age groups show a larger dispersion, which can be attributed to the fact that most OCs currently have a limited number of members – typically around three stars. This large scatter is also present in the youngest age groups of the HRS (Zhang et al. 2021a), and is likely due to the measurement uncertainties associated with young stars.
We determined the birth radii of our sample of OCs using the method presented by Lu et al. (2022a). Our findings suggest that the majority of the OCs near the Sun originate from the outer disk, indicating that the youngest OCs may have been destroyed in the inner disk as a result of tidal interaction. However, additional testing is required to confirm these findings.
![]() |
Fig. 18 Migration distance for the present sample of OCs as a function of their age. |
Data availability
Full Tables 1 and A.1 are available at the CDS via anonymous ftp to cdsarc.cds.unistra.fr (130.79.128.5) or via https://cdsarc.cds.unistra.fr/viz-bin/cat/J/A+A/692/A212
Acknowledgements
We would like to thank the anonymous referee for his/her very helpful and insightful comments and Ivan Minchev for sharing his evolutionary models. GJW acknowledges the Africa Europe Cluster of Research Excellence (CoRE-AI) fellowship. HFW is supported in this work by the Department of Physics and Astronomy of Padova University though the 2022 ARPE grant: Rediscovering our Galaxy with machines. JXC acknowledges the Support of the Postdoctoral Fellowship Program of CPSF under Grant Number GZC20240124. This work exploits the Simbad, Vizier, and NASA-ADS databases and the software TOPCAT (Taylor 2005). The Guo Shou Jing 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 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). This work is supported by the National Natural Science Foundation of China under program Nos. 12090040, 12090043, and 12003025, as well as the Basic Research Program of Yunnan Province (No. 202401AT070142). We also express our gratitude for the generous support from the International Centre of Supernovae, Yunnan Key Laboratory (No. 202302AN360001), and the Natural Science Foundation of Yunnan Province (No. 202201BC070003).
Appendix A The physics information for Individual OCs
Basic information of OCs in the present samples (extract). The full table is available at the CDS.
Appendix B The comparison between and chemo-dynamic model
![]() |
Fig. B1 Comparison between chemo-dynamic model (Minchev et al. 2014a,b) and observation for clusters younger than 3Gyr. The red lines and the blue lines are predictions from the MCM models for |z| < 0.3 kpc and 0.3 < |z| <0.8 kpc respectively. The colors in the cycle are the number of stars for each OC whose membership is larger than 3, and the cross symbol is the OC with a member less than 3. |
![]() |
Fig. B2 Same as Fig. B1 but for clusters between 1.5 and 4.5 Gyr. Note the paucity of OCs older than 2.5 Gyr |
References
- Abdurro’uf, Accetta, K., Aerts, C., et al. 2022, ApJS, 259, 35 [NASA ADS] [CrossRef] [Google Scholar]
- Anders, F., Gispert, P., Ratcliffe, B., et al. 2023, A&A, 678, A158 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Andreuzzi, G., Bragaglia, A., Tosi, M., et al. 2011, MNRAS, 412, 1265 [NASA ADS] [Google Scholar]
- Bai, Z.-R., Zhang, H.-T., Yuan, H.-L., et al. 2021, Res. Astron. Astrophys., 21, 249 [CrossRef] [Google Scholar]
- Bertelli Motta, C., Pasquali, A., Richer, J., et al. 2018, MNRAS, 478, 425 [NASA ADS] [CrossRef] [Google Scholar]
- Bird, J. C., Kazantzidis, S., & Weinberg, D. H. 2012, MNRAS, 420, 913 [NASA ADS] [CrossRef] [Google Scholar]
- Bovy, J. 2015, ApJS, 216, 29 [NASA ADS] [CrossRef] [Google Scholar]
- Buder, S., Sharma, S., Kos, J., et al. 2021, MNRAS, 506, 150 [NASA ADS] [CrossRef] [Google Scholar]
- Casamiquela, L., Carrera, R., Jordi, C., et al. 2016, MNRAS, 458, 3150 [NASA ADS] [CrossRef] [Google Scholar]
- Casamiquela, L., Blanco-Cuaresma, S., Carrera, R., et al. 2019, MNRAS, 490, 1821 [NASA ADS] [CrossRef] [Google Scholar]
- Casamiquela, L., Castro-Ginard, A., Anders, F., et al. 2021, A&A, 654, A151 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Casali, G., Magrini, L., Tognelli, E., et al. 2019, A&A, 629, A62 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Casali, G., Magrini, L., Frasca, A., et al. 2020, A&A, 643, A12 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Cantat-Gaudin, T., Anders, F., Castro-Ginard, A., et al. 2020, A&A, 640, A1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Cavallo, L., Spina, L., Carraro, G., et al. 2024, AJ, 167, 12 [NASA ADS] [CrossRef] [Google Scholar]
- Carrera, R., Bragaglia, A., Cantat-Gaudin, T., et al. 2019, A&A, 623, A80 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Carrera, R., Casamiquela, L., Bragaglia, A., et al. 2022, A&A, 663, A148 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Cayrel de Strobel, G., Bentolila, C., Hauck, B., et al. 1980, A&AS, 41, 405 [NASA ADS] [Google Scholar]
- Cayrel de Strobel, G., Bentolila, C., Hauck, B., et al. 1981, A&AS, 45, 97 [NASA ADS] [Google Scholar]
- Cayrel de Strobel, G., Bentolila, C., Hauck, B., et al. 1985, A&AS, 59, 145 [NASA ADS] [Google Scholar]
- Cayrel de Strobel, G., Hauck, B., Francois, P., et al. 1992, A&AS, 95, 273 [NASA ADS] [Google Scholar]
- Cayrel de Strobel, G., Soubiran, C., Friel, E. D., et al. 1997, A&AS, 124, 299 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Cayrel de Strobel, G., Soubiran, C., & Ralite, N. 2001, A&A, 373, 159 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Chaplin, W. J., & Miglio, A. 2013, ARA&A, 51, 353 [Google Scholar]
- Chaplin, W. J., Basu, S., Huber, D., et al. 2014, ApJS, 210, 1 [Google Scholar]
- Chaplin, W. J., Serenelli, A. M., Miglio, A., et al. 2020, Nat. Astron., 4, 382 [Google Scholar]
- Chen, Y. Q., & Zhao, G. 2020, MNRAS, 495, 2673 [NASA ADS] [CrossRef] [Google Scholar]
- Cropper, M., Katz, D., Sartoretti, P., et al. 2018, A&A, 616, A5 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Cui, X.-Q., Zhao, Y.-H., Chu, Y.-Q., et al. 2012, Res. Astron. Astrophys., 12, 1197 [Google Scholar]
- Deng, L.-C., Newberg, H. J., Liu, C., et al. 2012, Res. Astron. Astrophys., 12, 735 [Google Scholar]
- De Silva, G. M., Sneden, C., Paulson, D. B., et al. 2006, AJ, 131, 455 [NASA ADS] [CrossRef] [Google Scholar]
- De Silva, G. M., Freeman, K. C., Bland-Hawthorn, J., et al. 2015, MNRAS, 449, 2604 [NASA ADS] [CrossRef] [Google Scholar]
- Donor, J., Frinchaboy, P. M., Cunha, K., et al. 2020, AJ, 159, 199 [NASA ADS] [CrossRef] [Google Scholar]
- Du, B., Luo, A., Zhang, J., et al. 2012, Proc. SPIE, 8451, 845137 [NASA ADS] [CrossRef] [Google Scholar]
- D’Orazi, V., Oliva, E., Bragaglia, A., et al. 2020, A&A, 633, A38 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Frankel, N., Rix, H.-W., Ting, Y.-S., et al. 2018, ApJ, 865, 96 [NASA ADS] [CrossRef] [Google Scholar]
- Frasca, A., Alonso-Santiago, J., Catanzaro, G., et al. 2019, A&A, 632, A16 [NASA ADS] [EDP Sciences] [Google Scholar]
- Freeman, K., & Bland-Hawthorn, J. 2002, ARA&A, 40, 487 [Google Scholar]
- Friel, E. D. 1995, ARA&A, 33, 381 [Google Scholar]
- Friel, E. D., Janes, K. A., Tavarez, M., et al. 2002, AJ, 124, 2693 [NASA ADS] [CrossRef] [Google Scholar]
- Friel, E. D., Jacobson, H. R., & Pilachowski, C. A. 2010, AJ, 139, 1942 [NASA ADS] [CrossRef] [Google Scholar]
- Fu, X., Bragaglia, A., Liu, C., et al. 2022, A&A, 668, A4 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Fujii, M. S., & Portegies Zwart, S. 2016, ApJ, 817, 4 [NASA ADS] [CrossRef] [Google Scholar]
- Gaia Collaboration (Vallenari, A., et al.) 2023, A&A, 674, A1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Genovali, K., Lemasle, B., Bono, G., et al. 2014, A&A, 566, A37 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Grenon, M. 1972, IAU Colloq. 17: Age des Etoiles, 55 [Google Scholar]
- Grenon, M. 1989, Ap&SS, 156, 29 [NASA ADS] [CrossRef] [Google Scholar]
- Grand, R. J. J., Kawata, D., & Cropper, M. 2015, MNRAS, 447, 4018 [NASA ADS] [CrossRef] [Google Scholar]
- Gilmore, G., Randich, S., Worley, C. C., et al. 2022, A&A, 666, A120 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Hayden, M. R., Sharma, S., Bland-Hawthorn, J., et al. 2022, MNRAS, 517, 5325 [NASA ADS] [CrossRef] [Google Scholar]
- Hunt, E. L., & Reffert, S. 2021, A&A, 646, A104 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Hunt, E. L., & Reffert, S. 2023, A&A, 673, A114 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Jacobson, H. R., Friel, E. D., Jílková, L., et al. 2016, A&A, 591, A37 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Jofré, P., Heiter, U., & Soubiran, C. 2019, ARA&A, 57, 571 [Google Scholar]
- Katz, D., Sartoretti, P., Guerrier, A., et al. 2023, A&A, 674, A5 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Koleva, M., Prugniel, P., Bouchard, A., et al. 2009, A&A, 501, 1269 [CrossRef] [EDP Sciences] [Google Scholar]
- Kordopatis, G., Schultheis, M., McMillan, P. J., et al. 2023, A&A, 669, A104 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Kharchenko, N. V., Piskunov, A. E., Schilbach, E., et al. 2013, A&A, 558, A53 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Kubryk, M., Prantzos, N., & Athanassoula, E. 2013, MNRAS, 436, 1479 [NASA ADS] [CrossRef] [Google Scholar]
- Li, X., & Lin, B. 2023, MNRAS, 521, 6354 [NASA ADS] [CrossRef] [Google Scholar]
- Lépine, J. R. D., Cruz, P., Scarano, S., et al. 2011, MNRAS, 417, 698 [CrossRef] [Google Scholar]
- Liu, C., Fu, J., Shi, J., et al. 2020, arXiv e-prints [arXiv:2005.07210] [Google Scholar]
- Luo, A.-L., Zhao, Y.-H., Zhao, G., et al. 2015, RA&A, 15, 1095 [Google Scholar]
- Lu, Y., Minchev, I., Buck, T., et al. 2022, MNRAS, 535, 392 [Google Scholar]
- Mackereth, J. T., Bovy, J., Leung, H. W., et al. 2019, MNRAS, 489, 176 [Google Scholar]
- Martig, M., Minchev, I., & Flynn, C. 2014, MNRAS, 443, 2452 [Google Scholar]
- Marton, G., Ábrahám, P., Rimoldini, L., et al. 2023, A&A, 674, A21 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Minchev, I., & Famaey, B. 2010, ApJ, 722, 112 [Google Scholar]
- Minchev, I., Famaey, B., Combes, F., et al. 2011, A&A, 527, A147 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Minchev, I., Chiappini, C., & Martig, M. 2013, A&A, 558, A9 [CrossRef] [EDP Sciences] [Google Scholar]
- Minchev, I., Chiappini, C., Martig, M., et al. 2014a, ApJ, 781, L20 [Google Scholar]
- Minchev, I., Chiappini, C., & Martig, M. 2014b, A&A, 572, A92 [CrossRef] [EDP Sciences] [Google Scholar]
- Minchev, I., Anders, F., Recio-Blanco, A., et al. 2018, MNRAS, 481, 1645 [NASA ADS] [CrossRef] [Google Scholar]
- Moya, A., Sarro, L. M., Delgado-Mena, E., et al. 2022, A&A, 660, A15 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Myers, N., Donor, J., Spoo, T., et al. 2022, AJ, 164, 85 [NASA ADS] [CrossRef] [Google Scholar]
- Nair, V., & Hinton, G. E. 2010, Proc. of the 27th Int. Conf. on Machine Learning, eds. J. Fürnkranz, & T. Joachims (Madison, WI: Omnipress), 807 [Google Scholar]
- Ness, M. K., Wheeler, A. J., McKinnon, K., et al. 2022, ApJ, 926, 144 [NASA ADS] [CrossRef] [Google Scholar]
- Nguyen, C. T., Costa, G., Girardi, L., et al. 2022, A&A, 665, A126 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Niu, Z., Yuan, H., & Liu, J. 2023, ApJ, 950, 104 [NASA ADS] [CrossRef] [Google Scholar]
- Netopil, M., Paunzen, E., Heiter, U., et al. 2016, A&A, 585, A150 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Netopil, M., Oralhan, I. A., Çakmak, H., et al. 2022, MNRAS, 509, 421 [Google Scholar]
- Price-Whelan, A. M., Hogg, D. W., Johnston, K. V., et al. 2021, ApJ, 910, 17 [NASA ADS] [CrossRef] [Google Scholar]
- Queiroz, A. B. A., Anders, F., Santiago, B. X., et al. 2018, MNRAS, 476, 2556 [Google Scholar]
- Ratcliffe, B., Minchev, I., Anders, F., et al. 2023, MNRAS, 525, 2208 [NASA ADS] [CrossRef] [Google Scholar]
- Ren, J.-J., Wu, H., Wu, C.-J., et al. 2021, RA&A, 21, 051 [Google Scholar]
- Recio-Blanco, A., de Laverny, P., Palicio, P. A., et al. 2023, A&A, 674, A29 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Reddy, A. B. S., Lambert, D. L., & Giridhar, S. 2016, MNRAS, 463, 4366 [NASA ADS] [CrossRef] [Google Scholar]
- Riello, M., De Angeli, F., Evans, D. W., et al. 2021, A&A, 649, A3 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Roškar, R., Debattista, V. P., Stinson, G. S., et al. 2008, ApJ, 675, L65 [Google Scholar]
- Salaris, M., Weiss, A., & Percival, S. M. 2004, A&A, 414, 163 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Sanders, J. L., & Das, P. 2018, MNRAS, 481, 4093 [CrossRef] [Google Scholar]
- Salvatier, J., Wieckiâ, T. V., & Fonnesbeck, C. 2016, Astrophysics Source Code Library [record ascl:1610.016] [Google Scholar]
- Sellwood, J. A., & Binney, J. J. 2002, MNRAS, 336, 785 [Google Scholar]
- Schönrich, R., & Binney, J. 2009, MNRAS, 399, 1145 [Google Scholar]
- Schiappacasse-Ulloa, J., Pasquato, M., & Lucatello, S. 2023, arXiv e-prints [arXiv:2306.15487] [Google Scholar]
- Semenova, E., Bergemann, M., Deal, M., et al. 2020, A&A, 643, A164 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Sharma, S., Hayden, M. R., Bland-Hawthorn, J., et al. 2021, MNRAS, 506, 1761 [NASA ADS] [CrossRef] [Google Scholar]
- Soubiran, C., Le Campion, J.-F., Brouillet, N., et al. 2016, A&A, 591, A118 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Soubiran, C., Brouillet, N., & Casamiquela, L. 2022, A&A, 663, A4 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Spina, L., Ting, Y.-S., De Silva, G. M., et al. 2021, MNRAS, 503, 3279 [CrossRef] [Google Scholar]
- Sun, T., Ge, Z., Chen, X., et al. 2023, ApJS, 268, 29 [NASA ADS] [CrossRef] [Google Scholar]
- Taylor, M. B. 2005, Astron. Data Anal. Softw. Syst. XIV, 347, 29 [Google Scholar]
- Ting, Y.-S., & Rix, H.-W. 2019, ApJ, 878, 21 [NASA ADS] [CrossRef] [Google Scholar]
- Ting, Y.-S., & Weinberg, D. H. 2022, ApJ, 927, 209 [NASA ADS] [CrossRef] [Google Scholar]
- Tsantaki, M., Pancino, E., Marrese, P., et al. 2022, A&A, 659, A95 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- van den Bergh, S., & McClure, R. D. 1980, A&A, 88, 360 [NASA ADS] [Google Scholar]
- Viscasillas Vázquez, C., Magrini, L., Casali, G., et al. 2022, A&A, 660, A135 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Viscasillas Vázquez, C., Magrini, L., Spina, L., et al. 2023, A&A, 679, A122 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Wang, R., Luo, A.-L., Chen, J.-J., et al. 2019, ApJS, 244, 27 [NASA ADS] [CrossRef] [Google Scholar]
- Wang, G.-J., Li, S.-Y., & Xia, J.-Q. 2020a, ApJS, 249, 25 [NASA ADS] [CrossRef] [Google Scholar]
- Wang, H.-F., López-Corredoira, M., Huang, Y., et al. 2020b, MNRAS, 491, 2104 [NASA ADS] [CrossRef] [Google Scholar]
- Wang, H.-F., Huang, Y., Zhang, H.-W., et al. 2020c, ApJ, 902, 70 [NASA ADS] [CrossRef] [Google Scholar]
- Wang, G.-J., Cheng, C., Ma, Y.-Z., & Xia, J.-Q. 2022, ApJS, 262, 24 [NASA ADS] [CrossRef] [Google Scholar]
- Wang, R., Luo, A.-L., Zhang, S., et al. 2023a, ApJS, 266, 40 [NASA ADS] [CrossRef] [Google Scholar]
- Wang, G.-J., Cheng, C., Ma, Y.-Z., et al. 2023b, ApJS, 268, 7 [NASA ADS] [CrossRef] [Google Scholar]
- Wang, H.-F., Chrobáková, Ž., López-Corredoira, M., et al. 2023c, ApJ, 942, 12 [NASA ADS] [CrossRef] [Google Scholar]
- Wang, H.-F., Carraro, G., Li, X., et al. 2024, ApJ, in press [arXiv:2310.17196] [Google Scholar]
- Wu, Y., Singh, H. P., Prugniel, P., et al. 2011, A&A, 525, A71 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Wu, Y., Luo, A.-L., Li, H.-N., et al. 2011, RA&A, 11, 924 [Google Scholar]
- Wu, Y., Du, B., Luo, A., et al. 2014, Statistical Challenges in 21st Century Cosmology, 306, 340 [NASA ADS] [Google Scholar]
- Xiang, M., & Rix, H.-W. 2022, Nature, 603, 599 [NASA ADS] [CrossRef] [Google Scholar]
- Xiang, M., Ting, Y.-S., Rix, H.-W., et al. 2019, ApJS, 245, 34 [Google Scholar]
- Zhang, R., Lucatello, S., Bragaglia, A., et al. 2021a, A&A, 654, A77 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Zhang, H., Chen, Y., & Zhao, G. 2021b, ApJ, 919, 52 [NASA ADS] [CrossRef] [Google Scholar]
https://www.galah-survey.org/; GALAH has released its third set of data, which provides high-resolution (R=28 000) spectra and covers four wavelength ranges: 4713–4903, 5648–5873, 6478–6737, and 7585–7887 Å.
All Tables
Radial velocities for the OCs studied here compared to those from the literature.
Linear regression parameters and the 95% confidence intervals. The left panel shows the fitting results for the present sample and the left panel shows the linear regression parameters for selected samples from previous works.
Basic information of OCs in the present samples (extract). The full table is available at the CDS.
All Figures
![]() |
Fig. 1 Difference in RVs between 101 cross-matched targets from previous works, including those by Casali et al. (2019), Casamiquela et al. (2016), Spina et al. (2021), Frasca et al. (2019), Casali et al. (2020), D’Orazi et al. (2020), Zhang et al. (2021a), and Carrera et al. (2022), and the current study is shown. The gray lines indicate the mean offset value of 5.598 km/s with a standard deviation of 14.013 km/s for this comparison. The x-axis represents the RV of the present study, incorporating the error for the overlapping sample. |
In the text |
![]() |
Fig. 2 Difference distribution in RVs for 595 cross-matching targets between Gaia DR3 and the present study. The mean value of the offset is 6.075 with a standard deviation of 11.496 km s−1 as indicated by the gray lines in the plot. The x-axis represents the LAMOST RV and the error of the overlapped sample. The y-axis includes the RV uncertainties from Gaia DR3. |
In the text |
![]() |
Fig. 3 Structure of the ANN used in this work. The number of neurons in the hidden layer is 128. |
In the text |
![]() |
Fig. 4 Loss function of the training and validation sets. |
In the text |
![]() |
Fig. 5 Comparison of Metallicity Measurements from LAMOST DR8 and GALAH. The right panel shows the iron relation for 30 000 targets. The y-axis represents the metallicity from LAMOST dr8, while the x-axis represents the corresponding value from GALAH. The left plot displays the metallicity residuals of LAMOST and GALAH, both before and after correction. The blue dots represent the original data selected from LAMOST DR8, while the orange dot represents the LAMOST result after correction. |
In the text |
![]() |
Fig. 6 Histogram of the training set and the metallicity that needed to be corrected for the LASP result. The orange section represents the training data that were used to build the ANN model, while the blue section represents the metallicity that needed to be corrected for the selected samples. |
In the text |
![]() |
Fig. 7 Kiel diagram displaying the member distribution for six selected OCs, color-coded according to metallicity. The x-axis corresponds to the effective temperature, while the y-axis corresponds to the surface gravity. The gray dashed track represents the PARSEC isochrone, which is based on the OC parameters from Cavallo et al. (2024). |
In the text |
![]() |
Fig. 8 Kiel diagram and isochrone for Collinder 69, divided into six age bins. |
In the text |
![]() |
Fig. 9 [Fe/H] residual for 68 common stars between this work and that from the PASTEL catalog. |
In the text |
![]() |
Fig. 10 Residuals between APOGEE DR17 and LAMOST LRS calculated based on the effective temperature and surface gravity. The blue dots represent the difference between HRS and the original LRS result, while the orange curve shows the difference between HRS and LRS after correction. The gray-shaded region includes samples with an effective temperature of below 5000 K or above 6500 K and dwarfs with a surface gravity of around 4.5 dex. |
In the text |
![]() |
Fig. 11 Residuals for metallicities between HRS (Spina et al. 2021) and LRS with respect to effective temperature and surface gravity. The blue dots are the difference between HRS and the original result from LRS, and the orange circles show the difference between HRS and LRS after correction. The gray shadow highlights stars with Teff < 5000 K and log g > 4.5 dex. |
In the text |
![]() |
Fig. 12 Membership metallicity distribution vs. Teff and log g for six example OCs analyzed in the present study. The blue dots are the stars studied by Spina et al. (2021). The green points are from LASP DR8 and the orange points are the metallicity for the present study after correction. |
In the text |
![]() |
Fig. 13 Comparison of [Fe/H] Distribution between LAMOST and HRS Data. Left panel: distribution of [Fe/H] for both our results and results from the literature based on HRS. The metallicities from these two datasets follow the function [Fe/H]HRS = 0·729[Fe/H]LAMOST − 0.01. The dashed line represents a 1:1 linear relationship. Right panel: residual metallicities between the two datasets, with an average difference of −0.0023±0.114 dex. The red circles represent clusters with fewer than three members in the LAMOST data, while the green crosses represent HRS OCs with fewer than three members. |
In the text |
![]() |
Fig. 14 Linear regression relationship between OC metallicities determined with an ANN model (x-axis) and an MCMC method (y-axis). |
In the text |
![]() |
Fig. 15 Metallicity distribution as a function of galactocentric distance for OCs of all ages. The orange dots represent the sample from this work, while the blue points represent OCs selected from the literature where HRS was used to determine metallicity. The linear regression models for the present samples are shown in dark orange and orange, representing the 68% and 95% confidence intervals, respectively. Similarly, the confidence intervals for the literature samples are shown in dark blue and blue, representing the 68% and 95% confidence intervals, respectively. |
In the text |
![]() |
Fig. 16 Metallicity distribution as a function of guiding distance for our sample of OCs. The dark and light orange areas are the 68% and 95% confidence intervals for our sample, respectively. |
In the text |
![]() |
Fig. 17 Posterior distributions of parameters α and β obtained using Bayesian linear regression analyses of [Fe/H]−RGC and [Fe/H]−Rցuide. |
In the text |
![]() |
Fig. 18 Migration distance for the present sample of OCs as a function of their age. |
In the text |
![]() |
Fig. B1 Comparison between chemo-dynamic model (Minchev et al. 2014a,b) and observation for clusters younger than 3Gyr. The red lines and the blue lines are predictions from the MCM models for |z| < 0.3 kpc and 0.3 < |z| <0.8 kpc respectively. The colors in the cycle are the number of stars for each OC whose membership is larger than 3, and the cross symbol is the OC with a member less than 3. |
In the text |
![]() |
Fig. B2 Same as Fig. B1 but for clusters between 1.5 and 4.5 Gyr. Note the paucity of OCs older than 2.5 Gyr |
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.