Shu Yinan, Truhlar Donald G
Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States.
J Chem Theory Comput. 2024 Jun 11;20(11):4396-4426. doi: 10.1021/acs.jctc.4c00424. Epub 2024 May 31.
We reconsider recent methods by which direct dynamics calculations of electronically nonadiabatic processes can be carried out while requiring only adiabatic potential energies and their gradients. We show that these methods can be understood in terms of a new generalization of the well-known semiclassical Ehrenfest method. This is convenient because it eliminates the need to evaluate electronic wave functions and their matrix elements along the mixed quantum-classical trajectories. The new approximations and procedures enabling this advance are the curvature-driven approximation to the time-derivative coupling, the generalized semiclassical Ehrenfest method, and a new gradient correction scheme called the time-derivative matrix (TDM) scheme. When spin-orbit coupling is present, one can carry out dynamics calculations in the fully adiabatic basis using potential energies and gradients calculated without spin-orbit coupling plus the spin-orbit coupling matrix elements. Even when spin-orbit coupling is neglected, the method is useful because it allows calculations by electronic structure methods for which nonadiabatic coupling vectors are unavailable. In order to place the new considerations in context, the article starts out with a review of background material on trajectory surface hopping, the semiclassical Ehrenfest scheme, and methods for incorporating decoherence. We consider both internal conversion and intersystem crossing. We also review several examples from our group of successful applications of the curvature-driven approximation.
我们重新审视了近期的一些方法,通过这些方法可以进行电子非绝热过程的直接动力学计算,同时只需要绝热势能及其梯度。我们表明,这些方法可以从著名的半经典埃伦费斯特方法的一种新推广的角度来理解。这很方便,因为它消除了在混合量子 - 经典轨迹上评估电子波函数及其矩阵元的需要。实现这一进展的新近似方法和程序包括对时间导数耦合的曲率驱动近似、广义半经典埃伦费斯特方法以及一种称为时间导数矩阵(TDM)方案的新梯度校正方案。当存在自旋 - 轨道耦合时,可以在完全绝热基下进行动力学计算,使用在不考虑自旋 - 轨道耦合的情况下计算的势能和梯度加上自旋 - 轨道耦合矩阵元。即使忽略自旋 - 轨道耦合,该方法也很有用,因为它允许使用无法获得非绝热耦合矢量的电子结构方法进行计算。为了将新的考虑因素置于背景中,本文首先回顾了关于轨迹表面跳跃、半经典埃伦费斯特方案以及纳入退相干的方法的背景材料。我们考虑了内转换和系间窜越。我们还回顾了我们小组中曲率驱动近似成功应用的几个例子。