anti-Hermitian Contracted Schrödinger Equation Research Articles

Quantum-classical hybrid algorithm for the simulation of all-electron correlation.

While chemical systems containing hundreds to thousands of electrons remain beyond the reach of quantum devices, hybrid quantum-classical algorithms present a promising pathway toward a quantum advantage. Hybrid algorithms treat the exponentially scaling part of the calculation-the static correlation-on the quantum computer and the non-exponentially scaling part-the dynamic correlation-on the classical computer. While a variety of algorithms have been proposed, the dependence of many methods on the total wave function limits the development of easy-to-use classical post-processing implementations. Here, we present a novel combination of quantum and classical algorithms, which computes the all-electron energy of a strongly correlated molecular system on the classical computer from the 2-electron reduced density matrix (2-RDM) evaluated on the quantum device. Significantly, we circumvent the wave function in the all-electron calculations by using density matrix methods that only require input of the statically correlated 2-RDM. Although the algorithm is completely general, we test it with two classical density matrix methods, the anti-Hermitian contracted Schrödinger equation (ACSE) and multiconfiguration pair-density functional theories, using the recently developed quantum ACSE method for simulating the statically correlated 2-RDM. We obtain experimental accuracy for the relative energies of all three benzyne isomers and thereby demonstrate the ability of the developed algorithm to achieve chemically relevant and accurate results on noisy intermediate-scale quantum devices.

Accurate singlet-triplet gaps in biradicals via the spin averaged anti-Hermitian contracted Schrödinger equation.

The accurate description of biradical systems, and in particular the resolution of their singlet-triplet gaps, has long posed a major challenge to the development of electronic structure theories. Biradicaloid singlet ground states are often marked by strong correlation and, hence, may not be accurately treated by mainstream, single-reference methods such as density functional theory or coupled cluster theory. The anti-Hermitian contracted Schrödinger equation (ACSE), whose fundamental quantity is the two-electron reduced density matrix rather than the N-electron wave function, has previously been shown to account for both dynamic and strong correlations when seeded with a strongly correlated guess from a complete active space (CAS) calculation. Here, we develop a spin-averaged implementation of the ACSE, allowing it to treat higher multiplicity states from the CAS input without additional state preparation. We apply the spin-averaged ACSE to calculate the singlet-triplet gaps in a set of small main group biradicaloids, as well as the organic four-electron biradicals trimethylenemethane and cyclobutadiene, and naphthalene, benchmarking the results against other state-of-the-art methods reported in the literature.

Using reduced density matrix techniques to capture static and dynamic correlation in the energy landscape for the decomposition of the CH2CH2ONO radical and support a non-IRC pathway.

The unexpected abundance of HNO in the photodecomposition of the radical 2-nitrosooxy ethyl (CH2CH2ONO) is investigated through calculations of the potential energy surface by the anti-Hermitian contracted Schrödinger equation (ACSE) method, which directly generates the 2-electron reduced density matrix. The ACSE, which is able to balance single-reference (dynamic) and multi-reference (static) correlation effects, reveals some subtle correlation effects along the intrinsic reaction coordinate (IRC) en route to NO + oxirane, an IRC which offers a potential bifurcation to the HNO + vinoxy product channel. These effects were not fully captured by either single-reference techniques, such as coupled cluster, or multi-reference techniques, such as second-order multi-reference perturbation theory. These correlation effects reveal small to moderate energy changes in key transition states, which have implications for the reaction mechanism as related to the production of HNO.

Enhanced computational efficiency in the direct determination of the two-electron reduced density matrix from the anti-Hermitian contracted Schrödinger equation with application to ground and excited states of conjugated π-systems.

Determination of the two-electron reduced density matrix (2-RDM) from the solution of the anti-Hermitian contracted Schrödinger equation (ACSE) yields accurate energies and properties for both ground and excited states. Here, we develop a more efficient method to solving the ACSE that uses second-order information to select a more optimal step towards the solution. Calculations on the ground and excited states of water, hydrogen fluoride, and conjugated π systems show that the improved ACSE algorithm is 10-20 times faster than the previous ACSE algorithm. The ACSE can treat both single- and multi-reference electron correlation with the initial 2-RDM from a complete-active-space self-consistent-field (CASSCF) calculation. Using the improved algorithm, we explore the relationship between truncation of the active space in the CASSCF calculation and the accuracy of the energy and 2-RDM from the ACSE calculation. The accuracy of the ACSE, we find, is less sensitive to the size of the active space than the accuracy of other wavefunction methods, which is useful when large active space calculations are computationally infeasible.

Highly accurate excited-state energies from direct computation of the 2-electron reduced density matrix by the anti-Hermitian contracted Schrödinger equation

ABSTRACTDirectly solving for the 2-electron reduced density matrix (2-RDM) via the anti-Hermitian contracted Schrödinger equation (ACSE) enables computations for excited states energies without the N-electron wave function. Of particular interest are excitations and dissociation curves that exhibit strong multi-reference correlation effects. The ground and excited states of the molecules HF, H2O, and N2 are examined at both equilibrium and non-equilibrium geometries to compare the ability of the ACSE and widely used ab initio techniques to treat strong multi-reference electron correlation. Calculations are performed with double-ζ basis sets for calibration with full configuration interaction (FCI). Multi-reference second-order perturbation theory (MRPT2) and the ACSE both provide qualitative precision with respect to FCI data, although the ACSE's capability to include higher order correlation effects permits near-FCI accuracy for ground and excited states and excitation energies.

Molecule-Optimized Basis Sets and Hamiltonians for Accelerated Electronic Structure Calculations of Atoms and Molecules

Molecule-optimized basis sets, based on approximate natural orbitals, are developed for accelerating the convergence of quantum calculations with strongly correlated (multireferenced) electrons. We use a low-cost approximate solution of the anti-Hermitian contracted Schrödinger equation (ACSE) for the one- and two-electron reduced density matrices (RDMs) to generate an approximate set of natural orbitals for strongly correlated quantum systems. The natural-orbital basis set is truncated to generate a molecule-optimized basis set whose rank matches that of a standard correlation-consistent basis set optimized for the atoms. We show that basis-set truncation by approximate natural orbitals can be viewed as a one-electron unitary transformation of the Hamiltonian operator and suggest an extension of approximate natural-orbital truncations through two-electron unitary transformations of the Hamiltonian operator, such as those employed in the solution of the ACSE. The molecule-optimized basis set from the ACSE improves the accuracy of the equivalent standard atom-optimized basis set at little additional computational cost. We illustrate the method with the potential energy curves of hydrogen fluoride and diatomic nitrogen. Relative to the hydrogen fluoride potential energy curve from the ACSE in a polarized triple-ζ basis set, the ACSE curve in a molecule-optimized basis set, equivalent in size to a polarized double-ζ basis, has a nonparallelity error of 0.0154 au, which is significantly better than the nonparallelity error of 0.0252 au from the polarized double-ζ basis set.

Photoexcited tautomerization of vinyl alcohol to acetylaldehydevia a conical intersection from contracted Schrödinger theory

The photoexcited tautomerization of vinyl alcohol to acetylaldehyde via a conical intersection is analyzed through the direct calculation of two-electron reduced density matrices (2-RDMs) from solutions to the anti-Hermitian contracted Schrödinger equation (ACSE). The study utilizes the recent generalization of the ACSE method for the treatment of excited states [G. Gidofalvi and D. A. Mazziotti, Phys. Rev. A, 2009, 80, 022507]. We computed absolute energies of the critical points as well as various intermediate points along the ground- and excited-state potential energy surface of vinyl alcohol and acetylaldehyde. The ACSE, MCSCF, second-order multireference many-body perturbation theory (MRMP2), and various coupled cluster methods all demonstrate the existence of a family of pathways from vinyl alcohol to acetylaldehyde via a conical intersection that are monotonically decreasing in energy. The conical intersection, proposed for the first time in this paper, is both structurally and energetically similar to the ground-state transition state. We observe a relationship between conical intersections and transition states both in this paper and in our previous work on bicyclobutane's ring conical intersection [J. W. Snyder, Jr. and D. A. Mazziotti, J. Chem. Phys., 2011, 135, 024107]. To treat multireference correlation, we seeded the ACSE with an initial 2-RDM guess from a multiconfiguration self-consistent field (MCSCF) calculation. The ACSE recovers more single-reference correlation energy than MRMP2 and more multireference correlation energy than comparable single-reference wave function methods. The 2-RDMs from the ACSE nearly satisfy all necessary N-representability conditions.

Conical Intersection of the Ground and First Excited States of Water: Energies and Reduced Density Matrices from the Anti-Hermitian Contracted Schrödinger Equation

A conical intersection between the ground and first-excited states of water is computed through the direct calculation of two-electron reduced density matrices (2-RDMs) from solutions of the anti-Hermitian contracted Schrödinger equation (ACSE). This study is an extension of a previous study in which the ACSE was used to compute the energies around a conical intersection in the triplet excited states of methylene [Snyder, J. W., Jr.; Rothman, A. E.; Foley, J. J.; Mazziotti, D. A. J. Chem. Phys. 2010, 132, 154109]. We compute absolute energies of the 1(1)A' and 2(1)A' states of water (H(2)O) and the location of the conical intersection. The ACSE energies are compared to those from ab initio wave function methods. To treat multireference correlation, we seed the ACSE with an initial 2-RDM from a multiconfiguration self-consistent field (MCSCF) calculation. Unlike the situation for methylene, the two states in the vicinity of the conical intersection of water both have the same spatial symmetry. Hence, the study demonstrates the ability of the ACSE to resolve states of the same spatial symmetry that are nearly degenerate in energy. The 2-RDMs from the ACSE nearly satisfy necessary N-representability conditions. Comparison of the results from double-ζ and augmented double-ζ basis sets demonstrates the importance of augmented (or diffuse) functions for determining the location of the conical intersection.

Photoexcited conversion of gauche-1,3-butadiene to bicyclobutane via a conical intersection: Energies and reduced density matrices from the anti-Hermitian contracted Schrödinger equation

The photoexcited reaction pathway from gauche-1,3-butadiene to bicyclobutane via a conical intersection is analyzed through the direct calculation of two-electron reduced density matrices (2-RDMs) from solutions to the anti-Hermitian contracted Schrödinger equation (ACSE). The study utilizes the recent generalization of the ACSE method for the treatment of excited states [G. Gidofalvi and D. A. Mazziotti, Phys. Rev. A 80, 022507 (2009)]. We computed absolute energies of the critical points as well as various intermediate points along the ground-and excited-state potential energy surface of gauche-1,3-butadiene and bicyclobutane. To treat multi-reference correlation, we seeded the ACSE with an initial 2-RDM from a multi-configuration self-consistent field (MCSCF) calculation. The ACSE, MCSCF, and second-order multi-reference many-body perturbation theory (MRPT2) all demonstrate that there exists a family of pathways from gauche-1,3-butadiene to bicyclobutane via a conical intersection that are monotonically decreasing in energy, confirming a conjecture by Sicilia et al. [J. Phys. Chem. A 111, 2182 (2007)]. The ACSE recovers more single-reference correlation energy than MRPT2 and more multi-reference correlation energy than comparable single-reference wave function methods. The 2-RDMs from the ACSE nearly satisfy all necessary N-representability conditions.

Balancing single- and multi-reference correlation in the chemiluminescent reaction of dioxetanone using the anti-Hermitian contracted Schrödinger equation

Direct computation of energies and two-electron reduced density matrices (2-RDMs) from the anti-Hermitian contracted Schrödinger equation (ACSE) [D. A. Mazziotti, Phys. Rev. Lett. 97, 143002 (2006)], it is shown, recovers both single- and multi-reference electron correlation in the chemiluminescent reaction of dioxetanone especially in the vicinity of the conical intersection where strong correlation is important. Dioxetanone, the light-producing moiety of firefly luciferin, efficiently converts chemical energy into light by accessing its excited-state surface via a conical intersection. Our previous active-space 2-RDM study of dioxetanone [L. Greenman and D. A. Mazziotti, J. Chem. Phys. 133, 164110 (2010)] concluded that correlating 16 electrons in 13 (active) orbitals is required for realistic surfaces without correlating the remaining (inactive) orbitals. In this paper we pursue two complementary goals: (i) to correlate the inactive orbitals in 2-RDMs along dioxetanone's reaction coordinate and compare these results with those from multireference second-order perturbation theory (MRPT2) and (ii) to assess the size of the active space-the number of correlated electrons and orbitals-required by both MRPT2 and ACSE for accurate energies and surfaces. While MRPT2 recovers very different amounts of correlation with (4,4) and (16,13) active spaces, the ACSE obtains a similar amount of correlation energy with either active space. Nevertheless, subtle differences in excitation energies near the conical intersection suggest that the (16,13) active space is necessary to determine both energetic details and properties. Strong electron correlation is further assessed through several RDM-based metrics including (i) total and relative energies, (ii) the von Neumann entropy based on the 1-electron RDM, as well as the (iii) infinity and (iv) squared Frobenius norms based on the cumulant 2-RDM.

Convergence and computational efficiency enhancements in the iterative solution of the G‐particle‐hole hypervirial equation

AbstractThe G‐particle‐hole hypervirial (GHV) equation has been recently reported (Valdemoro et al., Sixth International Congress of the International Society for Theoretical Chemical Physics Vancouver: Canada, 2008. Alcoba et al., Int J Quantum Chem 2009, 109, 3178; Valdemoro et al., Int J Quantum Chem 2009, 109, 2622). This equation is the newest member of the family of equations which can be obtained by applying a matrix‐contracting mapping (Valdemoro, An R Soc Esp Fís 1983, 79, 106; Valdemoro, Phys Rev A 1985, 31, 2114; Valdemoro, in Density Matrices and Density Functionals, Reidel: Dordrecht, 1987; p 275.) to the matrix representation in the N‐electron space of the Schrödinger, Liouville and hypervirial equations. The procedure that we have applied in order to solve the GHV equation exploits the stationary property of the hypervirials (Hirschfelder, J Chem Phys 1960, 33, 1462; Hirschfelder and Epstein, Phys Rev 1961, 123, 1495) and follows the general lines of Mazziotti's variational approach for solving the anti‐Hermitian contracted Schrödinger equation (ACSE) (Mazziotti, Phys Rev Lett 2006, 97, 143002; Mazziotti, Phys Rev A 2007, 75, 022505; Mazziotti, J Chem Phys 2007, 126, 184101). In this article, we report how the method's convergence has been significantly enhanced and how its computational scaling has been considerably reduced (in both floating‐point operations and storage). The results for a variety of atomic and molecular calculations confirming these methodological improvements are reported here. © 2010 Wiley Periodicals, Inc. Int J Quantum Chem, 2011

Strongly correlated mechanisms of a photoexcited radical reaction from the anti-Hermitian contracted Schrödinger equation

Photoexcited radical reactions are critical to processes in both nature and materials, and yet they can be challenging for electronic structure methods due to the presence of strong electron correlation. Reduced-density-matrix (RDM) methods, based on solving the anti-Hermitian contracted Schrödinger equation (ACSE) for the two-electron RDM (2-RDM), are examined for studying the strongly correlated mechanisms of these reactions with application to the electrocyclic interconversion of allyl and cyclopropyl radicals. We combine recent extensions of the ACSE to excited states [G. Gidofalvi and D. A. Mazziotti, Phys. Rev. A 80, 022507 (2009)] and arbitrary spin states [A. E. Rothman, J. J. Foley IV, and D. A. Mazziotti, Phys. Rev. A 80, 052508 (2009)]. The ACSE predicts that the ground-state ring closure of the allyl radical has a high 52.5 kcal/mol activation energy that is consistent with experimental data, while the closure of an excited allyl radical can occur by disrotatory and conrotatory pathways whose transition states are essentially barrierless. Comparisons are made with multireference second- and third-order perturbation theories and multireference configuration interaction. While predicted energy differences do not vary greatly between methods, the ACSE appears to improve these differences when they involve a strongly and a weakly correlated radical by capturing a greater share of single-reference correlation that increases the stability of the weakly correlated radicals. For example, the ACSE predicts a -39.6 kcal/mol conversion of the excited allyl radical to the ground-state cyclopropyl radical in comparison to the -32.6 to -37.3 kcal/mol conversions predicted by multireference methods. In addition, the ACSE reduces the computational scaling with the number of strongly correlated orbitals from exponential (traditional multireference methods) to quadratic. Computed ground- and excited-state 2-RDMs are nearly N-representable.

Efficient geometry optimization by Hellmann–Feynman forces with the anti-Hermitian contracted Schrödinger equation

An efficient method for geometry optimization based on solving the anti-Hermitian contracted Schrödinger equation (ACSE) is presented. We formulate a reduced version of the Hellmann–Feynman theorem (HFT) in terms of the two-electron reduced Hamiltonian operator and the two-electron reduced density matrix (2-RDM). The HFT offers a considerable reduction in computational cost over methods which rely on numerical derivatives. While previous geometry optimizations with numerical gradients required 2M evaluations of the ACSE where M is the number of nuclear degrees of freedom, the HFT requires only a single ACSE calculation of the 2-RDM per gradient. Synthesizing geometry optimization techniques with recent extensions of the ACSE theory to arbitrary electronic and spin states provides an important suite of tools for accurately determining equilibrium and transition-state structures of ground- and excited-state molecules in closed- and open-shell configurations. The ability of the ACSE to balance single- and multi-reference correlation is particularly advantageous in the determination of excited-state geometries where the electronic configurations differ greatly from the ground-state reference. Applications are made to closed-shell molecules N2, CO, H2O, the open-shell molecules B2 and CH, and the excited state molecules N2, B2, and BH. We also study the HCN ↔ HNC isomerization and the geometry optimization of hydroxyurea, a molecule which has a significant role in the treatment of sickle-cell anaemia.

Conical intersections in triplet excited states of methylene from the anti-Hermitian contracted Schrödinger equation

A conical intersection in triplet excited states of methylene is computed through the direct calculation of two-electron reduced density matrices (2-RDMs) from solutions of the anti-Hermitian contracted Schrodinger equation (ACSE). The study synthesizes recent extensions of the ACSE method for the treatment of excited states [G. Gidofalvi and D. A. Mazziotti, Phys. Rev. A 80, 022507 (2009)] and arbitrary-spin states [A. E. Rothman, J. J. Foley, and D. A. Mazziotti, Phys. Rev. A 80, 052508 (2009)]. We compute absolute energies of the 1 (3)B(1), 1 (3)A(2), and 2 (3)B(1) states of methylene (CH(2)) and the location of the conical intersection along the 1 (3)A(2)-2 (3)B(1) potential-energy surfaces. To treat multireference correlation, we seed the ACSE with an initial 2-RDM from a multiconfiguration self-consistent field (MCSCF) calculation. The ACSE produces energies that significantly improve upon those from MCSCF and second-order multireference many-body perturbation theory, and the 2-RDMs from the ACSE nearly satisfy necessary N-representability conditions. Comparison of the results from augmented double-zeta and triple-zeta basis sets demonstrates the importance of augmented (or diffuse) functions for determining the location of the conical intersection.

Nonequilibrium, steady-state electron transport with N-representable density matrices from the anti-Hermitian contracted Schrödinger equation

We study molecular conductivity for a one-electron, bath-molecule-bath model Hamiltonian. The primary quantum-mechanical variable is the one-electron reduced density matrix (1-RDM). By identifying similarities between the steady-state Liouville equation and the anti-Hermitian contracted Schrödinger equation (ACSE) [D. A. Mazziotti, Phys. Rev. A 75, 022505 (2007)], we develop a way of enforcing nonequilibrium, steady-state behavior in a time-independent theory. Our results illustrate the relationship between current and voltage in molecular junctions assuming that the total number of electrons under consideration can be fixed across all driving potentials. The impetus for this work is a recent study by Subotnik et al. that also uses the 1-RDM to study molecular conductivity under different assumptions regarding the total number of electrons [J. E. Subotnik et al., J. Chem. Phys. 130, 144105 (2009)]. Unlike calculations in the previous study, our calculations result in 1-RDMs that are fully N-representable. The present work maintains N-representability through a bath-bath mixing that is related to a time-independent relaxation of the baths in the absence of the molecule, as governed by the ACSE. A lack of N-representability can be important since it corresponds to occupying energy states in the molecule or baths with more than one electron or hole (the absence of an electron) in violation of the Pauli principle. For this reason the present work may serve as an important, albeit preliminary, step in designing a 2-RDM/ACSE method for studying steady-state molecular conductivity with an explicit treatment of electron correlation.

Energy Barriers of Vinylidene Carbene Reactions from the Anti-Hermitian Contracted Schrödinger Equation

Computational studies of carbenes must take into account the possibility of multireference correlation because the highest occupied and lowest unoccupied molecular orbitals can be nearly energetically degenerate. We apply the anti-Hermitian contracted Schrodinger equation (ACSE) [Mazziotti, D. A. Phys. Rev. Lett. 2006, 97, 143002] to compute two-electron reduced density matrices (2-RDMs) and their energies for two carbene reactions: (i) the acetylene-vinylidene rearrangement and (ii) the rearrangement of pent-1-en-4-yn-3-one to acryloylvinylidene, which then cyclizes to cyclopenta-2,4-dienone. The ACSE has some unique advantages in the treatment of carbene reactions and more general families of reactions in which the importance of multireference correlation is not known a priori: (i) the ACSE is more reliable than single-reference methods for confirming the presence or absence of multireference correlation and (ii) in the absence of multireference correlation, unlike multireference second-order perturbation theory (MRPT2), the ACSE recovers more single-reference correlation energy than similarly scaling coupled-cluster methods. Because MRPT2 does not recover as much single-reference correlation as the coupled-cluster or ACSE methods, it tends to underestimate reaction barriers within the carbene reactions. For example, in the rearrangement of pent-1-en-4-yn-3-one, the ACSE and CCSD(T) methods produce cyclization barriers of 18.9 and 18.7 kcal/mol with the 6-31G(d) basis set, whereas MRPT2 predicts this barrier to be 12.1 kcal/mol; furthermore, both the ACSE and CCSD(T) determine the energy of the transition state for acryloylvinylidene formation to be 6.6-6.7 kcal/mol above that of the carbene, and yet, MRPT2 does not predict a transition state.

Energy Barriers in the Conversion of Bicyclobutane to gauche-1,3-Butadiene from the Anti-Hermitian Contracted Schrödinger Equation

In 1959 Charles Coulson popularized the idea of computing a molecule's ground-state energy as a functional of the two-electron reduced density matrix (2-RDM) without the many-electron wave function. For 50 years, however, a practical, direct calculation of the 2-RDM was not achieved because the 2-RDM must be constrained by N-representability conditions to represent an N-electron system. Recently, two general approaches to the direct calculation of the 2-RDM have emerged including the solution of the anti-Hermitian contracted Schrodinger equation (ACSE) [Mazziotti, Phys. Rev. Lett., 2006, 97, 143002]. In this article, after further extending the theoretical development of the ACSE method for multireference correlation, we apply the ACSE to studying an unresolved question regarding the opening of bicyclobutane to gauche-1,3-butadiene by conrotatory and disrotatory pathways. Previous theoretical values for the disrotatory energy barrier reveal a disagreement between correlation methods on the order of 10 kcal/mol. By capturing significantly more correlation energy than traditional multireference methods, the ACSE provides new insight into this discrepancy. The ACSE energy for the conrotatory energy barrier agrees with the 40.6 +/- 2.5 kcal/mol experimental value.

Determining the Energy Gap between the Cis and Trans Isomers of HO3-Using Geometry Optimization within the Anti-Hermitian Contracted Schrödinger and Coupled Cluster Methods

The cis and trans isomers of the HO3- anion, which are important in proposed mechanisms for ozonization, are studied computationally. Relative energies, geometries, and normal-mode frequencies are calculated with anti-Hermitian contracted Schrödinger equation (ACSE) and coupled cluster methods. Both the ACSE method and the coupled cluster method with single and double excitations (CCSD) are applied in a correlation-consistent polarized double-zeta basis set (cc-pVDZ). Using coupled cluster with singles, doubles, and perturbative triples (CCSD(T)), we treat the problem with larger basis sets than those in previous work, including correlation-consistent polarized quadruple-zeta basis sets with (aug-cc-pVQZ) and without (cc-pVQZ) diffuse functions, which permit extrapolation of the cis and trans energies to the complete-basis-set limit. The cis isomer is found to be lower in energy than the trans isomer by -3.5 kcal/mol, which is 50% larger in magnitude than the best previous result of -2.2 kcal/mol. The bond lengths between the O2 and OH fragments of the cis- and trans-HO3 are calculated to be 1.713 and 1.857 A, respectively, where both bond lengths are significantly longer than the 1.464 A O-O bond in hydrogen peroxide. In this paper, we extend the ACSE method [Mazziotti, D. A. J. Chem. Phys. 2007, 126, 184101], which computes the two-electron reduced density matrix directly, to include geometry optimization by a Newton's method with numerical derivatives. Calculation of the cis- and trans-HO3- isomers by the ACSE yields energies, geometries, and frequencies that are closer to those from CCSD(T) than those from CCSD.

Multireference many-electron correlation energies from two-electron reduced density matrices computed by solving the anti-Hermitian contracted Schrödinger equation

Two-electron reduced density matrices (2-RDMs) have recently been directly calculated by solving the anti-Hermitian contracted Schr\"odinger equation (ACSE) to obtain $95$--$100\phantom{\rule{0.3em}{0ex}}%$ of the ground-state correlation energy of atoms and molecules with the accuracy increasing with the size of the one-electron basis set [Mazziotti, Phys. Rev. Lett. 97, 143002 (2006).] In this paper, the ACSE method is extended to treat strong multireference correlation effects that are often important at nonequilibrium molecular geometries. While previous ACSE calculations have employed an initial 2-RDM from the Hartree-Fock method, we initialize the solution of the ACSE with a 2-RDM guess from a multiconfiguration self-consistent field calculation. Applications are made to multireference correlation in the potential energy surfaces of the molecules $\mathrm{HF}$, ${\mathrm{H}}_{2}\mathrm{O}$, and ${\mathrm{C}}_{2}$ in polarized valence double-zeta basis sets while ${\mathrm{N}}_{2}$ is treated in polarized valence double- and triple-zeta basis sets. Accurate ground-state energies and 1-RDM occupation numbers are obtained at both equilibrium and nonequilibrium geometries where the energies are within a few millihartrees of those from full configuration interaction.

Cumulant reconstruction of the three-electron reduced density matrix in the anti-Hermitian contracted Schrödinger equation

Differing perspectives on the accuracy of three-electron reduced-density-matrix (3-RDM) reconstruction in nonminimal basis sets exist in the literature. This paper demonstrates the accuracy of cumulant-based reconstructions, developed by Valdemoro (V) [F. Colmenero et al., Phys. Rev. A 47, 971 (1993)], Nakatsuji and Yasuda (NY) [Phys. Rev. Lett. 76, 1039 (1996)], Mazziotti (M) [Phys. Rev. A 60, 3618 (1999)], and Valdemoro-Tel-Perez-Romero (VTP) [Many-electron Densities and Density Matrices, edited by J. Cioslowski (Kluwer, Boston, 2000)]. Computationally, we extend previous investigations to study a variety of molecules, including LiH, HF, NH(3), H(2)O, and N(2), in Slater-type, double-zeta, and polarized double-zeta basis sets at both equilibrium and nonequilibrium geometries. The reconstructed 3-RDMs, compared with 3-RDMs from full configuration interaction, demonstrate in nonminimal basis sets the accuracy of the first-order expansion (V) as well as the important role of the second-order corrections (NY, M, and VTP). Calculations at nonequilibrium geometries further show that cumulant functionals can reconstruct the 3-RDM from a multireferenced 2-RDM with reasonable accuracy, which is relevant to recent multireferenced formulations of the anti-Hermitian contracted Schrodinger equation (ACSE) and canonical diagonalization. Theoretically, we perform a detailed perturbative analysis of the M functional to identify its second-order components. With these second-order components we connect the M, NY, and VTP reconstructions for the first time by deriving both the NY and VTP functionals from the M functional. Finally, these 3-RDM reconstructions are employed within the ACSE [D. Mazziotti, Phys. Rev. Lett. 97, 143002 (2006)] to compute ground-state energies which are compared with the energies from the contracted Schrodinger equation and several wave function methods.