Quantum gravity phenomenology via Lorentz violations
Abstract
The search for a quantum theory of gravity has been one of the main aims of theoretical physics for many years by now. However the efforts in this direction have been often hampered by the lack of experimental/observational tests able to select among, or at least constrain, the numerous quantum gravity models proposed so far. This situation has changed in the last decade thanks to the realization that some QG inspired violations of Lorentz symmetry could be constrained using current experiments and observations. This study it is not only allowing us to test at higher and higher energies a fundamental symmetry of spacetime but it is also providing us with hints and perspectives about the fundamental nature of gravity.
Quantum gravity phenomenology via Lorentz violations
Stefano Liberati^{†}^{†}thanks: Speaker. ^{†}^{†}thanks: liberati@sissa.it
SISSA/ISAS and INFN, Trieste, Italy
Email: liberati@sissa.it
\abstract@cs
1 Introduction
The quest for a quantum theory of gravity has been one of the main enterprise of modern physics for more than half a century. However this research has been systematically frustrated by the lack of experimental/observational test of the models so far proposed.
Albeit many intriguing and ingenious ideas have been explored, it seems safe to say that without both observing phenomena that depend on Quantum Gravity (QG), and extracting reliable predictions from candidate theories that can be compared with observations, the goal of a theory capable of incorporating quantum mechanics and general relativity will remain unattainable.
Of course one can claim that a primary tests to be passed by any candidate QG theory is to admit a suitable semiclassical limit, i.e. to recover GR at sufficiently low energies. However this is is more a consistency requirement than a true prediction to be tested for a given QG model. In this sense the only example of a prediction of a quantum gravity model that has received some support from observation is the spectrum of primordial cosmological perturbations where the quantized longitudinal linearized gravitational mode, albeit slave to the inflaton and not a dynamically independent degree of freedom, plays an essential role [1]. However also in this case one might object that somehow any theory of quantum gravity that admits a GR limit should admit as well a regime where gravitons are a valid concept.
Looking for empirical evidence the last decades have witnessed an increasing number of ideas about observable phenomena where QG could play a key role. A partial list includes: deviations from Newton’s law at very short distances [2, 3], Planckscale fuzziness of spacetime [4], possible production in TeVscale QG scenarios of miniblack holes at colliders [5] or in cosmic rays [6], QG induced violations of discrete symmetries of the Standard Model [7, 8] as well as spacetime symmetries [9]. This broad field of research goes under the general name of quantum gravity phenomenology.
We shall focus here on this last item, more specifically on the possibility that the Lorentz symmetry could be violated or deformed (we shall explain later what we exactly mean by this) by Plancksuppressed corrections. Given the necessary conciseness of these proceedings we shall not cover all the related issues, we direct the reader interested in deepening this subject to some of the extensive reviews available in the literature (see e.g. [10]). Discussion and results presented in this proceeding can be also be found (in an extended form) in the following works [11, 12, 13, 14, 15, 16, 17].
2 Lorentz symmetry at the Planck scale
Lorentz symmetry has been confirmed to ever greater precision, and it powerfully constrains theories in a way that has proved instrumental in discovering new laws of physics. It seems then natural to assume under these circumstances that Lorentz invariance is a symmetry of nature up to arbitrary boosts. Nevertheless, from a purely logical point of view, it is clear that such a conjecture is empirically not strongly motivated given that an infinite volume of the Lorentz group is (and will always be) experimentally untested being the Lorentz group noncompact, unlike the rotation group. Why should we assume that exact Lorentz invariance holds at all scales when this hypothesis cannot even in principle be tested? In this sense high energy tests of Lorentz symmetry are “per se” valuable being significant improvements of the precision with which an apparently fundamental symmetry of nature is tested.
While the above reasoning is logically sound, it is by itself not very encouraging as it is not providing any argument by which we should expect any departure from Lorentz symmetry due to QG. However, there are also several arguments that lead to suspect that there could be a failure of Lorentz symmetry in proximity of the Planck scale such as the ubiquity of UV divergences in quantum field theories or profound implications related to the very construction of a quantum theory of spacetime (see e.g. [18, 19, 20]).
Aside from general issues of principle, specific hints of Lorentz violation have come from tentative calculations in various approaches to QG hinting that this might be the case. Just to cite a few, we can recall string theory tensor VEVs [21], spacetime foam [22], semiclassical spinnetwork calculations in Loop QG [23], noncommutative geometry [24], some braneworld backgrounds [25], and condensed matter analogues of “emergent gravity” [12].
Of course, lacking a definitive theory of QG one cannot claim that there is a convincing prediction that some departure from Lorentz symmetry must be a feature of quantum gravity. However, taken together they do motivate the effort to characterize possible observable consequences of LV and to strengthen observational bounds. Moreover, although very different, these models have in common the fact that Lorentz violations express themselves through modified dispersion relations for elementary particles. Generically one can cast them in the form:
(2.0) 
where, , are the energy and momentum of the particle, is some particle physics mass scale (possibly the mass of the particle ) and denotes the mass scale at which the quantum gravity corrections become appreciable. Normally, one assumes that is of order the Planck mass: GeV.
We can expand the above dispersion relation in powers of the particle momentum^{1}^{1}1We assume, for simplicity, that rotational invariance is preserved and only boost invariance is affected by Planck scale corrections.
(2.0) 
The dependence on the mass terms and is now hidden in the dimensionfull coefficients . In general the idea that the above dispersion relations are a sort of mesoscopic effect of some Planckian physics leads to the natural expectation that the extra term in Eq. (2) contains suitable powers of the mass terms that make sufficiently Planck suppressed for being compatible with low energy observations. For the moment we simply follow the observational lead and insert at least one inverse power of in each term
(2.0) 
Later on, in characterizing the strength of a constraint we refer to the without the tilde, so we are comparing to what might be expected from Plancksuppressed LV. We shall also allow the LV parameters to dependent on the particle type .
It is remarkable that several significant constraints can be put on the intensity of the Lorentz violating term using current experiments and observations [10]. However before discussing some of such constraints we want to provide here an explicit toy model of a system which does show modified dispersion relations of the sort we just conjecture and which for this reason could in fact provide a test field for our ideas. Such a framework is provided by the so called analogue models of gravity [12].
3 Analogue models of gravity
Reduced to the bones most of the analogue models of gravity can be described as condensed matter systems where the dynamics of atoms or molecule provides an emergent spacetime geometry for the propagation of collective excitations of the background. In this sense analogue models of gravity are indeed analogue models of emergent gravity frameworks, i.e. of scenarios in which the metric and the affine connections are collective variables emerging from the dynamics of more fundamental objects (note that in this sense the metric or the connections would not be the right object to quantize the same way you would not expect to recover atomic interactions from the quantization of hydrodynamical variables).
The simplest example of an analog model is provided by a barotropic — , with and respectively equal to the density and pressure of the fluid — and inviscid fluid whose flow is irrotational (though possibly time dependent). For such a system the equation of motion of the velocity potential of an acoustic disturbance (defined as , where is the velocity of the flow) is identical to the d’Alembertian equation of motion for a minimally coupled massless scalar field propagating in a –dimensional Lorentzian geometry
(3.0) 
This signifies that under the above conditions, the propagation of sound is governed by an acoustic metric — . This acoustic metric describes a –dimensional Lorentzian (pseudo–Riemannian) geometry. The metric depends algebraically on the density, velocity of flow, and local speed of sound in the fluid. Specifically
(3.0) 
where is the speed of sound.
Correspondingly the dispersion relation for the quanta of sound (which can be meaningfully defined in some coherent systems like e.g. superfluids), the so caled “phonons”, will take the standard relativistic form
(3.0) 
where the role of the speed of light is now played by the speed of sound.
The interesting point for our discussion here is that in a realistic condensed matter system the underlying microscopic structure of the background (e.g. the fact that the fluid is made up of molecule) will generically show up in a breakdown of the acoustic Lorentz invariance of the phonon equation of motions which will lead to modified dispersion relations of the form (2) where now will correspond to the energy scale that marks the transition from an hydrodynamic description to a molecular/atomic one (e.g. intermolecular distance or coherence length). In order to explicitly see this we have however to fix a particular analogue model and a well known example in this sense is that of a BoseEinstein condensate [12, 26, 27, 28, 29].
Bose–Einstein condensates:
Let us start by very briefly reviewing the derivation of the acoustic metric for a BEC system, and show that the equations for the phonons of the condensate closely mimic the dynamics of a scalar field in a curved spacetime. In the dilute gas approximation, one can describe a Bose gas through a quantum field satisfying
(3.0) 
Here g parameterizes the strength of the interactions between the different bosons in the gas. It can be reexpressed in terms of the scattering length as
(3.0) 
The quantum field can be separated into a macroscopic (classical) condensate and a fluctuation: , with . After some suitable approximations (basically neglecting the backreaction of the excitations on the background, see e.g. [12, 30]), and adopting the so called Madelung representation for the wave function of the condensate
(3.0) 
the equation describing the evolution of the background (classical) field (Gross–Pitaevskii equation) can be rewritten as a continuity equation plus an Euler equation:
(3.0)  
(3.0) 
where we have defined the irrotational “velocity field” by .
Let us now note that these equations are completely equivalent to those of an irrotational and inviscid fluid apart from the existence of the socalled quantum potential
(3.0) 
From what we previously saw it should then be clear that if one neglects the quantum potential then the propagation of acoustic perturbations on the BEC background will be again characterized by an acoustic geometry of the form of (3).
More precisely if one defines as a quantum excitation of the phase of the background (which here plays the role of the velocity potential in the perfect fluid) then one finds that the equation for the field becomes that of a (massless minimally coupled) quantum scalar field over a curved background
(3.0) 
with an effective metric of the form
(3.0) 
where the magnitude represents the speed of the phonons in the medium:
(3.0) 
Let us now consider the case in which the above “hydrodynamical” approximation for BECs does not hold. In order to explore a regime where the contribution of the quantum potential cannot be neglected we can use the so called eikonal approximation, a highmomentum approximation where the phase fluctuation is itself treated as a slowlyvarying amplitude times a rapidly varying phase (see e.g. [12, 30] for further details). Specifically, we shall write . As a consequence of this assumption, gradients of the amplitude, and gradients of the background fields, are systematically ignored relative to gradients of . Then adopting the notation
(3.0) 
one can show that the dispersion relation for the excitations of the BEC (quasiparticles) in a flat background (i.e. when the fluid is at rest, ) takes the form:
(3.0) 
which is just a special case of the general form of dispersion relations (2) we consider in this work and corresponds to the so called Bogoliubov dispersion relation which was found in 1947 [31] via a diagonalization procedure for the Hamiltonian describing the system of bosons . Note also that the “quantum gravity scale” or more correctly the scale of the violation of the acoustic Lorentz invariance is the wavelength with respect to the “acoustic Compton wavelength” . In particular, for large wavelengths one gets a standard phonon dispersion relation . For wavelengths the quasiparticle energy tends to the kinetic energy of an individual gas particle and in fact (as it should be given that the BEC is in the end a system of many interacting nonrelativistic atoms as described by the nonlinear Schr̈odinger equation (3) ).
The above discussion clearly shows the potentiality of the analogues systems as an inspirational tool for understanding the interplay of microphysics and macrophysics in quantum gravity phenomenology of Lorentz violations, where the modified dispersion relations of the form (2) can be considered intrinsically a mesoscopic physics effect relevant at intermediate energies between our low energy world and the Planck scale. Of course one has to keep in mind that while the relativistic (low energy w.r.t. ) limit of the dispersion relation above is quite generic in many models (the microscopic parameters will in general only determine the value of the speed of sound ), the “high energy” deviations (i.e. the phenomenology at energies comparable to M) will not be generic but instead strongly dependent on the microphysics of the system. Nonetheless we shall see later on that sometimes useful lessons about viable mechanisms in nature for generating the dispersion relations under investigations can be learn and possibly exported to the more realistic scenarios.
4 Modified dispersion relations
Let us now move back to the problem of actually placing constraints on the departures from Lorentz invariance manifested in dispersion relations like (2). It was perhaps one of the greatest achievements of the last decade the realization that such constraints can be cast (see e.g. [10] for an extensive review). In fact it was realized that in some special situations even tiny corrections, like the one considered here, can be magnified to observable effects. A partial list of these “windows on quantum gravity” includes:

sidereal variation of Lorentz violation (LV) couplings as the lab moves with respect to a preferred frame or directions

cumulative effects: long baseline dispersion and vacuum birefringence (e.g. of signals from gamma ray bursts, active galactic nuclei, pulsars)

anomalous (normally forbidden) threshold reactions allowed by LV terms (e.g. photon decay, vacuum Čerenkov effect)

shifting of existing threshold reactions (e.g. photon annihilation from blazars, GZK reaction)

LV induced decays not characterized by a threshold (e.g. decay of a particle from one helicity to the other or photon splitting)

maximum velocity (e.g. synchrotron peak from supernova remnants)

dynamical effects of LV background fields (e.g. gravitational coupling and additional wave modes)
It has however to be stressed that not all of the above cited tests are in the same way robust against the underlying physical framework that one is choosing in order to justify the use of the modified dispersion relations (2). In fact while the above cited cumulative effects use exclusively the form of the modified dispersion relations basically all the other effects are dependent on the underlying dynamics of interacting particles and on the fact that the standard energymomentum conservation holds or not. Hence in order to cast most of the constrains on dispersion relations of the form (2) one needs to adopt a specific theoretical framework justifying the use of such deformed dispersion relations.
5 Theoretical frameworks for Lorentz violation
Dsr:
Quantum gravity effects seems to introduce a new dimensional fundamental scale given by the Planck length (or the corresponding area/volume). Deformed or Doubly Special Relativity (DSR) [32, 33, 34] can be understood as a tentative modification of Special Relativity (SR) in order to incorporate some QG scale (generally the Planck length) as a new invariant scale other than that provided by the speed of light but preserving at the same time the relativity principle.
This very ambitious program has so far found a partial realization only in momentum space where the DSR idea has been implemented by introducing a deformation of the Poincaré algebra in the boost sector [33, 35]. Specifically the Lorentz commutators among rotations and boost are left unchanged but the action of boosts on momenta is changed in a nontrivial way (see e.g. [33]) by corrections which are suppressed by .
Concretely one considers the standard Lorentz algebra of the generators of rotations, , and boosts, :
(5.0) 
(Latin indices run from 1 to 3) and supplements this with the following (deformed) commutators between the Lorentz generators and those of translations in spacetime (the momentum operators and ):
(5.0) 
(5.0) 
(5.0) 
Finally, assume
(5.0) 
The commutation relations (5)–(5) are given in terms of three unspecified, dimensionless structure functions , , and , and are sufficiently general to include all known DSR proposals — the DSR1 [32], DSR2 [33], and DSR3 [34]. Furthermore, in all the DSRs considered to date, the dimensionless arguments of these functions are specialized to
(5.0) 
so that rotational symmetry is completely unaffected. Furthermore, in order to recover ordinary special relativity in the limit, one has to demand that, in that limit, and tend to , and that tend to some finite value.
It was soon recognized [36] that such deformed boost algebra amounts to the assertion that physical energy and momentum of DSR can be always expressed as nonlinear functions of a fictitious pseudomomentum , whose components transform linearly under the action of the Lorentz group. More precisely one can assume the existence of an invertible map between two momentum spaces: the classical space , with coordinates where the Lorentz group acts linearly and the physical space , with coordinates , where the Lorentz group acts as the image of its action on . Also, must be such that for all elements on with and/or .
The main open issues in this momentum formulation of DSR are the so called multiplicity and saturation problems. The first is related to the fact that in principle there are many possible deformations (an infinite number, depending on the choice of an energy invariant scale, threemomentum scale or both [35]). This seems to suggest that the set of linear transformations (that is SR) is the only one that have a physical sense. Moreover the composition law for energy and momenta of DSR, being derived by imposing a standard composition law for the pseudo fourmomenta , is characterized by a saturation at the Planck scale apparently in open contrast with the everyday life observation of classical objects with transplanckian energies and momenta.
On the other hand, since DSR is not a formulation of QG, but gives a set of transformations with the typical QG scale, it might be plausible to consider it as a low energy limit of QG, that is, as some effective theory. (Indeed such point of view was taken in several works on the subject [37, 38, 39].) Along this line of though, it was recently proposed [13] that DSR could be interpreted as an effective theory of measurement for high energy particle properties. According to this framework, the relation between “true” energy and momentum of a particle (the classical variables of DSR) and observed quantities (the physical variables of DSR) acquires, at sufficiently high energies, Planck suppressed distortions induced by quantum gravity effects. These relations can be identified with DSRtype deformations. In [14, 15] it was further argued that this non linear nature might arise as a result of the unavoidable averaging over QG fluctuations of the metric around flat spacetime which is required in order to properly define energy and momentum in first place.
The main observational constraints on DSR type dispersion relations concerns limits on the delays between arrival times of sharp features of different energies observed in the intensities of radiation from very distant astrophysical sources. Not very much else is available at the moment given that DSR type dispersion relations do not generically show birefringence and that we cannot resort to anomalous threshold reactions as these phenomena are not allowed in DSR (the reason for this being simply that a kinematically forbidden reaction in the “classical” variables cannot be made viable just via a nonlinear redefinition of momenta). There have been attempts to consider constraints provided by shifts of normally allowed threshold reactions [40]. However for such reactions the possible constraints are strongly dependent not just on kinematical considerations, but also on reaction rates which require some working framework for their derivation, which presumably would take the form of some Effective Field Theory (EFT).
The problem in this case is the present lack of a DSR formulation in configuration space. In this sense the main stream approach in the DSR community has been oriented towards the possibility that the configuration space description of the just discussed deformed symmetry in momentum space could imply spacetime noncommutativity [41, 42]. While some promising results linking DSR to some special form of spacetime noncommutativity have been found in 2+1 dimensions [43], we are still lacking a consistent physical picture in 3+1 dimensions. Moreover attempts to develop a quantum field theory associated with different forms of noncommutativity, a much needed step in order to be able to effectively cast phenomenological constraints, led to highly nontrivial quantum field theories (possibly with problematic features such as IR/UVmixing [44]).
An alternative route it was recently proposed in [16] where it was shown that any modified dispersion relation of the form (2) can be associated with a Finsler geometry ^{2}^{2}2Finsler geometry is a generalization of Riemannian geometry: instead of defining an inner product structure over the tangent bundle, we define a norm . This norm will be a real function of a spacetime point and of a tangent vector , such that it satisfies the usual norm properties namely if , and , with . The Finsler metric is then defined as or equivalently .. While this approach can be naturally seen as a geometrical description of the violation of Lorentz symmetry (and in fact it naturally arises in analogue models of gravity [45]), it cannot be excluded at the moment that it might fit instead in a DSRlike framework if a non linear realization of the Lorentz algebra can be found leaving the Finsler norm invariant ^{3}^{3}3It was recently shown in [46] that this is not the case for the standard DSR algebra presented at the start of this section: while the equation of motions are indeed invariant the norm is so only modulo a total derivative. It is presently unclear the physical relevance of this point, given that it is not excluded that other realizations of the Lorentz algebra can leave the norm invariant, neither it is clear the observational relevance the Finsler norm once the invariance of the field equations is assured.. This seems an interesting perspective which probably deserve further attention in the next future.
In conclusion, missing a field theory implementation of the DSR idea we cannot safely pose most of the constraints listed before. We shall then move to consider the alternative route of an EFT with Lorentz violating terms.
EFT with Lorentz violation:
EFT has proven very effective and flexible in the past, it produces local energy and momentum conservation laws, and seems to require for its applicability just locality and local spacetime translation invariance above some length scale. It describes the standard model and general relativity (which are presumably not fundamental theories), a myriad of condensed matter systems at appropriate length and energy scales, and even string theory. Furthermore, it is at the moment the only framework within which we can compute reaction rates and in general fully describe the particle dynamics.
For what regards concrete realizations of this framework we can distinguish two main lines of research: one considers EFT with only renormalizable (i.e. mass dimension 3 and 4) LV operators [47, 48], the other considers instead EFT with nonrenormalizable (i.e. mass dimension 5 and higher) LV operators [49, 11].
Most of the research along the first direction has been carried out within the so called (minimal) standard model extension (SME) [47]. It consists of the standard model of particle physics plus all Lorentz violating renormalizable operators (i.e. of mass dimension ) that can be written without changing the field content or violating the gauge symmetry. For illustration, to lowest order in the Lorentz violating coefficients the dispersion relations respectively for electrons and photons in the, rotationally invariant, QED sector of the SME are [10, 47]
(5.0) 
where from here thereafter we shall identify with the electron mass keV.
The alternative approach is to study nonrenormalizable operators. All in all we nowadays consider the SM just an effective field theory and in this sense its renormalizability is seen as a consequence of neglecting some higher order operators which are suppressed in some appropriate mass scale. It is a short deviation from orthodoxy to imagine that such nonrenormalizable operators can be generated by quantum gravity effects (and hence be naturally suppressed by the Planck mass) and possibly associated to the violation of some fundamental spacetime symmetry like local Lorentz invariance.
A rigorous study of Lorentz violating EFT in the higher mass dimension sector was initiated in [49] through a classification of all LV dimension five operators that can be added to the QED Lagrangian and are quadratic in the same fields, rotation invariant, gauge invariant, not reducible to a combination of lower and/or higher dimension operators using the field equations, and contribute terms to the dispersion relation. Just three operators arise and all of the terms violate CPT symmetry as well as Lorentz invariance ^{4}^{4}4Note that while CPT violation implies LIV the opposite is not true [50].. In this case the QED sector dispersion relations come to be (in the limit of high energy )
(5.0)  
(5.0) 
were we have introduced the simplified notation and . The subscripts refer to helicity which can be shown to be a good quantum number in the presence of these LV terms [51]. Moreover, are the LV parameters of the two electron helicities, those for positrons can be shown to be [11, 51].
This is the framework in which most of the work on astrophysical constraints has been carried out (with or without taking into account explicitly the possible helicity dependence of the coefficients).
6 Constraints
We shall focus now on the QED sector with dimension 5 LV operators described by the dispersion relations Eq. (5) and (5). For these dispersion relations it is obvious that high energies are needed in order to cast constraints. How high, can be estimated considering, for example, threshold processes. In these cases sizable deviations start to appear when the LV term in the dispersion relation is of the order of the relevant mass term. For an electron this implies that one can cast a constraint of on when TeV. This energy is at the moment beyond the capacity of terrestrial experiments in particle accelerators however it is well within the realm of high energy astrophysics observations. We give here a very dry account, see e.g. [10, 11] for a detailed discussion.
Cumulative effects:
Constraints based on cumulative effects are not the strongest available but they do have a special status given that, being a purely of kinematical nature, they are the only ones independent on the choice between the above discussed theoretical frameworks.
They are cast by placing upper limits on the differences in the time of arrival at Earth of photons produced by a distant astrophysical events. In fact one can see from Eq. (5) that there is a LIV induced time delay given by
(6.0) 
which increases with the distance from the source and with the energy difference. Using single sources (generally gamma rays bursts (GRB) or active galactic nuclei (AGN)) the typical strength of the current constraints strength up to order [52].
Unfortunately these constraints suffer of large systematic errors, because it is hard to state if the photons, at different energies, are produced simultaneously in the source, even for a GRB. This problem might be avoided using the helicity dependence of (5). In fact one can consider the velocity difference of the two polarizations at a single energy [51]. This would lead to a constraint at least twice as large as the one arising from energy differences which is also independent of any intrinsic time lag between different energy photons. However would required the simultaneous detection of both photons’ polarizations and the possibility to exclude that some spurious helicity dependent mechanism (like for example the crossing of the light path of some birefringent medium) has affected the relative propagation of the two polarizations states.
Nowadays the most robust constraints on based on this mechanism are obtained via a statistical analysis on a large sample of GRBs with known redshifts. Looking at the arrival times of sharp features in the radiation intensity at different energies it is obtained a constraint of order [53].
The helicity dependence of (5) also implies that the direction of linear polarization is rotated through a frequencydependent angle, due to different phase velocities for opposite helicities. The difference in rotation angle between two wave vectors and after a propagation distance is
(6.0) 
Albeit small this vacuum birefringence can easily depolarize linearly polarized radiation composed of a spread of frequencies when this travels over astrophysical distances. Hence detection of polarized light from distant sources can cast a constraint on the strength of the LV coefficient . A reliable constraint on the dimension five term, , was deduced in [54] using UV light from distant galaxies. The much stronger constraint was derived [51, 55] from the report [56] of a high degree of polarization of MeV photons from the gammaray burst GRB021206. However, the data has been reanalyzed in two different studies and no statistically significant polarization was found [57]. Finally, using the same methodology, it was very recently derived the strongest up to date constraint by looking at polarization in the optical/UV afterglow of some gammaray bursts (GRB 020813 and GRB 021004). This limit is [58].
Anomalous threshold reactions:
These are typical phenomena allowed just in case of explicit LV however they are not very sensitive to the details of the dynamics given that the rates, once above threshold, are tremendously strong Thus we could tolerate huge modifications to the matrix element (the dynamics) and still be able to cast a strong constraint just assuming standard energy momentum conservation (i.e. no DSR).
The reactions used in this case are mainly those reactions related to the basic vertex of QED, i.e. gamma decay and vacuum Čerenkov (VC) . In both the cases the reaction happens so fast above threshold that the particle stops immediately to propagate [11] ^{5}^{5}5Above threshold the decay rate of photon goes as , while the electron rate of energy loss goes as [11].. Hence if photons or electrons are observed to propagate up to some given high energy one can infer that the latter should be still below the threshold for one of the above reactions to happen.
Taking ^{6}^{6}6Remember that from birefringence we already know that ., the photon decay threshold is
(6.0) 
In [11] was derived by requiring that that 50 TeV rays measured from the Crab Nebula had to be below threshold. Now, this constraint is somewhat better, , due to the fact that 80 TeV photons have been seen by HEGRA [59].
For the vacuum Čerenkov (taking again ), the threshold energy is given by
(6.0) 
Moreover just above threshold, this process has a time scale of the order of s, so it is extremely efficient. The strongest constraint up to now is again cast using the observation of high energy photons from the Crab nebula: the VHE ray emission of the Crab Nebula is usually interpreted as due to Inverse Compton (IC) scattering of accelerated electrons/positrons onto background photons. These leptons simply would not be able to produce the observed inverse Compton radiation if they would have been above the VC threshold (the VC rate above threshold is much higher than the IC scattering rate in the Crab).
In [11] it was used the fact that 50 TeV photons are observed from the Crab Nebula to infer that 50 TeV electrons or positrons in at least one helicity state must be propagating in the Crab nebula. Thus, the bound was deduced for one of the 4 fermion parameters. With the observation of 80 TeV photons by HEGRA [59] this is strengthened to .
Shifting of existing threshold reactions:
Two characteristic reactions of high energy astrophysics have been studied: photon pair production from TeV gamma rays hitting the infrared (IR) and cosmic microwave background (CMB) photons and the so called GZK reaction which describes the production of pions via the collision of ultra high energy protons^{7}^{7}7Note that we do not have at the moment a fundamental derivation of the dispersion relation for hadrons (i.e. for composite particles) and a form like Eq. (5) is just assumed. ( eV) on CMB, (see e.g. [60, 61, 62, 63]). In the presence of Lorentz violating dispersion relations the threshold for these processes is in general shifted. Moreover, it has been noticed [60, 63, 64], that in some cases there is an upper threshold beyond which the process does not occur. It is very difficult to cast robust constraints on reactions of this kind, in fact these are strongly dependent on the shape of the primary spectrum of the incident particles, on the backgrounds (which in the case of IR photons is not well known) and on the reaction rates. Moreover, for GZK reaction, we do have even an uncertainty about the very existence of the reaction hence it is not clear how present observations should be used. For all these reasons we do not discuss here the constraints and refer the reader to the reviews [10, 11].
LV induced decays not characterized by a threshold:
Characteristic examples in this case are decays of a particle from one helicity to the other or photon splitting. Although they are not characterized by a proper threshold energy they are strongly dependent on the decay rate which of course is totally negligible at low energies. Helicity decay can happen if the positive and negative helicity LV parameters for electrons are unequal. One can easily see [11] that the rate is nonnegligible for a given approximately at the same energy at which the threshold of the vacuum Čerenkov would be located. In this sense one might speak of an effective threshold. Unfortunately, we lack at the moment sufficiently precise observations for casting constraints using this reaction although the situation might change soon (see discussion in [11]). For what regards photon splitting we are still missing a complete calculation of the rate which would include the the helicity dependence of the photon dispersion, and the Lorentz violation in the electronpositron sector. However preliminary calculations seems to be promising casting a bound on the photon coefficient of order [65].
Maximum velocity:
Synchrotron radiation emission by electrons/positrons cycling in a magnetic field is strongly affected by LIV. In the Lorentz invariant case, as well as in the presence of Lorentz violation [66, 67], most part of the radiation due to an electron of energy is emitted around a critical frequency
(6.0) 
where , and is the electron group velocity.
The dispersion relation (5) implies that electrons (or positrons) with a negative value of will have a maximal group velocity smaller than the low energy speed of light. Consequently there will be for them a maximal synchrotron frequency that can be produced, regardless the energy of the radiating lepton. Thus for at least one electron or positron helicity must be greater than the maximum observed synchrotron emission frequency . This yields a constraint which is strongest for a system that has the smallest ratio [66]. Presently such a system is the Crab nebula, which emits synchrotron radiation up to 100 MeV and has a magnetic field no larger than 0.6 mG in the emitting region. Thus one can infer that at least one lepton population must have a LIV coefficient greater than (note this is not per se a constraint, see discussion below).
Combined constraints:
Putting all together we get the following picture where all the relevant information is given by the Crab nebula.
In summary we see that the coefficients are nowadays restricted to the region by birefringence and by photon decay. It is clear that while the constraint on the photon coefficient is remarkably strong not the same can be said about the LIV coefficients of the leptons. Of course we have a comparably strong estimate on the leptons coefficients, namely the synchrotron one, but this is not a double sided constraint and it only implies that the LIV coefficient of the population responsible for the Crab synchrotron emission cannot be more negative that . The same way the vacuum ČerenkovIC bound only tell us that at least one lepton population must satisfy it. These statements, although not void of physical significance, cannot be considered per se constraints on the parameters , since for each of them it is easy to realize that one of the two parameters will always satisfies the bound, as will do one of .
Something more can be said making partial use of the information provided by current modeling of the Crab nebula emission. In particular current reconstructions of the Crab emission fit very well the data by just assuming a single lepton population accounting for both the synchrotron and IC emission [11]. In this case it is possible to infer that at least one of the four pairs must lie in the narrow region bounded horizontally by the dashed lines of the synch and IC bounds and vertically by the birefringence constraint. However, we cannot a priori exclude that only one out of four populations is responsible for both the synchrotron and inverse Compton emission, so the electron sector is not yet strongly constrained (see again [11] for further details).
Nonetheless, it is clear that these simple arguments do not fully exploit the large amount of information we obtain from the Crab Nebula. A detailed comparison of the observations with the reconstructed spectrum in the LIV case, where all new reactions and modifications of classical processes are considered, could still provide us with strong constraints on both positive and negative for the four lepton populations. This study is currently in progress and there is reasonable expectation that it might settle this issue [68].
7 The naturalness problem: an analogue models lesson
Looking back at our ansatz Eq. (2) compatibility with observations induced us to assume that the lowest order coefficients and contain some appropriate power of the small ratio . Note however that we did not assumed any extra Planck suppression for the higher order dimensionless coefficients ( with ) which are then naturally of order one. Such an ansatz assures that for any the lowest higher order term allowed is always dominant on both the lowest oder ones as well as on all the other higher order ones. If CPT invariance is not a priori enforced, the dominant term would be the one cubic in the momentum in Eq. (2).
A naturalness problem arises because such a line of reasoning does not seem to be well justified within an EFT framework. In fact we implicitly assumed that there are no extra Planck suppressions hidden in the dimensionless coefficients with . However we cannot justify why only the dimensionless coefficients of the terms should be suppressed by powers of the small ratio . Furthermore it is easy to show [69] that, without some protecting symmetry, it is generic that radiative corrections due to particle interactions in an EFT with only Lorentz violations of order in (2) for the free particles, will generate and Lorentz violating terms in the dispersion relation which will then be dominant.
It has indeed been suggested in [70] that supersymmetry (SUSY) could play a protective role for the lowestorder operators: the dual requirements of supersymmetry and gauge invariance permit one to add to the SUSY standard model only those operators corresponding to terms in the dispersion relation. However SUSY is broken in the real world and when LV SUSY QED with softly broken SUSY was considered [71], it was found that, upon SUSY breaking, the dimension five SUSY operators generate dimension three operators large enough that the dimension five operators must be suppressed by a mass scale much greater than . In this sense, the naturalness problem is not completely solved in this setting. Alternatively one can try to identify LV operators that are protected against transmutation in lower order ones by a variety of other mechanisms like the irreducibility of the LV tensor structures,Tinvariance, and the lepton number conservation. This has been explicitly done for dimension 5 Lorentz Violating Interactions in the Standard Model showing that some LIV operators are indeed protected in this way [72] ^{8}^{8}8Renormalization group equations for QED with dimension five LIV terms have been also given in [72]. Solving such equations can be used to show that and are of order 1 at the TeV scale [68]. .
Finally also analogue models of gravity can be used as an inspiring tool in tackling the naturalness problem. We already saw the wave equation for phonons propagating on a Bose–Einstein condensate BEC can be described as that of a relativistic scalar field on a curved spacetime (the background is flat in the special case that the BEC is at rest and has constant density). Moreover the dispersion relations for such phonons (Eq. (3)) does indeed show high energy corrections. However, in order to address the naturalness problem one needs to consider a system where there are at least two kind of excitations/quasiparticles which moreover have to share the same Lorentzian geometry at low energies. In this sense a system of two coupled BEC is ideal because does satisfy this criteria when its microscopic quantities (background densities and couplings) are suitably tuned [17].
The basis of the model is an ultracold dilute atomic gas of bosons in two coupled singleparticle states and . There are three atomatom coupling constants, , , and , and an additional coupling that drives transitions between the two singleparticle states. Ignoring again back reaction effects of the quantum fluctuations one then obtains a pair of coupled Gross–Pitaevskii equations (GPE)
(7.0) 
where and again identify the classical wave function of each condensate specie.
Now consider small perturbations (sound waves) in the condensate cloud. The excitation spectrum is obtained by linearizing around some background, and after a straightforward analysis and a suitable tuning of the microscopic quantities it is possible to show that the dispersion relations at low energies are [17]
(7.0)  
(7.0) 
where and identified the two phonlike modes of the coupled system and where both the low energy speed of sound and the effective mass are a function of the microscopic parameters. In particular if the laser coupling .
One can then investigate if at high energies the interaction between these modes allows the primary quartic order (in particle momentum) LV to “percolate” to quadratic order (hence inducing nonzero terms) and make the latter the dominant contribution. In [17] it was shown that at high energies one indeed finds that the high energy dispersion relations become [17]
(7.0)  
(7.0) 
where plays the role of the LIV scale.
So that a percolation of the standard (for BEC) quartic LIV terms into quadratic ones seems to have taken place (at least for the massive quasiparticle). However looking in detail at the form of the and coefficients one finds that while the former one is actually a small ratio ^{9}^{9}9One can always tune the microscopic physics of the 2BEC system so that the latter are indeed coefficients of order one. So for any the higher order LIV terms will dominate over the quadratic one.
We can here propose a nice interpretation in terms of “emergent symmetry”: Nonzero simultaneously produces a nonzero mass for one of the phonons, and a corresponding nonzero LIV at order . Let us now drive but in such a ay that at low energies the Lorentzian geometry si still recovered (this can always be done [17]). In this case one gets an EFT which at low energies describes two noninteracting phonons propagating on a common background (both and go to 0 for ). Now this system possesses a symmetry and hence nonzero laser coupling softly breaks this , the mass degeneracy, and lowenergy Lorentz invariance. Such soft Lorentz violation is then characterized (as usual in EFT) by the ratio of the scale of the symmetry breaking , and that of the scale originating the LIV in first place . We stress that the symmetry is an “emergent/accidental symmetry” as it is not preserved beyond the hydrodynamic limit: the coefficients are in general different if , so is generically broken at high energies. Nevertheless this is enough for the protection of the lowestorder LIV operators.
The interesting lesson to be learn concerns the fact that in this case one can see that the symmetry that protects the lowest order operators (those associated with the term) is not an high energy one eventually broken at low energies, on the contrary it is an emergent (low energy) symmetry which being softly broken by interactions allows only for suppressed low order terms. It would be intriguing to explore the possibility that something along these lines could happen in the real world.
8 Conclusions
What we presented here is just a succinct account of an ongoing effort to constraints QG induced departures from Lorentz symmetry at high (but still subPlanckian) energies. It is remarkable that even without experiments specifically dedicated to the detection of LV effects in astrophysics so much could be done in constraining the EFT framework of [49]. In the future, experiments like AUGER, GLAST or OWL could help in improving our knowledge of the microscopic spacetime nature. There is also some hope that high energy neutrino astrophysics might play a pivotal role in the next decade advancing this branch of research. We direct again to [11] for further discussion. Moreover we wanted to stress here the interesting role which the interplay of different fields of physics — high energy physics, astrophysics, condensed matter physics — can have in this branch of research. All in all while a theory of quantum/emergent gravity is definitely in the realm of high energy physics, it is also true that our best option to test it is by describing its semiclassical limit via EFT scenarios, i.e. the “heart and soul” of condensed matter physics, and constrain their prediction using the highest energies observations so far available, i.e. the high energy astrophysics one. Who is writing strongly feels that the new challenges ahead us will not be overcome without further exploiting this rich mixing of ideas. In any case we hope that this brief account of the present achievements might have convinced the reader that the age of purely speculative QG models has ended and that the dawn of quantum gravity phenomenology has finally arrived.
Acknowledgments
References
 [1] V. F. Mukhanov, H. A. Feldman and R. H. Brandenberger, Phys. Rept. 215, 203 (1992).
 [2] C. D. Hoyle, D. J. Kapner, B. R. Heckel, E. G. Adelberger, J. H. Gundlach, U. Schmidt and H. E. Swanson, Phys. Rev. D 70 (2004) 042004
 [3] C. D. Hoyle, U. Schmidt, B. R. Heckel, E. G. Adelberger, J. H. Gundlach, D. J. Kapner and H. E. Swanson, Phys. Rev. Lett. 86 (2001) 1418
 [4] G. AmelinoCamelia, Phys. Rev. D 62, 024015 (2000)
 [5] S. Dimopoulos and G. Landsberg, Phys. Rev. Lett. 87 (2001) 161602
 [6] L. A. Anchordoqui, J. L. Feng, H. Goldberg and A. D. Shapere, Phys. Rev. D 65, 124027 (2002)
 [7] J. R. Ellis, J. L. Lopez, N. E. Mavromatos and D. V. Nanopoulos, Phys. Rev. D 53, 3846 (1996)
 [8] V. A. Kostelecky, “CPT and Lorentz symmetry”. Proceedings of the 3rd Meeting on CPT and Lorentz Symmetry (CPT 04), Bloomington, Indiana, 47 Aug 2004. Published in Hackensack, USA: World Scientific (2005).
 [9] L. J. Garay, Int. J. Mod. Phys. A 10, 145 (1995)
 [10] D. Mattingly, Living Rev. Rel. 8, 5 (2005)
 [11] T. Jacobson, S. Liberati and D. Mattingly, Annals Phys. 321, 150 (2006)
 [12] C. Barcelo, S. Liberati and M. Visser, “Analogue gravity,” Living Rev. Rel. 8, 12 (2005).
 [13] S. Liberati, S. Sonego and M. Visser, Phys. Rev. D 71, 045001 (2005).
 [14] R. Aloisio, A. Galante, A. Grillo, S. Liberati, E. Luzio and F. Mendez, Phys. Rev. D 73, 045020 (2006).
 [15] R. Aloisio, A. Galante, A. F. Grillo, S. Liberati, E. Luzio and F. Mendez, Phys. Rev. D 74, 085017 (2006)
 [16] F. Girelli, S. Liberati and L. Sindoni, Phys. Rev. D 75, 064015 (2007)
 [17] S. Liberati, M. Visser and S. Weinfurtner, Phys. Rev. Lett. 96, 151301 (2006). Class. Quant. Grav. 23, 3129 (2006).

[18]
C. J. Isham, “Structural issues in quantum gravity”,
[arXiv:grqc/9510063];
J. Butterfield and C. J. Isham, “On the emergence of time in quantum gravity”, [arXiv:grqc/9901024].  [19] K. V. Kuchar, “Time And Interpretations Of Quantum Gravity”, In “Winnipeg 1991, Proceedings, General relativity and relativistic astrophysics”, 211314.
 [20] L. J. Garay, Int. J. Mod. Phys. A 10, 145 (1995).
 [21] V. A. Kostelecky and S. Samuel, Phys. Rev. D 39, 683 (1989).

[22]
G. AmelinoCamelia, J. R. Ellis, N. E. Mavromatos and D. V. Nanopoulos,
Int. J. Mod. Phys. A 12, 607 (1997);
G. AmelinoCamelia, J. R. Ellis, N. E. Mavromatos, D. V. Nanopoulos and S. Sarkar, Nature 393, 763 (1998).  [23] R. Gambini and J. Pullin, Phys. Rev. D 59, 124021 (1999); J. Alfaro, H. A. MoralesTecotl and L. F. Urrutia, Phys. Rev. D 65, 103509 (2002).
 [24] S. M. Carroll, J. A. Harvey, V. A. Kostelecky, C. D. Lane and T. Okamoto, Phys. Rev. Lett. 87, 141601 (2001).
 [25] C. P. Burgess, J. Cline, E. Filotas, J. Matias and G. D. Moore, JHEP 0203, 043 (2002).
 [26] L. J. Garay, J. R. Anglin, J. I. Cirac and P. Zoller, Phys. Rev. Lett. 85, 4643 (2000)
 [27] L. J. Garay, J. R. Anglin, J. I. Cirac and P. Zoller, Phys. Rev. A 63, 023611 (2001)
 [28] C. Barcelo, S. Liberati and M. Visser, Class. Quant. Grav. 18, 1137 (2001)
 [29] M. Visser, C. Barcelo and S. Liberati, “Acoustics in BoseEinstein condensates as an example of broken Lorentz symmetry,” [arXiv:hepth/0109033].
 [30] C. Barcelo, S. Liberati and M. Visser, Phys. Rev. A 68, 053613 (2003) [arXiv:condmat/0307491].
 [31] N. Bogoliubov, J. Phys. (USSR) 11, 23 (1947).
 [32] G. AmelinoCamelia, Int. J. Mod. Phys. D 11, 35–59 (2002).

[33]
J. Magueijo and L. Smolin, Phys. Rev. Lett. 88, 190403 (2002). [arXiv:hepth/0112090];
scale,” Phys. Rev. D 67, 044017 (2003) [arXiv:grqc/0207085].  [34] G. AmelinoCamelia, Int. J. Mod. Phys. D 12, 1211–1226 (2003) [arXiv:astroph/0209232].
 [35] J. Lukierski and A. Nowicki, Int. J. Mod. Phys. A18, 7 (2003).
 [36] S. Judes and M. Visser, Phys. Rev. D 68, 045001 (2003).
 [37] L. Freidel, J. KowalskiGlikman and Lee Smolin, Phys. Rev. D69, 044001 (2004); J. KowalskiGlikman, 3rd International Sakharov Conference on Physics, Moscow, Russia, 2429 Jun 2002 [hepth/0209264].
 [38] J. Magueijo and L. Smolin, Class. Quant. Grav. 21, 1725 (2004).
 [39] F. Girelli, E. R. Livine and D. Oriti, Nucl. Phys. B 708, 411 (2005); F. Girelli and E. R. Livine, arXiv:grqc/0412004]; Braz. J. Phys. 35, 432 (2005).
 [40] D. Heyman, F. Hinteleitner and S. Major, Phys. Rev. D 69105016 (2004).
 [41] G. AmelinoCamelia, Int. J. Mod. Phys. D 11 (2002) 35 [arXiv:grqc/0012051].
 [42] J. KowalskiGlikman and S. Nowak, “Noncommutative spacetime of doubly special relativity theories,” Int. J. Mod. Phys. D 12 (2003) 299 [arXiv:hepth/0204245].
 [43] L. Freidel, J. KowalskiGlikman and L. Smolin, Phys. Rev. D 69, 044001 (2004) [arXiv:hepth/0307085].
 [44] G. AmelinoCamelia, G. Mandanici and K. Yoshida, “On the IR / UV mixing and experimental limits on the parameters of canonical noncommutative spacetimes,” JHEP 0401 (2004) 037 [arXiv:hepth/0209254].

[45]
C. Barcelo, S. Liberati and M. Visser,
Class. Quant. Grav. 19, 2961 (2002)
[arXiv:grqc/0111059];
S. Weinfurtner, S. Liberati and M. Visser, [arXiv:grqc/0605121].  [46] S. Mignemi, arXiv:0704.1728 [grqc].

[47]
D. Colladay and V. A. Kostelecky,
Phys. Rev. D 55, 6760 (1997);
Phys. Rev. D 58, 116002 (1998).  [48] S. R. Coleman and S. L. Glashow, Phys. Rev. D 59, 116008 (1999).
 [49] R. C. Myers and M. Pospelov, Phys. Rev. Lett. 90, 211601 (2003).
 [50] O. W. Greenberg, Phys. Rev. Lett. 89, 231602 (2002).
 [51] T. A. Jacobson, S. Liberati, D. Mattingly and F. W. Stecker, Phys. Rev. Lett. 93, 021101 (2004)
 [52] S. D. Biller et al., “Limits to quantum gravity effects from observations of TeV flares in active galaxies,” Phys. Rev. Lett. 83, 2108 (1999) [arXiv:grqc/9810044]; B. E. Schaefer, Phys. Rev. Lett. 82, 4964 (1999). P. Kaaret, “Pulsar radiation and quantum gravity,” Astronomy and Astrophysics, 345, L32L34 (1999) [arXiv:astroph/9903464]; S. E. Boggs, C. B. Wunderer, K. Hurley and W. Coburn, Astrophys. J. 611, L77 (2004).
 [53] J. R. Ellis, N. E. Mavromatos, D. V. Nanopoulos, A. S. Sakharov and E. K. G. Sarkisyan, Astropart. Phys. 25, 402 (2006).
 [54] R. J. Gleiser and C. N. Kozameh, Phys. Rev. D 64, 083007 (2001).
 [55] I. G. Mitrofanov, Nature 426, 139 (2003).
 [56] W. Coburn and S. E. Boggs, Nature 423, 415 (2003).
 [57] R. E. Rutledge and D. B. Fox, “ReAnalysis of Polarization in the Gammaray flux of GRB 021206,” [arXiv:astroph/0310385]; S. E. Boggs and W. Coburn, “Statistical Uncertainty in the ReAnalysis of Polarization in GRB021206,” [arXiv:astroph/0310515]. C. Wigger et al., “GammaRay Burst Polarization: Limits from RHESSI Measurements,” [arXiv:astroph/0405525].
 [58] Y. Z. Fan, D. M. Wei and D. Xu, “Gammaray Burst UV/optical afterglow polarimetry as a probe of Quantum Gravity,” [arXiv:astroph/0702006].
 [59] F. Aharonian et al. [The HEGRA Collaboration], Astrophys. J. 614, 897 (2004)
 [60] W. Kluzniak, Astropart. Phys. 11, 117 (1999).
 [61] R. Aloisio, P. Blasi, P. L. Ghia and A. F. Grillo, Phys. Rev. D 62, 053010 (2000).
 [62] G. AmelinoCamelia and T. Piran, Phys. Rev. D 64, 036005 (2001).
 [63] T. Jacobson, S. Liberati and D. Mattingly, Phys. Rev. D 66, 081302 (2002); Phys. Rev. D 67, 124011 (2003).
 [64] D. Mattingly, T. Jacobson and S. Liberati, Phys. Rev. D 67, 124012 (2003).
 [65] G. Gelmini, S. Nussinov and C. E. Yaguna, JCAP 0506, 012 (2005).
 [66] T. Jacobson, S. Liberati and D. Mattingly, Nature 424, 1019 (2003).
 [67] R. Montemayor and L. F. Urrutia, Phys. Rev. D 72, 045018 (2005) [arXiv:hepph/0505135].

[68]
L. Maccione, S. Liberati, A. Celotti, J. Kirk,
“The Crab Nebula constraints on Planckscale Lorentz Violation in QED”,
In preparation  [69] J. Collins, A. Perez, D. Sudarsky, L. Urrutia and H. Vucetich, Phys. Rev. Lett. 93, 191301 (2004).
 [70] S. G. Nibbelink and M. Pospelov, Phys. Rev. Lett. 94, 081601 (2005).
 [71] P. A. Bolokhov, S. G. Nibbelink and M. Pospelov, Phys. Rev. D 72, 015013 (2005).
 [72] P. A. Bolokhov and M. Pospelov, “Classification of dimension 5 Lorentz violating interactions in the standard model,” [arXiv:hepph/0703291].