Generalization of classical mechanics for nuclear motions on nonadiabatically coupled potential energy surfaces in chemical reactions

Classical trajectory study of nuclear motion on the Born-Oppenheimer potential energy surfaces is now one of the standard methods of chemical dynamics. In particular, this approach is inevitable in the studies of large molecular systems. However, as soon as more than a single potential energy surface is involved due to nonadiabatic coupling, such a naive application of classical mechanics loses its theoretical foundation. This is a classic and fundamental issue in the foundation of chemistry. To cope with this problem, we propose a generalization of classical mechanics that provides a path even in cases where multiple potential energy surfaces are involved in a single event and the Born-Oppenheimer approximation breaks down. This generalization is made by diagonalization of the matrix representation of nuclear forces in nonadiabatic dynamics, which is derived from a mixed quantum-classical representation of the electron-nucleus entangled Hamiltonian [Takatsuka, K. J. Chem. Phys. 2006, 124, 064111]. A manifestation of quantum fluctuation on a classical subsystem that directly contacts with a quantum subsystem is discussed. We also show that the Hamiltonian thus represented gives a theoretical foundation to examine the validity of the so-called semiclassical Ehrenfest theory (or mean-field theory) for electron quantum wavepacket dynamics, and indeed, it is pointed out that the electronic Hamiltonian to be used in this theory should be slightly modified.     

Nonadiabatic quantum-classical reaction rates with quantum equilibrium structure

Time correlation function expressions for quantum reaction-rate coefficients are computed in a quantum-classical limit. This form for the correlation function retains the full quantum equilibrium structure of the system in the spectral density function but approximates the time evolution of the operator by quantum-classical Liouville dynamics. Approximate analytical expressions for the spectral density function, which incorporate quantum effects in the many-body environment and reaction coordinate, are derived. The results of numerical simulations of the reaction rate are presented for a reaction model in which a two-level system is coupled to a bistable oscillator which is, in turn, coupled to a bath of harmonic oscillators. The nonadiabatic quantum-classical dynamics is simulated in terms of an ensemble of surface-hopping trajectories and the effects of the quantum equilibrium structure on the reaction rate are discussed.     

Bohmian mechanics with complex action: a new trajectory-based formulation of quantum mechanics

In recent years there has been a resurgence of interest in Bohmian mechanics as a numerical tool because of its local dynamics, which suggest the possibility of significant computational advantages for the simulation of large quantum systems. However, closer inspection of the Bohmian formulation reveals that the nonlocality of quantum mechanics has not disappeared-it has simply been swept under the rug into the quantum force. In this paper we present a new formulation of Bohmian mechanics in which the quantum action, S, is taken to be complex. This leads to a single equation for complex S, and ultimately complex x and p but there is a reward for this complexification-a significantly higher degree of localization. The quantum force in the new approach vanishes for Gaussian wave packet dynamics, and its effect on barrier tunneling processes is orders of magnitude lower than that of the classical force. In fact, the current method is shown to be a rigorous extension of generalized Gaussian wave packet dynamics to give exact quantum mechanics. We demonstrate tunneling probabilities that are in virtually perfect agreement with the exact quantum mechanics down to 10(-7) calculated from strictly localized quantum trajectories that do not communicate with their neighbors. The new formulation may have significant implications for fundamental quantum mechanics, ranging from the interpretation of non-locality to measures of quantum complexity.     

Predicting hydration free energies with a hybrid QM/MM approach: an evaluation of implicit and explicit solvation models in SAMPL4

The correct representation of solute-water interactions is essential for the accurate simulation of most biological phenomena. Several highly accurate quantum methods are available to deal with solvation by using both implicit and explicit solvents. So far, however, most evaluations of those methods were based on a single conformation, which neglects solute entropy. Here, we present the first test of a novel approach to determine hydration free energies that uses molecular mechanics (MM) to sample phase space and quantum mechanics (QM) to evaluate the potential energies. Free energies are determined by using re-weighting with the Non-Boltzmann Bennett (NBB) method. In this context, the method is referred to as QM-NBB. Based on snapshots from MM sampling and accounting for their correct Boltzmann weight, it is possible to obtain hydration free energies that incorporate the effect of solute entropy. We evaluate the performance of several QM implicit solvent models, as well as explicit solvent QM/MM for the blind subset of the SAMPL4 hydration free energy challenge. While classical free energy simulations with molecular dynamics give root mean square deviations (RMSD) of 2.8 and 2.3 kcal/mol, the hybrid approach yields an improved RMSD of 1.6 kcal/mol. By selecting an appropriate functional and basis set, the RMSD can be reduced to 1 kcal/mol for calculations based on a single conformation. Results for a selected set of challenging molecules imply that this RMSD can be further reduced by using NBB to reweight MM trajectories with the SMD implicit solvent model.     

Proton transfer in nanoconfined polar solvents. II. Adiabatic proton transfer dynamics

The reaction dynamics for a model phenol-amine proton transfer system in a confined methyl chloride solvent have been simulated by mixed quantum-classical molecular dynamics. In this approach, the proton vibration is treated quantum mechanically (and adiabatically), while the rest of the system is described classically. Nonequilibrium trajectories are used to determine the proton transfer reaction rate constant. The reaction complex and methyl chloride solvent are confined in a smooth, hydrophobic spherical cavity, and radii of 10, 12, and 15 A have been considered. The effects of the cavity radius and the heavy atom (hydrogen bond) distance on the reaction dynamics are considered, and the mechanism of the proton transfer is examined in detail by analysis of the trajectories.     

Theories and Applications of Mixed Quantum-Classical Non-adiabatic Dynamics

Electronically non-adiabatic processes are essential parts of photochemical process, collisions of excited species, electron transfer processes, and quantum information processing. Various non-adiabatic dynamics methods and their numerical implementation have been developed in the last decades. This review summarizes the most significant development of mixed quantum-classical methods and their applications which mainly include the Liouville equation, Ehrenfest mean-field, trajectory surface hopping, and multiple spawning methods. The recently developed quantum trajectory mean-field method that accounts for the decoherence corrections in a parameter-free fashion is discussed in more detail.     

Reconciling semiclassical and Bohmian mechanics. I. Stationary states

The semiclassical method is characterized by finite forces and smooth, well-behaved trajectories, but also by multivalued representational functions that are ill behaved at caustics. In contrast, quantum trajectory methods--based on Bohmian mechanics (quantum hydrodynamics)--are characterized by divergent forces and erratic trajectories near nodes, but also well-behaved, single-valued representational functions. In this paper, we unify these two approaches into a single method that captures the best features of both, and in addition, satisfies the correspondence principle. Stationary eigenstates in one degree of freedom are the primary focus, but more general applications are also anticipated.     

多原子体系非绝热动力学的近似理论方法

非绝热动力学普遍存在于光物理和光化学过程中。描述非绝热跃迁需要处理电子-原子核间的相互耦合运动。由于计算量随体系尺度增大剧烈增长,准确的量子动力学计算目前只适用于描述小分子体系。为了研究多原子分子体系的非绝热过程,近年来发展了一些基于量子-经典动力学近似方法。本文将对典型的这类方法包括经典Ehrenfest方法、面跳跃方法、基于Wigner表示的混合量子-经典方法进行简要的介绍,并讨论如何将量子-经典动力学方法与电子结构从头算手段相结合,模拟非绝热过程。重点阐明各种方法的基本思想和优缺点,并对该领域的发展进行展望。     

WavePacket: A Matlab package for numerical quantum dynamics. III. Quantum-classical simulations and surface hopping trajectories

WavePacket is an open-source program package for numerical simulations in quantum dynamics. Building on the previous Part I (Schmidt and Lorenz, Comput. Phys. Commun. 2017, 213, 223] and Part II (Schmidt and Hartmann, Comput. Phys. Commun. 2018, 228, 229] which dealt with quantum dynamics of closed and open systems, respectively, the present Part III adds fully classical and mixed quantum-classical propagation techniques to WavePacket. There classical phase-space densities are sampled by trajectories which follow (diabatic or adiabatic) potential energy surfaces. In the vicinity of (genuine or avoided) intersections of those surfaces, trajectories may switch between them. To model these transitions, two classes of stochastic algorithms have been implemented: (1) Tully's fewest switches surface hopping and (2) Landau–Zener-based single switch surface hopping. The latter one offers the advantage of being based on adiabatic energy gaps only, thus not requiring nonadiabatic coupling information any more. The present work describes the MATLAB version of WavePacket 6.1.0, which is essentially an object-oriented rewrite of previous versions, allowing to perform fully classical, quantum-classical and quantum-mechanical simulations on an equal footing, that is, for the same physical system described by the same WavePacket input. The software package is hosted and further developed at the Sourceforge platform, where also extensive Wiki-documentation as well as numerous worked-out demonstration examples with animated graphics are available. © 2019 Wiley Periodicals, Inc.     

Quantum trajectories in atom-surface scattering with single adsorbates: the role of quantum vortices

In this work, a full quantum study of the scattering of He atoms off single CO molecules, adsorbed onto the Pt(111) surface, is presented within the formalism of quantum trajectories provided by Bohmian mechanics. By means of this theory, it is shown that the underlying dynamics is strongly dominated by the existence of a transient vortitial trapping with measurable effects on the whole diffraction pattern. This kind of trapping emphasizes the key role played by quantum vortices in this scattering. Moreover, an analysis of the surface rainbow effect caused by the local corrugation that the CO molecule induces on the surface, and its manifestation in the corresponding intensity pattern, is also presented and discussed.     

Long-range electrostatic interactions in hybrid quantum and molecular mechanical dynamics using a lattice summation approach

A method based on a lattice summation technique for treating long-range electrostatic interactions in hybrid quantum mechanics/molecular mechanics simulations is presented in this article. The quantum subsystem is studied at the semiempirical level, whereas the solvent is described by a two-body potential of molecular mechanics. Molecular dynamics simulations of a (quantum) chloride ion in (classical) water have been performed to test this technique. It is observed that the application of the lattice summations to solvent-solvent interactions as well as on solute-solvent ones has a significant effect on solvation energy and diffusion coefficient. Moreover, two schemes for the computation of the long-range contribution to the electrostatic interaction energy are investigated. The first one replaces the exact charge distribution of the quantum solute by a Mulliken charge distribution. The long-range electrostatic interactions are then calculated for this charge distribution that interacts with the solvent molecule charges. The second one is more accurate and involves a modified Fock operator containing long-range electron-charge interactions. It is shown here that both schemes lead to similar results, the method using Mulliken charges for the evaluation of long-range interactions being, however, much more computationally efficient.     

Quantum wavepacket ab initio molecular dynamics: an approach for computing dynamically averaged vibrational spectra including critical nuclear quantum effects

We have introduced a computational methodology to study vibrational spectroscopy in clusters inclusive of critical nuclear quantum effects. This approach is based on the recently developed quantum wavepacket ab initio molecular dynamics method that combines quantum wavepacket dynamics with ab initio molecular dynamics. The computational efficiency of the dynamical procedure is drastically improved (by several orders of magnitude) through the utilization of wavelet-based techniques combined with the previously introduced time-dependent deterministic sampling procedure measure to achieve stable, picosecond length, quantum-classical dynamics of electrons and nuclei in clusters. The dynamical information is employed to construct a novel cumulative flux/velocity correlation function, where the wavepacket flux from the quantized particle is combined with classical nuclear velocities to obtain the vibrational density of states. The approach is demonstrated by computing the vibrational density of states of [Cl-H-Cl]-, inclusive of critical quantum nuclear effects, and our results are in good agreement with experiment. A general hierarchical procedure is also provided, based on electronic structure harmonic frequencies, classical ab initio molecular dynamics, computation of nuclear quantum-mechanical eigenstates, and employing quantum wavepacket ab initio dynamics to understand vibrational spectroscopy in hydrogen-bonded clusters that display large degrees of anharmonicities.