Electron Doping of $\mathrm{La_3Ni_2O_7}$ Thin Films: Candidate Metal Dopants and Their Potential Impact on Superconductivity

This study employs first-principles calculations to identify zirconium, hafnium, and thorium as effective electron dopants for $\mathrm{La_3Ni_2O_7}$ thin films that enhance interlayer hopping and potentially elevate superconducting $T_c$, while ruling out cerium as a viable candidate.

The Gist

Imagine a superconductor as a busy dance floor where electrons pair up and glide across the room without bumping into anything (resistance). For decades, scientists have been obsessed with a specific type of dance floor made of copper and oxygen, called cuprates. They figured out that if you add extra "holes" (missing dancers) to the floor, the music gets better, and the dancing becomes super-efficient at high temperatures.

Recently, a new dance floor made of nickel and oxygen, called La₃Ni₂O₇, was discovered. It's like a cousin to the copper dance floor, but it has a secret: it can superconduct at even higher temperatures (over 80 Kelvin). However, scientists are still arguing about why it works. Is it because of a specific type of dancer (an orbital called $\gamma$ or $d_{z^2}$) that needs to be on the floor? Or does the dance work even if that specific dancer leaves?

To solve this mystery, the researchers in this paper decided to try a different trick: instead of removing dancers (hole doping), they tried to add extra dancers (electron doping).

Here is what they found, explained simply:

1. The "Wrong Key" Attempt: Cerium (Ce)

In the old copper dance floors, scientists used Cerium to add extra electrons. It worked like a charm there. So, the researchers thought, "Let's try Cerium on this new nickel floor!"

The Result: It failed.
Think of Cerium as a guest who shows up to the party but decides to sit in the corner and not dance. Even when they added a lot of Cerium, the low-energy dance floor (where the superconducting magic happens) looked exactly the same as before. The extra electrons didn't make it onto the main floor; they got stuck in the VIP lounge (high-energy states). The nickel floor simply didn't accept Cerium as an electron donor.

2. The "Right Keys": Zirconium, Hafnium, and Thorium

Since Cerium didn't work, the team tried other guests: Zirconium (Zr), Hafnium (Hf), and Thorium (Th).

The Result: Success!
These three elements acted like enthusiastic guests who immediately jumped onto the dance floor. They successfully added extra electrons to the low-energy bands.

• Thorium was the most energetic, pushing the dancers down to a lower energy level, effectively "filling up" the floor with new electrons.
• Zirconium and Hafnium also worked well, though they behaved slightly differently than Thorium.

3. How the Dance Floor Changed (The Physics)

When these new guests arrived, they didn't just add numbers; they changed the shape of the room.

• The "Bridge" Strengthened: The nickel floor has two layers of dancers. For the superconducting magic to happen, the dancers on the top layer need to talk to the dancers on the bottom layer. The researchers found that adding Zr, Hf, or Th built a stronger "bridge" (called interlayer hopping) between these layers.
• The Connection: This stronger bridge means the dancers are more tightly coupled. In the world of superconductors, a stronger connection between layers often leads to a higher "temperature limit" for the superconducting state. It's like tightening the springs on a trampoline; the bounce becomes more powerful.

4. Why This Matters

The big debate in the scientific community is: Does the superconductivity depend on that specific $d_{z^2}$ dancer being present on the floor?

• Hole doping (removing dancers) hasn't been able to settle this argument yet.
• Electron doping (adding dancers) pushes that specific dancer off the main stage (below the energy level where the action happens).

By successfully adding electrons with Zr, Hf, and Th, the researchers have created a new way to test the theory. If the superconductivity disappears when these specific dancers are pushed off the stage, we know they were essential. If the dancing continues, we know the mechanism is different.

Summary

This paper is a "guest list" for a nickel-based superconductor.

• Cerium was invited but didn't show up to the dance (failed to dope).
• Zirconium, Hafnium, and Thorium showed up, brought extra energy, and strengthened the connection between the two layers of the material.
• This gives scientists a new tool to figure out the secret recipe for high-temperature superconductivity in nickel materials, potentially helping us understand how to make even better superconductors in the future.

The paper stops at identifying these candidates and explaining how they change the electronic structure. It does not claim to have built a working device or a commercial product yet; it's purely about understanding the fundamental rules of the dance.

Technical Summary

Technical Summary: Electron Doping of La3Ni2O7 Thin Films

Problem Statement
The bilayer Ruddlesden-Popper nickelate La$_3$Ni$_2$O$_7$ has emerged as a critical platform for investigating high-temperature superconductivity, sharing quasi-two-dimensional NiO$_2$ layers with cuprates but featuring a distinct electronic structure involving both Ni-3$d_{x^2-y^2}$ and Ni-3$d_{z^2}$ orbitals. A central theoretical debate concerns the role of the $\gamma$ band (derived from Ni-3$d_{z^2}$ orbitals) in the pairing mechanism: whether superconductivity requires the $\gamma$ band to cross the Fermi level or if it persists when the band lies below it. While hole doping via Strontium (Sr) substitution or oxygen reduction has been extensively studied, it remains unclear if these methods can definitively tune the $\gamma$ band position to resolve this debate. Consequently, electron doping remains a largely unexplored avenue for La$_3$Ni$_2$O$_7$, despite its success in cuprate systems (e.g., Nd$_{2-x}$Ce$_x$CuO$_4$).

Methodology
The authors employed first-principles density functional theory (DFT) calculations using the Vienna Ab initio Simulation Package (VASP). Key methodological features include:

• Models: Simulations of 1-unit-cell (1UC) thick La$_{3-x}$R$_x$Ni$_2$O$_7$ thin films, where R represents tetravalent dopants (Ce, Th, Zr, Hf) or trivalent Sr for comparison.
• Substrates: Calculations were performed on three representative substrates (LaAlO$_3$, NdGaO$_3$, SrTiO$_3$) to account for varying in-plane lattice constants (compressive to tensile strain).
• Electronic Structure: The DFT+U method was applied with an effective $U$ of 3.5 eV for Ni-3$d$ orbitals to account for localization. For Ce doping, $U$ was also applied to Ce-4$f$ orbitals.
• Analysis Tools:
  • Wannier90: Used to extract tight-binding parameters, specifically the interlayer hopping integral $t_\perp$ between $d_{z^2}$ orbitals.
  • Constrained Random Phase Approximation (cRPA): Used to calculate effective screened interaction parameters, including the Hubbard $U$ and Hund's coupling $J_H$.
  • Bader Charge Analysis: Employed to evaluate charge transfer and valence states.
  • Density of States (DOS): Integrated to determine electron filling numbers for specific orbitals ($d_{z^2}$, $d_{x^2-y^2}$, and associated O-2$p$ orbitals).

Key Contributions and Results

1. Identification of Effective Electron Dopants:
   The study systematically screened tetravalent elements. Contrary to expectations based on cuprate physics, Cerium (Ce) doping fails to effectively introduce electron carriers into the low-energy bands of La$_3$Ni$_2$O$_7$. Even at high concentrations ($x=1$), the low-energy bands remain nearly identical to the pristine system, suggesting Ce retains a Ce$^{3+}$ valence state rather than Ce$^{4+}$.
   In contrast, Zirconium (Zr), Hafnium (Hf), and Thorium (Th) are identified as efficient electron dopants. Substitution of La with these elements successfully shifts the $d_{z^2}$ and $d_{x^2-y^2}$ bands downward in energy, increasing electron occupancy in the active $e_g$ orbitals.

2. Structural and Electronic Modifications:

   • Bond Lengths: Electron doping generally elongates vertical Ni–O bonds (due to the expansion of the Ni$^{2+}$ electron cloud), whereas hole doping (Sr) shortens them. However, Zr and Hf doping induce a lattice contraction due to their smaller ionic radii compared to La, which competes with the electronic expansion.
   • Orbital Occupancy: DOS analysis confirms that Zr, Hf, and Th doping increases electron counts in both in-plane ($d_{x^2-y^2}$) and out-of-plane ($d_{z^2}$) orbitals. Notably, Th doping preferentially populates the $d_{z^2}$-derived bands, while Zr and Hf favor the $d_{x^2-y^2}$ bands.
   • Charge Transfer: Bader analysis shows that electron dopants supply additional charge to the system (approx. +0.57 $|e|$ per Th atom), while hole dopants reduce it.

3. Impact on Superconducting Parameters:

   • Interlayer Hopping ($t_\perp$): A crucial finding is that electron doping significantly increases the interlayer hopping integral $t_\perp$ between $d_{z^2}$ orbitals. The ratio $t_\perp/t_{1x}$ (interlayer vs. intralayer hopping) increases from $\sim1.7$ in pristine La$_3$Ni$_2$O$_7$ to $\sim2.5$ in electron-doped compounds. This contrasts with hole doping, which reduces this ratio.
   • Interaction Parameters: cRPA calculations reveal that while the average Hubbard $U$ decreases in electron-doped compounds, the relative strength of Hund's coupling ($J_H/U$) increases. This enhancement of the Hund's coupling regime may favor specific pairing mechanisms.

Significance and Claims
The paper claims to provide a viable route for achieving electron doping in La$_3$Ni$_2$O$_7$ thin films, filling a gap in the experimental and theoretical exploration of this system. By identifying Zr, Hf, and Th as effective dopants, the work offers a pathway to experimentally tune the position of the $\gamma$ band relative to the Fermi level. This capability is presented as essential for clarifying the ongoing debate regarding the pairing mechanism (e.g., distinguishing between scenarios where the $\gamma$ band crosses the Fermi level versus those where it does not). Furthermore, the predicted enhancement of interlayer coupling ($t_\perp$) and Hund's coupling suggests that electron-doped La$_3$Ni$_2$O$_7$ could exhibit distinct, and potentially elevated, superconducting properties compared to its hole-doped counterpart. The study serves as theoretical guidance for experimentalists aiming to synthesize these specific doped phases to probe the fundamental physics of nickelate superconductivity.