We combine all the available experimental information on B_{
s
} mixing, including the very recent tagged analyses of B_{
s
} → J/Ψϕ by the CDF and DØ collaborations. We find that the phase of the B_{
s
} mixing amplitude deviates more than 3σ from the Standard Model prediction. While no single measurement has a 3σ significance yet, all the constraints show a remarkable agreement with the combined result. This is a first evidence of physics beyond the Standard Model. This result disfavours New Physics models with Minimal Flavour Violation with the same significance.

PACS Codes: 12.15.Ff, 12.15.Hh, 14.40.Nb

1. Letter

In the Standard Model (SM), all flavour and CP violating phenomena in weak decays are described in terms of quark masses and the four independent parameters in the Cabibbo-Kobayashi-Maskawa (CKM) matrix [1, 2]. In particular, there is only one source of CP violation, which is connected to the area of the Unitarity Triangle (UT). A peculiar prediction of the SM, due to the hierarchy among CKM matrix elements, is that CP violation in B_{
s
} mixing should be tiny. This property is also valid in models of Minimal Flavour Violation (MFV) [3–8], where flavour and CP violation are still governed by the CKM matrix. Therefore, the experimental observation of sizable CP violation in B_{
s
} mixing is a clear (and clean) signal of New Physics (NP) and a violation of the MFV paradigm. In the past decade, B factories have collected an impressive amount of data on B_{
d
} flavour- and CP-violating processes. The CKM paradigm has passed unscathed all the tests performed at the B factories down to an accuracy just below 10% [9–11]. This has been often considered as an indication pointing to the MFV hypothesis, which has received considerable attention in recent years. The only possible hint of non-MFV NP is found in the penguin-dominated b → s non-leptonic decays. Indeed, in the SM, the coefficient of the time-dependent CP asymmetry in these channels is equal to the measured with decays, up to hadronic uncertainties related to subleading terms in the decay amplitudes. Present data show a systematic, although not statistically significant, downward shift of with respect to [12–21], while hadronic models predict a shift in the opposite direction in many cases [22–29].

From the theoretical point of view, the hierarchical structure of quark masses and mixing angles of the SM calls for an explanation in terms of flavour symmetries or of other dynamical mechanisms, such as, for example, fermion localization in models with extra dimensions. All such explanations depart from the MFV paradigm, and generically cause deviations from the SM in flavour violating processes. Models with localized fermions [30–32], and more generally models of Next-to-Minimal Flavour Violation [33], tend to produce too large effects in ε_{
K
} [34, 35]. On the contrary, flavour models based on nonabelian flavour symmetries, such as U(2) or SU(3), typically suppress NP contributions to s ↔ d and possibly also to b ↔ d transitions, but easily produce large NP contributions to b ↔ s processes. This is due to the large flavour symmetry breaking caused by the top quark Yukawa coupling. Thus, if (nonabelian) flavour symmetry models are relevant for the solution of the SM flavour problem, one expects on general grounds NP contributions to b ↔ s transitions. On the other hand, in the context of Grand Unified Theories (GUTs), there is a connection between leptonic and hadronic flavour violation. In particular, in a broad class of GUTs, the large mixing angle observed in neutrino oscillations corresponds to large NP contributions to b ↔ s transitions [36–39].

In this Letter, we show that present data give evidence of a B_{
s
} mixing phase much larger than expected in the SM, with a significance of more than 3σ. This result is obtained by combining all available experimental information with the method used by our collaboration for UT analyses and described in Ref. [40].

We perform a model-independent analysis of NP contributions to B_{
s
} mixing using the following parametrization [41–46]:

(1)

where is the effective Hamiltonian generated by both SM and NP, while only contains SM contributions. The angle β_{
s
} is defined as and it equals 0.018 ± 0.001 in the SM (we are using the usual CKM phase convention in which is real to a very good approximation).

We use the following experimental input: the CDF measurement of Δm_{
s
} [47], the semileptonic asymmetry in B_{
s
} decays [48], the dimuon charge asymmetry from DØ [49] and CDF [50], the measurement of the B_{
s
} lifetime from flavour-specific final states [51–59], the two-dimensional likelihood ratio for ΔΓ_{
s
} and from the time-dependent tagged angular analysis of B_{
s
} → J/ψϕ decays by CDF [60] and the correlated constraints on Γ_{
s
}, ΔΓ_{
s
} and ϕ_{
s
} from the same analysis performed by DØ [61]. For the latter, since the complete likelihood is not available yet, we start from the results of the 7-variable fit in the free-ϕ_{
s
} case from Table one of ref. [61]. We implement the 7 × 7 correlation matrix and integrate over the strong phases and decay amplitudes to obtain the reduced 3 × 3 correlation matrix used in our analysis. In the DØ analysis, the twofold ambiguity inherent in the measurement (ϕ_{
s
} → π - ϕ_{
s
}, ΔΓ_{
s
} → - ΔΓ_{
s
}, cos δ_{1,2} → - cos δ_{1,2}) for arbitrary strong phases was removed using a value for cos δ_{1,2} derived from the BaBar analysis of B_{
d
} → J/ΨK* using SU(3). However, the strong phases in B_{
d
} → J/ΨK* and B_{
s
} → J/Ψϕ cannot be exactly related in the SU(3) limit due to the singlet component of ϕ. Although the sign of cos δ_{1,2} obtained using SU(3) is consistent with the factorization estimate, to be conservative we reintroduce the ambiguity in the DØ measurement. To this end, we take the errors quoted by DØ as Gaussian and duplicate the likelihood at the point obtained by applying the discrete ambiguity. Indeed, looking at Fig. 2 of ref. [61], this seems a reasonable procedure. Hopefully DØ will present results without assumptions on the strong phases in the future, allowing for a more straightforward combination. Finally, for the CKM parameters we perform the UT analysis in the presence of arbitrary NP as described in ref. [34], obtaining , = 0.384 ± 0.035 and sin 2β_{
s
} = 0.0409 ± 0.0038. The new input parameters used in our analysis are summarized in Table 1, all the others are given in Ref. [34]. The relevant NLO formulae for ΔΓ_{
s
} and for the semileptonic asymmetries in the presence of NP have been already discussed in refs. [34, 62, 63].

We also show the correlation coefficients Cs of the measurements of ϕ_{
s
}, ΔΓ_{
s
} and from ref. [61].

The results of our analysis are summarized in Table 2. We see that the phase deviates from zero at 3.7σ. We comment below on the stability of this significance. In Fig. 1 we present the two-dimensional 68% and 95% probability regions for the NP parameters and , the corresponding regions for the parameters and , and the one-dimensional distributions for NP parameters. Notice that the ambiguity of the tagged analysis of B_{
s
} → J/Ψϕ is slightly broken by the presence of the CKM-subleading terms in the expression of Γ_{12}/M_{12} (see for example eq. (5) of ref. [63]). The solution around corresponds to and . The second solution is much more distant from the SM and it requires a dominant NP contribution (). In this case the NP phase is thus very well determined. The strong phase ambiguity affects the sign of cos ϕ_{
s
} and thus Re , while Im in any case.

Table 2

Fit results for NP parameters, semileptonic asymmetries and width differences.

Observable

68% Prob.

95% Prob.

-19.9 ± 5.6

[-30.45,-9.29]

-68.2 ± 4.9

[-78.45,-58.2]

1.07 ± 0.29

[0.62,1.93]

-51 ± 11

[-69,-27]

-79 ± 3

[-84,-71]

0.73 ± 0.35

[0.24,1.38]

1.87 ± 0.06

[1.50,2.47]

Im

-0.74 ± 0.26

[-1.54,-0.30]

Re

-0.13 ± 0.31

[-0.61,0.78]

-1.82 ± 0.28

[-2.68,-1.36]

-0.34 ± 0.21

[-0.75,0.03]

-2.1 ± 1.0

[-4.7,-0.3]

ΔΓ_{
s
}/Γ_{
s
}

0.105 ± 0.049

[0.02,0.20]

-0.098 ± 0.044

[-0.19,-0.02]

Whenever present, we list the two solutions due to the ambiguity of the measurements. The first line corresponds to the one closer to the SM.

Before concluding, we comment on our treatment of the DØ result for the tagged analysis and on the stability of the NP fit. Clearly, the procedure to reintroduce the strong phase ambiguity in the DØ result and to combine it with CDF is not unique given the available information. In particular, the Gaussian assumption can be questioned, given the likelihood profiles shown in Ref. [61]. Thus, we have tested the significance of the NP signal against different modeling of the probability density function (p.d.f.). First, we have used the 90% C.L. range for ϕ_{
s
} = [-0.06, 1.20]° given by DØ to estimate the standard deviation, obtaining ϕ_{
s
} = (0.57 ± 0.38)° as input for our Gaussian analysis. This is conservative since the likelihood has a visibly larger half-width on the side opposite to the SM expectation (see Fig. 2 of Ref. [61]). Second, we have implemented the likelihood profiles for ϕ_{
s
} and ΔΓ_{
s
} given by DØ, discarding the correlations but restoring the strong phase ambiguity. The likelihood profiles include the second minimum corresponding to ϕ_{
s
} → ϕ_{
s
}+π, ΔΓ → -ΔΓ, which is disfavoured by the oscillating terms present in the tagged analysis and is discarded in our Gaussian analysis. Also this approach is conservative since each one-dimensional profile likelihood is minimized with respect to the other variables relevant for our analysis. It is remarkable that both methods give a deviation of from zero of 3 σ (the 3 σ ranges for are [-88, -48]° ∪ [-41, 0]° and [-88, 0]° for the two methods respectively). We conclude that the combined analysis gives a stable evidence for NP, although the precise number of standard deviations depends on the procedure followed to combine presently available data.

To illustrate the impact of the experimental constraints, we show in Fig. 2 the p.d.f. for obtained without the tagged analysis of B_{
s
} → J/Ψϕ or including only CDF or DØ results. Including only the CDF tagged analysis, we obtain at 97.7% probability (2.3σ). For DØ, we show results obtained with the Gaussian and likelihood profile treatment of the errors. In the Gaussian case, the DØ tagged analysis gives at 98.0% probability (2.3σ), while using the likelihood profiles at 92.8% probability (1.8σ). Finally, it is remarkable that the different constraints in Fig. 2 are all consistent among themselves and with the combined result. We notice, however, that the top-left plot is dominated by the measurement of while favours positive , although with a very low significance. For completeness, in Table 2 we also quote the fit results for , and for ΔΓ_{
s
}/Γ_{
s
}.

In this Letter we have presented the combination of all available constraints on the B_{
s
} mixing amplitude leading to a first evidence of NP contributions to the CP-violating phase. With the procedure we followed to combine the available data, we obtain an evidence for NP at more than 3σ. To put this conclusion on firmer grounds, it would be advisable to combine the likelihoods of the tagged B_{
s
} → J/Ψϕ angular analyses obtained without theoretical assumptions. This should be feasible in the near future. We are eager to see updated measurements using larger data sets from both the Tevatron experiments in order to strengthen the present evidence, waiting for the advent of LHCb for a high-precision measurement of the NP phase.

It is remarkable that to explain the result obtained for ϕ_{
s
}, new sources of CP violation beyond the CKM phase are required, strongly disfavouring the MFV hypothesis. These new phases will in general produce correlated effects in ΔB = 2 processes and in b → s decays. These correlations cannot be studied in a model-independent way, but it will be interesting to analyse them in specific extensions of the SM. In this respect, improving the results on CP violation in b → s penguins at present and future experimental facilities is of the utmost importance.

Declarations

Authors’ Affiliations

(1)

CERN

(2)

INFN, Sezione di Roma Tre

(3)

INFN, Sezione di Roma

(4)

Dipartimento di Fisica, Università di Roma Tre

(5)

Dipartimento di Fisica, Università di Roma "La Sapienza"

(6)

Dipartimento di Fisica, Università di Genova and INFN

(7)

ETH Zurich

(8)

Laboratoire de l'Accélérateur Linéaire, IN2P3-CNRS et Université de Paris-Sud

Bona M, Ciuchini M, Franco E, Lubicz V, Martinelli G, Parodi F, Pierini M, Roudeau P, Schiavi C, Silvestrini L, Stocchi A, [UTfit Collaboration]: JHEP. 2005, 0507:028.View ArticleADS

Bona M, Ciuchini M, Franco E, Lubicz V, Martinelli G, Parodi F, Pierini M, Roudeau P, Schiavi C, Silvestrini L, Stocchi A, Vagnoni V, [UTfit Collaboration]: JHEP. 2006, 0610:081.

Charles J, Hocker A, Lacker H, Laplace S, Le Diberder FR, Malcles J, Ocariz J, Pivk M, Roos L, [CKMfitter Group]: Eur Phys J C. 2005, 41:1.View ArticleADS

Aubert B, [BABAR Collaboration], et al. arXiv:hep-ex/0607101

Aubert B, [BABAR Collaboration], et al.: Phys Rev Lett. 2007, 98:031801.View ArticleADS

Aubert B, [BABAR Collaboration], et al.: Phys Rev D. 2007, 76:071101.View ArticleADS

Aubert B, [BABAR Collaboration], et al.: Phys Rev D. 2007, 76:091101.View ArticleADS

Aubert B, [BABAR Collaboration], et al.: Phys Rev Lett. 2007, 99:161802.View ArticleADS

Aubert B, [BABAR Collaboration], et al.: Phys Rev D. 2008, 77:012003.View ArticleADS

Bona M, Ciuchini M, Franco E, Lubicz V, Martinelli G, Parodi F, Pierini M, Roudeau P, Schiavi C, Silvestrini L, Stocchi A, Vagnoni V, [UTfit Collaboration]: JHEP. 2006, 0603:080.

Bona M, Ciuchini M, Franco E, Lubicz V, Martinelli G, Parodi F, Pierini M, Roudeau P, Schiavi C, Silvestrini L, Stocchi A, Vagnoni V, [UTfit Collaboration]: Phys Rev Lett. 2006, 97:151803.View ArticleADS

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.