Two-Dimensional Topological Superconductors with Spin-Orbit Coupling and Zeeman Term ()
1. Introduction
Topological superconductors have attracted a lot of attention in the past few years, characterized by having gap in the bulk and topologically protected gapless edge states composed of Majorana fermions [1] - [9]. Interest in topological superconductors is closely related to the exotic Majorana fermions, which are self-Hermitian particles and have potential applications in quantum computing [6] [7] [10] [11] [12] [13] [14]. Edge states generated from p-wave superconductors can prevent decoherence topologically [11] [15] [16]. Fu and Kane proposed that topological insulators proximitized by s-wave superconductors that could lead to Majorana zero modes at vortices [17]. Although two-dimensional and three-dimensional topological insulators have been recently discovered [18] [19], finding suitable topological insulators for these experiments remains a challenge. Adding Zeeman field to semiconductors with Rashba spin orbit coupling and inducing s-wave pairing through proximity effect is another method to obtain chiral p-wave superconductors [20] [21] [22] [23] [24]. Majorana zero modes have been reported iron-based superconductors [25], which indicates that there is topological superconductivity in iron-based superconductors. One of the motivations for this paper is that iron-based superconductors have a property that can withstand large in-plane magnetic fields [26] [27].
The remainder of this paper is organized as follows. In Section 2, we describe our efficient model 8 × 8 Bogoliubov-de Gennes Hamiltonian and perform a symmetry analysis on it. Results and discussion are given in Section 3. In Section 3.1, under the transformation U, when
, the transformed Hamiltonian is block diagonalized, and each block is still particle-hole symmetric. Adding a Zeeman term to the original Hamiltonian, when the breaks the inversion symmetry while preserving the mirror symmetry, we plot a topological phase diagram about
plane. We can see from this topological phase diagram that the system is a topological superconductor. In Section 3.2, adding a Zeeman term to the original Hamiltonian, when breaking the mirror symmetry while preserving the inversion symmetry, we plot the band structures of the lattice model and the density of states distribution under the same parameters. From the band structure, it can be seen that the two nodes are connected by a flat band of zero-energy Andreev edge states. The topological number is non-zero
. Therefore, the system in this case is a topological nodal superconductor. A brief summary is given in Section 4.
2. Model
Under the Nambu basis, the electron operator
, where 1 and 2 denote the top and bottom surface states and
and
represent spin-up and spin-down states, respectively. In momentum space, this system is described by an 8 × 8
Bogoliubov-de Gennes Hamiltonian
, with
(1)
where
is the Fermi velocity to be set to unity,
, where
,
and
are Pauli matrices acting on the spin, surface and particle-hole space, respectively. t describes the tunneling effect between the top and bottom surface states, its specific form is
.
is the chemical potential.
is the spin singlet intra-surface pairing strength.
is the spin triplet pairing strength and
is the spin triplet pairing strength.
is the Zeeman term. The matrix form of the Hamiltonian (1) has a special symmetry, which is the particle-hole symmetry,
with
, where
is the complex conjugation operator. For the Hamiltonian (1) has inversion symmetry
with
. When we choose the Zeeman term
, the Hamiltonian (1) is time-reversal invariant. We define a time-reversal operator
that satisfies
. The Hamiltonian (1) also has mirror symmetry
, with
, this is mirror symmetry about the xy plane.
3. Results and Discussion
3.1. Breaking Inversion
While Preserving Mirror
We have
, with the transformation matrix
(2)
We can transform the model Hamiltonian (1) into the following, then
(3)
with the basis
. When
, we can turn the above Hamiltonian into block diagonal
with
(4)
(5)
The terms that break inversion while preserving mirror is
, adding a term
with a constant c to the original Hamiltonian (1),
then
(6)
For
, the transformed Hamiltonian is block diagonalized, and each block is still particle-hole symmetric with
. The physical meaning of the c term is that, on the untransformed basis,
(7)
which, in the real space, is an on-site x-Zeeman term but in opposite directions in the two layers, namely the staggered x-Zeeman term. Make the following substitution for the basis
, where
and
. While for the transformed Hamiltonian, this term for
reads
(8)
which is a spin-singlet s-wave pairing. In the presence of both spin-singlet s-wave pairing and spin-triplet p-wave pairing, the inversion is no longer a conserved quantity. The result for
is similar, one only needs to substitute
and
. We show the Topological phase diagram in the
plane in Figure 1 and Figure 2. As can be seen in Figure 1(c) in the middle-up and -down regime, the Chern number is a nontrivial zero, then it can also be written as
, while in the centered middle regime, Chern number is a trivial zero. As show in Figure 2(c) in the middle regime, the Chern number is a nontrivial zero, namely
.
When an x-Zeeman term
is added to the model Hamiltonian (1), x-Zeeman term breaks inversion symmetry while protecting mirror symmetry. And it can be known from Figure 1 and Figure 2 that this system is a topological superconductor at this time. Next, we discuss the case of adding a Zeeman term that preserves inversion symmetry while breaking mirror symmetry.
3.2. Breaking Mirror
While Preserving Inversion
The terms that break mirror while preserve inversion is
. Adding the x-Zeeman term
to the original Hamiltonian,
Figure 1. (color online) Topological phase diagrams of (a)
, (b)
, (c)
in
plane with
,
,
,
and
.
Figure 2. (color online) Topological phase diagrams of (a)
, (b)
(c)
in
plane with
,
,
,
and
.
under the transformation U,
(9)
which is not block diagonalized for
. Although the addition of the x-Zeeman term
to the Hamiltonian (1) breaks the time-reversal symmetry of the system, Hamiltonian respects an effective time-reversal symmetry
, where
. Combining time-reversal symmetry and particle-hole symmetry yields the chiral symmetry
, with
. As a result, Hamiltonian belongs to the class BDI [1]. For a superconducting system with chiral symmetry, an explicit expression of the topological number is given by [28],
(10)
with
denotes the momentum with periodic boundary condition and
the one with open boundary condition in the numerics. For the parameters chosen in Figure 3, namely,
,
,
,
,
,
and
, then
and
, which is consistent with Figure 3(a). In Figure 3(a), we plot the band structures of the lattice model. From Figure 3(a), the two nodes are connected by a flat band of zero-energy Andreev edge states. Figure 3(b) shows the wave function distribution for the flat band of the zero-energy Andreev edge states. Therefore, when breaking the mirror symmetry while preserving the inversion symmetry, the system is a topological nodal superconductor.
4. Conclusion
In summary, we show a system with 8 × 8 Bogoliubov-de Gennes Hamiltonian and add a Zeeman term to the system. Under the transformation U, the system is still particle-hole symmetry. When breaking the inversion symmetry while
Figure 3. (color online) The band structures of h with 50-sites along the y-direction with
,
,
,
,
,
and
, We also plot the spacial distributions of lowest positive energy states at
in the right-hand side accordingly.
preserving the mirror symmetry, according to the topological phase diagram of the system in the
plane in Figure 1 and Figure 2, the system is a topological superconductor. When breaking the mirror symmetry while preserving the inversion symmetry, the system has chiral symmetry. The two nodes are connected by a flat band of zero-energy Andreev edge states (Figure 3(a)) and the topological number is non-zero
. Therefore, the system is a topological nodal superconductor.