Two-dimensional time-reversal-invariant topological superconductors host helical Majorana fermions at their boundary. We study the fate of these edge states under the combined influence of strong interactions and disorder, using the effective one-dimensional (1D) lattice model for the edge introduced by Jones and Metlitski [Phys. Rev. B 104, 245130 (2021)]. We specifically develop a strong-disorder renormalization-group analysis of the lattice model and identify a regime in which time-reversal is broken spontaneously, creating random magnetic domains; Majorana fermions localize to domain walls and form an infinite-randomness fixed point, identical to that appearing in the random transverse-field Ising model. While this infinite-randomness fixed point describes a fine-tuned critical point in a purely 1D system, in our edge context there is no obvious time-reversal-preserving perturbation that destabilizes the fixed point. Our analysis thus suggests that the infinite-randomness fixed point emerges as a stable phase on the edge of two-dimensional topological superconductors when strong disorder and interactions are present.
The presence of a boundary enriches the nature of quantum phase transitions. However, the boundary critical phenomena in topological superconductors remain underexplored so far. Here, we investigate the boundary criticality in a two-dimensional correlated time-reversal-invariant topological superconductor tuned through a quantum phase transition into a trivial time-reversal-breaking superconductor. Using sign-problem-free determinant quantum Monte Carlo simulations, we chart the quantum phase diagram and reveal the boundary criticalities encompassing ordinary, special, and extraordinary transitions. Additionally, using renormalization group analysis, we compute the boundary critical exponent up to two loops. Remarkably, the simulations and two-loop renormalization group calculations consistently demonstrate that the presence of the boundary Majorana fermion at the special transition gives rise to a new type of boundary Gross-Neveu-Yukawa fixed point. We conclude with a discussion of possible experimental realizations in iron-based superconductors.
In this paper, we present an exactly solvable model for two-dimensional topological superconductors with helical Majorana edge modes protected by time-reversal symmetry. Our construction is based on the idea of decorated domain walls and makes use of the Kasteleyn orientation on a two-dimensional lattice, which was used for the construction of the symmetry protected fermion phase with
Z2 symmetry by Tarantino et al. and Ware et al. By decorating the time-reversal domain walls with spinful Majorana chains, we are able to construct a commuting projector Hamiltonian with zero correlation length ground state wave function that realizes a strongly interacting version of the two-dimensional topological superconductor. From our construction, it can be seen that the ${T}^{2}=\ensuremath{-}1$ transformation rule for the fermions is crucial for the existence of such a nontrivial phase; with
T2=1, our construction does not work.
In this paper, we discuss strong coupling limits of topological quantum critical points (TQCPs) where quantum phase transitions between two topological distinct superconducting states take place. We illustrate that while superconducting phases on both sides of TQCPs spontaneously break the same symmetries, universality classes of critical states can be identified only when global symmetries in topological states are further specified. In dimensions
d=2and3, we find that continuous
(d+1)th-order transitions at weakly interacting TQCPs that were pointed out previously in the presence of the emergent Lorentz symmetry can be terminated by strongly interacting fixed points of Majorana fields. For
2d time reversal symmetry breaking TQCPs, termination points are supersymmetric with
N=4Nf=1 (where
Nf is the number of four-component Dirac fermions and
N is the number of two-component real fermions) beyond which transitions are discontinuous first-order ones. For
2d time reversal symmetric TQCPs without other global symmetries, termination points of
(d+1)th-order continuous transition lines are generically conformal invariant without supersymmetry. Beyond these strong coupling fixed points, there are first-order discontinuous transitions as far as the protecting symmetry is not spontaneously broken but no direct transitions if the protecting symmetry is spontaneously broken in the presence of strong interactions. In this case, termination points can become supersymmetric with
N=4Nf=2 only when interacting Majorana fermions further display an emergent global
U(1) symmetry. In
3d, strong coupling termination points can be further effectively represented by new emergent gapless real bosons weakly coupled with free gapless Majorana fermions. However, in
1d, time reversal symmetric
(d+1)th continuous transition lines of TQCPs are terminated by simple free Majorana fermion fixed points.
Topological superconductors are believed to host exotic quasiparticle excitations known as Majorana zero modes (MZMs), with much of the evidence based on BCS mean-field theory. The direct application of mean-field arguments is tenuous in finite, isolated systems relevant in some experiments. Here, we develop a new correlation-based method for identifying MZMs in interacting, number-conserving systems. Using the density matrix renormalization group, we study fermion number-conserving models with long-range interactions, which under periodic boundary conditions exhibit robust topological and nontopological superconductivity, tuned by the strength of interaction [Ortiz et al., Phys. Rev. Lett. 113, 267002 (2014)PRLTAO0031-900710.1103/PhysRevLett.113.267002]. We find evidence that, on the topological side, Majorana edge modes appear in open chains, manifesting as the vanishing of the energy splitting between odd- and even-parity ground states with increasing system size. Additionally, off-diagonal two-point correlation functions show nonlocal, parity-dependent edge effects. These correlations reveal the spatial structure of Majorana modes in the many-body wave function. We show that the correlation diagnostic applies broadly, including to short-range interacting models, where topological superconductivity is more fragile due to the absence of a bulk excitation gap.
We study the low-energy physics of a chain of Majorana fermions in the presence of interaction and disorder, emphasizing the difference between Majoranas and conventional (complex) fermions. While in the noninteracting limit both models are equivalent (in particular, belong to the same symmetry class BDI and flow towards the same infinite-randomness critical fixed point), their behavior differs drastically once interaction is added. Our density-matrix renormalization group calculations show that the complex-fermion chain remains at the noninteracting fixed point. On the other hand, the Majorana fermion chain experiences a spontaneous symmetry breaking and localizes for repulsive interaction. To explain the instability of the critical Majorana chain with respect to a combined effect of interaction and disorder, we consider interaction as perturbation to the infinite-randomness fixed point and calculate numerically two-wave-function correlation functions that enter interaction matrix elements. The numerical results supported by analytical arguments exhibit a rich structure of critical eigenstate correlations. This allows us to identify a relevant interaction operator that drives the Majorana chain away from the infinite-randomness fixed point. For the case of complex fermions, the interaction is irrelevant.