{"id":"24e91522-2e5e-4049-94fa-8b9d25600661","arxiv_id":"2608.07934","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"First frequency-domain m-mode calculation of Schwarzschild metric perturbations in Lorenz gauge, solving ten coupled elliptic PDEs and matching known energy fluxes to about four digits.","lead":"This paper computes gravitational wave metric perturbations of a Schwarzschild black hole with a new frequency-domain approach that keeps only the m quantum number, and validates the resulting energy fluxes to about four significant digits. The method is a stepping stone toward second-order self-force calculations for spinning (Kerr) black holes, which LISA-era waveform models will need.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Near-horizon and gauge selection is load-bearing, but no interior Lorenz-gauge residual or local metric check is reported; asymptotic flux agreement cannot certify the full 10-component retarded field.","rationale":"I largely agree with the reader's weakest_assumption but sharpen it: the near-horizon boundary conditions are the suspected mechanism, and the missing check is the interior Lorenz-gauge residual and local field comparison. The external flux agreement with the Black Hole Perturbation Toolkit is genuine independent support, and the Richardson convergence in Figure 7 strengthens the flux numbers; I am not claiming the computed fluxes are wrong. The concern targets the stronger statement that a full Lorenz-gauge metric perturbation has been obtained. A gauge-residual computation is cheap and decisive: if the residuals converge to zero, the conditional verdict can be upgraded; if they do not, the method is not validated for the stated second-order self-force application. The paper's self-identified limitations, including the skipped m = 0 mode, the unreleased Mathematica code, and the explicit progress-report framing, reinforce but do not drive this verdict. I therefore keep the reader's CONDITIONAL verdict unchanged, with the condition made explicit.","tokens_in":23356,"tokens_out":16446,"duration_ms":191355,"concrete_test":"Compute the four Lorenz-gauge residuals G_nu = g^alpha beta nabla_alpha hbar_beta nu from the converged numerical solution for r0 = 6M, m = 1 and m = 2, on the same (r*, theta) grid, evaluating both inside and outside the worldtube, and repeat at two resolutions separated by a factor of two (e.g., Delta r* = 0.2M and 0.1M). Report max|G_nu| normalized by max|psi| and its Richardson convergence factor. If the residual does not converge toward zero at the expected second-order rate and fall below roughly 1e-4 of the field amplitude, the solution is not demonstrably in Lorenz gauge and the central claim for second-order self-force is unsupported.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section V C assumes that four Lorenz-gauge boundary conditions plus generalized Bayliss-Turkel conditions on psi4 through psi9 at r*_min approximately -6M eliminate all 14 unphysical homogeneous solutions and select the unique retarded field. This assumption is the load-bearing step for the central claim. The paper's quantitative validation, however, is an integrated asymptotic energy flux (Table II) plus qualitative plots against the time-domain code of Ref. [55]; no interior diagnostic is reported. In particular, Eq. (2.2) is imposed explicitly only at the near-horizon boundary, and the fourth-order puncture of Section III yields an effective source in Eq. (3.4) whose divergence need not vanish exactly. The reduced field equations therefore admit solutions with nonzero gauge residual that can still produce the correct radiative flux, because the flux is dominated by gauge-invariant radiative content. For the stated purpose of feeding second-order self-force calculations, the full metric perturbation in the correct Lorenz gauge is required, so the present evidence is necessary but not sufficient. This is an omitted validation, not an internal inconsistency; the paper itself describes the work as partly a progress report.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript develops a frequency-domain, m-mode method for computing first-order Lorenz-gauge metric perturbations of Schwarzschild. The authors reduce the ten coupled perturbation equations to elliptic PDEs in (r*, theta), discretize with second-order finite differences, regularize the particle singularity with a fourth-order effective-source puncture, and solve the resulting sparse linear systems. After diagnosing near-horizon pathologies through an eigenvalue analysis of the homogeneous solutions, they introduce improved near-horizon boundary conditions consisting of the four Lorenz-gauge conditions plus higher-order Bayliss-Turkel-type conditions on six components. They report outgoing energy fluxes for circular equatorial orbits at r0 = 6M, 8M, and 10M and find agreement with independent Black Hole Perturbation Toolkit Teukolsky data to about four significant digits. The paper also presents the components of the retarded metric perturbation and compares them qualitatively with time-domain results from Ref. [55].","tokens_in":23530,"tokens_out":4514,"duration_ms":55917,"significance":"If the full metric perturbation is correct, this is the first frequency-domain m-mode calculation of Lorenz-gauge Schwarzschild metric perturbations and an important step toward second-order Kerr self-force calculations. The flux validation is external and model-independent: the Black Hole Perturbation Toolkit Teukolsky results are not fitted, and no free parameter is tuned. The paper also contains a useful near-horizon eigenvalue analysis (Table I and Appendix B), novel higher-order near-horizon radiation boundary conditions (Section V C), and a new puncture expansion for a point mass in Kerr spacetime (Section III). The authors are appropriately cautious, explicitly describing the work as partly a progress report. However, the quantitative evidence for the central claim is currently limited to integrated fluxes and qualitative field plots; a full interior validation of the ten-component Lorenz-gauge field is missing.","major_comments":[{"comment":"The quantitative evidence that the full ten-component field is the retarded Lorenz-gauge solution is incomplete. Equation (2.2) is imposed only as a boundary condition at r_min*, and the effective source in Eq. (3.4) need not be divergence-free because the puncture in Sec. III is truncated at fourth order. The main validation in Table II is an integrated asymptotic flux, which is dominated by radiative, largely gauge-invariant content, together with qualitative agreement with Ref. [55] in Figs. 4 and 5. For the stated purpose of feeding second-order self-force calculations, the full metric perturbation in the correct Lorenz gauge is required, so this is an omitted validation rather than an internal inconsistency. Please report interior residuals of both the field equations and the Lorenz condition, and show convergence of individual psi^m_mu_nu components, not only the flux.","section":"Sec. V C and Sec. VI B"},{"comment":"The boundary-condition construction is not demonstrated to select the unique retarded solution. Table I identifies six homogeneous solutions with the correct eigenvalue lambda = -i m Omega that satisfy the Lorenz gauge condition; imposing the four gauge conditions and the Bayliss-Turkel-type conditions on psi_4 through psi_9 gives ten conditions, but the degeneracy means uniqueness does not follow automatically from counting conditions. The text asserts that Eq. (5.9c) at r_min* ~ -6M eliminates the problematic solutions, but no numerical experiment isolates this assumption. Please demonstrate uniqueness by, e.g., varying r_min* and the order of the boundary conditions, or by checking that the exponentially growing homogeneous modes are absent from the interior solution.","section":"Sec. V C"},{"comment":"The convergence evidence is limited to the m=2 flux (Fig. 7) and to total fluxes at three radii. Since Sec. V uses the sensitive m=1 mode as the primary test of the corrected near-horizon conditions, the paper should show m=1 (and at least one higher-m) flux convergence and per-mode Richardson-error estimates. Without this, the claimed 'approximately 4 significant digits' is not supported for all modes that contribute to the total flux.","section":"Sec. VI B and Table II"}],"minor_comments":[{"comment":"The word 'certiain' in the sentence 'One way of calculating h_mu_nu that has certiain advantages' should be 'certain'.","section":"Sec. II A"},{"comment":"In the sentence following Eq. (4.11), 'All three above conditions' should read 'All three of the above conditions'.","section":"Sec. IV B"},{"comment":"There is a duplicated word in 'satisfies the the larger boundary condition'; it should be 'satisfies the larger boundary condition'.","section":"Sec. V C"},{"comment":"The phrase 'A through explanation of the advantages' should read 'A thorough explanation of the advantages'.","section":"Sec. I"}],"recommendation":"major_revision","confidential_remarks":"This is a credible methods paper with a clear external flux validation and a transparent description of its limitations. The main barrier to acceptance is the missing interior validation of the Lorenz-gauge field; adding residual checks and per-mode convergence data should address it without changing the paper's scope. The paper also relies heavily on the authors' own prior work, which is appropriate for a methods-focused contribution."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Bottom line: this is the first frequency-domain m-mode Lorenz-gauge metric perturbation calculation on Schwarzschild, and the central numerical result—outgoing energy flux—agrees with independent Teukolsky data from the BHPToolkit to about four significant digits. That is real, externally validated progress, not a fitting exercise. The paper deserves a serious referee.\n\nWhat's genuinely new: the m-mode frequency-domain elliptic PDE approach for gravitational perturbations in Lorenz gauge. Prior m-mode work was scalar or time-domain. The near-horizon analysis is also worth reading: the paper identifies 20 homogeneous solutions, shows only 6 with Im(λ)<0 and Lorenz-gauge-compatible eigenvalues survive, and builds generalized Bayliss-Turkel boundary conditions plus the four gauge conditions at r*_min ≈ -6M. That strategy fixes the spurious exponential decay that plagued naive Sommerfeld conditions, and the m=1 mode—the sensitive one—now matches the time-domain code within its noise.\n\nThe soft spots are real but not fatal. The quantitative validation is an integrated asymptotic flux plus qualitative comparison with Dolan-Barack time-domain data. There is no reported check that the full 10-component residual field satisfies the Lorenz gauge condition inside the domain, and no local interior comparison of the metric perturbation itself. The flux is dominated by radiative, gauge-invariant content, so agreement at the 4-digit level does not by itself certify that the full retarded field is correct in the Lorenz gauge—which is exactly what a second-order self-force pipeline needs. The paper is honest about this: it calls itself partly a progress report. I'd also note the code and data are not released, the m=0 non-radiative mode is skipped, and there are no explicit error bars beyond the significant-digit claims. None of these are load-bearing flaws; they are missing evidence.\n\nThe citation pattern and related-work discussion look fair. The claim that this is a first is credible against the cited literature. The eigenvalue table and appendix are reproducible in principle, and the puncture expansion is described in enough detail to be useful.\n\nWho is this for? People working on self-force, Lorenz gauge metric perturbations, or numerical PDE methods in black hole perturbation theory. It's not a general-audience paper, but for that subfield it's a meaningful step.\n\nRecommendation: send it out. A good referee should push for an interior Lorenz-gauge residual diagnostic or local comparison, and for code/data release or at least a public version, but the core methodological advance and external flux validation justify referee time.","headline":"First frequency-domain m-mode Lorenz-gauge Schwarzschild solver with externally validated fluxes to ~4 digits; near-horizon BCs are load-bearing and the missing interior gauge-residual check is the main gap.","tokens_in":24117,"tokens_out":1888,"would_cite":true,"duration_ms":22000,"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":"The paper claims the first frequency-domain m-mode Schwarzschild metric perturbations, computed by solving ten coupled elliptic PDEs in the Lorenz gauge, with outgoing energy fluxes accurate to about four significant digits.","keywords":["black hole perturbation theory","Lorenz gauge","m-mode decomposition","elliptic PDE solver","effective source method","gravitational self-force","Schwarzschild metric perturbations","near-horizon boundary conditions"],"falsifier":"Compute the Lorenz-gauge residual of the numerical solution everywhere, especially near $r_*\\approx -6M$, at increasing resolution: if the boundary conditions select the true retarded solution, the residual should fall at the expected second-order rate as the grid is refined, while a contaminated solution would show a residual that stagnates or grows near the inner boundary. A second check is to move the inner boundary much closer to the horizon and verify that the flux and full field are unchanged; if the improved conditions only work at $-6M$, the selection is not robust.","tokens_in":23082,"feed_emoji":"🕳️","tokens_out":8313,"duration_ms":82263,"temperature":0.7,"pith_summary":"The paper tries to establish that first-order gravitational perturbations of a Schwarzschild black hole can be computed directly in the Lorenz gauge by an m-mode frequency-domain scheme, in which each azimuthal mode obeys a system of ten coupled elliptic partial differential equations in the tortoise radius and polar angle. This matters because the eventual goal is second-order perturbations of a Kerr black hole, where the field equations do not separate in the usual tensor-harmonic basis and mode coupling becomes severe; an m-mode elliptic approach sidesteps both problems. The authors report the first such Schwarzschild calculation, regularized by an effective source with a new puncture expansion of the singular field of a point mass in Kerr spacetime. They show that the resulting outgoing gravitational-wave energy fluxes for circular equatorial orbits agree with established independent values to roughly four significant digits, which they argue is enough for future second-order self-force calculations. A central part of the story is the diagnosis and cure of near-horizon boundary contamination, which initially produced unphysical decaying modes.","feed_headline":"Ten coupled PDEs solve Schwarzschild metric modes to four digits","feed_subtitle":"New elliptic m-mode solver computes the retarded metric perturbation and gravitational-wave energy flux accurately enough for second-order…","key_machinery":"The object that carries the argument is the ten-component vector $\\vec{\\psi}_m$ of rescaled Lorenz-gauge metric perturbation modes, satisfying an elliptic system of the form $(\\partial^2/\\partial r_*^2 + \\Delta/r^4\\, \\partial^2/\\partial\\theta^2 + A_m\\,\\partial/\\partial r_* + B_m\\,\\partial/\\partial\\theta + C_m)\\vec{\\psi}_m = \\vec{S}_m$. Three mechanisms make it work. First, the effective source method replaces the point-particle delta source with a smooth source on a worldtube, so the elliptic PDE can be solved for a residual field; the puncture expansion of the singular field, new here and written for Kerr spacetime, supplies the subtraction. Second, a second-order finite-difference stencil converts the PDE system into a sparse linear algebra problem for each m-mode. Third, the decisive ingredient is the near-horizon boundary treatment: an eigenvalue analysis of the homogeneous system near $r=2M$ shows that only six of twenty homogeneous solutions have the correct ingoing behavior $\\lambda=-im\\Omega$ and satisfy the Lorenz gauge, and the paper's boundary conditions combine the four gauge conditions with higher-order frequency-domain radiation conditions, evaluated at $r_*\\approx -6M$, to exclude the other fourteen.","core_discovery":"The central claim, stated by the authors in Section VII, is that they have 'calculated Schwarzschild metric perturbations via m-modes in the frequency domain for the first time,' by solving the ten coupled elliptic PDEs of the Lorenz gauge on a two-dimensional grid in r* and theta. The loading idea is that each m-mode of the trace-reversed metric perturbation is represented by a ten-component vector psi_m with carefully chosen prefactors, so that the field equations become a single elliptic system whose source vanishes outside a worldtube around the particle. Regularization is handled with the effective source method, and the paper introduces a puncture expansion written for a point mass in Kerr spacetime, evaluated here in the Schwarzschild limit. The main numerical obstacle was spurious near-horizon behavior: a naive ingoing Sommerfeld condition permits unphysical homogeneous solutions that decay exponentially toward the horizon, contaminating the global solution. The authors claim that imposing the four Lorenz gauge conditions at the inner boundary together with higher-order radiation boundary conditions derived from a 20-eigenvalue near-horizon analysis selects the physical retarded solution, and that the outgoing energy flux computed from it converges to approximately four significant digits.","pith_inferences":["If the boundary-condition selection is as robust as the flux agreement suggests, the m-mode elliptic approach could replace metric reconstruction for Kerr self-force calculations, since it produces the metric perturbation in the Lorenz gauge directly rather than reconstructing it from curvature scalars.","The failure of the naive ingoing Sommerfeld condition for gravitational perturbations, where the zero field also satisfies the condition, suggests that frequency-domain Lorenz-gauge PDE solvers generally need an explicit gauge-condition boundary treatment; scalar-field solvers may not be a reliable guide.","A natural next test is to apply the same near-horizon strategy to the omitted m=0 sector and to eccentric orbits; if the six-eigenvalue selection remains sufficient there, it would strengthen the claim that the method transfers to Kerr.","The Kerr puncture expansion derived here is presented in a form usable by time-domain m-mode codes as well, so it could serve other effective-source implementations beyond the frequency-domain solver presented."],"forward_implications":["The full retarded Lorenz-gauge metric perturbation is obtained directly, so downstream self-force quantities do not require a separate metric reconstruction step.","The flux accuracy of about four significant digits is sufficient for second-order dissipative self-force calculations, the stated target application.","Because all ingredients are written in m-mode form and the puncture is given for Kerr spacetime, the same elliptic strategy can in principle accept a second-order Kerr source once the source modes are built.","The near-horizon eigenvalue criterion identifies exactly which homogeneous modes are physical at the horizon, giving a transferable test for any future m-mode solver in the Lorenz gauge."],"supporting_citations":[{"why":"Supplies the second-order finite-difference elliptic-PDE method and the outer radiation boundary conditions that this work extends from the scalar field to Lorenz-gauge metric perturbations.","marker":"[58]"},{"why":"Provides the time-domain m-mode gravitational perturbation results in the Lorenz gauge used for comparison and motivates the prefactor scalings in Eq. (2.9).","marker":"[55]"},{"why":"Establishes the m-mode approach and its advantages for second-order Kerr sources, framing why the frequency-domain elliptic strategy is worth building.","marker":"[45]"},{"why":"Supplies the scalar puncture decomposition into m-modes via complete elliptic integrals that the gravitational puncture and effective source inherit.","marker":"[56]"},{"why":"Provides the Q-R effective-source scheme used to construct the regularized source and the worldtube treatment.","marker":"[61]"},{"why":"Provides the high-order singular-field expansion in Schwarzschild that underlies the puncture form.","marker":"[62]"},{"why":"Provides the corresponding high-order singular-field expansion in Kerr spacetime used for the new Kerr puncture.","marker":"[63]"},{"why":"Supplies the radiation boundary conditions whose higher-order frequency-domain forms are generalized to the near-horizon regime.","marker":"[65]"},{"why":"Provides the independent flux comparison values used to validate the total outgoing energy flux in Table II.","marker":"[67]"},{"why":"Provides the curvature-perturbation flux results used as the precision comparison in Table II.","marker":"[68]"}],"fun_headline_variants":["First Schwarzschild m-mode solution via Lorenz gauge elliptic equations","Elliptic m-mode solver computes Schwarzschild metric to 4 digits","Ten coupled elliptic equations solve Schwarzschild perturbations","New elliptic solver achieves four-digit accuracy for Schwarzschild m-modes"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The whole calculation stands on the assumption that the special conditions imposed at the inner edge of the grid (the four gauge conditions plus a set of improved wave-outgoing conditions at about six gravitational radii from the horizon) actually filter out all unphysical solutions and leave the one true retarded wave; if a wrong mode slips through, the computed energy fluxes could match the independent values by accident.","fun_headline_variants_meta":{"raw":{"variants":["First Schwarzschild m-mode solution via Lorenz gauge elliptic equations","Elliptic m-mode solver computes Schwarzschild metric to 4 digits","Ten coupled elliptic equations solve Schwarzschild perturbations","New elliptic solver achieves four-digit accuracy for Schwarzschild m-modes"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000732,"raw_usage":{"total_tokens":3308,"prompt_tokens":1009,"completion_tokens":2299,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":625,"completion_tokens_details":{"reasoning_tokens":2231}},"tokens_in":625,"tokens_out":2299,"duration_ms":16698,"temperature":1.0,"reasoning_tokens":2231,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T00:39:07.803734+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the Lorenz-gauge residual of the numerical solution everywhere, especially near $r_*\\approx -6M$, at increasing resolution: if the boundary conditions select the true retarded solution, the residual should fall at the expected second-order rate as the grid is refined, while a contaminated solution would show a residual that stagnates or grows near the inner boundary. A second check is to move the inner boundary much closer to the horizon and verify that the flux and full field are unchanged; if the improved conditions only work at $-6M$, the selection is not robust.","supporting_citations":[{"cited_title":"Wardell, N","cited_arxiv_id":null,"evidence_quote":"Provides the curvature-perturbation flux results used as the precision comparison in Table II."}],"review_version":1}