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

Hubbard--Heisenberg Thermodynamic Comparison at Half Filling in a Fixed Staggered Field

Tadahiro Miyao

Published May 28, 2026• 0 citations
arXivHTMLPDF

Abstract

We study the repulsive Hubbard model at half filling in the canonical ensemble and compare its thermodynamics with the spin-$\frac12$ antiferromagnetic Heisenberg model at superexchange coupling $J_0(U)=4t^2/U$. For fixed hopping $t\ne0$, bounded staggered fields $|h|\le h_0$, and inverse temperatures $βJ_0(U)\ge\ell_0>0$, we establish a volume-uniform pressure comparison with error $O(U^{-3}+|h|U^{-2})$. On the Heisenberg field scale $h=J_0(U)x$, the error improves to $O(U^{-3})$ for bounded $x$. Using the quantum Lee--Yang theorem of Suzuki and Fisher, we obtain a volume-uniform $O(U^{-1})$ comparison of staggered magnetizations on compact positive intervals of the scaled field $x$. By combining a recent fermionic Lee--Yang theorem with analytic interpolation, we improve this estimate to $O(U^{-2}\log U)$ in the thermodynamic limit. The latter analysis also establishes analyticity of the limiting Hubbard pressure at nonzero real fields. The improved comparison extends to vanishing positive fields satisfying $U^{-3}\ll h(U)\ll U^{-1}$, with error $O((1+\log(U^3h(U)))/(U^3h(U)))$. Consequently, sufficiently low-temperature Heisenberg Néel order yields a quantitative positive lower bound on the Hubbard staggered magnetization in this regime. Furthermore, we prove $O(U^{-2})$ bounds on the double-occupancy and charge-defect densities and the normalized charge-density-wave second moment. At zero staggered field, this decay rate is optimal for both densities, even in the thermodynamic limit. Our approach uses finite-order Schrieffer--Wolff transformations and does not establish spontaneous Néel order at fixed finite $U$.

DOI
10.48550/arxiv.2605.30663

Related papers

Néel order, ring exchange, and charge fluctuations in the half-filled Hubbard model

Jean-Yves P. Delannoy, Michel J. P. Gingras, P. C. W. Holdsworth +1 more

We investigate the ground state properties of the two-dimensional half-filled one band Hubbard model in the strong (large-$U$) to intermediate coupling limit (i.e., away from the strict Heisenberg limit) using an effective spin-only low-energy theory that includes nearest-neighbor exchange, ring exchange, and all other spin interactions to order $t{(t∕U)}^{3}$. We show that the operator for the staggered magnetization, transformed for use in the effective theory, differs from that for the order parameter of the spin model by a renormalization factor accounting for the increased charge fluctuations as $t∕U$ is increased from the $t∕U\ensuremath{\rightarrow}0$ Heisenberg limit. These charge fluctuations lead to an increase of the quantum fluctuations over and above those for an $S=1∕2$ antiferromagnet. The renormalization factor ensures that the zero temperature staggered moment for the Hubbard model is a monotonously decreasing function of $t∕U$, despite the fact that the moment of the spin Hamiltonian, which depends on transverse spin fluctuations only, in an increasing function of $t∕U$. We also comment on quantitative aspects of the $t∕U$ and $1∕S$ expansions.

arXiv:cond-mat/0412033•Physical Review B•Sep 16, 2005•53 citations

Temporal decay of Néel order in the one-dimensional Fermi-Hubbard model

A. Bauer, Florian Dorfner, Fabian Heidrich‐Meisner

Motivated by recent experiments with ultracold quantum gases in optical lattices we study the decay of the staggered moment in the one-dimensional Fermi-Hubbard model starting from a perfect N\'eel state using exact diagonalization and the infinite-system-size time-evolving-block-decimation method. This extends previous work in which the same problem has been addressed for pure spin Hamiltonians. As a main result, we show that the relaxation dynamics of the double occupancy and of the staggered moment are different. The former is controlled by the nearest-neighbor tunneling rate while the latter is much slower and strongly dependent on the interaction strength, indicating that spin excitations are important. This difference in characteristic energy scales for the fast charge dynamics and the much slower spin dynamics is also reflected in the real-time evolution of nearest-neighbor density and spin correlations. A very interesting time dependence emerges in the von Neumann entropy, which at short times increases linearly with a slope proportional to the tunneling matrix element while the long-time growth of entanglement is controlled by spin excitations. Our predictions for the different relaxation dynamics of the staggered moment and the double occupancy should be observable in state-of-the-art optical lattice experiments. We further compare time averages of the double occupancy to the expectation values in both the canonical and diagonal ensembles, which quantitatively disagree with each other on finite systems. We relate the question of thermalization to the eigenstate thermalization hypothesis.

arXiv:1503.02010•Physical Review A•May 27, 2015•21 citations

Phase diagram of the half-filled extended Hubbard model in two dimensions

Biplab Chattopadhyay, D. M. Gaitonde

We consider an extended Hubbard model of interacting fermions on a lattice. The fermion kinetic energy corresponds to a tight-binding Hamiltonian with nearest-neighbor (-t) and next-nearest-neighbor (${\mathrm{t}}^{|\mathrm{I}\mathrm{H}}$) hopping matrix elements. In addition to the on-site Hubbard interaction (U) we also consider a nearest-neighbor repulsion (V). We obtain the zero-temperature phase diagram of our model within the Hartree-Fock approximation. We consider ground states having charge- and spin-density-wave ordering as well as states with orbital antiferromagnetism or spin nematic order. The latter two states correspond to particle-hole binding with ${\mathrm{d}}_{{\mathrm{x}}^{2}\mathrm{\ensuremath{-}}{\mathrm{y}}^{2}}$ symmetry in the charge and spin channels, respectively. For ${\mathrm{t}}^{|\mathrm{I}\mathrm{H}}$=0, only the charge-density-wave and spin-density-wave states are energetically stable. For nonzero ${\mathrm{t}}^{|\mathrm{I}\mathrm{H}}$, we find that orbital antiferromagnetism (or spin nematic) order is stable over a finite portion of the phase diagram at weak coupling. This region of stability is seen to grow with increasing values of ${\mathrm{t}}^{|\mathrm{I}\mathrm{H}}$.

arXiv:cond-mat/9703213•Physical review. B, Condensed matter•Jun 15, 1997•24 citations

Onsager-reaction-field correction in the half-filled Hubbard model

Peter Kopietz

The Onsager-reaction-field correction to Hartree-Fock theory recently derived by Georges and Yedidia [Phys. Rev. B 43, 3475 (1991)] is evaluated analytically for the half-filled positive-U Hubbard model on a d-dimensional hypercubic lattice. It is shown that the Onsager correction is for all d\ensuremath{\ge}3 irrelevant in the weak-coupling limit, while for d=1 and d=2 it strongly renormalizes the Hartree-Fock result, although it does note completely destabilize the spin-density wave. We examine the perturbation theory for the energy change induced by a small staggered magnetization m to all orders in U, and show that one obtains an infinite series in powers of (U/t)ln[${\mathit{m}}^{\mathrm{\ensuremath{-}}1}$ln(${\mathit{m}}^{\mathrm{\ensuremath{-}}1}$)] for d=1, and in powers of (U/t)${\mathrm{ln}}^{2}$[${\mathit{m}}^{\mathrm{\ensuremath{-}}1}$${\mathrm{ln}}^{2}$(${\mathit{m}}^{\mathrm{\ensuremath{-}}1}$)] for d=1. The close similarity between d=1 and d=2 is interpreted as evidence that antiferromagnetism in the weak coupling regime is not stable in d=2.

Physical review. B, Condensed matter•Nov 1, 1993•15 citations

Three-dimensional repulsive Hubbard model

Alexei Sherman

The three-dimensional repulsive Hubbard model is investigated using the strong-coupling diagram technique. For half-filling, the boundary between paramagnetic and antiferromagnetic states is determined for the range of the Hubbard repulsion $4t\leq U\leq12t$, where $t$ is the hopping integral between neighboring sites. Along this boundary, the density of states is calculated, and it demonstrates the Mott transition at $U\approx9t$. For $U\geq6t$ and half-filling, the density of states has the shape inherent in the strong electron repulsion, while for $U=4t$, its shape points to weak coupling. The dependence of the Néel temperature $T_{\rm N}$ on the electron concentration $\bar{n}$ is investigated for the cases $U=4t$ and $12t$. In the former case, $T_{\rm N}$ decreases monotonously with $\bar{n}$, while in the latter case, there is a plateau in the dependence near $\bar{n}=0.87$. The plateau is connected to a reconstruction of the density of states caused by an effective weakening of electron coupling due to electron depopulation. For $U=12t$ and half-filling, the magnetic critical exponent $γ\approx1.4$, which is close to the value in the Heisenberg model. Some features resembling the first-order phase transition, revealing themselves in a finite crystal, are discussed.

arXiv:2606.17933•Jun 16, 2026•0 citations

Quantum critical behavior and thermodynamics of the repulsive one-dimensional Hubbard model in a magnetic field

Ovidiu I. Pâţu, Andreas Klümper, Angela Foerster

Even though the Hubbard model is one of the most fundamental models of highly correlated electrons, analytical and numerical data describing its thermodynamics at nonzero magnetization are relatively scarce. We present a detailed investigation of the thermodynamic properties for the one-dimensional repulsive Hubbard model in the presence of an arbitrary magnetic field for all values of the filling fraction and temperatures as low as $T\ensuremath{\sim}0.005t.$ Our analysis is based on the system of integral equations derived in the quantum transfer matrix framework. We determine the critical exponents of the quantum phase transitions and also provide analytical derivations for some of the universal functions characterizing the thermodynamics in the vicinities of the quantum critical points. Extensive numerical data for the specific heat, susceptibility, compressibility, and entropy are reported. The experimentally relevant double occupancy presents an interesting doubly nonmonotonic temperature dependence at intermediate values of the interaction strength and also at large repulsion and magnetic fields close to the critical value. The susceptibility in zero magnetic field has a logarithmic singularity at low temperatures for all filling factors similar to the behavior of the same quantity in the spin-1/2 isotropic Heisenberg model. We determine the density profiles for a harmonically trapped system and show that while the total density profile seems to depend mainly on the value of chemical potential at the center of the trap the distribution of phases in the inhomogeneous system changes dramatically as we increase the magnetic field.

arXiv:1910.09535•Physical review. B./Physical review. B•Jan 29, 2020•10 citations
We study the repulsive Hubbard model at half filling in the canonical ensemble and compare its thermodynamics with the spin-12\frac1221​ antiferromagnetic Heisenberg model at superexchange coupling J0(U)=4t2/UJ_0(U)=4t^2/UJ0​(U)=4t2/U. For fixed hopping t≠0t\ne0t=0, bounded staggered fields ∣h∣≤h0|h|\le h_0∣h∣≤h0​, and inverse temperatures βJ0(U)≥ℓ0>0βJ_0(U)\ge\ell_0>0βJ0​(U)≥ℓ0​>0, we establish a volume-uniform pressure comparison with error O(U−3+∣h∣U−2)O(U^{-3}+|h|U^{-2})O(U−3+∣h∣U−2). On the Heisenberg field scale h=J0(U)xh=J_0(U)xh=J0​(U)x, the error improves to O(U−3)O(U^{-3})O(U−3) for bounded xxx. Using the quantum Lee--Yang theorem of Suzuki and Fisher, we obtain a volume-uniform O(U−1)O(U^{-1})O(U−1) comparison of staggered magnetizations on compact positive intervals of the scaled field xxx. By combining a recent fermionic Lee--Yang theorem with analytic interpolation, we improve this estimate to O(U−2log⁡U)O(U^{-2}\log U)O(U−2logU) in the thermodynamic limit. The latter analysis also establishes analyticity of the limiting Hubbard pressure at nonzero real fields. The improved comparison extends to vanishing positive fields satisfying U−3≪h(U)≪U−1U^{-3}\ll h(U)\ll U^{-1}U−3≪h(U)≪U−1, with error O((1+log⁡(U3h(U)))/(U3h(U)))O((1+\log(U^3h(U)))/(U^3h(U)))O((1+log(U3h(U)))/(U3h(U))). Consequently, sufficiently low-temperature Heisenberg Néel order yields a quantitative positive lower bound on the Hubbard staggered magnetization in this regime. Furthermore, we prove O(U−2)O(U^{-2})O(U−2) bounds on the double-occupancy and charge-defect densities and the normalized charge-density-wave second moment. At zero staggered field, this decay rate is optimal for both densities, even in the thermodynamic limit. Our approach uses finite-order Schrieffer--Wolff transformations and does not establish spontaneous Néel order at fixed finite UUU.
We investigate the ground state properties of the two-dimensional half-filled one band Hubbard model in the strong (large-UUU) to intermediate coupling limit (i.e., away from the strict Heisenberg limit) using an effective spin-only low-energy theory that includes nearest-neighbor exchange, ring exchange, and all other spin interactions to order t(t∕U)3t{(t∕U)}^{3}t(t∕U)3. We show that the operator for the staggered magnetization, transformed for use in the effective theory, differs from that for the order parameter of the spin model by a renormalization factor accounting for the increased charge fluctuations as t∕Ut∕Ut∕U is increased from the $t∕U\ensuremath{\rightarrow}0$ Heisenberg limit. These charge fluctuations lead to an increase of the quantum fluctuations over and above those for an S=1∕2S=1∕2S=1∕2 antiferromagnet. The renormalization factor ensures that the zero temperature staggered moment for the Hubbard model is a monotonously decreasing function of t∕Ut∕Ut∕U, despite the fact that the moment of the spin Hamiltonian, which depends on transverse spin fluctuations only, in an increasing function of t∕Ut∕Ut∕U. We also comment on quantitative aspects of the t∕Ut∕Ut∕U and 1∕S1∕S1∕S expansions.
Motivated by recent experiments with ultracold quantum gases in optical lattices we study the decay of the staggered moment in the one-dimensional Fermi-Hubbard model starting from a perfect N\'eel state using exact diagonalization and the infinite-system-size time-evolving-block-decimation method. This extends previous work in which the same problem has been addressed for pure spin Hamiltonians. As a main result, we show that the relaxation dynamics of the double occupancy and of the staggered moment are different. The former is controlled by the nearest-neighbor tunneling rate while the latter is much slower and strongly dependent on the interaction strength, indicating that spin excitations are important. This difference in characteristic energy scales for the fast charge dynamics and the much slower spin dynamics is also reflected in the real-time evolution of nearest-neighbor density and spin correlations. A very interesting time dependence emerges in the von Neumann entropy, which at short times increases linearly with a slope proportional to the tunneling matrix element while the long-time growth of entanglement is controlled by spin excitations. Our predictions for the different relaxation dynamics of the staggered moment and the double occupancy should be observable in state-of-the-art optical lattice experiments. We further compare time averages of the double occupancy to the expectation values in both the canonical and diagonal ensembles, which quantitatively disagree with each other on finite systems. We relate the question of thermalization to the eigenstate thermalization hypothesis.
We consider an extended Hubbard model of interacting fermions on a lattice. The fermion kinetic energy corresponds to a tight-binding Hamiltonian with nearest-neighbor (-t) and next-nearest-neighbor (t∣IH{\mathrm{t}}^{|\mathrm{I}\mathrm{H}}t∣IH) hopping matrix elements. In addition to the on-site Hubbard interaction (U) we also consider a nearest-neighbor repulsion (V). We obtain the zero-temperature phase diagram of our model within the Hartree-Fock approximation. We consider ground states having charge- and spin-density-wave ordering as well as states with orbital antiferromagnetism or spin nematic order. The latter two states correspond to particle-hole binding with ${\mathrm{d}}_{{\mathrm{x}}^{2}\mathrm{\ensuremath{-}}{\mathrm{y}}^{2}}$ symmetry in the charge and spin channels, respectively. For t∣IH{\mathrm{t}}^{|\mathrm{I}\mathrm{H}}t∣IH=0, only the charge-density-wave and spin-density-wave states are energetically stable. For nonzero t∣IH{\mathrm{t}}^{|\mathrm{I}\mathrm{H}}t∣IH, we find that orbital antiferromagnetism (or spin nematic) order is stable over a finite portion of the phase diagram at weak coupling. This region of stability is seen to grow with increasing values of t∣IH{\mathrm{t}}^{|\mathrm{I}\mathrm{H}}t∣IH.
The Onsager-reaction-field correction to Hartree-Fock theory recently derived by Georges and Yedidia [Phys. Rev. B 43, 3475 (1991)] is evaluated analytically for the half-filled positive-U Hubbard model on a d-dimensional hypercubic lattice. It is shown that the Onsager correction is for all d\ensuremath{\ge}3 irrelevant in the weak-coupling limit, while for d=1 and d=2 it strongly renormalizes the Hartree-Fock result, although it does note completely destabilize the spin-density wave. We examine the perturbation theory for the energy change induced by a small staggered magnetization m to all orders in U, and show that one obtains an infinite series in powers of (U/t)ln[${\mathit{m}}^{\mathrm{\ensuremath{-}}1}$ln(${\mathit{m}}^{\mathrm{\ensuremath{-}}1}$)] for d=1, and in powers of (U/t)ln2{\mathrm{ln}}^{2}ln2[${\mathit{m}}^{\mathrm{\ensuremath{-}}1}$ln2{\mathrm{ln}}^{2}ln2(${\mathit{m}}^{\mathrm{\ensuremath{-}}1}$)] for d=1. The close similarity between d=1 and d=2 is interpreted as evidence that antiferromagnetism in the weak coupling regime is not stable in d=2.
The three-dimensional repulsive Hubbard model is investigated using the strong-coupling diagram technique. For half-filling, the boundary between paramagnetic and antiferromagnetic states is determined for the range of the Hubbard repulsion 4t≤U≤12t4t\leq U\leq12t4t≤U≤12t, where ttt is the hopping integral between neighboring sites. Along this boundary, the density of states is calculated, and it demonstrates the Mott transition at U≈9tU\approx9tU≈9t. For U≥6tU\geq6tU≥6t and half-filling, the density of states has the shape inherent in the strong electron repulsion, while for U=4tU=4tU=4t, its shape points to weak coupling. The dependence of the Néel temperature TNT_{\rm N}TN​ on the electron concentration nˉ\bar{n}nˉ is investigated for the cases U=4tU=4tU=4t and 12t12t12t. In the former case, TNT_{\rm N}TN​ decreases monotonously with nˉ\bar{n}nˉ, while in the latter case, there is a plateau in the dependence near nˉ=0.87\bar{n}=0.87nˉ=0.87. The plateau is connected to a reconstruction of the density of states caused by an effective weakening of electron coupling due to electron depopulation. For U=12tU=12tU=12t and half-filling, the magnetic critical exponent γ≈1.4γ\approx1.4γ≈1.4, which is close to the value in the Heisenberg model. Some features resembling the first-order phase transition, revealing themselves in a finite crystal, are discussed.
Even though the Hubbard model is one of the most fundamental models of highly correlated electrons, analytical and numerical data describing its thermodynamics at nonzero magnetization are relatively scarce. We present a detailed investigation of the thermodynamic properties for the one-dimensional repulsive Hubbard model in the presence of an arbitrary magnetic field for all values of the filling fraction and temperatures as low as $T\ensuremath{\sim}0.005t.$ Our analysis is based on the system of integral equations derived in the quantum transfer matrix framework. We determine the critical exponents of the quantum phase transitions and also provide analytical derivations for some of the universal functions characterizing the thermodynamics in the vicinities of the quantum critical points. Extensive numerical data for the specific heat, susceptibility, compressibility, and entropy are reported. The experimentally relevant double occupancy presents an interesting doubly nonmonotonic temperature dependence at intermediate values of the interaction strength and also at large repulsion and magnetic fields close to the critical value. The susceptibility in zero magnetic field has a logarithmic singularity at low temperatures for all filling factors similar to the behavior of the same quantity in the spin-1/2 isotropic Heisenberg model. We determine the density profiles for a harmonically trapped system and show that while the total density profile seems to depend mainly on the value of chemical potential at the center of the trap the distribution of phases in the inhomogeneous system changes dramatically as we increase the magnetic field.