L-9: Difference between revisions

From Disordered Systems Wiki
Jump to navigation Jump to search
Line 1: Line 1:
=Eigenstates=
=Eigenstates=
In absence of disorder the eigenstates are delocalized plane waves.
Without disorder, the eigenstates are delocalized plane waves.


In presence of disorder, three situations can occur and to distinguish them it is useful to introduce the inverse participation ratio, IPR
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.
 
To distinguish these regimes, it is useful to introduce the inverse participation ratio (IPR).
<center><math>
<center><math>
IPR(q)=\sum_n |\psi_n|^{2 q} \sim L^{-\tau_q}
IPR(q)=\sum_n |\psi_n|^{2 q} \sim L^{-\tau_q}
</math></center>
</math></center>
The normalization imposes <math>\tau_1 =0 </math>. For  <math>q=0</math>, <math>|\psi_n|^{2 q} =1 </math>, hence,  <math>\tau_0 =-d </math>.


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


* <Strong> Localized eigenstates</Strong> 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  
==Localized eigenstates==
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>
<center><math>
IPR(q)= \text{const},  \quad \tau_q=0
IPR(q)= \text{const},  \quad \tau_q=0
</math></center>
</math></center>


* <Strong> Multifractal eigenstates.</Strong>  At the transition(  the mobility edge) an anomalous scaling is observed:
<center><math>
IPR(q)=L^{D_q(1-q)}  \quad \tau_q=D_q(1-q) 
</math></center>
Here <math>D_q</math> is q-dependent multifractal dimension, smaller than <math>d</math> and larger than zero.


==Multifractal eigenstates==


 
The exponent <math>\tau_q</math> is called <Strong> multifractal exponent </Strong>. It is a non decreasing  function with q with some special points:
==Multifractality==
* <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.
 
The exponent <math>\tau_q</math> is called <Strong> multifractal exponent </Strong>:
 
* <math>\tau_1 =0 </math> imposed by normalization.
* <math>\tau_1 =0 </math> imposed by normalization.
* <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.
=== Delocalized eigenstates ===
In this case, <math>|\psi_n|^{2} \approx L^{-d} </math> for all the <math> L^{d} </math> sites. This gives
<center><math>
\tau_q^{\text{deloc}}=d(q-1)
</math></center>
* <Strong> Multifractal eigenstates.</Strong> 
This case correspond to more complex wave function for which
we expect


To have multifractal behaviour  we expect
<center><math>
<center><math>
|\psi_n|^{2} \approx L^{-\alpha} \quad  \text{for}\; L^{f(\alpha)} \; \text{sites}   
|\psi_n|^{2} \approx L^{-\alpha} \quad  \text{for}\; L^{f(\alpha)} \; \text{sites}   
</math></center>
</math></center>
The exponent <math>\alpha </math> is positive and <math>f(\alpha)</math> is called <Strong> multifractal spectrum </Strong>. It is a convex function and its maximum is the fractal dimension of the object, in our case d. We can determine the relation between multifractal spectrum and exponent  
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>
<center><math>
IPR(q)=\sum_n |\psi_n|^{2 q}\sim \int d \alpha L^{-\alpha q} L^{f(\alpha)}  
IPR(q)=\sum_n |\psi_n|^{2 q}\sim \int d \alpha L^{-\alpha q} L^{f(\alpha)}  
Line 59: Line 42:
This means that for <math>\alpha^*(q) </math> that verifies <math>
This means that for <math>\alpha^*(q) </math> that verifies <math>
f'(\alpha^*(q))  = q  
f'(\alpha^*(q))  = q  
</math> we have
</math>  
 
we have
<center><math>
<center><math>
\tau(q)= \alpha^*(q) q  -f(\alpha^*(q))   
\tau(q)= \alpha^*(q) q  -f(\alpha^*(q))   
Line 66: Line 49:




<Strong> A metal</Strong>  has a simple spectrum. Indeed, all sites have  <math>\alpha=d</math>, hence  <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.
<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.
<Strong> A multifractal </Strong> has 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>.
<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:
<center><math>
IPR(q)=L^{D_q(1-q)}  \quad \tau_q=D_q(1-q) 
</math></center>
Here <math>D_q</math> is q-dependent multifractal dimension, smaller than <math>d</math> and larger than zero.


=Larkin model=
=Larkin model=

Revision as of 15:05, 23 March 2025

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; 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.

To distinguish these regimes, it is useful to introduce the inverse participation ratio (IPR).

IPR(q)=n|ψn|2qLτq

Delocalized eigenstates

In this case, |ψn|2Ld. Hence, we expect

IPR(q)=Ld(1q)τq=d(1q)

Localized eigenstates

In this case, |ψn|21/ξloc1/d for ξlocd sites and almost zero elsewhere. Hence, we expect

IPR(q)=const,τq=0


Multifractal eigenstates

The exponent τq is called multifractal exponent . It is a non decreasing function with q with some special points:

  • τ0=d because the wave fuction is defined on all sites, in general τ0 is the fractal dimension of the object we are considering. It is simply a geometrical property.
  • τ1=0 imposed by normalization.

To have multifractal behaviour we expect

|ψn|2LαforLf(α)sites

The exponent α is positive and f(α) is called multifractal spectrum . Its maximum is the fractal dimension of the object, in our case d. We can determine the relation between multifractal spectrum f(α) and exponent τq

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

for large L

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

This means that for α*(q) that verifies f(α*(q))=q we have

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


Delocalized wave functions have a simple spectrum: For α=d, we have f(α=d)=d and f(αd)=. Then α*(q)=d becomes q independent. Multifractal wave functions smooth this edge dependence and display a smooth spectrum with a maximum at α0 with f(α0)=d. At q=1, f(α1)=1 and f(α1)=α1.

Sometimes one writes:

IPR(q)=LDq(1q)τq=Dq(1q)

Here Dq is q-dependent multifractal dimension, smaller than d and larger than zero.

Larkin model

In your homewoork 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, rf the length of correlation of the disorder along the h direction . This defines a Larkin length. Indeed from

(h(r)h(0))2=ddq(h(q)h(q)(1cos(qr)σ2r2ζL

You get

(h(L)h(0))2=rf2L=(rfσ)1/ζL

Above this scale, roguhness change and pinning starts with a crtical force

fc=σLd/(2ζL)

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