© 2026 Aqora Quantum S.A.S.

Terms Privacy Legal Notice

Physical Review B

Implications of the two nodal domains conjecture for ground state fermionic wave functions

Dario Bressanini

Published Sep 13, 2012• 27 citations

Abstract

The nodes of many-body wave functions are mathematical objects important in many different fields of physics. They are at the heart of the quantum Monte Carlo methods but outside this field their properties are neither widely known nor studied. In recent years a conjecture, already proven to be true in several important cases, has been put forward related to the nodes of the fermionic ground state of a many-body system, namely that there is a single nodal hypersurface that divides configuration space into only two connected domains. While this is obviously relevant to the fixed node diffusion Monte Carlo method, its repercussions have ramifications in various fields of physics as diverse as density functional theory or Feynman and Cohen's backflow wave function formulation. To illustrate this we explicitly show that, even if we knew the exact Kohn-Sham exchange correlation functional, there are systems for which we would obtain the exact ground state energy and density but a wave function quite different from the exact one. This paradox is only apparent since the Hohenberg-Kohn theorem relates the energy directly to the density and the wave function is not guaranteed to be close to the exact one. The aim of this paper is to stimulate the investigation of the properties of the nodes of many-body wave functions in different fields of physics. Furthermore, we explicitly show that this conjecture is related to the phenomenon of avoided nodal crossing but it is not necessarily caused by electron correlation, as sometimes has been suggested in the literature. We explicitly build a many-body uncorrelated example whose nodal structure shows the same phenomenon.

DOI: 10.1103/physrevb.86.115120

Related papers

Fermion nodes and nodal cells of noninteracting and interacting fermions

Luboš Mitáš

A fermion node is subset of fermionic configurations for which a real wave function vanishes due to the antisymmetry and the node divides the configurations space into compact nodal cells (domains). We analyze the properties of fermion nodes of fermionic ground state wave functions for a number of systems. For several models we demonstrate that noninteracting spin-polarized fermions in dimension two and higher have closed-shell ground state wave functions with the minimal two nodal cells for any system size and we formulate a theorem which sumarizes this result. The models include periodic fermion gas, fermions on the surface of a sphere, fermions in a box. We prove the same property for atomic states with up to 3d half-filled shells. Under rather general assumptions we then derive that the same is true for unpolarized systems with arbitrarily weak interactions using Bardeen-Cooper-Schrieffer (BCS) variational wave function. We further show that pair correlations included in the BCS wave function enable singlet pairs of particles to wind around the periodic box without crossing the node pointing towards the relationship of nodes to transport and many-body phases such as superconductivity. Finally, we point out that the arguments extend also to fermionic temperature dependent/imaginary-time density matrices. The results reveal fundamental properties of fermion nodal structures and provide new insights for accurate constructions of wave functions and density matrices in quantum and path integral Monte Carlo methods.

arXiv:cond-mat/0605550•May 23, 2006•5 citations

Fixed-node Monte Carlo study of the two-dimensional Hubbard model

Guozhong An, J. M. J. van Leeuwen

The fixed-node Monte Carlo method is extended to lattice fermion models. This method replaces the problem of finding the ground state of a many-fermion system by an effective eigenvalue problem of finding the lowest-energy wave function in a given region of the configuration space. It has previously only been applied to fermions moving in continuous space. The discreteness of the configuration space causes the algorithm to differ from that of the continuum. The method is tested against a known limiting case where exact results are available for comparison. Good agreement is found. Using the fixed-node Monte Carlo method, we study the domain-wall phase in the ground state of the two-dimensional Hubbard model. The existence of a domain-wall phase which has domains of antiferromagnetic phases separated by walls of holes is recently suggested by inhomogeneous Hartree-Fock (HF) and variational Monte Carlo (VMC) calculations. Large improvements of the energies are found. The domain walls are broader than those obtained by the HF and VMC calculations.

Physical review. B, Condensed matter•Nov 1, 1991•29 citations

Model Hamiltonians in density functional theory

Paola Gori‐Giorgi, Julien Toulouse, Andreas Savin B

The formalism of Kohn and Sham uses a specific (model) hamiltonian which highly simplifies the many-electron problem to that of noninteracting fermions. The theorem of Hohenberg and Kohn tells us that, for a given ground state density, this hamiltonian is unique. In principle, this density can be chosen as that of the real, interacting system. To obtain the energy, or other properties of the real system, approximations are needed. Working with non interacting fermions is an important simplification, but it may be easier to produce approximations with different choices of the model hamiltonian. The feature that the exact density is (ideally) reproduced can be kept in the newly defined fictitious systems. Using model hamiltonians having the same form as the physical one, that is, being built of one- and two-body operators, allows to approach the physical hamiltonian arbitrarily close, and thus a systematic reduction of the approximations.

arXiv:cond-mat/0605174•CRM proceedings & lecture notes•Jun 20, 2007•0 citations

Exact exchange-correlation potential of effectively interacting Kohn-Sham systems

Shunsuke A. Sato, Ángel Rubio

Aiming to combine density functional theory (DFT) and wave-function theory, we study a mapping from the many-body interacting system to an effectively interacting Kohn-Sham system instead of a noninteracting Kohn-Sham system. Because a ground state of effectively interacting systems requires having a solution for the correlated many-body wave functions, this provides a natural framework to many-body wave-function theories such as the configuration interaction and the coupled-cluster method in the formal theoretical framework of DFT. Employing simple one-dimensional two-electron systems---namely, the one-dimensional helium atom, the hydrogen molecule, and the heteronuclear diatomic molecule---we investigate properties of many-body wave functions and exact exchange-correlation potentials of effectively interacting Kohn-Sham systems. As a result, we find that the asymptotic behavior of the exact exchange-correlation potential can be controlled by optimizing that of the effective interaction. Furthermore, the typical features of the exact noninteracting Kohn-Sham system, namely, a spiky feature and a step feature in the exchange-correlation potential for the molecular dissociation limit, can be suppressed by a proper choice of the effective interaction. These findings open a possibility to construct numerically robust and efficient exchange-correlation potentials and functionals based on the effectively interacting Kohn-Sham scheme.

arXiv:1907.00183•Physical Review A•Jan 15, 2020•1 citation

Many-body energy density functional

Alejandro Kievsky, Giuseppina Orlandini, M. Gattobigio

The Hohenberg-Kohn theorem and the Kohn-Sham equations, which are at the basis of the density-functional theory, are reformulated in terms of a particular many-body density, which is translational and Galilean invariant and therefore is relevant for self-bound systems. In a similar way that there is a unique relation between the one-body density and the external potential that gives rise to it, we demonstrate that there is a unique relation between that particular many-body density and a definite many-body potential. The energy is then a functional of this density, and its minimization leads to the ground-state energy of the system. As a proof of principle, the analogous of the Kohn-Sham equation is solved in the specific case of $^{4}\mathrm{He}$ atomic clusters, to put in evidence the advantages of this formulation in terms of physical insights.

arXiv:2106.13582•Physical Review A•Sep 20, 2021•2 citations

Generalized Pairing Wave Functions and Nodal Properties for Electronic Structure Quantum Monte Carlo

Michal Bajdich

The quantum Monte Carlo (QMC) is one of the most promising many-body electronic structure approaches. It employs stochastic techniques for solving the stationary Schr\" odinger equation and for evaluation of expectation values. The key advantage of QMC is its capability to use the explicitly correlated wave functions, which allow the study of many-body effects beyond the reach of mean-field methods. The most important limit on QMC accuracy is the fixed-node approximation, which comes from necessity to circumvent the fermion sign problem. The size of resulting fixed-node errors depends on the quality of the nodes (the subset of position space where the wave function vanishes) of a used wave function. In this dissertation, we analyze the nodal properties of the existing fermionic wave functions and offer new types of variational wave functions with improved nodal structure.

arXiv:0712.3066•Dec 18, 2007•0 citations

The nodes of many-body wave functions are mathematical objects important in many different fields of physics. They are at the heart of the quantum Monte Carlo methods but outside this field their properties are neither widely known nor studied. In recent years a conjecture, already proven to be true in several important cases, has been put forward related to the nodes of the fermionic ground state of a many-body system, namely that there is a single nodal hypersurface that divides configuration space into only two connected domains. While this is obviously relevant to the fixed node diffusion Monte Carlo method, its repercussions have ramifications in various fields of physics as diverse as density functional theory or Feynman and Cohen's backflow wave function formulation. To illustrate this we explicitly show that, even if we knew the exact Kohn-Sham exchange correlation functional, there are systems for which we would obtain the exact ground state energy and density but a wave function quite different from the exact one. This paradox is only apparent since the Hohenberg-Kohn theorem relates the energy directly to the density and the wave function is not guaranteed to be close to the exact one. The aim of this paper is to stimulate the investigation of the properties of the nodes of many-body wave functions in different fields of physics. Furthermore, we explicitly show that this conjecture is related to the phenomenon of avoided nodal crossing but it is not necessarily caused by electron correlation, as sometimes has been suggested in the literature. We explicitly build a many-body uncorrelated example whose nodal structure shows the same phenomenon.

A fermion node is subset of fermionic configurations for which a real wave function vanishes due to the antisymmetry and the node divides the configurations space into compact nodal cells (domains). We analyze the properties of fermion nodes of fermionic ground state wave functions for a number of systems. For several models we demonstrate that noninteracting spin-polarized fermions in dimension two and higher have closed-shell ground state wave functions with the minimal two nodal cells for any system size and we formulate a theorem which sumarizes this result. The models include periodic fermion gas, fermions on the surface of a sphere, fermions in a box. We prove the same property for atomic states with up to 3d half-filled shells. Under rather general assumptions we then derive that the same is true for unpolarized systems with arbitrarily weak interactions using Bardeen-Cooper-Schrieffer (BCS) variational wave function. We further show that pair correlations included in the BCS wave function enable singlet pairs of particles to wind around the periodic box without crossing the node pointing towards the relationship of nodes to transport and many-body phases such as superconductivity. Finally, we point out that the arguments extend also to fermionic temperature dependent/imaginary-time density matrices. The results reveal fundamental properties of fermion nodal structures and provide new insights for accurate constructions of wave functions and density matrices in quantum and path integral Monte Carlo methods.

The fixed-node Monte Carlo method is extended to lattice fermion models. This method replaces the problem of finding the ground state of a many-fermion system by an effective eigenvalue problem of finding the lowest-energy wave function in a given region of the configuration space. It has previously only been applied to fermions moving in continuous space. The discreteness of the configuration space causes the algorithm to differ from that of the continuum. The method is tested against a known limiting case where exact results are available for comparison. Good agreement is found. Using the fixed-node Monte Carlo method, we study the domain-wall phase in the ground state of the two-dimensional Hubbard model. The existence of a domain-wall phase which has domains of antiferromagnetic phases separated by walls of holes is recently suggested by inhomogeneous Hartree-Fock (HF) and variational Monte Carlo (VMC) calculations. Large improvements of the energies are found. The domain walls are broader than those obtained by the HF and VMC calculations.

The formalism of Kohn and Sham uses a specific (model) hamiltonian which highly simplifies the many-electron problem to that of noninteracting fermions. The theorem of Hohenberg and Kohn tells us that, for a given ground state density, this hamiltonian is unique. In principle, this density can be chosen as that of the real, interacting system. To obtain the energy, or other properties of the real system, approximations are needed. Working with non interacting fermions is an important simplification, but it may be easier to produce approximations with different choices of the model hamiltonian. The feature that the exact density is (ideally) reproduced can be kept in the newly defined fictitious systems. Using model hamiltonians having the same form as the physical one, that is, being built of one- and two-body operators, allows to approach the physical hamiltonian arbitrarily close, and thus a systematic reduction of the approximations.

Aiming to combine density functional theory (DFT) and wave-function theory, we study a mapping from the many-body interacting system to an effectively interacting Kohn-Sham system instead of a noninteracting Kohn-Sham system. Because a ground state of effectively interacting systems requires having a solution for the correlated many-body wave functions, this provides a natural framework to many-body wave-function theories such as the configuration interaction and the coupled-cluster method in the formal theoretical framework of DFT. Employing simple one-dimensional two-electron systems---namely, the one-dimensional helium atom, the hydrogen molecule, and the heteronuclear diatomic molecule---we investigate properties of many-body wave functions and exact exchange-correlation potentials of effectively interacting Kohn-Sham systems. As a result, we find that the asymptotic behavior of the exact exchange-correlation potential can be controlled by optimizing that of the effective interaction. Furthermore, the typical features of the exact noninteracting Kohn-Sham system, namely, a spiky feature and a step feature in the exchange-correlation potential for the molecular dissociation limit, can be suppressed by a proper choice of the effective interaction. These findings open a possibility to construct numerically robust and efficient exchange-correlation potentials and functionals based on the effectively interacting Kohn-Sham scheme.

The Hohenberg-Kohn theorem and the Kohn-Sham equations, which are at the basis of the density-functional theory, are reformulated in terms of a particular many-body density, which is translational and Galilean invariant and therefore is relevant for self-bound systems. In a similar way that there is a unique relation between the one-body density and the external potential that gives rise to it, we demonstrate that there is a unique relation between that particular many-body density and a definite many-body potential. The energy is then a functional of this density, and its minimization leads to the ground-state energy of the system. As a proof of principle, the analogous of the Kohn-Sham equation is solved in the specific case of

^{4}\mathrm{He}

atomic clusters, to put in evidence the advantages of this formulation in terms of physical insights.

The quantum Monte Carlo (QMC) is one of the most promising many-body electronic structure approaches. It employs stochastic techniques for solving the stationary Schr\" odinger equation and for evaluation of expectation values. The key advantage of QMC is its capability to use the explicitly correlated wave functions, which allow the study of many-body effects beyond the reach of mean-field methods. The most important limit on QMC accuracy is the fixed-node approximation, which comes from necessity to circumvent the fermion sign problem. The size of resulting fixed-node errors depends on the quality of the nodes (the subset of position space where the wave function vanishes) of a used wave function. In this dissertation, we analyze the nodal properties of the existing fermionic wave functions and offer new types of variational wave functions with improved nodal structure.