Quantum gravity phenomenology induced in the propagation of UHECR, a kinematical solution in Finsler and generalized Finsler spacetime

It is well-known that the Universe is opaque to the propagation of Ultra-High-Energy Cosmic Rays (UHECRs) since these particles dissipate energy during their propagation interacting with the Cosmic Microwave Background (CMB) mainly in the so-called GZK cut-off phenomenon. Some experimental evidence seems to hint at the possibility of a dilation of the GZK predicted opacity sphere. It is well-known that kinematical perturbations caused by supposed quantum gravity (QG) effects can modify the foreseen GZK opacity horizon. The introduction of Lorentz Invariance Violation (LIV) can indeed reduce, in some cases making negligible, the CMB-UHECRs interaction probability. In this work we explore the effects induced by modified kinematics in the UHECRs phenomenology from the QG perspective. We explore the possibility of a geometrical description of the massive fermions interaction with the supposed quantum structure of spacetime in order to introduce a Lorentz covariance modifification. The kinematics is amended modifying the Dispersion Relations (DRs) of free particles in the context of a covariance-preserving framework. This spacetime description requires a more general geometry than the usual Riemannian one, indicating for instance the Finsler construction and the related generalized Finsler spacetime as ideal candidates. Finally we investigate the correlation between the magnitude of Lorentz covariance modification and the attenuation length of the photopion production process related to the GZK cut-off, demonstrating that the predicted opacity horizon can be dilated even in the context of a theory that does not require any privileged reference frame.

geometrical description of the massive fermions interaction with the supposed quantum structure of spacetime in order to introduce a Lorentz covariance modification. The kinematics is amended modifying the Dispersion Relations (DRs) of free particles in the context of a covariance-preserving framework. This spacetime description requires a more general geometry than the usual Riemannian one, indicating for instance the Finsler construction and the related generalized Finsler spacetime as ideal candidates. Finally we investigate the correlation

Introduction
The Universe is opaque to the propagation of radiation in particular of its most energetic component, that is Cosmic Rays (CR). This part of the extraterrestrial radiation is composed of charged particles, protons and bare nuclei, most likely of extragalactic origin, with energies above E ∼ 1 EeV. Since the Universe opacity is more pronounced at the highest energies, in this work we deal with Ultra High Energy Cosmic Rays (UHECRs) that is the most energetic fraction of CR with energies larger than ∼ 5 × 10 19 eV. UHECRs during their propagation may interact with the Cosmic Microwave Background (CMB) and may be attenuated in a way that depends on their energy and nature. Protons dissipate energy via pair production or photopion creation, instead bare nuclei can undergo photodissociation processes. This means that the UHECRs free propagation path is finite and if they propagate for long enough distances they can be detected on Earth only under a determined energy threshold. This phenomenon poses an upper limit on the detectable UHECRs energy and is named GZK cut-off from the name of the physicists Greisen, Zatsepin and Kuzmin [1,2].
Astroparticle physics can be a useful framework to conduct searches for possible departures from the Lorentz covariance since in some theoretical models the Quantum Gravity (QG) effects are supposed to be more visible in the highenergy limit as Lorentz Invariance Violation (LIV) effects. Currently there is no a definitive theory that unifies quantum physics with gravity. The greatest challenge in formulating such a unified theory is the actual impossibility of obtaining the energies needed to probe the Planck scale realm. It is commonly believed that the Planck energy E P = √ c 5 G −1 1.22 × 10 19 GeV represents the threshold that separates the classical formulation of physics, which is General Relativity (GR), from the quantum realm, represented by Quantum Field Theory (QFT). Nevertheless some QG signatures may become manifest in a low-energy regime as residual effects and these perturbations of standard physics may give access to the possible phenomenology induced by QG.
The physics correlated with the UHECR propagation can indeed be exploited to investigate new phenomena, such as the supposed quantum structure of the spacetime [3,5,4,6,7,8]. UHECRs propagate for cosmic distances and reach extremely high energies, therefore they can open a window on the Planck scale physics.
In this work we investigate the possibility to detect the supposed QG signatures as departures from the Lorentz covariance in the UHECRs sector in the context of a theoretical framework that does not require the introduction of a privileged reference frame. There are many theoretical models that investigate the possibility of modifications as QG-induced perturbations [9,10,3,11,12,13,14,15], but we conduct our investigation in the context of Homogeneously Modified Special Relativity (HMSR) [16]. In this theoretical framework the covariance is modified and not broken as in Doubly Special Relativity (DSR) [12,13,14,15], hence the kinematical symmetry is again valid in an amended formulation. In this context a minimal extension of the Standard Model (SM) of particle physics can be formulated preserving the modified covariance. The QG perturbations are introduced in the free particle kinematics modifying the dispersion relations.
This modification is motivated by the idea that the effects induced by the quantum structure of the background can be geometrized requiring a more general structure to set the theory than the usual Riemann geometry, the Finsler one.
The phenomenological analysis is limited to the lightest component of the UHE-CRs, that is protons. This component is interesting since it is the least deflected one by magnetic fields during propagation so proton trajectories are more straightforward than those of the heavier CRs. The introduction of QG kinematical perturbations can modify the allowed phase space for the interaction with the CMB and the pion creation process. The phenomenological effects are evaluated via numerical simulations conducted using an ad hoc version of the software SimProp [17], modified in order to include the QG kinematical perturbations [7]. The introduction of Lorentz covariance modifications induce an opacity horizon enlargement [6,7], as foreseen even in the case of pure LIV [3,5,4].
This work is structured as follows: first we introduce the theoretical framework motivated by the QG perspective, geometrizing the interaction of a free particle with the background. Then we introduce a minimal extension of the

Kinematical modifications in an isotropy preserving scenario
In this work we explore the possibility of investigating the presumed QG effects perturbing the free particle kinematics and we obtain a minimal extension of the particle SM in a framework that preserves an amended covariance formulation.

Modified Dispersion Relations
The QG-caused perturbations are introduced only in the free particle kinematics in order to prevent the introduction of any exotic particle or reaction.
In this work as a reasonable physical hypothesis only the dispersion relations of massive fermions are modified, because of the gravitational nature of these supposed effects. The QG perturbations are introduced modifying the Dispersion Relations (DRs) of massive particles, as done in the great part of theories [9,3,11,12,13,14,15,18,19] that tackle the Special Relativity (SR) modifications from a QG perspective. In this work every particle species is supposed to have its personal Modified Dispersion Relation (MDR): with peculiar functions h and k depending on the particle species. These functions have to satisfy the perturbation condition: with derivatives of first and second order limited with at least the same magnitude. One may consider, for instance, the linear function α(| p|/E) or the translated exponential exp(α| p|/E) − 1 with α a strongly limited coefficient that encodes the perturbation character of the selected function. In the following all the terms with comparable magnitude will be considered first-order The perturbation functions h and k are 0-degree homogeneous and depend on the 3-momentum magnitude | p| in order to be rotationally invariant, as conjectured in [16].  [16]: In Finsler geometry the norm can define a local metric via the equation: which requires that det g µν = 0 and g µν ∈ R ∀µ, ν. In this way it is possible to reobtain a vector norm: As in Riemann geometry it is possible to define a duality relation between vectors and dual forms. The definition of the Legendre transform can be posed using the metric and the resulting bijection is given by: The previous definition is valid for a Finsler structure with a definite signature of the metric, but it must be generalized in order to deal with the physical spacetime, which exhibits a Minkowski-type underlying metric of indefinite-signature.
A systematical introduction to this issue and the so-called pseudo-Finsler geometry can be found in [22,23,24,25,26,27]. Here we follow the definition given in [23,25] and pose a pseudo-Finsler geometric structure modifying the previous definition and requiring that: and T * (T M \ 0) structures, as in [28].
Considering T * (T M \ 0) one can introduce the vertical distribution generated by the derivative ∂/∂v as: Given the vertical distribution, it is possible to identify a complementary structure, the horizontal distribution, named nonlinear connection, for which the Whitney decomposition is valid [29]: The nonlinear connection is a collection of homogeneous functions of degree 1, locally defined on the manifold, such that: In this work we follow [28] and as nonlinear connection we choose the usual General Relativity (GR) form: where Γ α µν represents the GR affine connection. The choice for a nonlinear connection is not uniquely defined, but this particular formulation is the most convenient one in this context. Analogously, the cotangent space is spanned by the differential basis defined by: In the context of Finsler geometry it is possible to generalize the definition of the affine connection via the general Christoffel symbols: and the geodesic equations become: The study conducted in this work can be generalized with these prescriptions in a curved spacetime. In the following, we will consider a geometry with no matter induced curvature for the investigation of the phenomenology introduced in the propagation of UHECRs. Hence, the nonlinear connection eq. (12) is zero and the derivative eq. (11) reduces to the usual partial derivative.

Finsler geometry and HMSR
The MDR (1) can be interpreted as the Finsler pseudo-norm characterizing the spacetime (momentum space) structure. In HMSR the standard physics is modified in a CPT-even scenario, the MDRs, indeed, do not exhibit any dependence on helicities, for instance, and are supposed to be equal for particles and the related antiparticles. A well-known result about fundamental physics symmetry states that in a CPT-odd scenario the Lorentz symmetry as usually formulated must be violated, whereas the opposite statement is not automatically true and therefore it is possible to violate or modify covariance, while preserving the CPT symmetry [30,31]. A review on the implications of the CPT symmetry in particle and astroparticle physics and its relation to the Lorentz/Poincaré symmetry is presented, for instance, in [32].
Promoting the MDR eq. (1) to the role of a Finsler norm, it is possible to obtain the related metric of the momentum space as the Hessian of the squared norm: Since the function F 2 is 2-degree homogeneous, the computation produces a 0-homogeneous tensor that can be written as: with a diagonal matrix D: The matrix D is the associated generalized Hamilton's space metric that effectively contributes to the computation of the MDR eq. (1). Indeed, by direct calculation it is straightforward to demonstrate the relation: The complete Finsler metric tensor is written including the matrix A eq. (17) and the metric tensor defined in this way satisfies the Finsler requirement of a totally symmetric associated Cartan's tensor: The matrix A has both diagonal and nondiagonal entries and thanks to the explicit form acquired by the metric tensor g it is simple to find out that this matrix gives no contribution to the norm of a covector: p µ A µν p ν = 0 eq. (19). From the Finsler co-metric one can derive the metric associated with coordinate 2 See Appendix A space as the inverse one, via the defining equation: obtaining the explicit form of the coordinate space metric: and the diagonal matrix D * : where O(2) indicates a second-order perturbation.

Hamilton and Lagrangian Finsler geometry
The correlation between coordinate and momentum space must be further explored and it can be done starting from the action: where the Lagrangian L is defined using a Lagrangian multiplier λ to pose the constraint that the MDR, F 2 (p) = m 2 , must be satisfied [33]. The action S is obtained integrating the Lagrangian L with respect to the proper time τ .
The equations of motion can be obtained varying the Lagrangian L. Hence, one can obtain the partial derivatives with respect to the velocityẋ µ and the Lagrange multiplier λ: The variation calculated with respect to the momentum p µ gives: A simple application of the Euler's theorem on homogeneous functions implies that the second term on the right-hand side of the previous relation is zero thanks to the 0-homogeneity of the metric tensor. Therefore, the velocity can be written in the simple form:ẋ Inverting the previous relation, the momentum can be expressed as: Substituting the previous relation in the action eq. (24) it is possible to compute the Hamiltonian: which acquires this simple explicit form since the second term on the right-hand side in eq. (28) is orthogonal to the momentum p µ .
The momentum can be expressed as a function of the velocity solving eq. (28), for instance, perturbatively. As a consequence, the coordinate space metric g µν (p), defined as the inverse of the momentum space metric g µν (p), can be written as a function of coordinate x and velocityẋ: The Lagrangian can be computed substituting eq. (26) in the action eq. (24): where O(2) is a second-order perturbation with respect to the scale fixed by the functions h and k (2) again thanks to the orthogonality of p µ and the second term on the right-hand side of eq. (28). The variation of the previous relation with respect to the parameter λ gives: and then the Lagrangian can be written as: The relation between the Lagrangian and the Hamiltonian formulation has been established and it is possible to compute the norm for the velocityẋ. This function is associated with the momentum space norm F which is defined via the MDR eq. (1). The coordinate space norm can be determined starting from the metric eq. (22): As a final remark it is possible to state that the resulting structure is a Finsler geometry.

HMSR generalized Finsler spacetime
In the following we will consider a generalized Finsler geometry model to set the stage of the HMSR formulation. In this context we will obtain an extension of the Standard Model (SM) of particle physics and will conduct the computations related with the GZK cutoff phenomenon.
The generalized Finsler geometry is a less restrictive structure which does not require a totally symmetric Cartan's tensor associated to the metric eq. (20), therefore, the construction of the geometry is easier. In this context the space-time and the momentum metrics eqs. (17) and (22) can be simplified taking into account only the diagonal parts eqs. (18) and (23): In this way the coordinate and momentum metrics are reduced to the parts that really contribute in evaluating the squared norm of a vector or a covector respectively.
Starting from the defining equations: it is possible to compute the generalized associated vierbein, which can be written as: It is simple to determine the perturbative order O(1) of part of these matrix entries using the approximation The form of the vierbein eq. (38) will be used in all the following computations, for instance, when evaluating the kinematical invariants (the Mandelstam s variable of the reactions involved in the GZK cut-off phenomenon).

Generalized covariance
In the HMSR model every particle has its personal modified spacetime, which is parameterized by the particle momentum. This means that every physical quantity related with a given particle is generalized and lives in a spacetime that acquires an explicit dependence on the particle energy. It is therefore necessary to introduce an original formalism to correlate different local spaces, using the generalized vierbein elements as projectors from every spacetime to a common flat Minkowski support space. Here we report a scheme of how the correlation between different local spaces is established: The Greek indices refer to the local curved geometric structures, whereas the Latin ones refer to the common Minkowski support space. Referring to the previous scheme, using the vierbein elements as projectors, the generalized Lorentz transformations can be obtained: where the elements Λ b a belong to the Lorentz group and are defined in the flat Minkowski spacetime and the used vierbein is defined in eq. (38) and includes all terms up to the perturbative order under consideration. The Lorentz covariance is promoted to a diffeomorphism invariance and the introduced class of We point out that the previous construction about the generalization of covariance is valid if the spacetime is metrizable and admits a vierbein, hence the result is valid even in different contexts. Therefore, the present prescription can be used to pose analogous definitions and set the stage of the model in the context of the more restrictive Finsler geometry.

Affine and spinorial connections
The geometrical structure is characterized by the affine and the spinorial connections. In Finsler geometry the affine connection, that is the Christoffel symbol, can be defined using eq. (14). In the case of zero spacetime curvature, the nonlinear connection eq. (12) is N µν = 0 and therefore the derivative eq. (14) reduces to the ordinary one: Since in the case of absence of curvature, the metric depends only on the momentum (or equivalently on the velocity) and not on the coordinates, the Christoffel symbol eq. (14) becomes: The explicit form of the covariant derivative is a consequence of this result and it is equal to the ordinary derivative in flat spacetime: Now it is possible to define the spinorial connection using the covariant derivative: It is simple to demonstrate that all the connection coefficients vanish since the vierbein depends on the momentum but not on the coordinates. Finally we obtain the total geometric covariant derivative, which will be useful in defining the minimal extension of the Standard Model (SM) of particle physics: The resulting spacetime geometry is therefore a flat Finsler pseudo-structure [34,33,35,36,37,38], whereas the associated momentum space is asymptotically flat.

Modified Poincaré brackets
As a final result about the geometrical structure introduced by HMSR, we report the modified Poincaré brackets, computed for the local structure using the vierbein projectors: where { x µ } are the local coordinates and {x i , e j ν (p)} = 0 since the vierbein e j ν (p) is a function of the momentum p µ . Hence, time and space coordinates do not commute anymore [21]. Using the explicit form of the vierbein eq. (38) and admitting for the perturbation functions the approximations 3 : it is possible to obtain the following relations valid at the first perturbative order: with the antisymmetric matrix θ µν that satisfies the relations: in the high energy limit, proving that, if the parameters α and β are different, the spacetime coordinates do not commute anymore [39].
The relation existing between noncommutative field theory and LIV was analyzed in [40]; here we introduced a framework to investigate this point in the context of modified covariance.

Minimal extension of the Standard Model in an isotropic scenario
In the HMSR framework the SM of particle physics can be amended in order to include the QG-caused perturbations preserving covariance of the theory, even if in an amended formulation. Following, for instance, a strategy analogous to that used for the isotropic sector of the Standard Model Extension (SME) [9], the theory formulation requires the definition of the modified Dirac matrices with the related Clifford algebra and spinors.

Modified Clifford algebra and spinors
The Dirac matrices acquire an explicit dependence on the particle momentum and again using the vierbein projectors eq. (38) they can be written as: The γ 5 matrix is constant and this means that the chiral projectors are not affected.
The Dirac matrices modified via eq. (50) satisfy the defining relation of the Clifford algebra: The definition of spinor fields is now amended preserving the usual plane-wave formulation: The normalization of the spinors u r (p) and v r (p) is modified since these are defined using the newly introduced metric eq. (22) and the related internal product.
From the previous definitions the modified Dirac equation can be derived: An important consistency check of the new formalism can be obtained verifying that eq. (53) implies the MDR eq. (1).
As a final result, a minimal extension of the SM can be obtained from the formalism here introduced. Indeed using the vierbein elements to project to a common Minkowski support spacetime the physical quantities related to different interacting particles and using the explicit form of the total covariant derivative eq. (44), the minimal extension of the SM can be formulated for a flat spacetime. Here we illustrate the amended formulation of quantum electrodynamics (QED), whose Lagrangian can be written in the form: where the term | det [g]| is borrowed from the formulation of QFT in curved spacetime. The vierbein element e is related to the gauge field A a (x). The gauge field is supposed as Lorentz-covariant in the usual meaning, that is, the MAV of photons is the usual speed of light c. The gauge field therefore is set on a Minkowski spacetime (T M, η ab ) and the vierbein is given by: e µ a = δ µ a . The QG corrections can be introduced in the generic gauge boson sector modifying the definition of the vierbein related to the gauge field, as done for the massive fermion fields.
The interaction is governed by the conserved current.
In the low-energy scenario the covariance perturbations in the conserved current are negligible, whereas in the high-energy limit the formulation admits a constant form, since one can consider the incoming and outgoing momenta with the same constant high-energy limit.

Gauge symmetry
The SM minimal extension obtained in the context of HMSR preserves the classic internal gauge symmetry SU (3) ⊗ SU (2) ⊗ U (1). This result can be stated formulating an amended version of the Coleman-Mandula theorem [16].
The modified gauge symmetry group acquires the explicit form: where P({p}) is the kinematical symmetry group and is given by the direct product of the momentum-dependent Poincaré groups associated with the different particle species and G int is the internal gauge symmetry group, in this case the usual SM gauge group:

Modified kinematics
HMSR theory perturbs the kinematics geometrizing the interaction of free propagating particles with the supposed quantized structure of the spacetime.
In this work we investigate the phenomenological effects introduced by QG in the UHECRs propagation, where the main effects are caused by the modification of the kinematics. The detectable effects are caused by the interaction of different particle species, which modify in a proper way the related spacetime.
Introducing a generalized internal product in the momentum space for more than two different interacting particle species, it is possible to obtain a generalized formulation for the Mandelstam variables s, t and u, which are the dynamical invariant quantities linked to a reaction.
Considering p and q as the momenta of two interacting particles of different species, the internal product of their sum can be defined as: where the vierbeins e µ a (p) andẽ µ a (q)) are associated, respectively, with the two different particle species. The formulation of the modified internal product can be simplified in the form: using the generalized metric: The inner product defined in eq. (61) is invariant with respect to the Modified Lorentz Transformations (MLT) introduced in HMSR: The inner product remains invariant under the action of such generalized Lorentz transformations, indeed one can obtain: Λ · G · Λ t is the metric evaluated for the two particle momenta Λp,Λq.
The new formalism here introduced guarantees that HMSR theory can deal with the interaction of different particle species in the context of a QG-modified kinematics without the necessity of the introduction of a preferred reference frame.

Ultra-high-energy cosmic-ray propagation
Before investigating the QG-induced phenomenology in UHECRs physics, it is useful to illustrate the standard physics predictions for these highly energetic particles. During their propagation, UHECRs can interact with the CMB and depending on their nature and energy are attenuated. For instance, a propagating CR bare nucleus N with atomic number A can undergo a photo-dissociation process interacting with the CMB: where n represents a generic nucleon. The proton propagation is influenced by different interaction processes with the CMB, hence it can undergo, for instance, a pair production process: This process is the main interaction mechanism with the CMB for CR protons with an energy lower than a threshold E ∼ 5 × 10 19 eV. The dominant process for UHECR protons with an energy exceeding this threshold is the ∆ particle resonance photopion production process: photopion production is the main attenuation interaction of UHECRs and is the core mechanism for the so-called GZK cut-off phenomenon for protons [1,2]. This effect poses an upper limit on the energy of protons detected at ground and coming from distant sources. Since through this effect a particle dissipates energy but is not annihilated, a proton with enough energy can undergo the same interaction process again and can undergo what is a stochastic dissipation process. This way, it becomes possible to evaluate the attenuation length, defined as the average distance that the proton has to travel in order to reduce its energy by a factor of 1/e. The inverse of the attenuation length is given by [41]: where σ pγ (s) is the proton-photon interaction cross section as a function of the squared center of mass energy (the Mandelstam s variable), n( ) represents the background photon density per unit volume and photon energy , µ = cos θ is the impact parameter and th is the interaction threshold energy. K(s) represents the reaction inelasticity, that is, the energy fraction available for secondaryparticle production during the reaction. Complementary to the inelasticity is the elasticity function, defined as the energy fraction preserved by the primary particle after the interaction, η = E out /E in with the incoming particle energy E in and the residual energy E out . Elasticity and inelasticity are connected by the simple relation: The the Mandelstam s can be computed introducing the photon four-momentum ( , p γ ) defined in the rest frame of the nucleus. In the high-energy limit approximation for the proton velocity v p 1, with ds = −2E p d cos θ, the following relations hold: hence, the inverse of the attenuation length eq. (68) can be written as: where the primed quantities are defined in the proton rest frame, whereas the other quantities are defined in the laboratory rest frame. Since the n( ) distribution is a Planckian function of the temperature T , a further simplification of the previous relation of eq. (72) is possible, obtaining the explicit form for the inverse of the attenuation length [4]: The inelasticity computed for the standard physics case is given by the relation [41]:

QG introduced phenomenology in UHECR propagation
In this work, following the theoretical framework introduced by HMSR, the kinematics of free particles is modified in order to include the QG effects. The kinematical perturbations modify the allowed phase space for the reaction and may therefore influence the processes involved in UHECR propagation. The introduction of QG phenomenology can indeed affect the photopion production mechanism, the core reaction underlying the GZK phenomenon [4,5,6,7]. The QG-caused reduction of the allowed phase space can modify the inelasticity K eq. (74). The consequent reduction of the inelasticity means that the incoming UHECR proton dissipates less energy during the GZK process and therefore the resulting opacity horizon is enlarged with respect to the standard physics prediction. In this work we conduct our analysis in the HMSR framework [16], but in the context of a version of the Finsler geometry obtained by starting from a MDR written in a more general form of eq. (1) than the one given by eq. (3) and used in previous works [6,7].

Constraints from ∆ resonance creation
The photopion production requires a ∆ particle creation (eq. (67)) and can occur passing through a real ∆ particle, in the case of the dominant process, or alternatively through a virtual one. The introduction of QG phenomenology can presumably modify only the GZK cut-off phenomenon, enlarging the foreseen opacity sphere without a complete suppression. The Universe opacity to the propagation of CR seems to be confirmed by experimental evidence, the uncertainty concerns only the dimension of the horizon. The production of the real ∆ particle must therefore be preserved in order to foresee only small QG-induced deviations from the classical predictions for the GZK process. The four-momenta used in the following computations are considered as covariant vectors in order to simplify the use of the MLT for changing the reference frame.
In the proton-CMB interaction process the production of a ∆ particle is allowed if the proton-CMB interaction free energy, that is the Mandelstam √ s variable, exceeds the rest energy of the particle resonance. As a consequence of this threshold energy consideration, a first constraint on the magnitude of the QG perturbation parameters can be posed. Using the generalized internal product and the formalism developed for the modified kinematics in the eqs. (60) and (61) one can write the Mandelstam s variable as a function of the proton and the photon four-momenta respectively p (p) = (E (p) , p (p) ) and p (γ) = (E (γ) , p (γ) ): The vierbein eq. (38) has been used to project the momenta of different particles to the common support Minkowski spacetime and the symmetry of the mixed product of p p and p γ has been used. Photons are assumed to be Lorentzinvariant, hence the related vierbein is defined as usual: e a µ (p γ ) = δ a µ . Using the explicit form of the squared Finsler norm F 2 (p) given in eq. (34), the following inequality can be derived reordering the terms of the previous relation in an opportune manner: where h p (p p ) and k p (p p ) are the proton QG perturbation functions taken from the MDR eq. (38). The following inequality can be derived from the previous eq. (76) using the MDR eq. (1): where m ∆ 1232 MeV, m p 938 MeV. The covariant formulation of the model is fundamental since it allows the change of the reference frame. Indeed, in order to simplify the computation, it is possible to choose the more suitable reference frame, where the involved energy scales ratio allows the suppression of some perturbative terms. One can therefore consider the laboratory frame, where the proton energy is much higher than the CMB one. Considering that the UHECR proton energy has an upper limit E p 10 21 eV and taking into account the tiny average value of the CMB energy E γ 1.16 × 10 −3 eV, one can neglect the O(1) contributions in the second term obtaining: Defining the function f p = h p −k p the result obtained is comparable with the one presented in [6,7]. The derived inequality must be satisfied in order to produce a ∆ resonance, otherwise the QG perturbations totally suppress the GZK effect.
Imposing the validity of the previous relation the following constraint can be derived: Substituting the approximated average value of the CMB energy in the previous inequality one can obtain the approximated constraint: This result is comparable with the upper limit 4.5 · 10 −23 obtained numerically in [4] for the perturbation magnitude.
As a final remark from the relations eqs. (77) and (79) it follows that in order to generate visible effects on the GZK phenomenon, the QG perturbation function f p must satisfy the relation: Indeed from eq. (77) the high-energy limit excludes the possibility of any suppression of the GZK cut-off in the case of a negative function f p . We underline that the choice of positive values of the perturbation f p corresponds to the introduction of an effective MAV in the MDR of eq. (1) for every massive particle, which is lower than the standard speed of light.

Modified inelasticity
The introduction of QG perturbations can modify the phase space allowed for the photopion production process, determining a modification of the inelasticity function K(s) eq. (74). In this work we evaluate the new inelasticity in the context of the generalized Finsler geometry here introduced following the approach of [6,7], originally inspired by [4,5]. In the following, all the computations are conducted again considering the four-momenta as covariant vectors, in order to simplify the computations transforming the reference frame with the MLT. The introduction of the center-of-mass (CM) reference frame is useful and it is defined via the relation: where the momenta are defined on the common Minkowski support spacetime (T M, η ab ) and * labels the quantities related to the CM frame.
The next element necessary in evaluating the change from the CM to a generic reference frame is the γ CM Lorentz factor. Starting from the free energy available for the photopion production √ s = E * p + E * π one can obtain the relation: Using the CM definition eq. (82) the relation | p * p | = | p * π | follows and the fourmomentum of the photopion can be written in the CM reference frame as . From the latter, the free energy necessary to produce a photopion in the CM frame can be computed using the squared Finsler norm eq. (34): h p and k p are the perturbation functions of the proton MDR of eq. (1) and h π and k π are the analogous functions for the pion MDR, p and π label, respectively, the elements related to the proton and the pion.
We underline that the computation made here for the modifications induced on the inelasticity are still valid in the context of the more restrictive Finsler geometry. Indeed, the computation is conducted using the MDR eq. (1), that where K π (θ) is the pion inelasticity depending on the impact angle.
In the high-energy limit, where for ultra-relativistic particles | p| ∼ E and the velocity factor β ∼ 1 and γ CM is approximated by eq. (83), the following equation can be derived: Posing f p = (h p − k p ) in the previous relation, it is possible to obtain the following one: where the latter result was obtained in [6,7]. The previous equation can be numerically solved in order to evaluate the inelasticity as a function of the collision angle θ. The numerical computation is conducted in the high-energy limit, hence by using eq. (86) and considering √ s = E where E is the initial energy for the process.
The free energy s can be written as: introducing the energy of the photon defined in the proton rest frame p p = 0.
The inelasticity is finally averaged over the interval θ ∈ [0, π]: In the following, we plot the inelasticity for different choices of the parameter . The parameter f p is constrained to be positive in order to guarantee detectable QG effects. This hypothesis corresponds to the introduction of a MAV inferior to the classical speed of light c for every massive particle species.
The parameter is even constrained from above by the limit obtained in the previous section eq. (80), hence the plausible parameter must be included in the interval 0 < f p 10 −22 . The inelasticity is plotted as a function of the proton energy E p , defined in the laboratory reference frame, and the photon energy , defined in the proton rest frame. Under these hypotheses, the QG perturbations cause a dramatic drop of the inelasticity value, which is more visible for increasing values of the parameter f p . This implies a reduction of the allowed phase space for the photopion production process, and an enlargement of the foreseen opacity horizon for an increasing magnitude of the QG perturbation.

Simulated attenuation lenght
The impact of the QG perturbations on the UHECR free propagation can be evaluated computing the value of the attenuation length as a function of the modified inelasticity K in eq. (73). In this work we have obtained the attenuation length using an ad hoc modified version of the simulation software SimProp [17]. In this software version we substituted the inelasticity predicted by the standard physics with the modified one, being a function not only of the proton and the CMB energy but also of the QG parameter. In Figs

Conclusions
In this work we investigated the possibility of exploiting the UHECR physics to detect supposed signatures of the quantum structure of spacetime. Nowadays the Lorentz covariance stands at the base of our physics knowledge, but small departures can be induced by QG effects. We studied the induced QG phenomenology on UHECRs simulating their propagation in a framework of Lorentz covariance modification. The model is a generalization of the HMSR theory, that foresees a minimal extension of the particle SM, preserves the usual gauge symmetry SU (3) ⊗ SU (2) ⊗ U (1), and does not introduce exotic particles or reactions. The spacetime geometry is a generalized Finsler spacetime and all the GZK cutoff modification effects have been computed in this context. We underline that the methodology introduced in this work can be used to set the stage of the model in the context of the more restrictive Finsler geometry. Moreover, we have demonstrated that, for instance, the inelasticity modification that underlies the opacity sphere enlargement can be still computed in the Finsler geometric spacetime.
The conducted simulations foresee an enlargement of the attenuation length in accordance with other works [4,5,6,7,42,43], but are obtained in a covariancepreserving scenario as proposed in [44]. The kinematical solution proposed in our work is based on the phenomenological effects caused by the interaction of different particle species which are modified in different ways by their interactions with the QG structure of spacetime. Our proposal is analogous to that of [3], but is based on a covariant formulation, such as in DSR theories [45,46].
Moreover, in our proposal a way is outlined to generalize HMSR theory in curved spacetime and a formalism is introduced that defines the threshold energy of reactions in a modified covariant framework.
We emphasize the importance for this work of preserving covariance, even if in an amended formulation. The preserved covariance allows to define relativistic kinematical invariants. As a direct consequence, it is possible to choose the most suitable reference frame to simplify the computations, suppressing some perturbation terms thanks to the involved energy scale ratios. Furthermore, in the astroparticle sector covariance can be a great experimental advantage, since all the obtained data can be collected without the necessity to introduce any sidereal discrimination related to the orientation of the detector with respect to a fixed privileged reference frame.
In this work we set the stage for an analysis strategy that can be improved including a heavier UHECR component, which is a more realistic CR composition, as indicated by the Pierre Auger collaboration [47]. The QG perturbations are expected to be larger for a CR component heavier than protons. Indeed, the kinematical perturbations can affect the propagation of bare nuclei in a more significant way. The perturbative order of the matrix A follows from a direct computation. Indeed, the entries of the matrix A involve first and second-order derivatives of the perturbation functions h and k multiplied, respectively, by the momentum to the first or the second power.