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

We now focus on the self energies, since the criterion for localization is given in terms of these quantities. In this Problem we will determine for which values of parameters localization is stable, estimating the critical value of disorder where the transition to a delocalised phase occurs.

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 _{c=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 }{{ϵ}_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{t}^2\eta }{{ϵ}^2}}\mathrm{\Phi }\left(\frac{s{t}^2}{{ϵ}^2}\right)]}^k$



2. The stability analysis. We now wish to check the stability of our assumption to be in the localized phase, ${\mathrm{\Gamma }}_a\sim \eta \ll 1$, which led to the identity above for $\mathrm{\Phi }(s)$. Our assumption is that the typical value of ${\mathrm{\Gamma }}_a$ is small, except for cases in which one of the descendants $b$ is such that ${ϵ}_b$ is very small, in which case ${\mathrm{\Gamma }}_a\sim 1/{ϵ}_b^2$.
   • 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 fields 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 that this corresponds to $\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 }{t}^{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{t}{ϵ}\right)}^{2\beta }\equiv kI(\beta ).$



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

   $W_c=t2ke\log \left(\frac{W_c}{2t}\right)\sim t2ke\log \left(k\right)$

   Why the critical disorder increases with $k$?

References

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