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

Nonlinear Drude weight of the one-dimensional Hubbard model

Tetsuya Iwasaki, Hosho Katsura

Published Aug 20, 2026• 0 citations
arXivHTMLPDF

Abstract

We investigate nonlinear Drude weights (NLDWs) in the one-dimensional repulsive Hubbard model at zero temperature by combining exact Bethe-ansatz calculations with low-energy effective field theory. At quarter filling, we first derive the strong-coupling expansion of the NLDWs and confirm it numerically over a wide range of interaction strength. We then compare the numerical results with the prediction of the Tomonaga-Luttinger liquid (TLL) description including irrelevant perturbations. While band-curvature corrections yield finite contributions to higher-order NLDWs, the Umklapp interaction predicts divergent NLDWs when the order $n$ of the Drude weight exceeds a threshold determined by the TLL parameter. In contrast, finite-size scaling of the exact Bethe-ansatz results indicates that all calculated NLDWs remain finite in the thermodynamic limit, revealing a discrepancy between the exact results and the predictions of the low-energy effective field theory. At half filling, we analyze the finite-size scaling of the NLDWs across the Mott metal-insulator transition. We derive their asymptotic behavior in the insulating phase and propose a hyperscaling ansatz for NLDWs near the critical point, which is verified numerically. Our results clarify the interaction dependence and critical scaling of nonlinear transport coefficients in the one-dimensional Hubbard model and highlight limitations of the conventional low-energy effective description for higher-order transport.

DOI
10.48550/arxiv.2608.20269

Related papers

Charge and spin Drude weight of the one-dimensional extended Hubbard model at quarter filling

Tomonori Shirakawa, Eric Jeckelmann

We calculate the charge and spin Drude weight of the one-dimensional extended Hubbard model with on-site repulsion $U$ and nearest-neighbor repulsion $V$ at quarter filling using the density-matrix renormalization-group method combined with a variational principle. Our numerical results for the Hubbard model $(V=0)$ agree with exact results obtained from the Bethe ansatz solution. We obtain the contour map for both Drude weights in the $UV$-parameter space for repulsive interactions. We find that the charge Drude weight is discontinuous across the Kosterlitz-Thouless transition between the Luttinger liquid and the charge-density-wave insulator, while the spin Drude weight varies smoothly and remains finite in both phases. Our results can be generally understood using bosonization and renormalization-group results. The finite-size scaling of the charge Drude weight is well fitted by a polynomial function of the inverse system size in the metallic region. In the insulating region, we find an exponential decay of the finite-size corrections with the system size and a universal relation between the charge gap ${\ensuremath{\Delta}}_{c}$ and the correlation length $\ensuremath{\xi}$ which controls this exponential decay.

arXiv:0902.4139•Physical Review B•May 22, 2009•17 citations

Exact Drude weight for the one-dimensional Hubbard model at finite temperatures

Satoshi Fujimoto, Norio Kawakami

The Drude weight for the one-dimensional Hubbard model is investigated at finite temperatures by using the Bethe ansatz solution. Evaluating finite-size corrections to the thermodynamic Bethe ansatz equations, we obtain the formula for the Drude weight as the response of the system to an external gauge potential. We perform low-temperature expansions of the Drude weight in the case of half-filling as well as away from half-filling, which clearly distinguish the Mott-insulating state from the metallic state.

arXiv:cond-mat/9710313•Journal of Physics A Mathematical and General•Jan 16, 1998•85 citations

Fine structure of the nonlinear Drude weights in the spin- 12 XXZ chain

Yuhi Tanikawa, Hosho Katsura

We study nonlinear Drude weights (NLDWs) for the spin-$\frac{1}{2}$ XXZ chain in the critical regime at zero temperature. The NLDWs are generalizations of the linear Drude weight. Via the nonlinear extension of the Kohn formula, they can be read off from higher-order finite-size corrections to the ground-state energy in the presence of a $U(1)$ magnetic flux. The analysis of the ground-state energy based on the Bethe ansatz reveals that the NLDWs exhibit convergence, power-law, and logarithmic divergence, depending on the anisotropy parameter $\mathrm{\ensuremath{\Delta}}$. We determine the convergent and power-law divergent regions, which depend on the order of the response $n$. Then, we examine the behavior of the NLDWs at the boundary between the two regions and find that they converge for $n=0,1,2\phantom{\rule{4pt}{0ex}}(\mathrm{mod}\phantom{\rule{3.33333pt}{0ex}}4)$, while they show logarithmic divergence for $n=3\phantom{\rule{4pt}{0ex}}(\mathrm{mod}\phantom{\rule{3.33333pt}{0ex}}4)$. Furthermore, we identify particular anisotropies $\mathrm{\ensuremath{\Delta}}=cos(\ensuremath{\pi}r/(r+1))$ ($r=1,2,3,...$) at which the NLDW at any order $n$ converges to a finite value.

arXiv:2107.13784•Physical review. B./Physical review. B•Nov 15, 2021•9 citations

Transport and Temperature: Exact Spectrum and Drude Weight for the One-Dimensional Infinite- U Hubbard Model

Shuo Liu, Yuhao Ma, Hitesh J. Changlani +2 more

Understanding charge transport in strongly correlated systems remains a central challenge in condensed matter physics, while exact analytical results for finite-temperature transport remain rare. Here, we investigate charge transport in the one-dimensional Hubbard model in the infinite-interaction limit. We first construct the exact and explicit, i.e., beyond Bethe ansatz energy spectrum and eigenstates. Building on these explicit eigenstates, we derive a closed-form analytical expression for the charge Drude weight at arbitrary temperatures in the hole-dilute limit with a fixed number of doped holes, as well as analytical expressions in the low- and high-temperature limits at fixed hole density. We further show that, at low temperatures, the leading correction to the charge Drude weight is linear in T when the number of holes is fixed, whereas it is quadratic in T when the hole density is fixed. This analytically controlled example highlights the importance of identifying the underlying thermodynamic trajectory when interpreting finite-size numerical studies of charge transport in more general models.

arXiv:2604.02426•Physical Review Letters•Sep 9, 2026•0 citations

Exact results for nonlinear Drude weights in the spin- 12 XXZ chain

Yuhi Tanikawa, Kazuaki Takasan, Hosho Katsura

While the linear Drude weight in quantum many-body systems has been studied extensively, its nonlinear counterparts -- called nonlinear Drude weights (NLDWs) -- remain largely unexplored. Here, the authors present the first exact results for the NLDWs of the spin-$\frac{1}{2}$ XXZ chain, a paradigmatic quantum many-body system. Their exact results based on the Bethe ansatz reveal that the NLDWs diverge with the system size in certain parameter regions. Also, the analytical expressions for the NLDWs are obtained when they converge.

arXiv:2103.05838•Physical review. B./Physical review. B•May 27, 2021•15 citations

Exact numerical and analytical results for correlated lattice electrons in one dimension

Satoshi Ejima

Since the advent of high-$T_c$ cuprate superconductors in 1986, strongly correlated electron systems have attracted much attention. Since the cuprates are essentially two-dimensional, low-dimensional systems have moved into the focus of condensed-matter theory. From a theoretical point of view, one-dimensional systems are of particular interest because there are exact numerical and analytical methods which permit detailed studies and deep insights into the many-body problem. In the first part of this thesis, using the numerically exact methods Exact Diagonalization and the Density-Matrix Renormalization Group (DMRG), we investigate the properties of the Tomonaga–Luttinger liquid which is the generic metallic state of matter in one dimension. In particular, we concentrate on the investigation of the so-called Tomonaga–Luttinger liquid parameter which determines the critical exponent $\alpha$ for the density of states near the Fermi energy. Experimental results for some quasi one-dimensional materials report $\alpha \gtrsim 1$, which would imply $K_\rho \lesssim 0.17$, a value which cannot be reconciled with the bare Hubbard model where $K_\rho^H \geq 0.5$, i.e., $\alpha^H \lesssim 1/8$. We develop new accurate numerical methods to obtain $K_\rho$ and investigate how to obtain such small values for $K_\rho$ for slightly doped charge-density-wave insulators. In the second part of this thesis, using the Thermodynamic Bethe Ansatz (TBA) as exact analytical method, we investigate the one-dimensional Hubbard model in the spindisordered regime, which is characterized by the temperature being much larger than the magnetic energy scale but small compared to the Mott–Hubbard gap. Our study is motivated by the controversy about the Mott–Hubbard insulator in infinite dimensions whose ground state is also spin-disordered. In this system the determination of the precise value of the critical interaction strength Uc where the Mott–Hubbard gap closes is still unsolved. Therefore, we provide an example of a Hubbard-type model with a disordered spin background which can be solved exactly. The thermodynamics of our model can be understood in terms of gapped charged excitations with an effective dispersion which are decoupled form the spin degrees of freedom; the latter contribute only entropically. An interpretation of this regime in terms of a putative interacting-electron system at zero temperature leads to a metal-insulator transition at a finite interaction strength above which the gap opens linearly. Our exact results indicate that the strong-coupling expansion of the ground-state energy cannot be used to locate $U_c$. However, the strong-coupling expansion of the gap permits a reliable extrapolation of the critical interaction strength.

•May 3, 2006•0 citations
We investigate nonlinear Drude weights (NLDWs) in the one-dimensional repulsive Hubbard model at zero temperature by combining exact Bethe-ansatz calculations with low-energy effective field theory. At quarter filling, we first derive the strong-coupling expansion of the NLDWs and confirm it numerically over a wide range of interaction strength. We then compare the numerical results with the prediction of the Tomonaga-Luttinger liquid (TLL) description including irrelevant perturbations. While band-curvature corrections yield finite contributions to higher-order NLDWs, the Umklapp interaction predicts divergent NLDWs when the order nnn of the Drude weight exceeds a threshold determined by the TLL parameter. In contrast, finite-size scaling of the exact Bethe-ansatz results indicates that all calculated NLDWs remain finite in the thermodynamic limit, revealing a discrepancy between the exact results and the predictions of the low-energy effective field theory. At half filling, we analyze the finite-size scaling of the NLDWs across the Mott metal-insulator transition. We derive their asymptotic behavior in the insulating phase and propose a hyperscaling ansatz for NLDWs near the critical point, which is verified numerically. Our results clarify the interaction dependence and critical scaling of nonlinear transport coefficients in the one-dimensional Hubbard model and highlight limitations of the conventional low-energy effective description for higher-order transport.
We calculate the charge and spin Drude weight of the one-dimensional extended Hubbard model with on-site repulsion UUU and nearest-neighbor repulsion VVV at quarter filling using the density-matrix renormalization-group method combined with a variational principle. Our numerical results for the Hubbard model (V=0)(V=0)(V=0) agree with exact results obtained from the Bethe ansatz solution. We obtain the contour map for both Drude weights in the UVUVUV-parameter space for repulsive interactions. We find that the charge Drude weight is discontinuous across the Kosterlitz-Thouless transition between the Luttinger liquid and the charge-density-wave insulator, while the spin Drude weight varies smoothly and remains finite in both phases. Our results can be generally understood using bosonization and renormalization-group results. The finite-size scaling of the charge Drude weight is well fitted by a polynomial function of the inverse system size in the metallic region. In the insulating region, we find an exponential decay of the finite-size corrections with the system size and a universal relation between the charge gap ${\ensuremath{\Delta}}_{c}$ and the correlation length $\ensuremath{\xi}$ which controls this exponential decay.
The Drude weight for the one-dimensional Hubbard model is investigated at finite temperatures by using the Bethe ansatz solution. Evaluating finite-size corrections to the thermodynamic Bethe ansatz equations, we obtain the formula for the Drude weight as the response of the system to an external gauge potential. We perform low-temperature expansions of the Drude weight in the case of half-filling as well as away from half-filling, which clearly distinguish the Mott-insulating state from the metallic state.
We study nonlinear Drude weights (NLDWs) for the spin-12\frac{1}{2}21​ XXZ chain in the critical regime at zero temperature. The NLDWs are generalizations of the linear Drude weight. Via the nonlinear extension of the Kohn formula, they can be read off from higher-order finite-size corrections to the ground-state energy in the presence of a U(1)U(1)U(1) magnetic flux. The analysis of the ground-state energy based on the Bethe ansatz reveals that the NLDWs exhibit convergence, power-law, and logarithmic divergence, depending on the anisotropy parameter $\mathrm{\ensuremath{\Delta}}$. We determine the convergent and power-law divergent regions, which depend on the order of the response nnn. Then, we examine the behavior of the NLDWs at the boundary between the two regions and find that they converge for n=0,1,2(mod4)n=0,1,2\phantom{\rule{4pt}{0ex}}(\mathrm{mod}\phantom{\rule{3.33333pt}{0ex}}4)n=0,1,2(mod4), while they show logarithmic divergence for n=3(mod4)n=3\phantom{\rule{4pt}{0ex}}(\mathrm{mod}\phantom{\rule{3.33333pt}{0ex}}4)n=3(mod4). Furthermore, we identify particular anisotropies $\mathrm{\ensuremath{\Delta}}=cos(\ensuremath{\pi}r/(r+1))$ (r=1,2,3,...r=1,2,3,...r=1,2,3,...) at which the NLDW at any order nnn converges to a finite value.
Understanding charge transport in strongly correlated systems remains a central challenge in condensed matter physics, while exact analytical results for finite-temperature transport remain rare. Here, we investigate charge transport in the one-dimensional Hubbard model in the infinite-interaction limit. We first construct the exact and explicit, i.e., beyond Bethe ansatz energy spectrum and eigenstates. Building on these explicit eigenstates, we derive a closed-form analytical expression for the charge Drude weight at arbitrary temperatures in the hole-dilute limit with a fixed number of doped holes, as well as analytical expressions in the low- and high-temperature limits at fixed hole density. We further show that, at low temperatures, the leading correction to the charge Drude weight is linear in T when the number of holes is fixed, whereas it is quadratic in T when the hole density is fixed. This analytically controlled example highlights the importance of identifying the underlying thermodynamic trajectory when interpreting finite-size numerical studies of charge transport in more general models.
While the linear Drude weight in quantum many-body systems has been studied extensively, its nonlinear counterparts -- called nonlinear Drude weights (NLDWs) -- remain largely unexplored. Here, the authors present the first exact results for the NLDWs of the spin-12\frac{1}{2}21​ XXZ chain, a paradigmatic quantum many-body system. Their exact results based on the Bethe ansatz reveal that the NLDWs diverge with the system size in certain parameter regions. Also, the analytical expressions for the NLDWs are obtained when they converge.
Since the advent of high-TcT_cTc​ cuprate superconductors in 1986, strongly correlated electron systems have attracted much attention. Since the cuprates are essentially two-dimensional, low-dimensional systems have moved into the focus of condensed-matter theory. From a theoretical point of view, one-dimensional systems are of particular interest because there are exact numerical and analytical methods which permit detailed studies and deep insights into the many-body problem. In the first part of this thesis, using the numerically exact methods Exact Diagonalization and the Density-Matrix Renormalization Group (DMRG), we investigate the properties of the Tomonaga–Luttinger liquid which is the generic metallic state of matter in one dimension. In particular, we concentrate on the investigation of the so-called Tomonaga–Luttinger liquid parameter which determines the critical exponent α\alphaα for the density of states near the Fermi energy. Experimental results for some quasi one-dimensional materials report α≳1\alpha \gtrsim 1α≳1, which would imply Kρ≲0.17K_\rho \lesssim 0.17Kρ​≲0.17, a value which cannot be reconciled with the bare Hubbard model where KρH≥0.5K_\rho^H \geq 0.5KρH​≥0.5, i.e., αH≲1/8\alpha^H \lesssim 1/8αH≲1/8. We develop new accurate numerical methods to obtain KρK_\rhoKρ​ and investigate how to obtain such small values for KρK_\rhoKρ​ for slightly doped charge-density-wave insulators. In the second part of this thesis, using the Thermodynamic Bethe Ansatz (TBA) as exact analytical method, we investigate the one-dimensional Hubbard model in the spindisordered regime, which is characterized by the temperature being much larger than the magnetic energy scale but small compared to the Mott–Hubbard gap. Our study is motivated by the controversy about the Mott–Hubbard insulator in infinite dimensions whose ground state is also spin-disordered. In this system the determination of the precise value of the critical interaction strength Uc where the Mott–Hubbard gap closes is still unsolved. Therefore, we provide an example of a Hubbard-type model with a disordered spin background which can be solved exactly. The thermodynamics of our model can be understood in terms of gapped charged excitations with an effective dispersion which are decoupled form the spin degrees of freedom; the latter contribute only entropically. An interpretation of this regime in terms of a putative interacting-electron system at zero temperature leads to a metal-insulator transition at a finite interaction strength above which the gap opens linearly. Our exact results indicate that the strong-coupling expansion of the ground-state energy cannot be used to locate UcU_cUc​. However, the strong-coupling expansion of the gap permits a reliable extrapolation of the critical interaction strength.