{"id":"20b4f0cb-c34f-4e38-a3c3-a09fed9bc1b0","arxiv_id":"1908.06034","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"Quasi-maximal slicing, adding a maximal-slicing source term to 1+log slicing, removes the coordinate singularities that break moving puncture evolutions of near-critical Brill waves.","lead":"The authors introduce a new coordinate condition, 'quasi-maximal' slicing, that prevents the standard moving puncture gauge from crashing when simulating gravitational collapse of Brill waves. This enables longer and more complete numerical evolutions of black hole formation in axisymmetric vacuum spacetimes.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"At A=5, the main claim rests on a single filtered, time-extrapolated run with no resolution or filter-sensitivity check; the cure might be an artifact of the approximate solver rather than of quasi-maximal slicing.","rationale":"The reader's weakest_assumption already pointed to the approximate, filtered, time-interpolated Wa and the lack of A=5 filtering checks. I agree with that general direction, but I would put the weight on a sharper evidential gap: the central claim is made specifically for A=5, yet every quantitative validation of the quasi-maximal procedure is performed at A=0.1 or A=1. The A=5 section reports a single successful run at one resolution and one filter setting, with no convergence study and no filter-sensitivity study in the regime where the claim is being made. This does not make the claim false, and the paper is honest about the shift shock and the inability to evolve indefinitely, but it does mean the headline result is less secure than the presentation suggests. A positive outcome of the proposed resolution and filter-dependence test would substantially strengthen the paper; a negative outcome would show that the cure is an artifact of the numerical recipe. Because this concern is a request for evidence rather than a demonstrated contradiction, the reader's CONDITIONAL verdict remains appropriate.","tokens_in":14782,"tokens_out":7222,"duration_ms":81624,"concrete_test":"Run the A=5 case at two additional resolutions (for example 1024 and 4096 points per refinement level) and with at least one different filter exponent d (for example d=4 and d=8), keeping all other parameters fixed; compare K_max, the apparent-horizon mass, and the invariant χ in the equatorial plane. If these quantities differ by more than the stated ~10^-3 horizon-mass uncertainty, the A=5 cure is not a converged property of quasi-maximal slicing itself but of the specific numerical recipe used to approximate it.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The headline claim is specifically about A=5 Brill waves, but the convergence and sensitivity evidence in §4.2 is for A=0.1 and A=1. For A=5, the paper relies on one 2048-point-per-level zero-shift run that reaches t≈25. The quasi-maximal source term Wa is, by the authors' own description in §4.1, a 'very approximate' solution of the elliptic equation (21), built from time-extrapolated spectral coefficients and then heavily low-pass filtered via (25) with empirically chosen J=36, d=6. The switch back to 1+log at t0 in (26) is also empirically determined. A run that does not crash is not by itself evidence that the modified gauge cured a slicing singularity: the filter, the time extrapolation, or the switch schedule could be the actual stabilizers. Because the filtering and the proximity of Wa to the true quasi-maximal driver are only validated in the mild A=1 regime, there is no demonstrated connection between the A=5 result and the intended maximal-slicing behavior of equation (22).","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes a new 'quasi-maximal' slicing condition (Eq. 22), obtained by adding a source term W_a--an approximate solution of the elliptic equation (21)--to the 1+log slicing condition, in order to cure coordinate singularities that appear in moving-puncture evolutions of near-critical Brill wave data. The authors implement this in a BSSN code with the analytic Cartoon method, computing W_a with a pseudospectral elliptic solver, time extrapolation/interpolation, and a low-pass filter (25), and switching back to 1+log at an empirically chosen time t0 (26). They reproduce the known failure of 1+log slicing at A=5, show that max|K| tends to zero with the time step for A=0.1, that K is reduced at A=1, and that gauge-invariant quantities converge to fourth order under grid refinement at A=1. Their main result is a single A=5 run with zero shift that reaches t≈25, forms an apparent horizon, and allows construction of an event horizon and conformal diagrams; they also report long runs at A=4.6, 4.65, 4.8, 5.0, and 5.5. The paper acknowledges that with the Γ-driver the A=5 shift develops a shock and that the A=5 success is therefore limited to zero shift.","tokens_in":15042,"tokens_out":8431,"duration_ms":76098,"significance":"If the central claim is correct, the quasi-maximal slicing condition is a practically valuable and computationally cheap modification of 1+log slicing that enables long-term simulations of Brill wave collapse in a finite-difference moving-puncture code, and the gauge-invariant diagnostics (Kretschmann mass, circumferential radius, chi, conformal diagrams) are useful tools for studying the approach to critical collapse. The paper's strengths include a clearly motivated construction (W_a's origin in maximal slicing is stated explicitly, so the derivation is not circular), several nontrivial checks--Δt-convergence of K at A=0.1, reduced K at A=1, fourth-order convergence of invariants at A=1, and reproduction of the known A=5 failure--and the release of the simulation code. However, the headline A=5 claim rests on a single zero-shift run with no resolution or filter-sensitivity study at that amplitude, and the well-posedness of the modified gauge is justified only by linearization around flat space. The significance is therefore conditional on additional evidence at the amplitude where the claim is made.","major_comments":[{"comment":"The central claim that quasi-maximal slicing cures the A=5 singularities is supported by a single run. The convergence evidence in Sec. 4.2 and Figs. 3-4 is for A=0.1 and A=1 only; no resolution study is reported for A=5. Because W_a is, in the authors' words, a 'very approximate' solution of (21), and because the low-pass filter (25) with J=36, d=6 and the switch time t0 in (26) are empirically chosen, the run's survival could be an artifact of dissipation, of the time extrapolation, or of the switching schedule rather than of the intended maximal-slicing driver. A resolution study at A=5 (grid spacing, CFL factor, filter parameters J and d, and t0) is required to tie the observed stability to the quasi-maximal source term.","section":"Section 5, 'Near-critical case: A=5'"},{"comment":"The well-posedness argument is insufficient for the strong-field regime. The statement that the source term 'does not modify the principal part of the PDEs linearized around flat space' only addresses perturbations of Minkowski spacetime; near-critical A=5 Brill waves develop large K and strong lapse dynamics, and the actual W_a is time-interpolated and low-pass filtered. No argument is given that the nonlinear system with this approximate, time-dependent driver remains well-posed. At minimum, the A=5 run should be accompanied by a report of Hamiltonian and momentum constraint violations under resolution, or another diagnostic demonstrating that the modified gauge does not introduce an instability.","section":"Section 3, after Eq. (22)"},{"comment":"The demonstration is limited to zero shift. The authors state that with the Γ-driver the shift develops a shock profile around t≈10, so they set the shift to zero; the moving-puncture failure reported in [11] concerns the full gauge, including the Γ-driver. Thus the paper does not establish that quasi-maximal slicing cures the moving-puncture failure for A=5; it establishes that the modified slicing improves the lapse sector in a zero-shift run. This limitation should be stated explicitly, and if the claim is intended for the full moving-puncture gauge, further work on the shift sector is needed.","section":"Section 5, 'Near-critical case: A=5'"}],"minor_comments":[{"comment":"In the last paragraph, 'rougly' should be 'roughly'.","section":"Section 4.1, Quasi-maximal slicing"},{"comment":"The symbol K is used both for the trace of the extrinsic curvature (7b) and for the Kretschmann scalar (12); the same symbol in adjacent sections is confusing, and the Kretschmann invariant should be denoted differently (e.g., \\mathcal{K}).","section":"Section 2.3"},{"comment":"The y-axis label 'Ln \\bar\\rho^2/\\sigma' in the left panel is likely meant to denote the normal derivative L_n \\bar\\rho^2, but it reads like a natural logarithm; the subscript should be typeset clearly.","section":"Figure 2"},{"comment":"The description of the elliptic solver reuse (S=1024) does not state whether the LU decomposition is refreshed during the S solves; since the metric changes with time, the preconditioner quality will degrade, and a sentence on this would help reproducibility.","section":"Section 4.1, Quasi-maximal slicing"},{"comment":"The basis uses \\arccot without specifying the branch or range; for a spectral basis this convention matters and should be stated.","section":"Eq. (24a)"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is well within the scope of the journal and the numerical work appears careful, but the headline result is less solid than the abstract suggests because it rests on a single A=5 run with no resolution or filter-sensitivity study at that amplitude. The requested additional runs are standard and feasible; I would be happy to see a revised version."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: this is a real advance, probably right, and worth sending to referees. The genuinely new piece is the quasi-maximal slicing condition (22), with a source term W_a obtained from the time derivative of the maximal slicing equation. That is not in the prior literature. The paper shows that this condition keeps K smaller than 1+log slicing, and at A=5 it avoids the cusp-like slicing singularities that kill the moving puncture gauge. It also produces the first long-lived axisymmetric Brill-wave evolutions of this kind, including conformal diagrams and a claimed toroidal event horizon. That combination matters.\n\nStrengths: the evidence is honest and layered. At A=0.1 the maximum norm of K goes to zero with time step, which directly checks the intended target ∂_tt K = 0. At A=1, gauge-invariant quantities converge at fourth order under grid refinement. The reproduction of the known 1+log failure for A=5 is consistent with [11], and the apparent horizon masses agree with the independent pseudospectral results in [13]. The code is available. The authors also disclose the shift-shock problem and the limitations of their spectral solver rather than hiding them. Citation practice is fair: Hilditch et al. are credited and compared throughout.\n\nSoft spots, in proportion: the headline A=5 claim rests essentially on one filtered, time-extrapolated run, with no resolution study or filter-sensitivity check at that amplitude. The stress-test note has a real point here. The filter, the approximate elliptic solve, and the empirically chosen switch time could all be doing stabilizing work that is not exactly the intended quasi-maximal driver. That said, the same approximate driver demonstrably pushes K down at A=1, so the mechanism is at least connected. And the A=5 run reaches t≈25, forms an apparent horizon, and gives a mass compatible with [13]; I do not think the result is a numerical mirage, just less well supported than it could be.\n\nOther soft spots are minor. The toroidal event horizon is reported almost in passing for A=4.8, apparently without a dedicated convergence check; I would treat it as suggestive rather than established. Well-posedness of the modified gauge is argued only from the linearized principal part, which is acceptable for a numerical paper but not a proof. The zero-shift workaround is a clear limitation, and the authors say so.\n\nBottom line: anyone working on vacuum critical collapse or on gauge conditions in numerical relativity will get value from this. It deserves a serious referee, and I would engage with it myself.","headline":"Quasi-maximal slicing is a plausible new cure for Brill-wave gauge failures; the A=5 evidence is thinner than the abstract suggests, but the claim likely holds and deserves review.","tokens_in":15547,"tokens_out":1977,"would_cite":true,"duration_ms":21802,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Adding a maximal-slicing source term to the 1+log lapse cures the slicing blowups that stop near-critical Brill wave collapse evolutions.","keywords":["numerical relativity","Brill waves","quasi-maximal slicing","1+log slicing","moving puncture gauge","gravitational collapse","black hole formation","BSSN formulation"],"falsifier":"Run the $A=5$ Brill wave with the low-pass filter removed and the pseudo-spectral solver for $W_a$ made substantially more accurate; if the resolution-dependent kink in the normal derivative of $\\bar\\rho^2$ reappears before $t\\approx10$, the cure is an artifact of the approximate solver rather than a property of quasi-maximal slicing.","tokens_in":14594,"feed_emoji":"🕳️","tokens_out":9736,"duration_ms":88407,"temperature":0.7,"pith_summary":"This paper tries to establish that a small modification to the standard '1+log' lapse condition removes the coordinate singularity that has prevented numerical evolutions of near-critical Brill wave data from running to completion. The modification adds a source term derived from the maximal slicing condition, producing what the authors call quasi-maximal slicing. If the claim holds, existing moving-puncture codes can evolve strong gravitational wave packets long enough to see apparent horizons, event horizons, and the settling of the newly formed black hole to Schwarzschild, without excision or a different evolution formulation.","feed_headline":"Quasi-maximal slicing cures Brill wave collapse blowups","feed_subtitle":"Adding a maximal-slicing source term to the 1+log gauge lets near-critical collapse runs form a black hole","key_machinery":"The load-bearing object is the quasi-maximal slicing condition (22), $$(\\partial_t-\\mathcal{L}_\\$\\beta$)\\$\\alpha$=-2\\$\\alpha$ K+\\kappa W_a,$$ which appends to the 1+log slicing law a source term $W_a$ constructed as an approximate solution of the elliptic equation (21), derived by differentiating the maximal slicing condition in time. In practice $W_a$ is computed by a pseudo-spectral elliptic solver run once per coarse time step, linearly extrapolated and interpolated between solves, and multiplied by a low-pass filter (25) that damps basis functions once they are no longer resolved by the finite-difference grid. Because the auxiliary equation has trivial principal part $\\Delta W=0$, the added term does not change the principal part of the evolution system linearized around flat space, so hyperbolicity is inherited from 1+log slicing while the source term pushes $K$ toward zero.","core_discovery":"Brill wave data with amplitude $A=5$ evolve under the moving puncture gauge for only about $t\\approx4.9$ before the constant-time slices develop a cusp-like nonsmooth feature, visible as a resolution-dependent kink in the normal derivative of the circumferential radius $\\bar\\rho^2$; the simulation then fails. The central discovery is that replacing the 1+log slicing law $(\\partial_t-\\mathcal{L}_\\beta)\\alpha=-2\\alpha K$ with the quasi-maximal condition $(\\partial_t-\\mathcal{L}_\\beta)\\alpha=-2\\alpha K+\\kappa W_a$, where $W_a$ is a rough, low-pass-filtered, time-interpolated solution of the elliptic equation (21), removes this pathology. With zero shift the $A=5$ run passes $t=10$, an apparent horizon of mass $M_{\\mathrm{AH}}\\approx0.56$ appears, and the spacetime can be described by gauge-invariant quantities; for $A=5.5$ the $\\Gamma$-driver shift works and mass estimates converge to the apparent horizon mass as the hole settles to Schwarzschild.","pith_inferences":["The same 'quasi-' construction should extend to other elliptic gauge conditions, for example deriving a source term from harmonic or minimal-distortion conditions, to stabilize hyperbolic gauges in other strong-field regimes; the paper does not test this.","Because a low-pass filter confining $W_a$ to the central region suffices to cure the failure, the 1+log singularity appears to originate in lapse dynamics inside the strong-field core rather than in outer-boundary gauge behavior.","The shift shock seen at $A=5$ suggests the remaining obstacle to full moving-puncture capability for near-critical data sits in the shift condition, so an analogous quasi-maximal driver for the shift is the natural next step.","The dimensionless invariant $\\chi$, whose first negative region at Brill waves points 'the wrong way' relative to Schwarzschild, could be monitored over successive echo periods in a future critical-collapse search, since a discretely self-similar spacetime should repeat its range of values; the paper notes the possibility but does not pursue it."],"forward_implications":["Near-critical Brill wave data can be evolved with BSSN finite-difference moving-puncture codes far enough to form and track an apparent horizon, which 1+log slicing alone does not allow.","Quasi-maximal slicing with fixed solver parameters converges at fourth order in the time integrator to the same invariant spacetime, so the approximate source term does not alter the physics being simulated.","For supercritical amplitude $A=5.5$, late-time gauge-invariant mass estimates approach the final apparent horizon mass, showing the collapsed object settles to a Schwarzschild black hole.","Even away from criticality the gauge allows compactified conformal diagrams of the collapse, and reveals an early event horizon that is disk-like with a non-smooth rim, and transiently toroidal for $A=4.8$."],"supporting_citations":[{"why":"Defines the Brill wave initial data family and the seed function (2) that all simulations in the paper evolve.","marker":"[2]"},{"why":"Documents the failure of the moving puncture gauge for near-critical Brill waves and supplies the A=5 baseline this paper reproduces and cures.","marker":"[11]"},{"why":"Introduces the 1+log slicing condition whose lapse law is modified by the new source term.","marker":"[9]"},{"why":"Provides an earlier maximal-slicing treatment of gravitational wave collapse that motivates pushing the lapse toward maximal slicing.","marker":"[7]"},{"why":"Constrained axisymmetric evolutions of prolate Brill waves with maximal slicing, another precedent for the gauge choice.","marker":"[14]"},{"why":"Independent pseudospectral evolutions of centered Brill waves whose horizon masses and critical amplitude bounds are compared with the quasi-maximal runs.","marker":"[13]"},{"why":"Shows maximal slicing is only well posed when the K=0 constraint is enforced, justifying the approximate quasi-maximal route.","marker":"[18]"},{"why":"Provides the Gamma-driver shift condition used for the A=5.5 runs and relevant to the shift shock observed at A=5.","marker":"[19]"}],"fun_headline_variants":["New slicing tames Brill wave collapse","Quasi-maximal gauge beats 1+log fail","Gauge fix lets Brill collapse form black hole","Elliptic source term stabilizes collapse","Simple slice law saves critical runs"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The cure relies on the assumption that a rough, low-pass-filtered, time-interpolated solution $W_a$ of the auxiliary elliptic equation stays close enough to true maximal slicing to suppress the singularity at the near-critical amplitude where 1+log fails, a behaviour checked numerically at low amplitude but not proven generally.","fun_headline_variants_meta":{"raw":{"variants":["New slicing tames Brill wave collapse","Quasi-maximal gauge beats 1+log fail","Gauge fix lets Brill collapse form black hole","Elliptic source term stabilizes collapse","Simple slice law saves critical runs"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000167,"raw_usage":{"total_tokens":1249,"prompt_tokens":931,"completion_tokens":318,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":547,"completion_tokens_details":{"reasoning_tokens":252}},"tokens_in":547,"tokens_out":318,"duration_ms":3820,"temperature":1.0,"reasoning_tokens":252,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T12:58:28.153705+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run the $A=5$ Brill wave with the low-pass filter removed and the pseudo-spectral solver for $W_a$ made substantially more accurate; if the resolution-dependent kink in the normal derivative of $\\bar\\rho^2$ reappears before $t\\approx10$, the cure is an artifact of the approximate solver rather than a property of quasi-maximal slicing.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Defines the Brill wave initial data family and the seed function (2) that all simulations in the paper evolve."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Documents the failure of the moving puncture gauge for near-critical Brill waves and supplies the A=5 baseline this paper reproduces and cures."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Introduces the 1+log slicing condition whose lapse law is modified by the new source term."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides an earlier maximal-slicing treatment of gravitational wave collapse that motivates pushing the lapse toward maximal slicing."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Constrained axisymmetric evolutions of prolate Brill waves with maximal slicing, another precedent for the gauge choice."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Independent pseudospectral evolutions of centered Brill waves whose horizon masses and critical amplitude bounds are compared with the quasi-maximal runs."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Shows maximal slicing is only well posed when the K=0 constraint is enforced, justifying the approximate quasi-maximal route."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the Gamma-driver shift condition used for the A=5.5 runs and relevant to the shift shock observed at A=5."}],"review_version":1}