thanks: These two authors contributed equallythanks: These two authors contributed equally

Type II t-J model in charge transfer regime in bilayer La3Ni2O7 and trilayer La4Ni3O10

Hanbit Oh [email protected]    Boran Zhou    Ya-Hui Zhang [email protected] William H. Miller III Department of Physics and Astronomy, Johns Hopkins University, Baltimore, Maryland, 21218, USA
Abstract

Recent observations of an 80 K superconductor in La3Ni2O7 under high pressure have attracted significant attention. Recent experiments indicate that La3Ni2O7 may be in the charge transfer regime, challenging the previous models based purely on the Ni dx2y2subscript𝑑superscript𝑥2superscript𝑦2d_{x^{2}-y^{2}}italic_d start_POSTSUBSCRIPT italic_x start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_y start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_POSTSUBSCRIPT and dz2subscript𝑑superscript𝑧2d_{z^{2}}italic_d start_POSTSUBSCRIPT italic_z start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_POSTSUBSCRIPT orbitals. In this study, we propose a low energy model that incorporates doped holes in the oxygen p𝑝pitalic_p orbitals. Given that the parent nickel state is in the 3d83superscript𝑑83d^{8}3 italic_d start_POSTSUPERSCRIPT 8 end_POSTSUPERSCRIPT configuration with a spin-one moment, doped hole only screens it down to spin-half, in contrast to the Zhang-Rice singlet in cuprate. We dub the single hole state as Zhang-Rice spin-half and build an effective model which includes three spin-one states (d8superscript𝑑8d^{8}italic_d start_POSTSUPERSCRIPT 8 end_POSTSUPERSCRIPT) and two Zhang-Rice spin-half states (d8Lsuperscript𝑑8𝐿d^{8}Litalic_d start_POSTSUPERSCRIPT 8 end_POSTSUPERSCRIPT italic_L). At moderate pressure around 20202020 GPa, the dominated oxygen orbital is an in-plane Wannier orbital with the same lattice symmetry as the dx2y2subscript𝑑superscript𝑥2superscript𝑦2d_{x^{2}-y^{2}}italic_d start_POSTSUBSCRIPT italic_x start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_y start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_POSTSUBSCRIPT orbital. The resulting model reduces to the bilayer type II t-J model previously proposed in the Mott-Hubbard regime. Notably, the hopping between the in-plane p𝑝pitalic_p orbitals of the two layers is still suppressed. Density matrix renormalization group (DMRG) simulation reveals a pairing dome with the optimal hole doping level at x=0.40.5𝑥0.4similar-to0.5x=0.4\sim 0.5italic_x = 0.4 ∼ 0.5, distinct from the hole doped cuprate where optimal doping occurs around x=0.19𝑥0.19x=0.19italic_x = 0.19. Further increasing pressure initially raises the critical temperature (Tcsubscript𝑇𝑐T_{c}italic_T start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT) until reaching an optimal pressure beyond which the pzsubscript𝑝𝑧p_{z}italic_p start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT orbital of oxygen becomes favorable and superconductivity is diminished. This shift from in-plane p𝑝pitalic_p orbital to pzsubscript𝑝𝑧p_{z}italic_p start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT orbital may elucidate the experimentally observed superconducting dome with varying pressure. As an extension, we also suggest a trilayer version of the type II t-J model as the minimal model for pressured La4Ni3O10, which is distinct from the models in the Mott-Hubbard regime.

Introduction.— The recent discovery of a superconductor with a critical temperature of approximately 80 K in La3Ni2O7 under high pressureSun et al. (2023); Hou et al. (2023); Zhang et al. (2023a) has sparked considerable interestLiu et al. (2023a); Yang et al. (2023a); Zhang et al. (2023b); Luo et al. (2023); Zhang et al. (2023c); Yang et al. (2023b); Sakakibara et al. (2023); Gu et al. (2023); Shen et al. (2023); Wú et al. (2023); Christiansson et al. (2023); Liu et al. (2023b); Cao and Yang (2023); Oh and Zhang (2023); Lu et al. (2024a); Qu et al. (2023); Lu et al. (2023); Jiang et al. (2023a); Tian et al. (2023); Zhang et al. (2023d); Qin and Yang (2023); Huang et al. (2023); Zhang et al. (2023e); Jiang et al. (2023b); Yang et al. (2023c); Qin and Yang (2023); Zhang et al. (2023); Kitamine et al. (2023); Jiang et al. (2023); Lange et al. (2023a, b); Schlömer et al. (2023); Jiang et al. (2024a); Zhan et al. (2024); Lechermann et al. (2024); Wang et al. (2024); Xue and Wang (2024); Kaneko et al. (2024); Jiang et al. (2024b); Ryee et al. (2023); Liao et al. (2023); Wú et al. (2023). Additionally, there is emerging evidence of superconductivity in La4Ni3O10 under high pressure, exhibiting critical temperatures ranging from 20-30 K, which has further fueled research interest in this area Li et al. (2024a); Zhu et al. (2024); Sakakibara et al. (2024); Zhang et al. (2024a); Lu et al. (2024b); Zhang et al. (2024b); LaBollita et al. (2024); Yang et al. (2024). Identifying a minimal model that captures the essential physics is a crucial step forward.

According to density functional theory (DFT) Sun et al. (2023), the valence of the nickel (Ni) atom in the bilayer nickelate La3Ni2O7 is in the 3d8x3superscript𝑑8𝑥3d^{8-x}3 italic_d start_POSTSUPERSCRIPT 8 - italic_x end_POSTSUPERSCRIPT configuration. The dx2y2subscript𝑑superscript𝑥2superscript𝑦2d_{x^{2}-y^{2}}italic_d start_POSTSUBSCRIPT italic_x start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_y start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_POSTSUBSCRIPT orbital is at density n1=1xsubscript𝑛11𝑥n_{1}=1-xitalic_n start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT = 1 - italic_x per site with x=0.5𝑥0.5x=0.5italic_x = 0.5, while the dz2subscript𝑑superscript𝑧2d_{z^{2}}italic_d start_POSTSUBSCRIPT italic_z start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_POSTSUBSCRIPT orbital has the density n21subscript𝑛21n_{2}\approx 1italic_n start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ≈ 1 and is near Mott localization. Following this picture, many theoretical works propose a two-orbital Hubbard or t-J model in terms of the dx2y2subscript𝑑superscript𝑥2superscript𝑦2d_{x^{2}-y^{2}}italic_d start_POSTSUBSCRIPT italic_x start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_y start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_POSTSUBSCRIPT and the dz2subscript𝑑superscript𝑧2d_{z^{2}}italic_d start_POSTSUBSCRIPT italic_z start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_POSTSUBSCRIPT orbital. Especially, Refs.Oh and Zhang (2023); Lu et al. (2024a) highlight the role of Hund’s coupling JHsubscript𝐽𝐻J_{H}italic_J start_POSTSUBSCRIPT italic_H end_POSTSUBSCRIPT in transmitting a large inter-layer spin-spin coupling Jsubscript𝐽perpendicular-toJ_{\perp}italic_J start_POSTSUBSCRIPT ⟂ end_POSTSUBSCRIPT from the dz2subscript𝑑superscript𝑧2d_{z^{2}}italic_d start_POSTSUBSCRIPT italic_z start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_POSTSUBSCRIPT to the dx2y2subscript𝑑superscript𝑥2superscript𝑦2d_{x^{2}-y^{2}}italic_d start_POSTSUBSCRIPT italic_x start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_y start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_POSTSUBSCRIPT orbital, providing a plausible explanation for a high critical temperature (Tcsubscript𝑇𝑐T_{c}italic_T start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT) even at large doping x0.5𝑥0.5x\approx 0.5italic_x ≈ 0.5. In Ref.Lu et al. (2024a), the localized dz2subscript𝑑superscript𝑧2d_{z^{2}}italic_d start_POSTSUBSCRIPT italic_z start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_POSTSUBSCRIPT orbitals are simply integrated out, so the model is reduced to a one-orbital bilayer t-J model of dx2y2subscript𝑑superscript𝑥2superscript𝑦2d_{x^{2}-y^{2}}italic_d start_POSTSUBSCRIPT italic_x start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_y start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_POSTSUBSCRIPT with a negligible tsubscript𝑡perpendicular-tot_{\perp}italic_t start_POSTSUBSCRIPT ⟂ end_POSTSUBSCRIPT. But the integration of dz2subscript𝑑superscript𝑧2d_{z^{2}}italic_d start_POSTSUBSCRIPT italic_z start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_POSTSUBSCRIPT orbital is not appropriate in the large JHsubscript𝐽𝐻J_{H}italic_J start_POSTSUBSCRIPT italic_H end_POSTSUBSCRIPT limit Oh and Zhang (2023); Yang et al. (2023d). Instead, the minimal model is demonstrated to be a bilayer type-II t-J model Oh and Zhang (2023); Yang et al. (2023d), which includes the spin-one Ni2+ state (3d8superscript𝑑8d^{8}italic_d start_POSTSUPERSCRIPT 8 end_POSTSUPERSCRIPT) and spin-half Ni3+ state (3d7superscript𝑑7d^{7}italic_d start_POSTSUPERSCRIPT 7 end_POSTSUPERSCRIPT) at each site. The type II t-J model hosts unique physics such as a dome of pairing gap around x0.5𝑥0.5x\approx 0.5italic_x ≈ 0.5 due to doping induced Bardeen-Cooper-Schrieffer(BCS) to Bose-Einstein-condensate (BEC) crossoverYang et al. (2023d). There are also two different Fermi liquids above Tcsubscript𝑇𝑐T_{c}italic_T start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT with a jump of the Fermi surface volume by 1/2121/21 / 2 Brillouin zoneOh and Zhang (2023); Yang et al. (2023d); Wu et al. (2024). Clearly the physics is essentially different from the familiar hole doped cuprates. However, the model is derivedOh and Zhang (2023) assuming that the system is in the Mott-Hubbard regime of the Zaanen–Sawatzky–Allen classification scheme Zaanen et al. (1985). Depending on the energy splitting ΔΔ\Deltaroman_Δ between the oxygen p𝑝pitalic_p orbital and the Ni d𝑑ditalic_d orbital compared to Hubbard U, we have the Mott-Hubbard regime (Δ>UΔ𝑈\Delta>Uroman_Δ > italic_U) or the charge-transfer regime (Δ<UΔ𝑈\Delta<Uroman_Δ < italic_U). In the Mott-Hubbard regime, the doped holes enter the 3d orbitals of the Ni atom and we can ignore the oxygens. However, a recent experiment suggests that La3Ni2O7 might be within the charge transfer regime where holes enter the oxygen p𝑝pitalic_p orbitals while the nickel atom is pinned to be in the 3d8superscript𝑑8d^{8}italic_d start_POSTSUPERSCRIPT 8 end_POSTSUPERSCRIPT configuration with a localized spin-one moment Dong et al. (2023). This challenges the previous theoretical models. It is thus critical to derive an effective model for the charge transfer regime, akin to what has been done in hole doped cuprates Zhang and Rice (1988).

There is a notable distinction between bilayer nickelates and cuprates. The undoped Ni state is in the spin-one 3d83superscript𝑑83d^{8}3 italic_d start_POSTSUPERSCRIPT 8 end_POSTSUPERSCRIPT configuration with two electrons in the two egsubscript𝑒𝑔e_{g}italic_e start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT orbitals, (dx2y2,dz2)subscript𝑑superscript𝑥2superscript𝑦2subscript𝑑superscript𝑧2(d_{x^{2}-y^{2}},d_{z^{2}})( italic_d start_POSTSUBSCRIPT italic_x start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_y start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_POSTSUBSCRIPT , italic_d start_POSTSUBSCRIPT italic_z start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_POSTSUBSCRIPT ), strongly coupled together by a Hund coupling JHsubscript𝐽𝐻J_{H}italic_J start_POSTSUBSCRIPT italic_H end_POSTSUBSCRIPT. Thus doped hole in the oxygen can at best screen the spin-one moment down to spin-half, contrasting the Zhang-Rice singlet in cuprate. In this work, we demonstrate that the doped hole enters the in-plane oxygen orbital and forms a net spin-half moment together with the Ni spin-one moment at moderate pressure. The hole state state is dubbed as the Zhang-Rice spin-half. Then by taking the Zhang-Rice spin-half as the primary state of the single-hole d8Lsuperscript𝑑8𝐿d^{8}Litalic_d start_POSTSUPERSCRIPT 8 end_POSTSUPERSCRIPT italic_L state and keeping the spin-triplet parent state of the undoped d8superscript𝑑8d^{8}italic_d start_POSTSUPERSCRIPT 8 end_POSTSUPERSCRIPT state, the type-II t-J model Zhang and Vishwanath (2020) is shown to be the minimal model also in the charge-transfer regime of the doped bilayer nickelates. Besides, we reveal that the oxygen pzsubscript𝑝𝑧p_{z}italic_p start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT orbital living between the two layers can only mediate a very small inter-layer hopping tsubscript𝑡perpendicular-tot_{\perp}italic_t start_POSTSUBSCRIPT ⟂ end_POSTSUBSCRIPT due to symmetry. From the density matrix renormalization group (DMRG) simulation, we show that the pairing gradually decreases with tsubscript𝑡perpendicular-tot_{\perp}italic_t start_POSTSUBSCRIPT ⟂ end_POSTSUBSCRIPT but remains robust. If we further increase across an optimal value, the doped hole shifts from the intra-layer px,pysubscript𝑝𝑥subscript𝑝𝑦p_{x},p_{y}italic_p start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT , italic_p start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT oxygen orbital to the inter-layer pzsubscript𝑝𝑧p_{z}italic_p start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT oxygen and the pairing is suppressed. We finally point out that the trilayer La4Ni3O10 is also described by a trilayer version of type II t-J model in the charge transfer regime.

Charge transfer model.— We consider a model on bilayer square lattice (See Fig. 1), which includes three p𝑝pitalic_p-orbitals (2px,2py,2pz2subscript𝑝𝑥2subscript𝑝𝑦2subscript𝑝𝑧2p_{x},2p_{y},2p_{z}2 italic_p start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT , 2 italic_p start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT , 2 italic_p start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT) of oxygen as well as the two Egsubscript𝐸𝑔E_{g}italic_E start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT d𝑑ditalic_d-orbitals (3dx2y2,3dz23subscript𝑑superscript𝑥2superscript𝑦23subscript𝑑superscript𝑧23d_{x^{2}-y^{2}},3d_{z^{2}}3 italic_d start_POSTSUBSCRIPT italic_x start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_y start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_POSTSUBSCRIPT , 3 italic_d start_POSTSUBSCRIPT italic_z start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_POSTSUBSCRIPT) of nickel. In the hole picture, we first focus on one-unit cell with two Ni atoms and nine O atoms, the Hamiltonian is,

H𝐻\displaystyle Hitalic_H =\displaystyle== Hdp+i,l[U1n1;i;l;dn1;i;l;d+U2n2;i;l;dn2;i;l;d\displaystyle H_{dp}+\sum_{i,l}\left[U_{1}n_{1;i;l;\uparrow}^{d}n^{d}_{1;i;l;% \downarrow}+U_{2}n_{2;i;l;\uparrow}^{d}n_{2;i;l;\downarrow}^{d}\right.italic_H start_POSTSUBSCRIPT italic_d italic_p end_POSTSUBSCRIPT + ∑ start_POSTSUBSCRIPT italic_i , italic_l end_POSTSUBSCRIPT [ italic_U start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT italic_n start_POSTSUBSCRIPT 1 ; italic_i ; italic_l ; ↑ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_d end_POSTSUPERSCRIPT italic_n start_POSTSUPERSCRIPT italic_d end_POSTSUPERSCRIPT start_POSTSUBSCRIPT 1 ; italic_i ; italic_l ; ↓ end_POSTSUBSCRIPT + italic_U start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT italic_n start_POSTSUBSCRIPT 2 ; italic_i ; italic_l ; ↑ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_d end_POSTSUPERSCRIPT italic_n start_POSTSUBSCRIPT 2 ; italic_i ; italic_l ; ↓ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_d end_POSTSUPERSCRIPT (1)
+Un1;i;ldn2;i;ld2JH(s1;i;lds2;i;ld+14n1;i;ldn2;i;ld),superscript𝑈subscriptsuperscript𝑛𝑑1𝑖𝑙subscriptsuperscript𝑛𝑑2𝑖𝑙2subscript𝐽𝐻subscriptsuperscript𝑠𝑑1𝑖𝑙subscriptsuperscript𝑠𝑑2𝑖𝑙14subscriptsuperscript𝑛𝑑1𝑖𝑙subscriptsuperscript𝑛𝑑2𝑖𝑙\displaystyle+U^{\prime}n^{d}_{1;i;l}n^{d}_{2;i;l}-2J_{H}(\vec{s}^{d}_{1;i;l}% \cdot\vec{s}^{d}_{2;i;l}+\frac{1}{4}n^{d}_{1;i;l}n^{d}_{2;i;l}),+ italic_U start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT italic_n start_POSTSUPERSCRIPT italic_d end_POSTSUPERSCRIPT start_POSTSUBSCRIPT 1 ; italic_i ; italic_l end_POSTSUBSCRIPT italic_n start_POSTSUPERSCRIPT italic_d end_POSTSUPERSCRIPT start_POSTSUBSCRIPT 2 ; italic_i ; italic_l end_POSTSUBSCRIPT - 2 italic_J start_POSTSUBSCRIPT italic_H end_POSTSUBSCRIPT ( over→ start_ARG italic_s end_ARG start_POSTSUPERSCRIPT italic_d end_POSTSUPERSCRIPT start_POSTSUBSCRIPT 1 ; italic_i ; italic_l end_POSTSUBSCRIPT ⋅ over→ start_ARG italic_s end_ARG start_POSTSUPERSCRIPT italic_d end_POSTSUPERSCRIPT start_POSTSUBSCRIPT 2 ; italic_i ; italic_l end_POSTSUBSCRIPT + divide start_ARG 1 end_ARG start_ARG 4 end_ARG italic_n start_POSTSUPERSCRIPT italic_d end_POSTSUPERSCRIPT start_POSTSUBSCRIPT 1 ; italic_i ; italic_l end_POSTSUBSCRIPT italic_n start_POSTSUPERSCRIPT italic_d end_POSTSUPERSCRIPT start_POSTSUBSCRIPT 2 ; italic_i ; italic_l end_POSTSUBSCRIPT ) ,
+aϵd,ana;i,ld]+i,l,αϵpnα;i;lp+iϵp;znpz;i,\displaystyle\left.+\sum_{a}\epsilon_{d,a}n_{a;i,l}^{d}\right]+\sum_{i^{\prime% },l,\alpha}\epsilon_{p}n_{\alpha;i^{\prime};l}^{p}+\sum_{i}\epsilon_{p;z}n^{p}% _{z;i},+ ∑ start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT italic_ϵ start_POSTSUBSCRIPT italic_d , italic_a end_POSTSUBSCRIPT italic_n start_POSTSUBSCRIPT italic_a ; italic_i , italic_l end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_d end_POSTSUPERSCRIPT ] + ∑ start_POSTSUBSCRIPT italic_i start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT , italic_l , italic_α end_POSTSUBSCRIPT italic_ϵ start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT italic_n start_POSTSUBSCRIPT italic_α ; italic_i start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ; italic_l end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_p end_POSTSUPERSCRIPT + ∑ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT italic_ϵ start_POSTSUBSCRIPT italic_p ; italic_z end_POSTSUBSCRIPT italic_n start_POSTSUPERSCRIPT italic_p end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_z ; italic_i end_POSTSUBSCRIPT ,

and

Hdpsubscript𝐻𝑑𝑝\displaystyle H_{dp}italic_H start_POSTSUBSCRIPT italic_d italic_p end_POSTSUBSCRIPT =\displaystyle== i,l,a,σ[2tdp;ada;i;l;σpa;i;l;σ+H.c.]\displaystyle\sum_{i,l,a,\sigma}\left[2t_{dp;a}d^{\dagger}_{a;i;l;\sigma}p_{a;% i;l;\sigma}+\mathrm{H.c.}\right]∑ start_POSTSUBSCRIPT italic_i , italic_l , italic_a , italic_σ end_POSTSUBSCRIPT [ 2 italic_t start_POSTSUBSCRIPT italic_d italic_p ; italic_a end_POSTSUBSCRIPT italic_d start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_a ; italic_i ; italic_l ; italic_σ end_POSTSUBSCRIPT italic_p start_POSTSUBSCRIPT italic_a ; italic_i ; italic_l ; italic_σ end_POSTSUBSCRIPT + roman_H . roman_c . ]
+\displaystyle++ i,σ[tdp;z(d2;i;t;σd2;i;b;σ)pz;i;σ+H.c.]\displaystyle\sum_{i,\sigma}\left[t_{dp;z}(d^{\dagger}_{2;i;t;\sigma}-d^{% \dagger}_{2;i;b;\sigma})p_{z;i;\sigma}+\mathrm{H.c.}\right]∑ start_POSTSUBSCRIPT italic_i , italic_σ end_POSTSUBSCRIPT [ italic_t start_POSTSUBSCRIPT italic_d italic_p ; italic_z end_POSTSUBSCRIPT ( italic_d start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT start_POSTSUBSCRIPT 2 ; italic_i ; italic_t ; italic_σ end_POSTSUBSCRIPT - italic_d start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT start_POSTSUBSCRIPT 2 ; italic_i ; italic_b ; italic_σ end_POSTSUBSCRIPT ) italic_p start_POSTSUBSCRIPT italic_z ; italic_i ; italic_σ end_POSTSUBSCRIPT + roman_H . roman_c . ]

with

p1;i;l;σ=12[px;i+x^2;l;σpx;ix^2;l;σpy;i+y^2;l;σ+py;iy^2;l;σ],subscript𝑝1𝑖𝑙𝜎12delimited-[]subscript𝑝𝑥𝑖^𝑥2𝑙𝜎subscript𝑝𝑥𝑖^𝑥2𝑙𝜎subscript𝑝𝑦𝑖^𝑦2𝑙𝜎subscript𝑝𝑦𝑖^𝑦2𝑙𝜎\displaystyle p_{1;i;l;\sigma}=\frac{1}{2}\left[p_{x;i+\frac{{\widehat{x}}}{2}% ;l;\sigma}-p_{x;i-\frac{{\widehat{x}}}{2};l;\sigma}-p_{y;i+\frac{{\widehat{y}}% }{2};l;\sigma}+p_{y;i-\frac{{\widehat{y}}}{2};l;\sigma}\right],italic_p start_POSTSUBSCRIPT 1 ; italic_i ; italic_l ; italic_σ end_POSTSUBSCRIPT = divide start_ARG 1 end_ARG start_ARG 2 end_ARG [ italic_p start_POSTSUBSCRIPT italic_x ; italic_i + divide start_ARG over^ start_ARG italic_x end_ARG end_ARG start_ARG 2 end_ARG ; italic_l ; italic_σ end_POSTSUBSCRIPT - italic_p start_POSTSUBSCRIPT italic_x ; italic_i - divide start_ARG over^ start_ARG italic_x end_ARG end_ARG start_ARG 2 end_ARG ; italic_l ; italic_σ end_POSTSUBSCRIPT - italic_p start_POSTSUBSCRIPT italic_y ; italic_i + divide start_ARG over^ start_ARG italic_y end_ARG end_ARG start_ARG 2 end_ARG ; italic_l ; italic_σ end_POSTSUBSCRIPT + italic_p start_POSTSUBSCRIPT italic_y ; italic_i - divide start_ARG over^ start_ARG italic_y end_ARG end_ARG start_ARG 2 end_ARG ; italic_l ; italic_σ end_POSTSUBSCRIPT ] ,
p2;i;l;σ=12[px;i+x^2;l;σpx;ix^2;l;σ+py;i+y^2;l;σpy;iy^2;l;σ],subscript𝑝2𝑖𝑙𝜎12delimited-[]subscript𝑝𝑥𝑖^𝑥2𝑙𝜎subscript𝑝𝑥𝑖^𝑥2𝑙𝜎subscript𝑝𝑦𝑖^𝑦2𝑙𝜎subscript𝑝𝑦𝑖^𝑦2𝑙𝜎\displaystyle p_{2;i;l;\sigma}=\frac{1}{2}\left[p_{x;i+\frac{{\widehat{x}}}{2}% ;l;\sigma}-p_{x;i-\frac{{\widehat{x}}}{2};l;\sigma}+p_{y;i+\frac{{\widehat{y}}% }{2};l;\sigma}-p_{y;i-\frac{{\widehat{y}}}{2};l;\sigma}\right],italic_p start_POSTSUBSCRIPT 2 ; italic_i ; italic_l ; italic_σ end_POSTSUBSCRIPT = divide start_ARG 1 end_ARG start_ARG 2 end_ARG [ italic_p start_POSTSUBSCRIPT italic_x ; italic_i + divide start_ARG over^ start_ARG italic_x end_ARG end_ARG start_ARG 2 end_ARG ; italic_l ; italic_σ end_POSTSUBSCRIPT - italic_p start_POSTSUBSCRIPT italic_x ; italic_i - divide start_ARG over^ start_ARG italic_x end_ARG end_ARG start_ARG 2 end_ARG ; italic_l ; italic_σ end_POSTSUBSCRIPT + italic_p start_POSTSUBSCRIPT italic_y ; italic_i + divide start_ARG over^ start_ARG italic_y end_ARG end_ARG start_ARG 2 end_ARG ; italic_l ; italic_σ end_POSTSUBSCRIPT - italic_p start_POSTSUBSCRIPT italic_y ; italic_i - divide start_ARG over^ start_ARG italic_y end_ARG end_ARG start_ARG 2 end_ARG ; italic_l ; italic_σ end_POSTSUBSCRIPT ] ,

where a=1,2𝑎12a=1,2italic_a = 1 , 2 labels dx2y2subscript𝑑superscript𝑥2superscript𝑦2d_{x^{2}-y^{2}}italic_d start_POSTSUBSCRIPT italic_x start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_y start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_POSTSUBSCRIPT and dz2subscript𝑑superscript𝑧2d_{z^{2}}italic_d start_POSTSUBSCRIPT italic_z start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_POSTSUBSCRIPT orbital of Ni and l=t,b𝑙𝑡𝑏l=t,bitalic_l = italic_t , italic_b labels the top and bottom layer. σ=,𝜎\sigma=\uparrow,\downarrowitalic_σ = ↑ , ↓ labels the spin. Here da;i;l;σsubscriptsuperscript𝑑𝑎𝑖𝑙𝜎d^{\dagger}_{a;i;l;\sigma}italic_d start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_a ; italic_i ; italic_l ; italic_σ end_POSTSUBSCRIPT creates a hole on the nickel site i𝑖iitalic_i relative to the 3d103superscript𝑑103d^{10}3 italic_d start_POSTSUPERSCRIPT 10 end_POSTSUPERSCRIPT state. pa;i;l;σsubscriptsuperscript𝑝𝑎𝑖𝑙𝜎p^{\dagger}_{a;i;l;\sigma}italic_p start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_a ; italic_i ; italic_l ; italic_σ end_POSTSUBSCRIPT creates a hole in a superposition of occupying the four oxygen atom sites isuperscript𝑖i^{\prime}italic_i start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT around the Ni site i𝑖iitalic_i. The px,pysubscript𝑝𝑥subscript𝑝𝑦p_{x},p_{y}italic_p start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT , italic_p start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT orbitals are regrouped as p1,p2subscript𝑝1subscript𝑝2p_{1},p_{2}italic_p start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , italic_p start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT orbitals living on the Ni site according to the D4hsubscript𝐷4D_{4h}italic_D start_POSTSUBSCRIPT 4 italic_h end_POSTSUBSCRIPT symmetry. One can check that p1,p2subscript𝑝1subscript𝑝2p_{1},p_{2}italic_p start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , italic_p start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT have the same lattice symmetry as the d1,d2subscript𝑑1subscript𝑑2d_{1},d_{2}italic_d start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , italic_d start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT orbitals. In supplemental material (SM) Sec.I, we provide the details on the symmetry analysis. U𝑈Uitalic_U(Usuperscript𝑈U^{\prime}italic_U start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT) is the intra-(inter-) orbital onsite repulsion, and JH>0subscript𝐽𝐻0J_{H}>0italic_J start_POSTSUBSCRIPT italic_H end_POSTSUBSCRIPT > 0 is the Hund coupling. tdpsubscript𝑡𝑑𝑝t_{dp}italic_t start_POSTSUBSCRIPT italic_d italic_p end_POSTSUBSCRIPT is the nearest-neighbor hopping between dp𝑑𝑝d-pitalic_d - italic_p orbitals. There are three different hopping channels, (d1;l,p1;l)subscript𝑑1𝑙subscript𝑝1𝑙(d_{1;l},p_{1;l})( italic_d start_POSTSUBSCRIPT 1 ; italic_l end_POSTSUBSCRIPT , italic_p start_POSTSUBSCRIPT 1 ; italic_l end_POSTSUBSCRIPT ), (d2;l,p2;l)subscript𝑑2𝑙subscript𝑝2𝑙(d_{2;l},p_{2;l})( italic_d start_POSTSUBSCRIPT 2 ; italic_l end_POSTSUBSCRIPT , italic_p start_POSTSUBSCRIPT 2 ; italic_l end_POSTSUBSCRIPT ), and (d2;td2;b,pz)subscript𝑑2𝑡subscript𝑑2𝑏subscript𝑝𝑧(d_{2;t}-d_{2;b},p_{z})( italic_d start_POSTSUBSCRIPT 2 ; italic_t end_POSTSUBSCRIPT - italic_d start_POSTSUBSCRIPT 2 ; italic_b end_POSTSUBSCRIPT , italic_p start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT ) classified by D4hsubscript𝐷4D_{4h}italic_D start_POSTSUBSCRIPT 4 italic_h end_POSTSUBSCRIPT point group, respectively. The relative size of the onsite energies ϵdsubscriptitalic-ϵ𝑑\epsilon_{d}italic_ϵ start_POSTSUBSCRIPT italic_d end_POSTSUBSCRIPT and ϵpsubscriptitalic-ϵ𝑝\epsilon_{p}italic_ϵ start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT decides whether the system is in the charge transfer or Mott Hubbard regime. We consider the case that the undoped parent Ni state is in 3d83superscript𝑑83d^{8}3 italic_d start_POSTSUPERSCRIPT 8 end_POSTSUPERSCRIPT with two holes in the two egsubscript𝑒𝑔e_{g}italic_e start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT orbitals, forming a S=1𝑆1S=1italic_S = 1 localized moment. Then under further hole doping, additional holes enter the oxygen p𝑝pitalic_p orbitals. We note that the orbital p1subscript𝑝1p_{1}italic_p start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT and p2,pzsubscript𝑝2subscript𝑝𝑧p_{2},p_{z}italic_p start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT , italic_p start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT have different eigenvalues under the C4subscript𝐶4C_{4}italic_C start_POSTSUBSCRIPT 4 end_POSTSUBSCRIPT rotation and can not hybridize. So the doped hole enters either p1subscript𝑝1p_{1}italic_p start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT or p2,pzsubscript𝑝2subscript𝑝𝑧p_{2},p_{z}italic_p start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT , italic_p start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT depending on energetics.

Refer to caption
Figure 1: The lattice structure and p,d𝑝𝑑p,ditalic_p , italic_d orbitals of the bilayer Nickelates La3Ni2O7. The green (yellow) sphere denotes Ni (O) atom, respectively. At each Ni atom, there are two d𝑑ditalic_d orbitals, d1,d2subscript𝑑1subscript𝑑2d_{1},d_{2}italic_d start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , italic_d start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT (Left inset). Around each Ni atom at site i𝑖iitalic_i, the px,pysubscript𝑝𝑥subscript𝑝𝑦p_{x},p_{y}italic_p start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT , italic_p start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT orbitals from the four oxygen atoms around the Ni site in the same plane form two orbitals pi;1,pi;2subscript𝑝𝑖1subscript𝑝𝑖2p_{i;1},p_{i;2}italic_p start_POSTSUBSCRIPT italic_i ; 1 end_POSTSUBSCRIPT , italic_p start_POSTSUBSCRIPT italic_i ; 2 end_POSTSUBSCRIPT living on the same site i𝑖iitalic_i. Meanwhile, there is a pzsubscript𝑝𝑧p_{z}italic_p start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT orbital from the oxygen at the center of the z^^𝑧{\widehat{z}}over^ start_ARG italic_z end_ARG bond. At moderate pressure, the doped hole enters the Wannier orbital, p1;i=12[px;i+x^/2px;ix^/2py;i+y^/2+py;iy^/2]subscript𝑝1𝑖12delimited-[]subscript𝑝𝑥𝑖^𝑥2subscript𝑝𝑥𝑖^𝑥2subscript𝑝𝑦𝑖^𝑦2subscript𝑝𝑦𝑖^𝑦2p_{1;i}=\frac{1}{2}[p_{x;i+{\widehat{x}}/{2}}-p_{x;i-{\widehat{x}}/2}-p_{y;i+{% \widehat{y}}/2}+p_{y;i-{\widehat{y}}/2}]italic_p start_POSTSUBSCRIPT 1 ; italic_i end_POSTSUBSCRIPT = divide start_ARG 1 end_ARG start_ARG 2 end_ARG [ italic_p start_POSTSUBSCRIPT italic_x ; italic_i + over^ start_ARG italic_x end_ARG / 2 end_POSTSUBSCRIPT - italic_p start_POSTSUBSCRIPT italic_x ; italic_i - over^ start_ARG italic_x end_ARG / 2 end_POSTSUBSCRIPT - italic_p start_POSTSUBSCRIPT italic_y ; italic_i + over^ start_ARG italic_y end_ARG / 2 end_POSTSUBSCRIPT + italic_p start_POSTSUBSCRIPT italic_y ; italic_i - over^ start_ARG italic_y end_ARG / 2 end_POSTSUBSCRIPT ] which has the same symmetry as the d1subscript𝑑1d_{1}italic_d start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT orbital. The hole in this pi;1subscript𝑝𝑖1p_{i;1}italic_p start_POSTSUBSCRIPT italic_i ; 1 end_POSTSUBSCRIPT orbital couples to the S=1𝑆1S=1italic_S = 1 moment of the Ni 3d83superscript𝑑83d^{8}3 italic_d start_POSTSUPERSCRIPT 8 end_POSTSUPERSCRIPT state and forms a Zhang-Rice spin-half state.

Zhang-Rice Spin-half.— In the strong coupling limit, tdpU,Umuch-less-thansubscript𝑡𝑑𝑝𝑈superscript𝑈t_{dp}\ll U,U^{\prime}italic_t start_POSTSUBSCRIPT italic_d italic_p end_POSTSUBSCRIPT ≪ italic_U , italic_U start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT, the hole enters the oxygen p𝑝pitalic_p orbitals and interact with the localized spin-one moments from the Ni2+ state through a Kondo coupling,

HKsubscript𝐻𝐾\displaystyle H_{K}italic_H start_POSTSUBSCRIPT italic_K end_POSTSUBSCRIPT =\displaystyle== i,l,aJK;a[(pa;i;lσpa;i;l)sa;i;ld12na;i;lp]subscript𝑖𝑙𝑎subscript𝐽𝐾𝑎delimited-[]superscriptsubscript𝑝𝑎𝑖𝑙𝜎subscript𝑝𝑎𝑖𝑙subscriptsuperscript𝑠𝑑𝑎𝑖𝑙12subscriptsuperscript𝑛𝑝𝑎𝑖𝑙\displaystyle\sum_{i,l,a}J_{K;a}\big{[}(p_{a;i;l}^{\dagger}\vec{\sigma}p_{a;i;% l})\cdot\vec{s}^{d}_{a;i;l}-\frac{1}{2}n^{p}_{a;i;l}\big{]}∑ start_POSTSUBSCRIPT italic_i , italic_l , italic_a end_POSTSUBSCRIPT italic_J start_POSTSUBSCRIPT italic_K ; italic_a end_POSTSUBSCRIPT [ ( italic_p start_POSTSUBSCRIPT italic_a ; italic_i ; italic_l end_POSTSUBSCRIPT start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT over→ start_ARG italic_σ end_ARG italic_p start_POSTSUBSCRIPT italic_a ; italic_i ; italic_l end_POSTSUBSCRIPT ) ⋅ over→ start_ARG italic_s end_ARG start_POSTSUPERSCRIPT italic_d end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_a ; italic_i ; italic_l end_POSTSUBSCRIPT - divide start_ARG 1 end_ARG start_ARG 2 end_ARG italic_n start_POSTSUPERSCRIPT italic_p end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_a ; italic_i ; italic_l end_POSTSUBSCRIPT ]
+\displaystyle++ i,lJK;z[(pz;iσpz;i)s2;i;ld12nz;ip]subscript𝑖𝑙subscript𝐽𝐾𝑧delimited-[]superscriptsubscript𝑝𝑧𝑖𝜎subscript𝑝𝑧𝑖subscriptsuperscript𝑠𝑑2𝑖𝑙12subscriptsuperscript𝑛𝑝𝑧𝑖\displaystyle\sum_{i,l}J_{K;z}[(p_{z;i}^{\dagger}\vec{\sigma}p_{z;i})\cdot\vec% {s}^{d}_{2;i;l}-\frac{1}{2}n^{p}_{z;i}]∑ start_POSTSUBSCRIPT italic_i , italic_l end_POSTSUBSCRIPT italic_J start_POSTSUBSCRIPT italic_K ; italic_z end_POSTSUBSCRIPT [ ( italic_p start_POSTSUBSCRIPT italic_z ; italic_i end_POSTSUBSCRIPT start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT over→ start_ARG italic_σ end_ARG italic_p start_POSTSUBSCRIPT italic_z ; italic_i end_POSTSUBSCRIPT ) ⋅ over→ start_ARG italic_s end_ARG start_POSTSUPERSCRIPT italic_d end_POSTSUPERSCRIPT start_POSTSUBSCRIPT 2 ; italic_i ; italic_l end_POSTSUBSCRIPT - divide start_ARG 1 end_ARG start_ARG 2 end_ARG italic_n start_POSTSUPERSCRIPT italic_p end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_z ; italic_i end_POSTSUBSCRIPT ]
\displaystyle-- i,l2JH[s1;i;lds2;i;ld+14],subscript𝑖𝑙2subscript𝐽𝐻delimited-[]subscriptsuperscript𝑠𝑑1𝑖𝑙subscriptsuperscript𝑠𝑑2𝑖𝑙14\displaystyle\sum_{i,l}2J_{H}[\vec{s}^{d}_{1;i;l}\cdot\vec{s}^{d}_{2;i;l}+% \frac{1}{4}],∑ start_POSTSUBSCRIPT italic_i , italic_l end_POSTSUBSCRIPT 2 italic_J start_POSTSUBSCRIPT italic_H end_POSTSUBSCRIPT [ over→ start_ARG italic_s end_ARG start_POSTSUPERSCRIPT italic_d end_POSTSUPERSCRIPT start_POSTSUBSCRIPT 1 ; italic_i ; italic_l end_POSTSUBSCRIPT ⋅ over→ start_ARG italic_s end_ARG start_POSTSUPERSCRIPT italic_d end_POSTSUPERSCRIPT start_POSTSUBSCRIPT 2 ; italic_i ; italic_l end_POSTSUBSCRIPT + divide start_ARG 1 end_ARG start_ARG 4 end_ARG ] ,

where sa;i;ldsubscriptsuperscript𝑠𝑑𝑎𝑖𝑙\vec{s}^{d}_{a;i;l}over→ start_ARG italic_s end_ARG start_POSTSUPERSCRIPT italic_d end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_a ; italic_i ; italic_l end_POSTSUBSCRIPT is the spin 1/2 operator from the dasubscript𝑑𝑎d_{a}italic_d start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT orbital of the Ni at site i𝑖iitalic_i of layer l=t,b𝑙𝑡𝑏l=t,bitalic_l = italic_t , italic_b. We have

JK;asubscript𝐽𝐾𝑎\displaystyle J_{K;a}italic_J start_POSTSUBSCRIPT italic_K ; italic_a end_POSTSUBSCRIPT =\displaystyle== 4tdp;a2U1+U+JHΔp;a+4tdp;a2Δp;a+JHU,4subscriptsuperscript𝑡2𝑑𝑝𝑎subscript𝑈1superscript𝑈subscript𝐽𝐻subscriptΔ𝑝𝑎4subscriptsuperscript𝑡2𝑑𝑝𝑎subscriptΔ𝑝𝑎subscript𝐽𝐻superscript𝑈\displaystyle\frac{4t^{2}_{dp;a}}{U_{1}+U^{\prime}+J_{H}-\Delta_{p;a}}+\frac{4% t^{2}_{dp;a}}{\Delta_{p;a}+J_{H}-U^{\prime}},divide start_ARG 4 italic_t start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_d italic_p ; italic_a end_POSTSUBSCRIPT end_ARG start_ARG italic_U start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT + italic_U start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT + italic_J start_POSTSUBSCRIPT italic_H end_POSTSUBSCRIPT - roman_Δ start_POSTSUBSCRIPT italic_p ; italic_a end_POSTSUBSCRIPT end_ARG + divide start_ARG 4 italic_t start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_d italic_p ; italic_a end_POSTSUBSCRIPT end_ARG start_ARG roman_Δ start_POSTSUBSCRIPT italic_p ; italic_a end_POSTSUBSCRIPT + italic_J start_POSTSUBSCRIPT italic_H end_POSTSUBSCRIPT - italic_U start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_ARG ,
JK;zsubscript𝐽𝐾𝑧\displaystyle J_{K;z}italic_J start_POSTSUBSCRIPT italic_K ; italic_z end_POSTSUBSCRIPT =\displaystyle== tdp;z2U2+U+JHΔp;2z+tdp;z2Δp;2z+JHU.subscriptsuperscript𝑡2𝑑𝑝𝑧subscript𝑈2superscript𝑈subscript𝐽𝐻superscriptsubscriptΔ𝑝2𝑧subscriptsuperscript𝑡2𝑑𝑝𝑧superscriptsubscriptΔ𝑝2𝑧subscript𝐽𝐻superscript𝑈\displaystyle\frac{t^{2}_{dp;z}}{U_{2}+U^{\prime}+J_{H}-\Delta_{p;2}^{z}}+% \frac{t^{2}_{dp;z}}{\Delta_{p;2}^{z}+J_{H}-U^{\prime}}.divide start_ARG italic_t start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_d italic_p ; italic_z end_POSTSUBSCRIPT end_ARG start_ARG italic_U start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT + italic_U start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT + italic_J start_POSTSUBSCRIPT italic_H end_POSTSUBSCRIPT - roman_Δ start_POSTSUBSCRIPT italic_p ; 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_z end_POSTSUPERSCRIPT end_ARG + divide start_ARG italic_t start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_d italic_p ; italic_z end_POSTSUBSCRIPT end_ARG start_ARG roman_Δ start_POSTSUBSCRIPT italic_p ; 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_z end_POSTSUPERSCRIPT + italic_J start_POSTSUBSCRIPT italic_H end_POSTSUBSCRIPT - italic_U start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_ARG . (3)

with Δp;a=ϵpϵd;asubscriptΔ𝑝𝑎subscriptitalic-ϵ𝑝subscriptitalic-ϵ𝑑𝑎\Delta_{p;a}=\epsilon_{p}-\epsilon_{d;a}roman_Δ start_POSTSUBSCRIPT italic_p ; italic_a end_POSTSUBSCRIPT = italic_ϵ start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT - italic_ϵ start_POSTSUBSCRIPT italic_d ; italic_a end_POSTSUBSCRIPT, Δp;az=ϵp;zϵd;asuperscriptsubscriptΔ𝑝𝑎𝑧subscriptitalic-ϵ𝑝𝑧subscriptitalic-ϵ𝑑𝑎\Delta_{p;a}^{z}=\epsilon_{p;z}-\epsilon_{d;a}roman_Δ start_POSTSUBSCRIPT italic_p ; italic_a end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_z end_POSTSUPERSCRIPT = italic_ϵ start_POSTSUBSCRIPT italic_p ; italic_z end_POSTSUBSCRIPT - italic_ϵ start_POSTSUBSCRIPT italic_d ; italic_a end_POSTSUBSCRIPT. Here, we assume UJHΔp,a,Δp,azUa+Uformulae-sequencemuch-less-thansuperscript𝑈subscript𝐽𝐻subscriptΔ𝑝𝑎much-less-thansuperscriptsubscriptΔ𝑝𝑎𝑧subscript𝑈𝑎superscript𝑈U^{\prime}-J_{H}\ll\Delta_{p,a},\Delta_{p,a}^{z}\ll U_{a}+U^{\prime}italic_U start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT - italic_J start_POSTSUBSCRIPT italic_H end_POSTSUBSCRIPT ≪ roman_Δ start_POSTSUBSCRIPT italic_p , italic_a end_POSTSUBSCRIPT , roman_Δ start_POSTSUBSCRIPT italic_p , italic_a end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_z end_POSTSUPERSCRIPT ≪ italic_U start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT + italic_U start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT so the system is in the charge transfer regime. For illustration, we adopt the hopping parameters calculated at 29.5GPa by DFT, tdp;1=1.56subscript𝑡𝑑𝑝11.56t_{dp;1}=1.56italic_t start_POSTSUBSCRIPT italic_d italic_p ; 1 end_POSTSUBSCRIPT = 1.56eV, tdp;2=0.75subscript𝑡𝑑𝑝20.75t_{dp;2}=-0.75italic_t start_POSTSUBSCRIPT italic_d italic_p ; 2 end_POSTSUBSCRIPT = - 0.75eV, tdp;z=1.63subscript𝑡𝑑𝑝𝑧1.63t_{dp;z}=1.63italic_t start_POSTSUBSCRIPT italic_d italic_p ; italic_z end_POSTSUBSCRIPT = 1.63eV Wú et al. (2023). We assume the similar interaction strength, Ua=10subscript𝑈𝑎10U_{a}=10italic_U start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT = 10eV, U=6superscript𝑈6U^{\prime}=6italic_U start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT = 6eV, JH=2subscript𝐽𝐻2J_{H}=2italic_J start_POSTSUBSCRIPT italic_H end_POSTSUBSCRIPT = 2eV and as in the cuprates Ogata and Fukuyama (2008) and use Δp;a=Δp;az=9subscriptΔ𝑝𝑎superscriptsubscriptΔ𝑝𝑎𝑧9\Delta_{p;a}=\Delta_{p;a}^{z}=9roman_Δ start_POSTSUBSCRIPT italic_p ; italic_a end_POSTSUBSCRIPT = roman_Δ start_POSTSUBSCRIPT italic_p ; italic_a end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_z end_POSTSUPERSCRIPT = 9eV. Within those parameter sets, we found that JK;zJK;1/4similar-to-or-equalssubscript𝐽𝐾𝑧subscript𝐽𝐾14J_{K;z}\simeq J_{K;1}/4italic_J start_POSTSUBSCRIPT italic_K ; italic_z end_POSTSUBSCRIPT ≃ italic_J start_POSTSUBSCRIPT italic_K ; 1 end_POSTSUBSCRIPT / 4 with JK;1=3.03subscript𝐽𝐾13.03J_{K;1}=3.03italic_J start_POSTSUBSCRIPT italic_K ; 1 end_POSTSUBSCRIPT = 3.03eV and JK;z=0.83subscript𝐽𝐾𝑧0.83J_{K;z}=0.83italic_J start_POSTSUBSCRIPT italic_K ; italic_z end_POSTSUBSCRIPT = 0.83eV.

Refer to caption
Figure 2: The wave function of Zhang-Rice spin-1/2 states of the bilayer nickelates. The local d1subscript𝑑1d_{1}italic_d start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT, d2subscript𝑑2d_{2}italic_d start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT orbitals of Ni atoms surrounded by the four in-plane p𝑝pitalic_p orbitals. The yellow, green, and magenta arrow is for the spin of the p1subscript𝑝1p_{1}italic_p start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT, d1subscript𝑑1d_{1}italic_d start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT, and d2subscript𝑑2d_{2}italic_d start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT orbital, respectively. (a) In the JK;1JHmuch-greater-thansubscript𝐽𝐾1subscript𝐽𝐻J_{K;1}\gg J_{H}italic_J start_POSTSUBSCRIPT italic_K ; 1 end_POSTSUBSCRIPT ≫ italic_J start_POSTSUBSCRIPT italic_H end_POSTSUBSCRIPT limit, the Zhang-Rice spin-1/2 state is a simple product state of the Zhang-Rice singlet made by d1,p1subscript𝑑1subscript𝑝1d_{1},p_{1}italic_d start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , italic_p start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT and the decoupled d2subscript𝑑2d_{2}italic_d start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT orbital. (b) In the JK;1JHmuch-less-thansubscript𝐽𝐾1subscript𝐽𝐻J_{K;1}\ll J_{H}italic_J start_POSTSUBSCRIPT italic_K ; 1 end_POSTSUBSCRIPT ≪ italic_J start_POSTSUBSCRIPT italic_H end_POSTSUBSCRIPT limit, two d𝑑ditalic_d orbital forming a spin triplet forms a spin 1/2121/21 / 2 along with p1subscript𝑝1p_{1}italic_p start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT orbitals. Here, for simplicity, we illustrate only the |ket\ket{{\uparrow}}| start_ARG ↑ end_ARG ⟩, since |ket\ket{{\downarrow}}| start_ARG ↓ end_ARG ⟩ is just a its time-reversal partner. (c) The three orbital contributions in Stotz=12subscriptsuperscript𝑆𝑧tot12S^{z}_{\mathrm{tot}}=\frac{1}{2}italic_S start_POSTSUPERSCRIPT italic_z end_POSTSUPERSCRIPT start_POSTSUBSCRIPT roman_tot end_POSTSUBSCRIPT = divide start_ARG 1 end_ARG start_ARG 2 end_ARG of |ket\ket{\uparrow}| start_ARG ↑ end_ARG ⟩ state. As JK;1/JH+subscript𝐽𝐾1subscript𝐽𝐻J_{K;1}/J_{H}\rightarrow+\inftyitalic_J start_POSTSUBSCRIPT italic_K ; 1 end_POSTSUBSCRIPT / italic_J start_POSTSUBSCRIPT italic_H end_POSTSUBSCRIPT → + ∞, the spin-half state is dominated by the d2subscript𝑑2d_{2}italic_d start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT orbital. Generically the state is a combination of all three orbitals: p1subscript𝑝1p_{1}italic_p start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT, d1subscript𝑑1d_{1}italic_d start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT, and d2subscript𝑑2d_{2}italic_d start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT.

We should view p1,p2subscript𝑝1subscript𝑝2p_{1},p_{2}italic_p start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , italic_p start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT as Wannier orbitals centered on Ni atom. The doped hole is favored to enter one of them or the pzsubscript𝑝𝑧p_{z}italic_p start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT orbital by minimizing the Kondo coupling JK;1,JK;2,JK;zsubscript𝐽𝐾1subscript𝐽𝐾2subscript𝐽𝐾𝑧J_{K;1},J_{K;2},J_{K;z}italic_J start_POSTSUBSCRIPT italic_K ; 1 end_POSTSUBSCRIPT , italic_J start_POSTSUBSCRIPT italic_K ; 2 end_POSTSUBSCRIPT , italic_J start_POSTSUBSCRIPT italic_K ; italic_z end_POSTSUBSCRIPT. When the hole occupies the in-plane pasubscript𝑝𝑎p_{a}italic_p start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT orbital, it forms a net S=12𝑆12S=\frac{1}{2}italic_S = divide start_ARG 1 end_ARG start_ARG 2 end_ARG together with the original spin-one moment from the Ni2+. On the other hand, when the hole occupies the pzsubscript𝑝𝑧p_{z}italic_p start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT orbital, it couples to the spin-one moments of Ni from both layers. The ground state energies resulting from the Kondo interactions per unit cell summed over two layers for each case are

EGa=[JH+JH2+JK;aJH+JK;a2+JK;a],superscriptsubscript𝐸𝐺𝑎delimited-[]subscript𝐽𝐻superscriptsubscript𝐽𝐻2subscript𝐽𝐾𝑎subscript𝐽𝐻superscriptsubscript𝐽𝐾𝑎2subscript𝐽𝐾𝑎\displaystyle E_{G}^{a}=-\big{[}J_{H}+\sqrt{J_{H}^{2}+J_{K;a}J_{H}+J_{K;a}^{2}% }+J_{K;a}\big{]},italic_E start_POSTSUBSCRIPT italic_G end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_a end_POSTSUPERSCRIPT = - [ italic_J start_POSTSUBSCRIPT italic_H end_POSTSUBSCRIPT + square-root start_ARG italic_J start_POSTSUBSCRIPT italic_H end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + italic_J start_POSTSUBSCRIPT italic_K ; italic_a end_POSTSUBSCRIPT italic_J start_POSTSUBSCRIPT italic_H end_POSTSUBSCRIPT + italic_J start_POSTSUBSCRIPT italic_K ; italic_a end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG + italic_J start_POSTSUBSCRIPT italic_K ; italic_a end_POSTSUBSCRIPT ] ,
EGz=12[2JH+4JH2+8JK;zJH+9JK;z2+3JK;z].superscriptsubscript𝐸𝐺𝑧12delimited-[]2subscript𝐽𝐻4superscriptsubscript𝐽𝐻28subscript𝐽𝐾𝑧subscript𝐽𝐻9superscriptsubscript𝐽𝐾𝑧23subscript𝐽𝐾𝑧\displaystyle E_{G}^{z}=-\frac{1}{2}\big{[}2J_{H}+\sqrt{4J_{H}^{2}+8J_{K;z}J_{% H}+9{J_{K;z}}^{2}}+3J_{K;z}\big{]}.italic_E start_POSTSUBSCRIPT italic_G end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_z end_POSTSUPERSCRIPT = - divide start_ARG 1 end_ARG start_ARG 2 end_ARG [ 2 italic_J start_POSTSUBSCRIPT italic_H end_POSTSUBSCRIPT + square-root start_ARG 4 italic_J start_POSTSUBSCRIPT italic_H end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + 8 italic_J start_POSTSUBSCRIPT italic_K ; italic_z end_POSTSUBSCRIPT italic_J start_POSTSUBSCRIPT italic_H end_POSTSUBSCRIPT + 9 italic_J start_POSTSUBSCRIPT italic_K ; italic_z end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG + 3 italic_J start_POSTSUBSCRIPT italic_K ; italic_z end_POSTSUBSCRIPT ] .

Using the condition tdp;1tdp;2much-greater-thansubscript𝑡𝑑𝑝1subscript𝑡𝑑𝑝2t_{dp;1}\gg t_{dp;2}italic_t start_POSTSUBSCRIPT italic_d italic_p ; 1 end_POSTSUBSCRIPT ≫ italic_t start_POSTSUBSCRIPT italic_d italic_p ; 2 end_POSTSUBSCRIPT, we always have EG1(JK;1)EG2(JK;2)much-greater-thansuperscriptsubscript𝐸𝐺1subscript𝐽𝐾1superscriptsubscript𝐸𝐺2subscript𝐽𝐾2E_{G}^{1}(J_{K;1})\gg E_{G}^{2}(J_{K;2})italic_E start_POSTSUBSCRIPT italic_G end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 1 end_POSTSUPERSCRIPT ( italic_J start_POSTSUBSCRIPT italic_K ; 1 end_POSTSUBSCRIPT ) ≫ italic_E start_POSTSUBSCRIPT italic_G end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( italic_J start_POSTSUBSCRIPT italic_K ; 2 end_POSTSUBSCRIPT ) and the main competition is between the p1subscript𝑝1p_{1}italic_p start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT and pzsubscript𝑝𝑧p_{z}italic_p start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT orbitals. For the estimated parameter of JK;14JK;zsimilar-to-or-equalssubscript𝐽𝐾14subscript𝐽𝐾𝑧J_{K;1}\simeq 4J_{K;z}italic_J start_POSTSUBSCRIPT italic_K ; 1 end_POSTSUBSCRIPT ≃ 4 italic_J start_POSTSUBSCRIPT italic_K ; italic_z end_POSTSUBSCRIPT, we have EG1(JK;14JK;z)<EGz(JK;z)superscriptsubscript𝐸𝐺1similar-to-or-equalssubscript𝐽𝐾14subscript𝐽𝐾𝑧superscriptsubscript𝐸𝐺𝑧subscript𝐽𝐾𝑧E_{G}^{1}(J_{K;1}\simeq 4J_{K;z})<E_{G}^{z}(J_{K;z})italic_E start_POSTSUBSCRIPT italic_G end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 1 end_POSTSUPERSCRIPT ( italic_J start_POSTSUBSCRIPT italic_K ; 1 end_POSTSUBSCRIPT ≃ 4 italic_J start_POSTSUBSCRIPT italic_K ; italic_z end_POSTSUBSCRIPT ) < italic_E start_POSTSUBSCRIPT italic_G end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_z end_POSTSUPERSCRIPT ( italic_J start_POSTSUBSCRIPT italic_K ; italic_z end_POSTSUBSCRIPT ) and the p1subscript𝑝1p_{1}italic_p start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT orbital wins. In this regime the doped hole enters the in-plane oxygen orbital p1subscript𝑝1p_{1}italic_p start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT and screens the spin-one Ni2+ down to spin-half. We dub this hole state as Zhang-Rice spin-half state as analog to the Zhang-Rice singlet in hole doped cuprate. This Zhang-Rice spin-half consists of the spins from the p1subscript𝑝1p_{1}italic_p start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT orbital and the original d1subscript𝑑1d_{1}italic_d start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT and d2subscript𝑑2d_{2}italic_d start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT orbital of Ni. The wavefunction depends on the parameter r=JK;1/JH𝑟subscript𝐽𝐾1subscript𝐽𝐻r=J_{K;1}/J_{H}italic_r = italic_J start_POSTSUBSCRIPT italic_K ; 1 end_POSTSUBSCRIPT / italic_J start_POSTSUBSCRIPT italic_H end_POSTSUBSCRIPT:

|σl(1+r+α)|σ¯,σ,σ(r+α)|σ,σ¯,σ|σ,σ,σ¯,similar-tosubscriptket𝜎𝑙1𝑟𝛼ket¯𝜎𝜎𝜎𝑟𝛼ket𝜎¯𝜎𝜎ket𝜎𝜎¯𝜎\displaystyle\ket{\sigma}_{l}\sim(1+r+\alpha)\ket{\overline{\sigma},\sigma,% \sigma}-(r+\alpha)\ket{\sigma,\overline{\sigma},\sigma}-\ket{\sigma,\sigma,% \overline{\sigma}},\!| start_ARG italic_σ end_ARG ⟩ start_POSTSUBSCRIPT italic_l end_POSTSUBSCRIPT ∼ ( 1 + italic_r + italic_α ) | start_ARG over¯ start_ARG italic_σ end_ARG , italic_σ , italic_σ end_ARG ⟩ - ( italic_r + italic_α ) | start_ARG italic_σ , over¯ start_ARG italic_σ end_ARG , italic_σ end_ARG ⟩ - | start_ARG italic_σ , italic_σ , over¯ start_ARG italic_σ end_ARG end_ARG ⟩ , (4)

with σ=±12𝜎plus-or-minus12\sigma=\pm\frac{1}{2}italic_σ = ± divide start_ARG 1 end_ARG start_ARG 2 end_ARG. We used α(r)=[r2+r+1]1/2𝛼𝑟superscriptdelimited-[]superscript𝑟2𝑟112\alpha(r)=[r^{2}+r+1]^{1/2}italic_α ( italic_r ) = [ italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + italic_r + 1 ] start_POSTSUPERSCRIPT 1 / 2 end_POSTSUPERSCRIPT and here omitted the normalization factor, for simplicity (See SM Sec.II). Also, σ¯=σ¯𝜎𝜎\overline{\sigma}=-\sigmaover¯ start_ARG italic_σ end_ARG = - italic_σ, and |σ,σ1,σ2ket𝜎subscript𝜎1subscript𝜎2\ket{\sigma,\sigma_{1},\sigma_{2}}| start_ARG italic_σ , italic_σ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , italic_σ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_ARG ⟩ denotes the spin of p1;l,d1;l,d2;lsubscript𝑝1𝑙subscript𝑑1𝑙subscript𝑑2𝑙p_{1;l},d_{1;l},d_{2;l}italic_p start_POSTSUBSCRIPT 1 ; italic_l end_POSTSUBSCRIPT , italic_d start_POSTSUBSCRIPT 1 ; italic_l end_POSTSUBSCRIPT , italic_d start_POSTSUBSCRIPT 2 ; italic_l end_POSTSUBSCRIPT orbital respectively. A few remarks are as follows. In the limit of JHJK;1much-less-thansubscript𝐽𝐻subscript𝐽𝐾1J_{H}\ll J_{K;1}italic_J start_POSTSUBSCRIPT italic_H end_POSTSUBSCRIPT ≪ italic_J start_POSTSUBSCRIPT italic_K ; 1 end_POSTSUBSCRIPT, the state is simply the tensor product of the Zhang-Rice singlet from only d1,p1subscript𝑑1subscript𝑝1d_{1},p_{1}italic_d start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , italic_p start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT orbitals and a decoupled spin-half from the d2subscript𝑑2d_{2}italic_d start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT orbitals shown in Fig. 2(a). However, with a general JHsubscript𝐽𝐻J_{H}italic_J start_POSTSUBSCRIPT italic_H end_POSTSUBSCRIPT, the three orbitals are highly entangled and should be considered together. For instance, in Fig. 2(c), we plot the spin contributions of the net spin-1/2 from all three orbitals. The bipartite entropy between the d1,p1subscript𝑑1subscript𝑝1d_{1},p_{1}italic_d start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , italic_p start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT orbitals and the d2subscript𝑑2d_{2}italic_d start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT orbital increases with increasing JHsubscript𝐽𝐻J_{H}italic_J start_POSTSUBSCRIPT italic_H end_POSTSUBSCRIPT, reaching its maximum values in the limit where JHJK;1much-greater-thansubscript𝐽𝐻subscript𝐽𝐾1J_{H}\gg J_{K;1}italic_J start_POSTSUBSCRIPT italic_H end_POSTSUBSCRIPT ≫ italic_J start_POSTSUBSCRIPT italic_K ; 1 end_POSTSUBSCRIPT (See SM), where we illustrate the state in Fig. 2(b).

The state |σlsubscriptket𝜎𝑙|{\sigma}\rangle_{l}| italic_σ ⟩ start_POSTSUBSCRIPT italic_l end_POSTSUBSCRIPT is derived from solving the Hamiltonian in a single unit cell, next we should make the states at two different sites i,j𝑖𝑗i,jitalic_i , italic_j orthogonal. We just call this state as |σi;lsubscriptket𝜎𝑖𝑙|{\sigma}\rangle_{i;l}| italic_σ ⟩ start_POSTSUBSCRIPT italic_i ; italic_l end_POSTSUBSCRIPT. This can be simply done by replacing the p1subscript𝑝1p_{1}italic_p start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT orbital in Eq.(4) with a wannierized electron operator, ci;l;σ=𝒩jB(ij)p1;j;l;σsuperscriptsubscript𝑐𝑖𝑙𝜎superscript𝒩subscript𝑗𝐵𝑖𝑗subscript𝑝1𝑗𝑙𝜎c_{i;l;\sigma}^{\dagger}=\mathcal{N^{\prime}}\sum_{j}B(i-j)p_{1;j;l;\sigma}italic_c start_POSTSUBSCRIPT italic_i ; italic_l ; italic_σ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT = caligraphic_N start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ∑ start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT italic_B ( italic_i - italic_j ) italic_p start_POSTSUBSCRIPT 1 ; italic_j ; italic_l ; italic_σ end_POSTSUBSCRIPT where B(ij)=1Nkeik(ij)β(k)𝐵𝑖𝑗1𝑁subscript𝑘superscript𝑒𝑖𝑘𝑖𝑗𝛽𝑘B(i-j)=\frac{1}{N}\sum_{k}e^{ik\cdot(i-j)}\beta(k)italic_B ( italic_i - italic_j ) = divide start_ARG 1 end_ARG start_ARG italic_N end_ARG ∑ start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT italic_e start_POSTSUPERSCRIPT italic_i italic_k ⋅ ( italic_i - italic_j ) end_POSTSUPERSCRIPT italic_β ( italic_k ) and β(k)=[112(coskx+cosky)]1/2𝛽𝑘superscriptdelimited-[]112subscript𝑘𝑥subscript𝑘𝑦12\beta(k)=[1-\frac{1}{2}(\cos k_{x}+\cos k_{y})]^{-1/2}italic_β ( italic_k ) = [ 1 - divide start_ARG 1 end_ARG start_ARG 2 end_ARG ( roman_cos italic_k start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT + roman_cos italic_k start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT ) ] start_POSTSUPERSCRIPT - 1 / 2 end_POSTSUPERSCRIPT which form orthogonal and complete sets with {ci;l;σ,cj;l;σ}=δi,jδl,lδσ,σ𝒩2subscript𝑐𝑖𝑙𝜎superscriptsubscript𝑐𝑗superscript𝑙superscript𝜎subscript𝛿𝑖𝑗subscript𝛿𝑙superscript𝑙subscript𝛿𝜎superscript𝜎superscriptsuperscript𝒩2\{c_{i;l;\sigma},c_{j;l^{\prime};\sigma^{\prime}}^{\dagger}\}=\delta_{i,j}% \delta_{l,l^{\prime}}\delta_{\sigma,\sigma^{\prime}}\mathcal{N^{\prime}}^{2}{ italic_c start_POSTSUBSCRIPT italic_i ; italic_l ; italic_σ end_POSTSUBSCRIPT , italic_c start_POSTSUBSCRIPT italic_j ; italic_l start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ; italic_σ start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT } = italic_δ start_POSTSUBSCRIPT italic_i , italic_j end_POSTSUBSCRIPT italic_δ start_POSTSUBSCRIPT italic_l , italic_l start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT italic_δ start_POSTSUBSCRIPT italic_σ , italic_σ start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT caligraphic_N start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT. Here, we use 𝒩(r)=[2(r+α)2+2(r+α)+2]1/2/(r+α+1)superscript𝒩𝑟superscriptdelimited-[]2superscript𝑟𝛼22𝑟𝛼212𝑟𝛼1\mathcal{N}^{\prime}(r)=[2(r+\alpha)^{2}+2(r+\alpha)+2]^{1/2}/(r+\alpha+1)caligraphic_N start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ( italic_r ) = [ 2 ( italic_r + italic_α ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + 2 ( italic_r + italic_α ) + 2 ] start_POSTSUPERSCRIPT 1 / 2 end_POSTSUPERSCRIPT / ( italic_r + italic_α + 1 ), monotonic increasing function as increasing r𝑟ritalic_r, ranging [1.22,1.41]1.221.41[1.22,1.41][ 1.22 , 1.41 ]. Then, we can turn to the electron picture in the remaining part of the paper, which is more intuitive, through using ci;l;σsubscript𝑐𝑖𝑙𝜎c_{i;l;\sigma}italic_c start_POSTSUBSCRIPT italic_i ; italic_l ; italic_σ end_POSTSUBSCRIPT. In this context, we have ni;c=1xsubscript𝑛𝑖𝑐1𝑥n_{i;c}=1-xitalic_n start_POSTSUBSCRIPT italic_i ; italic_c end_POSTSUBSCRIPT = 1 - italic_x, which corresponds to the condition ni;p=xsubscript𝑛𝑖𝑝𝑥n_{i;p}=xitalic_n start_POSTSUBSCRIPT italic_i ; italic_p end_POSTSUBSCRIPT = italic_x in the hole representation.

Type-II t-J model.— Now we can see that the low energy model in the charge transfer regime is still the type II t-J model Zhang and Vishwanath (2020). The undoped Ni site is in the d8superscript𝑑8d^{8}italic_d start_POSTSUPERSCRIPT 8 end_POSTSUPERSCRIPT with spin-one, while the single hole d8Lsuperscript𝑑8𝐿d^{8}Litalic_d start_POSTSUPERSCRIPT 8 end_POSTSUPERSCRIPT italic_L state corresponds to the Zhang-Rice spin-half state above, which is an analog of the d7superscript𝑑7d^{7}italic_d start_POSTSUPERSCRIPT 7 end_POSTSUPERSCRIPT states in the Mott-Hubbard regime given that the p1subscript𝑝1p_{1}italic_p start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT orbital has the same crystal symmetry as the d1subscript𝑑1d_{1}italic_d start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT orbital. Then, the wannierized orbital ci;l;σsubscript𝑐𝑖𝑙𝜎c_{i;l;\sigma}italic_c start_POSTSUBSCRIPT italic_i ; italic_l ; italic_σ end_POSTSUBSCRIPT can be projected in the 2+3=52352+3=52 + 3 = 5 dimensional Hilbert space. We now turn to the electron picture for simplicity. For a general r=JK;1/JH𝑟subscript𝐽𝐾1subscript𝐽𝐻r=J_{K;1}/J_{H}italic_r = italic_J start_POSTSUBSCRIPT italic_K ; 1 end_POSTSUBSCRIPT / italic_J start_POSTSUBSCRIPT italic_H end_POSTSUBSCRIPT, the electron operator of the p1subscript𝑝1p_{1}italic_p start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT orbital is written as,

ci;l;subscript𝑐𝑖𝑙\displaystyle c_{i;l;\uparrow}italic_c start_POSTSUBSCRIPT italic_i ; italic_l ; ↑ end_POSTSUBSCRIPT =\displaystyle== +j<i(1)nj[|i;l1|i;l+12|i;l0|i;l],subscriptproduct𝑗𝑖superscript1subscript𝑛𝑗delimited-[]subscriptket𝑖𝑙subscriptbra1𝑖𝑙12subscriptket𝑖𝑙subscriptbra0𝑖𝑙\displaystyle+\prod_{j<i}{(-1)^{n_{j}}}\left[\ket{\downarrow}_{i;l}\bra{-1}_{i% ;l}+\frac{1}{\sqrt{2}}\ket{\uparrow}_{i;l}\bra{0}_{i;l}\right],+ ∏ start_POSTSUBSCRIPT italic_j < italic_i end_POSTSUBSCRIPT ( - 1 ) start_POSTSUPERSCRIPT italic_n start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT end_POSTSUPERSCRIPT [ | start_ARG ↓ end_ARG ⟩ start_POSTSUBSCRIPT italic_i ; italic_l end_POSTSUBSCRIPT ⟨ start_ARG - 1 end_ARG | start_POSTSUBSCRIPT italic_i ; italic_l end_POSTSUBSCRIPT + divide start_ARG 1 end_ARG start_ARG square-root start_ARG 2 end_ARG end_ARG | start_ARG ↑ end_ARG ⟩ start_POSTSUBSCRIPT italic_i ; italic_l end_POSTSUBSCRIPT ⟨ start_ARG 0 end_ARG | start_POSTSUBSCRIPT italic_i ; italic_l end_POSTSUBSCRIPT ] ,
ci;l;subscript𝑐𝑖𝑙\displaystyle c_{i;l;\downarrow}italic_c start_POSTSUBSCRIPT italic_i ; italic_l ; ↓ end_POSTSUBSCRIPT =\displaystyle== j<i(1)nj[|i;l1|i;l+12|i;l0|i;l],subscriptproduct𝑗𝑖superscript1subscript𝑛𝑗delimited-[]subscriptket𝑖𝑙subscriptbra1𝑖𝑙12subscriptket𝑖𝑙subscriptbra0𝑖𝑙\displaystyle-\prod_{j<i}{(-1)^{n_{j}}}\left[\ket{\uparrow}_{i;l}\bra{1}_{i;l}% +\frac{1}{\sqrt{2}}\ket{\downarrow}_{i;l}\bra{0}_{i;l}\right],- ∏ start_POSTSUBSCRIPT italic_j < italic_i end_POSTSUBSCRIPT ( - 1 ) start_POSTSUPERSCRIPT italic_n start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT end_POSTSUPERSCRIPT [ | start_ARG ↑ end_ARG ⟩ start_POSTSUBSCRIPT italic_i ; italic_l end_POSTSUBSCRIPT ⟨ start_ARG 1 end_ARG | start_POSTSUBSCRIPT italic_i ; italic_l end_POSTSUBSCRIPT + divide start_ARG 1 end_ARG start_ARG square-root start_ARG 2 end_ARG end_ARG | start_ARG ↓ end_ARG ⟩ start_POSTSUBSCRIPT italic_i ; italic_l end_POSTSUBSCRIPT ⟨ start_ARG 0 end_ARG | start_POSTSUBSCRIPT italic_i ; italic_l end_POSTSUBSCRIPT ] ,

with the Jordan-Wigner string, j<i(1)njsubscriptproduct𝑗𝑖superscript1subscript𝑛𝑗\prod_{j<i}(-1)^{n_{j}}∏ start_POSTSUBSCRIPT italic_j < italic_i end_POSTSUBSCRIPT ( - 1 ) start_POSTSUPERSCRIPT italic_n start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT end_POSTSUPERSCRIPT. The spin-triplet states of d8superscript𝑑8d^{8}italic_d start_POSTSUPERSCRIPT 8 end_POSTSUPERSCRIPT are |1l=d1;ld2;l;|Gsubscriptket1𝑙subscriptsuperscript𝑑1𝑙absentsubscriptsuperscript𝑑2𝑙ket𝐺\ket{-1}_{l}=d^{\dagger}_{1;l\downarrow}d^{\dagger}_{2;l;\downarrow}\ket{G}| start_ARG - 1 end_ARG ⟩ start_POSTSUBSCRIPT italic_l end_POSTSUBSCRIPT = italic_d start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT start_POSTSUBSCRIPT 1 ; italic_l ↓ end_POSTSUBSCRIPT italic_d start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT start_POSTSUBSCRIPT 2 ; italic_l ; ↓ end_POSTSUBSCRIPT | start_ARG italic_G end_ARG ⟩, |0l=12(d1;l;d2;l;+d1;l;d2;l;)|Gsubscriptket0𝑙12subscriptsuperscript𝑑1𝑙subscriptsuperscript𝑑2𝑙subscriptsuperscript𝑑1𝑙subscriptsuperscript𝑑2𝑙ket𝐺\ket{0}_{l}=\frac{1}{\sqrt{2}}(d^{\dagger}_{1;l;\uparrow}d^{\dagger}_{2;l;% \downarrow}+d^{\dagger}_{1;l;\downarrow}d^{\dagger}_{2;l;\uparrow})\ket{G}| start_ARG 0 end_ARG ⟩ start_POSTSUBSCRIPT italic_l end_POSTSUBSCRIPT = divide start_ARG 1 end_ARG start_ARG square-root start_ARG 2 end_ARG end_ARG ( italic_d start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT start_POSTSUBSCRIPT 1 ; italic_l ; ↑ end_POSTSUBSCRIPT italic_d start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT start_POSTSUBSCRIPT 2 ; italic_l ; ↓ end_POSTSUBSCRIPT + italic_d start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT start_POSTSUBSCRIPT 1 ; italic_l ; ↓ end_POSTSUBSCRIPT italic_d start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT start_POSTSUBSCRIPT 2 ; italic_l ; ↑ end_POSTSUBSCRIPT ) | start_ARG italic_G end_ARG ⟩ and |1l=d1;l;d2;l;|Gsubscriptket1𝑙subscriptsuperscript𝑑1𝑙subscriptsuperscript𝑑2𝑙ket𝐺\ket{1}_{l}=d^{\dagger}_{1;l;\uparrow}d^{\dagger}_{2;l;\uparrow}\ket{G}| start_ARG 1 end_ARG ⟩ start_POSTSUBSCRIPT italic_l end_POSTSUBSCRIPT = italic_d start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT start_POSTSUBSCRIPT 1 ; italic_l ; ↑ end_POSTSUBSCRIPT italic_d start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT start_POSTSUBSCRIPT 2 ; italic_l ; ↑ end_POSTSUBSCRIPT | start_ARG italic_G end_ARG ⟩. Note that da;l,σsubscript𝑑𝑎𝑙𝜎d_{a;l,\sigma}italic_d start_POSTSUBSCRIPT italic_a ; italic_l , italic_σ end_POSTSUBSCRIPT is a hole-operator of d𝑑ditalic_d orbital and |Gket𝐺\ket{G}| start_ARG italic_G end_ARG ⟩ is defined as vacuum as a 3d103superscript𝑑103d^{10}3 italic_d start_POSTSUPERSCRIPT 10 end_POSTSUPERSCRIPT configuration.

Then, the charge-transferred type II t-J model written in the electron picture is,

H=TK𝐻subscript𝑇𝐾\displaystyle H=T_{K}italic_H = italic_T start_POSTSUBSCRIPT italic_K end_POSTSUBSCRIPT +\displaystyle++ Jsslijsi;lsj;l+JddlijSi;lSj;lsuperscriptsubscript𝐽parallel-to𝑠𝑠subscript𝑙subscriptdelimited-⟨⟩𝑖𝑗subscript𝑠𝑖𝑙subscript𝑠𝑗𝑙subscriptsuperscript𝐽𝑑𝑑parallel-tosubscript𝑙subscriptdelimited-⟨⟩𝑖𝑗subscript𝑆𝑖𝑙subscript𝑆𝑗𝑙\displaystyle J_{\parallel}^{ss}\sum_{l}\sum_{\langle ij\rangle}\vec{s}_{i;l}% \cdot\vec{s}_{j;l}+J^{dd}_{\parallel}\sum_{l}\sum_{\langle ij\rangle}\vec{S}_{% i;l}\cdot\vec{S}_{j;l}italic_J start_POSTSUBSCRIPT ∥ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_s italic_s end_POSTSUPERSCRIPT ∑ start_POSTSUBSCRIPT italic_l end_POSTSUBSCRIPT ∑ start_POSTSUBSCRIPT ⟨ italic_i italic_j ⟩ end_POSTSUBSCRIPT over→ start_ARG italic_s end_ARG start_POSTSUBSCRIPT italic_i ; italic_l end_POSTSUBSCRIPT ⋅ over→ start_ARG italic_s end_ARG start_POSTSUBSCRIPT italic_j ; italic_l end_POSTSUBSCRIPT + italic_J start_POSTSUPERSCRIPT italic_d italic_d end_POSTSUPERSCRIPT start_POSTSUBSCRIPT ∥ end_POSTSUBSCRIPT ∑ start_POSTSUBSCRIPT italic_l end_POSTSUBSCRIPT ∑ start_POSTSUBSCRIPT ⟨ italic_i italic_j ⟩ end_POSTSUBSCRIPT over→ start_ARG italic_S end_ARG start_POSTSUBSCRIPT italic_i ; italic_l end_POSTSUBSCRIPT ⋅ over→ start_ARG italic_S end_ARG start_POSTSUBSCRIPT italic_j ; italic_l end_POSTSUBSCRIPT
+\displaystyle++ Jsdlij(si;lSj;l+Si;lsj;l)subscriptsuperscript𝐽𝑠𝑑parallel-tosubscript𝑙subscriptdelimited-⟨⟩𝑖𝑗subscript𝑠𝑖𝑙subscript𝑆𝑗𝑙subscript𝑆𝑖𝑙subscript𝑠𝑗𝑙\displaystyle J^{sd}_{\parallel}\sum_{l}\sum_{\langle ij\rangle}(\vec{s}_{i;l}% \cdot\vec{S}_{j;l}+\vec{S}_{i;l}\cdot\vec{s}_{j;l})italic_J start_POSTSUPERSCRIPT italic_s italic_d end_POSTSUPERSCRIPT start_POSTSUBSCRIPT ∥ end_POSTSUBSCRIPT ∑ start_POSTSUBSCRIPT italic_l end_POSTSUBSCRIPT ∑ start_POSTSUBSCRIPT ⟨ italic_i italic_j ⟩ end_POSTSUBSCRIPT ( over→ start_ARG italic_s end_ARG start_POSTSUBSCRIPT italic_i ; italic_l end_POSTSUBSCRIPT ⋅ over→ start_ARG italic_S end_ARG start_POSTSUBSCRIPT italic_j ; italic_l end_POSTSUBSCRIPT + over→ start_ARG italic_S end_ARG start_POSTSUBSCRIPT italic_i ; italic_l end_POSTSUBSCRIPT ⋅ over→ start_ARG italic_s end_ARG start_POSTSUBSCRIPT italic_j ; italic_l end_POSTSUBSCRIPT )
+\displaystyle++ Jssisi;tsi;b+JddiSi;tSi;bsuperscriptsubscript𝐽perpendicular-to𝑠𝑠subscript𝑖subscript𝑠𝑖𝑡subscript𝑠𝑖𝑏subscriptsuperscript𝐽𝑑𝑑perpendicular-tosubscript𝑖subscript𝑆𝑖𝑡subscript𝑆𝑖𝑏\displaystyle J_{\perp}^{ss}\sum_{i}\vec{s}_{i;t}\cdot\vec{s}_{i;b}+J^{dd}_{% \perp}\sum_{i}\vec{S}_{i;t}\cdot\vec{S}_{i;b}italic_J start_POSTSUBSCRIPT ⟂ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_s italic_s end_POSTSUPERSCRIPT ∑ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT over→ start_ARG italic_s end_ARG start_POSTSUBSCRIPT italic_i ; italic_t end_POSTSUBSCRIPT ⋅ over→ start_ARG italic_s end_ARG start_POSTSUBSCRIPT italic_i ; italic_b end_POSTSUBSCRIPT + italic_J start_POSTSUPERSCRIPT italic_d italic_d end_POSTSUPERSCRIPT start_POSTSUBSCRIPT ⟂ end_POSTSUBSCRIPT ∑ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT over→ start_ARG italic_S end_ARG start_POSTSUBSCRIPT italic_i ; italic_t end_POSTSUBSCRIPT ⋅ over→ start_ARG italic_S end_ARG start_POSTSUBSCRIPT italic_i ; italic_b end_POSTSUBSCRIPT
+\displaystyle++ Jsdi(si;tSi;b+Si;tsi;b)+Vini;tni;b,subscriptsuperscript𝐽𝑠𝑑perpendicular-tosubscript𝑖subscript𝑠𝑖𝑡subscript𝑆𝑖𝑏subscript𝑆𝑖𝑡subscript𝑠𝑖𝑏𝑉subscript𝑖subscript𝑛𝑖𝑡subscript𝑛𝑖𝑏\displaystyle J^{sd}_{\perp}\sum_{i}(\vec{s}_{i;t}\cdot\vec{S}_{i;b}+\vec{S}_{% i;t}\cdot\vec{s}_{i;b})+V\sum_{i}n_{i;t}n_{i;b},italic_J start_POSTSUPERSCRIPT italic_s italic_d end_POSTSUPERSCRIPT start_POSTSUBSCRIPT ⟂ end_POSTSUBSCRIPT ∑ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ( over→ start_ARG italic_s end_ARG start_POSTSUBSCRIPT italic_i ; italic_t end_POSTSUBSCRIPT ⋅ over→ start_ARG italic_S end_ARG start_POSTSUBSCRIPT italic_i ; italic_b end_POSTSUBSCRIPT + over→ start_ARG italic_S end_ARG start_POSTSUBSCRIPT italic_i ; italic_t end_POSTSUBSCRIPT ⋅ over→ start_ARG italic_s end_ARG start_POSTSUBSCRIPT italic_i ; italic_b end_POSTSUBSCRIPT ) + italic_V ∑ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT italic_n start_POSTSUBSCRIPT italic_i ; italic_t end_POSTSUBSCRIPT italic_n start_POSTSUBSCRIPT italic_i ; italic_b end_POSTSUBSCRIPT ,

with

TKsubscript𝑇𝐾\displaystyle T_{K}italic_T start_POSTSUBSCRIPT italic_K end_POSTSUBSCRIPT =\displaystyle== tl,σ,i,jci;l;σcj;l;σtσ,ici;t;σci;b;σ+H.c.,formulae-sequencesubscript𝑡parallel-tosubscript𝑙𝜎𝑖𝑗superscriptsubscript𝑐𝑖𝑙𝜎subscript𝑐𝑗𝑙𝜎subscript𝑡perpendicular-tosubscript𝜎𝑖superscriptsubscript𝑐𝑖𝑡𝜎subscript𝑐𝑖𝑏𝜎Hc\displaystyle-t_{\parallel}\sum_{l,\sigma,\langle i,j\rangle}c_{i;l;\sigma}^{% \dagger}c_{j;l;\sigma}-t_{\perp}\sum_{\sigma,i}c_{i;t;\sigma}^{\dagger}c_{i;b;% \sigma}+\mathrm{H.c.},- italic_t start_POSTSUBSCRIPT ∥ end_POSTSUBSCRIPT ∑ start_POSTSUBSCRIPT italic_l , italic_σ , ⟨ italic_i , italic_j ⟩ end_POSTSUBSCRIPT italic_c start_POSTSUBSCRIPT italic_i ; italic_l ; italic_σ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT italic_c start_POSTSUBSCRIPT italic_j ; italic_l ; italic_σ end_POSTSUBSCRIPT - italic_t start_POSTSUBSCRIPT ⟂ end_POSTSUBSCRIPT ∑ start_POSTSUBSCRIPT italic_σ , italic_i end_POSTSUBSCRIPT italic_c start_POSTSUBSCRIPT italic_i ; italic_t ; italic_σ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT italic_c start_POSTSUBSCRIPT italic_i ; italic_b ; italic_σ end_POSTSUBSCRIPT + roman_H . roman_c . ,

The spin operators for the spin-1/2 state are si;l=12σσ|σi;lσσσσ|i;lsubscript𝑠𝑖𝑙12subscript𝜎superscript𝜎subscriptket𝜎𝑖𝑙subscript𝜎𝜎superscript𝜎subscriptbrasuperscript𝜎𝑖𝑙\vec{s}_{i;l}=\frac{1}{2}\sum_{\sigma\sigma^{\prime}}\ket{\sigma}_{i;l}\vec{% \sigma}_{\sigma\sigma^{\prime}}\bra{\sigma^{\prime}}_{i;l}over→ start_ARG italic_s end_ARG start_POSTSUBSCRIPT italic_i ; italic_l end_POSTSUBSCRIPT = divide start_ARG 1 end_ARG start_ARG 2 end_ARG ∑ start_POSTSUBSCRIPT italic_σ italic_σ start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT | start_ARG italic_σ end_ARG ⟩ start_POSTSUBSCRIPT italic_i ; italic_l end_POSTSUBSCRIPT over→ start_ARG italic_σ end_ARG start_POSTSUBSCRIPT italic_σ italic_σ start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT ⟨ start_ARG italic_σ start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_ARG | start_POSTSUBSCRIPT italic_i ; italic_l end_POSTSUBSCRIPT. Meanwhile, the spin operators for the spin-one moment are written as Si;ld=α,β=1,0,1Tαβ]|αi;lβ|i;l\vec{S}_{i;l}^{d}=\sum_{\alpha,\beta=-1,0,1}\vec{T}_{\alpha\beta]}\ket{\alpha}% _{i;l}\bra{\beta}_{i;l}over→ start_ARG italic_S end_ARG start_POSTSUBSCRIPT italic_i ; italic_l end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_d end_POSTSUPERSCRIPT = ∑ start_POSTSUBSCRIPT italic_α , italic_β = - 1 , 0 , 1 end_POSTSUBSCRIPT over→ start_ARG italic_T end_ARG start_POSTSUBSCRIPT italic_α italic_β ] end_POSTSUBSCRIPT | start_ARG italic_α end_ARG ⟩ start_POSTSUBSCRIPT italic_i ; italic_l end_POSTSUBSCRIPT ⟨ start_ARG italic_β end_ARG | start_POSTSUBSCRIPT italic_i ; italic_l end_POSTSUBSCRIPT. We have Tx=12(010101010)superscript𝑇𝑥12matrix010101010T^{x}=\frac{1}{\sqrt{2}}\begin{pmatrix}0&1&0\\ 1&0&1\\ 0&1&0\end{pmatrix}italic_T start_POSTSUPERSCRIPT italic_x end_POSTSUPERSCRIPT = divide start_ARG 1 end_ARG start_ARG square-root start_ARG 2 end_ARG end_ARG ( start_ARG start_ROW start_CELL 0 end_CELL start_CELL 1 end_CELL start_CELL 0 end_CELL end_ROW start_ROW start_CELL 1 end_CELL start_CELL 0 end_CELL start_CELL 1 end_CELL end_ROW start_ROW start_CELL 0 end_CELL start_CELL 1 end_CELL start_CELL 0 end_CELL end_ROW end_ARG ) and Ty=12(0i0i0i0i0)superscript𝑇𝑦12matrix0𝑖0𝑖0𝑖0𝑖0T^{y}=\frac{1}{\sqrt{2}}\begin{pmatrix}0&-i&0\\ i&0&-i\\ 0&i&0\end{pmatrix}italic_T start_POSTSUPERSCRIPT italic_y end_POSTSUPERSCRIPT = divide start_ARG 1 end_ARG start_ARG square-root start_ARG 2 end_ARG end_ARG ( start_ARG start_ROW start_CELL 0 end_CELL start_CELL - italic_i end_CELL start_CELL 0 end_CELL end_ROW start_ROW start_CELL italic_i end_CELL start_CELL 0 end_CELL start_CELL - italic_i end_CELL end_ROW start_ROW start_CELL 0 end_CELL start_CELL italic_i end_CELL start_CELL 0 end_CELL end_ROW end_ARG ) in the Sz=1,0,1superscript𝑆𝑧101S^{z}=1,0,-1italic_S start_POSTSUPERSCRIPT italic_z end_POSTSUPERSCRIPT = 1 , 0 , - 1 basis. V𝑉Vitalic_V is the inter-layer repulsive density interaction. The expressions for Jsubscript𝐽perpendicular-toJ_{\perp}italic_J start_POSTSUBSCRIPT ⟂ end_POSTSUBSCRIPT, Jsubscript𝐽parallel-toJ_{\parallel}italic_J start_POSTSUBSCRIPT ∥ end_POSTSUBSCRIPT in terms of the microscopic parameters are provided with the derivation in SM Sec.III. Additionally, we establish the relationships Jss=𝒞(r)Jsd=𝒞(r)2Jddsubscriptsuperscript𝐽𝑠𝑠perpendicular-to𝒞𝑟subscriptsuperscript𝐽𝑠𝑑perpendicular-to𝒞superscript𝑟2subscriptsuperscript𝐽𝑑𝑑perpendicular-toJ^{ss}_{\perp}=\mathcal{C}(r)J^{sd}_{\perp}=\mathcal{C}(r)^{2}J^{dd}_{\perp}italic_J start_POSTSUPERSCRIPT italic_s italic_s end_POSTSUPERSCRIPT start_POSTSUBSCRIPT ⟂ end_POSTSUBSCRIPT = caligraphic_C ( italic_r ) italic_J start_POSTSUPERSCRIPT italic_s italic_d end_POSTSUPERSCRIPT start_POSTSUBSCRIPT ⟂ end_POSTSUBSCRIPT = caligraphic_C ( italic_r ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_J start_POSTSUPERSCRIPT italic_d italic_d end_POSTSUPERSCRIPT start_POSTSUBSCRIPT ⟂ end_POSTSUBSCRIPT and Jss=𝒞(r)𝒜Jsd=𝒞(r)2𝒜2Jddsubscriptsuperscript𝐽𝑠𝑠parallel-to𝒞𝑟𝒜subscriptsuperscript𝐽𝑠𝑑parallel-to𝒞superscript𝑟2superscript𝒜2subscriptsuperscript𝐽𝑑𝑑parallel-toJ^{ss}_{\parallel}=\mathcal{C}(r)\mathcal{A}J^{sd}_{\parallel}=\mathcal{C}(r)^% {2}\mathcal{A}^{2}J^{dd}_{\parallel}italic_J start_POSTSUPERSCRIPT italic_s italic_s end_POSTSUPERSCRIPT start_POSTSUBSCRIPT ∥ end_POSTSUBSCRIPT = caligraphic_C ( italic_r ) caligraphic_A italic_J start_POSTSUPERSCRIPT italic_s italic_d end_POSTSUPERSCRIPT start_POSTSUBSCRIPT ∥ end_POSTSUBSCRIPT = caligraphic_C ( italic_r ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT caligraphic_A start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_J start_POSTSUPERSCRIPT italic_d italic_d end_POSTSUPERSCRIPT start_POSTSUBSCRIPT ∥ end_POSTSUBSCRIPT, where 𝒜0.7705similar-to-or-equals𝒜0.7705\mathcal{A}\simeq 0.7705caligraphic_A ≃ 0.7705. Using the parameters mentioned above, we have Jdd=0.097superscriptsubscript𝐽parallel-to𝑑𝑑0.097J_{\parallel}^{dd}=0.097italic_J start_POSTSUBSCRIPT ∥ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_d italic_d end_POSTSUPERSCRIPT = 0.097eV, Jddsuperscriptsubscript𝐽perpendicular-to𝑑𝑑J_{\perp}^{dd}italic_J start_POSTSUBSCRIPT ⟂ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_d italic_d end_POSTSUPERSCRIPT=0.076eV, and tsubscript𝑡parallel-tot_{\parallel}italic_t start_POSTSUBSCRIPT ∥ end_POSTSUBSCRIPT is in the range of 0.320.490.320.490.32-0.490.32 - 0.49eV. 𝒞(r)𝒞𝑟\mathcal{C}(r)caligraphic_C ( italic_r ) ranges from 4/3434/34 / 3 at r0𝑟0r\rightarrow 0italic_r → 0 limit and 00 at r𝑟r\rightarrow\inftyitalic_r → ∞ limit. Here r=JK;1/JH𝑟subscript𝐽𝐾1subscript𝐽𝐻r=J_{K;1}/J_{H}italic_r = italic_J start_POSTSUBSCRIPT italic_K ; 1 end_POSTSUBSCRIPT / italic_J start_POSTSUBSCRIPT italic_H end_POSTSUBSCRIPT. The dependence of t(r)subscript𝑡parallel-to𝑟t_{\parallel}(r)italic_t start_POSTSUBSCRIPT ∥ end_POSTSUBSCRIPT ( italic_r ) and 𝒞(r)𝒞𝑟\mathcal{C}(r)caligraphic_C ( italic_r ) are provided in SM. In the Mott-Hubbard regime, we have Jss=2Jsd=4Jdd=4t2/Usubscriptsuperscript𝐽𝑠𝑠perpendicular-to2subscriptsuperscript𝐽𝑠𝑑perpendicular-to4subscriptsuperscript𝐽𝑑𝑑perpendicular-to4superscriptsubscript𝑡perpendicular-to2𝑈J^{ss}_{\perp}=2J^{sd}_{\perp}=4J^{dd}_{\perp}=4t_{\perp}^{2}/Uitalic_J start_POSTSUPERSCRIPT italic_s italic_s end_POSTSUPERSCRIPT start_POSTSUBSCRIPT ⟂ end_POSTSUBSCRIPT = 2 italic_J start_POSTSUPERSCRIPT italic_s italic_d end_POSTSUPERSCRIPT start_POSTSUBSCRIPT ⟂ end_POSTSUBSCRIPT = 4 italic_J start_POSTSUPERSCRIPT italic_d italic_d end_POSTSUPERSCRIPT start_POSTSUBSCRIPT ⟂ end_POSTSUBSCRIPT = 4 italic_t start_POSTSUBSCRIPT ⟂ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT / italic_U, and Jss=Jsd=0subscriptsuperscript𝐽𝑠𝑠parallel-tosubscriptsuperscript𝐽𝑠𝑑parallel-to0J^{ss}_{\parallel}=J^{sd}_{\parallel}=0italic_J start_POSTSUPERSCRIPT italic_s italic_s end_POSTSUPERSCRIPT start_POSTSUBSCRIPT ∥ end_POSTSUBSCRIPT = italic_J start_POSTSUPERSCRIPT italic_s italic_d end_POSTSUPERSCRIPT start_POSTSUBSCRIPT ∥ end_POSTSUBSCRIPT = 0, Jdd=t2/Usubscriptsuperscript𝐽𝑑𝑑parallel-tosuperscriptsubscript𝑡parallel-to2𝑈J^{dd}_{\parallel}=t_{\parallel}^{2}/Uitalic_J start_POSTSUPERSCRIPT italic_d italic_d end_POSTSUPERSCRIPT start_POSTSUBSCRIPT ∥ end_POSTSUBSCRIPT = italic_t start_POSTSUBSCRIPT ∥ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT / italic_U Yang et al. (2023d). We stress that the differences in their ratios do not alter the qualitative behavior significantly.

We highlight that even in the charge transfer regime, the type-II t-J model still shows suppressed hopping between the in-plane p1subscript𝑝1p_{1}italic_p start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT orbitals of the two layers, with t0similar-to-or-equalssubscript𝑡perpendicular-to0t_{\perp}\simeq 0italic_t start_POSTSUBSCRIPT ⟂ end_POSTSUBSCRIPT ≃ 0, akin to the behavior observed in the Mott-Hubbard regime. This can be validated by a simple symmetry argument that since p1subscript𝑝1p_{1}italic_p start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT and pzsubscript𝑝𝑧p_{z}italic_p start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT orbitals have different symmetries, there is no direct coupling channel between them at the lowest order. A finite tsubscript𝑡perpendicular-tot_{\perp}italic_t start_POSTSUBSCRIPT ⟂ end_POSTSUBSCRIPT can only be generated through virtual hopping to a pzsubscript𝑝𝑧p_{z}italic_p start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT orbital at a different site and its value should be small. In the following we will show that a small tsubscript𝑡perpendicular-tot_{\perp}italic_t start_POSTSUBSCRIPT ⟂ end_POSTSUBSCRIPT has no significant effect and thus as a good approximation we can set it to be zero.

Refer to caption
Figure 3: DMRG simulation results of the type-II t-J model, Eq.(Type II t-J model in charge transfer regime in bilayer La3Ni2O7 and trilayer La4Ni3O10) in the two-leg ladder at t=1subscript𝑡parallel-to1t_{\parallel}=1italic_t start_POSTSUBSCRIPT ∥ end_POSTSUBSCRIPT = 1, Jss=0.1subscriptsuperscript𝐽parallel-to𝑠𝑠0.1J^{\parallel}_{ss}=0.1italic_J start_POSTSUPERSCRIPT ∥ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_s italic_s end_POSTSUBSCRIPT = 0.1 Here, we impose Jss=43Jsd=169Jddsubscriptsuperscript𝐽𝑠𝑠perpendicular-to43subscriptsuperscript𝐽𝑠𝑑perpendicular-to169subscriptsuperscript𝐽𝑑𝑑perpendicular-toJ^{ss}_{\perp}=\frac{4}{3}J^{sd}_{\perp}=\frac{16}{9}J^{dd}_{\perp}italic_J start_POSTSUPERSCRIPT italic_s italic_s end_POSTSUPERSCRIPT start_POSTSUBSCRIPT ⟂ end_POSTSUBSCRIPT = divide start_ARG 4 end_ARG start_ARG 3 end_ARG italic_J start_POSTSUPERSCRIPT italic_s italic_d end_POSTSUPERSCRIPT start_POSTSUBSCRIPT ⟂ end_POSTSUBSCRIPT = divide start_ARG 16 end_ARG start_ARG 9 end_ARG italic_J start_POSTSUPERSCRIPT italic_d italic_d end_POSTSUPERSCRIPT start_POSTSUBSCRIPT ⟂ end_POSTSUBSCRIPT and Jss=43𝒜Jsd=169𝒜2Jddsubscriptsuperscript𝐽𝑠𝑠parallel-to43𝒜subscriptsuperscript𝐽𝑠𝑑parallel-to169superscript𝒜2subscriptsuperscript𝐽𝑑𝑑parallel-toJ^{ss}_{\parallel}=\frac{4}{3}\mathcal{A}J^{sd}_{\parallel}=\frac{16}{9}% \mathcal{A}^{2}J^{dd}_{\parallel}italic_J start_POSTSUPERSCRIPT italic_s italic_s end_POSTSUPERSCRIPT start_POSTSUBSCRIPT ∥ end_POSTSUBSCRIPT = divide start_ARG 4 end_ARG start_ARG 3 end_ARG caligraphic_A italic_J start_POSTSUPERSCRIPT italic_s italic_d end_POSTSUPERSCRIPT start_POSTSUBSCRIPT ∥ end_POSTSUBSCRIPT = divide start_ARG 16 end_ARG start_ARG 9 end_ARG caligraphic_A start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_J start_POSTSUPERSCRIPT italic_d italic_d end_POSTSUPERSCRIPT start_POSTSUBSCRIPT ∥ end_POSTSUBSCRIPT with 𝒜=0.7705𝒜0.7705\mathcal{A}=0.7705caligraphic_A = 0.7705. (a) The tsubscript𝑡perpendicular-tot_{\perp}italic_t start_POSTSUBSCRIPT ⟂ end_POSTSUBSCRIPT dependence of the spin gap, Δs=E(Sz=1)E(Sz=0)subscriptΔ𝑠𝐸subscript𝑆𝑧1𝐸subscript𝑆𝑧0\Delta_{s}=E(S_{z}=1)-E(S_{z}=0)roman_Δ start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT = italic_E ( italic_S start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT = 1 ) - italic_E ( italic_S start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT = 0 ), of the type-II t-J model at V=0𝑉0V=0italic_V = 0 and x=0.5𝑥0.5x=0.5italic_x = 0.5 with various Jss=1,1.5,2,3,3.5subscriptsuperscript𝐽perpendicular-to𝑠𝑠11.5233.5J^{\perp}_{ss}=1,1.5,2,3,3.5italic_J start_POSTSUPERSCRIPT ⟂ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_s italic_s end_POSTSUBSCRIPT = 1 , 1.5 , 2 , 3 , 3.5. Here, we use Lx=40subscript𝐿𝑥40L_{x}=40italic_L start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT = 40 and χ=2400𝜒2400\chi=2400italic_χ = 2400 for the simulation. The spin gap monotonically decreases as increasing in the tsubscript𝑡perpendicular-tot_{\perp}italic_t start_POSTSUBSCRIPT ⟂ end_POSTSUBSCRIPT in the large Jsubscript𝐽perpendicular-toJ_{\perp}italic_J start_POSTSUBSCRIPT ⟂ end_POSTSUBSCRIPT limit. (b) The doping dependence of the spin gap at t=0.1subscript𝑡perpendicular-to0.1t_{\perp}=0.1italic_t start_POSTSUBSCRIPT ⟂ end_POSTSUBSCRIPT = 0.1, Jss=1superscriptsubscript𝐽perpendicular-to𝑠𝑠1J_{\perp}^{ss}=1italic_J start_POSTSUBSCRIPT ⟂ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_s italic_s end_POSTSUPERSCRIPT = 1 with various V=0,1,2,2.5𝑉0122.5V=0,1,2,2.5italic_V = 0 , 1 , 2 , 2.5. The superconducting dome is exhibited in the BCS limit, where we increase V𝑉Vitalic_V, and the substantial spin gap remains even in x=0.5𝑥0.5x=0.5italic_x = 0.5. (c) The pair correlation at t=0.1subscript𝑡perpendicular-to0.1t_{\perp}=0.1italic_t start_POSTSUBSCRIPT ⟂ end_POSTSUBSCRIPT = 0.1 and x=0.5𝑥0.5x=0.5italic_x = 0.5 with various Jsssuperscriptsubscript𝐽perpendicular-to𝑠𝑠J_{\perp}^{ss}italic_J start_POSTSUBSCRIPT ⟂ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_s italic_s end_POSTSUPERSCRIPT values. It exhibits the power law decaying where the fitted power α𝛼\alphaitalic_α is denoted as the solid line. (d) The scaling of the entanglement entropy and the correlation length ξ𝜉\xiitalic_ξ at x=0.5𝑥0.5x=0.5italic_x = 0.5 within the infinite DMRG calculation. Here, we choose t=0,0.2,0.4subscript𝑡perpendicular-to00.20.4t_{\perp}=0,0.2,0.4italic_t start_POSTSUBSCRIPT ⟂ end_POSTSUBSCRIPT = 0 , 0.2 , 0.4 and J=2,3subscript𝐽perpendicular-to23J_{\perp}=2,3italic_J start_POSTSUBSCRIPT ⟂ end_POSTSUBSCRIPT = 2 , 3. The extracted central charge c𝑐citalic_c is carried by the relation S=c6logξ𝑆𝑐6𝜉S=\frac{c}{6}\log\xiitalic_S = divide start_ARG italic_c end_ARG start_ARG 6 end_ARG roman_log italic_ξ.

Numerical simulations.— We perform the density matrix renormalization group (DMRG) simulations applied to the type-II t-J model described by Eq.(Type II t-J model in charge transfer regime in bilayer La3Ni2O7 and trilayer La4Ni3O10) White (1992); Hauschild and Pollmann (2018). We consider the two-leg ladder configuration (Lz=2,Ly=1,Lxformulae-sequencesubscript𝐿𝑧2formulae-sequencesubscript𝐿𝑦1subscript𝐿𝑥L_{z}=2,L_{y}=1,L_{x}\rightarrow\inftyitalic_L start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT = 2 , italic_L start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT = 1 , italic_L start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT → ∞), rather than the full bilayer two-dimensional cubic lattice, due to the well-known limitation of the DMRG. In our simulations, we set t=1subscript𝑡parallel-to1t_{\parallel}=1italic_t start_POSTSUBSCRIPT ∥ end_POSTSUBSCRIPT = 1 and Jss=0.1superscriptsubscript𝐽parallel-to𝑠𝑠0.1J_{\parallel}^{ss}=0.1italic_J start_POSTSUBSCRIPT ∥ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_s italic_s end_POSTSUPERSCRIPT = 0.1. In Fig. 3, we present the spin gap varying the parameters Jsssuperscriptsubscript𝐽perpendicular-to𝑠𝑠J_{\perp}^{ss}italic_J start_POSTSUBSCRIPT ⟂ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_s italic_s end_POSTSUPERSCRIPT, tsubscript𝑡perpendicular-tot_{\perp}italic_t start_POSTSUBSCRIPT ⟂ end_POSTSUBSCRIPT, and the doping ratio. The convergence of the DMRG results under the system size is also checked an provided in the SM Sec.IV. The presence of a finite spin gap indicates the emergence of the Luther-Emily liquid phase Luther and Emery (1974) with power-law inter-layer pairing correlations. The pairing gap gradually decreases with tsubscript𝑡perpendicular-tot_{\perp}italic_t start_POSTSUBSCRIPT ⟂ end_POSTSUBSCRIPT because inter-layer pairing frustrates the tsubscript𝑡perpendicular-tot_{\perp}italic_t start_POSTSUBSCRIPT ⟂ end_POSTSUBSCRIPT term. But the pairing remains robust at small tsubscript𝑡perpendicular-tot_{\perp}italic_t start_POSTSUBSCRIPT ⟂ end_POSTSUBSCRIPT. In the bilayer nickelate we believe ttmuch-less-thansubscript𝑡perpendicular-to𝑡t_{\perp}\ll titalic_t start_POSTSUBSCRIPT ⟂ end_POSTSUBSCRIPT ≪ italic_t and hence its effect should be negligible. One of the remarkable features of the bilayer type-II t-J model is the doping dependence of the pairing gap. As presented in Fig. 3 (b), the pairing scale shows a dome centered near x=0.40.5𝑥0.40.5x=0.4-0.5italic_x = 0.4 - 0.5 in the presence of a repulsion. This is due to the doping induced BCS to BEC crossoverYang et al. (2023d). Finally, we check the key characteristics of Luther-Emily liquids, such as power-law pair correlation functions and a central charge c=1𝑐1c=1italic_c = 1, in Figs. 3 (c-d). In the main figure, we only illustrate the case with x=0.5𝑥0.5x=0.5italic_x = 0.5 relevant to the experiments, but the Luther-Emily liquid phases are manifested in the broad range of x𝑥xitalic_x (See SM Sec.IV).

Discussion.— Our theoretical proposal can provide a potential scenario on the dome of the Tcsubscript𝑇𝑐T_{c}italic_T start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT versus pressure in the experiment Li et al. (2024b). Increasing pressure should enhance tdp;zsubscript𝑡𝑑𝑝𝑧t_{dp;z}italic_t start_POSTSUBSCRIPT italic_d italic_p ; italic_z end_POSTSUBSCRIPT and thus the Jsubscript𝐽perpendicular-toJ_{\perp}italic_J start_POSTSUBSCRIPT ⟂ end_POSTSUBSCRIPT term. Initially, Tcsubscript𝑇𝑐T_{c}italic_T start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT increasing with the pressure due to the increase of the Jsubscript𝐽perpendicular-toJ_{\perp}italic_J start_POSTSUBSCRIPT ⟂ end_POSTSUBSCRIPT. However, when the pressure is larger than an optimal value Psubscript𝑃P_{*}italic_P start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT, the hole prefers to stay in the pzsubscript𝑝𝑧p_{z}italic_p start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT orbital and then the pairing is suppressed (see SM Sec.V). The shift from the in-plane p1subscript𝑝1p_{1}italic_p start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT orbital to the pzsubscript𝑝𝑧p_{z}italic_p start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT orbital is likely the origin of the dome with pressure.

Our analysis can also be generalized to the trilayer nickelates, La4Ni3O10 Zhu et al. (2024) (see SM Sec.VI) and a trilayer version of the type II t-J model is the minimal model. Now the model is very different from the models proposed assuming Mott-Hubbard regime. For the trilayer case both d1,d2subscript𝑑1subscript𝑑2d_{1},d_{2}italic_d start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , italic_d start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT orbitals are assumed to be at fractional filling and the mobile carriers are argued to be from both orbitals Sakakibara et al. (2024); Zhang et al. (2024b); Lu et al. (2024b). In contrast, in the charge transfer picture, both d1,d2subscript𝑑1subscript𝑑2d_{1},d_{2}italic_d start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , italic_d start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT are still localized which just provide a spin-one moment at each Ni site. Then the doped hole enters the p1subscript𝑝1p_{1}italic_p start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT orbital and the final model is still one-orbital like, similar to the bilayer case. We leave to future to analyze this trilayer model.

Conclusion.— In summary, we provide analytical and numerical study for bilayer nickelates, La3Ni2O7 within the charge transfer framework. Our primary discovery is the identification of the Zhang-Rice spin-half state as the dominant hole state. We emphasize that the type II t-J model again serves as a minimal model for both the bilayer and the trilayer nickelate. The physics is distinct from the hole doped cuprates due to the importance of strong inter-layer spin-spin coupling. This leads to an optimal doping as large as 4050%40percent5040-50\%40 - 50 % in the bilayer type II t-J model in contrast to 20%percent2020\%20 % in cuprate.

Acknowledgement. H.Oh thanks Hui Yang for valuable help on the density matrix renormalization group simulation codes. B.Zhou thanks Xinlong Liu for fruitful discussions on constructing Zhang-Rice spin-half states in the stage of starting the project. This work was supported by the National Science Foundation under Grant No. DMR-2237031.

References

  • Sun et al. (2023) H. Sun, M. Huo, X. Hu, J. Li, Z. Liu, Y. Han, L. Tang, Z. Mao, P. Yang, B. Wang, et al., Nature , 1 (2023).
  • Hou et al. (2023) J. Hou, P. T. Yang, Z. Y. Liu, J. Y. Li, P. F. Shan, L. Ma, G. Wang, N. N. Wang, H. Z. Guo, J. P. Sun, Y. Uwatoko, M. Wang, G. M. Zhang, B. S. Wang,  and J. G. Cheng, “Emergence of high-temperature superconducting phase in the pressurized la3ni2o7 crystals,”  (2023), arXiv:2307.09865 [cond-mat.supr-con] .
  • Zhang et al. (2023a) Y. Zhang, D. Su, Y. Huang, H. Sun, M. Huo, Z. Shan, K. Ye, Z. Yang, R. Li, M. Smidman, M. Wang, L. Jiao,  and H. Yuan, “High-temperature superconductivity with zero-resistance and strange metal behavior in la3ni2o7,”  (2023a), arXiv:2307.14819 [cond-mat.supr-con] .
  • Liu et al. (2023a) Z. Liu, M. Huo, J. Li, Q. Li, Y. Liu, Y. Dai, X. Zhou, J. Hao, Y. Lu, M. Wang, et al., arXiv preprint arXiv:2307.02950  (2023a).
  • Yang et al. (2023a) J. Yang, H. Sun, X. Hu, Y. Xie, T. Miao, H. Luo, H. Chen, B. Liang, W. Zhu, G. Qu, et al., arXiv preprint arXiv:2309.01148  (2023a).
  • Zhang et al. (2023b) M. Zhang, C. Pei, Q. Wang, Y. Zhao, C. Li, W. Cao, S. Zhu, J. Wu,  and Y. Qi, arXiv preprint arXiv:2309.01651  (2023b).
  • Luo et al. (2023) Z. Luo, X. Hu, M. Wang, W. Wú,  and D.-X. Yao, Phys. Rev. Lett. 131, 126001 (2023).
  • Zhang et al. (2023c) Y. Zhang, L.-F. Lin, A. Moreo,  and E. Dagotto, Physical Review B 108 (2023c), 10.1103/physrevb.108.l180510.
  • Yang et al. (2023b) Q.-G. Yang, H.-Y. Liu, D. Wang,  and Q.-H. Wang, arXiv preprint arXiv:2306.03706  (2023b).
  • Sakakibara et al. (2023) H. Sakakibara, N. Kitamine, M. Ochi,  and K. Kuroki, arXiv preprint arXiv:2306.06039  (2023).
  • Gu et al. (2023) Y. Gu, C. Le, Z. Yang, X. Wu,  and J. Hu, arXiv preprint arXiv:2306.07275  (2023).
  • Shen et al. (2023) Y. Shen, M. Qin,  and G.-M. Zhang, arXiv preprint arXiv:2306.07837  (2023).
  • Wú et al. (2023) W. Wú, Z. Luo, D.-X. Yao,  and M. Wang, arXiv preprint arXiv:2307.05662  (2023).
  • Christiansson et al. (2023) V. Christiansson, F. Petocchi,  and P. Werner, Phys. Rev. Lett. 131, 206501 (2023).
  • Liu et al. (2023b) Y.-B. Liu, J.-W. Mei, F. Ye, W.-Q. Chen,  and F. Yang, arXiv preprint arXiv:2307.10144  (2023b).
  • Cao and Yang (2023) Y. Cao and Y.-f. Yang, arXiv preprint arXiv:2307.06806  (2023).
  • Oh and Zhang (2023) H. Oh and Y.-H. Zhang, Phys. Rev. B 108, 174511 (2023).
  • Lu et al. (2024a) C. Lu, Z. Pan, F. Yang,  and C. Wu, Phys. Rev. Lett. 132, 146002 (2024a).
  • Qu et al. (2023) X.-Z. Qu, D.-W. Qu, J. Chen, C. Wu, F. Yang, W. Li,  and G. Su, arXiv preprint arXiv:2307.16873  (2023).
  • Lu et al. (2023) D.-C. Lu, M. Li, Z.-Y. Zeng, W. Hou, J. Wang, F. Yang,  and Y.-Z. You, arXiv preprint arXiv:2308.11195  (2023).
  • Jiang et al. (2023a) R. Jiang, J. Hou, Z. Fan, Z.-J. Lang,  and W. Ku, arXiv preprint arXiv:2308.11614  (2023a).
  • Tian et al. (2023) Y.-H. Tian, Y. Chen, J.-M. Wang, R.-Q. He,  and Z.-Y. Lu, arXiv preprint arXiv:2308.09698  (2023).
  • Zhang et al. (2023d) J.-X. Zhang, H.-K. Zhang, Y.-Z. You,  and Z.-Y. Weng, arXiv preprint arXiv:2309.05726  (2023d).
  • Qin and Yang (2023) Q. Qin and Y.-f. Yang, arXiv preprint arXiv:2308.09044  (2023).
  • Huang et al. (2023) J. Huang, Z. Wang,  and T. Zhou, arXiv preprint arXiv:2308.07651  (2023).
  • Zhang et al. (2023e) Y. Zhang, L.-F. Lin, A. Moreo, T. A. Maier,  and E. Dagotto, arXiv preprint arXiv:2308.07386  (2023e).
  • Jiang et al. (2023b) K. Jiang, Z. Wang,  and F.-C. Zhang, arXiv preprint arXiv:2308.06771  (2023b).
  • Yang et al. (2023c) Y.-f. Yang, G.-M. Zhang,  and F.-C. Zhang, arXiv preprint arXiv:2308.01176  (2023c).
  • Qin and Yang (2023) Q. Qin and Y.-f. Yang, arXiv e-prints , arXiv:2308.09044 (2023)arXiv:2308.09044 [cond-mat.supr-con] .
  • Zhang et al. (2023) Y. Zhang, L.-F. Lin, A. Moreo, T. A. Maier,  and E. Dagotto, arXiv e-prints , arXiv:2308.07386 (2023)arXiv:2308.07386 [cond-mat.supr-con] .
  • Kitamine et al. (2023) N. Kitamine, M. Ochi,  and K. Kuroki, arXiv e-prints , arXiv:2308.12750 (2023)arXiv:2308.12750 [cond-mat.supr-con] .
  • Jiang et al. (2023) R. Jiang, J. Hou, Z. Fan, Z.-J. Lang,  and W. Ku, arXiv e-prints , arXiv:2308.11614 (2023)arXiv:2308.11614 [cond-mat.supr-con] .
  • Lange et al. (2023a) H. Lange, L. Homeier, E. Demler, U. Schollwöck, A. Bohrdt,  and F. Grusdt, “Pairing dome from an emergent feshbach resonance in a strongly repulsive bilayer model,”  (2023a), arXiv:2309.13040 [cond-mat.str-el] .
  • Lange et al. (2023b) H. Lange, L. Homeier, E. Demler, U. Schollwöck, F. Grusdt,  and A. Bohrdt, “Feshbach resonance in a strongly repulsive bilayer model: a possible scenario for bilayer nickelate superconductors,”  (2023b), arXiv:2309.15843 [cond-mat.str-el] .
  • Schlömer et al. (2023) H. Schlömer, U. Schollwöck, F. Grusdt,  and A. Bohrdt, “Superconductivity in the pressurized nickelate la3ni2o7 in the vicinity of a bec-bcs crossover,”  (2023), arXiv:2311.03349 [cond-mat.str-el] .
  • Jiang et al. (2024a) R. Jiang, J. Hou, Z. Fan, Z.-J. Lang, W. Ku, et al., Physical Review Letters 132, 126503 (2024a).
  • Zhan et al. (2024) J. Zhan, Y. Gu, X. Wu,  and J. Hu, “Cooperation between electron-phonon coupling and electronic interaction in bilayer nickelates la3ni2o7,”  (2024), arXiv:2404.03638 [cond-mat.supr-con] .
  • Lechermann et al. (2024) F. Lechermann, S. Bötzel,  and I. M. Eremin, arXiv preprint arXiv:2403.12831  (2024).
  • Wang et al. (2024) G. Wang, N. Wang, X. Shen, J. Hou, L. Ma, L. Shi, Z. Ren, Y. Gu, H. Ma, P. Yang, et al., Physical Review X 14, 011040 (2024).
  • Xue and Wang (2024) J.-R. Xue and F. Wang, arXiv preprint arXiv:2402.07449  (2024).
  • Kaneko et al. (2024) T. Kaneko, H. Sakakibara, M. Ochi,  and K. Kuroki, Physical Review B 109, 045154 (2024).
  • Jiang et al. (2024b) K. Jiang, Z. Wang,  and F.-C. Zhang, Chinese Physics Letters 41, 017402 (2024b).
  • Ryee et al. (2023) S. Ryee, N. Witt,  and T. O. Wehling, arXiv preprint arXiv:2310.17465  (2023).
  • Liao et al. (2023) Z. Liao, L. Chen, G. Duan, Y. Wang, C. Liu, R. Yu,  and Q. Si, Phys. Rev. B 108, 214522 (2023).
  • Li et al. (2024a) Q. Li, Z. Ying-Jie, Z.-N. Xiang, Y. Zhang, X. Zhu,  and H.-H. Wen, Chinese Physics Letters 41, 017401 (2024a).
  • Zhu et al. (2024) Y. Zhu, E. Zhang, B. Pan, X. Chen, D. Peng, L. Chen, H. Ren, F. Liu, N. Li, Z. Xing, J. Han, J. Wang, D. Jia, H. Wo, Y. Gu, Y. Gu, L. Ji, W. Wang, H. Gou, Y. Shen, T. Ying, X. Chen, W. Yang, C. Zheng, Q. Zeng, J. gang Guo,  and J. Zhao, “Superconductivity in trilayer nickelate la4ni3o10 single crystals,”  (2024), arXiv:2311.07353 [cond-mat.supr-con] .
  • Sakakibara et al. (2024) H. Sakakibara, M. Ochi, H. Nagata, Y. Ueki, H. Sakurai, R. Matsumoto, K. Terashima, K. Hirose, H. Ohta, M. Kato, et al., Physical Review B 109, 144511 (2024).
  • Zhang et al. (2024a) M. Zhang, C. Pei, X. Du, W. Hu, Y. Cao, Q. Wang, J. Wu, Y. Li, H. Liu, C. Wen, Y. Zhao, C. Li, W. Cao, S. Zhu, Q. Zhang, N. Yu, P. Cheng, L. Zhang, Z. Li, J. Zhao, Y. Chen, H. Guo, C. Wu, F. Yang, S. Yan, L. Yang,  and Y. Qi, “Superconductivity in trilayer nickelate la4ni3o10 under pressure,”  (2024a), arXiv:2311.07423 [cond-mat.supr-con] .
  • Lu et al. (2024b) C. Lu, Z. Pan, F. Yang,  and C. Wu, arXiv preprint arXiv:2402.06450  (2024b).
  • Zhang et al. (2024b) Y. Zhang, L.-F. Lin, A. Moreo, T. A. Maier,  and E. Dagotto, arXiv preprint arXiv:2402.05285  (2024b).
  • LaBollita et al. (2024) H. LaBollita, J. Kapeghian, M. R. Norman,  and A. S. Botana, arXiv preprint arXiv:2402.05085  (2024).
  • Yang et al. (2024) Q.-G. Yang, K.-Y. Jiang, D. Wang, H.-Y. Lu,  and Q.-H. Wang, arXiv preprint arXiv:2402.05447  (2024).
  • Yang et al. (2023d) H. Yang, H. Oh,  and Y.-H. Zhang, “Strong pairing from small fermi surface beyond weak coupling: Application to la3ni2o7,”  (2023d), arXiv:2309.15095 [cond-mat.str-el] .
  • Wu et al. (2024) X. Wu, H. Yang,  and Y.-H. Zhang, “Deconfined fermi liquid to fermi liquid transition and superconducting instability,”  (2024), arXiv:2401.08753 [cond-mat.str-el] .
  • Zaanen et al. (1985) J. Zaanen, G. A. Sawatzky,  and J. W. Allen, Phys. Rev. Lett. 55, 418 (1985).
  • Dong et al. (2023) Z. Dong, M. Huo, J. Li, J. Li, P. Li, H. Sun, Y. Lu, M. Wang, Y. Wang,  and Z. Chen, arXiv preprint arXiv:2312.15727  (2023).
  • Zhang and Rice (1988) F. Zhang and T. Rice, Physical Review B 37, 3759 (1988).
  • Zhang and Vishwanath (2020) Y.-H. Zhang and A. Vishwanath, Physical Review Research 2, 023112 (2020).
  • Ogata and Fukuyama (2008) M. Ogata and H. Fukuyama, Reports on Progress in Physics 71, 036501 (2008).
  • White (1992) S. R. White, Phys. Rev. Lett. 69, 2863 (1992).
  • Hauschild and Pollmann (2018) J. Hauschild and F. Pollmann, SciPost Phys. Lect. Notes , 5 (2018).
  • Luther and Emery (1974) A. Luther and V. J. Emery, Phys. Rev. Lett. 33, 589 (1974).
  • Li et al. (2024b) J. Li, P. Ma, H. Zhang, X. Huang, C. Huang, M. Huo, D. Hu, Z. Dong, C. He, J. Liao, X. Chen, T. Xie, H. Sun,  and M. Wang, “Pressure-driven right-triangle shape superconductivity in bilayer nickelate la3ni2o7,”  (2024b), arXiv:2404.11369 [cond-mat.supr-con] .
  • Luo et al. (2024) Z. Luo, C.-Q. Chen, M. Wang, W. Wú,  and D.-X. Yao, “Trilayer multi-orbital models of La4Ni3O10subscriptLa4subscriptNi3subscriptO10\mathrm{La_{4}Ni_{3}O_{10}}roman_La start_POSTSUBSCRIPT 4 end_POSTSUBSCRIPT roman_Ni start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT roman_O start_POSTSUBSCRIPT 10 end_POSTSUBSCRIPT,”  (2024), arXiv:2402.07196 [cond-mat.supr-con] .

Supplemental Material for “Type II t-J model in charge transfer regime in bilayer La3Ni2O7 and trilayer La4Ni3O10

Hanbit Oh{}^{\ {\color[rgb]{1,0,0}\definecolor[named]{pgfstrokecolor}{rgb}{1,0,0}*}}start_FLOATSUPERSCRIPT ∗ end_FLOATSUPERSCRIPT, Boran Zhou, and Ya-Hui Zhang{}^{\ {\color[rgb]{1,0,0}\definecolor[named]{pgfstrokecolor}{rgb}{1,0,0}% \dagger}}start_FLOATSUPERSCRIPT † end_FLOATSUPERSCRIPT

William H. Miller III Department of Physics and Astronomy,
Johns Hopkins University, Baltimore, Maryland, 21218, USA

I I. Symmetry properties of the Charge transfer Hamiltonian

The charge transfer Hamiltonian enjoys D4h=D4subscript𝐷4tensor-productsubscript𝐷4D_{4h}=D_{4}\otimes\mathcal{I}italic_D start_POSTSUBSCRIPT 4 italic_h end_POSTSUBSCRIPT = italic_D start_POSTSUBSCRIPT 4 end_POSTSUBSCRIPT ⊗ caligraphic_I point group symmetry defined at Ni i𝑖iitalic_i site. While the two d𝑑ditalic_d orbitals, (d1,d2subscript𝑑1subscript𝑑2d_{1},d_{2}italic_d start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , italic_d start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT) are already one-dimensional irreducible representations of D4hsubscript𝐷4D_{4h}italic_D start_POSTSUBSCRIPT 4 italic_h end_POSTSUBSCRIPT, px,ysubscript𝑝𝑥𝑦p_{x,y}italic_p start_POSTSUBSCRIPT italic_x , italic_y end_POSTSUBSCRIPT orbitals living at the link of the cubic lattice are not. To classify them as an symmetry representation, we first identify the symmetry actions of the p𝑝pitalic_p orbitals. For example, under the D4={E,2C4,C2,2C2,2C2′′}subscript𝐷4𝐸2subscript𝐶4subscript𝐶22superscriptsubscript𝐶22superscriptsubscript𝐶2′′D_{4}=\{E,2C_{4},C_{2},2C_{2}^{\prime},2C_{2}^{\prime\prime}\}italic_D start_POSTSUBSCRIPT 4 end_POSTSUBSCRIPT = { italic_E , 2 italic_C start_POSTSUBSCRIPT 4 end_POSTSUBSCRIPT , italic_C start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT , 2 italic_C start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT , 2 italic_C start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ ′ end_POSTSUPERSCRIPT }, one can show that

C4:px;i±x2;lpy;i±y2;l,py;i±y2;lpx;ix2;l,\displaystyle C_{4}:\quad p_{x;i\pm\frac{x}{2};l}\rightarrow p_{y;i\pm\frac{y}% {2};l},\quad p_{y;i\pm\frac{y}{2};l}\rightarrow-p_{x;i\mp\frac{x}{2};l},italic_C start_POSTSUBSCRIPT 4 end_POSTSUBSCRIPT : italic_p start_POSTSUBSCRIPT italic_x ; italic_i ± divide start_ARG italic_x end_ARG start_ARG 2 end_ARG ; italic_l end_POSTSUBSCRIPT → italic_p start_POSTSUBSCRIPT italic_y ; italic_i ± divide start_ARG italic_y end_ARG start_ARG 2 end_ARG ; italic_l end_POSTSUBSCRIPT , italic_p start_POSTSUBSCRIPT italic_y ; italic_i ± divide start_ARG italic_y end_ARG start_ARG 2 end_ARG ; italic_l end_POSTSUBSCRIPT → - italic_p start_POSTSUBSCRIPT italic_x ; italic_i ∓ divide start_ARG italic_x end_ARG start_ARG 2 end_ARG ; italic_l end_POSTSUBSCRIPT ,
C2:px;i±x2;lpx;ix2;l,py;i±y2;lpy;iy2;l,\displaystyle C_{2}:\quad p_{x;i\pm\frac{x}{2};l}\rightarrow-p_{x;i\mp\frac{x}% {2};l},\quad p_{y;i\pm\frac{y}{2};l}\rightarrow-p_{y;i\mp\frac{y}{2};l},italic_C start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT : italic_p start_POSTSUBSCRIPT italic_x ; italic_i ± divide start_ARG italic_x end_ARG start_ARG 2 end_ARG ; italic_l end_POSTSUBSCRIPT → - italic_p start_POSTSUBSCRIPT italic_x ; italic_i ∓ divide start_ARG italic_x end_ARG start_ARG 2 end_ARG ; italic_l end_POSTSUBSCRIPT , italic_p start_POSTSUBSCRIPT italic_y ; italic_i ± divide start_ARG italic_y end_ARG start_ARG 2 end_ARG ; italic_l end_POSTSUBSCRIPT → - italic_p start_POSTSUBSCRIPT italic_y ; italic_i ∓ divide start_ARG italic_y end_ARG start_ARG 2 end_ARG ; italic_l end_POSTSUBSCRIPT ,
C2:px;i±x2;lpx;i±x2;l,py;i±y2;lpy;iy2;l,\displaystyle C_{2}^{\prime}:\quad p_{x;i\pm\frac{x}{2};l}\rightarrow p_{x;i% \pm\frac{x}{2};l},\quad p_{y;i\pm\frac{y}{2};l}\rightarrow-p_{y;i\mp\frac{y}{2% };l},italic_C start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT : italic_p start_POSTSUBSCRIPT italic_x ; italic_i ± divide start_ARG italic_x end_ARG start_ARG 2 end_ARG ; italic_l end_POSTSUBSCRIPT → italic_p start_POSTSUBSCRIPT italic_x ; italic_i ± divide start_ARG italic_x end_ARG start_ARG 2 end_ARG ; italic_l end_POSTSUBSCRIPT , italic_p start_POSTSUBSCRIPT italic_y ; italic_i ± divide start_ARG italic_y end_ARG start_ARG 2 end_ARG ; italic_l end_POSTSUBSCRIPT → - italic_p start_POSTSUBSCRIPT italic_y ; italic_i ∓ divide start_ARG italic_y end_ARG start_ARG 2 end_ARG ; italic_l end_POSTSUBSCRIPT ,
C2′′:px;i±x2;lpy;i±y2;l,py;i±y2;lpx;ix2;l,\displaystyle C_{2}^{\prime\prime}:\quad p_{x;i\pm\frac{x}{2};l}\rightarrow p_% {y;i\pm\frac{y}{2};l},\quad p_{y;i\pm\frac{y}{2};l}\rightarrow p_{x;i\mp\frac{% x}{2};l},italic_C start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ ′ end_POSTSUPERSCRIPT : italic_p start_POSTSUBSCRIPT italic_x ; italic_i ± divide start_ARG italic_x end_ARG start_ARG 2 end_ARG ; italic_l end_POSTSUBSCRIPT → italic_p start_POSTSUBSCRIPT italic_y ; italic_i ± divide start_ARG italic_y end_ARG start_ARG 2 end_ARG ; italic_l end_POSTSUBSCRIPT , italic_p start_POSTSUBSCRIPT italic_y ; italic_i ± divide start_ARG italic_y end_ARG start_ARG 2 end_ARG ; italic_l end_POSTSUBSCRIPT → italic_p start_POSTSUBSCRIPT italic_x ; italic_i ∓ divide start_ARG italic_x end_ARG start_ARG 2 end_ARG ; italic_l end_POSTSUBSCRIPT ,

Solving the eigenvalues of the above symmetry actions, we can find that the following linear combinations,

p1;i;l;σsubscript𝑝1𝑖𝑙𝜎\displaystyle p_{1;i;l;\sigma}italic_p start_POSTSUBSCRIPT 1 ; italic_i ; italic_l ; italic_σ end_POSTSUBSCRIPT =\displaystyle== 12[px;i+x^2;l;σpx;ix^2;l;σpy;i+y^2;l;σ+py;iy^2;l;σ],12delimited-[]subscript𝑝𝑥𝑖^𝑥2𝑙𝜎subscript𝑝𝑥𝑖^𝑥2𝑙𝜎subscript𝑝𝑦𝑖^𝑦2𝑙𝜎subscript𝑝𝑦𝑖^𝑦2𝑙𝜎\displaystyle\frac{1}{2}\left[p_{x;i+\frac{{\widehat{x}}}{2};l;\sigma}-p_{x;i-% \frac{{\widehat{x}}}{2};l;\sigma}-p_{y;i+\frac{{\widehat{y}}}{2};l;\sigma}+p_{% y;i-\frac{{\widehat{y}}}{2};l;\sigma}\right],divide start_ARG 1 end_ARG start_ARG 2 end_ARG [ italic_p start_POSTSUBSCRIPT italic_x ; italic_i + divide start_ARG over^ start_ARG italic_x end_ARG end_ARG start_ARG 2 end_ARG ; italic_l ; italic_σ end_POSTSUBSCRIPT - italic_p start_POSTSUBSCRIPT italic_x ; italic_i - divide start_ARG over^ start_ARG italic_x end_ARG end_ARG start_ARG 2 end_ARG ; italic_l ; italic_σ end_POSTSUBSCRIPT - italic_p start_POSTSUBSCRIPT italic_y ; italic_i + divide start_ARG over^ start_ARG italic_y end_ARG end_ARG start_ARG 2 end_ARG ; italic_l ; italic_σ end_POSTSUBSCRIPT + italic_p start_POSTSUBSCRIPT italic_y ; italic_i - divide start_ARG over^ start_ARG italic_y end_ARG end_ARG start_ARG 2 end_ARG ; italic_l ; italic_σ end_POSTSUBSCRIPT ] , (S1)
p2;i;l;σsubscript𝑝2𝑖𝑙𝜎\displaystyle p_{2;i;l;\sigma}italic_p start_POSTSUBSCRIPT 2 ; italic_i ; italic_l ; italic_σ end_POSTSUBSCRIPT =\displaystyle== 12[px;i+x^2;l;σpx;ix^2;l;σ+py;i+y^2;l;σpy;iy^2;l;σ],12delimited-[]subscript𝑝𝑥𝑖^𝑥2𝑙𝜎subscript𝑝𝑥𝑖^𝑥2𝑙𝜎subscript𝑝𝑦𝑖^𝑦2𝑙𝜎subscript𝑝𝑦𝑖^𝑦2𝑙𝜎\displaystyle\frac{1}{2}\left[p_{x;i+\frac{{\widehat{x}}}{2};l;\sigma}-p_{x;i-% \frac{{\widehat{x}}}{2};l;\sigma}+p_{y;i+\frac{{\widehat{y}}}{2};l;\sigma}-p_{% y;i-\frac{{\widehat{y}}}{2};l;\sigma}\right],divide start_ARG 1 end_ARG start_ARG 2 end_ARG [ italic_p start_POSTSUBSCRIPT italic_x ; italic_i + divide start_ARG over^ start_ARG italic_x end_ARG end_ARG start_ARG 2 end_ARG ; italic_l ; italic_σ end_POSTSUBSCRIPT - italic_p start_POSTSUBSCRIPT italic_x ; italic_i - divide start_ARG over^ start_ARG italic_x end_ARG end_ARG start_ARG 2 end_ARG ; italic_l ; italic_σ end_POSTSUBSCRIPT + italic_p start_POSTSUBSCRIPT italic_y ; italic_i + divide start_ARG over^ start_ARG italic_y end_ARG end_ARG start_ARG 2 end_ARG ; italic_l ; italic_σ end_POSTSUBSCRIPT - italic_p start_POSTSUBSCRIPT italic_y ; italic_i - divide start_ARG over^ start_ARG italic_y end_ARG end_ARG start_ARG 2 end_ARG ; italic_l ; italic_σ end_POSTSUBSCRIPT ] , (S2)

become one-dimensional representations with B1gsubscript𝐵1𝑔B_{1g}italic_B start_POSTSUBSCRIPT 1 italic_g end_POSTSUBSCRIPT, A1gsubscript𝐴1𝑔A_{1g}italic_A start_POSTSUBSCRIPT 1 italic_g end_POSTSUBSCRIPT representation. This can be checked by the character table in Table S1. We tabulate the characters of the each orbitals and the irreducible representation. Here, g(u)𝑔𝑢g(u)italic_g ( italic_u ) denotes the even (odd) under the inversion \mathcal{I}caligraphic_I.

The symmetry allowed nearest neighbor hopping between d,p𝑑𝑝d,pitalic_d , italic_p orbitals are given as

Hdpsubscript𝐻𝑑𝑝\displaystyle H_{dp}italic_H start_POSTSUBSCRIPT italic_d italic_p end_POSTSUBSCRIPT =\displaystyle== l,σa;i,α;i[tdpa;i,α;ida;i;l;σpα;i;l;σ+H.c.]+σi[tdp;z(d2;i;t;σd2;i;b;σ)pz;i;σ+H.c.]\displaystyle\sum_{l,\sigma}\sum_{\langle a;i,\alpha;i^{\prime}\rangle}\left[t% _{dp}^{a;i,\alpha;i^{\prime}}d^{\dagger}_{a;i;l;\sigma}p_{\alpha;i^{\prime};l;% \sigma}+\mathrm{H.c.}\right]+\sum_{\sigma}\sum_{i}\left[t_{dp;z}(d^{\dagger}_{% 2;i;t;\sigma}-d^{\dagger}_{2;i;b;\sigma})p_{z;i;\sigma}+\mathrm{H.c.}\right]∑ start_POSTSUBSCRIPT italic_l , italic_σ end_POSTSUBSCRIPT ∑ start_POSTSUBSCRIPT ⟨ italic_a ; italic_i , italic_α ; italic_i start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ⟩ end_POSTSUBSCRIPT [ italic_t start_POSTSUBSCRIPT italic_d italic_p end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_a ; italic_i , italic_α ; italic_i start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUPERSCRIPT italic_d start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_a ; italic_i ; italic_l ; italic_σ end_POSTSUBSCRIPT italic_p start_POSTSUBSCRIPT italic_α ; italic_i start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ; italic_l ; italic_σ end_POSTSUBSCRIPT + roman_H . roman_c . ] + ∑ start_POSTSUBSCRIPT italic_σ end_POSTSUBSCRIPT ∑ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT [ italic_t start_POSTSUBSCRIPT italic_d italic_p ; italic_z end_POSTSUBSCRIPT ( italic_d start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT start_POSTSUBSCRIPT 2 ; italic_i ; italic_t ; italic_σ end_POSTSUBSCRIPT - italic_d start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT start_POSTSUBSCRIPT 2 ; italic_i ; italic_b ; italic_σ end_POSTSUBSCRIPT ) italic_p start_POSTSUBSCRIPT italic_z ; italic_i ; italic_σ end_POSTSUBSCRIPT + roman_H . roman_c . ]

where i𝑖iitalic_i(isuperscript𝑖i^{\prime}italic_i start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT) are sites of nickel (oxygen) atoms. The phase factors should be defined accordingly on the hopping terms,

tdp1;i,α;isuperscriptsubscript𝑡𝑑𝑝1𝑖𝛼superscript𝑖\displaystyle t_{dp}^{1;i,\alpha;i^{\prime}}italic_t start_POSTSUBSCRIPT italic_d italic_p end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 1 ; italic_i , italic_α ; italic_i start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUPERSCRIPT =\displaystyle== tdp;1s={±1}(1)s[δα,xδi,i+sx^2δα,yδi,i+sy^2]subscript𝑡𝑑𝑝1subscript𝑠plus-or-minus1superscript1𝑠delimited-[]subscript𝛿𝛼𝑥subscript𝛿superscript𝑖𝑖𝑠^𝑥2subscript𝛿𝛼𝑦subscript𝛿superscript𝑖𝑖𝑠^𝑦2\displaystyle t_{dp;1}\sum_{s=\{\pm 1\}}(-1)^{s}[\delta_{\alpha,x}\delta_{i^{% \prime},i+s\frac{{\widehat{x}}}{2}}-\delta_{\alpha,y}\delta_{i^{\prime},i+s% \frac{{\widehat{y}}}{2}}]italic_t start_POSTSUBSCRIPT italic_d italic_p ; 1 end_POSTSUBSCRIPT ∑ start_POSTSUBSCRIPT italic_s = { ± 1 } end_POSTSUBSCRIPT ( - 1 ) start_POSTSUPERSCRIPT italic_s end_POSTSUPERSCRIPT [ italic_δ start_POSTSUBSCRIPT italic_α , italic_x end_POSTSUBSCRIPT italic_δ start_POSTSUBSCRIPT italic_i start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT , italic_i + italic_s divide start_ARG over^ start_ARG italic_x end_ARG end_ARG start_ARG 2 end_ARG end_POSTSUBSCRIPT - italic_δ start_POSTSUBSCRIPT italic_α , italic_y end_POSTSUBSCRIPT italic_δ start_POSTSUBSCRIPT italic_i start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT , italic_i + italic_s divide start_ARG over^ start_ARG italic_y end_ARG end_ARG start_ARG 2 end_ARG end_POSTSUBSCRIPT ]
tdp2;i,α;isuperscriptsubscript𝑡𝑑𝑝2𝑖𝛼superscript𝑖\displaystyle t_{dp}^{2;i,\alpha;i^{\prime}}italic_t start_POSTSUBSCRIPT italic_d italic_p end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 ; italic_i , italic_α ; italic_i start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUPERSCRIPT =\displaystyle== tdp;2s={±1}(1)s[δα,xδi,i+sx^2+δα,yδi,i+sy^2].subscript𝑡𝑑𝑝2subscript𝑠plus-or-minus1superscript1𝑠delimited-[]subscript𝛿𝛼𝑥subscript𝛿superscript𝑖𝑖𝑠^𝑥2subscript𝛿𝛼𝑦subscript𝛿superscript𝑖𝑖𝑠^𝑦2\displaystyle t_{dp;2}\sum_{s=\{\pm 1\}}(-1)^{s}[\delta_{\alpha,x}\delta_{i^{% \prime},i+s\frac{{\widehat{x}}}{2}}+\delta_{\alpha,y}\delta_{i^{\prime},i+s% \frac{{\widehat{y}}}{2}}].italic_t start_POSTSUBSCRIPT italic_d italic_p ; 2 end_POSTSUBSCRIPT ∑ start_POSTSUBSCRIPT italic_s = { ± 1 } end_POSTSUBSCRIPT ( - 1 ) start_POSTSUPERSCRIPT italic_s end_POSTSUPERSCRIPT [ italic_δ start_POSTSUBSCRIPT italic_α , italic_x end_POSTSUBSCRIPT italic_δ start_POSTSUBSCRIPT italic_i start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT , italic_i + italic_s divide start_ARG over^ start_ARG italic_x end_ARG end_ARG start_ARG 2 end_ARG end_POSTSUBSCRIPT + italic_δ start_POSTSUBSCRIPT italic_α , italic_y end_POSTSUBSCRIPT italic_δ start_POSTSUBSCRIPT italic_i start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT , italic_i + italic_s divide start_ARG over^ start_ARG italic_y end_ARG end_ARG start_ARG 2 end_ARG end_POSTSUBSCRIPT ] .

Substituting the Eqs.(S1-S2) into Hdpsubscript𝐻𝑑𝑝H_{dp}italic_H start_POSTSUBSCRIPT italic_d italic_p end_POSTSUBSCRIPT simply reduces to,

Hdpsubscript𝐻𝑑𝑝\displaystyle H_{dp}italic_H start_POSTSUBSCRIPT italic_d italic_p end_POSTSUBSCRIPT =\displaystyle== i,l,a,σ[2tdp;ada;i;l;σpa;i;l;σ+H.c.]+i,σ[tdp;z(d2;i;t;σd2;i;b;σ)pz;i;σ+H.c.]\displaystyle\sum_{i,l,a,\sigma}\left[2t_{dp;a}d^{\dagger}_{a;i;l;\sigma}p_{a;% i;l;\sigma}+\mathrm{H.c.}\right]+\sum_{i,\sigma}\left[t_{dp;z}(d^{\dagger}_{2;% i;t;\sigma}-d^{\dagger}_{2;i;b;\sigma})p_{z;i;\sigma}+\mathrm{H.c.}\right]∑ start_POSTSUBSCRIPT italic_i , italic_l , italic_a , italic_σ end_POSTSUBSCRIPT [ 2 italic_t start_POSTSUBSCRIPT italic_d italic_p ; italic_a end_POSTSUBSCRIPT italic_d start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_a ; italic_i ; italic_l ; italic_σ end_POSTSUBSCRIPT italic_p start_POSTSUBSCRIPT italic_a ; italic_i ; italic_l ; italic_σ end_POSTSUBSCRIPT + roman_H . roman_c . ] + ∑ start_POSTSUBSCRIPT italic_i , italic_σ end_POSTSUBSCRIPT [ italic_t start_POSTSUBSCRIPT italic_d italic_p ; italic_z end_POSTSUBSCRIPT ( italic_d start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT start_POSTSUBSCRIPT 2 ; italic_i ; italic_t ; italic_σ end_POSTSUBSCRIPT - italic_d start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT start_POSTSUBSCRIPT 2 ; italic_i ; italic_b ; italic_σ end_POSTSUBSCRIPT ) italic_p start_POSTSUBSCRIPT italic_z ; italic_i ; italic_σ end_POSTSUBSCRIPT + roman_H . roman_c . ]

which is introduced in the main-text.

Orbitals R E𝐸Eitalic_E 2C4(z)2subscript𝐶4𝑧2C_{4}(z)2 italic_C start_POSTSUBSCRIPT 4 end_POSTSUBSCRIPT ( italic_z ) C2subscript𝐶2C_{2}italic_C start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT 2C22superscriptsubscript𝐶22C_{2}^{\prime}2 italic_C start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT 2C2′′2superscriptsubscript𝐶2′′2C_{2}^{\prime\prime}2 italic_C start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ ′ end_POSTSUPERSCRIPT \mathcal{I}caligraphic_I
d1;l,p1;lsubscript𝑑1𝑙subscript𝑝1𝑙d_{1;l},p_{1;l}italic_d start_POSTSUBSCRIPT 1 ; italic_l end_POSTSUBSCRIPT , italic_p start_POSTSUBSCRIPT 1 ; italic_l end_POSTSUBSCRIPT B1gsubscript𝐵1𝑔B_{1g}italic_B start_POSTSUBSCRIPT 1 italic_g end_POSTSUBSCRIPT 1 -1 1 1 -1 1
d2;l,p2;lsubscript𝑑2𝑙subscript𝑝2𝑙d_{2;l},p_{2;l}italic_d start_POSTSUBSCRIPT 2 ; italic_l end_POSTSUBSCRIPT , italic_p start_POSTSUBSCRIPT 2 ; italic_l end_POSTSUBSCRIPT A1gsubscript𝐴1𝑔A_{1g}italic_A start_POSTSUBSCRIPT 1 italic_g end_POSTSUBSCRIPT 1 1 1 1 1 1
d2;td2;b,pzsubscript𝑑2𝑡subscript𝑑2𝑏subscript𝑝𝑧d_{2;t}-d_{2;b},p_{z}italic_d start_POSTSUBSCRIPT 2 ; italic_t end_POSTSUBSCRIPT - italic_d start_POSTSUBSCRIPT 2 ; italic_b end_POSTSUBSCRIPT , italic_p start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT A1usubscript𝐴1𝑢A_{1u}italic_A start_POSTSUBSCRIPT 1 italic_u end_POSTSUBSCRIPT 1 1 1 1 1 -1
Table S1: The irreducible representation (R) and the character table of D4hsubscript𝐷4D_{4h}italic_D start_POSTSUBSCRIPT 4 italic_h end_POSTSUBSCRIPT point group.

II II. Zhang-Rice spin-half state with general JK,JHsubscript𝐽𝐾subscript𝐽𝐻J_{K},J_{H}italic_J start_POSTSUBSCRIPT italic_K end_POSTSUBSCRIPT , italic_J start_POSTSUBSCRIPT italic_H end_POSTSUBSCRIPT

In this section, we provide the energy analysis of possible d8Lsuperscript𝑑8𝐿d^{8}Litalic_d start_POSTSUPERSCRIPT 8 end_POSTSUPERSCRIPT italic_L states with general values of JK,JHsubscript𝐽𝐾subscript𝐽𝐻J_{K},J_{H}italic_J start_POSTSUBSCRIPT italic_K end_POSTSUBSCRIPT , italic_J start_POSTSUBSCRIPT italic_H end_POSTSUBSCRIPT. In the strong coupling limit JK,JH1much-greater-thansubscript𝐽𝐾subscript𝐽𝐻1J_{K},J_{H}\gg 1italic_J start_POSTSUBSCRIPT italic_K end_POSTSUBSCRIPT , italic_J start_POSTSUBSCRIPT italic_H end_POSTSUBSCRIPT ≫ 1, we can consider the local Hamiltonian defined at each site, neglecting the hopping term,

H𝐻\displaystyle Hitalic_H =\displaystyle== l;a2JK;a[sl;adsl;ap14nl;ap]+l2JK;z[sl;2dszp14nzp]l2JH[sl;1dsl;2d+14],subscript𝑙𝑎2subscript𝐽𝐾𝑎delimited-[]subscriptsuperscript𝑠𝑑𝑙𝑎subscriptsuperscript𝑠𝑝𝑙𝑎14subscriptsuperscript𝑛𝑝𝑙𝑎subscript𝑙2subscript𝐽𝐾𝑧delimited-[]subscriptsuperscript𝑠𝑑𝑙2subscriptsuperscript𝑠𝑝𝑧14subscriptsuperscript𝑛𝑝𝑧subscript𝑙2subscript𝐽𝐻delimited-[]subscriptsuperscript𝑠𝑑𝑙1subscriptsuperscript𝑠𝑑𝑙214\displaystyle\sum_{l;a}2J_{K;a}[\vec{s}^{d}_{l;a}\cdot\vec{s}^{p}_{l;a}-\frac{% 1}{4}n^{p}_{l;a}]+\sum_{l}2J_{K;z}[\vec{s}^{d}_{l;2}\cdot\vec{s}^{p}_{z}-\frac% {1}{4}n^{p}_{z}]-\sum_{l}2J_{H}[\vec{s}^{d}_{l;1}\cdot\vec{s}^{d}_{l;2}+\frac{% 1}{4}],∑ start_POSTSUBSCRIPT italic_l ; italic_a end_POSTSUBSCRIPT 2 italic_J start_POSTSUBSCRIPT italic_K ; italic_a end_POSTSUBSCRIPT [ over→ start_ARG italic_s end_ARG start_POSTSUPERSCRIPT italic_d end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_l ; italic_a end_POSTSUBSCRIPT ⋅ over→ start_ARG italic_s end_ARG start_POSTSUPERSCRIPT italic_p end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_l ; italic_a end_POSTSUBSCRIPT - divide start_ARG 1 end_ARG start_ARG 4 end_ARG italic_n start_POSTSUPERSCRIPT italic_p end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_l ; italic_a end_POSTSUBSCRIPT ] + ∑ start_POSTSUBSCRIPT italic_l end_POSTSUBSCRIPT 2 italic_J start_POSTSUBSCRIPT italic_K ; italic_z end_POSTSUBSCRIPT [ over→ start_ARG italic_s end_ARG start_POSTSUPERSCRIPT italic_d end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_l ; 2 end_POSTSUBSCRIPT ⋅ over→ start_ARG italic_s end_ARG start_POSTSUPERSCRIPT italic_p end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT - divide start_ARG 1 end_ARG start_ARG 4 end_ARG italic_n start_POSTSUPERSCRIPT italic_p end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT ] - ∑ start_POSTSUBSCRIPT italic_l end_POSTSUBSCRIPT 2 italic_J start_POSTSUBSCRIPT italic_H end_POSTSUBSCRIPT [ over→ start_ARG italic_s end_ARG start_POSTSUPERSCRIPT italic_d end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_l ; 1 end_POSTSUBSCRIPT ⋅ over→ start_ARG italic_s end_ARG start_POSTSUPERSCRIPT italic_d end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_l ; 2 end_POSTSUBSCRIPT + divide start_ARG 1 end_ARG start_ARG 4 end_ARG ] , (S3)

where JK;a,JK;zsubscript𝐽𝐾𝑎subscript𝐽𝐾𝑧J_{K;a},J_{K;z}italic_J start_POSTSUBSCRIPT italic_K ; italic_a end_POSTSUBSCRIPT , italic_J start_POSTSUBSCRIPT italic_K ; italic_z end_POSTSUBSCRIPT is the Kondo coupling derived in main text, Eq.(3). Based on the parameters listed in Ref.Wú et al. (2023), we assume tdp;1tdp;z>tdp;2similar-to-or-equalssubscript𝑡𝑑𝑝1subscript𝑡𝑑𝑝𝑧subscript𝑡𝑑𝑝2t_{dp;1}\simeq t_{dp;z}>t_{dp;2}italic_t start_POSTSUBSCRIPT italic_d italic_p ; 1 end_POSTSUBSCRIPT ≃ italic_t start_POSTSUBSCRIPT italic_d italic_p ; italic_z end_POSTSUBSCRIPT > italic_t start_POSTSUBSCRIPT italic_d italic_p ; 2 end_POSTSUBSCRIPT and Δp;1Δp;2zΔp;2similar-to-or-equalssubscriptΔ𝑝1superscriptsubscriptΔ𝑝2𝑧similar-to-or-equalssubscriptΔ𝑝2\Delta_{p;1}\simeq\Delta_{p;2}^{z}\simeq\Delta_{p;2}roman_Δ start_POSTSUBSCRIPT italic_p ; 1 end_POSTSUBSCRIPT ≃ roman_Δ start_POSTSUBSCRIPT italic_p ; 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_z end_POSTSUPERSCRIPT ≃ roman_Δ start_POSTSUBSCRIPT italic_p ; 2 end_POSTSUBSCRIPT, leading to the relation JK;14JK;z>JK;2subscript𝐽𝐾14subscript𝐽𝐾𝑧subscript𝐽𝐾2J_{K;1}\approx 4J_{K;z}>J_{K;2}italic_J start_POSTSUBSCRIPT italic_K ; 1 end_POSTSUBSCRIPT ≈ 4 italic_J start_POSTSUBSCRIPT italic_K ; italic_z end_POSTSUBSCRIPT > italic_J start_POSTSUBSCRIPT italic_K ; 2 end_POSTSUBSCRIPT. There are three possible states depending on the p𝑝pitalic_p orbital occupancy : (i)𝑖(i)( italic_i ) nl;1p=1superscriptsubscript𝑛𝑙1𝑝1n_{l;1}^{p}=1italic_n start_POSTSUBSCRIPT italic_l ; 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_p end_POSTSUPERSCRIPT = 1, (ii)𝑖𝑖(ii)( italic_i italic_i ) nl;2p=1superscriptsubscript𝑛𝑙2𝑝1n_{l;2}^{p}=1italic_n start_POSTSUBSCRIPT italic_l ; 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_p end_POSTSUPERSCRIPT = 1, and (iii)𝑖𝑖𝑖(iii)( italic_i italic_i italic_i ) nl;zp=1superscriptsubscript𝑛𝑙𝑧𝑝1n_{l;z}^{p}=1italic_n start_POSTSUBSCRIPT italic_l ; italic_z end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_p end_POSTSUPERSCRIPT = 1.

  • nl;1p=1superscriptsubscript𝑛𝑙1𝑝1n_{l;1}^{p}=1italic_n start_POSTSUBSCRIPT italic_l ; 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_p end_POSTSUPERSCRIPT = 1 or nl;2p=1superscriptsubscript𝑛𝑙2𝑝1n_{l;2}^{p}=1italic_n start_POSTSUBSCRIPT italic_l ; 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_p end_POSTSUPERSCRIPT = 1 case: The local Hamiltonian is reduced into three-site Hamiltonian written in terms of (pa;l,d1;l,d2;l)subscript𝑝𝑎𝑙subscript𝑑1𝑙subscript𝑑2𝑙(p_{a;l},d_{1;l},d_{2;l})( italic_p start_POSTSUBSCRIPT italic_a ; italic_l end_POSTSUBSCRIPT , italic_d start_POSTSUBSCRIPT 1 ; italic_l end_POSTSUBSCRIPT , italic_d start_POSTSUBSCRIPT 2 ; italic_l end_POSTSUBSCRIPT ),

    H𝐻\displaystyle Hitalic_H =\displaystyle== 2JK;a[sl;adsl;ap14]l2JH[sl;1dsl;2d+14].2subscript𝐽𝐾𝑎delimited-[]subscriptsuperscript𝑠𝑑𝑙𝑎subscriptsuperscript𝑠𝑝𝑙𝑎14subscript𝑙2subscript𝐽𝐻delimited-[]subscriptsuperscript𝑠𝑑𝑙1subscriptsuperscript𝑠𝑑𝑙214\displaystyle 2J_{K;a}[\vec{s}^{d}_{l;a}\cdot\vec{s}^{p}_{l;a}-\frac{1}{4}]-% \sum_{l}2J_{H}[\vec{s}^{d}_{l;1}\cdot\vec{s}^{d}_{l;2}+\frac{1}{4}].2 italic_J start_POSTSUBSCRIPT italic_K ; italic_a end_POSTSUBSCRIPT [ over→ start_ARG italic_s end_ARG start_POSTSUPERSCRIPT italic_d end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_l ; italic_a end_POSTSUBSCRIPT ⋅ over→ start_ARG italic_s end_ARG start_POSTSUPERSCRIPT italic_p end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_l ; italic_a end_POSTSUBSCRIPT - divide start_ARG 1 end_ARG start_ARG 4 end_ARG ] - ∑ start_POSTSUBSCRIPT italic_l end_POSTSUBSCRIPT 2 italic_J start_POSTSUBSCRIPT italic_H end_POSTSUBSCRIPT [ over→ start_ARG italic_s end_ARG start_POSTSUPERSCRIPT italic_d end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_l ; 1 end_POSTSUBSCRIPT ⋅ over→ start_ARG italic_s end_ARG start_POSTSUPERSCRIPT italic_d end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_l ; 2 end_POSTSUBSCRIPT + divide start_ARG 1 end_ARG start_ARG 4 end_ARG ] .

    The ground states of the above Hamiltonian is the spin 1/2121/21 / 2 state whose wave function is given by,

    |lsubscriptket𝑙\displaystyle\ket{{\uparrow}}_{l}| start_ARG ↑ end_ARG ⟩ start_POSTSUBSCRIPT italic_l end_POSTSUBSCRIPT =\displaystyle== 12[(r+α)2+(r+α+1)][(1+r+α)|,,(r+α)|,,|,,],12delimited-[]superscript𝑟𝛼2𝑟𝛼1delimited-[]1𝑟𝛼ket𝑟𝛼ketket\displaystyle\frac{1}{\sqrt{2[(r+\alpha)^{2}+(r+\alpha+1)]}}[(1+r+\alpha)|% \downarrow,\uparrow,\uparrow\rangle-(r+\alpha)|\uparrow,\downarrow,\uparrow% \rangle-|\uparrow,\uparrow,\downarrow\rangle],divide start_ARG 1 end_ARG start_ARG square-root start_ARG 2 [ ( italic_r + italic_α ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + ( italic_r + italic_α + 1 ) ] end_ARG end_ARG [ ( 1 + italic_r + italic_α ) | ↓ , ↑ , ↑ ⟩ - ( italic_r + italic_α ) | ↑ , ↓ , ↑ ⟩ - | ↑ , ↑ , ↓ ⟩ ] , (S4)
    |lsubscriptket𝑙\displaystyle\ket{{\downarrow}}_{l}| start_ARG ↓ end_ARG ⟩ start_POSTSUBSCRIPT italic_l end_POSTSUBSCRIPT =\displaystyle== 12[(r+α)2+(r+α+1)][(1+r+α)|,,(r+α)|,,|,,].12delimited-[]superscript𝑟𝛼2𝑟𝛼1delimited-[]1𝑟𝛼ket𝑟𝛼ketket\displaystyle\frac{-1}{\sqrt{2[(r+\alpha)^{2}+(r+\alpha+1)]}}[(1+r+\alpha)|% \uparrow,\downarrow,\downarrow\rangle-(r+\alpha)|\downarrow,\uparrow,% \downarrow\rangle-|\downarrow,\uparrow,\downarrow\rangle].divide start_ARG - 1 end_ARG start_ARG square-root start_ARG 2 [ ( italic_r + italic_α ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + ( italic_r + italic_α + 1 ) ] end_ARG end_ARG [ ( 1 + italic_r + italic_α ) | ↑ , ↓ , ↓ ⟩ - ( italic_r + italic_α ) | ↓ , ↑ , ↓ ⟩ - | ↓ , ↑ , ↓ ⟩ ] . (S5)

    where |s0,s1,s2ketsubscript𝑠0subscript𝑠1subscript𝑠2\ket{s_{0},s_{1},s_{2}}| start_ARG italic_s start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT , italic_s start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , italic_s start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_ARG ⟩ implies the spin of the pa,d1,d2subscript𝑝𝑎subscript𝑑1subscript𝑑2p_{a},d_{1},d_{2}italic_p start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT , italic_d start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , italic_d start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT orbital respectively. Here, we use r=JK;a/JH𝑟subscript𝐽𝐾𝑎subscript𝐽𝐻r=J_{K;a}/J_{H}italic_r = italic_J start_POSTSUBSCRIPT italic_K ; italic_a end_POSTSUBSCRIPT / italic_J start_POSTSUBSCRIPT italic_H end_POSTSUBSCRIPT and α(r)=r2+r+1𝛼𝑟superscript𝑟2𝑟1\alpha(r)=\sqrt{r^{2}+r+1}italic_α ( italic_r ) = square-root start_ARG italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + italic_r + 1 end_ARG. These states are called as Zhang-Rice spin-half state, and the schematic illustrations are provided in Fig. S1. The wave function is obviously entangled by the three (p1,d1,d2subscript𝑝1subscript𝑑1subscript𝑑2p_{1},d_{1},d_{2}italic_p start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , italic_d start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , italic_d start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT) orbitals in the most range of r𝑟ritalic_r, except r=+𝑟r=+\inftyitalic_r = + ∞. In this specific point r=+𝑟r=+\inftyitalic_r = + ∞, the state becomes the tensor product of the Zhang-Rice singlet only made by (p1,d1subscript𝑝1subscript𝑑1p_{1},d_{1}italic_p start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , italic_d start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT) and decoupled d2subscript𝑑2d_{2}italic_d start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT orbitals. The subsystem entanglement entropy is obtained in Fig. S2(b), manifesting its values is finite due to the entanglement of (p1,d1subscript𝑝1subscript𝑑1p_{1},d_{1}italic_p start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , italic_d start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT) and d2subscript𝑑2d_{2}italic_d start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT. The associated energy is obtained as,

    EG=[1+r+r2+r+1]JH,subscript𝐸𝐺delimited-[]1𝑟superscript𝑟2𝑟1subscript𝐽𝐻\displaystyle E_{G}=-[1+r+\sqrt{r^{2}+r+1}]J_{H},italic_E start_POSTSUBSCRIPT italic_G end_POSTSUBSCRIPT = - [ 1 + italic_r + square-root start_ARG italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + italic_r + 1 end_ARG ] italic_J start_POSTSUBSCRIPT italic_H end_POSTSUBSCRIPT , (S6)

    introducing the dimensionless parameter, r=JK;a/JH𝑟subscript𝐽𝐾𝑎subscript𝐽𝐻r=J_{K;a}/J_{H}italic_r = italic_J start_POSTSUBSCRIPT italic_K ; italic_a end_POSTSUBSCRIPT / italic_J start_POSTSUBSCRIPT italic_H end_POSTSUBSCRIPT.

  • nzp=1superscriptsubscript𝑛𝑧𝑝1n_{z}^{p}=1italic_n start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_p end_POSTSUPERSCRIPT = 1 case: The local Hamiltonian is reduced into five-site Hamiltonian of (pz,d1;t,d1;b,d2;t,d2;b)subscript𝑝𝑧subscript𝑑1𝑡subscript𝑑1𝑏subscript𝑑2𝑡subscript𝑑2𝑏(p_{z},d_{1;t},d_{1;b},d_{2;t},d_{2;b})( italic_p start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT , italic_d start_POSTSUBSCRIPT 1 ; italic_t end_POSTSUBSCRIPT , italic_d start_POSTSUBSCRIPT 1 ; italic_b end_POSTSUBSCRIPT , italic_d start_POSTSUBSCRIPT 2 ; italic_t end_POSTSUBSCRIPT , italic_d start_POSTSUBSCRIPT 2 ; italic_b end_POSTSUBSCRIPT ),

    H𝐻\displaystyle Hitalic_H =\displaystyle== l2JK;z[sl;2dszp14]l2JH[sl;1dsl;2d+14].subscript𝑙2subscript𝐽𝐾𝑧delimited-[]subscriptsuperscript𝑠𝑑𝑙2subscriptsuperscript𝑠𝑝𝑧14subscript𝑙2subscript𝐽𝐻delimited-[]subscriptsuperscript𝑠𝑑𝑙1subscriptsuperscript𝑠𝑑𝑙214\displaystyle\sum_{l}2J_{K;z}[\vec{s}^{d}_{l;2}\cdot\vec{s}^{p}_{z}-\frac{1}{4% }]-\sum_{l}2J_{H}[\vec{s}^{d}_{l;1}\cdot\vec{s}^{d}_{l;2}+\frac{1}{4}].∑ start_POSTSUBSCRIPT italic_l end_POSTSUBSCRIPT 2 italic_J start_POSTSUBSCRIPT italic_K ; italic_z end_POSTSUBSCRIPT [ over→ start_ARG italic_s end_ARG start_POSTSUPERSCRIPT italic_d end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_l ; 2 end_POSTSUBSCRIPT ⋅ over→ start_ARG italic_s end_ARG start_POSTSUPERSCRIPT italic_p end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT - divide start_ARG 1 end_ARG start_ARG 4 end_ARG ] - ∑ start_POSTSUBSCRIPT italic_l end_POSTSUBSCRIPT 2 italic_J start_POSTSUBSCRIPT italic_H end_POSTSUBSCRIPT [ over→ start_ARG italic_s end_ARG start_POSTSUPERSCRIPT italic_d end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_l ; 1 end_POSTSUBSCRIPT ⋅ over→ start_ARG italic_s end_ARG start_POSTSUPERSCRIPT italic_d end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_l ; 2 end_POSTSUBSCRIPT + divide start_ARG 1 end_ARG start_ARG 4 end_ARG ] .

    The wave function |ψket𝜓\ket{\psi}| start_ARG italic_ψ end_ARG ⟩ of the ground states of the above Hamiltonian have four-fold degeneracy and are very complicated. However, the important thing to note is the total spin of the state is Stot=32subscript𝑆tot32S_{\mathrm{tot}}=\frac{3}{2}italic_S start_POSTSUBSCRIPT roman_tot end_POSTSUBSCRIPT = divide start_ARG 3 end_ARG start_ARG 2 end_ARG, i.e. Stot2|ψ=154|ψsuperscriptsubscript𝑆𝑡𝑜𝑡2ket𝜓154ket𝜓S_{tot}^{2}\ket{\psi}=\frac{15}{4}\ket{\psi}italic_S start_POSTSUBSCRIPT italic_t italic_o italic_t end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT | start_ARG italic_ψ end_ARG ⟩ = divide start_ARG 15 end_ARG start_ARG 4 end_ARG | start_ARG italic_ψ end_ARG ⟩, and Stotz|ψ=±12,±32|ψsuperscriptsubscript𝑆𝑡𝑜𝑡𝑧ket𝜓plus-or-minus12plus-or-minus32ket𝜓S_{tot}^{z}\ket{\psi}=\pm\frac{1}{2},\pm\frac{3}{2}\ket{\psi}italic_S start_POSTSUBSCRIPT italic_t italic_o italic_t end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_z end_POSTSUPERSCRIPT | start_ARG italic_ψ end_ARG ⟩ = ± divide start_ARG 1 end_ARG start_ARG 2 end_ARG , ± divide start_ARG 3 end_ARG start_ARG 2 end_ARG | start_ARG italic_ψ end_ARG ⟩, for any values of positive JK;zsubscript𝐽𝐾𝑧J_{K;z}italic_J start_POSTSUBSCRIPT italic_K ; italic_z end_POSTSUBSCRIPT and JHsubscript𝐽𝐻J_{H}italic_J start_POSTSUBSCRIPT italic_H end_POSTSUBSCRIPT. The associated ground state energy is obtained as

    EGz=12[23r4+8r+9r2]JH,superscriptsubscript𝐸𝐺𝑧12delimited-[]23𝑟48𝑟9superscript𝑟2subscript𝐽𝐻\displaystyle E_{G}^{z}=\frac{1}{2}[-2-3r-\sqrt{4+8r+9r^{2}}]J_{H},italic_E start_POSTSUBSCRIPT italic_G end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_z end_POSTSUPERSCRIPT = divide start_ARG 1 end_ARG start_ARG 2 end_ARG [ - 2 - 3 italic_r - square-root start_ARG 4 + 8 italic_r + 9 italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG ] italic_J start_POSTSUBSCRIPT italic_H end_POSTSUBSCRIPT , (S7)

    with a dimensionless parameter, r=JK;z/JH𝑟subscript𝐽𝐾𝑧subscript𝐽𝐻r=J_{K;z}/J_{H}italic_r = italic_J start_POSTSUBSCRIPT italic_K ; italic_z end_POSTSUBSCRIPT / italic_J start_POSTSUBSCRIPT italic_H end_POSTSUBSCRIPT,

The above-hybridized states are illustrated in Fig. S1. We now compare those three energies, Eqs.(S6-S7), using the fact JK;14JK;z>JK;2subscript𝐽𝐾14subscript𝐽𝐾𝑧subscript𝐽𝐾2J_{K;1}\approx 4J_{K;z}>J_{K;2}italic_J start_POSTSUBSCRIPT italic_K ; 1 end_POSTSUBSCRIPT ≈ 4 italic_J start_POSTSUBSCRIPT italic_K ; italic_z end_POSTSUBSCRIPT > italic_J start_POSTSUBSCRIPT italic_K ; 2 end_POSTSUBSCRIPT. Obviously, EG1(JK;1)<EG2(JK;2)superscriptsubscript𝐸𝐺1subscript𝐽𝐾1superscriptsubscript𝐸𝐺2subscript𝐽𝐾2E_{G}^{1}(J_{K;1})<E_{G}^{2}(J_{K;2})italic_E start_POSTSUBSCRIPT italic_G end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 1 end_POSTSUPERSCRIPT ( italic_J start_POSTSUBSCRIPT italic_K ; 1 end_POSTSUBSCRIPT ) < italic_E start_POSTSUBSCRIPT italic_G end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( italic_J start_POSTSUBSCRIPT italic_K ; 2 end_POSTSUBSCRIPT ) is established. Next, we compare the energy between EG1superscriptsubscript𝐸𝐺1E_{G}^{1}italic_E start_POSTSUBSCRIPT italic_G end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 1 end_POSTSUPERSCRIPT and EGzsuperscriptsubscript𝐸𝐺𝑧E_{G}^{z}italic_E start_POSTSUBSCRIPT italic_G end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_z end_POSTSUPERSCRIPT in Fig. S2 showing that EG1(JK;1)<EGz(JK;z4JK;1)superscriptsubscript𝐸𝐺1subscript𝐽𝐾1superscriptsubscript𝐸𝐺𝑧similar-to-or-equalssubscript𝐽𝐾𝑧4subscript𝐽𝐾1E_{G}^{1}(J_{K;1})<E_{G}^{z}(J_{K;z}\simeq 4J_{K;1})italic_E start_POSTSUBSCRIPT italic_G end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 1 end_POSTSUPERSCRIPT ( italic_J start_POSTSUBSCRIPT italic_K ; 1 end_POSTSUBSCRIPT ) < italic_E start_POSTSUBSCRIPT italic_G end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_z end_POSTSUPERSCRIPT ( italic_J start_POSTSUBSCRIPT italic_K ; italic_z end_POSTSUBSCRIPT ≃ 4 italic_J start_POSTSUBSCRIPT italic_K ; 1 end_POSTSUBSCRIPT ). Hence, we can conclude that the Zhang-Rice spin-1/2 state, especially with p1subscript𝑝1p_{1}italic_p start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT hole is the most stable state.

Refer to caption
Figure S1: Possible Zhang-Rice states of d8Lsuperscript𝑑8𝐿d^{8}Litalic_d start_POSTSUPERSCRIPT 8 end_POSTSUPERSCRIPT italic_L electron configurations of the bilayer nickelates. (a) The spin 1/2121/21 / 2 and (b) spin 3/2323/23 / 2 states are hybridized by intra-layer (p1,p2)subscript𝑝1subscript𝑝2(p_{1},p_{2})( italic_p start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , italic_p start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ) orbital or inter-layer pzsubscript𝑝𝑧p_{z}italic_p start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT orbital of oxygen atoms. Note that here all of the states are drawn in the hole description picture.
Refer to caption
Figure S2: (a) The JK/JH=rsubscript𝐽𝐾subscript𝐽𝐻𝑟J_{K}/J_{H}=ritalic_J start_POSTSUBSCRIPT italic_K end_POSTSUBSCRIPT / italic_J start_POSTSUBSCRIPT italic_H end_POSTSUBSCRIPT = italic_r dependence of EG1(r)superscriptsubscript𝐸𝐺1𝑟E_{G}^{1}(r)italic_E start_POSTSUBSCRIPT italic_G end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 1 end_POSTSUPERSCRIPT ( italic_r ) and EGz(r)superscriptsubscript𝐸𝐺𝑧𝑟E_{G}^{z}(r)italic_E start_POSTSUBSCRIPT italic_G end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_z end_POSTSUPERSCRIPT ( italic_r ). (b)Comparison of the ground states energy in (JK;1,JK;z)subscript𝐽𝐾1subscript𝐽𝐾𝑧(J_{K;1},J_{K;z})( italic_J start_POSTSUBSCRIPT italic_K ; 1 end_POSTSUBSCRIPT , italic_J start_POSTSUBSCRIPT italic_K ; italic_z end_POSTSUBSCRIPT ) plane. The blue (green) region indicates the phase space where p1subscript𝑝1p_{1}italic_p start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT (pzsubscript𝑝𝑧p_{z}italic_p start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT) orbitals energetically favored, so an additonal holes enter to p1subscript𝑝1p_{1}italic_p start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT (pzsubscript𝑝𝑧p_{z}italic_p start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT) orbitals. Using the estimation JK;14JK;zsimilar-tosubscript𝐽𝐾14subscript𝐽𝐾𝑧J_{K;1}\sim 4J_{K;z}italic_J start_POSTSUBSCRIPT italic_K ; 1 end_POSTSUBSCRIPT ∼ 4 italic_J start_POSTSUBSCRIPT italic_K ; italic_z end_POSTSUBSCRIPT (red dashed line), we conclude that the Zhang-Rice spin-1/2 state with p1subscript𝑝1p_{1}italic_p start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT hole is more stable. (c) The bipartite entanglement entropy of |ket\ket{\uparrow}| start_ARG ↑ end_ARG ⟩, Eqs.(S4-S5), as a function of r𝑟ritalic_r. Here, the entanglement entropy is evaluated by S=Trd1(ρ^logρ^)𝑆subscriptTrsubscript𝑑1^𝜌^𝜌S=\mathrm{Tr}_{d_{1}}({\widehat{\rho}}\log{\widehat{\rho}})italic_S = roman_Tr start_POSTSUBSCRIPT italic_d start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_POSTSUBSCRIPT ( over^ start_ARG italic_ρ end_ARG roman_log over^ start_ARG italic_ρ end_ARG ), with the reduced density matrix, ρ^=Trd1,p[||]^𝜌subscriptTrsubscript𝑑1𝑝delimited-[]ketbra{\widehat{\rho}}=\mathrm{Tr}_{d_{1},p}[\ket{\uparrow}\bra{\uparrow}]over^ start_ARG italic_ρ end_ARG = roman_Tr start_POSTSUBSCRIPT italic_d start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , italic_p end_POSTSUBSCRIPT [ | start_ARG ↑ end_ARG ⟩ ⟨ start_ARG ↑ end_ARG | ]. As r+𝑟r\rightarrow+\inftyitalic_r → + ∞ limit, Zhang-Rice spin half state is just a tensor product of the singlet formed by the (d1,p1subscript𝑑1subscript𝑝1d_{1},p_{1}italic_d start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , italic_p start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT) and spin 1/2 of d2subscript𝑑2d_{2}italic_d start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT orbital, and the entanglement between two subsystems goes to zero. However, in other value of r𝑟ritalic_r, the d2subscript𝑑2d_{2}italic_d start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT orbital does not be decoupled from d1,p1subscript𝑑1subscript𝑝1d_{1},p_{1}italic_d start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , italic_p start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT orbitals, forming the Zhang-Rice spin half state together.

III III. Derivation of type-II t-J model from the charge transfer Hamiltonian

In this section, we derive the type-II t-J model starting from the charge transfer Hamiltonian, Eq. (1). In particular, we show all the spin-exchange couplings of the type-II t-J model, in terms of the original Hamiltonian. The spin-exchange interaction of two nickel atoms can be derived by the fourth-order perturbation in the strong coupling limit. We divide the charge-transfer Hamiltonian into the kinetic part, V=Hdp𝑉subscript𝐻𝑑𝑝V=H_{dp}italic_V = italic_H start_POSTSUBSCRIPT italic_d italic_p end_POSTSUBSCRIPT, and remaining part, H0=HVsubscript𝐻0𝐻𝑉H_{0}=H-Vitalic_H start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = italic_H - italic_V.

III.1 A. Derivation of Jsubscript𝐽perpendicular-toJ_{\perp}italic_J start_POSTSUBSCRIPT ⟂ end_POSTSUBSCRIPT

First, we show how the Jsuperscript𝐽perpendicular-toJ^{\perp}italic_J start_POSTSUPERSCRIPT ⟂ end_POSTSUPERSCRIPT terms are derived, exemplified by Jsdsubscriptsuperscript𝐽perpendicular-to𝑠𝑑J^{\perp}_{sd}italic_J start_POSTSUPERSCRIPT ⟂ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_s italic_d end_POSTSUBSCRIPT, and the generalization into Jss,Jddsubscriptsuperscript𝐽perpendicular-to𝑠𝑠subscriptsuperscript𝐽perpendicular-to𝑑𝑑J^{\perp}_{ss},J^{\perp}_{dd}italic_J start_POSTSUPERSCRIPT ⟂ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_s italic_s end_POSTSUBSCRIPT , italic_J start_POSTSUPERSCRIPT ⟂ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_d italic_d end_POSTSUBSCRIPT is straightforward. Due to the SU(2) symmetry of the Hamiltonian, it is enough to evaluate the coefficient of one of the components of 12(stSb+)=12|st=12,sb=0st=12,sb=1|+12subscriptsuperscript𝑠𝑡subscriptsuperscript𝑆𝑏12ketformulae-sequencesubscript𝑠𝑡12subscript𝑠𝑏0braformulae-sequencesubscript𝑠𝑡12subscript𝑠𝑏1\frac{1}{2}(s^{-}_{t}S^{+}_{b})=\frac{1}{\sqrt{2}}\ket{s_{t}=-\frac{1}{2},s_{b% }=0}\bra{s_{t}=\frac{1}{2},s_{b}=-1}+\cdotsdivide start_ARG 1 end_ARG start_ARG 2 end_ARG ( italic_s start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT italic_S start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT ) = divide start_ARG 1 end_ARG start_ARG square-root start_ARG 2 end_ARG end_ARG | start_ARG italic_s start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT = - divide start_ARG 1 end_ARG start_ARG 2 end_ARG , italic_s start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT = 0 end_ARG ⟩ ⟨ start_ARG italic_s start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT = divide start_ARG 1 end_ARG start_ARG 2 end_ARG , italic_s start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT = - 1 end_ARG | + ⋯ piece, not all st;iSb;isubscript𝑠𝑡𝑖subscript𝑆𝑏𝑖\vec{s}_{t;i}\cdot\vec{S}_{b;i}over→ start_ARG italic_s end_ARG start_POSTSUBSCRIPT italic_t ; italic_i end_POSTSUBSCRIPT ⋅ over→ start_ARG italic_S end_ARG start_POSTSUBSCRIPT italic_b ; italic_i end_POSTSUBSCRIPT. There are six different processes (A)-(F) contributing the stSb+subscriptsuperscript𝑠𝑡subscriptsuperscript𝑆𝑏s^{-}_{t}S^{+}_{b}italic_s start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT italic_S start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT, as illustrated in Fig.S3 (a),

st=12,sb=0|VRVRVRV|st=12,sb=1braformulae-sequencesubscript𝑠𝑡12subscript𝑠𝑏0𝑉𝑅𝑉𝑅𝑉𝑅𝑉ketformulae-sequencesubscript𝑠𝑡12subscript𝑠𝑏1\displaystyle\bra{s_{t}=-\frac{1}{2},s_{b}=0}VRVRVRV\ket{s_{t}=\frac{1}{2},s_{% b}=-1}⟨ start_ARG italic_s start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT = - divide start_ARG 1 end_ARG start_ARG 2 end_ARG , italic_s start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT = 0 end_ARG | italic_V italic_R italic_V italic_R italic_V italic_R italic_V | start_ARG italic_s start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT = divide start_ARG 1 end_ARG start_ARG 2 end_ARG , italic_s start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT = - 1 end_ARG ⟩ =\displaystyle== (A)+(B)+(C)+(D)+(E)+(F),𝐴𝐵𝐶𝐷𝐸𝐹\displaystyle(A)+(B)+(C)+(D)+(E)+(F),( italic_A ) + ( italic_B ) + ( italic_C ) + ( italic_D ) + ( italic_E ) + ( italic_F ) , (S8)

where R=Q/(E0H0)𝑅𝑄subscript𝐸0subscript𝐻0R=Q/(E_{0}-H_{0})italic_R = italic_Q / ( italic_E start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT - italic_H start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) is the projection onto the subspace of all exited states. Here, Q=1P𝑄1𝑃Q=1-Pitalic_Q = 1 - italic_P is a complementary operator of the projection operator of the ground state of H0subscript𝐻0H_{0}italic_H start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT where the doubly occupied states are forbidden.

With some tedious calculations, we find each component is given by

(A)𝐴\displaystyle(A)( italic_A ) =\displaystyle== tdp,z(Δp;2z+J¯H)12|t0|bVRVRVpz;[(1+r+α)𝒩|sd1=12,sc=12t1𝒩|sd1=12,sc=12t]|1bsubscript𝑡𝑑𝑝𝑧superscriptsubscriptΔ𝑝2𝑧subscript¯𝐽𝐻subscriptbra12𝑡subscriptbra0𝑏𝑉𝑅𝑉𝑅𝑉superscriptsubscript𝑝𝑧delimited-[]1𝑟𝛼𝒩subscriptketformulae-sequencesubscript𝑠subscript𝑑112subscript𝑠𝑐12𝑡1𝒩subscriptketformulae-sequencesubscript𝑠subscript𝑑112subscript𝑠𝑐12𝑡subscriptket1𝑏\displaystyle\frac{t_{dp,z}}{(\Delta_{p;2}^{z}+\overline{J}_{H})}\bra{-\frac{1% }{2}}_{t}\bra{0}_{b}VRVRVp_{z;\uparrow}^{\dagger}[\frac{(1+r+\alpha)}{\mathcal% {N}}\ket{s_{d_{1}}=\frac{1}{2},s_{c}=-\frac{1}{2}}_{t}-\frac{1}{\mathcal{N}}% \ket{s_{d_{1}}=-\frac{1}{2},s_{c}=\frac{1}{2}}_{t}]\ket{-1}_{b}divide start_ARG italic_t start_POSTSUBSCRIPT italic_d italic_p , italic_z end_POSTSUBSCRIPT end_ARG start_ARG ( roman_Δ start_POSTSUBSCRIPT italic_p ; 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_z end_POSTSUPERSCRIPT + over¯ start_ARG italic_J end_ARG start_POSTSUBSCRIPT italic_H end_POSTSUBSCRIPT ) end_ARG ⟨ start_ARG - divide start_ARG 1 end_ARG start_ARG 2 end_ARG end_ARG | start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT ⟨ start_ARG 0 end_ARG | start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT italic_V italic_R italic_V italic_R italic_V italic_p start_POSTSUBSCRIPT italic_z ; ↑ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT [ divide start_ARG ( 1 + italic_r + italic_α ) end_ARG start_ARG caligraphic_N end_ARG | start_ARG italic_s start_POSTSUBSCRIPT italic_d start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_POSTSUBSCRIPT = divide start_ARG 1 end_ARG start_ARG 2 end_ARG , italic_s start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT = - divide start_ARG 1 end_ARG start_ARG 2 end_ARG end_ARG ⟩ start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT - divide start_ARG 1 end_ARG start_ARG caligraphic_N end_ARG | start_ARG italic_s start_POSTSUBSCRIPT italic_d start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_POSTSUBSCRIPT = - divide start_ARG 1 end_ARG start_ARG 2 end_ARG , italic_s start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT = divide start_ARG 1 end_ARG start_ARG 2 end_ARG end_ARG ⟩ start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT ] | start_ARG - 1 end_ARG ⟩ start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT
=\displaystyle== (tdp,z)2(Δp;2z+J¯H)(U+2JH)12|t0|bVRVd2;b;[(1+r+α)𝒩|sd1=12,sc=12t1𝒩|sd1=12,sc=12]t|1bsuperscriptsubscript𝑡𝑑𝑝𝑧2superscriptsubscriptΔ𝑝2𝑧subscript¯𝐽𝐻𝑈2subscript𝐽𝐻subscriptbra12𝑡subscriptbra0𝑏𝑉𝑅𝑉superscriptsubscript𝑑2𝑏subscriptdelimited-[]1𝑟𝛼𝒩subscriptketformulae-sequencesubscript𝑠subscript𝑑112subscript𝑠𝑐12𝑡1𝒩ketformulae-sequencesubscript𝑠subscript𝑑112subscript𝑠𝑐12𝑡subscriptket1𝑏\displaystyle\frac{(t_{dp,z})^{2}}{(\Delta_{p;2}^{z}+\overline{J}_{H})(U+2J_{H% })}\bra{-\frac{1}{2}}_{t}\bra{0}_{b}VRVd_{2;b;\uparrow}^{\dagger}[\frac{(1+r+% \alpha)}{\mathcal{N}}\ket{s_{d_{1}}=\frac{1}{2},s_{c}=-\frac{1}{2}}_{t}-\frac{% 1}{\mathcal{N}}\ket{s_{d_{1}}=-\frac{1}{2},s_{c}=\frac{1}{2}}]_{t}\ket{-1}_{b}divide start_ARG ( italic_t start_POSTSUBSCRIPT italic_d italic_p , italic_z end_POSTSUBSCRIPT ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG ( roman_Δ start_POSTSUBSCRIPT italic_p ; 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_z end_POSTSUPERSCRIPT + over¯ start_ARG italic_J end_ARG start_POSTSUBSCRIPT italic_H end_POSTSUBSCRIPT ) ( italic_U + 2 italic_J start_POSTSUBSCRIPT italic_H end_POSTSUBSCRIPT ) end_ARG ⟨ start_ARG - divide start_ARG 1 end_ARG start_ARG 2 end_ARG end_ARG | start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT ⟨ start_ARG 0 end_ARG | start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT italic_V italic_R italic_V italic_d start_POSTSUBSCRIPT 2 ; italic_b ; ↑ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT [ divide start_ARG ( 1 + italic_r + italic_α ) end_ARG start_ARG caligraphic_N end_ARG | start_ARG italic_s start_POSTSUBSCRIPT italic_d start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_POSTSUBSCRIPT = divide start_ARG 1 end_ARG start_ARG 2 end_ARG , italic_s start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT = - divide start_ARG 1 end_ARG start_ARG 2 end_ARG end_ARG ⟩ start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT - divide start_ARG 1 end_ARG start_ARG caligraphic_N end_ARG | start_ARG italic_s start_POSTSUBSCRIPT italic_d start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_POSTSUBSCRIPT = - divide start_ARG 1 end_ARG start_ARG 2 end_ARG , italic_s start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT = divide start_ARG 1 end_ARG start_ARG 2 end_ARG end_ARG ⟩ ] start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT | start_ARG - 1 end_ARG ⟩ start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT
=\displaystyle== (tdp,z)3(Δp;2z+J¯H)2(U+2JH)12|t0|bVpz;[(1+r+α)|sd1=12,sc=12t|sd1=12,sc=12t]|0b2superscriptsubscript𝑡𝑑𝑝𝑧3superscriptsuperscriptsubscriptΔ𝑝2𝑧subscript¯𝐽𝐻2𝑈2subscript𝐽𝐻subscriptbra12𝑡subscriptbra0𝑏𝑉superscriptsubscript𝑝𝑧delimited-[]1𝑟𝛼subscriptketformulae-sequencesubscript𝑠subscript𝑑112subscript𝑠𝑐12𝑡subscriptketformulae-sequencesubscript𝑠subscript𝑑112subscript𝑠𝑐12𝑡subscriptket0𝑏2\displaystyle\frac{-(t_{dp,z})^{3}}{(\Delta_{p;2}^{z}+\overline{J}_{H})^{2}(U+% 2J_{H})}\bra{-\frac{1}{2}}_{t}\bra{0}_{b}Vp_{z;\downarrow}^{\dagger}[(1+r+% \alpha)\ket{s_{d_{1}}=\frac{1}{2},s_{c}=-\frac{1}{2}}_{t}-\ket{s_{d_{1}}=-% \frac{1}{2},s_{c}=\frac{1}{2}}_{t}]\frac{\ket{0}_{b}}{\sqrt{2}}divide start_ARG - ( italic_t start_POSTSUBSCRIPT italic_d italic_p , italic_z end_POSTSUBSCRIPT ) start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT end_ARG start_ARG ( roman_Δ start_POSTSUBSCRIPT italic_p ; 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_z end_POSTSUPERSCRIPT + over¯ start_ARG italic_J end_ARG start_POSTSUBSCRIPT italic_H end_POSTSUBSCRIPT ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( italic_U + 2 italic_J start_POSTSUBSCRIPT italic_H end_POSTSUBSCRIPT ) end_ARG ⟨ start_ARG - divide start_ARG 1 end_ARG start_ARG 2 end_ARG end_ARG | start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT ⟨ start_ARG 0 end_ARG | start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT italic_V italic_p start_POSTSUBSCRIPT italic_z ; ↓ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT [ ( 1 + italic_r + italic_α ) | start_ARG italic_s start_POSTSUBSCRIPT italic_d start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_POSTSUBSCRIPT = divide start_ARG 1 end_ARG start_ARG 2 end_ARG , italic_s start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT = - divide start_ARG 1 end_ARG start_ARG 2 end_ARG end_ARG ⟩ start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT - | start_ARG italic_s start_POSTSUBSCRIPT italic_d start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_POSTSUBSCRIPT = - divide start_ARG 1 end_ARG start_ARG 2 end_ARG , italic_s start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT = divide start_ARG 1 end_ARG start_ARG 2 end_ARG end_ARG ⟩ start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT ] divide start_ARG | start_ARG 0 end_ARG ⟩ start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT end_ARG start_ARG square-root start_ARG 2 end_ARG end_ARG
=\displaystyle== (tdp,z)4(Δp;2z+J¯H)2(U+2JH)2(1+r+α)𝒩2=(D),superscriptsubscript𝑡𝑑𝑝𝑧4superscriptsuperscriptsubscriptΔ𝑝2𝑧subscript¯𝐽𝐻2𝑈2subscript𝐽𝐻21𝑟𝛼superscript𝒩2𝐷\displaystyle\frac{(t_{dp,z})^{4}}{(\Delta_{p;2}^{z}+\overline{J}_{H})^{2}(U+2% J_{H})}\frac{\sqrt{2}(1+r+\alpha)}{\mathcal{N}^{2}}=(D),divide start_ARG ( italic_t start_POSTSUBSCRIPT italic_d italic_p , italic_z end_POSTSUBSCRIPT ) start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT end_ARG start_ARG ( roman_Δ start_POSTSUBSCRIPT italic_p ; 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_z end_POSTSUPERSCRIPT + over¯ start_ARG italic_J end_ARG start_POSTSUBSCRIPT italic_H end_POSTSUBSCRIPT ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( italic_U + 2 italic_J start_POSTSUBSCRIPT italic_H end_POSTSUBSCRIPT ) end_ARG divide start_ARG square-root start_ARG 2 end_ARG ( 1 + italic_r + italic_α ) end_ARG start_ARG caligraphic_N start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG = ( italic_D ) ,

and

(B)𝐵\displaystyle(B)( italic_B ) =\displaystyle== (tdp,z)22(Δp;2z+J¯H)212|t0|bVRVpz;pz;[(1+r+α)𝒩|sd1=12,sc=12t1𝒩|sd1=12,sc=12t]|sd1=12bsuperscriptsubscript𝑡𝑑𝑝𝑧22superscriptsuperscriptsubscriptΔ𝑝2𝑧subscript¯𝐽𝐻2subscriptbra12𝑡subscriptbra0𝑏𝑉𝑅𝑉superscriptsubscript𝑝𝑧superscriptsubscript𝑝𝑧delimited-[]1𝑟𝛼𝒩subscriptketformulae-sequencesubscript𝑠subscript𝑑112subscript𝑠𝑐12𝑡1𝒩subscriptketformulae-sequencesubscript𝑠subscript𝑑112subscript𝑠𝑐12𝑡subscriptketsubscript𝑠subscript𝑑112𝑏\displaystyle\frac{(t_{dp,z})^{2}}{2(\Delta_{p;2}^{z}+\overline{J}_{H})^{2}}% \bra{-\frac{1}{2}}_{t}\bra{0}_{b}VRVp_{z;\downarrow}^{\dagger}p_{z;\uparrow}^{% \dagger}[\frac{(1+r+\alpha)}{\mathcal{N}}\ket{s_{d_{1}}=\frac{1}{2},s_{c}=-% \frac{1}{2}}_{t}-\frac{1}{\mathcal{N}}\ket{s_{d_{1}}=-\frac{1}{2},s_{c}=\frac{% 1}{2}}_{t}]\ket{s_{d_{1}}=-\frac{1}{2}}_{b}divide start_ARG ( italic_t start_POSTSUBSCRIPT italic_d italic_p , italic_z end_POSTSUBSCRIPT ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG 2 ( roman_Δ start_POSTSUBSCRIPT italic_p ; 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_z end_POSTSUPERSCRIPT + over¯ start_ARG italic_J end_ARG start_POSTSUBSCRIPT italic_H end_POSTSUBSCRIPT ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG ⟨ start_ARG - divide start_ARG 1 end_ARG start_ARG 2 end_ARG end_ARG | start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT ⟨ start_ARG 0 end_ARG | start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT italic_V italic_R italic_V italic_p start_POSTSUBSCRIPT italic_z ; ↓ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT italic_p start_POSTSUBSCRIPT italic_z ; ↑ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT [ divide start_ARG ( 1 + italic_r + italic_α ) end_ARG start_ARG caligraphic_N end_ARG | start_ARG italic_s start_POSTSUBSCRIPT italic_d start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_POSTSUBSCRIPT = divide start_ARG 1 end_ARG start_ARG 2 end_ARG , italic_s start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT = - divide start_ARG 1 end_ARG start_ARG 2 end_ARG end_ARG ⟩ start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT - divide start_ARG 1 end_ARG start_ARG caligraphic_N end_ARG | start_ARG italic_s start_POSTSUBSCRIPT italic_d start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_POSTSUBSCRIPT = - divide start_ARG 1 end_ARG start_ARG 2 end_ARG , italic_s start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT = divide start_ARG 1 end_ARG start_ARG 2 end_ARG end_ARG ⟩ start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT ] | start_ARG italic_s start_POSTSUBSCRIPT italic_d start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_POSTSUBSCRIPT = - divide start_ARG 1 end_ARG start_ARG 2 end_ARG end_ARG ⟩ start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT
=\displaystyle== (tdp,z)32(Δp;2z+J¯H)312|t0|bVpz;[(1+r+α)𝒩|sd1=12,sc=12t1𝒩|sd1=12,sc=12t]|spz=12|0b2superscriptsubscript𝑡𝑑𝑝𝑧32superscriptsuperscriptsubscriptΔ𝑝2𝑧subscript¯𝐽𝐻3subscriptbra12𝑡subscriptbra0𝑏𝑉superscriptsubscript𝑝𝑧delimited-[]1𝑟𝛼𝒩subscriptketformulae-sequencesubscript𝑠subscript𝑑112subscript𝑠𝑐12𝑡1𝒩subscriptketformulae-sequencesubscript𝑠subscript𝑑112subscript𝑠𝑐12𝑡ketsubscript𝑠subscript𝑝𝑧12subscriptket0𝑏2\displaystyle\frac{-(t_{dp,z})^{3}}{2(\Delta_{p;2}^{z}+\overline{J}_{H})^{3}}% \bra{-\frac{1}{2}}_{t}\bra{0}_{b}Vp_{z;\downarrow}^{\dagger}[\frac{(1+r+\alpha% )}{\mathcal{N}}\ket{s_{d_{1}}=\frac{1}{2},s_{c}=-\frac{1}{2}}_{t}-\frac{1}{% \mathcal{N}}\ket{s_{d_{1}}=-\frac{1}{2},s_{c}=\frac{1}{2}}_{t}]\ket{s_{p_{z}}=% -\frac{1}{2}}\frac{\ket{0}_{b}}{\sqrt{2}}divide start_ARG - ( italic_t start_POSTSUBSCRIPT italic_d italic_p , italic_z end_POSTSUBSCRIPT ) start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT end_ARG start_ARG 2 ( roman_Δ start_POSTSUBSCRIPT italic_p ; 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_z end_POSTSUPERSCRIPT + over¯ start_ARG italic_J end_ARG start_POSTSUBSCRIPT italic_H end_POSTSUBSCRIPT ) start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT end_ARG ⟨ start_ARG - divide start_ARG 1 end_ARG start_ARG 2 end_ARG end_ARG | start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT ⟨ start_ARG 0 end_ARG | start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT italic_V italic_p start_POSTSUBSCRIPT italic_z ; ↓ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT [ divide start_ARG ( 1 + italic_r + italic_α ) end_ARG start_ARG caligraphic_N end_ARG | start_ARG italic_s start_POSTSUBSCRIPT italic_d start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_POSTSUBSCRIPT = divide start_ARG 1 end_ARG start_ARG 2 end_ARG , italic_s start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT = - divide start_ARG 1 end_ARG start_ARG 2 end_ARG end_ARG ⟩ start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT - divide start_ARG 1 end_ARG start_ARG caligraphic_N end_ARG | start_ARG italic_s start_POSTSUBSCRIPT italic_d start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_POSTSUBSCRIPT = - divide start_ARG 1 end_ARG start_ARG 2 end_ARG , italic_s start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT = divide start_ARG 1 end_ARG start_ARG 2 end_ARG end_ARG ⟩ start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT ] | start_ARG italic_s start_POSTSUBSCRIPT italic_p start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT end_POSTSUBSCRIPT = - divide start_ARG 1 end_ARG start_ARG 2 end_ARG end_ARG ⟩ divide start_ARG | start_ARG 0 end_ARG ⟩ start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT end_ARG start_ARG square-root start_ARG 2 end_ARG end_ARG
=\displaystyle== (tdp,z)42(Δp;2z+J¯H)32(1+r+α)𝒩2=(C)=(E)=(F),superscriptsubscript𝑡𝑑𝑝𝑧42superscriptsuperscriptsubscriptΔ𝑝2𝑧subscript¯𝐽𝐻321𝑟𝛼superscript𝒩2𝐶𝐸𝐹\displaystyle\frac{(t_{dp,z})^{4}}{2(\Delta_{p;2}^{z}+\overline{J}_{H})^{3}}% \frac{\sqrt{2}(1+r+\alpha)}{\mathcal{N}^{2}}=(C)=(E)=(F),divide start_ARG ( italic_t start_POSTSUBSCRIPT italic_d italic_p , italic_z end_POSTSUBSCRIPT ) start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT end_ARG start_ARG 2 ( roman_Δ start_POSTSUBSCRIPT italic_p ; 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_z end_POSTSUPERSCRIPT + over¯ start_ARG italic_J end_ARG start_POSTSUBSCRIPT italic_H end_POSTSUBSCRIPT ) start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT end_ARG divide start_ARG square-root start_ARG 2 end_ARG ( 1 + italic_r + italic_α ) end_ARG start_ARG caligraphic_N start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG = ( italic_C ) = ( italic_E ) = ( italic_F ) ,

with α=r2+r+1𝛼superscript𝑟2𝑟1\alpha=\sqrt{r^{2}+r+1}italic_α = square-root start_ARG italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + italic_r + 1 end_ARG and 𝒩=2[(r+α)2+(r+α+1)]𝒩2delimited-[]superscript𝑟𝛼2𝑟𝛼1\mathcal{N}=\sqrt{2[(r+\alpha)^{2}+(r+\alpha+1)]}caligraphic_N = square-root start_ARG 2 [ ( italic_r + italic_α ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + ( italic_r + italic_α + 1 ) ] end_ARG. By putting everything together, Jsdsuperscriptsubscript𝐽𝑠𝑑perpendicular-toJ_{sd}^{\perp}italic_J start_POSTSUBSCRIPT italic_s italic_d end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ⟂ end_POSTSUPERSCRIPT is obtained as,

Jsdsuperscriptsubscript𝐽𝑠𝑑perpendicular-to\displaystyle J_{sd}^{\perp}italic_J start_POSTSUBSCRIPT italic_s italic_d end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ⟂ end_POSTSUPERSCRIPT =\displaystyle== 𝒞(r)tdp,z4(Δp;2z+J¯H)2[1U2+2JH+1Δp;2z+J¯H],𝒞𝑟superscriptsubscript𝑡𝑑𝑝𝑧4superscriptsuperscriptsubscriptΔ𝑝2𝑧subscript¯𝐽𝐻2delimited-[]1subscript𝑈22subscript𝐽𝐻1superscriptsubscriptΔ𝑝2𝑧subscript¯𝐽𝐻\displaystyle\mathcal{C}(r)\frac{t_{dp,z}^{4}}{(\Delta_{p;2}^{z}+\overline{J}_% {H})^{2}}\left[\frac{1}{U_{2}+2J_{H}}+\frac{1}{\Delta_{p;2}^{z}+\overline{J}_{% H}}\right],caligraphic_C ( italic_r ) divide start_ARG italic_t start_POSTSUBSCRIPT italic_d italic_p , italic_z end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT end_ARG start_ARG ( roman_Δ start_POSTSUBSCRIPT italic_p ; 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_z end_POSTSUPERSCRIPT + over¯ start_ARG italic_J end_ARG start_POSTSUBSCRIPT italic_H end_POSTSUBSCRIPT ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG [ divide start_ARG 1 end_ARG start_ARG italic_U start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT + 2 italic_J start_POSTSUBSCRIPT italic_H end_POSTSUBSCRIPT end_ARG + divide start_ARG 1 end_ARG start_ARG roman_Δ start_POSTSUBSCRIPT italic_p ; 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_z end_POSTSUPERSCRIPT + over¯ start_ARG italic_J end_ARG start_POSTSUBSCRIPT italic_H end_POSTSUBSCRIPT end_ARG ] , (S9)

with 𝒞(r)=4(1+r+α)/𝒩2𝒞𝑟41𝑟𝛼superscript𝒩2\mathcal{C}(r)=4(1+r+\alpha)/\mathcal{N}^{2}caligraphic_C ( italic_r ) = 4 ( 1 + italic_r + italic_α ) / caligraphic_N start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT, whose r𝑟ritalic_r dependence is illustrated in Fig. S3 (b). Similarly, we can show that

Jddsuperscriptsubscript𝐽𝑑𝑑perpendicular-to\displaystyle J_{dd}^{\perp}italic_J start_POSTSUBSCRIPT italic_d italic_d end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ⟂ end_POSTSUPERSCRIPT =\displaystyle== tdp,z4(Δp;2z+J¯H)2[1U2+2JH+1Δp;2z+J¯H],superscriptsubscript𝑡𝑑𝑝𝑧4superscriptsuperscriptsubscriptΔ𝑝2𝑧subscript¯𝐽𝐻2delimited-[]1subscript𝑈22subscript𝐽𝐻1superscriptsubscriptΔ𝑝2𝑧subscript¯𝐽𝐻\displaystyle\frac{t_{dp,z}^{4}}{(\Delta_{p;2}^{z}+\overline{J}_{H})^{2}}\left% [\frac{1}{U_{2}+2J_{H}}+\frac{1}{\Delta_{p;2}^{z}+\overline{J}_{H}}\right],divide start_ARG italic_t start_POSTSUBSCRIPT italic_d italic_p , italic_z end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT end_ARG start_ARG ( roman_Δ start_POSTSUBSCRIPT italic_p ; 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_z end_POSTSUPERSCRIPT + over¯ start_ARG italic_J end_ARG start_POSTSUBSCRIPT italic_H end_POSTSUBSCRIPT ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG [ divide start_ARG 1 end_ARG start_ARG italic_U start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT + 2 italic_J start_POSTSUBSCRIPT italic_H end_POSTSUBSCRIPT end_ARG + divide start_ARG 1 end_ARG start_ARG roman_Δ start_POSTSUBSCRIPT italic_p ; 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_z end_POSTSUPERSCRIPT + over¯ start_ARG italic_J end_ARG start_POSTSUBSCRIPT italic_H end_POSTSUBSCRIPT end_ARG ] , (S10)
Jsssuperscriptsubscript𝐽𝑠𝑠perpendicular-to\displaystyle J_{ss}^{\perp}italic_J start_POSTSUBSCRIPT italic_s italic_s end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ⟂ end_POSTSUPERSCRIPT =\displaystyle== 𝒞(r)2tdp,z4(Δp;2z+J¯H)2[1U2+2JH+1Δp;2z+J¯H],𝒞superscript𝑟2superscriptsubscript𝑡𝑑𝑝𝑧4superscriptsuperscriptsubscriptΔ𝑝2𝑧subscript¯𝐽𝐻2delimited-[]1subscript𝑈22subscript𝐽𝐻1superscriptsubscriptΔ𝑝2𝑧subscript¯𝐽𝐻\displaystyle\mathcal{C}(r)^{2}\frac{t_{dp,z}^{4}}{(\Delta_{p;2}^{z}+\overline% {J}_{H})^{2}}\left[\frac{1}{U_{2}+2J_{H}}+\frac{1}{\Delta_{p;2}^{z}+\overline{% J}_{H}}\right],caligraphic_C ( italic_r ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT divide start_ARG italic_t start_POSTSUBSCRIPT italic_d italic_p , italic_z end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT end_ARG start_ARG ( roman_Δ start_POSTSUBSCRIPT italic_p ; 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_z end_POSTSUPERSCRIPT + over¯ start_ARG italic_J end_ARG start_POSTSUBSCRIPT italic_H end_POSTSUBSCRIPT ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG [ divide start_ARG 1 end_ARG start_ARG italic_U start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT + 2 italic_J start_POSTSUBSCRIPT italic_H end_POSTSUBSCRIPT end_ARG + divide start_ARG 1 end_ARG start_ARG roman_Δ start_POSTSUBSCRIPT italic_p ; 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_z end_POSTSUPERSCRIPT + over¯ start_ARG italic_J end_ARG start_POSTSUBSCRIPT italic_H end_POSTSUBSCRIPT end_ARG ] , (S11)

where we have Jss=𝒞(r)Jsd=𝒞(r)2Jddsuperscriptsubscript𝐽𝑠𝑠perpendicular-to𝒞𝑟superscriptsubscript𝐽𝑠𝑑perpendicular-to𝒞superscript𝑟2superscriptsubscript𝐽𝑑𝑑perpendicular-toJ_{ss}^{\perp}=\mathcal{C}(r)J_{sd}^{\perp}=\mathcal{C}(r)^{2}J_{dd}^{\perp}italic_J start_POSTSUBSCRIPT italic_s italic_s end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ⟂ end_POSTSUPERSCRIPT = caligraphic_C ( italic_r ) italic_J start_POSTSUBSCRIPT italic_s italic_d end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ⟂ end_POSTSUPERSCRIPT = caligraphic_C ( italic_r ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_J start_POSTSUBSCRIPT italic_d italic_d end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ⟂ end_POSTSUPERSCRIPT. Here we use J¯H=JHUsubscript¯𝐽𝐻subscript𝐽𝐻superscript𝑈\overline{J}_{H}=J_{H}-U^{\prime}over¯ start_ARG italic_J end_ARG start_POSTSUBSCRIPT italic_H end_POSTSUBSCRIPT = italic_J start_POSTSUBSCRIPT italic_H end_POSTSUBSCRIPT - italic_U start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT, Δp;a=ϵpϵd;asubscriptΔ𝑝𝑎subscriptitalic-ϵ𝑝subscriptitalic-ϵ𝑑𝑎\Delta_{p;a}=\epsilon_{p}-\epsilon_{d;a}roman_Δ start_POSTSUBSCRIPT italic_p ; italic_a end_POSTSUBSCRIPT = italic_ϵ start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT - italic_ϵ start_POSTSUBSCRIPT italic_d ; italic_a end_POSTSUBSCRIPT, and Δp;az=ϵp;zϵd;asuperscriptsubscriptΔ𝑝𝑎𝑧subscriptitalic-ϵ𝑝𝑧subscriptitalic-ϵ𝑑𝑎\Delta_{p;a}^{z}=\epsilon_{p;z}-\epsilon_{d;a}roman_Δ start_POSTSUBSCRIPT italic_p ; italic_a end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_z end_POSTSUPERSCRIPT = italic_ϵ start_POSTSUBSCRIPT italic_p ; italic_z end_POSTSUBSCRIPT - italic_ϵ start_POSTSUBSCRIPT italic_d ; italic_a end_POSTSUBSCRIPT.

Refer to caption
Figure S3: (a) All possible processes contributing to the spin exchange interaction upon the fourth-order perturbation. The green (yellow) sphere represents an oxygen (nickel) atom. In (i,a,σ)𝑖𝑎𝜎(i,a,\sigma)( italic_i , italic_a , italic_σ ), i𝑖iitalic_i denotes the order of the process, a=1,2𝑎12a=1,2italic_a = 1 , 2 is for labeling dasubscript𝑑𝑎d_{a}italic_d start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT orbitals and σ𝜎\sigmaitalic_σ is for the spin. These processes specifically contribute to the sLsR+subscriptsuperscript𝑠𝐿subscriptsuperscript𝑠𝑅s^{-}_{L}s^{+}_{R}italic_s start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_L end_POSTSUBSCRIPT italic_s start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_R end_POSTSUBSCRIPT, sLSR+subscriptsuperscript𝑠𝐿subscriptsuperscript𝑆𝑅s^{-}_{L}S^{+}_{R}italic_s start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_L end_POSTSUBSCRIPT italic_S start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_R end_POSTSUBSCRIPT or SLSR+subscriptsuperscript𝑆𝐿subscriptsuperscript𝑆𝑅S^{-}_{L}S^{+}_{R}italic_S start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_L end_POSTSUBSCRIPT italic_S start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_R end_POSTSUBSCRIPT where L(R)𝐿𝑅L(R)italic_L ( italic_R ) denotes the left (right) nickel atom. (b) The JK/JH=rsubscript𝐽𝐾subscript𝐽𝐻𝑟J_{K}/J_{H}=ritalic_J start_POSTSUBSCRIPT italic_K end_POSTSUBSCRIPT / italic_J start_POSTSUBSCRIPT italic_H end_POSTSUBSCRIPT = italic_r dependence of 𝒞(r)𝒞𝑟\mathcal{C}(r)caligraphic_C ( italic_r ) in Eqs.(S9-S11). In particular, 𝒞(r)𝒞𝑟\mathcal{C}(r)caligraphic_C ( italic_r ) becomes 4343\frac{4}{3}divide start_ARG 4 end_ARG start_ARG 3 end_ARG in the r0𝑟0r\rightarrow 0italic_r → 0 limit.

III.2 B. Derivation of Jsubscript𝐽parallel-toJ_{\parallel}italic_J start_POSTSUBSCRIPT ∥ end_POSTSUBSCRIPT

Next, we show the derivation of Jsuperscript𝐽parallel-toJ^{\parallel}italic_J start_POSTSUPERSCRIPT ∥ end_POSTSUPERSCRIPT terms. There are two significant differences, compared to Jsuperscript𝐽perpendicular-toJ^{\perp}italic_J start_POSTSUPERSCRIPT ⟂ end_POSTSUPERSCRIPT. First, there are two kinds of hopping processes with tdp;1subscript𝑡𝑑𝑝1t_{dp;1}italic_t start_POSTSUBSCRIPT italic_d italic_p ; 1 end_POSTSUBSCRIPT and tdp;2subscript𝑡𝑑𝑝2t_{dp;2}italic_t start_POSTSUBSCRIPT italic_d italic_p ; 2 end_POSTSUBSCRIPT. Second, the p orbital of the Zhang-Rice spin-half state can be affected by the hoppings, since the intermediate-occupied oxygen is also shared by those states. For example, the hopping process consisting of px;i+x^2;σsuperscriptsubscript𝑝𝑥𝑖^𝑥2𝜎p_{x;i+\frac{{\widehat{x}}}{2};\sigma}^{\dagger}italic_p start_POSTSUBSCRIPT italic_x ; italic_i + divide start_ARG over^ start_ARG italic_x end_ARG end_ARG start_ARG 2 end_ARG ; italic_σ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT eliminates the piece of the same component in ϕi;σ=jB(ij)p1;j;σ=ci;σsubscriptitalic-ϕ𝑖𝜎subscript𝑗𝐵𝑖𝑗subscript𝑝1𝑗𝜎subscriptsuperscript𝑐𝑖𝜎\phi_{i;\sigma}=\sum_{j}B(i-j)p_{1;j;\sigma}=c^{\dagger}_{i;\sigma}italic_ϕ start_POSTSUBSCRIPT italic_i ; italic_σ end_POSTSUBSCRIPT = ∑ start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT italic_B ( italic_i - italic_j ) italic_p start_POSTSUBSCRIPT 1 ; italic_j ; italic_σ end_POSTSUBSCRIPT = italic_c start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_i ; italic_σ end_POSTSUBSCRIPT, where B(ij)=1Nkeik(ij)β(k)𝐵𝑖𝑗1𝑁subscript𝑘superscript𝑒𝑖𝑘𝑖𝑗𝛽𝑘B(i-j)=\frac{1}{N}\sum_{k}e^{ik\cdot(i-j)}\beta(k)italic_B ( italic_i - italic_j ) = divide start_ARG 1 end_ARG start_ARG italic_N end_ARG ∑ start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT italic_e start_POSTSUPERSCRIPT italic_i italic_k ⋅ ( italic_i - italic_j ) end_POSTSUPERSCRIPT italic_β ( italic_k ) and β(k)=[112(coskx+cosky)]1/2𝛽𝑘superscriptdelimited-[]112subscript𝑘𝑥subscript𝑘𝑦12\beta(k)=[1-\frac{1}{2}(\cos k_{x}+\cos k_{y})]^{-1/2}italic_β ( italic_k ) = [ 1 - divide start_ARG 1 end_ARG start_ARG 2 end_ARG ( roman_cos italic_k start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT + roman_cos italic_k start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT ) ] start_POSTSUPERSCRIPT - 1 / 2 end_POSTSUPERSCRIPT. To be specific, consider the following example,

px;i+x^2;σϕi,σsuperscriptsubscript𝑝𝑥𝑖^𝑥2𝜎superscriptsubscriptitalic-ϕ𝑖superscript𝜎\displaystyle p_{x;i+\frac{{\widehat{x}}}{2};\sigma}^{\dagger}\phi_{i,\sigma^{% \prime}}^{\dagger}italic_p start_POSTSUBSCRIPT italic_x ; italic_i + divide start_ARG over^ start_ARG italic_x end_ARG end_ARG start_ARG 2 end_ARG ; italic_σ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT italic_ϕ start_POSTSUBSCRIPT italic_i , italic_σ start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT =\displaystyle== px;i+x^2;σ[j{i,i+x^}cB(ij)p1;j,σ+j{i,i+x^}B(ij)p¯1;j,σ]=px;i+x^2;σϕ¯i,σ,superscriptsubscript𝑝𝑥𝑖^𝑥2𝜎delimited-[]subscript𝑗superscript𝑖𝑖^𝑥𝑐𝐵𝑖𝑗superscriptsubscript𝑝1𝑗superscript𝜎subscript𝑗𝑖𝑖^𝑥𝐵𝑖𝑗superscriptsubscript¯𝑝1𝑗superscript𝜎superscriptsubscript𝑝𝑥𝑖^𝑥2𝜎superscriptsubscript¯italic-ϕ𝑖superscript𝜎\displaystyle p_{x;i+\frac{{\widehat{x}}}{2};\sigma}^{\dagger}\left[\sum_{j\in% \{i,i+{\widehat{x}}\}^{c}}B(i-j)p_{1;j,\sigma^{\prime}}^{\dagger}+\sum_{j\in\{% i,i+{\widehat{x}}\}}B(i-j)\overline{p}_{1;j,\sigma^{\prime}}^{\dagger}\right]=% p_{x;i+\frac{{\widehat{x}}}{2};\sigma}^{\dagger}\overline{\phi}_{i,\sigma^{% \prime}}^{\dagger},italic_p start_POSTSUBSCRIPT italic_x ; italic_i + divide start_ARG over^ start_ARG italic_x end_ARG end_ARG start_ARG 2 end_ARG ; italic_σ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT [ ∑ start_POSTSUBSCRIPT italic_j ∈ { italic_i , italic_i + over^ start_ARG italic_x end_ARG } start_POSTSUPERSCRIPT italic_c end_POSTSUPERSCRIPT end_POSTSUBSCRIPT italic_B ( italic_i - italic_j ) italic_p start_POSTSUBSCRIPT 1 ; italic_j , italic_σ start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT + ∑ start_POSTSUBSCRIPT italic_j ∈ { italic_i , italic_i + over^ start_ARG italic_x end_ARG } end_POSTSUBSCRIPT italic_B ( italic_i - italic_j ) over¯ start_ARG italic_p end_ARG start_POSTSUBSCRIPT 1 ; italic_j , italic_σ start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT ] = italic_p start_POSTSUBSCRIPT italic_x ; italic_i + divide start_ARG over^ start_ARG italic_x end_ARG end_ARG start_ARG 2 end_ARG ; italic_σ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT over¯ start_ARG italic_ϕ end_ARG start_POSTSUBSCRIPT italic_i , italic_σ start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT ,

where we define the ϕ¯(p¯)¯italic-ϕ¯𝑝\overline{\phi}(\overline{p})over¯ start_ARG italic_ϕ end_ARG ( over¯ start_ARG italic_p end_ARG ) as the ϕ(p)italic-ϕ𝑝\phi(p)italic_ϕ ( italic_p ), but excluding the px;i+x^/2subscript𝑝𝑥𝑖^𝑥2p_{x;i+{\widehat{x}}/2}italic_p start_POSTSUBSCRIPT italic_x ; italic_i + over^ start_ARG italic_x end_ARG / 2 end_POSTSUBSCRIPT piece,

p¯1;i,σsuperscriptsubscript¯𝑝1𝑖superscript𝜎\displaystyle\overline{p}_{1;i,\sigma^{\prime}}^{\dagger}over¯ start_ARG italic_p end_ARG start_POSTSUBSCRIPT 1 ; italic_i , italic_σ start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT =\displaystyle== 12[px;ix^2;l;σpy;i+y^2;l;σ+py;iy^2;l;σ]=p1;i,σ12px;i+x^2;l;σ,12delimited-[]subscript𝑝𝑥𝑖^𝑥2𝑙𝜎subscript𝑝𝑦𝑖^𝑦2𝑙𝜎subscript𝑝𝑦𝑖^𝑦2𝑙𝜎superscriptsubscript𝑝1𝑖superscript𝜎12subscript𝑝𝑥𝑖^𝑥2𝑙𝜎\displaystyle\frac{1}{2}\left[-p_{x;i-\frac{{\widehat{x}}}{2};l;\sigma}-p_{y;i% +\frac{{\widehat{y}}}{2};l;\sigma}+p_{y;i-\frac{{\widehat{y}}}{2};l;\sigma}% \right]=p_{1;i,\sigma^{\prime}}^{\dagger}-\frac{1}{2}p_{x;i+\frac{{\widehat{x}% }}{2};l;\sigma},divide start_ARG 1 end_ARG start_ARG 2 end_ARG [ - italic_p start_POSTSUBSCRIPT italic_x ; italic_i - divide start_ARG over^ start_ARG italic_x end_ARG end_ARG start_ARG 2 end_ARG ; italic_l ; italic_σ end_POSTSUBSCRIPT - italic_p start_POSTSUBSCRIPT italic_y ; italic_i + divide start_ARG over^ start_ARG italic_y end_ARG end_ARG start_ARG 2 end_ARG ; italic_l ; italic_σ end_POSTSUBSCRIPT + italic_p start_POSTSUBSCRIPT italic_y ; italic_i - divide start_ARG over^ start_ARG italic_y end_ARG end_ARG start_ARG 2 end_ARG ; italic_l ; italic_σ end_POSTSUBSCRIPT ] = italic_p start_POSTSUBSCRIPT 1 ; italic_i , italic_σ start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT - divide start_ARG 1 end_ARG start_ARG 2 end_ARG italic_p start_POSTSUBSCRIPT italic_x ; italic_i + divide start_ARG over^ start_ARG italic_x end_ARG end_ARG start_ARG 2 end_ARG ; italic_l ; italic_σ end_POSTSUBSCRIPT ,
p¯1;i+x^,σsuperscriptsubscript¯𝑝1𝑖^𝑥superscript𝜎\displaystyle\overline{p}_{1;i+{\widehat{x}},\sigma^{\prime}}^{\dagger}over¯ start_ARG italic_p end_ARG start_POSTSUBSCRIPT 1 ; italic_i + over^ start_ARG italic_x end_ARG , italic_σ start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT =\displaystyle== 12[px;i+3x^2;l;σpy;i+x^+y^2;l;σ+py;i+x^y^2;l;σ]=p1;i,σ+12px;i+x^2;l;σ,12delimited-[]subscript𝑝𝑥𝑖3^𝑥2𝑙𝜎subscript𝑝𝑦𝑖^𝑥^𝑦2𝑙𝜎subscript𝑝𝑦𝑖^𝑥^𝑦2𝑙𝜎superscriptsubscript𝑝1𝑖superscript𝜎12subscript𝑝𝑥𝑖^𝑥2𝑙𝜎\displaystyle\frac{1}{2}\left[p_{x;i+\frac{3{\widehat{x}}}{2};l;\sigma}-p_{y;i% +{\widehat{x}}+\frac{{\widehat{y}}}{2};l;\sigma}+p_{y;i+{\widehat{x}}-\frac{{% \widehat{y}}}{2};l;\sigma}\right]=p_{1;i,\sigma^{\prime}}^{\dagger}+\frac{1}{2% }p_{x;i+\frac{{\widehat{x}}}{2};l;\sigma},divide start_ARG 1 end_ARG start_ARG 2 end_ARG [ italic_p start_POSTSUBSCRIPT italic_x ; italic_i + divide start_ARG 3 over^ start_ARG italic_x end_ARG end_ARG start_ARG 2 end_ARG ; italic_l ; italic_σ end_POSTSUBSCRIPT - italic_p start_POSTSUBSCRIPT italic_y ; italic_i + over^ start_ARG italic_x end_ARG + divide start_ARG over^ start_ARG italic_y end_ARG end_ARG start_ARG 2 end_ARG ; italic_l ; italic_σ end_POSTSUBSCRIPT + italic_p start_POSTSUBSCRIPT italic_y ; italic_i + over^ start_ARG italic_x end_ARG - divide start_ARG over^ start_ARG italic_y end_ARG end_ARG start_ARG 2 end_ARG ; italic_l ; italic_σ end_POSTSUBSCRIPT ] = italic_p start_POSTSUBSCRIPT 1 ; italic_i , italic_σ start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT + divide start_ARG 1 end_ARG start_ARG 2 end_ARG italic_p start_POSTSUBSCRIPT italic_x ; italic_i + divide start_ARG over^ start_ARG italic_x end_ARG end_ARG start_ARG 2 end_ARG ; italic_l ; italic_σ end_POSTSUBSCRIPT ,

and

ϕ¯i,σ=ϕi,σ[B(0)B(1)2]px;i+x2;σ.superscriptsubscript¯italic-ϕ𝑖superscript𝜎superscriptsubscriptitalic-ϕ𝑖superscript𝜎delimited-[]𝐵0𝐵12superscriptsubscript𝑝𝑥𝑖𝑥2𝜎\displaystyle\overline{\phi}_{i,\sigma^{\prime}}^{\dagger}=\phi_{i,\sigma^{% \prime}}^{\dagger}-[\frac{B(0)-B(1)}{2}]p_{x;i+\frac{x}{2};\sigma}^{\dagger}.over¯ start_ARG italic_ϕ end_ARG start_POSTSUBSCRIPT italic_i , italic_σ start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT = italic_ϕ start_POSTSUBSCRIPT italic_i , italic_σ start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT - [ divide start_ARG italic_B ( 0 ) - italic_B ( 1 ) end_ARG start_ARG 2 end_ARG ] italic_p start_POSTSUBSCRIPT italic_x ; italic_i + divide start_ARG italic_x end_ARG start_ARG 2 end_ARG ; italic_σ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT .

Note that ϕ¯|ϕ1[B(0)B(1)2]20.7705𝒜\bra{\overline{\phi}}\phi\rangle\simeq 1-[\frac{B(0)-B(1)}{2}]^{2}\simeq 0.770% 5\equiv\mathcal{A}⟨ start_ARG over¯ start_ARG italic_ϕ end_ARG end_ARG | italic_ϕ ⟩ ≃ 1 - [ divide start_ARG italic_B ( 0 ) - italic_B ( 1 ) end_ARG start_ARG 2 end_ARG ] start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ≃ 0.7705 ≡ caligraphic_A. Then, the next step is straightforward to consider all the possible processes depicted in Fig S3 (a). After finishing all the calculations with Δp;1=Δp;2subscriptΔ𝑝1subscriptΔ𝑝2\Delta_{p;1}=\Delta_{p;2}roman_Δ start_POSTSUBSCRIPT italic_p ; 1 end_POSTSUBSCRIPT = roman_Δ start_POSTSUBSCRIPT italic_p ; 2 end_POSTSUBSCRIPT, we obtained that

Jddsuperscriptsubscript𝐽𝑑𝑑parallel-to\displaystyle J_{dd}^{\parallel}italic_J start_POSTSUBSCRIPT italic_d italic_d end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ∥ end_POSTSUPERSCRIPT =\displaystyle== a,atdp;a2tdp;a2(Δp+J¯H)2[1Ua+2JH+1Δp+J¯H],subscript𝑎superscript𝑎superscriptsubscript𝑡𝑑𝑝𝑎2superscriptsubscript𝑡𝑑𝑝superscript𝑎2superscriptsubscriptΔ𝑝subscript¯𝐽𝐻2delimited-[]1subscript𝑈superscript𝑎2subscript𝐽𝐻1subscriptΔ𝑝subscript¯𝐽𝐻\displaystyle\sum_{a,a^{\prime}}\frac{t_{dp;a}^{2}t_{dp;a^{\prime}}^{2}}{(% \Delta_{p}+\overline{J}_{H})^{2}}\left[\frac{1}{U_{a^{\prime}}+2J_{H}}+\frac{1% }{\Delta_{p}+\overline{J}_{H}}\right],∑ start_POSTSUBSCRIPT italic_a , italic_a start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT divide start_ARG italic_t start_POSTSUBSCRIPT italic_d italic_p ; italic_a end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_t start_POSTSUBSCRIPT italic_d italic_p ; italic_a start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG ( roman_Δ start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT + over¯ start_ARG italic_J end_ARG start_POSTSUBSCRIPT italic_H end_POSTSUBSCRIPT ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG [ divide start_ARG 1 end_ARG start_ARG italic_U start_POSTSUBSCRIPT italic_a start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT + 2 italic_J start_POSTSUBSCRIPT italic_H end_POSTSUBSCRIPT end_ARG + divide start_ARG 1 end_ARG start_ARG roman_Δ start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT + over¯ start_ARG italic_J end_ARG start_POSTSUBSCRIPT italic_H end_POSTSUBSCRIPT end_ARG ] , (S12)
Jsdsuperscriptsubscript𝐽𝑠𝑑parallel-to\displaystyle J_{sd}^{\parallel}italic_J start_POSTSUBSCRIPT italic_s italic_d end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ∥ end_POSTSUPERSCRIPT =\displaystyle== 𝒞(r)𝒜a,atdp;a2tdp;a2(Δp+J¯H)2[1Ua+2JH+1Δp+J¯H],𝒞𝑟𝒜subscript𝑎superscript𝑎superscriptsubscript𝑡𝑑𝑝𝑎2superscriptsubscript𝑡𝑑𝑝superscript𝑎2superscriptsubscriptΔ𝑝subscript¯𝐽𝐻2delimited-[]1subscript𝑈superscript𝑎2subscript𝐽𝐻1subscriptΔ𝑝subscript¯𝐽𝐻\displaystyle\mathcal{C}(r)\mathcal{A}\sum_{a,a^{\prime}}\frac{t_{dp;a}^{2}t_{% dp;a^{\prime}}^{2}}{(\Delta_{p}+\overline{J}_{H})^{2}}\left[\frac{1}{U_{a^{% \prime}}+2J_{H}}+\frac{1}{\Delta_{p}+\overline{J}_{H}}\right],caligraphic_C ( italic_r ) caligraphic_A ∑ start_POSTSUBSCRIPT italic_a , italic_a start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT divide start_ARG italic_t start_POSTSUBSCRIPT italic_d italic_p ; italic_a end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_t start_POSTSUBSCRIPT italic_d italic_p ; italic_a start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG ( roman_Δ start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT + over¯ start_ARG italic_J end_ARG start_POSTSUBSCRIPT italic_H end_POSTSUBSCRIPT ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG [ divide start_ARG 1 end_ARG start_ARG italic_U start_POSTSUBSCRIPT italic_a start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT + 2 italic_J start_POSTSUBSCRIPT italic_H end_POSTSUBSCRIPT end_ARG + divide start_ARG 1 end_ARG start_ARG roman_Δ start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT + over¯ start_ARG italic_J end_ARG start_POSTSUBSCRIPT italic_H end_POSTSUBSCRIPT end_ARG ] , (S13)
Jsssuperscriptsubscript𝐽𝑠𝑠parallel-to\displaystyle J_{ss}^{\parallel}italic_J start_POSTSUBSCRIPT italic_s italic_s end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ∥ end_POSTSUPERSCRIPT =\displaystyle== (𝒞(r)𝒜)2a,atdp;a2tdp;a2(Δp+J¯H)2[1Ua+2JH+1Δp+J¯H]superscript𝒞𝑟𝒜2subscript𝑎superscript𝑎superscriptsubscript𝑡𝑑𝑝𝑎2superscriptsubscript𝑡𝑑𝑝superscript𝑎2superscriptsubscriptΔ𝑝subscript¯𝐽𝐻2delimited-[]1subscript𝑈superscript𝑎2subscript𝐽𝐻1subscriptΔ𝑝subscript¯𝐽𝐻\displaystyle(\mathcal{C}(r)\mathcal{A})^{2}\sum_{a,a^{\prime}}\frac{t_{dp;a}^% {2}t_{dp;a^{\prime}}^{2}}{(\Delta_{p}+\overline{J}_{H})^{2}}\left[\frac{1}{U_{% a^{\prime}}+2J_{H}}+\frac{1}{\Delta_{p}+\overline{J}_{H}}\right]( caligraphic_C ( italic_r ) caligraphic_A ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ∑ start_POSTSUBSCRIPT italic_a , italic_a start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT divide start_ARG italic_t start_POSTSUBSCRIPT italic_d italic_p ; italic_a end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_t start_POSTSUBSCRIPT italic_d italic_p ; italic_a start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG ( roman_Δ start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT + over¯ start_ARG italic_J end_ARG start_POSTSUBSCRIPT italic_H end_POSTSUBSCRIPT ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG [ divide start_ARG 1 end_ARG start_ARG italic_U start_POSTSUBSCRIPT italic_a start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT + 2 italic_J start_POSTSUBSCRIPT italic_H end_POSTSUBSCRIPT end_ARG + divide start_ARG 1 end_ARG start_ARG roman_Δ start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT + over¯ start_ARG italic_J end_ARG start_POSTSUBSCRIPT italic_H end_POSTSUBSCRIPT end_ARG ] (S14)

Now, we have Jss=𝒞(r)𝒜Jsd=𝒞(r)2𝒜2Jddsuperscriptsubscript𝐽𝑠𝑠parallel-to𝒞𝑟𝒜superscriptsubscript𝐽𝑠𝑑parallel-to𝒞superscript𝑟2superscript𝒜2superscriptsubscript𝐽𝑑𝑑parallel-toJ_{ss}^{\parallel}=\mathcal{C}(r)\mathcal{A}J_{sd}^{\parallel}=\mathcal{C}(r)^% {2}\mathcal{A}^{2}J_{dd}^{\parallel}italic_J start_POSTSUBSCRIPT italic_s italic_s end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ∥ end_POSTSUPERSCRIPT = caligraphic_C ( italic_r ) caligraphic_A italic_J start_POSTSUBSCRIPT italic_s italic_d end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ∥ end_POSTSUPERSCRIPT = caligraphic_C ( italic_r ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT caligraphic_A start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_J start_POSTSUBSCRIPT italic_d italic_d end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ∥ end_POSTSUPERSCRIPT. We stress this condition is generally satisfied without imposing Δp;1=Δp;2subscriptΔ𝑝1subscriptΔ𝑝2\Delta_{p;1}=\Delta_{p;2}roman_Δ start_POSTSUBSCRIPT italic_p ; 1 end_POSTSUBSCRIPT = roman_Δ start_POSTSUBSCRIPT italic_p ; 2 end_POSTSUBSCRIPT.

III.3 C. Derivation of tsubscript𝑡parallel-tot_{\parallel}italic_t start_POSTSUBSCRIPT ∥ end_POSTSUBSCRIPT

From the second order perturbation, we can derive the in-plane hopping term with and without spin dependence l,σ,i,jt;i,jci;l;σcj;l,σ\sum_{l,\sigma,i,j}-t_{\parallel;i,j}c^{\dagger}_{i;l;\sigma}c_{j;l,\sigma}∑ start_POSTSUBSCRIPT italic_l , italic_σ , italic_i , italic_j end_POSTSUBSCRIPT - italic_t start_POSTSUBSCRIPT ∥ ; italic_i , italic_j end_POSTSUBSCRIPT italic_c start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_i ; italic_l ; italic_σ end_POSTSUBSCRIPT italic_c start_POSTSUBSCRIPT italic_j ; italic_l , italic_σ end_POSTSUBSCRIPT. The main contribution can be classified as the effective O-hopping and spin exchange between the Ni and O holes, similar as the case in cuprate Zhang and Rice (1988). We have that:

t;i,j=\displaystyle t_{\parallel;i,j}=italic_t start_POSTSUBSCRIPT ∥ ; italic_i , italic_j end_POSTSUBSCRIPT = (2tdp;123(U+U+JHΔp)+JK;1122)δij,NN/𝒩2(r)4JK;1λ3N𝐤β𝐤1ei𝐤(𝐫𝐢𝐫𝐣)/𝒩4(r),2superscriptsubscript𝑡𝑑𝑝123𝑈superscript𝑈subscript𝐽𝐻subscriptΔ𝑝subscript𝐽𝐾1122subscript𝛿𝑖𝑗NNsuperscript𝒩2𝑟4subscript𝐽𝐾1𝜆3𝑁subscript𝐤subscriptsuperscript𝛽1𝐤superscript𝑒i𝐤subscript𝐫𝐢subscript𝐫𝐣superscript𝒩4𝑟\displaystyle\left(\frac{2t_{dp;1}^{2}}{3(U+U^{\prime}+J_{H}-\Delta_{p})}+% \frac{J_{K;1}}{12\sqrt{2}}\right)\delta_{ij,\mathrm{NN}}/\mathcal{N}^{2}(r)-% \frac{4J_{K;1}\lambda}{3N}\sum_{\mathbf{k}}\beta^{-1}_{\mathbf{k}}e^{-\mathrm{% i}\mathbf{k}\cdot(\mathbf{r_{i}-r_{j}})}/\mathcal{N}^{4}(r),( divide start_ARG 2 italic_t start_POSTSUBSCRIPT italic_d italic_p ; 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG 3 ( italic_U + italic_U start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT + italic_J start_POSTSUBSCRIPT italic_H end_POSTSUBSCRIPT - roman_Δ start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT ) end_ARG + divide start_ARG italic_J start_POSTSUBSCRIPT italic_K ; 1 end_POSTSUBSCRIPT end_ARG start_ARG 12 square-root start_ARG 2 end_ARG end_ARG ) italic_δ start_POSTSUBSCRIPT italic_i italic_j , roman_NN end_POSTSUBSCRIPT / caligraphic_N start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( italic_r ) - divide start_ARG 4 italic_J start_POSTSUBSCRIPT italic_K ; 1 end_POSTSUBSCRIPT italic_λ end_ARG start_ARG 3 italic_N end_ARG ∑ start_POSTSUBSCRIPT bold_k end_POSTSUBSCRIPT italic_β start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT bold_k end_POSTSUBSCRIPT italic_e start_POSTSUPERSCRIPT - roman_i bold_k ⋅ ( bold_r start_POSTSUBSCRIPT bold_i end_POSTSUBSCRIPT - bold_r start_POSTSUBSCRIPT bold_j end_POSTSUBSCRIPT ) end_POSTSUPERSCRIPT / caligraphic_N start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT ( italic_r ) , (S15)

in which δij,NNsubscript𝛿𝑖𝑗NN\delta_{ij,\mathrm{NN}}italic_δ start_POSTSUBSCRIPT italic_i italic_j , roman_NN end_POSTSUBSCRIPT means that only nearest neighbor term is non-zero. λ𝜆\lambdaitalic_λ is defined as 𝐤β𝐤10.96similar-to-or-equalssubscript𝐤subscriptsuperscript𝛽1𝐤0.96\sum_{\mathbf{k}}\beta^{-1}_{\mathbf{k}}\simeq 0.96∑ start_POSTSUBSCRIPT bold_k end_POSTSUBSCRIPT italic_β start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT bold_k end_POSTSUBSCRIPT ≃ 0.96. Here, we used 𝒩(r)=[2(r+α)2+2(r+α)+2]1/2/(r+α+1)𝒩𝑟superscriptdelimited-[]2superscript𝑟𝛼22𝑟𝛼212𝑟𝛼1\mathcal{N}(r)=[2(r+\alpha)^{2}+2(r+\alpha)+2]^{1/2}/(r+\alpha+1)caligraphic_N ( italic_r ) = [ 2 ( italic_r + italic_α ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + 2 ( italic_r + italic_α ) + 2 ] start_POSTSUPERSCRIPT 1 / 2 end_POSTSUPERSCRIPT / ( italic_r + italic_α + 1 ), a function as increasing r𝑟ritalic_r, which depends on r=JK;1/JH𝑟subscript𝐽𝐾1subscript𝐽𝐻r=J_{K;1}/J_{H}italic_r = italic_J start_POSTSUBSCRIPT italic_K ; 1 end_POSTSUBSCRIPT / italic_J start_POSTSUBSCRIPT italic_H end_POSTSUBSCRIPT. The function is monotonic increasing in the range of [32,2]322[\sqrt{\frac{3}{2}},\sqrt{2}][ square-root start_ARG divide start_ARG 3 end_ARG start_ARG 2 end_ARG end_ARG , square-root start_ARG 2 end_ARG ]. N𝑁Nitalic_N is the number of the site. In our main text, we only keep the nearest neighbor term of t;i,jt_{\parallel;i,j}italic_t start_POSTSUBSCRIPT ∥ ; italic_i , italic_j end_POSTSUBSCRIPT, and denoted by tsubscript𝑡parallel-tot_{\parallel}italic_t start_POSTSUBSCRIPT ∥ end_POSTSUBSCRIPT. The r𝑟ritalic_r dependnence is illustrated in Fig.S4.

Refer to caption
Figure S4: The r𝑟ritalic_r dependence of t(r)subscript𝑡parallel-to𝑟t_{\parallel}(r)italic_t start_POSTSUBSCRIPT ∥ end_POSTSUBSCRIPT ( italic_r ) and 𝒩(r)𝒩𝑟\mathcal{N}(r)caligraphic_N ( italic_r ) estimated by the parameters mentioned in the main-text. Note that tsubscript𝑡parallel-tot_{\parallel}italic_t start_POSTSUBSCRIPT ∥ end_POSTSUBSCRIPT is plotted in the eV unit, and 𝒩𝒩\mathcal{N}caligraphic_N is a dimensionless quantity.

IV IV. More detailed DMRG simulation results

In this section, we provide additional DMRG calculation results, which haven’t been shown in the main text. Here, we still use t=1subscript𝑡parallel-to1t_{\parallel}=1italic_t start_POSTSUBSCRIPT ∥ end_POSTSUBSCRIPT = 1 and Jss=0.1superscriptsubscript𝐽parallel-to𝑠𝑠0.1J_{\parallel}^{ss}=0.1italic_J start_POSTSUBSCRIPT ∥ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_s italic_s end_POSTSUPERSCRIPT = 0.1. In Fig. S5, we have shown the tsubscript𝑡perpendicular-tot_{\perp}italic_t start_POSTSUBSCRIPT ⟂ end_POSTSUBSCRIPT dependence with various interlayer repulsive interaction strength V𝑉Vitalic_V to demonstrate, that the non-monotonic decreasing region diminishes and eventually disappears increasing V𝑉Vitalic_V. In Fig. S6, we provide some DMRG simulation data at x=0.2𝑥0.2x=0.2italic_x = 0.2, for comparing the data provided in the main text at x=0.5𝑥0.5x=0.5italic_x = 0.5 and showing that the key characteristic behaviors are similar. In Fig. S7, we provide the convergence check under the system size of the spin-gap results through DMRG calculations.

Refer to caption
Figure S5: DMRG simulation results of the type-II t-J model with Lz=2,Ly=1formulae-sequencesubscript𝐿𝑧2subscript𝐿𝑦1L_{z}=2,L_{y}=1italic_L start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT = 2 , italic_L start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT = 1 at t=1subscript𝑡parallel-to1t_{\parallel}=1italic_t start_POSTSUBSCRIPT ∥ end_POSTSUBSCRIPT = 1, Jss=0.1superscriptsubscript𝐽parallel-to𝑠𝑠0.1J_{\parallel}^{ss}=0.1italic_J start_POSTSUBSCRIPT ∥ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_s italic_s end_POSTSUPERSCRIPT = 0.1. (a,b) The tsubscript𝑡perpendicular-tot_{\perp}italic_t start_POSTSUBSCRIPT ⟂ end_POSTSUBSCRIPT dependence of the spin gap with various V=0,0.5,1,1.5,2𝑉00.511.52V=0,0.5,1,1.5,2italic_V = 0 , 0.5 , 1 , 1.5 , 2 at (a) x=0.5𝑥0.5x=0.5italic_x = 0.5 and (b) x=0.2𝑥0.2x=0.2italic_x = 0.2. Here, we choose t=1subscript𝑡perpendicular-to1t_{\perp}=1italic_t start_POSTSUBSCRIPT ⟂ end_POSTSUBSCRIPT = 1, and use Lx=40subscript𝐿𝑥40L_{x}=40italic_L start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT = 40 and χ=2400𝜒2400\chi=2400italic_χ = 2400 for simulation. The non-monotonically decreasing region is decreasing as increasing the repulsive interaction where the system goes to the BCS limit.
Refer to caption
Figure S6: DMRG simulation results of the type-II t-J model, with Lz=2,Ly=1formulae-sequencesubscript𝐿𝑧2subscript𝐿𝑦1L_{z}=2,L_{y}=1italic_L start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT = 2 , italic_L start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT = 1 at t=1subscript𝑡parallel-to1t_{\parallel}=1italic_t start_POSTSUBSCRIPT ∥ end_POSTSUBSCRIPT = 1, Jss=0.1superscriptsubscript𝐽parallel-to𝑠𝑠0.1J_{\parallel}^{ss}=0.1italic_J start_POSTSUBSCRIPT ∥ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_s italic_s end_POSTSUPERSCRIPT = 0.1. The doping is fixed at x=0.2𝑥0.2x=0.2italic_x = 0.2 to compare the data at x=0.5𝑥0.5x=0.5italic_x = 0.5 provided in the main-text. At x=0.2𝑥0.2x=0.2italic_x = 0.2, the results exhibits the evidence of the Luther-Emily liquid phases. (a) The pair correlation function at t=0.1subscript𝑡perpendicular-to0.1t_{\perp}=0.1italic_t start_POSTSUBSCRIPT ⟂ end_POSTSUBSCRIPT = 0.1 and x=0.2𝑥0.2x=0.2italic_x = 0.2 shows the powe-law decaying behaviours. (b) The entanglement entropy and the correlation length at x=0.2𝑥0.2x=0.2italic_x = 0.2. The fitted central charge is nearly c=1𝑐1c=1italic_c = 1.
Refer to caption
Figure S7: The system size dependence of DMRG simulation results of the type-II t-J model with Lz=2,Ly=1formulae-sequencesubscript𝐿𝑧2subscript𝐿𝑦1L_{z}=2,L_{y}=1italic_L start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT = 2 , italic_L start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT = 1 at t=1subscript𝑡parallel-to1t_{\parallel}=1italic_t start_POSTSUBSCRIPT ∥ end_POSTSUBSCRIPT = 1, Jss=0.1superscriptsubscript𝐽parallel-to𝑠𝑠0.1J_{\parallel}^{ss}=0.1italic_J start_POSTSUBSCRIPT ∥ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_s italic_s end_POSTSUPERSCRIPT = 0.1 (a) tsubscript𝑡perpendicular-tot_{\perp}italic_t start_POSTSUBSCRIPT ⟂ end_POSTSUBSCRIPT dependence of the spin gap with various Jsubscript𝐽perpendicular-toJ_{\perp}italic_J start_POSTSUBSCRIPT ⟂ end_POSTSUBSCRIPT at x=0.2𝑥0.2x=0.2italic_x = 0.2. (b) Doping x𝑥xitalic_x dependence of the spin gap with various V𝑉Vitalic_V at J=1subscript𝐽perpendicular-to1J_{\perp}=1italic_J start_POSTSUBSCRIPT ⟂ end_POSTSUBSCRIPT = 1 (χ=2400𝜒2400\chi=2400italic_χ = 2400).

V V. The pressure dependence of critical temperature

In this section, we analyze the binding energy for predicting the critical temperature under changing the pressure. Applying pressure effectively increases the hopping between pzsubscript𝑝𝑧p_{z}italic_p start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT and inplane d𝑑ditalic_d orbitals. Hence, according to our energy estimations, we expect there is a crossover pressure, satisfying EG1=EGzsuperscriptsubscript𝐸𝐺1superscriptsubscript𝐸𝐺𝑧E_{G}^{1}=E_{G}^{z}italic_E start_POSTSUBSCRIPT italic_G end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 1 end_POSTSUPERSCRIPT = italic_E start_POSTSUBSCRIPT italic_G end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_z end_POSTSUPERSCRIPT, as illustrated in Fig.S8 (a). At P<P𝑃subscript𝑃P<P_{*}italic_P < italic_P start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT, the additional holes are occupied in p1subscript𝑝1p_{1}italic_p start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT orbital, while at P>P𝑃subscript𝑃P>P_{*}italic_P > italic_P start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT, the pzsubscript𝑝𝑧p_{z}italic_p start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT orbital becomes dominant.

Then, we analyze the binding energy the estimate the critical temperature. For each regime of P<P𝑃subscript𝑃P<P_{*}italic_P < italic_P start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT and P>P𝑃subscript𝑃P>P_{*}italic_P > italic_P start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT, the binding energy is differently estimated, since the additional hole is occupied by p1subscript𝑝1p_{1}italic_p start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT and pzsubscript𝑝𝑧p_{z}italic_p start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT respectively,

P<P:EB1=E[np1=1](E[np1=2]+E[np1=0])/2,\displaystyle P<P*:\quad E_{B}^{1}=E[n_{p_{1}}=1]-(E[n_{p_{1}}=2]+E[n_{p_{1}}=% 0])/2,italic_P < italic_P ∗ : italic_E start_POSTSUBSCRIPT italic_B end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 1 end_POSTSUPERSCRIPT = italic_E [ italic_n start_POSTSUBSCRIPT italic_p start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_POSTSUBSCRIPT = 1 ] - ( italic_E [ italic_n start_POSTSUBSCRIPT italic_p start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_POSTSUBSCRIPT = 2 ] + italic_E [ italic_n start_POSTSUBSCRIPT italic_p start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_POSTSUBSCRIPT = 0 ] ) / 2 ,
P>P:EBz=E[npz=1](E[npz=2]+E[npz=0])/2.\displaystyle P>P*:\quad E_{B}^{z}=E[n_{p_{z}}=1]-(E[n_{p_{z}}=2]+E[n_{p_{z}}=% 0])/2.italic_P > italic_P ∗ : italic_E start_POSTSUBSCRIPT italic_B end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_z end_POSTSUPERSCRIPT = italic_E [ italic_n start_POSTSUBSCRIPT italic_p start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT end_POSTSUBSCRIPT = 1 ] - ( italic_E [ italic_n start_POSTSUBSCRIPT italic_p start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT end_POSTSUBSCRIPT = 2 ] + italic_E [ italic_n start_POSTSUBSCRIPT italic_p start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT end_POSTSUBSCRIPT = 0 ] ) / 2 .

The pressure, (or equivalently tdp;zsubscript𝑡𝑑𝑝𝑧t_{dp;z}italic_t start_POSTSUBSCRIPT italic_d italic_p ; italic_z end_POSTSUBSCRIPT) dependence of the binding energy per unit cell, is illustrated in Fig. S8 (b). In the regime where p1subscript𝑝1p_{1}italic_p start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT dominates, the binding energy exhibits a positive trend, steadily increasing under pressure. This arises from the stabilization of the two-hole state, occupied by p1subscript𝑝1p_{1}italic_p start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT at each layer, which enhanced inter-layer spin coupling by larger tdp;zsubscript𝑡𝑑𝑝𝑧t_{dp;z}italic_t start_POSTSUBSCRIPT italic_d italic_p ; italic_z end_POSTSUBSCRIPT. Conversely, in the regime where pzsubscript𝑝𝑧p_{z}italic_p start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT dominates, the binding energy demonstrates a negative trend, decreasing under pressure. This stems from the fact that the two hole states occupied by one pzsubscript𝑝𝑧p_{z}italic_p start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT orbital blocks the hopping between two layers of nickel, potentially leading the inter-layer spin coupling as zero.

Refer to caption
Figure S8: (a) tdp;zsubscript𝑡𝑑𝑝𝑧t_{dp;z}italic_t start_POSTSUBSCRIPT italic_d italic_p ; italic_z end_POSTSUBSCRIPT dependence of EG1EGzsuperscriptsubscript𝐸𝐺1superscriptsubscript𝐸𝐺𝑧E_{G}^{1}-E_{G}^{z}italic_E start_POSTSUBSCRIPT italic_G end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 1 end_POSTSUPERSCRIPT - italic_E start_POSTSUBSCRIPT italic_G end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_z end_POSTSUPERSCRIPT per unit cell. There is a transition that hole prefers to stay from p1subscript𝑝1p_{1}italic_p start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT to pzsubscript𝑝𝑧p_{z}italic_p start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT at tdp;z=2.50superscriptsubscript𝑡𝑑𝑝𝑧2.50t_{dp;z}^{*}=2.50italic_t start_POSTSUBSCRIPT italic_d italic_p ; italic_z end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT = 2.50eV. Here, we used the parameter introduced in the main-text, tdp;1=1.56subscript𝑡𝑑𝑝11.56t_{dp;1}=1.56italic_t start_POSTSUBSCRIPT italic_d italic_p ; 1 end_POSTSUBSCRIPT = 1.56eV, tdp;2=0.75subscript𝑡𝑑𝑝20.75t_{dp;2}=-0.75italic_t start_POSTSUBSCRIPT italic_d italic_p ; 2 end_POSTSUBSCRIPT = - 0.75eV, tdp;z=1.63subscript𝑡𝑑𝑝𝑧1.63t_{dp;z}=1.63italic_t start_POSTSUBSCRIPT italic_d italic_p ; italic_z end_POSTSUBSCRIPT = 1.63eV, Ua=10subscript𝑈𝑎10U_{a}=10italic_U start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT = 10eV, U=6superscript𝑈6U^{\prime}=6italic_U start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT = 6eV, JH=2subscript𝐽𝐻2J_{H}=2italic_J start_POSTSUBSCRIPT italic_H end_POSTSUBSCRIPT = 2eV, and Δa=9subscriptΔ𝑎9\Delta_{a}=9roman_Δ start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT = 9 eV. (b) tdp;zsubscript𝑡𝑑𝑝𝑧t_{dp;z}italic_t start_POSTSUBSCRIPT italic_d italic_p ; italic_z end_POSTSUBSCRIPT dependence of the binding energy EBsubscript𝐸𝐵E_{B}italic_E start_POSTSUBSCRIPT italic_B end_POSTSUBSCRIPT per unit cell. In the p1subscript𝑝1p_{1}italic_p start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT dominant regime, the binding energy EB,1subscript𝐸𝐵1E_{B,1}italic_E start_POSTSUBSCRIPT italic_B , 1 end_POSTSUBSCRIPT is increasing, implying the rise of Tcsubscript𝑇𝑐T_{c}italic_T start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT under the pressure. Meanwhile in the pzsubscript𝑝𝑧p_{z}italic_p start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT dominant regime, the binding energy EB,zsubscript𝐸𝐵𝑧E_{B,z}italic_E start_POSTSUBSCRIPT italic_B , italic_z end_POSTSUBSCRIPT tends to decrease under the pressure.

VI VI. Discussion on trilayer La4Ni3O10

In this section, we expand the theoretical framework of trilayer nickelates La4Ni3O10 in the charge-transfer regime. According to the DFT data, the valence of the Ni atom of La4Ni3O10 is in the 3d7.333superscript𝑑7.333d^{7.33}3 italic_d start_POSTSUPERSCRIPT 7.33 end_POSTSUPERSCRIPT configuration (Ni+2.67) with the dx2y2subscript𝑑superscript𝑥2superscript𝑦2d_{x^{2}-y^{2}}italic_d start_POSTSUBSCRIPT italic_x start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_y start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_POSTSUBSCRIPT orbital is close to one-third occupied n1=1/3subscript𝑛113n_{1}=1/3italic_n start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT = 1 / 3 per site, while the dz2subscript𝑑superscript𝑧2d_{z^{2}}italic_d start_POSTSUBSCRIPT italic_z start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_POSTSUBSCRIPT orbital is nearly half-filled n21subscript𝑛21n_{2}\approx 1italic_n start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ≈ 1 Sakakibara et al. (2024); Zhang et al. (2024b); Lu et al. (2024b).

Refer to caption
Figure S9: The lattice structure and p,d𝑝𝑑p,ditalic_p , italic_d orbitals of the bilayer Nickelates La4Ni3O10. The green (yellow) sphere denotes Ni (O) atom, respectively. At each Ni atom, there are two d𝑑ditalic_d orbitals, d1,d2subscript𝑑1subscript𝑑2d_{1},d_{2}italic_d start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , italic_d start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT at top/middle/bottom layers. At each O atom at each layer, either pxsubscript𝑝𝑥p_{x}italic_p start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT,pysubscript𝑝𝑦p_{y}italic_p start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT orbitals reside and can hybridize with the Ni atom at each layer. Meanwhile, the O atoms living in the middle of the three layers have pzsubscript𝑝𝑧p_{z}italic_p start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT and should be shared by two Ni atoms at adjacent layers. In our effective charge-transfer model, we focus on the 17 atoms consisting of 3 Ni and 14 O atoms.

VI.1 A. Zhang-Rice spin-half state

We consider the trilayer cubic lattice with nickel dl;1,dl;2subscript𝑑𝑙1subscript𝑑𝑙2d_{l;1},d_{l;2}italic_d start_POSTSUBSCRIPT italic_l ; 1 end_POSTSUBSCRIPT , italic_d start_POSTSUBSCRIPT italic_l ; 2 end_POSTSUBSCRIPT orbitals and oxygen pl;1,pl;2subscript𝑝𝑙1subscript𝑝𝑙2p_{l;1},p_{l;2}italic_p start_POSTSUBSCRIPT italic_l ; 1 end_POSTSUBSCRIPT , italic_p start_POSTSUBSCRIPT italic_l ; 2 end_POSTSUBSCRIPT orbitals defined at each layer l=t,m,b𝑙𝑡𝑚𝑏l=t,m,bitalic_l = italic_t , italic_m , italic_b. The oxygen pl;zsubscript𝑝superscript𝑙𝑧p_{l^{\prime};z}italic_p start_POSTSUBSCRIPT italic_l start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ; italic_z end_POSTSUBSCRIPT orbitals are living between each layer with l=t,bsuperscript𝑙𝑡𝑏l^{\prime}=t,bitalic_l start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT = italic_t , italic_b (See Fig. S9). Now, in each unit cell, we have to consider the total 17 sites of atoms consisting of three Ni atoms and fourteen O atoms. Starting from the local trilayer version of charge-transfer Hamiltonian generalized from Eq.(1) in the main text, we again obtain the Kondo Hamiltonian,

H𝐻\displaystyle Hitalic_H =\displaystyle== l=t,m,ba2JK;a[sl;adsl;ap14nl;ap]+l=t,b2JK;z[(sl;2d+sm;2d)szp12nl;zp]l=t,m,b2JH[sl;1dsl;2d+14],subscript𝑙𝑡𝑚𝑏subscript𝑎2subscript𝐽𝐾𝑎delimited-[]subscriptsuperscript𝑠𝑑𝑙𝑎subscriptsuperscript𝑠𝑝𝑙𝑎14subscriptsuperscript𝑛𝑝𝑙𝑎subscript𝑙𝑡𝑏2subscript𝐽𝐾𝑧delimited-[]subscriptsuperscript𝑠𝑑𝑙2subscriptsuperscript𝑠𝑑𝑚2subscriptsuperscript𝑠𝑝𝑧12subscriptsuperscript𝑛𝑝𝑙𝑧subscript𝑙𝑡𝑚𝑏2subscript𝐽𝐻delimited-[]subscriptsuperscript𝑠𝑑𝑙1subscriptsuperscript𝑠𝑑𝑙214\displaystyle\sum_{l=t,m,b}\sum_{a}2J_{K;a}[\vec{s}^{d}_{l;a}\cdot\vec{s}^{p}_% {l;a}-\frac{1}{4}n^{p}_{l;a}]+\sum_{l=t,b}2J_{K;z}[(\vec{s}^{d}_{l;2}+\vec{s}^% {d}_{m;2})\cdot\vec{s}^{p}_{z}-\frac{1}{2}n^{p}_{l;z}]-\sum_{l=t,m,b}2J_{H}[% \vec{s}^{d}_{l;1}\cdot\vec{s}^{d}_{l;2}+\frac{1}{4}],∑ start_POSTSUBSCRIPT italic_l = italic_t , italic_m , italic_b end_POSTSUBSCRIPT ∑ start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT 2 italic_J start_POSTSUBSCRIPT italic_K ; italic_a end_POSTSUBSCRIPT [ over→ start_ARG italic_s end_ARG start_POSTSUPERSCRIPT italic_d end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_l ; italic_a end_POSTSUBSCRIPT ⋅ over→ start_ARG italic_s end_ARG start_POSTSUPERSCRIPT italic_p end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_l ; italic_a end_POSTSUBSCRIPT - divide start_ARG 1 end_ARG start_ARG 4 end_ARG italic_n start_POSTSUPERSCRIPT italic_p end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_l ; italic_a end_POSTSUBSCRIPT ] + ∑ start_POSTSUBSCRIPT italic_l = italic_t , italic_b end_POSTSUBSCRIPT 2 italic_J start_POSTSUBSCRIPT italic_K ; italic_z end_POSTSUBSCRIPT [ ( over→ start_ARG italic_s end_ARG start_POSTSUPERSCRIPT italic_d end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_l ; 2 end_POSTSUBSCRIPT + over→ start_ARG italic_s end_ARG start_POSTSUPERSCRIPT italic_d end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_m ; 2 end_POSTSUBSCRIPT ) ⋅ over→ start_ARG italic_s end_ARG start_POSTSUPERSCRIPT italic_p end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT - divide start_ARG 1 end_ARG start_ARG 2 end_ARG italic_n start_POSTSUPERSCRIPT italic_p end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_l ; italic_z end_POSTSUBSCRIPT ] - ∑ start_POSTSUBSCRIPT italic_l = italic_t , italic_m , italic_b end_POSTSUBSCRIPT 2 italic_J start_POSTSUBSCRIPT italic_H end_POSTSUBSCRIPT [ over→ start_ARG italic_s end_ARG start_POSTSUPERSCRIPT italic_d end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_l ; 1 end_POSTSUBSCRIPT ⋅ over→ start_ARG italic_s end_ARG start_POSTSUPERSCRIPT italic_d end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_l ; 2 end_POSTSUBSCRIPT + divide start_ARG 1 end_ARG start_ARG 4 end_ARG ] , (S16)

where JK;a,JK;zsubscript𝐽𝐾𝑎subscript𝐽𝐾𝑧J_{K;a},J_{K;z}italic_J start_POSTSUBSCRIPT italic_K ; italic_a end_POSTSUBSCRIPT , italic_J start_POSTSUBSCRIPT italic_K ; italic_z end_POSTSUBSCRIPT is expressed in Eq.(3) in the main-text. The only difference compared to the bilayer model is that pzsubscript𝑝𝑧p_{z}italic_p start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT orbital now has a layer index. There are still three possible one-hole states depending on the p𝑝pitalic_p orbital occupancy : (i)𝑖(i)( italic_i ) nl;1p=1superscriptsubscript𝑛𝑙1𝑝1n_{l;1}^{p}=1italic_n start_POSTSUBSCRIPT italic_l ; 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_p end_POSTSUPERSCRIPT = 1, (ii)𝑖𝑖(ii)( italic_i italic_i ) nl;2p=1superscriptsubscript𝑛𝑙2𝑝1n_{l;2}^{p}=1italic_n start_POSTSUBSCRIPT italic_l ; 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_p end_POSTSUPERSCRIPT = 1, and (iii)𝑖𝑖𝑖(iii)( italic_i italic_i italic_i ) nl;zp=1superscriptsubscript𝑛𝑙𝑧𝑝1n_{l;z}^{p}=1italic_n start_POSTSUBSCRIPT italic_l ; italic_z end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_p end_POSTSUPERSCRIPT = 1. For the first two cases (nl;1p=1superscriptsubscript𝑛𝑙1𝑝1n_{l;1}^{p}=1italic_n start_POSTSUBSCRIPT italic_l ; 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_p end_POSTSUPERSCRIPT = 1 or nl;2p=1superscriptsubscript𝑛𝑙2𝑝1n_{l;2}^{p}=1italic_n start_POSTSUBSCRIPT italic_l ; 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_p end_POSTSUPERSCRIPT = 1) with l=t,m,b𝑙𝑡𝑚𝑏l=t,m,bitalic_l = italic_t , italic_m , italic_b, the local Hamiltonian is

H𝐻\displaystyle Hitalic_H =\displaystyle== 2JK;a[sl;adsl;ap14]l=t,m,b2JH[sl;1dsl;2d+14],2subscript𝐽𝐾𝑎delimited-[]subscriptsuperscript𝑠𝑑𝑙𝑎subscriptsuperscript𝑠𝑝𝑙𝑎14subscriptsuperscript𝑙𝑡𝑚𝑏2subscript𝐽𝐻delimited-[]subscriptsuperscript𝑠𝑑superscript𝑙1subscriptsuperscript𝑠𝑑superscript𝑙214\displaystyle 2J_{K;a}[\vec{s}^{d}_{l;a}\cdot\vec{s}^{p}_{l;a}-\frac{1}{4}]-% \sum_{l^{\prime}=t,m,b}2J_{H}[\vec{s}^{d}_{l^{\prime};1}\cdot\vec{s}^{d}_{l^{% \prime};2}+\frac{1}{4}],2 italic_J start_POSTSUBSCRIPT italic_K ; italic_a end_POSTSUBSCRIPT [ over→ start_ARG italic_s end_ARG start_POSTSUPERSCRIPT italic_d end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_l ; italic_a end_POSTSUBSCRIPT ⋅ over→ start_ARG italic_s end_ARG start_POSTSUPERSCRIPT italic_p end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_l ; italic_a end_POSTSUBSCRIPT - divide start_ARG 1 end_ARG start_ARG 4 end_ARG ] - ∑ start_POSTSUBSCRIPT italic_l start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT = italic_t , italic_m , italic_b end_POSTSUBSCRIPT 2 italic_J start_POSTSUBSCRIPT italic_H end_POSTSUBSCRIPT [ over→ start_ARG italic_s end_ARG start_POSTSUPERSCRIPT italic_d end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_l start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ; 1 end_POSTSUBSCRIPT ⋅ over→ start_ARG italic_s end_ARG start_POSTSUPERSCRIPT italic_d end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_l start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ; 2 end_POSTSUBSCRIPT + divide start_ARG 1 end_ARG start_ARG 4 end_ARG ] ,

whose ground state energy is obtained as,

EG=[2+r+r2+r+1]JH,subscript𝐸𝐺delimited-[]2𝑟superscript𝑟2𝑟1subscript𝐽𝐻\displaystyle E_{G}=-[2+r+\sqrt{r^{2}+r+1}]J_{H},italic_E start_POSTSUBSCRIPT italic_G end_POSTSUBSCRIPT = - [ 2 + italic_r + square-root start_ARG italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + italic_r + 1 end_ARG ] italic_J start_POSTSUBSCRIPT italic_H end_POSTSUBSCRIPT , (S17)

with a dimensionless parameter, r=JK;a/JH𝑟subscript𝐽𝐾𝑎subscript𝐽𝐻r=J_{K;a}/J_{H}italic_r = italic_J start_POSTSUBSCRIPT italic_K ; italic_a end_POSTSUBSCRIPT / italic_J start_POSTSUBSCRIPT italic_H end_POSTSUBSCRIPT. For the last case (nl;zp=1superscriptsubscript𝑛𝑙𝑧𝑝1n_{l;z}^{p}=1italic_n start_POSTSUBSCRIPT italic_l ; italic_z end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_p end_POSTSUPERSCRIPT = 1) with l=t,b𝑙𝑡𝑏l=t,bitalic_l = italic_t , italic_b, the local Hamiltonian is,

H𝐻\displaystyle Hitalic_H =\displaystyle== 2JK;z[(sl;2d+sm;2d)sl;zp12]l=t,m,b2JH[sl;1dsl;2d+14],2subscript𝐽𝐾𝑧delimited-[]subscriptsuperscript𝑠𝑑𝑙2subscriptsuperscript𝑠𝑑𝑚2subscriptsuperscript𝑠𝑝𝑙𝑧12subscriptsuperscript𝑙𝑡𝑚𝑏2subscript𝐽𝐻delimited-[]subscriptsuperscript𝑠𝑑superscript𝑙1subscriptsuperscript𝑠𝑑superscript𝑙214\displaystyle 2J_{K;z}[(\vec{s}^{d}_{l;2}+\vec{s}^{d}_{m;2})\cdot\vec{s}^{p}_{% l;z}-\frac{1}{2}]-\sum_{l^{\prime}=t,m,b}2J_{H}[\vec{s}^{d}_{l^{\prime};1}% \cdot\vec{s}^{d}_{l^{\prime};2}+\frac{1}{4}],2 italic_J start_POSTSUBSCRIPT italic_K ; italic_z end_POSTSUBSCRIPT [ ( over→ start_ARG italic_s end_ARG start_POSTSUPERSCRIPT italic_d end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_l ; 2 end_POSTSUBSCRIPT + over→ start_ARG italic_s end_ARG start_POSTSUPERSCRIPT italic_d end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_m ; 2 end_POSTSUBSCRIPT ) ⋅ over→ start_ARG italic_s end_ARG start_POSTSUPERSCRIPT italic_p end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_l ; italic_z end_POSTSUBSCRIPT - divide start_ARG 1 end_ARG start_ARG 2 end_ARG ] - ∑ start_POSTSUBSCRIPT italic_l start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT = italic_t , italic_m , italic_b end_POSTSUBSCRIPT 2 italic_J start_POSTSUBSCRIPT italic_H end_POSTSUBSCRIPT [ over→ start_ARG italic_s end_ARG start_POSTSUPERSCRIPT italic_d end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_l start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ; 1 end_POSTSUBSCRIPT ⋅ over→ start_ARG italic_s end_ARG start_POSTSUPERSCRIPT italic_d end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_l start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ; 2 end_POSTSUBSCRIPT + divide start_ARG 1 end_ARG start_ARG 4 end_ARG ] ,

whose ground state energy is obtained as,

EGz=12[43r4+8r+9r2]JH,superscriptsubscript𝐸𝐺𝑧12delimited-[]43𝑟48𝑟9superscript𝑟2subscript𝐽𝐻\displaystyle E_{G}^{z}=\frac{1}{2}[-4-3r-\sqrt{4+8r+9r^{2}}]J_{H},italic_E start_POSTSUBSCRIPT italic_G end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_z end_POSTSUPERSCRIPT = divide start_ARG 1 end_ARG start_ARG 2 end_ARG [ - 4 - 3 italic_r - square-root start_ARG 4 + 8 italic_r + 9 italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG ] italic_J start_POSTSUBSCRIPT italic_H end_POSTSUBSCRIPT , (S18)

with a dimensionless parameter, r=JK;z/JH𝑟subscript𝐽𝐾𝑧subscript𝐽𝐻r=J_{K;z}/J_{H}italic_r = italic_J start_POSTSUBSCRIPT italic_K ; italic_z end_POSTSUBSCRIPT / italic_J start_POSTSUBSCRIPT italic_H end_POSTSUBSCRIPT. The above results are just a constant shift by JHsubscript𝐽𝐻-J_{H}- italic_J start_POSTSUBSCRIPT italic_H end_POSTSUBSCRIPT from the Eqs.(S6-S7), which leads to the same conclusion. Based on the parameters listed in Ref.Luo et al. (2024), we assume tdp;1tdp;z>tdp;2similar-to-or-equalssubscript𝑡𝑑𝑝1subscript𝑡𝑑𝑝𝑧subscript𝑡𝑑𝑝2t_{dp;1}\simeq t_{dp;z}>t_{dp;2}italic_t start_POSTSUBSCRIPT italic_d italic_p ; 1 end_POSTSUBSCRIPT ≃ italic_t start_POSTSUBSCRIPT italic_d italic_p ; italic_z end_POSTSUBSCRIPT > italic_t start_POSTSUBSCRIPT italic_d italic_p ; 2 end_POSTSUBSCRIPT and Δp;1Δp;2zΔp;2similar-to-or-equalssubscriptΔ𝑝1superscriptsubscriptΔ𝑝2𝑧similar-to-or-equalssubscriptΔ𝑝2\Delta_{p;1}\simeq\Delta_{p;2}^{z}\simeq\Delta_{p;2}roman_Δ start_POSTSUBSCRIPT italic_p ; 1 end_POSTSUBSCRIPT ≃ roman_Δ start_POSTSUBSCRIPT italic_p ; 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_z end_POSTSUPERSCRIPT ≃ roman_Δ start_POSTSUBSCRIPT italic_p ; 2 end_POSTSUBSCRIPT, leading to the relation JK;14JK;z>JK;2subscript𝐽𝐾14subscript𝐽𝐾𝑧subscript𝐽𝐾2J_{K;1}\approx 4J_{K;z}>J_{K;2}italic_J start_POSTSUBSCRIPT italic_K ; 1 end_POSTSUBSCRIPT ≈ 4 italic_J start_POSTSUBSCRIPT italic_K ; italic_z end_POSTSUBSCRIPT > italic_J start_POSTSUBSCRIPT italic_K ; 2 end_POSTSUBSCRIPT. Using the fact JK;14JK;z>JK;2subscript𝐽𝐾14subscript𝐽𝐾𝑧subscript𝐽𝐾2J_{K;1}\approx 4J_{K;z}>J_{K;2}italic_J start_POSTSUBSCRIPT italic_K ; 1 end_POSTSUBSCRIPT ≈ 4 italic_J start_POSTSUBSCRIPT italic_K ; italic_z end_POSTSUBSCRIPT > italic_J start_POSTSUBSCRIPT italic_K ; 2 end_POSTSUBSCRIPT, we again have EG1(JK;1)<EG2(JK;2)superscriptsubscript𝐸𝐺1subscript𝐽𝐾1superscriptsubscript𝐸𝐺2subscript𝐽𝐾2E_{G}^{1}(J_{K;1})<E_{G}^{2}(J_{K;2})italic_E start_POSTSUBSCRIPT italic_G end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 1 end_POSTSUPERSCRIPT ( italic_J start_POSTSUBSCRIPT italic_K ; 1 end_POSTSUBSCRIPT ) < italic_E start_POSTSUBSCRIPT italic_G end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( italic_J start_POSTSUBSCRIPT italic_K ; 2 end_POSTSUBSCRIPT ) and EG1(JK;1)<EGz(JK;z4JK;1)superscriptsubscript𝐸𝐺1subscript𝐽𝐾1superscriptsubscript𝐸𝐺𝑧similar-to-or-equalssubscript𝐽𝐾𝑧4subscript𝐽𝐾1E_{G}^{1}(J_{K;1})<E_{G}^{z}(J_{K;z}\simeq 4J_{K;1})italic_E start_POSTSUBSCRIPT italic_G end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 1 end_POSTSUPERSCRIPT ( italic_J start_POSTSUBSCRIPT italic_K ; 1 end_POSTSUBSCRIPT ) < italic_E start_POSTSUBSCRIPT italic_G end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_z end_POSTSUPERSCRIPT ( italic_J start_POSTSUBSCRIPT italic_K ; italic_z end_POSTSUBSCRIPT ≃ 4 italic_J start_POSTSUBSCRIPT italic_K ; 1 end_POSTSUBSCRIPT ). Hence, we can conclude that the Zhang-Rice spin-1/2 state, especially with p1subscript𝑝1p_{1}italic_p start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT hole is the most stable state even for the trilayer nickelates.

VI.2 B. Trilayer layer Type-II t-J model

Then by taking the Zhang-Rice spin-half as the primary state of the d8Lsuperscript𝑑8𝐿d^{8}Litalic_d start_POSTSUPERSCRIPT 8 end_POSTSUPERSCRIPT italic_L state and keeping a spin-triplet doublon state of the d8superscript𝑑8d^{8}italic_d start_POSTSUPERSCRIPT 8 end_POSTSUPERSCRIPT state, the minimal model of the hole-doped trilayer bilayer nickelates is the type-II t-J model Zhang and Vishwanath (2020); Oh and Zhang (2023). The trilayer type-II t-J model Hamiltonian is given by

H=HK𝐻subscript𝐻𝐾\displaystyle H=H_{K}italic_H = italic_H start_POSTSUBSCRIPT italic_K end_POSTSUBSCRIPT +\displaystyle++ l=t,m,bi,j[Jsssi;lsj;l+Jsd(si;lSj;l+Si;lsj;l)+JddSi;lSj;l]\displaystyle\sum_{l=t,m,b}\sum_{\langle i,j\rangle}\left[J_{ss}^{\parallel}% \vec{s}_{i;l}\cdot\vec{s}_{j;l}+J_{sd}^{\parallel}(\vec{s}_{i;l}\cdot\vec{S}_{% j;l}+\cdot\vec{S}_{i;l}\cdot\vec{s}_{j;l})+J_{dd}^{\parallel}\vec{S}_{i;l}% \cdot\vec{S}_{j;l}\right]∑ start_POSTSUBSCRIPT italic_l = italic_t , italic_m , italic_b end_POSTSUBSCRIPT ∑ start_POSTSUBSCRIPT ⟨ italic_i , italic_j ⟩ end_POSTSUBSCRIPT [ italic_J start_POSTSUBSCRIPT italic_s italic_s end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ∥ end_POSTSUPERSCRIPT over→ start_ARG italic_s end_ARG start_POSTSUBSCRIPT italic_i ; italic_l end_POSTSUBSCRIPT ⋅ over→ start_ARG italic_s end_ARG start_POSTSUBSCRIPT italic_j ; italic_l end_POSTSUBSCRIPT + italic_J start_POSTSUBSCRIPT italic_s italic_d end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ∥ end_POSTSUPERSCRIPT ( over→ start_ARG italic_s end_ARG start_POSTSUBSCRIPT italic_i ; italic_l end_POSTSUBSCRIPT ⋅ over→ start_ARG italic_S end_ARG start_POSTSUBSCRIPT italic_j ; italic_l end_POSTSUBSCRIPT + ⋅ over→ start_ARG italic_S end_ARG start_POSTSUBSCRIPT italic_i ; italic_l end_POSTSUBSCRIPT ⋅ over→ start_ARG italic_s end_ARG start_POSTSUBSCRIPT italic_j ; italic_l end_POSTSUBSCRIPT ) + italic_J start_POSTSUBSCRIPT italic_d italic_d end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ∥ end_POSTSUPERSCRIPT over→ start_ARG italic_S end_ARG start_POSTSUBSCRIPT italic_i ; italic_l end_POSTSUBSCRIPT ⋅ over→ start_ARG italic_S end_ARG start_POSTSUBSCRIPT italic_j ; italic_l end_POSTSUBSCRIPT ]
+\displaystyle++ l=t,bi[Jsssi;lsi;m+Jsd(si;lSi;m+Si;lsi;m)+JddSi;lSi;m]+Vl=t,bini;lni;m,subscript𝑙𝑡𝑏subscript𝑖delimited-[]superscriptsubscript𝐽𝑠𝑠perpendicular-tosubscript𝑠𝑖𝑙subscript𝑠𝑖𝑚superscriptsubscript𝐽𝑠𝑑perpendicular-tosubscript𝑠𝑖𝑙subscript𝑆𝑖𝑚subscript𝑆𝑖𝑙subscript𝑠𝑖𝑚superscriptsubscript𝐽𝑑𝑑perpendicular-tosubscript𝑆𝑖𝑙subscript𝑆𝑖𝑚𝑉subscript𝑙𝑡𝑏subscript𝑖subscript𝑛𝑖𝑙subscript𝑛𝑖𝑚\displaystyle\sum_{l=t,b}\sum_{i}\left[J_{ss}^{\perp}\vec{s}_{i;l}\cdot\vec{s}% _{i;m}+J_{sd}^{\perp}(\vec{s}_{i;l}\cdot\vec{S}_{i;m}+\vec{S}_{i;l}\cdot\vec{s% }_{i;m})+J_{dd}^{\perp}\vec{S}_{i;l}\cdot\vec{S}_{i;m}\right]+V\sum_{l=t,b}% \sum_{i}n_{i;l}n_{i;m},∑ start_POSTSUBSCRIPT italic_l = italic_t , italic_b end_POSTSUBSCRIPT ∑ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT [ italic_J start_POSTSUBSCRIPT italic_s italic_s end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ⟂ end_POSTSUPERSCRIPT over→ start_ARG italic_s end_ARG start_POSTSUBSCRIPT italic_i ; italic_l end_POSTSUBSCRIPT ⋅ over→ start_ARG italic_s end_ARG start_POSTSUBSCRIPT italic_i ; italic_m end_POSTSUBSCRIPT + italic_J start_POSTSUBSCRIPT italic_s italic_d end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ⟂ end_POSTSUPERSCRIPT ( over→ start_ARG italic_s end_ARG start_POSTSUBSCRIPT italic_i ; italic_l end_POSTSUBSCRIPT ⋅ over→ start_ARG italic_S end_ARG start_POSTSUBSCRIPT italic_i ; italic_m end_POSTSUBSCRIPT + over→ start_ARG italic_S end_ARG start_POSTSUBSCRIPT italic_i ; italic_l end_POSTSUBSCRIPT ⋅ over→ start_ARG italic_s end_ARG start_POSTSUBSCRIPT italic_i ; italic_m end_POSTSUBSCRIPT ) + italic_J start_POSTSUBSCRIPT italic_d italic_d end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ⟂ end_POSTSUPERSCRIPT over→ start_ARG italic_S end_ARG start_POSTSUBSCRIPT italic_i ; italic_l end_POSTSUBSCRIPT ⋅ over→ start_ARG italic_S end_ARG start_POSTSUBSCRIPT italic_i ; italic_m end_POSTSUBSCRIPT ] + italic_V ∑ start_POSTSUBSCRIPT italic_l = italic_t , italic_b end_POSTSUBSCRIPT ∑ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT italic_n start_POSTSUBSCRIPT italic_i ; italic_l end_POSTSUBSCRIPT italic_n start_POSTSUBSCRIPT italic_i ; italic_m end_POSTSUBSCRIPT ,

with

HKsubscript𝐻𝐾\displaystyle H_{K}italic_H start_POSTSUBSCRIPT italic_K end_POSTSUBSCRIPT =\displaystyle== tl=t,m,bσ,i,jci;l;σcj;l;σtl=t,bσ,ici;l;σci;m;σ+H.c.,formulae-sequencesubscript𝑡parallel-tosubscript𝑙𝑡𝑚𝑏subscript𝜎𝑖𝑗superscriptsubscript𝑐𝑖𝑙𝜎subscript𝑐𝑗𝑙𝜎subscript𝑡perpendicular-tosubscript𝑙𝑡𝑏subscript𝜎𝑖superscriptsubscript𝑐𝑖𝑙𝜎subscript𝑐𝑖𝑚𝜎Hc\displaystyle-t_{\parallel}\sum_{l=t,m,b}\sum_{\sigma,\langle i,j\rangle}c_{i;% l;\sigma}^{\dagger}c_{j;l;\sigma}-t_{\perp}\sum_{l=t,b}\sum_{\sigma,i}c_{i;l;% \sigma}^{\dagger}c_{i;m;\sigma}+\mathrm{H.c.},- italic_t start_POSTSUBSCRIPT ∥ end_POSTSUBSCRIPT ∑ start_POSTSUBSCRIPT italic_l = italic_t , italic_m , italic_b end_POSTSUBSCRIPT ∑ start_POSTSUBSCRIPT italic_σ , ⟨ italic_i , italic_j ⟩ end_POSTSUBSCRIPT italic_c start_POSTSUBSCRIPT italic_i ; italic_l ; italic_σ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT italic_c start_POSTSUBSCRIPT italic_j ; italic_l ; italic_σ end_POSTSUBSCRIPT - italic_t start_POSTSUBSCRIPT ⟂ end_POSTSUBSCRIPT ∑ start_POSTSUBSCRIPT italic_l = italic_t , italic_b end_POSTSUBSCRIPT ∑ start_POSTSUBSCRIPT italic_σ , italic_i end_POSTSUBSCRIPT italic_c start_POSTSUBSCRIPT italic_i ; italic_l ; italic_σ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT italic_c start_POSTSUBSCRIPT italic_i ; italic_m ; italic_σ end_POSTSUBSCRIPT + roman_H . roman_c . ,

where the coefficients J,Jsubscript𝐽parallel-tosubscript𝐽perpendicular-toJ_{\parallel},J_{\perp}italic_J start_POSTSUBSCRIPT ∥ end_POSTSUBSCRIPT , italic_J start_POSTSUBSCRIPT ⟂ end_POSTSUBSCRIPT are listed in Eqs.(S9-S14).