Pseudo-spectral methods in gravitational lensing: Spin precession in loop quantum black hole spacetimes

Technical Summary: Pseudo-spectral methods in gravitational lensing: Spin precession in loop quantum black hole spacetimes

Problem Statement
Gravitational lensing is conventionally treated as an initial value problem: given an impact parameter $b$, one computes the resulting deflection angle. This paper advocates an inverse formulation, treating lensing as a two-point boundary value problem (BVP). In this framework, a specific deflection angle (particularly for backward scattering) is prescribed, and the orbit is determined by boundary conditions at spatial infinity. Consequently, the impact parameter emerges not as an input, but as a nonlinear eigenvalue to be solved simultaneously with the orbit. The authors apply this framework to compute the spin precession of massive particles undergoing backward scattering in two specific black hole spacetimes derived from Loop Quantum Gravity (LQG): the covariant effective metric and the self-dual loop quantum black hole.

Methodology
The core methodological contribution is the application of pseudo-spectral (collocation) methods, specifically using Chebyshev-Gauss-Lobatto grids, to solve the orbital and transport equations.

1. Boundary Value Formulation:

   • The orbit equation is derived for a general static, spherically symmetric line element with three free functions ($f, g, h$).
   • For backward scattering with a total azimuthal angle $2\phi_0 = 2n\pi$, the problem is posed on the fixed interval $[0, 2n\pi]$.
   • The system consists of a second-order orbit equation with two boundary conditions ($u=0$ at both ends) and a supplementary condition derived from the first integral of motion at the boundary. This supplementary condition determines the squared impact parameter $\eta^2$, which acts as the nonlinear eigenvalue.
   • The authors exploit the symmetry of the orbit about the turning point to reduce the computational domain to half the interval.

2. Pseudo-spectral Implementation:

   • The continuous differential equations are discretized on a Chebyshev collocation grid. Derivatives are computed via multiplication by differentiation matrices ($D^{(1)}, D^{(2)}$).
   • The orbit equation becomes a nonlinear algebraic system solved via Newton's method with an analytically computed Jacobian.
   • Once the orbit $u(\phi)$ is determined, the parallel transport equations for the spin vector reduce to a single linear system on the same grid. This decoupling allows the spin precession to be computed with negligible additional cost.
   • The method achieves exponential (spectral) convergence, requiring modest grid sizes (e.g., $N=53$ to $101$ points) for high precision.

3. Spacetime Models:

   • Covariant Effective Metric: Retains the Schwarzschild $g_{tt}$ but features $g_{tt} \neq -1/g_{rr}$, controlled by a quantum parameter $\alpha = \zeta/M$.
   • Self-Dual Loop Quantum Black Hole: Obtained by polymerizing the Schwarzschild interior, featuring a minimum area gap and a polymeric parameter $P$. Here, the areal radius differs from the radial coordinate.

Key Results
The study computes the cosine of the spin precession angle, $\cos \xi$, for massive particles in backward scattering channels ($\Delta\phi = \pi, 3\pi, 5\pi, 7\pi$) across varying quantum parameters and particle velocities ($v$).

• Quantum Parameter Dependence:

  • In the non-relativistic limit ($v \ll 1$), quantum corrections are most pronounced.
  • For the covariant metric, deviations from Schwarzschild values are moderate and grow with the channel number.
  • For the self-dual black hole, the polymeric parameter $P$ acts as the dominant deformation, shifting the strong-field region outward and causing substantial drifts in $\cos \xi$. The minimum-area parameter $\hat{a}$ has a negligible effect.
  • In the $\pi$ channel (near-parabolic orbits), the $\pi$ channel remains close to the Schwarzschild value ($+1$) for the covariant metric and the self-dual metric (with respect to $\hat{a}$), but the self-dual metric shows sensitivity to $P$.

• Velocity Dependence:

  • As $v \to 1$ (ultra-relativistic limit), $\cos \xi$ universally approaches $-1$ for all channels and spacetimes, recovering the known result that massless particles undergo a complete spin flip in backward scattering.
  • The number of zeros in $\cos \xi$ as a function of $v$ in the interval $(0, 1)$ equals the channel number $n$.
  • The approach to the ultra-relativistic limit is not universal. The slope $d(\cos \xi)/dv$ at $v=1$ is a strong-field quantity that depends visibly on the quantum parameters of both spacetimes.

• Slope at $v=1$:

  • While $\cos \xi(1) = -1$ is universal, the rate of approach varies. For the covariant metric, the magnitude of the slope increases with $\alpha$. For the self-dual black hole, the polymeric parameter $P$ significantly reduces the slope magnitude (by up to a factor of three), whereas $\hat{a}$ has minimal impact.
  • This slope serves as a "clean imprint" of the underlying quantum geometry, even in channels (like $\pi$) where the low-velocity limit is nearly indistinguishable from Schwarzschild.

Significance and Claims
The paper claims three primary contributions:

1. Methodological: It demonstrates that pseudo-spectral methods are highly effective for gravitational lensing, transforming the problem into a well-defined BVP where the impact parameter is a nonlinear eigenvalue. This approach is computationally efficient and accurate, suitable for strong-field regimes where traditional integration might struggle or require fine-tuning.
2. Physical: It establishes that spin precession in backward scattering is a sensitive probe of quantum geometry for massive particles. While the ultra-relativistic limit washes out quantum corrections (universal spin flip), the non-relativistic regime and the specific rate of approach to the ultra-relativistic limit retain distinct signatures of the LQG parameters.
3. Generality: The framework is not restricted to the specific LQG models studied. It applies to any static, spherically symmetric spacetime defined by three functions $f, g, h$, allowing for the handling of diverse metrics (including those with mismatched $g_{tt}$ and $g_{rr}$ or non-standard areal radii) on the same mathematical footing.

The authors note that while the quantum parameters are expected to be tiny for astrophysical black holes, the scans performed amplify these signals to reveal the structural dependence of spin precession on the underlying geometry. The work lays the groundwork for future extensions to rotating black holes and bound orbits.