# Thomas–Fermi equation

Last updated

In mathematics, the Thomas–Fermi equation for the neutral atom is a second order non-linear ordinary differential equation, named after Llewellyn Thomas and Enrico Fermi, [1] [2] which can be derived by applying the Thomas–Fermi model to atoms. The equation reads

## Contents

${\displaystyle {\frac {d^{2}y}{dx^{2}}}={\frac {1}{\sqrt {x}}}y^{3/2}}$

subject to the boundary conditions

${\displaystyle y(0)=1\quad$ ;\quad y(+\infty )=0}

If ${\displaystyle y}$ approaches zero as ${\displaystyle x}$ becomes large, this equation models the charge distribution of a neutral atom as a function of radius ${\displaystyle x}$. Solutions where ${\displaystyle y}$ becomes zero at finite ${\displaystyle x}$ model positive ions. [3] For solutions where ${\displaystyle y}$ becomes large and positive as ${\displaystyle x}$ becomes large, it can be interpreted as a model of a compressed atom, where the charge is squeezed into a smaller space. In this case the atom ends at the value of ${\displaystyle x}$ for which ${\displaystyle dy/dx=y/x}$. [4] [5]

## Transformations

Introducing the transformation ${\displaystyle z=y/x}$ converts the equation to

${\displaystyle {\frac {1}{x^{2}}}{\frac {d}{dx}}\left(x^{2}{\frac {dz}{dx}}\right)-z^{3/2}=0}$

This equation is similar to Lane–Emden equation with polytropic index ${\displaystyle 3/2}$ except the sign difference. The original equation is invariant under the transformation ${\displaystyle x\rightarrow cx,\ y\rightarrow c^{-3}y}$. Hence, the equation can be made equidimensional by introducing ${\displaystyle y=x^{-3}u}$ into the equation, leading to

${\displaystyle x^{2}{\frac {d^{2}u}{dx^{2}}}-6x{\frac {du}{dx}}+12u=u^{3/2}}$

so that the substitution ${\displaystyle u=e^{t}}$ reduces the equation to

${\displaystyle {\frac {d^{2}u}{dt^{2}}}-7{\frac {du}{dt}}+12u=u^{3/2}.}$

If ${\displaystyle w(u)={\frac {du}{dt}}}$ then the above equation becomes

${\displaystyle w{\frac {dw}{du}}-7w+12u=u^{3/2}.}$

But this first order equation has no known explicit solution, hence, the approach turns to either numerical or approximate methods.

## Sommerfeld's approximation

The equation has a particular solution ${\displaystyle y_{p}(x)}$, which satisfies the boundary condition that ${\displaystyle y\rightarrow 0}$ as ${\displaystyle x\rightarrow \infty }$, but not the boundary condition y(0)=1. This particular solution is

${\displaystyle y_{p}(x)={\frac {144}{x^{3}}}.}$

Arnold Sommerfeld used this particular solution and provided an approximate solution which can satisfy the other boundary condition in 1932. [6] If the transformation ${\displaystyle x=1/t,\ w=yt}$ is introduced, the equation becomes

${\displaystyle t^{4}{\frac {d^{2}w}{dt^{2}}}=w^{3/2},\quad w(0)=0,\ w(\infty )\sim t.}$

The particular solution in the transformed variable is then ${\displaystyle w_{p}(t)=144t^{4}}$. So one assumes a solution of the form ${\displaystyle w=w_{p}(1+\alpha t^{\lambda })}$ and if this is substituted in the above equation and the coefficients of ${\displaystyle \alpha }$ are equated, one obtains the value for ${\displaystyle \lambda }$, which is given by the roots of the equation ${\displaystyle \lambda ^{2}+7\lambda -6=0}$. The two roots are ${\displaystyle \lambda _{1}=0.772,\ \lambda _{2}=-7.772}$. Since this solution already satisfies the second boundary condition, to satisfy the first boundary condition one writes

${\displaystyle W=w_{p}(1+\beta t^{\lambda })^{n}=[144t^{3}(1+\beta t^{\lambda })]t.}$

The first boundary condition will be satisfied if ${\displaystyle 144t^{3}(1+\beta t^{\lambda })^{n}=144t^{3}\beta ^{n}t^{\lambda n}(1+\beta ^{-1}t^{-\lambda })^{n}\sim 1}$ as ${\displaystyle t\rightarrow \infty }$. This condition is satisfied if ${\displaystyle \lambda n+3=0,\ 144\beta ^{n}=1}$ and since ${\displaystyle \lambda _{1}\lambda _{2}=-6}$, Sommerfeld found the approximation as ${\displaystyle \lambda =\lambda _{1},\ n=-3/\lambda _{1}=\lambda _{2}/2}$. Therefore, the approximate solution is

${\displaystyle y(x)=y_{p}(x)\{1+[y_{p}(x)]^{\lambda _{1}/3}\}^{\lambda _{2}/2}.}$

This solution predicts the correct solution accurately for large ${\displaystyle x}$, but still fails near the origin.

## Solution near origin

Enrico Fermi [7] provided the solution for ${\displaystyle x\ll 1}$ and later extended by Edward B. Baker. [8] Hence for ${\displaystyle x\ll 1}$,

{\displaystyle {\begin{aligned}y(x)={}&1-Bx+{\frac {1}{3}}x^{3}-{\frac {2B}{15}}x^{4}+\cdots {}\\[6pt]&\cdots +x^{3/2}\left[{\frac {4}{3}}-{\frac {2B}{5}}x+{\frac {3B^{2}}{70}}x^{2}+\left({\frac {2}{27}}+{\frac {B^{3}}{252}}\right)x^{3}+\cdots \right]\end{aligned}}}

where ${\displaystyle B\approx 1.588071}$. [9] [10]

It has been reported by Salvatore Esposito [11] that the Italian physicist Ettore Majorana found in 1928 a semi-analytical series solution to the Thomas–Fermi equation for the neutral atom, which however remained unpublished until 2001.

Using this approach it is possible to compute the constant B mentioned above to practically arbitrarily high accuracy; for example, its value to 100 digits is ${\displaystyle B=1.588071022611375312718684509423950109452746621674825616765677418166551961154309262332033970138428665}$.

## Related Research Articles

In physics, the Lorentz transformations are a six-parameter family of linear transformations from a coordinate frame in spacetime to another frame that moves at a constant velocity relative to the former. The respective inverse transformation is then parameterized by the negative of this velocity. The transformations are named after the Dutch physicist Hendrik Lorentz.

In mathematics, the Rayleigh quotient for a given complex Hermitian matrix M and nonzero vector x is defined as:

In mathematics, the classical orthogonal polynomials are the most widely used orthogonal polynomials: the Hermite polynomials, Laguerre polynomials, Jacobi polynomials.

In mathematics, separation of variables is any of several methods for solving ordinary and partial differential equations, in which algebra allows one to rewrite an equation so that each of two variables occurs on a different side of the equation.

In mathematics and its applications, classical Sturm–Liouville theory is the theory of real second-order linear ordinary differential equations of the form:

In mathematics, the matrix exponential is a matrix function on square matrices analogous to the ordinary exponential function. It is used to solve systems of linear differential equations. In the theory of Lie groups, the matrix exponential gives the exponential map between a matrix Lie algebra and the corresponding Lie group.

In mathematics, an Euler–Cauchy equation, or Cauchy–Euler equation, or simply Euler's equation is a linear homogeneous ordinary differential equation with variable coefficients. It is sometimes referred to as an equidimensional equation. Because of its particularly simple equidimensional structure, the differential equation can be solved explicitly.

In quantum mechanics, the results of the quantum particle in a box can be used to look at the equilibrium situation for a quantum ideal gas in a box which is a box containing a large number of molecules which do not interact with each other except for instantaneous thermalizing collisions. This simple model can be used to describe the classical ideal gas as well as the various quantum ideal gases such as the ideal massive Fermi gas, the ideal massive Bose gas as well as black body radiation which may be treated as a massless Bose gas, in which thermalization is usually assumed to be facilitated by the interaction of the photons with an equilibrated mass.

The Pearson distribution is a family of continuous probability distributions. It was first published by Karl Pearson in 1895 and subsequently extended by him in 1901 and 1916 in a series of articles on biostatistics.

In numerical analysis, the Crank–Nicolson method is a finite difference method used for numerically solving the heat equation and similar partial differential equations. It is a second-order method in time. It is implicit in time, can be written as an implicit Runge–Kutta method, and it is numerically stable. The method was developed by John Crank and Phyllis Nicolson in the mid 20th century.

In general relativity, a geodesic generalizes the notion of a "straight line" to curved spacetime. Importantly, the world line of a particle free from all external, non-gravitational forces is a particular type of geodesic. In other words, a freely moving or falling particle always moves along a geodesic.

A differential equation can be homogeneous in either of two respects.

A multi-compartment model is a type of mathematical model used for describing the way materials or energies are transmitted among the compartments of a system. Sometimes, the physical system that we try to model in equations is too complex, so it is much easier to discretize the problem and reduce the number of parameters. Each compartment is assumed to be a homogeneous entity within which the entities being modeled are equivalent. A multi-compartment model is classified as a lumped parameters model. Similar to more general mathematical models, multi-compartment models can treat variables as continuous, such as a differential equation, or as discrete, such as a Markov chain. Depending on the system being modeled, they can be treated as stochastic or deterministic.

In mathematics, the spectral theory of ordinary differential equations is the part of spectral theory concerned with the determination of the spectrum and eigenfunction expansion associated with a linear ordinary differential equation. In his dissertation Hermann Weyl generalized the classical Sturm–Liouville theory on a finite closed interval to second order differential operators with singularities at the endpoints of the interval, possibly semi-infinite or infinite. Unlike the classical case, the spectrum may no longer consist of just a countable set of eigenvalues, but may also contain a continuous part. In this case the eigenfunction expansion involves an integral over the continuous part with respect to a spectral measure, given by the Titchmarsh–Kodaira formula. The theory was put in its final simplified form for singular differential equations of even degree by Kodaira and others, using von Neumann's spectral theorem. It has had important applications in quantum mechanics, operator theory and harmonic analysis on semisimple Lie groups.

A differential equation is a mathematical equation for an unknown function of one or several variables that relates the values of the function itself and of its derivatives of various orders. A matrix differential equation contains more than one function stacked into vector form with a matrix relating the functions to their derivatives.

The vibration of plates is a special case of the more general problem of mechanical vibrations. The equations governing the motion of plates are simpler than those for general three-dimensional objects because one of the dimensions of a plate is much smaller than the other two. This suggests that a two-dimensional plate theory will give an excellent approximation to the actual three-dimensional motion of a plate-like object, and indeed that is found to be true.

Euler's critical load is the compressive load at which a slender column will suddenly bend or buckle. It is given by the formula:

In mathematics, a linear recurrence with constant coefficients sets equal to 0 a polynomial that is linear in the various iterates of a variable—that is, in the values of the elements of a sequence. The polynomial's linearity means that each of its terms has degree 0 or 1. A linear recurrence denotes the evolution of some variable over time, with the current time period or discrete moment in time denoted as t, one period earlier denoted as t − 1, one period later as t + 1, etc.

The Fokas method, or unified transform, is an algorithmic procedure for analysing boundary value problems for linear partial differential equations and for an important class of nonlinear PDEs belonging to the so-called integrable systems. It is named after Greek mathematician Athanassios S. Fokas.

The Kaniadakis Logistic distribution is a generalized version of the Logistic distribution associated with the Kaniadakis statistics. It is one example of a Kaniadakis distribution. The κ-Logistic probability distribution describes the population kinetics behavior of bosonic or fermionic character.

## References

1. Davis, Harold Thayer. Introduction to nonlinear differential and integral equations. Courier Corporation, 1962.
2. Bender, Carl M., and Steven A. Orszag. Advanced mathematical methods for scientists and engineers I: Asymptotic methods and perturbation theory. Springer Science & Business Media, 2013.
3. pp. 9-12, N. H. March (1983). "1. Origins – The Thomas–Fermi Theory". In S. Lundqvist and N. H. March. Theory of The Inhomogeneous Electron Gas. Plenum Press. ISBN   978-0-306-41207-3.
4. March 1983, p. 10, Figure 1.
5. p. 1562,Feynman, R. P.; Metropolis, N.; Teller, E. (1949-05-15). "Equations of State of Elements Based on the Generalized Fermi-Thomas Theory" (PDF). Physical Review. American Physical Society (APS). 75 (10): 1561–1573. Bibcode:1949PhRv...75.1561F. doi:10.1103/physrev.75.1561. ISSN   0031-899X.
6. Sommerfeld, A. "Integrazione asintotica dell’equazione differenziale di Thomas–Fermi." Rend. R. Accademia dei Lincei 15 (1932): 293.
7. Fermi, E. (1928). "Eine statistische Methode zur Bestimmung einiger Eigenschaften des Atoms und ihre Anwendung auf die Theorie des periodischen Systems der Elemente". Zeitschrift für Physik (in German). Springer Science and Business Media LLC. 48 (1–2): 73–79. Bibcode:1928ZPhy...48...73F. doi:10.1007/bf01351576. ISSN   1434-6001.
8. Baker, Edward B. (1930-08-15). "The Application of the Fermi-Thomas Statistical Model to the Calculation of Potential Distribution in Positive Ions". Physical Review. American Physical Society (APS). 36 (4): 630–647. Bibcode:1930PhRv...36..630B. doi:10.1103/physrev.36.630. ISSN   0031-899X.
9. Comment on: “Series solution to the Thomas–Fermi equation” [Phys. Lett. A 365 (2007) 111], Francisco M.Fernández, Physics Letters A372, 28 July 2008, 5258-5260, doi : 10.1016/j.physleta.2008.05.071.
10. The analytical solution of the Thomas-Fermi equation for a neutral atom, G I Plindov and S K Pogrebnya, Journal of Physics B: Atomic and Molecular Physics20 (1987), L547, doi : 10.1088/0022-3700/20/17/001.
11. Esposito, Salvatore (2002). "Majorana solution of the Thomas-Fermi equation". American Journal of Physics. 70 (8): 852–856. arXiv:. Bibcode:2002AmJPh..70..852E. doi:10.1119/1.1484144.