Relativistic Jets and Winds in Radio-Identified Supermassive Black Hole Binary Candidates

TL;DR

The work addresses identifying individual SMBHBs as sources of the nanohertz gravitational-wave background by presenting two radio-identified blazar candidates, PKS 2131-021 and PKS J0805-0111, with persistent multi-wavelength periodicity. It introduces a semi-analytic disk-wind–collimated jet model in which a sub-relativistic wind channels an ultra-relativistic jet, and orbital motion induces a conical helical jet that modulates flux through Doppler aberration; phase delays across frequencies arise from the wind‑driven geometry and optical depth. The model reproduces the observed light curves and phase shifts for both sources, predicts polarization and direct-imaging signatures, and forecasts emission-line variability, thereby offering concrete tests with future observations and simulations (e.g., RMHD, VLBI, and PTA campaigns). The results imply high jet powers and a wind–jet boundary speed around β_w ≈ 0.9, with similar parameters across the two candidates, supporting a common physical picture for jetted SMBHBs and their role as nHz GW sources with identifiable electromagnetic counterparts.

Abstract

Supermassive black hole binary systems (SMBHBs) are thought to emit the recently discovered nHz gravitational wave background; however, not a single individual nHz source has been confirmed to date. Long-term radio-monitoring at the Owens Valley Radio Observatory has revealed two potential SMBHB candidates: blazars PKS 2131-021 and PKS J0805-0111. These sources show periodic flux density variations across the electromagnetic spectrum, signaling the presence of a good clock. To explain the emission, we propose a generalizable jet model, where a mildly relativistic wind creates an outward-moving helical channel, along which the ultra-relativistic jet propagates. The observed flux variation from the jet is mostly due to aberration. The emission at lower frequency arises at larger radius and its variation is consequently delayed, as observed. Our model reproduces the main observable features of both sources and can be applied to other sources as they are discovered. We make predictions for radio polarization, direct imaging, and emission line variation, which can be tested with forthcoming observations. Our results motivate future numerical simulations of jetted SMBHB systems and have implications for the fueling, structure, and evolution of blazar jets.

Figures (8)

• Figure 1: A visualization of our model. (Top) An SMBHB composed of the primary black hole H and a secondary black hole S. H orbits with speed $v_H$, carries around a disk (D) and launches a jet (J) along the orbital angular momentum direction (L). The disk produces a wind (W) with speed $v_w<c$, that collimates the jet. (Middle) A geometric illustration of the conical jet along the $\hat{z}$ axis with opening angle $\theta_j$ confined by the wind. The wind-jet interface has velocity $\vec{\beta}_w=\vec{v}_w/c$. We show a unit vector $\hat{r}$ associated with angular position ($\chi$, $\psi$) as well as the observer viewing direction $\hat{n}$. The orbital motion with velocity $\vec{\beta}_H=\vec{v}_H/c$ will distort the jet into a conical helix. (Bottom) A 2D slice through the wind-collimated helical jet. The dashed line in the center denotes the sinuous path of relativistic flow moving through the center of the channel. The minimum radii $r_{\min}$ which contribute at two frequencies are labeled.
• Figure 2: The local fractional variation in Doppler factor $\delta D/D$ (top), PD $\delta \Pi/\Pi_0$ (middle), and EVPA $\delta$PA in degrees (bottom) for choices of $\beta_w$ at a position $(\chi, \psi)=(0,0)$ within the jet. We fix $\beta_H=0.02$, $\Gamma_j=10$, and $i=3.8^\circ$.
• Figure 3: PKS 2131$-$021 light curves at radio, mm, sub-mm, and optical wavelengths. (a) Haystack 2.7--31.4 GHz observations 1986AJ.....92.1262O together with the least-squares sine-wave fits to the data at each of the frequencies. (b) OVRO 15 GHz light curve plus ALMA light curves at 91.5, 104, and 345 GHz. (c) Comparison of the $i$, $r$, and $g$-band ZTF optical and OVRO radio light curves of PKS 2131. The curves show the least squares sine wave fits to the corresponding data. (d) The observed phase shifts, relative to the OVRO 15 GHz light curve. A positive phase shift indicates that the light curve is shifted to a later time than the 15 GHz light curve. The curve shows a quadratic polynomial fit to the phase offset. Reproduced from Fig. 6 in 2024Paper2.
• Figure 4: The emissivity as a function of radius $j_{\nu}(r)$ at three frequencies. The emissivity is assumed to be non-zero between $r_{\min}(\nu)$ to $r_{\max}(\nu)$. See text for details.
• Figure 5: Our model for PKS 2131. (Top) The predicted intensity map at 230 GHz. (Bottom left) The model flux density in a selection of observed bands compared with OVRO and ALMA light curves shifted down by a constant value to match the model. (Bottom Right) Model phase shifts relative to the 15 GHz light compared to the observed phase shifts. Parameters are $\beta_w=0.9$, $\beta_H=0.02$, $\theta_j=0.1$, $\Gamma_j=10$, $r_{\min}$ chosen so that $\tau_\nu(r_{\min})=20$, $r_{\max}=r_{\min}+60 \,(\nu (1+z_S) /15 \text{ GHz})^{-1}$ ly and $i=3.8^\circ$. We assume $L_e=10^{46} \text{ erg s}^{-1}\, (r/r_0)/(1+(r/r_0))$ (where $r_0=5$ ly) and $L_m=10^{46} \text{ erg s}^{-1}\, 1/(1+(r/r_0))$.