Open Access
Issue
A&A
Volume 654, October 2021
Article Number A122
Number of page(s) 15
Section Galactic structure, stellar clusters and populations
DOI https://doi.org/10.1051/0004-6361/202141366
Published online 22 October 2021

© P. A. B. Galli et al. 2021

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.

1. Introduction

The young stellar clusters of the Solar neighbourhood constitute a primary laboratory for the investigation of the early stages of star and planet formation, as well as to connect models to observations and visualise the local structure of the Galaxy. Significant progress has been made in recent years to discover and characterise nearby young stellar clusters thanks to the precise astrometry delivered by the Gaia satellite (see e.g. Cantat-Gaudin et al. 2018; Sim et al. 2019; Meingast et al. 2019; Castro-Ginard et al. 2020). One important contribution to the field has been the recognition of the χ1 Fornacis cluster (also known as Alessi 13) as one of a few open clusters of its distance (∼100 pc) and age (∼40 Myr), with more than a hundred members, making it a major constituent of the Solar neighbourhood (Zuckerman et al. 2019). In the following, we refer to the χ1 Fornacis cluster (XFOR) using the same terminology that is also used to name the young stellar groups in the Solar neighbourhood after either a constellation or the most prominent cluster member; we emphasise that both names used in the literature refer to the same cluster.

The XFOR cluster is listed in the catalogue of open clusters compiled by Dias et al. (2002), with only a few cluster members, as it was poorly studied in the pre-Gaia era. The cluster was reanalysed by Yen et al. (2018) based on data from the first data release of the Gaia catalogue (Gaia-DR1, Gaia Collaboration 2016) where the authors confirmed nine members. The subsequent study conducted by Cantat-Gaudin et al. (2018), using the second data release of the Gaia catalogue (Gaia-DR2, Gaia Collaboration 2018), identified 48 cluster members.

More recently, Zuckerman et al. (2019, hereafter ZKK2019) revisited the cluster using different selection criteria to search for cluster members based on Gaia-DR2 data and increased the sample of cluster members to 108 stars. The authors argue that the cluster is coeval and co-moving with the Tucana-Horologium (THA) and Columba (COL) young stellar associations by comparing their positions in the colour-magnitude diagram and assigned the age of about 40 Myr to the cluster. Thus, the cluster appears to be significantly younger than 525 Myr as reported previously by Kharchenko et al. (2013), but still in good agreement with the age of about 30 Myr estimated by Mamajek (2016) from the saturated X-ray emission of the stars. Another interesting finding reported by ZKK2019 is the abundance of M-dwarfs in the cluster with mid-infrared excess emission indicative of debris discs with warm dust that could be associated with rocky planet formation, making the XFOR cluster a promising target in the search for exoplanets.

The early installment of the third data release of the Gaia space mission (Gaia-EDR3, Gaia Collaboration 2021) represents a significant advance with respect to its predecessor delivering more precise astrometry and photometry, along with reduced systematic errors that allow for more precise studies. Given the close proximity and young age of the XFOR cluster, we believe that it will become an important cornerstone to many studies related to star formation from models to the physical properties of the stars (e.g. as done for the Pleiades). In this paper, we revisit the census of cluster members using Gaia-EDR3 data and combine our results with spectroscopic observations to revisit the cluster properties. This study is conducted as part of the Dynamical Analysis of Nearby Clusters (DANCe, Bouy et al. 2013) project and follows a similar methodology applied to other stellar clusters in this series (Olivares et al. 2019; Miret-Roig et al. 2019; Galli et al. 2020a,b, 2021).

The paper is organised as follows. In Sect. 2, we perform a new membership analysis of the cluster based on the methodology developed by our team for the DANCe project, which allows us to identify 164 cluster members and confirming most of the members previously identified by ZKK2019 as well as adding new ones to the list. Section 3 describes the results obtained from our observations to determine the spectral type of the stars from low-resolution spectroscopy. Then, in Sect. 4, we use our new sample of cluster members to revisit several properties of the cluster (distance, kinematics, age, spatial distribution, and initial mass function). Finally, we present our conclusions in Sect. 5.

2. Membership analysis

We selected the most likely members of the XFOR cluster based on Gaia-EDR3 data using the same methodology applied by our team to other young stellar groups (see e.g. Miret-Roig et al. 2019; Galli et al. 2020a,b, 2021). We summarise the main steps of our membership analysis in this section and refer the reader to the original papers (Sarro et al. 2014; Olivares et al. 2019) for further details on the methodology.

To begin with, we downloaded the Gaia-EDR3 catalogue in the region defined by 215°   ≤  l  ≤  260° and −68°   ≤  b  ≤   − 45° in Galactic coordinates. This input catalogue contains 2 057 639 sources and includes all cluster members previously identified by ZKK2019. The representation space (i.e. the set of observables) we use in the membership analysis includes both the astrometric and photometric features of the Gaia-EDR3 catalogue. It is defined by μα cos δ, μδ, ϖ, GRP, and G − GRP, and it also includes 1 633 156 sources of the input catalogue with complete data. We do not use the GBP magnitude in our analysis because the flux in this band is likely to be overestimated for red sources (see e.g. Fabricius et al. 2021). In addition, we corrected the published G-band photometry of the sources with a six-parameter astrometric solution in the Gaia-EDR3 catalogue fainter than G  =  13 mag, using the formula presented in Riello et al. (2021), as explained in Appendix A of Gaia Collaboration (2021).

We modelled the field population using Gaussian mixture models (GMM)1 with a random sample of 106 sources from the input catalogue. We tested the field model with different numbers of GMM components (i.e. 120, 140, 160, 180, 200, 220, and 240) and took the model with 160 components that returns the minimum Bayesian information criterion (BIC, Schwarz 1978) value. The cluster was modelled using GMM in the space of proper motions and a principal curve in the photometric space, with a spread at any point given by a multivariate Gaussian to define the cluster sequence. Both the cluster locus in the space of proper motions and the principal curve in the photometric space are initialised based on an initial list of members. We use the sample of cluster members identified by ZKK2019 as initial list for the first iteration of our membership analysis. The method computes membership probabilities for all sources in the catalogue and classifies them into members and non-members based on a probability threshold, pin, which is predefined by the user. The cluster members that result from this process serve as input for the next iteration of the algorithm and this procedure is repeated until the list of cluster members remains fixed after successive iterations. Then we use the field and cluster models inferred by our method to generate synthetic data and define the optimum probability threshold popt for our membership analysis based on the performance of our classifier (see Olivares et al. 2019). Finally, we classify the sources in our catalogue as members (p ≥ popt) and non-members.

In Table 1, we compare the results of our membership analysis using different values for the user-defined probability threshold pin. The true positive rate (TPR, i.e. the fraction of cluster members recovered by the method) and contamination rate (CR, i.e. the fraction of field stars classified as cluster member by the method) given for each solution are derived from the synthetic data (based on the cluster and field models inferred in the previous steps) and, therefore, these represent only estimates of the indicators that remain useful in the absence of the true distributions. We see little variation in the number of cluster members among the solutions obtained from different pin thresholds. However, we note that high pin values (pin = 0.8 and pin = 0.9) lead to low values of the optimum probability threshold, popt, and, consequently, to the inclusion of stars with lower membership probabilities in our solution, while the least conservative threshold of pin = 0.5 yields the highest popt value. The solution obtained with pin = 0.6 offers the best compromise between the two probability thresholds. We therefore use this solution to define the list of XFOR cluster members derived in this study (see more comments below).

Table 1.

Comparison of results for the membership analysis using different probability thresholds.

We note that only one component of three binary systems in that sample was directly identified by our methodology, so we add the missing components of these systems, namely: Gaia EDR3 5053156704532702720 (032317–343547), Gaia EDR3 4854879573671981056 (032215–354719), and Gaia EDR3 4854540344270707200 (032513–370909) to the list. These sources have proper motions and parallax that are consistent with membership in the cluster (see Fig. 1) and they were rejected by our method due to their photometry. Gaia EDR3 5053156704532702720 (032317–343547) is the faintest cluster member and it is located in a region of the colour-magnitude diagram where our model is less efficient because it is trained on only a few sources in this colour regime (i.e. candidate members from the literature). On the other hand, we see that Gaia EDR3 4854879573671981056 (032215–354719) and Gaia EDR3 4854540344270707200 (032513–370909) are significantly above the empirical isochrone of the cluster (see Fig. 2), which results from a less reliable photometry. Indeed, we verified that all transits in the GRP band for these two sources are affected by blending effects (as indicated by the fields phot_rp_n_obs and phot_rp_n_blended_transits in the Gaia-EDR3 catalogue).

thumbnail Fig. 1.

Proper motions and parallaxes of the 164 cluster members identified in our membership analysis. The stars are colour-coded based on their membership probability that is scaled from zero to one. The three sources marked with green triangles denote the components of binary systems that we added manually to the sample, as explained in the text of Sect. 2.

We also searched for other astronomical objects (among the rejected sources) that exhibit astrometric features consistent with membership in the cluster but are automatically discarded by our method due to their relative position with respect to the cluster isochrone in the colour-magnitude diagram. We find one white dwarf candidate, namely, Gaia EDR3 4856170984438879232 (see Gentile Fusillo et al. 2019), that has proper motion (μα cos δ = 30.218 ± 1.099 mas yr−1, μδ = −5.781 ± 1.339 mas yr−1) and parallax (ϖ = 8.270 ± 0.979 mas yr−1) consistent with the other cluster members in the sample. However, based on the Bergeron et al. (1995) evolutionary sequences for white dwarfs2, we conclude that this source is much older (about 9 Gyr) than the XFOR cluster (see discussion in Sect. 4.4) and it is unlikely to be a cluster member.

Our final list of cluster members consists of 164 stars (see Table A.1). We also provide in Table A.2 the membership probability for all sources in the input catalogue used for the membership analysis derived from the different pin values investigated in our study. This offers subsequent studies the opportunity to select cluster members with criteria that are more specific to their scientific objectives than the ones used in this paper. Our new sample of cluster members includes stars in the magnitude range from G = 5.1 to G = 19.6 mag. The absolute magnitude of the stars range from about Gabs = −0.3 to Gabs = 14.7 mag. Figures 1 and 2 show the proper motions, parallaxes, and colour-magnitude diagram of the selected sources. The empirical isochrone of the cluster derived in our membership analysis is given in Table A.3. The observed correlation between μα cos δ and ϖ illustrated in Fig. 1 can be explained by projection effects. Cluster members located at different sky positions (but with similar space velocities) display proper motions and parallaxes that vary in proportion to each other.

thumbnail Fig. 2.

Colour-magnitude diagram of the cluster members selected in our membership analysis. The stars are colour-coded based on their membership probability that is scaled from zero to one. The black solid line shows the empirical isochrone of the cluster (see Table A.3). The arrow indicates the extinction vector of AV = 1 mag converted to the Gaia bands using the relative extinction values given by Wang & Chen (2019). The three sources marked with green triangles denote the components of binary systems that we added manually to the sample, as explained in the text of Sect. 2.

The XFOR cluster is located in a sky region that overlaps with the nearby THA and COL young stellar associations. We compiled the lists of known members given in the literature for THA and COL (Kraus et al. 2014; Gagné et al. 2018a,b; Gagné & Faherty 2018), cross-matched them with the list of stars identified in our membership analysis, and confirmed that our sample of XFOR cluster members is not contaminated by known members of these two stellar groups. The published lists of stars in THA and COL suggest that these stellar groups are closer to Earth than the XFOR cluster by about 50 pc. In a recent study, Gagné et al. (2021) proposed that many of the nearby moving groups in the Solar neighbourhood are more extended than previously thought and may be fragments of dissolving open clusters. In particular, the authors associate THA and COL with the IC 2602 and Platais 8 open clusters, respectively. Thus, if this scenario is correct, then future studies may reveal the existence of additional THA and COL members in the vicinity of the XFOR cluster.

Our membership analysis confirms 103 stars (out of 107 sources, see Fig. 3) that were previously identified as cluster members by ZKK20193. Only four stars (Gaia EDR3 4860808316464009856, Gaia EDR3 5046516311431900160, Gaia EDR3 4852758100345682816, and Gaia EDR3 5048052294816437120) from that study cannot be recovered by our methodology. One of them, namely, Gaia EDR3 4860808316464009856 (CD-36 1289), was classified as a “core-member” of the cluster by ZKK2019 based on its proper motion (μα cos δ = 34.6 ± 1.0 mas yr−1 and μδ = −1.8 ± 1.0 mas yr−1), given in the UCAC4 catalogue (Zacharias et al. 2013). Indeed, the Gaia-EDR3 proper motion (μα cos δ = 25.255 ± 0.054 mas yr−1, μδ = −1.366 ± 0.080 mas yr−1) is not consistent with membership in the cluster, which explains the reason it was rejected by our method. The same conclusion holds if we use the Gaia-DR2 measurements in the membership analysis. Thus, we discard this star from our study for consistency with the rest of our membership analysis where we use for all stars the same source of astrometric data (i.e. Gaia-EDR3) for both proper motion and parallax. The other three sources rejected in this study are classified as “potential members” by ZKK2019. Our membership analysis shows that they exhibit low membership probabilities (p < 0.13) that are below the adopted probability threshold (see Table 1).

thumbnail Fig. 3.

Venn diagram comparing the sample of cluster members identified in our membership analysis with the results obtained by ZKK2019. The area of each circle is proportional to the number of stars in the samples.

In this study, we identified 61 new cluster members, which represents an increase of about 60% in the census of XFOR stars with respect to the number of confirmed cluster members from the literature. In the following, we use this new sample of stars to re-visit some properties of the cluster.

3. Spectroscopic observations

We conducted a low-resolution spectroscopy study of our new sample of cluster members with the EFOSC2 spectrograph (Buzzoni et al. 1984) mounted at the ESO/NTT telescope, operating at La Silla (ESO, Chile). The observations were taken with the 1.0″ slit, grism Gr#16, and filter OG530, yielding a spectral coverage from 6015 Å to 10 320 Å and resolution of about 16 Å. The exposure times ranged from a few seconds up to 55 min depending on the magnitude of the target, so that the signal-to-noise ratio (S/N) is often better than 50 for most spectra collected in our campaign. We performed our observations in designated visitor mode (i.e. remote observations) during the nights of November 23 to November 30, 2020 (programme: 106.D-0429, PI: Galli). We observed 91 stars from our sample of cluster members identified in this study. Two of the new cluster members identified in this paper, namely Gaia EDR3 4755198025592857088 (025251–435842) and Gaia EDR3 4755198025592857216 (025251–435845), form a close binary system and we observed both components in our campaign.

The spectra were reduced with the EFOSC2 pipeline using the EsoRex command line tool, which performs a bias subtraction, flat-fielding, along with wavelength and flux calibration. We observed five spectrophotometric standards (Feige110, CD-34 241, LTT1020, LTT2415, and LTT3218) during the nights allocated to our programme and took the closest one (in time and airmass) to each science frame to flux calibrate the spectrum. The reduced spectra were then used to derive the spectral type of the stars as explained below.

We take advantage of the PyHammer v2.0 facility (Roulston et al. 2020) to perform the spectral classification of the stars in our sample. This code is based on the IDL spectral typing suite called ‘The Hammer’ developed by Covey et al. (2007) and recently re-written in Python programming language. PyHammer automatically assigns spectral types to the input spectra by measuring spectral indices of atomic and molecular lines and comparing them to templates of an empirical stellar library (see Kesseli et al. 2017). The input spectrum is compared to each template by computing the chi-squared value between the two spectra. The best matching spectral type is defined as the one that minimises the chi-squared value between the spectral indices from the template and input spectrum.

To investigate the robustness of our spectral classification, we implemented a chi-squared minimisation routine that uses the entire spectrum (instead of spectral indices). We ran this comparison between input spectrum and templates using the same spectral library of the PyHammer code. The spectral types derived from the two methods are consistent among themselves within one subclass for M stars and two subclasses for the other stars in our sample. We adopted these numbers as uncertainties of our spectral classification and report the spectral types derived directly with PyHammer as our final results.

We complemented the spectral classification of cluster members with additional measurements that are available in the literature, mostly for the stars of spectral types A and F of our sample. Table A.1 lists the spectral classification of the stars and the references from the literature. By combining the spectral types derived in this study from EFOSC2 observations and literature results, we were able to provide the spectral classification for 100 stars that we use in the forthcoming analysis (see Sect. 4). The spectral type of the stars ranges from B6.5 to M8.

4. Properties of the cluster

4.1. Radial velocities

The scarcity of radial velocities for the stars in our sample is the main limitation to access the 3D space of motion of the XFOR cluster. We proceed as follows to complement the Gaia-EDR3 astrometry with radial velocity data. Firstly, we searched the CDS databases for published radial velocities in the literature for the stars in our sample. We found radial velocity information for 24 stars that have been collected from Gontcharov (2006), Kordopatis et al. (2013), Shkolnik et al. (2017), and the Gaia-EDR3 catalogue itself4. Secondly, we searched for high-resolution spectra in public archives and derived the radial velocities ourselves (see below). We found spectroscopic data for four stars in our sample (one of them with multiple spectra). The results of our search in public archives are summarised in Table 2.

Table 2.

Properties of the data downloaded from public archives.

We use the reduced spectra downloaded from the archives to derive the radial velocity of the stars with available data. We use the iSpec facility (Blanco-Cuaresma et al. 2014) to cross-correlate the spectrum of each star with masks of different spectral types (A0, F0, G2, K0, K5, and M5). The iSpec routines implement the cross-correlation algorithm to generate the velocity profile of the star (see e.g. Baranne et al. 1996), and fit a Gaussian/Voigt to the cross-correlation function (CCF) to derive its radial velocity with the associated uncertainty. We inspected the CCFs individually and discarded the radial velocity measurements that result from a poor fit to the CCF that can arise, for instance, due to the low S/N of the collected spectra or a mismatch between the spectral type of the stars and the adopted mask. We compute the radial velocity for each star from the closest mask available in iSpec to its spectral type. The scatter in the radial velocities that results from using masks of different spectral types ranges from 0.1 to 0.3 km s−1 for the stars with available spectra in our sample. We add the radial velocity scatter obtained from different masks in quadrature to the formal uncertainty on the radial velocity given by iSpec for each star. This procedure is likely to overestimate the uncertainties derived in this study, but it accounts for the observed fluctuation on radial velocity results obtained with masks of different spectral types (which often do not precisely match the spectral type of the star).

We note that two stars for which we derived the radial velocities ourselves, namely Gaia EDR3 4860643905115917312 (032807–355444) and Gaia EDR3 5056560315790459776 (033121–303059), have known radial velocities (see comparison in Table 3). The published radial velocities for these two sources are less precise than the ones obtained in this study, so we use the results derived by us in the upcoming analysis. The median precision of the radial velocities derived in this paper from high-resolution spectra in public archives is 0.4 km s−1, while the radial velocities published in the literature (including those from the Gaia-EDR3 catalogue) have a median precision of 1.9 km s−1.

Table 3.

Comparison of the radial velocities derived in this paper with the literature results.

By combining the radial velocities from the literature with the ones derived in this paper we end up with a sample of 28 stars with available radial velocity information. Two sources in this sample, namely Gaia EDR3 5056578663891978880 (033036–302006) and Gaia EDR3 5058417459651237504 (031426–300700), have poor radial velocities (Vr = 30.5 ± 6.2 km s−1 and Vr = −4.6 ± 17.5 km s−1) with uncertainties that greatly exceed the observed velocity dispersion of a few km s−1 in young stellar groups. We have therefore discarded these measurements from our analysis. We also note the existence of discrepant radial velocity measurements in this sample as illustrated in Fig. 4. We applied the interquartile range (IQR) criterion to identify the outliers in the distribution of radial velocities by rejecting the measurements that lie 1.5 × IQR below the first quartile or above the third quartile. We discard the radial velocity of two stars based on this method. However, we do not remove these sources from our list of cluster members because their proper motions, parallaxes, and photometry are consistent with membership in the cluster (p ≥ 0.9). The discrepant radial velocity is more likely to be affected by other reasons (e.g. undetected binarity) and this requires further investigation with a greater store of data. In addition, we verified that three stars (with radial velocity data) are binary systems and we decided to also discard these measurements from the forthcoming analysis (see Sect. 4.2). Altogether, this reduces the sample of cluster members with radial velocity information to 21 stars. The radial velocities used in this paper are given in Table A.1.

thumbnail Fig. 4.

Kernel density estimate of the distribution of radial velocities using a kernel bandwidth of 2 km s−1 that roughly corresponds to the median precision of the radial velocities in our sample. The tick marks in the horizontal axis indicate the individual radial velocity measurements of the stars.

4.2. Distance and spatial velocity

We converted the proper motion and parallax of the stars in our sample into distance and 2D tangential velocities using Bayesian inference following the online tutorials available in the Gaia archive (Luri et al. 2018). We used the Kalkayotl package (Olivares et al. 2020) to investigate the dependency of our results on the choice of the prior for the distance. The Kalkayotl code implements two families of priors for the distance that are either based on statistical probability density distributions (hereafter, statistical priors) or astrophysical assumptions (hereafter, astrophysical priors). We performed our calculations with two statistical priors that implement the Uniform and Gaussian probability density distributions, and another two astrophysical priors that are based on the surface brightness profile of young star clusters in the Large Magellanic Clouds (Elson et al. 1987) and the profile distribution of globular clusters observed by King (1962). Our analysis confirms that our distance results are independent of our choice of priors for the distance thanks to the good precision of the Gaia-EDR3 parallaxes. Thus, the individual distances listed in Table A.1 refer to the solution obtained with the Uniform prior, which is the least informative (i.e. most simple) prior investigated in our analysis.

Then we computed the 3D position of the stars from their distances in a XYZ reference system that has its origin at the Sun where X points to the Galactic centre, Y points to the direction of Galactic rotation, and Z points to the Galactic North pole. We converted the 2D tangential velocities and radial velocities of the stars into the UVW components of the Galactic velocity in the same reference system as before using the transformation outlined by Johnson & Soderblom (1987). The individual distances and spatial velocity of the stars are given in Table A.1.

The distance to the XFOR cluster determined from Bayesian inference in this study is 108.4 ± 0.3 pc and the distances of individual stars in the cluster range from 91.1 ± 0.2 pc to 125.7 ± 0.6 pc. To gain confidence in our results, we recomputed the distances (as previously explained) after correcting the Gaia-EDR3 parallaxes for the colour-dependent bias (i.e. zero-point correction) following the recipe outlined by Lindegren et al. (2021a), and we confirmed that our results are not significantly affected by this correction. The proposed functions to correct for the parallax bias are only indicative (as explained by the authors themselves in their study) and they do not bring about significant changes in our case due to the close proximity of the cluster. We therefore report our distance results based on the original parallaxes given in the Gaia-EDR3 catalogue (i.e. without applying any correction).

We derived the cluster space motion from the procedure described in Sect. 7.2 of Perryman et al. (1998) to account for the existence of correlated errors in the velocity components that result from the transformation of proper motions, parallaxes, and radial velocities into 3D velocities. The mean space motion is obtained by averaging the measured velocities and the uncertainty is computed from the mean of the covariance matrices of individual stars. This yields a mean space motion of U = −12.8 ± 0.6 km s−1, V = −21.6 ± 0.8 km s−1, W = −4.9 ± 1.4 km s−1, and a standard deviation of σU = 1.6 km s−1, σV = 1.7 km s−1, σW = 2.6 km s−1. The median uncertainties in the UVW velocity components are 0.3, 0.5, and 1.1 km s−1, implying that the internal velocity dispersion of the cluster is resolved.

The space motion of the XFOR cluster can now be compared with the space motion of the COL and THA associations. The relative motion of XFOR with respect to the space motions of COL and THA given by Gagné et al. (2018a) in the sense XFOR ‘minus’ COL or THA are (ΔU, ΔV, ΔW) = (−0.9, −0.3, 0.8)±(1.2, 1.5, 1.6) km s−1 and (ΔU, ΔV, ΔW) = (−3.0, −0.7, −3.9)±(1.1, 1.1, 1.6) km s−1, respectively. This comparison shows that the space motion of XFOR is consistent with COL and THA within 3σ of the reported uncertainties in the spatial velocities when the correlations are taken into account. We therefore confirm the previous finding of ZKK2019 based on the space motion of only four core members that the XFOR cluster is comoving with the COL and THA young stellar associations.

Finally, we performed a simplified traceback analysis by integrating the orbit of the XFOR cluster back in time and comparing its location in the past with the other young stellar groups of the Solar neighbourhood. Here, we follow a similar procedure as described in Miret-Roig et al. (2018) which consists of integrating the centroid of each stellar group. Thus, we use the UVW velocity derived in this section and the mean XYZ position of the stars to define the space motion and present-day location of the XFOR cluster. For the other young stellar groups, we use the values given in Table 9 of Gagné et al. (2018a). Our calculations use the heliocentric curvilinear coordinate system (ξ′,η′,ζ′) defined by Asiain et al. (1999) and the same 3D axisymmetric Milky Way potential as in Miret-Roig et al. (2020) to integrate the equations of motion. This model is based on the Allen & Santillan (1991) potential which consists of a spherical central bulge, a disc, and a massive spherical halo, but with updated parameters taken from Table 1 of Irrgang et al. (2013). As illustrated in Fig. 5, the minimum distance of the XFOR cluster to the THA and COL stellar groups occurred at t = −10.6 Myr (∼24 pc) and t = −19.1 Myr (∼18 pc), respectively. Interestingly, the closest approach of the XFOR cluster occurred with the Carina (CAR) young stellar group at t = −18.6 Myr with a separation of only ∼7 pc. The relative space motion of the XFOR cluster and CAR (see Gagné et al. 2018a) is (ΔU, ΔV, ΔW) = (−2.1, 0.3, 0.6)±(0.9, 1.3, 1.7) km s−1 indicating that these two young stellar groups are also comoving. Altogether, this suggests that XFOR, THA, COL, and CAR may have a common origin and that they could be the remnant of a much larger association.

thumbnail Fig. 5.

Orbits of the young stellar groups integrated back in time in a curvilinear coordinate system using a 3D axisymmetric galactic potential. The orbits are colour-coded based on the integration time used in our calculations. The thicker lines roughly indicate the assumed age (∼30−40 Myr) of these young stellar groups given in the literature (see also Sect. 4.4).

4.3. Internal motions

The 3D velocity vectors of the cluster members with useful radial velocity measurements (see Sect. 4.1) are mostly parallel to each other. The non-trivial question that arises at this stage is whether the observed velocities result from random or organised motions within the cluster and their quantitative importance. In this section, we investigate the existence of potential expansion and rotation effects in the XFOR cluster following the procedure applied by Rivera et al. (2015) and Galli et al. (2019) in the Taurus star-forming region.

First, we computed the unit position vector r ̂ = r / | r | $ {\hat{\boldsymbol{r}}_{\star}}={\boldsymbol{r}}_{\star}/|{\boldsymbol{r}}_{\star}| $ of each star that represents its distance to the centre of the cluster. Then we computed the relative velocity δv of each star with respect to the mean velocity of the cluster (see Sect. 4.2). The dot product ( r ̂ · δ v $ {\hat{\boldsymbol{r}}_{\star}} \cdot \delta {\boldsymbol{v}_{\star}} $) is large and positive (negative) if the cluster is undergoing expansion (contraction). Similarly, the cross product ( r ̂ × δ v $ {\hat{\boldsymbol{r}}_{\star}} \times \delta {\boldsymbol{v}_{\star}} $), which is related to angular momentum, is large (small) when the cluster has significant (negligible) rotation. It is important to mention that the dot and cross products are not exactly the expansion and rotation velocities of the cluster, but they can be used as a proxy for these quantities to access the existence (or absence) of these motions within the cluster.

We computed r ̂ · δ v $ {\hat{\boldsymbol{r}}_{\star}} \cdot \delta {\boldsymbol{v}_{\star}} $ and r ̂ × δ v $ {\hat{\boldsymbol{r}}_{\star}} \times \delta {\boldsymbol{v}_{\star}} $ for each star in our sample individually and then we took the average and standard error of the mean of the resulting values to characterise the cluster as a whole. In doing so, we find that r ̂ · δ v = 0.5 ± 0.4 $ {\hat{\boldsymbol{r}}_{\star}} \cdot \delta {\boldsymbol{v}_{\star}} =0.5\pm 0.4 $ km s−1 and r ̂ × δ v = ( 0.1 , 0.1 , 0.4 ) ± ( 0.3 , 0.4 , 0.3 ) $ {\hat{\boldsymbol{r}}_{\star}} \times \delta {\boldsymbol{v}_{\star}}=(0.1,-0.1,0.4)\pm(0.3,0.4,0.3) $ km s−1, showing that both quantities are consistent with zero within 1σ of the reported uncertainties. Thus, we conclude that the effects of a potential expansion and rotation of the cluster are negligible and that the observed velocity dispersion is most likely to be dominated by random motions within the cluster. However, we should note that we reached the present conclusion based on a small fraction of cluster members with complete data for this analysis (3D positions and 3D velocities); thus, it is clear that more study is clearly warranted when the radial velocity of a more significant number of cluster members becomes available in the future.

4.4. Isochronal ages

The first step of our analysis to revisit the age of the cluster consists of constructing the Hertzsprung-Russell (HR) diagram of the cluster. Thus, we fit the spectral energy distribution (SED) of the stars with the Virtual Observatory SED Analyzer (VOSA, Bayo et al. 2008) to derive the effective temperatures and bolometric luminosities. We used the BT-Settl grid of models (Allard et al. 2012) to fit the SEDs and leave the extinction AV as a free parameter in the fit varying in the range of 0 to 10 mag. We provide the input photometry to the VOSA facility to avoid erroneous cross-matches with the system interface. We cross-matched our list of cluster members with the 2MASS (Cutri et al. 2003), AllWISE (Wright et al. 2010), Pan-STARRS (Chambers et al. 2016), APASS (Henden et al. 2015), and Tycho2 (Høg et al. 2000) catalogues to retrieve the photometric data for our stars using the precomputed cross-matched tables available in the Gaia archive. In doing so, we derived effective temperatures and bolometric luminosities for 154 stars in our sample with available photometry for the SED analysis with VOSA.

We compared the effective temperatures obtained from the SED fit with the ones inferred from the spectral type of the sample of 100 stars with available spectral classification (see Sect. 3). We used the tabulated values from Pecaut & Mamajek (2013) to convert the spectral type of the stars to effective temperatures. The mean difference between the effective temperatures in the two data sets is 136 ± 19 K, which is consistent with the uncertainty obtained in our spectral classification. Similarly, we compared the bolometric luminosities derived from VOSA with the results we computed, using the J-band photometry from the 2MASS catalogue and the bolometric corrections given by Pecaut & Mamajek (2013) for a sample of 96 stars with available spectral classification and 2MASS photometry. The results of this comparison confirm the robustness of our results obtained with different methods. In the following, we restrict our discussion to the results obtained directly from the VOSA service because this method allows us to construct the HR diagram with a much more significant number of cluster members.

Figure 6 presents the resulting HR diagram. We compare the position of the stars in the HR diagram with four grids of pre-main sequence star models: Baraffe et al. (2015, BHAC15), Siess et al. (2000, SDF00), PARSEC v1.2S (Bressan et al. 2012), and MIST (Choi et al. 2016), with stellar rotation (ν = 0.4νcrit). Our sample covers the mass range from about 0.02 to 4 M, but we show in Fig. 6 only the mass domain that includes most cluster members and where the models overlap for a better comparison. The lower mass limit is consistent with a brown dwarf of spectral type ∼ M8 and magnitude of G ∼ 18 mag. We note from the HR diagram that many cluster members appear above the 10 Myr isochrone given by the BHAC15, SDF00 and MIST models. This suggests that the cluster can be somewhat younger than previously thought. However, the fraction of young disc-bearing stars in the cluster is rather low to support the hypothesis of the cluster being only a few Myr old (see Sect. 4.5) and some of these sources are probably binary or high-order multiple systems that will require further investigation in future studies. The observed spread in the HR diagram for the least massive stars in the sample (< 0.3 M) indeed translates into a broad range of ages, but the stars with higher masses give a much tighter constrain and clearly discard an age younger than 10 Myr for the cluster.

thumbnail Fig. 6.

HR diagram of the XFOR cluster compared to the grid of isochrones and tracks computed from different models. The solid red lines and black dashed lines indicate the isochrones and tracks for each model, with the age (in Myr) and mass (in M) values indicated in the panels.

From a visual inspection in the colour-magnitude diagram (see Fig. 7), it is apparent that the empirical isochrone of the cluster – derived in our analysis and converted to absolute magnitudes at the distance of the cluster (108.4 ± 0.3 pc, see Sect. 4.2) – is more consistent with an age of about 30 Myr when compared to the isochrones of the BHAC15 and PARSEC models. This age is consistent with the median age of 35 Myr inferred from the PARSEC models based on interpolating the stars in the HR diagram between the isochrones. It is also in good agreement with the age estimate given by Mamajek (2016) based on the X-ray emission of a few cluster members. As illustrated in Fig. 8 the age of 30 Myr is also still consistent with the assumed age for the COL (42 4 + 6 $ ^{+6}_{-4} $ Myr, Bell et al. 2015), THA (45 ± 4 Myr, Bell et al. 2015), and CAR (45 7 + 11 $ ^{+11}_{-7} $ Myr, Bell et al. 2015) stellar groups, supporting the idea that they are all coeval (see also ZKK2019).

thumbnail Fig. 7.

Absolute colour-magnitude diagram of the cluster members identified in our membership analysis (see Sect. 2). The lines indicate the empirical isochrone of the cluster and the 30 Myr isochrones computed from the BHAC15 and PARSEC models. The three stars located above the cluster sequence are Gaia EDR3 5056560315790459776 (033121–303059), Gaia EDR3 4854540344270707200 (032513–370909), and Gaia EDR3 4854879573671981056 (032215–354719).

thumbnail Fig. 8.

Absolute colour-magnitude diagram derived from Gaia-EDR3 photometry. The different colours and symbols indicate the XFOR cluster members identified in this study and known members of the COL, THA, and CAR stellar groups published in the literature (Gagné et al. 2018a,b; Gagné & Faherty 2018). The solid lines mark the 10 Myr and 40 Myr isochrones from the PARSEC models.

At this stage it seems premature to assign a definitive and unambiguous age for the cluster given the important age variation that we observe from the HR diagram depending on our choice of evolutionary model (see Fig. 6). In the remainder of this paper, we assume an age of 30 Myr for the cluster, but we believe that more study using alternative methods, such as the traceback (see Miret-Roig et al. 2020) and lithium depletion boundary (see Binks & Jeffries 2014) methods, is required to provide a more robust result for the age.

4.5. Spatial distribution

Figure 9 shows the spatial distribution of the stars from our sample of cluster members with available spectral classification (see Sect. 3). Using the spectral type of the stars as a proxy for the stellar mass, we can confirm from our observations that the cluster is mass segregated with the more massive stars located preferentially in the central portions of the cluster and the M dwarfs spread across multiple locations (see also ZKK2019). For example, we note that nine stars among the 15 cluster members with spectral types B, A, and F in our sample are located within 10 pc of the 3D median position of the cluster (X, Y, Z) = (−32, −49, −87) pc. The three most dispersed early-type cluster members are Gaia EDR3 4843613397714473088 (035629–385744), Gaia EDR3 4949081572411575552 (025352–401139), and Gaia EDR3 4949158198924394496 (025048–395556) which are located farther than 15 pc from the cluster centre. The scatter of XYZ positions is (σX, σY, σZ) = (5.6, 3.8, 3.0) pc for the B, A, and F stars in our sample (excluding the three most dispersed stars beyond 15 pc of the centre), and (σX, σY, σZ) = (10.0, 5.3, 4.1) pc for the remaining cluster members of spectral types G, K, and M.

thumbnail Fig. 9.

3D spatial distribution of the stars in the cluster and the effects of mass segregation. The size and colour of the symbols are scaled with the spectral type of the stars derived in this study (see Sect. 3) in the sense that small red symbols denote the M dwarfs in our sample.

We now compare the spatial distribution of the stars with different SED subclasses in our sample. We use the classification scheme proposed by Koenig & Leisawitz (2014) to identify Class I and Class II sources based on infrared colours from the AllWISE catalogue. We classified 142 stars in our sample of cluster members with available photometric data for this analysis (see Fig. 10) and we found only eight Class II stars (i.e. 5% of the sample) which are listed in Table A.1. These sources are late-type stars as confirmed from our observations and results derived from the SED fit. Five of them are listed in Tables 2 and 3 of ZKK2019. The remaining three sources, namely Gaia EDR3 4755745925980823680 (025902–423220), Gaia EDR3 5060663418310359424 (032443–273323), and Gaia EDR3 5081640790199917696 (033852–264615), are new members identified in this study. The remaining cluster members are likely to be Class III stars, given the young ages that we observe from the HR diagram and they are the dominant population in our sample (as expected for a 30 Myr old cluster). As illustrated in Fig. 11, the Class II stars are mostly concentrated in the central portions of the cluster.

thumbnail Fig. 10.

Colour-colour diagram of cluster members using infrared photometry from the AllWISE catalogue. This is one of the diagrams used in the classification scheme developed by Koenig & Leisawitz (2014) to distinguish the SED subclasses based on infrared excess emission.

thumbnail Fig. 11.

3D spatial distribution of the cluster members identified in our membership analysis with SED classification derived in this study (see Sect. 4.5). The different colours and symbols denote the SED subclasses.

4.6. Initial mass function

In this section, we use our new sample of XFOR stars to discuss the initial mass function of the cluster. The stellar masses inferred from evolutionary models can vary significantly depending on the grid of models that is used in the analysis and some models do not cover the entire mass domain of our sample (see Fig. 6). So, to overcome this difficulty, we used the distribution of spectral types in the sample derived from our observations (see Sect. 3) as a proxy for the initial mass function.

We estimate the completeness limits of our sample of XFOR stars from the completeness limits of the Gaia-EDR3 catalogue, which is the only source of data used in the membership analysis. The Gaia-EDR3 catalogue has a high completeness at G = 19 mag (see Sect. 2.2 of Fabricius et al. 2021) and provides proper motion and parallax (i.e. five and six-parameter solutions) for about 98.3% of the detected sources up to this magnitude (see Tables 4 and 5 of Lindegren et al. 2021b). We therefore consider our sample to be complete up to G = 19 mag. We used the 3D dust maps produced by Lallement et al. (2019) to estimate the extinction at the position of the stars in our sample. We verified that the maximum observed interstellar extinction for our sample (within the 1σ uncertainty of the reddening values delivered by the 3D maps) is AG = 0.04 mag. Thus, the completeness limit corrected for the maximum observed extinction at the distance of the cluster (d = 108.4 ± 0.3 pc, see Sect. 4.2) yields the absolute magnitude of Gabs = 13.79 mag. Now, if we assume that the cluster has an age of 30 Myr and we use the BHAC15 grid models to convert the absolute magnitude into mass, we may conclude that our sample is complete down to 0.04 M.

In Fig. 12, we illustrate the distribution of spectral types of the XFOR cluster and compare it to the samples of stars in IC 348 (Muench et al. 2007), Taurus (Esplin & Luhman 2019), and Upper Scorpius (Luhman et al. 2018). We show that both the XFOR cluster and the Upper Scorpius associations exhibit a surplus of early-type stars as compared to Taurus and IC 348. However, it is apparent that the distribution of the late-type stars in the XFOR cluster down to the completeness limit of our sample resembles the other young stellar groups. The similar distribution of spectral types that we observe among these young stellar groups suggests that the initial mass function most likely shows little variation. The Taurus and Upper Scorpius samples of stars include a more significant number of substellar objects, which are still missing in our sample of XFOR cluster members. These samples have been obtained from ancillary observations in the optical and infrared, in addition to the Gaia data. We therefore argue that future studies in the XFOR cluster using additional material to complement the Gaia catalogue will facilitate the detection of additional cluster members and potentially reveal a population of substellar objects in the cluster, extending the initial mass function down to the planetary-mass regime.

thumbnail Fig. 12.

Distribution of spectral types of the samples of stars in the XFOR cluster (this paper), IC 348 (Muench et al. 2007), Taurus (Esplin & Luhman 2019), and Upper Scorpius (Luhman et al. 2018). The dashed line indicates the completeness limit of the XFOR sample.

5. Conclusions

In this study, we performed a new membership analysis of the XFOR cluster based on the astrometric and photometric data of the Gaia-EDR3 catalogue. We applied a probabilistic method to infer membership probabilities of the stars over a field of more than 1000 deg2 and selected 164 cluster members. We confirmed most of the cluster members that have been previously identified in the literature based on Gaia-DR2 data and discovered new ones increasing the list of cluster members by about 60%. We conducted low-resolution spectroscopy of the cluster members, derive the spectral type of the stars and use our new sample of members to revisit some properties of the cluster.

Our results obtained in this study show that the XFOR cluster is located at 108.4 ± 0.3 pc and it is comoving with the THA, COL, and CAR young stellar associations. A dynamical traceback reveals that the cluster was significantly closer to these young stellar groups in the past, which suggests a common origin. We do not detect important expansion or rotation effects in the cluster and we argue that the observed velocity dispersion is most likely associated with random motions. The spatial distribution of the stars shows that the cluster is mass segregated and that the few disc-bearing stars in our sample are mostly located in the central portions of the cluster. The empirical isochrone of the cluster that we derive from our membership analysis is consistent with the age of 30 Myr. However, the observed scatter in the HR diagram, in particular for the least massive cluster members, along with the age differences among the various models, call for a revision using other age-dating techniques. We estimate our sample of cluster members to be complete down to about 0.04 M and argue that the initial mass function of the cluster (inferred from the distribution of the spectral types of the stars) shows little variation compared to other young stellar groups.

In this paper, we use the most recent data release delivered by the Gaia space mission combined with archival spectra and our own observations to deliver a complete picture of the stellar population and properties of the cluster. We believe that future studies using ancillary material to complement the Gaia catalogue have the potential to detect more substellar cluster members as done in other young stellar groups. The main limitation to the present study was the lack of high-resolution spectra. We therefore encourage astronomers to conduct high-resolution spectroscopy of the cluster that will allow us to investigate the kinematics of the stars in more detail, refine the age of the cluster, and prepare for exoplanet programmes. This will make the XFOR cluster a cornerstone for many studies related to star and planet formation at this young stage of stellar evolution.


1

A Gaussian Mixture Model is a model that describes the probability distribution of data points as a linear combination of a finite number of Gaussian distributions.

2

The evolutionary cooling sequences for white dwarfs used in this study were downloaded from http://www.astro.umontreal.ca/bergeron/CoolingModels.

3

ZKK2019 identified 107 candidate members with Gaia-DR2 data, but one of them, Gaia EDR3 4860808041586101632 (032559–355730), is not resolved by the Gaia satellite, so they refer to 108 cluster members throughout that study. Here, we count this system as one single entry in our tables.

4

The Gaia-EDR3 catalogue does not contain new radial velocity measurements. The radial velocities of Gaia-DR2 have been added to the Gaia-EDR3 catalogue for user convenience.

Acknowledgments

This research has received funding from the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation programme (grant agreement No 682903, P.I. H. Bouy), and from the French State in the framework of the “Investments for the future” Program, IdEx Bordeaux, reference ANR-10-IDEX-03-02. This research has made use of the SIMBAD database, operated at CDS, Strasbourg, France. This work has made use of data from the European Space Agency (ESA) mission Gaia (https://www.cosmos.esa.int/gaia), processed by the Gaia Data Processing and Analysis Consortium (DPAC, https://www.cosmos.esa.int/web/gaia/dpac/consortium). Funding for the DPAC has been provided by national institutions, in particular the institutions participating in the Gaia Multilateral Agreement. This publication makes use of data products from the Wide-field Infrared Survey Explorer, which is a joint project of the University of California, Los Angeles, and the Jet Propulsion Laboratory/California Institute of Technology, funded by the National Aeronautics and Space Administration. This publication makes use of VOSA, developed under the Spanish Virtual Observatory project supported by the Spanish MINECO through grant AyA2017-84089. VOSA has been partially updated by using funding from the European Union’s Horizon 2020 Research and Innovation Programme, under Grant Agreement no 776403 (EXOPLANETS-A).

References

  1. Allard, F., Homeier, D., & Freytag, B. 2012, Phil. Trans. R. Soc. London Ser. A, 370, 2765 [Google Scholar]
  2. Allen, C., & Santillan, A. 1991, Rev. Mex. Astron. Astrofis., 22, 255 [NASA ADS] [Google Scholar]
  3. Asiain, R., Figueras, F., & Torra, J. 1999, A&A, 350, 434 [NASA ADS] [Google Scholar]
  4. Baraffe, I., Homeier, D., Allard, F., & Chabrier, G. 2015, A&A, 577, A42 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  5. Baranne, A., Queloz, D., Mayor, M., et al. 1996, A&AS, 119, 373 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  6. Bayo, A., Rodrigo, C., Barrado Y Navascués, D., et al. 2008, A&A, 492, 277 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  7. Bell, C. P. M., Mamajek, E. E., & Naylor, T. 2015, MNRAS, 454, 593 [Google Scholar]
  8. Bergeron, P., Wesemael, F., & Beauchamp, A. 1995, PASP, 107, 1047 [Google Scholar]
  9. Binks, A. S., & Jeffries, R. D. 2014, MNRAS, 438, L11 [NASA ADS] [CrossRef] [Google Scholar]
  10. Blanco-Cuaresma, S., Soubiran, C., Heiter, U., & Jofré, P. 2014, A&A, 569, A111 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  11. Bouy, H., Bertin, E., Moraux, E., et al. 2013, A&A, 554, A101 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  12. Bressan, A., Marigo, P., Girardi, L., et al. 2012, MNRAS, 427, 127 [NASA ADS] [CrossRef] [Google Scholar]
  13. Buzzoni, B., Delabre, B., Dekker, H., et al. 1984, The Messenger, 38, 9 [NASA ADS] [Google Scholar]
  14. Cantat-Gaudin, T., Jordi, C., Vallenari, A., et al. 2018, A&A, 618, A93 [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. Chambers, K. C., Magnier, E. A., Metcalfe, N., et al. 2016, ArXiv e-prints [arXiv:1612.05560] [Google Scholar]
  17. Choi, J., Dotter, A., Conroy, C., et al. 2016, ApJ, 823, 102 [Google Scholar]
  18. Covey, K. R., Ivezić, Ž., Schlegel, D., et al. 2007, AJ, 134, 2398 [NASA ADS] [CrossRef] [Google Scholar]
  19. Cutri, R. M., Skrutskie, M. F., van Dyk, S., et al. 2003, VizieR Online Data Catalog: II/246 [Google Scholar]
  20. Dekker, H., D’Odorico, S., Kaufer, A., Delabre, B., & Kotzlowski, H. 2000, in Optical and IR Telescope Instrumentation and Detectors, eds. M. Iye, & A. F. Moorwood, SPIE Conf. Ser., 4008, 534 [Google Scholar]
  21. 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]
  22. Elson, R. A. W., Fall, S. M., & Freeman, K. C. 1987, ApJ, 323, 54 [Google Scholar]
  23. Esplin, T. L., & Luhman, K. L. 2019, AJ, 158, 54 [Google Scholar]
  24. Fabricius, C., Luri, X., Arenou, F., et al. 2021, A&A, 649, A5 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  25. Gagné, J., & Faherty, J. K. 2018, ApJ, 862, 138 [Google Scholar]
  26. Gagné, J., Mamajek, E. E., Malo, L., et al. 2018a, ApJ, 856, 23 [Google Scholar]
  27. Gagné, J., Roy-Loubier, O., Faherty, J. K., Doyon, R., & Malo, L. 2018b, ApJ, 860, 43 [Google Scholar]
  28. Gagné, J., Faherty, J. K., Moranta, L., & Popinchalk, M. 2021, ApJ, 915, L29 [CrossRef] [Google Scholar]
  29. Gaia Collaboration (Brown, A. G. A., et al.) 2016, A&A, 595, A2 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  30. Gaia Collaboration (Brown, A. G. A., et al.) 2018, A&A, 616, A1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  31. Gaia Collaboration (Brown, A. G. A., et al.) 2021, A&A, 649, A1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  32. Galli, P. A. B., Loinard, L., Bouy, H., et al. 2019, A&A, 630, A137 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  33. Galli, P. A. B., Bouy, H., Olivares, J., et al. 2020a, A&A, 634, A98 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  34. Galli, P. A. B., Bouy, H., Olivares, J., et al. 2020b, A&A, 643, A148 [EDP Sciences] [Google Scholar]
  35. Galli, P. A. B., Bouy, H., Olivares, J., et al. 2021, A&A, 646, A46 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  36. Gentile Fusillo, N. P., Tremblay, P.-E., Gänsicke, B. T., et al. 2019, MNRAS, 482, 4570 [Google Scholar]
  37. Gontcharov, G. A. 2006, Astron. Lett., 32, 759 [Google Scholar]
  38. Henden, A. A., Levine, S., Terrell, D., & Welch, D. L. 2015, Am. Astron. Soc. Meet. Abstr., 225, 336.16 [Google Scholar]
  39. Høg, E., Fabricius, C., Makarov, V. V., et al. 2000, A&A, 355, L27 [Google Scholar]
  40. Houk, N. 1982, Michigan Catalogue of Two-dimensional Spectral Types for the HD stars, 3 [Google Scholar]
  41. Irrgang, A., Wilcox, B., Tucker, E., & Schiefelbein, L. 2013, A&A, 549, A137 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  42. Jackson, J., & Stoy, R. H. 1955, Ann. Cape Obs., 18 [Google Scholar]
  43. Johnson, D. R. H., & Soderblom, D. R. 1987, AJ, 93, 864 [Google Scholar]
  44. Kaufer, A., Stahl, O., Tubbesing, S., et al. 1999, The Messenger, 95, 8 [Google Scholar]
  45. Kesseli, A. Y., West, A. A., Veyette, M., et al. 2017, ApJS, 230, 16 [NASA ADS] [CrossRef] [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. King, I. 1962, AJ, 67, 471 [Google Scholar]
  48. Koenig, X. P., & Leisawitz, D. T. 2014, ApJ, 791, 131 [Google Scholar]
  49. Kordopatis, G., Gilmore, G., Steinmetz, M., et al. 2013, AJ, 146, 134 [Google Scholar]
  50. Kraus, A. L., Shkolnik, E. L., Allers, K. N., & Liu, M. C. 2014, AJ, 147, 146 [Google Scholar]
  51. Lallement, R., Babusiaux, C., Vergely, J. L., et al. 2019, A&A, 625, A135 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  52. Lindegren, L., Bastian, U., Biermann, M., et al. 2021a, A&A, 649, A4 [EDP Sciences] [Google Scholar]
  53. Lindegren, L., Klioner, S. A., Hernández, J., et al. 2021b, A&A, 649, A2 [EDP Sciences] [Google Scholar]
  54. Luhman, K. L., Herrmann, K. A., Mamajek, E. E., Esplin, T. L., & Pecaut, M. J. 2018, AJ, 156, 76 [Google Scholar]
  55. Luri, X., Brown, A. G. A., Sarro, L. M., et al. 2018, A&A, 616, A9 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  56. Mamajek, E. E. 2016, in Young Stars& Planets Near the Sun, eds. J. H. Kastner, B. Stelzer, & S. A. Metchev, 314, 21 [NASA ADS] [Google Scholar]
  57. Meingast, S., Alves, J., & Fürnkranz, V. 2019, A&A, 622, L13 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  58. Miret-Roig, N., Antoja, T., Romero-Gómez, M., & Figueras, F. 2018, A&A, 615, A51 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  59. Miret-Roig, N., Bouy, H., Olivares, J., et al. 2019, A&A, 631, A57 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  60. Miret-Roig, N., Galli, P. A. B., Brandner, W., et al. 2020, A&A, 642, A179 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  61. Muench, A. A., Lada, C. J., Luhman, K. L., Muzerolle, J., & Young, E. 2007, AJ, 134, 411 [Google Scholar]
  62. Olivares, J., Bouy, H., Sarro, L. M., et al. 2019, A&A, 625, A115 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  63. Olivares, J., Sarro, L. M., Bouy, H., et al. 2020, A&A, 644, A7 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  64. Pecaut, M. J., & Mamajek, E. E. 2013, ApJS, 208, 9 [Google Scholar]
  65. Perryman, M. A. C., Brown, A. G. A., Lebreton, Y., et al. 1998, A&A, 331, 81 [NASA ADS] [Google Scholar]
  66. Riello, M., De Angeli, F., Evans, D. W., et al. 2021, A&A, 649, A3 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  67. Rivera, J. L., Loinard, L., Dzib, S. A., et al. 2015, ApJ, 807, 119 [NASA ADS] [CrossRef] [Google Scholar]
  68. Roulston, B. R., Green, P. J., & Kesseli, A. Y. 2020, ApJS, 249, 34 [NASA ADS] [CrossRef] [Google Scholar]
  69. Sarro, L. M., Bouy, H., Berihuete, A., et al. 2014, A&A, 563, A45 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  70. Schwarz, G. 1978, Ann. Stat., 6, 461 [Google Scholar]
  71. Shkolnik, E. L., Allers, K. N., Kraus, A. L., Liu, M. C., & Flagg, L. 2017, AJ, 154, 69 [Google Scholar]
  72. Siess, L., Dufour, E., & Forestini, M. 2000, A&A, 358, 593 [Google Scholar]
  73. Sim, G., Lee, S. H., Ann, H. B., & Kim, S. 2019, J. Korean Astron. Soc., 52, 145 [NASA ADS] [Google Scholar]
  74. Torres, C. A. O., Quast, G. R., da Silva, L., et al. 2006, A&A, 460, 695 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  75. Wang, S., & Chen, X. 2019, ApJ, 877, 116 [Google Scholar]
  76. Wright, E. L., Eisenhardt, P. R. M., Mainzer, A. K., et al. 2010, AJ, 140, 1868 [Google Scholar]
  77. Yen, S. X., Reffert, S., Schilbach, E., et al. 2018, A&A, 615, A12 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  78. Zacharias, N., Finch, C. T., Girard, T. M., et al. 2013, AJ, 145, 44 [Google Scholar]
  79. Zuckerman, B., Klein, B., & Kastner, J. 2019, ApJ, 887, 87 [NASA ADS] [CrossRef] [Google Scholar]

Appendix A: Tables

Table A.1.

Properties of the 164 cluster members selected from our membership analysis. This table will be available in its entirety in machine-readable form at the CDS.

Table A.2.

Membership probability of all sources in the field investigated in our membership analysis using different pin values. This table will be available in its entirety in machine-readable form at the CDS.

Table A.3.

Empirical isochrone inferred from our membership analysis. This table will be available in its entirety in machine-readable form at the CDS.

Table A.4.

Physical properties of the cluster member derived from the SED fit with VOSA. This table will be available in its entirety in machine-readable form at the CDS.

All Tables

Table 1.

Comparison of results for the membership analysis using different probability thresholds.

Table 2.

Properties of the data downloaded from public archives.

Table 3.

Comparison of the radial velocities derived in this paper with the literature results.

Table A.1.

Properties of the 164 cluster members selected from our membership analysis. This table will be available in its entirety in machine-readable form at the CDS.

Table A.2.

Membership probability of all sources in the field investigated in our membership analysis using different pin values. This table will be available in its entirety in machine-readable form at the CDS.

Table A.3.

Empirical isochrone inferred from our membership analysis. This table will be available in its entirety in machine-readable form at the CDS.

Table A.4.

Physical properties of the cluster member derived from the SED fit with VOSA. This table will be available in its entirety in machine-readable form at the CDS.

All Figures

thumbnail Fig. 1.

Proper motions and parallaxes of the 164 cluster members identified in our membership analysis. The stars are colour-coded based on their membership probability that is scaled from zero to one. The three sources marked with green triangles denote the components of binary systems that we added manually to the sample, as explained in the text of Sect. 2.

In the text
thumbnail Fig. 2.

Colour-magnitude diagram of the cluster members selected in our membership analysis. The stars are colour-coded based on their membership probability that is scaled from zero to one. The black solid line shows the empirical isochrone of the cluster (see Table A.3). The arrow indicates the extinction vector of AV = 1 mag converted to the Gaia bands using the relative extinction values given by Wang & Chen (2019). The three sources marked with green triangles denote the components of binary systems that we added manually to the sample, as explained in the text of Sect. 2.

In the text
thumbnail Fig. 3.

Venn diagram comparing the sample of cluster members identified in our membership analysis with the results obtained by ZKK2019. The area of each circle is proportional to the number of stars in the samples.

In the text
thumbnail Fig. 4.

Kernel density estimate of the distribution of radial velocities using a kernel bandwidth of 2 km s−1 that roughly corresponds to the median precision of the radial velocities in our sample. The tick marks in the horizontal axis indicate the individual radial velocity measurements of the stars.

In the text
thumbnail Fig. 5.

Orbits of the young stellar groups integrated back in time in a curvilinear coordinate system using a 3D axisymmetric galactic potential. The orbits are colour-coded based on the integration time used in our calculations. The thicker lines roughly indicate the assumed age (∼30−40 Myr) of these young stellar groups given in the literature (see also Sect. 4.4).

In the text
thumbnail Fig. 6.

HR diagram of the XFOR cluster compared to the grid of isochrones and tracks computed from different models. The solid red lines and black dashed lines indicate the isochrones and tracks for each model, with the age (in Myr) and mass (in M) values indicated in the panels.

In the text
thumbnail Fig. 7.

Absolute colour-magnitude diagram of the cluster members identified in our membership analysis (see Sect. 2). The lines indicate the empirical isochrone of the cluster and the 30 Myr isochrones computed from the BHAC15 and PARSEC models. The three stars located above the cluster sequence are Gaia EDR3 5056560315790459776 (033121–303059), Gaia EDR3 4854540344270707200 (032513–370909), and Gaia EDR3 4854879573671981056 (032215–354719).

In the text
thumbnail Fig. 8.

Absolute colour-magnitude diagram derived from Gaia-EDR3 photometry. The different colours and symbols indicate the XFOR cluster members identified in this study and known members of the COL, THA, and CAR stellar groups published in the literature (Gagné et al. 2018a,b; Gagné & Faherty 2018). The solid lines mark the 10 Myr and 40 Myr isochrones from the PARSEC models.

In the text
thumbnail Fig. 9.

3D spatial distribution of the stars in the cluster and the effects of mass segregation. The size and colour of the symbols are scaled with the spectral type of the stars derived in this study (see Sect. 3) in the sense that small red symbols denote the M dwarfs in our sample.

In the text
thumbnail Fig. 10.

Colour-colour diagram of cluster members using infrared photometry from the AllWISE catalogue. This is one of the diagrams used in the classification scheme developed by Koenig & Leisawitz (2014) to distinguish the SED subclasses based on infrared excess emission.

In the text
thumbnail Fig. 11.

3D spatial distribution of the cluster members identified in our membership analysis with SED classification derived in this study (see Sect. 4.5). The different colours and symbols denote the SED subclasses.

In the text
thumbnail Fig. 12.

Distribution of spectral types of the samples of stars in the XFOR cluster (this paper), IC 348 (Muench et al. 2007), Taurus (Esplin & Luhman 2019), and Upper Scorpius (Luhman et al. 2018). The dashed line indicates the completeness limit of the XFOR sample.

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.