{"id":"385e9b5c-c240-4382-817f-00f530dee28f","arxiv_id":"2508.20545","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"The displacement memory formula for Generalized Proca gravity is derived for a massive Lorentz-invariant branch and a massless Lorentz-violating branch, with the dispersive branch requiring a frequency-integrated treatment.","lead":"This paper derives the gravitational memory effect, a permanent shift in the distance between test masses after a gravitational wave passes, for Generalized Proca gravity, a broad modified-gravity theory with an extra vector field. It finds two families of memory signals, one with massive dispersive waves and one with Lorentz-violating speed-varying waves, and it outlines how future memory observations could constrain Lorentz violation.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Equation (84) integrates over dispersive-mode speeds without a spectral Jacobian or frequency-resolved flux, so the total case-(a) memory formula is not actually computable as written; the per-mode result (83) may be sound, but the headline 'full memory formula' is underdetermined.","rationale":"The reader's weakest assumption identifies exactly the same gap: Eq. (84) converts a sum over dispersive frequencies into an integral over group velocities without specifying the spectral measure or Jacobian. I agree that this is the most load-bearing weakness of the central claim, because the paper advertises Eq. (84) as the total tensor memory in the massive Lorentz-invariant branch, and that expression cannot be evaluated from the information given. I do not see a more fundamental objection. The per-mode derivation in Eqs. (69), (70), and (83) appears internally consistent: the mass terms drop out of the averaged stress tensor after imposing the Klein-Gordon equation, and the group-velocity retarded-time choice is physically motivated. The Lorentz-violating case (b) results in Eqs. (86)-(88) rest on the earlier Einstein-Æther analysis and are not obviously problematic beyond the acknowledged breakdown of perturbation theory at the critical direction, which the authors state explicitly. The observational constraint argument in Eq. (89) is qualitative but consistent with the existing literature. No code, data, or machine-checked proof accompanies the paper, so the lengthy algebra was not independently verified; however, the main structural issue remains the under-specified spectral integration in Eq. (84). Because the reader already issued a CONDITIONAL verdict on this basis, my stress-test does not move the verdict.","tokens_in":29430,"tokens_out":8300,"duration_ms":86587,"concrete_test":"Re-derive Eq. (84) from Eq. (83) with an explicit frequency-resolved flux Fχ(ω;u',Ω'), then change variables using βχ=sqrt(1−m²/ω²). For a concrete source such as a burst with Fχ(ω;u',Ω')=A(ω)δ(u'−u₀) and A(ω)=ω^{−n}, compute the memory both by the correct frequency integral and by Eq. (84). If the two disagree by the omitted Jacobian m²/(ω³βχ) or by a needed spectral weight, Eq. (84) is not self-contained and should be restated either with the explicit measure ∫dω or restricted to the per-mode result (83).","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim includes Eq. (84) as the total tensor displacement memory in the massive Lorentz-invariant case (a). The derivation preceding it gives, in Eq. (83), the memory sourced by a single monochromatic plane-wave mode with frequency ωχ, whose flux Fχ(u',Ω') is implicitly frequency-dependent. To obtain the full memory one must sum over frequencies: Δh ~ ∫ dω Fχ(ω;u',Ω') βχ(ω) [Vχ n' n']^TT. Changing variables to βχ=sqrt(1−m²/ω²) requires the Jacobian dω/dβ = ω³β/m² (equivalently dβ/dω = m²/(ω³β)). Equation (84) instead writes ∫₀¹ dβχ Fχ(u',Ω') with the same Fχ used for a single mode, supplying no spectral density, no Jacobian, and no statement of what Fχ means as a function of βχ. This is not a cosmetic omission: different radiation spectra produce different total memories, and the expression as written is dimensionally and operationally ambiguous. The paper itself acknowledges in Sec. IV that the full memory depends on an explicit integration over the frequency content, which confirms that Eq. (84) is incomplete rather than a fully determined prediction. The per-mode formula and the Lorentz-violating case (b) results are less affected by this issue, but the case-(a) total-memory claim is the one that needs repair before the headline result can be used.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents a computation of gravitational displacement memory in Generalized Proca (GP) gravity. Building on the authors' earlier unified framework for memory in scalar-tensor and Einstein-Æther theories, the paper constructs a gauge-invariant second-order action on an asymptotically flat background, identifies two nontrivial background branches: (a) a Lorentz-invariant massive branch with dispersive vector and scalar modes, and (b) a Lorentz-violating massless branch with non-dispersive modes. It derives the effective stress-energy tensor via the second-variation method and obtains per-mode memory formulas, giving Eq. (84) for the total tensor memory in case (a) and Eqs. (86)-(88) in case (b), where a superluminal source sector produces an unbounded memory enhancement at a critical angle. The paper also discusses polarization content and observational constraints from GW170817 and Cherenkov radiation.","tokens_in":29820,"tokens_out":16595,"duration_ms":147539,"significance":"This work is the first computation of gravitational displacement memory in GP gravity and a valuable extension of the Isaacson/second-variation framework to dispersive massive modes. The gauge-invariant SVT decomposition and the identification of the two distinct background branches are clearly presented, and the per-mode memory formulas and the case-(b) results closely follow the structure of the well-established Einstein-Æther analysis. The distinction between phase and group velocity in the memory integral is an important conceptual point. However, the headline total memory formula for the Lorentz-invariant branch, Eq. (84), is underdetermined as written, because the conversion of the sum over dispersive modes into an integral over group velocities lacks a spectral measure and Jacobian. This is a fixable but load-bearing gap in the paper's central claim.","major_comments":[{"comment":"The total case-(a) memory formula is not well-defined as written. The per-mode result (83) is derived for a monochromatic mode with fixed frequency ωχ and group velocity βχ = sqrt(1 - m²/ωχ²). Summing over frequencies requires an integral over ωχ (or over βχ with the Jacobian dωχ/dβχ = ωχ³βχ/m²) and a frequency-resolved flux dFχ/dωχ. Eq. (84) instead writes ∫₀¹ dβχ Fχ(u',Ω') with the same Fχ defined in Eq. (70) for a single mode, giving no spectral density, no Jacobian, and no statement of the frequency dependence of the flux. The expression is therefore dimensionally and operationally ambiguous, and the total memory depends on the radiation spectrum. The paper's own Sec. IV acknowledges that the full memory formula requires an explicit integration over the frequency content, confirming that Eq. (84) is incomplete. The authors should replace Eq. (84) with an explicit spectral integral.","section":"Section III C, Eq. (84)"},{"comment":"The retarded time u' used in the integrand of Eq. (84) is the same for all frequencies and for both χ = S and χ = V. For dispersive modes, the appropriate retarded time is u'_χ(ωχ) = t' - r'/βχ(ωχ), which depends on the frequency. A correct frequency-resolved formulation must evaluate Fχ on the frequency-dependent retarded time and integrate over ωχ (or βχ) with the associated Jacobian; otherwise modes of different frequencies are mis-weighted. This issue is directly tied to the previous comment and should be resolved together.","section":"Section III C, Eq. (84)"},{"comment":"The assertion that the memory becomes unbounded at the critical direction n'·n = βT/βψ is plausible but not quantitatively demonstrated. Unlike the luminal case, where the TT projection removes the divergence at n' = n, the critical direction here is at a finite angle from the observer direction, so the TT projection of n'_i n'_j is nonzero and the angular integral appears to diverge logarithmically. The authors should provide the explicit behavior of the integrand near the critical direction (or cite the analogous calculation in Ref. [41]) and discuss the cutoff provided by finite source size or wave-packet width, since this underpins the proposed observational constraints.","section":"Section III C, case (b), Eqs. (86)-(88)"}],"minor_comments":[{"comment":"The sentence 'This can effectively be obtain through an integration over all possible emission speeds' contains a typo; 'obtain' should be 'obtained'.","section":"Section III C, text after Eq. (83)"},{"comment":"The symbol Vχ in the integrand is defined in Eq. (80) as an angular factor, while Vψ elsewhere denotes phase velocities (e.g., Eq. (23)); this dual use is confusing, and a distinct symbol such as \\mathcal{V}_χ would improve clarity.","section":"Section III C, Eq. (84)"},{"comment":"The constraint on the tensor-mode speed is written as β_T^(b) = 1 + O(10^{-15}); since the GW170817 bound is on |c_T - 1|, it would be more precise to state the deviation can be either sign or to use an absolute value.","section":"Section III C, Eq. (89)"},{"comment":"The technical assumptions 1 - σ\\bar A² ≠ 0 and 3\\bar A² \\bar G_{3,X}² + 4\\bar G_{2,XX}\\bar G_4(1 - σ\\bar A²) ≠ 0 are introduced without physical motivation; a brief remark that these conditions ensure the corresponding modes are dynamical would help the reader.","section":"Section II B 2, Eqs. (30)-(31)"},{"comment":"The square-root factors in the vector polarization amplitudes assume the combination under the square root is positive; the authors should state the sign conditions required for real, physical polarizations.","section":"Section II C, Eqs. (47c)-(47d)"}],"recommendation":"major_revision","confidential_remarks":"The manuscript extends the authors' established framework to GP gravity and is generally careful. The main technical gap is the underdetermined spectral integral in Eq. (84), which directly affects the headline case-(a) memory formula. The case-(b) analysis is solid, but the divergence at the critical angle deserves a quantitative statement. The paper is suitable for this journal after these points are addressed."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"This is a serious paper and the first computation of gravitational displacement memory in Generalized Proca gravity. The authors do real work: they build the gauge-invariant second-order action on an asymptotically flat background, derive the effective stress-energy tensor, and produce per-mode memory formulas that are internally consistent and carefully derived. The distinction between phase and group velocity in the dispersive case is genuinely new, and the Lorentz-violating branch is a useful cross-check against the earlier Einstein-Æther result.\n\nThe main soft spot is exactly where the stress-test note lands. Eq. (84) claims to give the total tensor memory for the massive, Lorentz-invariant case by integrating over group velocities β_χ from 0 to 1, but it supplies no spectral density and no Jacobian relating frequencies to β_χ. The flux F_χ is defined for a single mode, so the integral as written is not computable. The paper itself concedes in Sec. IV that the full memory depends on an explicit integration over frequency content, which confirms the gap. This does not break the per-mode formula (83) or the case-(b) results, but it means the headline \"full memory formula\" for case (a) is not yet a well-defined prediction.\n\nThe Lorentz-violating critical-angle enhancement is interesting, but the paper's own caveat that perturbation theory breaks down near the critical direction is important. The quantitative constraint discussion is also qualitative—it leans on existing bounds rather than deriving new ones. That is acceptable for a first paper, but it should be flagged as a limitation rather than presented as a final observational statement. The heavy self-citation is natural given that the framework is the authors' own; it is not a flaw here.\n\nOn balance, the per-mode results and the second-order action are solid contributions, and the flaw in Eq. (84) is repairable. This paper deserves a serious referee, not a desk rejection. A good referee should ask for a spectral treatment of the dispersive sum, and the authors should either provide it or explicitly present Eq. (83) as the per-mode result and demote Eq. (84) to a schematic expression.","headline":"First GP memory formulas, but the total dispersive memory in Eq. (84) is under-specified—worth refereeing, not desk-rejecting.","tokens_in":30308,"tokens_out":1570,"would_cite":true,"duration_ms":16797,"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":"This paper computes the first displacement-memory formula for Generalized Proca gravity and finds an unbounded critical-angle enhancement in the Lorentz-violating branch.","keywords":["gravitational memory","displacement memory","Generalized Proca gravity","Lorentz violation","group velocity","second-order action","stress-energy tensor","superluminal propagation"],"falsifier":"Compute the memory from a concrete compact-binary model in case (a) using the paper's own flux formulas: if one needs a Jacobian $d\\omega/d\\beta$ to convert the frequency integral into the $\\beta$-integral of equation (84), and that Jacobian changes the angular pattern or amplitude, the formula as written is incomplete. Observationally, in case (b), a detector network should see a sharp enhancement of memory in the specific direction $\\mathbf{n}'\\cdot\\mathbf{n}=\\beta_T/\\beta_\\psi$ for a known source; the absence of such an angular feature, combined with existing bounds on subluminal and superluminal propagation speeds, would exclude the enhancement mechanism.","tokens_in":29217,"feed_emoji":"🌌","tokens_out":9393,"duration_ms":83337,"temperature":0.7,"pith_summary":"Generalized Proca gravity is the most general metric theory with a gravitational vector field and derivative self-interactions that still keeps second-order equations of motion. The paper derives, for the first time, the permanent displacement memory left in the spacetime metric after a burst of gravitational radiation in this theory. It identifies two physically distinct asymptotic backgrounds: a Lorentz-invariant massive case, whose scalar and vector modes propagate slower than light and disperse, and a Lorentz-violating massless case whose modes travel at fixed speeds. In the Lorentz-violating case, any mode faster than the tensor graviton focuses its energy onto a single critical direction, making the memory amplitude grow without bound at that angle. The resulting formulas turn gravitational memory into a direct probe of the Proca mass and of Lorentz-violating propagation speeds.","feed_headline":"Proca gravity predicts memory blowup at one angle","feed_subtitle":"Proca gravity's first memory formula; a faster-than-light mode makes memory diverge at one critical angle","key_machinery":"The load-bearing machinery is the gauge-invariant scalar-vector-tensor decomposition of the perturbed action, which isolates five dynamical degrees of freedom (two tensor, two vector, one scalar) and reduces the second-order action to a sum of independent wave sectors. From that action, the second-variation method produces the gauge-invariant asymptotic stress-energy tensor for each sector. The memory integral then follows from the Isaacson back-reaction equation with a Green's function adapted to the physical group velocity $\\beta_\\psi$ rather than the phase velocity, which introduces the angular factor $V_\\psi = (1 - \\beta_\\psi \\,\\mathbf{n}'\\cdot\\mathbf{n}/\\beta_T)^{-1}$; the sign change of this factor when $\\beta_\\psi>\\beta_T$ is what produces the critical-angle divergence.","core_discovery":"The central discovery is the complete displacement-memory formula for Generalized Proca gravity, obtained by reducing the theory to its five dynamical gauge-invariant degrees of freedom and computing the effective stress-energy tensor that sources the low-frequency metric perturbation. In the massive Lorentz-invariant branch, the tensor memory is the general-relativistic result plus an integral over subluminal, dispersive scalar and vector modes, with the integrand carrying the group velocity of each mode and a denominator $1 - (\\beta_\\chi/\\beta_T)\\,\\mathbf{n}'\\cdot\\mathbf{n}$; because these modes are dispersive, the total memory must be assembled from all emission speeds $\\beta_\\chi\\in[0,1]$. In the Lorentz-violating branch, the scalar, vector, and tensor modes are non-dispersive, and whenever a source mode is faster than the tensor mode the memory amplitude diverges at the critical angle $\\mathbf{n}'\\cdot\\mathbf{n}=\\beta_T/\\beta_\\psi$. The paper states this as the first computation of this memory in Generalized Proca gravity and ties the unbounded enhancement to the same causal-structure mechanism seen in the earlier Lorentz-violating vector-tensor analysis.","pith_inferences":["The $\\beta_\\chi$ integral in equation (84) implicitly assumes a spectral distribution; a natural extension is to derive the Jacobian $d\\omega/d\\beta$ from a specific emission model, which would turn the formula into a ready-made waveform template for matched-filtering searches.","The critical-angle divergence is a kinematic resonance of any theory with superluminal propagation, so the same angular memory feature should appear in other spontaneously Lorentz-violating vector theories; a null search at the predicted angle would constrain the whole class, not just Generalized Proca gravity.","The discreteness of the Lorentz-violating branch suggests the two backgrounds could coexist in different spacetime regions; if so, transitions between branches might generate stochastic or burst-like memory signatures that future gravitational-wave observatories could search for.","Because the memory integral depends on group velocity while wavefront arrival depends on phase velocity, multi-messenger timing of a gravitational-wave memory signal with an electromagnetic counterpart could, in principle, measure the dispersion relation of the additional modes directly."],"forward_implications":["In the Lorentz-invariant massive branch, the memory signal from the Proca scalar and vector modes arrives late and spread out in time, because each frequency travels at its own subluminal group velocity; the shape of that tail encodes the Proca mass.","In the Lorentz-violating branch, any mode with $\\beta_\\psi>\\beta_T$ produces a memory amplitude that grows without bound at the critical angle, so past null results in memory searches translate directly into upper bounds on the Proca coupling parameters.","The polarization content differs sharply between branches: the massive branch shows only the two tensor polarizations, whereas the Lorentz-violating branch generically activates all six metric polarizations, with the vector polarizations present only when the tensor speed deviates from luminal.","Combined with the multimessenger bound that the tensor speed is within $10^{-15}$ of the speed of light and with Cherenkov-radiation bounds on subluminal modes, the memory result restricts the Lorentz-violating parameter space to near-luminal propagation, or to the special luminality conditions.","If the additional GP modes are excited in compact binary coalescences, the nondetection of memory in existing catalogs bounds the amplitude of the critical-angle enhancement, making Generalized Proca gravity testable with current detectors."],"supporting_citations":[{"why":"Defines the Generalized Proca action whose perturbation theory is analyzed.","marker":"[49]"},{"why":"Establishes the gauge-invariant second-order-action method for computing asymptotic memory from the stress-energy tensor.","marker":"[39]"},{"why":"Provides the Horndeski-like treatment of dispersive scalar modes that the massive branch (a) mirrors.","marker":"[40]"},{"why":"Gives the Lorentz-violating vector-tensor memory formulas, the critical-angle divergence, and the observational constraints that case (b) parallels.","marker":"[41]"},{"why":"Earlier calculation of leading-order Isaacson equations in GP gravity that this paper completes and extends.","marker":"[50]"},{"why":"The second-variation method used to extract the gauge-invariant energy-momentum tensor from the second-order action.","marker":"[62]"},{"why":"Introduces the nonlinear high-frequency averaging and the effective stress-energy tensor that define the memory back-reaction equation.","marker":"[68]"},{"why":"Null search for gravitational-wave memory in binary coalescence catalogs, used to bound the memory enhancement amplitude.","marker":"[76]"},{"why":"Multimessenger measurement bounding the tensor propagation speed near the speed of light.","marker":"[78]"},{"why":"Cherenkov-radiation argument ruling out subluminal propagation modes.","marker":"[79]"}],"fun_headline_variants":["First Proca gravity memory formula predicts angle-dependent blowup","Proca memory diverges at critical angle when source beats tensor wave","Generalized Proca gravity: superluminal modes make memory signal diverge","Proca gravity memory blowup at critical angle from superluminal mode"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The fragile step is the claim that the total memory from dispersive scalar and vector modes can be obtained simply by integrating the flux over all group velocities $\\beta_\\chi$ from 0 to 1 without specifying the radiation spectrum that maps frequency to $\\beta_\\chi$; without that spectral density, equation (84) leaves the total memory amplitude underdetermined.","fun_headline_variants_meta":{"raw":{"variants":["First Proca gravity memory formula predicts angle-dependent blowup","Proca memory diverges at critical angle when source beats tensor wave","Generalized Proca gravity: superluminal modes make memory signal diverge","Proca gravity memory blowup at critical angle from superluminal mode"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00079,"raw_usage":{"total_tokens":3500,"prompt_tokens":982,"completion_tokens":2518,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":598,"completion_tokens_details":{"reasoning_tokens":2451}},"tokens_in":598,"tokens_out":2518,"duration_ms":17200,"temperature":1.0,"reasoning_tokens":2451,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T16:43:50.816197+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the memory from a concrete compact-binary model in case (a) using the paper's own flux formulas: if one needs a Jacobian $d\\omega/d\\beta$ to convert the frequency integral into the $\\beta$-integral of equation (84), and that Jacobian changes the angular pattern or amplitude, the formula as written is incomplete. Observationally, in case (b), a detector network should see a sharp enhancement of memory in the specific direction $\\mathbf{n}'\\cdot\\mathbf{n}=\\beta_T/\\beta_\\psi$ for a known source; the absence of such an angular feature, combined with existing bounds on subluminal and superluminal propagation speeds, would exclude the enhancement mechanism.","supporting_citations":[{"cited_title":"Spontaneous Breaking of Lorentz Symmetry in String Theory,","cited_arxiv_id":null,"evidence_quote":"The second-variation method used to extract the gauge-invariant energy-momentum tensor from the second-order action."}],"review_version":2}