We demonstrated in the previous chapter (Chapter 1) that the evaluation of matrix elements between states represented by the Slater determinant can be evaluated unambiguously using Slater’s rule [1, 28, 40]. While the determinantal representation and the evaluation matrix elements between the states are quite straightforward, this procedure can eventually be ineffective and cumbersome with increasing size of the system. The occupation number formalism or second quantization technique introduced by Dirac in quantum field theory, in this context has been found to be quite formidable and elegant. This formalism not only allows us to represent the states and the operators in a compact and convenient form but also enables us to deal with systems containing variable numbers of particles without explicitly specifying the number of particles. It implicitly assumes the existence of an unspecified number of functions in a single-particle basis. These one-particle functions, which constitute the Hilbert space for N particles, are generally chosen to be spinorbital because of their orthogonality. In the case of fermions, the many-particle basis is the anti-symmetric product of N one-particle basis functions (Slater determinant). The vector space consisting of all possible anti-symmetric Hilbert spaces is referred to as the Fock space, in which the operators and functions are represented without specifying the number of particles.
University of Colorado APPM preprint number 554. The wavefunction for the multiparticle Schrödinger equation is a function of many variables and satisfies an antisymmetry condition, so it is natural to approximate it as a sum of Slater determinants. Many current methods do so, but they impose additional structural constraints on the determinants, such as orthogonality between orbitals or an excitation pattern. We present a method without any such constraints, by which we hope to obtain much more efficient expansions, and insight into the inherent structure of the wavefunction. We use an integral formulation of the problem, a Green’s function iteration, and a fitting procedure based on the computational paradigm of separated representations. The core procedure is the construction and solution of a matrix-integral system derived from antisymmetric inner products involving the potential operators. We show how to construct and solve this system with computational complexity competitive with current methods.
The basic question is addressed of how the space dimension
d is encoded in the Hilbert space of
N identical fermions. There appears a finite number $N{!}^{d\ensuremath{-}1}$ of many-body wave functions, called shapes, which cannot be generated by trivial combinatorial extension of the one-dimensional ones. A general algorithm is given to list them all in terms of standard Slater determinants. Conversely, excitations which can be induced from the one-dimensional case are bosonized into a system of distinguishable bosons, called Euler bosons, much like the electromagnetic field is quantized in terms of photons distinguishable by their wave numbers. Their wave functions are given explicitly in terms of elementary symmetric functions, reflecting the fact that the fermion sign problem is trivial in one dimension. The shapes act as vacua for the Euler bosons. They are the natural generalization of the single-Slater-determinant form for the ground state to more than one dimension. In terms of algebraic invariant theory, the shapes are antisymmetric invariants which finitely generate the
N-fermion Hilbert space as a graded algebra over the ring of symmetric polynomials. Analogous results hold for identical bosons.
The wavefunction for the multiparticle Schrödinger equation is a function of many variables and satisfies an antisymmetry condition, so it is natural to approximate it as a sum of Slater determinants. Many current methods do so, but they impose additional structural constraints on the determinants, such as orthogonality between orbitals or an excitation pattern. We present a method without any such constraints, by which we hope to obtain much more efficient expansions and insight into the inherent structure of the wavefunction. We use an integral formulation of the problem, a Green’s function iteration, and a fitting procedure based on the computational paradigm of separated representations. The core procedure is the construction and solution of a matrix-integral system derived from antisymmetric inner products involving the potential operators. We show how to construct and solve this system with computational complexity competitive with current methods.
In order to calculate the average value of a physical quantity containing also many-particle interactions in a system of
N antisymmetric particles, a set of generalized density matrices are defined. In order to permit the investigation of the same physical situation in two complementary spaces, the Hermitean density matrix of order
k has two sets of indices of each
k variables, and it is further antisymmetric in each set of these indices.Every normalizable antisymmetric wave function may be expanded in a series of determinants of order
N over all ordered configurations formed from a basic complete set of one-particle functions ${\ensuremath{\psi}}_{k}$, which gives a representation of the wave function and its density matrices also in the discrete
k-space. The coefficients in an expansion of an eigenfunction to a particular operator may be determined by the variation principle, leading to the ordinary secular equation of the method of configurational interaction. It is shown that the first-order density matrix may be brought to diagonal form, which defines the "natural spin-orbitals" associated with the system. The situation is then partly characterized by the corresponding occupation numbers, which are shown to lie between 0 and 1 and to assume the value 1, only if the corresponding spin-orbital occurs in all configurations necessary for describing the situation. If the system has exactly
N spin-orbitals which are fully occupied, the total wave function may be reduced to a single Slater determinant. However, due to the mutual interaction between the particles, this limiting case is never physically realized, but the introduction of natural spin-orbitals leads then instead to a configurational expansion of most rapid convergence.In case the basic set is of finite order
M, the best choice of this set is determined by a form of extended Hartree-Fock equations. It is shown that, in this case, the natural spin-orbitals approximately fulfill some equations previously proposed by Slater.
We study separations between two fundamental models (or \emph{Ansätze}) of antisymmetric functions, that is, functions
f of the form
f(xσ(1),…,xσ(N))=sign(σ)f(x1,…,xN), where
σ is any permutation. These arise in the context of quantum chemistry, and are the basic modeling tool for wavefunctions of Fermionic systems. Specifically, we consider two popular antisymmetric Ansätze: the Slater representation, which leverages the alternating structure of determinants, and the Jastrow ansatz, which augments Slater determinants with a product by an arbitrary symmetric function. We construct an antisymmetric function in
N dimensions that can be efficiently expressed in Jastrow form, yet provably cannot be approximated by Slater determinants unless there are exponentially (in
N2) many terms. This represents the first explicit quantitative separation between these two Ansätze.