Sandpiles subjected to sinusoidal drive

This paper considers a sandpile model subjected to a sinusoidal external drive with the period T. We develop a theoretical model for the Green's function in a large T limit, which predicts that the avalanches are anisotropic and elongated in the oscillation direction. We track the problem numerically and show that the system additionally shows a regime where the avalanches are elongated in the perpendicular direction with respect to the oscillations. We find a crossover between these two regimes. The power spectrum of avalanche size and the grains wasted from the parallel and perpendicular directions are studied. These functions show power-law behavior in terms of the frequency with exponents, which run with T.

Self-similar but not conformally invariant traces obtained by modified Loewner forces

The two-dimensional Loewner exploration process is generalized to the case where the random force is self-similar with positively correlated increments. We model this random force by a fractional Brownian motion with Hurst exponent H≥1/2≡H_{BM}, where H_{BM} stands for the one-dimensional Brownian motion. By manipulating the deterministic force, we design a scale-invariant equation describing self-similar traces which lack conformal invariance. The model is investigated in terms of the "input diffusivity parameter" κ, which coincides with the one of the ordinary Schramm-Loewner evolution (SLE) at H=H_{BM}. In our numerical investigation, we focus on the scaling properties of the traces generated for κ=2,3, κ=4, and κ=6,8 as the representatives, respectively, of the dilute phase, the transition point, and the dense phase of the ordinary SLE. The resulting traces are shown to be scale invariant. Using two equivalent schemes, we extract the fractal dimension, D_{f}(H), of the traces which decrease monotonically with increasing H, reaching D_{f}=1 at H=1 for all κ values. The left passage probability (LPP) test demonstrates that, for H values not far from the uncorrelated case (small ε_{H}≡H-H_{BM}/H_{BM}), the prediction of the ordinary SLE is applicable with an effective diffusivity parameter κ_{eff}. Not surprisingly, the κ_{eff}'s do not fulfill the prediction of SLE for the relation between D_{f}(H) and the diffusivity parameter.

Invasion percolation in short-range and long-range disorder background

In the original invasion percolation model, a random number quantifies the role of necks, or generally the quality of pores, ignoring the structure of pores and impermeable regions (to which the invader cannot enter). In this paper, we investigate invasion percolation (IP), taking into account the impermeable regions, the configuration of which is modeled by ordinary and Ising-correlated site percolation (with short-range interactions, SRI), on top of which the IP dynamics is defined. We model the long-ranged correlations of pores by a random Coulomb potential (RCP). By examining various dynamical observables, we suggest that the critical exponents of Ising-correlated cases change considerably only in the vicinity of the critical point (critical temperature), while for the ordinary percolation case the exponents are robust against the occupancy parameter p. The properties of the model for the long-range interactions [LRI (RCP)] are completely different from the normal IP. In particular, the fractal dimension of the external frontier of the largest hole is nearly 4/3 for SRI far from the critical points, which is compatible with normal IP, while it converges to 1.099±0.04 for RCP. For the latter case, the time dependence of our observables is divided into three parts: the power law (short time), the logarithmic (moderate time), and the linear (long time) regimes. The second crossover time is shown to go to infinity in the thermodynamic limit, whereas the first crossover time is nearly unchanged, signaling the dominance of the logarithmic regime. The average gyration radius of the growing clusters, the length of their external perimeter, and the corresponding roughness are shown to be nearly constant for the long-time regime in the thermodynamic limit.

Self-organized criticality in cumulus clouds

The shape of clouds has proven to be essential for classifying them. Our analysis of images from fair weather cumulus clouds reveals that, in addition to turbulence, they are driven by self-organized criticality. Our observations yield exponents that support the fact the clouds, when projected to two dimensions, exhibit conformal symmetry compatible with c=-2 conformal field theory. By using a combination of the Navier-Stokes equation, diffusion equations, and a coupled map lattice, we successfully simulated cloud formation, and obtained the same exponents.

Superstatistical two-temperature Ising model

The previous approach of the nonequilibrium Ising model was based on the local temperature in which each site or part of the system has its own specific temperature. We introduce an approach of the two-temperature Ising model as a prototype of the superstatistic critical phenomena. The model is described by two temperatures (T_{1},T_{2}) in a zero magnetic field. To predict the phase diagram and numerically estimate the exponents, we develop the Metropolis and Swendsen-Wang Monte Carlo method. We observe that there is a nontrivial critical line, separating ordered and disordered phases. We propose an analytic equation for the critical line in the phase diagram. Our numerical estimation of the critical exponents illustrates that all points on the critical line belong to the ordinary Ising universality class.

Geometry-induced nonequilibrium phase transition in sandpiles

We study the sandpile model on three-dimensional spanning Ising clusters with the temperature T treated as the control parameter. By analyzing the three-dimensional avalanches and their two-dimensional projections (which show scale-invariant behavior for all temperatures), we uncover two universality classes with different exponents (an ordinary BTW class, and SOC_{T=∞}), along with a tricritical point (at T_{c}, the critical temperature of the host) between them. The transition between these two criticalities is induced by the transition in the support. The SOC_{T=∞} universality class is characterized by the exponent of the avalanche size distribution τ^{T=∞}=1.27±0.03, consistent with the exponent of the size distribution of the Barkhausen avalanches in amorphous ferromagnets Durin and Zapperi [Phys. Rev. Lett. 84, 4705 (2000)PRLTAO0031-900710.1103/PhysRevLett.84.4705]. The tricritical point is characterized by its own critical exponents. In addition to the avalanche exponents, some other quantities like the average height, the spanning avalanche probability (SAP), and the average coordination number of the Ising clusters change significantly the behavior at this point, and also exhibit power-law behavior in terms of ε≡T-T_{c}/T_{c}, defining further critical exponents. Importantly, the finite-size analysis for the activity (number of topplings) per site shows the scaling behavior with exponents β=0.19±0.02 and ν=0.75±0.05. A similar behavior is also seen for the SAP and the average avalanche height. The fractal dimension of the external perimeter of the two-dimensional projections of avalanches is shown to be robust against T with the numerical value D_{f}=1.25±0.01.

Elastic backbone phase transition in the Ising model

The two-dimensional (zero magnetic field) Ising model is known to undergo a second-order paraferromagnetic phase transition, which is accompanied by a correlated percolation transition for the Fortuin-Kasteleyn (FK) clusters. In this paper we uncover that there exists also a second temperature T_{eb}<T_{c} at which the elastic backbone of FK clusters undergoes a second-order phase transition to a dense phase. The corresponding universality class, which is characterized by determining various percolation exponents, is shown to be completely different from directed percolation, which leads us to propose a new anisotropic universality class with β=0.54±0.02, ν_{||}=1.86±0.01, ν_{⊥}=1.21±0.04, and d_{f}=1.53±0.03. All tested hyperscaling relations are shown to be valid.

Scaling properties of monolayer graphene away from the Dirac point

The statistical properties of the carrier density profile of graphene in the ground state in the presence of particle-particle interaction and random charged impurity in zero gate voltage has been recently obtained by Najafi et al. [Phys. Rev. E 95, 032112 (2017)2470-004510.1103/PhysRevE.95.032112]. The nonzero chemical potential (μ) in gated graphene has nontrivial effects on electron-hole puddles, since it generates mass in the Dirac action and destroys the scaling behaviors of the effective Thomas-Fermi-Dirac theory. We provide detailed analysis on the resulting spatially inhomogeneous system in the framework of the Thomas-Fermi-Dirac theory for the Gaussian (white noise) disorder potential. We show that the chemical potential in this system as a random surface destroys the self-similarity, and also the charge field is non-Gaussian. We find that the two-body correlation functions are factorized to two terms: a pure function of the chemical potential and a pure function of the distance. The spatial dependence of these correlation functions is double logarithmic, e.g., the two-point density correlation behaves like D_{2}(r,μ)∝μ^{2}exp[-(-a_{D}lnlnr^{β_{D}})^{α_{D}}] (α_{D}=1.82, β_{D}=0.263, and a_{D}=0.955). The Fourier power spectrum function also behaves like ln[S(q)]=-β_{S}^{-a_{S}}(lnq)^{a_{S}}+2lnμ (a_{S}=3.0±0.1 and β_{S}=2.08±0.03) in contrast to the ordinary Gaussian rough surfaces for which a_{S}=1 and β_{S}=1/2(1+α)^{-1} (α being the roughness exponent). The geometrical properties are, however, similar to the ungated (μ=0) case, with the exponents that are reported in the text.

Self-avoiding walk on a square lattice with correlated vacancies

The self-avoiding walk on the square site-diluted correlated percolation lattice is considered. The Ising model is employed to realize the spatial correlations of the metric space. As a well-accepted result, the (generalized) Flory's mean-field relation is tested to measure the effect of correlation. After exploring a perturbative Fokker-Planck-like equation, we apply an enriched Rosenbluth Monte Carlo method to study the problem. To be more precise, the winding angle analysis is also performed from which the diffusivity parameter of Schramm-Loewner evolution theory (κ) is extracted. We find that at the critical Ising (host) system, the exponents are in agreement with Flory's approximation. For the off-critical Ising system, we find also a behavior for the fractal dimension of the walker trace in terms of the correlation length of the Ising system ξ(T), i.e., D_{F}^{SAW}(T)-D_{F}^{SAW}(T_{c})∼1/sqrt[ξ(T)].

Mapping of the Bak, Tang, and Wiesenfeld sandpile model on a two-dimensional Ising-correlated percolation lattice to the two-dimensional self-avoiding random walk

The self-organized criticality on the random fractal networks has many motivations, like the movement pattern of fluid in the porous media. In addition to the randomness, introducing correlation between the neighboring portions of the porous media has some nontrivial effects. In this paper, we consider the Ising-like interactions between the active sites as the simplest method to bring correlations in the porous media, and we investigate the statistics of the BTW model in it. These correlations are controlled by the artificial "temperature" T and the sign of the Ising coupling. Based on our numerical results, we propose that at the Ising critical temperature T_{c} the model is compatible with the universality class of two-dimensional (2D) self-avoiding walk (SAW). Especially the fractal dimension of the loops, which are defined as the external frontier of the avalanches, is very close to D_{f}^{SAW}=4/3. Also, the corresponding open curves has conformal invariance with the root-mean-square distance R_{rms}∼t^{3/4} (t being the parametrization of the curve) in accordance with the 2D SAW. In the finite-size study, we observe that at T=T_{c} the model has some aspects compatible with the 2D BTW model (e.g., the 1/log(L)-dependence of the exponents of the distribution functions) and some in accordance with the Ising model (e.g., the 1/L-dependence of the fractal dimensions). The finite-size scaling theory is tested and shown to be fulfilled for all statistical observables in T=T_{c}. In the off-critical temperatures in the close vicinity of T_{c} the exponents show some additional power-law behaviors in terms of T-T_{c} with some exponents that are reported in the text. The spanning cluster probability at the critical temperature also scales with L^{1/2}, which is different from the regular 2D BTW model.