{"id":"266baf3e-6520-4b1d-b7de-57d59296d8b8","arxiv_id":"1908.10019","paper_version":1,"verdict":"ACCEPT","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"low","formal_verification":"none","parameter_count":0,"one_line_summary":"Using three small-eccentricity relativistic binary pulsars, this paper obtains order-of-magnitude upper limits on SME matter-gravity coefficients for neutrons, protons, and electrons.","lead":"This paper sets new upper limits on Lorentz-violating matter-gravity couplings in the Standard-Model Extension by analyzing timing data from three neutron star-white dwarf binary pulsars. The limits for neutron, proton, and electron coefficients complement existing constraints from lunar laser ranging and atomic clocks.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Main residual concern is the O(1) projection of unmeasured orbital orientation and proper-motion contamination in xp_dot; this is acknowledged and does not overturn the order-of-magnitude maximal-reach limits, but it is the weakest link.","rationale":"The paper's central claim is explicitly framed as maximal-reach limits, and the O(1) projection treatment is acknowledged in Sec. IV with the factor-of-a-few caveat. The order-of-magnitude arithmetic for Eq. (24) is consistent with the quoted limits, and the proper-motion scale for J1738 is comparable to, not far larger than, the observed xp_dot, so the no-conspiracy argument is plausible. The residual concern is a real soft spot for the precise value of the best-limit claim, but it does not change the overall verdict: the limits are defensible as order-of-magnitude, complementary constraints. The reader identified the same weakest assumption, and I agree that it is the most load-bearing caveat while remaining within the paper's stated scope.","tokens_in":11369,"tokens_out":24306,"duration_ms":267646,"concrete_test":"Run a Monte Carlo over the unknown Euler angles (Omega, i, omega) and the proper-motion position angle for PSR J1738+0333, using its timing parameters and Eq. (24); compute the 5th-percentile bound on c^n_jk. If the 5th percentile is within a factor of 3 of the quoted 1e-11, the O(1) treatment is adequate; if it is an order of magnitude weaker, the Table IV best-limit claim must be qualified as orientation-dependent.","verdict_should_be":"UNCHANGED","load_bearing_attack":"In Sec. IV, Eq. (24) converts the measured xp_dot of PSR J1738+0333 into the best limit on c^w_jk in Table IV. The coefficient of A_jl in Eq. (24) is a projection onto the orbital basis (a_hat, b_hat, c_hat), which depends on the unmeasured longitude of ascending node and the inclination. The author treats these projections as O(1) and argues that no Nature's conspiracy introduces more than a factor-of-a-few uncertainty. For a single nonzero Cartesian component, however, the projection is a product of sines and cosines of Euler angles; for a fixed orientation it can be much smaller than 1, and the factor-of-a-few is a prior, not a proven bound. The same applies to B_j in Eqs. (25)-(26), which set the (a_eff)_k and c_0k limits. Separately, the observed xp_dot for J1738 includes a proper-motion contribution whose natural scale is comparable to the measured value, so a partial cancellation could shift the derived limit by an O(1) factor without any fine-tuned conspiracy. These effects do not invalidate maximal-reach limits, but they mean the quoted Table IV numbers should be read as reach under a favorable-orientation assumption.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper derives new constraints on Lorentz-violating matter-gravity couplings in the Standard-Model Extension (SME) using three small-eccentricity neutron star–white dwarf binaries. It combines the osculating-element equations of Jennings et al. (2015) with published pulsar timing measurements of the projected semimajor axis derivative, ˙xp, and the Laplace–Lagrange parameters, ˙η and ˙κ. After a small-eccentricity expansion (Eqs. 21–26), the author converts these measurements into maximal-reach limits on the SME coefficients c^w_jk, c^w_0k, and (a^w_eff)_k for neutrons, protons, and electrons. The main results are collected in Table IV, with the best limits on c^w_jk coming from PSR J1738+0333 and the best limits on c^w_0k and (a^w_eff)_k from PSR J0751+1807. The analysis relies on treating the unknown orbital orientation (in particular the longitude of ascending node, Ω) as O(1) in the projections, and on using measured ˙xp values as upper limits despite a comparable proper-motion contribution. The paper explicitly acknowledges these simplifications and discusses caveats related to strong-field effects in neutron stars.","tokens_in":11648,"tokens_out":10497,"duration_ms":106036,"significance":"If the derived limits hold, they are competitive and complementary to existing bounds: c^n,jk limits reach about 1e-11, c^n,0k about 1e-8, and (a^n_eff)_k about 1e-8 GeV, with the neutron (a_eff)_k bound slightly stronger than the corresponding lunar laser ranging result. A notable strength is that the analysis uses externally published timing measurements rather than fitting the SME coefficients to the data, so there is no circularity in the central claims. The paper is transparent about its maximal-reach methodology and about the strong-field caveat for neutron stars, and it provides a clear scaling law (T_obs^{-1.5}) for future improvements. These features make the paper a useful and reproducible contribution to the experimental SME literature, even though the quoted limits are order-of-magnitude estimates under simplifying assumptions.","major_comments":[{"comment":"The Table IV limits assume that the projections of A_jl and B_j onto the orbital basis (a_hat, b_hat, c_hat) are O(1). For a single nonzero Cartesian component, the projection is a product of sines and cosines of the unknown angles Ω, i, and ω, and for a generic orientation it can be much smaller than unity. The text acknowledges this but does not quantify the resulting spread. I recommend adding a short Monte Carlo over Ω (and over i where not directly measured) to show the distribution of the resulting limits, or at least a clear statement that the quoted numbers are favorable-orientation maximal-reach estimates rather than orientation-independent bounds.","section":"Sec. IV, Eqs. (24)–(26) and Table IV"},{"comment":"The observed ˙xp for PSR J0751+1807 and PSR J1738+0333 is used directly as an upper limit on Lorentz-violating contributions, with the argument that proper-motion contamination is not conspiratorial. The proper-motion contribution to ˙xp can be of the same order as the measured values for these nearby pulsars, so a partial cancellation could shift the derived limit by an O(1) factor without any fine-tuning. This does not invalidate the maximal-reach approach, but it is the weakest link in the chain. Please provide the proper-motion values used and/or perform an Ω-randomization for the ˙xp analysis analogous to what was done for the η, κ analysis in Ref. [21], so that the systematic effect is quantified rather than asserted.","section":"Sec. IV, paragraph on using measured ˙xp as an upper limit"}],"minor_comments":[{"comment":"The symbol m in the definition (a_eff)_μ = a_μ − m e_μ is not explicitly identified as the fermion mass in the surrounding text; a brief parenthetical would improve readability.","section":"Sec. II, Eq. (4)"},{"comment":"The entry 'Ne ∼ 0' for neutron stars could be clarified as 'negligible' to avoid implying an exactly zero electron content.","section":"Table I"},{"comment":"The statement 'the more relativistic the binary (namely, the larger nb)' is slightly imprecise because nb is the orbital frequency; the intended meaning is the shorter the orbital period, the larger the orbital velocity, and rephrasing would avoid confusion.","section":"Sec. IV, first paragraph"},{"comment":"The phrase 'Nature's conspiracy' is informal for a journal article; consider replacing it with a more neutral formulation such as 'fine-tuning of the unknown longitude of ascending node'.","section":"Sec. IV, discussion of proper motion"},{"comment":"The numerical factor 0.0005 in the electron terms is the ratio m_e/m_p; stating this explicitly would help readers quickly verify the composition-weighted coefficients.","section":"Eq. (17)"}],"recommendation":"minor_revision","confidential_remarks":"The paper is a solid phenomenological contribution appropriate for the special issue. The central derivation is sound and the limitations are acknowledged in the text; the requested changes are local clarifications and quantifications rather than corrections of a fundamental error. I support publication after minor revision."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThe paper is a straightforward, honest application of existing SME machinery to three small-eccentricity NS-WD pulsars. What is actually new: the first binary-pulsar limits on matter-gravity SME coefficients for neutrons, protons, and electrons. The best numbers, c^w_jk ~ 1e-11 from J1738's xp_dot, and (a^n_eff)_k ~ 1e-8 GeV from J0751's eta/kappa, are competitive and complementary to LLR and clock bounds. The theory is from Kostelecky-Tasson and Jennings et al.; the author's contribution is the small-eccentricity reduction to Eqs. (21)-(26) and the limit-setting. That is a legitimate new application, not a breakthrough.\n\nThe paper does several things well. The expansions are clean, and the e<1e-6 condition is firmly met. The author is explicit about the simplifications: O(1) orientation projection, proper-motion contamination in xp_dot, and the strong-field caveat for NSs. The published timing measurements are external inputs, so there is no circularity. Table IV is easy to read, and the provenance of each limit is clear.\n\nThe soft spots are the ones the author flags, and they are real. Treating projections as O(1) operators is an assumption, not a bound. For a single nonzero Cartesian component, the orbital-frame coefficient is a product of sines and cosines of unknown Euler angles; for some orientations it is much smaller than one. The 'no conspiracy' argument is a plausible prior, but it does not make the limits conservative. Likewise, the measured xp_dot of J1738 is small enough that proper motion could be comparable, so a partial cancellation could shift that limit by an O(1) factor. These caveats mean Table IV should be read as favorable-orientation reach, not worst-case bounds. The author says as much in Sec. IV, but the abstract and conclusions could state it more prominently.\n\nIs this fatal? No. The paper explicitly frames the results as maximal-reach limits, and that is what they are. The assumptions are standard for first constraints of this type, and the order-of-magnitude claims stand. Minor improvements: explicit error ranges on Table IV and a sentence in the abstract about the favorable-orientation assumption.\n\nWho gets value from this: SME phenomenologists and pulsar timers. It deserves a serious referee; I would send it out with a minor-revision recommendation, mainly to sharpen the caveat language.","headline":"Solid, honest maximal-reach limits on matter-gravity SME coefficients from three small-eccentricity pulsars; the O(1) orientation projection is the acknowledged weak point but does not undercut the order-of-magnitude claims.","tokens_in":12152,"tokens_out":3298,"would_cite":true,"duration_ms":31020,"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":"Binary pulsar timing can set the tightest gravitational limits on Lorentz-violating matter–gravity couplings, down to $10^{-11}$ for spatial neutron and proton coefficients.","keywords":["Lorentz violation","Standard-Model Extension","binary pulsar timing","matter-gravity couplings","neutron star","white dwarf","maximal-reach limits","small-eccentricity binaries"],"falsifier":"Measure the longitude of the ascending node for PSR J0751+1807 or PSR J1738+0333 and recompute the projections of $B_j$ and $A_{jl}$ onto the orbital frame; if the true projection factors turn out to be much smaller than order one, the maximal-reach limits in Table IV would weaken by more than a factor of a few, contradicting the paper's order-one projection assumption for that pulsar.","tokens_in":11107,"feed_emoji":"🔭","tokens_out":13613,"duration_ms":119302,"temperature":0.7,"pith_summary":"This paper claims that three small-eccentricity neutron-star–white-dwarf binaries can act as precision laboratories for Lorentz violation in the way ordinary matter couples to gravity. It uses pulsar-timing measurements of the projected semimajor axis $\\dot{x}_p$ and of the Laplace–Lagrange parameters $\\dot{\\eta}$ and $\\dot{\\kappa}$ to put upper limits on the matter–gravity coefficients of the Standard-Model Extension for neutrons, protons, and electrons. The tightest spatial limits reach $\\sim 10^{-11}$ for $c^n_{jk}$ and $c^p_{jk}$, and the effective neutron vector coefficient $(\\bar a^n_{\\mathrm{eff}})_k$ is bounded near $\\sim 10^{-8}$ GeV. Because neutron stars are strongly self-gravitating, the limits are conservative strong-field tests, complementary to lunar laser ranging and laboratory clock experiments. If the analysis is right, binary pulsars become one of the best gravitational probes of matter–gravity Lorentz violation.","feed_headline":"Pulsar binaries pin matter–gravity Lorentz violation to 1e-11","feed_subtitle":"Three neutron-star–white-dwarf pairs put the tightest gravitational bounds yet on matter–gravity SME coefficients.","key_machinery":"The central object is the point-particle action $\\mathcal S_u = \\int d\\lambda\\,[-m\\sqrt{-(g_{\\mu\\nu}+2\\bar c_{\\mu\\nu})u^\\mu u^\\nu} - (\\bar a_{\\mathrm{eff}})_\\mu u^\\mu]$, in which $\\bar c_{\\mu\\nu}$ and $(\\bar a_{\\mathrm{eff}})_\\mu$ are the species-dependent coefficient fields for Lorentz violation. For a composite body, these coefficients are replaced by particle-number-weighted sums over neutron, proton, and electron constituents, giving the composite $A_{jl}$ and $B_j$ used in the orbital equations. The mechanism that carries the argument is the small-eccentricity limit of the osculating-element equations, where $\\langle d x_p/dt\\rangle$ and the Laplace–Lagrange parameter derivatives $\\langle d\\eta/dt\\rangle$, $\\langle d\\kappa/dt\\rangle$ become linear in $A_{jl}$ and $B_j$. The 'maximal-reach' strategy then assumes one SME coefficient is nonzero at a time and treats the unknown projection of the orbit onto the sky as an order-one number, converting the absence of measured secular changes into upper limits.","core_discovery":"On its own terms, the paper's central discovery is that in the small-eccentricity limit the secular orbital changes of a relativistic binary reduce to simple, clean signals of Lorentz violation. The averaged derivative of the projected semimajor axis, $\\langle d x_p/dt\\rangle$, is controlled by a combination of the tensor components $A_{\\hat a\\hat c}\\cos\\omega - A_{\\hat b\\hat c}\\sin\\omega$, while the averaged derivatives of the eccentricity-vector components, $\\langle d\\eta/dt\\rangle$ and $\\langle d\\kappa/dt\\rangle$, are controlled by the vector components $B_{\\hat a}$ and $B_{\\hat b}$. These $A$ and $B$ objects are particle-number-weighted sums of the SME coefficients $c^w_{(jk)}$, $c^w_{(0j)}$, and $(\\bar a^w_{\\mathrm{eff}})_j$ over neutrons, protons, and electrons, with the electron's $c^w$ weight suppressed by the electron–proton mass ratio. Using measured and estimated time derivatives for PSRs J0348+0432, J0751+1807, and J1738+0333, the paper derives maximal-reach bounds: $c^{n,p}_{jk}\\lesssim 10^{-11}$, $c^{n,p}_{0k}\\lesssim 10^{-8}$, and $(\\bar a^n_{\\mathrm{eff}})_k\\lesssim 10^{-8}$ GeV, with PSR J1738+0333 giving the best $c^w_{jk}$ limits and PSR J0751+1807 the best $c^w_{0k}$ and $(\\bar a^w_{\\mathrm{eff}})_k$ limits.","pith_inferences":["If the longitude of the ascending node were measured for PSR J0751+1807 or PSR J1738+0333, the order-one projection assumption could be replaced by an explicit three-dimensional fit, which would either confirm the quoted limits or sharpen them into correlated bounds on several coefficients at once.","A dedicated simultaneous fit of $\\dot{x}_p$, $\\dot{\\eta}$, and $\\dot{\\kappa}$ with proper motion modeled would separate a genuine Lorentz-violating residual from the astrophysical contaminants that the current maximal-reach approach absorbs into upper limits.","Applying the same small-eccentricity equations to newly discovered short-period neutron-star–white-dwarf binaries should tighten the limits faster than the $T^{-1.5}$ scaling alone, because the signal grows with the orbital frequency.","The current numbers should be read as order-of-magnitude reach: if a future signal appears, the maximal-reach approach cannot identify which coefficient is responsible, and a full orbital-orientation fit would be required."],"forward_implications":["The spatial coefficients $c^n_{jk}$ and $c^p_{jk}$ are pinned below $\\sim 10^{-11}$, with PSR J1738+0333 providing the tightest limit through its accurate $\\dot{x}_p$ measurement.","The mixed-index coefficients $c^n_{0k}$ and $c^p_{0k}$ are bounded below $\\sim 10^{-8}$, and the neutron vector coefficient $(\\bar a^n_{\\mathrm{eff}})_k$ below $\\sim 10^{-8}$ GeV, both best constrained by PSR J0751+1807 via $\\dot{\\eta}$ and $\\dot{\\kappa}$.","Because neutron stars are strongly self-gravitating, the limits are conservative strong-field versions of the SME coefficients; any strong-field enhancement would tighten them rather than weaken them.","Longer timing baselines improve the limits roughly as $T_{\\mathrm{obs}}^{-3/2}$ even without instrumental upgrades, and new radio telescopes should speed up the improvement.","The electron $c^w_{jk}$ and $c^w_{0k}$ limits are suppressed by the electron–proton mass ratio, while the electron $(\\bar a^w_{\\mathrm{eff}})_k$ limits are not, so both charged-lepton sectors receive independent bounds."],"supporting_citations":[{"why":"Supplies the point-particle action and the composite-coefficient prescription for Lorentz-violating matter–gravity couplings.","marker":"[11]"},{"why":"Supplies the osculating-element equations for a binary under Lorentz-violating matter–gravity couplings that the paper simplifies to small eccentricity.","marker":"[31]"},{"why":"Establishes the binary-pulsar SME methodology and the $\\sqrt{12}\\sigma/T_{\\mathrm{obs}}$ estimate used for unmeasured derivatives.","marker":"[21]"},{"why":"Provides PSR J0348+0432's masses and orbital parameters entering Table III and Table IV.","marker":"[32]"},{"why":"Provides PSR J0751+1807's timing parameters, including measured $\\dot{x}_p$, $\\eta$, and $\\kappa$, which yield the best limits on $c^w_{0k}$ and $(\\bar a^w_{\\mathrm{eff}})_k$.","marker":"[33]"},{"why":"Provides PSR J1738+0333's timing parameters and its precise $\\dot{x}_p$ measurement, which yields the best limits on $c^w_{jk}$.","marker":"[34]"},{"why":"Defines the maximal-reach approach of assuming one coefficient nonzero at a time, used for the bounds in Table IV.","marker":"[18]"},{"why":"Supplies the proper-motion contribution to $\\dot{x}_p$, the systematic effect the paper treats as an upper limit.","marker":"[40]"}],"fun_headline_variants":["Pulsars bound matter–gravity Lorentz violation to 1e-11","Matter–gravity Lorentz violation pinned by pulsar binaries","Small-eccentricity pulsars bind Lorentz violation in matter","Pulsar timing constrains matter–gravity Lorentz violation"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing assumption is that the unknown orientation of each orbit projects the Lorentz-violating coefficients onto the measured derivatives with order-one factors, and that the pulsar's proper motion does not systematically cancel the signal; if either fails, some quoted limits could weaken by more than a factor of a few.","fun_headline_variants_meta":{"raw":{"variants":["Pulsars bound matter–gravity Lorentz violation to 1e-11","Matter–gravity Lorentz violation pinned by pulsar binaries","Small-eccentricity pulsars bind Lorentz violation in matter","Pulsar timing constrains matter–gravity Lorentz violation"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00102,"raw_usage":{"total_tokens":4330,"prompt_tokens":995,"completion_tokens":3335,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":611,"completion_tokens_details":{"reasoning_tokens":3263}},"tokens_in":611,"tokens_out":3335,"duration_ms":24623,"temperature":1.0,"reasoning_tokens":3263,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T10:55:39.569761+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the longitude of the ascending node for PSR J0751+1807 or PSR J1738+0333 and recompute the projections of $B_j$ and $A_{jl}$ onto the orbital frame; if the true projection factors turn out to be much smaller than order one, the maximal-reach limits in Table IV would weaken by more than a factor of a few, contradicting the paper's order-one projection assumption for that pulsar.","supporting_citations":[{"cited_title":"New limits on the violation of local position invariance of gravity","cited_arxiv_id":"1307.2637","evidence_quote":"Provides PSR J0348+0432's masses and orbital parameters entering Table III and Table IV."},{"cited_title":"Maximal Tests in Minimal Gravity","cited_arxiv_id":"1907.08106","evidence_quote":"Defines the maximal-reach approach of assuming one coefficient nonzero at a time, used for the bounds in Table IV."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the proper-motion contribution to $\\dot{x}_p$, the systematic effect the paper treats as an upper limit."}],"review_version":1}