L-9: Difference between revisions

From Disordered Systems Wiki
Jump to navigation Jump to search
No edit summary
Line 3: Line 3:
Without disorder, the eigenstates are delocalized plane waves.
Without disorder, the eigenstates are delocalized plane waves.


In the presence of disorder, three scenarios can arise: delocalized eigenstates, where the wavefunction remains extended; localized eigenstates, where the wavefunction is exponentially confined to a finite region; multifractal eigenstates, occurring at the mobility edge, where the wavefunction exhibits a complex, scale-dependent structure.
In the presence of disorder, three scenarios can arise:
 
* '''Delocalized eigenstates''', where the wavefunction remains extended over the whole system.
* '''Localized eigenstates''', where the wavefunction is exponentially confined to a finite region.
* '''Multifractal eigenstates''', occurring at the '''mobility edge''' of the Anderson transition, where the wavefunction exhibits a complex, scale–dependent structure.
 
In three dimensions, increasing the disorder strength leads to a transition between a metallic phase with delocalized eigenstates and an insulating phase with localized eigenstates. The critical energy separating these two regimes is called the '''mobility edge'''. Eigenstates at the mobility edge display multifractal statistics.
 
To characterize these different regimes it is useful to introduce the '''inverse participation ratio''' (IPR)
 
<math display="block">
\mathrm{IPR}(q)=\sum_n |\psi_n|^{2 q}.
</math>
 
For a system of linear size <math>L</math>, one typically observes a scaling


To distinguish these regimes, it is useful to introduce the inverse participation ratio (IPR).
<math display="block">
<math display="block">
\mathrm{IPR}(q)=\sum_n |\psi_n|^{2 q} \sim L^{-\tau_q}.
\mathrm{IPR}(q) \sim L^{-\tau_q}.
</math>
</math>
The exponent <math>\tau_q</math> characterizes the spatial structure of the eigenstate.


== Delocalized eigenstates ==
== Delocalized eigenstates ==


In this case, <math>|\psi_n|^{2} \approx L^{-d}</math>. Hence, we expect
In a delocalized state the probability is spread uniformly over the system:
 
<math display="block">
|\psi_n|^{2} \approx L^{-d}.
</math>
 
Therefore
 
<math display="block">
\mathrm{IPR}(q)
=
\sum_n |\psi_n|^{2q}
\sim L^d (L^{-d})^q
=
L^{d(1-q)}.
</math>
 
Thus
 
<math display="block">
<math display="block">
\mathrm{IPR}(q)=L^{d(1-q)},
\quad
\tau_q=d(1-q).
\tau_q=d(1-q).
</math>
</math>
Line 21: Line 52:
== Localized eigenstates ==
== Localized eigenstates ==


In this case, <math>|\psi_n|^{2} \approx 1/\xi_{\text{loc}}^{d}</math> on <math>\xi_{\text{loc}}^{d}</math> sites and almost zero elsewhere. Hence, we expect
In a localized state the wavefunction is concentrated within a region of size <math>\xi_{\text{loc}}</math>.
 
Roughly
 
<math display="block">
<math display="block">
\mathrm{IPR}(q)=\text{const},
|\psi_n|^2 \sim \xi_{\text{loc}}^{-d}
\quad
</math>
 
on <math>\xi_{\text{loc}}^d</math> sites and negligible elsewhere. Hence
 
<math display="block">
\mathrm{IPR}(q)\sim \text{const},
\qquad
\tau_q=0.
\tau_q=0.
</math>
</math>
Thus localized states do not scale with system size.


== Multifractal eigenstates ==
== Multifractal eigenstates ==


The exponent <math>\tau_q</math> is called the '''multifractal exponent'''. It is a non-decreasing function of <math>q</math> with some special points:
At the mobility edge the eigenstates are neither extended nor localized. Instead, the amplitudes fluctuate strongly across the system.
 
The scaling of the IPR is characterized by a '''nonlinear function'''
 
<math display="block">
\tau_q.
</math>
 
The corresponding '''multifractal dimensions''' are defined as
 
<math display="block">
D_q=\frac{\tau_q}{q-1}.
</math>
 
Multifractality means that different moments of the wavefunction probe different effective dimensions.
 
== Multifractal spectrum ==
 
To describe this structure it is useful to introduce the '''multifractal spectrum''' <math>f(\alpha)</math>.


* <math>\tau_0 = -d</math>, since the wavefunction is defined on all sites. In general, <math>\tau_0</math> represents the fractal dimension of the object under consideration and is purely a geometric property.
We assume that the wavefunction amplitudes scale as
* <math>\tau_1 = 0</math>, imposed by normalization.


To observe multifractal behavior, we expect:
<math display="block">
<math display="block">
|\psi_n|^{2} \approx L^{-\alpha}
|\psi_n|^2 \sim L^{-\alpha}
\quad \text{for}\;
L^{f(\alpha)} \; \text{sites}.
</math>
</math>


The exponent <math>\alpha</math> is positive, and <math>f(\alpha)</math> is called the '''multifractal spectrum'''. Its maximum corresponds to the fractal dimension of the object, which in our case is <math>d</math>.
on approximately


The relation between the multifractal spectrum <math>f(\alpha)</math> and the exponent <math>\tau_q</math> is given by:
<math display="block">
<math display="block">
\mathrm{IPR}(q) = \sum_n |\psi_n|^{2q}
L^{f(\alpha)}
\sim \int d\alpha \, L^{-\alpha q} L^{f(\alpha)}
</math>
</math>
for large <math>L</math>. From this, we obtain:
 
sites.
 
The function <math>f(\alpha)</math> describes the fractal dimension of the set of sites where the wavefunction has exponent <math>\alpha</math>.
 
Using this representation,
 
<math display="block">
<math display="block">
\tau(q) = \min_{\alpha} (\alpha q - f(\alpha)).
\mathrm{IPR}(q)
=
\sum_n |\psi_n|^{2q}
\sim
\int d\alpha \,
L^{-\alpha q} L^{f(\alpha)} .
</math>
</math>


This implies that for <math>\alpha^*(q)</math>, which satisfies
For large <math>L</math>, the integral is dominated by the saddle point, giving
<math display="block">f'(\alpha^*(q)) = q,</math>
we have
<math display="block">\tau(q) = \alpha^*(q) q - f(\alpha^*(q)).</math>


'''Delocalized wavefunctions''' have a simple spectrum: for <math>\alpha = d</math>, we find <math>f(\alpha = d) = d</math> and <math>f(\alpha \neq d) = -\infty</math>. This means that <math>\alpha^*(q) = d</math> is independent of <math>q</math>.
<math display="block">
\tau(q)=\min_{\alpha}\left(\alpha q-f(\alpha)\right).
</math>
 
Thus <math>\tau(q)</math> and <math>f(\alpha)</math> are related by a '''Legendre transform'''.
 
If <math>\alpha^*(q)</math> satisfies
 
<math display="block">
f'(\alpha^*(q))=q,
</math>
 
then
 
<math display="block">
\tau(q)=\alpha^*(q)q-f(\alpha^*(q)).
</math>
 
Multifractal wavefunctions exhibit a smooth spectrum with maximum
 
<math display="block">
f(\alpha_0)=d,
</math>


'''Multifractal wavefunctions''' exhibit a smoother dependence, leading to a continuous spectrum with a maximum at <math>\alpha_0</math>, where <math>f(\alpha_0) = d</math>. At <math>q = 1</math>, we have <math>f'(\alpha_1) = 1</math> and <math>f(\alpha_1) = \alpha_1</math>.
indicating that the wavefunction explores the entire system but with strong spatial fluctuations.


= Larkin model =
= Larkin model =

Revision as of 17:08, 15 March 2026

Eigenstates

Without disorder, the eigenstates are delocalized plane waves.

In the presence of disorder, three scenarios can arise:

  • Delocalized eigenstates, where the wavefunction remains extended over the whole system.
  • Localized eigenstates, where the wavefunction is exponentially confined to a finite region.
  • Multifractal eigenstates, occurring at the mobility edge of the Anderson transition, where the wavefunction exhibits a complex, scale–dependent structure.

In three dimensions, increasing the disorder strength leads to a transition between a metallic phase with delocalized eigenstates and an insulating phase with localized eigenstates. The critical energy separating these two regimes is called the mobility edge. Eigenstates at the mobility edge display multifractal statistics.

To characterize these different regimes it is useful to introduce the inverse participation ratio (IPR)

IPR(q)=n|ψn|2q.

For a system of linear size L, one typically observes a scaling

IPR(q)Lτq.

The exponent τq characterizes the spatial structure of the eigenstate.

Delocalized eigenstates

In a delocalized state the probability is spread uniformly over the system:

|ψn|2Ld.

Therefore

IPR(q)=n|ψn|2qLd(Ld)q=Ld(1q).

Thus

τq=d(1q).

Localized eigenstates

In a localized state the wavefunction is concentrated within a region of size ξloc.

Roughly

|ψn|2ξlocd

on ξlocd sites and negligible elsewhere. Hence

IPR(q)const,τq=0.

Thus localized states do not scale with system size.

Multifractal eigenstates

At the mobility edge the eigenstates are neither extended nor localized. Instead, the amplitudes fluctuate strongly across the system.

The scaling of the IPR is characterized by a nonlinear function

τq.

The corresponding multifractal dimensions are defined as

Dq=τqq1.

Multifractality means that different moments of the wavefunction probe different effective dimensions.

Multifractal spectrum

To describe this structure it is useful to introduce the multifractal spectrum f(α).

We assume that the wavefunction amplitudes scale as

|ψn|2Lα

on approximately

Lf(α)

sites.

The function f(α) describes the fractal dimension of the set of sites where the wavefunction has exponent α.

Using this representation,

IPR(q)=n|ψn|2qdαLαqLf(α).

For large L, the integral is dominated by the saddle point, giving

τ(q)=minα(αqf(α)).

Thus τ(q) and f(α) are related by a Legendre transform.

If α*(q) satisfies

f(α*(q))=q,

then

τ(q)=α*(q)qf(α*(q)).

Multifractal wavefunctions exhibit a smooth spectrum with maximum

f(α0)=d,

indicating that the wavefunction explores the entire system but with strong spatial fluctuations.

Larkin model

In your homework you solved a toy model for the interface: th(r,t)=2h(r,t)+F(r). For simplicity, we assume Gaussian disorder F(r)=0, F(r)F(r)=σ2δd(rr).

You proved that:

  • the roughness exponent of this model is ζL=4d2 below dimension 4
  • the force per unit length acting on the center of the interface is f=σ/Ld
  • at long times the interface shape is

h(q)h(q)=σ2qd+2ζL.

In the real depinning model the disorder is, however, a non-linear function of h. The idea of Larkin is that this linearization is correct up to rf, the correlation length of the disorder along the h direction. This defines a Larkin length.

Indeed, from (h(r)h(0))2=ddqh(q)h(q)(1cos(qr))σ2r2ζL, you get (h(L)h(0))2=rf2,L=(rfσ)1/ζL.

Above this scale, roughness changes and pinning starts with a critical force fc=σLd/(2ζL).

In d=1 we have L=(rfσ)2/3.