Polarization-dependent observational signatures of Weyl-coupled photons around a black hole

Technical Summary: Polarization-dependent observational signatures of Weyl-coupled photons around a black hole

Problem Statement
The paper investigates the observational consequences of a non-minimal coupling between photons and the Weyl tensor in a curved spacetime, specifically around a Schwarzschild black hole. In this model, defined by an action term $\alpha C_{\mu\nu\rho\sigma}F^{\mu\nu}F^{\rho\sigma}$, the vacuum acts as a birefringent medium where the two linear polarization states (labeled PPL and PPM) propagate according to distinct optical metrics. While previous studies have addressed the geometric-optics limit (lensing, shadows) or specific wave-optics aspects separately, a unified, numerically controlled treatment connecting the partial-wave scattering matrix to the full range of observables (from absorption to photon rings) for both polarizations has been lacking. The primary goal is to quantify how the coupling constant $\alpha$ splits the photon sphere, critical impact parameters, and scattering observables, providing a framework to constrain $\alpha$ using both wave and geometric optics.

Methodology
The authors employ a dual approach combining wave optics and geometric optics on a Schwarzschild background ($f(r) = 1 - 2M/r$).

1. Wave Optics (Scattering Matrix):

   • The field equations are reduced to two decoupled Schrödinger-type radial equations for the two parity sectors (even/odd), characterized by polarization-dependent effective potentials $V^\pm_\ell(r)$.
   • The scattering matrix elements $S^\pm_\ell$ are computed using a Riccati–Hankel matching scheme. Solutions are integrated from the horizon to a matching radius using the Numerov method and matched to Riccati–Hankel functions at infinity. This specific basis is chosen to preserve the phase of the centrifugal tail, which is critical for phase-sensitive observables like the differential cross-section and glory.
   • Observables derived include greybody factors, absorption cross-sections, and differential cross-sections (separated into helicity-preserving $f$ and helicity-reversing $g$ amplitudes).
   • Numerical accuracy is validated against unitarity, free-potential limits, and literature benchmarks (e.g., Crispino et al. for the Schwarzschild glory).

2. Geometric Optics:

   • The high-frequency limit is analyzed via null geodesics on polarization-dependent optical metrics.
   • Key quantities calculated include the photon sphere radius $r_{ph}$, the critical impact parameter $b_c$, and the Lyapunov exponent $\gamma$ governing the decay of photon ring subrings.
   • Shadow and Ring Imaging: The authors utilize a four-color celestial-sphere renderer to visualize the "double shadow." The results are cross-validated using three independent ray tracers: an adaptive Runge-Kutta integrator, a GPU-based tracer (Mahakala style), and the Gyoto code.

Key Contributions and Results

• Unified Framework: The paper provides the first numerically controlled treatment that carries both polarization states through the entire chain of observables, from the scattering matrix to the shadow and photon ring.
• Polarization Splitting: The coupling $\alpha$ splits the critical impact parameter $b_c$ and the photon sphere radius.
  • For $\alpha/M^2 = 0.75$, the PPM (odd parity) branch contracts ($b_c \approx 3.22 M$) while the PPL (even parity) branch expands ($b_c \approx 6.46 M$).
  • The separation between the two shadow edges reaches 62% of the shadow radius at this coupling strength.
• Double Shadow and Photon Rings: The splitting manifests as a "double shadow" where the black disk represents the region where both polarizations are captured, surrounded by a colored ring where only one polarization is captured. Similarly, the photon ring splits into two interlaced systems with different Lyapunov exponents; the PPM ring decays much faster (fainter subrings) than the PPL ring.
• Wave-Optics Signatures:
  • Absorption: The absorption cross-sections oscillate around the geometric capture areas $\pi (b_c^\pm)^2$. The thresholds for different partial waves shift in opposite directions for the two polarizations.
  • Differential Scattering and Glory: The glory (backward scattering) patterns split, with fringe spacing determined by the respective $b_c^\pm$.
  • Backward Birefringence: A key finding is the behavior at the exact backward angle ($\theta = \pi$). In the minimal theory ($\alpha=0$), the backward flux is identically zero due to parity symmetry (helicity-reversing amplitude vanishes). For $\alpha \neq 0$, this symmetry is broken, generating a non-zero backward signal proportional to $|S_M - S_E|^2$. The interference nulls in this signal provide a direct interferometric readout of the splitting $\Delta b_c$.
• Stability Constraints: The authors note that dynamical stability restricts the admissible range of $\alpha$ to $-M^2/2 < \alpha < M^2$, and the values displayed in the paper lie within this safe window.

Significance and Claims
The paper claims that the backward birefringence signal is the "cleanest diagnostic" for this coupling. Because it vanishes identically in the standard Schwarzschild case and is protected by parity, it is free from the astrophysical foregrounds and degeneracies (such as spin or inclination) that plague total-intensity imaging.

However, the authors maintain a modest stance regarding immediate astrophysical detection:

• They explicitly state that for astrophysical black holes (e.g., M87*, Sgr A*), the QED-induced value of $\alpha$ is far too small ($\sim 10^{-35}$ to $10^{-55}$) to be observed.
• The analysis is presented as a phenomenological template and a "proof of principle" for constraining $\alpha$ as a generic effective field theory parameter.
• The paper argues that while current Event Horizon Telescope (EHT) precision is insufficient to resolve the predicted edge separations (requiring $\sim 8\%$ resolution for $\alpha/M^2 \sim 0.1$), the difference observables (polarized edge separation or backward scattering) offer the most promising route to breaking degeneracies in future next-generation horizon-scale polarimetry (ngEHT, space VLBI).

In summary, the work establishes a rigorous computational link between Weyl-coupled photon dynamics and multi-messenger observables, identifying specific polarization-dependent signatures that could, in principle, constrain non-minimal gravity-photon interactions if the coupling were phenomenologically larger than standard QED predictions.