Nodal Superconductors Research Articles

Unconventional self-similar Hofstadter superconductivity from repulsive interactions

Fractal Hofstadter bands have become widely accessible with the advent of moiré superlattices, opening the door to studies of the effect of interactions in these systems. In this work we employ a renormalization group (RG) analysis to demonstrate that the combination of repulsive interactions with the presence of a tunable manifold of Van Hove singularities provides a new mechanism for driving unconventional superconductivity in Hofstadter bands. Specifically, the number of Van Hove singularities at the Fermi energy can be controlled by varying the flux per unit cell and the electronic filling, leading to instabilities toward nodal superconductivity and chiral topological superconductivity with Chern number {{{{{{{mathcal{C}}}}}}}}=pm 6. The latter is characterized by a self-similar fixed trajectory of the RG flow and an emerging self-similarity symmetry of the order parameter. Our results establish Hofstadter quantum materials such as moiré heterostructures as promising platforms for realizing novel reentrant Hofstadter superconductors.

Topological superconductivity in multifold fermion metals

Recently, multifold fermions characterized by band crossings with multifold degeneracy and Fermi surfaces with higher Chern numbers have been discovered experimentally in AlPt (Schroter et al. in Nat Phys 15:759–765, 2019) and XSi(X = Rh,Co) (Sanchez et al. in Nature 567:500–505, 2019; Rao et al. in Nature 567:496–499, 2019; Takane et al. in Phys Rev Lett 122:076402, 2019). In this work, we largely expand the family of multifold fermion materials by pointing out that several well-studied noncentrosymmetric superconductors are indeed multifold fermion metals. Importantly, their normal state topological properties, which have been ignored in previous studies, play an important role in the superconducting properties. Taking mathrm{Li} _{2}mathrm{Pd} _{3}mathrm{B} and mathrm{Li} _{2}mathrm{Pt} _{3}mathrm{B} as examples, we found a large number of unconventional degenerate points, such as double spin-1, spin-3/2, Weyl and double Weyl topological band crossing points near the Fermi energy, which result in finite Chern numbers on Fermi surfaces. Long Fermi arc states in mathrm{Li} _{2}mathrm{Pd} _{3}mathrm{B}, originating from the nontrivial band topology were found. Importantly, it has been shown experimentally that mathrm{Li} _{2}mathrm{Pd} _{3}mathrm{B} and mathrm{Li} _{2}mathrm{Pt} _{3}mathrm{B} are fully gapped and gapless superconductors, respectively. By analyzing the possible pairing symmetries, we suggest that mathrm{Li} _{2}mathrm{Pd} _{3}mathrm{B} is a conventional s-wave superconductor or DIII class topological superconductor with Majorana surface states. mathrm{Li} _{2}mathrm{Pt} _{3}mathrm{B}, being gapless, is likely to be a nodal topological superconductor with dispersionless surface Majorana modes. We further identified that several noncentrosymmetric superconductors, such as mathrm{Mo} _{3}mathrm{Al} _{2}mathrm{C}, PdBiSe, mathrm{Y}_{2}mathrm{C}_{3} and mathrm{La} _{2}mathrm{C}_{3}, are multifold fermion superconductors. This work calls for a revisit for the study of noncentrosymmtric superconductors which provide platforms for investigating the interplay between superconductivity and topologically nontrivial Fermi surfaces.

Nodal and Nematic Superconducting Phases in NbSe_{2} Monolayers from Competing Superconducting Channels.

Transition metal dichalcogenides like 2H-NbSe_{2} in their two-dimensional (2D) form exhibit Ising superconductivity with the quasiparticle spins being firmly pinned in the direction perpendicular to the basal plane. This enables them to withstand exceptionally high magnetic fields beyond the Pauli limit for superconductivity. Using field-angle-resolved magnetoresistance experiments for fields rotated in the basal plane we investigate the field-angle dependence of the upper critical field (H_{c2}), which directly reflects the symmetry of the superconducting order parameter. We observe a sixfold nodal symmetry superposed on a twofold symmetry. This agrees with theoretical predictions of a nodal topological superconducting phase near H_{c2}, together with a nematic superconducting state. We demonstrate that in NbSe_{2} such unconventional superconducting states can arise from the presence of several competing superconducting channels.

Revealing the Thermal Properties of Superconducting Magic-Angle Twisted Bilayer Graphene.

The allegedly unconventional superconducting phase of magic-angle twisted bilayer graphene (MATBG) has been predicted to possess extraordinary thermal properties, as it is formed from a highly diluted electron ensemble with a record-low carrier density (n) of ∼1011 cm-2 and electronic heat capacity (Ce) of <100kB. While these attributes position MATBG as a ground-breaking material platform for revolutionary calorimetric applications, these properties have so far not been experimentally shown. Here, we reveal the thermal properties of superconducting MATBG by monitoring its temperature dependent critical current (Ic) under continuous laser heating at 1550 nm. From the bolometric effect, we extract the temperature dependence of the electronic thermal conductance (Gth), which has a value of Gth = 0.2 pW/K at 35 mK and in the low temperature limit is consistent with a power law dependence, as expected for nodal superconductors. Our work lays the foundation for future thermal transport studies on this system.

Twisted multilayer nodal superconductors

Twisted bilayers of nodal superconductors were recently proposed as a promising platform to host superconducting phases that spontaneously break time-reversal symmetry. Here we extend this analysis to twisted multilayers, focusing on two high-symmetry stackings with alternating ($\pm \theta$) and constant ($\theta$) twist angles. In analogy to alternating-twist multilayer graphene, the former can be mapped to twisted bilayers with renormalized interlayer couplings, along with a remnant gapless monolayer when the number of layers $L$ is odd. In contrast, the latter exhibits physics beyond twisted bilayers, including the occurrence of `magic angles' characterized by cubic band crossings when $L \mod 4 = 3$. Owing to their power-law divergent density of states, such multilayers are highly susceptible to secondary instabilities. Within a BCS mean-field theory, defined in the continuum and on a lattice, we find that both stackings host chiral topological superconductivity in extended regions of their phase diagrams.

Hot-lines topology and the fate of the spin resonance mode in three-dimensional unconventional superconductors

In the quasi-two-dimensional (quasi-2D) copper- and iron-based superconductors, the onset of superconductivity is accompanied by a prominent peak in the magnetic spectrum at momenta close to the wave-vector of the nearby antiferromagnetic state. Such a peak is well described in terms of a spin resonance mode, i.e., a spin-1 exciton theoretically predicted for quasi-2D superconductors with a sign-changing gap. The same theories, however, indicate that such a resonance mode should be absent in a three-dimensional (3D) system with a spherical Fermi surface. This raises the question of the fate of the spin resonance mode in layered unconventional superconductors that are not strongly anisotropic, such as certain heavy-fermion compounds and potentially the newly discovered nickelate superconductor NdNiO$_2$. Here, we use the random-phase-approximation to calculate the dynamical spin susceptibility of 3D superconductors with a $d_{x^2-y^2}$-wave gap symmetry and corrugated cylindrical-like Fermi surfaces. By varying the out-of-plane hopping anisotropy $t_z/t$, we demonstrate that the appearance of a spin resonance mode is determined by the topology of the hot lines -- i.e. lines on the Fermi surface that are connected by the magnetic wave-vector. For an in-plane antiferromagnetic wave-vector, the hot lines undergo a topological transition from open lines to closed loops at a critical $t_z/t$ value. The closed hot lines cross the nodal superconducting lines, making the spin resonance mode overdamped and incoherent. In contrast, for an out-of-plane antiferromagnetic wave-vector, the hot lines remain open and the spin resonance mode remains sharp. We discuss the experimental implications of our results for the out-of-plane dispersion of the spin resonance mode and, more generally, for inelastic neutron scattering experiments on unconventional superconductors.

High-Throughput Investigations of Topological and Nodal Superconductors.

The theory of symmetry indicators has enabled database searches for topological materials in normal conducting phases, which has led to several encyclopedic topological material databases. To date, such a database for topological superconductors is yet to be achieved because of the lack of information about pairing symmetries of realistic materials. In this Letter, sidestepping this issue, we tackle an alternative problem: the predictions of topological and nodal superconductivity in materials for each single-valued representation of point groups. Based on recently developed symmetry indicators for superconductors, we provide comprehensive mappings from pairing symmetries to the topological or nodal superconducting nature for nonmagnetic materials listed in the Inorganic Crystal Structure Database. We quantitatively show that around 90% of computed materials are topological or nodal superconductors when a pairing that belongs to a one-dimensional nontrivial representation of point groups is assumed. When materials are representation-enforced nodal superconductors, positions and shapes of the nodes are also identified. When combined with experiments, our results will help us understand the pairing mechanism and facilitate realizations of the long-sought Majorana fermions promising for topological quantum computations.

Higher-order topological semimetals and nodal superconductors with an order-two crystalline symmetry

Using a systematic relation between topological gapless phases in three dimensions and topological gapped phases in two dimensions, we identify four types of higher-order topological semimetals or nodal superconductors (HOTS), hosting (i) flat zero-energy "Fermi arcs" at crystal hinges, (ii) flat zero-energy hinge arcs coexisting with surface Dirac cones, (iii) chiral or helical hinge modes, or (iv) flat zero-energy hinge arcs connecting nodes only at finite momentum. Bulk-boundary correspondence relates the hinge states to the bulk topology protecting the nodal point or loop. We classify all HOTS for all tenfold-way classes with an order-two crystalline (anti-)symmetry, such as mirror, twofold rotation, or inversion.

Josephson junctions of topological nodal superconductors

Transition metal dichalcogenides (TMDs) offer a unique platform to study unconventional superconductivity, owing to the presence of strong spin-orbit coupling and a remarkable stability to an in-plane magnetic field. A recent study found that when an in-plane field applied to a superconducting monolayer TMD is increased beyond the Pauli critical limit, a quantum phase transition occurs into a topological nodal superconducting phase which hosts Majorana flat bands. We study the current-phase relation of this nodal superconductor in a Josephson junction geometry. We find that the nodal superconductivity is associated with an energy-phase relation that depends on the momentum transverse to the current direction, with a 4\pi4π periodicity in between pairs of nodal points. We interpret this response as a result of a series of quantum phase transitions, driven by the transverse momentum, which separate a trivial phase and two distinct topologically non-trivial phases characterized by different winding invariants. This analysis sheds light on the stability of the Majorana flat bands to symmetry-breaking perturbations.

Evidence for unconventional superconductivity in twisted trilayer graphene.

Magic-angle twisted trilayer graphene (MATTG) has emerged as a moiré material that exhibits strong electronic correlations and unconventional superconductivity1,2. However, local spectroscopic studies of this system are still lacking. Here we perform high-resolution scanning tunnelling microscopy and spectroscopy of MATTG that reveal extensive regions of atomic reconstruction favouring mirror-symmetric stacking. In these regions, we observe symmetry-breaking electronic transitions and doping-dependent band-structure deformations similar to those in magic-angle bilayers, as expected theoretically given the commonality of flat bands3,4. Most notably in a density window spanning two to three holes per moiré unit cell, the spectroscopic signatures of superconductivity are manifest as pronounced dips in the tunnelling conductance at the Fermi level accompanied by coherence peaks that become gradually suppressed at elevated temperatures and magnetic fields. The observed evolution of the conductance with doping is consistent with a gate-tunable transition from a gapped superconductor to a nodal superconductor, which is theoretically compatible with a sharp transition from a Bardeen-Cooper-Schrieffer superconductor to a Bose-Einstein-condensation superconductor with a nodal order parameter. Within this doping window, we also detect peak-dip-hump structures that suggest that superconductivity is driven by strong coupling to bosonic modes of MATTG. Our results will enable further understanding of superconductivity and correlated states in graphene-based moiré structures beyond twisted bilayers5.

Chiral symmetry protected topological nodal superconducting phase and Majorana Fermi arc

With an external in-plane magnetic field, we show the emergence of a topological nodal superconducting phase of the two-dimensional topological surface states. This nodal superconducting phase is protected by the chiral symmetry with a non-zero magnetic field, and there are corresponding Majorana Fermi arcs (also known as flat band Andreev bound states) connecting the two Majorana nodes along the edges, similar to the case of Weyl semimetal. The topological nodal superconductor is an intermediate phase between two different chiral superconductors, and is stable against the effects of substrates. The two-dimensional effective theory of the nodal superconducting phase also captures the low energy behavior of a three-dimensional lattice model which describes the iron-based superconductor with a thin film geometry. The localizations of the Majorana nodes can be manipulated through external in-plane magnetic fields, which may introduce a non-trivial topological Berry phase between them.

A new class of bilayer kagome lattice compounds with Dirac nodal lines and pressure-induced superconductivity

Kagome lattice composed of transition-metal ions provides a great opportunity to explore the intertwining between geometry, electronic orders and band topology. The discovery of multiple competing orders that connect intimately with the underlying topological band structure in nonmagnetic kagome metals AV3Sb5 (A = K, Rb, Cs) further pushes this topic to the quantum frontier. Here we report a new class of vanadium-based compounds with kagome bilayers, namely AV6Sb6 (A = K, Rb, Cs) and V6Sb4, which, together with AV3Sb5, compose a series of kagome compounds with a generic chemical formula (Am-1Sb2m)(V3Sb)n (m = 1, 2; n = 1, 2). Theoretical calculations combined with angle-resolved photoemission measurements reveal that these compounds feature Dirac nodal lines in close vicinity to the Fermi level. Pressure-induced superconductivity in AV6Sb6 further suggests promising emergent phenomena in these materials. The establishment of a new family of layered kagome materials paves the way for designer of fascinating kagome systems with diverse topological nontrivialities and collective ground states.

Majorana flat bands at structured surfaces of nodal noncentrosymmetric superconductors

Surfaces of nodal noncentrosymmetric superconductors can host flat bands of Majorana modes, which provide a promising platform for quantum computation if one can find methods for manipulating localized Majorana wave packets. We study the fate of such flat bands when part of the surface is subjected to an exchange field induced by a ferromagnetic insulator. We use exact diagonalization to find the eigenstates and eigenenergies of the Bogoliubov-de Gennes Hamiltonian of a model system, for which an exchange field is applied along a strip on the surface of a slab. We consider different orientations of the strip and the applied field. If the spin polarization of the field-free system along the field direction is sufficiently large perturbation theory predicts that energies of states which are mostly localized on the exchange-field strip are shifted away from zero energy by an amount proportional to the field strength. On the other hand, energies corresponding to states localized on the field-free strip are only weakly affected by the field. Exact diagonalization confirms this. Moreover, we discuss a setup with a small exchange field applied to the previously field-free strip with the goal of introducing a linear dispersion. By switching this dispersion on and off, a wave packet could be moved in a certain direction. We find that in our model system, a linear dispersion can indeed be achieved. The qualitative features of this dispersion can be predicted from the momentum-dependent spin polarization of the field-free surface.

Euler-obstructed Cooper pairing: Nodal superconductivity and hinge Majorana zero modes

Since the proposal of monopole Cooper pairing in [Phys. Rev. Lett. 120, 067003 (2018)], considerable research efforts have been dedicated to the study of Cooper pairing order parameters constrained (or obstructed) by the nontrivial normal-state band topology at Fermi surfaces in 3D systems. In the current work, we generalize the topologically obstructed pairings between Chern states (like the monopole Cooper pairing) by proposing Euler obstructed Cooper pairing in 3D systems. The Euler obstructed Cooper pairing widely exists between two Fermi surfaces with nontrivial band topology characterized by nonzero Euler numbers; such Fermi surfaces can exist in 3D PT-protected spinless-Dirac/nodal-line semimetals with negligible spin-orbit coupling, where PT is the space-time inversion symmetry. An Euler obstructed pairing channel must have pairing nodes on the pairing-relevant Fermi surfaces, and the total winding number of the pairing nodes is determined by the sum or difference of the Euler numbers on the Fermi surfaces. In particular, we find that when the normal state is time-reversal invariant and the pairing is weak, a sufficiently-dominant Euler obstructed pairing channel with zero total momentum leads to nodal superconductivity. If the Fermi surface splitting is small, the resultant nodal superconductor hosts hinge Majorana zero modes. The possible dominance of the Euler obstructed pairing channel near the superconducting transition and the robustness of the hinge Majorana zero modes against disorder are explicitly demonstrated using effective or tight-binding models. Our work presents the first class of higher-order nodal superconductivity originating from the topologically obstructed Cooper pairing.

Observation of the non-linear Meissner effect

A long-standing theoretical prediction is that in clean, nodal unconventional superconductors the magnetic penetration depth λ, at zero temperature, varies linearly with magnetic field. This non-linear Meissner effect is an equally important manifestation of the nodal state as the well studied linear-in-T dependence of λ, but has never been convincingly experimentally observed. Here we present measurements of the nodal superconductors CeCoIn5 and LaFePO which clearly show this non-linear Meissner effect. We further show how the effect of a small dc magnetic field on λ(T) can be used to distinguish gap nodes from non-nodal deep gap minima. Our measurements of KFe2As2 suggest that this material has such a non-nodal state.

Two-Dimensional Topological Superconductors with Spin-Orbit Coupling and Zeeman Term

Topological superconductors have attracted much attention for their potential applications in realizing topological quantum computing. In this paper, we show a system with 8 × 8 Bogoliubov-de Gennes Hamiltonian. The system has particle-hole symmetry. By adding a Zeeman term to the model, we discuss this system from two situations. In the first case of breaking the inversion symmetry while preserving the mirror symmetry, it is obtained from the topological phase diagram of the system that it is a topological superconductor. In another case where the mirror symmetry is broken while preserving the inversion symmetry, the system has two nodes connected by a flat band of zero-energy Andreev edge states and the topological number is non-zero . It can be concluded that the system is a topological nodal superconductor.

Observation of electronic structure and electron-boson coupling in the low-dimensional superconductor Ta4Pd3Te16

Ternary telluride ${\mathrm{Ta}}_{4}{\mathrm{Pd}}_{3}{\mathrm{Te}}_{16}$ is a layered superconductor with quasi-one-dimensional characteristics. It has attracted great research interests due to the possible nodal superconductivity and its interplay with adjacent charge-density wave as well as proposed nontrivial band topology. Here, we report low-energy electronic band structure of ${\mathrm{Ta}}_{4}{\mathrm{Pd}}_{3}{\mathrm{Te}}_{16}$ by means of high-resolution laser-based angle-resolved photoemission spectroscopy. We acquired the Fermi surface and detailed band dispersions, directly reflecting its multiband nature. Furthermore, dispersion kinks at \ensuremath{\sim}8 meV below the Fermi level are unambiguously observed, indicative of strong electron-boson coupling. We suggest these kinks are probably due to electron-phonon coupling with effective coupling constants quantified to be \ensuremath{\sim}0.72 and \ensuremath{\sim}0.48 for the ${\ensuremath{\gamma}}_{1}$ and ${\ensuremath{\gamma}}_{2}$ bands, respectively. Our results provide valuable information for future studies on the mechanism of superconductivity as well as its intimacy to charge-density wave in ${\mathrm{Ta}}_{4}{\mathrm{Pd}}_{3}{\mathrm{Te}}_{16}$.

Trapping Majorana zero modes in vortices of magnetic texture crystals coupled to nodal superconductors

We propose a mechanism for opening a full bulk energy gap and inducing vortex Majorana zero modes (MZMs) in nodal superconductors (SCs). We show that this becomes possible by coupling the nodal SC of interest to a magnetic texture crystal. The latter consists of superpositions of magnetic textures which repeat periodically in space according to suitable wave vectors that enable spin-flip scattering between all pairs of nodes of the SC, and thus open a full gap in its bulk energy spectrum. In this event, MZMs can be trapped in spin or shift vortices introduced in the magnetic texture crystal. Our approach is generic and applies to nodal SCs of spin-singlet, -triplet, or -mixed type of pairing. Therefore, it promises to find application in a variety of nodal SCs, where the magnetic textures appear either spontaneously due to electron-electron interactions, or are imposed by nanomagnets, or become induced by coupling to a lattice of localized magnetic moments.

Evidence of topological boundary modes with topological nodal-point superconductivity

The extension of the topological classification of band insulators to topological semimetals gave way to the topology classes of Dirac, Weyl, and nodal line semimetals with their unique Fermi arc and drum head boundary modes. Similarly, there are several suggestions to employ the classification of topological superconductors for topological nodal superconductors with Majorana boundary modes. Here, we show that the surface 1H termination of the transition metal dichalcogenide compound 4Hb-TaS$_2$, in which 1T-TaS$_2$ and 1H-TaS$_2$ layers are interleaved, has the phenomenology of a topological nodal point superconductor. We find in scanning tunneling spectroscopy a residual density of states within the superconducting gap. An exponentially decaying bound mode is imaged within the superconducting gap along the boundaries of the exposed 1H layer characteristic of a gapless Majorana edge mode. The anisotropic nature of the localization length of the edge mode aims towards topological nodal superconductivity. A zero-bias conductance peak is further imaged within fairly isotropic vortex cores. All our observations are accommodated by a theoretical model of a two-dimensional nodal Weyl-like superconducting state, which ensues from inter-orbital Cooper pairing. The observation of an intrinsic topological nodal superconductivity in a layered material will pave the way for further studies of Majorana edge modes and its applications in quantum information processing.

Evidence for unconventional superconductivity in twisted bilayer graphene.

The emergence of superconductivity and correlated insulators in magic-angle twisted bilayer graphene (MATBG) has raised the intriguing possibility that its pairing mechanism is distinct from that of conventional superconductors1-4, as described by the Bardeen-Cooper-Schrieffer (BCS) theory. However, recent studies have shown that superconductivity persists even when Coulomb interactions are partially screened5,6. This suggests that pairing in MATBG might be conventional in nature and a consequence of the large density of states of its flat bands. Here we combine tunnelling and Andreev reflection spectroscopy with a scanning tunnelling microscope to observe several key experimental signatures of unconventional superconductivity in MATBG. We show that the tunnelling spectra below the transition temperature Tc are inconsistent with those of a conventional s-wave superconductor, but rather resemble those of a nodal superconductor with an anisotropic pairing mechanism. We observe a large discrepancy between the tunnelling gap ΔT, which far exceeds the mean-field BCS ratio (with 2ΔT/kBTc ~ 25), and the gap ΔAR extracted from Andreev reflection spectroscopy (2ΔAR/kBTc ~ 6). The tunnelling gap persists even when superconductivity is suppressed, indicating its emergence from a pseudogap phase. Moreover, the pseudogap and superconductivity are both absent when MATBG is aligned with hexagonal boron nitride. These findings and other observations reported here provide a preponderance of evidence for a non-BCS mechanism for superconductivity in MATBG.