{"id":"1f291c49-a6ce-4235-a244-981a45692cc2","arxiv_id":"2507.22415","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":3.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"A proposed software algorithm for generating gravitational wave templates of dark-matter-admixed neutron stars, demonstrated for bosonic and fermionic dark matter models.","lead":"This paper presents Darksuite, a proposed software extension of the LIGO/Virgo analysis tools for generating gravitational wave templates of neutron stars that contain dark matter. The authors show such stars would produce coherent, possibly detectable phase shifts in merger waveforms, and their workflow chains stellar structure calculations into waveform generation.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Detectability claim rests on an inapplicable sinusoidal heuristic, and Figure 3 uses an inconsistent fermionic density scale; the claimed bosonic-versus-fermionic resolvability is therefore unverified.","rationale":"The reader's CONDITIONAL verdict is well calibrated: the algorithmic claim cannot be executed without code, and the detectability analysis is not a complete calculation. My stress-test identifies the single most checkable load-bearing weakness as the disconnect between the Section V sinusoidal heuristic and the actual monotonic tidal phase drift of the computed waveforms, compounded by the inconsistent fermionic density scale between the text and Figure 3. I do not elevate the two-fluid Love-number interface concern to the primary issue, because for these EOSs each fluid reaches zero pressure with zero density, so the y-equation with the summed source in Eq. (10) can be integrated without a singular interface; the absence of a derivation is a documentation gap rather than a demonstrated error. Still, because the code is not released and no benchmark against published dark-matter-admixed Λ values is given, the conditional status remains appropriate. If the matched-filter test shows the actual mismatch exceeds the threshold, the central significance claim would be supported; if not, the paper's motivation is weakened. Neither outcome moves the verdict away from CONDITIONAL at this stage.","tokens_in":15910,"tokens_out":16804,"duration_ms":223658,"concrete_test":"Regenerate the two waveforms of Figure 3 using the text value ρFℏ³ = 1.9×10⁻⁴ GeV⁴ (and also 1.4×10⁻⁴ for comparison), compute the frequency-domain phase difference ΔΨ(f) between the bosonic and fermionic templates at the GW170817 masses and SNR ≈ 30, and evaluate the noise-weighted mismatch relative to 1/(2SNR²) ≈ 5.6×10⁻⁴. Only if the mismatch exceeds this threshold does the Section V heuristic support the detectability claim; the same run resolves whether the caption/text density-scale discrepancy changes the conclusion.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section V derives a detectability threshold from a sinusoidal phase modulation model, Ψ(f) = Ψ0 sin(2πfT), giving Ψ0 ≳ 1/(SNR |sin(2πf0T)|). But the waveforms in Figure 3 are generated with NRTidalv3 tidal phase corrections, which produce a systematic, monotonic inspiral phase drift, not a sinusoidal modulation. The paper never extracts Ψ0 from Figure 3, never computes the noise-weighted mismatch between the bosonic and fermionic waveforms, and never shows the actual accumulated phase exceeds the stated threshold. The conclusion that 'phase modulations can accumulate coherently over the inspiral and may become detectable' is therefore not connected to the computed signals. Additionally, Figure 3's caption lists ρFℏ³ = 1.4×10⁻⁴ GeV⁴ while Section IV and Figure 1 use 1.9×10⁻⁴ GeV⁴; if Figure 3 was produced with the different scale, the bosonic-versus-fermionic phase difference could be a parameter artifact rather than a particle-statistics effect. The reader's weakest-assumption focus on the two-fluid Love-number interface treatment is less decisive: for the EOSs used, each fluid reaches p = ρ = 0 continuously, so a single y-integration with the summed source in Eq. (10) is a plausible standard extension; the real threat is the unreleased implementation plus the internal inconsistency in the one comparative result the paper highlights.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper introduces Darksuite, a proposed Python-based extension of LALSuite for modeling gravitational-wave signals from dark-matter-admixed neutron stars. The core algorithm solves the two-fluid Tolman-Oppenheimer-Volkoff equations for nuclear matter (BSk22) plus either self-interacting bosonic or ideal fermionic dark matter, computes the dimensionless tidal deformability via the Hinderer equation with a two-fluid modified source term, and generates waveforms using the NRTidalv3 tidal phase correction. It presents mass-radius and tidal deformability-mass sequences for several dark-matter fractions, a waveform comparison between bosonic and fermionic dark matter at f = 0.8, and a heuristic detectability discussion based on a sinusoidal phase-modulation model.","tokens_in":16184,"tokens_out":8359,"duration_ms":100762,"significance":"If the numerical implementation is correct, Darksuite would provide a useful template-generation tool for a class of exotic-compact-object searches, and the paper correctly identifies an existing gap in the LALSuite waveform infrastructure. The qualitative behavior of the M-R and Lambda-M curves is consistent with prior work on dark-matter-admixed neutron stars, and the authors are transparent about the exploratory nature of the proposal. However, the paper's central results are not yet fully validated: the two-fluid tidal deformability integration lacks a stated interface treatment, the one comparative waveform figure uses a dark-matter density scale inconsistent with the rest of the paper, and the detectability claim rests on a sinusoidal heuristic that is not connected to the computed monotonic tidal phase drift. The paper is better viewed as an algorithm description than as a demonstrated detection capability.","major_comments":[{"comment":"The two-fluid tidal deformability computation is not fully specified for configurations in which one fluid ends at a smaller radius than the other (the DM-core and DM-halo cases of Figure 1). The paper integrates a single y-perturbation equation with the summed two-fluid source term Q(r) to R = max(R_NM, R_DM), but it does not state how Q(r) is evaluated in the shell where one fluid has already reached zero pressure and density, nor whether any junction condition is imposed at that interface. Since the Love number k2 in Eq. (8) depends on y at the outer radius, and every downstream quantity (Lambda, phase shift in Eq. (3), and the comparison in Figure 3) inherits this dependence, the validity of the results for hybrid configurations is not established. Please specify the treatment of the interface region and, ideally, validate the two-fluid integration against a dedicated multi-fluid tidal code or release the Darksuite code so this step can be reproduced.","section":"Section II.C, Eq. (10); Section III (integration stopping criteria)"},{"comment":"There is an internal inconsistency in the dark-matter density scale used for the fermionic case. Section IV and Figure 1 quote rho_F hbar^3 = 1.9 x 10^-4 GeV^4, while the Figure 3 caption quotes rho_F hbar^3 = 1.4 x 10^-4 GeV^4. If Figure 3 was generated with the latter value, the bosonic-versus-fermionic phase difference may reflect the different density scales rather than the particle-statistics difference; if it was generated with 1.9 x 10^-4 GeV^4, the caption is wrong. In either case, the highlighted comparative result in Figure 3 is ambiguous and must be corrected, and the waveforms should be regenerated or the caption fixed and the analysis repeated with a consistent value.","section":"Section IV, Figure 3 caption vs Section IV/Figure 1"},{"comment":"The detectability argument is not connected to the waveforms actually computed. The sinusoidal phase-modulation model Psi(f) = Psi0 sin(2 pi f T) and the resulting threshold Psi0 >~ 1/(SNR |sin(2 pi f0 T)|) describe an oscillatory phase modulation, but the NRTidalv3 tidal phase corrections used in Figure 3 produce a monotonic, growing phase drift with frequency. The paper neither extracts Psi0 from the computed phase difference nor evaluates the noise-weighted mismatch between the bosonic and fermionic waveforms, so the claim in Section VI that 'phase modulations can accumulate coherently over the inspiral and may become detectable' is not supported by the presented quantitative analysis. Please replace the sinusoidal heuristic with a concrete overlap calculation for the actual waveforms, or substantially weaken the detectability conclusion.","section":"Section V and Section VI (detectability)"}],"minor_comments":[{"comment":"The caption contains typographical errors: 'fermonic' should be 'fermionic', and 'pannel' should be 'panel'.","section":"Figure 3 caption"},{"comment":"The notation for the density scales is inconsistent: 'rho_B hbar' and 'rho_F hbar' should be 'rho_B hbar^3' and 'rho_F hbar^3' to match the units GeV^4.","section":"Section IV, Figure 1 caption"},{"comment":"The phrase 'delivering ng a maximum mass' is a typo and should read 'delivering a maximum mass'.","section":"Section IV"},{"comment":"The phrase 'the inspiral phase phase of the gravitational waveform' contains a duplicated word and should be corrected.","section":"Section III"},{"comment":"Equation (12) is garbled in presentation; the denominator and exponential terms appear to be missing brackets (e.g., 'exp[a5(xi - a6) + 1}' should likely be '(exp[a5(xi - a6)] + 1)^-1'). Please re-typeset the BSk22 functional so it is unambiguous.","section":"Section II.D.1, Eq. (12)"},{"comment":"The statement that 'the integration ends when the pressure of either fluid approaches zero' is ambiguous for two-fluid configurations; clarify that after one fluid's pressure vanishes, the integration continues for the remaining fluid until its own pressure vanishes at its surface.","section":"Section III (integration stopping criteria)"},{"comment":"The manuscript describes Darksuite as a 'proposed extension' but does not provide a code repository or pseudocode. Given that the two-fluid tidal integration is a central and nonstandard step, releasing the code or providing detailed pseudocode would substantially improve reproducibility and verifiability.","section":"General"}],"recommendation":"major_revision","confidential_remarks":"The paper is essentially an algorithm/software description rather than a new physical result, and the authors are upfront about the proposal being exploratory. The main scientific value would be in providing a validated tool; currently the validation is incomplete due to the unspecified two-fluid interface treatment, the inconsistency in Figure 3, and the unquantified detectability claim. If the authors fix these points, the paper could become acceptable as a methods paper, but the bar for that should include either a code release or a detailed enough algorithmic specification to allow independent reproduction."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: the paper builds a plausible pipeline for dark-matter-admixed neutron star waveforms, but the pipeline is not shipped and the detectability argument does not actually apply to the waveforms it computes.\n\nWhat is genuinely new is the demonstration of a workflow: two-fluid TOV solutions, interpolation over compactness and dark matter fraction, and injection into NRTidalv3 via LALSuite to produce waveforms. The equations are all imported from prior work, but the combination is a practical contribution. The M-R and Λ-M curves look physically reasonable and reproduce known qualitative behavior. The authors are also honest in several places: they call the GW170817 comparison qualitative because spin is ignored, they concede the inverse-problem solver needs refinement, and they admit the detectability analysis is not a complete calculation.\n\nThe main problem is that Darksuite itself is not released, so the central claim — that you can generate and analyze these waveforms — cannot be checked. That alone would block acceptance in my view. The detectability section is also disconnected from the computed signals: it derives a threshold using sinusoidal phase modulation, but the waveforms in Figure 3 are generated with NRTidalv3, which produces a monotonic tidal phase drift, not a sinusoidal modulation. The paper never computes the actual mismatch or accumulated phase against its own threshold, so the conclusion that the bosonic versus fermionic difference is detectable is unverified. There is also an internal inconsistency in the one comparative result that matters: Figure 3's caption gives ρF ℏ³ = 1.4×10⁻⁴ GeV⁴, while Section IV and Figure 1 use 1.9×10⁻⁴ GeV⁴. If Figure 3 was made with a different scale, the displayed phase difference could be a parameter artifact rather than a particle-statistics effect. The reader's worry about the two-fluid tidal interface conditions is less decisive to me: for these EOSs each fluid reaches p = ρ = 0 continuously, so integrating the single y-equation with the summed source is a plausible extension, though without released code it remains unverified.\n\nThis deserves a serious referee but should not be accepted until the code is released, the parameter inconsistency is fixed, a 0% DM baseline is shown in Figure 3, and the detectability claim is replaced with a matched-filter or Fisher-matrix estimate on injected signals. The physics chain is standard and the authors are honest, so it is a reasonable conditional.","headline":"A useful-sounding pipeline for DM-admixed NS waveforms, but the code is absent, the detectability argument doesn't match the computed signals, and a key figure has an inconsistent parameter.","tokens_in":16850,"tokens_out":3487,"would_cite":false,"duration_ms":35791,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["04.30.-w","95.35.+d","97.60.Jd"],"model":"deepseek-v4-flash","headline":"The paper argues that adding a dark-matter component to neutron stars changes their tidal deformability enough to alter the gravitational-wave phase of inspiraling binaries, and that the difference between bosonic and fermionic dark…","keywords":["dark matter-admixed neutron stars","tidal deformability","gravitational waves","two-fluid TOV equations","bosonic dark matter","fermionic dark matter","LALSuite","waveform modeling"],"falsifier":"Recompute $k_2$ for the $f = 0.8$, $C = 0.08$ bosonic and fermionic configurations with a solver that enforces the proper junction conditions at each fluid surface (at $R_{\\rm NM}$ and $R_{\\rm DM}$ separately), and compare the resulting $\\Lambda$ values to the Darksuite-interpolated values; if the shifted $\\Lambda$ brings the bosonic and fermionic phase differences below the mismatch threshold of $1/(2 \\, \\mathrm{SNR}^2)$, the central claim of resolvability collapses.","tokens_in":15545,"feed_emoji":"🌌","tokens_out":5994,"duration_ms":61924,"temperature":0.7,"pith_summary":"The paper proposes Darksuite, a software extension that would let gravitational-wave analysis pipelines produce templates for neutron stars laced with dark matter. Its central claim is that the tidal deformability—and therefore the gravitational-wave phase evolution during inspiral—depends measurably on both the fraction of dark matter inside the star and on whether that dark matter is bosonic or fermionic. The authors integrate the two-fluid Tolman–Oppenheimer–Volkoff equations with a modified tidal perturbation equation, build an interpolated surface of Love numbers, and feed the result into standard LALSuite waveform models. They show waveforms for a GW170817-like binary whose bosonic and fermionic dark-matter versions accumulate a coherent, oscillatory phase difference over the inspiral. If the calculation is right, gravitational-wave observations of binary neutron stars could serve as an indirect probe of dark-matter microphysics.","feed_headline":"Neutron-star waves could betray dark matter's particle nature","feed_subtitle":"A proposed LALSuite extension computes how bosonic versus fermionic dark matter shifts the inspiral waveform.","key_machinery":"The load-bearing object is the two-fluid generalization of the relativistic stellar-structure and tidal-response equations: the two-fluid TOV system (Eqs. 4–6) plus the first-order $y$-perturbation equation (Eq. 9) with the two-fluid quadrupole source $Q(r)$ (Eq. 10). Integrating these outward to $R = \\max(R_{\\rm NM}, R_{\\rm DM})$ yields $y_R$, which feeds the Love-number formula (Eq. 8) and the dimensionless tidal deformability $\\Lambda = \\frac{2}{3} k_2 / C^5$. That $\\Lambda$ enters the NRTidalv3 phase correction (Eq. 3), and the resulting phase shift is embedded in IMRPhenomPv2 to produce waveforms. The interpolation of $(C, f, k_2)$ over a bank of TOV solutions is what lets the pipeline evaluate arbitrary dark-matter fractions without reintegrating the structure equations.","core_discovery":"The paper's central claim is that a user-supplied dark-matter fraction and particle statistics (bosonic or fermionic) produce a distinct tidal-deformability surface $\\Lambda(M)$ and, through the NRTidalv3 phase correction, gravitational waveforms that differ from pure-nuclear templates and from each other. The authors demonstrate this with a two-fluid TOV solver using the BSk22 equation of state for nuclear matter and either a self-interacting bosonic or an ideal fermionic dark-matter equation of state, with density scales chosen so pure dark-matter stars are solar-mass-scale. A key displayed result is the $f = 0.8$, $C = 0.08$ comparison in Figure 3, where the bosonic and fermionic waveforms go in and out of phase over the inspiral; the authors argue from a heuristic mismatch criterion that phase modulations above about $0.03$ radians would be resolvable at the GW170817 signal-to-noise ratio. The paper frames Darksuite as a first step toward including dark matter in standard gravitational-wave data analysis rather than as a finished parameter-estimation study.","pith_inferences":["Not claimed by the paper: a decisive validation would be to run the same pipeline with a multi-fluid tidal solver that imposes junction conditions at each fluid surface; if $k_2$ shifts, the central bosonic-versus-fermionic comparison may change.","The chosen density scales imply a maximum-mass difference between the bosonic and fermionic cases (about $2.7 M_\\odot$ versus $2.5 M_\\odot$), suggesting a mass-radius observable independent of tides, such as in merger remnant properties.","The paper's heuristic threshold implies that next-generation detectors with SNR well above 30 could resolve phase modulations below $0.03$ radians, extending sensitivity to lower dark-matter fractions.","The GW170817 comparison is qualitative because spins are ignored; a full parameter-estimation run with Darksuite templates would be needed to determine whether dark matter is actually favored over ordinary equation-of-state variation."],"forward_implications":["If a neutron star in a detected binary contains a non-negligible dark-matter admixture, standard templates will mis-estimate the tidal contribution to the phase, biasing recovered masses and radii.","The bosonic-versus-fermionic difference at $f = 0.8$ is large enough that, under the paper's SNR heuristic, current detectors could distinguish the two microphysical models from the inspiral alone.","The interpolated $(C, f, k_2)$ surfaces make it feasible to include dark-matter-admixed stars in parameter-estimation pipelines without rerunning TOV solvers for each sample.","The comparison of mass-radius curves with GW170817 posteriors shows dark-matter-admixed configurations can occupy the observationally allowed region, so such stars are not excluded by current constraints.","The framework extends to other exotic components or modified-gravity variants by adding new equation-of-state classes, so the same machinery generalizes beyond the two dark-matter models tested."],"supporting_citations":[{"why":"Supplies the two-fluid dark-matter-admixed neutron-star framework and the bosonic/fermionic density scales that make pure dark-matter stars solar-mass-scale.","marker":"[7]"},{"why":"Provides the closed-form tidal phase correction (Eq. 3) that Darksuite adopts for the gravitational-wave phase shift.","marker":"[14]"},{"why":"Defines the two-fluid stellar radius as the maximum of the dark-matter and nuclear-matter radii, the integration endpoint used here.","marker":"[18]"},{"why":"Gives the tidal Love number formula and the $y$-equation that the two-fluid modification extends.","marker":"[19]"},{"why":"Supplies the Brussels-Montreal BSk22 unified equation of state used for the nuclear-matter component.","marker":"[20]"},{"why":"Gives the self-interacting bosonic dark-matter equation of state used for the bosonic case.","marker":"[21]"},{"why":"Gives the parametric fermionic gas equation of state used for the fermionic case.","marker":"[25]"},{"why":"Provides the NRTidalv3 model whose tidal phase corrections are used to generate the waveforms.","marker":"[34]"},{"why":"Supplies the IMRPhenomPv2 waveform model used to embed the tidal phase shift into the strain.","marker":"[35]"}],"fun_headline_variants":["Dark matter's fingerprint in neutron-star gravitational waves","Algorithm reveals dark matter's particle type in neutron stars","Gravitational waves differentiate bosonic and fermionic dark matter","New framework models dark-matter-admixed neutron stars","Probing dark matter's nature via neutron-star inspiral signals"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the two-fluid tidal deformability is computed correctly as implemented, which requires that integrating the $y$-equation to the outer radius with no explicit interface conditions at the surface where one fluid ends is a valid prescription; if $k_2$ is wrong for these mixed configurations, every downstream $\\Lambda$ and the key bosonic-versus-fermionic waveform comparison are wrong.","fun_headline_variants_meta":{"raw":{"variants":["Dark matter's fingerprint in neutron-star gravitational waves","Algorithm reveals dark matter's particle type in neutron stars","Gravitational waves differentiate bosonic and fermionic dark matter","New framework models dark-matter-admixed neutron stars","Probing dark matter's nature via neutron-star inspiral signals"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00023,"raw_usage":{"total_tokens":1531,"prompt_tokens":1046,"completion_tokens":485,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":662,"completion_tokens_details":{"reasoning_tokens":406}},"tokens_in":662,"tokens_out":485,"duration_ms":5971,"temperature":1.0,"reasoning_tokens":406,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T11:45:43.039797+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Recompute $k_2$ for the $f = 0.8$, $C = 0.08$ bosonic and fermionic configurations with a solver that enforces the proper junction conditions at each fluid surface (at $R_{\\rm NM}$ and $R_{\\rm DM}$ separately), and compare the resulting $\\Lambda$ values to the Darksuite-interpolated values; if the shifted $\\Lambda$ brings the bosonic and fermionic phase differences below the mismatch threshold of $1/(2 \\, \\mathrm{SNR}^2)$, the central claim of resolvability collapses.","supporting_citations":[{"cited_title":"Two-fluid dark matter ad- mixed neutron stars","cited_arxiv_id":null,"evidence_quote":"Defines the two-fluid stellar radius as the maximum of the dark-matter and nuclear-matter radii, the integration endpoint used here."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the Brussels-Montreal BSk22 unified equation of state used for the nuclear-matter component."},{"cited_title":"Vásquez Flores, Alessandro Parisi, Chian-Shu Chen, and Germán Lugones","cited_arxiv_id":null,"evidence_quote":"Gives the self-interacting bosonic dark-matter equation of state used for the bosonic case."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the IMRPhenomPv2 waveform model used to embed the tidal phase shift into the strain."}],"review_version":1}