{"id":"0cceb7d4-ab2a-482b-b30d-cfe0b72692c8","arxiv_id":"2603.29125","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"A scalar singlet leptoquark that explains B-meson anomalies produces a ~0.7% decrease in Z→τ+τ−, which future Z-factory measurements could detect, while Z→μ+μ− is essentially unchanged.","lead":"This paper calculates how a hypothetical scalar leptoquark particle would subtly change how often Z bosons decay into tau pairs. The aim is to show that future 'Z factory' colliders could test a popular explanation of B-meson anomalies.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Central −0.7% result rests on unvalidated SloopS implementation; no renormalization scheme or independent cross-check is provided, and the top-quark loop, while consistent with λ1L_bτ, is not documented.","rationale":"The strongest claim is a −0.7% NLO effect in Z→τ+τ− from the scalar singlet leptoquark. I agree with the reader that the main vulnerability is the black-box SloopS implementation: no renormalization scheme is specified, no UV-finiteness check is shown, and no cross-check against the SMEFT matching (Ref. [26]) is given. However, the reader's specific charge of inconsistency about the 'internal top quark' is not correct in my reading: the Lagrangian in Eq. (2.1) has λ1L_{iα} \\bar{q}_i^c ε l_α, and for i=3, α=3 this couples both b_L to ν_τ and t_L to τ_L. Thus a top-quark triangle loop is expected and consistent with the stated couplings. The notation λ1L_bτ is just the (3,3) matrix element. That said, the paper does not explicitly state this, and the ambiguity itself is a minor clarity issue. The deeper concern is quantitative: for λ1L=1.4, λ²/(16π²)≈1.2%, so the claimed 0.7% is not a tiny correction; an unsuppressed counterterm or a γ5-scheme artifact in the automated code could produce a shift of the same order. The paper's decision to bypass the SMEFT route and use a private automated code—without input files or an independent test—makes the result unverifiable from the published text. Keeping the verdict CONDITIONAL (unchanged) is the right call, with the concrete test of an independent SMEFT-based recomputation.","tokens_in":13492,"tokens_out":21490,"duration_ms":224750,"concrete_test":"Recompute δ for BP0 (MS1=1 TeV, λ1L_bτ=1.4, λ1R_cτ=−0.1) using the one-loop SMEFT matching of Ref. [26] to extract the effective Zττ coupling modification and translate it into the Z→τ+τ− width shift, or equivalently generate the S1 model with the public FeynRules implementation of Ref. [37] and run it through an independent NLO framework (e.g., MadGraph5_aMC@NLO). If the resulting δ differs from the paper's −0.65% by more than 0.1 percentage points, the SloopS calculation is not trustworthy; agreement within integration errors would resolve the concern.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The single load-bearing assumption is that the SloopS automated one-loop calculation correctly implements the scalar singlet S1 model and its renormalization for Z→τ+τ−. The paper gives no renormalization conditions, no demonstration of UV-finite results after counterterms, and no comparison with the SMEFT one-loop matching of Ref. [26], which should reproduce the same Zττ vertex shift. The reader's specific worry about an 'internal top quark' is likely a misreading: because λ1L_bτ is the (3,3) element of the left-handed coupling matrix, it gives both b−ν_τ and t−τ vertices, so a top-quark loop is expected. The genuine issue is that the −0.7% central value is of the same order as the loop factor λ²/(16π²) for λ≈1.4, so any mishandled counterterm or γ5-scheme artifact could shift δ by the full size of the signal. Without an independent numerical cross-check or at least the explicit SloopS input and renormalization setup, the central claim depends on a black-box code.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript computes the one-loop electroweak corrections to Z->tau+tau-, Z->mu+mu-, and e+e- -> tau+tau-/mu+mu- induced by the scalar singlet leptoquark S1 in the minimal scenario with only lambda1L_btautau and lambda1R_ctautau non-zero. Using the SloopS automated framework, it finds that the mu-pair channel is negligible (O(10^-6)% ), while the tau-pair channel receives a correction of up to about -0.7%, governed mainly by lambda1L_btautau, for both 1 TeV and 2 TeV leptoquark masses. The paper also provides a fitted analytic function (Eq. 3.1), maps current and expected Z-factory constraints into the coupling plane, and argues that the e+e- result at the Z pole is essentially identical to the Z-decay result. The central phenomenological message is that a leptoquark addressing the charged-current B anomalies may be visible in Z->tau+tau- at a future Z factory despite the loop suppression.","tokens_in":13785,"tokens_out":7378,"duration_ms":80367,"significance":"If the numerical result is correct, the paper presents a useful and non-obvious target: a leptoquark motivated by B-physics anomalies can produce an observable ~0.1%-0.7% shift in Z->tau+tau- at FCC-ee/CEPC, with the large allowed left-handed coupling compensating the TeV-scale mass suppression. The paper's strengths include a genuine one-loop calculation rather than a leading-order effective-operator estimate, a clear parameter scan with current constraints, a convenient fitting formula for the coupling-mass dependence, and explicit numerical tables for the Z-pole collider. The main weakness is that the central result rests entirely on the SloopS automated implementation, with no renormalization conditions, no explicit UV-finiteness check, and no independent numerical cross-check; this makes the -0.7% claim currently impossible to verify from the manuscript alone.","major_comments":[{"comment":"The central numerical result, delta ~ -0.7%, is obtained solely from the SloopS automated one-loop implementation. The manuscript does not state the renormalization scheme or counterterm structure for the S1 model, does not demonstrate that the one-loop result is UV-finite after renormalization, and provides no independent cross-check. Since delta is itself a one-loop quantity, an error in the counterterms, the gamma5 scheme, or the treatment of the colored scalar could shift the result by its full size. Please provide the renormalization conditions, show the cancellation of UV poles or the residual scale dependence, and, if possible, compare at least one benchmark against an independent calculation or against the SMEFT one-loop matching of Ref. [26].","section":"Sections 2 and 3, Eq. (2.1)-(2.3)"},{"comment":"Eq. (3.1) is a fit to the same numerical points that are used to obtain the delta values; using this fitted function to interpret the future Z-factory precision in Fig. 3 is an interpolation-based inversion, not an independent analytic prediction. The fit quality and the validity range in (lambda1L_btautau, MS1) should be reported, and the function should not be used for extrapolation beyond the scanned region. The current presentation of Eq. (3.1) as an 'analytic function to quantify the LQ effects' is over-stated.","section":"Eq. (3.1), Fig. 5, Fig. 3"},{"comment":"There is a sign inconsistency in the central numerical output. Eq. (2.3) defines delta = (sigma_NLO_S1 - sigma_NLO_SM)/sigma_LO, so for negative Delta_sigma_NLO_S1 the value of delta should be negative. In Table 3, however, Delta_sigma_NLO_S1 is negative while delta_BP0 and delta_BP1 are quoted as positive percentages. Figures 3 and 7 also use positive-looking color/axis scales while the text and abstract state a maximum deviation of about -0.7%. If the plotted/tabulated quantity is |delta|, this must be stated explicitly; otherwise the signs should be corrected. This is not merely cosmetic, because the sign of the shift is physically relevant for asymmetry observables.","section":"Table 3 and Eq. (2.3)"},{"comment":"Eq. (4.1) states that the Z-pole collider correction equals the Z-decay correction. This ignores the fact that a shift in Gamma_tau_tau also shifts the total width Gamma_Z entering sigma_peak = 12 pi Gamma_ee Gamma_ff/(M_Z^2 Gamma_Z^2). The collider delta differs from the decay delta by an approximate factor (1 - 2 Gamma_ff/Gamma_Z) ~ 0.93 for f=tau. The manuscript should either present the exact expression or state this approximation and its numerical impact. This may partly explain why the quoted peak values (0.60-0.63% in Table 3) are slightly below the maximum -0.7% quoted for the decay.","section":"Section 4, Eq. (4.1)"}],"minor_comments":[{"comment":"The sentence 'The latter effects would be sensitive to the top quark mass by introducing the internal top quark in triangle loop' is unexplained. A top-quark loop is indeed expected because lambda1L_btautau is the (3,3) element of the left-handed coupling matrix and gives a t-tau vertex, but this should be stated explicitly so the loop-fermion content is clear.","section":"Section 3, Fig. 2"},{"comment":"The fit parameters K2, K1, Kd are given without units or statistical errors. State the units (e.g., TeV for masses) and the fit range/residuals so the formula can be used reliably.","section":"Eq. (3.1)"},{"comment":"The manuscript mixes positive and negative signs for delta between the text, figures, and tables. If the intended quantity is the absolute value (as the color bars in Figs. 3 and 7 suggest), label it |delta| and make the convention uniform.","section":"Throughout"},{"comment":"There are several typos: 'Bellec' and 'LHCba' in Table 1, 'the 2th-generation leptons' in Section 3, and 'Combin' in Ref. [14]. Also, the reference for the SloopS-based Higgs-strahlung paper [65] should include the final publication details if available.","section":"Minor text"},{"comment":"The figure caption says 'angular of final tau+' but the text says 'final state tau-'. Please make this consistent.","section":"Section 4, Fig. 8"}],"recommendation":"major_revision","confidential_remarks":"The paper is within the journal's scope and addresses an interesting phenomenological target. My main concern is that the -0.7% result is produced entirely by an automated code with no renormalization details and no independent cross-check; this should be addressed before publication. The sign inconsistency between Eq. (2.3) and Table 3 is easy to fix but suggests the numerical presentation needs careful proofreading. I am not treating the top-quark loop as an error, since lambda1L_btautau gives a t-tau vertex, but the manuscript should clarify this. Overall this is a defensible major-revision case: the central claim is plausible but currently rests on an unvalidated black-box calculation."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper computes NLO electroweak corrections from the scalar singlet leptoquark S1 in Z→τ+τ− and e+e−→τ+τ− at the Z pole, in the minimal charged-current anomaly scenario. The main result is that the τ-pair channel can receive a correction up to about −0.7%, driven by λ1L_bτ, while the μ-pair channel is negligible. That is a new numerical result, not present in the cited literature, and it is a plausible and useful extension of leptoquark phenomenology to future Z-factory precision measurements.\n\nWhat the paper does well: it follows an established automated pipeline (SloopS), correctly exploits the Z-pole relation between decay width and cross section, and includes the relevant constraints from R(D(*)), Bc, |gτ/gμ|, and LHC searches. The parameter scan is transparent, and the analytic fit in Eq. 3.1 is explicitly labeled as a fit, with curves checked against the numerical points for several masses. The bottom line—that the enlarged allowed parameter space at higher LQ mass compensates the mass suppression—is clearly demonstrated.\n\nSoft spots, in proportion: the biggest one is that the renormalization scheme is not described. The paper does not give the counterterm conditions, show UV finiteness, or provide any independent cross-check. Because the signal is of the same order as the loop factor for the benchmark couplings, a mishandled counterterm or γ5 scheme artifact could shift the central value by the full size of the effect. The authors should disclose the renormalization setup and ideally compare with the one-loop SMEFT matching of Ref. [26]. Second, Eq. 3.1 is a fit to the same numerical results that are then used to project sensitivities; that is not a fatal circularity, but it means the analytic function is a parametrization, not a derivation. Third, there are no systematic uncertainty estimates beyond the Monte Carlo integration errors shown in Table 3; for a 0.1%-level prediction, that would be worth addressing. The concern raised in the reader's report about an 'internal top quark' in the triangle loop is actually resolved when one remembers that λ1L_bτ is the (3,3) element of the left-handed coupling matrix: because the left-handed quark doublet contains both b and t, the top quark does couple through the same λ1L_bτ. The paper would be clearer if it stated this explicitly, but it is not an inconsistency.\n\nWho is this for? People working on leptoquark phenomenology and on Z-factory physics projections. It is not a breakthrough, but it is a clean, useful calculation that deserves a serious referee. My recommendation: send it to peer review, and ask the authors to add renormalization details and, if possible, an independent numerical check or input files.","headline":"Solid, incremental NLO EW calculation of scalar singlet leptoquark effects in Z→τ+τ−; worth refereeing if the authors disclose the renormalization setup.","tokens_in":14295,"tokens_out":2212,"would_cite":false,"duration_ms":24463,"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":"A scalar singlet leptoquark that explains B-meson anomalies leaves a −0.7% imprint in Z→τ+τ−, within reach of future Z factories.","keywords":["leptoquark","scalar singlet leptoquark","Z factory","Z→τ+τ−","NLO electroweak corrections","charged-current anomalies","lepton flavor universality","tau-pair production"],"falsifier":"An independent one-loop renormalization of Z→τ+τ− in the same two-coupling S1 model would settle the numerics: if the correction is not quadratic in λ1L_bτ, or not capped near −0.7% at the allowed couplings, the central claim fails. On the experimental side, measuring Γ(Z→τ+τ−)/Γ(Z→µ+µ−) at 0.1% precision at a future Z factory would expose a −0.7% tau-width shift if the scenario is right; a null result would rule it out.","tokens_in":13361,"feed_emoji":"⚛️","tokens_out":7883,"duration_ms":76789,"temperature":0.7,"pith_summary":"The paper asks whether a scalar singlet leptoquark built to explain the B-meson charged-current anomalies leaves a measurable trace in Z-boson decays and e+e− collisions at the Z pole. It finds that the tau-pair channel receives a next-to-leading-order correction of up to about −0.7% for leptoquark masses of 1 and 2 TeV, while the muon-pair channel is essentially unchanged. The effect is driven almost entirely by the left-handed coupling of the leptoquark to the bottom quark and tau lepton, so a precision Z factory could directly constrain that coupling. The paper also provides a simple analytic fit to the correction as a function of mass and coupling, making the prediction easy to test.","feed_headline":"Leptoquark shifts Z decays to tau pairs by 0.7%","feed_subtitle":"A future precision Z factory could turn the 0.7% shift into a direct constraint on the leptoquark's left-handed coupling.","key_machinery":"The scalar singlet leptoquark S1 — a color-triplet scalar coupling a quark to a lepton — enters the Z→τ+τ− amplitude at one loop through vector-boson self-energies, the tau self-energy, and the Zττ vertex. The paper's quantitative handle is the observable δ: the S1-induced NLO shift divided by the leading-order Standard Model rate or width. The effect is controlled by λ1L_bτ, with the fit function δ_fitted(λ1L_bτ, MS1) = (λ1L_bτ)^2 [K2/(MS1+Kd)^2 + K1/(MS1+Kd)] and fitted constants K2 = −0.5919, K1 = −0.03947, Kd = 0.4188.","core_discovery":"In the minimal scalar singlet leptoquark model built to explain the b→cτν anomalies, where only λ1L_bτ and λ1R_cτ are non-zero, the paper finds that one-loop electroweak corrections shift Z→τ+τ− by up to about −0.7% for both 1 TeV and 2 TeV leptoquark masses. The shift grows quadratically with λ1L_bτ, is almost independent of λ1R_cτ, and is identical for e+e−→τ+τ− at the Z pole. Muon-pair final states receive only negligible corrections. The paper provides an analytic fit, δ_fitted(λ1L_bτ, MS1), and uses it to translate projected 0.1–0.3% Z-factory precisions into limits on λ1L_bτ.","pith_inferences":["The same loops that shift Z→τ+τ− should also generate S1 corrections to Z→bb and Z→cc (or b/c-pair production at the Z pole), giving an independent experimental cross-check the paper does not compute.","The paper's remark that the vertex is 'sensitive to the top quark mass' does not match its stated λ1R_cτ and λ1L_bτ couplings to charm and bottom; if the internal quark is charm or bottom rather than top, the mass-sensitivity and the fitted coefficients could shift.","Because the correction is essentially independent of λ1R_cτ, a tau-pair precision measurement would not bound the right-handed coupling; combining with B-decay measurements would be needed to separate the two couplings.","The fitted analytic function turns any future measurement of the tau-to-muon width ratio into a direct bound on λ1L_bτ across S1 masses, effectively making the Z factory a parameter-light tester of the B-anomaly explanation."],"forward_implications":["If Z factories reach 0.1–0.2% precision on tau-pair observables, they will probe the currently allowed λ1L_bτ region for S1 masses up to at least 2 TeV.","The equality of the Z-decay and e+e− collision effects at the Z pole means the same calculation can be compared with both lineshape and forward–backward asymmetry measurements.","Muon-pair channels will look exactly like the Standard Model, so they will not constrain the scalar leptoquark in this scenario.","The effect's stability across pT, rapidity, and cosθ means integrated rate measurements, not just differential shapes, can capture the full signal.","The dominance of λ1L_bτ means the tau-pair channel is a direct handle on the left-handed coupling responsible for the B anomalies."],"fun_headline_variants":["Z factory's tau pairs reveal leptoquark's 0.7% signal","Leptoquark effect in tau pairs: a 0.7% window for new physics","Tau pair shifts at Z pole expose leptoquark coupling limits","Future Z factory can pin down leptoquark via tau pairs","0.7% tau shift at Z pole could constrain leptoquark model"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The load-bearing premise is that the automated one-loop electroweak calculation is correct — the renormalization scheme, the Zττ vertex, and the quark masses in the loops — and the paper's remark that the vertex is 'sensitive to the top quark mass' conflicts with its own non-zero couplings to charm and bottom, leaving the loop content not fully pinned down.","fun_headline_variants_meta":{"raw":{"variants":["Z factory's tau pairs reveal leptoquark's 0.7% signal","Leptoquark effect in tau pairs: a 0.7% window for new physics","Tau pair shifts at Z pole expose leptoquark coupling limits","Future Z factory can pin down leptoquark via tau pairs","0.7% tau shift at Z pole could constrain leptoquark model"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000721,"raw_usage":{"total_tokens":3105,"prompt_tokens":811,"completion_tokens":2294,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":555,"completion_tokens_details":{"reasoning_tokens":2193}},"tokens_in":555,"tokens_out":2294,"duration_ms":13863,"temperature":1.0,"reasoning_tokens":2193,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-02T17:04:04.526343+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"An independent one-loop renormalization of Z→τ+τ− in the same two-coupling S1 model would settle the numerics: if the correction is not quadratic in λ1L_bτ, or not capped near −0.7% at the allowed couplings, the central claim fails. On the experimental side, measuring Γ(Z→τ+τ−)/Γ(Z→µ+µ−) at 0.1% precision at a future Z factory would expose a −0.7% tau-width shift if the scenario is right; a null result would rule it out.","supporting_citations":[],"review_version":1}