Revealing new physics in B 0 s → D ∓ s K ± decays

The B 0 s → D ∓ s K ± systemoffersadetermination of the Unitarity Triangle angle γ . Intrigued by an LHCb analysis showing a surprisingly large result in tension with information on the Unitarity Triangle and other γ measurements, we make a transparent study of the measured observables, conﬁrming the LHCb picture. The corresponding γ puzzle at the 3 σ level would require CP-violating contributions of New Physics, which should also manifest themselves in the corresponding decay branching ratios. Indeed, we ﬁnd that the rates of the individual B 0 s → D ∓ s K ± channels show puzzling patterns, in accordance with similar decays, with tensions up to 4 . 8 σ , thereby making the situation much more exciting. We present a formalism to include New-Physics effects in a model-independent way and apply it to the data to constrain the corresponding parameters. Interestingly, new contributions of moderate size could accommodate the data. Utilising this formalism in the future high-precision B physics era may allow us to ﬁnally establish new sources of CP violation.


Introduction
Decays of B mesons provide a wide spectrum of probes for testing the quark-flavour sector of the Standard Model (SM). A particularly interesting aspect is given by the subtle difference between weak interactions of particles and their antiparticles, which is described by the phenomenon of CP violation. Here C stands for charge conjugation and P refers to parity, i.e. space inversion. In the SM, CP violation is described through a complex phase in the Cabibbo-Kobayashi-Maskawa (CKM) matrix. Since the SM cannot explain the matter-antimatter asymmetry of the Universe, failing by many order of magnitude, new sources of CP violation are suggested.
The decaysB 0 s → D + s K − and B 0 s → D + s K − (with their CP conjugates) offer a powerful probe for testing the SM description of CP violation [1][2][3]. Due to quantummechanical B 0 s -B 0 s oscillations, interference effects are generated between these decay channels, leading to the following time-dependent rate asymmetry: C cos(∆M s t) + S sin(∆M s t) cosh(y s t/τ Bs ) + A ∆Γ sinh(y s t/τ Bs ) ; (1) an analogous expression holds for the CP-conjugate D − s K + final state, where C, S and A ∆Γ are replaced by C, S and A ∆Γ , respectively.
TheB 0 s → D + s K − and B 0 s → D + s K − channels arise from b → cūs andb →ūcs quarklevel processes, respectively. In the SM, the observables of Eq. (1) and its CP conjugate allow a theoretically clean determination of the angle γ of the Unitarity Triangle (UT) of the CKM matrix [1][2][3]. Performing these measurements and assuming certain SM relations, LHCb [4] has reported the result γ = 128 +17 −22 • modulo 180 • . This finding is puzzling as global SM analyses of the UT give values around 70 • [5][6][7][8]. This regime is also consistent with a recent simultaneous analysis of various B-meson decays by the LHCb collaboration [9]. However, these modes have dynamics different from the B 0 s → D ∓ s K ± channels and show sensitivity on γ through very different interference effects.
Intrigued by the challenge of interpreting this result, we have a transparent look at the determination of γ and the associated parameters, confirming the LHCb picture and excluding the modulo 180 • ambiguity [10]. If confirmed through future more precise measurements, this result could only be explained through CP-violating contributions of New Physics (NP) at the B 0 s → D ∓ s K ± decay amplitude level, and should then also leave imprints on the corresponding branching ratios. Consequently, we move on to extract the individualB 0 s → D + s K − and B 0 s → D + s K − branching ratios and compare them with the SM picture, developing a method to minimise dependences on non-perturbative hadronic and CKM parameters. We in fact arrive at yet another puzzling situation, complementing the γ measurement. The situation is even more exciting, as we observe similar patterns also in branching ratios for other B decays with similar dynamics. Surprisingly small branching ratios in the latter modes have already received attention [11,12], also within NP interpretations [13][14][15].
In view of these findings, we present a generalized description of the B 0 s → D ∓ s K ± decays and corresponding analysis of CP violation to include NP contributions in a model-independent way. Applying this formalism to the current data, we calculate cor-relations between the NP parameters, involving in particular new sources of CP violation. This compact paper complements a much more detailed discussion given in Ref. [10].

Standard Model Framework
The interference effects between B 0 s orB 0 s mesons decaying into the final states D + s K − and D − s K + caused by the B 0 s -B 0 s oscillations are described by physical observables ξ and ξ, respectively. In the SM, their product provides a theoretically clean expression [2]: where φ s is the CP-violating B 0 s -B 0 s mixing phase. It should be noted in particular that non-perturbative hadronic matrix elements cancel in (2). Since ξ andξ can be determined from C, S, A ∆Γ and C, S, A ∆Γ , respectively [3], and φ s through measurements of CP violation in B 0 s → J/ψφ [16], this relation allows a clean extraction of γ. Furthermore, the following relations hold in the SM framework [2]: which were assumed by the LHCb collaboration [4]. Performing a sophisticated fit to their data, LHCb found |ξ| = 0.37 +0.10 −0.09 (4) with the following results [4]: were δ s describes the CP-conserving strong phase difference between theB 0 s → D + s K − and B 0 s → D + s K − amplitudes. Here we have used φ s = (−1.7±1.9) • , which was employed in the LHCb analysis, to convert γ into φ s + γ. Using the updated result φ s = −5 +1.6 which includes also penguin corrections in B 0 s → J/ψφ modes [16], we obtain In view of the tension of this value with the SM and the complex LHCb analysis, it is crucial to transparently understand the situation. How can we achieve that? Using (3) with the measured value of C, we find |ξ| = 0.40 ± 0.13, which is in excellent agreement with (4). Introducing the combinations we obtain the following relations [2,3,10]: where the numerical values correspond to the LHCb measurements of the corresponding observables, yielding φ s + γ = (125 +18 −22 ) • with δ s = (2 +34 −22 ) • . Here we have excluded the solutions modulo 180 • as they would be in huge conflict with factorization predicting δ s ∼ 0 • [10]. Consequently, we find excellent agreement between this simple -but transparent -analysis and the complex LHCb fit. The result for γ in (6) is much larger than the 70 • regime, and shows a discrepancy at the 3σ level. It is important to stress that this intriguing tension could not come from any long distance effects, as the γ determination is theoretically clean in SM. Should this puzzle remain once the experimental picture sharpens further, it would require new sources of CP violation. Since the experimental value of φ s used in the analysis includes possible CP-violating NP effects in B 0 s -B 0 s mixing, new contributions entering directly at the decay amplitudes of the B 0 s → D ∓ s K ± system would be required. Such NP effects should manifest themselves also in the corresponding branching ratios. Let us therefore focus on these observables next. Due to the time-dependent B 0 s -B 0 s oscillations, we have to distinguish between time-integrated "experimental" branching ratios [17] and their "theoretical" counterparts where such mixing effects are "switched off" [3,18]. Moreover, we have to disentangle the interference effects between the two decay paths that arise from the B 0 s -B 0 s oscillations. We obtain the following expressions for the theoretical branching ratios [10]: Analogous expressions hold for the branching ratios of the B 0 s andB 0 s decays into the final state D − s K + . Unfortunately, separate measurements of the experimental branching ratios for these final states are not available, just the average with B exp Σ = (2.27 ± 0.19) × 10 −4 [6]. Assuming SM expressions for the decay amplitudes, we have [2,3] with y s = 0.062 ± 0.004 [6] and A ∆Γ + = 0.35 ± 0.23 [10]. Finally, we obtain [10]: For the calculation of these branching ratios and the underlying decay amplitudes, factorization provides the theoretical framework. Here the corresponding hadronic matrix elements of four-quark operators entering the low-energy effective Hamiltonian are factorized into the product of the matrix elements of quark currents.
The decayB 0 s → D + s K − originating from b → cūs processes is a prime example where "QCD factorization" is expected to work excellently for the colour-allowed tree topologies [19][20][21][22][23]. We obtain the following amplitude: where G F is the Fermi constant, V * us V cb a factor of CKM matrix elements, f K the kaon decay constant and F Bs→Ds 0 (m 2 K ) a form factor parametrising the hadronic b → c quarkcurrent matrix element. The parameter describes the deviation from naive factorization: a DsK 1 characterises non-factorisable effects entering the colour-allowed tree amplitude T DsK , whereas E DsK denotes nonfactorisable exchange topologies.
The current state-of-the-art results within QCD factorization are found as |a 1 | ≈ 1.07 with a quasi-universal behaviour [12,24], with uncertainties at the percent level. Recently, even QED effects have been studied [25], which are small and fully included within the uncertainties.In Ref. [12], a theoretical analysis of theB 0 d → D + K − decay, which does not have an exchange topology, has been performed, yielding −0.0128 . The same numerical result is found forB 0 s → D + s π − , yet another pure colour-allowed tree decay which does not have an exchange contribution. A detailed discussion of the associated uncertainties is also given in [12]. TheB 0 s → D + s K − channel -the key player for our analysis -differs from theB 0 d → D + K − mode only through the spectator quark. We use the SU (3) flavour symmetry to relate the spectator quarks to each other and double the tiny error in view of possible SU (3)-breaking effects, and employ the following value for our analysis of the b → cūs transition: The exchange topology in theB 0 s → D + s K − decay is non-factorizable and cannot be reliably calculated from first principles. Consequently, we use experimental data to constrain this contribution. TheB 0 d → D + s K − transition originates only from an exchange topology E DsK , which differs from its counterpart E DsK inB 0 s → D + s K − through the down quark of the initialB 0 d meson. To be specific, we have the expression where τ B d and τ Bs denote the lifetimes of the B d and B s mesons, respectively, and Φ(x, y) is the usual phase-space function, depending on the meson masses. Using this expression with the measured branching ratio which offers direct access to the size of the exchange contribution. Another constraint is provided through the comparison of the branching ratio ofB 0 , yielding the following result [10]: which is remarkably consistent with (20). Due to the non-factorizable nature of the exchange amplitude, it may well have a large strong phase difference with respect to the colour-allowed tree amplitudes [3,11]. Interestingly, this feature is indicated by the comparison of (20) with (21), although the current uncertainties do not allow us to draw further conclusions. It is important to emphasize that no anomalous behaviour of the exchange topologies that could be caused by large rescattering or other non-factorizable effects is indicated by the data [10], as was also found in Ref. [11]. Our next step is to extract the |a 1 | parameters from the data in the cleanest possible way and to compare them with the theoretical SM predictions. In this respect, semileptonic decays provide a very useful tool [11,19,22]. In the case of theB 0 s → D + s K − channel, we have the partner decayB 0 s → D + s −ν , and introduce the ratio which takes the form with where Φ ph is a phase-space factor which is equal to 1 with excellent precision [10]. The CKM matrix element |V cb | cancels in R D + s K − . Moreover, due to the normalization condition F Bs→Ds 0 (0) = F Bs→Ds 1 (0), we have an essentially negligible impact of the nonperturbative hadronic form factors [26][27][28]. Using a recent LHCb measurement [29] of the differential rate of theB 0 s → D + We observe that this result is in tension with the theoretical prediction in (18) at the 2.2 σ level. Consequently, we indeed find another tension at the decay amplitude level with respect to the SM, as we would expect in view of the puzzling result for γ. This exciting situation is further strengthened by the fact that a similar pattern of the |a 1 | parameters -with surprisingly small values -arises also for other B (s) decays with similar dynamics. We have extracted these quantities from the data in an analogous way with semileptonic decay information [10], and show the results in the left panel of Fig. 1.
showing even a discrepancy of 4.8 σ. Puzzlingly small branching ratios for this channel and theB 0 d → D + d π − ,B 0 s → D + s π − modes were also pointed out in Refs. [11,12].
Let us next have a look at theB 0 s → K + D − s decay. The amplitude of this channel, which is caused by b → ucs processes, can be expressed in a way similar to the b → cūs case, as we discuss in detail in [10]. In analogy to theB 0 s → D + s K − mode, this transition is also a colour-allowed tree decay. However, as the roles of the heavy c and light u quarks are interchanged, the heavy-quark arguments which can be used to prove factorization up to tiny corrections for the b → c modes do not apply, and there may be larger nonfactorizable effects. In view of this less favourable theoretical situation, we will use the following range as a reference for our analysis: Here we have followed the QCD renormalisation group analysis in Ref. [30] as guidance, where |a 1 | = 1.01 ± 0.02 was found for the variation of the global |a 1 | parameter for colour-allowed tree decays with respect to variations of the renormalization scale and schemes. In Eq. (26), we allow for a five times larger uncertainty [10]. Interestingly, the experimental value of δ s in (5) is found in excellent agreement with factorization [2]. Since this strong phase difference characterizes the interference between the b → cūs and b → ucs decay paths, it supports factorisation -where such phases vanish -also in the b → ucs channel. In analogy Eq. (22), the partner decay for the clean extraction of the |a KDs 1 | parameter isB 0 s → K + ν . Although this channel has been observed by LHCb [31], a measurement of the differential rate has not yet been reported. Consequently, we have applied the SU (3) flavor symmetry and have utilised theB 0 d → π + −ν mode, for which we do have information from the BaBar and Belle collaborations [5,6]. We find as discussed in more detail in [10]. As in the b → cūs case, this result favours again a value smaller than our theoretical reference in Eq. (26). The result is illustrated in the panel on the right-hand side of Fig. 1, where we include also theB 0 d → π + D − s mode, which differs only through the spectator quark from theB 0 s → K + D − s channel. The current uncertainties are too large to draw further conclusions on these modes.
Within the SM, universal power-suppressed corrections of order Λ QCD /m b could in principle lead to a suppression of the |a 1 | parameters [24]. However, such effects would not allow us to accommodate the puzzling result for γ arising from the CP-violating observables of the B 0 s → D ∓ s K ± system. This phenomenon would require new sources for CP violation. The exciting possibility of NP effects in non-leptonic tree-level decays of B mesons was discussed in Refs. [32,33], and models for physics beyond the SM addressing the puzzles of the small branching ratios were studied in Refs. [13,14]. We consider these first models, which face challenges from direct NP searches at ATLAS and CMS [15], as interesting illustrations of specific scenarios. In the remainder of this paper, we will work with a model-independent parametrization of NP contributions.

New Physics Analysis
Let us now extend the analysis of the B 0 s → D ∓ s K ± system to include NP effects. As we have already noted, the NP contributions to accommodate the puzzling measurement of γ and the |a 1 | values would have to enter at the decay amplitude level. For the b → cūs quark-level transition, we generalise the decay amplitude as with the NP parameterρ whereφ andδ denote CP-violating and CP-conserving phases, respectively. A similar expression can be written for theb →ūcs transition, with parameters ρ, ϕ and δ.
In order to get first access to these NP pazs, we compare the SM predictions for the branching ratios with the corresponding experimental values. Let us first have a look at the b → cūs transition, where we introducē Here B(B 0 s → D + s K − ) th is the CP average of the theoretical branching ratios, converted from the experimental measurements. For the determination ofb, we utilize again the semileptonic ratio in (22), which is particularly clean with respect to form factor and CKM parameter uncertainties. In the presence of NP contributions with new sources of CP violation, we introduce the following generalized ratio [10]: which satisfies R DsK = R D + s K − in the case of vanishing direct CP asymmetries, as in the SM. Then we obtain the following expression forb in terms of The parameter |a DsK 1 eff | given in (17) is the product of the theoretical prediction of |a DsK 1 | in (18) obtained within QCD factorization and the parameter r DsK E in (21), which was constrained through experimental data. We obtain the value |a DsK 1 eff | = 1.07 ± 0.09, which we will use in the numerical analysis below. The small impact of the exchange topology following from the experimental data as reflected by Eq. (20) holds irrespectively of whether we have the SM or possible NP contributions to this topology. In analogy, for theb →ūcs transition, we introduce Making use of r KDs E = 1.00 ± 0.08, which follows from an analysis similar to the one for r DsK E given above [10], and the reference for |a KDs 1 | in Eq. (26), we find |a KDs For the observableξ, we obtain the generalisation Similarly, we rewrite ξ with ∆ϕ, where we interchange the NP parametersρ,δ,φ and ρ, δ, ϕ.
The product ξ ×ξ introduced in Eq. (2), which plays the key role for CP violation in the SM, is generalized as follows: In contrast to the SM, NP may generate non-vanishing direct CP asymmetries: An analogous expression holds for the CP-conjugate asymmetryĀ dir CP , involvingρ with δ andφ. We may then write the first ratio entering (37) as 1 + ρ e iδ e +iϕ 1 + ρ e iδ e −iϕ = e −i∆Φ 1 − A dir where tan ∆Φ = − 2ρ cos δ sin ϕ + ρ 2 sin 2ϕ 1 + 2ρ cos δ cos ϕ + ρ 2 cos 2ϕ ; the second ratio takes a similar form, involvingρ,φ andĀ dir CP with a phase ∆Φ. We then obtain where generalising the SM relations in (3). Finally, we arrive at which is a theoretical clean relation playing key role in our analysis. Here, the UT angle γ enters as the "effective" angle Consequently, (6) actually corresponds to γ eff . Note that in combined fits to the data for various B decays to extract γ, such as in Ref. [9], NP effects may average out to some extend, thereby yielding an effective angle with NP contributions which -in contrast to (44) -cannot transparently be quantified. It will rather be crucial to search for patterns in the individual γ determinations, aiming at the highest precision.
Let us now apply our formalism to the current data. In order to be consistent with the LHCb assumption, we set the strong phases δ andδ to 0 • . This implies vanishing direct CP asymmetries, in agreement with B → DK data within the uncertainties [6].As we have noted after Eq. (31), using these assumptions, we may identify the experimental values for R D + s K − = 0.05 ± 0.01 and R K + D − s = 3.64 ± 1.70, with R DsK and R KDs , respectively. Complementing these values with the other relevant parameters introduced above, we findb = 0.58 ± 0.16, b = 0.50 ± 0.26.
The numerical value refers to γ = (70 ± 7) • , which is consistent with UT analyses, and the result in (6). In order to convert the measured observables into constraints on the NP parameters, we first employb and b to determineρ and ρ as functions ofφ and ϕ, respectively. Using then (6), we may calculate ϕ as a function ofφ, fixing a contour in the ϕ-φ plane. Finally, using againρ(φ) and ρ(ϕ) allows us to calculate a correlation in theρ-ρ plane, where each point is linked withφ and ϕ.
In Fig. 2, we show the corresponding correlations for the central values of the current data. In Fig. 3, we show the impact of the uncertainties of the input quantities ∆ϕ, b andb, varying them separately. We can nicely see that the SM point corresponding to the origin in theρ-ρ plane is excluded, and notice that NP contributions with CPviolating phases are simultaneously required in the b → cūs andb →ūcs decay paths.  Interestingly, we can describe the measurements with NP contributions as small as about 30% of the SM amplitudes.
In view of the complexity of the strategy, we finally summarise the main steps: • Step 1: CP Violation Utilising C, S, A ∆Γ and their CP conjugates, we determine ξ andξ, respectively, unambiguously from the data. The product ξ ×ξ, generalised to include NP, allows a theoretically clean determination of γ eff ≡ γ + γ NP , where γ NP is a function of the NP parameters. Using information on γ from other processes, we extract γ NP .

• Step 2: Branching Ratio Information
To have a particularly clean setup, we combine the branching fractions of the nonleptonic decays with differential rate information from their semi-leptonic partners. We constrain the exchange topologies via other control channels. Complementing data with theoretical input for the |a 1 | parameters, we extract b andb.

• Step 3: NP Parameters Correlation
We make use of all three observables γ eff , b andb and explore the available space for NP by obtaining correlations between the NP parameters ρ(ϕ) andρ(φ).  Figure 4: Illustration of the two puzzles and the strategy we presented in this paper.

Concluding Remarks
TheB 0 s → D + s K − andB 0 s → D − s K + decays with their CP conjugates are key players in the testing of the SM. We have demonstrated that the intriguing picture arising from the measured CP-violating observables, which results in a tension with the SM at the 3σ level, is complemented by a puzzling pattern of the individual branching ratios. The latter finding is actually in accordance with measurements of rates of B (s) modes with similar dynamics, where we find tensions with up to 4.8σ significance. While the γ measurement cannot be explained through non-factorizable effects, the branching ratio puzzles could in principle be accommodated through such contributions. We would like to stress that the experimental result of the strong phase difference δ s between the b → cūs and b → ucs decay paths is in excellent agreement with factorization. It is exciting to reveal and link these new puzzles, complementing indications of NP in other corners of the flavour sector, where currently rare decays arising from b → s + − quark-level processes involving also leptons are in the spotlight [34].
We have presented a model-independent description of the B 0 s → D ∓ s K ± system to include NP effects and give a generalized expression for an effective angle γ in terms of NP parameters that can be determined in a theoretically clean way from the measured CPviolating observables. Utilizing furthermore the information from the decay branching ratios, the NP parameters can be constrained, as we have illustrated for the current data. In Fig. 4, we give a schematic overview. In the future high-precision B physics era, this formalism can be fully exploited through more sophisticated experimental analyses, searching for direct CP violation and measuring the experimental B s branching ratios for the separate D + s K − and D − s K + final states. The measurement of the differential B 0 s → K + ν rate would be another important ingredient to complement the analysis.
The NP parameters resulting from this analysis will serve as benchmarks for the model building community. It will be exciting to monitor the evolution of the data and the sharper picture emerging from the application of our strategy. The central question is whether these studies will finally allow us to establish the presence of new CP-violating contributions to B 0 s → D ∓ s K ± decays, thereby raising the question of whether such new sources of CP violation could open a window to understand the matter-antimatter asymmetry of the Universe.