{"id":"805da7b7-f7e2-4420-ad9b-fb9e18177d19","arxiv_id":"2507.06178","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"Dark matter accretion onto binary pulsars is far too weak to affect observed orbital decay, so existing pulsar timing data cannot probe dark matter microphysics.","lead":"Dark matter falling onto the neutron stars of a binary pulsar changes the stars' masses and could in principle alter the binary's orbital decay. The paper finds this effect is many orders of magnitude weaker than gravitational wave emission for pulsars near Earth, so current observations cannot constrain dark matter.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The derivation of Eq. (2) silently assumes captured DM adds no orbital energy or angular momentum; this does not threaten the local null result but could undermine the galactic-center projection.","rationale":"The paper presents a clean, mostly self-contained calculation: the multiscatter capture formalism is standard, the use of published pulsar timing data is appropriate, and the geometric-limit bound makes the local negative result very secure. The reader's CONDITIONAL verdict is well calibrated. The most load-bearing soft spot is indeed the unstated assumption that accreted DM does not inject orbital energy or angular momentum when deriving Eq. (2). This does not affect the local null result by more than an order-one factor on an already negligible term, but the galactic-center projection, which is the forward-looking part of the paper, depends on the same simplified formula. The paper also contains a secondary sign inconsistency: the text repeatedly says DM accretion 'enhances the orbital decay rate', while Eq. (5) shows a positive contribution that would make |Pdot/P| smaller. This should be corrected but does not change the quantitative comparison. Our concrete test would settle whether the neglected energy/angular-momentum injection materially changes the projected DM term at the galactic center; if it does, the paper should be revised to state the assumption and temper the projection. Since the reader's conditional verdict already demands such revisions, no change to the verdict is needed from this stress-test pass.","tokens_in":13473,"tokens_out":13255,"duration_ms":147420,"concrete_test":"Re-derive Eq. (2) including the accreted DM's energy and angular momentum. Concretely, compute the specific orbital energy and specific angular momentum carried by a captured DM particle at each pulsar's position, using the halo velocity distribution (u0 = 270 km/s) and the binary orbital velocity; add the resulting ˙E_acc and ˙L_acc to Eqs. (2) and (6), then re-evaluate the galactic-center benchmark (ρχ = 1e18 GeV/cm^3, mχ = 1e5 GeV, σχN = 1e-46 cm^2). If the corrected DM contribution to Pdot/P changes by more than 50% relative to Table II's quoted value, the projection is unreliable; if it remains within the same order, the concern is not load-bearing.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central null result is robust: even at the geometric capture limit, Fig. 2 and Table II give |Pdot/P|_DM ~ 1e-33 s^-1 versus |Pdot/P|_GW ~ 1e-16 s^-1, so local binary pulsars cannot constrain DM microphysics. The load-bearing weakness is the derivation of Eq. (2). Starting from P = 2π G M μ^{3/2} (-2E)^{-3/2}, the time derivative is taken and then ˙E is set exclusively to the GW quadrupole loss. This silently assumes that each captured DM particle contributes zero orbital energy and zero orbital angular momentum to the binary. If captured DM carries halo kinetic energy (mχ u^2/2 with u ~ 270 km/s) or specific angular momentum comparable to the pulsar's orbital velocity, there is an additional ˙E and a change in orbital angular momentum (which feeds into eccentricity and the GW term). For local pulsars the correction is at most comparable to the already-tiny DM term, so the 'no constraint' conclusion stands. However, the paper's forward-looking galactic-center projection uses the same simplified Eq. (5) and would become comparable to GW at ρ ~ 1e18 GeV/cm^3; an order-one correction from the neglected energy/angular-momentum injection could change the projected signal strength or even its sign. The assumption is never stated or justified, and it is not a standard limit unless accreting matter is added at rest with zero specific angular momentum. This is the most load-bearing weakness for the paper's broader claim that orbital decay of accreting pulsars can probe DM.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies whether dark matter (DM) accretion onto the neutron stars of a binary pulsar can measurably change the orbital period derivative. It derives an expression for Pdot/P from the Kepler relation, Eq. (1), adding a DM mass-accretion term to the Peters-Mathews gravitational-wave (GW) term, Eq. (6). Using the multiscatter DM capture formalism of Bramante et al. and the Asteria package, it computes the DM contribution for constant and velocity-dependent cross-sections. Comparing with four known binary pulsars, the paper finds the DM contribution is at least fifteen orders of magnitude below the GW contribution (Table II) and concludes that current pulsar timing data cannot constrain DM particle properties. It then suggests that binary pulsars near the galactic center, where the DM density is much higher, could probe DM in the future.","tokens_in":13675,"tokens_out":15349,"duration_ms":161422,"significance":"If the central result is correct, the paper provides a robust negative result: DM accretion onto the neutron stars of known binary pulsars is far too slow to affect the orbital decay rate, even when the capture rate is artificially pushed to the geometric limit. This is a useful null result because it shows that multiscatter capture saturates and does not yield the large effects one might guess from single-scatter capture. The paper also correctly emphasizes that the four pulsars studied are relatively far from the galactic center. The forward-looking galactic-center projection, however, is not yet firmly grounded because the orbital evolution formula used for the projection omits the energy and angular momentum carried by the accreted DM particles; the assumptions behind that omission are not stated or justified. With that caveat, the numerical null result for local pulsars is well supported.","major_comments":[{"comment":"The derivation of Eq. (2) sets dE/dt equal to the Peters-Mathews GW loss only, without stating the assumption that the accreted DM contributes zero orbital energy and zero orbital angular momentum. This assumption is not benign: for mχ = 10^5 GeV and u0 = 270 km/s, each captured DM particle carries kinetic energy about 0.04 GeV, so the associated injection rate dE/dt leads to a Pdot/P correction |(3/2) dE/dt / |E|| that is comparable to the quoted DM mass term (both of order 10^-33 s^-1 for the local pulsars). At the galactic-center density ρχ = 10^18 GeV/cm^3 used in Sec. IV C, the same correction becomes of order 10^-15 s^-1, comparable to the projected DM signal. The galactic-center projection is therefore not robust until the full energy and angular-momentum accounting is performed, or until the zero-injection assumption is explicitly stated and defended. The text in Sec. II even contradicts the assumption by claiming that DM accretion 'increases the orbital decay rate by transferring angular momentum to the binary,' although no angular-momentum term appears in Eq. (2).","section":"II, Eqs. (2) and (5)"},{"comment":"The abstract says DM accretion 'may also modify the orbital evolution by enhancing the orbital decay rate,' and Sec. II says DM accretion 'increases the orbital decay rate.' This is the opposite of the sign in Eq. (2): the DM mass term is positive, while the observed and GW terms are negative, so DM accretion makes Pdot less negative, i.e., it slows the orbital decay rather than enhancing it. The wording should be corrected to say that accretion opposes the GW-driven decay, or reduces the magnitude of the decay rate.","section":"Abstract and II"},{"comment":"The galactic-center projection is presented as a quantitative statement ('DM accretion will influence the rate of orbital period'), but it uses Eq. (5), which is derived under the same unstated zero-energy/zero-angular-momentum assumption as Eq. (2). In addition, the projection assumes ρχ = 10^18 GeV/cm^3 without discussing whether such a spike density is compatible with the 'missing pulsar problem' that the paper itself mentions in Sec. V. The projection should be framed as a speculative order-of-magnitude illustration, not as a firm prediction, until the accretion-injection terms are included.","section":"IV C"}],"minor_comments":[{"comment":"The distance listed for PSR B1913+16 is 0.34 kpc, while the published value is about 3.4 kpc (Weisberg & Huang 2016). The distance is not used in the calculation, but the entry should be corrected.","section":"Table I"},{"comment":"The text says 'σχN = 10−48 cm^2 for α = 1, and σχN = 10−48 cm^2 for α = 2 model', but Fig. 1 and the surrounding discussion indicate the quartic model uses σχN = 10−52 cm^2. One of these is a typo.","section":"IV A"},{"comment":"The sentence 'for lower cross sections, the capture rate saturates to a geometrical rate and is reduced for higher DM masses' is unclear; saturation occurs at high cross-section, not low cross-section. Please rephrase.","section":"IV A"},{"comment":"The sentence 'From Eq. (6), we find that gravitational wave emission decreases the orbital decay rate' is imprecise: GW emission makes Pdot more negative, i.e., it accelerates the decay. The intended meaning is that GW emission decreases the orbital period.","section":"II"},{"comment":"The notation 'N = (10, eτ)' in the restriction of the sum limit is ambiguous; please specify whether the maximum N is max(10, e^τ) or min(10, e^τ), and cite the corresponding comment [84] precisely.","section":"III"},{"comment":"The approximation σ(vrel) ≈ σχN (w0/uref)^{2α} with w0 = sqrt(u0^2+v_esc^2) replaces a velocity-dependent integrand with a fixed value and discards the velocity distribution f(u). Since the paper's claims about positive velocity dependence rest on this approximation, a validation against the full integration, or at least a statement of whether the approximation over- or under-estimates the capture rate, should be included.","section":"III"},{"comment":"The note added states that the expression for Pdot/P differs from that of Ref. [97], but the difference is not shown. A brief comparison would help readers understand the distinction between the single-scatter and multiscatter treatments.","section":"V, Note added"}],"recommendation":"major_revision","confidential_remarks":null},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Bottom line: this is a solid negative-result paper. The central calculation says that for the four binary pulsars with measured Pdot, even the geometric upper bound on DM accretion contributes |Pdot/P| ~ 10^-33 s^-1, fifteen orders of magnitude below the observed ~10^-17 to 10^-16 s^-1. That conclusion is robust, and the paper deserves serious consideration.\n\nWhat is genuinely new is the combination of multiscatter capture with the orbital-period derivative, together with the explicit contrast to the single-scatter work of [97]. The note added is candid and correct: single-scatter capture grows with sigma, so it can produce apparent constraints that disappear once you treat large cross-sections properly.\n\nThe paper does several things well. It uses published timing data for four well-measured pulsars, tabulates the contributions, and is explicit that the geometric limit is a safe upper bound. The limitations section lists the main simplifications: non-rotating neutron star, constant geometric cross-section, DM rest frame. Using the Asteria package is a plus for reproducibility.\n\nNow the soft spots, in proportion.\n\nFirst, the prose has a sign error that should not survive review. The abstract and Section II say DM accretion 'enhances the orbital decay rate,' but Eq (2) gives the DM term a positive sign: accretion increases M, which for fixed E increases P, i.e., it slows the GW-driven decay. The same verbal error appears in the discussion of torque balance. The equations themselves appear consistently positive on the DM terms, so this is a wording/interpretation issue, not a computational one, but it is confusing and needs fixing.\n\nSecond, and more substantive, the derivation of Eq (2) is the one-parameter mass-growth formula: it takes P(E, M1, M2), differentiates, and inserts only the GW energy loss for dE/dt. This silently assumes each captured DM particle contributes zero orbital energy and zero angular momentum to the binary. The stress-test note is right that this is unstated. For the local pulsars it does not matter: a factor of a few or even an order of magnitude leaves the 15-order gap intact. But the galactic-center projection uses the same formula at rho ~ 10^18 GeV/cm^3, where the DM term becomes comparable to the GW term. There, the sign and magnitude of the missing energy/angular-momentum terms could qualitatively change the projected signal. The paper should either include those terms or explicitly state the assumption and justify it (e.g., isotropy of the DM halo, spherical capture).\n\nThird, the GC projection is explicitly speculative and depends on DM spikes and on pulsars existing there; the paper acknowledges the missing-pulsar problem. That is fine, but it should be labeled as an order-of-magnitude illustration rather than a prediction.\n\nWho should read this: the compact-object DM community and anyone comparing single- vs multi-scatter capture in binaries. It is not a breakthrough, but it is a useful boundary on a proposed channel. I would send it to a serious referee; with the sign wording fixed and the accretion assumption stated, it can be published.","headline":"Solid null result for DM accretion in binary pulsars; a sign wording fix and an explicit accretion assumption are needed before publication.","tokens_in":14307,"tokens_out":5279,"would_cite":false,"duration_ms":51909,"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":"Dark matter accretion shifts binary pulsar periods by at least fifteen orders of magnitude less than gravitational waves, so current timing data cannot constrain dark matter.","keywords":["dark matter capture","binary pulsars","orbital period decay","pulsar timing","multiscatter capture","velocity-dependent cross sections","galactic center dark matter spike","neutron stars"],"falsifier":"Recompute the dark-matter term for the Table I pulsars under the paper's geometric-limit maximum capture rate and $\\rho_\\chi=0.4$ GeV cm$^{-3}$: if the result came out within even a few orders of magnitude of the gravitational-wave term rather than fifteen below it, the no-constraint conclusion would be falsified. Observationally, a binary pulsar whose residual $\\dot P/P$ after subtracting the gravitational-wave prediction grows with ambient dark-matter density at a level near $10^{-19}$ s$^{-1}$ or above would contradict the paper's claim that the effect is unobservably small.","tokens_in":13097,"feed_emoji":"🌌","tokens_out":9664,"duration_ms":106620,"temperature":0.7,"pith_summary":"This paper asks whether the rate at which binary pulsar orbits shrink can reveal dark matter. It derives how dark-matter accretion onto each neutron star enters the orbital period derivative, computes the accretion rate with multiscatter capture, and compares the result with four precisely timed binary pulsars near Earth. The answer is that the dark-matter term is at least fifteen orders of magnitude smaller than the gravitational-wave term, so existing observations cannot place meaningful constraints on dark-matter particle properties. The paper argues that a binary pulsar merging inside a high-density dark-matter spike near the galactic center could bring the effect into view, though such systems have not yet been found.","feed_headline":"Dark matter term in pulsar orbits is 15 orders too small to see","feed_subtitle":"Existing pulsar timing can't probe dark matter, but galactic-center mergers could.","key_machinery":"The carrying mechanism is the logarithmic derivative of the Kepler period formula, giving $\\dot P/P = -\\tfrac{3}{2}\\dot E/E + (1+\\tfrac{M_2}{2M})\\dot M_1/M_1 + (1+\\tfrac{M_1}{2M})\\dot M_2/M_2$, which converts two competing effects—energy loss to gravitational waves and mass gain from dark-matter capture—into a single period-change ratio. The paper feeds this formula with the standard quadrupole gravitational-wave energy-loss rate and a multiscatter dark-matter capture rate built from a Maxwellian halo velocity distribution, Poisson-distributed scattering probabilities, and optical-depth-dependent approximations, extended to velocity-dependent cross-sections $\\sigma\\propto v^{2\\alpha}$ with $\\alpha=1,2$. That machinery produces the numerical dark-matter entries that are compared with observed pulsar timing data.","core_discovery":"The central claim is that for the four binary pulsars with precise timing—B1534+12, B1913+16, J0737-3039, and J1757-1854—the contribution of dark-matter accretion to the normalized orbital period derivative is about $1.8\\times10^{-33}$ s$^{-1}$, while the gravitational-wave contribution is roughly $1\\times10^{-16}$ s$^{-1}$, a gap of at least fifteen orders of magnitude. Even at the geometric-limit maximum capture rate, and even with velocity-dependent cross-sections that enhance capture at low dark-matter mass and large cross-section, the dark-matter term remains far below the timing residuals. The paper concludes that current pulsar timing data cannot constrain dark-matter microphysics through this mechanism, but that a binary pulsar near the galactic center, where dark-matter density could reach $10^{18}$ GeV cm$^{-3}$, could in principle make the effect detectable.","pith_inferences":["The same logarithmic-derivative formula could in principle be applied to white-dwarf or black-hole binaries, where capture rates and dark-matter densities differ, a step the paper does not take.","The galactic-center projection assumes that the linear scaling of capture rate with dark-matter density holds at spike densities; if the velocity dispersion there differs from the local Maxwellian assumed, the enhancement could be larger or smaller than quoted.","If a binary pulsar is eventually found near the galactic center, the so-called missing pulsar problem means its absence or presence will be as informative about dark-matter spikes as the period-derivative measurement itself."],"forward_implications":["For the four nearby binary pulsars, the dark-matter accretion contribution to $\\dot P/P$ is at least fifteen orders of magnitude below the gravitational-wave contribution, so current timing residuals cannot resolve it even at the geometric-limit capture rate.","Positive velocity dependence of the dark-matter-baryon cross-section boosts capture at low dark-matter mass and large cross-section, but not enough to close the gap at the local dark-matter density $\\rho_\\chi=0.4$ GeV cm$^{-3}$.","Because the capture rate scales linearly with dark-matter density, placing the same binary in a galactic-center spike with $\\rho_\\chi\\sim10^{18}$ GeV cm$^{-3}$ lifts the dark-matter term by roughly the same factor, potentially making it comparable to the gravitational-wave term.","For cross-sections above the geometric limit $\\sigma_{\\rm geo}\\simeq2\\times10^{-45}$ cm$^2$, the single-scatter capture description breaks down, so constraints derived from single-scatter scaling should be treated as unreliable.","The observed orbital decay of the four binaries is consistent with gravitational-wave emission alone; the dark-matter term is far too small to be extracted from the difference between the measured and predicted period derivatives."],"supporting_citations":[{"why":"Supplies the review of dark matter in compact stars and the large-optical-depth capture approximations used in the calculation.","marker":"[5]"},{"why":"Provides the quadrupole gravitational-wave energy-loss formula that gives the dominant contribution to $\\dot P/P$.","marker":"[66]"},{"why":"Establishes the observed binary-pulsar orbital decay that calibrates the gravitational-wave prediction against timing data.","marker":"[70]"},{"why":"Provides the multiscatter stellar-capture formalism that determines the dark-matter accretion rate onto each neutron star.","marker":"[77]"},{"why":"Provides the numerical capture-rate package used to compute the accreted dark-matter mass across interaction regimes.","marker":"[79]"},{"why":"Supplies the timing parameters of PSR B1534+12 used in Table I and Table II.","marker":"[85]"},{"why":"Supplies the timing parameters of PSR B1913+16 used in Table I and Table II.","marker":"[86]"},{"why":"Supplies the timing parameters of PSR J0737-3039 used in Table I and Table II.","marker":"[87]"},{"why":"Supplies the timing parameters of PSR J1757-1854 used in Table I and Table II.","marker":"[88]"},{"why":"Provides the galactic-center dark-matter spike density used for the future-projection argument.","marker":"[89]"}],"fun_headline_variants":["Galactic center pulsars could reveal dark matter via orbital decay","Pulsar timing: dark matter effect 15 orders too faint for now","Dark matter's pull on pulsar orbits undetectable until galactic center","Pulsar orbital decay: dark matter term too small to probe now"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that accreted dark matter changes only the neutron-star masses in the Kepler formula and contributes nothing to the binary's orbital energy or angular momentum; if captured dark matter carried non-negligible orbital angular momentum, the derived dark-matter contribution to $\\dot P/P$ would differ.","fun_headline_variants_meta":{"raw":{"variants":["Galactic center pulsars could reveal dark matter via orbital decay","Pulsar timing: dark matter effect 15 orders too faint for now","Dark matter's pull on pulsar orbits undetectable until galactic center","Pulsar orbital decay: dark matter term too small to probe now"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00079,"raw_usage":{"total_tokens":3441,"prompt_tokens":863,"completion_tokens":2578,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":479,"completion_tokens_details":{"reasoning_tokens":2500}},"tokens_in":479,"tokens_out":2578,"duration_ms":18528,"temperature":1.0,"reasoning_tokens":2500,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T19:09:33.926108+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Recompute the dark-matter term for the Table I pulsars under the paper's geometric-limit maximum capture rate and $\\rho_\\chi=0.4$ GeV cm$^{-3}$: if the result came out within even a few orders of magnitude of the gravitational-wave term rather than fifteen below it, the no-constraint conclusion would be falsified. Observationally, a binary pulsar whose residual $\\dot P/P$ after subtracting the gravitational-wave prediction grows with ambient dark-matter density at a level near $10^{-19}$ s$^{-1}$ or above would contradict the paper's claim that the effect is unobservably small.","supporting_citations":[],"review_version":1}