d-wave superconductivity in quintuple-layer bulk superconductor La 6 Ni 5 O 12

arXiv:2610.22212 · cond-mat.supr-con, cond-mat.str-el · Submitted 2026-10-15 · Read on arXiv

Listen

Transcript

Introduction to the show: ident: Quantum Radio. Generated commentary on the latest quantum physics and condensed matter papers.

Kai: I'm Kai.

Mira: And I'm Mira, here with Lev.

Kai: Today's paper: "d-wave superconductivity in quintuple-layer bulk superconductor La 6 Ni 5 O 12".

Mira: The gist: Using an effective five-layer Hubbard model and variational Monte Carlo calculations, we find that a spin-singlet dx2−y2-wave state is energetically favored in La6Ni5O12,

Kai: First, who's behind it and why it matters.

Title and authors: Kai: We’ve just looked at the core findings of "d-wave superconductivity in quintuple-layer bulk superconductor La six Ni five O twelve" where they found a favored spin-singlet d x two-y squared wave state #pg1.

Mira: That preference for d-wave over extended s-wave is what drives the next part of the discussion, showing how this result fits into the broader picture of nickelate pairing symmetry #pg1.

Lev: Before we go deeper, I want to talk about what they actually built and measured to get those numbers #pg3.

Kai: Right, Lev, can you tell us a bit about the setup for these VMC calculations? What was the actual hardware used for this modeling?

Lev: They used an effective five-layer Hubbard model derived from DFT calculations using the projector augmented-wave method and generalized gradient approximation in the PBE parametrization #pg3.

Mira: That DFT framework gives them a tight-binding parameter set, which they then use to define the mean-field Hamiltonian H mf #pg3.

Kai: And that Hamiltonian has specific parameters like interlayer hopping amplitudes e', layer-dependent chemical potentials e mu, and intralayer pairing amplitudes eta s,d #pg3.

Mira: They then tested three different mean-field ansatzes: the d-wave state, the extended s-wave state, and an unpaired reference state #pg3.

Lev: The variational parameters characterizing that mean-field Hamiltonian are all real numbers, and they don't explicitly include an interlayer pairing field #pg3.

Kai: So the calculation is self-contained in terms of those variational parameters, which simplifies things for the model setup #pg1.

Mira: That simplicity is important because it keeps the focus squarely on how those specific pairing amplitudes eta s,d behave as we change layer structure #pg1.

Lev: If you were to run this on real hardware, you'd have to worry about the overhead of calculating those complex layer-resolved hopping terms t' (k x - k y) squared #pg3 <ref:2610.22212#pg1>.

Kai: The tight-binding parameters themselves were derived from DFT calculations using the PBE parametrization and a plane-wave kinetic-energy cutoff of five hundred eV #pg3.

Mira: That means the input physics for the model is based on standard Density Functional Theory calculations, specifically within VASP #pg3.

Lev: And they sampled their Brillouin zone with a seven times seven times seven k-point mesh for those DFT calculations #pg3 <ref:2610.22212#pg2>.

Kai: So to recap, they constructed this model using five layers and used VMC to find the favored spin-singlet d x two-y squared state #pg1 <ref:2610.22212#pg1>.

Mira: It sets up the baseline for what we expect to see in experimental probes like the density of states and superfluid density #pg1.

Lev: It’s a good starting point, but remember that the actual computational cost of running this level of detail on a real system is substantial #pg3.

The paper's summary: Kai: Now let’s look at what the paper actually summarizes regarding these results for "d-wave superconductivity in quintuple-layer bulk superconductor La six Ni five O twelve" #pg1.

Mira: The core summary is that the d-wave superconducting state is energetically favored over any non-superconducting ansatz they tested #pg3.

Lev: And crucially, the extended s-wave pairing is suppressed, meaning it's roughly two orders of magnitude weaker than the d-wave amplitudes #pg1.

Kai: That suppression isn't just a small correction; it’s a big difference in the energy landscape of the system #pg1.

Mira: They also stress that this finding is consistent with experimental signatures we see in bulk experiments, like the V-shaped density of states spectrum #pg1.

Lev: So they are linking the theoretical ground state directly to measurable spectroscopic features #pg3.

Kai: And they point out specific predictions for observables: the superfluid density shows a linear low-temperature depletion, and the spin-lattice relaxation rate follows one/T 1T proportional to T squared #pg3 <ref:2610.22212#pg1>.

Mira: These are all consistent with what we expect from a clean nodal superconductor, which is what this material seems to be behaving like #pg1.

Lev: The paper also explicitly states that the optimized extended s-wave pairing amplitude in the trilayer four hundred thirty-eight compound is vanishingly small #pg2 <ref:2610.22212#pg2>.

Kai: So they are using the bulk six thousand five hundred twelve as a key platform to test these pairing tendencies across different layer numbers #pg1 <ref:2610.22212#pg2>.

Mira: This is important because it helps establish a common d x two-y squared pairing tendency across the reduced RP nickelate family #pg1 <ref:2610.22212#pg1>.

Lev: The methodology itself, using the Gutzwiller projector on the mean-field ground state mf, is what allows them to capture that strong correlation physics #pg3.

Kai: It’s a solid way to handle the strong electron correlations inherent in these nickelates #pg1.

Mira: So, in short, the summary says: d-wave wins, s-wave loses big time, and the experimental signatures match a clean nodal superconductor profile #pg1.

Lev: It’s a very clear conclusion from their calculation that they are using standard correlated electron techniques to pinpoint the symmetry #pg3.

The paper's improvements: Kai: Now let’s talk about what improvements the authors suggest or what this work adds to the existing understanding of "d-wave superconductivity in quintuple-layer bulk superconductor La six Ni five O twelve" #pg1.

Mira: The main improvement is using this specific material, La six Ni five O twelve as a benchmark for probing the pairing symmetry of the reduced RP nickelate family #pg2.

Lev: It’s valuable because it connects superconductivity across different layer numbers, linking the infinite-layer LaNiO2 to finite-layer compounds like four hundred thirty-eight #pg2 <ref:2610.22212#pg2>.

Kai: They are using this quintuple-layer system to bridge that gap, showing how the pairing tendency evolves as you change the number of consecutive NiO2 layers #pg2.

Mira: They establish a robust d x two-y squared pairing tendency across the reduced RP nickelate family, suggesting it’s not just an artifact of one specific layer count #pg1 <ref:2610.22212#pg1>.

Lev: And they suggest that the strength of superconductivity is controlled by carrier concentration and layer structure #pg1.

Kai: That means if you want to understand why some nickelates are superconducting and others aren't, this paper gives you a way to look at those structural differences #pg2.

Mira: They also show that reducing the hole concentration by electron doping might be a promising route to stabilize d-wave superconductivity in four hundred thirty-eight #four hundred thirty-eight #pg1 <ref:2610.22212#pg2>.

Lev: That suggests we should investigate doping strategies to see if we can drive four hundred thirty-eight toward a more superconducting regime rather than just sitting near the boundary #pg2 <ref:2610.22212#pg2>.

Kai: It’s a practical direction for materials discovery, moving beyond just observing the pairing symmetry to influencing it #pg1.

Mira: So the improvement isn't just confirming d-wave in six thousand five hundred twelve; it’s creating a roadmap for understanding the whole nickelate family #pg2 <ref:2610.22212#pg2>.

Lev: It’s a useful piece of theoretical guidance for anyone trying to design correlated electron systems #pg3.

Conclusion: Kai: So we've walked through the title and authors, looked at the core summary of "d-wave superconductivity in quintuple-layer bulk superconductor La six Ni five O twelve" #pg1.

Mira: We’ve established that the VMC calculations strongly favor a spin-singlet d x two-y squared state and suppress the extended s-wave pairing significantly #pg3 <ref:2610.22212#pg1>.

Lev: And they’ve shown how this translates into measurable signatures like a V-shaped density of states and linear low-temperature depletion in the superfluid density #pg3.

Kai: This is important because it gives us concrete, testable predictions for bulk experiments, which is what we need to know #pg1.

Mira: The larger implication is that this work points toward a common d x two-y squared pairing tendency across the reduced RP nickelate family #pg1 <ref:2610.22212#pg1>.

Lev: It gives us a new reference point for how structural changes affect superconductivity in these materials #pg2.

Kai: So, to wrap up on "d-wave superconductivity in quintuple-layer bulk superconductor La six Ni five O twelve" it solidifies the picture of d-wave pairing being dominant here #pg1.

Mira: It's a foundation for exploring how carrier concentration and layer structure control the strength of that superconductivity #pg1.

Lev: For me, it confirms that theoretical models based on Gutzwiller projection can actually yield physically sensible results in this regime #pg3.

Kai: That’s it for this paper, a solid piece of work to keep us thinking about nickelates #pg1.

Zhan Wang, *, Yuxin Wang, *, Hui-Ke Jin, †Kμn Jiang, Jiangping Hu

Beijing National Laboratory for Condensed Matter Physics and Institute of Physics, Chinese Academy of Sciences, Beijing 100190, China · Kavli Institute for Theoretical Sciences, University of Chinese Academy of Sciences, Beijing, 100190, China · School of Physical Science and Technology, ShanghaiTech University · School of Physical Sciences, University of Chinese Academy of Sciences

cond-mat.supr-con, cond-mat.str-el

Submitted: 2026-10-15

Updated: 2026-10-15

Comments: 6+4 pages, 3+1 figures

License: http://creativecommons.org/licenses/by/4.0/

Importance score: 85/100

The gist: The gist: Using an effective five-layer Hubbard model and variational Monte Carlo calculations, we find that a spin-singlet dx2−y2-wave state is energetically favored in La6Ni5O12, while

Key concepts

Effective Five-Layer Hubbard Model: This is a simplified mathematical model used to describe the electronic behavior of the quintuple-layer nickelate. It focuses on the interactions between electrons in five consecutive NiO2 planes, treating them as interacting layers. This model helps researchers understand how electron correlations influence superconductivity in this material.
d_{x^2-y^2}-wave State: This describes a specific pattern of electron pairing where the superconducting pairs have a symmetry matching the d_{x^2-y^2} orbital shape. This type of pairing is characterized by nodes (points where the superconducting gap closes) in momentum space, which is a key feature for nodal superconductors.
Variational Monte Carlo (VMC): VMC is a computational technique used to find the lowest energy state of a complex quantum system, like La6Ni5O12. It uses statistical sampling to explore different possible electron configurations and determines the most stable superconducting state by minimizing the energy functional.
Clean Nodal Superconductor: This term describes a superconductor where the superconducting gap has nodes (zero values) in momentum space, rather than being fully gapped everywhere. The calculated experimental signatures—like a linear depletion of superfluid density and a specific temperature dependence for spin-lattice relaxation—confirm that La6Ni5O12 behaves as such a clean nodal superconductor.

Summary

The gist: Using an effective five-layer Hubbard model and variational Monte Carlo calculations, we find that a spin-singlet dx2−y2-wave state is energetically favored in La6Ni5O12, while extended-s-wave pairing is strongly suppressed, providing consistent signatures of a clean nodal superconductor through various experimental probes.

Modeling the System

The study utilizes an effective five-layer Hubbard model for the quintuple-layer nickelate La6Ni5O12 and employs variational Monte Carlo (VMC) calculations to determine the ground state Using an effective five-layer Hubbard model and variational Monte Carlo calculations, we find that a spin-singlet dx2−y2-wave state is energetically favored, while extended-s-wave pairing is strongly suppressed. The crystal structure of the quintuple-layer 6512 compound consists of five consecutive NiO2 planes formed by topotactic reduction

Superconducting State Determination

The VMC calculations investigated the ground state of the multilayer Hubbard model using a Gutzwiller-projected wave function Ψ⟩ = PGΦmf⟩ We find that the d-wave superconducting state is energetically favored over the non-superconducting ansatz, and the extended s-wave pairing is strongly suppressed, i.e., being roughly two orders of magnitude weaker than the d-wave amplitudes

Experimental Signatures

Based on this correlated superconducting state, several experimentally accessible observables were calculated to provide direct predictions for bulk experiments The calculated quasiparticle density of states exhibits a characteristic V-shaped low-energy spectrum, reflecting the linearly vanishing density of states associated with the d-wave nodes The superfluid density shows a linear low-temperature depletion, and the Knight shift decreases below Tc The spin-lattice relaxation rate follows 1/T1T ∝ T2, providing consistent signatures of a clean nodal superconductor

Comparison with Other Nickelates

The framework was extended to related reduced RP nickelates to compare pairing tendencies across the family For the trilayer 438 compound, the system lies close to the boundary between a d-wave superconducting state and a non-superconducting state, whereas the optimized extended-s-wave pairing amplitude is vanishingly small For the infinite-layer 112 compound, VMC calculations find the lowest-energy d-wave state below the best non-superconducting reference states, while the extended s-wave solution again exhibits a negligible pairing amplitude

Conclusion on Pairing Tendency

The results identify bulk La6Ni5O12 as a promising platform for determining the pairing symmetry of nickelate superconductors and suggest a common dx2−y2-wave pairing tendency across the reduced RP nickelate family Taken together, these results point to a common dx2−y2-wave pairing tendency across the reduced RP nickelate family, with the strength of superconductivity controlled by carrier concentration and layer structure The bulk 6512 therefore provides not only a new platform for establishing the pairing symmetry of nickelate superconductors, but also a bridge between superconducting infinite-layer 112 and the emerging superconducting regime of finite-layer compounds such as 438

How it works

The study constructs a two-dimensional five-band tight-binding (TB) model for the Ni dx2−y2 sector through Wannierization, treating the La-derived band as a self-doping reservoir that modifies the filling of the Ni dx2−y2 bands The full Hamiltonian is given by H = Ht + U X il nil↑nil↓

Superfluid Density Calculation

The in-plane superfluid stiffness ρs(T) is calculated from the static transverse electromagnetic response of the renormalized effective mean field Hamiltonian, which takes a BdG form The normalized low-temperature depletion [ρs(0) − ρs(T)]/ρs(0) is linear in T, arising from thermally excited quasiparticles near the gap nodes

NMR Signatures

NMR provides complementary information on the spin structure of the Cooper pairs by calculating both the spin Knight shift and the nuclear spin-lattice relaxation rate using a renormalized Hamiltonian The spin contribution to the Knight shift Ks decreases below Tc and follows an approximately linear temperature dependence at low T, consistent with spinsinglet pairing At low temperatures, relaxation is dominated by quasiparticles near the d-wave nodes, yielding 1/T1T ∝ T2

Model Details

The tight-binding parameters for the five-layer model were derived from DFT calculations using the projector augmented-wave (PAW) method and generalized gradient approximation (GGA) in the PBE parametrization The variational parameters characterizing Hmf include interlayer hopping amplitudes eχll′ between adjacent layers, layerdependent chemical potentials eµl, and intralayer pairing amplitudes ηs,dl; see Supplemental Material S2 for details The Gutzwiller projector multiplies each configuration amplitude by exp − P il glnil↑nil↓

Tight-Binding Parameters

The tight-binding parameters of the five-layer model derived from DFT calculations, in units of eV, are summarized in Table I The intralayer dispersion is written as εl,k = εl + X 3 m=1 tmlγmk The interlayer hopping is parameterized as τll′ (k) = −2tll′⊥ (cos kx − cos ky) 2 The layer-resolved filling nl obtained from DFT-derived TB parameters is written on the right side

Trilayer and Infinite-Layer Comparison

For the trilayer 438 compound, VMC calculations are performed on L × L × 3 lattices with periodic boundary conditions in the in-plane directions at the nominal hole concentration x = 1/3 For the infinite-layer 112 compound, VMC calculations find the lowest-energy d-wave state below the best non-superconducting reference states The results support a robust dx2−y2-wave pairing tendency in 112, consistent with that obtained for the quintuple-layer 6512 compound

Supplemental Information

The DFT computational method was performed within the framework of density functional theory (DFT), as implemented in the Vienna ab initio Simulation Package (VASP) The exchange–correlation energy was treated using the generalized gradient approximation (GGA) in the Perdew–Burke–Ernzerhof (PBE) parametrization The plane-wave kinetic-energy cutoff was set to 500 eV, and the Brillouin zone was sampled using a 7 × 7 × 7 k-point mesh The convergence criterion for the electronic self-consistent calculations was set to 10−8 eV The mean-field ansatz Φmf⟩ entering Eq. (5) is the ground state of the free-fermion Hamiltonian Hmf = X k,σ X 5 l=1 h ε0l(k) − eµl i c†klσcklσ + X k,σ X 4 l=1 eχl,l+1(k) c†klσck,l+1,σ + h.c.i The nearest-neighbor spin-singlet pairing fields are e∆dl(k) = 2ηdl(cos kx − cos ky), and the extended s-wave state is e∆sl(k) = 2ηsl(cos kx + cos ky) All variational parameters are real, and no explicit interlayer pairing field is included The layer-dependent factor is zl = √nlhnl↑ + √nl↓nld√ (1 − nl/2)(nl/2) The total carrier concentration is controlled by fixing the number of doped holes Nh The momentum integration is performed on a 2000 × 2000 mesh with a Lorentzian broadening Γ = 0.

Improvements for AI systems

A robust predictive model for pairing symmetry in nickelate superconductors can be developed by using variational Monte Carlo (VMC) calculations for a realistic five-layer model of 6512 to determine that a spin-singlet dx 2−y squared-wave state is energetically favored.
The improved system can accurately predict experimental signatures from bulk measurements, as the VMC framework allows calculation of observables such as the tunneling density of states exhibits a characteristic V-shaped low-energy spectrum, the superfluid density shows a linear low-temperature depletion, and the spin-lattice relaxation rate follows 1/T1T ∝ T squared.
The AI system can perform comparative analysis across different material families by using the results to determine how pairing tendencies differ across nickelate compounds with different layer numbers, specifically identifying a robust d-wave tendency in infinite-layer LaNiO2 versus the behavior of trilayer La4Ni3O8, which lies close to the boundary between d-wave superconductivity and a non-superconducting state.
The system can guide materials discovery by suggesting specific doping strategies, as the analysis indicates that reducing the hole concentration by electron doping may provide a promising route to stabilize d-wave superconductivity in 438.

Abstract

The recent realization of bulk superconductivity in the quintuple-layer nickelate La 6 Ni 5 O 12 opens new opportunities to probe the pairing symmetry of the reduced Ruddlesden--Popper (RP) nickelate family. Using an effective five-layer Hubbard model and variational Monte Carlo calculations, we find that a spin-singlet d x 2-y squared-wave state is energetically favored, while extended- s-wave pairing is strongly suppressed. Based on the optimized superconducting state, we further calculate several experimentally accessible observables. The tunneling density of states exhibits a characteristic V-shaped low-energy spectrum, the superfluid density shows a linear low-temperature depletion, the Knight shift decreases below T c, and the spin-lattice relaxation rate follows 1/T 1T proportional to T squared, together providing consistent signatures of a clean nodal superconductor. Extending the same framework to related reduced RP nickelates, we find a robust d-wave tendency in infinite-layer LaNiO 2, while trilayer La 4 Ni 3 O 8 at its nominal filling lies close to the boundary between d-wave superconductivity and a non-superconducting state. These results establish bulk La 6 Ni 5 O 12 as a promising platform for determining the pairing symmetry of nickelate superconductors and suggest a common d x 2-y squared-wave pairing tendency across the reduced RP nickelate family.