We propose an algorithm to obtain the ground-state energy of a many-electron system using the variational wave function of a linear combination of antisymmetrized geminal powers. We optimized this algorithm to obtain the energy and the other parameters of a many-electron system. Also we clarified the bottleneck of the total calculation in the tensor contraction and successfully reduced the computational time. As a result, we can use an extended number of geminal states to obtain the ground state of the water molecule and Hubbard models. The result for the water molecule with the Dunning double-zeta basis is of the sub-milihartree order above the energy of exact diagonalization. Further, we observe that the result for the one-dimensional Hubbard model with 14 sites shows good tendency to capture the right ground state and that for the two-dimensional Hubbard model still lacks some part of the energy reflecting the large size of the Hilbert space. We conclude that the required number of terms for geminal states for sufficiently accurate energy is only moderately affected by the problem size. We further show other technical details for the numerical algorithms of geminal states in the variation process. It is expected that with the use of more extended computing resources and larger sizes of electronic systems, our algorithm can provide improved results.
We propose an algorithm to obtain the ground-state energy of a many-electron system using the variational wave function of a linear combination of antisymmetrized geminal powers. We optimized this algorithm to obtain the energy and the other parameters of a many-electron system. Also we clarified the bottleneck of the total calculation in the tensor contraction and successfully reduced the computational time. As a result, we can use an extended number of geminal states to obtain the ground state of the water molecule and Hubbard models. The result for the water molecule with the Dunning double-zeta basis is of the sub-milihartree order above the energy of exact diagonalization. Further, we observe that the result for the one-dimensional Hubbard model with 14 sites shows good tendency to capture the right ground state and that for the two-dimensional Hubbard model still lacks some part of the energy reflecting the large size of the Hilbert space. We conclude that the required number of terms for geminal states for sufficiently accurate energy is only moderately affected by the problem size. We further show other technical details for the numerical algorithms of geminal states in the variation process. It is expected that with the use of more extended computing resources and larger sizes of electronic systems, our algorithm can provide improved results.
In previous research, we tested the wave-function format of a linear combination of several antisymmetrized geminal power states. A numerical problem in the geminal matrices was noted, which made the total energies of electronic systems with large numbers of electrons unstable. The underlying cause was found to be the large cancellation term in the geminal power series. We have obtained a format to resolve this problem for the case of total energies and partly for the first-order derivatives within the antisymmetrized geminal power states. By using this formalism, we have calculated the ground-state energies for several electronic systems, including the usage of a larger chemical basis set. The results are, in some cases, very close to the exact result, especially for one-dimensional Hubbard systems. Our result for a water molecule with the Dunning double zeta basis set is better than the configuration-interaction singles and doubles energy and approaches the coupled cluster singles and doubles energy.
To avoid the combinatorial computational cost of configuration interaction (CI), we previously introduced the symmetric tensor decomposition CI (STD-CI) method, which takes advantage of the antisymmetric nature of the electronic wave function and expresses the CI coefficients compactly as a series of Kronecker product states (STD series) [W. Uemura and O. Sugino, Phys. Rev. Lett. 109, 253001 (2012)]. Here we extend the variational degrees of freedom by using different molecular orbitals for different terms in the STD series. This scheme is equivalent to the linear combination of the Hartree-Fock-Bogoliubov state or the antisymmetrized geminal powers (AGPs). The total energy converges very rapidly within 0.72 $\ensuremath{\mu}\mathrm{hartree}$ taking only 10 terms for the water molecule, and the convergence is likewise fast for Hubbard tetramers. The computational cost scales as the fifth power of the number of electrons and the square of the number of terms in the STD series, indicating the promise of this AGP-based scheme for highly accurate and efficient computation of quantum systems.
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.
The energy of a molecule is not the sum of the energies of its atomic components. Our aim is to compute an accurate approximation of this quantity based on an electron pair model, that is to say by using an antisymmetric product of two-electron wave functions, called “geminals”. In this model, the total wave function can also be described in a more practical way, by a set of matrices, one for each geminal. However, without further restrictions, such a model has a factorial computational complexity with the number of electrons and its applicability is therefore limited to small systems.We have introduced generalized orthogonality constraints between geminals to reduce the computational effort to a polynomial cost, without sacrificing the indistinguishability of the electrons. In other words, the electron pairs corresponding to the geminals are not fully distinguishable, so that our model remains complex, in the sense that geminal products have still to be antisymmetrized according to the Pauli principle to form bona fide electronic wave functions.Our geometrical constraints translate, in terms of geminal matrices, into simple equations involving the traces of products of these matrices. In the simplest non-trivial model, a set of solutions is given by block-diagonal matrices where each block is of size 2x2 and consists of a Pauli matrix or the Identity, multiplied by a complex parameter to be optimized. With this simplified ansatz for geminals, the number of terms in the calculation of the matrix elements of quantum observables is considerably reduced.
Slater determinants have underpinned quantum chemistry for nearly a century, yet their full potential has remained challenging to exploit. In this work, we show that a variational wavefunction composed of a few hundred optimized non-orthogonal determinants can achieve energy accuracies comparable to the state of the art. This is obtained by introducing an optimization method that leverages the quadratic dependence of the variational energy on the orbitals of each determinant, enabling an exact iterative optimization, and uses an efficient tensor-contraction algorithm to evaluate the effective Hamiltonian with a computational cost that scales as the fourth power of the number of basis functions. We benchmark the accuracy of the proposed method with exact full-configuration interaction results where available, and we achieve lower variational energies than coupled cluster (CCSD(T)) for several molecules in the double-zeta basis. Accurate description of electronic wave functions typically requires large numbers of Slater determinants. Here, the authors develop an optimization method that achieves state-of-the-art accuracy for small molecules with only few hundred determinants.
Variational minimization of the ground-state energy as a function of the two-electron reduced density matrix (2-RDM), constrained by necessary N-representability conditions, provides a polynomial-scaling approach to studying strongly correlated molecules without computing the many-electron wave function. Here we introduce a route to enhancing necessary conditions for N representability through rank restriction of the 2-RDM. Rather than adding computationally more expensive N-representability conditions, we directly enhance the accuracy of two-particle (2-positivity) conditions through rank restriction, which removes degrees of freedom in the 2-RDM that are not sufficiently constrained. We select the rank of the particle-hole 2-RDM by deriving the ranks associated with model wave functions, including both mean-field and antisymmetrized geminal power (AGP) wave functions. Because the 2-positivity conditions are exact for quantum systems with AGP ground states, the rank of the particle-hole 2-RDM from the AGP ansatz provides a minimum for its value in variational 2-RDM calculations of general quantum systems. To implement the rank-restricted conditions, we extend a first-order algorithm for large-scale semidefinite programming. The rank-restricted conditions significantly improve the accuracy of the energies; for example, the percentages of correlation energies recovered for HF, CO, and N2 improve from 115.2%, 121.7%, and 121.5% without rank restriction to 97.8%, 101.1%, and 100.0% with rank restriction. Similar results are found at both equilibrium and nonequilibrium geometries. While more accurate, the rank-restricted N-representability conditions are less expensive computationally than the full-rank conditions.