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:2305.05683 · Physical Review A

Nagaoka ferromagnetism in doped Hubbard models in optical lattices

Rhine Samajdar, R. N. Bhatt

Published Aug 22, 2024• 14 citations
arXivHTMLPDF

Abstract

Finding ferromagnetic ground states of the Hubbard model, which usually features antiferromagnetic spin correlations, has been an outstanding challenge since Nagaoka's celebrated (but experimentally unrealistic) proposal in 1966. Here, the authors show how such elusive itinerant ferromagnetism can be enhanced by kinetic frustration or occupation-dependent hoppings and demonstrate the relevance of these mechanisms to recent optical-lattice experiments.

DOI
10.1103/physreva.110.l021303

Related papers

Evolution of Nagaoka phase with kinetic energy frustrating hopping

Franco T. Lisandrini, Barbara Bravo, A. E. Trumper +2 more

We investigate, using the density-matrix renormalization group, the evolution of the Nagaoka state with ${t}^{\ensuremath{'}}$ hopping that frustrates the hole kinetic energy in the $U=\ensuremath{\infty}$ Hubbard model on the square and anisotropic triangular lattices. We find that the Nagaoka ferromagnet survives up to a rather small ${t}_{c}^{\ensuremath{'}}/t\ensuremath{\sim}0.2.$ At this critical value, there is a transition to an antiferromagnetic phase that depends on the lattice: a $\mathbf{Q}=(Q,0)$ spiral order, which continuously evolves with ${t}^{\ensuremath{'}}$, for the triangular lattice and the usual $\mathbf{Q}=(\ensuremath{\pi},\ensuremath{\pi})$ N\'eel order for the square lattice. Remarkably, the local magnetization takes its classical value for all considered ${t}^{\ensuremath{'}}$ (${t}^{\ensuremath{'}}/t\ensuremath{\le}1$). Our results show that the recently found classical kinetic antiferromagnetism, a perfect counterpart of Nagaoka ferromagnetism, is a generic phenomenon in these kinetically frustrated electronic systems.

arXiv:1702.00417•Physical review. B./Physical review. B•May 2, 2017•14 citations

Quantum phases in Fermi Hubbard systems with tunable frustration

Xu, Muqing

Analog quantum simulation provides a unique toolkit to investigate quantum many-body problems of whose solutions may be incredibly challenging. In particular, using ultracold fermionic atoms in optical lattices it realizes the Fermi Hubbard model, which is a fundamental model in condensed matter physics exhibiting properties relevant to many intriguing strongly-correlated systems. However, most quantum simulators of the Hubbard model have so far restricted to the square geometries. In addition, access to low temperatures where the exotic quantum phases are predicted to develop has still remained elusive. In this thesis, we engineer a novel optical lattice with dynamically tunable geometries and significantly reduced technical noises. We extend for the first time the study of quantum magnetism into the exotic phases in a Fermi Hubbard system with tunable geometric frustration. As we continuously tune our system from a square to a triangular lattice geometry, we observe a transition from a Néel antiferromagnet to a 120-degree spin spiral state. We then study, again for the first time, how doping affects a frustrated quantum magnet. To our surprise, we found an emergent ferromagnetic state with strong particle doping, which is absent on the hole-doped side and is in stark contrast with the particle-hole symmetric doping dependence in a square lattice. To shed light on the role of itinerant dopants in the emergence of these new magnetic properties, we furthermore leverage the unique single-particle resolution capabilities of our platform to probe higher-order density and spin correlations. These measurements hint at a new type of magnetism induced by kinetic frustration, which manifests as local antiferromagnetic correlations around hole dopants and and ferromagnetic correlations around particle dopants. Unlike magnetism induced by exchange or superexchange interactions, the energy associated with this mechanism is the kinetic energy, closely related to the famous rigorous results by Nagaoka. Furthermore, the dynamical tunability of the lattice allows adibatic engineering of quantum states in optical lattices. We discuss the preliminery efforts to prepare strongly-correlated states with significantly reduced temperatures, which may pave the way towards low temperature states in the Hubbard model.

•Jan 1, 2024•0 citations

Finite-temperature kinetic ferromagnetism in the square-lattice Hubbard model

Robert C. Newby, Ehsan Khatami

While the exact phase diagram of the Fermi-Hubbard model remains poorly understood despite decades of progress, nearly 60 years ago Nagaoka proved that a single dopant in an otherwise half-filled Hubbard system can bring about ferromagnetism through kinetic means. The phenomenon was recently observed with ultracold atoms in triangular optical lattices. Here, we explore the kinetic ferromagnetism within the square lattice Hubbard model and its strong-coupling counterpart, the $t\ensuremath{-}J$ model, at finite temperatures in the thermodynamic limit via numerical linked-cluster expansions. We find evidence of ferromagnetic Nagaoka polarons at dopings up to $\ensuremath{\sim}30%$ away from half filling for a variety of interaction strengths and at temperatures as low as 0.2 of the hopping energy. We map out the boundaries of this phase through analyzing various correlation functions.

arXiv:2502.07252•Physical review. B./Physical review. B•Jun 6, 2025•6 citations

Geometric frustration assisted kinetic ferromagnetism in doped Mott insulators

Qianqian Chen, Shuai A. Chen, Zheng Zhu

Understanding ferromagnetism mechanism in doped Mott insulators on frustrated lattices remains challenging at intermediate coupling and finite doping. Here, we study the itinerant ferromagnetism and propose its mechanism in doped Mott insulators on a geometrically frustrated triangular lattice. Using large-scale density matrix renormalization group (DMRG) and unrestricted Hartree–Fock mean-field methods, we reveal that itinerant ferromagnetism appears at intermediate coupling (10 ≲ U ≪ ∞) near 50% electron doping in the triangular-lattice Hubbard model. By analyzing all microscopic hopping processes, we find that doublon-singlon exchange alone drives the fully polarized ferromagnetism and uncovers the particle-hole asymmetry. We also establish the magnetic phase diagram and compare local spin correlations with recent experiments. Random phase approximation and DMRG calculations consistently confirm that the ferromagnetism persists when SU(2) symmetry is explicitly broken by magnetic anisotropy. These results clarify a microscopic route to itinerant ferromagnetism at intermediate coupling and finite doping in doped Mott insulators. Understanding the mechanism of ferromagnetism in strongly correlated systems is an ongoing theoretical challenge and the Hubbard model is typically adopted to investigate such systems. Here, the authors numerically investigate the doped triangular Fermi-Hubbard model and show ferromagnetism at intermediate coupling and finite doping, which results from the interplay of kinetic doublon-singlon exchange and lattice geometry.

arXiv:2408.05971•Communications Physics•Oct 10, 2025•1 citation

Nagaoka ferromagnetism in large-spin fermionic and bosonic systems

Seiji Miyashita, Masao Ogata, Hans A. De Raedt

We study the magnetic properties of itinerant quantum magnetic particles, described by a generalized Hubbard model with large spin $(S>1/2)$, which may be realized in optical lattices of laser-cooled atom systems. In fermion systems (half-integer spins), an extended form of Nagaoka ferromagnetism may be realized. However, as novel aspects of the large-spin cases, we found that the condition on the lattice connectivity is more stringent than in the case of $S=1/2$ particles and that the system shows a peculiar degenerate structure of the ground state in which the ferromagnetic state is included. In contrast, it turns out that the ground state of itinerant bosonic systems (integer spins) has a degenerate structure similar to that of fermion system with $S>1/2$ regardless of the shape, connectivity, or filling of the lattice, and that the state with the maximum total spin is always one of the ground states. Because the system consists of $2S+1$ types of particles and we study a $\text{SU}(2S+1)$ invariant model, the degeneracy of the ground state is given by the multiplets of the fully symmetric Young tableau of $\text{SU}(2S+1)$ if the state with maximum total spin belongs to the ground state.

arXiv:0909.2082•Physical Review B•Nov 24, 2009•6 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
Finding ferromagnetic ground states of the Hubbard model, which usually features antiferromagnetic spin correlations, has been an outstanding challenge since Nagaoka's celebrated (but experimentally unrealistic) proposal in 1966. Here, the authors show how such elusive itinerant ferromagnetism can be enhanced by kinetic frustration or occupation-dependent hoppings and demonstrate the relevance of these mechanisms to recent optical-lattice experiments.
We investigate, using the density-matrix renormalization group, the evolution of the Nagaoka state with ${t}^{\ensuremath{'}}$ hopping that frustrates the hole kinetic energy in the $U=\ensuremath{\infty}$ Hubbard model on the square and anisotropic triangular lattices. We find that the Nagaoka ferromagnet survives up to a rather small ${t}_{c}^{\ensuremath{'}}/t\ensuremath{\sim}0.2.$ At this critical value, there is a transition to an antiferromagnetic phase that depends on the lattice: a Q=(Q,0)\mathbf{Q}=(Q,0)Q=(Q,0) spiral order, which continuously evolves with ${t}^{\ensuremath{'}}$, for the triangular lattice and the usual $\mathbf{Q}=(\ensuremath{\pi},\ensuremath{\pi})$ N\'eel order for the square lattice. Remarkably, the local magnetization takes its classical value for all considered ${t}^{\ensuremath{'}}$ (${t}^{\ensuremath{'}}/t\ensuremath{\le}1$). Our results show that the recently found classical kinetic antiferromagnetism, a perfect counterpart of Nagaoka ferromagnetism, is a generic phenomenon in these kinetically frustrated electronic systems.
Analog quantum simulation provides a unique toolkit to investigate quantum many-body problems of whose solutions may be incredibly challenging. In particular, using ultracold fermionic atoms in optical lattices it realizes the Fermi Hubbard model, which is a fundamental model in condensed matter physics exhibiting properties relevant to many intriguing strongly-correlated systems. However, most quantum simulators of the Hubbard model have so far restricted to the square geometries. In addition, access to low temperatures where the exotic quantum phases are predicted to develop has still remained elusive. In this thesis, we engineer a novel optical lattice with dynamically tunable geometries and significantly reduced technical noises. We extend for the first time the study of quantum magnetism into the exotic phases in a Fermi Hubbard system with tunable geometric frustration. As we continuously tune our system from a square to a triangular lattice geometry, we observe a transition from a Néel antiferromagnet to a 120-degree spin spiral state. We then study, again for the first time, how doping affects a frustrated quantum magnet. To our surprise, we found an emergent ferromagnetic state with strong particle doping, which is absent on the hole-doped side and is in stark contrast with the particle-hole symmetric doping dependence in a square lattice. To shed light on the role of itinerant dopants in the emergence of these new magnetic properties, we furthermore leverage the unique single-particle resolution capabilities of our platform to probe higher-order density and spin correlations. These measurements hint at a new type of magnetism induced by kinetic frustration, which manifests as local antiferromagnetic correlations around hole dopants and and ferromagnetic correlations around particle dopants. Unlike magnetism induced by exchange or superexchange interactions, the energy associated with this mechanism is the kinetic energy, closely related to the famous rigorous results by Nagaoka. Furthermore, the dynamical tunability of the lattice allows adibatic engineering of quantum states in optical lattices. We discuss the preliminery efforts to prepare strongly-correlated states with significantly reduced temperatures, which may pave the way towards low temperature states in the Hubbard model.
While the exact phase diagram of the Fermi-Hubbard model remains poorly understood despite decades of progress, nearly 60 years ago Nagaoka proved that a single dopant in an otherwise half-filled Hubbard system can bring about ferromagnetism through kinetic means. The phenomenon was recently observed with ultracold atoms in triangular optical lattices. Here, we explore the kinetic ferromagnetism within the square lattice Hubbard model and its strong-coupling counterpart, the $t\ensuremath{-}J$ model, at finite temperatures in the thermodynamic limit via numerical linked-cluster expansions. We find evidence of ferromagnetic Nagaoka polarons at dopings up to $\ensuremath{\sim}30%$ away from half filling for a variety of interaction strengths and at temperatures as low as 0.2 of the hopping energy. We map out the boundaries of this phase through analyzing various correlation functions.
Understanding ferromagnetism mechanism in doped Mott insulators on frustrated lattices remains challenging at intermediate coupling and finite doping. Here, we study the itinerant ferromagnetism and propose its mechanism in doped Mott insulators on a geometrically frustrated triangular lattice. Using large-scale density matrix renormalization group (DMRG) and unrestricted Hartree–Fock mean-field methods, we reveal that itinerant ferromagnetism appears at intermediate coupling (10 ≲ U ≪ ∞) near 50% electron doping in the triangular-lattice Hubbard model. By analyzing all microscopic hopping processes, we find that doublon-singlon exchange alone drives the fully polarized ferromagnetism and uncovers the particle-hole asymmetry. We also establish the magnetic phase diagram and compare local spin correlations with recent experiments. Random phase approximation and DMRG calculations consistently confirm that the ferromagnetism persists when SU(2) symmetry is explicitly broken by magnetic anisotropy. These results clarify a microscopic route to itinerant ferromagnetism at intermediate coupling and finite doping in doped Mott insulators. Understanding the mechanism of ferromagnetism in strongly correlated systems is an ongoing theoretical challenge and the Hubbard model is typically adopted to investigate such systems. Here, the authors numerically investigate the doped triangular Fermi-Hubbard model and show ferromagnetism at intermediate coupling and finite doping, which results from the interplay of kinetic doublon-singlon exchange and lattice geometry.
We study the magnetic properties of itinerant quantum magnetic particles, described by a generalized Hubbard model with large spin (S>1/2)(S>1/2)(S>1/2), which may be realized in optical lattices of laser-cooled atom systems. In fermion systems (half-integer spins), an extended form of Nagaoka ferromagnetism may be realized. However, as novel aspects of the large-spin cases, we found that the condition on the lattice connectivity is more stringent than in the case of S=1/2S=1/2S=1/2 particles and that the system shows a peculiar degenerate structure of the ground state in which the ferromagnetic state is included. In contrast, it turns out that the ground state of itinerant bosonic systems (integer spins) has a degenerate structure similar to that of fermion system with S>1/2S>1/2S>1/2 regardless of the shape, connectivity, or filling of the lattice, and that the state with the maximum total spin is always one of the ground states. Because the system consists of 2S+12S+12S+1 types of particles and we study a SU(2S+1)\text{SU}(2S+1)SU(2S+1) invariant model, the degeneracy of the ground state is given by the multiplets of the fully symmetric Young tableau of SU(2S+1)\text{SU}(2S+1)SU(2S+1) if the state with maximum total spin belongs to the ground state.
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.