{"id":"7f8d9f3c-e6e8-482c-9844-cbb9bf6d1a59","arxiv_id":"2509.03272","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"Dark matter capture onto neutron stars is amplified by about 4-5 times in circular binaries, tightening cross-section limits from GW170817 and capping the accreted DM fraction near 1e-3.","lead":"Dark matter falling onto neutron stars in binary pairs gets funneled more efficiently by the combined gravity of the two stars, boosting the capture rate by a factor of 4 to 5. The authors use this boost to tighten limits on dark matter interactions with ordinary matter from the neutron star merger GW170817.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The MC flux is measured at r<0.1 r_H and the claimed constancy of Rq down to the NS radius is asserted without shown data; this unsupported extrapolation carries the 4-5 factor into all physical constraints.","rationale":"I read the central claim as a numerical statement about Rq, with the GW170817 and fchi results following from it. The reader's circular-orbit concern is legitimate, but it attacks the application rather than the raw simulation result; the paper explicitly restricts its headline to circular binaries. My concern targets the unshown radius extrapolation, which is the only bridge between the simulated encounter rate at 0.1 r_H and the NS-surface flux used in every derived bound. I also checked the apparent tension between Table I and the text '~4 at Hubble time for vchi = 220 km/s': with a Maxwellian-weighted average, the low-velocity tail dominates the flux integral, so <Rq> ~ 4 is compatible with the fitted R(x); that is not a problem. The paper is transparent about its main caveats (eccentricity, gas accretion, late formation), but the constancy of Rq in r is the load-bearing assertion I could not verify from the text. I would keep the reader's CONDITIONAL verdict: the central result is plausible but should not be used for constraints until the r-dependence check or the numerical data are provided.","tokens_in":21352,"tokens_out":24805,"duration_ms":245217,"concrete_test":"Re-run the Sec. II A Monte Carlo with the same sampling and importance splitting, but record cumulative closest-approach fluxes at r/a = 1e-3, 1e-4, 1e-5, and 1e-6 in addition to 0.1 r_H/a, for q = 1 and q = 0.28 and for v/Vb values of 0.3, 0.6, and 1.0. If Rq changes by more than 10% over this range, the constant-Rq assumption fails and the GW170817 cross-section limits and accreted-fraction bounds in Figs. 4-5 need recomputation; if it is flat, the central numerical factor is supported.","verdict_should_be":"UNCHANGED","load_bearing_attack":"In Sec. II A the authors measure close encounters at r<0.1 r_H and then state that Rq is constant for multiple orders of magnitude in r, but no figure, table, or error estimate showing this convergence is provided. For an equal-mass NS binary, r_H ≈ 0.55 a, so 0.1 r_H ≈ 0.055 a, while the physical capture surface is r_NS/a ≈ 1e-6; the application therefore extrapolates over roughly five orders of magnitude in radius. The scale-free argument fixes only the dependence on v and Vb for fixed r/a; it does not remove a possible r/a dependence. If the close-encounter distribution has a bottleneck at intermediate radii, Rq at r_NS could differ from Rq at 0.1 r_H. Because Eq. (17), Eq. (20), and Figs. 4-5 all scale linearly with Rq, even a 20% error in this extrapolation changes the derived DM bounds and the accreted-fraction cap. The absence of released code or data makes this check currently irreproducible. This is an internal gap in the numerical support for the central factor, distinct from the astrophysical circular-orbit caveat that the authors do flag.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper argues that the flux of dark matter particles onto a stellar-mass compact object is amplified when the object is in a circular binary, because test-particle orbits in the restricted three-body problem are chaotic and can undergo repeated close encounters. The authors measure this amplification R_q(v/V_b) with a Monte Carlo simulation for two mass ratios (q=1 and q=0.28), fit the result to an empirical curve saturating at R_0 ~ 3.5-4.0, and then apply the factor to the DM capture rate of neutron stars in inspiraling binaries. Using this amplified capture rate, they derive updated upper bounds on the DM-nucleon scattering cross section from the non-collapse of the GW170817 binary, and they derive a cap on the accreted DM mass fraction enforced by dynamical friction from ambient DM. The main quantitative claim is that the amplification factor is approximately 4-5 for circular binaries with q ≳ 0.3 and small DM velocity dispersion, and that this tightens previous bosonic DM constraints by the same factor.","tokens_in":21610,"tokens_out":14260,"duration_ms":138430,"significance":"If the central amplification factor is correct, the paper adds a genuinely new physical ingredient to DM capture calculations: binary gravitational scattering can enhance the effective accretion flux beyond the isolated-star value. The simulation setup follows the established method of Quinlan (1993) and is described in enough detail that the procedure is credible. The amplification factor is computed independently of the DM-nucleon cross section, so there is no circularity in the resulting constraints. The GW170817-based bounds and the dynamical-friction mass-fraction cap are concrete, falsifiable predictions. The main risks are the unsupported small-radius extrapolation of R_q and the absence of statistical or systematic error estimates on the fitted amplification curve; both directly affect the claimed constraints. No code or data release is mentioned, which limits reproducibility of the central numerical result.","major_comments":[{"comment":"The constancy of R_q as a function of r/a is asserted but not demonstrated. The simulation records close encounters at r < 0.1 r_H, and for an equal-mass NS binary r_H = a (m_i/3M)^{1/3} ~ 0.55 a, so 0.1 r_H ~ 0.055 a, whereas the NS capture surface is at r_NS/a ~ 10^{-6}. The scale-free argument fixes only the dependence on v and V_b at fixed r/a; it does not exclude an additional r/a dependence at intermediate radii. Since Eqs. (17) and (20) and Figs. 4-5 all scale linearly with R_q, an unquantified 20% change in the extrapolation would directly shift the derived cross-section bounds and the DM-fraction cap. Please provide a figure or table showing R_q as a function of r/a for several v/V_b values down to the smallest resolved radius, together with statistical uncertainties, and state the range over which the plateau is actually verified.","section":"Sec. II A, Eq. (7)"},{"comment":"The empirical fit parameters R_0 and w are quoted without uncertainties, and the data points in Fig. 2 carry no error bars. Because the amplification factor enters as a linear prefactor in Eq. (17), the Monte Carlo statistical error in R_0 and the systematic uncertainty of the fitting form in Eq. (10) should be propagated into the cross-section bounds in Fig. 4. Additionally, the abstract's claim for q ≳ 0.3 rests on only two simulated mass ratios, q = 1 and q = 0.28; if that claim is retained, the authors should either simulate additional mass ratios or provide an explicit interpolation uncertainty for the q-dependence.","section":"Fig. 2, Table I"},{"comment":"The final Discussion correctly acknowledges that the calculations assume negligible eccentricity and that eccentricity can significantly change the GW inspiral time. However, the GW170817 cross-section constraint and the DM-fraction limit in Figs. 4-5 are derived under the circular-orbit, GW-only inspiral assumption. Since the title highlights the GW170817 constraint, the authors should quantify how the integrated amplification changes for moderate initial eccentricities, for example by integrating Eq. (20) using Peters' eccentric-orbit inspiral rate, or explicitly restrict the title claim to binaries that circularize early.","section":"Discussion (final paragraph), Sec. IV"}],"minor_comments":[{"comment":"In the sentence 'it is a topic of strong debate wether NS binaries', 'wether' should be 'whether'.","section":"Sec. IV A"},{"comment":"The figure would benefit from explicit axis labels for x = v/V_b and y = R_q, and a legend identifying the error bars or noting their absence.","section":"Fig. 2"},{"comment":"The table should report the uncertainties in R_0 and w obtained from the fit, and the text should state the reduced chi-squared or equivalent goodness-of-fit measure.","section":"Table I"},{"comment":"The sentence 'All the orbits are unstable and thus eventually become expelled except for a set of initial condition with a volume of zero in the initial parameter space' is too strong given that integrations are truncated at 10^6 steps; please clarify how the truncation affects the recorded encounter statistics.","section":"Sec. II A"},{"comment":"The abstract refers to mass ratios q ≳ 0.3, but Section IV restricts the constraints to binary NSs with q = 1; aligning the wording would avoid overstating the coverage.","section":"Sec. IV, abstract"}],"recommendation":"major_revision","confidential_remarks":"I recommend major revision rather than rejection because the central mechanism is plausible and the missing r/a-plateau check is likely straightforward to supply from the existing simulation. I would not accept the paper without a quantitative demonstration that R_q is constant down to the NS radius, with error bars, and without propagating the fit uncertainties through the constraints. The circular-orbit caveat is acknowledged by the authors, but it should be elevated in the presentation because the GW170817 constraint is the headline result."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The genuinely new result is the ~4-5 amplification factor for DM capture in circular compact binaries and the resulting boost to bosonic DM bounds from NS mergers. The paper is the first to quantify this for DM on neutron stars and to propagate it into cross-section constraints, and it does that cleanly. The MC setup is standard restricted three-body scattering in the spirit of Quinlan, described in enough detail to follow, and the empirical fit is simple and sensible. The timescale checks on ZKL, flybys, and Bondi accretion are useful, and the authors are honest about the circular-orbit limitation in the Discussion. That part is not hidden.\n\nThe soft spots are real but not necessarily fatal. The biggest is the r/a extrapolation. Encounters are recorded at r < 0.1 r_H, which for equal-mass NSs is about 0.05a, while the NS surface is roughly 10^-6 a. That is five orders of magnitude, and there is no plot, table, or error estimate showing that R_q has actually converged over that range. The scale-free argument fixes the dependence on v/V_b but not a possible residual dependence on r/a. Since Eq. (20) and all the derived constraints scale linearly with R_q, even a 20% error in that extrapolation shifts the bounds by 20%. The stress-test note is on target: the authors assert constancy without showing it, and the absence of code or data makes the check impossible right now.\n\nSmaller issues: Fig. 2 has no error bars, the fit parameters in Table I have no uncertainties, only two mass ratios are simulated, and the interpolation to q > 0.3 is an assumption. None of these is disqualifying, but they should be quantified. The circular-orbit caveat is a real astrophysical limitation, especially for late-formed binaries, but the authors flag it explicitly. The single-scattering capture model is standard for the mass range considered, and the GW170817 density choice is conservative. The citation pattern is fine; the relevant stellar-dynamics and TDE literature is cited, and the self-citations are not load-bearing.\n\nWho is this for? People working on DM capture in compact objects and on using NS mergers as DM probes. It deserves a serious referee, not a desk reject. I would send it out with a specific request: show the convergence of R_q over radius with error bars, or release the simulation code and data, and either add a third mass ratio or state clearly how the interpolation error affects the constraints.","headline":"Plausible new amplification factor for DM capture in binaries, worth refereeing, but the central factor rests on an unsupported five-order-of-magnitude extrapolation that needs to be shown.","tokens_in":700,"tokens_out":862,"would_cite":true,"duration_ms":37768,"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":"The paper shows that neutron stars in binaries capture dark matter four to five times faster than isolated stars, tightening bosonic dark-matter bounds from GW170817.","keywords":["dark matter capture","neutron star binaries","GW170817","bosonic dark matter","dynamical friction","restricted three-body problem","gravitational wave constraints","amplification factor"],"falsifier":"Measure the eccentricity distribution of merging neutron-star binaries from gravitational-wave waveforms: if a substantial fraction merge with eccentricity larger than about 0.1 at frequencies where the orbital velocity is comparable to the dark-matter velocity dispersion, the circular-orbit integration fails as a description of the accretion history. Equivalently, a single old neutron-star merger in a high-dark-matter-density host whose delay time and dark-matter profile are well measured would let the bound be recomputed from event-specific parameters rather than the assumed 0.1 GeV/cm$^3$ and 220 km/s values, checking whether the factor-of-four tightening survives.","tokens_in":21137,"feed_emoji":"🌌","tokens_out":6206,"duration_ms":57490,"temperature":0.7,"pith_summary":"The paper argues that dark matter accretion onto a compact object is not an isolated-star process: in a binary, dark-matter particles on chaotic three-body orbits repeatedly pass close to the neutron stars, so the effective capture flux is amplified by a factor of about 4–5 for circular binaries, mass ratios above roughly 0.3, and dark-matter velocity dispersions below the binary orbital speed. Using this amplification, the paper tightens the upper bound on the scattering cross section of non-annihilating bosonic dark matter with nucleons, using the observed binary neutron star merger GW170817 as the anchor. It further shows that dynamical friction from ambient dark matter limits how much dark matter neutron stars in binaries can accumulate before merger, capping the accreted mass fraction at roughly one part in a thousand. A sympathetic reader would care because most neutron stars that merge spent their lives in binaries, so the isolated-star rates used in earlier constraints may have underestimated the dark-matter exposure by about a factor of four.","feed_headline":"Neutron-star binaries capture 4-5 times more dark matter","feed_subtitle":"The boost tightens dark-matter limits from GW170817 and caps accreted mass at 0.1 percent.","key_machinery":"The load-bearing object is the restricted three-body scattering problem: test particles with a fixed velocity at infinity are injected toward a circular binary of masses $m_1,m_2$, and their orbits are integrated numerically, recording close encounters within the Hill sphere of each component. Because the test-particle orbits are non-integrable, a single particle can have multiple close encounters with the same star or with both stars, and the resulting flux ratio $R_q=\\Sigma_b/\\Sigma_s$ is the amplification factor. Gravitational dynamics is scale-free, so $R_q$ depends only on the dimensionless ratio $v/V_b$ of dark-matter velocity to binary orbital velocity, not on the individual values. The amplification is then folded into the standard capture-rate formula for neutron stars and integrated along the standard gravitational-wave inspiral for circular binaries, while the same scattering machinery supplies the dynamical-friction hardening rate that caps the total accreted mass.","core_discovery":"On the paper's own terms, the central discovery is that the accretion flux of dark matter on neutron stars in a circular binary is not the sum of two isolated fluxes but a larger flux described by an amplification factor $R_q(v/V_b)$, which saturates near 4.5 when the dark-matter velocity is small compared with the binary orbital velocity and drops to unity when $v\\gg V_b$. The factor is computed with few-body Monte Carlo scattering, fitted empirically as $R(x)=1+R_0[1+(x/w)^6]^{-1/2}$, and shown to be nearly constant over many orders of magnitude in encounter radius for $r\\ll a$. Integrated over a gravitational-wave-driven inspiral, this gives about a fourfold increase in the number of captured particles for binaries that merge in roughly a Hubble time with a standard velocity dispersion of 220 km/s. From that increase, the paper derives a tighter upper bound on the bosonic dark-matter–nucleon scattering cross section from GW170817, and derives an upper bound of $\\lesssim10^{-3}$ on the dark-matter mass fraction inside neutron stars at merger, enforced by the dynamical-friction drag of the ambient dark matter on the binary orbit.","pith_inferences":["The authors apply the amplification to neutron stars, but the mechanism is purely gravitational and scale-free, so the same factor should multiply dark-matter capture in other compact or stellar binaries, such as white-dwarf binaries and black-hole–neutron-star binaries, whenever the binary spends a Hubble time at separations where $v_\\chi\\lesssim V_b$.","They leave implicit that the amplification also raises the effective dark-matter exposure of old binaries in dense environments; a statistical sample of mergers with measured host-galaxy delay times could convert the single-event GW170817 bound into a population-level constraint.","A testable extension is to check eccentric binaries: if a population of neutron-star mergers is found with non-negligible eccentricity, the circular-orbit integration used here would need to be repeated, and the resulting integrated amplification could be either larger or smaller depending on how much time the binary spends at wide separations.","Because the amplification is geometric, it applies to any weakly interacting particle captured by neutron stars, not just bosonic ones; the authors note this for the mass-fraction cap, but the cross-section improvement would also propagate to other capture-based probes such as kinetic heating."],"forward_implications":["Dark-matter capture onto neutron stars in merging binaries is roughly four to five times the isolated-star rate over the inspiral, so earlier constraints built on isolated stars are conservative by about that factor.","The upper bound on the bosonic dark-matter–nucleon scattering cross section derived from the observed binary neutron star merger GW170817 is correspondingly tighter than the isolated-star bound for the same host-galaxy assumptions.","The accreted dark-matter mass fraction in binary neutron stars at merger is at most about $10^{-3}$, even in extreme dark-matter spikes, because dynamical friction drives the binary to merge before much mass can accumulate.","The dynamical-friction cap applies to fermionic as well as bosonic dark matter, since it does not rely on collapse or annihilation inside the star.","For the highest ambient dark-matter densities, the dynamical-friction dephasing of the gravitational-wave signal could become observable with upcoming space-based detectors."],"supporting_citations":[{"why":"Establishes the gravitational-capture formalism for dark matter by stars that the binary version extends.","marker":"[1]"},{"why":"Gives the isolated-neutron-star bosonic dark-matter constraints, the Bose–Einstein condensation thresholds, and the cross-section threshold $\\sigma_{th}$ that the binary calculation improves.","marker":"[11]"},{"why":"Demonstrates the same collapse-based constraint idea applied to GW190425, the comparison case for binary neutron star versus black hole mergers.","marker":"[12]"},{"why":"Identifies GW170817 as a binary neutron star merger through gravitational waves and electromagnetic counterparts, the event used for the new bound.","marker":"[80]"},{"why":"Supplies the GW170817 host-galaxy properties and the inferred merger delay time of 6.8–13.6 Gyr used in the constraint.","marker":"[91]"},{"why":"Provides the Monte Carlo scattering setup and the dynamical-friction hardening rate that the amplification simulations and the mass-fraction cap build on.","marker":"[93]"},{"why":"Gives the gravitational-wave inspiral time for circular binaries used to integrate the capture rate over the binary lifetime.","marker":"[107]"},{"why":"Provides the central dark-matter density of the GW170817 host galaxy that motivates the assumed 0.1 GeV/cm$^3$ density.","marker":"[121]"}],"fun_headline_variants":["Binary neutron stars amplify dark matter capture 4-5 fold","GW170817 tightens dark matter bounds via neutron star binary boost","Binary neutron stars capture 4x more dark matter, tighter bounds","Dark matter accretion boosted 4-5x in neutron star binaries","Neutron star binaries amplify dark matter flux, improving GW170817 bounds"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The calculation assumes the neutron stars spend essentially their whole lifetime as a circular binary whose separation shrinks only by gravitational waves and dark-matter dynamical friction; if real merger progenitors have significant eccentricity, form late from previously isolated neutron stars, or are hardened by other environmental processes, the integrated amplification and the derived constraints shift.","fun_headline_variants_meta":{"raw":{"variants":["Binary neutron stars amplify dark matter capture 4-5 fold","GW170817 tightens dark matter bounds via neutron star binary boost","Binary neutron stars capture 4x more dark matter, tighter bounds","Dark matter accretion boosted 4-5x in neutron star binaries","Neutron star binaries amplify dark matter flux, improving GW170817 bounds"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001477,"raw_usage":{"total_tokens":5955,"prompt_tokens":984,"completion_tokens":4971,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":600,"completion_tokens_details":{"reasoning_tokens":4879}},"tokens_in":600,"tokens_out":4971,"duration_ms":32838,"temperature":1.0,"reasoning_tokens":4879,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T16:33:03.294412+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the eccentricity distribution of merging neutron-star binaries from gravitational-wave waveforms: if a substantial fraction merge with eccentricity larger than about 0.1 at frequencies where the orbital velocity is comparable to the dark-matter velocity dispersion, the circular-orbit integration fails as a description of the accretion history. Equivalently, a single old neutron-star merger in a high-dark-matter-density host whose delay time and dark-matter profile are well measured would let the bound be recomputed from event-specific parameters rather than the assumed 0.1 GeV/cm$^3$ and 220 km/s values, checking whether the factor-of-four tightening survives.","supporting_citations":[{"cited_title":"NGC 4993, the shell galaxy host of GW170817: constraints on the recent galactic merger","cited_arxiv_id":"1801.01493","evidence_quote":"Provides the central dark-matter density of the GW170817 host galaxy that motivates the assumed 0.1 GeV/cm$^3$ density."}],"review_version":2}