Replica field theory for deterministic models. II. A non-random spin glass with glassy behaviour

Symmetry-based approach to species-rich ecological communities

Disordered systems theory provides powerful tools to analyze the generic behaviors of high-dimensional systems, such as species-rich ecological communities or neural networks. By assuming randomness in their interactions, universality ensures that many microscopic details are irrelevant to system-wide dynamics; but the choice of a random ensemble still limits the generality of results. We show here, in the context of ecological dynamics, that these analytical tools do not require a specific choice of ensemble and that solutions can be found based only on a fundamental rotational symmetry in the interactions, encoding the idea that traits can be recombined into new species without altering global features. Dynamical outcomes then depend on the spectrum of the interaction matrix as a free parameter, allowing us to bridge between results found in different models of interactions and extend beyond them to previously unidentified behaviors. The distinctive feature of ecological models is the possibility of species extinctions, which leads to an increased universality of dynamics as the fraction of extinct species increases. We expect these findings can inform new developments in theoretical ecology as well as other families of complex systems.

Equilibrium fluctuations in mean-field disordered models

Mean-field models of glasses that present a random first order transition exhibit highly nontrivial fluctuations. Building on previous studies that focused on the critical scaling regime, we here obtain a fully quantitative framework for all equilibrium conditions. By means of the replica method we evaluate Gaussian fluctuations of the overlaps around the thermodynamic limit, decomposing them in thermal fluctuations inside each state and heterogeneous fluctuations between different states. We first test and compare our analytical results with numerical simulation results for the p-spin spherical model and the random orthogonal model, and then analyze the random Lorentz gas. In all cases, a strong quantitative agreement is obtained. Our analysis thus provides a robust scheme for identifying the key finite-size (or finite-dimensional) corrections to the mean-field treatment of these paradigmatic glass models.

Self-Induced Glassy Phase in Multimodal Cavity Quantum Electrodynamics

We provide strong evidence that the effective spin-spin interaction in a multimodal confocal optical cavity gives rise to a self-induced glassy phase, which emerges exclusively from the peculiar Euclidean correlations and is not related to the presence of disorder as in standard spin glasses. As recently shown, this spin-spin effective interaction is both nonlocal and nontranslational invariant, and randomness in the atoms' positions produces a spin glass phase. Here we consider the simplest feasible disorder-free setting, where atoms form a one-dimensional regular chain and we study the thermodynamics of the resulting effective Ising model. We present extensive results showing that the system has a low-temperature glassy phase. The model depends on the adimensional parameter α=(a/w_{0})^{2}, a being a lattice spacing and w_{0} an interaction length scale. Notably, for rational values of α=p/q, the number of metastable states at low temperature grows exponentially with q and the problem of finding the ground state rapidly becomes computationally intractable, suggesting that the system develops high-energy barriers and ergodicity breaking occurs.

Exact solution to the random sequential dynamics of a message passing algorithm

We analyze the random sequential dynamics of a message passing algorithm for Ising models with random interactions in the large system limit. We derive exact results for the two-time correlation functions and the speed of convergence. The de Almedia-Thouless stability criterion of the static problem is found to be necessary and sufficient for the global convergence of the random sequential dynamics.

Optimal learning with excitatory and inhibitory synapses

Characterizing the relation between weight structure and input/output statistics is fundamental for understanding the computational capabilities of neural circuits. In this work, I study the problem of storing associations between analog signals in the presence of correlations, using methods from statistical mechanics. I characterize the typical learning performance in terms of the power spectrum of random input and output processes. I show that optimal synaptic weight configurations reach a capacity of 0.5 for any fraction of excitatory to inhibitory weights and have a peculiar synaptic distribution with a finite fraction of silent synapses. I further provide a link between typical learning performance and principal components analysis in single cases. These results may shed light on the synaptic profile of brain circuits, such as cerebellar structures, that are thought to engage in processing time-dependent signals and performing on-line prediction.

Memory-free dynamics for the Thouless-Anderson-Palmer equations of Ising models with arbitrary rotation-invariant ensembles of random coupling matrices

We propose an iterative algorithm for solving the Thouless-Anderson-Palmer equations of Ising models with arbitrary rotation-invariant (random) coupling matrices. In the thermodynamic limit, we prove by means of the dynamical functional method that the proposed algorithm converges when the so-called de Almeida Thouless criterion is fulfilled. Moreover, we give exact analytical expressions for the rate of the convergence.

Classical many-particle systems with unique disordered ground states

Classical ground states (global energy-minimizing configurations) of many-particle systems are typically unique crystalline structures, implying zero enumeration entropy of distinct patterns (aside from trivial symmetry operations). By contrast, the few previously known disordered classical ground states of many-particle systems are all high-entropy (highly degenerate) states. Here we show computationally that our recently proposed "perfect-glass" many-particle model [Sci. Rep. 6, 36963 (2016)10.1038/srep36963] possesses disordered classical ground states with a zero entropy: a highly counterintuitive situation . For all of the system sizes, parameters, and space dimensions that we have numerically investigated, the disordered ground states are unique such that they can always be superposed onto each other or their mirror image. At low energies, the density of states obtained from simulations matches those calculated from the harmonic approximation near a single ground state, further confirming ground-state uniqueness. Our discovery provides singular examples in which entropy and disorder are at odds with one another. The zero-entropy ground states provide a unique perspective on the celebrated Kauzmann-entropy crisis in which the extrapolated entropy of a supercooled liquid drops below that of the crystal. We expect that our disordered unique patterns to be of value in fields beyond glass physics, including applications in cryptography as pseudorandom functions with tunable computational complexity.

Large fluctuations at the lasing threshold of solid- and liquid-state dye lasers

Intensity fluctuations in lasers are commonly studied above threshold in some special configurations (especially when emission is fed back into the cavity or when two lasers are coupled) and related with their chaotic behaviour. Similar fluctuating instabilities are usually observed in random lasers, which are open systems with plenty of quasi-modes whose non orthogonality enables them to exchange energy and provides the sort of loss mechanism whose interplay with pumping leads to replica symmetry breaking. The latter however, had never been observed in plain cavity lasers where disorder is absent or not intentionally added. Here we show a fluctuating lasing behaviour at the lasing threshold both in solid and liquid dye lasers. Above and below a narrow range around the threshold the spectral line-shape is well correlated with the pump energy. At the threshold such correlation disappears, and the system enters a regime where emitted laser fluctuates between narrow, intense and broad, weak peaks. The immense number of modes and the reduced resonator quality favour the coupling of modes and prepares the system so that replica symmetry breaking occurs without added disorder.

Scaling collapse at the jamming transition

The jamming transition of particles with finite-range interactions is characterized by a variety of critical phenomena, including power-law distributions of marginal contacts. We numerically study a recently proposed simple model of jamming, which is conjectured to lie in the same universality class as the jamming of spheres in all dimensions. We extract numerical estimates of the critical exponents, θ=0.451±0.006 and γ=0.404±0.004, that match the exponents observed in sphere packing systems. We analyze finite-size scaling effects that manifest in a subcritical cutoff regime and size-independent but protocol-dependent scaling curves. Our results support the conjectured link with sphere jamming, provide more precise measurements of the critical exponents than previously reported, and shed light on the finite-size scaling behavior of continuous constraint satisfiability transitions.

Replica-symmetry-breaking transitions and off-equilibrium dynamics

I consider branches of replica-symmetry-breaking (RSB) solutions in glassy systems that display a dynamical transition at a temperature T_{d} characterized by a mode-coupling-theory dynamical behavior. Below T_{d} these branches of solutions are considered to be relevant to the system complexity and to off-equilibrium dynamics. Under general assumptions I argue that near T_{d} it is not possible to stabilize the one-step (1RSB) solution beyond the marginal point by making a full RSB (FRSB) ansatz. However, depending on the model, there may exist a temperature T strictly lower than T_{d} below which the 1RSB branch can be continued to a FRSB branch. Such a temperature certainly exists for models that display the so-called Gardner transition and in this case T_{G}<T_<T_{d}. An analytical study in the context of the truncated model reveals that the FRSB branch of solutions below T is characterized by a two-plateau structure and it ends where the first plateau disappears. These general features are confirmed in the context of the Ising p-spin model with p=3 by means of a numerical solution of the FRSB equations. The results are discussed in connection with off-equilibrium dynamics within Cugliandolo-Kurchan theory. In this context I assume that the RSB solution relevant for off-equilibrium dynamics is the 1RSB marginal solution in the whole range (T ,T_{d}) and it is the end point of the FRSB branch for T<T. Remarkably, under these assumptions it can be argued that T marks a qualitative change in off-equilibrium dynamics in the sense that the decay of various dynamical quantities changes from power law to logarithmic.

RKKY interaction and intrinsic frustration in non-Fermi-liquid metals

We study the RKKY interaction in non-Fermi-liquid metals. We find that the RKKY interaction mediated by some non-Fermi-liquid metals can be of much longer range than for a Fermi liquid. The oscillatory nature of the RKKY interaction thus becomes more important in such non-Fermi liquids, and gives rise to enhanced frustration when the spins form a lattice. Frustration suppresses the magnetic ordering temperature of the lattice spin system. Furthermore, we find that the spin system with a longer range RKKY interaction can be described by the Brazovskii model, where the ordering wave vector lies on a higher dimensional manifold. Strong fluctuations in such a model lead to a first-order phase transition and/or glassy phase. This may explain some recent experiments where glassy behavior was observed in stoichiometric heavy fermion material close to a ferromagnetic quantum critical point.

Metastable states and space-time phase transitions in a spin-glass model

We study large deviations of the dynamical activity in the random orthogonal model. This is a fully connected spin-glass model with one-step replica symmetry-breaking behavior, consistent with the random first-order transition scenario for structural glasses. We show that this model displays dynamical (space-time) phase transitions between active and inactive phases, as demonstrated by singularities in large deviation functions. We argue that such transitions are generic in systems with long-lived metastable states.

Statistical mechanical analysis of the Kronecker channel model for multiple-input multiple-output wireless communication

The Kronecker channel model of wireless communication is analyzed using statistical mechanics methods. In the model, spatial proximities among transmission/reception antennas are taken into account as certain correlation matrices, which generally yield nontrivial dependence among symbols to be estimated. This prevents accurate assessment of the communication performance by naively using a previously developed analytical scheme based on a matrix integration formula. In order to resolve this difficulty, we develop a formalism that can formally handle the correlations in Kronecker models based on the known scheme. Unfortunately, direct application of the developed scheme is, in general, practically difficult. However, the formalism is still useful, indicating that the effect of the correlations generally increase after the fourth order with respect to correlation strength. Therefore, the known analytical scheme offers a good approximation in performance evaluation when the correlation strength is sufficiently small. For a class of specific correlation, we show that the performance analysis can be mapped to the problem of one-dimensional spin systems in random fields, which can be investigated without approximation by the belief propagation algorithm.

Curvature-induced frustration in the XY model on hyperbolic surfaces

We study low-temperature properties of the XY spin model on a negatively curved surface. Geometric curvature of the surface gives rise to frustration in local spin configuration, which results in the formation of high-energy spin clusters scattered over the system. Asymptotic behavior of the spin-glass susceptibility suggests a zero-temperature glass transition, which is attributed to multiple optimal configurations of spin clusters due to nonzero surface curvature of the system. It implies that a constant ferromagnetic spin interaction on a regular lattice can exhibit glasslike behavior without possessing any disorder if the lattice is put on top of a negatively curved space such as a hyperbolic surface.

Mode-coupling theory for the slow collective dynamics of fluids adsorbed in disordered porous media

We derive a mode-coupling theory for the slow dynamics of fluids confined in disordered porous media represented by spherical particles randomly placed in space. Its equations display the usual nonlinear structure met in this theoretical framework, except for a linear contribution to the memory kernel which adds to the usual quadratic term. The coupling coefficients involve structural quantities which are specific of fluids evolving in random environments and have expressions which are consistent with those found in related problems. Numerical solutions for two simple models with pure hard core interactions lead to the prediction of a variety of glass transition scenarios, which are either continuous or discontinuous and include the possibility of higher-order singularities and glass-glass transitions. The main features of the dynamics in the two most generic cases are reviewed and illustrated with detailed computations. Moreover, a reentry phenomenon is predicted in the low fluid-high matrix density regime and is interpreted as the signature of a de-correlation mechanism by fluid-fluid collisions competing with the localization effect of the solid matrix.

Spin glasses and fragile glasses: statics, dynamics, and complexity

In this paper I will briefly review some theoretical results that have been obtained in recent years for spin glasses and fragile glasses. I will concentrate my attention on the predictions coming from the so called broken replica symmetry approach and on their experimental verifications. I will also mention the relevance or these results for other fields, and in general for complex systems.

Jamming percolation and glass transitions in lattice models

A new class of lattice gas models with trivial interactions but constrained dynamics is introduced. These models are proven to exhibit a dynamical glass transition: above a critical density rhoc ergodicity is broken due to the appearance of an infinite spanning cluster of jammed particles. The fraction of jammed particles is discontinuous at the transition, while in the unjammed phase dynamical correlation lengths and time scales diverge as exp[C(rhoc-rho)-mu]. Dynamic correlations display two-step relaxation similar to glass formers and jamming systems.

Shear-thickening and entropy-driven reentrance

We discuss a generic mechanism for shear thickening analogous to entropy-driven phase reentrance. We implement it in the context of nonrelaxational mean-field glassy systems: although very simple, the microscopic models we study present a dynamical phase diagram with second- and first-order stirring-induced jamming transitions leading to intermittency, metastability, and phase coexistence as seen in some experiments. The jammed state is fragile with respect to change in the stirring direction. Our approach provides a direct derivation of a mode-coupling theory of shear thickening.

Stable solution of the simplest spin model for inverse freezing

We analyze the Blume-Emery-Griffiths-Capel model with disordered interaction that displays the inverse freezing phenomenon. The behavior of this spin-1 model in crystal field is studied throughout the phase diagram, and the transition lines are computed using the full replica symmetry breaking ansatz. We compare the results both with the formulation of the same model in terms of Ising spins on lattice gas, where no reentrance takes place, and with the model with generalized spin variables recently introduced by Schupper and Shnerb [Phys. Rev. Lett. 93, 037202 (2004)], for which the reentrance is enhanced as the ratio between the degeneracy of full to empty sites increases. The simplest version of all these models, known as the Ghatak-Sherrington model, turns out to hold all the general features characterizing an inverse transition to an amorphous phase.

Mean-field glass transition in a model liquid

We investigate the liquid-glass phase transition in a system of pointlike particles interacting via a finite-range attractive potential in D -dimensional space. The phase transition is driven by an "entropy crisis" where the available phase space volume collapses dramatically at the transition. We describe the general strategy underlying the first-principles replica calculation for this type of transition; its application to our model system then allows for an analytic description of the liquid-glass phase transition within a mean-field approximation, provided the parameters are chosen suitably. We find a transition exhibiting all the features associated with an entropy crisis, including the characteristic finite jump of the order parameter at the transition while the free energy and its first derivative remain continuous.

Glass transition in a simple stochastic model with back-reaction

We formulate and solve a model of dynamical arrest in colloids. A particle is coupled to the bath of statistically identical particles. The dynamics is described by Langevin equation with stochastic external force described by telegraphic noise. The interaction with the bath is taken into account self-consistently through the back-reaction mechanism. Dynamically induced glass transition occurs for certain value of the coupling strength. Edwards-Anderson parameter jumps discontinuously at the transition. Another order parameter can be also defined, which vanishes continuously with exponent 1/2 at the critical point. Nonlinear response to harmonic perturbation is found.

Role of the interaction matrix in mean-field spin glass models

Mean-field models of two-spin Ising spin glasses with interaction matrices taken from ensembles that are invariant under O(N) transformations are studied. A general study shows that the nature of the spin glass transition can be deduced from the eigenvalue spectrum of the interaction matrix. A simple replica approach is derived to carry out the average over the O(N) disorder. The analytic results are confirmed by the extensive Monte Carlo simulations for large system sizes and by the exact enumeration for small system sizes.

Tensionless structure of a glassy phase

We study a class of homogeneous finite-dimensional Ising models which were recently shown to exhibit glassy properties. Monte Carlo simulations of a particular three-dimensional model in this class show that the glassy phase obtained under slow cooling is dominated by large-scale excitations whose energy E(l) scales with their size l as E(l) approximately l(Theta) with Theta approximately 1.33(20). Simulations suggest that in another model of this class, namely the four-spin model, energy is concentrated mainly in linear defects, making the domain walls tensionless in this case also. Two-dimensional variants of these models are trivial and most likely the energy of excitations scales with the exponent Theta=1.

Adaptive and self-averaging Thouless-Anderson-Palmer mean-field theory for probabilistic modeling

We develop a generalization of the Thouless-Anderson-Palmer (TAP) mean-field approach of disorder physics, which makes the method applicable to the computation of approximate averages in probabilistic models for real data. In contrast to the conventional TAP approach, where the knowledge of the distribution of couplings between the random variables is required, our method adapts to the concrete set of couplings. We show the significance of the approach in two ways: Our approach reproduces replica symmetric results for a wide class of toy models (assuming a nonglassy phase) with given disorder distributions in the thermodynamic limit. On the other hand, simulations on a real data model demonstrate that the method achieves more accurate predictions as compared to conventional TAP approaches.

Crystallization of a supercooled liquid and of a glass: Ising model approach

Using Monte Carlo simulations we study crystallization in the three-dimensional Ising model with four-spin interaction. We monitor the morphology of crystals which grow after placing crystallization seeds in a supercooled liquid. Defects in such crystals constitute an intricate and very stable network that separates various domains by tensionless domain walls. We also show that the crystallization which occurs during the continuous heating of the glassy phase takes place at a heating-rate-dependent temperature.

Toy model for the mean-field theory of hard-sphere liquids

We investigate a toy model of liquid, based on simplified hypernetted chain equations in very large spatial dimension D. The model does not exhibit a phase transition, but several regimes of behavior when D-->infinity can be observed in different intervals of the density.

Dynamics of the frustrated ising lattice gas

The dynamical properties of a three-dimensional model glass, the frustrated Ising lattice gas (FILG), are studied by Monte Carlo simulations. We present results of compression experiments, where the chemical potential is either slowly or abruptly changed, as well as simulations at constant density. One time quantities like density and two time quantities like correlations, responses, and mean square displacements are measured, and the departure from equilibrium clearly characterized. The aging scenario, particularly in the case of density autocorrelations, is reminiscent of spin glass phenomenology with violations of the fluctuation-dissipation theorem, typical of systems with one replica symmetry breaking. The FILG, as a valid on-lattice model of structural glasses, can be described with tools developed in spin glass theory and, being a finite-dimensional model, can open the way for a systematic study of activated processes in glasses.

Slow dynamics of ising models with energy barriers

Using Monte Carlo simulations we study the dynamics of three-dimensional Ising models with nearest-, next-nearest-, and four-spin (plaquette) interactions. During coarsening, such models develop growing energy barriers, which leads to very slow dynamics at low temperature. As already reported, the model with only the plaquette interaction exhibits some of the features characteristic of ordinary glasses: strong metastability of the supercooled liquid, a weak increase of the characteristic length under cooling, stretched-exponential relaxation, and aging. The addition of two-spin interactions, in general, destroys such behavior: the liquid phase loses metastability and the slow-dynamics regime terminates well below the melting transition, which is presumably related with a certain corner-rounding transition. However, for a particular choice of interaction constants, when the ground state is strongly degenerate, our simulations suggest that the slow-dynamics regime extends up to the melting transition. The analysis of these models leads us to the conjecture that in the four-spin Ising model domain walls lose their tension at the glassy transition and that they are basically tensionless in the glassy phase.

Discrete spin variables and critical temperature in deterministic models with glassy behavior

The problem of the existence of a glassy phase transition in deterministic spin models is reconsidered, examining an Ising model with general spin s and nontranslationally invariant interaction. The discrete nature of the spin variables is shown to allow the glass state.

Cryptographical properties of Ising spin systems

The relation between Ising spin systems and public-key cryptography is investigated using methods of statistical physics. The insight gained from the analysis is used for devising a matrix-based cryptosystem whereby the ciphertext comprises products of the original message bits; these are selected by employing two predetermined randomly constructed sparse matrices. The ciphertext is decrypted using methods of belief propagation. The analyzed properties of the suggested cryptosystem show robustness against various attacks and competitive performance to modern cryptographical methods.

Glassy dynamics and aging in an exactly solvable spin model

We introduce a simple two-dimensional spin model with short-range interactions which shows glassy behavior despite a Hamiltonian which is completely homogeneous and possesses no randomness. We solve exactly for both the static partition function of the model and the distribution of energy barriers, giving us the equilibration time scales at low temperature. Simulations of instantaneous quenches and of annealing of the model are in good agreement with the analytic calculations. We also measure the two-time spin correlation as a function of waiting time, and show that the model has aging behavior consistent with the distribution of barrier heights. The model appears to have no sharp glass transition. Instead, it falls out of equilibrium at a temperature which decreases logarithmically as a function of the cooling time.

Damage spreading transition in glasses: a probe for the ruggedness of the configurational landscape

We consider damage spreading transitions in the framework of mode-coupling theory. This theory describes relaxation processes in glasses in the mean-field approximation which are known to be characterized by the presence of an exponentially large number of metastable states. For systems evolving under identical but arbitrarily correlated noises, we demonstrate that there exists a critical temperature T0 which separates two different dynamical regimes depending on whether damage spreads or not in the asymptotic long-time limit. This transition exists for generic noise correlations such that the zero damage solution is stable at high temperatures, being minimal for maximal noise correlations. Although this dynamical transition depends on the type of noise correlations, we show that the asymptotic damage has the good properties of a dynamical order parameter, such as (i) independence of the initial damage; (ii) independence of the class of initial condition; and (iii) stability of the transition in the presence of asymmetric interactions which violate detailed balance. For maximally correlated noises we suggest that damage spreading occurs due to the presence of a divergent number of saddle points (as well as metastable states) in the thermodynamic limit consequence of the ruggedness of the free-energy landscape which characterizes the glassy state. These results are then compared to extensive numerical simulations of a mean-field glass model (the Bernasconi model) with Monte Carlo heat-bath dynamics. The freedom of choosing arbitrary noise correlations for Langevin dynamics makes damage spreading an interesting tool to probe the ruggedness of the configurational landscape.