Metastable states in spin glasses with short-ranged interactions

Free-energy barriers in the Sherrington-Kirkpatrick model

The free-energy landscape of the Sherrington-Kirkpatrick (SK) Ising spin glass is simple in the framework of the Thouless-Anderson-Palmer (TAP) equations as each solution (which are minima of the free energy) has associated with it a nearby index-one saddle point. The free-energy barrier to escape the minimum is just the difference between the saddle point free energy and that at its associated minimum. This difference is calculated for the states with free energies f>f_{c}. It is very small for these states, decreasing as 1/N^{2}, where N is the number of spins in the system. These states are not marginally stable. We argue that such small barriers are why numerical studies never find these states when N is large. Instead, the states that are found are those that have marginal stability. For them the barriers are at least of O(1). f_{c} is the free energy per spin below which the states develop broken replica-symmetry-like overlaps with each other. In the regime f<f_{c} we can only offer some possibilities based around scaling arguments. One of these suggest that the barriers might become as large as N^{1/3}. That might be consistent with recent numerical studies on the Viana-Bray model, which were at variance with the expectations of Cugliandolo and Kurchan for the SK model.

Realizable solutions of the Thouless-Anderson-Palmer equations

We show that the only solutions of the Thouless-Anderson-Palmer (TAP) equations for the Sherrington-Kirkpatrick model of Ising spin glasses which can be found by iteration are those whose free energy lies on the border between replica-symmetric and broken-replica-symmetric states, when the number of spins N is large. Convergence to this same borderline also happens in quenches from a high-temperature initial state to a locally stable state where each spin is parallel to its local field; both are examples of self-organized criticality. At this borderline the band of eigenvalues of the Hessian associated with a solution extends to zero, so the states reached have marginal stability. We have also investigated the factors which determine the free-energy difference between a stationary solution corresponding to a saddle point and its associated minimum, which is the barrier which has to be surmounted to escape from the vicinity of a TAP minimum or pure state.

Ising spin glasses in two dimensions: Universality and nonuniversality

Following numerous earlier studies, extensive simulations and analyses were made on the continuous interaction distribution Gaussian model and the discrete bimodal interaction distribution Ising spin glass (ISG) models in two dimensions [Lundow and Campbell, Phys. Rev. E 93, 022119 (2016)1539-375510.1103/PhysRevE.93.022119]. Here we further analyze the bimodal and Gaussian data together with data on two other continuous interaction distribution two-dimensional ISG models, the uniform and the Laplacian models, and three other discrete interaction distribution models, a diluted bimodal model, an "antidiluted" model, and a more exotic symmetric Poisson model. Comparisons between the three continuous distribution models show that not only do they share the same exponent η≡0 but that to within the present numerical precision they share the same critical exponent ν also, and so lie in a single universality class. On the other hand the critical exponents of the four discrete distribution models are not the same as those of the continuous distributions, and the present data strongly indicate that they differ from one discrete distribution model to another. This is evidence that discrete distribution ISG models in two dimensions have nonzero values of the critical exponent η and do not lie in a single universality class.

Metastable minima of the Heisenberg spin glass in a random magnetic field

We have studied zero-temperature metastable minima in classical m-vector component spin glasses in the presence of m-component random fields for two models, the Sherrington-Kirkpatrick (SK) model and the Viana-Bray (VB) model. For the SK model we have calculated analytically its complexity (the log of the number of minima) for both the annealed case where one averages the number of minima before taking the log and the quenched case where one averages the complexity itself, both for fields above and below the de Almeida-Thouless (AT) field, which is finite for m>2. We have done numerical quenches starting from a random initial state (infinite temperature state) by putting spins parallel to their local fields until there is no further decrease of the energy and found that in zero field it always produces minima that have zero overlap with each other. For the m=2 and m=3 cases in the SK model the final energy reached in the quench is very close to the energy E_{c} at which the overlap of the states would acquire replica symmetry-breaking features. These minima have marginal stability and will have long-range correlations between them. In the SK limit we have analytically studied the density of states ρ(λ) of the Hessian matrix in the annealed approximation. Despite the fact that in the presence of a random field there are no continuous symmetries, the spectrum extends down to zero with the usual sqrt[λ] form for the density of states for fields below the AT field. However, when the random field is larger than the AT field, there is a gap in the spectrum, which closes up as the AT field is approached. The VB model behaves differently and seems rather similar to studies of the three-dimensional Heisenberg spin glass in a random vector field.

Dealing with correlated choices: how a spin-glass model can help political parties select their policies

Starting from preferences on N proposed policies obtained via questionnaires from a sample of the electorate, an Ising spin-glass model in a field can be constructed from which a political party could find the subset of the proposed policies which would maximize its appeal, form a coherent choice in the eyes of the electorate, and have maximum overlap with the party's existing policies. We illustrate the application of the procedure by simulations of a spin glass in a random field on scale-free networks.

Counting metastable states of Ising spin glasses on arbitrary graphs

Using a field-theoretical representation of the Tanaka-Edwards integral, we develop a method to systematically compute the number Ns of one-spin stable states (local energy minima) of a glassy Ising system with nearest-neighbor interactions and random Gaussian couplings on an arbitrary graph. In particular, we use this method to determine Ns for K -regular random graphs and d -dimensional regular lattices for d=2,3 . The method works also for other graphs. Excellent accuracy of the results allows us to observe that the number of local energy minima depends mainly on local properties of the graph on which the spin glass is defined.

Conformational temperature characterizing the folding of a protein

The time sequences of the molecular dynamics simulation for the folding process of a protein is analyzed with the inherent structure landscape which focuses on the configurational dynamics of the system. Time-dependent energy and entropy for inherent structures are introduced, and from these quantities a conformational temperature is defined. The conformational temperature follows the time evolution of a slow relaxation process and reaches the bath temperature when the system is equilibrated. We show that the nonequilibrium system is described by two temperatures, one for fast vibration and the other for slow configurational relaxation, while the equilibrium system is described by one temperature. The proposed formalism is applicable widely for systems with many metastable states.

Modeling protein thermodynamics and fluctuations at the mesoscale

We use an extended Go model, in unfrustrated and frustrated variants, to study the energy landscape and the fluctuations of a model protein. The model exhibits two transitions, folding and dynamical transitions, when changing the temperature. The inherent structures corresponding to the minima of the landscape are analyzed and we show how their energy density can be obtained from simulations around the folding temperature. The scaling of this energy density is found to reflect the folding transition. Moreover, this approach allows us to build a reduced thermodynamics in the inherent structure landscape. Equilibrium studies, from full molecular dynamics (MD) simulations and from the reduced thermodynamics, detect the features of a dynamical transition at low temperature and we analyze the location and time scale of the fluctuations of the protein, showing the need of some frustration in the model to get realistic results. The frustrated model also shows the presence of a kinetic trap which strongly affects the dynamics of folding.

The inherent structure landscape of a protein

Using the Gō model of a real protein, we explore the landscape of its metastable structures. First, we show how the inherent structure energy density can be obtained from the probability density determined by sampling molecular dynamics trajectories and quenching. The analysis of the inherent structure landscape can characterize the folding transition. Then we show how thermodynamics of the inherent states can be established to study the equilibrium properties of proteins. Our work brings some elements into the current discussion about the protein dynamical transition. The study uses a simplified model to illustrate the ideas, but, as the inherent structure landscape is much simpler than the free energy surface of the protein, it appears to be accessible for an all-atom model of a small protein, at the expense of much longer calculations.

Metastable configurations of small-world networks

We calculate the number of metastable configurations of Ising small-world networks that are constructed upon superimposing sparse Poisson random graphs onto a one-dimensional chain. Our solution is based on replicated transfer-matrix techniques. We examine the denegeracy of the ground state and find a jump in the entropy of metastable configurations exactly at the crossover between the small-world and the Poisson random graph structures. We also examine the difference in entropy between metastable and all possible configurations, for both ferromagnetic and bond-disordered long-range couplings.

Exact occupation time distribution in a non-Markovian sequence and its relation to spin glass models

We compute exactly the distribution of the occupation time in a discrete non-Markovian toy sequence that appears in various physical contexts such as the diffusion processes and Ising spin glass chains. The non-Markovian property makes the results nontrivial even for this toy sequence. The distribution is shown to have non-Gaussian tails characterized by a nontrivial large deviation function that is computed explicitly. An exact mapping of this sequence to an Ising spin glass chain via a gauge transformation raises an interesting question for a generic finite sized spin glass model; at a given temperature, what is the distribution (over disorder) of the thermally averaged number of spins that are aligned to their local fields? We show that this distribution remains nontrivial even at infinite temperature and can be computed explicitly in few cases such as in the Sherrington-Kirkpatrick model with Gaussian disorder.

Metastable states in spin glasses and disordered ferromagnets

We study analytically M-spin-flip stable states in disordered short-ranged Ising models (spin glasses and ferromagnets) in all dimensions and for all M. Our approach is primarily dynamical, and is based on the convergence of sigma(t), a zero-temperature dynamical process with flips of lattice animals up to size M and starting from a deep quench, to a metastable limit sigma(infinity). The results (rigorous and nonrigorous, in infinite and finite volumes) concern many aspects of metastable states: their numbers, basins of attraction, energy densities, overlaps, remanent magnetizations, and relations to thermodynamic states. For example, we show that their overlap distribution is a delta function at zero. We also define a dynamics for M=infinity, which provides a potential tool for investigating ground state structure.