Open Access
Issue
A&A
Volume 695, March 2025
Article Number L1
Number of page(s) 9
Section Letters to the Editor
DOI https://doi.org/10.1051/0004-6361/202453535
Published online 25 February 2025

© The Authors 2025

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

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

1. Introduction

The Gaia BH3 system recently discovered from pre-release Gaia DR4 astrometry (Gaia Collaboration 2024) is composed of a low-mass giant star (0.76 ± 0.05 M) and a high-mass dark component (32.70 ± 0.82 M) that is most plausibly a black hole (BH). Gaia repeatedly measured the position of the photocenter of Gaia BH3 (LS II +14 13, Gaia DR3 4318465066420528000, 2MASS J19391872+1455542), that is, the position of the star assuming that the BH flux is negligible. This displacement effectively traces the orbital motion of the star around the center of mass of the system (see, e.g., Brandt 2024). Among the 1.5 million examined orbital solutions in the Gaia DR4 preliminary NSS data, Gaia BH3 has the largest astrometric mass function fM = 32.06 ± 0.64, indicating that this is a very exceptional object within several kiloparsecs from the Sun. The high mass of the dark component makes it the most massive of all known BHs in the Milky Way, apart from the supermassive Sagittarius A*. The age of the very low metallicity ([Fe/H] = −2.56 ± 0.11) giant star component is estimated to be ≈12 Ga1 (Gaia Collaboration 2024). At a distance of only 590 pc, it is the second nearest BH after Gaia BH1 (480 pc; El-Badry et al. 2023; Chakrabarti et al. 2023), which has a mass of ≈10 M. The discovery of Gaia BH3 offers an opportunity to search for emission from the close environment of this BH. Gaia BH3 is photometrically quiet, but the BH may be accreting a fraction of the stellar wind emitted by the giant companion as observed in symbiotic binaries (Lutovinov et al. 2020), potentially resulting in measurable emission. Remarkably, Gaia BH3 is similar to the BH merger progenitors detected through gravitational waves (GWs) that exhibit masses of 20–40 M. In particular, it is an analog of the individual 36 ± 4 and 31 ± 4 M BHs that produced the first GW-detected merger, GW150914 (Abbott et al. 2016). The high mass of Gaia BH3 also fits within the expected range for the early-Universe seeds of supermassive BHs (Aggarwal 2024).

We report on the first interferometric observations of Gaia BH3 collected with VLTI/GRAVITY. The observational setup and observing log are described in Sect. 2, and the search for an observational signature from the BH component is discussed in Sect. 3. Our conclusions are presented in Sect. 4.

2. Observations

We observed Gaia BH3 on 3 June, 23 June, and 18 July 2024 using the GRAVITY beam combiner (GRAVITY Collaboration 2017) of ESO’s Very Large Telescope Interferometer (VLTI; Mérand et al. 2014; Woillez et al. 2018; Haubois et al. 2020, 2022) at Cerro Paranal observatory (Chile). Although the apparent magnitude of the stellar component in Gaia BH3 is relatively bright in the infrared (mK = 9.08), we employed the four 8-meter Unit Telescopes (UTs) to reach a high S/N ratio and a correspondingly high sensitivity in terms of companion contrast. The UTs were equipped with the first generation MACAO adaptive optics system (Arsenault et al. 2003). The calibrator for the interferometric transfer function BD+08 4132 (Gaia DR3 4296148519431266176, 2MASS J19322054+0848235) was selected using the JMMC2SearchCal tool (Bonneau et al. 2011; Chelli et al. 2016) considering its proximity to Gaia BH3 (6.4°) and its similar magnitude in the K band. BD+08 4132 was observed immediately after Gaia BH3 using an identical spectral resolution and polarization setting of GRAVITY. The photospheres of both the giant star component in Gaia BH3 and the calibrator are unresolved by the interferometer (fringe visibility V > 99.8%). A summary of the properties of Gaia BH3 and BD+08 4132 is given in Table 1, and the log of the collected observations of Gaia BH3 is presented in Table 2. A description of the GRAVITY data sets at each observing epoch is provided in App. A. In the following analysis, we consider only the data sets collected on 23 June and 18 July 2024.

Table 1.

Summary of the properties of the giant star in LS II +14 13 (Gaia BH3) and the interferometric calibrator BD+08 4132.

Table 2.

Log of the observations of Gaia BH3 (LS II +14 13) and the interferometric calibrator (BD+08 4132) with GRAVITY.

3. Analysis

3.1. Predicted black hole properties

The orbital elements determined from Gaia astrometry and spectroscopic radial velocities by Gaia Collaboration (2024) provide a precise prediction of the positions of the BH and giant star components of Gaia BH3 (Fig. 1). On 23 June 2024, the positional offset from the giant star to the BH was (Δα, Δδ) = (+16.7, +6.6) mas, and it was ( + 16.4, +6.9) mas for the 18 July 2024 epoch. Knowing the position facilitates the search for infrared emission at the position of the BH in the GRAVITY data. The radial velocity of both the star and the BH is close to vR = −350 km s−1 (the systemic radial velocity of Gaia BH3). Gaia BH3 appears as a quiet BH in the UVES spectrum of the star, with no detectable emission lines (Gaia Collaboration 2024). While the radius of the event horizon of the BH itself is very small (≈96 km), its Roche lobe is very extended, thanks to its dominant mass. The Lagrangian L1 point of the system (i.e., the point between the two bodies where the gravitational forces from the star and the BH balance each other in the rotating frame) varies in position on the star-BH axis over the orbital cycle due to the eccentricity and varying orbital angular velocity. As shown in Fig. 1, the position of L1 in the Gaia BH3 system is approximately at 20% of the star-BH separation from the star toward the BH. Due to this peculiar geometrical configuration, a significant fraction of the wind emitted by the giant star is intercepted by the Roche lobe of the BH, and part of it may be accreted. The wind of the giant star may therefore ensure a steady flow toward the BH and a stable emission from accretion.

thumbnail Fig. 1.

Position of Gaia BH3 and its companion star relative to their center of mass at the epoch of GRAVITY observations. The position of the Lagrangian L1 point of the system is shown with a purple diamond symbol.

3.2. Search for continuum emission

We used the companion analysis tool CANDID3 (Gallenne et al. 2015) to systematically search for the signature of Gaia BH3 in the GRAVITY data. CANDID uses a two-dimensional position grid and a multi-parameter fit based on a least-squares minimization algorithm. We jointly considered the data from the 23 June and 18 July 2024 observing epochs in the search for BH emission, as the shift in position of the giant star is only 400 μas. We checked through a numerical simulation that this shift has a negligible impact on the sensitivity of the companion search. The presence of a secondary emission in addition to the giant star in the field of view of GRAVITY has an influence on several interferometric observables and on the closure phase (CP; or T3PHI) in particular. This observable has a low sensitivity to biases introduced by the atmospheric piston effect (time-dependent variation of the optical delay between the beams collected by each telescope). The CP is a very good tracer of any deviation from central symmetry of the light distribution in the field of view, which may be caused by a secondary source. To ensure a sufficiently complete coverage of the χ2 map, we conducted the search on a grid of starting positions from where the χ2 of the fit is minimized. The detection limit of CANDID is determined based on the injection and recovery of artificial companions over the search grid. At each search position (Δα, Δδ), CANDID outputs the relative flux with respect to the primary star of a secondary source compatible with the input data. Examples of recent applications of this algorithm to the detection of faint companions are presented in Evans et al. (2024) and Gallenne et al. (2025). The CANDID maps of detection probability and sensitivity are shown in Fig. 2. At the predicted location of the BH, we determined a contrast limit of Δm = 6.8 mag, that is, C = fBH/f = 0.19% over the K band (λ = 2.0 − 2.4 μm). The expected flux of the giant star from its 2MASS Ks magnitude (Skrutskie et al. 2006) is

f K , = [ 1.01 ± 0.01 ] × 10 13 W m 2 μ m 1 . $$ \begin{aligned} f_{K, \star } = [1.01 \pm 0.01]\times 10^{-13}\,\mathrm{W\,m^{-2}\,\upmu m}^{-1}. \end{aligned} $$(1)

thumbnail Fig. 2.

Top panel: CANDID map of companion detection probability expressed in the number of standard deviations of the noise level for the combination of 23 June and 18 July 2024 data in medium spectral resolution. The S/N color scale is from 0 to 2.1σ. The position of the star is marked with a star symbol, and the maximum S/N ratio of 2.3 is indicated with a red cross. The predicted position of the BH is framed with a white square. Bottom panel: Contrast map for a 3σ detection level around the expected position of the BH, where the color scale represents the contrast Δm in magnitudes.

The contrast limit translates into a BH flux upper limit of

f K , BH < 1.9 × 10 16 W m 2 μ m 1 . $$ \begin{aligned} f_{K, \mathrm{BH}} < 1.9 \times 10^{-16}\,{W\,\mathrm{m}^{-2}\,\upmu \mathrm{m}}^{-1}. \end{aligned} $$(2)

To validate this CANDID detection limit, we also applied the PMOIRED4 code (Mérand 2022) to the combined 23 June and 18 July 2024 data set. This library fits user-defined multi-component source models to model astronomical spectro-interferometric data. The minimum χ2 is determined using a classical gradient descent method starting from a grid, resulting in best-fit values, uncertainties, and a correlation matrix. The contrast map around the position of the BH is presented in the left panel of Fig. 3. The simulated CPs (T3PHI) and normalized visibilities (N|V|) corresponding to an artificial companion with a flux of 5% of that of the star at the predicted location of the BH are presented in Fig. A.5. These data are compared to the 23 June and 18 July 2024 combined data set to illustrate the absence of a significant detection of BH emission. The PMOIRED limiting contrast for the continuum BH emission is consistent with the Δm = 6.8 mag value obtained with CANDID (Fig. 3, right panel). Finally, we obtained a similar contrast limit compared to CANDID and PMOIRED by applying the exoGravity5 code (GRAVITY Collaboration 2020; Lacour et al. 2020).

thumbnail Fig. 3.

Left panel: PMOIRED map of the limiting magnitude (3σ, 99% confidence) of the BH around its predicted position for the 23 June and 18 July 2024 combined data set. The displayed field corresponds to the white square in Fig. 2. Right panel: statistical distribution of the magnitudes of the recovered 3σ companions. Here, the dashed line represents the CANDID detection limit.

3.3. Hydrogen Brackett γ line

The high spectral resolution setting of the GRAVITY science combiner offers the possibility to search for BH emission in the hydrogen Brγ line. At this wavelength, the contrast may be reduced due to the combination of absorption in the star’s spectrum and emission around the BH. It is, however, less efficient in terms of contrast limit for the detection of continuum emission. The 23 June 2024 high spectral resolution data do not exhibit any significant signature in CP and triple amplitude in the Brγ line blue shifted at 2.164 μm (vR = −350 km s−1; Fig. A.6). We determined the detection limit in the Brγ line from this data set by assuming the BH has a possible emission line with a velocity offset of ±30 km s−1 around the expected velocity and a velocity dispersion δv between 50 and 500 km s−1. The maximum emission line strength for each given model was then estimated by computing the amplitude equivalent to a 3σ detection. The 18 July 2024 medium resolution data also do not show any significant Brγ signal (Fig. A.7). As shown in Fig. A.8, the detection limit mostly depends on δv, and it reaches Δm = 4.5 mag for δv = 300 km s−1. This behavior can be intuitively explained as follows: The higher the velocity dispersion, the more spectral channels are involved in the potential detection.

4. Conclusions

We did not detect emission from the BH in Gaia BH3 in the K band at a level of fK, BH3 < 1.9 × 10−16 W m−2 μm−1. This flux limit is shown in Fig. 4 together with the upper limits established by Cappelluti et al. (2024) in the X-ray domain from Chandra ACIS-S observations and by Sjouwerman & Blanchard (2024a,b) in the radio domain (ν = 3, 10 and 22 GHz). As discussed by Cappelluti et al. (2024), accretion of the wind from the giant star (whose mass loss rate is estimated to be between 10−11 and 10−9M a−1) by the BH is likely to occur through an advection-dominated accretion flow (ADAF; Mahadevan 1997). In this configuration, the spectral energy distribution of the accreting BH including synchrotron, inverse Compton, and bremsstrahlung (Pesce et al. 2021) is determined by the BH mass and the Eddington ratio fEdd. This adimensional parameter is defined as the actual accretion rate divided by the Eddington rate (the accretion rate for which the BH radiates at the Eddington luminosity). As shown in Fig. 4, the GRAVITY limit is consistent with the upper limit of fEdd < 4.9 × 10−7 derived by Cappelluti et al. (2024).

thumbnail Fig. 4.

Upper limits on the spectral energy distribution of Gaia BH3 established from radio (August–September 2024; Sjouwerman & Blanchard 2024a,b) and X-ray (May 2024; Cappelluti et al. 2024) observations and the GRAVITY near-infrared observations (June–July 2024). The blue curve represents the accreting BH model presented by Cappelluti et al. (2024) for fEdd < 4.9 × 10−7.

Together with the non-detection in X-ray6 and radio domains of Gaia BH1 and BH2 recently reported by Rodriguez et al. (2024), the nearest stellar mass BHs appear to be particularly quiet. While the Gaia BH3 components are presently close to their maximum separation, the eccentricity of the orbit (Fig. 1) will bring them much closer together around 2030. This may result in a significant increase of the accretion efficiency of the giant star wind by the BH and an up to 50-fold boost in emission (Cappelluti et al. 2024). The new GRAVITY+ Adaptive Optics (GPAO) system (Bourdarot et al. 2024; GRAVITY Collaboration 2022; Woillez et al. 2024) will provide a threefold gain in photometric efficiency compared to the MACAO system for a bright target such as Gaia BH3. The reported contrast limit on the BH emission relies on only three hours of effective GRAVITY observing time, and deeper observations will realistically provide a gain in sensitivity by another factor of two. This combination of stronger emission at periastron and enhanced GRAVITY+ sensitivity by a factor of six opens up exciting perspectives to detect accretion-induced emission from Gaia BH3 toward the end of the 2020s.

Data availability

The GRAVITY calibrated data files are available at the CDS via anonymous ftp to cdsarc.cds.unistra.fr (130.79.128.5) or via https://cdsarc.cds.unistra.fr/viz-bin/cat/J/A+A/695/L1


1

We abbreviate the word “year” as “a” from the Latin “annum” following the recommendations of IAU (Wilkins 1990).

6

Available at http://www.astropy.org/

Acknowledgments

We thank Drs Nico Cappellutti and Fabio Pacucci for sharing the BH emission model presented in Fig. 4. Based on observations collected at the European Organisation for Astronomical Research in the Southern Hemisphere under ESO DDT programme 113.27R4. We are grateful to ESO’s Director General, Prof. Xavier Barcons, for the allocation of VLTI time to our program. This work has made use of data from the European Space Agency (ESA) mission Gaia (http://www.cosmos.esa.int/gaia), processed by the Gaia Data Processing and Analysis Consortium (DPAC, http://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. The research leading to these results has received funding from the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation program (project UniverScale, grant agreement 951549). MC acknowledges financial support from the Centre National d’Etudes Spatiales (CNES). This research has made use of Astropy, a community-developed core Python package for Astronomy (Astropy Collaboration 2013, 2018), the Numpy library (Harris et al. 2020), the Scipy library (Virtanen et al. 2020), the Astroquery library (Ginsburg et al. 2019) and the Matplotlib graphics environment (Hunter 2007). We used the SIMBAD and VizieR databases and catalog access tool at the CDS, Strasbourg (France), and NASA’s Astrophysics Data System Bibliographic Services. The original description of the VizieR service was published in Ochsenbein et al. (2000). This publication makes use of data products from the Two Micron All Sky Survey, which is a joint project of the University of Massachusetts and the Infrared Processing and Analysis Center/California Institute of Technology, funded by the National Aeronautics and Space Administration and the National Science Foundation.

References

  1. Abbott, B. P., Abbott, R., Abbott, T. D., et al. 2016, Phys. Rev. Lett., 116, 241103 [Google Scholar]
  2. Aggarwal, Y. 2024, MNRAS, 536, 3177 [Google Scholar]
  3. Arsenault, R., Alonso, J., Bonnet, H., et al. 2003, in Adaptive Optical System Technologies II, eds. P. L. Wizinowich, & D. Bonaccini, SPIE Conf. Ser., 4839, 174 [NASA ADS] [CrossRef] [Google Scholar]
  4. Astropy Collaboration (Robitaille, T. P., et al.) 2013, A&A, 558, A33 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  5. Astropy Collaboration (Price-Whelan, A. M., et al.) 2018, AJ, 156, 123 [Google Scholar]
  6. Bonneau, D., Delfosse, X., Mourard, D., et al. 2011, A&A, 535, A53 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  7. Bourdarot, G., Eisenhauer, F., Yazıcı, S., et al. 2024, ArXiv e-prints [arXiv:2409.08438] [Google Scholar]
  8. Brandt, T. D. 2024, PASP, 136, 073001 [NASA ADS] [CrossRef] [Google Scholar]
  9. Cappelluti, N., Pacucci, F., & Hasinger, G. 2024, ApJ, 973, 75 [NASA ADS] [CrossRef] [Google Scholar]
  10. Chakrabarti, S., Simon, J. D., Craig, P. A., et al. 2023, AJ, 166, 6 [NASA ADS] [CrossRef] [Google Scholar]
  11. Chelli, A., Duvert, G., Bourgès, L., et al. 2016, A&A, 589, A112 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  12. Claret, A., & Bloemen, S. 2011, A&A, 529, A75 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  13. El-Badry, K., Rix, H.-W., Quataert, E., et al. 2023, MNRAS, 518, 1057 [Google Scholar]
  14. Evans, N. R., Gallenne, A., Kervella, P., et al. 2024, ApJ, 972, 145 [NASA ADS] [CrossRef] [Google Scholar]
  15. Gaia Collaboration (Brown, A. G. A., et al.) 2021, A&A, 649, A1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  16. Gaia Collaboration (Panuzzo, P., et al.) 2024, A&A, 686, L2 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  17. Gallenne, A., Mérand, A., Kervella, P., et al. 2015, A&A, 579, A68 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  18. Gallenne, A., Evans, N. R., Kervella, P., et al. 2025, A&A, 693, A111 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  19. Ginsburg, A., Sipőcz, B. M., Brasseur, C. E., et al. 2019, AJ, 157, 98 [Google Scholar]
  20. GRAVITY Collaboration (Abuter, R., et al.) 2017, A&A, 602, A94 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  21. GRAVITY Collaboration (Nowak, M., et al.) 2020, A&A, 633, A110 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  22. GRAVITY Collaboration (Abuter, R., et al.) 2022, The Messenger, 189, 17 [NASA ADS] [Google Scholar]
  23. GRAVITY Collaboration (Widmann, F., et al.) 2024, A&A, 681, A115 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  24. Harris, C. R., Millman, K. J., van der Walt, S. J., et al. 2020, Nature, 585, 357 [NASA ADS] [CrossRef] [Google Scholar]
  25. Haubois, X., Abuter, R., Aller-Carpentier, E., et al. 2020, SPIE Conf. Ser., 11446, 1144606 [NASA ADS] [Google Scholar]
  26. Haubois, X., Mérand, A., Abuter, R., et al. 2022, in Optical and Infrared Interferometry and Imaging VIII, eds. A. Mérand, S. Sallum, & J. Sanchez-Bermudez, SPIE Conf. Ser., 12183, 1218306 [NASA ADS] [Google Scholar]
  27. Hunter, J. D. 2007, Computing in Science& Engineering, 9, 90 [NASA ADS] [CrossRef] [Google Scholar]
  28. Lacour, S., Wang, J. J., Nowak, M., et al. 2020, in Optical and Infrared Interferometry and Imaging VII, eds. P. G. Tuthill, A. Mérand, & S. Sallum, SPIE Conf. Ser., 11446, 114460O [NASA ADS] [Google Scholar]
  29. Lutovinov, A., Suleimanov, V., Manuel Luna, G. J., et al. 2020, New Astron. Rev., 91, 101547 [Google Scholar]
  30. Mahadevan, R. 1997, ApJ, 477, 585 [NASA ADS] [CrossRef] [Google Scholar]
  31. Mérand, A. 2022, in Optical and Infrared Interferometry and Imaging VIII, eds. A. Mérand, S. Sallum, & J. Sanchez-Bermudez, SPIE Conf. Ser., 12183, 121831N [Google Scholar]
  32. Mérand, A., Abuter, R., Aller-Carpentier, E., et al. 2014, in Optical and Infrared Interferometry IV, Proc. SPIE, 9146, 91460J [Google Scholar]
  33. Ochsenbein, F., Bauer, P., & Marcout, J. 2000, A&AS, 143, 23 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  34. Pesce, D. W., Palumbo, D. C. M., Narayan, R., et al. 2021, ApJ, 923, 260 [NASA ADS] [CrossRef] [Google Scholar]
  35. Rodriguez, A. C., Cendes, Y., El-Badry, K., & Berger, E. 2024, PASP, 136, 024203 [NASA ADS] [CrossRef] [Google Scholar]
  36. Sjouwerman, L. O., & Blanchard, J. M. 2024a, The Astronomer’s Telegram, 16832, 1 [Google Scholar]
  37. Sjouwerman, L. O., & Blanchard, J. M. 2024b, The Astronomer’s Telegram, 16677, 1 [Google Scholar]
  38. Skrutskie, M. F., Cutri, R. M., Stiening, R., et al. 2006, AJ, 131, 1163 [NASA ADS] [CrossRef] [Google Scholar]
  39. Virtanen, P., Gommers, R., Oliphant, T. E., et al. 2020, Nature Methods, 17, 261 [CrossRef] [Google Scholar]
  40. Wilkins, G. A. 1990, Transactions of the International Astronomical Union Series B, 20 [Google Scholar]
  41. Woillez, J. 2024, in Optical and Infrared Interferometry and Imaging IX, eds. J. Kammerer, S. Sallum, & J. Sanchez-Bermudez, International Society for Optics and Photonics (SPIE), 13095, 130950A [NASA ADS] [Google Scholar]
  42. Woillez, J., Darré, P., Egner, S., et al. 2018, in Optical and Infrared Interferometry and Imaging VI, eds. M. J. Creech-Eakman, P. G. Tuthill, & A. Mérand, SPIE Conf. Ser., 10701, 1070103 [NASA ADS] [Google Scholar]

Appendix A: GRAVITY data quality and contrast limits

We here discuss the interferometric observables of the GRAVITY data collected on 3 June, 23 June, and 18 July 2024 For all epochs, the calibrator BD+08 4132 was positioned on the optical axis of the fringe tracker (FT) and SC fibers, with the same spectral resolution and polarization setting as Gaia BH3.

A.1. 3 June 2024

The 3 June 2024 observations included a splitting of polarizations using a Wollaston prism inserted in the FT and SC beams (GRAVITY Collaboration 2024). Both the SC and FT fibers were initially positioned on the giant star component of Gaia BH3 (single-field, on-axis mode). Unfortunately, a first problem that we detected is that the input fiber of the FT and SC corresponding to the UT1 telescope drifted with respect to the target star during the observation. This resulted in a progressive loss of flux, and also in a displacement in phase of the fringes during the sequence. These symptoms are observed on all baselines that include UT1, and can be explained by an inaccurate sky tracking. We attribute the variation of the differential phase with wavelength to this drift. The drift in optical path difference (OPD) translates to a different phase offset depending on wavelength. Moreover, we see a differential behavior in polarization (blue and green curves in the U1-U3 panel of the top figure in Fig. A.1). This is likely also a symptom of the misalignment of the fibers on the target. A second problem is that the visibility is increasing with wavelength in the upper panel of Fig. A.1, while the star being unresolved we expect a constant value of 1.0 (as in the bottom panel of Fig. A.2). This indicates that the FT and SC channels were out of phase. The metrology was inactive during these observations, and although it is in principle not needed for single-field observations, this apparently caused an OPD offset between the FT (zero OPD) and the SC (non-zero OPD). Considering the biases that are created by these instrumental problems, we did not include the data set of 3 June 2024 is our search for BH emission (however, we present the corresponding radial contrast limit in App. A.4).

thumbnail Fig. A.1.

Interferometric observables collected on 3 June 2024 in medium spectral resolution with split polarizations (the blue and green points represent the two polarizations). In each panel, the three plots show the CPs (left), the spectral differential phases ϕ (center), and the squared visibilities V2 (right) measured with GRAVITY on Gaia BH3.

A.2. 23 June 2024

On 23 June 2024, we tested both the medium and high spectral settings of the SC, respectively providing spectral resolutions of R ≈ 500 and 4000, while stabilizing the fringes with the FT. The two sets of fibers (FT and SC) were on-axis, that is, both positioned on the giant star. As listed in Table 2, the first ’object’ exposure on Gaia BH3 in high spectral resolution mode starting at UT08:30:14 is affected by a large visibility bias and is discarded.

The feature present in Fig. A.2 around 2.16 μm is intriguing, as it could in principle indicate Brγ line emission from the BH. However, this feature is unlikely to be of astrophysical origin for two reasons. Firstly, no comparable signature is observed in CP. This is a consistency check as the differential phase amplitude reaches 5 to 10°, and it should be detectable in CP if caused by BH emission. Secondly, there is no apparent signal in the squared visibility V2 at 2.16 μm, while a secondary source creating such a differential phase offset should induce a detectable signature.

The CP and V2 signals that are used for companion search, however, behave correctly, with satisfactory quality control parameters. We therefore included the 23 June 2024 data in the estimate of the BH emission contrast limit in Sect. 3.2.

thumbnail Fig. A.2.

Same quantities as in Fig. A.1 but for the observation of 23 June 2024 in medium spectral resolution with combined polarizations.

thumbnail Fig. A.3.

Same quantities as in Fig. A.1 but for the observation of 23 June 2024 in high spectral resolution with combined polarizations.

A.3. 18 July 2024

For our last observation, the SC fiber was offset from the FT fiber (dual-field, on-axis mode) by the expected star-BH relative position to optimize the injection of flux coming from the BH position in the SC. Due to the 18 mas astrometric offset of the SC fibers with respect to the star’s position, a wavelength-dependent phase shift of the SC fringes with respect to the FT fringes is visible in the central panel of the lower part of Fig. A.4. The large oscillations in differential phase ϕ are caused by phase folding. The V2 and CP signals are consistent with the differential phases, and the photometric S/N level was high and stable during the whole observing sequence. The standard quality control parameters of this data set are all satisfactory. The data obtained on 18 July 2024 is therefore included together with the 23 June 2024 epoch in our contrast limit determination presented in Sect. 3.2.

thumbnail Fig. A.4.

Interferometric observables collected on 18 July 2024 in medium spectral resolution with combined polarizations. This observation has been obtained with the SC fiber centered on the predicted BH position. The astrometric offset with respect to the star’s position explains the large oscillations in differential phase ϕ as being due to phase folding. The plotted quantities are the same as in Fig. A.1.

A.4. Contrast limits

Fig. A.5 (left and central panels) shows a PMOIRED simulation of the interferometric CP (T3PHI) and normalized visibility (N|V|) signals that would be produced by a 5% emission from the expected location of the BH (red curves). The normalized visibility is defined as N|V|=V(λ)/L[V(λ)] where L[V(λ)] is a linear fit to V(λ) as a function of λ. The 18 July 2024 data set displayed as gray curves do not exhibit any signature of such a flux contribution. The radial average dependence of the contrast limit established at each epoch from the medium spectral resolution observations is presented in the right panel of Fig. A.5. The improvement in sensitivity for the 23 June and 18 July 2024 compared to the 3 June 2024 epoch is clearly visible. At the expected separation of the BH of 18 mas in Fig. A.5 (right panel), the radial average contrast is slightly shallower than the δm = 6.8 mag value reported in Sect. 3.2. This is due to the fact that the particular (u, v) plane coverage of the observations results in a deeper local sensitivity around the expected location of the BH (Fig. 2) than its radial average.

thumbnail Fig. A.5.

Left and central panels: Simulated CPs (T3PHI) and normalized visibility (N|V|) for a 5% flux companion (red curves) injected in the data collected on 18 July 2024 (gray curves). Right panel: Contrast limits from CANDID as a function of the distance from the star for the three epochs in medium resolution. The contrast with 90% confidence is shown in magnitude difference (left axis) and contrast (right axis).

A.5. High spectral resolution data and contrast limit

The CP and triple amplitude signals as a function of wavelength in the region of the Brγ line are shown in Fig. A.6 for the 23 June 2024 epoch (high and medium resolution) and in Fig. A.7 for the 18 July 2024 epoch (medium resolution). The upper limit on the Brγ emission from the BH based on the high spectral resolution data is shown in Fig. A.8 and discussed in Sect. 3.3.

thumbnail Fig. A.6.

GRAVITY CPs (top panel) and triple amplitudes (bottom panel) of Gaia BH3 measured on 24 June 2024 (high and medium spectral resolutions) around the hydrogen Brγ line (rest wavelength λ = 2.166 μm, Doppler-shifted in Gaia BH3 to λ = 2.164 μm).

thumbnail Fig. A.7.

Closure phases (top panel) and triple amplitudes (bottom panel) of Gaia BH3 in the region of the hydrogen Brγ line measured with GRAVITY on 18 July 2024 in medium spectral resolution mode with combined polarizations.

thumbnail Fig. A.8.

Upper limit of the equivalent width (left scale) and contrast limit (right scale) in the Brackett γ line emission at the astrometric position of the BH component of Gaia BH3 as a function of the assumed emission line velocity dispersion.

All Tables

Table 1.

Summary of the properties of the giant star in LS II +14 13 (Gaia BH3) and the interferometric calibrator BD+08 4132.

Table 2.

Log of the observations of Gaia BH3 (LS II +14 13) and the interferometric calibrator (BD+08 4132) with GRAVITY.

All Figures

thumbnail Fig. 1.

Position of Gaia BH3 and its companion star relative to their center of mass at the epoch of GRAVITY observations. The position of the Lagrangian L1 point of the system is shown with a purple diamond symbol.

In the text
thumbnail Fig. 2.

Top panel: CANDID map of companion detection probability expressed in the number of standard deviations of the noise level for the combination of 23 June and 18 July 2024 data in medium spectral resolution. The S/N color scale is from 0 to 2.1σ. The position of the star is marked with a star symbol, and the maximum S/N ratio of 2.3 is indicated with a red cross. The predicted position of the BH is framed with a white square. Bottom panel: Contrast map for a 3σ detection level around the expected position of the BH, where the color scale represents the contrast Δm in magnitudes.

In the text
thumbnail Fig. 3.

Left panel: PMOIRED map of the limiting magnitude (3σ, 99% confidence) of the BH around its predicted position for the 23 June and 18 July 2024 combined data set. The displayed field corresponds to the white square in Fig. 2. Right panel: statistical distribution of the magnitudes of the recovered 3σ companions. Here, the dashed line represents the CANDID detection limit.

In the text
thumbnail Fig. 4.

Upper limits on the spectral energy distribution of Gaia BH3 established from radio (August–September 2024; Sjouwerman & Blanchard 2024a,b) and X-ray (May 2024; Cappelluti et al. 2024) observations and the GRAVITY near-infrared observations (June–July 2024). The blue curve represents the accreting BH model presented by Cappelluti et al. (2024) for fEdd < 4.9 × 10−7.

In the text
thumbnail Fig. A.1.

Interferometric observables collected on 3 June 2024 in medium spectral resolution with split polarizations (the blue and green points represent the two polarizations). In each panel, the three plots show the CPs (left), the spectral differential phases ϕ (center), and the squared visibilities V2 (right) measured with GRAVITY on Gaia BH3.

In the text
thumbnail Fig. A.2.

Same quantities as in Fig. A.1 but for the observation of 23 June 2024 in medium spectral resolution with combined polarizations.

In the text
thumbnail Fig. A.3.

Same quantities as in Fig. A.1 but for the observation of 23 June 2024 in high spectral resolution with combined polarizations.

In the text
thumbnail Fig. A.4.

Interferometric observables collected on 18 July 2024 in medium spectral resolution with combined polarizations. This observation has been obtained with the SC fiber centered on the predicted BH position. The astrometric offset with respect to the star’s position explains the large oscillations in differential phase ϕ as being due to phase folding. The plotted quantities are the same as in Fig. A.1.

In the text
thumbnail Fig. A.5.

Left and central panels: Simulated CPs (T3PHI) and normalized visibility (N|V|) for a 5% flux companion (red curves) injected in the data collected on 18 July 2024 (gray curves). Right panel: Contrast limits from CANDID as a function of the distance from the star for the three epochs in medium resolution. The contrast with 90% confidence is shown in magnitude difference (left axis) and contrast (right axis).

In the text
thumbnail Fig. A.6.

GRAVITY CPs (top panel) and triple amplitudes (bottom panel) of Gaia BH3 measured on 24 June 2024 (high and medium spectral resolutions) around the hydrogen Brγ line (rest wavelength λ = 2.166 μm, Doppler-shifted in Gaia BH3 to λ = 2.164 μm).

In the text
thumbnail Fig. A.7.

Closure phases (top panel) and triple amplitudes (bottom panel) of Gaia BH3 in the region of the hydrogen Brγ line measured with GRAVITY on 18 July 2024 in medium spectral resolution mode with combined polarizations.

In the text
thumbnail Fig. A.8.

Upper limit of the equivalent width (left scale) and contrast limit (right scale) in the Brackett γ line emission at the astrometric position of the BH component of Gaia BH3 as a function of the assumed emission line velocity dispersion.

In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.