Periodic chain of disks in a magnetic field: bulk states and edge states
E.N. Grishanov1, D.A. Eremin1, D. A. Ivanov1, I.Yu. Popov2, P.I. Smirnov2
1 Department of Mathematics and IT, Ogarev Mordovia State University Bolshevistskaya Str. 68, Saransk, Russia 2Department of Higher Mathematics, ITMO University, Kroverkskiy pr. 49, St. Petersburg, 197101, Russia popov1955@gmail.com
PACS 73.23.Ad, 02.30.Tb DOI 10.17586/2220-8054-2015-6-5-637-643
An explicitly solvable model for periodic chain of coupled disks in orthogonal magnetic field is considered. The spectrum for the Hamiltonian is compared with the spectrum for the corresponding chain of circles. These models are used for the comparison of the bulk and edge states. It is found that for some range of the magnetic field values the lowest band for the circles system lies below the spectrum for the corresponding disks system, i.e. the edge band is below and is separated from the lowest bulk band. Keywords: nanostructure; magnetic field; explicitly solvable model. Received: 1 September 2015
1. Introduction
The difference between surface and bulk conductivities for materials and nanostructures is a subject of numerous experimental and theoretical investigations due to intriguing prospects (high-temperature superconductivity, topological insulator, etc.). In many cases, low-dimensional nanostructures have unusual electronic properties, see, e.g., results of recent experiments with highly ordered organometallic nanoribbons, whose intrinsic (defect-free) conductivity is found to be three orders of magnitude higher than that of macroscopic crystals [1]. It should be mentioned that in many cases there is no adequate mathematical model of the corresponding nanosystem, and the theoretical description is complex. Our goal is to construct a mathematical model for chain type periodic nanostructure and to describe its spectrum. Ideally, we would like to have a band (preferably, the lowest band) for the surface (edge) states be in the gap (or below the first band) for the bulk states. The corresponding system is described in the present paper. This case is, in some sense, similar to a topological insulator (TI). TIs are electronic materials that have a bulk band gap like an ordinary insulator, but have protected conducting states on their edge or surface [2-5]. We have similar phenomenon for surface state protection as in TI, although the physical nature is different.
We consider a periodic system of coupled disks as a basic geometry. As for the problem of edge states, our consideration is in the framework of a model of quantum graphs (periodic system of coupled circles). This model was used intensively in previous publications (starting from works [6, 7]) due to its advantages: in one way, it is quasi-one-dimensional and rather simple, however, despite this apparent simplicity, it preserves the properties of the corresponding physical system, particularly, its spectral properties [8-13]. As for the problem of bulk states, we deal with the model of quantum resonators coupled through point-like window. This model
is based on the theory of self-adjoint extensions of symmetric operators (see, e.g., [21]). We compare the spectra of the Landau operators for these two models and observe that for some values of the magnetic field, one has an edge state band below the first bulk band. The structure of the paper is rather natural. We describe, consequently, the models for chain of circles and for chain of disks, and then compare the corresponding spectra.
2. Chain of circles
Consider an infinite chain of circles Cn, n = 1, 2,... of radii a. The state space is
H = 0 L2(Cn). To
construct the Hamiltonian of the charge particle for the chain of coupled
n=1
circles, we start from self-adjoint operator Hcc, Hcc = HC. Here H is the Landau operator
n=1
for the circle Cn. Using the standard polar coordinates (r, p) on the circle, one can represent the Hamiltonian in the following form:
Hc 2ma2 ( 1 dp + $0
h
Here, m is an electron mass, $ = na B is a magnetic flux, $0 = — is the quantum of
2c h2
the magnetic flux. We use the unit system such that the factor -- is equal to 1 (i.e.
2ma2
h =1 ,m =1/2, a =1).
A rigorous mathematical model of our system of coupled circles is constructed by a conventional way by use of the operator extensions theory (see, e.g., [14-19]). We assume that the contacts are located at the opposite points (a, p1) and (a, p2), where p1 = 0, p2 = n (see
Fig. 1). We restrict the initial operator Hcc on the set D0 = {f e D(Hcc) : f(qn) = 0},
n=1
where qn is the contact point. The model Hamiltonian is given by a self-adjoint extension of this symmetric operator. We choose the extension corresponding to the following boundary conditions:
r+1(pn+1) = etqVn(pn),
^/n+1(pn+1 + 0) - ^'n+1(pn+1 - 0) = eiql[^/n(pn + 0) - ^/n(pn - 0)], ^n+1(pn+1 ) = ^n(pn)= ^
^/n+1(pn+1 + 0) - ^/n+1(pn+1 - 0) + ^/n(pn + 0) - ij/n (pn - 0)],
where l = 2a = 2, q is the quasi-momentum, ft is the coupling parameter for the contact points. Taking into account the expression for the electron wave function ^(p):
^(p) = a,1G(p,p1,E) + a2G(p,p2 ,E),
and (1), we obtain the following dispersion equation:
4Q21 cos(ql) + ft-1 det Q - 2Tr Q = 0, (2)
Fig. 1. Chain of circles (disks)
where Q is the Krein Q-matrix:
Qii(E) = Q22 (E ) = kl
sin nkr cos nkr
Q12 (E ) = Q2i(E) = kr
sin2 nn — sin2 nkr_ sin nkr cos nn
sin2 nn — sin2 nkr _ Solutions of (1) for various quasi-momenta give us the spectral bands.
3. Chain of disks
In this section, we consider a system with similar geometrical structure but circles in the chain are replaced by disks of the same radius Dn = {(r, : r < a}. Here, r, tp are the
standard polar coordinates. The initial Hamiltonian for our procedure is Hcd = Hn, where
n=1
Hn is the Landau operator for the disk Dn with the Dirichlet boundary conditions:
Hd = (inV-eBx. r
2m
We introduce the following notations:
| eB |
u =-, ^ =
cm
2c
hc \
jeBf )
1/2
Xo =
2^2
An eigenfunction of the operator Hcd with the Dirichlet boundary condition is as follows:
^(r, p) =
, 1/2
mu \ r
\ 2nhc,
in
(2^2)HI/2 exp(^ — 4^) K^111 + 1, ¿2
We denote by Eln, - = 0, ±1, ±2,... the corresponding eigenvalue:
Ein = hu ( ---
where e1n is the n—th root of the equation $(e, |1| + 1,x0) = 0, $ is the Kummer function. Here, c1n is the normalization constant:
xo
cin = exp(—x)x 11 1 $2(ein, |-| + 1,x)dX.
a
The Green function for a single disk has the following form [20]:
ru ' ' TP\ mu
G{r, p,r ,p; E) = —h exP
r2 + r'2
(rr')111
+L
exp(il(p — p' ))x
+L
£
n=1
9 /9
$(£.„, Ill + 1, ^|l| + 1, )
cln(Eln — E )
(3)
To construct the model for coupled disks, we use the operator extensions theory model, more precisely, the model of zero-width slit [21]. This model faces difficulties for the Dirichlet boundary condition. In this case, one has to extend the initial Hilbert space to the Pontryagin space with indefinite metrics. The mathematical structure of this model has been described previously, e.g., in [22]. Here, we present the resulting formulas for our case.
The wave function for the j-th disk has the form:
j 0G(r,p; r',p'; E)
rfi(r, P) = a1
dr'
+ dG(r,p; r',p'; E)
dr'
r'=1,(p' =n
The asymptotics of this function near the chosen (contact) points (r,p) = (1,0), (r,p) = (1,n) contain a singular term. To construct a self-adjoint extension in the Pontrya-gin space (i.e. the model Hamiltonian), one should introduce a linear relation between the main singular and regular terms of the asymptotics from the both sides of the contact points (i.e. in neighbor disks).
We assume the existence of the ^-potentials at the contact points. This gives us a rule for selection of the self-adjoint extension from the family of extensions.
As for the chain, we should take into account the periodicity. This can be made using the Bloch's theory (similarly to the case of circles) or by the transfer-matrix technique (see, e.g. [14]). The transfer-matrix M is determined as follows:
aj+i
b
j+i
M
bj
(4)
where:
j+i
aj+1 = a
bj = q12«2 + 911«1, bj+1
In this formula, the following notation is used:
(5)
j+1 j+1 q12<aJ2 + a\ qn.
911
d_ dr
dG(r, p; r',p'; E)
dr'
dG(r,p; r',p'; Eo)
r'=1,<p'=0
q12
j dG(r, p; r',p'; E)
dr
dr'
dr'
r'=1,<p'=0.
r'=1,<p'=0.
r=1,^=0
r=1,^=n
The above mentioned selection rule for the extension leads to the following relations:
a2+i + = a(a2qi2 + alqii), j i j j+i i j+i «2512 + aiqii = aj q2i + «2 Q22.
Here, a is the parameter describing the strength of the S—potential at the contact point.
a
j
a
j
Taking into account (5) and (6), we express aj+i, bj+i in terms of a, bj:
Oj+1
qi2 1 — aqn —aj--
qii qii
b
j+i
qii +
aqii — 1
q2i
522
+ bj — +
1 , (1 — aqii)(q2i— qy
1 — aqii qii J ' \qi2 qiiqi2
Using (4), one obtains from (7), (8) the formula for the transfer-matrix M:
( qi2 1 — qgii \
(7)
(8)
M=
V
qii qii
det Q 1 (1 — aqn)det Q
qii — 1----1--
qii qi2 qiiqi2
(9)
/
where det Q = q2i — ?22. One can check that det M =1. In order for E to belong to the continuous spectrum of the Hamiltonian, the modula of the eigenvalues of the matrix M = M(E) must be equal to 1 (see, e.g. [14]). In our case, it reduces to the following inequality:
|TrM(E)| < 2.
4. Results and discussion
The spectral bands of the Hamiltonian for the periodic chain of circles are obtained by solving equation (2) for various quasi-momenta q. The dependence of the dispersion equation
solution upon the quasi-momentum is shown in Fig. 2. The picture depends on the value of the
$
magnetic flux —. For the semi-integer flux (Fig. 2a), one has only levels (so-called, flat bands,
$o
$
i.e. infinitely degenerate eigenvalues). Fig. 2b shows the band structure for —- = 0.2.
$o
Fig. 2. The band structure for the circles chain for different values of the mag-
$$
netic flux: a - for — = 0.5; b - for — = 0.2.
$o $o
The dependence of the bands on the magnetic field for the disks system for different a parameter values is shown in Fig. 3.
The dependence of the magnetic field bands for both systems is shown in Fig. 4. Bands correspond to the dark areas. One can see that for chosen range of the magnetic field, the lowest band for the chain of circles lies below the lowest band for the chain of disks. Hence,
Fig. 3. The dependence of the bands on the magnetic field for the disks system for different a parameter values: a - for a = 0.5; b - for a =1; for a = 2.
for these magnetic field values, one has something similar to the topological insulator effect; the edge state lies below all bulk states.
B B
Fig. 4. The dependence of the energy bands on the magnetic field B: a - for the chain of disks; b - for the chain of circles.
Acknowledgements
This work was partially financially supported by the Government of the Russian Federation (grant 074-U01), by the Ministry of Science and Education of the Russian Federation (GOSZADANIE 2014/190, Project 14.Z50.31.0031), by grants of the President of Russia (state contracts 14.124.13.2045-MK and 14.124.13.1493-MK).
References
[1] C. Hermosa, J.V. Alvarez, et al. Intrinsic electrical conductivity of nanostructured metal-organic polymer chains. Nature Comm., 2013, 4, P. 1709.
[2] A. Medhi, V.B. Shenoy. Continuum theory of edge states of topological insulators: variational principle and boundary conditions. J. Phys.: Condens. Matter, 2012, 24 (35), P. 355001.
[3] M.Z. Hasan, C.L. Kane. Topological insulators. Rev. Mod. Phys., 2010, 82 P. 3045-3067.
[4] S. Mao, Y. Kuramoto, K.-I. Imura, A. Yamakage. Analytic theory of edge modes in topological insulators. J. Phys. Soc. Japan, 2010, 79, P. 124709.
[5] J.J. Cha, K.J. Koski, Yi. Cui. Topological insulator nanostructures. Phys. Status Solidi RRL, 2013, 7 (1-2), P. 15-25.
[6] N.I. Gerasimenko, B.S. Pavlov. Scattering problems on noncompact graphs. Theoret. Math. Phys., 1988, 74, P. 230-240.
[7] P. Exner, P. Seba. Free quantum motion on a branching graph. Rep. Math. Phys., 1989, 28, P. 7-26.
[8] P. Exner, et al. (Eds.) Analysis on Graphs and Its Applications. Proc. Symp. Pure Math., (AMS, Providence, Rhode Island, 2008).
[9] O. Post. Spectral Analysis on Graph-like Spaces. Springer, Berlin, 2010.
[10] I.S. Lobanov, I.Yu. Popov. Scattering by zigzag-armchair nanotube junction. Nanosystems: Physics, Chemistry, Mathematics, 2012, 3 (2), P. 6-28.
[11] I.Yu. Popov, A.N. Skorynina, I.V. Blinova. On the existence of point spectrum for branching strips quantum graph. J. Math. Phys., 2014, 55, P. 033504/1-20.
[12] V. Kostrykin, R. Schrader. Kirchhoff's rule for quantum wires. J. Physics A: Math. Gen., 1999, 32, P. 595-630.
[13] K. Pankrashkin. On the Robin eigenvalues of the Laplacian in the exterior of a convex polygon. Nanosystems: Phys., Chem., Math., 2015, 6 (1), P. 46-56.
[14] P. Duclos, P. Exner, O. Turek. On the spectrum of a bent chain graph. J. Phys. A: Math. Theor., 2008, 41, P. 415206.
[15] I.Yu. Popov, P.I. Smirnov. Spectral problem for branching chain quantum graph. Phys. Lett. A, 2013, 377, P. 439-442.
[16] I.S. Lobanov, E.S. Trifanova. Direct and inverse problems in the model of quantum graphs. Nanosystems: Phys., Chem., Math., 2012, 3 (5), P. 6-32.
[17] I.S. Lobanov, A.I. Trifanov, E.S. Trifanova. Genetic algorithm for constructing graphene nanoribbon with given electronic transport properties. Nanosystems: Phys., Chem., Math., 2013, 4, P. 512-523.
[18] D.A. Eremin, I.Yu. Popov. Quantum ring with wire: a model of two-particles problem. Nanosystems: Phys., Chem., Math., 2012, 2 (2), P. 6-32.
[19] M.A. Kokoreva, V.A. Margulis, M.A. Pyataev. Electron transport in a two-terminal AharonovBohm ring with impurities. Physica E, 2011, 43, P. 1610.
[20] B.S. Pavlov, I.Yu. Popov, V.A. Geyler, O.S. Pershenko. Possible construction of a quantum multiplexer. Europhys. Lett., 2000, 52, P. 196.
[21] I.Yu. Popov. The resonator with narrow slit and the model based on the operator extensions theory. J.Math. Phys., 1992, 33 (11), P. 3794-3801.
[22] I.Yu. Popov. Model of point-like window for electromagnetic Helmholtz resonator. Zeitschrift fur Analysis und ihre Anwendungen. 2013, 32 (2), P. 155-162.