REVIEW 4 major objections 4 minor 34 references
Angular observables and branching ratio for $B_s\to \phi \ell^+ \ell^-$ decay
T0 review · 4 major / 4 minor · reviewed 2026-08-04 · deepseek-v4-flash
Pith's one-line read The covariant confined quark model reproduces LHCb angular data for B_s→φμ+μ− once two-loop corrections to the Wilson coefficients are included.
desk verdict A useful CCQM cross-check undermined by a corrupted comparison table and an overclaimed agreement sentence. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The central object is the set of dimensionless form factors A0, A+, A−, V, a0, a+, g for the B_s→φ transition, parametrized by a dipole form F(q²)=F(0)/(1−a s + b s²) with s=q²/m_{Bs}². These feed the helicity amplitudes that define the observables. The second key ingredient is the effective Hamiltonian with Wilson coefficients, where the long-distance charm-loop effects and two-loop QCD corrections are absorbed into C7_eff and C9_eff; the two-loop corrections are included only in the q² ranges where their expansions are valid, 1.1–5.5 and 8.8–22 GeV². The combination of these ingredients produces predictions that match LHCb data.
What would settle it
A determination of B_s→φμ+μ− in the q² bins [5,8] and [6,8] GeV² with the same two-loop accuracy (or an experimental measurement with uncertainty below ~5%) would test the one-loop CCQM predictions directly; if the measured branching fraction in these bins deviates by more than the model's ~20% uncertainty while the two-loop-valid bins agree, the assumed treatment of the charm region would be falsified. Alternatively, a lattice QCD calculation of the B_s→φ form factors at q² = 0 that yields A0(0)=0.40 and V(0)=0.31 within errors would support the CCQM input.
Extended reading notes
Core claim
The paper computes the B_s→φ transition form factors within the CCQM and uses them together with the SM effective Hamiltonian to predict the differential branching ratio and angular observables for B_s→φℓ+ℓ−. It shows that after including the NNLL two-loop corrections to C7_eff and C9_eff, the predicted branching fraction B(B_s→φμ+μ−) and the longitudinal polarization F_L fall within the 2021 LHCb measurements in the bins [11,12.5], [15,17], [17,18.9], and [15,18.9] GeV², improving on the agreement with the earlier 2015 data. The paper concludes that the CCQM provides a consistent Standard Model description of this decay.
Load-bearing premise
The numerical predictions inherit the CCQM form factors and model parameters from an earlier paper without recomputation, so the agreement with LHCb data is conditional on those priors being correct; additionally, the long-distance charm-loop effects are modeled by a Breit-Wigner ansatz rather than computed from first principles.
Editorial extensions
If this is right
- If the CCQM description holds, the measured B_s→φμ+μ− observables are consistent with the Standard Model, and deviations seen in earlier bins may be attributed to missing higher-order corrections rather than new physics.
- The inclusion of NNLL corrections is necessary for agreement; one-loop predictions are insufficient, especially at low q².
- The optimized observable S7, which vanishes at leading order, becomes non-zero with the two-loop corrections, providing a testable SM prediction.
- For the bins [5,8] and [6,8] GeV², where two-loop corrections are not valid, the one-loop CCQM predictions are the best available and can be improved when full NNLO results for the charm region become available.
- The same framework yields predictions for the τ+τ− mode and B_s→φνν̄, which can be checked by future experiments.
Reading between the lines
- The paper's reliance on form factors from a previous calculation without recomputation means the numerical agreement is conditional on those model parameters; a public release or independent recalculation of the CCQM parameters would allow direct testing.
- Because the two-loop corrections are valid only in specific q² ranges, the 'agreement' claim is strongest in those bins; future high-precision data in the [5,8] and [6,8] bins would sharply discriminate the CCQM's one-loop treatment from other models.
- The near-equality A0(0)≈V(0) noted in the heavy-quark limit is not exact in the CCQM, which explains the smallness of P1; testing whether P1 remains small at higher q² in other models could probe the model's form-factor structure.
- If the charm-loop Breit-Wigner ansatz is inadequate, the agreement could be accidental; lattice QCD form factors combined with the same effective Hamiltonian would provide a cross-check.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. This paper presents a calculation of the B_s→φ transition form factors in the covariant confined quark model (CCQM), and uses them to predict branching fractions and angular observables (A_FB, F_L, A_i, S_i) for B_s→φℓ^+ℓ^- in various q² bins. The results are compared with LHCb data from 2015 and 2021 and with PQCD predictions from ref. [19]. The central claim is that the CCQM predictions agree with the latest LHCb data, particularly in specific high-q² bins, and that NNLL two-loop corrections to C7_eff and C9_eff are essential to this agreement. The manuscript includes tables of binned predictions, figures for differential distributions, and a comparison of integrated observables.
Significance. If the numerical results are correct, this would be a useful independent CCQM benchmark for b→sℓℓ transitions and a cross-check of the LHCb angular measurements. The paper's use of externally computed Wilson coefficients and published two-loop corrections is a strength, as is the direct comparison with two LHCb datasets. However, the main comparison table is substantially misaligned and the text overstates the level of agreement. The quantitative conclusions cannot be assessed until the table is corrected and the internal inconsistencies are resolved.
major comments (4)
- [Table VI] Table VI, the sole LHCb-comparison table, contains multiple row/column misalignments that invalidate the printed agreement claim. In the AFB block, row [1,6] lists CCQM = 0.69 ± 0.14 and PQCD = 0.777, identical to the FL values for the same bin, whereas the adjacent [1.1,6] AFB row is 0.034 ± 0.006. In the same AFB block, row [15,18.9] lists 0.394 ± 0.003, again an FL value. A5 [0.1,0.98] = 0.268 ± 0.054 while A5 [0.1,2] = 0.0031 ± 0.0006, an implausible jump for a smooth observable. A8 [1,6] = 0.17 ± 0.03 equals the S4 entry. Many S7 rows lack CCQM entries and instead list PQCD values around 0.45–0.99, inconsistent with the stated magnitude of S7. These errors must be corrected before the agreement claim can be evaluated.
- [Section V] The concluding sentence in Section V states that the integrated observables are in 'complete agreement' with PQCD [19]. This is contradicted by Table V: for the τ mode, <F_L> is 0.090 ± 0.02 (CCQM) versus 0.396+0.002/−0.003 (PQCD), a factor-of-4 difference. The μ-mode values are close, but the τ-mode F_L is not in agreement. This overstatement should be removed and the discrepancy discussed.
- [Section IV, Eq. (20)] The paper states that the [5,8] and [6,8] GeV² bins are omitted from NNLL corrections and therefore use one-loop Wilson coefficients, yet these bins are included in Table VI without any additional systematic uncertainty. Since the text credits NNLL corrections for the agreement with LHCb, the one-loop bins are not on the same footing as the rest of the table. The authors should either assign an uncertainty for the missing two-loop piece in these bins or explicitly exclude them from the agreement claim.
- [Eq. (20)] The manuscript gives two conflicting validity ranges for the low-q² NNLL corrections. Equation (20) states 1.1 ≤ q² ≤ 5.5 GeV², while the following paragraph says the low-q² two-loop results are reliable for 0.1 ≤ q² ≤ 6 GeV². Several bins in Table VI start at 0.1 or 0.98 GeV². Because the paper argues that NNLL corrections are essential to the agreement, the reader must know which bins actually received these corrections. Please reconcile the ranges and state the treatment of each bin.
minor comments (4)
- [Figures 1–6] The figure labels contain LaTeX artifacts such as '/LParen1GeV2/RParen1' and '/MiΝus'. The sentence before Fig. 1, 'The behavior of the differential branching B(Bs→φνν) is shown in Figs. 4', is incomplete.
- [Tables I–II] Table I gives A0(0)=0.40, while Table II lists the CCQM A0(0)=0.28±0.03. The relation between these quantities (via the BSW form-factor convention from ref. [5]) should be stated explicitly in the caption of Table II, otherwise the reader cannot reconcile the two tables.
- [Section IV] The text says 'the optimized observable P1 remains small across a large range of values' and that A0(0)−V(0)=0.09 leads to a 'truly small value of P1'. However, Table V reports <P1> = −0.52 ± 0.1 for the muon mode. This is not small in the conventional normalization of P1 and should be clarified.
- [Table IV] The branching ratios in Table IV are quoted without explicit reference to the q² integration range. For Bs→φµ+µ−, the total rate should be specified as the full q² range; otherwise the comparison with LHCb and other models is ambiguous.
Circularity Check
No significant circularity: the CCQM predictions are benchmarked against LHCb data using independent external inputs, not fitted to the comparison data.
full rationale
The derivation chain is not circular in any constructional sense. The form factors in Table I are explicitly inherited from the earlier CCQM paper [5] ('same as in our previous work'), but that is prior independent model output, not a quantity fitted to the LHCb data being compared; the same holds for the Wilson coefficients taken from [30], the two-loop corrections from [28,29], and the Breit-Wigner charm-loop ansatz from [27]. The dipole parametrization in Eq. (4) is a fit to the model's own numerical form factors with a stated <1% relative error, not a fit to the experimental target, so it does not turn the prediction into an input. The comparison with LHCb [1,2] and PQCD [19] is an external benchmark. The paper also honestly flags its own limitation that the [5,8] and [6,8] GeV^2 bins are one-loop because the two-loop validity intervals exclude them. The suspected Table VI row/column misalignments are a presentation and reproducibility concern about the claimed agreement, but they do not amount to a derivation that reduces by construction to its own inputs. The self-citation [5] is load-bearing for the numerical estimates but qualifies as independent prior support rather than circularity because its parameters were not tuned to the LHCb observables with which the paper compares.
Assumptions & free parameters
free parameters (3)
- CCQM model parameters (constituent quark masses, quark-meson couplings g_Bs, g_φ, vertex width, infrared cutoff) =
not given; fixed in [5]
- Dipole fit parameters F(0), a, b for each of the seven form factors =
Table I
- Breit-Wigner parameters for c-cbar resonances absorbed in C9_eff =
not specified in this paper
assumptions (4)
- domain assumption The Standard Model effective Hamiltonian (Eq. 5) with the ten-operator basis and Wilson coefficients taken from [30] is valid at the hadronic scale.
- domain assumption The CCQM confining vertex ansatz in Eqs. (2)-(3) reliably approximates the QCD matrix elements for the B_s→φ transition.
- domain assumption Charm-loop long-distance effects are represented by C9_eff with the Breit-Wigner ansatz of [27] plus the two-loop corrections of [28,29], and bins [5,8]/[6,8] may use one-loop coefficients.
- standard math The angular-observable definitions and helicity amplitudes from [31] and the optimized-observable relations from [33] are the correct ones for the B_s→φ(→K+K−)ℓ+ℓ− cascade.
Cite this review
Pith. "Pith review of Angular observables and branching ratio for $B_s\to \phi \ell^+ \ell^-$ decay." pith.science (2026). https://pith.science/paper/M63LVINN
@misc{pith2026251003739,
author = {Pith},
title = {Pith review of: Angular observables and branching ratio for $B_s\to \phi \ell^+ \ell^-$ decay},
year = {2026},
howpublished = {\url{https://pith.science/paper/M63LVINN}},
note = {Machine review of arXiv:2510.03739}
}
abstract
In this paper, an analysis of the $B_s\to \phi \ell^+ \ell^-$ rare decay is presented within the framework of the covariant confined quark model. The $B_s\to \phi$ transition form factors are calculated and then used to compute the branching fractions and angular observables in various $q^2$ bins, including the forward-backward asymmetry $A_{FB}$, the longitudinal polarization $F_L$, and the optimized observables $A_i$ and $S_i$. The results show agreement with the latest experimental data given by LHCb collaboration and compared with available theoretical predictions.
Figures
Figures from the paper (3 more)
Reference graph
Works this paper leans on
-
[19]
R. R. Horgan, Z. Liu, S. Meinel and M. Wingate, PoS LA TTICE2014 (2015), 372 doi:10.22323/1.214.0372 [arXiv:1501.00367 [hep-lat]]
arXiv 2015
-
[1]
R. Aaij et al. [LHCb], JHEP 09 (2015), 179 doi:10.1007/JHEP09(2015)179 [arXiv:1506.08 777 [hep-ex]]
-
[2]
R. Aaij et al. [LHCb], Phys. Rev. Lett. 127 (2021) no.15, 151801 doi:10.1103/PhysRevLett.127.151801 [arXiv:2105.14007 [hep-ex]]
arXiv 2021
-
[3]
T. Aaltonen et al. [CDF], Phys. Rev. Lett. 106 (2011), 161801 doi:10.1103/PhysRevLett.106.161801 [arXiv:1101.1028 [ hep-ex]]
arXiv 2011
-
[4]
R. Aaij et al. [LHCb], JHEP 07 (2013), 084 doi:10.1007/JHEP07(2013)084 [arXiv:1305.21 68 [hep-ex]]
-
[5]
S. Dubniˇ cka, A. Z. Dubniˇ ckov´ a, A. Issadykov, M. A. Ivanov, A. Liptaj and S. K. Sakhiyev, Phys. Rev. D 93 (2016) no.9, 094022 doi:10.1103/PhysRevD.93.094022 [arX iv:1602.07864 [hep-ph]]
arXiv 2016
-
[6]
0 5 10 15 /MiΝus 0.4 /MiΝus 0.2 0.0 0.2 q2 /LParen1GeV2/RParen1 A5 A6 A8 A9 0 5 10 15 /MiΝus 0.3 /MiΝus 0.2 /MiΝus 0.1 0.0 0.1 0.2 0.3 q2 /LParen1GeV2/RParen1 S3 S4 S7 FIG
(19) The q2-dependence of the Ai and Si for Bs →φµ+µ− decay is displayed in Fig.5. 0 5 10 15 /MiΝus 0.4 /MiΝus 0.2 0.0 0.2 q2 /LParen1GeV2/RParen1 A5 A6 A8 A9 0 5 10 15 /MiΝus 0.3 /MiΝus 0.2 /MiΝus 0.1 0.0 0.1 0.2 0.3 q2 /LParen1GeV2/RParen1 S3 S4 S7 FIG. 5: Observables Ai and Si for Bs →φµ+µ− decay. The q2-dependence of the Ai and Si for Bs →φτ +τ − deca...
2015
-
[7]
C. Q. Geng and C. C. Liu, J. Phys. G 29 (2003), 1103-1118 doi:10.1088/0954-3899/29/6/313 [arXiv:hep-ph/0303246 [hep-ph]]
arXiv 2003
Show all 34 references
-
[8]
Deandrea and A
A. Deandrea and A. D. Polosa, Phys. Rev. D 64 (2001), 074012 doi:10.1103/PhysRevD.64.074012 [arXiv:hep-ph/0105058 [hep-ph]]
2001 arXiv
-
[9]
Bobeth, G
C. Bobeth, G. Hiller and G. Piranishvili, JHEP 07 (2008), 106 doi:10.1088/1126- 6708/2008/07/106 [arXiv:0805.2525 [hep-ph]]
2008 arXiv
-
[10]
Ball and R
P. Ball and R. Zwicky, Phys. Rev. D 71 (2005), 014029 doi:10.1103/PhysRevD.71.014029 [arXiv:hep-ph/0412079 [hep-ph]]
2005 arXiv
-
[11]
Y. L. Wu, M. Zhong and Y. B. Zuo, Int. J. Mod. Phys. A 21 (2006), 6125-6172 doi:10.1142/S0217751X06033209 [arXiv:hep-ph/0604007 [ hep-ph]]
2006 arXiv
-
[12]
Altmannshofer and D
W. Altmannshofer and D. M. Straub, Eur. Phys. J. C 75 (2015) no.8, 382 doi:10.1140/epjc/s10052-015-3602-7 [arXiv:1411.3161 [ hep-ph]]
2015 arXiv
-
[13]
Bharucha, D
A. Bharucha, D. M. Straub and R. Zwicky, JHEP 08 (2016), 098 doi:10.1007/JHEP08(2016)098 [arXiv:1503.05534 [hep-ph ]]
2016 arXiv
-
[14]
Descotes-Genon, L
S. Descotes-Genon, L. Hofer, J. Matias and J. Virto, JHE P 06 (2016), 092 doi:10.1007/JHEP06(2016)092 [arXiv:1510.04239 [hep-ph ]]
2016 arXiv
-
[15]
Li and J
Y. Li and J. Hua, Eur. Phys. J. C 71 (2011), 1764 doi:10.1140/epjc/s10052-011-1764-5 [arXiv:1105.3031 [hep-ph]]. 19
2011 arXiv
-
[16]
Y. G. Xu, L. H. Zhou, B. Z. Li and R. M. Wang, Chin. Phys. C 37 (2013), 063104 doi:10.1088/1674-1137/37/6/063104 [arXiv:1305.5010 [h ep-ph]]
2013 arXiv
-
[17]
R. R. Horgan, Z. Liu, S. Meinel and M. Wingate, Phys. Rev. D 89 (2014) no.9, 094501 doi:10.1103/PhysRevD.89.094501 [arXiv:1310.3722 [hep- lat]]
2014 arXiv
-
[18]
R. R. Horgan, Z. Liu, S. Meinel and M. Wingate, Phys. Rev. Lett. 112 (2014), 212003 doi:10.1103/PhysRevLett.112.212003 [arXiv:1310.3887 [ hep-ph]]
2014 arXiv
-
[20]
S. P. Jin and Z. J. Xiao, Adv. High Energy Phys. 2021 (2021), 3840623 doi:10.1155/2021/3840623 [arXiv:2011.11409 [hep-ph]]
2021 arXiv
-
[21]
M. A. Ivanov, J. G. Korner, S. G. Kovalenko, P. Santorell i and G. G. Saidullaeva, Phys. Rev. D 85 (2012), 034004 doi:10.1103/PhysRevD.85.034004 [arXiv:1 112.3536 [hep-ph]]
2012 doi
-
[22]
R. N. Faustov and V. O. Galkin, Eur. Phys. J. C 73 (2013) no.10, 2593 doi:10.1140/epjc/s10052-013-2593-5 [arXiv:1309.2160 [ hep-ph]]
2013 arXiv
-
[23]
A. Ali, G. Kramer, Y. Li, C. D. Lu, Y. L. Shen, W. Wang and Y. M. Wang, Phys. Rev. D 76 (2007), 074018 doi:10.1103/PhysRevD.76.074018 [arXiv:h ep-ph/0703162 [hep-ph]]
2007
-
[24]
Melikhov and B
D. Melikhov and B. Stech, Phys. Rev. D 62 (2000), 014006 doi:10.1103/PhysRevD.62.014006 [arXiv:hep-ph/0001113 [hep-ph]]
2000 arXiv
-
[25]
R. H. Li, C. D. Lu and W. Wang, Phys. Rev. D 79 (2009), 034014 doi:10.1103/PhysRevD.79.034014 [arXiv:0901.0307 [hep- ph]]
2009 arXiv
-
[26]
C. D. Lu, W. Wang and Z. T. Wei, Phys. Rev. D 76 (2007), 014013 doi:10.1103/PhysRevD.76.014013 [arXiv:hep-ph/0701265 [hep-ph]]
2007 arXiv
-
[27]
Buchalla, A
G. Buchalla, A. J. Buras and M. E. Lautenbacher, Rev. Mod . Phys. 68 (1996), 1125-1144 doi:10.1103/RevModPhys.68.1125 [arXiv:hep-ph/9512380 [hep-ph]]
1996 arXiv
-
[28]
A. Ali, T. Mannel and T. Morozumi, Phys. Lett. B 273 (1991), 505-512 doi:10.1016/0370- 2693(91)90306-B
1991 doi
-
[29]
H. H. Asatryan, H. M. Asatrian, C. Greub and M. Walker, Ph ys. Rev. D 65 (2002), 074004 doi:10.1103/PhysRevD.65.074004 [arXiv:hep-ph/0109140 [hep-ph]]
2002 arXiv
-
[30]
Greub, V
C. Greub, V. Pilipp and C. Schupbach, JHEP 12 (2008), 040 doi:10.1088/1126- 6708/2008/12/040 [arXiv:0810.4077 [hep-ph]]
2008 arXiv
-
[31]
Descotes-Genon, T
S. Descotes-Genon, T. Hurth, J. Matias and J. Virto, JHE P 05 (2013), 137 20 doi:10.1007/JHEP05(2013)137 [arXiv:1303.5794 [hep-ph] ]
2013 arXiv
-
[32]
Faessler, T
A. Faessler, T. Gutsche, M. A. Ivanov, J. G. Korner and V. E. Lyubovitskij, Eur. Phys. J. direct 4 (2002) no.1, 18 doi:10.1007/s1010502c0018 [arXiv:hep-ph /0205287 [hep-ph]]
2002 doi
-
[33]
K. A. Olive et al. [Particle Data Group], Chin. Phys. C 38 (2014), 090001 doi:10.1088/1674- 1137/38/9/090001
2014 doi
-
[34]
Descotes-Genon, J
S. Descotes-Genon, J. Matias, M. Ramon and J. Virto, JHE P 01 (2013), 048 doi:10.1007/JHEP01(2013)048 [arXiv:1207.2753 [hep-ph] ]. 21
2013 arXiv
Reviewed August 4, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.