Arxiv:2101.05282V3 [Astro-Ph.SR] 2 Feb 2021 at Angular Separations Greater Than About One Arcsecond, Wide Et Al
Total Page:16
File Type:pdf, Size:1020Kb
MNRAS 000,1–26 (2021) Preprint 3 February 2021 Compiled using MNRAS LATEX style file v3.0 A million binaries from Gaia eDR3: sample selection and validation of Gaia parallax uncertainties Kareem El-Badry,1,2¢ Hans-Walter Rix,2 and Tyler M. Heintz3 1Department of Astronomy and Theoretical Astrophysics Center, University of California Berkeley, Berkeley, CA 94720, USA 2Max-Planck Institute for Astronomy, Königstuhl 17, D-69117 Heidelberg, Germany 3Department of Astronomy, Boston University, 725 Commonwealth Ave., Boston, MA 02215, US Accepted to MNRAS ABSTRACT We construct from Gaia eDR3 an extensive catalog of spatially resolved binary stars within ≈ 1 kpc of the Sun, with projected separations ranging from a few au to 1 pc. We estimate the probability that each pair is a chance alignment empirically, using the Gaia catalog itself to calculate the rate of chance alignments as a function of observables. The catalog contains 1.3 (1.1) million binaries with >90% (>99%) probability of being bound, including 16,000 white dwarf – main sequence (WD+MS) binaries and 1,400 WD+WD binaries. We make the full catalog publicly available, as well as the queries and code to produce it. We then use this sample to calibrate the published Gaia DR3 parallax uncertainties, making use of the binary components’ near-identical parallaxes. We show that these uncertainties are generally reliable for faint stars (퐺 & 18), but are underestimated significantly for brighter stars. The underestimates are generally ≤ 30% for isolated sources with well-behaved astrometry, but are larger (up to ∼80%) for apparently well-behaved sources with a companion within . 4 arcsec, and much larger for sources with poor astrometric fits. We provide an empirical fitting function to inflate published fs values for isolated sources. The public catalog offers wide ranging follow-up opportunities: from calibrating spectroscopic surveys, to precisely constraining ages of field stars, to the masses and the initial-final mass relation of white dwarfs, to dynamically probing the Galactic tidal field. Key words: binaries: visual – stars: evolution – methods: statistical – catalogues – parallaxes 1 INTRODUCTION genuine binaries, in which the two stars have nearly identical proper motions (e.g. Luyten 1971, 1979; Salim & Gould 2003; Chanamé & About half of all solar-type stars are members of binary systems, and Gould 2004; Dhital et al. 2010). Many wide binary searches have a majority of these are so widely separated that the two components specifically targeted high-proper motion stars, which have fewer never interact (e.g. Moe & Di Stefano 2017). With orbital periods phase-space neighbors that can be mistaken for binary companions. ranging from ∼10 to ∼ 108 years, most of these binaries can in If available, parallaxes and radial velocities are also useful for dis- some sense be viewed as clusters of two: the components formed tinguishing binaries from chance alignments (e.g. Close et al. 1990; from the same gas cloud and have orbited one another ever since. Andrews et al. 2017). They thus have essentially the same age, initial composition, and distance, but generally different masses and occasionally different Prior to the Gaia mission (Gaia Collaboration et al. 2016), useful evolutionary phases. This makes wide binaries useful for calibrating parallaxes for this purpose were only available for (relatively) small stellar models as well as spectroscopic and astrometric surveys. samples of nearby and bright stars. Gaia DR2 (Gaia Collaboration arXiv:2101.05282v3 [astro-ph.SR] 2 Feb 2021 At angular separations greater than about one arcsecond, wide et al. 2018) dramatically expanded the sample of stars with well- binaries are easily resolvable as two point sources. Distinguishing measured parallaxes and proper motions, enabling the construction of physically bound binary stars from chance alignments (“optical dou- unprecedentedly pure and extensive wide binary samples. El-Badry bles”) has been a long-standing challenge for binary star astronomy. & Rix(2018, hereafter ER18) searched Gaia DR2 for pairs of stars Indeed, the first systematic binary star catalog was constructed under within 200 pc of the Sun with parallaxes and proper motions consis- the assumption that all close pairs were chance alignments (Her- tent with being gravitationally bound, and projected separations of schel 1782), an assumption that was only shown to be incorrect two up to 50,000 AU (0.24 pc). Their catalog prioritized purity over com- decades later (Herschel 1803). pleteness, and thus imposed relatively strict cuts on astrometric and For bright binaries at close angular separations, chance alignments photometric quality and signal-to-noise ratio (SNR). This resulted in can be excluded probabilistically. However, the contamination rate a catalog of ∼ 55, 000 binaries with an estimated contamination rate from chance alignments increases at wider separations and fainter of ∼ 0.1%. magnitudes. Inclusion of proper motion data can aid the selection of Using the same basic strategy but less stringent cuts on astrometric SNR, Tian et al.(2020) extended the ER18 binary search to larger distances (3 < 4 kpc) and wider separations (B < 1 pc). This pro- ¢ E-mail: [email protected] duced a substantially larger sample of ∼ 800, 000 binary candidates, © 2021 The Authors 2 El-Badry et al. but with a higher contamination rate: chance-alignments dominate and parallax_over_error > 5 their catalog at B & 20, 000 au, though higher-purity subsamples can and parallax_error < 2 be selected by imposing stricter cuts on astrometric SNR. Their cat- and phot_g_mean_mag is not null alog contained 325,000 binaries with B < 20, 000 au; it is expected to be reasonably pure in this separation regime. The query returns a total of # = 64, 407, 853 sources, corresponding Another binary catalog was produced by Hartman & Lépine to # ¹# ¸ 1º/2 ≈ 2 × 1015 possible pairs. Of these, we consider as (2020), who combined Gaia DR2 astrometry with a catalog of high- initial binary candidates all pairs that satisfy the following: proper motion stars not contained in Gaia DR2. In order to reduce • Projected separation less than 1 parsec: the angular separation contamination from chance alignments, they limited their search to between the two stars, \, must satisfy binaries with proper motions larger than 40 mas yr−1; this translates \ s to a distance limit of order 200 pc for typical stars in the Galactic ≤ 206.265 × 1 , (1) disk, but to a larger search volume for stars on halo-like orbits with arcsec mas large tangential velocities. Their primary binary catalog contains where s1 is the parallax of the star with the brighter 퐺 magnitude. ∼ 100, 000 binary candidates, with projected separations as large as The maximum search radius of 1 pc (corresponding to an orbital 10 pc. period of ∼ 108 years) is chosen because a vanishingly small number In this paper, we use Gaia eDR3 data (Gaia Collaboration et al. of bound binaries are expected to exist at separations wider than this, 2020a) to further expand the sample of known wide binaries. Com- where the Galactic tidal field becomes comparable to the gravita- pared to DR2, eDR3 astrometry is based on a ∼1.5 × longer time tional attraction of the two stars (e.g. Binney & Tremaine 2008). The baseline. This yields significant improvements in parallax and espe- separation beyond which the Galactic tidal field dominates a binary’s cially proper motion uncertainties. For example, the median uncer- internal acceleration is called the Jacobi radius. In the Solar neigh- parallax pmra tainties in and (one-dimenstional proper motion) at borhood, it is given by A퐽 ≈ 1.35 pc × ¹"tot/" º, where "tot is 퐺 = 18 have respectively improved from 0.165 mas to 0.120 mas, the total mass of the binary (Jiang & Tremaine 2010). At separations −1 −1 and from 0.280 mas yr to 0.123 mas yr (Lindegren et al. 2018, slightly below A퐽 , binaries are efficiently disrupted by gravitational 2020b). Improvements at the bright end are more significant, due to perturbations from objects such as other stars and molecular clouds better handling of systematics; i.e., at 퐺 = 13, the median parallax (e.g. Weinberg et al. 1987). uncertainty decreased from 0.029 mas to 0.015 mas. This improved • Parallaxes consistent within 3 (or 6) sigma: the parallaxes of astrometric precision allows us to distinguish bound binaries from the two components, s1 and s2, must satisfy chance alignments to larger distances and wider separations than q was possible with DR2. Our approach is a compromise between the j − j 2 ¸ 2 s1 s2 < 1 fs,1 fs,2, (2) strategies adopted by Tian et al.(2020) and ER18. Like Tian et al. (2020), we search out to wide separations and relatively low astro- where fs,8 is the parallax uncertainty of the 8-th component, and metric SNR, so that at wide separations and faint magnitudes, the 1 = 3 for pairs with \¡ 4 arcsec, or 1 = 6 for pairs with \ < 4 full catalog is dominated by chance alignments. However, we also arcsec. The less stringent cut at \ < 4 arcsec is adopted because empirically estimate and assign each binary a probability that it is a the chance alignment rate there is low and parallax uncertainties are chance alignment, making it straightforward to select pure subsets of significantly underestimated at close angular separations (Section5). the catalog. • Proper motions consistent with a Keplerian orbit: The two stars The remainder of the paper is organized as follows. In Section2, we in a wide binary will have proper motions that are similar, but, due describe how our binary candidate sample is selected and cleaned. to orbital motion, not identical. We require Section3 describes how we quantify the contamination rate from Δ` ≤ Δ` ¸ 2f , (3) chance alignments and estimate the probability that each binary can- orbit Δ` didate is bound. Section4 details basic properties of the catalog and where Δ` is the observed scalar proper motion difference, fΔ` its a cross-match with the LAMOST survey.