The signal of ill-defined CPT weakening entanglement in the B d system.

In the presence of quantum-gravity fluctuations (space-time foam), the CPT operator may be ill-defined. Its perturbative treatment leads to a modification of the Einstein-Podolsky-Rosen correlation of the neutral meson system by adding an entanglement-weakening term of the wrong exchange symmetry, the ω -effect. In the current paper we identify how to probe the complex ω in the entangled B d -system using the flavour (f)-CP(g) eigenstate decay channels: the connection between the intensities for the two time-ordered decays (f, g) and (g, f) is lost. Appropriate observables are constructed allowing independent experimental determinations of Re( ω ) and Im( ω ), disentangled from CPT violation in the evolution Hamiltonian Re( θ ) and Im( θ ). 2 σ tensions for both Re( θ ) and Im( ω ) are shown to be uncorrelated.

Introduction

The physics of discrete symmetries in particle and nuclear physics has always been a fascinating subject, since the observation of CP violation in the neutral kaon system [1], which was a clear experimental surprise and set the scene for subsequent precision tests of such discrete symmetries in other systems, including entangled neutral meson factories. Today CP violation in the K and systems, as well as T violation with entangled ’s [2], have been demonstrated experimentally to great accuracy. However, their combination CPT remains unbroken. This is believed to be due to one of the crucial theorems of modern physics, ensuring CPT Invariance of quantum field theory models that are Lorentz invariant, local (in their interactions) and unitary (that is, they conserve probability) [3]. This is basically a theorem of flat space-time. Quantum gravity or in general deviations from any of the three assumptions may lead to (independent) violations of CPT, which, if observed in nature, would undoubtedly constitute an indication of completely novel physics. Having mentioned quantum gravity, it is worth recalling a corollary by Wald [4], according to which a potential decoherence induced during observations in local scattering experiments in which the experimenter has no access to microscopic quantum gravity degrees of freedom, may lead to an effectively ill-defined CPT quantum mechanical operator. This observation prompted the authors of [5] to introduce a different observable for this kind of decoherence-induced CPT violation, termed -effect. The -effect is different from the situation where CPT violation is violated in the effective Hamiltonian, parameterized by the complex parameter. Among other possible sources, can be due to, e.g., Lorentz violation [6, 7] as a result of propagation in some Lorentz violating space-time (or otherwise) backgrounds. In the latter case the quantum mechanical operator that implements CPT symmetry is well defined but simply does not commute with the Hamiltonian. The -effect, if observed, points to an observation of a phenomenon that is exclusively linked to ill-defined nature of the CPT operator, which to date is theoretically linked only to fundamental decoherence [4], independently of any violation of CPT in the Hamiltonian.

Recently, a study for separate direct evidence of T, CP, CPT symmetry violation was accomplished [8]. It was based on the precise identification of genuine asymmetry parameters in the time evolution of intensities between the two decays in a B-Factory of entangled neutral -meson states. Their values were obtained from the BaBar measurements [2] of the different flavour-CP eigenstate decay channels. The concept, put forward in [9, 10], uses the entangled character of the initial state as the crucial ingredient to (i) connect experimental double decay rates with specific meson transition probabilities and (ii) identify the transformed transition to that taken as a reference [11, 12]. Possible fake effects [8] were demonstrated to be well under control by measurements in the same experiment. The methodology, discussed in [9, 10], appears to be [13] crucially dependent on the assumed maximal entanglement between and , or between two orthogonal superpositions of them, as given by the Einstein–Podolsky–Rosen (EPR) correlation [14] imposed by their decay from the -state with . The corresponding antisymmetric state of the system has two important implications: (i) the program of using entanglement and the decays as filtering measurements to prepare and detect the meson states can be implemented at any time for the first decay, even in the presence of mixing during the previous entangled evolution; (ii) the coefficients of the different time-dependent terms in the double decay rate intensities for the time-ordered decays to (g, f) are related to those for the time-ordered decays to (f, g). The antisymmetry of the entangled state is kept for any two independent states of the neutral mesons, so its evolution leads to a trivial time dependence with definite symmetry under the combined exchange . As a consequence, the double decay rate intensity (see Eq. 14 below) satisfies for the coefficients of its time dependence with ,where the time-ordered decays (f, g) and (g, f) are, in general, not connected by any symmetry transformation. At this level, they can be considered as two different experimental ways of measuring the same quantity when .

In the application to definite flavour or CP eigenstates decay products, the preparation by maximal entanglement of the initial state of a single neutral meson is usually referred to as “flavour tagging” , , or “CP tagging” , . The underlying assumption considers , as two states of the same field, in order to impose Bose statistics with charge conjugation C and permutation with C , and it may be invalidated if the CPT operator cannot be intrinsically well defined, as mentioned above. This latter circumstance may occur, for example, in the context of an extended class of quantum-gravity models, where the structure of quantum space-time at Planckian scales ( m) may actually be fuzzy, characterized by a “foamy” nature (space-time foam) [5, 15, 16]. Let us emphasize once more that this kind of CPT breaking is different from an explicit CPT violation in the Hamiltonian dynamics such that [CPT, H] , as conventionally introduced, in the context of the Weisskopf–Wigner approach [17-19] for the neutral meson system, in the mass matrix. This last CPT violation does not invalidate the analysis followed in [8] and, in fact, genuine observables for CPT violation were found with their values obtained from experiment. However, the CPT breaking associated to “ill-defined” particle–antiparticle states modifies the EPR correlation, producing the aforementioned -effect [5, 20, 21]. Treating it in perturbation theory, in such a way that we still talk the language of , , the perturbed two-particle state will contain a component of the “wrong” symmetry at the instant of their production by the decay of :where is a complex CPT-breaking parameter [5, 20], associated with the non-identical particle nature of the neutral meson and antimeson states. The presence of an -effect weakens the entanglement of the initial state (2), as follows from the fact that when the state simply reduces to a product state, whilst when the state is fully entangled.

We emphasize that the modification in Eq. (2) is due to the loss of indistinguishability of and and not due to violation of symmetries in the production process. Evidently, the probabilities for the two states connected by a permutation are different due to the presence of . This modification of the initial state vector has far-reaching consequences for the concept of meson tagging and for the relation of the time-dependent intensities between the decays to time-ordered (f, g) and (g, f) channel.1

In what follows we will study the non-trivial time evolution of Eq. (2), in the simplified but physically relevant case of a time-independent , in order to (i) establish the appearance of terms of the (previously forbidden) type and , and (ii) introduce a set of observables, which actually serve as a direct way for measuring , based on the violation of the relations of Eq. (1), i.e. using as observables for :and checking experimentally the robustness of the correlation between the two states assumed during the tagging. This paper demonstrates that the comparison between the double decay rate Intensities for time-ordered ( flavour, CP) eigenstate decay products and (g, f) is sensitive to both and .

Time evolution

Double decay rates, time-dependent intensities

The eigenstates of the effective Hamiltonian are2 In terms of thesewhere is a CP and CPT violating complex parameter given by . The time evolution of two-meson flavour states iswhere , ,andThe , matrices are In Eqs. (10)–(12), is the usual meson mixing quantity given by . Before addressing actual observables, it is worth noting that, considering Eq. (7), it is clear that the presence of the symmetric state in Eq. (2) induces the appearance of and states.

The transition amplitude for the decay of the first state into at time , and then the second state into at time is . Squaring and integrating over , the double decay rate I(f, g; t) is obtained:Expanding to first order in , and taking , I(f, g; t) has the following form for generic f and g decay channels3:with where, in terms of the decays amplitudes and , the following parameters are used4: In addition, and . It is worth recalling that, for the flavour-specific decay channels (“” for short in the following), we have , .

Sensitivity to

Coming back to the transition amplitude , it has the following structure:The prefactor gives a global dependence in . One can readily observe that the -dependent terms, even for (i.e. already for the leading contribution), do introduce an additional non-trivial dependence. Ignoring that prefactor, it is clear that combining the transformations and , the first contribution, the standard one, just receives a sign. This implies that, in the absence of , in the t-dependence of I(f, g; t),we necessarily have [8]: , and .

In the presence of the situation changes drastically. From the remaining contributions in Eq. (20), the ones induced by the evolution of the -dependent term in Eq. (2), the situation is more involved: the first one, proportional to and -independent, is clearly invariant under the combination of and . The last two terms are separately invariant under , but have no well-defined transformation under ; moreover, contrary to the previous contributions, they depend on , the time elapsed between production of the pair and the first decay.5 Out of those properties, the simple assignment of symmetry/antisymmetry under to the t-even/t-odd terms in , possible when , does not apply when . This simple remark provides the first understanding of the potential sensitivity to the presence of : while in the absence of , the measurement of intensities for decays into f and g with the two different orderings (i) first f then g and (ii) first g then f, provides two experimentally independent measurements of the same theoretical quantities, in the presence of the situation has changed. Deviations from the standard symmetry properties are a gateway to probe for .

The BaBar collaboration performed separate analyses [2] for the two different time orderings of the two B meson decays. Previous studies, like [21], exploited the use of two flavour specific decay channels to obtain bounds on through the appearance of and states for . Equation (14) shows that, using flavour specific channels alone, there is no sensitivity to : since , the terms in would be absent.6 Fortunately enough, besides addressing the two different time orderings, in [2], one decay is flavour specific (labeled ), while the other is CP specific (decays into , labeled for short): sensitivity to both and is thus expected.

Experimental observables

In order to reduce experimental uncertainties in the different channels, the BaBar collaboration, in Ref. [2], fixed the constant term and measured the coefficients and of the decay intensityusing for the f and g states one flavour specific channel, or , and one CP eigenstate, or . Obviously we should havewhere one should remember that in the coefficients and , the ordering of f and g means that f corresponds to the first (in time) decay product of the entangled state evolved in time, and g corresponds to the second (in time) decay product. In the case under consideration, the flavour specific decays simplify significantly the expressions, which are, at linear order in , , In the presence of , the time ordering definite symmetry is not valid anymore and therefore it is relevant to write the completely different coefficients As anticipated, and are linear in , and thus the fact that the BaBar collaboration distinguished the different decay time orderings in [2], now reveals crucial to disentangle the effect: These combinations are linearly sensitive not only to but also to when . The sensitivity to depends critically on the use of a CP eigenstate channel with large , as is the case with and .

Before we present our numerical results it is important to summarize the assumptions we have made in our analysis. In the limit, as explained after Eq. (21), we also have the set of equationstherefore it is clear that any deviation from those zeros, as identified by Eqs. (28) and (29), will give a positive signal for . This result is independent of any additional assumption as regards the decay amplitudes, it is valid even if we consider any pair of decay channels f and g. This should not come as a surprise because we compare two double decay rates to the same final states: these observables only reverse the order in time of the “same” decay channels and the results are dictated by the symmetry properties of the initial state entering Eq. (20).

Independently of the previous remark it is important to stress that in the present analysis there are no assumptions on the decay amplitudes in the channels . We follow here the parameterisation used in Ref. [8]. In that work, in order to point out genuine T violating observables independently of any assumption as regards CP or CPT, one is forced to allow for CPT violation both in the – mixing and in the decay amplitudes. One also needs to have a completely general parameterisation in order to measure any possible contamination coming from these contributions and others explained in [8]. Therefore we have parameterized and in terms of two arbitrary complex numbers that later on become four real arbitrary parameters , , and to be extracted from the fit itself ( and only add two signs). This parameterisation fits any New Physics including Standard Model penguins, New Physics models or CPT violation in the decay amplitudes. Nevertheless we have worked under the following approximations, which we deem to be physically reasonable: (i) we have made the approximation in order to use the Babar data, (ii) we have not allowed for wrong sign decays in flavour specific channels with the following implications: , as indicated at the end of Sect. 2.1. Note that this last assumption does not invalidate the fact that a difference between and (or between and ) would be an unambiguous signal of .

Results

We are now ready to present the results obtained from a global fit to available BaBar experimental data, following the same statistical treatment as in Ref. [8]. We use the 16 experimental observables measured by BaBar in [2]: , , and . Taking into account full covariance information on statistical and systematic uncertainties, we perform a fit in terms of the set of parameters , , , , , , , , , with the known constraints . Therefore we generalize the corresponding fit presented in Ref. [8] to the actual situation where deviations from EPR entanglement are present due to the -effect [5]. A more restricted fit is also done in the case where no wrong sign flavour decays are allowed in the decays, that is with .

Fig. 3

Correlations among and . (Colour coding as indicated in Fig. 1.) a vs. , b vs. , c vs. , d vs.

Fig. 4

and vs. ; the presence of affects the extracted values of and . (Colour coding as indicated in Fig. 1.) a vs. , b vs.

Fig. 1

vs. in the general fit (blue regions with solid contours), and in the fit with (red regions with dashed contours); darker to lighter regions correspond to two-dimensional 68, 95 and 99% CL. Figures 3 and 4 obey the same colour coding for the two fits and the CL regions

In Table 1(I) we present the general result of the fit, whose most salient features are the following:In Table 1(II) we present the results of the same fit with the additional requirement of not having wrong flavour decays, . No significant differences were noticed with respect to the conclusions discussed above for the general case. For completeness we show, when relevant, both analyses together in the same plots without further comments.

Table 1

Summary of results

(I) Parameters – general analysis
Reθ	±(6.11±3.45)10^-2	Imθ	(0.99±1.98)10^-2
Reω	(1.09±1.60)10^-2	Imω	±(6.40±2.80)10^-2
S_K_S	-0.624±0.030	R_K_S	±(0.781±0.024)
C_K_S	-1.44±3.2810^-2
S_K_S+S_K_L	(3.7±4.9)10^-2	R_K_S+R_K_L	(-3.27±4.3)10^-2
C_K_S-C_K_L	(-6.8±6.3)10^-2
(II) Parameters – λ_K_S+λ_K_L=0 analysis
Reθ	±(3.10±1.51)10^-2	Imθ	(0.14±1.67)10^-2
Reω	(1.17±1.59)10^-2	Imω	±(5.46±2.70)10^-2
S_K_S	-0.640±0.025	R_K_S	±(0.769±0.022)
C_K_S	(1.61±1.88)10^-2

Experimental data – more precisely the BaBar measurements in [2] – are sensitive for the first time to , revealing a tantalizing deviation from . These observables are also sensitive to , but they do not show any significant deviation from , and the previous determination [21] – using semileptonic channels – is still better than the present one.

The results of the fit for the CPT violating parameter – in the evolution Hamiltonian – are compatible with the previous determination in [8] and the one performed by the BaBar collaboration in Ref. [29]. An exciting effect in is still present.

The parameters that measure the presence of wrong flavour decays in , i.e. , and , do not show any significant deviation from zero and the results are consistent with [8].

In the case of and we observe that they differ by more than with respect to the determination in [8] without including the effect. Should this persist in the future, it could affect the precise determination of the unitarity triangle angle .

In Fig. 1 is shown the result for the new parameters not previously considered in the analyses where EPR entangled initial states where assumed. A deviation of the complex number from zero is found at 95% confidence level. This deviation comes essentially from and it represents a measurement of this parameter for the first time; the measurement of does not improve on the value obtained previously [21] from flavour specific decays.

The stability of the fitted value of the complex CPT violating parameter is shown in Fig. 2a, and b, where it is clear that the results for and do not change from the constrained case to the general case with arbitrary . Cross correlations among the different components of and are shown in Fig. 3. For example, Fig. 3c shows the independence of and : furthermore one can see in that figure that the point (0, 0) in this projection is at more than from the best fit values (or even at in the constrained analysis).

Fig. 2

Comparison of the results with and without ; vs. regions are shown in different scenarios: blue regions with solid contours correspond to fits including the -effect while red regions with dashed contours correspond fits without , i.e. with . Panel a shows the results for the general analyses while panel b shows the results for the analyses with . a Complete fit. b Restricted

Finally, in Fig. 4, one can see the near linear correlation among and . This explains why the presence of affects both and .

Before closing this section we would like to comment on an important issue that allows for the genuine CPTV -effect to be disentangled from ‘fake’ effects that mimick and exist already within conventional Standard Model physics. The latter are background contributions associated with processes given in Fig. 5a–c, in which the pair is in a configuration. However, they have a invariant mass dependence that is distinct from that of the genuine effect, which involves an amplitude coherent with the state, and therefore corresponds to a resonant peak in the invariant mass distribution. A similar situation characterizes the entangled neutral kaon system [5].

Fig. 5

Potential background contributions to the effect

In particular, the process of Fig. 5a associated with an intermediate state, has been discussed by Colangelo et al. [23], who performed a detailed calculation obtaining that the fake effect is of order , and concluded that a CPT test can be done at B factories [24].

The other two processes in Fig. 5b, c, are also expected to be very small; this can be seen, for example, for the absorptive part of Fig. 5b, by scaling the analog calculation in [25] to the – system. Similar processes have been analysed in factories for the – system [5, 26–28]. In general, these contributions are expected to be far below the present sensitivity.

Conclusions

In the present article we have discussed the possibility of probing the entanglement-weakening CPT violating parameter , that potentially signifies the breakdown of CPT operation as a result of quantum decoherence of matter in some models of quantum gravity, by means of identifying appropriate asymmetry parameters in the time evolution of intensities (14) between the two decays in a B factory, based on observables that have already been used in previous studies [8] probing independently T, CP and CPT symmetries in the absence of . In the current analysis we have included, simultaneously with the , also the conventional CPT parameter , already considered in [8], which parameterizes CPT violation in the case of a well-defined CPT operator which, however, does not commute with the Hamiltonian of the system, indicating a violation of CPT parameterized within the framework of effective field theories (e.g. due to Lorentz symmetry violation by a space-time background), in contrast to the parameter , which goes beyond that framework.

As we have demonstrated in the present article the set of observables of the B system (26), (27), (28) and (29) allow for a simultaneous determination (bounds) of the CPT violating parameters and , which can thus be disentangled. The results obtained from the experimental data from the BaBar measurements [2] (see Table 1(I)) are sensitive for the first time to , pointing towards a deviation from , which we interpret as an upper bound. The observables (28), (29) are also sensitive to , but they do not show any significant deviation from , and in this sense they are inferior to the previous analyses [21] using equal sign semileptonic decay asymmetries of the B system, which yield . The results (26) and (27) also allow a fit to the CPT violating parameter , and they are compatible with the previous determination in [8] and the one performed by the BaBar collaboration in [29], pointing towards a effect in , also interpreted as an upper bound for the corresponding parameter.

Moreover, the parameters that measure the presence of wrong flavour decays in , i.e. , and , do not show any significant deviation from zero and the results are consistent with [8]. In the case of and we observe that they differ by more than with respect to the determination in [8] without including the effect. Should this persist in the future, it could affect the precise determination of the unitarity triangle angle .

Before closing we stress once more that a quantum-gravity-decoherence-induced CPT violating and entanglement-weakening parameter may not only characterize the initial state of an entangled (neutral) meson system, but may also be generated as a result of a decoherening time evolution that goes beyond the local effective field theory framework [22]. A full analysis of that case will appear in a forthcoming publication.