Perpendicular electric field induced $s^\pm$-wave to $d$-wave superconducting transition in thin film La$_3$Ni$_2$O$_7$

Using dynamical cluster quantum Monte Carlo calculations on a two-orbital bilayer Hubbard model, this study demonstrates that a perpendicular electric field induces a superconducting transition from $s^\pm$-wave to $d$-wave symmetry in thin-film La$_3$Ni$_2$O$_7$ by suppressing $d_{z^2}$-orbital pairing and enhancing $d_{x^2-y^2}$-orbital contributions, with the resulting $d$-wave pairing exhibiting a dome-like dependence on field strength.

The Gist

Imagine you have a very special, two-story sandwich made of a rare material called La₃Ni₂O₇. In the world of physics, this sandwich is famous because it can conduct electricity with zero resistance (superconductivity) when you squeeze it tight with high pressure. Scientists are desperate to figure out how to make it superconduct at normal pressure and even higher temperatures, because that would revolutionize our power grids and electronics.

This paper is like a recipe experiment. The scientists asked: "What happens if we zap this sandwich with an electric field from the top, instead of just squeezing it?"

Here is the story of what they found, explained simply:

1. The Two "Flavors" of Electrons

Inside this sandwich, the electrons (the tiny particles carrying electricity) live in two different "rooms" or orbitals:

• The $d_{z^2}$ Room: Think of this as the interior decorator. It's great at connecting the top floor to the bottom floor. When the sandwich is just sitting there (no electric field), these electrons hold hands across the layers in a specific way called $s_{\pm}$-wave. It's a stable, cozy handshake that allows superconductivity to happen.
• The $d_{x^2-y^2}$ Room: Think of this as the acrobatic gymnast. It's very good at moving quickly within a single floor but isn't as interested in connecting the two floors.

2. The Electric Field "Zap"

The researchers applied a perpendicular electric field (a zap from top to bottom). Imagine this like tilting the sandwich so gravity pulls everything toward the bottom floor.

• What happened to the Interior Decorator ($d_{z^2}$)?
  The electric field messed up the balance between the top and bottom floors. The "handshake" between the layers got shaky and weak. The original superconducting style ($s_{\pm}$-wave) started to fade away.

• What happened to the Gymnast ($d_{x^2-y^2}$)?
  As the electric field pushed electrons from the bottom floor up to the top, the "gymnast" electrons got more active. Suddenly, they started forming a new kind of superconducting handshake called $d$-wave. This is a different style of dancing, one that happens mostly on the top floor.

3. The "Goldilocks" Zone (The Dome Shape)

Here is the most interesting part. The scientists found that the new $d$-wave superconductivity didn't just get stronger and stronger as they increased the electric field.

Instead, it behaved like a dome or a hill:

• Too little zap: The gymnasts aren't active enough. No superconductivity.
• Just the right zap: The electrons are perfectly balanced. The superconductivity peaks! This is the "Goldilocks" zone.
• Too much zap: You push too many electrons around, and the system gets chaotic. The superconductivity collapses again.

It's like tuning a radio: if you turn the dial too far left or too far right, you get static. But right in the middle, the music is crystal clear.

4. The Twist: Doping (Adding Extra Ingredients)

The scientists also tested what happens if they add extra electrons (electron-doping) or remove some (hole-doping).

• Hole-doping (removing electrons): The "Goldilocks" zone shifted. You needed a stronger electric zap to get the peak performance, but the superconductivity was still there.
• Electron-doping (adding electrons): The gymnasts got too crowded. The new $d$-wave style didn't work at all, and the original style was already weak. The sandwich just wouldn't superconduct well with this method.

The Big Picture

This paper tells us that electric fields are a powerful switch for this material. By applying the right amount of electric "zap," you can force the electrons to change their dance style from the old, layer-connecting style to a new, high-energy style.

Why does this matter?
It suggests that we might not need massive, expensive pressure machines to make these materials superconduct. Instead, we might be able to use simple electric fields (like in a transistor) to tune the material and make it superconduct at higher temperatures. It's like realizing you don't need to squeeze a sponge to get water out; sometimes, you just need to give it the right little push.

In short: The scientists discovered a way to "tune" a superconductor with electricity, finding a sweet spot where a new type of superconductivity shines brighter than the old one, but only if you don't push it too hard.

Technical Summary

1. Problem Statement

The discovery of high-temperature superconductivity (SC) in pressurized Ruddlesden–Popper (RP) bilayer nickelate La$_3$Ni$_2$O$_7$ has sparked intense debate regarding its microscopic mechanism. Two primary theories exist:

1. $d_{z^2}$-dominated: Pairing originates from interlayer spin fluctuations in the $d_{z^2}$ orbitals, with coherence established via $d_{x^2-y^2}$–$d_{z^2}$ hybridization.
2. $d_{x^2-y^2}$-dominated: Hund's coupling transfers magnetic correlations to the $d_{x^2-y^2}$ orbitals, driving $s_{\pm}$-wave pairing.

A critical open question is the nature of the pairing symmetry ($s$-wave vs. $d$-wave) and how external parameters, specifically a perpendicular electric field ($E$), influence the superconducting state in thin films. While weak-coupling methods (e.g., slave-boson mean-field, RPA) have predicted a transition from $s_{\pm}$ to $d$-wave under an electric field, strong-coupling many-body effects in the presence of an electric field had not been rigorously verified.

2. Methodology

The authors employed Dynamical Cluster Approximation (DCA) combined with Continuous-Time Auxiliary-Field (CT-AUX) Quantum Monte Carlo (QMC) simulations.

• Model: An imbalanced bilayer two-orbital Hubbard model on a 2D square lattice.
  • Orbitals: $d_{x^2-y^2}$ (denoted $x$) and $d_{z^2}$ (denoted $z$).
  • Layers: Top ($t$) and Bottom ($b$) layers.
  • Electric Field Simulation: Modeled by introducing an on-site energy difference ($\epsilon_b \neq 0$) on the bottom layer relative to the top layer ($\epsilon_t = 0$), creating a potential drop $E = (\epsilon_b - \epsilon_t)/e$.
  • Parameters: Includes nearest-neighbor, next-nearest-neighbor, and inter-layer hoppings. Coulomb interactions ($U=4.0$ eV, $U'=2.4$ eV, Hund's coupling $J=0.8$ eV) are set to typical values for nickelates.
• Cluster Size: Primarily used an $N_c = 8$ ($4 \times 2$) cluster to capture $d$-wave pairing symmetry and access lower temperatures. A smaller $N_c = 2$ cluster was used for pair-field susceptibility analysis to reach even lower temperatures.
• Superconducting Instability: Determined by solving the Bethe-Salpeter Equation (BSE) in the particle-particle channel. The leading eigenvalue $\lambda_\alpha(T)$ (where $\alpha = s_{\pm}, d$) was extrapolated to find the critical temperature $T_c$.
  • $s_{\pm}$-wave form factor: $g_{s\pm}(k) = \cos k_z$.
  • $d$-wave form factor: $g_d(k) = \cos k_x - \cos k_y$.

3. Key Contributions

• Strong-Coupling Verification: Provided the first large-scale many-body QMC confirmation of the electric-field-induced pairing symmetry transition from $s_{\pm}$-wave to $d$-wave, validating previous weak-coupling predictions.
• Orbital Selectivity: Demonstrated that the electric field suppresses the $d_{z^2}$-driven interlayer $s_{\pm}$ pairing while enhancing $d_{x^2-y^2}$-driven intralayer $d$-wave pairing.
• Dome Behavior: Identified a dome-shaped dependence of the $d$-wave $T_c$ on the electric field strength, linked to the optimal filling of the top-layer $d_{x^2-y^2}$ orbital.
• Doping Dependence: Systematically compared undoped, hole-doped, and electron-doped regimes, revealing distinct responses to the electric field.

4. Key Results

A. Undoped Case ($n=3.0$)

• Symmetry Transition: At $E=0$, the system exhibits dominant $s_{\pm}$-wave pairing driven by interlayer $d_{z^2}$ correlations. As $E$ increases, the interlayer $s_{\pm}$ pairing is suppressed.
• Emergence of $d$-wave: At intermediate fields ($E \approx 0.5 - 1.0$ V), a transition occurs to $d$-wave pairing dominated by the $d_{x^2-y^2}$ orbital (specifically in the top layer).
• $T_c$ Evolution:
  • $s_{\pm}$-$T_c$ decreases monotonically with increasing $E$.
  • $d$-wave $T_c$ exhibits a dome-like behavior, peaking at $E \approx 0.5$ V, where it exceeds the $s_{\pm}$-$T_c$.
• Mechanism: The electric field drives electrons from the bottom layer to the top layer. This creates a density imbalance that weakens interlayer $d_{z^2}$ coupling (suppressing $s_{\pm}$) but optimizes the filling of the top-layer $d_{x^2-y^2}$ orbital for $d$-wave pairing.

B. Hole-Doped Case ($n=2.8$)

• Enhanced $T_c$: Hole doping generally increases the $T_c$ of the $s_{\pm}$ state compared to the undoped case.
• Shifted Optimal Field: The $d$-wave $T_c$ dome persists but shifts to a higher optimal electric field ($E \approx 1.0$ V).
• Mechanism: The shift is attributed to the need for a specific top-layer $d_{x^2-y^2}$ filling ($n^t_x \approx 0.68$) to maximize $d$-wave pairing, which requires a stronger field to achieve at this doping level.

C. Electron-Doped Case ($n=3.2$)

• Suppressed $d$-wave: No significant $d$-wave pairing is observed. The $T_c$ remains low and nearly constant regardless of the electric field.
• Stable $s_{\pm}$: The $s_{\pm}$-$T_c$ decreases very slowly with increasing $E$.
• Reasoning: The electron doping prevents the system from reaching the optimal filling for the top-layer $d_{x^2-y^2}$ orbital required for $d$-wave instability, even under strong electric fields.

D. Orbital-Resolved Susceptibility (Low-Temperature Analysis)

Using a smaller cluster ($N_c=2$) to access lower temperatures, the authors analyzed the pair-field susceptibility.

• They found that the $d_{x^2-y^2}$ orbital actually possesses a stronger pairing instability than the $d_{z^2}$ orbital in both undoped and doped cases at low temperatures.
• This corrects the high-temperature observation where $d_{z^2}$ appeared dominant, suggesting that the $d_{x^2-y^2}$ orbital is the primary driver of superconductivity in the low-$T$ limit, consistent with the $d$-wave transition.

5. Significance and Outlook

• Mechanism Insight: The study provides strong evidence that the superconducting mechanism in RP nickelates is not static; it is highly tunable via external fields. The dominant orbital can shift from $d_{z^2}$ to $d_{x^2-y^2}$ depending on the electric field and doping.
• Experimental Relevance: The predicted dome-like behavior of $T_c$ versus electric field offers a concrete experimental signature for thin-film La$_3$Ni$_2$O$_7$ devices. It suggests that applying a gate voltage could potentially optimize $T_c$ or switch the pairing symmetry.
• Theoretical Validation: The results bridge the gap between weak-coupling theories and strong-coupling many-body physics, confirming that the interplay between interlayer mismatch and orbital filling is crucial for the SC state.
• Limitations: The study is limited by the fermionic sign problem in QMC, restricting the lowest accessible temperatures. The authors note that further investigation with advanced methods is needed to definitively determine the exact $T_c$ values and the precise role of the $d_{x^2-y^2}$ orbital at the lowest temperatures.

In conclusion, the paper establishes that a perpendicular electric field acts as a powerful control knob in La$_3$Ni$_2$O$_7$ thin films, capable of inducing a transition from $s_{\pm}$-wave to $d$-wave superconductivity, with the $d$-wave state exhibiting a non-trivial dome-shaped dependence on field strength driven by orbital-selective electron redistribution.