{"id":"b72f81ab-9a17-4fb8-a075-6c77232ad81a","arxiv_id":"1909.00403","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 leptoquark S1 near 1 TeV can raise B to mu nu to the Belle central value, but required Yukawa couplings are large and several LHC and EDM constraints already cut into that parameter space.","lead":"Researchers show that a hypothetical 1 TeV particle called a scalar leptoquark could boost the rate of the rare B meson decay to a muon and a neutrino, bringing it in line with a slightly high Belle measurement. The paper maps how large the leptoquark's couplings must be and which LHC and low-energy experiments currently leave those couplings unconstrained.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The claimed surviving parameter space for S1 in B→μν depends on an approximate ATLAS recast whose correlated branching fractions and lack of nuisance parameters make the exclusion boundary insecure.","rationale":"The reader's CONDITIONAL verdict is well matched to the paper. The matching calculation is internally consistent, the two-phase fits reproduce the quoted Wilson-coefficient values, and the paper honestly labels its simplifying assumptions. The reader's weakest_assumption concerns the truncation of S1 Yukawa couplings to a small set; that is a scope limitation rather than a threat to the central claim, since the paper explicitly explores only a slice of parameter space and does not claim full coverage. By contrast, the approximate treatment of the LHC constraints directly supports the assertion that part of the parameter space remains viable, which is part of the strongest claim as summarized. The recast limitation is acknowledged in the text (Sec. IV F, footnote 3), and the correlated branching fractions make the naive comparison to the ATLAS limit potentially unreliable. Still, the qualitative conclusion that an S1 of mass around 1 TeV can significantly modify B→μν is robust: the required couplings are small, and some low-coupling region will likely survive even after a more careful recast. Therefore the reader's CONDITIONAL verdict should stand, and no further verdict change is needed. The concrete test proposed would sharpen the quantitative exclusion boundary and resolve whether the surviving-region claim is correct.","tokens_in":15718,"tokens_out":16718,"duration_ms":157385,"concrete_test":"Recompute the purple exclusion region of Fig. 3 (left) with a full simulated recast of S1 pair production at sqrt(s)=13 TeV for m_S1=1.2 TeV, including all decay modes (S1→uμ, tμ, bν_μ) with branching fractions determined by the same y^R_12 and y^L_32 couplings, applying the ATLAS [16] event selection and its likelihood with nuisance parameters; then overlay the Belle 1σ and central-value contours for B(B→μν) to see whether any of the central-value contour survives. A minimal version of this check is to compute the single-muon pair-production yield 2 B(μ)[1−B(μ)] σ_pair ε for each point and compare it directly with the observed 95% CL upper limit from [16], rather than comparing B(S1→uμ) to the quoted branching-fraction limit.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The load-bearing step is not the matching calculation but the statement that after LHC constraints a part of the S1 Yukawa parameter space remains for B→μν. This rests on two recasts in Sec. IV F: the ATLAS pair-production search [16] used for the purple regions in Fig. 3, and the ditau search [48] used for Fig. 4. Both are explicitly approximate: the test statistic in Eq. (25) uses ΔN_i = sqrt(N_obs) with no nuisance parameters (footnote 3), and the pair-production constraint is applied by comparing B(S1→uμ) directly to the 95% CL upper limit on B(LQ→qμ) from [16]. But for the two-coupling scenario (y^R_12, y^L_32), S1 has several decay modes—uμ, tμ, bν_μ—whose branching fractions are correlated with the same couplings that set the B→μν rate. These modes have different jet multiplicities, b-tag content, and muon/MET acceptances, so the signal rate is not simply proportional to B(S1→uμ). The Belle central-value contour for φ=π corresponds to a fairly sizable B(S1→uμ) (order 0.4), so a factor-of-two shift in the effective limit could remove or significantly shrink the surviving region. Thus the quantitative claim 'a part of the parameter space remains unexcluded' is not yet secure, even though the core observation that S1 can modify B→μν is sound.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper studies the impact of the scalar leptoquark S1 on the purely leptonic decays B→μν̄ and B→τν̄. Starting from the S1 Yukawa Lagrangian, the authors derive the effective Hamiltonian and Wilson coefficients at the S1 mass scale, including RGE running to the B-meson scale. Using the Belle measurement of B(B→μν) and the world average for B(B→τν), they extract the required magnitudes and phases of scalar Wilson coefficients, highlighting the m_B^2/(m_l m_b) enhancement of scalar over vector operators. They then translate these requirements into the S1 Yukawa couplings under the assumption that all couplings except selected products vanish or are CKM-suppressed, and overlay constraints from b→cℓν observables, (g−2)_μ, τ→μγ, neutron EDM, and LHC direct searches. The main quantitative claim is that for m_S1 ≈ 1.2 TeV, a part of the S1 Yukawa parameter space that can explain the Belle central value remains unexcluded after LHC and EDM constraints.","tokens_in":16055,"tokens_out":6260,"duration_ms":54772,"significance":"The analytic matching in Eqs. (5)–(8) is standard and internally consistent, and the paper correctly identifies the scalar enhancement factor m_B^2/(m_l m_b) and its implications for B→μν. The compilation of flavor and collider constraints, including an explicit MC recast for the ditau search, is useful. The paper's central claim that sizable S1 couplings remain viable, however, rests on approximate recasts: the ATLAS pair-production constraint is applied as a simple scaling by B(S1→uμ), and the ditau recast uses a simplified test statistic without nuisance parameters. These limitations are partially self-flagged, but they directly affect the surviving-region conclusion. The results are nonetheless a reasonable first map of the S1 parameter space for B→μν, provided the approximations are treated as indicative rather than definitive.","major_comments":[{"comment":"The ATLAS pair-production constraint is implemented by comparing B(S1→uμ) directly with the 95% CL upper limit on B(LQ→qμ) from Ref. [16]. In the two-coupling scenario (y^R_12, y^L_32), S1 also decays to tμ and bν_μ, with branching fractions correlated to the same couplings; these modes have different jet multiplicities, b-tag content, and lepton/MET acceptances, so the signal rate is not simply proportional to B(S1→uμ). For φ=π, the Belle central-value contour corresponds to B(S1→uμ) ~ 0.4, and a factor-of-two change in the effective limit could remove or substantially shrink the surviving region. The exclusion boundaries in Fig. 3 should therefore be treated as approximate; a full recast including all decay modes is needed to secure the claim that part of the parameter space remains unexcluded.","section":"Sec. IV.F, Fig. 3"},{"comment":"The ditau recast uses the test statistic with ΔN_i = sqrt(N_obs) and no nuisance parameters, as acknowledged in footnote 3. This ignores correlated systematic uncertainties and can bias the exclusion boundary; the light-gray excluded region in Fig. 4 is correspondingly not robust at the quoted 2σ level. Since this constraint is used to limit the τ-related parameter space, the quantitative boundary should be revisited with a more complete likelihood, or the paper should state a conservative lower bound rather than a precise exclusion contour.","section":"Sec. IV.F, footnote 3, Eq. (25)"},{"comment":"The analysis assumes all S1 Yukawa couplings other than the plotted products are zero (y^L(R)_i1 = 0) or do not overpower CKM suppression. Under this truncation, the Wilson coefficients in Eq. (5) reduce to y^R_1l y^L_3l' terms, and the LHC decay modes are restricted to uℓ, tℓ, bνℓ. If, for instance, y^L_2i or y^L_3i couplings are non-negligible, both the Wilson coefficients entering Eq. (7) and the S1 branching fractions used in Sec. IV.F change, so the quoted coupling ranges in Figs. 3 and 4 are not a fully general characterization of the S1 parameter space. This limitation is stated, but its impact on the surviving-region claim should be quantified or made more prominent.","section":"After Eq. (5) and Sec. IV.A"}],"minor_comments":[{"comment":"The word 'in prinicple' should be corrected to 'in principle'.","section":"Abstract"},{"comment":"The heading 'YUKA W A COUPLINGS' appears to be a typo; it should read 'YUKAWA COUPLINGS'.","section":"Section IV heading"},{"comment":"The reference to 'Sec. 4.6' should be updated to 'Sec. IV.F' for consistency with the section numbering.","section":"Sec. IV.F"},{"comment":"The sentence 'the black colors does not correspond to any phase labeling' contains an agreement error; it should be 'the black color does not correspond to any phase labeling.'","section":"Figure 3 caption"},{"comment":"The extraction of the required Wilson coefficients is a fit to the Belle central value rather than a prediction; the wording 'can modify the B→μν rate significantly' in the abstract should be qualified as 'can accommodate the measured central value' to avoid overstatement.","section":"Eq. (11) and Sec. III.A"}],"recommendation":"major_revision","confidential_remarks":"The paper is within the journal's scope and the central mechanism is sound. The main weakness is the reliance on approximate LHC recasts for the claim of surviving parameter space; this should be addressed by either a more complete recast or by clearly downgrading the claim. I do not see grounds for rejection, but the quantitative conclusions need revision."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The useful core here is the application of the S1 leptoquark to the muon channel specifically, including complex-phase dependence. The matching calculation is standard and internally consistent, and the paper correctly emphasizes the m_B^2/(m_l m_b) enhancement of scalar over vector operators. It also honestly flags its own simplifications: the y_i1 couplings set to zero, the approximate ditau recast, and the fact that the central coupling values are fit to Belle's central branching ratio rather than predicted. That candor is worth something.\n\nThe main soft spot is the one the stress test identifies: the claim that part of the Yukawa parameter space survives LHC constraints rests on two recasts that are not robust at the factor-of-two level. For the two-coupling scenario, S1 decays to uμ, tμ, and bν, and the signal rate is not simply proportional to B(S1→uμ). The ATLAS pair-production limit is applied by direct comparison of that single branching fraction, and the ditau recast uses ΔN = sqrt(N_obs) with no nuisance parameters. For the φ=π contour, where B(S1→uμ) is sizable, a factor-of-two shift in the effective limit could remove the surviving region. So the qualitative statement that S1 can modify B→μν is sound, but the quantitative \"parameter space remains\" claim is not yet secure.\n\nMinor points: the truncation to only two nonzero Yukawas at a time is stated explicitly, and the bounds from b→c transitions and tau→mu gamma are checked and found weak in the plotted ranges—that is done honestly. The paper also correctly notes that prior Ref. [11] already covered B→tau nu, so the novelty is the muon application and the LHC/EDM interplay, not the framework itself.\n\nWho is this for? People working on leptoquark phenomenology who want a quick map of what S1 can do for B→mu nu before Belle II. It is a useful target-setting paper, but the exclusion boundaries should be treated as indicative, not precise. A serious referee would want the recasts redone with proper correlated branching fractions and, ideally, validated against the ATLAS public limits.\n\nMy recommendation: send it to peer review. The central idea is legitimate, the math is clean, and the flaws are explicitly disclosed approximate recasts rather than hidden errors. With modest revision—tone down the exclusion claims, or replace the recast with a more careful one—it would be a solid contribution.","headline":"A competent, honest application of S1 leptoquark to B→μν with a standard matching calculation, but the LHC recast that secures the surviving parameter space is too crude to trust for quantitative exclusion regions.","tokens_in":16579,"tokens_out":634,"would_cite":false,"duration_ms":8071,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["13.20.He","12.60.-i","14.80.Sv"],"model":"deepseek-v4-flash","headline":"A TeV-scale scalar leptoquark can modify the B→μν decay rate enough to match the measured central value.","keywords":["scalar leptoquark","S1 leptoquark","B to mu nu decay","purely leptonic B decays","Wilson coefficients","helicity suppression","neutron electric dipole moment","future collider searches"],"falsifier":"A future measurement at a B factory that finds $\\mathcal{B}(B \\to \\mu \\bar{\\nu})$ within $1\\sigma$ of the Standard Model value $3.92 \\times 10^{-7}$ with an uncertainty below roughly $0.5 \\times 10^{-7}$ would remove the empirical motivation for the S1 explanation; a future neutron-EDM bound $|d_n/e| < 10^{-28}$ cm would specifically kill the CP-violating $\\phi_{\\tau\\tau} = \\pi/2$ branch of the $B \\to \\tau \\bar{\\nu}$ solution.","tokens_in":15479,"feed_emoji":"⚛️","tokens_out":9551,"duration_ms":82131,"temperature":0.7,"pith_summary":"Purely leptonic B decays are a clean window onto new physics, and this paper argues that the scalar leptoquark S1 at roughly 1 TeV can open that window wide in the muon channel. Although the Standard Model rate for B→μν̄ is helicity-suppressed, S1 generates a scalar Wilson coefficient that interferes with the Standard Model amplitude, and the scalar operator is enhanced by $m_B^2/(m_\\mu m_b) \\sim 60$ relative to the vector operator. Matching the measured branching ratio of $5.3 \\times 10^{-7}$ requires $|C_{\\mathrm{SL}}^{\\mu\\mu}| \\simeq 0.0344$ at relative phase $0$, or $\\simeq 0.0026$ at phase $\\pi$; a tau antineutrino in the final state needs $|C_{\\mathrm{SL}}^{\\mu\\tau}| = 0.009$. If the paper is right, a TeV-scale S1 with Yukawa couplings that are still allowed by direct searches and electric dipole moments is a concrete explanation of the current $2.8\\sigma$ excess, and the next round of B-factory data should be able to distinguish it from the Standard Model.","feed_headline":"A scalar leptoquark can explain the high B→μν rate","feed_subtitle":"A 1-TeV S1 shifts the decay through a helicity-enhanced scalar operator; EDM and collider limits still leave room.","key_machinery":"The load-bearing object is the scalar Wilson coefficient $C_{\\mathrm{SL}}^{\\ell\\ell'}$ generated by the TeV-scale leptoquark in the effective Hamiltonian $\\mathcal{H}_{\\mathrm{eff}} = \\frac{4G_F}{\\sqrt{2}} V_{ub}(C_{VL}O_{VL} + C_{SL}O_{SL})$. The coefficient is the product $y^{R*}_{1\\ell} y^L_{3\\ell'}/m_{S1}^2$, run down to the B-meson scale, and it works because the scalar operator $O_{\\mathrm{SL}}$ is helicity-enhanced in $B \\to \\ell \\bar{\\nu}$ by $m_B^2/(m_\\ell m_b)$; the phase between this coefficient and $V_{ub}$ controls whether the S1 contribution interferes constructively or destructively with the Standard Model amplitude. This one coefficient, rather than the full leptoquark Lagrangian, is what lets the paper translate the measured rate into required Yukawa coupling ranges.","core_discovery":"The paper's central claim is that the scalar leptoquark S1, with mass around 1.2 TeV, can reproduce the observed B→μν̄ rate through the scalar operator $O_{\\mathrm{SL}}^{\\ell\\ell'} = (\\bar{u}_L b)(\\bar{\\ell}_L \\nu_{\\ell'})$, not through the vector operator. The S1 contribution to the Wilson coefficient is $C_{\\mathrm{SL}}^{\\ell\\ell'} = -\\frac{\\sqrt{2}}{8G_F V_{ub}}\\frac{y^{R*}_{1\\ell} y^L_{3\\ell'}}{m_{S1}^2}$, and it enters the branching ratio with a relative phase and with the helicity enhancement $m_B^2/(m_\\ell m_b)$. For the muon channel this enhancement is about 60, so the required coefficient is an order of magnitude smaller than the vector-coefficient alternative. The paper shows quantitatively that $|C_{\\mathrm{SL}}^{\\mu\\mu}| \\lesssim 0.0344$ at phase $0$ or $\\lesssim 0.0026$ at phase $\\pi$ matches the measured central value, and that $|C_{\\mathrm{SL}}^{\\mu\\tau}| = 0.009$ does the same by adding in quadrature when the unseen antineutrino is a $\\tau$. After imposing S1 pair-production searches and neutron-EDM constraints, part of the corresponding Yukawa parameter space remains.","pith_inferences":["Editorial extension: the ratio $R^{\\mu/\\tau}_B = \\mathcal{B}(B \\to \\mu \\bar{\\nu})/\\mathcal{B}(B \\to \\tau \\bar{\\nu})$, which the paper notes is 0.0045 in the SM, would be a nearly parameter-free cross-check: both decays are set by the same $y^R_{1\\ell} y^L_{3\\ell'}$ products, and their ratio is almost insensitive to $f_B$ and $|V_{ub}|$.","Editorial extension: the paper truncates to $y^R_{1\\ell}$ and $y^L_{3\\ell'}$; allowing non-negligible $y^L_{2i}$ or $y^L_{3i}$ would change the Wilson coefficients, the inferred coupling ranges, and the LHC branching-fraction interpretation, so a global fit with those couplings switched on is a natural next step.","Editorial extension: if the R(D^{(*)}) anomalies persist, the same $y^L_{33}$ that drives $b \\to c\\tau\\nu$ would also, when paired with $y^R_{12}$, feed $B \\to \\mu \\bar{\\nu}$ through $C_{\\mathrm{SL}}^{\\mu\\tau}$; a combined flavor and dijet analysis could connect the two anomalies."],"forward_implications":["If the measured B→μν̄ rate remains above the Standard Model, an S1 leptoquark near 1.2 TeV with the plotted Yukawa products is a viable explanation; Belle II should confirm the rate with early data.","A future neutron-EDM bound near $|d_n/e| < 10^{-28}$ cm would exclude the CP-violating $\\phi_{\\tau\\tau} = \\pi/2$ region that supports the higher $B \\to \\tau \\bar{\\nu}$ rate, while leaving the $\\phi = 0, \\pi$ solutions.","HL-LHC searches for $S_1 \\to t\\mu$ and heavy resonances in $\\tau\\tau$, $\\tau\\nu$, and $\\mu\\nu$ final states can cover much of the still-open Yukawa parameter space, especially if $|y^L_{33}|$ is as large as the R(D^{(*)}) tension favors.","The paper's comparison shows that the $B \\to \\mu \\bar{\\nu}$ channel is the sharper probe among purely leptonic B decays: the same S1 couplings need much larger values to move $B \\to \\tau \\bar{\\nu}$, because the scalar enhancement is only about 4 there."],"supporting_citations":[{"why":"Supplies the exclusive |Vub| and particle data inputs that set the SM branching ratio.","marker":"[2]"},{"why":"Provides the SM expectation 3.92e-7 and the Belle-style Wilson-coefficient presentation this paper extends to the S1 leptoquark.","marker":"[5]"},{"why":"Reports the B→μν̄ measurement that motivates studying the decay.","marker":"[6]"},{"why":"Gives the updated central value 5.3e-7 and the 1σ/2σ ranges that the paper reproduces with S1.","marker":"[7]"},{"why":"Supplies the neutron-EDM effective Hamiltonian and S1 loop contributions used to constrain complex phases and motivate setting first-generation couplings to zero.","marker":"[11]"},{"why":"Provides the S1 pair-production exclusion that cuts the allowed y^R_{12}-y^L_{32} plane.","marker":"[16]"},{"why":"Electron EDM limit that motivates setting the first-generation couplings to zero.","marker":"[19]"},{"why":"FLAG value f_B = 190 MeV fixing the SM B→μν̄ normalization.","marker":"[20]"},{"why":"Ditau resonance search that provides the 2σ exclusion on |y^R_{13}| shown in the τν̄ plots.","marker":"[48]"}],"fun_headline_variants":["Scalar leptoquark tweaks B→μν rate at 1 TeV","Helicity boost lets S1 match B→μν data","S1 leptoquark passes EDM and LHC limits","B→μν anomaly explained by scalar operator","1-TeV leptoquark fits B→μν without vector term"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The analysis assumes that all S1 Yukawa couplings except the plotted $y^R_{1\\ell}$ and $y^L_{3\\ell'}$ products are zero (to satisfy electron EDM) or CKM-suppressed, so the Wilson coefficients simplify to those products; if those other couplings are not small, the quoted ranges and exclusion plots no longer describe the full S1 parameter space.","fun_headline_variants_meta":{"raw":{"variants":["Scalar leptoquark tweaks B→μν rate at 1 TeV","Helicity boost lets S1 match B→μν data","S1 leptoquark passes EDM and LHC limits","B→μν anomaly explained by scalar operator","1-TeV leptoquark fits B→μν without vector term"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000394,"raw_usage":{"total_tokens":2082,"prompt_tokens":971,"completion_tokens":1111,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":587,"completion_tokens_details":{"reasoning_tokens":1021}},"tokens_in":587,"tokens_out":1111,"duration_ms":8446,"temperature":1.0,"reasoning_tokens":1021,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T05:54:02.463657+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A future measurement at a B factory that finds $\\mathcal{B}(B \\to \\mu \\bar{\\nu})$ within $1\\sigma$ of the Standard Model value $3.92 \\times 10^{-7}$ with an uncertainty below roughly $0.5 \\times 10^{-7}$ would remove the empirical motivation for the S1 explanation; a future neutron-EDM bound $|d_n/e| < 10^{-28}$ cm would specifically kill the CP-violating $\\phi_{\\tau\\tau} = \\pi/2$ branch of the $B \\to \\tau \\bar{\\nu}$ solution.","supporting_citations":[{"cited_title":"Enhanced $B \\to \\mu \\bar\\nu$ Decay at Tree Level as Probe of Extra Yukawa Couplings","cited_arxiv_id":"1903.03016","evidence_quote":"Provides the SM expectation 3.92e-7 and the Belle-style Wilson-coefficient presentation this paper extends to the S1 leptoquark."},{"cited_title":"Sibidanov et al","cited_arxiv_id":null,"evidence_quote":"Reports the B→μν̄ measurement that motivates studying the decay."},{"cited_title":"Correlating Tauonic B Decays to the Neutron EDM via a Scalar Leptoquark","cited_arxiv_id":"1905.08257","evidence_quote":"Supplies the neutron-EDM effective Hamiltonian and S1 loop contributions used to constrain complex phases and motivate setting first-generation couplings to zero."},{"cited_title":"B¨ uchmuller, R","cited_arxiv_id":null,"evidence_quote":"Provides the S1 pair-production exclusion that cuts the allowed y^R_{12}-y^L_{32} plane."},{"cited_title":"Hiller, D","cited_arxiv_id":null,"evidence_quote":"Electron EDM limit that motivates setting the first-generation couplings to zero."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"FLAG value f_B = 190 MeV fixing the SM B→μν̄ normalization."},{"cited_title":"Raj, Phys","cited_arxiv_id":null,"evidence_quote":"Ditau resonance search that provides the 2σ exclusion on |y^R_{13}| shown in the τν̄ plots."}],"review_version":1}