Variational reduced density matrix method in the doubly occupied configuration interaction space using three-particleN-representability conditions

Determination of electronic excitation energies within the doubly occupied configuration interaction space by means of the Hermitian operator method

In this work, we formulate the equations of motion corresponding to the Hermitian operator method in the framework of the doubly occupied configuration interaction space. The resulting algorithms turn out to be considerably simpler than the equations provided by that method in more conventional spaces, enabling the determination of excitation energies in N-electron systems under an affordable polynomial computational cost. The implementation of this technique only requires to know the elements of low-order reduced density matrices of an N-electron reference state, which can be obtained from any approximate method. We contrast our procedure against the reduced Bardeen-Cooper-Schrieffer and Richardson-Gaudin-Kitaev integrable models, pointing out the reliability of our proposal.

Variational determination of the two-electron reduced density matrix within the doubly occupied configuration interaction framework: Treatments of triplet N-electron systems

This work performs variational determinations of two-electron reduced density matrices corresponding to open-shell N-electron systems within the framework of the doubly occupied configuration interaction treatment, traditionally limited to studies of closed-shell systems. The procedure has allowed us to describe satisfactorily molecular systems in triplet states following two methods. One of them adds hydrogen atoms at an infinite distance of the triplet system studied, constituting a singlet supersystem. Energies and reduced density matrices of the triplet system are obtained by removing the contributions of the added atoms from the singlet supersystem results. The second procedure determines variationally the two-electron reduced density matrices corresponding to the triplet systems by means of adequate couplings of basis-set functions. Both models have been managed by imposing N-representability conditions on the reduced density matrix calculations. Results obtained from these methods for molecular systems in triplet ground states are reported and compared with those provided by benchmark methods.

Variational determination of ground and excited-state two-electron reduced density matrices in the doubly occupied configuration space: A dispersion operator approach

This work implements a variational determination of the elements of two-electron reduced density matrices corresponding to the ground and excited states of N-electron interacting systems based on the dispersion operator technique. The procedure extends the previously reported proposal [Nakata et al., J. Chem. Phys. 125, 244109 (2006)] to two-particle interaction Hamiltonians and N-representability conditions for the two-, three-, and four-particle reduced density matrices in the doubly occupied configuration interaction space. The treatment has been applied to describe electronic spectra using two benchmark exactly solvable pairing models: reduced Bardeen-Cooper-Schrieffer and Richardson-Gaudin-Kitaev Hamiltonians. The dispersion operator combined with N-representability conditions up to the four-particle reduced density matrices provides excellent results.

Quantum Information and Algorithms for Correlated Quantum Matter

Discoveries in quantum materials, which are characterized by the strongly quantum-mechanical nature of electrons and atoms, have revealed exotic properties that arise from correlations. It is the promise of quantum materials for quantum information science superimposed with the potential of new computational quantum algorithms to discover new quantum materials that inspires this Review. We anticipate that quantum materials to be discovered and developed in the next years will transform the areas of quantum information processing including communication, storage, and computing. Simultaneously, efforts toward developing new quantum algorithmic approaches for quantum simulation and advanced calculation methods for many-body quantum systems enable major advances toward functional quantum materials and their deployment. The advent of quantum computing brings new possibilities for eliminating the exponential complexity that has stymied simulation of correlated quantum systems on high-performance classical computers. Here, we review new algorithms and computational approaches to predict and understand the behavior of correlated quantum matter. The strongly interdisciplinary nature of the topics covered necessitates a common language to integrate ideas from these fields. We aim to provide this common language while weaving together fields across electronic structure theory, quantum electrodynamics, algorithm design, and open quantum systems. Our Review is timely in presenting the state-of-the-art in the field toward algorithms with nonexponential complexity for correlated quantum matter with applications in grand-challenge problems. Looking to the future, at the intersection of quantum information science and algorithms for correlated quantum matter, we envision seminal advances in predicting many-body quantum states and describing excitonic quantum matter and large-scale entangled states, a better understanding of high-temperature superconductivity, and quantifying open quantum system dynamics.

An adiabatic connection for doubly-occupied configuration interaction wave functions

An adiabatic connection (AC) is developed as an electron correlation correction for doubly occupied configuration interaction (DOCI) wave functions. Following the work of Pernal [Phys. Rev. Lett. 120, 013001 (2018)], the working equations of the approach, termed AC-DOCI, are rooted in the extended random phase approximation (ERPA) and require knowledge of only the ground-state two-electron reduced density matrix (2RDM) from the DOCI. As such, the AC is naturally suited to modeling electron correlation in variational 2RDM (v2RDM)-based approximations to the DOCI. The v2RDM-driven AC-DOCI is applied to the dissociation of molecular nitrogen and the double dissociation of water; the approach yields energies that are similar in quality to those from second-order multireference perturbation theory near equilibrium, but the quality of the AC-DOCI energy degrades at stretched geometries. The exact adiabatic connection path suggests the assumption that the one-electron reduced-density matrix is constant along the AC path is invalid at stretched geometries, but this deficiency alone cannot explain the observed behavior. Rather, it appears that the ERPA's single-particle-transition ansatz cannot, in general, provide good approximations to the 2RDM along the AC path. The AC-DOCI is also applied to a set of 45 reaction energies; for these systems, the approach has an average accuracy that is comparable to that of single-reference second-order many-body perturbation theory.

Variational reduced density matrix method in the doubly-occupied configuration interaction space using four-particle N-representability conditions: Application to the XXZ model of quantum magnetism

This work deals with the variational determination of the two-particle reduced density matrix (2-RDM) and the energy corresponding to the ground state of N-particle systems within the doubly occupied configuration interaction (DOCI) space. Here, we impose for the first time up to four-particle N-representability constraint conditions in the variational determination of the 2-RDM matrix elements using the standard semidefinite programming algorithms. The energies and 2-RDMs obtained from this treatment and the corresponding computational costs are compared with those arisen from previously reported less restrictive variational methods [D. R. Alcoba et al., J. Chem. Phys. 149, 194105 (2018)] as well as with the exact DOCI values. We apply the different approximations to the one-dimensional XXZ model of quantum magnetism, which has a rich phase diagram with one critical phase and constitutes a stringent test for the method. The numerical results show the usefulness of our treatment to achieve a high degree of accuracy.

Unrestricted treatment for the direct variational determination of the two-electron reduced density matrix for doubly occupied-configuration-interaction wave functions

This work extends to the unrestricted orbital approach the procedure described in our previous report [Alcoba et al., J. Chem. Phys. 148, 024105 (2018)] for determining variationally the two-electron reduced density matrix arising from doubly occupied-configuration-interaction wave functions by imposing two- and three-index N-representability conditions. An analysis of the numerical results obtained in selected systems, from both restricted and unrestricted treatments, allows one to assess the performance of these methodologies as well as to show the influence of the P, Q, G, T1, and T2 positivity conditions. We highlight the satisfactory results obtained within the unrestricted scheme.