{"id":"244435e8-adc8-412f-9ac3-ee64a9b6e739","arxiv_id":"2504.20825","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"First consistent numerical simulations of dark-matter-admixed neutron star mergers show that dark matter cores favor black hole collapse, halos form common envelopes, and standard tidal deformability calculations fail for halos.","lead":"For the first time, researchers simulated neutron star mergers where the stars contain dark matter, using fully consistent initial conditions. They found that dark matter arranged as a central core or a surrounding halo changes the merger outcome, and that standard calculations of the stars' tidal deformation fail for halo cases.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Halo tidal-deformability claim rests on a fitted Lambda_est, not an independent effective-tidal calculation; the conclusion that prior GW DM constraints may be invalid is underdetermined.","rationale":"The paper is a genuinely useful first step: it constructs constraint-solved two-fluid initial data and performs new NR evolutions, and it is transparent about its limitations. The reader's conditional verdict is appropriate. My primary concern differs from the reader's stated weakest assumption (the non-interacting, zero-temperature fermionic DM model). Instead, the weakest load-bearing step is the tidal-deformability interpretation in Sec. III F and the conclusion that previous GW-based DM exclusions may be invalid. That conclusion depends on the unverified assertion that the standard two-fluid Lambda calculation overestimates the effective tidal response for halo configurations. The paper demonstrates only that choosing a smaller Lambda improves agreement with a phenomenological waveform model; it does not derive that smaller value from first principles. Because Lambda is a tunable waveform parameter, the improved agreement is not independent evidence. This is a correctness risk rather than an internal inconsistency, and the paper explicitly calls for refined tidal-deformability calculations, so the appropriate outcome remains conditional acceptance. The concern also reinforces the reader's point (i) about the hand-picked Lambda_est and is consistent with the absence of a clean convergence order in Appendix C. A focused two-fluid tidal-Love calculation, or an independent NR-based Lambda measurement, would settle whether the halo result is a real physical effect or an artifact of comparing with an inappropriate effective model. For these reasons I keep the reader's conditional verdict unchanged.","tokens_in":32131,"tokens_out":6367,"duration_ms":77836,"concrete_test":"Independently compute the effective tidal deformability of the M2405H and M2805H initial configurations by solving the coupled two-fluid relativistic Love equations with proper junction conditions at the BM-DM interface and at the DM surface, or by extracting the tidal phase from the early-inspiral NR waveforms with Lambda treated as a free matched-filter parameter. If a correct calculation yields Lambda_eff close to Lambda_out (~2908 and ~901), the paper's proposed resolution fails and the claim that prior GW constraints may need revisiting is unsupported; if it yields values near Lambda_est (~810 and ~340), the concern is resolved. Additionally, run an R3 halo simulation to verify that the dephasing and f2 results are not resolution artifacts.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The most externally consequential claim is the GW/tidal-deformability conclusion in Sec. III F: halo configurations show large phase disagreement with IMRPhenomXAS NRTidalv3 when the two-fluid Love number is integrated to the halo edge, while using a hand-selected Lambda_est roughly three times smaller than Lambda_out restores agreement. The paper interprets this as evidence that the standard two-fluid tidal-deformability calculation is inadequate and that previously excluded fDM-mDM parameter space may be allowed. The load-bearing step is not the simulation itself but the assertion that the correct effective tidal deformability of an isolated halo-admixed star is far below Lambda_out. No independent computation of this effective tidal deformability is provided; Lambda_est is chosen to reduce dephasing. Because Lambda is a free parameter in the waveform model, tuning it to improve agreement does not validate the physical interpretation. If the true effective tidal deformability is close to Lambda_out, the phase mismatch instead indicates that a single-Lambda quasi-circular model is inadequate for two-fluid extended halos, which would not imply that previous GW DM constraints are invalid. The conclusion is also not protected by convergence: only R1 and R2 are available for M2405H, and Appendix C reports no clear convergence order, with the adopted error band taken from R2 versus R1.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents numerical-relativity simulations of dark-matter-admixed binary neutron star mergers, using constraint-solved initial data from the sgrid code and dynamical evolutions with BAM. Dark matter is modeled as a non-interacting, zero-temperature fermionic gas coupled to baryonic matter only through gravity, with baryonic matter described by the SLy4 equation of state. Six configurations are simulated: two total masses (2.4 and 2.8 solar masses), with DM-free, DM-core (3% DM fraction), and DM-halo (0.5% DM fraction) morphologies. The main reported results are that DM-core systems retain a compact central DM structure and can lead to more compact remnants or prompt collapse, while DM-halo systems develop a common DM envelope embedding the binary; that halo configurations suppress baryonic ejecta; and that gravitational-wave dephasing comparisons with IMRPhenomXAS NRTidalv3 show large disagreement for halo configurations when the standard two-fluid tidal deformability Lambda_out is used, whereas a hand-chosen Lambda_est roughly three times smaller restores agreement. Based on this, the authors suggest that the standard two-fluid tidal-deformability calculation is inadequate for extended halos and that previously excluded dark-matter parameter space may be allowed by GW170817 and GW190425.","tokens_in":32369,"tokens_out":5634,"duration_ms":65245,"significance":"If the central claims hold, this would be a meaningful step: it is, to my knowledge, the first set of BNS merger simulations with DM-admixed initial data that satisfy the Einstein constraint equations and use a two-fluid treatment with tabulated microphysical equations of state. The morphological findings—DM halos forming a common envelope, DM cores remaining distinct and tidally deformed, and the preservation of DM morphology in the remnant—come directly from the simulations and are the paper's strongest contribution. The paper also provides useful quantitative data on DM ejecta masses, post-merger angular-velocity profiles, and the gravitational-wave l=2,|m|=1 mode, and it makes waveform data available on Zenodo. However, the externally consequential claim about tidal deformability and previous GW-based DM constraints is underdetermined, and the numerical convergence support for that claim is incomplete. The paper is therefore best viewed as a promising initial exploration whose quantitative and phenomenological conclusions need further substantiation.","major_comments":[{"comment":"The conclusion that halo configurations invalidate the standard two-fluid tidal-deformability calculation and that previously excluded DM parameter space may be allowed rests on the ad hoc choices Lambda_est = 810 and 340, which are selected to reduce the NR-vs-model dephasing rather than derived from an independent computation. Because Lambda is a free parameter of IMRPhenomXAS NRTidalv3, tuning it to improve agreement does not validate the physical interpretation; if the true effective tidal deformability is close to Lambda_out, the mismatch instead signals that a single-Lambda quasi-circular model is inadequate for two-fluid extended halos. The claim needs support from an independent effective-tidal calculation or an explicit demonstration of why integrating the two-fluid Love equation to R_out overestimates the tidal response. As written, the abstract statement that 'previous conclusions' may be invalid is not justified by the evidence presented.","section":"Section III F, Fig. 9"},{"comment":"The convergence analysis shows no well-defined convergence order, and the adopted error band is the R2-R1 phase difference. For the two halo runs that carry the main tidal claim, Table II shows only R1 and R2 resolutions (M2405H and M2805H have no R3 entries), so the 'error band' in Fig. 9 for these configurations is a single two-resolution difference without an R3 check. Similarly, Table IV lists f2 for M2405H only at R2. The quantitative dephasing comparison for the halo configurations is therefore not validated at the same level as the DM-core runs, and the error estimate should be treated as provisional until a third resolution is available.","section":"Appendix C and Fig. 9"},{"comment":"The DM-free baselines have substantially larger initial separations than the DM-admixed runs: d_in = 53.05 km for M2400 and 56.02 km for M2800, compared with approximately 47 km for the DM-admixed configurations. Because a longer inspiral increases numerical diffusion before merger, the quantitative ejecta suppression factors ('factor of approximately 10' and 'suppression by factor 100') and the BH mass differences cannot be cleanly attributed to dark matter. The authors acknowledge this issue for the BH mass, but the ejecta claims are stated without the same caveat. Matched-separation baselines or a quantitative estimate of the diffusion-driven mass loss are needed before these factors can be taken at face value.","section":"Table I and Sections III A, III E"},{"comment":"The statement that scenarios with a dark matter core 'tend to exhibit a higher probability of prompt collapse' is not supported by the two simulated total masses. Only the M28 series collapses, and within that series all configurations except the low-resolution DM-free R1 run collapse regardless of DM morphology at the higher resolutions. With two mass points and one resolution-dependent survivor, the data cannot establish a probability trend; the conclusion should be restricted to these specific configurations or softened to a qualitative statement about the simulated cases.","section":"Abstract and Section III C"}],"minor_comments":[{"comment":"The text contains a typo: 'more accuratly' should read 'more accurately'.","section":"Section III F"},{"comment":"The f2 values for M243C at R2 (2.990 kHz) and R3 (3.156 kHz) differ by approximately 5%; the statement of 'consistent peak frequencies' should either quantify this spread or explain why it does not affect the conclusions, especially given the later caveat about resolution dependence.","section":"Table IV, Section III F"},{"comment":"The meaning of 'adapted coordinates of the XCTS system' for the initial separation d_in is unclear; a definition of how this coordinate separation relates to the actual orbital separation would help readers interpret Table I.","section":"Section III A"},{"comment":"The caption uses the word 'Lambda' in the text for Λ; the notation should be unified with the symbol used elsewhere in the paper.","section":"Figure 9 caption"},{"comment":"The atmosphere parameters f_atm and f_th are specified, but no test of sensitivity to these choices is presented; a brief justification or a reference to a convergence test would strengthen the ejecta analysis, which depends on the density floor.","section":"Section II E"}],"recommendation":"major_revision","confidential_remarks":"The paper is within scope for the journal and the novelty claim of constraint-solved DM-admixed BNS initial data appears well supported by prior work by Rüter et al. The main hazard is the underdetermined tidal-deformability conclusion and the uneven resolution coverage for halo configurations. I would ask the authors to provide either an independent effective-tidal calculation or an explicit argument for why Lambda_out is the wrong quantity, and to add an R3 run for at least one halo configuration before the GW-constraint implications are stated in the abstract."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"You should know two things. First, this is exactly what it claims to be: the first set of numerical-relativity binary neutron star merger simulations with two-fluid dark-matter admixed stars, using constraint-solved initial data, tabulated EOSs for both fluids, and DM halos that extend beyond the baryonic star. That is a real technical advance. Second, the paper's most consequential claim—that halo configurations invalidate the standard two-fluid tidal deformability calculation and could reopen DM parameter space excluded by GW170817/GW190425—is not supported by the evidence here. The agreement is restored by choosing Lambda_est by hand, roughly three times smaller than Lambda_out. No independent effective tidal deformability is computed. So the observation is qualitatively interesting but qualified: a single-Lambda quasi-circular waveform model does not capture an extended two-fluid halo, yet the paper does not show that the correct effective Lambda is far below Lambda_out.\n\nWhat it does well: the ID setup is careful, the runs are expensive and reported with transparency, and the appendices are honest—the convergence study openly says there is no clear convergence order. The morphological results—DM cores give more compact remnants, DM halos form a common envelope, cores stay core-like, halos stay diffuse—are new outputs and look qualitatively robust across the resolutions available. The angular-velocity analysis, despite the center-of-rotation uncertainty, is a nice addition. The ejecta trends, especially the factor-of-ten suppression for halos, are plausible.\n\nSoft spots, in proportion. The main one is the fitted Lambda_est: because Lambda is a free parameter in IMRPhenomXAS NRTidalv3, tuning it to reduce the dephasing does not validate the physical interpretation. The stress-test note has this right. The DM-free baselines also use a larger initial separation, so the comparison is not apples-to-apples for merger time and mass diffusion; the authors acknowledge it, but it weakens the quantitative BH mass and ejecta comparisons. The convergence order is absent, as they say, and the f2 comparison rests on resolution-dependent numbers. All of this is addressable; none of it undermines the core capability claim.\n\nWho is this for: anyone working on DM in compact objects and NR waveform modeling. It deserves serious peer review as a technical first step. If I were handling it, I would ask for a revision that either computes an independent effective tidal deformability for halos or drops the GW-constraint-revival claim, and that sharpens the DM-free comparison caveats. My own verdict would be accept after major revision.","headline":"First constraint-solved NR simulations of DM-admixed BNS mergers, with a real technical advance and a robust core/halo phenomenology; but the tidal-deformability conclusion rests on a hand-fitted Lambda_est and should not be used to reopen DM parameter space yet.","tokens_in":32971,"tokens_out":3115,"would_cite":true,"duration_ms":30940,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Dark matter's spatial layout—core or halo—decides how a neutron star merger ends and whether standard gravitational-wave models break.","keywords":["binary neutron star mergers","dark matter admixed neutron stars","numerical relativity","gravitational waves","tidal deformability","fermionic dark matter","prompt collapse","common envelope"],"falsifier":"Rerun the halo configuration (0.5% dark matter, 0.17 GeV particle mass) with a self-interacting or warm dark-matter equation of state and check whether the common envelope still forms and whether the standard two-fluid tidal deformability still overpredicts the gravitational-wave phase; if the mismatch disappears, the morphology-driven conclusion is an artifact of the non-interacting model.","tokens_in":31938,"feed_emoji":"🌌","tokens_out":8220,"duration_ms":83792,"temperature":0.7,"pith_summary":"This paper reports the first numerical-relativity simulations of binary neutron star mergers that contain dark matter and begin from constraint-solved initial data, meaning the two stars are set up as a consistent two-fluid gravitational system rather than superposed single-star snapshots. Modeling dark matter as a non-interacting fermionic gas, the paper finds that the spatial arrangement of dark matter—a compact core inside each star or a diffuse halo surrounding it—changes the merger outcome. Cores make the remnant more compact and promote prompt collapse to a black hole; halos merge early into a common envelope embedding the whole binary. The gravitational waves from halo mergers also disagree with analytical waveform models when the tidal deformability is computed in the standard two-fluid way, suggesting that earlier gravitational-wave-based exclusions of dark matter parameter space may need revisiting.","feed_headline":"Dark matter's shape decides neutron star merger fate","feed_subtitle":"First constraint-solved simulations show cores prompt collapse while halos form common envelopes","key_machinery":"The argument rests on a two-fluid general-relativistic framework in which baryonic matter and dark matter are separate perfect fluids whose energy-momentum tensors are conserved independently and coupled only through gravity. The initial data are generated as quasi-equilibrium, constraint-solved binaries using an extended conformal thin-sandwich formulation of the Einstein constraint equations, with the dark matter treated as a non-interacting, zero-temperature Fermi gas of spin-1/2 particles; two particle masses (1 GeV and 0.17 GeV) reproduce the two morphologies, a dense core and an extended halo. The evolution follows both fluids with ideal general-relativistic hydrodynamics, and the tidal deformability is computed by integrating Love's equation to the outermost radius. The machinery's job is to allow the two fluids to interact self-consistently through spacetime curvature throughout inspiral, merger, and post-merger, so that morphology-driven differences in dynamics and gravitational waves can be attributed to the dark matter structure rather than to inconsistent initial data.","core_discovery":"The paper's central claim is that dark matter morphology, not just its mass fraction, controls the merger dynamics of dark-matter-admixed neutron star binaries. With a 3% dark matter core, the baryonic stars are more compact and the post-merger remnant has higher central density; in the heavier 2.8 solar mass case this leads to prompt collapse to a black hole with a slightly more massive black hole than in dark-matter-free runs. With a 0.5% dilute halo, the two halos come into contact before the baryonic stars, forming a common dark matter envelope that embeds the binary, and the final remnant keeps a halo-like distribution. The tidal-deformability calculation in a standard two-fluid framework, integrating to the outermost radius, gives values roughly three times larger than the waveform-consistent estimate for halos, and with the lower estimate the numerical-relativity waveforms match an analytical waveform model within the error band. The paper concludes that previous constraints on fermionic dark matter from gravitational-wave observations may need to be revisited.","pith_inferences":["If the halo tidal-deformability mismatch is real, gravitational-wave searches that use two-fluid tidal deformabilities to exclude dark-matter parameter regions may be excluding configurations that are actually consistent with observed events; the exclusion regions would need to be recomputed with a halo-aware tidal deformability.","The common-envelope phase formed by halos could leave an observable imprint in the pre-merger gravitational-wave signal that is not captured by current waveform models, because the envelope changes the effective quadrupole moment and tidal response before the baryonic stars touch.","Dark matter ejected during the merger could later be re-accreted by surrounding objects, a 'dark matter recycling' channel that would modify the inferred accumulation history of neutron stars in dense dark-matter environments.","A testable extension would be to build initial data that already contain a common dark-matter envelope, which the current initial-data solver cannot represent, and check whether the envelope forms even earlier and strengthens the gravitational-wave dephasing."],"forward_implications":["If dark matter forms a dilute halo around each neutron star, mergers should show an early common-envelope phase and a suppression of baryonic ejecta by roughly an order of magnitude compared with dark-matter-free binaries.","If dark matter forms a dense core, heavier binaries are more likely to collapse promptly to a black hole, and the resulting black hole can be slightly more massive than in the dark-matter-free case.","Standard two-fluid tidal deformabilities computed to the outermost halo radius overestimate the halo's tidal effect by about a factor of three, and using a lower, waveform-consistent value brings the gravitational-wave phase into agreement with analytical models.","Dark matter ejecta masses in these mergers lie in the range 10^-6 to 10^-4 solar masses, with halos ejecting more dark matter than cores.","Post-merger angular-velocity profiles differ by dark matter morphology: dark matter cores rotate faster than the baryonic component, while halo remnants show a central plateau in baryonic angular velocity."],"supporting_citations":[{"why":"Supplies the constraint-solved quasi-equilibrium initial-data construction for two-fluid dark-matter-admixed binaries, the paper's central methodological upgrade.","marker":"[41]"},{"why":"Earlier numerical simulations of dark-matter-admixed binaries with superimposed initial data, the baseline this work improves on and compares with.","marker":"[40]"},{"why":"Provides the non-interacting fermionic dark matter model and the halo/core morphology framework used throughout.","marker":"[33]"},{"why":"Earlier constraints on fermionic dark matter from neutron star observations that the halo results suggest need revisiting.","marker":"[88]"},{"why":"Two-fluid tidal deformability calculation that the paper shows fails for extended halos, motivating the corrected estimate.","marker":"[97]"},{"why":"Quasi-universal relation for the post-merger f2 frequency that yields roughly 30% deviation for halo configurations with the standard two-fluid tidal deformability.","marker":"[109]"},{"why":"Analytical waveform model used for the dephasing comparison that quantifies the halo tidal-deformability mismatch.","marker":"[100]"}],"fun_headline_variants":["Dark matter core or halo flips merger outcome","Dark matter morphology sets merger collapse path","Dark matter cores prompt collapse, halos form envelopes","Neutron star mergers reveal dark matter's shape matters"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that dark matter inside these stars is a cold, non-interacting fermionic gas that couples to ordinary matter only through gravity; if dark matter self-interacts, has finite temperature, or has other statistics, the core-versus-halo outcomes could change or vanish.","fun_headline_variants_meta":{"raw":{"variants":["Dark matter core or halo flips merger outcome","Dark matter morphology sets merger collapse path","Dark matter cores prompt collapse, halos form envelopes","Neutron star mergers reveal dark matter's shape matters"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000865,"raw_usage":{"total_tokens":3771,"prompt_tokens":984,"completion_tokens":2787,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":600,"completion_tokens_details":{"reasoning_tokens":2727}},"tokens_in":600,"tokens_out":2787,"duration_ms":20238,"temperature":1.0,"reasoning_tokens":2727,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-16T05:19:10.915802+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Rerun the halo configuration (0.5% dark matter, 0.17 GeV particle mass) with a self-interacting or warm dark-matter equation of state and check whether the common envelope still forms and whether the standard two-fluid tidal deformability still overpredicts the gravitational-wave phase; if the mismatch disappears, the morphology-driven conclusion is an artifact of the non-interacting model.","supporting_citations":[],"review_version":1}