Linear-Quadratic GUP and Thermodynamic Dimensional Reduction
Abstract
In this paper we investigate the statistical mechanics within the Linear-Quadratic GUP (LQGUP, i.e, GUP with linear and quadratic terms in momentum) models in the semiclassical regime. Then, some thermodynamic properties of a system of 3-dimensional harmonic oscillators are investigated by calculating the deformed partition functions. According to the equipartition theorem, we show that the number of accessible microstates decreases sharply in the very high temperatures regime. When the thermal de Broglie wavelength is of the order of the Planck length, three degrees of freedom are frozen in this setup. In other words, it is observed that there is an effective reduction of the degrees of freedom from to for a system of 3D harmonic oscillators in this framework. The calculations are carried out using both approximate analytical and exact numerical methods. The results of the analytical method are also presented in the form of thermal wavelengths for better understanding. Finally, the case of a 2-dimensional harmonic is treated as another example to comprehend the results, leading to a reduction of the degrees of freedom from to .
- PACS numbers
-
04.60.Bc
- Key Words
-
Quantum Gravity Phenomenology, Thermodynamics, Generalized Uncertainty Principle, Dimensional Reduction
1 Introduction
Since Planck announced his solution to the black body radiation, a fundamental and very important quantity, namely the Planck constant () entered physics. Today, this fundamental constant plays a very crucial and consequential role in science, especially in theoretical high energy physics [1, 2, 3]. For example, in quantum gravity phenomenology theories (in fact, approaches to quantum gravity proposal) such as string theory, loop quantum gravity, doubly special relativity (DSR), and also Noncommutative geometry, new quantities were introduced into physics that represent Planck’s constant dependence of the structure of spacetime at high energy level regime [4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19]. One of these new quantities is the fundamental length scale called the Planck Length , which is defined by combining with the gravitational constant and the constant speed of light , that is ( m). Planck mass, that is, kg, is another such a quantity. Combining quantum field theory with gravity leads naturally to an effective cutoff (a minimal measurable length preferably of the order of the Planck length) in the ultraviolet regime. Indeed, the high energies used to probe small distances significantly disturb the space-time structure by their strong gravitational effects. Thus one common feature of all approaches to Quantum Gravity (QG) proposal is the existence of a minimal measurable length. The existence of a minimal measurable length modifies the Heisenberg uncertainty principle (HUP) to the so called generalized uncertainty principle (GUP) (and also extended uncertainty principle (EUP)). GUP can be understood in the context of string theory and supports the existence of a minimal length/maximal momentum (energy) (UV cutoff) [20, 21, 22, 23, 24, 25].
In the HUP framework, there isn’t any essential limitation on the accuracy of measuring the position of test particles. Therefore, the minimal uncertainty in position measurement () can be arbitrarily reduced even to a value of zero. On the other hand, there is a limitation in the GUP framework due to the minimal uncertainty in position measurement.
This mentioned phenomenological model, as a flat limit of quantum gravity proposal, represents the deformation for the density of states of the system at very high temperature regimes, which typically leads to self-exclusive measurements over microstates. This property leads to the reduction of the degrees of freedom which may reflect in essence a reduction of the dimensions of the space for the desired system [26, 27, 28, 29, 30]. For a recent and elegant work on dimensional reduction in quantum gravity, see [31, 32]. There are several versions of GUP. The first one, quadratic GUP, modifies the standard Heisenberg Uncertainty Principle (HUP) by including an additional quadratic term in momentum. This suggests that gravity may exhibit different behavior at the minimal length scale than it does under general relativity [33]. Another version, the linear and quadratic GUP (LQGUP ), incorporates both the maximum momentum and minimum length, includes both linear and quadratic terms for momentum, and is compatible with Doubly Special Relativity (DSR) theories [34, 35, 36, 37, 38, 39, 40].
In this interesting and challenging streamline, we are going to study the issue of thermodynamic dimensional reduction through the statistical mechanics of a -D harmonic oscillator system within a version of the GUP model; the Linear-quadratic GUP, in the semiclassical regime. We use the deformed partition functions and study some thermodynamic properties of a system of three-dimensional harmonic oscillators. We’ll see that the number of microstates decreases strongly in the high temperature regime (or equivalently at the low thermal de Broglie wavelength), and this may be effectively related to the reduction of the space dimensions at high energy regime. As an important result, a trace of the fractional degrees of freedom (and possibly space dimensions) will be observed in this setup for a system of -D harmonic oscillators.
In this paper, in Section 2, the statistical mechanics of the LQGUP framework is studied and the corresponding partition function is formulated. In section 3, the thermodynamics of a system of D harmonic oscillators is studied analytically (approximated solutions) and also numerically (exact solutions) within the framework of LQGUP. Section 4 provides an analysis of a 2-dimensional harmonic oscillator thermodynamics in LQGUP setup. Section 5 is devoted to the summary and conclusions.
2 Statistical Mechanics in LQGUP Phase Space
The classical system’s kinematics and dynamics on the phase space offer a framework for formulating statistical mechanics in the semiclassical regime. The Liouville volume is the crucial factor that determines the density of states, which can be utilized to derive all of a system’s thermodynamic properties.
2.1 Kinematics and Dynamics
The generalized uncertainty relation that implies the existence of a nonzero minimal uncertainty in position measurement has the following form
(1) |
Inspired by doubly special relativity (DSR), there exists a limit on a test particle’s momentum, which prevents it from being imprecise. Consequently, this may result in a maximum measurable momentum [17, 41, 42]. Within this framework, the LQGUP that forecasts a minimum observable length and a maximum momentum can be expressed as follows, [34, 43, 44, 45, 46, 47],
(2) |
We set For simplicity ([22, 23, 24, 25]). It should be mentioned that the constants and are related through dimensional analysis as [] = []. In the given relation, represents the LQGUP parameter, which supports both the minimal length and the maximal momentum, and it’s defined as that GeV is the 4-dimensional fundamental scale. is a dimensionless constant that can be determined through experimentats [48] and it’s close to the unity, that is, . This modified relation suggests that the uncertainty in position measurement, , is always greater than , and the uncertainty in momentum measurement, , is always smaller than [24]. If tends to zero (), the Heisenberg uncertainty relation (HUP) is restored. Using this approach, we can derive the following algebraic structure for position and momentum operators
(3) |
(4) |
In ordinary quantum mechanics, the coordinates and momenta satisfy the Heisenberg algebra , whose classical counterpart is the Poisson algebra . In the LQGUP framework, these relations are written as
(5) |
(6) |
The Darboux theorem states that it is always possible to find variables and in terms of and such that the following relation holds
(7) |
Therefore, the Eq. (5) is a specific instance of the Darboux theorem and in the LQGUP framework . On the other hand, Eq. (6) which ensures via the Jacobi identity does not constrain our study. Therefore, we can conclude that the Eqs. (5) and (6) are generally valid simultaneously in the presence of space noncommutativity.
The Jacobian of transformation for 2D-dimensional LQGUP phase space is as follows [34].
(8) |
with the Jacobian of the transformation (8) in hand, one is led to the LQGUP Liouville volume as
(9) |
Where represents the canonical Liouville volume, we assume that the system consisting of particles is kinematically separated.
2.2 Number of Microstates
In the LQGUP phase space, the microstate count is affected by the modified commutation relations between position and momentum operators. These deformations introduce noncommutativity in the phase space, leading to modifications in the distribution of microstates and consequently impacting thermodynamic properties. Understanding this qualitative aspect is useful before delving into the quantitative analysis provided by the partition function.
The number of microstates for a phase space containing a particle is given as follows111Throughout this article, we work in units , where and are the Boltzmann constant and speed of light in vacuum, respectively. Also, the figures are plotted in units of and the mass is taken to be in appropriate units..
(10) |
We consider the domain where momentum is restricted, i.e. where (with and being constant). In this scenario, it is reasonable to assume that the number of microstates is finite. As per equation (10), the number of microstates for a non-deformed phase space is , where represents the spatial volume of the system. In the context of the LQGUP phase space, when substituting the Liouville volume (9) into equation (10), it is formulated as
(11) |
In the low temperature regime, i.e., (), equation (11) behaves as for the number of microstates. This aligns with the standard result of no deformation in this regime. In the high temperature regime (), the number of microstates (11) is independent of and behaves as indicating a three-order reduction. This is illustrated in Fig. 1, where the number of microstates (11) is plotted versus and compared to the standard state without deformation. This result demonstrates that in the high-energy regime where quantum gravity effects (minimum length and maximal momentum) prevail, the number of microstates will significantly decrease. Even though this outcome is qualitatively derived for a generic system (without specifying a Hamiltonian function and ensemble density) in this subsection, we will clarify this outcome explicitly in the following section for the specific scenario of a 3D-harmonic oscillator in canonical ensemble based on the renowned equipartition theorem of energy.

2.3 Partition Function
The Hamiltonian equations for each of the phase space coordinates are and , where is the Hamiltonian and determines the dynamics of the whole system. Hence, for a system such as a three-dimensional harmonic oscillator, the partition function for the LQGUP framework is given by
(13) |
in which the Gibbs factor is omitted, as it is assumed that the oscillators are localized. As evident from the relation (13), the fundamental volume element of the phase space, , is larger than the usual fundamental volume element for large momentum, where is the Planck constant. Now, let’s assume that the particles in the system are also dynamically decoupled, represented by , where and denote the position and momentum of the -th particle. In this case, the total partition function for this type of system is also decoupled as , where
(15) |
where is the total volume of phase space.
3 Thermodynamics of a System of 3D Harmonic Oscillators in LQGUP Phase Space
In this section, we study the thermodynamics of a system of 3-dimensional harmonic oscillators in the LQGUP framework, which accounts for minimal length and maximal momentum effects. The Hamiltonian function for a single harmonic oscillator is then , and by substituting it in Eq. (15), the corresponding partition function for a single particle in semi-classical regime can be derived as follows
(16) |
The partition function obtained above cannot be solved analytically. Therefore, we must rely on numerical calculation methods to determine this function and other related thermodynamic variables. However, we can approximate the answer by expanding the deformed phase space volume as follows
(17) |
Thus, the partition function (16) can be approximated as follows
(18) |
After integration of the above relation, we’ll have
(19) |
In this relation, “erf“ is the error function, defined by: . To better understand the qualitative features of the partition function (19), it is convenient to rewrite it in terms of the thermal de Broglie wavelength , defined by: as
(20) |
Where is the characteristic length of the harmonic oscillator system. In this formula, is the thermal de Broglie wavelength in the LQGUP framework, which is of the order of the Planck scale thermal de Broglie wavelength,
(21) |
By expanding the partition function (20) for very high and low temperature regimes, we find
(22) |
where is the non-deformed partition function for a three
dimensional single-harmonic oscillator. As observed from the above relation, the influence of quantum gravity (minimal length and maximal momentum) can be disregarded at the low temperature limit , and the classical outcome for the harmonic oscillator partition function is restored. Additionally, it is evident from Eq.(22) that three degrees of freedom will freeze () at the Planck length scale in the LQGUP framework, and the partition function behaves as if it were effectively a partition function of a D harmonic oscillator. Therefore, it is possible to conclude an effective reduction in the degrees of freedom from to and therefore in the space dimensions from to in the GUP framework for harmonic oscillator systems. This feature evidently signals the existence of fractional degrees of freedom and possibly space dimensions for such a system [49, 50].
It is important to note that in these processes of decreasing number of degrees of freedom, at intermediate temperatures in all the figures presented in this paper, we observe a fractional nature of the number of the degrees of freedom. In fact, it would be more accurate to state that the number of the degrees of freedom of the deformed models mentioned are not completely integers at high temperatures; rather, they are indeed fractional (see also figure 3).
From Eq.(20), the total partition function of the system of harmonic oscillators in the LQGUP framework can be expressed as follows
(23) |
Now, utilizing Eq. , we can derive all thermodynamic quantities using the standard definitions. First, we begin with the Helmholtz free energy , which is expressed as follows
(24) |
Having equipped with Helmholtz free energy, we can obtain other thermodynamic quantities as follows.
3.1 Internal energy
With the modified Helmholtz free energy as given by Eq. (24), the internal energy of the system can be modified as follows
(27) |
For a standard 3D-harmonic oscillator, the internal energy linearly depends on the temperature through the well-known relation . Any high energy scale is accessible for the system just by a sufficiently large increment in temperature. In figure (2), the internal energy diagram versus temperature (approximated and exact solutions) and also versus thermal de Broglie wavelength in the LQGUP framework are shown. As can be seen from the figure, as the temperature increases, the deviation from the classical state (quantum gravity effects that appear as a minimal length and maximal momentum) becomes more and more appreciable. To expand the relation (27) for both of the high and low temperature regimes, we have.
(28) |
As the above relations clearly show, the system has a reduction of the degrees of freedom in the LQGUP phase space model. In fact, the degrees of freedom of the system, i.e. are here reduced from to . The degrees of freedom are plotted in figure (3). As one can see in this figure, the reduction of the degrees of freedom disappears at lower temperatures for systems with greater than the fixed . As we have mentioned earlier, when is of the order of the , the system behaves exactly like a D harmonic oscillator or equivalently with degrees of freedom. So, as we have seen from the figures, there are fractal dimensions for intermediate temperatures, especially the fractal dimension of for a system of harmonic oscillators in the LQGUP because three degrees of freedom are frozen for this model.




3.2 Specific heat
The high temperature behavior of the internal energy in the LQGUP framework could be also more precisely understood from another important thermodynamical quantity named as specific heat, that is defined from the relation (27) by the standard definition as follows
(29) |
The temperature-dependent behavior of this quantity is plotted in figure (4). Expansion of the specific heat (29) for the low and high temperatures regimes gives
(30) |
which shows that the specific heat approaches at the high temperature limit (see also figure (4)). As it can be observed, the equation (29) is a lengthy equation that gives the precise value of under a limiting operation for extremely high temperatures or equivalently for very short wavelengths, and must be taken into account to obtain the accurate value. If, for example, the Planck wavelength is neglected from the fourth order onwards, an incorrect and highly inaccurate value will be result.
According to the classical equipartition theorem, each of the degrees of freedom of the system has a contribution of to its internal energy and to its specific heat too. As it is clear from figure (4), the specific heat is bounded as in this setup.
As we have stated previously, when is of the order of , the reduction of the degrees of freedom will not be temperature dependent, and the system behaves completely like a dimension of of harmonic oscillators, or equivalently, with degrees of freedom. So, as we have seen from the figures, there are fractal dimensions for intermediate temperatures, especially the fractal dimension of for a system of harmonic oscillators in the LQGUP model because three degrees of freedom are frozen for this model.


3.3 Entropy
From Eq. (24), the entropy of the system in LQGUP framework will be as follows
(31) |
According to Eq. , the entropy of each system is directly proportional to the number of accessible microstates. Therefore, it is reasonable to expect that the number of microstates in the GUP framework is lower than in the classical case. This can be observed in Figure (5), where the entropy of the system in the LQGUP model increases at a slower rate compared to its classical counterpart. This result was anticipated, as we have previously observed a reduction in the degrees of freedom in this model, from to .


3.4 Numerical Results
The partition function obtained in equation (16) can be solved exactly by numerical techniques using Mathematica. Subsequently, other related thermodynamic variables can be determined exactly and their plots can be generated. In this subsection, in addition to numerical analysis encoded in the plots, we present some numerical results in the table (1). It will be observed that the approximate analytical results from the previous section closely align with the numeric findings. In the diagrams (2,3,4,5), we have plotted the temperature-dependent rate of change for various thermodynamic parameters such as internal energy, number of degrees of freedom, specific heat and entropy. As can be seen from the diagrams and the table, the analytical values determined in the previous section are very close to the exact numerical result.
T | 20 | 40 | 60 | 80 | 100 |
(analytical) | 61.6745 | 101.695 | 135.125 | 166.8 | 197.786 |
(numerical) | 60.0699 | 99.6916 | 133.094 | 164.78 | 195.778 |
T | 20 | 40 | 60 | 80 | 100 |
(analytical) | 6.16745 | 5.08475 | 4.50417 | 4.17 | 3.95572 |
(numerical) | 6.00699 | 4.98458 | 4.43648 | 4.11949 | 3.91557 |
T | 20 | 40 | 60 | 80 | 100 |
(analytical) | 2.42667 | 1.75936 | 1.61284 | 1.56211 | 1.53909 |
(numerical) | 2.37137 | 1.75497 | 1.61298 | 1.56276 | 1.53973 |
T | 20 | 40 | 60 | 80 | 100 |
(analytical) | 5.45347 | 6.8666 | 7.54631 | 8.00228 | 8.34809 |
(numerical) | 5.28132 | 6.6787 | 7.3578 | 7.8139 | 8.15986 |
4 The Case of 2D Harmonic Oscillator
In this section we consider a 2D harmonic oscillator to see the reduction of the degrees of freedom in this system and whether the results are as expected or not. The 2D LQGUP Liouville volume, , is obtained from equation (8) as . Hence, for 2D harmonic oscillators in LQGUP phase space, the partition function of a single particle in the semi-classical regime can be derived as follows
(33) |
Expanding the deformed phase space volume yields
(34) |
So, the partition function (33) can be approximated as follows
(35) |
Upon integration of the above relation, we obtain
(36) |
Expressing it in terms of the thermal de Broglie wavelength , we get
(37) |
Expanding the partition function (37) at very high and low temperatures gives
(38) |
It is evident from Eq. (38) that two degrees of freedom will freeze () at the Planck length scale in the LQGUP framework, and the partition function behaves as if it were effectively a partition function of a D harmonic oscillator.
The internal energy of the system can be modified as follows
(41) |
Expanding the relation (41) for both of the high and low temperature regimes, we have.
(42) |
Finally, expanding the specific heat which is defined through the relation (42) is as follows
(43) |
These results for a system of 2-dimensional harmonic oscillator are in agreement with the general prescription of the reduction of the degrees of freedom as expected.


5 Summary and Conclusions
In this paper, we have investigated the thermostatistics of a harmonic oscillator system in the Linear Quadratic GUP phase space model. Our aim was to gain further insight into the topic of the reduction of the degrees of freedom in high energy regimes. The calculations have been based on the existence of an ultraviolet cutoff (or minimal measurable length and maximal measurable momentum). Referring to the equipartition theorem of energy, we have clearly shown that as the temperature increases towards the Planck temperature in the mentioned framework, the three degrees of freedom freeze out. Such freezing (reduction) of the degrees of freedom has been observed for the Snyder Phase Space framework in a previous work by our research group [27]. The most important observation here is the fact that, while in the Snyder phase space scheme, the reduction of the number of degrees of freedom at high temperatures is -folds, in the GUP framework, this reduction is -folds.
In this direction, we note the following important points:
First, we note that the reduction of the degrees of freedom is essentially continuous, which means the possibility of having a non-integer degrees of freedom (see also figures 3). This may signal a fractal space dimension in the presence of the phenomenological quantum gravity effect encoded in the presence of a minimal measurable length.
Second, we considered the classical equipartition theorem for our treatment. It is in fact necessary to consider a modified equipartition theorem and even a Planck scale thermodynamics that may be in dramatic digression with the standard thermodynamics. But these are the issues of the yet-unsolved quantum gravity problem.
In the absence of a well-formulated quantum gravity theory, such phenomenological-based studies may shed light on the issue of the final quantum gravity proposal. So, we emphasize that in general, the deformed spaces in different phenomenological approaches to quantum gravity proposal all address a reduction of the degrees of freedom in high temperature regimes from a thermodynamic point of view. All of these theories coincide with their classical counterparts at low temperatures, as required. The possibility of realizing the fractal degrees of freedom in these setups is also an important outcome of this study.
The behavior shown in the diagrams arises in the field of quantum gravity, where the effects of the generalized uncertainty principle (GUP) lead to changes in the structure of phase space. The study of the LQGUP deformation scenario reveals a remarkable aspect: a significant deviation from classical behavior in the number of microstates at high energies. This deviation, characterized by the blocking of degrees of freedom, underlines a complex interplay between the basic principles of quantum mechanics and gravity. The shift towards a low-dimensional phase space at high energies emphasizes the intricate correlation between the dynamics of microstates and the inherent quantum essence of spacetime. These insights provide valuable perspectives on the nature of spacetime in a realm where the influences of quantum gravity are significant and provide a deeper understanding of the fundamental structure of the universe.
Finally we have checked the case with 2-dimensional harmonic oscillator to see the results which are in agreement with the general prescription: a reduction of the degrees of freedom from to , as expected.
Acknowledgement
We appreciate the referee for his/her insightful and constructive comments.
References
- [1] D. Kleppner and R. Jackiw, Science 289 (2000), 893-898.
- [2] A. Zeilinger, A. Nature vol. 408 (2000), 639-41.
- [3] G. Amelino-Camelia, Nature 408 (2000), 661-664.
- [4] D. J. Gross and P. F. Mende, Nucl. Phys. B 303 (1988) 407.
- [5] D. Amati, M. Ciafaloni and G. Veneziano, Phys. Lett. B 216 (1989) 41.
- [6] M. Kato, Phys. Lett. B 245 (1990) 43.
- [7] L. Garay, Int. J. Mod. Phys. A 10 (1995) 145.
- [8] C. Rovelli and L. Smolin, Nucl. Phys. B 331 (1990) 80.
- [9] C. Rovelli and L. Smolin, Nucl. Phys. B 442 (1995)593.
- [10] A. Ashtekar and J. Lewandowski, Class. Quantum Grav. 14 (1997) A55.
- [11] A. Ashtekar and J. Lewandowski, Class. Quantum Grav. 21 (2004) R53.
- [12] H. Snyder, Phys. Rev. 71 (1947) 38.
- [13] S. Mignemi, Phys. Rev. D 84 (2011) 025021.
- [14] S. Mignemi, Class. Quantum Grav. 29 (2012) 215019.
- [15] G. Amelino-Camelia, Int. J. Mod. Phys. D 11 (2000) 35.
- [16] G. Amelino-Camelia, Nature 418 (2002) 34.
- [17] J. Magueijo and L. Smolin, Phys. Rev. Lett. 88 (2002) 190403.
- [18] J. Kowalski-Glikman, Lect. Notes Phys. 669 (2005) 131.
- [19] G. Amelino-Camelia, L. Freidel, J. Kowalski-Glikman and L. Smolin, Phys. Rev. D 84 (2011) 084010.
- [20] M. Maggiore, Phys. Lett. B 319 (1993) 83.
- [21] L. Garay, Int. J. Mod. Phys. A 10 (1995) 145.
- [22] A. Kempf, G. Mangano and R. B. Mann, Phys. Rev. D 52 (1995) 1108.
- [23] A. Kempf and G. Mangano, Phys. Rev. D 55 (1997) 7909.
- [24] K. Nozari and A. Etemadi, Phys. Rev. D 85 (2012) 104029.
- [25] F. Scardigli, [arXiv:2212.13776].
- [26] T. Fityo, Phys. Lett. A 372 (2008) 5872.
- [27] K. Nozari, V. Hosseinzadeh and M. A. Gorji, Phys. Lett. B 750 (2015) 218.
- [28] K. Nozari, H. Ramezani, Int. J. Mod. Phys. D 26 (2017) no.11, 1750126.
- [29] S. Das, S. Ghosh, J. Phys. A: Math. Theor. 42 (2009) 345401, [arXiv:0812.3512].
- [30] S. Das, D. Roychowdhury, Phys. Rev. D 81 (2010) 085039.
- [31] S. Carlip, Class. Quantum Grav. 34 (2017) 193001.
- [32] S. Carlip, Universe 5 (2019) 83.
- [33] L. N. Chang, D. Minic, N. Okamura and T. Takeuchi, Phys. Rev. D 65 (2002), 125028.
- [34] E. C. Vagenas, A. F. Ali, M. Hemeda and H. Alshal, Eur. Phys. J. C 79 (2019) no.5, 398.
- [35] S. Bhattacharyya, S. Sen and S. Gangopadhyay, Eur. Phys. J. C 84 (2024) 425.
- [36] J. P. R. Bernaldez, A. G. Abac and R. E. S. Otadoy, Annals Phys. 457 (2023) 169402.
- [37] A. Rahaman, Int .J. Mod. Phys. A 38 (2023) 2350066
- [38] E. C. Vagenas, A. Farag Ali, M. Hemeda and H. Alshal, Annals Phys. 432 (2021) 168574.
- [39] H. Barman, A. Rahaman and S. K. Jha, Mod. Phys. Lett. A 36 (2021) 2150130.
- [40] M. Hemeda, H. Alshal, A. Farag Ali and Elias C. Vagenas, Nucl. Phys. B 1000 (2024) 116456.
- [41] J. Magueijo and L. Smolin, Phys. Rev. D 67 (2003), 044017.
- [42] J. Magueijo and L. Smolin, Phys. Rev. D 71 (2005), 026010.
- [43] S. Das and E. C. Vagenas, Phys. Rev. Lett. 101 (2008), 221301.
- [44] A. F. Ali, S. Das and E. C. Vagenas, Phys. Lett. B 678 (2009), 497-499.
- [45] S. Das and E. C. Vagenas, Can. J. Phys. 87 (2009), 233-240.
- [46] S. Das, E. C. Vagenas and A. F. Ali, Phys. Lett. B 690 (2010), 407-412
- [47] S. Basilakos, S. Das and E. C. Vagenas, JCAP 09 (2010), 027.
- [48] S. Das and E. C. Vagenas, Phys. Rev. Lett. 101 (2008) 221301.
- [49] P. S. Addison, Fractals and Chaos: An Illustrated Course. Institute of Physics Publishing 19. ISBN 0-7503-0400-6, UK (1997).
- [50] H. Lauwerier, Fractals: Endlessly Repeated Geometrical Figures. p. 36. ISBN 0-691-02445-6, Princton University Press, 1991, USA.