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

Ferromagnetic domains in the large-$U$ Hubbard model with a few holes: an FCIQMC study

Yun, Sujun, Werner Dobrautz, Luo Hongjun, Vamshi M. Katukuri, Niklas Libermann, Ali Alavi

Published Oct 21, 2022• 1 citation
arXivHTMLPDF

Abstract

Two-dimensional Hubbard lattices with two or three holes are investigated as a function of $U$ in the large-$U$ limit. In the so-called Nagaoka limit (one-hole system at infinite $U$), it is known that the Hubbard model exhibits a ferromagnetic ground-state. Here, by means of exact FCIQMC simulations applied to periodic lattices up to 24 sites, we compute spin-spin correlation functions as a function of increasing $U$. The correlation functions clearly demonstrate the onset of ferromagnetic domains, centred on individual holes. The overall total spin of the wavefunctions remain the lowest possible (0, or $\frac{1}{2}$, depending on the number of holes). The ferromagnetic domains appear at interaction strengths comparable to the critical interaction strengths of the Nagaoka transition in finite systems with strictly one hole. The existence of such ferromagnetic domains is the signature of Nagaoka physics in Hubbard systems with a small (but greater than one) number of holes.

DOI
10.48550/arxiv.2210.12049

Related papers

Ferromagnetic domains in the large- U Hubbard model with a few holes: A full configuration interaction quantum Monte Carlo study

Su-Jun Yun, Werner Dobrautz, Hongjun Luo +3 more

Two-dimensional Hubbard lattices with two or three holes are investigated as a function of $U$ in the large-$U$ limit. In the so-called Nagaoka limit (one-hole system at infinite $U$), it is known that the Hubbard model exhibits a ferromagnetic ground state. Here, by means of exact full configuration interaction quantum Monte Carlo simulations applied to periodic lattices up to 24 sites, we compute spin-spin correlation functions as a function of increasing $U$. The correlation functions clearly demonstrate the onset of ferromagnetic domains, centered on individual holes. The overall total spin of the wave functions remains the lowest possible (0 or $\frac{1}{2}$, depending on the number of holes). The ferromagnetic domains appear at interaction strengths comparable to the critical interaction strengths of the Nagaoka transition in finite systems with strictly one hole. The existence of such ferromagnetic domains is the signature of Nagaoka physics in Hubbard systems with a small (but greater than 1) number of holes.

Physical review. B./Physical review. B•Feb 7, 2023•8 citations

Benchmark study of Nagaoka ferromagnetism by spin-adapted full configuration interaction quantum Monte Carlo

Sujun Yun, Werner Dobrautz, Hongjun Luo +1 more

We investigate Nagaoka ferromagnetism in the two-dimensional Hubbard model with one hole using the spin-adapted [$\text{SU}(2)$ conserving] full configuration interaction quantum Monte Carlo method. This methodology gives us access to the ground-state energies of all possible spin states $S$ of finite Hubbard lattices, here obtained for lattices up to 26 sites for various interaction strengths ($U$). The critical interaction strength, ${U}_{c}$, at which the Nagaoka transition occurs is determined for each lattice and is found to be proportional to the lattice size for the larger lattices. Below ${U}_{c}$, the overall ground states are found to favour the minimal total spin ($S=\frac{1}{2}$), and no intermediate spin state is found to be the overall ground state on lattices larger than 16 sites. However, at ${U}_{c}$, the energies of all the spin states are found to be nearly degenerate, implying that large fluctuations in total spin can be expected in the vicinity of the Nagaoka transition.

arXiv:2105.06802•Physical review. B./Physical review. B•Dec 2, 2021•13 citations

Holes in the infinite-UHubbard model: Instability of the Nagaoka state

Yan Yan Fang, Andrei E. Ruckenstein, Elbio R. Dagotto +1 more

The solution of a three-body Schr\"odinger equation (for two holes and a spin flip) shows that in the U=\ensuremath{\infty} Hubbard model the saturated ferromagnetic state with two holes is unstable (in a finite volume and for periodic boundary conditions). The ground state carries a finite momentum, and becomes degenerate with the Nagaoka state in the thermodynamic limit. This is in agreement with exact diagonalization results for small lattices where, in addition, we find that the true two-hole ground state is a singlet. As a result, our notions of ferromagnetism in the Hubbard model might need revision.

Physical review. B, Condensed matter•Oct 1, 1989•79 citations

Infinite-U Hubbard model in the large-spin regime: exact diagonalization study

J C Anglès d'Auriac, Benoît Douçot, R. Rammal

The infinite-U Hubbard model on finite two-dimensional square lattices is investigated by the use of exact numerical diagonalization. Configurations with one spin flip away from the ferromagnetic state and up to three holes are studied for lattice sizes ranging from 4*4 to 12*12. For all these situations, the state with one spin flip is found to be more stable than the Nagaoka ferromagnetic state. The crystal momentum of the ground state in this spin sector depends on the system size. Systematic comparison with the Gutzwiller wavefunction shows that it reproduces qualitatively the crystal momentum dependence of the energy, but not the exact correlation between the spin flip and the holes.

Journal of Physics Condensed Matter•Jun 3, 1991•4 citations

Nagaoka Ferromagnetism in the Two-Dimensional Infinite- U Hubbard Model

Federico Becca, Sandro Sorella

We present different numerical calculations based on variational quantum Monte Carlo simulations supporting a ferromagnetic ground state for finite and small hole densities in the two-dimensional infinite-U Hubbard model. Moreover, by studying the energies of different total spin sectors, these calculations strongly suggest that the paramagnetic phase is unstable against a phase with a partial polarization for large hole densities delta approximately 0.40 with evidence for a second-order transition to the paramagnetic large doping phase.

arXiv:cond-mat/0102249•Physical Review Letters•Apr 9, 2001•58 citations

Hole dynamics in the two-dimensional strong-coupling Hubbard Hamiltonian

Efthimios Kaxiras, Efstratios Manousakis

The effective Hamiltonian, obtained from the Hubbard model in the strong-coupling limit, is diagonalized exactly for a periodic two-dimensional square lattice of ten sites. We obtain the ground state of the system for any filling of the lattice and the $K$-space excitation spectrum for a single hole. Within our finite-size system, it is found that for very small coupling ratio $\frac{t}{U}$, each hole creates a local ferromagnetic environment at least of the size of our system, whereas for increasing $\frac{t}{U}$, the holes tend to form clusters, induced by antiferromagnetic spin correlations. A range of $\frac{t}{U}$ exists for which pairing of the holes may be possible.

Physical review. B, Condensed matter•Jul 1, 1988•66 citations
Ferromagnetic domains in the large-UUU Hubbard model with a few holes: an FCIQMC study
Two-dimensional Hubbard lattices with two or three holes are investigated as a function of UUU in the large-UUU limit. In the so-called Nagaoka limit (one-hole system at infinite UUU), it is known that the Hubbard model exhibits a ferromagnetic ground-state. Here, by means of exact FCIQMC simulations applied to periodic lattices up to 24 sites, we compute spin-spin correlation functions as a function of increasing UUU. The correlation functions clearly demonstrate the onset of ferromagnetic domains, centred on individual holes. The overall total spin of the wavefunctions remain the lowest possible (0, or 12\frac{1}{2}21​, depending on the number of holes). The ferromagnetic domains appear at interaction strengths comparable to the critical interaction strengths of the Nagaoka transition in finite systems with strictly one hole. The existence of such ferromagnetic domains is the signature of Nagaoka physics in Hubbard systems with a small (but greater than one) number of holes.
Two-dimensional Hubbard lattices with two or three holes are investigated as a function of UUU in the large-UUU limit. In the so-called Nagaoka limit (one-hole system at infinite UUU), it is known that the Hubbard model exhibits a ferromagnetic ground state. Here, by means of exact full configuration interaction quantum Monte Carlo simulations applied to periodic lattices up to 24 sites, we compute spin-spin correlation functions as a function of increasing UUU. The correlation functions clearly demonstrate the onset of ferromagnetic domains, centered on individual holes. The overall total spin of the wave functions remains the lowest possible (0 or 12\frac{1}{2}21​, depending on the number of holes). The ferromagnetic domains appear at interaction strengths comparable to the critical interaction strengths of the Nagaoka transition in finite systems with strictly one hole. The existence of such ferromagnetic domains is the signature of Nagaoka physics in Hubbard systems with a small (but greater than 1) number of holes.
We investigate Nagaoka ferromagnetism in the two-dimensional Hubbard model with one hole using the spin-adapted [SU(2)\text{SU}(2)SU(2) conserving] full configuration interaction quantum Monte Carlo method. This methodology gives us access to the ground-state energies of all possible spin states SSS of finite Hubbard lattices, here obtained for lattices up to 26 sites for various interaction strengths (UUU). The critical interaction strength, Uc{U}_{c}Uc​, at which the Nagaoka transition occurs is determined for each lattice and is found to be proportional to the lattice size for the larger lattices. Below Uc{U}_{c}Uc​, the overall ground states are found to favour the minimal total spin (S=12S=\frac{1}{2}S=21​), and no intermediate spin state is found to be the overall ground state on lattices larger than 16 sites. However, at Uc{U}_{c}Uc​, the energies of all the spin states are found to be nearly degenerate, implying that large fluctuations in total spin can be expected in the vicinity of the Nagaoka transition.
The solution of a three-body Schr\"odinger equation (for two holes and a spin flip) shows that in the U=\ensuremath{\infty} Hubbard model the saturated ferromagnetic state with two holes is unstable (in a finite volume and for periodic boundary conditions). The ground state carries a finite momentum, and becomes degenerate with the Nagaoka state in the thermodynamic limit. This is in agreement with exact diagonalization results for small lattices where, in addition, we find that the true two-hole ground state is a singlet. As a result, our notions of ferromagnetism in the Hubbard model might need revision.
The infinite-U Hubbard model on finite two-dimensional square lattices is investigated by the use of exact numerical diagonalization. Configurations with one spin flip away from the ferromagnetic state and up to three holes are studied for lattice sizes ranging from 4*4 to 12*12. For all these situations, the state with one spin flip is found to be more stable than the Nagaoka ferromagnetic state. The crystal momentum of the ground state in this spin sector depends on the system size. Systematic comparison with the Gutzwiller wavefunction shows that it reproduces qualitatively the crystal momentum dependence of the energy, but not the exact correlation between the spin flip and the holes.
We present different numerical calculations based on variational quantum Monte Carlo simulations supporting a ferromagnetic ground state for finite and small hole densities in the two-dimensional infinite-U Hubbard model. Moreover, by studying the energies of different total spin sectors, these calculations strongly suggest that the paramagnetic phase is unstable against a phase with a partial polarization for large hole densities delta approximately 0.40 with evidence for a second-order transition to the paramagnetic large doping phase.
The effective Hamiltonian, obtained from the Hubbard model in the strong-coupling limit, is diagonalized exactly for a periodic two-dimensional square lattice of ten sites. We obtain the ground state of the system for any filling of the lattice and the KKK-space excitation spectrum for a single hole. Within our finite-size system, it is found that for very small coupling ratio tU\frac{t}{U}Ut​, each hole creates a local ferromagnetic environment at least of the size of our system, whereas for increasing tU\frac{t}{U}Ut​, the holes tend to form clusters, induced by antiferromagnetic spin correlations. A range of tU\frac{t}{U}Ut​ exists for which pairing of the holes may be possible.