We propose an approach to the electronic structure problem based on noninteracting electron pairs that has similar computational cost to conventional methods based on noninteracting electrons. In stark contrast to other approaches, the wave function is an antisymmetric product of nonorthogonal geminals, but the geminals are structured so the projected Schrödinger equation can be solved very efficiently. We focus on an approach where, in each geminal, only one of the orbitals in a reference Slater determinant is occupied. The resulting method gives good results for atoms and small molecules. It also performs well for a prototypical example of strongly correlated electronic systems, the hydrogen atom chain.
We propose an approach to the electronic structure problem based on noninteracting electron pairs that has similar computational cost to conventional methods based on noninteracting electrons. In stark contrast to other approaches, the wave function is an antisymmetric product of nonorthogonal geminals, but the geminals are structured so the projected Schrödinger equation can be solved very efficiently. We focus on an approach where, in each geminal, only one of the orbitals in a reference Slater determinant is occupied. The resulting method gives good results for atoms and small molecules. It also performs well for a prototypical example of strongly correlated electronic systems, the hydrogen atom chain.
The Lewis electron-pair picture pervades all of chemistry and a substantialdomain of materials science. “Geminals” are electron-pair quantum states whichare used to express approximate solutions of the Schrödinger equation of polyelec-tronic systems. In quantum chemistry, geminal-based methods are potentially moreeffective than the traditional ones based on the orbital picture: electrons occupyorbitals and experience only the average effect of each other. This is especially truefor strongly correlated-systems. However, the computational cost of using arbitrarygeminal wave functions is non-polynomial in the system size.It is therefore an important challenge to find additional constraints on geminalmodels which preserve the accuracy of the wavefunction description while reducingthe calculations to polynomial cost. The Ayers’ group at Mc Master has developedone such successful geminal model called AP1roG because the wavefunction is anantisymmetric product of geminals with one distinct reference orbital occupied ineach geminal [1]. This ansatz relaxes the so-called “strong-orthogonality” constraintbetween geminals. The latter makes wavefunction matrix element computationsvery easy but appears too drastic from the physical point of view, since the geminalelectron pairs are then distinguishable.We have followed a different path to relax the strong-orthogonality constraint. Wehave defined the geometrical concept of n-orthogonality [2] which can be viewed asa graded indistinguishability measure for electronic states: 1-orthogonality coincideswith strong orthogonality, the larger n the less distinguishable the n-orthogonal elec-tronic states will be. We have studied geminal models constrained by n-orthogonalityrelations for different values of n, imposed to different combinations of geminal prod-uct states [3]. However, in this work the scaling of matrix element computations wasstill unsatisfactory. We will present new ideas developed since then to remedy thisproblem.References[1] P. A Limacher, P. W. Ayers, P. A. Johnson, S. De Baerdemacker, D. Van Neck and P. Bultinck, J. of Chem. Theory and Computation 9, 1394-1401 (2013).[2] P. Cassam-Chenaı̈, Phys. Rev. A77, 032103 (2008).[3] P. Cassam-Chenaı̈, V. Rassolov, Chem. Phys. Lett. 487, 147-152 (2010).
Wave functions constructed from electron-pair states can accurately model strong electron correlation effects and are promising approaches especially for larger many-body systems. In this article, we analyze the nature and the type of electron correlation effects that can be captured by wave functions restricted to electron-pair states. We focus on the pair-coupled-cluster doubles (pCCD) ansatz also called the antisymmetric product of the 1-reference orbital geminal (AP1roG) method, combined with an orbital optimization protocol presented in Boguslawski et al. [Phys. Rev. B 89, 201106(R) (2014)], whose performance is assessed against electronic structures obtained form density-matrix renormalization-group reference data. Our numerical analysis covers model systems for strong correlation: the one-dimensional Hubbard model with a periodic boundary condition as well as metallic and molecular hydrogen rings. Specifically, the accuracy of pCCD/AP1roG is benchmarked using the single-orbital entropy, the orbital-pair mutual information, as well as the eigenvalue spectrum of the one-orbital and two-orbital reduced density matrices. Our study indicates that contributions from singly occupied states become important in the strong correlation regime which highlights the limitations of the pCCD/AP1roG method. Furthermore, we examine the effect of orbital rotations within the pCCD/AP1roG model on correlations between orbital pairs.
In this review, we discuss the recent progress in developing geminal-based theories for challenging problems in quantum chemistry. Specifically, we focus on the antisymmetrized geminal power, generalized valence bond, antisymmetrized product of strongly orthogonal geminals, singlet-type orthogonal geminals, the antisymmetric product of 1-reference orbital geminal, also known as the pair coupled cluster doubles ansatz, and geminals constructed from Richardson-Gaudin states. Furthermore, we review various corrections to account for the missing dynamical correlation effects in geminal models and possible extensions to target electronically excited states and open-shell species. Finally, we discuss some numerical examples and present-day challenges for geminal-based models, including a quantitative and qualitative analysis of wave functions, and software availability.
Strong electronic correlations are at the heart of many interesting phenomena. For the theoretical description of these materials, a proper treatment of the local atomic physics is required. We propose a novel approach combining functionals of the electron density and functionals of the one-particle reduced density matrix to improve this description. The proposed method has a solid foundation in reduced density-matrix functional theory and no double-counting problem arises. It employs a decomposition of the electron-electron interaction in real space. The interaction close to the correlated orbitals, for example, the partially filled 3d-orbitals of transition-metal ions, is described with a density-matrix functional and otherwise with a local or semi-local density functional. We propose to evaluate the density-matrix functional from Levy’s constrained search problem, i.e., via a constrained minimization over an ensemble of many-particle wave functions. In contrast to approximate parametrized density-matrix functionals, this evaluation allows us to systematically improve the functional towards the exact result. In situations where the one-particle basis is too large to evaluate the density-matrix functional from Levy’s constrained minimization problem, we apply a series of approximations that each can be converged to the exact result: the first approximation step is the local approximation of the density-matrix functional proposed by Blöchl, Walther and Pruschke. For the density-matrix functionals within the local interactions, we propose the adaptive cluster approximation (ACA) that systematically truncates non-interacting one-particle states and drastically reduces the computational effort. The resulting density-matrix functional for a local interaction and a small number of non-interacting one-particle states is then evaluated with the constrained minimization problem. We explore different parametrizations of many-particle wave functions in this constrained minimization problem. A parametrization based on a configuration-interaction-like ansatz is shown to converge rapidly if suitable selection criterion for the Slater determinants is chosen. An impurity-bath-separation ansatz is shown to be suitable for single-impurity Anderson models. It is shown that the constrained minimization can be solved for matrix product states with a DMRG-like iterative minimization. Furthermore, we show that Gutzwiller-Jastrow-correlated wave functions can be used with a quantum Monte Carlo procedure as many-particle wave functions. Finally, we formulate an algorithm for the evaluation of the density-matrix functional on near-term quantum computers. Results from the execution of the algorithm on an existing quantum computer with transmon qubits are presented. The proposed approach combining density functionals and density-matrix functionals is implemented in the CP-PAW code based on the projector augmented-wave formalism. We present results for the dissociation curve of the hydrogen molecule as the prototypical case of strong static correlation and show that static correlation is well described with the new approach. We show results for the nonmagnetic state of the transition-metaoxide NiO that is described qualitatively wrong with the DFT+U method. The proposed method properly describes the nonmagnetic state of NiO as an insulator and predicts a qualitatively correct spectral function.
We introduce an efficient method to construct optimal and system adaptive basis sets for use in electronic structure and quantum Monte Carlo calculations. The method is based on an embedding scheme in which a reference atom is singled out from its environment, while the entire system (atom and environment) is described by a Slater determinant or its antisymmetrized geminal power (AGP) extension. The embedding procedure described here allows for the systematic and consistent contraction of the primitive basis set into geminal embedded orbitals (GEOs), with a dramatic reduction of the number of variational parameters necessary to represent the many-body wave function, for a chosen target accuracy. Within the variational Monte Carlo method, the Slater or AGP part is determined by a variational minimization of the energy of the whole system in presence of a flexible and accurate Jastrow factor, representing most of the dynamical electronic correlation. The resulting GEO basis set opens the way for a fully controlled optimization of many-body wave functions in electronic structure calculation of bulk materials, namely, containing a large number of electrons and atoms. We present applications on the water molecule, the volume collapse transition in cerium, and the high-pressure liquid hydrogen.
We investigated some variational methods to compute a wavefunction based on antisymmetric product of geminals (APG). The Waring decomposition on the APG wavefunction leads a finite sum of antisymmetrized geminal power (AGP) wavefunctions, each for which the variational principle can be applied. We call this as AGP-CI method which provides a variational solution of the APG wavefunction efficiently. However, number of AGP wavefunctions in the exact AGP-CI formalism become exponentially large in case of many-electron systems. Therefore, we also investigate the low-rank APG wavefunction, in which the coefficient matrices of geminals are factorized by the Schur decomposition. Interestingly, only a few non-zero eigenvalues (up to half number of electrons) were found from the Schur decomposition on the APG wavefunction. We developed some methods to approximate the APG wavefunction by lowering the ranks of coefficient matrices of geminals, and demonstrate their performance. Our new geminal method based on the low-rank decomposition can drastically reduce the number of variational parameters, although no efficient formalism has been elaborated so far, due to the difficulty of optimization caused by the inability to compute analytical derivatives.