Matter-antimatter asymmetry and non-inertial effects

We investigate non-inertial effects on CP-violating processes using a model, based on the framework of quantum field theory in curved spacetimes, devised to account for the decay of accelerated particles. We show that the CP violation parameter for the decay of accelerated kaons into two pions decreases very slightly as very high accelerations are achieved, implying decreased asymmetry between matter and antimatter in this regime. We discuss the relationship between these results and cosmological processes surrounding matter-antimatter asymmetry and argue that, due to the connection between non-inertial and thermal phenomena established by the Unruh effect, this kind of computation may prove useful in furthering the understanding of thermodynamical effects in curved spacetimes.


Introduction
One of the great modern scientific mysteries is the abundance of matter over antimatter. Cosmic matter-antimatter asymmetry arises in the framework of elementary particle physics through baryogenesis models, which offer mechanisms for obtaining matterantimatter asymmetry from an initially symmetric Universe [1]. It is well established that baryogenesis requires three crucial ingredients [2]: (i) baryon number violation, (ii) violation of C (particle-antiparticle) symmetry and the combination CP of C and P (left-right or parity) symmetries and (iii) departure from thermal equilibrium.
In the Standard Model (SM) the matter-antimatter asymmetry of the Universe is credited to CP violation (CP v), although experimental evidence has systematically shown that it may not be sufficient to explain this imbalance [3,4]. CP v has been observed in various weak decays involving strange and beauty quarks, being recently confirmed by the LHCb collaboration for the charmed D meson [5]. The main sources of CP v in the SM are (i) the quark sector, involving the Cabibbo-Kobayashi-Maskawa (CKM) matrix, (ii) the strong interaction and, (iii) the Pontecorvo-Maki-Nakagawa-Sakata (PMNS) matrix in the lepton sector [6]. The CKM matrix has been observed experimentally and can only account for a small portion of the CP v required to explain the matter-antimatter asymmetry. The failure, until now, to observe the electric dipole moment of the neutron in experiments suggests that any CP v in the strong sector is also too small to account for the necessary matter-antimatter asymmetry in the early Universe. In the case of neutrinos being Majorana fermions, for example, the PMNS matrix could have two additional CP -violating Majorana phases, which would lead to a new source of CP v in the framework of the SM in the lepton sector. Alternatively, CP v in the lepton sector could result, experimentally, to JHEP03(2021)285 be too small to account for matter-antimatter asymmetry, but additional sources of CP v could arise from some still unknown physics beyond the SM.
In general, CP -violating observables are computed from decay amplitudes, decay rates, or quantities derived from these, such as a particle's lifetime, which is obtained from the decay rates and is regarded as one of its inherent and characteristic properties. It is well known since the seminal work of Fulling, Davies and Unruh [7][8][9] in quantum field theory in curved spacetimes (QFTCS) that a uniformly accelerated detector moving through the usual flat spacetime vacuum of a conventional quantum field theory responds as though it were in a thermal bath of temperature where a is the acceleration (this is often referred to as the Unruh effect). One may, thus, expect that the acceleration would causes a modification of particle lifetimes with respect to their own proper time, i.e., in their accelerated rest frame, much as a thermal bath would do to a particle in an inertial frame. This is indeed the case -as was shown by Müller [10] and Vanzella and Matsas [11,12] in their investigation of the decay rates of accelerated particles -and the objective of this paper is to investigate the impact of these non-inertial effects on CP -violating observables, i.e., we are interested in the non-inertial factor of the following heuristic equation, Specifically, we investigate these effects on the decay rates and in a CP v parameter in neutral kaons. We may envisage a connection between the cosmological processes responsible for the existence of matter-antimatter asymmetry and non-inertial effects by noting the similarities between different effects arising from QFTCS, particularly those concerning the thermodynamical phenomena stemming from the existence of event horizons in certain spacetimes [13,14], along with the observation that the temperature of the Universe is tied to particle creation due to its expansion [15,16]. In what follows, we argue that investigations of mechanical phenomena such as the ones alluded above can play an important role in understanding thermodynamical effects, given this connection.
Throughout this work we use natural units 8πG = c = = k B = 1, unless stated otherwise, and the (−, +, +, +) convention for the metric signature.

CP violation in the K-system
We start by briefly reviewing the basic ideas related to CP v in the kaon system [17,18]. In order to understand the unusual properties of neutral kaons, we first observe that, since both K 0 and K 0 decay in two pions,

JHEP03(2021)285
they are not independent particles with respect to the weak interaction. Charge conjugation C and parity P operations on kaon and pion states result in From these equations we see that the CP operation does not leave the kaon states invariant, and K 0 → 2π decays appear to be prohibited for CP -symmetric processes.
If CP symmetry is to be restored, we may consider, alternatively, linear combinations of theses states resulting in CP eigenstates, Using eq. (2.2) one obtains According to eq. (2.2) and (2.4), K 1 (K 2 ) can only decay into a state with eigenvalue CP = +1 (CP = −1). Hence, we have only two possible decay channels: K 1 → 2π and K 2 → 3π. Furthermore, a CP projection operator can be defined, To accommodate CP v in weak processes, the following mass eigenstates are introduced: where the (complex) mixing parameter q characterizes the strength of the CP v. There is a measurable difference in the lifetimes (with τ S < τ L ) and masses (with m K S < m K L ) of these eigenstates, justifying the "S" (referring to the short-lived or small state) and "L" (referring to the long-lived or large state) labels. We note that which shows that the states in eq. (2.7) are mixtures of CP = +1 and CP = −1 eigenstates. If q = 1, the system is CP -symmetric and the projections in eq. (2.8) reduce to those of eq. (2.6). The lifetimes of the states in question are [6] τ (2.11) In figure 1 we see a schematic representation of the Cronin-Fitch experiment [19], the first experimental result to provide clear evidence, from kaon decays, that CP symmetry could be broken. An initially mixed beam of K S and K L is injected into a tube. The K S particles traverse only their characteristic distance cτ S (region 1), decaying rapidly to 2π. The beam that emerges in region 2 is composed solely of K L particles. Therefore, if CP is a conserved symmetry (i.e., if q = 1) one should only observe the decay to 3π in this region. CP v occurs because decays of type K L → 2π, are also observed, implying that the weak eigenstates do not correspond to the CP eigenstates.
The CP v that occurs in the transition K → 2π can be studied by defining an observable η that relates the transition amplitudes A of the K S and K L species: Then, putting η = η 00 = η +− and using eq. (2.8) one obtains It must be noted that the experimental data suggests a difference between the absolute values of η 00 and η +− [6]: (2.14) In fact, setting η +− = + and η 00 = − 2 , the value of is nonzero but three orders of magnitude smaller then the value of [18]. Alternatively, |η| 2 is given in terms of the decay rates Γ by which we shall use in the subsequent analysis.

JHEP03(2021)285 3 Acceleration and CP -violating decays
The Unruh effect (see [20] for a review) states that accelerated observers perceive the inertial vacuum as a thermal bath of particles. 1 It seems, then, that accelerated particles should have increased proper decay rates, as they would if immersed in a thermal bath. This is in fact the case, as was show for scalar fields by R. Müller in [10] and for spinorial fields by D.A.T. Vanzella and G.E.A. Matsas in [11,12]. These results also show a dependence of the rate of the increase on the mass of the decaying particle. Given the relationship between CP v and the kaon decays described above, we seek to investigate the impact that non-inertial effects have on CP v in the K-system with a model for the decay of accelerated kaons (related investigations, probing CP v in the lepton sector, can be found in [21][22][23]). Such a model, based on the ones introduced by Müller, is presented below and used to evaluate the behavior of the decay rates and of the CP v parameter η, the latter of which being possible due the difference in the masses of K S and K L .

Model for accelerated decays
We consider the following interaction Lagrangian, where Φ is a scalar field of mass M , φ 1 and φ 2 are scalar fields of mass m and G Γ is the coupling parameter of the interaction. Although kaons and pions are described by pseudoscalar fields, it is assumed that scalar fields provide approximate descriptions of their behavior. The decay rate may then be obtained from the decay amplitude, namely the transition amplitude for the process Φ → φ 1 φ 2 , given, up to first order in G Γ , by where the final state consists of two particles with momenta k 1 and k 2 . The decay probability can be computed from the amplitude, where f (x) is the mode associated to the initial state |i and 0|φ † j (x)φ j (x )|0 is the Wightman function of φ j .

JHEP03(2021)285
Choosing a frame in Minkowski spacetime where x = x and assuming t > t (with x = (t, x) and x = (t , x )), the Wightman function of a complex scalar field φ with mass m, for timelike separations of x and x , reads It is straightforward to show that the function above reduces to 1 (m∆s) ∆s , where H 1 is a Hankel function of the second kind and ∆s is the spacetime interval of the timelike separated events x and x .
Assuming that f is peaked over a trajectory x(τ ) parameterized by its proper time τ , i.e., over the trajectory of a particle, we can write it as (3.6) in the instantaneous rest frame. Assuming, furthermore, that the decay products do not deviate much from this trajectory, the decay probability can be written as where κ is given by κ = | d 3 x h(x)| 2 . A change of variables of the form v := τ − τ leads to Since the Wightman function depends only on the difference x − x , the integral over τ is trivial (and infinite). We make use, thus, of the proper decay rate i.e., the decay probability by unit of proper time.  Table 1. Values for the masses of the pions π 0 , π ± and the neutral kaon K 0 and the mass difference between K L and K S [6].
where a is the magnitude of the proper acceleration. The squared spacetime interval ∆s 2 at certain values τ and τ of the proper times is given by If we substitute ∆s 2 from (3.11) into (3.9), after introducing the variable u = a(τ − τ )/2 = av/2, we are led to , (3.12) which is the decay rate for uniformly accelerated scalar particles.

Results and analysis
The experimental data relevant to the computation of eq. the parameters c i , i = 1, . . . , 5, were obtained by fitting the curve to the numerical data. The form of Γ fit is based on the singular parts of the integral in eq. (3.12) (given in eqs. (A.19)). These results clearly indicate that the decay rate for this process grows as the acceleration increases, as expected from the results of previous analyses. It is also evident that this JHEP03(2021)285 There is also the problem of accelerating neutral particles, for which there is no obvious solution (see section 4 for a discussion on why a cosmological setting may be better suited for the observation of effects of this kind).
To investigate the behavior of η it is sufficient, having eq. (2.15) in mind, to compute Γ for the two values of M corresponding to the K S and K L masses and take their ratio. Since the difference between these values is very small, we take m K S ≈ m K 0 and m K L ≈ m K 0 + (m K L − m K S ). A small inconvenience for this computation is the fact that G Γ (K S → ππ) and G Γ (K L → ππ) must be different to account for the considerable difference in the magnitudes of Γ(K S → ππ) and Γ(K L → ππ). This may be remedied by considering a rescaled CP v parameter η , given by The numerical results for |η | 2 − 1, as well as a plot of an approximation of it obtained from the singular parts of eq. (3.12) (see appendix A for more details), can be found in figure 3. While they point to a variation of |η| 2 with increasing acceleration, which implies that CP v is sensitive to non-inertial effects, it is a decrease -which indicates that the contribution of these kaon decays to CP v is smaller at very high acceleration scales -of the order of 10 −13 |η| 2 (currently impossible to detect, even at extremely high accelerations, given that the experimental uncertainty for η is of at least 0.4%). The tiny fluctuations of the numerical data at small acceleration scales can be attributed to errors in the numerical calculation of the integral for the decay rate.  Further investigation of the behavior of η using the singular parts of the decay rate integrals leads to an interesting relationship between the amplitude of the variation of η and the value of the relative mass difference between the two species modeled by the field Φ.

JHEP03(2021)285
As an example of the kind of physics this would entail, consider, in place of K S and K L , two species Φ (1) and Φ (2) with masses (in units of m K 0 ) M 1 = 1 + (M 1 − M 2 ) > 1 and M 2 = 1, respectively. Using the singular integral approximation mentioned above, the computation of i are scalar fields of mass m = m π 0 /m K 0 , for different values of M 1 − M 2 yields the results show in figure 4. They indicate a roughly linear relation between the mass difference and the amplitude of the fluctuation of |η | 2 − 1, from which one may infer that the scale of the effect is given by the relative mass difference (M 1 − M 2 )/M 2 (or, in the case of the kaon decays, (m K L − m K S )/m K 0 ). It seems reasonable to infer that for a species with M 1 < M 2 there should be an increase in the magnitude of η, still proportional to the mass difference M 2 − M 1 . This is the case for K S → 3π processes, which are also CP -violating [6], though a model for decays of this kind would be slightly different. The nature of the results on the K S → 3π channel is, nevertheless, expected to be the same as for the K L → 2π channel (see [10] for a comparison of computations using the method discussed above for different models).

Cosmological effects
An intimate relationship between non-inertial and cosmological effects exists in QFTCS: the presence of an event horizon is correlated to the existence of non-Hilbert-Schmidt Bogoliubov transformations, i.e., QFT constructions that are not unitarily equivalent, which manifest themselves as correspondences between vacuum and thermal (KMS) states (see [25] JHEP03(2021)285 for an overview). This is the case for the aforementioned Unruh effect, the well known Hawking effect for black holes [13] and the Gibbons-Hawking effect in de Sitter spacetime [14], all of which present Killing horizons. A similar effect is that of particle creation in expanding universes, first described by Parker in [15,16], which predicts a thermal signature for particles created by the expansion of a Friedmann-Lemaître-Robertson-Walker (FLRW) universe. This section seeks to draw parallels between these effects and to encourage further research on the impact of phenomena stemming from QFTCS on other areas of physics by discussing the extension of the predictions of section 3 to cosmological settings.
A direct connection between the Unruh effect and the Gibbons-Hawking effect is established by noting that both are due to the presence of a Killing horizon. The Killing fields generating these horizons are the Lorentz boost generators for the Unruh effect and boostlike generators for the Gibbons-Hawking effect (in fact, these vector fields can be lifted to Lorentz boost generators on a 5-dimensional Minkowski spacetime where the de Sitter spacetime can be embedded). The surface gravity κ of these horizons is related to the temperature of the KMS state by with κ = a for the Unruh effect and κ = Λ/3 for the Gibbons-Hawking effect, where Λ is the cosmological constant. One may then wonder if CP -violating processes in de Sitter spacetime would behave like those of accelerated particles: would an observer also measure the difference in the behaviour of the observables discussed above when compared to measurements by an inertial observer in Minkowski spacetime? This seems reasonable (so long as a = Λ/3), given that both the accelerated observer in Minkowski spacetime and the observer in de Sitter spacetime perceive KMS states of the same temperature. A computation of decay rates of particles in a thermal bath is discussed in [26,27], as is the relationship between it and the results obtained in [11,12], which firmly establishes the complementarity of the non-inertial and thermal effects. With respect to particle creation in expanding universes (i.e., in FLRW universes), the connection is more tenuous: though there exists a kind of horizon in solutions in this class, it is not (in general) a Killing horizon but a particle horizon. Nevertheless, the prediction

JHEP03(2021)285
of QFTCS is that the expansion leads to particle creation at a temperature where S(t) is the scale factor and S 1 is some lower bound for the scale factor (see [16] for more details). This allows one to draw the same conclusion as above (so long as a ≈ √ 2πS 1 /S). More realistic cosmological models, such as the ΛCDM model, present both particle and event horizons (the latter being closely related to the de Sitter event horizon, given the presence of the cosmological constant), leading to a combination of effects, though the Gibbons-Hawking component is very small for our Universe, since Λ = 7.15(19) × 10 −121 [28] (compare to a 2 TPL ≈ 10 −60 ). In the context of the thermal history of the Universe, kaons freeze-out near the QCD crossover, at a temperature T QCD ≈ 155 MeV/k B [29], the moment at which the CP -violating processes discussed above are expected to contribute to the observed matterantimatter asymmetry. Inverting the expression for the Unruh temperature allows the determination of the acceleration scale corresponding to the crossover temperature, a QCD = 2πT QCD ≈ 2m K 0 . Figure 3 then allows us to infer that, at the crossover, η would be near to its lowest value, implying a lower amplitude of CP v when compared to current values and diminished matter-antimatter asymmetry.

Conclusions
The model introduced above for the K → 2π decays allows for the computation of the decay rate over accelerated trajectories, corroborating the results of ref. [10]. The subsequent analysis of the behaviour of the CP v parameter η leads to the conclusion that its squared magnitude |η| 2 decreases very slightly with increasing acceleration (around three parts in 10 14 for an acceleration a ≈ 2m K 0 ≈ 4 × 10 32 m/s 2 ). We have shown that the amplitude of the decrease is proportional to the mass difference between the weak eigenstates (for the kaons, (m K L − m K L )/m K 0 ≈ 10 −14 ).
The discussion on the relationship between the Unruh effect, the Gibbons-Hawking effect and the phenomenon of particle creation in accelerating universes led to the argument that, given the complementarity between non-inertial and thermal effects, a similar conclusion may hold for fields in de Sitter and FLRW universes (though the proper computations still need to be executed to fully justify this affirmation): CP v contributions of this nature are smaller at higher temperatures in these spacetimes.
Further research on the connection between the results presented in this work and those on the thermal dependence of CP v in the quark sector is in progress. A decrease in the effectiveness of CP -violating processes in this sector with increasing temperature has been shown in [30], but the relation between it and our result is not immediately clear.

A Numerical treatment of the decay rate
In this appendix, the numerical treatment of the decay rate, eq. (3.12), is presented. The integral appearing in this expression is 1 [2m sinh(u)/a] 2 sinh 2 (u) .
It is difficult to tackle this expression analytically, so the numerical approach is favored. Problems in the implementation of the numerical methods appear due to singularities in this expression, which means that the singular parts of the integral must be separated into an integral I 0 and treated analytically, so that I can be computed as where γ is the Euler-Mascheroni constant. Since this function appears squared in the expression for the decay rate and also involves a division by a second order polynomial on the argument of function, 2m sinh(u)/a, the following expression is computed: