Lagrangian reconstruction of Barrow holographic dark energy in interacting tachyon model

We consider a correspondence between the tachyon dark energy model and Barrow holographic dark energy (BHDE). The latter is a modified scenario based on the application of the holographic principle with Barrow entropy instead of the usual Bekenstein-Hawking one. We reconstruct the dynamics of the tachyon scalar field $T$ in a curved Friedmann-Robertson-Walker Universe both in the presence and absence of interactions between dark energy and matter. As a result, we show that the tachyon field exhibits a non-trivial dynamics. In a flat Universe, $\dot T^2$ must always be vanishing, independently of the existence of interaction. This implies $\omega_D=-1$ for the equation-of-state parameter, which in turn can be used for modeling the cosmological constant behavior. On the other hand, for a non-flat Universe and various values of Barrow parameter, we find that $\dot T^2$ decreases monotonically for increasing $\cos(R_h/a)$ and $\cosh(R_h/a)$, where $R_h$ and $a$ are the future event horizon and the scale factor, respectively. Specifically, $\dot T^2\ge0$ for a closed Universe, while $\dot T^2<0$ for an open one, which is physically not allowed. We finally comment on the inflation mechanism and Trans-Planckian Censorship Conjecture in BHDE and discuss observational consistency of our model.


I. INTRODUCTION
Experimental evidences from Supernova SNIa, Baryon acoustic oscillations and gravitational waves have definitely proved that our Universe is expanding at an accelerated rate.In spite of enormous effort, a fully consistent explanation for the origin of this behavior is missing.Among the various mechanisms, the existence of an unknown form of energy (Dark Energy, DE) affecting the Universe on large scales is the most widely accepted proposal.But yet, the nature of DE remains quite elusive.The possibility that DE is modeled by the cosmological constant acting as source of vacuum energy has been originally considered as natural way out of the DE puzzle [1,2].However, this scenario is at odds with our field theoretical understanding of the quantum properties of vacuum, thus requiring further investigation.Along this line, a plethora of DE models have been put forward over the years [3][4][5][6][7][8][9][10][11][12][13][14][15][16][17][18][19][20][21][22].
An interesting model to account for the nature of DE is the so called Holographic Dark Energy (HDE) [23][24][25][26][27][28][29], which emerges within quantum gravity framework.The main ingredient of this approach is the holographic principle, according to which the description of a volume of space can be thought of as encoded on a lower-dimensional boundary surface to the region.In [30,31] it has been pointed out that effective local quantum field theories over-count the number of independent degrees of freedom, predicting that entropy scales extensively (S ∼ L 3 ) for systems of size L with UV cutoff Λ.Later on, a solution to this problem has been provided in [23], where it has been argued that the total energy of a system with size L should not exceed that of an equally sized black hole, i.e.L 3 ρ Λ ≤ LM 2 p .Here, M p = (8πG) −1/2 is the reduced Planck mass, while ρ Λ denotes the quantum zero-point energy density caused by the UV cutoff Λ (we are working in natural units = c = 1).The inequality is saturated for the largest value of L. In this context, the holographic dark energy density is obtained as where c is a dimensionless constant and the factor 3 has been introduced for convenience.Cosmological applications of the holographic principle and HDE have been largely considered in literature.As an example, it was analyzed by [32] that the consequence of excluding those degrees of freedom of the system that will never be observed by the effective field theory results into an IR cutoff L at the future event horizon.In a DE dominated Universe, such an horizon is then predicted to tend toward a constant value of the order H −1 0 , with H 0 being the present Hubble rate [33].Furthermore, the issue of assuming the apparent (Hubble) horizon R A = 1/H as IR cutoff in a flat Universe has been examined in [34,35].Despite the intensive study, the shortcomings of the HDE in describing the history of a flat Friedmann-Robertson-Walker (FRW) Universe have prompted tentative changes to this approach.For instance, HDE has been used to address the DE problem in Brans-Dicke Cosmology [36][37][38][39][40][41][42] by considering different IR cutoffs [28,29,41] or and/or generalized entropies [43][44][45][46][47][48][49][50][51][52][53][54][55][56][57][58].In particular, the latter path has led to promising models, such as Tsallis [43][44][45][46]59], Barrow [47,51,53] and Kaniadakis [49,50,52] holographic dark energy, the latter two being motivated by quantum gravitational and relativistic considerations, respectively.
While exhibiting a richer phenomenology comparing to standard Cosmology, HDE with generalized entropies is often based on ad hoc deformations of the entropy-area law, which might somehow question the relevance of this model in improving our knowledge of the Universe at very fundamental level.In the absence of solid empirical guidelines, valuable hints can be gained by looking at Noether symmetries of the underlying theoretical framework, which are notoriously linked to physical conserved quantities (see also [60,61] for other recent usage of Noether symmetries in cosmology and gravitation).In order to export Noether's theorem to the present analysis, a reformulation of extended HDE is needed using the Lagrangian language.Toward this end, a pathway is to consider reconstruction scenarios, i.e. to compare the relative energy density of extended HDE and other solid cosmological models to find the reconstructed action which reproduces the whole cosmic history of the Universe.Reconstruction paradigm could be particularly inspiring for Barrow Holographic Dark Energy as a preliminary attempt to formulate the effective action of cosmological model in a quantum gravity-oriented picture.
Along this line and motivated by the analysis of [33], in [62] it has been shown that Tsallis holographic description of DE (THDE) is non-trivially intertwined with tachyon dark energy model [62].In this regard, we would like to remark that the tachyon field has been proposed as a possible candidate for dark energy.In particular, a rolling tachyon has a peculiar equation of state parameter that interpolates between the values −1 and 0 [63].Therefore, the tachyon can be realized as a suitable candidate for the high energy inflation [64] and at the same time as a source of dark energy depending on the form of the tachyon potential [65].Not least, such a model is mathematically less cumbersome than scalar-tensor theory, thus allowing for a more direct and intuitive interpretation of results.In [62] a correspondence between the tachyon field and THDE has been established based on the reconstruction of the dynamics of the tachyon field in THDE.In recent years, more applications have been analyzed in f (R) [66], f (R, T ) [67], f (G, T ) [68], teleparallel [69], Brans-Dicke [70], logarithmic Brans-Dicke [71] and Saez-Ballester [72,73] theories, among others.
Starting from the above premises, in this work we explore more in-depth the connection between the tachyon dark energy model and HDE.We frame our analysis in the context of HDE based on Barrow entropy [47].The ensuing scenario is typically named Barrow Holographic Dark Energy (BHDE) and arises from the application of the holographic principle at a cosmological framework, but employing Barrow entropy [74] instead of the standard Bekenstein-Hawking one.We analyze the case of a non-flat FRW Universe for interacting dark energy.Since scalar fields are generally conjectured to have driven inflation in the very early Universe, we then study the inflation mechanism in our BHDE model.We find an analytical solution for the slow-roll parameters, the scalar spectral index and the tensor-to-scalar ratio.We also compare our findings with recent results in the literature.
The remainder of the work is structured as follows: in the next Section we introduce BHDE.Section III is devoted to analyze the correspondence between the tachyon dark energy and BHDE in a non-flat FRW Universe.In Sec.IV we discuss inflation in BHDE, while conclusions and outlook are summarized in Sec.V.

II. BARROW HOLOGRAPHIC DARK ENERGY
Let us briefly review the basics of BHDE.We consider the four-dimensional Friedmann-Robertson-Walker (FRW) metric of scale factor a(t) and spatial curvature k = 0, 1, −1 for a flat, closed and open Universe, respectively.We use the definition (1) for the holographic dark energy in standard Cosmology and assume [33] where r(t) is the (time-dependent) radius that is relevant to the future event horizon of the Universe.Since we easily obtain where R h is the future event horizon given by [33] HDE relies on the holographic principle, which asserts that the number of degrees of freedom describing the physics of any quantum gravity system i ) scales as the bounding surface (rather than the volume) of the system and ii ) should be constrained by an infrared cutoff [30,31].This is in tune with Bekenstein-Hawking (BH) relation S BH = A/A 0 for black holes, where S BH and A denote the entropy and area of the black hole, respectively, while A 0 = 4G its Planck area.Recently, deformations of this relation have been proposed to take account of quantum [74][75][76] and/or relativistic [77] effects.In particular, in [74] it has been argued that quantum gravity may introduce intricate, fractal features on the black hole horizon, leading to the modified area law Deviations from BH entropy are quantified by the exponent 0 ≤ ∆ ≤ 1, with ∆ = 0 giving the BH limit, while ∆ = 1 corresponding to the maximal horizon deformation.We emphasize that although this relation resembles Tsallis entropy in non-extensive statistical thermodynamics [75,76], the origin and motivation underlying Eq. ( 7) are completely different.Cosmological implications of Barrow entropy have been recently studied in the context of Big Bang Nucleosynthesis [78], Baryogenesis [79] and tests of gravity theories from observations of Sagittarius A* [80], among others.The possibility of a running ∆ has also been considered in [81].Strictly speaking, Eq. ( 7) has been formulated for black holes.However, it is known that in any gravity theory one can consider the entropy for the Universe horizon in the same form as the black hole entropy, the only adjustment being the replacement of the black hole horizon radius with the apparent horizon radius.This is at the heart of the various generalizations of HDE with modified entropy laws (see, e.g.[43,44,[47][48][49]). Now, in [23] Cohen et al. have proposed the following inequality between the entropy, the IR (L) and UV (Λ) cutoffs for a given system in an effective local quantum field theory If we use for the entropy the modified expression (7), we have where Λ 4 denotes the vacuum energy density, i.e. the energy density of DE (ρ D ) in the HDE hypothesis [82].By using the above inequality, Barrow holographic dark energy density can be proposed as where C is an unknown parameter with dimensions [L] −2−∆ .It is worth noticing that for ∆ = 0, the above relation reduces to the standard HDE (1), provided that C = 3c 2 M 2 p .On the other hand, in the case where deformations effects switch on (∆ = 0), BHDE departs from the standard HDE, leading to different cosmological scenarios [47].
Following the standard literature, we now define the critical energy density ρ cr and the curvature energy density ρ k as We also introduce the three fractional energy densities where ρ m is the matter energy density.In particular, by setting L = H −1 as in [35,43,[83][84][85][86][87] we obtain 1 From Eqs. ( 3) and ( 4), one can derive the following expression for the time derivative of L [62] where we have defined [33] 1 Now, for a flat FRW Universe filled by non-interacting BHDE and pressureless DM, the first Friedmann equation takes the form which, by use of Eqs. ( 12) and ( 13), can be rewritten as where Since BHDE does not interact with other parts of cosmos (DM), the conservation equations of dust and THDE read 1 There are several choices for the IR cutoff L. Following [43], here we resort to the simplest one L = H −1 .Other possible choices are the particle horizon, the future event horizon, the GO cutoff [88] or combination thereof.However, in these cases one must generally resort to numerical evaluation to study the cosmological evolution of the model [44].Since we are interested in extracting analytical solutions and given the degree of arbitrariness in the selection of the most reliable description of dark energy, we leave the analysis of dark energy models with different IR cutoffs for future investigation.
where we have denoted by ω D = p D /ρ D and p D the equation of state parameter and pressure of THDE, respectively.From Eq. ( 21), we obtain ρ m = ρ m,0 /a 3 , where ρ m,0 is the present matter energy density.Deriving Eq. ( 18) respect to time and using the continuity equations ( 21) and ( 22), after some algebra we are led to Likewise, by plugging Eq. ( 10) into (22), we find which gives, by comparison with Eq. ( 23) With the aid of Eq. ( 19), this finally yields On the other hand, if there exists an interaction of the type between BHDE and matter, the continuity equations ( 21) and ( 22) become Following similar calculations as above, one can show that Barrow holographic energy equation of state takes the form where b is the coupling parameter that quantifies the interaction.

A. The age of the Universe
Let us now consider the following integral where Ω D,i ≡ Ω D (t = 0) and Ω D,0 is the current value of the DE density.Here, we have used Eqs.( 24) and ( 26) along with By integrating the above relation, it follows that where 2 F 1 (a, b; c; d) is the hypergeometric function of first kind.This equation can be used to estimate the order of the age of the current Universe (z = 0) in our model.Specifically, by using the relation between the Hubble and EoS parameters, we can approximately write For ω D (z = 0) = −2/3, we then have t = 1/H 0 for ∆ = 1, corresponding to the maximal deformation of the Bekenstein-Hawking area law.As observed in [43], further corrections to Eq. ( 34) may arise due to either different modifications of the horizon entropy or other IR cutoffs.

III. TACHYON SCALAR FIELD AS BARROW HOLOGRAPHIC DARK ENERGY IN A NON-FLAT FRW UNIVERSE
In this Section we analyze the correspondence between the tachyon dark energy model and the BHDE scenario in a non-flat FRW Universe.Toward this end, we recall that in [33] it has been shown that the energy density ρ T and pressure p T for the tachyon scalar field take the form where V (T ) is the tachyon potential energy.From these relations, we derive the equation of state parameter (EoS) for the tachyon as which will be later equated to the EoS parameter of BHDE in our reconstructed scenario.We also remind that the cosmological model based on the effective Lagrangian of tachyonic matter with V (T ) = const.coincides with Chaplygin gas model3 [89].
Let us explicitly explore the connection between BHDE and tachyon dark energy model.Toward this end, we consider the time derivative of BHDE (10) and use Eq. ( 16) to get where x ≡ R h /a.By means of the continuity equation ( 29), this can be cast as We can now resort to Eq. ( 37) to obtain where we have made use of the required correspondence ω T = ω D , as explained above.This relation can be further manipulated by dividing both sides by 3Hρ D and using Eq. ( 20) to give After employing Eq. ( 10) and the condition L = H −1 , we finally reach which can be equivalently written as Notice that, once known Ṫ 2 , one can express the tachyon potential V (T ) in terms of H and the parameter ∆ by simply equating Eqs. ( 10) and (35).
Let us now study the cases k = 0, 1, −1 separately.The k = 0 framework corresponds to a flat FRW Universe.In this case it is easy to show that where we have used Eqs.( 17) and (20).Thus, if we require Ṫ to be real, then b = 0 in flat space, which allows us to conclude that Ṫ 2 must always be vanishing in flat space, independently of the existence of interaction.In turn, from Eqs. (35) and (36) this implies We conclude that in this case, the equation of state is always ω D = −1, reproducing a cosmological constant-like behavior.On the other hand, for k = 1 (closed FRW Universe) we get from the definition ( 17) In order for Ṫ to be zero (i.e.ω T = −1), we must have which clearly admits non-trivial solution, since ∆ < 1.This provides the condition for which the tachyon model of BHDE reproduces the cosmological constant behavior in a closed FRW Universe.By contrast, in [62] it is argued that Ṫ 2 cannot be zero in Tsallis holographic dark energy in a non-flat Universe.The evolution trajectory of Ṫ 2 in Eq. ( 47) is plotted in Fig. 1 for fixed b and u and various allowed ∆.One can see that Ṫ 2 decreases monotonically for increasing cos(x) and Ṫ 2 ≥ 0, so that ω T ≥ −1 (quintessence or cosmological constant-like behavior), while tending to negative values for cos(x) → 1, which becomes physically invalid.The non-interacting case can be simply derived by setting b 2 = 0 in Eqs.(47) and (48).The ensuing behavior is similar to that described above, the only difference being that Ṫ 2 ≥ 0 throughout the whole evolution in this case.
Similarly, we can consider the dynamics of the tachyon field in an open (k = −1) Universe.Following the same reasoning as above, we obtain The evolution of Ṫ 2 in Eq. ( 49) is plotted in Fig. 2 for different values of ∆.As before, we notice that Ṫ 2 decreases monotonically for increasing cosh(x), but in this case it is always negative, which is not a physically valid situation.We also see that Eq. ( 49) vanishes, provided that However, since cosh(x) > 1 (we remind that x = 0), we infer Ṫ 2 can never be zero, in agreement with the result of [62].The same behavior is exhibited in the absence of interactions (b 2 = 0).

A. Observational studies
This section is devoted to explore some observational implications of BHDE.For simplicity, we focus on the case where there is no interaction between the dark sectors of the cosmos and positive values of curvature.First, we notice that the use of Eq. ( 23) in non-flat Universe allows us to write down the deceleration parameter q as From the above definition, one can check that q > 0 (q < 0) corresponds to a decelerated (accelerated) expansion of the Universe, since ä < 0 (ä > 0).By computing the time derivative of Eq. ( 10), we then get while the further usage of the definition (12) leads to where . By plugging into (22), the BHDE EoS parameter and the fractional BHDE density take the form The latter equation is solved numerically and plotted in Fig. 3, which shows a monotonic increasing of dark energy for decreasing redshift, implying a DE dominated Universe in the far future.On the other hand, the evolution of the EoS parameter is displayed in Fig. 4, which indicates that BHDE behaves like quintessence or (asymptotically) cosmological constant, consistently with the discussion below Fig. 2. Specifically, for the considered ∆'s, we can see that the present value of this parameter lies in the range −0.7 ω D,0 −0.4, which only slightly deviates from recent Planck+WP+BAO predictions [90].
Similarly, from Eq. ( 51) we can derive the evolution of the deceleration parameter, whose behavior is plotted in Fig. 5.We observe that the present model allows to explain the sequence of an early decelerating (i.e.q > 0) expansion of the Universe, followed by an accelerated (i.e.q < 0) epoch.Also, it is consistent with the description of the current acceleration, although the predicted value −0.3 q 0 −0.01 is slightly higher than q 0 −0.5 obtained in the standard ΛCDM model [90].
FIG. 4: Evolution trajectories of ωD versus z (we considered the same initial conditions as in Fig. 3).The dashed vertical line marks the value at present time.FIG.5: Evolution trajectories of q versus z (we considered the same initial conditions as in Fig. 3).The dashed vertical line marks the value at present time.
To further check the phenomenological consistency of our model, one could also study the evolution of some other relevant quantities, such as the Hubble rate or the effects of Barrow model on the growth of cosmological perturbations and structure formation.The analysis of these aspects is more challenging and requires further effort.Results will be presented in a future work.

IV. INFLATION IN BARROW HOLOGRAPHIC DARK ENERGY
In this Section we discuss inflation in BHDE.For reasons that will appear clear below and following [91], here we consider the more general expression for the length scale L −2 = αH 2 + β Ḣ, where α and β are dimensionless constant.Assuming that the expansion of the Universe is driven by BHDE (10) and neglecting the matter contribution due to the rapid inflationary expansion, Eq. ( 18) becomes from which we infer From this relation, it is clear that setting the IR cutoff L H −1 (i.e.β = 0) as in the previous study would give rise to technical issues in the present framework.
To simplify the resolution of Eq. ( 57), we introduce the e-folds variable N = log (a/a i ), where a i is the initial value of the scale factor a. By observing that dN = Hdt and Ḣ = 1 2 dH 2 dN , integration of Eq. ( 57) gives where H = H/M p is the dimensionless Hubble parameter and Here we have denoted the Hubble parameter at the end of inflation by Hf .From Eq. ( 56) we can now compute the characteristic parameters of slow-roll inflation.Specifically, the first slow-roll parameter is given by The other slow-roll parameters can be derived by using the definition n+1 = d log( n )/dN .For the second parameter Let us now evaluate the Hubble parameter at the end of inflation.This phase is characterized by 1 = 1.By straightforward calculations, we obtain On the other hand, at the beginning of inflation (including the horizon crossing time) Eq. ( 58) gives which can be used to calculate the slow-roll parameters for earlier time by direct substitution in Eqs. ( 60) and (61).
In order to derive the scalar spectral index n s − 1 and the tensor-to-scalar ratio r, we follow [91,92] and make use of the usual perturbation procedure.We are led to Clearly, a full perturbation analysis is needed to obtain the exact expressions of n s − 1 and r.Two comments are in order here: first, we notice that the constant γ does not intervene in the calculation of the slow-roll parameters at the horizon crossing time, which means that neither n s − 1 nor r depend on it.As explained in [91], this constant can be estimated by considering the amplitude of the scalar perturbation.Furthermore, it is worth mentioning that a similar analysis of inflation and correspondence between BHDE and tachyon field has been proposed in [93].However, in that case the authors consider values of Barrow parameter ∆ higher than unity, which are actually forbidden in Barrow model.This somehow questions the results exhibited in [93].

A. Trans-Planckian Censorship Conjecture
The large-scale structures we currently see in the Universe originated from matter and energy quantum fluctuations produced during inflation.Such fluctuations cross the Hubble radius during the early phase, are stretched out and classicalize, and finally re-enter the Hubble horizon to produce the CMB anisotropies.The key point is that if inflation lasted longer than the supposed minimal period, then it would be possible to observe length scales originated from modes smaller than the Planck length at inflation [94].This problem is usually referred to as "Trans-Planckian problem".To avoid inconsistencies, it has been conjectured that this problem cannot arise in any consistent model of quantum gravity ("Trans-Planckian Censorship Conjecture", TCC) [95].
The TCC states that no length scales which cross the Hubble horizon could ever have had a wavelength smaller than the Planck length.This is imposed by requiring that where L p = 1/M p is the Planck length and we have denoted by a f the scale factor at the end of inflation.By using Eq. ( 62) for the Hubble parameter at the final time, the TCC (65) becomes the validity of which can be examined by comparison with observational data.This aspect is under active investigation and will be addressed in more detail in a future work.

V. CONCLUSIONS AND OUTLOOK
The origin of the accelerated expansion of the Universe is an open problem in modern Cosmology.To date, the most reliable explanation is provided by the existence of an enigmatic form of energy -the Dark Energy -affecting the Universe on large scales.Several candidates have been considered to account for this phenomenon.In particular, the holographic dark energy has been largely studied, also in connection with different real scalar field theories, such as quintessence [3][4][5][6], K-essence [7], phantom [9][10][11], interacting models [16] and tachyon [33] (see also [96,97] for arguments that theoretically rule out quintessence and K-essence).Recently, the interest has been extended to the Tsallis holographic dark energy [43,44,[47][48][49] and the possibility of using it to describe the dynamics of the tachyon field [62].
In this work we have considered the further scenario of tachyon model as Barrow holographic dark energy.Barrow entropy arises from the effort to include quantum gravity effects on the black hole horizon.In this sense, the present analysis must be intended as a preliminary step toward a fully quantum gravity extension of [62].In the absence of empirical guidelines, we have exploited the powerful tool of Lagrangian formalism and, in perspective, the precious clues that may be provided by the related Noether's theorem.In particular, we have established a correspondence between BHDE and the tachyon field model in a FRW Universe, both in the presence and absence of interactions between dark energy and matter.In this regard, we would like to stress that the aim of the present analysis is to reconstruct BHDE in tachyon model, rather than showing the consistency between the two frameworks.In other terms, in the absence of a lagrangian formulation of BHDE, the question we have asked ourselves is how BHDE would appear if one requires that the evolution of its energy density can be described in terms of that of tachyon field.As a result, we have shown that the tachyon field should exhibit a non-trivial dynamics.In particular, in a flat Universe, Ṫ 2 must always be vanishing, independently of the existence of interaction, which implies ω D = −1 for the equation-of-state parameter.On the other hand, for a non-flat Universe and various values of Barrow parameter, we have found that Ṫ 2 decreases monotonically for increasing cos(R h /a) and cosh(R h /a).Specifically, Ṫ 2 ≥ 0 for a closed Universe, while Ṫ 2 < 0 for an open one, which is physically not allowed.We have finally investigated an inflationary scenario described by a Universe filled with BHDE and commented on the Trans-Planckian Censorship Conjecture.
We have shown that the tachyon field exhibits a non-trivial dynamics.In a flat Universe, Ṫ 2 must always be vanishing, independently of the existence of interaction, which implies ω D = −1 for the equation-of-state parameter.On the other hand, for a non-flat Universe and various values of Barrow parameter, we have found that Ṫ 2 decreases monotonically for increasing cos(R h /a) and cosh(R h /a).Specifically, Ṫ 2 ≥ 0 for a closed Universe, while Ṫ 2 < 0 for an open one, which is physically not allowed.We have finally discussed observational consistency of our model and investigated an inflationary scenario described by a Universe filled with BHDE, and commented on the Trans-Planckian Censorship Conjecture.
Further aspects remain to be addressed.For instance, we can look at the correspondence between the tachyon field and other dark energy scenarios, in particular stable dark energy models.Furthermore, it would be interesting to disclose the effects of Barrow entropic corrections on the growth of perturbations and structure formation in the present model.Preliminary studies along this direction have been recently conducted in [98].One more suggestive perspective concerns the extension of the above framework by using different IR cutoffs (such as the future event horizon or Grand-Oliveros cutoff) and/or other deformed entropies, such as Kaniadakis entropy [49], which is based on a relativistic self-consistent generalization of the classical Boltzmann-Gibbs entropy [77] (see also [99] for a recent review of Kaniadakis entropy applications in gravity and cosmology).Finally, it is essential to examine to what extent our effective model reconciles with predictions of more fundamental candidate theories of quantum gravity, such as String Theory and Loop Quantum Gravity, or more phenomenological approaches, such as generalizations of the Heisenberg relation at Planck scale [100][101][102][103]. Work along these directions requires further investigation and will be presented elsewhere.