Recognition: unknown
Axial w-modes of anisotropic neutron stars
Pith reviewed 2026-05-09 17:18 UTC · model grok-4.3
The pith
Axial w-mode frequencies of anisotropic neutron stars depend approximately linearly on compactness, with anisotropy modifying the slope and intercept.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The axial w-modes of anisotropic neutron stars exhibit an approximately linear dependence of their frequencies on stellar compactness. For fixed anisotropy strength the frequency decreases monotonically with mass along the stable branch. Configurations with dominant radial pressure have higher frequencies than those with dominant tangential pressure at low mass, but the ordering reverses near the upper end of the stable branch. Damping times increase with mass and become shorter when tangential pressure dominates at fixed mass. The Bowers-Liang ansatz produces a wider spread in both frequencies and damping times than the Horvat ansatz. Empirical expressions are given for frequency and for ln
What carries the argument
The continued-fraction method for solving the linearized axial perturbation equations on static anisotropic neutron-star backgrounds constructed with the Horvat and Bowers-Liang pressure-anisotropy ansatzes.
If this is right
- For any fixed compactness the axial w-mode frequency can be estimated from the supplied empirical fit once the anisotropy parameter is specified.
- Damping times shorten when tangential pressure exceeds radial pressure at the same mass, with the effect strongest near the maximum stable mass.
- The Bowers-Liang ansatz generates a broader range of frequencies and longer damping times than the Horvat ansatz at equal compactness.
- At low masses, stars with radial pressure larger than tangential pressure oscillate at higher frequencies, while the reverse holds near the maximum mass.
- The sensitivity of damping time to anisotropy grows markedly for more compact configurations.
Where Pith is reading between the lines
- If axial w-modes are detected in gravitational waves from neutron-star events, the reported relations could be inverted to constrain the level of pressure anisotropy inside the star.
- The dominance of compactness in the linear fit suggests that anisotropy enters mainly as a correction term that could be calibrated against future observations.
- The same modeling approach could be applied to polar w-modes or to rotating anisotropic stars to check whether the linear compactness trend persists.
- Direct comparison of these empirical formulas against full numerical evolutions of oscillating anisotropic stars would test the accuracy of the continued-fraction results.
Load-bearing premise
The chosen functional forms for pressure anisotropy together with the BSk21 and SLy4 equations of state produce stable stellar configurations whose axial w-mode spectra are reliably obtained by the continued-fraction technique.
What would settle it
A numerical computation or gravitational-wave observation of an axial w-mode frequency for a known compactness and anisotropy strength that lies significantly off the reported linear frequency-compactness line would falsify the central relation.
Figures
read the original abstract
We investigate the axial $w$-mode oscillations of anisotropic neutron stars. Stellar configurations are constructed using two realistic equations of state, BSk21 and SLy4, with two prescriptions for pressure anisotropy, the Horvat ansatz and the Bowers-Liang ansatz. The axial $w$-mode frequencies are computed by solving the linearized perturbation equations using a continued-fraction method. For each fixed anisotropy strength, the axial $w$-mode frequency decreases monotonically with increasing stellar mass along the stable branch, with its magnitude depending on both the equation of state and the nature of the anisotropy. At low stellar masses, configurations with dominant radial pressure ($p_r>p_t$) exhibit higher frequencies than those with dominant tangential pressure, whereas toward the upper end of the stable branch this ordering is reversed, and configurations with $p_t>p_r$ attain higher frequencies at the same mass. The axial $w$-mode frequency displays an approximately linear dependence on compactness, with anisotropy modifying both the slope and the intercept. The Bowers-Liang ansatz produces a wider spread in the frequency values compared to the Horvat ansatz. We also analyze the damping times associated with the axial $w$-modes and find that they increase with stellar mass, with a rapid rise toward the upper end of the stable branch. At a fixed mass, increasing the tangential pressure relative to the radial pressure leads to shorter damping times, while configurations with dominant radial pressure exhibit longer damping times. The sensitivity of the damping time to anisotropy is more pronounced for more compact stars, and the Bowers-Liang ansatz yields systematically larger damping times than the Horvat ansatz. Finally, we provide empirical expressions for the axial $w$-mode frequency and damping time as functions of stellar compactness and anisotropy strength.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript investigates axial w-modes of anisotropic neutron stars constructed with BSk21 and SLy4 equations of state and Horvat/Bowers-Liang anisotropy prescriptions. It solves the linearized axial perturbation equations via the continued-fraction method, reporting that w-mode frequency decreases monotonically with mass on the stable branch, exhibits approximate linear dependence on compactness (with anisotropy modifying slope and intercept), and provides empirical fitting expressions for frequency and damping time as functions of compactness and anisotropy strength. Damping times increase with mass and show greater sensitivity to anisotropy at higher compactness, with Bowers-Liang yielding larger spreads than Horvat.
Significance. If the numerical trends hold, the work extends asteroseismology to anisotropic neutron-star models, potentially aiding interpretation of gravitational-wave signals from oscillating or merging compact objects. The use of realistic EOS, standard methods, and provision of empirical relations are positive features, though the results remain specific to the chosen ansatzes and may require further validation for broader applicability.
major comments (2)
- The empirical fitting expressions for frequency and damping time (described in the abstract and results) are constructed from the computed data points; the manuscript should detail the functional form chosen, the fitting procedure, the exact ranges of compactness and anisotropy strength, and quantitative fit metrics (e.g., residuals or R²) to substantiate their utility beyond interpolation.
- Application of the continued-fraction method is load-bearing for all reported frequencies, damping times, and linearity claims, yet the text provides limited information on convergence criteria, truncation errors, or direct validation against the isotropic (zero-anisotropy) limit; explicit error budgets or comparison tables would be required to confirm the monotonic trends and ordering reversals with anisotropy.
minor comments (3)
- The abstract could briefly indicate the range of anisotropy parameters explored to contextualize the reported ordering reversals between radial- and tangential-dominant cases.
- Figures displaying frequency versus mass or compactness should label the stable branch explicitly and, if possible, indicate numerical resolution or precision.
- Ensure the introduction cites key prior isotropic w-mode studies to clearly frame the extension to anisotropic configurations.
Simulated Author's Rebuttal
We thank the referee for the positive evaluation and constructive feedback on our manuscript. The suggestions will help strengthen the presentation of the numerical results and empirical relations. We address each major comment below.
read point-by-point responses
-
Referee: The empirical fitting expressions for frequency and damping time (described in the abstract and results) are constructed from the computed data points; the manuscript should detail the functional form chosen, the fitting procedure, the exact ranges of compactness and anisotropy strength, and quantitative fit metrics (e.g., residuals or R²) to substantiate their utility beyond interpolation.
Authors: We agree that additional details on the empirical fits are warranted to allow readers to assess their robustness. In the revised manuscript we will add a dedicated paragraph (or short subsection) in the results section that specifies: (i) the exact functional forms adopted (linear in compactness with additive and multiplicative terms linear in the anisotropy parameter), (ii) the least-squares fitting procedure, (iii) the precise ranges of compactness (approximately 0.10–0.32) and anisotropy strengths explored for each EOS and ansatz, and (iv) quantitative goodness-of-fit measures including R² values and maximum relative residuals. These additions will be placed immediately after the presentation of the fitting expressions. revision: yes
-
Referee: Application of the continued-fraction method is load-bearing for all reported frequencies, damping times, and linearity claims, yet the text provides limited information on convergence criteria, truncation errors, or direct validation against the isotropic (zero-anisotropy) limit; explicit error budgets or comparison tables would be required to confirm the monotonic trends and ordering reversals with anisotropy.
Authors: We acknowledge that the current description of the continued-fraction implementation is concise and would benefit from more explicit numerical validation. In the revised version we will expand the methods section to include: (i) the convergence criteria employed (number of continued-fraction terms and tolerance on the frequency), (ii) estimates of truncation error obtained by increasing the truncation order, and (iii) a direct comparison table (or figure) of our zero-anisotropy results against published isotropic w-mode frequencies for BSk21 and SLy4. This will provide an explicit error budget supporting the reported monotonic trends and anisotropy-induced ordering reversals. revision: yes
Circularity Check
No significant circularity
full rationale
The paper constructs stellar models from two standard EOS (BSk21, SLy4) and two established anisotropy ansatzes (Horvat, Bowers-Liang), then solves the linearized axial perturbation equations via the continued-fraction method to obtain w-mode frequencies and damping times. All reported trends (monotonic decrease with mass, approximate linearity with compactness, ordering reversals, sensitivity to anisotropy) and differences between ansatzes are direct numerical outputs of this procedure. The empirical expressions are explicitly described as fits to the computed data rather than first-principles derivations or predictions, and the central claims do not depend on these fits for their validity. No step reduces by construction to its inputs, no self-citation is load-bearing, and the numerical chain is self-contained against external benchmarks.
Axiom & Free-Parameter Ledger
free parameters (1)
- anisotropy strength parameter
axioms (2)
- domain assumption The background metric and hydrostatic equilibrium are described by the anisotropic Tolman-Oppenheimer-Volkoff equations.
- domain assumption The continued-fraction method accurately solves the linearized axial perturbation equations with appropriate boundary conditions at the center and at infinity.
Forward citations
Cited by 1 Pith paper
-
On the non-radial oscillations of realistic anisotropic neutron stars: Axial modes
Axial modes of anisotropic neutron stars show mass-scaled frequency and damping time with nearly universal quadratic dependence on compactness, insensitive to EOS and anisotropy model.
Reference graph
Works this paper leans on
-
[1]
In the present work, we derive the governing equations for the axialw-modes of anisotropic neutron stars and compute their complex frequencies
computed the polarw-modes of anisotropic neutron stars. In the present work, we derive the governing equations for the axialw-modes of anisotropic neutron stars and compute their complex frequencies. To solve these equa- tions and obtain the corresponding oscillation frequen- cies and damping times, we consider two realistic nuclear equations of state, na...
1999
- [2]
-
[3]
L. D. Landau, On the theory of stars, Phys. Z. Sowjetu- nion1, 285 (1932)
1932
-
[4]
Baade and F
W. Baade and F. Zwicky, Cosmic Rays from Super- Novae, Proc. Nat. Acad. Sci.20, 259 (1934). 15
1934
-
[5]
Baade and F
W. Baade and F. Zwicky, On Super-Novae, Proc. Nat. Acad. Sci.20, 254 (1934)
1934
-
[6]
R. C. Tolman, Static solutions of Einstein’s field equa- tions for spheres of fluid, Phys. Rev.55, 364 (1939)
1939
-
[7]
J. R. Oppenheimer and G. M. Volkoff, On massive neu- tron cores, Phys. Rev.55, 374 (1939)
1939
-
[8]
Hewish, S
A. Hewish, S. J. Bell, J. D. H. Pilkington, P. F. Scott, and R. A. Collins, Observation of a rapidly pulsating radio source, Nature217, 709 (1968)
1968
-
[9]
R. N. Manchesteret al., The Australia Telescope Na- tional Facility Pulsar Catalogue, Astron. J.129, 1993 (2005), arXiv:astro-ph/0412641 [astro-ph]
work page Pith review arXiv 1993
-
[10]
B. P. Abbottet al.(LIGO Scientific, Virgo), Observation of Gravitational Waves from a Binary Black Hole Merger, Phys. Rev. Lett.116, 061102 (2016), arXiv:1602.03837 [gr-qc]
work page internal anchor Pith review arXiv 2016
-
[11]
B. P. Abbottet al., Multi-messenger Observations of a Binary Neutron Star Merger, Astrophys. J. Lett.848, L12 (2017), arXiv:1710.05833
work page Pith review arXiv 2017
- [12]
-
[13]
B. Haskell and D. I. Jones, Glitching pulsars as gravita- tional wave sources, Astropart. Phys.157, 102921 (2024), arXiv:2311.04586 [astro-ph.HE]
-
[14]
B. Haskell and K. Schwenzer, Isolated Neutron Stars (2022) arXiv:2104.03137 [gr-qc]
- [15]
- [16]
-
[17]
A. Bauswein, N. Stergioulas, and H.-T. Janka, Explor- ing properties of high-density matter through remnants of neutron-star mergers, Eur. Phys. J. A52, 56 (2016), arXiv:1508.05493 [astro-ph.HE]
-
[18]
K. S. Thorne and A. Campolattaro, Non-Radial Pulsa- tion of General-Relativistic Stellar Models. I. Analytic Analysis forL >= 2, Astrophys. J.149, 591 (1967)
1967
-
[19]
L. Gualtieri, E. M. Kantor, M. E. Gusakov, and A. I. Chugunov, Quasinormal modes of superfluid neutron stars, Phys. Rev. D90, 024010 (2014), arXiv:1404.7512 [gr-qc]
-
[20]
N. Andersson, G. L. Comer, and D. Langlois, Oscillations of general relativistic superfluid neutron stars, Phys. Rev. D66, 104002 (2002), arXiv:gr-qc/0203039
- [21]
-
[22]
G. Miniutti, J. A. Pons, E. Berti, L. Gualtieri, and V. Fer- rari, Non-radial oscillation modes as a probe of density discontinuities in neutron stars, Mon. Not. Roy. Astron. Soc.338, 389 (2003), arXiv:astro-ph/0206142
-
[23]
A. Kunjipurayil, T. Zhao, B. Kumar, B. K. Agrawal, and M. Prakash, Impact of the equation of state on f- and p- mode oscillations of neutron stars, Phys. Rev. D106, 063005 (2022), arXiv:2205.02081 [nucl-th]
- [24]
-
[25]
H. Sotani and A. Kumar, Universal relation involv- ing fundamental modes in two-fluid dark matter ad- mixed neutron stars, Eur. Phys. J. C85, 1438 (2025), arXiv:2512.07105 [astro-ph.HE]
- [26]
-
[27]
K. D. Kokkotas and B. F. Schutz, Normal modes of a model radiating system, Gen. Rel. Grav.18, 913 (1986)
1986
-
[28]
K. D. Kokkotas and B. F. Schutz, W-modes - A new fam- ily of normal modes of pulsating relativistic stars, Mon. Not. Roy. Astron. Soc.255, 119 (1992)
1992
-
[29]
Chandrasekhar and V
S. Chandrasekhar and V. Ferrari, On the non-radial os- cillations of a star. III - A reconsideration of the axial modes, Proceedings of the Royal Society of London Se- ries A434, 449 (1991)
1991
-
[30]
K. D. Kokkotas, Axial Modes for Relativistic Stars, Mon. Not. Roy. Astron. Soc.268, 1015 (1994)
1994
-
[31]
K. D. Kokkotas and B. G. Schmidt, Quasinormal modes of stars and black holes, Living Rev. Rel.2, 2 (1999), arXiv:gr-qc/9909058
work page internal anchor Pith review arXiv 1999
- [32]
-
[33]
D. Chatterjee and D. Bandyopadhyay, Probing dense matter in neutron stars with axialwmodes, Phys. Rev. D80, 023011 (2009), arXiv:0904.1949
-
[34]
H. Sotani and K. D. Kokkotas, Stellar oscillations in scalar-tensor theory of gravity, Phys. Rev. D71, 124038 (2005), arXiv:gr-qc/0506060
-
[35]
Z. Altaha Motahar, J. L. Bl´ azquez-Salcedo, B. Kleihaus, and J. Kunz, Axial quasinormal modes of scalarized neu- tron stars with realistic equations of state, Phys. Rev. D 98, 044032 (2018), arXiv:1807.02598
- [36]
- [37]
-
[38]
Lemaitre, The expanding universe, Annales Soc
G. Lemaitre, The expanding universe, Annales Soc. Sci. Bruxelles A53, 51 (1933)
1933
-
[39]
R. L. Bowers and E. P. T. Liang, Anisotropic Spheres in General Relativity, Astrophys. J.188, 657 (1974)
1974
-
[40]
Ruderman, Pulsars: structure and dynamics, Ann
M. Ruderman, Pulsars: structure and dynamics, Ann. Rev. Astron. Astrophys.10, 427 (1972)
1972
-
[41]
Hoffberg, A
M. Hoffberg, A. E. Glassgold, R. W. Richardson, and M. Ruderman, Anisotropic Superfluidity in Neutron Star Matter, Phys. Rev. Lett.24, 775 (1970)
1970
-
[42]
R. F. Sawyer, Condensed pi- phase in neutron star mat- ter, Phys. Rev. Lett.29, 382 (1972)
1972
-
[43]
S. Nelmes and B. M. A. G. Piette, Phase transition and anisotropic deformations of neutron star matter, Phys. Rev. D85, 123004 (2012), arXiv:1204.0910 [astro-ph.SR]
-
[44]
Barreto and S
W. Barreto and S. Rojas, An equation of state for radiat- ing dissipative spheres in general relativity, Astrophysics and space science193, 201 (1992)
1992
-
[45]
Barreto, Exploding radiating viscous spheres in gen- eral relativity, Astrophysics and space science201, 191 (1993)
W. Barreto, Exploding radiating viscous spheres in gen- eral relativity, Astrophysics and space science201, 191 (1993). 16
1993
-
[46]
Yazadjiev, Relativistic models of magnetars: Nonper- turbative analytical approach, Phys
S. Yazadjiev, Relativistic models of magnetars: Nonper- turbative analytical approach, Phys. Rev. D85, 044030 (2012), arXiv:1111.3536 [gr-qc]
-
[47]
M. Karlovini and L. Samuelsson, Elastic stars in general relativity. 1. Foundations and equilibrium models, Class. Quant. Grav.20, 3613 (2003), arXiv:gr-qc/0211026
-
[48]
A. Alho, J. Nat´ ario, P. Pani, and G. Raposo, Compact elastic objects in general relativity, Physical Review D 105, 044025 (2022)
2022
-
[49]
Herrera and N
L. Herrera and N. O. Santos, Local anisotropy in self- gravitating systems, Phys. Rept.286, 53 (1997)
1997
-
[50]
S. Mondal and M. Bagchi, f-mode oscillations of anisotropic neutron stars in full general relativity, Phys. Rev. D110, 123011 (2024), arXiv:2309.00439 [gr-qc]
-
[51]
S. Mondal and M. Bagchi, Quasinormal f-modes of anisotropic quark stars in full general relativity, Phys. Rev. D111, 103035 (2025), arXiv:2504.20589 [gr-qc]
- [52]
- [53]
- [54]
-
[55]
F. Douchin and P. Haensel, A unified equation of state of dense matter and neutron star structure, Astron. As- trophys.380, 151 (2001), arXiv:astro-ph/0111092
-
[56]
Chabanat, P
E. Chabanat, P. Bonche, P. Haensel, J. Meyer, and R. Schaeffer, A Skyrme parametrization from subnuclear to neutron star densities Part II. Nuclei far from stabili- ties, Nucl. Phys. A635, 231 (1998)
1998
- [57]
- [58]
-
[59]
Folomeev, Anisotropic neutron stars inR 2 gravity, Phys
V. Folomeev, Anisotropic neutron stars inR 2 gravity, Phys. Rev. D97, 124009 (2018), arXiv:1802.01801 [gr- qc]
- [60]
- [61]
- [62]
- [63]
-
[64]
R. Abbottet al.(LIGO Scientific, Virgo), GW190814: Gravitational Waves from the Coalescence of a 23 Solar Mass Black Hole with a 2.6 Solar Mass Compact Object, Astrophys. J. Lett.896, L44 (2020), arXiv:2006.12611 [astro-ph.HE]
-
[65]
N. Andersson and G. L. Comer, Relativistic fluid dynam- ics: Physics for many different scales, Living Rev. Rel. 10, 1 (2007), arXiv:gr-qc/0605010
-
[66]
Kojima, Equations governing the nonradial oscilla- tions of a slowly rotating relativistic star, Phys
Y. Kojima, Equations governing the nonradial oscilla- tions of a slowly rotating relativistic star, Phys. Rev. D 46, 4289 (1992)
1992
-
[67]
E. W. Leaver, An Analytic Representation for the Quasi- Normal Modes of Kerr Black Holes, Proc. Roy. Soc. Lond. A402, 285 (1985)
1985
-
[68]
Leins, H
M. Leins, H. P. Nollert, and M. H. Soffel, Nonradial os- cillations of neutron stars: A New branch of strongly damped normal modes, Phys. Rev. D48, 3467 (1993)
1993
-
[69]
Akmal, V
A. Akmal, V. R. Pandharipande, and D. G. Ravenhall, Equation of state of nucleon matter and neutron star structure, Phys. Rev. C58, 1804 (1998)
1998
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.