Center for Superfunctional Materials, Department of Chemistry, Pohang University of Science and Technology, Namgu, Pohang, Korea.
J Chem Phys. 2011 Dec 28;135(24):244112. doi: 10.1063/1.3671952.
The electron dynamics with complex third-order Suzuki-Trotter propagator (ST(3)) has been implemented into a planewave (PW) based density functional theory program, and several applications including linear absorption spectra and coupled electron-nucleus dynamics have been calculated. Since the ST(3) reduces the number of Fourier transforms to less than half compared to the fourth-order Suzuki-Trotter propagator (ST(4)), more than twice faster calculations are possible by exploiting the ST(3). We analyzed numerical errors of both the ST(3) and the ST(4) in the presence∕absence of an external field for several molecules such as Al(2), N(2), and C(2)H(4). We obtained that the ST(3) gives the same order of numerical errors (10(-5) Ry after 100 fs) as the ST(4). Also, the time evolution of dipole moments, hence the absorption spectrum, is equivalent for both ST(3) and ST(4). As applications, the linear absorption spectrum for an ethylene molecule was studied. From the density difference analysis, we showed that the absorption peaks at 6.10 eV and 7.65 eV correspond to the π → 4a(g) and π → π* excitation bands, respectively. We also investigated the molecular vibrational effect to the absorption spectra of an ethylene molecule and the dynamics of a hydrogen molecule after the σ → σ* transition by formulating coupled electron-nucleus dynamics within the Ehrenfest regime. The trajectory of nuclei follows the excited state potential energy curve exactly.
已经将具有复杂三阶 Suzuki-Trotter 传播子(ST(3))的电子动力学实现到基于平面波(PW)的密度泛函理论程序中,并计算了包括线性吸收光谱和耦合电子-核动力学在内的几个应用。由于 ST(3)与第四阶 Suzuki-Trotter 传播子(ST(4))相比,将傅里叶变换的数量减少到不到一半,因此通过利用 ST(3)可以实现两倍以上的快速计算。我们分析了在存在/不存在外部场的情况下,ST(3)和 ST(4)的数值误差对于 Al(2)、N(2)和 C(2)H(4)等几种分子。我们得到 ST(3)给出的数值误差相同(10(-5) Ry 在 100 fs 之后)作为 ST(4)。此外,偶极矩的时间演化,因此吸收光谱,对于 ST(3)和 ST(4)都是等效的。作为应用,研究了乙烯分子的线性吸收光谱。从密度差分析,我们表明 6.10 eV 和 7.65 eV 处的吸收峰分别对应于π→4a(g)和π→π激发带。我们还通过在 Ehrenfest 范围内制定耦合电子-核动力学,研究了乙烯分子吸收光谱和 σ→σ跃迁后氢分子的分子振动效应。原子核的轨迹完全遵循激发态势能曲线。