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:2602.03771

Spin and Charge Conductivity in the Square Lattice Fermi-Hubbard Model

Linh Pham, Ehsan Khatami

Published Feb 3, 2026• 0 citations
arXivHTMLPDF

Abstract

Dynamical properties are notoriously difficult to compute in numerical treatments of the Fermi-Hubbard model, especially in two spatial dimensions. However, they are essential in providing us with insight into some of the most important and less well-understood phases of the model, such as the pseudogap and strange metal phases at relatively high temperatures, or unconventional superconductivity at lower temperatures, away from the commensurate filling. Here, we use the numerical linked-cluster expansions to compute spin and charge optical conductivities of the model at different temperatures and strong interaction strengths via the exact real-time-dependent correlation functions of the current operators. We mitigate systematic errors associated with having a limited access to the long-time behavior of the correlators by introducing fits and allowing for non-zero Drude weights when appropriate. We compare our results to available data from optical lattice experiments and find that the Drude contributions can account for the theory-experiment gap in the DC spin conductivity of the model at half filling in the strong-coupling region. Our method helps paint a more complete picture of the conductivity in the two-dimensional Hubbard model and opens the door to studying dynamical properties of quantum lattice models in the thermodynamic limit.

DOI
10.48550/arxiv.2602.03771

Related papers

Transport properties of the one-dimensional Hubbard model at finite temperature

Christoph Karrasch, Dante Marvin Kennes, Joel E. Moore

We study finite-temperature transport properties of the one-dimensional Hubbard model using the density matrix renormalization group. Our aim is twofold: First, we compute both the charge- and the spin-current correlation function of the integrable model at half filling. The former decays rapidly, implying that the corresponding Drude weight is either zero or very small. Second, we calculate the optical charge conductivity ${\ensuremath{\sigma}}_{r}eg(\ensuremath{\omega})$ in the presence of small integrability-breaking next-nearest-neighbor interactions (the extended Hubbard model). The dc conductivity is finite and diverges as the temperature is decreased below the gap. Our results thus suggest that the half-filled, gapped Hubbard model is a normal charge conductor at finite temperatures. As a test bed for our numerics, we compute ${\ensuremath{\sigma}}_{r}eg(\ensuremath{\omega})$ for the integrable XXZ spin chain in its gapped phase.

arXiv:1408.1891•Physical Review B•Oct 3, 2014•56 citations

Cluster dynamical mean field theory study of antiferromagnetic transition in the square-lattice Hubbard model: Optical conductivity and electronic structure

Toshihiro Sato, Hirokazu Tsunetsugu

We numerically study optical conductivity $\ensuremath{\sigma}(\ensuremath{\omega})$ near the ``antiferromagnetic'' phase transition in the square-lattice Hubbard model at half filling. We use a cluster dynamical mean field theory and calculate conductivity including vertex corrections and, to this end, we have reformulated the vertex corrections in the antiferromagnetic phase. We find that the vertex corrections change various important details in temperature and $\ensuremath{\omega}$ dependencies of conductivity in the square lattice, and this contrasts sharply the case of the Mott transition in the frustrated triangular lattice. Generally, the vertex corrections enhance variations in the $\ensuremath{\omega}$ dependence, and sharpen the Drude peak and a high-$\ensuremath{\omega}$ incoherent peak in the paramagnetic phase. They also enhance the dip in $\ensuremath{\sigma}(\ensuremath{\omega})$ at $\ensuremath{\omega}=0$ in the antiferromagnetic phase. Therefore, the dc conductivity is enhanced in the paramagnetic phase and suppressed in the antiferromagnetic phase, but this change occurs slightly below the transition temperature. We also find a temperature region above the transition temperature in which the dc conductivity shows an insulating behavior but $\ensuremath{\sigma}(\ensuremath{\omega})$ retains the Drude peak, and this region is stabilized by the vertex corrections. We also investigate which fluctuations are important in the vertex corrections and analyze momentum dependence of the vertex function in detail.

Physical review. B./Physical review. B•Aug 8, 2016•18 citations

Publisher's Note: Cluster dynamical mean field theory study of antiferromagnetic transition in the square-lattice Hubbard model: Optical conductivity and electronic structure [Phys. Rev. B94, 085110 (2016)]

Toshihiro Sato, Hirokazu Tsunetsugu

We numerically study optical conductivity $σ(ω)$ near the "antiferromagnetic" phase transition in the square-lattice Hubbard model at half filling. We use a cluster dynamical mean field theory and calculate conductivity including vertex corrections, and to this end, we have reformulated the vertex corrections in the antiferromagnetic phase. We find that the vertex corrections change various important details in temperature- and $ω$-dependencies of conductivity in the square lattice, and this contrasts sharply the case of the Mott transition in the frustrated triangular lattice. Generally, the vertex corrections enhance variations in the $ω$-dependence, and sharpen the Drude peak and a high-$ω$ incoherent peak in the paramagnetic phase. They also enhance the dip in $σ(ω)$ at $ω$=0 in the antiferromagnetic phase. Therefore, the dc conductivity is enhanced in the paramagnetic phase and suppressed in the antiferromagnetic phase, but this change occurs slightly below the transition temperature. We also find a temperature region above the transition temperature in which the dc conductivity shows an insulating behavior but $σ(ω)$ retains the Drude peak, and this region is stabilized by the vertex corrections. We also investigate which fluctuations are important in the vertex corrections and analyze momentum dependence of the vertex function in detail.

arXiv:1605.00387•Physical review. B./Physical review. B•Aug 22, 2016•4 citations

Dynamic conductivity of strongly correlated electrons: The Hubbard model on a cubic lattice

L. Tan, J. Callaway

The frequency-dependent (optical) conductivity of the Hubbard model on a three-dimensional cubic lattice is studied through calculations on a small cluster. The interaction strength is varied between the weak- and strong-coupling limits. Results are reported for the half-filled band, and for small hole dopings away from half filling. The formation of a Hubbard gap is observed. Results are related to aspects of the current discussion of the metal-insulator transition in bulk systems.

Physical review. B, Condensed matter•Sep 1, 1992•12 citations

Effect of frustration on charge dynamics for a doped two-dimensional triangular Hubbard lattice: Comparison with a square lattice

Takami Tohyama

We examine the optical conductivity $\ensuremath{\sigma}(\ensuremath{\omega})$ and the chemical potential $\ensuremath{\mu}$, together with the spin correlation, in the strong-coupling limit of a hole-doped two-dimensional triangular Hubbard model near half filling by using an exact diagonalization technique. In contrast to the case of a square lattice without frustration, the doping dependences of $\ensuremath{\mu}$ and the Drude weight indicate that the charge degree of freedom is weakly coupled to the spin degree of freedom. However, we find that $\ensuremath{\sigma}(\ensuremath{\omega})$ shows strong incoherent excitations extended to a higher energy region. This implies that geometrical frustration in strongly correlated electron systems influences incoherent charge dynamics.

arXiv:cond-mat/0605007•Physical Review B•Sep 21, 2006•6 citations

Charge fluctuations, hydrodynamics, and transport in the square-lattice Hubbard model

J. Vučičević, Sonja Predin, Michel Ferrero

Recent experimental results suggest that a particular hydrodynamic theory describes charge fluctuations at long wavelengths in the square-lattice Hubbard model. Due to the continuity equation, the correlation functions for the charge and the current are directly connected: the parameters of the effective hydrodynamic model thus determine the optical conductivity. Here we investigate the validity of the proposed hydrodynamic theory in the full range of parameters of the Hubbard model. In the noninteracting case, there is no effective hydrodynamics, and the charge fluctuations present a rich variety of nonuniversal behaviors. At weak coupling, the optical conductivity is consistent with the hydrodynamic theory: at low frequency one observes a Lorentzian-shaped Drude peak, but the high-frequency asymptotics is necessarily different; the high-temperature limit for the product of the two hydrodynamic model parameters is also in agreement with numerical data. At strong coupling, we find that a generalization of the proposed hydrodynamic law is consistent with our quantum Monte Carlo, as well as the finite-temperature Lanczos results from literature. Most importantly, the temperature dependence of the hydrodynamic parameters as well as the dc resistivity are found to be very similar in the weak- and the strong-coupling regimes.

arXiv:2208.04047•Physical review. B./Physical review. B•Apr 25, 2023•14 citations
Dynamical properties are notoriously difficult to compute in numerical treatments of the Fermi-Hubbard model, especially in two spatial dimensions. However, they are essential in providing us with insight into some of the most important and less well-understood phases of the model, such as the pseudogap and strange metal phases at relatively high temperatures, or unconventional superconductivity at lower temperatures, away from the commensurate filling. Here, we use the numerical linked-cluster expansions to compute spin and charge optical conductivities of the model at different temperatures and strong interaction strengths via the exact real-time-dependent correlation functions of the current operators. We mitigate systematic errors associated with having a limited access to the long-time behavior of the correlators by introducing fits and allowing for non-zero Drude weights when appropriate. We compare our results to available data from optical lattice experiments and find that the Drude contributions can account for the theory-experiment gap in the DC spin conductivity of the model at half filling in the strong-coupling region. Our method helps paint a more complete picture of the conductivity in the two-dimensional Hubbard model and opens the door to studying dynamical properties of quantum lattice models in the thermodynamic limit.
We study finite-temperature transport properties of the one-dimensional Hubbard model using the density matrix renormalization group. Our aim is twofold: First, we compute both the charge- and the spin-current correlation function of the integrable model at half filling. The former decays rapidly, implying that the corresponding Drude weight is either zero or very small. Second, we calculate the optical charge conductivity ${\ensuremath{\sigma}}_{r}eg(\ensuremath{\omega})$ in the presence of small integrability-breaking next-nearest-neighbor interactions (the extended Hubbard model). The dc conductivity is finite and diverges as the temperature is decreased below the gap. Our results thus suggest that the half-filled, gapped Hubbard model is a normal charge conductor at finite temperatures. As a test bed for our numerics, we compute ${\ensuremath{\sigma}}_{r}eg(\ensuremath{\omega})$ for the integrable XXZ spin chain in its gapped phase.
We numerically study optical conductivity $\ensuremath{\sigma}(\ensuremath{\omega})$ near the ``antiferromagnetic'' phase transition in the square-lattice Hubbard model at half filling. We use a cluster dynamical mean field theory and calculate conductivity including vertex corrections and, to this end, we have reformulated the vertex corrections in the antiferromagnetic phase. We find that the vertex corrections change various important details in temperature and $\ensuremath{\omega}$ dependencies of conductivity in the square lattice, and this contrasts sharply the case of the Mott transition in the frustrated triangular lattice. Generally, the vertex corrections enhance variations in the $\ensuremath{\omega}$ dependence, and sharpen the Drude peak and a high-$\ensuremath{\omega}$ incoherent peak in the paramagnetic phase. They also enhance the dip in $\ensuremath{\sigma}(\ensuremath{\omega})$ at $\ensuremath{\omega}=0$ in the antiferromagnetic phase. Therefore, the dc conductivity is enhanced in the paramagnetic phase and suppressed in the antiferromagnetic phase, but this change occurs slightly below the transition temperature. We also find a temperature region above the transition temperature in which the dc conductivity shows an insulating behavior but $\ensuremath{\sigma}(\ensuremath{\omega})$ retains the Drude peak, and this region is stabilized by the vertex corrections. We also investigate which fluctuations are important in the vertex corrections and analyze momentum dependence of the vertex function in detail.
We numerically study optical conductivity σ(ω)σ(ω)σ(ω) near the "antiferromagnetic" phase transition in the square-lattice Hubbard model at half filling. We use a cluster dynamical mean field theory and calculate conductivity including vertex corrections, and to this end, we have reformulated the vertex corrections in the antiferromagnetic phase. We find that the vertex corrections change various important details in temperature- and ωωω-dependencies of conductivity in the square lattice, and this contrasts sharply the case of the Mott transition in the frustrated triangular lattice. Generally, the vertex corrections enhance variations in the ωωω-dependence, and sharpen the Drude peak and a high-ωωω incoherent peak in the paramagnetic phase. They also enhance the dip in σ(ω)σ(ω)σ(ω) at ωωω=0 in the antiferromagnetic phase. Therefore, the dc conductivity is enhanced in the paramagnetic phase and suppressed in the antiferromagnetic phase, but this change occurs slightly below the transition temperature. We also find a temperature region above the transition temperature in which the dc conductivity shows an insulating behavior but σ(ω)σ(ω)σ(ω) retains the Drude peak, and this region is stabilized by the vertex corrections. We also investigate which fluctuations are important in the vertex corrections and analyze momentum dependence of the vertex function in detail.
We examine the optical conductivity $\ensuremath{\sigma}(\ensuremath{\omega})$ and the chemical potential $\ensuremath{\mu}$, together with the spin correlation, in the strong-coupling limit of a hole-doped two-dimensional triangular Hubbard model near half filling by using an exact diagonalization technique. In contrast to the case of a square lattice without frustration, the doping dependences of $\ensuremath{\mu}$ and the Drude weight indicate that the charge degree of freedom is weakly coupled to the spin degree of freedom. However, we find that $\ensuremath{\sigma}(\ensuremath{\omega})$ shows strong incoherent excitations extended to a higher energy region. This implies that geometrical frustration in strongly correlated electron systems influences incoherent charge dynamics.
Recent experimental results suggest that a particular hydrodynamic theory describes charge fluctuations at long wavelengths in the square-lattice Hubbard model. Due to the continuity equation, the correlation functions for the charge and the current are directly connected: the parameters of the effective hydrodynamic model thus determine the optical conductivity. Here we investigate the validity of the proposed hydrodynamic theory in the full range of parameters of the Hubbard model. In the noninteracting case, there is no effective hydrodynamics, and the charge fluctuations present a rich variety of nonuniversal behaviors. At weak coupling, the optical conductivity is consistent with the hydrodynamic theory: at low frequency one observes a Lorentzian-shaped Drude peak, but the high-frequency asymptotics is necessarily different; the high-temperature limit for the product of the two hydrodynamic model parameters is also in agreement with numerical data. At strong coupling, we find that a generalization of the proposed hydrodynamic law is consistent with our quantum Monte Carlo, as well as the finite-temperature Lanczos results from literature. Most importantly, the temperature dependence of the hydrodynamic parameters as well as the dc resistivity are found to be very similar in the weak- and the strong-coupling regimes.