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

Three-dimensional Hubbard model in the thermodynamic limit

Ehsan Khatami

Published Sep 9, 2016• 26 citations
arXivHTMLPDF

Abstract

We employ the numerical linked-cluster expansion to study finite-temperature properties of the uniform cubic lattice Hubbard model in the thermodynamic limit for a wide range of interaction strengths and densities. We carry out the expansion to the 9th order and find that the convergence of the series extends to lower temperatures as the strength of the interaction increases, giving us access to regions of the parameter space that are difficult to reach by most other numerical methods. We study the precise trends in the specific heat, the double occupancy, and magnetic correlations at temperatures as low as 0.2 of the hopping amplitude in the strong-coupling regime. We show that in this regime, accurate estimates for transition temperatures to the N\'eel ordered phase, in agreement with the predicted asymptotic behavior, can be deduced from the low-temperature magnetic structure factor. We also find evidence for possible instability to the magnetically ordered phase away from, but close to, half filling. Our results have important implications for parametrizing fermionic systems in optical lattice experiments and for benchmarking other numerical methods.

DOI
10.1103/physrevb.94.125114

Related papers

Thermodynamics of the disordered Hubbard model studied via numerical linked-cluster expansions

Jacob Park, Ehsan Khatami

The interplay of disorder and strong correlations in quantum many-body systems remains an open question. That is despite much progress made in recent years with ultracold atoms in optical lattices to better understand phenomena such as many-body localization or the effect of disorder on Mott metal-insulator transitions. Here, we utilize the numerical linked-cluster expansion technique, extended to treat disordered quantum lattice models, and study exact thermodynamic properties of the disordered Fermi-Hubbard model on the square and cubic geometries. We consider box distributions for the disorder in the onsite energy, the interaction strength, as well as the hopping amplitude and explore how energy, double occupancy, entropy, heat capacity, and magnetic correlations of the system in the thermodynamic limit evolve as the strength of disorder changes. We compare our findings with those obtained from determinant quantum Monte Carlo simulations and discuss the relevance of our results to experiments with cold fermionic atoms in optical lattices.

arXiv:2101.12721•Physical review. B./Physical review. B•Oct 4, 2021•6 citations

Finite-temperature superconducting correlations of the Hubbard model

Ehsan Khatami, Richard T. Scalettar, Rajiv R. P. Singh

We utilize numerical linked-cluster expansions (NLCEs) and the determinantal quantum Monte Carlo algorithm to study pairing correlations in the square-lattice Hubbard model. To benchmark the NLCE, we first locate the finite-temperature phase transition of the attractive model to a superconducting state away from half filling. We then explore the superconducting properties of the repulsive model for the $d$-wave and extended $s$-wave pairing symmetries. The pairing structure factor shows a strong tendency to $d$-wave pairing and peaks at an interaction strength comparable to the bandwidth. The extended $s$-wave structure factor and correlation length are larger at higher temperatures but clearly saturate as temperature is lowered, whereas the $d$-wave counterparts, which start off lower at high temperatures, continue to rise near half filling. This rise is even more dramatic in the $d$-wave susceptibility. The convergence of NLCEs breaks down as the susceptibilities and correlation lengths become large, so we are unable to determine the onset of long-range order. However, our results extend the conclusion, previously restricted to only magnetic and charge correlations, that NLCEs offer a unique window into pairing in the Hubbard model at strong coupling.

arXiv:1503.06213•Physical Review B•Jun 18, 2015•24 citations

Thermodynamics of strongly interacting fermions in two-dimensional optical lattices

Ehsan Khatami, Marcos Rigol

We study finite-temperature properties of strongly correlated fermions in two-dimensional optical lattices by means of numerical linked cluster expansions, a computational technique that allows one to obtain exact results in the thermodynamic limit. We focus our analysis on the strongly interacting regime, where the on-site repulsion is of the order of or greater than the band width. We compute the equation of state, double occupancy, entropy, uniform susceptibility, and spin correlations for temperatures that are similar to or below the ones achieved in current optical lattice experiments. We provide a quantitative analysis of adiabatic cooling of trapped fermions in two dimensions, by means of both flattening the trapping potential and increasing the interaction strength.

arXiv:1104.5494•Physical Review A•Nov 14, 2011•91 citations

Finite-Temperature Ferromagnetic Correlations of the Kagome Lattice Hubbard Model

Francisco Correia, Kyle Corbett, Ehsan Khatami

The Kagome lattice Fermi-Hubbard model is one of the most physically rich, and at the same time most challenging, models to study in strongly-correlated physics. Among its special features are geometric frustration and a flat energy band that create conditions favorable to ferromagnetism near the band insulating limit. Here, we utilize two exact finite-temperature methods, the numerical linked-cluster expansion and the determinant quantum Monte Carlo, to study the extent as well as doping and interaction dependence of ferromagnetic correlations and other thermodynamic properties of the model. We find that repulsive interactions enhance ferromagnetic correlations at high electron densities and that increasing the interaction strength, extends the region with strong ferromagnetic correlations towards half filling, smoothly connecting it to Nagaoka ferromagnetism near the Mott insulating region. We also use the charge compressibility to obtain an accurate estimate for the critical interaction strength for the metal-insulator transition at half filling. These results improve our understanding of the magnetic tendencies of the model away from half filling and pave the way for further studies, including with ultracold atoms in optical lattices.

arXiv:2604.26827•Apr 29, 2026•0 citations

Thermodynamic properties of one-dimensional hubbard model at finite temperatures

Yang C; Kocharian AN; Chiang YL

[[abstract]]The phase diagram and thermodynamic properties of the Hubbard model in one dimension are calculated numerically within the generalized self-consistent approach over a wide range of temperature T, magnetic field h, interaction strength U/t and electron concentration n. The temperature variation of the transverse magnetization, double occupied sites and chemical potential at various n, U/t and h is also analyzed. The wave vector q for the magnetic order parameter in the vicinity of the half-filling in the metallic phase shows a crossover from incommensurate magnetic phase (0 < q < pi) into the state with strong antiferromagnetic correlations (q = pi) as temperature is increased. Overall our numerical results in the approximate theory at finite temperatures are in good quantitative and qualitative agreement with the Bethe-ansatz results.

•Jan 1, 2003•0 citations

Thermodynamic Properties of One-Dimensional Hubbard Model at Finite Temperatures

Chi Yang, ARMEN N. KOCHARIAN, You-Ling Chiang

The phase diagram and thermodynamic properties of the Hubbard model in one dimension are calculated numerically within the generalized self-consistent approach over a wide range of temperature T, magnetic field h, interaction strength U/t and electron concentration n. The temperature variation of the transverse magnetization, double occupied sites and chemical potential at various n, U/t and h is also analyzed. The wave vector q for the magnetic order parameter in the vicinity of the half-filling in the metallic phase shows a crossover from incommensurate magnetic phase (0<q<π) into the state with strong antiferromagnetic correlations (q=π) as temperature is increased. Overall our numerical results in the approximate theory at finite temperatures are in good quantitative and qualitative agreement with the Bethe-ansatz results.

International Journal of Modern Physics B•Aug 10, 2003•0 citations
We employ the numerical linked-cluster expansion to study finite-temperature properties of the uniform cubic lattice Hubbard model in the thermodynamic limit for a wide range of interaction strengths and densities. We carry out the expansion to the 9th order and find that the convergence of the series extends to lower temperatures as the strength of the interaction increases, giving us access to regions of the parameter space that are difficult to reach by most other numerical methods. We study the precise trends in the specific heat, the double occupancy, and magnetic correlations at temperatures as low as 0.2 of the hopping amplitude in the strong-coupling regime. We show that in this regime, accurate estimates for transition temperatures to the N'eel ordered phase, in agreement with the predicted asymptotic behavior, can be deduced from the low-temperature magnetic structure factor. We also find evidence for possible instability to the magnetically ordered phase away from, but close to, half filling. Our results have important implications for parametrizing fermionic systems in optical lattice experiments and for benchmarking other numerical methods.
The interplay of disorder and strong correlations in quantum many-body systems remains an open question. That is despite much progress made in recent years with ultracold atoms in optical lattices to better understand phenomena such as many-body localization or the effect of disorder on Mott metal-insulator transitions. Here, we utilize the numerical linked-cluster expansion technique, extended to treat disordered quantum lattice models, and study exact thermodynamic properties of the disordered Fermi-Hubbard model on the square and cubic geometries. We consider box distributions for the disorder in the onsite energy, the interaction strength, as well as the hopping amplitude and explore how energy, double occupancy, entropy, heat capacity, and magnetic correlations of the system in the thermodynamic limit evolve as the strength of disorder changes. We compare our findings with those obtained from determinant quantum Monte Carlo simulations and discuss the relevance of our results to experiments with cold fermionic atoms in optical lattices.
We utilize numerical linked-cluster expansions (NLCEs) and the determinantal quantum Monte Carlo algorithm to study pairing correlations in the square-lattice Hubbard model. To benchmark the NLCE, we first locate the finite-temperature phase transition of the attractive model to a superconducting state away from half filling. We then explore the superconducting properties of the repulsive model for the ddd-wave and extended sss-wave pairing symmetries. The pairing structure factor shows a strong tendency to ddd-wave pairing and peaks at an interaction strength comparable to the bandwidth. The extended sss-wave structure factor and correlation length are larger at higher temperatures but clearly saturate as temperature is lowered, whereas the ddd-wave counterparts, which start off lower at high temperatures, continue to rise near half filling. This rise is even more dramatic in the ddd-wave susceptibility. The convergence of NLCEs breaks down as the susceptibilities and correlation lengths become large, so we are unable to determine the onset of long-range order. However, our results extend the conclusion, previously restricted to only magnetic and charge correlations, that NLCEs offer a unique window into pairing in the Hubbard model at strong coupling.
We study finite-temperature properties of strongly correlated fermions in two-dimensional optical lattices by means of numerical linked cluster expansions, a computational technique that allows one to obtain exact results in the thermodynamic limit. We focus our analysis on the strongly interacting regime, where the on-site repulsion is of the order of or greater than the band width. We compute the equation of state, double occupancy, entropy, uniform susceptibility, and spin correlations for temperatures that are similar to or below the ones achieved in current optical lattice experiments. We provide a quantitative analysis of adiabatic cooling of trapped fermions in two dimensions, by means of both flattening the trapping potential and increasing the interaction strength.
The Kagome lattice Fermi-Hubbard model is one of the most physically rich, and at the same time most challenging, models to study in strongly-correlated physics. Among its special features are geometric frustration and a flat energy band that create conditions favorable to ferromagnetism near the band insulating limit. Here, we utilize two exact finite-temperature methods, the numerical linked-cluster expansion and the determinant quantum Monte Carlo, to study the extent as well as doping and interaction dependence of ferromagnetic correlations and other thermodynamic properties of the model. We find that repulsive interactions enhance ferromagnetic correlations at high electron densities and that increasing the interaction strength, extends the region with strong ferromagnetic correlations towards half filling, smoothly connecting it to Nagaoka ferromagnetism near the Mott insulating region. We also use the charge compressibility to obtain an accurate estimate for the critical interaction strength for the metal-insulator transition at half filling. These results improve our understanding of the magnetic tendencies of the model away from half filling and pave the way for further studies, including with ultracold atoms in optical lattices.
[[abstract]]The phase diagram and thermodynamic properties of the Hubbard model in one dimension are calculated numerically within the generalized self-consistent approach over a wide range of temperature T, magnetic field h, interaction strength U/t and electron concentration n. The temperature variation of the transverse magnetization, double occupied sites and chemical potential at various n, U/t and h is also analyzed. The wave vector q for the magnetic order parameter in the vicinity of the half-filling in the metallic phase shows a crossover from incommensurate magnetic phase (0 < q < pi) into the state with strong antiferromagnetic correlations (q = pi) as temperature is increased. Overall our numerical results in the approximate theory at finite temperatures are in good quantitative and qualitative agreement with the Bethe-ansatz results.
The phase diagram and thermodynamic properties of the Hubbard model in one dimension are calculated numerically within the generalized self-consistent approach over a wide range of temperature T, magnetic field h, interaction strength U/t and electron concentration n. The temperature variation of the transverse magnetization, double occupied sites and chemical potential at various n, U/t and h is also analyzed. The wave vector q for the magnetic order parameter in the vicinity of the half-filling in the metallic phase shows a crossover from incommensurate magnetic phase (0<q<π) into the state with strong antiferromagnetic correlations (q=π) as temperature is increased. Overall our numerical results in the approximate theory at finite temperatures are in good quantitative and qualitative agreement with the Bethe-ansatz results.