Alexander Yu. Sokolov

A time-dependent formulation of multi-reference perturbation theory

We discuss the time-dependent formulation of perturbation theory in the context of the interacting zeroth-order Hamiltonians that appear in multi-reference situations. As an example, we present a time-dependent formulation and implementation of second-order n-electron valence perturbation theory. The resulting time-dependent n-electron valence second-order perturbation theory (t-NEVPT2) method yields the fully uncontracted n-electron valence perturbation wavefunction and energy, but has a lower computational scaling than the usual contracted variants, and also avoids the construction of high-order density matrices and the diagonalization of metrics. We present results of t-NEVPT2 for the water, nitrogen, carbon, and chromium molecules and outline directions for the future.

Orbital-optimized density cumulant functional theory

In density cumulant functional theory (DCFT) the electronic energy is evaluated from the one-particle density matrix and two-particle density cumulant, circumventing the computation of the wavefunction. To achieve this, the one-particle density matrix is decomposed exactly into the mean-field (idempotent) and correlation components. While the latter can be entirely derived from the density cumulant, the former must be obtained by choosing a specific set of orbitals. In the original DCFT formulation [W. Kutzelnigg, J. Chem. Phys. 125, 171101 (2006)] the orbitals were determined by diagonalizing the effective Fock operator, which introduces partial orbital relaxation. Here we present a new orbital-optimized formulation of DCFT where the energy is variationally minimized with respect to orbital rotations. This introduces important energy contributions and significantly improves the description of the dynamic correlation. In addition, it greatly simplifies the computation of analytic gradients, for which expressions are also presented. We offer a perturbative analysis of the new orbital stationarity conditions and benchmark their performance for a variety of chemical systems.

Free Cyclooctatetraene Dianion: Planarity, Aromaticity, and Theoretical Challenges

The planarity and 10 π-electron aromaticity of the free cyclooctatetraene dianion (C8H82–, COT2–) have been questioned recently on the basis of conflicting density functional and second-order Moller–Plesset perturbation computations. Rigorous coupled-cluster methods are employed here to establish the structure and properties of COT2–. Like many multiply charged anions, COT2– exists in isolation only as a short-lived resonance state lying above neutral COT. Wave function stability analysis demonstrates that predictions of nonplanar COT2– rings are artifacts of using overly diffuse basis sets. The resulting broken-symmetry wave functions are not characteristic of COT2– but mainly describe COT in a continuum of free electrons. All-electron coupled cluster theory extended through triple excitations [AE-CCSD(T)] yields a planar D8h symmetry COT2– structure. Final focal point analyses place the COT2– resonance state 61.6 kcal mol–1 above neutral COT. Nonetheless, COT2– exhibits structural, magnetic, and energet...

Analytic gradients for density cumulant functional theory: The DCFT-06 model

Density cumulant functional theory (DCFT) is one of a number of nascent electron correlation methods that are derived from reduced density matrices and cumulants thereof, instead of the wavefunction. Deriving properties from the density cumulant naturally yields methods that are size extensive and size consistent. In this work, we derive expressions for the analytic gradient, with respect to an external perturbation, for the DCFT-06 variant of density cumulant functional theory. Despite the fact that the DCFT-06 energy functional is stationary with respect to the density cumulant, the analytic gradients of the energy require the solution of perturbation-independent equations for both orbital and cumulant response. These two sets of linear response equations are coupled in nature and are solved iteratively with the solution of orbital and cumulant response equations each macroiteration, exhibiting rapid convergence. The gradients are implemented and benchmarked against coupled cluster theory with single and double excitations (CCSD) and CCSD with perturbative triple excitations [CCSD(T)], as well as accurate empirically corrected experimental data, for a test set comprising 15 small molecules. For most of the test cases, results from DCFT-06 are closer to CCSD(T) and empirical data than those from CCSD. Although the total energy and analytic gradient have the same asymptotic scaling, the present experience shows that the computational cost of the gradient is significantly lower.

Benchmark Study of Density Cumulant Functional Theory: Thermochemistry and Kinetics

We present an extensive benchmark study of density cumulant functional theory (DCFT) for thermochemistry and kinetics of closed- and open-shell molecules. The performance of DCFT methods (DC-06, DC-12, ODC-06, and ODC-12) is compared to that of coupled-electron pair methods (CEPA0 and OCEPA0) and coupled-cluster theory (CCSD and CCSD(T)) for the description of noncovalent interactions (A24 database), barrier heights of hydrogen-transfer reactions (HTBH38), radical stabilization energies (RSE30), adiabatic ionization energies (AIE), and covalent bond stretching in diatomic molecules. Our results indicate that out of four DCFT methods the ODC-12 method is the most reliable and accurate DCFT formulation to date. Compared to CCSD, ODC-12 shows superior results for all benchmark tests employed in our study. With respect to coupled-pair theories, ODC-12 outperforms CEPA0 and shows similar accuracy to the orbital-optimized CEPA0 variant (OCEPA0) for systems at equilibrium geometries. For covalent bond stretching, ODC-12 is found to be more reliable than OCEPA0. For the RSE30 and AIE data sets, ODC-12 shows competitive performance with CCSD(T). In addition to benchmark results, we report new reference values for the RSE30 data set computed using coupled cluster theory with up to perturbative quadruple excitations.

Density cumulant functional theory: the DC-12 method, an improved description of the one-particle density matrix.

Density cumulant functional theory (DCFT) is a theory that, in principle, can compute energies and properties exactly without a wavefunction. To accomplish this, the energy is expressed as an exact, known functional of the one-particle density matrix and two-particle density cumulant. The correlation contribution to the one-particle density matrix is obtained from the cumulant, to eliminate redundancy in the equations. The previous formulation of DCFT introduced this relationship in an approximate way, to obtain tractable equations. In this research, it is demonstrated that the correlation contribution to the one-particle density matrix can be extracted exactly from the cumulant, with minimal computational overhead and no increase in the asymptotic cost of the theory. We present numerical results, showing the improvements resulting from this reformulation (DC-12), and offer a perturbative analysis of the new equations to compare them to their predecessors.

Time-dependent N-electron valence perturbation theory with matrix product state reference wavefunctions for large active spaces and basis sets: Applications to the chromium dimer and all-trans polyenes

In earlier work [A. Y. Sokolov and G. K.-L. Chan, J. Chem. Phys. 144, 064102 (2016)], we introduced a time-dependent formulation of the second-order N-electron valence perturbation theory (t-NEVPT2) which (i) had a lower computational scaling than the usual internally contracted perturbation formulation and (ii) yielded the fully uncontracted NEVPT2 energy. Here, we present a combination of t-NEVPT2 with a matrix product state (MPS) reference wavefunction (t-MPS-NEVPT2) that allows us to compute uncontracted dynamic correlation energies for large active spaces and basis sets, using the time-dependent density matrix renormalization group algorithm. In addition, we report a low-scaling MPS-based implementation of strongly contracted NEVPT2 (sc-MPS-NEVPT2) that avoids computation of the four-particle reduced density matrix. We use these new methods to compute the dissociation energy of the chromium dimer and to study the low-lying excited states in all-trans polyenes (C4H6 to C24H26), incorporating dynamic correlation for reference wavefunctions with up to 24 active electrons and orbitals.

Can Density Cumulant Functional Theory Describe Static Correlation Effects

We evaluate the performance of density cumulant functional theory (DCT) for capturing static correlation effects. In particular, we examine systems with significant multideterminant character of the electronic wave function, such as the beryllium dimer, diatomic carbon, m-benzyne, 2,6-pyridyne, twisted ethylene, as well as the barrier for double-bond migration in cyclobutadiene. We compute molecular properties of these systems using the ODC-12 and DC-12 variants of DCT and compare these results to multireference configuration interaction and multireference coupled-cluster theories, as well as single-reference coupled-cluster theory with single, double (CCSD), and perturbative triple excitations [CCSD(T)]. For all systems the DCT methods show intermediate performance between that of CCSD and CCSD(T), with significant improvement over the former method. In particular, for the beryllium dimer, m-benzyne, and 2,6-pyridyne, the ODC-12 method along with CCSD(T) correctly predict the global minimum structures, while CCSD predictions fail qualitatively, underestimating the multireference effects. Our results suggest that the DC-12 and ODC-12 methods are capable of describing emerging static correlation effects but should be used cautiously when highly accurate results are required. Conveniently, the appearance of multireference effects in DCT can be diagnosed by analyzing the DCT natural orbital occupations, which are readily available at the end of the energy computation.

A transformed framework for dynamic correlation in multireference problems

We describe how multireference dynamic correlation theories can be naturally obtained as single-reference correlation theories in a canonically transformed frame. Such canonically transformed correlation theories are very simple and involve identical expressions to their single-reference counterparts. The corresponding excitations involve quasiparticles rather than the bare particles of the system. High-order density matrices (or their approximations) and the numerical metric instabilities common to multireference correlation theories do not appear. As an example, we formulate the Bogoliubov canonically transformed version of second-order Møller-Plesset perturbation theory and demonstrate its performance in H2, H2O, N2, and BeH2 bond dissociation.

Density cumulant functional theory from a unitary transformation: N-representability, three-particle correlation effects, and application to O4+

A new approach to density cumulant functional theory is developed that derives density cumulant N-representability conditions from an approximate Fock space unitary transformation. We present explicit equations for the third- and fourth-order two-particle cumulant N-representability, as well as the second-order contributions that depend on the connected three-particle density cumulant. These conditions are used to formulate the ODC-13 method and the non-iterative (λ3) correction that employ an incomplete description of the fourth-order two-particle cumulant N-representability and the second-order three-particle correlation effects, respectively. We perform an analysis of the ODC-13 N-representability description for the dissociation of H2 and apply the ODC-13 method and the (λ3) correction to diatomic molecules with multiple bond character and the symmetry-breaking tetraoxygen cation (O4(+)). For the O4(+) molecule, the vibrational frequencies of the ODC-13(λ3) method do not exhibit spatial symmetry breaking and are in a good agreement with the recent infrared photodissociation experiment. We report the O4(+) equilibrium structure, harmonic frequencies, and dissociation energy computed using ODC-13(λ3) with a diffuse, core-correlated aug-cc-pCVTZ basis set.