{"id":"c833cb17-530a-4b88-81dc-90782f771b16","arxiv_id":"1908.11155","paper_version":2,"verdict":"ACCEPT","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"low","formal_verification":"none","parameter_count":1,"one_line_summary":"Scalar leptoquark couplings are constrained by a comprehensive set of upper bounds from tau, muon, electron, and kaon decay data, including standard model interference in rare K to pi l+ l- modes.","lead":"This paper works out how tightly experiments on rare lepton and kaon decays restrict hypothetical particles called scalar leptoquarks, which could link quarks and leptons. It provides a comprehensive map of the strongest limits and shows where future measurements can improve them.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Table 9's K_L→π0ℓ+ℓ− ranges are chirality/Wilson-sign dependent (Eq. 70) but the grouped column (˜R2, ˜S1, 4×S3) holds three distinct (C, s_Y) classes; the quoted ranges, incl. the '(for ˜S1)' entry, match only one class, so some bounds are misassigned.","rationale":"This is a careful, transparent constraints catalog. Before forming the critique, I verified the internal normalization of several bounds: from Eq. (57) and the stated f_K, m_K, τ_L inputs, the K_L→e+e− and K_L→μ+μ− entries in Table 8 (2.0×10^-3 and 4.7×10^-5) follow from allowing the leptoquark amplitude to saturate the quoted experimental uncertainty, and the K+→π+νν̄ structure in Eq. (80) is dimensionally consistent with Eqs. (75)–(76). The one-at-a-time extraction flagged by the reader is explicitly stated in Section 4.1 and is standard for a first catalog, so I do not treat it as an error. The load-bearing concern is narrower and sits in the paper's most intricate new ingredient, the SM–NP interference treatment for K_L→π0ℓ+ℓ−. Eq. (70)'s interference terms depend on the operator chirality (s_Y) and on the overall Wilson-coefficient sign, and the four leptoquark types that contribute to these modes fall into three different (C, s_Y) classes. Table 9 quotes only two ranges per row (one for R2, one for the grouped ˜R2/˜S1/S3 column), and the values are mutually consistent only if the grouped range is the S3-class interval, contradicting its '(for ˜S1)' annotation. Since these rows are among the strongest constraints of the whole catalog (10^-4 level on imaginary parts of coupling products), a misassignment would change the central catalog entries, not merely their interpretation. If the check I propose confirms the concern, the appropriate remedy is to split the grouped column into per-chirality ranges, or to state explicitly which chirality each quoted range corresponds to; the paper's other results are unaffected. If the check shows that the published table already contains correct per-leptoquark assignments, the ACCEPT verdict stands unchanged. Given the non-trivial probability of a genuine misassignment in these two rows, I recommend CONDITIONAL rather than UNCHANGED.","tokens_in":125,"tokens_out":49191,"duration_ms":626206,"concrete_test":"Recompute the 90% C.L. allowed ranges for Im(y11 y*12) and Im(y21 y*22) (R2) and for Im(y*11 y21) and Im(y*12 y22) (˜R2, ˜S1, S3^{4/3}) in K_L→π0e+e− and K_L→π0μ+μ− by evaluating Eq. (70) separately for each (C, s_Y) class: R2 (−,+), ˜S1 (+,+), S3^{4/3} (+,−), ˜R2 (−,−). Then compare the four intervals per mode with the two quoted intervals in Table 9. If the grouped interval equals the S3-class interval and the '(for ˜S1)' annotation is not explained, Table 9 needs correction; if each leptoquark's interval matches its assigned quoted entry, the concern is resolved.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The catalog's strongest kaon constraints include the K_L→π0ℓ+ℓ− rows of Table 9, extracted from the SM–NP interference distribution of Eq. (70). The interference terms carry two operator-dependent signs: the lepton-chirality sign s_Y, explicit in Eq. (70), and the overall sign C of the Wilson coefficient g^{XY}_{V,d} (C>0 for ˜S1 and S3^{4/3}; C<0 for R2 and ˜R2). The three leptoquark types in Table 9's grouped column (˜R2, ˜S1, 4×S3) span three distinct interference classes: (C>0, s_Y=+1) for ˜S1, (C>0, s_Y=−1) for S3^{4/3}, and (C<0, s_Y=−1) for ˜R2. Each class gives a different allowed range for Im(coupling), because the interference drive is C(v0 + s_Y·a0), with v0 and a0 the effective vector and axial direct-CPV SM amplitudes. A consistency read of the two quoted ranges per row: for the μμ mode, the asymmetries (−6.5 vs +5.1) and (−5.8 vs +5.7), in units of 10^-4, fix C(v0+s_Y·a0) of opposite signs, so the grouped entry matches the (C>0, s_Y=−1) class, i.e., the S3 interval. But the table annotates this entry '(for ˜S1)', which is right-handed and would need a range near the mirror of the R2 entry; ˜R2 would need a third interval. As rendered, Table 9 misassigns semileptonic-kaon bounds for at least two of the three grouped leptoquark types, and the Section 5.2 claim that SM–NP interference is 'fully taken into account' is not reflected in the table structure. These rows are among the strongest catalog constraints (≈10^-4 on Im of coupling products), so the affected entries are load-bearing. If the published table contains per-chirality sub-rows lost in this rendering, the concern dissolves; the check below settles it.","agreement_with_reader":"disagree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript constructs the low-energy effective Lagrangian induced by tree-level exchange of the five scalar leptoquarks, lists the resulting Wilson coefficients, runs them to hadronic scales with leading-log QCD, and uses existing upper bounds on rare lepton and kaon processes to derive a catalog of 90% C.L. limits on products of leptoquark Yukawa couplings. It also derives SM differential distributions for K→πℓ+ℓ− and K_L→π0ℓ+ℓ−, including the interference of SM and leptoquark amplitudes, and applies them to extract the bounds in Table 9. The numerical results are summarized in Tables 1-11 at a reference leptoquark mass of 1 TeV.","tokens_in":37424,"tokens_out":14873,"duration_ms":165035,"significance":"If correct, this is a useful model-independent reference for scalar leptoquark searches: the operator basis is complete, the RG treatment is standard, and all numerical inputs are tabulated, making the constraints reproducible. The treatment of SM-NP interference in K_L→π0ℓ+ℓ− is a genuine improvement over earlier kaon analyses. The main caveat is the one-coupling-at-a-time assumption, which the authors explicitly acknowledge but which should be kept in mind when using individual entries of the catalog.","major_comments":[{"comment":"The K_L→π0e+e− and K_L→π0μ+μ− rows in Table 9 are not correctly rendered for the grouped column '˜R2, ˜S1, 4×S3'. Equation (70) contains linear SM-NP interference terms whose signs depend on both the lepton-chirality factor sY (sR=+1, sL=−1) and the overall sign C of the Wilson coefficient g^{XY}_{V,d}. The three types in that grouped column sit in three different classes: ˜S1 is (C>0, sY=+1), S3^{4/3} is (C>0, sY=−1), and ˜R2 is (C<0, sY=−1). The two printed ranges, one labelled '(for ˜S1)' and one unlabelled, cannot cover all three classes: for example, flipping C while keeping sY fixed mirrors the interval through the V0(z) interference term, so the ˜R2 range should differ from the S3 range. As published, the table either misassigns or omits bounds for at least two of the three leptoquark types, and the Section 5.2 statement that the interference is 'fully taken into account' is not reflected in the table structure. Please recompute and present separate ranges for each (C, sY) class, or state explicitly which class each printed range belongs to.","section":"§5.2, Table 9"}],"minor_comments":[{"comment":"The column heading '4×S3' is used without a definition; please state explicitly that it refers to the factor 4 appearing in the relevant S3 operator combination, so that readers do not confuse it with the number of leptoquark fields.","section":"§4.1, Tables 1 and 2"},{"comment":"The meaning of the parenthetical second line of numerical coefficients and the superscripts on a^{1,2,3}_1 and a^{1,2,3}_2 is not explained in the text; please define the convention for the two leptoquark charge states and for the bracketed entries.","section":"§4.2.4, Table 6"},{"comment":"Calling the procedure of neglecting the SM contribution in order to derive the K+→π+ℓ+ℓ− bounds 'conservative' is misleading: since the measured rates are dominated by the SM amplitude, a full treatment that allows destructive interference would generally permit larger leptoquark couplings, so the quoted limits may be stronger rather than weaker than a complete fit. Please reword or justify the statement.","section":"§5.2, Eqs. (62)-(63)"},{"comment":"The statement that the limits assume all other contributions to be absent is easy to miss because the same single-operator assumption underlies most of Tables 1-11; I suggest restating this caveat prominently in the Summary so that the catalog is not over-interpreted as simultaneous bounds on all couplings.","section":"§4.1 and Summary"}],"recommendation":"major_revision","confidential_remarks":"The paper is well within the journal's scope and the overall derivation is sound. The main issue is the Table 9 K_L→π0ℓ+ℓ− ranges for the grouped leptoquark column, which affects some of the strongest kaon constraints and needs to be corrected before publication. Once that is addressed, I would be willing to accept the manuscript."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"You should know two things about this paper. First, it is a genuinely useful catalog: five scalar leptoquarks, matched to a four-fermion effective Lagrangian, run down to low energies, and translated into 90% C.L. bounds from tau, muon, electron, and kaon data. The genuinely new piece is the differential treatment of K→πℓ+ℓ− and K_L→π0ℓ+ℓ− that keeps the SM–NP interference terms, including the s_Y chirality factor of Eq. (70). The numerical inputs are tabulated, the assumptions are stated, and the formulas are explicit enough that a model builder can reproduce the limits. The one-coupling-at-a-time extraction is acknowledged in Section 4.1, and for a catalog of this type that is an acceptable limitation, not a hidden fit.\n\nThe bigger soft spot is Table 9. The K_L→π0ℓ+ℓ− bounds are chirality- and Wilson-sign-dependent because the interference drive in Eq. (70) is C(v0 + s_Y·a0). The grouped column \"R̃2, S̃1, 4×S3\" lumps together three distinct classes: S̃1 has (C>0, s_Y=+1), S3^{4/3} has (C>0, s_Y=−1), and R̃2 has (C<0, s_Y=−1). The rendered table gives two ranges per mode, labels one \"(for S̃1)\", and repeats the same coupling expression Im(y*_{1m} y_{2m}). That expression is actually R̃2's combination, not S̃1's, and the quoted ranges appear to match the S3 class rather than S̃1. In other words, at least two of the three grouped leptoquark types are assigned the wrong interval. I can't fully verify the numerics without rerunning the fit, but the internal inconsistency is visible from the text alone. If the published table contains per-chirality sub-rows lost in this extraction, then the concern dissolves; otherwise those entries are misassigned and need correction.\n\nThis issue does not sink the catalog. Most of the constraints in Tables 1–8 and 10–11 are absolute values or loop-dominated bounds where the sign convention does not matter. But the K_L→π0ℓ+ℓ− rows are among the strongest kaon constraints, so a model builder using those rows for S̃1 or R̃2 would get the wrong allowed region. The g−2 \"explanation\" ranges in Eqs. (29)–(30) are explicitly fits, not predictions, and the KOTO discussion is appropriately cautious.\n\nWho gets value: anyone filtering scalar leptoquark scenarios against low-energy flavor data, and kaon phenomenologists interested in the SM–NP interference formalism. It deserves a serious referee, and the misassignment, if real, is a minor-revision fix rather than a deep flaw. I would accept after the authors clarify or correct Table 9.","headline":"A solid, transparent constraint catalog for scalar leptoquarks with a genuinely new SM–NP interference treatment for K→πℓ+ℓ−, but Table 9's chirality-dependent ranges look misassigned in the rendered version.","tokens_in":38002,"tokens_out":13349,"would_cite":true,"duration_ms":116646,"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":"Tree-level scalar leptoquark exchange generates a definite set of four-fermion operators, and current rare-decay data convert that into 90% C.L.","keywords":["scalar leptoquarks","effective four-fermion Lagrangian","rare kaon decays","lepton flavour violation","muon-electron conversion","Yukawa couplings","Wilson coefficients","low-energy constraints"],"falsifier":"Measure the $q^2$ distribution of $K^+\\to\\pi^+\\mu^+\\mu^-$ with the current form-factor inputs; the paper predicts a definite interference term between the Standard Model vector amplitude and the leptoquark contribution, so a shape that excludes that interference would rule out the scalar-leptoquark explanation. Alternatively, a confirmed $K_L\\to\\pi^0\\nu\\bar\\nu$ rate above about $7.8\\times10^{-10}$ would falsify the lepton-flavour-conserving scalar-leptoquark picture and force the lepton-flavour-violating neutrino channel.","tokens_in":36765,"feed_emoji":"⚛️","tokens_out":11940,"duration_ms":108920,"temperature":0.7,"pith_summary":"The paper tries to establish that rare decays of leptons and kaons already put sharp, model-independent limits on scalar leptoquarks—hypothetical bosons that can convert a quark into a lepton. It constructs the full low-energy four-fermion Lagrangian produced by tree-level scalar leptoquark exchange, writes the Wilson coefficients as products of Yukawa couplings, and converts current experimental upper bounds into a catalog of 90% C.L. limits on those products. The strongest constraints come from $K_L\\to\\mu^+\\mu^-$, $K_L\\to\\pi^0\\ell^+\\ell^-$, $K_L\\to\\pi^0\\nu\\bar\\nu$, $\\mu\\to e\\gamma$, and $\\mu$-$e$ conversion in gold. If the catalog is right, it tells model builders which leptoquark couplings remain viable at the TeV scale and which future measurements will tighten the limits first.","feed_headline":"Kaon decays pin down scalar leptoquark couplings","feed_subtitle":"A new catalog of 90% C.L. limits shows rare decays squeeze the Yukawa couplings of TeV-scale leptoquarks.","key_machinery":"The central object is the low-energy effective Lagrangian generated by tree-level scalar leptoquark exchange. Scalar leptoquarks are hypothetical bosons carrying both quark and lepton flavour; after integrating them out, all new physics is encoded in the coefficients (Wilson coefficients) that multiply a fixed set of four-fermion operators: charged-current operators, neutral-current operators with charged leptons, and neutral-current operators with neutrinos. The five scalar leptoquarks—the $SU(2)_L$ singlets $S_1,\\tilde S_1$, the doublets $R_2,\\tilde R_2$, and the triplet $S_3$—each map onto a distinct subset of these coefficients; for instance, $\\tilde S_1$ produces exactly one operator, while $R_2$ and $S_1$ generate both left- and right-handed couplings whose interference is enhanced by the heavy-quark mass in loops. These coefficients, after QCD running of the scalar and tensor currents, feed into the branching-ratio formulas that are compared with data.","core_discovery":"The central claim is that tree-level scalar leptoquark exchange generates a specific pattern of vector, scalar, and tensor four-fermion operators, and that the present upper bounds on rare lepton and kaon decays can be turned into quantitative limits on the corresponding products of Yukawa couplings. The paper works this out for all five scalar leptoquark types, including loop-induced transitions with only leptons or quarks as external states. It finds that $K_L\\to\\mu^+\\mu^-$ and the $K_L\\to\\pi^0\\ell^+\\ell^-$ modes constrain the relevant couplings to order $10^{-5}$ and below (times $(M_\\mathrm{LQ}/\\mathrm{TeV})^2$), while $\\mu\\to e\\gamma$ forces certain left-right top-quark combinations below $10^{-15}$ in the squared coupling, and gold-nucleus $\\mu$-$e$ conversion bounds combinations at the $10^{-11}$ level. The kaon modes are complementary: $K_S$ decays probe the real parts, $K_L$ decays the imaginary parts, and $K^+\\to\\pi^+\\ell^+\\ell^-$ the absolute values of the same coupling combinations.","pith_inferences":["The tabulated bounds are derived one coupling combination at a time, so they should be read as ceilings rather than as a global fit; switching on several Yukawa entries at once could, in principle, produce interference or cancellations that relax the individual limits.","The operator catalog itself is independent of the leptoquark hypothesis: any new physics that generates the same four-fermion Wilson coefficients inherits the same bounds, provided the matching and running are done at the same scales.","The pattern of limits suggests a natural next step: a global fit to all five leptoquark types with simultaneous flavour-structure assumptions, which would convert this catalog into a likelihood and identify which specific Yukawa textures are most constrained."],"forward_implications":["A TeV-mass scalar leptoquark cannot have order-one Yukawa entries in the combinations probed by these modes; the catalog gives the allowed ceiling for each product, so any proposed leptoquark explanation of flavour anomalies must pass these numbers first.","The tightest constraints come from $K_L\\to\\mu^+\\mu^-$, $K_L\\to\\pi^0\\ell^+\\ell^-$, $K_L\\to\\pi^0\\nu\\bar\\nu$, $\\mu\\to e\\gamma$, and $\\mu$-$e$ conversion, so improvements in those channels will most directly shrink the allowed parameter space.","Because the Standard Model and leptoquark amplitudes interfere in $K\\to\\pi\\ell^+\\ell^-$, the differential distributions predicted here can reveal new physics through the $q^2$ shape even before branching-ratio measurements improve.","The limit from $K^+\\to\\pi^+\\nu\\bar\\nu$ forces $K_L\\to\\pi^0\\nu\\bar\\nu$ below about $7.8\\times10^{-10}$ when neutrino flavours are conserved; if the recent candidate events survive, the only scalar-leptoquark way to accommodate them is through neutrino-flavour-violating final states.","The limits scale as $(M_\\mathrm{LQ}/\\mathrm{TeV})^4$ for tree-level processes and $(M_\\mathrm{LQ}/\\mathrm{TeV})^2$ for loop processes, so for $M_\\mathrm{LQ}$ of a few TeV the low-energy bounds weaken but remain competitive with collider searches."],"supporting_citations":[{"why":"supplies the five-leptoquark classification and several earlier decay-mode constraints that this analysis updates","marker":"[7]"},{"why":"gives the interaction Lagrangian conventions and the genuine-leptoquark distinction used in Eq. (1)","marker":"[9]"},{"why":"compiles the experimental branching ratios, decay constants, and hadronic inputs used for every numerical bound","marker":"[61]"},{"why":"sets the $\\mu\\to e\\gamma$ upper bound that dominates the radiative-decay constraints in Table 5","marker":"[82]"},{"why":"provides the gold-nucleus $\\mu$-$e$ conversion upper bound used in Table 7","marker":"[88]"},{"why":"supplies the nuclear conversion-rate formula with overlap integrals used for the $\\mu$-$e$ bounds","marker":"[90]"},{"why":"provides the Standard Model kaon-decay amplitudes and long-distance contributions included in the differential distributions","marker":"[92]"},{"why":"reports the $K^+\\to\\pi^+\\mu^+e^-$ upper bound used in Table 9","marker":"[107]"},{"why":"reports the recent $K^+\\to\\pi^+\\nu\\bar\\nu$ upper bound that drives the constraints in Table 10","marker":"[120]"}],"fun_headline_variants":["Scalar leptoquarks face tight kaon and lepton constraints","Lepton and kaon data squeeze scalar leptoquark couplings","Kaon modes probe scalar leptoquark Yukawa couplings","Rare kaon and lepton decays constrain leptoquark Yukawas","Leptoquark couplings bounded by kaon and muon data"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"Every numerical bound is obtained by letting one leptoquark coupling combination contribute to a given process at a time; if several couplings are present simultaneously, interference or cancellation could weaken or change those limits.","fun_headline_variants_meta":{"raw":{"variants":["Scalar leptoquarks face tight kaon and lepton constraints","Lepton and kaon data squeeze scalar leptoquark couplings","Kaon modes probe scalar leptoquark Yukawa couplings","Rare kaon and lepton decays constrain leptoquark Yukawas","Leptoquark couplings bounded by kaon and muon data"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000775,"raw_usage":{"total_tokens":3423,"prompt_tokens":936,"completion_tokens":2487,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":552,"completion_tokens_details":{"reasoning_tokens":2397}},"tokens_in":552,"tokens_out":2487,"duration_ms":20960,"temperature":1.0,"reasoning_tokens":2397,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T10:22:36.612493+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the $q^2$ distribution of $K^+\\to\\pi^+\\mu^+\\mu^-$ with the current form-factor inputs; the paper predicts a definite interference term between the Standard Model vector amplitude and the leptoquark contribution, so a shape that excludes that interference would rule out the scalar-leptoquark explanation. Alternatively, a confirmed $K_L\\to\\pi^0\\nu\\bar\\nu$ rate above about $7.8\\times10^{-10}$ would falsify the lepton-flavour-conserving scalar-leptoquark picture and force the lepton-flavour-violating neutrino channel.","supporting_citations":[],"review_version":1}