Next Article in Journal
Spin Polarization of Electrons in Two-Color XUV + Optical Photoionization of Atoms
Previous Article in Journal
Ultradilute Quantum Droplets in the Presence of Higher-Order Quantum Fluctuations
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

On the Wigner-Kirkwood Expansion of the Free Energy and the Evaluation of the Quantum Correction

by
Luc Kazandjian
1,
François Soubiran
1,2 and
Jean-Christophe Pain
1,2,*
1
CEA, DAM, DIF, F-91297 Arpajon, France
2
Laboratoire Matière en Conditions Extrêmes, Université Paris-Saclay, CEA, F-91680 Bruyères-le-Châtel, France
*
Author to whom correspondence should be addressed.
Atoms 2022, 10(2), 65; https://doi.org/10.3390/atoms10020065
Submission received: 20 February 2022 / Revised: 12 May 2022 / Accepted: 13 June 2022 / Published: 20 June 2022
(This article belongs to the Section Atomic, Molecular and Nuclear Spectroscopy and Collisions)

Abstract

:
The Wigner-Kirkwood expansion of the quantum correction to the classical free energy is generally said to be in powers of 2 and only its first few terms are presented. In this work, we use the Bloch differential equation to obtain a general description of all terms in a dimensionless form. The first corrective term turns out to be proportional to the product of λ 2 / a 2 , where λ is the thermal de Broglie wavelength and a 3 is the volume per particle, by an effective coupling constant. This dimensionless parameter can be used to assess the magnitude of the quantum correction. Using the one-component plasma as an illustration we highlight the importance of the magnitude of the potential on the quantum correction. The results presented are not formally new; the emphasis is placed on a simple and didactic presentation.

1. Introduction

Describing the physical properties of matter under a wide range of thermodynamic conditions requires a very diverse set of tools. At low density one can easily use classical approximations, even though ionization effects have to be properly included as temperature increases. But as the density increases one needs to include N-body interactions as well as quantum effects. For low temperatures degeneracy effects become also important. The most complicated systems are probably the ones involving condensed and warm dense matter for which electrons are degenerate and even the nuclei can have a quantum behavior.
Quantum Monte-Carlo (QMC) simulations [1] or path-integral molecular dynamics (PIMD) [2] are able to tackle such systems but with a very high computer cost. At higher temperatures one can advantageously consider nuclei as classical and use Kohn-Sham-density-functional-theory-based molecular dynamics (see references in [3]), path-integral Monte-Carlo [4,5] or even pseudo-atom molecular dynamics [6]. But these techniques are still very computer-time-consuming and it is important to determine what is the proper physics to include and to determine the extent of N-body and quantum effects so as to choose the appropriate tool. The N-body effects can be advantageously anticipated by the consideration of coupling parameters. For the quantum behaviors it is a bit more difficult to properly anticipate their importance.
The thermodynamic properties of a system can be obtained from the partition function Q defined in the canonical ensemble by [7]
Q = Tr ( e β H ^ ) ,
where H ^ is the Hamiltonian operator of the system, Tr denotes the trace of the argument, and β = 1 / k B T , where T is the temperature of the system and k B is Boltzmann’s constant. In the classical limit of Boltzmann statistics, the partition function can be simplified and taken equal to [7]
Q cl = 1 h 3 N N ! exp ( β H ) d r d p ,
where h is Planck’s constant, N is the number of atoms of the system, m is the mass of the atoms (they are assumed to be identical for the sake of simplicity), H = p 2 / 2 m + V ( r ) , r = ( r 1 , , r N ) is the 3N-dimensional position vector in configuration space, V ( r ) is the potential energy of the system, and p = ( p 1 , . . . , p N ) is the classical 3N-dimensional momentum. V ( r ) is a regular function of the position coordinates.
For a nearly classical system, if the fermion or boson exchange effects between quantum atoms are neglected, Q can be expressed through the so-called Wigner-Kirkwood expansion [7,8,9,10,11]. The latter is usually said to be in powers of 2 [7,10], where is h / ( 2 π ) , or equivalently of λ 2 where λ = β / m is the thermal de Broglie wavelength [11,12]. It is generally stated that the quantum effects are characterized by the dimensionless ratio λ / l where l is a typical length of the classical particle system [11] or the average radius of an atom [7]. Besides, the Wigner-Kirkwood expansion is obtained with tedious calculations. The aim of the present work is to use simple calculations to deduce the general form of the Wigner-Kirkwood expansion and of its first corrective term. The interest of the latter term is that it enables one to assess the quantum correction.

2. General Form of the 2 Expansion of the Quantum Correction

Let = ( 1 , . . . . , N ) be the gradient in the 3N-dimensional configuration space. The system is assumed to be confined to a cubic volume V . We define a and r’ by
a 3 = V N ,
and
r = r a .
In the framework of Boltzmann statistics, the partition function Q can be expressed as [11]
Q = 1 N ! r | e β H ^ | r d r ,
where we have used the bra-ket notation and where H ^ is the Hamiltonian operator which is expressed in | r representation by
H ^ = 2 2 m 2 + V ( r ) .
By using the closure relation in | p representation in Equation (5), we find
Q = 1 N ! p | r r | e β H ^ | p d p d r ,
or equivalently
Q = 1 N ! h 3 N e ( i / ) p · r g ( p , r ) d p d r ,
with
g ( p , r ) = e β H ^ e ( i / ) p · r .
In Equation (9), H ^ is the Hamiltonian operator expressed in | r representation and defined in Equation (6); exp ( β H ^ ) is applied to exp ( i p · r / ) regarded as a scalar function of r so that g ( p , r ) is a scalar function of r . We define w ( p , r , β ) by [7]
g ( p , r ) = e β H e ( i / ) p · r w ( p , r , β ) ,
where H is defined in the introduction. g ( p , r ) verifies the Bloch differential equation [7]
g β = H ^ [ g ] .
Equations (10) and (11) imply
w β = 2 2 m e β V 2 ( w e β V ) + i m e β V p · ( w e β V ) .
If Q is assumed to be an analytic function of , it is possible to look for w as an expansion in powers of [7]:
w ( p , r , β ) = l = 0 l w l ( p , r , β ) .
Note that the analyticity of Q as a function of does not hold for any potential [13]. We define
p = β m p .
We substitute Equation (13) into Equation (12) and equate the coefficients of the same powers of . With the obvious condition w ( β = 0 ) = 1 , we get w 0 = 1 , w 1 = i ( λ / 2 a ) p · β V , where the superscript ’ indicates that the differentiation is carried out with respect to r ; we also obtain the following relation
β β ( w l + 2 l + 2 ) = 1 2 λ 2 a 2 e β V 2 ( w l l e β V ) + i λ a e β V p · ( w l + 1 l + 1 e β V ) .
It can be shown with Equation (15) that each w l l is a finite linear combination of terms which are proportional to β V differentiated l and only l times with respect to r , to ( λ / a ) l , and possibly to products of l components of p . Note that, for l 1 , all terms of this finite linear combination depend on powers of β with positive integer as exponent. Note also that the number of such terms depends on N. Let us assume that, up to some integer l + 1 , w l + 1 l + 1 verifies the above property. Then the right-hand side of Equation (15) also verifies this property and, since each term depends on β as an integer power of β , so does w l + 2 l + 2 obtained by integration of the right-hand side of Equation (15) under the constraint w l + 2 ( β = 0 ) = 0 . The property stated about w l l is therefore shown by induction and we can write
w l l = λ a l F l ( β V , p , N ) ,
where F l is a dimensionless function which depends on β V , p , and N only, with the dependence on β V and p described above. By substituting Equations (10), (13) and (16) into the right-hand side of Equation (8), we find
Q Q c l = 1 + L = 1 A 2 L [ β V ; N ] λ 2 a 2 L ,
with
A 2 L [ β V ; N ] = F 2 L ( β V , p , N ) c ,
where the subscript c designates the canonical average [7], and the result w 0 = 1 has been taken into account. The terms in ( λ / a ) l with l odd are equal to 0 because Q / Q c l is an even function of . A 2 L is a dimensionless functional of β V ( r ) depending on N only; the dependence of A 2 L on N comes from both the dependence of F 2 L on N and from the bounds of the integrals on r in the calculation of the partition function (all coordinates x j vary between 0 and N 1 / 3 ). Besides, as F 2 L is a linear combination of separable functions of p and r , A 2 L depends linearly on canonical averages of 2 L differentiations of β V with respect to r . Thus, for instance, A 2 can only be of the form
A 2 [ β V ; N ] = α 1 2 β V c + α 2 ( β V ) 2 c ,
where α 1 and α 2 are dimensionless real constants.
In practice, we are interested in the system thermodynamic properties, that is, in the canonical ensemble, in the free energy F, which is equal to [7]
F = k B T ln Q .
In view of Equations (17) and (20), F is equal to
F = F c l k B T ln 1 + L = 1 A 2 L [ β V ; N ] λ 2 a 2 L ,
where F c l = k B T ln Q c l is the free energy in the classical limit of Boltzmann statistics. Equation (21) shows that F can be expressed as an expansion in powers of λ 2 / a 2
F = F c l + k B T L = 1 B 2 L [ β V ; N ] λ 2 a 2 L ,
where B 2 L is a dimensionless functional of β V ( r ) depending on N only. It remains to be found how B 2 L varies with N. As F is an extensive quantity and λ 2 / a 2 is an intensive quantity, the B 2 L ’s must be extensive quantities. As a result, Equation (22) implies
F N k B T = F c l N k B T + L = 1 C 2 L [ β V ] λ 2 a 2 L ,
where C 2 L is a dimensionless functional of β V ( r ) depending on no parameter.

3. Discussion

In Equation (23), the quantum correction depends on through λ 2 / a 2 only. Therefore the 0 limit is in fact the λ 2 / a 2 0 limit. This limit can be obtained by decreasing mass density or by increasing temperature. It can be noted that, in order for the terms in the series in Equation (23) to decrease with increasing L in the λ 2 / a 2 0 limit, the decrease of ( λ 2 / a 2 ) L with L must not be compensated or dominated by the possible increase of C L ( β V ) in this limit. In the present discussion, we consider systems such that, in the λ 2 / a 2 0 limit, the terms in the series decrease with increasing L. It is difficult to find all the sufficient conditions for the existence of such systems. An obvious sufficient condition is for instance the independence of the function V ( r ) with respect to temperature and mass density.
The term of first order in λ 2 / a 2 in Equation (23) is of practical interest. If this term is of the same order or greater than the term of order 0, F c l / N k B T , the λ 2 / a 2 0 limit is not obtained; the quantum correction is large and cannot be evaluated by a perturbative expansion. If the term of first order is small compared to F c l / N k B T , it can indicate that the quantum correction is negligible or can possibly be used as an evaluation of the quantum correction.
We consider the term of first order in λ 2 / a 2 in Equation (23). By noting that, for any regular function f ( r ) ,
2 e β f = β e β f 2 f + β 2 e β f ( f ) 2 ,
and that 2 e β f gives no volume contribution when integrated over r , A 2 can also be written
A 2 [ β V ; N ] = α 3 2 β V c ,
where α 3 is a dimensionless real constant. It can be found with tedious calculations that α 3 = 1 / 24 [11]. As all particles are equivalent, A 2 can also be expressed by
A 2 [ β V ; N ] = α 3 N j 2 β V c ,
where j designates any arbitrary particle. Equations (21)–(23) show that C 2 = B 2 / N = A 2 / N ; therefore Equations (23) and (26) show that, at order 1 in λ 2 / a 2 , F / N k B T is equal to
F N k B T = F c l N k B T + ξ ,
with
ξ = α 3 λ 2 a 2 Γ ,
and
Γ = β a 2 j 2 V c .
In Equation (29), j 2 has been replaced by j 2 by multipying by a 2 . Note that the parameter ξ is an adequate evaluation of the quantum correction for thermodynamics. Let C v be the heat capacity of the system. If, as a first approximation, j 2 V c is assumed to be independent of T, it can be simply deduced from Equation (27) that the quantum correction to C v / N k B is of the same order as ξ .
We now interpret the parameter Γ . If, as a first approximation, the particles are regarded as having an average oscillation length equal to a under the influence of a potential V which is mainly a quadratic function of the coordinates of the particles, j 2 V c a 2 can be regarded as a measure of the potential energy of the individual particle j. Then β j 2 V c a 2 , or equivalently Γ in Equation (29), can be interpreted as a coupling constant. This coupling constant depends in particular on the typical length of variation of V. Note that, in view of Equation (24), Γ can also be found proportional to ( j V ) 2 c = F j 2 c , where F j is the force exerted on particle j. This latter form allows the practical calculation of Γ , for instance in a molecular dynamics simulation. The fact that the coupling constant influences the magnitude of the quantum correction is not a surprise because the latter stems from the non-zero commutator of the operators β K ^ and β V ^ , where K ^ and V ^ are respectively the kinetic energy operator and the potential energy operator [7].

4. Application to the One-Component Plasma

As an illustration, we consider the specific case of the one-component plasma (OCP) which consists of a single species of point ions with charge Z e , where Z is the atomic number and e is the proton charge, immersed in a uniform neutralizing electronic background [14]. The potential energy of the system comes from the Coulombic interaction among ions and from the Coulombic interaction between ions and the neutralizing background; it can therefore be expressed by
V = 1 2 i j = 1 N Z 2 e 2 4 π ϵ 0 r i r j V n 0 i = 1 N Z e 2 4 π ϵ 0 r i r e d r e ,
where n 0 = N Z / V and ϵ 0 is the dielectric permittivity of vacuum. From Equation (30), it follows that [11]
j 2 V c = n 0 Z e 2 ϵ 0 .
Let ξ O C P be ξ for the OCP. By using Equation (31) in Equation (28), we find
ξ O C P = π 6 λ 2 a 2 Γ O C P ,
where α 3 = 1 / 24 has been taken into account and where
Γ O C P = β Z 2 e 2 4 π ϵ 0 a ,
is the usual definition of the coupling constant in the OCP. As interpreted above, ξ O C P is proportional to the coupling constant. This expression of ξ O C P is identical to that found in the literature and expressed as [14,15,16]
ξ O C P = 1 24 ( β ω p ) 2 ,
where ω p is the plasma frequency defined by
ω p 2 = N V Z 2 e 2 ϵ 0 m .
It is worth noting that, in Ref. [15], the authors also derive the λ 4 / a 4 term in an explicit form for both the liquid and solid one-component plasmas. The quantum correction is dependent on λ 2 / a 2 but also on the relative intensity of the potential energy compared to the kinetic energy. To illustrate the importance of the potential, we plotted in Figure 1 the evolution of the quantum correction to the free energy for a fully ionized hydrogen plasma (an OCP) as a function of density and temperature. The quantum corrections are of course very small at low density and high temperature but become the dominant contribution to the free energy at low temperature and high density. There, they cannot be considered as “corrections” anymore. On this graph it is clear that the evolution of the quantum correction is the same as that of the ξ OCP parameter at high temperature and low density. However, the λ / a ratio does not reproduce the density-temperature dependence of the quantum correction. It is thus clear that the “ 0 limit” generally evoked is actually a “ ξ 0 limit”.

5. Conclusions

In the framework of Boltzmann statistics, we have revisited the 2 expansion of the partition function and of the free energy by using simple calculations. We have shown that, in a dimensionless form, the expansion of the free energy is equivalent to an expansion in λ 2 / a 2 whose coefficients are dimensionless functionals of the potential energy (divided by k B T ) depending on no parameter. These coefficients depend in particular (but not only) on the typical length of variation of the potential.
By retaining the first-order corrective term of the expansion in λ 2 / a 2 to evaluate the quantum correction, the latter turns out to be characterized by a dimensionless parameter proportional to the product of λ 2 / a 2 by an effective coupling constant, as shown in Equation (28). As an illustration, the characteristic parameter of the OCP is considered; then the effective coupling constant is proportional to the coupling constant usually defined, as shown in Equation (32). In order to properly evaluate the intensity of the quantum correction, it is thus necessary to consider not only temperature and density but also the magnitude of the potential in which the particles evolve.
An extreme case is that of the ideal gas. Since there is no interaction among atoms, the operators K ^ and V ^ trivially commute so that there is no quantum correction under the assumptions of this work, and this for any thermal de Broglie wavelength. However, the fermion or boson exchange effects between particles, which have been assumed to be negligible in the present work, must now be taken into account in the calculation of the quantum correction; for an ideal gas, the latter varies like λ 3 / a 3 for small values of λ / a [7].
As indicated in the introduction, the main purpose of this work is of a didactic nature. Other approaches allow to obtain the same results. One of them is to deduce w from the identity of the right-hand sides of Equations (9) and (10). Another approach is to consider the Laplace transform of the matrix elements involved in the calculation of Q and to obtain an expansion of the resulting resolvant of the operator H ^ [11]. For obtaining the general form of the expansion of Q and of its coefficients, we have found these two approaches to be much more tedious than the one presented in this work.

Author Contributions

Conceptualization, L.K., F.S. and J.-C.P.; methodology, L.K., F.S. and J.-C.P.; writing—original draft preparation, L.K. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author. The data are not publicly available due to privacy reasons.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
OCPOne-Component Plasma

References

  1. Becca, F.; Sorella, S. Quantum Monte Carlo Approaches for Correlated Systems; Cambridge University Press: Cambridge, UK, 2017. [Google Scholar]
  2. Geneste, G.; Torrent, M.; Bottin, F.; Loubeyre, P. Strong Isotope Effect in Phase II of Dense Solid Hydrogen and Deuterium. Phys. Rev. Lett. 2012, 109, 155303. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Clérouin, J.; Recoules, V.; Soubiran, F. The advent of ab initio simulations of dense plasmas. Contrib. Plasma Phys. 2021, 61, e202100095. [Google Scholar] [CrossRef]
  4. Ceperley, D.M. Path Integral Monte Carlo Methods for Fermions. In Monte Carlo and Molecular Dynamics of Condensed Matter Systems; Binder, K., Ciccotti, G., Eds.; Editrice Compositori: Bologna, Italy, 1996. [Google Scholar]
  5. Militzer, B.; González-Cataldo, F.; Zhang, S.; Driver, K.; Soubiran, F. First-principles equation of state database for warm dense matter computation. Phys. Rev. E 2021, 103, 013203. [Google Scholar] [CrossRef] [PubMed]
  6. Saumon, D.; Starrett, C.E. Pseudo-atom molecular dynamics: A model for warm and hot dense matter. AIP Conf. Proc. 2020, 2272, 090002. [Google Scholar]
  7. McQuarrie, D.A. Statistical Mechanics; University Science Books: Sausalito, CA, USA, 2000. [Google Scholar]
  8. Wigner, E. On the Quantum Correction For Thermodynamic Equilibrium. Phys. Rev. 1932, 40, 749. [Google Scholar] [CrossRef]
  9. Kirkwood, J.G. Quantum Statistics of Almost Classical Assemblies. Phys. Rev. 1933, 44, 31. [Google Scholar] [CrossRef]
  10. Brack, M.; Bhaduri, R.K. Semiclassical Physics; CRC Press: Boca Raton, FL, USA, 2003. [Google Scholar]
  11. Samaj, L.; Jancovici, B. Wigner Kirkwood expansion for semi-infinite quantum fluids. J. Stat. Mech. 2007, 2007, P02002. [Google Scholar] [CrossRef] [Green Version]
  12. Jizba, P.; Zatloukal, V. Path-integral approach to the Wigner-Kirkwood expansion. Phys. Rev. E 2014, 89, 012135. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. DeWitt, H.E. Analytic properties of the quantum corrections to the second virial coefficient. J. Math. Phys. 1962, 3, 1003. [Google Scholar] [CrossRef]
  14. Hansen, J.-P. Statistical Mechanics of Dense Ionized Matter. I. Equilibrium Properties of the Classical One-Component Plasma. Phys. Rev. A 1973, 8, 3096. [Google Scholar] [CrossRef]
  15. Hansen, J.P.; Vieillefosse, P. Quantum corrections in dense ionized matter. Phys. Lett. 1975, 53, 188. [Google Scholar] [CrossRef]
  16. Chabrier, G.; Mazevet, S.; Soubiran, F. A New Equation of State for Dense Hydrogen–Helium Mixtures. Astrophys. J. 2019, 51, 872. [Google Scholar] [CrossRef] [Green Version]
  17. Chabrier, G.; Potekhin, A.Y. Equation of state of fully ionized electron-ion plasmas. Phys. Rev. B 1998, 58, 4941. [Google Scholar] [CrossRef] [Green Version]
  18. Potekhin, A.Y.; Chabrier, G. Thermodynamic Functions of Dense Plasmas: Analytic Approximations for Astrophysical Applications. Contrib. Plasma Phys. 2010, 50, 82. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Temperature-density dependence of the first corrective term Δ F O C P in the Wigner-Kirkwood expansion of the free energy in the OCP approximation for a fully ionized hydrogen plasma. The OCP free energy F OCP corresponds to the ionic classical contribution as defined in Refs. [17,18]. The full lines correspond to iso-values of the parameter ξ OCP and the dashed lines represent iso-values of λ / a .
Figure 1. Temperature-density dependence of the first corrective term Δ F O C P in the Wigner-Kirkwood expansion of the free energy in the OCP approximation for a fully ionized hydrogen plasma. The OCP free energy F OCP corresponds to the ionic classical contribution as defined in Refs. [17,18]. The full lines correspond to iso-values of the parameter ξ OCP and the dashed lines represent iso-values of λ / a .
Atoms 10 00065 g001
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Kazandjian, L.; Soubiran, F.; Pain, J.-C. On the Wigner-Kirkwood Expansion of the Free Energy and the Evaluation of the Quantum Correction. Atoms 2022, 10, 65. https://doi.org/10.3390/atoms10020065

AMA Style

Kazandjian L, Soubiran F, Pain J-C. On the Wigner-Kirkwood Expansion of the Free Energy and the Evaluation of the Quantum Correction. Atoms. 2022; 10(2):65. https://doi.org/10.3390/atoms10020065

Chicago/Turabian Style

Kazandjian, Luc, François Soubiran, and Jean-Christophe Pain. 2022. "On the Wigner-Kirkwood Expansion of the Free Energy and the Evaluation of the Quantum Correction" Atoms 10, no. 2: 65. https://doi.org/10.3390/atoms10020065

APA Style

Kazandjian, L., Soubiran, F., & Pain, J. -C. (2022). On the Wigner-Kirkwood Expansion of the Free Energy and the Evaluation of the Quantum Correction. Atoms, 10(2), 65. https://doi.org/10.3390/atoms10020065

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop