A mechanism for nonmonotonic $T_{c,max}(n)$ in multilayer cuprates

This paper proposes a preformed pair mechanism to explain the nonmonotonic dependence of the maximal critical temperature on the number of conducting layers in cuprates, attributing the peak to an optimal balance between pair compactness and lightness that is disrupted by rising kinetic energy in thicker multilayers.

The Gist

The Big Picture: The "Goldilocks" Problem of Superconductors

Imagine you are trying to build the ultimate super-highway for electricity. In a normal wire, electricity flows like cars in rush hour traffic—lots of bumps, stops, and wasted energy (heat). In a superconductor, the cars (electrons) pair up and glide perfectly without any friction.

The holy grail of physics is to make these super-highways work at room temperature (around 20°C or 68°F). Right now, the best we have works at about -135°C (138 K). That's cold, but not room temperature.

Scientists have noticed a weird pattern in these materials (called cuprates). If you stack more layers of the "highway" on top of each other, the temperature at which superconductivity works goes up, hits a peak, and then starts to go back down.

• 1 layer: Good.
• 2 layers: Better.
• 3 layers: The BEST.
• 4+ layers: It gets worse again.

Why does adding more layers eventually hurt the performance? This paper by Pavel Kornilovitch tries to solve that mystery.

---

The Core Idea: The "Dancing Pair" Analogy

To understand this, imagine the electrons aren't just cars; they are dancers.

• The Goal: For superconductivity to happen, these dancers must form pairs and move in perfect unison (this is called a "condensate").
• The Rule: To get the most dancers on the floor at once (which creates the strongest superconductivity), the dancers need to be light (easy to move) and compact (taking up little space).

The paper argues that the number of layers ($n$) acts like a tuning knob that changes how heavy and how big these dancing pairs are.

1. From 1 to 3 Layers: Getting Lighter

Imagine you are trying to dance in a narrow hallway (1 layer). It's hard to move sideways because the walls are close. You feel heavy and slow.

• Adding a second layer is like opening a second hallway right next to the first. Now you can hop between them easily. You feel lighter and faster.
• Adding a third layer makes it even easier to move up and down. The dancers become very "light" (low mass).
• Result: Because they are lighter, they can condense into a super-state at a higher temperature. $T_c$ goes up.

2. From 3 to 16 Layers: Getting Too Big

Now, imagine you keep adding layers until you have a massive skyscraper of hallways (16 layers).

• The Problem: The "attraction" that holds the dancers together (the glue) is fighting against their natural desire to run around (kinetic energy).
• As you add more layers, the dancers get so much energy from running around that the "glue" starts to stretch. The dancers stop holding hands tightly and start drifting apart.
• The Result: The pairs become huge and fluffy (inflated volume). Even though they are light, they are so big that you can't fit many of them on the dance floor. They bump into each other and break the rhythm.
• $T_c$ goes down.

The "Sweet Spot" (The Goldilocks Zone)

The paper explains that there is a perfect balance between being light and being compact.

• Too Anisotropic (1 layer): The dancers are stuck in one lane. They are heavy and slow.
• Too Isotropic (Many layers): The dancers are free to run everywhere, but they get so excited they drift apart and become huge balloons.
• Just Right (3 layers): The dancers are light enough to move fast, but the "glue" is still strong enough to keep them tight and compact.

This is why the peak is usually at $n=3$.

The "Apical Oxygen" Twist: Why 3 is Special

The author adds a clever detail about the architecture of these materials.

• In a 3-layer stack, the outer layers are close to the "ceiling" and "floor" of the building (called apical oxygens).
• These ceiling/floor atoms act like magnets that pull the dancers down, making them feel even lighter and helping them stay tight.
• In a 4 or 5-layer stack, you have inner layers that are far away from the magnets. They don't get the same "help." The system gets messy, and the balance is lost.

What Does This Mean for the Future?

The paper isn't just explaining the past; it's a recipe book for the future.

If we want to break the current record (138 K) and reach room temperature, we need to:

1. Keep the layers thin: Don't just keep stacking them forever; find the perfect number (maybe 3, maybe 4, depending on the material).
2. Tweak the "Ceiling": We need to adjust the atoms above and below the layers to make the dancers even lighter without making them drift apart.
3. Squeeze the Glue: We need to make the attraction between dancers stronger so they stay compact even when they are running fast.

Summary in One Sentence

Superconductivity works best when the electron pairs are light enough to move fast but small enough to fit together; adding too many layers makes the pairs get too big and fluffy, ruining the party, so the "Goldilocks" number of layers is usually three.

Technical Summary

1. Problem Statement

The paper addresses a long-standing puzzle in high-temperature superconductivity: the non-monotonic dependence of the maximal critical temperature ($T_{c,max}$) on the number of conducting $\text{CuO}_2$ layers ($n$) in the unit cell.

• Observation: In most copper-oxide families (e.g., Hg, Tl, Bi), $T_{c,max}$ increases from $n=1$ to a peak at $n=3$ or $n=4$, then decreases for $n > 4$.
• The Anomaly: The "infinite-layer" compound $(\text{Sr,Ca})\text{CuO}_2$ ($n=\infty$) has a high $T_c \approx 110$ K. However, inserting charge reservoir layers to create finite-$n$ compounds generally lowers $T_c$, yet specific $n=3$ and $n=4$ members (like $\text{HgBa}_2\text{Ca}_2\text{Cu}_3\text{O}_{8+\delta}$) exceed this baseline, reaching up to 138 K.
• Gap: Previous theories (phenomenological coupling, Coulomb energy reduction, magnetic exchange, or phonon-induced pairing) failed to consistently explain the non-monotonic peak and the subsequent drop.

2. Methodology

The author adopts the preformed pair mechanism (also known as real-space pairing or Bose-Einstein Condensation (BEC) of local pairs). The approach combines qualitative physical arguments with rigorous exact solutions of a two-body quantum mechanical problem.

• Theoretical Framework:
  • $T_{c,max}$ is modeled using the anisotropic BEC formula for a gas of preformed pairs:
    $$T_{c,max} \propto \frac{1}{(m_x^* m_y^* m_z^*)^{1/3} \Omega_p^{2/3}}$$
    where $m_i^*$ are the effective masses of the pair and $\Omega_p$ is the pair volume.
  • High $T_c$ requires pairs that are simultaneously light (low effective mass) and compact (small volume).
• The Model:
  • A multilayer attractive Hubbard model is constructed for $n$ layers.
  • Hamiltonian: Includes in-plane hopping ($t$), inter-stack hopping ($t'$), and on-site attractive interaction ($-|U|$).
  • Two Models Tested:
    1. Uniform Model: All layers have identical on-site energy ($\epsilon_\alpha = 0$).
    2. "Low-Outer-Planes" Model: Outer layers have a lower on-site energy ($\epsilon_1 = \epsilon_n = \Delta < 0$) due to polaron shifts from apical oxygens, while inner layers have $\epsilon = 0$.
• Calculation Technique:
  • The two-hole Schrödinger equation is solved exactly in real space for $n=1$ to $5$.
  • The solution involves transforming to momentum space, reducing the problem to a system of linear equations for auxiliary integrals (Green's functions), and solving for bound-state energy, effective masses ($m^*$), and pair radii ($r^*$).
  • Pair volume $\Omega_p$ is regularized to prevent unphysical collapse at strong coupling.

3. Key Contributions and Physical Mechanisms

The paper identifies the competition between kinetic energy and pairing strength as the driver of the $T_{c,max}(n)$ curve.

• Role of Anisotropy:
  • $n=1$ to $n=2$: Adding a second layer reduces lattice anisotropy (increases $t'$ relative to the stack). This lowers the out-of-plane effective mass ($m_z^*$), making pairs lighter and increasing $T_c$. This rise is universal across cuprate families.
  • $n \geq 3$: As $n$ increases further, the kinetic energy ($K$) of the holes increases, while the attractive strength ($V$) remains relatively constant. The ratio $V/K$ decreases, approaching the pairing threshold.
• The Trade-off:
  • As the system approaches the pairing threshold (weak binding), the pair volume ($\Omega_p$) balloons (pairs become diffuse).
  • Simultaneously, the out-of-plane mass ($m_z^*$) behavior becomes complex. In the uniform model, $m_z^*$ increases for $n \geq 3$ due to the localization of the pair wavefunction on the inner planes (which are separated by tunneling barriers).
• The "Low-Outer-Planes" Solution:
  • The uniform model predicts a peak at $n=2$, contradicting experiments.
  • Introducing the energy shift $\Delta < 0$ (mimicking apical oxygen polaron effects) pulls the pair wavefunction toward the outer planes.
  • This delocalizes the pair across the stack, significantly reducing $m_z^*$ for $n \geq 3$.
  • This reduction in mass counteracts the increase in pair volume, shifting the peak of $T_{c,max}$ to $n=3$ (or potentially higher depending on parameters), matching experimental data.

4. Results

• Quantitative Agreement: The "Low-Outer-Planes" model successfully reproduces the non-monotonic $T_{c,max}(n)$ curve with a peak at $n=3$, consistent with the Hg and Tl families (Fig. 11 vs. Fig. 1).
• Mechanism Validation:
  • $n=1 \to 2$: $T_c$ rises because $m_z^*$ drops (anisotropy reduction).
  • $n=2 \to 3$: $T_c$ continues to rise in the "Low-Outer-Planes" model because the reduction in $m_z^*$ (due to $\Delta$) outweighs the increase in pair volume.
  • $n > 3$: $T_c$ falls because the kinetic energy dominates, causing the pair volume $\Omega_p$ to expand rapidly (pairs become too large), reducing the close-packing density.
• Parameter Sensitivity: The position of the peak ($n_{max}$) is tunable. With specific parameter sets, the model can predict peaks at $n=2$, $n=3$, or even $n>3$, suggesting that different cuprate families are tuned to different points on this curve.

5. Significance and Implications

• Unification: The paper provides a unified explanation for the $T_c(n)$ dependence within the preformed pair framework, resolving the contradiction between the high $T_c$ of infinite-layer compounds and the peak at $n=3$.
• Design Rules for Higher $T_c$: The authors propose strategies to boost $T_c$ beyond the current record (138 K):
  1. Reduce $m_z^*$: Optimize charge reservoir layers to minimize out-of-plane mass without destroying attraction (e.g., via apical oxygen engineering).
  2. Compact Pairs: Increase the attractive interaction $|U|$ to reduce pair volume $\Omega_p$, thereby increasing the close-packing density, provided phase separation is avoided.
  3. Polaron Engineering: Utilize the "polaron undressing" effect (e.g., via pressure or THz excitation) to reduce effective mass dynamically.
• Theoretical Advance: The work demonstrates the necessity of exact two-body solutions in multiband systems to capture subtle effects like layer-dependent localization and mass renormalization, which mean-field approximations might miss.

In conclusion, the paper argues that the optimal $T_c$ in multilayer cuprates arises from a "Goldilocks" zone of anisotropy and interaction strength, where the copper-oxygen planes fine-tune the balance between kinetic energy and attractive forces to produce pairs that are both light enough to condense and compact enough to pack densely.