{"id":"304f35e4-b511-4049-b0af-b827ce83538a","arxiv_id":"2506.20736","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"A Bayesian analysis of dark matter admixed neutron stars finds that astrophysical data constrain only the dark matter fraction, which is governed by the hadronic equation of state stiffness.","lead":"This paper fits three dark matter parameters in neutron stars that also contain fermionic dark matter in their cores, using Bayesian inference against NICER and gravitational wave data. It concludes that current data only constrain the dark matter mass fraction, and that this fraction is set mainly by the stiffness of the hadronic equation of state.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The inferred fDM constraints depend on the core-confinement assumption, but the MD prior overlaps the halo branch and the transition also depends on Cvd and fDM; the central claim is therefore not established for the full two-fluid configuration space.","rationale":"The reader's weakest-assumption analysis identified the core-confinement assumption as the key vulnerability, and the paper itself flags the halo possibility in Section I. I agree that this is the most load-bearing concern. The central claim is not incorrect for the core-confined branch, but the manuscript does not establish that the inferred fDM constraints survive when the physically allowed halo branch is included, and part of the MD prior overlaps the stated halo regime. This does not require rejection, because the analysis is coherent within its stated assumption, but it does require the authors to either justify that the relevant prior volume is core-confined or to extend the analysis to both branches. The reader's CONDITIONAL verdict is therefore appropriate, and my stress-test does not change it. I considered other potential issues, such as the overstated 'well constrained' language and the lack of code/data, but these are secondary: the fDM posteriors are indeed narrower than the priors, and reproducibility concerns do not bear directly on the physical argument. I also considered whether the 'fDM insensitive to observations' claim is contradicted by Table II, where Case I vs Case II medians differ by up to ~30% for some EoSs; this is a quantitative overstatement but is within the reported 1-sigma uncertainties and does not threaten the main conclusion. The halo issue is the one that could, if real, change the interpretation of every reported posterior.","tokens_in":17239,"tokens_out":3919,"duration_ms":51837,"concrete_test":"Repeat the Bayesian analysis with a two-fluid TOV solver that admits both branch types. Specifically, integrate Eq. (5) with independent boundary conditions P_HM(R_HM)=0 and P_DM(R_DM)=0, and accept a configuration if it is stable, regardless of whether R_DM <= R_HM (core-confined) or R_DM > R_HM (halo). Then rerun the Case II likelihood for the same priors. If the fDM posterior median or 1-sigma interval shifts by more than the quoted errors, or if MD becomes constrained, the central claim is branch-dependent. A cheaper complementary check is to map the core/halo phase boundary across the MD–Cvd–fDM prior using the criterion in Refs. [19,39,55] and report the fraction of prior volume that yields halos; if that fraction is non-negligible, the omitted branch must be included before the claim can stand.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central claim is that current observations constrain only the DM mass fraction fDM, not MD or Cvd, and that fDM is set mainly by hadronic EoS stiffness. This claim is derived entirely from two-fluid TOV solutions in which DM is confined to the stellar core. However, Section I explicitly states that DM particles with masses of a few hundred MeV form an extended halo around the NS, while more massive particles are core-confined, and the transition depends on MD, the coupling, and the DM fraction (Refs. [19,39,55]). The Bayesian prior is MD = 500–3000 MeV, so the lower part of the prior lies in or near the halo regime. Moreover, because the transition is not a hard mass cutoff, some configurations with MD above 500 MeV and small Cvd will also produce halos. For halo configurations, the two-fluid TOV equations (Eq. 5) must be solved with the DM fluid extending beyond the hadronic surface, and the resulting M–R and Lambda–M relations differ from core-confined solutions. The reported posteriors, e.g., fDM ~ 2–5% in Cases I–III, are therefore conditional on one branch of the solution space. Since the halo branch is not explored, the conclusion that observations constrain fDM rather than MD or Cvd, and that fDM is governed by hadronic stiffness, is not shown to hold for the full model class. This is a scope limitation, not an internal inconsistency, but it is load-bearing because the paper's main message is stated without this caveat.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies dark matter admixed neutron stars (DMANS) with fermionic dark matter interacting only gravitationally with hadronic matter. Using RMF equations of state for both sectors and the two-fluid TOV equations, the authors perform a Bayesian analysis to constrain the dark matter particle mass MD, the vector coupling-to-mass ratio Cvd, and the dark matter mass fraction fDM. Three treatments of astrophysical data are considered: Case I uses a maximum-mass cut plus radii from PSR J0030+0451 and GW170817 tidal deformability; Case II uses full NICER mass-radius posteriors for PSR J0740+6620 and PSR J0030+0451 through KDEs; Case III uses mock data with reduced uncertainties. The analysis is repeated for two realistic EoSs (BITSH-E and BITSH-I), with and without high-density speed-of-sound uncertainties, and for stiffer EoSs obtained by setting the vector self-interaction term to zero. The central claim is that current observations mainly constrain fDM, while MD and Cvd remain poorly constrained, and that fDM is largely insensitive to the likelihood treatment and high-density EoS uncertainties, being instead determined primarily by the stiffness of the hadronic EoS.","tokens_in":17535,"tokens_out":4952,"duration_ms":56529,"significance":"If the central claim holds, the paper provides a useful, observationally grounded statement about what current neutron-star data can tell us about fermionic dark matter: essentially only the DM mass fraction, not the particle mass or coupling, and with a strong degeneracy with the hadronic EoS stiffness. The study's strengths include the use of recently developed nuclear EoSs constrained by finite nuclei, heavy-ion collisions, and astrophysical data; a nested-sampling Bayesian framework; and a systematic comparison of three likelihood constructions, multiple EoSs, and high-density uncertainty treatments. The result that fDM below roughly 5% is favored for realistic EoSs is a falsifiable statement relevant to future NICER and gravitational-wave observations. The main caveat, discussed in the report, is that the analysis is restricted to core-confined dark matter configurations, so the quantitative posteriors apply only to that branch of the model space.","major_comments":[{"comment":"The manuscript explicitly assumes that DM is entirely confined to the neutron star core (Section I), but the Bayesian prior on MD is 500-3000 MeV, and the text states that particles of order a few hundred MeV form an extended halo, with the core/halo transition depending on MD, Cvd, and fDM (Refs. [19,39,55]). Because the lower part of the MD prior lies in or near the halo regime, the two-fluid TOV solutions (Eq. 5) and the resulting posteriors in Tables II and III are conditional on one branch of the configuration space. The conclusion that current observations constrain only fDM, and that fDM is set primarily by hadronic stiffness, is therefore not established over the full model class; I request either a restriction of the prior to the core-confined regime or an explicit treatment or discussion of the halo branch.","section":"Sections I and III"},{"comment":"The likelihood in Eq. (13) is written as a symmetric Gaussian, and the text states that all data are assumed to follow a symmetric Gaussian distribution, despite the reported measurements being asymmetric: R1.34=12.71+1.14-1.19 km, R1.44=13.02+1.24-1.06 km, and Lambda_1.4=190+390-120. The paper does not specify how the sigma in Eq. (13) is chosen from these asymmetric errors. Since Case I is one of the three observational treatments used to support the claim that the DM fraction is insensitive to the treatment of observations, this misspecification should be corrected by using an asymmetric likelihood or by convolving with the actual posterior/KDE samples.","section":"Section III, Eq. (13)"},{"comment":"The maximum-mass constraint from PSR J0740+6620 is described as a 'stringent cut' at 2.073 +/- 0.069 M_sun, but no explicit likelihood or prior term is given for it. A hard cut with no uncertainty, a cut using only the lower bound, or a Gaussian term with the quoted 1-sigma error would lead to different accepted parameter regions, and the current text does not allow the reader to reproduce this part of the analysis. Please state explicitly how the Mmax constraint enters the likelihood or posterior.","section":"Section III, Case I"},{"comment":"The statement that 'the DM fraction is largely insensitive to how astrophysical observations are integrated' is stronger than Table II supports. For example, for BITSH-E with C=0 and HDU, the median fDM changes from 8.12% in Case I to 4.82% in Case II and 3.58% in Case III, a factor of about 2.3 between Cases I and III; for BITSH-I with C not equal 0 and HDU, the ordering of BITSH-E and BITSH-I reverses between Case I and Case II. These median shifts are within the broad 1-sigma intervals, so the qualitative conclusion may survive, but the wording should be qualified by the actual numerical spread rather than presented as a strong insensitivity.","section":"Section IV.A and Table II"}],"minor_comments":[{"comment":"The title contains a typo: 'th e' should be 'the'.","section":"Title"},{"comment":"In the text describing the finite-nuclei constraints, 'chare radii' should be 'charge radii'.","section":"Section II.A"},{"comment":"In the description of Case III, 'PSR J740+6620' should be 'PSR J0740+6620'.","section":"Section IV.A, Case III"},{"comment":"The definition of P(m|Theta) appears after Eq. (15) although it is used in Eq. (14); it should be introduced before Eq. (14).","section":"Section III, Eqs. (14)-(15)"},{"comment":"Several references have incomplete bibliographic information, for example Refs. [7] and [30] lack publication years; please standardize the reference list.","section":"References"}],"recommendation":"major_revision","confidential_remarks":"The stress-test concern about core confinement is valid and is the main reason for major revision. The paper's central message is a posterior observation about which parameter is best constrained, not a prediction from the model, and the authors are mostly careful about that distinction in the text. The symmetrized likelihood and the unspecified maximum-mass cut are technical issues that can be fixed locally. I would not reject: the scope limitation is clearly acknowledged in Section I, and the authors can respond by restricting the prior, exploring the halo branch, or substantially qualifying the conclusions."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"You should know this paper is a standard Bayesian two-fluid TOV analysis of fermionic DM-admixed neutron stars, using the authors' own RMF EoSs (BITSH-E and BITSH-I). The main message—that current astrophysical data constrain only the DM mass fraction, not the particle mass or coupling, and that the fraction tracks hadronic EoS stiffness—is consistent with prior work (Refs. [24] and [32]) and is probably right. The paper is honest about citing those antecedents, so the novelty is incremental: new EoSs, a systematic comparison of three likelihood treatments (including a mock future dataset), and a useful scan over several other RMF EoSs in Table III. What the paper does well: the Bayesian machinery is standard but applied carefully, and the demonstration that fDM is insensitive to how the observational likelihood is constructed while sensitive to hadronic stiffness is a clean, believable result. The authors also correctly note that their conclusion is about which fitted parameter is best constrained, not a prediction. The soft spots, in order of importance. First, the core-confinement assumption is load-bearing. The paper states that DM particles with masses of a few hundred MeV form an extended halo, yet the prior on MD runs from 500 MeV up to 3000 MeV. So a non-negligible part of the prior overlaps the halo regime, and the transition also depends on Cvd and fDM. The posteriors—and the claim that only fDM is constrained—are conditional on one branch of the two-fluid solution space. This is a scope limitation rather than an internal inconsistency, but the abstract states the conclusion without that caveat, which overstates the generality. Second, the language \"well constrained\" for fDM goes beyond what the numbers support: medians like 2.46(+1.49/-1.65)% have relative uncertainties of order 50–100%. \"Least poorly constrained\" would be accurate. Third, reproducibility: no code or data are released, and the speed-of-sound priors are only referenced to [61], not listed. The Gaussian symmetrization of clearly asymmetric NICER errors is a minor issue I would not press hard. Who this is for: people working on DMANS constraints, especially those wondering how the BITSH EoSs behave in this context. It is not a breakthrough, but it is a useful data point. I would send it to peer review rather than desk-reject it, but the authors should be asked to state the core-confined scope in the abstract, temper the \"well constrained\" language, and provide either code/data or at least a complete list of priors.","headline":"Plausible, competent, and incremental DMANS constraints on the authors' own EoSs; the central claim holds for core-confined DM but the abstract overstates the scope.","tokens_in":690,"tokens_out":1885,"would_cite":false,"duration_ms":46546,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Current neutron-star observations constrain the dark-matter mass fraction inside the star's core, but they leave the dark-matter particle mass and its vector coupling essentially unconstrained.","keywords":["dark matter admixed neutron stars","fermionic dark matter","two-fluid TOV equations","relativistic mean-field equation of state","Bayesian inference","NICER mass-radius constraints","tidal deformability","dark matter mass fraction"],"falsifier":"Measure the mass, radius, and tidal deformability of a neutron star whose hadronic EoS is already constrained (for example, a future NICER radius measurement of a ~1.4 $M_\\odot$ pulsar matched to a gravitational-wave event). If the measured radius falls outside the 95% band that the paper's core-confined two-fluid model predicts for the posterior $f_{\\rm DM}$ values, the claimed constraint would be contradicted; conversely, a radius that cannot be reproduced without dark matter would support it. A more direct test is to fit the same data with and without a halo component: if a halo model and a core-confined model fit the NICER and GW170817 data equally well, the paper's core-confined fraction is not uniquely determined.","tokens_in":17023,"feed_emoji":"🌌","tokens_out":13851,"duration_ms":148155,"temperature":0.7,"pith_summary":"This paper asks whether current neutron-star observations can reveal anything about dark matter that has accreted into a star and settled in its core. The authors model the star as two fluids—hadronic matter and fermionic dark matter—that interact only through gravity, and they fit the dark sector's particle mass, vector coupling, and mass fraction to pulsar mass–radius data, NICER posterior distributions, and the GW170817 tidal-deformability constraint. Their central finding is that these observations constrain the dark-matter mass fraction but leave the particle mass and coupling essentially undetermined. They also find that the inferred fraction is controlled mainly by the stiffness of the hadronic equation of state at high densities, with stiffer hadronic models allowing larger dark-matter fractions, while the way observations are included and the high-density hadronic uncertainty play a secondary role. As a result, the dark-matter mass fraction is the one dark-sector quantity that current neutron-star observations can meaningfully constrain.","feed_headline":"Neutron star data pin dark matter fraction, not its mass","feed_subtitle":"The hadronic equation of state, not the observations, decides how much dark matter a neutron star can hold.","key_machinery":"The central object is the two-fluid Tolman–Oppenheimer–Volkoff system, in which hadronic matter and dark matter each have their own conserved energy-momentum tensor and pressure gradient but share the same metric, so their only interaction is gravitational. Each fluid is described by a relativistic mean-field equation of state: the hadronic side by the BITSH-E and BITSH-I parametrizations, calibrated to finite-nuclei, heavy-ion, and neutron-star data, and the dark side by a Fermi gas of mass $M_D$ with a repulsive dark-vector coupling $C_{\\rm vd}=g_{\\rm vd}/m_{\\rm vd}$. Varying the ratio of central energy densities fixes the dark-matter mass fraction $f_{\\rm DM}$, and the two-fluid TOV solution yields mass–radius curves and tidal deformabilities that are compared with NICER and GW170817 data through a Bayesian likelihood. A speed-of-sound parametrization above $2\\rho_0$ supplies the high-density hadronic uncertainty used in a second set of fits.","core_discovery":"The paper's central claim, stated in its own terms, is that a core-confined fermionic dark-matter component in a neutron star is currently constrained only through its mass fraction $f_{\\rm DM}$, not through the particle mass $M_D$ or the dark-vector coupling ratio $C_{\\rm vd}$. Bayesian fits with the realistic BITSH-E and BITSH-I hadronic EoSs yield median $f_{\\rm DM}$ values around 2–4% (1.7–3.4% depending on likelihood case), with 2σ upper limits near 10–14%, while $M_D$ and $C_{\\rm vd}$ posteriors remain nearly flat over their priors. The authors also claim that the preferred DM fraction is set by the stiffness of the hadronic equation of state at high densities: removing the vector-meson self-interaction makes the EoS stiffer and raises median $f_{\\rm DM}$ to 4–8% for BITSH-E, and the very stiff NL3 model gives a median of about 15%, whereas the softer models stay below 5%. The way astrophysical data are inserted (mass-cut plus radii, full NICER KDE posteriors, or mock future measurements) changes $f_{\\rm DM}$ little, and the high-density speed-of-sound uncertainty has only a weak effect.","pith_inferences":["Inference: If the true dark-matter particle is light enough to form a halo rather than a core component, the core-confined fractions reported here should be read as upper limits on core dark matter; the same two-fluid machinery would need an extended component to capture the actual signature.","Inference: Because the preferred $f_{\\rm DM}$ moves with the assumed hadronic stiffness, dark-matter constraints from neutron stars should be reported jointly with the high-density hadronic EoS; a fully simultaneous Bayesian marginalization over both sectors would likely widen the $f_{\\rm DM}$ intervals beyond those in Tables II and III.","Inference: A multi-messenger test is possible: across several binary neutron-star events, the model predicts a correlation between the inferred dark-matter fraction and the tidal deformability that is different from what pure hadronic EoS variation would produce, so a gravitational-wave catalog could distinguish the two."],"forward_implications":["The median dark-matter mass fraction for the realistic BITSH-E and BITSH-I hadronic EoSs sits near 2–4%, with 2σ upper limits around 10–14%, while the fermion mass $M_D$ and coupling ratio $C_{\\rm vd}$ remain essentially unconstrained.","Switching between likelihood treatments (mass-cut plus radii, full NICER KDE posteriors, or mock future precision data) changes the inferred $f_{\\rm DM}$ only mildly, so current data are already close to what these observations can say about the dark-matter fraction.","Removing the vector-meson self-interaction makes the hadronic EoS stiffer and raises the preferred $f_{\\rm DM}$, and the very stiff NL3 model pushes the median fraction to roughly 15%, so the allowed dark-matter fraction is tied to the high-density hadronic stiffness.","Including high-density hadronic uncertainties through the speed-of-sound parametrization has a weak effect on $f_{\\rm DM}$, indicating that the fraction constraint is not driven by the high-density part of the hadronic EoS.","The star's maximum mass, canonical radius, and tidal deformability correlate mainly with $f_{\\rm DM}$, not with $M_D$ or $C_{\\rm vd}$, so the mass fraction is the dark-sector input that shapes observable properties of dark-matter admixed neutron stars."],"supporting_citations":[{"why":"supplies the relativistic mean-field hadronic EoS (BITSH-E) built from finite-nuclei, heavy-ion, and neutron-star constraints.","marker":"[53]"},{"why":"supplies the companion hadronic EoS (BITSH-I) with the alternative implicit finite-nuclei treatment.","marker":"[54]"},{"why":"provides the finite-nuclei and heavy-ion experimental constraints used to fix the low-density hadronic EoS family.","marker":"[57]"},{"why":"gives the speed-of-sound parametrization used to vary the hadronic EoS above $2\\rho_0$ in the high-density uncertainty analyses.","marker":"[61]"},{"why":"provides the GW170817 tidal-deformability constraint $\\Lambda_{1.4}=190^{+390}_{-120}$ entering the likelihood as $L_{\\rm GW}$.","marker":"[63]"},{"why":"provides the $2.073\\pm0.069\\,M_\\odot$ mass of PSR J0740+6620 used as the Case I maximum-mass cut.","marker":"[64]"},{"why":"supplies the PSR J0030+0451 radius measurement $R_{1.34}$ used in the Case I likelihood.","marker":"[65]"},{"why":"supplies the PSR J0030+0451 radius measurement $R_{1.44}$ and the NICER mass-radius posterior used in Case II.","marker":"[66]"},{"why":"provides the PSR J0740+6620 NICER mass-radius posterior used in the Case II KDE likelihood.","marker":"[67]"},{"why":"provides the simultaneous NICER mass-radius measurement of PSR J0740+6620 used alongside [67] in Case II.","marker":"[68]"}],"fun_headline_variants":["Stiff hadronic EoS raises neutron star dark matter fraction","Dark matter fraction in neutron stars set by EoS stiffness","Astrophysical data constrain only dark matter fraction","Neutron star dark matter: fraction constrained, mass free"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing assumption is that all dark matter is confined to the neutron star core, whereas the paper itself notes that sufficiently light dark-matter particles would instead form an extended halo, so a halo configuration would invalidate the core-confined fraction constraints.","fun_headline_variants_meta":{"raw":{"variants":["Stiff hadronic EoS raises neutron star dark matter fraction","Dark matter fraction in neutron stars set by EoS stiffness","Astrophysical data constrain only dark matter fraction","Neutron star dark matter: fraction constrained, mass free"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000645,"raw_usage":{"total_tokens":3061,"prompt_tokens":1138,"completion_tokens":1923,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":754,"completion_tokens_details":{"reasoning_tokens":1856}},"tokens_in":754,"tokens_out":1923,"duration_ms":13617,"temperature":1.0,"reasoning_tokens":1856,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T22:42:37.642567+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the mass, radius, and tidal deformability of a neutron star whose hadronic EoS is already constrained (for example, a future NICER radius measurement of a ~1.4 $M_\\odot$ pulsar matched to a gravitational-wave event). If the measured radius falls outside the 95% band that the paper's core-confined two-fluid model predicts for the posterior $f_{\\rm DM}$ values, the claimed constraint would be contradicted; conversely, a radius that cannot be reproduced without dark matter would support it. A more direct test is to fit the same data with and without a halo component: if a halo model and a core-confined model fit the NICER and GW170817 data equally well, the paper's core-confined fraction is not uniquely determined.","supporting_citations":[{"cited_title":"Thakur, A","cited_arxiv_id":null,"evidence_quote":"supplies the companion hadronic EoS (BITSH-I) with the alternative implicit finite-nuclei treatment."},{"cited_title":"Malik, K","cited_arxiv_id":null,"evidence_quote":"gives the speed-of-sound parametrization used to vary the hadronic EoS above $2\\rho_0$ in the high-density uncertainty analyses."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"provides the GW170817 tidal-deformability constraint $\\Lambda_{1.4}=190^{+390}_{-120}$ entering the likelihood as $L_{\\rm GW}$."},{"cited_title":"Stuart and J","cited_arxiv_id":null,"evidence_quote":"provides the $2.073\\pm0.069\\,M_\\odot$ mass of PSR J0740+6620 used as the Case I maximum-mass cut."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"supplies the PSR J0030+0451 radius measurement $R_{1.34}$ used in the Case I likelihood."},{"cited_title":"Salmi et al","cited_arxiv_id":null,"evidence_quote":"supplies the PSR J0030+0451 radius measurement $R_{1.44}$ and the NICER mass-radius posterior used in Case II."}],"review_version":1}