高精度相对论密度泛函计算方法

对适用于含重元素体系的高精度相对论密度泛函计算方法作简要的评述.结合本实验室的研究工作,重点介绍严格处理相对论效应的四分量相对论密度泛函计算方法和近似处理相对论效应的两分量和标量相对论密度泛函计算方法,包括零级规则展开近似(ZORA)方法及其改进和排除奇点的近似展开(SEAX)方法,以及适合处理局部包含重元素大体系的接合两分量-标量相对论(或非相对论)计算方法.     

对含重元素体系的接合二分量-标量相对论密度泛函计算方法

在相对论密度泛函ZORA方法的基础上,提出一种用于含重元素体系的接合二分量-标量相对论密度泛函计算方法.对于只含少数几个重元素的较大体系,仅对其中旋轨耦合作用强的重元素作二分量相对论计算,而对体系的其余部分则作标量相对论计算,通过对动能矩阵元的近似处理实现两种计算的接合.对一系列含6p区重元素分子进行计算的结果表明,当非重元素是第三周期以前的元素时,此方法与二分量ZO-RA方法的计算结果吻合得很好.当非重元素为第四周期元素时,计算结果有一定偏差,表明在后一种情况下旋轨耦合作用已比较显著,但误差仍在目前近似密度泛函计算的精度范围内.此方法可以有效地节省计算量,而且避免了Dyall方法的缺点.     

几种密度泛函理论公式用于镧系硫属化合物计算的比较

以镧系元素La、Gd、Lu的硫属化合物为对象,系统考察几种密度泛函理论公式对镧系化合物计算的适用情况,考虑了相对论效应的影响. 计算结果显示,相对论效应引起的键长变化在十2Pm到-3Pm之间,引起的键能减小为0. 4～0. 6eV,与采用的密度泛函公式关系不明显. 不同的密度泛函公式对键长的计算结果影响也不太大,但对键能有显著影响,其中LDA(VWN)＋PW86X公式给出最好结果. 交换能梯度校正明显改善键能计算结果,而相关能梯度校正反而使之变差. 简单的X_α公式给出相当好的键能计算结果. 在考虑相对论效应和梯度校正以后,密度泛函理论方法给出比较可靠的键长数值,键能则仍然偏高,但不超过20％.     

元素电负性和硬度的密度泛函理论研究

应用密度泛函理论的DFT LDA、DFT LDA/NL和改进的Slater过渡态方法,把元素的电离能和电子亲合能的计算扩展到周期表的103种元素.并用有限差分方法计算了这103种元素的电负性和硬度.计算中考虑了相对论效应.计算结果比以前Robles等用密度泛函理论的XGL和Xα近似的交换相关泛函的计算结果有所改进,更接近实验值.     

密度泛函理论及其数值方法新进展

综述了密度泛函理论及其数值方法的最新进展.密度泛函理论的发展以寻找合适 的交换相关近似为主线,从最初的局域密度近似、广义梯度近似到现在的非局域泛函、自相 互作用修正,多种泛函形式的相继出现使得密度泛函理论可以提供越来越精确的计算结果. 除了交换相关近似的发展,近年来密度泛函理论向含时理论、相对论等方面的扩展也很活跃 .另外,在密度泛函理论体系发展的同时,相应的数值计算方法的发展也非常迅速.从古老 的有限差分、有限元到新兴的小波分析都被用来实现密度泛函理论的数值计算.与此同时, 线性标度的密度泛函理论算法日趋成熟,使得通过密度泛函理论研究诸如生物大分子之类的 体系成为可能.随着密度泛函理论本身及其数值方法的发展,它的应用也越来越广泛,一些新的应用领域和研究方向不断涌现.     

Yb、YbO电子激发态的相对论含时密度泛函理论研究

本文用基于精确二分量哈密顿（exact two—component Hamiltonian）的相对论含时密度泛函理论（time-dependent relativistic density functional theory）计算了Yb和YbO的电子激发态,并利用对称性、自然原子轨道对激发态性质和归属进行了详细分析,所得结果支持实验对YbO基态与激发态的指认．     

相对论量子化学新进展*

本文详细阐述了相对论量子化学的基本概念和原理,在此基础上评述了相对论量子化学领域的最新进展。指出,不靠数学技巧,而仅凭"用原子（分子片）合成分子"这一思想就可以大大简化分子的相对论计算,使四分量完全相对论和二分量准相对论方法在简洁性、计算精度、计算效率诸方面达到完全一致;作者发展的新一代准相对论方法XQR（exact matrix quasi-relativistic theory）不仅准确、简单,而且是联系相对论Dirac方程和非相对论Schrodinger方程的"无缝桥梁"。这是概念上的一大突破。可以说,化学（和普通物理）中的相对论问题已经得到解决。本文还展望了相对论量子化学未来的发展方向。     

密度泛函理论研究镥二聚体低电子能态的相对稳定性、结构以及成键性质

用密度泛函理论方法研究了镥二聚体（Lu2）低能量电子态的性质,计算了电子态相对能量、平衡键长、振动频率以及基态解离能,考察了密度泛函性质、相对论有效势种类以及Hartree—Fock交换作用大小对计算结果的影响．结果表明,无论采用何种密度泛函和相对论有效势,体系的基态都为三重态,与其他一些基于分子轨道理论的从头计算方法得到的结论是一致的．另外,与分子轨道从头计算结果以及实验结果比较发现,采用杂化密度泛函理论和Stuttgart小核有效势计算得到的结果总体吻合最好．最后,特别分析研究了B3LYP计算中Hartree—Fock交换作用大小对基态键长和基态解离能的影响,发现随着交换作用的增大,键长增长,解离能减小,这是由于5d轨道杂化导致的共价成键作用减弱造成的．     

镱硫属化合物的密度泛函理论研究

用密度泛函理论(DFT)研究镱硫属化合物的电子结构和性质,通过与实验比较考察了现有的几种近似密度泛函公式对镧系元素化合物的适用程度和相对论效应的影响.结果表明,用DFT计算的YbO键长对实验值的偏差约为0.002nm;但得到的键能即使在考虑梯度校正和相对论效应之后,仍比实验值高,在定域密度近似基础上引入交换梯度校正使键能计算值减小,其中PW86x使键能计算值减小稍多些,结果更接近实验值;相关梯度校正使键能计算值升高.相对论效应使键长缩短0.004～0.006nm,键能减小约0.5eV.计算结果的分析表明,Yb的5d轨道和配体的np轨道间形成σ键和π键.在所研究的分子体系中,配体原子从O到Te、Yb原子的5d轨道布后数依次减少,同键能减弱的顺序一致.相对论效应使键能减小的主要原因是在成键过程中发生了Yb的6s电子向5d轨道的转移,而相对论效应使该过程能量增加.偶极矩和电荷分布的计算表明,Yb-L键以共价性为主,相对论效应使共价性成份增加.     

惰性元素化合物相对论效应的研究 I. 相对论效应的计算

分子形成中的相对论效应定义。用相对论电荷迭代EHT程序ITEREX计算了惰性元素氟化物XeF2, XeF4, XeF6, KrF2和RnF2的相对论效应, 分别为-72.91,是-160.34,-281.82, -25.44和-220.7kJ·mol^-^1。惰性元素化合物的相对论效应都是负值, 使分子趋于稳定。还计算了IF和CsF的相对论效应, 分别是10.68和17.42kJ·mol6-^1, 均为正值, 使分子能量升高, 证实了相对论效应在惰性元素化合物中的重要作用。     

Analytical energy gradient evaluation in relativistic and nonrelativistic density functional calculations

The expressions of analytical energy gradients in density functional theory and their implementation in programs are reported. The evaluation of analytical energy gradients can be carried out in the fully 4-component relativistic, approximate relativistic, and nonrelativistic density functional calculations under local density approximation or general gradient approximation with or without frozen core approximation using different basis sets in our programs. The translational invariance condition and the fact that the one-center terms do not contribute to the energy gradients are utilized to improve the calculation accuracy and to reduce the computational effort. The calculated results of energy gradients and optimized geometry as well as atomization energies of some molecules by the analytical gradient method are in very good agreement with results obtained by the numerical derivative method.     

Relativistic effects on the electronic structure and bonding of [Ir(CN)5]3−

Four-component relativistic and nonrelativistic molecular orbital calculations were performed for the covalent paramagnetic complex [Ir(CN)₅]³⁻, employing the self-consistent discrete variational method, in the framework of density functional theory. Relativistic effects on the electronic structure and chemical bonding are discussed by comparison of relativistic and nonrelativistic one-electron energy levels, populations, and bond orders. The influence of relativistic effects on calculated absorption energies of the electronic spectrum is briefly assessed. © 1996 John Wiley & Sons, Inc.     

Four-component relativistic Kohn-Sham theory

A four-component relativistic implementation of Kohn-Sham theory for molecular systems is presented. The implementation is based on a nonredundant exponential parametrization of the Kohn-Sham energy, well suited to studies of molecular static and dynamic properties as well as of total electronic energies. Calculations are presented of the bond lengths and the harmonic and anharmonic vibrational frequencies of Au(2), Hg(2+)(2), HgAu(+), HgPt, and AuH. All calculations are based on the full four-component Dirac-Coulomb Hamiltonian, employing nonrelativistic local, gradient-corrected, and hybrid density functionals. The relevance of the Coulomb and Breit operators for the construction of relativistic functionals is discussed; it is argued that, at the relativistic level of density-functional theory and in the absence of a vector potential, the neglect of current functionals follows from the neglect of the Breit operator.     

Ab initio and density functional theory study of structures and energies for dimethoxymethane as a model for the anomeric effect

Ab initio molecular orbital theory and density functional theory calculations have been carried out on dimethoxymethane as a model for the anomeric effect. We optimized various conformations of dimethoxymethane using Gaussian 92 at the MP2/6-311 + + G**, MP2/DZP + Diffuse, MP2/6-31G**, and Becke3LYP/6-31G** levels of theory. These methods were evaluated based on their performance in reproducing structures and energies of dimethoxymethane when compared to experiment. This study also examined the structure and energy of dimethoxymethane as a function of dihedral angles for examining the anomeric effect at the MP2/6-31G** and Becke3LYP/6-31G** levels of theory. These calculations are qualitatively consistent with the anomeric effect observations in carbohydrates and with earlier calculations. Quantitative comparisons with earlier results reveal that dimethoxymethane has lower total energies, smaller rotational barriers, and shorter bond lengths than was previously determined. The Becke3LYP calculations were also compared to the MP2 results. The density functional theory findings show that the minimum energy structures correspond well with experimental and MP2 data. The total and relative energies from molecular orbital theory and density functional theory vary to some extent. Contour plots of the relative energies of dimethoxymethane were evaluated and compared to a relative energy contour plot determined by MM3. The contour plots were similar, showing slightly larger changes in energies for the MP2 results than for the Becke3LYP results, which in turn were slightly larger than the MM3 results. Density functional theory calculations are an excellent alternative method of calculation due to increased speed and reliable accuracy of the density functional calculations. These results will serve as a benchmark for modelling the anomeric effect in carbohydrates. © 1996 John Wiley & Sons, Inc.     

Time-dependent four-component relativistic density functional theory for excitation energies

Time-dependent four-component relativistic density functional theory within the linear response regime is developed for calculating excitation energies of heavy element containing systems. Since spin is no longer a good quantum number in this context, we resort to time-reversal adapted Kramers basis when deriving the coupled Dirac-Kohn-Sham equation. The particular implementation of the formalism into the Beijing density functional program package utilizes the multipolar expansion of the induced density to facilitate the construction of the induced Coulomb potential. As the first application, pilot calculations on the valence excitation energies and fine structures of the rare gas (Ne to Rn) and Group 12 (Zn to Hg) atoms are reported. To the best of our knowledge, it is the first time to be able to account for spin-orbit coupling within time-dependent density functional theory for excitation energies.     

Implementation of the SCC-DFTB method for hybrid QM/MM simulations within the amber molecular dynamics package

Self-consistent charge density functional tight-binding (SCC-DFTB) is a semiempirical method based on density functional theory and has in many cases been shown to provide relative energies and geometries comparable in accuracy to full DFT or ab initio MP2 calculations using large basis sets. This article shows an implementation of the SCC-DFTB method as part of the new QM/MM support in the AMBER 9 molecular dynamics program suite. Details of the implementation and examples of applications are shown.     

The electronic structure of alkali aurides. A four-component Dirac-Kohn-Sham study

Spectroscopic constants, including dissociation energies, harmonic and anharmonic vibrational frequencies, and dipole moments, are calculated for the complete alkali auride series (LiAu, NaAu, KAu, RbAu, CsAu). The four-component formulation of relativistic density functional theory has been employed in this study, using the G-spinor basis sets implemented recently in the program BERTHA. The performance of four standard nonrelativistic density functionals employed is investigated by comparing the results with the best available theoretical and experimental data. The present work provides the first theoretical predictions on the molecular properties of RbAu. The intermetallic bond that occurs in the alkali auride series is highly polar and is characterized by a large charge transfer from the alkali metals to gold. The extent of this electron transfer has been investigated using several different charge analysis methods, enabling us to reach some general conclusions on their relative performance. We further report a detailed analysis of the topological properties of relativistic electron density in the bonding region, discussing the features of this approach which characterize the nature of the chemical bond. We have also computed the fully relativistic density for the alkali halides MBr and MI (M = Li, Na, K, Rb, and Cs). The comparative study shows that, on the basis of several topological properties and the variation in bond lengths, the gold atom behaves similarly to a halogen intermediate between Br and I.     

Direct and indirect causes of Fermi level pinning at the SiO/GaAs interface

The correlation between atomic bonding sites and the electronic structure of SiO on GaAs(001)-c(2x8)/(2x4) was investigated using scanning tunneling microscopy (STM), scanning tunneling spectroscopy (STS), and density functional theory (DFT). At low coverage, STM images reveal that SiO molecules bond Si end down; this is consistent with Si being undercoordinated and O being fully coordinated in molecular SiO. At approximately 5% ML (monolayer) coverage, multiple bonding geometries were observed. To confirm the site assignments from STM images, DFT calculations were used to estimate the total adsorption energies of the different bonding geometries as a function of SiO coverage. STS measurements indicated that SiO pins the Fermi level midgap at approximately 5% ML coverage. DFT calculations reveal that the direct causes of Fermi level pinning at the SiO GaAs(001)-(2x4) interface are a result of either local charge buildups or the generation of partially filled dangling bonds on Si atoms.     

Accelerating Kohn‐Sham response theory using density fitting and the auxiliary‐density‐matrix method

An extension of the formulation of the atomic‐orbital‐based response theory of Larsen et al., JCP 113, 8909 (2000) is presented. This new framework has been implemented in LSDalton and allows for the use of Kohn‐Sham density‐functional theory with approximate treatment of the Coulomb and Exchange contributions to the response equations via the popular resolution‐of‐the‐identity approximation as well as the auxiliary‐density matrix method (ADMM). We present benchmark calculations of ground‐state energies as well as the linear and quadratic response properties: vertical excitation energies, polarizabilities, and hyperpolarizabilities. The quality of these approximations in a range of basis sets is assessed against reference calculations in a large aug‐pcseg‐4 basis. Our results confirm that density fitting of the Coulomb contribution can be used without hesitation for all the studied properties. The ADMM treatment of exchange is shown to yield high accuracy for ground‐state and excitation energies, whereas for polarizabilities and hyperpolarizabilities the performance gain comes at a cost of accuracy. Excitation energies of a tetrameric model consisting of units of the P700 special pigment of photosystem I have been studied to demonstrate the applicability of the code for a large system.