Gravitational wave signatures from periodic orbits around a Schwarzschild-Bertotti-Robinson black hole

In this paper, we investigate periodic bound orbits and gravitational wave (GW) emission in the Schwarzschild-Bertotti-Robinson (Schwarzschild-BR) spacetime-an exact electrovacuum solution describing a static black hole (BH) immersed in a uniform magnetic field. We explore how the background magnetic field qualitatively alters the BH’s gravitational dynamics, affecting timelike geodesics such as the marginally bound orbit (MBO) and the innermost stable circular orbit (ISCO). We then analyze periodic bound orbits using the frequency ratio ${ω_φ}/{ω_{r}}$, which characterizes the orbits by their azimuthal and radial motions. Based on the numerical kludge method we further compute the gravitational waveforms emitted from periodic orbits around a supermassive Schwarzschild-BR BH. We show that the background magnetic field significantly changes orbital frequencies, resonance conditions, zoom-whirl structures, and the resulting waveforms. Finally, we examine the frequency spectra in the mHz range and the detectability of these GW signals by computing the characteristic strain via a discrete Fourier transform on the time-domain waveforms, comparing the results with the sensitivity curves of space-based GW detectors such as LISA, Taiji, and TianQin. Our results show that intrinsically magnetic fields modify spacetime and leave observable imprints on extreme mass-ratio inspiral GWs, which may be tested by future observations.

💡 Research Summary

This paper presents a comprehensive study of periodic bound orbits and the associated gravitational‑wave (GW) emission in the Schwarzschild‑Bertotti‑Robinson (Schwarzschild‑BR) spacetime, an exact electrovacuum solution describing a non‑rotating black hole immersed in a uniform magnetic field. The authors first introduce the metric, which depends on a dimensionless magnetic‑field parameter B, and derive the conserved energy E and angular momentum L for a neutral test particle moving in the equatorial plane. By constructing the effective potential V_eff, they show that increasing B pushes the location of stable circular orbits inward and raises the potential barrier. Analytic expressions for the innermost stable circular orbit (ISCO) radius, angular momentum, and energy are obtained (Eqs. 10‑12), revealing a clear B‑dependence: r_ISCO = 6M/(1 − B²M²), etc. Numerical scans of the marginally bound orbit (MBO) confirm that as B grows, r_MBO and r_ISCO increase while L_MBO decreases, a trend summarized in Table I.

The core of the work focuses on periodic orbits, defined by a rational ratio of azimuthal to radial frequencies, q = ω_φ/ω_r − 1 = w + v/z, where (z,w,v) are the zoom, whirl, and vertex integers. Using the geodesic equations, the authors compute q as a function of (E, L) for several B values (Fig. 3). They find that for a given B, q rises sharply as E approaches the MBO limit and that larger B shifts the q‑versus‑E curve upward, meaning higher energies are required to achieve the same resonance. Conversely, q exhibits a peak near the minimum allowed L and then declines for larger L, with stronger magnetic fields moving the peak to lower L. Table II lists specific (E, L) pairs for selected (z,w,v) combinations, and Figures 4–5 illustrate the resulting orbital shapes for B = 0.02, highlighting how increasing z produces more intricate zoom‑whirl patterns while larger w adds extra revolutions between successive apastra.

To assess the GW signatures, the authors adopt the “numerical kludge” approach. They first solve the geodesic equations numerically to obtain the particle trajectory, then feed this trajectory into the quadrupole formula h_{ij}=2G/(c⁴D_L) · ¨I_{ij} (Eq. 15) to generate approximate waveforms. This method assumes the adiabatic regime of EMRIs, where radiation reaction is slow compared to the orbital period, allowing the motion to be treated as a sequence of geodesics. Waveforms are computed for B = 0, 0.02, 0.04, and the authors observe that stronger magnetic fields suppress high‑frequency components associated with the whirl phase and shift the overall spectral power toward lower frequencies in the millihertz band. By performing a discrete Fourier transform, they calculate the characteristic strain h_c(f)=2f|ĥ(f)| and compare it against the sensitivity curves of LISA, Taiji, and TianQin. Even for B ≈ 0.04, the strain remains above LISA’s noise floor, especially for resonant orbits with (z,w,v) = (3–4, 2–3, *), where the signal‑to‑noise ratio is markedly enhanced.

The paper concludes that (i) the Schwarzschild‑BR spacetime provides a clean, exact laboratory for studying the interplay between magnetic fields and strong‑gravity dynamics; (ii) the background magnetic field significantly modifies key orbital parameters (ISCO, MBO), the rational frequency ratio, and consequently the zoom‑whirl structure; (iii) numerical‑kludge waveforms capture these effects efficiently, offering a practical tool for future data‑analysis pipelines; and (iv) the resulting GW signals lie comfortably within the detection band of upcoming space‑based interferometers, opening the possibility of probing intrinsic magnetic fields around supermassive black holes through EMRI observations. The authors suggest extensions to the rotating Kerr‑BR case and inclusion of radiation‑reaction corrections as natural next steps.