L-1: Difference between revisions

From Disordered Systems Wiki
Jump to navigation Jump to search
No edit summary
Line 43: Line 43:


The solution of the SK is difficult. To make progress we first study the radnom energy model (REM) introduced by B. Derrida.  
The solution of the SK is difficult. To make progress we first study the radnom energy model (REM) introduced by B. Derrida.  
 
This model neglects the correlations between the <math> M=2^N </math> configurations and assumes the <math> E_{\alpha} </math> as iid variables.
=== Derivation of the model===
The REM neglects the correlations between the <math> 2^N </math> configurations and assumes the <math> E_{\alpha} </math> as iid variables.
* Show that the energy distribution is  
* Show that the energy distribution is  
<center><math> p(E_\alpha) =\frac{1}{\sqrt{2 \pi \sigma^2}}e^{-\frac{E_{\alpha}^2}{2 \sigma^2}}</math></center>
<center><math> p(E_\alpha) =\frac{1}{\sqrt{2 \pi \sigma^2}}e^{-\frac{E_{\alpha}^2}{2 \sigma^2}}</math></center>
and determine <math>\sigma^2</math>
and determine <math>\sigma^2</math>


=== Solution of the Random Energy model  ===


We provide different solutions of the Random Energy Model (REM).  The first one focus on the statistics of the smallest energies among the ones associated to the <math>M=2^N</math> configurations.


Consider the  <math>M=2^N</math>  energies: <math>(E_1,...,E_M)</math>. They are i.i.d. variables, drawn from the Gaussian distribution <math>p(E)</math>.
We provide different solutions of the Random Energy Model (REM). The first one focus on the statistics of the smallest energies among the ones associated to the <math>M=2^N</math> configurations. For this, we need to become familiar with the main results of extreme value statistic of iid variables.
 
===  Extreme value statistics  ===
 
Consider the  <math>M=2^N</math>  energies: <math>(E_1,...,E_M)</math>. They are iid variables, drawn from the distribution <math>p(E)</math>.
It is useful to  use the following notations:
It is useful to  use the following notations:
* <math>P^<(E)=\int_{-\infty}^E dx p(x)  \sim \frac{\sigma}{\sqrt{2 \pi}|E|}e^{-\frac{E^2}{2 \sigma^2}} \; </math> for  <math>x \to -\infty</math>. It  represents the probability to find an energy smaller than ''E''.  
* <math>P^<(E)=\int_{-\infty}^E dx p(x)  \sim \frac{\sigma}{\sqrt{2 \pi}|E|}e^{-\frac{E^2}{2 \sigma^2}} \; </math> for  <math>x \to -\infty</math>. It  represents the probability to find an energy smaller than ''E''.  
* <math> P^>(E)=\int_E^{+\infty} dx p(x) = 1- P^<(E) </math>. It represents the probability to dfind an energy  larger than ''E''.
* <math> P^>(E)=\int_E^{+\infty} dx p(x) = 1- P^<(E) </math>. It represents the probability to dfind an energy  larger than ''E''.
At the end we will discuss the case where <math>p(E)</math> is Gaussian, but we can remain for general for this section.




====  Extreme value statistics for iid  ====
We denote
We denote
<center><math>E_{\min}=\min(E_1,...,E_M)</math></center>
<center><math>E_{\min}=\min(E_1,...,E_M)</math></center>
Line 75: Line 76:
Moreover, if we define <math>  P^<(\epsilon)=\exp(-A(\epsilon)) </math> we recover the famous Gumbel distribution:
Moreover, if we define <math>  P^<(\epsilon)=\exp(-A(\epsilon)) </math> we recover the famous Gumbel distribution:
<center><math>Q_M(\epsilon) \sim \exp\left(-e^{- A'(a_M) (\epsilon-a_M)}\right)  </math> </center>
<center><math>Q_M(\epsilon) \sim \exp\left(-e^{- A'(a_M) (\epsilon-a_M)}\right)  </math> </center>
From these results we conclude that:
* the minimum is typically located around <math> a_M  </math>
* the fluctuations of the minimum (i.e. its standard deviation) scale as <math>1/A'(a_M)</math>
===== Exercise L1-A: the Gaussian case =====
===== Exercise L1-A: the Gaussian case =====
Specify these results to the Guassian case and find
Specify these results to the Guassian case and find
Line 84: Line 89:


===Density of states above the minimum===
===Density of states above the minimum===
For a given disorder realization, we compute <math> d(x) </math>, the number of configurations above the minimum, but with an energy smaller than <math> E_{\min}+x</math>.
For a given disorder realization, we compute <math> d(x) </math>, the number of configurations above the minimum with an energy smaller than <math> E_{\min}+x</math>. The key relation for this quantity is:
<center><math>  \text{Prob}(d(x) = k) = M \binom{M-1}{k}\int  dE \; p(E) [P^>(E) -  P^>(E+x)  ]^{k} P^>(E+x)^{M - k - 1}
<center><math>  \text{Prob}(d(x) = k) = M \binom{M-1}{k}\int  dE \; p(E) [P^>(E) -  P^>(E+x)  ]^{k} P^>(E+x)^{M - k - 1}
     </math></center>
     </math></center>
Line 91: Line 96:
     \overline{d(x)} = M (M-1) \int  dE \; p(E) \left[P^>(E) -  P^>(E+x)  \right] P^>(E)^{M-2}   
     \overline{d(x)} = M (M-1) \int  dE \; p(E) \left[P^>(E) -  P^>(E+x)  \right] P^>(E)^{M-2}   
</math></center>
</math></center>
 
Now, in the integral <math> E </math> is the energy of the minimum, hence we can use
=== Number ===
 
 


==Bibliography==
==Bibliography==

Revision as of 11:11, 28 November 2023

Spin glass Transition

Experiments

Parlare dei campioni di rame dopati con il magnesio, marino o no: trovare due figure una di suscettivita e una di calore specifico, prova della transizione termodinamica.

Edwards Anderson model

We consider for simplicity the Ising version of this model.

Ising spins takes two values Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \sigma=\pm 1} and live on a lattice of Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle N } sitees . The enregy is writteen as a sum between the nearest neighbours <i,j>:

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle E= - \sum_{ <i, j> } J_{ij} \sigma_i \sigma_j }

Edwards and Anderson proposed to study this model for couplings that are i.i.d. random variables with zero mean. We set Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \pi(J)} the coupling distribution indicate the avergage over the couplings called disorder average, with an overline:

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \overline{J} \equiv \int d J \, J \, \pi(J)=0 }

It is crucial to assume Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \bar{ J}=0 } , otherwise the model displays ferro/antiferro order. We sill discuss two distributions:

  • Gaussian couplings:
  • Coin toss couplings, Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle J= \pm 1 } , selected with probability Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle 1/2 } .

Edwards Anderson order parameter

The SK model

Sherrington and Kirkpatrik considered the fully connected version of the model with Gaussian couplings:

At the inverse temperature Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \beta } , the partion function of the model is

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle Z= \sum_{\alpha=1}^{2^N} z_{\alpha}, \quad \text{with}\; z_{\alpha}= e^{-\beta E_\alpha} }

Here is the energy associated to the configuration Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \alpha } . This model presents a thermodynamic transition at Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \beta_c=?? } .

Random energy model

The solution of the SK is difficult. To make progress we first study the radnom energy model (REM) introduced by B. Derrida. This model neglects the correlations between the configurations and assumes the Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle E_{\alpha} } as iid variables.

  • Show that the energy distribution is
Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle p(E_\alpha) =\frac{1}{\sqrt{2 \pi \sigma^2}}e^{-\frac{E_{\alpha}^2}{2 \sigma^2}}}

and determine


We provide different solutions of the Random Energy Model (REM). The first one focus on the statistics of the smallest energies among the ones associated to the Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle M=2^N} configurations. For this, we need to become familiar with the main results of extreme value statistic of iid variables.

Extreme value statistics

Consider the Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle M=2^N} energies: . They are iid variables, drawn from the distribution Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle p(E)} . It is useful to use the following notations:

  • Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle P^<(E)=\int_{-\infty}^E dx p(x) \sim \frac{\sigma}{\sqrt{2 \pi}|E|}e^{-\frac{E^2}{2 \sigma^2}} \; } for . It represents the probability to find an energy smaller than E.
  • Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle P^>(E)=\int_E^{+\infty} dx p(x) = 1- P^<(E) } . It represents the probability to dfind an energy larger than E.

At the end we will discuss the case where Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle p(E)} is Gaussian, but we can remain for general for this section.


We denote

Our goal is to compute the cumulative distribution Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle Q_M(\epsilon)\equiv\text{Prob}(E_{\min}> \epsilon)} for large M and iid variables.

We need to understand two key relations:

  • The first relation is exact:
Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle Q_M(\epsilon) = \left(P^>(\epsilon)\right)^M }
  • The second relation identifies the typical value of the minimum, namely :
Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle P^<(a_M) = \frac1 M }

.

Close to Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle a_M } , we expect . Hence, from the limit Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \lim_{M\to \infty} (1-\frac{k}{M})^M =\exp(-k)} we re-write the first relation:

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle Q_M(\epsilon) \sim \exp\left(-M P^<(\epsilon)\right)}

Moreover, if we define we recover the famous Gumbel distribution:

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle Q_M(\epsilon) \sim \exp\left(-e^{- A'(a_M) (\epsilon-a_M)}\right) }

From these results we conclude that:

  • the minimum is typically located around Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle a_M }
  • the fluctuations of the minimum (i.e. its standard deviation) scale as Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle 1/A'(a_M)}
Exercise L1-A: the Gaussian case

Specify these results to the Guassian case and find

  • the typical value of the minimum

%

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle a_M = \sigma \sqrt{2 \log M}-\frac{1}{2}\sqrt{\log(\log M)} +O(1) }
  • The expression Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle A(\epsilon) =\frac{\epsilon^2}{2\sigma^2} -\frac{\sqrt{2 \pi}}{\sigma} \log|\epsilon|+\ldots }
  • The expression of the Gumbel distribution for the Gaussian case
Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle Q_M(\epsilon) \sim \exp\left(-e^{- \frac{\sqrt{2 \log M}}{\sigma} (\epsilon-a_M)}\right) }

Density of states above the minimum

For a given disorder realization, we compute Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle d(x) } , the number of configurations above the minimum with an energy smaller than Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle E_{\min}+x} . The key relation for this quantity is:

Taking the average Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \overline{d(x)} = \sum_k k \text{Prob}(d(x) = k) } , we derive

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \overline{d(x)} = M (M-1) \int dE \; p(E) \left[P^>(E) - P^>(E+x) \right] P^>(E)^{M-2} }

Now, in the integral is the energy of the minimum, hence we can use

Bibliography

Bibliography

  • Theory of spin glasses, S. F. Edwards and P. W. Anderson, J. Phys. F: Met. Phys. 5 965, 1975