Origin of superconductivity in bilayer nickelates: a Quantum Monte Carlo study for a sign-problem-free effective model

arXiv:2512.09025 · cond-mat.str-el, cond-mat.supr-con · Submitted 2025-12-09 · Read on arXiv

Listen

Transcript

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

Kai: Today's paper: "Origin of superconductivity in bilayer nickelates".

Mira: Determinant Quantum Monte Carlo simulations investigate how doping, interlayer tunneling, and onsite Hund’s coupling stabilize superconductivity in bilayer Nickelate La3Ni2O7 by analyzing a sign-problem-free effective model.

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

Paper summary: Kai: So to wrap up our discussion on "Origin of superconductivity in bilayer nickelates: a Quantum Monte Carlo study for a sign-problem-free effective model," the paper essentially provides a path through the complexity of this system using specific symmetries.

Mira: Exactly, Kai; the authors are showing that by utilizing Kramers anti-unitary symmetries, they can create an effective model where they can avoid the sign problem in Quantum Monte Carlo simulations.

Lev: And from a quantum information perspective, if we can map these symmetry constraints onto physical Hamiltonians, it gives us a blueprint for constructing effective low-energy models that might be tractable for error correction simulations on real hardware.

Kai: The main implication is clarifying the role of inter-layer tunneling in deciding whether the system settles into superconductivity or exciton condensation at certain momentum points.

Mira: They also point toward directions on how we might enhance the superconducting transition temperature and stabilize that SC phase by modifying parameters within their framework.

Lev: If they suggest reducing intersubspace tunneling to modify the behavior, then that gives us a concrete theoretical target for experimentalists to aim for when tuning material parameters in synthesis. This paper, "Origin of superconductivity in bilayer nickelates: a Quantum Monte Carlo study for a sign-problem-free effective model," offers clarity on the mechanisms driving SC and competing instabilities in these complex materials.

Conclusion: Kai: So we've been looking at how they managed to tame that sign problem in their quantum Monte Carlo setup for bilayer nickelates, and now we get to talk about what this whole paper is actually about.

Mira: Exactly, Kai; the title itself tells us they focused on finding the origin of superconductivity through a method that avoids those pesky sign problems.

Lev: From a computational standpoint, avoiding the sign problem is huge because it means we can't just stop at theoretical predictions; it means we can actually run these models on hardware that has more than a few qubits.

Kai: Right, so if they successfully mapped the physics onto an effective model like this, what does that actually tell us about how these materials behave in a lab?

Mira: It tells us that by focusing on specific symmetries—those anti-unitary ones—they can distill the complex reality of the nickelate into something mathematically manageable.

Lev: And for error correction, if they can define a sign-problem-free effective Hamiltonian, it gives us a much clearer target for building quantum circuits that mimic these condensed matter systems.

Kai: So it seems like this work isn't just about finding another material property; it's about developing a new way to model complex strongly correlated systems using better computational tools.

Mira: It points toward a general strategy where exploiting underlying symmetries can unlock physical insights that were previously hidden behind mathematical intractability.

Lev: And the implication for real hardware is that if we can build these simulators, we get a direct testbed for how quantum algorithms handle these kinds of realistic electronic structures.

Kai: It's wild to think about the scale of what this means for materials science and computation when you consider all those layers of complexity they managed to untangle.

Mira: Indeed, the real impact here lies in showing that even with immense complexity, structure-preserving transformations can guide us toward a stable solution.

Lev: So we've seen how it works conceptually; now I'm curious if we can actually translate these symmetry constraints into something you could program onto a superconducting processor.

Department of Physics and Astronomy, Ghent University

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

Submitted: 2025-12-09

Updated: 2026-10-07

Comments: 7 pages, 4 figures

Journal ref: Phys. Rev. B 114, L231102 (2026)

DOI: 10.1103/h16z-tqcl

License: http://arxiv.org/licenses/nonexclusive-distrib/1.0/

Importance score: 83/100

The gist: Determinant Quantum Monte Carlo simulations investigate how doping, interlayer tunneling, and onsite Hund’s coupling stabilize superconductivity in bilayer Nickelate La3Ni2O7 by analyzing a

Key concepts

Quantum Monte Carlo Simulations: This is a computational method used to solve complex many-body problems in physics. In this study, it was used to investigate the electronic interactions in the bilayer nickelate material by simulating the system's behavior and finding stable superconducting states.
Sign-Problem-Free Effective Model: The original problem with Monte Carlo simulations is a 'sign problem,' which makes calculations impossible. This paper uses a specific mathematical transformation to create an effective model where this sign problem is absent, allowing for reliable simulation of the material's properties.
Competing Orders (SC vs. EC): The study analyzed different physical states the material could adopt, such as superconductivity (SC) and exciton condensation (EC). The results showed that SC is preferred over EC in certain conditions, but EC instability at specific momentum points $(\pi, \pi)$ is a strong candidate for the density wave seen in experiments.

Terminology

Quantum Monte Carlo — 19× in this paper
Hubbard Model — 4× in this paper · explained in 3 episodes
DFT Calculations — 2× in this paper · explained in 2 episodes
Trace — 2× in this paper · explained in 2 episodes
Antiferromagnetism — 1× in this paper · explained in 1 episode
Cooperation — 1× in this paper · explained in 2 episodes
Effective Hamiltonian — 1× in this paper · explained in 2 episodes
Entanglement — 1× in this paper · explained in 3 episodes
Hilbert Space — 1× in this paper · explained in 2 episodes
Phase Transitions — 1× in this paper · explained in 2 episodes
Quantum Circuits — 1× in this paper · explained in 2 episodes
Symmetry Breaking — 1× in this paper · explained in 3 episodes

Summary

Determinant Quantum Monte Carlo simulations investigate how doping, interlayer tunneling, and onsite Hund’s coupling stabilize superconductivity in bilayer Nickelate La3Ni2O7 by analyzing a sign-problem-free effective model.

The gist

Using a realistic band dispersion and doping for the bilayer Nickelate La3Ni2O7 [82], we find a SC phase with transition temperature Tc ≲ 0.06 coming from the dx2−y2 orbital, generated by a simplified interaction with Jzz = Jxz = 0.4.

Model Setup and Simulation Method

The study employs a Hamiltonian H = H0 + HI, where H0 describes the kinetic energy terms including nearest and next-nearest neighbor hopping between layers (l=1, 2) and orbitals (x for dx2−y2, z for dz2), while HI represents the interaction terms. The parameters in H0 are derived from DFT calculations [82], specifying hopping values such as t∥xx = −0.483 and t⊥zz = −0.635. The interaction term HI is approximated using a Hubbard-Stratonovich (HS) transformation, leading to an effective form amenable to simulation, where the parameters are tuned by fixing Jzz = 0.4 and varying Jxz from 0.4 to 0.1 in the simulations.

Symmetries and Sign Problem Resolution

A crucial aspect of this study is that the same anti-unitary symmetries which guarantee the absence of the sign problem also imply that the SC pairing is maximal in the ∆† = c†l xσ ym0K channel [199]. The paper details several Kramers anti-unitary symmetries, such as T1, T2, and T3, which are invariant under specific transformations. These symmetries ensure that the trace over fermionic degrees of freedom is positive definite for any auxiliary field configuration and hence gives rise to the well-defined probability weight for Monte Carlo sampling.

Competing Orders: SC vs. EC and SDW

The analysis of order parameters reveals a competition between different instabilities. The exchange contribution to the correlation function is maximal for the spin-singlet layer-triplet SC MSC = ⟨c†l xσ ym0/zc†⟩ due to the T1 symmetry, while for inter-layer EC MEC = ⟨c†l x/yσ 0m0/zc⟩ because of the T2, T3 symmetries. The direct contribution to the SC/EC order parameter correlation function is zero due to charge conservation. Adding back kinetic terms shows that the correlations of MSC are strictly larger than those of MEC, which indicates that SC is preferred over exciton condensation. Furthermore, when inter-layer tunneling t⊥zz is reinstated, the phase diagram is determined by the competition between the direct and exchange contributions, suggesting a transition from MSC = ⟨c†l xσ ym0 c†⟩ zero momentum order to MEC = ⟨c†l xσ 0m0 c⟩ order at momentum (π, π).

Filling Dependence and Conclusion

The study examines the filling dependence by tuning the chemical potential. It is found that superconductivity prefers small-doping region, and a large doping in dx2−y2 orbital will suppress and even eliminate the superconductivity. Turning off inter-layer tunneling (t⊥zz = 0) increases hole density on the dz2 orbitals, making Tc insensitive to a decrease in Hund’s coupling strength Jxz. Conversely, when reinstating t⊥zz to its non-zero (realistic) value, we find that decreasing the Hund’s coupling to Jxz = 0.2 eliminates SC. The findings suggest that the EC instability at (π, π) is a promising candidate for the unknown density wave observed in experiment and points toward avenues for enhancing superconductivity by strongly breaking particle-hole symmetry or reducing the intersubspace tunneling.

Spectral Features and Instabilities

Single-particle excitations show a clear single-particle gap developing around the X points in Fig. 1(c), which is reflected in the small value of G(k, τ = β/2) in Fig. 1(d). The static structure factors indicate that the EC structure factor is maximal at (π, π), in agreement with our analysis of the dominant ordering instabilities based on the Kramers anti-unitary symmetries, and a peak in the static spin structure factor indicates fluctuating SDW order near momentum close to (π2, π2). This suggests that fluctuating EC order is accompanied by fluctuating SDW order near momentum (π2, π2). The overall conclusion is that SC emerges from both types of orbitals, and the competition between direct and exchange contributions in the presence of inter-layer tunneling dictates the final phase diagram.

Improvements for AI systems

As a fastidious and diligent AI researcher, I have analyzed this paper, How a bilayer Nickelate superconducts: a Quantum Monte Carlo study. The core contribution lies in using determinant Quantum Monte Carlo (QMC) to map the phase diagram of La3Ni2O7 bilayer nickelates by studying doping, inter-layer tunneling, and Hund’s coupling.

Here are the specific improvements I can propose for AI systems based on this research:

The improved AI system will be a high-fidelity predictive model capable of simulating and interpreting strongly correlated electron systems, specifically focusing on emergent phenomena in transition metal oxides.

The system will incorporate a specialized QMC engine optimized for sign-problem-free fermion Hamiltonians (as established by the paper's use of Kramers anti-unitary symmetries).
The system will be trained on the derived two-orbital Hamiltonian, including realistic hopping parameters from DFT and interaction terms for Hund’s coupling and inter-layer exchange.

The improved AI system can perform the following specific tasks:

Predict the superconducting phase diagram (SC vs. EC vs. SDW) of nickelate materials by inputting varying doping levels and inter-layer tunneling strengths, allowing it to identify optimal material compositions for high-Tc superconductivity (e.g., predicting that large hole-doping in the Ni-3dx2−y2 orbital eliminates SC).
Determine the critical parameters for stabilizing specific pairing symmetries, such as identifying the competition between spin-singlet layer-triplet SC and competing exciton condensation (EC) order at momentum vectors like (π, π).
Predict the phase boundary behavior regarding inter-layer coupling: it can determine if reducing inter-layer tunneling is a viable strategy to suppress competing instabilities like EC, or conversely, how realistic tunneling values interact with Hund's coupling to either stabilize or destroy SC.
Analyze spectral functions and single-particle excitations (as shown in Fig. 1), enabling the AI to diagnose the presence of superconducting gaps and identify characteristic Fermi surface features under different doping conditions.
Identify the relationship between microscopic symmetries (Kramers anti-unitary symmetries) and macroscopic order parameters, allowing the AI to predict which pairing channels (e.g., spin-singlet layer-triplet channel vs. EC) are dominant in a given electronic environment.