{"id":"2a5e93ab-3f1a-4b51-b295-ff3a0f1af340","arxiv_id":"2510.17905","paper_version":3,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":3.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"Dark-matter-admixed neutron star properties are computed across DM particle masses 0.2–1 GeV and mass fractions 0–0.8, showing that only small accumulated DM fractions (f ≤ 0.2) stay consistent with NICER and GW170817 constraints.","lead":"This paper computes how adding fermionic dark matter changes the mass, radius, and tidal properties of neutron stars, and finds that only small dark-matter fractions survive current observations. It is a parametric survey that extends earlier two-fluid dark-matter neutron-star studies to a wider grid of particle masses and two representative equations of state.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Tidal Love number calculation ignores two-fluid structure: Eq. (20) uses a single-fluid dp/dρ, which is undefined for two independent DM/BM fluids, so the Λ–M curves and f ≤ 0.2 exclusion may be invalid.","rationale":"The paper's central contribution is a systematic parameter-space survey of DM-admixed neutron stars, and the most valuable quantitative output is the exclusion of large DM fractions using NICER, pulsar masses, and GW170817 tidal constraints. The mass-radius part of the analysis uses the standard two-fluid TOV equations and is likely sound. However, the tidal deformability part, which is essential for the f ≤ 0.2 conclusion, applies the single-fluid Hinderer formalism to a two-fluid star by simply summing pressures and densities. That step is not justified: in the two-fluid model the constituent fluids have independent EOSs, so the total pressure is not a barotropic function of the total energy density, and the adiabatic index dp/dρ in Eq. (20) is not well-defined. This is an internal formal issue rather than a question of whether the zero-temperature ideal Fermi gas model is astrophysically realistic; even granting all the authors' stated assumptions, the Love number calculation may be incorrect. The reader's verdict was already CONDITIONAL, and the reader did flag the tidal-formalism issue as secondary, but the reader's weakest_assumption focused on the physical DM model assumption. My read agrees partially: the DM-model assumption is a limitation, but the tidal-formalism gap is more directly load-bearing because it threatens the numerical values that drive the main exclusion. The recommended verdict remains CONDITIONAL: the paper should be accepted only after the two-fluid tidal equations are properly derived (or shown to reduce to the single-fluid equations) and the Λ–M curves are recomputed. I do not recommend REJECT because the mass-radius results and the qualitative DM-core/halo behavior may survive, and the correct tidal calculation could plausibly reproduce the same qualitative constraints.","tokens_in":13131,"tokens_out":8743,"duration_ms":83447,"concrete_test":"Recompute the k2/Λ curves in Figs. 11 and 12 using the two-fluid tidal perturbation formalism (e.g., as in Leung, Chu & Lin 2022 [23]) with the same SLy/MPA1 EOSs, µ ∈ {0.5, 0.6, 0.7} GeV, and f ∈ {0.2, 0.4}, rather than substituting p_tot and ρ_tot into the single-fluid ODE. If Λ_1.4 changes by more than ~10% for any of these models, or if the f=0.4 curves cross the GW170817 band differently than in the paper, the f ≤ 0.2 claim is not robust. A minimal analytical version is to re-derive Eq. (22) from linearized Einstein plus separate conservation laws for the two fluids and check for extra terms proportional to (c_{s,nm}^2 − c_{s,dm}^2).","verdict_should_be":"UNCHANGED","load_bearing_attack":"Sec. 2.2 solves the standard single-fluid Hinderer equation (Eqs. 17–22) for k2 and then states that the parameters are obtained by summing the nm and dm components. This is not a valid reduction for two gravitationally coupled fluids. Eq. (20) contains (ρ+p)/(dp/dρ), i.e., the adiabatic sound speed of a single barotropic fluid. In the two-fluid model, ρ_nm and ρ_dm are independent functions with different EOSs; total pressure is not a function of total energy density, so dp_tot/dρ_tot is not defined without an additional assumption. The linearized perturbation equations require separate displacement fields and two sound speeds; the single-fluid master equation follows only after using δp = c_s^2 δρ, which does not hold when δp = c_{s,nm}^2 δρ_nm + c_{s,dm}^2 δρ_dm. The paper neither derives the two-fluid tidal equations nor cites a source justifying the application of the single-fluid formula to summed variables. The headline claim that astrophysical constraints force f ≤ 0.2 is read from the Λ–M curves in Figs. 11–12, so if the Love number calculation is incorrect, the central exclusion is not established. The paper cites Leung, Chu, and Lin (2022) [23], who do derive DM-admixed tidal deformability; using that formalism, or a fresh derivation, is the necessary check.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies neutron stars admixed with non-interacting fermionic dark matter (DM) using a two-fluid Tolman-Oppenheimer-Volkhoff (TOV) formalism. For DM particle masses μ = 0.2–1 GeV and DM mass fractions f = 0–0.8, and for the SLy and MPA1 nuclear equations of state, the authors compute mass–radius relations, density profiles, tidal Love numbers k2, and dimensionless tidal deformabilities Λ. These results are compared with constraints from ~2 M⊙ pulsars, NICER radius measurements, and the GW170817 tidal-deformability bounds. The paper concludes that (i) the DM distribution passes from a core to a halo as f increases, (ii) larger μ and f soften the effective equation of state and lower the maximum mass, and (iii) observational constraints exclude large DM fractions, with f ≤ 0.2 remaining consistent with current data.","tokens_in":13454,"tokens_out":11211,"duration_ms":82222,"significance":"If the computational framework is valid, the paper would provide a useful systematic scan of the fermionic-DM parameter space (μ, f) for DM-admixed neutron stars, connecting two nuclear equations of state to NICER and gravitational-wave constraints. Its qualitative findings—smaller maximum masses, smaller radii, and lower Λ for larger f—are plausible and consistent with earlier literature. The paper does not introduce new formalism or provide reproducible code, but it does present a clear set of figures and correctly quotes the standard ideal Fermi gas EOS and single-fluid TOV equations. The value of the work depends on the correctness of the two-fluid hydrostatic and tidal calculations; the central quantitative claim (the f ≤ 0.2 exclusion) rests on the tidal-deformability curves, where the manuscript currently has a load-bearing technical gap.","major_comments":[{"comment":"The tidal Love number calculation is not valid for a two-fluid star as written. Equation (20) contains (ρ+p)/(dp/dρ), the single-fluid barotropic sound speed. In the two-fluid model p_tot = p_nm(ρ_nm) + p_dm(ρ_dm), with ρ_tot = ρ_nm + ρ_dm, but ρ_nm(r) and ρ_dm(r) are independent functions obtained from the separate TOV equations (7)–(10). Hence p_tot is not a function of ρ_tot alone and dp_tot/dρ_tot is undefined without an additional assumption. The linearized perturbation equations require separate displacement fields for the two fluids and two sound speeds; the single-fluid Hinderer equation cannot be applied by simply summing masses, densities, and pressures. The manuscript neither derives the two-fluid tidal equations nor cites a source for the reduction. Since Figs. 11 and 12, and the resulting statement that f ≥ 0.4 is excluded by GW170817 (and f ≤ 0.2 is consistent), are based o","section":"Sec. 2.2, Eq. (20) and final paragraph"},{"comment":"There is a factor-of-2 inconsistency in the definition of dν/dr. Equation (4) gives dν/dr = 2(m + 4πr^3 p)/(r(r−2m)), which is the standard relation for the metric exponent e^ν. In the two-fluid system, Eq. (11) gives dν/dr = (m_tot + 4πr^3 p_tot)/(r(r−2m_tot)) without the factor of 2, and Eq. (7) uses dP_nm/dr = −(P_nm + ρ_nm) dν/dr. If the equations are taken literally, the hydrostatic equilibrium equations are not the standard two-fluid TOV equations, and the ν′^2 term in Eq. (20) is inconsistent by a factor of 4. This must be corrected or explicitly clarified; all numerical results, including the mass–radius curves, are affected if the code follows the equations as written.","section":"Sec. 2.1, Eqs. (4), (7), and (11)"}],"minor_comments":[{"comment":"Typo: 'IN GeV' should read 'in GeV'.","section":"Sec. 2.5"},{"comment":"The panels are labelled (a), (b), (c) in the captions, but the text refers to them as 'fig.5a' and similar. Please use standard subfigure references (e.g., 'Fig. 5a').","section":"Figs. 5 and 6"},{"comment":"The text states that four nuclear EOSs (SLy, MPA1, ENG, AP3) are used, but the admixed-star results are shown only for SLy and MPA1. Clarify why ENG and AP3 are not used in the DM-admixed analysis (or state that they are only used for the pure neutron-star comparison).","section":"Sec. 2.3"},{"comment":"The caption mentions a 'black horizontal line' indicating the 1.4 M⊙ radius span 11–13 km, but the figure appears to show a horizontal or vertical band. Please make the description of the plotted constraint explicit.","section":"Fig. 2 caption"},{"comment":"The conclusion states that 'for f = 0.2, all of the Λ−M curves match the observational limitations of GW170817', but Figs. 11–12 display only μ = 0.5, 0.6, 0.7 GeV. Indicate whether this statement holds for the full scanned μ range (0.2–1 GeV) or only for the plotted cases.","section":"Sec. 4"}],"recommendation":"major_revision","confidential_remarks":"The paper is a systematic but essentially incremental application of existing two-fluid DM-admixed neutron-star models. The principal concern is technical: the tidal deformability calculation appears to be the single-fluid Hinderer equation applied to summed two-fluid variables, which is not justified and directly affects the paper's main constraint on f. There is also a factor-of-2 inconsistency in the definition of the metric derivative between the single-fluid and two-fluid sections that needs resolution. Both issues are potentially fixable, so I recommend major revision rather than rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Two things you should know. The mass-radius half of this paper is a competent, unsurprising scan of dark-matter-admixed neutron stars; the tidal half is built on an unjustified step that likely invalidates the headline exclusion. Specifically, Sec. 2.2 takes the single-fluid Hinderer equation and says it applies to the two-fluid star by summing nm and dm contributions to mass, pressure, and energy density. That only works if total pressure is a barotropic function of total energy density, which it is not when you have two independent fluids with different EOSs. The perturbation needs two displacement fields and two sound speeds. The authors cite Leung-Chu-Lin 2022, who derive the correct two-fluid tidal formalism, but they don't use it. So the Λ-M curves in Figs. 11 and 12 are not trustworthy, and the conclusion that GW170817 forces f ≤ 0.2 is not established by this paper.\n\nWhat the paper does well: the two-fluid TOV integration is standard and the mass-radius curves are consistent with earlier work. The systematic grid over µ = 0.5–0.7 GeV and f = 0.2–0.8 for SLy and MPA1 is a useful map, and the core-to-halo transition as f grows is clearly demonstrated. The authors also deserve credit for being explicit that they use a non-interacting zero-temperature Fermi gas, even though that is a strong physical assumption that any realistic DM model would depart from.\n\nSoft spots, in order. The tidal issue is the big one. There's also a smaller mismatch: the abstract advertises µ = 0.2–1 GeV, but the admixed-star figures only show 0.5, 0.6, 0.7. And no code or data are provided, so the numerics can't be independently checked. That last one is minor; many papers in this area don't share code.\n\nWho this is for: model builders who want a quick look at how much non-interacting fermionic DM a NS can hold while staying consistent with mass and radius constraints. The M-R part of the paper could survive even if the tidal part is redone. I'd send it to a referee, but with a clear request: either derive the two-fluid tidal equations properly or use the existing formalism, and fix the abstract. If that gets done, it's a modest but usable contribution. As is, the headline exclusion should be treated as unsupported.","headline":"The M-R half is a competent scan of a known model; the tidal half is built on an unjustified single-fluid application and the f≤0.2 exclusion from Λ is not supported.","tokens_in":14056,"tokens_out":3475,"would_cite":false,"duration_ms":31096,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["04.40.Dg","95.35.+d","97.60.Jd"],"model":"deepseek-v4-flash","headline":"A two-fluid model shows that fermionic dark matter mixed into a neutron star can account for only about 20 percent of its mass before the star's radius and tidal deformability fall outside observed ranges.","keywords":["fermionic dark matter","neutron stars","two-fluid TOV","mass-radius relation","tidal deformability","Love number","ideal Fermi gas EOS","dark matter admixed neutron stars"],"falsifier":"Recompute the same mass-radius and tidal-deformability curves using any interacting fermionic dark-matter equation of state (for instance, with a repulsive self-interaction). If the f ≈ 0.2 compatibility boundary shifts by more than the observational uncertainty, the paper's exclusion is an artifact of the ideal-gas assumption; if it stays fixed, the constraint is robust.","tokens_in":12946,"feed_emoji":"🌌","tokens_out":6737,"duration_ms":57163,"temperature":0.7,"pith_summary":"This paper asks whether fermionic dark matter can hide inside neutron stars without breaking what we know about them. It models a neutron star as two fluids — ordinary nuclear matter and a zero-temperature ideal Fermi gas of dark matter — coupled only by gravity, and scans dark matter particle mass and mass fraction. The central result is that the star's mass-radius curve and its tidal deformability shift strongly with the dark matter fraction, and only small dark fractions (about f ≤ 0.2) remain compatible with current neutron star mass, radius, and tidal measurements. This matters because it turns neutron stars into a probe that can exclude whole regions of dark matter parameter space, and because the star's response depends on whether the dark matter sits as a compact core or an extended halo.","feed_headline":"Neutron star data cap dark matter share near 20 percent","feed_subtitle":"Only a small dark-matter fraction keeps neutron stars inside the observed mass, radius, and tidal-deformability windows.","key_machinery":"The central tool is the relativistic two-fluid stellar-structure (TOV) system: ordinary matter and dark matter each carry their own pressure, energy density, and mass, coupled only through a shared gravitational potential. The dark matter fluid is described by a zero-temperature ideal Fermi gas equation of state, parametrized by particle mass μ and number density. The star's total mass is split into a fraction f of dark matter. The system is completed by the standard differential equation for the quadrupolar tidal Love number, from which the dimensionless tidal deformability Λ is derived. This machinery lets the authors compute mass-radius and Λ–mass curves for each (μ, f) pair and compare t","core_discovery":"In the two-fluid picture, dark matter changes neutron star structure in a way controlled mainly by the dark matter particle mass μ and mass fraction f. For fixed compactness, small f produces a dark-matter core (dark-matter radius smaller than the nuclear-matter radius), while larger f produces a dark-matter halo surrounding the baryonic core. As μ grows, the dark-matter equation of state softens, lowering the maximum mass and shifting tidal-deformability curves toward lower masses. When matched against observations — the existence of roughly two-solar-mass pulsars, radius constraints near 1.4 solar masses, and the gravitational-wave tidal bound — only f ≈ 0.2 or less survives for the tested","pith_inferences":["If the ideal Fermi gas is replaced with a self-interacting fermionic dark-matter equation of state, the two-fluid machinery still applies, so the f ≈ 0.2 boundary can be recomputed; the direction of any shift would reveal whether the exclusion is an artifact of the non-interacting assumption.","The paper sets f by hand rather than deriving it from accretion. Combining the excluded (μ, f) region with a capture model for dark matter over the neutron star's lifetime would translate the bound into a limit on the dark-matter–nucleon cross-section — a connection the authors leave implicit.","The core-to-halo transition suggests a potentially non-monotonic response in the tidal Love number: a concentrated dark-matter core may affect tidal deformability differently than a diffuse halo even at the same total f. The paper's figures hint at this but do not isolate the two contributions."],"forward_implications":["A neutron star that accumulates more than roughly a fifth of its mass in non-interacting fermionic dark matter would be measurably different: its maximum mass drops, its radius at 1.4 solar masses shrinks, and its tidal deformability falls outside the observed band.","For a fixed overall compactness, the dark-matter component switches from a compact core (small f) to an extended halo (larger f), and the two geometries leave different fingerprints on the tidal Love number.","Because the maximum mass falls as μ increases, the existence of massive pulsars removes the high-μ corner of the (μ, f) plane.","The boundary f ≈ 0.2 is the paper's central quantitative result: below it the admixed star looks almost normal, above it the star begins to resemble a dark-matter star."],"fun_headline_variants":["Neutron stars allow only tiny dark matter shares","Dark matter fraction in neutron stars limited to 20%","Neutron star data shrink allowed dark matter fraction","Fermionic dark matter must be sparse in neutron stars","Observations restrict dark matter in neutron stars to under 20%"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The entire result depends on treating the dark matter inside the neutron star as a zero-temperature, non-interacting ideal Fermi gas that feels only gravity and has no accretion history — change that assumption and the f ≤ 0.2 ceiling may move.","fun_headline_variants_meta":{"raw":{"variants":["Neutron stars allow only tiny dark matter shares","Dark matter fraction in neutron stars limited to 20%","Neutron star data shrink allowed dark matter fraction","Fermionic dark matter must be sparse in neutron stars","Observations restrict dark matter in neutron stars to under 20%"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000805,"raw_usage":{"total_tokens":3362,"prompt_tokens":724,"completion_tokens":2638,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":468,"completion_tokens_details":{"reasoning_tokens":2557}},"tokens_in":468,"tokens_out":2638,"duration_ms":15523,"temperature":1.0,"reasoning_tokens":2557,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-04T09:08:50.935730+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Recompute the same mass-radius and tidal-deformability curves using any interacting fermionic dark-matter equation of state (for instance, with a repulsive self-interaction). If the f ≈ 0.2 compatibility boundary shifts by more than the observational uncertainty, the paper's exclusion is an artifact of the ideal-gas assumption; if it stays fixed, the constraint is robust.","supporting_citations":[],"review_version":1}