arXiv:2012.05283 · Physical Review A

Calculating the distance from an electronic wave function to the manifold of Slater determinants through the geometry of Grassmannians

Yuri Alexandre Aoto, Márcio Fabiano da Silva

Published Nov 6, 2020• 13 citations

Abstract

The set of all electronic states that can be expressed as a single Slater determinant forms a submanifold, isomorphic to the Grassmannian, of the projective Hilbert space of wave functions. We explored this fact by using tools of Riemannian geometry of Grassmannians as described by [Absil et al., Acta Appl. Math. 80, 199 (2004)] to propose an algorithm that converges to a Slater determinant that is a critical point of the overlap function with a correlated wave function. This algorithm can be applied to quantify the entanglement or correlation of a wave function. We show that this algorithm is equivalent to the Newton method using the standard parametrization of Slater determinants by orbital rotations, but it can be more efficiently implemented because the orbital basis used to express the correlated wave function is kept fixed throughout the iterations. We present the equations of this method for a general configuration interaction wave function and for a wave function with up to double excitations over a reference determinant. Applications of this algorithm to selected electronic systems are also presented and discussed.

DOI: 10.1103/physreva.102.052803

Related papers

Precise Quantum Chemistry calculations with few Slater Determinants

Clemens Giuliani, Jannes Nys, Rocco Martinazzo +2 more

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.

arXiv:2503.14502Nature CommunicationsMar 16, 20260 citations

An Optimisation on the Grassmannian with Applications toQuantum Chemistry

Yuri Alexandre Aoto, Márcio Fabiano da Silva

In this work we propose an algorithm to find critical points of the inner product between an element of the Grassmannian and a fixed point of the projective space of the exterior algebra where the Grassmannian is embedded. This has interesting applications to electronic structure theory, where the wave functions are represented by elements of the exterior algebra. This method is exemplified for a Grassmannian that is a model for the hydrogen molecule, H 2 .

Proceeding Series of the Brazilian Society of Computational and Applied MathematicsDec 20, 20211 citation

Influence of single particle orbital sets and configuration selection on multideterminant wavefunctions in quantum Monte Carlo

Raymond C. Clay, Miguel A. Morales

Multideterminant wavefunctions, while having a long history in quantum chemistry, are increasingly being used in highly accurate quantum Monte Carlo calculations. Since the accuracy of QMC is ultimately limited by the quality of the trial wavefunction, multi-Slater determinants wavefunctions offer an attractive alternative to Slater-Jastrow and more sophisticated wavefunction ansatz for several reasons. They can be efficiently calculated, straightforwardly optimized, and systematically improved by increasing the number of included determinants. In spite of their potential, however, the convergence properties of multi-Slater determinant wavefunctions with respect to orbital set choice and excited determinant selection are poorly understood, which hinders the application of these wavefunctions to large systems and solids. In this paper, by performing QMC calculations on the equilibrium and stretched carbon dimer, we find that convergence of the recovered correlation energy with respect to number of determinants can depend quite strongly on basis set and determinant selection methods, especially where there is strong correlation. We demonstrate that properly chosen orbital sets and determinant selection techniques from quantum chemistry methods can dramatically reduce the required number of determinants (and thus the computational cost) to reach a given accuracy, which we argue shows clear need for an automatic QMC-only method for selecting determinants and generating optimal orbital sets.

The Journal of Chemical PhysicsJun 16, 201523 citations

Excited states from quantum Monte Carlo in the basis of Slater determinants

Alexander Humeniuk, Roland Mitrić

Building on the full configuration interaction quantum Monte Carlo (FCIQMC) algorithm introduced recently by Booth et al. [J. Chem. Phys. 131, 054106 (2009)] to compute the ground state of correlated many-electron systems, an extension to the computation of excited states (exFCIQMC) is presented. The Hilbert space is divided into a large part consisting of pure Slater determinants and a much smaller orthogonal part (the size of which is controlled by a cut-off threshold), from which the lowest eigenstates can be removed efficiently. In this way, the quantum Monte Carlo algorithm is restricted to the orthogonal complement of the lower excited states and projects out the next highest excited state. Starting from the ground state, higher excited states can be found one after the other. The Schrödinger equation in imaginary time is solved by the same population dynamics as in the ground state algorithm with modified probabilities and matrix elements, for which working formulae are provided. As a proof of principle, the method is applied to lithium hydride in the 3-21G basis set and to the helium dimer in the aug-cc-pVDZ basis set. It is shown to give the correct electronic structure for all bond lengths. Much more testing will be required before the applicability of this method to electron correlation problems of interesting size can be assessed.

The Journal of Chemical PhysicsNov 17, 201417 citations

Matrix Product State representation for Slater Determinants and Configuration Interaction States

Pietro Silvi, Davide Rossini, Giuseppe E. Santoro +1 more

Slater determinants are product states of filled quantum fermionic orbitals. When they are expressed in a configuration space basis chosen a priori, their entanglement is bound and controlled. This suggests that an exact representation of Slater determinants as finitely-correlated states is possible. In this paper we analyze this issue and provide an exact Matrix Product representation for Slater determinant states. We also argue possible meaningful extensions that embed more complex configuration interaction states into the description.

Jan 1, 20147 citations

Model Wavefunction Forms to Describe Strong Correlation in Quantum Chemistry

Paul A. Johnson

Standard approaches to approximating solutions to the Schroedinger equation rely upon the assumption that there is a single dominant electronic configuration with a small number of low-lying excited states. As such, they are referred to as single reference approaches and are accurate only when the band gap is small in comparison with the inter-electronic repulsion of the valence electrons. Systems for which this is the case are weakly-correlated. If the inter-electronic repulsion larger than the band gap, then a single electronic reference is not an appropriate ansatz for the physical wavefunction. As the strength of the interaction is increased many Slater determinants become near-degenerate, each providing non-negligible contributions. Such systems are strongly-correlated. The goal of the thesis was to employ wavefunctions which are: non-trivial, tractable, size-consistent and extensive, and not given in an active space. We have outlined approaches based on wavefunction forms more suited to describe strong correlation. Built on the Lie algebra su(2), our first trial wavefunction was the Bethe ansatz solution to the reduced BCS Hamiltonian derived by Richardson. This wavefunction is expressed as a product of quasi-particles each parametrized in terms of quasi-momenta. Consistency among the quasi-momenta is ensured by a set of non-linear equations to be solved numerically. This approach is feasible in both variational and projected flavours. We next considered a wavefunction of the same form without requiring it to be an eigenvector of the reduced BCS Hamiltonian. A variational approach in this case is only feasible provided there are two equivalent descriptions in terms of particles and holes. In any case, a bivariational principle and a projected approach are both feasible. Also based on su(2), there is AP1roG (coupled cluster pair doubles). Open-shell singlets and triplets were constructed with the Lie algebras sp(N) and so(4N). The resulting formulae appear to be infeasible unless it is possible to efficiently compute generalizations of determinants to three-dimensional objects. Treating each spatial orbital individually, all four possible occupations close the Lie superalgebra gl(2|2).

Jan 6, 20153 citations