{"id":"5526ccd0-b913-4494-90e9-9515ad4cc179","arxiv_id":"1908.10440","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"The existence of PSR J1748-2446ad and PSR B1937+21 excludes ultralight scalars coupled to neutrons with masses near 10^-11 eV at couplings as small as 10^-6 times gravity.","lead":"This paper works out how rapidly spinning neutron stars (millisecond pulsars) can lose energy by emitting ultralight particles through superradiance, and uses the two fastest known pulsars to place new bounds on such particles. If correct, it improves existing fifth-force limits by up to three orders of magnitude and may explain why no pulsar spins faster than about 700 Hz.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"PSR J1748-2446ad constraint hinges on an unmeasured 3×10^8 yr age; the measured lower bound of 2.5×10^7 yr weakens the excluded coupling by a factor of ~3.5 and undercuts the headline 3-order improvement.","rationale":"The reader's weakest assumption identifies the same load-bearing parameter: the adopted pulsar lifetime τ = 3×10^8 yr for PSR J1748-2446ad. This is the most direct place where an input to the exponential-growth criterion is both poorly measured and quantitatively significant. The paper notes the measured lower bound of 2.5×10^7 yr and then substitutes an indirect estimate of ~10^9 yr based on typical millisecond-pulsar magnetic fields, calling the resulting 3×10^8 yr choice 'conservative' even though a larger τ strengthens the bound. Since the excluded coupling scales as τ^(-1/2), the factor-12 difference shifts the boundary by about 3.5 in ϵ. That is comparable to the claimed three-order improvement, so the headline number is not stable under the paper's own quoted uncertainty. This does not invalidate the overall argument: even with the shorter age, the constraints may still improve on torsion balance by roughly two orders of magnitude, and the underlying superradiance mechanism is physically plausible. But the conditional verdict is appropriate until the age or an equivalent upper bound on the spin-down time is established, or until the analysis is recast as an integral over a realistic spin-up/spin-down history.","tokens_in":20739,"tokens_out":37003,"duration_ms":420431,"concrete_test":"Recompute the ψ322 exclusion curve in Fig. 2 for PSR J1748-2446ad, replacing τ = 3×10^8 yr with the measured spin-down lower bound τ = 2.5×10^7 yr while keeping all other inputs fixed. If the lower edge of the excluded region rises above 5×10^-23 at any scalar mass, the claimed three-order improvement over torsion balance is not robust; if it remains below 5×10^-23, the age assumption is not the limiting uncertainty and the conditional verdict can be relaxed.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central exclusion contours in Figs. 2 and 3 are set by requiring the superradiant growth factor Γτ to saturate the available angular momentum, Γτ ≲ ln(L_s/ℏ) ≈ 176. The paper takes τ = 3×10^8 yr for both pulsars, but for PSR J1748-2446ad the direct spin-down measurement gives only a lower bound of 2.5×10^7 yr, as the paper itself states in Section IV B. Because Γ is proportional to ϵ², replacing 3×10^8 yr with the measured lower bound raises the minimum excluded coupling by a factor of √(12) ≈ 3.5. The abstract's headline claim of improving torsion-balance bounds by three orders of magnitude is therefore not a robust consequence of measurement: it depends on an indirect age estimate inferred from a typical magnetic-field strength rather than on the pulsar's observed spin-down. The same fragility affects other hand-picked parameters (T = 10^7 K, Γ1 = 10 Hz, the phonon overlap integral), but the age is the clearest case because the paper itself supplies a conflicting, directly measured number.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes that rotational superradiance of light scalar fields gravitationally bound to millisecond pulsars can be an efficient spin-down mechanism if the scalar has a Yukawa coupling to neutrons. The authors derive a superradiance rate by modeling the stellar medium as a harmonic chain of neutrons, estimating the phonon-excitation absorption coefficient, and then converting it to a growth rate under the superradiance condition. Using the existence of the two fastest known pulsars, PSR J1748-2446ad (716 Hz) and PSR B1937+21 (642 Hz), they place upper bounds on the Yukawa coupling for scalar masses near 10^-11 eV, claiming an improvement over torsion-balance fifth-force bounds by up to three orders of magnitude. The paper also analyzes mode stability against mixing from free precession, equatorial ellipticity, stellar quakes, and companion stars, and speculates that a new particle could explain the apparent absence of pulsars above about 700 Hz.","tokens_in":21045,"tokens_out":10811,"duration_ms":114841,"significance":"If the rate estimate is reliable, the proposal is significant: pulsar spin frequencies are measured with extraordinary precision, so a superradiance-based probe avoids the systematic uncertainties of black-hole spin measurements. The paper's strengths are the explicit physical mechanism, the detailed treatment of damping via mixing with absorptive modes, and the use of independent pulsar observations as the test data. The bounds would probe ultra-light scalars with nucleon couplings far below current laboratory limits and would give a falsifiable prediction for the pulsar frequency distribution. The main caveat, reflected in the major comments below, is that the quantitative exclusion contours inherit order-of-magnitude uncertainties from the phonon-model calculation and from the assumed pulsar age; these issues do not invalidate the concept but they do affect the advertised factor-of-1000 improvement.","major_comments":[{"comment":"The text states that the measured spin-down of PSR J1748-2446ad gives only a lower bound on its characteristic lifetime of 2.5e7 years, but then says 'we conservatively take the stellar lifetime to be tau = 3e8 years for each pulsar in setting our constraints.' This wording is not conservative for exclusion purposes: the superradiance constraint is Gamma*tau <~ ln(L_s/hbar) ~ 176, so the minimum excluded coupling scales as epsilon ~ tau^{-1/2}. Replacing the measured lower bound of 2.5e7 years with the assumed 3e8 years weakens the J1748 bound by a factor of sqrt(12) ~ 3.5 in epsilon. Because the abstract's headline claim of a three-order-of-magnitude improvement rests on combining the two fastest pulsars, the authors should recompute the J1748 contours with tau = 2.5e7 years and quote the weakened bound, or explicitly state which exclusion regions depend on the indirect age estimate rather than on measured spin-down. If the B1937+21 constraint alone preserves the full three-order improvement, that should be stated and supported by its measured characteristic age of about 2e8 years.","section":"Section III A 1, Eqs. (27)-(28)"},{"comment":"The central rate formula is obtained by starting from a 1D harmonic chain of N neutrons and then 'extrapolating to 3D' by taking N to be the total number of neutrons and writing the integral with n(r). This step is not justified. In the 1D chain the phonon frequencies are omega_j ~ (j/N) omega, so the lowest mode frequency is set by the chain length in a way that does not map directly onto the global l = 0, 1, 2 stellar oscillation frequencies (omega_1 ~ 2pi x 2-4 kHz) used later in the paper. The relation between the microscopic chain normalization y_js ~ N^{-1/2} and the normalized displacement eigenfunctions of a realistic neutron-star oscillation mode is not established. In addition, the reduction from Eq. (27) to Eq. (28) assumes that the phonon wavefunction has the same angular structure as the scalar mode; the text acknowledges that otherwise the integral vanishes for a spherical star. For l = 1 the authors introduce an additional (R/a0)^2 suppression by hand. Since the predicted superradiance rate is proportional to the squared overlap integral in Eq. (28), all exclusion contours in Figs. 2-4 inherit this uncontrolled factor. The authors should either provide a direct matching of the toy-model matrix element to realistic neutron-star phonon/oscillation modes or give an explicit uncertainty estimate for the overlap integral. The same concern applies to the l = 0 and l = 1 damping rates, which the text itself describes as 'admittedly rough' and which enter the psi211 bounds and the mixing cutoffs.","section":"Section III A 1, Eqs. (27)-(28)"},{"comment":"The reported exclusion contours depend on several unmeasured or loosely constrained astrophysical inputs: the temperature T = 1e7 K, the radius R = 12 km, the mass of PSR B1937+21 (taken as 1.4 M_sun without a measurement), the phonon frequencies, and the damping rates. The superradiance rate in Eq. (28) scales linearly with T and quadratically with the overlap integral, so these choices propagate directly into the excluded coupling. The paper does not provide a sensitivity analysis, and the claimed three-order-of-magnitude improvement over torsion-balance limits is a quantitative statement that requires such an analysis. At minimum, the authors should show how the contours in Figs. 2-4 shift when T is varied over the quoted 5e5-1e8 K range and when R and the unmeasured mass of B1937+21 are varied within their plausible ranges.","section":"Section IV B (parameter choices and uncertainty propagation)"}],"minor_comments":[{"comment":"The caption refers to 'PSR B1937-21' while the text and Figure 2 use 'PSR B1937+21'; make the notation consistent.","section":"Figure 3 caption"},{"comment":"There are several typographical errors: 'Millsecond', 'equillibrium', 'absorpative', and 'disovered' should be corrected.","section":"General"},{"comment":"Reference [21] appears to be a webpage on sodium atomic data that is unrelated to the surrounding discussion; it should be removed or replaced with a relevant citation.","section":"References"},{"comment":"The conclusions state that the constraints improve current bounds by 'two to four orders of magnitude', while the abstract and Section IV B say 'three orders of magnitude'; harmonize these quantitative claims.","section":"Section V vs. abstract"},{"comment":"The suggestion that a scalar with mass ~2pi x 1500-3000 Hz would produce a cutoff near 700 Hz requires high-m superradiant modes, whose overlap with the star is strongly suppressed; the paper should either provide an estimate showing this is plausible or clearly label the suggestion as an unsupported speculation.","section":"Conclusions (pulsar cutoff speculation)"}],"recommendation":"major_revision","confidential_remarks":"The manuscript proposes an interesting and potentially important new probe, but the quantitative bounds need revision. The most urgent issue is the age assumption for PSR J1748-2446ad in Section IV B: the paper uses 3e8 years while itself reporting a measured lower bound of 2.5e7 years, and this choice directly affects the headline improvement factor. The 1D-to-3D extrapolation in Eqs. (27)-(28) also needs better justification or an explicit uncertainty estimate. If the authors can show that the B1937+21 constraint alone preserves the claimed three-order improvement, or if they recompute the J1748 contours with the measured age, the paper would be much stronger."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Here's my read. The genuinely new content is the phonon-mediated absorption calculation for a Yukawa-coupled scalar and the first real treatment of mixing with absorptive modes. Earlier pulsar-superradiance papers only considered stellar conductivity; this paper adds a concrete dissipation mechanism and thinks carefully about precession, ellipticity, quakes, and companions that could kill the instability. That is a real step forward, and the paper is honest about what is rough.\n\nThe central physics is plausible. The superradiance rate comes from an absorption calculation with a 1D toy model extrapolated to 3D, and the l=0,1 phonon damping rates are admittedly rough. The constraints should be read as order-of-magnitude, which the paper mostly does, but it doesn't propagate any uncertainties—everything is point estimates. That is my main technical complaint.\n\nThe biggest soft spot is the pulsar age. For PSR B1937+21, τ=3×10^8 yr is consistent with the measured characteristic age, so the bound there is fine. For PSR J1748-2446ad, the measured lower bound is 2.5×10^7 yr, and the paper adopts 3×10^8 yr via a typical-field estimate, calling it conservative. That's wrong in direction: a longer age gives a stronger exclusion, so the choice strengthens the headline result. Using the measured lower bound raises the excluded coupling by about sqrt(12)≈3.5. The torsion-balance improvement is still maybe two orders of magnitude, not three. This is addressable—show the bounds with both ages—but it should be fixed before the paper is taken as a precision result.\n\nThe 700 Hz cutoff possibility is clearly speculative, and the paper says so; it's not load-bearing.\n\nNo circularity; the derivation is independent of the pulsar observations used for the bounds. The citation pattern looks fair, including to the prior pieces it builds on.\n\nRecommendation: send it to peer review. The idea is interesting, the formalism is a real contribution, and the soft spots are fixable in revision. Just make sure the age dependence is presented honestly.","headline":"A genuinely new phonon-mediated superradiance calculation for millisecond pulsars, worth refereeing, but the headline bound rests on an optimistic age choice for the fastest pulsar.","tokens_in":21553,"tokens_out":2907,"would_cite":true,"duration_ms":29249,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Rotational superradiance in millisecond pulsars can turn the existence of the two fastest pulsars into bounds on ultralight scalars that beat torsion-balance experiments by three orders of magnitude.","keywords":["superradiance","millisecond pulsars","ultralight scalar bosons","fifth-force constraints","QCD axion","neutron star phonons","Yukawa coupling","pulsar spin-down"],"falsifier":"Measure the spin-down age of PSR J1748-2446ad directly; if it is close to 25 million years rather than the assumed 300 million years, the derived superradiance limit weakens by roughly an order of magnitude and the excluded region shrinks.","tokens_in":20521,"feed_emoji":"🌀","tokens_out":11425,"duration_ms":93631,"temperature":0.7,"pith_summary":"The paper argues that rotational superradiance—the same instability that can drain angular momentum from spinning black holes—can operate efficiently in millisecond pulsars, because the dense neutron medium absorbs the emitted particles through phonon excitations. It then uses the existence of the two fastest known pulsars, PSR J1748-2446ad and PSR B1937+21, as an exclusion device: had a scalar of mass near the rotation frequency (about $10^{-11}$ eV) coupled to neutrons through $\\epsilon\\,\\Psi\\,\\bar{n} n$ with strength above roughly $10^{-6}$ times gravity, it would have spun these pulsars down within their lifetime. The resulting bounds beat torsion-balance fifth-force limits by up to three orders of magnitude. If correct, the argument turns ordinary pulsar timing into a particle-physics laboratory and offers a new-physics explanation for the puzzling absence of pulsars spinning above about 700 Hz.","feed_headline":"Two pulsars rule out scalar forces a million times weaker than gravity","feed_subtitle":"Superradiance in two millisecond pulsars beats torsion-balance fifth-force bounds by up to 3 orders of magnitude.","key_machinery":"The engine is rotational superradiance of gravitationally bound scalar modes around a rotating neutron star, with absorption supplied by phonon excitation through the Yukawa interaction $\\epsilon\\,\\Psi\\,\\bar{n} n$. The growth rate is $\\Gamma_{nlm}=C_{nlm}(\\mu-m\\Omega)/\\mu$, so when $\\mu<m\\Omega$ the absorptive term becomes emissive; efficiency peaks at $\\mu\\sim\\Omega$ because higher angular-momentum modes have suppressed overlap with the star. The authors compute the absorption coefficient $C_{nlm}$ from a scalar-to-phonon conversion rate in the stellar medium, then check that astrophysical asymmetries—free precession, equatorial ellipticity, stellar quakes, and companion tides—do not mix the superradiant mode into absorptive modes fast enough to stop its growth, except in the regions where the exclusion plots are cut off.","core_discovery":"The central claim is that a light scalar $\\Psi$ with a Yukawa coupling $\\epsilon\\,\\Psi\\,\\bar{n} n$ to neutrons can be superradiantly amplified in the gravitational bound states around a rapidly rotating neutron star. When the scalar mass $\\mu$ is close to the stellar rotation frequency $\\Omega$, the lowest angular-momentum modes ($\\psi_{211}$ and $\\psi_{322}$) overlap strongly with the star, and the coupling excites phonons that dissipate energy; in the rotating frame the absorption term flips sign into an emissive term. The paper shows that for $\\mu\\sim 10^{-11}$ eV the growth rate is fast enough that the observed 716 Hz and 642 Hz pulsars, taken with ages around $3\\times10^8$ years, exclude couplings down to roughly $10^{-6}$ of gravitational strength, improving on torsion-balance fifth-force bounds by up to three orders of magnitude. For certain neutron-star equations of state, the same argument would also rule out the QCD axion with a Planck-scale decay constant in a mass window near $10^{-12}$ eV.","pith_inferences":["If the excluded band is real, laboratory searches for axion-like particles should concentrate on the $\\sim10^{-11}$ eV mass window, where the astrophysical bound is strongest; a direct detection there would turn the two pulsars into a calibrated probe of scalar-neutron couplings.","The argument could be inverted as an age diagnostic: for a scalar with coupling just below the exclusion boundary, the observed spin frequency of a millisecond pulsar would encode how long it has been spinning, making superradiance a clock rather than only a constraint.","Because the low-mass edge of the PSR J1748-2446ad exclusion is set by tidal mixing from its companion, a precise measurement of that companion's orbit would sharpen or weaken the boundary; better measurements of the low-order phonon damping rates would similarly tighten the $\\psi_{211}$ bounds."],"forward_implications":["Scalars with masses between about $2\\times10^{-12}$ and $6\\times10^{-12}$ eV and Yukawa couplings to neutrons near $10^{-6}$ of gravity are excluded by the two pulsars, improving torsion-balance bounds by up to three orders of magnitude.","If neutron-star equations of state produce $\\theta_{\\rm eff}\\sim1$, the QCD axion with mass between $5\\times10^{-13}$ and $3\\times10^{-12}$ eV and a Planck-scale decay constant is ruled out.","The observed absence of pulsars above about 700 Hz, despite equations of state allowing rotation up to about 1500 Hz, could be explained by a scalar of mass near $10^{-11}$ eV coupled to nucleons.","Superradiant braking would predict a pile-up of pulsar spin frequencies near half the particle mass, a signature that distinguishes it from gravitational-wave or r-mode spin-down.","A hypothetical isolated pulsar rotating at 1200 Hz would extend the excluded region to larger masses, showing that faster pulsar discoveries sharpen the particle-physics reach."],"supporting_citations":[{"why":"Supplies the observed 716 Hz rotation of PSR J1748-2446ad, its measured mass, and its companion parameters used in the constraints.","marker":"[17]"},{"why":"Supplies PSR B1937+21, the accretion-based formation picture for millisecond pulsars, and the nominal neutron-star parameters used in the bounds.","marker":"[18]"},{"why":"Supplies the torsion-balance fifth-force bounds that the superradiance constraints are compared against and exceed.","marker":"[5]"},{"why":"Provides the phonon-mode frequencies and gravitational-wave damping rates that set the scalar absorption rate in the star.","marker":"[25]"},{"why":"Provides the neutron-star equations of state used to argue for a pseudoscalar condensate with $\\theta_{\\rm eff}\\sim1$, connecting the bounds to the QCD axion.","marker":"[31]"},{"why":"Supplies the maximum equatorial ellipticity that bounds mixing-induced damping of the superradiant mode.","marker":"[22]"},{"why":"Establishes superradiance as a general instability of rotating absorptive systems, the mechanism the paper applies to pulsars.","marker":"[14]"}],"fun_headline_variants":["Superradiant pulsars rule out weak couplings down to 1e-6 of gravity","Pulsar superradiance beats fifth-force bounds by 1000x","Fast pulsars constrain axions and ultralight scalars","Millisecond pulsars probe forces 1e6 times weaker than gravity","Superradiant pulsars exclude new forces 1000x better than lab tests"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The constraints assume both pulsars have been spinning at their observed frequencies for about 300 million years; the faster pulsar's measured lower age bound is only 25 million years, so if its true age is nearer the lower bound the excluded region shrinks.","fun_headline_variants_meta":{"raw":{"variants":["Superradiant pulsars rule out weak couplings down to 1e-6 of gravity","Pulsar superradiance beats fifth-force bounds by 1000x","Fast pulsars constrain axions and ultralight scalars","Millisecond pulsars probe forces 1e6 times weaker than gravity","Superradiant pulsars exclude new forces 1000x better than lab tests"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000869,"raw_usage":{"total_tokens":3766,"prompt_tokens":949,"completion_tokens":2817,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":565,"completion_tokens_details":{"reasoning_tokens":2716}},"tokens_in":565,"tokens_out":2817,"duration_ms":19307,"temperature":1.0,"reasoning_tokens":2716,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T10:43:16.489000+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the spin-down age of PSR J1748-2446ad directly; if it is close to 25 million years rather than the assumed 300 million years, the derived superradiance limit weakens by roughly an order of magnitude and the excluded region shrinks.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the torsion-balance fifth-force bounds that the superradiance constraints are compared against and exceed."}],"review_version":1}