Raman response in superconducting multiorbital systems with application to nickelates

This paper investigates the electronic Raman scattering response in superconducting nickelates by analyzing various multiorbital models and pairing symmetries to identify characteristic fingerprints that can help determine the minimal model and pairing mechanism underlying high-$T_c$ superconductivity in these materials.

The Gist

Imagine you are a detective trying to solve a mystery inside a tiny, super-conductive city called a Nickelate. Recently, scientists discovered that this city can conduct electricity with zero resistance (superconductivity) at surprisingly high temperatures, especially when squeezed (pressure) or made very thin. But there's a catch: we don't know how the citizens (electrons) are holding hands to create this super-state.

This paper is like a forensic toolkit designed to help us figure out the "hand-holding" pattern (the superconducting gap symmetry) by listening to the city's "echoes."

Here is the breakdown of the paper using simple analogies:

1. The Mystery: What is the "Hand-Holding" Pattern?

In a superconductor, electrons pair up to move without friction. Think of these pairs as dance partners.

• The Gap: The energy required to break these dance partners apart.
• The Symmetry: The shape of the dance floor. Are they dancing in a circle (s-wave)? In a four-leaf clover shape (d-wave)? Or something more complex?

The scientists in the paper are trying to figure out which dance floor shape the Nickelate electrons are using.

2. The Tool: Raman Scattering (The "Echo Chamber")

To solve the mystery, the researchers use a technique called Raman Scattering.

• The Analogy: Imagine shining a flashlight (light) into a dark room full of bouncing balls (electrons). When the light hits the balls, they bounce back with a slightly different color (energy).
• The Echo: By analyzing the color of the light that bounces back, we can hear the "echo" of the electrons breaking apart. If the electrons are holding hands tightly, the echo happens at a specific energy level. If they are holding hands loosely, the echo is different.
• The Goal: The paper maps out exactly what these echoes should look like for different dance patterns.

3. The Suspects: Different Models of the City

The researchers didn't just look at one version of the Nickelate city. They tested three different "blueprints" to see which one matches reality:

• Blueprint A: The Single-Layer City (One Floor)
  Imagine a city with just one floor. The electrons live in two different "neighborhoods" (orbitals: $d_{x^2-y^2}$ and $d_{z^2}$).

  • The Finding: If the electrons dance in a clover shape (d-wave), the echo shows two distinct peaks (like two different drumbeats). If they dance in a circle (s-wave), there is one sharp peak.

• Blueprint B: The Double-Layer City (Two Floors)
  This is the more exciting one! Recent discoveries suggest the real Nickelate superconductors have two floors stacked on top of each other, and the electrons can jump between floors.

  • The Finding: When you add the second floor, the echoes get more complex. The "drumbeats" split apart. If the floors are tightly connected, the echoes look very different than if they are loosely connected.

• Blueprint C: The "One-Orbital" City (Simplified)
  Some scientists think only one neighborhood matters. The paper checks if this simple view works. They found that while it's a good start, it misses the subtle "whispers" that the two-neighborhood models catch.

4. The Big Twist: The "Additive" Trap

Here is the most important lesson from the paper.

For a long time, scientists used a shortcut to predict these echoes. They thought: "If I listen to Floor 1's echo and Floor 2's echo separately, and just add them together, I'll get the total sound."

• The Analogy: Imagine listening to a choir. The shortcut assumes the sound of the whole choir is just the sound of the tenors plus the sound of the sopranos added up.
• The Reality: The paper shows this is wrong for these complex cities. In the real Nickelate city, the floors talk to each other. The electrons on Floor 1 affect the electrons on Floor 2 in a way that creates new sounds that don't exist if you just add them up.
• The Warning: If you use the "add them up" shortcut, you might miss the most important clues or hear a fake clue. You need to listen to the whole city at once (the "Full Multiorbital" calculation).

5. Why Does This Matter?

Right now, there is a debate in the scientific community. Some say the Nickelate electrons are dancing in a circle (s-wave), others say a clover (d-wave), and others say a mix.

This paper provides a legend for the map.

• If you do the Raman experiment and see two peaks at low energy, it suggests a specific "clover" dance.
• If you see one sharp peak, it suggests a "circle" dance.
• If you see the peaks split in a specific way, it tells you the two floors are tightly connected.

The Takeaway

The authors are essentially saying: "We have built a detailed soundboard for the Nickelate superconductor. When experimentalists finally shine their light on these materials (which is happening now!), they can compare the echoes they hear against our soundboard. This will tell us exactly how the electrons are dancing, which is the key to understanding how to make even better superconductors for the future."

It's like giving the detectives a perfect audio recording of the crime scene so they can finally identify the culprit.

Technical Summary

1. Problem Statement

The recent discovery of high-temperature superconductivity ($T_c$) in pressurized bilayer nickelates ($La_3Ni_2O_7$) and thin-film infinite-layer nickelates has sparked intense debate regarding the underlying pairing mechanism. Key unresolved issues include:

• Gap Symmetry: Whether the pairing is $d$-wave, $s$-wave, or $s_{\pm}$-wave.
• Minimal Model: Whether superconductivity arises from a single active orbital ($d_{x^2-y^2}$) or requires a multiorbital description involving both $d_{x^2-y^2}$ and $d_{z^2}$ orbitals.
• Experimental Limitations: While Electronic Raman Scattering (ERS) is a powerful tool for determining gap size and symmetry, theoretical predictions for multiorbital systems are complex. A common approximation (Additive Response) treats bands independently, but its validity in strongly coupled multiorbital systems like nickelates is unverified.

The paper aims to provide a theoretical framework to interpret future Raman experiments on nickelates by calculating the Raman response for various minimal models and pairing symmetries, comparing full multiorbital calculations against the additive approximation.

2. Methodology

The authors employ a microscopic theoretical approach based on the Random Phase Approximation (RPA) and Nambu-Gorkov formalism.

• Models Investigated:

  1. Bilayer One-Orbital Model: Two coupled $d_{x^2-y^2}$ layers (relevant for infinite-layer nickelates or weakly coupled bilayers).
  2. Single-Layer Two-Orbital Model: A single layer with active $d_{x^2-y^2}$ and $d_{z^2}$ orbitals, including inter-orbital hopping.
  3. Bilayer Two-Orbital Model: Two layers of the above system, with strong interlayer hopping specifically between $d_{z^2}$ orbitals (relevant for $La_3Ni_2O_7$).

• Pairing Symmetries:
  The study analyzes various gap structures:

  • Isotropic $s$-wave ($s$).
  • $d$-wave pairing (both intra-orbital on nearest-neighbor bonds and inter-orbital).
  • $s_{\pm}$-wave pairing.

• Calculation Approaches:

  1. Full Multiorbital (MO) Calculation: The Raman response is computed using the full Green's function in the band basis, explicitly including inter-band and intra-band scattering processes and the mixing of orbital characters via the rotation matrix $U$. The Raman vertices are derived from the second derivatives of the full Hamiltonian.
  2. Additive (Isolated Band - IB) Approximation: The response is calculated as the sum of independent contributions from each band, ignoring inter-band coherence and mixing. This is a common approximation used in cuprates.
  3. Screening: For $A_{1g}$ symmetry, Coulomb screening is included to account for charge conservation constraints.

3. Key Contributions

• Validation of Approximations: The paper rigorously compares the Full MO approach with the IB approximation. It demonstrates that while the IB approximation works well for simple bilayer one-orbital models (like cuprates), it fails qualitatively in multiorbital systems where inter-orbital hopping and orbital mixing are significant.
• Fingerprint Identification: The authors identify specific "fingerprints" in the Raman spectra (peak positions, intensities, and low-energy power laws) that distinguish between different pairing symmetries and orbital models.
• Power Law Analysis: They confirm that for $d$-wave pairing, the low-energy Raman response follows characteristic power laws ($\sim \omega$ for $A_{1g}$ and $\sim \omega^3$ for $B_{1g}$) even in multiorbital systems, provided the analysis is done correctly.
• Screening Mechanism: The study elucidates how interlayer coupling ($t_{\perp}$) affects the efficiency of Coulomb screening in the $A_{1g}$ channel, showing that strong coupling reduces screening efficiency, making the signal more observable.

4. Key Results

A. Bilayer One-Orbital Model

• Weak Coupling ($t_{\perp}/t \approx 0.05$): Resembles infinite-layer nickelates. The two Fermi surface (FS) sheets are close, leading to overlapping peaks at $2\Delta_0$. Screening in the $A_{1g}$ channel is highly efficient (signal nearly vanishes), while $B_{1g}$ remains strong.
• Strong Coupling ($t_{\perp}/t \approx 0.85$): Resembles bilayer $La_3Ni_2O_7$. The FS sheets are well-separated. Two distinct peaks appear for $d$-wave pairing. Screening in $A_{1g}$ is much less efficient, allowing a measurable signal.

B. Single-Layer Two-Orbital Model

• Discrepancy in Approximations: For intra-orbital $d$-wave pairing on nearest-neighbor bonds, the MO calculation predicts two distinct peaks (low energy $\sim \Delta_0/2$ and high energy $\sim 1.5\Delta_0$). The IB approximation suppresses the high-energy peak and broadens the low-energy one, leading to qualitatively different predictions.
• Inter-orbital $d$-wave: The MO method shows that the low-energy peak shifts to higher energies compared to the intra-orbital case. The IB approximation fails to capture the correct spectral weight distribution, particularly suppressing the $B_{1g}$ response significantly compared to MO.
• Conclusion: The IB approximation is unreliable for multiorbital nickelates; full MO calculations are necessary.

C. Bilayer Two-Orbital Model (Most Relevant for $La_3Ni_2O_7$)

• $s$-wave Pairing: Both $A_{1g}$ and $B_{1g}$ show sharp peaks at $2\Delta_0$. Distinguishing between intra-layer and inter-layer $s$-wave is difficult based solely on peak position.
• $d$-wave Pairing:
  • Intra-layer: Similar to the single-layer case, two peaks are observed. The second peak (higher energy) is more intense than the first.
  • Inter-orbital: A sharp peak appears at $\sim 2\Delta_0$ in both channels. Crucially, the $A_{1g}$ channel shows additional low-energy spectral weight (peaks A and B) that is absent in the $B_{1g}$ channel due to the momentum dependence of the Raman vertex.
• $s_{\pm}$-wave: Shows a dominant peak at $2\Delta_0$ and sub-leading peaks at lower energies in $A_{1g}$, consistent with more complex DFT-based models found in recent literature.
• IB vs. MO: In the bilayer two-orbital case with strong interlayer coupling, the IB approximation reproduces the MO results better than in the single-layer case, but discrepancies in peak intensity and low-energy behavior remain.

5. Significance

• Experimental Guidance: The paper provides a roadmap for experimentalists analyzing Raman data in nickelates. It warns that observing a peak at a specific energy does not uniquely determine the pairing symmetry without considering the full multiorbital nature of the material.
• Resolution of Controversies: The results suggest that if Raman experiments show a strong $A_{1g}$ signal with low-energy features, it could point towards specific inter-orbital pairing mechanisms or strong interlayer coupling, helping to resolve conflicting reports on gap symmetry in $La_3Ni_2O_7$.
• General Applicability: The formalism and findings extend beyond nickelates to other multiorbital superconductors (e.g., iron-based superconductors), emphasizing that the "additive" approach is often insufficient for systems with strong orbital mixing.
• Future Directions: The authors suggest that the large gap and multiband nature of nickelates make them ideal platforms for studying collective modes (Higgs modes) in optical responses.

In summary, this work establishes that full multiorbital calculations are essential for interpreting Raman scattering in nickelates, as simplified additive models can lead to incorrect conclusions regarding gap symmetry and magnitude.