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:1608.01600 · Physical review. B./Physical review. B

Transport and optical conductivity in the Hubbard model: A high-temperature expansion perspective

Edward Perepelitsky, Andrew Galatas, Jernej Mravlje, Rok Žitko, Ehsan Khatami, B. Sriram Shastry, Antoine Georges

Published Dec 5, 2016• 57 citations
arXivHTMLPDF

Abstract

We derive analytical expressions for the spectral moments of the dynamical response functions of the Hubbard model using the high-temperature series expansion. We consider generic dimension $d$ as well as the infinite-$d$ limit, arbitrary electron density $n$, and both finite and infinite repulsion $U$. We use moment-reconstruction methods to obtain the one-electron spectral function, the self-energy, and the optical conductivity. They are all smooth functions at high temperature and, at large $U$, they are featureless with characteristic widths of the order of the lattice hopping parameter $t$. In the infinite-$d$ limit, we compare the series expansion results with accurate numerical renormalization group and interaction expansion quantum Monte Carlo results. We find excellent agreement down to surprisingly low temperatures, throughout most of the bad-metal regime, which applies for $T\ensuremath{\gtrsim}(1\ensuremath{-}n)D$, the Brinkman-Rice scale. The resistivity increases linearly in $T$ at high temperature without saturation. This results from the $1/T$ behavior of the compressibility or kinetic energy, which play the role of the effective carrier number. In contrast, the scattering time (or diffusion constant) saturates at high $T$. We find that $\ensuremath{\sigma}(n,T)\ensuremath{\approx}(1\ensuremath{-}n)\ensuremath{\sigma}(n=0,T)$ to a very good approximation for all $n$, with $\ensuremath{\sigma}(n=0,T)\ensuremath{\propto}t/T$ at high temperatures. The saturation at small $n$ occurs due to a compensation between the density dependence of the effective number of carriers and that of the scattering time. The $T$ dependence of the resistivity displays a kneelike feature which signals a crossover to the intermediate-temperature regime where the diffusion constant (or scattering time) starts increasing with decreasing $T$. At high temperatures, the thermopower obeys the Heikes formula, while the Wiedemann-Franz law is violated with the Lorenz number vanishing as $1/{T}^{2}$. The relevance of our calculations to experiments probing high-temperature transport in materials with strong electronic correlations or ultracold atomic gases in an optical lattice is briefly discussed.

DOI
10.1103/physrevb.94.235115

Related papers

Conductivity in the Bose-Hubbard model at high temperatures

Ivana Vasić, J. Vučičević

Quantum simulations are quickly becoming an indispensable tool for studying particle transport in correlated lattice models. One of the central topics in the study of transport is bad-metal behavior, characterized by the direct-current resistivity linear in temperature. In the fermionic Hubbard model, optical conductivity has been studied extensively, and a recent optical lattice experiment demonstrated bad-metal behavior in qualitative agreement with theory. Far less is known about transport in the bosonic Hubbard model. We investigate the conductivity in the Bose-Hubbard model and focus on the regime of strong interactions and high temperatures. We use numerically exact calculations for small lattice sizes. At weak tunneling, we find multiple peaks in the optical conductivity that stem from the Hubbard bands present in the many-body spectrum. This feature slowly washes out as the tunneling rate gets stronger. At high temperature, we identify a regime of $T$-linear resistivity, as expected. When the interactions are very strong, the leading inverse-temperature coefficient in conductivity is proportional to the tunneling amplitude. As the tunneling becomes stronger, this dependence takes quadratic form. At very strong coupling and half filling, we identify a separate linear resistivity regime at lower temperature, corresponding to the hard-core boson regime. Additionally, we unexpectedly observe that at half filling, in a big part of the phase diagram, conductivity is an increasing function of the coupling constant before it saturates at the hard-core-boson result. We explain this feature based on the analysis of the many-body energy spectrum and the contributions to conductivity of individual eigenstates of the system.

arXiv:2404.16559•Physical review. B./Physical review. B•Aug 1, 2024•3 citations

The Wiedemann-Franz law in doped Mott insulators without quasiparticles

Wen O. Wang, Jixun K. Ding, Yoni Schattner +3 more

Many metallic quantum materials display anomalous transport phenomena that defy a Fermi liquid description. Here, we use numerical methods to calculate thermal and charge transport in the doped Hubbard model and observe a crossover separating high- and low-temperature behaviors. Distinct from the behavior at high temperatures, the Lorenz number [Formula: see text] becomes weakly doping dependent and less sensitive to parameters at low temperatures. At the lowest numerically accessible temperatures, [Formula: see text] roughly approaches the Wiedemann-Franz constant [Formula: see text], even in a doped Mott insulator that lacks well-defined quasiparticles. Decomposing the energy current operator indicates a compensation between kinetic and potential contributions, which may help to clarify the interpretation of transport experiments beyond Boltzmann theory in strongly correlated metals.

arXiv:2208.09144•Science•Nov 30, 2023•15 citations

Transport Properties of the Infinite-Dimensional Hubbard Model

Thomas Pruschke, Daniel Lee Cox, Mark Jarrell

Results for the optical conductivity and resistivity of the Hubbard model in infinite spatial dimensions are presented. At half-filling we observe a gradual crossover from a metal with a Drude peak at ω = 0 in the optical conductivity to an insulator as a function of U for temperatures above the antiferromagnetic phase transition. When doped, the "insulator" becomes a Fermi liquid with a corresponding temperature dependence of the optical conductivity and resistivity. We find a T 2 -coefficient in the low-temperature resistivity which suggests that the carriers in the system acquire a considerable mass-enhancement due to the strong local correlations. At high temperatures, a crossover into a semi-metallic regime takes place.

arXiv:cond-mat/9208002•Europhysics Letters (EPL)•Feb 10, 1993•60 citations

Charge transport in the Hubbard model at high temperatures: Triangular versus square lattice

Ana Vranić, J. Vučičević, J. Kokalj +4 more

High-temperature bad-metal transport has been recently studied both theoretically and in experiments as one of the key signatures of strong electronic correlations. Here we use the dynamical mean field theory and its cluster extensions, as well as the finite-temperature Lanczos method to explore the influence of lattice frustration on the thermodynamic and transport properties of the Hubbard model at high temperatures. We consider the triangular and the square lattices at half-filling and at 15% hole doping. We find that for $T\ensuremath{\gtrsim}1.5t$ the self-energy becomes practically local, while the finite-size effects become small at lattice size $4\ifmmode\times\else\texttimes\fi{}4$ for both lattice types and doping levels. The vertex corrections to optical conductivity, which are significant on the square lattice even at high temperatures, contribute less on the triangular lattice. We find approximately linear temperature dependence of dc resistivity in doped Mott insulator for both types of lattices.

arXiv:2006.01707•Physical review. B./Physical review. B•Sep 21, 2020•42 citations

Spectral properties of the Hubbard bands

Henk J. Eskes, Andrzej M. Oleś, Marcel B. J. Meinders +1 more

Using strong-coupling perturbation theory for the Hubbard model, explicit expressions are obtained for the integrated weights, energy positions, and widths of the upper and lower Hubbard bands separately, both for the optical and photoemission spectrum. In one dimension all expressions can be explicitly evaluated using the large-U Bethe-ansatz wave function. The k\ensuremath{\rightarrow}-dependent moments for the one-particle spectrum are compared with the closely related two-pole ansatz. The strong momentum dependence of the spectra demonstrates the importance of the k\ensuremath{\rightarrow}-dependent correlation functions in the spectral moments, which are often neglected. In order to estimate corrections due to the neglect of higher-order terms (in t/U), all results are compared with numerical-diagonalization data for intermediate-U values. The sum-rule results show a rapid spectral-weight transfer from the high- to the low-energy regime upon doping, similar to what is observed experimentally in one-particle and optical spectra of Cu (high-${\mathit{T}}_{\mathit{c}}$) and Ni compounds. Such fast weight redistributions as a function of the carrier density are a natural consequence of the strong correlations in these materials. The redistribution of intensity in the optical conductivity is connected with the dilution of the spin system by the added holes. For the one-particle spectrum the decrease of weight in the upper Hubbard band away from half filling can be understood by state counting at U=\ensuremath{\infty}, and the effect is strongly enhanced for small doping by both the first- and second-order terms in t/U. Altogether, the one-particle and optical spectra show the importance of the three-site hopping term. Therefore the t-J model does not represent well the spectral properties of the Hubbard model at large U. The local spin order is essential and determines both the k\ensuremath{\rightarrow} dependence of the one-particle spectrum and the intensities in the optical spectrum.

Physical review. B, Condensed matter•Dec 15, 1994•158 citations

Repulsive versus attractive Hubbard model: Transport properties and spin-lattice relaxation rate

Rok Žitko, Žiga Osolin, P. Jeglič

We contrast the transport properties (dc resistivity, Seebeck coefficient), optical conductivity, spectral functions, dynamical magnetic susceptibility, and the nuclear magnetic resonance $1/{T}_{1}$ spin-lattice relaxation rate of the repulsive and attractive infinite-dimensional Hubbard models in the paramagnetic phase for a generic band filling. The calculations are performed in a wide temperature interval using the dynamical mean-field theory with the numerical renormalization group as the impurity solver. The attractive case exhibits significantly more complex temperature dependencies which can be explained by the behavior of the half-filled Hubbard model in external magnetic field with constant magnetization, to which the attractive Hubbard model maps through the partial particle-hole transformation. The resistivity is nonmonotonous for the strongly attractive case: it peaks significantly above the Mott-Ioffe-Regel value at a temperature ${T}_{\mathrm{max}}$ where the quasiparticle band disappears. For both signs of $U$ we find particle-hole asymmetry in the self-energy at low energies, but with the opposite kind of excitations having longer lifetime. This leads to a strong suppression of the slope of the Seebeck coefficient in the attractive case rather than an enhancement as in the repulsive case. The spin-lattice relaxation rate in the strongly attractive case has a nonmonotonic temperature dependence, thereby revealing the pairing fluctuations.

arXiv:1412.4506•Physical Review B•Apr 8, 2015•15 citations
We derive analytical expressions for the spectral moments of the dynamical response functions of the Hubbard model using the high-temperature series expansion. We consider generic dimension ddd as well as the infinite-ddd limit, arbitrary electron density nnn, and both finite and infinite repulsion UUU. We use moment-reconstruction methods to obtain the one-electron spectral function, the self-energy, and the optical conductivity. They are all smooth functions at high temperature and, at large UUU, they are featureless with characteristic widths of the order of the lattice hopping parameter ttt. In the infinite-ddd limit, we compare the series expansion results with accurate numerical renormalization group and interaction expansion quantum Monte Carlo results. We find excellent agreement down to surprisingly low temperatures, throughout most of the bad-metal regime, which applies for T\ensuremath≳(1\ensuremath−n)DT\ensuremath{\gtrsim}(1\ensuremath{-}n)DT\ensuremath≳(1\ensuremath−n)D, the Brinkman-Rice scale. The resistivity increases linearly in TTT at high temperature without saturation. This results from the 1/T1/T1/T behavior of the compressibility or kinetic energy, which play the role of the effective carrier number. In contrast, the scattering time (or diffusion constant) saturates at high TTT. We find that \ensuremathσ(n,T)\ensuremath≈(1\ensuremath−n)\ensuremathσ(n=0,T)\ensuremath{\sigma}(n,T)\ensuremath{\approx}(1\ensuremath{-}n)\ensuremath{\sigma}(n=0,T)\ensuremathσ(n,T)\ensuremath≈(1\ensuremath−n)\ensuremathσ(n=0,T) to a very good approximation for all nnn, with \ensuremathσ(n=0,T)\ensuremath∝t/T\ensuremath{\sigma}(n=0,T)\ensuremath{\propto}t/T\ensuremathσ(n=0,T)\ensuremath∝t/T at high temperatures. The saturation at small nnn occurs due to a compensation between the density dependence of the effective number of carriers and that of the scattering time. The TTT dependence of the resistivity displays a kneelike feature which signals a crossover to the intermediate-temperature regime where the diffusion constant (or scattering time) starts increasing with decreasing TTT. At high temperatures, the thermopower obeys the Heikes formula, while the Wiedemann-Franz law is violated with the Lorenz number vanishing as 1/T21/{T}^{2}1/T2. The relevance of our calculations to experiments probing high-temperature transport in materials with strong electronic correlations or ultracold atomic gases in an optical lattice is briefly discussed.
Quantum simulations are quickly becoming an indispensable tool for studying particle transport in correlated lattice models. One of the central topics in the study of transport is bad-metal behavior, characterized by the direct-current resistivity linear in temperature. In the fermionic Hubbard model, optical conductivity has been studied extensively, and a recent optical lattice experiment demonstrated bad-metal behavior in qualitative agreement with theory. Far less is known about transport in the bosonic Hubbard model. We investigate the conductivity in the Bose-Hubbard model and focus on the regime of strong interactions and high temperatures. We use numerically exact calculations for small lattice sizes. At weak tunneling, we find multiple peaks in the optical conductivity that stem from the Hubbard bands present in the many-body spectrum. This feature slowly washes out as the tunneling rate gets stronger. At high temperature, we identify a regime of TTT-linear resistivity, as expected. When the interactions are very strong, the leading inverse-temperature coefficient in conductivity is proportional to the tunneling amplitude. As the tunneling becomes stronger, this dependence takes quadratic form. At very strong coupling and half filling, we identify a separate linear resistivity regime at lower temperature, corresponding to the hard-core boson regime. Additionally, we unexpectedly observe that at half filling, in a big part of the phase diagram, conductivity is an increasing function of the coupling constant before it saturates at the hard-core-boson result. We explain this feature based on the analysis of the many-body energy spectrum and the contributions to conductivity of individual eigenstates of the system.
Many metallic quantum materials display anomalous transport phenomena that defy a Fermi liquid description. Here, we use numerical methods to calculate thermal and charge transport in the doped Hubbard model and observe a crossover separating high- and low-temperature behaviors. Distinct from the behavior at high temperatures, the Lorenz number [Formula: see text] becomes weakly doping dependent and less sensitive to parameters at low temperatures. At the lowest numerically accessible temperatures, [Formula: see text] roughly approaches the Wiedemann-Franz constant [Formula: see text], even in a doped Mott insulator that lacks well-defined quasiparticles. Decomposing the energy current operator indicates a compensation between kinetic and potential contributions, which may help to clarify the interpretation of transport experiments beyond Boltzmann theory in strongly correlated metals.
Results for the optical conductivity and resistivity of the Hubbard model in infinite spatial dimensions are presented. At half-filling we observe a gradual crossover from a metal with a Drude peak at ω = 0 in the optical conductivity to an insulator as a function of U for temperatures above the antiferromagnetic phase transition. When doped, the "insulator" becomes a Fermi liquid with a corresponding temperature dependence of the optical conductivity and resistivity. We find a T 2 -coefficient in the low-temperature resistivity which suggests that the carriers in the system acquire a considerable mass-enhancement due to the strong local correlations. At high temperatures, a crossover into a semi-metallic regime takes place.
High-temperature bad-metal transport has been recently studied both theoretically and in experiments as one of the key signatures of strong electronic correlations. Here we use the dynamical mean field theory and its cluster extensions, as well as the finite-temperature Lanczos method to explore the influence of lattice frustration on the thermodynamic and transport properties of the Hubbard model at high temperatures. We consider the triangular and the square lattices at half-filling and at 15% hole doping. We find that for $T\ensuremath{\gtrsim}1.5t$ the self-energy becomes practically local, while the finite-size effects become small at lattice size $4\ifmmode\times\else\texttimes\fi{}4$ for both lattice types and doping levels. The vertex corrections to optical conductivity, which are significant on the square lattice even at high temperatures, contribute less on the triangular lattice. We find approximately linear temperature dependence of dc resistivity in doped Mott insulator for both types of lattices.
Using strong-coupling perturbation theory for the Hubbard model, explicit expressions are obtained for the integrated weights, energy positions, and widths of the upper and lower Hubbard bands separately, both for the optical and photoemission spectrum. In one dimension all expressions can be explicitly evaluated using the large-U Bethe-ansatz wave function. The k\ensuremath{\rightarrow}-dependent moments for the one-particle spectrum are compared with the closely related two-pole ansatz. The strong momentum dependence of the spectra demonstrates the importance of the k\ensuremath{\rightarrow}-dependent correlation functions in the spectral moments, which are often neglected. In order to estimate corrections due to the neglect of higher-order terms (in t/U), all results are compared with numerical-diagonalization data for intermediate-U values. The sum-rule results show a rapid spectral-weight transfer from the high- to the low-energy regime upon doping, similar to what is observed experimentally in one-particle and optical spectra of Cu (high-Tc{\mathit{T}}_{\mathit{c}}Tc​) and Ni compounds. Such fast weight redistributions as a function of the carrier density are a natural consequence of the strong correlations in these materials. The redistribution of intensity in the optical conductivity is connected with the dilution of the spin system by the added holes. For the one-particle spectrum the decrease of weight in the upper Hubbard band away from half filling can be understood by state counting at U=\ensuremath{\infty}, and the effect is strongly enhanced for small doping by both the first- and second-order terms in t/U. Altogether, the one-particle and optical spectra show the importance of the three-site hopping term. Therefore the t-J model does not represent well the spectral properties of the Hubbard model at large U. The local spin order is essential and determines both the k\ensuremath{\rightarrow} dependence of the one-particle spectrum and the intensities in the optical spectrum.
We contrast the transport properties (dc resistivity, Seebeck coefficient), optical conductivity, spectral functions, dynamical magnetic susceptibility, and the nuclear magnetic resonance 1/T11/{T}_{1}1/T1​ spin-lattice relaxation rate of the repulsive and attractive infinite-dimensional Hubbard models in the paramagnetic phase for a generic band filling. The calculations are performed in a wide temperature interval using the dynamical mean-field theory with the numerical renormalization group as the impurity solver. The attractive case exhibits significantly more complex temperature dependencies which can be explained by the behavior of the half-filled Hubbard model in external magnetic field with constant magnetization, to which the attractive Hubbard model maps through the partial particle-hole transformation. The resistivity is nonmonotonous for the strongly attractive case: it peaks significantly above the Mott-Ioffe-Regel value at a temperature Tmax{T}_{\mathrm{max}}Tmax​ where the quasiparticle band disappears. For both signs of UUU we find particle-hole asymmetry in the self-energy at low energies, but with the opposite kind of excitations having longer lifetime. This leads to a strong suppression of the slope of the Seebeck coefficient in the attractive case rather than an enhancement as in the repulsive case. The spin-lattice relaxation rate in the strongly attractive case has a nonmonotonic temperature dependence, thereby revealing the pairing fluctuations.