Search

The presence of a random magnetic field in ferromagnetic systems leads, in the broken phase, to an anomalous $O(\sqrt{1/N})$ convergence of some thermodynamic quantities to their asymptotic limits. Here we show a general method, based on the replica trick, to compute analytically the $O(\sqrt{1/N})$ finite size correction to the average free energy. We apply this method to two mean field Ising models, fully connected and random regular graphs, and compare the results to exact numerical algorithms. We argue that this behaviour is present in finite dimensional models as well.

We propose a new algorithm for inferring the state of hidden spins and reconstructing the connections in a synchronous kinetic Ising model, given the observed history. Focusing on the case in which the hidden spins are conditionally independent of each other given the state of observable spins, we show that calculating the likelihood of the data can be simplified by introducing a set of replicated auxiliary spins. Belief propagation (BP) and susceptibility propagation (SusP) can then be used to infer the states of hidden variables and to learn the couplings. We study the convergence and performance of this algorithm for networks with both Gaussian-distributed and binary bonds. We also study how the algorithm behaves as the fraction of hidden nodes and the amount of data are changed, showing that it outperforms the Thouless-Anderson-Palmer (TAP) equations for reconstructing the connections. QC 20150625

In many mean-field glassy systems, the low-temperature Gibbs measure is dominated by exponentially many metastable states. We analyze the evolution of the metastable states as temperature changes adiabatically in the solvable case of the spherical $s+p$-spin glass model, extending the work of Barrat, Franz and Parisi J. Phys. A 30, 5593 (1997). We confirm the presence of level crossings, bifurcations, and temperature chaos. For the states that are at equilibrium close to the so-called dynamical temperature $T_d$, we find, however, that the following state method (and the dynamical solution of the model as well) is intrinsically limited by the vanishing of solutions with non-zero overlap at low temperature.

An efficient microcanonical dynamics has been recently introduced for Ising spin models embedded in a generic connected graph even in the presence of disorder i.e. with the spin couplings chosen from a random distribution. Such a dynamics allows a coherent definition of local temperatures also when open boundaries are coupled to thermostats, imposing an energy flow. Within this framework, here we introduce a consistent definition for local energy currents and we study their dependence on the disorder. In the linear response regime, when the global gradient between thermostats is small, we also define local conductivities following a Fourier dicretized picture. Then, we work out a linearized "mean-field approximation", where local conductivities are supposed to depend on local couplings and temperatures only. We compare the approximated currents with the exact results of the nonlinear system, showing the reliability range of the mean-field approach, which proves very good at high temperatures and not so efficient in the critical region. In the numerical studies we focus on the disordered cylinder but our results could be extended to an arbitrary, disordered spin model on a generic discrete structures., Comment: 12 pages, 6 figures

The aim of this work is to try to bridge over theoretical immunology and disordered statistical mechanics. Our long term hope is to contribute to the development of a quantitative theoretical immunology from which practical applications may stem. In order to make theoretical immunology appealing to the statistical physicist audience we are going to work out a research article which, from one side, may hopefully act as a benchmark for future improvements and developments, from the other side, it is written in a very pedagogical way both from a theoretical physics viewpoint as well as from the theoretical immunology one. Furthermore, we have chosen to test our model describing a wide range of features of the adaptive immune response in only a paper: this has been necessary in order to emphasize the benefit available when using disordered statistical mechanics as a tool for the investigation. However, as a consequence, each section is not at all exhaustive and would deserve deep investigation: for the sake of completeness, we restricted details in the analysis of each feature with the aim of introducing a self-consistent model., Comment: 22 pages, 14 figure

For a non-interacting Bose gas on a lattice we compute the shift of the critical temperature for condensation when random-bond and onsite disorder are present. We evidence that the shift depends on the space dimensionality D and the filling fraction f. For D -> infinity (infinite-range model), using results from the theory of random matrices, we show that the shift of the critical temperature is negative, depends on f, and vanishes only for large f. The connections with analogous results obtained for the spherical model are discussed. For D=3 we find that, for large f, the critical temperature Tc is enhanced by disorder and that the relative shift does not sensibly depend on f; at variance, for small f, Tc decreases in agreement with the results obtained for a Bose gas in the continuum. We also provide numerical estimates for the shift of the critical temperature due to disorder induced on a non-interacting Bose gas by a bichromatic incommensurate potential., 18 pages, 8 figures; Fig. 8 improved adding results for another value of q (q=830/1076)

We present a finite-size scaling analysis of high-statistics Monte Carlo simulations of the three-dimensional randomly site-diluted and bond-diluted Ising model. The critical behavior of these systems is affected by slowly-decaying scaling corrections which make the accurate determination of their universal asymptotic behavior quite hard, requiring an effective control of the scaling corrections. For this purpose we exploit improved Hamiltonians, for which the leading scaling corrections are suppressed for any thermodynamic quantity, and improved observables, for which the leading scaling corrections are suppressed for any model belonging to the same universality class. The results of the finite-size scaling analysis provide strong numerical evidence that phase transitions in three-dimensional randomly site-diluted and bond-diluted Ising models belong to the same randomly dilute Ising universality class. We obtain accurate estimates of the critical exponents, $\nu=0.683(2)$, $\eta=0.036(1)$, $\alpha=-0.049(6)$, $\gamma=1.341(4)$, $\beta=0.354(1)$, $\delta=4.792(6)$, and of the leading and next-to-leading correction-to-scaling exponents, $\omega=0.33(3)$ and $\omega_2=0.82(8)$., Comment: 45 pages, 22 figs, revised estimate of nu

We present results of a Monte Carlo study of the equilibrium dynamics of the one dimensional long-range Ising spin glass model. By tuning a parameter $\sigma$, this model interpolates between the mean field Sherrington-Kirkpatrick model and a proxy of the finite dimensional Edward-Anderson model. Activated scaling fits for the behavior of the relaxation time $\tau$ as a function of the number of spins $N$ (Namely $\ln(\tau)\propto N^{\psi}$) give values of $\psi$ that are not stable against inclusion of subleading corrections. Critical scaling ($\tau\propto N^{\rho}$) gives more stable fits, at least in the non mean field region. We also present results on the scaling of the time decay of the critical remanent magnetization of the Sherrington-Kirkpatrick model, a case where the simulation can be done with quite large systems and that shows the difficulties in obtaining precise values for dynamical exponents in spin glass models.

We calculate the optimal solutions of the fully heterogeneous Von Neumann expansion problem with $N$ processes and $P$ goods in the limit $N\to\infty$. This model provides an elementary description of the growth of a production economy in the long run. The system turns from a contracting to an expanding phase as $N$ increases beyond $P$. The solution is characterized by a universal behavior, independent of the parameters of the disorder statistics. Associating technological innovation to an increase of $N$, we find that while such an increase has a large positive impact on long term growth when $N\ll P$, its effect on technologically advanced economies ($N\gg P$) is very weak., Comment: 8 pages, 1 figure

Crystalline plasticity is strongly interlinked with dislocation mechanics and nowadays is relatively well understood. Concepts and physical models of plastic deformation in amorphous materials on the other hand - where the concept of linear lattice defects is not applicable - still are lagging behind. We introduce an eigenstrain-based finite element lattice model for simulations of shear band formation and strain avalanches. Our model allows us to study the influence of surfaces and finite size effects on the statistics of avalanches. We find that even with relatively complex loading conditions and open boundary conditions, critical exponents describing avalanche statistics are unchanged, which validates the use of simpler scalar lattice-based models to study these phenomena., Comment: Journal of Statistical Mechanics: Theory and Experiment, 2015, P02011

We present a method derived by cavity arguments to compute the spin-glass and higher-order susceptibilities in diluted mean-field spin-glass models. The divergence of the spin-glass susceptibility is associated to the existence of a non-zero solution of a homogeneous linear integral equation. Higher order susceptibilities, relevant for critical dynamics through the parameter exponent $\lambda$, can be expressed at criticality as integrals involving the critical eigenvector. The numerical evaluation of the corresponding analytic expressions is discussed. The method is illustrated in the context of the de Almeida-Thouless line for a spin-glass on a Bethe lattice but can be generalized straightforwardly to more complex situations., Comment: v3: typo error in eq. (10) (and similar equations) corrected

We study the performance of different message passing algorithms in the two dimensional Edwards Anderson model. We show that the standard Belief Propagation (BP) algorithm converges only at high temperature to a paramagnetic solution. Then, we test a Generalized Belief Propagation (GBP) algorithm, derived from a Cluster Variational Method (CVM) at the plaquette level. We compare its performance with BP and with other algorithms derived under the same approximation: Double Loop (DL) and a two-ways message passing algorithm (HAK). The plaquette-CVM approximation improves BP in at least three ways: the quality of the paramagnetic solution at high temperatures, a better estimate (lower) for the critical temperature, and the fact that the GBP message passing algorithm converges also to non paramagnetic solutions. The lack of convergence of the standard GBP message passing algorithm at low temperatures seems to be related to the implementation details and not to the appearance of long range order. In fact, we prove that a gauge invariance of the constrained CVM free energy can be exploited to derive a new message passing algorithm which converges at even lower temperatures. In all its region of convergence this new algorithm is faster than HAK and DL by some orders of magnitude., Comment: 19 pages, 13 figures

During the last years, through the combined effort of the insight, coming from physical intuition and computer simulation, and the exploitation of rigorous mathematical methods, the main features of the mean field Sherrington-Kirkpatrick spin glass model have been firmly established. In particular, it has been possible to prove the existence and uniqueness of the infinite volume limit for the free energy, and its Parisi expression, in terms of a variational principle, involving a functional order parameter. Even the expected property of ultrametricity, for the infinite volume states, seems to be near to a complete proof. The main structural feature of this model, and related models, is the deep phenomenon of spontaneous replica symmetry breaking (RSB), discovered by Parisi many years ago. By expanding on our previous work, the aim of this paper is to investigate a general frame, where replica symmetry breaking is embedded in a kind of mechanical scheme of the Hamilton-Jacobi type. Here, the analog of the "time" variable is a parameter characterizing the strength of the interaction, while the "space" variables rule out quantitatively the broken replica symmetry pattern. Starting from the simple cases, where annealing is assumed, or replica symmetry, we build up a progression of dynamical systems, with an increasing number of space variables, which allow to weaken the effect of the potential in the Hamilton-Jacobi equation, as the level of symmetry braking is increased. This new machinery allows to work out mechanically the general K-step RSB solutions, in a different interpretation with respect to the replica trick, and lightens easily their properties as existence or uniqueness., Comment: 24 pages, no figures

What is the origin of the sharp slowdown displayed by glassy systems? Physical common sense suggests there must be a concomitant growing correlation length, but finding this length has been nontrivial. In random first-order theory, it is given by the size of amorphous excitations, which depends on a balance between their mutual interfacial energy and their configurational entropy. But how these excitations disappear when crossing over to the normal high temperature phase is unclear, chiefly due to lack of data about the surface tension. We measure the energy cost for creating amorphous excitations in a model glass-former, and discover that the surface tension vanishes at a well-defined spinodal energy, above which amorphous excitations cannot be sustained. This spinodal therefore marks the true onset of glassiness., Facultad de Ciencias Exactas, Instituto de Investigaciones Fisicoquímicas Teóricas y Aplicadas

In this paper we obtain two results for the Sherrington-Kirkpatrick (SK) model, and we show that they both emerge from a single approach. First, we prove that the average of the overlap takes positive values when it is non zero. More specificly, the average of the overlap, which is naively expected to take values in the whole interval $[-1,+1]$, becomes positive if we ``first'' apply an external field, so to destroy the gauge invariance of the model, and ``then'' remove it in the thermodynamic limit. This phenomenon emerges at the critical point. This first result is weaker that the one obtained by Talagrand (not limited to the average of the overlap), but we show here that, at least in average, the overlap is proven to be non-negative with no use of the Ghirlanda-Guerra identities. The latter are instead needed to obtain the second result, which is the control the behavior of the overlap at the critical point: we find the critical exponents of all the overlap correlation functions., Comment: 19 pages