Some Consequences of a Simple Approach for Constructing a Theory of a Relativistic Fermi Gas ()

Roberto Lopez-Boada^{1}, Luis Grave de Peralta^{2,3*}

^{1}Science and Wellness Department, Broward College, Pembroke Pines, FL, USA.

^{2}Department of Physics and Astronomy, Texas Tech University, Lubbock, TX, USA.

^{3}Nano Tech Center, Texas Tech University, Lubbock, TX, USA.

**DOI: **10.4236/jmp.2021.1214113
PDF
HTML XML
222
Downloads
761
Views
Citations

It is shown how to build a simple but exact theory of a relativistic Fermi gas at 0 °K, which is based in a recently reported analytic formula for the energies of a relativistic spin-0 particle in a box. A white dwarf star is then simulated as a sphere filled with a relativistic Fermi gas. The Chandrasekhar mass limit is simply obtained using this model. We then discuss, using the proposed approach to relativistic quantum mechanics, how the interplay between the special theory of relativity, quantum mechanics, and gravity determines the stability of the matter.

Share and Cite:

Lopez-Boada, R. and Grave de Peralta, L. (2021) Some Consequences of a Simple Approach for Constructing a Theory of a Relativistic Fermi Gas. *Journal of Modern Physics*, **12**, 1966-1974. doi: 10.4236/jmp.2021.1214113.

1. Introduction

The theory of a non-relativistic Fermi gas is commonly discussed in the context of the theory of metals in solid state physics [1]. The theory of a non-relativistic Fermi gas is built from a well-known analytical expression, which gives the energies of a non-relativistic spin-0 particle in a three-dimensional infinite well (particle in a box) [2]. The Pauli exclusion principle is used, in the non-relativistic theory of a Fermi gas at 0 ˚K, for taking care of the fermion nature of the particles forming a Fermi gas [1] [2]. However, this relatively simple theory cannot be directly extended to the relativistic domain [3] [4]. This is because no such analytical expression was known, until recently, for the energies of a relativistic spin-0 particle in a box. In this work, we show, for the first time, how a simple but precise theory of a relativistic Fermi gas can be constructed. We use for this a recently reported analytical expression for the energies of a relativistic particle in a box [5] [6] [7] [8]. As it is shown in the Appendix, this analytical expression was obtained by solving a surprising Schrödinger-like but quasi-relativistic wave equation; therefore, using the same mathematical techniques required for obtaining the non-relativistic formula [2] [5] [6] [7] [8] [9].

In what follows, first, in Section 2, some results of such a simple theory of a relativistic Fermi gas, at a 0 ˚K temperature, are presented. While the consequences for solid state physics of this theory will be presented elsewhere, here, in Section 3, these results are used, for the first time, for obtaining in an alternative way, simple but precise, the Chandrasekhar mass limit of Fermi gas stars. Then, in Section 4 is presented a discussion about how the interplay between quantum mechanics, Newtonian gravity, and special relativity determines the stability of the matter. Finally, the conclusions of this work are given in Section 5.

2. The Approach

The theory of a non-relativistic Fermi gas, formed by *N* non-interacting fermions with spin-1/2, is based on the knowledge of the analytic expression of the energies corresponding to a 0-spin particle trapped in an infinite well (particle in a box), which can be easily calculated using the Schrödinger equation [1] [2]. The fermion nature of the particles is included by considering the Pauli exclusion principle, which is purely quantum in nature, and implies that only two fermions can simultaneously be at the same quantum state [1] [2].

Only recently, an analytic expression giving the energies of a spin-0 particle trapped in an infinite well was reported [5] [6] [7] [8]. This result was obtained using the GPPP (Grave de Peralta-Poirier-Poveda) method [5] [6] [7] [8] [9]. A fast review of the GPPP method is described for self-suffice reasons in the Appendix. It is worth noting that, first, there is an excellent agreement between the energies calculated using the analytic formula [5] [6] [7] [8], and the energies numerically calculated using the Dirac equation [10]. Second, it is well understood the relationship between the GPPP method and the Dirac equation [11]. Consequently, we can now construct a simple but exact theory of Fermi gases, which is valid from the non-relativistic to the ultrarelativistic regime. For doing this, we just need to repeat the same steps often follows for constructing the non-relativistic theory of a Fermi gas but using the new energy formula. We should start by substituting the non-relativistic formula for the energies of a quantum particle of mass *m* in a cubic box of size *L* [1] [2]:

${E}_{\left|n\right|}=\frac{{\pi}^{2}{\hslash}^{2}}{2m{L}^{2}}{\left|n\right|}^{2},\text{\hspace{0.17em}}\text{\hspace{0.17em}}n=\left({n}_{x},{n}_{y},{n}_{z}\right),\text{\hspace{0.17em}}\text{\hspace{0.17em}}{n}_{x,y,z}=1,2,\cdots $ (1)

by the corresponding formula for a relativistic particle [5] [6] [7] [8]:

${E}_{\left|n\right|}=\frac{{\pi}^{2}{\hslash}^{2}}{\left(\gamma +1\right)m{L}^{2}}{\left|n\right|}^{2}=\left(\gamma -1\right)m{c}^{2}.$ (2)

In Equations (1) and (2), *ħ* is the reduced Plank constant [1] [2], and *γ* is the especial relativity Lorentz factor [5] - [12]:

$\gamma =\sqrt{1+\frac{{p}_{\left|n\right|}^{2}}{{m}^{2}{c}^{2}}}=\sqrt{1+{\left|n\right|}^{2}{\pi}^{2}{\left(\frac{{\u019b}_{C}}{L}\right)}^{2}},\text{\hspace{0.17em}}\text{\hspace{0.17em}}{\u019b}_{C}=\frac{\hslash}{mc}.$ (3)

In Equation (3), *p* is the magnitude of the linear momentum of the particle in the box, and (*c*) is the speed of the light in the vacuum. In the non-relativistic limit,
$L\gg \hslash /mc$ ; therefore *γ* ≈ 1 and then Equation (2) coincides with Equation (1). In the ultrarelativistic limit,
$L\ll \hslash /mc$ ; therefore
$\gamma \gg \text{1}$, thus Equation (3) also gives the correct relation between *E* and *p* for an ultrarelativistic particle [3] [4]:

${E}_{\left|n\right|}={p}_{\left|n\right|}c,\text{\hspace{0.17em}}\text{\hspace{0.17em}}{p}_{\left|n\right|}=\frac{\left|n\right|}{L}h.$ (4)

In the thermodynamic limit,
$N\gg \text{1}$, each quantum state corresponds to a point in the “*n*-space” with energy given by Equation (2) [1]. In the ground state of the Fermi gas at 0 ˚K, all the energy levels up to the Fermi energy (*E _{F}*) level are occupied, and all the higher levels are empty [1]. The ground state is then represented by a three-dimensional isotropic and uniform Fermi sphere. Therefore, the number of states in the Fermi sphere and its radius are related by the following equation [1]:

$N=2\times \frac{1}{8}\times \frac{4}{3}\pi {\left|{n}_{F}\right|}^{3}\Rightarrow \left|{n}_{F}\right|={\left(\frac{3N}{\pi}\right)}^{1/3}.$ (5)

In Equation (5), the factor of two expresses the two spin states, and the factor of 1/8 expresses the fraction of the sphere that lies in the region where all *n _{x,y,z} *are positive. Substituting Equation (5) in Equation (2), and replacing

${E}_{F}=\frac{{\hslash}^{2}{\left(3{\pi}^{2}\right)}^{2/3}}{\left(\gamma +1\right)m}{\left(\frac{N}{V}\right)}^{2/3}=\left(\gamma -1\right)m{c}^{2},\text{\hspace{0.17em}}\text{\hspace{0.17em}}\gamma =\sqrt{1+{\left(3{\pi}^{2}\right)}^{2/3}{\u019b}_{C}^{2}{\left(\frac{N}{V}\right)}^{2/3}}.$ (6)

It is straightforward to show that Equation (6) gives the correct values of *E _{F}* in both the non-relativistic and the ultrarelativistic limits [1] [3] [4]. Moreover, Equation (6) is valid in all this range. The rest of the theory of a Fermi gas is constructed from the analytic formula for

${E}_{T}=\frac{{\hslash}^{2}}{m}{\left(\frac{3{\pi}^{2}}{V}\right)}^{2/3}{\displaystyle {\int}_{0}^{N}\frac{{\left({N}^{\prime}\right)}^{2/3}}{\gamma +1}\text{d}{N}^{\prime}}\approx N\theta {E}_{F}.$ (7)

The integral in Equation (7) can be easily calculated in the non-relativistic (*θ* = 3/5) and the ultrarelativistic (*θ* = 3/4) limits. Consequently, Equation (7) with *θ* slowly changing from 3/5 to 3/4 gives the exact value of *E _{T}* in the relativistic and ultrarelativistic limits [1] [3] [4], and it is a good approximation to

$P=-\frac{\partial {E}_{T}}{\partial V}\propto \{\begin{array}{l}{\left(\frac{N}{V}\right)}^{\frac{5}{3}},\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{non-relat}\\ {\left(\frac{N}{V}\right)}^{\frac{4}{3}},\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{ultrarelat}\end{array}$ (8)

We should emphasize here that the existence of the degenerate pressure in a Fermi gas is a purely quantum effect.

3. The Chandrasekhar Mass Limit

Using the Fermi gas as a model, it is possible to calculate the Chandrasekhar limit, *i.e.*, the maximum mass any white dwarf star may have (without significant thermally generated pressure) without collapsing into a black hole or a neutron star. The latter is a star mainly composed of neutrons, where the collapse is also avoided by neutron degeneracy pressure [13] [14].

A simple model for a Fermi gas star, formed by
$N\gg \text{1}$ spin-1/2 particles of mass *m*, is a spherical star of radius *r* formed by a Fermi gas of total mass *M _{S}* =

${E}_{s}\left(r\right)={E}_{T}\left(r\right)-{\alpha}_{G}\frac{G{M}_{S}^{2}}{r}.$ (9)

In Equation (9), *G* is the gravitational constant, and the parameter *α _{G}* is equal to 3/5 = 0.6 for a constant density sphere, but it takes slightly different values depending on the details of de internal structure of the star (

${E}_{s}\left(r\right)=\theta a\frac{{M}_{S}^{5/3}}{\left(\gamma +1\right){r}^{2}}-{\alpha}_{G}\frac{G{M}_{S}^{2}}{r},\text{\hspace{0.17em}}\text{\hspace{0.17em}}a={\left(\frac{3\pi}{4}\right)}^{2/3}\frac{{\hslash}^{2}}{{m}^{8/3}},\text{\hspace{0.17em}}\text{\hspace{0.17em}}\gamma =\sqrt{1+\frac{a}{{c}^{2}}\frac{{M}_{S}^{2/3}}{{r}^{2}}}.$ (10)

Or alternatively:

${E}_{s}\left(r\right)=\theta \left(\sqrt{1+\frac{a}{{c}^{2}}\frac{{M}_{S}^{2/3}}{{r}^{2}}}-1\right){M}_{S}{c}^{2}-{\alpha}_{G}\frac{G{M}_{S}^{2}}{r}.$ (11)

The radius of the Fermi gas star can be then estimated as the values of *r* for which *E _{S}*(

$r=\frac{\theta a}{{\alpha}_{G}G{M}_{S}^{1/3}}=\left[{\left(\frac{3\pi}{4}\right)}^{2/3}\frac{\theta}{{\alpha}_{G}}\right]{l}_{P}{\left(\frac{{m}_{P}}{m}\right)}^{3}{\left(\frac{m}{{M}_{S}}\right)}^{1/3}.$ (12)

For obtaining Equation (12), we used the following relations involving the Plank’s length (*l _{P}*) and mass (

$\frac{{\hslash}^{2}}{G{m}^{8/3}}=\frac{{\hslash}^{2}}{G{m}^{3}}{m}^{1/3},\text{\hspace{0.17em}}\text{\hspace{0.17em}}\frac{{\hslash}^{2}}{G{m}^{3}}={l}_{P}{\left(\frac{{m}_{P}}{m}\right)}^{3},\text{\hspace{0.17em}}\text{\hspace{0.17em}}{l}_{P}=\sqrt{\frac{\hslash G}{{c}^{3}}},\text{\hspace{0.17em}}\text{\hspace{0.17em}}{m}_{P}=\sqrt{\frac{\hslash c}{G}}.$ (13)

Therefore, a non-null radio, which decreases monotonically when *M _{S}* increases, exists for any star mass. This indicates the existence of an equilibrium between the centrifugal influence of the degeneracy pressure (a purely quantum mechanical effect) and the crushing gravitational force in the massive star. In contrast, when the particles forming the gas move at ultrarelativistic speeds, we can use either one of Equations (10) and (11) with:

$\gamma \approx \frac{\sqrt{a}}{c}\frac{{M}_{S}^{1/3}}{r}\gg 1.$ (14)

For obtaining:

${E}_{S}\left(r\right)\approx \left(\theta c\sqrt{a}{M}_{S}^{4/3}-{\alpha}_{G}G{M}_{S}^{2}\right)\frac{1}{r}.$ (15)

Consequently,
${E}_{S}\to -\infty $ when
$r\to 0$ ; *i.e.*, the star collapses, when:

${M}_{S}>{M}_{Ch}={\left(\frac{\theta}{{\alpha}_{G}}\right)}^{3/2}\sqrt{\frac{3\pi}{4}}{\left(\frac{{m}_{P}}{m}\right)}^{2}{m}_{P}.$ (16)

We obtained an excellent agreement with the Chandrasekhar mass limit (*M _{Ch}* ≈ 1.43 solar mases) [3] [11] [12], when taking

$r=\beta \left[{\left(\frac{3\pi}{4}\right)}^{2/3}\frac{\theta}{{\alpha}_{G}}\right]{l}_{P}{\left(\frac{{m}_{P}}{m}\right)}^{3}{\left(\frac{m}{{M}_{S}}\right)}^{1/3}\sqrt{1-{\left(\frac{{M}_{S}}{{M}_{Ch}}\right)}^{4/3}}.$ (17)

The dependence of *r* (in solar radii) on *M _{S}* (in solar mases) corresponding to Equation (17) is shown as a continuous red curve in Figure 1. We used here the same values of

4. Stability of the Matter

The curves shown in Figure 1 describe the stability of huge cosmological objects formed by numerous quantum particles. However, these curves streakily resemble other curves that were obtained while describing the stability of elemental

Figure 1. Radius-mass relations in a Fermi gas star model, relativistic (continuous red curve) vs non-relativistic (discontinuous blue curve).

particles and atoms [5]. It has been shown that a good estimate of the size of a Hydrogen-like atom, which is formed by an electron of mass *m _{e}* bounded to a nucleus of charge

$E\left(r\right)\approx \frac{{\hslash}^{2}}{\left(\gamma +1\right){m}_{e}{r}^{2}}-\frac{Z{e}^{2}}{4\pi {\u03f5}_{0}r},\text{\hspace{0.17em}}\text{\hspace{0.17em}}\gamma =\sqrt{1+{\left(\frac{{\u019b}_{C}}{r}\right)}^{2}}.$ (18)

Equation (18) is like Equation (10) but with the relativistic quantum mechanics term, *E _{T}*(

It is not then surprising that the dependence of the size of the atom (*r*) on *Z*, as shown in Figure 2(a), resembles the curves plotted in Figure 1. In this case, a non-relativistic description (continuous blue curve) predicts atoms with any value of *Z* should all be stable. However, when special relativity in included in the model (red points), and *Z* < 137, the electrical force between the nucleus and the electron tends to collapse the atom but pure quantum mechanical effects stabilize it. In contrast, atoms with *Z* > 137 are unstable (the electron collapses to the nucleus,
$r\to 0$ ) and thus these atoms should not exist [5]. Indeed, no atom with *Z *> 118 has ever been observed.

The following hypothesis, explaining why there are not elemental quantum particles with a mass larger than the Plank mass [5] [16], is much more closely related with the previous discussion about the fate of Fermi gas stars. If an elemental particle of mass *m* were able to interact gravitationally with itself, due the spread of its mass density through its wave function [5] [16] [17], the energy of the free particle could be estimated using the following equation [5] [16]:

$E\left(r\right)\approx \frac{{\hslash}^{2}}{\left(\gamma +1\right)m{r}^{2}}-\frac{G{m}^{2}}{r}.$ (19)

With *γ* given by Equation (18). The value of *r* that minimizes Equation (19) was found solving the equation d*E*/d*r* = 0 [5] [16]:

$r={a}_{D}\sqrt{1-{\left(\frac{m}{{m}_{P}}\right)}^{4}},\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.05em}}{a}_{D}={l}_{P}{\left(\frac{{m}_{P}}{m}\right)}^{3}.$ (20)

Figure 2. (a) Calculated radii of hydrogen-like atoms in
${\u019b}_{C}$ (reduced Compton wavelength) units as a function of *Z* for (continuous blue curve) *γ* ≈ 1 and (red points)
$\gamma \gg \text{1}$. (b) Calculated radii of the radius of a quantum particle, in Plank units, as a function of *m* in Plank units for (discontinuous blue curve) *γ* ≈ 1 and (continuous red curve)
$\gamma \gg \text{1}$.

In contrast, for the non-relativistic case; *i.e.*, using Equation (19) with *γ* = 1, it was found that *r* = *a _{D}* [5] [16]. As shown in Figure 2(b), a non-relativistic description (discontinuous blue curve) predicts elemental particles with any mass are possible. The dependence of

The stability of atoms, the so-called stability of the first kind [18], disappears when the combined effects of special relativity and electrostatic overcome the stability provides by quantum mechanics effects, thus producing the collapse of superheavy atoms. The stability of single quantum particles, a kind of zero-order stability, disappears when the combined effects of special relativity and gravity overcome the stability provided by quantum mechanics effects, thus producing the collapse of elemental particles. While elemental quantum particles with a mass larger than the Plank mass may not exist [5], massive cosmological bodies formed by an extremely large number of quantum particles do exist. The stability of the white dwarfs and neutron stars is an instance of the so-called stability of the second kind [18]. This stability also disappears when the combined effects of special relativity and gravity overcome the quantum effects, associated to the Pauli exclusion principle, that makes the stability possible.

5. Conclusion

We have shown that the analytical expression of the energies of a relativistic spin-0 particle trapped in a cubic box, which can be obtained using the GPPP approach (as shown in the Appendix), can be used as the foundation of a simple but exact theory of a Fermi gas, which is valid from the non-relativistic to the ultrarelativistic regimes. For instance, we were able to obtain the Chandrasekar mass limit using this approach and, due to the simplicity of our approach, we were able to illustrate the interplay between gravity, quantum mechanics, and special theory of relativity that determines the stability of matter.

Appendix

The Poirier-Grave de Peralta equation [5] [19]:

$i\hslash \frac{\partial}{\partial t}\Psi =\left[\frac{{\stackrel{^}{p}}^{2}}{\left(\stackrel{^}{\gamma}+1\right)m}+V\right]\Psi ,\text{\hspace{0.17em}}\text{\hspace{0.17em}}\stackrel{^}{\gamma}=\sqrt{1+\frac{{\stackrel{^}{p}}^{2}}{{m}^{2}{c}^{2}}},\text{\hspace{0.17em}}\text{\hspace{0.17em}}\stackrel{^}{p}=-i\hslash \nabla .$ (A1)

is fully Lorentz-covariant [5]. It can be shown that for a particle in the box, Equation (A1) exactly reduces to solving the following Schrödinger-like equation [5] [19]:

$i\hslash \frac{\partial}{\partial t}\Psi =-\frac{{\hslash}^{2}}{\left(\gamma +1\right)m}{\nabla}^{2}\Psi +V\Psi .$ (A2)

In Equation (A2), *γ* is not an operator but the parameter (GPPP approach) [5]:

$\gamma =\sqrt{1+\frac{2}{m{c}^{2}}{E}_{Sch}},\text{\hspace{0.17em}}\text{\hspace{0.17em}}{E}_{Sch}=\langle {\psi}_{Sch}|\frac{{\stackrel{^}{p}}^{2}}{2m}|{\psi}_{Sch}\rangle .$ (A3)

In Equation (A3), *E _{Sch}* and

${E}_{\left|n\right|}=\frac{{\pi}^{2}{\hslash}^{2}}{\left(\gamma +1\right)m{L}^{2}}{\left|n\right|}^{2},\text{\hspace{0.17em}}\text{\hspace{0.17em}}\gamma =\sqrt{1+\frac{2}{m{c}^{2}}{E}_{Sch}}.$ (A4)

Equations (2) and (3) follows directly from the substitution of *E _{Sch}* =

Conflicts of Interest

The authors declare no conflicts of interest regarding the publication of this paper.

[1] | Kittel, C. (2005) Introduction to Solid State Physics. 8th Edition, J. Wiley & Sons, New York. |

[2] |
Griffiths, D.J. (2018) Introduction to Quantum Mechanics. 3rd Edition, Prentice Hall, New York. https://doi.org/10.1017/9781316995433 |

[3] | Greiner, W., Neise, L. and Stöcker, H. (1997) Thermodynamic and Statistical Mechanics. Springer-Verlag, New York. |

[4] | Landau, L. and Lifshits, E.M. (1970) Statistical Physics. 2nd Edition, Pergamon Press, Oxford. |

[5] |
Grave de Peralta, L., Poveda, L. and Poirier, B. (2021) European Journal of Physics, 42, Article ID: 055404. https://doi.org/10.1088/1361-6404/ac0ecc |

[6] |
Grave de Peralta, L. (2020) The European Journal of Physics, 41, Article ID: 065404. https://doi.org/10.1088/1361-6404/aba7dc |

[7] |
Grave de Peralta, L. (2020) Results in Physics, 18, Article ID: 103318. https://doi.org/10.1016/j.rinp.2020.103318 |

[8] |
Grave de Peralta, L. (2020) Journal of Modern Physics, 11, 196-213. https://doi.org/10.4236/jmp.2020.112012 |

[9] |
Grave de Peralta, L. (2020) Scientific Reports, 10, Article No. 14925. https://doi.org/10.1038/s41598-020-71505-w |

[10] |
Ruiz-Columbié, A., Farooq, H. and Grave de Peralta, L. (2021) Journal of Modern Physics, 12, 1358-1374. https://doi.org/10.4236/jmp.2021.1210085 |

[11] |
Grave de Peralta, L. and Farooq, H. (2021) Journal of Modern Physics, 12, 1145-1159. https://doi.org/10.4236/jmp.2021.128068 |

[12] | Christodeulides, C. (2016) The Special Theory of Relativity. Springer, New York. |

[13] | Morison, I. (2008) Introduction to Astronomy and Cosmology. Wiley, London. |

[14] |
Kutner, M.L. (2003) Astronomy, a Physical Perspective. 2nd Edition, Cambridge University Press, Cambridge. https://doi.org/10.1017/CBO9780511802195 |

[15] | Walker, J., Resnick, R. and Halladay, D. (2014) Fundamentals of Physics. 10th Edition, J. Wiley & Sons, New York. |

[16] |
Poveda, L.A., Ruiz-Columbié, A. and Grave de Peralta, L. (2021) Relativistic Corrections to the Diósi-Penrose Model. http://arxiv.org/abs/2103.00994 |

[17] |
Diósi, L. (1984) Physics Letters A, 105, 199. https://doi.org/10.1016/0375-9601(84)90397-9 |

[18] |
Lieb, H.E. and Seiringer, R. (2010) The Stability of Matter in Quantum Mechanics. Cambridge University Press, New York. https://doi.org/10.1017/CBO9780511819681 |

[19] | Grave de Peralta, L., Poirier, B. and Poveda, L.A. (2021) APS March Meeting U71-00234. |

Journals Menu

Contact us

+1 323-425-8868 | |

customer@scirp.org | |

+86 18163351462(WhatsApp) | |

1655362766 | |

Paper Publishing WeChat |

Copyright © 2024 by authors and Scientific Research Publishing Inc.

This work and the related PDF file are licensed under a Creative Commons Attribution 4.0 International License.