{"id":"c8db07ab-13ba-4049-8cfa-8636487998a4","arxiv_id":"2506.03078","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"A tailored gauge, constraint addition and rescaling recipe makes the strongest hyperboloidal compactification (n=2) converge in a spherical dual-foliation generalized harmonic formulation, recovering scalar quasinormal modes and tails at null infinity.","lead":"This numerical relativity paper shows that the strongest hyperboloidal compactification parameter (n=2), previously inaccessible in the dual-foliation generalized harmonic gauge formulation, can be made to run stably with a carefully chosen gauge, constraint addition, and field rescalings. The payoff is direct extraction of signals, here quasinormal ringing and power-law tails of a scalar field, at future null infinity.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Boundedness of the second-order compactification coefficient hat_C+ is assumed rather than established; the reported runs do not adequately test it, so the n=2 claim remains conditional.","rationale":"The reader's weakest-assumption analysis identifies exactly the load-bearing concern: the regularity of the n=2 scheme relies on hat_C+ remaining bounded, and the paper explicitly states this is not guaranteed by [20]. This is the right point to stress because the whole rescaling strategy, and the regularity of the Jacobian in particular, depends on it. I agree with the reader's CONDITIONAL verdict: the numerical evidence is credible and the recovered physics (QNM, Price tail, Bondi-mass behavior) is a positive, falsifiable check, but the absence of an explicit evolution code and the lack of a direct test of the hat_C+ ansatz prevent a stronger verdict. I would keep the verdict unchanged rather than moving it, because the paper itself is transparent about the test nature of the run and the results support the claim for the presented configuration. The concrete test I propose would settle the concern by directly monitoring the quantity whose boundedness is in question, over a broader range of amplitudes and times, and by checking whether the Jacobian factor remains regular with resolution.","tokens_in":13646,"tokens_out":4886,"duration_ms":64936,"concrete_test":"Reconstruct hat_C+ from the evolved C+ field and Eq. (33) (equivalently from the reduction constraint (34)) at grid points near r_I, and plot its time dependence for the existing run. Then repeat the evolution with amplitudes 10^-3 and 10^-2 and with integration time extended by at least a factor of 5, monitoring the Jacobian factor R'(1 - H'C+) in Eq. (32). If |hat_C+| shows secular growth or the Jacobian factor drifts toward zero or diverges with increasing resolution, the n=2 rescaling is not robust beyond the single tested case.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim that n=2 compactification yields regular, self-convergent evolutions depends on the Jacobian factor (32) staying O(1). To make this factor explicitly regular, the paper introduces the second-order coefficient hat_C+ in Eq. (33). The paper itself states in Section III.C that the Good-Bad-Ugly-F heuristic results of [20] do not guarantee that hat_C+ remains bounded during evolution, and that the presented run should be interpreted as a test. If |hat_C+| grows in time, the denominator in Eq. (32) can become formally singular, and the n=2 rescaling loses its regularity control. The single reported numerical run uses one amplitude (10^-4), one value of r_I (20), and a limited integration time; no time series or monitor of hat_C+ (or of the Jacobian factor) is reported, and the convergence study in Fig. 2 uses a global norm that need not expose growth of a boundary coefficient. Thus the headline conclusion is conditioned on an unproven asymptotic ansatz, exactly as the authors flag. The physical content recovered (QNM ringing, t^-2 tail, Bondi-mass monotonicity) is encouraging and shows the scheme works for the tested configuration, but it does not establish that hat_C+ remains bounded generally.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies a spherically symmetric reduction of the Dual-Foliation Generalized Harmonic Gauge (DF-GHG) formulation of general relativity, with the goal of using the strongest hyperboloidal compactification parameter n=2. The authors choose gauge source functions, constraint additions, and field rescalings to minimize formally singular terms, and present numerical evolutions of small perturbations of Schwarzschild with future null infinity included. They report self-convergence with second-order accuracy, and for constraint-satisfying initial data they extract a Bondi mass that decreases monotonically and a scalar field at I+ that exhibits quasinormal ringing followed by a t^{-2} tail, matching linear theory.","tokens_in":13921,"tokens_out":6877,"duration_ms":80789,"significance":"If the n=2 compactification is indeed regular for the DF-GHG system, this is a useful technical advance: it makes the metric-based hyperboloidal approach directly analogous to conformal compactification and opens a path toward extracting radiation at I+ in 3D. The paper is careful in constructing the rescaled system, provides the reduced equations and the reduction constraints, and ships Mathematica notebooks as auxiliary material. The numerical validation is meaningful: the convergence factors trend toward the expected value 2 with resolution, and the extracted physical outputs (quasinormal frequency, Price tail exponent, Bondi-mass monotonicity) are the right targets. However, the headline claim rests on an assumption that the paper itself identifies as not guaranteed by the underlying asymptotic analysis.","major_comments":[{"comment":"The central claim that the n=2 system is regular depends on the coefficient \\hat{C}_+ in the rescaling (33) remaining bounded, so that the Jacobian factor (32) stays O(1). The paper explicitly states in Section III.C that the Good-Bad-Ugly-F heuristics of [20] do not guarantee this and that the run should be interpreted as a test. As a referee, I flag this as a self-identified missing support: the paper never reports the outcome of that test. No time series of \\hat{C}_+ or of the Jacobian denominator is shown, and the global norm (41) used in the convergence study need not expose growth localized at I+. The reported run also uses a single amplitude (10^{-4}), a single r_I=20, and the convergence study covers roughly 500 time units. Please add a diagnostic of \\hat{C}_+ (or the Jacobian factor) over the full evolution, or otherwise provide evidence that it remains bounded; without this, the n=2 regularity conclusion is conditional.","section":"Section III.C (Eqs. 32-33) and Section IV"},{"comment":"The statement that the rescaling 'should be interpreted as a test of whether \\hat{C}_+ remains bounded' is not made operational: no criterion is given for what would constitute failure of the test (e.g., a threshold on |\\hat{C}_+| or on the Jacobian factor), and no failure mode is analyzed. Consequently the reader cannot tell whether the reported run would distinguish between a bounded and an unbounded \\hat{C}_+. A concrete monitoring prescription and a statement of what was observed in the presented evolution should be added.","section":"Section III.C"}],"minor_comments":[{"comment":"The height function is written as H(R) = R - m_{C+} ln R - r, but r is the compactified coordinate; since H should be a function of R only, this expression mixes the two coordinate systems. Please clarify, presumably by writing r = r(R).","section":"Section III.B, Eq. (25)"},{"comment":"The symbol '˜C□' appears instead of the tilde-C-minus variable (\\tilde{C}_-); please fix.","section":"Figure 1 caption"},{"comment":"In the text, 'as can bee seen' should read 'as can be seen'.","section":"Section IV.A"},{"comment":"The sentence 'We then computed the norm ... by use of the norm' repeats the word 'norm'; please rephrase.","section":"Section IV.A"},{"comment":"The norm notation uses Z generically; please state explicitly that the norm is summed over all evolved fields, or define Z as a vector of fields.","section":"Section IV.A, Eq. (41)"},{"comment":"The axis label ' at +' should read 'Ψ at I^{+}'; also specify the time intervals used for the quasinormal-mode fit and the tail fit.","section":"Section IV.B, Figure 4"},{"comment":"The y-axis label '1e 5+1' is confusing; it presumably means 1 + 10^{-5} times the plotted quantity. Please relabel.","section":"Section IV.B, Figure 3"}],"recommendation":"major_revision","confidential_remarks":"The paper is a reasonable incremental advance, but the authors themselves flag the main gap: the boundedness of \\hat{C}_+ is assumed and only informally tested. I would urge the editor to request the additional diagnostics described in the major comments. The auxiliary Mathematica notebooks are a useful strength, and the physical validation targets are well chosen."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"First thing you should know: this paper achieves something the authors' previous attempt couldn't — numerically stable, self-convergent evolutions at the n=2 compactification in the spherical DF-GHG system. The new ingredients are a specific constraint addition, a new reduction field for C+, and a second-order term in the C+ rescaling. The numerical results are credible: the convergence factors head to 2, the Bondi mass decays monotonically and settles, and the scalar at I+ shows the right QNM frequency and a t^-2 tail. That is a real technical step forward, and the physics checks are exactly the right ones.\n\nThe soft spot is the one the authors themselves point out: the whole n=2 setup depends on the second-order coefficient hat_C+ in the C+ rescaling staying bounded. They say the GBUF heuristics from [20] don't guarantee this, and the run is meant as a test. But the test is thin — one amplitude, one r_I, one integration time, no monitor of hat_C+ or the Jacobian factor. The self-convergence norm is computed over the whole domain and could easily miss growth localized near I+. So the claim \"n=2 works\" is, strictly speaking, \"n=2 works in this one configuration and under an unproven asymptotic assumption.\" That's not a fatal flaw — the authors flag it — but it should be front and center in any review.\n\nI also notice the paper doesn't give a well-posedness argument for the final reduced system with the new constraint addition, and the reproducibility is partial: Mathematica notebooks for the symbolics are available, but the evolution code itself isn't posted. That's not disqualifying, but it would help.\n\nThe citation pattern is honest — heavy use of the authors' own prior work on the GBUF model is natural here, since the paper is a direct continuation. The new physics checks against linear theory are external anchors, so the circularity concern is mild.\n\nOverall: this is a solid, honest technical contribution. The open question about hat_C+ is real and should be tested more aggressively — different amplitudes, larger r_I, longer runs, and a direct diagnostic. But the paper deserves a serious referee and likely publication after some strengthening. I'd take it to the reading group.","headline":"Real technical advance on n=2 hyperboloidal compactification with credible numerics, but the regularity claim rests on a single test of an unproven asymptotic coefficient.","tokens_in":14474,"tokens_out":2683,"would_cite":true,"duration_ms":29746,"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":"In spherical symmetry, the paper shows that the DF-GHG formulation of GR can produce self-convergent numerical evolutions that include future null infinity, with the strongest compactification parameter $n=2$, and that the waves extracted…","keywords":["numerical relativity","hyperboloidal compactification","future null infinity","generalized harmonic gauge","dual foliation","quasinormal modes","Price tail","constraint addition"],"falsifier":"Evolve the identical initial data at $n=2$ while tracking $\\hat C_+$ at $\\mathcal{I}^+$: if this coefficient grows without bound, or if the Jacobian denominator $R'(1-H'C_+)$ vanishes on the grid, the regularization is only apparent. Independently, if the $n=2$ results do not converge, under resolution increase, to the same extracted waveform obtained with $n<2$ (where formally singular terms have trivial limits), the claim that $n=2$ captures the correct continuum solution would be refuted.","tokens_in":13412,"feed_emoji":"📡","tokens_out":12110,"duration_ms":114658,"temperature":0.7,"pith_summary":"The paper asks whether the dual-foliation generalized harmonic gauge (DF-GHG) formulation of general relativity can be evolved numerically all the way to future null infinity ($\\mathcal{I}^+$) with the strongest allowed hyperboloidal compactification, the parameter value $n=2$ that makes the setup directly analogous to conformal compactification. It answers yes in spherical symmetry, provided three ingredients are chosen together: a specific pair of gauge source functions, a specific constraint addition that converts formally singular terms into regular ones, and reduction fields rescaled with the correct powers of the areal radius. With those choices the dangerous denominator $R'(1-H'C_+)$ is made regular by writing $C_+ = 1 + m_{C_+}/R + \\hat C_+/R^2$. The numerical evolutions are self-convergent at second order, and for small scalar-field perturbations of Schwarzschild the field extracted directly at $\\mathcal{I}^+$ shows quasinormal ringing followed by a $t^{-2}$ tail, with the Bondi mass decreasing monotonically and settling to a constant. A reader should care because this is evidence that the strongest compactification, previously out of reach in this formulation, can be used in practice, opening a direct comparison with conformal methods.","feed_headline":"Null infinity reached with the strongest compactification","feed_subtitle":"Black-hole runs now converge with the strongest compactification and reproduce ringdown and tail decay.","key_machinery":"The central mechanism is the hyperboloidal compactification $R(r)=r/(1-r^2/r_I^2)$, the $n=2$ member of the family $R(r)=r/\\Omega(r)^{1/(n-1)}$, combined with the height function $H(R)=R-m_{C_+}\\ln R - r$ that lifts Cauchy slices so that they intersect future null infinity. What carries the argument is the Jacobian denominator $R'(1-H'C_+)$: it becomes formally singular at $\\mathcal{I}^+$ unless $C_+$ is expanded to second order and its incoming null derivative is rescaled with a full power $R^2$. Equally load-bearing are the Good-Bad-Ugly-F decay rates, a model classification of wave equations by their falloff at null infinity, which dictate which fields need which rescaling, and the constraint addition that substitutes regular combinations for formally singular terms in the evolution equations for $\\Theta_-$ and $\\Delta_-$.","core_discovery":"The central discovery is that the $n=2$ compactified hyperboloidal evolution of the DF-GHG system becomes numerically regular when the formally singular terms are eliminated by a specific combination of gauge source functions, constraint addition, and reduction fields. The key novelty relative to the previous spherical work is the rescaling $\\Theta_+ \\equiv R^2 D_\\sigma C_+/\\kappa$ and its companion, together with the second-order expansion $C_+ = 1 + m_{C_+}/R + \\hat C_+/R^2$, which makes the Jacobian factor $R'(1-H'C_+)$ finite at $\\mathcal{I}^+$ instead of formally singular. With this setup the equations still contain a small number of formally singular terms, handled at the single boundary point by l'Hôpital's rule, and the rescaled fields are $O(1)$ at $\\mathcal{I}^+$. The paper reports self-convergence factors approaching 2 for both constraint-violating and constraint-satisfying perturbations of Schwarzschild, and extraction of the scalar field at $\\mathcal{I}^+$ reproduces the fundamental quasinormal mode and the $t^{-2}$ tail from linear theory.","pith_inferences":["If $\\hat C_+$ stays bounded for a wider family of data and gauges, the same regularization should transfer to 3D, because the argument relies on asymptotic decay rates rather than on spherical symmetry; the spherical reduction, however, hides angular derivative terms that could introduce new formally singular combinations.","The $n=2$ DF-GHG evolutions are now directly comparable to conformal compactification codes on identical initial data, so a head-to-head waveform or Bondi mass comparison would be a sharper test of both approaches than matching linear theory alone.","The paper excludes constraint damping terms because they conflict with the regularizing constraint addition; a smooth blending of damping near the strong-field region with the asymptotically regularizing terms is an obvious extension, and its absence may set the practical accuracy limit at late times.","Only the spherically symmetric scalar mode is tested, so recovering the gravitational quasinormal mode spectrum and the corresponding tail for nonspherical perturbations remains an open check that the 3D extension will need to address."],"forward_implications":["The strongest compactification parameter $n=2$, the only one for which the DF-GHG metric fields admit a conformal-type compactification, is usable in practice rather than only at the continuum level.","Radiation can be read off directly at $\\mathcal{I}^+$ instead of being extrapolated from finite radius: the extracted scalar field shows the fundamental quasinormal mode followed by a $t^{-2}$ tail.","The Bondi mass computed at $\\mathcal{I}^+$ is positive, monotonically decreasing, and settles to a constant slightly above the initial mass, consistent with the black hole accreting part of the scalar field.","Both constraint-violating and constraint-satisfying initial data evolve stably, with self-convergence factors approaching the expected value 2 for second-order finite differences.","The construction is stated to extend to full 3D, where the same formal singularities appear, and the authors identify that extension as the next step."],"supporting_citations":[{"why":"Establishes the spherical DF-GHG setup and the previous compactification parameter range that the present paper extends to n=2.","marker":"[23]"},{"why":"Supplies the Good-Bad-Ugly-F asymptotic decay rates and the gauge-driver construction that justify the field rescalings.","marker":"[20]"},{"why":"Shows constraint addition can force the improved R^{-2} decay of incoming null derivatives of C+ and epsilon that the regularization uses.","marker":"[25]"},{"why":"Introduces the dual-foliation height-function construction and the n=2 conformal compactification analogy.","marker":"[19]"},{"why":"Provides the scri-fixing gauge whose leading form the gauge source functions (35) parallel.","marker":"[16]"},{"why":"Determines when an extra power of R can be used in rescaling null derivatives, informing the choice of the R^2 rescaling for C+.","marker":"[27]"},{"why":"Gives the truncation-error-matching one-sided stencil used to evaluate the remaining formally singular terms at the outer gridpoint.","marker":"[28]"},{"why":"Provides the quasinormal mode frequency omega M = 0.11 + 0.10i used to fit the ringing segment.","marker":"[30]"},{"why":"Predicts the t^{-2} Price tail that the late-time extracted field is compared against.","marker":"[31]"}],"fun_headline_variants":["Ringdown and tail extracted directly at null infinity","Stronger compactification enables regular evolution at null infinity","Black hole perturbations: ringdown and tail extracted at null infinity","Direct null infinity extraction with stronger compactification","No more singular terms: compactified GR reaches null infinity"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The scheme rests on the assumption that the asymptotic decay rates derived for a model wave system carry over to the full nonlinear Einstein equations, and in particular that the second-order coefficient $\\hat C_+$ in the expansion of $C_+$ stays bounded during the whole evolution; the paper states that the model analysis does not guarantee this bound and that the present run is a test of it.","fun_headline_variants_meta":{"raw":{"variants":["Ringdown and tail extracted directly at null infinity","Stronger compactification enables regular evolution at null infinity","Black hole perturbations: ringdown and tail extracted at null infinity","Direct null infinity extraction with stronger compactification","No more singular terms: compactified GR reaches null infinity"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000962,"raw_usage":{"total_tokens":4096,"prompt_tokens":942,"completion_tokens":3154,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":558,"completion_tokens_details":{"reasoning_tokens":3078}},"tokens_in":558,"tokens_out":3154,"duration_ms":23814,"temperature":1.0,"reasoning_tokens":3078,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T11:10:03.665547+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Evolve the identical initial data at $n=2$ while tracking $\\hat C_+$ at $\\mathcal{I}^+$: if this coefficient grows without bound, or if the Jacobian denominator $R'(1-H'C_+)$ vanishes on the grid, the regularization is only apparent. Independently, if the $n=2$ results do not converge, under resolution increase, to the same extracted waveform obtained with $n<2$ (where formally singular terms have trivial limits), the claim that $n=2$ captures the correct continuum solution would be refuted.","supporting_citations":[{"cited_title":"Price, Nonspherical perturbations of relativistic grav- itational collapse","cited_arxiv_id":null,"evidence_quote":"Predicts the t^{-2} Price tail that the late-time extracted field is compared against."}],"review_version":1}