Virtual Orbital Many-Body Expansions: A Possible Route towards the Full Configuration Interaction Limit
In the present letter, it is demonstrated how full configuration interaction (FCI) results in extended basis sets may be obtained to within sub-kJ/mol accuracy by decomposing the energy in terms of many-body expansions in the virtual orbitals of the molecular system at hand. This extension of the FCI application range lends itself to two unique features of the current approach, namely that the total energy calculation can be performed entirely within considerably reduced orbital subspaces and may be so by means of embarrassingly parallel programming. Facilitated by a rigorous and methodical screening protocol and further aided by expansion points different from the Hartree-Fock solution, al…
Communication: The performance of non-iterative coupled cluster quadruples models
We compare the numerical performance of various non-iterative coupled cluster (CC) quadruples models. The results collectively show how approaches that attempt to correct the CC singles and doubles energy for the combined effect of triple and quadruple excitations all fail at recovering the correlation energy of the full CC singles, doubles, triples, and quadruples (CCSDTQ) model to within sufficient accuracy. Such a level of accuracy is only achieved by models that make corrections to the full CC singles, doubles, and triples (CCSDT) energy for the isolated effect of quadruple excitations of which the CCSDT(Q–3) and CCSDT(Q–4) models of the Lagrangian-based CCSDT(Q–n) perturbation series a…
On the convergence of perturbative coupled cluster triples expansions:Error cancellations in the CCSD(T) model and the importance of amplitude relaxation
Recently, we proposed a novel Lagrangian-based perturbation series-the CCSD(T-n) series-which systematically corrects the coupled cluster singles and doubles (CCSD) energy in orders of the Møller-Plesset fluctuation potential for effects due to triple excitations. In the present study, we report numerical results for the CCSD(T-n) series up through fourth order which show the predicted convergence trend throughout the series towards the energy of its target, the coupled cluster singles, doubles, and triples (CCSDT) model. Since effects due to the relaxation of the CCSD singles and doubles amplitudes enter the CCSD(T-n) series at fourth order (the CCSD(T-4) model), we are able to separate th…
Incremental Treatments of the Full Configuration Interaction Problem
The recent many-body expanded full configuration interaction (MBE-FCI) method is reviewed by critically assessing its advantages and drawbacks in the context of contemporary near-exact electronic structure theory. Besides providing a succinct summary of the history of MBE-FCI to date within a generalized and unified theoretical setting, its finer algorithmic details are discussed alongside our optimized computational implementation of the theory. A selected few of the most recent applications of MBE-FCI are revisited, before we close by outlining its future research directions as well as its place among modern near-exact wave function-based methods.
A Lagrangian framework for deriving triples and quadruples corrections to the CCSD energy.
Using the coupled cluster Lagrangian technique, we have determined perturbative corrections to the coupled cluster singles and doubles (CCSD) energy that converge towards the coupled cluster singles, doubles, and triples (CCSDT) and coupled cluster singles, doubles, triples, and quadruples (CCSDTQ) energies, considering the CCSD state as the unperturbed reference state and the fluctua- tion potential as the perturbation. Since the Lagrangian technique is utilized, the energy corrections satisfy Wigner’s 2n + 1 rule for the cluster amplitudes and the 2n + 2 rule for the Lagrange multi- pliers. The energy corrections define the CCSD perturbation series, CCSD(T–n) and CCSD(TQ–n), which are ter…
Generalized Many-Body Expanded Full Configuration Interaction Theory
Facilitated by a rigorous partitioning of a molecular system's orbital basis into two fundamental subspaces - a reference and an expansion space, both with orbitals of unspecified occupancy - we generalize our recently introduced many-body expanded full configuration interaction (MBE-FCI) method to allow for electron-rich model and molecular systems dominated by both weak and strong correlation to be addressed. By employing minimal or even empty reference spaces, we show through calculations on the one-dimensional Hubbard model with up to 46 lattice sites, the chromium dimer, and the benzene molecule how near-exact results may be obtained in a entirely unbiased manner for chemical and physi…
Equation-of-motion coupled cluster perturbation theory revisited
The equation-of-motion coupled cluster (EOM-CC) framework has been used for deriving a novel series of perturbative corrections to the coupled cluster singles and doubles energy that formally con- verges towards the full configuration interaction energy limit. The series is based on a Møller-Plesset partitioning of the Hamiltonian and thus size extensive at any order in the perturbation, thereby rem- edying the major deficiency inherent to previous perturbation series based on the EOM-CC ansatz. © 2014 AIP Publishing LLC. [http://dx.doi.org/10.1063/1.4873138]
Foreword: Prof. Gauss Festschrift
As guest editors, we are excited to present the Molecular Physics Festschrift in honour of Jurgen Gauss, professor of theoretical chemistry at the Johannes Gutenberg-Universitat Mainz, Germany, on ...
Assessment of the accuracy of coupled cluster perturbation theory for open-shell systems. I. Triples expansions
The accuracy at which total energies of open-shell atoms and organic radicals may be calculated is assessed for selected coupled cluster perturbative triples expansions, all of which augment the coupled cluster singles and doubles (CCSD) energy by a non-iterative correction for the effect of triple excitations. Namely, the second- through sixth-order models of the recently proposed CCSD(T-n) triples series [J. Chem. Phys. 140, 064108 (2014)] are compared to the acclaimed CCSD(T) model for both unrestricted as well as restricted open-shell Hartree-Fock (UHF/ROHF) reference determinants. By comparing UHF- and ROHF-based statistical results for a test set of 18 modest-sized open-shell species …
The Dalton quantum chemistry program system
Dalton is a powerful general-purpose program system for the study of molecular electronic structure at the Hartree-Fock, Kohn-Sham, multiconfigurational self-consistent-field, MOller-Plesset, confi ...
Convergence of coupled cluster perturbation theory.
The convergence of a recently proposed coupled cluster (CC) family of perturbation series [Eriksen, J. J. et al., J. Chem. Phys. 140, 064108 (2014)], in which the energetic difference between two CC models - a low-level parent and a high-level target model - is expanded in orders of the M{\o}ller-Plesset (MP) fluctuation potential, is investigated for four prototypical closed-shell systems (Ne, singlet methylene, distorted HF, and the fluoride anion) in standard and augmented basis sets. In these investigations, energy corrections of the various series have been calculated to high orders and their convergence radii determined by probing for possible front- and back-door intruder states, the…
Ground and Excited State First-Order Properties in Many-Body Expanded Full Configuration Interaction Theory
The recently proposed many-body expanded full configuration interaction (MBE-FCI) method is extended to excited states and static first-order properties different from total, ground state correlation energies. Results are presented for excitation energies and (transition) dipole moments of two prototypical, heteronuclear diatomics---LiH and MgO---in augmented correlation consistent basis sets of up to quadruple-$\zeta$ quality. Given that MBE-FCI properties are evaluated without recourse to a sampled wave function and the storage of corresponding reduced density matrices, the memory overhead associated with the calculation of general first-order properties only scales with the dimension of …
Assessment of the accuracy of coupled cluster perturbation theory for open-shell systems. II. Quadruples expansions
We extend our assessment of the potential of perturbative coupled cluster (CC) expansions for a test set of open-shell atoms and organic radicals to the description of quadruple excitations. Namely, the second- through sixth-order models of the recently proposed CCSDT(Q-n) quadruples series [J. Chem. Phys. 140, 064108 (2014)] are compared to the prominent CCSDT(Q) and lambda-CCSDT(Q) models. From a comparison of the models in terms of their recovery of total CC singles, doubles, triples, and quadruples (CCSDTQ) energies, we find that the performance of the CCSDT(Q-n) models is independent of the reference used (unrestricted or restricted (open-shell) Hartree-Fock), in contrast to the CCSDT(…
Efficient and portable acceleration of quantum chemical many-body methods in mixed floating point precision using OpenACC compiler directives
It is demonstrated how the non-proprietary OpenACC standard of compiler directives may be used to compactly and efficiently accelerate the rate-determining steps of two of the most routinely applied many-body methods of electronic structure theory, namely the second-order M{\o}ller-Plesset (MP2) model in its resolution-of-the-identity (RI) approximated form and the (T) triples correction to the coupled cluster singles and doubles model (CCSD(T)). By means of compute directives as well as the use of optimized device math libraries, the operations involved in the energy kernels have been ported to graphics processing unit (GPU) accelerators, and the associated data transfers correspondingly o…
The Ground State Electronic Energy of Benzene.
We report on the findings of a blind challenge devoted to determining the frozen-core, full configuration interaction (FCI) ground state energy of the benzene molecule in a standard correlation-consistent basis set of double-$\zeta$ quality. As a broad international endeavour, our suite of wave function-based correlation methods collectively represents a diverse view of the high-accuracy repertoire offered by modern electronic structure theory. In our assessment, the evaluated high-level methods are all found to qualitatively agree on a final correlation energy, with most methods yielding an estimate of the FCI value around $-863$ m$E_{\text{H}}$. However, we find the root-mean-square devia…
Many-Body Expanded Full Configuration Interaction. II. Strongly Correlated Regime
In this second part of our series on the recently proposed many-body expanded full configuration interaction (MBE-FCI) method, we introduce the concept of multideterminantal expansion references. Through theoretical arguments and numerical validations, the use of this class of starting points is shown to result in a focussed compression of the MBE decomposition of the FCI energy, thus allowing chemical problems dominated by strong correlation to be addressed by the method. The general applicability and performance enhancements of MBE-FCI are verified for standard stress tests such as the bond dissociations in H$_2$O, N$_2$, C$_2$, and a linear H$_{10}$ chain. Furthermore, the benefits of em…