{"id":"6925f9b2-0be9-4e15-bdd8-e991a980248b","arxiv_id":"2607.08570","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"Direct-wave envelopes from black-hole mergers are predicted to decay at about 62% of the Kerr surface-gravity rate because a redshift-suppressed near-horizon source is convolved with the screened black-hole response, matching GW250114 and GW231226 residuals.","lead":"After a black-hole merger, the brief 'direct wave' fades more slowly than the black hole's surface gravity predicts; this paper argues the slowdown comes from gravitational redshift near the horizon shaping the signal, not from a different horizon property. The authors test this with a Teukolsky calculation and find their predicted damping agrees with residual data from two gravitational-wave events.","discovery_kind":"first_principles","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Source prescription λ_s≈κ is the load-bearing, unvalidated input: γ_eff/κ≈0.62 shifts if the true plunge-source decay differs, and the event fit grid-tunes source-decay time.","rationale":"The reader’s CONDITIONAL verdict is appropriate. The paper is methodologically sound in its complex-frequency Teukolsky kernel and the mathematical convolution argument is correct: two comparable decay rates naturally produce a finite-window logarithmic slope below κ, and this part is robust for a wide range of λ_s and γ_imp. However, the event-specific claim γ_eff/κ≈0.62, and the claimed consistency with GW250114 and GW231226, depend on a source model that is asserted, partly grid-tuned, and not reproducible because the code and data are not yet public. The near-horizon scaling argument motivates an exponential source but does not derive the actual radiation amplitude reaching infinity from a full inhomogeneous Teukolsky source term. The paper itself flags conditionality, so this is not a hidden flaw; it is the known weak point that justifies a conditional rather than an unconditional acceptance. No verdict change is needed beyond the reader’s CONDITIONAL.","tokens_in":9885,"tokens_out":9887,"duration_ms":97126,"concrete_test":"Derive S(t) from a full inhomogeneous spin-2 Teukolsky calculation using a plunging trajectory in Kerr obtained from the public GW250114 posterior samples, compute the convolution with the same screened kernel and Eq. (4) with no free source parameters, and compare the resulting γ_eff/κ to 0.62 and to the best-fit source-decay time. If the derived effective λ_s differs from κ by more than 20%, or if γ_eff/κ moves outside 0.59–0.66, the central numerical claim is not supported.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central numerical result γ_eff/κ≈0.62 is computed from Eq. (4) using α_s=λ_s/κ=1, an assumption justified only by the near-horizon scaling dx/dt≈−2κx and lapse √x∝e^{−κt}. That argument gives the redshift of a test frequency, not the amplitude of the outgoing radiation emitted by a plunging source as seen at infinity; in a full inhomogeneous Teukolsky calculation, S(t) is a source integral involving the trajectory and the homogeneous solution, and is not guaranteed to be e^{−κt}. The “near-horizon test-particle source” that supposedly provides independent confirmation in Table 1 is never specified in Methods, so it cannot be checked. Moreover, the GW250114 comparison in Methods performs a grid search over source-decay times (0.8–2.0 ms) and chooses the value that minimizes the whitened residual; the best fit is 2 ms, whereas κ=0.6306 ms⁻¹ corresponds to about 1.59 ms. Unless this parameter has a different meaning, the event-level fit implies α_s≈0.79, not the α_s=1 used in Eq. (4). Thus the event-level agreement is partly a fit, not a blind prediction, and the theoretical curve in Fig. 5 may not use the same source parameters as the data comparison. The paper’s own Discussion concedes the ratio is “conditional on the source prescription and window.” This does not refute the convolution mechanism, but it is the load-bearing assumption on which the headline quantitative claim rests.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper argues that the observed envelope damping of direct gravitational-wave radiation from a black-hole merger is not the bare Kerr surface gravity κ. The authors compute a spin-2, ℓ=m=2 Teukolsky response at complex frequency, define a screened response kernel with a zero at the complex horizon frequency, and convolve it with a finite-duration, near-horizon source whose amplitude decays as e^{-λ_s t} with λ_s≃κ. The resulting two-exponential envelope is then projected onto a finite log-envelope window, yielding the central result γ_eff/κ≈0.62. For GW250114 this gives γ_eff≈0.40 ms^{-1}, and for GW231226 γ_eff≈0.31 ms^{-1}, both presented as consistent with public residual data. The paper concludes that direct-wave damping is a horizon-redshift transfer observable rather than a direct measurement of κ.","tokens_in":10318,"tokens_out":5437,"duration_ms":50285,"significance":"If the central mechanism is correct, the paper offers a new and important way to read measured direct-wave envelopes: a source–response convolution with comparable decay rates naturally produces a finite-window slope below κ, so observing γ_eff<κ does not contradict Kerr. The analytic formula in Eq. (4) is clean, parameter-free once the source rate, kernel damping, and window are fixed, and it makes a specific, falsifiable prediction for the ratio γ_eff/κ across spins and carrier frequencies. The use of a published complex-frequency Teukolsky solver and the attempt at two-event residual consistency checks are strengths. However, the quantitative prediction and the event-level consistency depend critically on the phenomenological source decay rate λ_s=κ, which is asserted from near-horizon scaling rather than derived from a sourced Teukolsky calculation, and the GW250114 comparison grid-selects a source-decay time that appears inconsistent with α_s=1. These issues must be resolved before the central claim can be accepted.","major_comments":[{"comment":"The load-bearing input is the source prescription S(t)=e^{-λ_s t} with λ_s=κ. The justification (dx/dt≈−2κx, lapse ∝√x) gives the redshift of a test frequency near the horizon, but not the amplitude of outgoing radiation emitted by a plunging source and propagated to infinity in an inhomogeneous Teukolsky calculation. In the latter, S(t) is a source integral involving the orbit and the homogeneous solution, and need not be a pure exponential with rate κ. Since Eq. (4) returns γ_eff/κ≈0.62 only for α_s=1, the headline number is conditional on an unvalidated source model. The paper acknowledges this conditionality in the Discussion, but a specific derivation or a quantitative robustness scan over λ_s is needed. For illustration, inserting α_s=0.79 in Eq. (4) gives γ_eff/κ≈0.49, quite different from 0.62.","section":"Physical mechanism of the reduced damping, Eqs. (2) and (4)"},{"comment":"The grid search over source-decay times includes τ_s=0.8, 1.2, 1.5, 2.0 ms and selects 2.0 ms as the minimum. If 'source-decay time' is the e-folding time, then with κ=0.6306 ms^{-1} the implied α_s=(κτ_s)^{-1}=(0.6306×2.0)^{-1}=0.79, not α_s=1. The theoretical curve used in Fig. 2 and the value quoted from Eq. (4) use α_s=1. The paper does not explain this discrepancy. If τ_s has a different meaning (e.g., a rise-plus-decay timescale), it must be defined explicitly. As written, the event-level agreement is partly a fit to a different source rate than the one used in the central prediction.","section":"Methods, 'GW250114 residual scan'"},{"comment":"The 'near-horizon test-particle source' is listed as an independent finite-duration realization giving γ_eff≈0.40 ms^{-1}, but no equations, trajectory, or numerical ingredients are given anywhere in Methods. This source cannot be checked, and if it shares the same λ_s=κ assumption it is not independent of the phenomenological plunge source. Please specify the model and show that it does not reduce to the same ansatz.","section":"Table 1 and Methods"},{"comment":"The fitting window [u1,u2]=[1.26,5.04] is inherited from the GW250114 residual analysis (2–8 ms with κ=0.6306 ms^{-1}). For other remnants it is rescaled as u=κt, which changes the physical window. Because γ_eff/κ is a window average, the paper should quantify sensitivity to the window choice: e.g., a small plot or table of γ_eff/κ versus (u1,u2) or versus physical window would show how much of the claimed 0.62 is an artifact of the chosen interval. Similarly, the GW231226 profile peaks at η≈0.56 while the prediction is η≈0.63; the paper says the prediction lies within the 'broad maximum' but does not quantify the profile likelihood at the predicted point. A numerical likelihood ratio would turn this visual consistency into a test.","section":"Finite-window damping projection and Fig. 5b"}],"minor_comments":[{"comment":"Typos: in the Abstract, 'spin-$-2$' should be 'spin−2'; in the Introduction, 'GW250114 and GW231226' residuals are described as QNM-subtracted, but the precise waveform model used for subtraction is only given later. Also, '0.18 below the maximum' in the text near Fig. 2 lacks units or a definition of the profile quantity.","section":"Abstract and Introduction"},{"comment":"The quoted interval z=0.230^{+0.041}_{-0.063} has the upper and lower errors in nonstandard order relative to the other quoted intervals; please check consistency.","section":"Eq. (5)"},{"comment":"Ref. [24] contains 'doi: 10.1103/kkmt-fbjb', which is not a valid DOI and appears to be a placeholder. Ref. [13] also has a volume-less DOI format. Please update these.","section":"References"},{"comment":"The statement 'will be deposited in a versioned public repository upon publication' is not verifiable during review. For a paper whose validation relies on public residual data and numerical code, please provide the code/data repository or a detailed description of the analysis scripts.","section":"Data availability"}],"recommendation":"major_revision","confidential_remarks":"I agree with the stress-test assessment: the source prescription λ_s=κ is the central vulnerability. The convolution mechanism in Eq. (4) is sound, but the quantitative claim is not established without either a derivation of the source amplitude from a full inhomogeneous Teukolsky calculation (or NR-informed model) or a demonstration that the event-level results are insensitive to α_s. The apparent mismatch between the grid-selected source-decay time (2 ms) and the α_s=1 used in the theory is the most concrete fixable issue. If the authors can reconcile this and specify the test-particle source, the paper would be a solid contribution. As it stands, it is too conditional for acceptance."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper deserves a serious referee and a read, but the quantitative ratio γ_eff/κ≈0.62 is a model output resting on an asserted source decay law, not a measured or derived fact. The conceptual mechanism is new and plausible; the event claims are looser than they look.\n\nWhat is new: Han and Jiang split the carrier-frequency question from the envelope-damping question and argue that direct-wave envelope damping should be viewed as the convolution of a redshift-suppressed near-horizon source with a screened Kerr response. The analytic two-rate projection, Eq. (4), is clean and makes the mechanism transparent. The complex-frequency Teukolsky kernel comes from their published solver, so that piece is reproducible. The GW231226 residual profile is a genuine second-event consistency check, and the model locus intersects the 68% region in both events.\n\nWhat is soft: The load-bearing input is λ_s≈κ. The heuristic dx/dt≈−2κx and lapse √x∝e^{−κt} describe the redshift of a test frequency, not the amplitude of the outgoing radiation; in a full inhomogeneous Teukolsky calculation, the effective source integral could decay differently. If λ_s differs, the ratio shifts, and in the impulsive limit you recover γ_imp≈0.94κ. The paper itself concedes the ratio is conditional on the source prescription and window, which is honest but undercuts the idea of a sharp prediction.\n\nThere is also a tension between the theory and the event fit. The Methods grid-search over source-decay times picks 2 ms, whereas α_s=1 gives about 1.59 ms. That suggests the data comparison is using α_s≈0.79, not 1, and the paper does not reconcile this with Eq. (4). The 'near-horizon test-particle source' in Table 1 is never specified, so it cannot be checked. No code or data is public yet.\n\nThese are not fatal to the mechanism, but they mean the agreement with GW250114 is partly a fit, not a blind prediction. The stress-test note is right on this.\n\nFor whom: This is for ringdown spectroscopists and black-hole perturbation people. It reframes the direct-wave debate and gives a concrete mechanism for sub-κ damping. I'd bring it to a reading group, with the source-model caveat flagged in advance.\n\nBottom line: Send it to peer review. A good referee should push for a derivation of the source decay from an inhomogeneous Teukolsky or NR calculation, and for the event comparison to use the same source parameters as the prediction. But the core idea is solid enough to deserve referee time.","headline":"A plausible two-rate convolution mechanism for direct-wave damping, with a clean analytic core, but the headline γ_eff/κ≈0.62 rests on an asserted source decay law and an event fit that may be using different source parameters.","tokens_in":10789,"tokens_out":3718,"would_cite":true,"duration_ms":34144,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["83C57","83C35"],"pacs":["04.30.-w","04.70.-s"],"model":"deepseek-v4-flash","headline":"The observed envelope of direct waves from a black-hole merger is set by the convolution of a redshifted plunge source with the screened Kerr response, not by surface gravity alone.","keywords":["gravitational waves","black-hole ringdown","direct waves","Kerr surface gravity","Teukolsky equation","horizon redshift","effective damping","GW250114"],"falsifier":"Run a full inhomogeneous Teukolsky or numerical-relativity simulation that supplies the actual plunge source amplitude, apply the same 2–8 ms log-envelope fit to the resulting direct-wave strain, and check the fitted slope; if the slope is close to κ (≈0.63 ms^-1) rather than ≈0.40 ms^-1 for a GW250114-scale remnant, the λ_s ≈ κ source prescription is wrong. Alternatively, use a third merger with the same final spin χ ≈ 0.67 but a different remnant mass: the model predicts γ_eff ∝ 1/M, so the measured rate in ms^-1 should move accordingly, while a constant rate would refute it.","tokens_in":9740,"feed_emoji":"🕳️","tokens_out":6985,"duration_ms":61329,"temperature":0.7,"pith_summary":"This paper argues that the observed envelope damping of direct gravitational waves from black-hole mergers is not a direct reading of the Kerr surface gravity κ. Because a finite-duration near-horizon source has its outgoing amplitude suppressed by gravitational redshift, and because that source decays at a rate comparable to the screened response, the product is a non-exponential envelope with an effective damping γ_eff below κ. Convolving a complex-frequency Teukolsky response kernel with such a source gives γ_eff/κ≈0.62, corresponding to about 0.40 ms^-1 for GW250114 and 0.31 ms^-1 for GW231226, both consistent with residual analyses of those events. If right, this turns direct-wave damping into an observable of horizon-redshift transfer and requires any inference of κ from envelopes to account for the source's finite duration.","feed_headline":"Direct-wave damping falls to 62% of the Kerr surface-gravity rate","feed_subtitle":"A redshifted plunge source, not the horizon alone, sets the damping seen in GW250114 and GW231226.","key_machinery":"The load-bearing object is a two-rate convolution: the screened complex-frequency Teukolsky response kernel I(t) and a redshift-suppressed plunge source S(t). The kernel is built from the spin-2, ℓ=m=2 Teukolsky equation solved at complex frequency along a contour with Imω = −κ, with a screening zero placed at the complex horizon frequency mΩ_H − iκ and a Gaussian window centered on the real carrier ω_c. The source is a finite-duration exponential S(t) = e^{-λ_s t} e^{-iω_c t} whose rate λ_s ≈ κ is motivated by near-horizon redshift. Convolving the two turns the product into a difference of exponentials, and a finite-window log-envelope projection converts that difference into the dimensionl","core_discovery":"The paper's central claim is that a black-hole merger's direct-wave strain, viewed as a carrier times a slow envelope, obeys h_DW(t) ∝ (e^{-γ_imp t} − e^{-λ_s t})/(λ_s − γ_imp) e^{-iω_c t} rather than a single e^{-κt}. Here I(t) = e^{-γ_imp t} e^{-iω_c t} is the screened Teukolsky response kernel (with γ_imp ≈ 0.94κ from the complex-frequency contour and Gaussian window), and S(t) = e^{-λ_s t} e^{-iω_c t} is the finite-duration near-horizon source, with λ_s ≈ κ following from the near-horizon redshift scaling dx/dt ≈ −2κx and √x ∝ e^{-κt}. The difference of two exponentials produces a broader envelope; projecting its log slope over the 2–8 ms post-peak window gives γ_eff/κ ≈ 0.62. For the re","pith_inferences":["If the γ_eff/κ ≈ 0.62 ratio persists, it offers a simple explanation for why some numerical-relativity analyses find filtered-strain damping that drifts away from κ: the measured slope is a window-dependent convolution of two comparable rates, not a clean κ.","A testable extension: at fixed remnant spin, γ_eff should scale with κ ∝ 1/M, so a third event with χ ≈ 0.67 but a different remnant mass would discriminate the model from a fixed damping rate; the two current events already test mass rescaling but not spin dependence.","The source prescription λ_s ≈ κ is heuristic; a full inhomogeneous Teukolsky computation or a numerical-relativity plunge trajectory could directly measure λ_s and either confirm the ≈0.62 prediction or show that the real ratio is closer to the impulsive limit 0.94.","Because the paper decouples the real carrier frequency from the horizon-frequency reference, it implies that carrier-phase matching is not required for the envelope test; future searches could target envelope damping alone, sidestepping the debate over accidental carrier-frequency crossings."],"forward_implications":["Envelope damping measured in direct-wave residuals should be compared with the source-convolved effective rate γ_eff, not with the bare Kerr rate κ; single damped-sinusoid fits that identify the envelope slope with κ will misread the signal.","For remnant spins 0.50–0.95 and carrier detuning −2 ≤ (ω_c − mΩ_H)/κ ≤ 1, the model predicts γ_eff/κ between 0.589 and 0.663, with the ratio near 0.63 for spins around 0.67.","GW250114's residual profile places the best-fit source-convolved model at f_c ≈ 171 Hz and γ_eff ≈ 0.410 ms^-1, while the bare-κ reference at f_H lies outside the nominal 95% region; GW231226's damping profile peaks at η ≈ 0.56, with the predicted η ≈ 0.63 inside its broad maximum.","In the impulsive-source limit, the fitting projection recovers γ_imp ≈ 0.94κ, so the reduced damping is specifically a finite-duration source effect rather than an artifact of the fitting window.","Joint inference of remnant parameters and source duration could ultimately allow κ to be extracted from measured damping without assuming γ_eff = κ."],"fun_headline_variants":["Horizon-redshift clocks down direct-wave damping","Direct-wave damping set by redshifted plunge, not bare Kerr","Direct-wave decay reveals redshift transfer, not surface gravity","Black-hole direct waves decay slower than Kerr predicts"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The quantitative prediction rests on the source model S(t) ≈ e^{-κt} for the near-horizon plunge amplitude; if the actual source decays at a different rate or shuts off before the measurement window, the two-rate convolution no longer gives γ_eff/κ ≈ 0.62, and the impulsive limit instead approaches 0.94κ.","fun_headline_variants_meta":{"raw":{"variants":["Horizon-redshift clocks down direct-wave damping","Direct-wave damping set by redshifted plunge, not bare Kerr","Direct-wave decay reveals redshift transfer, not surface gravity","Black-hole direct waves decay slower than Kerr predicts"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000778,"raw_usage":{"total_tokens":3290,"prompt_tokens":773,"completion_tokens":2517,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":517,"completion_tokens_details":{"reasoning_tokens":2454}},"tokens_in":517,"tokens_out":2517,"duration_ms":17088,"temperature":1.0,"reasoning_tokens":2454,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-04T04:23:46.008023+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run a full inhomogeneous Teukolsky or numerical-relativity simulation that supplies the actual plunge source amplitude, apply the same 2–8 ms log-envelope fit to the resulting direct-wave strain, and check the fitted slope; if the slope is close to κ (≈0.63 ms^-1) rather than ≈0.40 ms^-1 for a GW250114-scale remnant, the λ_s ≈ κ source prescription is wrong. Alternatively, use a third merger with the same final spin χ ≈ 0.67 but a different remnant mass: the model predicts γ_eff ∝ 1/M, so the measured rate in ms^-1 should move accordingly, while a constant rate would refute it.","supporting_citations":[],"review_version":2}