Next Article in Journal
A Hybrid MCDM Model to Evaluate and Classify Outsourcing Providers in Manufacturing
Next Article in Special Issue
High-Energy Neutrino Astronomy—Baikal-GVD Neutrino Telescope in Lake Baikal
Previous Article in Journal
Fluctuating Asymmetry and Population Dynamics of the Common Shrew, Sorex araneus, in Central Siberia under Climate Change Conditions
Previous Article in Special Issue
Exploring Finite-Sized Scale Invariance in Stochastic Variability with Toy Models: The Ornstein–Uhlenbeck Model
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Predictions of Ultra-High Energy Cosmic Ray Propagation in the Context of Homogeneously Modified Special Relativity

by
Marco Danilo Claudio Torri
1,2,*,
Lorenzo Caccianiga
1,2,
Armando di Matteo
3,
Andrea Maino
1,2 and
Lino Miramonti
1,2
1
Dipartimento di Fisica, Università degli Studi di Milano, via Celoria 16, I-20133 Milano, Italy
2
INFN, Sezione di Milano, via Celoria 16, I-20133 Milano, Italy
3
INFN, Sezione di Torino, via Pietro Giuria 1, 10125 Torino, Italy
*
Author to whom correspondence should be addressed.
Symmetry 2020, 12(12), 1961; https://doi.org/10.3390/sym12121961
Submission received: 4 November 2020 / Revised: 23 November 2020 / Accepted: 25 November 2020 / Published: 27 November 2020
(This article belongs to the Special Issue Cosmic Rays)

Abstract

:
Ultra high energy cosmic rays (UHECRs) may interact with photon backgrounds and thus the universe is opaque to their propagation. Many Lorentz Invariance Violation (LIV) theories predict a dilation of the expected horizon from which UHECRs can arrive to Earth, in some case even making the interaction probability negligible. In this work, we investigate this effect in the context of the LIV theory that goes by the name of Homogeneously Modified Special Relativity (HMSR). In this work, making use of a specifically modified version of the SimProp simulation program in order to account for the modifications introduced by the theory to the propagation of particles, the radius of the proton opacity horizon (GZK sphere), and the attenuation length for the photopion production process are simulated and the modifications of these quantities introduced by the theory are studied.

1. Introduction

Cosmic rays (CRs) are the part of the radiation of extraterrestrial origin made up by charged particles. In this work, we deal only with ultra-high energy cosmic rays (UHECRs), which are made up by protons and other atomic nuclei with energies above 1 EeV most likely of extraterrestrial (most likely extragalactic [1,2,3,4,5]) origin. They can be classified based on their composition in three main groups: the light component principally composed by protons (e.g., protons and He bare nuclei), an intermediate group (e.g., C, N and O bare nuclei), and the heavy component composed of iron type bare nuclei (e.g., Fe and Ni).
Nowadays, it is recognized that the Universe is opaque to the propagation of cosmic rays (CRs) with energy larger than ∼ 5 × 10 19 eV ; this means that such high energy particles can have only a finite free propagation path. Indeed, the propagation of these particles is influenced by their interactions with the Cosmic Microwave Background (CMB). Through these processes, UHECRs dissipate energy during their propagation so that they are attenuated in a way that depends on their energy and nature. This means that UHECRs have a finite free path of order 100 Mpc for 10 20 eV protons and even shorter for higher energies or heavier nuclei, hence, if they propagate for longer distances, they can be detected on Earth only under a determined energy threshold. This phenomenon is named GZK-cutoff after the physicists Greisen, Zatsepin, and Kuzmin, who first theorized it [6,7], posing an upper limit on the energy of UHECR (which are supposed of extragalactic origin) that can be detected on Earth.
The UHECR sources are still unknown; nevertheless, the physics correlated with their propagation can be exploited to investigate new phenomena, such as the hypothetical ones induced by quantum gravity and described by LIV theories. These perturbations can be considered as Planck scale physics residual effects in the low energy limit; therefore, they may be an open window on possible phenomenology induced by quantum gravity. These effects can be considered caused by the interaction of particles with the presumed quantum structure of the background spacetime.
In this work, we introduce a new methodology to investigate this sector exploiting the peculiar features of UHECR physics. Indeed, the CR opacity horizon may be enlarged by new physics, such as the introduction of LIV extensions to the particle Standard Model (SM) [8,9,10,11,12].
In this study, to investigate LIV effects on UHECR propagation, we limit our analysis to the lightest component of UHECRs, which is protons, since the lightest component of UHECRs is the least deflected by magnetic fields and propagates on more straightforward trajectories than heavier nuclei. The LIV effects are expected to be more visible for the intermediate component of CRs and composition data show a mix of light and intermediate elements (except possibly at the highest energies) so all the constraints on LIV obtained for the light component are still valid for the real physical case, for which they are expected to be even more stringent.
There are various LIV theoretical models [8,13,14,15,16,17,18], but we resort to Homogeneously Modified Special Relativity (HMSR) [19] since, for this kind of analysis, it is useful that the theory preserves spacetime isotropy and essential that it can foresee an extension of the particle Standard Model. The introduction of LIV in UHECR physics can determine a kinematical perturbation that modifies the allowed phase space of the photopion production process. The perturbation is then introduced in the SimProp [20] software that becomes able to simulate the LIV modified propagation of light cosmic rays in order to evaluate the GZK sphere radius modifications and the UHECR attenuation length.
This work is structured as follows: in Section 2, we present the LIV model employed here (HMSR) [19]; in Section 3, we introduce UHECR propagation; in Section 4, we illustrate the modifications introduced in the UHECR physics by LIV in order to obtain a modified SimProp software version able to simulate their free propagation in a Lorentz violating scenario; finally, in Section 5, we summarize the obtained results and, in Section 6, we present the future perspectives to extend the analysis to a more realistic CR composition scenario.

2. Homogeneously Modified Special Relativity

In this work, we explore the possibility to introduce LIV to justify a possible GZK phenomenon modification, in order to explore quantum gravity effects. The purpose of our study implies the necessity to use a LIV model that can both modify particle physics phenomenology as required and preserve the isotropy of spacetime. We decided to set our study in the HMSR theoretical framework [19], since this model was developed in order to investigate the opportunity to extend the Standard Model (SM) of elementary particles preserving spacetime isotropy and homogeneity in a LIV scenario.

2.1. Geometric Structure: Hamilton/Finsler Geometry

In order to obtain an extension of the particle SM that preserves spacetime isotropy and does not introduce any exotic reaction, only the free propagation of massive fermions is modified geometrizing their interaction with the assumed background quantum structure. This approach is exploited modifying the dispersion relations (DR) of massive particles, as in the majority of LIV theories [8,13,14,15,16,17,18,21,22].
Every particle species is supposed to have a personal modified dispersion relation (MDR) with the form:
MDR F 2 ( p ) = E 2 p 2 1 f p E = m 2
where the bold letters indicate four-vectors; the perturbation function f is related to the particle species and depends only on the magnitude of the three-momentum vector | p | in order to be rotationally invariant. The perturbation function is chosen homogeneous of degree 0 to preserve the geometrical origin of the MDR. In this way, the F 2 ( p ) function is homogeneous of degree 2 and can be a candidate Finsler pseudo-norm. Here, we are dealing with pseudo-Finsler structures, since the local underlying geometry is the Minkowski one, with metric η μ ν = diag ( 1 , 1 , 1 , 1 ) that is not positive definite. On the contrary, a Finsler geometry is constructed starting from a positive defined metric. Some issues can emerge in defining a pseudo-Finsler structure [23,24,25,26,27], but, in the case of HMSR, the function f is supposed to be of a tiny magnitude compared to the other MDR quantities, so it represents a perturbation; therefore, as a result, it is possible to deal with these issues. In the actual HMSR formulation, MDR are supposed CPT even, since they do not present an explicit dependence on particle helicity or spin. Indeed as a result, it is possible to introduce LIV preserving the CPT symmetry. It is well known that Lorentz violation does not imply the CPT one [28,29,30], whereas the opposite statement is widely debated in literature [31,32,33,34,35].
In HMSR, every particle species has its own personal MDR, with a specific perturbation f. This function depends on the ratio | p | / E which admits a finite limit:
| p | E ( 1 + δ ) for E with 0 < δ 1 f | p | E ϵ 1
This last relation implies that every massive particle admits a different personal maximum attainable velocity (MAV):
c = 1 f | p | E ( 1 ϵ )
as assumed in the first formulation of Very Special Relativity [8]. This HMSR feature is of a very general nature and allows for obtaining a rich phenomenology even in the neutrino oscillation physics [36,37].
The MDR (1) can be promoted to the role of momentum space norm, making it possible to obtain the associated momentum space Finsler one, using the equation:
g ˜ ( p ) μ ν = 1 2 p μ p ν F 2 ( E , p )
The obtained metric presents a non diagonal part that does not contribute to the computation of the MDR, so it can be neglected, thanks to the polarization identity, and the final form of the obtained metric can be reduced to a pure diagonal form, without a changing of the basis:
g ˜ μ ν ( p ) = 1 0 0 t 1 f p / E I 3 × 3
Considering the peculiar form of the perturbation function f that depends on the ratio | p | / E (2), it is possible to correlate at least at the leading order the momentum space with the coordinate space via the Legendre transform. Indeed, the homogeneity of the perturbation f allows for neglecting the derivatives with respect to the momentum, since these terms are a second order perturbation [19]. HMSR introduces therefore a simple way of inverting the correspondence between four-momentum and four-velocity (the demonstration that derivatives of this homogeneous function are negligible is given in Appendix A.). The final form of the coordinate space metric is therefore given by:
g ( x , x ˙ ( p ) ) μ ν = 1 0 0 t 1 + f p / E I 3 × 3
This allows us to correlate the Finsler co-metric resulting from the modified dispersion relation and the Finsler metric describing the modified geometry in coordinate space. The pseudo-Finsler norm can indeed be written as a function of the coordinates:
G ( x ˙ ( p ) ) = F ( p )
and even the related metric in the coordinate space can be computed via the relation:
g ( x , x ˙ ( p ) ) μ ν = 1 2 2 G x ˙ μ x ˙ ν
The metric in coordinate space constitutes the inverse of the one in momentum space:
g ˜ μ α g α ν = δ ν μ
To complete the HMSR modified geometry description, we give the Hamiltonian of the free massive particle as a function of the proper time:
H = g ˜ μ ν ( p ) p μ p ν = F ( p )
consistently with standard Special Relativity. The velocity correlated with the momentum can be computed resorting to the Legendre transformation:
x ˙ μ = p μ F ( p ) g ˜ μ ν ( p ) p ν g ˜ μ ν ( p ) p μ p ν = g ˜ μ ν ( p ) p ν m
The Lagrangian can be evaluated starting from the Hamiltonian:
L = p · x ˙ H = x ˙ μ p μ = m g μ ν ( p ) x ˙ μ x ˙ ν
Starting from the obtained metrics (Equations (5) and (6)), it is possible to introduce the generalized vierbein, in order to satisfy the definitional equations:
g μ ν ( x ˙ ) = e μ a ( p ( x ˙ ) ) η a b e ν b ( p ( x ˙ ) ) g ˜ μ ν ( p ) = e a μ ( p ) η a b e b ν ( p )
obtaining the explicit forms:
e a μ ( p ) = 1 0 0 t 1 f p / E I 3 × 3 e μ a ( p ) = 1 0 0 t 1 + f p / E I 3 × 3
The peculiar feature of HMSR is that every particle feels a specific personal modified spacetime, parameterized by the particle momentum. Every particle lives in a personal curved spacetime and therefore every physical quantity is generalized, acquiring an explicit dependence on the momentum. It is therefore necessary to introduce an original formalism to correlate different local curved spaces, using the generalized vierbein as a projector from the curved spacetime to a common flat Minkowski support space. Here, we report a scheme of how the original mathematical formalism of HMSR works to relate different local curved spaces:
Symmetry 12 01961 i001
where the Greek indices refer to the local curved geometric structures, whereas the Latin ones refer to the common support Minkowski space. Using again the vierbein as projectors from the common support space to the local curved one, it is simple to obtain the generalized Lorentz transformations:
Λ μ ν ( p ) = e μ a ( Λ p ) Λ a b e ν b ( p )
where Λ a b are the classic Lorentz ones valid for the flat Minkowski spacetime. These Modified Lorentz Transformations (MLTs) are the isometries of the MDR (1), so every particle species presents its personal MLTs, which preserve its specific MDR.
The geometrical structure is defined by the affine and the spinorial connections and here we give their characterization. It is simple to evaluate the explicit forms of the affine connection components, starting from the metric tensor. The affine connection is defined via the Christoffel symbols:
Γ μ ν α = 1 2 g α β μ g β ν + ν g μ β β g μ ν
The following Christoffel symbols are identically equal to zero:
Γ μ 0 0 = Γ 00 i = Γ μ ν i = 0 μ ν
but even the remaining elements are null at the leading order, since they depend on derivatives of the 0 degree homogeneous perturbation function f:
Γ i i 0 = 1 2 0 f ( p ) 0 Γ 0 i 0 = Γ i 0 0 = 1 2 ( 1 + f ( p ) ) 0 f ( p ) 0 Γ i i i = 1 2 ( 1 + f ( p ) ) i f ( p ) 0 Γ j j i = 1 2 ( 1 + f ( p ) ) i f ( p ) 0 i j Γ i j i = Γ j i i = 1 2 ( 1 + f ( p ) ) i f ( p ) 0 i j
The Cartan or spinorial connection is defined as:
ω μ a b = e a ν μ e b ν e a ν μ e b ν
Even the Cartan or spinorial connection is negligible at the leading order, since every element is again proportional to derivatives of the 0 degree homogeneous perturbation function f. Indeed, from the first Cartan structural equations:
d e = e ω
applied to the external forms
e 0 μ = d x μ e i μ = 1 f ( p ) d x μ ,
it follows that the spinorial connection elements are given by:
1 2 ϵ i j k ω i j = 1 2 1 1 f ( p ) ϵ i j k ( i f ( p ) d x j j f ( p ) d x i ) 0
Since the two fundamental connections can be neglected, the total geometric covariant derivative becomes:
D μ v a ν = μ v a ν + Γ μ α ν v a α ω μ ν a v b ν μ v b ν
The resulting geometric structure is an asymptotically flat Finsler pseudo-structure [38,39,40,41,42,43].
Even the Poincaré brackets of the local coordinates are modified:
{ x ˜ μ , x ˜ ν } = { x i e i μ ( p ) , x j e j ν ( p ) } = { x i , e j ν ( p ) } e i μ ( p ) x j + { e i μ ( p ) , x j } x i e j ν ( p ) = 0 { x ˜ μ , p ˜ ν } = { x i e i μ ( p ) , p j e ν j ( p ) } = { x i , p j } e i μ ( p ) e ν j ( p ) + { x i , e ν j ( p ) } e i μ ( p ) p j = δ ν μ + { x i , e ν j ( p ) } e i μ ( p ) p j { p ˜ μ , p ˜ ν } = { p i e μ i ( p ) , p j e ν j ( p ) } = 0
where, in the second equation, the parentheses { x i , e ν j ( p ) } 0 since the vierbein e ν j ( p ) is a function of the derivative p i = x i . Time and space therefore do not commute anymore [37].

2.2. HMSR Standard Model Extension

In the HMSR scenario, as a result, it is possible to amend the formulation of the SM of particle physics following a procedure analogous to that used in the Standard Model Extension (SME) [13]. In order to modify the theory formulation, it is necessary to start from the definition of the amended Dirac matrices, introducing the usual momentum dependence:
Γ μ = e a μ ( p ) γ a Γ 5 = ϵ μ ν α β 4 ! Γ μ Γ ν Γ α Γ β = γ 5
Next, it is possible to modify the Clifford Algebra, starting from the anticommutator of the Γ matrices:
{ Γ μ , Γ ν } = 2 g μ ν ( p ) = 2 e μ a ( p ) η a b e ν b ( p )
Finally it is possible to introduce the modified spinor fields, preserving their plane wave formulation:
ψ + ( x ) = u r ( p ) e i p μ x μ ψ ( x ) = v r ( p ) e i p μ x μ
where the spinors u r ( p ) and v r ( p ) normalization results modified. It can be useful to notice that the modified Dirac equation:
( i Γ μ μ m ) ψ = 0
implies the MDR (1):
( i Γ μ μ + m ) ( i Γ μ μ m ) ψ + = 0 ( Γ μ p μ + m ) ( Γ μ p μ m ) u r ( p ) = 0 1 2 { Γ μ , Γ ν } p μ p ν m 2 u r ( p ) = 0 ( p μ g μ ν p ν m 2 ) u r ( p ) = 0
As a final result, it is possible to formulate the amended particle SM, resorting to the vierbein formalism to project every physical quantity from the modified curved spacetime to the common support flat Minkowski spacetime. For instance, the quantum electrodynamics (QED) Lagrangian can be written as:
L = det [ g ] ψ ¯ ( i Γ μ μ m ) ψ + e det [ g ˜ ] ψ ¯ Γ μ ( p , p ) ψ e ¯ ν μ A ν
where the e ¯ vierbein has been introduced related to the gauge field, and the index μ represents a coordinate of the Minkowski space-time ( T M , η μ ν ) . In order to simplify the theory formulation, the gauge bosons can be supposed to be Lorentz invariant. | det [ g ] | is the term borrowed from curved space-time QFT.
While in the low energy scenario the LIV perturbation can be neglected, in the high energy limit, it is possible to consider incoming and outgoing momenta with approximately the same magnitude, even after interaction. The conserved currents therefore admit a high energy constant limit not dependent on the momentum. The conserved current reduces then to the form:
J μ = e det 1 2 { Γ μ , Γ ν } ψ ¯ Γ μ ψ = e det [ g ] ψ ¯ Γ μ ψ
The modified SM formulation preserves the classical gauge S U ( 3 ) × S U ( 2 ) × U ( 1 ) . The Coleman–Mandula theorem is still valid in a modified formulation that provides an amended symmetry group given by P ( { p } ) G i n t , where P ( { p } ) is the direct product of modified—particle and momentum depending—Poincaré groups:
P { p } = i P ( i ) ( p ( i ) )
and G i n t is the internal symmetry group (in this case, S U ( 3 ) × S U ( 2 ) × U ( 1 ) ).

2.3. HMSR Modified Kinematics

In HMSR framework, the most interesting phenomenological effects caused by the introduction of LIV emerge in the interaction of different particle species, since every particle type physics is modified in a different way by LIV. It is therefore necessary to determine how the reaction invariants—that is, the Mandelstam s, t, and u variables—are modified (Figure 1), considering the p and q momenta as belonging in general to different particle species. It is therefore necessary to generalize the internal product definition for the sum of two or more particle species momenta as:
p + q | p + q = ( p μ e a μ ( p ) + q μ e ˜ a μ ( q ) ) η a b ( p ν e b ν ( p ) + q ν e ˜ b ν ( q ) )
where the vierbeins e a μ ( p ) and e ˜ a μ ( q ) ) —associated respectively with the two different species—are used. Now, one can generalize the definition of the Mandelstam variables s, t, and u resorting to the modified internal product definition. The new internal product can be written introducing the generalized metric:
G = g μ ν ( p ) e a μ ( p ) e ˜ a β ( q ) e ˜ a α ( q ) e a ν ( p ) g ˜ α β ( q )
The internal product (33) can be written as:
p + q | p + q = p q g μ ν ( p ) e a μ ( p ) e ˜ a β ( q ) e ˜ a α ( q ) e a ν ( p ) g ˜ α β ( q ) p q = = p μ g μ ν ( p ) p ν + p μ e a μ e ˜ a β ( q ) q β + q α e ˜ a α ( q ) e a ν ( p ) p ν + q α g ˜ α β ( q ) q β
This generalization presents the feature to be invariant respect to the Modified Lorentz Transfomations MLTs, that is, HMSR presents a covariant formulation, with respect to the MLTs, with the explicit form:
Λ = Λ μ μ 0 0 Λ ˜ α α
Using the MLTs of the two particle species:
p + q | p + q = p q t · G · p q = = Λ p q t · Λ · G · Λ t · Λ p q = Λ ( p + q ) | Λ ( p + q ) =
where Λ · G · Λ t represents the metric evaluated for the two particles momenta given, respectively, by Λ p , Λ ˜ q . It is therefore clear that HMSR does not require the introduction of a preferred reference frame, in contrast with most of the other LIV models.

3. Ultra High Energy Cosmic Rays Propagation and GZK Cut-Off Effect

UHECRs interact with the backgrounds photons, and they are attenuated in a way that depends on their composition and energy; for instance, nuclei can undergo a photo-dissociation process
A + γ ( A 1 ) + n
where A is the atomic number of the considered CR bare nucleus, while protons interacting with the background photons can undergo a pair production:
p + γ p + e + e +
or a photopion process, through a Δ particle resonance:
p + γ Δ p + π 0 p + γ Δ n + π +
The first process is dominant for protons with energy lower than the threshold energy E = 5 × 10 19 eV , whereas the second process becomes the main attenuation mechanism for more energetic protons.
The photopion production is the core mechanism of the GZK phenomenon [6,7] for protons. This process gives an upper limit to the energy of CRs coming from distant sources, since through this effect a proton dissipates energy, but it is not annihilated. If the proton has enough energy, it can undergo the same interaction process again. This way, it becomes possible to evaluate the attenuation length 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: [44]:
1 l p γ = ϵ t h + n ( ϵ ) d ϵ 1 + 1 1 2 s ( 1 v p μ ) σ p γ ( s ) K ( s ) d μ = 1 l p γ = ϵ t h + n ( ϵ ) d ϵ 1 + 1 1 2 s ( 1 v p cos θ ) σ p γ ( s ) K ( s ) d cos θ
where σ p γ ( s ) is the cross section for the proton–photon interaction as a function of the squared energy in the Mandelstam s variable (the squared center of mass energy), μ is the impact parameter μ = cos θ , n ( ϵ ) is the density of background photons per unit volume and photon energy ϵ and ϵ t h is the threshold energy for the interaction. K ( s ) is the reaction inelasticity, defined as the energy fraction available for the production of secondary particles during the reaction. K = ( 1 η ) and η = E out / E in is the elasticity, defined as the energy fraction preserved by the primary particle after the interaction, with E in the incoming particle energy and E out the residual energy. Using the following relations:
s = ( m p + ϵ ) 2 p γ 2 = m p 2 + 2 m p ϵ ϵ = γ ϵ ( 1 v p cos θ )
where the Mandelstam s is computed introducing the photon four momentum ( ϵ , p γ ) defined in the rest frame of the nucleus. In the high energy limit approximation, it is possible to consider the proton velocity v p 1 and with d s = 2 E p ϵ d cos θ the Equation (41) becomes:
1 l p γ = 1 2 γ 2 ϵ t h / 2 γ + d ϵ n ( ϵ ) ϵ 2 ϵ t h ϵ m a x = 2 γ ϵ 1 2 ϵ σ p γ ( ϵ ) K ( ϵ ) d ϵ
where all the standard quantities are defined in the laboratory rest frame and the primed ones are defined in the proton rest frame. Using the fact that the background most relevant to UHECR propagation is the CMB, whose distribution n ( ϵ ) is a black-body spectrum with T 2.7 K , the previous relation (43) can be simplified in the following form, obtaining the inverse of the attenuation length [9]:
1 l p γ = k B T 2 π 2 γ 2 ϵ t h + ϵ σ p γ ( ϵ ) K ( ϵ ) ln 1 e ϵ / 2 k B T γ d ϵ
In the classical Lorentz invariant physics scenario, the inelasticity is given by the relation [44]:
K ( s ) = 1 2 1 m p 2 m π 2 s

4. LIV Introduced Phenomenology in UHECR Propagation

The LIV modifications introduced by the HMSR model are conceived in order to modify the kinematics, avoiding the introduction of exotic particles or reactions. The kinematical perturbations modify the allowed phase space for the reactions and may therefore influence the processes involved in UHECR propagation, such as the GZK effect. The introduction of LIV can indeed affect the photopion production, the core reaction underlying the GZK phenomenon for protons.

4.1. Δ Resonance Creation Constraint

The photopion production passes mainly through a Δ particle creation (38). The process can occur passing through a real Δ particle (dominant process) or a virtual one. The production of the real Δ particle must be preserved in order to foresee only small LIV-induced deviations from the classic GZK process. Indeed, if the real Δ particle production is kinematically forbidden, the photopion production effect and, as a consequence, even the GZK cutoff are strongly suppressed by the LIV introduction. The free energy—that is, the squared root of the Mandelstam s variable—of the proton interacting with the CMB photon must be larger than the rest energy of the Δ resonance, in order to produce one. This threshold energy consideration allows us to pose a first constraint on the LIV magnitude. Using the modified internal product (33), the Mandelstam s variable can be written as a function of the proton and the photon four momenta, respectively, ( E p , p p ) and ( ϵ , p γ ) :
s = ( E p + ϵ ) 2 ( p p [ e p 1 ] + p γ ) 2 m Δ 2 E p 2 p p ( 1 f p ( p p ) ) 2 f p ( p p ) p p + 2 E p ϵ + 2 p p · p γ 1 + 1 2 f p ( p p ) m Δ 2
where the approximation [ e p 1 ] ( p ) = 1 + 1 2 f p ( p ) has been used for the vierbein (14), and f p ( p ) represents the proton LIV perturbation function. From the previous relation (46), it is possible to derive the following relation:
2 f p ( p p ) E p 2 E p ( 4 ϵ + f p ( p p ) ϵ ) + m Δ 2 m p 2 0
This second grade inequality must be satisfied in order to produce a Δ resonance; otherwise, the GZK effect is suppressed by the introduction of LIV. It is therefore simple to derive the following constraint:
f p ( p p ) < Δ M 2 4 ϵ ( Δ M ) 2 8 E γ Δ M 2 4 E γ 2
where Δ M 2 = ( m Δ 2 m p 2 ) . Admitting the average values E γ 7.0 × 10 4 eV , m Δ 1232 MeV , m p 938 MeV , we obtain the disequality:
f p ( p p ) < 6 · 10 23
to guarantee the existence of the GZK effect, comparable with the superior limit 4.5 × 10 23 obtained numerically in [9] assuming that the observed cutoff is entirely due to the propagation. (Scenarios in which the sources themselves have an injection cutoff can be compatible with even stronger LIV in which GZK interactions are disabled altogether [45,46,47].)

4.2. Reduced Phase Space and Modified Inelasticity

The introduction of LIV determines a change of the phase space allowed for the photopion reaction, determining a modification of the inelasticity K (45). The modified inelasticity can be evaluated starting from the modified kinematics. Following the approach of [9,10], in order to simplify the computation, it is useful to define the center of momenta reference frame, where the following relation is satisfied:
p p * + p π * = 0
In this last equation, the vectors are defined in ( T M , η a b ) . The γ C M factor, involved in the change of reference frame from the center of momenta to a generic one, can be computed considering the free energy of the photo-pion production:
s = ( E p * + E π * )
obtaining the equation:
γ C M ( E p * + E π * ) = γ C M s = ( E p + E π ) γ C M = E p + E π s = E t o t s
The free energy necessary for the creation of a photo-pion in the CM frame of reference can be computed, using the CM definition (50) p p * = p π * :
( s E p * ) 2 ( [ e π 1 ] p p * ) 2 = m π 2 ( s 2 s E p * ) + E p * 2 p p * 2 ( 1 + f p ) + p p * 2 f p p p * 2 f π = m π 2
where f p and f π represent respectively the proton and pion LIV correction functions.
From the previous equation follows the equality:
E p * = s + m p 2 + f p p p * 2 m π 2 f π p p * 2 2 s = s + m p 2 m π 2 + f p π p * 2 2 s = F ( s )
where f p π = f p f π is the LIV parameter.
From this relation, it emerges immediately that the residual proton energy after the photopion reaction is increased by the proton correction and reduced by the pion one. This means that the proton correction must be larger than the pion one in order to produce a GZK sphere enlargement. In this work, we make the assumption that every particle has a maximum attainable velocity lower than that of light c. Moreover, we admit that the heavier particles have a bigger LIV induced modification, since the supposed quantum space perturbations have a gravitational origin. This allows us to consider the correction factor of the pion negligible and to pose f p π f p .
Now, it is possible to use the approximations:
p p * E p * = ( 1 k π ( θ ) ) s p π * E π * = k π ( θ ) s
In these equations, s represents the initial free total energy and E p and E π are the final energies of the proton and the pion, respectively.
At this point using the Lorentz changing frame equation, as a result, it is possible to write:
E p = γ C M ( E p * + β cos θ p p * ) E p * = ( 1 k π ( θ ) ) s
using the pion inelasticity.
Approximating the three-momentum magnitude with the energy and the velocity factor β with 1, in the hypothesis of ultra-relativistic particles and substituting γ C M with the value computed in Equation (52), it is possible to obtain the following equation:
( 1 k π ( θ ) ) = 1 s F ( s ) + cos θ F ( s ) 2 m p 2 + 2 f p p 2
from which it is possible to compute the inelasticity in function of the collision angle θ . The equation is solved introducing the energy of the photon ϵ defined in the reference frame where the proton momentum is null p p = 0 , so the s free energy can be written as:
s = ( E p + ϵ ) 2 p γ 2 = 2 m p ϵ + m p 2
We can now average the inelasticity on the interval θ [ 0 , π ] :
k π = 1 π 0 π k π ( θ ) d θ
In Figure 2, Figure 3 and Figure 4, we plot the inelasticity computed in the case of different LIV parameter f p π f p values as a function of the interacting proton E p and photon ϵ energies in the laboratory rest frame. In Figure 2, we report the inelasticity obtained in a Lorentz invariant scenario and given by Equation (45)—LIV parameter posed equal to f p π f p = 0 . In Figure 3, we plot the inelasticity computed in the case of a LIV parameter f p π f p = 9 × 10 23 . Finally, in Figure 4, the inelasticity is plotted for f p π f p = 3 × 10 24 . The LIV caused effects determining a dramatic drop in the inelasticity value (shown in the E p ϵ space where the proton energy E p is defined in the laboratory reference frame and the photon one ϵ is defined in the proton rest frame) for the different values of the LIV parameter. This implies a drastic reduction of the allowed phase space for the creation of a photopion during the interaction of a proton with a CMB photon (54).

5. Simulations on LIV Modified Propagation

Making use of the modified shape of the inelasticity K function derived from the theory, we can now compute the modifications experienced by an UHECR proton for the GZK opacity horizon and the attenuation length, assuming that their production rate is constant in space and time. To define the GZK opacity sphere, we consider the region in which 50 % of the protons that reach Earth with a residual energy higher (or equal to) the threshold of the photopion production process, which is 5 × 10 19 eV , has originated.
In order to evaluate the modifications introduced by the LIV theory, we simulated the UHECR protons’ free propagation using a specifically modified version of the simulation program SimProp [20]. In the modified version of the software, we substituted the classical physics predicted inelasticity with the amended one given by an analytic function of the energy of the proton, the CMB photon interacting energy defined in the proton rest frame and of the LIV parameter. This modification determines that the GZK produced photopion energy is quickly reduced in the case the energies of the interacting particles belong to a determined range in the phase space, since the value of the inelasticity in Figure 2, Figure 3 and Figure 4 drops quickly to negligible values.
The LIV-caused GZK sphere modifications are evaluated simulating the propagation of 1000 protons for various LIV parameters magnitudes with various fixed energies. The results of the simulations are then used to plot the GZK opacity sphere radius as a function of the energy and the LIV parameters (see Figure 5). Every plot line is obtained fixing the LIV parameter difference magnitude from a minimum value of 8 × 10 24 to a maximum of 4 × 10 22 and is compared with the result obtained for null LIV. For fixed LIV parameters, the 1000 particles are initially simulated with an acceleration energy E = 5 × 10 19 eV . The simulations are then repeated increasing the production energy by steps of 5 × 10 18 eV reaching a maximum value of E = 2 × 10 20 eV , in order to obtain every plot line. We made the conservative assumption that, in a CR group with varying production energy, the dominant component is the lowest one. Indeed, the energy loss rate increases steeply with increasing CR energy, so it is possible to simplify the analysis considering every detected proton as generated with the inferior energy value for every CR group. The particles are simulated assuming a uniform distribution inside a sphere of radius ∼1270 Mpc centered on Earth that is inside a portion of Universe included between redshift parameters z = 0 and z = 0.25 .
From Figure 5, one can clearly notice how the opacity sphere for UHECR protons is modified by the LIV introduction; indeed, the GZK proton horizon is LIV parameter dependent. Figure 5 has visible inflections points where the steepness of the GZK radius plot increases as a function of the energy. These inflection points depend on the LIV parameter magnitude; indeed, they are present at lower energies for increasing LIV values. This behaviour is correlated with analogous inflections that are present in the plot of the attenuation length (Figure 6) as a function of energy and LIV as explained below.
Following the same simulation procedure, we can also plot the attenuation length for the reaction. The results of the simulations for different LIV parameters are shown in Figure 6. In the LIV-less scenario, the attenuation length decreases with increasing energy values. In the presence of LIV perturbations, it is clear that the simulated energy loss length undergoes a recovery that becomes more and more important for increasing values of the LIV parameter. Indeed, the attenuation length at first decreases with increasing energy; then, there is an inflection point and, after this, the increases for increasing energy. This effect is caused by the average energy lost in every photopion production process that is reduced by LIV. Indeed, up to LIV parameters of 10 22 , the variation from the Lorentz invariant (LI) case is noticeable starting from production energy values of ∼ 10 20 eV , while, for LIV parameters at least equal to 2 × 10 22 , the modification is appreciable starting already from ∼ 6 × 10 19 eV . The inflection points depends therefore on the LIV magnitude as in the case of GZK-cutoff plot. Note that we are neglecting the effects of LIV on electron–positron pair production interactions, which for large values of the LIV parameter would presumably further increase the energy loss lengths.
As a result, one can state that the GZK sphere enlargement and the attenuation length increase imply that the introduction of LIV determines a reduction of the Universe opacity for UHECR propagation.

6. Conclusions

Lorentz covariance underlies nowadays our knowledge of physics; nevertheless, small departures from this fundamental symmetry are supposed as possible residual evidence of the phenomenology induced by quantum gravity. Our study is conducted in the framework of HMSR theory a model that foresees a LIV minimal extension of the SM of particle physics preserving the covariance of the formulation in a context able to predict a rich phenomenology. HMSR theory foresees a minimal extension of the SM of particle physics that preserves the usual gauge invariance S U ( 3 ) × S U ( 2 ) × U ( 1 ) and does not introduce exotic particles or reactions. Moreover, the theory preserves the covariance in an amended formulation that does not require the introduction of a preferred reference frame. This last aspect is of great experimental advantage since all the obtained observations concerning UHECR can be collected without discriminations related to the orientation of the observatory respect to a privileged reference frame. This model can therefore be used to conduct studies of correlation between the anisotropy of the arrival direction of UHECR with the distribution of candidate sources thanks to the preservation of the spacetime isotropy.
Our study is conducted by simulating the propagation of UHECR protons in the presence of LIV using an ad hoc modified version of the SimProp software. The obtained result foresees an enlargement of the GZK opacity sphere and of the correlated attenuation lengths. This predictions are in accordance with other works [9,10,11,48,49], but, with the difference of being obtained in a covariant preserving framework. The accordance of these results with previous works indicates a great accuracy of the SimProp modifications introduced to take into account the LIV generated effects on UHECR protons and guarantees a solid basis for the research framework introduced here.
The theoretical predictions about the Δ resonance creation suggest a LIV parameter value of ∼ 6 × 10 23 (49), similar to values found in previous works [9,50] by fitting the observed UHECR spectrum assuming a pure proton injection with no cutoff at sources. To obtain concrete results from the comparison with experimental data about UHECR, it may be necessary to improve our study. The new physics introduced by LIV leaves room for an analysis refinement studying a more realistic scenario. Future developments could implement the study of heavier UHECRs (for instance, the intermediate component), in order to follow the indications given by the Pierre Auger Collaboration that privileges a mixed UHECR composition [51]. Indeed, LIV effects are expected to be larger for intermediate nuclei due to the nature of the kinematical modifications, which could affect their propagation in a more relevant way than that of the lighter or heavier UHECR components. We made an excessive conservative hypothesis considering only protons that is the lighter UHECR component, but, to support this first approximation, it is important to underline that all the results obtained for this scenario are still valid for the more realistic one, for which they are expected to be even more important. The results reported in this paper can be used to compare predictions with experimental observations. In particular, one may try and compute how observables such as spectrum, composition, and secondary photon flux are expected to change at different LIV intensities, as was done in [45,46,47,52,53,54]. Moreover, since the HMSR has the peculiarity of preserving spacetime isotropy, it is particularly fit to make predictions on how UHECR anisotropy evolves with LIV in different sources scenarios, as we plan to do this in a forthcoming work.

Author Contributions

Writing—original draft: M.D.C.T., L.C., A.d.M., A.M. and L.M.; review and editing: M.D.C.T. All authors contributed equally to this work. All authors have read and agreed to the published version of the manuscript.

Funding

This work was partially supported by a post doctoral fellowship financed by the Fondazione Fratelli Giuseppe Vitaliano, Tullio e Mario Confalonieri, Milano—Italy.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
UHECRUltra High Energy Ccosmic Ray
CRCosmic Rays
LIVLorentz Invariance Violation
HMSRHomogeneously Modified Special Relativity
SMStandard Model
SMEStandard Model extension
VSRVery Special Relativity
DRDispersion Relation
MDRModified Dispersion Relation
MLTModified Lorentz Transformation

Appendix A. Derivative of the 0 Degree Homogeneous Function

One fundamental feature of HMSR is the possibility to establish a direct relation between the momentum and the coordinate space, at least at the leading order. This peculiar characteristic follows from the possibility to neglect the derivatives of the co-metric by momentum. The negligibility of these derivatives implies that the geometric structure is asymptotically flat (that is approximately flat but not trivial) and the total covariant derivative can be approximated with the standard partial one. All the previous considerations follow from the homogeneity of the perturbation function f in the MDR (1). The homogeneity of the perturbation function f is fundamental to neglect the derivatives by the momentum, which is proportional to terms with the form:
p μ f | p | E
Finally, taking the Taylor expansion series of the f function and the high energy limit, it follows the relation:
p j f ( p ) = p j k α k | p | k E k p j k α k | p | k ( | p | 2 + m 2 ) k = = k α k k | p | k 2 p j ( | p | 2 + m 2 ) k α k k | p | k p j ( | p | 2 + m 2 ) k + 2 0
where the approximate equivalence E | p | 2 + m 2 has been used.

References

  1. Pierre Auger Collaboration. Observation of a Large-scale Anisotropy in the Arrival Directions of Cosmic Rays above 8 × 1018 eV. Science 2017, 357, 1266–1270. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Pierre Auger Collaboration. Constraints on the origin of cosmic rays above 1018 eV from large scale anisotropy searches in data of the Pierre Auger Observatory. Astrophys. J. Lett. 2012, 762, L13. [Google Scholar]
  3. Tinyakov, P.G.; Urban, F.R.; Ivanov, D.; Thomson, G.B.; Tirone, A.H. A signature of EeV protons of Galactic origin. Mon. Not. R. Astron. Soc. 2016, 460, 3479–3487. [Google Scholar] [CrossRef] [Green Version]
  4. Abbasi, R.U.; Abe, M.; Abu-Zayyad, T.; Allen, M.; Azuma, R.; Barcikowski, E.; Belz, J.W.; Bergman, D.R.; Blake, S.A.; Cady, R.; et al. Search for EeV Protons of Galactic Origin. Astropart. Phys. 2017, 86, 21–26. [Google Scholar] [CrossRef] [Green Version]
  5. Blasi, P. Origin of very high- and ultra-high-energy cosmic rays. C. R. Physique 2014, 15, 329–338. [Google Scholar] [CrossRef] [Green Version]
  6. Greisen, K. End to the cosmic ray spectrum? Phys. Rev. Lett. 1966, 16, 748–750. [Google Scholar] [CrossRef]
  7. Zatsepin, G.; Kuzmin, V. Upper limit of the spectrum of cosmic rays. JETP Lett. 1966, 4, 78–80. [Google Scholar]
  8. Coleman, S.R.; Glashow, S.L. High-energy tests of Lorentz invariance. Phys. Rev. D 1999, 59, 116008. [Google Scholar] [CrossRef] [Green Version]
  9. Stecker, F.W.; Scully, S.T. Searching for New Physics with Ultrahigh Energy Cosmic Rays. New J. Phys. 2009, 11, 085003. [Google Scholar] [CrossRef]
  10. Scully, S.; Stecker, F. Lorentz Invariance Violation and the Observed Spectrum of Ultrahigh Energy Cosmic Rays. Astropart. Phys. 2009, 31, 220–225. [Google Scholar] [CrossRef] [Green Version]
  11. Torri, M.D.C.; Bertini, S.; Giammarchi, M.; Miramonti, L. Lorentz Invariance Violation effects on UHECR propagation: A geometrized approach. J. High Energy Astrophys. 2018, 18, 5–14. [Google Scholar] [CrossRef]
  12. Torri, M.D.C. Lorentz Invariance Violation Effects on Ultra High Energy Cosmic Rays Propagation, a Geometrical Approach. Ph.D. Thesis, Milan University, Milan, Italy, 2019. [Google Scholar]
  13. Colladay, D.; Kostelecký, V.A. Lorentz violating extension of the standard model. Phys. Rev. D 1998, 58, 116002. [Google Scholar] [CrossRef] [Green Version]
  14. Cohen, A.G.; Glashow, S.L. Very special relativity. Phys. Rev. Lett. 2006, 97, 021601. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Amelino-Camelia, G. Doubly special relativity. Nature 2002, 418, 34–35. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Amelino-Camelia, G. Doubly special relativity: First, results and key open problems. Int. J. Mod. Phys. D 2002, 11, 1643–1669. [Google Scholar] [CrossRef]
  17. Amelino-Camelia, G.; Freidel, L.; Kowalski-Glikman, J.; Smolin, L. The principle of relative locality. Phys. Rev. D 2011, 84, 084010. [Google Scholar] [CrossRef] [Green Version]
  18. Amelino-Camelia, G.; Bianco, S.; Rosati, G. Planck-Scale-Deformed Relativistic Symmetries and Diffeomorphisms on Momentum Space. Phys. Rev. D 2020, 101, 026018. [Google Scholar] [CrossRef] [Green Version]
  19. Torri, M.; Antonelli, V.; Miramonti, L. Homogeneously Modified Special relativity (HMSR). Eur. Phys. J. C 2019, 79, 808. [Google Scholar] [CrossRef] [Green Version]
  20. Aloisio, R.; Boncioli, D.; Grillo, A.; Petrera, S.; Salamida, F. SimProp: A Simulation Code for Ultra High Energy Cosmic Ray Propagation. J. Cosmol. Astropart. Phys. 2012, 10, 7. [Google Scholar] [CrossRef] [Green Version]
  21. Magueijo, J.; Smolin, L. Gravity’s rainbow. Class. Quant. Grav. 2004, 21, 1725. [Google Scholar] [CrossRef] [Green Version]
  22. Magueijo, J.; Smolin, L. Lorentz invariance with an invariant energy scale. Phys. Rev. Lett. 2002, 88, 190403. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Pfeifer, C.; Wohlfarth, M.N.R. Finsler geometric extension of Einstein gravity. Phys. Rev. D 2012, 85, 064009. [Google Scholar] [CrossRef] [Green Version]
  24. Hohmann, M.; Pfeifer, C.; Voicu, N. Finsler gravity action from variational completion. Phys. Rev. D 2019, 100, 064035. [Google Scholar] [CrossRef] [Green Version]
  25. Pfeifer, C. Finsler spacetime geometry in Physics. Int. J. Geom. Meth. Mod. Phys. 2019, 16, 1941004. [Google Scholar] [CrossRef] [Green Version]
  26. Javaloyes, M.A.; Sánchez, M. On the definition and examples of cones and Finsler spacetimes. arXiv 2018, arXiv:1805.06978. [Google Scholar] [CrossRef] [Green Version]
  27. Bernal, A.; Javaloyes, M.Á.; Sánchez, M. Foundations of Finsler spacetimes from the Observers’ Viewpoint. Universe 2020, 6, 55. [Google Scholar] [CrossRef] [Green Version]
  28. Greenberg, O. Why is CPT fundamental? Found. Phys. 2006, 36, 1535–1553. [Google Scholar] [CrossRef] [Green Version]
  29. Greenberg, O. CPT violation implies violation of Lorentz invariance. Phys. Rev. Lett. 2002, 89, 231602. [Google Scholar] [CrossRef] [Green Version]
  30. Antonelli, V.; Miramonti, L.; Torri, M.D.C. Phenomenological Effects of CPT and Lorentz Invariance Violation in Particle and Astroparticle Physics. Symmetry 2020, 12, 1821. [Google Scholar] [CrossRef]
  31. Chaichian, M.; Dolgov, A.D.; Novikov, V.A.; Tureanu, A. CPT Violation Does Not Lead to Violation of Lorentz Invariance and Vice Versa. Phys. Lett. B 2011, 699, 177–180. [Google Scholar] [CrossRef] [Green Version]
  32. Tureanu, A. CPT and Lorentz Invariance: Their Relation and Violation. J. Phys. Conf. Ser. 2013, 474, 2031. [Google Scholar] [CrossRef] [Green Version]
  33. Chaichian, M.; Fujikawa, K.; Tureanu, A. Electromagnetic interaction in theory with Lorentz invariant CPT violation. Phys. Lett. B 2013, 718, 1500–1504. [Google Scholar] [CrossRef] [Green Version]
  34. Duetsch, M.; Gracia-Bondia, J.M. On the assertion that PCT violation implies Lorentz non-invariance. Phys. Lett. B 2012, 711, 428–433. [Google Scholar] [CrossRef] [Green Version]
  35. Greenberg, O. Remarks on a Challenge to the Relation between CPT and Lorentz Violation. arXiv 2011, arXiv:1105.0927. [Google Scholar]
  36. Antonelli, V.; Miramonti, L.; Torri, M.D.C. Neutrino oscillations and Lorentz Invariance Violation in a Finslerian Geometrical model. Eur. Phys. J. C 2018, 78, 667. [Google Scholar] [CrossRef] [Green Version]
  37. Torri, M.D.C. Neutrino Oscillations and Lorentz Invariance Violation. Universe 2020, 6, 37. [Google Scholar] [CrossRef] [Green Version]
  38. Kostelecky, A. Riemann-Finsler geometry and Lorentz-violating kinematics. Phys. Lett. B 2011, 701, 137. [Google Scholar] [CrossRef] [Green Version]
  39. Girelli, F.; Liberati, S.; Sindoni, L. Planck-scale modified dispersion relations and Finsler geometry. Phys. Rev. D 2007, 75, 064015. [Google Scholar] [CrossRef] [Green Version]
  40. Edwards, B.R.; Kostelecky, V.A. Riemann-Finsler geometry and Lorentz-violating scalar fields. Phys. Lett. B 2018, 786, 319. [Google Scholar] [CrossRef]
  41. Lämmerzahl, C.; Perlick, V. Finsler geometry as a model for relativistic gravity. Int. J. Geom. Meth. Mod. Phys. 2018, 15, 1850166. [Google Scholar] [CrossRef] [Green Version]
  42. Lämmerzahl, C.; Perlick, V. Axiomatic formulations of modied gravity theories with nonlinear dispersion relations and Finsler Lagrange Hamilton geometry. Eur. Phys. J. C 2018, 78, 969. [Google Scholar]
  43. Schreck, M. Classical Lagrangians and Finsler structures for the nonminimal fermion sector of the Standard-Model Extension. Phys. Rev. D 2016, 93, 105017. [Google Scholar] [CrossRef] [Green Version]
  44. Stecker, F.W. Photodisintegration of ultrahigh-energy cosmic rays by the universal radiation field. Phys. Rev. 1969, 180, 1264–1266. [Google Scholar] [CrossRef]
  45. Boncioli, D.; di Matteo, A.; Salamida, F.; Aloisio, R.; Blasi, P.; Ghia, P.L.; Grillo, A.; Petrera, S.; Pierog, T. Future prospects of testing Lorentz invariance with UHECRs. arXiv 2015, arXiv:1509.01046. [Google Scholar]
  46. Boncioli, D.; Pierre Auger Collaboration. Probing Lorentz symmetry with the Pierre Auger Observatory. In Proceedings of the ICRC17, Busan, Korea, 10–20 July 2017; Volume 561. [Google Scholar]
  47. Lang, R.G.; Pierre Auger Collaboration. Testing Lorentz Invariance Violation at the Pierre Auger Observatory. In Proceedings of the ICRC19, Madison, WI, USA, 24 July–1 August 2019; p. 327. [Google Scholar]
  48. Bietenholz, W. Cosmic Rays and the Search for a Lorentz Invariance Violation. Phys. Rep. 2011, 505, 145–185. [Google Scholar] [CrossRef] [Green Version]
  49. Mattingly, D. Modern tests of Lorentz invariance. Living Rev. Rel. 2005, 8, 5. [Google Scholar] [CrossRef] [Green Version]
  50. Kostelecky, V.A.; Russell, N. Data Tables for Lorentz and CPT Violation. Rev. Mod. Phys. 2011, 83, 11–31. [Google Scholar] [CrossRef] [Green Version]
  51. Pierre Auger Collaboration. Depth of maximum of air-shower profiles at the Pierre Auger Observatory. II. Composition implications. Phys. Rev. D 2014, 90, 122006. [Google Scholar] [CrossRef] [Green Version]
  52. Diaz, J.S.; Klinkhamer, F.R.; Risse, M. Changes in extensive air showers from isotropic Lorentz violation in the photon sector. Phys. Rev. D 2016, 94, 085025. [Google Scholar] [CrossRef] [Green Version]
  53. Klinkhamer, F.R.; Niechciol, M.; Risse, M. Improved bound on isotropic Lorentz violation in the photon sector from extensive air showers. Phys. Rev. D 2017, 96, 116011. [Google Scholar] [CrossRef] [Green Version]
  54. Lang, R.G.; Martínez-Huerta, H.; de Souza, V. Limits on the Lorentz Invariance Violation from UHECR astrophysics. Astrophys. J. 2018, 853, 23. [Google Scholar] [CrossRef]
Figure 1. Feynman diagrams related to Mandelstam variables.
Figure 1. Feynman diagrams related to Mandelstam variables.
Symmetry 12 01961 g001
Figure 2. Inelasticity obtained in the case of LIV parameter f p π f p = 0 as a function of the proton energy E p and of the photon energy ϵ defined in the proton rest frame.
Figure 2. Inelasticity obtained in the case of LIV parameter f p π f p = 0 as a function of the proton energy E p and of the photon energy ϵ defined in the proton rest frame.
Symmetry 12 01961 g002aSymmetry 12 01961 g002b
Figure 3. Inelasticity obtained in the case of LIV parameter f p π f p = 9 × 10 23 as a function of the proton energy E p and of the photon energy ϵ defined in the proton rest frame.
Figure 3. Inelasticity obtained in the case of LIV parameter f p π f p = 9 × 10 23 as a function of the proton energy E p and of the photon energy ϵ defined in the proton rest frame.
Symmetry 12 01961 g003aSymmetry 12 01961 g003b
Figure 4. Inelasticity obtained in the case of LIV parameter f p π f p = 3 × 10 24 as a function of the proton energy E p and of the photon energy ϵ defined in the proton rest frame.
Figure 4. Inelasticity obtained in the case of LIV parameter f p π f p = 3 × 10 24 as a function of the proton energy E p and of the photon energy ϵ defined in the proton rest frame.
Symmetry 12 01961 g004aSymmetry 12 01961 g004b
Figure 5. GZK sphere radius as a function of production energy E simulated for 1000 protons generated randomly inside of a sphere with a radius of ∼ 1270 Mpc , centered on Earth. The simulations for ten different values of the LIV parameter are shown in the energy range 0 eV 2 × 10 20 eV .
Figure 5. GZK sphere radius as a function of production energy E simulated for 1000 protons generated randomly inside of a sphere with a radius of ∼ 1270 Mpc , centered on Earth. The simulations for ten different values of the LIV parameter are shown in the energy range 0 eV 2 × 10 20 eV .
Symmetry 12 01961 g005
Figure 6. Attenuation length as a function of energy E, plotted for ten different values of the LIV parameter.
Figure 6. Attenuation length as a function of energy E, plotted for ten different values of the LIV parameter.
Symmetry 12 01961 g006
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Torri, M.D.C.; Caccianiga, L.; di Matteo, A.; Maino, A.; Miramonti, L. Predictions of Ultra-High Energy Cosmic Ray Propagation in the Context of Homogeneously Modified Special Relativity. Symmetry 2020, 12, 1961. https://doi.org/10.3390/sym12121961

AMA Style

Torri MDC, Caccianiga L, di Matteo A, Maino A, Miramonti L. Predictions of Ultra-High Energy Cosmic Ray Propagation in the Context of Homogeneously Modified Special Relativity. Symmetry. 2020; 12(12):1961. https://doi.org/10.3390/sym12121961

Chicago/Turabian Style

Torri, Marco Danilo Claudio, Lorenzo Caccianiga, Armando di Matteo, Andrea Maino, and Lino Miramonti. 2020. "Predictions of Ultra-High Energy Cosmic Ray Propagation in the Context of Homogeneously Modified Special Relativity" Symmetry 12, no. 12: 1961. https://doi.org/10.3390/sym12121961

APA Style

Torri, M. D. C., Caccianiga, L., di Matteo, A., Maino, A., & Miramonti, L. (2020). Predictions of Ultra-High Energy Cosmic Ray Propagation in the Context of Homogeneously Modified Special Relativity. Symmetry, 12(12), 1961. https://doi.org/10.3390/sym12121961

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop