Phenomenological Effects of CPT and Lorentz Invariance Violation in Particle and Astroparticle Physics

It is well known that a fundamental theorem of Quantum Field Theory (QFT) set in at spacetime ensures the CPT invariance of the theory. This symmetry is strictly connected to the Lorentz covariance, and consequently to the fundamental structure of spacetime. Therefore it may be interesting to investigate the possibility of departure from this fundamental symmetry, since it can furnish a window to observe possible effects of a more fundamental quantum gravity theory in a"lower energy limit". Moreover, in the past, the inquiry of symmetry violations provided a starting point for new physics discoveries. A useful physical framework for this kind of search is provided by astroparticle physics, thanks to the high energy involved and to the long path travelled by particles accelerated by an astrophysical object and then revealed on Earth. Astrophysical messengers are therefore very important probes for investigating this sector, involving high energy photons, charged particles, and neutrinos of cosmic origin. In addition, one can also study artificial neutrino beams, investigated at accelerator experiments. Here we discuss the state of art for all these topics and some interesting new proposals, both from a theoretical and phenomenological point of view.

invariance under CP T symmetry, the combined transformation given by the product of C, P , and T in any order, is guaranteed in flat spacetime by a fundamental result: The so-called CP T theorem that was implicit in ' 51 Schwinger's work [1] and was more formally proved by Lúders [2] and by Pauli [3] a few years later. For every symmetry, it is interesting to test at which level of accuracy we can verify it and which are the conditions required to prove formally its validity.
Moreover, symmetries can even guide to the discovery of new physics, through the study of the consequences of their potential violation, or, at least, modification. In the specific case of the Lorentz and CPT invariance the search for modifications of these fundamental principles is motivated by the idea that small violations of space-time symmetries could provide an experimental window on the Planck scale effects caused by a more fundamental quantum gravity theory.
In this work we will give an introductory demonstration of the CP T theorem, showing that this result is fundamental since it is strictly related to LI and therefore to the space-time intrinsic structure. We will start introducing the algebraic demonstration, with a review of the original Jost proof [4], based on axiomatic field theory formulation. Then we will consider the proof based on the Lagrangian field theory, following by Bell [5], Lúders [2], Pauli [3], and even Schwinger [1,6]. After the introduction of the CP T theorem, we will discuss the connection between Lorentz and CP T invariance in QFT, giving a demonstration of the Greenberg theorem, which states that CP T violation implies Lorentz Invariance Violation (LIV) [7]. Then we discuss the implications of CP T symmetry and expose the motivations for conducing the search for possible violations of this fundamental symmetry. In the following section, we explain how it may be possible to generalize the CP T theorem to curved spacetime and the implications of its validity or violation in the gravity sector.
Finally we discuss the phenomenology that can be induced by CP T violation and by LIV in the astrophysical sector, in particular in gamma-ray bursts, Cosmic Rays (CRs), and the neutrino sector.

CP T Theorem
In this section we will give a review of the CP T theorem by first introducing the axiomatic demonstration and then the Lagrangian or operator-based formulation, as discussed in [8,9,10]. The CP T theorem states that a Lorentz invariant QFT that preserves unitarity and locality and is defined in a flat spacetime is invariant under the action of the operator θ = CP T , where the product order of C, P , and T does not matter. Under the action of such an operator, a generic physical state transforms as: so the theory must describe particles and antiparticles and: Since two combined charge conjugations, parity, and time reversal transforms leave the final state unaffected, it follows that the unitary property: Now it is important to underline that T is antilinear, while C and P are both linear transformations, so the total transformation θ = CP T is antilinear, therefore θ results antiunitary. Hence for a generic hermitian operator A, the CP T correlated one is given by: therefore: thanks to the antiunitarity of θ.

Wightman Axioms
The axiomatic argumentation is important since it can give an inside view on the importance of CP T symmetry caused by its strict correlation with Lorentz invariance and can underline that complete CP T symmetry presents a more fundamental nature than either C, P , or T alone. Before introducing the Jost argumentation [4], it is first necessary to review the Wightman axioms background required to define a QFT on a flat spacetime, using for instance the version present in [11]: 1. Poincaré covariance of the Hilbert space where the theory is set.
This means that unitary operators U (Λ, a) exist that implement Lorentz transformations and space-time translations; 2. Existence of a vacuum state |0 .
There is a unique state, the vacuum, that is unaffected by Poincaré transformations up to a phase: |U (Λ, a)|0 = e iφ |0 . This implies that this state can only have null four-momentum and angular momentum, since these quantities change under Lorentz transformations. The vacuum state must have even the lowest allowed energy and must be cyclic, which is acting on it via the creation operators every Hilbert space state can be constructed; 3. Fields constructed via operators.
All physical quantities can be constructed using polynomials of fields acting on the Hilbert space. The fields transform under Poincaré symmetry as scalars, spinors, and tensors. The fields are defined in such a way that each one corresponds to a definite physical state, i.e., a particle with defined physical quantities and quantum numbers, such as mass, spin etc.; 4. Energy positivity.
The Hamiltonian operator is supposed to have non negative eigenvalues.
This property together with covariance under the action of Lorentz group implies that the physical four-momentum belongs to the light cone;

Microscopic causality.
Causality is imposed in the meaning of locality, that is, the field operators can commute/anticommute only if they are defined on points separated by space-like vectors: [φ(x), φ(y)] ± = 0 for all space-time points such that A QFT is a Lorentz invariant if and only if it is covariant both in and out of the light cone. "In cone" means that the momenta are physical, i.e., they belong to the forward light cone in momenta space. Analougously "out of cone" means that the considered momenta are not physical. Under these axioms, QFT results totally defined by the vacuum expectation values of the product of field operators. These functions are called Wightman functions and are defined as: where ∆x i = x i−1 − x i . They depend only on the space-time points differences because of the translation invariance postulated in axiom 1 and 3.

Complex Lorentz Group
In this part of the work, illustrating the Jost argumentation, we follow the exposure of [12], which is based on the fact that the space-time inversion is included in the connected component of the complexified Lorentz group. Hence it is necessary to review its properties.
The real Lorentz group is the real set of the 4 × 4 matrices that preserve the Minkowski metric: that is SO (1,3). This group splits in four disconnected components L ↑ + , L ↓ + , L ↑ − , and L ↓ − classified according to the sign of the time changing coefficient Λ 0 0 and the determinant of the matrix.
Since it is not possible to continuously change the determinant of a transformation, the space inversion matrix is not connected to the identity, and from this there is the necessity to employ the complexified Lorentz group, in particular its covering group: L ↑ + = SL(2, C). It is well known that SL(2, C) has two fundamental representations, given by the spinors transformation rules: where spinors are divided in dotted and undotted. Every SL(2, C) representation uses a spinor with m undotted and n dotted indices. The covering group of the complex Lorentz L + subgroup is given by the direct product SL(2, C) ⊗ SL(2, C). In this context it is possible to identify a class of transformations that continuously change from the identity to the space-time inversion:

Axiomatic Cpt Theorem Demonstration
In this section we illustrate the core of the Jost argumentation. Considering a ψ (m, n) (x) spinor field, it transforms under the action of the Poincaré group as: where S mn (Λ) −1 is an opportune matrix determined by the Λ transformation.
The generic Wightman function becomes: Due to the translation invariance, the Wightman function can be written as: Now it is necessary to prove the existence of a complex analytic continuation of these functions. One must consider the action of translation on a Wightman function: where the identity expansion I = pr p r p r | has been introduced. From the previous relation it follows that the momenta p j and the space-time differences ∆x j are conjugate coordinates. As a result, one can write a generic Wightman function as a Fourier transform.
where the p j momenta are physical, that is they belong to the forward light cone, as it is stated by the axiom 4. The exponential present in this Fourier decomposition must provide a decaying factor in order to guarantee that this distribution behaves as an analytic function. This condition is satisfied when the real part of the exponent is negative: Re(−i p j ∆x j ) < 0. This is possible only if Im∆x j = 0 and in particular Im∆x j is defined on the backward light cone, and this domain is called T k−1 . Now one can use this result to generalize the Lorentz symmetry. Equation (12) remains valid for determined complex valued points ∆x j . Now it is possible to use a theorem by Bargman, Hall, and Wightman [13], which states that if Equation (12) is still valid, the function W (k) has a defined and unique analytic continuation on the domain ΛT k−1 = T ′ k−1 . T k−1 is the domain of validity for Equation (12) and Λ are elements of the complex Lorentz L + group. The {x i } points belonging to this set are called Jost points and are characterized by the feature that all sums i λ i (x i −x i+1 ) are space-like vectors for every λ i ≥ 0 and i λ i > 0. From this it follows that i λ i x i ∼ 0 implies x i ∼ 0. Since the complex Lorentz transformations include space-time inversion, from Equation (12) one can obtain the relation: It can be shown that the factor (−1) L is the correct form of the matrix S mini (Λ) for the space-time inversion. The next step of the demonstration considers the k points vacuum expectation functions. From the characteristic of fermionic fields, it is possible to demonstrate the following equality: if N fermionic fields are involved. The N coefficient derives from the spinstatistic theorem, that is from the T transpositions of the N fermionic fields This condition is called weak local commutativity at Jost point. Using Equation (15) written for the fields: together with relation (16) one obtains: In terms of Wightman functions the previous equality becomes: This last result is obtained considering weak local commutativity with spacetime reflection. From Equation (18) using the complex conjugation, it is possible to obtain the relation: where N = i N i and: is the generic transformation rule for fermionic fields.
From the equality: and the vacuum invariance under CP T , the theorem follows. Therefore weak local commutativity is a necessary and sufficient condition together with Lorentz invariance to guarantee CP T symmetry preservation.
A proof based on a more complete mathematical formulation can be found in [14].
The Lagrangian approach is based on constructive method: It is shown that all the physical meaningful Lagrangian terms must be CP T even. We will give a review of this result for the most significative fields: • Scalar: • Fermionic field with spin = 1 2 : • Bosonic fields with spin = 1 In previous relations, the coordinate inversion follows from P and T . The operator C instead acts on the fields, causing the hermitian conjugation.
The four-gradient under the θ = CP T action remains unchanged: For a generalized tensor field T µ1...µn (x), the CP T transform is: as it can be shown for the Faraday tensor field F µν = ∂ µ A ν (x) − ∂ ν A µ (x) using the previous result (25): Now to proceed with the proof one has to notice that the Lagrangian density must be a Lorentz scalar, so all the spinorial indices must be contracted. This means that spinors must be always paired to form bilinears: ψψ, ψγ µ ψ, ψσ µν ψ,

etc.
For instance, in the case of a spinor bilinear χψ, the θ = CP T action gives: where the anticommutation of spinor fields and of γ 0 and γ 5 matrices has been used and the transposition on the last equality can be neglected since it is applied to a scalar. The last equality shows that the general CP T transformation rule for spinor bilinears is given by (23). With a simple generalization the following relation follows for generic Lorentz invariant Lagrangian densities: where the contraction of Lorentz indices and even number of fermionic fields has been used. A physics meaningful Lagrangian must be not only Lorentz invariant, but even hermitian L(x) = L † (x). This last property is motivated by the necessity to guarantee the theory of unitarity. This feature used together with Equation (30) justifies the relation: Finally this last result, with the assumption that the fields interact in the same space-time point leads to the relation for the action: hence the theorem assertion is proved.

CPT Violation Implies Lorentz Invariance Violation
In this section we resume the Greenberg argumentation [7], which states that a CP T odd QFT necessarily also violates the Lorentz covariance. The proof uses the previously introduced algebraic QFT principles.
First it is necessary to define the general τ function as: where σ is a permutation of a k indices set and Θ is the time ordering operator. The proof is based on the idea that if CP T symmetry is violated for every Wightman function, than the related τ functions and therefore even the complete QFT theory violate Lorentz covariance. The demonstration proceeds using the idea that if the τ functions are covariant, the time reversing Lorentz transformation that leaves invariant the Wightman functions must not affect even the τ functions. This statement must be true even for space-like points.
Therefore there must be a complex Lorentz transformation (time inversion) that does not change the Wightman functions but makes negative the time differences between successive points. From the τ invariance it follows that Wightman functions are equal if their fields order is totally inverted, as stated before. Hence to preserve the τ functions covariance it is necessary to have weak local commutativity of the Wightman functions, a sufficient and necessary condition for the CP T invariance and this demonstrates the assertion. However the idea that CPT violation automatically implies LIV was confuted in [15] and the argument has been widely debated in the literature [16,17,18,19].

Consequences of CPT Symmetry
In this section we briefly analyze the consequences of CP T symmetry, focusing in particular on the relation between particles and antiparticles.
The first part of this section is devoted to the properties of the general physical eigenstates. It is not necessary to deal with the exact form of these states, but only to consider the Lorentz invariant quantum numbers. The quantum numbers are supposed to arise from continuous symmetries, which are associated to the conserved currents by the Noether theorem. The conserved currents are composed of the product of fields and fields derivatives and can be used to form hermitian operators whose functional form is space-time independent.
Using the property of CP T symmetry of a generic n indices tensor, it is easy to obtain the following relations for a conserved current: Considering the charge: it follows: The energy-momentum tensor transforms under CP T as: From this relation and from the definition of momentum: it is possible to obtain: In the same way, starting from the CP T transformation relation valid for the angular momentum tensor: The Lorentz group generator is assigned: it is possible to obtain: Considering now a generic one-particle state |p µ , j, j p , q , where p represents the four-momentum, j is the total angular momentum, j p is the angular momentum projection along the momentum direction, and q is the particle charge, it is possible to obtain: The CP T transformations for p µ and q were obtained before. The total angular momentum is given by the sum of angular momentum and spin J µ = L µ +S µ . The total angular momentum is represented by the three indices tensor J µνα (x), so the j p projection transforms as θj p θ † = −j p . On the contrary, the total angular momentum is not affected by the CP T action since it is associated to the J 2 squared operator, that is invariant. Now we give a brief review of the effects of CP T on coupling constants present in the Hamiltonian. There are a lot of issues in defining couplings, but regardless of these ambiguities, one can for instance consider the interaction Hamiltonian: If A CP T (x) = A(x) then the coupling constants λ and λ may differ. A CP T (x) = θA(x)θ † is the CP T conjugate and the A operator is supposed to not contribute to the anti-particle matrix element and viceversa: The interaction Hamiltonian transforms therefore under θ = CP T as: The generalized Hamiltonian, which must be CP T even, assumes the explicit form: with H 0 CP T = H 0 . The CP T the invariance of H implying therefore that λ = λ, so the coupling constants for particles and antiparticles must be equal.
Another observation can be made on masses, indeed from the dispersion relation: and from the conjugation property of p µ it results that every particle must have the same mass of its related antiparticle. For the same reason a particle and its antiparticle must have the same dispersion relation in every CP T even theory.
To conclude this section we underline how it is possible to construct CP T odd terms for an Hamiltonian: It is possible by introducing different masses or dispersion relations or coupling constants for particles and the respective antiparticles. Finally it is always necessary to violate the LI, even if this is only a necessary but not sufficient condition. More extensive analysis can be found in [11].

CPT Violation Motivations
The first motivation to investigate eventual violations of CP T symmetry was given by Hawking [20,21]. He argued that gravitational effects can influence the predictability of future and with this the unitarity of a theory. Indeed black holes with their event horizons can limit the observer knowledge of physics. This argument evolved in the space-time foam idea [22,23], where a pure state evolves in a mixed state because of the quantum background effects causing decoherence.
This foam is supposed to be made of fluctuations of the spacetime that can manifest such as Planck scale black holes. The CP T violations can manifest as unitarity violations, since part of the physical information can disappear inside these black holes event horizons. Hence even in the low energy limit and in the flat space-time, approximation gravitational effects can spoil the validity of CP T symmetry. Defining the S scattering matrix as the operator thank whom an in state evolves to an asymptotic out state: and the density state matrix as: the effects of the quantum foam on the states evolution can induce decoherence and unitarity violation in the following way: whereS is a non unitary operator and thereforeS = SS † . This ensures a violation of strong CP T invariance intended as the standard formulation of this symmetry. However it is not possible to exclude the validity of the weak version of this symmetry. In this last case, CP T is preserved in the state probabilities of the asymptotic states obtained from the evolution of pure initial quantum states. This hypothesis can be supported by ideas muted by black hole thermodynamics. From black hole evaporation the disappeared information can be re-obtained and unitarity can be restored. A more detailed discussion can be found in [23]. The OPE is introduced giving the collection of the series coefficients: where I is an index set suitably defined for this case.
The OPE coefficients must satisfy the following axiomatic requirements: 1. Locality and covariance.
The indices must transform under the action of a generic ρ transformation as: where M ′ is the transformed manifold; 2. Identity element..
where i stands for the identity I in the k-place; for an opportune π 0 permutation; 4. Commutativity-anticommutativity.
with F (i k ) and F (i k+1 ) indices related to the bosonic or fermionic nature of the field involved; 5. Scaling degree. The singularities on the field product must have a positive frequency;
An opportunely defined notion of associativity is required; 9. Analytic dependence upon the metric.
The OPE coefficients must be regular functionals of the space-time metric.
In this work we avoid discussing in more detail mathematical issues that can arise from previous requirements. We limit the discussion to underline that the Wightman functions can be written as: and resorting to the OPE an axiomatic construction of a QFT can be obtained. ∀a i and ∀f i must be true that:

Construction of a QFT from OPE Coefficients
3. Relations arising from OPE. 4. Anticommutation Relation.
for a properly defined index K(i 1 i 2 ); 5. Positivity.
The previously cited spectrum condition can be written in the following form: where W F stands for wave front and Γ n (M ) is a properly defined set of points belonging to M .
In this QFT scenario, the following lemma is true:

Lemma
The map α L : where ψ (i) is an embedding isometry that preserves causality, spin structure, and orientation.
Resorting to this lemma, one can prove the CP T theorem generalized for curved spacetime, which can be formulated in the following way: and for odd n: where F (i) is the number of spinor indices and F (i) = U (i) + P (i) mod 2 with U (i) the unprimed (that is the unchanged) and P (i) the primed (that is the transformed) indices, and θ CP T M is the antilinear isomorfism already introduced.
The key idea of the demonstration consists in choosing the function ψ (i) of the previous Lemma equal to the multiplication of the antilinear immersion: It can be demonstrated that the * -homomorphism constructed in this way is an isomorphism proving the assertion.

CPT and Gravity
Since the discovery of antimatter, many physicists have wondered about antimatter behavior in a gravitational field. Many theoretical arguments sustain that particles and antiparticles have the same gravitational properties. Indeed CP T symmetry affects a particle's internal quantum numbers but presumably leaves unaffected the gravity sector. This statement is supported by the fact that the motion in a gravitational field seems to be symmetric under the time reflection, so thanks to the crossing symmetry of QFT a particle falling follows the same geodesic as for an antiparticle motion, with reversed time. Obviously gravitational behavior can be more complex if the mass is different for particles and antiparticles. What can happen is that antimatter is self-attractive as matter, but the gravitational mutual interaction may be either attractive or repulsive. Arguments against any gravitational repulsive interactions are present in [29,30,31], but were questioned and criticized in [32,33,34,35], so this issue is still under debate. In some proposed interpretations of antimatter gravity, the gravitational mass of antimatter is supposed with negative charge, meaning antimatter is self attractive but interacts repulsively with ordinary matter [35,36].
Instead other authors have argued that antimatter can be self repulsive [37,38].
The issue arises from the interpretation of the charge present in the formulation of General Relativity (GR) equations of motion. Here we discuss the idea presented in [39], based on the consideration that in GR, the real charge is not the gravitational mass m g , but the four momentum p µ = m g The Einstein equation results therefore invariant under the action of CP T , as it follows from the transformation rule for even-indices tensor fields. Antimatter behaves therefore like ordinary matter and must be gravitationally self attractive. This result can be shown even more directly resorting to the geodesic equation of motion: where the inertial mass m i and gravitational mass m g are explicitly written.
From the form of Equation (69), it is simple to notice that the charge is given by the term m g Γ µ αβ dx β dτ . The Christoffel symbol: is a 3-indices object, so it results to CP T odd. The differential form dx µ is CP T odd, therefore the terms: are odd under CP T action. The CP T operator applied on the test body acts on the following terms of Equation (69): The CP T action on the gravitational field generator transforms the term: The action of the CP T operator on the test body and gravitational field generator leaves unaffected the equation of motion (69), so antimatter is self attractive like ordinary matter. It is interesting to notice what happens instead for antimatter in a gravitational field generated by ordinary matter, or vice versa. The motion of antimatter in an ordinary gravitational field is described by the equation: The opposite situation, that is the motion of matter in a gravitational field generated by antimatter, is described by: In both cases the equation of motion becomes: so the gravitational interaction between matter and antimatter is repulsive.

CPT Violation and LIV Research
The

Very Special Relativity
The first of these effective theories is the Very Special Relativity (VSR) theory of Coleman and Glashow [40]. In this framework the supposed effects induced by the quantum background modifies the Maximal Attainable Velocity (MAV) of every particle species. This effect implies a modification of a particle's dispersion relations: where ǫ represents the correction factor for the MAV of the particle species taken into account. This modification implies the new propagator form: .
The change introduced in the dispersion relation and in the particle's propagator is obviously Lorentz violating, but CP T even, since there is no differences between the Modified Dispersion Relation (MDR) of the particle and the related antiparticle. This theory is conceived in order to be isotropic, but in a preferred reference frame. This means that the Lorentz symmetry is broken.
A following formulation of VSR [41] is The DISIM b (2) generalization [42]. The implementation of this second formulation of VSR leads to the coupling of fields with fixed background tensors and we refer to the next section for a more complete exposition of this approach.

Standard Model Extension
The second framework used to study LIV and even CP T violation effects is the Standard Model Extension (SME) [43]. This model is based on the idea that in some string theories, particular operators can assume vacuum expectation values different from zero and can therefore can introduce a breaking of space-time isotropy [24]. In this theory, the background effects are studied by supplementing all the possible LIV or CP T odd operators, preserving the gauge symmetry and microcausality with respect to the positive energy and four-momentum conservation law. Moreover, the SME formulation is conceived in such a way so as to guarantee the existence of relativistic Dirac and non- where the CP T odd terms are those with an odd number of contracted Lorentz indices.
Even in SME, the Lorentz covariance is broken modifying the particles dispersion relations. Nowadays SME appears to be the most complete framework to study LIV and CP T violation, since its rich phenomenological predictions are capable of testing the validity of these fundamental symmetries. As a final remark, the SME framework standard physics is perturbed in every regime and presumed phenomenological effects could be detectable in low-energy experiments. Indeed many of the tightest bounds on SME coefficients are set by low-energy atomic and optic experiments [44].

Theories Preserving Covariance
Another approach to the study of LIV or CP T violations consists in attempting the construction of complete physical theories. This is the case of Doubly Special Relativity (DSR) [45,46,47,48], in which LIV is studied formulating a modification of Special Relativity (SR) that attempts to include in the theory formulation another invariant quantity, the Planck length, in addition to the light speed. In this model the momentum space and not the spacetime is supposed to be the fundamental structure describing physics. The concept of absolute locality results relaxed and different observers feel a personal spacetime structure, which presents an energy dependence. Space-time description is constructed by every observer in a local way, losing its universality and becomes, therefore, a personal and auxiliary concept that emerges from the fundamental momentum space, where the real dynamics takes place. The dispersion relations result again modified and are implemented in the action formulation: with the particle's MDR M (p): where f (| p|, E) is the perturbation function introduced by the quantum gravity phenomenology. Moreover to determine the momentum space connection, the interaction processes kinematics are modified. More in detail, a modified composition rule for momenta is defined as: where g(p, q) represents a perturbation of the usual momenta sum. In an interaction process, a particle is therefore supposed by influencing the other incoming particles proportionally to its energy. In this theory the Lorentz symmetry is modified and not simply broken.
Finally, another scenario where the LIV effects are analyzed is the one denoted as Homogeneously Modified Special Relativity (HMSR) [49,51,50]. This theory was developed in order to preserve, in addition to the SME symmetries, and is also a covariant formulation. Indeed in this model the Lorentz symmetry is modified and the Lorentz group results amended in order to preserve covariance. In this last theory the interaction with the background is geometrized, indeed the dispersion relations are modified perturbing the kinematic. This implies a modification of the underlying geometry of SR. This model introduces a peculiar form of MDRs that present a perturbation homogeneous of zero degree: In this way the dispersion relation satisfies the requirements act to be a The dispersion relations are supposed to be different for every particle species, as in VSR [40]. The model built in this way is an extension of the SM, which preserves the internal symmetry SU (3) × SU (2) × U (1) and does not introduce any exotic particle or interaction and is CP T even. Moreover the isotropic LIV corrections are conceived in order to obtain a modified Lorentz covariant formulation that implies the preservation of space-time isotropy.

CPT Violation and LIV Geometry Framework
A theoretical model including the possibility of LIV or CPT violation usually also aims to find an appropriate geometry to formulate an extended version of the "usual" theory, making possible some kind of unification between quantum field theory and GR, i.e., QFT and gravity [52]. Indeed it is well known that the conventional Riemann geometry, even in the extended Riemann/Cartan formulation is incompatible with Lorentz or CP T violation. The natural setting of such theories appears to be the Riemann/Finsler geometry [53,54,55,56]. This idea has been investigated in the SME framework [43], in which the free particle's kinematic is described by the geodesic motion in a spacetime that emerges due to the coupling of the SME background tensors with the fields describing particles. Finsler geometry is used even in SR modification theories, such as DSR [45,46,47,48], VSR [40,41], and HMSR [49,51], where the Lorentz symmetry is only modified and not merely broken as in SME. In DSR formulation local relativity states that the momentum space is the fundamental structure at the basis of the physical processes description, indeed the momentum space is supposed curved and the underlying geometry appears to be the Finsler one, in which the framework presents an explicit dependence on momenta [57,58]. Even VSR [59] and HMSR [49] are settled in the Finsler geometry, which naturally emerges in this case as the appropriate one from the kinematic perturbation introduced via the MDRs [60].

Search for CPT and LIV Violation in Astroparticle Physics
Since in some theoretical models the supposed quantum gravity effects are expected to be more visible in the high energy limit as LIV perturbations, astroparticle physics can be a useful framework to conduct search on possible departures from CP T symmetry or eventually from the Lorentz covariance.
Indeed the tiny effects caused by quantum gravity are supposed to be more visible at the highest energies reached for instance in the astrophysics sector.
Moreover these effects are supposed to add up during the particle's propagation for the cosmic distances involved in astroparticle physics. In this section, the presumed LIV effects could manifest as modifications of the photons and the Ultra High Energy Cosmic Rays (UHECRs) propagation and finally LIV can influence the neutrino oscillations.

Ultra High Energy Cosmic Rays
In this section we give a review of the expected effects related to the Cosmic Rays (CRs) physics, in particular the enlargement/modification of the Greisen Zatsepin Kutzmin (GZK) opacity sphere [61,62] and via a photopion production, that through a ∆ particle resonance is dominant for the high energy limit: Instead heavy nuclei CR can interact via a photodissocitation process: where A is the atomic number of the considered CR bare nucleus.
The creation energy of UHECRs is limited by the dimensions of the astrophysical object that accelerated these particles. As a consequence of this upper limit and of the dissipation mechanism, the UHECRs originated outside a sphere centered on Earth and of a determined radius can be detected only under a determined energetic threshold. For this reason, free space is revealed to be opaque to the propagation of massive particles such as UHECRs.
There are some experimental hints that the GZK opacity sphere could be modified in respect to theoretical predictions, as discussed in the work [63], where UHECR are correlated with cosmic neutrino sources located farther than the classic opacity length. This modification can be justified by the introduction of LIV perturbations in the standard particle physics. This possibility is investigated for instance in works as [64,65,66,51], in which the propagation of the light massive cosmic rays (protons and similar particles) is analyzed. The energy dissipation mechanism caused by the interaction with the CMB determines a modification of the attenuation length or the mean free path of a proton, defined as the average distance that the particle has to travel in order to reduce its energy by a factor of 1 e . The inverse of the attenuation length is given by [64]: where µ is the impact parameter µ = cos θ and n(E) is the density function of the CMB: A Planck's black body energy distribution, σ pγ (E) is the cross section for the proton-photon interaction and E th is the threshold energy for this physical process. K(s) is the inelasticity of this reaction, defined as the energy fraction available for the production of secondary particles. One can define the elasticity as the energy fraction preserved by the primary particle after the interaction where E in is the primary incoming energy and E out is its energy after the interaction. From this relation the definition of the inelasticity Using the following relations: where the Mandelstam variable s is computed using the photon four momentum defined in the rest frame of the cosmic ray (ǫ ′ , p ′ γ ). In the high energy limit the proton velocity can be approximated as v p ≃ 1, with ds = −2E p ǫ d cos θ and using the explicit form of the n(E) Planck distribution, the previous relation (88) can be simplified, obtaining: In the classical Lorentz invariant physics scenario, the inelasticity is given by the relation [64]: The introduction of LIV perturbations in the standard physics determines kinematic modifications, which reduce the phase space part available for the photopion production, amending the inelasticity K (91) in Equation (88). These kinematic modifications determine a reduction of the inelasticity as a function of the proton energy, implying a dilatation of the attenuation length [64,65,66].
The introduction of LIV can influence even the propagation of heavy CR, modifying the photodissociation process they undergo. For its peculiar feature photodissociation determines a change in the chemical composition of the CR considered, and this means that after every interaction the CR changes. However this process can be described as a dissipation effect grouping the CR in families based on similar chemical composition. Since the photodissociation is a relatively rare process and in the hypothesis only few nucleons are lost in every process, one can suppose that a UHECR belongs to the same family even after some interactions, so one can study this phenomenon as the proton energy dissipation and can see that even in this case the attenuation length is enlarged [67].
As already underlined, a CP T odd theory is also a LIV theory, but in the presumed GZK modification the introduction of explicit CP T violation is not directly required.

Time Delays
Among other theoretical and experimental evidences that SR is modified or breaks down in the ultra-high energy scenario, time delay in the propagation of high energetic particles plays an important role.
It is well known that intense gravitational fields can modify the propagation of photons traveling nearby massive objects because of space-time dilation and this effect is a classic test of GR [68]. However some experimental results hint the possibility that photons propagating with different energies and in the absence of intense gravitational fields can present an analogous time delay.
In this section, the best constraints on the quantum gravity induced effects magnitude can be posed using astrophysical sources with distant, bright and very energetic astrophysical sources such as Gamma-Ray Bursts (GRB), flaring Active Galactic Nuclei (AGN), and Pulsars (PSR). For instance this effect was observed between gamma-ray flares, with different energies, originated from the Markarian 501 galaxy center [69]. This phenomenon can be explained introducing modifications to SR that imply a departure from the Lorentz/Poincaré symmetry perturbing the free particle's kinematic both in a Lorentz breaking [70,71,72] or in a SR modification scenario [73,74,75,76]. In the following explanation of this phenomenon, we briefly review the argumentation of the SR kinematic symmetry modification approach.
As shown before modifying the dispersion relations, the resulting geometry acquires an explicit dependence on the particle's momentum, which is propagating particles probe of different spacetimes as a function of their energy. In DSR or desitter projective relativity the line element acquires the explicit form: where n(E) is the function expressing the metric dependence on the particle's energy that is a departure from the standard Minkowski geometry, and is supposedly increasing with the energy, hence it can be written as: where E P ≃ 10 19 GeV is the Planck energy, which is the characteristic suppression scale for the LIV phenomenon, and k n are positive adimensional coefficients.
From the previous relation (92), one can compute the photon propagation velocity as a function of its energy, considering that this particle propagates along the null direction: finally one can obtain: where the metric correction factor n(E) acquires the role of a refractive index and n(E) > 1. Taking into account that the velocity of a photon depends on its energy, one can obtain the time difference ∆t between two gamma ray burst with different energies to travel the same distance L: where the standard light speed is posed c = 1. This result is in accordance with the experimental evidence that a higher energy gamma-ray burst is delayed with respect to a lower energetic one, since n(E 2 ) > n(E 1 ) for E 2 > E 1 . A complete phenomenological discussion on this topic can be found in [77].
An analogous time delay effect is theoretically predicted even in massive particle's propagation. Indeed starting from a massive particle's MDR written as: where f ( p, E) is the LIV-induced perturbation as a function of the particle's four-momentum. Now it is possible to compute the particle's four-velocity: obtaining a result that depends on the particle's momentum.

CPT and LIV in Neutrino Physics
The consequences of presumed LIV on neutrino phenomenology are illustrated, for instance, in [50,51,78]. The idea of introducing LIV in neutrino physics was first investigated in 1999 by Coleman and Glashow [40].
The presence of tiny perturbative corrections, breaking Lorentz invariance determines a modification of the dispersion relations, introducing the idea of a MAV (that can differ from one particle to the other) modifying the neutrino propagation and, consequently, the oscillation phenomenon. This aspect is common to different studies of LIV, with and without the violation of CP T invariance, which followed over the years.
Some of these works consider also very exotic scenarios that are difficult to reconcile with present phenomenological data or predicting phenomenogical effects that are difficult to distinguish from the standard scenario. This is the case, for instance, of the already cited VSR model, in which the "magnitude" of Lorentz violating effects is limited by the connection between LIV and the violation of some internal symmetry of the theory (CP symmetry in this case), as underlined in [41]. In some of the papers adopting this model [40,79], the authors also consider a possible "natural origin" of neutrino oscillation, different from the usual one, which would require neither the introduction of additional states, nor the lepton number violation, but would imply ultra-relativistic neutrinos (γ >> 1). In this case, VSR and conventional neutrino masses would be indistinguishable, but VSR effects could be significant near the beta decay endpoint, where neutrinos are non ultra-relativistic [79].
The possibility of ultra-luminal neutrinos and the limits on LIV magnitude coming from high energy cosmic neutrino phenomenology has also been investigated in papers like [80], in which the Lorentz violation hypothesis is introduced in a more general quantum gravity framework, following an approach based on EFT formalism.
We have already explained in previous sections that the reference work for this approach is the SME [43], propagation. Similarly to the VST case, the presence of these LIV operators introduces deviations from the standard physics predicted c for every particle's MAV.
For a particle of species i and mass m, the MDRs can be rewritten in the form: wherem i is a sort of effective mass term, containing the corrections generated by LIV operators: The energy dependent dimensionless (here and in the rest of the paper we adopt, for simplicity, the natural units = c = 1) δ i coefficients are functions of the lagrangian parameters, which can be different for different particle species and determine the MAV for every particle. The The coefficient corresponding to [d] = 4 dimension operator can be supposed identically zero, or at least suppressed with respect to the ones of the [d] = 5 and [d] = 6 operators (see e.g, [83]). In any case the coefficient of this operator entering in processes in which neutrinos and electrons are involved has been seriously constrained by the analysis [84] of the IceCube data [82], which fixed k νx,e,0 < 5.2 × 10 −21 .
Concerning the remaining LIV operators, the ones with odd mass dimensions (for which n is an odd number) are also CPT -odd, while the other operators (for which n is an even number) are associated to terms in the Lagrangian that violate the Lorentz invariance but preserve CP T symmetry, as discussed in [81], in which all gauge-invariant Lorentz and CP T -violating terms with dimensions up to nine in the quadratic Lagrangian density are classified.
As already underlined, one can distinguish two different classes of theories: The ones assuming the violation of CP T and Lorentz symmetry and the other, which predict LIV but are CP T -even. For both classes of models, the possible constraints on the different operator coefficients coming from experiments with fermions have been analyzed e.g. in [85] and, in the specific case of high energy astrophysical neutrino data and CP T -conservation, in [86].
Using, the notation (adopted also by [80]) of Equation (100) In principle one should introduce also a helicity dependence in the coefficient of the n = 1 operator, because in this case (due CP T -violation) one has to distinguish the LIV contributions associated in the Lagrangian to terms (allowed in Standard Model extensions) involving right handed neutrinos. However the potential contributions of these terms to electron coefficients are negligible because they are very strongly constrained by experimental data, mainly by Crab nebula observation [87,88,89,90,91].
The possibility of using neutrinos from astrophysical sources to test the hypothesis of theories predicting LIV and eventually also CP T -violation has been widely exploited, not only in the studies that we just discussed (most of which rely on the investigation of 2014 Ice Cube data [82]) and that will also continue in future (concerning these kind of analyses relying mainly on future Ice Cube data, see also [92]), but also with the use of other sets of data and considering different astrophysical sources. For instance, the case of the study of the propagation of neutrinos from a Core-Collapse Supernova (looking for variations in time of the order of a few milliseconds [93]). This includes even the search for a possible LIV-induced spread out of the time profile signal characteristic of a ν e neutronization burst from the next galactic supernova [94] and also of the LIV and CPT violation tests by the oscillation pattern observation in experiments measuring atmospheric neutrinos [95,96].
The most promising opportunity is probably offered by the recently born and increasingly relevant multimessenger astronomy, looking for very high energy neutrinos emitted by astrophysical objects (like blazars, neutron stars, and so on) together with other electromagnetic and/or gravitational signals. Examples of this kind are offered by [97,98].
In the last few years, the investigation of astrophysical neutrinos as probes for potential signals of Lorentz and CP T symmetry violation has been flanked with analogous studies performed with artificial neutrino sources, produced mainly by accelerators and also partially by reactors.
A first interesting example is given by [99], in which the authors discuss the possibility of using a generic Short Baseline (SBL) accelerator neutrino experiment to look for a possible anomalous energy dependence in the spectrum (differing from the one predicted by usual massive neutrino oscillation) and for a dependence on the direction of neutrino propagation. They apply this idea in particular to the LSND data, obtaining a possible non zero value of the order of (3 ± 1) × 10 −19 GeV, for a combination of LIV coefficients. A similar study has been performed by the MiniBooNE Collaboration in [100] using both the ν e and ν e appearance data of this experiment and looking for sidereal time variation effects that could be induced by CPT and Lorentz violations. These effects do not seem to be present in most of the data, but one can find a competitive upper limit of the order of 10 −20 GeV (similar to that derived in [99]), for some combination of the SME model coefficients, which would eventually contribute to the ν µ → ν e oscillation channel. A search for such sidereal modulation effects in the neutrino interactions rate has also been performed more recently [101], using the T2K on-axis near detector (therefore for a short value of the baseline, even if in the case of a near detector built for a long baseline experiment (LBL)).
The analysis gave no evidence of such a signal.
An even more interesting possibility is offered by the LBL accelerator neutrino and the medium baseline reactor antineutrino experiments, as discussed since a few years ago in [102] and in [103], respectively for the first and sec-  [109,110], and NOνA [111]. The LBL experiments have already given fundamental contributions to our understanding of neutrino properties, from the determination of the "atmospheric" oscillation parameters θ 32 and ∆m 2 32(31) , to the indication of θ 31 = 0 (complementary to the studies performed by reactor antineutrino experiments). The new generation of LBL is designed to attack some of the main open questions of neutrino physics of great impact for all elementary particle physics and astrophysics, like the search for CP violation in the leptonic sector and the determination of neutrino mass ordering. A very remarkable result has been obtained recently [112] by T2K, which found important constraints on the possible CP violation phase (δ CP = 0 and δ CP = π at almost 3 σ), confirming the hints derived already by other LBL and mainly by NOνA. [113].
The LBL accelerator experiments have automatically made the possibility of studying both the neutrino and antineutrino channels and, therefore, they are the ideal framework for CP T studies. Moreover this new generation of experi-ments is characterized by high values of the energy E and of the baseline L, and, therefore, they also offer the opportunity to search for signals of LIV (which is expected to introduce a subleading effect in the oscillation pattern proportional to L × E). In addition to this, the high intensity of the beams guarantees high statistics, essential for every investigation looking for tiny effects.
A recent interesting example of such kinds of studies is given by [114].
The authors perform their analysis in the already cited framework of SME [43], restricting their analysis to the possibility of LIV together with CPT -violation, like in [115,116]. Following the line of [117] and of [118], in which the possible operators appearing in the SME Lagrangian are classified according to their mass dimension (as already discussed for [80]), and keeping into account only renormalizable terms (that is contributions of mass dimension [d] ≤ 4), the LIV part of the neutrino lagrangian density can be written as: The different contributions can be reorganized, considering the fact that only the left handed part of the neutrino wave functions takes part in the interactions and keeping only the terms that violate CPT, in addition to LI. In this way the CPT violating and LIV corrections can be parametrized in terms of just two parameters: (a L ) µ αβ = (p+q) µ αβ and (c L ) µν αβ = (r+s) µ αβ , which are represented by hermitian matrices in the flavor space. The choice adopted in [117] and pursued also in [114] is that of choosing µ = ν = 0 in order to realize in the simplest way (but not the only possibility) isotropy for what concern space-time rotations in a fixed preferred reference frame with respect to the isotropy preserving coefficient αβ . An important note is that the term (a L ) µ αβ explicitly violates Lorentz invariance introducing a preferred direction in spacetime. In this way one can write the LIV part of the Hamiltonian using the matrices a 0 αβ and (c) 00 αβ : with the formal correspondence a 0 αβ = √ 2G F N e ǫ αβ . The values of the non diagonal matrix elements, that would contribute to lepton flavor violation, are of course limited by the data of oscillation experiments, like the SuperKamiokande atmospheric neutrino data [119], which put the following limits (in units of GeV): |a 0 eµ | ≤ 2.5 × 10 −23 ; |a 0 eτ | ≤ 5 × 10 −23 ; |a 0 µτ | ≤ 8.3 × 10 −24 . In [114] the impact of the different CPT and LIV coefficients on the neutrino flavor oscillation probabilities and consequently on the ν µ → ν e and ν µ → ν µ channels for the two LBL experiments T2K and NoνA and the experiments' sensitivity to these coefficients are analyzed. The main outcome of the analysis is that the a 0 eµ and a 0 eτ , together with the matrix diagonal element a 0 ee are relevant for the ν µ → ν e appearance channel, whereas the ν µ survival probability is influenced (but in a weaker way) by a 0 µτ . Generally speaking the upper limits obtained on these coefficients from the analysis of the T2K experiment are almost one order of magnitude lower than the corresponding ones derived by the atmospheric neutrino data. The limits become more stringent in the NOνA case and even more if one considers a "combined T2K + NOνA" analysis, but in any case they do not reach the level of SuperKamiokande limits. The study performed in [114], for the NOνA case also shows that the CPT violating LIV corrections can affect quite significantly the CP violation and mass ordering studies at LBL experiments, increasing or decreasing their sensitivity according to the values of the LIV parameters phases.
A similar study [120] has also been performed for the case of the DUNE experiment [109].The authors investigated the impact of potential CP T -violating LIV corrections on the experiment sensitivity for the measurement of the θ 23 octant and CP phases. They found that the presence of a non null single LIV coefficient, even with a relatively small value (|a eµ | or |a eτ | = 5 × 10 −24 ) could reduce the DUNE discrimination power, and in addition due to the destructive interference between the phenomenological parameters one wants to extract from the analysis (that is the mixing angle, the usual Dirac CP violation phase δ and the new CP phase of dynamical origin φ eµ /φ eτ ). The reduction of the θ 23 octant resolution would, instead be less significant in case one considers both the LIV parameters (a eµ and a eτ ) together in the analysis.
As already discussed, there are models in literature in which the modification of the traditional Lorentz invariance symmetry is introduced without breaking the fundamental CP T invariance. They also predict interesting phenomenological consequences that can be tested by using various experiments, both with natural and artificial neutrino sources. We refer in particular to the HMSR model [49,50,51,78,121], in which the LIV has a geometric kinematic origin and does not break space-time isotropy. This feature is very important for phenomenological analyses because it guarantees the possibility of looking for isotropic corrections to the "traditional" experimental signature, without the need of fixing any preferred reference frame.
The MDRs cause a variation in the neutrino propagation and, consequently, a modification in the neutrino oscillation probabilities, which has been explicitly computed for different neutrino flavors in [50,51,78]. The LIV-induced corrections present an energy dependence different from the one of the "usual" leading term (they are proportional to L × E, instead of L E ). This introduces a partial deformation of the oscillating neutrino transition probability and of the corresponding experimental spectrum for various neutrino interaction channels.
This effect could in principle be observed in experiments involving high values of the factor L × E, that is the product of the neutrino energy E and of the baseline L between the neutrino production and interaction points. The amplitudes of these corrections in the model considered would be proportional to δf ij , which is the difference between the LIV coefficients f i and f j , correcting the dispersion relations for the two neutrino species ν i and ν j involved in the process under investigation.
As shown explicitly in [50,51,78], opportunities has been performed [50,51,78,121,122]. Further analyses can be done about the neutrino telescope (ANTARES [123], KM3NeT [124], and Ice-Cube [92,125,126]) detection of high-energy neutrinos (with E in the range between TeV and PeV), the Auger study of cosmic neutrinos with E > EeV [127], and the potential analysis of high-energy atmospheric neutrinos at the JUNO experiment [128,129]. Some of these studies are in progress and will be published soon.
Finally, to complete the theoretical and phenomenological framework, it is worthwhile to recall that an interesting general analysis of the limits that can be derived for LIV from the study of high energy cosmic neutrinos can also be found in [130]. Conflicts of interest The authors declare no conflict of interest. References