Avalanche distributions in the two-dimensional nonequilibrium zero-temperature random field Ising model

Hysteresis-loop phenomena in disordered ferromagnets with demagnetizing field and finite temperature

Using numerical simulations, we investigate the impact of the demagnetization field and finite temperature on the hysteresis phenomena in disordered ferromagnetics systems. We model the behavior of thin systems employing the thermal nonequilibrium random field Ising model driven by a finite-driving rate protocol to study the shape of the hysteresis loop and demagnetization line and the magnetization fluctuations for varied parameters. Our results reveal a significant interplay of the disorder, the demagnetizing fields, and thermal fluctuations. In particular, at a fixed disorder and temperature, the increasing demagnetization coefficient gradually prologues the magnetization reversal process, altering the multifractal nature of the magnetization fluctuations. The process accompanies the appearance of extended linear segments in the hysteresis loop and changes in the demagnetization line while practically preserving the value of the coercive field and slightly changing the remanent magnetization. On the other hand, increasing temperature expands the system's response fluctuations and narrows the loop, affecting the coercive field and remanent magnetization. The interplay of thermal fluctuations and demagnetizing fields fully manifests in limiting the large-scale magnetization fluctuations, revealed by the multifractal spectra and the scaling functions of the avalanches. Our research offers new modeling perspectives with a more realistic scenario and provides insight into new hysteresis loop phenomena relevant to theoretical analysis and applications.

Barkhausen noise in disordered striplike ferromagnets: Experiment versus simulations

In this work, we present a systematic comparison of the results obtained from the low-frequency Barkhausen noise recordings in nanocrystalline samples with those from the numerical simulations of the random-field Ising model systems. We performed measurements at room temperature on a field-driven metallic glass stripe made of VITROPERM 800 R, a nanocrystalline iron-based material with an excellent combination of soft and magnetic properties, making it a cutting-edge material for a wide range of applications. Given that the Barkhausen noise emissions emerging along a hysteresis curve are stochastic and depend in general on a variety of factors (such as distribution of disorder due to impurities or defects, varied size of crystal grains, type of domain structure, driving rate of the external magnetic field, sample shape and temperature, etc.), adequate theoretical modeling is essential for their interpretation and prediction. Here the Random field Ising model, specifically its athermal nonequilibrium version with the finite driving rate, stands out as an appropriate choice due to the material's nanocrystalline structure and high Curie temperature. We performed a systematic analysis of the signal properties and magnetization avalanches comparing the outcomes of the numerical model and experiments carried out in a two-decade-wide range of the external magnetic field driving rates. Our results reveal that with a suitable choice of parameters, a considerable match with the experimental results is achieved, indicating that this model can accurately describe the Barkhausen noise features in nanocrystalline samples.

Interplay of disorder and type of driving in disordered ferromagnetic systems

We investigate the effects of adiabatic, quasistatic, and finite-rate types of driving on the evolution of disordered three-dimensional ferromagnetic systems, studied within the frame of the nonequilibrium athermal random field Ising model. The effects were examined in all three domains of disorder (low, high, and transitional) for all types of driving, and in a wide range of driving rates for quasistatic and finite-rate driving, providing an extensive overview and comparison of the joint effects that the disorder, type of driving, and rate regime have on the system's behavior.

Linking space-time correlations for a class of self-organized critical systems

The hypothesis of self-organized criticality explains the existence of long-range "space-time" correlations, observed inseparably in many natural dynamical systems. A simple link between these correlations is yet unclear, particularly in fluctuations at an "external drive" timescale. As an example, we consider a class of sandpile models displaying nontrivial correlations. We apply the scaling method and determine spatial cross-correlation by establishing a relationship between local and global temporal correlations. We find that the spatial cross-correlation decays in a power-law manner with an exponent γ=1-δ, where δ characterizes a scaling of the total power of the global temporal process with the system size.

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.

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.

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.

Identifying epidemic threshold by temporal profile of outbreaks on networks

Identifying epidemic threshold is of great significance in preventing and controlling disease spreading on real-world networks. Previous studies have proposed different theoretical and numerical approaches to determine the epidemic threshold for the susceptible-infected-recovered (SIR) model, but the numerical study of the critical points on networks by utilizing temporal characteristics of epidemic outbreaks is still lacking. Here, we study the temporal profile of epidemic outbreaks, i.e., the average avalanche shapes of a fixed duration. At the critical point, the rescaled average terminating and nonterminating avalanche shapes for different durations collapse onto two universal curves, respectively, while the average number of subsequent events essentially remains constant. We propose two numerical measures to determine the epidemic threshold by analyzing the convergence of the rescaled average nonterminating avalanche shapes for varying durations and the stability of the average number of subsequent events, respectively. Extensive numerical simulations demonstrate that our methods can accurately identify the numerical threshold for the SIR dynamics on synthetic and empirical networks. Compared with traditional numerical measures, our methods are more efficient due to the constriction of observation duration and thus are more applicable to large-scale networks. This work helps one to understand the temporal profile of disease propagation and would promote further studies on the phase transition of epidemic dynamics.

The critical Barkhausen avalanches in thin random-field ferromagnets with an open boundary

The interplay between the critical fluctuations and the sample geometry is investigated numerically using thin random-field ferromagnets exhibiting the field-driven magnetisation reversal on the hysteresis loop. The system is studied along the theoretical critical line in the plane of random-field disorder and thickness. The thickness is varied to consider samples of various geometry between a two-dimensional plane and a complete three-dimensional lattice with an open boundary in the direction of the growing thickness. We perform a multi-fractal analysis of the Barkhausen noise signals and scaling of the critical avalanches of the domain wall motion. Our results reveal that, for sufficiently small thickness, the sample geometry profoundly affects the dynamics by modifying the spectral segments that represent small fluctuations and promoting the time-scale dependent multi-fractality. Meanwhile, the avalanche distributions display two distinct power-law regions, in contrast to those in the two-dimensional limit, and the average avalanche shapes are asymmetric. With increasing thickness, the scaling characteristics and the multi-fractal spectrum in thicker samples gradually approach the hysteresis loop criticality in three-dimensional systems. Thin ferromagnetic films are growing in importance technologically, and our results illustrate some new features of the domain wall dynamics induced by magnetisation reversal in these systems.

Barkhausen Noise from Precessional Domain Wall Motion

The jerky dynamics of domain walls driven by applied magnetic fields in disordered ferromagnets-the Barkhausen effect-is a paradigmatic example of crackling noise. We study Barkhausen noise in disordered Pt/Co/Pt thin films due to precessional motion of domain walls using full micromagnetic simulations, allowing for a detailed description of the domain wall internal structure. In this regime the domain walls contain topological defects known as Bloch lines which repeatedly nucleate, propagate, and annihilate within the domain wall during the Barkhausen jumps. In addition to bursts of domain wall propagation, the in-plane Bloch line dynamics within the domain wall exhibits crackling noise and constitutes the majority of the overall spin rotation activity.

Hysteresis in the Ising model with Glauber dynamics

We use Glauber dynamics to study time and temperature dependence of hysteresis in the pure (without quenched disorder) Ising model on cubic, square, honeycomb lattices as well as random graphs. Results are discussed in the context of more extensive studies of hysteresis in the random field Ising model.

Universality of Critically Pinned Interfaces in Two-Dimensional Isotropic Random Media

Based on extensive simulations, we conjecture that critically pinned interfaces in two-dimensional isotropic random media with short-range correlations are always in the universality class of ordinary percolation. Thus, in contrast to interfaces in >2 dimensions, there is no distinction between fractal (i.e., percolative) and rough but nonfractal interfaces. Our claim includes interfaces in zero-temperature random field Ising models (both with and without spontaneous nucleation), in heterogeneous bootstrap percolation, and in susceptible-weakened-infected-removed epidemics. It does not include models with long-range correlations in the randomness and models where overhangs are explicitly forbidden (which would imply nonisotropy of the medium).

Crossover from three-dimensional to two-dimensional systems in the nonequilibrium zero-temperature random-field Ising model

We present extensive numerical studies of the crossover from three-dimensional to two-dimensional systems in the nonequilibrium zero-temperature random-field Ising model with metastable dynamics. Bivariate finite-size scaling hypotheses are presented for systems with sizes L×L×l which explain the size-driven critical crossover from two dimensions (l=const, L→∞) to three dimensions (l∝L→∞). A model of effective critical disorder R_{c}^{eff}(l,L) with a unique fitting parameter and no free parameters in the R_{c}^{eff}(l,L→∞) limit is proposed, together with expressions for the scaling of avalanche distributions bringing important implications for related experimental data analysis, especially in the case of thin three-dimensional systems.

Threshold-induced correlations in the Random Field Ising Model

We present a numerical study of the correlations in the occurrence times of consecutive crackling noise events in the nonequilibrium zero-temperature Random Field Ising model in three dimensions. The critical behavior of the system is portrayed by the intermittent bursts of activity known as avalanches with scale-invariant properties which are power-law distributed. Our findings, based on the scaling analysis and collapse of data collected in extensive simulations show that the observed correlations emerge upon applying a finite threshold to the pertaining signals when defining events of interest. Such events are called subavalanches and are obtained by separation of original avalanches in the thresholding process. The correlations are evidenced by power law distributed waiting times and are present in the system even when the original avalanche triggerings are described by a random uncorrelated process.

Critical behavior of the two-dimensional nonequilibrium zero-temperature random field Ising model on a triangular lattice

We present a numerical study of the critical behavior of the nonequilibrium zero-temperature random field Ising model in two dimensions on a triangular lattice. Our findings, based on the scaling analysis and collapse of data collected in extensive simulations of systems with linear sizes up to L=65536, show that the model is in a different universality class than the same model on a quadratic lattice, which is relevant for a better understanding of model universality and the analysis of experimental data.

Criteria for infinite avalanches in the zero-temperature nonequilibrium random-field Ising model on a Bethe lattice

We present general criteria for the occurrence of infinite avalanches and critical hysteresis in the zero-temperature nonequilibrium random-field Ising model on a Bethe lattice. Drawing upon extant results as well as a result on a dilute four-coordinated (z=4) lattice, we show that diverging avalanches can occur if an arbitrarily small fraction of sites on a spanning cluster have connectivity z≥4.

Influence of the aspect ratio and boundary conditions on universal finite-size scaling functions in the athermal metastable two-dimensional random field Ising model

This work studies universal finite size scaling functions for the number of one-dimensional spanning avalanches in a two-dimensional (2D) disordered system with boundary conditions of different nature and different aspect ratios. To this end, we will consider the 2D random field Ising model at T=0 driven by the external field H with athermal dynamics implemented with periodic and forced boundary conditions. We have chosen a convenient scaling variable z that accounts for the deformation of the distance to the critical point caused by the aspect ratio. In addition, assuming that the dependence of the finite size scaling functions on the aspect ratio can be accounted for by an additional multiplicative factor, we have been able to collapse data for different system sizes, different aspect ratios, and different types of the boundary conditions into a single scaling function Q̂.

Nonequilibrium random-field Ising model on a diluted triangular lattice

We study critical hysteresis in the random-field Ising model on a two-dimensional periodic lattice with a variable coordination number z(eff) in the range 3≤z(eff)≤6. We find that the model supports critical behavior in the range 4<z(eff)≤6, but the critical exponents are independent of z(eff). The result is discussed in the context of the universality of nonequilibrium critical phenomena and extant results in the field.

Analysis of spanning avalanches in the two-dimensional nonequilibrium zero-temperature random-field Ising model

We present a numerical analysis of spanning avalanches in a two-dimensional (2D) nonequilibrium zero-temperature random field Ising model. Finite-size scaling analysis, performed for distribution of the average number of spanning avalanches per single run, spanning avalanche size distribution, average size of spanning avalanche, and contribution of spanning avalanches to magnetization jump, is augmented by analysis of spanning field (i.e., field triggering spanning avalanche), which enabled us to collapse averaged magnetization curves below critical disorder. Our study, based on extensive simulations of sufficiently large systems, reveals the dominant role of subcritical 2D-spanning avalanches in model behavior below and at the critical disorder. Other types of avalanches influence finite systems, but their contribution for large systems remains small or vanish.

Effect of coordination number on the nonequilibrium critical point

We study the nonequilibrium critical point of the zero-temperature random-field Ising model on a triangular lattice and compare it with known results on honeycomb, square, and simple cubic lattices. We suggest that the coordination number of the lattice rather than its dimension plays the key role in determining the universality class of the nonequilibrium critical behavior. This is discussed in the context of numerical evidence that equilibrium and nonequilibrium critical points of the zero-temperature random-field Ising model belong to the same universality class. The physics of this curious result is not fully understood.

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.

Scaling behavior in probabilistic neuronal cellular automata

We study a neural network model of interacting stochastic discrete two-state cellular automata on a regular lattice. The system is externally tuned to a critical point which varies with the degree of stochasticity (or the effective temperature). There are avalanches of neuronal activity, namely, spatially and temporally contiguous sites of activity; a detailed numerical study of these activity avalanches is presented, and single, joint, and marginal probability distributions are computed. At the critical point, we find that the scaling exponents for the variables are in good agreement with a mean-field theory.