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

: The Wigner-Kirkwood expansion of the quantum correction to the classical free energy is generally said to be in powers of ¯ h 2 and only its ﬁrst 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 ﬁrst 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.


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-Shamdensity-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] whereĤ 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] 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 /2m + 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 ofh 2 [7,10], whereh is h/(2π), or equivalently of λ 2 where λ =h β/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.

General Form of theh 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 and In the framework of Boltzmann statistics, the partition function Q can be expressed as [11] where we have used the bra-ket notation and whereĤ is the Hamiltonian operator which is expressed in |r representation bŷ By using the closure relation in |p representation in Equation (5), we find or equivalently e −(i/h)p·r g(p, r) dp dr, with g(p, r) = e −βĤ e (i/h)p·r .
In Equation (9),Ĥ is the Hamiltonian operator expressed in |r representation and defined in Equation (6); exp(−βĤ) is applied to exp(ip · r/h) 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/h)p·r w(p, r, β), where H is defined in the introduction. g(p, r) verifies the Bloch differential equation [7] ∂g ∂β Equations (10) and (11) imply If Q is assumed to be an analytic function ofh, it is possible to look for w as an expansion in powers ofh [7]: Note that the analyticity of Q as a function ofh does not hold for any potential [13]. We define We substitute Equation (13) into Equation (12) and equate the coefficients of the same powers ofh. With the obvious condition w( where the superscript ' indicates that the differentiation is carried out with respect to r ; we also obtain the following relation It can be shown with Equation (15) that each w lh 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+1h 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+2h l+2 obtained by integration of the right-hand side of Equation (15) under The property stated about w lh l is therefore shown by induction and we can write 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 with 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 cl is an even function ofh. A 2L is a dimensionless functional of βV(r ) depending on N only; the dependence of A 2L on N comes from both the dependence of F 2L 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 2L is a linear combination of separable functions of p and r , A 2L depends linearly on canonical averages of 2L differentiations of βV with respect to r . Thus, for instance, A 2 can only be of the form 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] In view of Equations (17) and (20), F is equal to where F cl = −k B T ln Q cl 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 where B 2L is a dimensionless functional of βV(r ) depending on N only. It remains to be found how B 2L varies with N. As F is an extensive quantity and λ 2 /a 2 is an intensive quantity, the B 2L 's must be extensive quantities. As a result, Equation (22) implies where C 2L is a dimensionless functional of βV(r ) depending on no parameter.

Discussion
In Equation (23), the quantum correction depends onh through λ 2 /a 2 only. Therefore theh → 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 cl /Nk 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 cl /Nk 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), and that ∇ 2 e −β f gives no volume contribution when integrated over r, A 2 can also be written 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 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/Nk B T is equal to with and In Equation (29), ∇ 2 j has been replaced by ∇ 2 j 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, ∇ 2 j V c is assumed to be independent of T, it can be simply deduced from Equation (27) that the quantum correction to C v /Nk 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, ∇ 2 j V c a 2 can be regarded as a measure of the potential energy of the individual particle j. Then β ∇ 2 j 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 2 j 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, whereK andV are respectively the kinetic energy operator and the potential energy operator [7].

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 Ze, 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 where n 0 = NZ/V and 0 is the dielectric permittivity of vacuum. From Equation (30), it follows that [11] ∇ 2 j V c = n 0 Ze 2 0 . (31) Let ξ OCP be ξ for the OCP. By using Equation (31) in Equation (28), we find where α 3 = −1/24 has been taken into account and where is the usual definition of the coupling constant in the OCP. As interpreted above, ξ OCP is proportional to the coupling constant. This expression of ξ OCP is identical to that found in the literature and expressed as [14][15][16] where ω p is the plasma frequency defined by 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 "h → 0 limit" generally evoked is actually a "ξ → 0 limit".  [17,18]. The full lines correspond to iso-values of the parameter ξ OCP and the dashed lines represent iso-values of λ/a.

Conclusions
In the framework of Boltzmann statistics, we have revisited theh 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 operatorsK andV 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 anh expansion of the resulting resolvant of the operatorĤ [11]. For obtaining the general form of theh 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.