General Commentary ARTICLE
Commentary on “Force-field functor theory: classical force-fields which reproduce equilibrium quantum distributions”
- Istituto dei Sistemi Complessi - Consiglio Nazionale delle Ricerche and Istituto Nazionale di Fisica Nucleare - Sezione di Firenze, Sesto Fiorentino, Italy
by Babbush, R., Parkhill, J., and Aspuru-Guzik, A. (2013). Front. Chem. 1:26. doi: 10.3389/fchem.2013.00026
This comment regards a recent paper by Babbush, Parkhill, and Aspuru-Guzik (henceforth BPA) (Babbush et al., 2013). The subject is the equilibrium thermodynamics of a system of many quantum particles with Hamiltonian , where and the commutator of the coordinate- and momentum operators is . In BPA it is correctly observed (see, e.g., Hillery et al., 1984; Cuccoli et al., 1995) that the coordinate distribution function at temperature T = β−1, namely
can be used to define the effective potential W(q) (the kinetic part of the partition function Z is omitted in BPA) in terms of which the exact quantum equilibrium average of any operator takes the classical form of a configuration integral,
It is convincingly proven in BPA that the mapping η(q) ↔ W(q) (i.e., the exponential function) is bijective, as well as V(q) ↔ η(q), and then it follows that V(q) ↔ W(q) is one-to-one. Furthermore, it is correctly pointed out that the Giachetti–Tognetti–Feynman–Kleinert (Giachetti and Tognetti, 1985; Feynman and Kleinert, 1986) (GTFK) effective potential Veff(q) differs from W(q). Indeed, Veff(q), which accounts exactly for any quadratic potential, entails that for approximating one has to include a further Gaussian average accounting for purely-quantum fluctuations, as shown, e.g., in Vaia and Tognetti (1990); Cuccoli et al. (1995): there W(q) is also introduced and dubbed the local effective potential. However, in BPA it is not shown that Equation (29), the paper's main result obtained as the Jensen's approximation W(q) ≈ WBPA(q) to the exact formula (26), can be calculated explicitly as
i.e., the convolution between the potential V(q) and a Gaussian with variance σ2 = βħ2/(6m) proportional to the squared de-Broglie wavelength. This is in agreement with the Wigner series (Wigner, 1932) up to lowest order, but lacks the nonlinear contributions to it. How accurate is the approximation made in Equation (29) of BPA? One can estimate this by considering a single (N = 1) quantum harmonic oscillator, , whose frequency is . Expanding V(q + ξ) in Equation (3) one finds
which does not improve upon the classical result using Equation (2). From the known density for the quantum harmonic oscillator the exact functor for the class of harmonic potentials can be easily derived:
For a linear functor, this expression should be proportional to mω2: evidently this is true only in the classical limit, βħω « 1 or T » ħω, where Equation (4) is recovered. However, the mapping V → W can surely be locally linear, namely V + εδV → W + εδW with δW independent of the small parameter ε. Hence, WBPA(q) is reliable only when the temperature overcomes the typical quantum energy scale ħω; for instance, taking ω2 ~ V″(qm)/m [qm being the minimum of V(q)], a pair of hydrogen molecules has typically ħω ~ 102 K (Vaia and Tognetti, 1990 and references therein) and using WBPA(q) would only be reliable at very high T » 102 K, i.e., just in the classical limit. Such an approximation is indeed used in the high-T propagator of path-integral Monte Carlo algorithms in order to improve convergence in the Trotter number (Takahashi and Imada, 1984). Hence, the approximation (29) of BPA can “reproduce quantum distributions” just when these are almost classical.
On the other hand, the use of the exact effective pair-potential, rather than that obtained from Equation (29) of BPA, is a good starting point for treating a not too dense quantum fluid by means of a classical-like simulation, as shown in the last section of BPA and as noted by several authors (see, e.g., Thirumalai et al., 1984 and many references cited in BPA). At variance with the procedure of BPA, based on the heavy calculation of a (locally) “linear functor” at fixed T, it would be more practical to directly obtain the exact pair-potential W(q) for the chosen , a task that can easily be carried out at any T.
Babbush, R., Parkhill, J. A., and Aspuru-Guzik, A. (2013). Force-field functor theory: classical force-fields which reproduce equilibrium quantum distributions. Front. Chem. 1:26. doi: 10.3389/fchem.2013.00026
Cuccoli, A., Giachetti, R., Tognetti, V., Vaia, R., and Verrucchi, P. (1995). The effective potential and effective Hamiltonian in quantum statistical mechanics. J. Phys. 7, 7891–7938. doi: 10.1088/0953-8984/7/41/003
Giachetti, R., and Tognetti, V. (1985). Variational approach to quantum statistical mechanics of nonlinear systems with application to sine-Gordon chains. Phys. Rev. Lett. 55, 912–915. doi: 10.1103/PhysRevLett.55.912
Keywords: quantum thermodynamics, effective classical potential, path integration, Monte Carlo simulations, force-field functor
Citation: Vaia R (2013) Commentary on “Force-field functor theory: classical force-fields which reproduce equilibrium quantum distributions” Front. Chem. 1:34. doi: 10.3389/fchem.2013.00034
Received: 29 October 2013; Accepted: 23 November 2013;
Published online: 24 December 2013.
Edited by:Francesco Aquilante, Uppsala University, Sweden
Copyright © 2013 Vaia. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.