-
The test of CP violation in the Standard Model, which is described by the complex phase in the Cabibbo-Kobayashi-Maskawa (CKM) matrix, is the central goal of the heavy flavor physics programme. Using B decays to determine the three angles
α ,β andγ of the usual non-squashed unitarity triangle of the CKM matrix, and thus to test the closure of the unitarity triangle, is a straightforward and promising way to accomplish this goal. Any discrepancies would suggest possible new sources of CP violation beyond the Standard Model.In principle,
α ,β andγ can be determined via the measurements of CP violating asymmetry in neutral B decays to CP eigenstates. If a single CKM amplitude or different amplitudes with the same CKM phase contribute to the decay ofB0 meson, then the mixing-induced CP asymmetry is a pure function of CKM parameters, which are derived from the two neutral B meson mixing andB0 decay with no strong phase uncertainties. As is well known, the angleβ can be determined in a reliable way with the help of the mixing-induced CP violation in the single "gold plated" modeB0→J/ψKS . Likewise,α can be extracted using the neutral B decayB0→π+π− , using the isospin symmetry to separate the tree and penguin contributions in theB→ππ decays. Similarly to the measurement ofβ andα , a straightforward way to obtainγ theoretically might be to use the CKM suppressedB0s decayB0s→ρKS , or to analyze the decaysB0s→D0ϕ,ˉD0ϕ andD01ϕ , as in [1]. However, the observed mixing-induced CP asymmetries are expected to be strongly diluted by the largeBs−¯Bs mixing, so that determiningγ in this way is considerably more involved than it is forβ andα .The third angle
γ is currently one of the least known. It depends on the difference between the strong phases of the amplitude contributions in B decays, which is difficult to calculate reliably. One of the theoretically best ways of determiningγ is to use the interference between theb→cˉus andb→uˉcs decay amplitudes, with the intermediate statesD0 and¯D0 mesons subsequently decaying into common final states. Depending on the common final states, the methods can be divided into: the GLW method [2] , where D mesons decay into CP eigenstates; the ADS method [3], where the final state is not a CP eigenstate but uses the doubly Cabibbo-suppressed decays to enhance the CP violation effect; the GGSZmethod [4], which exploits the three-body D decays to self-conjugate modes, such asD0→KSK+(π+)K−(π−) . The world average values are obtained by combining the above methods in order to reduce the statistical uncertainties, and areγ=(71.1+4.6−5.3)∘ [5],γ=(73.5+4.2−5.1)∘ [6] andγ=(70.0±4.2)∘ [7] fitted by the HFLAV, CKMFitter and UTfit collaborations, respectively. The latest measurement ofγ by the LHCb collaboration yielded(74.0+5.0−5.8)∘ [8]. As of mid 2019 , the preliminary result by CKMfitter forγ from all sources isγ=(65.80+0.94−1.29)∘ [9]. As there is no penguin diagram pollution in the charmed B decays, and as the theoretical uncertainty onγ determined from higher-order electroweak corrections is very small, a shiftδγ≲O(10−7) calculated in [10], the uncertainty of approximately5∘ onγ is statistically limited. The reason for the relatively large statistical uncertainty is due to the small branching fraction Br(B−→¯D0K− )∼O(10−6) .The two-body charmless B meson decays receive contributions from the tree and penguin diagrams with relatively large branching ratios, of the order of
O(10−6− 10−5) . The branching ratios and CP asymmetry parameters depend strongly on the interference of the tree and penguin diagrams with different weak and strong phases. This provides a possible way to measure the CKM angleγ . The only problem here is how to calculate or extract reliably the different strong phases of the tree and penguin diagrams of charmless B decays. The methods proposed in Refs. [11–13] extract the strong phases in theB→ππ ,B→πK andB→KK decays by using the flavor SU(3) symmetry. Fleischer proposed a method to use the decaysBd→π+π− andBs→K+K− via the U spin flavor symmetry of strong interactions [14]. All these methods require a number of experimental measurements but do not depend on the non-perturbative QCD calculations. However, the precision of determiningγ is limited by the theoretical uncertainties of the flavor SU(3) breaking effects or U spin breaking corrections. Some of them can only provide a bound onγ , and serve as a complementary and indirect constraint on the unitarity triangle. Recently, three-body charmless B decays, whose amplitudes are related by the flavor SU(3) symmetry [15–17] or U spin flavor symmetry [18], were used to extract the CKM angleγ . The uncertainty of this fit was of the order of10∘ with six possible solutions found in the latest paper [17].In order to improve the precision of
γ angle measurements, one has to deal with the flavor SU(3) breaking effect. Recently, the factorization assisted topological amplitude approach was proposed in [19–23] to parameterize all contributions to charmless B decays by topological diagrams, but keeping most of the SU(3) breaking effects. Similarly to the previous version of the topological diagram approach [24], most of the hadronic decay amplitudes of the weak diagrams are fitted using the experimental measurements instead of the perturbative QCD calculations. Thus, the fit is model independent. We also take into account the flavor SU(3) breaking effects in each flavor topological diagram characterized by decay constants and weak transition form factors. As a result, we can reduce the number of unknown hadronic parameters by fitting all charmless two-bodyB→PP ,PV decays together [19], while the study in [24] fitted theB→PP andB→PV decay modes separately with two sets of parameters. In the present work, we use all experimentally measuredB→PP ,PV decay observables to perform a global fit, but leave the weak phaseγ to be fitted from the abundant experimental data together with the hadronic parameters. There are also a number of measured B decay channels with two vector meson final states. Since this kind of decay is more complicated as it includes transverse polarization, we do not include them in our current study to avoid introducing more free parameters. Therefore, in the present work, we fit 15 parameters using 37 experimentally measured branching fractions and 11 CP asymmetry parameters of theB→PP ,PV decays.We begin in Sec. 2 with a summary of the parametrization of the tree and penguin topological amplitudes of the charmless
B→PP PV decays, which leave the weak phaseγ as a free parameter, to be fitted together with the hadronic parameters. The fit result of the CKM angleγ with the experimental and theoretical uncertainties is presented in Sec. 3, and the conclusions are given in Sec. 4. -
The charmless two-body B meson decays are induced by weak interactions via the tree diagram and the penguin diagram at the quark level. The tree level weak decay diagram contributes via the so-called color-favored tree emission diagram T, color-suppressed tree emission diagram C, W-exchange tree diagrams E, and the W annihilation tree diagrams A, which are shown in Fig. 1. The 1-loop corrections of the QCD penguin diagrams are not suppressed due to the larger CKM matrix element compared with the tree diagram. They are also grouped into four categories: (a) the QCD-penguin emission diagram P, (b) the flavor-singlet QCD-penguin diagram
PC orEW-penguin diagram PEW , (c) the time-like QCD-penguin diagramPE , and (d) the space-like QCD-penguin annihilation diagramPA , shown in Fig. 2.Figure 1. Topological tree diagrams contributing to
B→PP andB→PV decays: (a) the color-favored tree emission diagram T, (b) the color-suppressed tree emission diagram C, (c) the W-exchange diagram E , and (d) the W-annihilation diagram A.Figure 2. Topological penguin diagrams contributing to
B→PP andB→PV decays: (a) the QCD-penguin emission diagram P, (b) the flavor-singlet QCD-penguin diagramPC orEW-penguin diagram PEW , (c) the time-like QCD-penguin diagramPE , and (d) the space-like QCD-penguin annihilation diagramPA .The color-favored tree (T) topology shown in Fig. 1(a)was proven to be factorizable to all orders of
αs in the QCD factorization approach [25], the perturbative QCD factorization approach [26] and the soft-collinear effective theory [27], and it is easily parametrized asTP1P2=iGF√2VubV∗uq′a1(μ)fp2(m2B−m2p1)FBP10(m2p2),TPV=√2GFVubV∗uq′a1(μ)fVmVFB−P1(m2V)(ε∗V⋅pB),TVP=√2GFVubV∗uq′a1(μ)fPmVAB−V0(m2P)(ε∗V⋅pB).
(1) The superscript in
TP1P2 denotes that the final mesons are two pseudoscalar mesons, and inTPV(VP) that the recoiling mesons are pseudoscalar meson (vector meson).q′=d,s quark in theΔS=0,1 transition, respectively.a1(μ) is the effective Wilson coefficient from the short distance QCD corrections,a1(μ)=C2(μ)+C1(μ)/3=1.05 at next-to-leading order [28], with factorization scaleμ , insensitive to final state mesons, usually chosen atmb/2=2.1GeV . The decay constantsfP ,fV and the form factorsFBP10 ,FB−P1 andAB−V0 characterize theSU(3) breaking effects.ε∗V is the polarization vector of vector meson, andpB is the 4-momentum of B meson.We parametrize the non-factorization dominant color suppressed tree diagrams as
CP1P2=iGF√2VubV∗uq′χCeiϕCfp2(m2B−m2p1)FBP10(m2p2),CPV=√2GFVubV∗uq′χC′eiϕC′fVmVFB−P1(m2V)(ε∗V⋅pB),CVP=√2GFVubV∗uq′χCeiϕCfPmVAB−V0(m2P)(ε∗V⋅pB),
(2) where
χC ,ϕC represent the magnitude and phase of the pseudo-scalar meson emitted decaysB→PP ,VP . The prime inχC′ ,ϕC′ denotes the difference with respect to the hadronic parameters in the amplitude of the vector meson emitted decayB→PV . Similarly, we parametrize the W exchange diagrams asEP1P2=iGF√2VubV∗uq′χEeiϕEfBm2B(fp1fp2f2π),EPV,VP=√2GFVubV∗uq′χEeiϕEfBmV(fPfVf2π)(ε∗V⋅pB),
(3) where
χE ,ϕE represent the magnitude and strong phase. We ignore the W annihilation topology, as its contribution is negligible, as discussed in [24].The penguin emission diagram was also proven to be factorizable to all orders of
αS in various QCD inspired approaches and the soft-collinear effective theory. Thus, its amplitudes are describe as follows:PPP=−iGF√2VtbV∗tq′[a4(μ)+χPeiϕPrχ]fp2(m2B−m2p1)×FBP10(m2p2),PPV=−√2GFVtbV∗tq′a4(μ)fVmVFB−P1m2V(ε∗V⋅pB),
PVP=−√2GFVtbV∗tq′[a4(μ)−χPeiϕPrχ]×fPmVAB−V0(m2P)(ε∗V⋅pB),
(4) where the second terms denote the “chiral enhanced” penguin contributions, with the chiral factor of pseudo-scalar meson
rχ . The Wilson coefficienta4(μ) of effective penguin operators is calculated to the next-to-leading order [28]. We parametrize the flavor-singlet penguin diagram asPPPC=−iGF√2VtbV∗tq′χPCeiϕPCfp2(m2B−m2p1)FBP10(m2p2),PPVC=−√2GFVtbV∗tq′χP′CeiϕP′CfVmVFB−P1(m2V)(ε∗V⋅pB),PVPC=−√2GFVtbV∗tq′χPCeiϕPCfPmVAB−V0(m2P)(ε∗V⋅pB).
(5) Looking carefully at Fig. 2(a) and 2(d), one can see that the only difference between them is the hard gluon attached to a light quark pair. Since we do not calculate the QCD corrections but fit them to the experimental measurements, these two diagrams give the same contribution. Since the contribution of pseudo-scalar meson emission
PPP,VPA is already encoded in the parametersrχχP ,eiϕP in Eq. (4) for diagram Fig. 2(a), we only have the vector meson contribution left for the space-like penguin diagrams,PPVA=−√2GFVtbV∗tq′χPAeiϕPAfBmV(fPfVf2π)(ε∗V⋅pB).
(6) The contribution from the time-like penguin (
PE ) diagram is arguably small, and can be ignored in the decay modes where it is not dominant [19].The electroweak penguin topology (
PEW ) is similar to the tree and penguin emission diagrams, which can be proven by factorization. We calculate this diagram in QCD asPPPEW=−iGF√2VtbV∗tq′eq32a9(μ)fp2(m2B−m2p1)FBP10(m2p2),PPVEW=−√2GFVtbV∗tq′eq32a9(μ)fVmVFB−P1(m2V)(ε∗V⋅pB),PVPEW=−√2GFVtbV∗tq′eq32a9(μ)fPmVAB−V0(m2P)(ε∗V⋅pB),
(7) where
a9(μ) is the effective Wilson coefficient at the next-to-leading order. -
From Eqs. (1-7), one notices that all tree amplitudes are proportional to
VubV∗uq′ while the penguin amplitudes are proportional toVtbV∗tq′=−(VubV∗uq′+VcbV∗cq′) . ExceptVub≡|Vub|e−iγ , all other CKM matrix elementsV∗uq′ ,VcbV∗cq′ are approximately real numbers without the electroweak phase. The phase angleγ is a free parameter to be fitted from the experimental data. The magnitudes of CKM matrix elements are input parameters from Ref. [29]|Vud|=0.97420±0.00021,|Vus|=0.2243±0.0005,|Vub|=0.00394±0.00036,|Vcd|=0.218±0.004,|Vcs|=0.997±0.017,|Vcb|=0.0422±0.0008.
(8) The remaining parameters in the decay amplitude formulas are the meson decay constants and transition form factors. The meson decay constants are measured by experiments or calculated by theoretical approaches, such as the covariant light front approach [30], light-front holographic QCD [31], light-cone sum rules [32, 33], and QCD sum rules [34, 35]. We show their values in Table 1as an average from PDG [29].
fπ fK fB fρ fK∗ fω fϕ 130.2±1.7 155.6±0.4 190.9±4.1 213±11 220±11 192±10 225±11 Table 1. The decay constants of light pseudo-scalar mesons and vector mesons (in unit of MeV).
The transition form factors of B meson decays are usually measured via the semileptonic B decay modes together with the CKM matrix elements. Theoretically, they are calculated in various approaches: the constitute quark model and light cone quark model [36–41], covariant light front approach (LFQM) [30, 42, 43], light-cone sum rules [33, 44–66], PQCD [67–76] and lattice QCD [77–79]. We combine these results and use the average of the transition form factors of B meson decays at
q2 = 0, shown in Table 2. Theq2 dependence of the transition form factors of B meson decays are described with the dipole parametrization,FB→π0 FB→K0 FB→ηq0 FB→π1 FB→K1 FB→ηq1 AB→ρ0 AB→ω0 AB→K∗0 Fi(0) 0.28±0.03 0.31±0.03 0.21±0.02 0.28±0.03 0.31±0.03 0.21 \pm 0.02 0.36 \pm 0.04 0.32 \pm 0.03 0.39\pm 0.04 \alpha_1 0.50 0.53 0.52 0.52 0.54 1.43 1.56 1.60 1.51 \alpha_2 -0.13 -0.13 0 0.45 0.50 0.41 0.17 0.22 0.14 Table 2. The transition form factors of B meson decays at
q^{2} = 0 and dipole model parameters.\begin{split} F_{i}(q^{2}) = {F_{i}(0)\over 1-\alpha_{1}\displaystyle{q^{2}\over M_{\rm pole}^{2}}+ \alpha_{2}\displaystyle{q^{4}\over M_{\rm pole}^{4}}}, \end{split}
(9) where
F_{i} denotes the form factorsF_{0} ,F_{1} , orA_{0} , whileM_{\rm pole} is the mass of the corresponding pole state, such asB_{(s)} forA_{0} , andB^{*}_{(s)} forF_{0,1} .q^2 of charmless B meson decays is not far from zero, thus the uncertainties of the dipole model parameters are neglected in our calculations. The dipole model parameters are also listed in Table 2.To minimize the statistical uncertainties, we use the maximum number of experimental observables of the
B \to PP ,PV decays. However, some of them were measured with very poor precision. Hence, we do not use in our fit those measurements with less than a3\sigma significance. As a result, we have a total of 37 branching ratios and 11 CP violation observations for theB \to PP ,PV decays shown in Table 3 and Table 4, respectively, taken from the experimental data in Ref. [29], and the 2019 by averaging the measurements from LHCb [80–82], BELLE [83] and BABAR [84] and so on.Mode Amplitudes {\rm {Exp}} Mode Amplitudes {\rm {Exp}} \pi^{-} \pi^{0} T,C,P_{EW} 5.5\pm0.4 \pi^{+} \rho^{-} T,E,P,(P_E),P_A 14.6\pm1.6 \pi^{-} \eta T,C,P,P_C,P_{EW} 4.02\pm0.27 \pi^{-} \rho^{+} T,E,P,(P_E) 8.4\pm1.1 \pi^{-} \eta^{'} T,C,P,P_C,P_{EW} 2.7\pm0.9 \pi^{0} \rho^{0} C,C^{\prime},E,P,P_A,(P_E),P_{EW} 2\pm0.5 \pi^{+} \pi^{-} T,E,(P_E),P 5.12\pm0.19 \pi^{-} \bar{K^{*0}} P,P_A 10.1\pm0.9 \pi^{0} \pi^{0} C,E,P,(P_E),P_{EW} 1.59\pm0.26 \pi^{0} K^{*-} T,C,P,P_A,P_{EW} 8.2\pm1.9 K^{-} K^{0} P 1.31\pm0.17 \eta K^{*-} T,C,P,P_C,P_A,P_{EW} 19.3\pm1.6 K^{0} \bar{K^{0}} P 1.21\pm0.16 K^{-} \rho^{0} T,C^{\prime},P,P_{EW} 3.7\pm0.5 \pi^{-} \bar{K^{0}} P 23.7\pm0.8 K^{-} \omega T,C^{\prime},P,P_C^{\prime},P_{EW} 6.5\pm0.4 \pi^{0} K^{-} T,C,P,P_{EW} 12.9\pm0.5 K^{-} \phi P,P_C^{\prime},P_A,P_{EW} 8.8^{+0.7}_{-0.6} \eta K^{-} T,C,P,P_C,P_{EW} 2.4\pm0.4 \bar{K^{0}} \rho^{-} P 7.3^{+1.0}_{-1.2} \eta^{'} K^{-} T,C,P,P_C,P_{EW} 70.6\pm2.5 \pi^{+} K^{*-} T,P,P_A 8.4\pm0.8 \pi^{+} K^{-} T,P 19.6\pm0.5 \pi^{0} \bar{K^{*0}} C,P,P_A,P_{EW} 3.3\pm0.6 \pi^{0} \bar{K^{0}} C,P,P_{EW} 9.9\pm0.5 \eta \bar{K^{*0}} C,P,P_C,P_A,P_{EW} 15.9\pm1 \eta \bar{K^{0}} C,P,P_C,P_{EW} 1.23^{+0.27}_{-0.24} \eta^{'} \bar{K^{*0}} C,P,P_C,P_C^{\prime},P_A,P_{EW} 2.8\pm0.6 \eta^{'} \bar{K^{0}} C,P,P_C,P_{EW} 66\pm4 K^{-} \rho^{+} T,P 7\pm0.9 \pi^{-} \rho^{0} T,C^{\prime},P,P_A,P_{EW} 8.3\pm1.2 \bar{K^{0}} \rho^{0} C^{\prime},P,P_{EW} 4.7\pm0.6 \pi^{-} \omega T,C^{\prime},P,P_C^{\prime},P_A,P_{EW} 6.9\pm0.5 \bar{K^{0}} \omega C^{\prime},P,P_C^{\prime},P_{EW} 4.8\pm0.4 \pi^{0} \rho^{-} T,C,P,P_A,P_{EW} 10.9\pm1.4 \bar{K^{0}} \phi P,P_C^{\prime},P_A,P_{EW} 7.3\pm0.7 \eta^{'} \rho^{-} T,C,P,P_C,P_A,P_{EW} 9.7\pm2.2 Table 3. Branching fractions
(\times10^{-6}) of the various\bar B \to PP, PV decay modes with the decay amplitudes listed in the second column.Mode {{\cal{A}}_{\rm {exp}} } \mathcal{S_{\rm {exp}} } Mode {{\cal{A}}_{\rm {exp}} } \mathcal{S_{\rm {exp}} } \pi^{+} \pi^{-} 0.31±0.05 −0.67±0.06 K_s \omega 0.7±0.21 \pi^{0} K_s 0.58±0.17 K_s \phi 0.59±0.14 \eta^{'} K_s 0.63±0.06 K^{-} \rho^{0} 0.37±0.10 \eta K^{-} −0.37±0.08 \pi^{+} K^{*-} −0.22±0.06 \pi^{+} K^{-} −0.082±0.006 \eta \bar{K^{*0}} 0.19±0.05 Table 4. The direct CP asymmetries (
{{\cal{A}}_{\rm {exp}} } ) and the mixing-induced CP asymmetries (\mathcal{S_{\rm {exp}} } ) in the\bar B\to PP, PV decays.We use the
\chi^2 fit method of the Minuit program [85], where the\chi^2 function is defined in terms of n experimental observablesx_i \pm \Delta x_i and the corresponding theoretical valuesx_i^{\mathrm{th}} as\begin{split} \chi^{2} = \sum\limits_{i = 1}^{n}\left(\frac{x_i^{\rm th}-x_i}{\Delta x_i}\right)^2 . \end{split}
(10) The theoretical values are written as functions of the 15 theoretical parameters of the topological amplitudes. The best fit parameters are
\begin{array}{clcl} \gamma &\hspace{-3mm}= (69.8 \pm 2.1)^{ \circ} , \\[1mm] \chi^{C} &\hspace{-3mm}= 0.41 \pm 0.06,\; \; \; &\phi^{C} & \hspace{-3mm}= -1.74 \pm 0.09, \\[1mm] \chi^{C^{\prime}} &\hspace{-3mm}= 0.40 \pm 0.17,\; \; \; & \phi^{C^{\prime}} &\hspace{-3mm} = 1.78\pm 0.10, \\[1mm] \chi^{E} &\hspace{-3mm}= 0.06\pm0.006,\; \; \; &\phi^{E} &\hspace{-3mm} = 2.76\pm 0.13,\\[1mm] \chi^{P} &\hspace{-3mm}= 0.09\pm0.003,\; \; \;& \phi^{P} &\hspace{-3mm} = 2.55\pm 0.03 ,\\[1mm] \chi^{P_C}&\hspace{-3mm}= 0.045 \pm 0.003,\; \; \; &\phi^{P_C} &\hspace{-3mm} = 1.53 \pm 0.08, \\[1mm] \chi^{P_C^{\prime}} &\hspace{-3mm}= 0.037\pm 0.003,\; \; \; & \phi^{P_C^{\prime}} &\hspace{-3mm} = 0.67 \pm 0.08, \\[1mm] \chi^{P_A} &\hspace{-3mm}= 0.006\pm0.0008,\; \; \; & \phi^{P_A} &\hspace{-3mm} = 1.49\pm 0.09, \end{array}
(11) with
\chi^{2}/\text{d.o.f} = 45.4/33 = 1.4 . The uncertainties of the parameters shown above are given by the\chi^2 fit program Minuit, and are mainly due to the experimental statistical and systematic uncertainties. There are also theoretical uncertainties of the extracted\gamma values. The major source of theoretical uncertainties in our calculations are the uncertainties of the input parameters:|V_{ub}| ,|V_{cb}| , hadronic form factors and decay constants. We repeated the fit by varying the input parameters one by one assuming their Gaussian distribution. The theoretical uncertainties of\gamma were then assessed using the distribution of the central values of\gamma form these fits. The values of the theoretical uncertainty\sigma_{\mathrm{(T.)}} obtained are0.2{^\circ} ,0.2^{\circ} ,0.9^{\circ} and0.1^{\circ} from the uncertainties of|V_{ub}| ,|V_{cb}| , form factors and decay constants, respectively. The total theoretical uncertainty is0.9^{\circ} . Our final result for\gamma is then(69.8 \pm 2.1 \pm 0.9) ^{\circ} , which is in good agreement with the current world averages:\gamma = (71.1^{+4.6}_{-5.3})^{\circ} [5],\gamma = (73.5^{+4.2}_{-5.1})^{\circ} [6] and\gamma = (70.0 \pm 4.2)^{\circ} [7], and the latest measurement of(74.0^{+5.0}_{-5.8})^{\circ} by the LHCb collaboration [8]. It is obvious that our uncertainties of\gamma are reduced to roughly a half of the uncertainties of the world-average values. -
The charmless B meson decays receive contributions from the tree amplitudes and the loop penguin amplitudes. The interference between these amplitudes makes the branching ratios of these decay modes sensitive to the CKM angles, where large direct CP asymmetries are observed. Since non-perturbative dynamics is involved, the hadronic matrix elements of these decays are difficult to calculate precisely. We parametrized the decay amplitudes of the topological diagrams, which can be fitted by the experimentally measured quantities, such as the branching ratios and CP asymmetry parameters. To improve the precision of the global fit, we factorized the corresponding decay constants and form factors in order to characterize the flavor SU(3) breaking effect. We extracted the CKM weak angle
\gamma using the measured two-body charmlessB\to PP ,PV decays in the factorization assisted topological amplitude approach. The value obtained is(69.8 \pm 2.1 \pm 0.9) ^{\circ} , where the first uncertainty is from the experimental errors and the second from the hadronic parameters and CKM matrix elements. This result is compatible with the current world average values and the latest measurement of(74.0^{+5.0}_{-5.8})^{\circ} by the LHCb collaboration, but has a smaller uncertainty.We are grateful to Wen-Bin Qian and Xiao-Rui Lyu for useful discussion.
