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=1}^L|{\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.

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.