REVIEW 3 major objections 5 minor 31 references
The linearized Einstein equations in Kerr can be directly separated, yielding a unique set of decoupled mode functions for the spin-2 perturbation, with all ten metric components expressed explicitly in terms of the Weyl scalars f0 and f4.
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 →
T0 review · deepseek-v4-flash
2026-08-01 20:32 UTC pith:4T6JMT3M
load-bearing objection A serious and readable attack on a decades-old problem, but the claimed uniqueness rests on an admitted ansatz and a shaky trace argument; worth refereeing, not yet a settled result. the 3 major comments →
Separating the linearized Einstein equations in Kerr
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The central claim is that the linearized Einstein equations in a Kerr background can be separated directly, without the usual detour through a Hertz potential or radiation gauge. The paper shows that the trace h of the metric perturbation satisfies the spin-0 wave equation (13), so the spin-2 sector is traceless; imposing that condition together with the de Donder gauge (2), the Killing-Yano eigen-equation (7), the parity condition (17), and the definitions of the perturbed Weyl scalars f0 and f4 reduces the coupled system first to 26 and then to a single 12th-order equation (29) for the mode function fxx. Equation (29) is solved by the rational ansatz (31), producing explicit formulas (33)
What carries the argument
The central machinery is the rational ansatz (31), fxx = (HY)^(-n) sum_i [p_i(x) f_i + q_i(x) f'_i], with f_i running over f0, f4 and their parity conjugates, H = r^2 + a^2 x^2, Y = 1 - x^2, and n tested up to 6. This ansatz is what makes the last reduced equation (29) algebraically solvable: the unknown polynomial coefficients p_i, q_i are determined by linear equations after lowering derivatives with (22), and once fxx is known the other nine components follow sequentially. The matching companion object is the Teukolsky-Starobinsky-type identity (37), whose constant Q fixes the 'unique' relation between f0 and f4 and is identified with the Teukolsky-Starobinsky constant.
Load-bearing premise
The load-bearing premise is that the rational ansatz (31) — fxx written as a finite combination of f0, f4, their parity conjugates and first derivatives over (HY)^n — covers every solution of the reduced equation (29); the paper itself concedes that some generality may have been lost, and if a valid solution lies outside this rational family the claimed uniqueness collapses.
What would settle it
Take the Schwarzschild limit a -> 0 in the explicit formulas (33)/(39)-(41): the components must reduce (up to pure gauge) to the standard odd- and even-parity metric perturbation for Schwarzschild. A term-by-term mismatch would falsify the claim that (31) captures the full spin-2 solution space. Alternatively, solve (29) numerically for fixed M, a, omega, m, lambda without imposing the rational ansatz: finding any solution not of the form (31) would refute the uniqueness.
If this is right
- The full linearized metric around Kerr is obtained directly from f0/f4, eliminating the Hertz-potential step and the radiation-gauge singularities that complicate metric reconstruction.
- The two independent solution branches found in earlier symmetry-based constructions are recovered as a byproduct of the derivation, now derived rather than assumed, so the mode-function space has a first-principles derivation.
- The traceless condition (14) follows from the spin-0 equation (13) plus the de Donder gauge, so spin-0 and spin-2 sectors decouple; spin-2 perturbations carry zero trace.
- For Q ≠ 0, the relation (37) forces ψ0 = 0 to imply ψ4 = 0 (and vice versa), matching the expectation for physically regular perturbations; algebraically special modes are excluded by fixed parity when Q ≠ 0.
- Because the master functions satisfy separable ordinary differential equations, the resulting mode functions are amenable to standard Teukolsky-mode numerical routines; this sets up direct metric assembly for gravitational-wave templates.
Where Pith is reading between the lines
- I would test completeness head-on: search for solutions of (29) with a denominator more general than (HY)^n. If any exists, the 'unique' set of mode functions is unique only inside the ansatz family, a weaker statement than the abstract suggests.
- The h = 0 step is presented as following from (13), but (13) only says the trace obeys a scalar wave equation; treating the trace as pure spin-0 and dropping it may be a gauge-fixing choice. A worthwhile test is whether a residual gauge transformation can reinstate a nonzero trace without changing physical observables.
- The explicit formulas suggest a practical numerical pipeline: solve the ODEs for f4, use (37) for f0, and evaluate (40)-(41). Comparing the result in the Schwarzschild (a -> 0) limit against the standard odd- and even-parity metric perturbation would be a cheap, decisive validation and could expose missing pure-gauge terms.
- If the ansatz completeness is confirmed, the method likely extends to other Petrov-D backgrounds with a Killing-Yano tensor; the same reduction logic might separate metric perturbations in that broader class of spacetimes.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper claims a direct separation of the linearized Einstein equations (LEEs) in Kerr spacetime. Using the de Donder gauge, the Killing-Yano symmetry operator, a fixed-parity condition, the Weyl-scalar constraints, and the tracelessness condition h=0, the author reduces the ten metric-perturbation components to explicit expressions in terms of Teukolsky master functions f0, f4 and their parity conjugates. The central results are Eqs. (33), (40), and (41), together with the Teukolsky-Starobinsky-type relation (37). The paper argues that the resulting mode functions are unique and that they reproduce the two solutions found earlier in [23].
Significance. If the derivation is valid, this would be a substantial technical step: explicit metric perturbations in Kerr expressed directly through the separable Teukolsky functions would be useful for self-force calculations, gravitational-wave modeling, and studies of second-order perturbations. The paper is also commendably explicit about the assumptions it makes, and it provides a concrete consistency check by recovering the earlier solutions of [23] and the standard Teukolsky-Starobinsky constant. However, the advertised uniqueness and the claim that the trace is irrelevant are not established by the argument as written. The result is best read, at this stage, as the construction of a particular traceless, parity-fixed, ansatz-restricted family of separated mode functions, not as a proof of the uniqueness of separated modes.
major comments (3)
- [Section II, Eqs. (13)-(14)] The step from g^{μν} e2[hμν] = -1/2 ∇^μ∇_μ h = 0 to h=0 is a non sequitur. Any solution of the scalar wave equation satisfies the same equation; the trace is not forced to vanish. Moreover, under the de Donder gauge there remains residual gauge freedom, and h can be changed by such transformations, so h=0 is an additional gauge-type restriction rather than a consequence of the field equations. This restriction is then used algebraically in Eq. (25) to eliminate ftϕ. The derivation therefore covers only the traceless sector. The abstract's statement that the trace 'is irrelevant' and the conclusion h=0 in Eq. (14) need to be replaced either by a proof that h can be set to zero by residual gauge freedom for all perturbations under consideration, or by an explicit statement that tracelessness is being imposed as an assumption and that the final claims are restricted accordingly.
- [Section III.B.4, Eqs. (29)-(31)] The uniqueness claim rests on an unproven ansatz. Equation (29) is a single linear ODE containing derivatives of fxx up to order 12. It is solved by postulating the rational form (31), fxx = (HY)^{-n} Σ_i [p_i(x)f_i + q_i(x)f_i'], with n chosen by trial up to 6 and p_i, q_i finite-degree polynomials. Substituting this ansatz and matching coefficients yields linear equations and a unique solution within the ansatz family, but nothing in the derivation shows that every solution of (29) lies in this family, nor that homogeneous solutions of (29) are absent or pure gauge. The paper itself concedes in Section V that 'some generality might have been lost.' A second completeness gap appears when the r-dependence and parameter dependence are restored by assuming polynomial coefficient structures and matching into (26b). Consequently the word 'unique' in the abstract and in Section IV is not supp
- [Section IV, Eqs. (33)-(41)] The formulas are not independently checkable from the main text: all polynomials Ai, Bi, Ci, Di, Ei are relegated to separate supplemental files, and the derivation of Eq. (29) itself is described only schematically. More importantly, the input data include the Teukolsky equations (22)-(23), the Killing-Yano eigen-equation (7), and the parity-conjugated constraints (18), and the final mode functions are found to be identical to those of [23]. This does not by itself invalidate the calculation, but it changes the status of the claimed 'direct separation': the paper should state clearly which parts of the result are derived from the LEEs and which are assumed through the Teukolsky/KY input, and it should discuss whether the construction is independent of [23] or a rederivation of it within a more systematic scheme.
minor comments (5)
- [Section III.B.4] The ansatz (31) is motivated by the structure of the Kerr metric, but the statement that 'by experimenting with n as large as 6, a unique solution can be found' needs more detail: how was n chosen, how was uniqueness within the ansatz verified, and were higher n checked to ensure that no additional solutions appear?
- [Eq. (25)] The tracelessness condition is solved for ftϕ with a denominator M a r. The special cases M=0, a=0, or r=0 are not discussed; at least a comment on the Schwarzschild and static limits is needed.
- [Eqs. (40)-(41)] The notation f^{(4)}_{μν} and f^{(0)}_{μν} is potentially confusing, since 0 and 4 also label the two Teukolsky functions. A brief explanation of the superscript convention would improve readability.
- [Section V] The paragraph on algebraically special solutions with Q=0 is interesting but appears only as a discussion item. If Q=0 is indeed a necessary condition for such solutions in this construction, that claim should be stated more precisely and, ideally, demonstrated within the derivation.
- [General] There are several typographical issues in the schematic equations, e.g., the stray comma in Eq. (26a) and the inconsistent use of f••, f•, f∗∗ in Eqs. (26)-(29). These should be cleaned up before publication.
Circularity Check
No significant circularity: the ansatz completeness caveat is a limitation, not a circular reduction; K4 from [23] is explicit and independently checkable.
full rationale
The claimed derivation chain is not circular. Section II introduces the KY operator K4 explicitly (Eq. 6), so the load-bearing algebraic object is not imported as a black box; the commutation claim is an independently checkable mathematical assertion. The paper uses the Teukolsky definitions (10), parity, tracelessness, and gauge constraints as inputs, reduces the system to Eq. (29), and solves it by the explicit rational ansatz (31). This is a standard solution method, not a fit to a target: the polynomial coefficients in (31) are determined by requiring the residual equation (29) to vanish, and the final step checks the previously unused equations (26a), leading to the nontrivial f0-f4 relation (35)/(37). The later statement that the resulting f(0), f(4) are "identical to the two solutions found in [23]" is a post hoc cross-check, not an input to the derivation. The main caveat is a completeness gap, not circularity: Section V concedes "some generality might have been lost" because (29) was solved by an ansatz. That undermines the abstract's unqualified "unique set," but it is a correctness/justification limitation, not a reduction of the output to the input by construction. The self-citation [23] is not load-bearing in a circular way; the operator and equations are stated in the present paper, and the solution is not assumed from [23]. Hence no significant circularity.
Axiom & Free-Parameter Ledger
free parameters (2)
- ansatz denominator power n =
positive integer, found by experiment as large as 6
- degree bounds for the ansatz polynomials p_i(x), q_i(x) =
not stated in text
axioms (7)
- domain assumption Teukolsky's master equation describes the perturbed Weyl scalars ψ0 and ψ4 and is separable in Kerr.
- domain assumption The symmetry operator K4 commutes with the linearized Einstein operator e2 off-shell, and all solutions can be classified by eigenmodes K4[h_μν]=λ h_μν.
- ad hoc to paper The trace of the perturbed metric is set to zero: h=0.
- domain assumption Mode functions have fixed parity under x→−x, Eq. (17).
- ad hoc to paper The final unknown mode function fxx admits the rational ansatz (31): fxx = 1/(HY)^n Σ_i [p_i(x) f_i + q_i(x) f'_i] with finite-order polynomials.
- domain assumption The (t,ϕ) dependence factors as e^{-i(wt-mϕ)}.
- domain assumption The de Donder gauge condition (2) is imposed and does not restrict generality.
read the original abstract
A direct separation of the linearized Einstein equations in Kerr is presented. The trace of the perturbed metric is found to obey the spin-0 perturbation equation and so is irrelevant for the spin-2 perturbations. By combining the traceless condition, the Killing-Yano symmetry, a parity requirement, and the de Donder gauge condition, it has been found possible to reduce the linearized Einstein equations in Kerr to a manageable level, thereby enabling the direct separation of the equations and the derivation of a unique set of decoupled mode functions for the spin-2 perturbation of the Kerr black hole.
Reference graph
Works this paper leans on
-
[1]
On a computer, this can be achieved by iterating the DCS procedure without much difficulty
Reducing the order ofr-derivatives To reduce the number of equations, one can firstly try to lower the order of derivatives with respect to one of the coordinates, sayr. On a computer, this can be achieved by iterating the DCS procedure without much difficulty. The results can be written schematically as ˙fµν =F (10) µν h f , f′ , f′′ , f′′′ i ,(20a) 0 =E...
-
[2]
At this stage, however, each it- eration of the DCS procedure often takes a too long time to finish
Reducing the order ofx-derivatives The equations in (20b) can be further reduced by low- ering the order of derivatives with respect tox, again by 4 using the DCS procedure. At this stage, however, each it- eration of the DCS procedure often takes a too long time to finish. To accelerate the calculation, one can give nu- merical values to the physical par...
-
[3]
Applying the traceless condition The traceless condition (14) is an algebraic equation of the mode functions and can be solved as ftϕ =− 1 4M ar nh (r2 +a 2)2 −a 2XY i ftt −X 2frr −XY fxx − X Y −a 2 fϕϕ o .(25) This eliminates one mode function from all the equations
-
[4]
Reducing the order ofr-derivatives For the 28 partial differential equations in (9), (12) and (18), one can repeat what have been done for (9) and (12) in the previous subsection and reduce the order ofr-derivatives by iterating the DCS procedure. In this way, the 28 equations can be reduced to ˙fµν =F (9) µν h f••, f′ ••, f′′ ••, f′′′ •• i , ,(26a) 0 =E ...
-
[5]
So, there are six (three) equations havingx-derivatives up to the first (second) order
Reducing the order ofx-derivatives By giving numerical values to the physical parame- ters, one can further reduce (26b) to 13 equations, among which nine are of the form: f ′ µν =F (6) µν h f••, f′ ••, f•, f′ • i ,(27a) f ′′ µν =F (3) µν h f••, f′ ••, f•, f′ • i ,(27b) where for (27a),f µν ∈ {ftt,f tx,f rx,f rϕ,f xϕ,f ϕϕ}, and for (27b),f µν ∈ {ftr,f rr,...
-
[6]
Firstly, only three mode functions in (27) have second- 5 order derivatives, while the remaining 6 only have first- order derivatives
Solving for the mode functions There are significant improvements of (27) over (21). Firstly, only three mode functions in (27) have second- 5 order derivatives, while the remaining 6 only have first- order derivatives. Secondly, the overall expressions of (27) are much simpler than those of (21). Consequently, all mode functions can now be solved sequent...
-
[7]
C. Shi, Q. Zhang, and J. Mei, Phys. Rev. D110, 124007 (2024), arXiv:2407.13110 [gr-qc]
Pith/arXiv arXiv 2024
-
[8]
J. Luoet al., Living Rev. Rel.29, 1 (2026), arXiv:2502.20138 [gr-qc]
Pith/arXiv arXiv 2026
-
[9]
K. G. Arunet al.(LISA), Living Rev. Rel.25, 4 (2022), arXiv:2205.01597 [gr-qc]
Pith/arXiv arXiv 2022
-
[10]
Y. Yang, C. Shi, and Y.-M. Hu, (2025), arXiv:2510.16903 [gr-qc]
arXiv 2025
-
[11]
Y.-F. Wang, S. Ma, N. Khera, and H. Yang, (2026), arXiv:2601.05734 [gr-qc]
arXiv 2026
-
[12]
E. Poisson, A. Pound, and I. Vega, Living Rev. Rel.14, 7 (2011), arXiv:1102.0529 [gr-qc]
Pith/arXiv arXiv 2011
-
[13]
L. Barack and A. Pound, Rept. Prog. Phys.82, 016904 (2019), arXiv:1805.10385 [gr-qc]
Pith/arXiv arXiv 2019
-
[14]
A. Pound and B. Wardell, (2021), 10.1007/978-981-15- 4702-7 38-1, arXiv:2101.04592 [gr-qc]
Pith/arXiv arXiv 2021
- [15]
-
[16]
S. Hollands and V. Toomani, Class. Quant. Grav.43, 055001 (2026), arXiv:2405.18604 [gr-qc]
Pith/arXiv arXiv 2026
-
[17]
S. A. Teukolsky, Phys. Rev. Lett.29, 1114 (1972)
1972
-
[18]
S. A. Teukolsky, Astrophys. J.185, 635 (1973)
1973
-
[19]
P. L. Chrzanowski, Phys. Rev.D11, 2042 (1975)
2042
-
[20]
L. S. Kegeles and J. M. Cohen, Phys. Rev.D19, 1641 (1979)
1979
-
[21]
R. M. Wald, Phys. Rev. Lett.41, 203 (1978)
1978
-
[22]
J. M. Stewart, Proc. Roy. Soc. Lond.A367, 527 (1979)
1979
-
[23]
A. K.-W. Chung, P. Wagle, and N. Yunes, Phys. Rev. D 107, 124032 (2023), arXiv:2302.11624 [gr-qc]
Pith/arXiv arXiv 2023
-
[24]
A. K.-W. Chung, P. Wagle, and N. Yunes, Phys. Rev. D 109, 044072 (2024), arXiv:2312.08435 [gr-qc]
Pith/arXiv arXiv 2024
-
[25]
A. K.-W. Chung and N. Yunes, Phys. Rev. D110, 064019 (2024), arXiv:2406.11986 [gr-qc]
Pith/arXiv arXiv 2024
-
[26]
A. K.-W. Chung and N. Yunes, Phys. Rev. Lett.133, 181401 (2024), arXiv:2405.12280 [gr-qc]
Pith/arXiv arXiv 2024
-
[27]
B. Chen and L. C. Stein, Phys. Rev. D96, 064017 (2017), arXiv:1707.05319 [gr-qc]
Pith/arXiv arXiv 2017
-
[28]
N. Franchini, Phys. Rev. D108, 044079 (2023), arXiv:2305.19313 [gr-qc]
Pith/arXiv arXiv 2023
-
[29]
J. Mei, Sci. China Phys. Mech. Astron.69, 270415 (2026), arXiv:2311.18409 [gr-qc]
Pith/arXiv arXiv 2026
-
[30]
R. M. Wald, J. Math. Phys.14, 1453 (1973)
1973
-
[31]
Chandrasekhar, Proc
S. Chandrasekhar, Proc. Roy. Soc. Lond. A392, 1 (1984)
1984
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.