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:1402.2968 · Physical Review B

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

Published May 8, 2014• 8 citations
arXivHTMLPDF

Abstract

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.

DOI
10.1103/physrevb.89.195109

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

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

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 +1 more

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.

arXiv:1204.2006•Physical Review B•Jun 26, 2012•48 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

Exact and variational investigations of Hubbard rings

William Hodge

We study the one-dimensional Hubbard model for interacting electrons in a crystal lattice, using two methods: solution of the Lieb-Wu equations in the atomic limit of no intersite hopping, and variational techniques to obtain a lower bound for the ground state energy. From the Lieb-Wu equations we derive an exact expression for the ground state energy in the atomic limit. We obtain our result by exploiting the spin symmetry of the Hubbard Hamiltonian. Our energy is lower than that obtained by previous published results, and we are able to explain why our results differ. Our method also gives information about the magnetic ordering of the ground state as a function of the number of electrons. Our results are valid for arbitrary electron density values. The variational techniques we use are applied to the two-particle reduced-density matrix (2-RDM), because this quantity suffices to obtain expectation values of all one- and two-particle operators, and because it is less computationally intensive than obtaining the complete N -particle wave function. Because the 2-RDM must be positive semidefinite, we use an efficient programming algorithm known as the Semidefinite Programming Algorithm. A naive search through the space of two-particle matrices encounters a difficult problem, known as the " N -representability problem", which is to assure that the trial density matrices actually come from a wave function that has the required exchange symmetry. This problem is ameliorated by imposing constraints on the variational search that eliminate trial matrices with the incorrect exchange symmetry. These constraints express known physical properties of the ground state (for example its total spin and its translational symmetry) and the so-called D, Q , and G two-positivity conditions that express the required exchange symmetry. To calculate matrix elements, we use the two-particle eigenstates of the reduced Hubbard Hamiltonian; this choice of basis significantly decreases the number of matrix elements needed for energy calculations. After an optimal two-particle density matrix has been obtained, we calculate additional properties of the ground state, including spin correlation functions, chemical potentials, and the momentum distribution. We compare all of these properties to values obtained by exact numerical solution of the Schrödinger equation.

•Jan 1, 2008•0 citations

Visualization of quantum fluctuations by superposition of optimized nonorthogonal Slater determinants

Norikazu Tomita

Electronic states in the one-dimensional (1D) doped Hubbard model are described by superposition of optimized nonorthogonal Slater determinants (S-dets). Analysis on the S-dets allows us to visualize quantum fluctuations. In the weak and intermediate interaction regimes, quantum fluctuations are described by translation and breathing motions of spin-charge coupled defects called polarons as well as spin-charge decoupled defects called holons and spinons. In the strong-interaction regime, on the other hand, spin and charge fluctuations are mostly separated as spinons and holons, especially in the lightly doped systems $(\ensuremath{\delta}=0.08)$. In the highly doped systems $(\ensuremath{\delta}=0.24)$, polarons also play an important role. Finally, it is shown that in the 1/3-filled systems, the concepts of holon, spinon, and polaron do not work anymore. The electronic structure is qualitatively described by mixtures of two different density waves having triple periodicity to the lattice. The domains of these different density waves make quantum fluctuations.

Physical Review B•Feb 13, 2009•8 citations
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 Ne\ifmmode±\else\textpm\fi1{N}_{e}\ifmmode\pm\else\textpm\fi{}1Ne​\ifmmode±\else\textpm\fi1 electron systems, display features beyond a simple quasiparticle distribution, as well as spin-charge separation trends.
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.
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.
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.
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.
We study the one-dimensional Hubbard model for interacting electrons in a crystal lattice, using two methods: solution of the Lieb-Wu equations in the atomic limit of no intersite hopping, and variational techniques to obtain a lower bound for the ground state energy. From the Lieb-Wu equations we derive an exact expression for the ground state energy in the atomic limit. We obtain our result by exploiting the spin symmetry of the Hubbard Hamiltonian. Our energy is lower than that obtained by previous published results, and we are able to explain why our results differ. Our method also gives information about the magnetic ordering of the ground state as a function of the number of electrons. Our results are valid for arbitrary electron density values. The variational techniques we use are applied to the two-particle reduced-density matrix (2-RDM), because this quantity suffices to obtain expectation values of all one- and two-particle operators, and because it is less computationally intensive than obtaining the complete N -particle wave function. Because the 2-RDM must be positive semidefinite, we use an efficient programming algorithm known as the Semidefinite Programming Algorithm. A naive search through the space of two-particle matrices encounters a difficult problem, known as the " N -representability problem", which is to assure that the trial density matrices actually come from a wave function that has the required exchange symmetry. This problem is ameliorated by imposing constraints on the variational search that eliminate trial matrices with the incorrect exchange symmetry. These constraints express known physical properties of the ground state (for example its total spin and its translational symmetry) and the so-called D, Q , and G two-positivity conditions that express the required exchange symmetry. To calculate matrix elements, we use the two-particle eigenstates of the reduced Hubbard Hamiltonian; this choice of basis significantly decreases the number of matrix elements needed for energy calculations. After an optimal two-particle density matrix has been obtained, we calculate additional properties of the ground state, including spin correlation functions, chemical potentials, and the momentum distribution. We compare all of these properties to values obtained by exact numerical solution of the Schrödinger equation.
Electronic states in the one-dimensional (1D) doped Hubbard model are described by superposition of optimized nonorthogonal Slater determinants (S-dets). Analysis on the S-dets allows us to visualize quantum fluctuations. In the weak and intermediate interaction regimes, quantum fluctuations are described by translation and breathing motions of spin-charge coupled defects called polarons as well as spin-charge decoupled defects called holons and spinons. In the strong-interaction regime, on the other hand, spin and charge fluctuations are mostly separated as spinons and holons, especially in the lightly doped systems $(\ensuremath{\delta}=0.08)$. In the highly doped systems $(\ensuremath{\delta}=0.24)$, polarons also play an important role. Finally, it is shown that in the 1/3-filled systems, the concepts of holon, spinon, and polaron do not work anymore. The electronic structure is qualitatively described by mixtures of two different density waves having triple periodicity to the lattice. The domains of these different density waves make quantum fluctuations.