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

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

化学动力学中的非绝热过程及其理论研究

伴随着原子核的运动出现的电子态之间的跃迁,被称为非绝热过程.近年来,实验结果表明,非绝热过程广泛存在于多种化学体系中.从理论上理解这一过程极具挑战性,因为该过程中电子和原子核的运动耦合在一起,量子化学的基本假设"波恩-奥本海默近似"被打破.本文介绍通过理论化学方法模拟非绝热过程的现状,讨论所涉及的非绝热动力学(量子和半经典动力学)和量子化学计算方法(高精度电子相关方法、半经验方法和TDDFT方法),同时也介绍一些典型的非绝热过程(生物分子的光稳定性、光异构化、自旋翻转、激发态电子和能量转移).尽管目前非绝热动力学领域已经取得了很大进展,但处理复杂体系仍然面临诸多挑战.解决这一问题需要动力学理论、量子化学方法和计算机技术等领域的共同发展.     

轨线面跳跃方法模拟有机过氧化物化学发光原理及关键问题

化学发光反应是一类特殊的激发态化学反应,在化学分析和生活中具有广泛应用.本文介绍了轨线面跳跃非绝热动力学方法模拟化学发光反应的一般原理和关键的技术,包括初始采样、势能和力的计算方法以及含时和不含时策略中跃迁概率的计算以及动力学模拟过程.     

分子激发态非绝热动力学的理论模拟

<正>光诱导反应过程是今天物理化学研究的中心问题之一。光吸收导致分子跃迁到电子激发态,发生一系列复杂的光诱导反应,如荧光、磷光、电子和能量转移、光化学反应等。当分子体系处于激发态上时,伴随着原子核的运动,不同电子态之间可能会出现势能面交叉。在交叉区域体系会从一个电子态"跳跃"到另外一个电子态,发生非绝热跃迁~(1,2)。该过程广泛的存在于自然界中。从简单的     

防晒化合物咖啡酸的直接-非绝热动力学（英文）

作为防晒霜中一种广泛使用的成分,咖啡酸具有很强的紫外线吸收能力.本文通过研究咖啡酸的光诱导动力学,探讨其光保护机制.在半经验电子结构理论水平上使用直接面跳跃模拟方法研究最稳定的咖啡酸异构体非绝热动力学.从不同电子态开始进行动力学模拟,探讨光诱导反应通道对激发波长的依赖性.结果表明,氢原子解离和环变形驱动的通道主导了咖啡酸的非绝热动力学.同时改变激发波长会调控通道的分支比例,本研究揭示了咖啡酸化合物光保护原理的分子机制.     

化学反应的高精度从头算势能面

势能面是化学反应动力学研究的基础。近年来随着理论方法的发展与计算技术的进步,不但含三、四个原子反应体系的电子基态势能面的构建精度进一步提高,一些反应体系的多电子态耦合势能面的构建和含六个原子以上反应体系的高维从头算势能面的构建也取得了重要进展。本文结合若干典型体系势能面的构建工作,主要介绍了高精度电子基态势能面,包括Renner-Teller、旋轨耦合等非绝热效应的耦合势能面以及高维势能面方面的研究进展。     

H+BrCH_3体系反应动力学的准经典轨迹研究

本文在推广LEPS势能面上用准经典轨迹理论方法计算了20000条碰撞轨迹,得到各种散射通道的几率和速度常数,讨论了振动增强、通道竞争、振动绝热性、质量组合及角动量耦合等重要动力学效应,并建立了它们与势能面性质之间的联系。据此判断该反应按直接型回弹机理进行并可满意地解释交叉分子束实验观测到的不对称性产物角分布。     

改进的半经典动力学模拟二苯乙烯光致顺反异构化反应

发展了一种改进的半经典动力学模拟方法, 并将其程序化用于气相二苯乙烯光致顺反异构化反应的机理研究.新的方法不仅采用e 指数模型改进了原有Zhu-Nakamura 理论中计算电子非绝热跃迁几率的计算方法, 而且将约束哈密顿方法用于限制性分子动力学模拟过程中. 计算结果表明, 采用此方法得到的统计平均的量子产率及反应机理与以前的实验与理论结果吻合较好, 从而可以应用于全量子动力学方法无法进行的大分子体系的动力学研究.     

基本不变量神经网络解析梯度方法的研究

提升势能面的运行速度对于动力学模拟至关重要. 相对于计算简单、 但耗时更长的数值梯度计算, 直接求解势能面梯度的解析公式能够大幅提高势能面的运行效率. 本文发展了基本不变量神经网络解析梯度的生成方法. 计算解析梯度的代码可以通过程序自动生成. 对大量数据点进行测试后, 证明了该方法可以得到正确的势能面梯度输出结果. 通过测试不同势能面的调用时间, 发现采用解析梯度方法能够带来10倍以上的性能提升. 随着体系的增大, 这种性能提升也会越明显.     

基于新一代密度泛函和神经网络势能面的量子反应动力学计算

介绍了近年来发展起来的新一代密度泛函XYG3及利用神经网络构造分子体系势能面的最新进展。以H₃和CH₅等体系为实例,表明基于高效准确的密度泛函电子结构计算,与神经网络高精度势能面构造的理想结合,可以得到可靠的化学动力学结果,并有望用于更大更复杂的体系。     

Adiabatic states derived from a spin-coupled diabatic transformation: semiclassical trajectory study of photodissociation of HBr and the construction of potential curves for LiBr+

The development of spin-coupled diabatic representations for theoretical semiclassical treatments of photodissociation dynamics is an important practical goal, and some of the assumptions required to carry this out may be validated by applications to simple systems. With this objective, we report here a study of the photodissociation dynamics of the prototypical HBr system using semiclassical trajectory methods. The valence (spin-free) potential energy curves and the permanent and transition dipole moments were computed using high-level ab initio methods and were transformed to a spin-coupled diabatic representation. The spin-orbit coupling used in the transformation was taken as that of atomic bromine at all internuclear distances. Adiabatic potential energy curves, nonadiabatic couplings and transition dipole moments were then obtained from the diabatic ones and were used in all the dynamics calculations. Nonadiabatic photodissociation probabilities were computed using three semiclassical trajectory methods, namely, coherent switching with decay of mixing (CSDM), fewest switches with time uncertainty (FSTU), and its recently developed variant with stochastic decoherence (FTSU/SD), each combined with semiclassical sampling of the initial vibrational state. The calculated branching fraction to the higher fine-structure level of the bromine atom is in good agreement with experiment and with more complete theoretical treatments. The present study, by comparing our new calculations to wave packet calculations with distance-dependent ab initio spin-orbit coupling, validates the semiclassical trajectory methods, the semiclassical initial state sample scheme, and the use of a distance-independent spin-orbit coupling for future applications to polyatomic photodissociation. Finally, using LiBr(+) as a model system, it is shown that accurate spin-coupled potential curves can also be constructed for odd-electron systems using the same strategy as for HBr.     

Simulating quantum dynamics in classical environments

Methods for simulating the dynamics of composite systems, where part of the system is treated quantum mechanically and its environment is treated classically, are discussed. Such quantum–classical systems arise in many physical contexts where certain degrees of freedom have an essential quantum character while the other degrees of freedom to which they are coupled may be treated classically to a good approximation. The dynamics of these composite systems are governed by a quantum–classical Liouville equation for either the density matrix or the dynamical variables which are operators in the Hilbert space of the quantum subsystem and functions of the classical phase space variables of the classical environment. Solutions of the evolution equations may be formulated in terms of surface-hopping dynamics involving ensembles of trajectory segments interspersed with quantum transitions. The surface-hopping schemes incorporate quantum coherence and account for energy exchanges between the quantum and classical degrees of freedom. Various simulation algorithms are discussed and illustrated with calculations on simple spin-boson models but the methods described here are applicable to realistic many-body environments.     

Trajectory-based solution of the nonadiabatic quantum dynamics equations: an on-the-fly approach for molecular dynamics simulations

The non-relativistic quantum dynamics of nuclei and electrons is solved within the framework of quantum hydrodynamics using the adiabatic representation of the electronic states. An on-the-fly trajectory-based nonadiabatic molecular dynamics algorithm is derived, which is also able to capture nuclear quantum effects that are missing in the traditional trajectory surface hopping approach based on the independent trajectory approximation. The use of correlated trajectories produces quantum dynamics, which is in principle exact and computationally very efficient. The method is first tested on a series of model potentials and then applied to study the molecular collision of H with H(2) using on-the-fly TDDFT potential energy surfaces and nonadiabatic coupling vectors.     

Coherent switching with decay of mixing: an improved treatment of electronic coherence for non-Born-Oppenheimer trajectories

The self-consistent decay-of-mixing (SCDM) semiclassical trajectory method for electronically nonadiabatic dynamics is improved by modifying the switching probability that determines the instantaneous electronic state toward which the system decoheres. This method is called coherent switching with decay of mixing (CSDM), and it differs from the previously presented SCDM method in that the electronic amplitudes controlling the switching of the decoherent state are treated fully coherently in the electronic equations of motion for each complete passage through a strong interaction region. It is tested against accurate quantum mechanical calculations for 12 atom-diatom scattering test cases. Also tested are the SCDM method and the trajectory surface hopping method of Parlant and Gislason that requires coherent passages through each strong interaction region, and which we call the "exact complete passage" trajectory surface hopping (ECP-TSH) method. The results are compared with previously presented results for the fewest switches with time uncertainty and Tully's fewest switches (TFS) surface hopping methods and the semiclassical Ehrenfest method. We find that the CSDM method is the most accurate of the semiclassical trajectory methods tested. Including coherent passages improves the accuracy of the SCDM method (i.e., the CSDM method is more accurate than the SCDM method) but not of the trajectory surface hopping method (i.e., the ECP-TSH method is not more accurate on average than the TFS method).     

Mixed quantum/classical investigation of the photodissociation of NH3(A) and a practical method for maintaining zero-point energy in classical trajectories

The photodissociation dynamics of ammonia upon excitation of the out-of-plane bending mode (mode nu(2) with n(2)=0,[ellipsis (horizontal)],6 quanta of vibration) in the A electronic state is investigated by means of several mixed quantum/classical methods, and the calculated final-state properties are compared to experiments. Five mixed quantum/classical methods are tested: one mean-field approach (the coherent switching with decay of mixing method), two surface-hopping methods [the fewest switches with time uncertainty (FSTU) and FSTU with stochastic decay (FSTU/SD) methods], and two surface-hopping methods with zero-point energy (ZPE) maintenance [the FSTUSD+trajectory projection onto ZPE orbit (TRAPZ) and FSTUSD+minimal TRAPZ (mTRAPZ) methods]. We found a qualitative difference between final NH(2) internal energy distributions obtained for n(2)=0 and n(2)>1, as observed in experiments. Distributions obtained for n(2)=1 present an intermediate behavior between distributions obtained for smaller and larger n(2) values. The dynamics is found to be highly electronically nonadiabatic with all these methods. NH(2) internal energy distributions may have a negative energy tail when the ZPE is not maintained throughout the dynamics. The original TRAPZ method was designed to maintain ZPE in classical trajectories, but we find that it leads to unphysically high internal vibrational energies. The mTRAPZ method, which is new in this work and provides a general method for maintaining ZPE in either single-surface or multisurface trajectories, does not lead to unphysical results and is much less time consuming. The effect of maintaining ZPE in mixed quantum/classical dynamics is discussed in terms of agreement with experimental findings. The dynamics for n(2)=0 and n(2)=6 are also analyzed to reveal details not available from experiment, in particular, the time required for quenching of electronic excitation and the adiabatic energy gap and geometry at the time of quenching.     

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.     

Quantum-classical dynamics of wave fields

An approach to the quantum-classical mechanics of phase space dependent operators, which has been proposed recently, is remodeled as a formalism for wave fields. Such wave fields obey a system of coupled nonlinear equations that can be written by means of a suitable non-Hamiltonian bracket. As an example, the theory is applied to the relaxation dynamics of the spin-boson model. In the adiabatic limit, a good agreement with calculations performed by the operator approach is obtained. Moreover, the theory proposed in this paper can take nonadiabatic effects into account without resorting to surface-hopping approximations. Hence, the results obtained follow qualitatively those of previous surface-hopping calculations and increase by a factor of (at least) 2, the time length over which nonadiabatic dynamics can be propagated with small statistical errors. Moreover, it is worth to note that the dynamics of quantum-classical wave fields proposed here is a straightforward non-Hamiltonian generalization of the formalism for nonlinear quantum mechanics that Weinberg introduced recently.