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

arXiv:1204.2006 · Physical Review B

Symmetry-projected variational approach for ground and excited states of the two-dimensional Hubbard model

R. R. Rodríguez-Guzmán, K.W. Schmid, Carlos A. Jiménez-Hoyos, Gustavo E. Scuseria

Published Jun 26, 2012• 48 citations
arXivHTMLPDF

Abstract

We present a symmetry-projected configuration mixing scheme to describe ground and excited states, with well defined quantum numbers, of the two-dimensional Hubbard model with nearest-neighbor hopping and periodic boundary conditions. Results for the half-filled $2\ifmmode\times\else\texttimes\fi{}4$, $4\ifmmode\times\else\texttimes\fi{}4$, and $6\ifmmode\times\else\texttimes\fi{}6$ lattices as well as doped $4\ifmmode\times\else\texttimes\fi{}4$ systems, compare well with available results, both exact and from other state-of-the-art approximations. We report spectral functions and density of states obtained from a well-controlled ansatz for the $({N}_{e}\ifmmode\pm\else\textpm\fi{}1)$-electron system. Symmetry projected methods have been widely used for the many-body nuclear physics problem but have received little attention in the solid state community. Given their relatively low (mean-field) computational cost and the high quality of results here reported, we believe that they deserve further scrutiny.

DOI
10.1103/physrevb.85.245130

Related papers

Multireference symmetry-projected variational approaches for ground and excited states of the one-dimensional Hubbard model

R. R. Rodríguez-Guzmán, Carlos A. Jiménez-Hoyos, Roman Schutski +1 more

We present a multireference configuration mixing scheme for describing ground and excited states, with well-defined spin and space-group symmetry quantum numbers, of the one-dimensional Hubbard model with nearest-neighbor hopping and periodic boundary conditions. Within this scheme, each state is expanded in terms of nonorthogonal and variationally determined symmetry-projected configurations. The results for lattices up to 30 and 50 sites compare well with the exact Lieb-Wu solutions as well as with results from other state-of-the-art approximations. In addition to spin-spin correlation functions in real space and magnetic structure factors, we present results for spectral functions and density of states computed with an Ansatz whose quality can be well controlled by the number of symmetry-projected configurations used to approximate the systems with ${N}_{e}$ and ${N}_{e}\ifmmode\pm\else\textpm\fi{}1$ electrons. The intrinsic symmetry-broken determinants resulting from the variational calculations have rich structures in terms of defects that can be regarded as basic units of quantum fluctuations. Given the quality of the results here reported, as well as the parallelization properties of the considered scheme, we believe that symmetry-projection techniques, which have found ample applications in nuclear structure physics, deserve further attention in the study of low-dimensional correlated many-electron systems.

arXiv:1304.4192•Physical Review B•Jun 24, 2013•52 citations

Symmetry-projected variational approach to the one-dimensional Hubbard model

K.W. Schmid, T. Dahm, Jérôme Margueron +1 more

We apply a variational method devised for the nuclear many-body problem to the one-dimensional Hubbard model with nearest neighbor hopping and periodic boundary conditions. The test wave function consist for each state out of a single Hartree--Fock determinant mixing all the sites (or momenta) as well as the spin projections of the electrons. Total spin and linear momentum are restored by projection methods before the variation. It is demonstrated that this approach reproduces the results of exact diagonalizations for half-filled $N=12$ and $N=14$ lattices not only for the energies and occupation numbers of the ground but also of the lowest excited states rather well. Furthermore, a system of ten electrons in an $N=12$ lattice is investigated and, finally, an $N=30$ lattice is studied. In addition to energies and occupation numbers we present the spectral functions computed with the help of the symmetry-projected wave functions as well.

arXiv:cond-mat/0409691•Physical Review B•Aug 12, 2005•19 citations

Multireference symmetry-projected variational approximation for the ground state of the doped one-dimensional Hubbard model

R. R. Rodríguez-Guzmán, Carlos A. Jiménez-Hoyos, Gustavo E. Scuseria

The few determinant (FED) approximation introduced in our previous work [Phys. Rev. B 87, 235129 (2013)] is used to describe the ground state, characterized by well-defined spin and space group symmetry quantum numbers as well as doping fractions ${N}_{e}/{N}_{\mathrm{sites}}$, of one-dimensional Hubbard lattices with nearest-neighbor hopping and periodic boundary conditions. Within this multireference scheme, each ground state is expanded in a given number of nonorthogonal and variationally determined symmetry-projected configurations. The results obtained for the ground-state and correlation energies of half-filled and doped lattices with 30, 34, and 50 sites compare well with the exact Lieb-Wu solutions as well as with those obtained with other state-of-the-art approximations. The structure of the intrinsic symmetry-broken determinants resulting from the variational procedure is interpreted in terms of solitons whose translational and breathing motions can be regarded as basic units of quantum fluctuations. It is also shown that in the case of doped one-dimensional lattices, a part of such fluctuations can also be interpreted in terms of polarons. In addition to momentum distributions, both spin-spin and density-density correlation functions are studied as functions of doping. The spectral functions and density of states, computed with an ansatz whose quality can be well controlled by the number of symmetry-projected configurations used to approximate the ${N}_{e}\ifmmode\pm\else\textpm\fi{}1$ electron systems, display features beyond a simple quasiparticle distribution, as well as spin-charge separation trends.

arXiv:1402.2968•Physical Review B•May 8, 2014•8 citations

Variational description of the ground state of the repulsive two-dimensional Hubbard model in terms of nonorthogonal symmetry-projected Slater determinants

R. R. Rodríguez-Guzmán, Carlos A. Jiménez-Hoyos, Gustavo E. Scuseria

The few-determinant (FED) methodology, introduced in our previous works [R. Rodr\'{\i}guez-Guzm\'an et al., Phys. Rev. B 87, 235129 (2013); Phys. Rev. B 89, 195109 (2014)] for one-dimensional (1D) lattices, is here adapted for the repulsive two-dimensional Hubbard model at half filling and with finite doping fractions. Within this configuration mixing scheme, a given ground state with well-defined spin and space group quantum numbers is expanded in terms of a nonorthogonal symmetry-projected basis determined through chains of variation-after-projection calculations. The results obtained for the ground-state and correlation energies of half-filled and doped $4\ifmmode\times\else\texttimes\fi{}4,\phantom{\rule{0.16em}{0ex}}6\ifmmode\times\else\texttimes\fi{}6,\phantom{\rule{0.16em}{0ex}}8\ifmmode\times\else\texttimes\fi{}8$, and $10\ifmmode\times\else\texttimes\fi{}10$ lattices, as well as momentum distributions and spin-spin correlation functions in small lattices, compare well with those obtained using other state-of-the-art approximations. The structure of the intrinsic determinants resulting from the variational strategy is interpreted in terms of defects that encode information on the basic units of quantum fluctuations in the considered 2D systems. The varying nature of the underlying quantum fluctuations, reflected in a transition to a stripe regime for increasing on-site repulsions, is discussed using the intrinsic determinants belonging to a $16\ifmmode\times\else\texttimes\fi{}4$ lattice with 56 electrons. Such a transition is further illustrated by computing spin-spin and charge-charge correlation functions with the corresponding multireference FED wave functions. In good agreement with previous studies, the analysis of the pairing correlation functions reveals a weak enhancement of the extended $s$-wave and ${d}_{{x}^{2}\ensuremath{-}{y}^{2}}$ pairing modes. Given the quality of results here reported together with those previously obtained for 1D lattices and the parallelization properties of the FED scheme, we believe that symmetry projection techniques are very well suited for building ground-state wave functions of correlated electronic systems, regardless of their dimensionality.

arXiv:1408.2567•Physical Review B•Nov 7, 2014•15 citations

Variational ground states of the two-dimensional Hubbard model

D. Baeriswyl, David Eichenberger, M. Z. Menteshashvili

Recent refinements of analytical and numerical methods have improved our understanding of the ground-state phase diagram of the two-dimensional (2D) Hubbard model. Here, we focus on variational approaches, but comparisons with both quantum cluster and Gaussian Monte Carlo methods are also made. Our own ansatz leads to an antiferromagnetic ground state at half filling with a slightly reduced staggered order parameter (as compared to simple mean-field theory). Away from half filling, we find d-wave superconductivity, but confined to densities where the Fermi surface passes through the antiferromagnetic zone boundary (if hopping between both nearest-neighbour and next-nearest-neighbour sites is considered). Our results agree surprisingly well with recent numerical studies using the quantum cluster method. An interesting trend is found by comparing gap parameters Δ (antiferromagnetic or superconducting) obtained with different variational wave functions. Δ varies by an order of magnitude and thus cannot be taken as a characteristic energy scale. In contrast, the order parameter is much less sensitive to the degree of sophistication of the variational schemes, at least at and near half filling.

arXiv:0907.1593•New Journal of Physics•Jul 23, 2009•45 citations

Spectral functions in the hubbard model with half-filling

Sergei G. Ovchinnikov, Elena I. Shneyder

Under the assumption of long-range antiferromagnetic order at low temperatures, the spectral functions and the density of states are calculated in the two-dimensional Hubbard model with half-filling in the Hubbard-I approximation. The results are compared with the data obtained using an exact numerical technique, namely, the quantum Monte Carlo method. The influence of hopping to the next-to-nearest neighbor on the formation of the electronic structure is considered.

Physics of the Solid State•Aug 1, 2004•12 citations
We present a symmetry-projected configuration mixing scheme to describe ground and excited states, with well defined quantum numbers, of the two-dimensional Hubbard model with nearest-neighbor hopping and periodic boundary conditions. Results for the half-filled 2\ifmmode×\else\texttimes\fi42\ifmmode\times\else\texttimes\fi{}42\ifmmode×\else\texttimes\fi4, 4\ifmmode×\else\texttimes\fi44\ifmmode\times\else\texttimes\fi{}44\ifmmode×\else\texttimes\fi4, and 6\ifmmode×\else\texttimes\fi66\ifmmode\times\else\texttimes\fi{}66\ifmmode×\else\texttimes\fi6 lattices as well as doped 4\ifmmode×\else\texttimes\fi44\ifmmode\times\else\texttimes\fi{}44\ifmmode×\else\texttimes\fi4 systems, compare well with available results, both exact and from other state-of-the-art approximations. We report spectral functions and density of states obtained from a well-controlled ansatz for the (Ne\ifmmode±\else\textpm\fi1)({N}_{e}\ifmmode\pm\else\textpm\fi{}1)(Ne​\ifmmode±\else\textpm\fi1)-electron system. Symmetry projected methods have been widely used for the many-body nuclear physics problem but have received little attention in the solid state community. Given their relatively low (mean-field) computational cost and the high quality of results here reported, we believe that they deserve further scrutiny.
We present a multireference configuration mixing scheme for describing ground and excited states, with well-defined spin and space-group symmetry quantum numbers, of the one-dimensional Hubbard model with nearest-neighbor hopping and periodic boundary conditions. Within this scheme, each state is expanded in terms of nonorthogonal and variationally determined symmetry-projected configurations. The results for lattices up to 30 and 50 sites compare well with the exact Lieb-Wu solutions as well as with results from other state-of-the-art approximations. In addition to spin-spin correlation functions in real space and magnetic structure factors, we present results for spectral functions and density of states computed with an Ansatz whose quality can be well controlled by the number of symmetry-projected configurations used to approximate the systems with Ne{N}_{e}Ne​ and ${N}_{e}\ifmmode\pm\else\textpm\fi{}1$ electrons. The intrinsic symmetry-broken determinants resulting from the variational calculations have rich structures in terms of defects that can be regarded as basic units of quantum fluctuations. Given the quality of the results here reported, as well as the parallelization properties of the considered scheme, we believe that symmetry-projection techniques, which have found ample applications in nuclear structure physics, deserve further attention in the study of low-dimensional correlated many-electron systems.
We apply a variational method devised for the nuclear many-body problem to the one-dimensional Hubbard model with nearest neighbor hopping and periodic boundary conditions. The test wave function consist for each state out of a single Hartree--Fock determinant mixing all the sites (or momenta) as well as the spin projections of the electrons. Total spin and linear momentum are restored by projection methods before the variation. It is demonstrated that this approach reproduces the results of exact diagonalizations for half-filled N=12N=12N=12 and N=14N=14N=14 lattices not only for the energies and occupation numbers of the ground but also of the lowest excited states rather well. Furthermore, a system of ten electrons in an N=12N=12N=12 lattice is investigated and, finally, an N=30N=30N=30 lattice is studied. In addition to energies and occupation numbers we present the spectral functions computed with the help of the symmetry-projected wave functions as well.
The few determinant (FED) approximation introduced in our previous work [Phys. Rev. B 87, 235129 (2013)] is used to describe the ground state, characterized by well-defined spin and space group symmetry quantum numbers as well as doping fractions Ne/Nsites{N}_{e}/{N}_{\mathrm{sites}}Ne​/Nsites​, of one-dimensional Hubbard lattices with nearest-neighbor hopping and periodic boundary conditions. Within this multireference scheme, each ground state is expanded in a given number of nonorthogonal and variationally determined symmetry-projected configurations. The results obtained for the ground-state and correlation energies of half-filled and doped lattices with 30, 34, and 50 sites compare well with the exact Lieb-Wu solutions as well as with those obtained with other state-of-the-art approximations. The structure of the intrinsic symmetry-broken determinants resulting from the variational procedure is interpreted in terms of solitons whose translational and breathing motions can be regarded as basic units of quantum fluctuations. It is also shown that in the case of doped one-dimensional lattices, a part of such fluctuations can also be interpreted in terms of polarons. In addition to momentum distributions, both spin-spin and density-density correlation functions are studied as functions of doping. The spectral functions and density of states, computed with an ansatz whose quality can be well controlled by the number of symmetry-projected configurations used to approximate the ${N}_{e}\ifmmode\pm\else\textpm\fi{}1$ electron systems, display features beyond a simple quasiparticle distribution, as well as spin-charge separation trends.
The few-determinant (FED) methodology, introduced in our previous works [R. Rodr\'{\i}guez-Guzm\'an et al., Phys. Rev. B 87, 235129 (2013); Phys. Rev. B 89, 195109 (2014)] for one-dimensional (1D) lattices, is here adapted for the repulsive two-dimensional Hubbard model at half filling and with finite doping fractions. Within this configuration mixing scheme, a given ground state with well-defined spin and space group quantum numbers is expanded in terms of a nonorthogonal symmetry-projected basis determined through chains of variation-after-projection calculations. The results obtained for the ground-state and correlation energies of half-filled and doped $4\ifmmode\times\else\texttimes\fi{}4,\phantom{\rule{0.16em}{0ex}}6\ifmmode\times\else\texttimes\fi{}6,\phantom{\rule{0.16em}{0ex}}8\ifmmode\times\else\texttimes\fi{}8$, and $10\ifmmode\times\else\texttimes\fi{}10$ lattices, as well as momentum distributions and spin-spin correlation functions in small lattices, compare well with those obtained using other state-of-the-art approximations. The structure of the intrinsic determinants resulting from the variational strategy is interpreted in terms of defects that encode information on the basic units of quantum fluctuations in the considered 2D systems. The varying nature of the underlying quantum fluctuations, reflected in a transition to a stripe regime for increasing on-site repulsions, is discussed using the intrinsic determinants belonging to a $16\ifmmode\times\else\texttimes\fi{}4$ lattice with 56 electrons. Such a transition is further illustrated by computing spin-spin and charge-charge correlation functions with the corresponding multireference FED wave functions. In good agreement with previous studies, the analysis of the pairing correlation functions reveals a weak enhancement of the extended sss-wave and ${d}_{{x}^{2}\ensuremath{-}{y}^{2}}$ pairing modes. Given the quality of results here reported together with those previously obtained for 1D lattices and the parallelization properties of the FED scheme, we believe that symmetry projection techniques are very well suited for building ground-state wave functions of correlated electronic systems, regardless of their dimensionality.
Recent refinements of analytical and numerical methods have improved our understanding of the ground-state phase diagram of the two-dimensional (2D) Hubbard model. Here, we focus on variational approaches, but comparisons with both quantum cluster and Gaussian Monte Carlo methods are also made. Our own ansatz leads to an antiferromagnetic ground state at half filling with a slightly reduced staggered order parameter (as compared to simple mean-field theory). Away from half filling, we find d-wave superconductivity, but confined to densities where the Fermi surface passes through the antiferromagnetic zone boundary (if hopping between both nearest-neighbour and next-nearest-neighbour sites is considered). Our results agree surprisingly well with recent numerical studies using the quantum cluster method. An interesting trend is found by comparing gap parameters Δ (antiferromagnetic or superconducting) obtained with different variational wave functions. Δ varies by an order of magnitude and thus cannot be taken as a characteristic energy scale. In contrast, the order parameter is much less sensitive to the degree of sophistication of the variational schemes, at least at and near half filling.