Pairing and charge distribution in Emery ladders preserving the ratio of Cu to O atoms

Using the density-matrix renormalization group method, this study demonstrates that Emery ladders preserving the Cu-to-O ratio transition from charge-transfer insulators to Luther-Emery liquids with enhanced pairing upon doping, successfully reproducing the relationships between charge distribution, pairing strength, and interactions observed in both two-dimensional models and experiments.

The Gist

Imagine you are trying to understand how a complex machine works, like a high-performance race car engine. The engine is a giant, flat sheet of metal (a 2D surface) with thousands of tiny, interacting parts. Studying the whole flat sheet at once is incredibly difficult because there are too many variables, and our current computers can't handle the math.

So, physicists often take a shortcut: they study a ladder.

Think of a ladder as a narrow strip cut out of that giant engine sheet. If the strip is narrow enough, we can simulate it perfectly on a computer. The hope is that if we understand how the ladder behaves, we can predict how the whole engine works. This is the core idea of this paper.

The Problem with Old Ladders

For decades, scientists have been building these "ladders" to study cuprates—a family of materials that can conduct electricity with zero resistance (superconductivity) at surprisingly high temperatures.

However, the authors of this paper, Gökmen Polat and Eric Jeckelmann, realized that the ladders everyone else was using were broken models.

Here is the analogy: Imagine the real engine (the cuprate) is made of a specific recipe: 1 part Copper and 2 parts Oxygen. It's like a cake that needs exactly 1 cup of flour and 2 cups of sugar.

The old ladders scientists were using were like cakes made with the wrong recipe:

• Some ladders had a ratio of 2 parts Copper to 3 parts Oxygen.
• Others had 2 parts Copper to 5 parts Oxygen.

It's like trying to bake a cake with 2 cups of flour and 3 cups of sugar, then claiming you've learned how to bake the original cake. The flavor (the physics) might be okay, but the chemical balance (the charge distribution) is wrong. Because the ratio was off, scientists couldn't accurately predict how the "flour" (electrons on Copper) and "sugar" (electrons on Oxygen) would share the load when the cake was baked (doped with extra ingredients).

The New Solution: Perfectly Balanced Ladders

In this paper, the authors built three new types of ladders. They carefully designed these strips so that they perfectly preserve the 1:2 ratio of Copper to Oxygen, just like the real 2D material.

Think of these new ladders as miniature, perfect replicas of the real engine, just sliced into a long, narrow strip. Even though these new strips look a bit more "jagged" or less symmetrical than the old smooth ones, they have the correct ingredients.

What They Found

Using a powerful computer technique called DMRG (which is like a super-smart way to solve puzzles by breaking them into smaller, manageable pieces), they tested these new ladders.

Here are the key discoveries, translated into everyday language:

1. The "Off" Switch (Undoped State): When the material has no extra electricity flowing through it (undoped), it acts like a perfect insulator. The electrons are stuck in place, like cars stuck in a traffic jam. The authors confirmed their new ladders behave exactly like this, acting as "charge-transfer insulators."
2. The "On" Switch (Doped State): When they added a few extra "cars" (holes/electrons) to the traffic, the jam broke. The material transformed into a Luther-Emery liquid.
   • Analogy: Imagine the traffic jam suddenly turning into a synchronized dance. The cars (electrons) start pairing up and moving together in perfect harmony. This "pairing" is the secret sauce of superconductivity.
3. The Perfect Ratio Matters: Because their new ladders had the correct 1:2 Copper-to-Oxygen ratio, the authors could finally see exactly where the extra electricity was going.
   • When they added "holes" (missing electrons), the holes mostly went to the Oxygen atoms.
   • When they added "electrons," they mostly removed them from the Copper atoms.
   • This specific behavior matches what experiments see in real life. The old, broken-ratio ladders couldn't show this clearly.

Why This Is a Big Deal

The authors found that these new, "imperfectly shaped" but "perfectly balanced" ladders show the exact same superconducting behavior as the old ladders, but with a crucial advantage: they tell the truth about the ingredients.

They discovered that the strength of the superconducting "dance" (pairing) is directly linked to how the electricity is shared between the Copper and Oxygen atoms.

• If you tweak the "recipe" (change the energy levels of the atoms), the dance gets stronger or weaker.
• They found a "sweet spot" where the dance is strongest, which matches what scientists see in real-world experiments.

The Bottom Line

This paper is like a mechanic realizing that to fix a race car, they need to stop testing it on a model with the wrong fuel mixture. By building a model with the exact right ratio of ingredients, they can finally understand how the fuel (electrons) distributes itself to create the super-speed (superconductivity).

This new approach gives scientists a much clearer window into the mysterious world of high-temperature superconductors, potentially helping us design better materials for lossless power grids, faster computers, and maglev trains in the future.

Technical Summary

Here is a detailed technical summary of the paper "Pairing and charge distribution in Emery ladders preserving the ratio of Cu to O atoms" by Gökmen Polat and Eric Jeckelmann.

1. Problem Statement

Understanding the mechanism of high-temperature superconductivity in cuprates remains a major theoretical challenge. The Emery model (three-band Hubbard model), which explicitly includes copper ($d$) and oxygen ($p$) orbitals, is the standard framework for describing the low-energy physics of $\text{CuO}_2$ planes.

While studying 2D systems directly is computationally prohibitive for strongly correlated electron systems, researchers often use quasi-one-dimensional (1D) ladder subsystems as proxies. However, existing ladder models used in the literature suffer from two critical flaws when attempting to model the 2D $\text{CuO}_2$ lattice:

1. Incorrect Stoichiometry: Common ladder structures (e.g., 3-chain and 5-chain ladders) have Copper-to-Oxygen (Cu:O) ratios of 2:3 and 2:5, respectively, deviating significantly from the 1:2 ratio of the actual 2D lattice. This prevents accurate study of charge redistribution between Cu and O atoms.
2. Symmetry Constraints: The "4-chain tube" model has the correct 1:2 ratio but only exhibits a superconducting (Luther-Emery) phase if the hopping between Cu and O orbitals is anisotropic (stronger along the legs than the rungs), which contradicts the rotational symmetry of the 2D lattice.

The authors aim to resolve these issues by proposing new ladder structures that preserve the correct Cu:O ratio and investigating whether they can reproduce the essential physics of doped cuprates, specifically the Luther-Emery phase (characterized by enhanced pairing correlations) and the charge distribution between orbitals.

2. Methodology

• Model: The authors study the Emery model on two-leg ladders. The Hamiltonian includes:
  • Hopping terms: $t_{dp_x}$ (Cu-O along legs), $t_{dp_y}$ (Cu-O along rungs), and $t_{pp}$ (O-O).
  • On-site energies: $\epsilon$ (energy difference between $p$ and $d$ orbitals).
  • Coulomb repulsion: $U_d$ (on Cu) and $U_p$ (on O).
• Lattice Structures: Three novel ladder structures are introduced, each acting as a supercell of the 2D $\text{CuO}_2$ lattice, ensuring a Cu:O ratio of 1:2:
  1. t-ladder: Preserves translational symmetry but breaks reflection symmetry.
  2. g-ladder: Possesses glide symmetry (unit cell is twice as large).
  3. r-ladder: Preserves reflection symmetry but has a doubled periodicity along the legs.
  • Comparison: These are compared against common structures (3-chain, 5-chain, 4-chain tube).
• Numerical Technique: The Density Matrix Renormalization Group (DMRG) method is used to calculate ground-state properties.
  • System sizes up to $L=40$ rungs were simulated.
  • Open boundary conditions were used with oxygen sites at the ends to minimize boundary effects.
  • Bond dimensions up to 9000 were employed to ensure high accuracy.
• Parameters: Parameters were optimized to maximize the Pair Binding Energy (PBE) for two doped holes while ensuring the undoped system is a charge-transfer insulator with a realistic charge ratio ($\approx$ 3:7 for O:Cu).

3. Key Contributions

1. Proposal of New Ladder Geometries: The authors define and validate three new ladder structures (t, g, r) that strictly preserve the 1:2 Cu:O atomic ratio, a prerequisite for studying charge redistribution phenomena accurately.
2. Validation of Luther-Emery Phase: They demonstrate that these new ladders, despite having lower symmetries than common models, robustly exhibit a Luther-Emery liquid phase upon doping. This phase is characterized by a gapless charge sector and a gapped spin sector, along with enhanced pairing correlations.
3. Isotropic Hopping Verification: Crucially, the authors show that the Luther-Emery phase persists even when hopping is isotropic ($t_{dp_x} = t_{dp_y}$), a condition required for these ladders to be valid supercells of the 2D lattice.
4. Charge Distribution Analysis: The preservation of the Cu:O ratio allows for a systematic investigation of how charge distributes between Cu ($d$) and O ($p$) orbitals as a function of doping and Hamiltonian parameters, directly linking microscopic parameters to macroscopic superconducting properties.

4. Key Results

A. Undoped State (Insulating Phase)

• For parameters listed in Table II, all three new ladders are charge-transfer insulators with a finite charge gap ($E_c$) and a smaller but finite spin gap ($E_s$).
• The ground state exhibits short-range antiferromagnetic order on the Cu sites.
• The charge gap scales correctly with system size, confirming the insulating nature in the thermodynamic limit.

B. Doped State (Luther-Emery Phase)

• Gap Behavior: Upon doping (both hole and electron), the charge gap closes ($E_c \to 0$), while the spin gap remains finite. This is the signature of a Luther-Emery liquid.
• Pair Binding Energy (PBE): The PBE is positive at low doping levels, indicating that doped carriers form bound pairs. The PBE converges to finite values as system size increases.
• Pairing Correlations: Pairing correlation functions decay as $r^{-\alpha}$ with $\alpha \le 1$ (slower than the $r^{-2}$ decay of non-interacting systems), indicating quasi-long-range pairing order.
• Isotropy: Positive PBE and enhanced correlations are observed even when $t_{dp_x} = t_{dp_y}$, confirming the models' validity as 2D supercells.

C. Charge Distribution and Pairing Strength

• Doping Asymmetry:
  • Hole doping: Holes primarily occupy Oxygen ($p$) orbitals.
  • Electron doping: Holes are removed primarily from Copper ($d$) orbitals.
  • This behavior confirms the charge-transfer nature of the insulator.
• Correlation with PBE: The authors find a systematic relationship between the oxygen occupation ($2n_p$) and the Pair Binding Energy.
  • The PBE exhibits a maximum in the range $0.3 \lesssim 2n_p \lesssim 0.4$.
  • This parabolic-like dependence mirrors the experimental "dome" shape of the critical temperature ($T_c$) versus oxygen content in real cuprates (e.g., $\text{YBa}_2\text{Cu}_3\text{O}_{6+\delta}$).
• Parameter Sensitivity:
  • Increasing $U_d$ increases PBE for hole doping but decreases it for electron doping.
  • Decreasing the energy difference $\epsilon$ increases the maximal PBE but shifts the optimal oxygen occupation to higher values.

5. Significance

This work provides a superior theoretical framework for studying high-temperature superconductivity using 1D ladder approximations.

• Resolution of Stoichiometry Issues: By correcting the Cu:O ratio, the new models allow for a direct and accurate comparison with experimental data regarding charge distribution, which was previously impossible with common ladder models.
• Bridging Theory and Experiment: The ability to map the calculated charge distribution ($n_d$ vs. $n_p$) and pairing strength to experimental observations (like the $T_c$ dome) suggests that material-specific parameters for the Emery model can be determined more reliably.
• Validation of 2D Physics: The confirmation that these lower-symmetry ladders still capture the essential Luther-Emery physics (even with isotropic hopping) strengthens the validity of using 1D DMRG to infer properties of 2D cuprates.

In conclusion, the authors argue that the t, g, and r-ladder structures constitute a more robust approach than previously used models for investigating the interplay between charge redistribution, doping, and pairing in the Emery model.