Pith. sign in

REVIEW 4 major objections 4 minor 1 cited by

This paper establishes that for a star made of two weakly interacting fluids—neutron matter and dark matter—the appearance of a zero-frequency radial oscillation mode is exactly equivalent to the static critical-curve condition, so the stab

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 →

For two-fluid neutron–dark-matter stars, the paper claims a rigorous equivalence between the static critical-curve stability criterion and the vanishing of a radial oscillation mode, and maps the resulting stable configurations, including mass–radius twin stars.

T0 review reviewed 2026-08-03 challenge →

load-bearing objection The numerical catalog is worth having, but the claimed rigorous proof of the critical-curve/zero-mode equivalence has a real gap that the paper does not close. the 4 major comments →

arxiv 2512.01641 v2 pith:5IEDJ5S3 submitted 2025-12-01 astro-ph.HE gr-qchep-phhep-thmath-phmath.MP

Stability Boundary of Neutron-Dark Matter Mixed Stars

classification astro-ph.HE gr-qchep-phhep-thmath-phmath.MP MSC 85A1583C55 PACS 97.60.Jd95.35.+d04.40.Dg
keywords two-fluid starsmixed starsdark matter admixed neutron starsradial oscillation stabilitycritical curve criterionzero-frequency modetwin starsmulti-fluid TOV equations
verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

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

The reading

The paper proves that two ways of finding the stability boundary of a star made of two weakly interacting fluids—neutron matter and dark matter—give the same answer. One way solves the relativistic radial oscillation equations and asks when the lowest frequency reaches zero. The other way inspects equilibrium sequences and locates where gradients of total mass and particle numbers become parallel in the central-pressure plane. The paper shows these two criteria coincide, so the static method, which is far cheaper, can be trusted to delimit the stable region. For the worked dark-matter-admixed neutron star models, the stable set is a two-dimensional surface in the mass-radius-central-pressure space, and it contains twin stars with identical mass and radius but different interiors.

Core claim

On its own terms, the paper's central claim is the equivalence stated as Eq. (42): for a two-fluid mixed star, any one of the three pairwise parallelisms among the gradients of total mass, neutron-matter particle number, and dark-matter particle number in the central-pressure plane is equivalent to the existence of a radial normal mode with zero frequency. The authors interpret this as making rigorous the longtime implicit assumption that the static critical-curve criterion locates genuine dynamical stability boundaries. Their proof proceeds in two steps: a zero eigenmode yields a static perturbation connecting neighboring equilibria while conserving mass and particle numbers, giving extrema

What carries the argument

The load-bearing object is the triplet of gradient fields ∇M, ∇N_NM, and ∇N_DM on the two-dimensional central-pressure plane of an equilibrium family. The paper's equivalence theorem asserts that parallelism of any two of these gradients is equivalent to the existence of a zero-frequency radial mode. The proof machinery is the variational extremum principle: equilibrium configurations extremize total gravitational mass M subject to fixed particle numbers N_NM and N_DM, derived through a Lagrange-multiplier scheme in the multi-fluid TOV setting; together with covariant conservation of mass and particle number, this turns static zero modes into simultaneous extrema and back into gradient paral

Load-bearing premise

The argument assumes that whenever two of the three gradients—of total mass, neutron-matter particle number, and dark-matter particle number—are parallel in the pressure plane, the third gradient must also be parallel, an assertion that does not hold for arbitrary three vectors in a two-dimensional plane.

What would settle it

Inspect a two-fluid configuration where, say, ∇M and ∇N_NM are parallel but ∇N_DM is not, by computing the three normalized cross products on a fine grid in the central-pressure plane. If the set where any two gradients are parallel is strictly larger than the set where all three are parallel, and if the fundamental radial mode vanishes only on the latter, then Eq. (42)'s OR-condition is not equivalent to the zero-mode condition.

Watch this falsifier. Get emailed when new claim-graph text bears on it.

If this is right

  • The critical-curve criterion, which costs only equilibrium TOV solutions, can be used as a reliable substitute for full radial pulsation calculations when locating the fundamental-mode stability boundary of a two-fluid mixed star.
  • Every configuration on the plotted critical curves in Figs. 5–6 is a genuine boundary where some radial mode has zero frequency; the stable region is the interior mapped in Figs. 7–8.
  • Mixed stars do not obey the classical BTM turning-point rule: the last stable configuration along a fixed-dark-matter-pressure sequence is not the maximum-mass point, and the deviation depends on the dark matter equation of state.
  • Stable mixed stars occupy a two-dimensional surface in the mass-radius-central-pressure space, and its projection onto the mass-radius plane contains overlapping twin regions with identical mass and radius but different internal compositions.
  • The framework extends to general multi-fluid systems, since the conservation and variational arguments are not special to the two-fluid case.

Where Pith is reading between the lines

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

  • If the equivalence is accepted, one can invert the logic: an observed stability boundary in the mass-radius data of a candidate mixed star would probe the dark matter equation of state, because the paper's numerics show the boundary's shape is sensitive to dark matter parameters.
  • The twin-star overlap implies a degeneracy that gravitational-wave and X-ray observations alone may not break; combining two different observables, such as tidal deformability and radius, could distinguish a dark-core twin from a dark-halo twin.
  • A direct extension would be to test the equivalence for slowly rotating or finite-temperature two-fluid stars, where the equilibrium sequence is no longer one-dimensional and the variational principle needs modification.
  • The paper's own proof route suggests a sharper check: verify that the three pairwise parallelism sets coincide on an open neighborhood, not just along the ∇M–∇N_DM curve shown; if they do not, the OR-condition in Eq. (42) is not what the numerical boundary computes.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit.

Referee Report

4 major / 4 minor

Summary. The paper derives and numerically applies a stability criterion for two-fluid neutron-star/dark-matter mixed stars. The central claim is Eq. (42): the OR of three pairwise parallelisms among the gradients of total mass and the two particle numbers in the central-pressure plane is equivalent to the existence of a zero-frequency radial pulsation mode. The proof is split into Theorem III.1 (zero mode implies mutual parallelism of the three gradients) and Theorem III.2 (equilibrium is an extremum of M under N-conserving variations), with a Lagrange-multiplier derivation in Appendix A. The authors then use the ∇M×∇N_DM cross product to plot critical curves (Figs. 5–6) and map stable regions (Figs. 7–8), highlighting twin-star configurations. The paper also discusses deviations from the single-fluid BTM criterion and potential observational constraints.

Significance. If the claimed equivalence were correct, it would provide a computationally cheap, rigorous stability diagnostic for multi-fluid stars and would justify the plotted stability boundaries and the twin-star predictions. The paper has several genuine strengths: Appendix A correctly reproduces the two-fluid TOV equation from a variational principle (Eq. A28 → Eq. 7), the conservation argument supporting the forward direction of Theorem III.1 is physically sound, and the numerical exploration of stable surfaces and twin configurations is a useful contribution to the dark-matter admixed neutron-star literature. However, the central theorem as stated is not established by the proof given, and the numerical validation tests only one of the three gradient parallelisms. The significance of the paper therefore rests on an assertion that currently lacks both a valid proof and numerical support in its full generality.

major comments (4)
  1. [§III.C, Eq. (42)] The claimed OR-equivalence is not proven. The proof establishes at most the AND-condition (43) in the forward direction, and even that only partially. Theorem III.2 proves that equilibrium configurations are extrema of M under N-conserving variations, which is true at every equilibrium point and does not imply that if any two of ∇M, ∇N_NM, ∇N_DM are parallel, the third must also be parallel. In a two-dimensional pressure plane, three vectors can satisfy the OR condition without satisfying the AND condition (e.g., ∇M=(1,0), ∇N_NM=(1,0), ∇N_DM=(1,1)). The permutation-symmetry argument following Theorem III.2 is therefore invalid, and Eq. (42) is not a logical consequence of the given derivations.
  2. [Theorem III.1 (converse direction)] The statement of Theorem III.1 is 'if and only if', but the proof only argues that a zero-frequency mode implies a static displacement conserving M and N_I, hence parallel gradients. The converse — parallel gradients imply the existence of a zero eigenvalue of the second variation — is never established. The remark 'or almost zero, i.e., ω_n^2 → 0^+' does not supply the missing spectral or turning-point argument. Thus even the AND-equivalence is incomplete.
  3. [§IV.A.2, Fig. 5 caption] The numerical validation computes only the normalized cross product of ∇M and ∇N_DM (the background heatmap in Fig. 5 and the curves in Fig. 6). It does not check ∇M×∇N_NM or ∇N_NM×∇N_DM along the plotted boundaries. Consequently, the numerical comparison with the fundamental-mode boundary in Fig. 4b/5a tests only one of the three pairwise parallelisms and cannot distinguish the OR condition (42) from the AND condition (43), nor can it validate the claim that the plotted curves are the fundamental-mode stability boundary.
  4. [Abstract and Conclusion] The abstract states that the paper 'rigorously prove[s] the equivalence' between the zero-frequency radial mode and the static critical-curve criterion, and the conclusion repeats this as a 'formal equivalence'. Given that the proof does not establish Eq. (42) and that the numerical tests cover only one cross-product, these statements overstate the demonstrated result. At minimum, the claims should be weakened to the provable subset and the missing steps explicitly flagged.
minor comments (4)
  1. [§III.A.2, Eq. (38)–(39)] The reduction of the fluid-boundary condition to Δp_I = 0 is stated after asserting that the coupling term contains a factor p_I0 that vanishes at R_I. This is plausible but should be shown explicitly; also, the sign of the dp_I0/dr term in Eq. (38) is not obviously consistent with the preceding expression for δp_I.
  2. [§III.A.2, Eq. (37)] The symbol R in the coupling terms appears to denote both the star radius and the combination 4πe^{2λ+2φ}/r; this dual use is confusing and should be renamed.
  3. [General] Typos and notational inconsistencies: 'purturbation', 'steller', the duplicated 'N_F' in the critical-curve criterion, and the undefined symbol F in Eq. (45) after the conservation-law discussion. References [82] and [83] are the same work and should be merged.
  4. [Abstract vs. §III.C] The abstract calls the relation a 'sufficient-condition relation' while the rest of the paper claims an equivalence. This inconsistency should be resolved in a revision.

Circularity Check

2 steps flagged

The central equivalence in Eq. (42) is partly circular: the converse of Theorem III.1 is assumed, and the OR-to-AND reduction is asserted via 'permutation symmetry.'

specific steps
  1. self definitional [§III.C, Theorem III.1, Eq. (44)]
    "If for a certain static equilibrium configuration, there exists n∈N such that ω_n^2=0 (or almost zero, i.e., ω_n^2→0^+) under radial perturbation, the corresponding eigenmode becomes quasi-static, implying the existence of a perturbation normal mode that is almost time-independent, called a static solution: ξI_n(r,t)=ξI_n(r)e^{-iω_nt}=ξI_n(r) (44). This perturbation mode signifies an infinitely slow translation between configurations σ and σ+δσ."

    This proves only that a zero mode gives a static, charge-conserving displacement of the equilibrium family. The converse needed for Eq. (42) is that a direction in (p_c^DM, p_c^NM) along which dM=dN_NM=dN_DM=0 actually solves the linearized pulsation equations at ω=0. The text bridges this gap by saying the mode 'signifies an infinitely slow translation between configurations,' i.e., it identifies the critical-curve configuration with the zero-mode configuration by assertion rather than by constructing the mode. The equivalence to be proved is therefore placed in the characterization of the critical point.

  2. other [§III.C, Step 2, after Theorem III.2]
    "Because the arbitrary time-independent variations considered in the theorem apparently include the static perturbations corresponding to radial oscillation modes, this permutation symmetry implies that if any two of the three gradients (∇M,∇N_NM, and ∇N_DM) are parallel (indicating simultaneous extrema), the third must also be parallel to them. This reduces the requirement of mutual parallelism among all three gradients to the parallelism of any two gradients."

    This is the step that identifies the OR condition in Eq. (42) with the AND condition in Eq. (43) that the proof actually targets. Theorem III.2 establishes that an equilibrium is an extremum of M under variations conserving both particle numbers; it does not by itself imply a relation between the three gradient fields at a given point of the two-dimensional pressure plane. The asserted 'permutation symmetry' is doing the work of converting (42) into (43), so the equivalence is assumed in the reduction rather than derived from the variational principle.

full rationale

The paper's numerical comparison (Figs. 4-5) is an independent check for one EoS combination and one cross-product branch, so the stability maps are not entirely fitted inputs. Likewise, the self-citations ([17], [20], [88]) are not the sole support for the main theorem. However the claimed derivation of Eq. (42) is not fully self-contained: the converse of Theorem III.1 is not proved, and the OR/AND identification relies on an asserted symmetry. These are construction-level reductions of the target equivalence to the paper's own characterization of critical points, meriting a partial-circularity score. The stability sign relation in Eq. (56) is also imported without derivation, although external references accompany the self-citation [88].

Axiom & Free-Parameter Ledger

4 free parameters · 6 axioms · 0 invented entities

The central pressures (p_c^NM, p_c^DM) are scanned independent variables of the two-fluid TOV family, not fitted parameters. The listed EoS constants (ℓ, B, m_f) and the mixing ratio are inputs chosen by hand from the cited model literature; none is fitted to observational data. No new entities are postulated — the DM components are standard literature models and the twin-star/halo configurations are numerical outputs of the two-fluid TOV system. The claims load most heavily on the domain assumptions and the unproven zero-mode/equilibrium-family correspondence, listed under axioms.

free parameters (4)
  • Holographic EoS brane separation ℓ = ℓ^{−7} = 10300 in Figs. 5a/5c and 7-9
    Sole adjustable parameter of the AdS/QCD double-polytropic NM EoS (Eq. 8, from [74]); hand-selected as a reference EoS (M_max ≈ 1.85 M⊙), not fitted to the stability result.
  • Bosonic DM coupling B = 0.08√λ4 (m/GeV)² = 0.1 (Figs. 3a,4,5a,7-9), 0.03 (Fig. 5b), 0.0012/0.0022 (Table I)
    Free parameter of the SIDM EoS (Eq. 10, from [34,79]); varied by hand to span halo masses 55-101 M⊙.
  • Fermionic DM mass m_f = 0.5 in Figs. 5c/5d (units not stated)
    Free parameter of the free Fermi gas DM EoS (Eq. 9); chosen for comparison with Kain [19].
  • Mixing ratio r = p_c^DM / p_c^NM = 0.1 and 0.9 in Table I and Fig. 10
    Hand-chosen central-pressure ratio defining the DM admixture in the 'prediction prospects' demonstration; M_max varies from ~55 to ~101 M⊙ across the table.
axioms (6)
  • domain assumption Mixed-star decoupling ansatz: the two fluids interact only gravitationally; each energy-momentum tensor is conserved separately; each fluid is barotropic and zero-temperature (p_I = p_I(ε_I)).
    §II.A-B: defines scenario II (mixed core) and yields the multi-fluid TOV system (Eq. 7) and the decoupled pulsation equations (Eqs. 32-37); the entire paper operates within it.
  • domain assumption Spectral theory of the two-fluid pulsation operator: real eigenvalues, complete mode expansion, existence of a fundamental mode, and the energy-sign correspondence sgn(δ²M) = sgn(ω₀²).
    §III.C 'Mathematical Foundations' and Eq. (56): explicitly 'adopted without further proof' from [82,87,88]; this is what turns the constrained-extremum picture into a dynamic stability statement.
  • domain assumption Mode continuity: ω²_n vary continuously in the central-pressure plane, so a zero-mode curve separates stable and unstable regions.
    §III.C Preliminary Remark 2, justified by time-reversal invariance [86]; required for the critical curve to be the stability boundary rather than an isolated locus.
  • domain assumption Zero-mode/equilibrium-family correspondence: a vanishing-frequency radial mode is a static perturbation connecting two neighboring equilibrium configurations with the same M, N_NM, N_DM.
    §III.C Theorem III.1 (Eqs. 44-49): asserted via the conservation argument; this is the multi-fluid analog of the turning-point principle proven for one fluid in [84,85] and is precisely the step the proof does not fully establish (see weakest_assumption).
  • standard math Isentropic first law used in the variational derivation: δε_I = (ε_I+p_I)δn_I/n_I.
    Appendix A, Eq. (A6); standard zero-temperature thermodynamics, consistent with the barotropic EoS assumption.
  • standard math General relativity with standard TOV construction: spherical symmetry, Schwarzschild-like coordinates, ideal-fluid energy-momentum tensor.
    §II.B Eqs. (1)-(6); background framework of the paper.

reviewed 2026-08-03 · how reviews work

0 comments
Cite this review

Pith. "Pith review of Stability Boundary of Neutron-Dark Matter Mixed Stars." pith.science (2026). https://pith.science/paper/5IEDJ5S3

@misc{pith2026251201641,
  author       = {Pith},
  title        = {Pith review of: Stability Boundary of Neutron-Dark Matter Mixed Stars},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/5IEDJ5S3}},
  note         = {Machine review of arXiv:2512.01641}
}
Share X Bluesky LinkedIn Reddit HN
read the original abstract

Regarding the stability of two-fluid star models, we rigorously prove the equivalence between the emergence of a zero-frequency radial oscillation mode and the static critical-curve criterion for mixed stars, after briefly reviewing the hybrid star case. This establishes a sufficient-condition relation between two independently developed stability criteria. Although this connection has often been implicitly assumed in previous studies of mixed stars, it has rarely been demonstrated explicitly. Our derivation can be extended to general multi-fluid systems. As an illustrative example, we consider dark matter-admixed neutron star models and show that their stability boundary differs from that of single-fluid stars. In this case, stable configurations form a surface in the three-dimensional parameter space spanned by central pressure, mass, and radius, giving rise to a class of stable mixed stars. This class includes "twin stars" with identical masses and radii but distinct internal compositions and structures. These results provide a useful framework for interpreting compact star observations and for constraining dark matter properties through astrophysical measurements.

Figures

Figures reproduced from arXiv: 2512.01641 by Kilar Zhang, Si-Man Wu, Tian-Shun Chen, Xiao-Ding Zhou.

Figure 1
Figure 1. Figure 1: FIG. 1: Schematic diagram: two different possible [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2: Schematic diagram: the stability turning point [PITH_FULL_IMAGE:figures/full_fig_p003_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3: Single fluid stability analysis. Left panel (Fig.3a) shows fundamental eigenvalue curves as functions of [PITH_FULL_IMAGE:figures/full_fig_p013_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4: Radial oscillation analysis results for the mixed star model with EoS of Holographic NM ( [PITH_FULL_IMAGE:figures/full_fig_p013_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5: Critical curves for different EoS combinations. The black curves indicate multiple stability boundaries. The [PITH_FULL_IMAGE:figures/full_fig_p014_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6: The critical curves corresponding to the [PITH_FULL_IMAGE:figures/full_fig_p015_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: FIG. 7: The two-dimensional stable region from the pressure parameter space is mapped onto a three-dimensional [PITH_FULL_IMAGE:figures/full_fig_p016_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: FIG. 8: In both panels, [PITH_FULL_IMAGE:figures/full_fig_p017_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: FIG. 9: Deviation from the BTM criteria in mixed stars. [PITH_FULL_IMAGE:figures/full_fig_p017_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: FIG. 10: Prediction prospects for admixed star models displayed through the [PITH_FULL_IMAGE:figures/full_fig_p018_10.png] view at source ↗

discussion (0)

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

Forward citations

Cited by 1 Pith paper

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

  1. Rotating Fermion-Boson Stars in $R$-squared Gravity

    gr-qc 2026-07 conditional novelty 6.0

    R-squared gravity enlarges the equilibrium domain of rotating fermion-boson stars and raises static and Keplerian maximum masses relative to GR while remaining compatible with current compact-object constraints.

Reference graph

Works this paper leans on

101 extracted references · 7 canonical work pages · cited by 1 Pith paper

  1. [1]

    Numerator

    Single Fluid Chandrasekhar, in his 1964 paper, rigorously estab- lished the pulsation equation for radial oscillations in single-fluid stars based on radial perturbations and vari- ational principles [2], and it is widely applied in many works such as [81]. For the need of unifying notation, we briefly review the main steps and key equations in the deriva...

  2. [2]

    The oscillation equations for such a system have been previously derived in the literature and applied to other mixed models (see, e.g., [[18][19][83]])

    Multi-Fluid For a compact star composed of multiple ideal fluids that are in the absence of any non-gravitational coupling, the pulsation equations for a single fluid can be general- ized to the multi-fluid case. The oscillation equations for such a system have been previously derived in the literature and applied to other mixed models (see, e.g., [[18][1...

  3. [3]

    BTM Criteria for Single-Fluid The BTM criteria is a method for assessing the radial stability of general relativistic star models. The key prin- ciples of this criterion are (we apply its generalized ver- sion in [54], that when passing through each extremum of theM-Rcurve, the direction can go along either in- creasing or decreasing the core pressure): •...

  4. [4]

    Critical Curve Criteria for Multi-Fluid The stability analysis of multi-fluid compact stars can be significantly simplified through the application of the critical curve method. This approach, originally devel- oped in [8] for determining the stability of boson-fermion stars, is in fact applicable to any spherically symmetric ideal multi-fluid star model ...

  5. [5]

    The connection between the existence of static solu- tions and the critical curve relies on conserved quan- tities

    corresponds toξ In (r). The connection between the existence of static solu- tions and the critical curve relies on conserved quan- tities. Recalling the derivation of the two-fluid equi- librium and oscillation equations, the definitions of the physical quantities and their perturbations in terms of 11 Lagrangian displacements originate from the covarian...

  6. [6]

    (b) CorrespondingM-Rrelation curves

    Single Fluid As shown in Fig.3, we show the eigenvalue-central pres- sure relations for single-fluid models with three different 13 (a) Single fluid eigenvalue curves as functions of central configuration for different EoS. (b) CorrespondingM-Rrelation curves. FIG. 3: Single fluid stability analysis. Left panel (Fig.3a) shows fundamental eigenvalue curves...

  7. [7]

    To numerically identify the stability boundary stati- cally, we need to precisely identify where the gradients become parallel

    Multi-Fluid The numerical results in Fig.4a provides precise squared frequencies of fundamental radial oscillation us- ing Bosonic DM (B 4 = 0.1) and Holographic NM (ℓ −7 = 10300), thus in addition to identifying stable regions, it also reveals the trend of fundamental frequency variation, showing non-monotonic behavior. To numerically identify the stabil...

  8. [8]

    sta- ble group

    These re- lations obtained through single-parameter variations in mixed star models no longer satisfy the BTM criteria. Furthermore, the deviation from the BTM criteria can be quantitatively analyzed, as demonstrated in Fig.9b. This panel plots the relative separation between the cen- tral pressure of the last stable point,p last-stable, and that of the m...

  9. [9]

    Our mixed star model demonstrates the ability to cover a wide range ofM-Rarea by adjusting the EoS parameters

    Prediction prospects Collecting the latestM-Rcredible window from ob- servational data of compact stars on theM-Rstable groups can impose constraints on DM EoS. Our mixed star model demonstrates the ability to cover a wide range ofM-Rarea by adjusting the EoS parameters. For in- stance, Fig.8a covers the mass and radius values of data PSR J0437-4715 from ...

  10. [10]

    The maximum mass point (circle marker) and the last stable configuration point (star marker), determined by the critical curve method, do not coincide. (b) The deviation is described as the relative error between the central pressure of the last stable point and that of the maximum mass point, and it is a function of the DM central pressurep c

  11. [11]

    sta- ble group

    This explicitly shows that the turning point of stability deviates from the mass peak, invalidating the simple BTM criteria for these sequences. FIG. 9: Deviation from the BTM criteria in mixed stars. V. CONCLUSION In this work, we have conducted a comprehensive in- vestigation into the stability properties and macroscopic characteristics of DM-NM mixed s...

  12. [12]

    B. P. Abbott, R. Abbott, and e. a. Abbott (LIGO Scien- tific Collaboration and Virgo Collaboration), Gw170817: Observation of gravitational waves from a binary neutron star inspiral, Phys. Rev. Lett.119, 161101 (2017)

  13. [13]

    Chandrasekhar, The Dynamical Instability of Gaseous Masses Approaching the Schwarzschild Limit in General Relativity, Astrophys

    S. Chandrasekhar, The Dynamical Instability of Gaseous Masses Approaching the Schwarzschild Limit in General Relativity, Astrophys. J.140, 417 (1964)

  14. [14]

    J. M. Bardeen, K. S. Thorne, and D. W. Meltzer, A cat- alog of methods for studying the normal modes of radial pulsation of general-relativistic stellar models, The As- trophysical Journal145, 505 (1966)

  15. [15]

    MISHRA, S

    H. MISHRA, S. MISRA, P. PANDA, and B. PARIDA, Neutron matter – quark matter phase transition and quark star, International Journal of Modern Physics E 02, 547–563 (1993)

  16. [16]

    S. K. Ghosh, S. C. Phatak, and P. K. Sahu, Hybrid stars and quark hadron phase transition in chiral colour dielec- tric model, Zeitschrift f¨ ur Physik A Hadrons and Nuclei 352, 457–466 (1995)

  17. [17]

    KHADKIKAR, A

    S. KHADKIKAR, A. MISHRA, and H. MISHRA, Con- finement, quark matter equation of state and hybrid stars, Modern Physics Letters A10, 2651–2663 (1995)

  18. [18]

    Di Clemente, M

    F. Di Clemente, M. Mannarelli, and F. Tonelli, Reliable description of the radial oscillations of compact stars, Physical Review D101, 10.1103/physrevd.101.103003 (2020)

  19. [19]

    Henriques, A

    A. Henriques, A. R. Liddle, and R. Moorhouse, Stabil- ity of boson-fermion stars, Physics Letters B251, 511 (1990)

  20. [20]

    Zhang, B.-A

    N. Zhang, B.-A. Li, J. Zhang, W. Shen, and H. Zhang, Il- luminating dark matter admixed in neutron stars with si- multaneous mass–radius constraints, Symmetry17, 1669 (2025)

  21. [21]

    Panotopoulos, How compact stars challenge our view about dark matter (2019), arXiv:1910.02279 [gr-qc]

    G. Panotopoulos, How compact stars challenge our view about dark matter (2019), arXiv:1910.02279 [gr-qc]

  22. [22]

    Rezaei, Double dark matter admixed neutron star, In- ternational Journal of Modern Physics D27, 1950002 (2018)

    Z. Rezaei, Double dark matter admixed neutron star, In- ternational Journal of Modern Physics D27, 1950002 (2018)

  23. [23]

    X. D. Wang, B. Qi, G. L. Yang, N. B. Zhang, and S. Y. Wang, Possible maximum mass of dark matter existing in compact stars based on the self-interacting fermionic model, International Journal of Modern Physics D28, 1950148 (2019)

  24. [24]

    Mukhopadhyay and J

    P. Mukhopadhyay and J. Schaffner-Bielich, Quark stars admixed with dark matter, Physical Review D93, 10.1103/physrevd.93.083009 (2016)

  25. [25]

    Xiang, W.-Z

    Q.-F. Xiang, W.-Z. Jiang, D.-R. Zhang, and R.-Y. Yang, Effects of fermionic dark matter on properties of neutron stars, Physical Review C89, 10.1103/physrevc.89.025803 (2014)

  26. [26]

    Leung, M.-C

    S.-C. Leung, M.-C. Chu, and L.-M. Lin, Dark- matter admixed neutron stars, Physical Review D84, 10.1103/physrevd.84.107301 (2011)

  27. [28]

    Zhang, G.-Z

    K. Zhang, G.-Z. Huang, J.-S. Tsao, and F.-L. Lin, Gw170817 and gw190425 as hybrid stars of dark and nu- clear matter, The European Physical Journal C82, 1 (2022)

  28. [29]

    Kain, Radial oscillations and stability of multiple- fluid compact stars, Phys

    B. Kain, Radial oscillations and stability of multiple- fluid compact stars, Phys. Rev. D102, 023001 (2020), arXiv:2007.04311 [gr-qc]

  29. [30]

    Kain, Dark matter admixed neutron stars, Phys

    B. Kain, Dark matter admixed neutron stars, Phys. Rev. D103, 043009 (2021), arXiv:2102.08257 [gr-qc]

  30. [31]

    Chen, X.-D

    T.-S. Chen, X.-D. Zhou, and K. Zhang, Equivalence of stability criteria for multi-fluid stars, to appear (2025)

  31. [32]

    R. C. Tolman, Static solutions of einstein’s field equations for spheres of fluid, Phys. Rev.55, 364 (1939)

  32. [33]

    J. R. Oppenheimer and G. M. Volkoff, On massive neu- tron cores, Physical Review55, 374 (1939)

  33. [34]

    Panotopoulos and I

    G. Panotopoulos and I. Lopes, Dark matter effect on re- alistic equation of state in neutron stars, Physical Review D96, 083004 (2017)

  34. [35]

    Leung, M.-C

    S.-C. Leung, M.-C. Chu, and L.-M. Lin, Dark-matter admixed neutron stars, Physical Review D—Particles, Fields, Gravitation, and Cosmology84, 107301 (2011)

  35. [36]

    Kolb,The early universe(CRC press, 2018)

    E. Kolb,The early universe(CRC press, 2018)

  36. [37]

    F. D. Steffen, Dark-matter candidates: Axions, neutrali- 22 nos, gravitinos, and axinos, The European Physical Jour- nal C59, 557 (2009)

  37. [38]

    G¨ uver, A

    T. G¨ uver, A. E. Erkoca, M. H. Reno, and I. Sarcevic, On the capture of dark matter by neutron stars, Journal of Cosmology and Astroparticle Physics2014(05), 013

  38. [39]

    Davoudiasl, Gravitationally induced dark matter asymmetry and dark nucleon decay, Phys

    H. Davoudiasl, Gravitationally induced dark matter asymmetry and dark nucleon decay, Phys. Rev. D88, 095004 (2013)

  39. [40]

    A. R. Zentner and A. P. Hearin, Asymmetric dark matter may alter the evolution of very low-mass stars and brown dwarfs, Phys. Rev. D84, 101302 (2011)

  40. [41]

    Tulin and H.-B

    S. Tulin and H.-B. Yu, Dark matter self-interactions and small scale structure, Physics Reports730, 1–57 (2018)

  41. [42]

    L. G. van den Aarssen, T. Bringmann, and C. Pfrom- mer, Is dark matter with long-range interactions a solu- tion to all small-scale problems of Λ cold dark matter cosmology?, Physical Review Letters109, 10.1103/phys- revlett.109.231301 (2012)

  42. [43]

    Tulin, H.-B

    S. Tulin, H.-B. Yu, and K. M. Zurek, Resonant dark forces and small-scale structure, Phys. Rev. Lett.110, 111301 (2013)

  43. [44]

    J. L. Feng, M. Kaplinghat, and H.-B. Yu, Halo-shape and relic-density exclusions of sommerfeld-enhanced dark matter explanations of cosmic ray excesses, Physi- cal Review Letters104, 10.1103/physrevlett.104.151301 (2010)

  44. [45]

    Zhang, L.-W

    K. Zhang, L.-W. Luo, J.-S. Tsao, C.-S. Chen, and F.-L. Lin, Dark stars and gravitational waves: Topical review, Results in Physics , 106967 (2023)

  45. [46]

    Aprile, K

    E. Aprile, K. Arisaka, F. Arneodo,et al.(XENON100 Collaboration), First dark matter results from the xenon100 experiment, Phys. Rev. Lett.105, 131302 (2010)

  46. [47]

    Ahmed, D

    Z. Ahmed, D. S. Akerib, S. Arrenberg,et al.(CDMS Collaboration), Results from a low-energy analysis of the cdms ii germanium data, Phys. Rev. Lett.106, 131302 (2011)

  47. [48]

    C. E. Aalseth, P. Barbeau, N. Bowden, B. Cabrera- Palmer, J. Colaresi, J. Collar, S. Dazeley, P. De Lurgio, J. E. Fast, N. Fields,et al., Results from a search for light-mass dark matter with a p-type point contact ger- manium detector, Physical Review Letters106, 131301 (2011)

  48. [49]

    Bernabei, P

    R. Bernabei, P. Belli, F. Cappella, R. Cerulli, C. Dai, A. d’Angelo, H. He, A. Incicchitti, H. Kuang, J. Ma, et al., First results from dama/libra and the combined results with dama/nai, The European Physical Journal C56, 333 (2008)

  49. [50]

    Leung, M.-C

    S.-C. Leung, M.-C. Chu, and L.-M. Lin, Dark-matter ad- mixed neutron stars, Phys. Rev. D84, 107301 (2011)

  50. [51]

    Goldman and S

    I. Goldman and S. Nussinov, Weakly interacting mas- sive particles and neutron stars, Phys. Rev. D40, 3221 (1989)

  51. [52]

    Kouvaris, Wimp annihilation and cooling of neutron stars, Phys

    C. Kouvaris, Wimp annihilation and cooling of neutron stars, Phys. Rev. D77, 023006 (2008)

  52. [53]

    Bertone and M

    G. Bertone and M. Fairbairn, Compact stars as dark mat- ter probes, Phys. Rev. D77, 043515 (2008)

  53. [54]

    de Lavallaz and M

    A. de Lavallaz and M. Fairbairn, Neutron stars as dark matter probes, Phys. Rev. D81, 123521 (2010)

  54. [55]

    Kouvaris and P

    C. Kouvaris and P. Tinyakov, Can neutron stars con- strain dark matter?, Phys. Rev. D82, 063531 (2010)

  55. [56]

    Brito, V

    R. Brito, V. Cardoso, and H. Okawa, Accretion of dark matter by stars, Phys. Rev. Lett.115, 111301 (2015)

  56. [57]

    Cermeno, M

    M. Cermeno, M. ´A. P´ erez-Garc ´ ıa, and J. Silk, Fermionic light dark matter particles and the new physics of neu- tron stars, Publications of the Astronomical Society of Australia34, e043 (2017)

  57. [58]

    M. I. Gresham and K. M. Zurek, Asymmetric dark stars and neutron star stability, Physical Review D99, 083008 (2019)

  58. [59]

    Ellis, G

    J. Ellis, G. H¨ utsi, K. Kannike, L. Marzola, M. Raidal, and V. Vaskonen, Dark matter effects on neutron star properties, Phys. Rev. D97, 123007 (2018)

  59. [60]

    Deliyergiyev, A

    M. Deliyergiyev, A. Del Popolo, L. Tolos, M. Le Delliou, X. Lee, and F. Burgio, Dark compact objects: An exten- sive overview, Phys. Rev. D99, 063015 (2019)

  60. [61]

    A. E. Nelson, S. Reddy, and D. Zhou, Dark halos around neutron stars and gravitational waves, Journal of Cos- mology and Astroparticle Physics2019(07), 012

  61. [62]

    Bramante and N

    J. Bramante and N. Raj, Dark matter in compact stars (2025), arXiv:2307.14435 [hep-ph]

  62. [63]

    Alford, M

    M. Alford, M. Braby, M. Paris, and S. Reddy, Hybrid stars that masquerade as neutron stars, The Astrophys- ical Journal629, 969 (2005)

  63. [64]

    Issifu, A

    A. Issifu, A. Konstantinou, P. Thakur, and T. Frederico, Rotating proto-neutron stars admixed with mirror dark matter: A two fluid approach (2025), arXiv:2507.20823 [astro-ph.HE]

  64. [65]

    Zhang and F.-L

    K. Zhang and F.-L. Lin, Constraint on hybrid stars with gravitational wave events, Universe6, 231 (2020)

  65. [66]

    S. L. Pitz and J. Schaffner-Bielich, Generating ultra- compact neutron stars with bosonic dark matter (2024), arXiv:2408.13157 [astro-ph.HE]

  66. [67]

    Kumar and H

    A. Kumar and H. Sotani, Stability analysis of two-fluid neutron stars featuring twin star and ultradense config- urations (2025), arXiv:2509.03862 [astro-ph.HE]

  67. [68]

    Ivanytskyi, V

    O. Ivanytskyi, V. Sagun, and I. Lopes, Neutron stars: New constraints on asymmetric dark matter, Physical Review D102, 10.1103/physrevd.102.063028 (2020)

  68. [69]

    M. I. Gresham and K. M. Zurek, Asymmetric dark stars and neutron star stability, Physical Review D99, 10.1103/physrevd.99.083008 (2019)

  69. [70]

    J. M. Bardeen, K. S. Thorne, and D. W. Meltzer, A Cat- alogue of Methods for Studying the Normal Modes of Radial Pulsation of General-Relativistic Stellar Models, Astrophys. J.145, 505 (1966)

  70. [71]

    Shahrbaf, D

    M. Shahrbaf, D. R. Karkevandi, A. Ayriyan, and S. Typel, Observational probes of the neutron star equa- tion of state with hyperons, bosonic dark matter, and quark matter, arXiv preprint arXiv:2402.18686 (2024)

  71. [72]

    Shahrbaf, P

    M. Shahrbaf, P. Thakur, and D. R. Karkevandi, Prob- ing strange dark matter throughf-mode oscillations of neutron stars with hyperons and quark matter (2025), arXiv:2510.08115 [nucl-th]

  72. [73]

    Mukhopadhyay, D

    S. Mukhopadhyay, D. Atta, K. Imam, D. N. Basu, and C. Samanta, Compact bifluid hybrid stars: hadronic matter mixed with self-interacting fermionic asymmet- ric dark matter, The European Physical Journal C77, 10.1140/epjc/s10052-017-5006-3 (2017)

  73. [74]

    V. R. Pandharipande, C. J. Pethick, and D. G. Raven- hall, Equation of state and neutron star models, Nuclear Physics AA498, 313c (1989), proceedings of the Inter- national Nuclear Physics Conference, Harrogate, UK

  74. [75]

    Akmal, V

    A. Akmal, V. Pandharipande, and D. Ravenhall, Equa- tion of state of nucleon matter and neutron star struc- ture, Physical Review C58, 1804 (1998)

  75. [76]

    R. B. Wiringa, V. Fiks, and A. Fabrocini, Equation of 23 state for dense nucleon matter, Physical Review C38, 1010 (1988)

  76. [77]

    M¨ uther, M

    H. M¨ uther, M. Prakash, and T. Ainsworth, The nuclear symmetry energy in relativistic brueckner-hartree-fock calculations, Physics Letters B199, 469 (1987)

  77. [78]

    Mueller and B

    H. Mueller and B. D. Serot, Relativistic mean-field theory and the high-density nuclear equation of state, Nuclear Physics A606, 508 (1996)

  78. [79]

    B. D. Lackey, M. Nayyar, and B. J. Owen, Observational constraints on hyperons in neutron stars, Phys. Rev. D 73, 024021 (2006)

  79. [80]

    Goriely, N

    S. Goriely, N. Chamel, and J. Pearson, Skyrme-hartree- fock-bogoliubov nuclear mass formulas: Crossing the 0.6 mev accuracy threshold with microscopically deduced pairing, Physical review letters102, 152503 (2009)

  80. [81]

    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, Nuclear Physics A635, 231 (1998)

Showing first 80 references.

This paper was first reviewed by deepseek-v4-flash on August 3, 2026.