Pith. sign in

REVIEW 2 major objections 6 minor 3 cited by

Black hole Green function factored into angular and radial parts

Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

T0 review · glm-5.2

2026-07-05 01:27 UTC pith:7CLQ47DG

load-bearing objection Paper extends scalar Green function regularization to spin-±1 and ±2 in Schwarzschild; the S² factorization is the main new result, but one step in its derivation is asserted rather than shown. the 2 major comments →

arxiv 2604.21218 v2 pith:7CLQ47DG submitted 2026-04-23 astro-ph.GA

Early metal-enriched baryon cycling before the midpoint of cosmic reionization

Yongda Zhu , Zhiyuan Ji , George D. Becker , Jiani Ding , Eiichi Egami , Xiaohui Fan , Xiangyu Jin , Weizhe Liu
show 8 more authors
Jianwei Lyu Zheng Ma Suprabhas Narisetty George H. Rieke Yunjing Wu Minghao Yue Junyu Zhang Marcia J. Rieke
This is my paper
classification astro-ph.GA PACS 04.70.-s04.30.-w04.25.Nj
keywords beforecosmicgalaxiesreionizationbaryoncyclingenrichmentionic
verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

This paper obtains exact analytic expressions for the direct (singular) part of the retarded Green function for the Bardeen-Press-Teukolsky equation on Schwarzschild spacetime, covering electromagnetic (spin ±1) and gravitational (spin ±2) perturbations. The key structural result is a separable factorization made possible by working in a spacetime conformal to Schwarzschild that decomposes as a direct product of a 2-dimensional Lorentzian spacetime M² and a 2-sphere S². The S² factor is obtained in closed form in terms of Euler angles, exploiting a connection between geodesics on the sphere, spin-weighted spherical harmonics, and rigid rotations. The M² factor splits into the square root of the van Vleck determinant (a standard geometric quantity measuring geodesic focusing) times a spin-dependent transport factor that the authors evaluate exactly via elliptic integrals for constant-radius orbits. With these factors in hand, the authors derive analytic multipolar (ℓ-mode) expressions for the direct part of the Green function. Subtracting these modes from the full Green function modes yields a non-direct remainder that provides a more accurate representation of the Green function near—but not exactly at—coincidence of the two spacetime points. This matters because practical calculations of quantities like the gravitational self-force on particles orbiting black holes require accurate Green function values in precisely this near-coincidence regime, where raw finite mode-sums are contaminated by the smearing of delta-distribution singularities.

Core claim

The central result is that the direct part of the Hadamard form of the Bardeen-Press-Teukolsky retarded Green function on conformal Schwarzschild spacetime admits a clean separable factorization Ũ = Ū·Ŭ, where the S² angular factor Ŭ = exp(−is(ᾱ+β̄))·(γ/sinγ)^{1/2} is given in closed form via Euler angles (Eq. 4.36), and the M² factor Ū = exp(2sλ)·Δ̄^{1/2} (Eq. 4.37) is expressible exactly through elliptic integrals for constant-radius worldlines. This factorization yields analytic ℓ-modes of the direct part (Eq. 5.20) for spins s = ±1 and ±2, which when subtracted from the full Green function ℓ-modes produce a non-direct part that better represents the Green function near coincidence than a

What carries the argument

The argument is carried by three structures: (1) the 2+2 conformal decomposition of Schwarzschild into M² × S², which enables variable separation in the Hadamard transport equation; (2) the Euler-angle parametrization of geodesics on S², which connects spin-weighted spherical harmonics to the angular factor via the addition theorem and yields the closed-form exp(−is(ᾱ+β̄)) prefactor; and (3) the (τ, ξ) coordinate system on M² — proper time and initial radial velocity along geodesics — in which the van Vleck determinant and the spin-dependent transport integral reduce to elliptic integrals for constant-radius orbits.

Load-bearing premise

The general-case formulas (beyond constant-radius worldlines) depend on an unpublished reference that provides the (τ, ξ) coordinate system on M² and the associated inverse Jacobian matrix elements. Additionally, the non-direct Green function calculation relies on ℓ-modes of the full Green function from a companion paper by two of the same authors. If either dependency contains errors, the central results propagate them.

What would settle it

Compute the ℓ-modes of the direct part numerically for a constant-radius orbit at r = 6M using an independent method (e.g., direct numerical integration of the transport equation in Schwarzschild coordinates without the 2+2 decomposition) and compare against Eq. 5.20. If the closed-form Euler-angle angular factor or the elliptic-integral M² factor produces ℓ-modes that disagree with the independent numerical calculation beyond expected discretization error, the factorization or the coordinate construction underlying it would be incorrect.

Watch this falsifier — get emailed when new claim-graph text bears on it.

If this is right

  • The analytic ℓ-modes of the direct part can be used to model the early-time (direct-wave) portion of black hole ringdown waveforms, complementing late-time quasinormal-mode fits.
  • The non-direct Green function calculation narrows the regime where separate Hadamard tail (V) expansions are needed, reducing the computational cost of self-force calculations for extreme-mass-ratio inspirals.
  • The separable factorization structure may extend to Kerr spacetime if an analogous conformal decomposition or approximate separation can be found, though Kerr lacks the direct-product structure exploited here.
  • The closed-form S² factor involving Euler angles and spin-weighted spherical harmonics reveals a geometric relationship that could inform calculations of spin-weighted fields on other spherically symmetric backgrounds.
  • The coincidence limit [Ṽ] = −s²·2M/r³ (Eq. 2.14) provides a direct check on numerical Green function codes for electromagnetic and gravitational perturbations.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit.

Referee Report

2 major / 6 minor

Summary. This manuscript derives a separable factorization of the direct (Hadamard) part of the retarded Green function for the Bardeen-Press-Teukolsky (BPT) equation on Schwarzschild spacetime, working in the conformally related M² × S² product spacetime. The S² factor is obtained in closed form in terms of Euler angles (Eq. 4.36), and the M² factor is expressed via elliptic integrals for constant-radius worldlines (Eqs. 4.49, 4.51, 4.53). These results yield analytic ℓ-modes of the direct part for spins s = ±1 and s = ±2 (Eq. 5.20), which are then subtracted from the full Green function ℓ-modes to produce a non-direct part that better represents the Green function near coincidence. The derivations are cross-validated through multiple independent methods: mollified integrals (App. A), elliptic integral representations (App. B), small-coordinate expansions (App. C), and near-coincidence expansions, all showing agreement (Fig. 2).

Significance. The paper provides genuinely new results: the closed-form S² factor involving Euler angles (Eq. 4.36) and the analytic ℓ-modes of the direct part for s ≠ 0 (Eq. 5.20) are not available in the prior literature. The cross-validation of four independent computational methods (Fig. 2, App. A–C) is a notable strength, as is the provision of a Mathematica notebook [29] for the small-coordinate expansions. The coincidence limit [V] = −s²(2M/r³) (Eq. 2.14) is derived via a covariant expansion (App. D) and is a useful independent check. The practical motivation—improving Green function representations for self-force calculations—is well articulated.

major comments (2)
  1. [§IV.A, between Eqs. (4.35) and (4.36)] The derivation of the S² factor Ŭ solves the transport equation (4.9) in a restricted coordinate system (θ' = π/2, φ' = 0, φ ∈ (0, π/2)) and obtains ζ_s = −is(ᾱ + β̄) in that restricted case (Eq. 4.35). The text then states: 'We now drop the restrictions on φ, γ, ᾱ and β̄ and use a direct calculation to show that this formula holds in all cases.' This direct calculation is not shown. The extension is non-trivial because spin-weighting breaks spherical symmetry (as the authors themselves note in Section II), so the WLOG argument that works for individual geodesics does not automatically apply to two-point functions with a chosen spin axis. Since Eq. (4.36) and all downstream results—the ℓ-modes (5.20) and the non-direct Green function (Figs. 8–9)—depend on this extension, the authors should either include the calculation or provide a rigorous argument for why the restricted-case result is
  2. [§IV.B, Eqs. (4.44) and (4.53)] The general-case formulas (beyond constant-radius worldlines) depend on reference [24] (Nolan, 'Coordinates on the geodesic spray and two-point functions on conformal Schwarzschild spacetime,' In preparation, 2026), which provides the (τ, ξ) coordinate system on M² and the inverse Jacobian matrix elements used in Eq. (4.44) and (4.53). Since this reference is unpublished, the general-case results cannot be independently verified by readers. The constant-radius case is self-contained (Eqs. 4.49, 4.51, 4.53 with App. A–B), but the claim in the abstract and conclusions that the expressions 'can be generalised to essentially all pairs of points' rests on [24]. The authors should clarify which results are fully self-contained in this paper versus which depend on [24], and ideally include the key coordinate construction in an appendix.
minor comments (6)
  1. [§IV.A, Eq. (4.9)] The separation constant is set to zero with the phrase 'where an arbitrary separation constant has been set equal to zero.' While the stress-test note correctly identifies that this is forced by the coincidence condition [Ū] = 1 (evaluating the M²-side equation at coincidence yields C = 0), the manuscript does not show this reasoning. A brief sentence explaining why C = 0 is forced would improve clarity.
  2. [§V.B, Figs. 8–9] The improvement from subtracting the direct part (pushing validity from Δt ≈ 4M to ≈ 1.6M in the circular case) is modest, as the authors acknowledge. The discussion of why the non-direct part still fails very near coincidence is clear, but it would help to state more precisely what practical gain this represents for self-force calculations—e.g., a quantitative estimate of the error reduction at intermediate Δt.
  3. [Figures 6–7] The plots show visual agreement between exact and approximate ℓ-modes, but no quantitative error analysis is provided for these figures (unlike Fig. 3 for the van Vleck determinant). A brief comment on the relative error would strengthen the validation.
  4. [§II.A, Eq. (2.19)] The statement that the spin-dependent prefactor 'increases or decreases exponentially with time depending on the sign of s(r−3M)' is correct but could note that this behavior is specific to circular geodesics and does not generalize to arbitrary worldlines, to avoid potential confusion.
  5. [References] Several references show placeholder URLs (e.g., [1], [3], [4], [7], [13], [26], [27]) rendered as repeated 'URL' strings. These should be corrected to actual URLs or DOIs.
  6. [Appendix C, Eq. (C8)] The small-coordinate expansion for Gd_ℓ is given to O(Δt⁶/M⁶) but the coefficients become quite unwieldy. The statement that higher-order coefficients are 'extremely unwieldy' is understandable, but it would help to note the convergence radius or the Δt range where this expansion is reliable.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for a careful reading and for identifying two substantive gaps in our presentation. Both points are well-taken, and we will address them in a revised manuscript.

read point-by-point responses

  1. Referee: The derivation of the S² factor Ŭ solves the transport equation in a restricted coordinate system and then claims the result holds generally without showing the calculation. The extension is non-trivial because spin-weighting breaks spherical symmetry.

    Authors: The referee is correct that the extension from the restricted case (θ'=π/2, φ'=0, φ∈(0,π/2)) to the general case is non-trivial and that the argument as currently stated is incomplete. We acknowledge that the standard 'without loss of generality' reasoning that applies to individual geodesics does not automatically extend to two-point functions with a fixed spin axis, precisely because spin-weighting breaks spherical symmetry, as we ourselves note in Section II. In the revised manuscript, we will include the direct calculation that verifies Eq. (4.35) for general (θ, φ, θ', φ'). The calculation proceeds by substituting the general-coordinate expressions for the Euler angles (4.14)–(4.15) and the world function ˚σ into the transport equation (4.22) and checking that ζ_s = −is(ᾱ + β̄) satisfies it identically. This is a straightforward but lengthy verification that involves differentiating the arctan expressions for ᾱ and β̄ with respect to θ and φ and confirming that the transport equation is satisfied. We will present this calculation, or at minimum a detailed sketch with the key steps, in the revised Section IV.A. We agree that without this, the derivation of Eq. (4.36) and all downstream results rests on an unjustified step. revision: yes

  2. Referee: The general-case formulas depend on unpublished reference [24] (Nolan, 'In preparation, 2026'), so the general-case results cannot be independently verified. The authors should clarify which results are self-contained versus which depend on [24], and ideally include the key coordinate construction in an appendix.

    Authors: The referee is correct that the general-case results (beyond constant-radius worldlines) depend on the (τ, ξ) coordinate construction and the inverse Jacobian matrix elements from [24], which is unpublished. We will address this in two ways in the revised manuscript. First, we will add a clear statement delineating which results are fully self-contained in this paper: the constant-radius case (Eqs. 4.49, 4.51, 4.53, together with Appendices A–C) is self-contained and does not depend on [24]. The general-case formula (4.44) and the claim that the expressions 'can be generalised to essentially all pairs of points' do depend on [24] for the coordinate construction and Jacobian elements. Second, we will include the key elements of the (τ, ξ) coordinate construction—specifically, the definition of the coordinates (Eqs. 4.40–4.41), the form of the metric (4.42), the expression for ˉσ (4.43), and the relevant inverse Jacobian matrix elements—in a new appendix. This will make the general-case results verifiable without requiring access to [24]. We note that the constant-radius results, which are the ones used in all numerical calculations and figures in the paper (Figs. 2–9), are already fully self-contained. revision: yes

Circularity Check

0 steps flagged

No significant circularity; self-citations are legitimate building blocks, not definitional restatements

full rationale

The paper derives new results (the S² factor Ŭ in Eq. 4.36, the M² factor Ū via elliptic integrals, and the ℓ-modes of the direct part in Eq. 5.20) that are not restatements of the cited prior work. The self-citations [8], [9], [13], and [24] serve as legitimate building blocks: [8] provides the scalar-case framework that the current paper generalizes to s≠0; [9] provides the 2+2 decomposition; [13] provides the full GF ℓ-modes that are subtracted from the newly derived direct-part ℓ-modes; and [24] (unpublished) provides the (τ,ξ) coordinate system used in the M² calculations. None of these citations define the target result in terms of itself. The separation constant being set to zero (Eq. 4.9) is mathematically forced by the coincidence condition [Ū]=1, as the reader's context correctly notes. The extension of the S² factor from the restricted equatorial case to general point pairs (the sentence after Eq. 4.35) is asserted via a 'direct calculation' that is not shown, but this is a gap in verification (a correctness/completeness concern), not a circularity: the result is not defined in terms of itself. The cross-checks between independent methods (elliptic integrals vs. small-coordinate expansions vs. near-coincidence expansions, shown in Fig. 2 and Appendix C) provide genuine internal validation. The dependence on unpublished [24] for the general (non-constant-radius) case is a verifiability concern, not a circularity concern. No step in the derivation chain reduces to its inputs by construction. The minor self-citations are standard in a multi-paper program and do not raise the circularity score above 2.

Axiom & Free-Parameter Ledger

0 free parameters · 5 axioms · 0 invented entities

No new physical entities, particles, forces, or dimensions are introduced. The Euler angles (ᾱ, γ, β̄) are standard mathematical constructs, not invented entities.

axioms (5)
  • standard math The Hadamard form (Eq. 2.7) provides a valid representation of the retarded GF in normal neighbourhoods, with the direct biscalary Ũ satisfying the transport equation (2.8).
    Invoked in Section II; this is a standard result from Friedlander [6] and Hadamard [5].
  • standard math Schwarzschild spacetime is conformal to M² × S² with the direct product metric (3.1), and the conformal Teukolsky operator (3.5) governs the rescaled field ψ̂ = rψ.
    Invoked in Section III; this is a standard conformal decomposition used in [8, 9].
  • domain assumption The (τ, ξ) coordinate system on M² and the inverse Jacobian matrix elements used in Eq. (4.44) are well-defined and correctly computed.
    Invoked in Section IV.B; depends on reference [24] (Nolan, in preparation, 2026), which is unpublished.
  • domain assumption The ℓ-modes of the full BPT GF for s=-2 from [13] (Aruquipa & Casals, 2026) are correct.
    Used in Section V.B to compute the non-direct part via subtraction (Eq. 5.10).
  • ad hoc to paper The separation constant in Eq. (4.7)-(4.9) is set to zero, allowing the factorization Ũ = Ū·Ŭ.
    Stated in Section IV: 'an arbitrary separation constant has been set equal to zero.' No independent justification is given beyond the coincidence limits [Ū]=[Ŭ]=1.

pith-pipeline@v1.1.0-glm · 28926 in / 4745 out tokens · 218198 ms · 2026-07-05T01:27:16.724620+00:00 · methodology

0 comments
read the original abstract

Models predict that chemical enrichment and gas redistribution should begin rapidly once star formation starts, but direct constraints at the earliest epochs have been scarce. Here we show that metal-enriched gas in multiple ionic phases was already present around galaxies before the midpoint of cosmic reionization. Using JWST/NIRSpec rest-frame ultraviolet spectroscopy from SPURS, we detect blueshifted metal absorption in three galaxies at $7.2<z<9.3$. The detected transitions span neutral, low-ionization, and high-ionization species, including O I, Si II, C II, Si IV, and C IV, with velocity offsets of $|\Delta v|\sim 50$--$250\,\mathrm{km\,s^{-1}}$ relative to nebular systemic redshifts. The ionic coexistence, overlapping velocity structure, and equivalent-width ratios are consistent with outflowing or otherwise kinematically disturbed galaxy-associated gas, implying rapid metal enrichment. These results show that key conditions for baryon cycling were established in at least a subset of luminous galaxies within the first several hundred million years of cosmic time, well before the completion of reionization.

Review history (2 revisions) →

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Forward citations

Cited by 3 Pith papers

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

  1. JWST absorption line spectroscopy with SPURS: ISM covering fractions and kinematics in individual galaxies at $z=5-9$

    astro-ph.GA 2026-06 unverdicted novelty 7.0

    JWST spectra of six z=5-9 galaxies show low-ionization covering fractions of 0.2-0.9 and diverse kinematics including blueshifted outflows, indicating heterogeneous multiphase ISM.

  2. Chemical signatures from the first stars embedded in metal-poor gas in galaxies at cosmic dawn

    astro-ph.GA 2026-06 unverdicted novelty 6.0

    JWST absorption spectra of galaxies at z~8-9 show metal-poor gas with high [C/O] suggesting enrichment by Population III supernovae.

  3. JWST Absorption-Line Analysis of UV-Bright Galaxies at $z=7.2-10.6$: Early Chemical Enrichment Traced by C, O, Mg, Al, Si, and Fe

    astro-ph.GA 2026-06 unverdicted novelty 5.0

    JWST spectra reveal that two z~7 galaxies already show near-solar iron-to-silicon ratios with no strong odd-even effect, favoring early Type Ia supernovae over pair-instability supernovae as the source of iron enrichment.

Reference graph

Works this paper leans on

32 extracted references · 32 canonical work pages · cited by 3 Pith papers · 1 internal anchor

  1. [1]

    Hence, Eq

    for general null geodesics of Schwarzschild spacetime. Hence, Eq. (2.20) yields �U(x, x ′) = � r′ r � �(3�+1)�2 � r′ −2M r−2M � �(1−�)�2 along radial null geodesics.(2.21) Further progress taking this 4-dimensional approach seems difficult, and as we have seen, the results above are limited by the assumption thatθ=θ ′ =π/2. Thus we take an alternative app...

    work page

  2. [2]

    There is an interesting corollary of the calculation above which we state as the following proposition

    An aside on Euler angles. There is an interesting corollary of the calculation above which we state as the following proposition. This relates geodesics and Euler angles onS 2. The result follows immediately from the fact thatκis constant along geodesics, and from the formula (4.33). ����������� �LetCbe a great circle ofS 2 and fixP ′ ∈C. LetPbe any other...

    work page

  3. [3]

    TwodLocalExpan- sions arxiv.nb

    fors= 0 (even if the values ofℓare also different). B. Calculation of the non-direct GF In the previous subsection we have calculated theℓ-modes −2Gd � of the direct part of the BPT GF. In [13], two of the authors calculated the BPT GF −2Gret and itsℓ-modes −2G�; we refer the reader to that paper for details of those calculations. Equipped with these, we ...

    work page

  4. [4]

    J. J. Duistermaat and L. H¨ ormander, Acta Mathematica128, 183 (1972), URL������������������������ ����������

    work page 1972

  5. [5]

    Poisson, A

    E. Poisson, A. Pound, and I. Vega, Living Reviews in Relativity14, 7 (2011)

    work page 2011

  6. [6]

    Casals, S

    M. Casals, S. Dolan, A. C. Ottewill, and B. Wardell, Phys. Rev. D88, 044022 (2013), URL�������������������� ������������������������������

    work page 2013

  7. [7]

    Wardell, C

    B. Wardell, C. R. Galley, A. Zengino˘ glu, M. Casals, S. R. Dolan, and A. C. Ottewill, Phys. Rev. D89, 084021 (2014), URL��������������������������������������������������

    work page 2014

  8. [8]

    Hadamard,Lectures on Cauchy’s Problem in Linear Partial Differential Equations(Dover Publications, 1923), ISBN 978-0486495491

    J. Hadamard,Lectures on Cauchy’s Problem in Linear Partial Differential Equations(Dover Publications, 1923), ISBN 978-0486495491

    work page 1923

  9. [9]

    F. G. Friedlander,The wave equation on a curved space-time, vol. 2 (Cambridge university press, 1975)

    work page 1975

  10. [10]

    R. H. Jonsson, D. Q. Aruquipa, M. Casals, A. Kempf, and E. Mart´ ın-Mart´ ınez, Phys. Rev. D101, 125005 (2020), URL����������������������������������������������������

    work page 2020

  11. [11]

    Casals, B

    M. Casals, B. C. Nolan, A. C. Ottewill, and B. Wardell, Physical Review D100, 104037 (2019)

    work page 2019

  12. [12]

    Casals and B

    M. Casals and B. C. Nolan, Physical Review D108, 044033 (2023)

    work page 2023

  13. [13]

    J. M. Bardeen and W. H. Press, Journal of Mathematical Physics14, 7 (1973), https://doi.org/10.1063/1.1666175, URL���������������������������������

    work page doi:10.1063/1.1666175 1973

  14. [14]

    S. A. Teukolsky, Phys. Rev. Lett.29, 1114 (1972)

    work page 1972

  15. [15]

    S. A. Teukolsky, Astrophys. J.185, 635 (1973)

    work page 1973

  16. [16]

    D. Q. Aruquipa and M. Casals,Green functions of the Regge-Wheeler and Teukolsky equations in Schwarzschild spacetime(2026), 2603.07747, URL��������������������������������

    work page arXiv 2026

  17. [17]

    Wald,General Relativity(Chicago, University of Chicago Press, 1984)

    R. Wald,General Relativity(Chicago, University of Chicago Press, 1984)

    work page 1984

  18. [18]

    Iuliano and J

    C. Iuliano and J. Zahn, Phys. Rev. D108, 125017 (2023), 2307.07467

    work page arXiv 2023

  19. [19]

    Goldberg, A

    J. Goldberg, A. Macfarlane, E. T. Newman, F. Rohrlich, and E. Sudarshan, Journal of Mathematical Physics8, 2155 (1967)

    work page 1967

  20. [20]

    T. J. Hollowood, G. M. Shore, and R. J. Stanley, Journal of High Energy Physics2009, 089 (2009)

    work page 2009

  21. [21]

    Casals and B

    M. Casals and B. C. Nolan, Physical Review D92, 104030 (2015)

    work page 2015

  22. [22]

    Casals and B

    M. Casals and B. C. Nolan, Physical Review D—Particles, Fields, Gravitation, and Cosmology86, 024038 (2012)

    work page 2012

  23. [23]

    D. A. Varshalovich, A. N. Moskalev, and V. K. Khersonskii,Quantum theory of angular momentum(World Scientific, 1988)

    work page 1988

  24. [24]

    Monteverdi and E

    A. Monteverdi and E. Winstanley, arXiv preprint arXiv:2410.23201 (2024)

    work page arXiv 2024

  25. [25]

    Michel and K

    V. Michel and K. Seibert, inMathematische Geod¨ asie/Mathematical Geodesy: Handbuch der Geod¨ asie, herausgegeben von Willi Freeden und Reiner Rummel(Springer, 2020), pp. 195–307

    work page 2020

  26. [26]

    Allen and T

    B. Allen and T. Jacobson, Communications in Mathematical Physics103, 669 (1986). 28

    work page 1986

  27. [27]

    B. C. Nolan,Coordinates on the geodesic spray and two-point functions on conformal Schwarzschild spacetime, In preparation. (2026)

    work page 2026

  28. [28]

    Pade Approximants of the Green Function in Spherically Symmetric Spacetimes

    M. Casals, S. Dolan, A. C. Ottewill, and B. Wardell, Phys. Rev.D79, 124044 (2009), 0903.5319

    work page internal anchor Pith review Pith/arXiv arXiv 2009

  29. [29]

    Oshita, S

    N. Oshita, S. Ma, Y. Chen, and H. Yang,Probing direct waves in black hole ringdowns(2025), 2509.09165, URL ��������������������������������

    work page arXiv 2025

  30. [30]

    J. Su, N. Khera, M. Casals, S. Ma, A. Chowdhuri, and H. Yang, arXiv e-prints arXiv:2601.22015 (2026), 2601.22015

    work page arXiv 2026

  31. [31]

    P. F. Byrd and M. D. Friedman,Handbook of elliptic integrals for engineers and physicists(Springer, 2013)

    work page 2013

  32. [32]

    Q. D. Aruquipa,�������������������������������������������

    work page