Restoration of dimensional reduction in the random-field Ising model at five dimensions

Fluids in Random Media and Dimensional Augmentation

We propose a solution to the puzzle of dimensional reduction in the random field Ising model, asking the following: To what random problem in D=d+2 dimensions does a pure system in d dimensions correspond? For a continuum binary fluid and an Ising lattice gas, we prove that the mean density and other observables equal those of a similar model in D dimensions, but with infinite range interactions and correlated disorder in the extra two dimensions. There is no conflict with rigorous results that the finite range model orders in D=3. Our arguments avoid the use of replicas and perturbative field theory, being based on convergent cluster expansions, which, for the lattice gas, may be extended to the critical point by the Lee-Yang theorem. Although our results may be derived using supersymmetry, they follow more directly from the matrix-tree theorem.

Finite-size scaling of the random-field Ising model above the upper critical dimension

Finite-size scaling above the upper critical dimension is a long-standing puzzle in the field of statistical physics. Even for pure systems various scaling theories have been suggested, partially corroborated by numerical simulations. In the present manuscript we address this problem in the even more complicated case of disordered systems. In particular, we investigate the scaling behavior of the random-field Ising model at dimension D=7, i.e., above its upper critical dimension D_{u}=6, by employing extensive ground-state numerical simulations. Our results confirm the hypothesis that at dimensions D>D_{u}, linear length scale L should be replaced in finite-size scaling expressions by the effective scale L_{eff}=L^{D/D_{u}}. Via a fitted version of the quotients method that takes this modification, but also subleading scaling corrections into account, we compute the critical point of the transition for Gaussian random fields and provide estimates for the full set of critical exponents. Thus, our analysis indicates that this modified version of finite-size scaling is successful also in the context of the random-field problem.

Random-bond antiferromagnetic Ising model in a field

Using combinatorial optimization techniques we study the critical properties of the two- and three-dimensional Ising models with uniformly distributed random antiferromagnetic couplings (1≤J_{i}≤2) in the presence of a homogeneous longitudinal field, h, at zero temperature. In finite systems of linear size, L, we measure the average correlation function, C_{L}(ℓ,h), when the sites are either on the same sublattice, or they belong to different sublattices. The phase transition, which is of first order in the pure system, turns to mixed order in two dimensions with critical exponents 1/ν≈0.5 and η≈0.7. In three dimensions we obtain 1/ν≈0.7, which is compatible with the value of the random-field Ising model, but we cannot discriminate between second-order and mixed-order transitions.

Theory and experiments for disordered elastic manifolds, depinning, avalanches, and sandpiles

Domain walls in magnets, vortex lattices in superconductors, contact lines at depinning, and many other systems can be modeled as an elastic system subject to quenched disorder. The ensuing field theory possesses a well-controlled perturbative expansion around its upper critical dimension. Contrary to standard field theory, the renormalization group (RG) flow involves a function, the disorder correlator Δ(w), and is therefore termed the functional RG. Δ(w) is a physical observable, the auto-correlation function of the center of mass of the elastic manifold. In this review, we give a pedagogical introduction into its phenomenology and techniques. This allows us to treat both equilibrium (statics), and depinning (dynamics). Building on these techniques, avalanche observables are accessible: distributions of size, duration, and velocity, as well as the spatial and temporal shape. Various equivalences between disordered elastic manifolds, and sandpile models exist: an elastic string driven at a point and the Oslo model; disordered elastic manifolds and Manna sandpiles; charge density waves and Abelian sandpiles or loop-erased random walks. Each of the mappings between these systems requires specific techniques, which we develop, including modeling of discrete stochastic systems via coarse-grained stochastic equations of motion, super-symmetry techniques, and cellular automata. Stronger than quadratic nearest-neighbor interactions lead to directed percolation, and non-linear surface growth with additional Kardar-Parisi-Zhang (KPZ) terms. On the other hand, KPZ without disorder can be mapped back to disordered elastic manifolds, either on the directed polymer for its steady state, or a single particle for its decay. Other topics covered are the relation between functional RG and replica symmetry breaking, and random-field magnets. Emphasis is given to numerical and experimental tests of the theory.

Parisi-Sourlas Supersymmetry in Random Field Models

By the Parisi-Sourlas conjecture, the critical point of a theory with random field (RF) disorder is described by a supersymmeric (SUSY) conformal field theory (CFT), related to a d-2 dimensional CFT without SUSY. Numerical studies indicate that this is true for the RF ϕ^{3} model but not for the RF ϕ^{4} model in d<5 dimensions. Here we argue that the SUSY fixed point is not reached because of new relevant SUSY-breaking interactions. We use a perturbative renormalization group in a judiciously chosen field basis, allowing systematic exploration of the space of interactions. Our computations agree with the numerical results for both cubic and quartic potential.

Universality in the two-dimensional dilute Baxter-Wu model

We study the question of universality in the two-dimensional spin-1 Baxter-Wu model in the presence of a crystal field Δ. We employ extensive numerical simulations of two types, providing us with complementary results: Wang-Landau sampling at fixed values of Δ and a parallelized variant of the multicanonical approach performed at constant temperature T. A detailed finite-size scaling analysis in the regime of second-order phase transitions in the (Δ,T) phase diagram indicates that the transition belongs to the universality class of the four-state Potts model. Previous controversies with respect to the nature of the transition are discussed and attributed to the presence of strong finite-size effects, especially as one approaches the pentacritical point of the model.

Mechanism of subcritical avalanche propagation in three-dimensional disordered systems

We present a numerical study on necessary conditions for the appearance of infinite avalanche below the critical point in disordered systems that evolve throughout metastable states. The representative of those systems is the nonequilibrium athermal random-field Ising model. We investigate the impact on propagation of infinite avalanche of both the interface of flipped spins at the avalanche's starting point and the number of independent islands of flipped spins in the system at the moment when the avalanche starts. To deduce what effects are originated due to finite system's size, and to distinguish them from the real necessary conditions for the appearance of the infinite avalanche, we examined lattices of different sizes as well as other key parameters for the avalanche propagation.

Effects of external noise on threshold-induced correlations in ferromagnetic systems

In the present paper we investigate the impact of the external noise and detection threshold level on the simulation data for the systems that evolve through metastable states. As a representative model of such systems we chose the nonequilibrium athermal random-field Ising model with two types of the external noise, uniform white noise and Gaussian white noise with various different standard deviations, imposed on the original response signal obtained in model simulations. We applied a wide range of detection threshold levels in analysis of the signal and show how these quantities affect the values of exponent γ_{S/T} (describing the scaling of the average avalanche size with duration), the shift of waiting time between the avalanches, and finally the collapses of the waiting time distributions. The results are obtained via extensive numerical simulations on the equilateral three-dimensional cubic lattices of various sizes and disorders.

Nonequilibrium athermal random-field Ising model on hexagonal lattices

We present the results of a study providing numerical evidence for the absence of critical behavior of the nonequilibrium athermal random-field Ising model in adiabatic regime on the hexagonal two-dimensional lattice. The results are obtained on the systems containing up to 32768×32768 spins and are the averages of up to 1700 runs with different random-field configurations per each value of disorder. We analyzed regular systems as well as the systems with different preset conditions to capture behavior in thermodynamic limit. The superficial insight to the avalanche propagation in this type of lattice is given as a stimulus for further research on the topic of avalanche evolution. With obtained data we may conclude that there is no critical behavior of random-field Ising model on hexagonal lattice which is a result that differs from the ones found for the square and for the triangular lattices supporting the recent conjecture that the number of nearest neighbors affects the model criticality.

Ising universality in the two-dimensional Blume-Capel model with quenched random crystal field

Using high-precision Monte Carlo simulations based on a parallel version of the Wang-Landau algorithm and finite-size scaling techniques, we study the effect of quenched disorder in the crystal-field coupling of the Blume-Capel model on a square lattice. We mainly focus on the part of the phase diagram where the pure model undergoes a continuous transition, known to fall into the universality class of a pure Ising ferromagnet. A dedicated scaling analysis reveals concrete evidence in favor of the strong universality hypothesis with the presence of additional logarithmic corrections in the scaling of the specific heat. Our results are in agreement with an early real-space renormalization-group study of the model as well as a very recent numerical work where quenched randomness was introduced in the energy exchange coupling. Finally, by properly fine tuning the control parameters of the randomness distribution we also qualitatively investigate the part of the phase diagram where the pure model undergoes a first-order phase transition. For this region, preliminary evidence indicate a smoothing of the transition to second-order with the presence of strong scaling corrections.

Dimensional reduction breakdown and correction to scaling in the random-field Ising model

We provide a theoretical analysis by means of the nonperturbative functional renormalization group (NP-FRG) of the corrections to scaling in the critical behavior of the random-field Ising model (RFIM) near the dimension d_{DR}≈5.1 that separates a region where the renormalized theory at the fixed point is supersymmetric and critical scaling satisfies the d→d-2 dimensional reduction property (d>d_{DR}) from a region where both supersymmetry and dimensional reduction break down at criticality (d<d_{DR}). We show that the NP-FRG results are in very good agreement with recent large-scale lattice simulations of the RFIM in d=5 and we detail the consequences for the leading correction-to-scaling exponent of the peculiar boundary-layer mechanism by which the dimensional-reduction fixed point disappears and the dimensional-reduction-broken fixed point emerges in d_{DR}.

Avalanche properties in striplike ferromagnetic systems

We present numerical findings on the behavior of the athermal nonequilibrium random-field Ising model of spins at the thin striplike L_{1}×L_{2}×L_{3} cubic lattices with L_{1}<L_{2}<L_{3}. Changing of system sizes highly influences the evolution and shape of avalanches. The smallest avalanches [classified as three-dimension- (3D) like] are unaffected by the system boundaries, the larger are sandwiched between the top and bottom system faces so are 2D-like, while the largest are extended over the system lateral cross section and propagate along the length L_{3} like in 1D systems. Such a structure of avalanches causes double power-law distributions of their size, duration, and energy with larger effective critical exponent corresponding to 3D-like and smaller to 2D-like avalanches. The distributions scale with thickness L_{1} and are collapsible following the proposed scaling predictions which, together with the distributions' shape, might be important for analysis of the Barkhausen noise experimental data for striplike samples. Finally, the impact of system size on external field that triggers the largest avalanche for a given disorder is presented and discussed.

Behavior of the random-field XY model on simple cubic lattices at h_{r}=1.5

We have performed studies of the three-dimensional random-field XY model on 32 samples of L×L×L simple cubic lattices with periodic boundary conditions, with a random field strength of h_{r} = 1.5, for L= 128, using a parallelized Monte Carlo algorithm. We present results for the sample-averaged magnetic structure factor S(k[over ⃗]) over a range of temperature, using both random hot start and ferromagnetic cold start initial states, and k[over ⃗] along the [1,0,0] and [1,1,1] directions. At T= 1.875, S(k[over ⃗]) shows a broad peak near |k[over ⃗]|=0, with a correlation length which is limited by thermal fluctuations, rather than the lattice size. As T is lowered, this peak grows and sharpens. By T= 1.5, it is clear that the correlation length is larger than L= 128. The lowest temperature for which S(k[over ⃗]) was calculated is T= 1.421875, where the hot start and cold start initial conditions usually do not find the same local minimum in the phase space. Our results are consistent with the idea that there is a finite value of T below which S(k[over ⃗]) diverges slowly as |k[over ⃗]| goes to zero. This divergence would imply that the relaxation time of the spins is also diverging. That is the signature of an ergodicity-breaking phase transition.

Directed Polymers and Interfaces in Disordered Media

We consider field theory formulation for directed polymers and interfaces in the presence of quenched disorder. We write a series representation for the averaged free energy, where all the integer moments of the partition function of the model contribute. The structure of field space is analysed for polymers and interfaces at finite temperature using the saddle-point equations derived from each integer moments of the partition function. For the case of an interface we obtain the wandering exponent ξ = ( 4 - d ) / 2 , also obtained by the conventional replica method for the replica symmetric scenario.

Critical disorder and critical magnetic field of the nonequilibrium athermal random-field Ising model in thin systems

In the present study of the nonequilibrium athermal random-field Ising model we focus on the behavior of the critical disorder R_{c}(l) and the critical magnetic field H_{c}(l) under different boundary conditions when the system thickness l varies. We propose expressions for R_{c}(l) and H_{c}(l) as well as for the effective critical disorder R_{c}^{eff}(l,L) and effective critical magnetic field H_{c}^{eff}(l,L) playing the role of the effective critical parameters for the L×L×l lattices of finite lateral size L. We support these expressions by the scaling collapses of the magnetization and susceptibility curves obtained in extensive simulations. The collapses are achieved with the two-dimensional (2D) exponents for l below some characteristic value, providing thus a numerical evidence that the thin systems exhibit a 2D-like criticality which should be relevant for the experimental analyses of thin ferromagnetic samples.

Evidence for Supersymmetry in the Random-Field Ising Model at D=5

We provide a nontrivial test of supersymmetry in the random-field Ising model at five spatial dimensions, by means of extensive zero-temperature numerical simulations. Indeed, supersymmetry relates correlation functions in a D-dimensional disordered system with some other correlation functions in a D-2 clean system. We first show how to check these relationships in a finite-size scaling calculation and then perform a high-accuracy test. While the supersymmetric predictions are satisfied even to our high accuracy at D=5, they fail to describe our results at D=4.

Universality from disorder in the random-bond Blume-Capel model

Using high-precision Monte Carlo simulations and finite-size scaling we study the effect of quenched disorder in the exchange couplings on the Blume-Capel model on the square lattice. The first-order transition for large crystal-field coupling is softened to become continuous, with a divergent correlation length. An analysis of the scaling of the correlation length as well as the susceptibility and specific heat reveals that it belongs to the universality class of the Ising model with additional logarithmic corrections which is also observed for the Ising model itself if coupled to weak disorder. While the leading scaling behavior of the disordered system is therefore identical between the second-order and first-order segments of the phase diagram of the pure model, the finite-size scaling in the ex-first-order regime is affected by strong transient effects with a crossover length scale L^{*}≈32 for the chosen parameters.