{"id":"0829516b-9d12-421c-a456-e86035ff30c0","arxiv_id":"2502.09883","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"Nonfactorizable one-loop QED corrections enhance Bs to Ds(*) l nu branching ratios and reduce the ratios R(Ds(*)), especially for the electron mode.","lead":"This paper calculates how QED corrections that do not factorize into separate hadron and lepton pieces change the predicted rates of Bs meson decays to Ds mesons plus a lepton and a neutrino. The corrections raise the predicted branching ratios and lower the lepton-universality ratios R(Ds(*)), but the underlying form-factor uncertainties remain very large.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Virtual-only QED corrections contain uncanceled double logarithms; the large electron-mode enhancement and R(Ds) reduction are not physical without real photon emission.","rationale":"The reader’s weakest_assumption is the neglect of spectator scattering, which the paper explicitly acknowledges. That is a valid concern, but it is not the single most load-bearing issue. The omission of real photon emission is more fundamental: it affects the infrared and collinear safety of the entire calculation. The large electron-mode effect—a 24% branching-ratio enhancement—has the characteristic size of double Sudakov logarithms α/π ln^2(mℓ/mB) present in the virtual correction but cancelled by soft/collinear real emission. Without including b→cℓνγ and specifying the photon energy cut, the predicted branching ratios and R(Ds(*)) have no clear physical meaning. This is not a small correction that could be added later; it is an essential part of the one-loop QED correction. The paper’s summary even states that verification is impracticable, but more importantly the calculation does not yet define an observable. The reader’s conditional acceptance was based on mending the spectator-scattering approximation and error budget; our concern requires a substantially different calculation. Hence the verdict should move to REJECT: the central claim as stated is not supported by the presented calculation.","tokens_in":12345,"tokens_out":12955,"duration_ms":145099,"concrete_test":"Compute the full O(α) decay rate for Bs→Ds e−ν̄e by adding the real-photon bremsstrahlung diagram b→c e−ν̄eγ to the virtual corrections of Eq. (12), using a realistic photon-energy resolution (e.g., ΔE = 20 MeV in the Bs rest frame), and recompute the branching ratios and R(Ds)e. If the electron-mode branching ratio changes by more than 10% relative to the virtual-only value, or if R(Ds)e moves significantly back toward the ηEW value, then the claimed enhancement and reduction are artifacts of the omitted real-photon contribution.","verdict_should_be":"REJECT","load_bearing_attack":"The central numerical claim rests on the factor ˜ηEW = 1 + α(ηb+ηc) in Eq. (12), whose ηb,c are the QED vertex corrections of Fig. 1(b),(c) only. No real-photon bremsstrahlung diagrams are introduced anywhere in the paper, yet physical branching ratios require the sum of virtual and real photon contributions to be infrared- and collinear-finite. The expressions in Eqs. (13)–(14) contain double logarithms of the charged-lepton mass, e.g., terms like ln(tb) ln((tb−sb)/(1−tb)) and analogous Dilog terms; these are precisely the contributions that are cancelled or converted to single logarithms with a photon-energy cut when b→cℓνγ is added. The paper’s numerical results show a 24% enhancement of the electron-mode branching ratio (2.23%→2.77% in Table II) and a 20% reduction of R(Ds)e (0.298→0.240 in Table III), effects that are driven by these uncanceled double logs. Because no photon-energy acceptance is specified and no real-emission contribution is included, the branching ratios and ratios in Tables II and III are not well-defined physical observables, and the abstract’s claim that QED contributions enhance branching ratios and reduce R(Ds(*)) is not established.","agreement_with_reader":"disagree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This manuscript reevaluates the semileptonic decays \\bar{B}_s \\to D_s^{(*)} \\ell \\bar{\\nu}_\\ell including a one-loop QED vertex correction. The authors take the analytic expressions for the correction factors \\eta_b and \\eta_c from their earlier paper (Ref. [6]), replace the universal short-distance factor \\eta_EW in Eq. (11) with the lepton-flavor-dependent factor \\tilde{\\eta}_EW = 1 + \\alpha(\\eta_b+\\eta_c), and combine this with HPQCD lattice form factors and PDG inputs to compute branching ratios, R(D_s^{(*)}), helicity fractions, and differential distributions. The main numerical claims are that the QED factor enhances the electron-mode branching ratio by about 24% (2.23% to 2.77% in Table II), leaves the tau mode almost unchanged, and reduces R(D_s)_e and R(D_s^*)_e by about 20% (Table III), with much smaller shifts in the muon modes.","tokens_in":12549,"tokens_out":6465,"duration_ms":72710,"significance":"If the calculation defined a physical O(\\alpha) observable, the result would be relevant to lepton-flavor-universality tests at LHCb, Belle II, and future Z factories, and would provide an interesting comparison with the B \\to D^{(*)} analysis of Ref. [6]. The paper has clear strengths: it uses published lattice form factors, quotes and propagates form-factor uncertainties, introduces no free parameters, and presents useful tables and figures including a comparison with existing LHCb data. The SU(3) comparison in Fig. 3 is a valuable consistency check. However, the central quantitative claim is currently tied to a set of virtual corrections that are not made infrared- and collinear-safe by a corresponding real-photon treatment, so the physical interpretation of the numbers in Tables II and III is not established as it stands.","major_comments":[{"comment":"The calculation contains no real-photon bremsstrahlung contribution and no specification of a photon-energy acceptance. The analytic expressions for \\eta_b and \\eta_c contain double-logarithmic and Dilogarithmic functions of the charged-lepton mass, e.g. \\ln(t_b)\\ln((t_b-s_b)/(1-t_b)) and Li_2 terms, which in QED must be combined with real emission (or a photon-energy cut) to produce an infrared- and collinear-finite observable. The statement that ultraviolet and infrared divergences have been subtracted does not remove this problem: after such a subtraction the residual numerical value depends on the subtraction scheme unless a physical observable is defined. Consequently, the branching ratios in Table II and the ratios R(D_s^{(*)}) in Table III are not well-defined physical observables, and the abstract's claim that the QED contributions enhance branching ratios and reduce R(D_s^{(*)}) is not established by the calculation as presented. The revision should include a full O(\\alpha) treatment with real photon emission, or should reframe the quantity as a scheme-dependent ingredient of such a treatment.","section":"Section II, Eqs. (12)-(14); Section III, Tables II and III"},{"comment":"The spectator-scattering QED corrections, in which the photon couples to the spectator s quark and the charged lepton, are discarded with the statement that they are left out 'at a first approximation for the time being.' No quantitative estimate or symmetry argument is given. These are O(\\alpha) corrections of the same type as the ones computed, and their lepton-mass dependence could in principle alter the flavor ordering responsible for the R(D_s^{(*)}) shifts. The manuscript needs at least a power-counting estimate, a numerical bound, or a demonstration that such contributions cancel for the B_s system before the central claim can be regarded as a complete QED result.","section":"Section II, paragraph after Eq. (11)"},{"comment":"The analytic core of the correction is imported from the authors' own Ref. [6] without derivation or an independent cross-check within this manuscript. Since all quantitative conclusions flow from the functions \\eta_b and \\eta_c, the revision should either reproduce the derivation (for example, in an appendix) or provide a check against a known limit or against a numerical evaluation, so that the sign and magnitude of the double-logarithmic terms can be assessed independently of the infrared-safety concern raised above.","section":"Section II, Eqs. (13)-(14)"}],"minor_comments":[{"comment":"The decay-rate formula in Eq. (19) is written with the universal factor |\\eta_EW|^2, but the numerical results in Tables II and III use the lepton-flavor-dependent \\tilde{\\eta}_EW. Unless the symbol in Eq. (19) is intended generically, this is an inconsistency that should be corrected.","section":"Eq. (19)"},{"comment":"The variables s_b, t_b, s_c, and t_c are introduced only through the relations in Eqs. (15)-(18). An explicit definition at first use would improve readability and help the reader check the arguments of the logarithms and dilogarithms.","section":"Eqs. (13)-(18)"},{"comment":"There are several typographical artifacts, including 'bran ching' in the abstract and 'Tara-Z' in the introduction, which should presumably read 'Tera-Z'.","section":"Introduction and Abstract"},{"comment":"The quark-mass approximations m_b \\approx m_{B_s} and m_c \\approx m_{D_s^{(*)}}, together with \\mu_MS = m_b, should be accompanied by an estimate of the induced uncertainty, because the renormalization logarithms in Eqs. (13)-(14) depend on this choice.","section":"Section II, after Eq. (14)"},{"comment":"The axis labels in Figs. 2 and 5 appear garbled in the manuscript text; these figures should be regenerated with correct notation.","section":"Figures 2 and 5"}],"recommendation":"major_revision","confidential_remarks":null},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The short version: this is a clean, transparent extension of the authors' earlier QED vertex-correction calculation to the Bs -> Ds(*) l nu family, and the numerics are internally consistent. But the headline claim is not physically well-defined because the calculation includes only virtual photon corrections. The formulas in Eqs. (13)-(14) contain double logarithms of the charged-lepton mass (e.g., ln(t_b) ln((t_b - s_b)/(1 - t_b))), which are mass singularities that cancel only when real photon emission is added with a specified photon-energy cut. No real-emission diagram appears anywhere in the paper. So the 24% enhancement of the electron-mode branching ratio and the 20% reduction of R(D_s)_e in Tables II and III are not predictions for any measurable branching ratio or ratio. They are artifacts of an infrared-incomplete calculation.\n\nThe paper does several things well. It takes HPQCD lattice form factors for both Bs -> Ds and Bs -> Ds*, propagates form-factor uncertainties, and is honest about the large errors on the Ds* modes. The SU(3) comparison with B -> D(*) in Fig. 3 is a nice cross-check. The paper also states explicitly that verification is impracticable at present, which is commendable.\n\nThe soft spots beyond the missing real emission: the QED correction factor is imported from the authors' own Ref. [6] without derivation, which is acceptable as self-citation but means the reader must trust that earlier paper. Spectator scattering is neglected without a quantitative estimate; the paper admits this. And the Ds* form-factor uncertainties are so large that the claimed R(Ds*) reduction is not statistically significant, as the authors note.\n\nThe missing real-photon piece is the load-bearing flaw. It is fixable, and the fix is standard: add b -> c l nu gamma with a photon veto or energy cut, compute the infrared-finite combination, and present results for a specified cut. Until then, the numerical predictions are not observables.\n\nFor a reader working on semileptonic B decays and radiative corrections, this paper is worth engaging with as a draft, but not worth citing yet. I would send it to a competent referee, because the flaw is subtle and fixable and the topic is relevant for LHCb and Belle II. The authors are thinking seriously about a real issue, but the central claim does not hold as written.","headline":"Transparent extension to Bs decays, but the virtual-only QED correction leaves uncanceled lepton-mass logarithms, so the headline effect is not a physical observable yet.","tokens_in":13137,"tokens_out":2838,"would_cite":false,"duration_ms":28508,"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":"Nonfactorizable one-loop QED corrections enhance $\\overline{B}_s \\to D_s^{(*)} \\ell \\bar{\\nu}_\\ell$ branching ratios and reduce $R(D_s^{(*)})$ in a lepton-flavor-dependent way.","keywords":["semileptonic B_s decays","nonfactorizable QED corrections","lepton flavor universality","R(D_s) ratios","branching ratios","lattice QCD form factors","CKM matrix element V_cb"],"falsifier":"Measure the ratio $\\mathcal{B}(\\overline{B}_s\\to D_s e\\bar{\\nu}_e)/\\mathcal{B}(\\overline{B}_s\\to D_s\\mu\\bar{\\nu}_\\mu)$ precisely; the paper's $\\tilde{\\eta}_{EW}$ prediction is about 1.20, while the lepton-flavor-universal $\\eta_{EW}$ prediction is about 1.00. A value near 1.00 with small uncertainty would falsify the claimed nonfactorizable QED enhancement.","tokens_in":12073,"feed_emoji":"⚛️","tokens_out":8945,"duration_ms":76757,"temperature":0.7,"pith_summary":"This paper tries to establish that virtual photon exchange between the $b$ or $c$ quark and the charged lepton in $\\overline{B}_s \\to D_s^{(*)} \\ell \\bar{\\nu}_\\ell$ decays is not a flavor-blind correction. The one-loop QED vertex terms are collected into a factor $\\tilde{\\eta}_{EW}$ that depends on the lepton mass, so the branching ratios and their ratios $R(D_s^{(*)})$ shift by different amounts for $e$, $\\mu$, and $\\tau$. The paper finds these nonfactorizable corrections raise the branching ratios and lower $R(D_s^{(*)})$, most strongly for electrons, and that the semimuonic rates agree better with measured data. It also claims the SU(3) flavor symmetry is well preserved across the $R(D)$--$R(D^*)$ plane for the charmed semileptonic $B_{u,d,s}$ decays. The stakes are whether a lepton-flavor-dependent Standard Model effect of this kind can be identified or must be mimicked by new physics.","feed_headline":"QED corrections shrink R(Ds) ratios and boost Bs decay rates","feed_subtitle":"One-loop photon exchange between quarks and leptons makes the rates lepton-dependent and pulls the electron ratios down.","key_machinery":"The central object is the one-loop QED correction factor $\\tilde{\\eta}_{EW} = 1 + \\alpha_{em}(\\eta_b + \\eta_c)$, where $\\eta_b$ and $\\eta_c$ are the analytic vertex-correction functions in Eqs. (13) and (14). These functions are built from kinematic ratios $s_b,t_b,s_c,t_c$ involving the quark and lepton momenta, dilogarithms, and logarithms of $m_\\ell/m_{b,c}$, and they carry the lepton-mass dependence that the standard $\\eta_{EW}$ lacks. The factor multiplies the leading-order decay amplitude before the helicity amplitudes and phase-space integration are performed, so it acts directly on the differential decay rate.","core_discovery":"Within the Standard Model, the one-loop QED correction factor $\\tilde{\\eta}_{EW}=1+\\alpha_{em}(\\eta_b+\\eta_c)$ replaces the universal short-distance factor $\\eta_{EW}$. The functions $\\eta_b$ and $\\eta_c$ encode photon exchange between the $b$/$c$ quarks and the charged lepton and depend on $q^2$, the lepton mass, and the scattering angle. Using lattice QCD form factors, the paper obtains branching ratios that grow relative to the $\\eta_{EW}$ results, while the ratios $R(D_s)_e$ and $R(D_s^*)_e$ drop from $0.298$ to $0.240$ and from $0.248$ to $0.199$, respectively; the corresponding muon ratios move only slightly. The authors read this as evidence that nonfactorizable QED corrections introduce a lepton-flavor dependence in the effective weak coupling, which can move the predicted ratios away from the measured $R(D)$--$R(D^*)$ region and sharpen, rather than resolve, the lepton flavor universality tension.","pith_inferences":["If the same nonfactorizable QED mechanism is extrapolated to $B\\to D^{(*)} \\ell\\bar{\\nu}_\\ell$ decays, the electron ratios there should also be suppressed relative to the universal-$\\eta_{EW}$ prediction, which would change how much room remains for new physics in the R(D) tension.","A direct test is the ratio $\\mathcal{B}(\\overline{B}_s\\to D_s e\\bar{\\nu}_e)/\\mathcal{B}(\\overline{B}_s\\to D_s\\mu\\bar{\\nu}_\\mu)$: the calculation predicts about 1.20, while the flavor-blind factor predicts about 1.00, so a precise measurement near 1.00 would rule out the claimed lepton-flavor dependence.","The spectator-scattering corrections the paper sets aside could be estimated with the same helicity machinery; if they are not small or are strongly flavor-dependent, the central pattern of enhancement and ratio reduction could change."],"forward_implications":["Branching ratios for $\\overline{B}_s\\to D_s^{(*)}e\\bar{\\nu}_e$ and $\\overline{B}_s\\to D_s^{(*)}\\mu\\bar{\\nu}_\\mu$ are predicted to increase relative to the universal-$\\eta_{EW}$ calculation, while the $\\tau$ channels barely move.","The ratios $R(D_s)_e$ and $R(D_s^*)_e$ are predicted to fall by roughly 19--20 percent, making the electron-versus-muon pattern of the ratios more pronounced.","The semimuonic branching ratios $\\mathcal{B}(\\overline{B}_s\\to D_s^{(*)} \\mu\\bar{\\nu}_\\mu)$ move closer to the available measurements.","The ratios $R(D)$--$R(D^*)$ for $B_{u,d,s}$ decays remain consistent with SU(3) flavor symmetry under the same lepton-flavor-dependent corrections.","For $\\overline{B}_s\\to D_s^*\\ell\\bar{\\nu}_\\ell$, current form-factor uncertainties are large enough to mask the QED effect and complicate $V_{cb}$ extraction."],"supporting_citations":[{"why":"Supplies the QED vertex-correction formalism and analytic expressions this paper extends from B to Bs decays.","marker":"[6]"},{"why":"Defines the short-distance electroweak correction factor $\\eta_{EW}$ that the paper compares against and replaces with $\\tilde{\\eta}_{EW}$.","marker":"[8]"},{"why":"Provides the $B_s\\to D_s$ lattice form factors used in the $B_s\\to D_s$ helicity amplitudes.","marker":"[11]"},{"why":"Provides the $B_s\\to D_s^*$ lattice form factors used in the $B_s\\to D_s^*$ helicity amplitudes.","marker":"[12]"},{"why":"Supplies the particle masses, lifetime, $|V_{cb}|$, and measured branching fractions used as inputs and comparisons.","marker":"[1]"},{"why":"Supplies the measured $B_s\\to D_s^{(*)} \\mu\\bar{\\nu}_\\mu$ branching fractions used to compare the semimuonic rates.","marker":"[3]"},{"why":"Supplies the estimated $R(D_s^*)_\\mu$ value quoted in Table III for comparison.","marker":"[13]"},{"why":"Provides the experimental $R(D)$ and $R(D^*)$ averages used in the SU(3) flavor-symmetry correlation plot.","marker":"[2]"}],"fun_headline_variants":["Loop QED lowers R(Ds) and lifts Bs branching ratios","Nonfactorizable QED cuts R(Ds) and boosts Bs rates","One-loop photon exchange shifts R(Ds) down in Bs decays","Lepton-dependent QED correction shrinks R(Ds)"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The calculation assumes that photon-exchange corrections involving the spectator $s$ quark are negligible; if they are not, the predicted lepton-flavor pattern of the QED effect could change.","fun_headline_variants_meta":{"raw":{"variants":["Loop QED lowers R(Ds) and lifts Bs branching ratios","Nonfactorizable QED cuts R(Ds) and boosts Bs rates","One-loop photon exchange shifts R(Ds) down in Bs decays","Lepton-dependent QED correction shrinks R(Ds)"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000597,"raw_usage":{"total_tokens":2801,"prompt_tokens":961,"completion_tokens":1840,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":577,"completion_tokens_details":{"reasoning_tokens":1763}},"tokens_in":577,"tokens_out":1840,"duration_ms":13698,"temperature":1.0,"reasoning_tokens":1763,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T20:09:51.507320+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the ratio $\\mathcal{B}(\\overline{B}_s\\to D_s e\\bar{\\nu}_e)/\\mathcal{B}(\\overline{B}_s\\to D_s\\mu\\bar{\\nu}_\\mu)$ precisely; the paper's $\\tilde{\\eta}_{EW}$ prediction is about 1.20, while the lepton-flavor-universal $\\eta_{EW}$ prediction is about 1.00. A value near 1.00 with small uncertainty would falsify the claimed nonfactorizable QED enhancement.","supporting_citations":[{"cited_title":"Yang , L","cited_arxiv_id":null,"evidence_quote":"Supplies the QED vertex-correction formalism and analytic expressions this paper extends from B to Bs decays."},{"cited_title":"Sirlin, Large mW , mZ behaviour of the O(α) corrections to semileptonic processes mediated by W , Nucl","cited_arxiv_id":null,"evidence_quote":"Defines the short-distance electroweak correction factor $\\eta_{EW}$ that the paper compares against and replaces with $\\tilde{\\eta}_{EW}$."},{"cited_title":"McLean, C","cited_arxiv_id":null,"evidence_quote":"Provides the $B_s\\to D_s$ lattice form factors used in the $B_s\\to D_s$ helicity amplitudes."},{"cited_title":"Harrison, C","cited_arxiv_id":null,"evidence_quote":"Provides the $B_s\\to D_s^*$ lattice form factors used in the $B_s\\to D_s^*$ helicity amplitudes."},{"cited_title":"Navas, C","cited_arxiv_id":null,"evidence_quote":"Supplies the particle masses, lifetime, $|V_{cb}|$, and measured branching fractions used as inputs and comparisons."},{"cited_title":"Aaij et al","cited_arxiv_id":null,"evidence_quote":"Supplies the measured $B_s\\to D_s^{(*)} \\mu\\bar{\\nu}_\\mu$ branching fractions used to compare the semimuonic rates."},{"cited_title":"Paolucci (LHCb Collaboration), Study of the measure ment of the ratio R(D∗ s ) at LHCb, Nuovo Cim","cited_arxiv_id":null,"evidence_quote":"Supplies the estimated $R(D_s^*)_\\mu$ value quoted in Table III for comparison."}],"review_version":1}