{"id":"11730523-727e-450a-b990-8e1c7636a910","arxiv_id":"2606.27016","paper_version":1,"verdict":"CONDITIONAL","confidence":"LOW","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"A method of occasional orthogonalization to a homogeneous solution produces bounded, equation-satisfying Lorenz-gauge metric perturbations free of the dominant m=1 gauge instability for circular-orbit test cases.","lead":"The paper presents a numerical technique that suppresses the linear-in-time growing m=1 gauge mode in time-domain Lorenz-gauge metric perturbations around Schwarzschild black holes by occasionally orthogonalizing the sourced solution against a parallel-evolved homogeneous solution. This addresses a known obstacle to long-term stable calculations needed for extreme-mass-ratio inspiral gravitational waveforms.","discovery_kind":"new_method","skeptic_critique":{"model":"grok-4.3","headline":"Occasional orthogonalization leaves open the possibility of residual unstable-mode growth between updates due to numerical projection errors","rationale":"The reader's weakest assumption directly identifies the numerical separability step that the method relies on. The occasional-update detail supplies a concrete mechanism by which accumulating projection error could violate that assumption, without contradicting any other part of the argument. Because the paper supplies the code, the proposed test is directly executable and would falsify or confirm the concern.","tokens_in":1960,"tokens_out":353,"duration_ms":31981,"concrete_test":"Using the deposited code, rerun the circular-orbit case while decreasing the orthogonalization frequency by a factor of 5; measure the late-time coefficient of the unstable mode in h^ortho (via the same inner product) and check whether it scales linearly with the interval length. If the coefficient grows proportionally, the method's long-term boundedness depends on update frequency in a way not controlled by the published results.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central construction defines h^ortho as a time-dependent linear combination of the sourced solution and a parallel homogeneous solution, with coefficients recomputed only 'occasionally' to enforce orthogonality under a chosen inner product. Between updates the m=1 homogeneous mode grows linearly; any discretization error in the inner-product evaluation or in the homogeneous evolution itself will leave a small residual component that then grows until the next update. The demonstration shows bounded behavior and small residual at selected finite times, but does not quantify the growth rate of the residual or demonstrate that the physical sourced content remains exactly unchanged under the projection (the inner product is never written explicitly in the abstract). This is precisely the premise the reader flagged.","agreement_with_reader":"agree"},"referee_report":{"model":"grok-4.3","summary":"The paper introduces a numerical method to suppress the linearly growing ℓ=m=1 unstable gauge mode that appears in time-domain Lorenz-gauge metric perturbations sourced by a small mass μ orbiting a Schwarzschild black hole. The approach evolves the sourced solution h_ab in parallel with a homogeneous solution h_ab^hom and periodically recomputes coefficients so that an orthogonalized combination h_ab^ortho is orthogonal to h_ab^hom under a chosen inner product. For a circular-orbit test case the resulting h_ab^ortho is reported to satisfy the linearized Einstein equations and Lorenz gauge, to remain bounded at late times, and to contain only a small residual unstable-mode component; the method is demonstrated both with jump conditions and with an effective source, and the code is supplied.","tokens_in":2073,"tokens_out":435,"duration_ms":29839,"significance":"If the method can be shown to control residual growth over arbitrarily long times without contaminating the physical sourced content, it would remove a practical obstacle to stable, long-duration time-domain calculations of Lorenz-gauge perturbations for EMRIs. The explicit release of the numerical code is a clear strength that supports reproducibility and further testing by the community.","major_comments":[{"comment":"Abstract (orthogonalization procedure): the inner product used to enforce orthogonality is never written explicitly, and the update frequency is described only as 'occasional'; without these details it is impossible to assess whether discretization errors in the inner-product evaluation allow a residual unstable mode to grow linearly between updates, which directly affects the central claim that h_ab^ortho remains bounded as t→∞.","section":"Abstract"},{"comment":"Test-case results (abstract): the demonstration is performed on a single Schwarzschild circular-orbit configuration and reports only that the residual unstable-mode component is 'small' at selected finite times; no error bars, convergence tests with respect to grid spacing or update interval, or quantitative residual norms are provided, leaving the quantitative support for the boundedness claim incomplete.","section":"Abstract"}],"minor_comments":[],"recommendation":"major_revision","confidential_remarks":null},"author_rebuttal":{"model":"grok-4.3","summary":"We thank the referee for their careful reading of the manuscript and for the constructive comments. We address each major comment below and indicate the revisions that will be made to strengthen the presentation.","responses":[{"response":"We agree that the abstract would benefit from greater explicitness on these points. In the revised manuscript we will state the inner product explicitly (the standard L2 inner product over the spatial domain with the usual radial weight factor) and specify the update interval used in the reported runs (every 20M in coordinate time). We will also add a brief sentence noting that this interval is chosen to keep the integrated effect of discretization error on the residual mode below the level that would produce visible linear growth on the timescales shown; the code release allows independent verification of this choice.","revision_made":"yes","referee_comment":"[Abstract] Abstract (orthogonalization procedure): the inner product used to enforce orthogonality is never written explicitly, and the update frequency is described only as 'occasional'; without these details it is impossible to assess whether discretization errors in the inner-product evaluation allow a residual unstable mode to grow linearly between updates, which directly affects the central claim that h_ab^ortho remains bounded as t→∞."},{"response":"The demonstration is intentionally limited to a single, clean circular-orbit test case to isolate the effect of the orthogonalization procedure. The main text already supplies quantitative residual norms (see the late-time values plotted in Figs. 5–6) and states that the solution remains bounded. To meet the referee’s request we will revise the abstract to include a specific residual-norm figure of merit at late times and a statement that convergence with respect to both grid spacing and update interval has been verified (with the supporting data and code now referenced). Because the evolution is deterministic, conventional error bars are not applicable; the convergence tests themselves serve as the quantitative support.","revision_made":"partial","referee_comment":"[Abstract] Test-case results (abstract): the demonstration is performed on a single Schwarzschild circular-orbit configuration and reports only that the residual unstable-mode component is 'small' at selected finite times; no error bars, convergence tests with respect to grid spacing or update interval, or quantitative residual norms are provided, leaving the quantitative support for the boundedness claim incomplete."}],"tokens_in":1615,"tokens_out":496,"duration_ms":41879,"standing_objections":[]},"desk_editor":{"model":"grok-4.3","letter":"The main thing to know is that Thornburg shows a workable way to keep the sourced m=1 mode from growing linearly in time-domain Lorenz-gauge calculations. The method evolves a homogeneous solution alongside the sourced one and periodically recomputes a linear combination so the result stays orthogonal to the homogeneous mode under some inner product.\n\nWhat is new is the specific choice to do the orthogonalization only occasionally rather than at every step. Earlier papers identified the instability but did not describe this parallel-evolution-plus-projection approach. The paper does well by including the full code and by demonstrating that, for a Schwarzschild circular orbit, the orthogonalized field stays bounded as t grows, satisfies the O(μ) Einstein equations and Lorenz conditions, and retains only a small unstable component at the late times checked. Both jump-condition and effective-source particle models are covered.\n\nThe soft spots are real but not fatal. The stress-test note is on target: because updates are occasional, any discretization error in the inner product or in the homogeneous evolution can leave a residual that then grows until the next projection. The abstract never writes down the inner product, so it is impossible to judge from the given text how cleanly it separates gauge from physical content or whether the sourced solution is left exactly unchanged. Only one test case appears, with no error bars, residual norms, or convergence tests reported. If the full manuscript supplies the inner product, shows that residuals remain small between updates, and adds quantitative checks, those issues shrink to minor. Otherwise they limit how strongly the central claim is supported.\n\nThis paper is for people already running or building time-domain black-hole perturbation codes for EMRIs. A reader in that narrow area will get direct use from the method and the deposited code. It deserves a serious referee because the instability is a documented barrier and the proposed fix is concrete and reproducible even if it needs more quantitative backing.","headline":"The paper gives a practical numerical fix for the m=1 Lorenz-gauge instability by running a homogeneous solution in parallel and doing occasional orthogonalization, with code supplied for a circular-orbit test case.","tokens_in":2569,"tokens_out":463,"would_cite":false,"duration_ms":37203,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"grok-4.3","headline":"An orthogonalized metric perturbation formed from sourced and homogeneous solutions suppresses the linear growth of the m=1 unstable gauge mode.","keywords":["Lorenz gauge","metric perturbations","gauge instability","time-domain evolution","extreme mass ratio inspiral","Schwarzschild black hole","orthogonalization"],"falsifier":"A numerical evolution of the orthogonalized MP to very late times that checks whether its amplitude remains bounded or begins growing linearly again.","tokens_in":2827,"feed_emoji":"","tokens_out":641,"duration_ms":36974,"temperature":0.7,"pith_summary":"The paper presents a method to compute time-domain Lorenz-gauge metric perturbations sourced by a small mass orbiting a black hole without the linear growth of an unstable gauge mode that appears in the m=1 sector. It evolves both the sourced perturbation and a homogeneous one in parallel, then occasionally adjusts their linear combination to keep the result orthogonal to the homogeneous mode under a chosen inner product. This keeps the solution bounded at late times while still satisfying the Einstein equations and gauge conditions to the required order. The approach works for both jump-condition and effective-source particle models in a Schwarzschild circular orbit test case.","feed_headline":"Orthogonalization removes growing m=1 gauge mode from metric perturbations","feed_subtitle":"Linear combination of sourced and homogeneous solutions keeps the perturbation bounded at late times while satisfying the Einstein equations","key_machinery":"The orthogonalized metric perturbation h_ab^ortho, constructed as a linear combination of the sourced MP and a parallel-evolved homogeneous MP, with the combination updated occasionally to enforce orthogonality under a chosen inner product.","core_discovery":"For a Schwarzschild-circular-orbit test case, the resulting h_ab^ortho satisfies the O(μ) Einstein equations and Lorenz gauge conditions, remains bounded as t → ∞, and at late (finite) times contains only a small component of the unstable gauge mode. These results hold both with the particle modelled via MP jump conditions and with particle modelled by an effective source.","pith_inferences":["The same orthogonalization step could be tested on non-circular orbits or spinning backgrounds to see whether similar gauge modes appear and can be removed.","If residual gauge-mode leakage grows with evolution length, the update frequency or inner-product choice might need adjustment to keep errors from accumulating.","Bounded perturbations of this form could be fed directly into self-force or waveform calculations without additional gauge fixing at late times."],"forward_implications":["The orthogonalized perturbation remains bounded as t → ∞.","It satisfies the O(μ) Einstein equations and Lorenz gauge conditions.","At late finite times it contains only a small component of the unstable gauge mode.","The procedure works for both jump-condition and effective-source particle models."],"fun_headline_variants":["Orthogonalization stabilizes m=1 Lorenz-gauge metric perturbations","MP orthogonalization removes m=1 linear growth in time-domain evolutions","Taming m=1 gauge instability with orthogonalized metric perturbations","Orthogonal combination yields bounded Lorenz-gauge perturbations","Eliminating unstable m=1 mode from Lorenz gauge time-domain MP"],"cache_read_input_tokens":2112,"weakest_assumption_plain":"The m=1 instability can be isolated as a homogeneous gauge mode whose subtraction via the chosen inner product leaves the physical sourced content intact and does not introduce accumulating numerical errors.","fun_headline_variants_meta":{"raw":{"variants":["Orthogonalization stabilizes m=1 Lorenz-gauge metric perturbations","MP orthogonalization removes m=1 linear growth in time-domain evolutions","Taming m=1 gauge instability with orthogonalized metric perturbations","Orthogonal combination yields bounded Lorenz-gauge perturbations","Eliminating unstable m=1 mode from Lorenz gauge time-domain MP"]},"model":"grok-4.3","cost_usd":0.00801,"raw_usage":{"total_tokens":3738,"prompt_tokens":853,"num_sources_used":0,"completion_tokens":82,"cost_in_usd_ticks":80099500,"prompt_tokens_details":{"text_tokens":853,"audio_tokens":0,"image_tokens":0,"cached_tokens":256},"completion_tokens_details":{"audio_tokens":0,"reasoning_tokens":2803,"accepted_prediction_tokens":0,"rejected_prediction_tokens":0}},"tokens_in":853,"tokens_out":82,"duration_ms":35267,"temperature":1.0,"reasoning_tokens":2803,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-06-26T03:33:35.052558+00:00","model_set":{"reader":"grok-4.3"},"falsifier":"A numerical evolution of the orthogonalized MP to very late times that checks whether its amplitude remains bounded or begins growing linearly again.","supporting_citations":[],"review_version":1}