Optimal array geometries for kinetic magnetism and Nagaoka polarons
Phys. Rev. B 113, 075156 – Published 26 February, 2026
DOI: https://doi.org/10.1103/51bm-vxbn
Abstract
Quantum dot platforms have enabled the direct observation of Nagaoka ferromagnetism (NFM) in small arrays and noninfinite interaction strength. However, optimizing the cluster connectivity characteristics that yield a ground state with maximal spin and their robustness against magnetic fields remains unexplored. Employing exact diagonalization of the Hubbard Hamiltonian, we find a connection between the existence of kinetic ferromagnetism and graph theory descriptions. Algebraic connectivity () and Katz centrality (KC) are shown to be related to the spin correlation over the system. In square arrays, the onset of NFM is found to be . In optimal cluster geometries, large and low KC fluctuation per site are found to enhance , extending the NFM phase while diminishing the strength of spin-correlation clouds. A perpendicular magnetic field introduces Aharonov-Bohm phases, and a critical flux for which NFM is destroyed. We further find that tuning the flux phase to results in a ground state that exhibits antiferromagnetic correlations (counter-Nagaoka state). Our results illustrate how NFM and polaron formation can be predicted from the array's connectivity ( and KC), and how the introduction of flux results in the counterintuitive destruction of kinetic ferromagnetism in the system.