Recognition: unknown
Static Tidal Perturbations of Relativistic Stars: Corrected Center Expansion and Love Numbers-I
Pith reviewed 2026-05-10 10:17 UTC · model grok-4.3
The pith
A correction to the subleading term in the center expansion of static even-parity perturbations of relativistic stars does not change the extracted tidal Love number k2.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The regular-center Frobenius expansion of the interior even-parity master function has a corrected subleading coefficient that differs from the expression commonly used in the literature, yet numerical integrations for polytropic equations of state show that the correction affects only subleading initial data and leaves the extracted Love number k2 unchanged within numerical accuracy.
What carries the argument
The even-parity master function in Regge-Wheeler gauge and its regular Frobenius series expansion at the stellar center.
If this is right
- The standard quadrupolar equation can continue to be used for Love-number calculations once supplied with the corrected center expansion.
- Static even-parity perturbations admit a consistent formulation on Schwarzschild-de Sitter backgrounds that include a cosmological horizon.
- The full interior system in Regge-Wheeler gauge reduces exactly to the usual master equation employed in tidal-deformability studies.
- Polytropic models confirm that k2 is insensitive to the subleading center coefficient at current numerical precision.
Where Pith is reading between the lines
- Similar coefficient corrections could be derived and tested for higher multipoles to ensure consistency in gravitational-wave template modeling.
- For realistic neutron-star equations of state, the numerical insensitivity of k2 should be rechecked with independent codes at higher resolution.
- The Schwarzschild-de Sitter extension provides a controlled setting to study how a positive cosmological constant modifies tidal Love numbers.
- Gauge-invariant matching at the stellar surface might be examined to confirm that no hidden dependence on the center coefficient reappears.
Load-bearing premise
The standard quadrupolar even-parity master equation remains an accurate reduction of the full system, and surface boundary conditions do not reintroduce dependence on the corrected center coefficient into the final Love number.
What would settle it
A higher-precision numerical integration or analytic matching across the stellar surface that produces a statistically significant difference in the extracted k2 when the corrected center coefficient is used instead of the standard one.
read the original abstract
We revisit static tidal perturbations of relativistic stars with emphasis on two technical issues in the standard quadrupolar formulation. First, we derive the regular-center Frobenius expansion of the interior even-parity master function and obtain a corrected subleading coefficient, which differs from the expression commonly used in the literature. Second, we derive the static even-parity master equation on a Schwarzschild-de Sitter background, extending the usual asymptotically flat problem to a two-horizon geometry. To place these results on a common footing, we also show how the general interior even-parity system in Regge-Wheeler gauge reduces to the standard quadrupolar equation used in Love-number calculations. Numerical integrations for polytropic equations of state show that the corrected center coefficient affects only subleading initial data and leaves the extracted Love number $k_2$ unchanged within numerical accuracy. Taken together, these results fix the regular-center input to the standard quadrupolar problem and extend the static even-parity formalism to Schwarzschild-de Sitter backgrounds.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper revisits static tidal perturbations of relativistic stars, focusing on the regular-center Frobenius expansion of the interior even-parity master function, deriving a corrected subleading coefficient that differs from the literature. It also derives the static even-parity master equation on a Schwarzschild-de Sitter background and demonstrates the reduction of the general interior even-parity system in Regge-Wheeler gauge to the standard quadrupolar equation. Numerical integrations using polytropic equations of state show that this correction only affects subleading initial data and leaves the Love number k2 unchanged within numerical accuracy.
Significance. If the central results hold, this manuscript provides a valuable technical clarification to the standard formalism for computing tidal Love numbers of relativistic stars by correcting the center expansion coefficient and showing it has no effect on the extracted k2. The extension to Schwarzschild-de Sitter backgrounds broadens the applicability of the formalism. Strengths include the derivation of the reduction from the full system and the numerical verification that the correction is inconsequential for observables. This work helps ensure the accuracy of initial conditions in such calculations without changing physical predictions.
major comments (1)
- In the section deriving the reduction from the general interior even-parity system in Regge-Wheeler gauge to the standard quadrupolar master equation, the manuscript states that the reduction is shown. To confirm that the corrected subleading center coefficient remains invisible to the surface matching conditions (and thus to k2) in the unreduced system, it would be helpful to explicitly demonstrate that no additional coupling terms arise from the gauge choice or boundary conditions at the stellar surface when the full system is restored.
minor comments (2)
- The explicit algebraic form of the corrected subleading Frobenius coefficient and its direct comparison to the commonly used literature expression should be stated clearly in the derivation section to allow immediate verification by readers.
- The numerical methods section would benefit from additional details on grid resolution, surface boundary condition implementation, and convergence tests to support the claim that k2 is unchanged within numerical accuracy and to address reproducibility concerns.
Simulated Author's Rebuttal
We thank the referee for their positive evaluation of our manuscript and for the constructive suggestion that can further clarify the implications of our results. We address the major comment below and will incorporate the requested clarification in the revised version.
read point-by-point responses
-
Referee: In the section deriving the reduction from the general interior even-parity system in Regge-Wheeler gauge to the standard quadrupolar master equation, the manuscript states that the reduction is shown. To confirm that the corrected subleading center coefficient remains invisible to the surface matching conditions (and thus to k2) in the unreduced system, it would be helpful to explicitly demonstrate that no additional coupling terms arise from the gauge choice or boundary conditions at the stellar surface when the full system is restored.
Authors: We agree that an explicit demonstration at the surface would strengthen the presentation. In the revised manuscript we will add a concise paragraph (or short appendix) showing that the algebraic reduction from the full Regge-Wheeler-gauge system to the master equation remains valid up to the stellar surface. Because the reduction follows directly from the linearized Einstein equations without additional assumptions, the metric perturbations and their derivatives that enter the surface matching conditions map one-to-one between the two formulations; no extra coupling terms are generated by the gauge choice or by the boundary conditions themselves. Consequently, the regular-center Frobenius series for the master function translates directly into the appropriate initial data for the full system, and the corrected subleading coefficient continues to affect only terms that vanish at the surface values relevant for exterior matching. This equivalence guarantees that k2 extracted from the unreduced system is identical to the value obtained from the master equation, consistent with our existing numerical results. revision: yes
Circularity Check
No significant circularity; derivations are self-contained from linearized Einstein equations
full rationale
The paper begins from the standard linearized Einstein equations in Regge-Wheeler gauge for static even-parity perturbations, derives the regular-center Frobenius series for the interior master function by direct substitution into the master equation, obtains the corrected subleading coefficient analytically, demonstrates the reduction of the general interior system to the quadrupolar master equation, and extends the master equation to the Schwarzschild-de Sitter background. The numerical claim that the corrected coefficient leaves k2 unchanged is obtained by integrating the resulting standard equation with the new initial data; no fitted parameters are renamed as predictions, no load-bearing uniqueness theorems or ansatze are imported via self-citation, and none of the central results reduce to their own inputs by construction. The derivation chain is therefore independent and externally falsifiable against the underlying field equations.
Axiom & Free-Parameter Ledger
axioms (2)
- domain assumption Linearized Einstein equations hold for small static even-parity perturbations on a spherically symmetric background
- domain assumption The stellar interior is described by a perfect fluid with a barotropic equation of state
Forward citations
Cited by 1 Pith paper
-
Stationary solutions in the small-$c$ expansion of GR
Exact Lense-Thirring-type, C-metric-type, and Hartle-Thorne-type stationary vacuum solutions are constructed in the NLO and NNLO small-c expansion of GR, revealing a richer sector than magnetic Carroll gravity.
Reference graph
Works this paper leans on
-
[1]
Odd-parity equations Settingh tr = 0, thetθcomponent of the linearized Ricci tensor becomes 2δRtθ = 1 r2 sin2 θ ∂ϕ∂θhtϕ + −f ∂2 r − 4M r3 − 1 r2 sin2 θ ∂2 ϕ htθ.(63) We now substitute the odd-parity harmonic expansions. Since the perturbations are static, the radial amplitudes depend only onr: htθ = X lm −hlm t (r) 1 sinθ ∂ϕYlm(θ, ϕ),(64) htϕ = X lm hlm t...
-
[2]
Even-parity equations We next turn to the even-parity sector of the static exterior perturbation. In the vacuum regionr > R, the field equations reduce to δRµν = 0.(75) In addition to the constraints (62), we impose the angular identity hϕϕ = sin2 θ hθθ.(76) Starting from the general static expression forδRtt, given in Appendix B1, and using the Schwarzsc...
-
[3]
Thetθcomponent of Eq
Odd-parity sector We begin with the axial (odd-parity) sector. Thetθcomponent of Eq. (108), together with¯gtθ = 0andδT tθ = −¯ρ htθ, yields δRtθ + 4π(3¯p+ ¯ρ)htθ = 0.(118) For the relevant components ofδTµν, see Appendix H. Using the coordinate expression forδRtθ given in Appendix I, together with the odd-parity harmonic expansions htθ = X lm −hlm t (r) 1...
-
[4]
Thettcomponent of Eq
Even-parity sector We now turn to the polar (even-parity) sector. Thettcomponent of Eq. (108) reads δRtt + 1 2 e2ψδR−4π(¯ρ−3¯p)htt = 8π e2ψδρ−¯ρhtt ,(127) which may be rearranged as δRtt + 1 2 e2ψδR+ 4π(3¯p+ ¯ρ)htt −8πe 2ψδρ= 0.(128) To eliminateδR, we take the trace of Eq. (107) and use the background field equations, obtaining δR= 8π hµν ¯Tµν −¯gµνδTµν ...
-
[5]
Tidal love numbers of neutron stars.Astrophysical Journal, 677(2):1216–1220, 2008
Tanja Hinderer. Tidal love numbers of neutron stars.Astrophysical Journal, 677(2):1216–1220, 2008
2008
-
[6]
Flanagan and Tanja Hinderer
Éanna E. Flanagan and Tanja Hinderer. Constraining neutron-star tidal love numbers with gravitational-wave detectors. Physical Review D, 77(2):021502, 2008
2008
-
[7]
Relativistic theory of tidal love numbers.Physical Review D, 80(8):084018, 2009
Taylor Binnington and Eric Poisson. Relativistic theory of tidal love numbers.Physical Review D, 80(8):084018, 2009
2009
-
[8]
Relativistic tidal properties of neutron stars.Physical Review D, 80(8):084035, 2009
Thibault Damour and Alessandro Nagar. Relativistic tidal properties of neutron stars.Physical Review D, 80(8):084035, 2009
2009
-
[9]
Tidal Response of Compact Objects
Sumanta Chakraborty and Paolo Pani. Tidal Response of Compact Objects.arXiv preprint arXiv:2604.08679, 4 2026
work page internal anchor Pith review Pith/arXiv arXiv 2026
-
[10]
Love numbers of black holes and compact objects
María J. Rodríguez, Luca Santoni, and Adam R. Solomon. Love numbers of black holes and compact objects.arXiv preprint arXiv:2604.08653, 4 2026
work page internal anchor Pith review Pith/arXiv arXiv 2026
-
[11]
Frank J. Zerilli. Effective potential for even-parity regge-wheeler gravitational perturbation equations.Phys. Rev. Lett., 24:737–738, Mar 1970
1970
-
[12]
Magnetic tidal love numbers clarified
Paolo Pani, Leonardo Gualtieri, Tiziano Abdelsalhin, and Xisco Jiménez-Forteza. Magnetic tidal love numbers clarified. Physical Review D, 98(12):124023, 2018
2018
-
[13]
Gw170817: Observation of gravitational waves from a binary neutron star inspiral.Physical Review Letters, 119(16):161101, 2017
LIGO Scientific Collaboration and Virgo Collaboration. Gw170817: Observation of gravitational waves from a binary neutron star inspiral.Physical Review Letters, 119(16):161101, 2017
2017
-
[14]
Properties of the binary neutron star merger gw170817.Physical Review X, 9(1):011001, 2019
LIGO Scientific Collaboration and Virgo Collaboration. Properties of the binary neutron star merger gw170817.Physical Review X, 9(1):011001, 2019
2019
-
[15]
Tullio Regge and John A. Wheeler. Stability of a schwarzschild singularity.Physical Review, 108(4):1063–1069, 1957
1957
-
[16]
Non-radial pulsation of general-relativistic stellar models
Kip S Thorne and Alfonso Campolattaro. Non-radial pulsation of general-relativistic stellar models. i. analytic analysis for l>= 2.Astrophysical Journal, vol. 149, p. 591, 149:591, 1967
1967
-
[17]
Tidal deformability of neutron stars in unimodular gravity.Universe, 8(11):576, 2022
Rui-Xin Yang, Fei Xie, and Dao-Jun Liu. Tidal deformability of neutron stars in unimodular gravity.Universe, 8(11):576, 2022
2022
-
[18]
Thorne and James B
Kip S. Thorne and James B. Hartle. Laws of motion and precession for black holes and other bodies.Phys. Rev. D, 31:1815–1837, Apr 1985
1985
-
[19]
Kip S. Thorne. Tidal stabilization of rigidly rotating, fully relativistic neutron stars.Physical Review D, 58(12), November 1998
1998
-
[20]
Robert F. Tooper. General Relativistic Polytropic Fluid Spheres.Astrophys. J., 140:434, August 1964
1964
-
[21]
Robert F. Tooper. Adiabatic Fluid Spheres in General Relativity.Astrophys. J., 142:1541, November 1965
1965
-
[22]
Courier Corporation, 2013
Richard P Brent.Algorithms for minimization without derivatives. Courier Corporation, 2013
2013
-
[23]
On some static solutions of Einstein’s gravitational field equations in a spherically symmetric case.Sci
Hidekazu Nariai. On some static solutions of Einstein’s gravitational field equations in a spherically symmetric case.Sci. Rep. Tohoku University Eighth Ser., 34:160, 1950
1950
-
[24]
On a New Cosmological Solution of Einstein’s Field Equations of Gravitation.Gen
Hidekazu Nariai. On a New Cosmological Solution of Einstein’s Field Equations of Gravitation.Gen. Rel. Grav., 31(6):963– 971, 1999
1999
-
[25]
Static Tidal Perturbations of Schwarzschild–de Sitter Spacetime Near the Nariai Limit-II.to appear
Emel Altas, Onur Oktay, Ercan Kilicarslan, and Bayram Tekin. Static Tidal Perturbations of Schwarzschild–de Sitter Spacetime Near the Nariai Limit-II.to appear. Appendix A: Background geometry: connection and curvature components This appendix summarizes the Levi-Civita connection and curvature for the background metric (9). 27
-
[26]
Nonzero Christoffel symbols For the line element (9), the non-vanishing connection coefficients are: ¯Γt tr =ψ ′ , ¯Γr rr =− f ′ 2f , ¯Γr tt =fe 2ψψ′ ,(A1) ¯Γθ rθ = 1 r , ¯Γϕ rϕ = 1 r , ¯Γr θθ =−rf ,(A2) ¯Γr ϕϕ =−rfsin 2 θ , ¯Γθ ϕϕ =−sinθcosθ , ¯Γϕ θϕ = cotθ .(A3)
-
[27]
Ricci tensor and scalar The independent nonzero Ricci components can be written as ¯Rtt = e2ψf ψ′′ + (ψ′)2 + ψ′f ′ 2f + 2ψ′ r ,(A4) ¯Rrr =−ψ ′′ −(ψ ′)2 − ψ′f ′ 2f − f ′ rf ,(A5) ¯Rθθ =f − rf ′ 2f −rψ ′ −1 + 1, ¯Rϕϕ = sin2 θ ¯Rθθ .(A6) The scalar curvature ¯R= ¯gµν ¯Rµν reduces to ¯R=−2f(ψ ′)2 −2f ψ ′′ −f ′ψ′ − 4f r ψ′ − 2f ′ r − 2f r2 + 2 r2 .(A7) Appendi...
-
[28]
Diagonal components a.ttcomponent 2δRtt =f 2e2ψ −ψ′∂r − 2f ′ f ψ′ −2(ψ ′)2 −2ψ ′′ − 4 r ψ′ hrr + −f ∂2 r + 2f ψ′ − f ′ 2 − 2f r ∂r −2f(ψ ′)2 − 1 r2 D2 htt + 2f e2ψ r2 ψ′ (−∂θ −cotθ)h rθ − 2f e2ψ r2 sin2 θ ψ′ ∂ϕhrϕ + f e2ψ r2 ψ′ ∂r − 2 r hθθ + f e2ψ r2 sin2 θ ψ′ ∂r − 2 r hϕϕ.(B1) 28 b.rrcomponent 2δRrr = f ψ′ + 2f r ∂r +ψ ′f ′ + 2f ′ r − 1 r2 D2 hrr +e −2ψ...
-
[29]
Off-diagonal components a.tθcomponent 2δRtθ = f ∂r∂θ + f ψ′ + f ′ 2 ∂θ htr + 1 r2 sin2 θ ∂ϕ∂θhtϕ + −f ∂2 r + f ψ′ − f ′ 2 ∂r − 4f r ψ′ − 1 r2 sin2 θ ∂2 ϕ htθ.(B5) b.tϕcomponent 2δRtϕ = f ∂r∂ϕ + f ψ′ + f ′ 2 ∂ϕ htr + 1 r2 ∂θ∂ϕ − 1 r2 cotθ ∂ ϕ htθ + −f ∂2 r + f ψ′ − f ′ 2 ∂r − 4f r ψ′ − 1 r2 ∂2 θ + 1 r2 cotθ ∂ θ htϕ.(B6) c.rθcomponent 2δRrθ = −2f r ψ′ − f ′...
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.