L-9: Difference between revisions

From Disordered Systems Wiki
Jump to navigation Jump to search
 
(11 intermediate revisions by the same user not shown)
Line 1: Line 1:
=Eigenstates=
= Eigenstates =
 
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:


To distinguish these regimes, it is useful to introduce the inverse participation ratio (IPR).
* '''Delocalized eigenstates''', where the wavefunction remains extended over the whole system.
<center><math>
* '''Localized eigenstates''', where the wavefunction is exponentially confined to a finite region.
IPR(q)=\sum_n |\psi_n|^{2 q} \sim L^{-\tau_q}
* '''Multifractal eigenstates''', occurring at the '''mobility edge''' of the Anderson transition, where the wavefunction exhibits a complex, scale–dependent structure.
</math></center>


== Delocalized eigenstates==
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.
In this case, <math>|\psi_n|^{2} \approx L^{-d} </math>. Hence, we expect
<center><math>
IPR(q)=L^{d(1-q)}  \quad \tau_q=d(1-q)
</math></center>


==Localized eigenstates==
To characterize these different regimes it is useful to introduce the inverse participation ratio (IPR)
In this case, <math>|\psi_n|^{2} \approx 1/\xi_{\text{loc}}^{1/d} </math> for <math>\xi_{\text{loc}}^{d}</math> sites and almost zero elsewhere. Hence, we expect
<center><math>
IPR(q)= \text{const},  \quad \tau_q=0
</math></center>


<math display="block">
\mathrm{IPR}(q)=\sum_n |\psi_n|^{2 q}.
</math>


==Multifractal eigenstates==
For a system of linear size <math>L</math>, one typically observes a scaling


The exponent <math>\tau_q</math> is called <Strong> multifractal exponent </Strong>. It is a non decreasing  function with q with some special points:
<math display="block">
* <math>\tau_0 =-d </math> because the wave fuction is defined  on all sites, in general <math>\tau_0 </math> is the fractal dimension of the object we are considering. It is simply a geometrical property.
\mathrm{IPR}(q) \sim L^{-\tau_q}.
* <math>\tau_1 =0 </math> imposed by normalization.
</math>


To have multifractal behaviour  we expect
The exponent <math>\tau_q</math> characterizes the spatial structure of the eigenstate.
<center><math>
|\psi_n|^{2} \approx L^{-\alpha} \quad  \text{for}\; L^{f(\alpha)} \; \text{sites} 
</math></center>
The exponent <math>\alpha </math> is positive and <math>f(\alpha)</math> is called <Strong> multifractal spectrum </Strong>. Its maximum is the fractal dimension of the object, in our case d. We can determine the relation between multifractal spectrum <math>f(\alpha)</math> and  exponent <math>\tau_q</math>
<center><math>
IPR(q)=\sum_n |\psi_n|^{2 q}\sim \int d \alpha L^{-\alpha q} L^{f(\alpha)}
</math></center>
for large L
<center><math>
\tau(q)= \min_{\alpha}{(\alpha q -f(\alpha))} 
</math></center>
This means that for <math>\alpha^*(q) </math> that verifies <math>
f'(\alpha^*(q))  = q
</math>
we have
<center><math>
\tau(q)= \alpha^*(q) q  -f(\alpha^*(q)) 
</math></center>


== Delocalized eigenstates ==


<Strong> Delocalized wave functions</Strong>  have a simple spectrum: For <math>\alpha=d</math>, we have  <math>f(\alpha=d)=d</math> and  <math>f(\alpha\ne d ) =-\infty</math>. Then <math>\alpha^*(q)=d </math> becomes <math>q</math> independent.
In a delocalized state the probability is spread uniformly over the system:
<Strong> Multifractal wave functions </Strong>  smooth this edge  dependence and display a smooth spectrum with a maximum at <math>\alpha_0</math> with <math>f(\alpha_0)=d</math>. At <math>q=1</math>, <math>f'(\alpha_1)=1</math> and <math>f(\alpha_1)=\alpha_1</math>.


Sometimes one writes:
<math display="block">
<center><math>
|\psi_n|^{2} \approx L^{-d}.
IPR(q)=L^{D_q(1-q)} \quad \tau_q=D_q(1-q) 
</math>
</math></center>
Here <math>D_q</math> is q-dependent multifractal dimension, smaller than <math>d</math> and larger than zero.


=Larkin model=
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">
\tau_q=d(q-1).
</math>
 
== Localized eigenstates ==
 
In a localized state the wavefunction is concentrated within a region of size <math>\xi_{\text{loc}}</math>.
 
Roughly
 
<math display="block">
|\psi_n|^2 \sim \xi_{\text{loc}}^{-d}
</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.
</math>
 
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
 
<math display="block">
\tau_q.
</math>
 
The corresponding multifractal dimensions are defined as
 
<math display="block">
D_q=\frac{\tau_q}{q-1}.
</math>
 
It is useful to contrast fractal and multifractal scaling.
 
* In a '''fractal object''', a single exponent describes the scaling of the measure. If the wavefunction amplitudes scale with a single exponent <math>D</math>, then
 
<math display="block">
\mathrm{IPR}(q)\sim L^{-D(q-1)} .
</math>
 
All moments are controlled by the same dimension <math>D</math>.
 
* In a '''multifractal object''', different regions of the system scale with different exponents. In this case there is no single fractal dimension: different moments probe different effective dimensions of the wavefunction.
 
== Multifractal spectrum ==
 
To describe this structure it is useful to introduce the '''multifractal spectrum''' <math>f(\alpha)</math>.
 
We assume that the wavefunction amplitudes scale as
 
<math display="block">
|\psi_n|^2 \sim L^{-\alpha}
</math>
 
on approximately
 
<math display="block">
L^{f(\alpha)}
</math>
 
sites.
 
The function <math>f(\alpha)</math> therefore describes the fractal dimension of the set of sites where the wavefunction has exponent <math>\alpha</math>.
 
Using this representation,
 
<math display="block">
\mathrm{IPR}(q)
=
\sum_n |\psi_n|^{2q}
\sim
\int d\alpha \,
L^{-\alpha q} L^{f(\alpha)} .
</math>
 
For large <math>L</math>, the integral is dominated by the saddle point, giving
 
<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>
 
which corresponds to the most typical amplitude of the wavefunction.
 
It is useful to contrast this situation with the case of a simple fractal.
 
In a fractal object, the measure has a single scaling exponent. In our notation this means that the wavefunction amplitudes scale with a single value <math>\alpha_0</math> on a set of fractal dimension <math>D</math>:
 
<math display="block">
|\psi_n|^2 \sim L^{-\alpha_0}
</math>
 
on approximately
 
<math display="block">
L^{D}
</math>
 
sites, and are essentially negligible elsewhere.
 
In terms of the multifractal spectrum this corresponds to a trivial spectrum
 
<math display="block">
f(\alpha_0)=D,
\qquad
f(\alpha\neq\alpha_0)=-\infty.
</math>
 
By contrast, in a multifractal state the amplitudes are broadly distributed and many values of <math>\alpha</math> contribute. The function <math>f(\alpha)</math> then becomes a smooth curve describing the distribution of scaling exponents of the wavefunction.
 
<!-- = Larkin model =
 
In your homework you solved a toy model for the interface
 
<math display="block">
\partial_t h(r,t) = \nabla^2 h(r,t) + F(r).
</math>


In your homewoork you solved a toy model for the interface:
<center><math>
\partial_t h(r,t) =  \nabla^2 h(r,t)  + F(r)
</math></center>
For simplicity, we assume Gaussian disorder
For simplicity, we assume Gaussian disorder
<math>\overline{F(r)}=0</math>, <math> \overline{F(r)F(r')}=\sigma^2 \delta^d(r-r') </math>.
 
<math display="block">
\overline{F(r)}=0,
\qquad
\overline{F(r)F(r')}=\sigma^2 \delta^d(r-r').
</math>


You proved that:
You proved that:
* the roughness exponent of this model is  <math>\zeta_L=\frac{4-d}{2}</math> below dimension 4
* The force per unit length acting on the center of the interface is <math> f= \sigma/\sqrt{L^d}</math>
* at long times the  interface shape is
<center><math>
\overline{h(q)h(-q)} =  \frac{\sigma^2}{q^{d+2\zeta_L}}
</math></center>


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, <math> r_f</math> the length of correlation of the disorder <Strong> along the h direction </Strong>. This defines a Larkin length. Indeed from
* the roughness exponent of this model is <math>\zeta_L=\frac{4-d}{2}</math> below dimension 4
<center><math>
* the force per unit length acting on the center of the interface is <math>f= \sigma/\sqrt{L^d}</math>
\overline{(h(r)-h(0))^2}= \int _d^dq (\overline{h(q)h(-q)}(1-\cos(q r) \sim \sigma^2 r^{2 \zeta_L}
* at long times the interface shape is
</math></center>
 
You get
<math display="block">
<center><math>
\overline{h(q)h(-q)} = \frac{\sigma^2}{q^{d+2\zeta_L}}.
\overline{(h(\ell_L)-h(0))^2}= r_f^2 \quad \ell_L=\left(\frac{r_f}{\sigma} \right)^{1/\zeta_L}
</math>
</math></center>
 
Above this scale, roguhness change and pinning starts with a crtical force
In the real depinning model the disorder is, however, a non-linear function of <math>h</math>. The idea of Larkin is that this linearization is valid only up to a scale <math>r_f</math>, the correlation length of the disorder along the <math>h</math> direction. This defines the '''Larkin length'''.
<center><math>
 
f_c= \frac{\sigma}{\ell_L^{d/(2 \zeta_L)}}
Indeed, starting from
</math></center>
 
In <math>d=1</math> we have <math>\ell_L=\left(\frac{r_f}{\sigma} \right)^{2/3}</math>
<math display="block">
\overline{(h(r)-h(0))^2}
= \int d^d q \,\overline{h(q)h(-q)}\,(1-\cos(q r))
\sim \sigma^2 r^{2\zeta_L},
</math>
 
we obtain
 
<math display="block">
\overline{(h(\ell_L)-h(0))^2}= r_f^2,
\qquad
\ell_L=\left(\frac{r_f}{\sigma}\right)^{1/\zeta_L}.
</math>
 
Above this scale the roughness changes and pinning sets in with a critical force
 
<math display="block">
f_c= \frac{\sigma}{\ell_L^{d/(2 \zeta_L)}}.
</math>
 
In <math>d=1</math> we therefore obtain
 
<math display="block">
\ell_L=\left(\frac{r_f}{\sigma}\right)^{2/3}.
</math>
-->

Latest revision as of 22:07, 23 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(q1).

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.

It is useful to contrast fractal and multifractal scaling.

  • In a fractal object, a single exponent describes the scaling of the measure. If the wavefunction amplitudes scale with a single exponent D, then

IPR(q)LD(q1).

All moments are controlled by the same dimension D.

  • In a multifractal object, different regions of the system scale with different exponents. In this case there is no single fractal dimension: different moments probe different effective dimensions of the wavefunction.

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(α) therefore 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,

which corresponds to the most typical amplitude of the wavefunction.

It is useful to contrast this situation with the case of a simple fractal.

In a fractal object, the measure has a single scaling exponent. In our notation this means that the wavefunction amplitudes scale with a single value α0 on a set of fractal dimension D:

|ψn|2Lα0

on approximately

LD

sites, and are essentially negligible elsewhere.

In terms of the multifractal spectrum this corresponds to a trivial spectrum

f(α0)=D,f(αα0)=.

By contrast, in a multifractal state the amplitudes are broadly distributed and many values of α contribute. The function f(α) then becomes a smooth curve describing the distribution of scaling exponents of the wavefunction.