Perturbation Expansion of Internally Contracted Coupled-Cluster Theory up to Third Order

Author(s)

Aoto, Yuri Alexandre ^{[1, 2]} ; Bargholz, Arne ^[1] ; Kats, Daniel ^{[3, 1]} ; Werner, Hans-Joachim ^[1] ; Koehn, Andreas ^[1]

Affiliation

^[1] Univ Stuttgart, Inst Theoret Chem, Pfaffenwaldring 55, D-70569 Stuttgart - Germany

^[2] Fed Univ ABC UFABC, Ctr Math Comp & Cognit, Ave Estados 5001, Santo Andre - Brazil

^[3] Max Planck Inst Solid State Res, Heisenbergstr 1, D-70569 Stuttgart - Germany

• Document type: Journal article
• Source: JOURNAL OF CHEMICAL THEORY AND COMPUTATION; v. 15, n. 4, p. 2291-2305, APR 2019.
• Web of Science Citations: 4

Abstract

The internally contracted multireference coupled-cluster (icMRCC) method is analyzed through third order in perturbation theory. Up to second order, the icMRCC perturbation expansion is equivalent to that of the standard Rayleigh- Schrodinger perturbation theory, which is based on a linear ansatz for the wave function, and the resulting theory is, depending on the employed zeroth-order Hamiltonian, equivalent to either second-order complete active space perturbation theory (CASPT2), N-electron valence perturbation theory (NEVPT2), or Fink's retention of the excitation degree perturbation theory (REPT2). At third order, the icMRCC perturbation expansion features additional terms in comparison to the Rayleigh- Schrodinger perturbation theory, but these are shown to be nearly negligibly small by both analytic arguments and numerical examples. Considering these systematic cancellations, however, may be important in future work on approximations to icMRCC theory. In addition, we provide an extensive set of tests of the second and third-order perturbation theories based on three different zeroth-order Hamiltonians, namely, the projected effective Fock operator as used for CASPT, the Dyall Hamiltonian as used for NEVPT, and the Fink Hamiltonian used for REPT. While the third-order variant of REPT often gives absolute energies that are rather close to values from higher level calculations, the results for relative energies and spectroscopic constants such as harmonic frequencies, give a less clear picture and a general conclusion about any best zeroth-order Hamiltonian does not emerge from our data. For small active spaces, REPT is rather prone to intruder state problems. (AU)

• FAPESP's process: 18/04617-6 - The differentiable manifolds of the electronic structure theory
• Grantee: Yuri Alexandre Aoto
• Support Opportunities: Scholarships in Brazil - Young Researchers
• FAPESP's process: 17/21199-0 - The differentiable manifolds of the electronic structure theory
• Grantee: Yuri Alexandre Aoto
• Support Opportunities: Research Grants - Young Investigators Grants