Topological gapless phase in Kitaev model on square lattice.

We study the topological feature of gapless states in the fermionic Kitaev model on a square lattice. There are two types of gapless states which are topologically trivial and nontrivial. We show that the topological gapless phase lives in a wide two-dimensional parameter region and are characterized by two vertices of an auxiliary vector field de-fined in the two-dimensional momentum space, with opposite winding numbers. The isolated band touching points, as the topological defects of the field, move, emerge, and disappear as the parameters vary. The band gap starts to open only at the merg-ing points, associated with topologically trivial gapless states. The symmetry protect-ing the topological gapless phase and the robustness under perturbations are also discussed.

Introduction

The concept of topology in condensed matter physics has received great attention[1-4] as it not only predicts new physical phenomena with potential technological applications, but is also closely related to fundamental physics, such as the discovery of fermionic particles and phenomena predicted in high-energy physics, including Majorana[5-10], Dirac[11-17] and Weyl fermions[18-26]. These concepts relate to an interesting topic, topological gapless phase. System in the topological gapless phase exhibits band structures with band-touching points in the momentum space, where these kinds of nodal points appear as topological defects of an auxiliary vector field. Then these points are unremovable due to the symmetry protection, until a pair of them meets and annihilates together. On the other hand, a gapful phase can be topologically non-trivial, commonly referred to as topological insulators and superconductors, the band structure of which is characterized by nontrivial topology.

In this paper, we study the topological gapless phase in the Kitaev model on a square lattice based on analytical solutions. It has been shown that a large class of two-dimensional spinless fermion models exhibit topological superconducting phases[27]. Two different phases are separated by the gapless state as quantum phase boundary. We will show that there are two types of gapless states which are topologically trivial and nontrivial. The topological gapless states are characterized by two vertices with opposite winding numbers in the two-dimensional momentum space. As parameters vary, the isolated topological band touching points, as topological defects, move, emerge, and disappear as the band gap opens. In general, a gapless state only lives at a boundary line, which separates two gapful phases. We show that the topological gapless phase in the present model exist in a wide two-dimensional parameter region. We also analyze the symmetry which protects the topological gapless states the robustness under perturbations. In contrast to previous study, this symmetry does not involve any anti-unitary operation.

Model and phase diagram

We consider the Kitaev model on a square lattice which is employed to depict 2D p-wave superconductors. The Hamiltonian of the tight-binding model on a square lattice takes the following formawhere r is the coordinates of lattice sites and is the fermion annihilation operators at site r. Vectors a = a i, a j, are the lattice vectors in the x and y directions with unitary vectors i and j. The hopping between (pair operator of) neighboring sites is described by the hopping amplitude t (the real order parameter Δ). The last term gives the chemical potential.

Imposing periodic boundary conditions on both directions, the Hamiltonian can be exactly diagonalized. Taking the Fourier transformationwe getwhere the summation of k = (k , k ) is . The Hamiltonian can be diagonalized asby introducing the Bogoliubov operatorwhich satisfies the commutation relations of fermion

Here the spectrum isand the normalized amplitudes arewhere θ k is determined by

The ground state can be constructed aswith the groundstate energywhere is the vacuum state of c k, satisfying for all k. We note that the gapless ground state appears when ε k has a zero point, or band touching point of single γ k-particle spectrum.

In this paper, we are interested in the gapless state arising from the band touching point of the spectrum. The band degenerate point k 0 = (k 0, k 0) is determined by

For different parameters, two bands touch at three types of configurations, single point, double points, and curves in the k -k plane. There are two typical cases:

(i) For Δ ≠ 0, we havein the region . It indicates that there are two nodal points for μ ≠ 0 and |μ/t| ≠ 2. The two points move along the line represented by the equation k 0 = − k 0, and merge at k 0 = (π, −π) when μ/t = ±2. In the case of μ = 0, the nodal points become two nodal lines represented by the equations k 0y = ±π + k 0.

(ii) For Δ = 0, the system becomes interaction free. In this case, Eq. (12) becomesin the regionwhich indicates that the band gap close at the nodal line represented by Eq. (14). In the case of |μ/t| = 0, Eq. (12) represents four linesand

In the case of |μ/t| = 2, the nodal points reduce to a single point.

The phase diagram is illustrated in Fig. 1, depending on the values of μ and Δ (compared with the hopping strength t). We select several points in the phase diagram of Fig. 1, which are representative of all possible typical cases. In addition, we plot the band structures in Fig. 2 for these cases. The configurations of nodal points in the figures are in agreement with our analysis.

Figure 1

Phase diagram of the Kitaev model on a square lattice system on the parameter μ − Δ plane (in units of t). The blue lines indicate the boundary, which separate the gapful phases (yellow) and gapless phase (green). Several points (a–g) at typical positions are indicated and the same letter represents the situations with the similar band structures. The corresponding band structures and the topology of the nodal point in the momentum space are given in Figs 2 and 3, respectively.

Figure 2

Energy spectra from Eq. (7) at points (a–g) marked in the phase diagram in Fig. 1. We see that (i) in the gapless phase region with Δμ ≠ 0 (b,c,e,f), the zero energy points are isolated points; (ii) in the gapless phase region with Δμ = 0 (d,h,i), the zero energy points become lines; (iii) At the phase boundary (b,f), two isolated points merge; and (iv) the zero energy points disappear in the gapful regions (a,g). In the cases of (d,h,i), we only plot the lower bands in order to display the entire nodal lines. In the momentum plane, both k and k range from −π to π.

Figure 3

Isolated zero energy points as topological defects. The planar vector fields defined in Eq. (24) in the momentum space for several typical cases indicated in Fig. (1). It shows that the isolated degeneracy points correspond to the vortices in the momentum space with winding numbers 0 and ±1. The red point represents winding number −1 and the bule one 1 in Fig. (3b). The nontrivial winding number indicates that the isolated nodal point is topologically protected and thus cannot be removed in adiabatic procedures. Two zero energy points (red and blue dots) can merge to a single vortex (black empty circle) with zero winding number. Red lines denote the nodal lines. Abscissa and ordinate represent k and k respectively, ranging from −π to π.

Topological nodal points

In this section, we will show that zero-gap systems are in the topological gapless phase. We demonstrate this point by rewriting the Hamiltonian in the formwhere

The core matrix can be expressed aswhere the components of the auxiliary field (k) = (B , B , B ) are

The Pauli matrices σ are taken as the form

In condensed matter physics, the Dirac or Weyl point acts like a singularity of the Berry curvature in the Brillouin zone, or a magnetic monopole in k-space. When a degenerate point is isolated, it should be a vortex of the vector field B(k), which is the topological defect of the field. We will show that the appearance of gapless states in the present model and their existence regions are fully determined by the topological configurations of the field defects. And these defects can be classified as topologically trivial and non-trivial using topological number, or invariant. The topological invariant of a defect is the winding number

where the unit vector and ∇ = ∂/∂ is the nabla operator in k-space. On the other hand, the winding number of the vortices (the isolated degenerate points) is equivalent to the chirality of Weyl fermion in the context of quantum field theory.

Actually, in the vicinity of the defects, we havewhere q = k − k 0, k 0 = (k 0, k 0) and (k 0, k 0) satisfy Eq. (13), is the momentum in another frame. Around these degenerate points, the Hamiltonian h can be linearized as the formwhich is equivalent to the Hamiltonian for two-dimensional massless relativistic fermions. The corresponding chirality for these particle is defined as

Then we havewhich leads to w = ±1 for two nodal points. The chiral relativistic fermions serve as two-dimensional Weyl fermions. Two Weyl nodes located at two separated degenerate points have opposite chirality. We note that for Δ|μ/t|(|μ/t| − 2) = 0, we have w = 0. At this situation, two Weyl nodes merge at (0, 0) and (±π, ∓π). The topology of the nodal point becomes trivial, and a perturbation hence can open up the energy gap. We illustrate the vortex structure of the degeneracy point in k -k plane in Fig. 3. As shown in figures, we find three types of topological configurations: pair of vortices with opposite chirality, single trivial vortex (or degeneracy lines), and no vortex, corresponding to topological gapless, trivial gapless and gapped phases, respectively.

Symmetry protection of nodal points

In this section, we focus on the symmetry, which protects the nodal point or gapless state. We begin with the translational symmetry of the system. The Hamiltonian is invariant via a translational transformation, i.e., [T , H] = 0 and [T , H] = 0, where T and T are the shift operators defined as

In k-space, the actions of the operators are

Combining with the above two operators we havewhere . In general, for an excited state , which is also the eigenstate of operators , , and , we have

The two-fold zero-energy degenerate states can be constructed aswhere k 0 = |k 0|. It is easy to check that is not the eigenstate of , and

This indicates that the nodal points are protected by the symmetry related to . Furthermore, the zero-energy excited state is a degenerate state, which breaks the symmetry. In contrast to previous studies[28], the symmetry involved in our model does not contain an anti-unitary operator.

Summary

In this paper we have studied the topological gapless state and edge modes of the Kitaev model on a square lattice. The advantage of studying the Kitaev model is that it is the minimal model in two dimensions where one can derive a number of analytical results for the topological gapless phase. It is shown that the topological gapless phase is characterized by two topological vortices with opposite chirality in the momentum plane. These two defects are unremovable until they get together some fixed points. We find that the two topological vortices do not meet when the system parameters drop in a large area of the parameter plane. Then it enhance the possibility to acquire the topological gapless phase in practice. Furthemore, we also analyze the symmetry which protects the topological gapless states the robustness under perturbations. In contrast to previous study, this symmetry does not involve any anti-unitary operation. This is may due to the fact that the magnetic flux is not necessary for the existence of the topological gapless phase in the present model.

In general, a quantum phase is a gapped phase, which exists in a region of parameter space. This region is very robust to parameter variations since the gap cannot close suddenly. Two different quantum phases are separated by gapless phase, which usually lives in the interface of two regions as quantum phase boundary. In this work, we find that the gapless phase also exists in a region of parameter space. And the region is also robust to parameter variations due to the topological nodual points, which cannot be removed suddenly. It indicates that, a transition from a gapped phase to a topological gapless phase should be regarded as a new type of topological quantum phase transition.

Method

The essential nature of the topological boundary is that the band touching points are unavoidable under certain perturbation. We demonstrate this point by consider a perturbations with an extra diagonal hopping term. The perturbed Hamiltonian can be written aswhere t D denotes the diagonal hopping amplitude and vector b = a i + a j. Employing the Fourier transformations in Eq. (2), we still havewhere

The spectrum iswhich only has a shift on t, i.e., μ → μ + t D cos (k  + k ), from the spectrum ε in Eq. (7). Thus the zero point can be obtained directly as following. We are only interested in the non-trivial case with nonzero Δ. In this case, the positions of two vertices are determined by the equation

The existence of the solution (k 0, k 0) requires that

It indicates that t D− term cannot destroy the topological gapless phase but just shifts the region along μ− axis by t D. Figure 4 illustrates this point. Then the gapless state is topologically invariant under the perturbation from the t D− term.

Figure 4

Phase diagram of the Kitaev model with t D− term on a square lattice on the parameter μ − Δ plane (in units of t). The blue lines indicate the original boundary, while the red ones denote the new boundary. It demonstrates that the topological gapless phase is robust under the t D− term perturbation.