The Unanticipated Phenomenology of the Blazar PKS 2131-021: A Unique Supermassive Black Hole Binary Candidate

Most large galaxies host supermassive black holes in their nuclei and are subject to mergers, which can produce a supermassive black hole binary (SMBHB), and hence periodic signatures due to orbital motion. We report unique periodic radio flux density variations in the blazar PKS 2131-021, which strongly suggest an SMBHB with an orbital separation of similar to 0.001-0.01 pc. Our 45.1 yr radio light curve shows two epochs of strong sinusoidal variation with the same period and phase to within less than or similar to 2% and similar to 10%, respectively, straddling a 20 yr period when this variation was absent. Our simulated light curves accurately reproduce the "red noise" of this object, and Lomb-Scargle, weighted wavelet Z-transform and least-squares sine-wave analyses demonstrate conclusively, at the 4.6 sigma significance level, that the periodicity in this object is not due to random fluctuations in flux density. The observed period translates to 2.082 +/- 0.003 yr in the rest frame at the z = 1.285 redshift of PKS 2131-021. The periodic variation in PKS 2131-021 is remarkably sinusoidal. We present a model in which orbital motion, combined with the strong Doppler boosting of the approaching relativistic jet, produces a sine-wave modulation in the flux density that easily fits the observations. Given the rapidly developing field of gravitational-wave experiments with pulsar timing arrays, closer counterparts to PKS 2131-021 and searches using the techniques we have developed are strongly motivated. These results constitute a compelling demonstration that the phenomenology, not the theory, must provide the lead in this field.