Interplay of two $E_{g}$ orbitals in Superconducting La${}_{3}$Ni${}_{2}$O${}_{7}$ Under Pressure

The different orientations of two $E_{g}$ orbitals result in the strong inter-layer hopping $t^{z}_{\perp}$ between two $3d_{z^{2}}$ orbitals along the rung, and also the relatively strong intra-layer hopping $t^{x}_{\parallel}$ between the $3d_{x^{2}-y^{2}}$ orbitals. The $3d_{z^{2}}$ and $3d_{x^{2}-y^{2}}$ orbitals within the intra-layer nearest neighbor (NN) bond further hybridize, but those along the rung do not. The corresponding band structure is formulated in the following tight-binding Hamiltonian,

	$\displaystyle H_{0}$	$\displaystyle=-t^{x}_{\parallel}\sum_{\langle i,j\rangle\alpha\sigma}\big{(}d_{x^{2}\alpha\sigma}^{\dagger}(i)d_{x^{2}\alpha\sigma}(j)+\text{h.c.}\big{)}+\Delta_{g}\sum_{i\alpha\sigma}d_{x^{2}\alpha\sigma}^{\dagger}(i)d_{x^{2}\alpha\sigma}(i)$		(1)
		$\displaystyle-t^{z}_{\perp}\sum_{i\sigma}\big{(}d_{z^{2}1\sigma}^{\dagger}(i)d_{z^{2}2\sigma}(i)+\text{h.c.}\big{)}-t^{z}_{\parallel}\sum_{\langle i,j\rangle\alpha\sigma}\big{(}d_{z^{2}\alpha\sigma}^{\dagger}(i)d_{z^{2}\alpha\sigma}(j)+\text{h.c.}\big{)}-\sum_{\langle i,j\rangle\alpha\sigma}t_{\parallel,j-i}^{xz}\big{(}d_{x^{2}\alpha\sigma}^{\dagger}(i)d_{z^{2}\alpha\sigma}(j)+\text{h.c.}\big{)}$

where $i$ is the lattice site index; $\alpha$ represents the layer index; $\sigma$ is the spin component index. $d_{a\alpha\sigma}^{\dagger}(i)$ represents the electron creation operator for the $3d_{x^{2}-y^{2}}$ orbital ($a=x^{2}$) or the $3d_{z^{2}}$ orbital ($a=z^{2}$) with $a$ referring to the orbital index. The $3d_{z^{2}}$ orbital has a smaller on-site energy as described by the $\Delta_{g}$ term. The intra-layer hybridization $t_{\parallel,j-i}^{xz}$ between $3d_{x^{2}-y^{2}}$ and $3d_{z^{2}}$ orbitals on neighboring sites is comparable to hopping strengthes $t^{x}_{\parallel}$ and $t^{z}_{\perp}$. The NN hybridization takes opposite signs along the $x$ and $y$ direction, respectively, i.e., $t^{xz}_{\parallel,x}=-t^{xz}_{\parallel,y}\equiv t^{xz}_{\parallel}$, due to the different symmetries of orbital wavefunctions. A finite but small intra-layer hopping $t^{z}_{\parallel}$ for the $3d_{z^{2}}$ orbital is also added into Eq. (1).

The multi-orbital Hubbard interaction reads,

		$\displaystyle H_{1}=U\sum_{ia}n_{a\alpha\uparrow}(i)n_{a\alpha\downarrow}(i)+V\sum_{i\alpha}n_{z^{2}\alpha}(i)n_{x^{2}\alpha}(i)$		(2)
		$\displaystyle+P\sum_{i\alpha}\left(d^{\dagger}_{x^{2}\alpha\uparrow}(i)d^{\dagger}_{x^{2}\alpha\downarrow}(i)d_{z^{2}\alpha\downarrow}(i)d_{z^{2}\alpha\uparrow}(i)+\text{h.c}\right)$
		$\displaystyle-J_{H}\sum_{i\alpha}\bm{S}_{z^{2}\alpha}(i)\cdot\bm{S}_{x^{2}\alpha}(i),$

where the $U$- and $V$-terms describe the intra-orbital and inter-orbital repulsions, respectively; the $P$-term describes the singlet pair hopping between two orbitals; the $J_{H}$-term is the Hund’s rule coupling aligning electron spins in two onsite $E_{g}$ orbitals; $n_{s\alpha\sigma}(i)=d_{s\alpha\sigma}^{\dagger}(i)d_{s\alpha\sigma}(i)$ represents the particle number operator for an electron with spin-$\sigma$ in orbital $s$ with $s=x^{2}$ oder $z^{2}$; $n_{s\alpha}(i)=n_{s\alpha\uparrow}(i)+n_{s\alpha\downarrow}(i)$; $\bm{S}_{s\alpha}(i)=\frac{1}{2}d_{s\alpha}^{\dagger}(i)[\bm{\sigma}]d_{s\alpha}(i)$ is the spin operator for the electron in orbital $s$, where Pauli matrices $\bm{\sigma}=(\sigma_{x},\sigma_{y},\sigma_{z})$.

In the strong coupling limit, the Hubbard-$U$ interaction generates an effective inter-layer AFM spin superexchange $J_{\perp}$ for the $3d_{z^{2}}$ orbitals and an intra-layer one $J_{\parallel}$ for the $3d_{x^{2}-y^{2}}$ orbitals, respectively. Two electrons in the same site but different $E_{g}$ orbitals are aligned in their spin direction due to the strong Hund’s coupling. Consequently, the inter-layer superexchange $J_{\perp}$ of $3d_{z^{2}}$ orbital is transmitted to the $3d_{x^{2}-y^{2}}$ orbital [25]. Please note that this interlayer AFM coupling does not arise from the virtual hopping due to the inter-layer tunneling between two $3d_{x^{2}-y^{2}}$ orbitals, which is negligible.

We arrive at an effective bilayer $t$-$J$ two-orbital model as depicted in Fig. 1: The band part, or, the $t$-term, looks the same as Eq. (1) but under the constraint of no double-occupancy in each orbital. The interaction part reads

	$\displaystyle H_{J,x^{2}}$	$\displaystyle=J_{\parallel}\sum_{\langle i,j\rangle\alpha}\bm{S}_{x^{2}\alpha}(i)\cdot\bm{S}_{x^{2}\alpha}(j)$		(3)
		$\displaystyle+J_{\perp}\sum_{i}\bm{S}_{x^{2}1}(i)\cdot\bm{S}_{x^{2}2}(i),$
	$\displaystyle H_{J,z^{2}}$	$\displaystyle=J_{\perp}\sum_{i}\bm{S}_{z^{2}1}(i)\cdot\bm{S}_{z^{2}2}(i).$

The $3d_{x^{2}-y^{2}}$-orbital $H_{J,x^{2}}$ takes the form of a bilayer $t$-$J_{\parallel}$-$J_{\perp}$ model [72, 25]. Similar single-orbital bilayer $t$-$J_{\parallel}$-$J_{\perp}$ models have also been considered in the context of cuprates  [63, 64, 65, 66, 67, 68, 69, 70, 71, 72] and in cold atom systems [95, 96, 97]. For the $3d_{z^{2}}$-orbital $H_{z^{2}}$, only the inter-layer superexchange $J_{\perp}$ is relevant.

We adopt the hopping parameters from the DFT calculation provided by Ref.  [32] for the LNO material. The relevant parameters are taken as, $t^{x}_{\parallel}=0.526$ eV, $t^{z}_{\parallel}=0.113$ eV, $t^{z}_{\perp}=0.676$ eV, $t^{xz}_{\parallel}=0.25$ eV, $\Delta_{g}=0.528$ eV. The superexchange coupling parameters $J_{\parallel}$ and $J_{\perp}$ should in principle be determined by various interaction parameters in the strong-coupling limit. For simplicity, we use the simple relation $J\approx 4t^{2}/U$ by neglecting other interaction parameters. Since the accurate value of the Hubbard $U$ is hard to obtain, $J_{\parallel}$ and $J_{\perp}$ are set as tuning parameters. Since $J_{\perp}/J_{\parallel}\approx(t^{z}_{\perp}/t^{x}_{\parallel})^{2}$, from the tight-binding parameters provided in Ref. [32] and Ref. [12], the ratio between two superexchange strengths is approximated as $J_{\perp}/J_{\parallel}\approx 1.7-1.8$. Similarly, $J_{\perp}/J_{\parallel}\approx 2$ as given in Ref. [21]. In our calculation, their values are tuned with their ratio fixed as $J_{\perp}/J_{\parallel}=2$ for convenience. An estimation in Ref. [84] provides $U=4.75$ eV, leading to $J_{\parallel}\approx 0.23$ eV. On the other hand, a smaller value of $J_{\parallel}=0.1$eV is used Ref. [21]. Considering the different values of parameters in literature, $J_{\parallel}$ is varied from 0.1 to 0.3eV in our calculation to cover these variations.

In the slave boson approach  [91, 94], the electron operators of two $E_{g}$ orbitals are represented as $d_{a\alpha\sigma}^{\dagger}(i)=f_{a\alpha\sigma}^{\dagger}(i)b_{a\alpha}(i)$, where $f_{a\alpha\sigma}^{\dagger}(i)$ is the spinon creation operator and $b_{a\alpha}(i)$ is the holon annihilation operator for the two orbitals ($a=x^{2},z^{2}$) respectively. The following no-double-occupancy constraint is imposed on each orbital at each site,

$\displaystyle\sum_{\sigma}f_{a\alpha\sigma}^{\dagger}(i)f_{a\alpha\sigma}(i)+b_{a\alpha}^{\dagger}(i)b_{a\alpha}(i)=1.$

(4)

These constraints are imposed through Lagrange multipliers $\lambda_{a\alpha}(i)$ and in the mean-field (MF) theory, $\lambda_{a\alpha}(i)\equiv\lambda_{a}$ on average irrespective of the layer and site.

The total electron filling $x_{\text{tot}}$ of the $E_{g}$ orbital is determined from the material, which yeids $x_{\text{tot}}=0.75$ for LNO. The NN hybridization causes the density fluctuation between two orbitals while keeps the total holon density $\delta_{\text{tot}}=2-2x_{\text{tot}}$ in the $E_{g}$ orbitals per Ni site fixed,

$\displaystyle\frac{1}{N}\sum_{i\alpha}\big{(}b_{x^{2}\alpha}^{\dagger}(i)b_{x^{2}\alpha}(i)+b_{z^{2}\alpha}^{\dagger}(i)b_{z^{2}\alpha}(i)\big{)}=2\delta_{\text{tot}}.$

(5)

For the case of undoped LNO, $\delta_{\text{tot}}=0.5$. Upon chemical doping, $\delta_{\text{tot}}$ varies within $0\sim 1$, and is controlled by $\lambda_{\text{tot}}$. $\lambda_{\text{tot}}$ plays the role of chemical potential. In addition to $\lambda_{\text{tot}}$, $\lambda_{x^{2}}$, $\lambda_{z^{2}}$ are also introduced. Their combinations together with the onsite $\Delta_{g}$ act as 4 chemical potentials for the four degrees of freedom, i.e., spinons and holons in two $E_{g}$ orbitals, respectively.

After the decomposition, the neareast neighboring AFM superexchange $J_{\parallel}$ and $J_{\perp}$ terms give rise to spinon’s intra-orbital bondings and singlet pairings. The $\chi_{a}^{(\alpha)}(i,j)$ ($\alpha=1,2$) and $\chi_{a}^{\perp}(i)$ represent spinon’s intra- and inter-layer hopping, respectively. $\Delta_{a}^{(\alpha)}(i,j)$ and $\Delta_{a}^{\perp}(i)$ represent the intra- and inter-layer singlet pairings, respectively. The electron hopping terms of $t^{a}_{\parallel}$ and $t^{z}_{\perp}$ are decoupled as the intra-orbital spinon hopping $\chi_{a}^{(\alpha)}(i,j)$ and $\chi_{z}^{\perp}(i)$ as well as holon hopping $\kappa_{a}^{(\alpha)}(i,j)$ and $\kappa_{z}^{\perp}(i)$, respectively. The inter-orbital electron hybridization $t_{\parallel,j-i}^{xz}$ term is decoupled as the inter-orbital spinon hybridization $\chi_{xz}^{(\alpha)}(i,j)$ and holon hybridization $\kappa_{xz}^{(\alpha)}(i,j)$ ($\alpha=1,2$), respectively.

The various parameters defined above are expressed in terms of the expectation values of the corresponding fermion bi-linear operators at thermal equilibrium. At the mean-field level, order parameters are assumed to be site/layer independent. The spinon hopping parameters are expressed as

	$\displaystyle\chi_{\parallel}^{a}$	$\displaystyle=\frac{3}{8}J_{\parallel}\langle f_{a\alpha\uparrow}^{\dagger}(j)f_{a\alpha\uparrow}(i)+f_{a\alpha\downarrow}^{\dagger}(j)f_{a\alpha\downarrow}(i)\rangle,$		(6)
	$\displaystyle\chi^{a}_{\perp}$	$\displaystyle=\frac{3}{8}J_{\perp}\langle f_{a2\uparrow}^{\dagger}(i)f_{a1\uparrow}(i)+f_{a2\downarrow}^{\dagger}(i)f_{a1\downarrow}(i)\rangle,$
	$\displaystyle\chi_{\parallel}^{xz}$	$\displaystyle=t_{\parallel}^{xz}\langle f_{z^{2}\alpha\uparrow}^{\dagger}(j)f_{x^{2}\alpha\uparrow}(i)+f_{z^{2}j\alpha\downarrow}^{\dagger}(j)f_{x^{2}\alpha\downarrow}(i)\rangle,$

respectively, where $\langle\cdots\rangle$ means the thermal average. The spinon pairing order parameters are expressed as

	$\displaystyle\Delta_{\parallel}^{a}$	$\displaystyle=\frac{3}{8}J_{\parallel}\langle f_{a\alpha\downarrow}(j)f_{a\alpha\uparrow}(i)-f_{a\alpha\uparrow}(j)f_{a\alpha\downarrow}(i)\rangle,$		(7)
	$\displaystyle\Delta^{a}_{\perp}$	$\displaystyle=\frac{3}{8}J_{\perp}\langle f_{a2\downarrow}(i)f_{a1\uparrow}(i)-f_{a2\uparrow}(i)f_{a1\downarrow}(i)\rangle,$

respectively. Similarly for the holon part, their hopping parameters are defined as follows as

	$\displaystyle\kappa_{\parallel}^{a}=$	$\displaystyle t_{\parallel}^{a}\langle b_{a\alpha}^{\dagger}(j)b_{a\alpha}(i)\rangle,$		(8)
	$\displaystyle\kappa^{a}_{\perp}=$	$\displaystyle t_{\perp}^{a}\langle b_{a2}^{\dagger}(i)b_{a1}(i)\rangle,$
	$\displaystyle\kappa_{\parallel}^{xz}=$	$\displaystyle t_{\parallel}^{xz}\langle b_{x^{2}\alpha}^{\dagger}(j)b_{z^{2}\alpha}(i)\rangle.$