Internally contracted multireference coupled-cluster theory in a multistate framework

Performance Tests of the Second-Order Approximate Internally Contracted Multireference Coupled-Cluster Singles and Doubles Method icMRCC2

Benchmark results are presented for the second-order approximation of the internally contracted multireference coupled-cluster method with single and double excitations, icMRCC2 [Köhn, Bargholz, J. Chem. Phys. 2019, 151, 041106], which was designed as a multireference analogue of the single-reference second-order approximate coupled-cluster method CC2 [Christiansen, Koch, Jørgensen, Chem. Phys. Lett. 1995, 243, 409-418]. Vertical excitation energies of various small to medium-sized organic molecules are investigated based on established test sets from the literature. Additionally, the spectroscopic constants of ground and excited states of diatomics and the geometric parameters of excited triatomic molecules were determined and compared to the experimental data. The results show that the method clearly extends the applicability of single-reference CC2, including doubly excited states, and also artifacts of CC2 like too low Rydberg excitations and too weak multiple bonds are eliminated. The method is computationally more demanding than standard multireference second-order perturbation theories but improves significantly in accuracy, as shown by the benchmark results. In addition, it is demonstrated that small active spaces are often sufficient to obtain accurate energies with icMRCC2. Example applications like the automerization of cyclobutadiene, the deactivation pathway of ethylene, and the excited states of an iron complex with a noninnocent nitrosyl ligand demonstrate the potential of icMRCC2 in cases with strong multireference character.

Assessment of State-Averaged Driven Similarity Renormalization Group on Vertical Excitation Energies: Optimal Flow Parameters and Applications to Nucleobases

We present a comprehensive excited-state benchmark for the state-averaged (SA) driven similarity renormalization group (DSRG) [Li, C.; Evangelista, F. A. J. Chem. Phys. 2018, 148, 124106]. Following the QUEST database [Véril, M.; Scemama, A.; Caffarel, M.; Lipparini, F.; Boggio-Pasqua, M.; Jacquemin, D.; Loos, P.-F. Wiley Interdiscip. Rev. Comput. Mol. Sci. 2021, 11, e1517], 280 vertical transition energies of 35 medium-sized molecules are computed using the SA-DSRG derived second- and third-order perturbation theories (PT2/PT3) along with a nonperturbative approach [sq-LDSRG(2)]. Comparing to the theoretical best estimates, the optimal flow parameter is found to be 0.35 and 2.0 Eh-2 for SA-DSRG-PT2 and SA-DSRG-PT3, respectively. For SA-sq-LDSRG(2), a flow parameter of 1.5 Eh-2 provides converged equations without compromising the accuracy. We then assess the accuracy of the SA-DSRG hierarchy using these parameters. The SA-DSRG-PT2 scheme outperforms the level-shifted CASPT2 by 0.10 eV in mean absolute error (MAE), yet this accuracy is slightly inferior than that of CASPT2 with the ionization-potential-electron-affinity shift. Both SA-DSRG-PT3 and SA-sq-LDSRG(2) yield a MAE of 0.10 eV, which is comparable to that of CASPT3 (0.09 eV). Finally, we compute vertical excitation energies of several low-lying singlet states of nucleobases. The SA-sq-LDSRG(2) approach provides highly accurate results for π → π* excitations, while n → π* transitions are better described by SA-DSRG-PT3.

Automatic derivation of many-body theories based on general Fermi vacua

This paper describes Wick&d, an implementation of the algebra of second-quantized operators normal ordered with respect to general correlated references and the corresponding Wick theorem [D. Mukherjee, Chem. Phys. Lett. 274, 561 (1997) and W. Kutzelnigg and D. Mukherjee, J. Chem. Phys. 107, 432 (1997)]. Wick&d employs a compact representation of operators and a backtracking algorithm to efficiently evaluate Wick contractions. Since Wick&d can handle both fully and partially contracted terms, it can be applied to both projective and Fock-space many-body formalisms. To demonstrate the usefulness of  Wick&d, we use it to evaluate the single-reference coupled cluster equations up to octuple excitations and report an automated derivation and implementation of the second-order driven similarity renormalization group multireference perturbation theory.

Geometric interpretation for coupled-cluster theory. A comparison of accuracy with the corresponding configuration interaction model

Although coupled-cluster theory is well-known for its accuracy, the geometry associated to the manifold of wave functions reached by the coupled-cluster ansatz has not been deeply explored. In this article we look for an interpretation for the high accuracy of coupled-cluster theory based on how the manifold of coupled-cluster wave functions is embedded within the space of n-electron wave functions. We define the coupled-cluster and configuration interaction manifolds and measure the distances from the full configuration interaction (FCI) wave function to these manifolds. We clearly observe that the FCI wave function is closer to the coupled-cluster manifold, that is curved, than to the configuration interaction manifold, that is flat, for the selected systems studied in this work. Furthermore, the decomposition of the distances among these manifolds and wave functions into excitation ranks gives insights on the failure of the coupled-cluster approach for multireference systems. The present results show a new interpretation for the quality of the coupled-cluster method, as contrasted to the truncated configuration interaction approach, besides the well-established argument based on size-extensivity. Furthermore, we show how a geometric description of wave function methods can be used in electronic structure theory.

A generalized hybrid scheme for multireference methods

A generalization of the hybrid scheme for multireference methods as recently put forward by Saitow and Yanai [J. Chem. Phys. 152, 114 111 (2020)] is presented. The hybrid methods are constructed by defining internal and external excitation spaces and evaluating these two subsets of excitations at different levels of theory. New hybrids that use the mix of internally contracted multireference coupled-cluster, unshifted multireference coupled electron pair, and multireference perturbation methods are derived and benchmarked. A new separation of the excitation space, which combines all singles and doubles excitations to the virtual orbitals into the external space, is also presented and tested. In general, the hybrid methods improve upon their non-hybrid parent method and offer a good compromise between computational complexity and numerical accuracy.

Second-Order CASSCF Algorithm with the Cholesky Decomposition of the Two-Electron Integrals

In this contribution, we present the implementation of a second-order complete active space–self-consistent field (CASSCF) algorithm in conjunction with the Cholesky decomposition of the two-electron repulsion integrals. The algorithm, called norm-extended optimization, guarantees convergence of the optimization, but it involves the full Hessian and is therefore computationally expensive. Coupling the second-order procedure with the Cholesky decomposition leads to a significant reduction in the computational cost, reduced memory requirements, and an improved parallel performance. As a result, CASSCF calculations of larger molecular systems become possible as a routine task. The performance of the new implementation is illustrated by means of benchmark calculations on molecules of increasing size, with up to about 3000 basis functions and 14 active orbitals.

The Molpro quantum chemistry package

Molpro is a general purpose quantum chemistry software package with a long development history. It was originally focused on accurate wavefunction calculations for small molecules but now has many additional distinctive capabilities that include, inter alia, local correlation approximations combined with explicit correlation, highly efficient implementations of single-reference correlation methods, robust and efficient multireference methods for large molecules, projection embedding, and anharmonic vibrational spectra. In addition to conventional input-file specification of calculations, Molpro calculations can now be specified and analyzed via a new graphical user interface and through a Python framework.

Benchmarks for Electronically Excited States with CASSCF Methods

The accuracy of three different complete active space (CAS) self-consistent field (CASSCF) methods is investigated for the electronically excited-state benchmark set of Schreiber , M. ; et al. J. Chem. Phys. 2008 , 128 , 134110 . Comparison of the CASSCF linear response (LR) methods MC-RPA and MC-TDA and the state-averaged (SA) CASSCF method is made for 122 singlet excitation energies and 69 oscillator strengths. Of all CASSCF methods, when considering the complete test set, MC-RPA performs best for both excitation energies and oscillator strengths with a mean absolute error (MAE) of 0.74 eV and 51%, respectively. MC-TDA and SA-CASSCF show a similar accuracy for the excitation energies with a MAE of ∼1 eV with respect to more accurate coupled cluster (CC3) excitation energies. The opposite trend is observed for the subset of n → π* excitation energies for which SA-CASSCF exhibits the least deviations (MAE 0.65 eV). By looking at s-tetrazine in more detail, we conclude that better performance for the n → π* SA-CASSCF excitation energies can be attributed to a fortunate error compensation. For oscillator strengths, SA-CASSCF performs worst for the complete test set (MAE 100%) as well as for the subsets of n → π* (MAE 192%) and π → π* excitations (MAE 84.9%). In general, CASSCF gives the worst performance for excitation energies of all excited-state ab initio methods considered so far due to lacking the major part of dynamic electron correlation, though MC-RPA and TD-DFT (BP86) show similar performance. Among all LR-type methods, LR-CASSCF oscillator strengths are the ones with the least accuracy for the same reason. As state-specific orbital relaxation effects are accounted for in LR-CASSCF, oscillator strengths are significantly more accurate than those of MS-CASPT2. Our findings should encourage further developments of response theory-based multireference methods with higher accuracy and feasibility.

CASSCF linear response calculations for large open-shell molecules

The complete active space self-consistent-field (CASSCF) linear response method for the simulation of ultraviolet-visible (UV/Vis) absorption and electronic circular dichroism (ECD) spectra of large open-shell molecules is presented. By using a one-index transformed Hamiltonian, the computation of the most time-consuming intermediates can be pursued in an integral-direct fashion, which allows us to employ the efficient resolution-of-the-identity and overlap-fitted chain-of-spheres approximation. For the iterative diagonalization, pairs of Hermitian and anti-Hermitian trial vectors are used which facilitate, on the one hand, an efficient solution of the pair-structured generalized eigenvalue problem in the reduced space, and on the other hand, make the full multiconfigurational random phase approximation as efficient as the corresponding Tamm-Dancoff approximation. Electronic transitions are analyzed and characterized in the particle-hole picture by natural transition orbitals that are introduced for CASSCF linear response theory. For a small organic radical, we can show that the accuracy of simulated UV/Vis absorption spectra with the CASSCF linear response approach is significantly improved compared to the popular state-averaged CASSCF method. To demonstrate the efficiency of the implementation, the 50 lowest roots of a large Ni triazole complex with 231 atoms are computed for the simulated UV/Vis and ECD spectra.

Revisiting the F + HCl → HF + Cl reaction using a multireference coupled-cluster method

A potential energy surface for the title reaction is constructed using a multireference coupled-cluster method, giving rate constant in excellent agreement with experiments.