Gravity, entanglement and CPT-symmetry violation in particle mixing

TL;DR

This paper investigates how self-gravitational interactions affect particle mixing, revealing CPT-symmetry violation linked to entanglement growth, with potential implications for early universe conditions and experimental simulations using Rydberg atoms.

Contribution

It introduces a novel mechanism where gravitational interactions induce CPT-symmetry violation through entanglement in particle mixing systems.

Findings

CPT-symmetry is broken due to T- and CP-symmetry dynamics.

Entanglement increases with system size, affecting particle oscillations.

Proposes experimental tests using Rydberg atom setups.

Abstract

We study the probability oscillations of mixed particles in the presence of self-gravitational interaction. We show a breaking of the CPT-symmetry due to the contemporary violation of the T-symmetry and the CP-symmetry preservation. This violation is directly associated to the rising of the entanglement among the elements of the system that can be seen as a pure many-body effect scaling with the number of the elements in the system. This effect could have played a relevant role in the first stages of the Universe or in core of very dense systems. Experiments based on Rydberg atoms confined in microtraps can simulate the mixing and the mutual interaction and could allow to test the mechanism here presented.

Equations

$$|n_A⟩$$

$$|n_B⟩$$

$$H^{(1)}=E+\frac{c^2}{2E}(m_1^2|m_1⟩⟨m_1|+m_2^2|m_2⟩⟨m_2|).$$

$$H^{(1)}={\omega }_0{\sigma }^z;{\omega }_0=\frac{c^2}{4E}(m_1^2-m_2^2).$$

$$P_{n_A\to n_B}=P_{n_B\to n_A}={\sin }^2(2\theta ){\sin }^2({\omega }_{0}t).$$

$$H^{(2)}$$

$$H^{(2)}=\omega ({\sigma }_i^z+{\sigma }_j^z)+\Omega {\sigma }_i^z\cdot {\sigma }_j^z,$$

$$\mathcal{P}({\rho }_i(t))=1-\frac{1}{2}{\sin }^4(2\theta ){\sin }^2(2t\Omega ).$$

$$P_{n_A\to n_B}$$

$$P_{n_B\to n_A}$$

$${\Delta }_T$$

$$H^{(N)}=\sum _i{\omega }_i{\sigma }_i+\frac{1}{2}\sum _{i,j}{\Omega }_{i,j}{\sigma }_i\cdot {\sigma }_j.$$

$${\rho }_k(t)=\frac{1}{2}\left(1+\sum _{\alpha }⟨{\psi }^{(N)}(0)|U^{†}(t){\sigma }_k^{\alpha }U(t)|{\psi }^{(N)}(0)⟩{\sigma }_k^{\alpha }\right),$$

$$|\psi (0)⟩=\sum _{\{l\}}{R}_l\left(a_k|m_1,l⟩+b_k|m_2,l⟩\right),$$

$$|\tilde{\psi }(t)⟩=\sum _{\{l\}}{R}_l\left(a_k{e}^{-ı(\omega +{\Gamma }_k)t}|1,l⟩+b_k{e}^{ı(\omega +{\Gamma }_k)t}|2,l⟩\right),$$

$$\!\rho_{k}(t)\!=\!\frac{1}{2}\!\left(\begin{array}[]{cc}1+\zeta_{k}\cos(2\theta)&\zeta_{k}e^{-\imath\phi}\sin(2\theta)a_{k}^{*}(t)\\ \zeta_{k}e^{\imath\phi}\sin(2\theta)a_{k}(t)&1-\zeta_{k}\cos(2\theta)\end{array}\right),\!\!\!\!\!\!\!$$

$$a_k(t)=e^{ı2{\omega }_{k}t}\prod _{j=1}^N(\cos (2{\Omega }_{k,j}t)+ı{\zeta }_k\cos (2\theta )\sin (2{\Omega }_{k,j}t)),$$

$$\mathcal{P}({\rho }_k(t))=1-{\sin }^2(2\theta )\left(1-|a_k(t){|}^2\right).$$

$$P_{n_A\to n_B}$$

$$P_{n_B\to n_A}$$

$$a_k(t)\simeq e^{ı2{\omega }_{k}t}\left(1\pm 2ı\cos (2\theta )\sum _{j=1}^N{\Omega }_{k,j}t\right),$$

$${\Delta }_T$$

$$H^{(N)}(t)=\sum _k{\omega }_k(t){\sigma }_k+\frac{1}{2}\sum _{k,j}{\Omega }_{k,j}(t){\sigma }_k\cdot {\sigma }_j.$$

$$U(t)=\exp \left(\sum _l\frac{1}{l!}{\Lambda }_l\right),$$

$${\Lambda }_1$$

$${\Lambda }_2$$

$$U(t)$$

$$\tilde{H}(t)$$

$${\tilde{\omega }}_k$$

$${\tilde{\Omega }}_{k,j}$$

Full text

I introduction

Particle mixing and oscillations have provided some of the most direct and robust indications of physics beyond the standard model Wu1957 ; Giunti2007 ; Griffiths2008 . We have several examples of such phenomenon both in the bosonic and fermionic sectors. In the first sector we have mixing among axion–photon Raffelt1996 ; Sikivie2011 ; Capolupo2019_1 , $\eta$–$\eta^{\prime}$ Pham2015 , neutral kaons Christenson1964 and $B$ meson Abashian2001 . In the second one we can find the neutrino flavor oscillations Bilenky1978 ; Gonzalez2008 , the neutron–antineutron oscillations that could be observed in the next generation of experiments using slow neutrons with kinetic energies of a few $meV$ Phillips2016 , and the quark mixing Nakamura2010 . Apart from the last one which involves particles confined inside hadrons, all the other mixing phenomena concern only neutral particles. All of them are characterized by the fact that the physical fields, called flavor fields, are superpositions of free fields with definite different masses.

Since the difference between the masses is very small, also weak perturbations can produce measurable deviations from vacuum oscillation frequencies. The extreme sensibility to a wide set of perturbations is at the basis of different experiment proposals. For instance, many studies in recent years have been devoted to the possibility to test the quantum nature of gravity using concepts of quantum information theory Bose2017 ; Marletto2017 . These proposals are based on the idea of using a system in which the intensity of the gravitational interaction depends on some internal degrees of freedom Giampaolo2018 . The sensibility of the oscillation frequency of neutral particles such as neutrino provides a natural system where to analyze these effects Marletto2018 . Indeed, it is well-known that neutrinos interact exclusively via gravity and weak interaction. This last interaction is stronger than gravity but has an extremely short range Greiner2009 (about $d=10^{-16}/10^{-18}m$) and, hence, it can be neglected for distances bigger than $d$ for which gravity survives.

The effects of gravity on the oscillation of neutral particles are not limited to a change in the frequency of flavor oscillations. Gravity is also considered as one of the possible sources of decoherence in flavor oscillation Ellis1984 ; Amelino-Camelia2000 ; Rovelli2004 that leads to many interesting effects like the $CPT$-symmetry violation in particle mixing Gago2001 ; Guzzo2016 ; Lisi2000 ; Benatti2000 ; Benatti2001 ; Capolupo2018 ; Capolupo2019 . In all these papers the non–unitary evolution is introduced by considering a dissipator that generates a completely positive dynamics Benatti1997 ; Benatti1998 . This dissipator summarizes the effects of all possible sources of decoherence and does not allow to analyze the origin and the relative weight of the different sources of decoherence.

In the present paper we adopt a different approach to analyze the role of gravity in the particle mixing phenomena. Instead to consider a single particle as an open system affected by several uncontrolled phenomena, we consider an ensemble of self-interacting particles as a closed system where all internal physical quantities are under control. Therefore, we consider a system of $N$ mixed neutral particles evolving under the self-gravity and neglect all other possible interactions acting between the system and the environment. We prove that, because of the difference in mass of the free fields, the self-gravity induces a violation of the $T$-symmetry whereas the $CP$-symmetry is preserved. This implies a violation of the $CPT$-symmetry. Such a violation represents an emerging many-body effect associated to the rising of a non-vanishing entanglement among the different fields.

Both the entanglement and the $CPT$ violation are extremely small and difficult to be observed in laboratories. However, being many-body effects, they are related to the number $N$ of elements of the system and to the density $n$ of the particles. Therefore, the $CPT$ violation induced by the gravity could play an important role in very dense astrophysical objects and it could have affected the early stages of the Universe Hamada2016 .

The paper is organized as follows. In Sec. II, for the reader’s convenience, we resume the main results of the quantum mechanical approach to the particle mixing. In Sec. III we consider the effects of gravitational interaction in a simple model made of two interacting fields, and show the $CPT$-symmetry breaking and the role played by the entanglement. In Sec. IV we generalize the formalism to the case of many interacting particles. In Sec. V we discuss how to generalize our results to the case in which the distance among the particles changes in time and in Sec. VI we draw our conclusions.

II Neutral particle oscillation in vacuum

For sake of completeness, in this section we review some aspects of neutral particle oscillations. We consider the very general case represented by the mixing of two flavor fields named $n_{A}$ and $n_{B}$. To fix the ideas, in case of neutrino oscillations, the two flavor fields coincide with the leptonic flavors as $n_{A}=\nu_{e}$ and $n_{B}=\nu_{\mu}$ while, in the case of neutron–antineutron oscillations we have $n_{A}=n$ and $n_{B}=\bar{n}$. Since the flavor fields do not coincide with those of definite masses, the mixing relations can be written as

[TABLE]

where $\theta$ is the mixing angle and $\phi$ is the Majorana phase which is zero in case of Dirac fermions Giunti2007 and $\ket{m_{i}}$ are the states with definite masses $m_{i}$.

Neglecting any interaction with the rest of the Universe, the particle is a closed system which travels through space with its energy $E$. Assuming that the masses $m_{i}$ are much smaller than $E$, we can write the Hamiltonian of mixed fields as

[TABLE]

Introducing the Pauli operator $\sigma^{z}$ which discriminates between the mass eigenstates, $\sigma^{z}\!=\!\ket{m_{1}}\!\!\bra{m_{1}}\!-\!\ket{m_{2}}\!\!\bra{m_{2}}$, and neglecting state-independent terms proportional to the identity operator in the Hamiltonian, eq. (2) becomes

[TABLE]

The single-particle state, that at $t=0$ is in one of the two flavor states of eq. (II), for $t>0$ evolves under the effect of $H^{(1)}$. Due to this evolution we have, for $t>0$, a non-vanishing probability to observe a change in the flavor state of the particle. The time-dependent expression of the flavor oscillation formula is than given by $P_{n_{A}\rightarrow n_{B}}=|\bra{n_{B}}\exp(-\imath Ht)\ket{n_{A}}|^{2}$ and is invariant under the exchange of flavors states, i.e. $P_{n_{A}\rightarrow n_{B}}=P_{n_{B}\rightarrow n_{A}}$. Explicitly, we have

[TABLE]

This is the well-known Pontecorvo formula Bilenky1978 that describes the oscillation of a neutral particle in the vacuum.

III Oscillation of two interacting neutral particles

Let us now consider the case in which the system that we analyze is composed not of a single particle but of two mixing particles interacting gravitationally.

Before starting our analysis, let us discuss the basic hypothesis we will use along the section. At first, we assume the validity of the equivalence principle between inertial and gravitational mass. Moreover, we represent the gravitational interaction with the standard Newtonian potential. Replacing it with ghost-free theories of gravity as the one showed in Ref. Biswas2012 would induce quantitative but not qualitative changes in the physical behavior. Furthermore, for the sake of simplicity, we assume that the particles travel in space with the same energy along the same direction, hence, keeping their relative distance, that we denote with $d$, a time-independent parameter. It is worth to note that this last assumption is made only to simplify the explanation of our results. Indeed, as we will show in Sec. V, it is possible to extend our results to the case of time-dependent $d$.

Within the above assumptions the Hamiltonian of the system made of two mixed particles that interact gravitationally becomes

[TABLE]

where $G$ is the gravitational constant, $d$ the distance between the two particles and the indices in the two particles states refer, respectively, to the $i$-th particle (the first) and to the $j$-th particle (the second).

It is useful to rewrite the Hamiltonian in eq. (5) in a more compact form. By recalling the definition of $\sigma^{z}_{i}$ and neglecting all terms proportional to the identity, the Hamiltonian of the system can be written as

[TABLE]

where $\omega=\omega_{0}+g(m_{1}^{2}-m_{2}^{2})$, $\Omega=g(m_{1}-m_{2})^{2}$, and $g=-\frac{G}{4d}$. Comparing the above Hamiltonian with the one in eq. (3) we can see that the presence of the gravitational interaction has two different effects. The first is that the interaction changes the value of $\omega$ from $\omega_{0}$ to $\omega_{0}+g(m_{1}^{2}-m_{2}^{2})$ while the second one is the appearance of a new term, with amplitude equal to $\Omega=g(m_{1}-m_{2})^{2}$, involving operators defined on both the fields. It is worth to note that both of them will disappear in the case in which $m_{1}$ and $m_{2}$ coincide.

Before going further, let us recall two basic results of quantum information that can be found in all quantum information books as, for example, that in Ref. Nielsen2011 . The first result is that, given a bipartite system, the evolution induced by Hamiltonian terms acting only on a single part of the system can never affect the value of the entanglement between the two parts. On the contrary, terms that act simultaneously on both the parts usually modify the entanglement. The second result is that when we project a pure state defined on a bipartite system into one of its parts, the projection obtained in such a way is still pure if and only if the state was not entangled.

In our case, the natural bipartition of the system under analysis is obtained by considering each part coinciding with one of the two particles. With respect to this partition the terms proportional to $\omega$ are local, since each one of them acts on one single particle. Therefore, they cannot create, or destroy, entanglement inside the system. On the contrary, the term proportional to $\Omega$, induced by the presence of the state-dependent gravitational interaction, is non-local respect to the natural bipartition. Hence, it can increase, or decrease, the value of the entanglement Wang2002 .

At $t\!=\!0$, i.e. when the two particles were created, we can assume that there is no entanglement between them. Thus, the initial state of the whole system is a two–body fully separable (i.e. without any entanglement between the two particles) pure flavor state of the form $\ket{\psi(0)}=\ket{n_{\eta}}_{1}\ket{n_{\chi}}_{2}$, where $\eta$ and $\chi$ could assume all possible combinations of $A$ and $B$.

Once fixed $\ket{\psi(0)}$, the state at $t>0$ can be obtained as $\ket{\psi(t)}=U(t)\ket{\psi(0)}$, where the time evolution operator is $U(t)=\exp(-\imath tH^{(2)})$. The operator $U(t)$ is unitary because we assume that the system under analysis is closed, i.e. does not interact with the surrounding world. Therefore, for any time $t\geq 0$, the state $\ket{\psi(t)}$ is still a pure state exactly as at $t=0$. But, if $\Omega\neq 0$ the state $\ket{\psi(t)}$ holds, in general, a non-vanishing entanglement between the two particles. This implies that the projection of $\ket{\psi(t)}$ on any of the two particles would be, in general, a mixed state.

It is possible to quantify how much a projection is pure using a quantity called purity defined as $\mathcal{P}(\rho_{i}(t))=\mathrm{Tr}(\rho^{2}_{i}(t))$, where $\rho_{i}(t)=\mathrm{Tr}_{j}(\ket{\psi(t)}\bra{\psi(t)})$ is the projection over the $i$-th particle of the state $\ket{\psi(t)}$ Nielsen2011 . The purity holds a relevant role in the theory of entanglement for pure states defined in bipartite systems. Indeed, it is also associated to the $2$–Renyi entropy, defined as $S_{2}=-\ln(\mathcal{P}(\rho_{i}(t)))$, that represents a proper measure of the entanglement between a particle and the rest of the system Horodecki2002 ; Plenio2007 ; Giampaolo2013 . Projection of fully separable states holds a purity equal to 1 and a vanishing Renyi entropy while entangled states are characterized by $\mathcal{P}(\rho_{i}(t))<1$ and a non-vanishing entropy.

In our case it is easy to verify that the time-dependent expression of purity is, independently on the initial state, equal to

[TABLE]

In the presence of flavor mixing, i.e. for $\theta\neq n\frac{\pi}{2}$, and for any time $t\neq n\frac{\pi}{2\Omega}$ (with $n$ integer), we have $\mathcal{P}(\rho_{i}(t))<1$ and, thus, the single particle state is not a pure state, hence, implying that $\ket{\psi(t)}$ is entangled. The expression of purity is also a proof of the fact that the entanglement between two flavor fields is a direct consequence of the presence of the gravitational interaction and vanishes if $\Omega$ is neglected, i.e. if $m_{1}=m_{2}$. This fact can be also considered a proof of the quantum nature of gravity since quantum correlations and entanglement, can be created only through a quantum channel Bennet1999 .

Let us now show that, in the system we are analyzing, the presence of the entanglement induces also a violation of $T$- and $CPT$-symmetry. In order to provide this proof, we consider a simple conceptual experiment. We take into account two copies of the system already described. We assume that the two copies are identical except for the fact that, in the initial state of the first, both particles are in the state $\ket{n_{A}}$, i.e. $\ket{\psi_{0}}=\ket{n_{A}}\ket{n_{A}}$ while in the second one are both in $\ket{n_{B}}$ and, hence, $\ket{\psi_{0}}=\ket{n_{B}}\ket{n_{B}}$. At the same time $t>0$, we observe, in both copies, one of the two particles and analyze the oscillation probabilities.

Let us name $\rho_{A}(t)$ the projection obtained in the first copy of our system and $\rho_{B}(t)$ the one obtained in the second copy. In these two cases, the probability of flavor transitions are given, respectively, by $P_{n_{A}\rightarrow n_{B}}=\mathrm{Tr}(\rho_{A}(t)\ket{n_{B}}\!\!\bra{n_{B}})$ and $P_{n_{B}\rightarrow n_{A}}=\mathrm{Tr}(\rho_{B}(t)\ket{n_{A}}\!\!\bra{n_{A}})$. Explicitly, we obtain

[TABLE]

It is easy to note that the probabilities in eqs. (III) are independent on the Majorana phase $\phi$, and, hence, the $CP$-symmetry is preserved, $\Delta_{CP}=0$. On the contrary, since the probability is not invariant under the exchange of the two flavors, we have a violation of the time-reversal symmetry

[TABLE]

For $m_{1}\neq m_{2}$ we have $\Delta_{T}=0$ only if $t=\frac{k\pi}{2\Omega}$, or $\theta=\frac{k\pi}{4},\,k\in\mathbb{Z}$. Since $\Delta_{T}\neq\Delta_{CP}$ we also have the violation of the $CPT$-symmetry. Therefore, the entanglement between the two particles induces a $CPT$-symmetry breaking.

It is worth to note that, even if in the system that we have analyzed, the presence of the entanglement induces a violation of the $T$–symmetry, this is not a general result. Indeed, it is possible to find several Hamiltonians that can induce entanglement without breaking the $T$-symmetry, or breaking both the $T$- and the $CP$-symmetry etc. To provide a simple example, if we consider an interaction that can be summarized by a Hamiltonian as $H=\Omega\,\sigma^{z}_{i}\cdot\sigma^{z}_{j}$, that can be obtained from the one in eq. (6) assuming $\omega=0$, we can immediately recover from eq. (7) and eq. (9) that $\Delta_{T}=0$ even in the presence of a non-vanishing entanglement.

A simple numerical analysis of the above model shows that, for many mixed particle systems, the non-unitary evolution effect is negligible. However, as we will show in the next section, this effect is a many-body effect and, hence, its relevance increases proportionally to the number of particles in the system.

IV Oscillation of $N$ interacting neutral particles

We now generalize the scheme presented above to the case where each copy of the analyzed system is made of a large number $N$ of particles. We consider the same assumptions used in the previous section: 1) we assume that the system is closed; 2) we take into account only the gravitational interaction among the particles; 3) we assume the identity between inertial and gravitational masses; 4) we consider the Newtonian potential valid; 5) we assume the invariance of the relative distances among the fields during the time evolution.

Within the above hypothesis, the system evolves with time-independent Hamiltonian which generalize that in eq. (6), i.e.

[TABLE]

The main difference between eq. (6) and eq. (10) is that, now, all parameters of the Hamiltonian depend on the index running on the set of particles. Indeed, $\omega_{i}$ and $\Omega_{i,j}$ are now given by $\omega_{i}=\omega_{0}+\sum_{j}g_{i,j}(m_{1}^{2}-m_{2}^{2})$ and $\Omega_{i,j}=g_{i,j}(m_{1}-m_{2})^{2}$, where $g_{i,j}=\frac{G}{4d_{i,j}}$ and $d_{i,j}$ is the relative distance between the $i$-th and the $j$-th fields. Despite this loss of symmetry, the Hamiltonian in eq. (10) still holds the fundamental property that it can be seen as a sum of mutual commuting terms. This property plays a key role in the rest of our paper. Indeed, usually the dynamic of a quantum many–body system is extremely complex to be analyzed exactly and numerical and/or approximate approach are needed. However this is not the case. In fact, exploiting such a property, we have that the time evolution operator can be written as the product of several operators each one of them accounts for the evolution induced by a single term of the Hamiltonian in eq. (10). Hence, collecting all the terms it is possible to obtain an exact expression of the state at a time $t>0$ and, more important for our analysis, of its projection into the Hilbert space defined on a single particle.

As in the previous section we consider, at time $t\!=\!0$, that the system is described by a fully separable state. We assume that the first $M$ particles are created in the state $\ket{n_{A}}$ and the rest is in the state $\ket{n_{B}}$, so that the initial state is $\ket{\psi^{(N)}(0)}=\bigotimes_{\alpha=1}^{M}\ket{n_{A}}_{\alpha}\bigotimes_{\beta=M+1}^{N}\ket{n_{B}}_{\beta}$. Soon after $t\!=\!0$ the system will start to evolve under the influence of the self-gravity and for any $t>0$, the whole system is represented by the pure state (because we are assuming that the system is closed) $\ket{\psi^{(N)}(t)}=U(t)\ket{\psi^{(N)}(0)}$, where the unitary time evolution operator is $U(t)=\exp(-\imath tH^{(N)})$. Knowing the initial state, the reduced density matrix on the selected $k$-th particle can be obtained in terms of the Pauli matrix Osborne2002 as

[TABLE]

where $\alpha$ runs over the ensemble $\{x,y,z\}$. Since all terms in the Hamiltonian commutes with each other, the operator $U(t)$ can be arranged as the product of three different terms, i.e. $U(t)=u_{k}(t)u_{k,r}(t)u_{r}(t)$. Here $u_{k}(t)=\exp(-\imath\omega_{k}\sigma^{z}_{k}t)$ is the part of the unitary evolution that acts only on the selected $k$-th particle, $u_{k,r}(t)=\exp(-\imath t\sum_{j}\Omega_{k,j}\sigma_{k}\cdot\sigma_{j})$ while $u_{r}(t)$ includes all the other Hamiltonian terms that do not involve directly the $k$-th field.

Taking into account the fact that Pauli operators on different particles commute with each other, we have that in the evaluation of $\rho_{k}(t)$ the operator $u_{r}(t)$ can be neglected. Hence, $\rho_{k}(t)$ depends only on $u_{k}(t)$ and $u_{k,r}(t)$. Moreover, since both $u_{k}(t)$ and $u_{k,r}(t)$ depend only on $\sigma^{z}_{k}$, we have that the coefficient of $\sigma_{k}^{z}$, and, hence, the elements on the diagonal of the reduced density matrix are time-independent. On the contrary, the coefficients of $\sigma_{k}^{x}$ and $\sigma_{k}^{y}$ depend on time and their derivation is long but straightforward. Substituting the expression of the flavor fields in eq. (II) in $\ket{\psi(0)}$, we have that the initial state can be written as

[TABLE]

where $\ket{m_{1},l}$ $(\ket{m_{2},l})$ is a generic tensor product of mass states in which the state in the $k$-th field is equal to $m_{1}$ ($m_{2}$). For the different parameters, $a_{k}$ ($b_{k}$) is equal to $\cos(\theta)$ $(e^{\imath\phi}\sin(\theta))$ for $k\leq M$ and to $e^{-\imath\phi}\sin(\theta)$ $(\cos(\theta))$ for $k>M$. On the other hand, $R_{l}=\prod_{s}c_{l,s}$, where $c_{l,s}$ is equal to $a_{l}$ ($b_{l}$) if in $\ket{l}$ the $s$-th field in the mass state $m_{1}$ ($m_{2}$).

From this expression, it is immediate to obtain the expression of $\ket{\tilde{\psi}(t)}=u_{k}(t)u_{k,r}(t)\ket{\psi(0)}$,

[TABLE]

where $\Gamma_{k}=\sum_{s}(-1)^{\lambda_{s}}\Omega_{s,k}$ with $\lambda_{s}=1$ ($\lambda_{s}=2$) if in $\ket{l}$ the mass state of the $s$-th particle is $m_{1}$ $(m_{2})$.

The knowledge of $\ket{\tilde{\psi}(t)}$ allows us to construct the reduced density matrix taking into account that $\bra{\psi^{(N)}(0)}U^{\dagger}(t)\sigma^{\alpha}_{k}U(t)\ket{\psi^{(N)}(0)}=\bra{\tilde{\psi}(t)}\sigma^{\alpha}_{k}\ket{\tilde{\psi}(t)}$. After some algebras we obtain the following general exact expression for the reduced density matrix,

[TABLE]

where $\zeta_{k}$ is a function that is equal to $+1$ for $k\leq M$, and to $-1$ for $k>M$, $a_{k}(t)$ is given by

[TABLE]

and we assume, as definition, that $\Omega_{k,k}=0$. It is worth to note that the expression of the reduced density matrix in eq. (14) is exact and obtained, once given the set of the relative distances, without any approximation and without the necessity to use any master equation approach.

As we have already said, since we are neglecting any interaction among the elements of the system and the surrounding world, the time evolution is unitary. As a consequence, $\ket{\psi^{(N)}(t)}$ is always a pure state. Therefore, it is possible to use the $2$–Renyi entropy defined as $S_{2}=-\ln(\mathcal{P}(\rho_{i}(t)))$ to quantify the total entanglement that any single particle shares with the rest of the system. From eq. (14), we obtain for the time-dependent purity

[TABLE]

Now the couplings $\Omega_{k,j}$ are not invariant under the change of fields. Thus, we have $|a_{k}(t)|^{2}<1,\,\forall t>0$, which reduces to 1 only at $t=0$. Therefore, for $t>0$, any single particle is entangled with the rest of the system.

By means of eq. (14) we can generalize the result presented in eqs. (III) and eq. (9) by analyzing the oscillation probability in two copies of the system in which the first one has $M=N$ and the second one has $M=0$. Differently from the previous case, now the reduced density matrices at $t>0$, and, hence, also the oscillation probabilities, are site-dependent. Thus, we have to consider the average over all elements of the system. Explicitly we obtain

[TABLE]

where $\mathrm{Re}(a_{k}^{(A)}(t))$ ($\mathrm{Re}(a_{k}^{(B)}(t))$) is the real part of $a_{k}^{(A)}(t)$ ($a_{k}^{(B)}(t)$) that are the functions $a_{k}(t)$ when $M=N$ ($M=0$). As well as in eq. (III), also the transition probability in eq. (IV) does not depend on the $CP$-violating Majorana phase, so that $\Delta_{CP}=0$. On the other hand, $P_{n_{A}\rightarrow n_{B}}\neq P_{n_{B}\rightarrow n_{A}}$ because of $a_{k}^{(A)}(t)\neq a_{k}^{(B)}(t)$ and, hence, $\Delta_{T}\neq 0$. In order to make this violation more evident, let us assume that $\Omega_{k,j}t\ll 1$. In this case $a_{k}(t)$ becomes

[TABLE]

where the sign $+$ is for the system composed at $t=0$ only of $n_{A}$-particles, i.e. $a_{k}^{(A)}(t)$, and the sign $-$ is for systems composed at $t=0$ only of $n_{B}$-particles, i.e. $a_{k}^{(B)}(t)$. Substituting eq. (18) in eqs. (IV), we have the time-reversal symmetry violation becomes

[TABLE]

Since $\Delta_{CP}\neq\Delta_{T}$, the $CPT$-symmetry is broken.

The exact value of the violation of the time-reversal symmetry in eq. (19)depends on the whole set of relative distances among the particles in the system. For system with large $N$ this set is not known but we can express $\Delta_{T}$ in terms of average values. Indeed, since the gravity has a very long range, then the sum in eq. (19) contains $N(N-1)$ non–zero terms (we assumed $\Omega_{i,i}=0\;\forall i$). Moreover, all these terms have the same sign. In fact, since gravity is attractive, all $\Omega_{k,j}$ are negative regardless of the particular choice of $k$ and $j$. On the other hand, the sign of $\omega_{k}$ depends on the difference between $m_{1}$ and $m_{2}$ and, hence, the sign does not depend on $k$. Therefore, inside the sum, for time short enough such that $\max(\omega_{k}t)<\pi/4$, all the terms have the same sign. In other words, by defining $f_{k}=\frac{\sin(2\omega_{k}t)}{N}\sum_{j}\Omega_{k,j}$ we have that $f_{k}$ is of the order of unity and all $f_{k}\,\forall k$ have the same sign. Therefore, defining $F$ as the average of $f_{k}$, i.e. $F=\frac{1}{N}\sum_{k=1}^{N}f_{k}$ we have

[TABLE]

where we see explicitly that $\Delta_{T}$ is proportional to the number of particles of the system.

Similar $CPT$ violation can be obtained for all configurations in which the difference $M$ and $N-M$ is of the same order of magnitude of $N$. Indeed, when this does not happen, as in the case in which at $t=0$ we have $N/2$ particles in the flavor state $\ket{n_{A}}$ and $N/2$ in $\ket{n_{B}}$, it is possible to show, using the Lindeberg–Lévy theorem VanderVaart1998 , that $\Delta_{T}$ is proportional to $\sqrt{N}$.

V Time-dependent relative distances

All the results obtained up to now were derived assuming constant the relative distances among the particles of the system. However, this assumption is not crucial. In this section we extend our analysis also to the more realistic case in which the distances change in time.

Removing the constraints of the independence of the distances on time, for any $t>0$, the Hamiltonian of the gravitational self-interacting system can be written as

[TABLE]

Differently from eq. (10), now the parameters $\omega_{k}(t)$ and $\Omega_{k,j}(t)$ depend explicitly on time. Indeed, they are, respectively, $\omega_{k}(t)=\omega_{0}+\sum_{j}g_{k,j}(t)(m_{1}^{2}-m_{2}^{2})$ and $\Omega_{k,j}(t)=g_{k,j}(t)(m_{1}-m_{2})^{2}$, where $g_{k,j}(t)=-\frac{G}{4d_{k,j}(t)}$ and $d_{k,j}(t)$ is the relative time-dependent distance between the $k$-th and the $j$-th fields.

The time dependence of the Hamiltonian affects the evaluation of the time evolution unitary operator. Indeed, this operator, that is obtained as solution of the Schrödinger equation can be, in general, written in terms of Magnus expansion Blanes2009 ; Blanes2010 as

[TABLE]

where the first terms of the expansion are

[TABLE]

where the square brackets denote the commutator between the Hamiltonian at diffferent times. Moreover, all the other terms of the Magnus expansion depend on a combination of commutators between the Hamiltonian in eq. (21) at different times Blanes2010 .

However, in the case that we are analyzing, it is easy to check that $[H(\tau_{1}),H(\tau_{2})]=0$ $\forall\tau_{1},\,\tau_{2}$ and, hence, we have that $\Lambda_{l}=0\,\forall\,l\geq 2$. Therefore, from eq. (22) we obtain that

[TABLE]

where

[TABLE]

The time evolution operator in eq. (24) is, formally, equivalent to that obtained in Sec. IV. Therefore, independently on $N$ and $M$, we can use the relations in eq. (V) to generalize the results obtained in Sec. IV to the case of time-dependent relative distances.

Before to conclude, it is worth to underline that this surprising result holds only because $[H(\tau_{1}),H(\tau_{2})]=0$ $\forall\tau_{1},\,\tau_{2}$. In the general case this is not true and, hence, the generalization to time-dependent relative distances cannot be evaluated exactly. In these cases we are forced to use different approaches, such as master equations, Dyson series expansions, etc.

VI Conclusions

We have shown that the gravity in a self-interacting particles mixing system leads to the $CPT$ violation. This violation is related to the emergence of a non-zero entanglement among the elements of the system induced by a difference of mass of the free fields. Moreover, since the gravitational interaction is additive, the $CPT$ violation is proportional to the number of elements of the system and its density. Therefore, this effect could play a crucial role in galactic objects and in the first stage of the Universe where the densities and the number of particles are very high.

The $CPT$-symmetry violation presented in this paper is not the first one discovered in the context of neutral particle oscillations. In neutrino physics, several studies have been devoted to the analysis of symmetry violations induced by dissipative dynamics Benatti2000 ; Benatti2001 ; Capolupo2018 . However, our work presents several aspects of novelty. In fact, instead to consider a single particle as an open system affected by several uncontrolled phenomena, we consider an ensemble of self-interacting particles as a closed system with all the physical quantities being under control. As a consequence we have a difference $CPT$ violation. In fact, in the previous works, $CPT$ violation was generated by a $CP$-symmetry breaking and not, as in our case, by a violation of the $T$-symmetry.

However, our results must not be considered in contrast with the ones presented in Refs. Benatti2000 ; Benatti2001 ; Capolupo2018 . Indeed, the non-unitary dynamics includes a wide family of physical sources of decoherence. On the contrary, we have limited our analysis to the effects due to the self-gravitational interaction so neglecting all other possible sources of decoherence Kostelecky1995 ; Kostelecky2004 ; Colladay1997 ; Amelino-Camelia2010 ; Bernabeu2017 . Nevertheless, our work paves the way to several other works, in which a detailed analysis of each individual contribution to decoherence can be realized.

Moreover, it is worth to note that the mechanism here presented is not only limited to gravitational interaction. In fact, the two main requirements are: 1) the presence of neutral particles whose flavor states are superpositions of the eigenstates of a free field Hamiltonian; 2) The presence of an interaction depending on the eigenstates of the free Hamiltonian. Within this hypothesis the interaction, not necessarily of gravitational origin, among two or more of these particles will generate entanglement and, hence, induce a $CPT$-symmetry violation in the flavors oscillations.

Table-top experiments, based on Rydberg atoms confined in microtraps and optically manipulated Pillet2016 , can simulate the mixing and the mutual interaction. In this system, the two internal states, i.e. the ground state and the excited Rydberg level can represent the mass eigenstates, whereas two particular orthonormal superpositions can simulate the flavor states, and the dipole-dipole interaction can play the role of the gravity Jau2016 ; Labuhn2016 ; Ravets2014 . Thus, next experiments on atomic physics could allow to test the fundamental laws and symmetries of nature.

Acknowledgements

K.S. is grateful to the Austrian Science Fund (FWF-P30821, FWF-P26783). A.C. thanks partial financial support from MIUR and INFN and the COST Action CA1511 Cosmology and Astrophysics Network for Theoretical Advances and Training Actions (CANTATA). SMG acknowledges support from the European Regional Development Fund for the Competitiveness and Cohesion Operational Programme (KK.01.1.1.06–RBI TWIN SIN) and from the Croatian Science Fund Project No. IP-2016–6–3347. SMG also acknowledges the QuantiXLie Center of Excellence, a project co–financed by the Croatian Government and European Union through the European Regional Development Fund–the Competitiveness and Cohesion Operational Programme (Grant KK.01.1.1.01.0004).