Adaptive Configuration Interaction for Computing Challenging Electronic Excited States with Tunable Accuracy

Quantum Equation-of-Motion Method with Single, Double, and Triple Excitations

The quantum equation-of-motion (qEOM) method with single and double excitations (qEOM-SD) has been proposed to study electronically excited states, but it fails to handle states dominated by double excitations. In this work, we reformulate the qEOM method within the effective Hamiltonian framework that satisfies the killer condition, and then present an efficient implementation incorporating single, double, and triple excitations. To reduce computational complexity, we employ point-group symmetry and perturbation theory to screen triple excitations, effectively reducing the scaling from No6Nv6 to No5Nv5, where No and Nv are the numbers of occupied and virtual spin orbitals, respectively. Furthermore, we account for the effect of neglected triple excitations by introducing a perturbative correction to the excitation energy. We apply this method to challenging cases where the qEOM-SD method exhibits significant errors, such as the 2 1Δ state of CH+ and the 2 1Σ state of HF. Our new method achieves energy errors below 0.18 eV while incorporating less than 8.2% of triple excitations. Additionally, we extend the operator screening technique to the quantum subspace expansion method for the efficient inclusion of selected triple excitations.

Compactification of determinant expansions via transcorrelation

Although selected configuration interaction (SCI) algorithms can tackle much larger Hilbert spaces than the conventional full CI method, the scaling of their computational cost with respect to the system size remains inherently exponential. In addition, inaccuracies in describing the correlation hole at small interelectronic distances lead to the slow convergence of the electronic energy relative to the size of the one-electron basis set. To alleviate these effects, we show that the non-Hermitian, transcorrelated (TC) version of SCI significantly compactifies the determinant space, allowing us to reach a given accuracy with a much smaller number of determinants. Furthermore, we note a significant acceleration in the convergence of the TC-SCI energy as the basis set size increases. The extent of this compression and the energy convergence rate are closely linked to the accuracy of the correlation factor used for the similarity transformation of the Coulombic Hamiltonian. Our systematic investigation of small molecular systems in increasingly large basis sets illustrates the magnitude of these effects.

Selected Configuration Interaction for Resonances

Electronic resonances are metastable states that can decay by electron loss. They are ubiquitous across various fields of science, such as chemistry, physics, and biology. However, current theoretical and computational models for resonances cannot yet rival the level of accuracy achieved by bound-state methodologies. Here, we generalize selected configuration interaction (SCI) to treat resonances by using the complex absorbing potential (CAP) technique. By modifying the selection procedure and the extrapolation protocol of standard SCI, the resulting CAP-SCI method yields resonance positions and widths of full configuration interaction quality. Initial results for the shape resonances of N2- and CO- reveal the important effect of high-order correlation, which shifts the values obtained with CAP-augmented equation-of-motion coupled-cluster with singles and doubles by more than 0.1 eV. The present CAP-SCI approach represents a cornerstone in the development of highly accurate methodologies for resonances.

Forte: A suite of advanced multireference quantum chemistry methods

Forte is an open-source library specialized in multireference electronic structure theories for molecular systems and the rapid prototyping of new methods. This paper gives an overview of the capabilities of Forte, its software architecture, and examples of applications enabled by the methods it implements.

Reference Energies for Double Excitations: Improvement and Extension

In the realm of photochemistry, the significance of double excitations (also known as doubly excited states), where two electrons are concurrently elevated to higher energy levels, lies in their involvement in key electronic transitions essential in light-induced chemical reactions as well as their challenging nature from the computational theoretical chemistry point of view. Based on state-of-the-art electronic structure methods (such as high-order coupled-cluster, selected configuration interaction, and multiconfigurational methods), we improve and expand our prior set of accurate reference excitation energies for electronic states exhibiting a substantial amount of double excitations [Loos et al. J. Chem. Theory Comput. 2019, 15, 1939]. This extended collection encompasses 47 electronic transitions across 26 molecular systems that we separate into two distinct subsets: (i) 28 "genuine" doubly excited states where the transitions almost exclusively involve doubly excited configurations and (ii) 19 "partial" doubly excited states which exhibit a more balanced character between singly and doubly excited configurations. For each subset, we assess the performance of high-order coupled-cluster (CC3, CCSDT, CC4, and CCSDTQ) and multiconfigurational methods (CASPT2, CASPT3, PC-NEVPT2, and SC-NEVPT2). Using as a probe the percentage of single excitations involved in a given transition (%T1) computed at the CC3 level, we also propose a simple correction that reduces the errors of CC3 by a factor of 3, for both sets of excitations. We hope that this more complete and diverse compilation of double excitations will help future developments of electronic excited-state methodologies.

Coupling molecular density functional theory with converged selected configuration interaction methods to study excited states in aqueous solution

This paper presents the first implementation of a coupling between advanced wavefunction theories and molecular density functional theory (MDFT). This method enables the modeling of solvent effect into quantum mechanical (QM) calculations by incorporating an electrostatic potential generated by solvent charges into the electronic Hamiltonian. Solvent charges are deduced from the spatially and angularly dependent solvent particle density. Such a density is obtained through the minimization of the functional associated with the molecular mechanics (MM) Hamiltonian describing the interaction between the fluid particles. The introduced QM/MDFT framework belongs to QM/MM family of methods, but its originality lies in the use of MDFT as the MM solver, offering two main advantages. First, its functional formulation makes it competitive with respect to sampling-based molecular mechanics. Second, it preserves a molecular-level description lost in macroscopic continuum approaches. The excited state properties of water and formaldehyde molecules solvated into water have been computed at the selected configuration interaction (SCI) level. The excitation energies and dipole moments have been compared with experimental data and previous theoretical work. A key finding is that using the Hartree-Fock method to describe the solute allows for predicting the solvent charge around the ground state with sufficient precision for the subsequent SCI calculations of excited states. This significantly reduces the computational cost of the described procedure, paving the way for the study of more complex molecules.

Using a multistate mapping approach to surface hopping to predict the ultrafast electron diffraction signal of gas-phase cyclobutanone

Using the recently developed multistate mapping approach to surface hopping (multistate MASH) method combined with SA(3)-CASSCF(12,12)/aug-cc-pVDZ electronic structure calculations, the gas-phase isotropic ultrafast electron diffraction (UED) of cyclobutanone is predicted and analyzed. After excitation into the n-3s Rydberg state (S2), cyclobutanone can relax through two S2/S1 conical intersections, one characterized by compression of the CO bond and the other by dissociation of the α-CC bond. Subsequent transfer into the ground state (S0) is then achieved via two additional S1/S0 conical intersections that lead to three reaction pathways: α ring-opening, ethene/ketene production, and CO liberation. The isotropic gas-phase UED signal is predicted from the multistate MASH simulations, allowing for a direct comparison to the experimental data. This work, which is a contribution to the cyclobutanone prediction challenge, facilitates the identification of the main photoproducts in the UED signal and thereby emphasizes the importance of dynamics simulations for the interpretation of ultrafast experiments.

Distributed Implementation of Full Configuration Interaction for One Trillion Determinants

Full configuration interaction (FCI) can provide an exact molecular ground-state energy within a given basis set and serve as a benchmark for approximate methods in quantum chemical calculations, including the emerging variational quantum eigensolver. However, its exponential computational and memory requirements easily exceed the capability of a single server and limit its applicability to large molecules. In this paper, we present a distributed FCI implementation employing a hybrid parallelization scheme with multithreading and multiprocessing to expand FCI's applicability. We optimize this scheme to minimize the bottlenecks arising from interprocess communications and interthread data management. Our implementation achieves higher scalability than the naive combination of prior works and successfully calculates the exact energy of C3H8/STO-3G with 1.31 trillion determinants, which is the largest FCI calculation to the best of our knowledge. Furthermore, we provide a comprehensive list of FCI results with 136 combinations of molecules and basis sets for future evaluation and development of approximate methods.

Modular Approach to Selected Configuration Interaction in an Arbitrary Spin Basis: Implementation and Comparison of Approaches

A modular selected configuration interaction (SCI) code has been developed that is based on the existing Monte-Carlo configuration interaction code (MCCI). The modularity allows various selection protocols to be implemented with ease and allows for fair comparison between wave functions built via different criteria. We have initially implemented adaptations of existing SCI theories, which are based on either energy- or coefficient-driven selection schemes. These codes have been implemented not only in the basis of Slater determinants (SDs) but also in the basis of configuration state functions (CSFs) and extended to state-averaged regimes. This allows one to take advantage of the reduced dimensionality of the wave function in the CSF basis and also the guarantee of pure spin states. All SCI methods were found to be able to predict potential energy surfaces to high accuracy, producing compact wave functions, when compared to full configuration interaction (FCI) for a variety of bond-breaking potential energy surfaces. The compactness of the error-controlled adaptive configuration interaction approach, particularly in the CSF basis, was apparent with nonparallelity errors within chemical accuracy while containing as little as 0.02% of the FCI CSF space. The size-to-accuracy was also extended to FCI spaces approaching one billion configurations.

Analytic Non-adiabatic Couplings for Selected Configuration Interaction via Approximate Degenerate Coupled Perturbed Hartree-Fock

We use degenerate perturbation theory and assume that for degenerate pairs of orbitals, the coupled perturbed Hartree-Fock coefficients are symmetric in the degenerate basis to show [Formula: see text] is the only modification needed in the original molecular orbital basis. This enables us to develop efficient and accurate analytic nonadiabatic couplings between electronic states for selected configuration interactions (CIs). Even when the states belong to different irreducible representations, degenerate orbital pairs cannot be excluded by symmetry. For various excited states of carbon monoxide and trigonal planar ammonia, we benchmark the method against the full CI and find it to be accurate. We create a semi-numerical approach and use it to show that the analytic approach is correct even when a high-symmetry structure is distorted to break symmetry so that near degeneracies in orbitals occur. For a range of geometries of trigonal planar ammonia, we find that the analytic non-adiabatic couplings for selected CI can achieve sufficient accuracy using a small fraction of the full CI space.

Coordinate Descent Full Configuration Interaction for Excited States

An efficient excited state method, named xCDFCI, in the configuration interaction framework is proposed. xCDFCI extends the unconstrained nonconvex optimization problem in coordinate descent full configuration interaction (CDFCI) to a multicolumn version for low-lying excited states computation. The optimization problem is addressed via a tailored coordinate descent method. In each iteration, a determinant is selected based on an approximated gradient, and coefficients of all states associated with the selected determinant are updated. A deterministic compression is applied to limit memory usage. We test xCDFCI applied to H2O and N2 molecules under the cc-pVDZ basis set. For both systems, five low-lying excited states in the same symmetry sector are calculated, together with the ground state. xCDFCI also produces accurate binding curves of the carbon dimer in the cc-pVDZ basis with chemical accuracy, where the ground state and four excited states in the same symmetry sector are benchmarked.

Approximate Excited-State Geometry Optimization with the State-Averaged Adaptive Sampling Configuration Interaction Algorithm with Large Active Spaces

The selected configuration interaction (SCI) wave function is a useful approximation to the full configuration interaction (FCI) one. The adaptive sampling CI (ASCI) method is a deterministic SCI method. By combining ASCI and orbital optimization, the ASCI self-consistent field (ASCI-SCF) method, which is an approximation of the complete active space self-consistent field (CASSCF) method, can be formulated as well. However, their applicability has been tested mainly on the systems in their electronically ground states. In this work, we implement the state-average (SA) ansatz in ASCI-SCF calculations to calculate excited states. We also derive expressions for the approximate analytical gradient and implement them as a computer program. We demonstrate the applicability of the current method for calculating vertical and adiabatic excitation energies and optimizing the molecular geometries of thermally activated delayed fluorescence (TADF) molecules.

Intruder-free cumulant-truncated driven similarity renormalization group second-order multireference perturbation theory

Accurate multireference electronic structure calculations are important for constructing potential energy surfaces. Still, even in the case of low-scaling methods, their routine use is limited by the steep growth of the computational and storage costs as the active space grows. This is primarily due to the occurrence of three- and higher-body density matrices or, equivalently, their cumulants. This work examines the effect of various cumulant truncation schemes on the accuracy of the driven similarity renormalization group second-order multireference perturbation theory. We test four different levels of three-body reduced density cumulant truncations that set different classes of cumulant elements to zero. Our test cases include the singlet-triplet gap of CH2, the potential energy curves of the XΣg+1 and AΣu+3 states of N2, and the singlet-triplet splittings of oligoacenes. Our results show that both relative and absolute errors introduced by these cumulant truncations can be as small as 0.5 kcal mol-1 or less. At the same time, the amount of memory required is reduced from O(NA6) to O(NA5), where NA is the number of active orbitals. No additional regularization is needed to prevent the intruder state problem in the cumulant-truncated second-order driven similarity renormalization group multireference perturbation theory methods.

Converging high-level coupled-cluster energetics via adaptive selection of excitation manifolds driven by moment expansions

A novel approach to rapidly converging high-level coupled-cluster (CC) energetics in an automated fashion is proposed. The key idea is an adaptive selection of excitation manifolds defining higher--than--two-body components of the cluster operator inspired by CC(P;Q) moment expansions. The usefulness of the resulting methodology is illustrated by molecular examples where the goal is to recover the electronic energies obtained using the CC method with a full treatment of singly, doubly, and triply excited clusters (CCSDT) when the noniterative triples corrections to CCSD fail.

Corrigendum: Coupled cluster theory on modern heterogeneous supercomputers

[This corrects the article DOI: 10.3389/fchem.2023.1154526.].

Unitary Coupled Cluster: Seizing the Quantum Moment

Shallow, CNOT-efficient quantum circuits are crucial for performing accurate computational chemistry simulations on current noisy quantum hardware. Here, we explore the usefulness of noniterative energy corrections, based on the method of moments of coupled-cluster theory, for accelerating convergence toward full configuration interaction. Our preliminary numerical results relying on iteratively constructed ansätze suggest that chemically accurate energies can be obtained with substantially more compact circuits, implying enhanced resilience to gate and decoherence noise.

Analytic Gradients for Selected Configuration Interaction

We develop analytic gradients for selected configuration interaction wave functions. Despite all pairs of molecular orbitals now potentially having to be considered for the coupled perturbed Hartree-Fock equations, we show that degenerate orbital pairs belonging to different irreducible representations in the largest abelian subgroup do not need to be included and instabilities due to degeneracies are avoided. We introduce seminumerical gradients and use them to validate the analytic approach even when near degeneracies are present due to high-symmetry geometries being slightly distorted to break symmetry. The method is applied to carbon monoxide, ammonia, square planar H₄, hexagonal planar H₆, and methane for a range of bond lengths where we demonstrate that analytic gradients for selected configuration interaction can approach the quality of full configuration interaction yet only use a very small fraction of its determinants.

Advances in Parallel Heat Bath Configuration Interaction

Heat-bath configuration interaction (HCI) is a deterministic method that approaches the full CI limit at greatly reduced computational cost. In this work, computational improvements to the HCI algorithm are introduced targeting speed, parallel efficiency, and memory requirements. The new implementation introduces a hash function to distribute determinants and takes advantage of MPI and OpenMP for parallelism allowing for a (22e,168o) active space to be studied, which explicitly includes 2.39 × 10⁷ variational determinants and 8.95 × 10¹⁰ perturbative determinants. Benchmarks show up to 86% parallel efficiency of the perturbative step on 32 nodes (4096 cores) and a total efficiency of 74%. The new HCI implementation is benchmarked for accuracy against prior results and applied to study the triplet-quintet gap in the challenging [FeO(NH₃)₅]²⁺ complex.

Good Vibrations: Calculating Excited-State Frequencies Using Ground-State Self-Consistent Field Models

The use of Δ-self-consistent field (SCF) approaches for studying excited electronic states has received a renewed interest in recent years. In this work, the use of this scheme for calculating excited-state vibrational frequencies is examined. Results from Δ-SCF calculations for a set of representative molecules are compared with those obtained using configuration interaction with single substitutions (CIS) and time-dependent density functional theory (TD-DFT) methods. The use of an approximate spin purification model is also considered for cases where the excited-state SCF solution is spin-contaminated. The results of this work demonstrate that an SCF-based description of an excited-state potential energy surface can be an accurate and cost-effective alternative to CIS and TD-DFT methods.

Quantum Orbital Minimization Method for Excited States Calculation on a Quantum Computer

We propose a quantum-classical hybrid variational algorithm, the quantum orbital minimization method (qOMM), for obtaining the ground state and low-lying excited states of a Hermitian operator. Given parametrized ansatz circuits representing eigenstates, qOMM implements quantum circuits to represent the objective function in the orbital minimization method and adopts a classical optimizer to minimize the objective function with respect to the parameters in ansatz circuits. The objective function has an orthogonality constraint implicitly embedded, which allows qOMM to apply a different ansatz circuit to each input reference state. We carry out numerical simulations that seek to find excited states of H₂, LiH, and a toy model consisting of four hydrogen atoms arranged in a square lattice in the STO-3G basis with UCCSD ansatz circuits. Comparing the numerical results with existing excited states methods, qOMM is less prone to getting stuck in local minima and can achieve convergence with more shallow ansatz circuits.

Relativistic nonorthogonal configuration interaction: application to L2,3-edge X-ray spectroscopy

In this article, we develop a relativistic exact-two-component nonorthogonal configuration interaction (X2C-NOCI) for computing L-edge X-ray spectra. This article to our knowledge is the first time NOCI has been used for relativistic wave functions. A set of molecular complexes, including SF₆, SiCl₄ and [FeCl₆]^3-, are used to demonstrate the accuracy and computational scaling of the X2C-NOCI method. Our results suggest that X2C-NOCI is able to satisfactorily capture the main features of the L_2,3-edge X-ray absorption spectra. Excitations from the core require a large amount of orbital relaxation to yield reasonable energies and X2C-NOCI allows us to treat orbital optimization explicitly. However, the cost of computing the nonorthogonal coupling is higher than in conventional CI. Here, we propose an improved integral screening using overlap-scaled density combined with a continuous measure of the generalized Slater-Condon rules that allows us to estimate if an element is zero before attempting a two-electron integral contraction.

SOiCI and iCISO: combining iterative configuration interaction with spin-orbit coupling in two ways

The near-exact iCIPT2 approach for strongly correlated systems of electrons, which stems from the combination of iterative configuration interaction (iCI, an exact solver of full CI) with configuration selection for static correlation and second-order perturbation theory (PT2) for dynamic correlation, is extended to the relativistic domain. In the spirit of spin separation, relativistic effects are treated in two steps: scalar relativity is treated by the infinite-order, spin-free part of the exact two-component (X2C) relativistic Hamiltonian, whereas spin-orbit coupling (SOC) is treated by the first-order, Douglas-Kroll-Hess-like SOC operator derived from the same X2C Hamiltonian. Two possible combinations of iCIPT2 with SOC are considered, i.e., SOiCI and iCISO. The former treats SOC and electron correlation on an equal footing, whereas the latter treats SOC in the spirit of state interaction, by constructing and diagonalizing an effective spin-orbit Hamiltonian matrix in a small number of correlated scalar states. Both double group and time reversal symmetries are incorporated to simplify the computation. Pilot applications reveal that SOiCI is very accurate for the spin-orbit splitting (SOS) of heavy atoms, whereas the computationally very cheap iCISO can safely be applied to the SOS of light atoms and even of systems containing heavy atoms when SOC is largely quenched by ligand fields.

Multicomponent heat-bath configuration interaction with the perturbative correction for the calculation of protonic excited states

In this study, we extend the multicomponent heat-bath configuration interaction (HCI) method to excited states. Previous multicomponent HCI studies have been performed using only the variational stage of the HCI algorithm as they have largely focused on the calculation of protonic densities. Because this study focuses on energetic quantities, a second-order perturbative correction after the variational stage is essential. Therefore, this study implements the second-order Epstein-Nesbet correction to the variational stage of multicomponent HCI for the first time. Additionally, this study introduces a new procedure for calculating reference excitation energies for multicomponent methods using the Fourier-grid Hamiltonian (FGH) method, which should allow the one-particle electronic basis set errors to be better isolated from errors arising from an incomplete description of electron-proton correlation. The excited-state multicomponent HCI method is benchmarked by computing protonic excitations of the HCN and FHF^- molecules and is shown to be of similar accuracy to previous excited-state multicomponent methods such as the multicomponent time-dependent density-functional theory and equation-of-motion coupled-cluster theory relative to the new FGH reference values.

Removing the Deadwood from DFT/MRCI Wave Functions: The p-DFT/MRCI Method

The combined density functional theory and multireference configuration interaction (DFT/MRCI) method is a powerful tool for the calculation of excited electronic states of large molecules. There exists, however, a large amount of superfluous configurations in a typical DFT/MRCI wave function. We show that this deadwood may be effectively removed using a simple configuration pruning algorithm based on second-order Epstein-Nesbet perturbation theory. The resulting method, which we denote p-DFT/MRCI, is shown to result in orders of magnitude saving in computational timings, while retaining the accuracy of the original DFT/MRCI method.

High-level coupled-cluster energetics by merging moment expansions with selected configuration interaction

Inspired by our earlier semi-stochastic work aimed at converging high-level coupled-cluster (CC) energetics [J. E. Deustua, J. Shen, and P. Piecuch, Phys. Rev. Lett. 119, 223003 (2017) and J. E. Deustua, J. Shen, and P. Piecuch, J. Chem. Phys. 154, 124103 (2021)], we propose a novel form of the CC(P; Q) theory in which the stochastic Quantum Monte Carlo propagations, used to identify dominant higher-than-doubly excited determinants, are replaced by the selected configuration interaction (CI) approach using the perturbative selection made iteratively (CIPSI) algorithm. The advantages of the resulting CIPSI-driven CC(P; Q) methodology are illustrated by a few molecular examples, including the dissociation of F₂ and the automerization of cyclobutadiene, where we recover the electronic energies corresponding to the CC calculations with a full treatment of singles, doubles, and triples based on the information extracted from compact CI wave functions originating from relatively inexpensive Hamiltonian diagonalizations.

The Static-Dynamic-Static Family of Methods for Strongly Correlated Electrons: Methodology and Benchmarking

A series of methods (SDSCI, SDSPT2, iCI, iCIPT2, iCISCF(2), iVI, and iCAS) is introduced to accurately describe strongly correlated systems of electrons. Born from the (restricted) static-dynamic-static (SDS) framework for designing many-electron wave functions, SDSCI is a minimal multireference (MR) configuration interaction (CI) approach that constructs and diagonalizes a [Formula: see text] matrix for [Formula: see text] states, regardless of the numbers of orbitals and electrons to be correlated. If the full molecular Hamiltonian H in the QHQ block (which describes couplings between functions of the first-order interaction space Q) of the SDSCI CI matrix is replaced with a zeroth-order Hamiltonian [Formula: see text] before the diagonalization is taken, we obtain SDSPT2, a CI-like second-order perturbation theory (PT2). Unlike most variants of MRPT2, SDSPT2 treats single and multiple states in the same way and is particularly advantageous in the presence of near degeneracy. On the other hand, if the SDSCI procedure is repeated until convergence, we will have iterative CI (iCI), which can converge quickly from the above to the exact solutions (full CI) even when starting with a poor guess. When further combined with the selection of important configurations followed by a PT2 treatment of dynamic correlation, iCI becomes iCIPT2, which is a near-exact theory for medium-sized systems. The microiterations of iCI for relaxing the coefficients of contracted many-electron functions can be generalized to an iterative vector interaction (iVI) approach for finding exterior or interior roots of a given matrix, in which the dimension of the search subspace is fixed by either the number of target roots or the user-specified energy window. Naturally, iCIPT2 can be employed as the active space solver of the complete active space (CAS) self-consistent field, leading to iCISCF(2), which can further be combined with iCAS for automated selection of active orbitals and assurance of the same CAS for all states and all geometries. The methods are calibrated by taking the Thiel set of benchmark systems as examples. Results for the corresponding cations, a new set of benchmark systems, are also reported.

Reinforcement Learning Configuration Interaction

Selected configuration interaction (sCI) methods exploit the sparsity of the full configuration interaction (FCI) wave function, yielding significant computational savings and wave function compression without sacrificing the accuracy. Despite recent advances in sCI methods, the selection of important determinants remains an open problem. We explore the possibility of utilizing reinforcement learning approaches to solve the sCI problem. By mapping the configuration interaction problem onto a sequential decision-making process, the agent learns on-the-fly which determinants to include and which to ignore, yielding a compressed wave function at near-FCI accuracy. This method, which we call reinforcement-learned configuration interaction, adds another weapon to the sCI arsenal and highlights how reinforcement learning approaches can potentially help solve challenging problems in electronic structure theory.

Software for the frontiers of quantum chemistry: An overview of developments in the Q-Chem 5 package

This article summarizes technical advances contained in the fifth major release of the Q-Chem quantum chemistry program package, covering developments since 2015. A comprehensive library of exchange-correlation functionals, along with a suite of correlated many-body methods, continues to be a hallmark of the Q-Chem software. The many-body methods include novel variants of both coupled-cluster and configuration-interaction approaches along with methods based on the algebraic diagrammatic construction and variational reduced density-matrix methods. Methods highlighted in Q-Chem 5 include a suite of tools for modeling core-level spectroscopy, methods for describing metastable resonances, methods for computing vibronic spectra, the nuclear-electronic orbital method, and several different energy decomposition analysis techniques. High-performance capabilities including multithreaded parallelism and support for calculations on graphics processing units are described. Q-Chem boasts a community of well over 100 active academic developers, and the continuing evolution of the software is supported by an "open teamware" model and an increasingly modular design.

Toward an Accurate Ab Initio Description of Low-Lying Singlet Excited States of Polyenes

The low-lying excited states of carotenoids play a crucial role in many important biophysical processes such as photosynthesis. Most of these excited states are strongly correlated, which makes them both challenging for a qualitative ab initio description and an engaging model system for trying out emerging multireference methods. Among these methods, driven similarity renormalization group (DSRG) and its perturbative version (DSRG-MRPT2) are especially attractive in terms of both accuracy and moderate numerical complexity. In this paper, we applied density matrix renormalization group (DMRG) followed by DSRG-MRPT2 for the calculation of vertical and adiabatic excitation energies into the 2A_g^-, 1B_u^-, and 1B_u⁺ electronic states of polyenes containing from 8 to 13 conjugating double bonds acting as a model for natural carotenoids. It was shown that the DSRG flow parameter should be adjusted to ensure both the energy convergence with respect to it and the agreement with the experimental data. With the increased flow parameter, the proposed combination of methods provides a reasonable agreement with the experiment. The deviations of the adiabatic excitation energies are less than 1000 cm^-1 for the 2A_g^- and less than 3000 cm^-1 for the excited states of the B_u symmetry, which in terms of accuracy significantly outperforms the N-electron valence state perturbation theory. At the same time, DSRG-MRPT2 is shown to be robust with respect to variation of quality of the DMRG reference wave function such as the orbital optimization or the number of electronic states in the averaging.

Is Externally Corrected Coupled Cluster Always Better Than the Underlying Truncated Configuration Interaction?

The short answer to the question in the title is "no". We identify classes of truncated configuration interaction (CI) wave functions for which the externally corrected coupled-cluster (ec-CC) approach using the three-body (T₃) and four-body (T₄) components of the cluster operator extracted from CI does not improve the results of the underlying CI calculations. Implications of our analysis, illustrated by numerical examples, for the ec-CC computations using truncated and selected CI methods are discussed. We also introduce a novel ec-CC approach using the T₃ and T₄ amplitudes obtained with the selected CI scheme abbreviated as CIPSI, correcting the resulting energies for the missing T₃ correlations not captured by CIPSI with the help of moment expansions similar to those employed in the completely renormalized CC methods.

High-level coupled-cluster energetics by Monte Carlo sampling and moment expansions: Further details and comparisons

We recently proposed a novel approach to converging electronic energies equivalent to high-level coupled-cluster (CC) computations by combining the deterministic CC(P;Q) formalism with the stochastic configuration interaction (CI) and CC Quantum Monte Carlo (QMC) propagations. This article extends our initial study [J. E. Deustua, J. Shen, and P. Piecuch, Phys. Rev. Lett. 119, 223003 (2017)], which focused on recovering the energies obtained with the CC method with singles, doubles, and triples (CCSDT) using the information extracted from full CI QMC and CCSDT-MC, to the CIQMC approaches truncated at triples and quadruples. It also reports our first semi-stochastic CC(P;Q) calculations aimed at converging the energies that correspond to the CC method with singles, doubles, triples, and quadruples (CCSDTQ). The ability of the semi-stochastic CC(P;Q) formalism to recover the CCSDT and CCSDTQ energies, even when electronic quasi-degeneracies and triply and quadruply excited clusters become substantial, is illustrated by a few numerical examples, including the F-F bond breaking in F₂, the automerization of cyclobutadiene, and the double dissociation of the water molecule.

Further Development of iCIPT2 for Strongly Correlated Electrons

The efficiency of the recently proposed iCIPT2 [iterative configuration interaction (iCI) with selection and second-order perturbation theory (PT2); J. Chem. Theory Comput. 2020, 16, 2296] for strongly correlated electrons is further enhanced (by up to 20×) by using (1) a new ranking criterion for configuration selection, (2) a new particle-hole algorithm for Hamiltonian construction over randomly selected configuration state functions (CSF), and (3) a new data structure for the quick sorting of the variational and first-order interaction spaces. Meanwhile, the memory requirement is also significantly reduced. As a result, this improved implementation of iCIPT2 can handle 1 order of magnitude more CSFs than the previous version, as revealed by taking the chromium dimer and an iron-sulfur cluster, [Fe₂S₂(SCH₃)]₄^2-, as examples.

Multicomponent CASSCF Revisited: Large Active Spaces Are Needed for Qualitatively Accurate Protonic Densities

Multicomponent methods seek to treat select nuclei, typically protons, fully quantum mechanically and equivalent to the electrons of a chemical system. In such methods, it is well-known that due to the neglect of electron-proton correlation, a Hartree-Fock (HF) description of the electron-proton interaction catastrophically fails leading to qualitatively incorrect protonic properties. In single-component quantum chemistry, the qualitative failure of HF is normally indicative of the need for multireference methods such as complete active space self-consistent field (CASSCF). While a multicomponent CASSCF method was implemented nearly 20 years ago, it is only able to perform calculations with very small active spaces (∼10⁵ multicomponent configurations). Therefore, in order to extend the realm of applicability of the multicomponent CASSCF method, this study derives and implements a new two-step multicomponent CASSCF method that uses multicomponent heat-bath configuration interaction for the configuration interaction step, enabling calculations with very large active spaces (up to 16 electrons in 48 orbitals). We find that large electronic active spaces are needed to obtain qualitatively accurate protonic densities for the HCN and FHF^- molecules. Additionally, the multicomponent CASSCF method implemented here should have further applications for double-well protonic potentials and systems that are inherently electronically multireference.

The Shape of Full Configuration Interaction to Come

We present a Perspective on what the future holds for full configuration interaction (FCI) theory, with an emphasis on conceptual rather than technical details. Upon revisiting the early history of FCI, a number of its key contemporary approximations are compared on as equal a footing as possible, using a recent blind challenge on the benzene molecule as a testbed [Eriksen et al., J. Phys. Chem. Lett., 2020 11, 8922]. In the process, we review the scope of applications for which FCI continues to prove indispensable, and the required traits in terms of robustness, efficacy, and reliability its modern approximations must satisfy are discussed. We close by conveying a number of general observations on the merits offered by the state-of-the-art alongside some of the challenges still faced to this day. While the field has altogether seen immense progress over the years-the past decade, in particular-it remains clear that our community as a whole has a substantial way to go in enhancing the overall applicability of near-exact electronic structure theory for systems of general composition and increasing size.

Coupled Cluster Externally Corrected by Adaptive Configuration Interaction

An externally corrected coupled cluster (CC) method, where an adaptive configuration interaction (ACI) wave function provides the external cluster amplitudes, named ACI-CC, is presented. By exploiting the connection between configuration interaction and CC through cluster analysis, the higher-order T₃ and T₄ terms obtained from ACI are used to augment the T₁ and T₂ amplitude equations from traditional CC. These higher-order contributions are kept frozen during the CC iterations and do not contribute to an increased cost with respect to coupled cluster including the single and double excitations (CCSD). We have benchmarked this method on three closed-shell systems: beryllium dimer, carbonyl oxide, and cyclobutadiene, with good results compared to other corrected CC methods. In all cases, the inclusion of these external corrections improved upon the "gold standard" CCSD(T) results, indicating that ACI-CCSD(T) can be used to assess strong correlation effects in a system and as an inexpensive starting point for more complex external corrections.

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-ζ quality. As a broad international endeavor, 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 mE_H. However, we find the root-mean-square deviation of the energies from the studied methods to be considerable (1.3 mE_H), which in light of the acclaimed performance of each of the methods for smaller molecular systems clearly displays the challenges faced in extending reliable, near-exact correlation methods to larger systems. While the discrepancies exposed by our study thus emphasize the fact that the current state-of-the-art approaches leave room for improvement, we still expect the present assessment to provide a valuable community resource for benchmark and calibration purposes going forward.

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 an entirely unbiased manner for chemical and physical problems of not only academic but also applied chemical interest. Given the massive parallelism and overall accuracy of the resulting method, we argue that generalized MBE-FCI theory possesses an immense potential to yield near-exact correlation energies for molecular systems of unprecedented size, composition, and complexity in the years to come.