T-5: Difference between revisions
No edit summary |
|||
Line 1: | Line 1: | ||
<strong>Goal: </strong> | <strong>Goal: </strong> | ||
So far we have discussed the equilibrium properties of disordered systems, that are encoded in their partition function | So far we have discussed the equilibrium properties of disordered systems, that are encoded in their partition function/free energy. When a system (following Langevin, Monte Carlo dynamics) equilibrates at sufficiently large times, its long-time properties are captured by these equilibrium calculations. In glassy systems the equilibration timescales are extremely large: for very large timescales the system does not visit equilibrium configurations, but rather metastable states. In this set of problems, we characterize the energy landscape of the spherical <math>p</math>-spin by studying its metastable states (local minima). | ||
<br> | <br> | ||
<strong>Techniques: </strong> differential geometry (just a tiny bit!), random matrix theory. | <strong>Techniques: </strong> differential geometry (just a tiny bit!), random matrix theory. | ||
Line 25: | Line 25: | ||
<br> | <br> | ||
<li> '''Stationary points and complexity.''' To guess where gradient descent dynamics (or | <li> '''Stationary points and complexity.''' To guess where gradient descent dynamics (or <ins> Langevin dynamics </ins>) are expected to converge, it is useful to understand the distribution of the stationary points, i.e. the number <math> \mathcal{N}(\epsilon)</math> of such configuration having a given energy density <math> \epsilon </math>. In fully-connected models, this quantity has an exponential scaling, <math> \mathcal{N}(\epsilon) \sim \text{exp}\left(N \Sigma(\epsilon) \right)</math>, where <math> \Sigma(\epsilon)</math> is the landscape’s <ins>complexity</ins>. <sup>[[#Notes|[*] ]]</sup>. Stationary points can be stable (local minima), or unstable (saddles or local maxima): their stability is encoded in the spectrum of the <ins> Hessian matrix </ins> <math>\nabla_{\perp}^2 E(\vec{\sigma})</math>: when all the eigenvalues of the Hessian are positive, the point is a local minimum (and a saddle otherwise). | ||
</li> | </li> | ||
</ul> | </ul> | ||
Line 44: | Line 44: | ||
=== Problem 5.1: the Kac-Rice formula and the complexity === | === Problem 5.1: the Kac-Rice formula and the complexity === | ||
First, a few notions of geometry: we define with <math> \hat \Pi(\vec{ | First, a few notions of geometry: we define with <math> \hat \Pi(\vec{\sigma}) </math> the projector on the tangent plane to the sphere at <math> \vec{\sigma}</math>: this is the plane orthogonal to the vector <math> \vec{\sigma}</math>. The gradient <math> \nabla_\perp E(\vec{\sigma}) </math> is a <math> (N-1)</math>-dimensional vector that is obtained projecting the gradient <math> [\nabla E(\vec{\sigma})]_i=\partial E/\partial \sigma_i </math> on the tangent plane, <math> \nabla_\perp E(\vec{\sigma})=\hat \Pi \nabla E(\vec{\sigma})</math>. The Hessian <math> \nabla^2_\perp E(\vec{\sigma}) </math> is a <math> (N-1) \times (N-1)</math>-dimensional matrix that is obtained from the Hessian <math> [\nabla^2 E(\vec{\sigma})]_{ij}=\partial^2 E/\partial \sigma_i \partial \sigma_j </math> as <math> \nabla^2_\perp E(\vec{\sigma})= \hat \Pi(\vec{\sigma}) \, \nabla^2 E(\vec{\sigma}) \, \hat\Pi(\vec{\sigma}) - N^{-1}\nabla E(\vec{\sigma}) \cdot \vec{\sigma} \mathbb{I}</math> where <math> \mathbb{I}</math> is the identity matrix. <br> | ||
Revision as of 18:51, 15 February 2024
Goal:
So far we have discussed the equilibrium properties of disordered systems, that are encoded in their partition function/free energy. When a system (following Langevin, Monte Carlo dynamics) equilibrates at sufficiently large times, its long-time properties are captured by these equilibrium calculations. In glassy systems the equilibration timescales are extremely large: for very large timescales the system does not visit equilibrium configurations, but rather metastable states. In this set of problems, we characterize the energy landscape of the spherical -spin by studying its metastable states (local minima).
Techniques: differential geometry (just a tiny bit!), random matrix theory.
Dynamics, optimization, trapping local minima
- Rugged landscapes. Consider the spherical -spin model for concreteness: is an energy landscape . It is a random function on configuration space (in our case, the surface of the sphere in dimension 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 } ). This landscape has its global minimum(a) at the ground state configuration(s): the energy density of the ground state(s) can be obtained studying the partition 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 Z } in the limit . Besides the ground state(s), the energy landscape can have other local minima; fully-connected models of glasses are characterized by the fact that there are plenty of these local minima: the energy landscape is rugged, see the sketch.
- Optimization by gradient descent. Suppose that we are interested in finding the configurations of minimal energy, starting from an arbitrary 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 \vec{\sigma}_0}
: we can implement a dynamics in which we progressively update the configuration moving towards lower and lower values of the energy, hoping to eventually converge to the ground state(s). The simplest dynamics of this sort is gradient descent,
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 \frac{d \vec{\sigma}(t)}{dt}=- \nabla_{\perp} E(\vec{\sigma}) } where is the gradient of the landscape restricted to the sphere. The dynamics stops when it reaches a stationary point , a configuration 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 \nabla_\perp E(\vec{\sigma})=0} . If the landscape has a convex structure, this will be the ground state; if the energy landscape is very non-convex like in glasses, the end point of this algorithm will be a local minimum at energies much higher than the ground state (see sketch).
- Stationary points and complexity. To guess where gradient descent dynamics (or Langevin dynamics ) are expected to converge, it is useful to understand the distribution of the stationary points, i.e. the number 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 \mathcal{N}(\epsilon)} of such configuration having a given energy density . In fully-connected models, this quantity has an exponential scaling, 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 \mathcal{N}(\epsilon) \sim \text{exp}\left(N \Sigma(\epsilon) \right)} , 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 \Sigma(\epsilon)} is the landscape’s complexity. [*] . Stationary points can be stable (local minima), or unstable (saddles or local maxima): their stability is encoded in the spectrum of the Hessian matrix : when all the eigenvalues of the Hessian are positive, the point is a local minimum (and a saddle otherwise).
- [*] - This quantity looks similar to the entropy 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 S(\epsilon) } we computed for the REM in Problem 1. However, while the entropy counts all configurations at a given energy density, the complexity 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(\epsilon) } accounts only for the stationary points.
Problems
In these problems, we discuss the computation of the annealed complexity of the spherical -spin model, which is defined by
Problem 5.1: the Kac-Rice formula and the complexity
First, a few notions of geometry: we define 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 \hat \Pi(\vec{\sigma}) }
the projector on the tangent plane to the sphere at : this is the plane orthogonal to the vector 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 \vec{\sigma}}
. The gradient 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 \nabla_\perp E(\vec{\sigma}) }
is a -dimensional vector that is obtained projecting the gradient 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 [\nabla E(\vec{\sigma})]_i=\partial E/\partial \sigma_i }
on the tangent plane, 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 \nabla_\perp E(\vec{\sigma})=\hat \Pi \nabla E(\vec{\sigma})}
. The Hessian is a 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-1) \times (N-1)}
-dimensional matrix that is obtained from the Hessian 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 [\nabla^2 E(\vec{\sigma})]_{ij}=\partial^2 E/\partial \sigma_i \partial \sigma_j }
as 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 \mathbb{I}}
is the identity matrix.
- The Kac-Rice formula. Consider first a random function of one 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 f(x)}
defined on an interval , and let 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 \mathcal{N}}
be the number of points 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 x }
such that . Justify why
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}}= \int_a^b dx \,p_0(x) , \quad \quad p_0(x)=\overline{\delta(f(x)) |f'(x)|} }
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_0(x) } is the probability density that is a zero of the function. In particular, why is the derivative of the function appearing in this formula? Consider now the number of stationary points 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 \mathcal{N}(\epsilon)} 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 p} -spin energy landscape, which satisfy . Justify why the generalization of the formula above gives
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}(\epsilon)}= \int_{\mathcal{S}_N} d \vec{S} \,p_{\epsilon}(\vec{S}) , \quad \quad p_{\epsilon}(\vec{S})=\overline{|\text{det} \nabla_\perp^2 E (\vec{S})|\,\, \delta(\nabla_\perp E(\vec{S})=0) \, \,\delta(E(\vec{S})- N \epsilon)} }
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_{\epsilon}(\vec{S})} is the probability density that is a stationary point of 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 } , 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 \nabla_\perp^2 E (\vec{S}) } is the Hessian matrix of the function restricted to the sphere.
- Statistical rotational invariance. Recall the expression of the correlations of the energy landscape 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 p}
-spin computed in Problem 2.1: in which sense the correlation function is rotationally invariant? Justify why rotational invariance implies 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 \overline{\mathcal{N}(\epsilon)}= (2 \pi e)^{\frac{N}{2}} \, p_{\epsilon}(\vec{1}) }
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 \vec{1}=(1,1,1, \cdots, 1) } is one fixed vector belonging to the surface of the sphere. Where does the prefactor arise from?
- Gaussianity and correlations. Determine the distribution of the quantity 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 (\vec{1})}
. Show that the components of the vector 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 \nabla E (\vec{1})}
are Gaussian random variables with zero mean and covariances
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{(\nabla E)_i \, (\nabla E)_j}= \frac{p}{2} \, \delta_{ij}+O\left(\frac{1}{N} \right) }
The quantity 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 \nabla_\perp E (\vec{1})} can be shown to be uncorrelated 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 E (\vec{1}), \nabla^2_\perp E (\vec{1})} . The entries 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 (N-1)\times (N-1) } matrix 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 \nabla_\perp^2 E (\vec{\sigma}) } are also Gaussian variables. Computing their correlation, one finds that the matrix conditioned to the fact 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 E(\vec 1)=N \epsilon } can be written 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 [\nabla_\perp^2 E(\vec{1})]_{\alpha \beta}= \left[\hat \Pi \nabla^2 E (\vec{1}) \hat \Pi - N^{-1}p \, E (\vec{\sigma}) \mathbb{I} \right]_{\alpha \beta}= M_{\alpha \beta}- p \epsilon\, \delta_{\alpha \beta}, }
where the matrix 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 } has random entries with zero average and correlations
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{{M}_{\alpha \beta} \, {M}_{\gamma \delta}}= \frac{p (p-1)}{2 N} \left( \delta_{\alpha \gamma} \delta_{\beta \delta}+ \delta_{\alpha \delta} \delta_{\beta \gamma}\right) }
Combining everything, show that this implies
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}(\epsilon)}= (2 \pi e)^{\frac{N}{2}} \,\frac{1}{(\pi \, p)^{\frac{N-1}{2}}}\; \sqrt{\frac{N}{\pi}} e^{-N \epsilon^2}\;\overline{|\text{det} \left( M- p \epsilon \mathbb{I} \right)|} }
Problem 5.2: the Hessian and random matrix theory
To get the complexity, it remains to compute the expectation value of the determinant of the Hessian matrix: this is the goal of this problem. We will do this exploiting results from random matrix theory.
- Gaussian Random matrices. Show that the matrix 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 } is a GOE matrix, i.e. a matrix taken from the Gaussian Orthogonal Ensemble, meaning that it is a symmetric matrix with 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(M)= Z_N^{-1}e^{-\frac{N}{4 \sigma^2} \text{Tr} M^2}. } What is the value 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 \sigma^2 } ?
- Eigenvalue density and concentration. Let 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 \lambda_\alpha }
be the eigenvalues of the matrix 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 }
. Show that the following identity holds:
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{|\text{det} \left(M - p \epsilon \mathbb{I} \right)|}= \overline{\text{exp} \left[(N-1) \left( \int d \lambda \, \rho_N(\lambda) \, \log |\lambda - p \epsilon|\right) \right]}, \quad \quad \rho_{N}(\lambda)= \frac{1}{N-1} \sum_{\alpha=1}^{N-1} \delta (\lambda- \lambda_\alpha) }
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 \rho_{N}(\lambda)} is the empirical eigenvalue density. It can be shown that if 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 } is a GOE matrix, the distribution of the empirical density has a large deviation form (recall TD1) with speed 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^2 } , meaning 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 P_N[\rho] = e^{-N^2 \, g[\rho]} } where now 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 g[\cdot] } is a functional (a function of a function). Using a saddle point argument, show that this implies
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{\text{exp} \left[(N-1) \left( \int d \lambda \, \rho_N(\lambda) \, \log |\lambda - p \epsilon|\right) \right]}=\text{exp} \left[N \left( \int d \lambda \, \rho_{\text{ty}}(\lambda+p \epsilon) \, \log |\lambda|\right)+ o(N) \right] }
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 \rho_{\text{ty}}(\lambda) } is the typical value of the eigenvalue density, which satisfies 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 g[\rho_{\text{ty}}]=0 } .
- The semicircle, the threshold and the ground state. The eigenvalue density of GOE matrices is self-averaging, and it equals 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 \lim_{N \to \infty}\rho_N (\lambda)=\lim_{N \to \infty} \overline{\rho_N}(\lambda)= \rho_{\text{ty}}(\lambda)= \frac{1}{2 \pi \sigma^2}\sqrt{4 \sigma^2-\lambda^2 } }
- Check this numerically: generate matrices for various values 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 } , plot their empirical eigenvalue density and compare with the asymptotic curve. Is the convergence faster in the bulk, or in the edges of the eigenvalue density, where it vanishes?
- Sketch 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 \rho_{\text{ty}}(\lambda+p \epsilon) }
for different values 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 \epsilon }
; recalling that the Hessian encodes for the stability of the stationary points, show that there is a transition in the stability of the stationary points at a critical value of the 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_{\text{th}}= -\sqrt{\frac{2(p-1)}{p}} }
When are the critical point stable local minima? When are they saddles? Why the stationary points 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 \epsilon= \epsilon_{\text{th}}} are called marginally stable ?
- Combining all the results, show that the annealed complexity 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 \Sigma_{\text{a}}(\epsilon)= \frac{1}{2}\log [4 e (p-1)]- \epsilon^2+ I_p(\epsilon), \quad \quad I_p(\epsilon)= \frac{2}{\pi}\int d x \sqrt{1-\left(x- \frac{\epsilon}{ \epsilon_{\text{th}}}\right)^2}\, \log |x| , \quad \quad \epsilon_{\text{th}}= -\sqrt{\frac{2(p-1)}{p}}. } The integral 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 I_p(\epsilon)} can be computed explicitly, and one finds:
Plot the annealed complexity, and determine numerically where it vanishes: why is this a lower bound or the ground state energy density?
Check out: key concepts
Gradient descent, rugged landscapes, metastable states, Hessian matrices, random matrix theory, landscape’s complexity.