{"id":"223bfde3-3bd0-4dc2-b339-ce06d0e0be68","arxiv_id":"1908.05150","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"Rotating white dwarfs with dark matter cores produce gravitational-wave bursts whose peak-amplitude ratios can reveal the dark matter mass and rotation, allowing dark matter cores above 0.03 solar masses to be inferred.","lead":"This paper simulates the collapse of rotating white dwarfs that contain a compact core of dark matter, and shows the dark matter leaves a distinctive mark in the gravitational waves emitted at the moment of bounce. The authors find that measuring the ratios of the wave's three biggest peaks could reveal both the dark matter mass and the inner core's rotation rate.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"EOS-robustness test in Fig. 11 / Appendix B was run only for MDM=0, so the claimed 0.03 Msun inference threshold is not yet supported for DM-admixed collapse.","rationale":"The reader's stationary-DM-assumption concern is legitimate, but the EOS test gap is more directly load-bearing for the precise headline claim. The conclusion 'Even within the uncertainties of nuclear matter equation of state...' is a quantitative robustness statement; its evidence base is Fig. 11's light shaded regions, which Appendix B does not produce for DM-admixed models. The full text contains an explicit limitation in Section 3.3 that h2 can change by ~30% with EOS and electron-capture rate variations, so the separation used to set 0.03 Msun could be eroded by the very uncertainties the abstract claims to tolerate. I therefore partially agree with the reader: the stationary-DM issue is a real modeling caveat that could shift the calibrated relations, but the missing DM x EOS grid is the single check that would settle the stated threshold. The proposed test is inexpensive relative to a full two-fluid DM treatment and directly targets the claim. Since the concern is a support gap rather than a demonstrated contradiction, the appropriate verdict remains CONDITIONAL, as in the reader's assessment.","tokens_in":21470,"tokens_out":5871,"duration_ms":65836,"concrete_test":"Re-run the fixed-beta series beta3-DM0 through beta3-DM4 and beta6-DM0 through beta6-DM4 with the HShen and SFHo EOSs (the same EOSs used in Appendix B), using the corresponding GR1D Ye(rho) parametrizations for each EOS. Compute h2/h1 and h3/h1 for every model and overplot them on Fig. 11. If the MDM >= 0.03 Msun points remain outside the MDM=0 EOS+beta spread for all three EOSs, the threshold claim is confirmed; if any overlap occurs, the 0.03 Msun robustness statement must be weakened or the threshold increased.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The abstract's central claim ('Even within the uncertainties of nuclear matter equation of state, a DM core can be inferred if its mass is greater than 0.03 M_sun') rests on the light-shaded EOS uncertainty bands in Fig. 11, but those bands come from Appendix B, which explicitly tests three EOSs (HShen, LS220, SFHo) only for rotating AIC models without DM admixture. DM changes the collapse trajectory qualitatively and quantitatively: Table E1 shows rho_c,b decreases from 3.85e14 to 2.73e14 g/cm3 and Mic,b changes across the R5 series, and bounce is delayed from 33 ms to 93 ms. There is no reason to assume the EOS sensitivity of h2/h1 and h3/h1 is unaffected by a compact DM core, because the bounce occurs at different densities, entropies, and inner-core masses. Furthermore, the paper itself notes h2 changes by ~30% among EOSs and by ~30% under electron-capture rate variations (Sec. 3.3). If those ~30% shifts are comparable to the h2/h1 difference between MDM=0 and MDM=0.03 at fixed beta_ic,b, the threshold is not established. This is an evidentiary gap rather than a demonstrated failure; the claim may survive, but the current manuscript does not provide the necessary DM+EOS runs.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents axisymmetric, Newtonian hydrodynamical simulations of accretion-induced collapse (AIC) of rotating white dwarfs that contain a central, static, non-rotating dark matter core, modeled as 1 GeV fermionic DM. The authors vary the admixed DM mass MDM and the initial rotation rate, and extract the burst gravitational-wave signals using the quadrupole formula. They find that DM delays the plunge and bounce, lowers the central density and inner-core mass at bounce, and reduces the proto-neutron star mass. The gravitational waveforms exhibit three characteristic spikes around bounce; the amplitude ratios h2/h1 and h3/h1 depend on both MDM and the inner-core rotation parameter beta_ic,b, which is used to break the degeneracy between these two parameters. The paper also proposes rescaling relations for h1 and h2 to retrieve MDM and beta_ic,b, discusses the detectability with advanced LIGO, and estimates the effect of DM on PNS g-mode frequencies.","tokens_in":21776,"tokens_out":5781,"duration_ms":60469,"significance":"If the results hold, this is a useful and original contribution to the gravitational-wave literature on AIC: it adds DM admixture as a new physical degree of freedom, proposes a concrete observable (bounce-peak amplitude ratios) for disentangling DM mass from rotation, and connects the result to a possible indirect DM detection channel. The paper includes several checks that strengthen confidence in the core dynamics: convergence tests in Appendix D, comparison of three nuclear EOSs in Appendix B, central-density variation tests in Appendix C, and a waveform comparison with Abdikamalov et al. (2010) in Section 4.1. The peak-amplitude ratios are raw simulation outputs rather than fitted quantities, so the central observable is not circularly defined. The main significance, if the DM+EOS robustness can be established, is a potential new probe of DM in stellar collapse.","major_comments":[{"comment":"The central claim in the abstract that a DM core with MDM >= 0.03 Msun can be inferred 'even within the uncertainties of nuclear matter equation of state' is not directly supported by the simulations as presented. The light-shaded EOS uncertainty bands in Fig. 11 are taken from Appendix B, which compares HShen, LS220, and SFHo only for rotating AIC models without DM admixture. DM changes the bounce conditions substantially (Table E1: for the R5 series, rho_c,b decreases from 3.85e14 to 2.73e14 g/cm3 and tb increases from 33.2 to 93.1 ms as MDM goes from 0 to 0.04 Msun), so there is no demonstrated basis to assume that the EOS sensitivity of h2/h1 and h3/h1 is the same for DM-admixed collapses. The same caveat applies to the central-density robustness test in Appendix C, which appears to vary rho_c only for MDM=0. This is an evidentiary gap rather than a demonstrated failure, but the 0.03 Msun threshold requires either additional DM+EOS simulations or a more limited claim restricted to a fixed EOS.","section":"Section 2.2 / Section 3.3"},{"comment":"The two-fluid model assumes that the DM core is stationary and non-rotating, affecting the baryons only through a fixed external gravitational potential. Because the DM core's gravitational pull is the physical mechanism that changes the bounce density and the GW peak ratios used for MDM inference, this assumption is load-bearing. A dynamical response of the DM core to the collapsing baryons, or even a different assumed DM density profile, could shift rho_c,b, Mic,b, and hence h1, h2, and h3. I would like to see a concrete robustness test, such as a time-dependent DM potential or an explicit check with an altered DM profile, or at minimum a clear statement in the abstract and conclusions that the inference applies only within this modeling assumption.","section":"Section 3.3 / Eqs. (17)-(18)"},{"comment":"The retrieval procedure for MDM and beta_ic,b is calibrated and tested on the same set of simulations from which Eqs. (17) and (18) are fitted. This demonstrates internal consistency but not predictive power: an independent set of simulations, for example with different Ye parameterizations, rotation laws, or mass grids, is needed to establish that the h1* and h2* relations are universal enough for an actual GW observation. As written, the 'in principle' retrieval claim is somewhat stronger than the evidence provided.","section":"Section 3.3 / Richers et al. (2017)"},{"comment":"The paper cites roughly 30% variations in h2 (and in h1 - h2) with electron-capture rate variations, yet the >0.03 Msun threshold is framed only against EOS uncertainties. Since the DM indicator is built from h2/h1, the electron-capture uncertainty could also overlap the MDM=0 and MDM=0.03 signal separation. A quantitative statement of how much of the h2/h1 difference between MDM=0 and MDM>=0.03 remains after combined EOS and electron-capture uncertainties is needed to support the headline claim.","section":"Section 3.3 / Richers et al. (2017)"}],"minor_comments":[{"comment":"The fitted coefficient in Eq. (5) should specify the fit range and the number of models used, and the sentence containing 'W e found' in Section 2.1 contains a typographical spacing error.","section":"Section 2.1 / Eq. (5)"},{"comment":"The model labeling in Tables E1 and E2 should be cross-checked; in particular, the footnote 'aAlso Rmax-DM4' is unclear because R5-DM4 appears in the R5 block and Rmax-DM4 is not listed as a separate row.","section":"Table E1 / Table E2"},{"comment":"The caption of Fig. C1 states that GW amplitudes are multiplied by a constant to match h1, but the constant is not given; please state explicitly that this rescaling does not affect the amplitude ratios used in the analysis.","section":"Appendix C / Fig. C1"}],"recommendation":"major_revision","confidential_remarks":"The paper is within the journal's scope and the simulations appear carefully executed. My main reservation is that the headline EOS-robustness claim is not backed by DM+EOS simulations; this is fixable either by adding those runs or by softening the abstract and conclusions. The static-DM assumption is also worth an explicit robustness test. I would support publication after a major revision that addresses these points."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Bottom line: this is a solid extension of the group's DM-admixed white dwarf program, and it earns peer review, but the headline 0.03 Msun inference threshold is not yet demonstrated. The EOS uncertainty band in Fig. 11 comes from Appendix B, which varies the nuclear EOS only for MDM=0 models. DM changes the collapse trajectory enough (bounce delayed from 33 to 93 ms, central bounce density down by ~30%) that you cannot simply assume the peak-ratio sensitivity is EOS-independent. That is an evidentiary gap, not a demonstrated failure, but it should be fixed before the abstract claim is taken at face value.\n\nWhat is actually new and good: these are the first axisymmetric simulations of rotating DM-admixed WDs through AIC with GW extraction. The peak-ratio observables h2/h1 and h3/h1 are raw simulation outputs, not fits, so they are not circularly defined. The convergence tests (Appendix D), the initial central-density test (Appendix C), and the comparison with Abdikamalov et al. (2010) using their Ye profile all support the core dynamics. The rescaling relations for h1 and h2 are calibrated and tested on the same simulations, which is a weak test, but it is standard forward-model calibration, and the authors are appropriately cautious about microphysics uncertainties elsewhere.\n\nSoft spots, in proportion: the stationary, non-rotating DM core is a real simplification. If the collapsing baryons compress, heat, or mix with the DM, the bounce density and PNS mass shift, and the calibrated h1, h2 relations could move. The authors flag it as future work, which is honest, but it means the quantitative inference pipeline is model-dependent in a way not yet quantified. Also, the detection context is sobering: Galactic AIC rates around 1e-4 to 1e-3 per year put LIGO detection at the margin, and the paper says so itself. That does not undercut the physics, but it does temper the \"new indirect DM channel\" framing.\n\nThe citation pattern is fine: the paper builds on its own previous work, cites the relevant AIC and CCSN GW literature, and the comparison to Abdikamalov et al. is a genuine check, not a strawman.\n\nWho is this for: anyone working on DM effects in compact objects, AIC nucleosynthesis or transients, or GW bursts from stellar collapse. A serious referee should engage with it. The main request should be DM-plus-EOS runs for the threshold claim, or a clearly stated downgrade of that claim.\n\nRecommendation: send to peer review with an expectation of revision. The core simulation work is careful and reproducible enough to warrant referee time; the inference claim needs strengthening or softening.","headline":"Careful, genuinely new numerical study whose central EOS-robustness claim is not yet supported because the EOS tests were run only without dark matter.","tokens_in":22349,"tokens_out":1328,"would_cite":true,"duration_ms":17226,"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":"Gravitational-wave peak ratios can reveal a dark matter core inside a collapsing white dwarf.","keywords":["dark matter","gravitational waves","white dwarfs","accretion-induced collapse","proto-neutron star","core bounce","nuclear equation of state","numerical hydrodynamics"],"falsifier":"A detected accretion-induced-collapse burst whose peak-amplitude ratios match the dark-matter-free values $h_2/h_1=-2.14\\pm0.14$ and $h_3/h_1=1.37\\pm0.04$ while the absolute amplitudes imply a high $\\beta_\\mathrm{ic,b}$ would contradict the predicted DM shift; and an event yielding rescaled amplitudes ($h_1^*,h_2^*$) that do not fall on a single universal curve would falsify the parameter-retrieval scheme.","tokens_in":21266,"feed_emoji":"🌌","tokens_out":10542,"duration_ms":87463,"temperature":0.7,"pith_summary":"This paper argues that the burst of gravitational waves emitted when a rotating, dark-matter-admixed white dwarf collapses can reveal the presence and mass of the dark matter core. Around the moment of core bounce, the waveform has three characteristic amplitude spikes; the ratios between those spikes are nearly invariant for ordinary white dwarfs but shift systematically as the dark-matter core mass grows. Using axisymmetric hydrodynamical simulations, the paper shows that these ratios break the degeneracy between rotation rate and dark-matter mass, and that a core heavier than 0.03 solar masses can be identified even without pinning down the nuclear equation of state. This matters because stellar-collapse observations could then probe dark matter in a way direct-search experiments cannot.","feed_headline":"Gravitational waves can unmask dark matter cores in white dwarfs","feed_subtitle":"A dark-matter core over 0.03 solar masses shifts the bounce signal's peak ratios even under nuclear physics uncertainty.","key_machinery":"The central machinery is the set of amplitude ratios and rescalings derived from the bounce gravitational-wave spikes. The collapse is simulated with a two-dimensional (axisymmetric) Newtonian hydrodynamics code with an effective general-relativistic potential, a parameterized electron-capture profile, and a degenerate-electron-gas white dwarf; the dark matter is treated as a non-rotating, stationary fluid whose gravity is added to the Poisson equation. From the waveforms, the paper defines $h_1$, $h_2$, $h_3$ as the amplitudes of the three spikes around bounce and shows that $h_2/h_1$ and $h_3/h_1$ depend strongly on $M_\\mathrm{DM}$ while being only mildly sensitive to $\\beta_\\mathrm{ic,b}$; the rescaled quantities $h_1^* \\equiv h_1/[1-15.36(M_\\mathrm{DM}/M_\\odot)]$ and $h_2^* \\equiv h_2/[3.24-(1-11.6(M_\\mathrm{DM}/M_\\odot))^{-1}]$ then each collapse onto a universal curve against $\\beta_\\mathrm{ic,b}$. This two-step procedure — ratios for detection, rescalings for parameter extraction — is what carries the argument.","core_discovery":"The paper's central claim is that the collapse-bounce gravitational-wave signal of an accreting white dwarf carries a measurable imprint of an admixed dark-matter core. In the simulations, the quantities that govern the collapse—bounce time, central density at bounce, and proto-neutron-star mass—depend on the dark-matter mass $M_\\mathrm{DM}$ and the inner-core rotation parameter $\\beta_\\mathrm{ic,b}$ in a factorised way, and the amplitudes $h_1$, $h_2$, $h_3$ of the three dominant waveform spikes inherit this dependence. The ratios $h_2/h_1$ and $h_3/h_1$, which are nearly invariant for ordinary white dwarfs (about $-2.14$ and $1.37$), deviate systematically with $M_\\mathrm{DM}$, and the rescaled amplitudes $h_1^*$ and $h_2^*$ each follow a universal monotonic relation with $\\beta_\\mathrm{ic,b}$. The paper concludes that a dark-matter core mass $M_\\mathrm{DM} \\geq 0.03\\,M_\\odot$ can be inferred even within nuclear-equation-of-state uncertainties, and that smaller cores could be inferred if the equation of state is better constrained.","pith_inferences":["If the stationary-core assumption is relaxed, the dark-matter core itself would be compressed and heated during collapse, likely raising the bounce density and shrinking the amplitude ratios; two-fluid dynamical simulations would bound this systematic.","The universal rescaling curves suggest a direct parameter-estimation recipe for future detectors: measure $h_1$ and $h_2$, invert on $h_1^*(\\beta)$ and $h_2^*(\\beta)$, and the crossing point gives $(M_\\mathrm{DM}, \\beta_\\mathrm{ic,b})$; the recipe's practicality depends on calibrating the electron-capture profile, which neutrino-transport simulations can provide.","The same ratio-based fingerprint could be searched for in ordinary core-collapse supernova waveforms, where the 'extra central mass' would be the compact remnant itself rather than dark matter; the detection logic transfers.","The paper's claim that the 0.03 solar-mass threshold is insensitive to equation-of-state choice rests on only three nuclear equations of state; a wider ensemble, including exotic compositions, would map how far the threshold moves."],"forward_implications":["A gravitational-wave detection of an accretion-induced collapse could double as an indirect dark-matter detection from a stellar-scale object, something direct searches cannot probe.","If the inner-core rotation parameter can be pinned down from the absolute wave amplitude, the peak-amplitude ratios give the dark-matter core mass, breaking the rotation–dark-matter degeneracy.","The proto-neutron star left behind is lighter when the dark-matter core is heavier, shifting the frequency of post-bounce g-mode gravitational waves by up to roughly 20%; combining bounce and g-mode signals strengthens the inference.","For dark-matter cores above 0.03 solar masses the identification survives nuclear-equation-of-state uncertainty; below that, it hinges on a better-constrained equation of state."],"supporting_citations":[{"why":"Provides the rotating-AIC gravitational-wave waveforms and the beta_ic,b–amplitude correlations that this paper extends to dark-matter-admixed models.","marker":"Abdikamalov et al. 2010"},{"why":"Establishes the strong correlation between gravitational-wave amplitude and inner-core rotation parameter used to motivate the degeneracy analysis.","marker":"Dimmelmeier et al. 2008"},{"why":"Shows that h1 and h1-h2 increase monotonically with beta_ic,b for core-collapse supernovae, the basis for the ratio-based parameter extraction.","marker":"Richers et al. 2017"},{"why":"Supplies the two-fluid hydrostatic construction of dark-matter-admixed white dwarfs and the DM density profiles used as initial models.","marker":"Leung et al. 2015a"},{"why":"Gives the iterative method for constructing self-consistently rotating white-dwarf initial models.","marker":"Hachisu 1986"},{"why":"Provides the LS220 nuclear equation of state used for the standard simulations.","marker":"Lattimer & Swesty 1991"},{"why":"Provides the parameterized temperature profile and electron-capture treatment used to build AIC initial conditions.","marker":"Dessart et al. 2006"},{"why":"Supplies the GR1D code whose two-moment neutrino transport yields the electron-fraction profile parameterization.","marker":"O'Connor 2015"}],"fun_headline_variants":["Gravitational wave peak ratios unmask hidden dark matter cores","White dwarf collapse waves could expose dark matter cores","Dark matter cores in white dwarfs leave gravitational wave fingerprints","Peak ratios in AIC gravitational waves break dark matter degeneracy","Collapsing white dwarfs give away dark matter via wave peaks"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The dark matter core is assumed to stay fixed and rigid while the white dwarf collapses around it, pulling on the ordinary matter only through a constant gravitational field; if the collapse compresses, heats, or mixes with that core, the calibrated relations between dark-matter mass, rotation, and wave amplitudes no longer hold.","fun_headline_variants_meta":{"raw":{"variants":["Gravitational wave peak ratios unmask hidden dark matter cores","White dwarf collapse waves could expose dark matter cores","Dark matter cores in white dwarfs leave gravitational wave fingerprints","Peak ratios in AIC gravitational waves break dark matter degeneracy","Collapsing white dwarfs give away dark matter via wave peaks"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000396,"raw_usage":{"total_tokens":2147,"prompt_tokens":1091,"completion_tokens":1056,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":707,"completion_tokens_details":{"reasoning_tokens":973}},"tokens_in":707,"tokens_out":1056,"duration_ms":8386,"temperature":1.0,"reasoning_tokens":973,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T13:22:01.807095+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A detected accretion-induced-collapse burst whose peak-amplitude ratios match the dark-matter-free values $h_2/h_1=-2.14\\pm0.14$ and $h_3/h_1=1.37\\pm0.04$ while the absolute amplitudes imply a high $\\beta_\\mathrm{ic,b}$ would contradict the predicted DM shift; and an event yielding rescaled amplitudes ($h_1^*,h_2^*$) that do not fall on a single universal curve would falsify the parameter-retrieval scheme.","supporting_citations":[{"cited_title":"B., Ott, C","cited_arxiv_id":null,"evidence_quote":"Provides the rotating-AIC gravitational-wave waveforms and the beta_ic,b–amplitude correlations that this paper extends to dark-matter-admixed models."},{"cited_title":"D., Marek, A., & Janka, H.-T","cited_arxiv_id":null,"evidence_quote":"Establishes the strong correlation between gravitational-wave amplitude and inner-core rotation parameter used to motivate the degeneracy analysis."},{"cited_title":"1986, ApJS, 61, 479","cited_arxiv_id":null,"evidence_quote":"Gives the iterative method for constructing self-consistently rotating white-dwarf initial models."},{"cited_title":"M., & Swesty, F","cited_arxiv_id":null,"evidence_quote":"Provides the LS220 nuclear equation of state used for the standard simulations."},{"cited_title":"D., et al","cited_arxiv_id":null,"evidence_quote":"Provides the parameterized temperature profile and electron-capture treatment used to build AIC initial conditions."}],"review_version":1}