L-8

Goal: we will introduce the Anderson model, discuss the behaviour as a function of the dimension. In 1d localization can be connected to the product of random matrices.

Anderson model (tight binding model)

We consider disordered non-interacting particles hopping between nearest neighbors sites on a lattice. The hamiltonian reads:

$H=-t\sum _{<i,j>}(c_i^{†}{c}_j+c_j^{†}{c}_i)\sum _i{V}_i{c}_i^{†}{c}_i$

The single particle hamiltonian in 1d reads

$H=\left[\begin{array}{cccccc}{V}_1& -t& 0& 0& 0& 0\\ -t& V_2& -t& 0& 0& 0\\ 0& -t& V_3& -t& 0& 0\\ 0& 0& -t& \dots & -t& 0\\ 0& 0& 0& -t& \dots & -t\\ 0& 0& 0& 0& -t& V_L\end{array}\right]$

For simplicity we set the hopping $t=1$. The disorder are iid random variables drawn, uniformly from the box $(-\frac{W}{2},\frac{W}{2})$.

The final goal is to study the statistical properties of eigensystem

$H\psi =ϵ\psi ,\text{with}\sum _n|{\psi }_n{|}^2=1$

Density of states (DOS)

In 1d and in absence of disorder, the dispersion relation is $ϵ(k)=-2\cos k,k\in (-\pi ,\pi ),-2<ϵ(k)<2$. From the dispersion relation, we compute the density of states (DOS) :

$\rho (ϵ)=\underset{-\pi }{\overset{\pi }{\int }}\frac{dk}{2\pi }\delta (ϵ-ϵ(k))=\frac{1}{\pi }\frac{1}{\sqrt{4-{ϵ}^2}}\text{for }ϵ\in (-2,2)$

In presence of disorder the DOS becomes larger, and display sample to sample fluctuations. One can consider its mean value, avergaed over disorder realization.

Eigenstates

In absence of disorder the eigenstates are plane waves delocalized along all the system. In presence of disorder, three situations can occur and to distinguish them it is useful to introduce the inverse participation ratio, IPR

$IPR(q)=\sum _n|{\psi }_n{|}^{2q}\sim L^{-{\tau }_q}$

The normalization imposes ${\tau }_1=0$. For $q=0$, $|{\psi }_n{|}^2=1$, hence, ${\tau }_0=-d$.

• Delocalized eigenstates In this case, $|{\psi }_n{|}^2\approx L^{-d}$. Hence, we expect

$IPR(q)=L^{d(1-q)}{\tau }_q=d(1-q)$

• Localized eigenstates In this case, $|{\psi }_n{|}^2\approx 1/{\xi }_{\text{loc}}^{1/d}$ for ${\xi }_{\text{loc}}^d$ sites and almost zero elsewhere. Hence, we expect

$IPR(q)=\text{const},{\tau }_q=0$

• Multifractal eigenstates. At the transition( the mobility edge) an anomalous scaling is observed:

$IPR(q)=L^{D_q(1-q)}{\tau }_q=D_q(1-q)$

Here $D_q$ is q-dependent multifractal dimension, smaller than $d$ and larger than zero.

Transfer matrices and Lyapunov exponents

Product of random variables and Central limit theorem

Consider a set of positive iid random variables $x_1,x_2,\dots x_N$ with finite mean and variance and compute their product

${\mathrm{\Pi }}_N=\prod _{n=1}^N{x}_i,\text{or}\ln {\mathrm{\Pi }}_N=\sum _{n=1}^N\ln x_i$

For large N, the Central Limit Theorem predicts:

$\log {\mathrm{\Pi }}_N=\gamma N+{\gamma }_2\sqrt{N}\chi$

• $\chi$ is a Gaussian number of zero mean and unit variance
• $\gamma ,{\gamma }_2$ are N indepent and can be written as

$\gamma =\overline{\ln x}<\ln \bar{x},{\gamma }_2=\sqrt{\overline{(\ln x{)}^2}-(\overline{\ln x}{)}^2}$

Log-normal distribution

The distribution of ${\mathrm{\Pi }}_N$ is log-normal

$P({\mathrm{\Pi }}_N)d{\mathrm{\Pi }}_N=\frac{1}{\sqrt{2\pi {\gamma }_2^{2}N}}\exp [-\frac{(\ln ({\mathrm{\Pi }}_N)-\gamma N{)}^2}{2{\gamma }_2^{2}N}]\frac{d{\mathrm{\Pi }}_N}{{\mathrm{\Pi }}_N}$

Quenched and Annealed averages

To compute the moments of the log-normal distribution, it is convenient to introduce the variable

$X\equiv \ln ({\mathrm{\Pi }}_N)$

which is Gaussian distributed:

$p(X)=\frac{1}{\sqrt{2\pi {\sigma }^2}}\exp [-\frac{(X-\mu {)}^2}{2{\sigma }^2}]$

with $\mu =\gamma N$ and ${\sigma }^2={\gamma }_2^{2}N$.

The moments of

${\mathrm{\Pi }}_N$

can be easily computed:

$\overline{{\mathrm{\Pi }}_N^n}=\int dX{e}^{nX}p(X)=\exp [\mu n+{\sigma }^2\frac{n^2}{2}]=\exp [(\gamma n+{\gamma }_2^2\frac{n^2}{2})N]$

The variable ${\mathrm{\Pi }}_N$ is therefore not self-averaging (see Valentina's lecture 1) since its fluctuations grow with $N$ faster than its mean:

$\frac{\overline{{\mathrm{\Pi }}_N^2}}{(\overline{{\mathrm{\Pi }}_N}{)}^2}=\exp \left[{\gamma }_2^{2}N\right]$

Hence, ${\mathrm{\Pi }}_N$ is not self-averaging, while $\ln {\mathrm{\Pi }}_N$ is self-averaging.

In particular, the mean $\overline{{\mathrm{\Pi }}_N}=\exp [(\gamma +{\gamma }_2^2/2)N]$ grows much faster than the typical value ${\mathrm{\Pi }}_N^{\text{typ}}\equiv \exp (\gamma N)$.

Product of random matrices

Let's consider again the Anderson Model in 1d. The eigensystem is well defined in a box of size L with Dirichelet boundary condition on the extremeties of the box.

Here we will solve the second order differential equation imposing instead Cauchy boundaries on one side of the box. Let's rewrite the previous eigensystem in the following form

$\left[\begin{array}{c}{\psi }_{n+1}\\ {\psi }_n\end{array}\right]=\left[\begin{array}{cc}{V}_n-ϵ& -1\\ 1& 0\end{array}\right]\left[\begin{array}{c}{\psi }_n\\ {\psi }_{n-1}\end{array}\right]$

We can continue the recursion

$\left[\begin{array}{c}{\psi }_{n+1}\\ {\psi }_n\end{array}\right]=\left[\begin{array}{cc}{V}_n-ϵ& -1\\ 1& 0\end{array}\right]\left[\begin{array}{cc}{V}_{n-1}-ϵ& -1\\ 1& 0\end{array}\right]\left[\begin{array}{c}{\psi }_{n-1}\\ {\psi }_{n-2}\end{array}\right]$

It is useful to introduce the transfer matrix and their product

$T_n=\left[\begin{array}{cc}{V}_n-ϵ& -1\\ 1& 0\end{array}\right],{\mathrm{\Pi }}_n=T_n\cdot T_{n-1}\cdot \dots T_1$

The Schrodinger equation can be written as

$\left[\begin{array}{c}{\psi }_{n+1}\\ {\psi }_n\end{array}\right]={\mathrm{\Pi }}_n\left[\begin{array}{c}{\psi }_1\\ {\psi }_0\end{array}\right]=\left[\begin{array}{cc}{\pi }_{11}& {\pi }_{12}\\ {\pi }_{21}& {\pi }_{22}\end{array}\right]\left[\begin{array}{c}{\psi }_1\\ {\psi }_0\end{array}\right]$

Fustenberg Theorem

We define the norm of a 2x2 matrix:

$\Vert {\mathrm{\Pi }}_n{\Vert }^2=\frac{{\pi }_{11}^2+{\pi }_{21}^2+{\pi }_{12}^2+{\pi }_{22}^2}{2}$

For large N, the Fustenberg theorem ensures the existence of the non-negative Lyapunov exponent, namely

$\underset{n\to \mathrm{\infty }}{\lim }\frac{\ln \Vert {\mathrm{\Pi }}_n\Vert }{n}=\gamma \ge 0$

In absence of disorder $\gamma =0$ for $ϵ\in (-2,2)$. Generically the Lyapunov is positive, $\gamma >0$, and depends on $ϵ$ and on the distribution of $V_i$.

Consequences

Localization length

Together with the norm, also $|{\psi }_n{|}^2$ grows exponentially with n. We can write

$\ln |{\psi }_n|\sim \gamma n+{\gamma }_2\chi \sqrt{n}$

which means that $\ln |{\psi }_n|$ is performing a random walk with a drift.

However, our initial goal is a properly normalized eigenstate at energy $ϵ$. We need to switch from Cauchy, where you set the initial condition, to Dirichelet or vonNeuman, where you set the behaviour at the two boundaries. The true eigenstate is obtained by matching two "Cauchy" solutions on the half box and imposing the normalization. Hence, we obtain a localized eigenstate and we can identify

${\xi }_{\text{loc}}(ϵ)\equiv {\gamma }^{-1}(ϵ)$

Fluctuations

We expect strong fluctuations on quantites like $|{\psi }_n|,\Vert {\mathrm{\Pi }}_n\Vert ,G,\dots$, while their logarithm is self averaging.