Search

ABC-stacked multilayer graphene (ABC-MLG) exhibits topological surface flat bands with a divergent density of states, leading to many-body instabilities at charge neutrality. Here, we explore electronic ordering within a mean-field approach with full generic treatment of all spin-isotropic, two-site charge density and spin interactions up to next-nearest neighbor (NNN) sites. We find that surface superconductivity and magnetism are significantly enhanced over bulk values. We find spin-singlet $s$ wave and unconventional NNN bond spin-triplet $f$ wave to be the dominant superconducting pairing symmetries, both with a full energy gap. By establishing the existence of ferromagnetic intra-sublattice interaction, $(J_2<0)$ we conclude that the $f$-wave state is favored in ABC-MLG, in sharp contrast to bulk ABC-graphite where chiral $d$- or $p$-wave states, together with s-wave states, display stronger ordering tendencies albeit not achievable at charge neutrality. We trace this distinctive surface behavior to the strong sublattice polarization of the surface flat bands. We also find competing ferrimagnetic order, fully consistent with density functional theory (DFT) calculations. The magnetic order interpolates between sublattice ferromagnetism and antiferromagnetism, but only with the ratio of the sublattice magnetic moments ($R$) being insensitive to the DFT exchange correlation functional. We finally establish the full phase diagram by constraining the interactions to the $R$-value identified by DFT. We find $f$-wave superconductivity being favored for all weak to moderately strong couplings $J_2$ and as long as $J_2$ is a sufficiently large part of the full interaction mix. Gating ABC-MLG away from charge neutrality further enhances the $f$-wave state over the ferrimagnetic state, establishing ABC-MLG as a strong candidate for $f$-wave superconductivity., Comment: 20 pages, 11 figures + supplementary (1 page). Minor corrections implemented

We investigate the effects of single, multiple, and extended defects in the form of non-magnetic impurities and vacancies in twisted bilayer graphene (TBG) at and away from the magic angle, using a fully atomistic model and focusing on the behavior of the flat low-energy moir\'e bands. For strong impurities and vacancies in the $AA$ region we find a complete removal of one of the four moir\'e bands, resulting in a significant depletion of the charge density in the $AA$ regions even at extremely low defect concentrations. We find similar results for other defect locations, with the exception of the least coordinated sites in the $AB$ region, where defects instead result in a peculiar band replacement process within the moir\'e bands. In the vacancy limit, this process yields a band structure misleadingly similar to the pristine case. Moreover, we show that triple point fermions (TPFs), which are the crossing of the Dirac point by a flat band, appearing for single, periodic, defects, are generally not preserved when adding extended or multiple defects, and thus likely not experimentally relevant. We further identify two universal length scales for defects, consisting of charge modulations on the atomic scale and on the moir\'e scale, illustrating the importance of both the atomic and moir\'e structures for understanding TBG. We show that our conclusions hold beyond the magic angle and for fully isolated defects. In summary, our results demonstrate that the normal state of TBG and its moir\'e flat bands are extremely sensitive to both the location and strength of non-magnetic impurities and vacancies, which should have significant implications for any emergent ordered state., Comment: 15 pages, 10 figures

Twisted bilayer graphene (TBG) develops large moir\'e patterns at small twist angles with flat energy bands hosting domes of superconductivity. The large system size and intricate band structure have however hampered investigations into the superconducting state. Here, using full-scale atomistic modelling with local electronic interactions, we find at and above experimentally relevant temperatures a highly inhomogeneous superconducting state with nematic ordering on both atomic and moir\'e length scales. The nematic state has a locally anisotropic real-valued d-wave pairing, with a nematic vector winding throughout the moir\'e pattern, and is three-fold degenerate. Although d-wave symmetric, the superconducting state has a full energy gap, which we tie to a {\pi}-phase interlayer coupling. The superconducting nematicity is further directly detectable in the local density of states. Our results show that atomistic modeling is essential and also that very similar local interactions produce very different superconducting states in TBG and the high-temperature cuprate superconductors., Comment: 4 figures

In this work we develop EPOCH: Equilibrium Propagator by Orthogonal polynomial CHain, a computationally efficient method to calculate the time-dependent equilibrium Green's functions, including the anomalous Green's functions of superconductors, to capture the time-evolution in large inhomogeneous systems. The EPOCH method generalizes the Chebyshev wave-packet propagation method from quantum chemistry and efficiently incorporates the Fermi-Dirac statistics that is needed for equilibrium quantum condensed matter systems. The computational cost of EPOCH scales only linearly in the system degrees of freedom, generating an extremely efficient algorithm also for very large systems. We demonstrate the power of the EPOCH method by calculating the time-evolution of an excitation near a superconductor-normal metal interface in two and three dimensions, capturing transmission as well as normal and Andreev reflections., Comment: 17 pages, 5 figures. Published version

We investigate the induced superconducting pair correlations in junctions between a conventional spin-singlet $s$-wave superconductor and a disordered normal metal. Decomposing the pair amplitude based on its symmetries in the time domain, we demonstrate that the odd-time, or equivalently odd-frequency, spin-singlet $p$-wave correlations are both significant in size and entirely robust against random non-magnetic disorder. We find that these odd-frequency correlations can even be generated by disorder. Our results show that anisotropic odd-frequency pairing represent an important fraction of the proximity-induced correlations in disordered superconducting hybrid structures., Comment: 7 pages, 4 figures + 7 pages supplemental material, 4 figures. Published version

Condensed matter systems with flat bands close to the Fermi level generally exhibit, due to their very large density of states, extraordinary high critical ordering temperatures of symmetry breaking orders, such as superconductivity and magnetism. Here we show that the critical temperatures follow one of two universal curves with doping away from a flat band depending on the ordering channel, which completely dictates both the general order competition and the phase diagram. Notably, we find that orders in the particle-particle channel (superconducting orders) survive decisively further than orders in the particle-hole channel (magnetic or charge orders) because the channels have fundamentally different polarizabilities. Thus, even if a magnetic or charge order initially dominates, superconducting domes are still likely to exist on the flanks of flat bands. We apply these general results to both the topological surface flat bands of rhombohedral ABC-stacked graphite and to the van Hove singularity of graphene.

A chiral time-reversal symmetry breaking (d+id)-wave superconducting state is likely to emerge in graphene doped close to the Van Hove singularity. As heavy doping procedures are expected to introduce defects, we here investigate the effects of microscopic defects on the (d+id)-wave superconducting state at the Van Hove singularity. We find that, while the superconducting order is reduced near a defect, the (d+id)-wave state remains intact and recovers in an exponential manner away from the defect. The recovery length is found to be on the order of one lattice constant for weak couplings. This is comparable to the recovery length of a conventional s-wave state, demonstrating that the unconventional (d+id)-wave state is quite resilient to defects. Moreover, we find no significant changes between a single site defect and more extended defects, such as a bivacancy. While the (d+id)-wave state is fully gapped, we also show that defects introduce localized midgap states with non-zero energies, which should be accessible via scanning probe experiments., Comment: 7 pages, 8 figures