The CKM parameters in the SMEFT

Extraction of the Cabibbo-Kobayashi-Maskawa (CKM) matrix from flavour observables can be affected by physics beyond the Standard Model (SM). We provide a general roadmap to take this into account, which we apply to the case of the Standard Model Effective Field Theory (SMEFT). We choose a set of four input observables that determine the four Wolfenstein parameters, and discuss how the effects of dimension-six operators can be included in their definition. We provide numerical values and confidence intervals for the CKM parameters, and compare them with the results of CKM fits obtained in the SM context. Our approach allows one to perform general SMEFT analyses in a consistent fashion, independently of any assumptions about the way new physics affects flavour observables. We discuss a few examples illustrating how our approach can be implemented in practice.


Introduction
Quark flavour-changing transitions have been used to probe high-energy scales since several decades, with many successes along the way, such as predicting the existence and properties of the charm and top quarks. They are particularly well suited for this purpose in the context of the Standard Model (SM), where flavour transitions are controlled by the unitary Cabibbo-Kobayashi-Maskawa (CKM) matrix defined by only 4 real parameters. Over the last decades, the improvement in measurements has been matched by theoretical progress in computing accurately SM contributions, from the high-energy side (electroweak and perturbative QCD contributions) but also from the lowenergy side (hadronic matrix elements involving QCD in the non-perturbative regime). Global fits to the wealth of experimental data on flavour transitions show an overall excellent agreement with the SM picture, leading also to a precise determination of the CKM parameters [1].
Such a precise determination of the CKM parameters is essential for precise predictions of many flavour observables, used to set bounds on beyond-the-SM (BSM) physics. Some of these bounds are among the most stringent BSM constraints available, and can be translated to lower bounds on the BSM scale much above the reach of present and near-future colliders. However, these bounds must be extracted with great care, because the presence of new physics may well invalidate the assumptions implicit in the extraction of the CKM parameters in a SM-based analysis. The goal of this article is to address, in a systematic and model-independent way, how new physics (NP) will affect the global CKM fit and what is the best way to fix the full CKM matrix in a generic BSM context.
Direct searches for new heavy particles at the LHC have not been conclusive, which suggests a significant gap between the electroweak (EW) and BSM scales. In this context, the Standard Model Effective Field Theory (SMEFT) [2,3] represents an appropriate theoretical framework to analyse flavour data. Such data analysis is done once and for all in the SMEFT and the subsequent correlated bounds on the relevant Wilson coefficients can be applied to a plethora of NP models in a simple way. This is true whether one is addressing anomalies in the data, or just setting bounds on BSM physics. Moreover, the SMEFT allows one to account for non-trivial correlations between different classes of observables, such as quark-flavour transitions, leptonic processes, and EW precision measurements. In addition it embeds resummations that are needed to tame large logarithms in the presence of large scale hierarchies through Renormalisation Group Evolutions (RGEs).
The SMEFT is the effective theory of any fundamental theory that contains the SM supplemented by a set of heavy particles with masses M ∼ Λ m Z , and in which the EW symmetry is linearly realised. Technically speaking, it extends the SM Lagrangian with higher-dimensional operators built from SM fields, with the assumption that these operators also obey the gauge symmetry of the SM. The leading NP effects are typically encoded in the Wilson coefficients of operators of dimension six, and thus the primary goal is to establish confidence intervals for these parameters. However, it is important to realise that the free parameters of the SMEFT are not only the Wilson coefficients, but also the parameters already present in the SM Lagrangian (or the "SM parameters" with a slight abuse of language): the gauge and Yukawa couplings, the Higgs mass and its vacuum expectation value (VEV). In a consistent analysis one needs to take into account the presence of the NP contributions affecting the input observables from which the SM parameters are extracted.
The issue of NP "contamination" has been carefully studied in the context of EW precision observables (see e.g. Ref. [4]), where the relevant SM parameters are the SU (2) L × U (1) Y gauge couplings g L , g Y , and the Higgs VEV v. In the SM their numerical values are determined from three very precisely measured experimental inputs (see e.g. Ref. [5]): the fine-structure constant α em (0) (from the Rydberg constant), the Z-boson mass (from the Z lineshape in LEP-1), and the Fermi constant G F (from the muon lifetime). In the SMEFT, the relation between these observables and g L , g Y , and v is affected by dimension-six operators, and a straightforward extraction of the SM parameters is thus not possible. For example, the Fermi constant is given by G F = ( √ 2v 2 ) −1 (1 + δG F /G F ) where δG F is a linear combination of several dimension-six Wilson coefficients (see Section 3.2 for details). One option would be to perform a global fit simultaneously to the Wilson coefficients and the SM parameters. A more convenient and practical approach, however, consists in absorbing δG F into a redefinition of the SM parameters, in a procedure akin to the renormalisation of the SM at one loop. Namely, one can define the "tilde VEV" parameter v = v (1 + δG F /G F ) −1/2 such that it relates to the Fermi constant in the same way as the Higgs VEV in the SM: G F = ( √ 2ṽ 2 ) −1 , and has a well-defined numerical value,ṽ = 246.21965(6) GeV. Although δG F has been redefined away here, it does not disappear. Instead, NP corrections proportional to δG F emerge in SMEFT predictions of other electroweak observables (such as the W boson mass) that depend on the Higgs VEV in the SM limit, once v is traded forṽ. This kind of approach was followed in some previous global analyses of EW precision observables within the SMEFT, see e.g. Ref. [6][7][8][9][10].
We want to develop an analogous approach to deal with the CKM parameters in the SMEFT consistently. This task is much more challenging than in the EW sector where a small set of precisely measured and theoretically clean input observables can be distinguished. On the contrary, SM CKM fits rely on many distinct observables, often measured in elaborate experimental set-ups and displaying a complicated dependence on non-perturbative hadronic inputs. These require more involved theoretical approaches, sometimes relying on assumptions only justified within the SM. For this reason the current global CKM fits developed within the SM cannot be used without a careful adaptation in a general BSM framework such as the SMEFT. Perhaps for these reasons there is, in fact, no complete SMEFT analysis of flavour data available in the literature to this day. 1 Model-independent analyses available in the literature (and also many model-dependent ones) involve different classes of assumptions. A usual approach, motivated by setups with a very high NP scale and an arbitrary flavour structure, is to use ∆F = 2 processes to extract NP bounds [14][15][16][17][18][19], with the implicit or explicit assumption that the extraction of the CKM parameters themselves (which are used to calculate the SM prediction of these ∆F = 2 processes) is not affected by NP [16]. Moreover, it is sometimes assumed that the "PDG values" of the CKM parameters (obtained from the global SM CKM fits) can still be used in such BSM setups. Such assumptions are highly non-trivial and they greatly reduce the model-independent nature of these studies. The goal of this work is to provide the necessary results so that such assumptions are not needed.
In this article we propose a consistent framework to use the CKM parameters in the SMEFT. To this end, we will select a number of input observables, and identify specific combinations of the CKM parameters and Wilson coefficients that are determined by these observables. These combinations will define the "tilde Wolfenstein parameters" which, in analogy to the tilde VEV, can be used to predict numerical values and NP dependence of other flavour observables. The outline of the article is the following. Section 2 introduces the theory framework and notation. In Section 3 we describe our strategy to extract the CKM parameters in the general context of the SMEFT, and justify our choice of the input observables. In Section 4 we extract numerical values for the tilde Wolfenstein parameters from the input observables, and give the necessary formulas to apply our results in phenomenological SMEFT applications. In Section 5 with discuss some examples of applications of our formalism. Section 6 contains our conclusions and some future perspectives.
2 Theoretical framework 2.1 Fermion masses and CKM matrix beyond the SM We assume that there is a hierarchy between the EW and NP scales (µ EW Λ NP ), and that EW symmetry breaking is linearly realised. In that case the physics at the EW scale is described by the SMEFT [2,3]: and C i are respectively the dimension-six effective operators and their Wilson coefficients, and the dots include operators that violate lepton or baryon number and operators of dimension larger than six, which we will not consider. We will use the Warsaw basis and the notation and conventions in Ref. [20] except for the fact that we include 1/Λ 2 in the coefficients C i .
In the broken phase, the Lagrangian at the EW scale contains the fermion mass terms: where the mass matrices include contributions from the SM Yukawa couplings as well as from dimension-six operators: where v denotes the VEV of the Higgs doublet in the presence of dimension-six operators. It is always possible to define the Lagrangian in a "weak" basis for the fermion fields where the mass matrices are given by with V a unitary matrix. We will adopt this convention, in line with Refs. [21,22]. Thus, the right-and left-handed lepton and up-quark fields, as well as the right-handed down-quark fields are the same in the weak and mass eigenstate bases: e L,1 = e L , e R,2 = µ R , u L,2 = c L , d R,3 = b R , etc., while the translation from weak to mass eigenstate flavour indices for the left-handed down quarks is given by the V matrix: Formally, V ix has a weak index i = {1, 2, 3} and a mass-eigenstate index x = {d, s, b} [22]. Since in our convention the weak and mass bases for up-type quarks are the same, it holds that V 1x = V ux , V 2x = V cx and V 3x = V tx , and from now on we can use V rx with both mass-eigenstates indices r = u, c, t and x = d, s, b. In this article we use the Wolfenstein parameterization for V : We refer to the unitary matrix V as the CKM matrix. Its definition is affected by the presence of certain dimension-six operators, cf. Eq. (2.3). Moreover, and contrary to the SM, the flavour structure of charged currents is not uniquely determined by the CKM matrix, but is also affected by the presence of dimension-six operators with generic flavour structure. 2 In the following we discuss the consistent extraction of V from flavour observables within the general context of the SMEFT.

Effective theory below the EW scale
While it is possible that, in the future, precision high-energy measurements at the EW scale might be used to extract the parameters of the CKM matrix (see e.g. [23]), low-energy flavour-violating observables remain currently the best window to CKM physics. These observables are calculated in an effective theory where particles with EW-scale masses have been integrated out [22,[24][25][26]. Low-energy flavour observables probe directly the Wilson coefficients of the operators in this Lowenergy EFT (LEFT) at the appropriate hadronic scale, which can be related to the SMEFT through RGE together with a matching at the EW scale. In this article we will use the LEFT basis and notation of Ref. [22]:  quarks and leptons except the top quark [26], while for physics at lower energies one may integrate out additional fields such as the b quark. Anticipating the relevant observables that will be chosen in Section 3.3 to fix the CKM matrix, we focus on the semileptonic and ∆F = 2 operators (mediating d → uµ −ν µ , s → uµ −ν µ , b → uτ −ν τ transitions as well as B d and B s mixings). The relevant LEFT operators are collected for convenience in Table 1.
The tree-level matching conditions for the full set of Wilson coefficients L i in terms of the SMEFT Wilson coefficients can be found in Ref. [22]. The matching conditions in the SM are known to much higher orders (see e.g. [24,27]), while some one-loop contributions from non-SM operators are also known [28][29][30]. We will consider the state-of-the-art SM matching conditions but only tree-level matching from dimension-six operators in the SMEFT, as given in Ref. [22]. A typical matching condition has the structure: where F SM i ( g, m, µ EW ) are the SM matching conditions as functions of the set of couplings and masses in the SMEFT (collectively called g and m), and the product F (6) ij ( g, m, µ EW ) C j (µ EW ) denotes the contribution from the dimension-six SMEFT operator Q (6) j . The relevant expressions for the specific SMEFT operators needed in the analysis of Section 4 are given in Appendix A.
The low-energy amplitudes used to compute the processes of interest for the CKM parameters are given by default in terms of LEFT Wilson coefficients at a low, hadronic scale, where non-perturbative matrix elements are computed: L i (µ b ∼ 4.3 GeV) in the case of B physics and L i (µ s ∼ 2 GeV) in the case of K physics. In order to relate the coefficients at these scales with the matching conditions at the EW scale without generating large logarithms one needs to use RGEs. As in the case of the matching conditions, the anomalous dimensions of the SM operators are known to high orders in α s . For some sets of BSM operators, two-or three-loop anomalous dimensions in QCD are also known, including the operators that will be relevant in Section 4 [31,32]. One-loop anomalous dimensions in QED+QCD are known for the full set of LEFT operators [26,33], which are implemented in publicly available software tools [10,34]. The RG evolution, which can be implemented matricially as is included in the relevant formulas given in Appendix B.
3 Model-independent determination of the CKM matrix

Basics of the CKM fit in the SM
In the SM, the numerical values of the Wolfenstein parameters W i ≡ {λ, A,ρ,η} are extracted from a global fit to a long list of experimental observables (see e.g. [1]) that are accurately measured and whose SM predictions are well understood. They can be separated in four broad categories: • Leptonic decays (∆F = 1 branching ratios): the branching ratios provide information on the modulus of a CKM matrix element, provided that one knows the corresponding decay constant, i.e. the coupling between the axial current and the relevant meson. Currently, the main observables with accurate theoretical and experimental inputs are π → µν , K → eν , K → µν , τ → Kν , τ → πν , • Semileptonic decays (∆F = 1 branching ratios): these measurements provide information on the modulus of a CKM matrix element, provided that one knows the corresponding form factors, i.e. the couplings between the vector/axial and scalar/pseudoscalar currents and the relevant mesons. Currently, observables with accurate theoretical and experimental inputs are: Also in this category we can include the superallowed nuclear transitions, which provide a precise value for |V ud | and rely on a detailed description of nuclei and their weak transitions. There are also determinations of |V ub | and |V cb | from inclusive b → c ν and b → u ν transitions, which are however not fully compatible with the determinations from exclusive transitions (hinting at underestimated systematics). We can also include in this category |V us | determinations from inclusive τ →ūsν decays.
• CP-asymmetries (∆F = 1 CP-violating observables): These measurements allow one to extract CP-violating phases (α, β and γ, see Ref. [1] for their definition). They typically combine information from different channels or exploit time-dependent asymmetries involving the same hadronic matrix elements, so that the latter can be determined from the data or cancel out in ratios depending only on the CKM elements. The presence of the same hadronic matrix elements may hinge on SM flavour symmetries (isospin symmetry for the angle α), the hierarchy of CKM contributions (neglect of penguins for the angles β and β s ) or the knowledge of hadronic matrix elements from other sources (hadronic D decays for the angle γ). The main current processes of interest are: • Neutral-meson mixing (∆F = 2 observables): these measurements deal with the time evolution of the system composed by a neutral-meson flavour state and its CP-conjugate. They measure properties of the transition from one mass eigenstate of a neutral-meson system to the other (difference of masses, CP violation). They rely on the knowledge of the matrix element of the relevant ∆F = 2 operator in the SM between both neutral mesons. The main current observables of interest are All these measurements show a remarkable agreement with the CKM picture for the quark-flavour transitions embedded in the SM, leading to an accurate determination of the four CKM parameters [1]. However, these results have to be reassessed in the presence of NP, as non-SM contributions may not respect the assumptions implicit in their derivation. In the remainder of this article we propose an algorithm for extracting the CKM parameters in a general SMEFT framework, and discuss how to translate the measurements of flavour observables into constraints on NP in a consistent way.

Interlude: Higgs VEV in the SMEFT
Before embarking on the extraction of the CKM parameters in the SMEFT, it is worth recalling how parameters in the EW sector can be defined to illustrate this strategy. We take as an example the Higgs VEV in the SMEFT. In the SM, v is related to the Fermi constant G F , which in turn can be defined as a coefficient of the 4-fermion interaction between muons, electrons and neutrinos in the effective theory at a scale µ ∼ m µ : Integrating out the tree-level W exchange in the SM one finds G F = ( √ 2v 2 ) −1 . Given the experimental value G F = 1.1663787(6) × 10 −5 GeV −2 [5] precisely measured in muon decay, one can assign the numerical value to the Higgs VEV, v = 246.21965(6) GeV. However, this logic is perturbed if the SMEFT (and not the SM) is the relevant theory at µ m W . In that case one finds that dimension-six operators affect the Fermi constant as 3 where v is the VEV of the Higgs field in the presence of dimension-six operators, and C , C are Wilson coefficients of the corresponding operators in the Warsaw basis [20]. At this stage we cannot assign a numerical value to v without knowing the Wilson coefficients. Instead, it is convenient to define the tilde VEV parameterṽ via the relatioñ With this definition we recover G F = ( √ 2ṽ 2 ) −1 , and we can assign a numerical value toṽ, which is equal to that of v in the SM context,ṽ = 246.21965(6) GeV. In fact, this procedure is similar to the renormalisation of the SM at one loop. Let us however stress that we are dealing with finite tree-level corrections in the present situation.
At this point the dependence of the muon decay width on the SMEFT Wilson coefficients has been absorbed into the definition ofṽ, hence this observable alone does not constrain NP. However, the physical effect of δG F is not void. Using Eq. (3.7), we should replace v withṽ in the expression for any other EW observable sensitive to the Higgs VEV in the SM limit, in order to isolate the SM prediction for that observable. This way, δG F will modify the linear combination of Wilson coefficients to which the observable is sensitive: The "direct" contribution comes from the computation using the initial SMEFT parameters, whereas the "indirect" part comes from the redefinition of the Higgs VEV. We remark that the separation between direct and indirect NP effects is semantic. Both effects are in general equally large and physical, namely O(Λ −2 ) in the SMEFT expansion.

Strategy for the extraction of the CKM parameters in the SMEFT
We turn to the determination of CKM parameters in the general context of the SMEFT. There are two distinct strategies one could envisage here. One could aim at performing a global fit to all available flavour observables where not only the dimension-six Wilson coefficients but also the 4 independent CKM parameters are considered free parameters. Treating the unknown Wilson coefficients as nuisance parameters would return confidence intervals for the Wolfenstein parameters. This is a formidable task, and a much greater challenge than the SM CKM fit considering the number of parameters involved. An illustration of this strategy in the more limited case of NP only in the ∆F = 2 sector can be found in Refs. [17][18][19].
In this article we opt for a simpler strategy. We will identify a minimal set of four optimal observables constraining specific combinations of CKM and SMEFT parameters. These observables will define the 4 Wolfenstein parameters to which we will assign numerical values and errors. These in turn can be used to predict numerical values of other flavour observables, which can be compared with the experimental values in order to constrain NP. Note that our strategy can be embedded in the former one at a later stage, by using the results obtained in this article as "pseudo-observables" in a global fit.
In line with the discussion in Section 3.2, we will denote the combinations of Wolfenstein parameters and SMEFT Wilson coefficients extracted from the selected observables by (denoting ρ ≡ρ andη ≡η) with the understanding that in the SM limit W j → {λ, A,ρ,η}. We will often refer to these quantities as tilde Wolfenstein parameters, or simply tilde parameters.
In order to determine W j one should choose a quartet of observables from the pool of observables listed in Section 3.1. Ideally, we want the input observables to satisfy the following rather constraining conditions: 1. The set of observables must have a good sensitivity to all the four Wolfenstein parameters; 2. Each observable should be accurately measured and its theoretical prediction (in the general BSM case) should be well understood and non-controversial; 3. The general SMEFT expression for the observables should involve as few SMEFT Wilson coefficients as possible, in order to minimize the number of correlated observables needed in phenomenological applications.
This greatly reduces the available choices from the list in Section 3.1. Condition #1 is obviously mandatory and can be checked in the SM limit. Condition #2 rules out b → c ν ( = e, µ) transitions due to the tensions observed between exclusive vs. inclusive determinations [1]. Condition #3 is needed to select observables that can be used in the general SMEFT framework while depending only on a limited set of theoretical inputs and unknown parameters.
Let us now discuss some classes of the observables from the CKM fit in the SM in Section 3.1 in more detail. Observables from non-leptonic decays in Eq. (3.3) involve a limited set of hadronic matrix elements in the SM, which can be determined or eliminated thanks to additional observables and symmetries. Beyond the SM, however, these observables involve a much wider set of hadronic matrix elements that are currently not known and, in a general SMEFT context, cannot be related to other hadronic quantities through flavour symmetries. A similar issue affects K , which can be extracted from K → ππ decays only under specific assumptions about the weak amplitudes.
Concerning the semileptonic decays such as K → π ν, D → K ν, or B → π ν, the rates depend on form factors whose momentum dependence is usually extracted from the measurement of the differential distributions, which are themselves modified by BSM effects. Thus in order to use this information, a new BSM analysis of both differential distribution and rate is required (see e.g. Ref. [11]). This is in contrast to the leptonic decays, whose hadronic input is limited to decay constants, well known from lattice QCD. In addition, semileptonic decays are often sensitive to a larger set of BSM operators than leptonic decays, disfavouring semileptonic decays on the basis of condition #3. Overall these arguments favor using leptonic as opposed to semileptonic decays as our input observables.
We can now determine the most appropriate observables for the determination of the CKM parameters. Concerning observables sensitive (only) to λ, condition #2 suggests to disfavour D and D s meson decays compared to K decays. The latter are measured with a better accuracy and thus exhibit better sensitivity to λ. One technical complication, however, arises due to the dependence of the leptonic K decays on the decay constant f K + , as its most recent determinations rely on the "experimental" value of f π from π → µν to set the reference scale in the lattice QCD calculations [35]. This reintroduces an SM assumption (i.e., that the pion leptonic decay is completely dominated by SM contributions) that is not appropriate for a general analysis in the SMEFT setup [11]. To avoid this complication, we take the ratio Γ(K → µν) to Γ(π → µν) as our input observable, as the lattice determinations of f K + /f π + are free from this problem (and known with higher accuracy). Concerning the parameter A, we may consider observables sensitive to V ub , V cb , V td , or V ts , while the highest sensitivity toρ andη comes from V ub and V td . All in all, the remaining observables satisfying our criteria and sensitive to these three CKM parameters are B → τ ν (for V ub ), ∆M d (for V td ), and ∆M s (for V ts ).
This leaves us with the following set of input observables that we consider optimal: These four observables indeed obey the criteria listed above. In Section 4 we will show that they provide an accurate determination of the four Wolfenstein parameters W j in the generic SMEFT case, with only a moderate loss of accuracy compared the SM case. One should stress that our choice is not set in stone, and some variations on the input observables are of course possible, similarly to different input schemes used in EW precision physics. Furthermore, we emphasise that the "optimal choice" may vary over time. For example, if the inclusive-vs-exclusive tensions for b → c or b → u transitions disappear, or (theoretical or experimental) progress is achieved in some of the flavour transitions that we dismissed, our input observables may need to be appropriately reconsidered.
Summarizing, our approach to constraining NP in the SMEFT using flavour observables adheres to the following algorithm: 1. We identify the dependence of the input observables O input i in Eq. (3.10) on the LEFT Wilson coefficients L k and, given the LEFT-to-SMEFT map, on the Wilson coefficients of dimensionsix operators in the SMEFT C k : where we keep the dependence of the f and g functions on the CKM parameters W j ≡ {λ, A,ρ,η} as implicit.
2. We define the parameters W j ≡ {λ,Ã,ρ,η} by where δW j /W j are functions of the LEFT (or SMEFT) Wilson coefficients. They are defined such that the input observables depend explicitly only on W j in a way similar to the SM expression and involving no additional LEFT or SMEFT Wilson coefficients: 3. We extract the numerical value of the tilde parameters W j using Eq. (3.13) along with the necessary experimental and theoretical inputs (as described above), keeping full track of correlations. These values, their correlated uncertainties, and the contribution from BSM operators, given in Eq. (3.12), will be the main result of this work.
4. At this point we can translate any other flavour measurement, O α , into a model-independent NP constraint: where δO direct α,NP stands for a combination of Wilson coefficients contributing directly to the observable, and the indirect contribution is 4 Eq. (3.14) is the flavour analogue of Eq. (3.8) relevant for EW precision observables.
Once the tilde Wolfenstein parameters have been determined, it is convenient to introduce the tilde CKM matrix V . Given the SM expression V (λ, A,ρ,η) in Eq. (2.6), we define V by The elements of this matrix can be used to calculate the numerical SM predictions for observables depending on the CKM parameters. The NP effects included in them should be taken into account through the method described above. The V matrix defined above is unitary by construction. This does not entail any loss of generality, because we do not define the nine elements of V as the elements extracted from nine different measurements (such matrix would not be unitary in the SMEFT). Unitarity is a key and necessary ingredient, since we only have four independent CKM parameters to fix, and thus we only need to "lose" four measurements (and not nine) to fix them. Any additional observable becomes in this way a NP probe, as it should be. A last comment is in order concerning the choice of the hadronic inputs related to these observables. Lattice QCD provides a self-consistent theoretical framework to compute these inputs in global CKM analysis, but it still requires phenomenological inputs to determine the values of the parameters of the Lagrangian: the quark masses and the strong coupling constant (i.e. the lattice scale in physical units). However, as discussed above for f K , the "experimental" value of f π from π → µν is often used to set the scale in the lattice QCD calculations. This reintroduces an SM assumption (i.e., the pion leptonic decay is completely dominated by SM contributions), which will propagate in all dimensionful lattice QCD inputs and which is thus not appropriate for a general analysis in the SMEFT setup. From this point of view it is thus better to use determinations of the scale where an observable dominated by strong dynamics is used to set the scale (for instance the masses of hadrons, or the quark-antiquark potential).

Analysis and Results
4.1 K → µν µ , π → µν µ and B → τν τ We start with the leptonic decay rates Γ(P + → + ν ), with P = {π, K, B}, and = {µ, τ }. For an exhaustive study of V us from K → µν µ along the lines of the present article, we refer to Ref. [11]. The decay rate for the process P − → −ν can be written as where f P is the decay constant defined by 0|qγ µ γ 5 u|P + (k) = ik µ f P , and the quantity δ P accounts for all electromagnetic corrections in the SM (see e.g. [36]), 5 as well as isospin-breaking corrections if they are not already included in the decay constants [37]. The full and linearized expressions for the NP corrections are given by Factoring out the SM corrections induces (tiny) NP × QED corrections, which may not be the correct ones but which are beyond the current theoretical precision.
The definitions ofṽ and δv/v can be found in Eq. (3.7). The parameters uq A are O(Λ −2 ) in the SMEFT expansion, and they are connected to the LEFT Wilson coefficients at the hadronic scale by: Using the input values collected in Table 2, from the input observable Γ(B → τ ν) we obtain where we neglected the electromagnetic correction (i.e., δ B 0) since it induces an effect much smaller that the current experimental sensitivity (see Ref. [38] for further detail on this issue). The error in | V ub | (12%) is dominated by the experimental uncertainty on B(B → τ ν τ ). The value used for the decay constant f B ± [39] and quoted in Table 2 does not rely on the experimental value for f π , c.f. the discussion in Section 3.3.
In a similar fashion we could determine | V us | from the observable Γ(K → µν), which would translate directly into a determination of the tilde Wolfenstein parameterλ. However, as discussed in Sec. 3.3, this choice would lead to a relatively large uncertainty onλ, which come in particular come from the lattice input for f K : the MILC09 calculation [40] is the most precise lattice determination of f K not relying on the pion leptonic width to set the lattice scale, with fairly large uncertainties. A more precise value ofλ can be obtained by considering instead the ratio Γ(K → µν µ )/Γ(π → µν µ ), which allows one to extract the ratio | V us / V ud | given the ratio of decay constants f K /f π . The latter can be consistently taken from the FLAG average [35], which combines several lattice determinations for this ratio of decay constants without introducing any uncontrollable dependence on NP via the pion leptonic width [35,[41][42][43].
In this case we have: with and Given the inputs in Table 2, we find with a relative error of 0.2%, dominated by the uncertainty on the lattice determination on the decay constant ratio.

∆M d and ∆M s
The mass differences ∆M q of neutral B q mesons (q = {d, s}) are given by [44,45] where B q i are the so-called bag parameters, defined e.g. in Ref. [46] 6 . The quantities V tb and V tq have been defined such that: with the BSM effects contained in where we have defined C  Table 2 of Ref. [46]. It is well known that the combination The bag parameters B q i are the matrix elements B q |O i (µ)|B q up to a normalization factor. With our conventions, B q i denote the bag parameters in the renormalisation scheme of Ref. [31] at a scale µ b = 4.3 GeV, as given Table 2 of Ref. [46], and in agreement with the RGE factors used in Appendix B. It is worth mentioning that f π has been used to set the scale in Ref. [46]; however, this choice induces a subleading effect on the determination of the (dimensionless) bag parameters, such that we can safely use this computation given the current uncertainties. Γ(K + → µ + ν µ ) = 3.3793(79) · 10 −8 eV [5] δ K/π = −0.0069(17) [47] Γ(π + → µ + ν µ ) = 2.5281(5) · 10 −8 eV [5] f K ± /f π ± = 1.1932 (19)   is more precisely determined in the lattice than numerator and denominator separately due to the presence of parametric correlations. Often this is exploited by trading one mass difference by the observable ∆M s /∆M d . Instead, we will take into account these parametric correlations by writing 2) only the tree-level matching conditions between the LEFT and the SMEFT, but higher order corrections can be included trivially, once they are known (e.g. [29]). This might be relevant since non-log-enhanced GIM-violating contributions may be numerically large in observables such as ∆M q . This is an important issue to keep in mind.
Using the numerical inputs from Table 2, we obtain | V tb V td | = 0.00851 ± 0.00025 , and | V tb V ts | = 0.0414 ± 0.0010 . The errors in | V tb V td | (2.9%) and | V tb V ts | (2.5%) are both dominated by the uncertainties of the f 2 Bq B q 1 combinations. Again, the f Bs value [39] quoted in Table 2 does not rely on the experimental value for f π , according to the discussion in Section 3.3.

Summary and Results
To summarise, we have obtained the following numerical constraints on the tilde CKM elements: with negligible correlations, except for the last two elements that present a correlation of +87%. Moreover, we have identified the new physics contributions to the above quantities, which can be found in Eqs. (4.2), (4.7), and (4.11).
It is now possible to write our results in terms of the tilde Wolfenstein parameters W i , making use of the following relations: Given these definitions, one can translate Eq. (4.14) into correlated constraints on the CKM parameters W i . We find the following results:  where ∆ K/π , ∆ Bτ 2 , ∆ ∆M d , and ∆ ∆Ms are the (linearized) NP contributions to the four chosen observables, which can be found in Eqs.
The numerical value of M is given by  Table 3 summarises our results for the Wolfenstein parameters in the presence of NP, and compares them to the results of the canonical SM fits. As could be expected, our procedure of using only four input observables to determine the four Wolfenstein parameter leads to some loss of accuracy in the limit where BSM corrections are absent, compared to the SM fits using a much larger set of observables. 7 Nevertheless, in most physical applications the error bars of our tilde parameters will be anyway a subleading effect compared to other sources of experimental or theoretical uncertainties (as illustrated in a few concrete examples in the next section). On the other hand our tilde parameters can be consistently used for generic NP frameworks described by the SMEFT, unlike the results of the SM fits. We remark that our input observables allow other solutions than the one displayed in Eq. (4.16); in particular, there is another solution with the opposite sign ofη. This discrete ambiguity will lead to "mirror solutions" in global fits, where the CKM parameters differ significantly from the ones obtained in the SM context, but the resulting shift of precisely measured flavour observables is canceled by a relatively large (and 7 The difference of precision between our determination ofλ and the SM CKMfitter determination of λ is partially due to the use of a larger set of constraints in the latter case, but mainly due to current internal tensions between some of the constraints on λ (within the SM), which generate a smaller error on λ in the CKMfitter statistical approach. Compared to UTfit, we obtain a more preciseλ because we use the new FLAG average for f K /f π [35]. fine-tuned) contribution from SMEFT Wilson coefficients. In this article we will not discuss the mirror solutions any further, and focus on the SM-like solution in Eq. (4.16) where the NP effects are subleading compared to the SM contributions.
With the likelihood function in Eq. (4.16) we find the following 1σ intervals for the elements of the tilde CKM matrix defined in Eq. We do not give here the (non-trivial) correlations between the various tilde CKM elements, but they are encoded in the likelihood in Eq. (4.16). The NP effects absorbed in these CKM elements should be taken into account through the method described in Section 3.3; see Section 5.2 for an example. The numerical form of the matrix V is given here for illustration purposes only: ideally one should always express all CKM input in terms of Wolfenstein parameters if one wishes to use the approach and results of this paper appropriately, including correlations.

Applications
In this section we discuss through a few examples how to use the tilde parameters to analyze, in a consistent fashion, other flavour processes and set bounds on NP.

Leptonic decays of pions and D mesons
Consider the pion decay π → µν. The goal is to compare the precisely measured branching fraction to the SM predictions so as to place constraints on effective interactions beyond the SM. The SM prediction is proportional to |V ud |. In the presence of generic NP we cannot use |V ud | determined by fits performed in the SM context, such as the ones provided by CKMfitter or the PDG, as the observables used in those analyses may themselves be affected by new physics, which might even have the same underlying effective operators. Instead, we can use our results in Eq. (4.16) where the NP effects have been absorbed into the definition of the tilde parameters. All we need to do is to express the theoretical prediction for B(π → µν) in terms ofλ defined in Eq. (4.17). The π → µν decay width can thus be written as where the decay constant is f π ± = 130.2(8) MeV (from the average of Ref. [35] with a lattice scale set using QCD observables [50][51][52]), the electromagnetic radiative corrections are encoded in δ πµ = 0.0176(21) [36], and ∆ πµ2 is given by The NP quantities µud A,P , δv/v and δλ have been defined in Eqs. (4.3), (3.7), (4.17). The terms proportional to δλ are due to NP affecting the observables used to determine the CKM parameters in our approach; note that they depend on the same Wilson coefficients that also enter into µud A,P . Their effect is to change the linear combination of Wilson coefficients∆ πµ2 probed by π → µν decays. Given the current experimental measurement, B(π → µν) = 0.9998770(4) [5], combined with τ π = 2.6033(5) · 10 −8 s, we obtain the following constraints on the linear combinations of Wilson coefficients in Eq. (5.2): ∆ πµ2 = 0.004 ± 0.013.
The error is totally dominated by the lattice uncertainty on f π ± . The error of our determination ofλ in Eq. (4.16) is completely negligible for this constraint. Up to small O(λ 4 ) corrections, the CKM elements V ud,us,cd,cs are only functions of λ in the Wolfenstein parameterization. Thus, besides pion decays, there is a long list of observables which are only functions of λ and NP parameters. A global fit to d → u ν and d → u ν transitions was performed Ref. [11], where simultaneous constraints were derived on λ and the relevant LEFT Wilson coefficients. We note that such a global approach obviates the need to define the tilde CKM parameters; 8 however, so far it was realised only for observables depending on λ, and extending it to the full set of the CKM parameters will involve considerable technical difficulties. The approach in this article bypasses this problem when setting constraints from individual observables, as has been exemplified here.
Going beyond the analysis of Ref. [11], we consider the decay D → ν. Analogously to π → ν, we write where f D ± = 212.7 (6) MeV [43], and and ∆ D 2 is given by where the NP quantities cd A,P can be found again in Eq. (4.3). The experimental measurement B(D → µν) = 3.74 (17) · 10 −4 [5] combined with the lifetime τ D ± = 1.040(7) · 10 −12 s, and δ Dµ = 0.007 (6) [43], results in the following constraint: showing a small 2.1 σ tension with the SM. Our analysis affects the NP interpretation of this result, adding the δλ term in Eq. (5.5). In this case this correction is enhanced byλ (unlike for π → ν where it is suppressed). One comment about the lattice input for f D ± is in order in connection with the discussion in Section 3.3. The value of f D ± obtained in [43] is normalised to the PDG value f exp π , which in principle propagates the NP contribution in π → µν into the D → µν constraints. However, the comparison of f lattice π and f exp π limits these NP effects to the 1% level, cf. Eq. (5.3), which makes them subdominant compared to the 4% error in Eq. (5.6) mainly due to the experimental uncertainty on B(D → µν). To avoid this issue altogether, lattice collaborations should quote f D ± setting the scale using a QCD-dominated observable free of NP. These considerations may be relevant in the future, when the experimental error on B(D → µν) is reduced.

Exclusive hadronic W decays
We now consider the processes W → u j d k , where u j = (u, c) and d k = (d, s, b) denote particular quark flavours. We assume that flavour tagging allows one to separate the distinct exclusive hadronic W decays, and that it is possible to measure the partial widths with a good precision. In the SM, the measurement of Γ(W → u j d k ) can be interpreted as an alternative probe of the CKM element V jk . Beyond the SM, the coupling strength of the W boson to quarks may be affected by new physics. In the SMEFT, the leading order effects can be described by the vertex corrections δg W q L to the couplings between the left-handed quarks and W : The right-handed vertex correction δg W q R does not affect W decays at O(Λ −2 ), and will be neglected in this discussion.g L ,g Y are the SU (2) L × U (1) L gauge couplings extracted from the EW input observables α, m Z in the presence of dimension-6 operators, in analogy toṽ extracted from G F discussed in Section 3.2. The left-handed vertex correction is related to the parameters in the Warsaw basis as [9] δg W q L jk = [C Naively, from Eq. (5.7) one could conclude that each exclusive W decays probes simply δg W q L jk and thus constrains the particular combination of the Wilson coefficients given in Eq. (5.8). However, to constrain new physics, the experimentally measured partial width has to be compared with the corresponding SM prediction. The exclusive decay widths predicted in the SM depend on the numerical value for the specific CKM matrix element V jk extracted from experiment, which in turn may be affected by NP. This effect can be disentangled in our scheme by trading V jk in Eqs. (5.7) and (5.8) for the tilde CKM element V jk defined in Eq. (3.16), i.e. V jk →Ṽ jk − δV jk . We then have where Γ(W → u j d k ) SM is the SM prediction calculated using the numerical values of the tilde CKM elements in Eq. (4.21). In this approach, the combination probed by the exclusive decay The CKM shifts relevant for W decays, up to O(Λ −4 ) and at leading order inλ, are given by δV cb = 2Ãλ δλ +λ 2 δA + O(λ 6 ) . (5.10) The corrections to the Wolfenstein parameters δλ, δA, δη and δρ in terms of EFT Wilson coefficients are defined in Eq. (4.17). If needed, (more lengthy) O(λ 4 ) corrections to the expressions in Eq. (5.10) can be easily calculated, but we do not list them here for the sake of simplicity. The data on exclusive W decays are presently very limited. We are only aware of the Delphi measurement of Γ(W → cs) [53] with the relative precision of ∼ 40%. However, progress in flavour tagging should enable more precise measurements of Γ(W → u j d k ) in the near future. For example, Ref. [23] argues that the measurement of Γ(W → cb) with the relative precision of order ∼ 15% should be possible in ATLAS or CMS using the existing data sets. In the scheme proposed in the present article, an LHC measurement of Γ(W → cb) probes not only the vertex correction δg W q L cb , but also 4-quark operators affecting the B s meson mass difference ∆M s . Our formalism allows for a consistent interpretation of these measurements as model-independent constraints in the SMEFT.

A Z model for b → s anomalies
To close this section, we discuss an application of our formalism in the context of a toy model addressing b → s ( = e, µ) anomalies, and in particular the ratios R K and R K * violating leptonflavour universality [54][55][56][57]. Let us consider a simple BSM toy model featuring a massive Z boson coupled in an SU (3)×SU (2)×U (1) invariant way to left-handed b and s quarks and to left-handed muons: Here 2 = (ν µ , µ L ) is the second-generation lepton doublet, and q 2 = (V † 2x u L,x , s L ), q 3 = (V † 3x u L,x , b L ) are the second-and third-generation quark doublets in the down-type basis. Integrating out Z yields new contributions to four-fermion contact interactions in the effective theory below the scale m Z . In particular, we generate a new contribution to the effective interaction (b L γ ρ s L )(μ L γ ρ µ L ) that adds to the SM loop-level contribution and may help to explain the R K ( * ) anomalies, as pointed out by several independent analyses [58][59][60][61][62][63][64][65][66] We have used the results of Ref. [58] for the best fit to b → s flavour observables. The small correction ρ V,LL ed takes into account the running from the Z mass to the b mass.
This constrains one combination of the three toy-model parameters g bs , g µµ , and m Z . In addition, other operators generated by integrating out the Z boson lead to further constraints on these parameters. First, the low-energy theory contains new contributions to four-lepton interactions that can be probed by the trident muon production in neutrino scattering [67][68][69]. Namely where ρ V,LL νe accounts for the running from m Z to the EW scale. The numerical value in the RHS was taken from the global fit in Ref. [9].
Furthermore, we generate the contribution to the ∆F = 2 operator responsible for B s mixing C where ρ (s) 1 accounts for RG running from m Z to the EW scale. One may be tempted to derive a constraint on the combination g 2 bs /m 2 Z using the measured value of ∆M s [5] together with its SM prediction calculated using the CKM elements extracted within the SM [70], namely ∆ ∆Ms = ∆M exp There is however one conceptual difficulty. The SM prediction for ∆M s crucially depends on the numerical value of the CKM matrix element V ts . In the SM context this value is extracted from a global fit to multiple observables (including ∆M s ) and it may be shifted in the presence of NP contributions. Our approach allows us to solve this conundrum. Since the ∆M s measurement has been selected as one of our input observables, it serves as an input to fix the tilde Wolfenstein parameters and by itself it does not constrain new physics. Of course, this does not mean that there is no possibility to probe C (s) 1 , but we need to use other observables than ∆M s such as the B → D ( * ) ν ( = e, µ) decays. Following the analysis of Ref. [71] in the SM limit, B → D * ν yields the 68% CL constraint on the CKM parameter |V cb | = (3.90 ± 0.07) × 10 −2 , while from B → D ν one obtains |V cb | = (3.96 ± 0.09) × 10 −2 . These extractions happen to be valid also in our model, since it does not introduce NP contributions to b → c ν transitions. Following the scheme proposed in this article, we can then relate the extracted |V cb | to its tilde value, which gives us the following constraint on the model: where we used Eq. (4.17) to relate δA to ∆ ∆Ms . We used as well that ∆ K/π (and thus also δλ), ∆ Bτ 2 , and ∆ ∆M d are zero in this model. Illustrative limits on our model are shown in Fig. 1 in the case m Z = 100 GeV. Together with the trident constraints, B → D ( * ) ν leaves a corner of allowed parameter space in the g µµ -g bs plane where the Z couplings are small enough. In our approach, B → D * ν provides the strongest constraint on g bs , which is however weaker than what would be found if we (incorrectly) used ∆M s to set limits. The outcome of the two approaches could differ even more significantly in more general situations, for example when an extra gauge boson also generates bc ν effective operators at low energies.

Conclusions and outlook
In the present article we have discussed the role played by the CKM matrix in the search for NP in the model-independent context of the SMEFT. We recall that the determination of the CKM parameters themselves are then affected by the presence of dimension-6 operators, and that the results from SM global fits combining all available observables performed by CKMfitter [14] or UTfit [15] cannot be used directly to exploit other flavour constraints involving the CKM matrix.
We have identified a set of four observables: which are deemed appropriate to determine the CKM parameters in the context of the SMEFT, based on the accuracy of the measurements, the theoretical understanding of their computation, the precision reached on the hadronic inputs, and the simplicity of their calculation within the SMEFT. We have determined the NP corrections to these four processes in the low-energy EFT (below the weak scale), and expressed these corrections in terms of SMEFT contributions by running from the low hadronic scale to the weak scale and performing a matching at the latter. The corrections from dimension-six operators can then be included in the definition of the "tilde parameters" W j ≡ {λ,Ã,ρ,η} (in a procedure similar to one-loop renormalisation), and constraints on W j can be extracted. Our results for the tilde Wolfenstein parameters are: In this exercise, particular attention must be paid to hadronic inputs from lattice QCD simulations: we select data with a lattice scale set by pure QCD quantities (such as bound-state masses) and not by quantities involving the weak interaction (such as decay constants) and potentially modified by NP.
We have also discussed several applications (leptonic meson decays, hadronic W decays, and constraints on a Z model) to illustrate how our approach leads to a clear interpretation of the measurements of quark-flavour observables and a separation between NP contributions coming from the determination of the CKM parameters and those linked to the process of interest.
The result of our analysis is a set of observables that differs from more usual choices. Indeed, the generality of our approach prevents us from using a rather common option, i.e. using only tree-level processes to extract the CKM parameters. It is often advocated that SM loop-level transitions are much more sensitive than SM tree-level transitions to NP effects, so that tree-level transitions should be used preferentially to determine the CKM parameters. One cannot assume this premise in the general SMEFT set-up. In addition, a hierarchy of NP effects between SM treeand loop-level transitions is not supported by the current B-anomalies -where the (potential) BSM contribution relative to the SM is of the same order in b → cτ ν (SM tree) and b → sµµ (SM loop) transitions (see e.g. [57]). It is also not present in well-known theoretical frameworks such as Minimal Flavour Violation [73][74][75].
The present article outlines the procedure to determine the CKM parameters in SMEFT analyses where only a subset of all flavour observables is taken into account. We have proposed a choice of input observables that we consider optimal at this point. This might change in the future if, e.g., new theory developments appear and/or if experimental measurements improve the consistency or the accuracy of some of the other constraints. Moreover, this choice among observables is not needed in a fit that includes all the measurements that are most sensitive to the CKM parameters, taking into account all correlations. In this case, all observables contribute to the bounds on all the parameters of the fit, and the separation between the processes (mainly) used to extract the CKM parameters and processes (mainly) used to set NP bounds is only useful in order to illustrate their different sensitivities to each type of parameters. This approach was illustrated in studies of NP restricted to ∆F = 2 transitions performed by UTfit [15] and CKMfitter [14,[17][18][19]. But such analysis was possible because of the very simple structure of the NP scenario considered. In more general settings, and in particular in the full SMEFT case, it is at the moment not possible to proceed in the same way. Indeed, in global SMEFT analyses so far the CKM parameters are not treated as free variables, and NP effects affecting their extraction are not taken into account, as discussed in Refs. [76,77]. Our work allows one to overcome this limitation, providing an appropriate framework to consistently include such NP effects as well as the uncertainty on the CKM parameters.
We have discarded here some of the observables because in the general SMEFT they have complicated expressions involving unknown BSM hadronic matrix elements that cannot be easily computed or connected to other hadronic inputs through symmetries. There are however more specific cases of NP where the expressions of these observables are simple and could provide interesting alternatives to the subset chosen here. The simplest example consists in the case of NP contributions with an SM-like structure for the additional operators so that only (V − A) × (V − A) charged currents are generated (see, e.g., Ref. [78]). In practical terms, it would then be useful to include these simplified expressions to be used only if the necessary assumptions are satisfied by the underlying NP operators. One would thus be able to recover in a trivial algorithmic way the setup described above where tree-level processes are approximately free of NP, or to recover the most accurate results from the SM global fit when the SM case is considered.
In connection with the flavour anomalies currently observed, the constraints on NP coming for flavour physics should be assessed with a particular attention. These deviations constitute essential probes of the physics at play at energies beyond the current LHC frontier, and one should aim at exploiting current and forthcoming data in global analyses combining large sets of observables from different sectors of particle physics. A careful determination of the CKM parameters in the SMEFT will thus play an important role in this strategy, which should ultimately provide us with new insights in the structure of the physics at higher energies, beyond the Standard Model. A Matching the LEFT to the SMEFT The relevant SMEFT operators in the Warsaw basis [20] onto which the LEFT operators from Table 1 match are collected in Table 4. Here, ϕ is a Higgs field, (e) is a left-(right-)handed lepton field, q is a left-handed quark field and u(d) is an up-(down-)type right-handed quark field, with {i, j, k, l} family indices, D µ is the covariant derivative, and T A are the Gell-Mann matrices.
We consider first the matching conditions for the LEFT operators relevant for the semileptonic charged-current transitions collected in the left column in Table 1. The tree-level matching conditions to the SMEFT at the EW scale are given by [22,25,26]: where we assume all SMEFT parameters given at the EW scale by default. These matching conditions refer to the coefficients of the LEFT operators in the mass basis (where x = {d, s, b}, i = {e, µ, τ } and k = {u, c}), which are related to the ones in the weak basis by the CKM matrix, by virtue of Eq. (2.5). The LEFT operators relevant forB q − B q mixing are collected in the second column of Table 1. The tree-level matching to the SMEFT at order 1/Λ 2 is given by [22]: and analogously for d → s. Some 1/Λ 4 terms are known [22,28], but have been dropped here. In addition, one-loop matching corrections give non-zero contributions to L S1,RR dd and L S8,RR dd , which we have also ignored, for simplicity. Some of these can be found in Ref. [28].
In the computation of the mass differences we will use the the traditional "SUSY basis" for the ∆F = 2 operators [44], for which the matrix elements are explicitly known. These are denoted by O

B Renormalisation Group Evolution
We start with the leptonic decays P → ν discussed in Section 4.1. Using three-loop QCD running plus one-loop QED running [26,32], the parameters lxy X defined in Eq. and corresponds to what is traditionally denoted byη B S 0 (x t ) 0.83798 × 2.36853 = 1.985, where the factorη B encodes NLO QCD matching corrections and running effects simultaneously (see e.g. Ref. [70]). In this way we have made the separation between EW matching and RGE running explicit. Note that in our conventions the sign of the Wilson coefficients C i is opposite than that in Ref. [44].
Hq ij = (ϕ † i ← → D I µ ϕ)(q i σ I γ µ q j ) ∆F = 2 Q Hud ij = i( ϕ † D µ ϕ)(ū i γ µ d j ) + h.c. Q (1) qq ijkl = (q i γ µ q j )(q k γ µ q l ) qq ijkl = (q i γ µ σ I q j )(q k γ µ σ I q l ) equ ijkl = (¯ m i σ µν e j ) mn (q n k σ µν u l ) Q qd ijkl = (q i γ µ q j )(d k γ µ d l ) Q edq ijkl = (l i e j )(d k q l ) Q (8) qd ijkl = (q i γ µ T A q j )(d k γ µ T A d l )  To describeB q − B q mixing we use the Wilson coefficients in the SUSY basis (see Eq. (A.2)) in the MS scheme of Ref. [31], at the scale µ b = 4.3 GeV, in accordance with the matrix elements of the operators provided in [46]. We call these Wilson coefficients C (q) 1,...,5 , C (q) 1,2,3 . In order to relate these coefficients to the ones at the EW scale, we use the NLO evolution matrix given in Ref. [31] (in the same scheme but different basis -see Section 3 of Ref. [81] for detail), running α s at four loops [82]. In the end we find, for µ EW = M Z and q = {d, s}:

C NP shifts to the CKM parameters beyond linear order
The relation between the NP shifts in the Wolfenstein parameters, δW = {δλ, δA, δη, δρ}, and the NP contributions to the chosen observables, denoted by ∆ ≡ {∆ K/π , ∆ Bτ 2 , ∆ ∆M d , ∆ ∆Ms }, is in general a complicated non-linear equation. However, assuming that NP is a small perturbation one can expand that equation around the SM point to any given order, and obtain its unique solution. Let us note that, by construction, this approach discards possible additional solutions where the NP correction is comparable to or larger than the SM contribution. The generalization of Eq. (4.17) to include quadratic NP terms is the following where the NP corrections to the observables, ∆, should include quadratic corrections. We have introduced the quantities where M corresponds to the matrix in Eq.
All quantities above are implicitly evaluated at W = W .