---
title: Three-Orbital Emery Model in Cuprates
url: https://www.emergentmind.com/topics/three-orbital-emery-model
type: topic
---

# Three-Orbital Emery Model in Cuprates

The three-orbital Emery model, also called the three-band Hubbard model, is the canonical multiorbital lattice model for a CuO\(_2\) plane in which one explicitly retains the Cu \(d_{x^2-y^2}\) orbital and the two planar oxygen \(p_\sigma\) orbitals \(p_x\) and \(p_y\). Its defining purpose is to encode Cu–O hybridization, O–O hopping, charge-transfer physics, and strong local correlations on copper within a single low-energy framework. Relative to a one-band Hubbard reduction, it keeps the oxygen degrees of freedom explicit and therefore treats orbital occupancies, Zhang–Rice singlets, particle–hole asymmetry, pseudogap phenomenology, and transport asymmetry in a more microscopic way [2403.11218][2409.20497]. Subsequent work has also clarified that the conventional three-hopping version is often qualitatively effective but can be quantitatively insufficient for realistic cuprate superconductivity and transport unless longer-range processes are restored [2605.07739][2604.08085].

## 1. Orbital content and conceptual role

The model is defined on the CuO\(_2\) geometry, equivalently a Lieb-lattice arrangement, with Cu orbitals on the vertices of a square lattice and oxygen orbitals on the horizontal and vertical links. The active local basis is
\[
(d_{x^2-y^2},\,p_x,\,p_y),
\]
or, in momentum space,
\[
\psi^{\dagger}_{\mathbf{k},\sigma}=(d^{\dagger}_{\mathbf{k},\sigma},p^{\dagger}_{x,\mathbf{k},\sigma},p^{\dagger}_{y,\mathbf{k},\sigma}).
\]
This orbital structure is the minimal one that directly distinguishes copper and oxygen charge sectors and thereby permits a model-internal discussion of where doped carriers reside [2403.11218][2605.07739].

A central reason the Emery model matters is that cuprates are charge-transfer materials rather than simple single-band Mott systems. In the Zaanen–Sawatzky–Allen language, the relation between the local Cu repulsion and the Cu–O level separation distinguishes Mott-Hubbard-like and charge-transfer-insulator-like regimes. One determinant quantum Monte Carlo study states this explicitly as \(U_{dd}<\Delta\) for a Mott-Hubbard-like regime and \(U_{dd}>\Delta\) for a charge-transfer-insulator-like regime, with \(\Delta=\epsilon_p-\epsilon_d\) in its hole-language convention [2509.11028]. A cellular dynamical mean-field formulation uses an electron-language convention with
\[
\Delta=\epsilon_d+U_d-\tilde{\epsilon}_p,
\]
showing that the precise definition depends on representation and on how the oxygen level is renormalized [2607.06462]. This suggests that comparisons across the literature require care not only about parameters, but also about conventions for particles, holes, and double counting.

The model is also the natural setting for discussing the asymmetry between hole-doped and electron-doped cuprates. Several studies report that holes predominantly enter oxygen-derived states while electrons predominantly enter copper-derived states, an asymmetry that is unavailable in a one-orbital description [2407.19545][2503.07810]. At the same time, recent comparisons between the one-band Hubbard and three-orbital Emery models show that the two can share broad trends while differing strongly in spectral gaps, resistivity, effective masses, and the coupling strength required to match experiment [2604.08085]. The model is therefore used both as a direct theory of CuO\(_2\) planes and as the reference from which one judges the adequacy of one-band reductions.

## 2. Hamiltonian, conventions, and parameterizations

A widely used hole-language form is
\[
H = K_0 + K_{pd} + K_{pp} + V_{dd} + V_{pp},
\]
with
\[
K_0 = (\epsilon_d - \mu)\sum_{i\sigma} n^d_{i\sigma} + (\epsilon_p - \mu)\sum_{j\alpha\sigma} n^p_{j\alpha\sigma},
\]
\[
K_{pd} = \sum_{\langle ij \rangle \alpha \sigma} t_{pd}^{i,j,\alpha} \left(d^\dagger_{i\sigma} p_{j\alpha\sigma} + p^\dagger_{j\alpha\sigma} d_{i\sigma}\right),
\]
\[
K_{pp} = \sum_{\langle jj' \rangle \alpha\alpha' \sigma} t_{pp}^{j,j',\alpha,\alpha'} \left(p^\dagger_{j\alpha\sigma} p_{j'\alpha'\sigma} + p^\dagger_{j'\alpha'\sigma} p_{j\alpha\sigma}\right),
\]
\[
V_{dd} = U_{dd}\sum_i n^d_{i\uparrow}n^d_{i\downarrow}, \qquad V_{pp} = U_{pp}\sum_{j,\alpha} n^p_{j\alpha\uparrow}n^p_{j\alpha\downarrow}.
\]
Here \(\epsilon_d,\epsilon_p\) are Cu and O on-site energies, \(\mu\) is the chemical potential, \(t_{pd}\) and \(t_{pp}\) are Cu–O and O–O hoppings, and \(U_{dd},U_{pp}\) are local Coulomb repulsions on Cu and O orbitals [2409.20497][2403.11218].

Many modern formulations instead use a \(3\times 3\) one-particle matrix \(\mathbf h_0(\mathbf k)\) with an explicit next-nearest-neighbor oxygen–oxygen hopping \(t'_{pp}\), together with a local interaction only on Cu:
\[
H=\sum_{\mathbf{k}\sigma} C_{\mathbf{k}\sigma}^\dagger \big[\mathbf{h}_0(\mathbf{k})-\mu \mathbf{I}\big] C_{\mathbf{k}\sigma} +U_d\sum_{\mathbf{R}_i} n_{d\mathbf{R}_i\uparrow}n_{d\mathbf{R}_i\downarrow}.
\]
In this family of parameterizations,
\[
V_{dp_\alpha}=t_{pd}(1-e^{-ik_\alpha}), \qquad W_{p_\alpha p_\alpha}=2t_{pp}'\cos k_\alpha, \qquad W_{p_xp_y}=t_{pp}(1-e^{ik_x})(1-e^{-ik_y}),
\]
and the oxygen level often appears as the renormalized \(\tilde{\epsilon}_p=\epsilon_p-2t_{pp}\) [2607.06462][2311.09023]. Other works use equivalent matrix forms with different gauge choices, and one DQMC study emphasizes explicitly that the hopping-sign convention is gauge-dependent and not unique [2509.11028].

Filling conventions vary. In a hole representation, one work sets \(\rho=1\) at half-filling and defines
\[
\delta \equiv |\rho-1|.
\]
In an electron-language convention used in several CDMFT studies, the half-filled three-orbital unit cell has
\[
n_d+2n_p=5,
\]
and doping is written
\[
\delta=5-n_{\rm tot}.
\]
Thus \(\delta>0\) may denote hole doping in one convention and a magnitude of deviation from \(\rho=1\) in another [2409.20497][2407.19545][2503.07810]. For technical work this is not a semantic issue; it changes the meaning of “half-filling,” the sign of the dopant, and the definition of the charge-transfer scale.

Representative parameter sets depend on the physical regime being targeted. For cuprate-relevant DCA studies, one set is
\[
U_{dd}=8.5~{\rm eV},\quad t_{pd}=1.13~{\rm eV},\quad t_{pp}=0.49~{\rm eV},\quad \epsilon_d=0,\quad \epsilon_p=3.24~{\rm eV},
\]
with comparisons between \(U_{pp}=0\) and \(U_{pp}=4.1~{\rm eV}\) [2409.20497]. Another DCA study keeps
\[
U_{dd}=7.5~{\rm eV},\quad U_{pp}=0,\quad t_{pd}=1.13~{\rm eV},\quad t_{pp}=0.49~{\rm eV},\quad \epsilon_d=0,
\]
and varies \(\epsilon_p\) between \(3.24\) and \(6.0\) to interpolate between cuprate-like and nickelate-like behavior [2403.11218]. A DQMC study instead fixes
\[
t_{pd}=1.0,\quad t_{pp}=0.4,\quad \epsilon_d=0,
\]
and varies \(U_{dd}=4.0\text{ to }8.0\) and \(\epsilon_p=2.0\text{ to }6.0\), emphasizing systematic scans through interaction and charge-transfer regimes [2509.11028].

The “conventional Emery model” is often taken to include only three hopping parameters, \(t_{pd}\), \(t_{pp}\), and \(t'_{pp}\). A 2026 study argues that this conventional three-parameter truncation is not quantitatively sufficient and replaces it by the full \(d\)-\(p\) Wannier Hamiltonian, retaining all significant longer-range in-plane hoppings \(\bar h_r(\mathbf R-\mathbf R')\) obtained from a DFT-based downfolding [2605.07739]. This has become a major point of contemporary discussion because it bears directly on realistic phase diagrams and on whether the model should be regarded as merely three-orbital or also explicitly longer ranged.

## 3. Solution strategies and diagnostic observables

The three-orbital Emery model has been treated with an unusually broad range of many-body techniques. Finite-temperature normal-state calculations include the Dynamical Cluster Approximation with a continuous-time auxiliary-field quantum Monte Carlo solver, determinant quantum Monte Carlo, cellular dynamical mean-field theory with CT-HYB, single-site and cluster DMFT, DMFT with NRG for transport, the dynamical vertex approximation in its \(\lambda\)-corrected ladder formulation, and a multiorbital TPSC+DMFT scheme [2403.11218][2509.11028][2607.06462][2604.08085][2308.14091]. Ordered and quasi-one-dimensional problems have additionally been studied by variational Monte Carlo, exact diagonalization, density-matrix renormalization group, and channel-decomposed functional renormalization group [1305.7275][2603.10755][1310.3136].

Cluster approaches are prominent because they access momentum differentiation. In DCA calculations, particular emphasis is placed on the antinodal momentum \(\mathbf K=(\pi,0)\) and on the nodal point \((\pi/2,\pi/2)\). One study mainly uses \(N_c=16\) for momentum-resolved scattering rates and \(N_c=4\) for lower-temperature access when the sign problem becomes severe, while another focuses on a small \(N_c=4\) DCA cluster specifically to extract the doping dependence of the \(T\)-linear scattering slope [2403.11218][2409.20497]. CDMFT work on superconductivity and ambipolar doping uses a 12-site cluster containing \(N_d=4\) Cu sites and \(N_p=8\) O sites, which is the minimal cluster capable of supporting \(d\)-wave superconductivity in that framework [2607.06462][2407.19545].

Several observables recur across the literature. The electronic scattering rate is commonly defined from the zero-frequency self-energy as
\[
\gamma_{\mathbf K}\equiv -\operatorname{Im}\Sigma(\mathbf K,\omega=0),
\]
or, in one implementation,
\[
\gamma_k \equiv -\mathrm{Im}\,\Sigma^{(2)}(\mathbf{K},\omega=0).
\]
Because \(\omega=0\) is not directly available on the Matsubara axis, the low-frequency Matsubara self-energy is fitted and extrapolated to zero frequency, typically with a second-order polynomial [2403.11218][2409.20497]. The quasiparticle scattering rate is then defined by
\[
\frac{1}{\tau_{\mathbf K}} = Z_{\mathbf K}\gamma_{\mathbf K},
\qquad
Z_{\mathbf K}\approx \left[1-\frac{\operatorname{Im}\Sigma(\mathbf K,i\omega_0)}{\omega_0}\right]^{-1}.
\]
For transport-motivated analysis, one work further uses the empirical relation
\[
\alpha_1(\delta)\propto \delta\,A_1(\delta)
\]
to estimate a linear-in-\(T\) resistivity coefficient \(A_1\) from computed scattering slopes [2409.20497].

Other diagnostics are tailored to the charge-transfer problem. In CDMFT, the charge-transfer gap \(\Delta_{\rm I}\) can be extracted from the compressibility plateau at \(n_{\rm tot}=5\), with the plateau operationally defined by
\[
\frac{dn_{\rm tot}}{d\mu}<0.005.
\]
The oxygen hole content is measured as
\[
2p_p = 2(2-n_p).
\]
These quantities are used as state variables for organizing superconducting trends in the Zaanen–Sawatzky–Allen plane [2607.06462]. In DMFT studies of magnetic response, the uniform static spin susceptibility is computed as
\[
\chi_m(\mathbf q=0,\Omega=0)=\left.\frac{\partial m}{\partial H}\right|_{H=0},
\]
with \(m=n_{d\uparrow}-n_{d\downarrow}\), enabling direct comparison to Knight-shift measurements [2311.09023].

Methodological work has also clarified how multiorbital self-consistency should be organized. In the TPSC+DMFT formulation for the Emery model, interacting orbital densities rather than noninteracting ones must be used in the sum rules, and the effective spin vertex \(U_{sp}\) decreases rapidly with filling at fixed bare \(U_d\) [2308.14091]. This does not merely modify a technical detail; it alters the interpretation of effective correlation strength across the phase diagram.

## 4. Linear-\(T\) scattering, transport, and non-Fermi-liquid regimes

A major contemporary use of the Emery model is the analysis of \(T\)-linear scattering and strange-metal transport. In DCA+CT-AUX calculations with cuprate-like \(\epsilon_p=3.24\), the electronic scattering rate exhibits a nearly linear-in-\(T\) behavior at low temperature for intermediate densities, roughly \(\rho\sim 1.2-1.5\), while low densities show strong increases associated with proximity to the charge-transfer insulating state and pseudogap-like momentum differentiation [2403.11218]. At larger \(\epsilon_p=6.0\), intended as a nickelate-like regime, the scattering is generally larger, more isotropic, and develops a downturn below \(T\sim 0.1\), with evidence from \(N_c=4\) calculations for two consecutive nearly linear-\(T\) regimes separated by a smooth crossover around \(T\sim 0.08\text{--}0.1\) [2403.11218].

The doping dependence of the linear-in-\(T\) slope \(\alpha_1\) has been studied directly in the three-orbital model. For the electronic scattering rate, the slope grows approximately linearly with electron doping, whereas on the hole-doped side it is approximately inversely proportional to doping at intermediate doping, \(\alpha_1\sim 1/\delta\), before crossing over around \(\delta\sim 0.3\) to a more linear-like dependence at larger hole doping [2409.20497]. This electron–hole asymmetry survives the inclusion of finite oxygen interaction: \(U_{pp}\) has little effect on the electron-doped side, enhances the magnitude and slope on the hole-doped side, but does not qualitatively change the trend [2409.20497].

The same work shows that the quasiparticle scattering slope is less clean than the electronic one. On the hole-doped side it retains an inverse-like trend, but on the electron-doped side it is more sensitive to the approximation used for \(Z_{\mathbf K}\). The estimated resistivity coefficient \(A_1\), inferred through \(\alpha_1(\delta)\propto \delta A_1(\delta)\), comes out roughly inverse proportional to doping on both sides; agreement with experiment is better for hole doping than for electron doping, leading the authors to treat the electronic scattering rate as the more reliable comparator to transport data within that calculation [2409.20497].

These results also constrain “Planckian” interpretations. In the temperature-dependent DCA study, the quasiparticle scattering rate generically departs from the unity slope predicted by a strict Planckian dissipation picture; approximate unity occurs only in narrow cases near \(\rho\sim 1.4\) for \(\epsilon_p=3.24\) and \(\rho\sim 1.5\) for \(\epsilon_p=6.0\) [2403.11218]. The Emery model therefore supports extended non-Fermi-liquid-like linear-\(T\) regimes, but not a universal or parameter-insensitive Planckian slope.

A broader transport comparison between the single-band Hubbard and three-orbital Emery models reaches a related conclusion from a different direction. Within DMFT+NRG, the two models yield a similar broad picture but not the same quantitative transport. The Emery model generally gives a larger resistivity than the Hubbard model at comparable gap size, and no single mapping of coupling strengths works across all observables: matching gap size suggests \(U^{\mathrm{Hubbard}}\approx \tfrac23 U^{\mathrm{Emery}}\), matching low-energy spectra suggests \(U^{\mathrm{Hubbard}}\approx U^{\mathrm{Emery}}\), and matching dc resistivity suggests \(U^{\mathrm{Hubbard}}\approx \tfrac54 U^{\mathrm{Emery}}\) [2604.08085]. This is consistent with the view that transport is especially sensitive to the explicit oxygen sector rather than just to a renormalized single-band coupling.

## 5. Charge-transfer physics, pseudogap behavior, and orbital differentiation

The most robust normal-state result across methods is the orbital asymmetry of doping. In CDMFT, ambipolar doping of a charge-transfer insulator shows that doped electrons mostly enter copper orbitals while doped holes mostly enter oxygen orbitals [2407.19545]. DQMC studies recover the same asymmetry and further show that it strengthens as \(\epsilon_p\) and \(U_{dd}\) increase on the hole-doped side [2509.11028]. A zero-temperature CDMFT study of the \(n_d\)-\(n_p\) relation similarly finds that larger \(U\) sharpens the asymmetry across half-filling, with holes becoming more oxygen-like and electrons more copper-like [2503.07810]. This pattern is often taken as a hallmark of the charge-transfer-insulator regime.

Yet the nature of the doping-driven transition is not simply “electron side metallic, hole side pseudogapped.” In a 12-site CDMFT study restricted to the normal state, both electron and hole doping produce a two-stage evolution:
\[
\text{charge-transfer insulator} \rightarrow \text{pseudogap} \rightarrow \text{correlated metal}.
\]
The CTI-to-pseudogap change is continuous at \(\delta=0\), while the pseudogap-to-correlated-metal transition is first order at low temperature and terminates at a finite-\(T\) critical end point, with a Widom line above it [2407.19545]. The concrete implication drawn there is that merely doping the same charge-transfer insulator with opposite carrier signs is not sufficient to explain the full electron–hole asymmetry of the cuprate normal-state phase diagram.

Pseudogap structure in the Emery model is also strongly controlled by the charge-transfer scale. At smaller \(\epsilon_p\), momentum-resolved DQMC spectra show stronger nodal–antinodal differentiation, while at \(\epsilon_p=6\) the pseudogap-like suppression is largely absent [2509.11028]. The same study interprets the larger-\(\epsilon_p\) regime as more relevant to infinite-layer nickelates and the smaller-\(\epsilon_p\) regime as more cuprate-like. This matches the DCA finding that larger \(\epsilon_p\) yields more isotropic scattering and a low-\(T\) downturn rather than a straightforward extension of the cuprate-like linear-\(T\) regime [2403.11218].

The model is also widely used to discuss Zhang–Rice physics. DQMC local and momentum-resolved spectra identify a ZRS-like low-energy feature near the Fermi level at low and moderate hole doping, but find that in the heavily overdoped regime the low-energy local density of states can split into multiple peaks, with redistributed oxygen weight, which is interpreted as a possible Zhang–Rice singlet breakdown [2509.11028]. Finite \(U_{pp}\), although often neglected to reduce the sign problem, tends to stabilize the ZRS against hole doping and to push the onset of low-energy splitting to higher doping [2509.11028]. A separate DMFT study likewise identifies the low-energy quasiparticle peak as ZRS-like, but emphasizes that single-site DMFT cannot produce the momentum-selective spectral pseudogap even though it does reproduce the non-Curie-like Knight-shift downturn through emerging oxygen–copper singlet fluctuations [2311.09023]. This is a common point of clarification: in the Emery model, one may capture singlet-induced magnetic anomalies without yet reproducing the full pseudogap spectrum.

Charge distribution observables have also been tied directly to experiment. By matching CDMFT \(n_d\)-\(n_p\) curves to NMR estimates of Cu and O occupations, approximate values \(U-\epsilon_p\approx 10\) for LCO, \(6\) for YBCO, and \(5\) for NCCO were inferred within that parameterization [2503.07810]. The same study reports that the Emery model works reasonably well for hole-doped LCO and YBCO but less satisfactorily for electron-doped NCCO, suggesting either underestimated correlation strength or limitations of that CDMFT treatment for the electron-doped material [2503.07810].

## 6. Magnetism, loop currents, and superconductivity

The ordered-state literature on the three-orbital Emery model is structurally diverse. In a variational Monte Carlo and exact-diagonalization study of an Emery-type Hamiltonian with Cu–O and O–O hoppings, onsite \(U_d,U_p\), charge-transfer energy \(\Delta\), and nearest-neighbor \(V_{dp}\), a loop-current state of the \(\theta_2\) type is stabilized for parameter choices compatible with a typical hole-doped Fermi surface, provided the oxygen–oxygen transfer signs are modified by indirect Cu-4s–mediated processes [1305.7275]. For \(t_{pp}=-0.3\) eV and \(\epsilon_p-\epsilon_d\approx 1\) eV, the resulting sequence is AFM for \(\delta<5\%\), loop current for roughly \(5\%\lesssim\delta\lesssim 12\%\), \(d\)-wave superconductivity for \(12\%<\delta<20\%\), and Fermi liquid at higher doping [1305.7275]. That work interprets the loop-current state as lowering interaction energy at a kinetic-energy cost and argues that the relevant sign structure of O–O transfers is central to its stability.

A weak-coupling channel-decomposed fRG study instead finds the dominant competition to be between antiferromagnetism and \(d_{x^2-y^2}\)-wave superconductivity. Around generic hole-doped parameters, the flow exhibits a strong AFM–\(d\)-wave competition with a smooth doping dependence, and increasing incommensurability in the magnetic channel deforms the \(d\)-wave gap near its maxima [1310.3136]. That analysis concludes that multiorbital effects do not radically alter the qualitative instability structure relative to a one-band model, but they do matter quantitatively through orbital makeup and longer-range effective hopping.

In cluster and diagrammatic strong-coupling approaches, the central superconducting question has become what controls the dome and the maximal pairing scale. A 2026 DΓA study argues that the conventional three-hopping Emery model is qualitatively correct but quantitatively insufficient: for CaCuO\(_2\), the full \(d\)-\(p\) model with long-range hopping gives a superconducting dome from about \(7\%\) to \(22\%\) hole doping, whereas the conventional model yields only \(9\%\) to \(14\%\) [2605.07739]. In the conventional model, the van Hove singularity is too close to the Fermi level, pseudogap effects turn on too early, antinodal spectral weight is suppressed, and the superconducting eigenvector becomes strongly modulated instead of exhibiting a more conventional cuprate \(d\)-wave form [2605.07739]. The argument is not that long-range hoppings pair carriers directly, but that they correct the Fermi surface and effective one-band \(t'\) into the range where \(d\)-wave superconductivity is optimized.

A complementary CDMFT study formulates the superconducting problem in terms of parent-state observables. It finds that the maximum superconducting critical temperature \(T_c^{\rm max}\) increases monotonically as the oxygen hole content \(2p_p\) increases and the charge-transfer gap size \(\Delta_{\rm I}\) decreases, with the oxygen hole content being the dominant variable in the explored regime [2607.06462]. In that work, optimal superconductivity occurs both near the metal–insulator boundary and deep in the charge-transfer regime, showing that these are distinct mechanisms rather than equivalent descriptions. This gives the Emery model an internal materials-design language: parent-state covalency and parent-state gap size are separate levers.

Quasi-one-dimensional calculations preserve the same basic theme. DMRG studies on ladder supercells that keep the physical Cu:O ratio \(1:2\) find that the undoped ladders are charge-transfer insulators and that doping turns them into Luther–Emery liquids with enhanced pairing correlations [2603.10755]. Because these ladders preserve the Cu/O stoichiometry, they allow a direct study of how pairing strength varies with the Cu–O charge distribution; the pairing strength is reported to peak around
\[
0.3 \lesssim 2n_p \lesssim 0.4,
\]
resembling the experimentally observed nonmonotonic relation between superconductivity and oxygen-hole content [2603.10755]. This reinforces the broader CDMFT conclusion that oxygen participation is not a secondary detail but a controlling variable for pairing.

## 7. Extensions, reductions, and non-solid-state realizations

A persistent question is whether the Emery model should be viewed as a starting point to be downfolded or as the minimal faithful model in its own right. Weak-coupling fRG comparisons find qualitative agreement between the Emery model and an effective one-band Hubbard model but report that the one-band description underestimates the instability scale and misses quantitative multiorbital corrections from orbital makeup and longer-range hopping [1310.3136]. DMFT+NRG comparisons reach a similar conclusion from the strong-coupling side: both models reproduce broad trends, but there is no observable-independent mapping of \(U\), and transport in particular often requires a larger effective coupling in the one-band model than expected [2604.08085]. A plausible implication is that downfolding can remain phenomenologically useful while still being non-universal across observables.

The model has also been extended beyond cuprates. In nickelate-oriented work, raising \(\epsilon_p\) toward \(6.0\) suppresses pseudogap behavior and yields more isotropic scattering, already suggesting a deformation of the cuprate problem within the same three-orbital framework [2403.11218][2509.11028]. A further extension adds an interstitial \(s\) orbital to form a four-orbital \(d\)-\(p\)-\(s\) model for infinite-layer nickelates. DQMC then finds that strong correlations substantially reduce the \(s\)-derived electron pocket, push the \(\Gamma\)-centered pocket above the Fermi level while leaving an \(A\)-centered pocket of experimentally comparable size, flatten the \(d_{x^2-y^2}\) \(k_z\) dispersion, and enhance short-range antiferromagnetic correlations relative to the conventional three-orbital \(d\)-\(p\) model [2603.20705]. This indicates that the three-orbital Emery model is not merely a cuprate-specific construct; it is also the baseline from which more realistic nickelate models are built.

A distinct line of work proposes direct quantum simulation of the Emery model with ultracold atoms. One proposal uses a bichromatic optical superlattice to realize the three-band geometry, retain tunable \(U_d\), \(U_p\), and \(\Delta_{pd}\), benchmark the tight-binding model through single-particle quantum walks, and infer effective single-band Hubbard parameters by Hamiltonian learning; for the studied finite system the learned values are reported around
\[
U/t\approx 16.7,\qquad t'/t\approx -0.1
\]
[2604.22955]. Another proposal implements the model in a phase-stable optical lattice with a controllable interference term, achieving accessible parameters around \(\Delta_{pd}/t_{pd}\approx 2\) to \(8\), \(t_{pp}/t_{pd}\approx 0.1\), and \(t_{pp}'/t_{pd}\approx -0.2\), explicitly targeting both cuprate-like and nickelate-like regimes [2603.11037]. These developments suggest that the Emery model has become not only a theoretical framework for correlated oxides but also a candidate Hamiltonian for controlled analog simulation.

In its present usage, the “three-orbital Emery model” therefore denotes more than a single fixed Hamiltonian. It is a family of closely related Cu–O multiband models whose common core is the explicit retention of \(d_{x^2-y^2}\), \(p_x\), and \(p_y\) orbitals. Around that core, the current literature debates how much realism must be restored—oxygen interactions, \(t'_{pp}\), full long-range \(d\)-\(p\) hoppings, three-dimensional spectator bands, or experimentally calibrated charge-transfer scales—to obtain the correct pseudogap, transport, and superconducting phenomenology. The broad consensus across these studies is that the explicit oxygen sector is indispensable whenever the questions concern charge transfer, orbital occupancies, electron–hole asymmetry, or the quantitative structure of strange-metal and superconducting behavior [2409.20497][2605.07739].

Source: https://www.emergentmind.com/topics/three-orbital-emery-model