L-7: Difference between revisions

From Disordered Systems Wiki
Jump to navigation Jump to search
No edit summary
Line 1: Line 1:
<Strong> Goal </Strong>: This is the first lecture about the localization. Localization is a <Strong> wave phenomenon induced by disorder.</Strong>
'''Goal.''' This lecture introduces the phenomenon of localization. Localization is a '''wave phenomenon induced by disorder''' that suppresses transport in a system.


== Free particles and ballistic behaviour ==
== Free particles and ballistic behaviour ==
The Schrodinger equation govern the evolution of the quantum state of a particle in 1D:
<center><math> i\hbar \, \partial_t \psi(x,t) = H \psi(x,t), \quad H= -\frac{\hbar^2}{2 m} \frac{d^2}{dx^2} +V(x)  </math></center>
Here <math> H  </math> is the Hamiltonian. For a free particle the potential is absent, <math> V(x) =0</math>.


One first look for separable solutions, or eigenstates or stationary solutions. These solutions are stationary because if the particle is in an eigenstate, its physical properties are time independent. For a free particle, the separable solutions have the form
The Schrödinger equation governs the evolution of the quantum state of a particle in one dimension:
<center><math> \psi_k(x,t) =\frac{1}{\sqrt{2 \pi}} e^{i k x} e^{- i E_k t/\hbar},\quad \text{with}  \; E_k =\frac{\hbar^2 k^2}{2 m}  </math></center>
Here <math>k </math> is a real number, hence the spectrum is continuum. These solutions are delocalized on the entire real axis and they cannot be normalized.


The physical states of the particle can be decomposed on the eigenstates of the Hamiltonian. Indeed, For the superposition principle, any linear combination of separable solution is also a solution of the Schrodinger equation. However this superposed solution are not separable and the particle's properties will evolve in time. Hencewe can costruct a localized wave packet, with  the correct normalization:  
<math display="block">
<center><math> \psi(x,t)=\int_{-\infty}^\infty d k c(k) \psi_k(x,t), \quad  \text{with}  \;  \int_{-\infty}^\infty d k  |c(k)|^2 =1</math></center>
i\hbar \partial_t \psi(x,t)=H\psi(x,t),
\qquad
H=-\frac{\hbar^2}{2m}\frac{d^2}{dx^2}+V(x).
</math>
 
Here <math>H</math> is the Hamiltonian. For a free particle the potential vanishes, <math>V(x)=0</math>.
 
One first looks for separable solutions, also called eigenstates or stationary solutions. If the particle is in an eigenstate, all physical observables are time independent.
 
For a free particle the stationary solutions are plane waves
 
<math display="block">
\psi_k(x,t)=\frac{1}{\sqrt{2\pi}}e^{ikx}e^{-iE_k t/\hbar},
\qquad
E_k=\frac{\hbar^2k^2}{2m}.
</math>
 
Here <math>k</math> is a real number and the spectrum is continuous. These solutions are completely delocalized over the real axis and cannot be normalized.
 
Physical states are obtained as superpositions of eigenstates. By the superposition principle any linear combination of stationary solutions is again a solution of the Schrödinger equation. Hence we can construct a localized wave packet:
 
<math display="block">
\psi(x,t)=\int_{-\infty}^{\infty}dk\,c(k)\psi_k(x,t),
\qquad
\int_{-\infty}^{\infty}dk\,|c(k)|^2=1.
</math>
 
=== Evolution of a Gaussian wave packet ===
 
* '''Initial state.'''  
At time <math>t=0</math> consider the Gaussian packet
 
<math display="block">
\psi(x,0)=\frac{e^{-x^2/(4a^2)}}{(2\pi a^2)^{1/4}},
\qquad
|\psi(x,0)|^2=\frac{e^{-x^2/(2a^2)}}{\sqrt{2\pi a^2}}.
</math>
 
Show that the coefficients of the plane wave decomposition are
 
<math display="block">
c(k)=\left(\frac{2a^2}{\pi}\right)^{1/4}e^{-a^2k^2}.
</math>
 
* '''Time evolution.'''
 
Define the spreading velocity
 
<math display="block">
v_s=\frac{\hbar}{2ma^2}.
</math>
 
Show that the time evolution of the packet is
 
<math display="block">
\psi(x,t)=
\frac{e^{-x^2/(4a^2(1+i v_s t))}}
{[2\pi a^2(1+i v_s t)]^{1/4}}.
</math>
 
* '''Ballistic spreading.'''
 
The probability density becomes
 
<math display="block">
|\psi(x,t)|^2=
\frac{e^{-x^2/(2a^2(1+v_s^2 t^2))}}
{\sqrt{2\pi a^2(1+v_s^2 t^2)}}.
</math>
 
Hence
 
<math display="block">
\sqrt{\langle x^2\rangle}
=
\left(\int dx\,x^2|\psi(x,t)|^2\right)^{1/2}
=
a\sqrt{1+v_s^2 t^2}.
</math>
 
At long times
 
<math display="block">
\sqrt{\langle x^2\rangle}\sim v_s t.
</math>
 
This behaviour is called '''ballistic spreading'''.


===Evolution of the Gaussian packet===
It should be contrasted with the two other possible transport regimes:


* '''Diffusive motion''' (random walk): 
<math>\sqrt{\langle x^2\rangle}\sim \sqrt{t}</math>


* <Strong> Initial state </Strong>. At time <math> t=0  </math>, the wave packet and its probability density function are
* '''Localized regime''':  
<center><math> \psi(x,t=0) = \frac{e^{-x^2/(4a^2)}}{(2 \pi a^2)^{1/4}}, \quad  |\psi(x,t=0)|^2=\frac{e^{-x^2/(2a^2))}}{\sqrt{2 \pi a^2}}  </math></center>
<math>\sqrt{\langle x^2\rangle}</math> saturates at long times.
Show that the coefficent <math> c(k)  </math> are
<center><math> c_k=  \left(\frac{2 a^2}{\pi}\right)^{1/4} e^{-a^2 k^2}</math></center>


Understanding how disorder changes ballistic motion into diffusion and eventually localization is the main goal of the following lectures.


== Localization of the packet: general idea and experiment ==


* <Strong>Time evolution </Strong>. We define the spreading speed <math> v_s=\hbar/(2 m a^2) </math>. Show that the time  evolution of the Gaussian wave packet has the form:
[[File:Localization1DB.png|thumb|left|x240px|BEC condensate expanding in a 1D disordered potential. Billy et al., Nature (2008).]]
<center><math> \psi(x,t) = \frac{e^{-x^2/(4 a^2 (1+i v_s t))}}{[2 \pi a^2(1+ iv_s t)]^{1/4}}  </math></center>


* <Strong> Ballistic Spreading</Strong>. Show that the Probability density function of the position is
In 1958, P. Anderson proposed that disorder can suppress transport in quantum systems. This phenomenon, called '''Anderson localization''', has since been observed both numerically and experimentally.
<center><math> |\psi(x,t)|^2 = \frac{e^{-x^2/(2 a^2 (1 + v_s^2 t^2))}}{[2 \pi a^2(1+ v_s^2 t^2)]^{1/2}}  </math></center>
Hence, the particle spreads on a distance
<center><math> \sqrt{\langle x^2 \rangle} \equiv \left( \int_{-\infty}^\infty d x x^2 |\psi(x,t)|^2\right)^{1/2} = a \sqrt{1+ v_s^2 t^2}
</math></center>


===Localization of the packet: General idea and experiment===
In the experiment shown here, a Gaussian packet of a Bose–Einstein condensate is prepared in a harmonic trap. When the trap is removed the cloud initially starts to expand but quickly stops spreading and remains localized.


 
To understand this behaviour we must solve the eigenstate problem of the Hamiltonian in the presence of disorder.
[[File:Localization1DB.png|thumb|left|x240px|BCE condensate confined with an harmonic trap in a 1D disordered geometry. Billy  et al, Nature 2008 ]]


In 1958, P. Anderson proposed that disorder can quench the dynamics of non ineracting particles. This phenomenon, called Anderson localization, is very well observed in numerical simulations. Recenlty, it has been clearly observed in experiments. Here, a Gaussian packet of a BCE is prepared with an harmonic trap. When the trap is removed,  particles start to spread, but remain quickly localized in space.
[[File:Localization1dB.png|thumb|left|x140px|Semilog plot of the particle density. Billy et al., Nature (2008).]]


To explain how this is possible we need to solve the eigenstate problem of the hamiltonian with a disorder potential <math> V(x)</math>, e.g. a Gaussian white noise. We will study this problem in different dimensions and see that eigenstates  can be localized in space.
In a disordered potential an eigenstate of energy <math>E_k</math> has the form


<math display="block">
\psi_k(x,t)=\psi_k(x)e^{-iE_k t/\hbar}.
</math>


[[File:Localization1dB.png|thumb|left|x140px|Semilog plot of the particle concentration.
The spatial part of the wavefunction is localized around some position <math>\bar{x}</math> and decays exponentially:
Billy  et al, Nature 2008 ]]
 
<math display="block">
\psi_k(x)\sim
e^{-|x-\bar{x}|/\xi_{\text{loc}}(E_k)}.
</math>
 
Here <math>\xi_{\text{loc}}</math> is the localization length.
 
Because the eigenstates are localized, each of them has support only in a finite spatial region. The wave packet is therefore built from eigenstates localized near its initial position.
 
Two important observations follow:
 
* eigenstates localized far from the packet contribute exponentially little to the wavefunction near the packet,
* eigenstates contributing to the packet decay exponentially far from their center.
 
As a consequence, transport far away from the initial position of the particle is '''exponentially suppressed'''.
 
Localization therefore does not mean that the particle is frozen. The wave packet can still evolve in time, but it cannot propagate arbitrarily far.
 
== Conductance and diffusive transport ==
 
In most materials the effect of weak disorder is not localization but diffusion.
 
In the Drude picture electrons scatter randomly on impurities. Beyond the mean free path the motion of electrons becomes diffusive.
 
In this regime Ohm's laws hold.
 
* First law
 
<math display="block">
\frac{V}{I}=R,
\qquad
\frac{I}{V}=G.
</math>


In this case an eigenstate <math> \psi_k(x,t) </math> of energy <math>E_k</math>  is concentrated around some position <math> \overline{x} </math>. Far from its center, the wave function decays exponentially with a finite localization length:  <math>\xi_{\text{loc}}(E_k)</math> (this length has a mild dependence on <math>E_k</math> that we  neglect for the moment).
Here <math>R</math> is the resistance and <math>G</math> the conductance.
<center><math> \psi_k(x,t) = \psi_k(x) e^{- i E_k t/\hbar},\quad \text{with}  \; \psi_k(x) \sim e^{-\frac{|x-\overline{x}|}{\xi_{\text{loc}}(E_k)}}  </math></center>
In the experiment an exponential decay of the Gaussian packet has been clearly shown.
How it is related with the exponential decay of the eigenstates?  To describe the evolution of the packet, we  proceed as before:
* First, decompose the wave packet as a superposition on the new eigenstates. The coefficients <math> c(k) </math> are concentrated only on the eigenstates which are localized around the packet.
* Second, evolve the packet in a given point as the superposition of the time evolved eigenstates in that point.
* It is time for two crucial observations: Eigenstates localized far away from the packet provide an  exponetially small contribution it. Viceversa, the eigenstates on which the packet is decomposed provide an exponentially small contribution far away from it.
As a consequence, transport far away from the initial location of the particle is <Strong> exponentially</Strong> suppressed.


== The conductance and the diffusive behaviour  ==
The  main effect of disorder in real material is not localization, but is that transport is diffusive. Indeed, from the time of Drude model, we know  that electrons scatter on the impurities with a random velocity. Beyond the mean free path an electron perform a random walk in the material. We will no go into that, but we will assume that if the effect of the disorder is weak, the material is a <Strong> conductor</Strong> and obeys Ohm's laws:
* First law:
<center><math> \frac{V}{I}= R, \quad \text{or} \quad \frac{I}{V}= G  </math></center>
Here <math>R </math> is the resistence of the sample and <math> G </math> is its conductance.
* Second law
* Second law
<center><math>  R =\rho \frac{L}{S}\sim \rho L^{2-d} , \quad \text{or} \quad  G = \sigma \frac{S}{L}\sim \sigma L^{d-2}  </math></center>
Here <math>  \rho, \sigma</math> are the resistivity and the conductivity. These are material properties, independent of the geometry of the sample. Let us stress once again that these  two phenomenological laws are a macroscpic manifestation of the diffusive motion of the electrons.


==Conductance in presence of a localized spectrum==
<math display="block">
R=\rho\frac{L}{S}\sim \rho L^{2-d},
\qquad
G=\sigma\frac{S}{L}\sim \sigma L^{d-2}.
</math>
 
Here <math>\rho</math> and <math>\sigma</math> are the resistivity and conductivity. These are material properties independent of the geometry of the sample.


If the effect of the disorder is strong also diffusion is suppressed. The conductance in presence of localization decay exponentially
These phenomenological laws are the macroscopic manifestation of diffusive transport.
<center><math>  G \sim e^{-2 L/\xi_{\text{loc}}}  </math></center>
There are different ways to derive this result here we use the Landauer model.


We present the idea <math>d=1</math> and consider a problem in which disorder is present only in a box of size <math>(0,L)</math>.
== Conductance in the localized regime ==


* Construct the <Strong> scattering states </Strong> of energy <math>E=\frac{\hbar^2 k^2}{2 m}</math>. Far from the box the states are plane waves coming from the left:
When disorder is strong, diffusion is suppressed and the system becomes insulating.
<center><math>
 
\psi_{k,L}(x)=\begin{cases}(
In the localized phase the conductance decays exponentially with the system size:
e^{ikx} + r e^{-ikx})\quad & \text{for} \;\; x \to - \infty\\
 
t e^{ikx} \quad  & \text{for} \;\; x \to \infty
<math display="block">
\end{cases}
G\sim e^{-2L/\xi_{\text{loc}}}.
</math></center>
</math>
or coming from the right
 
<center><math>
The localization length <math>\xi_{\text{loc}}</math> characterizes the spatial decay of the eigenstates.
\psi_{k,R}(x)=\begin{cases}(
 
e^{-ikx} + r e^{ikx})\quad & \text{for} \;\; x \to \infty\\
== The “Gang of Four” scaling theory ==
t e^{-ikx} \quad  & \text{for} \;\; x \to - \infty
 
\end{cases}
In a famous PRL (1979), Abrahams, Anderson, Licciardello and Ramakrishnan proposed a scaling theory of localization.
</math></center>


* The current probability of a wave function in 1D writes
It is based on the idea that the relevant quantity is the dimensionless conductance
<center><math>
J = \frac{\hbar}{2 m i}\left( \psi^* \frac{d\psi}{dx} - \psi \frac{d\psi^*}{dx}\right)
</math></center>
Hence you can derive
<center><math>
J_{k, L}=-J_{-k,R}=\frac{\hbar k}{m} |t|^2
</math></center>


* In the Laundauer model we apply a voltage <math> V</math> and meassure the current as the sum of <math> J_{k, L}</math>  up to the energy    <math> E_F+ e V</math> and of    <math> J_{k, L}</math> up to <math> E_F</math>. Pay attention to the Jacobian <math> \frac{\hbar^2 k}{  m} dk = dE </math> and get:
<math display="block">
<center><math>
g=\frac{G\hbar}{e^2}.
I= e \int_0^{E_F+ e V} \frac{ m d E}{\hbar^2 k} J_{k,L} + e \int_0^{E_F} \frac{ m d E}{\hbar^2 k} J_{-k,R} \sim \frac{e^2}{\hbar} |t|^2 V
</math>
</math></center>
Conclude
<center><math>
G = \frac{e^2}{\hbar} |t|^2 
</math></center>
* If the sample is in the localized phase  we can expect
<center><math>
|t|^2 = e^{-\frac{2 L}{\xi_{\text{loc}}}} 
</math></center>


== The gang of the four==
The scaling equation reads


In a famous PRL (1979) four physicists, Abrahams, Anderson, Licciardello, Ramakrishnan, proposed a scaling theory for localization. It provides a valid argument on the fate of localization at different dimensions. It is based on two well justified assumptions:
<math display="block">
* <Strong> Single parameter scaling</Strong> (due to Thouless): the RG equation of the dimensionless conductance <math> g= G \hbar/ e^2 </math> is a function of the dimension and of <math> g </math> only:
\frac{d\ln g}{d\ln L}=\beta(g).
<center>
<math>
\frac{d \ln g}{d \ln L}= \beta(g)
</math>
</math>
</center>
 
The beta function is unknown and the <math> \ln g </math> mean we are interested in the typical behaviour instead of the mean behaviour.
The function <math>\beta(g)</math> depends only on <math>g</math> and the spatial dimension.
* <Strong> Asymptotic behaviour</Strong>: For large conductance we can assume that the material (at these scales) is a conductor, for very small conductance we take instead the insulating behaviour. Hence
 
<center><math>
The asymptotic behaviours are
\beta(g)=\begin{cases}
 
d-2 \quad & \text{for} \;\; \ln g \to \infty\\
<math display="block">
-\frac{2 L}{\xi_{\text{loc}}} =\ln g \quad  & \text{for} \;\; \ln g \to -\infty
\beta(g)=
\begin{cases}
d-2 & g\to\infty \\
\sim \ln g & g\to 0
\end{cases}
\end{cases}
</math></center>
</math>
 
The second relation reflects the exponential suppression of conductance in the localized regime.
 
If the beta function is monotonic, the scaling theory predicts:


The simplest scenario is that the beta function is monotonous. Within this picture, we expect a single metal-insulator transition for <math>d>2</math> and a fully localized phase at <math>d=2</math> and below. More complex scenarios, with multiple transitions go against physical intuition.
* a metal–insulator transition for <math>d>2</math>,
* complete localization for <math>d\le2</math>.

Revision as of 19:43, 5 March 2026

Goal. This lecture introduces the phenomenon of localization. Localization is a wave phenomenon induced by disorder that suppresses transport in a system.

Free particles and ballistic behaviour

The Schrödinger equation governs the evolution of the quantum state of a particle in one dimension:

itψ(x,t)=Hψ(x,t),H=22md2dx2+V(x).

Here H is the Hamiltonian. For a free particle the potential vanishes, V(x)=0.

One first looks for separable solutions, also called eigenstates or stationary solutions. If the particle is in an eigenstate, all physical observables are time independent.

For a free particle the stationary solutions are plane waves

ψk(x,t)=12πeikxeiEkt/,Ek=2k22m.

Here k is a real number and the spectrum is continuous. These solutions are completely delocalized over the real axis and cannot be normalized.

Physical states are obtained as superpositions of eigenstates. By the superposition principle any linear combination of stationary solutions is again a solution of the Schrödinger equation. Hence we can construct a localized wave packet:

ψ(x,t)=dkc(k)ψk(x,t),dk|c(k)|2=1.

Evolution of a Gaussian wave packet

  • Initial state.

At time t=0 consider the Gaussian packet

ψ(x,0)=ex2/(4a2)(2πa2)1/4,|ψ(x,0)|2=ex2/(2a2)2πa2.

Show that the coefficients of the plane wave decomposition are

c(k)=(2a2π)1/4ea2k2.

  • Time evolution.

Define the spreading velocity

vs=2ma2.

Show that the time evolution of the packet is

ψ(x,t)=ex2/(4a2(1+ivst))[2πa2(1+ivst)]1/4.

  • Ballistic spreading.

The probability density becomes

|ψ(x,t)|2=ex2/(2a2(1+vs2t2))2πa2(1+vs2t2).

Hence

x2=(dxx2|ψ(x,t)|2)1/2=a1+vs2t2.

At long times

x2vst.

This behaviour is called ballistic spreading.

It should be contrasted with the two other possible transport regimes:

  • Diffusive motion (random walk):

x2t

  • Localized regime:

x2 saturates at long times.

Understanding how disorder changes ballistic motion into diffusion and eventually localization is the main goal of the following lectures.

Localization of the packet: general idea and experiment

BEC condensate expanding in a 1D disordered potential. Billy et al., Nature (2008).

In 1958, P. Anderson proposed that disorder can suppress transport in quantum systems. This phenomenon, called Anderson localization, has since been observed both numerically and experimentally.

In the experiment shown here, a Gaussian packet of a Bose–Einstein condensate is prepared in a harmonic trap. When the trap is removed the cloud initially starts to expand but quickly stops spreading and remains localized.

To understand this behaviour we must solve the eigenstate problem of the Hamiltonian in the presence of disorder.

Semilog plot of the particle density. Billy et al., Nature (2008).

In a disordered potential an eigenstate of energy Ek has the form

ψk(x,t)=ψk(x)eiEkt/.

The spatial part of the wavefunction is localized around some position x¯ and decays exponentially:

ψk(x)e|xx¯|/ξloc(Ek).

Here ξloc is the localization length.

Because the eigenstates are localized, each of them has support only in a finite spatial region. The wave packet is therefore built from eigenstates localized near its initial position.

Two important observations follow:

  • eigenstates localized far from the packet contribute exponentially little to the wavefunction near the packet,
  • eigenstates contributing to the packet decay exponentially far from their center.

As a consequence, transport far away from the initial position of the particle is exponentially suppressed.

Localization therefore does not mean that the particle is frozen. The wave packet can still evolve in time, but it cannot propagate arbitrarily far.

Conductance and diffusive transport

In most materials the effect of weak disorder is not localization but diffusion.

In the Drude picture electrons scatter randomly on impurities. Beyond the mean free path the motion of electrons becomes diffusive.

In this regime Ohm's laws hold.

  • First law

VI=R,IV=G.

Here R is the resistance and G the conductance.

  • Second law

R=ρLSρL2d,G=σSLσLd2.

Here ρ and σ are the resistivity and conductivity. These are material properties independent of the geometry of the sample.

These phenomenological laws are the macroscopic manifestation of diffusive transport.

Conductance in the localized regime

When disorder is strong, diffusion is suppressed and the system becomes insulating.

In the localized phase the conductance decays exponentially with the system size:

Ge2L/ξloc.

The localization length ξloc characterizes the spatial decay of the eigenstates.

The “Gang of Four” scaling theory

In a famous PRL (1979), Abrahams, Anderson, Licciardello and Ramakrishnan proposed a scaling theory of localization.

It is based on the idea that the relevant quantity is the dimensionless conductance

g=Ge2.

The scaling equation reads

dlngdlnL=β(g).

The function β(g) depends only on g and the spatial dimension.

The asymptotic behaviours are

β(g)={d2glngg0

The second relation reflects the exponential suppression of conductance in the localized regime.

If the beta function is monotonic, the scaling theory predicts:

  • a metal–insulator transition for d>2,
  • complete localization for d2.