pith. sign in

arxiv: 2604.15975 · v1 · submitted 2026-04-17 · 🌀 gr-qc

The double Schwarzschild solution in bispherical coordinates

Christian Klein , El Mehdi Zejly This is my paper

Pith reviewed 2026-05-10 07:58 UTC · model grok-4.3

classification 🌀 gr-qc
keywords bispherical coordinatesdouble Schwarzschild solutionconformal transformationelliptic functionsspectral methodsWeyl coordinatesblack hole spacetimes
0
0 comments X

The pith

An explicit elliptic-function map transforms the equal-mass double Schwarzschild solution into bispherical coordinates.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper constructs a conformal transformation that takes the standard cylindrical Weyl coordinates of the equal-mass double Schwarzschild solution to bispherical coordinates, with the map written in closed form using elliptic functions. It then introduces a multi-domain spectral method adapted to bispherical coordinates and applies the method to recover the solution numerically. The coordinate system matches the natural geometry of two distinct black holes separated along an axis. This supplies both an analytic expression for the coordinate change and a practical numerical scheme for the spacetime.

Core claim

The equal-mass double Schwarzschild solution is studied in bispherical coordinates. An explicit conformal transformation from cylindrical Weyl coordinates to bispherical coordinates is given in terms of elliptic functions. A multi-domain spectral method for spacetimes in bispherical coordinates is presented to numerically reconstruct this solution.

What carries the argument

The conformal transformation from Weyl to bispherical coordinates expressed via elliptic functions, which carries the vacuum metric while adapting to the two-center geometry.

If this is right

  • The vacuum Einstein equations remain satisfied after the coordinate change.
  • The multi-domain spectral method produces a faithful numerical representation of the solution in the new coordinates.
  • The bispherical system avoids artificial coordinate singularities that appear in Weyl coordinates for two separated sources.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The same elliptic-function technique may apply to other axisymmetric vacuum solutions written in Weyl form.
  • Numerical evolution schemes could be built on this coordinate system to study time-dependent binary black hole dynamics.
  • Comparison of the reconstructed metric against known asymptotic properties of the double Schwarzschild solution offers a direct test of the method's accuracy.

Load-bearing premise

The elliptic-function map is truly conformal and reproduces the original double Schwarzschild geometry exactly, without adding singularities or breaking the vacuum equations.

What would settle it

A pointwise check showing that the pulled-back metric fails to satisfy the vacuum Einstein equations at a regular point in the domain.

Figures

Figures reproduced from arXiv: 2604.15975 by Christian Klein, El Mehdi Zejly.

Figure 1
Figure 1. Figure 1: The metric functions (7) for the double Schwarzschild solution for m1 = m2 = 1 and R0 = 4, on the left f, on the right e 2k . There are two horizons located on the symmetry axis between R0/2±m1 and −R0/2±m2. Both functions, f and e 2k , vanish there. The solution is asymptotically flat. On the regular part of the axis, here for |z| > 3, the metric function k vanishes. However, this is not the case between … view at source ↗
Figure 2
Figure 2. Figure 2: Illustration of the angle θ at point P in the x-z plane. Here, d(P, F1) and d(P, F2) denote the Euclidean distances from P to F1 and F2, respectively. If ψ ∈ [0, 2π) denotes the azimuthal angle around the z-axis, then (η, θ, ψ) are called the bispherical coordinates of P. These coordinates are described in more detail in [36]. The surfaces {η = const.} are nested spheres with center (0, 0, a coth η) and ra… view at source ↗
Figure 3
Figure 3. Figure 3: The surfaces {η = const.} for η = 0.25, 0.5, 1, 2, and 3 (shown in red), together with the corresponding symmetric surfaces for negative values of η (shown in blue), with a = 1 [PITH_FULL_IMAGE:figures/full_fig_p006_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: The sections of the surfaces {η = const.} in the plane spanned by the x- and z-axes, for the same values as in [PITH_FULL_IMAGE:figures/full_fig_p006_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: Weyl coordinates as functions of the bispherical coordinates: on the left, ρ, and on the right, z. The pole at u = 0 is clearly visible. In order to better observe how the Weyl coordinates vary with respect to the bispherical coor￾dinates, it is convenient to normalize them by the factor Q defined in (9), since it vanishes at the pole with the same order, as shown in Appendix C [PITH_FULL_IMAGE:figures/fu… view at source ↗
Figure 6
Figure 6. Figure 6: Normalized Weyl coordinates as functions of the bispherical coordinates: on the left, Qρ, and on the right, Qz. Denoting by (13) habdx a dx b := e 2k (dρ 2 + dz 2 ) + ρ 2 dϕ 2 the spatial part of the metric, up to the conformal factor f −1 , in Weyl coordinates, one finds that the metric in bispherical coordinates can be written as hηη = hθθ = |w ′ (η + iθ)| 2 e 2k , hϕϕ = ρ 2 (η, θ) [PITH_FULL_IMAGE:figu… view at source ↗
Figure 7
Figure 7. Figure 7: Weyl potentials as functions of the bispherical coordinates: on the left, f, and on the right, e 2k . Note that the boundary of the domains in [PITH_FULL_IMAGE:figures/full_fig_p009_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: The functions W (left) and U (right) as defined in (17) and (20) respec￾tively [PITH_FULL_IMAGE:figures/full_fig_p009_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: The function W (17) for the double Schwarzschild solution for m1 = m2 = 1 and R0 = 5 on the left and the spectral coefficients on the right. The corresponding solution for ρQ is shown on the left of [PITH_FULL_IMAGE:figures/full_fig_p011_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: The function ρQ (17) for the double Schwarzschild solution for m1 = m2 = 1 and R0 = 5 on the left and the difference between numerical and exact solution on the right. 6. Multi-domain spectral approach In this section we will summarise basic concepts of the multi-domain spectral approach which is also the basis of Kadath [34]. 6.1. Domains. Spectral methods as discussed in the previous section are very ef… view at source ↗
Figure 11
Figure 11. Figure 11: The five numerical domains. implemented with a τ -method. At the rectangle cut out near infinity, the exact solution will be imposed. 6.2. The Ernst potential f. In order to take care of the horizons where f vanishes, we make the ansatz (20) f = (1 − η 2 /η2 0 ) 2 e U and get for the Ernst equation (21) Uηη + (ln ρ)ηUη − 4 η 2 0 − η 2 − 8η 2 (η 2 0 − η 2 ) 2 − 4η η 2 0 − η 2 (ln ρ)η + Uθθ + (ln ρ)θUθ = 0.… view at source ↗
Figure 12
Figure 12. Figure 12: The spectral coefficients for U in the domains I (left), II (middle) and IV (right), for m1 = m2 = 1 and R0 = 5. 6.3. Numerical solution for U. The numerical solution with the above numerical parameters can be seen in [PITH_FULL_IMAGE:figures/full_fig_p014_12.png] view at source ↗
Figure 13
Figure 13. Figure 13: The function U (20) for the double Schwarzschild solution for m1 = m2 = 1 and R0 = 5 on the left, and the difference between numerical and exact solution on the right. The solution for the Ernst equation itself is shown in [PITH_FULL_IMAGE:figures/full_fig_p014_13.png] view at source ↗
Figure 14
Figure 14. Figure 14: The function f (20) for the double Schwarzschild solution for m1 = m2 = 1 and R0 = 5 on the left, and the difference between numerical and exact solution on the right. with a vanishing condition on the boundary of the computational domain and the Euler-Darboux equation for the logarithm of the (real) Ernst potential. If the behavior of these two functions at the horizons and the axis is addressed by an ex… view at source ↗
Figure 15
Figure 15. Figure 15: Mapping from the (half) fundamental domain D to the Weyl coordinate domain Dw through the bispherical domain Db [PITH_FULL_IMAGE:figures/full_fig_p017_15.png] view at source ↗
Figure 16
Figure 16. Figure 16: Domain-coloring plot of the conformal map w in the domain Db. The θ axis is reversed in this figure. Appendix B. Uniqueness of the conformal map The map w is in fact the unique conformal map that preserves the structure of the horizons in Weyl coordinates. Indeed, suppose that there exists another conformal map w˜ : Db → Dw such that w˜(η0) = i R0 2 + m  , w˜(η0 + iπ) = i R0 2 − m  , w˜(−η0) = −i [PITH… view at source ↗
read the original abstract

The double Schwarzschild solution in the equal mass case is studied in bispherical coordinates. An explicit conformal transformation from cylindrical Weyl coordinates to bispherical coordinates is given in terms of elliptic functions. A multi-domain spectral method for spacetimes in bispherical coordinates is presented to numerically reconstruct this solution.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The manuscript studies the equal-mass double Schwarzschild solution in bispherical coordinates. It claims to provide an explicit conformal transformation from cylindrical Weyl coordinates to bispherical coordinates expressed in terms of elliptic functions, and presents a multi-domain spectral method to numerically reconstruct the solution in the new coordinates.

Significance. If the elliptic-function map is confirmed to be conformal, to satisfy the vacuum Einstein equations exactly, and to map the rod singularities without introducing extraneous curvature or coordinate singularities, the work would supply a concrete coordinate representation of a known binary black-hole geometry that could facilitate analytic or numerical studies of equal-mass systems. The multi-domain spectral approach, if shown to converge, would also constitute a reusable numerical tool for other vacuum spacetimes in bispherical coordinates.

major comments (2)
  1. [the section presenting the coordinate transformation] The central claim that the elliptic-function transformation is conformal and exactly preserves the double-Schwarzschild geometry rests on an unverified assertion. The manuscript must demonstrate explicitly that the given map satisfies the two-dimensional Cauchy-Riemann conditions (or the equivalent Jacobian relation for conformality) throughout the domain excluding the axis and horizons, and that the branch structure maps the two rod singularities precisely onto the bispherical coordinate singularities without extra zeros, poles, or curvature.
  2. [the section describing the multi-domain spectral method] The numerical reconstruction via the multi-domain spectral method lacks reported error estimates, convergence tests with respect to the number of domains or spectral order, and direct comparison against the known Weyl-coordinate values of the metric functions. Without these, it is impossible to assess whether the numerical solution faithfully reproduces the analytic double-Schwarzschild geometry.
minor comments (2)
  1. [the coordinate transformation] The precise definition of the elliptic functions (including the choice of branch cuts and integration constants) should be stated with an explicit equation or reference so that the map can be reproduced independently.
  2. [introduction or coordinate section] Notation for the bispherical coordinates and the conformal factor should be introduced once and used consistently; currently the relation between the cylindrical (ρ,z) and bispherical variables is not summarized in a single equation.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading and constructive comments on our manuscript. We address each major comment below and indicate the revisions we will make to strengthen the presentation.

read point-by-point responses

  1. Referee: [the section presenting the coordinate transformation] The central claim that the elliptic-function transformation is conformal and exactly preserves the double-Schwarzschild geometry rests on an unverified assertion. The manuscript must demonstrate explicitly that the given map satisfies the two-dimensional Cauchy-Riemann conditions (or the equivalent Jacobian relation for conformality) throughout the domain excluding the axis and horizons, and that the branch structure maps the two rod singularities precisely onto the bispherical coordinate singularities without extra zeros, poles, or curvature.

    Authors: We agree that an explicit verification strengthens the central claim. The transformation was constructed from the known conformal properties of the elliptic functions that solve the relevant Laplace equation, but the manuscript did not include a direct check of the Cauchy-Riemann conditions or the Jacobian determinant. In the revised manuscript we will add a short subsection that computes the Jacobian relation throughout the domain (away from the axis and horizons), verifies that it equals the required conformal factor, and analyzes the branch cuts to confirm that the two rod singularities map precisely onto the bispherical coordinate singularities without introducing extraneous zeros, poles, or curvature singularities. revision: yes

  2. Referee: [the section describing the multi-domain spectral method] The numerical reconstruction via the multi-domain spectral method lacks reported error estimates, convergence tests with respect to the number of domains or spectral order, and direct comparison against the known Weyl-coordinate values of the metric functions. Without these, it is impossible to assess whether the numerical solution faithfully reproduces the analytic double-Schwarzschild geometry.

    Authors: We accept that quantitative validation is required. The original manuscript described the multi-domain spectral method and showed representative solutions but did not report error norms, convergence studies, or direct comparisons with the Weyl-coordinate data. In the revised version we will add these diagnostics: L2 and maximum-norm errors for increasing spectral order and number of domains, together with pointwise comparisons of the metric functions against the known Weyl expressions at selected interior points. revision: yes

Circularity Check

0 steps flagged

No significant circularity detected in the derivation chain.

full rationale

The paper starts from the established Weyl form of the equal-mass double Schwarzschild solution in cylindrical coordinates and supplies an explicit conformal map to bispherical coordinates expressed via elliptic functions, followed by a multi-domain spectral reconstruction. This constitutes a direct coordinate transformation and numerical verification rather than any self-definitional loop, fitted parameter renamed as prediction, or load-bearing self-citation chain. The central claim remains independent of its own outputs, with the elliptic map asserted as a construction from known inputs and the spectral method serving as external reconstruction.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The central claims rest on the existence of the double Schwarzschild solution in Weyl coordinates and on standard properties of conformal transformations in vacuum GR; no free parameters or new entities are introduced in the abstract.

axioms (2)
  • domain assumption The equal-mass double Schwarzschild solution satisfies the vacuum Einstein equations in Weyl coordinates.
    Standard background fact in general relativity invoked to justify the starting point of the transformation.
  • standard math A conformal rescaling preserves the vacuum character of the metric.
    Basic property of conformal transformations in four-dimensional Lorentzian geometry used to map the solution.

pith-pipeline@v0.9.0 · 5339 in / 1246 out tokens · 36444 ms · 2026-05-10T07:58:34.218332+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

39 extracted references · 39 canonical work pages

  1. [1]

    Observation of gravitational waves from a binary black hole merger,

    B. P. A. et al., “Observation of gravitational waves from a binary black hole merger,” Physical Review Letters , vol. 116, p. 061102, 2016

    work page 2016

  2. [2]

    Detweiler in Frontiers in Numerical Relativity (C

    S. Detweiler in Frontiers in Numerical Relativity (C. R. Evans, L. S. Finn, and D. W. Hobill, eds.), p. 43, Cambridge, England: Cambridge University Press, 1989

    work page 1989

  3. [3]

    Blackburn and S

    J. Blackburn and S. Detweiler Physical Review D , vol. 46, p. 2318, 1992

    work page 1992

  4. [4]

    Detweiler Physical Review D , vol

    S. Detweiler Physical Review D , vol. 50, p. 4929, 1994

    work page 1994

  5. [5]

    Schönberg Physical Review, vol

    M. Schönberg Physical Review, vol. 69, p. 211, 1946

    work page 1946

  6. [6]

    Schild Physical Review, vol

    A. Schild Physical Review, vol. 131, p. 2762, 1963

    work page 1963

  7. [7]

    Friedman, K

    J. Friedman, K. Uryu, and M. Shibata Physical Review D , vol. 65, p. 064035, 2002

    work page 2002

  8. [8]

    Gourgoulhon, P

    E. Gourgoulhon, P. Grandclément, and S. Bonazzola Physical Review D , vol. 65, p. 044020, 2002

    work page 2002

  9. [9]

    Grandclément, E

    P. Grandclément, E. Gourgoulhon, and S. Bonazzola Physical Review D , vol. 65, p. 044021, 2002

    work page 2002

  10. [10]

    Andrade, C

    Z. Andrade, C. Beetle, A. Blinov, B. Bromley, L. Burko, M. Cranor, R. Owen, and R. Price Physical Review D, vol. 70, p. 064001, 2003

    work page 2003

  11. [11]

    Caudill, G

    M. Caudill, G. B. Cook, J. D. Grigsby, and H. P. Pfeiffer Physical Review D , vol. 74, p. 064011, 2006

    work page 2006

  12. [12]

    K. Uryu, F. Limousin, J. L. Friedman, E. Gourgoulhon, and M. Shibata Physical Review D , vol. 80, p. 124004, 2009

    work page 2009

  13. [13]

    J. L. Friedman and K. Uryu Physical Review D , vol. 73, p. 104039, 2006

    work page 2006

  14. [14]

    N. T. Bishop, R. Gómez, L. Lehner, M. Maharaj, and J. Winicour Physical Review D , vol. 72, p. 024002, 2005

    work page 2005

  15. [15]

    K. Uryu, E. Gourgoulhon, and C. Markakis Physical Review D , vol. 82, p. 104054, 2010

    work page 2010

  16. [16]

    Yoshida, B

    S. Yoshida, B. C. Bromley, J. S. Read, K. Uryu, and J. L. Friedman Classical and Quantum Gravity , vol. 23, no. 16, 2006

    work page 2006

  17. [17]

    Beetle, B

    C. Beetle, B. Bromley, and R. H. Price Physical Review D , vol. 74, p. 024013, 2006

    work page 2006

  18. [18]

    Bromley, R

    B. Bromley, R. Owen, and R. H. Price Physical Review D , vol. 71, p. 104017, 2005

    work page 2005

  19. [19]

    Beetle, B

    C. Beetle, B. Bromley, N. Hernández, and R. H. Price Physical Review D , vol. 76, p. 084016, 2007

    work page 2007

  20. [20]

    S. R. Lau and R. H. Price Journal of Computational Physics , vol. 227, no. 2, p. 1126, 2007

    work page 2007

  21. [21]

    S. R. Lau and R. H. Price Journal of Computational Physics , vol. 231, no. 22, p. 7695, 2012

    work page 2012

  22. [22]

    Bonazzola, J

    S. Bonazzola, J. L. Jaramillo, and J. Novak Classical and Quantum Gravity , vol. 24, no. 16, 2007

    work page 2007

  23. [23]

    Klein Physical Review D , vol

    C. Klein Physical Review D , vol. 70, p. 124026, 2004

    work page 2004

  24. [24]

    Ehlers, Konstruktion und Charakterisierungen von Lösungen der Einstein ’schen Gravitationsgleichungen

    J. Ehlers, Konstruktion und Charakterisierungen von Lösungen der Einstein ’schen Gravitationsgleichungen . PhD thesis, University of Hamburg, 1957. Dissertation (in German)

    work page 1957

  25. [25]

    Geroch Journal of Mathematical Physics , vol

    R. Geroch Journal of Mathematical Physics , vol. 12, p. 918, 1971

    work page 1971

  26. [26]

    F. J. Ernst Physical Review, vol. 167, p. 1175, 1968

    work page 1968

  27. [27]

    Klein and O

    C. Klein and O. Richter, Ernst Equation and Riemann Surfaces , vol. 685 of Lecture Notes in Physics . Berlin: Springer, 2005

    work page 2005

  28. [28]

    Absence of asymptotically flat solutions of einstein’s equations which are periodic and empty near infinity,

    G. Gibbons and J. Stewart, “Absence of asymptotically flat solutions of einstein’s equations which are periodic and empty near infinity,” in Classical General Relativity (W. Bonnor, J. Islam, and M. MacCallum, eds.), p. 77, Cambridge, England: Cambridge University Press, 1984

    work page 1984

  29. [29]

    Ashtekar and B

    A. Ashtekar and B. Xanthopoulos Journal of Mathematical Physics , vol. 19, p. 2216, 1978

    work page 1978

  30. [30]

    S. Bai, G. Izquierdo, and C. Klein Physical Review D , vol. 93, p. 124001, 2016

    work page 2016

  31. [31]

    Kramer and G

    D. Kramer and G. Neugebauer Physics Letters A , vol. 75, p. 259, 1980

    work page 1980

  32. [32]

    Kramer, H

    D. Kramer, H. Stephani, E. Herlt, and M. MacCallum, Exact Solutions of Einstein ’s Field Equations . Cam- bridge: Cambridge University Press, 1980. THE DOUBLE SCHW ARZSCHILD SOLUTION IN BISPHERICAL COORDINATES 23

    work page 1980

  33. [33]

    de Leon, J

    E. de Leon, J. Frauendiener, and C. Klein, 2025

    work page 2025

  34. [34]

    Grandclément Journal of Computational Physics , vol

    P. Grandclément Journal of Computational Physics , vol. 220, p. 3334, 2010

    work page 2010

  35. [35]

    V. S. Manko and E. Ruiz Physics Letters B , vol. 794, pp. 36–40, 2019

    work page 2019

  36. [36]

    Parry Moon, Domina Eberle Spencer, Field theory handbook , pp. 110–112. Springer Verlag, 1961

    work page 1961

  37. [37]

    Lawden, Elliptic functions and applications , vol

    D. Lawden, Elliptic functions and applications , vol. 80 of Applied Mathematical Sciences. New York: Springer, 1989

    work page 1989

  38. [38]

    L. N. Trefethen, Spectral Methods in MATLAB . Philadelphia, PA: SIAM, 2000

    work page 2000

  39. [39]

    Crespo, M

    S. Crespo, M. Fasondini, C. Klein, N. Stoilov, and C. Vallée Numerical Algorithms, vol. 84(1), pp. 1–35, 2019. (C. Klein) Université Bourgogne Europe, CNRS, IMB UMR 5584, 21000 Dijon, France, Institut Universitaire de France Email address : christian.klein@u-bourgogne.fr (M. Zejly) Université Bourgogne Europe, CNRS, IMB UMR 5584, 21000 Dijon, France Email...

    work page 2019