A time-reversal invariant Kitaev-type model is introduced in which spins (Dirac matrices) on the square lattice interact via anisotropic nearest-neighbor and next-nearest-neighbor exchange interactions. The model is exactly solved by mapping it onto a tight-binding model of free Majorana fermions coupled with static
Z2 gauge fields. The Majorana fermion model can be viewed as a model of time-reversal-invariant superconductor and is classified as a member of symmetry class DIII in the Altland-Zirnbauer classification. The ground-state phase diagram has two topologically distinct gapped phases which are distinguished by a
Z2 topological invariant. The topologically nontrivial phase supports both a Kramers' pair of gapless Majorana edge modes at the boundary and a Kramers' pair of zero-energy Majorana states bound to a 0-flux vortex in the $\ensuremath{\pi}$-flux background. Power-law decaying correlation functions of spins along the edge are obtained by taking the gapless Majorana edge modes into account. The model is also defined on the one-dimension ladder, in which case again the ground-state phase diagram has
Z2 trivial and nontrivial phases.
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.
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 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 2D topological superconductors when strong disorder and interactions are present.
The two-dimensional Majorana-Hubbard model on the
π-flux lattice involves the most local interaction possible for a lattice with one relativistic Majorana fermion per site. It is predicted to host a Majorana semimetal, dimerized phases, and quantum-critical points with an emergent
U(1) symmetry. Nonperturbative lattice calculations for this model remain a major open challenge, thus far limited to zero hopping. We present quantum Monte Carlo simulations of the model at finite hopping using Hamiltonian fermion bags, which are based on a continuous-time expansion with Boltzmann weights that are Pfaffians of a matrix with dimension set by the expansion order, and which give the sign of the weight. Sign problem severity is related to the location of zero modes in the free hopping matrix, and for couplings
g≲0.7 and
β=L/4 it is fully removed for lattices up to at least
256 sites. The interaction strength
g>0 enhances the
U(1)-charged dimerization doublet and suppresses the neutral channel, as the renormalization group predicts. Increasing the interactions also restores
U(1) symmetry in the component of the dimerization correlator that is even under bond orientation exchange, but not its odd component. There is a persistent anisotropy in the odd component with the symmetry of a lattice nematic. The data is consistent with a gapless semimetal throughout the sign-problem-minimal window.
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.
A fundamental open issue in physics is whether and how the fermion sign problem in quantum Monte Carlo (QMC) simulations can be solved generically. Here, we show that Majorana-time-reversal (MTR) symmetries can provide a unifying principle to solve the fermion sign problem in interacting fermionic models. By systematically classifying Majorana-bilinear operators according to the anticommuting MTR symmetries they respect, we rigorously prove that there are two and only two fundamental symmetry classes which are sign-problem-free and which we call the "Majorana class" and "Kramers class," respectively. Novel sign-problem-free models in the Majorana class include interacting topological superconductors and interacting models of charge-4e superconductors. We believe that our MTR unifying principle could shed new light on sign-problem-free QMC simulation on strongly correlated systems and interacting topological matters.