Realizing the Emery Model in Optical Lattices for Quantum Simulation of Cuprates and Nickelates

This paper proposes a quantum simulation scheme using ultracold atoms in optical lattices to realize the three-band Emery model, enabling the study of high-temperature superconductivity in cuprates and nickelates on system sizes that are currently inaccessible to numerical methods.

The Gist

Imagine you are trying to understand why a specific type of metal (called a cuprate) conducts electricity with zero resistance at surprisingly high temperatures. This is the "holy grail" of physics: High-Temperature Superconductivity.

For decades, scientists have been looking at this problem through a simplified lens, like trying to understand a complex orchestra by only listening to the lead violinist. This simplified model is called the Fermi-Hubbard model. It's been very useful, but recent evidence suggests it's missing crucial details. The real "music" of these materials involves a whole section of instruments playing together, not just one.

This paper proposes a way to build a quantum simulator to hear the whole orchestra.

The Problem: The "Three-Player" Game

The real materials (cuprates and a newer family called nickelates) are made of layers of Copper and Oxygen atoms.

• The Old View: Scientists used to pretend the Oxygen atoms were just background noise and focused only on the Copper atoms.
• The New View: The Oxygen atoms are actually active players! They interact with the Copper in a specific way. This requires a more complex model called the Emery Model, which tracks three types of "players" (orbitals) at once: Copper and two types of Oxygen.

The problem? This three-player game is so mathematically messy that even the world's fastest supercomputers struggle to solve it for large groups of atoms.

The Solution: Building a "Lego" Universe with Light

The authors propose using ultracold atoms (atoms cooled to near absolute zero) trapped in a grid of laser beams (an optical lattice). Think of this as building a miniature, controllable universe out of light.

Here is how they plan to build it, using a creative analogy:

1. The Stage (The Lattice)

Imagine a checkerboard made of laser light. Usually, every square on the board is identical. But to simulate the Emery model, we need a special pattern.

• The Setup: They use a single laser beam that bounces back and forth.
• The Trick: By inserting a special mirror (a half-wave plate) in the middle of the path, they can twist the "polarization" (the direction the light waves vibrate) of the beam on the return trip.
• The Result: This creates a subtle difference in energy between the "Copper" spots and the "Oxygen" spots on the board. It's like painting every other square of the checkerboard a slightly different shade of blue, creating an energy gap that mimics the real material.

2. The Actors (The Atoms)

They drop fermionic atoms (like Lithium-6) into this laser grid. These atoms act as the "electrons" (or more accurately, "holes" in the material).

• The Magic: By adjusting the angle of that special mirror, they can dial in the exact energy difference between the Copper and Oxygen spots. This allows them to simulate both Cuprates (where the energy gap is small) and Nickelates (where the gap is huge) just by turning a knob.

3. The Goal: Finding the "Secret Sauce"

Why go through all this trouble? Because they want to see if a specific phenomenon called the Zhang-Rice Singlet forms.

• The Analogy: Imagine a Copper atom is a lonely dancer. When an extra "hole" (a missing electron) enters the system, it doesn't just sit on the Copper. It grabs the Copper's hand and forms a tight, inseparable dance pair with the surrounding Oxygen atoms.
• Why it matters: This "dance pair" (the singlet) is believed to be the secret ingredient that allows electricity to flow without resistance.
• The Prediction: The authors ran simulations on a computer to predict what would happen in their laser setup. They found that:
  • In Cuprates, these dance pairs form easily, and the magnetic order breaks down in a specific way.
  • In Nickelates, the energy gap is so large that the dance pairs are much harder to form, leading to different physics.

Why This Matters

Currently, we can't build a supercomputer big enough to simulate a whole sheet of this material to see if it becomes a superconductor. But we can build this laser setup.

• The "Real" Material: Instead of guessing with math, they can literally watch the atoms dance in the laser grid.
• The Outcome: If they see the "dance pairs" (Zhang-Rice singlets) forming and behaving exactly as their simulations predict, it confirms that the three-band Emery model is the correct way to understand high-temperature superconductivity.

In a Nutshell

This paper is a blueprint for building a quantum playground using lasers and cold atoms. It's designed to let scientists play with the complex rules of copper-oxide materials in a controlled environment, finally allowing them to solve the 40-year-old mystery of how these materials become superconductors. It's like moving from trying to solve a puzzle in the dark to turning on the lights and seeing the pieces fit together.

Technical Summary

Here is a detailed technical summary of the paper "Realizing the Emery Model in Optical Lattices for Quantum Simulation of Cuprates and Nickelates."

1. Problem Statement

The microscopic origin of high-temperature superconductivity in cuprates remains a central unsolved problem in condensed matter physics. While the single-band Fermi-Hubbard model has been extensively studied, growing evidence suggests it is insufficient to fully capture the physics of cuprates and infinite-layer nickelates. Specifically, the single-band model may miss crucial multi-band effects, such as charge-transfer physics and the formation of Zhang-Rice singlets.

The canonical three-band Emery model, which describes copper-oxide planes with copper ($d$) and oxygen ($p$) orbitals on a Lieb lattice geometry, offers a more accurate description. However, simulating this model numerically is extremely challenging due to the large system sizes required to capture low-temperature phenomena and the complexity of the three-band Hilbert space. Current numerical methods (DMRG, Monte Carlo, DMFT) struggle with the necessary system sizes and parameter regimes.

2. Methodology

The authors propose a quantum simulation scheme using ultracold fermionic atoms in optical lattices to realize the Emery model.

• Lattice Geometry: The proposal utilizes a Lieb lattice geometry, previously realized experimentally, consisting of a central site ($d$-orbital, representing Copper) and two surrounding sites ($p_x, p_y$ orbitals, representing Oxygen) per unit cell.
• Tunable Charge-Transfer Energy ($\Delta_{pd}$): A key innovation is the method to introduce the energy offset $\Delta_{pd}$ (the charge-transfer gap) between $d$ and $p$ sites.
  • The setup uses a single-beam, passively phase-stable optical lattice in a bow-tie geometry.
  • By inserting a Half-Wave Plate (HWP) between the forward and retro-reflected passes of the laser, the polarization of the beam along the $y$-axis is rotated by an angle $\theta$ relative to the $x$-axis.
  • This polarization difference creates an interference contrast that results in a tunable potential offset between nearest-neighbor sites, effectively controlling $\Delta_{pd}$.
• Site Selectivity: To achieve the specific Lieb lattice connectivity (where $p$ sites connect only to $d$ sites, not each other directly in the tight-binding limit, though $p$-$p$ hopping is induced), the authors propose superimposing a pattern of repulsive potentials using a Digital Micromirror Device (DMD). This selectively lifts the energy of every fourth site (the "copper" sites in the raw lattice) to create the desired unit cell structure.
• Particle-Hole Transformation: The simulation operates in the hole picture. By transforming the energy axis ($E_p \to E_h$), the system maps to a scenario with half-filled $d$-bands and empty $p$-bands. This allows for a low filling fraction ($\langle n \rangle = 1/3$), which is experimentally and numerically more tractable than the particle picture ($\langle n \rangle = 5/3$).
• Hamiltonian: The realized Hamiltonian is:
  $$\hat{H} = -t_{pd} \sum_{\langle ij \rangle, \sigma} (\hat{d}^\dagger_{i\sigma} \hat{p}_{j\sigma} + h.c.) - \sum_{jj', \sigma} t_{pp}^{(j,j')} (\hat{p}^\dagger_{j\sigma} \hat{p}_{j'\sigma} + h.c.) + \Delta_{pd} \sum_{j, \sigma} \hat{p}^\dagger_{j\sigma} \hat{p}_{j\sigma} + \sum_{\nu \in \{p,d\}} U_\nu \sum_i \hat{n}_{i\uparrow}^\nu \hat{n}_{i\downarrow}^\nu$$
  where $t_{pd}$ is the $d$-$p$ hopping, $t_{pp}$ and $t'_{pp}$ are $p$-$p$ hoppings, and $U_\nu$ are on-site interactions.

3. Key Contributions

• Experimental Protocol: The paper provides a concrete, experimentally feasible protocol to realize the Emery model using current ultracold atom technology, specifically leveraging polarization control to tune $\Delta_{pd}$ without complex multi-beam phase-locking.
• Parameter Accessibility: The authors demonstrate that the proposed setup can access parameter regimes relevant to both cuprates ($\Delta_{pd}/t_{pd} \approx 3.5$) and infinite-layer nickelates ($\Delta_{pd}/t_{pd} \approx 9$).
• Bridging Theory and Experiment: By combining the experimental proposal with numerical benchmarks, the work establishes a pathway to study multi-band physics on system sizes (e.g., cylinders of $12 \times 2$ unit cells) that are currently inaccessible to classical numerical methods.

4. Results

The authors performed numerical simulations (using Density Matrix Renormalization Group - DMRG) on cylinders with $12 \times 2$ unit cells to predict observables for the proposed experiment.

• Band Structure Fitting: The optical lattice band structure was fitted to the Emery model, confirming that the tunable polarization angle $\theta$ allows continuous control of $\Delta_{pd}/t_{pd}$ from $\approx 2$ to $\approx 8$, covering the regimes of both material classes.
• Spin Correlations:
  • Antiferromagnetism: Nearest-neighbor spin correlations between $d$-$d$ sites decrease with doping, similar to the single-band Hubbard model.
  • Asymmetry: A distinct particle-hole asymmetry was observed. Hole-doping leads to strong correlations between $d$ and $p$ sites, while electron-doping does not. This asymmetry is a hallmark of the multi-band Emery model not captured by single-band models.
• Zhang-Rice Singlets:
  • In the cuprate regime (hole-doped), the number of occupied $d$-$p$-$d$ bonds scales linearly with doping ($N_{bond} \approx N_h/2$), indicating the formation of Zhang-Rice singlets (a bound state of a hole on a $d$-site and a symmetric superposition of holes on surrounding $p$-sites).
  • In the nickelate regime (large $\Delta_{pd}$), Zhang-Rice singlet formation is suppressed, and additional dopants predominantly reside on $d$-sites, leading to different magnetic behavior.
• Magnetic Order: The simulations show that while antiferromagnetic order exists in both regimes, it decays faster in nickelates due to the competition between the large charge-transfer energy and other energy scales.

5. Significance

• Validation of Multi-Band Physics: This work provides a route to experimentally verify that the single-band Hubbard model is an approximation that may fail to capture essential physics (like Zhang-Rice singlets) in real materials.
• Material Simulation: It enables the direct simulation of infinite-layer nickelates, a recently discovered class of superconductors, allowing researchers to understand why their physics differs from cuprates (specifically regarding the role of $\Delta_{pd}$).
• Scalability: By using quantum simulators, the field can move beyond the limitations of numerical methods to study larger systems and lower temperatures, potentially uncovering the mechanism of unconventional superconductivity.
• Effective Models: The platform allows for the systematic derivation of effective single-band models from multi-band physics by measuring single-particle coherences and constructing many-body Wannier orbitals, bridging the gap between microscopic models and effective theories.

In conclusion, this paper presents a robust experimental blueprint for realizing the Emery model, offering a powerful tool to resolve long-standing questions about the nature of high-temperature superconductivity in copper-oxide and nickelate materials.