The Journal of Chemical Physics

Classical Wigner model based on a Feynman path integral open polymer

S. Karl-Mikael Svensson, Jens Aage Poulsen, Gunnar Nyman

Published Mar 5, 2020• 3 citations

Abstract

The classical Wigner model is one way to approximate the quantum dynamics of atomic nuclei. Here, a new method is presented for sampling the initial quantum mechanical distribution that is required in the classical Wigner model. The new method is tested for the position, position-squared, momentum, and momentum-squared autocorrelation functions for a one-dimensional quartic oscillator and double well potential as well as a quartic oscillator coupled to harmonic baths of different sizes. Two versions of the new method are tested and shown to possibly be useful. Both versions always converge toward the classical Wigner limit. For the one-dimensional cases, some results that are essentially converged to the classical Wigner limit are acquired and others are not far off. For the multi-dimensional systems, the convergence is slower, but approximating the sampling of the harmonic bath with classical mechanics was found to greatly improve the numerical performance. For the double well, the new method is noticeably better than the Feynman-Kleinert linearized path integral method at reproducing the exact classical Wigner results, but they are equally good at reproducing exact quantum mechanics. The new method is suggested as being interesting for future tests on other correlation functions and systems.

DOI: 10.1063/1.5126183

Related papers

Approximate quantum dynamics methods for time correlation functions

Kyle K. G. Smith

The dynamic structure factor of liquid para-hydrogen and ortho-deuterium in corresponding thermodynamic states, (T = 20.0 K, n = 21.24 nm⁻³) and (T = 23.0 K, n = 24.61 nm⁻³) respectively, has been computed by both the Feynman-Kleinert linearized path-integral (FK-LPI) and Ring-Polymer Molecular Dynamics (RPMD) methods and compared with Inelastic X-ray Scattering spectra. The combined use of computational and experimental methods enables a reduction in experimental uncertainties for the determination of the true sample spectrum. Furthermore, the refined experimental spectrum of para-hydrogen and ortho-deuterium is consistently reproduced by both FK-LPI and RPMD at momentum transfers lower than 12.8nm⁻¹. At larger momentum transfers the F K - LP I results agree with experiment much better for ortho-deuterium than for para-hydrogen. More specifically we found that for k ~ 20.0 nm⁻¹ para-hydrogen provides a test case for improved approximations to quantum dynamics. We meet this demand for an improved approximate quantum dynamics method by developing two classes of quasi-classical dynamics that are shown to conserve the initial quantum ensemble when used in conjunction with the Feynman-Kleinert approximation of the density operator. As shown, both classes of dynamics are able to recover the exact classical and high temperature limits of the quantum time correlation function, while a subset is able to recover the exact harmonic limit. A comparison of the approximate quantum time correlation functions obtained from both classes of dynamics are made with the exact results for the challenging model problems of the quartic and double-well potentials. It is found that this new Feynman-Kleinert Quasi-Classical Wigner (FK-QCW) method provides a great improvement over the Feynman-Kleinert implementation of the classical Wigner approximation, also known as FK-LPI, in which purely classical dynamics are used. Furthermore, it is shown that the first class of dynamics reduces to Centroid Molecular Dynamics (CMD) when used within the framework of the classical Wigner approximation for the Kubo transformed time correlation function. Finally, we apply the Feynman-Kleinert Quasi-Classical Wigner (FK- QCW) method to the same liquid para-hydrogen and ortho-deuterium system, previously studied using FK-LPI and RPMD. When applied to this challenging system, it is shown that this new FK-QCW method consistently reproduces the experimental dynamic structure factor for all momentum transfers considered. This shows that FK-QCW provides a great improvement over FK-LPI for not only model problems, but also realistic systems. Furthermore, for small momentum transfers, where RPMD is applicable, it is shown that FK-QCW provides nearly the same results as RPMD, thus suggesting that FK-QCW provides a potentially more appealing algorithm than RPMD since one is not limited to correlation functions involving linear operators. This then suggests that the FK-QCW method is a top contender in the realm of approximate quantum dynamics methods which allow for the practical evaluation of time correlation functions.

May 1, 20140 citations

Semiclassical propagation of Wigner functions

T. Dittrich, E. A. Gómez, Leonardo A. Pachón

We present a comprehensive study of semiclassical phase-space propagation in the Wigner representation, emphasizing numerical applications, in particular as an initial-value representation. Two semiclassical approximation schemes are discussed. The propagator of the Wigner function based on van Vleck's approximation replaces the Liouville propagator by a quantum spot with an oscillatory pattern reflecting the interference between pairs of classical trajectories. Employing phase-space path integration instead, caustics in the quantum spot are resolved in terms of Airy functions. We apply both to two benchmark models of nonlinear molecular potentials, the Morse oscillator and the quartic double well, to test them in standard tasks such as computing autocorrelation functions and propagating coherent states. The performance of semiclassical Wigner propagation is very good even in the presence of marked quantum effects, e.g., in coherent tunneling and in propagating Schrodinger cat states, and of classical chaos in four-dimensional phase space. We suggest options for an effective numerical implementation of our method and for integrating it in Monte-Carlo-Metropolis algorithms suitable for high-dimensional systems.

arXiv:0911.3871The Journal of Chemical PhysicsJun 2, 201036 citations

Momentum distribution functions of strongly correlated systems of particles: Wigner approach and path integrals

Alexander S. Larkin, Vladimir S. Filinov, В. Е. Фортов

The new numerical version of the Wigner approach to quantum mechanics for treatment thermodynamic properties of the strongly interacting systems of particles has been developed for extreme conditions, when there are no small physical parameters and analytical approximations used in different kind of perturbation theories can not be applied. The new path integral representation of the quantum Wigner function in the phase space has been developed for canonical ensemble. Explicit analytical expression of the Wigner function has been obtained in linear and harmonic approximations. The new quantum Monte-Carlo method for calculations of average values of arbitrary quantum operators has been proposed. Preliminary calculations of the momentum distribution function of the Coulomb systems of particles have been carried out. Comparison with classical Maxwell-Boltzmann distribution shows the significant influence of quantum effects on the high energy asymptotics ("tails') of the calculated momentum distribution functions, which resulted in appearance of sharp oscillations.

Journal of Physics Conference SeriesNov 1, 20161 citation

Path integral approach to the Wigner representation of canonical density operators for discrete systems coupled to harmonic baths

Andrés Montoya−Castillo, David R. Reichman

We derive a semi-analytical form for the Wigner transform for the canonical density operator of a discrete system coupled to a harmonic bath based on the path integral expansion of the Boltzmann factor. The introduction of this simple and controllable approach allows for the exact rendering of the canonical distribution and permits systematic convergence of static properties with respect to the number of path integral steps. In addition, the expressions derived here provide an exact and facile interface with quasi- and semi-classical dynamical methods, which enables the direct calculation of equilibrium time correlation functions within a wide array of approaches. We demonstrate that the present method represents a practical path for the calculation of thermodynamic data for the spin-boson and related systems. We illustrate the power of the present approach by detailing the improvement of the quality of Ehrenfest theory for the correlation function Czz(t)=Re⟨σz(0)σz(t)⟩ for the spin-boson model with systematic convergence to the exact sampling function. Importantly, the numerically exact nature of the scheme presented here and its compatibility with semiclassical methods allows for the systematic testing of commonly used approximations for the Wigner-transformed canonical density.

arXiv:1610.04243The Journal of Chemical PhysicsJan 12, 201717 citations

Polymer density functional approach to efficient evaluation of path integrals

Andrei Broukhno, P. N. Vorontsov‐Velyaminov, Henrik Georg Bohr

A polymer density functional theory (P-DFT) has been extended to the case of quantum statistics within the framework of Feynman path integrals. We start with the exact P-DFT formalism for an ideal open chain and adapt its efficient numerical solution to the case of a ring. We show that, similarly, the path integral problem can, in principle, be solved exactly by making use of the two-particle pair correlation function (2p-PCF) for the ends of an open polymer, half of the original. This way the exact data for one-dimensional quantum harmonic oscillator are reproduced in a wide range of temperatures. The exact solution is not, though, reachable in three dimensions (3D) because of a vast amount of storage required for 2p-PCF. In order to treat closed paths in 3D, we introduce a so-called "open ring" approximation which proves to be rather accurate in the limit of long chains. We also employ a simple self-consistent iteration so as to correctly account for the interparticle interactions. The algorithm is speeded up by taking convolutions with the aid of fast Fourier transforms. We apply this approximate path integral DFT (PI-DFT) method to systems within spherical symmetry: 3D harmonic oscillator, atoms of hydrogen and helium, and ions of He and Li. Our results compare rather well to the known data, while the computational effort (some seconds or minutes) is about 100 times less than that with Monte Carlo simulations. Moreover, the well-known "sign problem" is expected to be considerably reduced within the reported PI-DFT, since it allows for a direct estimate of the corresponding partition functions.

Physical Review EOct 3, 20054 citations

Stochastic Solution to Quantum Dynamics

Sarah John, Walter John

The quantum Liouville equation in the Wigner representation is solved numerically by using Monte Carlo methods. For incremental time steps, the propagation is implemented as a classical evolution in phase space modified by a quantum correction. The correction, which is a momentum jump function, is simulated in the quasi-classical approximation via a stochastic process. In this paper the technique, which is developed and validated in two- and three-dimensional momentum space, extends an earlier one-dimensional work. Also, by developing a new algorithm, the application to bound state motion in an anharmonic quartic potential shows better agreement with exact solutions in two-dimensional phase space. I. Introduction Nuclear interaction theory is formulated in the language of quantum mechanics, and hence the development of general methods of solutions to quantum dynamics will provide useful tools for application to a large class of problems in nuclear many-body theory. Different approaches...

NASA STI/Recon Technical Report NFeb 1, 20030 citations