The role of the apical oxygen in cuprate high-temperature superconductors

By combining density-functional theory and cluster dynamical mean-field theory, this study demonstrates that while apical oxygen displacement modulates the superconducting order parameter in cuprates, this effect is primarily driven by changes in effective hole doping rather than the charge-transfer gap, urging caution in interpreting correlations between $T_c$ and apical oxygen distance across different compounds.

The Gist

Imagine a high-temperature superconductor (a material that conducts electricity with zero resistance) as a bustling, crowded dance floor. The dancers are electrons, and for the music (electricity) to flow perfectly without anyone tripping or bumping into each other, the dancers need to pair up and move in perfect sync. These pairs are called Cooper pairs.

For decades, scientists have been trying to figure out exactly what makes the dance floor "perfect." One specific feature of the building's architecture has been a major suspect: the Apical Oxygen.

The Suspect: The "Ceiling Light"

In these copper-oxide materials, the dance floor is a flat layer of copper and oxygen atoms. But sticking up from this floor, like a chandelier hanging from the ceiling, is an extra oxygen atom. This is the Apical Oxygen.

Scientists noticed something strange: when you move this "chandelier" up or down (changing its distance from the floor), the quality of the dancing changes.

• The Old Theory: Some researchers thought moving the chandelier changed the "energy gap" (the difficulty of getting a new dancer onto the floor). They believed that moving the chandelier up made it easier for dancers to pair up, leading to better superconductivity.
• The New Discovery: This paper says, "Hold on, let's look closer."

The Investigation: A Digital Simulation

The authors of this paper are like master architects who built a perfect, virtual 3D model of these materials. They didn't just guess; they used powerful computers to simulate what happens when they physically move that "chandelier" (the apical oxygen) up and down, while keeping everything else exactly the same.

They tested three different types of superconducting buildings (Bi-2201, Bi-2212, and Hg-1201) to see if the results were consistent.

The Big Reveal: It's Not the Energy Gap, It's the Crowd Density

Here is the twist they found, explained simply:

1. The "Energy Gap" Myth: They proved that moving the chandelier does not significantly change the energy gap. The old idea that the chandelier controls the "difficulty of entry" for electrons was wrong.
2. The Real Culprit: The Crowd (Doping): What actually changes is the number of dancers on the floor.
   • Think of the superconductor as a dance floor that is slightly too crowded (overdoped). The dancers are bumping into each other, and the pairs aren't forming well.
   • When the authors moved the apical oxygen (the chandelier), it acted like a subtle adjustment to the building's plumbing. It changed how many electrons were flowing into the dance floor from the rest of the building.
   • In the Bi-based materials: Moving the chandelier up actually reduced the number of dancers on the floor, bringing the crowd closer to the "perfect density" where pairing is easiest.
   • In the Hg-based material: The effect was the opposite because the starting crowd density was different, but the mechanism was the same: the chandelier movement changed the crowd size.

The Analogy: The Concert Hall

Imagine a concert hall where the acoustics (superconductivity) are perfect only if the audience size is exactly 1,000 people.

• The Old View: Scientists thought moving the stage lights (apical oxygen) changed the sound quality of the room itself.
• The New View: This paper shows that moving the lights actually changed the number of people who decided to buy tickets and enter the hall.
  • If the hall was too full (overcrowded), moving the lights up let some people leave, bringing the crowd down to the sweet spot of 1,000.
  • If the hall was too empty, moving the lights might have encouraged more people to enter.

Why Does This Matter?

For a long time, scientists looked at different superconducting materials and said, "Oh, this one has a taller chandelier and a higher superconducting temperature, so the chandelier height must be the secret!"

This paper says: Be careful.
The change in the "chandelier" height only causes a small improvement in superconductivity (about a 20-30% change in the pairing strength). It is not the magic bullet. The real magic is that the chandelier height subtly controls how many electrons are on the copper layer.

The Takeaway

The role of the apical oxygen isn't to change the rules of the dance; it's to act as a dimmer switch for the crowd size. By understanding this, scientists can stop looking for a single "magic distance" and start focusing on how to precisely control the number of electrons (doping) in these materials to create even better superconductors for the future.

In short: It's not about where the ceiling light is; it's about how many people are dancing underneath it.

Technical Summary

1. Problem Statement

The paper addresses a long-standing unresolved question in condensed matter physics: What is the causal relationship between the apical oxygen distance ($\delta_{api}$) and superconductivity in cuprates?

• Context: Recent Scanning Tunneling Microscopy (STM) experiments on optimally doped $\text{Bi}_2\text{Sr}_2\text{CaCu}_2\text{O}_{8+\delta}$ (Bi-2212) revealed a correlation between the natural superstructure modulation and the superconducting order parameter ($m_{SC}$).
• The Hypothesis: Previous interpretations (e.g., O'Mahony et al., 2022) suggested that the modulation of $\delta_{api}$ alters the Charge-Transfer Gap (CTG). Specifically, it was hypothesized that increasing $\delta_{api}$ reduces the CTG, thereby enhancing superconductivity.
• The Gap: While correlations between $\delta_{api}$ and the critical temperature ($T_c$) have been observed across different cuprate families, it remains unclear if $\delta_{api}$ is the direct driver of $T_c$ variations or if other structural/chemical factors are dominant. Furthermore, previous theoretical attempts to isolate the effect of $\delta_{api}$ lacked quantitative agreement with experiment or relied on questionable approximations (e.g., ignoring hole doping or using single-site DMFT).

2. Methodology

The authors employ a rigorous first-principles ($ab$ $initio$) framework combining Density Functional Theory (DFT) with Cluster Dynamical Mean-Field Theory (CDMFT).

• Materials Studied:
  • $\text{Bi}_2\text{Sr}_2\text{CuO}_{6+\delta}$ (Bi-2201, single-layer)
  • $\text{Bi}_2\text{Sr}_2\text{CaCu}_2\text{O}_{8+\delta}$ (Bi-2212, bilayer)
  • $\text{HgBa}_2\text{CuO}_{4+\delta}$ (Hg-1201, single-layer)
• Computational Strategy:
  • Structural Modeling: Instead of simulating the complex long-range supermodulation (which is computationally intractable), the authors constructed homogeneous tetragonal unit cells with four distinct values of $\delta_{api}$, matching experimental ranges.
  • Optimization: Lattice parameters ($a, c$) and internal atomic positions were optimized using the eDMFT approach while keeping $\delta_{api}$ fixed.
  • Electronic Structure: Calculations were performed using the DFT+CDMFT framework (extending Wien2k with PyQCM).
    • Correlated subspace: $2 \times 2$ plaquette of Cu-$d_{x^2-y^2}$ orbitals.
    • Interaction strength: $U = 9$ eV (with $J_{Hund} = 1$ eV).
    • Solver: Exact diagonalization.
  • Self-Doping: The natural self-doping of Bi-based cuprates was treated explicitly without using the Virtual Crystal Approximation (VCA). For Hg-1201, VCA was used to fix doping in the underdoped regime.
• Analysis Tools:
  • Wannier Functions: Used to downfold the DFT bands into effective three-band (Emery) and single-band (Hubbard) models to analyze hopping parameters ($t_{pd}, t_{pp}$), on-site energies, and the CTG.
  • Superexchange ($J$): Calculated via exact diagonalization of $\text{Cu}_2\text{O}_2$ clusters.

3. Key Contributions & Results

A. Quantitative Reproduction of Experimental Data

The authors successfully reproduced the experimental STM observations of O'Mahony et al.

• Result: The calculated relative variations of the superfluid density ($|m_{SC}|^2$) as a function of $\delta_{api}$ show remarkable quantitative agreement (approx. 30% variation) with experimental data for Bi-2212.
• Magnitude: While the relative variation is significant, the absolute variation in $m_{SC}$ induced solely by $\delta_{api}$ is modest (approx. $0.5 \times 10^{-2}$). This is significantly smaller than the differences observed between different cuprate families (e.g., Bi-2201 vs. Bi-2212).

B. Refutation of the Charge-Transfer Gap (CTG) Mechanism

The paper definitively rules out the hypothesis that $\delta_{api}$ modulates superconductivity via the CTG.

• Finding: As $\delta_{api}$ increases, the CTG increases (or remains stable), rather than decreasing.
  • In the effective three-band model, the energy separation between Cu-$d$ and O-$p$ orbitals increases with $\delta_{api}$.
  • The superexchange interaction $J$ was found to decrease with increasing $\delta_{api}$. Since $J \propto t^2/\Delta_{CTG}$, a decreasing $J$ is incompatible with the observed increase in $m_{SC}$ if the CTG were the primary driver.
• Conclusion: The scenario proposed in Ref. [20] (that increasing $\delta_{api}$ reduces CTG to enhance $T_c$) is not viable.

C. Identification of the True Mechanism: Effective Hole Doping

The authors identify modulation of effective hole doping in the $\text{CuO}_2$ planes as the dominant mechanism.

• Mechanism: Increasing $\delta_{api}$ alters the charge transfer between the $\text{CuO}_2$ planes and the charge reservoir layers.
  • In Bi-2201/2212 (Overdoped): These materials are naturally overdoped. Increasing $\delta_{api}$ reduces the effective hole doping, driving the system closer to the optimal doping regime, thereby increasing $m_{SC}$.
  • In Hg-1201 (Underdoped): Calculations on underdoped Hg-1201 show the opposite trend: increasing $\delta_{api}$ increases effective hole doping, which also increases $m_{SC}$ (as it moves the system toward optimal doping).
• Evidence: The variation in electron occupation ($\Delta n_{\text{CuO}_2}$) correlates perfectly with the changes in $m_{SC}$. The spectral weight changes observed in STM are attributed to shifts in the Fermi level and occupation, not the CTG.

D. Distinction Between "Within-Compound" and "Between-Compound" Trends

The study highlights a critical distinction:

• Within a compound: Variations in $\delta_{api}$ modulate $m_{SC}$ primarily by tuning the effective doping toward the optimal point.
• Between compounds: Correlations between $\delta_{api}$ and $T_c$ across different families (e.g., La-based vs. Bi-based) cannot be attributed solely to $\delta_{api}$. The absolute magnitude of the effect of $\delta_{api}$ is small compared to other structural differences (e.g., single-layer vs. bilayer).

4. Significance

• Resolution of a Controversy: The paper resolves the debate regarding the O'Mahony et al. STM findings, attributing the observed correlations to doping modulation rather than CTG changes.
• Validation of $Ab$ $Initio$ Frameworks: It demonstrates that modern DFT+CDMFT approaches can quantitatively predict superconducting order parameters and resolve specific structural degrees of freedom in correlated oxides.
• Guidance for Material Design: The work emphasizes that optimizing $T_c$ requires careful control of effective hole doping. While apical oxygen distance is a lever to tune doping, it is not a direct "knob" for the CTG.
• Caveat on Cross-Compound Comparisons: The authors warn against inferring direct causal links between $\delta_{api}$ and $T_c$ when comparing different cuprate families, as the effect of $\delta_{api}$ is secondary to other structural and chemical factors.

In summary, the paper establishes that the apical oxygen distance influences superconductivity in cuprates indirectly by modulating the effective hole doping of the $\text{CuO}_2$ planes, rather than by directly altering the charge-transfer gap.