T-9

Goal: the goal of this problem is to determine when the solution of the distributional equations corresponding to localization is unstable, providing an estimate of thee mobility edge on the Bethe lattice.
Techniques: stability analysis, Laplace transforms.

Problems

In this Problem we determine for which values of parameters localization is stable, estimating the critical value of disorder where the transition to a delocalised phase occurs. Recall the results of Problem 8: the real and imaginary parts of the local self energy satisfy the self-consistent equations:

${\mathrm{\Gamma }}_a=\sum _{b\in \partial a}{t}_{ab}^2\frac{{\mathrm{\Gamma }}_b+\eta }{(E-W{V}_b-R_b{)}^2+({\mathrm{\Gamma }}_b+\eta {)}^2},R_a=\sum _{b\in \partial a}{t}_{ab}^2\frac{E-W{V}_b-R_b}{(E-W{V}_b-R_b{)}^2+({\mathrm{\Gamma }}_b+\eta {)}^2}$

These equations admit the solution ${\mathrm{\Gamma }}_a={\mathrm{\Gamma }}_b=0$ when $\eta =0$, which corresponds to localization. We now determine when this solution becomes unstable.

Problem 9: an estimate of the mobility edge

1. Imaginary approximation and distributional equation. We consider the equations for ${\mathrm{\Gamma }}_a$ and neglect the terms $R_b$ in the denominators, which couple the equations to those for the real parts of the self energies (“imaginary” approximation). Moreover, we assume to be in the localized phase, where ${\mathrm{\Gamma }}_a\sim \eta \ll 1$. Finally, we set $t_{ab}\equiv t$ and $E=0$ for simplicity. Show that under these assumptions the probability density for the imaginary part, $P_{\mathrm{\Gamma }}(\mathrm{\Gamma })$, satisfies

   $P_{\mathrm{\Gamma }}(\mathrm{\Gamma })=\int \prod _{b=1}^{k}d{ϵ}_{b}p({ϵ}_b)\int \prod _{b=1}^{k}d{\mathrm{\Gamma }}_b{P}_{\mathrm{\Gamma }}({\mathrm{\Gamma }}_b)\delta (\mathrm{\Gamma }-t^2\sum _{b\in \partial a}\frac{{\mathrm{\Gamma }}_b+\eta }{W^2{V}_b^2})$

   Show that the Laplace transform of this distribution, $\mathrm{\Phi }(s)=\underset{0}{\overset{\mathrm{\infty }}{\int }}d\mathrm{\Gamma }{e}^{-s\mathrm{\Gamma }}{P}_{\mathrm{\Gamma }}(\mathrm{\Gamma })$, satisfies

   $\mathrm{\Phi }(s)={[\int dϵp(ϵ)e^{-\frac{s{V}^2\eta }{{ϵ}^2}}\mathrm{\Phi }\left(\frac{s{V}^2}{{ϵ}^2}\right)]}^k$



2. The stability analysis. We now assume to be in the localized phase, when for $\eta \to 0$ the distribution $P_{\mathrm{\Gamma }}(\mathrm{\Gamma })\to \delta (\mathrm{\Gamma })$. We wish to check the stability of our assumption. This is done by controlling the tails of the distribution $P_{\mathrm{\Gamma }}(\mathrm{\Gamma })$ for finite $\eta$.
   • For finite $\eta$, we expect that typically ${\mathrm{\Gamma }}_a\sim \eta \ll 1$, and thus $P_{\mathrm{\Gamma }}(\mathrm{\Gamma })$ should have a peak at this scale; however, we also expect some power law decay $P_{\mathrm{\Gamma }}(\mathrm{\Gamma })\sim {\mathrm{\Gamma }}^{-\alpha }$ for large $\mathrm{\Gamma }$. These tails are contributed by the events in which one of the descendants $b$ is such that ${ϵ}_b$ is very small, in which case ${\mathrm{\Gamma }}_a\sim 1/{ϵ}_b^2\gg 1$. Show that if $\mathrm{\Gamma }\sim 1/{ϵ}^2$ and $p(ϵ)$ is not gapped around zero, then $P_{\mathrm{\Gamma }}(\mathrm{\Gamma })\sim {\mathrm{\Gamma }}^{-3/2}$, i.e. the distribution has tails contributed by these events in which the local random potential happen to be very small.
   • Assume more generally that $P_{\mathrm{\Gamma }}(\mathrm{\Gamma })\sim {\mathrm{\Gamma }}^{-\alpha }$ for large $\mathrm{\Gamma }$ and $\alpha \in (1,3/2]$. Show, using a dimensional analysis argument, that this corresponds to a non-analytic behaviour of the Laplace transform, $\mathrm{\Phi }(s)\sim 1-A|s{|}^{\beta }$ for $s$ small, with $\beta =\alpha -1\in (0,1/2]$.
   • Show that the equation for $\mathrm{\Phi }(s)$ gives for $s$ small $1-A{s}^{\beta }=1-Ak\int dϵp(ϵ)\frac{s^{\beta }{V}^{2\beta }}{{ϵ}^{2\beta }}+o(s^{\beta })$, and therefore this is consistent provided that there exists a $\beta \in (0,1/2]$ solving

     $1=k\int dϵp(ϵ){\left(\frac{V}{|ϵ|}\right)}^{2\beta }\equiv kI(\beta ).$






Behaviour of the integral $I(\beta )$ in the case of uniformily distributed disorder, for $W<W_c$ .

3. The critical disorder. Consider now local fields $ϵ$ taken from a uniform distribution in $[-W/2,W/2]$. Compute $I(\beta )$ and show that it is non monotonic, with a local minimum ${\beta }^{*}$ in the interval $[0,1/2]$. Show that $I({\beta }^{*})$ increases as the disorder is made weaker and weaker, and thus the transition to delocalisation occurs at the critical value of disorder when $I({\beta }^{*})=k^{-1}$. Show that this gives the following estimate for the critical disorder at which the transition to delocalisation occurs:

   $\frac{W_c}{V}=2ke\log \left(\frac{W_c}{2V}\right)\sim 2ek\log \left(k\right)$

   Why the critical disorder increases with $k$?

Check out: key concepts of this TD

Stability analysis, critical disorder, mobility edge.

References

• Abou-Chacra, Thouless, Anderson. A selfconsistent theory of localization. Journal of Physics C: Solid State Physics 6.10 (1973)