Bytautas Laimutis, Scuseria Gustavo E, Ruedenberg Klaus
Department of Chemistry, Galveston College, 4015 Ave. Q, Galveston, Texas 77550, USA.
Department of Chemistry, Rice University, Houston, Texas 77005, USA.
J Chem Phys. 2015 Sep 7;143(9):094105. doi: 10.1063/1.4929904.
The present study further explores the concept of the seniority number (Ω) by examining different configuration interaction (CI) truncation strategies in generating compact wave functions in a systematic way. While the role of Ω in addressing static (strong) correlation problem has been addressed in numerous previous studies, the usefulness of seniority number in describing weak (dynamic) correlation has not been investigated in a systematic way. Thus, the overall objective in the present work is to investigate the role of Ω in addressing also dynamic electron correlation in addition to the static correlation. Two systematic CI truncation strategies are compared beyond minimal basis sets and full valence active spaces. One approach is based on the seniority number (defined as the total number of singly occupied orbitals in a determinant) and another is based on an excitation-level limitation. In addition, molecular orbitals are energy-optimized using multiconfigurational-self-consistent-field procedure for all these wave functions. The test cases include the symmetric dissociation of water (6-31G), N2 (6-31G), C2 (6-31G), and Be2 (cc-pVTZ). We find that the potential energy profile for H2O dissociation can be reasonably well described using only the Ω = 0 sector of the CI wave function. For the Be2 case, we show that the full CI potential energy curve (cc-pVTZ) is almost exactly reproduced using either Ω-based (including configurations having up to Ω = 2 in the virtual-orbital-space) or excitation-based (up to single-plus-double-substitutions) selection methods, both out of a full-valence-reference function. Finally, in dissociation cases of N2 and C2, we shall also consider novel hybrid wave functions obtained by a union of a set of CI configurations representing the full valence space and a set of CI configurations where seniority-number restriction is imposed for a complete set (full-valence-space and virtual) of correlated molecular orbitals, simultaneously. We discuss the usefulness of the seniority number concept in addressing both static and dynamic electron correlation problems along dissociation paths.
本研究通过系统地研究不同的组态相互作用(CI)截断策略来生成紧凑波函数,进一步探索了 seniority 数(Ω)的概念。虽然在许多先前的研究中已经探讨了 Ω 在解决静态(强)相关问题中的作用,但尚未系统地研究 seniority 数在描述弱(动态)相关方面的有用性。因此,本工作的总体目标是研究 Ω 在解决除静态相关之外的动态电子相关中的作用。除了最小基组和全价活性空间之外,还比较了两种系统的 CI 截断策略。一种方法基于 seniority 数(定义为行列式中单占据轨道的总数),另一种基于激发能级限制。此外,对所有这些波函数使用多组态自洽场程序对分子轨道进行能量优化。测试案例包括水(6 - 31G)、N₂(6 - 31G)、C₂(6 - 31G)和 Be₂(cc - pVTZ)的对称解离。我们发现,仅使用 CI 波函数的 Ω = 0 扇区就能相当好地描述 H₂O 解离的势能曲线。对于 Be₂ 的情况,我们表明,使用基于 Ω 的(包括在虚拟轨道空间中具有高达 Ω = 2 的组态)或基于激发的(高达单加双取代)选择方法,都能几乎精确地重现全 CI 势能曲线(cc - pVTZ),这两种方法均基于全价参考函数。最后,在 N₂ 和 C₂ 的解离案例中,我们还将考虑通过将一组表示全价空间的 CI 组态与一组对相关分子轨道的完整集合(全价空间和虚拟空间)施加 seniority 数限制的 CI 组态相结合而获得的新型混合波函数。我们讨论了 seniority 数概念在解决解离路径上的静态和动态电子相关问题中的有用性。