We numerically investigate critically delocalized wave functions in models of two-dimensional Dirac fermions, subject to vector potential disorder. These describe the surface states of three-dimensional topological superconductors, and can also be realized through long-range correlated bond randomness in artificial materials like molecular graphene. A ``frozen'' regime can occur for strong disorder in these systems, wherein a single wave function presents a few localized peaks separated by macroscopic distances. Despite this rarefied spatial structure, we find robust correlations between eigenstates at different energies, at both weak and strong disorder. The associated level statistics are always approximately Wigner-Dyson. The system shows generalized Chalker (quantum critical) scaling, even when individual states are quasilocalized in space. We confirm analytical predictions for the density of states and multifractal spectra. For a single Dirac valley, we establish that finite energy states show universal multifractal spectra consistent with the integer quantum Hall plateau transition. A single Dirac fermion at finite energy can therefore behave as a ``Quantum Hall critical metal.'' For the case of two valleys and non-Abelian disorder, we verify predictions of conformal field theory. Our results for the non-Abelian case imply that both delocalization and conformal invariance are topologically protected for multivalley topological superconductor surface states.
Two-dimensional (2D) disordered superconductor (SC) in class D exhibits a disorder-induced quantum multicritical phenomenon among diffusive thermal metal (DTM), topological superconductor (TS), and conventional localized (AI) phases. To characterize the quantum tricritical point where these three phases meet, we carry out a two-loop renormalization group (RG) analysis for 2D Dirac fermion with random mass in terms of the $\ensuremath{\epsilon}$-expansion in the spatial dimension $d=2\ensuremath{-}\ensuremath{\epsilon}$. In two dimensions $(\ensuremath{\epsilon}=0)$, the random mass is marginally irrelevant around a clean-limit fixed point of the gapless Dirac fermion, while there exists an IR unstable fixed point at finite disorder strength that corresponds to the tricritical point. The critical exponent, dynamical exponent, and scaling dimension of the (uniform) mass term are evaluated around the tricritical point by the two-loop RG analysis. Using a mapping between an effective theory for the 2D random-mass Dirac fermion and the
(1+1)-dimensional Gross-Neveu model, we further deduce the four-loop evaluation of the critical exponent, and the scaling dimension of the uniform mass around the tricritical point. Both the two-loop and four-loop results suggest that criticalities of a AI-DTM transition line as well as TS-DTM transition line are controlled by other saddle-point fixed point(s) at finite uniform mass.
We discuss the effects of disorder in time-reversal-invariant topological insulators and superconductors in three spatial dimensions. For three-dimensional topological insulator in the symplectic (AII) symmetry class, the phase diagram in the presence of disorder and a mass term, which drives a transition between trivial and topological insulator phases, is computed numerically by the transfer matrix method. The numerics is supplemented by a field-theory analysis (the large-
Nf expansion, where
Nf is the number of valleys or Dirac cones), from which we obtain the correlation length exponent, and several anomalous dimensions at a nontrivial critical point separating a metallic phase and a Dirac semimetal. A similar-field theory approach is developed for disorder-driven transitions in symmetry classes AIII, CI, and DIII. For these three symmetry classes, where topological superconductors are characterized by integer topological invariant, a complementary description is given in terms of the nonlinear sigma model supplemented with a topological term, which is a three-dimensional analog of the Pruisken term in the integer quantum Hall effect.