Approximate rational solutions to the Thomas–Fermi equation based on dynamic consistency
Introduction
The major purpose of this paper is to construct a new rational approximation to the Thomas–Fermi equation with boundary conditions
This equation arises in the modeling of certain phenomena in the atomic structure of matter [1], [2].
Since it is expected that no exact analytical solution formula exists for the general solution in terms of a finite combination of elementary functions, a broad range of approximate solutions have been derived [3], [4], [5], [6], [7], [8], [9]. The mathematical procedures used to obtain these expressions include variational methods [3], [6], the use of iteration techniques [4], application of the Adomian decomposition method [5], homotopy analysis [7], and rational approximations [8]. The indicated references [3], [4], [5], [6], [7], [8], [9] are but a small subset of the hundreds of publications on this subject. Note that the existence of such a substantial research literature is due, in large part, to the fact that this thorny boundary-value problem has drawn out a significant number of ideas as to how approximations to its solutions should be done.
The purpose of this paper is to use a new methodology to construct approximate solutions to the Thomas–Fermi (TF) equation. We do this within the framework of the concept or principle of dynamic consistency [10] and the restrictions that it imposes on any mathematical representation of the forms used for rational approximations to a TF equation solution. We demonstrate that this is rather easy to do.
In outline, Section 2 provides a summary of several of the general, exact, and for the most part qualitative properties of the TF equation [11]. This section includes several new results of the authors. Section 3 introduces and defines the concept of dynamic consistency [10] and it is used in Section 4 to construct two elementary rational approximations to the TF solution. Finally, in Section 5, we provide a summary of our results and discuss possible extensions to this work.
Section snippets
Preliminaries: Exact results
It turns out that even in the absence of knowledge of an explicit, exact solution for the TF equation, many of the essential properties of such a solution can be derived. Below, we give a concise summary of these items and refer the reader to the paper of Hille [11] for proofs of some of these statements.
(a) An exact solution to the TF equation is Note that it contains no arbitrary parameters and thus is not a special case of the (unknown) general solution. We denote it as, and
Dynamic consistency
Consider two systems and . Let system have the property . If system also has property , then is said to be dynamic consistent to , with respect to property [10].
Dynamic consistency (DC) has served as one of the fundamental principles of the nonstandard, finite difference methodology for constructing improved discrete models for the numerical integration of differential equations [10]. It serves as an assessment of the “closeness” or “fidelity” of two systems based on their
Approximate solutions
We take the following two ansatzes to represent approximations to the solution of the TF equation where , and as a consequence of (33), .
Close inspection of and shows that they are DC with the five conditions listed in Section 3. Observe that is a rational function of , while is a rational function of .
Fig. 2 gives plots of and vs. . The two curves lie close to each other, but
Discussion
In summary, we have constructed two representations of approximate solutions to the TF equation. However, it must be stated that this equation is itself an approximation to an (unknown) equation describing atomic quantum phenomena [1], [2]. In addition, all numerical procedures have errors associated with their algorithms and the hardware on which they are implemented.
The methodology used to generate the ansatzes, given in Eqs. (31), (32), (33), is based on constructing for , the exact
References (15)
Solution to the Thomas–Fermi equation
Appl. Math. Lett.
(1998)Un methodo statistico per la determinazione di alcunepriorieta dell’atomé
Rend. Acad. Naz. Lincei
(1927)The calculation of atomic fields
Proc. Camb. Phil. Soc.
(1927)- et al.
Variational solutions to the Thomas–Fermi equation
Nuovo Cimento B Ser. 10
(1968) - et al.
The analytical solution of the Thomas–Fermi equation for a neutral atom
J. Phys. B: At. Mol. Phys.
(1987) - et al.
Variational approach to the Thomas–Fermi equation
Eur. J. Phys.
(2004) Analytic approximation for Thomas–Fermi equation
Acta Phys. Polon. A
(2008)
Cited by (5)
Dynamical analysis of a generalized hepatitis B epidemic model and its dynamically consistent discrete model
2023, Mathematics and Computers in SimulationCitation Excerpt :Hence, NSFD schemes are effective to simulate dynamics of dynamical differential models over long time periods. Because of this, nowadays NSFD schemes have been recognized as one of the effective approaches to solve differential equation equations arising in theory and practice [1,5,10,14,26,28,40,54,57–59,63,64]. In recent works [16–21,30,32], we have successfully developed the Mickens’s methodology to construct NSFD schemes for some mathematical models arising in biology and epidemiology.
A dynamically consistent nonstandard finite difference scheme for a generalized SEIR epidemic model
2023, Journal of Difference Equations and ApplicationsA novel second-order nonstandard finite difference method preserving dynamical properties of a general single-species model
2023, International Journal of Computer Mathematics