A semiclassical hybrid approach to many particle quantum dynamics

We analytically derive a correlated approach for a mixed semiclassical many particle dynamics, treating a fraction of the degrees of freedom by the multitrajectory semiclassical initial value method of Herman and Kluk [Chem. Phys. 91, 27 (1984)] while approximately treating the dynamics of the remaining degrees of freedom with fixed initial phase space variables, analogously to the thawed Gaussian wave packet dynamics of Heller [J. Chem. Phys. 62, 1544 (1975)]. A first application of this hybrid approach to the well studied Secrest-Johnson [J. Chem. Phys. 45, 4556 (1966)] model of atom-diatomic collisions is promising. Results close to the quantum ones for correlation functions as well as scattering probabilities could be gained with considerably reduced numerical effort as compared to the full semiclassical Herman-Kluk approach. Furthermore, the harmonic nature of the different degrees of freedom can be determined a posteriori by comparing results with and without the additional approximation.     

A prefactor free semiclassical initial value series representation of the propagator

A new class of prefactor free semiclassical initial value representations (SCIVR) of the quantum propagator is presented. The derivation is based on the physically motivated demand, that on the average in phase space and in time, the propagator obey the exact quantum equation of motion. The resulting SCIVR series representation of the exact quantum propagator is also free of prefactors. When using a constant width parameter, the prefactor free SCIVR propagator is identical to the frozen Gaussian propagator of Heller [J. Chem. Phys. 75, 2923 (1981)]. A numerical study of the prefactor free SCIVR series is presented for scattering through a double slit potential, a system studied extensively previously by Gelabert et al. [J. Chem. Phys. 114, 2572 (2001)]. As a basis for comparison, the SCIVR series is also computed using the optimized Herman-Kluk SCIVR. We find that the sum of the zeroth order and the first order terms in the series suffice for an accurate determination of the diffraction pattern. The same exercise, but using the prefactor free propagator series needs also the second order term in the series, however the numerical effort is not greater than that needed for the Herman-Kluk propagator, since one does not need to compute the monodromy matrix elements at each point in time. The numerical advantage of the prefactor free propagator grows with increasing dimensionality of the problem.     

A justification for a nonadiabatic surface hopping Herman-Kluk semiclassical initial value representation of the time evolution operator

A justification is given for the validity of a nonadiabatic surface hopping Herman-Kluk (HK) semiclassical initial value representation (SC-IVR) method. The method is based on a propagator that combines the single surface HK SC-IVR method [J. Chem. Phys. 84, 326 (1986)] and Herman's nonadiabatic semiclassical surface hopping theory [J. Chem. Phys. 103, 8081 (1995)], which was originally developed using the primitive semiclassical Van Vleck propagator. We show that the nonadiabatic HK SC-IVR propagator satisfies the time-dependent Schrodinger equation to the first order of variant Planck's over 2pi and the error is O(variant Planck's over 2pi(2)). As a required lemma, we show that the stationary phase approximation, under current assumptions, has an error term variant Planck's over 2pi(1) order higher than the leading term. Our derivation suggests some changes to the previous development, and it is shown that the numerical accuracy in applications to Tully's three model systems in low energies is improved.     

A new time evolving Gaussian series representation of the imaginary time propagator

Frantsuzov and Mandelshtam [J. Chem. Phys. 121, 9247 (2004)] have recently demonstrated that a time evolving Gaussian approximation (TEGA) to the imaginary time propagator exp(-betaH) is useful for numerical computations of anharmonically coupled systems with many degrees of freedom. In this paper we derive a new exact series representation for the imaginary time propagator whose leading order term is the TEGA. One can thus use the TEGA not only as an approximation but also to obtain the exact imaginary time propagator. We also show how the TEGA may be generalized to provide a family of TEGA's. Finally, we find that the equations of motion governing the evolution of the center and width of the Gaussian may be thought of as introducing a quantum friction term to the classical evolution equations.     

Renormalization of the frozen Gaussian approximation to the quantum propagator

The frozen Gaussian approximation to the quantum propagator may be a viable method for obtaining "on the fly" quantum dynamical information on systems with many degrees of freedom. However, it has two severe limitations, it rapidly loses normalization and one needs to know the Gaussian averaged potential, hence it is not a purely local theory in the force field. These limitations are in principle remedied by using the Herman-Kluk (HK) form for the semiclassical propagator. The HK propagator approximately conserves unitarity for relatively long times and depends only locally on the bare potential and its second derivatives. However, the HK propagator involves a much more expensive computation due to the need for evaluating the monodromy matrix elements. In this paper, we (a) derive a new formula for the normalization integral based on a prefactor free HK propagator which is amenable to "on the fly" computations; (b) show that a frozen Gaussian version of the normalization integral is not readily computable "on the fly"; (c) provide a new insight into how the HK prefactor leads to approximate unitarity; and (d) how one may construct a prefactor free approximation which combines the advantages of the frozen Gaussian and the HK propagators. The theoretical developments are backed by numerical examples on a Morse oscillator and a quartic double well potential.     

New coherent state representation for the imaginary time propagator with applications to forward-backward semiclassical initial value representations of correlation functions

There have been quite a few attempts in recent years to provide an initial value coherent state representation for the imaginary time propagator exp(-betaH). The most notable is the recent time evolving Gaussian approximation of Frantsuzov and Mandelshtam [J. Chem. Phys. 121, 9247 (2004)] which may be considered as an expansion of the imaginary time propagator in terms of coherent states whose momentum is zero. In this paper, a similar but different expression is developed in which exp(-betaH) is represented in a series whose terms are weighted phase space averages of coherent states. Such a representation allows for the formulation of a new and simplified forward-backward semiclassical initial value representation expression for thermal correlation functions.     

Herman-Kluk semiclassical dynamics of molecular rotations in laser fields

The action-angle mapping algorithm [R. Saha and M. Ovchinnikov, J. Chem. Phys. 124, 204112 (2006)] is utilized to provide a Herman-Kluk semiclassical initial value representation (SC-IVR) treatment of quantum dynamics of systems with non-Cartesian degrees of freedom. The non-Cartesian system under investigation is a linear rotor molecule in static electric and pulsed laser field. The results demonstrate that the SC-IVR procedure described in this work provides an accurate representation of quantum rotational dynamics of the system.     

Nonadiabatic surface hopping Herman-Kluk semiclassical initial value representation method revisited: applications to Tully's three model systems

The nonadiabatic surface hopping Herman-Kluk (HK) semiclassical initial value representation (SC-IVR) method for nonadiabatic problems is reformulated. The method has the same spirit as Tully's surface hopping technique [J. Chem. Phys. 93, 1061 (1990)] and almost keeps the same structure as the original single-surface HK SC-IVR method except that trajectories can hop to other surfaces according to the hopping probabilities and phases, which can be easily integrated along the paths. The method is based on a rather general nonadiabatic semiclassical surface hopping theory developed by Herman [J. Chem. Phys. 103, 8081 (1995)], which has been shown to be accurate to the first order in h and through all the orders of the nonadiabatic coupling amplitude. Our simulation studies on the three model systems suggested by Tully demonstrate that this method is practical and capable of describing nonadiabatic quantum dynamics for various coupling situations in very good agreement with benchmark calculations.     

A proof of Jarzynski's nonequilibrium work theorem for dynamical systems that conserve the canonical distribution

We present a derivation of the Jarzynski [Phys. Rev. Lett. 78, 2690 (1997)] identity and the Crooks [J. Stat. Phys. 90, 1481 (1998)] fluctuation theorem for systems governed by deterministic dynamics that conserves the canonical distribution such as Hamiltonian dynamics, Nose-Hoover dynamics, Nose-Hoover chains, and Gaussian isokinetic dynamics. The proof is based on a relation between the heat absorbed by the system during the nonequilibrium process and the Jacobian of the phase flow generated by the dynamics.     

Using the thermal Gaussian approximation for the Boltzmann operator in semiclassical initial value time correlation functions

The thermal Gaussian approximation (TGA) recently developed by Frantsuzov et al. [Chem. Phys. Lett. 381, 117 (2003)] has been demonstrated to be a practical way for approximating the Boltzmann operator exp(-betaH) for multidimensional systems. In this paper the TGA is combined with semiclassical (SC) initial value representations (IVRs) for thermal time correlation functions. Specifically, it is used with the linearized SC-IVR (LSC-IVR, equivalent to the classical Wigner model), and the "forward-backward semiclassical dynamics" approximation developed by Shao and Makri [J. Phys. Chem. A 103, 7753 (1999); 103, 9749 (1999)]. Use of the TGA with both of these approximate SC-IVRs allows the oscillatory part of the IVR to be integrated out explicitly, providing an extremely simple result that is readily applicable to large molecular systems. Calculation of the force-force autocorrelation for a strongly anharmonic oscillator demonstrates its accuracy, and calculation of the velocity autocorrelation function (and thus the diffusion coefficient) of liquid neon demonstrates its applicability.     

Interference and quantization in semiclassical response functions

Application of the Herman-Kluk semiclassical propagator to the calculation of spectroscopic response functions for anharmonic oscillators has demonstrated the quantitative accuracy of these approximate dynamics. In this approach, spectroscopic response functions are expressed as multiple phase-space integrals over pairs of classical trajectories and their associated stability matrices. Here we analyze the Herman-Kluk semiclassical approximation to a linear response function and determine the origin of the capacity of this method to reproduce quantum effects in a response function from classical dynamical information. Our analysis identifies those classical trajectories that contribute most significantly to the response function on different time scales. This finding motivates a procedure for computing the linear response function in which the interference between pairs of classical trajectories is treated approximately, resulting in an integral over a single average trajectory, as in a purely classical calculation.     

Operator splitting algorithm for isokinetic SLLOD molecular dynamics

We apply an operator splitting method to develop a simulation algorithm that has complete analytical solutions for the Gaussian thermostated SLLOD equations of motion [D. J. Evans and G. P. Morriss, Phys. Rev. A 30, 1528 (1984)] for a system under shear. This leads to a homogeneous algorithm for performing both equilibrium and nonequilibrium isokinetic molecular dynamics simulation. The resulting algorithm is computationally efficient. In particular, larger integration time steps can be used compared to simulations with regular Gaussian thermostated SLLOD equations of motion. The utility and accuracy of the algorithm are demonstrated through application to the Weeks-Chandler-Anderson fluid. Although strict conservation of the kinetic energy suppresses thermal fluctuations in the system, this algorithm does not allow simulations at lower shear rates than those normally afforded by older nonequilibrium molecular dynamics simulations.     

Finite temperature application of the corrected propagator method to reactive dynamics in a condensed-phase environment

The recently proposed mixed quantum-classical method is extended to applications at finite temperatures. The method is designed to treat complex systems consisting of a low-dimensional quantum part (the primary system) coupled to a dissipative bath described classically. The method is based on a formalism showing how to systematically correct the approximate zeroth-order evolution rule. The corrections are defined in terms of the total quantum Hamiltonian and are taken to the classical limit by introducing the frozen Gaussian approximation for the bath degrees of freedom. The evolution of the primary system is governed by the corrected propagator yielding the exact quantum dynamics. The method has been tested on a standard model system describing proton transfer in a condensed-phase environment: a symmetric double-well potential bilinearly coupled to a bath of harmonic oscillators. Flux correlation functions and thermal rate constants have been calculated at two different temperatures for a range of coupling strengths. The results have been compared to the fully quantum simulations of Topaler and Makri [J. Chem. Phys. 101, 7500 (1994)] with the real path integral method.     

On the properties of a primitive semiclassical surface hopping propagator for nonadiabatic quantum dynamics

A previously developed nonadiabatic semiclassical surface hopping propagator [M. F. Herman J. Chem. Phys. 103, 8081 (1995)] is further studied. The propagator has been shown to satisfy the time-dependent Schrodinger equation (TDSE) through order h, and the O(h2) terms are treated as small errors, consistent with standard semiclassical analysis. Energy is conserved at each hopping point and the change in momentum accompanying each hop is parallel to the direction of the nonadiabatic coupling vector resulting in both transmission and reflection types of hops. Quantum mechanical analysis and numerical calculations presented in this paper show that the h2 terms involving the interstate coupling functions have significant effects on the quantum transition probabilities. Motivated by these data, the h2 terms are analyzed for the nonadiabatic semiclassical propagator. It is shown that the propagator can satisfy the TDSE for multidimensional systems by including another type of nonclassical trajectories that reflect on the same surfaces. This h2 analysis gives three conditions for these three types of trajectories so that their coefficients are uniquely determined. Besides the nonadiabatic semiclassical propagator, a numerically useful quantum propagator in the adiabatic representation is developed to describe nonadiabatic transitions.     

Near-edge x-ray absorption and natural circular dichroism spectra of L-alanine: a theoretical study based on the complex polarization propagator approach

The complex polarization propagator method [J. Chem. Phys. 123, 194103 (2005)] has been employed in conjunction with density functional theory and gauge-including atomic orbitals in order to determine the near-edge x-ray absorption and natural circular dichroism spectra of L-alanine in its neutral and zwitterionic forms. Results are presented for the K-edges of carbon, nitrogen, and oxygen. In contrast to traditional methods, the proposed approach enables a direct determination of the spectra at an arbitrary frequency instead of focusing on the rotatory strengths for individual electronic transitions. The propagator includes a complete set of nonredundant electron-transfer operators and allows for full core-hole relaxation. The theoretical spectrum at the nitrogen K-edge of the zwitterion compares well with the experimental spectrum.     

Imaginary time Gaussian dynamics of the Ar3 cluster

Semiclassical Gaussian approximations to the Boltzmann operator have become an important tool for the investigation of thermodynamic properties of clusters of atoms at low temperatures. Usually, numerically expensive thawed Gaussian variants are applied. In this article, we introduce a numerically much cheaper frozen Gaussian approximation to the imaginary time propagator with a width matrix especially suited for the dynamics of clusters. The quality of the results is comparable to that of thawed Gaussian methods based on the single-particle ansatz. We apply the method to the argon trimer and investigate the dissociation process of the cluster. The results clearly show a classical-like transition from a bounded moiety to three free particles at a temperature T ≈ 20 K, whereas previous studies of the system were not able to resolve this transition. Quantum effects, i.e., differences with the purely classical case manifest themselves in the low-temperature behavior of the mean energy and specific heat as well as in a slight shift of the transition temperature. We also discuss the influence of an artificial confinement of the atoms usually introduced to converge numerical computations. The results show that restrictive confinements often implemented in studies of clusters can influence the thermodynamic properties drastically. This finding may have implications on other studies of atomic clusters.     

Wavepacket propagation using time-sliced semiclassical initial value methods

A new semiclassical initial value representation (SC-IVR) propagator and a SC-IVR propagator originally introduced by Kay [J. Chem. Phys. 100, 4432 (1994)], are investigated for use in the split-operator method for solving the time-dependent Schrodinger equation. It is shown that the SC-IVR propagators can be derived from a procedure involving modified Filinov filtering of the Van Vleck expression for the semiclassical propagator. The two SC-IVR propagators have been selected for investigation because they avoid the need to perform a coherent state basis set expansion that is necessary in other time-slicing propagation schemes. An efficient scheme for solving the propagators is introduced and can be considered to be a semiclassical form of the effective propagators of Makri [Chem. Phys. Lett. 159, 489 (1989)]. Results from applications to a one-dimensional, two-dimensional, and three-dimensional Hamiltonian for a double-well potential are presented.     

Density functional theory of liquid crystals and surface anchoring: hard Gaussian overlap-sphere and hard Gaussian overlap-surface potentials

This article applies the density functional theory to confined liquid crystals, comprised of ellipsoidal shaped particles interacting through the hard Gaussian overlap (HGO) potential. The extended restricted orientation model proposed by Moradi and co-workers [J. Phys.: Condens. Matter 17, 5625 (2005)] is used to study the surface anchoring. The excess free energy is calculated as a functional expansion of density around a reference homogeneous fluid. The pair direct correlation function (DCF) of a homogeneous HGO fluid is approximated, based on the optimized sum of Percus-Yevick and Roth DCF for hard spheres; the anisotropy introduced by means of the closest approach parameter, the expression proposed by Marko [Physica B 392, 242 (2007)] for DCF of HGO, and hard ellipsoids were used. In this study we extend an our previous work [Phys. Rev. E 72, 061706 (2005)] on the anchoring behavior of hard particle liquid crystal model, by studying the effect of changing the particle-substrate contact function instead of hard needle-wall potentials. We use the two particle-surface potentials: the HGO-sphere and the HGO-surface potentials. The average number density and order parameter profiles of a confined HGO fluid are obtained using the two particle-wall potentials. For bulk isotropic liquid, the results are in agreement with the Monte Carlo simulation of Barmes and Cleaver [Phys. Rev. E 71, 021705 (2005)]. Also, for the bulk nematic phase, the theory gives the correct density profile and order parameter between the walls.     

Accurate and efficient method for the treatment of exchange in a plane-wave basis

We describe an accurate and efficient extension of Chawla and Voth's [J. Chem. Phys. 108, 4697 (1998)] plane-wave based algorithm for calculating exchange energies, exchange energy densities, and exchange energy gradients with respect to wave-function parameters in systems of electrons subject to periodic boundary conditions. The theory and numerical results show that the computational effort scales almost linearly with the number of plane waves and quadratically with the number of k vectors. To obtain high accuracy with relatively few k vectors, we use an adaptation of Gygi and Baldereschi's [Phys. Rev. B 34, 4405 (1986)] method for reducing Brillouin-zone integration errors.