Co-operating multiorbital and nonlocal correlations in bilayer nickelate

This study employs an advanced many-body framework to reveal how the interplay between multiorbital physics and nonlocal correlations in bilayer nickelate La$_3$Ni$_2$O$_7$ determines the position of the flat $\gamma$ band and drives spin-polaron formation, offering a potential explanation for conflicting angle-resolved photoemission experimental results.

The Gist

Imagine a superconductor as a busy highway where electrons are cars. Usually, these cars drive smoothly, but in certain materials, they get stuck in traffic jams or form strange, slow-moving convoys. The material in this paper, La₃Ni₂O₇ (a bilayer nickelate), is a very special highway that scientists are trying to understand because it might hold the key to superconductivity (electricity with zero resistance) at higher temperatures.

Here is the story of what the researchers found, explained without the heavy math.

The Setup: A Three-Lane Highway

Think of the electrons in this material as living in a three-story apartment building (the "three-orbital model").

1. Two busy floors (Orbitals $\alpha$ and $\beta$): These are like the main living areas where most of the action happens. The electrons here are active and moving around.
2. One quiet, flat floor (Orbital $\gamma$): This is the "flat band." Imagine a perfectly flat, featureless meadow. In physics, when electrons are on a flat band, they don't move much; they just sit there, heavy and slow. This flat floor is the star of the show.

The Problem: The "Flat Floor" is in the Wrong Place

The researchers wanted to know: Is this flat meadow ($\gamma$ band) above the ground (Fermi level) or below it?

• Below the ground: The electrons are stuck in the basement. They are occupied but not really participating in the high-energy party.
• Above the ground: The electrons are free to roam the main floor, ready to interact with everything else.

Recent experiments (like taking photos of the electrons with a camera called ARPES) gave conflicting answers. Some said the flat floor was below ground; others said it was above. The scientists in this paper asked: Why are the photos different?

The Solution: A "Social Distance" Knob

The researchers built a super-computer simulation to act as a control room. They found that the position of that flat floor depends on a "knob" they call the interorbital interaction ($J$). You can think of this knob as a social distance regulator between the electrons in different apartments.

• Turning the knob down (Low $J$): The electrons in the flat floor stay put, below the ground level. The system is calm. The electrons behave like normal, predictable cars. The "photos" taken here look like the flat floor is hidden underground.
• Turning the knob up (High $J$): The electrons get more agitated. The flat floor gets pushed up until it crosses the ground level. Now, the electrons are free to roam the main floor.

The Big Surprise: The "Shadow Ghost"

Here is where it gets magical. When the researchers turned the knob high enough to push the flat floor above the ground, something weird happened. The flat floor didn't just stay flat; it split in two.

1. The Real Band: One part stayed where it was, moving normally.
2. The Shadow Band: A "ghost" version of the band appeared below the ground level.

The Analogy: Imagine a loudspeaker playing a song. Suddenly, you hear the main song, but you also hear a faint, distorted echo coming from the basement. That echo is the Shadow Band.

Why Does This Happen? (The Spin-Polaron)

Why did the ghost appear? It's because of magnetic fluctuations.

• When the flat floor crosses the ground, the electrons there start interacting with "magnetic waves" (called paramagnons) that are rippling through the material.
• Imagine a electron trying to run through a crowd of people who are all waving their arms (magnetic waves). The electron gets tangled up with the crowd, slowing down and forming a heavy, sluggish bundle.
• In physics, this bundle is called a Spin-Polaron. It's a "Frankenstein" creature made of an electron and a magnetic wave stuck together.
• Because this creature is heavy and confused, it creates that "shadow" signal below the ground level.

The "Aha!" Moment

The researchers realized that the conflicting photos from the real world weren't mistakes. They were just snapshots taken at different settings of the "social distance knob."

• If the knob was set low, the flat floor was underground (no shadow).
• If the knob was set high, the flat floor was above ground, creating a shadow below.

The Takeaway

This paper solves a mystery. It tells us that the bilayer nickelate is a shape-shifter. Depending on how the electrons interact with each other, the material can look completely different.

• It explains why some experiments see the flat band and others don't.
• It predicts that if we can control this "knob" (perhaps by changing pressure or doping the material), we might be able to create these "shadow bands" on purpose.
• These shadow bands (spin-polarons) might be the secret ingredient that helps the material become a superconductor.

In short: The material is like a stage with a hidden trapdoor. Depending on how you tune the lights (the interactions), the actors (electrons) either stay in the basement, come up to the stage, or create a spooky ghost reflection below the stage. Understanding this trick helps us figure out how to make better superconductors.

Technical Summary

1. Problem Statement

The paper addresses the complex electronic structure of the high-pressure superconducting bilayer nickelate La$_3$Ni$_2$O$_7$. While Density Functional Theory (DFT) reveals a low-energy electronic structure involving three active orbitals (two $d_{x^2-y^2}$-dominated antibonding bands, one $d_{z^2}$-dominated antibonding band, and a non-bonding flat $d_{z^2}$-dominated $\gamma$ band), the role of nonlocal electronic correlations remains unclear.

Previous studies have largely relied on local approximations (like Dynamical Mean-Field Theory, DMFT) or reduced orbital models. However, the interplay between multiorbital physics and nonlocal self-energy effects is critical for understanding:

• The precise location of the flat $\gamma$ band relative to the Fermi level ($E_F$).
• The nature of low-energy excitations and potential competing orders (e.g., charge ordering vs. superconductivity).
• The controversy in recent Angle-Resolved Photoemission Spectroscopy (ARPES) experiments, which report conflicting results regarding whether the $\gamma$ pocket lies above or below $E_F$.

2. Methodology

The authors employ an advanced many-body framework to solve an effective three-orbital Hubbard model derived from DFT.

• Model Construction:

  • The system is mapped onto maximally localized Wannier orbitals: $w_1$ (derived from non-bonding $Ni-d_{z^2}$ across the bilayer, representing the $\gamma$ band) and $w_2, w_3$ (derived from $Ni-d_{x^2-y^2}$ hybridized with $O-2p$, representing the $\alpha$ and $\beta$ bands).
  • The Hamiltonian includes intra-orbital Coulomb repulsion ($U$), inter-orbital repulsion ($U'$), and a specific Hund's coupling term ($J$) that accounts for competition between intra-layer Hund's coupling and inter-layer antiferromagnetic coupling.
  • The model is filled with three electrons per unit cell.

• Computational Approach: D-TRILEX

  • The model is solved using the Dual Triply Irreducible Local Expansion (D-TRILEX) method.
  • Key Advantage: Unlike standard DMFT, which treats local correlations exactly but neglects nonlocal spatial fluctuations, D-TRILEX performs a diagrammatic expansion beyond DMFT. It consistently incorporates momentum-dependent ($k$-dependent) collective fluctuations in charge, spin, and orbital channels into the electronic self-energy.
  • The impurity problem is solved using Continuous-Time Quantum Monte Carlo (CT-QMC) within the w2dynamics package.
  • Calculations are performed on a $36 \times 36$ $k$-grid in the Brillouin zone, with analytical continuation via the Maximum Entropy Method.

3. Key Contributions

• Orbital-Dependent Fluctuation Mechanism: The study demonstrates that the nature of electronic fluctuations in La$_3$Ni$_2$O$_7$ is highly sensitive to the occupation of the $w_1$ orbital (controlled by the interorbital interaction parameter $J$).
• Spin-Polaron Formation: The authors identify a mechanism where itinerant electrons scatter off ferromagnetic (FM) paramagnon excitations originating from the flat $\gamma$ band, leading to the formation of spin-polaron bound states.
• Resolution of ARPES Controversy: The paper provides a theoretical explanation for conflicting experimental observations regarding the $\gamma$ band's position, attributing the differences to the formation of a "shadow band" (incoherent spectral weight) when the flat band crosses the Fermi level.

4. Results

A. Regime 1: Flat $\gamma$ Band Below $E_F$ ($J = 0.10$ eV)

• Band Position: The flat part of the $\gamma$ band lies just below the Fermi level.
• Fluctuations: Non-local charge and spin fluctuations are weak. The leading fluctuations originate from the $w_2$ and $w_3$ orbitals ($\alpha$ and $\beta$ sheets).
• Self-Energy: The electronic self-energy is nearly momentum-independent, resembling standard DMFT results. The system behaves as a conventional metal with weak correlations.
• Susceptibility: Charge fluctuations are dominated by small momenta (near $\Gamma$), while spin fluctuations are incommensurate.

B. Regime 2: Flat $\gamma$ Band Crossing/Above $E_F$ ($J \geq 0.13$ eV)

• Band Position: Increasing $J$ reduces the occupancy of the $w_1$ orbital, shifting the flat $\gamma$ band upward across $E_F$.
• Critical Instability: At $J \approx 0.12$ eV, the system approaches a Charge Ordering (CO) instability.
• Spin-Polaron Formation: Once the flat band crosses $E_F$ (e.g., at $J = 0.15$ eV):
  • Magnetic Fluctuations: Strong, low-momentum ferromagnetic (FM) spin fluctuations emerge, driven by the large spectral weight of the $w_1$ orbital at $E_F$.
  • Self-Energy: The self-energy becomes strongly momentum-dependent, exhibiting non-Fermi-liquid behavior specifically at the $M$ point (where the $\gamma$ band is located).
  • Shadow Band: The scattering of electrons with these FM paramagnons creates a bound state (spin polaron). This manifests in the spectral function as a splitting of the flat band: a main peak above $E_F$ and a weaker, incoherent "shadow band" below $E_F$.
• Spectral Weight: The incoherent spectral weight of the shadow band diminishes as $J$ increases further, eventually vanishing as the band moves deeper into the unoccupied region.

C. Temperature Dependence

• The formation of the spin-polaron band is robust. While lowering the temperature enhances the spectral weight of the main peak (driving the system toward CO instability below ~150 K), the position and intensity of the spin-polaron shadow band remain largely unchanged in the normal state ($T > 150$ K).

5. Significance and Implications

• Theoretical Insight: The work establishes that nonlocal correlations are not merely corrections but are fundamental to the low-energy physics of bilayer nickelates. It highlights a regime where multiorbital physics and nonlocal self-energy effects cooperate to generate exotic composite excitations (spin polarons).
• Experimental Validation:
  • ARPES: The predicted "shadow band" below the Fermi level offers a direct explanation for recent ARPES data that observed spectral weight in the $\gamma$ pocket below $E_F$ even when the main band appeared above it.
  • Proposed Tests: The authors suggest experimental verification via:
    1. Measuring a strong temperature-dependent increase in uniform spin susceptibility (indicating FM tendencies).
    2. Optical measurements to detect composite excitations.
    3. Observing the effect of hole doping, which should trigger the polaronic band formation once the $\gamma$ pocket crosses $E_F$.
• Broader Context: The findings challenge the view that nonlocal correlations are negligible in multiorbital systems (unlike in some iron-based superconductors) and suggest that the interplay of flat bands and nonlocal fluctuations is a key ingredient in the phase diagram of high-$T_c$ nickelates.

In summary, the paper utilizes the D-TRILEX method to reveal that the interplay between multiorbital interactions and nonlocal spin fluctuations in La$_3$Ni$_2$O$_7$ leads to the formation of spin-polaron bound states, providing a unified theoretical framework to explain conflicting experimental observations of the electronic structure.