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. B, Condensed matter

Electron correlations in the two-dimensional Hubbard model: A group-theoretical and numerical study

G. Fano, F. Ortolani, Alberto Parola

Published Jul 1, 1992• 53 citations
Publisher

Abstract

The Lanczos algorithm has been applied to the Hubbard model on a 4\ifmmode\times\else\texttimes\fi{}4 square lattice and several equal-time correlation functions have been evaluated. This computation has been made possible by a classification of all the irreducible representations of the space symmetry group, which allowed a considerable reduction of the dimension of the Hilbert space of the system. Such an analysis is quite general and might be useful in other physical applications. Our results shed light on the effects of doping on the magnetic ordering, suggest a spin-charge decoupling at large U, and show a remarkable reduction of the quasiparticle weight at strong coupling.

DOI
10.1103/physrevb.46.1048

Related papers

Hybrid-space density matrix renormalization group study of the doped two-dimensional Hubbard model

G. Ehlers, Steven R. White, R. M. Noack

The performance of the density matrix renormalization group (DMRG) is strongly influenced by the choice of the local basis of the underlying physical lattice. We demonstrate that, for the two-dimensional Hubbard model, the hybrid--real-momentum-space formulation of the DMRG is computationally more efficient than the standard real-space formulation. In particular, we show that the computational cost for fixed bond dimension of the hybrid-space DMRG is approximately independent of the width of the lattice, in contrast to the real-space DMRG, for which it is proportional to the width squared. We apply the hybrid-space algorithm to calculate the ground state of the doped two-dimensional Hubbard model on cylinders of width four and six sites; at $n=0.875$ filling, the ground state exhibits a striped charge-density distribution with a wavelength of eight sites for both $U/t=4.0$ and $8.0$. We find that the strength of the charge ordering depends on $U/t$ and on the boundary conditions. Furthermore, we investigate the magnetic ordering as well as the decay of the static spin, charge, and pair-field correlation functions.

arXiv:1701.03690•Physical review. B./Physical review. B•Mar 21, 2017•90 citations

Two-dimensional Hubbard model at finite temperature: Weak, strong, and long correlation regimes

Fedor Šimkovic, Riccardo Rossi, Michel Ferrero

We investigate the momentum-resolved spin and charge susceptibilities, as well as the chemical potential and double occupancy in the two-dimensional Hubbard model as functions of doping, temperature, and interaction strength. Through these quantities, we identify a weak-coupling regime, a strong-coupling regime with short-range correlations and an intermediate-coupling regime with long magnetic correlation lengths. In the spin channel, we observe an additional crossover from commensurate to incommensurate correlations. In contrast, we find charge correlations to be only short ranged for all studied temperatures, which suggests that the spin and charge responses are decoupled. These findings were obtained by a connected determinant diagrammatic Monte Carlo algorithm for the computation of double expansions, which we introduce in this paper. This permits us to obtain numerically exact results at temperatures as low as $T=0.067$ and interactions up to $U=7$, while working on arbitrarily large lattices. Our method also allows us to gain physical insights from investigating the analytic structure of perturbative series. We connect to previous work by studying smaller lattice geometries and report substantial finite-size effects.

arXiv:2110.05863•Physical Review Research•Dec 20, 2022•45 citations

Ground state of the strong-coupling Hubbard Hamiltonian: A numerical diagonalization study

Efthimios Kaxiras, Efstratios Manousakis

We exactly diagonalize the effective Hamiltonian obtained from the Hubbard model in the strong-coupling limit for a two-dimensional $\sqrt{10}\ifmmode\times\else\texttimes\fi{}\sqrt{10}$ size square lattice. The effective Hamiltonian operates on a restricted Hilbert space containing only states with singly occupied sites, which makes the diagonalization possible within reasonable computational time for the above lattice. The ground-state energy and wave function are obtained for several values of the coupling ratio $\frac{t}{U}$ and the doping fraction $x$. We find three different phases in the ($x,\frac{t}{U}$) phase diagram. The first is characterized by antiferromagnetic order which extends to the longest possible distance in the 10-site system. The second and third phases are characterized by antiferromagnetic and ferromagnetic short-range correlations, respectively. We comment on the possible relevance of the results to the recently discovered high-temperature copper-oxide superconductors.

Physical review. B, Condensed matter•Jan 1, 1988•53 citations

The Study of Model Hamiltonians on Finite Clusters.

Duan Chen

We present two methods to study model Hamiltonians defined on clusters. One is the symmetrization of the Hilbert space basis according to irreducible representations of the symmetry group of the Hamiltonian. The other one adapts the Lanczos diagonalization method directly to the second quantized form of the Hamiltonian studied. By the symmetrization of the basis, one takes advantage of the symmetries of the Hamiltonian, and uses the group theory to find the invariant subspaces of the Hamiltonian in the complete Hilbert space. Then one needs only to work within the invariant subspaces. This method cuts down substantially the size of the matrix representation of the Hamiltonian to be diagonalized. Therefore, one could study relatively larger systems. As a result, all the eigenvectors and eigenvalues are obtained. Consequently, one could calculate thermodynamical quantities and various correlation functions, such as, the spin correlation, the charge correlation, the superconductivity pairing correlation, and the photoemission spectra. One the other hand, the Lanczos method, mentioned here in particular, is a specially developed computer program for the study of the model Hamiltonians. This program uses the key feature of the Lanczos method in diagonalization, allows one to use the second quantized form of the Hamiltonian directly without pregeneration or storing the matrix elements of the Hamiltonian. Two methods could be combined to take the full advantage of the symmetries of the Hamiltonian by input of the symmetrized wavefunction as the starting trial wavefunction in the Lanczos method. Here we illustrate the above two methods by applying to the Hubbard, and extended Hubbard models on a simple cubic cluster, and the lattice Anderson model on the two dimensional square lattice. Other interesting results for different systems are also briefly discussed here.

Civil War Book Review•Jan 1, 1989•0 citations

Exact-diagonalization method for correlated-electron models

Canio Noce, Mario Cuoco

We present a method to compute eigenstates and eigenvalues of strongly correlated models on clusters. The calculation is performed applying the SU(2) symmetry, which allows a considerable reduction of the dimension of Hilbert space of the system. We illustrate the method classifying the states of the four-sites Hubbard model on a chain. This classification has been made using also the translational symmetry and the S\ifmmode \tilde{}\else \~{}\fi{}\ifmmode \tilde{U}\else \~{U}\fi{}(2) pseudospin symmetry. \textcopyright{} 1996 The American Physical Society.

Physical review. B, Condensed matter•Nov 1, 1996•22 citations

Properties of strongly correlated one dimensional chains beyond the Hubbard model

Samuel Milner

The Hubbard model is a simple but profound theoretical tool for understanding the physics underlying strongly correlated systems. This phenomenology includes, but is not limited to, metal-insulating transitions, charge and spin order (charge density/spin density waves), superconducting phases, non-Fermi liquid physics, among others. Examining one dimensional physical models described by the Hubbard Hamiltonian can teach us much about the related physics of higher dimensional systems. Applying advanced numerical computational techniques such as density matrix numerical renormalization group (DMRG) to one-dimensional manypbody model Hamiltonians allows for the calculation of ground state energies and dynamical spectral functions with extremely high precision and accuracy. This thesis explores the application of DMRG numerical techniques to one-dimensional physical systems described by Hamiltonians beyond the conventional Hubbard model, and the nuanced physical effects observed for these systems across various regimes and phases. A Wannier Hubbard chain composed of non-orthogonal overlapping atomic orbitals was transformed to an orthogonal atomic basis using a non-unitary Lowdin transformation. The calculation of energies, dynamical spectral functions and the photoemission spectrum for this Wannier Hubbard chain revealed phases not previously observed for Hubbard chains at half-filling. Additional analysis and calculations were performed for a CuO chain described by the Luther-Emery Hamiltonian that epitomizes charge-transfer physics. The motivation for studying this CuO chain was previous literature on transition metal oxides (TMOs) and the mapping of multiband Hamiltonians to single band models in certain energy regimes. Finally, we performed computational analysis for a spin-3/2 generalized Hubbard chain. By spanning the full phase space of our spin-3/2 Hubbard chain in terms of chemical potential $\mu$ and external magnetic field $h$ we observed non-typical Cooper pairing states and novel phases for different symmetries of our Hamiltonian.--Author's abstract

•Jan 1, 2025•0 citations
The Lanczos algorithm has been applied to the Hubbard model on a 4\ifmmode\times\else\texttimes\fi{}4 square lattice and several equal-time correlation functions have been evaluated. This computation has been made possible by a classification of all the irreducible representations of the space symmetry group, which allowed a considerable reduction of the dimension of the Hilbert space of the system. Such an analysis is quite general and might be useful in other physical applications. Our results shed light on the effects of doping on the magnetic ordering, suggest a spin-charge decoupling at large U, and show a remarkable reduction of the quasiparticle weight at strong coupling.
The performance of the density matrix renormalization group (DMRG) is strongly influenced by the choice of the local basis of the underlying physical lattice. We demonstrate that, for the two-dimensional Hubbard model, the hybrid--real-momentum-space formulation of the DMRG is computationally more efficient than the standard real-space formulation. In particular, we show that the computational cost for fixed bond dimension of the hybrid-space DMRG is approximately independent of the width of the lattice, in contrast to the real-space DMRG, for which it is proportional to the width squared. We apply the hybrid-space algorithm to calculate the ground state of the doped two-dimensional Hubbard model on cylinders of width four and six sites; at n=0.875n=0.875n=0.875 filling, the ground state exhibits a striped charge-density distribution with a wavelength of eight sites for both U/t=4.0U/t=4.0U/t=4.0 and 8.08.08.0. We find that the strength of the charge ordering depends on U/tU/tU/t and on the boundary conditions. Furthermore, we investigate the magnetic ordering as well as the decay of the static spin, charge, and pair-field correlation functions.
We investigate the momentum-resolved spin and charge susceptibilities, as well as the chemical potential and double occupancy in the two-dimensional Hubbard model as functions of doping, temperature, and interaction strength. Through these quantities, we identify a weak-coupling regime, a strong-coupling regime with short-range correlations and an intermediate-coupling regime with long magnetic correlation lengths. In the spin channel, we observe an additional crossover from commensurate to incommensurate correlations. In contrast, we find charge correlations to be only short ranged for all studied temperatures, which suggests that the spin and charge responses are decoupled. These findings were obtained by a connected determinant diagrammatic Monte Carlo algorithm for the computation of double expansions, which we introduce in this paper. This permits us to obtain numerically exact results at temperatures as low as T=0.067T=0.067T=0.067 and interactions up to U=7U=7U=7, while working on arbitrarily large lattices. Our method also allows us to gain physical insights from investigating the analytic structure of perturbative series. We connect to previous work by studying smaller lattice geometries and report substantial finite-size effects.
We exactly diagonalize the effective Hamiltonian obtained from the Hubbard model in the strong-coupling limit for a two-dimensional $\sqrt{10}\ifmmode\times\else\texttimes\fi{}\sqrt{10}$ size square lattice. The effective Hamiltonian operates on a restricted Hilbert space containing only states with singly occupied sites, which makes the diagonalization possible within reasonable computational time for the above lattice. The ground-state energy and wave function are obtained for several values of the coupling ratio tU\frac{t}{U}Ut​ and the doping fraction xxx. We find three different phases in the (x,tUx,\frac{t}{U}x,Ut​) phase diagram. The first is characterized by antiferromagnetic order which extends to the longest possible distance in the 10-site system. The second and third phases are characterized by antiferromagnetic and ferromagnetic short-range correlations, respectively. We comment on the possible relevance of the results to the recently discovered high-temperature copper-oxide superconductors.
We present two methods to study model Hamiltonians defined on clusters. One is the symmetrization of the Hilbert space basis according to irreducible representations of the symmetry group of the Hamiltonian. The other one adapts the Lanczos diagonalization method directly to the second quantized form of the Hamiltonian studied. By the symmetrization of the basis, one takes advantage of the symmetries of the Hamiltonian, and uses the group theory to find the invariant subspaces of the Hamiltonian in the complete Hilbert space. Then one needs only to work within the invariant subspaces. This method cuts down substantially the size of the matrix representation of the Hamiltonian to be diagonalized. Therefore, one could study relatively larger systems. As a result, all the eigenvectors and eigenvalues are obtained. Consequently, one could calculate thermodynamical quantities and various correlation functions, such as, the spin correlation, the charge correlation, the superconductivity pairing correlation, and the photoemission spectra. One the other hand, the Lanczos method, mentioned here in particular, is a specially developed computer program for the study of the model Hamiltonians. This program uses the key feature of the Lanczos method in diagonalization, allows one to use the second quantized form of the Hamiltonian directly without pregeneration or storing the matrix elements of the Hamiltonian. Two methods could be combined to take the full advantage of the symmetries of the Hamiltonian by input of the symmetrized wavefunction as the starting trial wavefunction in the Lanczos method. Here we illustrate the above two methods by applying to the Hubbard, and extended Hubbard models on a simple cubic cluster, and the lattice Anderson model on the two dimensional square lattice. Other interesting results for different systems are also briefly discussed here.
We present a method to compute eigenstates and eigenvalues of strongly correlated models on clusters. The calculation is performed applying the SU(2) symmetry, which allows a considerable reduction of the dimension of Hilbert space of the system. We illustrate the method classifying the states of the four-sites Hubbard model on a chain. This classification has been made using also the translational symmetry and the S\ifmmode \tilde{}\else \~{}\fi{}\ifmmode \tilde{U}\else \~{U}\fi{}(2) pseudospin symmetry. \textcopyright{} 1996 The American Physical Society.
The Hubbard model is a simple but profound theoretical tool for understanding the physics underlying strongly correlated systems. This phenomenology includes, but is not limited to, metal-insulating transitions, charge and spin order (charge density/spin density waves), superconducting phases, non-Fermi liquid physics, among others. Examining one dimensional physical models described by the Hubbard Hamiltonian can teach us much about the related physics of higher dimensional systems. Applying advanced numerical computational techniques such as density matrix numerical renormalization group (DMRG) to one-dimensional manypbody model Hamiltonians allows for the calculation of ground state energies and dynamical spectral functions with extremely high precision and accuracy. This thesis explores the application of DMRG numerical techniques to one-dimensional physical systems described by Hamiltonians beyond the conventional Hubbard model, and the nuanced physical effects observed for these systems across various regimes and phases. A Wannier Hubbard chain composed of non-orthogonal overlapping atomic orbitals was transformed to an orthogonal atomic basis using a non-unitary Lowdin transformation. The calculation of energies, dynamical spectral functions and the photoemission spectrum for this Wannier Hubbard chain revealed phases not previously observed for Hubbard chains at half-filling. Additional analysis and calculations were performed for a CuO chain described by the Luther-Emery Hamiltonian that epitomizes charge-transfer physics. The motivation for studying this CuO chain was previous literature on transition metal oxides (TMOs) and the mapping of multiband Hamiltonians to single band models in certain energy regimes. Finally, we performed computational analysis for a spin-3/2 generalized Hubbard chain. By spanning the full phase space of our spin-3/2 Hubbard chain in terms of chemical potential μ\muμ and external magnetic field hhh we observed non-typical Cooper pairing states and novel phases for different symmetries of our Hamiltonian.--Author's abstract