Bao Junwei Lucas, Odoh Samuel O, Gagliardi Laura, Truhlar Donald G
Department of Chemistry, Chemical Theory Center, and Supercomputing Institute, University of Minnesota , 207 Pleasant Street SE, Minneapolis, Minnesota 55455-0431, United States.
J Chem Theory Comput. 2017 Feb 14;13(2):616-626. doi: 10.1021/acs.jctc.6b01102. Epub 2017 Jan 4.
We study the performance of multiconfiguration pair-density functional theory (MC-PDFT) and multireference perturbation theory for the computation of the bond dissociation energies in 12 transition-metal-containing diatomic molecules and three small transition-metal-containing polyatomic molecules and in two transition-metal dimers. The first step is a multiconfiguration self-consistent-field calculation, for which two choices must be made: (i) the active space and (ii) its partition into subspaces, if the generalized active space formulation is used. In the present work, the active space is chosen systematically by using three correlated-participating-orbitals (CPO) schemes, and the partition is chosen by using the separated-pair (SP) approximation. Our calculations show that MC-PDFT generally has similar accuracy to CASPT2, and the active-space dependence of MC-PDFT is not very great for transition-metal-ligand bond dissociation energies. We also find that the SP approximation works very well, and in particular SP with the fully translated BLYP functional SP-ftBLYP is more accurate than CASPT2. SP greatly reduces the number of configuration state functions relative to CASSCF. For the cases of FeO and NiO with extended-CPO active space, for which complete active space calculations are unaffordable, SP calculations are not only affordable but also of satisfactory accuracy. All of the MC-PDFT results are significantly better than the corresponding results with broken-symmetry spin-unrestricted Kohn-Sham density functional theory. Finally we test a perturbation theory method based on the SP reference and find that it performs slightly worse than CASPT2 calculations, and for most cases of the nominal-CPO active space, the approximate SP perturbation theory calculations are less accurate than the much less expensive SP-PDFT calculations.
我们研究了多组态对密度泛函理论(MC-PDFT)和多参考扰动理论在计算12个含过渡金属的双原子分子、3个含过渡金属的小分子多原子分子以及2个过渡金属二聚体的键解离能方面的性能。第一步是进行多组态自洽场计算,为此必须做出两个选择:(i)活性空间;(ii)如果使用广义活性空间公式,将其划分为子空间。在本工作中,通过使用三种相关参与轨道(CPO)方案系统地选择活性空间,并通过使用分离对(SP)近似来选择划分。我们的计算表明,MC-PDFT通常具有与CASPT2相似的精度,并且对于过渡金属-配体键解离能,MC-PDFT对活性空间的依赖性不是很大。我们还发现SP近似效果很好,特别是使用完全平移的BLYP泛函的SP-ftBLYP比CASPT2更准确。相对于CASSCF,SP大大减少了组态态函数的数量。对于具有扩展CPO活性空间的FeO和NiO的情况,由于完全活性空间计算难以承受,SP计算不仅可行,而且精度令人满意。所有MC-PDFT结果都明显优于具有破缺对称性的自旋非限制Kohn-Sham密度泛函理论的相应结果。最后,我们测试了一种基于SP参考的扰动理论方法,发现它的性能略逊于CASPT2计算,并且对于大多数标称CPO活性空间的情况,近似SP扰动理论计算的精度低于成本低得多的SP-PDFT计算。