1 Introduction
Topological phases of matter have been the central issues of condensed matter study in the past decades [
1–
4], which is recently extended to higher-order band topology [
5]. Initially addressed as quantized electric multiple insulators [
6], higher-order topological insulators (HOTIs) are described by new topological invariants [
7,
8] and have different bulk-boundary correspondence. In 2-dimensional (2D) systems, the second-order topological phases are characterized by 0D in-gap corner states, while a 3D HOTI may host 1D gapless hinge states or in-gap corner states [
5]. Up to now, HOTIs have been reported in various systems [
9–
17], including quasicrystals [
18–
20] and Anderson insulators [
21,
22].
Honeycomb lattice materials, such as graphene and MoS
2, are ideal platforms for exploring topological phases [
2,
23–
28]. The famous Kane–Mele model of quantum spin Hall effect is first theoretically proposed on graphehe [
29], which is later demonstrated to have an extremely small band gap [
30]. Honeycomb lattice materials supporting the Kane–Mele model are continuously found, for instance, silicene [
31–
34], germanium [
32], and stanene [
35–
37], etc. First-principles calculation suggested jacutingaite Pt
2HgSe
3 as a new Kane–Mele quantum spin Hall insulator [
38–
40], where a 0.1 eV band gap is achieved in experiment very recently [
41]. A modified Kane–Mele model with in-plane Zeeman field is shown to be second-order topological insulator [
42], where transport induced dimer state is formed during the resonant tunneling of two topological corner states [
43]. The topological corner states are embedded in a large energy gap and localized in real-space obtuse corners, which makes them hard to detect via transport means.
To increase the probability of observing these corner states, a possible solution is creating appropriate locations to host more corner states in one single system. In this work, based on the modified Kane–Mele model [
42,
43], we propose a multihollow structure through precise engineering of 2D honeycomb lattice. The existence of multiple corner states is revealed in the energy spectrum as well as eigenfunction distributions, and their transport features are numerically investigated in a two-probe setup. Both the transmission function and partial local density of states confirm that, multiple corner states contribute to resonant tunneling and result in global resonant states in a large range of system parameters. We suggest a singleelectron source via adiabatic pumping of multiple corner states, and the pumped charge is found to be well quantized.
The rest of the paper is organized as follows. In Section 2, the model Hamiltonian for topological corner states in honeycomb lattice is introduced, and quantum transport formalism is briefly reviewed. Section 3 includes numerical results and relevant discussion on transport features of the multihollow structure. In addition, a singleelectron source based on this multihollow system is proposed. A summary is finally given in Section 4.
2 Model and formalism
We start with the modified Kane–Mele model in the presence of an in-plane Zeeman field, where second-order topological corner states reside in real-space obtuse intersections of the honeycomb lattice [
42]. The tight-binding Hamiltonian of this model is written as [
42]
Here
is the creation operator at site
i for both spin up (↑) and down (↓).
t in the first term is the nearest-neighbor hopping energy, which is taken as the energy unit. The second term denotes intrinsic spin–orbit coupling (SOC) with strength
tSO, which is represented by next-nearest-neighbor hopping. The first and the second terms together give the conventional Kane–Mele model for quantum spin Hall effect.
νij = ±1 if the electron takes a clockwise or anticlockwise turn from site
j to
i. The in-plane Zeeman field is introduced in the third term, with
B = (
Bx, By, 0) and
λ the Zeeman field strength. This Zeeman field breaks the time-reversal symmetry and hence destroys quantum spin Hall effect. For a certain field strength, such as
λ ≥ 0.2 [
42], a pair of zero-energy corner states appear at the obtuse intersections of two different zigzag boundaries in a diamond-shaped honeycomb lattice, as shown in Fig. 1(a), which are the second-order topological phase protected by the mirror symmetry.
To study the transport features of topological corner states, the diamond-shaped region hosting corner states is connected to two semi-infinite leads. The lead is chosen to be free of magnetic field, and its Hamiltonian follows the conventional Kane–Mele model:
where is the creation/annihilation operator in this lead and n/m labels its lattice site.
Influence of the
α lead is treated as the self-energy,
, which is conveniently evaluated through the transfer matrix approach [
44,
45]. As one of the most important transport properties, the transmission function of an electron with energy
E from lead
β to lead
α is expressed in terms of the nonequilibrium Green’s functions as [
46]
where
Gr/
Ga is the retarded/advanced Green’s function of the central region,
and
Ga =
Gr†. The linewidth function of lead
α is defined as
, and its eigenvector form can effectively speed up the calculation [
47].
The partial local density of states (PLDOS) [
48–
50] is another useful tool in transport study, which describes the real-space distribution of a scattering state:
“Partial” means only electrons incident from lead α contribute to dosα, and i represents the real-space lattice site. Transport properties including transmission and PLDOS of multiple corner states are numerically investigated and thoroughly discussed in the following section.
3 Numerical results and discussion
In the calculation, we set
tSO = 0.1
t, and the Zeeman field is along y direction, i.e.,
B = (0,
By, 0). Without loss of generality, we choose a typical diamond-shaped flake with length and width
Lx =
Ly = 30
a, where
a is the lattice constant. After diagonalizing the Hamiltonian of this isolated system, the eigenfunction distribution in Fig. 1(a) clearly addresses the existence of second-order topological corner states, whose energy level highlighted in blue is isolated in a large gap. Such topological corner states have fractional charge distribution, i.e.,
e/2 charge at each corner, which is independent of the Zeeman field orientation [
42].
Next we introduce hollows into the diamond-shaped flake fulfilling the requirement of topological corner states. The hollow is also diamond-shaped and has obtuse interior intersections. In Fig. 1(b), a single hollow structure is analyzed. It is found that two pairs of zero-energy states arise in the energy gap, whose eigenfunctions are localized at corners displayed in blue and red, respectively. The interior corner states in red color are induced by this single hollow, and its energy level deviates slightly from
E = 0. Note that the energy levels of in-gap states are symmetric about
E = 0, and we only show the
E > 0 states in the following to highlight the eigenfunction distributions. When the isolated system has two identical hollows, three pairs of corner states are shown in Fig. 1(c), where the two red states are degenerate in energy and real spaces. For the four-hollow structure exhibited in Fig. 1(d), five near zero-energy states appear in the large gap, with distinguishable energy difference between each other. These numerical evidence proves that: by introducing hollows into the honeycomb lattice, multiple corner states emerge in the energy spectrum as well as real space intersections, which greatly enlarges the probability of detecting them. The multihollow structure of honeycomb lattice can be realized through fine structure engineering. For instance, the interior edges of graphene can be exactly engineered as armchair, zigzag, or hosting periodic line defects [
51–
53]. Patterned graphene structures such as triangles and hexagons with zigzag boundaries have been fabricated via hydrogen-assisted anisotropic etching [
54,
55], which provides a solid experimental foundation for the proposed multihollow structure.
We take the four-hollow structure as the example, and introduce two transport setups for the detection and utilization of multiple corner states.
3.1 Global resonant state in bulk insulator
To study the transport features of multiple corner states, we build a two-probe system by connecting two narrow leads to the four-hollow structure, as shown in Fig. 2(a). In this case, all corner states remain intact in the diamond region. For instance, the eigenfunction distribution of zero-energy states locates at the exterior corners of this isolated central region. In Fig. 2(b), transmission spectrums with respect to the Zeeman field strength are presented. When
λ = 0, quantum spin Hall edge states provide two perfect conducting channels, leading to
T = 2. In the presence of the Zeeman field, these edge states are gradually destroyed and topological corner states form, which causes the transmission quickly drops to
T = 0. With the increasing of
λ, multiple
T = 1 transmission peaks emerge for energies ranging from
E = 10
−13 to
E = 10
−2. We demonstrate in the following that these peaks correspond to the resonant tunneling of corner states. Compared with the transport of one or two corner states [
43], the energy scope of detecting corner states is greatly enlarged in this four-hollow system.
Partial local density of states (PLDOS) can effectively reveal the dynamic property of resonant tunneling process. We select one transmission curve from Fig. 2(b), and reproduce it in Fig. 3(a). As shown in red arrows, five typical Zeeman field strengths are chosen to observe their PLDOS, which are listed in Figs. 3(b)–(f). For a weak field λ1, the transmission is approximately T = 1, and Fig. 3(b) shows that the transport is dominated by weakened edge states. When λ is increased to λ2, T ≈ 0 and corner states emerge in the PLDOS distribution. At λ3 prior to the resonance, corner states are dressed by incident electrons and delocalize from corners. A perfect resonant peak T = 1 occurs at λ4, where the corresponding PLDOS occupies almost entire bulk of the system in Fig. 3(e). For a narrower transmission peak at λ5, the resonant state in Fig. 3(f) even extends to the λ = 0 region. Notice that the central region is a bulk insulator with localized corner states, and Figs. 3(e) and (f) show that global resonant states spread in this bulk insulating region through resonant tunneling process. Here PLDOS plays an important role in revealing the dynamic distribution of scattering states.
Similar extended states have been reported in Anderson localized systems [
56], which are known as the necklace states [
57–
60]. In a disordered Anderson insulator, local states trapped between potential barriers widely exist but hardly contribute to transport. As a rare event, the necklace state passes through many localized states through the multiresonance process and results in huge density of states (DOS) [
61,
62]. Here we identify another extended state: a global resonant state in bulk topological insulator, in which incident electrons combine multiple corner states and the scattering wavefunction occupies the entire scattering region. Although originating from different mechanisms, both the necklace state and this global resonant state lead to large transmission and DOS.
To summarize this subsection, we provide a 2D diagram of the transmission coefficient T as a function of both Fermi energy E and Zeeman field strength λ for the twoterminal system with four hollows. By varying E and λ, we see that there is a -shaped region with large T in the upper-right corner of Fig. 4. Meanwhile, the bright orange region indicates that the transmission approaches 1 under those circumstances. Besides, there are also many discrete points with high transmission, which correspond to global resonant events in bulk topological insulator. Numerical results displayed in Fig. 4 further confirm that, the probability of detecting topological corner states through transport measurement is greatly increased in the multihollow system.
3.2 Single-electron source via single-parameter pumping
Single-electron source is an on-demand nanoscale device in quantum computing and quantum information [
63], which reveals the dynamic aspect of coherent transport. Quantum parametric pump provides an effective way to the precise control of single-electron emitting or charge quantization [
64], which can be operated in both adiabatic and nonadiabatic regimes [
65–
67]. Recently, a silicon single-electron source with internal electron motion up to 250 GHz has been experimentally realized [
69]. Based on the well-separated energy levels of topological corner states in multihollow structures, we propose a singleelectron source through single-parameter pumping. The system setup is shown in the inset of Fig. 5(b), where a periodic gate voltage is applied on the two-probe multihollow system and accurately controls its internal potential,
V(
t) =
Vp cos(
ωt). Here
Vp is the amplitude of pumping potential and
ω is its frequency.
We focus on the adiabatic regime. In adiabatic pumping, the internal potential changes slowly compared with the relaxation time, and the system is under equilibrium at any moment of the whole pumping period. Due to the time-dependent pumping potential
V(
t), the instantaneous pumped current detected in the left lead is expressed as [
68] (
q = 1)
f is the Fermi distribution function. At zero temperature, the current can be further simplified as
where the retarded Green’s function Gr = (E−HC− is time-dependent. In the calculation, we set Vp = 0.01 and ω = 1. Note that here ω is irrelevant to the energy scale and simply stands for a long pumping period T = 2π/ω. The average pumped current in one period is given by , which is zero in this single-parameter pumping setup.
We first plot the transmission function of the two-probe multihollow system versus the Fermi energy at λ = 0.3 in the absence of the pumping potential. Figure 5(a) shows that there are four sharp resonant peaks (T = 1) in the transmission spectrum, which are well separated from one another. These discrete and perfect transmission peaks are ideal for single-electron manipulation. We further calculate the instantaneous pumped current ip(t) as a function of the pumping potential Vp cos(t). The Fermi energy is fixed at E = 0 in both leads. From Fig. 5(b), we see that ip(t) behaves as extremely sharp spikes when V (t) varies with time. It is obvious that the current spikes are closely related to the transmission peaks, where a narrower transmission peak gives rise to a higher current spike.
The discrete nature of ip(t) allows us to investigate the quantization of pumped charge. In this single-parameter adiabatic pump, the total pumped charge or current is zero in one period, since the current flowing in the central region in the first half-period will be pumped out in the second half-period. But we can integrate the current spike one by one to obtain the pumped charge, i.e., ∆Qp = ∫∆tip(t)dt. Numerical result gives that:
These quantized charges are strong evidences of singleelectron pumping. Therefore, this two-probe system with multiple topological corner states can serve as a singleelectron source via single-parameter pumping mechanism. Given the diagram in Fig. 4, we expect that this singleelectron source can operate in a large range of system parameters.
4 Conclusion
In summary, we have investigated the quantum transport properties of a diamond-shaped honeycomb lattice with multihollow structure, where multiple second-order topological states reside in the energy spectrum as well as real-space corners. Particularly, for a system with four hollows, five pairs of corner states exist. After connecting to two narrow leads, multiple resonant peaks are found in the transmission spectrum for a relatively large energy range. Partial local density of states shows that the resonant behavior corresponds to a global resonant state, in which the scattering state expands to the whole bulk insulator. These global resonant states are induced by incident electron dressing, and reveal the dynamic transport features of topological corner states. Based on the discrete energy levels of multiple corner states in the multihollow system, we propose a single-parameter pumping setup to realize single-electron source, and numerical results confirm that the pumped charge is perfectly quantized and well separated from one another.