绝热表象下非绝热分子动力学方法  

Nonadiabatic molecular dynamics under adiabatic representation

作　　者：孙震 吕项 李盛 安忠[3] Sun Zhen;LüXiang;Li Sheng;An Zhong(Department of Physics,Zhejiang Normal University,Jinhua 321004,China;Zhejiang Institute of Photoelectronics,Jinhua 321004,China;College of Physics,Hebei Normal University,Shijiazhuang 050024,China)

机构地区：[1]浙江师范大学物理与电子信息工程学院,金华321004 [2]浙江光电子研究院,金华321004 [3]河北师范大学物理学院,石家庄050024

出　　处：《物理学报》2024年第14期16-23,共8页Acta Physica Sinica

基　　金：浙江省自然科学基金(批准号:LY19A040007)资助的课题.

摘　　要：基于SSH(Su-Schriffer-Heeger)哈密顿的非绝热分子动力学方法广泛应用于模拟有机共轭聚合物中光激发过程和极化子运动.目前,该方法中电子波函数的演化是在透热表象中进行的,本文对该方法进行扩展,让电子波函数的演化在绝热表象下进行,给出了详细的公式推导过程.分别利用新、旧方法模拟了一条共轭聚合物链中光激发动力学过程,两种方法得到的数值计算结果相符.新方法可以加深对非绝热分子动力学方法的理解,提供激发态弛豫过程中不同分子轨道之间非绝热耦合强度等重要信息.In this paper,we develop a nonadiabatic molecular dynamics method based on Su-Schriffer-Heeger(SSH)Hamiltonian,and this method is widely used to study the photoexcitation dynamics and polaron motion in conjugated polymers.However,in this method,the time-dependent Schrödinger equation has so far been solved in a diabatic representation,also known as site representation.In order to provide a deeper insight into the nonadiabatic molecular dynamics method,we solve the time-dependent Schrödinger equation in an adiabatic representation.The new method can directly provide the important information about the strength of nonadiabatic couplings between different molecular orbitals in the excited-state relaxation process,helping us to predict the electron and energy transfer within or between polymer chains.Solving the time-dependent Schrödinger equation in an adiabatic representation is much more complicated,it is mainly because we need to calculate the nonadiabatic couplings between different molecular orbitals.In this paper,the detailed formula derivation and actual calculation process of the nonadiabatic molecular dynamics method in an adiabatic representation are given.Using this new method,we simulate three photoexcitation processes in a conjugated polymer chain,HOMO→LUMO,HOMO–1→LUMO+1 and HOMO–2→LUMO+2.We analyze in detail the time evolutions of lattice configuration for these three photoexcitation processes,and compare these results with those obtained by diabatic representation(site representation)showing that the results obtained from these two representations are consistent with each other.

关 键 词：非绝热分子动力学 绝热表象 共轭聚合物 Su-Schriffer-Heeger 模型 

分 类 号：O561[理学—原子与分子物理]