{"id":"8cf18670-f14e-4452-ad56-7df81ecedfb8","arxiv_id":"2511.05657","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"Lambda-antilambda oscillations are already bounded far more strongly by Super-Kamiokande's pp→K+K+ limit (δmΛ ≲ 10^-32 GeV) than by BESIII's direct search, making dedicated oscillation searches non-competitive in the considered models.","lead":"The paper studies Λ–Λ̄ oscillations, baryon-number-violating transitions driven by six-quark operators with strangeness, and derives indirect bounds from existing neutron and dinucleon experiments. It concludes that the pp→K+K+ process already constrains these oscillations far more strongly than dedicated BESIII searches, redirecting experimental effort.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"pp→K+K+ bound on δmΛ assumes the two-hadron matrix element equals the single-baryon one; a counterterm or lattice ratio could change the 14-order claim.","rationale":"The reader's weakest_assumption correctly points to hadronic matrix element uncertainties, but the more load-bearing gap is the relation between the pp→K+K+ amplitude and δmΛ itself. The paper's Eq. (3.27) is a pole-diagram estimate that omits the full propagator numerator and ignores potential counterterms; Eq. (3.32) then equates a two-hadron matrix element to δmΛ without justification. If this ratio is not O(1), the quantitative hierarchy — the paper's central assertion — shifts. However, the concern does not fully undermine the qualitative ordering: even a factor of 100 error in the amplitude would still leave pp→K+K+ much stronger than BESIII, and the paper's own order-of-magnitude caveats are visible. Thus the appropriate verdict remains CONDITIONAL, not REJECT or ACCEPT. The proposed lattice ratio is the specific check that would settle whether the fourteen-order superiority is real or an artifact of the hadronization assumption.","tokens_in":28270,"tokens_out":21554,"duration_ms":187054,"concrete_test":"Compute on the lattice the ratio R = ⟨KK|O_{(uds)^2}|NN⟩ / ⟨Λ|O_{(uds)^2}|Λ⟩ with the same operator O and physical quark masses, using the framework of Refs. [38,39]. If R differs from 1 by more than an order of magnitude, the dinucleon bound δmΛ ≲ 10^-32 GeV must be rescaled by R, and the claimed fourteen-order gap over BESIII would shrink accordingly. Until such a calculation exists, treat Eq. (3.30) as an assumption, not a bound.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim that Super-K's pp→K+K+ limit gives δmΛ ≲ 10^-32 GeV (Eq. 3.30) — fourteen orders below BESIII — is derived from a tree-level BχPT diagram in which an off-shell Λ propagator connects two NN→ΛK vertices, with δmΛ inserted at the Λ-Λbar transition. This calculation (Eq. 3.27) drops the /q+mΛ numerator of the fermion propagator and replaces the full t+u amplitude by 2× the t-channel evaluated at one kinematic point. More importantly, it implicitly assumes that the short-distance NN→KK amplitude is saturated by the δmΛ pole diagram. In EFT, δmΛ is a low-energy constant for the single-baryon transition; the two-baryon operator ⟨KK|O|NN⟩ has independent counterterms not fixed by δmΛ. The alternative nuclear estimate (Eqs. 3.31–3.35) makes this explicit by setting C_i⟨KK|O|pp⟩ ≡ δmΛ — i.e., equating a two-meson matrix element to the Λ-Λbar mass parameter. If ⟨KK|O|NN⟩ differs from ⟨Λ|O|Λ⟩ substantially (there is no lattice calculation for either), the bound shifts by exactly that ratio, and the claimed 'strongest constraint' could be weakened by orders of magnitude. The paper itself labels the nuclear treatment 'an order-of-magnitude estimate', but the headline conclusion quotes no uncertainty.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper extends the effective-field-theory treatment of |ΔB|=2 baryon-number violation from neutron–antineutron oscillations to strangeness-violating Λ–Λbar oscillations. It classifies dimension-9 LEFT and SMEFT operators relevant for the (uds)^2 sector, enumerates scalar-mediated UV completions (trilinear, quartic, and fermionic topologies), and derives indirect bounds on the mass-mixing parameter δmΛ from n–nbar searches, from Super-Kamiokande pp→K+K+ searches, and from BESIII. In a simplified model with the scalars Sbar1 and Ω4, the paper finds that n–nbar is induced only at two loops while Λ–Λbar appears at tree level, and that pp→K+K+ gives the strongest lower bound on the effective BNV scale, Λ_BNV ≳ 300 TeV. The central claim is that current Λ–Λbar oscillation searches are far from competitive with dinucleon-decay bounds.","tokens_in":28708,"tokens_out":19981,"duration_ms":179012,"significance":"If the hadronic matrix-element assumptions are controlled, this is a useful and original phenomenological guide. The operator classification goes beyond prior n–nbar studies, with machine-checked enumeration via Sym2Int and a systematic survey of UV completions; the appendices on fermionic and quartic completions are valuable. The simplified model is explicit and falsifiable: it predicts that pp→K+K+ at Super-K, not BESIII, is the best probe of strangeness-violating |ΔB|=2 physics. The paper also usefully identifies the need for lattice calculations of ⟨Λ|O|Λ⟩ and two-meson nuclear matrix elements. I find no circularity: the bounds are derived from external experimental limits, with the simplified-model bounds being legitimate parameter translations. The main weakness is that the headline numerical hierarchy relies on hadronic matrix-element identifications that are assumed rather than derived, and the paper gives no uncertainty budget.","major_comments":[{"comment":"The headline bound δmΛ ≲ 1.8×10^-32 GeV from pp→K+K+ is not a model-independent constraint. The BχPT amplitude in Eq. (3.27) inserts the single-baryon mass-mixing operator δmΛ into a long-distance pole diagram. In the hadronic EFT, the same six-quark operator also matches onto local NN→KK operators whose coefficients are not fixed by δmΛ. Eq. (3.32) makes the implicit identification explicit: C_i⟨K+K+|O_i|pp⟩ ≡ δmΛ. If the two-hadron matrix element differs from ⟨Λ|O|Λ⟩, the quoted bound shifts by exactly that ratio. Since the paper presents this as the central result and uses it in Sec. 4 and Table 4, please provide an estimate or a conservative range for R = ⟨K+K+|O|pp⟩/⟨Λ|O|Λ⟩, and show how the claimed fourteen-order separation from BESIII changes under plausible variations. The label 'order-of-magnitude estimate' is not sufficient when the numerical separation is the main quantitative","section":"Sec. 3.2, Eqs. (3.27)–(3.30)"},{"comment":"The alternative nuclear estimate is not internally consistent as written. In Eq. (3.32), M_A = δmΛ ρpp(0) with ρpp(0) = κ|ρ_p(0)|^2: with densities in fm^-6 and δmΛ in GeV, this amplitude has dimension mass^7. The phase-space integral is quoted in Eq. (3.34) as 4×10^-7 GeV^-1, whereas the standard three-body phase space for this decay has dimension mass^2, and no fm↔GeV conversion factors are displayed in the rate formula. Consequently Eq. (3.35) cannot be reproduced from Eq. (3.31) by a dimensionally consistent calculation. This is load-bearing because Table 5 presents the nuclear estimate as independent confirmation of the chiral estimate. Please rewrite with explicit dimensions and numerical conversions, or clearly label Eq. (3.35) as a parametric dimensional estimate.","section":"Sec. 3.2, Eqs. (3.31)–(3.35)"},{"comment":"The numerical bounds are quoted without any uncertainty budget. The unknown Λ matrix element is replaced by the lattice neutron value in Eq. (2.13), with only a symbolic O(m_u,d/m_s) correction; the one-loop integrals in Eqs. (3.13) and (3.21) are regulated with an ad hoc 5 GeV cutoff; and the chiral couplings a,b in Eq. (3.8) carry fit errors that are not propagated. Since the paper's conclusion is quantitative—fourteen orders of magnitude—please provide a conservative error range for each bound, or at least show how the bounds in Tables 4 and 5 depend on the key inputs. This pass should also fix the dimension/step inconsistency in Eq. (3.22): the quantity δm_n^K is written as 9.8×10^-16 GeV, but the subsequent bound δm_Λ^K ≲ 1.4×10^-18 GeV requires that this be interpreted as the dimensionless coefficient 9.8×10^-16 times δmΛ.","section":"Secs. 2.2, 3.1 and Tables 4–5"}],"minor_comments":[{"comment":"The amplitude in Eq. (3.27) drops the /q+mΛ numerator of the Λ propagator. This may be an O(1) approximation for nonrelativistic kinematics, but it should be stated. The factor of 2 multiplying 2mp for the u-channel contribution is also not explained in the text.","section":"Eq. (3.29) and Fig. 5"},{"comment":"The Λ operator basis is not displayed; the reader is told that the number of independent operators is 52 (LEFT) and 20 (SMEFT), with the derivation left to a forthcoming publication. A representative subset or an ancillary file would make the counting reproducible and easier to check.","section":"Sec. 2.2"},{"comment":"The figure caption does not identify the colored regions or the axes beyond 'parameter space'. Please add a legend and define the plotted quantity (e.g., the effective scale Λ_BNV versus scalar mass).","section":"Figure 9"},{"comment":"The hard-cutoff dependence of the one-loop bounds is not discussed. Varying Λχ between ~1 and 5 GeV changes the coefficients in Eqs. (3.14) and (3.22), and the resulting δmΛ bounds are close enough to the BESIII limit that the comparison should include this variation.","section":"Sec. 3.1"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is within the scope of the journal and contains useful new classification and model-building material. The main risk is the hadronic matrix-element identification behind the pp→K+K+ bound; this is a physics assumption, not a calculation error, but it is central to the paper's strongest claim. The dimensional inconsistency in the nuclear estimate strengthens the need for a careful revision. I would recommend sending the revised version to a referee with expertise in chiral perturbation theory and nuclear matrix elements. The paper's own caveats and the 'forthcoming publication' for the operator basis should be addressed during revision."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Here's the short version: this is a genuinely useful paper for the BNV subfield. It gives the first systematic EFT/chiral treatment of Λ–Λ̄ oscillations, classifies the (uds)^2 dimension-9 operators (52 LEFT / 20 SMEFT), adds a previously neglected class of trilinear scalar UV completions, and—most importantly—shows that Super-K's pp→K+K+ limit already constrains δmΛ far more strongly than BESIII's direct oscillation search. The qualitative hierarchy is almost certainly right; the fourteen-order gap to BESIII is enormous.\n\nWhat's new and good: the operator counting, the µX1X2X3 completions, and the simplified model where Λ–Λ̄ is tree-level and n–n̄ is two-loop. The chiral translation of n–n̄, BESIII, and pp→K+K+ limits into common bounds on δmΛ is exactly the kind of apples-to-apples comparison the field needs. The paper is careful about the proton-decay constraints on scalar LQs and about the GIM-like suppression in the two-loop n–n̄ amplitude.\n\nSoft spots, in order of importance. First, the pp→K+K+ bound rests on the assumption that the two-nucleon→two-kaon amplitude is saturated by the δmΛ pole diagram. That is not forced by EFT: the short-distance operator has independent matrix elements between two-nucleon and two-kaon states. The alternative nuclear estimate (Eq. 3.32) makes this explicit by setting C_i⟨KK|O|pp⟩ ≡ δmΛ, and if that matrix element differs from the single-baryon one, the bound shifts by exactly that ratio. It would take a large suppression to matter given the fourteen-order gap, so I would not bet against the ordering, but the number 10^-32 GeV has no uncertainty attached and is an order-of-magnitude estimate by the authors' own admission. Second, the Λ matrix element is just set equal to the neutron lattice value (Eq. 2.13); there is no lattice input for the strange system. Third, the operator basis is asserted via Sym2Int but not actually listed—the 52/20 counts are not checkable from the paper. Fourth, the chiral loops are regulated with a hard cutoff at 5 GeV, which is above the usual BχPT validity range, and Eq. (3.22) quotes a dimensionless coefficient with a GeV unit.\n\nI am not worried about circularity or data-manipulation; the bounds are derived from external limits and standard LECs. The paper would benefit from a referee insisting that the authors (1) present the explicit operator lists or deposit them, (2) attach an honest error bar to the δmΛ bounds, and (3) fix the typo and explain the matching assumption more carefully.\n\nVerdict: send to peer review. It's a solid subfield contribution that will be cited and should be on the radar of anyone working on baryon-number violation or planning dinucleon searches. I'd bring it to a reading group.","headline":"Useful, carefully scoped EFT/chiral paper: the pp→K+K+ bound on δmΛ is the strongest current probe and the ordering is likely robust, but the headline number is an order-of-magnitude estimate with unquantified hadronic uncertainties.","tokens_in":29225,"tokens_out":6517,"would_cite":true,"duration_ms":48812,"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":"Dinucleon decay, not hyperon oscillations, sets the tightest bound on strangeness-changing baryon-number violation.","keywords":["baryon number violation","Λ–Λ̄ oscillations","neutron–antineutron oscillations","dinucleon decay","six-quark operators","SMEFT","chiral perturbation theory","dimension-9 operators"],"falsifier":"A lattice QCD computation of the Λ→Λ̄ six-quark matrix element that found it to be many orders of magnitude smaller than the neutron value would invalidate the claim that pp→K+K+ constrains δm_Λ at the 10^-32 GeV level; conversely, an improved bound or a first positive signal in pp→K+K+ would directly test the prediction.","tokens_in":28163,"feed_emoji":"⚛️","tokens_out":4563,"duration_ms":43227,"temperature":0.7,"pith_summary":"The paper extends baryon-number-violating (ΔB=2) physics beyond neutron–antineutron oscillations to the strange baryon Λ, showing that Λ–Λ̄ mixing can be generated at tree level by six-quark operators of the form (uds)^2 that are independent of the (udd)^2 operators behind n–n̄ oscillations. Using the Standard Model Effective Field Theory and scalar-mediated UV completions, the authors derive indirect bounds on the effective Λ mass-mixing parameter δm_Λ from existing limits on n–n̄ oscillations and from the dinucleon decay p p → K+ K+. They find that the dinucleon channel dominates by far, constraining δm_Λ below about 10^-32 GeV — roughly fourteen orders of magnitude stronger than the direct bound from J/ψ→Λ Λ̄ oscillations. The paper thereby establishes pp→K+K+ searches as the most powerful current probe of strangeness-violating ΔB=2 transitions, sensitive to new-physics scales up to ~300 TeV, and argues that Λ–Λ̄ oscillation experiments would need roughly four orders of magnitude improvement to become competitive.","feed_headline":"pp→K+K+ bound beats Λ–Λ̄ oscillation limits","feed_subtitle":"Strangeness-changing ΔB=2 effects probe new physics up to ~300 TeV, far beyond direct hyperon limits.","key_machinery":"The central object is the effective mass mixing δm_Λ, the off-diagonal element in the two-state Hamiltonian for Λ and its antiparticle, generated by six-quark operators (uds)^2 with Wilson coefficients suppressed by the fifth power of a new-physics scale Λ_BNV. The argument is carried by chiral effective theory, which connects δm_Λ to observables: a tree-level diagram converts δm_Λ into n–n̄ mixing via weak vertices, pion and kaon loops give subleading contributions, and a tree-level t-channel diagram produces the amplitude for p p→K+ K+. That amplitude, combined with a nuclear-density estimate of the intranuclear rate, is what turns the dinucleon lifetime limit into the sharp bound on δm_Λ.","core_discovery":"The paper's central claim is that strangeness-violating |ΔB|=2 baryon-number violation is currently best probed by the dinucleon decay p p→K+ K+, not by Λ–Λ̄ oscillations. Starting from the Standard Model Effective Field Theory, the authors identify 52 six-quark operators of the (uds)^2 type (up from 14 for (udd)^2), classify their tree-level scalar-mediated UV completions, and identify models in which Λ–Λ̄ mixing is generated at tree level while n–n̄ mixing appears only at two loops. Using chiral effective theory to bridge quark operators and baryon observables, they derive indirect bounds on δm_Λ from the measured n–n̄ oscillation limit and from the water-Cherenkov detector's limit on p p→","pith_inferences":["If a lattice QCD computation of the Λ six-quark matrix element deviates from the neutron value by more than an order of magnitude, the quantitative hierarchy among bounds would shift, though the qualitative ordering would likely survive because the dinucleon bound is so much stronger.","The same chiral machinery could be applied to Ξ− or Ω− hyperon oscillations, or to ΔB=2 processes with charm quarks, extending the operator catalogue to (ucs)^2 or mixed-flavour combinations.","A dedicated re-analysis of existing water-Cherenkov data on pp→K+K+ using modern nuclear matrix elements would be the cheapest experimental test of the paper's central claim.","The two-loop, GIM-like suppression of n–n̄ relative to Λ–Λ̄ in the simplified model suggests a generic way to hide n–n̄ oscillations while leaving Λ–Λ̄ or dinucleon channels observable."],"forward_implications":["If correct, pp→K+K+ searches at water-Cherenkov detectors currently probe |ΔB|=2 at effective scales up to about 300 TeV, far beyond direct collider reach.","Models with (uds)^2 operators can yield Λ–Λ̄ oscillations at tree level while n–n̄ mixing is loop-suppressed, so the two channels probe genuinely different operator directions.","Improving the pp→K+K+ limit by an order of magnitude would push the new-physics scale into the PeV region.","Λ–Λ̄ oscillation experiments would need to improve δm_Λ sensitivity by roughly four orders of magnitude to compete with collider mass bounds, which appears infeasible in the foreseeable future.","The classification of 52 LEFT operators and 20 SMEFT operators provides a working basis for future lattice computations of the hadronic matrix elements."],"fun_headline_variants":["pp→K+K+ beats Lambda-LambdaBar for ΔB=2","Λ–Λ̄ oscillations? Try pp→K+K+ decay","Strangeness flips baryon number: (uds)^2 to 300 TeV","Baryon violation with strangeness: pp→K+K+ wins"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The quantitative bound on δm_Λ assumes the six-quark matrix element for the Λ equals the lattice neutron value up to O(1) factors, and that logarithmically divergent one-loop chiral integrals are regulated with a hard cutoff at 5 GeV; if those hadronic estimates are severely wrong, the numerical hierarchy among bounds shifts, though the ordering probably survives.","fun_headline_variants_meta":{"raw":{"variants":["pp→K+K+ beats Lambda-LambdaBar for ΔB=2","Λ–Λ̄ oscillations? Try pp→K+K+ decay","Strangeness flips baryon number: (uds)^2 to 300 TeV","Baryon violation with strangeness: pp→K+K+ wins"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001353,"raw_usage":{"total_tokens":5379,"prompt_tokens":845,"completion_tokens":4534,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":589,"completion_tokens_details":{"reasoning_tokens":4451}},"tokens_in":589,"tokens_out":4534,"duration_ms":29961,"temperature":1.0,"reasoning_tokens":4451,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-03T23:26:33.807490+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A lattice QCD computation of the Λ→Λ̄ six-quark matrix element that found it to be many orders of magnitude smaller than the neutron value would invalidate the claim that pp→K+K+ constrains δm_Λ at the 10^-32 GeV level; conversely, an improved bound or a first positive signal in pp→K+K+ would directly test the prediction.","supporting_citations":[],"review_version":1}