LBan-1: Difference between revisions
| Line 160: | Line 160: | ||
In the Gaussian case, we start from the third relation introduced earlier and expand <math>A(E)</math> around <math>a_M</math>: | In the Gaussian case, we start from the third relation introduced earlier and expand <math>A(E)</math> around <math>a_M</math>: | ||
<center> <math>Q_M(E) \sim \exp\left(-M P(E)\right) \sim \exp\left(-M e^{A(a_M) +A'(a_M)E-a_M)}\right) </math> </center> | <center> <math>Q_M(E) \sim \exp\left(-M P(E)\right) \sim \exp\left(-M e^{A(a_M) +A'(a_M)\cdot (E-a_M)}\right) </math> </center> | ||
By setting | By setting | ||
<center> <math>a_M = E_{\min}^{\text{typ}}=-\sigma \sqrt{2 \log M} + \ldots \quad \text{and} \quad b_M = \frac{1}{A'(a_M)}= \frac{ \sigma}{\sqrt{2 \log M}}</math> </center> | <center> <math>a_M = E_{\min}^{\text{typ}}=-\sigma \sqrt{2 \log M} + \ldots \quad \text{and} \quad b_M = \frac{1}{A'(a_M)}= \frac{ \sigma}{\sqrt{2 \log M}}</math> </center> | ||
Revision as of 16:36, 6 August 2025
Overview
This lesson is structured in three parts:
- Self-averaging and disorder in statistical systems
Disordered systems are characterized by a random energy landscape, however, in the thermodynamic limit, physical observables become deterministic. This property, known as self-averaging, does not always hold for the partition function which is the quantity that we can compute. When it holds the annealed average and the quenched average coincides otherwiese we have
- The Random Energy Model
We study the Random Energy Model (REM) introduced by Bernard Derrida. In this model at each configuration is assigned an independent energy drawn from a Gaussian distribution of extensive variance. The model exhibits a freezing transition at a critical temperature, below which the free energy becomes dominated by the lowest energy states.
- Extreme value statistics and saddle-point analysis
The results obtained from a saddle-point approximation can be recovered using the tools of extreme value statistics.
Part I
Random energy landascape
In a system with degrees of freedom, the number of configurations grows exponentially with . For simplicity, consider Ising spins that take two values, , located on a lattice of size in dimensions. In this 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 N = L^d} and the number of configurations 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 M = 2^N = e^{N \log 2}} .
In the presence of disorder, the energy associated with a given configuration becomes a random quantity. For instance, in the Edwards-Anderson model:
where the sum runs over nearest neighbors 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 \langle i, j \rangle} , and the 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_{ij}} are independent and identically distributed (i.i.d.) Gaussian random variables with zero mean and unit variance.
The energy of a given configuration is a random quantity because each system corresponds to a different realization of the disorder. In an experiment, this means that each of us has a different physical sample; in a numerical simulation, it means that each of us has generated a different set of 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_{ij}} .
To illustrate this, consider a single configuration, for example the one where all spins are up. The energy of this configuration is given by the sum of all the couplings between neighboring spins:
Since the the couplings are random, the energy associated with this particular configuration is itself a Gaussian random variable, with zero mean and a variance proportional to the number of terms in the sum — that is, of order 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} . The same reasoning applies to each of 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. So, in a disordered system, the entire energy landscape is random and sample-dependent.
Self-averaging observables
A crucial question is whether the macroscopic properties measured on a given sample are themselves random or not. Our everyday experience suggests that they are not: materials like glass, ceramics, or bronze have well-defined, reproducible physical properties that can be reliably controlled for industrial applications.
From a more mathematical point of view, it means that the free energy 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 F_N(\beta)=N f_N(\beta)} and its derivatives (magnetization, specific heat, susceptibility, etc.), in 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 N \to \infty } , these random quantities concentrates around a well defined value. These observables are called self-averaging. This means that,
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_{N \to \infty} f_N (\beta)= \lim_{N \to \infty} f_N^{\text{typ}}(\beta) =\lim_{N \to \infty} \overline{f_N(\beta)} =f_\infty(\beta) }
Hence 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 f_N(\beta) } becomes effectively deterministic and its sample-to sample fluctuations vanish in relative terms:
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_{N \to \infty} \frac{\overline{f_N^2(\beta)}}{\overline{f_N(\beta)}^2}=1. }
The partition function
The partition function
is itself a random variable in disordered systems. Analytical methods can capture the statistical properties of this variable. We can define to average over the disorder realizations:
- The annealed average corresponds to the calculation of the moments of the partition function. The annealed free energy is
- the quenched average corresponds to the average of the logarithm of the partition function, which is self-averaging for sure.
Do these two averages coincide?
If the partition function is self-averaging in the thermodynamic limit, then
As a consequence, the annealed and the quenched averages coincide.
If the partition function is not self-averaging, only typical partition function concentrates, but extremely rare configurations contribute disproportionately to its moments:
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_{N \to \infty} Z_N^{\text{typ}}(\beta) = e^{-\beta N f_\infty(\beta)} < \lim_{N \to \infty} \overline{Z_N(\beta)} = e^{-\beta N f^{\text{ann.}}(\beta)} }
There are then two main strategies to determine the deterministic value of the observable :
- Compute directly the quenched average using methods such as the replica trick and the Parisi solution.
- Determine the typical value 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_N^{\text{typ}}(\beta) } and evaluate
Part II
Random Energy Model
The Random energy model (REM) neglects the correlations between 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. The energy associated to each configuration is an independent Gaussian variable with zero mean and variance 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} . The simplest solution of the model is with the microcanonical ensemble.
Microcanonical calculation
Step 1: Number of states .
Let the number of states of energy in the interval (E,E+dE). It is a random number and we use the representation
with 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 \chi_\alpha(E)=1} if and 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 \chi_\alpha(E)=0} otherwise. We can cumpute its average
Here 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 \epsilon =E/N } is the energy density and the annealed entropy density in the thermodynamic limit is
Step 2: Self-averaging.
Let compute now the second moment
We can then check the self averaging condition:
A critical energy density 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 \epsilon^* = \sqrt{2 \ln 2}} separates a self-averaging regime for 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 |\epsilon| < \epsilon^*} and a non self-averaging regime where for . In the first regime, 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{\mathcal{N}_N(E)}} is exponentially large and its value is determinstic (average, typical, median are the same). In the secon regime, 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{\mathcal{N}_N(E)}} is exponentially small but nonzero. The typical value instead is exactly zero, : for most disorder realizations, there are no configurations with energy below 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 - \epsilon^* N} and only a vanishingly small fraction of rare samples gives a positive contribution to the average. As a result, the quenched average on the entropy density is:
Back to canonical ensemble: the freezing transition
The annealed partition function is the average of the partition function over the disorder:
Using the saddle point for large N we find and thus
The quenched partition function is obtained replacing the mean with the typical value:
Using the saddle point for large N we find a critical 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_c = \sqrt{2 \ln 2}} separating two phases:
- For 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 < \beta_c } , 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 \epsilon_\min =-\beta} and the annealed calculation works
- For 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 > \beta_c } , 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 \epsilon_\min =-\beta_c} and the free energy freezes to a temperature independent value. As a result, the quenched average on the free energy density is:
Part III
Detour: Extreme Value Statistics
Consider the REM spectrum 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 M} energies 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_1, \dots, E_M} drawn from a 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 introduce the cumulative probability of finding an energy smaller than E
We also define:
The statistical properties 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 E_{\min}} are derived using three key relations:
- First relation:
This relation is exact but depends on M and the precise form 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 p(E)} .
- Second relation:
This is an estimation of the typical value of the minimum. It is a crucial relation that will be used frequently.
- Third relation
This is an approximation valid for large M and around the typical value of the minimum energy. It allows to extract the universal scaling.
Back to REM
Let us analyze in detail the case of a Gaussian distribution with zero mean and variance 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^2} . Using integration by parts, we can write :
The asymptotic expansion for is :
In this case it is convenient to introduce the function 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(E)} defined 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 P(E) = \exp(A(E))} . hence we have:
From the second relation we impose 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_{\text{gauss}}(E_{\min}^{\text{typ}})=- \ln M} . For large 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} we get:
We look for a scaling form of the random variable 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} = a_M + b_M z}
such that 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}
, 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 b_M}
are deterministic and M dependent while 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}
is random and M independent.
In the Gaussian case, we start from the third relation introduced earlier and expand 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(E)}
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}
:
By setting