JWST Confirmation of a Runaway Supermassive Black Hole via its Supersonic Bow Shock

TL;DR

This study confirms RBH-1 as a runaway SMBH at z=0.96 by revealing a kiloparsec-scale supersonic bow shock at the wake tip through JWST/NIRSpec IFU data. A simple bow-shock model ties a velocity of $v_{BH} \approx 954$ km s$^{-1}$ and an inclination $i \approx 29^{\circ}$ to the observed ~600 km s$^{-1}$ gradient over ~1 kpc, with line-ratio diagnostics from fast radiative shocks corroborating the shock velocity and preshock conditions. Energy conservation arguments yield a SMBH mass $M_{BH} \gtrsim \text{a few} \times 10^{7}$ M$_{\odot}$, consistent with bulge-dynamical expectations and supporting a GW recoil or multi-body ejection origin. The wake results from turbulent entrainment and CGM mixing, explaining the downstream kinematics and the wake’s longevity, while UV data offer tentative localization of the SMBH vicinity. Overall, the work provides compelling empirical validation for runaway SMBHs and highlights wakes as a morphological pathway to identify them in future wide-field surveys.

Abstract

We present JWST/NIRSpec IFU observations of a candidate runaway supermassive black hole at the tip of a 62 kpc-long linear feature at z=0.96. The JWST data show a sharp kinematic discontinuity at the tip, with a radial velocity change of $\approx 600$ km/s across 0.1'' (1 kpc). The velocity gradient, together with the projected post-shock flow velocity of $\approx 300$ km/s, is well described by a simple shock-compression model of a supersonic object, with a velocity of $v_{BH} = 954^{+110}_{-126}$ km/s and an inclination $i=29^{+6}_{-3}$ deg. The previously puzzling kinematics along the linear feature, with the observed radial velocity decreasing from $\approx 300$ km/s near the tip to $\approx 100$ km/s closer to the former host galaxy, are naturally explained as gradual downstream mixing of shocked gas with the circumgalactic medium through turbulent entrainment. The runaway black hole interpretation is further supported by the morphology of the gas at the tip of the wake and an analysis of the [OIII]/H$α$, [NII]/H$α$, [SII]/H$α$, and [SIII]/[SII] line ratios. The line ratios are consistent with fast radiative shocks and rapid cooling, with best-fit shock velocities that are in good agreement with expectations from the black hole velocity and the shock geometry. Energy conservation over the lifetime of the wake suggests a SMBH mass of $M_{BH} \gtrsim 10^7$ M$_{\odot}$. These results confirm that the wake is powered by a supersonic runaway supermassive black hole, a long-predicted consequence of gravitational-wave recoil or multi-body ejection from galactic nuclei.

Figures (13)

• Figure 1: Top: HST/WFC3 UVIS imaging of the candidate $z=0.96$ runaway black hole wake RBH-1. The displayed image spans $13\farcs 2 \times 7\farcs 7$. It is the sum of two long pass filters, F200LP and F350LP, with a combined integration time of 29,898 s. The two JWST NIRSpec pointings are indicated in yellow. Bottom: summed flux profile along the feature. There is continuous emission all the way from the galaxy to the tip at $r\approx 62$ kpc, followed by a sudden drop of a factor of $>40$. The profile does not resemble the exponential fall-off expected for an edge-on spiral galaxy sanchez:23.
• Figure 2: Overview of the NIRSpec IFU data. Insets show a map of the redshifted [O III] emission, with three characteristic locations indicated by vertical lines. The main panels show the spectra at those locations. There is a distinct, very bright knot at the tip of RBH-1, with a high [O III]/H$\alpha$ ratio. This ratio decreases gradually behind the tip, while the [N II]/H$\alpha$ and [S II]/H$\alpha$ ratios increase. In all panels [O III] and H$\alpha$ are the strongest lines; these are used for the dynamical analysis in § \ref{['kinematics.sec']}.
• Figure 3: Left: Summed [O III] + H$\alpha$ line emission at the tip of the linear feature. The tip is clearly resolved, and resembles a classic bow shock with a thin turbulent wake behind it. The thick black line is a restricted parabolic fit to a faint isophote at $63\leq r \leq 65$ kpc. Right: Illustration of a bow shock, at approximately the same scale. Upon encountering the shock gas flows tangentially along the shock surface. The low pressure behind the shock causes flow lines to bend toward the symmetry axis, creating a narrow wake.
• Figure 4: Illustration of a flow surface for a shock that is projected at an intermediate angle ($50^{\circ}$) toward the observer. Gas on the near limb is blueshifted and gas on the far limb is redshifted with respect to the wake velocity. The blueshifted gas is projected further downstream than the redshifted gas, leading to a velocity gradient.
• Figure 5: Middle: Wavelength slices through the data cube, near the wavelength of the redshifted [O III] $\lambda 5007$ line (left) and H$\alpha$ (right). Fluxes have been scaled to the peak in each slice. There is a striking pattern, with the emission near the tip systematically shifting further upstream for higher velocities. Left and right: Position-velocity diagrams for the regions in between the white lines, obtained by averaging the flux in the direction perpendicular to the wake. These diagrams show an unambiguous gradient of $\sim 600$ km s$^{-1}$ over $\sim 1$ kpc. The gradient is quantified by the white data points.