Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

Physical Review

Variational Approach to the Many-Body Problem

R. U. Ayres

Published Sep 15, 1958• 39 citations
Publisher

Abstract

Certain problems connected with the choice of trial variational functions for the quantum many-body problem are discussed. It is convenient to make use of the reduced-density-matrix formulation developed by L\"owdin and Mayer, since a system whose Hamiltonian contains only 1- and 2-body interactions is characterized by the two simplest reduced density matrices ${\ensuremath{\Gamma}}^{(1)}$ and ${\ensuremath{\Gamma}}^{(2)}$. It is known that symmetry, Hermiticity, and normalizability conditions are not sufficient to ensure a physically realizable choice of the functions ${\ensuremath{\Gamma}}^{(1)}$ and ${\ensuremath{\Gamma}}^{(2)}$. However an additional nontrivial restriction is imposed by the fact that such functions must be possible variational extrema. The implications of this condition are investigated in some detail and it turns out that all the nonphysical choices of ${\ensuremath{\Gamma}}^{(1)}$ and ${\ensuremath{\Gamma}}^{(2)}$ which have been exhibited by various authors are eliminated. The condition is essentially nonquantum inasmuch as an analogous argument should be possible for a classical system. Moreover one of the consequences of this condition is that the system (classical or quantum) may exhibit macroscopic behavior, i.e., order-disorder transitions. The analysis leads to apparently reasonable choices of trial forms for ${\ensuremath{\Gamma}}^{(1)}$ and ${\ensuremath{\Gamma}}^{(2)}$, although it is still not known whether any further restriction must be imposed to complete the sufficiency argument.

DOI
10.1103/physrev.111.1453

Related papers

Reduction of the N -Particle Variational Problem

Claude Garrod, Jerome K. Percus

A variational method is presented which is applicable to N-particle boson or fermion systems with two-body interactions. For these systems the energy may be expressed in terms of the two-particle density matrix: Γ(1, 2 | 1′, 2′)=(Ψ |a2+a1+a1′a2′| Ψ). In order to have the variational equation: δE/δΓ = 0 yield the correct ground-state density matrix one must restrict Γ to the set of density matrices which are actually derivable from N-particle boson (or fermion) systems. Subsidiary conditions are presented which are necessary and sufficient to insure that Γ is so derivable. These conditions are of a form which render them unsuited for practical application. However the following necessary (but not sufficient) conditions are shown by some applications to yield good results: It is proven that if Γ(1, 2 | 1′, 2′) and γ(1 | 1′) are the two-particle and one-particle density matrices of an N-particle system [normalized by tr Γ = N(N − 1) and trγ = N] then the associated operator: G(1,2 | 1′,2′)=δ(1−1′)γ(2 | 2′)+σ Γ(1′,2 | 1,2′)−γ(2 | 1)γ(1′ | 2′) is a nonnegative operator. [Here σ is + 1 or − 1 for bosons or fermions respectively.]

Journal of Mathematical Physics•Dec 1, 1964•477 citations

Variational principles for the determination of single-particle density matrices

Viraht Sahni, Joseph B. Krieger

A variational formalism for the determination of the spinless single-particle density matrix $\ensuremath{\gamma}(\stackrel{\ensuremath{\rightarrow}}{\mathrm{r}}, {\stackrel{\ensuremath{\rightarrow}}{\mathrm{r}}}^{\ensuremath{'}})$ for an interacting many-particle system, correct to second order, is developed. The method is based on variational principles for determining single-particle expectation values, correct to second order, without requiring highly accurate wave functions in order to represent the interacting system. The expectation value of the operator $X=Y+iZ$, where $Y=(\frac{1}{2}){\ensuremath{\Sigma}}_{j}[W({\stackrel{\ensuremath{\rightarrow}}{\mathrm{r}}}_{i}){T}_{j}(\stackrel{\ensuremath{\rightarrow}}{\mathrm{a}})+W({\stackrel{\ensuremath{\rightarrow}}{\mathrm{r}}}_{j}^{\ensuremath{'}}){T}_{j}(\ensuremath{-}\stackrel{\ensuremath{\rightarrow}}{\mathrm{a}})]$, $Z=\ensuremath{-}(\frac{i}{2}){\ensuremath{\Sigma}}_{j}[W({\stackrel{\ensuremath{\rightarrow}}{\mathrm{r}}}_{j}){T}_{j}(\stackrel{\ensuremath{\rightarrow}}{\mathrm{a}})\ensuremath{-}W({\stackrel{\ensuremath{\rightarrow}}{\mathrm{r}}}_{j}^{\ensuremath{'}}){T}_{j}(\ensuremath{-}\stackrel{\ensuremath{\rightarrow}}{\mathrm{a}})]$, $W({\stackrel{\ensuremath{\rightarrow}}{\mathrm{r}}}_{j})=\ensuremath{\delta}({\stackrel{\ensuremath{\rightarrow}}{\mathrm{r}}}_{j}\ensuremath{-}\stackrel{\ensuremath{\rightarrow}}{\mathrm{r}})$, ${T}_{j}(\stackrel{\ensuremath{\rightarrow}}{\mathrm{a}})$ is a translation operator, and $\stackrel{\ensuremath{\rightarrow}}{\mathrm{a}}={\stackrel{\ensuremath{\rightarrow}}{\mathrm{r}}}^{\ensuremath{'}}\ensuremath{-}\stackrel{\ensuremath{\rightarrow}}{\mathrm{r}}$, is the reduced single-particle density matrix. For a given choice of a Slater-determinant-type trial wave function it is shown that the density matrix may be obtained from the charge density of the system where the latter too is determined via the same formalism to second order as the expectation value of the operator $W$ employing the same system trial wave function. The applicability of the technique is then demonstrated by showing that for the model hydrogen-atom problem it leads to highly accurate results for the momentum density and the Compton profile in the impulse approximation, and that the virial theorem is closely satisfied, even when the calculations of the various properties employing an approximate trial wave function alone are substantially in error.

Physical Review A•Feb 1, 1975•19 citations

Classical simulation of quantum many-body systems

Yichen Huang

Classical simulation of quantum many-body systems is in general a challenging problem for the simple reason that the dimension of the Hilbert space grows exponentially with the system size. In particular, merely encoding a generic quantum many-body state requires an exponential number of bits. However, condensed matter physicists are mostly interested in local Hamiltonians and especially their ground states, which are highly non-generic. Thus, we might hope that at least some physical systems allow efficient classical simulation.Starting with one-dimensional (1D) quantum systems (i.e., the simplest nontrivial case), the first basic question is: Which classes of states have efficient classical representations? It turns out that this question is quantitatively related to the amount of entanglement in the state, for states with ``little entanglement'' are well approximated by matrix product states (a data structure that can be manipulated efficiently on a classical computer). At a technical level, the mathematical notion for ``little entanglement'' is area law, which has been proved for unique ground states in 1D gapped systems. We establish an area law for constant-fold degenerate ground states in 1D gapped systems and thus explain the effectiveness of matrix-product-state methods in (e.g.) symmetry breaking phases. This result might not be intuitively trivial as degenerate ground states in gapped systems can be long-range correlated.Suppose an efficient classical representation exists. How can one find it efficiently? The density matrix renormalization group is the leading numerical method for computing ground states in 1D quantum systems. However, it is a heuristic algorithm and the possibility that it may fail in some cases cannot be completely ruled out. Recently, a provably efficient variant of the density matrix renormalization group has been developed for frustration-free 1D gapped systems. We generalize this algorithm to all (i.e., possibly frustrated) 1D gapped systems. Note that the ground-state energy of 1D gapless Hamiltonians is computationally intractable even in the presence of translational invariance.It is tempting to extend methods and tools in 1D to two and higher dimensions (2+D), e.g., matrix product states are generalized to tensor network states. Since an area law for entanglement (if formulated properly) implies efficient matrix product state representations in 1D, an interesting question is whether a similar implication holds in 2+D. Roughly speaking, we show that an area law for entanglement (in any reasonable formulation) does not always imply efficient tensor network representations of the ground states of 2+D local Hamiltonians even in the presence of translational invariance. It should be emphasized that this result does not contradict with the common sense that in practice quantum states with more entanglement usually require more space to be stored classically; rather, it demonstrates that the relationship between entanglement and efficient classical representations is still far from being well understood.Excited eigenstates participate in the dynamics of quantum systems and are particularly relevant to the phenomenon of many-body localization (absence of transport at finite temperature in strongly correlated systems). We study the entanglement of excited eigenstates in random spin chains and expect that its singularities coincide with dynamical quantum phase transitions. This expectation is confirmed in the disordered quantum Ising chain using both analytical and numerical methods.Finally, we study the problem of generating ground states (possibly with topological order) in 1D gapped systems using quantum circuits. This is an interesting problem both in theory and in practice. It not only characterizes the essential difference between the entanglement patterns that give rise to trivial and nontrivial topological order, but also quantifies the difficulty of preparing quantum states with a quantum computer (in experiments).

•Jan 1, 2015•19 citations

Variational determination of the two-particle density matrix as a quantum many-body technique

Brecht Verstichel

In this thesis the variational optimisation of the density matrix is discussed as a method in many-body quantum mechanics. This is a relatively unknown technique in which one tries to obtain the two-particle reduced density matrix directly in a variational approach. In the first Chapter we introduce the subject in the broader context of many-body quantum mechanics, and briefly sketch the history of the field. The second Chapter tries to summarise what is known about N-representability of reduced density matrices, and derives these results in a unified framework. The optimisation problem can be formulated as a semidefinite program. Three different algorithms are discussed in Chapter three, and their performance is compared. In Chapter four we show how symmetry can be exploited to reduce the computational cost considerably. Several applications of the method are discussed in Chapter five. We show that the standard two-index constraints fail in the strong-correlation limit of the one-dimensional Hubbard model. In Chapter six we identify the spin-adapted lifting conditions as the most compact constraints that can correctly describe this limit.

arXiv:1203.5659•Mar 26, 2012•4 citations

Remarks on Mayer's Reduced Density Matrix Method

Yukio Mizuno, Takeo Izuyama

Mayer has proposed a variation method of finding the two-body reduced density matrix ρ2 for the ground state of a many fermion system. In this paper it is shown explicitly that his variation method is imperfect and sometimes leads to quite unsatisfactory results. Further, his trial form of ρ2 for electron gas is translated into that for the π electron system of benzene molecule and it is shown that such a trial form is unsatisfactory in several respects.

Progress of Theoretical Physics•Jul 1, 1957•24 citations

Reexamining the variational two-particle reduced density matrix for nuclear systems

J. G. Li, N. Michel, Zuo Wei +1 more

In most nuclear many-body methods, observables are calculated using many-body wave functions explicitly. The variational two-particle reduced density matrix method is one of the few exceptions to the rule. Ground-state energies of both closed-shell and open-shell nuclear systems can indeed be evaluated by minimizing a constrained linear functional of the two-particle reduced density matrix. However, it has virtually never been used in nuclear theory, because nuclear ground states were found to be well overbound, contrary to those of atoms and molecules. Consequently, we introduced new constraints in the nuclear variational two-particle reduced density matrix method, developed recently for atomic and molecular systems. Our calculations then show that this approach can provide a proper description of nuclear systems where only valence neutrons are included. For the nuclear systems where both neutrons and protons are active, however, the energies obtained with the variational two-particle reduced density matrix method are still overbound. The possible reasons for the noticed discrepancies and solutions to this problem will be discussed.

arXiv:2106.04103•Physical Review C•Jun 23, 2021•2 citations
Certain problems connected with the choice of trial variational functions for the quantum many-body problem are discussed. It is convenient to make use of the reduced-density-matrix formulation developed by L"owdin and Mayer, since a system whose Hamiltonian contains only 1- and 2-body interactions is characterized by the two simplest reduced density matrices \ensuremathΓ(1){\ensuremath{\Gamma}}^{(1)}\ensuremathΓ(1) and \ensuremathΓ(2){\ensuremath{\Gamma}}^{(2)}\ensuremathΓ(2). It is known that symmetry, Hermiticity, and normalizability conditions are not sufficient to ensure a physically realizable choice of the functions \ensuremathΓ(1){\ensuremath{\Gamma}}^{(1)}\ensuremathΓ(1) and \ensuremathΓ(2){\ensuremath{\Gamma}}^{(2)}\ensuremathΓ(2). However an additional nontrivial restriction is imposed by the fact that such functions must be possible variational extrema. The implications of this condition are investigated in some detail and it turns out that all the nonphysical choices of \ensuremathΓ(1){\ensuremath{\Gamma}}^{(1)}\ensuremathΓ(1) and \ensuremathΓ(2){\ensuremath{\Gamma}}^{(2)}\ensuremathΓ(2) which have been exhibited by various authors are eliminated. The condition is essentially nonquantum inasmuch as an analogous argument should be possible for a classical system. Moreover one of the consequences of this condition is that the system (classical or quantum) may exhibit macroscopic behavior, i.e., order-disorder transitions. The analysis leads to apparently reasonable choices of trial forms for \ensuremathΓ(1){\ensuremath{\Gamma}}^{(1)}\ensuremathΓ(1) and \ensuremathΓ(2){\ensuremath{\Gamma}}^{(2)}\ensuremathΓ(2), although it is still not known whether any further restriction must be imposed to complete the sufficiency argument.
A variational method is presented which is applicable to N-particle boson or fermion systems with two-body interactions. For these systems the energy may be expressed in terms of the two-particle density matrix: Γ(1, 2 | 1′, 2′)=(Ψ |a2+a1+a1′a2′| Ψ). In order to have the variational equation: δE/δΓ = 0 yield the correct ground-state density matrix one must restrict Γ to the set of density matrices which are actually derivable from N-particle boson (or fermion) systems. Subsidiary conditions are presented which are necessary and sufficient to insure that Γ is so derivable. These conditions are of a form which render them unsuited for practical application. However the following necessary (but not sufficient) conditions are shown by some applications to yield good results: It is proven that if Γ(1, 2 | 1′, 2′) and γ(1 | 1′) are the two-particle and one-particle density matrices of an N-particle system [normalized by tr Γ = N(N − 1) and trγ = N] then the associated operator: G(1,2 | 1′,2′)=δ(1−1′)γ(2 | 2′)+σ Γ(1′,2 | 1,2′)−γ(2 | 1)γ(1′ | 2′) is a nonnegative operator. [Here σ is + 1 or − 1 for bosons or fermions respectively.]
A variational formalism for the determination of the spinless single-particle density matrix $\ensuremath{\gamma}(\stackrel{\ensuremath{\rightarrow}}{\mathrm{r}}, {\stackrel{\ensuremath{\rightarrow}}{\mathrm{r}}}^{\ensuremath{'}})$ for an interacting many-particle system, correct to second order, is developed. The method is based on variational principles for determining single-particle expectation values, correct to second order, without requiring highly accurate wave functions in order to represent the interacting system. The expectation value of the operator X=Y+iZX=Y+iZX=Y+iZ, where $Y=(\frac{1}{2}){\ensuremath{\Sigma}}_{j}[W({\stackrel{\ensuremath{\rightarrow}}{\mathrm{r}}}_{i}){T}_{j}(\stackrel{\ensuremath{\rightarrow}}{\mathrm{a}})+W({\stackrel{\ensuremath{\rightarrow}}{\mathrm{r}}}_{j}^{\ensuremath{'}}){T}_{j}(\ensuremath{-}\stackrel{\ensuremath{\rightarrow}}{\mathrm{a}})]$, $Z=\ensuremath{-}(\frac{i}{2}){\ensuremath{\Sigma}}_{j}[W({\stackrel{\ensuremath{\rightarrow}}{\mathrm{r}}}_{j}){T}_{j}(\stackrel{\ensuremath{\rightarrow}}{\mathrm{a}})\ensuremath{-}W({\stackrel{\ensuremath{\rightarrow}}{\mathrm{r}}}_{j}^{\ensuremath{'}}){T}_{j}(\ensuremath{-}\stackrel{\ensuremath{\rightarrow}}{\mathrm{a}})]$, $W({\stackrel{\ensuremath{\rightarrow}}{\mathrm{r}}}_{j})=\ensuremath{\delta}({\stackrel{\ensuremath{\rightarrow}}{\mathrm{r}}}_{j}\ensuremath{-}\stackrel{\ensuremath{\rightarrow}}{\mathrm{r}})$, ${T}_{j}(\stackrel{\ensuremath{\rightarrow}}{\mathrm{a}})$ is a translation operator, and $\stackrel{\ensuremath{\rightarrow}}{\mathrm{a}}={\stackrel{\ensuremath{\rightarrow}}{\mathrm{r}}}^{\ensuremath{'}}\ensuremath{-}\stackrel{\ensuremath{\rightarrow}}{\mathrm{r}}$, is the reduced single-particle density matrix. For a given choice of a Slater-determinant-type trial wave function it is shown that the density matrix may be obtained from the charge density of the system where the latter too is determined via the same formalism to second order as the expectation value of the operator WWW employing the same system trial wave function. The applicability of the technique is then demonstrated by showing that for the model hydrogen-atom problem it leads to highly accurate results for the momentum density and the Compton profile in the impulse approximation, and that the virial theorem is closely satisfied, even when the calculations of the various properties employing an approximate trial wave function alone are substantially in error.
Classical simulation of quantum many-body systems is in general a challenging problem for the simple reason that the dimension of the Hilbert space grows exponentially with the system size. In particular, merely encoding a generic quantum many-body state requires an exponential number of bits. However, condensed matter physicists are mostly interested in local Hamiltonians and especially their ground states, which are highly non-generic. Thus, we might hope that at least some physical systems allow efficient classical simulation.Starting with one-dimensional (1D) quantum systems (i.e., the simplest nontrivial case), the first basic question is: Which classes of states have efficient classical representations? It turns out that this question is quantitatively related to the amount of entanglement in the state, for states with ``little entanglement'' are well approximated by matrix product states (a data structure that can be manipulated efficiently on a classical computer). At a technical level, the mathematical notion for ``little entanglement'' is area law, which has been proved for unique ground states in 1D gapped systems. We establish an area law for constant-fold degenerate ground states in 1D gapped systems and thus explain the effectiveness of matrix-product-state methods in (e.g.) symmetry breaking phases. This result might not be intuitively trivial as degenerate ground states in gapped systems can be long-range correlated.Suppose an efficient classical representation exists. How can one find it efficiently? The density matrix renormalization group is the leading numerical method for computing ground states in 1D quantum systems. However, it is a heuristic algorithm and the possibility that it may fail in some cases cannot be completely ruled out. Recently, a provably efficient variant of the density matrix renormalization group has been developed for frustration-free 1D gapped systems. We generalize this algorithm to all (i.e., possibly frustrated) 1D gapped systems. Note that the ground-state energy of 1D gapless Hamiltonians is computationally intractable even in the presence of translational invariance.It is tempting to extend methods and tools in 1D to two and higher dimensions (2+D), e.g., matrix product states are generalized to tensor network states. Since an area law for entanglement (if formulated properly) implies efficient matrix product state representations in 1D, an interesting question is whether a similar implication holds in 2+D. Roughly speaking, we show that an area law for entanglement (in any reasonable formulation) does not always imply efficient tensor network representations of the ground states of 2+D local Hamiltonians even in the presence of translational invariance. It should be emphasized that this result does not contradict with the common sense that in practice quantum states with more entanglement usually require more space to be stored classically; rather, it demonstrates that the relationship between entanglement and efficient classical representations is still far from being well understood.Excited eigenstates participate in the dynamics of quantum systems and are particularly relevant to the phenomenon of many-body localization (absence of transport at finite temperature in strongly correlated systems). We study the entanglement of excited eigenstates in random spin chains and expect that its singularities coincide with dynamical quantum phase transitions. This expectation is confirmed in the disordered quantum Ising chain using both analytical and numerical methods.Finally, we study the problem of generating ground states (possibly with topological order) in 1D gapped systems using quantum circuits. This is an interesting problem both in theory and in practice. It not only characterizes the essential difference between the entanglement patterns that give rise to trivial and nontrivial topological order, but also quantifies the difficulty of preparing quantum states with a quantum computer (in experiments).
In this thesis the variational optimisation of the density matrix is discussed as a method in many-body quantum mechanics. This is a relatively unknown technique in which one tries to obtain the two-particle reduced density matrix directly in a variational approach. In the first Chapter we introduce the subject in the broader context of many-body quantum mechanics, and briefly sketch the history of the field. The second Chapter tries to summarise what is known about N-representability of reduced density matrices, and derives these results in a unified framework. The optimisation problem can be formulated as a semidefinite program. Three different algorithms are discussed in Chapter three, and their performance is compared. In Chapter four we show how symmetry can be exploited to reduce the computational cost considerably. Several applications of the method are discussed in Chapter five. We show that the standard two-index constraints fail in the strong-correlation limit of the one-dimensional Hubbard model. In Chapter six we identify the spin-adapted lifting conditions as the most compact constraints that can correctly describe this limit.
In most nuclear many-body methods, observables are calculated using many-body wave functions explicitly. The variational two-particle reduced density matrix method is one of the few exceptions to the rule. Ground-state energies of both closed-shell and open-shell nuclear systems can indeed be evaluated by minimizing a constrained linear functional of the two-particle reduced density matrix. However, it has virtually never been used in nuclear theory, because nuclear ground states were found to be well overbound, contrary to those of atoms and molecules. Consequently, we introduced new constraints in the nuclear variational two-particle reduced density matrix method, developed recently for atomic and molecular systems. Our calculations then show that this approach can provide a proper description of nuclear systems where only valence neutrons are included. For the nuclear systems where both neutrons and protons are active, however, the energies obtained with the variational two-particle reduced density matrix method are still overbound. The possible reasons for the noticed discrepancies and solutions to this problem will be discussed.