The quest for binary and dual supermassive black holes (SMBHs) at the dawn of the multi-messenger era is compelling. Detecting dual active galactic nuclei (AGN)-active SMBHs at projected separations larger than several parsecs-and binary AGN-probing the scale where SMBHs are bound in a Keplerian binary-is an observational challenge. The study of AGN pairs (either dual or binary) also represents an overarching theoretical problem in cosmology and astrophysics. The AGN triggering calls for detailed knowledge of the hydrodynamical conditions of gas in the imminent surroundings of the SMBHs and, at the same time, their duality calls for detailed knowledge on how galaxies assemble through major and minor mergers and grow fed by matter along the filaments of the cosmic web. This review describes the techniques used across the electromagnetic spectrum to detect dual and binary AGN candidates and proposes new avenues for their search. The current observational status is compared with the state-of-the-art numerical simulations and models for formation of dual and binary AGN. Binary SMBHs are among the loudest sources of gravitational waves (GWs) in the Universe. The search for a background of GWs at nHz frequencies from inspiralling SMBHs at low redshifts, and the direct detection of signals from their coalescence by the Laser Interferometer Space Antenna in the next decade, make this a theme of major interest for multi-messenger astrophysics. This review discusses the future facilities and observational strategies that are likely to significantly advance this fascinating field.

Very Long Baseline Interferometry (VLBI) allows for high-resolution and highsensitivity observations of relativistic jets, that can reveal periodicities of several years in their structure. We perform an analysis of long-term VLBI data of the quasar S5 1928+738 in terms of a geometric model of a helical structure projected onto the plane of the sky. We monitor the direction of the jet axis through its inclination and position angles. We decompose the variation of the inclination of the inner 2 milliarcseconds of the jet of S5 1928+738 into a periodic term with amplitude of ∼ 0.89 • and a linear decreasing trend with rate of ∼ 0.05 • yr −1 . We also decompose the variation of the position angle into a periodic term with amplitude of ∼ 3.39 • and a linear increasing trend with rate of ∼ 0.24 • yr −1 . We interpret the periodic components as arising from the orbital motion of a binary black hole inspiraling at the jet base and derive corrected values of the mass ratio and separation from the accumulated 18 years of VLBI data. Then we identify the linear trends in the variations as due to the slow reorientation of the spin of the jet emitter black hole induced by the spin-orbit precession and we determine the precession period T SO = 4852 ± 646 yr of the more massive black hole, acting as the jet emitter. Our study provides indications, for the first time from VLBI jet kinematics, for the spinning nature of the jet-emitting black hole.