Transition between extended and localized states in a one ...
Transition between extended and localized states in a one ...
Transition between extended and localized states in a one ...
Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
TRANSITION BETWEEN EXTENDED AND LOCALIZED . . . PHYSICAL REVIEW A 64 033416<br />
FIG. 3. Energy b<strong>and</strong>s for V 00.5 <strong>and</strong> 0 as a function of<br />
rational p/q; values with pq30 are shown. Each po<strong>in</strong>t <strong>in</strong><br />
this plot corresponds to a full narrow energy b<strong>and</strong>.<br />
structure with large energy gaps, rem<strong>in</strong>iscent of the Hofstadter<br />
butterfly 16. This occurs at energies of the order of<br />
V 0.<br />
When is an irrational number Bloch’s theorem is no<br />
longer applicable <strong>and</strong> the calculation of the spectrum is nontrivial.<br />
This <strong>in</strong>commensurate case can be studied by tak<strong>in</strong>g a<br />
sequence of rational numbers np n /q n such that it converges<br />
to as n→. The potential with n is called a periodic<br />
approximant to the <strong>in</strong>commensurate potential. One important<br />
irrational value for is the <strong>in</strong>verse of the golden<br />
mean, (51)/2, which characterizes the quasiperiodicity<br />
<strong>in</strong> some quasicrystals 17 <strong>and</strong> possesses <strong>in</strong>terest<strong>in</strong>g<br />
number-theoretical properties. The sequence of approximants<br />
used is obta<strong>in</strong>ed by truncat<strong>in</strong>g the expansion of <strong>in</strong> terms of<br />
a cont<strong>in</strong>ued fraction, result<strong>in</strong>g <strong>in</strong> the recursion p n<br />
q n1 ,q np n1q n1 with <strong>in</strong>itial values p 0q 01. In<br />
the rest of the paper we will study this particular case, although<br />
the results are similar for other irrational values of .<br />
A. Label<strong>in</strong>g the gaps<br />
As we <strong>in</strong>crease the order n of the approximation, qn <strong>in</strong>creases<br />
<strong>and</strong> thus the size of the Brillou<strong>in</strong> z<strong>one</strong> shr<strong>in</strong>ks, converg<strong>in</strong>g<br />
to zero. With gaps open<strong>in</strong>g up at both the center <strong>and</strong><br />
the edges of the Brillou<strong>in</strong> z<strong>one</strong>, as n→ gaps open up almost<br />
everywhere <strong>in</strong> the spectrum; the spectrum becomes a<br />
Cantor set, which is a fractal 18. It is of <strong>in</strong>terest to f<strong>in</strong>d<br />
properties of the system that rema<strong>in</strong> cont<strong>in</strong>uous <strong>in</strong> this limit.<br />
As shown by Fig. 3, <strong>one</strong> such property is the location of the<br />
ma<strong>in</strong> gaps <strong>in</strong> the spectrum.<br />
The total number of gaps <strong>in</strong> the spectrum is <strong>in</strong>f<strong>in</strong>ite, but<br />
countable. They can be <strong>in</strong>dexed us<strong>in</strong>g the gap-labell<strong>in</strong>g theorem<br />
19. Given the density of <strong>states</strong> per unit length (E),<br />
let us def<strong>in</strong>e the <strong>in</strong>tegrated density of <strong>states</strong> I(E)<br />
<br />
E (E)dE. When evaluated at an energy ly<strong>in</strong>g on an<br />
energy gap E gap there exists a unique pair of <strong>in</strong>tegers n 1 <strong>and</strong><br />
n 2 such that<br />
IE gapn 11/2n 21/2. 3.2<br />
This pair (n 1 ,n 2) identifies the gap.<br />
It is <strong>in</strong>structive to plot the <strong>in</strong>tegrated b<strong>and</strong>width as a function<br />
of the energy. The plots for different rational approxi-<br />
033416-3<br />
FIG. 4. Integrated b<strong>and</strong>width for V 00.068, show<strong>in</strong>g the labels<br />
for the ma<strong>in</strong> gaps <strong>in</strong> the spectrum.<br />
mants quickly converge to a curve <strong>in</strong> the region support<strong>in</strong>g<br />
<strong>extended</strong> <strong>states</strong>; the ma<strong>in</strong> gaps are displayed as horizontal<br />
segments. In Fig. 4 we show the result for V 00.068, for<br />
which the whole spectrum supports <strong>extended</strong> <strong>states</strong>. The<br />
largest gaps correspond to comb<strong>in</strong>ations of small <strong>in</strong>tegers n 1<br />
<strong>and</strong> n 2. If we th<strong>in</strong>k of the potential as a perturbation of the<br />
free-particle Hamiltonian, then n 1n 2 is the order of the<br />
correction that gives rise to the gap. Thus, <strong>in</strong>tegers with<br />
larger magnitudes give smaller gaps.<br />
B. Quasimomentum space distributions<br />
For a periodic potential we can calculate the energy as a<br />
function of the quasimomentum for the different energy<br />
b<strong>and</strong>s. In the quasiperiodic case the quasimomentum is not a<br />
well def<strong>in</strong>ed quantity, so there is no dispersion relation. In<br />
order to study a closely related concept, we plot the eigenfunctions<br />
of the Hamiltonian as both a function of momentum<br />
p <strong>and</strong> energy E.<br />
Let us consider first the periodic potential,<br />
VxV 0cosx1cosx/2. 3.3<br />
In the top panel of Fig. 5 we have plotted the energy E as a<br />
function of the momentum p for V 00.1. The gray scale <strong>and</strong><br />
size of the po<strong>in</strong>t represents the probability of measur<strong>in</strong>g a<br />
value of p for a given eigenvector of the Hamiltonian with<br />
energy E. In the absence of a potential, the plot would just be<br />
the free-particle dispersion relation,<br />
E 1<br />
2 p2 , 3.4<br />
shown as the thick parabola <strong>in</strong> the lower panel. In the presence<br />
of the weak periodic potential 3.3, Bragg scatter<strong>in</strong>g<br />
with p 1<br />
2 ,1 becomes possible. Consider<strong>in</strong>g multiplescatter<strong>in</strong>g<br />
events the momentum can change by <strong>in</strong>teger multiples<br />
of 1<br />
2 . We <strong>in</strong>clude <strong>in</strong> the lower panel the free-particle<br />
parabola displaced <strong>in</strong> momentum by such amounts.<br />
We can clearly see <strong>in</strong> the top panel of Fig 5 the remnants<br />
of the free particle parabola, as well as some of the secondary<br />
<strong>one</strong>s. The gaps <strong>in</strong> the spectrum open up where the energy<br />
eigen<strong>states</strong> are degenerate, which occurs where two of the<br />
parabolas meet. Interpret<strong>in</strong>g the secondary parabolas as<br />
Bragg scattered <strong>states</strong>, the gray scale of the po<strong>in</strong>t signals the