Search

We show that spin-singlet $s$-wave multi-band superconductors have a topological phase protected by rotation symmetry and time-reversal symmetry without spin-orbit coupling in two and three dimensions. This topological phase, an example of a representation-protected topological phase, has a $\mathbb{Z}_2$ topological index and is stable as long as the bands at the Fermi energy are formed by a doublet of orbital states with finite angular momenta. In the limit of weak superconducting pair potential, the $\mathbb{Z}_2$ index gives a Fermi-surface formula and is related to the winding number of three-dimensional strong topological superconductors of class CI. We present a model of a topological $s_{\pm}$-wave superconductor that has gapless surface states with a quadratic dispersion and suggest a connection with iron-based superconductors., Comment: 7+10 pages, 2+1 figures

We show that the ground-state expectation value of twisting operator is a topological order parameter for $\text{U}(1)$- and $\mathbb{Z}_{N}$-symmetric symmetry-protected topological (SPT) phases in one-dimensional ``spin'' systems -- it is quantized in the thermodynamic limit and can be used to identify different SPT phases and to diagnose phase transitions among them. We prove that this (non-local) order parameter must take values in $N$-th roots of unity, and its value can be changed by a generalized lattice translation acting as an $N$-ality transformation connecting distinct phases. This result also implies the Lieb-Schultz-Mattis ingappability for SU($N$) spins if we further impose a general translation symmetry. Furthermore, our exact result for the order parameter of SPT phases can predict a large number of LSM ingappabilities by the general lattice translation. We also apply the $N$-ality property to provide an efficient way to construct possible multi-critical phase transitions starting from a single Hamiltonian with a unique gapped ground state., Comment: 6 pages

We discuss the topology of the parameter space of invertible phases with an onsite symmetry $G$, i.e., quantum many-body ground states that have neither fractionalization nor spontaneous breaking of the symmetry. The classification of invertible phases is known to be obtained by counting the connected components in the parameter space of the invertible phases. We consider its generalization -- the deformation classes of the mappings from $n$-dimensional spheres $S^n$ to this parameter space for arbitrary integer $n$. We argue a direct one-to-one correspondence in the framework of lattice models between the non-contractibility of $S^n$ and (i) the classification of invertible phases in $d$ dimensions when $d\geq n$; or (ii) zero-dimensional invertible Hamiltonians parametrized by $S^{n-d}$ when $d

We propose from first-principles calculations that spin-orbit-free materials in the $RE_8\mathrm{Co}X_3$ group ($RE$ = rare earth elements, $X$ = Al, Ga, or In) are ideal spinless Dirac semimetals whose Fermi surfaces are fourfold degenerate band-crossing points (without including spin degeneracy). Despite the lack of space inversion symmetry in these materials, Dirac points are formed on the rotation-symmetry axis due to accidental degeneracies of two bands corresponding to different 2-dimensional irreducible representations of $C_{6v}$ group. The surface states have two midgap bands emanating from the projection in the surface Brillouin zone of the bulk Dirac points, in accordance with nontrivial Zak phases for each glide sector. We also investigate, through first-principles calculations and effective model analysis, various phase transitions caused by lattice distortion or elemental substitutions from the Dirac semimetal phase to distinct topological semimetallic phases such as nonmagnetic linked-nodal-line and Weyl semimetals (characterized by the second Stiefel-Whitney class) and ferromagnetic Weyl semimetals.

Lieb-Schultz-Mattis (LSM) theorems impose non-perturbative constraints on the zero-temperature phase diagrams of quantum lattice Hamiltonians (always assumed to be local in this paper). LSM theorems have recently been interpreted as the lattice counterparts to mixed 't Hooft anomalies in quantum field theories that arise from a combination of crystalline and global internal symmetry groups. Accordingly, LSM theorems have been reinterpreted as LSM anomalies. In this work, we provide a systematic diagnostic for LSM anomalies in one spatial dimension. We show that gauging subgroups of the global internal symmetry group of a quantum lattice model obeying an LSM anomaly delivers a dual quantum lattice Hamiltonian such that its internal and crystalline symmetries mix non-trivially through a group extension. This mixing of crystalline and internal symmetries after gauging is a direct consequence of the LSM anomaly, i.e., it can be used as a diagnostic of an LSM anomaly. We exemplify this procedure for a quantum spin-1/2 chain obeying an LSM anomaly resulting from combining a global internal $\mathbb{Z}^{\,}_{2}\times\mathbb{Z}^{\,}_{2}$ symmetry with translation or reflection symmetry. We establish a triality of models by gauging a $\mathbb{Z}^{\,}_{2}\subset\mathbb{Z}^{\,}_{2}\times\mathbb{Z}^{\,}_{2}$ symmetry in two ways, one of which amounts to performing a Kramers-Wannier duality, while the other implements a Jordan-Wigner duality. We discuss the mapping of the phase diagram of the quantum spin-1/2 $XYZ$ chains under such a triality. We show that the deconfined quantum critical transitions between Neel and dimer orders are mapped to either topological or conventional Landau-Ginzburg transitions. Finally, we extend our results to $\mathbb{Z}^{\,}_{n}$ clock models and provide a reinterpretation of the dual internal symmetries in terms of $\mathbb{Z}^{\,}_{n}$ charge and dipole symmetries., Comment: 88 pages, 6 figures

Vortex lines in superconducting Dirac semimetals realize crystal-symmetry-protected gapless vortex-line phases in which gapless excitations propagate inside a vortex line, in the presence of appropriate crystal symmetry, spin-orbit coupling, and multi-band structures. Here we present a general scheme to classify possible gapless vortex-line phases in superconducting Dirac semimetals with rotation (or screw) symmetry and inversion symmetry, assuming that the rotation (screw) axis is parallel to the vortex line. The rotation (screw)-symmetry-protected gapless modes are stable as long as they have different rotation (screw) eigenvalues. The underlying mechanism for the formation of gapless vortex bound states depends on irreducible representations of rotation (screw) symmetry subject to a vortex field and is classified into three types: (i) accidental band crossing of two vortex bound-state modes under rotation symmetry; (ii) accidental and (iii) enforced band crossing of four vortex bound-state modes under screw symmetry. We present a tight-binding model of screw-symmetry-protected Dirac semimetal with an $s$-wave pair potential, demonstrating a gapless vortex-line phase of type (ii). We obtain four gapless modes of vortex bound states whose gapless points (Majorana zero modes) pinned at a time-reversal invariant momentum (TRIM) when the Fermi energy is close to the Dirac points. As the Fermi energy is moved away from the Dirac points, the four gapless modes are split into a pair of two gapless modes with vanishing excitation energy at non-TRIMs. In closing, we discuss Nb$_3$Pt as a candidate material with the four-fold screw-symmetry-protected Dirac cones that can host a gapless vortex-line phase., Comment: 19 pages, 6 figures

Three-dimensional (3D) gapped topological phases with fractional excitations are divided into two subclasses: one has topological order with point-like and loop-like excitations fully mobile in the 3D space, and the other has fracton order with point-like excitations constrained in lower-dimensional subspaces. These exotic phases are often studied by exactly solvable Hamiltonians made of commuting projectors, which, however, are not capable of describing those exhibiting surface states with gapless chiral dispersion. Here we introduce a systematic way, based on cellular construction recently proposed for 3D topological phases, to construct another type of exactly solvable models in terms of coupled quantum wires with given inputs of cellular structure, two-dimensional Abelian topological order, and their gapped interfaces. We show that our models can describe both 3D topological and fracton orders (and even their hybrid) and study their universal properties such as quasiparticle statistics and topological ground-state degeneracy. We also apply this construction to two-dimensional coupled-wire models with ordinary topological orders and translation-symmetry-enriched topological orders. Our results pave the way for effective quantum field theory descriptions or microscopic model realizations of fracton orders with chiral gapless surface states., Comment: 58+64 pages, 40+1 figures, 1+12 tables, v3: References added. Published version

We propose an index $\mathcal{I}_G$ which characterizes the degree of ingappability, namely the difficulty to induce a unique ground state with a nonvanishing excitation gap, in the presence of a symmetry $G$. $\mathcal{I}_G$ represents the dimension of the subspace of ambient uniquely-gapped in the entire $G$-invariant "theory space". The celebrated Lieb-Schultz-Mattis theorem corresponds, in our formulation, to the case $\mathcal{I}_G=0$ (completely ingappable) for the symmetry $G$ including the lattice translation symmetry. We illustrate the usefulness of the index by discussing the phase diagram of spin-$1/2$ antiferromagnets in various dimensions, which do not necessarily have the translation symmetry., Comment: 3 figures

We propose a geometric {approach to Lieb-Schultz-Mattis theorem for} quantum many-body systems with discrete spin-rotation symmetries and lattice inversion or rotation symmetry, but without translation symmetry assumed. Under symmetry-twisting on a $(d-1)$-dimensional plane, we find that any $d$-dimensional inversion-symmetric spin system possesses a doubly degenerate spectrum when it hosts a half-integer spin at the inversion-symmetric point. We also show that any rotation-symmetric generalized spin model with a projective representation at the rotation center has a similar degeneracy under symmetry-twisting. We argue that these degeneracies imply that {a unique symmetric gapped ground state that is smoothly connected to product states} is forbidden in the original untwisted systems -- generalized inversional/rotational Lieb-Schultz-Mattis theorems without lattice translation symmetry imposed. The traditional Lieb-Schultz-Mattis theorems with translations also fit in the proposed framework., Comment: 6 pages, 7 figures

The fractional Josephson effect is known to be a characteristic phenomenon of topological Josephson junctions hosting Majorana zero modes (MZMs), where the Josephson current has a $4\pi$ (rather than a $2\pi$) periodicity in the phase difference between the two topological superconductors. We introduce a one-dimensional model of a topological superconductor/normal-metal/superconductor (SNS) junction with the normal-metal (N) region of finite length, which is intermediate regime between the short- and long-junction limits. Assuming weak tunneling at the SN interfaces, we investigate resonance and finite-size effects on the fractional Josephson effect due to the existence of several discrete energy levels in the N region in which wavefunctions have oscillating nodal structure. Through careful analysis of the sign change in the transmission amplitudes through the junction and the fermion parity of the two MZMs, we find that the fractional Josephson current is proportional to the parity of total fermion numbers including both filled normal levels and two MZMs. Furthermore, we elucidate drastic enhancement of the Josephson current due to the resonance between a discrete level in the N region and MZMs., Comment: 15 pages, 6 figures; published version

We present a series of models of three-dimensional rotation-symmetric fragile topological insulators in class AI (time-reversal symmetric and spin-orbit-free systems), which have gapless surface states protected by time-reversal ($T$) and $n$-fold rotation ($C_n$) symmetries ($n=2,4,6$). Our models are generalizations of Fu's model of a spinless topological crystalline insulator, in which orbital degrees of freedom play the role of pseudo-spins. We consider minimal surface Hamiltonian with $C_n$ symmetry in class AI and discuss possible symmetry-protected gapless surface states, i.e., a quadratic band touching and multiple Dirac cones with linear dispersion. We characterize topological structure of bulk wave functions in terms of two kinds of topological invariants obtained from Wilson loops: $\mathbb{Z}_2$ invariants protected by $C_n$ ($n=4,6$) and time-reversal symmetries, and $C_2T$-symmetry-protected $\mathbb{Z}$ invariants (the Euler class) when the number of occupied bands is two. Accordingly, our models realize two kinds of fragile topological insulators. One is a fragile $\mathbb{Z}$ topological insulator whose only nontrivial topological index is the Euler class that specifies the number of surface Dirac cones. The other is a fragile $\mathbb{Z}_2$ topological insulator having gapless surface states with either a quadratic band touching or four (six) Dirac cones, which are protected by time-reversal and $C_4$ ($C_6$) symmetries. Finally, we discuss the instability of gapless surface states against the addition of $s$-orbital bands and demonstrate that surface states are gapped out through hybridization with surface-localized $s$-orbital bands., Comment: 20 pages, 7 figures

We study the stability of topological crystalline superconductors in the symmetry class DIIIR and in two-dimensional space when perturbed by quartic contact interactions. It is known that no less than eight copies of helical pairs of Majorana edge modes can be gapped out by an appropriate interaction without spontaneously breaking any one of the protecting symmetries. Hence, the noninteracting classification $\mathbb{Z}$ reduces to $\mathbb{Z}^{\,}_{8}$ when these interactions are present. It is also known that the stability when there are less than eight modes can be understood in terms of the presence of topological obstructions in the low-energy bosonic effective theories, which prevent opening of a gap. Here, we investigate the stability of the edge theories with four, two, and one edge modes, respectively. We give an analytical derivation of the topological term for the first case, because of which the edge theory remains gapless. For two edge modes, we employ bosonization methods to derive an effective bosonic action. When gapped, this bosonic theory is necessarily associated to the spontaneous symmetry breaking of either one of time-reversal or reflection symmetry whenever translation symmetry remains on the boundary. For one edge mode, stability is explicitly established in the Majorana representation of the edge theory., Comment: 16 pages with 1 figure

We formulate a $\mathbb{Z}_k$-parafermionization/bosonization scheme for one-dimensional lattice models and field theories on a torus, starting from a generalized Jordan-Wigner transformation on a lattice, which extends the Majorana-Ising duality at $k=2$. The $\mathbb{Z}_k$-parafermionization enables us to investigate the critical theories of parafermionic chains whose fundamental degrees of freedom are parafermionic, and we find that their criticality cannot be described by any existing conformal field theory. The modular transformations of these parafermionic low-energy critical theories as general consistency conditions are found to be unconventional in that their partition functions on a torus transform differently from any conformal field theory when $k>2$. Explicit forms of partition functions are obtained by the developed parafermionization for a large class of critical $\mathbb{Z}_k$-parafermionic chains, whose operator contents are intrinsically distinct from any bosonic or fermionic model in terms of conformal spins and statistics. We also use the parafermionization to exhaust all the $\mathbb{Z}_k$-parafermionic minimal models, complementing earlier works on fermionic cases., Comment: 27 pages, 3 figures

The interplay of time-reversal and $n$-fold rotation symmetries ($n=2,4,6$) is known to bring a new class of topological crystalline insulators (TCIs) having $n$ surface Dirac cones due to surface rotation anomaly. We show that the proximity-induced $s$-wave superconductivity on the surface of these TCIs yields a topological superconducting phase in which two Majorana zero modes are bound to a vortex, and that $n$-fold rotation symmetry ($n=2,4,6$) enriches the topological classification of a superconducting vortex from $\mathbb{Z}_2$ to $\mathbb{Z}_2\times\mathbb{Z}_2$. Using a model of a three-dimensional high-spin topological insulator with $s$-wave superconductivity and two-fold rotation symmetry, we show that, with increasing chemical potential, the number of Majorana zero modes at one end of a vortex changes as $2\to1\to0$ through two topological vortex phase transitions. In addition, we show that additional magnetic-mirror symmetry further enhances the topological classification to $\mathbb{Z} \times \mathbb{Z}$, Comment: 10 pages, 4 figures

We propose and investigate a simple one-dimensional model for a single-channel quantum wire hosting electrons that interact repulsively and are subject to a significant spin-orbit interaction. We show that an external Zeeman magnetic field, applied at the right angle to the Rashba spin-orbit axis, drives the wire into a correlated spin-density wave state with gapped spin and gapless charge excitations. By computing the ground-state degeneracies of the model with either (anti-)periodic or open boundary conditions, we conclude that the correlated spin-density state realizes a gapless symmetry-protected topological phase, as the ground state is unique in the ring geometry while it is two-fold degenerate in the wire with open boundaries. Microscopically the two-fold degeneracy is found to be protected by the conservation of the magnetization parity. Open boundaries induce localized zero-energy (midgap) states which are described, at the special Luther-Emery point of the model, by Majorana fermions. We find that spin densities at the open ends of the wire exhibit unusual long-ranged correlations despite the fact that all correlations in the bulk of the wire decay in a power-law or exponential fashion. Our study exposes the crucial importance of the long-ranged string operator needed to implement the correct commutation relations between spin densities at different points in the wire. Along the way we rederive the low-energy theory of Galilean-invariant electron systems in terms of current operators.

The phase diagram of the quantum spin-1/2 antiferromagnetic $J^{\,}_{1}$-$J^{\,}_{2}$ XXZ chain was obtained by Haldane using bosonization techniques. It supports three distinct phases for $0\leq J^{\,}_{2}/J^{\,}_{1}<1/2$, i.e., a gapless algebraic spin liquid phase, a gapped long-range ordered Neel phase, and a gapped long-range ordered dimer phase. Even though the Neel and dimer phases are not related hierarchically by a pattern of symmetry breaking, it was shown that they meet along a line of quantum critical points with a U(1) symmetry and central charge $c=1$. Here, we extend the analysis made by Haldane on the quantum spin-1/2 antiferromagnetic $J^{\,}_{1}$-$J^{\,}_{2}$ XYZ chain using both bosonization and numerical techniques. We show that there are three Neel phases and the dimer phase that are separated from each other by six planes of phase boundaries realizing U(1) criticality when $0\leq J^{\,}_{2}/J^{\,}_{1}<1/2$. We also show that each long-range ordered phase harbors topological point defects (domain walls) that are dual to those across the phase boundary in that a defect in one ordered phase locally binds the other type of order around its core. By using the bosonization approach, we identify the critical theory that describes simultaneous proliferation of these dual point defects, and show that it supports an emergent U(1) symmetry that originates from the discrete symmetries of the XYZ model. To confirm this numerically, we perform DMRG calculation and show that the critical theory is characterized by the central charge $c=1$ with critical exponents that are consistent with those obtained from the bosonization approach. Furthermore, we generalize the field theoretic description of direct continuous phase transition to higher dimensions, especially in $d=3$, by using a non-linear sigma model (NLSM) with a topological term., Comment: 25 pages with 14 figures

We establish the correspondence between the fractional charge bound to a vortex in a textured lattice and the relevant bulk band topology in two-dimensional (2D) topological crystalline insulators. As a representative example, we consider the Kekule textured graphene whose bulk band topology is characterized by a 2D $\mathbb{Z}_{2}$ topological invariant $\nu_{\rm 2D}$ protected by inversion symmetry. The fractional charge localized at a vortex in the Kekule texture is shown to be related to the change in the bulk topological invariant $\nu_{\rm 2D}$ around the vortex, as in the case of the Su-Schriefer-Heeger model in which the fractional charge localized at a domain wall is related to the change in the bulk charge polarization between degenerate ground states. We show that the effective three-dimensional (3D) Hamiltonian, where the angle $\theta$ around a vortex in Kekule-textured graphene is a third coordinate, describes a 3D axion insulator with a quantized magnetoelectric polarization. The spectral flow during the adiabatic variation of $\theta$ corresponds to the chiral hinge modes of an axion insulator and determines the accumulated charge localized at the vortex, which is half-quantized when chiral symmetry exists. When chiral symmetry is absent, electric charge localized at the vortex is no longer quantized, but the vortex always carries a half-quantized Wannier charge as long as inversion symmetry exists. For the cases when magnetoelectric polarization is quantized due to the presence of symmetry that reverses the space-time orientation, we classify all possible topological crystalline insulators whose vortex defect carries a fractional charge., Comment: 5+8 pages, 3+6 figures, 1+1 tables

We propose a systematic approach to constructing microscopic models with fractional excitations in three-dimensional (3D) space. Building blocks are quantum wires described by the (1+1)-dimensional conformal field theory (CFT) associated with a current algebra $\mathfrak{g}$. The wires are coupled with each other to form a 3D network through the current-current interactions of $\mathfrak{g}_1$ and $\mathfrak{g}_2$ CFTs that are related to the $\mathfrak{g}$ CFT by a nontrivial conformal embedding $\mathfrak{g} \supset \mathfrak{g}_1 \times \mathfrak{g}_2$. The resulting model can be viewed as a layer construction of a 3D topologically ordered state, in which the conformal embedding in each wire implements the anyon condensation between adjacent layers. Local operators acting on the ground state create point-like or loop-like deconfined excitations depending on the branching rule. We demonstrate our construction for a simple solvable model based on the conformal embedding $SU(2)_1 \times SU(2)_1 \supset U(1)_4 \times U(1)_4$. We show that the model possesses extensively degenerate ground states on a torus with deconfined quasiparticles, and that appropriate local perturbations lift the degeneracy and yield a 3D $Z_2$ gauge theory with a fermionic $Z_2$ charge., Comment: 6+10 pages, 3 figures. v2: References added, published version

The coupled-wire construction provides a useful way to obtain microscopic Hamiltonians for various two-dimensional topological phases, among which fractional quantum Hall states are paradigmatic examples. Using the recently introduced flux attachment and vortex duality transformations for coupled wires, we show that this construction is remarkably versatile to encapsulate phenomenologies of hierarchical quantum Hall states: the Jain-type hierarchy states of composite fermions filling Landau levels and the Haldane-Halperin hierarchy states of quasiparticle condensation. The particle-hole conjugate transformation for coupled-wire models is also given as a special case of the hierarchy construction. We also propose coupled-wire models for the composite Fermi liquid, which turn out to be compatible with a sort of the particle-hole symmetry implemented in a nonlocal way at $\nu=1/2$. Furthermore, our approach shows explicitly the connection between the Moore-Read Pfaffian state and a chiral $p$-wave pairing of the composite fermions. This composite fermion picture is also generalized to a family of the Pfaffian state, including the anti-Pfaffian state and Bonderson-Slingerland hierarchy states., Comment: 38 pages, 11 figures, v2: published version

Correlation functions of dimer operators, the product operators of spins on two adjacent sites, are studied in the spin-$\frac{1}{2}$ XXZ chain in the critical regime. The amplitudes of the leading oscillating terms in the dimer correlation functions are determined with high accuracy as functions of the exchange anisotropy parameter and the external magnetic field, through the combined use of bosonization and density-matrix renormalization group methods. In particular, for the antiferromagnetic Heisenberg model with SU(2) symmetry, logarithmic corrections to the dimer correlations due to the marginally-irrelevant operator are studied, and the asymptotic form of the dimer correlation function is obtained. The asymptotic form of the spin-Peierls excitation gap including logarithmic corrections is also derived., Comment: 11 pages, 8 figures, 1 table + Supplemental Material (7 pages, 6 tables). v2: accepted version, numerical data for finite magnetic field (Figs. 7, 8 and Tables in Supplemental Material) are modified slightly, some references are added

In three-dimensional noncentrosymmetric materials two-fold screw rotation symmetry forces electron's energy bands to have Weyl points at which two bands touch. This is illustrated for space groups No. 19 ($P2_12_12_1$) and No. 198 ($P2_13$), which have three orthogonal screw rotation axes. In the case of space groups No. 61 ($Pbca$) and No. 205 ($P$a-3) that have extra inversion symmetry, Weyl points are promoted to four-fold degenerate line nodes in glide-invariant planes. The three-fold rotation symmetry present in the space groups No. 198 and No. 205 allows Weyl and Dirac points, respectively, to appear along its rotation axes in the Brillouin zone and generates four-fold and six-fold degeneracy at the $\Gamma$ point and R point, respectively., Comment: 8 pages, v2 added references; v3 corrected typos, sec.IIB slightly expanded

Topological semimetals have energy bands near the Fermi energy sticking together at isolated points/lines/planes in the momentum space, which are often accompanied by stable surface states and intriguing bulk topological responses. Although it has been known that certain crystalline symmetries play an important role in protecting band degeneracy, a general recipe for stabilizing the degeneracy, especially in the presence of spin-orbit coupling, is still lacking. Here we show that a class of novel topological semimetals with point/line nodes can emerge in the presence of an off-centered rotation/mirror symmetry whose symmetry line/plane is displaced from the center of other symmorphic symmetries in nonsymmorphic crystals. Due to the partial translation perpendicular to the rotation axis/mirror plane, an off-centered rotation/mirror symmetry always forces two energy bands to stick together and form a doublet pair in the relevant invariant line/plane in momentum space. Such a doublet pair provides a basic building block for emerging topological semimetals with point/line nodes in systems with strong spin-orbit coupling., Comment: 37 pages, 9 figures, 2 tables

We study two-dimensional bosonic symmetry protected topological (SPT) phases which are protected by reflection symmetry and local symmetry [$Z_N\rtimes R$, $Z_N\times R$, U(1)$\rtimes R$, or U(1)$\times R$], in the search for two-dimensional bosonic analogs of topological crystalline insulators in integer-$S$ spin systems with reflection and spin-rotation symmetries. To classify them, we employ a Chern-Simons approach and examine the stability of edge states against perturbations that preserve the assumed symmetries. We find that SPT phases protected by $Z_N\rtimes R$ symmetry are classified as $\mathbb{Z}_2\times\mathbb{Z}_2$ for even $N$ and 0 (no SPT phase) for odd $N$ while those protected by U(1)$\rtimes R$ symmetry are $\mathbb{Z}_2$. We point out that the two-dimensional Affleck-Kennedy-Lieb-Tasaki state of $S=2$ spins on the square lattice is a $\mathbb{Z}_2$ SPT phase protected by reflection and $\pi$-rotation symmetries., Comment: 11+epsilon pages, 2 figures; v3: typos corrected

In general, the stability of a band crossing point indicates the presence of a quantized topological number associated with it. In particular, the recent discovery of three-dimensional Dirac semimetals in Na$_{3}$Bi and Cd$_{3}$As$_{2}$ demonstrates that a Dirac point with four-fold degeneracy can be stable as long as certain crystalline symmetries are supplemented in addition to the time-reversal and inversion symmetries. However, the topological charges associated with Na$_{3}$Bi and Cd$_{3}$As$_{2}$ are not clarified yet. In this work, we identify the topological charge of three-dimensional Dirac points. It is found that although the simultaneous presence of the time-reversal and inversion symmetries forces the net chiral charge to vanish, a Dirac point can carry another quantized topological charge when an additional rotation symmetry is considered. Two different classes of Dirac semimetals are identified depending on the nature of the rotation symmetries. First, the conventional symmorphic rotational symmetry which commutes with the inversion gives rise to the class I Dirac semimetals having a pair of Dirac points on the rotation axes. Since the topological charges of each pair of Dirac points have the opposite sign, a pair-creation or a pair-annihilation is required to change the number of Dirac points in the momentum space. On the other hand, the class II Dirac semimetals possess a single isolated Dirac point at a time-reversal invariant momentum, which is protected by a screw rotation. The non-symmorphic nature of screw rotations allows the anti-commutation relation between the rotation and inversion symmetries, which enables to circumvent the doubling of the number of Dirac points and create a single Dirac point at the Brillouin zone boundary., Comment: 23 pages, 13 figures, 4 tables

We study interaction effects on the topological crystalline insulators protected by time-reversal ($T$) and reflection symmetry ($R$) in two and three spatial dimensions. From the stability analysis of the edge states with bosonization, we find that the classification of the two-dimensional SPT phases protected by $Z_2\times[\mbox{U(1)}\rtimes T]$ symmetry is reduced from $\mathbb{Z}$ to $\mathbb{Z}_4$ by interactions, where the $Z_2$ symmetry denotes the reflection whose mirror plane is the two-dimensional plane itself. By extending the approach recently proposed by Isobe and Fu, we show that the classification of the three-dimensional SPT phases (i.e., topological crystalline insulators) protected by $R\times[\mbox{U(1)}\rtimes T]$ symmetry is reduced from $\mathbb{Z}$ to $\mathbb{Z}_8$ by interactions., Comment: 9 pages, 1 figures, v2: typos corrected

The conditions for both the stability and the breakdown of the topological classification of gapped ground states of noninteracting fermions, the tenfold way, in the presence of quartic fermion-fermion interactions are given for any dimension of space. This is achieved by encoding the effects of interactions on the boundary gapless modes in terms of boundary dynamical masses. Breakdown of the noninteracting topological classification occurs when the quantum nonlinear sigma models for the boundary dynamical masses favor quantum disordered phases. For the tenfold way, we find that (i) the noninteracting topological classification $\mathbb{Z}^{\,}_{2}$ is always stable, (ii) the noninteracting topological classification $\mathbb{Z}$ in even dimensions is always stable, (iii) the noninteracting topological classification $\mathbb{Z}$ in odd dimensions is unstable and reduces to $\mathbb{Z}^{\,}_{N}$ that can be identified explicitly for any dimension and any defining symmetries. We also apply our method to the three-dimensional topological crystalline insulator SnTe from the symmetry class AII$+R$, for which we establish the reduction $\mathbb{Z}\to\mathbb{Z}^{\,}_{8}$ of the noninteracting topological classification., Comment: 29 pages

We propose that the topological magneto-electric (ME) effect, a hallmark of topological insulators (TIs), can be realized in thin films of TIs in the $\nu=0$ quantum Hall state under magnetic field or by doping two magnetic ions with opposite signs of exchange coupling. These setups have the advantage compared to previously proposed setups that a uniform configuration of magnetic field or magnetization is sufficient for the realization of the topological ME effect. To verify our proposal, we numerically calculate ME response of TI thin films in the cylinder geometry and that of effective 2D models of surface Dirac fermions. The ME response is shown to converge to the quantized value corresponding to the axion angle $\theta=\pm \pi$ in the limit of the large top and bottom surface area of TI films, where non-topological contributions from the bulk and the side surface are negligible., Comment: 8 pages, 6 figures

We construct the generic phase diagrams encoding the topologically distinct localized and delocalized phases of noninteracting fermionic quasiparticles for any symmetry class from the tenfold way in one, two, and three dimensions. To this end, we start from a massive Dirac Hamiltonian perturbed by a generic disorder for any dimension of space and for any one of the ten symmetry classes from the tenfold way. The physics of Anderson localization is then encoded by a two-dimensional phase diagram that we deduce from the topology of the space of normalized Dirac masses. This approach agrees with previously known results and gives an alternative explanation for the even-odd effect in the one-dimensional chiral symmetry classes. We also give a qualitative explanation for the Gade singularity and Griffiths effects in the density of states using the first homotopy group of the normalized Dirac masses in two dimensions. Finally, this approach is used to analyze the stability of massless Dirac fermions on the surface of three-dimensional topological crystalline insulators., Comment: 45 pages, 6 figures

Three-dimensional topological insulators of finite thickness can show the quantum Hall effect (QHE) at the filling factor $\nu=0$ under an external magnetic field if there is a finite potential difference between the top and bottom surfaces. We calculate energy spectra of surface Weyl fermions in the $\nu=0$ QHE and find that gapped edge states with helical spin structure are formed from Weyl fermions on the side surfaces under certain conditions. These edge channels account for the nonlocal charge transport in the $\nu=0$ QHE which is observed in a recent experiment on (Bi$_{1-x}$Sb$_x$)$_2$Te$_3$ films. The edge channels also support spin transport due to the spin-momentum locking. We propose an experimental setup to observe various spintronics functions such as spin transport and spin conversion., Comment: 7 pages, 5 figures

We study $\mathbb{Z}_3$ symmetry-protected topological (SPT) phases in one-dimensional spin systems with $Z_3 \times Z_3$ symmetry. We construct ground-state wave functions of the matrix product form for nontrivial $\mathbb{Z}_3$ phases and their parent Hamiltonian from a cocycle of the group cohomology $H^2(Z_3\times Z_3,U(1))$. The Hamiltonian is an SU(3) version of the Affleck-Kennedy-Lieb-Tasaki (AKLT) model, consisting of bilinear and biquadratic terms of su(3) generators in the adjoint representation. A generalization to the SU($N$) case, the SU($N$) AKLT Hamiltonian, is also presented which realizes nontrivial $\mathbb{Z}_N$ SPT phases. We use the infinite-size variant of the density matrix renormalization group (iDMRG) method to determine the ground-state phase diagram of the SU(3) bilinear-biquadratic model as a function of the parameter $\theta$ controlling the ratio of the bilinear and biquadratic coupling constants. The nontrivial $\mathbb{Z}_3$ SPT phase is found for a range of the parameter $\theta$ including the point of vanishing biquadratic term ($\theta=0$) as well as the SU(3) AKLT point [$\theta=\arctan(2/9)$]. A continuous phase transition to the SU(3) dimer phase takes place at $\theta \approx -0.027\pi$, with a central charge $c\approx3.2$. For SU(3) symmetric cases we define string order parameters for the $\mathbb{Z}_3$ SPT phases in a similar way to the conventional Haldane phase. We propose simple spin models that effectively realize the SU(3) and SU(4) AKLT models., Comment: 23 pages, 6 figures

We study the stability of gap-closing (Weyl or Dirac) points in the three-dimensional Brillouin zone of semimetals using Clifford algebras and their representation theory. We show that a pair of Weyl points with $\mathbb{Z}_2$ topological charge are stable in a semimetal with time-reversal and reflection symmetries when the square of the product of the two symmetry transformations equals minus identity. We present toy models of $\mathbb{Z}_2$ Weyl semimetals which have surface modes forming helical Fermi arcs. We also show that Dirac points with $\mathbb{Z}_2$ topological charge are stable in a semimetal with time-reversal, inversion, and SU(2) spin rotation symmetries when the square of the product of time-reversal and inversion equals plus identity. Furthermore, we briefly discuss the topological stability of point nodes in superconductors using Clifford algebras., Comment: 14 pages, 2 figures

We study the stability against disorder of surface states of weak $\mathbb{Z}_2$ topological insulators (superconductors) which are stacks of strong $\mathbb{Z}_2$ topological insulators (superconductors), considering representative Dirac Hamiltonians in the Altland-Zirnbauer symmetry classes in various spatial dimensions. We show that, in the absence of disorder, surface Dirac fermions of weak $\mathbb{Z}_2$ topological insulators (superconductors) can be gapped out by a Dirac mass term which couples surface Dirac cones and leads to breaking of a translation symmetry (dimerization). The dimerization mass is a unique Dirac mass term in the surface Dirac Hamiltonian, and the two dimerized gapped phases which differ in the sign of the Dirac mass are distinguished by a $\mathbb{Z}_2$ index. In other words the dimerized surfaces can be regarded as a strong $\mathbb{Z}_2$ topological insulator (superconductor). We argue that the surface states are not localized by disorder when the ensemble average of the Dirac mass term vanishes., Comment: 8 pages

A two-dimensional spin-directed $\mathbb{Z}^{\,}_{2}$ network model is constructed that describes the combined effects of dimerization and disorder for the surface states of a weak three-dimensional $\mathbb{Z}^{\,}_{2}$ topological insulator. The network model consists of helical edge states of two-dimensional layers of $\mathbb{Z}^{\,}_{2}$ topological insulators which are coupled by time-reversal symmetric interlayer tunneling. It is argued that, without dimerization of interlayer couplings, the network model has no insulating phase for any disorder strength. However, a sufficiently strong dimerization induces a transition from a metallic phase to an insulating phase. The critical exponent $\nu$ for the diverging localization length at metal-insulator transition points is obtained by finite-size scaling analysis of numerical data from simulations of this network model. It is shown that the phase transition belongs to the two-dimensional symplectic universality class of Anderson transition., Comment: 36 pages and 27 figures, plus Supplemental Material

We classify topological insulators and superconductors in the presence of additional symmetries such as reflection or mirror symmetries. For each member of the 10 Altland-Zirnbauer symmetry classes, we have a Clifford algebra defined by operators of the generic (time-reversal, particle-hole, or chiral) symmetries and additional symmetries, together with gamma matrices in Dirac Hamiltonians representing topological insulators and superconductors. Following Kitaev's approach, we classify gapped phases of non-interacting fermions under additional symmetries by examining all possible distinct Dirac mass terms which can be added to the set of generators of the Clifford algebra. We find that imposing additional symmetries in effect changes symmetry classes and causes shifts in the periodic table of topological insulators and superconductors. Our results are in agreement with the classification under reflection symmetry recently reported by Chiu et al. Several examples are discussed including a topological crystalline insulator with mirror Chern numbers and mirror superconductors., Comment: 18 pages

We study the ground-state properties of one-dimensional mixtures of bosonic and fermionic atoms resonantly coupled to fermionic Feshbach molecules. When the particle densities of fermionic atoms and Feshbach molecules differ, the system undergoes various depletion transitions between binary and ternary mixtures, as a function of the detuning parameter. However, when the particle densities of fermionic atoms and Feshbach molecules are identical, the molecular conversion and disassociation processes induce a gap in a sector of low-energy excitations, and the remaining system can be described by a two-component Tomonaga-Luttinger liquid. Using a bosonization scheme, we derive the effective low-energy Hamiltonian for the system, which has a similar form as that of the two-chain problem of coupled Tomonaga-Luttinger liquids. With the help of improved perturbative renormalization group analysis of the latter problem, we determine the ground-state phase diagram and find that it contains a phase dominated by composite superfluid or pairing correlations between the open and closed resonant channels., Comment: 14 pages, 4 figures, published version

We discuss the thermal (or gravitational) responses in topological superconductors and in topological phases in general. Such thermal responses (as well as electromagnetic responses for conserved charge) provide a definition of topological insulators and superconductors beyond the single-particle picture. In two-dimensional topological phases, the Str\v{e}da formula for the electric Hall conductivity is generalized to the thermal Hall conductivity. Applying this formula to the Majorana surface states of three-dimensional topological superconductors predicts cross-correlated responses between the angular momentum and thermal polarization (entropy polarization). We also discuss a use of D-branes in string theory as a systematic tool to derive all such topological terms and topological responses. In particular, we relate the $\mathbb{Z}_2$ index of topological insulators introduced by Kane and Mele (and its generalization to other symmetry classes and to arbitrary dimensions) to the K-theory charge of non-BPS D-branes, and vice versa. We thus establish a link between the stability of non-BPS D-branes and the topological stability of topological insulators., Comment: 16 pages, 2 figures; Submitted to a topical issue of the Comptes Rendus de l Academie des Sciences (CRAS)

The ground-state phase diagram of a spin-1/2 XXZ chain with competing ferromagnetic nearest-neighbor (J_1<0) and antiferromagnetic second-neighbor (J_2>0) exchange couplings is studied by means of the infinite time evolving block decimation algorithm and effective field theories. For the SU(2)-symmetric (Heisenberg) case, we show that the nonmagnetic phase in the range -4

Anderson localization is studied for two flavors of massless Dirac fermions in 2D space perturbed by static disorder that is invariant under a chiral symmetry (chS) and a time-reversal symmetry (TRS) operation which, when squared, is equal either to plus or minus the identity. The former TRS (symmetry class BDI) can for example be realized when the Dirac fermions emerge from spinless fermions hopping on a 2D lattice with a linear energy dispersion such as the honeycomb lattice (graphene) or the square lattice with $\pi$-flux per plaquette. The latter TRS is realized by the surface states of 3D $\mathbb{Z}_{2}$-topological band insulators in symmetry class CII. In the phase diagram parametrized by the disorder strengths, there is an infrared stable line of critical points for both symmetry classes BDI and CII. Here we discuss a "global phase diagram" in which disordered Dirac fermion systems in all three chiral symmetry classes, AIII, CII, and BDI, occur in 4 quadrants, sharing one corner which represents the clean Dirac fermion limit. This phase diagram also includes symmetry classes AII [e.g., appearing at the surface of a disordered 3D $\mathbb{Z}_2$-topological band insulator in the spin-orbit (symplectic) symmetry class] and D (e.g., the random bond Ising model in two dimensions) as boundaries separating regions of the phase diagram belonging to the three chS classes AIII, BDI, and CII. Moreover, we argue that physics of Anderson localization in the CII phase can be presented in terms of a non-linear-sigma model (NLsM) with a $\mathbb{Z}_{2}$-topological term. We thereby complete the derivation of topological or Wess-Zumino-Novikov-Witten terms in the NLsM description of disordered fermionic models in all 10 symmetry classes relevant to Anderson localization in two spatial dimensions., Comment: 28 pages

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 Z_2 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 Z_2 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 \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 Z_2 trivial and non-trivial phases., Comment: 17 pages, 9 figures

We study nontrivial responses of topological superconductors and superfluids to the temperature gradient and rotation of the system. In two-dimensional gapped systems, the Str\v{e}da formula for the electric Hall conductivity is generalized to the thermal Hall conductivity. Applying this formula to the Majorana surface states of three-dimensional topological superconductors predicts cross-correlated responses between the orbital angular momentum and thermal polarization (entropy polarization). These results can be naturally related to the gravitoelectromagnetism description of three-dimensional topological superconductors and superfluids, analogous to the topological magnetoelectric effect in ${\mathbb{Z}}_2$ topological insulators., Comment: 7 pages, 1 figure

We summarize our theoretical findings on the ground-state phase diagram of the spin-1/2 XXZ chain having competing nearest-neighbor (J1) and antiferromagnetic next-nearest-neighbor (J2) couplings. Our study is mainly concerned with the case of ferromagnetic J1, and the case of antiferromagnetic J1 is briefly reviewed for comparison. The phase diagram contains a rich variety of phases in the plane of J1/J2 versus the XXZ anisotropy Delta: vector-chiral phases, N\'eel phases, several dimer phases, and Tomonaga-Luttinger liquid phases. We discuss the vector-chiral order that appears for a remarkably wide parameter space, successive N\'eel-dimer phase transitions, and an emergent nonlocal string order in a narrow region of ferromagnetic J1 side., Comment: 8 pages (1 column), 3 figures, published version

We study theoretically a strongly type-II s-wave superconducting state of two-dimensional Dirac fermions in proximity to a ferromagnet having in-plane magnetization. It is shown that a magnetic domain wall can host a chain of equally spaced vortices in the superconducting order parameter, each of which binds a Majorana fermion state. The overlap integral of neighboring Majorana states is sensitive to the position of the chemical potential of the Dirac fermions. This leads to a characteristic V-shaped dependence of thermal conductivity of Majorana fermions on the chemical potential., Comment: 4+ pages, 2 figures

We show that Majorana edge modes of two-dimensional spin-triplet topological superconductors have Ising-like spin density whose direction is determined by the d-vector characterizing the spin-triplet pairing symmetry. Exchange coupling between an impurity spin (S=1/2) and Majorana edge modes is thus Ising-type. Under external magnetic field perpendicular to the Ising axis, the system can be mapped to a two-level system with Ohmic dissipation, which is equivalent to the anisotropic Kondo model. The magnetic response of the impurity spin can serve as a local experimental probe for the order parameter., Comment: 4+ pages, 2 figures

We study the ground-state phase diagram of a one-dimensional spin-1/2 easy-plane XXZ model with a ferromagnetic nearest-neighbor (NN) coupling $J_1$ and a competing next-nearest-neighbor (NNN) antiferromagnetic coupling $J_2$ in the parameter range $0

Many classes of non-linear sigma models (NLSMs) are known to contain composite operators with an arbitrary number 2s of derivatives ("high-gradient operators") which appear to become strongly relevant within RG calculations at one (or fixed higher) loop order, when the number 2s of derivatives becomes large. This occurs at many conventional fixed points of NLSMs which are perturbatively accessible within the usual epsilon-expansion in d=2+\epsilon dimensions. Since such operators are not prohibited from occurring in the action, they appear to threaten the very existence of such fixed points. At the same time, for NLSMs describing metal-insulator transitions of Anderson localization in electronic conductors, the strong RG-relevance of these operators has been previously related to statistical properties of the conductance of samples of large finite size ("conductance fluctuations"). In this paper, we analyze this question, not for perturbative RG treatments of NLSMs, but for 2d Wess-Zumino-Witten (WZW) models at level k, perturbatively in the current-current interaction of the Noether current. WZW models are special ("Principal Chiral") NLSMs on a Lie Group G with a WZW term at level k. In these models the role of high-gradient operators is played by homogeneous polynomials of order 2s in the Noether currents, whose scaling dimensions we analyze. For the Lie Supergroup G=GL(2N|2N) and k=1, this corresponds to time-reversal invariant problems of Anderson localization in the so-called chiral symmetry classes, and the strength of the current-current interaction, a measure of the strength of disorder, is known to be completely marginal (for any k). We find that all high-gradient (polynomial) operators are, to one loop order, irrelevant or relevant depending on the sign of that interaction., Comment: 22 pages

The quantum spin Hall effect shares many similarities (and some important differences) with the quantum Hall effect for the electric charge. As with the quantum (electric charge) Hall effect, there exists a correspondence between bulk and boundary physics that allows to characterize the quantum spin Hall effect in diverse and complementary ways. In this paper, we derive from the network model that encodes the quantum spin Hall effect, the so-called Z_2 network model, a Dirac Hamiltonian in two dimensions. In the clean limit of this Dirac Hamiltonian, we show that the bulk Kane-Mele Z_2 invariant is nothing but the SU(2) Wilson loop constructed from the SU(2) Berry connection of the occupied Dirac-Bloch single-particle states. In the presence of disorder, the non-linear sigma model (NLSM) that is derived from this Dirac Hamiltonian describes a metal-insulator transition in the standard two-dimensional symplectic universality class. In particular, we show that the fermion doubling prevents the presence of a topological term in the NLSM that would change the universality class of the ordinary two-dimensional symplectic metal-insulator transition. This analytical result is fully consistent with our previous numerical studies of the bulk critical exponents at the metal-insulator transition encoded by the Z_2 network model. Finally, we improve the quality and extend the numerical study of boundary multifractality in the Z_2 topological insulator. We show that the hypothesis of two-dimensional conformal invariance at the metal-insulator transition is verified within the accuracy of our numerical results., Comment: 31 pages, 9 figures, submitted to New J Phys

It has recently been shown that in every spatial dimension there exist precisely five distinct classes of topological insulators or superconductors. Within a given class, the different topological sectors can be distinguished, depending on the case, by a Z or a Z_2 topological invariant. This is an exhaustive classification. Here we construct representatives of topological insulators and superconductors for all five classes and in arbitrary spatial dimension d, in terms of Dirac Hamiltonians. Using these representatives we demonstrate how topological insulators (superconductors) in different dimensions and different classes can be related via dimensional reduction by compactifying one or more spatial dimensions (in Kaluza-Klein-like fashion). For Z-topological insulators (superconductors) this proceeds by descending by one dimension at a time into a different class. The Z_2-topological insulators (superconductors), on the other hand, are shown to be lower-dimensional descendants of parent Z-topological insulators in the same class, from which they inherit their topological properties. The 8-fold periodicity in dimension d that exists for topological insulators (superconductors) with Hamiltonians satisfying at least one reality condition (arising from time-reversal or charge-conjugation/particle-hole symmetries) is a reflection of the 8-fold periodicity of the spinor representations of the orthogonal groups SO(N) (a form of Bott periodicity). We derive a relation between the topological invariant that characterizes topological insulators/superconductors with chiral symmetry and the Chern-Simons invariant: it relates the invariant to the electric polarization (d=1), or to the magnetoelectric polarizability (d=3). Finally, we discuss topological field theories describing the space time theory of linear responses, and study how the presence of inversion symmetry modifies the classification., Comment: 66 pages, 5 figures, updated version of manuscript to appear in New Journal of Physics (copy of proofs submitted to Journal on May 12, 2010)

We generalize universal relations between the multifractal exponent \alpha_0 for the scaling of the typical wave function magnitude at a (Anderson) localization-delocalization transition in two dimensions and the corresponding critical finite size scaling (FSS) amplitude \Lambda_c of the typical localization length in quasi-one-dimensional (Q1D) geometry: (i) When open boundary conditions are imposed in the transverse direction of Q1D samples (strip geometry), we show that the corresponding critical FSS amplitude \Lambda_c^o is universally related to the boundary multifractal exponent \alpha_0^s for the typical wave function amplitude along a straight boundary (surface). (ii) We further propose a generalization of these universal relations to those symmetry classes whose density of states vanishes at the transition. (iii) We verify our generalized relations [Eqs. (6) and (7)] numerically for the following four types of two-dimensional Anderson transitions: (a) the metal-to-(ordinary insulator) transition in the spin-orbit (symplectic) symmetry class, (b) the metal-to-(Z_2 topological insulator) transition which is also in the spin-orbit (symplectic) class, (c) the integer quantum Hall plateau transition, and (d) the spin quantum Hall plateau transition., Comment: 13 pages, 7 figures