Open Access
Issue
A&A
Volume 694, February 2025
Article Number A129
Number of page(s) 12
Section Stellar structure and evolution
DOI https://doi.org/10.1051/0004-6361/202450311
Published online 11 February 2025

© The Authors 2025

Licence Creative CommonsOpen Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

This article is published in open access under the Subscribe to Open model. Subscribe to A&A to support open access publication.

1. Introduction

The multiplicity of stars is essential in several areas of astronomy, such as stellar evolution, dynamics, and gravitational theories. Binary stars provide critical information that can be used to test and refine theoretical stellar models (Popper 1980; Andersen 1991; Southworth et al. 2005; Torres et al. 2010). The study of multiple systems sheds light on star formation processes and the dynamical effects, bringing on new insights into star formation and galactic evolution through their evolution over time. Binaries are key systems for measuring stellar masses accurately, even including radii, in cases of eclipsing binaries when applying Kepler’s third law. This information is fundamental to infer the stellar properties of each component and provides key inputs to theoretical models of stellar structure and evolution (Lodieu et al. 2020). The observational determination of masses, temperatures, and luminosities of stars allow for the validation or refining of existing models.

Metal-poor stars, also known as sub-dwarfs, are astronomical objects with abundances of elements heavier than hydrogen and helium in their atmospheres at levels that are at least ten times lower compared to the Sun and other Population I stars in the Milky Way. These stars have ages around 10–13 billion years, making them among the oldest objects in the Universe, having formed early on in the history of the Milky Way. They exhibit high proper motions and heliocentric velocities with kinematic properties typical of the thick disk or halo of the Galaxy with velocities typical larger than 60–100 km/s (Gizis 1997; Lépine et al. 2007; Jao et al. 2008). Subdwarfs typically exhibit lower luminosities than their solar counterparts and higher gravities, hence, smaller radii as well, compared to giants and subgiants. Metal-poor dwarfs usually bluer colours mainly due to enhanced calcium hydrides (CaH) in the optical and stronger collision-induced absorption-affected near-infrared (NIR) wavelengths. Sub-dwarfs provide valuable information about the conditions characterising the early universe and the processes that led to the formation of the first generation of stars in the Milky Way and other galaxies (Buser 2000; Bland-Hawthorn & Gerhard 2016; Gallart et al. 2019).

Stellar multiplicity is the result of a complex interplay of many physical processes, including metallicity, initial formation conditions, stellar evolution, and environment. Metal-poor environments tend to have lower opacities and cooling rates, which can affect the fragmentation process during star formation. Lower metallicity gas is less efficient at radiating away heat, leading to the formation of fewer but more massive fragments (Bate 2019; Myers 2014; Chon et al. 2021; Guszejnov et al. 2022; Bate 2023). Comprehensive studies of the multiplicity of metal-poor stars, including observational surveys, statistical analyses, and theoretical insights can be found in several key papers and references therein; for instance Abt & Cardona (1983), Carney (1983b,a), Latham et al. (2002), Raghavan et al. (2010), Moe & Di Stefano (2017), Gao et al. (2017), Badenes et al. (2018), El-Badry & Rix (2019), Liu (2019), Mazzola et al. (2020), Hwang et al. (2021), and Niu et al. (2022).

Several surveys have specifically looked into the close binary fraction of stars at different metallicity, with some controversial outcomes (e.g. Stryker et al. 1985; Carney et al. 1994; Latham et al. 2002), as well as differences settled in the work by Moe et al. (2018) where the metallicity does play a role in the frequency of close-in binaries with projected separations below 10 au. These authors showed that stellar interactions are an important factor that links binary frequency with metallicity. The overall frequency of F- and G-type sub-dwarfs (30–40%; Latham et al. 2002; Goldberg et al. 2002; Carney et al. 2003; Jao et al. 2009; Moe & Di Stefano 2017; El-Badry & Rix 2019) seems slightly lower than for their solar-metallicity counterparts (50–60% Duquennoy & Mayor 1991; Raghavan et al. 2010). Similarly, the global frequency of M sub-dwarf binaries is lower than their solar-metallicity counterparts (15% vs. 30–35%; Jao et al. 2009; Riaz et al. 2008) in addition to a lower number of wide co-moving pairs (Riaz et al. 2008; Jao et al. 2009; Lodieu et al. 2009; Ziegler et al. 2015); this is a trend that was previously observed for FGK sub-dwarfs as well (Hwang et al. 2021; Niu et al. 2022). The separation between components of M sub-dwarf multiples is lower and the mass ratio closer to unity than at solar-metallicity (Fischer & Marcy 1992; Reid & Gizis 1997; Delfosse et al. 2004; Dhital et al. 2010; Bergfors et al. 2010).

This paper presents a preliminary study of the multiplicity of metal-poor stars with metallicities below −1.5 dex over a large range of projected physical separations accessible to Gaia and seeing-limited images. In Section 2 we describe the sample extracted from the Gaia latest data release. In Section 3 we present several independent searches for wide companions: one based on Gaia data only, another one cross-matching the VISTA Hemisphere Survey (VHS; McMahon 2012) and Wide Field Infrared Survey Explorer (AllWISE; Wright et al. 2010), along with a lucky imaging survey of a sub-sample of stars. In Section 5 we put the numbers of common proper motion pairs identified in our search into context to determine the multiplicity of metal-poor stars and make comparisons with previous studies in the literature.

2. Sample selection

We launched an Astronomical Data Query Language (ADQL; Appendix A) query in the Gaia science archive1 to retrieve all metal-poor stars in the Gaia Data Release 3 (DR3; Gaia Collaboration 2016, 2023). We set the following constraints to identify good quality sources with physical parameters derived by Gaia (Creevey et al. 2023; Fouesneau et al. 2023), namely: distances below 400 pc with 3σ detections, metallicities below −1.5 dex, and transversal velocities higher than 60 km/s (Tian et al. 2020). The choice of the upper limit on the metallicity is our own choice to look at truly metal-poor dwarfs with little contamination from solar-metallicity sources due to the current uncertainties in the determination of metallicities that can be up to 1 dex. The query returned 15 432 sources (hereafter “our sample”), which can be easily downloaded from the Gaia archive by running the ADQL query provided in Appendix A.

The Gaia metallicities suffer from shortcomings (e.g. systematics and outliers) that other teams have looked into to derive metallicities with higher fidelity and fewer outliers (left-hand side plot in Fig. 2), particularly below −2.5 dex. As a consequence, we have cross-matched our sample with the catalogues published by Andrae et al. (2023)2 and Zhang et al. (2023)3. We note that we imposed a constraint on the quality flag in this catalogue: quality_flags ≤ 8 that contain good quality 175 and 220 million sources with Gaia spectra and more reliable metallicities. If we restrict these catalogues to the parameters set for our sample, we end up with 6500 and 1589 metal-poor ([Fe/H] ≤ −1.5 dex) sources with parallaxes larger than 2.5 mas with a 3σ confidence. After cross-matching our sample with both catalogues using a matching radius of 3 arcseconds (arcsecs), we identified 610 sources in common that constitute our main best-confidence sample of study in the rest of this work. The main difference with the original catalogue is the disappearance of the sources with high extinction (Fig. 1). The key properties of the samples are plotted in Fig. 1. Most of the stars in this sample have effective temperatures and gravities in the range 5000–6500 K (approximately mid-F to early-K; Pecaut & Mamajek 2013) and 3.8–4.6 dex, respectively (Andrae et al. 2023).

thumbnail Fig. 1.

Main parameters of the original (light gray) and final (gray) samples. Sources with high RUWE (greater or equal to 1.4) are highlighted in red. Multiple systems identified in the Gaia and VHS-AllWISE catalogues are marked with large black and blue symbols, respectively.

We also added sources from other catalogues of metal-poor stars freely available in the literature. In particular, we considered the Jinabase database of metal-poor stars (Abohalima & Frebel 2018)4 as well as the sample of 385 very metal-poor stars, including the spectroscopic follow-up (Li et al. 2022a). To ensure a reliable identification of common proper motion candidates in both samples, we limited our search to objects with total Gaia proper motions larger or equal to 200 mas/yr. Only half of the stars in the Jinabase catalogue satisfy this astrometric criterion, while just ten stars remain in the sample of Li et al. (2022a), most of them outside the VHS coverage. We kept both samples because of the trustworthy physical parameters and low metallicities, but we caution that the statistics on both catalogues is extremely limited.

3. Search for wide faint companions

3.1. Search for common proper motion companions

We cross-correlated internally our main sample of 610 sources varying the matching radius up to 300 arcsec. This search did not returned any wide pair. Gaia is expected to be complete for binaries with separation of a few arcsec down to its full depth. However, our cut in metallicity is strict at −1.5 dex and some true common proper motion candidates to sources in the main sample may exist with slightly higher metallicities or no [Fe/H] measurement at all.

We cross-matched these 610 sources with Gaia DR3 up to a separation of 25 arcsec, yielding 4309 candidates. After removing the primary sources themselves and keeping only common proper motion candidates at similar distances within 3σ of the error bars in parallax, μα cos δ, and μδ, we are left with nine pairs (Table B.1 in Appendix B; large black symbols in Fig. 1). The Gaia metallicity of the secondaries are either above the −1.5 dex limit or absent from the Gaia DR3 catalogue; this indicated the reason why they escaped the original selection. None of the primaries have large renormalised unit weight error (RUWE; Lindegren et al. 2018), whereas two secondaries have values larger than 1.25, suggesting that they might be binaries with separations below 0.02 arcsec that are unresolved by Gaia (Penoyre et al. 2020; Belokurov et al. 2020; Gandhi et al. 2022). The chance alignment is extremely small at these separations is displayed in Fig. 4 of El-Badry et al. (2021). Therefore, these systems might actually be higher order multiple. These systems seem to suffer little extinction (AG) as computed by Gaia in the G-band along the line-of-sight (histogram at the bottom-right panel in Fig. 1).

We found an optical spectrum of the source 00:49:20.33 +06:25:09.7 in the latest data releases of LAMOST5. The spectral type derived by LAMOST is A9–F0 with a temperature of 6075–6133 ± 34 K, a metallicity between −1.85 ± 0.01 and −1.77 ± 0.02 dex and a large radial velocity of about −188 ± 5 km/s. We found an APOGEE spectrum of 14:53:10.23 +48:14:43.8 in the latest release of the SDSS-III6 survey. The inferred parameters are Teff = 5759 ± 59 K, [Fe/H] = −1.83 ± 0.02 dex, and RV = −156.8 ± 0.2 km/s. Thus, the primaries of both systems have been spectroscopically confirmed as metal-poor dwarfs with metallicities in agreement with the ones in the Gaia catalogue to better than 0.3 dex.

3.2. Search for faint wide companions

The VISTA Hemisphere Survey (VHS; McMahon 2012) is a NIR survey of the entire Southern Sky (20 000 square degrees) using the 4m Visible and Infrared Survey Telescope for Astronomy (VISTA; Emerson 2001; Emerson et al. 2004) equipped with a camera made of sixteen 2048 × 2048 pixel detectors covering 1.65 square degrees, after dithering, with a mean pixel scale of 0.34 arcsec (Dalton et al. 2006). Half of the total area is covered in J and Ks while 5000 deg2 are covered in JHKs and another 5000 deg2 in YJHKs. The mean 5σ point source depth of the combined images are about Y = 20.0 mag, J = 19.5 mag, H = 18.6 mag, and Ks = 17.9 mag in the Vega system.

The All-Sky Wide-Field Infrared Survey Explorer (AllWISE) space-based mission covered the entire sky in four mid-infrared (MIR) passbands: 3.4 (W1), 4.6 (W2), 12 (W3), and 22 (W4) microns down to 17.1, 15.7, 11.7, 7.7 mag (95% completeness), respectively (Wright et al. 2010). It was a follow-up to the WISE mission, which operated from 2009 to 2011, with improved data processing and coverage. The total number of entries in the AllWISE source catalogue with W1 detections is close to 750 million sources, number decreasing with longer wavelength.

We cross-matched VHS and AllWISE to look for companions to metal-poor stars in the final common sample selection. We find ∼44% of these objects (268/610) covered by VHS. We prepared lists of objects from both catalogues that are located at a distance of less than 180 arcsec from each source in the selection, discarding objects that do not show motion signatures between VHS and AllWISE. We remained sensitive to companions about 1.5 times the full width at half maximum (FWHM) of WISE in the W1 and W2, which is about 10 arcsec. Then, we prepared lists of pair of sources with a separation between VHS and AllWISE less than 60 arcsec. This places an upper limit in proper motion of about 10 arcsec/yr. We computed the proper motion for each pair in the lists rejecting pairs with a difference in proper motion with the object in the Gaia sample larger than a 15%. The output lists contains mostly primaries from the Gaia sample along with contaminants and background extended objects, which were discarded after visual inspection. This search returned eight potential pairs with separations larger than 80 arcsec and J-band magnitude in the 17.6–20.3 mag range after excluding the primaries themselves, placing an upper limit of 8/268 = 3.0 ± 1.0% at the 1σ level on the frequency of wide companions. We list these companion candidates in Table C.1 in Appendix C and plot them as blue crosses in Fig. 1.

We repeated the search for faint companions around the individual catalogues to evaluate whether or not we could set a more stringent constraint on the multiplicity. We found 6580 stars of the 15 432 covered by VHS, whose distances range from 50 to 400 pc. In this sample, we identified one bright companion detected in Gaia DR3 at 4.85 arcsec from L 313-33 (Schneider et al. 2016), whose RUWE parameter is around 1.9 , thusindicating a potential multiple system. We did not find this secondary in our original Gaia sample because its metallicity is −0.56 dex in Gaia DR3, while the primary has a metallicity of −1.88 dex although proper motion and parallax are in agreement for both components. The metallicity of the pair is most likely not so low because it is not in the final sample of metal-poor stars discussed in this work. At larger separations beyond 20 arcsec from their primaries, we spotted another 12 companions also detected in Gaia DR3, all with metallicities higher than the lower limit of −1.5 dex and outside our final sample. Three out of 12 secondaries have RUWE parameters significantly larger, pointing towards multiple systems. At fainter magnitudes (J > 18 mag), we identified two wide companion candidates at 31.1 and 23.3 arcsec from TYC 7650-1108-1 with [Fe/H] = −2.38 dex) and CD-24 17504 with [Fe/H] = −4.04 dex (Ishigaki et al. 2012; Yoon et al. 2016), respectively. The latter is very metal-poor primary whose companion has been identified as part of the sample of Li et al. (2022a) but rejected with additional GTC observations (see below). We inferred an upper limit of 0.21 ± 0.06% ((1+12+2−1)/6580; 1σ) on the multiplicity rate for projected separations larger 500–72 000 au for the range of validity of the VHS-WISE search (10–180 arcsec) and distance interval.

We reiterated the cross-match process for the 6500 and 1589 metal-poor stars in the catalogues of Andrae et al. (2023) and Zhang et al. (2023), yielding 2671 and 665 counterparts in VHS. We identified 39 and 38 potential wide companions in Andrae et al. (2023) and Zhang et al. (2023), respectively. After further inspection and removal of sources detected in Gaia revealing distinct proper motions, we ended up with two and one companions confirmed with Gaia for the samples of Andrae et al. (2023) and Zhang et al. (2023). The system in common to both catalogues and our original selection is the LP 850-23/LP 850-24 system. The additional wide pair in the catalogue of Andrae et al. (2023) is L 158-59/L 158-60 with a separation of about 63 arcsec. Once again, the issue comes from the Gaia metallicities of one of the components, which is above our original cut in metallicity. Hence, the success rate of truly co-moving pairs of the VHS/AllWISE search is below 2/39 = 5.1% and 1/38 = 2.6% for the catalogues of Andrae et al. (2023) and Zhang et al. (2023), respectively.

We also identified 5 (17) and 3 (15) companion candidates without astrometric information in Gaia DR3 (no Gaia DR3 counterpart) in Andrae et al. (2023) and Zhang et al. (2023), respectively. All are faint in the NIR, with VHS J-band magnitudes typically fainter than 17.5 mag, implying G − J colours redder than 3.5 mag typically of ultracool dwarfs (Reylé 2018; Reylé et al. 2021). Considering the low success rate identifying companions in VHS-AllWISE above, we did not expect to find many bona fide companions in those lists without Gaia counterparts, a few at most. Therefore, we inferred a binary frequency of metal-poor dwarfs in the sample of Andrae et al. (2023) in the range of 1/2671 = 0.04% and (1+5+17)/2671 = 0.86%, while those percentages are around 0.15–2.9% for the sample from Zhang et al. (2023). We conclude that the overall binary frequency is at most 3%, most likely below 5% for projected separations larger than 1500–7200 au because of the success rate discussed above.

We repeated the same procedure to look for wide companions among the metal-poor stars in the catalogue of Li et al. (2022a) and the Jinabase sample covered by VHS. We identified a potential wide companion to 02:13:40.24 −00:05:18.6 in Li et al. (2022a); however we rejected it because we spotted an optical counterpart in the optical and no clear motion between the Sloan digital sky survey (SDSS; York et al. 2000) and Panoramic Survey Telescope and Rapid Response System (Pan-STARRS; Chambers et al. 2016) detections. We did detect one possible companion at a separation of 23.7 arcsec (23:07:22.25 +23.52:35.0) to CEMP-no CD-2417504 (Ishigaki et al. 2012; Yoon et al. 2016) in the Jinabase sample, respectively. We obtained CCD images of this companion candidate on 21 July 2023 with Optical System for Imaging and Low Resolution Integrated Spectroscopy (OSIRIS) camera-spectrograph (Cepa et al. 2000 in its new upgraded OSIRIS+ configuration7 on the Gran Telescopio de Canarias (GTC).

The OSIRIS+ detector is a blue ensitive monolithic 4k × 4k CCD, operated in 2 × 2 binned mode giving a plate scale of 0.25″/pix and providing a unvignetted FOV of 7.8′ × 7.8′. The standard operation mode of the instrument uses a readout speed of 233 kHz, which gives a gain of 1.9 electrons/ADU and readout noise of 4.3 electrons. We obtained a serie of 5 images using the Sloan g, r, i, and z filters with individual exposure times of 120 s with a small offset of 10 arcsec between them, giving a total on-source time of 600 s per filter. We reduced the images using an in-house pipeline for OSIRIS+ in imaging mode (SAUSERO), which comprises the standard procedures including bias, flat field, cosmic ray removal and a final combined image calibrated astrometrically in each passband with an average seeing of 0.9″. We did not detect any obvious motion between the OSIRIS and VHS epoch of the proposed companion to CD-2417504, yielding to the rejection of this candidate. The statistical multiplicity results are consistent with the former search due to the much lower number statistics.

3.3. Lucky imaging follow-up

We carried out high spatial resolution imaging of a sub-sample of binary system candidates from the original Gaia sample of 15 432 objects with FastCam (Oscoz et al. 2008), the “lucky imaging” instrument on the 1.52 m Carlos Sanchez telescope at Teide Observatory, Tenerife, Canary Islands in visitor mode on the nights of 18 October (observers EGE & RC) and 14 December 2023 (observer JYZ & RC). We note that none of these sources are common to the final sample of 610 metal-poor stars. We did not see any bias in the effective temperature and metallicity distribution of this sub-sample with respect to the full sample of 15 432 objects. The night on 18 October was affected by clouds whose thickness and presence became more significant towards the end of the night. The seeing was variable between 1.5 and 2.5 arcsec. The night on 14 December was dark and clear with no moon and no clouds but the peak of Geminides. The seeing was below 1.5 arcsec all night.

The upgraded FastCam camera is equipped with a Andor iXon DU-888U3-CSO#BV back-illuminated system8. The camera has 1024 × 1024 pixels with a size of 13 microns, low read-out noise (< 1 electron), and more than 95% quantum efficiency. We conducted observations in the broad-band Johnson I filter with a reduced field of view of 400 pixels aside to reduce the individual on-source integration times to about 15 ms. Typically, we selected short exposures of 15ms repeated 30 000 times for each target, yielding a total on-source integration of 7.5 min.

We targeted a total of 48 sub-dwarf candidates: 17 on 18 October and 31 on 14 December (Table D.1). We also targeted the Cyg 22 region and the Trapezium Cluster with several objects in the field to measure accurately the pixel scale on 18 October and 14 December 2023, respectively. We measured a pixel scale of 0.03493 and 0.03498 arcsec per pixel and rotation angle of 1.033 and 1.288 degree on 18 October and 14 December 2023, respectively.

We also collected high resolution images of nine targets with FastCam on the Nordic Optical Telescope (NOT) in service mode on the night of 22 and 23 January 2024 (observer RC). Among these, four were also observed with FastCam on the Carlos Sanchez (see above), while the other five had no first epoch with FastCam. Using the Trapezium Cluster as astrometric reference, we measured a pixel scale of 0.02521 arcsec per pixel and a rotation angle of 89.48 deg.

The data reduction of the raw images was carried out with a homemade pipeline developed specifically for FastCam. This pipeline is based on the “lucky images” technique, as described in Baldwin et al. (2001), Tubbs et al. (2002) and references therein. Images are chosen based on their Strehl ratio, as it can be regarded as a metric for assessing the impact of atmospheric turbulences on the image. The higher the Strehl ratio, the better the quality of the image. We selected the best 10% images to generate the final combined image shown in Appendix D. When the pipeline detects two stars in the image with similar magnitudes, it can utilise either of them for recentering, resulting in an image with a central star surrounded by two ghost stars. To address this issue, we devised an algorithm named “two stars”, which forces the pipeline to identify the pair of brightest speckles with a fixed separation and angular orientation, using them as a reference to recenter using always the same star. In particular, we note that 0219+30 appears as a false triple in the TCS/FastCam images; we confirm that this system is a binary, using NOT data with improved flux and resolution.

We measured the separation and angle of the binary systems detected in our sample with the ruler option on the ds9 fits viewer. The angle was calculated taking the peak of the primary as reference going from west to east. The typical uncertainty on the angle is a one degree, while the uncertainty on the separation is two tens of a pixel. We also give the separation in au by multiplying by the pixel scale of FastCam on the night of observations. We also report the ratio of the peak of the flux between the primary and the secondary in the I-band, which is a good proxy for the flux ratio. The flux of the primary is the taken from the pixel with the higher number of counts while the peak of the secondary comes either from the most illuminated pixel or the average of the four pixel with the largest fluxes. The typical error bar on the fluxes of the secondaries range from 20 to 200 for the best and worse cases. The targets with the largest uncertainties in all three parameters are 0927−08, 2007+13, and 2307+21. We note that 0426+15 might not be single while 0047−01 is either a binary or a triple.

4. Discussion

The final sample is the most complete sample presented in this work and relies on Gaia metallicities derived from three distinct techniques, that have indicated these sources are metal-poor with a high level of confidence. This sample contains mid-F to early-K stars with temperatures intervals of 4970–6470 K, 5370–6660 K, and 5380–6770 K from Andrae et al. (2023), Zhang et al. (2023), and Gaia DR3, respectively. The mean and median effective temperatures are 5773/5780 K, 5998/6002 K, and 6202/6271 K, respectively, showing that there are significant offsets between the determination of Andrae et al. (2023) and the other two. Naturally, the differences in temperatures impact the gravity and metallicity solutions. The lowest metallicities in Andrae et al. (2023), Zhang et al. (2023), and Gaia DR3 samples are −2.05, −1.85, and −4.0 dex with very comparable median values of −2.05, −1.85, and −1.98 dex, respectively. The significant difference in the lowest metallicities is obvious in Fig. 2. The gravities inferred in all catalogue are in agreement, with minimum and maximum values of log(g) of 3.6–4.9 dex. Hence, we did not divide the sample into sub-classes.

thumbnail Fig. 2.

Comparison of the Gaia metallicities (x-axis) with the determinations from Andrae et al. (2023) in red and from Zhang et al. (2023) in black.

We cross-matched the sample of 610 sub-dwarfs with several catalogues with spectroscopically derived physical parameters, with a special focus on the metallicity (the key parameters in this work). First, we considered the 400 very metal-poor stars studied with LAMOST and Subaru (Li et al. 2022a), which includes 11 objects in our sample (red triangles in Fig. 4). Second, we identified two sources in the catalogue of 1.2 million giants with abundances from LAMOST DR8 (Li et al. 2022b), which are plotted as green squares in Fig. 4. Third, we added the 89 sources with metallicities from APOGEE and PASTEL in Table 1 of Wang et al. (2022) and plotted as black crosses and dots in Fig. 4, respectively. We observe that the spectroscopic metallicities from APOGEE (Jönsson et al. 2020) tend to be lower by about 0.3 dex on average compared to the three estimates based on Gaia data, while the metallicities derived using the bibliographical catalogue PASTEL (Soubiran et al. 2010, 2016) appear higher (Wang et al. 2022). The spectroscopic determination of Li et al. (2022a) agree with the photometric estimates. Overall, the spectroscopic metallicities seem to point towards lower values. We conclude that the metallicities of our sample are accurate to ±0.5 dex in absolute terms, suggesting that the subsequent analysis considers truly metal-poor stars as those with metallicities below −1.0 dex with a high-level of confidence.

5. Multiplicity of metal-poor stars

In this section, we discuss the multiplicity of metal-poor stars based on the samples studied in Section 3.1. Hence, we aim to offer a global view on the multiplicity of metal-poor solar-type stars.

At separations closer than about 20 mas, Gaia is unable to separate stars. Assuming a distance interval of 50–400 pc for our final sample, this translates into projected physical separations in the ∼1–8 au range. The RUWE parameter is close to 1.0 for sources well fit by a single-star model. If the RUWE is larger than 1.25 in Gaia DR3 (Castro-Ginard et al. 2024), the astrometric solution is less compatible with a single-star model, implying that the probability that a system is, in fact, a multiple would increase. Assuming that all objects with RUWE ≥ 1.25 are binaries or higher-order multiples in Gaia DR3 (which is not necessarily the case; see below) as described in Castro-Ginard et al. (2024), we found 113 candidates out of 610. This would imply a multiplicity fraction of 113/610 = 18.5 ± 1.7% for separations below 1–8 au, assuming Poisson statistics. Our sample exhibits RUWE values between 1.15 and 1.36, implying an interval for the frequency of multiple systems of around 13.8–26.0% (Fig. 3). If we increase the lower limit of the RUWE parameter to 2.0, the number of sources decreases to 6.2% (38 sources). Hence, we conclude that the frequency of multiple systems increases swiftly for separations closer than ∼8 au (20 mas at a maximum distance of 400 pc), yielding a minimum multiplicity rate of 20 ± 6%. This number might be twice as high due to the Gaia incompleteness on the true fraction of binaries. If we focus on unresolved binaries (< 20 mas), which are better sampled by the percent of successful image parameter determination windows IPDfmp parameter (Lindegren et al. 2021) with more than one peak (IPDfmp > 2) of Gaia (Tokovinin 2023), we can thereby infer a fraction of 9/610 = 2.0 ± 0.5%.

thumbnail Fig. 3.

Histogram of the RUWE parameter across the sky for our sample of sub-dwarfs following the simulation from Castro-Ginard et al. (2024).

thumbnail Fig. 4.

Comparison of metallicities from different catalogues. Left: comparison of the Gaia metallicities (x-axis) with the spectroscopic determinations of Li et al. (2022a,b), and Wang et al. (2022), shown as red triangles, green squares, and grey symbols, respectively. Middle: same as middle panel but comparing the spectroscopic metallicities with the determination from Andrae et al. (2023) on the x-axis. Right: same as other panels but with metallicities from Zhang et al. (2023) on the x-axis.

We have double-checked the multiplicity cross-matching the original sample with the catalogues of Andrae et al. (2023) and Zhang et al. (2023) independently. We found 1129 and 695 sources, respectively, and RUWE (IPDfmp) percentages of 18.0 ± 1.3% (2.2 ± 0.4%) and 20.7 ± 1.7% (3.5 ± 0.7%), which is consistent with the former derivation within the Poisson errors. These multiplicities agree with the frequency of solar-type stars of about 20% for separations closer than 8 au derived from Figure 13 of Raghavan et al. (2010).

At intermediate separations of a few to hundreds of AU, the fraction has been found to remain constant and independent of metallicity (Köhler et al. 2000; Zinnecker et al. 2004). We identified nine wide co-moving pairs with separations between 1.5 and 22.5 arcsec, corresponding to projected physical separations in the 500–8700 au considering the distance of the primaries. We note that Gaia can detect approximately 60% of all binaries down to G = 20 mag closer than 250 pc to the Sun. This efficiency drops to 35% out to 1000 pc. If we interpolate linearly to the maximum distance of our sample (400 pc), the frequency would be complete to about 55%. Taking into account all these uncertainties, we derive a frequency of wide companions of 9/610 = 1.48 ± 0.49% for projected separations between 8 and 10 000 au (Table 5). Applying a correction for the Gaia incompleteness, this rate could be as high as 3.0% at 3σ. This rate is significantly lower than frequency of close-in binaries among solar-type FGK dwarfs (27–30% for separations of 8–2000 au; Raghavan et al. 2010), denoting a sharp decrease in the multiplicity of sub-dwarfs with separation.

thumbnail Fig. 5.

Properties of the sub-dwarf systems resolved with FastCam on the 1.5-m Carlos Sánchez telescope.

We did find nine co-moving pairs that could be resolved on seeing-limited images because the smaller separation is 1.5 arcsec (corresponding to about 500 au). This is at odds with the seeing-limited deep imaging surveys of metal-poor stars (Allen et al. 2000; Zapatero Osorio & Martín 2004) suggesting a fraction of ∼15% and solar-type stars (10.8% for a > 1000 au; Raghavan et al. 2010). One possible explanation for this is that we have specifically focussed on stars with metallicities below −1.5 dex. Allen et al. (2000) identified 122 wide binaries among about 1200 metal-poor stars (Schuster & Nissen 1988, 1989; Schuster et al. 1993) but only seven binaries have metallicities below −1.5 dex (see Table 2 of their appendix). Among these, the separations of the three most metal-poor stars with [Fe/H] below −2.0 dex lie between 32 000 and 125 000 au while the remaining ones are closer, with separations between 530 and 8600 au. Zapatero Osorio & Martín (2004) searched for companions up to 25 arcsec from a sample 473 metal-poor stars, with about half of them displaying metallicities below −1.5 dex. These authors reported six previously known companions and five new companions with separations between 37 and ∼56 500 au, yielding a frequency of about 4.5%. Many stars with [Fe/H] below −2.0 dex also show separations larger than 20 000 au, as in the sample of Allen et al. (2000). The wide binary frequency studied for a wide range of metallicities from −1.5 dex to +0.5 dex at separations of 10 000–100 000 au combining LAMOST spectra with astrometry from Gaia DR2 concluded that there must be a decrease in the multiplicity fraction with metallicity, reaching about 1% at [Fe/H] = −1.5 dex (Hwang et al. 2021), which is consistent with our results. Our sample is the largest homogeneous sample of metal-poor stars with metallicities below −1.5 dex with a dedicated analysis of the multiplicity. One possible explanation for the difference with earlier studies may come from the Gaia metallicities that may differ between the primary and secondary, as in the case of three of the nine pairs listed in Table B.1 in Appendix B, combined with the limited sensitivity of Gaia to very cool companions.

Overall, we find a sharp decrease in the frequency of multiple systems beyond 8 au (assuming the maximum distance of 400 pc for our sample). On the one hand, this decrease appears sharper than for solar-type stars. On the other hand, the frequency is significantly lower for metal-poor solar-type stars than their solar-metallicity counterparts (3.0% at 3σ confidence vs. 10–30% depending on the separation). However, the FastCam follow-up of 48 sub-dwarfs shows that a high value of the Gaia RUWE parameter does not always correlate with multiplicity. We find that 62.5 ± 2% of 48 sub-dwarfs are multiples in a range of projected separations of 105–840 au. Looking at the full sample of targets imaged with FastCam, we find that the numbers of resolved systems among low and high RUWE is comparable (48% vs. 52%), demonstrating that dedicated high-resolution imaging campaigns are needed to shed more light on the multiplicity as a function of metallicity and projected separations. These numbers suggest that the field population might rather come from low density environments rather than dense regions, as seen in open clusters and star-forming regions (Deacon & Kraus 2020).

Table 1.

Multiplicity statistics of our samples.

Among the sample of 49 sub-dwarfs targeted with FastCam, we resolved 30 systems with separations ranging between 0.7 and 2.1 arcsec, corresponding to projected separations between 21 and 840 au considering the range of distances of this FastCam sample (30–400 pc). We note that 13 of the 30 resolved systems had a RUWE higher than 1.4 (43.3%), while only 5 of the 18 unresolved systems had RUWE below 1.4 (27.8%). We emphasise that this lower limit for the RUWE parameter applied for Gaia DR2 prior to the updated analysis of Gaia DR3 by Castro-Ginard et al. (2024). The derived frequency is high with 63.2% of multiple systems with an uncertainty of about 2% because the multiplicity of one target (2307+21) is dubious. We emphasise that this FastCam sample is highly biased because we focussed mainly on candidates with a high probability of being multiples. Therefore, the derived frequency must be interpreted as a strict upper limit in an incomplete sample.

We plot several parameters in Fig. 5, including the Gaia RUWE and IPDfmp versus the physical separation in pixels and the ratio of the peak fluxes of the secondary divided by the primary. We do not observe any significant trend in the physical separation or flux ratios as a function of metallicity. The majority of systems seem to be concentrated around separation of 30–40 pixels (1.05–1.40 arcsec) and flux ratios of 0.2–0.4. However, we do observe a decline in the measured separation with increasing IDPfmp (lower right panel in Fig. 5). We cross-matched the list of systems observed with FastCam with Gaia DR3 to figure out how many pairs our original sample based on a strict limit on the metallicity may have lost. We find that 19 had no companion reported in Gaia DR3, 16 have a companion without effective temperature and metallicity estimates, 7 have a nearby object without parallax nor proper motion, 1 has two companions, 6 have more than one companion, whose metallicities are greater than −1.5 dex, but lower than −0.5 dex; this indicates that they are metal-poor, but with a metallicity discrepant between the primary and secondary. Among the pairs resolved with FastCam, most of them were indeed possible binaries in Gaia DR3, except for a handful of exceptions marked in Table D.1. We conclude that the multiplicity derived from the Gaia sample for separations below 25 arcsec is a lower limit because we are missing pairs due to various incompleteness factors; in particular, discrepant metallicities between both components of the system as well as secondaries with no parallax and proper motion measurements.

The multiplicity statistics set from the catalogues of Li et al. (2022a) and Jinabase sample agree with those extracted from the Gaia only sample, but they have a lower level of significance due to the smaller numbers of sources in the input catalogues (10 and 143, respectively). We inferred upper limits on the former of 10 ± 5% and of 0.70 ± 0.06% for projected separations of 50–4000 au and 80–1600 au, respectively.

We compare the multiplicity of the most statistically-significant sample with previous studies of metal-poor stars published in the literature. At close separations (< 8 au) based on the RUWE, our fractions of about 20% agree with the results of long-term radial velocity surveys of large samples of truly metal-poor stars (15%; Partridge 1967; Stryker et al. 1985; Carney et al. 1994; Latham et al. 2002). We conclude that initial conditions in the molecular clouds do not seem to evolve significantly over time.

6. Conclusions

We present dedicated searches for wide companions to three samples of stars with metallicities below [Fe/H] = −1.5 dex to look for low-mass companions and infer the multiplicity statistics. The main results of our analysis are summarised below:

  • Our sample involves 610 stars with metallicities below −1.5 dex and distances between 70 and 400 pc in common to three catalogues with independent determinations of the metallicities.

  • For our sample of sub-dwarfs, we infer a minimum binary fraction of 20 ± 6% based on the RUWE parameter for projected separations closer than ≤8 au, assuming a maximum distance of 400 pc.

  • Nine wide pairs separated by 1.5–25 arcsec are detected in the Gaia DR3 sample (500–8700 au), implying a fraction below 1.5 ± 0.5% with an upper limit of 3% at 3σ confidence level at odds with former studies. However, this fraction is a very likely lower limit because of various incompleteness factors that could be up to a factor of two due the completeness and the determination of metallicities from Gaia.

  • No pair with Gaia astrometry was found with separations of 25–300 arcsec (projected separations of 3750–120 000 au), yielding a very low fraction below 1% much lower with previous seeing-limited searches (4–5%).

  • Several wide common proper motion companion candidates have been identified in the 10–180 arcsec range in the VHS-AllWISE cross-match for three different samples of metal-poor stars with precise astrometry from Gaia. However, they require an additional astrometric and spectroscopic confirmation. We infer a multiplicity frequency below 3.0±1.0% (1σ confidence).

  • A sub-sample of systems with indication of multiple pairs targeted with lucky imaging revealed that 62.5% of them are pairs in the 105–840 au separation range, with statistics relatively independent of the RUWE parameter. There is a visible trend whereby the physical separation decreases with an increasing Gaia parameter IPDfmp.

  • High-resolution surveys are needed to complement the Gaia data and further constrain the impact of key astrometric parameters that are most sensitive to binarity.

The main conclusions of our study reveals that short-period binaries do not seem to be impacted by metallicity or by any of the specific conditions present in the original molecular clouds. The fraction of binaries with separations beyond 8 au appear less frequent at metallicities below −1.5 dex than at solar metallicity, pointing towards a significant disruption of the systems with less binding energy that is due to their old(er) age and increased probability of interactions. Complementary higher resolution imaging studies are needed to further constrain the multiplicity of this metal-poor population.

Data availability

Images of the targets observed with the FastCam instrument are available on zenodo at the URL https://zenodo.org/records/14558611.


Acknowledgments

NL and JYZ acknowledge support from the Agencia Estatal de Investigación del Ministerio de Ciencia e Innovación (AEI-MCINN) under grants PID2019-109522GB-C53 and PID2022-137241NB-C41. APG acknowledges support from the grant PID2020-120052GB-I00 financed by MCIN/AEI/10.13039/501100011033. ELM is supported by the European Research Council Advanced grant SUBSTELLAR, project number 101054354. We would like to thank the referee, Dr. Haijun Tian, for his comments on the original manuscript that led to a stronger and compelling analysis. We thank Alfred Castro-Ginard for sharing his latest GUMS simulation with us to look at the distribution across the sky (Castro-Ginard et al. 2024). This article is based on observations made with the 1.52m Carlos Sánchez telescope operated on the island of Tenerife by the IAC in the Spanish Observatorio del Teide (Canary Islands). This research has made use of the Simbad and Vizier databases, and the Aladin sky atlas operated at the centre de Données Astronomiques de Strasbourg (CDS), and of NASA’s Astrophysics Data System Bibliographic Services (ADS). This work presents results from the European Space Agency (ESA) space mission Gaia. Gaia data are being processed by the Gaia Data Processing and Analysis Consortium (DPAC). Funding for the DPAC is provided by national institutions, in particular the institutions participating in the Gaia MultiLateral Agreement (MLA). The Gaia mission website is https://www.cosmos.esa.int/gaia. The Gaia archive website is https://archives.esac.esa.int/gaia. Based on observations obtained as part of the VISTA Hemisphere Survey, ESO Progam, 179.A-2010 (PI: McMahon) Based on observations made with the Gran Telescopio Canarias (GTC), installed at the Spanish Observatorio del Roque de los Muchachos of the Instituto de Astrofísica de Canarias, on the island of La Palma. This work is (partly) based on data obtained with the instrument OSIRIS+, built by a Consortium led by the Instituto de Astrofísica de Canarias in collaboration with the Instituto de Astronomía of the Universidad Autónoma de México. OSIRIS was funded by GRANTECAN and the National Plan of Astronomy and Astrophysics of the Spanish Government. Based on observations made with the Nordic Optical Telescope, owned in collaboration by the University of Turku and Aarhus University, and operated jointly by Aarhus University, the University of Turku and the University of Oslo, representing Denmark, Finland and Norway, the University of Iceland and Stockholm University at the Observatorio del Roque de los Muchachos, La Palma, Spain, of the Instituto de Astrofisica de Canarias. 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, and NEOWISE, whicjklh is a project of the Jet Propulsion Laboratory/California Institute of Technology. WISE and NEOWISE are funded by the National Aeronautics and Space Administration. Guoshoujing Telescope (the Large Sky Area Multi-Object Fiber Spectroscopic Telescope LAMOST) is a National Major Scientific Project built by the Chinese Academy of Sciences. Funding for the project has been provided by the National Development and Reform Commission. LAMOST is operated and managed by the National Astronomical Observatories, Chinese Academy of Sciences.

References

  1. Abohalima, A., & Frebel, A. 2018, ApJS, 238, 36 [NASA ADS] [CrossRef] [Google Scholar]
  2. Abt, H. A., & Cardona, O. 1983, ApJ, 272, 182 [NASA ADS] [CrossRef] [Google Scholar]
  3. Allen, C., Poveda, A., & Herrera, M. A. 2000, A&A, 356, 529 [NASA ADS] [Google Scholar]
  4. Andersen, J. 1991, A&ARv, 3, 91 [Google Scholar]
  5. Andrae, R., Rix, H.-W., & Chandra, V. 2023, ApJS, 267, 8 [NASA ADS] [CrossRef] [Google Scholar]
  6. Badenes, C., Mazzola, C., Thompson, T. A., et al. 2018, ApJ, 854, 147 [NASA ADS] [CrossRef] [Google Scholar]
  7. Baldwin, J. E., Tubbs, R. N., Cox, G. C., et al. 2001, A&A, 368, L1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  8. Bate, M. R. 2019, MNRAS, 484, 2341 [NASA ADS] [CrossRef] [Google Scholar]
  9. Bate, M. R. 2023, MNRAS, 519, 688 [Google Scholar]
  10. Belokurov, V., Penoyre, Z., Oh, S., et al. 2020, MNRAS, 496, 1922 [Google Scholar]
  11. Bergfors, C., Brandner, W., Janson, M., et al. 2010, A&A, 520, A54 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  12. Bland-Hawthorn, J., & Gerhard, O. 2016, ARA&A, 54, 529 [Google Scholar]
  13. Buser, R. 2000, Science, 287, 69 [CrossRef] [Google Scholar]
  14. Carney, B. W. 1983a, AJ, 88, 623 [NASA ADS] [CrossRef] [Google Scholar]
  15. Carney, B. W. 1983b, AJ, 88, 610 [Google Scholar]
  16. Carney, B. W., Latham, D. W., Laird, J. B., & Aguilar, L. A. 1994, AJ, 107, 2240 [Google Scholar]
  17. Carney, B. W., Latham, D. W., Stefanik, R. P., Laird, J. B., & Morse, J. A. 2003, AJ, 125, 293 [NASA ADS] [CrossRef] [Google Scholar]
  18. Castro-Ginard, A., Penoyre, Z., Casey, A. R., et al. 2024, A&A, 688, A1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  19. Cepa, J., Aguiar, M., Escalera, V. G., et al. 2000, in SPIE Conf. Ser., eds. M. Iye, & A. F. Moorwood, 4008, 623 [Google Scholar]
  20. Chambers, K. C., Magnier, E. A., Metcalfe, N., et al. 2016, ArXiv e-prints [arXiv:1612.05560] [Google Scholar]
  21. Chon, S., Omukai, K., & Schneider, R. 2021, MNRAS, 508, 4175 [NASA ADS] [CrossRef] [Google Scholar]
  22. Creevey, O. L., Sordo, R., Pailler, F., et al. 2023, A&A, 674, A26 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  23. Dalton, G. B., Caldwell, M., Ward, A. K., et al. 2006, SPIE Conf. Ser., 6269 [Google Scholar]
  24. Deacon, N. R., & Kraus, A. L. 2020, MNRAS, 496, 5176 [NASA ADS] [CrossRef] [Google Scholar]
  25. Delfosse, X., Beuzit, J., Marchal, L., et al. 2004, in Spectroscopically and Spatially Resolving the Components of the Close Binary Stars, eds. R. W. Hilditch, H. Hensberge, & K. Pavlovski., ASP Conf. Ser., 318, 166 [Google Scholar]
  26. Dhital, S., West, A. A., Stassun, K. G., & Bochanski, J. J. 2010, AJ, 139, 2566 [NASA ADS] [CrossRef] [Google Scholar]
  27. Duquennoy, A., & Mayor, M. 1991, A&A, 248, 485 [NASA ADS] [Google Scholar]
  28. El-Badry, K., & Rix, H.-W. 2019, ApJ, 482, L139 [Google Scholar]
  29. El-Badry, K., Rix, H.-W., & Heintz, T. M. 2021, MNRAS, 506, 2269 [NASA ADS] [CrossRef] [Google Scholar]
  30. Emerson, J. P. 2001, in The New Era of Wide Field Astronomy, eds. R. Clowes, A. Adamson, & G. Bromage, ASP Conf. Ser., 232, 339 [NASA ADS] [Google Scholar]
  31. Emerson, J. P., Sutherland, W. J., McPherson, A. M., et al. 2004, The Messenger, 117, 27 [NASA ADS] [Google Scholar]
  32. Fischer, D. A., & Marcy, G. W. 1992, ApJ, 396, 178 [NASA ADS] [CrossRef] [Google Scholar]
  33. Fouesneau, M., Frémat, Y., Andrae, R., et al. 2023, A&A, 674, A28 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  34. Gaia Collaboration (Brown, A. G. A., et al.) 2016, A&A, 595, A2 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  35. Gaia Collaboration (Vallenari, A., et al.) 2023, A&A, 674, A1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  36. Gallart, C., Bernard, E. J., Brook, C. B., et al. 2019, Nature Astronomy, 3, 932 [NASA ADS] [CrossRef] [Google Scholar]
  37. Gandhi, P., Buckley, D. A. H., Charles, P. A., et al. 2022, MNRAS, 510, 3885 [CrossRef] [Google Scholar]
  38. Gao, S., Zhao, H., Yang, H., & Gao, R. 2017, ApJ, 469, L68 [Google Scholar]
  39. Gizis, J. E. 1997, AJ, 113, 806 [NASA ADS] [CrossRef] [Google Scholar]
  40. Goldberg, D., Mazeh, T., Latham, D. W., et al. 2002, AJ, 124, 1132 [NASA ADS] [CrossRef] [Google Scholar]
  41. Guszejnov, D., Grudić, M. Y., Offner, S. S. R., et al. 2022, MNRAS, 515, 4929 [Google Scholar]
  42. Hwang, H.-C., Ting, Y.-S., Schlaufman, K. C., Zakamska, N. L., & Wyse, R. F. G. 2021, MNRAS, 501, 4329 [NASA ADS] [CrossRef] [Google Scholar]
  43. Ishigaki, M. N., Chiba, M., & Aoki, W. 2012, ApJ, 753, 64 [NASA ADS] [CrossRef] [Google Scholar]
  44. Jao, W.-C., Henry, T. J., Beaulieu, T. D., & Subasavage, J. P. 2008, AJ, 136, 840 [NASA ADS] [CrossRef] [Google Scholar]
  45. Jao, W.-C., Mason, B. D., Hartkopf, W. I., Henry, T. J., & Ramos, S. N. 2009, AJ, 137, 3800 [NASA ADS] [CrossRef] [Google Scholar]
  46. Jönsson, H., Holtzman, J. A., Allende Prieto, C., et al. 2020, AJ, 160, 120 [Google Scholar]
  47. Köhler, R., Kunkel, M., Leinert, C., & Zinnecker, H. 2000, A&A, 356, 541 [NASA ADS] [Google Scholar]
  48. Latham, D. W., Stefanik, R. P., Torres, G., et al. 2002, AJ, 124, 1144 [NASA ADS] [CrossRef] [Google Scholar]
  49. Lépine, S., Rich, R. M., & Shara, M. M. 2007, ApJ, 669, 1235 [CrossRef] [Google Scholar]
  50. Li, H., Aoki, W., Matsuno, T., et al. 2022a, ApJ, 931, 147 [NASA ADS] [CrossRef] [Google Scholar]
  51. Li, Z., Zhao, G., Chen, Y., Liang, X., & Zhao, J. 2022b, MNRAS, 517, 4875 [NASA ADS] [CrossRef] [Google Scholar]
  52. Lindegren, L., Hernández, J., Bombrun, A., et al. 2018, A&A, 616, A2 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  53. Lindegren, L., Klioner, S. A., Hernández, J., et al. 2021, A&A, 649, A2 [EDP Sciences] [Google Scholar]
  54. Liu, C. 2019, MNRAS, 490, 550 [NASA ADS] [Google Scholar]
  55. Lodieu, N., Zapatero Osorio, M. R., & Martín, E. L. 2009, A&A, 499, 729 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  56. Lodieu, N., Paunzen, E., & Zejda, M. 2020, Low-Mass and Sub-stellar Eclipsing Binaries in Stellar Clusters, 213 [Google Scholar]
  57. Mazzola, C. N., Badenes, C., Moe, M., et al. 2020, MNRAS, 499, 1607 [CrossRef] [Google Scholar]
  58. McMahon, R. 2012, The Messenger, 150, 67 [NASA ADS] [Google Scholar]
  59. Moe, M., & Di Stefano, R. 2017, ApJS, 230, 15 [Google Scholar]
  60. Moe, M., Kratter, K. M., & Badenes, C. 2018, ApJ, 875, 61 [Google Scholar]
  61. Myers, P. C. 2014, ApJ, 781, 33 [NASA ADS] [CrossRef] [Google Scholar]
  62. Niu, Z., Yuan, H., Wang, Y., & Liu, J. 2022, ApJ, 931, 124 [NASA ADS] [CrossRef] [Google Scholar]
  63. Oscoz, A., Rebolo, R., López, R., et al. 2008, SPIE Conf. Ser., 7014 [Google Scholar]
  64. Partridge, R. B. 1967, AJ, 72, 713 [NASA ADS] [CrossRef] [Google Scholar]
  65. Pecaut, M. J., & Mamajek, E. E. 2013, ApJS, 208, 9 [Google Scholar]
  66. Penoyre, Z., Belokurov, V., Wyn Evans, N., Everall, A., & Koposov, S. E. 2020, MNRAS, 495, 321 [NASA ADS] [CrossRef] [Google Scholar]
  67. Popper, D. M. 1980, ARA&A, 18, 115 [NASA ADS] [CrossRef] [Google Scholar]
  68. Raghavan, D., McAlister, H. A., Henry, T. J., et al. 2010, ApJS, 190, 1 [Google Scholar]
  69. Reid, I. N., & Gizis, J. E. 1997, AJ, 113, 2246 [NASA ADS] [CrossRef] [Google Scholar]
  70. Reylé, C. 2018, ApJ, 619, L8 [Google Scholar]
  71. Reylé, C., Jardine, K., Fouqué, P., et al. 2021, A&A, 650, A201 [Google Scholar]
  72. Riaz, B., Gizis, J. E., & Samaddar, D. 2008, ApJ, 672, 1153 [NASA ADS] [CrossRef] [Google Scholar]
  73. Schneider, A. C., Greco, J., Cushing, M. C., et al. 2016, ApJ, 817, 112 [Google Scholar]
  74. Schuster, W. J., & Nissen, P. E. 1988, A&AS, 73, 225 [NASA ADS] [Google Scholar]
  75. Schuster, W. J., & Nissen, P. E. 1989, A&A, 221, 65 [NASA ADS] [Google Scholar]
  76. Schuster, W. J., Parrao, L., & Contreras de Martinez, M. E. 1993, A&AS, 97, 951 [NASA ADS] [Google Scholar]
  77. Soubiran, C., Le Campion, J. F., Cayrel Strobel, G., & Caillo, A. 2010, A&A, 515, A111 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  78. Soubiran, C., Le Campion, J.-F., Brouillet, N., & Chemin, L. 2016, A&A, 591, A118 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  79. Southworth, J., Maxted, P. F. L., & Smalley, B. 2005, A&A, 429, 645 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  80. Stryker, L. L., Hesser, J. E., Hill, G., Garlick, G. S., & Okeefe, L. M. 1985, PASP, 97, 247 [NASA ADS] [CrossRef] [Google Scholar]
  81. Tian, H.-J., El-Badry, K., Rix, H.-W., & Gould, A. 2020, ApJS, 246, 4 [Google Scholar]
  82. Tokovinin, A. 2023, AJ, 165, 180 [NASA ADS] [CrossRef] [Google Scholar]
  83. Torres, G., Andersen, J., & Giménez, A. 2010, A&ARv, 18, 67 [Google Scholar]
  84. Tubbs, R. N., Baldwin, J. E., Mackay, C. D., & Cox, G. C. 2002, ApJ, 387, L21 [Google Scholar]
  85. Wang, C., Huang, Y., Yuan, H., et al. 2022, ApJS, 259, 51 [NASA ADS] [CrossRef] [Google Scholar]
  86. Wright, E. L., Eisenhardt, P. R. M., Mainzer, A. K., et al. 2010, AJ, 140, 1868 [Google Scholar]
  87. Yoon, J., Beers, T. C., Placco, V. M., et al. 2016, ApJ, 833, 20 [Google Scholar]
  88. York, D. G., Adelman, J., Anderson, J. E., et al. 2000, AJ, 120, 1579 [Google Scholar]
  89. Zapatero Osorio, M. R., & Martín, E. L. 2004, MNRAS, 419, 167 [Google Scholar]
  90. Zhang, X., Green, G. M., & Rix, H.-W. 2023, MNRAS, 524, 1855 [NASA ADS] [CrossRef] [Google Scholar]
  91. Ziegler, C., Law, N. M., Baranec, C., Riddle, R. L., & Fuchs, J. T. 2015, ApJ, 804, 30 [NASA ADS] [CrossRef] [Google Scholar]
  92. Zinnecker, H., Köhler, R., & Jahreiß, H. 2004, in Revista Mexicana de Astronomia y Astrofisica Conference Series, eds. C. Allen, & C. Scarfe, 27, 33 [Google Scholar]

Appendix A: ADQL query in Gaia archive

SELECT ra, dec, l, b, parallax, parallax_error, pm, pmra, pmdec, ruwe, astrometric_excess_noise, astrometric_excess_noise_sig, phot_g_mean_mag, phot_bp_mean_mag, phot_rp_mean_mag, radial_velocity, teff_gspphot, logg_gspphot, mh_gspphot FROM gaiadr3.gaia_source WHERE parallax ≥ 2.5 AND parallax_over_error ≥ 3. AND sqrt(pmra*pmra+pmdec*pmdec)*4.74/parallax > = 60.0 AND mh_gspphot BETWEEN −10.0 AND −1.5 ORDER BY ra,dec

Appendix B: List of wide pairs up to 300 arcsec with some Gaia parameters

Table B.1.

List of wide pairs up to 300 arcsec with some Gaia parameters.

Appendix C: List of wide pairs in VHS and AllWISE.

Table C.1.

List of wide pairs up to 180 arcsec identified in VHS and AllWISE.

Appendix D: FastCam data

Table D.1.

List of systems followed-up with lucky imaging.

All Tables

Table 1.

Multiplicity statistics of our samples.

Table B.1.

List of wide pairs up to 300 arcsec with some Gaia parameters.

Table C.1.

List of wide pairs up to 180 arcsec identified in VHS and AllWISE.

Table D.1.

List of systems followed-up with lucky imaging.

All Figures

thumbnail Fig. 1.

Main parameters of the original (light gray) and final (gray) samples. Sources with high RUWE (greater or equal to 1.4) are highlighted in red. Multiple systems identified in the Gaia and VHS-AllWISE catalogues are marked with large black and blue symbols, respectively.

In the text
thumbnail Fig. 2.

Comparison of the Gaia metallicities (x-axis) with the determinations from Andrae et al. (2023) in red and from Zhang et al. (2023) in black.

In the text
thumbnail Fig. 3.

Histogram of the RUWE parameter across the sky for our sample of sub-dwarfs following the simulation from Castro-Ginard et al. (2024).

In the text
thumbnail Fig. 4.

Comparison of metallicities from different catalogues. Left: comparison of the Gaia metallicities (x-axis) with the spectroscopic determinations of Li et al. (2022a,b), and Wang et al. (2022), shown as red triangles, green squares, and grey symbols, respectively. Middle: same as middle panel but comparing the spectroscopic metallicities with the determination from Andrae et al. (2023) on the x-axis. Right: same as other panels but with metallicities from Zhang et al. (2023) on the x-axis.

In the text
thumbnail Fig. 5.

Properties of the sub-dwarf systems resolved with FastCam on the 1.5-m Carlos Sánchez telescope.

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.