Variational Monte Carlo Optimization of Topological Chiral Superconductors

Abstract

We perform the variational Monte Carlo calculation for recently proposed chiral superconducting states driven by strong Coulomb interactions. We compare the resulting energetics of these electronic phases for the electron dispersion relation $E_k = c_2 k^2+c_4 k^4$. Motivated by the recent discovery of chiral superconductivity in rhombohedral graphene systems, we apply our analysis to relevant parameter regimes. We demonstrate that topological chiral superconducting phases (including a spin-unpolarized state) can be energetically favored over the spin-valley polarized Fermi liquid above the density of Wigner crystal phase. Our results show that the preference for chiral superconductivity is strongest when $c_2$ lies between zero and a negative value, corresponding to a system on the verge of forming a hole pocket around $k=0$. This finding suggests that superconductivity can arise from pure repulsive Coulomb interactions in systems with an almost flat band bottom, without relying on the pairing instability of a Fermi surface. This mechanism opens a new pathway to superconductivity beyond the conventional BCS mechanism.