Loss of Ergodicity in the Transition from Annealed to Quenched Disorder in a Finite Kinetic Ising Model

We consider a kinetic Ising model which represents a generic agent-based model for various types of socio-economic systems. We study the case of a finite (and not necessarily large) number of agents N as well as the asymptotic case when the number of agents tends to infinity. The main ingredient are individual decision thresholds which are either fixed over time (corresponding to quenched disorder in the Ising model, leading to nonlinear deterministic dynamics which are generically non-ergodic) or which may change randomly over time (corresponding to annealed disorder, leading to ergodic dynamics). We address the question how increasing the strength of annealed disorder relative to quenched disorder drives the system from non-ergodic behavior to ergodicity. Mathematically rigorous analysis provides an explicit and detailed picture for arbitrary realizations of the quenched initial thresholds, revealing an intriguing ""jumpy"" transition from non-ergodicity with many absorbing sets to ergodicity. For large N we find a critical strength of annealed randomness, above which the system becomes asymptotically ergodic. Our theoretical results suggests how to drive a system from an undesired socio-economic equilibrium (e. g. high level of corruption) to a desirable one (low level of corruption).

Anomalous roughening of Hele-Shaw flows with quenched disorder

The kinetic roughening of a stable oil-air interface moving in a Hele-Shaw cell that contains a quenched columnar disorder (tracks) has been studied. A capillary effect is responsible for the dynamic evolution of the resulting rough interface, which exhibits anomalous scaling. The three independent exponents needed to characterize the anomalous scaling are determined experimentally. The anomalous scaling is explained in terms of the initial acceleration and subsequent deceleration of the interface tips in the tracks coupled by mass conservation. A phenomenological model that reproduces the measured global and local exponents is introduced.

Fluctuations in Saffman-Taylor fingers with quenched disorder

We make an experimental characterization of the effect that static disorder has on the shape of a normal Saffman-Taylor finger. We find that static noise induces a small amplitude and long wavelength instability on the sides of the finger. Fluctuations on the finger sides have a dominant wavelength, indicating that the system acts as a selective amplifier of static noise. The dominant wavelength does not seem to be very sensitive to the intensity of static noise present in the system. On the other hand, at a given flow rate, rms fluctuations of the finger width, decrease with decreasing intensity of static noise. This might explain why the sides of the fingers are flat for typical Saffman-Taylor experiments. Comparison with previous numerical studies of the effect that temporal noise has on the Saffman-Taylor finger, leads to conclude that the effect of temporal noise and static noise are similar. The behavior of fluctuations of the finger width found in our experiments, is qualitatively similar to one recently reported, in the sense that, the magnitude of the width fluctuations decays as a power law of the capillary number, at low flow rates, and increases with capillary number for larger flow rates.

Universality in models for disorder induced phase transitions

We have systematically analyzed six different reticular models with quenched disorder and no thermal fluctuations exhibiting a field-driven first-order phase transition. We have studied the nonequilibrium transition, appearing when varying the amount of disorder, characterized by the change from a discontinuous hysteresis cycle (with one or more large avalanches) to a smooth one (with only tiny avalanches). We have computed critical exponents using finite size scaling techniques and shown that they are consistent with universal values depending only on the space dimensionality d.

Wave propagation in a medium with disordered excitability

The effect of quenched disorder on the propagation of autowaves in excitable media is studied both experimentally and numerically in relation to the light-sensitive Belousov-Zhabotinsky reaction. The spatial disorder is introduced through a random distribution with two different levels of transmittance. In one dimension the (time-averaged) wave speed is smaller than the corresponding to a homogeneous medium with the mean excitability. Contrarily, in two dimensions the velocity increases due to the roughening of the front. Results are interpreted using kinematic and scaling arguments. In particular, for d = 2 we verify a theoretical prediction of a power-law dependence for the relative change of the propagation speed on the disorder amplitude.

Growth of unstable interfaces in disordered media

The effects of a disordered medium in the growth of unstable interfaces are studied by means of two local models with multiplicative and additive quenched disorder, respectively. For short times and large pushing the multiplicative quenched disorder is equivalent to a time-dependent noise. In this regime, the linear dispersion relation contains a destabilizing contribution introduced by the noise. For long times, the interface always gets pinned. We model the systematics of the pinned shapes by means of an effective nonlinear model. These results show good agreement with numerical simulations. For the additive noise we find numerically that a depinning transition occurs.

Wave propagation in a medium with disordered excitability

The effect of quenched disorder on the propagation of autowaves in excitable media is studied both experimentally and numerically in relation to the light-sensitive Belousov-Zhabotinsky reaction. The spatial disorder is introduced through a random distribution with two different levels of transmittance. In one dimension the (time-averaged) wave speed is smaller than the corresponding to a homogeneous medium with the mean excitability. Contrarily, in two dimensions the velocity increases due to the roughening of the front. Results are interpreted using kinematic and scaling arguments. In particular, for d = 2 we verify a theoretical prediction of a power-law dependence for the relative change of the propagation speed on the disorder amplitude.

Critical behavior of a contact process with aperiodic transition rates

We performed Monte Carlo simulations to investigate the steady-state critical behavior of a one-dimensional contact process with an aperiodic distribution of rates of transition. As in the presence of randomness, spatial fluctuations can lead to changes of critical behavior. For sufficiently weak fluctuations, we give numerical evidence to show that there is no departure from the universal critical behavior of the underlying uniform model. For strong spatial fluctuations, the analysis of the data indicates a change of critical universality class.

Rounding of a first-order quantum phase transition to a strong-coupling critical point

We investigate the effects of quenched disorder on first-order quantum phase transitions on the example of the N-color quantum Ashkin-Teller model. By means of a strong-disorder renormalization group, we demonstrate that quenched disorder rounds the first-order quantum phase transition to a continuous one for both weak and strong coupling between the colors. In the strong-coupling case, we find a distinct type of infinite-randomness critical point characterized by additional internal degrees of freedom. We investigate its critical properties in detail and find stronger thermodynamic singularities than in the random transverse field Ising chain. We also discuss the implications for higher spatial dimensions as well as unusual aspects of our renormalization-group scheme. DOI: 10.1103/PhysRevB.86.214204

First-order transition in a three-dimensional disordered system

We present the first detailed numerical study in three dimensions of a first-order phase transition that remains first order in the presence of quenched disorder (specifically, the ferromagnetic-paramagnetic transition of the site-diluted four states Potts model). A tricritical point, which lies surprisingly near the pure-system limit and is studied by means of finite-size scaling, separates the first-order and second-order parts of the critical line. This investigation has been made possible by a new definition of the disorder average that avoids the diverging-variance probability distributions that plague the standard approach. Entropy, rather than free energy, is the basic object in this approach that exploits a recently introduced microcanonical Monte Carlo method.

First Order Phase Transition in a 3D disordered system

We present a detailed numerical study on the effects of adding quenched impurities to a three dimensional system which in the pure case undergoes a strong first order phase transition (specifically, the ferromagnetic/paramagnetic transition of the site-diluted four states Potts model). We can state that the transition remains first-order in the presence of quenched disorder (a small amount of it) but it turns out to be second order as more impurities are added. A tricritical point, which is studied by means of Finite-Size Scaling, separates the first-order and second-order parts of the critical line. The results were made possible by a new definition of the disorder average that avoids the diverging-variance probability distributions that arise using the standard methodology. We also made use of a recently proposed microcanonical Monte Carlo method in which entropy, instead of free energy, is the basic quantity.

Phase transitions and memory effects in the dynamics of Boolean networks

The generating functional method is employed to investigate the synchronous dynamics of Boolean networks, providing an exact result for the system dynamics via a set of macroscopic order parameters. The topology of the networks studied and its constituent Boolean functions represent the system's quenched disorder and are sampled from a given distribution. The framework accommodates a variety of topologies and Boolean function distributions and can be used to study both the noisy and noiseless regimes; it enables one to calculate correlation functions at different times that are inaccessible via commonly used approximations. It is also used to determine conditions for the annealed approximation to be valid, explore phases of the system under different levels of noise and obtain results for models with strong memory effects, where existing approximations break down. Links between Boolean networks and general Boolean formulas are identified and results common to both system types are highlighted. © 2012 Copyright Taylor and Francis Group, LLC.

Statistical physics of error-correcting codes

In this thesis we use statistical physics techniques to study the typical performance of four families of error-correcting codes based on very sparse linear transformations: Sourlas codes, Gallager codes, MacKay-Neal codes and Kanter-Saad codes. We map the decoding problem onto an Ising spin system with many-spins interactions. We then employ the replica method to calculate averages over the quenched disorder represented by the code constructions, the arbitrary messages and the random noise vectors. We find, as the noise level increases, a phase transition between successful decoding and failure phases. This phase transition coincides with upper bounds derived in the information theory literature in most of the cases. We connect the practical decoding algorithm known as probability propagation with the task of finding local minima of the related Bethe free-energy. We show that the practical decoding thresholds correspond to noise levels where suboptimal minima of the free-energy emerge. Simulations of practical decoding scenarios using probability propagation agree with theoretical predictions of the replica symmetric theory. The typical performance predicted by the thermodynamic phase transitions is shown to be attainable in computation times that grow exponentially with the system size. We use the insights obtained to design a method to calculate the performance and optimise parameters of the high performance codes proposed by Kanter and Saad.