{"id":"a8db19ec-1952-4ff4-9083-a13f1fd88481","arxiv_id":"2501.08964","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"Infalling and outgoing gravitational radiation around a perturbed Schwarzschild horizon share the same quasi-normal mode spectrum, with an explicit transfer formula between the horizon and null infinity.","lead":"The paper derives, in linear perturbation theory around a Schwarzschild black hole, that the radiation falling into the horizon and the radiation escaping to infinity ring with the same quasi-normal mode frequencies, and it gives explicit formulas connecting the two. This provides the first analytical support for 'black hole tomography', the program of inferring horizon geometry from gravitational wave observations.","discovery_kind":"first_principles","skeptic_critique":{"model":"deepseek-v4-flash","headline":"No-incoming-radiation check at I− is deferred and likely ill-posed because the QNM solution diverges at I−; Eq. (138) may relate only formal modes.","rationale":"The reader's weakest_assumption is exactly the deferred check of Ψ1 peeling at I−. I agree, but the concern is sharper: the paper's own asymptotic expressions show the QNM solution diverges at past null infinity, so the claimed check is not merely a lengthy computation—it is a check of whether a well-defined I− limit exists at all. The paper flags this in Sec. V ('A complete discussion of this behavior is beyond the scope... tentative proposal'). This makes the no-incoming-radiation condition the hinge of the paper: if it fails, the analyticity/stability selection is not the standard QNM boundary condition, and Eq. (138) does not relate detected ringdown amplitudes to horizon infall. The rest of the paper—Teukolsky separation, horizon algebra, Starobinsky consistency—is substantial and appears internally coherent; those parts support the formal mode relation. My verdict remains CONDITIONAL, matching the reader, because the gap is explicitly acknowledged and a concrete calculation can settle it. I am not raising an external-consensus objection; the issue is the paper's own unverified assertion.","tokens_in":49601,"tokens_out":7194,"duration_ms":76434,"concrete_test":"Perform the asymptotic analysis at I−: take the full solutions for eΨ0 and eΨ4 (Eqs. (97)/(119)) with k^{(2)}=0, express all spin coefficients via the radial Bianchi identities and frame equations (App. A), and compute the leading 1/r behavior of Ψ1 along v=const null generators as r→∞, specifying whether the overtone sum or the limit is taken first. Accept the boundary condition only if Ψ1 = Ψ^{(0)}_1 Ω^2 + O(Ω^3) with Ω=-1/r as a well-defined asymptotic series. If the limit is divergent or the leading power differs, the no-incoming-radiation condition fails and the QNM identification collapses.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The weakest point is the asserted verification of the no-incoming-radiation condition from I−. In Sec. VB the authors state that the second sufficient condition of [70]—the peeling behavior of Ψ1 at I−—'has been checked' but the calculation is deferred. This is not a cosmetic omission: it is the step that identifies the analyticity/stability-selected solution as the physical ringdown with no incoming radiation, and hence licenses Eq. (138) as a statement about detected amplitudes. The paper itself shows the difficulty: eΨ0 behaves as e^{iω(-v+2r)}/r^5 (Eq. (120)) and 'the exponential terms in fact diverge at i0 and at past-null infinity'; the authors then 'tentatively' take polynomial fall-off as implementing no-incoming radiation. The Walker-Will conditions, however, assume existence and regularity of I−, which a diverging QNM solution does not obviously possess. Asserting the Ψ1 peeling without presenting the asymptotic limit—which must be taken in the r→∞, v=const sector where the mode sum diverges—leaves the physical boundary condition unestablished. Consequently, the central tomography formula may be a relation among formal modes only.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper develops a characteristic initial-value formulation of linear perturbations of a slowly spinning isolated horizon, with free data on the horizon encoded in the perturbed Weyl scalar Ψ0 and outgoing radiation registered on a transverse null surface N that is eventually identified with future null infinity. The horizon transport equations (41) are solved explicitly in Eq. (54), the Teukolsky equations for Ψ0 and Ψ4 are separated in the adapted horizon-penetrating coordinates, and the requirement of analyticity plus future stability is shown to select the standard Schwarzschild QNM frequencies through the Leaver continued-fraction condition (84). The central result is the transfer relation (138), with explicit coefficients in Eqs. (139)-(140) and inverses in Eqs. (C4)-(C5), linking Ψ4 modes at I+ to Ψ0 modes at the horizon; the paper also provides an explicit expression for Ψ2 in Sec. VII. The presentation is constructive and the core algebraic structure is transparent, but the physical boundary-condition step that identifies the formal QNM solution with the no-incoming-radiation ringdown is not fully established.","tokens_in":49707,"tokens_out":6459,"duration_ms":75045,"significance":"If the construction is sound, this is a valuable analytic explanation of the numerically observed correlation between infalling horizon radiation and outgoing ringdown radiation, and it gives a concrete, falsifiable transfer function for black-hole tomography. Strengths include the explicit solution of the horizon equations, the demonstration that the same QNM frequencies arise for Ψ0 and Ψ4 through the recurrence identities (116), and the explicit, convergent numerical estimates for the transfer coefficients F_lmn. The paper is also careful to flag the slicing dependence of the transfer function. The main weakness is that the no-incoming-radiation boundary condition, which is load-bearing for the physical interpretation of Eq. (138), is asserted without presenting the required calculation and is in tension with the divergent behavior of the QNM solution at I-.","major_comments":[{"comment":"The identification of the analytic/stability-selected solution with the physical ringdown having no incoming radiation from I- rests on the second Walker-Will sufficient condition, the peeling behavior of Ψ1 at I-. This check is deferred: the text states that the procedure is straightforward but lengthy and 'will be presented explicitly elsewhere,' and that the authors have 'indeed checked' the condition. This is a load-bearing omission: the no-incoming-radiation condition is what licenses Eq. (138) as a relation between detected amplitudes at I+ and horizon amplitudes, rather than a relation among formal modes. The calculation should be presented in the paper, or the physical claim should be weakened accordingly.","section":"Sec. V B, paragraph following Eq. (128)"},{"comment":"The statement that 'as a tentative proposal for this paper, we shall take the polynomial fall-off of eΨ0 as implementing a no-incoming radiation from I-' conflicts with the definitive assertion in Sec. V B. As the authors note, the QNM solution contains factors e^{iω(-v+2r)} with Im ω < 0, which diverge as v → -∞; consequently the mode sum does not define a regular field on I-, and the Walker-Will peeling argument, which assumes existence and regularity of I-, cannot be applied without a further limiting prescription. The paper should specify the precise asymptotic procedure under which the Ψ1 peeling is evaluated and show that the result is independent of the order of limits.","section":"Sec. V A, paragraph after Eq. (120)"},{"comment":"The conclusion states that analyticity and future stability select the 'unique solution' describing the ringdown. This is stronger than what is demonstrated: the explicit solution (54) is obtained after discarding the growing shear mode c1 e^{κ v}, setting the static tidal solution to zero, restricting to the Fourier ansatz (52), and later imposing k(2)=0 by stability. The paper does not prove that no other solutions in the stated solution class survive these conditions. Either a precise uniqueness statement with proof should be supplied, or the wording should be softened to describe the solution that satisfies the stated additional assumptions.","section":"Sec. VIII, Conclusions"}],"minor_comments":[{"comment":"The displayed equalities appear to contain index typos: the last three terms use l0n where the context requires lmn (c−(2)_{l0n}, b+(2)_{l0n}, c+(2)_{l0n} should be c−(2)_{lmn}, b+(2)_{lmn}, c+(2)_{lmn}).","section":"Eq. (109)"},{"comment":"The denominator in the m≠0 expression uses ω_{l0n} in the factors (κ² + ω²) and (2κ - iω); consistency with Eq. (C5a) indicates these should be ω_{lmn}.","section":"Eq. (140a)"},{"comment":"The symbol eΨ0 is used in the abstract before its definition as the perturbation symbol in Sec. III; a brief parenthetical definition would help the reader.","section":"Abstract and Introduction"},{"comment":"There are several typos, including 'analicity' for 'analyticity' in the heading of Sec. V A, 'perturation' for 'perturbation' near the end of the Introduction, and 'Ilustration' in the caption of Fig. 3.","section":"Sec. V A, heading and Introduction"},{"comment":"The numerical values of F_lmn in Table I are presented as averages over N ∈ [30, 200], and the text notes numerical instability for large N. If these constants are intended for actual tomographic estimates, error bars from the truncation and from the approximated QNM frequencies should be provided or at least quantified in order of magnitude.","section":"Sec. VI, Table I"}],"recommendation":"major_revision","confidential_remarks":"The core derivation is largely convincing and the transfer formula (138) is an explicit, checkable result. The decisive issue is the deferred verification of the no-incoming-radiation boundary condition at I-; because the QNM solution diverges at I-, this is not a cosmetic omission but a point that can affect whether Eq. (138) applies to the physically detected ringdown. I would be willing to recommend acceptance once that calculation (or a precise weakened statement of the boundary condition) is included. I do not see grounds for rejection, since the structural relation between Ψ0 at the horizon and Ψ4 at I+ is derived explicitly from the field equations."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: this paper deserves a serious referee. The core analytic machinery is substantial: a first-principles 4D characteristic formulation on a perturbed isolated horizon, explicit solutions of the horizon equations, and a derivation that infalling Psi0 and outgoing Psi4 share the same QNM spectrum, with concrete transfer formulas (Eqs. 138-140). The QNM frequencies are recovered from the Leaver continued fraction, not fitted, and the connection to the Teukolsky-Starobinsky identities is made cleanly. The authors also honestly flag the slicing dependence of their F coefficients and give only order-of-magnitude numerical values. That is a lot of solid work, and the central Psi0-to-Psi4 link is structurally coherent.\n\nThe softest spot is exactly where the reader and the stress-test point: Section VB. The claim that the solution has no incoming radiation from I- is load-bearing—it is what upgrades a formal mode relation into a statement about detectable amplitudes at I+ and physics at H. The authors defer the check of the Psi1 peeling condition and also acknowledge that the QNM exponential blows up at i0 and I-. That is not a cosmetic gap. If the peeling check cannot be performed because I- is not regular for QNMs, then the Walker-Will criteria may not apply, and the boundary condition is, as the stress-test says, possibly ill-posed. I would not call it a demonstrated error; it is a missing step that needs to be supplied or qualified.\n\nThe \"unique solution\" claim in Sec. VIII is also stronger than what is shown: it is uniqueness within a Fourier-mode, separable, analytic ansatz, not in the full space of perturbations. That is a minor overreach. The F_lmn values in Table I are explicitly order-of-magnitude and slicing-dependent; acceptable as long as the tomography claim is read accordingly.\n\nBottom line: this is a paper for the gr-qc community working on ringdown, horizon physics, and perturbation theory. It will be cited, and it deserves conditional acceptance: send it to referees, and have the authors provide the deferred I- calculation or rewrite the boundary-condition discussion with a precise sense in which no-incoming-radiation holds. I would bring it to a reading group.","headline":"A genuinely new analytic bridge between horizon influx and outgoing ringdown modes, with a load-bearing no-incoming-radiation check that is asserted rather than shown.","tokens_in":50396,"tokens_out":5444,"would_cite":true,"duration_ms":59679,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["83C57","83C35","83C05","83C60"],"pacs":["04.30.-w","04.70.-s"],"model":"deepseek-v4-flash","headline":"An analytic transfer formula now connects the ringdown waves a detector sees to the radiation falling into the remnant black hole.","keywords":["black hole perturbation theory","quasi-normal modes","isolated horizons","gravitational wave ringdown","Newman-Penrose formalism","Teukolsky equation","black hole tomography","characteristic initial value formulation"],"falsifier":"Carry out the deferred calculation: expand the radial solution for $\\widetilde{\\Psi}_1$ at $\\mathcal{I}^-$ and check whether it decays as $\\Omega^2$ (equivalently as $r^{-3}$ in the Bondi frame); if it does not, the no-incoming-radiation condition fails. A numerical counterpart would be to measure the ratio of horizon to infinity QNM amplitudes in a binary black-hole ringdown simulation and compare with $F_{\\ell mn}$ of Table I; a disagreement beyond the reported convergence error would invalidate the transfer formula.","tokens_in":49313,"feed_emoji":"🕳️","tokens_out":8658,"duration_ms":85904,"temperature":0.7,"pith_summary":"This paper establishes, in linear perturbation theory around a Schwarzschild (slowly spinning) black hole, what numerical simulations have long suggested: the radiation falling into a remnant black hole during ringdown is a superposition of the same quasi-normal modes as the outgoing gravitational-wave signal. The authors reformulate the ringdown problem on null surfaces—the horizon and a transverse outgoing null surface—rather than on spacelike slices, and show that demanding analyticity and future stability of the perturbation selects the standard Schwarzschild quasi-normal frequencies without imposing the usual no-incoming-radiation boundary condition. The central result is an explicit mode-by-mode transfer formula, Eq. (138), linking the amplitude of $\\Psi_4$ at future null infinity to the amplitude of $\\Psi_4$ at the horizon, and through Eqs. (139)–(140), to the horizon data $\\widetilde{\\Psi}_0$. If correct, this gives an analytic foundation for “black hole tomography”: the horizon geometry and infalling flux could be inferred from observed ringdown amplitudes.","feed_headline":"Ringdown waves at infinity fix the infalling radiation at the horizon","feed_subtitle":"One coefficient links detected ringdown waves to horizon infall, making black hole tomography concrete.","key_machinery":"The load-bearing mechanism is the characteristic initial value formulation tailored to a perturbed isolated horizon: the horizon is a null surface $\\mathcal{H}$ on which the free data is the spin-weight-2 Weyl scalar $\\widetilde{\\Psi}_0$ (equivalently the shear), and a transverse null surface $N$ registers the outgoing $\\widetilde{\\Psi}_4$. The Newman–Penrose field equations at the horizon reduce to a closed system, Eqs. (41), that propagates $\\widetilde{\\Psi}_0$ through the spin coefficients and Weyl scalars up to $\\widetilde{\\Psi}_4$; the Teukolsky equation then carries each mode radially. Separability plus analyticity of the confluent-Heun series pins the frequencies to the roots of the continued fraction (84), i.e. the standard QNM frequencies, and the transfer coefficient $F_{\\ell mn} = \\sum_k a_k$ is the value of the radial series at $\\mathcal{I}^+$ relative to the horizon. The hyperboloidal coordinate transformation (132) makes the horizon-to-infinity map explicit and shows its slicing dependence.","core_discovery":"The paper's central discovery is that the correlation between horizon and infinity is not heuristic: it is encoded in a transfer coefficient. Working in a characteristic initial value problem with data on a perturbed isolated horizon and a transverse null surface, the authors solve the horizon field equations sourced by $\\widetilde{\\Psi}_0$, then solve the Teukolsky equations radially, and impose analyticity and stability toward the future. This selects the QNM frequencies $\\omega_{\\ell mn}$ of Schwarzschild and leaves the radiation at $\\mathcal{I}^+$ written as $\\Psi_4^{\\mathcal{I}^+,\\pm}_{\\ell mn} = F_{\\ell mn}(\\omega_{\\ell mn},c)\\, \\Psi_4^{\\mathcal{H},\\pm}_{\\ell mn}$ (Eq. (138)), where $F_{\\ell mn}$ is a convergent sum over the series coefficients of the radial solution. Inverting this relation gives the horizon amplitudes $\\Psi_{0,\\ell mn}^{\\mathcal{H},\\pm}$ in terms of the observed $\\Psi_4^{\\mathcal{I}^+,\\pm}$ (Appendix C), completing the analytic “tomography” map. The paper also shows that the same QNM frequencies appear in $\\widetilde{\\Psi}_0$, $\\widetilde{\\Psi}_4$, and all intermediate Weyl scalars, and that the solution satisfies the no-incoming-radiation condition at $\\mathcal{I}^-$; the verification of the $\\Psi_1$ peeling condition is stated but its calculation is deferred.","pith_inferences":["A direct test of the paper's transfer formula would be to extract horizon and infinity QNM amplitudes from numerical relativity simulations of a head-on merger or ringdown and compare their ratio to $F_{\\ell mn}$; the paper reports $F$ values but does not perform this comparison.","The deferred $\\Psi_1$ peeling check is the natural spot to probe: until it is published, the claim that analyticity and stability alone select the physical no-incoming-radiation solution rests on an assertion, and a failure there would not break Eq. (138) but would change which solution it describes.","Because $F_{\\ell mn}$ depends on slicing, the notion of “the” horizon amplitude inferred from a waveform is gauge-dependent; a future extension might identify a gauge-invariant combination, such as energy flux ratios, that removes this ambiguity.","The same mode-correlation idea may extend to extreme-mass-ratio inspirals, where the horizon flux is driven by the orbiting companion; the paper announces this as forthcoming work, but if the transfer map survives, tidal heating rates could be measured from the inspiral waveform."],"forward_implications":["A detected ringdown mode at $\\mathcal{I}^+$ fixes the amplitude of the corresponding infalling mode at the horizon through $F_{\\ell mn}$, so horizon flux and shear could be inferred from gravitational-wave observations.","The QNM frequencies arise from analyticity and future stability alone, which explains why the late-time ringdown and the horizon flux share the same frequencies and damping times.","In the minimal gauge, $|F_{\\ell 00}| < 1$ for the fundamental $\\ell=2,3,4$ modes, so the corresponding horizon amplitudes are larger than the observed amplitudes at infinity; this is a quantitative prediction that numerical-relativity simulations can test.","The map between horizon and infinity is slicing dependent, so any observational tomography claim must specify the height function; the same frequencies survive a change of slicing, but the amplitude relation does not.","The same construction applies, with minor changes, to a general isolated-horizon background, indicating that the correlation is a feature of horizon boundary data rather than of Schwarzschild-specific coordinates."],"supporting_citations":[{"why":"Companion paper that sets up the perturbed isolated horizon tetrad and gauge used here; supplies the horizon data construction.","marker":"[27]"},{"why":"Establishes separability of the Teukolsky equation, which the mode decomposition relies on.","marker":"[58]"},{"why":"Supplies the continued-fraction method and tabulated QNM frequencies that the paper's eigenfrequency condition reproduces.","marker":"[64]"},{"why":"Provides the series representations and three-term recurrence relations used to build radial solutions and to identify QNM eigenfrequencies.","marker":"[69]"},{"why":"Gives the peeling fall-off conditions at $\\mathcal{I}^-$ used to argue there is no incoming radiation; the paper defers the explicit $\\Psi_1$ check.","marker":"[70]"},{"why":"Numerical evidence that horizon fluxes are QNM superpositions; the empirical observation the paper explains analytically.","marker":"[21]"},{"why":"Recent numerical study correlating horizon and outgoing radiation; another target of the analytic explanation.","marker":"[22]"},{"why":"Characteristic initial value formulation framework that justifies prescribing data on null surfaces.","marker":"[23]"}],"fun_headline_variants":["Analytic map links horizon infall to ringdown waves at infinity","Horizon tomography made exact: detected ringdown fixes infall","QNM frequencies proven shared between horizon and null infinity","One transfer coefficient predicts horizon infall from ringdown data","Gravitational wave observations now yield analytic horizon snapshot"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the solution selected by analyticity and future stability really has no incoming radiation from past null infinity; the paper states that the required peeling fall-off of $\\Psi_1$ has been checked but defers the calculation, so if that check fails the physical interpretation changes.","fun_headline_variants_meta":{"raw":{"variants":["Analytic map links horizon infall to ringdown waves at infinity","Horizon tomography made exact: detected ringdown fixes infall","QNM frequencies proven shared between horizon and null infinity","One transfer coefficient predicts horizon infall from ringdown data","Gravitational wave observations now yield analytic horizon snapshot"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000169,"raw_usage":{"total_tokens":1299,"prompt_tokens":1014,"completion_tokens":285,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":630,"completion_tokens_details":{"reasoning_tokens":203}},"tokens_in":630,"tokens_out":285,"duration_ms":3672,"temperature":1.0,"reasoning_tokens":203,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T20:13:07.645259+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Carry out the deferred calculation: expand the radial solution for $\\widetilde{\\Psi}_1$ at $\\mathcal{I}^-$ and check whether it decays as $\\Omega^2$ (equivalently as $r^{-3}$ in the Bondi frame); if it does not, the no-incoming-radiation condition fails. A numerical counterpart would be to measure the ratio of horizon to infinity QNM amplitudes in a binary black-hole ringdown simulation and compare with $F_{\\ell mn}$ of Table I; a disagreement beyond the reported convergence error would invalidate the transfer formula.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the continued-fraction method and tabulated QNM frequencies that the paper's eigenfrequency condition reproduces."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the series representations and three-term recurrence relations used to build radial solutions and to identify QNM eigenfrequencies."},{"cited_title":"Walker and C","cited_arxiv_id":null,"evidence_quote":"Gives the peeling fall-off conditions at $\\mathcal{I}^-$ used to argue there is no incoming radiation; the paper defers the explicit $\\Psi_1$ check."}],"review_version":1}