NANOSYSTEMS: PHYSICS, CHEMISTRY, MATHEMATICS, 2017, 8 (3), P. 305-309
Quantum graphs with the Bethe-Sommerfeld property
P. Exner1'2, O. Turek2'3'4
1Doppler Institute for Mathematical Physics and Applied Mathematics, Czech Technical University, Brehova 7, 11519 Prague, Czechia
2 Department of Theoretical Physics, Nuclear Physics Institute CAS, 25068 ReZ near Prague, Czech Republic
3Bogoliubov Laboratory of Theoretical Physics, Joint Institute for Nuclear Research, 141980 Dubna, Russia
4Laboratory for Unified Quantum Devices, Kochi University of Technology, Kochi 782-8502, Japan
[email protected], [email protected] PACS 03.65.-w, 02.30.Tb, 02.10.Db, 73.63.Nm DOI 10.17586/2220-8054-2017-8-3-305-309
In contrast to the usual quantum systems which have at most a finite number of open spectral gaps if they are periodic in more than one direction, periodic quantum graphs may have gaps arbitrarily high in the spectrum. This property of graph Hamiltonians, being generic in a sense, inspires the question about the existence of graphs with a finite and nonzero number of spectral gaps. We show that the answer depends on the vertex couplings together with commensurability of the graph edges. A finite and nonzero number of gaps is excluded for graphs with scale invariant couplings; on the other hand, we demonstrate that graphs featuring a finite nonzero number of gaps do exist, illustrating the claim on the example of a rectangular lattice with a suitably tuned ¿-coupling at the vertices.
Keywords: periodic quantum graphs, gap number, ¿-coupling, rectangular lattice graph, scale-invariant coupling, Bethe-Sommerfeld conjecture, golden mean.
Quantum graphs [1] have attracted much attention both from the practical point of view as models of nanos-tructures as well as a tool to study the properties of quantum systems with a nontrivial topology of the configuration space. The topological richness of quantum graphs allows them to exhibit properties different from those of the 'usual' quantum Hamiltonians; examples are well known, for instance, the existence of compactly supported eigen-functions on infinite graphs [1, Sec. 3.4] or the possibility of having flat bands only as is the case for magnetic chain graphs with a half-of-the-quantum flux through each chain element [2].
In this letter, we are going to consider another situation where quantum graphs are known to behave unusually. Our problem concerns the gap structure of the spectrum of periodic quantum graphs. Recall that the finiteness of the open gap number for periodic quantum systems in dimension two or more was conjectured in the early days of quantum theory by Bethe and Sommerfeld [3]. The validity of the conjecture was taken for granted even if its proof turned out to pose a mathematically difficult problem. It took a long time before it was rigorously established for the 'usual' periodic Schrodinger operators [4-8]. Nevertheless, the situation appears to be different for quantum graphs, as we will see below.
The traditional reasoning behind the Bethe-Sommerfeld conjecture relies on the behavior of the spectral bands identified with the ranges of the dispersion curves or surfaces which, in contrast to the one-dimensional situation, typically overlap, making opening of gaps more and more difficult as we proceed to higher energies. The situation with graphs might be similar [1, Sec. 4.7] but the spectral behavior need not be the same, one reason being the possibility of resonant gaps. The existence of gaps coming from a graph decoration was first observed in the discrete graph context [9] and the effect is present for metric graphs as well [1, Sec. 5.1]. In addition, the recently discovered universality property of periodic graphs [10], valid in the generic situation when the graph edges are incommensurate and the vertex coupling is the simplest possible, usually called Kirchhoff, shows that the occurrence of infinitely many open gaps is quite typical.
This prompts one to ask whether there are quantum graphs with the band spectrum similar to that of the 'usual' multidimensional periodic systems, i.e. a nonzero and finite number of open gaps; for the sake of brevity, we shall speak of the Bethe-Sommerfeld property. With this question in mind, our aim in this letter is twofold. On the one hand we will show that the answer depends on the vertex coupling and there are classes of couplings for which such a behavior is excluded. On the other hand, using a simple example we shall demonstrate that periodic graphs with Bethe-Sommerfeld property do exist.
Consider an infinite graph r periodic in v directions, with a slight abuse of notation we will speak of a Zv-periodicity, v > 2. The Hamiltonian is supposed to act as —d2/dx2 on each edge; to make it a self-adjoint operator, one has to impose appropriate coupling conditions at each vertex. The most general form of them [11,12] is (U — I)^ + i(U + I= 0, where are vectors of the function and derivative values at the vertex,
respectively, and U is a unitary n x n matrix for n denoting the degree of the vertex. According to the eigenvalues of U one can split the coupling into the Dirichlet, Neumann, and Robin parts [1, Sec. 1.4]. Equivalently, the coupling can be written in so-called ST-form [13],
( I(r) T \ ^ ( S 0 \
^ 0 0 J = y —T* I(n_r) J ( )
for certain r, S, and T, where I(r) is the identity matrix of order r and S is a Hermitian matrix that refers to the Robin part. A coupling is called scale-invariant if the matrix U has no eigenvalues other than ±1; it is easy to see that this happens if and only if S = 0 [14]. Recall that the on-shell scattering matrix S(k) for the vertex in question is in the ST-formalism given by
S(k) = —I(n) + 2 ^ y? j (j(r) + TT* — 1 S^j (l(r) t)
and it is obvious that S(k) is independent of k iff S = 0.
The spectrum is obtained using the Bloch-Floquet theory [1, Sec. 4.2]. We cut from r its elementary cell rper which is assumed to be a finite graph with a family of pairs of 'antipodal' vertices related mutually by the conditions ^(v+) = e1^^(v_) and (v+) = e1^^'(v_) with some G (—n, n], where l = 1,..., v with v being the dimension of translation group associated with graph periodicity; the pair of edges with the endpoints v± can be turned into a single edge by identifying these endpoints, and the acquired phase $i coming from the Bloch conditions can be also regarded as being induced by a suitable magnetic potential.
The spectral problem can be solved in the usual way, cf. [1, Sec. 2.1] or [10]. Assuming that rper has E edges, we consider three 2E x 2E matrices. The diagonal matrix L is determined by the lengths of the directed edges (bonds) of rper, the diagonal matrix A has the entries e1^ or e-1^ at the positions corresponding to the edges created by the mentioned vertex identification, and all its other entries are zero, and finally, the matrix S is the bond scattering matrix, which contains directed edge-to-edge scattering coefficients. Using them, we define
F(k;$) := det (i — e1(A+fcL)S(k)) ; (2)
then k2 g <r(H) holds iff there is a $ G (—n, n]v such that the secular equation F(k; $) = 0 is satisfied.
Suppose now that the couplings at all the vertices of r are scale-invariant. This, in particular, means the matrix S entering formula (2) is independent of k, hence the value F(k; $) depends on the vectors $ and ki0, ki^..., kid, where {io,^,... , id}, d < E — 1, is the set of mutually different edge lengths of r, and F(k; $) is obviously 2n-periodic in the terms k^0, k^1,..., The secular equation can be then written as
F ({kM(2n), {kM(2n), ..., |kM(2n); $) = 0 , (3)
where {x}(2n) := 2n{x/(2n)}, for {•} denoting the difference between the number and its nearest integer. This allows us to prove that
(i) if <r(H) has a gap, then it has infinitely many gaps,
(ii) the gaps can be divided into series with asymptotically constant lengths with respect to k, and
(iii) in particular, if all the edge lengths are commensurate, the momentum spectrum is periodic.
The easiest part to verify is (iii). In that case there is an elementary length L > 0 and integers mj G N such that ij = mjL holds for j = 0,1,..., d, and consequently, the left-hand side of (3) is 2n/L-periodic with respect to k. Parts (i) and (ii) are more involved and we just sketch the argument referring to [15] for the full proof.
To prove (i) we consider a k > 0 satisfying k2 G ^(H) and use it to prove the existence, for any given C > 0, a k' > C such that (k')2 G ^(H). Due to the continuity of F, it is sufficient to find a k' so that the values k'ij are arbitrarily close to kij up to an integer multiple of 2n. To this aim, we denote a.j = ij /i0 and employ the simultaneous version of the Dirichlet's approximation theorem by which for any N g N there are integers p1,... ,pd, q G Z, 1 < q < N, such that
Pj
jq
< qN/d. (4)
Choosing integers m > i0C/2n and N > (2nm/6) , and putting k's := k + 2nmq/£0, it is straightforward to see that k's > C and, using (4), to check that | {k's j — kj}(2n) | <6 for all j. Moreover, the latter inequality in combination with the argument used to prove (i) yields the claim (ii).
In fact, one can exclude the Bethe-Sommerfeld property for a wider class of graphs. Given a vertex coupling described by condition (1), we consider the associated scale-invariant one obtained by replacing the Robin part S by zero. The vertex scattering matrix can be then written as S(k) = S0 + S1 (k)/k, where S0 is the scattering matrix of the associated scale-invariant vertex coupling. Then the function F(k; ■) in the secular equation is
Fo({Mo}(2п), {kM(2n^
,{k£d}(2n); 0) + 1 Fi(k; 0),
where the subscript zero at F0 refers to the Hamiltonian H0 of the graph in which all the couplings have been replaced by the associated scale-invariant ones. Using the fact that the leading behavior of F(•; ■) at high energies comes from the scale-invariant term, it is not difficult to see that [15]
(i) if <t(Ho) has an open gap, then a(H) has infinitely many gaps,
(ii) if all the edge lengths of r are commensurate, then the gaps of a(H) asymptotically coincide with those of *(Ho).
Let us turn to our second topic, the existence of graphs with the Bethe-Sommerfeld property. To this goal, we revisit the lattice Kronig-Penney model introduced in [16] and further discussed in [17,18]. In this case, r is a rectangular lattice graph in the plane with edges of lengths a and b, cf. Fig. 1.
Fig. 1. The rectangular-lattice graph
The Hamiltonian H = Ha,a/b is the negative Laplacian with the S coupling condition in each vertex, i.e. the functions are continuous there and satisfy J24=1 ^'(v) = a^(v) with a fixed parameter a e R.
According to [17], a number k2 > 0 belongs to a gap if and only if k > 0 satisfies the gap condition, which reads
( ka n tan---
V 2 2
ka n
kb n
+ tan---
+ 1 2 2
kb n
<
2k
(5)
for a > 0 and the analogous one with tan replaced by cot if a < 0; we neglect the case a = 0 where the spectrum is trivial, a(H) = [0, to). For a < 0 the spectrum extends to the negative part of the real axis and may have a gap there, but since such a gap always has a positive part [18], we may restrict our attention to examining gaps in the positive spectrum.
The crucial quantity is the ratio 6 = a/b. It is obvious that a(H) has infinitely many gaps once a = 0 and 6 is rational, and the same is true for the 'well approximable' irrationals [17]. We thus focus on the other irrationals, called badly approximable, i.e those to which there is a c > 0 such that
B-P
> — q2
for all p,q G Z with q = 0. These numbers form a set of zero Lebesgue measure. Alternatively they can be characterized as irrationals for which the sequence in the continued-fraction representation, 6 = [c0; c1, c2,...], is bounded [19], or as numbers whose Markov constant p(6), defined [20] for 6 G R as
M(B) = inf{c > 0 | (3TO(p,q) e N2)(
B-P
< — q2
,
(6)
is strictly positive. It is convenient to introduce a one-sided analog v(6) of the Markov constant, with the last inequality in (6) replaced by 0 <6 — p/q < c/q2. We have ^(6) = min{u(6), u(6-1)}; the number u(6) may or may not coincide with the Markov constant [21].
a
c
q
c
q
To obtain the existence claim, we focus on the situation where 9 is the 'worst approximable' irrational, the golden mean, ^ = (%/5 + 1)/2 . Our result about the golden-mean lattice is the following:
(i) If a > n2/(%/5a) or a < — n2/(%/5a), the spectrum has infinitely many gaps.
(ii) If —2n/a tan ((3 — V5)n/4) < a < n2/(V5a) there are no gaps in the spectrum.
(iii) If — n2/(V5a) < a < —2n/a tan ((3 — V5)n/4), there is a nonzero and finite number of gaps in the spectrum.
(iv) Moreover, put Aj := 2n (^2j — ^_2j)/%/5 tan (n^_2j/2), then there are exactly N gaps in the spectrum if
—An+i < a < — An .
Note that the window in which Bethe-Sommerfeld property occurs in this example (statement (iii)), is rather narrow, roughly 4.298 < —aa < 4.414.
The proof of these claims is rather involved and we limit ourselves with mentioning its key elements referring to [15] for the full exposition. The central notion is that of the Diophantine approximation of third type from below (from above, respectively). The former is a number p/q with p, q G Z such that
0 < q(q9 — p) < q'(q'9 — p') (7)
holds for all p'/q' > 9 with p'/q' = p/q, p', q' G Z and 0 < q' < q; the approximation from above has (7) replaced by 0 < q(p — q9) < q'(p' — q'9). We note that u(9) is the infimum of those q(q9 — p) for which p/q is a best approximation from below to 9. These approximations are also closely related to convergents obtained from truncated continued-fraction representation of 9. Specifically, every best approximation of the third kind from below to a 9 g R is a convergent of 9, and on the other hand, every best approximation from above is either [9] or a convergent of 9, where [•] is the ceiling function.
The described Diophantine approximation in combination with the gap condition allows us to estimate the number of gaps for a given ratio 9 = a/b and coupling parameter a. This has to be done for each sign of a separately. If a > 0, condition (5) yields that
• if a < n2 • min {u(9)/b, u(9-1)/^, the spectrum has at most finitely many gaps. If the opposite (sharp) inequality holds true, the number of gaps is infinite;
• if a < y+ for y+ given by
Y+ := min inf < 2nmw — tan ( — (mn — Imnl )] >, (8)
n=e,e-i men\ y ab V2v ' L /JVJ' w
the spectrum has no gaps. If a > y+, there are gaps in the spectrum;
• in particular, if y+ < a < n2 • min {u(0)/6, , there is a nonzero and finite number of gaps in
the spectrum.
The golden mean ^ in our example has continued-fraction representation ^ = [1; 1,1,... ], therefore, its convergents are ratios Fn+1/Fn of the Fibonnaci numbers Fn = - (-^)-n)/\/5, and we have = u(^) = 1/%/5 in view of the Hurwitz theorem [22]. Hence, we get
and n2 -min'
Y+
■\/5a [ b ' a J A/5a'
which implies for all positive a either infinite number of spectral gaps or none at all. Alternatively, for a < 0 the gap condition implies that
• if |a| < n2 • min {u(9)/a, u(9-1)/^, the number of gaps in the positive spectrum is at most finite. By contrast, for | a| greater than the right-hand side of the above inequality, there are infinitely many spectral gaps, while
• if |a| < y_ for y_ given by the relation analogous to (8) in which mn — |_mnj is replaced by [mn] — mn, the spectrum has no gaps. If |a| > y_, there are gaps in the spectrum;
• in particular, if y_ < |a| < n2 • mi^u(9)/a, u(9_1)/^, there is a nonzero and finite number of gaps in the spectrum.
The last of these results, together with the easily verifiable formula
2n (3 — V5)n Y_ = — tan ■
a 4
proves the claim (iii) above, demonstrating thus the existence of graphs with the Bethe-Sommerfeld property. To prove (iv), the first step consists in showing that the number of gaps in the golden-mean lattice graph is equal to the number of solutions m G N of
2пт ( п \
-tan — ( I тф I — тф) < |a| .
a V 2 /
If — An +1 < a < —An, one can verify that the inequality is satisfied only for m = Fn with n = 2,4,6,..., 2N; this implies the existence of exactly N gaps.
Since a finite nonzero number of gaps occurred in the above example only for attractive vertex couplings, it is natural to ask whether the attractivity of the coupling is always a necessary condition for the Bethe-Sommerfeld property. It appears that it is not the case, a more thorough analysis of the gap condition [15] shows that, for instance, the edge ratio
л 2t3 — 2t2 — 1 + V5 , _
в = —-3-^———- with t G N, t > 3 ,
2(t4 — t3 + t2 — t +1) , > ,
which has the continued-fraction representation [0; t,t, 1,1,1,1,...], yields the lattice graph spectrum with the said
property for a certain a > 0 and for a certain a < 0 as well. This observation can be stated more generally [15]
and allows one to explicitly construct ratios to achieve the Bethe-Sommerfeld property.
In conclusion, we have proved and demonstrated on concrete examples that there are periodic quantum graphs,
the spectrum of which contains a nonzero and finite number of open gaps. We have also described how this
property depends on the type of the vertex coupling; in particular, we showed that a quantum graph cannot be of
the Bethe-Sommerfeld type if its couplings are scale invariant or associated with scale-invariant ones.
Acknowledgements
The research was supported by the Czech Science Foundation (GACR) within the project 17-01706S. References
[8 [9 [10 [11 [12 [13
[14 [15 [16 [17 [18 [19 [20 [21 [22
Berkolaiko G., Kuchment P. Introduction to Quantum Graphs. AMS, Providence, R.I., 2013.
Exner P., Manko S. S. Spectra of magnetic chain graphs: coupling constant perturbations. J. Phys. A: Math. Theor., 2015, 48, P. 125302. Sommerfeld A., Bethe H. Electronentheorie der Metalle. 2nd ed., Handbuch der Physik, Springer Verlag, Berlin, 1933. Skriganov M. M. Proof of the Bethe-Sommerfeld conjecture in dimension two. Soviet Math. Dokl., 1979, 20, P. 956-959. Dahlberg J., Trubowitz E. A remark on two dimensional periodic potentials. Comment. Math. Helvetici, 1982, 57, P. 130-134. Skriganov M.M. The spectrum band structure of the threedimensional Schrödinger operator with periodic potential. Invent. Math., 1985, 80, P. 107-121.
Helffer B., Mohamed A. Asymptotic of the density of states for the Schrödinger operator with periodic electric potential. Duke Math. J., 1998, 92, P. 1-60.
Parnovski L. Bethe-Sommerfeld conjecture. Ann. Henri Poincare, 2008, 9, P. 457-508.
Schenker J. H., Aizenman M. The creation of spectral gaps by graph decoration. Lett. Math. Phys., 2000, 53, P. 253-262.
Band R., Berkolaiko G. Universality of the momentum band density of periodic networks. Phys. Rev. Lett., 2013, 111, P. 130404.
Kostrykin V., Schrader R. Kirchhoff's rule for quantum wires. J. Phys. A: Math. Gen., 1999, 32, P. 595-620.
Harmer M. Hermitian symplectic geometry and extension theory. J. Phys. A: Math. Gen., 2000, 33, P. 9193-9203.
Cheon T., Exner P., Turek O. Approximation of a general singular vertex coupling in quantum graphs. Ann. Phys. (NY), 2010, 325,
P. 548-578.
Cheon T., Exner P., Turek O. Tripartite connection condition for a quantum graph vertex. Phys. Lett. a, 2010, 375, P. 113-118.
Exner P., Turek O. Periodic quantum graphs from the Bethe-Sommerfeld perspective, arXiv: 1705.07306 [math-ph]
Exner P. Lattice Kronig-Penney models. Phys. Rev. Lett. 1995, 74, P. 3503-3506.
Exner P. Contact interactions on graph superlattices. J. Phys. A: Math. Gen., 1996, 29, P. 87-102.
Exner P., Gawlista R. Band spectra of rectangular graph superlattices. Phys. Rev. B, 1996, 53, P. 7275-7286.
Khinchin Y. Ya. Continued Fractions. University of Chicago Press, 1964.
Cassels J. W. An introduction to Diophantine Approximation. Cambridge University Press, 1957.
Pelantova E., Starosta S., Znojil M. Markov constant and quantum instabilities. J. Phys. A: Math. Theor., 2016, 49, P. 155201. Hurwitz A. Uber die angenaherte Darstellung der Irrationalzahlen durch rationale Brüche. Math. Ann., 1891, 39, P. 279-284.
[7