{"id":"a5f939ec-cb43-4211-a35a-cb2ea5b4be7c","arxiv_id":"2412.19262","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"Head-on collisions of fuzzy dark matter solitons leave an oscillating merged soliton that emits gravitational waves with periods from a few years to tens of years for particle masses of 1e-17 to 1e-18 eV/c^2.","lead":"This paper simulates head-on collisions of fuzzy dark matter solitons and calculates the gravitational waves emitted by the merged, oscillating soliton. A correct result would give a new way to connect low-frequency gravitational wave detections to the mass of the fuzzy dark matter particle.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The quoted GW periods are read from three undocumented PyUltraLight runs with no convergence, timestep, duration, or spectral analysis; the frequency claim is not yet established and is also λ-dependent.","rationale":"The reader's conditional verdict is appropriate and I do not see a reason to move away from it. The quadrupole computation itself is standard, the back-reaction ratio in Section IV is reassuring, and the C3 case is correctly identified as compromised by the unstable excited state. The load-bearing gap is that the central frequency numbers are numerical observables extracted from simulations that are not documented at the level needed for independent verification: no timestep, no total duration, no convergence test, and no spectral analysis. Energy conservation alone cannot rule out underresolution or periodic-box contamination. The additional lambda dependence, via Eq. (5) and the chosen lambda in Table II, means the abstract's mass-only frequency statement is also over-generalized unless the corresponding halo mass is specified. A focused resolution/box-size and spectral-peak test would settle whether the quoted few-year and few-ten-year periods are robust or are artifacts of the specific numerical setup.","tokens_in":10056,"tokens_out":16547,"duration_ms":170847,"concrete_test":"Rerun C1 (and C2) with (N, l) = (256, 2e4), (512, 2e4), (256, 4e4), and (512, 4e4), recording timestep and total duration; compute the power spectral density of h_+(t) from Eq. (18) over at least 10 cycles and require the dominant peak to be stable to better than 10% across resolution and box size. Also rerun with lambda=5.78e-8 from Eq. (8) for M_halo=1e12 M_sun and compare the resulting peak frequency to the lambda-scaled expectation. If the peak shifts with resolution/box size, or if the canonical halo gives a period about 30 times longer, the abstract frequency claim must be qualified.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central prediction in the abstract and Section IV is a frequency read off three PyUltraLight runs (Section III, Table II) with N=256, l=2e4, and unreported timestep and total duration. The post-collision signal is at dimensionless frequencies around 10^-6 to 10^-8, i.e. periods of 10^6 to 10^8 code-time units; no convergence study, spectral estimate, or statement of how many cycles were captured is given. Energy conservation (Figs. 2-4) does not establish that the low-frequency quadrupole oscillation is physical, since a pseudo-spectral solver can conserve discrete energy while underresolving field gradients, and periodic images at distance 2e4 (with soliton radius roughly 2e3 for lambda=1.66e-6) can contaminate the quadrupole moment in Eq. (16). In addition, by the scaling symmetry (5), the dimensionless oscillation frequency scales with the chosen lambda, so the mass-only statement in the abstract is only for the specific lambda=1.66e-6, which via Eq. (8) corresponds to M_halo ~ 1.5e14 M_sun; for a typical 1e12 M_sun halo, lambda=5.78e-8 and the periods would be about 30 times longer.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies gravitational-wave emission from the oscillating merged object that forms after head-on collisions of fuzzy dark matter solitons. It introduces a dimensionless unit system, computes ground and first-excited soliton profiles, simulates three collision setups (C1, C2, C3) with PyUltraLight at N=256 in a box of length l=2e4, and applies the quadrupole formula to the post-collision density field to obtain h_+ and h_× waveforms. The central quantitative claim is that post-collision GWs have periods of 'few ten years' for FDM mass m=10^-18 eV/c^2 and 'few years' for m=10^-17 eV/c^2, and that GW back-reaction is negligible because Eg/ΔEp is of order 10^-12.","tokens_in":10259,"tokens_out":14114,"duration_ms":145913,"significance":"If established, the mass-dependent frequency prediction would be a falsifiable, novel signature of FDM soliton mergers and would give a concrete target for GW searches. The paper uses the standard Schrödinger-Poisson system, a clean dimensionless rescaling, and an explicit back-reaction check (Eq. 19 and the quoted ratios), and I see no circularity in the construction of the signal: the frequency-mass relation is a scaling consequence of the simulation, not fitted to a target. The claim is not yet established, however, because the headline frequencies are read from three runs with no reported timestep, duration, spectral analysis, or convergence study, and because the quoted periods depend on the halo mass through λ even though the abstract presents them as functions of m alone. The manuscript has the structure of a proof-of-principle study; with additional numerical documentation and a corrected parameter statement, the central idea is worth publishing.","major_comments":[{"comment":"The paper never states the timestep Δt, the total simulation time, or how many post-merger cycles are used to read the GW frequency, and it provides no resolution or box-size convergence study. The periods quoted in the abstract correspond to dimensionless frequencies of roughly 2–7×10^-7 for m=10^-18 eV/c^2 and 10^-17 eV/c^2, yet the only frequency quoted in the text is '~10^-8' (Section III and Section V); if that quoted value were the GW frequency, the periods would be of order 10^3 yr, not 'few ten years.' No Fourier spectrum of h_+ is given, so the frequency underlying the abstract cannot be checked from the paper. This matters because the box has l=2×10^4 while the C1/C2 solitons have λ=1.66×10^-6 and radii of order 2×10^3 code units; with periodic boundary conditions, image contamination of the quadrupole moment in Eq. (16) is not excluded. Energy conservation in Figs. 2-4 does not discipline this low-frequency, long-wavelength mode, since a pseudo-spectral solver can conserve discrete energy while underresolving phase gradients. I request the timestep and total duration, a convergence study (e.g., N=128/256/512 or l=4×10^4), and a spectrum of h_+.","section":"Section III, Table II; Section IV, Figs. 5-7"},{"comment":"The headline frequency is presented as a function of the FDM particle mass m alone, but in these simulations the dimensionless frequency is set by the soliton scale λ. Runs C1 and C2 use λ=1.66×10^-6, which via Eq. (8) corresponds to M_halo≈1.5×10^14 M_sun, not a typical galaxy-scale halo. For M_halo=10^12 M_sun, Eq. (8) gives λ=5.78×10^-8, and the scaling symmetry of the SP system (which requires t̃→λ^-1 t̃ in addition to the listed transformations) shifts all dimensionless frequencies down by a factor of roughly 29 relative to the quoted runs. The abstract should therefore either quote the halo mass or λ at which the period statements hold, or state the frequency as an explicit function f(m, M_halo). As written, the mass-only statement in the abstract is incomplete and cannot be reproduced from the simulation setup described in Table II.","section":"Abstract, Section IV, Eq. (8), Table II"},{"comment":"Run C3 is described by the authors themselves as not reaching a collision: the soliton with the first excited-state profile is unstable and is 'dismembered before collision.' Consequently, the waveform in Fig. 7 is not a post-collision signal from two merging solitons and should not be counted as support for the paper's central claim. This run should either be removed from Section IV or explicitly labeled as a failed/stability test, with its waveform not presented as a post-collision GW result. The remaining evidence for the headline frequency therefore rests on C1 and C2 alone, which strengthens the need for the numerical documentation requested above.","section":"Section III, Table II (C3); Section IV, Fig. 7"}],"minor_comments":[{"comment":"There are repeated typos: 'sollitons' and 'FDM soltions' should be 'solitons', and 'ring-dwon' in Section V should be 'ring-down'.","section":"Throughout, especially Section IV and figure captions"},{"comment":"The paper asserts that the GW frequency is 'mainly determined by the FDM soliton size and mass but not sensitive to the angular momentum' without a calculation or test. A brief argument or a run with non-zero impact parameter would make this claim credible.","section":"Section V"},{"comment":"All simulations set the initial phases δ_1=δ_2=0 and the same initial speed |v|=10^-4. A sentence on the expected sensitivity of the post-collision frequency to these choices would help the reader gauge the generality of the quoted periods.","section":"Table II"},{"comment":"The scaling symmetry is written for the static reduced system (4) and omits the transformation of time; stating the full symmetry, including t̃→λ^-1 t̃, would remove ambiguity in the frequency-scaling argument and in the derivation of Eq. (8) from the soliton-halo relation.","section":"Eq. (5)"}],"recommendation":"major_revision","confidential_remarks":"The manuscript fits the journal's scope and the central idea is not circular, but the numerical evidence for the headline frequency is currently incomplete. I would ask the authors to supply the missing simulation metadata (timestep, duration, number of cycles), a convergence test, and a spectral estimate of h_+, and to revise the abstract so that the halo-mass dependence through λ is explicit. With those changes the paper could be acceptable; without them, the abstract's mass-only frequency claim is not supported by the reported runs."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThe paper gives a concrete GW frequency prediction from fuzzy dark matter soliton post-merger oscillations, and that specific result is genuinely new relative to the cited literature. The method is routine — a PyUltraLight simulation plus the quadrupole formula — but applying it to this post-collision oscillation is the new step. The back-reaction check (Eg/DeltaEp ~ 1e-12) is a sensible internal consistency test, and the paper is explicit about the linearized-theory approximation. The frequency claim is, in principle, falsifiable by pulsar timing arrays.\n\nThe soft spots are real, though. The central frequency claim is read from three simulations with N=256, box length 2e4, no stated timestep, total duration, or resolution convergence study. The post-collision GW frequencies are around 1e-6 to 1e-8 in code units, meaning periods of 1e6 to 1e8 code-time units; without a statement of how many cycles were captured, or a spectral estimate, the numbers could be contaminated by underresolution or periodic-box artifacts. Energy conservation alone does not establish the physicality of the low-frequency oscillation, since a pseudo-spectral solver can conserve energy while misrepresenting spectral content.\n\nThe C3 run uses a first-excited-state soliton that the paper itself says dismembers before collision; that run should not count as a clean post-collision signal. And the abstract's mass-frequency statement is incomplete: the quoted periods correspond to lambda = 1.66e-6, which through Eq. (8) means a halo of about 1.5e14 solar masses, not a typical galaxy halo. For a 1e12 solar mass halo the periods would be tens of times longer. That qualification belongs in the abstract, not as a detail in the setup table.\n\nThe paper also over-sells the new adimensional unit system; it is the standard reduced Compton length/time/mass, and the advantage over earlier units is bookkeeping.\n\nWhat is here is a plausible mechanism and a clean calculation once you accept the simulation outputs. But the manuscript does not yet establish the quantitative prediction. The authors need to add a resolution study, report timestep and run duration, give a spectral estimate of the oscillation, and present the halo-mass dependence explicitly. With those additions it could be a useful contribution to the FDM-GW literature.\n\nFor a reading group: maybe, if the group wants to discuss how numerical papers should document convergence. I would not cite it yet in my own work.\n\nThe paper deserves a serious referee. I would send it to review with a request for major revisions — not because the idea is wrong, but because the numerical support for the central claim is missing.","headline":"Plausible new FDM soliton GW source, but the frequency claim rests on underdocumented simulations and a halo-mass choice the abstract omits.","tokens_in":10864,"tokens_out":5464,"would_cite":false,"duration_ms":49082,"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":"This paper shows that the merged remnant of two colliding fuzzy dark matter solitons oscillates and emits gravitational waves whose period is set by the dark matter particle mass—tens of years for $10^{-18}\\,{\\rm eV}/c^2$ and a few years…","keywords":["fuzzy dark matter","soliton collision","gravitational waves","Schrodinger-Poisson","quadrupole formula","ultralight scalar field","pulsar timing arrays","dark matter cores"],"falsifier":"Re-run the three collision setups at resolutions $N=128$, $256$, and $512$ and with a larger box at fixed resolution, then compare the frequency and amplitude of the post-merger $h_+$ waveform; if the dominant frequency shifts by more than a few percent or depends on the box size, the claimed period prediction is not settled. A direct pulsar-timing-array detection of a periodic signal with the predicted period from a nearby soliton merger would confirm the claim.","tokens_in":9779,"feed_emoji":"🌊","tokens_out":7457,"duration_ms":65670,"temperature":0.7,"pith_summary":"This paper claims that the remnant of two colliding fuzzy dark matter solitons—the lumpy cores that ultralight dark matter forms—oscillates and emits gravitational waves with a period set by the dark matter particle mass. Simulating head-on collisions with the Schrödinger–Poisson equations and applying the quadrupole formula, the authors find wave periods of tens of years for a particle mass of $10^{-18}\\,\\mathrm{eV}/c^2$ and a few years for $10^{-17}\\,\\mathrm{eV}/c^2$. If correct, the predicted frequency gives a concrete, mass-dependent target for future low-frequency gravitational-wave searches, and a detection would constrain the fuzzy dark matter particle mass and the properties of solitons. The paper also introduces a new set of dimensionless units that make the simulations independent of particle mass, so results can be rescaled to any FDM mass without rerunning.","feed_headline":"Dark matter soliton collisions may ring every few to tens of years","feed_subtitle":"Ultralight dark matter's mass fixes the wave period, giving pulsar timing arrays a specific target to chase.","key_machinery":"The key machinery is the Schrödinger–Poisson system solved with the pseudo-spectral code PyUltraLight under periodic boundary conditions, combined with the quadrupole formula for gravitational waves (Eqs. 15–18). The SP system evolves the FDM wavefunction and its self-gravitational potential; the three head-on collision setups—equal-mass ground-state, unequal-mass ground-state, and ground-state versus first-excited-state—produce a merged soliton whose energy exchange between kinetic and potential forms gives the oscillating quadrupole moment. A new scale system (Eq. 6) expresses length, time, and mass in terms of $\\hbar/(mc)$, $\\hbar/(mc^2)$, and $\\hbar c/(Gm)$, making the dimensionless dynamics independent of the FDM mass. The gravitational-wave strain is then computed from the second time derivative of the quadrupole moment of the density field.","core_discovery":"The central claim is that gravitational waves from the post-collision stage of FDM soliton mergers are not negligible and their frequency is determined by the soliton size and the FDM mass. With the linearized quadrupole formula applied to the simulated density evolution, the authors predict GWs with a period of tens of years for $m=10^{-18}\\,\\mathrm{eV}/c^2$ and a few years for $m=10^{-17}\\,\\mathrm{eV}/c^2$. The waves originate from the irregular spherically asymmetric oscillation of the merged soliton after a head-on collision, with the $h_\\times$ polarization vanishing by symmetry. The paper further asserts that gravitational-wave back reaction is negligible, since the energy radiated is only about $10^{-12}$ of the kinetic–potential energy exchange in the simulated period. The new dimensionless units proposed here decouple the simulation from the particle mass, so all runs are mass-independent and only the physical rescaling changes with $m$.","pith_inferences":["If the paper's picture is right, pulsar timing arrays could search for ultra-low-frequency signals from individual nearby soliton mergers, stacking many events to compensate for the small strain.","A natural extension is to simulate off-axis collisions; breaking the head-on symmetry would generate $h_\\times$ and likely change the amplitude, while the frequency may remain set by the soliton size as the paper argues.","The predicted periods depend on the post-merger oscillation being physical rather than a grid artifact; a resolution-convergence study would settle whether the dimensionless frequency near $10^{-5}$ is stable as the grid is refined.","If soliton mergers are common across cosmic history, their incoherent sum could form a stochastic gravitational-wave background in the nHz band, giving pulsar timing arrays an independent way to constrain the FDM model."],"forward_implications":["The gravitational-wave period scales directly with the FDM particle mass: $10^{-18}\\,\\mathrm{eV}/c^2$ gives periods of tens of years and $10^{-17}\\,\\mathrm{eV}/c^2$ gives a few years, so a detection fixes a narrow mass window.","Because the dimensionless simulation is mass-independent, the same waveforms can be rescaled to any FDM particle mass without additional runs, covering a continuous parameter space from one set of collision simulations.","The vanishing $h_\\times$ polarization in head-on collisions gives a symmetry check: any observed $h_\\times$ would indicate the collision was not head-on or that angular momentum plays a role.","The tiny radiated energy compared with the internal energy transfer justifies ignoring gravitational-wave back reaction for these systems, but this approximation would need rechecking for more massive or faster collisions.","Future detection of such waves would constrain the FDM particle mass and the soliton–halo mass relation used to set the initial soliton masses."],"supporting_citations":[{"why":"Defines fuzzy dark matter as an ultralight scalar field with kpc-scale de Broglie wavelength, setting the physical context for solitons.","marker":"[23]"},{"why":"Establishes soliton solutions and evolution of the Schrödinger–Newton/SP system, providing the equilibrium configurations used as initial conditions.","marker":"[24]"},{"why":"Describes the PyUltraLight pseudo-spectral solver used for all collision simulations in this paper.","marker":"[27]"},{"why":"Supports the identification of the SP system as the weak-field limit of the Einstein–Klein–Gordon system, justifying the use of linearized theory for GW emission.","marker":"[29]"},{"why":"Supplies the soliton–halo mass relation used to set the rescaling factors $\\lambda_n$ for the colliding solitons.","marker":"[32]"},{"why":"Provides the analogy of gravitational waves from post-collision oscillations of Proca stars, motivating the quadrupole formula application here.","marker":"[33]"}],"fun_headline_variants":["Fuzzy soliton crashes may ring gravitational waves every few years","Dark matter soliton mergers could produce decade-period waves","Soliton collisions in fuzzy dark matter might create GW pulses","Post-collision dark matter solitons could emit low-frequency waves","Fuzzy dark matter smashups may signal via pulsar timing"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that three computer runs, each with 256 grid points per side in a box of side length 20,000 simulation units, truly capture the sloshing of the merged dark matter clump that sets the wave period; the paper reports no test that a finer grid or a larger box would give the same period.","fun_headline_variants_meta":{"raw":{"variants":["Fuzzy soliton crashes may ring gravitational waves every few years","Dark matter soliton mergers could produce decade-period waves","Soliton collisions in fuzzy dark matter might create GW pulses","Post-collision dark matter solitons could emit low-frequency waves","Fuzzy dark matter smashups may signal via pulsar timing"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000208,"raw_usage":{"total_tokens":1404,"prompt_tokens":945,"completion_tokens":459,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":561,"completion_tokens_details":{"reasoning_tokens":374}},"tokens_in":561,"tokens_out":459,"duration_ms":4845,"temperature":1.0,"reasoning_tokens":374,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T00:48:25.000143+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Re-run the three collision setups at resolutions $N=128$, $256$, and $512$ and with a larger box at fixed resolution, then compare the frequency and amplitude of the post-merger $h_+$ waveform; if the dominant frequency shifts by more than a few percent or depends on the box size, the claimed period prediction is not settled. A direct pulsar-timing-array detection of a periodic signal with the predicted period from a nearby soliton merger would confirm the claim.","supporting_citations":[{"cited_title":"Evolution of the Schrodinger-Newton system for a selfgravitating 9 scalar field,","cited_arxiv_id":null,"evidence_quote":"Establishes soliton solutions and evolution of the Schrödinger–Newton/SP system, providing the equilibrium configurations used as initial conditions."},{"cited_title":"Klein-Gordon Geon,","cited_arxiv_id":null,"evidence_quote":"Supports the identification of the SP system as the weak-field limit of the Einstein–Klein–Gordon system, justifying the use of linearized theory for GW emission."}],"review_version":1}