---
title: Bilayer Two-Orbital Kanamori-Hubbard Model
url: https://www.emergentmind.com/topics/bilayer-two-orbital-kanamori-hubbard-model
type: topic
---

# Bilayer Two-Orbital Kanamori-Hubbard Model

The bilayer two-orbital Kanamori-Hubbard model denotes a family of correlated-electron lattice theories in which bilayer geometry and two active orbital flavors are treated on comparable footing, while local interactions are organized in the Kanamori form or in controlled truncations of it. In current usage, the term covers several distinct but related constructions: explicit bilayer square-lattice \(e_g\) models with layer labels \(\alpha=1,2\) and orbitals \(d_{x^2-y^2}\) and \(d_{z^2}\); strong-coupling descendants such as bilayer \(t\)-\(J\) and Kugel-Khomskii models; and bilayer-derived moiré theories in which the microscopic layer degree of freedom has been integrated into effective Wannier orbitals, so that the final Hamiltonian is no longer an explicit bilayer lattice model even though its origin is bilayer physics [2601.06440][2508.04554][2303.02305].

## 1. Scope of the model class

In the explicit nickelate-oriented formulations, each in-plane lattice site carries two layers and two orbitals per layer, so the local flavor structure before spin is \((x_1,z_1,x_2,z_2)\), with \(x\equiv d_{x^2-y^2}\) and \(z\equiv d_{z^2}\) or \(d_{3z^2-r^2}\). This is the setting adopted in bilayer square-lattice models motivated by \(\mathrm{La_3Ni_2O_7}\), where interlayer \(d_{z^2}\)-\(d_{z^2}\) hopping and intralayer \(d_{x^2-y^2}\)-\(d_{z^2}\) hybridization are central one-body ingredients [2312.03605][2506.07741].

A narrower usage reserves the phrase for models with a fully rotationally invariant onsite Kanamori interaction. In that strict sense, the quarter-hole-filled bilayer \(e_g\) model studied in a molecular-orbital basis is an explicit bilayer two-orbital Kanamori-Hubbard model: it contains two layers \(T,B\), two orbitals \(x,z\), crystal-field splitting, and onsite \(U\), \(U'\), Hund exchange, spin-flip, and pair-hopping terms with \(U'=U-2J_{\rm H}\) [2601.06440]. The same strict classification applies to the bilayer nickelate model used to study two intertwined \(s_{\pm}\)-wave superconductivities, which employs the onsite multiorbital Kanamori form with \(U=7t\) and \(J_H\) scanned over \(0.1U\sim 0.3U\) [2508.04554].

A broader usage includes density-density or longitudinal-Hund truncations. This is the case for the bilayer two-orbital Hubbard model studied by cluster DMFT for nickelate superconductivity, where the main calculations retain only \(U_z\) on the \(d_{z^2}\) orbital and later add \(U_x\), \(U_{xz}\), and a longitudinal Hund coupling \(J_H\), but no spin-flip or pair-hopping [2312.03605]. It is also the case for realistic DCA and DMRG studies of bilayer nickelates that keep \(U\), \(U'\), and density-density Hund terms but omit the full rotational structure for sign-problem or tractability reasons [2506.07741][2511.01801].

An important limiting case is the bilayer-derived moiré construction. In twisted \(\Gamma\)-valley TMD homobilayers, Wannierization of the continuum model produces a two-orbital honeycomb theory with \(p_x\)- and \(p_y\)-like Wannier states on the moiré sublattices. The microscopic origin is bilayer, but the low-energy Hamiltonian is not an explicit bilayer lattice Hamiltonian; instead, it is an effective honeycomb model with two orbitals per site and a local sector that resembles Kanamori interactions while also containing substantial nonlocal Coulomb, exchange, pair-hopping, and interaction-assisted hopping terms. For that reason the authors term it a Kanamori-Moiré-Hubbard model rather than a standard bilayer Kanamori-Hubbard model [2303.02305].

## 2. Interaction structure and the meaning of “Kanamori”

A canonical two-orbital Kanamori interaction is written as
\[
H_{\rm int}^{\rm Kanamori} = U\sum_m n_{m\uparrow}n_{m\downarrow} + U'\sum_{m<n,\sigma,\sigma'} n_{m\sigma}n_{n\sigma'} - J_H\sum_{m<n} 2\,\mathbf S_m\cdot\mathbf S_n + J_H\sum_{m<n} (P^\dagger_m P_n+\text{h.c.}),
\]
often together with the rotationally invariant constraint \(U'=U-2J_H\). In the current literature, explicit bilayer two-orbital models span the full range from this form to drastic truncations [2303.02305].

| Model family | Bilayer representation | Interaction sector |
|---|---|---|
| Explicit bilayer \(e_g\) Kanamori models | Layers retained explicitly | Full local Kanamori with spin-flip and pair-hopping |
| Bilayer two-orbital density-density models | Layers retained explicitly | \(U\), \(U'\), longitudinal or Ising-like Hund terms |
| Kanamori-Moiré-Hubbard models | Bilayer encoded implicitly after Wannierization | Local Kanamori-like terms plus strong nonlocal terms |

The fully rotationally invariant case is represented by the quarter-hole-filled bilayer model in which
\[
H_{\rm int} = \sum_{i\alpha\eta}\epsilon_{\alpha} n_{i\alpha}^{\eta} +U\sum_{i,\alpha,\eta}n_{i\alpha\uparrow}^{\eta} n_{i\alpha\downarrow}^{\eta} +\sum_{i,\alpha<\beta,\sigma,\eta} \Big[ U^\prime n_{i\alpha\sigma}^{\eta}n_{i\beta\bar\sigma}^{\eta} +(U^\prime-J_{\rm H})n_{i\alpha\sigma}^{\eta}n_{i\beta\sigma}^{\eta}
\]
together with explicit spin-flip and pair-hopping terms, again with \(U'=U-2J_{\rm H}\) [2601.06440]. The Kanamori model for \(\mathrm{La_3Ni_2O_7}\) used in CDMFT likewise includes intraorbital Hubbard repulsion, interorbital density-density interactions, Hund exchange, spin-flip, and pair-hopping [2508.04554].

By contrast, several influential nickelate studies are only Kanamori-like. One adopts
\[
H_U= U_{z}\sum_{i\alpha}n_{i\alpha\uparrow}^{z}n_{i\alpha\downarrow}^{z}
\]
for the main calculations, and later extends it to
\[
H^{\prime}_U= U_{x}\sum_{i\alpha}n_{i\alpha\uparrow}^{x}n_{i\alpha\downarrow}^{x} + U_{xz}\sum_{i\alpha}n_{i\alpha \sigma }^{x}n_{i\alpha \overline{\sigma} }^{z} +(U_{xz}-J_H)\sum_{i\alpha, \sigma }n_{i\alpha\sigma}^{x}n_{i\alpha \sigma}^{z},
\]
with \(U_x=U_z=7\), \(J_H=1.4\), and \(U_{xz}=U_z-2J_H\), but explicitly without spin-flip and pair-hopping [2312.03605]. A realistic DCA study of pressurized \(\mathrm{La_3Ni_2O_7}\) similarly uses a density-density two-orbital Hubbard-Hund interaction with \(U=3~{\rm eV}\), \(U'=2~{\rm eV}\), and \(J=0.5~{\rm eV}\), consistent with \(U'=U-2J\), while omitting spin-flip and pair-hopping [2506.07741]. The zero-temperature CDMFT+NORG work on La327 goes further toward an Ising-Hund representation by separating interorbital couplings into \(U'_a\) and \(U'_p\) and identifying the local interorbital spin coupling \(\Delta U' = U'_a-U'_p\) as the density-density proxy for Hund physics [2407.13737].

The moiré case is more general still. Its intra-unit-cell interaction \(H_i\) contains the usual intra-orbital Hubbard term \(U_0\), interorbital density-density terms \(U_{\alpha\beta}\), Hund/direct-exchange terms \(J_{\alpha\beta}\), and pair-hopping, but also interaction-assisted hoppings, spin-flip hopping accompanied by local spin flip, and doublon scattering processes parameterized by \(\tilde A_{\alpha\beta\gamma}\), \(A_{\alpha\beta\gamma}\), and \(J_{\beta\alpha\gamma}\). The authors state that these additional terms are precisely what place the model beyond the standard multi-orbital Hubbard description [2303.02305].

A persistent source of confusion is therefore terminological. “Bilayer two-orbital Kanamori-Hubbard model” may refer either to a fully SU(2)-invariant local interaction with explicit layer operators, or to a broader family of bilayer two-orbital Hubbard-Hund models in which only the density-density subset is kept. The literature represented here contains both usages [2601.06440][2312.03605].

## 3. Bilayer geometry, orbital selectivity, and one-body structure

The defining one-body feature of the nickelate models is strongly orbital-selective bilayer hopping. In the minimal bilayer square-lattice formulation motivated by \(\mathrm{La_3Ni_2O_7}\), the noninteracting Hamiltonian contains intralayer \(d_{x^2-y^2}\) and \(d_{z^2}\) hoppings, interlayer hopping \(t_\perp\) only between \(d_{z^2}\) orbitals on the same in-plane site, and intralayer interorbital hybridization \(V\). One representative parametrization fixes \(t_x^1=1\), \(t_x^2=0.15\), \(t_z^1=-0.25\), \(t_z^2=-0.05\), with total filling \(n=0.75\) and typical orbital occupancies \(n_z\approx 0.45\text{--}0.47\) and \(n_x\approx 0.28\text{--}0.30\) [2312.03605].

The realistic DCA model for pressurized \(\mathrm{La_3Ni_2O_7}\) sharpens this hierarchy. Using parameters for 25 GPa, it takes \(t_{xx}=-0.515\) eV, \(t_{zz}=-0.110\) eV, \(t_{xz}=0.243\) eV, \(t^\perp_{zz}=0.666\) eV, and crystal-field splitting \(\epsilon_x=0.506\) eV at filling \(\langle n\rangle=1.5\). The largest hopping is \(t^\perp_{zz}\), so the dominant bonding-antibonding splitting is in the \(d_{3z^2-r^2}\) sector, and the low-energy superconducting physics becomes overwhelmingly \(z\)-orbital dominated despite the presence of both \(e_g\) orbitals [2506.07741].

Pressure-dependent DFT+Wannier downfolding reaches the same conclusion from a different direction. In the bilayer ladder model derived for \(\mathrm{La_3Ni_2O_7}\), \(t^{xz}_{[100]}\) increases by about \(30\%\) and \(t^{zz}_{[001]}\) by about \(37.6\%\) between 21.6 and 100 GPa, while the effective interaction-to-hopping ratio \(U/|t_\parallel^{xx}|\) decreases from \(8.0\) to \(6.6\) by 80.3 GPa. This explicitly identifies pressure as simultaneously enhancing bilayer/interorbital kinetic couplings and reducing relative correlation strength [2511.01801].

The strong-coupling bilayer Kanamori model at quarter-hole filling isolates the same orbital selectivity in molecular-orbital language. There the relevant hierarchy is
\[
t_{\perp}^{zz},\, U,\, J_H,\, |\epsilon_x-\epsilon_z| \gg t_{\parallel}^{xx},
\]
and the \(d_{z^2}\) sector is reorganized into bonding and antibonding molecular orbitals. The low-energy manifold on each rung consists of a doubly occupied bonding \(d_{z^2}\) singlet plus one remaining \(d_{x^2-y^2}\) electron, giving four states \(|S^z,\tau^z\rangle\) with \(S^z=\pm \tfrac12\) and a layer pseudospin \(\tau^z=\pm \tfrac12\) that tracks whether the \(d_{x^2-y^2}\) electron sits on the top or bottom layer [2601.06440].

The moiré realization changes the geometry but preserves the multiorbital logic. Starting from a continuum moiré Hamiltonian
\[
H = -\frac{\hbar k^2}{2m^*} + \Delta(\mathbf r),\qquad \Delta(\mathbf r)=\sum_{s=1}^{6} V_s e^{i(\mathbf g_s\cdot \mathbf r + \phi_s)},
\]
Wannierization yields four well-localized \(p_x\)- and \(p_y\)-like Wannier functions on a moiré honeycomb lattice. The resulting tight-binding description lives on a single effective honeycomb lattice with two orbitals per site and hoppings up to third nearest neighbors; the bilayer origin survives only implicitly through the continuum bands and the Wannier orbitals [2303.02305].

A plausible implication is that “bilayer” is not a unique microscopic category but a projection-dependent one. When interlayer dimerization is retained explicitly, the natural degrees of freedom are layers, orbitals, and sometimes bonding-antibonding combinations. When the same physics is downfolded into extended Wannier states, the low-energy theory may instead look like a two-orbital single-layer model, even though its parameters remain bilayer controlled [2601.06440][2303.02305].

## 4. Theoretical methods and effective reductions

The model class has been studied with a correspondingly heterogeneous methodological toolkit. Weak-to-intermediate-coupling superconductivity in explicit bilayer two-orbital lattice models has mainly been addressed by cluster DMFT and DCA. The nickelate study built around \(U_z\), \(V\), and \(t_\perp\) uses single-site DMFT, four-site CDMFT, and eight-site DCA with CTQMC and HFQMC impurity solvers, plus high-temperature DQMC benchmarking on \(12\times 12\times 4\) orbitals [2312.03605]. The realistic bilayer two-orbital Hubbard-Hund study of pressurized \(\mathrm{La_3Ni_2O_7}\) employs DCA with a continuous-time auxiliary-field QMC solver on clusters \(N_c=2\), \(8\), and \(16\) [2506.07741].

Zero-temperature normal-state correlation physics has been treated by CDMFT with a natural orbitals renormalization group impurity solver. In the La327 study this choice is tailored to a two-site bilayer cluster containing the two layers \(A,B\) and both orbitals on each site, thereby resolving interlayer bonding/antibonding structure and Hund-mediated magnetic transfer exactly at the impurity level [2407.13737]. Sign-problem-free determinant QMC becomes available when the interaction sector is simplified sufficiently; this route was used both in the half-filled bilayer two-orbital Hubbard benchmark and in a related bilayer \(U(2)\times U(2)\)-symmetric model with an Ising-like ferromagnetic interlayer coupling [2408.12042][2304.13083].

Strong-coupling reductions have been central to the subject. The quarter-hole-filled bilayer Kanamori model is reduced by a Schrieffer-Wolff transformation to an anisotropic Kugel-Khomskii Hamiltonian,
\[
H_{\text{eff}} = \sum_{ij}J_{s}\,\hat{\mathbf S}_i\cdot \hat{\mathbf S}_j +\sum_{ij,\alpha} K_{s}^{\alpha\alpha} \Big(\frac14-\hat{\mathbf S}_i\cdot \hat{\mathbf S}_j\Big) \hat\tau_i^\alpha \hat\tau_j^\alpha +\sum_{ij,\alpha} K_{t}^{\alpha\alpha} \Big(\hat{\mathbf S}_i\cdot \hat{\mathbf S}_j+\frac34\Big) \hat\tau_i^\alpha \hat\tau_j^\alpha,
\]
and its ordered phases are then analyzed by Weiss mean-field theory and generalized flavor-wave theory [2601.06440]. A related derivation from a two-orbital Hubbard model with \(U\), \(U'\), crystal-field splitting, and Hund coupling yields a projected multicomponent \(t\)-\(J\) theory with local spin-\(\tfrac12\) states, spin-1 triplets, and additional orbital-conversion and Kugel-Khomskii-type terms; this was proposed as the appropriate low-energy model for bilayer nickelates [2504.10114].

An alternative strong-coupling route starts directly from a two-orbital bilayer \(t\)-\(J\) model with onsite ferromagnetic Hund’s-rule coupling and intersite interorbital hybridization. In that framework, \(d_{z^2}\) behaves as a localized bilayer-rung subsystem with large \(t_\perp\) and \(J_\perp\), whereas \(d_{x^2-y^2}\) is the more itinerant pairing-active band [2311.12769]. The moiré counterpart combines projection-based Wannierization, unrestricted Hartree-Fock, and Lanczos exact diagonalization to expose how broad moiré Wannier functions invalidate a purely local Kanamori approximation [2303.02305].

Methodological pluralism is therefore not incidental. It reflects the fact that the same nominal model class supports weak-coupling Fermi-surface competition, strong-coupling spin-orbital-layer exchange, and bilayer-to-effective-single-orbital reductions, depending on filling, hopping anisotropy, and the degree of interaction truncation [2312.03605][2601.06440][2504.10114].

## 5. Correlated phases outside superconductivity

At half filling in the bilayer-derived Kanamori-Moiré-Hubbard model, the phase diagram in dielectric constant \(\varepsilon\) and twist angle \(\theta\) contains three principal phases: antiferromagnetism, an \(S=1\) ferromagnetic insulator, and a charge density wave. For \(\theta<1.75^\circ\) and lower screening the dominant phase is the \(S=1\) ferromagnetic insulator, stabilized by competition between nonlocal ferromagnetic direct exchange and antiferromagnetic superexchange; the AFM state appears only for rather large dielectric constant, \(\varepsilon>20\). For \(\theta\ge 1.75^\circ\), lower \(\varepsilon\) instead favors a CDW due primarily to nonlocal density-density repulsion. A simplified interaction retaining only density-density, direct-exchange, and pair-hopping terms reproduces nearly the same three-phase structure, indicating that nonlocal density repulsion and nonlocal direct exchange are the essential half-filled ingredients [2303.02305].

The explicit bilayer Kanamori model in the quarter-hole-filled molecular-orbital regime yields a different but equally rich phase structure. With couplings \(J_s=0.2\), \(K_t^{zz}=1.5\), \(K_t^{xx}=K_t^{yy}=1.2\), and variable \(K_s^{xx}=K_s^{yy}=\Delta_s K_s^{zz}\), Weiss mean-field theory finds four phases: FM-LS, AFM-LS, AFM-LC, and SLE. The LS states exhibit layer-staggered pseudospin order, the AFM-LC phase has spontaneous interlayer coherence with \(\langle \hat\tau^\pm\rangle\neq 0\), and the SLE phase has no ordinary dipolar spin or layer order but instead a composite order parameter \(\hat{\mathcal Q}_i^{\beta\alpha}=\hat S_i^\beta \hat\tau_i^\alpha\). The SLE phase is described as arising from emergent \(O(4)\to O(3)\) symmetry breaking and supports three entangled gapless Goldstone modes [2601.06440].

Normal-state bilayer nickelate models reveal yet another sequence. In zero-temperature CDMFT+NORG, the realistic La327-motivated density-density Kanamori-Hubbard model displays
\[
\text{Mott} \rightarrow \text{PG} \rightarrow \text{NFL} \rightarrow \text{FL}
\]
with increasing hole doping. The nominal La327 point lies at \(p\approx 0.25\), \(\mu=-1.59\), \(n_x=0.59\), \(n_z=0.91\), inside the NFL regime. A central result is that Hund spin correlation transmits strong interlayer AFM correlations from the \(z\) orbital sector into the \(x\) orbital sector; when the bonding \(z+\) band is no longer fully filled, these interlayer AFM correlations weaken rapidly [2407.13737].

At half filling in a sign-problem-free simplified bilayer two-orbital Hubbard model, DQMC reveals a nonmagnetic weakly insulating phase at weak \(U\), followed by a transition to an antiferromagnetic Mott insulator at
\[
U/t_1^x \approx 4.15.
\]
The transition is consistent with the 3D \(O(3)\) Heisenberg universality class. The weakly insulating phase shows a full single-particle gap, a small but finite spin gap, no extrapolated CDW, PDW, excitonic, or bond-order parameter, and pronounced interlayer \(d_{3z^2-r^2}\) correlations. The same calculations show a downward shift and extended flatness of the \(\gamma\) band, mirroring ARPES observations [2408.12042].

These results collectively undercut a simplified expectation that bilayer two-orbital models generically interpolate only between weak itinerant antiferromagnetism and a standard Mott insulator. Nonlocal exchange can stabilize ferromagnetic or CDW phases in moiré realizations; strong interlayer dimerization can generate nonmagnetic weak insulators; and Hund-coupled bilayer geometries can produce composite spin-layer order not visible in conventional dipolar observables [2303.02305][2408.12042][2601.06440].

## 6. Superconductivity, pairing mechanisms, and model reductions

The superconducting literature on bilayer two-orbital Hubbard and Kanamori-Hubbard models is dominated by \(s_{\pm}\)-type interlayer pairing, but the microscopic meaning of that state varies substantially across models. In the cluster-DMFT study of the bilayer two-orbital Hubbard model with \(U_z=7\), superconductivity is diagnosed from the interlayer \(d_{z^2}\) spin-singlet susceptibility
\[
\chi_{z_1 z_2} = \frac{1}{N}\sum_{i,j}\int_0^\beta \left\langle \mathcal{T}_\tau p_j^\dagger(\tau) p_i(0)\right\rangle d\tau,
\qquad
p_i = \frac{1}{\sqrt{2}} \left( d_{i,1,\uparrow} d_{i,2,\downarrow} - d_{i,1,\downarrow} d_{i,2,\uparrow} \right),
\]
and the resulting state is identified as \(s_{\pm}\)-wave. Its \(T_c\) is nonmonotonic in both \(t_\perp\) and \(V\): a threshold \(V_c\sim 0.4\) is required, the optimal regime is around \(V/t\approx 0.55\text{--}0.60\), and for \(t_\perp=1.3\), \(V=0.5\), \(n_z=0.45\), the calculated \(T_c\) becomes almost independent of \(U_z\) for \(U_z\gtrsim 6\), corresponding to about \(100\) K when \(t_x^1=0.48\) eV is used for conversion. The physical interpretation is a two-component mechanism in which \(t_\perp\) provides pairing glue through \(J_\perp\sim t_\perp^2/U_z\), while \(V\) supplies phase coherence via the itinerant \(d_{x^2-y^2}\) sector [2312.03605].

A more explicitly Kanamori-based CDMFT study argues that bilayer \(\mathrm{La_3Ni_2O_7}\) supports two intertwined \(s_{\pm}\)-wave superconductivities. SC I is strongest at \(\delta_z=0\), with the \(d_{z^2}\) bonding band below \(E_F\), and is strongly enhanced by \(J_H\) but only weakly dependent on \(V\); at \(J_H=0.2U\) and \(V=0.5t\), it gives \(T_c\approx 0.017t\), and it survives even at \(V=0\). SC II becomes dominant at finite \(d_{z^2}\) hole doping, is relatively insensitive to \(J_H\), depends critically on hybridization with \(V_c\sim 0.35t\), and peaks near \(\delta_z=0.08\) with \(T_c\approx 0.02t\). Because SC I weakens while SC II strengthens over comparable \(\delta_z\) intervals, the combined \(T_c(\delta_z)\) evolves smoothly even though the underlying pairing mechanism changes [2508.04554].

The realistic DCA study of the two-orbital Hubbard-Hund model for pressurized \(\mathrm{La_3Ni_2O_7}\) finds a leading \(s_{\pm}\) instability with form factor
\[
g_{s^\pm}^{\ell\ell'}({\bm k}) = \cos k_z \,\delta_{\ell\ell'},
\]
so the sign change is between bonding and antibonding bilayer sectors rather than between distinct in-plane pockets. Real-space eigenvector analysis shows that the leading pair is primarily interlayer, intraorbital, local in-plane, and overwhelmingly \(d_{3z^2-r^2}\)-dominated; interorbital pair correlations do not control the leading instability. The same study identifies interlayer antiferromagnetic spin fluctuations in the \(d_{3z^2-r^2}\) orbital as the pairing driver and argues that a single-orbital bilayer Hubbard model for \(d_{3z^2-r^2}\) is therefore an excellent effective low-energy description of the superconducting behavior [2506.07741].

The DMRG study of a DFT-parameterized bilayer two-orbital Hubbard ladder adds a pressure axis. Near \(\eta_e=3/8\), the dominant instability is an SDW with \(Q_S=\pi/2\); at \(\eta_e=21/64\), superconducting correlations appear with interlayer singlet pairing in both \(d_{z^2}\) and \(d_{x^2-y^2}\), plus in-plane \(d_{z^2}\) singlet pairing. At 21.6 GPa the exponents are approximately \(K_{SC}^z\approx 1\) and \(K_{SC}^x\approx 1.5\), while increasing pressure weakens superconductivity and drives an SC-to-Luttinger-liquid crossover near \(80\) GPa. The comparative analysis in that work attributes the suppression primarily to the reduction of \(U/t\) and \(J_H/t\), rather than to large changes in hopping anisotropy ratios [2511.01801].

Auxiliary models clarify what is generic and what is not. In the single-orbital bilayer Hubbard model, DCA finds a crossover from \(B_{1g}\) \(d_{x^2-y^2}\)-wave pairing at small \(t_\perp/t\) to \(A_{1g}\) \(s^\pm\)-wave pairing at larger \(t_\perp/t\), with interlayer spin fluctuations favoring the latter [1107.0401]. In a sign-problem-free bilayer \(U(2)\times U(2)\)-symmetric Hubbard model with an Ising-like ferromagnetic interlayer interaction, the favored superconducting state is not singlet \(s_{\pm}\) but equal-spin interlayer pairing \(\Delta_{\perp\sigma}=c_{1i\sigma}c_{2i\sigma}\), again coexisting with stripe order and forming a dome in doping [2304.13083]. These analogues suggest that bilayer geometry robustly promotes interlayer pairing channels, while the spin structure of the condensate depends sensitively on whether the interorbital coupling is antiferromagnetic, ferromagnetic, SU(2)-invariant, or Ising-like.

A recurrent conceptual outcome is model reduction. Realistic two-orbital bilayer models often flow toward smaller effective descriptions: to a \(d_{3z^2-r^2}\)-only bilayer Hubbard model in DCA [2506.07741], to a \(d_{x^2-y^2}\)-based bilayer \(t\)-\(J\)-\(J_\perp\) model after Hund-mediated transfer of interlayer AF correlations from a localized \(d_{z^2}\) sector [2311.12769], or to a projected \(t\)-\(J\) theory with local spin-\(\tfrac12\) and spin-1 states once the on-site two-orbital multiplets are treated exactly in strong coupling [2504.10114]. This suggests that the encyclopedic subject is best understood not as a single Hamiltonian but as a hierarchy of related models whose common core is the competition between bilayer bonding, orbital differentiation, and multiorbital local interactions.

Source: https://www.emergentmind.com/topics/bilayer-two-orbital-kanamori-hubbard-model