JMPJournal of Modern Physics2153-1196Scientific Research Publishing10.4236/jmp.2021.122008JMP-106572ArticlesPhysics&Mathematics Brownian Motion in an External Field Revisited AngeloPlastino1*MarioCarlos Rocca2DianaMonteoliva3AlbertoHernando4Departamento de Física, Universidad Nacional de La Plata, La Plata, ArgentinaKido Dynamics SA, Lausanne, SwitzerlandConsejo Nacional de Investigaciones Científicas y Tecnológicas, (IFLP-CCT-CONICET)-C. C. 727, La Plata, ArgentinaDepartamento de Mateática, Universidad Nacional de La Plata, La Plata, Argentina13012021120282901, December 202015, January 2021 18, January 2021© Copyright 2014 by authors and Scientific Research Publishing Inc. 2014This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/

In many interesting physical examples, the partition function is divergent, as first pointed out in 1924 by Fermi (for the hydrogen-atom case). Thus, the usual toolbox of statistical mechanics becomes unavailable, notwithstanding the well-known fact that the pertinent system may appear to be in a thermal steady state. We tackle and overcome these difficulties hereby appeal to firmly established but not too well-known mathematical recipes and obtain finite values for a typical divergent partition function, that of a Brownian particle in an external field. This allows not only for calculating thermodynamic observables of interest, but for also instantiating other kinds of statistical mechanics’ novelties.

Divergent Partition Functions Statistical Mechanics Fisher Information1. Introduction

In many interesting physical examples, the partition function is divergent [1] [2] [3] [4]. Thus, the usual toolbox of statistical mechanics becomes unavailable, notwithstanding the well-known fact that the pertinent system may appear to be in a thermal steady state (see, for instance [5] [6] [7] [8] [9] ) and references therein]. Our goal here is to deal with a specific divergent partition function, and obtain a finite value for it. This permits to compute new observables of interest and also to develop some hopefully new statistical mechanics’ insights.

2. The Central Issue2.1. Partition Function

We will consider here the partition function for Brownian motion in an external field, given by [4]

Z = ∫ − ∞ ∞     e β U 0 1 + x 2 d x , (2.1)

with β = 1 / ( k B T ) and k B = Boltzmann’s constant. Change now variables to y = 1 + x 2 . Taking advantage now of well-known features of Schwartz’ theory of distributions [10], we can recast the integral that defines Z in the fashion

Z = ∫ 1 ∞ ( y − 1 ) − 1 2 e β U 0 y d y ≡ lim ν → 1 ∫ 1 ∞     y ν − 1 ( y − 1 ) − 1 2 e β U 0 y d y , (2.2)

and remember that the limit of an integral equals the integral of the limit. We consult then the Tableof Ref. [11] and find that our current integral is a special case of the more general one

W = ∫ u ∞     x ν − 1 ( x − u ) μ − 1 e β x d x = B ( 1 − μ − ν , μ ) u μ + ν − 1 ϕ ( 1 − μ − ν ; 1 − ν , β u ) . (2.3)

Here B is the well-known beta function and ϕ the confluent hypergeometric function, that reads, appealing to the Gamma function Γ ,

B = Γ ( 1 − μ − ν ) Γ ( μ ) / Γ ( 1 − ν ) . (2.4)

Comparing integrals, we see at this stage that the right hand side of (2.2) will coincide with W in (2.3) by setting

μ = 1 / 2 ;             ν = 1 ;             u = 1 , (2.5)

so that these special values are to be inserted in

W = Γ ( 1 − μ − ν ) [ Γ ( μ ) / Γ ( 1 − ν ) ] u μ + ν − 1 ϕ ( 1 − μ − ν ; 1 − ν , β u ) . (2.6)

Note also that

Γ ( 1 / 2 ) = π ;             Γ ( − 1 / 2 ) = − 2 π . (2.7)

We have a Γ ( 0 ) in a denominator now. This induces us to appeal once again to [11] to employ the useful relation

lim γ → 0 ϕ ( α ; γ ; s ) = z α ϕ ( α + 1 ; 2 ; z ) , (2.8)

so that we can finally arrive at the result

Z = π β U 0 ϕ ( 1 2 ; 2 ; β U 0 ) , (2.9)

our desired finite form. We see that we arrive at Z via a straightforward path. The essential step here is that of consulting an appropriate table of integrals and performing adequate manipulations. Note that at very low temperatures quantum effects raise their head and our treatment becomes invalid. Below it will be shown that one also encounters problems or exceedingly high temperatures. We have found a finite partition function for our Brownian problem and proceed to calculate with it, below, important quantifiers of statistical mechanics.

2.2. Units for Our Graphs

We find it convenient to plot our thermal quantities versus y = k B T / U 0 in the range 0 ≤ y ≤ 1 . Given the smallness of k B , this encompasses an immense T-range, since k B is of the order of 10−23 in its appropriate units. In particular, we plot the logarithm of the partition function in Figure 1. We appreciate the fact that it converges to a definite value as T grows.

3. Other Thermal Quantities3.1. Mean Energy

One has

〈 U 〉 = − ∂ ln Z ∂ β , (3.1)

so that

〈 U 〉 = − 1 Z [ π U 0 ϕ ( 1 2 ; 2 ; β U 0 ) + π β U 0 2 4 ϕ ( 3 2 ; 3 ; β U 0 ) ] . (3.2)

Note that at very low temperatures quantum effects raise their head and our classical treatment becomes invalid.

3.2. Entropy S

We have

S = ∂ ( k B T ln Z ) ∂ T , (3.3)

so that

S = ln [ π β U 0 ϕ ( 1 2 ; 2 ; β U 0 ) ] − β Z [ π U 0 ϕ ( 1 2 ; 2 ; β U 0 ) + π β U 0 2 4 ϕ ( 3 2 ; 3 ; β U 0 ) ] , (3.4)

that is plotted in Figure 2. Note that at very low temperatures, quantum effects raise their head and our treatment becomes invalid. This is evident whenever S becomes negative at low T. A new effect is observed at very large T. Whenever T ≥ 10 22 , the treatment becomes invalid as well. Such high-T outcome is typical of classical self-gravitating systems [12] [13] [14].

3.3. Specific Heat C

One defines it as

C = − β T ∂ 〈 U 〉 ∂ β , (3.5)

so that

C = − 1 Z 2 [ π U 0 ϕ ( 1 2 ; 2 ; β U 0 ) + π β U 0 2 4 ϕ ( 3 2 ; 3 ; β U 0 ) ]     × [ π β U 0 T ϕ ( 1 2 ; 2 ; β U 0 ) + π β 2 U 0 2 4 T ϕ ( 3 2 ; 3 ; β U 0 ) ]     + 1 Z [ π β U 0 2 2 T ϕ ( 3 2 ; 3 ; β U 0 ) + π β 2 U 0 3 8 T ϕ ( 5 2 ; 4 ; β U 0 ) ] , (3.6)

depicted in Figure 3. Note that at very low temperatures, quantum effects raise their head and our treatment becomes invalid. Thus, the third thermodynamics’ law is violated here. Interestingly enough there is a Schottky anomaly. This is an effect typical of solid-state physics: the specific heat at low temperature exhibits a peak. When T is high, the specific heat decreases. A new effect is observed at very large T. The specific heat becomes negative. Such outcome is typical of classical self-gravitating systems [12] [13] [15].

4. Moment Generating Functions

We pass to the moment generating function for our extant probability distribution function (PDF) f ( x ) [consult (2.1)]

f ( x ) = e β U 0 1 + x 2 Z , (4.1)

where Z is given by (2.1). In the naive traditional treatment, these moments diverge. The mean value for x 2 n + 1 , ( n = 1 , 2 , 3 , ⋯ ) vanishes by parity. That of x 2 n becomes

〈 x 2 n 〉 = 1 Z ∫ − ∞ ∞     x 2 n e β U 0 1 + x 2 d x . (4.2)

Appeal again to the variables change y = 1 + x 2 and face

〈 x 2 n 〉 = 1 Z ∫ 1 ∞ ( y − 1 ) n − 1 2 e β U 0 y d y , (4.3)

so that, proceeding in a fashion similar to that above we find

〈 x 2 n 〉 = β U 0 Z Γ ( − n + 1 2 ) Γ ( n + 1 2 ) ϕ ( 1 2 − n ; 2 ; β U 0 ) , (4.4)

Thus, we get for the moment generating function M 1 ( t )

M 1 ( t ) = β U 0 Z ∑ n = 0 ∞ t 2 n ( 2 n ) ! Γ ( 1 2 − n ) Γ ( 1 2 + n ) ϕ ( 1 2 − n ; 2 ; β U 0 ) . (4.5)

As particular cases, we obtain the values

〈 x 2 〉 = − π β U 0 Z ϕ ( − 1 2 ; 2 ; β U 0 ) , (4.6)

and

〈 x 4 〉 = π β U 0 Z ϕ ( − 3 2 ; 2 ; β U 0 ) . (4.7)

The first one is plotted in Figure 4. We encounter again here the high temperature effect already reported in [2] [12] [13] (and references therein) and in precedent graphs: a high temperature upper bound, beyond which our treatment becomes invalid. Such bound manifests itself in making negative these types of expectation values at temperatures of the order of 1022 Kelvin. For reference, 100 seconds after the Big Bang it is estimated that the temperature is of a billion K-degrees, and 0.0001 seconds after the Big Bang it is of about T = 1013 K [16].

5. Fisher Information Measure (FIM)

Given a continuous probability distribution function (PDF) f ( x ) with x ∈ Δ ⊂ ℝ and ∫ Δ     f ( x ) d x = 1 , its associated Shannon Entropy S is, as we saw above,

S ( f ) = − ∫ Δ     f ln ( f ) d x (5.1)

a quantifier of global nature that it is not very sensitive to strong changes in the distribution that may take place in a small-sized region. This is not the case for Fisher’s Information Measure (FIM) F [17] [18], which constitutes a quantifier of the gradient content of f ( x ) , being accordingly quite sensitive even to small localized perturbations. One writes

F ( f ) = ∫ Δ 1 f ( x ) [ d f ( x ) d x ] 2 d x = 4 ∫ Δ [ d ψ ( x ) d x ] 2 (5.2)

FIM can be interpreted in variegated fashions. 1) As a quantifier of the ability

to estimate a parameter. 2) As the amount of information that can be extracted from a set of measurements. 3) A quantifier of the state of disorder of a system or phenomenon [18], and finally, at more recent times 4) As a strict measure of order [19] [20] [21]. In the above definition of FIM the division by f ( x ) is not desirable if f ( x ) → 0 at certain x-values. We bypass this issue by working with a real probability amplitudes f ( x ) = ψ 2 ( x ) [17] [18], which is a simpler form (no divisors), while showing that F simply measures the gradient content of ψ ( x ) . The gradient operator significantly influences the contribution of minute local f-changes in FIM’s values. Thus, this quantifier is called a local measure [18].

For the f of (3.5) one has

F ( f ) = 2 Z ∫ 0 ∞     e β U 0 1 + x 2 [ de β U 0 1 + x 2 d x ] 2 d x , (5.3)

or

F ( f ) = 8 β 2 U 0 2 Z ∫ 0 ∞ x 2 ( 1 + x 2 ) 2 e β U 0 1 + x 2 d x . (5.4)

Changing variables in the fashion y = 1 + x 2 we get

F ( f ) = 4 β 2 U 0 2 Z ∫ 1 ∞     y − 2 ( y − 1 ) 1 2 e β U 0 y d y , (5.5)

that after evaluation yields for the Fisher information measure the value

F ( f ) = 2 β U 0 , (5.6)

clearly a very large positive number, given the smallness of the Boltzmann constant entering the denominator. Let us look for the Cramer-Rao (CR) product [ F ( f ) 〈 x 2 〉 f ] , that is always ≥ 1 [18]. The CR relation has been linked to the Heisenberg uncertainty relation (HUR) for the D-dimensional quantum central problem [22]. Still further, Frieden has shown that all UHRs can be derived from the CR relation [18].

We need a value for 〈 x 2 〉 , that we take from (4.4). The Cramer-Rao product 〈 x 2 〉 F is then

F ( f ) 〈 x 2 〉 = − 2 π β 2 U 0 2 Z ϕ ( − 1 2 ; 2 ; β U 0 ) . (5.7)

The CR product is plotted in Figure 5. We see that it is indeed ≥ 1 till we reach a very high temperature, of the order of 1022 Kelvin, at which our probability distribution no longer makes sense. We have already encountered above this effect, in connection with 〈 x 2 〉 -graph, the entropy, and the specific heat.

6. Conclusions

In deceptively simple fashion, we have regularized the partition function for Brownian functions moving in an external potential, thus solving a very old

problem. Some other special cases were already treated by the present authors. One is that of the Z-expression in the case of Newton’s gravity [12], where the divergences are of a different nature from the ones here discussed. A second case is that of Fermi’s problem, cited in the Introduction [23]. Our treatment displays two noticeable features.

· Being of a classical nature, it fails at very low temperatures, where quantum effects become predominant.

· At extremely high temperatures, of the order of 1022 Kelvin, we face a T-upper bound. This fact has already been reported, in another context, by Refs. [12] [13]. Our partition function is saying to us that the system can not exist at such high temperatures.

Summing up: We were here tackling partition function’ divergences, a physically-motivated mathematical problem, that we indeed solved. As for applications, the most we can say at this stage is that we have at our disposal a new canonical probability distribution. Can one use the concomitant partition function Z in a concrete problem? To answer this question, more research is needed. We guess that with this Z some density distribution might be constructed that could describe a quasi-stationary solution in some suitable scenario.

Conflicts of Interest

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

Cite this paper

Plastino, A., Rocca, M.C., Monteoliva, D. and Hernando, A. (2021) Brownian Motion in an External Field Revisited. Journal of Modern Physics, 12, 82-90. https://doi.org/10.4236/jmp.2021.122008

ReferencesFermi, E. (1924) Zeitschrift fur Physik, 26, 54. https://doi.org/10.1007/BF01327311Plastino, A., Rocca, M.C. and Ferri, G.L. (2019) Physica A, 534, 15. https://doi.org/10.1016/j.physa.2019.122169McQuarrie, D.A. (1975) Statistical Mechanics. Harper and Row, New York.Defaveri, L., Anteneodo, C., Kessler, D.A. and Barkai, E. (2020) Physical Review Research, 2, Article ID: 043088.Aghion, E., Kessler, D.A. and Barkai, E. (2019) Physical Review Letters, 122, Article ID: 010601. https://doi.org/10.1103/PhysRevLett.122.010601Aghion, E., Kessler, D.A. and Barkai, E. (2020) Chaos Solitons & Fractals, 138, Article ID: 109890.Ryabov, A., Holubec, V. and Berestneva, E. (2019) Journal of Statistical Mechanics, 2019, Article ID: 084014. https://doi.org/10.1088/1742-5468/ab333fSiler, M., Ornigotti, L., Brzobohat’y, O., J’akl, P., Ryabov, A., Holubec, V. and Filip, R. (2018) Physical Review Letters, 121, Article ID: 230601. https://doi.org/10.1103/PhysRevLett.121.230601Neri, I. (2020) Physical Review Letters, 124, Article ID: 040601. https://doi.org/10.1103/PhysRevLett.124.040601Schwartz, L. (1966) Théorie des distributions. Hermann, Paris.Gradshteyn, I.S. and Rizhik, I.M. (1965) Table of Integrals Series and Products. Academic Press, Cambridge.Plastino, A. and Rocca, M.C. (2018) Physica A, 503, 793. https://doi.org/10.1016/j.physa.2018.03.019Pennini, F., Plastino, A., Rocca, M.C. and Ferri, G.L. (2019) Entropy, 21, 677. https://doi.org/10.3390/e21070677Bountis, T. and Helleman, R.H.G. (1978) Journal of Mathematical Physics, 19, 477. https://doi.org/10.1063/1.523669Lynden-Bell, D. and Lynden-Bell, R.M. (1977) Monthly Notices of the Royal Astronomical Society, 181, 405. https://doi.org/10.1093/mnras/181.3.405Brief History of the Universe. See Web Page. http://www.astro.ucla.edu/wright/BBhistory.htmlFisher, A. (1922) Philosophical Transactions of the Royal Society of London, 222, 309. https://doi.org/10.1098/rsta.1922.0009Roy Frieden, B. (2004) Science from Fisher information: A Unification. Cambridge University Press, Cambridge. https://doi.org/10.1017/CBO9780511616907Pennini, F., Plastino, A. and Ferri, G.L. (2012) Entropy, 14, 2081. https://doi.org/10.3390/e14112081Frieden, B.R. and Hawkins, R.J. (2010) Physical Review E, 82, Article ID: 066117. https://doi.org/10.1103/PhysRevE.82.066117Frieden, B.R., Plastino, A. and Plastino, A.R. (2011) Physica A, 391, 410. https://doi.org/10.1016/j.physa.2011.08.005Dehesa, J.S., Gonzalez-Ferez, R. and Sanchez-Moreno, P. (2007) Journal of Physics A: Mathematical and Theoretical, 40, 1845. https://doi.org/10.1088/1751-8113/40/8/011Plastino, A. and Rocca, M.C. (2019) Physica A, 543, Article ID: 122169. https://doi.org/10.1016/j.physa.2019.122169