Variational Monte Carlo Optimization of Topological Chiral Superconductors

Using variational Monte Carlo calculations, this study demonstrates that topological chiral superconducting states driven by strong repulsive Coulomb interactions are energetically favored over Fermi liquid phases in rhombohedral graphene systems with specific dispersion relations, revealing a new mechanism for superconductivity that does not rely on conventional Fermi surface pairing instabilities.

The Gist

The Big Picture: A New Kind of Dance

Imagine a crowded dance floor where the dancers are electrons. Usually, when these dancers get too close, they push each other away because they all have the same electric charge (like trying to hug someone who hates you). This is called repulsion.

In most superconductors (materials that conduct electricity with zero resistance), the dancers pair up like couples holding hands to glide smoothly across the floor. This usually happens because of a specific "glue" (often vibrations in the material) that helps them overcome their natural urge to push apart. This is the famous BCS theory.

However, scientists recently discovered something strange in a special type of stacked graphene (a form of carbon). They saw a Chiral Superconductor. "Chiral" means the dance has a specific direction, like everyone spinning clockwise. Even stranger, this dance seems to happen without the usual "glue." It seems to be driven purely by the dancers' intense desire to stay apart from each other.

The Question: Is this new dance actually possible? Does it cost less energy than the dancers just standing around in a chaotic crowd (a "Fermi liquid") or forming a rigid crystal (a "Wigner crystal")?

The Experiment: The Virtual Dance Floor

The authors of this paper, Minho Luke Kim, Abigail Timmel, and Xiao-Gang Wen, didn't build a physical lab. Instead, they built a super-computer simulation (Variational Monte Carlo) to act as a virtual dance floor.

They wanted to test three specific types of "chiral dances" (superconducting states) against the standard "chaotic crowd" (the Fermi liquid) to see which one is the most energy-efficient.

The Three Dancers (The States)

1. The Pfaffian (The Solo Spin): Imagine a dance where everyone is forced to spin in the exact same direction (spin-polarized). They move in a complex, coordinated pattern that avoids bumping into each other.
2. The K2a (The Mixed Spin): Imagine a dance where half the crowd spins clockwise and half spins counter-clockwise, but they still manage to coordinate a perfect, flowing dance together.
3. The K2b (The Exotic Spin): A more complex, rare dance pattern that is very different from anything we usually see.

The Terrain: The "Flat" Floor

The dance floor isn't perfectly flat. It has a specific shape defined by a mathematical curve. The authors found that the shape of the floor matters immensely.

• If the floor is a gentle hill, the dancers prefer to just stand around (the Fermi liquid).
• The Sweet Spot: If the floor is shaped like a shallow bowl that is almost flat at the bottom, the dancers suddenly prefer to start their special chiral dance. It's like the floor is so flat that the only way to move without bumping into each other is to spin in a coordinated circle.

The Results: Who Wins the Dance-Off?

The computer simulation calculated the "energy cost" of every scenario. In physics, the state with the lowest energy is the one that actually happens in nature.

1. The Surprise: The authors found that the Chiral Superconducting dances (Pfaffian and K2a) actually cost less energy than the chaotic crowd (Fermi liquid) in certain conditions.
2. The "Almost Flat" Key: This only happens when the "floor" (the electron energy band) is almost flat at the bottom. This is a crucial finding because it means you don't need a mysterious "glue" to make superconductors. You just need a specific shape of the floor and strong repulsion.
3. The Magnetic Field Test: They also tested what happens if you put a giant magnet near the dance floor.
   • The "Solo Spin" dance (Pfaffian) is very sensitive to magnets.
   • The "Mixed Spin" dance (K2a) is very robust. It keeps dancing even when a strong magnet is applied. This matches real-world experiments where the superconductivity didn't disappear even under strong magnetic fields.

The "Aha!" Moment: Why This Matters

This paper suggests a new pathway to superconductivity.

• Old Way (BCS): You need a "glue" to pair up electrons.
• New Way (This Paper): You don't need a glue. If you have a system where electrons really hate each other (strong repulsion) and the energy landscape is "flat" just right, the electrons will spontaneously organize into a superconducting dance to avoid each other.

It's like a group of people who hate being crowded. Instead of standing still (which is boring and takes energy), they decide to run in a perfect circle. By moving together in a specific pattern, they actually avoid collisions better than if they were just standing still.

Summary in One Sentence

The authors used a super-computer to prove that electrons, driven purely by their mutual dislike of each other, can spontaneously form a special, spinning superconducting dance on a nearly flat energy floor, offering a brand new way to create superconductors without the usual "glue."

Technical Summary

1. Problem Statement

The paper addresses the origin of chiral superconductivity recently observed in rhombohedral stacked graphene systems (specifically multilayer graphene).

• Context: Experiments have identified superconducting states between the spin-valley polarized Fermi liquid (Quarter Fermi Liquid, QFL) and the Wigner crystal phase. These states exhibit broken time-reversal symmetry (chirality) and robustness against magnetic fields.
• Theoretical Challenge: While many theories attribute this to BCS-like pairing instabilities on a Fermi surface, the observed short coherence length and the fact that the neighboring phases (QFL and Wigner crystal) are driven by strong repulsive Coulomb interactions suggest a different mechanism.
• Core Question: Can purely repulsive Coulomb interactions, without relying on Fermi surface pairing instabilities, energetically stabilize topological chiral superconducting states over the standard Fermi liquid? Previous works (e.g., Ref. 17) proposed such states but lacked a rigorous energetic comparison with an optimized Fermi liquid wavefunction.

2. Methodology

The authors employ Variational Monte Carlo (VMC) calculations to compare the ground state energies of proposed chiral superconducting states against a highly optimized Quarter Fermi Liquid (QFL).

• Model System:

  • Dispersion Relation: $E_k = c_2 k^2 + c_4 k^4$. This quartic dispersion is crucial for modeling the flat-band bottom and the tunable band curvature in multilayer graphene via displacement fields.
  • Interactions: Strong repulsive Coulomb interactions ($e^2/\epsilon r$).
  • Parameters: Electron density ($n_e$) and the quadratic coefficient ($c_2$), which is tuned experimentally.

• Trial Wavefunctions:

  1. Chiral Superconductors:
     • Single-species (Spin-polarized): A modified Pfaffian wavefunction. The authors introduce a new ansatz (Eq. 7) that replaces the Gaussian decay of the standard Laughlin-related ansatz with algebraic long-range density correlations. This allows for better optimization of potential energy while maintaining the topological structure.
     • Two-species (Spin-unpolarized): $K_{2a}$ and $K_{2b}$ states. These are based on Laughlin-type wavefunctions with two species of electrons (spin/valley). $K_{2a}$ corresponds to a spin-triplet $p+ip$ phase, while $K_{2b}$ carries fractional statistics.
  2. Quarter Fermi Liquid (QFL):
     • The authors construct a Slater-Jastrow wavefunction ($\Psi = D(R)e^{J(R)}$) for the QFL.
     • Unlike previous studies limited to quadratic dispersion, they optimize the Jastrow factor for the general quartic dispersion ($c_4 \neq 0$) to accurately capture correlation energy.

• Optimization & Calculation:

  • Variational parameters (e.g., Jastrow coefficients, matrix elements of the $K$-matrix, decay lengths) are continuously tuned to minimize the total energy.
  • Kinetic Energy: Calculated via numerical derivatives of the wavefunction. For the quartic term ($k^4$), the authors utilize integration by parts on a torus to sample $|\nabla^2 \psi / \psi|^2$, avoiding direct computation of the biharmonic operator.
  • Potential Energy: Computed by numerically extracting the pair distribution function $g(z)$ and integrating with the Coulomb potential, including Ewald summation for periodic boundary conditions and a neutralizing background.

3. Key Contributions

• Energetic Comparison: This is the first study to perform a fully optimized VMC comparison between topological chiral superconductors and a correlated Fermi liquid with quartic dispersion.
• New Pfaffian Ansatz: The authors propose a superior trial wavefunction for the single-species Pfaffian state that incorporates algebraic long-range correlations, significantly lowering the energy compared to standard Laughlin-based ansatzes.
• QFL Optimization: They provide a rigorous calculation of the QFL ground state energy for general $E_k = c_2 k^2 + c_4 k^4$, showing that optimizing the Jastrow factor significantly lowers the Fermi liquid energy, making the competition with superconductors much tighter and more realistic.
• Mechanism Identification: They demonstrate that superconductivity can emerge from pure repulsive interactions in systems with an almost flat band bottom, bypassing the conventional BCS pairing instability mechanism.

4. Key Results

• Phase Diagram (Fig. 1): The authors map the ground state phase diagram as a function of electron density ($n_e$) and band curvature ($c_2$).
  • Favorable Regime: Topological chiral superconducting phases (both Pfaffian and $K_{2a}$) are energetically favored over the QFL when $c_2$ is slightly negative (approaching the formation of a hole pocket at $k=0$) and at densities relevant to experiments ($0.2 \sim 0.7 \times 10^{12} \text{cm}^{-2}$).
  • Critical Condition: The preference for chiral superconductivity is strongest when the system is on the verge of a Fermi surface transition (where a hole pocket forms).
• Energy Gains:
  • The superconducting condensation energy is approximately 1 meV (roughly 5% of the Coulomb energy scale $\sqrt{n_e}e^2/\epsilon$).
  • The Pfaffian state is generally more favorable at lower densities, while the $K_{2a}$ state (spin-unpolarized) remains competitive at higher densities.
  • The energy difference between the Pfaffian and $K_{2a}$ states is very small, suggesting both could be realized depending on specific material parameters.
• Magnetic Field Robustness:
  • For the single-species Pfaffian and QFL, the Zeeman energy shifts are identical, so the energy difference remains unchanged to first order.
  • For the two-species $K_{2a}$ state, the Zeeman shifts cancel out (half spin up, half spin down), while the QFL energy lowers. However, the authors find that the $K_{2a}$ state remains robust against magnetic fields up to 5 Tesla, consistent with experimental observations in rhombohedral graphene.

5. Significance

• Beyond BCS: The results provide strong theoretical evidence that superconductivity in flat-band graphene systems can arise from strongly correlated repulsive interactions rather than weak-coupling BCS pairing. This establishes a new pathway to superconductivity.
• Topological Order: The study validates the existence of topological chiral superconducting states (with Majorana or complex fermion edge modes) as the true ground state in specific regimes of multilayer graphene.
• Experimental Guidance: The findings explain why chiral superconductivity is observed near the Fermi surface transition line in rhombohedral graphene. The prediction that the state is robust against magnetic fields supports the experimental observation of superconductivity persisting up to 5T.
• Future Directions: The paper highlights the need to include Berry curvature effects (which may further lower the energy of chiral states) and to study the energetics of the Wigner crystal phase to complete the full phase diagram.

In summary, the paper successfully demonstrates that repulsive Coulomb interactions alone can stabilize topological chiral superconductivity in multilayer graphene, offering a compelling explanation for recent experimental discoveries in rhombohedral graphene systems.