Chiral d-wave superconductivity in doped graphene

Topological Superconductivity in Heavily Doped Single-Layer Graphene

The existence of superconductivity (SC) appears to be established in both twisted and nontwisted graphene multilayers. However, whether their building block, single-layer graphene (SLG), can also host SC remains an open question. Earlier theoretical works predicted that SLG could become a chiral d-wave superconductor driven by electronic interactions when doped to its van Hove singularity, but questions such as whether the d-wave SC survives the strong band renormalizations seen in experiments, its robustness against the source of doping, or if it will occur at any reasonable critical temperature (Tc) have remained difficult to answer, in part due to uncertainties in model parameters. Furthermore, doping of graphene beyond its van Hove singularity remained experimentally challenging and was not demonstrated until recently. In this study, we n dope SLG past the van Hove singularity by employing Tb intercalation and derive structural models from angle-resolved photoemission spectroscopy measurements. We adopt a reliable numerical framework based on a random-phase approximation technique to investigate the emergence of unconventional SC in the heavily doped monolayer. We predict that robust d + id topological SC could arise in SLG doped by Tb, with a Tc up to 600 mK. We also employ first-principles calculations to investigate the possibility of realizing d-wave SC with other dopants, such as Li or Cs. We find that dopants that change the lattice symmetry of SLG are detrimental to the d-wave state. The stability of the d-wave SC predicted here in Tb-doped SLG could provide a valuable insight for guiding future experimental efforts aimed at exploring topological superconductivity in monolayer graphene.

Van Hove Singularity and Enhanced Superconductivity in Ca-Intercalated Bilayer Graphene Induced by Confinement Epitaxy

We demonstrate the impact of high-density calcium introduction into Ca-intercalated bilayer graphene on a SiC substrate, wherein a metallic layer of Ca has been identified at the interface. We have discerned that the additional Ca layer engenders a free-electron-like band, which subsequently hybridizes with a Dirac band, leading to the emergence of a van Hove singularity. Coinciding with this, there is an increase in the critical temperature for superconductivity. These findings allude to the manifestation of Ca-driven confinement epitaxy, augmenting superconductivity through the enhancement of attractive interactions in a pair of electron and hole bands with flat dispersion around the Fermi level.

Chiral-Flux-Phase-Based Topological Superconductivity in Kagome Systems with Mixed Edge Chiralities

Recent studies have attracted intense attention on the quasi-2D kagome superconductors AV_{3}Sb_{5} (A=K, Rb, and Cs) where the unexpected chiral flux phase (CFP) associates with the spontaneous time-reversal symmetry breaking in charge density wave states. Here, commencing from the 2-by-2 charge density wave phases, we bridge the gap between topological superconductivity and time-reversal asymmetric CFP in kagome systems. Several chiral topological superconductor (TSC) states featuring distinct Chern numbers emerge for an s-wave or a d-wave superconducting pairing symmetry. Importantly, these CFP-based TSC phases possess unique gapless edge modes with mixed chiralities (i.e., both positive and negative chiralities), but with the net chiralities consistent with the Bogoliubov-de Gennes Chern numbers. We further study the transport properties of a two-terminal junction, using Chern insulator or normal metal leads via atomic Green's function method with Landauer-Büttiker formalism. In both cases, the normal electron tunneling and the crossed Andreev reflection oscillate as the chemical potential changes, but together contribute to plateau transmissions (1 and 3/2, respectively) that exhibit robustness against disorder. These behaviors can be regarded as the signature of a TSC hosting edge states with mixed chiralities.

Andreev Reflection and Klein Tunneling in High-Temperature Superconductor-Graphene Junctions

Scattering processes in quantum materials emerge as resonances in electronic transport, including confined modes, Andreev states, and Yu-Shiba-Rusinov states. However, in most instances, these resonances are driven by a single scattering mechanism. Here, we show the appearance of resonances due to the combination of two simultaneous scattering mechanisms, one from superconductivity and the other from graphene p-n junctions. These resonances stem from Andreev reflection and Klein tunneling that occur at two different interfaces of a hole-doped region of graphene formed at the boundary with superconducting graphene due to proximity effects from Bi_{2}Sr_{2}Ca_{1}Cu_{2}O_{8+δ}. The resonances persist with gating from p^{+}-p and p-n configurations. The suppression of the oscillation amplitude above the bias energy which is comparable to the induced superconducting gap indicates the contribution from Andreev reflection. Our experimental measurements are supported by quantum transport calculations in such interfaces, leading to analogous resonances. Our results put forward a hybrid scattering mechanism in graphene-high-temperature superconductor heterojunctions of potential impact for graphene-based Josephson junctions.

Pair-Density-Wave and Chiral Superconductivity in Twisted Bilayer Transition Metal Dichalcogenides

We theoretically explore possible orders induced by weak repulsive interactions in twisted bilayer transition metal dichalcogenides (e.g., WSe_{2}) in the presence of an out-of-plane electric field. Using renormalization group analysis, we show that superconductivity survives even with the conventional van Hove singularities. We find that topological chiral superconducting states with Chern number N=1, 2, 4 (namely, p+ip, d+id, and g+ig) appear over a large parameter region with a moiré filling factor around n=1. At some special values of applied electric field and in the presence of a weak out-of-plane Zeeman field, spin-polarized pair-density-wave (PDW) superconductivity can emerge. This spin-polarized PDW state can be probed by experiments such as spin-polarized STM measuring spin-resolved pairing gap and quasiparticle interference. Moreover, the spin-polarized PDW could lead to a spin-polarized superconducting diode effect.

Inter-valley coherent order and isospin fluctuation mediated superconductivity in rhombohedral trilayer graphene

Superconductivity was recently discovered in rhombohedral trilayer graphene (RTG) in the absence of a moiré potential. Superconductivity is observed proximate to a metallic state with reduced isospin symmetry, but it remains unknown whether this is a coincidence or a key ingredient for superconductivity. Using a Hartree-Fock analysis and constraints from experiments, we argue that the symmetry breaking is inter-valley coherent (IVC) in nature. We evaluate IVC fluctuations as a possible pairing glue, and find that they lead to chiral unconventional superconductivity when the fluctuations are strong. We further elucidate how the inter-valley Hund’s coupling determines the spin-structure of the IVC ground state and breaks the degeneracy between spin-singlet and triplet superconductivity. Remarkably, if the normal state is spin-unpolarized, we find that a ferromagnetic Hund’s coupling favors spin-singlet superconductivity, in agreement with experiments. Instead, if the normal state is spin-polarized, then IVC fluctuations lead to spin-triplet pairing.

Trilayer graphene was recently shown to exhibit superconductivity without a Moire pattern that had proved important in tuning superconductivity in bilayer graphene. Here, the authors explore correlated metallic phases and the pairing mechanism of superconductivity in trilayer graphene, and show that intervalley coherent fluctuations can act as a pairing glue, giving rise to chiral unconventional superconductivity.

Triplet Superconductivity from Nonlocal Coulomb Repulsion in an Atomic Sn Layer Deposited onto a Si(111) Substrate

Atomic layers deposited on semiconductor substrates introduce a platform for the realization of the extended electronic Hubbard model, where the consideration of electronic repulsion beyond the on-site term is paramount. Recently, the onset of superconductivity at 4.7 K has been reported in the hole-doped triangular lattice of tin atoms on a silicon substrate. Through renormalization group methods designed for weak and intermediate coupling, we investigate the nature of the superconducting instability in hole-doped Sn/Si(111). We find that the extended Hubbard nature of interactions is crucial to yield triplet pairing, which is f-wave (p-wave) for moderate (higher) hole doping. In light of persisting challenges to tailor triplet pairing in an electronic material, our finding promises to pave unprecedented ways for engineering unconventional triplet superconductivity.

Two-Dimensional Superconductivity of Ca-Intercalated Graphene on SiC: Vital Role of the Interface between Monolayer Graphene and the Substrate

Ca-intercalation has enabled superconductivity in graphene on SiC. However, the atomic and electronic structures that are critical for superconductivity are still under discussion. We find an essential role of the interface between monolayer graphene and the SiC substrate for superconductivity. In the Ca-intercalation process, at the interface a carbon layer terminating SiC changes to graphene by Ca-termination of SiC (monolayer graphene becomes a bilayer), inducing more electrons than a free-standing model. Then, Ca is intercalated in between the graphene layers, which shows superconductivity with the updated critical temperature (T_C) of up to 5.7 K. In addition, the relation between T_C and the normal-state conductivity is unusual, "dome-shaped". These findings are beyond the simple C₆CaC₆ model in which s-wave BCS superconductivity is theoretically predicted. This work proposes a general picture of the intercalation-induced superconductivity in graphene on SiC and indicates the potential for superconductivity induced by other intercalants.

Inducing chiral superconductivity on honeycomb lattice systems

Superconductivity in graphene-based systems has recently attracted much attention, as either intrinsic behavior or induced by proximity to a superconductor may lead to interesting topological phases and symmetries of the pairing function. A prominent system considers the pairing to have chiral symmetry. The question arises as to the effect of possible spin-orbit coupling on the resulting superconducting quasiparticle (QP) spectrum. Utilizing a Bogolyubov-de Gennes (BdG) Hamiltonian, we explore the interplay of different interaction terms in the system, and their role in generating complex Berry curvatures in the QP spectrum, as well as non-trivial topological behavior. We demonstrate that the topology of the BdG Hamiltonian in these systems may result in the appearance of edge states along the zigzag edges of nanoribbons in the appropriate regime. For suitable chemical potential and superconducting pairing strength, we find the appearance of robust midgap states at zigzag edges, well protected by large excitation gaps and momentum transfer.

Unconventional Superconductivity in Systems with Annular Fermi Surfaces: Application to Rhombohedral Trilayer Graphene

We show that in a two-dimensional electron gas with an annular Fermi surface, long-range Coulomb interactions can lead to unconventional superconductivity by the Kohn-Luttinger mechanism. Superconductivity is strongly enhanced when the inner and outer Fermi surfaces are close to each other. The most prevalent state has chiral p-wave symmetry, but d-wave and extended s-wave pairing are also possible. We discuss these results in the context of rhombohedral trilayer graphene, where superconductivity was recently discovered in regimes where the normal state has an annular Fermi surface. Using realistic parameters, our mechanism can account for the order of magnitude of T_{c}, as well as its trends as a function of electron density and perpendicular displacement field. Moreover, it naturally explains some of the outstanding puzzles in this material, that include the weak temperature dependence of the resistivity above T_{c}, and the proximity of spin singlet superconductivity to the ferromagnetic phase.

Nematicity Arising from a Chiral Superconducting Ground State in Magic-Angle Twisted Bilayer Graphene under In-Plane Magnetic Fields

Recent measurements of the resistivity in magic-angle twisted bilayer graphene near the superconducting transition temperature show twofold anisotropy, or nematicity, when changing the direction of an in-plane magnetic field [Cao et al., Science 372, 264 (2021)SCIEAS0036-807510.1126/science.abc2836]. This was interpreted as strong evidence for exotic nematic superconductivity instead of the widely proposed chiral superconductivity. Counterintuitively, we demonstrate that in two-dimensional chiral superconductors the in-plane magnetic field can hybridize the two chiral superconducting order parameters to induce a phase that shows nematicity in the transport response. Its paraconductivity is modulated as cos(2θ_{B}), with θ_{B} being the direction of the in-plane magnetic field, consistent with experiment in twisted bilayer graphene. We therefore suggest that the nematic response reported by Cao et al. does not rule out a chiral superconducting ground state.

Solid-state quantum chemistry with ΘΦ (ThetaPhi): Spin-liquids, superconductors, and magnetic superstructures made computationally available

We present a standalone ΘΦ (ThetaPhi) package capable to read the results of ab initio DFT/PAW quantum-chemical solid-state calculations processed through various tools projecting them to the atomic basis states as an input and to perform on top of this an analysis of so derived electronic structure which includes (among other options) the possibility to obtain a superconducting (Bardeen-Cooper-Schrieffer, BCS), spin-liquid (resonating valence bond, RVB) states/phases as solutions of the electronic structure problem along with the magnetically ordered phases with an arbitrary pitch (magnetic superstructure) vector. Remarkably, different solutions of electronic-structure problems come out as temperature-dependent (exemplified by various superconducting and spin-liquid phases) which feature is as well implemented. All that is exemplified by model calculations on 1D chain, 2D square lattice as well as on more realistic superconducting doped graphene, magnetic phases of iron, and spin-liquid and magnetically ordered states of a simplest nitrogen-based copper pseudo-oxide, CuNCN, resembling socalled metal-oxide framework (MOF) phases by the atomic interlinkage.

Berry Curvature Effects on Quasiparticle Dynamics in Superconductors

We construct a theory for the semiclassical dynamics of superconducting quasiparticles by following their wave packet motion and reveal rich contents of Berry curvature effects in the phase space spanned by position and momentum. These Berry curvatures are traced back to the characteristics of superconductivity, including the nontrivial momentum-space geometry of superconducting pairing, the real-space supercurrent, and the charge dipole of quasiparticles. The Berry-curvature effects strongly influence the spectroscopic and transport properties of superconductors, such as the local density of states and the thermal Hall conductivity. As a model illustration, we apply the theory to study the twisted bilayer graphene with a d_{x^{2}+y^{2}}+id_{xy} superconducting gap function and demonstrate Berry-curvature induced effects.

Overdoping Graphene beyond the van Hove Singularity

At very high doping levels the van Hove singularity in the π^{*} band of graphene becomes occupied and exotic ground states possibly emerge, driven by many-body interactions. Employing a combination of ytterbium intercalation and potassium adsorption, we n dope epitaxial graphene on silicon carbide past the π^{*} van Hove singularity, up to a charge carrier density of 5.5×10^{14}  cm^{-2}. This regime marks the unambiguous completion of a Lifshitz transition in which the Fermi surface topology has evolved from two electron pockets into a giant hole pocket. Angle-resolved photoelectron spectroscopy confirms these changes to be driven by electronic structure renormalizations rather than a rigid band shift. Our results open up the previously unreachable beyond-van-Hove regime in the phase diagram of epitaxial graphene, thereby accessing an unexplored landscape of potential exotic phases in this prototype two-dimensional material.

Geometric Induction in Chiral Superconductors

We consider a number of effects due to the interplay of superconductivity, electromagnetism, and elasticity, which are unique for thin membranes of layered chiral superconductors. Some of them should be within the reach of present technology, and could be useful for characterizing materials. More speculatively, the enriched control of Josephson junctions they afford might find useful applications.

Charge Order in the Holstein Model on a Honeycomb Lattice

The effect of electron-electron interactions on Dirac fermions, and the possibility of an intervening spin-liquid phase between the semimetal and antiferromagnetic (AF) regimes, has been a focus of intense quantum simulation effort over the last five years. We use determinant quantum Monte Carlo simulations to study the Holstein model on a honeycomb lattice and explore the role of electron- phonon interactions on Dirac fermions. We show that they give rise to charge-density-wave (CDW) order and present evidence that this occurs only above a finite critical interaction strength. We evaluate the temperature for the transition into the CDW which, unlike the AF transition, can occur at finite values owing to the discrete nature of the broken symmetry.

Proposed Spontaneous Generation of Magnetic Fields by Curved Layers of a Chiral Superconductor

We demonstrate that two-dimensional chiral superconductors on curved surfaces spontaneously develop magnetic flux. This geometric Meissner effect provides an unequivocal signature of chiral superconductivity, which could be observed in layered materials under stress. We also employ the effect to explain some puzzling questions related to the location of zero-energy Majorana modes.

Bogoliubov Fermi Surfaces in Superconductors with Broken Time-Reversal Symmetry

It is commonly believed that, in the absence of disorder or an external magnetic field, there are three possible types of superconducting excitation gaps: The gap is nodeless, it has point nodes, or it has line nodes. Here, we show that, for an even-parity nodal superconducting state which spontaneously breaks time-reversal symmetry, the low-energy excitation spectrum generally does not belong to any of these categories; instead, it has extended Bogoliubov Fermi surfaces. These Fermi surfaces can be visualized as two-dimensional surfaces generated by "inflating" point or line nodes into spheroids or tori, respectively. These inflated nodes are topologically protected from being gapped by a Z_{2} invariant, which we give in terms of a Pfaffian. We also show that superconducting states possessing these Fermi surfaces can be energetically stable. A crucial ingredient in our theory is that more than one band is involved in the pairing; since all candidate materials for even-parity superconductivity with broken time-reversal symmetry are multiband systems, we expect these Z_{2}-protected Fermi surfaces to be ubiquitous.

Helical Majorana fermions in d(x2-y2) + id(xy)-wave topological superconductivity of doped correlated quantum spin Hall insulators

There has been growing interest in searching for exotic self-conjugate, charge-neutral low-energy fermionic quasi-particles, known as Majorana fermions (MFs) in solid state systems. Their signatures have been proposed and potentially observed at edges of topological superconcuctors with non-trivial topological invariant in the bulk electronic band structure. Much effort have been focused on realizing MFs in odd-parity superconductors made of strong spin-orbit coupled materials in proximity to conventional superconductors. In this paper, we propose a novel mechanism for realizing MFs in 2D spin-singlet topological superconducting state induced by doping a correlated quantum spin Hall (Kane-Mele) insulator. Via a renormalized mean-field approach, the system is found to exhibits time-reversal symmetry (TRS) breaking d(x2-y2) + id(xy)-wave (chiral d-wave) superconductivity near half-filling in the limit of large on-site repulsion. Surprisingly, however, at large spin-orbit coupling, the system undergoes a topological phase transition and enter into a new topological phase protected by a pseudo-spin Chern number, which can be viewed as a persistent extension of the quantum spin Hall phase upon doping. From bulk-edge correspondence, this phase is featured by the presence of two pairs of counter-propagating helical Majorana modes per edge, instead of two chiral propagating edge modes in the d + id' superconductors.

Accessing topological superconductivity via a combined STM and renormalization group analysis

The search for topological superconductors has recently become a key issue in condensed matter physics, because of their possible relevance to provide a platform for Majorana bound states, non-Abelian statistics, and quantum computing. Here we propose a new scheme which links as directly as possible the experimental search to a material-based microscopic theory for topological superconductivity. For this, the analysis of scanning tunnelling microscopy, which typically uses a phenomenological ansatz for the superconductor gap functions, is elevated to a theory, where a multi-orbital functional renormalization group analysis allows for an unbiased microscopic determination of the material-dependent pairing potentials. The combined approach is highlighted for paradigmatic hexagonal systems, such as doped graphene and water-intercalated sodium cobaltates, where lattice symmetry and electronic correlations yield a propensity for a chiral singlet topological superconductor. We demonstrate that our microscopic material-oriented procedure is necessary to uniquely resolve a topological superconductor state.