Peng Daoling, Zou Wenli, Liu Wenjian
Institute of Theoretical and Computational Chemistry, State Key Laboratory of Rare Earth Materials Chemistry and Applications, Peking University, Beijing, People's Republic of China.
J Chem Phys. 2005 Oct 8;123(14):144101. doi: 10.1063/1.2047554.
A time-dependent quasirelativistic density-functional theory for excitation energies of systems containing heavy elements is developed, which is based on the zeroth-order regular approximation (ZORA) for the relativistic Hamiltonian and a noncollinear form for the adiabatic exchange-correlation kernel. To avoid the gauge dependence of the ZORA Hamiltonian a model atomic potential, instead of the full molecular potential, is used to construct the ZORA kinetic operator in ground-state calculations. As such, the ZORA kinetic operator no longer responds to changes in the density in response calculations. In addition, it is shown that, for closed-shell ground states, time-reversal symmetry can be employed to simplify the eigenvalue equation into an approximate form that is similar to that of time-dependent nonrelativistic density-functional theory. This is achieved by invoking an independent-particle approximation for the induced density matrix. The resulting theory is applied to investigate the global potential-energy curves of low-lying LambdaS- and omega omega-coupled electronic states of the AuH molecule. The derived spectroscopic parameters, including the adiabatic and vertical excitation energies, equilibrium bond lengths, harmonic and anharmonic vibrational constants, fundamental frequencies, and dissociation energies, are in good agreement with those of time-dependent four-component relativistic density-functional theory and ab initio multireference second-order perturbation theory. Nonetheless, this two-component relativistic version of time-dependent density-functional theory is only moderately advantageous over the four-component one as far as computational efforts are concerned.
我们发展了一种用于含重元素体系激发能的含时准相对论密度泛函理论,它基于相对论哈密顿量的零阶正则近似(ZORA)以及绝热交换关联核的非共线形式。为避免ZORA哈密顿量的规范依赖性,在基态计算中使用模型原子势而非全分子势来构建ZORA动能算符。这样一来,在响应计算中ZORA动能算符不再对密度变化作出响应。此外,研究表明,对于闭壳层基态,可利用时间反演对称性将本征值方程简化为一种近似形式,该形式类似于含时非相对论密度泛函理论的形式。这是通过对诱导密度矩阵采用独立粒子近似来实现的。所得理论被用于研究AuH分子低激发态LambdaS和ωω耦合电子态的全局势能曲线。推导得到的光谱参数,包括绝热和垂直激发能、平衡键长、简谐和非简谐振动常数、基频以及解离能,与含时四分量相对论密度泛函理论和从头算多参考二阶微扰理论的结果吻合良好。然而,就计算量而言,这种含时密度泛函理论的两分量相对论版本相较于四分量版本仅具有一定程度的优势。