Searching for Nearby Supermassive Black Hole Binaries with GRAVITY

We present a systematic interferometric search for close supermassive black hole (SMBH) binaries in 20 nearby type 1 active galactic nuclei (AGN) observed with VLTI/. We investigate parsec-scale asymmetric emission by modeling the squared visibilities and closure phases with a single circular Gaussian component plus a background point source for the central AGN, and searching for additional compact companions through grid-based binary models. The robustness of the search is validated using injection-recovery simulations and point-spread-function diagnostics. Several features are identified in the initial significance maps, but only one candidate out of 20 shows a persistent compact asymmetry, corresponding to an observed incidence of 5%. If NGC 3783 is interpreted as a binary AGN, we derive a Poissonian one-sided 1σ upper limit of f bin<16.5% for GRAVITY-detectable binaries among bright nearby AGN (K<12.5) at separations of about1 - 100 pc. Conversely, if the observed asymmetry in NGC 3783 is associated with asymmetric circumnuclear dust emission, the corresponding upper limit on the binary fraction is f bin<9.2%. We compare these results with predictions from the Horizon-AGN simulation, finding comparable binary fractions. We also noted that the expected fraction strongly depends on the assumed distribution of individual SMBH accretion rates within binaries, which is currently unresolved in cosmological models. Our results demonstrate the potential of infrared interferometry to identify parsec-scale binary AGN candidates and motivate future GRAVITY+ observations to expand the sample and follow up known candidates, helping to confirm their nature and distinguish genuine binaries from asymmetric nuclear structures.