We report spectroscopic and photometric followup of a dormant black hole BH candidate from Gaia DR3. The system which we call Gaia BH2 contains a 1 MSUBSUB red giant and a dark companion with mass M2  8.9pm 0.3 rm Modot  that is very likely a BH. The orbital period PSUBorbSUB  1277 d is much longer than that of any previously studied BH binary. Our radial velocity RV followup over a 7month period spans gt90 per cent of the orbits RV range and is in excellent agreement with the Gaia solution. UV imaging and highresolution optical spectra rule out plausible luminous companions that could explain the orbit. The star is a bright G  12.3 slightly metalpoor rm FeH0.22 lowluminosity giant Trm eff4600 rm K R  7.8 Rodot  log leftgleftrm cm s2rightright  2.6. The binarys orbit is moderately eccentric e  0.52. The giant is enhanced in elements with rm alpha Fe  0.26 but the systems Galactocentric orbit is typical of the thin disc. We obtained Xray and radio nondetections of the source near periastron which support BH accretion models in which the net accretion rate at the horizon is much lower than the BondiHoyleLyttleton rate. At a distance of 1.16 kpc Gaia BH2 is the secondnearest known BH after Gaia BH1. Its orbit  like that of Gaia BH1  seems too wide to have formed through common envelope evolution. Gaia BH1 and BH2 have orbital periods at opposite edges of the Gaia DR3 sensitivity curve perhaps hinting at a bimodal intrinsic period distribution for wide BH binaries. Dormant BH binaries like Gaia BH1 and Gaia BH2 significantly outnumber their close Xray bright cousins but their formation pathways remain uncertain.
['binaries: spectroscopic', 'stars: black holes', 'Astrophysics - Solar and Stellar Astrophysics', 'Astrophysics - Astrophysics of Galaxies']
A red giant orbiting a black hole
{'A red giant orbiting a black hole': 'Kareem El-Badry 1 GLYPH<150> 2 GLYPH<150> 3 GLYPH<150> 4 ★ , Hans-Walter Rix 3 , Yvette Cendes 1 , Antonio C. Rodriguez 4 , Charlie Conroy 1 , Eliot Quataert 5 , Keith Hawkins 6 , Eleonora Zari 3 , Melissa Hobson 3 , Katelyn Breivik 7 , Arne Rau 8 , Edo Berger 1 , Sahar Shahaf 9 , Rhys Seeburger 3 , Kevin B. Burdge 10 , David W. Latham 1 , Lars A. Buchhave 11 , Allyson Bieryla 1 , Dolev Bashi 12 , Tsevi Mazeh 12 , Simchon Faigler 12  \n1 Center for Astrophysics j Harvard & Smithsonian, 60 Garden Street, Cambridge, MA 02138, USA  \n- 2 Harvard Society of Fellows, 78 Mount Auburn Street, Cambridge, MA 02138\n- 3 Max-Planck Institute for Astronomy, Königstuhl 17, D-69117 Heidelberg, Germany\n- 4 Department of Astronomy, California Institute of Technology, Pasadena, CA 91125, USA  \n5 Department of Astrophysical Sciences, Princeton University, Princeton, NJ 08544, USA  \n- 6 Department of Astronomy, The University of Texas at Austin, 2515 Speedway Boulevard, Austin, TX 78712, USA\n- 7 Center for Computational Astrophysics, Flatiron Institute, 162 Fifth Ave, New York, NY, 10010, USA\n- 8\n- Max-Planck-Institut für extraterrestrische Physik, Gießenbachstraße 1, 85748 Garching, Germany  \n9 Department of Particle Physics and Astrophysics, Weizmann Institute of Science, Rehovot 7610001, Israel  \n- 10 MIT-Kavli Institute for Astrophysics and Space Research 77 Massachusetts Ave. Cambridge, MA 02139, USA\n- 11 DTU Space, National Space Institute, Technical University of Denmark, Elektrovej 328, DK-2800 Kgs. Lyngby, Denmark\n- 12 School of Physics and Astronomy, Tel Aviv University, Tel Aviv, 6997801, Israel', 'ABSTRACT': "We report spectroscopic and photometric follow-up of a dormant black hole (BH) candidate from Gaia DR3. The system, which we call Gaia BH2, contains a GLYPH<24> 1 𝑀 GLYPH<12> red giant and a dark companion with mass 𝑀 2 = 8 GLYPH<149> 9 GLYPH<6> 0 GLYPH<149> 3 𝑀 GLYPH<12> that is very likely a BH. The orbital period, 𝑃 orb = 1277 days, is much longer than that of any previously studied BH binary. Our radial velocity (RV) follow-up over a 7-month period spans ¡ 90% of the orbit's RV range and is in excellent agreement with the Gaia solution. UV imaging and high-resolution optical spectra rule out plausible luminous companions that could explain the orbit. The star is a bright ( 𝐺 = 12 GLYPH<149> 3), slightly metal-poor ( » Fe GLYPH<157> H … = GLYPH<0> 0 GLYPH<149> 22) low-luminosity giant ( 𝑇 eff = 4600 K; 𝑅 = 7 GLYPH<149> 8 𝑅 GLYPH<12> ; log GLYPH<2> 𝑔 GLYPH<157> GLYPH<0> cms GLYPH<0> 2 GLYPH<1> GLYPH<3> = 2 GLYPH<149> 6). The binary's orbit is moderately eccentric ( 𝑒 = 0 GLYPH<149> 52). The giant is enhanced in 𝛼 GLYPH<0> elements, with » 𝛼 GLYPH<157> Fe … = , 0 GLYPH<149> 26, but the system's Galactocentric orbit is typical of the thin disk. We obtained X-ray and radio nondetections of the source near periastron, which support BH accretion models in which the net accretion rate at the horizon is much lower than the Bondi-Hoyle-Lyttleton rate. At a distance of 1.16 kpc, Gaia BH2 is the second-nearest known BH, after Gaia BH1. Its orbit - like that of Gaia BH1 seems too wide to have formed through common envelope evolution. Gaia BH1 and BH2 have orbital periods at opposite edges of the Gaia DR3 sensitivity curve, perhaps hinting at a bimodal intrinsic period distribution for wide BH binaries. Dormant BH binaries like Gaia BH1 and Gaia BH2 significantly outnumber their close, X-ray bright cousins, but their formation pathways remain uncertain.  \nKey words: binaries: spectroscopic - stars: black holes", '1 INTRODUCTION': "The Milky Way is very likely teeming with stellar-mass black holes (BHs). We can see today that it contains a few GLYPH<2> 10 4 O stars with masses 𝑀 & 20 𝑀 GLYPH<12> (Garmany et al. 1982; Reed 2003), a majority of which will likely leave behind BHs when they die. The lifetime of these massive stars is of order 10 Myr, or 0.1% the age of the Milky Way. Notwithstanding variations in the Galactic star formation rate with cosmic time, this implies that the Milky Way should contain GLYPH<24> 1000 BHs for every luminous BH progenitor alive today; i.e., a few GLYPH<2> 10 7 BHs.  \nOnly a tiny fraction of this BH population has been observed. Observations of X-ray binaries have dynamically confirmed the presence of a BH in GLYPH<24> 20 systems (Remillard & McClintock 2006), and another  \nGLYPH<24> 50 X-ray transients are suspected to contain BHs (Corral-Santana et al. 2016). In all these systems, there is ongoing mass transfer from a luminous star onto a BH, giving rise to an accretion disk that shines brightly in X-rays. BH X-ray binaries are an intrinsically rare outcome of binary evolution (e.g. Portegies Zwart et al. 1997), and empirical population models suggest that only about 1000 exist in the Milky Way (Corral-Santana et al. 2016). They nevertheless make up the bulk of the observed BH population, because they are easier to find than BHs that are not accreting.  \nThe Gaia mission has opened a new window on the Milky Way's binary star population. By precisely measuring the positions of GLYPH<24> 2 billion sources over time, Gaia can detect subtle wobbles caused by the gravitational effects of binary companions - even companions that do not emit any light themselves. As of the mission's 3rd data release ('DR3'; Gaia Collaboration et al. 2022c), nearly 500 million sources have been observed in more than 20 well-separated epochs. These  \ndata are ultimately expected to yield astrometric binary solutions for of order 10 7 binaries (e.g. Söderhjelm 2004), but the quality cuts applied to solutions included in DR3 were stringent, so that only 1 GLYPH<149> 7 GLYPH<2> 10 5 full astrometric orbital solutions have been published so far. The DR3 binary catalog nevertheless represents more than an order of magnitude increase in sample size over previous samples of binary orbits and is therefore a promising dataset to search for rare objects.  \nGaia DR3 has thus far resulted in the identification of one unambiguous BH binary (El-Badry et al. 2023; Chakrabarti et al. 2022). That system, Gaia BH1, consists of a Sun-like star in a wide ( 𝑃 orb = 186 day; 𝑎 = 1 GLYPH<149> 4 au) orbit with a GLYPH<24> 10 𝑀 GLYPH<12> BH companion. The binary's wide orbit is not easily explained by standard binary evolution models. At a distance of only 𝑑 = 480 pc, Gaia BH1 is a factor of GLYPH<24> 3 closer than the next-nearest known BH, and so the Copernican principle suggests that widely separated BH binaries similar to Gaia BH1 are likely common compared to BHs in X-ray binaries. Most theoretical models also predict that wide binaries with non-accreting BHs should significantly outnumber their accreting cousins (e.g. Breivik et al. 2017; Langer et al. 2020).  \nGaia data are well suited to identify such systems. While Gaia BH1 was the most obvious BH binary candidate published in DR3, El-Badry et al. (2023) also identified a second promising candidate, whose longer orbital period demanded longer-term spectroscopic follow-up than was necessary for Gaia BH1. This paper is focused on that second object, which we refer to as Gaia BH2. Similar to Gaia BH1, the system contains a solar-mass star in a wide orbit around an object we suspect to be a BH. Unlike in Gaia BH1, the luminous star has left the main sequence and is ascending the giant branch. At a distance of 1.16 kpc, the object is likely the second-nearest known BH, after Gaia BH1.  \nThe remainder of this paper is organized as follows. Section 2 describes the system's discovery and Gaia orbital solution. Most of the observational data and modeling are described in Section 3, including light curves (Section 3.1), extinction (Section 3.2), UV photometry (Section 3.3), SED modeling (Section 3.4), mass inference of the giant (Section 3.5), spectroscopy (Section 3.6), orbit modeling (Section 3.7), spectral analysis (Section 3.8), Galactic orbit (Section 3.9), X-ray observations (Section 3.10), and radio observations (Section 3.11). We discuss the nature of the unseen companion in Section 4 and its possible formation history in Section 5. Prospects for detecting and characterizing additional BHs in binaries with Gaia are discussed in Section 6, and we conclude in Section 7. The appendices provide additional details about several aspects of the data and modeling.", '2 DISCOVERY': "Being a bright optical source ( 𝐺 = 12 GLYPH<149> 3), Gaia BH2 appears in many archival photometric catalogs. To our knowledge, no spectroscopic observations of the source occurred prior to the Gaia observations. An orbital solution with a best-fit period of 1352 days - somewhat longer than the GLYPH<24> 1000 day time baseline for observations included in Gaia DR3-was published in Gaia DR3(Gaia Collaboration et al. 2022a) with source id 5870569352746779008. We discuss the cadence and phase coverage of the Gaia observations in Appendix A, and the astrometric gooodness of fit in Appendix B. Both are unproblematic.  \nUnlike Gaia BH1, which has a purely astrometric solution ( nss\\_solution\\_type = Orbital ), Gaia BH2 has a joint astrometric + single-lined spectroscopic solution ( nss\\_solution\\_type  \n= AstroSpectroSB1 ). The reason for this difference is that individual-epoch RVs are only derived for sources brighter than 𝐺 RVS = 12, where 𝐺 RVS is a source's magnitude in the Gaia RVS bandpass. Gaia BH2 has 𝐺 RVS = 11 GLYPH<149> 11, and RVs were measured (but not published) in 34 unique transits. Gaia BH1 has 𝐺 RVS = 12 GLYPH<149> 83, and thus no individual-epoch RVs.  \nWeidentified Gaia BH2 as a promising BH candidate in the course of a broader program to spectroscopically follow-up Gaia DR3 binaries suspected to contain compact objects. The astrometric selection that first brought the system to our attention is described in Appendix E of El-Badry et al. (2023). The Gaia solution has a period of 𝑃 orb = 1352 GLYPH<6> 46 days. The astrometric mass function implied by the astrometric data alone is 𝑓 ' 𝑀 2 ' ast = ' 𝑎 0 GLYPH<157> 𝜛 ' 3 ' 𝑃 orb GLYPH<157> yr ' GLYPH<0> 2 = 6 GLYPH<149> 72 GLYPH<6> 0 GLYPH<149> 5 𝑀 GLYPH<12> ; here 𝑎 0 and 𝜛 represent the semi-major axis of the photocenter and the parallax. Assuming a mass of 1 𝑀 GLYPH<12> for the luminous star (Section 3.5), this implied an unseen companion mass of 𝑀 2 GLYPH<25> 8 GLYPH<149> 4 GLYPH<6> 0 GLYPH<149> 5 𝑀 GLYPH<12> . Since this far exceeds the apparent photometric mass of the luminous primary, we initiated a spectroscopic follow-up campaign, obtaining our first spectrum of the source in August 2022.  \nIn El-Badry et al. (2023), we reported that early follow-up data were consistent with the Gaia solution. However, since the data inhand then covered only a small fraction of the orbit's predicted RV range, we cautioned that continued RV monitoring was required to better assess the reliability of the Gaia solution. We have now carried out the necessary observations.  \nGaia data for the source were also analyzed by Tanikawa et al. (2022), who were the first to draw attention to the source as a BH candidate in a September 2022 preprint. They concluded that the Gaia orbital solution makes the source the most promising BH candidate among the GLYPH<24> 64000 binary solutions in Gaia DR3 for which both time-resolved spectroscopic and astrometric data are available. Their work showed convincingly that the Gaia data suggest the presence of an unseen massive companion, which would be most simply explained by a BH. These authors did not obtain any follow-up data, so their conclusions depend entirely on the reliability of the Gaia solution.", '3 PROPERTIES OF THE SOURCE': "Gaia BH2 is a source in Centaurus ( 𝛼 = 13:50:16.7; 𝛿 = GLYPH<0> 59:14:20), near the Galactic plane ( 𝑙 = 310 GLYPH<149> 4 GLYPH<150> 𝑏 = 2 GLYPH<149> 8), with 𝐺 = 12 GLYPH<149> 3. Its basic properties are summarized in Figure 1. The field is relatively crowded, with 17 Gaia -detected neighbors within 10 arcsec, though all are at least 3.5 magnitudes fainter than the source itself. The nearest Gaia -detected neighbor is a 𝐺 = 17 GLYPH<149> 95 source at a distance of 1.77 arcsec. This source is unresolved in the postage stamp shown in Figure 1, which is from the DECaPS survey (Schlafly et al. 2018; Saydjari et al. 2022). Its astrometry is inconsistent with being bound to Gaia BH2.  \nThe Gaia astrometric solution places Gaia BH2 at a distance 𝑑 = 1 GLYPH<149> 16 GLYPH<6> 0 GLYPH<149> 02 kpc. We do not attempt to correct this distance for the Gaia parallax zeropoint (Lindegren et al. 2021b), since there is little reason to expect that this zeropoint - inferred from single-star astrometric solutions - applies to the binary solutions. On the color-magnitude diagram, Gaia BH2 appears on the lower giant branch, just below the red clump. Our constraints on the star's temperature and radius (Section 3.4) imply a mass near 1 𝑀 GLYPH<12> , as shown in the right panel of Figure 1.  \n<!-- image -->  \n<!-- image -->  \n-  \nFigure 1. Properties of the luminous star. Left panel shows a 1.1-arcmin 𝑔𝑟𝑧 -band postage stamp from the DECaPS survey. Gaia BH2 is the bright (saturated) source at the center. Middle panel shows the source's position on the dereddened Gaia color-magnitude diagram, with the full Gaia DR3 astrometric binary sample shown for context. Right panel compares the star's inferred temperature and radius to single-star MIST evolutionary models. The optical source appears to be an evolved GLYPH<24> 1 𝑀 GLYPH<12> star on the lower giant branch.  \n<!-- image -->", '3.1 Light curves': "WeinspectedtheASAS-SN 𝑉 GLYPH<0> and 𝑔 GLYPH<0> bandlight curves of the source (Kochanek et al. 2017), which contain more than 3000 photometric epochs over a 7-year period, with a typical uncertainty of 0.02 mag. This did not reveal any significant periodic or long-term photometric variability: the optical source is constant at the & 0 GLYPH<149> 01 mag level.  \nGaiaBH2wasobservedby TESS (Ricker et al. 2015) during sectors 11 and 38. The source is sufficiently bright ( 𝑇 GLYPH<25> 11 GLYPH<149> 5) that it may be possible to measure asteroseismic parameters from the light curve (Stello et al. 2022). We extracted light curves from the TESS full frame images using eleanor (Feinstein et al. 2019) and analyzed their power spectra following the methods described by Stello et al. (2022); see Appendix C for details. We found a marginally significant detection of 𝜈 max = 61 𝜇 Hz in both sectors. If reliable, this value can be translated into a constraint on the giant's surface gravity using asteroiseismic scaling relations (e.g. Chaplin & Miglio 2013). Together with the spectroscopically-measured effective temperature, it implies a surface gravity log h 𝑔 GLYPH<157> GLYPH<16> cms GLYPH<0> 2 GLYPH<17> i GLYPH<25> 2 GLYPH<149> 67, consistent with the spectroscopic log 𝑔 (Section 3.8). Given the marginal nature of the detection, we did not include it in our subsequent analysis.", '3.2 Extinction': "We estimate the extinction to Gaia BH2 using the 3D dust map from Lallement et al. (2022), which predicts a modest integrated extinction of 𝐴 550 nm = 0 GLYPH<149> 58 mag at 𝑑 = 1 GLYPH<149> 16 kpc, corresponding to 𝐸 ' 𝐵 GLYPH<0> 𝑉 ' GLYPH<25> 0 GLYPH<149> 2 mag. The 𝐺 GLYPH<0> band extinction inferred by the Gaia GSP-Phot pipeline (Andrae et al. 2022) is 𝐴 𝐺 = 0 GLYPH<149> 70, corresponding to 𝐸 ' 𝐵 GLYPH<0> 𝑉 ' GLYPH<25> 0 GLYPH<149> 26, but we expect the value from the 3D dust map to be more reliable for stars in this evolutionary phase. We adopt 𝐸 ' 𝐵 GLYPH<0> 𝑉 ' = 0 GLYPH<149> 2 GLYPH<6> 0 GLYPH<149> 03 mag in our SED modeling.", '3.3 Ultraviolet observation': "Gaia BH2 does not have a published UV magnitude and is outside the published GALEX footprint. Measuring a UV magnitude or upper limit for the source is important for constraining possible flux contributions from a luminous secondary, so we obtained a 1625s observation of the source with the UVOT photometer (Roming et al.  \n2005) on board the Neil Gehrels Swift Observatory in September 2022 (ToO ID 17935). We used the UVM2 band, with an effective wavelength of 2246 Å, and also obtained simultaneous X-ray observations using XRT (Burrows et al. 2005).  \nWe analyzed the UVM2 image using the uvotsource routine distributed with Heasoft , using a 4 arcsec aperture centered on the Gaia source and an 8 arcsec background aperture in a nearby region with no obvious sources in the UVM2 image (ra = 207.558463, dec = -59.231955). Gaia BH2 is detected in the UVM2 image only with 2.7 𝜎 significance. A low-significance detection turns out to be sufficient for our purposes, since the main goal was to rule out a UV-bright companion. The AB magnitude in the UVM2 band is 22 GLYPH<149> 20 GLYPH<6> 0 GLYPH<149> 41 mag, corresponding to a flux density of ' 2 GLYPH<149> 9 GLYPH<6> 1 GLYPH<149> 1 ' GLYPH<2> 10 GLYPH<0> 17 erg s GLYPH<0> 1 cm GLYPH<0> 2 .  \nThe source was not detected in the simultaneous XRT observation. We did not analyze the XRT data further, since our Chandra observations (Section 3.10) provide significantly deeper X-ray constraints.", '3.4 SED model': "We constructed the source's broadband spectral energy distribution (SED) by combining the Swift UVM2 magnitude, synthetic SDSS 𝑢𝑔𝑟𝑖𝑧 photometry constructed from Gaia BP/RP spectra (Gaia Collaboration et al. 2022b), 2MASS 𝐽𝐻𝐾 photometry (Skrutskie et al. 2006), and WISE 𝑊 1 𝑊 2 𝑊 3 photometry (Wright et al. 2010). We set an uncertainty floor of 0.02 mag in all bands to allow for photometric calibration issues and imperfect models. We then fit the SED with single-star models to infer the temperature and radius of the red giant. We predict bandpass-integrated magnitudes using empirically calibrated model spectral from the BaSeL library (v2.2; Lejeune et al. 1997, 1998). We assume a Cardelli et al. (1989) extinction law with 𝑅 𝑉 = 3 GLYPH<149> 1 and adopt a prior on the reddening 𝐸 ' 𝐵 GLYPH<0> 𝑉 ' = 0 GLYPH<149> 20 GLYPH<6> 0 GLYPH<149> 03. We use pystellibs 1 to interpolate between model SEDs, and pyphot 2 to calculate synthetic photometry. We then fit the SED using emcee (Foreman-Mackey et al. 2013) to sample from the posterior, with the temperature, radius, metallicity, and reddening sampled as free parameters.  \nFigure 2. Spectral energy distribution of the Gaia BH2 system. In all panels, red points show observed photometry. Black line and open squares show the continuous SED and model photometry for the best-fit single-star model, whose parameters are listed in the top panel. 2nd panel shows the SED of an 8 GLYPH<149> 9 𝑀 GLYPH<12> main-sequence star; this would be 200 times more luminous than the observed source, and is thus ruled out. In the third and fourth panels, we consider a scenario where the Gaia astrometric inclination constraint (which implies 𝑀 2 GLYPH<24> 8 GLYPH<149> 9 𝑀 GLYPH<12> ) is spurious. In this case, the minimum plausible companion mass is 𝑀 2 GLYPH<24> 2 GLYPH<149> 6 𝑀 GLYPH<12> , corresponding to an edge-on orbit (Figure 5). Main-sequence companions of this mass would still be brighter than the observed source in the UV, both for the case of a single 2 GLYPH<149> 6 𝑀 GLYPH<12> companion (3rd panel) and for a 1 GLYPH<149> 3 , 1 GLYPH<149> 3 𝑀 GLYPH<12> inner binary (4th panel). These scenarios would also conflict with limits on the companion flux from the optical spectra (Section 3.8.1) which are shown with a blue upper limit.  \n<!-- image -->  \nThe results are shown in the top panel of Figure 2. A single-star model yields a reasonably good fit to the data, with 𝜒 2 GLYPH<157> 𝑁 data = 1 GLYPH<149> 16, where 𝑁 data is the number of photometric points. The best-fit SED model has 𝑇 eff = 4615 K and 𝑅 = 7 GLYPH<149> 77 𝑅 GLYPH<12> , with a bolometric luminosity of 𝐿 GLYPH<25> 25 𝐿 GLYPH<12> . Modeling of the source's high-resolution spectrum (Section 3.8) leads to a very similar constraint on 𝑇 eff . The high-precision of the Gaia parallax ( 𝜛 GLYPH<157> 𝜎 𝜛 = 47) allows us to constrain the star's radius within GLYPH<25> 3% and the luminosity within GLYPH<25> 6%.  \nThe lack of significant UV excess leads to stringent limits on any possible main-sequence companions: if the companion were a normal star, it would be impossible to obtain a good fit to the SED with a single-star model (e.g. El-Badry & Rix 2022). Luminous stars that are massive enough to explain the observed astrometry and RVs are predicted to outshine the giant in the UV and contribute significantly to the optical spectrum. These scenarios are explored in the lower panels of Figure 2 and discussed further in Section 4.", '3.5 Luminous star mass and evolutionary state': "We compare the measured radius and effective temperature of the luminous star to MIST single-star evolutionary models (Choi et al. 2016). We use solar-scaled models with 𝑍 = 0 GLYPH<149> 012 to account approximately for the star's measured sub-solar iron abundance and enhancement of 𝛼 GLYPH<0> elements (see Fu et al. 2018; Gaia BH2 has » Fe GLYPH<157> H … = GLYPH<0> 0 GLYPH<149> 22 and » 𝛼 GLYPH<157> Fe … = , 0 GLYPH<149> 26, as described in Section 3.8). These models are shown in the right panel of Figure 1.  \nFitting the measured temperature and radius implies a luminous star mass of 1 GLYPH<149> 07 GLYPH<6> 0 GLYPH<149> 19 𝑀 GLYPH<12> . The uncertainty is dominated by uncertainty in the observed effective temperature, though uncertainties in the stellar models likely contribute at a similar level (e.g. Joyce et al. 2022). The corresponding age range is about 5-13 Gyr. In all plausible models, the star is on the lower red giant branch and has recently begun to expand following core hydrogen exhaustion; it is not a core helium burning red clump star that has already reached the tip of the giant branch and experienced a helium flash.", '3.5.1 Could the star be a low-mass stripped giant?': "We considered whether the giant might be a product of recent envelope stripping, in which case it might have significantly lower mass than implied by single-star evolutionary models (e.g. El-Badry et al. 2022; El-Badry & Rix 2022). We consider such a scenario unlikely given the binary's wide and non-circular orbit, and the lack of evidence for ongoing mass transfer. In the current orbit, the star would overflow its Roche lobe at periastron only if it had a radius 𝑅 & 100 𝑅 GLYPH<12> . This essentially precludes a low-mass ( 𝑀 . 0 GLYPH<149> 5 𝑀 GLYPH<12> ) stripped giant, because a giant of that mass would not have reached a sufficiently large radius to overflow its Roche lobe in the Gaia BH2 system (e.g. Rappaport et al. 1995). A past period of mass transfer from the giant to the companion would also be expected to have circularized the orbit. Our analysis of the star's spectrum reveals no excess absorption or emission in H 𝛼 (Section 3.8), as has been commonly found in giants with ongoing or recently-terminated mass transfer (Jayasinghe et al. 2021, 2022). Finally, the spectroscopic log 𝑔 (Section 3.8) and the radius constraint from the SED fit imply a giant mass of 𝑀 ★ = 1 GLYPH<149> 17 , 0 GLYPH<149> 88 GLYPH<0> 0 GLYPH<149> 50 𝑀 GLYPH<12> ; i.e., a 1 𝜎 lower limit of 0 GLYPH<149> 67 𝑀 GLYPH<12> .", '3.6 Spectroscopy': 'To measure atmospheric parameters, abundances, and radial velocities, we obtained multi-epoch spectra of Gaia BH2 using two instruments.', '3.6.1 FEROS': "We observed Gaia BH2 43 times with the Fiberfed Extended Range Optical Spectrograph (FEROS; Kaufer et al. 1999) on the 2.2m ESO/MPG telescope at La Silla Observatory (programs P109.A-9001 and P110.9014). The first several observations used 2 GLYPH<2> 2 binning to reduce readout noise at the expense of spectral resolution; the remainder used 1 GLYPH<2> 1 binning. The resulting spectra have resolution 𝑅 GLYPH<25> 40 GLYPH<150> 000 (2 GLYPH<2> 2 binning) and 𝑅 GLYPH<25> 50 GLYPH<150> 000 (1 GLYPH<2> 1 binning). Most of our observations used 1200s or 900s exposures. The typical SNR at 5800 Å is 15-30.  \nWe reduced the data using the CERES pipeline (Brahm et al. 2017), which performs bias-subtraction, flat fielding, wavelength calibration, and optimal extraction. The pipeline measures and corrects for small shifts in the wavelength solution during the course a night via simultaneous observations of a ThAr lamp with a second fiber. Several minor modifications to the pipeline were necessary to obtain good performance with 2 GLYPH<2> 2 binning. We first calculate RVs by cross-correlating a template spectrum with each order individually and then report the mean RV across 15 orders with wavelengths between 4500 and 6700 Å. We calculate the uncertainty on this mean RV from the dispersion between orders; i.e., 𝜎 RV GLYPH<25> std ' RVs ' GLYPH<157> p 15.", '3.6.2 UVES': 'We observed Gaia BH2 5 times with the Ultraviolet-Visible Echelle Spectrograph (UVES; Dekker et al. 2000) mounted on the 8.2m UT2 telescope at the VLT on Cerro Paranal (program 2110.D-5024). We used Dichroic 2 with central wavelengths of 437 and 760 nm on the blue and red arms, providing spectral coverage of the wavelength ranges from 3730-4990 and 5670-9450 Å. We took 600s exposures and used 0.4 and 0.3 arcsec slits on the blue and red arms, respectively, which yielded resolutions of 80,000 and 110,000 with SNR GLYPH<24> 50 at 6500 Å.  \nWereducedthedatausing the ESO Reflex pipeline (Freudling et al. 2013) with standard calibrations. This performs bias-subtraction, flat fielding, wavelength calibration using afternoon ThAr arcs, optimal extraction, and order merging. We verified the stability of the wavelength solution using telluric absorption lines, measured RVs via cross-correlation of the merged red-arm spectra with a template, and adopted a per-epoch RV uncertainty of 0.1 km s GLYPH<0> 1 to account for drifts in the wavelength solution between science exposures and afternoon arcs. Near-simultaneous FEROS and UVES observations allowed us to verify good agreement between the two RV scales at the GLYPH<25> 0 GLYPH<149> 1 km s GLYPH<0> 1 level.', '3.6.3 Radial velocities': 'Our follow-up RVs are shown in Figure 3 and listed in Table 1. In total, we obtained 48 RVs between August 2022 and March 2023. All our observations have a precision better than 0.1 km s GLYPH<0> 1 . The precision of the FEROS RVs is somewhat higher than that of the UVESdata because of the precise wavelength calibration enabled by simultaneous arcs.  \nTable 1. Radial velocities.', '3.7 Orbit modeling': "We jointly fit our follow-up RVs and the astrometry+RV constraints from Gaia with a model that has 14 free parameters: the five standard astrometric parameters ( 𝛼GLYPH<150> 𝛿GLYPH<150> 𝜇 GLYPH<3> 𝛼 GLYPH<150> 𝜇 𝛿 GLYPH<150> and 𝜛 ); and the binary parameters period, eccentricity, inclination, angle of the ascending node Ω , argument of periastron 𝜔 , periastron time, center-of-mass velocity, luminous star mass 𝑀 ★ , and companion mass 𝑀 2 . For each call to the likelihood function, we predict the semi-major axis in physical units,  \n𝑎 = h 𝑃 2 orb 𝐺 ' 𝑀 ★ , 𝑀 2 ' GLYPH<157> GLYPH<16> 4 𝜋 2 GLYPH<17> i 1 GLYPH<157> 3 . We then calculate the angular semi-major axis of the photocenter ( 𝑎 0 ), and that of the primary ( 𝑎 1 ;  \nTable 2. Physical parameters and 1 𝜎 uncertainties for both components of Gaia BH2. We compare constraints on the orbit based on both Gaia and our RVs (3rd block), Gaia alone (4th block), and our RVs alone (5th) block. The 'our RVs alone' constraints adopt a prior of 𝑃 orb GLYPH<159> 2000 days. Flat and broad priors are used for all other constraints.  \ni.e., the star whose RVs are being measured):  \n𝑎 1 = 𝑎 𝑑 GLYPH<18> 𝑞 1 , 𝑞 GLYPH<19> (1)  \n𝑎 0 = 𝑎 𝑑 GLYPH<18> 𝑞 1 , 𝑞 GLYPH<0> 𝜖 1 , 𝜖 GLYPH<19> GLYPH<149> (2)  \nHere, 𝑞 = 𝑀 2 GLYPH<157> 𝑀 ★ is the mass ratio and 𝜖 is the 𝐺 GLYPH<0> band flux ratio of the companion to the primary. 𝑑 is the distance, which we compute as 𝑑 » kpc … = 1 GLYPH<157> 𝜛 » mas … . For a dark companion ( 𝜖 = 0), 𝑎 1 = 𝑎 0 . We then predict the six Thiele-Innes parameters:  \n𝐴 = 𝑎 0 ' cos 𝜔 cos Ω GLYPH<0> sin 𝜔 sin Ω cos 𝑖 ' (3)  \n𝐵  \n=  \n𝑎  \n0  \n'  \ncos  \n𝜔  \nsin  \nΩ  \n,  \nsin  \n𝜔  \ncos  \nΩ  \ncos  \n𝑖  \n'  \n(4)  \n𝐹 = GLYPH<0> 𝑎 0 ' sin 𝜔 cos Ω , cos 𝜔 sin Ω cos 𝑖 ' (5)  \n𝐺 = GLYPH<0> 𝑎 0 ' sin 𝜔 sin Ω GLYPH<0> cos 𝜔 cos Ω cos 𝑖 ' (6)  \n𝐶 = 𝑎 1 𝜛 sin 𝜔 sin 𝑖 (7)  \n𝐻 = 𝑎 1 𝜛 cos 𝜔 sin 𝑖GLYPH<149> (8)  \nNote that in the convention used here and in the Gaia archive, the astrometric parameters 𝐴GLYPH<150> 𝐵GLYPH<150> 𝐹 , and 𝐺 have angular units (mas), while the spectroscopic parameters 𝐶 and 𝐻 have physical units (au). In the Gaia data processing, 𝐶 and 𝐻 are constrained by the measured RVs, while 𝐴GLYPH<150> 𝐵GLYPH<150> 𝐹 , and 𝐺 are constrained by the astrometry (e.g. Pourbaix et al. 2022).  \nFor each call to the likelihood function, we construct the predicted vector of Gaia -constrained parameters:  \n𝜽 Gaia = GLYPH<2> 𝛼GLYPH<150> 𝛿GLYPH<150> 𝜛GLYPH<150> 𝜇 GLYPH<3> 𝛼 GLYPH<150> 𝜇 𝛿 GLYPH<150> 𝐴GLYPH<150> 𝐵GLYPH<150> 𝐹GLYPH<150> 𝐺GLYPH<150> 𝐶GLYPH<150> 𝐻GLYPH<150> 𝛾GLYPH<150> 𝑒GLYPH<150> 𝑃 orb GLYPH<150> 𝑇 𝑝 GLYPH<3> GLYPH<149> (9)  \nWe then calculate a likelihood term that quantifies the difference between these quantities and the Gaia constraints:  \nln 𝐿 Gaia = GLYPH<0> 1 2 GLYPH<0> 𝜽 Gaia GLYPH<0> 𝝁 Gaia GLYPH<1> | 𝚺 GLYPH<0> 1 Gaia GLYPH<0> 𝜽 Gaia GLYPH<0> 𝝁 Gaia GLYPH<1> GLYPH<149> (10)  \nHere 𝝁 Gaia and 𝚺 Gaia represent the vector of best-fit parameters constrained by Gaia and their covariance matrix, which we construct from the corr\\_vec parameter reported in the Gaia archive.  \nWeadditionally predict the RVs of the luminous star at the array of times 𝑡 𝑖 at which we obtained spectra. This leads to a radial velocity term in the likelihood,  \nln 𝐿 RVs = GLYPH<0> 1 2 ∑︁ 𝑖 GLYPH<16> RV pred ' 𝑡 𝑖 ' GLYPH<0> RV 𝑖 GLYPH<17> 2 𝜎 2 RV GLYPH<150> 𝑖 GLYPH<150> (11)  \nwhere RV 𝑖 and RV pred ' 𝑡 𝑖 ' are the measured and predicted RVs, with their uncertainties 𝜎 RV GLYPH<150> i . The predicted RVs are calculated by solving Kepler's equation using standard methods. The full likelihood is then given by  \nln 𝐿 = ln 𝐿 Gaia , ln 𝐿 RVs GLYPH<149> (12)  \nWe use flat priors on all parameters except 𝑀 ★ , for which we use a normal distribution informed by isochrones, 𝑀 ★ GLYPH<157> 𝑀 GLYPH<12> GLYPH<24> N' 1 GLYPH<149> 07 GLYPH<150> 0 GLYPH<149> 19 ' . We sample from the posterior using emcee (ForemanMackey et al. 2013) with 64 walkers, taking 3000 steps per walker after a burn-in period of 3000 steps. The results of the fit are summarized in Table 2 (3rd block).  \nWe initially considered including the flux ratio 𝜖 as a free parameter. But in this case we found the posterior for 𝜖 to pile up against 0, with a 1 𝜎 upper limit of 𝜖 GLYPH<159> 0 GLYPH<149> 11. That is, a luminous companion is  \ndisfavored by the astrometry because it would produce too small of a photocenter semi-major axis, independent of any arguments about the expected mass-luminosity relation of plausible luminous secondaries. In subsequent analysis, we therefore fixed 𝜖 = 0 in order to avoid introducing a bias due to unrealistically large flux contributions from the secondary, which we suspect is truly dark.", '3.7.1 Gaia-only constraints': "In order to explore the relative importance of our follow-up RVs and the Gaia data for our final constraints, we also repeat the fit without the follow-up RVs; i.e., simply removing the ln 𝐿 RVs term from Equation 12. We list the resulting constraints in Table 2 (4th block).  \nFigure 3 (top panel) compares the measured RVs to the predictions of the Gaia -only solution. There is significant ( GLYPH<6> 3 months) uncertainty in the predicted RV turnover time, presumably because the Gaia data did not cover a full period and were obtained two orbital cycles before our follow-up (see Appendix A). The follow-up RVs are fully consistent with the Gaia -only solution.  \nPerhaps counterintuitively, the best-fit period obtained from the Gaia -only fit is 𝑃 orb = 1300 GLYPH<6> 26 days, which is different at the GLYPH<24> 1 𝜎 level from the period constraint in the Gaia archive, 𝑃 orb = 1352 GLYPH<6> 45 days. Several other parameters are also slightly discrepant. The reason for this difference is that our modeling explicitly links the spectroscopic and astrometric Thiele-Innes parameters via the requirement that 𝑎 0 = 𝑎 1 , and by the requirement that a single pair of ( 𝜔GLYPH<150> 𝑖 ) can explain all six Thiele-Innes parameters. In contrast, the Gaia data processing only requires that the spectroscopic and astrometric orbits have the same period, eccentricity, and eccentric anomaly at fixed time; 𝐴GLYPH<150> 𝐵GLYPH<150> 𝐹GLYPH<150> 𝐺GLYPH<150> 𝐶 , and 𝐻 are fit to the astrometric and RV time series without enforcing constraints between them (see Pourbaix et al. 2022). Inspection of Equations 3-8 reveals that the six Thiele-Innes parameters are functions of four Campbell elements ' 𝑎 0 GLYPH<150> 𝜔GLYPH<150> Ω GLYPH<150> 𝑖 ), or of five if 𝑎 1 is fit as a separate parameter. This makes the problem over-constrained. In the presence of noise, there is thus no guarantee that a single set of orbital parameters can actually match all the Thiele-Innes parameters, even for a perfectly well-behaved orbit. Our approach minimizes the difference between the predicted and measured parameters. As we show in Table 3, there does indeed exist a single set of orbital parameters that can reproduce all the Gaia constraints and our follow-up RVs.", '3.7.2 Joint constraints': "Themiddle and bottom panels of Figure 3 compare the measured RVs to predictions of the joint Gaia +RVfit. With the RVs included, the fit includes data taken over the course of GLYPH<24> 3 orbital periods, resulting in a tightly constrained orbit. The follow-up RVs cover most of the best-fit orbit's dynamic range in RV and thus constrain the RV mass function more tightly than the pureGaia constraints (Section 3.7.3). Weshowtheresiduals of the measured RVs with respect to the best-fit orbit in the top panel of Figure 5. The fit is good, with 𝜒 2 GLYPH<157> 𝑁 RVs GLYPH<25> 0 GLYPH<149> 7, suggesting that the typical RV uncertainty is slightly overestimated.  \nIn Figure 4, we compare the constraints from joint modeling to those based only on the Gaia constraints. Reassuringly, the two sets of constraints are consistent, but the constraints that include our RVs are tighter.  \nFigure 3. Radial velocities (RVs). In the top panel, we compare predictions of the Gaia orbital solution (cyan lines) to our measured RVs. There is uncertainty in these predictions, which are based on Gaia observations taken GLYPH<24> 6 years earlier, but the RVs are consistent with them. In the middle and bottom panels, we show predictions from a joint fit of the RVs and Gaia constraints. Insets highlight periods with high-cadence data. The RVs have a typical uncertainty of 30 m s GLYPH<0> 1 and span most of the orbit's predicted RV range, so they tightly constrain the orbit when combined with the longer-baseline Gaia constraints.  \n<!-- image -->  \n-  \n-", '3.7.3 Companion mass and inclination': "Our joint fit constrains the mass of the unseen companion to 𝑀 2 = 8 GLYPH<149> 9 GLYPH<6> 0 GLYPH<149> 3 𝑀 GLYPH<12> . The dominant source of uncertainty in this estimate is the uncertainty in the red giant mass.  \nThe astrometric constraint on the eccentricity plays a critical role in constraining the mass, because the inclination, 𝑖 GLYPH<25> 35 deg, is relatively face-on. This means that the component of the star's orbital velocity projected onto our line-of-sight is only sin ' 35 deg ' GLYPH<25> 57% of the total orbital velocity. We show how the implied companion mass varies with inclination in Figure 5, where we plot the implied companion mass as a function of assumed inclination, giving the RV  \nmass function, 𝑓 ' 𝑀 2 ' = 𝑃 orb 𝐾 3 ★ GLYPH<16> 1 GLYPH<0> 𝑒 2 GLYPH<17> 3 GLYPH<157> 2 GLYPH<157> ' 2 𝜋𝐺 ' GLYPH<25> 1 GLYPH<149> 33 𝑀 GLYPH<12> . If  \nthe orbit were edge-on, the companion mass implied by the period, RV semi-amplitude, and eccentricity would be 𝑀 2 GLYPH<25> 2 GLYPH<149> 6 𝑀 GLYPH<12> . While this is much lower than the fiducial mass constraint when the inclination measurement is taken into account, it is still too massive for the unseen companion to be any plausible luminous star, whether a single star or a binary itself (Section 4).  \nOnly 18% of randomly oriented orbits are expected to have an inclination as low as 35 degrees. We note that in contrast to RV surveys, astrometric binary surveys have a selection bias in favor of low inclinations (e.g. Gaia Collaboration et al. 2022a). The good agreement between the RV prediction from the Gaia solution and our measured RVs, as well as the fact that RVs alone rule out a stellar companion even for an edge-on orbit (Figure 2), suggests that the Gaia inclination constraint is reliable.  \n0.25  \n/circledot  \n/circledot  \nFigure 4. Comparison between constraints obtained from joint fitting of our RVs and the Gaia data (black) and constraints from Gaia alone (cyan). The two sets of constraints are generally consistent, but those conditioned on the RVs are tighter. Numbers on the diagonals reflect the joint constraints. Both sets of constraints are listed in Table 2.  \n<!-- image -->  \nFigure 5. Top: comparison between the best-fit Gaia + RV solution and the measured RVs. The median absolute residual is 0.02 km s GLYPH<0> 1 , consistent with the RV uncertainties. Bottom: mass of the unseen companion implied by the RV orbital solution as a function of the assumed inclination. Shaded region shows the GLYPH<6> 1 𝜎 mass range for the luminous star. Hatched vertical band shows the Gaia inclination constraint, which implies a companion mass of 𝑀 2 GLYPH<25> 8 GLYPH<149> 9 GLYPH<6> 0 GLYPH<149> 3 𝑀 GLYPH<12> . If the orbit were edge-on (in conflict with the astrometric solution), the minimum companion mass would be GLYPH<25> 2 GLYPH<149> 6 𝑀 GLYPH<12> .  \n<!-- image -->", '3.7.4 Pure RV constraints': "Recognizing that the Gaia orbital solution could be wrong, we now explore constraints that can be obtained from our follow-up RVs alone.  \nWe fit the FEROS and UVES RVs with the Joker (PriceWhelan et al. 2017), which uses rejection sampling to obtain unbiased posterior samples consistent with a set of RVs. We used a 𝑝 ' 𝑃 orb ' / 1 GLYPH<157> 𝑃 orb prior between 100 and 2000 days, and broad uniform priors on other orbital parameters. Even with 10 10 prior samples, the Joker returned only a few samples that were marginally consistent with the RVs, and all of these samples were concentrated in a single posterior mode. This reflects the fact that the RVs strongly constrain the orbit, such that only a tiny fraction of orbital parameter space is consistent with the data. To better sample the posterior, we initialized standard MCMC chains near the maximum-likelihood sample and drew more samples using emcee . The results are reported in the bottom block of Table 2 and are shown in Figure 6. In the top panel, gray lines show RV predictions from random posterior draws; red line shows the prediction of the best-fit joint Gaia + RV solution. In the bottom panels, contours show constraints from the RVs alone, while the red square marks the best-fit joint solution.  \nThe RV-only constraints are fully consistent with those from the joint fit. Because our RVs cover only a small fraction of the period, the orbit is not fully constrained by RVs alone. In particular, the RVs allow for orbits with arbitrarily long periods combined with increas-  \nFigure 6. Constraints on the Gaia BH2 orbit from RVs alone, ignoring the Gaia constraints entirely. Black lines show predicted RV curves (top) and posterior constraints (bottom) from fitting our follow-up RVs. Red lines show the best-fit joint Gaia + RV solution. The two sets of constraints are fully consistent. Although they cover only a small fraction of the orbit, the RVs alone rule out periods below GLYPH<24> 1150 days and RV mass functions below GLYPH<24> 1 GLYPH<149> 25 𝑀 GLYPH<12> . The RVs alone do not provide a meaningful upper limit on the orbital period (here we adopted a prior with an upper limit of 𝑃 orb GLYPH<159> 2000 days), but any longer period would imply an even larger companion mass function.  \n<!-- image -->  \n/circledot  \ningly large eccentricities. For this reason, there are many posterior samples with periods near our adopted upper limit of 2000 days. However, with a prior of 𝑃 orb GLYPH<159> 2000 days, the orbit becomes wellconstrained. Most importantly, the RVs alone provide a well-defined minimum period and mass function, ruling out any 𝑃 orb . 1150 days and any RV mass function 𝑓 ' 𝑀 2 ' RVs . 1 GLYPH<149> 25 𝑀 GLYPH<12> . This means that independent of the Gaia data - the RV orbit and luminous star mass imply a companion mass of at least 2 GLYPH<149> 5 𝑀 GLYPH<12> .", '3.8 Spectral analysis': "We analyzed the UVES spectra to measure atmospheric parameters and abundances of the red giant and place further constraints on the light contributions of the secondary, using a combination of abinitio and empirical modeling. Most of our analysis focused on the spectrum obtained on JD 2459935, but we verified that there are no significant differences between the rest-frame UVES spectra taken on different nights.  \nTable 3. Comparison of the Gaia constraints on Gaia BH2's orbit to constraints from the joint RVs+astrometry fit. All parameters are consistent at the GLYPH<24> 2 𝜎 level. The most discrepant parameter is h\\_thiele\\_innes , for which there is 2.2 𝜎 tension. Overall, this comparison suggests that the Gaia solution is reliable, though its uncertainties could be underestimated slightly.  \nWe first fit the UVES spectra using the Brussels Automatic Code for Characterizing High accUracy Spectra (BACCHUS; Masseron et al. 2016) with the same set up as in El-Badry et al. (2023). BACCHUS enables us to derive the stellar atmospheric parameters, including the effective temperature ( 𝑇 eff ), surface gravity (log 𝑔 ), metallicity ([Fe/H]) and microturblent velocity ( 𝑣 micro ) by assuming Fe excitation/ionization balance; i.e., the requirement that lines with different excitation potentials all imply the same abundances. We use the fifth version of the Gaia-ESO atomic linelist (Heiter et al. 2021). Hyperfine structure splitting is included for Sc I, V I Mn I, Co I, Cu I, Ba II, Eu II, La II, Pr II, Nd II, Sm II (see more details in Heiter et al. 2021). We also include molecular line lists for the following species: CH (Masseron et al. 2014), and CN, NH, OH, MgH and C 2 (T. Masseron, private communication). Finally, we also include the SiH molecular line list from the Kurucz linelists 3 . Spectral synthesis for BACCHUS is done using the TURBOSPECTRUM (Alvarez & Plez 1998; Plez 2012) code along with the line lists listed above and the MARCS model atmosphere grid (Gustafsson et al. 2008).  \nOnce the stellar parameters were determined, the model atmosphere is fixed and individual abundances were derived using BACCHUS' 'abund' module. For each spectral absorption feature, this module creates a set of synthetic spectra, which range between -0.6 GLYPH<159> [X/Fe] GLYPH<159> +0.6 dex, and performs a 𝜒 2 minimization between the observed and synthetic spectra. The reported atmospheric [X/Fe] abundances are the median of derived [X/Fe] across all lines for a given species. The uncertainty in the atmospheric [X/Fe] is defined as the dispersion of the [X/H] abundance measured across all lines for a species. If only 1 absorption line is used, we conservatively assume a [X/Fe] uncertainty of 0.10 dex. For a more detailed discussion of the BACCHUS code, we refer the reader to Section 3 of both Hawkins et al. (2020a,b). We ran BACCHUS on the UVES spectrum after merging individual de-blazed orders and performing a preliminary continuum normalization using a polynomial spline. Further normalization is performed by BACCHUS. Spectral regions affected by telluric absorption were masked during fitting, and we only fit the portion of the spectrum with 𝜆 GLYPH<159> 7600 Å, where tellurics are less severe.  \nThe resulting stellar parameters and abundances are listed in Table 4. Figure 7 compares portions of the observed spectrum to the best-fit BACCHUS model. The fit is overall quite good: while there are some deviations between the model and data, the quality of the fit is typical for ab-initio (as opposed to data-driven) fits to high-resolution spectra. The remaining differences between data and  \nmodel can be attributed to a combination of imperfect continuum normalization and imperfections in the line list.  \nThe most important conclusion from this analysis is that the star's iron abundance is subsolar ( » Fe GLYPH<157> H … = GLYPH<0> 0 GLYPH<149> 22), but the abundances of 𝛼 GLYPH<0> elements (i.e.; O, Mg, Si, Ca, and Ti) are enhanced relative to solar, with » 𝛼 GLYPH<157> Fe … = 0 GLYPH<149> 26 GLYPH<6> 0 GLYPH<149> 05. This abundance pattern is characteristic of old stars in the Galactic thick disk or metal-rich halo, which formed from material that was enriched by type II SNe, with less enhancement of iron-peak elements by SNe Ia than typical stars in the thin disk. If Gaia BH2 formed in the thick disk, this would imply an age of 8-12 Gyr (e.g Nissen et al. 2020; Xiang & Rix 2022) and (assuming mass transfer was negligible) a stellar mass of 𝑀 ★ GLYPH<25> 1 GLYPH<149> 0 GLYPH<6> 0 GLYPH<149> 07 𝑀 GLYPH<12> for the giant.  \nBeyond the strong 𝛼 GLYPH<0> enhancement, the abundances of the giant appear normal. We compared the star's position in the [X/Fe] vs. [Fe/H] plane to the population of stars in the solar neighborhood with abundances measured by Adibekyan et al. (2012), Bensby et al. (2014), and Battistini & Bensby (2016), finding that all elements fall within the GLYPH<24> 2 𝜎 scatter of the observed local population. Lithium is not detected, as expected for a star on the giant branch that has undergone first dredge-up. There is no evidence of enhancement in neutron-capture elements relative to typical thick-disk stars.", '3.8.1 Empirical comparison to other red giant spectra and constraints on luminous companions': 'To search for possible anomalous spectral features - which could arise, for example, due to a luminous companion, emission from an accretion disk, or absorption by an accretion stream - we compared the spectrum of Gaia BH2 to spectra of stars with similar stellar parameters and abundances observed by the GALAH survey (De Silva et al. 2015; Buder et al. 2021). To this end, we degraded the observed UVES spectrum to 𝑅 = 28 GLYPH<150> 000, shifted it to rest frame, and identified its nearest neighbor in pixel space among stars observed by GALAH with SNR > 100. The closest match was Gaia DR3 source 5909513401713377280, which in GALAH DR3 has 𝑇 eff = 4591 GLYPH<6> 74 K, log 𝑔 = 2 GLYPH<149> 45 GLYPH<6> 0 GLYPH<149> 19, » Fe GLYPH<157> H … = GLYPH<0> 0 GLYPH<149> 36 GLYPH<6> 0 GLYPH<149> 05, and » 𝛼 GLYPH<157> Fe … = 0 GLYPH<149> 28 GLYPH<6> 0 GLYPH<149> 01. These parameters are all similar to those we find for Gaia BH2, providing independent validation of our BACCHUS-inferred parameters. In Figure 8, we compare the normalized and resolution-matched spectra of the two sources in the 1st and 3rd GALAH wavelength windows, which respectively contain H 𝛽 and H 𝛼 . The spectra are nearly indistinguishable.  \nThesimilarity of the two observed spectra, and the good agreement between the data and the BACCHUS model spectrum in Figure 7, speaks against the presence of a luminous secondary. The quantitative  \nTable 4. Parameters of the red giant inferred by BACCHUS. 𝑁 lines is the number of independent absorption lines used to infer each element. 𝛼 is an average of O, Mg, Si, Ca, and Ti.  \nupper limit depends on the spectral type of the putative companion: a rapidly rotating secondary would be harder to detect. However, the depth of the observed absorption lines rules out even a companion that contributes only continuum. The total flux in the cores of the deepest lines is only GLYPH<25> 2% of the continuum value (Figure 7). A companion contributing more than GLYPH<25> 2% of the total light at 40004300 Å would dilute the absorption lines beyond this value and is thus ruled out. We conservatively adopt a GLYPH<159> 3%fluxlimit on the total flux; this is the origin of the blue upper limit in Figure 2.', '3.9 Galactic orbit': 'Toinvestigate the Galactic orbit of Gaia BH2, we used its parallax and proper motion from the Gaia astrometric binary solution, as well as the center-of-mass RV inferred from the joint fit, as starting points to compute its orbit backward in time for 500 Myr using galpy (Bovy 2015). We used the Milky Way potential from McMillan (2017). The result is shown in Figure 9; for comparison, we also show the orbit of the Sun. The orbit is typical of a thin -disk star, with modest eccentricity and excursions above the disk midplane limited to GLYPH<6> 80 pc.', '3.9.1 Interpretation of the abundances and kinematics': "The combination of thin-disk kinematics with 𝛼 GLYPH<0> rich chemistry is unusual. This is illustrated in Figure 10, where we compare Gaia BH2's chemistry and kinematics to the solar neighborhood sample from Bensby et al. (2014). The left panel shows the Toomre diagram, while the right panel shows the abundance of titanium (a representative 𝛼 GLYPH<0> element) relative to iron. Most of the 𝛼 GLYPH<0> enhanced stars in the Bensby et al. (2014) sample have thick-disk kinematics. In contrast,  \nGaia BH2 is strongly enhanced in all 𝛼 elements but has a Galactic orbit typical of the thin disk.  \nGiven its thin-disk-like orbit, it is tempting to attribute the 𝛼 GLYPH<0> enhancement of Gaia BH2 to pollution from the BH progenitor. Indeed, similar 𝛼 GLYPH<0> enhancement in the donor stars of BH X-ray binaries has been interpreted as evidence for pollution during the supernova event. For example, 0.2-0.5 dex enhancement of some 𝛼 GLYPH<0> elements in the donor stars in GRO J1655-40 and V404 Cyg - the BH X-ray binaries whose detailed abundance patterns can be measured most readily - has been interpreted by Israelian et al. (1999) and González Hernández et al. (2011) as evidence that the BHs formed in a supernova. In those systems, the observed abundance patterns can plausibly be explained by pollution from a supernova with yields similar to those predicted for stars with initial masses of 25 GLYPH<0> 40 𝑀 GLYPH<12> .  \nIt is more challenging to apply a similar a similar explanation to Gaia BH2. At the periastron distance of GLYPH<25> 2 GLYPH<149> 4 au, a star with 𝑅 ★ GLYPH<25> 1 𝑅 GLYPH<12> subtends only GLYPH<25> 10 GLYPH<0> 6 of the sky as seen from the BH progenitor. The supernova ejecta are expected to consist mostly ( GLYPH<25> 70%) of oxygen (Thielemann et al. 1996) and must have had a total mass of . 10 𝑀 GLYPH<12> in order for the binary to remain bound. For a spherically symmetric explosion, we thus expect at most 7 GLYPH<2> 10 GLYPH<0> 6 𝑀 GLYPH<12> of oxygen to have been deposited on the companion during a supernova. During first dredge-up, this would have been mixed with the GLYPH<25> 0 GLYPH<149> 005 𝑀 GLYPH<12> of oxygen already present in the star's convective envelope, and would not significantly increase the surface oxygen abundance observable today. Similar calculations apply to the predicted yields of other 𝛼 GLYPH<0> elements. It thus seems unlikely that the observed 𝛼 GLYPH<0> enhancement is a consequence of pollution by supernova ejecta, if these ejecta escaped the BH progenitor at high velocity and were deposited on the secondary directly.  \nIt is possible, however, that some of the ejecta escaped the BH progenitor at low velocity, reached a separation of several au, and remained bound to the binary. In this case, a significant amount of the material may have eventually been accreted by the companion star, and it is possible that this supplied the GLYPH<25> 0 GLYPH<149> 005 𝑀 GLYPH<12> of 𝛼 GLYPH<0> elements required to explain the observed abundances. The expected total accretion onto the secondary in this scenario will depend on the energetics of the explosion and on the structure of the BH progenitor before its death; detailed calculations are required to explore this further.", '3.10 X-ray observations': 'We observed Gaia BH2 for 20ks using the Advanced CCD Imaging Spectrometer (ACIS; Garmire et al. 2003) on board the Chandra X-ray telescope on 2023 January 25 (proposal ID 23208881; PI: ElBadry). We used the ACIS-S configuration, with a spatial resolution of about 1 arcsec. The observations were timed to occur near the periastron passage, when the separation between the BH and the star was GLYPH<25> 2 GLYPH<149> 47 au.  \nThe Chandra data are shown in the upper left panel of Figure 11. There is no obvious X-ray source coincident with the optical source, whose position is marked with a red circle. We used the Chandra Interactive Analysis of Observations (Fruscione et al. 2006) software to place upper limits on the X-ray flux. We first ran the chandra\\_repro tool to reprocess the observation; this creates a new bad pixel file and de-streaks the event file. We performed aperture photometry at the optical position of Gaia BH2 in the reprocessed event file using the srcflux tool. We detect no significant flux at the location of Gaia BH2 and obtain a background count rate of 1 GLYPH<149> 45 GLYPH<2> 10 GLYPH<0> 4 counts s GLYPH<0> 1 . We assume a foreground hydrogen column density of 𝑁 𝐻 GLYPH<25> 1 GLYPH<149> 36 GLYPH<2> 10 21 cm GLYPH<0> 2 based on the optical extinction  \nFigure 7. UVES spectral cutouts of Gaia BH2 (red) compared to BACCHUS spectral model (black; parameters are in Table 4). Top panel shows a 20 Å-wide region containing the Ca I resonance line at 4227 Å, as well as many other metal lines. Bottom panels highlight lines from 5 different 𝛼 -elements, whose abundances are all enhanced relative to the solar abundance pattern. The deepest metal lines reach a depth of GLYPH<25> 2% of the continuum; this rules out luminous companions that contribute more than GLYPH<25> 2% of the light at these wavelengths.  \n<!-- image -->  \n(e.g. Güver & Özel 2009) and use the Portable, Interactive MultiMission Simulator (PIMMS) tool to calculate the corresponding upper limit on the unabsorbed source flux. Based on the background count rate and 𝑁 𝐻 , and assuming a power law spectrum with photon index 2, we obtain a 2 𝜎 flux limit in the 0.5-7 keV energy range of 𝐹 𝑋 GLYPH<159> 5 GLYPH<149> 1 GLYPH<2> 10 GLYPH<0> 15 erg s GLYPH<0> 1 cm GLYPH<0> 2 . This limit is shown as a horizontal dashed line in Figure 11.', '3.11 Radio observations': "WeobservedGaiaBH2for4hourswiththeMeerKATradiotelescope in L-band (1.28 GHz) on 2023 January 13 (DDT-20230103-YC-01, PI: Cendes), when the separation between the BH and the star was GLYPH<25> 2 GLYPH<149> 54 au. We used the flux calibrator J1939-6342 and the gain calibrator J1424-4913, and used the calibrated images obtained via the SARAO Science Data Processor (SDP) 4 for our analysis.  \nTheMeerKATdataareshowninthebottomleftpanel of Figure 11. There is no detectable radio source coincident with Gaia BH2. We measured the flux density using the imtool package within pwkit (Williams et al. 2017) at the location of Gaia BH2. The RMS at the source's position is 17 𝜇 Jy, so we report a 2 𝜎 (3 𝜎 ) nondetection of GLYPH<159> 34 𝜇 Jy ( GLYPH<159> 51 𝜇 Jy). The 2 𝜎 limit is shown with a horizontal dashed line in Figure 11.", '3.12.1 Bondi-Hoyle-Lyttleton accretion': "A rough estimate of the expected accretion rate onto the BH can be obtained under the assumption that a spherically-symmetric wind from the giant is accreted at the Bondi-Hoyle-Lyttleton (BHL) rate:  \n/ 𝑀 BHL = 𝐺 2 𝑀 2 BH / 𝑀 wind 𝑣 4 wind 𝑑 2 p = 2 GLYPH<2> 10 GLYPH<0> 13 𝑀 GLYPH<12> yr GLYPH<0> 1 GLYPH<18> 𝑀 BH 10 𝑀 GLYPH<12> GLYPH<19> 2 GLYPH<18> / 𝑀 wind 10 GLYPH<0> 11 𝑀 GLYPH<12> yr GLYPH<0> 1 GLYPH<19> GLYPH<2> GLYPH<18> 𝑣 wind 150 km s GLYPH<0> 1 GLYPH<19> GLYPH<0> 4 GLYPH<18> 𝑑 p 2 GLYPH<149> 5 au GLYPH<19> GLYPH<0> 2 GLYPH<149> (13)  \nThe giant's mass-loss rate, / 𝑀 wind , is uncertain. We consider the prescription from Reimers (1975):  \n/ 𝑀 wind = 4 GLYPH<2> 10 GLYPH<0> 13 𝛽 𝑅 𝑀 GLYPH<12> yr GLYPH<0> 1 GLYPH<18> 𝐿 ★ 𝐿 GLYPH<12> GLYPH<19> GLYPH<18> 𝑅 ★ 𝑅 GLYPH<12> GLYPH<19> GLYPH<18> 𝑀 ★ 𝑀 GLYPH<12> GLYPH<19> GLYPH<0> 1 GLYPH<150> (14)  \nwhere 𝛽 𝑅 is a dimensionless constant; we adopt 𝛽 𝑅 = 0 GLYPH<149> 1 following Choi et al. (2016). For the wind velocity, a reasonable approximation is that it scales with the escape velocity at the stellar surface,  \n𝑣 wind = 600 km s GLYPH<0> 1 𝛽 wind GLYPH<18> 𝑀 ★ 𝑀 GLYPH<12> GLYPH<19> 1 GLYPH<157> 2 GLYPH<18> 𝑅 ★ 𝑅 GLYPH<12> GLYPH<19> GLYPH<0> 1 GLYPH<157> 2 GLYPH<150> (15)  \nwhere 𝛽 wind is another dimensionless constant.  \nFigure 8. Spectral cutouts of Gaia BH2 (red) compared to a reference star (black), a thick-disk star with similar stellar parameters and abundances observed by the GALAH survey. The two spectra are very similar. This rules out significant light contributions from a companion and speaks against ongoing mass transfer. It also indicates that the abundance pattern of the giant is not grossly unusual.  \n<!-- image -->  \nFigure 9. Galactic orbit of Gaia BH2, calculated backward for 500 Myr from the measured proper motion and center-of-mass RV. For comparison, we show the Sun's orbit calculated over the same period. The orbit is typical of a thin-disk star. This suggests that the BH was not born with a large natal kick and also likely did not form in a classical globular cluster.  \n<!-- image -->  \nWe then assume that the X-ray flux is emitted isotropically with a  \nradiative efficiency 𝜂 𝑋 , such that 𝐹 X GLYPH<150> BHL = 𝜂 𝑋 / 𝑀 BHL 𝑐 2 4 𝜋𝑑 2 . This leads to a predicted X-ray flux for accretion at the BHL rate:  \n𝐹 X GLYPH<150> BHL = 10 GLYPH<0> 14 erg s GLYPH<0> 1 cm GLYPH<0> 2 GLYPH<18> 𝐿 ★ 25 𝐿 GLYPH<12> GLYPH<19> GLYPH<18> 𝑅 ★ 8 𝑅 GLYPH<12> GLYPH<19> 3 GLYPH<18> 𝑀 ★ 𝑀 GLYPH<12> GLYPH<19> GLYPH<0> 3 GLYPH<18> 𝑀 BH 9 𝑀 GLYPH<12> GLYPH<19> 2 GLYPH<2> GLYPH<18> 𝑑 p 2 GLYPH<149> 5 au GLYPH<19> GLYPH<0> 2 GLYPH<18> 𝑑 1 kpc GLYPH<19> GLYPH<0> 2 GLYPH<18> 𝜂 𝑋 10 GLYPH<0> 4 GLYPH<19> GLYPH<18> 𝛽 𝑅 0 GLYPH<149> 1 GLYPH<19> GLYPH<18> 𝛽 w 1 GLYPH<19> 4 GLYPH<150>  \n<!-- image -->  \nFigure 10. Toomre diagram (left) and [Ti/Fe] vs. [Fe/H] diagram (right; Ti is a represenative 𝛼 -element). Small points show stars in the solar neighborhood from Bensby et al. (2014) that are 𝛼 GLYPH<0> rich (red) and 𝛼 GLYPH<0> poor (black). Large red and cyan stars show Gaia BH2 and BH1, respectively. Gaia BH2 is kinematically part of the thin disk. Its 𝛼 GLYPH<0> enhanced abundances are, however, more characteristic of the thick disk. The opposite is true for Gaia BH1.  \n<!-- image -->  \n(16)  \nwhere we have scaled to parameters similar to Gaia BH2. Given the manyapproximationsthatenterthis scaling, it is only expected to hold to order of magnitude. The most important prediction of Equation 16 is the strong scaling with the radius and luminosity of the luminous star, 𝐹 X GLYPH<150> BHL / 𝑅 3 ★ 𝐿 ★ . This means that at fixed mass, a giant leads to a much larger predicted X-ray flux than a dwarf. Physically, this occurs because (a) the giant has a significantly higher mass loss rate, / 𝑀 wind , and (b), the wind from a giant is slower, allowing the BH to capture a larger fraction of it. As a result, the predicted X-ray flux from Gaia BH2 is GLYPH<25> 1 GLYPH<150> 000 times larger than that from Gaia BH1, even though Gaia BH2 is farther away and has a wider orbit.  \nWe can also predict the radio flux assuming the source falls on the canonical X-ray - radio flux correlation. We take the version of this relation from Gallo et al. (2006):  \nlog 𝐿 𝑅 = 0 GLYPH<149> 6 log 𝐿 𝑋 , 0 GLYPH<149> 8 log ' 𝑀 BH GLYPH<157> 𝑀 GLYPH<12> ' , 7 GLYPH<149> 3 GLYPH<150> (17)  \nwhere 𝐿 𝑅 and 𝐿 𝑋 represent the radio and X-ray luminosity. Following Gallo et al. (2006), we translate this into a radio flux density under the assumption of a flat radio spectrum ( 𝑆 𝜈 / 𝜈 0 ) up to a maximum frequency 𝜈 max = 8 GLYPH<149> 4 GHz.  \nWe show the expected X-ray and radio fluxes from Gaia BH2 in Figure 11 with red-dashed lines. For each value of / 𝑀 BHL , we predict the radiative efficiency using the fitting function from Xie & Yuan (2012), who used simulations of hot accretion flows to calculate these efficiencies as a function of the accretion rate at the event horizon in Eddington units, for several values of the electron heating parameter, 𝛿 . We consider their predictions for 𝛿 = 0 GLYPH<149> 5. For a plausible wind speed of 𝑣 wind GLYPH<25> 150 km s GLYPH<0> 1 (somewhat lower than the giant's escape velocity, since the wind must escape the giant's gravitational potential before reaching the BH), the predicted X-ray and radio fluxes are in tension with the observed non-detections.", '3.12.2 Reduction in / 𝑀 near the BH': "The calculations above implicitly assume that the accretion rate near the event horizon is equal to the BHL rate. This may not be realized  \nin practice if, for example, a significant fraction of the inflowing material is lost to winds (e.g. Blandford & Begelman 1999). Numerous simulations of radiatively inefficient accretion flows predict that / 𝑀 decreases toward the event horizon roughly as / 𝑀 / 𝑟 0 GLYPH<149> 5 , where 𝑟 is the distance from the BH (Pen et al. 2003; Pang et al. 2011; Yuan et al. 2012; Ressler et al. 2021). There are several proposed physical explanations for this scaling. For concreteness, we consider the model from Yuan et al. (2012), in which / 𝑀 / 𝑟 𝑠 , where 𝑠 GLYPH<25> 0 GLYPH<149> 5 down to 𝑟 GLYPH<24> 10 𝑅 𝑠 , and 𝑠 GLYPH<25> 0 at smaller radii, where 𝑅 𝑠 is the Schwarzschild radius.  \nThe BHL accretion rate corresponds to the accretion rate on scales of the BH's accretion radius, 𝑅 acc = 𝐺𝑀 BH GLYPH<157> 𝑣 2 wind GLYPH<25> 0 GLYPH<149> 3 au. The model described above thus predicts that the accretion rate at the event horizon (and all radii interior to GLYPH<24> 10 𝑅 𝑠 ) is  \n/ 𝑀 GLYPH<25> GLYPH<18> 10 𝑅 𝑠 𝑅 acc GLYPH<19> 0 GLYPH<149> 5 / 𝑀 BHL = p 20 𝑣 wind 𝑐 / 𝑀 BHL (18) GLYPH<25> 0 GLYPH<149> 002 / 𝑀 BHL GLYPH<150> (19)  \ni.e., a factor of GLYPH<24> 500 reduction in / 𝑀 relative to the BHL rate!  \nA lower / 𝑀 also leads to a lower predicted radiative efficiency. The expected accretion rate at the horizon is / 𝑀 GLYPH<25> 4 GLYPH<2> 10 GLYPH<0> 16 GLYPH<25> 2 GLYPH<2> 10 GLYPH<0> 9 / 𝑀 edd , and the corresponding radiative efficiency from Xie & Yuan (2012) is predicted to be 𝜂 𝑋 GLYPH<25> 7 GLYPH<149> 5 GLYPH<2> 10 GLYPH<0> 5 , as opposed to 𝜂 𝑋 GLYPH<25> 0 GLYPH<149> 004fortheBHLcase.Substitutingthis value into Equation 16 and finally multiplying by 0.002 (Equation 19), we obtain a predicted 𝐹 𝑋 GLYPH<25> 10 GLYPH<0> 17 erg s GLYPH<0> 1 cm GLYPH<0> 2 for Gaia BH2, which is more than 2 orders of magnitude below the Chandra limit. The corresponding radio flux is GLYPH<25> 0 GLYPH<149> 2 𝜇 Jy, more than 2 orders of magnitude below the MeerKAT limit. The predicted X-ray and radio fluxes are shown with solid black lines in Figure 11.  \nTo summarize: if we take the accretion rate onto the BH to be the BHL rate predicted by Equation 13, then the predicted radiative efficiency is 𝜂 𝑋 = 0 GLYPH<149> 004, and the X-ray flux predicted by Equation 16 is GLYPH<25> 4 GLYPH<2> 10 GLYPH<0> 13 erg s GLYPH<0> 1 cm GLYPH<0> 2 , which would have been easily detected. On the other hand, if we reduce the predicted accretion rate to account for a reduction in / 𝑀 between the accretion radius and the horizon (Equation 18), the predicted accretion rate and radiative efficiency both fall significantly, and the expected X-ray and radio fluxes would  \nnot be detectable. The nondetection of Gaia BH2 thus supports models in which winds, convection, or other processes cause a significant reduction in the horizon accretion rate relative to the BHL rate. Such models have also been supported by polarization measurements of emission from Sgr A* and M87*, which show that the accretion rate near the event horizon is 2-3 orders of magnitude less than the Bondi value (Quataert & Gruzinov 2000; Kuo et al. 2014).", '4 NATURE OF THE COMPANION': "Since we have not detected the companion directly, we can infer its nature only from (a) dynamical constraints on its mass and (b) the fact that it emits little if any detectable light. The combination of Gaia astrometry and our RVs yields a robust and precise constraint on its mass: 8 GLYPH<149> 9 GLYPH<6> 0 GLYPH<149> 3 𝑀 GLYPH<12> . This relies on astrometric constraints on the orbital inclination, which are not accessible with RVs alone.  \nIt is instructive to disregard the astrometric inclination constraint and consider only the RVs. The orbital solution from our joint fit has a spectroscopic mass function, 𝑓 ' 𝑀 2 ' = 𝑃 orb 𝐾 3 ★ GLYPH<16> 1 GLYPH<0> 𝑒 2 GLYPH<17> 3 GLYPH<157> 2 GLYPH<157> ' 2 𝜋𝐺 ' = 1 GLYPH<149> 33 GLYPH<6> 0 GLYPH<149> 01 𝑀 GLYPH<12> . This represents the absolute minimum companion mass compatible with the RV orbit, in the limit of an edge-on orbit and a luminous star with zero mass. Adopting 𝑀 ★ = 1 GLYPH<149> 07 GLYPH<6> 0 GLYPH<149> 19 𝑀 GLYPH<12> and assuming an edge-on orbit (in conflict with the astrometric solution), this implies an minimum companion mass of 2 GLYPH<149> 64 GLYPH<6> 0 GLYPH<149> 16 𝑀 GLYPH<12> (Figure 5).  \nWe now consider what kinds of companions could have the dynamically implied mass. As demonstrated in Figure 2, a 9 𝑀 GLYPH<12> mainsequence star would outshine the giant at all wavelengths. It would have a bolometric luminosity of 𝐿 2 & 5 GLYPH<150> 000 𝐿 GLYPH<12> , at least 200 times brighter than the giant. An evolved 9 𝑀 GLYPH<12> star would be even more luminous, so any 9 𝑀 GLYPH<12> star is ruled out by the non-detection of a second component in the UVES spectra and by the faintness of the source in the Swift UVM2 data. A 9 𝑀 GLYPH<12> companion is also too massive to be a white dwarf or neutron star. While a small cluster of white dwarfs or neutron stars could escape detection, such a configuration would be difficult to assemble and would rapidly become unstable. We conclude that if the Gaia astrometric solution is correct, the companion must be a BH, or a close binary containing a BH.  \nWhat if the astrometric constraint on the inclination is spurious? In that case - if the period and eccentricity from the Gaia solution are still reliable - the companion mass could be as low as GLYPH<25> 2 GLYPH<149> 6 𝑀 GLYPH<12> . This is still too massive to be a white dwarf or a neutron star, and a 2 GLYPH<149> 6 𝑀 GLYPH<12> main-sequence companion is inconsistent with the SED and spectra (3rd panel of Figure 2). We also consider the possibility that the companion is an unresolved binary consisting of two stars with total mass 2 GLYPH<149> 6 𝑀 GLYPH<12> . An equal-mass inner binary will yield the highest mass-luminosity ratio, and the bottom panel of Figure 2 shows that this scenario is also ruled out by the faintness of the source in the UV and the spectral limit in the blue optical.  \nThere is also an evolutionary problem for scenarios in which the companion contains one or more main-sequence stars: in all cases, the stars would be more massive than the giant itself, and should thus be more evolved than it is. But one or more red giant companions would outshine the giant in the optical.", '4.1 Could the Gaia orbit be spurious?': "Since our RV follow-up has covered only a limited fraction of Gaia BH2's orbit, it is in principle possible that the Gaia orbit is simply wrong. In practice, our follow-up RVs alone constrain the orbit  \nfairly precisely (Figure 6), excluding any orbit with 𝑃 orb . 1150 days or 𝑓 ' 𝑀 2 ' RVs GLYPH<159> 1 GLYPH<149> 25 𝑀 GLYPH<12> . While we do not have access to the epoch-level astrometric data on which the Gaia solution is based, it is very unlikely that the observed RVs would just happen to match the predictions of the Gaia solution if that solution were incorrect. The Gaia orbital solution also appears unproblematic based on its goodness of fit metrics (Appendix B), and the source was observed by Gaia almost 90 times over a 1000-day period, sampling most of the orbital phase (Appendix A). Gaia BH2 has a AstroSpectroSB1 solution. In order to obtain an AstroSpectroSB1 solution in Gaia DR3, a source must first have independent spectroscopic and astrometric solutions, which can be combined in a manner that improves the goodness of fit relative to the purely astrometric solution (Pourbaix et al. 2022). This means that if the solution were spurious, the spectroscopic and astrometric orbits would independently have to be spurious in compatible ways - a scenario that is unlikely to occur for a source with such good astrometric phase overage. And finally, the minimum companion mass implied by the RVs alone - 2 GLYPH<149> 5 𝑀 GLYPH<12> - is still too high to be a white dwarf, or a main-sequence star that could escape spectroscopic and photometric detection.  \nMore broadly, the fact that the measured RVs are in good agreementwiththepredictions of the Gaia solution, and the Gaia fit quality metrics are unremarkable, leaves us with little reason to doubt the reliability of the Gaia solution, including the inclination constraint. In this case, the inclination is robustly constrained to 34 GLYPH<149> 9 GLYPH<6> 0 GLYPH<149> 4 degrees, and the companion mass to 8 GLYPH<149> 9 GLYPH<6> 0 GLYPH<149> 3 𝑀 GLYPH<12> . A companion of this mass must contain at least one BH. A single BH with 𝑀 GLYPH<24> 8 GLYPH<149> 9 𝑀 GLYPH<12> thus seems to be the simplest explanation, although the data do not rule out an inner binary containing 2 BHs or a BH and another compact object or low-mass star.", '4.2 Comparison to other known BHs': "In Figures 12-14, we compare Gaia BH2 and BH1 to other known BHs. Red and blue points show low- and high-mass X-ray binaries, whose parameters we take from Remillard & McClintock (2006) and the BlackCAT catalog of X-ray transients introduced by CorralSantana et al. (2016). We also show the binaries VFTS 243 (in the LMC; Shenar et al. 2022) and HD 130298 (in the Milky Way; Mahy et al. 2022), which both are single-lined binaries containing GLYPH<24> 25 𝑀 GLYPH<12> O stars and unseen companions suspected to be BHs. In magenta, we show two binaries with suspected BH companions in the globular cluster NGC 3201, discovered with MUSE (Giesers et al. 2018, 2019).  \nFor the BHs in the Milky Way - which includes all the LMXBs but only 2 of the HMXBs - we additionally collected distance estimates and quiescent 𝐺 -band magnitudes. 5 systems - A0620-00, Cyg X1, V4641 Sgr, GRO J1655-40, and HD 130298 - have reasonably precise Gaia DR3pallaxes, with pallax\\_over\\_error >5.Forthese sources, we use the geometric distance estimates from Bailer-Jones et al. (2021), which are informed by the parallax and a Galactic model prior. For the other sources, which are fainter and/or more distant, the parallaxes are not very constraining, and so we collect distance estimates from the literature that are based on the orbital period, donor spectral type, and apparent magnitude in quiescence. Where available, we take these estimates from Jonker & Nelemans (2004). For the objects not included in that work, we take the estimates from individual-object papers: Casares et al. (2009) for GS 1354-64, Heida et al. (2017) for GX 339-4, Mikołajewska et al. (2022) for MAXIJ1820+070, and Mata Sánchez et al. (2021) for MAXI J1305704. A few sources are too faint in the optical (in some cases due extinction) to be detected by Gaia ; these are not shown in Figure 14.  \n<!-- image -->  \n<!-- image -->  \nFigure 11. Left: MeerKAT and Chandra observations of Gaia BH2. Images are 2.7 arcmin wide. Red circles are centered on Gaia BH2 and have a radius of 3.5 arcsec, which is roughly the size of the MeerKAT beam and slightly larger than the Chandra PSF. The source is not detected in either image. Right: predicted X-ray and radio flux of Gaia BH2 under 2 scenarios: (1; red-dashed) the BH accretes the giant's wind at the Bondi-Hoyle-Lyttleton (BHL) rate, with the radiative efficiency calculated following models for radiatively ineffiecient accretion flows, and (2; solid black) the mass flux toward the BH decreases at smaller radii due to winds (Section 3.12.2 and Equation 18), such that the accretion rate at the event horizon is much lower than at the Bondi radius. Dashed lines show the observed 2 𝜎 upper limits. The data are inconsistent with accretion at the BHL rate but consistent with models in which most of the would-be-accreted mass is lost to winds.  \n<!-- image -->  \nGaia BH2 has the longest orbital period of all currently known and suspected BHs. Its orbital period is a factor of seven longer than that of Gaia BH1, and GLYPH<25> 40 times longer than that of the longest-period X-ray binary, GRS 1915+105. The masses of Gaia BH1 and BH2 are similar to the bulk of the BH X-ray binary population; only GRS 1009-45 and GRO J0422+32 host BHs of significantly lower mass.  \nBoth Gaia BH1 and BH2 are closer, and have longer orbital periods, than any other BH binaries. They are brighter than any of the BH X-ray binaries with low-mass donors, but fainter than the HMXB Cyg X-1 (which hosts a GLYPH<24> 400 GLYPH<150> 000 𝐿 GLYPH<12> luminous star). Given that their companion stars are cool and slowly rotating, their abundances and RVs can be measured with higher fidelity than in other BH companions. The two systems' close distances imply that widely separated BH binaries significantly outnumber close and accreting BH X-ray binaries.", '5 FORMATION HISTORY': 'As with Gaia BH1 (El-Badry et al. 2023), it is unclear whether Gaia BH2 formed from a primordial binary or via a more exotic channel. Wefirst consider formation from a primordial binary and then discuss alternative channels in Section 5.5.', '5.1 Natal kick constraints': "The formation of a compact object after core collapse is expected to impart a natal kick to the compact object, modifying the orbital size and eccentricity. Given the wide current orbit of Gaia BH2, it is unlikely that tides have appreciably circularized the binary following the formation of the BH. Assuming that the system formed from a pre-existing binary, its semi-major axis and eccentricity may have been changed by an asymmetric supernova kick and by instantaneous mass loss (e.g. Blaauw 1961). We used COSMIC (Breivik et al. 2020) and emcee (Foreman-Mackey et al. 2013) to model the combined binary properties at core collapse and natal kick strength based on the present day orbital period and eccentricity of Gaia BH2 using the 1 𝜎 measurement errors. We followed a similar procedure to the one described in Wong et al. (2022) and El-Badry et al. (2023). But instead of modeling the primordial binary properties, we initialized our model with the BH progenitor mass 12 𝑀 GLYPH<12> just before core collapse, secondary mass 1 GLYPH<149> 07 𝑀 GLYPH<12> , and zero eccentricity. This effectively marginalizes over all previous evolution up to the point of core collapse of Gaia BH2's progenitor. Our choice of BH progenitor mass is defined by the Fryer et al. (2012) 'delayed' model and produces a BH with mass 9 𝑀 GLYPH<12> . The Blaauw kick associated with the instantaneous mass loss imparts an eccentricity of GLYPH<24> 0 GLYPH<149> 3, so a  \nFigure 12. Comparison of Gaia BH1 and BH2 (black points) to other known BHs. Red and blue symbols correspond to accreting BHs with low- and high-mass companions. Magenta symbols show detached binaries in the globular cluster NGC 3201, and cyan points show detached binaries in which the luminous star is a high-mass ( & 20 𝑀 GLYPH<12> ) star. Gaia BH2 stands out from the rest of the population most strongly due to its orbital period, which is 7 GLYPH<2> longer than that of Gaia BH1 and GLYPH<24> 1000 times longer than the periods of typical X-ray binaries.  \n<!-- image -->  \nFigure 13. Comparison of Gaia BH1 and BH2 (black points) to known Galactic BHs in the plane of distance and orbital period. Color scheme is as in Figure 12. The Gaia -discovered systems are in entirely different part of this parameter space from other known BHs: they both have longer periods and are closer to Earth than any other known BHs.  \n<!-- image -->  \nFigure 14. Comparison of Gaia BH1 and BH2 (black points) to known Galactic BHs in the plane of distance and quiescent optical magnitude. Color scheme is as in Figure 12. Gaia BH1 and BH2 are nearer than any of the other objects. They are significantly brighter than the BH LMXBs (red), but not the HMXBs, whose donors are intrinsically much more luminous.  \n<!-- image -->  \nmoderate natal kick is required to produce the observed eccentricity, assuming that the orbit was previously circularized by tides.  \nThe strength of the natal kick is degenerate with the orbital separation at core collapse, 𝑎 cc, which we therefore include as a parameter in our model. We used uniform priors for the orbital separation at core-collapse between 1 and 5 au, the natal kick velocity between 0 and 100 km s GLYPH<0> 1 , and natal kick angles that are uniform on a sphere. We ran 1024 walkers for 50 GLYPH<150> 000 steps, thinned the chains by a factor of 10, and retained only the last 2000 steps to ensure a proper burn in. We find that the present-day orbital period and eccentricity prefer a core-collapse semimajor axis 𝑎 cc = 2 GLYPH<149> 9 GLYPH<6> 0 GLYPH<149> 93 au and natal kick speed 𝑣 kick = 36 , 21 GLYPH<0> 11 kms GLYPH<0> 1 . Assuming Gaia BH2 formed as a primordial binary, this suggests that the semi-major axis was likely smaller in the past than its present-day value. The exact constraints depend on the assumed pre-supernova mass of the progenitor, which is uncertain. However, a generic result of this modeling is that the observed wide orbit and moderate eccentricity disfavor strong natal kicks to the BH.", '5.2 BH progenitor mass': "Given the BH's mass of GLYPH<24> 9 𝑀 GLYPH<12> , its progenitor likely had an initial mass 𝑀 init & 25 𝑀 GLYPH<12> (e.g. Sukhbold et al. 2016; Raithel et al. 2018). Because more massive stars are predicted to have stronger winds, the relation between initial mass and BH mass for stars of near-solar metallicity is predicted to be non-monotonic, and - given current observational constraints - essentially any initial mass 𝑀 init & 25 𝑀 GLYPH<12> could plausibly have produced a 9 𝑀 GLYPH<12> BH.", '5.3 Problems with formation through common envelope evolution': "If Gaia BH2 formed as a primordial binary with a & 25 𝑀 GLYPH<12> primary and a GLYPH<24> 1 𝑀 GLYPH<12> secondary, mass transfer through Roche lobe overflow of the primary would almost certainly lead to an episode of common envelope evolution due to the extreme mass ratio. This is predicted to reduce the orbital separation by a factor of 100 GLYPH<0> 1000, depending on the structure of the primary when mass transfer began (see the discussion in El-Badry et al. 2023). Since Roche lobe overflow could not have begun at a separation significantly larger than 10 au, this leads to a predicted post-common envelope separation at least a factor of 50 closer than observed. Given the system's low space velocity and modest eccentricity, we also cannot invoke a fine-tuned natal kick to the BH to widen the orbit again after a common-envelope episode. The same considerations also disfavored a common envelope for Gaia BH1, and Gaia BH2's wider orbit exacerbates the tension.", '5.4 Formation from a star too massive to become a red supergiant': "We next consider whether Gaia BH2 could have formed from a massive star that never became a red supergiant and thus avoided mass transfer altogether. We modeled the evolution of the progenitor with single massive star evolutionary models of solar metallicity, calculated with MESA by Klencki et al. (2020). The orbital separation at the time of the progenitor's death is unknown, since a natal kick of unknown magnitude and direction could have modified the orbit in a variety of ways. But as discussed above, the modest eccentricity suggests that the final pre-supernova separation was likely not too different from the current semi-major axis ( GLYPH<24> 5 au). We adopt a value of 3 au.  \nThe massive star models undergo significant mass loss through  \nwinds. Assuming this occurs isotropically on a timescale that is long compared to the orbital period, the orbit will expand such that the product of the semi-major axis and total binary mass is conserved. This means that the orbit must have been tighter in the past in the absence of mass transfer. For each massive star model, we assume a final separation of 3 au, and then calculate the separation at each timestep based on the mass at that time. We also calculate the Roche lobe radius, using the fitting formula from Eggleton (1983) and assuming a 1 𝑀 GLYPH<12> companion. Finally, we calculate the ratio of the massive star model's radius to the Roche lobe. We plot the maximum value of this quantity across all timesteps in Figure 15. If this ratio is larger than 1, it means that the BH progenitor would have overflowed its Roche lobe at some point if it were in the Gaia BH2 system, leading to a common envelope inspiral.  \nAll models with 𝑀 init GLYPH<20> 65 𝑀 GLYPH<12> would overflow their Roche lobes in the Gaia BH2 system by more than a factor of 2, meaning that they could not have avoided a common envelope. However, the most massive models behave quite differently: these models lose their hydrogen envelopes to strong winds while still on the main sequence and never become red giants (see e.g. Smith & Conti 2008). If the progenitor of Gaia BH2 followed such an evolutionary channel, the system could have formed as, e.g., a 70 , 1 𝑀 GLYPH<12> primordial binary with separation GLYPH<24> 0 GLYPH<149> 5 au, completely avoided mass transfer through Roche lobe overflow, and died as a Wolf Rayet + 1 𝑀 GLYPH<12> binary with separation GLYPH<24> 3 au. The mass threshold above which stars avoid a red supergiant phase varies significantly between models, and is as low as GLYPH<24> 40 𝑀 GLYPH<12> in some models (e.g. Agrawal et al. 2020; Bavera et al. 2022).  \nAn important caveat to these calculation is that the maximum radius of massive stars depends sensitively on parameters that are uncertain, particularly their wind mass loss rates. Some recent studies (e.g. Beasor et al. 2020) suggest that winds of red supergiants with 𝑀 & 25 𝑀 GLYPH<12> are weaker than predicted by classical models. This may imply that the mass loss rates of more massive stars are also lower, in whichcaseenvelope loss while on the main sequence might not occur. The mass threshold above which models avoid a red supergiant phase also depends significantly on assumptions about rotation (Maeder 1987). We conclude that there are significant uncertainties in the evolutionary models for massive stars that avoid expansion, but it is possible that Gaia BH2 formed through such a channel.", '5.5 Other formation channels': "It is also possible that Gaia BH2 did not form from a primordial binary. One possibility is that the system formed dynamically through an exchange interaction in a dense cluster. The binary's thin disklike Galactic orbit precludes formation in a classic globular cluster, but the system may have formed in an open cluster that has since dissolved. This scenario does not leave any smoking gun signatures and is thus difficult to test definitively for any single binary. N-body simulations of binaries in clusters (e.g. Kremer et al. 2018) that are tailored to systems like Gaia BH1 and BH2 are required to determine whether such dynamical exchange+ejection events in clusters are common enough to explain the observed rate of wide BH + normal star binaries. If the observed 𝛼 GLYPH<0> enhancement of the giant is a result of pollution from the BH progenitor (see Section 3.9.1), this would disfavor a dynamical formation channel.  \nIt is also possible that Gaia BH2 formed as a triple, with the red giant an outer tertiary to an inner binary containing two massive stars. Interactions between the two stars could have prevented either one from becoming a red supergiant. This scenario is discussed further by El-Badry et al. (2023) for Gaia BH1; a similar scenario could  \nFigure 15. Predictions of single massive star evolutionary models for different initial masses of the BH progenitor. Top panel shows the maximum radius reached by the model, and gray shaded region shows the periastronto-apastron separation of Gaia BH2's current orbit. Bottom panel shows the maximum ratio of the star's radius to the Roche lobe radius, assuming a presupernova separation of 3 au and adiabatic orbital evolution due to mass loss. For all progenitor masses below 65 𝑀 GLYPH<12> , the model would overflow its Roche lobe in the Gaia BH2 system, presumably leading to an episode of common envelope evolution. Progenitors with initial masses 𝑀 init ¡ 65 𝑀 GLYPH<12> , however, lose their envelopes to strong winds at the end of their main-sequence evolution and never become red supergiants. Such models could have avoided a common envelope episode in the Gaia BH2 system.  \n<!-- image -->  \n/circledot  \nalso work for Gaia BH2. Indeed, a triple scenario likely requires less fine-tuning for Gaia BH2 than for BH1, because the wider orbit makes it possible to accommodate a wider range of inner binaries without the system becoming unstable. As with Gaia BH1, precision RVmonitoring of the luminous star may eventually detect or rule out a dark inner binary, which would cause a short-period, low-amplitude modulation of the overall RV curve.", '5.6 Future evolution': "Over the next GLYPH<25> 100 Myr, the luminous star is expected to continue ascending the red giant branch, until it ignites core helium burning. The MIST models shown in Figure 1 predict that if it evolved in isolation, the giant would reach a radius of about 170 𝑅 GLYPH<12> just before the helium flash. However, in the current binary orbit, the giant would overflow its Roche lobe at periastron when it reached 𝑅 GLYPH<25> 110 R GLYPH<12> . The orbit will likely circularize as the giant expands, and there will likely be a brief period of mass transfer through stable Roche lobe overflow near the tip of the giant branch, when the system will be observable as a symbiotic X-ray binary. After helium ignition, the giant will contract and the binary will appear as a wide BH binary with a red clump or hot subdwarf star. Another episode of stable Roche lobe overflow is possible after core helium burning, depending  \non how much of the star's envelope mass remains. Either way, the binary is expected to terminate its evolution as a wide BH + white dwarf binary, probably with a low-eccentricity orbit.", '6.1 Completeness of follow-up': 'Weidentified both Gaia BH1 and BH2 as BH candidates using simple astrometric cuts designed to find massive and dark companions to low-mass stars (see Appendix E of El-Badry et al. 2023). We have now followed-up all of the 6 candidates those cuts yielded. Two have been validated, and the remaining four turned out to have spurious Gaia orbital solutions.  \nAbroader search for compact object companions among the Gaia DR3 astrometric solutions was carried out by Shahaf et al. (2023). Gaia BH2 was excluded from their sample both because the luminous source is not on the main sequence and because its orbital period is longer than 1000 days. Besides Gaia BH1, Shahaf et al. (2023) identified an additional seven BH candidates, whose Gaia orbital solutions appeared to imply a dark companion with 𝑀 2 ¡ 2 𝑀 GLYPH<12> . We are carrying out RV follow-up of these sources, which we will describe in detail elsewhere. Only for one of them, Gaia DR3 6328149636482597888, do the RVs inspire confidence in the Gaia solution. However, that source is a low-metallicity halo star, and so the masses of both the luminous star and the unseen companion are lower than inferred when the system is interpreted with solar-metallicity models (see El-Badry et al. 2023). The best-fit companion mass is 𝑀 2 GLYPH<25> 2 GLYPH<149> 2 𝑀 GLYPH<12> . This is consistent with a variety of astrophysical scenarios; we defer more detailed analysis of this source to future work.  \nIn summary, we suspect that Gaia BH1 and BH2 are the only BHs with DR3 astrometric orbital solutions, at least among sources with implied companion masses 𝑀 2 & 2 GLYPH<149> 5 𝑀 GLYPH<12> . There are, however, a few dozen candidates with 1 . 𝑀 2 . 2 𝑀 GLYPH<12> - most readily interpreted as neutron stars and white dwarfs - still to be followed-up.', '6.2 Is there a BH - neutron star mass gap?': "It is conspicuous that there are now two vetted Gaia compact object binary solutions with 8 GLYPH<159> 𝑀 2 GLYPH<157> 𝑀 GLYPH<12> GLYPH<159> 10, a few dozen candidates with 1 . 𝑀 2 . 2 GLYPH<149> 5 𝑀 GLYPH<12> , and no reliable solutions with 3 GLYPH<159> 𝑀 2 GLYPH<157> 𝑀 GLYPH<12> GLYPH<159> 8. The implied BH/NS companion mass distribution is reminiscent of the mass gap observed for BHs in X-ray binaries (e.g. Kreidberg et al. 2012), though the mean BH mass inferred for X-ray binaries is somewhat lower ( GLYPH<24> 7 𝑀 GLYPH<12> ). While Gaia is somewhat more sensitive to high-mass companions, which produce a larger photocenter wobble at fixed orbital period, both Gaia BH1 and BH2 would plausibly have been detectable if the BH had a mass of, say, 5 𝑀 GLYPH<12> . The data are thus already in tension with a scenario in which lowmass BHs significantly outnumber GLYPH<24> 10 𝑀 GLYPH<12> BHs, at least among BHs with stellar companions in wide orbits. The fact that Shahaf et al. (2023) identified several dozen wide neutron star candidates with dynamically-implied masses of 1 GLYPH<0> 2 𝑀 GLYPH<12> also suggests that Gaia would be sensitive to lower-mass BHs, but further follow-up is required to vet the neutron star candidates.  \nThenext decade will likely see the discovery of both more wide BH +normal star astrometric binaries and isolated single BHs discovered via astrometric microlensing (e.g. Sahu et al. 2022; Lam et al. 2022; Mroz et al. 2022). Comparison of the mass distributions of BHs  \nFigure 16. Orbital period distribution of all astrometric binaries in Gaia DR3 (black) and Gaia BH1 and BH2 (red). In the bottom panel, we divide the observed period distribution of all astrometric binaries by the expected intrinsic binary period distribution for solar-type stars. This provides an approximate sensitivity function of Gaia DR3 as a function of orbital period, marginalized over other observables. The two BHs have periods at the longand short-period edges of the full observed period distribution, where DR3's sensitivity is GLYPH<24> 10 GLYPH<2> lower than at its peak at 1.5 years. This is suggestive of a bimodal intrinsic period distribution, with fewer BHs at periods of 1-3 years.  \n<!-- image -->  \nwithin these populations may shed light on the formation pathways of the BH binaries.", '6.3 Orbital period distribution': 'In Figure 16, we compare the orbital periods of Gaia BH2 and BH1 to the period distribution of all GLYPH<25> 186 GLYPH<150> 000 astrometric binaries published in DR3. Most of these systems are solar-type binaries, whose intrinsic period distribution is well-approximated by a lognormal that peaks at a period of 10 5 days with a dispersion of 2.3 dex (Raghavan et al. 2010). To estimate the relative sensitivity of Gaia DR3 to binaries at different periods, we divide the observed period distribution by this lognormal. The latter decreases toward longer periods (for linear period intervals), so this rescaling gives more weight to binaries with long periods.  \nFigure 16 shows that Gaia DR3 was most sensitive to binaries with orbital periods between 1.25 and 2.5 years, since shorter periods produce smaller astrometric wobbles at fixed mass, and longer periods could not be fully sampled during the GLYPH<24> 1000 day observing baseline for data published in DR3. The strongly reduced sensitivity to periods of one year, due to confusion of orbital and parallax motion, is also evident. The fact that the two BHs discovered from Gaia DR3 have orbital periods on the long- and short-period edges of the observational sensitivity curve hints that the BH companion period distribution may be bimodal, with a valley at 1-3 year periods. The strength of the conclusions we can draw is of course limited by the  \nfact that we have detected only two BHs. We can nevertheless draw some conclusions. If the intrinsic period distribution of BH binaries were flat in linear orbital period, 80% the detected binaries would have periods between 1 and 3 years. There is thus only a 0 GLYPH<149> 2 2 = 4% chance of detecting two BHs outside this range and none inside it if the underlying BH period distribution is flat.  \nAs discussed by El-Badry et al. (2023), Gaia DR4 and DR5 will likely enable the discovery of dozens of BH + normal star wide binaries, both because they will be based on a longer observing baseline than DR3 and because epoch-level astrometric data will be published for all sources, without stringent cuts on signal-to-noise ratio. DR4 will be based on GLYPH<24> 2000 days of observations, and DR5 will likely be based on GLYPH<24> 4000 days of observations. Gaia will thus eventually enable the discovery and characterization of BH binaries with periods up to GLYPH<24> 15 years (and perhaps longer if partial orbits can be fit uniquely; e.g., Andrews et al. 2021). The intrinsic period distribution of BH + normal star binaries is quite uncertain, but some models predict that it will rise steeply at 𝑃 orb & 10 years, where the BH binaries will have formed without the BH progenitor ever having overflowed its Roche lobe (e.g. Breivik et al. 2017; Chawla et al. 2022). This prediction depends sensitively on the strength of BH natal kicks, which can easily unbind the widest binaries. In any case, the empirical period and mass distributions of BH binaries will become much better sampled in future Gaia releases.', '7 CONCLUSIONS': "We have carried out follow-up observations of a newly identified black hole (BH) + luminous star binary discovered via Gaia astrometry and spectroscopy. The system, which we call Gaia BH2, consists of an apparently unremarkable lower red giant branch star in a 3.5year orbit with an unseen massive companion. Basic properties of the system are summarized in Table 2. Our main conclusions are as follows:  \n(i) The luminous source : the luminous source is a bright ( 𝐺 = 12 GLYPH<149> 3), nearby ( 𝑑 = 1 GLYPH<149> 16 kpc) star on the lower red giant branch, just below the red clump (Figure 1). From evolutionary models, we infer a mass 𝑀 ★ = 1 GLYPH<149> 07 GLYPH<6> 0 GLYPH<149> 19 𝑀 GLYPH<12> . From a high-resolution spectrum (Figure 7), we measured moderately subsolar metallicity ( » Fe GLYPH<157> H … = GLYPH<0> 0 GLYPH<149> 22) but strong enhancement of 𝛼 GLYPH<0> elements ( » 𝛼 GLYPH<157> Fe … = , 0 GLYPH<149> 26). The detailed abundance pattern is otherwise unremarkable (Table 4). Because the system's orbit is wide, the 𝛼 GLYPH<0> enhancement is most likely primordial, but it is also possible that it is a result of pollution from low-velocity ejecta during a failed supernova.  \n(ii) Orbit and companion mass : our RV follow-up over a 7-month period validates the Gaia orbital solution (Figure 3) and yields improved constraints on the companion mass (Figure 4). The orbit has a 1277-day period - longer by a factor of 7 than any other known BH binary - and is moderately eccentric, with 𝑒 = 0 GLYPH<149> 52. The inclination is GLYPH<25> 35 deg; i.e., relatively face-on. The dynamically implied mass of the companion is 𝑀 2 = 8 GLYPH<149> 9 GLYPH<6> 0 GLYPH<149> 3 𝑀 GLYPH<12> . This constraint depends only weakly on the assumed mass of the luminous star, but significantly on the inclination (Figure 5).  \n(iii) Limits on luminous companions : the SED from the UV to the near-IR (Figure 2) and high-resolution spectrum (Figures 7 and 8) are both well fitted by single-star models. A luminous companion contributing more than GLYPH<25> 3% of the light at 4000 Å is ruled out spectroscopically. The giant is cool ( 𝑇 eff = 4600 K ' , and so the fact that the source is faint in the UV rules out any scenario in which the companion is one or more main-sequence stars (Figure 2). These limits rule out any plausible luminous stellar companion if the Gaia  \nsolution is correct. Even if the inclination inferred from the Gaia orbit were spurious - in which case the RVs would allow a companion mass as low as 2 GLYPH<149> 6 𝑀 GLYPH<12> - a stellar companion would be detectable.  \n(iv) Nature of the companion : the dynamically implied mass of the companion - 8 GLYPH<149> 9 GLYPH<6> 0 GLYPH<149> 3 𝑀 GLYPH<12> - is too large to be a neutron star or white dwarf, or any plausible combination of white dwarfs, neutron stars, and faint stars. The simplest explanation is a single BH. As with Gaia BH1, it is possible that the companion is an unresolved binary containing at least one BH.  \n(v) Formation history : the current orbit of Gaia BH2, with a periastron distance of 2 GLYPH<149> 4 au, is too close to accommodate BH progenitors with initial masses 20 . 𝑀 init GLYPH<157> 𝑀 GLYPH<12> . 65, which are predicted to become red supergiants with radii larger than the current orbit (Figure 15). One possibility is that the BH progenitor had a sufficiently high mass ( 𝑀 init & 65 𝑀 GLYPH<12> , though the precise value varies between models) to lose its envelope while on the main sequence, avoiding a red supergiant phase. It is also possible that the system was formed dynamically through an exchange interaction in a dense cluster, or through triple evolution.  \n(vi) X-ray and radio observations : we observed Gaia BH2 with both X-ray and radio facilities but did not detect it (Figures 11). Expressed in terms of the Bondi-Hoyle-Lyttleton (BHL) accretion rate, the non-detection implies a radiative efficiency . 10 GLYPH<0> 4 . The non-detection is consistent with model predictions for radiativelyinefficient accretion flows in which most of the inflowing material is lost to winds or accretion is suppressed by convection, such that only a small fraction of the BHL accretion rate actually reaches the event horizon. For typical radiative efficiencies, the non-detection is inconsistent with models in which a mass flux equal to the BHL rate reaches the event horizon.  \n(vii) The Gaia BHs in context : Gaia DR3 has thus far resulted in the identification of two widely separated BH + normal star binaries, whose properties are quite different from previously known BH binaries (Figures 12-14). Both have longer periods and are closer to Earth than any previously known BHs. This suggests that widely-separated BH binaries significantly outnumber close systems with ongoing mass transfer. Both have masses between 8 and 10 𝑀 GLYPH<12> . There do not appear to be any BHs with astrometric solutions published in DR3 that have masses between 3 and 8 𝑀 GLYPH<12> . The two BHs have periods at the short- and long-period edges of the Gaia sensitivity curve (Figure 16), perhaps hinting at a bimodal intrinsic period distribution for BH binaries.", 'ACKNOWLEDGEMENTS': "We thank Frédéric Arenou, Jim Fuller, and Andy Gould for useful discussions, and Maren Hempel, Paul Eigenthaler, and Régis Lachaume for carrying out the FEROS observations. We are grateful to the ESO, MeerKAT, Chandra , and Swift directorial offices and support staff for prompt assistance with DDT observations.  \nHWR acknowledges the European Research Council for the ERC Advanced Grant [101054731]. This research used pystrometry, an open source Python package for astrometry timeseries analysis (Sahlmann 2019). This work used Astropy, 5 a community-developed core Python package and an ecosystem of tools and resources for astronomy (Astropy Collaboration et al. 2022), and the CIAO software (Fruscione et al. 2006) provided by the Chandra X-ray Center (CXC). We acknowledge the use of public data from the Swift data archive.  \nThis project was developed in part at the Gaia Fête, held at the Flatiron Institute's Center for Computational Astrophysics in June 2022, and in part at the Gaia Hike, held at the University of British Columbia in June 2022. 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.  \nThe MeerKAT telescope is operated by the South African Radio Astronomy Observatory, which is a facility of the National Research Foundation, an agency of the Department of Science and Innovation. The scientific results reported in this article are based in part on observations made by the Chandra X-ray Observatory, which is operated by the Smithsonian Astrophysical Observatory for and on behalf of the National Aeronautics Space Administration under contract NAS8-03060.", 'DATA AVAILABILITY': 'Data used in this study are available upon request from the corresponding author.', 'A1 Predicted astrometry': 'Figure A1 shows the plane-of-the-sky motion of the photocenter of Gaia BH2 as it would be observed by Gaia . In the left panel, the dashed blue line shows the motion expected due to parallax and proper motion alone, while the black and gray lines show the combined effects of orbital motion, parallax, and proper motion. Lines are colored darker in the time interval covered by Gaia observations included in DR3, between July 2014 and May 2017. The interplay between parallax and proper motion leads to rather complex motion on the plane of the sky, which - given the limited time baseline of DR3-doesnot yet appear periodic. Nonetheless, the Gaia data alone constrain the orbit fairly precisely.  \nIn the right panel, individual cyan lines show draws from the posterior of the fit to the pure Gaia solution, with proper motion and parallax motion removed. Although there is noticeable dispersion in these predictions owing to observational uncertainties, the size and shape of the orbital ellipse are well-determined. The middle panel shows draws from the posterior of the joint fit of the Gaia data and RVs. As expected, the resulting orbital ellipse is consistent with the Gaia -only fit, but its parameters are more tightly constrained.', 'A2 Observation times': "The AstroSpectroSB1 solution for Gaia BH2 is based on both astrometric and RV measurements. Single-epoch astrometric and RV data is not published in DR3, but some insight into when a source was observed can be provided by the Gaia observation scheduling tool (GOST). 6 Given a source position, GOST provides a list of observation times when the scanning law predicts that a source will transit across the Gaia focal plane. For Gaia BH2, GOST predicts 97 CCDs transits during the observation window that contributed to DR3, distributed over 27 visibility periods (i.e., groups of observations separated by at least 4 days). Meanwhile, the gaia\\_source table reports matched\\_transits = 93 and visibility\\_periods\\_used = 26 . This suggests that a majority of the planned transits indeed contributed to the DR3 solution, but a few were omitted, as is common. Fewer transits contributed RVs, with rv\\_nb\\_transits = 32 and rv\\_visibility\\_periods\\_used = 12. This is also expected, since only a fraction of focal plane transits cover the RVS spectrometer (Gaia Collaboration et al. 2016).  \nIn Figure A2, we show the predicted astrometric positions and RVs at observation times predicted by GOST. It is important to note that (a) epoch-level data is not published in DR3, so there is no guarantee that the actual observed positions and RVs match the model predictions, (b) we have not simulated noise, (c) actual Gaia astrometric measurements are 1D, with the scan angle varying from visit to visit, and (d) while a majority (93/97) of the planned transits actually contributed data to the astrometric solution, only about one third (32/93) contributed an RV.  \nWecannevertheless draw several conclusions from Figure A2. The predicted astrometric measurement times span 1000 days - 79 % of the binary's orbital period - and cover most of the predicted plane-ofthe sky motion. The typical actually-realized along-scan uncertainty  \nFigure A1. Left: plane-of-the-sky photocenter motion as inferred from our joint RV + Gaia fit. Blue line shows the motion due to parallax and proper motion alone. Black and grey lines show the full motion, including wobble due to the star's orbit. Black line shows the time period covered by Gaia DR3 astrometry; gray line extrapolates beyond that period. Middle and right panels show constraints on the orbital motion alone, with parallax and proper motion subtracted. Right panel shows the constraints from the Gaia data alone, while middle panel shows constraints from the joint fit.  \n<!-- image -->  \nfor each astrometric measurement for a source with 𝐺 = 12 GLYPH<149> 3 is 0.2 mas (Lindegren et al. 2021a), which is comparable to the size of the red points in Figure A2. The expected per-epoch RV uncertainty for a source of Gaia BH2's apparent magnitude and spectral type is also about 1.5 km s GLYPH<0> 1 (see e.g. El-Badry et al. 2019, their Figure F1), also comparable to the size of the points in the lower right panel. The reported rv\\_amplitude\\_robust =37kms GLYPH<0> 1 (representing the total range in measured epoch RVs, after clipping of outliers) suggests that RVs were measured over most of the orbit's dynamic range in RV, though the two transits predicted near the periastron passage in 2016 likely did not result in an RV measurement.  \nWe conclude that the Gaia measurements likely are adequate to constrain the orbit robustly, though it would obviously be preferable to have access to the actual epoch-level data.", 'APPENDIX B: GAIA GOODNESS OF FIT': 'One of the quality diagnostics calculated for Gaia orbital solutions is the goodness\\_offit estimator, 𝐹 2 , which is defined (e.g. Halbwachs et al. 2022) as  \n𝐹 2 = √︂ 9 𝜈 2 " GLYPH<18> 𝜒 2 𝜈 GLYPH<19> 1 GLYPH<157> 3 , 2 9 𝜈 GLYPH<0> 1 # GLYPH<149> (B1)  \nHere 𝜈 is the number of degrees of freedom and 𝜒 2 is the sum of the uncertainty-normalized squared residuals between model and data. When there is good agreement between model and data and the uncertainties are reliably estimated, one expects 𝐹 2 to be normally distributed; i.e., 𝐹 2 GLYPH<24> N\' 0 GLYPH<150> 1 \' . Values of 𝐹 2 GLYPH<29> 1 thus indicate a formally poor fit. Gaia BH2 has 𝐹 2 = 3 GLYPH<149> 07, raising potential concern about the quality of the Gaia solution.  \nIn Figure B1, we compare the 𝐹 2 values of Gaia BH2 and BH1 to those of all astrometric binary solutions published in DR3 (including both Orbital and AstroSpectroSB1 solutions). This comparison reveals that Gaia BH2 actually has a somewhat lower 𝐹 2 value (indicative of a better fit) than typical sources at its apparent magnitude: almost all sources with 𝐺 GLYPH<159> 13 have 𝐹 2 ¡ 1. This most likely is a  \nFigure A2. Predicted observation times of Gaia BH2 from the Gaia observation scheduling tool (GOST). Black line shows the best-fit orbit from our combined fit. Red points show the predicted photocenter positions at the times when GOST predicts Gaia would have observed the source. Note that we do not have access to the actual measured Δ RA, Δ Dec and RV values; only to the predicted scan times. Further caveats are discussed in the text.  \n<!-- image -->  \nresult of the change in Gaia window class at 𝐺 = 13 (Rowell et al. 2021): epoch astrometry for sources with 𝐺 GLYPH<159> 13 is derived from a 2D window fitted with a point spread function, while epoch astrometry for fainter sources is calculated from 1D data. This leads to discontinuities in many astrometric quantities at 𝐺 = 13 (e.g. El-Badry et al. 2021; Lindegren et al. 2021b).  \nThe fact that Gaia BH2 has 𝐹 2 ¡ 1 does not necessarily imply that the uncertainties on the AstroSpectroSB1 solution - which also enter into our joint fit - are underestimated. In calculating the  \nFigure B1. Reported goodness\\_offit ( 𝐹 2; Equation B1) for all Gaia astrometric binaries, with Gaia BH1 and BH2 highlighted. Large 𝐹 2 values indicate a poor formal fit. There is a discontinuity in the typical 𝐹 2 values at 𝐺 = 13 owing to a quirk of the Gaia data processing: brighter sources typically have larger 𝐹 2 values because the uncertainties in their individualepoch astrometric data are underestimated. Both Gaia BH1 and BH2 have typical 𝐹 2 values for their apparent magnitude, indicative of an unremarkable astrometric solution.  \n<!-- image -->  \nastrometric solution, the uncertainties of the epoch-level data are inflated by multiplying by a constant, 𝑐 = √︂ 𝜒 2 𝜈 GLYPH<0> 1 GLYPH<16> 1 GLYPH<0> 2 9 𝜈 GLYPH<17> GLYPH<0> 3 , after which 𝐹 2 = 0 (Halbwachs et al. 2022; Pourbaix et al. 2022).', 'APPENDIX C: TESS LIGHT CURVES': 'The eleanor light curves of Gaia BH2 (TIC 207885877) are shown in Figure C1. We extract these using a 15 GLYPH<2> 15 pixel target pixel file, with the background estimated in a 31 GLYPH<2> 31 pixel region. We use the systematics-corrected flux in the default aperture chosen by eleanor ; i.e., the aperture with the lowest combined differential photometric precision (CDPP) on one-hour timescales. The expected flux contamination from background stars in the aperture is of order 10%.  \nThe source is clearly not strongly variable on short timescales. Much of the apparent variability is noise; some variability on longer timescales may be the result of imperfect background subtraction and removal of systematics. The light curve from sector 38 is noisier than the one from sector 11, mainly because sector 38 used 3 GLYPH<2> shorter exposure times.  \nPower spectra for both sectors are shown in the right panels of Figure C1. The red line shows the result of smoothing the power spectra with a Gaussian kernel with 𝜎 = 5 𝜇 Hz. Several potentially significant peaks are evident in the smoothed power spectra. We tentatively identify the peak with 𝜈 = 61 𝜇 Hzas the global asteroseismic parameter 𝜈 max, since it is the strongest feature evident in data from both sectors. We emphasize that the identification is tentative since the time baseline and signal-to-noise ratio of the TESS data is still insufficient to unambiguously constrain 𝜈 max in a source of this apparent magnitude (e.g. Stello et al. 2022). The implied surface gravity is log 𝑔 GLYPH<25> log 𝑔 GLYPH<12> , log GLYPH<0> 𝜈 max GLYPH<157> 𝜈 max GLYPH<150> GLYPH<12> GLYPH<1> , 1 2 log GLYPH<0> 𝑇 eff GLYPH<157> 𝑇 eff GLYPH<150> GLYPH<12> GLYPH<1> GLYPH<25> 2 GLYPH<149> 68, where 𝑔 is in cgs units and log 𝑔 GLYPH<12> GLYPH<25> 4 GLYPH<149> 43, 𝜈 max GLYPH<150> GLYPH<12> = 3090 𝜇 Hz, and 𝑇 eff GLYPH<150> GLYPH<12> = 5777 K.  \n6  \nFigure C1. TESS data for Gaia BH2. Left panels show normalized light curves from the full frame images, with an exposure time of 1800s for sector 11 and 600s for sector 38. Right panels show power spectra, and red lines show power spectra smoothed with a Gaussian kernel with 𝜎 = 5 𝜇 Hz. Vertical dashed line shows a tentative detection of 𝜈 max.  \n<!-- image -->'}