Range of biquadratic and triquadratic Heisenberg effective couplings deduced from multiorbital Hubbard models
Abstract
We studied a multi-orbital Hubbard model at half-filling for two and three orbitals per site on a two-site cluster via full exact diagonalization, in a wide range for the onsite repulsion U, from weak to strong coupling, and multiple ratios of the Hund coupling JH to U. The hopping matrix elements among the orbitals were also varied extensively. At intermediate and large U, we mapped the results into a Heisenberg model. For two orbitals per site, the mapping is into a S=1 Heisenberg model where by symmetry both nearest-neighbor (Si·Sj) and (Si·Sj)2 are allowed, with respective couplings J1 and J2. For the case of three orbitals per site, the mappping is into a S=3/2 Heisenberg model with (Si·Sj), (Si·Sj)2, and (Si·Sj)3 terms, and respective couplings J1, J2, and J3. The strength of these coupling constants in the Heisenberg models depend on the U, JH, and hopping amplitudes of the underlying Hubbard model. Our study allows to establish bounds on how large the ratios J2/J1 and J3/J1 can be. We show that those ratios are severely limited and, as a crude guidance, we conclude that J2/J1 is less than 0.4 and J3/J1 is less than 0.2, establishing bounds on effective models for strongly correlated Hubbard systems.
Turn this paper into a full lesson
ArcXiv compiles a staged curriculum from this paper: 8-12 lessons across beginner → advanced, synthesised section guides, visuals, flashcards, a quiz, exercises, and on-demand deep dives per section. Grounded in the abstract, never invented.