{"id":"6d4e7b54-c3b7-4375-baf1-5c83c82bc4c6","arxiv_id":"gr-qc/0507014","paper_version":1,"verdict":"ACCEPT","confidence":"HIGH","novelty_score":9.0,"correctness_risk":"low","formal_verification":"none","parameter_count":0,"one_line_summary":"A generalized-harmonic-coordinate code carries an equal-mass, non-spinning black hole binary through plunge, merger, and ringdown, leaving a Kerr remnant with spin parameter about 0.70.","lead":"A numerical code evolves two equal-mass black holes through their final orbit, merger, and ringdown in full three-dimensional general relativity, the first time a non-axisymmetric binary has been carried through coalescence. The remnant settles to a spinning Kerr black hole, and the run emits a gravitational waveform of the kind detectors like LIGO would need to match against data.","discovery_kind":"unclear","skeptic_critique":{"model":"selfhosted/deepseek-v4.1-flash","headline":"The remnant parameters (a≈0.70, M≈1.9M0) and the 5% energy estimate rest on three runs whose convergence is asserted but not shown; the convergence tests are deferred to an unpublished 'elsewhere,' and the resolution spread is large enough that the high-resolution numbers are not demonstrably conver","rationale":"The reader's weakest_assumption — horizon-based extraction of a and Mf — is a genuine concern, and the paper's own concession about applying Kerr formulas too early is a red flag. But the more load-bearing issue is one step upstream: the three resolutions are not shown to be converged, and the convergence tests that would establish trustworthiness are deferred to an unpublished 'elsewhere.' The reader notes this deferral in the rationale but does not elevate it to the weakest assumption. I would keep the verdict at ACCEPT because the qualitative existence claim — the code completes the evolution and produces a physically plausible merger waveform — does not depend on convergence, and the paper hedges the quantitative claims ('indications are,' 'rough calculation'). The quantitative numbers should be read as preliminary, which the reader already says about the 5% figure; the same caveat applies to a≈0.70 and M≈1.9M0. The shared axisymmetry assumption behind the two horizon methods is a further reason not to treat their agreement as independent validation.","tokens_in":8540,"tokens_out":11357,"duration_ms":183117,"concrete_test":"Perform a three-level convergence analysis on a robust scalar — e.g., the final BH mass from Table I, or the Ψ4 waveform at a fixed extraction radius — and report the observed convergence order and Richardson-extrapolated error bars. If the observed order is below first order, or if the extrapolated uncertainty on a and Mf exceeds the 6-7% spread across the three runs, the high-resolution remnant parameters are not converged and the a≈0.70 / M≈1.9M0 headline should be downgraded to a rough estimate.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's qualitative claim — that the generalized-harmonic code can carry an equal-mass binary through orbit, merger, and ringdown — is well supported: the runs complete, the horizon tracks look physical, and the waveform has the expected shape. The load-bearing soft spot is the quantitative layer built on top of it. The abstract and conclusion state a≈0.70 and M≈1.9M0 as results, and the abstract quotes ~5% radiated energy. These rest on three resolutions whose final BH masses are 1.77, 1.85, 1.90 M0 and spins 0.74, 0.73, 0.70, with radiated energies at r=25M0 of 7.1%, 5.1%, 4.7% (Table I and footnote [20]). The paper says 'convergence tests suggest that with sufficient resolution the code can evolve the system for as long as needed' (Sec. I) but presents no convergence data; the tests are deferred to 'elsewhere' (Sec. I). Without a measured convergence order, the reader cannot tell whether the high-resolution numbers are approaching a continuum limit or whether the apparent stability is resolution-dependent. The horizon extraction compounds this: the two J/a methods share the same axisymmetry assumption (rotation axis orthogonal to z=0, azimuthal orbits in z=constant), so their agreement is not independent validation of a≈0.70. The paper concedes the timing is questionable ('one might have expected the black hole to still be too far from its stationary state') but not the shared assumption.","agreement_with_reader":"partial"},"referee_report":{"model":"selfhosted/deepseek-v4.1-flash","summary":"The manuscript reports the first fully three-dimensional numerical evolution of a non-axisymmetric binary black hole spacetime through inspiral, plunge, merger and ringdown. The code solves Einstein's equations in a generalized harmonic formulation with constraint-damping terms, uses compactified coordinates with outer boundaries at spatial infinity, adaptive mesh refinement, and dynamical excision. Initial data are produced by two Lorentz-boosted scalar-field profiles that collapse to form the two black holes. A single equal-mass, non-spinning example is evolved with three grid hierarchies; the binary completes roughly one orbit before merging, and the remnant is characterized as a Kerr black hole with a≈0.70 and mass ≈1.9 M0. Gravitational waves are extracted via the Newman-Penrose scalar Ψ4, and the radiated energy is estimated at roughly 5% of 2M0, with the caveat that the estimate is resolution-limited in the wave zone.","tokens_in":8858,"tokens_out":4181,"duration_ms":65849,"significance":"If the result holds, this is a genuine milestone: the last uncomputed segment of the binary black hole waveform — late inspiral, plunge, merger and ringdown — becomes accessible to numerical solution, which is precisely the regime where post-Newtonian and perturbative methods are expected to fail. The paper is candid about what is and is not established, and the qualitative claim (the code survives through coalescence and produces a physically sensible horizon track and waveform) is directly demonstrated by the runs shown in Figs. 1–3 rather than merely asserted. The methodological ingredients — generalized harmonic evolution with constraint damping, excision tracking, compactified outer boundaries — are described concretely enough to be reproduced, and the paper is appropriately framed as a Letter reporting early success rather than a precision measurement. The quantitative layer (remnant parameters, radiated energy) is weaker and is largely flagged as such by the authors themselves.","major_comments":[{"comment":"The two determinations of the remnant angular momentum are presented as mutually corroborating, but they are not independent in the way the text implies. The dynamical-horizon estimate explicitly assumes that the rotation axis is orthogonal to the z=0 orbital plane and that each closed orbit of the azimuthal vector field lies in a z=constant surface; the circumference-ratio method of Eq. (5) likewise assumes the polar axis of the horizon coincides with z. Both methods therefore share the same axisymmetry and axis-alignment assumptions, so their agreement constrains the magnitude of J but not the orientation or the axisymmetry of the remnant. A concrete test would be to repeat the Cr-based estimate using polar circles defined about a tilted axis and report the sensitivity of a to that choice; if a is stable at the few-percent level under such variation, the agreement becomes meaningful evidence rather than a consistency check of a shared premise.","section":"Sec. III, Table I and Fig. 2","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"Equation (5) is a fitting formula calibrated to stationary Kerr horizons, so applying it to a highly distorted apparent horizon shortly after merger and recovering a Kerr-like value is partly built into the procedure. The authors acknowledge the timing concern in the sentence beginning 'Remarkably, the dynamical horizon estimate...', but the concern is not resolved by the agreement of the two numbers, since both are calibrated against stationarity. A more probative test would be to show that the recovered a is constant in time over the ringdown once the horizon has settled, and to quantify the residual deviation of the horizon geometry from Kerr (for example, the size of the higher multipole moments of the horizon mass or curvature) as a function of time. If those residuals decay on the ringdown timescale, the claim that the remnant is Kerr with a≈0.70 is well supported; as presented, the claim rests on a formula whose domain of validity is exactly what is in question.","section":"Sec. III, Eq. (5) and Fig. 2 (right panel)","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"The radiated-energy estimates decrease monotonically with extraction radius — 4.7%, 3.2%, 2.7%, 2.3% at r = 25, 50, 75, 100 M0 in the high-resolution run — and the paper attributes the decrease to numerical dissipation in the low-resolution wave zone. This means the r→∞ extrapolation is not under control, and the abstract's 'on the order of 5%' is in fact the smallest-radius, most contaminated number. The independent horizon-mass-difference estimate (≈5%) is a useful cross-check, but it carries its own systematic from the same numerical loss of mass after merger that the paper notes in the Fig. 2 caption. The abstract should either quote the full range across extraction radii or attribute the 5% explicitly to the horizon-mass method, and the text should state what fraction of the discrepancy between the two methods is attributable to dissipation rather than to the finite-radius validity of Eq. (6).","section":"Sec. III, Eq. (6) and footnote [20]","authors_acknowledge":null,"authors_acknowledge_p":null}],"minor_comments":[{"comment":"The outline states 'We conclude in Sec. V', but the manuscript contains no Sec. IV; the sections run I, II, III, then V. Either a section was removed in preparation or the numbering should be corrected.","section":"Sec. I, final paragraph of the outline","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"The caption defines M0 as the mass of one initial black hole 'measured from the higher resolution simulation', so the high-resolution entry of exactly M0 is definitional rather than a measurement, and only the 0.97M0 and 0.99M0 entries carry information. This should be stated in the caption so the row is not read as three independent determinations.","section":"Table I, 'initial BH masses' row","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"The orbital eccentricity is quoted only as a range, 0–0.2, inferred by overlaying reference ellipses. Since the paper's stated goal is to extract orbital information from the waveform, the residual eccentricity matters for the phasing; the authors should state whether the range reflects genuine uncertainty in the initial data parameters or a fit degeneracy, and whether the measured orbit is consistent with zero eccentricity at the available resolution.","section":"Sec. III, Fig. 1 and Table I","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"The symmetrization notation in the term 2H(α,β) is not defined, and the index structure of the constraint-damping term (in particular the placement of n^γ C_γ relative to g_{αβ}) would benefit from one explicit sentence. The same applies to Eq. (3), where the dimensions of the ξ1 term and the meaning of the η offset in the lapse sector are not explained.","section":"Sec. II, Eq. (1)","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"The caption notes that the waveforms must be shifted by more than the light travel time between extraction radii, which the authors interpret as evidence that the extraction points are too close to the source. It would help to state the implied systematic on the phase of the extracted waveform, since phase accuracy is what matters for matched filtering.","section":"Sec. III, Fig. 3","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"The 5% figure is quoted relative to 2M0, but the ADM mass of the spacetime is 2.39M0 and roughly 15% of the scalar-field energy does not collapse. The abstract's phrase 'initial rest mass of the system' is therefore ambiguous; the denominator should be stated explicitly wherever the percentage appears.","section":"Abstract and Sec. V","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"No comparison is made between the extracted waveform and a post-Newtonian or perturbative prediction in the regime where those should still be approximately valid. Even a qualitative overlay during the early inspiral would strengthen the claim that the numerical waveform is physically faithful rather than merely plausible in shape.","section":"Sec. III, waveform discussion","authors_acknowledge":null,"authors_acknowledge_p":null}],"recommendation":"accept","confidential_remarks":"This is a landmark manuscript and the qualitative claim — that a generalized-harmonic code can carry an equal-mass binary through merger and ringdown — is directly demonstrated rather than merely asserted. The quantitative layer is thinner than the abstract suggests: the convergence evidence is deferred to a follow-up paper, the two spin estimates share an axisymmetry assumption, and the radiated-energy numbers decrease monotonically with extraction radius. None of these undermines the central result, and the authors flag most of them, but the abstract's a≈0.70 and 5% would read better with the resolution spread attached. I would not hold publication for this; the follow-up paper is the right place for the convergence study."},"author_rebuttal":null,"desk_editor":{"model":"selfhosted/deepseek-v4.1-flash","letter":"You already know the headline: this is the paper that first got a full 3D, non-axisymmetric binary black hole through plunge, merger, and ringdown. That claim has held up. The ingredients are borrowed — generalized harmonic formulation from Garfinkle, constraint damping from Gundlach et al., dynamical excision, compactified outer boundaries, AMR — but the combination is the contribution, and it works. The paper is also unusually honest about what it can't yet do: it says outright that resolution far from the binary prevents an accurate energy estimate, that the amplitude falloff with radius is numerical dissipation, that the extraction radii may be too close, and that the horizon formulas may be applied before the remnant is stationary. That candor is worth crediting, not penalizing. The stress-test's central worry is fair, though: the abstract's 'with sufficient resolution' is backed by convergence tests that are deferred to 'elsewhere' (Sec. I). For a letter that's normal, but it means the resolution study is asserted, not shown, and the spread across the three runs — final masses 1.77/1.85/1.90 M0, spins 0.74/0.73/0.70, radiated energy at r=25M0 of 7.1/5.1/4.7% — is wide enough that the high-resolution numbers aren't demonstrably converged. The 5% figure should be read as order-of-magnitude; the paper's own numbers span 2.1 to 7.1% with no quoted uncertainty. On the spin: the stress-test is right that the two J estimates share the axisymmetry assumption (rotation axis orthogonal to z=0, azimuthal orbits in z=constant), so their agreement is consistency, not independent validation. The paper does flag the assumption as needing confirmation, so it's acknowledged rather than hidden, but the point stands. The Kerr-calibrated Eq. (5) applied to a still-ringing horizon is the weakest single step, partly offset by the dynamical-horizon estimate. None of this touches the existence result. Who should read it: numerical relativists, and anyone building waveform models who needs to know the merger segment is now computable. It deserves a serious referee — accept.","headline":"The paper that first carried a non-axisymmetric binary through merger and ringdown; the qualitative claim has held up for twenty years, the quantitative numbers are rougher than the abstract's tone suggests.","tokens_in":9512,"tokens_out":1437,"would_cite":true,"duration_ms":22995,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"pith_extraction":{"msc":[],"pacs":["04.25.Dm","04.30.Db","04.70.Bw"],"model":"selfhosted/deepseek-v4.1-flash","headline":"A numerical code evolves two black holes through their final orbit, merger, and ringdown, producing the first complete strong-field waveform from a non-axisymmetric binary.","keywords":["numerical relativity","binary black holes","generalized harmonic coordinates","constraint damping","gravitational waves","black hole merger","apparent horizon","adaptive mesh refinement"],"falsifier":"A higher-resolution run with the wave zone properly resolved, extracting the radiated energy at several large radii and showing the answers converge to a single value, would confirm the ~5% figure; if instead the extracted energy keeps falling with radius and never settles, or if a better-resolved evolution fails before merger, the central claim of a usable merger waveform weakens accordingly.","tokens_in":9827,"feed_emoji":"🕳️","tokens_out":2288,"duration_ms":42730,"temperature":0.7,"pith_summary":"The last uncomputed stretch of a binary black hole's life is the one that matters most for gravitational wave detection: the final orbits, the plunge, the merger, and the ringdown, where perturbation theory stops working and only a full numerical solution of Einstein's equations will do. This paper reports that a code built on generalized harmonic coordinates can evolve an equal-mass, non-spinning binary through exactly that stretch without the simulation falling apart. The author runs one configuration through roughly a single orbit before coalescence, continues past merger, and reads off a remnant black hole with mass about $1.9M_0$ and spin parameter $a$ near $0.70$, plus a rough estimate that about 5% of the initial rest mass leaves as gravitational waves. The point is not the particular numbers, which the paper treats as provisional, but the demonstration that the scheme survives long enough, at achievable resolution, to hand over the waveform segment that analytic methods cannot supply.","feed_headline":"Binary black holes evolved through merger and ringdown","feed_subtitle":"A generalized harmonic code follows an equal-mass pair through a plunge orbit, leaving a spinning remnant and a first waveform.","key_machinery":"The load-bearing mechanism is the generalized harmonic formulation with constraint damping. Writing the field equations so that the coordinates satisfy a wave-like condition sourced by freely specifiable functions $H_μ$ turns the gauge choice into part of the evolution system rather than an external constraint; the residual constraint violation $C_μ$ is then damped by adding terms proportional to a positive constant $κ$ to the field equations, following the damping scheme the author credits to a Z4 treatment. Around this sit the practical supports that make a binary survivable: a compactified domain with physically correct outer boundary conditions, mesh refinement concentrated near the holes, excision boundaries that move with the apparent horizons, a slicing condition that slows the collapse of the lapse, and numerical dissipation. The paper's own assessment is that the damping terms had little effect without black holes present but a significant effect on how long a black hole evolution stays accurate.","core_discovery":"The claim is that a specific numerical scheme can carry a genuinely three-dimensional, non-axisymmetric binary black hole spacetime through coalescence and ring down, which the paper states had not been done before. The scheme combines several ingredients: field equations written in a generalized harmonic form with freely chosen source functions encoding the gauge, a damping term added to the constraints that the author reports makes a large difference to how long a black hole simulation runs accurately, a spatially compactified grid whose outer boundary sits at infinity, adaptive mesh refinement, and dynamical excision that follows each hole through the grid. The demonstration case is a binary of two equal-mass, non-spinning holes built from collapsing scalar field pulses, with initial proper separation about $16.6M_0$ and orbital angular velocity about $0.023/M_0$, evolved at three resolutions. The remnant is characterised two ways: from horizon area and angular momentum through the Smarr formula, and from the ratio of polar to equatorial horizon circumference through a Kerr fitting formula, with the angular momentum itself taken from a dynamical horizon calculation. Both routes give $a$ between about 0.70 and 0.74, and the paper notes their agreement even shortly after merger, when neither method would obviously be expected to apply.","pith_inferences":["The paper's own energy estimate varies by roughly a factor of two across extraction radii (4.7% down to 2.3%), which is the signature of dissipation eating the wave amplitude with distance; a converged answer would need the wave zone resolved as well as the holes, and the author says as much.","If the scheme generalises, the natural next tests are unequal masses and non-zero spins, since spin-orbit and spin-spin couplings change both the waveform phasing and the remnant spin, and neither is exercised here.","The agreement between the two independent spin estimators so soon after merger is a check worth pushing harder: if they continue to agree in better-resolved runs, the circumference-ratio formula becomes a cheap diagnostic for near-stationary horizons in dynamical settings.","The scalar field construction of initial data is a workaround for the lack of a clean binary black hole initial data prescription; a reader could reasonably expect the method's real test to come once puncture-style or excision-compatible data are fed to it."],"forward_implications":["Gravitational wave searches gain access to the merger and ringdown segment of the waveform, the part matched filtering needs and perturbation theory cannot supply.","The radiated energy and final spin of an equal-mass non-spinning merger become quantities a numerical code can compute and other methods can be checked against.","The generalized harmonic formulation with constraint damping becomes a candidate platform for a wider family of binary configurations rather than a one-off.","Because the outer boundary is placed at spatial infinity, wave extraction is not contaminated by an artificial outer edge, in principle.","Constraint damping is identified as a practical lever on simulation lifetime, separate from resolution."],"supporting_citations":[{"why":"The earlier description of this code and its method, which the present letter extends with constraint damping and reports results from.","marker":"[9]"},{"why":"The suggestion to formulate the field equations in harmonic coordinates, which the evolution system is built on.","marker":"[10]"},{"why":"The constraint-damping scheme adapted here, which the author identifies as decisive for simulation lifetime with black holes present.","marker":"[13]"},{"why":"The dynamical horizon framework used to estimate the remnant angular momentum.","marker":"[15]"},{"why":"The polar-to-equatorial circumference ratio formula used as the second, independent estimate of the Kerr spin parameter.","marker":"[16]"},{"why":"The flux formula used to convert the extracted Newman-Penrose scalar into a radiated energy estimate.","marker":"[17]"},{"why":"The standard initial data techniques on which the conformally flat, maximal, harmonic starting slice is built.","marker":"[14]"}],"fun_headline_variants":["A generalized harmonic code carries equal-mass black holes through merger","An equal-mass black hole binary evolved through plunge, merger, and ringdown","The remnant spin is near a=0.70 from horizon geometry and the Smarr formula","Far-field resolution limits the gravitational-wave energy estimate to roughly 5%"],"cache_read_input_tokens":0,"weakest_assumption_plain":"The remnant's mass and spin are read off a highly distorted, still-changing apparent horizon using formulas derived for stationary black holes, together with an angular momentum estimate that assumes the rotation axis is perpendicular to the orbital plane and that the horizon's rotational orbits lie in flat slices.","fun_headline_variants_meta":{"raw":{"variants":["A generalized harmonic code carries equal-mass black holes through merger","An equal-mass black hole binary evolved through plunge, merger, and ringdown","The remnant spin is near a=0.70 from horizon geometry and the Smarr formula","Far-field resolution limits the gravitational-wave energy estimate to roughly 5%"]},"model":"selfhosted/deepseek-v4.1-flash","effort":"low","attempts":1,"cost_usd":0.0,"raw_usage":{"total_tokens":2966,"prompt_tokens":920,"completion_tokens":2046,"prompt_tokens_details":null,"completion_tokens_details":{"reasoning_tokens":1976}},"tokens_in":920,"tokens_out":2046,"duration_ms":38902,"image_pages":[],"temperature":1.0,"reasoning_budget":50,"reasoning_tokens":1976,"transport_provider":"selfhosted","cache_read_input_tokens":0,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-10-29T16:12:37.828081+00:00","model_set":{"reader":"selfhosted/deepseek-v4.1-flash"},"falsifier":"A higher-resolution run with the wave zone properly resolved, extracting the radiated energy at several large radii and showing the answers converge to a single value, would confirm the ~5% figure; if instead the extracted energy keeps falling with radius and never settles, or if a better-resolved evolution fails before merger, the central claim of a usable merger waveform weakens accordingly.","supporting_citations":[{"cited_title":"Class. Quant. Grav. 22 425,","cited_arxiv_id":null,"evidence_quote":"The earlier description of this code and its method, which the present letter extends with constraint damping and reports results from."},{"cited_title":"Phys.Rev. D65, 044029","cited_arxiv_id":null,"evidence_quote":"The suggestion to formulate the field equations in harmonic coordinates, which the evolution system is built on."},{"cited_title":"Constraint damping in the Z4 formulation and harmonic gauge","cited_arxiv_id":"gr-qc/0504114","evidence_quote":"The constraint-damping scheme adapted here, which the author identifies as decisive for simulation lifetime with black holes present."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"The dynamical horizon framework used to estimate the remnant angular momentum."},{"cited_title":"Phys.Rev. D52, 870","cited_arxiv_id":null,"evidence_quote":"The polar-to-equatorial circumference ratio formula used as the second, independent estimate of the Kerr spin parameter."},{"cited_title":"in Sources of Gravitational Radiation , ed. L. Smarr, Seattle, Cambridge University Press","cited_arxiv_id":null,"evidence_quote":"The flux formula used to convert the extracted Newman-Penrose scalar into a radiated energy estimate."},{"cited_title":"Living Rev.Rel. 3, 5","cited_arxiv_id":null,"evidence_quote":"The standard initial data techniques on which the conformally flat, maximal, harmonic starting slice is built."}],"review_version":1}