© 2026 Aqora Quantum S.A.S.

Terms Privacy Legal Notice

arXiv:2112.11346 · New Journal of Physics

Estimation of biquadratic and bicubic Heisenberg effective couplings from multiorbital Hubbard models

Rahul Soni, Nitin Kaushal, Cengiz Şen, Fernando Agustín Reboredo, Adriana Moreo, Elbio R. Dagotto

Published Jun 23, 2022• 16 citations

Abstract

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 J H 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 (S i ⋅ S j ) and ( S i ⋅ S j ) 2 are allowed, with respective couplings J 1 and J 2. For the case of three orbitals per site, the mapping is into a S = 3/2 Heisenberg model with (S i ⋅ S j ), ( S i ⋅ S j ) 2 , and ( S i ⋅ S j ) 3 terms, and respective couplings J 1, J 2, and J 3. The strength of these coupling constants in the Heisenberg models depend on the U, J H, and hopping amplitudes of the underlying Hubbard model. Our study provides a first crude estimate to establish bounds on how large the ratios J 2/J 1 and J 3/J 1 can be. We show that those ratios appear rather limited and, as a qualitative guidance, we conclude that J 2/J 1 is less than 0.4 and J 3/J 1 is less than 0.2, establishing bounds on effective models for strongly correlated Hubbard systems. Moreover, the intermediate Hubbard U regime was found to be the most promising to enhance J 2/J 1 and J 3/J 1.

DOI: 10.1088/1367-2630/ac7b9c

Related papers

Exact solutions of frustrated ordinary and chiral eight-site Hubbard models

James K. Freericks, L. M. Falicov, Daniel S. Rokhsar

This contribution addresses the question of whether or not the ground state of a frustrated spin-1/2 Heisenberg model can be smoothly related to the ground state of a simple tight-binding model at half-filling in an appropriately chosen magnetic field. This continuity is considered explicitly for an eight-site square-lattice Hubbard model with nearest- and next-nearest-neighbor hopping, which approaches a frustrated Heisenberg model with nearest- and next-nearest-neighbor antiferromagnetic exchange coupling in the limit of large Hubbard U. In addition, a counterexample to the widely held belief that the half-filled band of the ordinary Hubbard model is nonmagnetic has been found in the regime where the hopping parameters and Hubbard interaction are all the same order of magnitude.

Physical review. B, Condensed matter•Jul 15, 1991•17 citations

Two-orbital Hubbard model vs spin S=1 Heisenberg model: studies on clusters

R. Lemański, Matysiak

We perform exact numeric calculations for the two-orbital Hubbard model on the four-site cluster.In the limit of large on-site coupling the model becomes equivalent to the spin S " 1 Heisenberg model.By comparing energy spectra of these two models, we quantified the range of interaction parameters for which the Heisenberg model satisfactorily reproduces the two-orbital Hubbard model.Then we examined how the spectrum evolves when we are outside of this region, focusing especially on checking of how it is modified when various ways of interatomic hoppings of electrons between different orbitals are taken into account.We finally show how these modifications affect the dependence of specific heat on temperature.

arXiv:1809.09483•Condensed Matter Physics•Sep 1, 2018•3 citations

Multisite versus multiorbital Coulomb correlations studied within finite-temperature exact diagonalization dynamical mean-field theory

A. Liebsch, H. Ishida, Jaime Merino

The influence of short-range Coulomb correlations on the Mott transition in the single-band Hubbard model at half filling is studied within cellular dynamical mean-field theory for square and triangular lattices. Finite-temperature exact diagonalization is used to investigate correlations within two-, three-, and four-site clusters. Transforming the nonlocal self-energy from a site basis to a molecular-orbital basis, we focus on the interorbital charge transfer between these cluster molecular orbitals in the vicinity of the Mott transition. In all cases studied, the charge transfer is found to be small, indicating weak Coulomb-induced orbital polarization despite sizable level splitting between orbitals. These results demonstrate that all cluster molecular orbitals take part in the Mott transition and that the insulating gap opens simultaneously across the entire Fermi surface. Thus, at half filling we do not find orbital-selective Mott transitions or a combination of band filling and Mott transition in different orbitals. Nevertheless, the approach toward the transition differs greatly between cluster orbitals, giving rise to a pronounced momentum variation along the Fermi surface, in agreement with previous works. The near absence of Coulomb-induced orbital polarization in these clusters differs qualitatively from single-site multiorbital studies of several transition-metal oxides, where the Mott phase exhibits nearly complete orbital polarization as a result of a correlation driven enhancement of the crystal-field splitting. The strong single-particle coupling among cluster orbitals in the single-band case is identified as the source of this difference.

arXiv:0810.4837•Physical Review B•Oct 28, 2008•50 citations

Comparison of the Multi-Band Hubbard and Heisenberg Models on Small Finite Clusters

J. Matysiak, R. Lemański

We perform exact numeric calculations for a few versions of the two-orbital Hubbard model on the four site cluster.We show that the conventional spin 1 Heisenberg Hamiltonian can be obtained in the strong coupling limit, but only for some versions of the two-orbital Hubbard model.Then we propose a modified version of the Heisenberg Hamiltonian, that we call the two-orbital Heisenberg model, which is relevant in those cases, where the conventional spin 1 Heisenberg model fails.

Acta Physica Polonica A•Mar 1, 2018•1 citation

Magnetism in the Hubbard model: An effective spin Hamiltonian approach

Michael A. Tusch, Yolande H. Szczech, David E. Logan

We present an approach to the magnetic properties of the half-filled Hubbard model, based on an approximate mapping of its low-energy transverse spin excitations on to those of an effective underlying Heisenberg model, but with effective spin interactions which are self-consistently determined and not confined solely to nearest-neighbor couplings. The mapping is exact in strong-coupling and is found to be accurate over a very wide range of interaction strengths, down to weak coupling. At zero temperature, it permits ready evaluation at finite U of the one-loop effects of zero-point spin fluctuations on, e.g., the sublattice magnetization. At finite temperatures, thermodynamic properties of the system in the thermal paramagnet are studied via a physically transparent Onsager reaction field approach, which amounts to a self-consistent treatment of paramagnetic spin correlations. This is central not only in recovering the correct dimensionality dependence of antiferromagnetic long-ranged order, but also for the d=3 case of primary interest here yields a N\'eel temperature in close agreement with known strong- and weak-coupling limits. Spin correlation functions and magnetic susceptibilities also show very good agreement with quantum Monte Carlo calculations over an appreciable temperature range in which the low-lying transverse spin excitations are thermally dominant. \textcopyright{} 1996 The American Physical Society.

Physical review. B, Condensed matter•Mar 1, 1996•16 citations

Configuration-interaction approach to hole pairing in the two-dimensional Hubbard model

Enrique Louis, F. Guinea, María P. López-Sancho +1 more

The interactions between holes in the Hubbard model, in the low density, intermediate- to strong-coupling limit, are investigated by systematically improving mean field calculations. The configuration-interaction basis set is constructed by applying to local unrestricted Hartree-Fock configurations all lattice translations and rotations. It is shown that this technique reproduces, correctly, the properties of the Heisenberg model, in the limit of large U. Upon doping, dressed spin polarons in neighboring sites have an increased kinetic energy and an enhanced hopping rate. Both effects are of the order of the hopping integral and lead to an effective attraction at intermediate couplings. The numerical results also show that when more than two holes are added to the system, they do not tend to cluster, but rather hole pairs remain far apart. Hole-hole correlations are also calculated and shown to be in qualitative agreement with exact calculations for $4\ifmmode\times\else\texttimes\fi{}4$ clusters. In particular our results indicate that for intermediate coupling the hole-hole correlation attains a maximum when the holes are in the same sublattice at a distance of $\sqrt{2}$ times the lattice spacing, in agreement with exact results and the $t\ensuremath{-}J$ model. The method is also used to derive known properties of the quasiparticle band structure of isolated spin polarons.

arXiv:cond-mat/9811298•Physical review. B, Condensed matter•Jun 1, 1999•21 citations

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 J H 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 (S i ⋅ S j ) and ( S i ⋅ S j ) 2 are allowed, with respective couplings J 1 and J 2. For the case of three orbitals per site, the mapping is into a S = 3/2 Heisenberg model with (S i ⋅ S j ), ( S i ⋅ S j ) 2 , and ( S i ⋅ S j ) 3 terms, and respective couplings J 1, J 2, and J 3. The strength of these coupling constants in the Heisenberg models depend on the U, J H, and hopping amplitudes of the underlying Hubbard model. Our study provides a first crude estimate to establish bounds on how large the ratios J 2/J 1 and J 3/J 1 can be. We show that those ratios appear rather limited and, as a qualitative guidance, we conclude that J 2/J 1 is less than 0.4 and J 3/J 1 is less than 0.2, establishing bounds on effective models for strongly correlated Hubbard systems. Moreover, the intermediate Hubbard U regime was found to be the most promising to enhance J 2/J 1 and J 3/J 1.

This contribution addresses the question of whether or not the ground state of a frustrated spin-1/2 Heisenberg model can be smoothly related to the ground state of a simple tight-binding model at half-filling in an appropriately chosen magnetic field. This continuity is considered explicitly for an eight-site square-lattice Hubbard model with nearest- and next-nearest-neighbor hopping, which approaches a frustrated Heisenberg model with nearest- and next-nearest-neighbor antiferromagnetic exchange coupling in the limit of large Hubbard U. In addition, a counterexample to the widely held belief that the half-filled band of the ordinary Hubbard model is nonmagnetic has been found in the regime where the hopping parameters and Hubbard interaction are all the same order of magnitude.

We perform exact numeric calculations for the two-orbital Hubbard model on the four-site cluster.In the limit of large on-site coupling the model becomes equivalent to the spin S " 1 Heisenberg model.By comparing energy spectra of these two models, we quantified the range of interaction parameters for which the Heisenberg model satisfactorily reproduces the two-orbital Hubbard model.Then we examined how the spectrum evolves when we are outside of this region, focusing especially on checking of how it is modified when various ways of interatomic hoppings of electrons between different orbitals are taken into account.We finally show how these modifications affect the dependence of specific heat on temperature.

The influence of short-range Coulomb correlations on the Mott transition in the single-band Hubbard model at half filling is studied within cellular dynamical mean-field theory for square and triangular lattices. Finite-temperature exact diagonalization is used to investigate correlations within two-, three-, and four-site clusters. Transforming the nonlocal self-energy from a site basis to a molecular-orbital basis, we focus on the interorbital charge transfer between these cluster molecular orbitals in the vicinity of the Mott transition. In all cases studied, the charge transfer is found to be small, indicating weak Coulomb-induced orbital polarization despite sizable level splitting between orbitals. These results demonstrate that all cluster molecular orbitals take part in the Mott transition and that the insulating gap opens simultaneously across the entire Fermi surface. Thus, at half filling we do not find orbital-selective Mott transitions or a combination of band filling and Mott transition in different orbitals. Nevertheless, the approach toward the transition differs greatly between cluster orbitals, giving rise to a pronounced momentum variation along the Fermi surface, in agreement with previous works. The near absence of Coulomb-induced orbital polarization in these clusters differs qualitatively from single-site multiorbital studies of several transition-metal oxides, where the Mott phase exhibits nearly complete orbital polarization as a result of a correlation driven enhancement of the crystal-field splitting. The strong single-particle coupling among cluster orbitals in the single-band case is identified as the source of this difference.

We present an approach to the magnetic properties of the half-filled Hubbard model, based on an approximate mapping of its low-energy transverse spin excitations on to those of an effective underlying Heisenberg model, but with effective spin interactions which are self-consistently determined and not confined solely to nearest-neighbor couplings. The mapping is exact in strong-coupling and is found to be accurate over a very wide range of interaction strengths, down to weak coupling. At zero temperature, it permits ready evaluation at finite U of the one-loop effects of zero-point spin fluctuations on, e.g., the sublattice magnetization. At finite temperatures, thermodynamic properties of the system in the thermal paramagnet are studied via a physically transparent Onsager reaction field approach, which amounts to a self-consistent treatment of paramagnetic spin correlations. This is central not only in recovering the correct dimensionality dependence of antiferromagnetic long-ranged order, but also for the d=3 case of primary interest here yields a N\'eel temperature in close agreement with known strong- and weak-coupling limits. Spin correlation functions and magnetic susceptibilities also show very good agreement with quantum Monte Carlo calculations over an appreciable temperature range in which the low-lying transverse spin excitations are thermally dominant. \textcopyright{} 1996 The American Physical Society.

The interactions between holes in the Hubbard model, in the low density, intermediate- to strong-coupling limit, are investigated by systematically improving mean field calculations. The configuration-interaction basis set is constructed by applying to local unrestricted Hartree-Fock configurations all lattice translations and rotations. It is shown that this technique reproduces, correctly, the properties of the Heisenberg model, in the limit of large U. Upon doping, dressed spin polarons in neighboring sites have an increased kinetic energy and an enhanced hopping rate. Both effects are of the order of the hopping integral and lead to an effective attraction at intermediate couplings. The numerical results also show that when more than two holes are added to the system, they do not tend to cluster, but rather hole pairs remain far apart. Hole-hole correlations are also calculated and shown to be in qualitative agreement with exact calculations for $4\ifmmode\times\else\texttimes\fi{}4$ clusters. In particular our results indicate that for intermediate coupling the hole-hole correlation attains a maximum when the holes are in the same sublattice at a distance of

\sqrt{2}

times the lattice spacing, in agreement with exact results and the $t\ensuremath{-}J$ model. The method is also used to derive known properties of the quasiparticle band structure of isolated spin polarons.