利用基于直接动力学的轨线面跳跃方法处理非绝热过程

势能面交叉引起的非绝热过程广泛存在于光化学和光物理中。对这一过程进行描述是理论化学的重要挑战之一。非绝热过程涉及原子核与电子之间的耦合运动,因此量子化学的基本假设之一"玻恩-奥本海默"近似被打破,所以对其进行描述需要发展新的动力学理论方法。在这些方法中,Tully发展的最少轨线面跳跃方法凭借易于程序化、便于计算等优点已经发展成为处理非绝热问题的主要动力学方法之一。其中原子核以经典的方式在单一势能面上进行演化,电子以量子的方式沿着同一轨线进行演化。在整个演化过程中,非绝热跃迁通过轨线在不同势能面间的跃迁来描述,其中跳跃发生的几率与电子的演化有关。如果将该方法与从头算直接动力学相结合,可以在全原子水平上研究实际分子体系的非绝热动力学,给出其激发态寿命、非绝热动力学中分子的主要运动方式、反应通道以及分支比等重要信息。本文旨在讨论最少面跳跃直接动力学方法研究非绝热问题的一些进展,包括动力学基本理论,特别关注将最少面跳跃方法和直接动力学结合的数值实现细节,同时讨论该方法在研究实际体系当中的一些应用,并对轨线面跳跃方法下一步发展的一些方向进行合理的展望。     

Quantum manifestations of classical trajectories in molecular systems

Quantum chaos, understood as the effect of the underlying classical dynamics on the stationary quantum properties in classically chaotic systems, is examined in two molecular floppy systems. Realistic models of two degrees of freedom for HO₂ and HCN/HNC are considered. The structure of the classical phase space is studied using Poincaré surfaces of section and the dynamical characteristics of the corresponding wave functions analyzed also in phase space with the aid of Husimi functions. Some wave functions show strong localization along periodic orbits. © 2002 John Wiley & Sons, Inc. Int J Quantum Chem, 2001     

Non-adiabatic molecular dynamics with quantum solvent effects

Three novel approaches extending quantum-classical non-adiabatic (NA) molecular dynamics (MD) to include quantum effects of solvent environments are described. In a standard NA-MD the solute subsystem is treated quantum mechanically, while the larger solvent part of a system is treated classically. The three novel approaches presented here are based on the Bohmian formulation of quantum mechanics, the stochastic Schrödinger equation for the evolution of open quantum systems and the quantized Hamilton dynamics generalization of classical mechanics. The approaches extend the standard NA-MD to incorporate the following quantum effects of the solvent. (1) Branching, i.e. the ability of solvent quantum wave packets to split and follow asymptotically diverging trajectories correlated with different quantum states of the solute. (2) Decoherence, i.e. loss of quantum interference within the solute subsystem induced by the diverging solvent trajectories. (3) Zero point energy that contributes to NA coupling and must be preserved during the energy exchange between solvent and solute degrees of freedom. The Bohmian quantum-classical mechanics, stochastic mean-field and quantized mean-field approximations incorporate the quantum solvent effects into the standard quantum-classical NA-MD in a straightforward and efficient way that can be easily applied to quantum dynamics of condensed phase chemical systems.     

Reduced density matrix hybrid approach: an efficient and accurate method for adiabatic and non-adiabatic quantum dynamics

We present a new approach to calculate real-time quantum dynamics in complex systems. The formalism is based on the partitioning of a system's environment into "core" and "reservoir" modes with the former to be treated quantum mechanically and the latter classically. The presented method only requires the calculation of the system's reduced density matrix averaged over the quantum core degrees of freedom which is then coupled to a classically evolved reservoir to treat the remaining modes. We demonstrate our approach by applying it to the spin-boson problem using the noninteracting blip approximation to treat the system and core, and Ehrenfest dynamics to treat the reservoir. The resulting hybrid methodology is accurate for both fast and slow baths, since it naturally reduces to its composite methods in their respective regimes of validity. In addition, our combined method is shown to yield good results in intermediate regimes where neither approximation alone is accurate and to perform equally well for both strong and weak system-bath coupling. Our approach therefore provides an accurate and efficient methodology for calculating quantum dynamics in complex systems.     

A structurally driven spectral transform Lanczos method for mixed quantum—classical canonical simulations (MQCC): the thermodynamics of Kr10–H and the energies of Krn–H, (n:1→9)

We develop and test an approximate approach for canonical simulations of weakly bound atomic or molecular systems for which some degrees of freedom can be treated separately by quantum mechanics. The system chosen for testing is Kr₁₀–H, for which the adiabatic approximation applied to separate the hydrogen degrees of freedom works reasonably well. The hydrogen atom is bound to the Kr clusters at cold temperatures and we calculate several bound states for clusters in the n=1–9 range, in the global minimum configuration. The structural character of the mixed quantum classical simulation is substantially different than the classical simulation for Kr₁₀–H as a result of zero point energy effects. When quantum effects are included, the low temperature dynamics of Kr₁₀–H are dominated by a significant well to well hopping about an incomplete icosahedral krypton core.     

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

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

Independent-cluster methods as mappings of quantum theory into classical mechanics

Summary The algebraic structures of theconfiguration interaction, normal coupled cluster, andextended coupled cluster methods are reviewed and developed. These methods are pointed out to perform a mapping of the quantum mechanical problem into a classical phase space, where in each case the classical canonical coordinates have characteristically different cluster and locality properties. Special focus is given to the extended coupled cluster method (ECCM), which alone is based on an entirely additively separable coordinate system. The general principles are formulated for systems with both bosonic and fermionic degrees of freedom, allowing both commutative and anticommutative (Grassmann) cluster amplitudes. The properties of the classical images are briefly discussed. It is proposed that phase spaces may exist which are fixed points of quantization.Based on a talk given at theWorkshop on Coupled-Cluster Theory at the Interface of Atomic Physics and Quantum Chemistry, Harvard-Smithsonian Institute for Theoretical Atomic and Molecular Physics, Cambridge, MA, August 7–11, 1990     

Surface-hopping dynamics and decoherence with quantum equilibrium structure

In open quantum systems, decoherence occurs through interaction of a quantum subsystem with its environment. The computation of expectation values requires a knowledge of the quantum dynamics of operators and sampling from initial states of the density matrix describing the subsystem and bath. We consider situations where the quantum evolution can be approximated by quantum-classical Liouville dynamics and examine the circumstances under which the evolution can be reduced to surface-hopping dynamics, where the evolution consists of trajectory segments exclusively evolving on single adiabatic surfaces, with probabilistic hops between these surfaces. The justification for the reduction depends on the validity of a Markovian approximation on a bath averaged memory kernel that accounts for quantum coherence in the system. We show that such a reduction is often possible when initial sampling is from either the quantum or classical bath initial distributions. If the average is taken only over the quantum dispersion that broadens the classical distribution, then such a reduction is not always possible.     

Quantized Hamilton Dynamics

The paper describes the quantized Hamilton dynamics (QHD) approach that extends classical Hamiltonian dynamics and captures quantum effects, such as zero point energy, tunneling, decoherence, branching, and state-specific dynamics. The approximations are made by closures of the hierarchy of Heisenberg equations for quantum observables with the higher order observables decomposed into products of the lower order ones. The technique is applied to the vibrational energy exchange in a water molecule, the tunneling escape from a metastable state, the double-slit interference, the population transfer, dephasing and vibrational coherence transfer in a two-level system coupled to a phonon, and the scattering of a light particle off a surface phonon, where QHD is coupled to quantum mechanics in the Schrödinger representation. Generation of thermal ensembles in the extended space of QHD variables is discussed. QHD reduces to classical mechanics at the first order, closely resembles classical mechanics at the higher orders, and requires little computational effort, providing an efficient tool for treatment of the quantum effects in large systems.     

The infrared-driven cis–trans isomerization of nitrous acid HONO III: A mixed quantum–classical simulation

A mixed quantum–classical simulation of the IR-driven cis–trans isomerization of HONO in a Kr matrix at 30 K is presented, treating the hydrogen atom as quantum particle and the Kr matrix as well as intermolecular degrees of freedom of the ONO-body as classical. A new method is presented to time-propagate the coupled set of equations in a DVR basis in internal spherical coordinates, rather than in laboratory frame fixed cartesian coordinates. In spherical coordinates, a much more precise computation of the weak vibrational couplings is possible using a still manageable basis size. Good qualitative agreement between simulation and experiment is obtained, underestimating relaxation and isomerization rates by a modest factor ≈5. Upon matrix fluctuations, frequent curve crossings occur between the initially excited OH-stretch vibration and a closely lying combination mode of torsional and bending coordinate that lead to a transfer of population. The subsequent pathway of energy flow is deduced and discussed within a tier model, where trans-states, that belong to the second tier, are populated through a first tier of states that is composed of combinations of bending and torsional excitations. No specific energy pathway is revealed that would funnel the hydrogen atom directly towards the trans-side, hence the experimentally observed high cis → trans quantum yield of close to one probably has to be explained in a statistical scenario on a timescale much longer than that of the present simulation.     

Incorporation of quantum effects for selected degrees of freedom into the trajectory-based dynamics using spatial domains

The approach of defining quantum corrections on nuclear dynamics of molecular systems incorporated approximately into selected degrees of freedom, is described. The approach is based on the Madelung-de-Broglie-Bohm formulation of time-dependent quantum mechanics which represents a wavefunction in terms of an ensemble of trajectories. The trajectories follow classical laws of motion except that the quantum potential, dependent on the wavefunction amplitude and its derivatives, is added to the external, classical potential. In this framework the quantum potential, determined approximately for practical reasons, is included only into the "quantum" degrees of freedom describing light particles such as protons, while neglecting with the quantum force for the heavy, nearly classical nuclei. The entire system comprised of light and heavy particles is described by a single wavefunction of full dimensionality. The coordinate space of heavy particles is divided into spatial domains or subspaces. The quantum force acting on the light particles is determined for each domain of similar configurations of the heavy nuclei. This approach effectively introduces parametric dependence of the reduced dimensionality quantum force, on classical degrees of freedom. This strategy improves accuracy of the quantum force and does not restrict interaction between the domains. The concept is illustrated for two-dimensional scattering systems, where the quantum force is required to reproduce vibrational energy of the quantum degree of freedom.     

Quantum simulation of solution phase intramolecular electron transfer rates in betaine-30

Mixed quantum-classical atomistic simulations have been carried out to investigate the mechanistic details of excited state intramolecular electron transfer in a betaine-30 molecule in acetonitrile. The key electronic degrees of freedom of the solute molecule are treated quantum mechanically using the semiempirical Pariser-Parr-Pople Hamiltonian, including the solvent influence on electronic structure. The intramolecular vibrational modes are also treated explicitly at a quantum level, with the remaining elements treated classically using empirical potentials. The electron-transfer rate, corresponding to S1 --> S0 relaxation, is evaluated via time-dependent perturbation theory with the explicit inclusion of the dynamics of solvation and intramolecular conformation. The calculations reveal that, while solvation dynamics is critical to the rate, the intramolecular torsional dynamics also plays an important role. The importance of the use of multiple high-frequency quantum modes is also discussed.     

Mixed time slicing in path integral simulations

A simple and efficient scheme is presented for using different time slices for different degrees of freedom in path integral calculations. This method bridges the gap between full quantization and the standard mixed quantum-classical (MQC) scheme and, therefore, still provides quantum mechanical effects in the less-quantized variables. Underlying the algorithm is the notion that time slices (beads) may be "collapsed" in a manner that preserves quantization in the less quantum mechanical degrees of freedom. The method is shown to be analogous to multiple-time step integration techniques in classical molecular dynamics. The algorithm and its associated error are demonstrated on model systems containing coupled high- and low-frequency modes; results indicate that convergence of quantum mechanical observables can be achieved with disparate bead numbers in the different modes. Cost estimates indicate that this procedure, much like the MQC method, is most efficient for only a relatively few quantum mechanical degrees of freedom, such as proton transfer. In this regime, however, the cost of a fully quantum mechanical simulation is determined by the quantization of the least quantum mechanical degrees of freedom.     

On a semiclassical approach to energy transfer in polyatomic molecules

A semiclassical method for energy transfer in polyatomic molecules is suggested. The method is based on a partial quantization of the vibrational degrees of freedom. The remaining degrees of freedom are treated classically. As an example a system with three quantum coordinates is considered. The Coriolis coupling terms are taken into account in the quantum mechanical part of the system and discussed for the N₂ + CO₂ system.     

Quantum dynamics through conical intersections in macrosystems: Combining effective modes and time-dependent Hartree

We address the nonadiabatic quantum dynamics of (macro)systems involving a vast number of nuclear degrees of freedom (modes) in the presence of conical intersections. The macrosystem is first decomposed into a system part carrying a few, strongly coupled modes, and an environment, comprising the remaining modes. By successively transforming the modes of the environment, a hierarchy of effective Hamiltonians for the environment can be constructed. Each effective Hamiltonian depends on a reduced number of effective modes, which carry cumulative effects. The environment is described by a few effective modes augmented by a residual environment. In practice, the effective modes can be added to the system’s modes and the quantum dynamics of the entire macrosystem can be accurately calculated on a limited time-interval. For longer times, however, the residual environment plays a role. We investigate the possibility to treat fully quantum mechanically the system plus a few effective environmental modes, augmented by the dynamics of the residual environment treated by the time-dependent Hartree (TDH) approximation. While the TDH approximation is known to fail to correctly reproduce the dynamics in the presence of conical intersections, it is shown that its use on top of the effective-mode formalism leads to much better results. Two numerical examples are presented and discussed; one of them is known to be a critical case for the TDH approximation.     

Quantum wave packet ab initio molecular dynamics: an approach to study quantum dynamics in large systems

A methodology to efficiently conduct simultaneous dynamics of electrons and nuclei is presented. The approach involves quantum wave packet dynamics using an accurate banded, sparse and Toeplitz representation for the discrete free propagator, in conjunction with ab initio molecular dynamics treatment of the electronic and classical nuclear degree of freedom. The latter may be achieved either by using atom-centered density-matrix propagation or by using Born-Oppenheimer dynamics. The two components of the methodology, namely, quantum dynamics and ab initio molecular dynamics, are harnessed together using a time-dependent self-consistent field-like coupling procedure. The quantum wave packet dynamics is made computationally robust by using adaptive grids to achieve optimized sampling. One notable feature of the approach is that important quantum dynamical effects including zero-point effects, tunneling, as well as over-barrier reflections are treated accurately. The electronic degrees of freedom are simultaneously handled at accurate levels of density functional theory, including hybrid or gradient corrected approximations. Benchmark calculations are provided for proton transfer systems and the dynamics results are compared with exact calculations to determine the accuracy of the approach.     

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.     

Quantum wavepacket ab initio molecular dynamics: generalizations using an extended Lagrangian treatment of diabatic states coupled through multireference electronic structure

We present a generalization to our previously developed quantum wavepacket ab initio molecular dynamics (QWAIMD) method by using multiple diabatic electronic reduced single particle density matrices, propagated within an extended Lagrangian paradigm. The Slater determinantal wavefunctions associated with the density matrices utilized may be orthogonal or nonorthogonal with respect to each other. This generalization directly results from an analysis of the variance in electronic structure with quantum nuclear degrees of freedom. The diabatic electronic states are treated here as classical parametric variables and propagated simultaneously along with the quantum wavepacket and classical nuclei. Each electronic density matrix is constrained to be N-representable. Consequently two sets of new methods are derived: extended Lagrangian-QWAIMD (xLag-QWAIMD) and diabatic extended Lagrangian-QWAIMD (DxLag-QWAIMD). In both cases, the instantaneous potential energy surface for the quantum nuclear degrees of freedom is constructed from the diabatic states using an on-the-fly nonorthogonal multireference formalism. By introducing generalized grid-based electronic basis functions, we eliminate the basis set dependence on the quantum nucleus. Subsequent reuse of the two-electron integrals during the on-the-fly potential energy surface computation stage yields a substantial reduction in computational costs. Specifically, both xLag-QWAIMD and DxLag-QWAIMD turn out to be about two orders of magnitude faster than our previously developed time-dependent deterministic sampling implementation of QWAIMD. Energy conservation properties, accuracy of the associated potential surfaces, and vibrational properties are analyzed for a family of hydrogen bonded systems.