# Thomas-Fermi model in the presence of natural cutoffs.

1. Introduction

According to equivalence principal in general relativity, gravitational field is coupled to everything. It has been characterized that gravity in very small length scales causes serious change in the structure of space-time. It causes minimal uncertainty in positions of atomic and subatomic particles [1-15]. In fact, there is an absolutely small uncertainty in position measurement of any quantum mechanical system and this feature leads nontrivially to the existence of a minimal measurable length on the order of Planck length. Existence of this natural cutoff requires deformation of the standard Heisenberg uncertainty principle to the so-called generalized uncertainty principle (GUP) (see, for instance, [13,14,16-21]).

n one dimension for positions and momentum operators, the deformed Heisenberg algebra can be represented as

[X,P] = i[??](1 + [beta][P.sup.2]). (1)

In general, for two symmetric operators A and B, we have

[DELTA]A[DELTA]B [greater than or equal to] [absolute value of <[A, B]>]. (2)

So the generalized uncertainty principle can be deduced as

[DELTA]X[DELTA]P [greater than or equal to] [[??]/2][1 + [beta][([DELTA]P).sup.2]]. (3)

While in ordinary quantum mechanics [DELTA]X can be made arbitrarily small by letting [DELTA]P grow correspondingly, this is no longer the case if (3) holds. If, for decreasing [DELTA]X, [DELTA]P increases, the new term [beta][([DELTA]P).sup.2] on the right hand side of (3) will eventually grow faster than the left hand side. Hence [DELTA]X can no longer be made arbitrarily small [16-20]. To obtain this minimal uncertainty, we saturate inequality in (3) and solve the resulting equation for [DELTA]P:

[DELTA]P = [DELTA]X [+ or -] [square root of [([DELTA]X).sup.2] - [[??].sup.2][beta])]/[??][beta]. (4)

The reality of solutions requires positivity of the term in square root, leading to

[([DELTA]X).sub.0] = [??][square root of [beta]]. (5)

This smallest uncertainty in position measurement leads nontrivially to the existence of a minimal measurable length. On the other hand, in the context of the doubly special relativity (DSR) theories (for review see [22-27]), one can show that a test particle's momentum cannot be arbitrarily imprecise. In fact there is an upper bound for momentum fluctuations [28-31]. As a nontrivial assumption, this may lead to a maximal measurable momentum for a test particle (see [20, 32-34]). In this framework, the GUP that predicts both a minimal observable length and a maximal momentum can be written (with [??] = 1) as follows 32, 33]:

[DELTA]X[DELTA]P

[greater than or equal to] [1/2][1 + ([alpha]/[square root of <[P.sup.2]>] + 4[[alpha].sup.2])[([DELTA]P).sup.2] + 4[[alpha].sup.2][&lt;P&gt;.sup.2] - 2[alpha][square root of <[P.sup.2]>]]. (6)

Since [([DELTA]P).sup.2] = <[P.sup.2]> - [&lt;P&gt;.sup.2], by setting &lt;P&gt; = 0 for simplicity, we find

[DELTA]X[DELTA]P [greater than or equal to] [1/2] [1 - [alpha] ([DELTA]P) + 4[[alpha].sup.2] [([DELTA]P).sup.2]]. (7)

This GUP contains both a minimal length and a maximal momentum. To see how a maximal momentum arises in this setup (see [20], for details), we note that with GUP (7) the absolute minimal measurable length is given by [DELTA][X.sub.min] (&lt;P&gt; = 0) [equivalent to] [DELTA][X.sub.0] = 3[alpha]/2. Due to duality of position and momentum operators, it is reasonable to assume [DELTA][X.sub.min] [varies] [DELTA][P.sub.max]. By saturating the inequality in relation (7), we find

2([DELTA]X[DELTA]P) = (1 - [alpha]([DELTA]P) + 4[[alpha].sup.2][([DELTA]P).sup.2]). (8)

This results in

[([DELTA]P).sup.2] - [(2[DELTA]X + [alpha])/4[[alpha].sup.2]][DELTA]P + [1/4[[alpha].sup.2]] = 0. (9)

So we obtain

[([DELTA][P.sub.max]).sup.2] - [(2[DELTA][X.sub.min] + [alpha])/4[[alpha].sup.2]] [DELTA][P.sub.max] + [1/4[[alpha].sup.2]] = 0. (10)

Now by using the value of [DELTA][X.sub.min], we find

[([DELTA][P.sub.max]).sup.2] - [1/[alpha]] [DELTA][P.sub.max] + [1/4[[alpha].sup.2]] = 0. (11)

The solution of this equation is

[DELTA][P.sub.max] = 1/2[alpha]. (12)

So, there is an upper bound on particle's momentum uncertainty. As a nontrivial assumption, we assume that this maximal uncertainty in particle's momentum is indeed the maximal measurable momentum. This is of the order of Planck momentum.

The presence of these strong natural constraints in measurement of position and momentum of particles has been ignored in studies of many-body interactions of condensed matter physics. Based on these preliminaries as our main motivation, we are going to study Thomas-Fermi theory in the presence of these natural cutoffs.

Thomas-Fermi theory is a basic theory in many-body interactions. For a system that consists of N particles each with mass m, the deformed algebra (1) is as follows:

[[X.sub.j], [P.sub.k]] = i[??][[delta].sub.jk] (1 + [beta][P.sup.2.sub.k]), [[X.sub.j], [X.sub.k]] = [[P.sub.j], [P.sub.k]] = 0, (13)

where j, k = 1, ..., N. If the particles of system are interacting via the pairwise potential V, the Hamiltonian of our N-body system is given by [24]

[H.sup.N] = [N.summation over (j=1)][[P.sup.2.sub.j]/2m] + [N.summation over (j<k=1)]V([X.sub.j] - [X.sub.k]). (14)

Generally, Hamiltonian of a many-body system with variety of interactions between particles is too difficult to handle. Nevertheless, someprogresscan be made if we neglect mutual interactions of particles. In this first step analysis the total Hamiltonian of the many-body system transforms to a set of independent single-particle Hamiltonians and these are usually easy to solve. A solid consists of array of atoms in close proximity. This array may be periodic, as in a crystal lattice, but this is not essential. Each atom consists of a positive nucleus surrounded by neutralizing set of electrons and preforms a many-body system. The Hamiltonian of a solid can be written as follows [35]:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (15)

where the nuclear positions are given by ([[??].sub.i], [[??].sub.j]) and the electron positions are given by ([[??].sub.k], [[??].sub.l]); Z is nuclear charge and (M, m) are the masses of the nucleus and electron, respectively. The primes on the summations mean that we must exclude the terms i = j, ... to prevent self-interactions. Within the Born-Oppenheimer approximation, we have (see [35], for details)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (16)

The second term of this relation shows the interaction between electrons that is the only obstacle to use single-particle wave function. Within Hartree's approximation we replace the interaction term

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (17)

by an approximate form

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (18)

Now we focus on correlation effects and the Thomas-Fermi model in a solid in equilibrium. Consider a system that includes electronic gas in equilibrium. In this system the local electron density is determined by the Fermi level which is a property of the solid as a whole. Suppose we apply a local electrostatic force that disturbs the system. In this situation, there is a tendency for electrons to migrate until the density is again consistent with the Fermi level. Suppose that [U.sub.0](r) is the potential of the applied electrostatic force and V(r) is the electrostatic potential due to movement of charges. Then the total potential of the system U(r) can be obtained from the following relation:

U(r) = [U.sub.0] (r) + V(r). (19)

The electrostatic potential satisfies Poisson's equation with the change in the charge density [rho](r):

[[nabla].sup.2]V (r) = 4[pi][rho] (r), (20)

where [rho](r) is given by

[rho](r) = eU (r)N([E.sub.f]) (21)

and N([E.sub.f]) is the density of states at the Fermi energy. Thus we have

[[nabla].sup.2] (U(r) - [U.sub.0](r)) = [[lambda].sup.2]U (r), (22)

where

[[lambda].sup.2] [equivalent to] 4[pi]eN ([E.sub.f]). (23)

Now (22) can be rewritten as

([[nabla].sup.2] - [[lambda].sup.2]) U (r) = [[nabla].sup.2][U.sub.0] (r). (24)

The corresponding homogeneous Thomas-Fermi equation is

([[nabla].sup.2] - [[lambda].sup.2])U(r) = 0. (25)

By using this equation we calculate the total potential U(r) and then we are able to obtain a relation for interaction potential. In what follows we are going to investigate Thomas-Fermi model firstly in the presence of a minimal measurable length encoded in a generalized uncertainty principle. We obtain total potential, interaction potential, and the correlation energy that is the rate of required energy for putting electron in the correlation hole and at the end we will find a relation for generalized dielectric function. In each step we compare our results with the standard case to see the effects of quantum gravity corrections. We note that the GUP as a manifestation of quantum gravity effects reflects the universality of quantum gravity corrections [36]. So, all quantum mechanical systems are affected by the quantum gravity effects and Thomas-Fermi model for many-body quantum interactions is not an exception in this regard. On the other hand, quantum gravity effects may have considerable effects on many-body quantum systems and one cannot neglect these effects in principle. It is important to note also that incorporation of the GUP effects in many-body problems may open new window to test quantum gravity in the lab.

2. Thomas-Fermi Equation in the Presence of Minimal Length

In one-dimensional case the deformed Heisenberg algebra reads

[X, P] = i[??](1 + [beta][P.sup.2]). (26)

One of the possible representations of this algebra can be written as

P = p(1 + [[beta]/3][p.sup.2]), X = x, (27)

where x and p are the conventional position and momentum operators. Note that while the existence of the natural cutoffs such as the minimal length is a common feature of all quantum gravity candidates, the form of the GUP containing this minimal length is model-dependent. In principle, these representations are not unique. For instance, we have supposed &lt;p&gt; = 0 in these representations and we have taken into account only terms of the first order in [beta] in a perturbative expansion. Since [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], so the deformed gradient operator ([??]') will be as follows:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (28)

Up to the first order in [beta] we have

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (29)

Now in deformed space (25) is replaced with

([[nabla]'.sup.2] - [[lambda].sup.2]) U' (r) = 0, (30)

where U'(r) is the total potential in deformed space. Therefore we find the following differential equation to first order in [beta]:

(-[2[beta][[??].sup.2]/3][[nabla].sup.4] + [[nabla].sup.2] - [[lambda].sup.2]) U'(r) = 0. (31)

We suppose the system is spherically symmetric so that there is no angular dependence and therefore

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (32)

Now (31) takes the following form:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (33)

This differential equation has the following solutions:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (34)

For a given [beta] with [beta] < 9/24[[??].sup.2][[lambda].sup.2], the terms with coefficients [C.sub.2] and [C.sub.4] are not bounded in the limit of r [right arrow] [infinity]. So, to have correct limiting case we set [C.sub.2] = [C.sub.4] = 0. On the other hand, only the first term gives the physically appropriate limiting result (note also that the term containing [C.sub.3] has no suitable limit; that is, in the limit of [beta] [right arrow] 0, it cannot recover the standard result in the absence of GUP). So we consider just this term and we set [C.sub.1] = e, where e is the electron charge. We define

[lambda]' [equivalent to] 1/2 [square root of (3 - [square root of (9 - 24[[??].sup.2][beta][[lambda].sup.2])]/[[??].sup.2][beta])]. (35)

So, we find

U'(r) = e[e.sup.-[lambda]'r]/r. (36)

This relation gives the effective interaction potential between electrons in Thomas-Fermi model in the presence of minimal measurable length. It is obvious that, in the limit of [beta] [right arrow] 0, [lambda]' should be transformed into the [lambda]. In fact, the limit of [beta] [right arrow] 0 can be calculated easily. The square root

[1/2] [square root of (3 - [square root of (9 - 24[[??].sup.2][beta][[lambda].sup.2])]/[[??].sup.2][beta])] (37)

can be expanded as

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (38)

The limit of this expression for [beta] [right arrow] 0 gives A as required. It could then be introduced, instead of the normal interaction potential equation (17), into the Schrodinger equation

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (39)

This induced potential causes the correlation energy which an electron needs in order to be located at the correlation hole. In solid state physics, to find this energy we use the following relation [35]:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (40)

where

[V.sub.IND] (r) [equivalent to] U'(r) - U(r [right arrow] 0). (41)

Using (36) we have

[V.sub.IND] (r) = -[e/2] (1 - [e.sup.-[lambda]'r]). (42)

Therefore, we find

[E.sub.CORR] = -[lambda]' [[e.sup.2]/2]. (43)

Note that the effect of GUP is hidden in the definition of [lambda]'. This energy appears as an additive constant to the energy of each electron.

3. The Generalized Dielectric Function in Deformed Space

In Thomas-Fermi model, the total potentialU(r) is obtained from a known external potential [U.sub.0](r). In classical electrostatics, we can introduce a dielectric constant that expresses the relationship between these two potentials:

U [approximately equal to] [U.sub.0]/[epsilon]. (44)

It would be correct also to describe this potential in terms of a dielectric function [epsilon](r):

U(r) = [U.sub.0](r)/[epsilon](r). (45)

Now we come back to (24) and write it in deformed space as

([[nabla]'.sup.2] - [[lambda].sup.2]) U(r) = [[nabla]'.sup.2][U.sub.0] (r), (46)

where

[[nabla]'.sup.2] = [[nabla].sup.2] - [2[beta][[??].sup.2]/3] [[nabla].sup.4]. (47)

So we find

[[nabla].sup.2]U(r) - [2[beta][[??].sup.2]/3] [[nabla].sup.4]U(r) - [[lambda].sup.2]U(r) = [[nabla].sup.2][U.sub.0](r) - [2[beta][[??].sup.2]/3][[nabla].sup.4][U.sub.0] (r). (48)

Taking Fourier transform and using a wave vector-dependent dielectric function [epsilon]([??]) we can write

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (49)

By substituting these relations in (48) we find

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (50)

So, we find

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (51)

Therefore,

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (52)

From relation [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], the generalized dielectric function in deformed space is obtained as

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (53)

To first order in [beta] and with [??] = 1, we can conclude that [??] in deformed space is defined as

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (54)

Therefore,

[epsilon]([??]) = 1 + [[[lambda].sup.2]/[q.sup.2] (1 + (2[beta]/3)[q.sup.2])]. (55)

It is important to note that corrections to the dielectric function in deformed space are wave number dependent. By the Fourier transform of U([??]), we can calculate the total potential U(r) to find

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (56)

Since

[U.sub.0](r) = e/[absolute value of [??]], (57)

the Fourier transform of this potential is

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (58)

thus we find

[U.sub.0] ([??]') = -4[pi]e/[q'.sup.2]. (59)

By substituting [??]' = [??](1 + ([beta]/3)[q.sup.2]) in this relation we find

[U.sub.0] ([??]) = -4[pi]e/[q.sup.2][(1 + ([beta]/3) 2)[q.sup.2]).sup.2]. (60)

To first order in [beta] this equation can be written as

[U.sub.0] ([??]) = -4[pi]e/[q.sup.2](1 + (2[beta]/3) 2)[q.sup.2]). (61)

From (55) and (61) we can obtain the total potential as follows:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (62)

thus,

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (63)

Figure 1 shows the real part and also the norm of U(r) for different values of [beta]. We note that, since GUP comes from quantum gravity effects via a perturbational scheme, it is natural to think that the GUP parameter is a very small parameter. The case with [beta] = 0 is the standard case and any deviation from this value incorporates the quantum gravitational effects. We have chosen the mentioned values of [beta] in figure just as examples and there is no reason to prevent selection of other nonvanishing values for [beta].

The correlation energy can be expressed in the form of

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (64)

where

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (65)

The Fourier transform of the [V.sub.IND]([??]) is as follows:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (66)

Thus,

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (67)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (68)

This integral can be calculated numerically to obtain [E.sub.CORR]. Figure 2 shows the real part and also the norm of [E.sub.CORR] for different values of [beta].

4. Thomas-Fermi Model in the Presence of Minimal Length and Maximal Momentum

Magueijo and Smolin have shown that in the context of the doubly special relativity a test particle's momentum cannot be arbitrarily imprecise and therefore there is an upper bound for momentum fluctuation [28-30] (see also [31]). Then it has been shown that this may lead nontrivially to a maximal measurable momentum for a test particle [32, 33]. In this framework, the GUP that predicts both a minimal observable length and a maximal momentum for a test particle can be written as follows [20, 32-34]:

[DELTA]X[DELTA]P [greater than or equal to] [[??]/2](1 - 2[beta] &lt;P&gt; + 4[[beta].sup.2] <[P.sup.2]>), (68)

or since [DELTA]P [equivalent to] P - &lt;P&gt;,

[DELTA]X[DELTA]P [greater than or equal to] [[??]/2] [1 - [beta]([DELTA]P) + 2[[beta].sup.2][([DELTA]P).sup.2]]. (69)

In this framework the following algebraic structure can be deduced (see 32, 33]):

[X, P] = i[??](1 - [beta]P + 2[[beta].sup.2][P.sup.2]), (70)

where [beta] is the GUP parameter in the presence of both a minimal length and a maximal momentum. Similar to the minimal length case, we can define [20, 34]

X = x,

P = p(1 - [beta]p + [[[beta].sup.2]/3] [p.sup.2]), (71)

P = P(1-PP+JP),

where, as before, x and p satisfy the canonical commutation relations via the Jacobi identity and X and P satisfy the generalized commutation relations in the presence of minimal length and maximal momentum [20, 32, 33]:

[X, P] = i[??](1 - [beta]p + [[[beta].sup.2]/3] [p.sup.2]). (72)

In this section we want to calculate the effects of both minimal length and maximal momentum on the total potential and generalized dielectric function in Thomas-Fermi model. In this case the deformed momentum is defined as

P [right arrow] P' = p(1 - [beta]p + [[[beta].sup.2]/3] [p.sup.2]). (73)

So the gradient operator will change as follows:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (74)

To first order in [beta], [[nabla]'.sup.2] can be written as

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (75)

Now, the Thomas-Fermi equation ([[nabla].sup.2] - [[lambda].sup.2])U(r) = 0 in this case transforms to the following equation:

(-[2[[beta].sup.2][[??].sup.2]/3][[nabla].sup.4] + 2i[[beta].sup.3][??][[nabla].sup.3] + [[nabla].sup.2] - [[lambda].sup.2])U(r) = 0. (76)

To solve this equation we assume [lambda] to be a complex quantity so that

[lambda] = [[lambda].sub.1] + i[[lambda].sub.2], [[lambda].sup.2] = [[lambda].sup.2.sub.1] - [[lambda].sup.2.sub.2] + 2i[[lambda].sup.*], (77)

where [[lambda].sup.*] [equivalent to] [[lambda].sub.1][[lambda].sub.2]. Then the total potential will be a complex function as

U(r) = [U.sub.R] + i[U.sub.I]. (78)

In spherical coordinate

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (79)

and using (32), we find

[-[2[[beta].sup.4][[??].sup.2]/3] [[nabla].sup.4] + [[nabla].sup.2] - ([[lambda].sup.2.sub.1] - [[lambda].sup.2.sub.2]) + i(-2[[lambda].sup.*] + 2[[beta].sup.3][??][[nabla].sup.3])]U(r) = 0; (80)

therefore,

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (81)

The real part of this complex differential equation can be solved to obtain

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (82)

where by definition

[lambda]' [equivalent to] [[lambda].sup.2.sub.1] - [[lambda].sup.2.sub.2]. (83)

Only the first term of the solution (82) satisfies the required boundary conditions and gives correct limiting results

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (84)

We suppose [C.sub.1] = e to find

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (85)

We note that the limit of [beta] [right arrow] 0 can be obtained easily much similar to the analysis that has been done just after (36).

In this case the effective interaction potential between electrons in the Thomas-Fermi model takes the following form:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (86)

and the correlation energy is

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (87)

so that

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (88)

To calculate the generalized dielectric function in the presence of both a minimal length and a maximal momentum, we proceed in the same manner as in the previous section. Starting with

([[nabla].sup.2] - [[lambda].sup.2]) U(r) = [[nabla].sup.2][U.sub.0](r) (89)

now the gradient operator changes as follows:

[[nabla].sup.2] [right arrow] [[nabla]'.sup.2] = [[nabla].sup.2] + 2i[[beta].sup.3][??][[nabla].sup.3] - [2[[beta].sup.4][[??].sup.2]/3][[nabla].sup.4] (90)

so we find

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (910

where

[[lambda].sup.2] = [[lambda].sup.2.sub.1] - [[lambda].sup.2.sub.2] + 2i[[lambda].sup.*]. (92)

Thus, we have

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (93)

Taking the Fourier transforms of U(r) and [U.sub.0](r), we find

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (94)

which gives

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (95)

So we find

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (96)

This gives finally

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (97)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (98)

From the relation [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], we can conclude that the generalized dielectric function in the presence of a minimal length and a maximal momentum is as follows:

[epsilon](q) = 1 + [[[lambda].sup.2]/[q.sup.2](1 + (2[beta]/3)[[??].sup.2][q.sup.2] - 2[beta][??]q)]. (99)

By setting [??] = 1 and considering only the terms that are first order in [beta], we find

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (100)

Thus,

[epsilon] ([??]) = 1 + [[[lambda].sup.2]/[q.sup.2](1 + (2[[beta].sup.2]/3)[q.sup.2] - 2[beta]q)]. (101)

To calculate the total potential U(r), we preform the Fourier transform of U([??]):

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (102)

As before we have

[U.sub.0] ([??]) = -4[pi]e/[q'.sup.2], (103)

where q' is given by (100). Therefore, taking only the terms that are first order in [beta], we find

[U.sub.0] ([??]) = -[4[pi]e/[q.sup.2](1 + (2[[beta].sup.2]/3)[q.sup.2] - 2[beta]q]. (104)

Now, the total potential can be written as

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (105)

So we have

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (106)

Figure 3 shows the real part and also norm of U(r) for different values of [beta].

To estimate correlation energy between electrons defined as

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (107)

we take the Fourier transform of [V.sub.IND]([??]) to find

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (108)

thus,

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (109)

Note that, due to the presence of the maximal momentum of the order of Planck momentum, this integral should be calculated for q ranging from -[q.sub.pl] to + [q.sub.pl] that [q.sub.pl] = [p.sub.pl]/[??]. Figure 4 shows the real part and also the norm of ECORR for different values of [beta].

Finally we note that it is possible to derive numerical bound on the value of the GUP parameter [beta] or the parameter [[beta].sub.0] defined as [beta] = [[beta].sub.0]/[M.sub.pl]c. As we have shown in [34], inspection of a bouncing particle in a gravitational field reveals the bound [[beta].sub.0] < [10.sup.29]. However, authors of [37] obtained the stronger bound of [[beta].sub.0] < [10.sup.10] by treating the Lamb shift in hydrogen atom.

5. Summary and Conclusions

In this paper the one-dimensional quantum N-body problem has been studied in condensed matter physics within the framework of a modified Heisenberg algebra admitting a minimal measurable length and maximal momentum. We have introduced within some approximation the effective Hamiltonian describing a solid and the interaction potential between electrons in this system. The formalism developed has been explicitly applied to the case of Thomas-Fermi theory that is a basic theory in many-body interactions. The system under study is an electronic gas in equilibrium that is disturbed with a local electrostatic force. At the first step we investigated this system in the presence of a minimal length and obtained total potential and interaction potential between electrons and then correlation energy and generalized dielectric function in this setup. Then we generalized our study to the case where there are both minimal length and maximal momentum as natural cutoffs. We have shown that, in contrast to the undeformed space, in this deformed space there are some new contributions in total potential and correlation energy that are determined by the deformation parameter. As an important result, we have shown that corrections to the dielectric function in deformed space are wave number dependent. This feature may provide a clue to see these tiny effects in the lab in future experiments. In conclusion we would like to stress that the formulation proposed here could be very useful in studying properties of particles in many-body systems in condensed matter physics and their various associated observables in the context of the deformed Heisenberg algebra.

http://dx.doi.org/10.1155/2014/418342

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgment

This work is financially supported by the Research Council of the Islamic Azad University, Sari Branch, Sari, Iran.

References

[1] G. Veneziano, "A stringy nature needs just two constants," Europhysics Letters, vol. 2, no. 3, article 199, 1986.

[2] D. Amati, M. Cialfaloni, and G. Veneziano, "Superstring collisions at planckian energies," Physics Letters B, vol. 197, no. 1-2, pp. 81-88, 1987

[3] D. Amati, M. Cialfaloni, and G. Veneziano, "Can spacetime be probed below the string size?" Physics Letters B, vol. 216, no. 1-2, pp. 41-47, 1989.

[4] D. J. Gross and P. F. Mende, "The high-energy behavior of string scattering amplitudes," Physics Letters B, vol. 197, no. 1-2, pp. 129-134, 1987.

[5] K. Konishi, G. Paffuti, and P. Provero, "Minimum physical length and the generalized uncertainty principle in string theory," Physics Letters B, vol. 234, no. 3, pp. 276-284, 1990.

[6] R. Guida, K. Konishi, and P. Provero, "On the short distance behavior of string theories," Modern Physics Letters A, vol. 6, no. 16, pp. 1487-1504, 1991.

[7] M. Kato, "Particle theories with minimum observable length and open string theory," Physics Letters B, vol. 245, no. 1, pp. 43-47, 1990.

[8] L. J. Garay, "Quantum gravity and minimal length," International Journal of Modern Physics A, vol. 10, no. 2, pp. 145-166, 1995.

[9] S. Capozziello, G. Lambiase, and G. Scarpetta, "Generalized uncertainty principle from quantum geometry," International Journal of Theoretical Physics, vol. 39, no. 1, pp. 15-22, 2000.

[10] M. Maggiore, "A generalized uncertainty principle in quantum gravity" Physics Letters B, vol. 304, no. 1-2, pp. 65-69, 1993.

[11] M. Maggiore, "The algebraic structure of the generalized uncertainty principle," Physics Letters B, vol. 319, no. 1-3, pp. 83-86, 1993.

[12] M. Maggiore, "Quantum groups, gravity, and the generalized uncertainty principle," Physical Review D, vol. 49, no. 10, pp. 5182-5187, 1994.

[13] S. Hossenfelder, "Minimal length scale scenarios for quantum gravity," Living Reviews in Relativity, vol. 16, 2013.

[14] S. Hossenfelder, "Can we measure structures to a precision better than the Planck length?" Classical and Quantum Gravity, vol. 29, no. 11, Article ID 115011, 2012.

[15] F. Scardigli, "Generalized uncertainty principle in quantum gravity from micro-black hole gedanken experiment," Physics Letters B, vol. 452, no. 1-2, pp. 39-44, 1999.

[16] A. Kempf, G. Mangano, and R. B. Mann, "Hilbert space representation of the minimal length uncertainty relation," Physical Review D, vol. 52, no. 2, pp. 1108-1118, 1995.

[17] H. Hinrichsen and A. Kempf, "Maximal localization in the presence of minimal uncertainties in positions and in momenta," Journal of Mathematical Physics, vol. 37, no. 5, pp. 2121-2137, 1996.

[18] A. Kempf, "On quantum field theory with nonzero minimal uncertainties in positions and momenta," Journal of Mathematical Physics, vol. 38, no. 3, pp. 1347-1372, 1997

[19] M. Bojowald and A. Kempf, "Generalized uncertainty principles and localization of a particle in discrete space," Physical Review D, vol. 86, no. 8, Article ID 085017, 15 pages, 2012.

[20] K. Nozari and A. Etemadi, "Minimal length, maximal momentum, and Hilbert space representation of quantum mechanics," Physical Review D, vol. 85, no. 10, Article ID 104029, 12 pages, 2012.

[21] R. J. Adler, "Six easy roads to the Planck scale," The American Journal of Physics, vol. 78, article 925, 2010.

[22] G. Amelino-Camelia, "Relativity in spacetimes with short-distance structure governed by an observer-independent (Planckian) length scale," International Journal of Modern Physics D, vol. 11, no. 1, pp. 35-59, 2002.

[23] G. Amelino-Camelia, "Special treatment," Nature, vol. 418, no. 6893, pp. 34-35, 2002.

[24] G. Amelino-Camelia, "Doubly-special relativity: first results and key open problems," International Journal of Modern Physics D, vol. 11, no. 10, pp. 1643-1669, 2002.

[25] J. Kowalski-Glikman, "Introduction to doubly special relativity," in Planck Scale Effects in Astrophysics and Cosmology, vol. 669 of Lecture Notes in Physics, pp. 131-159, Springer, Berlin, Germany, 2005.

[26] G. Amelino-Camelia, G. Mandanici, A. Procaccini, and J. Kowalski-Glikman, "Phenomenology of doubly special relativity," International Journal of Modern Physics A, vol. 20, no. 26, pp. 6007-6037, 2005.

[27] K. Imilkowska and J. Kowalski-Glikman, "Doubly special relativity as a limit of gravity," in Special Relativity, vol. 702 of Lecture Notes in Physics, pp. 279-298, Springer, Berlin, Germany, 2006.

[28] J. Magueijo and L. Smolin, "Lorentz invariance with an invariant energy scale," Physical Review Letters, vol. 88, no. 19, Article ID 190403, 4 pages, 2002.

[29] J. Magueijo and L. Smolin, "Generalized Lorentz invariance with an invariant energy scale," Physical Review D, vol. 67, no. 4, Article ID 044017, 12 pages, 2003.

[30] J. Magueijo and L. Smolin, "String theories with deformed energy-momentum relations, and a possible nontachyonic bosonic string," Physical Review D, vol. 71, no. 2, Article ID 026010, p. 6, 2005.

[31] J. L. Cortes and J. Gamboa, "Quantum uncertainty in doubly special relativity," Physical Review D, vol. 71, no. 6, Article ID 065015, p. 4, 2005.

[32] A. F. Ali, S. Das, and E. C. Vagenas, "Discreteness of space from the generalized uncertainty principle," Physics Letters B, vol. 678, no. 5, pp. 497-499, 2009.

[33] S. Das, E. C. Vagenas, and A. Farag Ali, "Discreteness of space from GUP II: relativistic wave equations," Physics Letters B, vol. 690, pp. 407-412, 2010.

[34] P. Pedram, K. Nozari, and S. H. Taheri, "The effects of minimal length and maximal momentum on the transition rate of ultra cold neutrons in gravitational field," Journal of High Energy Physics, vol. 2011, article 93, 2011.

[35] J. C. Inkson, Many-Body Theory of Solids, Plenum Press, New York, NY, USA, 1984.

[36] S. Das and E. C. Vagenas, "Universality of quantum gravity corrections," Physical Review Letters, vol. 101, Article ID 221301, 4 pages, 2008.

[37] A. Farag Ali, S. Das, and E. C. Vagenas, "Proposal for testing quantum gravity in the lab," Physical Review D, vol. 84, Article ID 044013, 2011.

Kourosh Nozari, Z. Haghani, and J. Vahedi

Department of Physics, Islamic Azad University, Sari Branch, Sari, Iran

Correspondence should be addressed to Kourosh Nozari; knozari@umz.ac.ir

Received 5 November 2013; Revised 26 January 2014; Accepted 29 January 2014; Published 28 April 2014