25.07.2013 Views

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 ...

SHOW MORE
SHOW LESS

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

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!