Path-integral centroid dynamics for general initial conditions: A nonequilibrium projection operator formulation

Operator Formulation of Feynman Path Centroid Dynamics for Rotations

An operator formulation of centroid molecular dynamics (CMD) for rotational degrees of freedom is presented. The quasi-density operator concept was introduced by Jang and Voth [J. Chem. Phys 111, 2357 (1999)] and is used to obtain a phase-space mapping without the need for discretized path integrals. The approach allows the calculation of approximate Kubo-transformed time correlation functions. The particle on a ring is chosen as an illustrative example. Numerical results demonstrate that the proposed approach leads to accurate results when compared with exact diagonalization calculations for linear operators. At very low temperatures, it is found that rotational CMD yields results that are in exact agreement with the quantum dynamics of a spin-1 system.

Simple and General Unitarity Conserving Numerical Real-Time Propagators of the Time-Dependent Schrödinger Equation Based on Magnus Expansion

Magnus expansion (ME) provides a general way to expand the real-time propagator of a time-dependent Hamiltonian within the exponential such that the unitarity is satisfied at any order. We use this property and explicit integration of Lagrange interpolation formulas for the time-dependent Hamiltonian within each time interval and derive approximations that preserve unitarity for the differential time evolution operators of general time-dependent Hamiltonians. The resulting second-order approximation is the same as using the average of Hamiltonians for two end points of time. We identify three fourth-order approximations involving commutators of Hamiltonians at different times and also derive a sixth-order expression. A test of these approximations along with other available expressions for a two-state time-dependent Hamiltonian with sinusoidal time dependences provides information on the relative performance of these approximations and suggests that the derived expressions can serve as useful numerical tools for time evolution in time-resolved spectroscopy, quantum control, quantum sensing, real-time ab initio quantum dynamics, and open system quantum dynamics.

Formulation of state projected centroid molecular dynamics: Microcanonical ensemble and connection to the Wigner distribution

A derivation of quantum statistical mechanics based on the concept of a Feynman path centroid is presented for the case of generalized density operators using the projected density operator formalism of Blinov and Roy [J. Chem. Phys. 115, 7822-7831 (2001)]. The resulting centroid densities, centroid symbols, and centroid correlation functions are formulated and analyzed in the context of the canonical equilibrium picture of Jang and Voth [J. Chem. Phys. 111, 2357-2370 (1999)]. The case where the density operator projects onto a particular energy eigenstate of the system is discussed, and it is shown that one can extract microcanonical dynamical information from double Kubo transformed correlation functions. It is also shown that the proposed projection operator approach can be used to formally connect the centroid and Wigner phase-space distributions in the zero reciprocal temperature β limit. A Centroid Molecular Dynamics (CMD) approximation to the state-projected exact quantum dynamics is proposed and proven to be exact in the harmonic limit. The state projected CMD method is also tested numerically for a quartic oscillator and a double-well potential and found to be more accurate than canonical CMD. In the case of a ground state projection, this method can resolve tunnelling splittings of the double well problem in the higher barrier regime where canonical CMD fails. Finally, the state-projected CMD framework is cast in a path integral form.

The properties of water: insights from quantum simulations

The properties of water play a central role in many phenomena of relevance to different areas of science, including physics, chemistry, biology, geology, and climate research. Although well studied for decades, the behavior of water under different conditions and in different environments still remains mysterious and often surprising. In this article, various efforts aimed at providing a comprehensive representation of the water properties at a molecular level through computer modeling and simulation will be described. In particular, the unique role played by the hydrogen-bond network will be examined, first in liquid water, then in the solvation of model biological compounds, and finally in ice, especially highlighting the important effects related to the quantization of the nuclear motion.

Inclusion of inversion symmetry in centroid molecular dynamics: a possible avenue to recover quantum coherence

Inversion symmetry is included in the operator formulation of the centroid molecular dynamics (CMD). This work involves the development of a symmetry-adapted CMD (SA-CMD), here particularly for symmetrization and antisymmetrization projections. A symmetry-adapted quasidensity operator, as defined by Blinov and Roy [J. Chem. Phys. 115, 7822 (2001)], is employed to obtain the centroid representation of quantum mechanical operators. Numerical examples are given for a single particle confined to one-dimensional symmetric quartic and symmetric double-well potentials. Two SA-CMD simulations are performed separately for both projections, and centroid position autocorrelation functions are obtained. For each projection, the quality of the approximation as well as the accuracy are similar to those of regular CMD. It is shown that individual trajectories from two separate SA-CMD simulations can be properly combined to recover trajectories for Boltzmann statistics. Position autocorrelation functions are compared to the exact quantum mechanical ones. This explicit account of inversion symmetry provides a qualitative improvement on the conventional CMD approach and allows the recovery of some quantum coherence.