Pith. Sign in

REVIEW

Accurate Evolutions of Orbiting Black-Hole Binaries Without Excision

PaperPDF from arXiv

data files

ReviewThe objections and verdict, as JSON Paper breakdownThe authors' claims, conditions and measurements, as JSON Signed receiptProof this review has not been altered
X Bluesky LinkedIn Reddit Hacker News
M. Campanelli, C. O. Lousto, P. Marronetti, Y. Zlochower

General Relativity and Quantum Cosmology · Submitted 9 Nov 2005 · 29 references

Pith's verdict: The paper that made moving punctures the default: a real method with a real demonstration, though the puncture-regularity argument is power counting plus a gauge chosen by trial rather than a proof.

T0 · AI review

How this review was made: An AI referee (selfhosted/deepseek-v4.1-flash) read the paper on October 29, 2026 and answered a fixed list of questions about it. T0 is the starting level every paper gets. The levels describe how deeply Pith has checked the paper, never the authors.

  1. T0 AI review An AI read the paper and answered a fixed list of questions.
  2. T1 Audit Automatic checks that the references and files the paper relies on exist and match.
  3. T2 Compute The paper's code or calculations re-run with the same results.
  4. T3 Formal statement The claims written out in a formal math language a computer can read.
  5. T4 Proof Those claims proved, and the proof checked by computer.

How Pith checks papers →

· 2 major objections · 10 minor · Ask about this paper · The compiled record · Read the full review →
Challenge this review Comment (0) Read on arXiv ↗
Cite ▾
arXiv IDgr-qc/0511048 (version 2)
Pith IDDY52LRQE, Pith's permanent record of this paper and its review
CitationC. O. Lousto, M. Campanelli, P. Marronetti, Y. Zlochower. "Accurate Evolutions of Orbiting Black-Hole Binaries Without Excision." arXiv:gr-qc/0511048 (2005). https://pith.science/paper/DY52LRQE
Cite this reviewBibTeX for Pith's review of this paper
Claim paper

Why Pith says this

Highlight any text to comment on it.

Referee report

How Pith checks papers → · Full audit of this paper →

The paper

The paper claims

The paper argues that orbiting black-hole binaries can be evolved to merger without cutting the singularities out of the grid, by swapping the divergent conformal factor for a regular variable and letting the punctures move with the shift.

The reading

Numerical relativity's long-standing obstacle was the singularity: the conformal factor of puncture initial data diverges at each black hole, so codes either excised the singular region from the grid or froze the punctures in place, and both choices tended to destroy a long binary run. This paper proposes instead to evolve $¥chi = e^{-4¥varphi}$ in place of the BSSN conformal exponent $¥varphi$, since $¥chi$ is smooth where $¥varphi$ diverges, and to start from a 'pre-collapsed' lapse $¥alpha = ¥psi_{BL}^{-2}$ that vanishes like $r^2$ at each puncture. The punctures are then allowed to move, tracked by $¥partial_t ¥vec{x}_{punct} = -¥vec{¥beta}$, so no corotating shift and no excision are needed. The test case is a full evolution of an equal-mass binary started near the innermost stable circular orbit, carried through common-horizon formation to a ringing Kerr remnant, with waveforms that converge at fourth order under resolution and radiated energy that matches the measured drop in horizon mass. If the argument holds, the cheapest kind of black-hole initial data becomes directly evolvable, which removes a major piece of machinery from the binary problem.

What carries the argument

The load-bearing device is the variable substitution $¥chi = e^{-4¥varphi}$ together with a lapse that vanishes quadratically at the puncture. In the standard BSSN system the conformal exponent $¥varphi$ diverges logarithmically at a puncture, and the $¥partial_i¥varphi$ and $¥partial_i¥partial_j¥varphi$ terms in the Ricci tensor and in the $¥tilde{¥Gamma}^i$ equation inherit that divergence. Writing the equations in $¥chi$ turns those terms into $¥partial_i¥chi/¥chi$ and $¥partial_i¥partial_j¥chi/¥chi$, which are individually singular but always appear multiplied by factors of $¥chi$ or by $¥alpha ¥sim r^2$, so the products are finite and in fact smooth. The second device is to let the punctures move rather than pinning them: at $¥chi = 0$ the evolution equation for $¥chi$ reduces to $¥partial_t ¥vec{x}_{punct} = -¥vec{¥beta}$, giving a cheap tracker for the hole positions throughout the inspiral. A multiple-transition Fisheye coordinate supplies the resolution range from horizon scale to wave zone on a single grid.

Extended reading notes

Core claim

The central claim is that puncture data can be evolved stably and accurately with the singularities left on the grid. The authors replace the BSSN variable $¥varphi$ with $¥chi = e^{-4¥varphi}$, which is $C^4$ at a puncture, choose an initial lapse $¥alpha = ¥psi_{BL}^{-2}$ that is $O(r^2)$ there, and adopt the 1+log lapse $¥partial_0 ¥alpha = -2¥alpha K$ together with a Gamma-driver shift $¥partial_t ¥beta^a = B^a$, $¥partial_t B^a = ¥frac{3}{4}¥partial_t ¥tilde{¥Gamma}^a - ¥eta B^a$. With this combination the terms that blow up individually, $¥chi D_i D_j ¥alpha$, $¥alpha¥chi ¥partial_i¥partial_j¥varphi$ and $¥alpha ¥tilde{A}_{ij}¥partial_i¥varphi$, are each suppressed by enough powers of $r$ that $¥tilde{A}_{ij}$ stays $C^4$ and $¥tilde{¥Gamma}^i$ is driven from $C^0$ to $C^2$. Evolving an equal-mass binary from near ISCO, they find a common apparent horizon at $t = 19M$ with irreducible mass $0.9055M$, a polar-to-equatorial circumference ratio settling to $0.898 ¥pm 0.002$ (corresponding to spin $¥tilde{a} = 0.683 ¥pm 0.006$ and horizon mass $0.973M ¥pm 0.002M$), radiated energy $2.8¥% ¥pm 0.2¥%$ of $M$ and radiated angular momentum $J_z = -0.11M^2 ¥pm 0.01M^2$, with waveforms fourth-order convergent across $h = M/16, M/24, M/36$ and in good agreement with the Lazarus estimates.

Load-bearing premise

The claim that the evolution stays regular rests on counting powers of $r$ at the puncture to argue that every apparently divergent term is multiplied by enough of $¥chi$ or $¥alpha$ to be smooth; the paper does not prove this for the discretized equations, does not say how the code evaluates $¥partial_i¥chi/¥chi$ where $¥chi$ vanishes, and does not show a numerical convergence test of the constraint at the puncture itself.

Editorial extensions

If this is right

  • Puncture initial data, which is comparatively easy to construct, becomes usable for full merger evolutions without excision, corotating shifts, or special boundary treatment at the holes.
  • Radiated energy and angular momentum can be read off a complete inspiral-to-ringdown simulation and checked against the remnant's mass and spin, as the paper does with $2.8Â¥% Â¥pm 0.2Â¥%$ of $M$ against a horizon mass deficit of $0.028M Â¥pm 0.002M$.
  • The approach extends naturally to unequal-mass binaries, spinning holes, and more astrophysically motivated initial data such as thin-sandwich or post-Newtonian sets, which the authors list as next steps.
  • Fourth-order convergence at resolutions up to $h = M/36$ on a unigrid code with Fisheye coordinates suggests that the accuracy needed for waveform templates is reachable without adaptive mesh refinement.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • The regularity argument is a leading-order power count at a single point, so the technique's real reach will be decided by how the discretized right-hand side behaves at grid points where $Â¥chi$ is tiny but nonzero; a code that evaluates $Â¥partial_iÂ¥chi/Â¥chi$ naively could reintroduce the divergence the substitution was meant to remove.
  • Because the run reported here stops at about $Â¥sim 85M$ for computational reasons, the paper establishes stability through merger but says nothing yet about the many-orbit inspirals that matched-filtering templates actually need; whether the moving-puncture gauge remains well behaved over tens of orbits is an open question this result sets up rather than answers.
  • The trick of absorbing a singular conformal factor into a variable that is smooth where the original diverges looks portable to other elliptic-plus-evolution problems with point sources, such as boosted or spinning puncture data and possibly higher-dimensional or matter-coupled punctures.
  • The authors' own note that Hamiltonian constraint violation shows cubic convergence near the former puncture positions and only quadratic outside, which they attribute to second-order-accurate initial data, is a concrete loose end: a higher-order elliptic solver for the initial data would test whether the residual constraint growth is an initial-data artifact or a property of the evolution scheme.
Original abstract

We present a new algorithm for evolving orbiting black-hole binaries that does not require excision or a corotating shift. Our algorithm is based on a novel technique to handle the singular puncture conformal factor. This system, based on the BSSN formulation of Einstein's equations, when used with a `pre-collapsed' initial lapse, is non-singular at the start of the evolution, and remains non-singular and stable provided that a good choice is made for the gauge. As a test case, we use this technique to fully evolve orbiting black-hole binaries from near the Innermost Stable Circular Orbit (ISCO) regime. We show fourth order convergence of waveforms and compute the radiated gravitational energy and angular momentum from the plunge. These results are in good agreement with those predicted by the Lazarus approach.

The compiled record

The paper's own sentences, arranged, kept beside the review and separate from it. Each label is drawn from a closed vocabulary; the number beside it is the probability the classifier gave that label. Numbers in the tables are the authors', checked by code. No sentence in this record was written by a model.

Compiled by classifier-1.13-20260917 via OpenRouter, question bank abstract-p1-20260923, 2026-09-23: 1 request, 4,912 input tokens, $0.0002 reported.

What the abstract itself states, as probabilities: claims first 0.08 · states limitation 0.66 · human participants 0.01 · controlled comparison 0.20 · randomized assignment 0.01 · states quantitative result 0.46 · compares with prior results 0.97 · mentions code or data availability 0.02.

The abstract, sentence by sentence

arXiv listing abstract v2. Results and the main claim are underlined; .

We present a new algorithm for evolving orbiting black-hole binaries that does not require excision or a corotating shift. Our algorithm is based on a novel technique to handle the singular puncture conformal factor. This system, based on the BSSN formulation of Einstein's equations, when used with a `pre-collapsed' initial lapse, is non-singular at the start of the evolution, and remains non-singular and stable provided that a good choice is made for the gauge. As a test case, we use this technique to fully evolve orbiting black-hole binaries from near the Innermost Stable Circular Orbit (ISCO) regime. We show fourth order convergence of waveforms and compute the radiated gravitational energy and angular momentum from the plunge. These results are in good agreement with those predicted by the Lazarus approach.

background contribution method result interpretation main claim

What the classifier reports for each sentence. The number is the probability it gave the label.
#sentencerolethe authors' own workstates the main conclusionstrength
s0 We present a new algorithm for evolving orbiting black-hole binaries that does not require excision or a corotating shift. contribution 1.00 0.93 0.44 asserted 1.00
s1 Our algorithm is based on a novel technique to handle the singular puncture conformal factor. method 0.76 0.84 0.08 asserted 1.00
s2 This system, based on the BSSN formulation of Einstein's equations, when used with a `pre-collapsed' initial lapse, is non-singular at the start of the evolution, and remains non-singular and stable provided that a good choice is made for the gauge. method 0.89 0.68 0.21 asserted 0.98
s3 As a test case, we use this technique to fully evolve orbiting black-hole binaries from near the Innermost Stable Circular Orbit (ISCO) regime. method 0.66 0.94 0.09 asserted 0.96
s4 We show fourth order convergence of waveforms and compute the radiated gravitational energy and angular momentum from the plunge. result 0.99 0.96 0.20 asserted 1.00
s5 These results are in good agreement with those predicted by the Lazarus approach. result 0.72 0.70 0.16 asserted 0.90

What kind of paper the classifier reports this is: method or tool 1.00 · evidence from simulation or computation 1.00.

What was examined, and what was not

sources read
  • arXiv listing abstract v2
not examined
  • full text
  • figures and tables
  • references and citing sentences
  • relations to other papers
levels run
  • P abstract
not run
  • S sections
  • T body sentences
  • C citations
  • R cross-paper relations

Every quoted sentence on this record is the authors'. Every label is from Pith's closed vocabulary, with the probability shown beside it. Every number in the measurements table is the authors', checked by code. Record as JSON

Figures
Figure 1. FIG. 1: The trajectories of the punctures (the circles cor
Figure 2. FIG. 2: QC0 waveform. The real and imaginary parts of the

Figures reproduced from arXiv:gr-qc/0511048 by the authors.

References
How later papers used it

414 later papers cite it

Evidence: substantial. 363 of the citing papers share no author with it, and 131 of those carry a confident label from Pith's citation classifier.

Pith read 39 of the 576 citing sentences, choosing only those whose citation markers point at this paper alone, and wrote the paraphrases from them. The summary, the counts, the evidence level, the links and the quotes come from code; the summary counts the papers listed below. A later paper is listed under a label when Pith's citation classifier gives one of its sentences that label at probability 0.8 or more (0.95 for support), Pith's reading finds that the citing paper itself did it, and a second check confirms the paraphrase against the sentence. This section lists uses only. It does not list criticism or disagreement, so their absence here says nothing about how the paper has been received. Papers that share an author with this one are listed apart, under "By the same authors". A quote is the citing authors' own sentence, shown when their paper's licence is CC BY, CC BY-SA or CC0. JSON

Uses its method or tool 29 of 156 papers with this label · 16 by the same authors

By the same authors

Uses its data or values 1 of 5 papers with this label · 1 by the same authors

By the same authors

  • By contrast, as a result of the merger, the time-dependent quadrupole moment of the gravitational potential disappears, leaving only its monopolar component, which is slightly weaker due to the loss of mass associated with the energy radiated in gravitational waves ($∼3%$ of the total for equal-mass, nonspinning binaries [ref] ).

    Relativistic gas accretion onto supermassive black Hole binaries from inspiral through merger (2025) (same authors) · Luminosity · 0.93

How the literature uses it

333 arXiv papers cite this in 486 sentences

Each citing sentence is taken from the citing paper's LaTeX source and typed by Pith's classifier. A label appears at probability 0.8 or more (0.95 for support), a contest at 0.7 or more; the quoted words are the citing authors'. How each label was checked · JSON

Papers citing it · 152 sentences without a confident label

Uses its method or tool 125 papers
Uses its data or values 5 papers
  • We use a catalog of NR waveforms [ref] , produced with the [ref] that describe BBHs with mass-ratios \(1≤q ≤10\), and eccentricities up to \(e_0≤0.18\) fifteen cycles before merger.

    Fusing numerical relativity and deep learning to detect higher-order multipole waveforms from eccentric binary black hole mergers · Numerical relativity catalog · 0.88

  • The kick velocity could reach thousands of $km s^-1$ for special cases [ref] , and in more ordinary conditions the velocity is on average $400 km s^-1$ (e.g. [ref] ).

    Phase shift of gravitational waves induced by aberration · Introduction · 0.90

  • Different astrophysical models predict a variety of velocities for GW sources, where we consider the following three particularly interesting cases: (i) gravitational kicks which in average have a velocity of $400 km s^-1$ [ref] while extreme kicks of up to $5000 km s^-1$ can occur for some special configurations [ref] , (ii) BBHs, BNSs, and BHNSBs orbiting a SMBH in the center of a galaxy [ref] , where the velocity depends on the mass of the central SMBH and the radius of the orbit – for a SMBH with a typical mass of $10^7 M_⊙$ and for an orbital radius of $5 mpc$ the orbital velocity is around $3000 km s^-1$ while for an orbital radius of $100 R_s$ ($R_s = 2GM/c^2$ is the Schwarzschild radius of a BH with mass $M$) we have a velocity of around 7 % the speed of light, and (iii) the peculiar velocity of galaxies induced by the gravitational interaction of the host galaxy with its environ

    Moving gravitational wave sources at cosmological distances: Impact on the measurement of the Hubble constant · Observables for high redshift and fast moving sources · 0.81

Compares with it 4 papers
Cites it as support 2 papers
  • A physical motivation is that the results of numerical relativity [ref] show that the pulse of radiation coming from the merger phase is very short and seems to connect smoothly to the previous inspiral and subsequent ringdown phases.

    The Close-Limit Approximation for Black Hole Binaries with Post-Newtonian Initial Conditions · Introduction · 0.99

  • An encouraging motivation for hoping so is that in other physical settings, like for example the merger of two equal mass Black Holes, the perturbative analysis of the early 70s [ref] in the extreme mass ratio limit was able to single out the important physical elements (i.e., the presence of QNMs), sketching a picture that has been later refined, but substantially confirmed, by numerical relativity simulations [ref] .

    Dynamical excitation of space-time modes of compact objects · Conclusions · 0.98

Cites it as background 149 papers

Discussion

Disagree with the review, or know something it missed? Leave a comment, highlight any passage above to comment on it, or challenge this review formally.