2-electron Reduced Density Matrix 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.

Assessing the orbital-optimized unitary Ansatz for density cumulant theory.

The previously proposed Ansatz for density cumulant theory that combines orbital-optimization and a parameterization of the 2-electron reduced density matrix cumulant in terms of unitary coupled cluster amplitudes (OUDCT) is carefully examined. Formally, we elucidate the relationship between OUDCT and orbital-optimized unitary coupled cluster theory and show the existence of near-zero denominators in the stationarity conditions for both the exact and some approximate OUDCT methods. We implement methods of the OUDCT Ansatz restricted to double excitations for numerical study, up to the fifth commutator in the Baker-Campbell-Hausdorff expansion. We find that methods derived from the Ansatz beyond the previously known ODC-12 method tend to be less accurate for equilibrium properties and less reliable when attempting to describe H2 dissociation. New developments are needed to formulate more accurate density cumulant theory variants.

Active-Space Pair Two-Electron Reduced Density Matrix Theory for Strong Correlation

An active-space variational calculation of the two-electron reduced density matrix (2-RDM) is derived and implemented where the active orbitals are correlated within the pair approximation. The pair approximation considers only doubly occupied configurations of the wave function, which enables the calculation of the 2-RDM at a computational cost of [Formula: see text]. Calculations were performed both with the pair active-space configuration interaction (PASCI) method and the pair active-space self-consistent field (PASSCF) method. The latter includes a mixing of the active and inactive orbitals through unitary transformations. The active-space pair 2-RDM method is applied to the nitrogen molecule, the p-benzyne diradical, a newly synthesized biscobalt complex, and the nitrogenase cofactor FeMoco. The FeMoco molecule is treated in a [120,120] active space. Fractional occupations are recovered in each of these systems, indicating the presence of strong electron correlation.

Redox, transmetalation, and stacking properties of tetrathiafulvalene-2,3,6,7-tetrathiolate bridged tin, nickel, and palladium compounds.

Here we report that capping the molecule TTFtt (TTFtt = tetrathiafulvalene-2,3,6,7-tetrathiolate) with dialkyl tin groups enables the isolation of a stable series of redox congeners and facile transmetalation to Ni and Pd. TTFtt has been proposed as an attractive building block for molecular materials for two decades as it combines the redox chemistry of TTF and dithiolene units. TTFttH4, however, is inherently unstable and the incorporation of TTFtt units into complexes or materials typically proceeds through the in situ generation of the tetraanion TTFtt4−. Capping of TTFtt4− with Bu2Sn2+ units dramatically improves the stability of the TTFtt moiety and furthermore enables the isolation of a redox series where the TTF core carries the formal charges of 0, +1, and +2. All of these redox congeners show efficient and clean transmetalation to Ni and Pd resulting in an analogous series of bimetallic complexes capped by 1,2-bis(diphenylphosphino)ethane (dppe) ligands. Furthermore, by using the same transmetalation method, we synthesized analogous palladium complexes capped by 1,1′-bis(diphenylphosphino)ferrocene (dppf) which had been previously reported. All of these species have been thoroughly characterized through a systematic survey of chemical and electronic properties by techniques including cyclic voltammetry (CV), ultraviolet-visible-near infrared spectroscopy (UV-vis-NIR), electron paramagnetic resonance spectroscopy (EPR), nuclear magnetic resonance spectroscopy (NMR) and X-ray diffraction (XRD). These detailed synthetic and spectroscopic studies highlight important differences between the transmetalation strategy presented here and previously reported synthetic methods for the installation of TTFtt. In addition, the utility of this stabilization strategy can be illustrated by the observation of unusual TTF radical–radical packing in the solid state and dimerization in the solution state. Theoretical calculations based on variational 2-electron reduced density matrix methods have been used to investigate these unusual interactions and illustrate fundamentally different levels of covalency and overlap depending on the orientations of the TTF cores. Taken together, this work demonstrates that tin-capped TTFtt units are ideal reagents for the installation of redox-tunable TTFtt ligands enabling the generation of entirely new geometric and electronic structures.

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.

Method for making 2-electron response reduced density matrices approximately N-representable.

In methods like geminal-based approaches or coupled cluster that are solved using the projected Schrödinger equation, direct computation of the 2-electron reduced density matrix (2-RDM) is impractical and one falls back to a 2-RDM based on response theory. However, the 2-RDMs from response theory are not N-representable. That is, the response 2-RDM does not correspond to an actual physical N-electron wave function. We present a new algorithm for making these non-N-representable 2-RDMs approximately N-representable, i.e., it has the right symmetry and normalization and it fulfills the P-, Q-, and G-conditions. Next to an algorithm which can be applied to any 2-RDM, we have also developed a 2-RDM optimization procedure specifically for seniority-zero 2-RDMs. We aim to find the 2-RDM with the right properties which is the closest (in the sense of the Frobenius norm) to the non-N-representable 2-RDM by minimizing the square norm of the difference between this initial response 2-RDM and the targeted 2-RDM under the constraint that the trace is normalized and the 2-RDM, Q-matrix, and G-matrix are positive semidefinite, i.e., their eigenvalues are non-negative. Our method is suitable for fixing non-N-representable 2-RDMs which are close to being N-representable. Through the N-representability optimization algorithm we add a small correction to the initial 2-RDM such that it fulfills the most important N-representability conditions.

Development and application of a 2-electron reduced density matrix approach to electron transport via molecular junctions.

Can an electronic device be constructed using only a single molecule? Since this question was first asked by Aviram and Ratner in the 1970s [Chem. Phys. Lett. 29, 277 (1974)], the field of molecular electronics has exploded with significant experimental advancements in the understanding of the charge transport properties of single molecule devices. Efforts to explain the results of these experiments and identify promising new candidate molecules for molecular devices have led to the development of numerous new theoretical methods including the current standard theoretical approach for studying single molecule charge transport, i.e., the non-equilibrium Green's function formalism (NEGF). By pairing this formalism with density functional theory (DFT), a wide variety of transport problems in molecular junctions have been successfully treated. For some systems though, the conductance and current-voltage curves predicted by common DFT functionals can be several orders of magnitude above experimental results. In addition, since density functional theory relies on approximations to the exact exchange-correlation functional, the predicted transport properties can show significant variation depending on the functional chosen. As a first step to addressing this issue, the authors have replaced density functional theory in the NEGF formalism with a 2-electron reduced density matrix (2-RDM) method, creating a new approach known as the NEGF-RDM method. 2-RDM methods provide a more accurate description of electron correlation compared to density functional theory, and they have lower computational scaling compared to wavefunction based methods of similar accuracy. Additionally, 2-RDM methods are capable of capturing static electron correlation which is untreatable by existing NEGF-DFT methods. When studying dithiol alkane chains and dithiol benzene in model junctions, the authors found that the NEGF-RDM predicts conductances and currents that are 1-2 orders of magnitude below those of B3LYP and M06 DFT functionals. This suggests that the NEGF-RDM method could be a viable alternative to NEGF-DFT for molecular junction calculations.

Pair 2-electron reduced density matrix theory using localized orbitals.

Full configuration interaction (FCI) restricted to a pairing space yields size-extensive correlation energies but its cost scales exponentially with molecular size. Restricting the variational two-electron reduced-density-matrix (2-RDM) method to represent the same pairing space yields an accurate lower bound to the pair FCI energy at a mean-field-like computational scaling of O(r3) where r is the number of orbitals. In this paper, we show that localized molecular orbitals can be employed to generate an efficient, approximately size-extensive pair 2-RDM method. The use of localized orbitals eliminates the substantial cost of optimizing iteratively the orbitals defining the pairing space without compromising accuracy. In contrast to the localized orbitals, the use of canonical Hartree-Fock molecular orbitals is shown to be both inaccurate and non-size-extensive. The pair 2-RDM has the flexibility to describe the spectra of one-electron RDM occupation numbers from all quantum states that are invariant to time-reversal symmetry. Applications are made to hydrogen chains and their dissociation, n-acene from naphthalene through octacene, and cadmium telluride 2-, 3-, and 4-unit polymers. For the hydrogen chains, the pair 2-RDM method recovers the majority of the energy obtained from similar calculations that iteratively optimize the orbitals. The localized-orbital pair 2-RDM method with its mean-field-like computational scaling and its ability to describe multi-reference correlation has important applications to a range of strongly correlated phenomena in chemistry and physics.

Ligand non-innocence and strong correlation in manganese superoxide dismutase mimics.

We examine the 1-electron reduction of manganese porphyrin complexes Mn(iii) porphyrin and Mn(iii) TDE-2-ImP5+, which have attracted recent interest due to their properties as superoxide dismutase mimics. We perform a series of electronic structure calculations using the variational 2-electron reduced density matrix (2-RDM) method with a large [30,30] active space that represents a wavefunction with 1019 variables, as well as the more traditional complete active space self-consistent field (CASSCF) method with a [14,14] active space. We show that the larger 2-RDM calculation, intractable with CASSCF, is required to capture the full effects of electron correlation in the molecule and predict the non-innocence of the porphyrin ligand during the reduction. The CASSCF method predicts single-reference systems exhibiting a metal-centered reduction, but the 2-RDM method predicts a strongly correlated system exhibiting a ligand-centered reduction. Based on these results, we find that the porphyrin ligand is reduced rather than the manganese, and suggest that the electron correlation plays a role in driving the ligand non-innocence.

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.

The tensor hypercontracted parametric reduced density matrix algorithm: Coupled-cluster accuracy with O(r4) scaling

Tensor hypercontraction is a method that allows the representation of a high-rank tensor as a product of lower-rank tensors. In this paper, we show how tensor hypercontraction can be applied to both the electron repulsion integral tensor and the two-particle excitation amplitudes used in the parametric 2-electron reduced density matrix (p2RDM) algorithm. Because only O(r) auxiliary functions are needed in both of these approximations, our overall algorithm can be shown to scale as O(r4), where r is the number of single-particle basis functions. We apply our algorithm to several small molecules, hydrogen chains, and alkanes to demonstrate its low formal scaling and practical utility. Provided we use enough auxiliary functions, we obtain accuracy similar to that of the standard p2RDM algorithm, somewhere between that of CCSD and CCSD(T).

Comparison of low-rank tensor expansions for the acceleration of quantum chemistry computations

Low-rank spectral expansion and tensor hypercontraction are two promising techniques for reducing the size of the two-electron excitation tensor by factorizing it into products of smaller tensors. Both methods can potentially realize an O(r(4)) quantum chemistry method where r is the number of one-electron orbitals. We compare the two factorizations in this paper by applying them to the parametric 2-electron reduced density matrix method with the M functional [D. A. Mazziotti, Phys. Rev. Lett. 101, 253002 (2008)]. We study several inorganic molecules, alkane chains, and potential curves as well as reaction and dissociation energies. The low-rank spectral expansion, we find, is typically more efficient than tensor hypercontraction due to a faster convergence of the energy and a smaller constant prefactor in the energy optimization. Both factorizations are applicable to the acceleration of a wide range of wavefunction and reduced-density-matrix methods.

Relative Energies and Geometries of the cis- and trans-HO3 Radicals from the Parametric 2-Electron Density Matrix Method

The parametric 2-electron reduced density matrix (2-RDM) method employing the M functional [Mazziotti, D. A. Phys. Rev. Lett. 2008, 101, 253002], also known as the 2-RDM(M) method, improves on the accuracy of coupled electron-pair theories including coupled cluster with single-double excitations at the computational cost of configuration interaction with single-double excitations. The cis- and trans-HO(3) isomers along with their isomerization transition state were examined using the recent extension of 2-RDM(M) to nonsinglet open-shell states [Schwerdtfeger, C. A.; Mazziotti, D. A. J. Chem. Phys. 2012, 137, 034107] and several coupled cluster methods. We report the calculated energies, geometries, natural-orbital occupation numbers, and reaction barriers for the HO(3) isomers. We find that the 2-RDM(M) method predicts that the trans isomer of HO(3) is lower in energy than the cis isomer by 1.71 kcal/mol in the correlation-consistent polarized valence quadruple-ζ (cc-pVQZ) basis set and 1.84 kcal/mol in the augmented correlation-consistent polarized valence quadruple-ζ (aug-cc-pVQZ) basis set. Results include the harmonic zero-point vibrational energies calculated in the correlation-consistent polarized valence double-ζ basis set. On the basis of the results of a geometry optimization in the augmented correlation consistent polarized valence triple-ζ basis set, the parametric 2-RDM(M) method predicts a central oxygen-oxygen bond of 1.6187 Å. We compare these energies and geometries to those predicted by three single-reference coupled cluster methods and experimental results and find that the inclusion of multireference correlation is important to describe properly the relative energies of the cis- and trans-HO(3) isomers and improve agreement with experimental geometries.

Populations of Carbonic Acid Isomers at 210 K from a Fast Two-Electron Reduced-Density Matrix Theory

Parametrization of the 2-electron reduced density matrix (2-RDM) rather than the many-electron wave function yields a new family of electronic-structure methods that are faster and more accurate than traditional coupled electron-pair methods including coupled cluster with single and double excitations. Deriving the parametrization from N-representability conditions generates a 2-RDM that captures significant correlation from triple and higher-order excitations at the cost of double excitations. We apply the parametric 2-RDM method to confirm recent experiments determining the relative thermodynamic populations of the cis-cis and cis-trans isomers of carbonic acid. In 2010 Bernard et al. showed by infrared spectroscopy that the populations of cis-cis and cis-trans isomers have a 10:1 ratio at 210 K. By use of the parametric 2-RDM method, we predict a 8:1 ratio at 210 K. Comparable ab initio methods overestimate the stability of the cis-cis isomer with 24:1 and 21:1 ratios. These 2-RDM-based methods promise to have significant applications throughout chemistry.