Finite-size scaling in Ising-like systems with quenched random fields: evidence of hyperscaling violation

Random Field Ising Model Criticality in a Complex Binary Liquid System

While Ising criticality in classical liquids has been firmly established both theoretically and experimentally, much less is known about criticality in liquids in which the growth of the correlation length is frustrated by finite-size effects. A theoretical approach for dealing with this issue is the random-field Ising model (RFIM). While experimental critical-exponent values have been reported for magnetic samples (here, we consider γ, ν and η), little experimental information is available for critical fluctuations in corresponding liquid systems. In this paper, we present a study on a binary liquid consisting of 3-methyl pyridine and heavy water in a very light-weight porous gel. We find that the experimental results are in agreement with the theoretical predictions from the RFIM.

Breakdown of Ergodicity and Self-Averaging in Polar Flocks with Quenched Disorder

We show that spatial quenched disorder affects polar active matter in ways more complex and far reaching than heretofore believed. Using simulations of the 2D Vicsek model subjected to random couplings or a disordered scattering field, we find in particular that ergodicity is lost in the ordered phase, the nature of which we show to depend qualitatively on the type of quenched disorder: for random couplings, it remains long-range ordered, but qualitatively different from the pure (disorderless) case. For random scatterers, polar order varies with system size but we find strong non-self-averaging, with sample-to-sample fluctuations dominating asymptotically, which prevents us from elucidating the asymptotic status of order.

Random-field Ising model criticality in a glass-forming liquid

We use computer simulations to investigate the extended phase diagram of a supercooled liquid linearly coupled to a quenched reference configuration. An extensive finite-size scaling analysis demonstrates the existence of a random-field Ising model (RFIM) critical point and of a first-order transition line, in agreement with recent field-theoretical approaches. The dynamics in the vicinity of this critical point resembles the peculiar activated scaling of RFIM-like systems, and the overlap autocorrelation displays a logarithmic stretching. Our study demonstrates RFIM criticality in the thermodynamic limit for a three-dimensional supercooled liquid at equilibrium.

Elastic Backbone Defines a New Transition in the Percolation Model

The elastic backbone is the set of all shortest paths. We found a new phase transition at p_{eb} above the classical percolation threshold at which the elastic backbone becomes dense. At this transition in 2D, its fractal dimension is 1.750±0.003, and one obtains a novel set of critical exponents β_{eb}=0.50±0.02, γ_{eb}=1.97±0.05, and ν_{eb}=2.00±0.02, fulfilling consistent critical scaling laws. Interestingly, however, the hyperscaling relation is violated. Using Binder's cumulant, we determine, with high precision, the critical probabilities p_{eb} for the triangular and tilted square lattice for site and bond percolation. This transition describes a sudden rigidification as a function of density when stretching a damaged tissue.

Determination of the conformal-field-theory central charge by the Wang-Landau algorithm

We present a simple method to estimate the central charge of the conformal field theory corresponding to a critical point of a two-dimensional lattice model from Monte Carlo simulations. The main idea is to use the Wang-Landau flat-histogram algorithm, which allows us to obtain the free energy of a lattice model on a torus as a function of torus radii. The central charge is calculated with good precision from a free-energy scaling at the critical point. We apply the method to the Ising, tricritical Ising (Blume-Capel), Potts, and site-diluted Ising models, and we also discuss an estimation of the conformal weights.

Numerical evidence of hyperscaling violation in wetting transitions of the random-bond Ising model in d=2 dimensions

We performed extensive simulations of the random-bond Ising model confined between walls where competitive surface fields act. By properly taking the thermodynamic limit we unambiguously determined wetting transition points of the system, as extrapolation of localization-delocalization transitions of the interface between domains of different orientation driven by the respective fields. The finite-size scaling theory for wetting with short-range fields [E. V. Albano and K. Binder, Phys. Rev. Lett. 109, 036101 (2012)PRLTAO0031-900710.1103/PhysRevLett.109.036101] establishes that the average magnetization of the sample, with critical exponent β, is the proper order parameter for the study of wetting. While the hyperscaling relationship given by γ+2β=ν_{∥}+ν_{⊥} requires β=1/2 (γ=4, ν_{∥}=3, and ν_{⊥}=2), the thermodynamic scaling establishes that Δ_{s}=γ+β, which in contrast requires β=0 (Δ_{s}=4), where γ, ν_{∥}, ν_{⊥}, and Δ_{s} are the critical exponents of the susceptibility, the correlation lengths parallel and perpendicular to the interface, and the gap exponent, respectively. So, we formulate a finite-size scaling theory for wetting without hyperscaling and perform numerical simulations that provide strong evidence of hyperscaling violation (i.e., β=0) and a direct measurement of the susceptibility critical exponent γ/ν_{⊥}=2.0±0.2, in agreement with theoretical results for the strong fluctuation regime of wetting transitions with quenched noise.

Efficient numerical methods for the random-field Ising model: Finite-size scaling, reweighting extrapolation, and computation of response functions

It was recently shown [Phys. Rev. Lett. 110, 227201 (2013)PRLTAO0031-900710.1103/PhysRevLett.110.227201] that the critical behavior of the random-field Ising model in three dimensions is ruled by a single universality class. This conclusion was reached only after a proper taming of the large scaling corrections of the model by applying a combined approach of various techniques, coming from the zero- and positive-temperature toolboxes of statistical physics. In the present contribution we provide a detailed description of this combined scheme, explaining in detail the zero-temperature numerical scheme and developing the generalized fluctuation-dissipation formula that allowed us to compute connected and disconnected correlation functions of the model. We discuss the error evolution of our method and we illustrate the infinite limit-size extrapolation of several observables within phenomenological renormalization. We present an extension of the quotients method that allows us to obtain estimates of the critical exponent α of the specific heat of the model via the scaling of the bond energy and we discuss the self-averaging properties of the system and the algorithmic aspects of the maximum-flow algorithm used.

Phase Transition for Quenched Coupled Replicas in a Plaquette Spin Model of Glasses

We study a three-dimensional plaquette spin model whose low temperature dynamics is glassy, due to localized defects and effective kinetic constraints. The thermodynamics of this system is smooth at all temperatures. We show that coupling it to a second system with a fixed (quenched) configuration leads to a phase transition, at finite coupling. The order parameter is the overlap between the copies, and the transition is between phases of low and high overlap. We find critical points whose properties are consistent with random-field Ising universality. We analyze the interfacial free energy cost between the high- and low-overlap states that coexist at (and below) the critical point, and we use this cost as the basis for a finite-size scaling analysis. We discuss these results in the context of mean-field and dynamical facilitation theories of the glass transition.

First-order phase transition and tricritical scaling behavior of the Blume-Capel model: A Wang-Landau sampling approach

We investigate the tricritical scaling behavior of the two-dimensional spin-1 Blume-Capel model by using the Wang-Landau method of measuring the joint density of states for lattice sizes up to 48×48 sites. We find that the specific heat deep in the first-order area of the phase diagram exhibits a double-peak structure of the Schottky-like anomaly appearing with the transition peak. The first-order transition curve is systematically determined by employing the method of field mixing in conjunction with finite-size scaling, showing a significant deviation from the previous data points. At the tricritical point, we characterize the tricritical exponents through finite-size-scaling analysis including the phenomenological finite-size scaling with thermodynamic variables. Our estimation of the tricritical eigenvalue exponents, yt=1.804(5), yg=0.80(1), and yh=1.925(3), provides the first Wang-Landau verification of the conjectured exact values, demonstrating the effectiveness of the density-of-states-based approach in finite-size scaling study of multicritical phenomena.

Finite-size scaling, dynamic fluctuations, and hyperscaling relation in the Kuramoto model

We revisit the Kuramoto model to explore the finite-size scaling (FSS) of the order parameter and its dynamic fluctuations near the onset of the synchronization transition, paying particular attention to effects induced by the randomness of the intrinsic frequencies of oscillators. For a population of size N, we study two ways of sampling the intrinsic frequencies according to the same given unimodal distribution g(ω). In the "random" case, frequencies are generated independently in accordance with g(ω), which gives rise to oscillator number fluctuation within any given frequency interval. In the "regular" case, the N frequencies are generated in a deterministic manner that minimizes the oscillator number fluctuations, leading to quasiuniformly spaced frequencies in the population. We find that the two samplings yield substantially different finite-size properties with clearly distinct scaling exponents. Moreover, the hyperscaling relation between the order parameter and its fluctuations is valid in the regular case, but it is violated in the random case. In this last case, a self-consistent mean-field theory that completely ignores dynamic fluctuations correctly predicts the FSS exponent of the order parameter but not its critical amplitude.

Fluids in porous media: the case of neutral walls

The bulk phase behavior of a fluid is typically altered when the fluid is brought into confinement by the walls of a random porous medium. Inside the porous medium, phase-transition points are shifted, or may disappear altogether. A crucial determinant is how the walls interact with the fluid particles. In this work, we consider the situation whereby the walls are neutral with respect to the liquid and vapor phases. In order to realize the condition of strict neutrality, we use a symmetric binary mixture inside a porous medium that interacts identically with mixture species. Monte Carlo simulations are then used to obtain the phase behavior. Our main finding is that, in the presence of the porous medium, a liquid-vapor critical point still exists. At the critical point, the distribution of the order parameter remains scale invariant, but self-averaging is violated. These findings provide further evidence that random confinement by neutral walls induces critical behavior of the random Ising model (i.e., Ising models with dilution type disorder, where the disorder couples to the energy).

Dynamical properties of random-field Ising model

Extensive Monte Carlo simulations are performed on a two-dimensional random field Ising model. The purpose of the present work is to study the disorder-induced changes in the properties of disordered spin systems. The time evolution of the domain growth, the order parameter, and the spin-spin correlation functions are studied in the nonequilibrium regime. The dynamical evolution of the order parameter and the domain growth shows a power law scaling with disorder-dependent exponents. It is observed that for weak random fields, the two-dimensional random field Ising model possesses long-range order. Except for weak disorder, exchange interaction never wins over pinning interaction to establish long-range order in the system.

Random pinning limits the size of membrane adhesion domains

Theoretical models describing specific adhesion of membranes predict (for certain parameters) a macroscopic phase separation of bonds into adhesion domains. We show that this behavior is fundamentally altered if the membrane is pinned randomly due to, e.g., proteins that anchor the membrane to the cytoskeleton. Perturbations which locally restrict membrane height fluctuations induce quenched disorder of the random-field type. This rigorously prevents the formation of macroscopic adhesion domains following the Imry-Ma argument [Imry and Ma, Phys. Rev. Lett. 35, 1399 (1975)]. Our prediction of random-field disorder follows from analytical calculations and is strikingly confirmed in large-scale Monte Carlo simulations. These simulations are based on an efficient composite Monte Carlo move, whereby membrane height and bond degrees of freedom are updated simultaneously in a single move. The application of this move should prove rewarding for other systems also.

Membrane lateral structure: the influence of immobilized particles on domain size

In experiments on model membranes, formation of large domains of different lipid composition is readily observed. However, no such phase separation is observed in the membranes of intact cells. Instead, small transient inhomogeneities called lipid rafts are expected in these systems. One of the numerous attempts to explain small domains refers to the coupling of the membrane to its surroundings, which leads to the immobilization of some of the membrane molecules. These immobilized molecules then act as static obstacles for the remaining mobile ones. We present detailed Molecular Dynamics simulations demonstrating that this can indeed account for small domains. This confirms previous Monte Carlo studies based on simplified models. Furthermore, by directly comparing domain structures obtained using Molecular Dynamics to Monte Carlo simulations of the Ising model, we demonstrate that domain formation in the presence of obstacles is remarkably insensitive to the details of the molecular interactions.

Fluids with quenched disorder: scaling of the free energy barrier near critical points

In the context of Monte Carlo simulations, the analysis of the probability distribution P(L)(m) of the order parameter m, as obtained in simulation boxes of finite linear extension L, allows for an easy estimation of the location of the critical point and the critical exponents. For Ising-like systems without quenched disorder, P(L)(m) becomes scale-invariant at the critical point, where it assumes a characteristic bimodal shape featuring two overlapping peaks. In particular, the ratio between the value of P(L)(m) at the peaks (P(L, max)) and the value at the minimum in between (P(L, min)) becomes L-independent at criticality. However, for Ising-like systems with quenched random fields, we argue that instead ΔF(L) := ln(P(L, max)/P(L, min)) proportional to L(θ) should be observed, where θ > 0 is the 'violation of hyperscaling' exponent. Since θ is substantially non-zero, the scaling of ΔF(L) with system size should be easily detectable in simulations. For two fluid models with quenched disorder, ΔF(L) versus L was measured and the expected scaling was confirmed. This provides further evidence that fluids with quenched disorder belong to the universality class of the random field Ising model.

Universality in disordered systems: the case of the three-dimensional random-bond Ising model

We study the critical behavior of the d=3 Ising model with bond randomness through extensive Monte Carlo simulations and finite-size scaling techniques. Our results indicate that the critical behavior of the random-bond model is governed by the same universality class as the site- and bond-diluted models, clearly distinct from that of the pure model, thus providing a complete set of universality in disordered systems.