Semiclassical molecular dynamics of wave packets

Development of a new quantum trajectory molecular dynamics framework

An extension to the wave packet description of quantum plasmas is presented, where the wave packet can be elongated in arbitrary directions. A generalized Ewald summation is constructed for the wave packet models accounting for long-range Coulomb interactions and fermionic effects are approximated by purpose-built Pauli potentials, self-consistent with the wave packets used. We demonstrate its numerical implementation with good parallel support and close to linear scaling in particle number, used for comparisons with the more common wave packet employing isotropic states. Ground state and thermal properties are compared between the models with differences occurring primarily in the electronic subsystem. Especially, the electrical conductivity of dense hydrogen is investigated where a 15% increase in DC conductivity can be seen in our wave packet model compared with other models. This article is part of the theme issue 'Dynamic and transient processes in warm dense matter'.

Implementation of the self-consistent phonons method with ab initio potentials (AI-SCP)

The self-consistent phonon (SCP) method allows one to include anharmonic effects when treating a many-body quantum system at thermal equilibrium. The system is then described by an effective temperature-dependent harmonic Hamiltonian, which can be used to estimate its various dynamic and static properties. In this paper, we combine SCP with ab initio (AI) potential energy evaluation in which case the numerical bottleneck of AI-SCP is the evaluation of Gaussian averages of the AI potential energy and its derivatives. These averages are computed efficiently by the quasi-Monte Carlo method utilizing low-discrepancy sequences leading to a fast convergence with respect to the number, S, of the AI energy evaluations. Moreover, a further substantial (an-order-of-magnitude) improvement in efficiency is achieved once a numerically cheap approximation of the AI potential is available. This is based on using a perturbation theory-like (the two-grid) approach in which it is the average of the difference between the AI and the approximate potential that is computed. The corresponding codes and scripts are provided.

Adiabatic transfer of population in a dense fluid: The role of dephasing statistics

We report the results of simulation studies of the statistics of vibrational dephasing of a YCl (Y=H, D, T, and I) diatom in dense fluid Ar at two temperatures, including the effect of strong field driving on the energy level modulation statistics. The distribution of energy level modulations is found to be non-Gaussian with a high energy tail. Aspects of stimulated Raman adiabatic passage (STIRAP) between the vibrational levels of HCl in dense fluid Ar have been investigated. For HCl with nearly degenerate v=0-->v=1 and v=1-->v=2 transitions, the combined effect of modulation and power broadening reduces the STIRAP efficiency for population transfer from v=0 to v=2 of the order of 30%. However, if the transitions used have very different frequencies, as in the original model studied by Demirplak and Rice [J. Chem. Phys. 116, 8028 (2002)], the STIRAP efficiency for population transfer remains high, of the order of 80%, even with non-Gaussian modulation of energy levels.

Thermodynamics and equilibrium structure of Ne38 cluster: Quantum mechanics versus classical

The equilibrium properties of classical Lennard-Jones (LJ38) versus quantum Ne38 Lennard-Jones clusters are investigated. The quantum simulations use both the path-integral Monte Carlo (PIMC) and the recently developed variational-Gaussian wave packet Monte Carlo (VGW-MC) methods. The PIMC and the classical MC simulations are implemented in the parallel tempering framework. The classical heat capacity Cv(T) curve agrees well with that of Neirotti et al. [J. Chem. Phys. 112, 10340 (2000)], although a much larger confining sphere is used in the present work. The classical Cv(T) shows a peak at about 6 K, interpreted as a solid-liquid transition, and a shoulder at approximately 4 K, attributed to a solid-solid transition involving structures from the global octahedral (Oh) minimum and the main icosahedral (C5v) minimum. The VGW method is used to locate and characterize the low energy states of Ne38, which are then further refined by PIMC calculations. Unlike the classical case, the ground state of Ne38 is a liquidlike structure. Among the several liquidlike states with energies below the two symmetric states (Oh and C5v), the lowest two exhibit strong delocalization over basins associated with at least two classical local minima. Because the symmetric structures do not play an essential role in the thermodynamics of Ne38, the quantum heat capacity is a featureless curve indicative of the absence of any structural transformations. Good agreement between the two methods, VGW and PIMC, is obtained. The present results are also consistent with the predictions by Calvo et al. [J. Chem. Phys. 114, 7312 (2001)] based on the quantum superposition method within the harmonic approximation. However, because of its approximate nature, the latter method leads to an incorrect assignment of the Ne38 ground state as well as to a significant underestimation of the heat capacity.

Quantum statistical mechanics with Gaussians: Equilibrium properties of van der Waals clusters

The variational Gaussian wave-packet method for computation of equilibrium density matrices of quantum many-body systems is further developed. The density matrix is expressed in terms of Gaussian resolution, in which each Gaussian is propagated independently in imaginary time beta=(k(B)T)(-1) starting at the classical limit beta=0. For an N-particle system a Gaussian exp[(r-q)(T)G(r-q)+gamma] is represented by its center qinR(3N), the width matrix GinR(3Nx3N), and the scale gammainR, all treated as dynamical variables. Evaluation of observables is done by Monte Carlo sampling of the initial Gaussian positions. As demonstrated previously at not-very-low temperatures the method is surprisingly accurate for a range of model systems including the case of double-well potential. Ideally, a single Gaussian propagation requires numerical effort comparable to the propagation of a single classical trajectory for a system with 9(N(2)+N)/2 degrees of freedom. Furthermore, an approximation based on a direct product of single-particle Gaussians, rather than a fully coupled Gaussian, reduces the number of dynamical variables to 9N. The success of the methodology depends on whether various Gaussian integrals needed for calculation of, e.g., the potential matrix elements or pair correlation functions could be evaluated efficiently. We present techniques to accomplish these goals and apply the method to compute the heat capacity and radial pair correlation function of Ne(13) Lennard-Jones cluster. Our results agree very well with the available path-integral Monte Carlo calculations.

Semiquantal time-dependent Hartree approach to condensed phase chemical dynamics: application to the system-bath model

A semiquantal analysis of condensed phase chemical dynamics, outlined recently for a double-well linearly coupled to dissipative harmonic bath, is formulated in detail to clarify its general features as well as the specifics of the linear and quadratic coupling cases. The theory may be called a "semiquantal time-dependent Hartree (SQTDH)" approach, as it assumes a factorized product of the squeezed coherent state wave packets for the variational subspace of the many-dimensional time-dependent wave function. Due to this assumption, it straightforwardly satisfies the canonicity condition introduced by Marumori et al. is described by a set of Hamilton equations of motion in an extended phase space that includes auxiliary coordinates representing the wave packet widths. The potential in the extended phase space provides a pictorial understanding of the quantum effects affected due to the bath coupling, e.g., suppression of the wave packet spreading in terms of the potential wall developing along the auxiliary coordinates. The idea is illustrated by prototypical models of quartic double-well and cubic metastable potentials linearly and quadratically coupled to the bath. Further applications and extensions, where the SQTDH method will offer a practical approach for introducing quantum effects into realistic molecular dynamics simulations, are also discussed.

Calculation of infrared absorption spectra using Gaussian variational wave packets

Formulas are developed for calculation of infrared absorption spectra with the help of Gaussian variational wave packets. The excitation is viewed as occurring to a linear combination of basis states corresponding to a Gaussian multiplied by Hermite polynomials of coordinates. Application of the formulas is examplified on HDO isolated in a D2O matrix. A simple approximation is proposed for the fundamental transitions.

A canonical averaging in the second-order quantized Hamilton dynamics

Quantized Hamilton dynamics (QHD) is a simple and elegant extension of classical Hamilton dynamics that accurately includes zero-point energy, tunneling, dephasing, and other quantum effects. Formulated as a hierarchy of approximations to exact quantum dynamics in the Heisenberg formulation, QHD has been used to study evolution of observables subject to a single initial condition. In present, we develop a practical solution for generating canonical ensembles in the second-order QHD for position and momentum operators, which can be mapped onto classical phase space in doubled dimensionality and which in certain limits is equivalent to thawed Gaussian. We define a thermal distribution in the space of the QHD-2 variables and show that the standard beta=1/kT relationship becomes beta'=2/kT in the high temperature limit due to an overcounting of states in the extended phase space, and a more complicated function at low temperatures. The QHD thermal distribution is used to compute total energy, kinetic energy, heat capacity, and other canonical averages for a series of quartic potentials, showing good agreement with the quantum results.