Mathematical Problems in Engineering
Volume 2015 (2015), Article ID 721637, 8 pages
http://dx.doi.org/10.1155/2015/721637

Research Article

Closed-Form Solutions of the Thomas-Fermi in Heavy Atoms and the Langmuir-Blodgett in Current Flow ODEs in Mathematical Physics

Efstathios E. Theotokoglou, Theodoros I. Zarmpoutis, and Ioannis H. Stampouloglou

School of Applied Mathematical and Physical Sciences, Department of Mechanics, Laboratory of Testing and Materials, National Technical University of Athens, Zographou Campus, Theocaris Building, 5 Heroes of Polytechniou Avenue, 157-73 Athens, Greece

Received 25 March 2015; Revised 2 July 2015; Accepted 7 July 2015

Academic Editor: Hang Xu

Copyright © 2015 Efstathios E. Theotokoglou et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Two kinds of second-order nonlinear, ordinary differential equations (ODEs) appearing in mathematical physics are analyzed in this paper. The first one concerns the Thomas-Fermi (TF) equation, while the second concerns the Langmuir-Blodgett (LB) equation in current flow. According to a mathematical methodology recently developed, the exact analytic solutions of both TF and LB ODEs are proposed. Both of these are nonlinear of the second order and by a series of admissible functional transformations are reduced to Abel’s equations of the second kind of the normal form. The closed form solutions of the TF and LB equations in the phase and physical plane are given. Finally a new interesting result has been obtained related to the derivative of the TF function at the limit.

The paper is dedicated to the memory of the authors’ Professor in National Technical University of Athens Dr. D. E. Panayotounakos who contributed a lot to the study of Thomas-Fermi and Langmuir-Blodgett equations

1. Introduction

An equation of considerable interest, called the Langmuir-Blodgett (LB) equation [1, 2], appeared in connection with the theory of flow of a current from a hot cathode to a positively charged anode in a high vacuum. The cathode and anode are long coaxial cylinders. It was proved [2, page 409] that the LB equation has an analytic expansion in the neighborhood of a fixed value of the independent variable, which assumes arbitrarily given values of the dependent variable and its derivative provided the zero value. However, the above equation does not admit exact analytic solution in terms of known tabulated functions [3, 4]. A lot of applications have been proposed for Langmuir equation. Horváth et al. [5] used Martin-Synge algorithm or anti-Langmuir in liquid chromatography, Shang and Zheng [6] proposed the system of Zakharov equations which involves the interaction between Langmuir and ion-acoustic waves in plasma, Graham and Cairns [7] posed the constraints on the formation and structure of Langmuir eigenmodes in the solar wind, and so forth.

On the other hand, the Thomas-Fermi (TF) equation appears in the problem of determining the effective nuclear charge in heavy atoms [2, 8, 9]. Numerous applications have also been proposed for TF equation. Tsurumi and Wadati [10] applied the TF approximation to obtain the dynamics of magnetically trapped boson-fermion, Liao [11] proposed an analytic technique for TF equation, Abbasbandy and Bervillier [12] used analytic continuation of Taylor series to confront the TF equation, Ourabah and Tribeche [13] presented a TF model based on thermal nonextensive relativistic effects, and so forth.

Both LB and TF equations do not admit exact analytic solution in terms of known (tabulated) functions. In this paper using series of admissible functional transformations, at first both LB and TF nonlinear ODEs are reduced to equivalent Abel’s equations of the second kind of the normal form. According to a mathematical methodology recently developed [14–16], we extract closed-form solution of the above Abel equations in both phase and physical planes under initial data in accordance with the physical problems.

The mathematical technique introduced in our study is general and can be applied to a very large class of yet unsolvable ODEs in nonlinear mechanics and generally in mathematical physics. The method employed is exact in accordance with the boundary conditions of each problem under consideration and gives us the possibility not to consider the a priori approximate construction of the solutions in the form of power series, some coefficient of which must be estimated. In addition with our analysis the exact value of the first order derivative of Thomas Fermi equation at is derived in a different way from the methods already existing in the literature [17–21].

2. Preliminaries: Notations

A nonlinear ODE of considerable interest, called the Langmuir-Blodgett (LB) equation [1], is the following:

Here the notations , are used for the total derivatives.

This equation of constant coefficients appeared in connection with the theory of flow of a current from a hot cathode to a positive charged anode in a high vacuum. The cathode and anode are long coaxial cylinders and the independent variable is defined by the equation: , where is the radius of the anode enclosing a cathode of radius . The independent variable is given by .

Equation (1) is converted into a somewhat more tractable form by means of the transformation proposed by Davis [2] and thus becomes

Numerical techniques concerning the solution of (1) are included in [6].

The TF nonlinear ODE

appears in the problem of determining the effective nuclear charge in heavy atoms [8, 9]. The solution is defined for the boundary values

Interesting approximate and numerical integrations methods have been proposed by Sommerfeld [22], Bush and Caldwell [17], Kobayashi et al. [20], Feynman et al. [18], and Coulson and March [19].

The differential equation (3) belongs to equations of Emden-Fowler type [3, 4]. Taking into account the known exact parametric solutions of this class of nonlinear ODEs [4, pages 241–250], the above mentioned equation can not solve analytically in terms of known functions.

Many investigations have studied the TF equation obeying to the boundary conditions (4) (see [18–21]) and they carried out the outward numerical integration (, ) for a sequence of values of initial slope converging step by step to that curve which should approach the -axis asymptotically. Finally for a semianalytical solution methodology we must refer to Sommerfeld [22].

The Emden-Fowler nonlinear ODE of the normal form is

where and , are arbitrary parameters. We note the following regarding the reduction of (5) (see [4]).

For and the admissible functional transformations

reduce the nonlinear ODE (5) to the following Abel equation of the second kind of the normal form

In what follows, by a series of admissible functional transformations, we prove that both of (1) and (3) can be exactly reduced to Abel’s equations of the second kind of the normal form . According to the following reduction procedure the exact analytic solutions of the above Abel equations are constructed, avoiding the approximate expression of the solutions in power series.

3. The Reduction Procedure

3.1. The LB Equation

The LB equation (1) is a second order nonlinear ODE of constant coefficients. Thus, by the substitution

it can be reduced to the following Abel equation of the second kind:

Furthermore, the well-known transformation [4, page 50]

giving the expressions

transforms (9) into the following Abel equation of the second kind:

Finally, the substitution

reduces (12) to the following Abel equation of the second kind of the normal form:

or, setting

to the form

The solution of (16), and thus the solutions of (14), (12), and (9), constitutes the intermediate integral of the LB equation (1) in the phase plane. In other words, when obtaining the solution of the transformed equation (16) in the form , = first integration constant, the solution of (14) becomes , the solution of (12) is , and finally the solution of (8) is . Thus, the solution to the original Langmuir equation (1) in the physical plane can be obtained by the integration by parts of the following equation:

Note that the reduced Abel equation (16) does not admit an exact analytic solution in terms of known (tabulated) functions [4, pages 29–45].

3.2. The TF Equation

Equation (3) is a typical Emden-Fowler nonlinear equation with , , and . A thorough examination concerning approximate expansion solution of (3) is included in [20] as well as all the references cited there. By now, with the aid of the transformations (6) one reduces the original TF equation (3) to the following Abel equation of the second kind of the normal form:

where

In what follows, based on a mathematical construction recently developed in [14, 15] concerning the closed-form analytic solution of an Abel equation of the second kind of the normal form , we provide the closed-form solutions of the reduced Abel equation (18), that is to say, the construction of the intermediate integral of the original TF equation (3) in the phase plane, as well as the final solution in the physical plane in accordance with the given boundary conditions.

4. Closed-Form Solutions of the TF Equation in the Phase Plane-Final Solutions

We consider the reduced Abel equation (18), namely, the equation

obeying (as the first of the boundary conditions (4) and the well-known result [2, 8, 17–20]) the transformed (19) and (20) equivalent condition:

It was recently proved [14, 15] that an Abel equation of the second kind of the normal form admits an exact analytic solution in terms of known (tabulated) functions. Hence, we perform the solution of (21) which constitutes the solution of the TF equation in the phase plane as follows:

where

Here is integration constant, while the function is given implicitly in terms of the subsidiary function , which can be determined from (23b), and the known member of (21) as in the following three cases.

Case 1 ( ()). Consider

Case 2 (). Consider

Case 3 (). Consider

In all these formulae the quantities ,  , and are given by

In order to define the type of the function () which must be selected among the closed-form solutions ((23a), (23b), and (23c)) to (26), as well as the value of the subsidiary function at and the value of the integration constant , we must combine the above Abel solutions with the boundary data.

We symbolize by the substitution while by any function at . Thus functions and given in (23a) perform

On the other hand, the Abel nonlinear ODE (21) at becomes

Combination of the above equations results in the following quadratic to equation:

including the value of and the constant on integration .

Similarly, by (23b) one gets

that constitute a second equation including the value of and the constant of integration . Also in this equation by means of (23c) we have

Summarizing, from (19) and (20) one concludes that, since by the first boundary conditions for then for or , where is a fixed value. Thus, combining (23a), (23b), and (23c) together with the first of (28) and the Abel nonlinear ODE (21), one may write that, for , , . Therefore, by making use of the first boundary condition given in (4) (, ,) as well as of all the above observations, the solutions of (31) and (32) perform the values for and and they can be written in the following algebraic forms:

where

Thus, we are now able to estimate , , and given by (27) in terms of the known values and . This estimation furnishes the value of and therefore the type of the solution () near the origin among the formulae (24) to (26).

From the above evaluation, it is concluded that the closed form solutions of the Abel nonlinear ODE (21) given as

have been completely determined.

Since we obtained the above closed- form solutions of the intermediate Abel nonlinear ODE (21) the solutions to the original TF equation (3) can be determined as follows. By (19) and (21) we perform

while taking the total differential in (19) we extract

where denotes the known solution of the Abel nonlinear ODE (18). Using (34a), (34b), and (37) one writes

that is, the Bernoulli equation

The last step of the above analysis concerns the calculation of the second constant of integration . Thus from the first of the boundary conditions (4) and the parametric solution (40), one writes

and estimates the constant of integration . By the second of the boundary conditions (4) one writes , , and . Again through the parametric solution (40) one writes

and evaluates the value of . The evaluation of the second constant of integration can be obtained again through the boundary condition ; that is,

Equations (41) and (42a), (42b) complete the solutions of the problem under consideration that is the construction of a closed-form parametric solution of the TF equation (3).

It is worthwhile to remark that since the solution is given in a closed form the mathematical methodology being developed results in the ad hoc definition of the TF function in approximate power series with ( = a factor that must be determined) and with boundary conditions given by ; as . The value can be evaluated by way of the already constructed closed-form solution. This is a very interesting result in accordance with previous results presented by Bush and Caldwell [17], Feynman et al. [18], Coulson and March [19], Kobayashi et al. [20], and Kobayashi [21].

5. Closed-Form Solutions of the LB Equation in the Physical Plane

It was already proved (Section 3) that the LB nonlinear ODE (1) has been reduced to an Abel ODE of the second kind of the normal form (15). It was also proved that this Abel equation admits exact analytic solutions given, as in case of the TF nonlinear ODE, by (24) to (26). Thus, the solution of the LB equation (1) in the phase plane is given by the formulae (24) to (26), if instead of the parenthesis , the quantity is introduced. The solution of the original equation (1) in the physical plane is given through the combination of the above prescribed formulae in combination with (17). As in case of the TF equation, the whole problem includes two constants of integration and , which are to be determined by using the convenient boundary data of the problem under consideration. In order to define simultaneously the type of the function () which must be selected in accordance with the modified previously developed formulae (24) to (26) and thus define the modified function and the constant of integration , we must combine the above Abel solutions with the prescribed convenient boundary data.

Based on what was previously developed in Section 3 for the original (LB) equation (1), the following boundary conditions hold true:

where and are arbitrary given values provided .

Based on these boundary conditions, the substitution (8) furnishes , while transformation (10) furnishes too. Thus, the substitutions (13), (15) and the Abel nonlinear ODE (16) perform the following relations:

We are now able to evaluate all the before mentioned quantities. Indeed, (23a) and (23b) for , , furnish the following two equations:

where and are the values of and at , respectively. On the other hand, the formula for in (23a) at results in

while by the first of (23a) and (44) we have

or equivalently

The last relation (48) in combination together with (44) permits us to write

The above three equations (46) to (49) constitute a nonlinear (transcendental) system for , , and by means of which one calculates these unknowns in terms of the parameters and . Thus, the discriminant ( at ) given in (27) is also known, the fact that permits us to define the type of the function for (24) to (26). This completes the solution of the problem under consideration, which is the construction of the solution of the LB equation in the phase-plane (and thus in the physical plane) in the main interval . If this solution is not unique inside , then we follow step by step the methodology developed in [15].

6. Conclusions

By a series of admissible functional transformations we reduce the nonlinear TF and LB equations to Abel’s equations of the second kind of the normal form of (16) and (18), respectively. These equations do not admit closed-form solutions in terms of known (tabulated) functions. This unsolvability is due to the fact that only very special forms of this kind of (16) and (18) can be solved in parametric form [4]. Our goal is the development of the construction of the exact analytic solutions of the above equations based on a mathematical technique leading to the derivation of close form solutions for the Abel equation of the second kind of the normal form [3, 4]. The proposed methodology constitutes the intermediate integral of the TF and LB equations in the phase and physical planes. The reduction procedure in the paper and the constructed solutions are very general and can be applied to a large number of nonlinear ODEs in mathematical physics and nonlinear mechanics. Also the method employed is exact in accordance with the boundary conditions of each problem under consideration and guides us to the a priori approximate construction of the solutions in power series, some coefficients of which must be evaluated.

With the proposed procedure the exact value of first order derivative of TF equation is derived at . In particular, the solution for the TF equation is finally constructed by calculating 4 constants from the specific boundary conditions (4) (see (34a), (34b), (41), (42a), and (42b)). Equations (40) constitute the parametric solutions of the TF problem. With our method essentially a different way of calculating (Section 4) is proposed with reference to the methods proposed in [17–21], where the value of has been calculated having kept the solution from approximate series in which their coefficients had been calculated from the boundary conditions. In conclusion, in this paper an analytic method is developed that does not contradict previous results but it is a different way of calculating them.

Conflict of Interests

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

References

I. Langmuir and K. B. Blodgett, “Currents limited by space charge between coaxial cylinders,” Physical Review, vol. 22, no. 4, pp. 347–356, 1923.

H. T. Davis, Introduction to Nonlinear Differential and Integral Equations, Dover, New York, NY, USA, 1962.

E. Kamke, Differentialgleichungen I, Gewöhnliche Differentialgleichungen and II, Partielle Differentialgleichungen, Akademische Verlagsgesellschaft, Leipzig, Germany, 1962.

A. D. Polyanin and V. F. Zaitsev, Handbook of Exact Solutions for Ordinary Differential Equations, CRC Press, Boca Raton, Fla, USA, 1999.

K. Horváth, J. N. Fairchild, K. Kaczmarski, and G. Guiochon, “Martin-Synge algorithm for the solution of equilibrium-dispersive model of liquid chromatography,” Journal of Chromatography A, vol. 1217, no. 52, pp. 8127–8135, 2010.

Y. Shang and X. Zheng, “The first-integral method and abundant explicit exact solutions to the zakharov equations,” Journal of Applied Mathematics, vol. 2012, Article ID 818345, 16 pages, 2012.

D. B. Graham and I. H. Cairns, “Constraints on the formation and structure of langmuir eigenmodes in the solar wind,” Physical Review Letters, vol. 111, no. 12, Article ID 121101, 2013.

L. H. Thomas, “The calculation of atomic fields,” Mathematical Proceedings of the Cambridge Philosophical Society, vol. 23, no. 5, pp. 542–548, 1927.

E. Fermi, “Statistical method of investigating electrons in atoms,” Zeitschrift für Physik, vol. 48, pp. 73–79, 1927.

T. Tsurumi and M. Wadati, “Dynamics of magnetically trapped boson-fermion mixtures,” Journal of the Physical Society of Japan, vol. 69, no. 1, pp. 97–103, 2000.

S. Liao, “An explicit analytic solution to the Thomas-Fermi equation,” Applied Mathematics and Computation, vol. 144, no. 2-3, pp. 495–506, 2003.

S. Abbasbandy and C. Bervillier, “Analytic continuation of Taylor series and the boundary value problems of some nonlinear ordinary differential equations,” Applied Mathematics and Computation, vol. 218, no. 5, pp. 2178–2199, 2011.

K. Ourabah and M. Tribeche, “Relativistic formulation of the generalized nonextensive Thomas-Fermi model,” Physica A, vol. 393, pp. 470–474, 2014.

D. E. Panayotounakos, “Exact analytic solutions of unsolvable classes of first and second order nonlinear ODEs (part I: Abel's equations),” Applied Mathematics Letters, vol. 18, no. 2, pp. 155–162, 2005.

D. E. Panayotounakos, N. B. Sotiropoulos, A. B. Sotiropoulou, and N. D. Panayotounakou, “Exact analytic solutions of nonlinear boundary value problems in fluid mechanics (Blasius equations),” Journal of Mathematical Physics, vol. 46, no. 3, Article ID 033101, 2005.

D. E. Panayotounakos, E. E. Theotokoglou, and M. P. Markakis, “Exact analytic solutions for the unforded damped duffing nonlinear oscillator,” Comptes Rendus Mechanics Journal, vol. 334, pp. 311–316, 2006.

V. Bush and S. H. Caldwell, “Thomas-Fermi equation solution by the differential analyzer,” Physical Review, vol. 38, no. 10, pp. 1898–1901, 1931.

R. P. Feynman, N. Metropolis, and E. Teller, “Equations of state of elements based on the generalized fermi-thomas theory,” Physical Review, vol. 75, no. 10, pp. 1561–1573, 1949.

C. A. Coulson and N. H. March, “Momenta in atoms using the Thomas-Fermi method,” Proceedings of the Physical Society Section A, vol. 63, no. 4, pp. 367–374, 1950.

S. Kobayashi, T. Matsukuma, S. Nagai, and K. Umeda, “Accurate value of the initial slope of the ordinary TF function,” Journal of the Physical Society of Japan, vol. 10, no. 9, pp. 759–762, 1955.

S. Kobayashi, “Some coefficients of the series expansion of the TFD function,” Journal of the Physical Society of Japan, vol. 10, no. 9, pp. 824–825, 1955.

S. A. Sommerfeld, “Integrazione Asintotica dell'Equazione Differenziali di Thomas-Fermi,” Accademia dei Lincei, Atti-Rendiconte, vol. 4, pp. 85–110, 1932.