Wave Optics Effects from Gravitational Wave Propagation Through Dark Matter Halos

Technical Summary: Wave Optics Effects from Gravitational Wave Propagation Through Dark Matter Halos

Problem Statement
Gravitational wave (GW) propagation is conventionally analyzed using the geometric optics approximation, which assumes the GW wavelength is negligible compared to the characteristic scale of the gravitational lens. However, this approximation breaks down when the Schwarzschild radius of the lens is comparable to or smaller than the GW wavelength. In such regimes, wave optics effects—specifically diffraction and interference—become significant, modifying GW signals in a frequency-dependent manner. While analytical methods exist, they are often restricted to idealized geometries (e.g., point lenses) and simple perturbative expansions, failing to capture the complex, inhomogeneous matter distributions found in realistic astrophysical environments like dark matter halos. There is a need for numerical frameworks that can evolve GWs through realistic potentials to understand how wave optics effects influence signal interpretation and lens property inference.

Methodology
The authors present a numerical framework for evolving linearized plane GWs through stationary, spherical astrophysical structures in both weak and strong gravity regimes. The core of the methodology involves:

• Full Einstein Equations: The authors evolve the full set of 10 components of the linearized Einstein equations for trace-reversed metric perturbations ($\bar{h}_{\mu\nu}$) in a curved background, rather than relying solely on simplified scalar wave equations. The background metric is defined by a gravitational potential $\phi$ derived from specific density profiles.
• Potential Models: Three distinct density profiles are simulated:
  1. Gaussian: Mimics a standard spherical mass distribution.
  2. Navarro–Frenk–White (NFW): A common model for cold dark matter halos.
  3. Burkert: A model often used for dwarf galaxies.
     For NFW and Burkert profiles, where densities diverge at the origin, the authors solve a "density-free" system ($T_{\mu\nu} = 0$) to isolate curvature effects, while the Gaussian profile is simulated with and without the source term to quantify the impact of the matter coupling.
• Numerical Implementation: The simulations utilize the CarpetX framework within the Einstein Toolkit, employing Adaptive Mesh Refinement (AMR) to resolve steep curvature gradients and the propagating wave packet simultaneously. The equations are discretized using a fourth-order central finite difference stencil and evolved via a fourth-order Runge-Kutta integrator.
• Initial Conditions: Circularly polarized monochromatic waves are initialized in a flat Minkowski region and propagated through the curved spacetime. To handle boundary conditions and avoid spurious gauge modes, a quintic spline is used to smoothly transition the potential to zero at the simulation boundaries.
• Regimes: Simulations are conducted in both weak gravity ($|\phi_{max}| \approx 0.1$) and strong gravity ($|\phi_{max}| \approx 0.49$) settings. The GW wavelength is set to be less than the Schwarzschild radius of the structure to suppress pure diffraction effects and focus on interference and curvature-induced scattering.

Key Contributions and Results

• Deviation from Null Geodesics: The simulations demonstrate that GWs do not strictly propagate along null geodesics when passing through a Gaussian density distribution, even in weak gravity. In strong gravity, this deviation is more pronounced, with parts of the wavefront scattering in directions not orthogonal to the geodesics.
• Role of Potential Shape: The scattering and amplification of GWs depend not only on the total mass of the lens but strongly on the specific shape of the gravitational potential. The Burkert potential in the strong gravity regime produces the most significant amplification (nearly a factor of 30) and exhibits unique behavior, such as the flipping of the plane wave's convexity upon exiting the structure.
• Comparison with Simplified Models:
  • Inside the Lens: When the wave is inside the central potential, differences between the full Einstein equation evolution and linearized scalar wave predictions can be of order unity (significant).
  • Outside the Lens: Once the wave exits the structure, the discrepancy between the full solution and simplified models reduces to a few percent for most cases.
  • Density Coupling: In the strong gravity Gaussian case, the presence of the density source ($T_{\mu\nu} \neq 0$) induces a weak but non-negligible scattering of the GW into multiple directions (amplitude $\sim$ few percent), a feature absent in the null geodesic description.
• Limitations of Simplified Equations: The authors find that simplified evolution equations fail to capture multidirectional scattering patterns and large transient amplifications in strong curvature regimes. Notably, for the Burkert potential in strong gravity, the simplified model differs from the full Einstein result by more than 35% even after the wave exits the high-curvature region.

Significance and Claims
The paper claims that while geometric optics and simplified analytical treatments are sufficient for weak gravity scenarios, they are inadequate for accurately modeling GW propagation in strong curvature regimes or for lenses with complex potential shapes. The authors assert that:

1. Precision Modeling: Future GW detectors and Pulsar Timing Arrays (PTAs) may be sensitive to these wave optics signatures. Accurate signal modeling for inferring lens properties or cosmological parameters requires moving beyond simplified analytical treatments to full numerical relativity simulations.
2. Dark Matter Constraints: These wave optics effects could potentially help constrain the structure of dark matter spikes or halos, as the signal modifications carry information about the mass distribution and potential shape.
3. Machine Learning Training: The authors highlight a risk in using simple analytical models to train machine learning algorithms for GW data analysis; such models may produce biased results compared to training on realistic simulations that include full wave optics effects.
4. Methodological Validation: The study provides a benchmark for testing simplified analytical methods against more complex systems, demonstrating that while simplified models capture the main features in weak fields, they miss critical transient and scattering effects in strong fields.

The work concludes that full Einstein equation simulations are essential for interpreting GW signals in regimes where finite-wavelength effects and detailed lens structures produce non-trivial waveform modifications not captured by geometric optics alone.