Emerging $\Delta M_{d}$-Anomaly from Tree-Level Determinations of $|V_{cb}|$ and the Angle $\gamma$

TL;DR

Recent measurements of the CKM angle gamma and |V_{cb}| suggest a 30% enhancement in ΔM_d, challenging existing flavor models and indicating potential new physics contributions.

Contribution

The paper highlights a tension between experimental CKM parameters and ΔM_{d,s} predictions, proposing a new approach using gamma, beta, |V_{cb}|, and |V_{us}| as fundamental parameters.

Findings

ΔM_d shows a 30% enhancement over data.

Tension with CMFV and U(2)^3 models.

Implications for |V_{ub}|, |V_{td}|, and rare decays.

Abstract

We point out that the recently increased value of the angle $\gamma$ in the Unitarity Triangle (UT), determined in tree-level decays to be $\gamma=(74.0^{+5.0}_{-5.8})^\circ$ by the LHCb collaboration, combined with the most recent value of $|V_{cb}|$ implies an enhancement of $\Delta M_{d}$ over the data in the ballpark of $30\%$. This is roughly by a factor of two larger than the enhancement of $\Delta M_{s}$ that is independent of $\gamma$. This disparity of enhancements is problematic for models with Constrained Minimal Flavour Violation (CMFV) and also for $U(2)^3$ models. In view of the prospects of measuring $\gamma$ with the precision of $\pm 1^\circ$ by Belle II and LHCb in the coming years, we propose to use the angles $\gamma$ and $\beta$ together with $|V_{cb}|$ and $|V_{us}|$ as the fundamental parameters of the CKM matrix until $|V_{ub}|$ from tree-level decays will be known precisely. Displaying $\Delta M_{s,d}$ as functions of $\gamma$ clearly demonstrates the tension between the value of $\gamma$ from tree-level decays, free from new physics (NP) contributions, and $\Delta M_{s,d}$ calculated in CMFV and $U(2)^3$ models and thus exhibits the presence of NP contributions to $\Delta M_{s,d}$ beyond these frameworks. We calculate the values of $|V_{ub}|$ and $|V_{td}|$ as functions of $\gamma$ and $|V_{cb}|$ and discuss the implications of our results for $\varepsilon_K$ and rare $K$ and $B$ decays. We also briefly discuss a future strategy in which $\beta$, possibly affected by NP, is replaced by $|V_{ub}|$.

Tables

Table 1. Table 1: Values of the experimental and theoretical quantities used as input parameters. For future updates see PDG [ 22 ] and HFLAV [ 6 ] .

$m_{B_s}=5366.8\left(2\right)\mathrm{MeV}$ [22]	$m_{B_d}=5279.58\left(17\right)\mathrm{MeV}$ [22]
$\Delta M_s=17.757\left(21\right){\text{ps}}^{-1}$ [6]	$\Delta M_d=0.5055\left(20\right){\text{ps}}^{-1}$ [6]
$S_{\psi K_S}=0.691\left(17\right)$ [6]	$S_{\psi \varphi }=0.015\left(35\right)$ [6]
$\left|V_{us}\right|=0.2253\left(8\right)$ [22]	$\left|{\epsilon }_K\right|=2.228\left(11\right)\cdot {10}^{-3}$ [22]
$F_{B_s}$ = $228.6\left(3.8\right)\mathrm{MeV}$ [23]	$F_{B_d}$ = $193.6\left(4.2\right)\mathrm{MeV}$ [23]
$m_t\left(m_t\right)=163.53\left(85\right)\mathrm{GeV}$	$S_0\left(x_t\right)=2.322\left(18\right)$
${\eta }_{cc}=1.87\left(76\right)$ [24]	${\eta }_{ct}=0.496\left(47\right)$ [25]
${\eta }_{tt}=0.5765\left(65\right)$ [21]	${\eta }_B=0.55\left(1\right)$ [21, 26]
${\tau }_{B_s}=1.510\left(5\right)\text{ps}$ [6]	$\Delta {\Gamma }_s/{\Gamma }_s=0.124\left(9\right)$ [6]
${\tau }_{B_d}=1.520\left(4\right)\text{ps}$ [6]	${\kappa }_{\epsilon }=0.94\left(2\right)$ [27]

Table 2. Table 2: Left: Central values for the branching ratio ℬ ​ ( K + → π + ​ ν ​ ν ¯ ) ℬ → superscript 𝐾 superscript 𝜋 𝜈 ¯ 𝜈 \mathcal{B}(K^{+}\to\pi^{+}\nu\bar{\nu}) for various values of γ 𝛾 \gamma and | V c ​ b | subscript 𝑉 𝑐 𝑏 |V_{cb}| . The angle β 𝛽 \beta is fixed to β = 21.85 ∘ 𝛽 superscript 21.85 \beta=21.85^{\circ} determined from S ψ ​ K S subscript 𝑆 𝜓 subscript 𝐾 𝑆 S_{\psi K_{S}} . Right: Central values for the branching ratio ℬ ​ ( K L → π 0 ​ ν ​ ν ¯ ) ℬ → subscript 𝐾 𝐿 superscript 𝜋 0 𝜈 ¯ 𝜈 \mathcal{B}(K_{L}\to\pi^{0}\nu\bar{\nu}) for various values of β 𝛽 \beta , γ 𝛾 \gamma and | V c ​ b | subscript 𝑉 𝑐 𝑏 |V_{cb}| .

$\gamma {[}^{\circ }]$	${10}^3\cdot \left|V_{cb}\right|$
	39	40	41	42	43
64	6.7	7.2	7.8	8.3	8.9
66	6.9	7.4	7.9	8.5	9.1
68	7.1	7.6	8.1	8.7	9.3
70	7.2	7.7	8.3	8.9	9.5
72	7.4	7.9	8.5	9.1	9.7
74	7.5	8.1	8.6	9.2	9.9
76	7.7	8.2	8.8	9.4	10.0
78	7.8	8.4	9.0	9.6	10.3

Equations

$$\Delta M_s,\Delta M_d,S_{\psi K_S},S_{\psi \varphi },{\epsilon }_K$$

$$|V_{us}|,|V_{cb}|,|V_{ub}|,\gamma ,\beta ,$$

$$(\beta ,\gamma ),(R_b,\gamma ),$$

$$\gamma =(74.0_{-5.8}^{+5.0}{)}^{\circ },.$$

$$|V_{cb}|=(42.0\pm 0.6)\times 10^{-3},$$

$$|V_{us}|,|V_{cb}|,\gamma ,\beta .$$

$$F_{B_s}\sqrt{{\hat{B}}_{B_s}},F_{B_d}\sqrt{{\hat{B}}_{B_d}}.$$

$$F_{B_s}\sqrt{{\hat{B}}_{B_s}}=(274.6\pm 8.8)\mathrm{MeV},F_{B_d}\sqrt{{\hat{B}}_{B_d}}=(227.7\pm 9.8)\mathrm{MeV},$$

$$\xi =\frac{F_{B_s}\sqrt{{\hat{B}}_{B_s}}}{F_{B_d}\sqrt{{\hat{B}}_{B_d}}}=1.206\pm 0.019.$$

$$\gamma =(63.0\pm 2.1{)}^{\circ },$$

$$\beta =(21.85\pm 0.67{)}^{\circ },\gamma =(74.0_{-5.8}^{+5.0}{)}^{\circ },$$

$$S_{\psi K_S}=\sin 2\beta =0.691\pm 0.017.$$

$$R_b=\frac{\sin (\beta )}{\sin (\gamma +\beta )}=0.374\pm 0.012,R_t=\frac{\sin (\gamma )}{\sin (\gamma +\beta )}=0.964\pm 0.035.$$

$$V_{td}=|V_{td}|e^{-i\beta },V_{ub}=|V_{ub}|e^{-i\gamma }.$$

$$|V_{td}|=|V_{us}||V_{cb}|R_t,|V_{ts}|={\eta }_R|V_{cb}|$$

$${\eta }_R=1-|V_{us}|\xi \sqrt{\frac{\Delta M_d}{\Delta M_s}}\sqrt{\frac{m_{B_s}}{m_{B_d}}}\cos \beta +\frac{{\lambda }^2}{2}+\mathcal{O}({\lambda }^4)=0.9825,$$

$$|V_{ub}|=\lambda |V_{cb}|\frac{R_b}{1-\frac{{\lambda }^2}{2}}.$$

$$\Delta M_d$$

$$\Delta M_s$$

$$S_0(x_t)=\frac{4x_t-11x_t^2+x_t^3}{4(1-x_t{)}^2}-\frac{3x_t^2\log x_t}{2(1-x_t{)}^3}=2.322{\left[\frac{{\overline{m}}_{\mathrm{t}}(m_{\mathrm{t}})}{163.5\mathrm{GeV}}\right]}^{1.52},$$

$$(\Delta M_d{)}_{\text{SM}}=(0.648\pm 0.077){\text{ps}}^{-1},(\Delta M_s{)}_{\text{SM}}=(19.8\pm 1.4){\text{ps}}^{-1},$$

$$S(v)\ge S_0(x_t)=2.322.$$

$$\frac{|V_{td}|}{|V_{ts}|}=\xi \sqrt{\frac{m_{B_s}}{m_{B_d}}}\sqrt{\frac{\Delta M_d}{\Delta M_s}}$$

$$\frac{|V_{td}|}{|V_{ts}|}=\xi \sqrt{\frac{m_{B_s}}{m_{B_d}}}\sqrt{\frac{(\Delta M_d{)}_{\exp }}{(\Delta M_s{)}_{\exp }}}=0.2052\pm 0.0033,$$

$$|{\epsilon }_K{|}_{\mathrm{SM}}=(2.26\pm 0.27)\cdot 10^{-3},$$

$$R_b=\left(1-\frac{{\lambda }^2}{2}\right)\frac{1}{\lambda }\left|\frac{V_{ub}}{V_{cb}}\right|.$$

$$R_t=\sqrt{1+R_b^2-2R_b\cos \gamma },\cot \beta =\frac{1-R_b\cos \gamma }{R_b\sin \gamma },$$

$$S_{\psi K_S}=\sin (2\beta +2{\phi }_{\mathrm{new}})=0.691\pm 0.017,$$

$${\phi }_{\mathrm{new}}\simeq -2.2^{\circ }.$$

$$S_i=S_0(x_t)-\Delta S_i{e}^{i{\delta }_i}(i=K,s,d),$$

$$\Delta S_s<\Delta S_d.$$

Full text

TTP18-044

AJB-18-10

**Emerging $\Delta M_{d}$-Anomaly from Tree-Level

Determinations of $|V_{cb}|$ and the Angle $\gamma$**

Monika Blankea,b and Andrzej J. Burasc,d a Institut für Kernphysik, Karlsruhe Institute of Technology, Hermann-von-Helmholtz-Platz 1, D-76344 Eggenstein-Leopoldshafen, Germany

b Institut für Theoretische Teilchenphysik, Karlsruhe Institute of Technology, Engesserstraße 7, D-76128 Karlsruhe, Germany

c TUM-IAS, Lichtenbergstr. 2a, D-85748 Garching, Germany

dPhysik Department, TUM, D-85748 Garching, Germany

**Abstract

**

We point out that the recently increased value of the angle $\gamma$ in the Unitarity Triangle (UT), determined in tree-level decays to be $\gamma=(74.0^{+5.0}_{-5.8})^{\circ}$ by the LHCb collaboration, combined with the most recent value of $|V_{cb}|$ implies an enhancement of $\Delta M_{d}$ over the data in the ballpark of $30\%$. This is larger by roughly a factor of two than the enhancement of $\Delta M_{s}$ that is independent of $\gamma$. This disparity of enhancements is problematic for models with Constrained Minimal Flavour Violation (CMFV) and also for $U(2)^{3}$ models. In view of the prospects of measuring $\gamma$ with the precision of $\pm 1^{\circ}$ by Belle II and LHCb in the coming years, we propose to use the angles $\gamma$ and $\beta$ together with $|V_{cb}|$ and $|V_{us}|$ as the fundamental parameters of the CKM matrix until $|V_{ub}|$ from tree-level decays will be known precisely. Displaying $\Delta M_{s,d}$ as functions of $\gamma$ clearly demonstrates the tension between the value of $\gamma$ from tree-level decays, free from new physics (NP) contributions, and $\Delta M_{s,d}$ calculated in CMFV and $U(2)^{3}$ models and thus exhibits the presence of NP contributions to $\Delta M_{s,d}$ beyond these frameworks. We calculate the values of $|V_{ub}|$ and $|V_{td}|$ as functions of $\gamma$ and $|V_{cb}|$ and discuss the implications of our results for $\varepsilon_{K}$ and rare $K$ and $B$ decays. We also briefly discuss a future strategy in which $\beta$, possibly affected by NP, is replaced by $|V_{ub}|$.

1 Introduction

The $\Delta F=2$ transitions in the down-quark sector, that is $B^{0}_{s,d}-\bar{B}^{0}_{s,d}$ and $K^{0}-\bar{K}^{0}$ mixings, have been vital in constraining the Standard Model (SM) and in the search for new physics (NP) for several decades[1, 2]. However, theoretical uncertainties related to the hadronic matrix elements entering these transitions and their large sensitivity to the CKM parameters made clear cut conclusions about the presence of NP impossible. As we demonstrate in this paper, this could change in the near future.

Among the most important flavour observables we have at our disposal are

[TABLE]

with $\Delta M_{s,d}$ being the mass differences in $B^{0}_{s,d}-\bar{B}^{0}_{s,d}$ mixings and $S_{\psi K_{S}}$ and $S_{\psi\phi}$ the corresponding mixing induced CP-asymmetries. $\varepsilon_{K}$ describes the magnitude of indirect CP-violation in $K^{0}-\bar{K}^{0}$ mixing. $\Delta M_{s,d}$ and $\varepsilon_{K}$ are already known experimentally with impressive precision. The asymmetries $S_{\psi K_{S}}$ and $S_{\psi\phi}$ are less precisely measured but have the advantage of being subject to only very small hadronic uncertainties.

On the other hand the CKM parameters of particular interest are

[TABLE]

with the first three being the moduli of the most intensively studied elements of the CKM matrix, and $\gamma$ and $\beta$ being two angles in the Unitarity Triangle (UT). The angle $\gamma$ is to an excellent approximation equal to the sole complex phase in the standard parametrization of the CKM matrix.

Now, as elaborated in [3], there are many ways to construct the rescaled UT. They all involve only two inputs, but as quantified in the latter paper, some pairs are particularly suited for the determination of the apex ($\bar{\rho},\bar{\eta}$) of this triangle, as only moderate precision on them is required to obtain a satisfactory determination of $\bar{\rho}$ and $\bar{\eta}$. The clear winners from this study are the pairs

[TABLE]

with $R_{b}$ being the length of one side in the UT related to the ratio $|V_{ub}|/|V_{cb}|$.

Ideally, one would like to use the second pair which allows to construct the so-called reference unitarity triangle (RUT) [4] that is supposed to be free of NP contributions. Unfortunately, the persistent discrepancy between inclusive and exclusive determinations of $|V_{ub}|$ from tree-level decays precludes a satisfactory determination of the RUT at present.

On the other hand the tree-level determination of the angle $\gamma$ has significantly been improved in the last years by various measurements of the LHCb collaboration, with the latest average being [5]111The HFLAV average $\gamma=(73.5^{+4.2}_{-5.1})^{\circ}$ [6] does not include the latest LHCb result.

[TABLE]

Moreover, the prospects of LHCb and Belle II [7, 8] to decrease the error down to $\pm 1^{\circ}$ are promising. In view of this situation and significant recent progress in the determination of $|V_{cb}|$, giving [9]

[TABLE]

we will choose as the four fundamental CKM parameters

[TABLE]

Within the SM and CMFV models [10, 11, 12], the hadronic uncertainties in $\Delta M_{s,d}$ reside within a good approximation in the parameters

[TABLE]

Fortunately, during the last years their uncertainties decreased significantly. In particular, an impressive progress has been made by the Fermilab Lattice and MILC Collaborations (Fermilab-MILC) that find [13]

[TABLE]

with uncertainties of $3\%$ and $4\%$, respectively. An even higher precision is achieved for the ratio

[TABLE]

Based on the results in (8) and (9) we have performed in [14] a detailed analysis of $\Delta F=2$ processes in CMFV models, finding a significant tension between $\Delta M_{s,d}$ and $\varepsilon_{K}$ in these models with the pattern of the tension strongly dependent on the value of $|V_{cb}|$. Moreover, constructing the universal unitarity triangle (UUT) [10] via $R_{t}$ and $\beta$ we could predict, independently of $|V_{cb}|$, the value of $\gamma$ to be

[TABLE]

significantly below the value in (4). This number has not changed with respect to our 2016 analysis, and now displays a $1.8\sigma$ tension with the improved tree level measurement in (4).222Using instead of (9) the very recent RBC-UKQCD result $\xi=1.1853\pm 0.0054^{+0.0116}_{-0.0156}$ [15], one obtains $\gamma=(60.7\pm 1.5)^{\circ}$, increasing the tension to $2.2\sigma$. As we have discussed in [14] this problem arises not only in the SM and more generally in CMFV models but also in minimally broken $U(2)^{3}$ models, where NP contributions in the $B_{d}$ and $B_{s}$ systems are universal and hence cancel in the ratio.

As the present paper deals again with the tensions between $\Delta F=2$ observables in CMFV models, it is mandatory for us to state what is new in our paper:

• •

In [14], we have considered two strategies. One in which $\varepsilon_{K}$ has been used to determine $|V_{cb}|$, implying a value consistent with the inclusive determination as well as $\Delta M_{s,d}$ values well above the data. In the second strategy, $|V_{cb}|$ has been determined from $\Delta M_{s}$ resulting in a low value of $|V_{cb}|$ consistent with the exclusive determination at that time. The predicted $\varepsilon_{K}$ then turned out to be well below its experimental value. The recent improvements in the determinations of $|V_{cb}|$ [16, 17] disfavours the second strategy and also the recent claim in [18] that there is a $4\sigma$ anomaly in $\varepsilon_{K}$.

• •

More importantly, in view of the improved value of $\gamma$, we decided to use it as an input in the present analysis, instead of the usual determination of the UT in CMFV models through $S_{\psi K_{S}}~{}(\beta)$ and the side $R_{t}$ of the UT determined from the ratio $\Delta M_{d}/\Delta M_{s}$ and $\xi$ in (9).

• •

The most recent discussions, see in particular [19, 20], dealt exclusively with the implications of the enhanced value of $\Delta M_{s}$ and not $\Delta M_{d}$, for which in addition to the increased value of $|V_{cb}|$ also the increased value of $\gamma$ matters.

In the context of the second item we remark that the $(R_{t},\beta)$ strategy for the determination of the UT has been found in [3] to be less powerful than the $(\beta,\gamma)$ strategy used here. Moreover, as NP now is expected in $\Delta M_{s,d}$, it appears as a better strategy to replace their ratio by the angle $\gamma$ and instead treat $\Delta M_{s,d}$ as outputs being functions of $\gamma$, $\beta$ and $|V_{cb}|$.

One could wonder why the emerging $\Delta M_{d}$ anomaly pointed out by us has not been noticed in the global fits performed by the CKMfitter and UTfit collaborations. In our view such global fits, involving simultaneously many quantities, are likely to miss NP effects present in only a subset of observables, in particular when the significance has not reached the discovery level. We are optimistic that the findings of this paper pointing towards NP in the $B_{d}$ system will motivate both theorists and experimentalists to intensify the search for NP in $b\to d$ transitions, after the last five years being dominated by the study of $b\to s$ and $b\to c$ transitions.

Our paper is organized as follows. In Section 2 we present the determination of the UT and of the CKM matrix using the $(\beta,\gamma)$ strategy. In Section 3 we evaluate $\Delta M_{d}$ and $\Delta M_{s}$ as functions of $\gamma$, finding their values to disagree with the data. The new result relative to [14] and other recent papers [19, 20] is the disagreement of $\Delta M_{d}$ and the ratio $\Delta M_{d}/\Delta M_{s}$ with the data, a direct consequence of the increased value of $\gamma$. On the other hand $\varepsilon_{K}$ agrees well with the data. We therefore provide the SM predictions for the branching ratios of $K^{+}\rightarrow\pi^{+}\nu\bar{\nu}$ and $K_{L}\rightarrow\pi^{0}\nu\bar{\nu}$ for different values of $\gamma$, $\beta$ and $|V_{cb}|$. In Section 4 we have a look at the $(R_{b},\gamma)$ strategy, which could become favourable in the next decade, once the tree-level determination of $|V_{ub}|$ is settled. In Section 5 we briefly investigate what kind of NP could be responsible for the $\Delta M_{s,d}$ anomalies found in Section 3 and what are the implications for NP in $\Delta F=1$ transitions. We conclude in Section 6.

2 Deriving the UT and the CKM matrix

Our determination of the UT and of the CKM matrix proceeds in two steps:

Step 1:

We use as input parameters

[TABLE]

with $\beta$ obtained from

[TABLE]

This allows us to determine the two sides $R_{b}$ and $R_{t}$ of the UT shown in Fig. 1, that are given in terms of $\beta$ and $\gamma$ as follows [3]

[TABLE]

The angles $\beta$ and $\gamma$ of the unitarity triangle are directly related to the complex phases of the CKM-elements $V_{td}$ and $V_{ub}$, respectively, through

[TABLE]

Step 2:

Including $\lambda\equiv|V_{us}|$ and $|V_{cb}|$ as the remaining input parameters we determine $|V_{td}|$ and $|V_{ts}|$ through

[TABLE]

with

[TABLE]

where we have used $\beta$ in (11). $\Delta M_{d,s}$ are taken from experiment as given in Table 1 but using our $|V_{cb}|$ and $\gamma$ dependent values would change the result by less than $1\%$.

Finally, we find

[TABLE]

In Fig. 2 we show the constraints on the UT from the tree-level measurement of $\gamma$, from $\beta$ extracted from $S_{\psi K_{S}}$, and $R_{t}$ from $\Delta M_{d}/\Delta M_{s}$. The advantage of the $(\gamma,\beta)$ strategy over the $(R_{t},\beta)$ strategy is not seen yet because of a significant error in $\gamma$. With the future uncertainty on $\gamma$ of $\pm 1^{\circ}$ represented by the black area, the power of the $(\gamma,\beta)$ strategy in determining the UT is clearly visible. However, already now we observe that the apex of the UT obtained from the $(\gamma,\beta)$ strategy disagrees with the one from the $(R_{t},\beta)$ one. This tension indicates the presence of some NP contributions.

In Fig. 3 we show $|V_{td}|$ as a function of $\gamma$ and $|V_{ub}|$ as a function of $\beta$ for different values of $|V_{cb}|$. The dependences of $|V_{td}|$ on $\beta$ and of $|V_{ub}|$ on $\gamma$ are very small. These plots will allow to monitor the values of $|V_{td}|$ and $|V_{ub}|$ that enter various observables as the uncertainties of $\gamma$, $\beta$ and $|V_{cb}|$ will shrink with time.

3 Calculating observables

For the mass differences in the $B^{0}_{s,d}-\bar{B}^{0}_{s,d}$ systems we have the very accurate expressions [14]

[TABLE]

Here $S(v)$ is the box-function in CMFV models with $v$ denoting parameters of a given model including $x_{t}=m_{t}^{2}/M_{W}^{2}$. The value $2.322$ in the normalization of $S(v)$ is its SM value for $m_{t}(m_{t})=163.5\,{\rm GeV}$ obtained from

[TABLE]

and $\eta_{B}$ is the perturbative QCD correction [21]. Our input parameters, equal to the ones used in [13], are collected in Table 1. We find

[TABLE]

which differ from the experimental values by $1.9\sigma$ and $1.4\sigma$, respectively. As the correlation matrix of the relevant lattice parameters entering these predictions is unknown to us, we do not attempt to derive a global significance for the anomaly in $B_{d,s}-\bar{B}_{d,s}$ mixing.

Now, the overall factors in (18) and (19) are the central experimental values, and in CMFV models $S(v)$ is bounded from below by its SM value [28]

[TABLE]

Consequently, with the values of $|V_{td}|$ found in the previous section, that are significantly larger than its nominal value in (18), it is evident that CMFV models have difficulties in describing the data for $\Delta M_{d}$. In addition, with the value of $|V_{cb}|$ in (5) also $|V_{ts}|$ is significantly larger than its nominal value in (19). Therefore $\Delta M_{s}$ in CMFV models is enhanced over its experimental value as already pointed out in [13, 14] and recently analysed in [19, 20]. Yet the latter enhancement is not as large as for $\Delta M_{d}$ because $\Delta M_{s}$ does not depend on $\gamma$.

In Fig. 4 we show in the left panel $\Delta M_{d}$ and $\Delta M_{s}$ normalized to their experimental values. Evidently, for central values of all parameters, $\Delta M_{d}$ differs by roughly $30\%$ from the data while in the case of $\Delta M_{s}$ the corresponding difference amounts only to $12\%$. But the uncertainties in other parameters like $|V_{cb}|$ and the hadronic parameters in (8) are still significant. However, we expect that in the coming years these uncertainties will significantly be reduced.

In the right panel of Fig. 4 we show the ratio $\Delta M_{s}/\Delta M_{d}$ as a function of $\gamma$. The dependence on $|V_{cb}|$ cancels in this ratio and the error on $\xi$ in (9) is much smaller than the errors in (8). Consequently the disagreement of the ratio in question with the data, shown as a horizontal line at $35.1$, is clearly visible and expresses the problem of CMFV models and those based on the $U(2)^{3}$ symmetry.

Of interest is also the ratio

[TABLE]

with $\Delta M_{s,d}$ predicted here to be compared with

[TABLE]

with $\Delta M_{s,d}$ taken from experiment. In CMFV models and those with $U(2)^{3}$ symmetry this ratio depends only on the angle $\gamma$. We show this in Fig. 5. A significant enhancement of $|V_{td}|/|V_{ts}|$ over the value in (24) is observed.

As far as $\varepsilon_{K}$ is concerned, using the standard expression as given e. g. in [2] and all input parameters collected in Table 1, we find the SM value of $\varepsilon_{K}$ to be fully consistent with the data:

[TABLE]

with higher values for the remaining CMFV models due to the bound in (22).

Despite the fact that the SM fails to describe the data for $\Delta M_{s,d}$, having determined the CKM parameters, the agreement of the SM with the experimental value for $\varepsilon_{K}$ invites us to calculate the branching ratios for $K^{+}\rightarrow\pi^{+}\nu\bar{\nu}$ and $K_{L}\rightarrow\pi^{0}\nu\bar{\nu}$ in the SM. This is of interest in view of the NA62 and KOTO experiments that should provide results for these decays in the coming years. Using the parametric formulae of [29] we find the central values of $\mathcal{B}(K^{+}\rightarrow\pi^{+}\nu\bar{\nu})$ and $\mathcal{B}(K_{L}\rightarrow\pi^{0}\nu\bar{\nu})$ given in table 2 for different values of $\gamma$, $\beta$ and $|V_{cb}|$.

4 $(R_{b},\gamma)$ strategy

It is likely that in the next decade the $(\beta,\gamma)$ strategy will be replaced by the $(R_{b},\gamma)$ strategy. This could turn out to be even necessary if the value of $|V_{ub}|$ determined from tree-level processes turned out to be very different from the one determined in the previous section. Therefore for completeness we want to give the relevant formulae for this strategy.

Knowing $|V_{ub}|$ determined in tree-level decays, one finds $R_{b}$ using

[TABLE]

Together with $\gamma$, this result allows to determine $R_{t}$ and $\beta$ by means of

[TABLE]

so that the RUT is completely fixed.

If the resulting value of $\beta$ differs from the one in (11), then the expression in (12) will have to be replaced by

[TABLE]

with $\varphi_{\rm new}$ being a new CP-violating phase. For instance for $|V_{ub}|=4.0\cdot 10^{-3}$ we find

[TABLE]

5 Going beyond CMFV

Our analysis signals the violation of flavour universality in the function $S(v)$, characteristic for CMFV models. It hints for the presence of new sources of flavour and CP-violation and/or new operators contributing to $\Delta F=2$ transitions beyond the SM $(V-A)\otimes(V-A)$ ones.333In a more general formulation of MFV new operators could be present [30]. For simplicity we restrict first our discussion of NP scenarios to the ones in which only SM operators are present.

A fully general and very convenient solution in this case is just to consider instead of the flavour universal function $S(v)$ three functions

[TABLE]

with $\Delta S_{i}$ being real and positive definite quantities, and the minus sign required to suppress $\Delta M_{d,s}$ below their SM values. It is evident that with two free parameters in each meson system it is always possible to obtain an agreement with the data on $\Delta F=2$ observables. Our analysis indicates the following pattern of these parameters:

• •

A clear breakdown of the universality of $S(v)$ with

[TABLE]

• •

The new phases

[TABLE]

in order not to spoil the good agreement of the SM with the experimental values of $S_{\psi K_{S}}$ and $S_{\psi\phi}$.

• •

In the case of $K^{0}-\bar{K}^{0}$ mixing, the good agreement of $\varepsilon_{K}$ with its measured value implies a small imaginary part of the NP contribution. This can either be achieved by a small value of $\Delta S_{K}$, or by an appropriately chosen value of the new phase $\delta_{K}$.

Note that the fate of $\delta_{d}$ and to a lesser extend of $\delta_{K}$ will depend on the future value of $|V_{ub}|$ as remarked in connection with (28).

This pattern cannot be explained in models with a minimally broken $U(2)^{3}$ flavour symmetry [31, 32, 33] in which the equality $\Delta S_{s}=\Delta S_{d}$ is predicted, although the near equality of $\delta_{s}$ and $\delta_{d}$ is a property of these models. This could change if for instance $|V_{ub}|$ was found significantly different from the value followed from our strategy, as shown in Fig. 6. But as these models fail anyway we will not consider them further.

The simplest models beyond the CMFV and $U(2)^{3}$ frameworks one could consider are models with tree-level $Z^{\prime}$ and $Z$ exchanges. While in [34, 35, 36, 37] general studies of such scenarios have been considered, specific examples are models with vector-like quarks [38] and 331 models [39]. These models have sufficient numbers of parameters to obtain an agreement with the data for $\Delta F=2$ processes. This is explicitly shown for the case of 331 models in [39].

The minus sign in (30) has been introduced by us by hand. Strictly speaking, as already discussed in the context of $\Delta M_{s}$ in [19], in the presence of only left-handed currents the minus sign in (30) truly requires the NP phases to be $\pi+\delta_{d,s}$. Following the reasoning in [40], this implies the CP-violating phases in the corresponding $\Delta F=1$ $b\to d,s$ transitions to be close to $\pi/2$, i. e. maximal. We hence conclude that, within models with only left-handed currents, the observed suppression of $\Delta M_{d}$ and to a lesser extent $\Delta M_{s}$ implies significant deviations from the SM in CP-asymmetries of radiative and rare $b\to d$ and $b\to s$ decays. As quantitative predictions for these observables are model-dependent, we leave their thorough analysis for future work.

In the presence of both left- and right-handed couplings, on the other hand, the suppression of $\Delta M_{d}$ is much easier to achieve without introducing large CP-violating phases. In this context probably most interesting are models in which the SMEFT operator $\mathcal{O}_{Hd}$ involving right-handed flavour violating couplings to down-quarks is generated at the NP scale. As demonstrated in [37], the renormalisation group evolution to low-energy scales involving also left-handed currents present already within the SM generates left-right $\Delta F=2$ operators representing FCNCs mediated by the $Z$ boson. At NLO this effect has also been discussed in [36]. An explicit realization of such a NP scenario is provided by models with vector-like quarks with an additional $U(1)$ gauge symmetry so that both tree-level $Z$ and $Z^{\prime}$ exchanges are present, and in some models of this type also box diagram contributions with vector-like quarks, Higgs and other scalar and pseudoscalar exchanges are important [38]. The test of these scenarios is then mainly offered through the correlations of $\Delta M_{d,s}$ with $\Delta F=1$ processes, that is rare $K$ or $B_{s,d}$ decays, the ratio $\varepsilon^{\prime}/\varepsilon$ and other observables. This is evident from the analyses in [38, 37] and once the data on $\gamma$, $|V_{cb}|$ and $|V_{ub}|$ improve, could be an arena for further investigation of the implications of the $\Delta M_{d}$ anomaly pointed out here.

6 Summary

The main message of our paper is the emerging $\Delta M_{d}$ anomaly which is significantly larger than the $\Delta M_{s}$ one discussed in [14, 19, 20]. Its fate will depend strongly on the improved values of $\gamma$ and $|V_{cb}|$ from tree-level decays and, to a lesser extent, on $|V_{ub}|$, which is more relevant for the prediction of $\sin 2\beta$ in the SM. This anomaly, if confirmed, will have implications for observables sensitive to $b\to d$ transitions like $b\to d\ell^{+}\ell^{-}$ and $b\to d\nu\bar{\nu}$ which will be explored by Belle II. It will open a new oasis of NP, analogous to the one related to the recent anomalies in $b\to s\ell^{+}\ell^{-}$ and their implications for $b\to s\nu\bar{\nu}$ transitions. Depending on the NP flavour structure, it could also have implications for $K^{+}\rightarrow\pi^{+}\nu\bar{\nu}$ and $K_{L}\rightarrow\pi^{0}\nu\bar{\nu}$.

Acknowledgements

The research of AJB was fully supported by the DFG cluster of excellence “Origin and Structure of the Universe”.