Open Access
Issue
A&A
Volume 668, December 2022
Article Number A13
Number of page(s) 13
Section Catalogs and data
DOI https://doi.org/10.1051/0004-6361/202244570
Published online 30 November 2022

© C. J. Hao 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

Classical Cepheids are intermediate-mass stars that can exhibit regular pulsations in fundamental and overtone modes, and possess ages that overlap with those of open clusters (OCs; see Table A.1). Therefore, classical Cepheids are expected to be in OCs, and such pairs are significant objects in astronomy.

Cepheids have a famous period-luminosity relationship (PLR), also known as the Leavitt law, securing them as standard candles for constituting the cosmic distance scale and as essential indicators, empowering firmer constraints on Hubble’s law (e.g., Leavitt & Pickering 1912; Feast 1999; Madore & Freedman 1991). However, one of the critical debates in cosmology is the Hubble constant (H0) tension. Based on the current standard ∧ cold dark matter (∧CDM) model, estimates of H0 based on measurements of fluctuations in the temperature and polarization of the cosmic microwave background (CMB) from Planck and ACT+WMAP consistently yield values of 67.4 ± 0.5 and 67.6 ± 1.1 km s−1 Mpc−1, respectively (Planck Collaboration VI 2020; Aiola et al. 2020). In contrast, Riess et al. (2019, 2021) argued that the Cepheid distance scale implies a larger value of H0 = 73.2 ± 1.3 km s−1 Mpc−1, in tension with the CMB result. Using the tip of the red giant branch in the Large Magellanic Cloud, the Carnegie-Chicago Hubble Program put forward an H0 = 69.6 ± 1.9 km s−1 Mpc−1, close to the CMB prediction considering the uncertainty (Freedman et al. 2019, 2020). Likewise, Majaess (2020) pointed out that the accuracy of Cepheid distances may be swayed by neglected or inaccurate blending corrections for remote targets, which would result in an overestimated H0. Therefore, it is not clear so far whether the H0 tension means new physics are needed for the ACDM model or that unknown systematic errors exist in the above measurements. Open cluster–Cepheid pairs can help to solve this problem. As a large number of cluster members makes it possible to derive more accurate astronomical parameters than those of individual sources, OCs harboring Cepheids can optimize the zero point of the PLR for Cepheids, and ultimately contribute to reducing the H0 uncertainty by providing a Galactic calibration of the Leavitt law (e.g., Turner & Burke 2002; Chen et al. 2017; Breuval et al. 2020). Recently, using only 17 reported cluster Cepheids with Hubble Space Telescope photometry and Gaia Early Data Release 3 (EDR3) parallax, Riess et al. (2022) constrained the PLR of Cepheids and obtained a 5–7% reduction in the H0 uncertainty.

Open cluster–Cepheid pairs are important in additional ways. First, cluster Cepheids can be used to benchmark Gaia parallax, especially for distant Cepheids. Gaia parallax is affected by a systematic error, the so-called “parallax zero point”, which was first clearly identified in Data Release 2 using quasars (Gaia Collaboration 2018; Lindegren et al. 2018). Also, OC–Cepheid pairs could be an important diagnostic tool for examining the Gaia parallax offset, whereby the cluster distance is tied to more stars of perhaps a different mean color than the Cepheid, and offsets can be examined as a function of direction, magnitude, color, and distance. For an ongoing investigation using cluster Cepheids to assess the systematic errors in the Gaia data, we recommend Riess et al. (2022). Second, cluster Cepheids are expected to benchmark the stellar evolution models and theoretical isochrones (e.g., Bono et al. 2005; Turner et al. 2006). Third, OC–Cepheid pairs can be of service in elucidating the pulsating nature and dynamical evolution of Cepheids (e.g., Fry & Carney 1997; Lemasle et al. 2017; Dinnbier et al. 2022). In addition, Cepheids adhere to a period–age relation (PAR), whereby longer period Cepheids are more luminous and massive (e.g., Efremov 2003; Bono et al. 2005; Turner et al. 2012a). Therefore, the ages of OC–Cepheid pairs based on isochrone fitting to the cluster colour–magnitude diagrams (CMDs) can be used to calibrate the PAR of Cepheids (e.g., Efremov 2003; Anderson et al. 2016; De Somma et al. 2020). As both OCs and Cepheids are young objects, they can also be used as good tracers for studying spiral arms and the detailed disk structure of the Milky Way (e.g., Chen et al. 2019; Hao et al. 2021; Poggio et al. 2021). In particular, Minniti et al. (2020) used classical Cepheids to trace the gradient of metallicity on both sides of the Galactic disk.

Taking advantage of the Gaia data, Hao et al. (2020, 2022) found more than 700 new OC candidates. Although over 2000 new OCs have been identified, the Galactic OC census is considered to be far from complete and this particularly true along heavily obscured sightlines (e.g., Cantat-Gaudin et al. 2020; Castro-Ginard et al. 2022). Thanks to the VISTA variables in the Vía Láctea survey and other longer wavelength surveys, classical Cepheids are being discovered on the other side of the Galaxy (e.g., Dékány et al. 2015b; Minniti et al. 2021). At present, more than 3000 classical Cepheids have been discovered in the Milky Way (e.g., Chen et al. 2020; Pietrukowicz et al. 2021).

Despite the significance of OC–Cepheid pairs and the discovery of thousands of OCs and classical Cepheids, only a few OCs harboring Cepheids have been reported so far. The presence of Cepheids in Galactic OCs was first discovered by Doig (1925, 1926) and probably others, as noted by Fernie (1969). By cross-matching the known OCs in Gaia with a list of Cepheids, Zhou & Chen (2021) reported the largest OC–Cepheid sample to date, in which only 33 Cepheids exist in 29 OCs. Therefore, the identification of more OC–Cepheid pairs is an indispensable undertaking with great importance.

The latest Gaia data release, DR3, includes astrometric and photometric measurements of 1.8 billion objects, and the determination of the mean radial velocity (RV) of 33 million objects (Gaia Collaboration 2022b). Meanwhile, the data quality of Gaia has been further improved. Previous studies identifying OC–Cepheid pairs mainly cross-matched Cepheids with known OCs. However, as mentioned above, the OC census of the Milky Way is incomplete. Hence, based on the aforementioned comprehensive datasets, this work aims to find OC–Cepheid pairs in the regions of Cepheids using the cluster-search method employed in our previous studies (Hao et al. 2020, 2022).

We organize this paper as follows. The data used in this work are described in Sect. 2. Section 3 presents the method applied to the search for OCs holding classical Cepheids. The results are displayed in Sect. 4 and we summarize this work in Sect. 5.

2 Data

The classical Cepheids adopted in this work come from the catalog compiled by Pietrukowicz et al. (2021), who carefully inspected candidate Cepheids from many surveys based on their long-term experience in the field of variable stars. The surveys related to the sources in this catalog are as follows: General catalog of Variable Stars, All Sky Automated Survey, Northern Sky Variability Survey, Optical Gravitational Lensing Experiment, Asteroid Terrestrial-impact Last Alert System, All-Sky Automated Survey for Supernovae, Wide-field Infrared Survey Explorer, Gaia astrometric mission, VISTA Variables in the Via Lactea, Zwicky Transient Facility, and so on (for more details, see Pietrukowicz et al. 2021, and references within). This catalog, with a purity exceeding 97%, contains a total of 3 352 classical Cepheid variables, which are available at the OGLE Internet Data Archive1. In addition, Pietrukowicz et al. (2021) not only presented detailed information on these classical Cepheids but also cross-matched them with the Gaia EDR3 catalog. In this work, considering that the present-day uncertainty on parallax in the Gaia area can reach 0.02−0.03 mas (Lindegren et al. 2021), we decided to mainly concentrate on the objects within 5 kpc of the Sun, resulting in a subsample that contains 1 085 classical Cepheids. For each of these Cepheids, we have obtained detailed information such as astrometric and photometric parameters in the Gaia DR3 dataset.

We used the dataset from the Gaia DR3 catalog, which can be found at the Gaia archive2. The astrometric and photometric parameters of stars in Gaia DR3 are the same as those in Gaia EDR3 (Gaia Collaboration 2020), which contain stellar celestial positions, parallaxes, proper motions (l, b, ϖ, µα*, and µδ), and three photometric bands (G, GBP, and GRP). Meanwhile, Gaia DR3 includes 33 million objects with new determinations on their mean radial velocities. As in many previous works searching for OCs (e.g., Cantat-Gaudin et al. 2018; Hao et al. 2022; Castro-Ginard et al. 2022), we only chose sources brighter than G = 18. For sources with five-parameter solutions at this magnitude, the median parallax uncertainty on them is 0.120 mas, and the median uncertainties on proper motions in µα* and µδ are 0.123 and 0.111 mas yr−1, respectively (Lindegren et al. 2021). We also rejected sources that present negative parallaxes or large proper motions with |µα*|, |µδ| > 30 mas yr−1. As the aim of this work is to find OCs containing classical Cepheids, we only focus on the 1085 spatial regions where selected Cepheids are present; our experimental setup is described in the following section.

3 Methods

The sample-based clustering search method was successfully adopted in our previous studies and we found more than 700 new OC candidates in the Gaia area (Hao et al. 2020, 2022). This method consists primarily of four steps: determining spatial regions, searching for stellar clusters in the five-dimensional (5D) parametric space (l, b, ϖ, µα*, and µδ), selecting targets based on proper-motion dispersion, and visually inspecting the multiple distributions. In this section, we describe the method and process of searching for OCs containing classical Cepheids.

Sample determination

Under the selection criteria of stars described in Sect. 2, we first extracted stars around the spatial locations of classical Cepheids. The key aspect of this step is to decide the spatial size and distance boundaries of the different regions where Cepheids exist. In this work, by referencing our previous experience in searching for OCs (Hao et al. 2020, 2022), the celestial sizes and distance boundaries of different spatial regions are set in the range of [1°, 10°] and [0.5, 5] kpc, respectively, mainly according to the distances of classical Cepheids and the spatial density of regions. Thus, we obtained 1085 original stellar samples. In addition, if in the following search there are star clusters located at the borders of regions, we change the sizes of the regions and re-detect these star clusters.

Clustering algorithm

The powerful unsupervised clustering algorithm, known as the density-based spatial clustering of applications with noise (DBSCAN, Ester et al. 1996), was applied in the multidimensional space (i.e., l, b, ϖ, µα*, and µδ) to find spatial over-density structures in the regions. For DBSCAN, the parameters ϵ and minPts are needed, which specify the radius of a neighborhood with respect to some point and the desired minimum cluster size, respectively. We let the parameter minPts vary in the range of six to ten stars. The corresponding parameters ϵ of each sample were then automatically calculated using the z-score standardization, a Gaussian kernel-density-estimation method (Lampe et al. 2011), and the k-nearest (k = minPts – 1) neighbors algorithm (Altman 1992); see Hao et al. (2020) for elaborate details on the determination of ϵ. After preparing the parameter pairs for all samples, we conducted the search for statistical stellar clusters in all regions of existing classical Cepheids.

We noticed that there are some classical Cepheids located close to each other in spatial locations, which means regional overlaps exist among some initial samples, resulting in several stellar clusters being found repeatedly. We resolved this case by identifying those stellar clusters that have comparable mean parameters within 3σi (σ is the standard deviation, i = l, b, ϖ, µα*, and µδ) in the 5D parametric space. Subsequently, we visually inspected these stellar clusters and filtered out the repetitions.

Selection based on proper-motion dispersion

The apparent proper-motion dispersion of a stellar cluster can serve as an effective means to distinguish reasonable OCs from implausible ones, and this technique has been used in some work to select potentially real OCs from spatial over-density structures (e.g., Hunt & Reffert 2021; Hao et al. 2022). The gravitationally bound OC systems generally have small internal velocity dispersions. Here, as in our previous work (Hao et al. 2022), we selected the possibly genuine OCs using the following proper-motion criterion: σμa*2+σμδ2{ 0.5mas yr1if ϖ < 1mas22ϖ4.7404mas yr1if ϖ 1mas. $\sqrt {\sigma _{{{\rm{\mu }}_{a*}}}^2 + \sigma _{{{\rm{\mu }}_\delta }}^2} \le \left\{ {\matrix{{0.5\,{\rm{mas y}}{{\rm{r}}^{ - 1}}} \hfill &amp; {{\rm{if}}\,\varpi \,{\rm{ &lt; }}\,{\rm{1mas}}} \hfill \cr {2\sqrt 2 {\varpi \over {4.7404}}{\rm{mas}}\,{\rm{y}}{{\rm{r}}^{ - 1}}} \hfill &amp; {{\rm{if}}\,\varpi \, \ge \,{\rm{1mas}}} \hfill \cr } .} \right.$(1)

Here, σμα*${\sigma _{{{\rm{\mu }}_{\alpha *}}}}$, σμδ${\sigma _{{{\rm{\mu }}_\delta }}}$, and ϖ are the standard deviations of two proper motions and the mean parallax of statistical stellar clusters. For the property of proper-motion dispersions of OCs identified in the Gaia area, we recommend Cantat-Gaudin & Anders (2020).

Visual inspection and confirmation

For the potentially real OCs, we visually inspected their multiple characters, such as their members, the density relative to the background stellar distribution, the distribution of radial velocities when available, and especially the CMDs. As the member stars of an OC are almost simultaneously formed in the same molecular cloud, their distribution on the CMD can exhibit the feature of some empirical isochrone. These measures were mainly taken for the potentially genuine OCs that have not been reported to select the more reliable OC candidates.

We computed the mean RV, and the weighted standard deviation of RV for each OC candidate using the method described in Soubiran et al. (2018). We also estimated the ages, line-of-sight extinction (absorption, AG), and distance modulus (DM) of the proposed OC candidates using the photometric parameters (G, GBP, and GRP) provided in Gaia to construct the theoretical isochrones in their CMDs. It should be noted that the differential reddening seen in some OCs can complicate matters considerably, in that a tight obvious main sequence and evolved pattern may not be seen. In this case, differential reddening procedures need to be employed, potentially using a color–color diagram (e.g., Turner 1994, see Figs. 2 and 3 for NGC 1545). The isochrones we adopted are derived from the PARSEC library (Bressan et al. 2012), which has been updated for the Gaia passbands with the photometric calibration of Evans et al. (2018) and contains logarithmic ages (i.e., log(age/year)) ranging from 5.92 to 10.13 and metal fractions (z) ranging from 0.015 to 0.029, as well as a series of isochrones that are generated with steps of Δlog(age/year) = 0.02 and Δz = 0.001. We used the least square fitting method to obtain an isochrone that is as reliable as possible for each OC candidate. Meanwhile, the extinction and reddening that we adopted are corrected using an extinction law of Rv = 3.1 (O'Donnell 1994) for the Milky Way galaxy, and the applied approximate relation is E(GBP – GRP) = 0.50 Av as recommended by Andrae et al. (2018). We refer the reader to Sect. 3.4 in Hao et al. (2022) for more details on the determinations.

4 Results

After implementing the processes described in Sect. 3, we cross-matched the stellar clusters found in this work with reported OCs to identify the OCs that are already known and those that are new findings. In total, 635 known OCs were re-detected, and we proposed 38 stellar clusters as newly found OC candidates. The classical Cepheids used in this work were cross-matched with the above OCs, and we found 45 OCs harboring a total of 50 Cepheids. For the newly found OCs, we determined their parameters and investigated their characteristics.

4.1 Cross-matching with known OCs

The stellar clusters found in this work were first cross-matched with previously reported OCs in the Gaia area: Cantat-Gaudin et al. (2018, 2019), Castro-Ginard et al. (2018, 2019, 2020, 2022), Sim et al. (2019), Liu & Pang (2019), Ferreira et al. (2019, 2020, 2021), Hao et al. (2020, 2022), He et al. (2021), and Hunt & Reffert (2021). In general, two stellar clusters can be considered a match if their mean parameters are within 3σi (i = l, b, ϖ, µα* and µδ) in the 5D parametric space. Thus, in this step, we obtained 594 previously known OCs in the Gaia data.

Subsequently, stellar clusters obtained by us were cross-matched with OCs that were known before Gaia, which mainly derive from the catalogs compiled by Dias et al. (2002, V3.5, 2015 edition) and Kharchenko et al. (2013), which contain approximately 2000 and 3000 objects, respectively, while they have repetitive sources. The criterion used in this step is the same as in our previous work (Hao et al. 2022). We first selected star clusters that fall within a circle of 0.5° radii of objects listed in these two catalogs. If the mean proper motions (µα* and µδ) of the listed objects are within three times the standard deviations of the stellar clusters found in this work, and their distances are compatible, they can be considered matched. Here, the distance of a star cluster found in this work is the inverse of the mean parallax, and the corresponding error comes from three times the standard deviation of the parallax. Cantat-Gaudin et al. (2018) took into account many objects in these two catalogs, but we re-detected 24 known OCs in Dias et al. (2002) and 21 known OCs in Kharchenko et al. (2013), respectively, and 4 OCs among them are repetitive. Similar to the OC catalogs, we also cross-matched our findings to the globular cluster catalog reported by Baumgardt et al. (2019), but we did not explore any coincidences.

For the remaining stellar clusters found in this work that were not yet cross-matched with known OCs, we carefully inspected their multidimensional nature as described in Sect. 3 and proposed 38 reliable ones to be OC candidates, which are numbered from OC–0705 in order to be the continuation of our previous work (Hao et al. 2022). Combining these 38 newly found OC candidates with 594 OCs in Gaia and 41 OCs in the catalogs before Gaia, we searched for a total of 673 OCs in this work. Subsequently, we investigated these OCs to identify OCs hosting classical Cepheids.

4.2 Open clusters housing classical Cepheids

We regarded OCs and classical Cepheids as combinations in two scenarios: Cepheids are member stars of OCs; or the 5D astrometric parameters (l, b, ϖ, µα*, and µδ) of Cepheids and OCs are concordant within three times the corresponding standard deviations. As described in Sect. 2, we obtained Gaia DR3 source_id for each classical Cepheid. For the first case, we cross-matched the list of classical Cepheids with the member stars of the 673 OCs found in this work. The cross-match was performed using Gaia DR3 source_id and returned 30 OCs holding 32 Cepheids. The 5D astrometric parameters of the remaining classical Cepheids were then compared with all found OCs, resulting in 18 Cepheids associated with 18 OCs, of which 3 OCs are the same ones as above. Thus, we identified 45 OCs housing 50 classical Cepheids, and the names of these OC–Cepheid pairs are listed in Table A.1.

We cross-matched the OC–Cepheid pairs found in this work with those in Zhou & Chen (2021), who reported the discovery of 29 OCs containing 33 Cepheids. Finally, we were able to re-detect 27 OCs holding 30 Cepheids. We noticed that the differences in proper motion in µδ between the remaining two OCs (Ruprecht 97 and UBC 375) and their two associated classical Cepheids proposed by Zhou & Chen (2021) are significantly large (10.9σ and 7.3σ, respectively), resulting in the two pairs not being identified in this work. There is another Cepheid (V379 Cas) not being matched by us because its spatial location is about one degree away from the target OC (NGC 129) in the sky. The dynamical evolution proposed by Turner (1985) suggests that instances of Cepheids may be possible members of cluster coronae. Coulson & Caldwell (1985) and Majaess et al. (2013) both reported that the Cepheid QZ Nor is a coronal member of the OC NGC 6067. Therefore, the Cepheid V379 Cas may be an ejected member of the NGC 129 cluster. Although V379 Cas is not followed up here, it is still worth investigating. We also recovered some interesting OCs that contain more than one classical Cepheid: for example, the widely known Galactic OC NGC 7790 hosting three Cepheids, and the OC Kronberger 84 holding two Cepheids reported by Zhou & Chen (2021). The cross-match using Gaia DR3 source_id shows that the Cepheid V733 Cyg is a non-member of Kronberger 84, but that the 5D astrometric parameters of V733 Cyg are consistent with the cluster, within three times the corresponding standard deviations.

For the 45 OCs containing classical Cepheids, 35 are previously known OCs compiled by Cantat-Gaudin et al. (2020), two OCs (OC–0125 and OC–0675) were reported by Hao et al. (2022), one OC (Alessi 95) was reported by Turner et al. (2012b) and Majaess et al. (2012), one OC (UBC 1424) was reported by Castro-Ginard et al. (2022), and the last six are newly discovered in this work. Based on the member stars of the 45 OCs in Gaia DR3, we determined their mean astrometric parameters (α, δ, ϖ, µα*, and µδ) and the corresponding standard deviations, as presented in Table A.1. The apparent angular radius (θ) of each OC was estimated as the square root of the quadratic sum of σl and σb (e.g., Castro-Ginard et al. 2020, 2022; Hao et al. 2022). The cluster apparent diameters are then determined as two times the θ. The ages and uncertainties of the known OCs are from the previous determinations in Sana et al. (2010), Turner et al. (2012b), Cantat-Gaudin et al. (2020), Hao et al. (2022), and Castro-Ginard et al. (2022), and those of newly found OCs were estimated by the method described in Sect. 3. Table A.1 also presents the astrometric parameters (α, δ, ϖ, µα*, and µδ) of Classical Cepheids, which are from the dataset of Gaia DR3, and the listed periods of Cepheids are from Pietrukowicz et al. (2021). The Cepheid ages listed in Table A.1 are taken from the mean calculation of the Cepheid PAR reported in following three works: Efremov (2003), Bono et al. (2005), and Turner et al. (2012a), and the deviation is cited as the uncertainty. We note that Efremov (2003) and Bono et al. (2005) only present the PAR for fundamental and first-overtone classical Cepheids, and so the second-overtone pulsating mode of the Cepheid OGLE-GD-CEP-1673 means that its age is unavailable from the two works, and its age uncertainty is not listed in Table A.1. However, it should be noted that the obtained ages of Cepheids are uncertain parameters, and are affected by reddening, distance, photometric calibrations, and so on, as pointed out by Bono et al. (2005). We also estimated the angular separation between the classical Cepheid and its host OC, where the cluster centre is derived from the mean astrometric positions of member stars.

We noticed that some OCs may exhibit pre-main sequence (PMS) characteristics, and therefore would be unrelated to Cepheids, particularly the lower period ones. For example, the readily discernible PMS of Trumpler 14 was studied by Sana et al. (2010), and the presence of O-type stars in NGC 6193 was reported by Skinner et al. (2005), suggesting that the two OCs are not related to Cepheids. Indeed, OCs with ages younger than ten million years are unlikely to be related to Cepheids, and those with ages near ten million years also tend to be related to very long-period Cepheids (e.g., Anderson et al. 2016). We investigated all the identified OC–Cepheid pairs by considering the OC ages and the periods of Cepheids, and found that 11 Cepheids may be unrelated to OCs. However, it is interesting that the astrometric matches of these Cepheids and their host OCs still emerge, and so these Cepheids are worthy of follow-up investigation in order to come to a firm conclusion on their nature. In Table A.1, class A Cepheids are related to OCs and have been identified as cluster members, Cepheids in class B are related to OCs but have not been identified as cluster members, and class C Cepheids may be unrelated to OCs but they possess similar astrometric parameters. Table A.1 can be found online at the CDS, together with the table reporting the member stars of OCs harboring classical Cepheids.

Figure 1 shows three examples of OCs holding classical Cepheids, including the members of OCs and Cepheids in the astrometric spaces as well as CMDs. The remaining cases are presented in Fig. A.1. As intermediate-mass stars, classical Cepheids have evolved away from the main sequence and are in the stage of shell hydrogen burning or core helium burning (e.g., Turner et al. 2006). Hence, in the CMDs of bona fide OCs harboring classical Cepheids, Cepheids are generally located on the classical instability strip and are brighter than the main sequence stars. We visually inspected the CMDs of the OC–Cepheid pairs identified in this work and found 38 Cepheids with the above features. For the last 12 classical Cepheids that may not be on the classical instability strip of their host OCs, these classical Cepheids require further investigation, and we associated asterisks with the names of these classical Cepheids in Table A.1.

Figure 2 presents the distribution of the 45 OCs hosting classical Cepheids in Galactic l and b coordinates. The vast majority (~ 93%) of them are located within the Galactic disk at latitudes of |b| ≤ 5°. We also made a comparison between the parallaxes of OCs containing classical Cepheids found in this work and those previously reported in Cantat-Gaudin et al. (2020), Castro Ginard et al. (2022), and Hao et al. (2022), as shown in Fig. 3. The mean parallaxes of the known OCs re-detected in this work are in good agreement with previously published values. However, compared with the previous results, the smaller errors of parallaxes obtained here demonstrate the improved data quality of Gaia DR3 over Gaia data release 2. We also noticed that 43 of these OCs have RV measurements, of which 40 are based on more than one member star, and 35 are based on more than two member stars. The OCs marked with asterisks in Table A.1 have an RV measurement. Therefore, in addition to providing useful indications for optimizing the zero-point of the Cepheid period–luminosity relationship and period–age relation, it is hopeful that these OCs can be used to unveil the nature and especially the dynamical evolution of classical Cepheids.

thumbnail Fig. 1

Examples of OCs (blue dots) harboring classical Cepheids (red hexagram). The columns of each panel represent the distributions of the member stars of OCs and classical Cepheids for position in RA (α) and Dec (δ), proper motions in µα* and µδ, and the CMD in G vs. GBPGRP, as well as the mean parallax and total proper-motion dispersion of OCs. Here, the listed OCs are Berkeley 58, NGC 103, and NGC 129.

thumbnail Fig. 2

Distributions in Galactic coordinates of the OCs housing classical Cepheids (red dots), newly found OCs (blue triangles) and previously known OCs (black plus) obtained in this work.

4.3 Newly found OCs

As described in Sect. 4.1, aside from re-detecting 635 previously known OCs, we are fortunate to propose 38 reliable ones as new OC candidates in this work. Figure 4 displays the distributions of some examples of newly found OCs. Here we present the features of their member stars in astrometric spaces and CMDs, indicating that these proposed new OCs show high concentrations of member stars in the five astrometric parameters, and the members in the CMDs are concentrated on potentially theoretical isochrones. Table A.2 contains a list of the parameters of the new OC candidates, including their mean astrometric parameters (α, δ, l, b, ϖ, µα*, and µδ) and standard deviations, RVs, and the weighted standard deviations of RVs when available, as well as the estimated ages, line-of-sight extinction (AG), and DM. In addition, the apparent angular sizes of these new OC candidates were computed as the square root of the quadratic sum of σl and σb following Castro-Ginard et al. (2020, 2022), and Hao et al. (2022). The table can be found online at the CDS, including the table of member stars.

The new OC candidates are generally at larger distances, of which ~95% are further than 1 kpc, and only two OCs are located within 1 kpc. For the newly found OCs, 31 of them have member stars with RV s in Gaia DR3, more than half of which contain more than one star with an RV measurement, and the names of these OCs have been placed with asterisks in Table A.2. The distributions of these OCs in the Galactic plane are also shown in Fig. 2. The newly found OCs are mostly found at low latitudes within the Galactic disk, with the vast majority of new OC candidates (~97%) located at |b| ≤ 4°, and only one at b ≈ 7°.

thumbnail Fig. 3

Parallactic comparison between OCs containing classical Cepheids found in this work and those previously reported in Cantat-Gaudin et al. (2020, TCG20), Castro-Ginard et al. (2022, ACG22), and Hao et al. (2022, HAO22).

thumbnail Fig. 4

Examples of newly found OCs in this work. The columns of each panel are the same as Fig. 1. Here, the listed OCs are OC–0707, OC–0709, OC–0719, OC–0724, OC–0727, OC–0734, OC–0737, OC–0738, and OC–0739.

5 Summary

Combining the high-precision data provided by Gaia DR3 and known classical Cepheids with our previously practiced clustering search method, we conducted a search for OCs containing classical Cepheids. We identified 50 classical Cepheids associated with 45 OCs in 5D astrometric space. After investigating the ages of the OCs and the periods of the Cepheids, we found 39 OC–Cepheid pairs to be probable, and a remaining 11 Cepheids to be improbable and possibly unrelated to OCs but worthy of follow-up investigation. For all of these OCz–Cepheid pairs, we provide their astrometric and astrophysical parameters. In addition to re-detecting 635 known OCs in the literature, we also propose 38 newly found OC candidates within the Galactic disk based on a detailed inspection of the multidimensional distributions of their member stars. For each of these new OC candidates, we present its mean astrometric values and radial velocity when available, together with estimations of age, line-of-sight extinction, and distance modulus. The new OCs are generally located at greater distances from the Sun, which can help to improve the completeness of the Galactic OC census.

Our current search mainly focuses on the regions within 5 kpc of the Sun. Hence, it is expected that more OC–classical Cepheid pairs will be detected in our Galaxy along with the improvements in parallax and proper motion precision of the future data releases of Gaia.

Acknowledgements

We sincerely appreciate the anonymous referee for the instructive comments which help us to improve the paper. This work was funded by NSFC Grants 11933011, 11873019, 12203104, and the Key Laboratory for Radio Astronomy. Y.J.L. thanks support from the Natural Science Foundation of Jiangsu Province (grant number BK20210999) and the Entrepreneurship and Innovation Program of Jiangsu Province. Our 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.

Appendix A Additional tables and figures

Table A.1

Parameters of the OC-Cepheid pairs identified in this work.

Table A.2

Parameters of the proposed newly found OCs ordered by increasing l.

thumbnail Fig. A.1

Open clusters (blue dots) harboring classical Cepheids (red hexagram). The columns of each panel represent the distributions of the member stars of OCs and classical Cepheids for position in RA (α) and Dec (δ), proper motions in µα* and µδ, and the CMD in G vs. GBPGRP, as well as the mean parallax and total proper-motion dispersion of OCs. Here, the listed OCs are Alessi 95, FSR 0951, vdBergh 1, OC-0717, UBC 229, UBC 1424, Ruprecht 79, NGC 5662, NGC 6067, NGC 6087, OC-0675, and IC 4725.

thumbnail Fig. A.1

continued. Here, the listed OCs are NGC 6649, NGC 6664, FSR 0158, UBC 129, UBC 135, Berkeley 51, OC-0125, Berkeley 55, Kronberger 84, NGC 7790, OC-0718, and OC-0721.

thumbnail Fig. A.1

continued. Here, the listed OCs are OC-0723, UBC 491, UBC 290, BH 222, UBC 106, Trumpler 35, UBC 130, FSR 0172, OC-0715, OC-0720, IC 2395, and Trumpler 14.

thumbnail Fig. A.1

ccontinued. Here, the listed OCs are UBC 266, BH 121, Lynga 6, NGC 6193, UBC 323, and Czernik 41.

References

  1. Aiola, S., Calabrese, E., Maurin, L., et al. 2020, JCAP, 2020, 047 [Google Scholar]
  2. Altman, N. S. 1992, BAN, 46, 175 [Google Scholar]
  3. Anderson, R. I., Saio, H., Ekström, S., Georgy, C., & Meynet, G. 2016, A&A, 591, A8 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  4. Andrae, R., Fouesneau, M., Creevey, O., et al. 2018, A&A, 616, A8 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  5. Baumgardt, H., Hilker, M., Sollima, A., & Bellini, A. 2019, MNRAS, 482, 5138 [Google Scholar]
  6. Bono, G., Marconi, M., Cassisi, S., et al. 2005, ApJ, 621, 966 [NASA ADS] [CrossRef] [Google Scholar]
  7. Bressan, A., Marigo, P., Girardi, L., et al. 2012, MNRAS, 427, 127 [NASA ADS] [CrossRef] [Google Scholar]
  8. Breuval, L., Kervella, P., Anderson, R. I., et al. 2020, A&A, 643, A115 [EDP Sciences] [Google Scholar]
  9. Cantat-Gaudin, T., & Anders, F. 2020, A&A, 633, A99 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  10. Cantat-Gaudin, T., Jordi, C., Vallenari, A., et al. 2018, A&A, 618, A93 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  11. Cantat-Gaudin, T., Krone-Martins, A., Sedaghat, N., et al. 2019, A&A, 624, A126 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  12. Cantat-Gaudin, T., Anders, F., Castro-Ginard, A., et al. 2020, A&A, 640, A1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  13. Castro-Ginard, A., Jordi, C., Luri, X., et al. 2018, A&A, 618, A59 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  14. Castro-Ginard, A., Jordi, C., Luri, X., Cantat-Gaudin, T., & Balaguer-Nunez, L. 2019, A&A, 627, A35 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  15. Castro-Ginard, A., Jordi, C., Luri, X., et al. 2020, A&A, 635, A45 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  16. Castro-Ginard, A., Jordi, C., Luri, X., et al. 2022, A&A, 661, A118 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  17. Chen, X., de Grijs, R., & Deng, L., 2017, MNRAS, 464, 1119 [NASA ADS] [CrossRef] [Google Scholar]
  18. Chen, X., Wang, S., Deng, L., et al. 2019, NatAs, 3, 320 [Google Scholar]
  19. Chen, X., Wang, S., Deng, L., et al. 2020, ApJS, 249, 18 [NASA ADS] [CrossRef] [Google Scholar]
  20. Coulson, I. M., & Caldwell, J. A. R. 1985, MNRAS, 216, 671 [NASA ADS] [Google Scholar]
  21. De Somma, G., Marconi, M., Cassisi, S., et al. 2020, MNRAS, 496, 5039 [CrossRef] [Google Scholar]
  22. Dékány, I., Minniti, D., Majaess, D., et al. 2015b, ApJ, 812, L29 [CrossRef] [Google Scholar]
  23. 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]
  24. Dinnbier, F., Anderson, R. I., & Kroupa, P. 2022, A&A, 659, A169 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  25. Doig, P. 1925, JBAA, 35, 202 [Google Scholar]
  26. Doig, P. 1926, JBAA, 36, 60 [Google Scholar]
  27. Efremov, Y. N. 2003, ARep, 47, 1000 [NASA ADS] [Google Scholar]
  28. Ester, M., Kriegel, H.-P., Sander, J., & Xu X. 1996, 2nd International Conference on Knowledge Discovery and Data Mining (AAAI Press), 96, 226 [Google Scholar]
  29. Evans, D. W., Riello, M., De Angeli, F., et al. 2018, A&A, 616, A4 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  30. Feast, M. W. 1999, PASP, 111, 775 [NASA ADS] [CrossRef] [Google Scholar]
  31. Ferreira, A. F., Santos, J. F. C., Corradi, W. J. B., Maia, F. F. S., & Angelo, M. S. 2019, MNRAS, 483, 5508 [NASA ADS] [CrossRef] [Google Scholar]
  32. Ferreira, A. F., Corradi, W. J. B., Maia, F. F. S., Angelo, M. S., & Santos, J. F. C. 2020, MNRAS, 496, 2021 [NASA ADS] [CrossRef] [Google Scholar]
  33. Ferreira, A. F., Corradi, W. J. B., Maia, F. F. S., Angelo, M. S., & Santos, J. F. C. 2021, MNRAS, 502, L90 [NASA ADS] [CrossRef] [Google Scholar]
  34. Fernie, J. D. 1969, PASP, 81, 707 [NASA ADS] [CrossRef] [Google Scholar]
  35. Freedman, W. L., Madore, B. F., Hatt, D., et al. 2019, ApJ, 882, 34 [Google Scholar]
  36. Freedman, W. L., Madore, B. F., Hoyt, T., et al. 2020, ApJ, 891, 57 [Google Scholar]
  37. Fry, A. M., & Carney, B. W. 1997, AJ, 113, 1073 [NASA ADS] [CrossRef] [Google Scholar]
  38. Gaia Collaboration (Brown, A. G. A., et al.) 2018, A&A, 616, A1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  39. Gaia Collaboration (Brown, A. G. A., et al.) 2020, A&A, 649, A1 [Google Scholar]
  40. Gaia Collaboration (Vallenari, A., et al.) 2022, A&A, in press https://doi.org/10.1051/0004-6361/202243940 [Google Scholar]
  41. Hao, C. J., Xu, Y., Wu, Z. Y., et al. 2020, PASP, 132, 034502 [NASA ADS] [CrossRef] [Google Scholar]
  42. Hao, C. J., Xu, Y., Hou, L. G., et al. 2021, A&A, 652, A102 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  43. Hao, C. J., Xu, Y., Wu, Z. Y., et al. 2022, A&A, 660, A4 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  44. He, Z. H., Xu, Y., Hao, C. J., Wu, Z. Y., & Li, J. J. 2021, Res. Astron. Astrophys. 21, 093 [Google Scholar]
  45. Hunt, E. L., & Reffert, S. 2021, A&A, 646, A104 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  46. 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]
  47. Lampe, O. D., & Hauser, H. 2011, IEEE Pacific Visualization Symposium (Piscataway, NJ: IEEE), 171 [Google Scholar]
  48. Leavitt, H. S., & Pickering, E. C. 1912, Harv. Coll. Obs. Circ., 173, 1 [Google Scholar]
  49. Lemasle, B., Groenewegen, M. A. T., Grebel, E. K., et al. 2017, A&A, 608, A85 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  50. Lindegren, L., Hernandez, J., Bombrun, A., et al. 2018, A&A, 616, A2 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  51. Lindegren, L., Klioner, S. A., Hernández, J. et al. 2021, A&A, 649, A2 [EDP Sciences] [Google Scholar]
  52. Liu, L., & Pang, X. 2019, ApJS, 245, 32 [CrossRef] [Google Scholar]
  53. Majaess, D. J., Turner, D. G., Gallo, L., et al. 2012, ApJ, 753, 144 [NASA ADS] [CrossRef] [Google Scholar]
  54. Majaess, D. J., Sturch, L., Moni Bidin, C., et al. 2013, Ap&SS, 347, 61 [NASA ADS] [CrossRef] [Google Scholar]
  55. Majaess, D. J. 2020, ApJ, 897, 13 [NASA ADS] [CrossRef] [Google Scholar]
  56. Madore, B. F., & Freedman, W. L. 1991, PASP, 103, 933 [Google Scholar]
  57. Minniti, J. H., Sbordone, L., Rojas-Arriagada, A., et al. 2020, A&A, 640, A92 [EDP Sciences] [Google Scholar]
  58. Minniti, J. H., Zoccali, M., Rojas-Arriagada, A., et al. 2021, A&A, 654, A138 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  59. O’Donnell, J. E. 1994, ApJ, 422, 158 [Google Scholar]
  60. Pietrukowicz, P., Soszyński, I., & Udalski, A. 2021, Acta Astron., 71, 205 [NASA ADS] [Google Scholar]
  61. Planck Collaboration VI. 2020, A&A, 641, A6 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  62. Poggio, E., Drimmel, R., Cantat-Gaudin, T., et al. 2021, A&A, 651, A104 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  63. Riess, A. G., Casertano, S., Yuan, W., Macri, L. M., & Scolnic, D. 2019, ApJ, 876, 85 [Google Scholar]
  64. Riess, A. G., Casertano, S., Yuan, W., et al. 2021, ApJ, 908, L6 [NASA ADS] [CrossRef] [Google Scholar]
  65. Riess, A. G., Breuval, L., Yuan, W., et al. 2022, ApJ, 938, 36 [NASA ADS] [CrossRef] [Google Scholar]
  66. Sana, H., Momany, Y., Gieles, M., et al. 2010, A&A, 515, A26 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  67. Sim, G., Lee, S. H., Ann, H. B., & Kim, S. 2019, J. Korean Astron. Soc., 52, 145 [NASA ADS] [Google Scholar]
  68. Skinner, S. L., Zhekov, S. A., Palla, F., & Barbosa, C. L. D. R. 2005, MNRAS, 361, 191 [NASA ADS] [CrossRef] [Google Scholar]
  69. Soubiran, C., Cantat-Gaudin, T., Romero-Gómez, M., et al. 2018, A&A, 619, A155 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  70. Turner, D. 1985, in IAU Colloq. 82, Cepheids: Theory and Observation, ed. B. F. Madore (New York: Cambridge University Press), 209 [Google Scholar]
  71. Turner, D. G. 1994, J. Roy. Astron. Soc. Can., 88, 176 [NASA ADS] [Google Scholar]
  72. Turner, D. G., & Burke, J. F. 2002, AJ, 124, 2931 [NASA ADS] [CrossRef] [Google Scholar]
  73. Turner, D. G., Abdel-Sabour Abdel-Latif, M., & Berdnikov, L. N. 2006, PASP, 118, 410 [Google Scholar]
  74. Turner, D. G., Forbes, D., English, D., et al. 2008, MNRAS, 388, 444 [NASA ADS] [CrossRef] [Google Scholar]
  75. Turner, D. G. 2012a, JAAVSO, 40, 502 [NASA ADS] [Google Scholar]
  76. Turner, D. G., Majaess, D. J., Lane, D. J., et al. 2012b, MNRAS, 422, 2501 [NASA ADS] [CrossRef] [Google Scholar]
  77. Zhou, X., & Chen, X. 2021, MNRAS, 504, 4768 [NASA ADS] [CrossRef] [Google Scholar]

All Tables

Table A.1

Parameters of the OC-Cepheid pairs identified in this work.

Table A.2

Parameters of the proposed newly found OCs ordered by increasing l.

All Figures

thumbnail Fig. 1

Examples of OCs (blue dots) harboring classical Cepheids (red hexagram). The columns of each panel represent the distributions of the member stars of OCs and classical Cepheids for position in RA (α) and Dec (δ), proper motions in µα* and µδ, and the CMD in G vs. GBPGRP, as well as the mean parallax and total proper-motion dispersion of OCs. Here, the listed OCs are Berkeley 58, NGC 103, and NGC 129.

In the text
thumbnail Fig. 2

Distributions in Galactic coordinates of the OCs housing classical Cepheids (red dots), newly found OCs (blue triangles) and previously known OCs (black plus) obtained in this work.

In the text
thumbnail Fig. 3

Parallactic comparison between OCs containing classical Cepheids found in this work and those previously reported in Cantat-Gaudin et al. (2020, TCG20), Castro-Ginard et al. (2022, ACG22), and Hao et al. (2022, HAO22).

In the text
thumbnail Fig. 4

Examples of newly found OCs in this work. The columns of each panel are the same as Fig. 1. Here, the listed OCs are OC–0707, OC–0709, OC–0719, OC–0724, OC–0727, OC–0734, OC–0737, OC–0738, and OC–0739.

In the text
thumbnail Fig. A.1

Open clusters (blue dots) harboring classical Cepheids (red hexagram). The columns of each panel represent the distributions of the member stars of OCs and classical Cepheids for position in RA (α) and Dec (δ), proper motions in µα* and µδ, and the CMD in G vs. GBPGRP, as well as the mean parallax and total proper-motion dispersion of OCs. Here, the listed OCs are Alessi 95, FSR 0951, vdBergh 1, OC-0717, UBC 229, UBC 1424, Ruprecht 79, NGC 5662, NGC 6067, NGC 6087, OC-0675, and IC 4725.

In the text
thumbnail Fig. A.1

continued. Here, the listed OCs are NGC 6649, NGC 6664, FSR 0158, UBC 129, UBC 135, Berkeley 51, OC-0125, Berkeley 55, Kronberger 84, NGC 7790, OC-0718, and OC-0721.

In the text
thumbnail Fig. A.1

continued. Here, the listed OCs are OC-0723, UBC 491, UBC 290, BH 222, UBC 106, Trumpler 35, UBC 130, FSR 0172, OC-0715, OC-0720, IC 2395, and Trumpler 14.

In the text
thumbnail Fig. A.1

ccontinued. Here, the listed OCs are UBC 266, BH 121, Lynga 6, NGC 6193, UBC 323, and Czernik 41.

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.