×
近期发现有不法分子冒充我刊与作者联系,借此进行欺诈等不法行为,请广大作者加以鉴别,如遇诈骗行为,请第一时间与我刊编辑部联系确认(《中国物理C》(英文)编辑部电话:010-88235947,010-88236950),并作报警处理。
本刊再次郑重声明:
(1)本刊官方网址为cpc.ihep.ac.cn和https://iopscience.iop.org/journal/1674-1137
(2)本刊采编系统作者中心是投稿的唯一路径,该系统为ScholarOne远程稿件采编系统,仅在本刊投稿网网址(https://mc03.manuscriptcentral.com/cpc)设有登录入口。本刊不接受其他方式的投稿,如打印稿投稿、E-mail信箱投稿等,若以此种方式接收投稿均为假冒。
(3)所有投稿均需经过严格的同行评议、编辑加工后方可发表,本刊不存在所谓的“编辑部内部征稿”。如果有人以“编辑部内部人员”名义帮助作者发稿,并收取发表费用,均为假冒。
                  
《中国物理C》(英文)编辑部
2024年10月30日

Discriminating 1D new physics solutions in bsℓℓ decays

  • The recent measurements of $ R_{K^+} $, $ R_{K_S^0} $, $ R_{K^{*+}} $, $ B_s\to\mu^+\mu^- $, a set of CP-averaged angular observables for the $ B^0\to K^{*0}\mu^+\mu^- $ decay and its isospin partner $ B^+\to K^{*+}\mu^+\mu^- $ by the LHCb Collaboration consistently hint at lepton universality violation in the $ b\to s\ell\ell $ transitions. In this work, we first perform global fits to the $ b\to s\ell\ell $ data and show that five one-dimensional scenarios, i.e, $ \delta C_9^{\mu} $, $ \delta C_{10}^{\mu} $, $ \delta C_L^{\mu} $, $ \delta C_9^{\mu}=C_{10}^{\mu\prime} $, and $ \delta C_9^{\mu}=-C_9^{\mu\prime} $ can best explain the so-called B anamolies. Furthermore, we explore how these scenarios can be distinguished from each other. For this purpose, we first study the combinations of four angular asymmetries $ A_i ~~(i=3,4,5,9) $ and find that they cannot distinguish the five new physics scenarios. We then show that a newly constructed ratio $ R_{S} $ can uniquely discriminate the five new physics scenarios in proper intervals of $ q^2 $ if it can be measured with percent-level precision.
  • 加载中
  • [1] S. L. Glashow, J. Iliopoulos, and L. Maiani, Phys. Rev. D 2, 1285 (1970) doi: 10.1103/PhysRevD.2.1285
    [2] Y. Li and C.-D. Lü, Sci. Bull. 63, 267 (2018), arXiv:1808.02990[hep-ph doi: 10.1016/j.scib.2018.02.003
    [3] S. Bifani, S. Descotes-Genon, A. Romero Vidal et al., J. Phys. G 46, 023001 (2019), arXiv:1809.06229[hep-ex
    [4] S. Chen, Y. Li, W. Qian et al., (2021), arXiv: 2111.14360[hep-ex]
    [5] P. Gambino et al., Eur. Phys. J. C 80, 966 (2020), arXiv:2006.07287[hep-ph
    [6] J. Albrecht, D. van Dyk, and C. Langenbruch, Prog. Part. Nucl. Phys. 120, 103885 (2021), arXiv:2107.04822[hep-ex doi: 10.1016/j.ppnp.2021.103885
    [7] D. London and J. Matias, (2021), 10.1146/annurev-nucl-102020-090209, arXiv: 2110.13270[hep-ph]
    [8] J. M. Camalich and M. Patel, Science Bulletin 67, 1 (2022) doi: 10.1016/j.scib.2021.09.012
    [9] R. Aaij et al., Nature Phys. 18(3), 277-282 (2022) doi: 10.1038/s41567-021-01478-8
    [10] L.-S. Geng, B. Grinstein, S. Jäger et al., Phys. Rev. D 96, 093006 (2017), arXiv:1704.05446[hep-ph
    [11] R. Aaij et al. (LHCb), Phys. Rev. Lett. 113, 151601 (2014), arXiv:1406.6482[hep-ex doi: 10.1103/PhysRevLett.113.151601
    [12] R. Aaij et al., Phys. Rev. D 105(1), 012010 (2022) doi: 10.1103/PhysRevD.105.012010
    [13] R. Aaij et al., Phys. Rev. Lett. 128(4), 041801 (2022) doi: 10.1103/PhysRevLett.128.041801
    [14] L.-S. Geng, B. Grinstein, S. Jäger et al., Phys. Rev. D 104, 035029 (2021), arXiv:2103.12738[hep-ph
    [15] M. Beneke, C. Bobeth, and R. Szafron, JHEP 10, 232 (2019), arXiv:1908.07011[hep-ph
    [16] A. M. Sirunyan et al. (CMS), JHEP 04, 188 (2020), arXiv:1910.12127[hep-ex
    [17] M. Aaboud et al. (ATLAS), JHEP 04, 098 (2019), arXiv:1812.03017[hep-ex
    [18] R. Aaij et al. (LHCb), (2021), arXiv: 2110.09501[hep-ex]
    [19] A. Abdesselam et al. (Belle), Phys. Rev. Lett. 126, 161801 (2021), arXiv:1904.02440[hep-ex doi: 10.1103/PhysRevLett.126.161801
    [20] S. Choudhury et al. (BELLE), JHEP 03, 105 (2021), arXiv:1908.01848[hep-ex
    [21] R. Aaij et al. (LHCb), Phys. Rev. Lett. 125, 011802 (2020), arXiv:2003.04831[hep-ex doi: 10.1103/PhysRevLett.125.011802
    [22] R. Aaij et al. (LHCb), Phys. Rev. Lett. 126, 161802 (2021), arXiv:2012.13241[hep-ex doi: 10.1103/PhysRevLett.126.161802
    [23] R. Aaij et al. (LHCb), JHEP 02, 104 (2016), arXiv:1512.04442[hep-ex
    [24] W. Altmannshofer and P. Stangl, Eur. Phys. J. C 81, 952 (2021), arXiv:2103.13370[hep-ph doi: 10.1140/epjc/s10052-021-09725-1
    [25] C. Cornella, D. A. Faroughy, J. Fuentes-Martin et al., JHEP 08, 050 (2021), arXiv:2103.16558[hep-ph
    [26] J. Kriewald, C. Hati, J. Orloff et al., in 55th Rencontres de Moriond on Electroweak Interactions and Unified Theorie (2021) arXiv: 2104.00015[hep-ph]
    [27] M. Algueró, B. Capdevila, S. Descotes-Genon et al., Eur. Phys. J. C 82(4), 326 (2022) doi: 10.1140/epjc/s10052-022-10231-1
    [28] T. Hurth, F. Mahmoudi, D. M. Santos et al., Phys. Lett. B 824, 136838 (2022) doi: 10.1016/j.physletb.2021.136838
    [29] A. Carvunis, F. Dettori, S. Gangal et al., JHEP 12, 078 (2021) doi: 10.1007/JHEP12(2021)078
    [30] A. Angelescu, D. Bečirević, D. A. Faroughy et al., Phys. Rev. D 104, 055017 (2021), arXiv:2103.12504[hep-ph doi: 10.1103/PhysRevD.104.055017
    [31] G. Isidori, D. Lancierini, P. Owen et al., Phys. Lett. B 822, 136644 (2021), arXiv:2104.05631[hep-ph doi: 10.1016/j.physletb.2021.136644
    [32] J. Alda, J. Guasch, and S. Penaranda, PoS EPS-HEP2021, 494 (2022) doi: 10.22323/1.398.0494
    [33] G. Isidori, D. Lancierini, A. Mathad et al., (2021), arXiv: 2110.09882[hepph]
    [34] J.-P. Lee, J. Korean Phys. Soc. 80(1), 13-19 (2022) doi: 10.1007/s40042-021-00350-5
    [35] R. Bause, H. Gisbert, M. Golz et al., JHEP 12, 061 (2021) doi: 10.1007/JHEP12(2021)061
    [36] M. Ciuchini, M. Fedele, E. Franco et al., (2021), arXiv: 2110.10126[hep-ph]
    [37] J.-Y. Cen, Y. Cheng, X.-G. He et al., Nucl. Phys. B 978, 115762 (2022) doi: 10.1016/j.nuclphysb.2022.115762
    [38] J. Kawamura and S. Raby, Phys. Rev. D 104, 035007 (2021), arXiv:2104.04461[hep-ph doi: 10.1103/PhysRevD.104.035007
    [39] B. C. Allanach, J. M. Butterworth, and T. Corbett, Eur. Phys. J. C 81(12), 1126 (2021) doi: 10.1140/epjc/s10052-021-09919-7
    [40] A. K. Alok, N. R. S. Chundawat, and D. Kumar, Eur. Phys. J. C 82(1), 30 (2022) doi: 10.1140/epjc/s10052-021-09975-z
    [41] P. Ko, T. Nomura, and H. Okada, (2021), arXiv: 2110.10513[hep-ph]
    [42] M. F. Navarro and S. F. King, Phys. Rev. D 105(3), 035015 (2022) doi: 10.1103/PhysRevD.105.035015
    [43] R. Bause, G. Hiller, T. Höhne et al., Eur. Phys. J. C 82(1), 42 (2022) doi: 10.1140/epjc/s10052-021-09957-1
    [44] A. K. Alok, A. Dighe, S. Gangal et al., (2021), arXiv: 2108.05614[hep-ph]
    [45] X. Wang, (2021), arXiv: 2108.01279[hep-ph]
    [46] J. Davighi, JHEP 08, 101 (2021), arXiv:2105.06918[hep-ph
    [47] A. Greljo, P. Stangl, and A. E. Thomsen, Phys. Lett. B 820, 136554 (2021), arXiv:2103.13991[hep-ph doi: 10.1016/j.physletb.2021.136554
    [48] T. Nomura and H. Okada, Phys. Rev. D 104, 035042 (2021), arXiv:2104.03248[hep-ph doi: 10.1103/PhysRevD.104.035042
    [49] M. Du, J. Liang, Z. Liu et al., (2021), arXiv: 2104.05685[hep-ph]
    [50] P. F. Perez, C. Murgui, and A. D. Plascencia, Phys. Rev. D 104, 035041 (2021), arXiv:2104.11229[hep-ph doi: 10.1103/PhysRevD.104.035041
    [51] A. Greljo, Y. Soreq, P. Stangl et al., JHEP 04, 151 (2022) doi: 10.1007/JHEP04(2022)151
    [52] S. F. King, JHEP 11, 161 (2021) doi: 10.1007/JHEP11(2021)161
    [53] S. Singirala, S. Sahoo, and R. Mohanta, Phys. Rev. D 105(1), 015033 (2022) doi: 10.1103/PhysRevD.105.015033
    [54] W.-F. Chang, JHEP 09, 043 (2021), arXiv:2105.06917[hep-ph
    [55] K. Ban, Y. Jho, Y. Kwon et al., (2021), arXiv: 2104.06656[hep-ph]
    [56] D. Marzocca and S. Trifinopoulos, Phys. Rev. Lett. 127, 061803 (2021), arXiv:2104.05730[hep-ph doi: 10.1103/PhysRevLett.127.061803
    [57] H. M. Lee, Phys. Rev. D 104, 015007 (2021), arXiv:2104.02982[hep-ph doi: 10.1103/PhysRevD.104.015007
    [58] J. Chen, Q. Wen, F. Xu et al., (2021), arXiv: 2104.03699[hep-ph]
    [59] W.-F. Duan, S.-P. Li, X.-Q. Li et al., (2021), arXiv: 2111.05178[hep-ph]
    [60] K. S. Babu, P. S. B. Dev, S. Jana et al., JHEP 03, 179 (2021), arXiv:2009.01771[hep-ph
    [61] A. K. Alok, A. Dighe, S. Gangal et al., JHEP 06, 089 (2019), arXiv:1903.09617[hep-ph
    [62] F. Munir Bhutta, Z.-R. Huang, C.-D. Lü et al., Nucl. Phys. B 979, 115763 (2022) doi: 10.1016/j.nuclphysb.2022.115763
    [63] A. K. Alok, S. Kumbhakar, J. Saini et al., Nucl. Phys. B 967, 115419 (2021), arXiv:2011.14668[hep-ph doi: 10.1016/j.nuclphysb.2021.115419
    [64] J. Kumar and D. London, Phys. Rev. D 99, 073008 (2019), arXiv:1901.04516[hep-ph
    [65] A. K. Alok, B. Bhattacharya, D. Kumar et al., Phys. Rev. D 96, 015034 (2017), arXiv:1703.09247[hep-ph doi: 10.1103/PhysRevD.96.015034
    [66] B. Capdevila, A. Crivellin, S. Descotes-Genon et al., JHEP 01, 093 (2018), arXiv:1704.05340[hep-ph
    [67] M. Algueró, P. A. Cartelle, A. M. Marshall et al., JHEP 12, 085 (2021) doi: 10.1007/JHEP12(2021)085
    [68] T. Hurth, F. Mahmoudi, D. Martinez Santos et al., Phys. Rev. D 96, 095034 (2017), arXiv:1705.06274[hep-ph
    [69] B. Capdevila, S. Descotes-Genon, J. Matias et al., JHEP 10, 075 (2016), arXiv:1605.03156[hep-ph
    [70] W. Altmannshofer, P. Stangl, and D. M. Straub, Phys. Rev. D 96, 055008 (2017), arXiv:1704.05435[hep-ph doi: 10.1103/PhysRevD.96.055008
    [71] S. Descotes-Genon, J. Matias, and J. Virto, Phys. Rev. D 88, 074002 (2013), arXiv:1307.5683[hep-ph doi: 10.1103/PhysRevD.88.074002
    [72] S. Descotes-Genon, L. Hofer, J. Matias et al., JHEP 06, 092 (2016), arXiv:1510.04239[hep-ph
    [73] W. Altmannshofer, P. Ball, A. Bharucha et al., JHEP 01, 019 (2009), arXiv:0811.1214[hep-ph
    [74] S. Jäger and J. Martin Camalich, JHEP 05, 043 (2013), arXiv:1212.2263[hep-ph
    [75] R. Mandal, R. Sinha, and D. Das, Phys. Rev. D 90, 096006 (2014), arXiv:1409.3088[hep-ph
    [76] A. K. Alok, S. Kumbhakar, and S. Uma Sankar, Eur. Phys. J. Plus 137(3), 354 (2022) doi: 10.1140/epjp/s13360-022-02576-z
    [77] T. Hurth, F. Mahmoudi, and S. Neshatpour, Nucl. Phys. B 909, 737 (2016), arXiv:1603.00865[hep-ph
    [78] S. Bhattacharya, A. Biswas, S. Nandi et al., Phys. Rev. D 101, 055025 (2020), arXiv:1908.04835[hep-ph
    [79] S. Jäger and J. Martin Camalich, Phys. Rev. D 93, 014028 (2016), arXiv:1412.3183[hep-ph
    [80] W. Altmannshofer, S. Gori, M. Pospelov et al., Phys. Rev. D 89, 095033 (2014), arXiv:1403.1269[hep-ph doi: 10.1103/PhysRevD.89.095033
    [81] M. Ciuchini, A. M. Coutinho, M. Fedele et al., Eur. Phys. J. C 77, 688 (2017), arXiv:1704.05447[hep-ph
    [82] M. Ciuchini, A. M. Coutinho, M. Fedele et al., Eur. Phys. J. C 79, 719 (2019), arXiv:1903.09632[hep-ph
    [83] M. Algueró, B. Capdevila, S. Descotes-Genon et al., Phys. Rev. D 99, 075017 (2019), arXiv:1809.08447[hep-ph
    [84] M. Algueró, B. Capdevila, S. Descotes-Genon et al., JHEP 07, 096 (2019), arXiv:1902.04900[hep-ph
    [85] R. Aaij et al. (LHCb), JHEP 05, 159 (2013), arXiv:1304.3035[hep-ex
    [86] R. Aaij et al. (LHCb), JHEP 04, 064 (2015), arXiv:1501.03038[hep-ex
    [87] R. Aaij et al. (LHCb), JHEP 12, 081 (2020), arXiv:2010.06011[hep-ex
  • 加载中

Figures(6) / Tables(4)

Get Citation
Shuang-Yi Li, Rui-Xiang Shi and Li-Sheng Geng. Discriminating 1D new physics solutions in bsℓℓ decays[J]. Chinese Physics C. doi: 10.1088/1674-1137/ac5a3a
Shuang-Yi Li, Rui-Xiang Shi and Li-Sheng Geng. Discriminating 1D new physics solutions in bsℓℓ decays[J]. Chinese Physics C.  doi: 10.1088/1674-1137/ac5a3a shu
Milestone
Received: 2022-01-22
Article Metric

Article Views(1209)
PDF Downloads(49)
Cited by(0)
Policy on re-use
To reuse of Open Access content published by CPC, for content published under the terms of the Creative Commons Attribution 3.0 license (“CC CY”), the users don’t need to request permission to copy, distribute and display the final published version of the article and to create derivative works, subject to appropriate attribution.
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Email This Article

Title:
Email:

Discriminating 1D new physics solutions in bsℓℓ decays

    Corresponding author: Li-Sheng Geng, lisheng.geng@buaa.edu.cn, Corresponding author
  • 1. School of Physics, Beihang University, Beijing 102206, China
  • 2. School of Space and Environment, Beihang University, Beijing 102206, China
  • 3. School of Physics & Beijing Key Laboratory of Advanced Nuclear Materials and Physics, Beihang University, Beijing 102206, China
  • 4. School of Physics and Microelectronics, Zhengzhou University, Zhengzhou 450001, China

Abstract: The recent measurements of $ R_{K^+} $, $ R_{K_S^0} $, $ R_{K^{*+}} $, $ B_s\to\mu^+\mu^- $, a set of CP-averaged angular observables for the $ B^0\to K^{*0}\mu^+\mu^- $ decay and its isospin partner $ B^+\to K^{*+}\mu^+\mu^- $ by the LHCb Collaboration consistently hint at lepton universality violation in the $ b\to s\ell\ell $ transitions. In this work, we first perform global fits to the $ b\to s\ell\ell $ data and show that five one-dimensional scenarios, i.e, $ \delta C_9^{\mu} $, $ \delta C_{10}^{\mu} $, $ \delta C_L^{\mu} $, $ \delta C_9^{\mu}=C_{10}^{\mu\prime} $, and $ \delta C_9^{\mu}=-C_9^{\mu\prime} $ can best explain the so-called B anamolies. Furthermore, we explore how these scenarios can be distinguished from each other. For this purpose, we first study the combinations of four angular asymmetries $ A_i ~~(i=3,4,5,9) $ and find that they cannot distinguish the five new physics scenarios. We then show that a newly constructed ratio $ R_{S} $ can uniquely discriminate the five new physics scenarios in proper intervals of $ q^2 $ if it can be measured with percent-level precision.

    HTML

    I.   INTRODUCTION
    • In the standard model (SM), the rare decays of B mesons induced by flavour-changing neutral-current (FCNC) $ b\to s\ell\ell $ transitions at the quark level are suppressed by the GIM mechanism [1]. Therefore they provide an ideal laboratory to indirectly probe new physics (NP) beyond the SM. Interestingly, the measurements of several observables yield results in tension with the SM expectations (see Refs. [28] for recent reviews).

      Recently, the LHCb Collaboration reported the most precise measurement of $ R_{K^+}=\Gamma(B^+\to K^+\mu^+\mu^-)/ \Gamma(B^+\to $$ K^+e^+e^-) $ in the bin $ [1.1,6]\; {\rm{GeV}}^2 $ [9]

      $ \begin{eqnarray} R_{K^+}=0.846^{+0.042 +0.013}_{-0.039 -0.012}, \end{eqnarray} $

      (1)

      showing that the significance of deviation from the SM prediction (see Table 1 in Ref. [10]) is at $ 3.1\sigma $ confidence level. Compared to the 2014 measurement [11], the tension with respect to the SM prediction has significantly increased. At the same time, a new result for the $ B_s\to\mu^+\mu^- $ branching fraction,

      Coeff. Best fit $ \chi^2_{\rm{min}} $ P-value $ {\rm{Pull}}_{\rm{SM}} $ 1σ range 3σ range
      $ \delta C_9^{\mu}=C_{10}^{\mu\prime} $ −0.49 118.75 [95 d.o.f.] 0.05 3.64 [−0.62, −0.35] [−0.90, −0.09]
      $ \delta C_9^{\mu}=-C_9^{\mu\prime} $ −1.03 117.99 [95 d.o.f.] 0.06 3.74 [−1.31, −0.76] [−1.91, −0.20]
      $ \delta C_9^{\mu} $ −0.90 108.05 [95 d.o.f.] 0.17 4.89 [−1.11, −0.70] [−1.55, −0.33]
      $ \delta C_{10}^{\mu} $ 0.57 109.64 [95 d.o.f.] 0.14 4.73 [0.45, 0.70] [0.20, 0.98]
      $ \delta C_L^{\mu} $ $ -0.42 $ 104.00 [95 d.o.f.] 0.25 5.29 $ [-0.50,-0.34] $ $ [-0.68,-0.18] $
      $ (\delta C_9^{\mu},\delta C_{10}^{\rm{\mu}}) $ $ (-0.58, 0.32) $ $ 103.51 $ [94 d.o.f.] $ 0.24 $ $ 4.97 $ $ \delta C_9^\mu \in ~[-0.81,-0.34] $ $ \delta C_{10}^\mu \in ~[0.17, 0.48] $

      Table 1.  Best fit values, $ \chi_{\rm{min}}^2 $, p-value, $ {\rm{Pull}}_{\rm{SM}} $ and confidence intervals of the WCs in the five one dimensional scenarios.

      $ \begin{eqnarray} {\rm{BR}}(B^0_s\to\mu^+\mu^-)=(3.09^{+0.46 +0.15}_{-0.43 -0.11})\times10^{-9}, \end{eqnarray} $

      (2)

      was also published by the LHCb Collaboration [12, 13], which allows for a better constraint on the Wilson coefficient $ C_{10} $. In Ref. [14], together with the theoretical prediction of Ref. [15] including the double-logarithmic QED and QCD corrections, the CMS measurement [16], and the ATLAS measurement [17], the following ratio is obtained

      $ \begin{eqnarray} R=\frac{{\rm{BR}}(B_s^0\to\mu^+\mu^-)_{\rm{exp}}}{{\rm{BR}}(B_s^0\to\mu^+\mu^-)_{\rm{SM}}}=0.78(9). \end{eqnarray} $

      (3)

      Very recently, the LHCb Collaboration measured two new lepton-universality ratios $ R_{K_S^0}=\Gamma(B^0\to K_S^0\mu^+\mu^-)/ $$ \Gamma(B^0\to K_S^0 e^+ e^-) $ and $ R_{K^{*+}}=\Gamma(B^+\to K^{*+}\mu^+\mu^-)/ \Gamma(B^+\to $$ K^{*+} e^+ e^-) $ in the $ q^2 $ ranges $ [1.1,6.0]\; \;{\rm{GeV}}^2 $ and $ [0.045,6.0] $$ {\rm{GeV}}^2 $, respectively, using proton–proton collision data corresponding to an integrated luminosity of 9 fb$ ^{-1} $ [18]. The two ratios are

      $ R_{K_S^0}=0.66_{\; -0.14\; -0.04}^{\; +0.20\; +0.02},\quad R_{K^{*+}}=0.70_{\; -0.13\; -0.04}^{\; +0.18\; +0.03}, $

      (4)

      where the first error is statistical and the second is systematic. We note that the central values are close to the Belle measurements [19, 20], but the uncertainties are smaller, resulting in deviations from the SM predictions with a significance of $ 1.5\sigma $ and $ 1.4\sigma $, respectively.

      Also, in 2020 the LHCb Collaboration reported angular analyses of the $ B^0\to K^{*0}\mu^+\mu^- $ decay and its isospin partner $ B^+\to K^{*+}\mu^+\mu^- $ decay [21, 22]. The results reconfirm the global tension with respect to the SM predictions previously reported for the decay of $ B^0\to K^{*0}\mu^+\mu^- $ [23].

      These new measurements have attracted much attention and led to many model-independent global analyses [14, 2436] assuming the presence of NP only in the $ b\to s\mu^+\mu^- $ mode, and it is shown that the significance of the SM exclusion in the global fits is about $ 4\sim6\sigma $. These data have also been explained in NP models involving tree-level exchanges of new particles such as the neutral gauge boson $ Z^\prime $ [3746], leptoquarks [31, 33, 34, 4757] or scalar Higgs [58, 59]. It is interesting to note that the leptoquark models can simultaneously explain the $ b\to s\ell^+\ell^- $ (with $ \ell=e,\mu $) and charged current (CC) $ b\to c\ell^-\bar{\nu}_\ell $ (with $ \ell=\mu,\tau $) flavor anomalies [31, 33, 34, 49, 53, 55, 60].

      The so-called B anamolies can be best explained in five one-dimensional scenarios, i.e, $ \delta C_9^{\mu} $, $ \delta C_{10}^{\mu} $, $ \delta C_L^{\mu} $, $ \delta C_9^{\mu}=C_{10}^{\mu\prime} $, and $ \delta C_9^{\mu}=-C_9^{\mu\prime} $, as demonstrated in the model independent analyses [14, 2428, 30, 31, 34, 6172]. In this work, assuming that NP appears in the muon mode, we first present global fits to the $ b\to s\ell^+\ell^- $ data for the SM and five one-dimensional NP scenarios $ \delta C_9^{\mu} $, $ \delta C_{10}^{\mu} $, $ \delta C_L^{\mu} $, $ \delta C_9^{\mu}=C_{10}^{\mu\prime} $, and $ \delta C_9^{\mu}=-C_9^{\mu\prime} $; the latter two were not considered in Ref. [14]. Next, we study the four angular asymmetries $ A_i $ ($ i=3,4,5,9 $) and find that they cannot discriminate the five new physics scenarios. Finally, we show that a set of ratios $ R_i $ are helpful to uniquely discriminate the NP scenarios. In addition, we present predictions in the SM and different NP scenarios for binned observables $ A_i $ and $ R_i $ in certain ranges of bins, which might be relevant for future experiments.

      This work is organized as follows. We introduce the observables $ A_i $ and $ R_i $ in Sec. III. Results and discussions are given in Sec. III, followed by a short summary and outlook in Sec. IV.

    II.   ANGULAR OBSERVABLES $ {\boldsymbol A_i} $ AND $ {\boldsymbol R_i} $ FOR THE ${\boldsymbol B\to {\boldsymbol K^*}{\boldsymbol \ell^+}}{\boldsymbol \ell^-}$ DECAY
    • For details about our theoretical framework, the low-energy effective Hamiltonian, the origin of theoretical uncertainties and how they are parameterized, and the statistical methods used to perform global fits to all the relevant experimental data, we refer to Ref. [14]. In this section, we introduce some angular observables which can potentially discriminate different NP scenarios. For this, we focus on the $ B\to K^*\ell^+\ell^- $ decay. The differential decay rate for the four-body $ B(p)\to K^*(k)(\to K\pi)\ell^+(q_1) $$ \ell^-(q_2) $ decay is of the following form [73, 74]

      $ \begin{eqnarray} &&\frac{{\rm d}^{(4)}\Gamma}{{\rm d}q^2 {\rm d}\cos\theta_{\ell}{\rm d}\cos\theta_K {\rm d}\phi} = \frac{9}{32\pi}I^{(\ell)}(q^2,\theta_\ell,\theta_K,\phi), \end{eqnarray} $

      (5)

      with

      $ \begin{aligned}[b] I^{(\ell)}(q^2,\theta_\ell,\theta_K,\phi)=&\left(I_1^s\sin^2\theta_K + I_1^c\cos^2\theta_K \right.\\ &+ (I_2^s\sin^2\theta_K + I_2^c\cos^2\theta_K)\cos2\theta_{\ell} \\ &+ I_3\sin^2\theta_K\sin^2 \theta_{\ell}\cos2\phi \\ &+ I_4\sin2\theta_K\sin2\theta_{\ell}\cos\phi\\ &+ I_5\sin2\theta_K\sin\theta_{\ell}\cos\phi \\ &+ (I_6^s\sin^2\theta_K + I_6^c\cos^2\theta_K)\cos\theta_{\ell} \\ &+ I_7\sin2\theta_K\sin\theta_{\ell}\sin\phi\\ &+ I_8\sin2\theta_K\sin2\theta_{\ell}\sin\phi\\ &\left.+ I_9\sin^2\theta_K\sin^2 \theta_{\ell}\sin2\phi\right), \end{aligned} $

      (6)

      where the kinematic variables $ q^2 $, $ \theta_{\ell} $, $ \theta_K $ and ϕ are defined, respectively, as follows: (i) $ q^2=(p-k)^2 $ is the square of the dilepton invariant mass, (ii) $ \theta_{\ell} $ is the angle between the flight direction of the B meson and the $ \ell^- $ lepton in the dilepton rest frame, (iii) $ \theta_K $ is the angle between momenta of the B meson and the K meson in the dimeson ($ K\pi $) rest frame, (iv) ϕ is the angle between the dimeson ($ K\pi $) and dilepton rest frames. As shown in Fig. 1, both $ \theta_{\ell} $ and $ \theta_K $ are defined in the interval $ [0,\pi] $ while the range of ϕ is $ [0,2\pi] $. The angular coefficients $ I_i $ expressed in terms of the helicity amplitudes can be found in Ref. [74], while the expressions of the CP-conjugate decay $ \bar{B}\to\bar{K}^*\ell^+\ell^- $ can be obtained by the following replacements

      Figure 1.  (color online) Kinematics of the $ B\to K^*\ell^+\ell^- $ decay.

      $ \begin{align} I_{1s(c),2s(c),3,4,7}\rightarrow \bar{I}_{1s(c),2s(c),3,4,7},\qquad I_{5,6,8,9}\rightarrow -\bar{I}_{5,6,8,9}, \end{align} $

      (7)

      where $ \bar{I}_i=I_i^{*} $.

      Following the strategy outlined in Refs. [75, 76], we can "construct" four angular asymmetries $ A_i $ ($ i=3,4,5,9 $) by integrating over different intervals in the angles $ \theta_\ell $, $ \theta_K $, and ϕ for the $ B\to K^*\mu^+\mu^- $ decay as follows:

      $ \begin{aligned}[b] A_3(q^2)=&\frac{\left[\int_{-\pi/4}^{\pi/4}-\int_{\pi/4}^{3\pi/4}+\int_{-\pi/4}^{3\pi/4}-\int_{3\pi/4}^{7\pi/4}\right]{\rm d}\phi \int_{-1}^{1}{\rm d}\cos\theta_K\int_{-1}^{1}{\rm d}\cos\theta_{\ell}\cdot I_{\rm{sum}}(q^2,\theta_\ell,\theta_K,\phi)} {\int_{0}^{2\pi}{\rm d}\phi\int_{-1}^{1}{\rm d}\cos\theta_K\int_{-1}^{1}{\rm d}\cos\theta_{\ell}\cdot I_{\rm{sum}}(q^2,\theta_\ell,\theta_K,\phi)},\\ A_4(q^2)=&\frac{\left[\int_0^{\frac{\pi}{2}}-\int_{\frac{\pi}{2}}^{\pi}\right]{\rm d}\phi \left[\int_{0}^{1}-\int_{-1}^{0}\right]{\rm d}\cos\theta_K\left[\int_{0}^{1}-\int_{-1}^{0}\right]{\rm d}\cos\theta_{\ell}\cdot I_{\rm{sum}}(q^2,\theta_\ell,\theta_K,\phi)} {\int_{0}^{2\pi}{\rm d}\phi\int_{-1}^{1}{\rm d}\cos\theta_K\int_{-1}^{1}{\rm d}\cos\theta_{\ell}\cdot I_{\rm{sum}}(q^2,\theta_\ell,\theta_K,\phi)},\\ A_5(q^2)=&\frac{\left[\int_0^{\frac{\pi}{2}}-\int_{\frac{\pi}{2}}^{\pi}\right]{\rm d}\phi \left[\int_{0}^{1}-\int_{-1}^{0}\right]{\rm d}\cos\theta_K\int_{-1}^{1}{\rm d}\cos\theta_{\ell}\cdot I_{\rm{diff}}(q^2,\theta_\ell,\theta_K,\phi)} {\int_{0}^{2\pi}{\rm d}\phi\int_{-1}^{1}{\rm d}\cos\theta_K\int_{-1}^{1}{\rm d}\cos\theta_{\ell}\cdot I_{\rm{sum}}(q^2,\theta_\ell,\theta_K,\phi)},\\ A_9(q^2)=&\frac{\left[\int_{0}^{\pi/2}-\int_{\pi/2}^{\pi}+\int_{0}^{\pi}-\int_{\pi}^{2\pi}\right]{\rm d}\phi \int_{-1}^{1}{\rm d}\cos\theta_K\int_{-1}^{1}{\rm d}\cos\theta_{\ell}\cdot I_{\rm{diff}}(q^2,\theta_\ell,\theta_K,\phi)} {\int_{0}^{2\pi}{\rm d}\phi\int_{-1}^{1}{\rm d}\cos\theta_K\int_{-1}^{1}{\rm d}\cos\theta_{\ell}\cdot I_{\rm{sum}}(q^2,\theta_\ell,\theta_K,\phi)}, \end{aligned} $

      (8)

      with

      $ \begin{aligned}[b] I_{\rm{sum}}(q^2,\theta_\ell,\theta_K,\phi)=&I^{(\ell)}(q^2,\theta_\ell,\theta_K,\phi)+\bar{I}^{(\ell)}(q^2,\theta_\ell,\theta_K,\phi),\\ I_{\rm{diff}}(q^2,\theta_\ell,\theta_K,\phi)=&I^{(\ell)}(q^2,\theta_\ell,\theta_K,\phi)-\bar{I}^{(\ell)}(q^2,\theta_\ell,\theta_K,\phi), \end{aligned} $

      (9)

      where the angular asymmetries $ A_i $ are related to the CP-averaged observables $ S_i $ introduced in Ref. [73], namely,

      $ \begin{eqnarray} A_3=\frac{1}{\pi}S_3,\quad A_4=\frac{1}{\pi}S_4,\quad A_5=\frac{3}{8}S_5,\quad A_9=\frac{1}{\pi}S_9. \end{eqnarray} $

      (10)

      Clearly, analogous to $ S_i $ [73, 77, 78], the observables $ A_i $ are also sensitive to NP. As a matter of fact, the relative uncertainties of observables $ A_i $ and those of $ S_i $ are equal. On the experimental side, these observables can be measured either by fitting to angular distributions, or by performing angular integration of $ I_i $. However, this may result in different statistical uncertainties because of different angular coverage. It should be noted that the observables $ A_i $ are ratios of combinations of the well-known angular coefficients $ I_i $. Therefore it is feasible to measure them in the current experiments.

      In addition to the $ A_i $ observables, following Refs. [10, 79], we revisit the ratios of angular observables $ R_i $,

      $ \begin{eqnarray} R_i(q^2)=\frac{I_i^{(\mu)}(q^2)+\bar{I}_i^{(\mu)}(q^2)}{I_i^{(e)}(q^2)+\bar{I}_i^{(e)}(q^2)}, \end{eqnarray} $

      (11)

      which are defined with the angular coefficients in Eq. (6). As pointed out in Ref. [79], the observables $ R_i $ are sensitive to new physics and their hadronic uncertainties are almost exactly cancelled. It should be emphasized that in this work the lepton mass $ m_\ell $ is non-vanishing and the angular coefficient $ I_6^c $ is absent because we do not consider scalar operators.

    III.   RESULTS AND DISCUSSIONS
    • In this section, we first update the global fits of Ref. [14] by considering the latest measurements $ R_{K_S^0} $ and $ R_{K^{*+}} $ [18] and perform global fits for two more NP scenarios, in comparison with Ref. [14]. Next, we study whether with the four asymmetry observables $ A_i $ and the ratios $ R_i $ introduced in Refs. [10, 79], one can distinguish different one-dimensional NP scenarios.

    • A.   Updated global fits

    • In the Appendix of Ref. [14], we have updated the clean fits insensitive to hadronic uncertainties by considering the latest measurements of $ R_{K_S^0} $ and $ R_{K^{*+}} $. It is shown that the three one-dimensional scenarios $ \delta C_9^{\mu} $, $ \delta C_{10}^{\mu} $ and $ \delta C_L^{\mu}=\delta C_9^{\mu}=-\delta C_{10}^{\mu} $ and a two-dimensional scenario $ (\delta C_9^\mu,\delta C_{10}^\mu) $ can well describe the current experimental data. Interestingly, a likelihood ratio test favours the $ \delta C_L^\mu $ scenario over the SM at $ 5\sigma $ significance. Additionally, although the angular observables of the $ B\to K^*\mu^+\mu^- $ decay suffer from large hadronic uncertainties, they can still provide constraints on new physics. Therefore, it is interesting to study how the conclusions of the clean fits change when including these angular observables. For this purpose, we update the global fits of Ref. [14]. We should note that the total number of fitted data for the updated global fits is 96. We obtain $ \chi_{\rm{SM,min}}^2=131.98 $ with 96 degrees of freedom (d.o.f.), corresponding to a p-value of $ 0.009 $. Compared to the global fits of Ref. [14], the updated fits given in Table 1 and Fig. 2 show that the confidence level of the exclusion of the SM predictions increases by $ 0.3\sigma - 0.38\sigma $ for all the NP scenarios. In particular, the significance of deviation from SM in the $ \delta C_L^\mu $ scenario is more than $ 5\sigma $ as well. In addition, we find that the updated global fit constrains better $ \delta C_9^\mu $ in the two-dimensional $ (\delta C_9^\mu,\delta C_{10}^\mu) $ scenario and excludes positive values at $ 1\sigma $ confidence level, as shown in Fig. 2.

      Figure 2.  (color online) Contour plots at the $ 1\sigma $ and $ 3\sigma $ confidence level for the $ (\delta C_9^{\mu},\delta C_{10}^{\mu}) $ scenario. The updated global fit with all observables for Gaussian (regions in light red and orange) is compared with the corresponding global fit of Ref. [14] (dashed lines in blue).

      We note that several other groups [24, 2629, 36] have also performed global fits considering the latest LHCb data. These fits show that three one-dimensional scenarios $ \delta C_9^{\mu} $, $ \delta C_{10}^{\mu} $ and $ \delta C_L^{\mu} $ can best explain the current experimental data, which is consistent with our conclusion. However, the significance of the SM exclusion is different among the various global fits, as shown in Table 2. The different significances mainly originate from two factors, i.e., which experimental data are fitted and how the uncertainties of non-perturbative effects and experimental data are dealt with. It is interesting to compare the results of Refs. [27, 36] and ours because only these three works considered the latest $ R_{K_S} $ and $ R_{K^{*+}} $ data. In Ref. [36] a "data driven" approach is used to conservatively estimate the theoretical uncertainties, resulting in the smallest $ {\rm{Pull}}_{\rm{SM}} $ values for the scenarios of $ \delta C^\mu_9 $ and $ \delta C^\mu_L $, but somehow a larger $ {\rm{Pull}}_{\rm{SM}} $ than ours for the scenario of $ \delta C^\mu_{10} $. On the other hand, the $ {\rm{Pull}}_{\rm{SM}} $ values of Ref. [27] are larger than ours.

      Coeff. [24] [26] [27] [28] [29] [36] This work
      $ \delta C_9^{\mu} $ $ 5.2\sigma $ $ 6.3\sigma $ $ 7.0\sigma $ $ 7.6\sigma $ $ 6.7\sigma $ $ 2.1\sigma $ $ 4.89\sigma $
      $ \delta C_{10}^{\mu} $ $ 4.7\sigma $ $ 4.5\sigma $ $ - $ $ 4.5\sigma $ $ 6.0\sigma $ $ 5.1\sigma $ $ 4.73\sigma $
      $ \delta C_L^{\mu} $ $ 5.6\sigma $ $ 5.9\sigma $ $ 6.5\sigma $ $ 6.7\sigma $ $ 7.0\sigma $ $ 4.3\sigma $ $ 5.29\sigma $

      Table 2.  Significance of the SM exclusion in the three NP scenarios from various global fits. Note that only Refs. [27, 36] and our present work considered the latest $ R_{K_S} $ and $ R_{K^{*+}} $ data [18].

      In fact, in the case of NP involving tree-level exchanges of new particles such as the neutral gauge boson $ Z^\prime $, leptoquarks or scalar Higgs, further correlations can be induced between the Wilson coefficients $ C_9^{\mu(\prime)} $ and $ C_{10}^{\mu(\prime)} $. As demonstrated in some model independent analyses [27, 61, 6466, 69, 72], two additional one-dimensional scenarios $ \delta C_9^{\mu}= C_{10}^{\mu\prime} $ and $ \delta C_9^{\mu}=- C_9^{\mu\prime} $ cannot be completely excluded by the current $ b\to s\ell\ell $ data. The former scenario can be realized in the $ U_1 $ leptoquark model in the case of $ \texttt{g}_{\ell q}^{\mu b}\left(\texttt{g}_{\ell q}^{\mu s}\right)^*=\texttt{g}_{ed}^{\mu b}\left(\texttt{g}_{ed}^{\mu s}\right)^* $ [49]. The latter one could be realized in $ Z^\prime $ models with vector-like fermions and $ L_\mu-L_\tau $ symmetry [80]. Therefore, we would like to investigate these two NP scenarios in our theoretical framework. This can be done following the same strategy as that of Ref. [14] and the corresponding results are shown in Table 1. We note that the $ \delta C_9^{\mu}=C_{10}^{\mu\prime} $ and $ \delta C_9^{\mu}=-C_9^{\mu\prime} $ scenarios can well explain the new $ b\to s\ell\ell $ data and the deviation from the SM has a significance of more than $ 3\sigma $, smaller than the three one-dimensional $ \delta C_9^{\mu} $, $ \delta C_{10}^{\mu} $ and $ \delta C_L^{\mu} $ scenarios. Indeed, we also tested the other one-dimensional scenarios induced by NP models with the new data studied in the present work and found that they are excluded at $ 3\sigma $ confidence level, consistent with the conclusions of Refs. [27, 61, 6466, 69, 72].

      As discussed above, among the new NP scenarios studied, all of the five one-dimensional cases, i.e., $ \delta C_9^{\mu} $, $ \delta C_{10}^{\mu} $, $ \delta C_L^{\mu} $, $ \delta C_9^{\mu}= C_{10}^{\mu\prime} $ and $ \delta C_9^{\mu}=-C_9^{\mu\prime} $, can provide a good description of the latest data. Thus, a natural question is which observables can discriminate between them. For such a purpose, we study the four angular asymmetries $ A_i $ ($ i=3,4,5,9 $) and the ratios $ R_i $ introduced in Sec. II.

    • B.   Sensitivity of angular asymmetries ${\boldsymbol A_i}$ to new physics

    • In Fig. 3, we plot $ A_{3,5,4,9} $ as functions of $ q^2 $, where q is the dilepton invariant mass, with the best fitted Wilson coefficients determined for the SM and the five one-dimensional NP scenarios. The shaded bands reflect the uncertainties originating from form factors [10, 74, 79]. The asymmetries integrated over [15, 19] GeV$ ^2 $ are given in Table 3.

      Figure 3.  (color online) Dilepton mass squared ($ q^2 $) spectra of observables $ A_i $ in the SM (solid lines in red) and different NP scenarios (dashed lines). Shaded regions denote the uncertainties of SM and five NP scenarios.

      Observables SM $ \delta C_9^{\mu} $ $ \delta C_{10}^{\mu} $ $ \delta C_9^{\mu}= C_{10}^{\mu\prime} $ $ \delta C_L^{\mu} $ $ \delta C_9^{\mu}=- C_9^{\mu\prime} $
      $ A_3 $ [15, 19] GeV$ ^2 $ $ -0.067^{+0.006}_{-0.006} $ $ -0.067^{+0.006}_{-0.006} $ $ -0.067^{+0.006}_{-0.006} $ $ -0.060^{+0.007}_{-0.007} $ $ -0.067^{+0.006}_{-0.006} $ $ -0.054^{+0.007}_{-0.007} $
      $ A_4 $ [15, 19] GeV$ ^2 $ $ -0.096^{+0.002}_{-0.002} $ $ -0.096^{+0.002}_{-0.002} $ $ -0.096^{+0.002}_{-0.002} $ $ -0.093^{+0.003}_{-0.002} $ $ -0.096^{+0.002}_{-0.002} $ $ -0.090^{+0.003}_{-0.003} $
      $ A_5 $ [15, 19] GeV$ ^2 $ $ -0.104^{+0.008}_{-0.008} $ $ -0.101^{+0.008}_{-0.008} $ $ -0.104^{+0.008}_{-0.008} $ $ -0.113^{+0.008}_{-0.008} $ $ -0.104^{+0.008}_{-0.008} $ $ -0.115^{+0.008}_{-0.008} $
      $ A_9 $ [15, 19] GeV$ ^2 $ $ -0.0003^{+0.0001}_{-0.0002} $ $ -0.0003^{+0.0001}_{-0.0002} $ $ -0.0003^{+0.0002}_{-0.0002} $ $ -0.0003^{+0.0002}_{-0.0002} $ $ -0.0003^{+0.0002}_{-0.0002} $ $ -0.0013^{+0.0004}_{-0.0002} $

      Table 3.  Predictions in the SM and five NP scenarios for binned observables.

      By examining Fig. 3 and Table 3, we can see that these observables can provide two possible strategies to distinguish between different NP structures. The first strategy is to analyse the $ q^2 $ spectra of certain observables. It is clear that the discriminating power of these observables $ A_i(q^2) $ is very limited. For instance, to distinguish the NP scenario $ \delta C_9^{\mu}=- C_9^{\mu\prime} $, an experimental uncertainty of sub-percent level or lower is needed. The second strategy is to study binned observables in proper ranges of $ q^2 $. As shown in Table 3, only binned $ A_3 $ and $ A_9 $ have the ability to discriminate the scenarios $ \delta C_9^{\mu}=-C_9^{\mu\prime} $ and $ \delta C_9^{\mu}=C_{10}^{\mu\prime} $ at $ 3\sigma $ confidence level if they can be measured with an uncertainty of 0.001. However, the other NP scenarios cannot be discriminated by these observables $ A_i $.

      In conclusion, the angular asymmetries $ A_i $ cannot uniquely distinguish between the five different NP solutions in the one-dimensional scenario.

    • C.   Sensitivity of ratios $ {\boldsymbol R_i} $ to new physics

    • In this subsection, we study whether the ratios $ R_i $ introduced in Refs. [10, 79] can be used to discriminate the five NP scenarios.

      First we plot the $ q^2 $ spectra of $ R_i(q^2) $ in the low bin of [0.045,6] GeV$ ^2 $ in Fig. 4. One can see that the two observables $ R_{1s}(q^2) $ and $ R_{2s}(q^2) $ could discriminate all the five different NP scenarios, while the observables $ R_{1c}(q^2) $ or $ R_{2c}(q^2) $ can only discriminate the NP scenario $ \delta C_9^{\mu}=- C_9^{\mu\prime} $. As explicitly shown in Appendix B, the other ratios suffer from uncertainties induced by power corrections and charm loop contributions and have no ability to distinguish the NP scenarios. We stress that it is possible to distinguish the five different NP scenarios if future experimental statistics are high enough that an uncertainty of percent level or lower can be achieved.

      Figure 4.  (color online) Same as Fig. 3 but for observables $ R_{1c,1s,2c,2s,C,S} .$

      Next, we study binned observables in proper ranges of $ q^2 $. As shown in Table 4, once $ R_{1s} $ or $ R_{2s} $ is measured with the same accuracy as that of theoretical predictions, the five different NP scenarios can be either confirmed or excluded at $ 2\sigma $ confidence level. Therefore, the binned observables $ R_{1s} $ and $ R_{2s} $ are useful to distinguish between different NP scenarios, though less constraining than the $ q^2 $ spectra. We therefore conclude that the ratios $ R_{1s} $ and $ R_{2s} $ could uniquely distinguish between the five different NP solutions in the one-dimensional scenario if they can be measured with a precision of percent level.

      Observables SM $ \delta C_9^{\mu} $ $ \delta C_{10}^{\mu} $ $ \delta C_9^{\mu}= C_{10}^{\mu\prime} $ $ \delta C_L^{\mu} $ $ \delta C_9^{\mu}=- C_9^{\mu\prime} $
      $ R_{1s} $ [1, 4] GeV$ ^2 $ $ 0.982^{+0.002}_{-0.002} $ $ 1.14^{+0.02}_{-0.03} $ $ 0.84^{+0.03}_{-0.02} $ $ 1.06^{+0.03}_{-0.03} $ $ 0.94^{+0.03}_{-0.03} $ $ 1.19^{+0.02}_{-0.03} $
      $ R_{2s} $ [1, 4] GeV$ ^2 $ $ 0.966^{+0.001}_{-0.001} $ $ 1.12^{+0.01}_{-0.02} $ $ 0.83^{+0.02}_{-0.02} $ $ 1.04^{+0.02}_{-0.03} $ $ 0.92^{+0.02}_{-0.03} $ $ 1.17^{+0.02}_{-0.02} $
      $ R_{1c} $ [1, 6] GeV$ ^2 $ $ 1.007^{+0.001}_{-0.001} $ $ 0.83^{+0.02}_{-0.01} $ $ 0.85^{+0.02}_{-0.02} $ $ 0.77^{+0.02}_{-0.01} $ $ 0.80^{+0.01}_{-0.01} $ $ 0.68^{+0.06}_{-0.03} $
      $ R_{2c} $ [1, 6] GeV$ ^2 $ $ 0.977^{+0.001}_{-0.001} $ $ 0.81^{+0.02}_{-0.01} $ $ 0.83^{+0.02}_{-0.02} $ $ 0.75^{+0.02}_{-0.01} $ $ 0.78^{+0.01}_{-0.01} $ $ 0.66^{+0.05}_{-0.03} $
      $ R_{S} $ [1, 4] GeV$ ^2 $ $ 0.985^{+0.003}_{-0.002} $ $ 1.14^{+0.02}_{-0.03} $ $ 0.84^{+0.03}_{-0.02} $ $ 1.07^{+0.03}_{-0.03} $ $ 0.95^{+0.03}_{-0.03} $ $ 1.19^{+0.02}_{-0.03} $
      $ R_{C} $ [1, 6] GeV$ ^2 $ $ 0.9994^{+0.0007}_{-0.0008} $ $ 0.82^{+0.02}_{-0.01} $ $ 0.85^{+0.02}_{-0.02} $ $ 0.77^{+0.02}_{-0.01} $ $ 0.80^{+0.01}_{-0.01} $ $ 0.68^{+0.06}_{-0.03} $

      Table 4.  Same as Table 3 but for observables $ R_{1c,1s,2c,2s,C,S} $.

      On the experimental side, the angular coefficients $ I_{1s} $ and $ I_{2s} $ are not directly measured. Instead, the measured quantities are the branching ratio $ {\rm{BR}} $ and the longitudinal polarization $ F_L $ as combinations of $ I_{1c} $, $ I_{1s} $, $ I_{2c} $, and $ I_{2s} $, which read

      $ \begin{aligned}[b] \frac{{\rm d}{\rm{BR}}}{{\rm d}q^2}=&\frac{\tau_B}{2}\frac{{\rm d}\left(\Gamma+\bar{\Gamma}\right)}{{\rm d}q^2}=\frac{\tau_B}{4}\left[\left(3\Sigma_{1c}(q^2)-\Sigma_{2c}(q^2)\right)\right. \\&\left.+2\left(3\Sigma_{1s}(q^2)-\Sigma_{2s}(q^2)\right)\right], \end{aligned} $

      (12)

      $ \begin{eqnarray} F_L=\frac{1}{2}\frac{3\Sigma_{1c}(q^2)-\Sigma_{2c}(q^2)}{{\rm d}\left(\Gamma+\bar{\Gamma}\right)/{\rm d}q^2}, \end{eqnarray} $

      (13)

      where

      $ \begin{eqnarray} \Sigma_{1s,1c,2s,2c}(q^2)=I_{1s,1c,2s,2c}(q^2)+\bar{I}_{1s,1c,2s,2c}(q^2). \end{eqnarray} $

      (14)

      Clearly, the four angular coefficients $ I_{1c} $, $ I_{1s} $, $ I_{2c} $, and $ I_{2s} $ cannot be uniquely determined by two observables $ {\rm{BR}} $ and $ F_L $. However, we note that two combinations of angular coefficients $ 3\Sigma_{1c}-\Sigma_{2c} $ and $ 3\Sigma_{1s}-\Sigma_{2s} $ can be extracted from $ {\rm{BR}} $ and $ F_L $. Therefore, we define two new ratios

      $ \begin{eqnarray} &&R_C(q^2)=\frac{3\Sigma_{1c}^\mu(q^2)-\Sigma_{2c}^\mu(q^2)}{3\Sigma_{1c}^e(q^2)-\Sigma_{2c}^e(q^2)}, \end{eqnarray} $

      (15)

      $ \begin{eqnarray} &&R_S(q^2)=\frac{3\Sigma_{1s}^\mu(q^2)-\Sigma_{2s}^\mu(q^2)}{3\Sigma_{1s}^e(q^2)-\Sigma_{2s}^e(q^2)}. \end{eqnarray} $

      (16)

      It should be feasible to extract them experimentally. Therefore, we study the $ q^2 $ spectra and binned values of $ R_{C,S}(q^2) $ in the SM and different NP scenarios. The relevant results are shown in Table 4 and Fig. 4. It is clear that the ratio $ R_S $ could uniquely distinguish between the five different NP solutions in the one-dimensional scenario if it can be measured with a precision of percent level.

    IV.   SUMMARY AND OUTLOOK
    • In conclusion, in the one-dimensional NP case only five scenarios, i.e., $ \delta C_9^{\mu} $, $ \delta C_{10}^{\mu} $, $ \delta C_L^{\mu} $, $ \delta C_9^{\mu}= C_{10}^{\mu\prime} $ and $ \delta C_9^{\mu}=-C_9^{\mu\prime} $, can well describe the latest LHCb, Belle, ATLAS, and CMS data. In order to discriminate the five different NP scenarios, we first "construct" four angular asymmetries $ A_i $ and study the sensitivity of $ A_i $ to new physics. We find that these observables $ A_i $ cannot discriminate different NP scenarios. Further, we construct a new ratio $ R_S $ by revisiting a set of ratios $ R_i $ and show that the observable $ R_S $ can uniquely discriminate the five NP solutions in the one-dimensional case once they are measured in future experiments, by analysing the $ q^2 $ spectra of these observables and their binned results.

      In the next few years, with the collection of more data at the LHCb, Belle II experiments and improvement of experimental precision, such global fits should be continually updated. In addition, new ideas should be explored to further reduce hadronic uncertainties. Furthermore, it will be interesting to study the interplay between semi-leptonic baryon decays, such as $ \Lambda_b\to\Lambda $, $ \Xi_b\to\Xi $, and semi-leptonic meson decays. Because of rich helicity structures, if measured to the same precision, semi-leptonic baryon decays could constrain some relevant Wilson coefficients more tightly than their mesonic counterparts.

    ACKNOWLEDGMENTS
    • We thank Jorge Martin Camalich for providing the code for the calculation of the $ b\to s \ell\ell $ amplitudes.

    APPENDIX

      A.   Discussion of NP in the electron and muon channels

    • In our current work and most works performed so far, NP is assumed to only originate from the muon mode. On the other hand, one can also assume that NP exists in the electron mode, or a combination of the electron and muon modes. In the electron mode, the signs of the fitted Wilson coefficients should be opposite to those of the muon mode if only the observables $ R_K $ and $ R_{K^*} $ are taken into account. The case of allowing for NP in the electron channel has been discussed quite extensively in Refs. [68, 81] (for some recent studies see, e.g., Refs. [61, 64, 8284]). Here, we investigate a more general scenario that NP comes from a combination of electron and muon channels. As a result, more experimental data need to be taken into account, such as the branching fraction [85] as well as eight angular observables [86, 87] for the $ B \to K^* e^+ e^- $ decay. Thus, the total number of data fitted becomes 103.

      In Fig. A1, we show the two-parameter fits for the ($ \delta C_L^\mu $, $ \delta C_L^e $) and ($ \delta C_L^\mu $, $ \delta C_9^{\rm{univ}} $) cases including the full 103 data set. Here, one requires $ \delta C_9^{\rm{univ}}=\delta C_9^\mu =\delta C_9^e $. Compared to the results in 2017, all the NP scenarios are better constrained as shown by the smaller contours. In addition, the pure $ \delta C_L^{\rm{univ}} $ scenario does not change significantly and is still in agreement with the SM within the $ 1\sigma $ confidence level while other scenarios are not. More specifically, the significance of the SM exclusion in the pure $ \delta C_L^\mu $ scenario becomes larger. Also, the new data favor a zero $ \delta C_9^e $ contribution. It should be mentioned that these variations from 2017 to 2021 are mainly caused by the most precise $ R_K $ from the LHCb experiment. The results indicate that new physics is most likely to exist in the muon mode.

      Figure A1.  (color online) Contours for the two-parameter fits ($ \delta C_L^\mu $, $ \delta C_L^e $) and ($ \delta C_L^\mu $, $ \delta C_9^{\rm{univ}} $) including the full data set (103 data). The regions in orange and light red represent the 1σ and 3σ bounds. The lines in red correspond to those in Fig. 7 of Ref. [10].

    • B.   $ {\boldsymbol q}^{\boldsymbol 2} $ spectra of ${\boldsymbol R_{\boldsymbol {3-9}}}$

    • In Fig. B1, we show the $ q^2 $ spectra of $ R_{3-9} $. Clearly, these ratios cannot discriminate the five NP scenarios.

      Figure B1.  (color online) Same as Fig. 4 but for observables $ R_{3-9} $.

Reference (87)

目录

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return