arxiv: 2605.03305 · v1 · submitted 2026-05-05 · 🌀 gr-qc · astro-ph.HE· hep-th· nucl-th

Recognition: unknown

Non-radial pulsations of gravitationally coupled two-fluid neutron stars in general relativity

Ankit Kumar , Daniel A. Caballero , Hajime Sotani , Nicol\'as Yunes

Authors on Pith no claims yet

Pith reviewed 2026-05-07 14:41 UTC · model grok-4.3

classification 🌀 gr-qc astro-ph.HEhep-thnucl-th

keywords two-fluid neutron starsnon-radial pulsationspolar perturbationsgeneral relativityasteroseismologygravitational wavesmode spectraneutron star oscillations

0 comments

The pith

A fully relativistic framework now exists to compute polar pulsations of two-fluid neutron stars coupled only by gravity.

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

The paper develops the complete set of linear perturbation equations that describe non-radial polar oscillations in neutron stars composed of two fluids whose only interaction is the shared spacetime in general relativity. It derives the coupled equations for the metric and the two fluid variables, then supplies the regularity conditions at the center, the surface continuity requirements, and the exterior vacuum matching conditions needed to turn the system into a well-posed eigenvalue problem. The authors implement the resulting equations numerically and obtain example spectra for representative two-fluid stellar models. If the construction holds, the computed eigenfunctions reveal which fluid component dominates each mode through their radial node structure. This matters because multi-fluid interiors are expected in many compact objects, and their distinct oscillation signatures could be detectable in gravitational waves.

Core claim

We develop a fully relativistic framework for polar perturbations of gravitationally coupled two-fluid neutron stars, assuming that the two fluids interact only through the common spacetime and are not coupled by entrainment or direct microphysical interactions. We derive the coupled linear perturbation equations governing the metric and both fluid components, and complete the formulation by establishing the regularity, surface, and exterior matching conditions required for a well-posed oscillation eigenvalue problem. We then implement the resulting system numerically and compute representative polar mode spectra for gravitationally coupled two-fluid stellar models, allowing the fundamental,

What carries the argument

The coupled linear perturbation equations for the metric and the two independent fluid displacement fields, closed by regularity at the center, continuity across each fluid surface, and matching to the exterior vacuum spacetime.

If this is right

The fundamental and pressure mode branches can be classified according to their dominant inner- or outer-fluid character through the associated eigenfunctions and their node structure.
The implementation supplies a practical method for mode identification in gravitationally coupled two-fluid stars.
The formalism provides a foundation for extending relativistic asteroseismology to multi-fluid compact stars.
It enables systematic exploration of potential gravitational-wave signatures from such oscillations in a fully general relativistic setting.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

The same matching conditions could be reused if entrainment is later introduced, provided the relative velocity between fluids is incorporated into the surface jump relations.
Applying the framework to tabulated equations of state that include both nuclear and superfluid phases would produce concrete predictions for frequency shifts relative to single-fluid models.
The derivation steps can be repeated for axial perturbations or for slowly rotating two-fluid configurations without changing the core structure of the boundary conditions.

Load-bearing premise

The two fluids interact solely through the shared gravitational field with no direct microphysical coupling or entrainment between them.

What would settle it

A numerical calculation in which the density of one fluid is reduced continuously to zero and the resulting mode frequencies and eigenfunctions are checked for smooth recovery of the known single-fluid polar-mode spectrum.

Figures

Figures reproduced from arXiv: 2605.03305 by Ankit Kumar, Daniel A. Caballero, Hajime Sotani, Nicol\'as Yunes.

**Figure 1.** Figure 1: FIG. 1. Flowchart of the numerical procedure used to construct the non-radial mode spectrum of a gravitationally coupled two-fluid star. view at source ↗

**Figure 2.** Figure 2: FIG. 2. Gravitational mass contours of two-fluid neutron stars with view at source ↗

**Figure 3.** Figure 3: FIG. 3. Boundary-condition residual (error function) view at source ↗

**Figure 4.** Figure 4: FIG. 4. Normalized radial-displacement eigenfunctions for polar non-radial modes in the QHC21-BT model at fixed gravitational mass view at source ↗

**Figure 5.** Figure 5: shows that, in the single-fluid case, the fundamental-mode frequency increases smoothly with stellar mass for neutron stars constructed with the QHC21-BT equation of state, thereby providing a useful reference against which the two-fluid trends can be assessed. A clear feature of view at source ↗

**Figure 6.** Figure 6: depicts the p1–mode frequency as a function of gravitational mass for the same single-fluid and two-fluid mirror dark matter sequences considered in view at source ↗

**Figure 7.** Figure 7: FIG. 7. Mass-scaled fundamental-mode relation for mirror dark matter view at source ↗

**Figure 8.** Figure 8: FIG. 8. Mass-scaled fundamental-mode relation plotted against view at source ↗

read the original abstract

Non-radial oscillations of neutron stars provide a powerful probe of stellar structure and relativistic gravity, but a fully general relativistic treatment for gravitationally coupled two-fluid stars with independently conserved currents has so far been lacking. In this work, we develop a fully relativistic framework for polar perturbations of gravitationally coupled two-fluid neutron stars, assuming that the two fluids interact only through the common spacetime and are not coupled by entrainment or direct microphysical interactions. We derive the coupled linear perturbation equations governing the metric and both fluid components, and complete the formulation by establishing the regularity, surface, and exterior matching conditions required for a well-posed oscillation eigenvalue problem. We then implement the resulting system numerically and compute representative polar mode spectra for gravitationally coupled two-fluid stellar models. This implementation provides a practical way to address mode identification in gravitationally coupled two-fluid stars, allowing the fundamental ($\mathsf{f}$) and pressure ($\mathsf{p}$) mode branches of the spectrum to be classified according to their dominant inner- or outer-fluid character through the associated eigenfunctions and their node structure. The formalism developed here provides a foundation for extending relativistic asteroseismology to multi-fluid compact stars and for exploring their potential gravitational-wave signatures in a fully general relativistic setting.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Desk Editor's Note private letter to a colleague

This paper supplies the first complete GR linear perturbation setup for polar modes in two-fluid neutron stars coupled only by gravity, with boundary conditions and sample spectra.

read the letter

The core advance is closing the gap on relativistic two-fluid asteroseismology. They start from the linearized Einstein equations plus two independent fluid conservation laws, derive the coupled system for metric and fluid variables under the no-entrainment assumption, and supply the full set of regularity conditions at the center, surface matching, and exterior vacuum conditions to make it a proper eigenvalue problem. They then code it up and extract representative f- and p-mode branches, tagging each by which fluid dominates the eigenfunction and node count. That classification step is practical and new for the two-fluid case in GR. The work is technically clean on its own terms and directly usable for follow-on calculations of gravitational-wave frequencies from multi-component stars. The main limitation is the explicit decoupling: the fluids interact only through the metric, so entrainment or any microphysical coupling is left out. That is a reasonable first cut and they state it up front, but it means the results are a baseline rather than a direct model of real superfluid neutron-star matter. The numerics are shown for a few example configurations; more extensive convergence checks or recovery of the single-fluid limit would strengthen the claim, though nothing in the stated scope suggests an inconsistency. This is aimed at researchers who already work on neutron-star oscillations or gravitational-wave asteroseismology and need the relativistic two-fluid extension. It is solid enough to deserve referee time; the derivation is self-contained and the numerical implementation turns the formalism into something others can test and extend.

Referee Report

0 major / 2 minor

Summary. The manuscript develops a fully relativistic framework for polar non-radial pulsations of gravitationally coupled two-fluid neutron stars, deriving the coupled linear perturbation equations for the metric and both fluid components under the explicit assumption of no entrainment or direct microphysical interactions (only gravitational coupling via the shared spacetime). It completes the formulation with regularity, surface, and exterior matching conditions to yield a well-posed eigenvalue problem, implements the system numerically, and computes representative polar mode spectra while classifying f- and p-mode branches according to their dominant inner- or outer-fluid character via eigenfunctions and node structure.

Significance. If the derivation and implementation hold, this supplies the first complete GR treatment of polar perturbations for two-fluid stars with independently conserved currents, filling a clear gap in relativistic asteroseismology. The provision of a closed linear system, explicit boundary conditions, and numerical mode spectra offers a practical foundation for studying multi-fluid compact objects and their potential gravitational-wave signatures, with the mode-classification approach providing a concrete tool for interpretation.

minor comments (2)

The numerical implementation and resulting spectra would benefit from explicit statements of the discretization scheme, convergence tests, and error estimates (e.g., in the section describing the code or results), as these are essential for reproducibility of the eigenvalue spectra.
Notation for the two fluids (e.g., labels for inner/outer or fluid 1/2) should be introduced consistently at first use and maintained throughout the equations and figures to avoid ambiguity when discussing dominant character of modes.

Simulated Author's Rebuttal

0 responses · 0 unresolved

We thank the referee for their positive and constructive assessment of our manuscript. We are pleased that the work is recognized as supplying the first complete general-relativistic treatment of polar perturbations for gravitationally coupled two-fluid neutron stars with independently conserved currents. The recommendation for minor revision is noted; however, no specific major comments were raised in the report.

Circularity Check

0 steps flagged

No significant circularity; derivation self-contained from standard GR linearization

full rationale

The paper constructs the polar perturbation equations by applying the standard linearized Einstein equations and baryon current conservation laws to two perfect fluids that share only the metric, with the no-entrainment assumption stated explicitly as an input rather than derived. Regularity, surface, and exterior matching conditions are obtained from the same linearized system and asymptotic flatness, without reduction to fitted parameters or prior self-citations that carry the central claim. Numerical implementation and mode classification follow directly from the resulting eigenvalue problem. No step equates a prediction to its own input by construction.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The framework rests on standard general-relativistic hydrodynamics for multi-fluid systems together with the explicit modeling choice that the fluids couple solely through the metric.

axioms (2)

standard math Linearized Einstein equations and conservation laws for two independently conserved fluid currents govern the perturbations.
Invoked to obtain the coupled system of perturbation equations.
domain assumption The two fluids interact only gravitationally, with no entrainment or direct microphysical coupling.
Stated explicitly as the modeling assumption that defines the class of stars under study.

pith-pipeline@v0.9.0 · 5534 in / 1288 out tokens · 105495 ms · 2026-05-07T14:41:28.811095+00:00 · methodology

Review history (2 revisions) →

discussion (0)

Reference graph

Works this paper leans on

75 extracted references · 13 canonical work pages · 1 internal anchor

[1]

Vertical dashed lines mark the radius of the inner-fluid component. Insets zoom into the neighborhood of Wnorm ≃0 to make node structure and sign changes visible, enabling an unambiguous identification of which spectral feature corresponds to an outer- or inner-led mode. 16 this mode does not belong to the outer-fluid-led branch. This also clarifies why t...
[2]

Andersson and K

N. Andersson and K. D. Kokkotas, Phys. Rev. Lett.77, 4134 (1996)

1996
[3]

B. F. Schutz, Classical and Quantum Gravity6, 1761 (1989)

1989
[4]

T. G. Cowling, Monthly Notices of the Royal Astronomical Society101, 367 (1941)

1941
[5]

K. D. Kokkotas and B. G. Schmidt, Living Reviews in Relativity 2, 2 (1999)

1999
[6]

Lindblom and S

L. Lindblom and S. L. Detweiler, apjs53, 73 (1983)

1983
[7]

K. S. Thorne and A. Campolattaro, Non-Radial Pulsation of General-Relativistic Stellar Models. I. Analytic Analysis for L >=2 (1967)

1967
[8]

Campolattaro and K

A. Campolattaro and K. S. Thorne, Astrophys. J.159, 847 (1970)

1970
[9]

K. D. Kokkotas, mnras268, 1015 (1994)

1994
[10]

Price and K

R. Price and K. S. Thorne, Astrophys. J.155, 163 (1969)

1969
[11]

Detweiler and L

S. Detweiler and L. Lindblom, Astrophys. J.292, 12 (1985)

1985
[12]

Andersson and K

N. Andersson and K. D. Kokkotas, Monthly Notices of the Royal Astronomical Society299, 1059 (1998)

1998
[13]

L. S. Finn, Monthly Notices of the Royal Astronomical Society 222, 393 (1986)

1986
[14]

Torres-Forn´e, P

A. Torres-Forn´e, P. Cerd´a-Dur´an, A. Passamonti, and J. A. Font, Monthly Notices of the Royal Astronomical Society474, 5272 (2017)

2017
[15]

Benhar, E

O. Benhar, E. Berti, and V . Ferrari, Monthly Notices of the Royal Astronomical Society310, 797 (1999)

1999
[16]

Benhar, V

O. Benhar, V . Ferrari, and L. Gualtieri, Phys. Rev. D70, 124015 (2004)

2004
[17]

Jaikumar, A

P. Jaikumar, A. Semposki, M. Prakash, and C. Constantinou, Phys. Rev. D103, 123009 (2021), arXiv:2101.06349 [nucl-th]

work page arXiv 2021
[18]

Sotani, N

H. Sotani, N. Yasutake, T. Maruyama, and T. Tatsumi, Phys. Rev. D83, 024014 (2011)

2011
[19]

Chatterjee and D

D. Chatterjee and D. Bandyopadhyay, Phys. Rev. D80, 023011 (2009)

2009
[20]

B. K. Pradhan, D. Chatterjee, M. Lanoye, and P. Jaikumar, Phys. 23 Rev. C106, 015805 (2022)

2022
[21]

L. S. Finn, Monthly Notices of the Royal Astronomical Society 227, 265 (1987)

1987
[22]

Sotani, K

H. Sotani, K. Tominaga, and K.-i. Maeda, Phys. Rev. D65, 024010 (2001)

2001
[23]

K. D. Kokkotas and B. F. Schutz, Monthly Notices of the Royal Astronomical Society255, 119 (1992)

1992
[24]

Leins, H

M. Leins, H. P. Nollert, and M. H. Soffel, Phys. Rev. D48, 3467 (1993)

1993
[25]

A New class of unstable modes of rotatin g relativistic stars

N. Andersson, Astrophys. J.502, 708 (1998), arXiv:gr-qc/9706075 [gr-qc]

work page arXiv 1998
[26]

Kojima and M

Y . Kojima and M. Hosonuma, The Astrophysical Journal520, 788 (1999)

1999
[27]

R. I. Epstein, Astrophys. J.333, 880 (1988)

1988
[28]

G. BAYM, C. PETHICK, and D. PINES, Nature224, 673 (1969)

1969
[29]

Lombardo and H.-J

U. Lombardo and H.-J. Schulze, Superfluidity in neutron star matter, in Physics of Neutron Star Interiors, edited by D. Blaschke, A. Sedrakian, and N. K. Glendenning (Springer Berlin Heidelberg, Berlin, Heidelberg, 2001) pp. 30–53

2001
[30]

Chamel, Monthly Notices of the Royal Astronomical Society 388, 737 (2008)

N. Chamel, Monthly Notices of the Royal Astronomical Society 388, 737 (2008)

2008
[31]

Lindblom and G

L. Lindblom and G. Mendell, Astrophys. J.421, 689 (1994)

1994
[32]

G. L. Comer, D. Langlois, and L. M. Lin, Phys. Rev. D60, 104025 (1999)

1999
[33]

Andersson, G

N. Andersson, G. L. Comer, and D. Langlois, Phys. Rev. D66, 104002 (2002)

2002
[34]

Leung, M.-C

S.-C. Leung, M.-C. Chu, and L.-M. Lin, Phys. Rev. D84, 107301 (2011)

2011
[35]

L.-M. Lin, N. Andersson, and G. L. Comer, Phys. Rev. D78, 083008 (2008)

2008
[36]

Biesdorf, J

C. Biesdorf, J. Schaffner-Bielich, and L. Tolos, Phys. Rev. D 111, 083038 (2025)

2025
[37]

Goldman, R

I. Goldman, R. Mohapatra, S. Nussinov, D. Rosenbaum, and V . Teplitz, Physics Letters B725, 200 (2013)

2013
[38]

Kain,Dark matter admixed neutron stars,Phys

B. Kain, Phys. Rev. D103, 043009 (2021), arXiv:2102.08257 [gr-qc]

work page arXiv 2021
[39]

Shawqi and S

S. Shawqi and S. M. Morsink, The Astrophysical Journal975, 123 (2024)

2024
[40]

Sandin and P

F. Sandin and P. Ciarcelluti, Astroparticle Physics32, 278 (2009)

2009
[41]

Hippert, E

M. Hippert, E. Dillingham, H. Tan, D. Curtin, J. Noronha-Hostler, and N. Yunes, Phys. Rev. D107, 115028 (2023)

2023
[42]

P. N. McDermott, H. M. van Horn, and C. J. Hansen, Astrophys. J.325, 725 (1988)

1988
[43]

and Lee, U., A&A395, 201 (2002)

Yoshida, S. and Lee, U., A&A395, 201 (2002)

2002
[44]

Stavridis, A

A. Stavridis, A. Passamonti, and K. Kokkotas, Phys. Rev. D75, 064019 (2007), arXiv:gr-qc/0701122

work page arXiv 2007
[45]

Pratten, P

G. Pratten, P. Schmidt, and T. Hinderer, Nature Communications 11, 2553 (2020)

2020
[46]

W. C. G. Ho, D. I. Jones, N. Andersson, and C. M. Espinoza, Phys. Rev. D101, 103009 (2020)

2020
[47]

Sotani, Phys

H. Sotani, Phys. Rev. D103, 123015 (2021)

2021
[48]

Zhao and J

T. Zhao and J. M. Lattimer, Phys. Rev. D106, 123002 (2022)

2022
[49]

D. D. Doneva, E. Gaertig, K. D. Kokkotas, and C. Kr¨uger, Phys. Rev. D88, 044052 (2013), arXiv:1305.7197 [astro-ph.SR]

work page arXiv 2013
[50]

V . Tran, S. Ghosh, N. Lozano, D. Chatterjee, and P. Jaikumar, Phys. Rev. C108, 015803 (2023)

2023
[51]

Sotani, K

H. Sotani, K. Nakazato, K. Iida, and K. Oyamatsu, Phys. Rev. Lett.108, 201101 (2012), arXiv:1202.6242 [astro-ph.HE]

work page arXiv 2012
[52]

Andersson, K

N. Andersson, K. Glampedakis, and L. Samuelsson, Monthly Notices of the Royal Astronomical Society396, 894 (2009)

2009
[53]

Passamonti and N

A. Passamonti and N. Andersson, Monthly Notices of the Royal Astronomical Society419, 638 (2011)

2011
[54]

Samuelsson and N

L. Samuelsson and N. Andersson, Classical and Quantum Gravity26, 155016 (2009)

2009
[55]

Kumar, S

A. Kumar, S. Girmohanta, and H. Sotani, Eur. Phys. J. C85, 1109 (2025), arXiv:2501.16829 [astro-ph.HE]

work page arXiv 2025
[56]

D. A. Caballero, J. L. Ripley, and N. Yunes, Phys. Rev. D110, 103038 (2024), arXiv:2408.04701 [gr-qc]

work page arXiv 2024
[57]

Kumar and H

A. Kumar and H. Sotani, Phys. Rev. D112, 063030 (2025), arXiv:2509.03862 [astro-ph.HE]

work page arXiv 2025
[58]

Gleason, B

T. Gleason, B. Brown, and B. Kain, Phys. Rev. D105, 023010 (2022), arXiv:2201.02274 [gr-qc]

work page arXiv 2022
[59]

Rutherford, C

N. Rutherford, C. Prescod-Weinstein, and A. Watts, Phys. Rev. D111, 123034 (2025)

2025
[60]

Ellis, G

J. Ellis, G. H ¨utsi, K. Kannike, L. Marzola, M. Raidal, and V . Vaskonen, Phys. Rev. D97, 123007 (2018)

2018
[61]

Sotani and A

H. Sotani and A. Kumar, Eur. Phys. J. C85, 1438 (2025), arXiv:2512.07105 [astro-ph.HE]

work page arXiv 2025
[62]

Thakur, A

P. Thakur, A. Kumar, V . B. Thapa, V . Parmar, and M. Sinha, Journal of Cosmology and Astroparticle Physics2024(12), 042
[63]

Sotani and A

H. Sotani and A. Kumar, Phys. Rev. D111, 123013 (2025), arXiv:2505.18800 [astro-ph.HE]

work page arXiv 2025
[64]

P. N. McDermott, H. M. van Horn, and J. F. Scholl, Astrophys. J.268, 837 (1983)

1983
[65]

Kumar and H

A. Kumar and H. Sotani, Phys. Rev. D113, 063055 (2026)

2026
[66]

J. L. Friedman, Communications in Mathematical Physics62, 247 (1978)

1978
[67]

Sotani and T

H. Sotani and T. Takiwaki, Phys. Rev. D102, 063025 (2020)

2020
[68]

F. J. Zerilli, Phys. Rev. Lett.24, 737 (1970)

1970
[69]

Chandrasekhar and S

S. Chandrasekhar and S. Detweiler, Proceedings of the Royal Society of London Series A344, 441 (1975)

1975
[70]

Regge and J

T. Regge and J. A. Wheeler, Phys. Rev.108, 1063 (1957)

1957
[71]

Beyond Three Terms: Continued Fractions for Rotating Black Holes in Modified Gravity

G. Karikos, J. A. Saes, P. Wagle, and N. Yunes, Beyond three terms: Continued fractions for rotating black holes in modified gravity (2026), arXiv:2604.18680 [gr-qc]

work page internal anchor Pith review Pith/arXiv arXiv 2026
[72]

E. W. Leaver, Phys. Rev. D41, 2986 (1990)

1990
[73]

T. Kojo, G. Baym, and T. Hatsuda, The Astrophysical Journal 934, 46 (2022)

2022
[74]

Togashi, K

H. Togashi, K. Nakazato, Y . Takehara, S. Yamamuro, H. Suzuki, and M. Takano, Nuclear Physics A961, 78 (2017)

2017
[75]

Sotani, The European Physical Journal C82, 477 (2022)

H. Sotani, The European Physical Journal C82, 477 (2022)

2022