Center for Molecular Modeling, Ghent University, Technologiepark, 903, B-9052 Zwijnaarde, Belgium.
Phys Chem Chem Phys. 2012 Aug 14;14(30):10690-704. doi: 10.1039/c2cp41086a. Epub 2012 Jun 28.
In this paper an overview is presented of several approximations within Density Functional Theory (DFT) to calculate g-tensors in transition metal containing systems and a new accurate description of the spin-other-orbit contribution for high spin systems is suggested. Various implementations in a broad variety of software packages (ORCA, ADF, Gaussian, CP2K, GIPAW and BAND) are critically assessed on various aspects including (i) non-relativistic versus relativistic Hamiltonians, (ii) spin-orbit coupling contributions and (iii) the gauge. Particular attention is given to the level of accuracy that can be achieved for codes that allow g-tensor calculations under periodic boundary conditions, as these are ideally suited to efficiently describe extended condensed-phase systems containing transition metals. In periodic codes like CP2K and GIPAW, the g-tensor calculation schemes currently suffer from an incorrect treatment of the exchange spin-orbit interaction and a deficient description of the spin-other-orbit term. In this paper a protocol is proposed, making the predictions of the exchange part to the g-tensor shift more plausible. Focus is also put on the influence of the spin-other-orbit interaction which becomes of higher importance for high-spin systems. In a revisited derivation of the various terms arising from the two-electron spin-orbit and spin-other-orbit interaction (SOO), new insight has been obtained revealing amongst other issues new terms for the SOO contribution. The periodic CP2K code has been adapted in view of this new development. One of the objectives of this study is indeed a serious enhancement of the performance of periodic codes in predicting g-tensors in transition metal containing systems at the same level of accuracy as the most advanced but time consuming spin-orbit mean-field approach. The methods are first applied on rhodium carbide but afterwards extended to a broad test set of molecules containing transition metals from the fourth, fifth and sixth row of the periodic table. The set contains doublets as well as high-spin molecules.
本文综述了几种密度泛函理论(DFT)中的近似方法,用于计算含过渡金属体系中的 g 张量,并提出了一种新的精确描述高自旋体系中自旋-轨道贡献的方法。在各种软件包(ORCA、ADF、Gaussian、CP2K、GIPAW 和 BAND)中,各种实现方法在多个方面进行了批判性评估,包括(i)非相对论与相对论哈密顿量,(ii)自旋-轨道耦合贡献,以及(iii)规范。特别关注那些允许在周期性边界条件下计算 g 张量的代码所能达到的精度水平,因为这些代码非常适合有效地描述含有过渡金属的扩展凝聚相体系。在 CP2K 和 GIPAW 等周期性代码中,g 张量计算方案目前存在交换自旋-轨道相互作用处理不当和自旋-轨道以外项描述不足的问题。本文提出了一种方案,使得 g 张量中交换部分的预测更加合理。同时还关注了自旋-轨道以外相互作用的影响,它对于高自旋体系变得更为重要。在对来源于双电子自旋-轨道和自旋-轨道以外相互作用(SOO)的各种项的重新推导中,获得了新的见解,揭示了 SOO 贡献的新项等问题。针对这一新发展,对周期性 CP2K 代码进行了调整。本研究的目标之一是在预测含过渡金属体系中的 g 张量方面,显著提高周期性代码的性能,达到最先进但耗时的自旋-轨道平均场方法的精度水平。这些方法首先应用于碳化铑,但随后扩展到一个包含过渡金属的广泛测试集的分子,这些分子来自元素周期表的第四、五和六周期。该集合包含双态和高自旋分子。