pith. sign in

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

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

Ankit Kumar , Daniel A. Caballero , Hajime Sotani , Nicol\'as Yunes This is my paper

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
0 comments X

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. 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. Figure 2: FIG. 2. Gravitational mass contours of two-fluid neutron stars with view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. Boundary-condition residual (error function) view at source ↗
Figure 4
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. 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. 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. Figure 7: FIG. 7. Mass-scaled fundamental-mode relation for mirror dark matter view at source ↗
Figure 8
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.

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)
  1. 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.
  2. 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)

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

Reference graph

Works this paper leans on

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

  1. [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...

    work page

  2. [2]

    Andersson and K

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

    work page 1996

  3. [3]

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

    work page 1989

  4. [4]

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

    work page 1941

  5. [5]

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

    work page 1999

  6. [6]

    Lindblom and S

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

    work page 1983

  7. [7]

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

    work page 1967

  8. [8]

    Campolattaro and K

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

    work page 1970

  9. [9]

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

    work page 1994

  10. [10]

    Price and K

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

    work page 1969

  11. [11]

    Detweiler and L

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

    work page 1985

  12. [12]

    Andersson and K

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

    work page 1998

  13. [13]

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

    work page 1986

  14. [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)

    work page 2017

  15. [15]

    Benhar, E

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

    work page 1999

  16. [16]

    Benhar, V

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

    work page 2004

  17. [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. [18]

    Sotani, N

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

    work page 2011

  19. [19]

    Chatterjee and D

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

    work page 2009

  20. [20]

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

    work page 2022

  21. [21]

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

    work page 1987

  22. [22]

    Sotani, K

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

    work page 2001

  23. [23]

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

    work page 1992

  24. [24]

    Leins, H

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

    work page 1993

  25. [25]

    A new class of unstable modes of rotating relativistic stars

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

    work page Pith review arXiv 1998

  26. [26]

    Kojima and M

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

    work page 1999

  27. [27]

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

    work page 1988

  28. [28]

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

    work page 1969

  29. [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

    work page 2001

  30. [30]

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

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

    work page 2008

  31. [31]

    Lindblom and G

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

    work page 1994

  32. [32]

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

    work page 1999

  33. [33]

    Andersson, G

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

    work page 2002

  34. [34]

    Leung, M.-C

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

    work page 2011

  35. [35]

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

    work page 2008

  36. [36]

    Biesdorf, J

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

    work page 2025

  37. [37]

    Goldman, R

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

    work page 2013

  38. [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. [39]

    Shawqi and S

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

    work page 2024

  40. [40]

    Sandin and P

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

    work page 2009

  41. [41]

    Hippert, E

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

    work page 2023

  42. [42]

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

    work page 1988

  43. [43]

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

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

    work page 2002

  44. [44]

    Stavridis, A

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

    work page arXiv 2007

  45. [45]

    Pratten, P

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

    work page 2020

  46. [46]

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

    work page 2020

  47. [47]

    Sotani, Phys

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

    work page 2021

  48. [48]

    Zhao and J

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

    work page 2022

  49. [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. [50]

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

    work page 2023

  51. [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. [52]

    Andersson, K

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

    work page 2009

  53. [53]

    Passamonti and N

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

    work page 2011

  54. [54]

    Samuelsson and N

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

    work page 2009

  55. [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. [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. [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. [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. [59]

    Rutherford, C

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

    work page 2025

  60. [60]

    Ellis, G

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

    work page 2018

  61. [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. [62]

    Thakur, A

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

    work page

  63. [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. [64]

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

    work page 1983

  65. [65]

    Kumar and H

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

    work page 2026

  66. [66]

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

    work page 1978

  67. [67]

    Sotani and T

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

    work page 2020

  68. [68]

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

    work page 1970

  69. [69]

    Chandrasekhar and S

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

    work page 1975

  70. [70]

    Regge and J

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

    work page 1957

  71. [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. [72]

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

    work page 1990

  73. [73]

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

    work page 2022

  74. [74]

    Togashi, K

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

    work page 2017

  75. [75]

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

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

    work page 2022