The double Schwarzschild solution in bispherical coordinates

This paper studies the equal-mass double Schwarzschild solution in bispherical coordinates by providing an explicit conformal transformation from Weyl coordinates using elliptic functions and presenting a multi-domain spectral method to numerically reconstruct the solution.

The Gist

The Big Picture: Two Black Holes Holding Hands (Sort Of)

Imagine two massive black holes floating in space, staring at each other. In reality, they would spiral inward and crash together, creating a massive explosion of gravitational waves (like the ripples in a pond when you drop a stone).

However, this paper isn't about the crash. It's about a very specific, frozen moment in time: What if they were perfectly balanced and not moving?

To keep them from crashing, imagine there is an invisible, unbreakable rod (called a "Weyl strut") pushing them apart. This isn't a real physical rod; it's a mathematical trick to keep the system static so scientists can study it. The authors of this paper wanted to map out the shape of space around these two frozen black holes with extreme precision.

The Problem: The Wrong Map

To draw a map of this space, scientists usually use a standard grid system called Weyl coordinates. Think of this like a standard city map with straight streets running North-South and East-West.

The problem is that black holes are round spheres. Trying to map a round sphere onto a grid of straight lines is like trying to wrap a flat piece of graph paper perfectly around a basketball. You get wrinkles, tears, and it's very hard to calculate the exact shape of the ball. In math terms, the equations get "singular" (they break or become infinite) right where the black holes are.

The Solution: A New Shape-Shifting Lens

The authors, Christian Klein and El Mehdi Zejly, decided to change the map entirely. Instead of a flat grid, they used a coordinate system called Bispherical Coordinates.

The Analogy:
Imagine you have a pair of glasses that can reshape reality.

• Old Glasses (Weyl): You see the black holes as long, flat strips on a line. It's messy to calculate the curves.
• New Glasses (Bispherical): You put on the new glasses, and suddenly, the space around the black holes transforms. The black holes now look like perfect, smooth bubbles sitting on a coordinate grid. The "rod" holding them apart is still there, but the grid lines wrap around the bubbles perfectly.

This new system is called "bispherical" because it naturally fits two spheres (the two black holes) side-by-side.

The Magic Trick: Elliptic Functions

How did they switch from the old map to the new one? They used a mathematical "magic wand" involving Jacobi Elliptic Functions.

Think of these functions as a complex, twisting rubber sheet. The authors found a specific way to stretch and twist this sheet so that the messy, broken points of the old map (where the black holes are) became smooth, clean points on the new map. They wrote down the exact formula for this transformation, which is a major achievement because it had never been done this clearly before.

The Computer Challenge: The "Pixel" Problem

Once they had the perfect map, they wanted to use a computer to solve Einstein's equations (the rules of gravity) to see if the computer could recreate this known solution.

The Challenge:
Computers solve problems by breaking them into tiny pieces (like pixels on a screen).

• The Smooth Parts: Most of space is smooth. A computer can solve this easily, like filling in a blue sky.
• The "Rod" Problem: The invisible rod holding the black holes apart creates a sharp "kink" or a "cusp" in the math. It's like trying to draw a sharp corner with a soft brush. If you try to draw the whole picture with one big brush, the corner looks blurry and jagged (this is called the Gibbs phenomenon).

The Fix: Multi-Domain Spectral Method
To fix this, the authors didn't use one big brush. They used a patchwork quilt approach (called a multi-domain method).

1. They cut the map into five different zones.
2. In the smooth zones, they used high-resolution math to get a perfect picture.
3. In the tricky zones (near the black holes and the "rod"), they used smaller, specialized patches to handle the sharp corners without blurring.

The Result: Machine Precision

When they ran their computer simulation, the result was stunning.

• They took the "exact" mathematical solution they knew was correct.
• They let their computer try to guess it using their new method.
• The computer's guess matched the exact solution almost perfectly—down to 15 decimal places (machine precision).

It's as if you asked a computer to draw a perfect circle, and it drew one so perfect that if you zoomed in with a microscope, you couldn't tell the difference between the drawing and a real circle.

Why Does This Matter?

You might ask, "Why study two black holes that aren't moving?"

1. Testing the Tools: Before scientists can simulate black holes crashing (which is what we actually detect with gravitational wave detectors), they need to make sure their computer codes work perfectly. This paper proves their new "Bispherical Glasses" and "Patchwork Quilt" method works flawlessly.
2. Future Simulations: The ultimate goal is to simulate black holes that are moving and spinning, held in a "helical" balance (like a dance). This paper is the training ground. If they can perfectly simulate the static version, they are one step closer to simulating the dynamic, crashing version that creates the gravitational waves we hear in the universe.

Summary

The authors invented a new way to look at two black holes (Bispherical Coordinates) that makes the math much easier. They proved that their computer method can solve the equations for this setup with incredible accuracy, paving the way for better simulations of black hole collisions in the future.

Technical Summary

1. Problem Statement

The paper addresses the numerical and analytical challenges associated with modeling binary black hole systems, specifically focusing on the equal-mass static double Schwarzschild solution.

• Context: Binary black holes are primary sources of gravitational waves. While dynamic evolution requires full numerical relativity, early-stage evolution is often approximated by spacetimes with a helical Killing vector (quasi-stationary phase).
• The Challenge: Solving the Einstein equations for such systems is difficult due to:
  • Singularities: The presence of Killing horizons and a "light cylinder" where the Killing vector changes causal nature.
  • Coordinate Issues: The standard Weyl (cylindrical) coordinates represent black hole horizons as intervals on the symmetry axis. This is numerically inconvenient because the horizons are not constant coordinate surfaces, and the axis contains a conical singularity (the Weyl strut) required to keep the black holes static.
  • Infinity: Standard coordinates do not naturally compactify spatial infinity, making boundary conditions at infinity difficult to handle numerically.
• Goal: To construct an explicit representation of the double Schwarzschild solution in bispherical coordinates, where horizons are constant coordinate surfaces, and to demonstrate a high-precision numerical reconstruction of this solution as a testbed for future helical Killing vector simulations.

2. Methodology

A. Analytical Transformation (Weyl to Bispherical)

The authors derive an exact conformal transformation from Weyl coordinates $(\rho, z)$ to bispherical coordinates $(\eta, \theta)$.

• Bispherical Coordinates: Defined by two focal points $F_1, F_2$ on the $z$-axis. Surfaces of constant $\eta$ are spheres, and surfaces of constant $\theta$ are toroidal surfaces.
• The Conformal Map: The transformation is constructed using Jacobi elliptic functions. The complex coordinate $w = \rho + iz$ is mapped to $u = \eta + i\theta$ via:
  $$w(u) = i \left(\frac{R_0}{2} + m\right) \text{ns}\left(\frac{K}{\eta_0}u, \mu\right)$$
  where $\text{ns}(x) = 1/\text{sn}(x)$, $K$ is the complete elliptic integral of the first kind, and parameters $\mu$ and $\eta_0$ depend on the black hole mass $m$ and separation $R_0$.
• Key Features of the Map:
  • The horizons (intervals in Weyl) become the spherical surfaces $\eta = \pm \eta_0$.
  • Spatial infinity (the origin in Weyl $\rho=0, z=0$ is not infinity; infinity is $\rho^2+z^2 \to \infty$) corresponds to the pole $u=0$ in the bispherical domain.
  • The symmetry axis is split: $\theta=0$ represents the regular axis outside the horizons, and $\theta=\pi$ represents the axis between the horizons (containing the Weyl strut).

B. Numerical Approach

The authors employ a multi-domain spectral method (similar to the Kadath code) to solve the Einstein equations in these new coordinates.

• Equations: The vacuum Einstein equations reduce to:
  1. A harmonic equation for the metric function related to $\rho$ (specifically a function $W$).
  2. The Ernst equation (a non-linear elliptic PDE) for the metric potential $f$.
• Handling Singularities:
  • Horizons and Axis: The authors use explicit ansatzes to factor out the known singular behavior of the metric functions at the horizons ($\eta = \pm \eta_0$) and the axis ($\theta = 0, \pi$). This reduces the problem to solving for regular functions.
  • Infinity: To handle the cusp-like behavior of the Ernst potential $f$ at infinity (where $f \sim 1 - 2M/r$), the computational domain is modified. A small rectangular region near the pole $u=0$ (infinity) is excised, and exact boundary conditions are imposed on the cut boundary.
• Discretization:
  • Chebyshev Collocation: Functions are expanded in Chebyshev polynomials.
  • Multi-Domain Decomposition: The domain is split into five sub-domains (I–V) to handle the non-smoothness of the metric function $e^{2k}$ (caused by the Weyl strut) and the cusp at infinity. This avoids the Gibbs phenomenon and improves conditioning.
  • Boundary Conditions: Implemented via the Lanczos $\tau$-method, replacing differential equation rows in the differentiation matrices with boundary constraints.

3. Key Contributions

1. Explicit Conformal Map: The paper provides, for the first time, an explicit formula for the double Schwarzschild solution in bispherical coordinates using Jacobi elliptic functions.
2. Coordinate System Optimization: It demonstrates that bispherical coordinates naturally compactify infinity to a single point and align the horizons with constant coordinate surfaces, significantly simplifying the numerical boundary conditions.
3. Numerical Validation: The authors successfully reconstruct the exact solution numerically.
   • They solved the linearized harmonic equation for the spatial metric component.
   • They solved the non-linear Ernst equation for the potential $f$.
4. Uniqueness Proof: Appendix B provides a mathematical proof that the constructed conformal map is the unique map preserving the horizon structure.

4. Results

• Accuracy: The numerical reconstruction achieves machine precision (order $10^{-12}$ to $10^{-16}$) for the regular parts of the metric.
  • The function $W$ (related to $\rho$) was reproduced with errors $\sim 10^{-13}$.
  • The Ernst potential $U$ (related to $f$) showed errors $\sim 10^{-9}$ to $10^{-12}$, with the largest errors occurring near the horizon where the potential vanishes.
• Spectral Convergence: The spectral coefficients of the solution decay exponentially, confirming the analyticity of the solution within the sub-domains and the efficacy of the spectral method.
• Metric Behavior: The study confirmed the expected behavior of the metric near the horizons (vanishing $f$) and the axis (conical singularity/Weyl strut where $e^{2k} < 1$).
• Limitation: The metric function $e^{2k}$, which is discontinuous due to the Weyl strut, was not numerically resolved to high precision (as spectral methods struggle with discontinuities), but the authors argue this is expected and not the primary focus for future helical Killing vector studies.

5. Significance and Outlook

• Preparatory Step: This work serves as a critical "test case" for developing numerical algorithms for binary black holes with helical Killing vectors. The double Schwarzschild solution is static (a special case of helical symmetry with zero rotation), making it an ideal benchmark.
• Future Applications: The techniques developed here (multi-domain spectral methods, handling Fuchsian singularities at horizons, and excising infinity) will be directly applied to the more complex, rotating binary black hole systems.
• Mathematical Insight: The explicit connection between Weyl coordinates and bispherical coordinates via elliptic functions provides a deeper mathematical understanding of the geometry of multi-black hole spacetimes.

In summary, the paper successfully bridges the gap between the analytical exact solution of the double Schwarzschild metric and high-precision numerical relativity techniques, establishing a robust framework for tackling the more complex problem of rotating binary black holes.