Abstract:Nagaoka ferromagnetism (NF) is a long-predicted example of itinerant ferromagnetism (IF) in the Hubbard model that has been studied theoretically for many years. The condition for NF, an infinite on-site Coulomb repulsion and a single hole in a half-filled band, does not arise naturally in materials. NF was only realized recently for the first time in experiments on a $2\times 2$ array of gated quantum dots. Dopant arrays and gated quantum dots in Si allow for engineering controllable systems with complex geometries. This makes quantum dot arrays good candidates to study NF in different geometries through analog quantum simulation. Here we present theoretical simulations done for $3\times 3$ arrays and larger $N\times N$ arrays, and predict the emergence of different forms of ferromagnetism in different geometries. We find NF in perfect $3\times 3$ arrays, and in $N\times N$ arrays for one hole doping of a half-filled band. The ratio of the hopping $t$ to Hubbard on-site repulsion $U$ that defines the onset of NF scales as $1/N^{3.5}$ as $N$ increases, approaching the bulk limit of infinite $U$ for large $N$. Additional simulations are done for geometries made by removing sites from $N\times N$ arrays. Different forms of ferromagnetism are found for different geometries. Loops show ferromagnetism, but only for three electrons. For loops, the critical $t/U$ for the onset of ferromagnetism scales as $N$ as the loop length increases. We show that the different dependences on size for loops and $N\times N$ arrays can be understood by scaling arguments that highlight the different energy contributions to each form of ferromagnetism. Our results show how analog quantum simulation with small arrays can elucidate the role of effects including wavefunction connectivity; system geometry, size and symmetry; bulk and edge sites; and kinetic energy in determining quantum magnetism of small systems.

Stability Of The Nagaoka State In The One-Band Hubbard-Model

The Nagaoka state and its stability in the one-band Hubbard model

Lattice dependence of saturated ferromagnetism in the Hubbard model

Low-Density Ferromagnetism In The Hubbard Model

The Nagaoka State in the One-Hand Hubbard Model with Two and More Holes

THE NAGAOKA STATE IN THE ONE-BAND HUBBARD-MODEL WITH 2 AND MORE HOLES

Nagaoka Ferromagnetism in Two-Band Hubbard Model

Nagaoka ferromagnetism in doped Hubbard models in optical lattices

Stability of ferromagnetism in the Hubbard model on the kagomé lattice

Polaronic mechanism of Nagaoka ferromagnetism in Hubbard models

A finite temperature version of the Nagaoka--Thouless theorem in the $\mathrm{SU}(n)$ Hubbard model

Ferrimagnetic Long-Range Order of the Hubbard Model

Superconductivity in the Two-Orbital Hubbard Model of Infinite-Layer Nickelates

Existence of Saturated Ferromagnetic and Spiral States in 1D Lieb-Ferrimagnetic Models away from Half-Filling

Numerical observation of $\mathrm{SU}(N)$ Nagaoka ferromagnetism

Stability of Ferromagnetism in Hubbard Models with Nearly-Flat Bands

Stability of the 1/3 magnetization plateau of the $J_1-J_2$ kagome Heisenberg model

Transition from Band Insulator to Mott Insulator and Formation of Local Moment in the Half-Filled Ionic SU(N) Hubbard Model

Ferromagnetism in Hubbard Models: Low Density Route

Hartree-Fock variational bounds for ground state energy of chargeless fermions with finite magnetic moment in the presence of a hard core potential: A stable ferromagnetic state

Nagaoka Ferromagnetism in $3 \times 3$ Arrays and Beyond