Self-organized critical state of sandpile automaton models

Average height for Abelian sandpiles and the looping constant on Sierpiński graphs

For the Abelian sandpile model on Sierpiński graphs, we investigate several statistics such as average height, height probabilities and looping constant. In particular, we calculate the expected average height of a recurrent sandpile on the finite iterations of the Sierpiński gasket and we also give an algorithmic approach for calculating the height probabilities of recurrent sandpiles under stationarity by using the connection between recurrent configurations of the Abelian sandpile Markov chain and uniform spanning trees. We also calculate the expected fraction of vertices of height i for i ∈ { 0 , 1 , 2 , 3 } of sandpiles under stationarity and relate the bulk average height to the looping constant on the Sierpiński gasket.

Anomalous relaxation and hyperuniform fluctuations in center-of-mass conserving systems with broken time-reversal symmetry

We study the Oslo model, a paradigm for absorbing-phase transition, on a one-dimensional ring of L sites with a fixed global density ρ[over ¯]; we consider the system strictly above critical density ρ_{c}. Notably, microscopic dynamics conserve both mass and center of mass (CoM), but lack time-reversal symmetry. We show that, despite having highly constrained dynamics due to CoM conservation, the system exhibits diffusive relaxation away from criticality and superdiffusive relaxation near criticality. Furthermore, the CoM conservation severely restricts particle movement, causing the mobility-a transport coefficient analogous to the conductivity for charged particles-to vanish exactly. Indeed, the steady-state temporal growth of current fluctuation is qualitatively different from that observed in diffusive systems with a single conservation law. Remarkably, far from criticality where the relative density Δ=ρ[over ¯]-ρ_{c}≫ρ_{c}, the second cumulant, or the variance, 〈Q_{i}^{2}(T,Δ)〉_{c}, of current Q_{i} across the ith bond up to time T in the steady-state saturates as 〈Q_{i}^{2}〉_{c}≃Σ_{Q}^{2}(Δ)-constT^{-1/2}; near criticality, it grows subdiffusively as 〈Q_{i}^{2}〉_{c}∼T^{α}, with 0<α<1/2, and eventually saturates to Σ_{Q}^{2}(Δ). Interestingly, the asymptotic current fluctuation Σ_{Q}^{2}(Δ) is a nonmonotonic function of Δ: It diverges as Σ_{Q}^{2}(Δ)∼Δ^{2} for Δ≫ρ_{c} and Σ_{Q}^{2}(Δ)∼Δ^{-δ}, with δ>0, for Δ→0^{+}. Using a mass-conservation principle, we exactly determine the exponents δ=2(1-1/ν_{⊥})/ν_{⊥} and α=δ/zν_{⊥} via the correlation-length and dynamic exponents, ν_{⊥} and z, respectively. Finally, we show that in the steady state the self-diffusion coefficient D_{s}(ρ[over ¯]) of tagged particles is connected to activity through the relation D_{s}(ρ[over ¯])=a(ρ[over ¯])/ρ[over ¯].

Self-organized criticality of magnetic avalanches in disordered ferrimagnetic material

We observe multiple steplike jumps in a Dy-Fe-Ga-based ferrimagnetic alloy in its magnetic hysteresis curve at 2 K. The observed jumps are found to have a stochastic character with respect to their magnitude and the field position, and the jumps do not correlate with the duration of the field. The distribution of jump size follows a power law variation indicating the scale invariance nature of the jumps. We have invoked a simple two-dimensional random bond Ising-type spin system to model the dynamics. Our computational model can qualitatively reproduce the jumps and their scale-invariant character. It also elucidates that the flipping of antiferromagnetically coupled Dy and Fe clusters is responsible for the observed jumps in the hysteresis loop. These features are described in terms of the self-organized criticality.

Dynamic correlations in the conserved Manna sandpile

We study dynamic correlations for current and mass, as well as the associated power spectra, in the one-dimensional conserved Manna sandpile. We show that, in the thermodynamic limit, the variance of cumulative bond current up to time T grows subdiffusively as T^{1/2-μ} with the exponent μ≥0 depending on the density regimes considered and, likewise, the power spectra of current and mass at low frequency f varies as f^{1/2+μ} and f^{-3/2+μ}, respectively. Our theory predicts that, far from criticality, μ=0 and, near criticality, μ=(β+1)/2ν_{⊥}z>0 with β, ν_{⊥}, and z being the order parameter, correlation length, and dynamic exponents, respectively. The anomalous suppression of fluctuations near criticality signifies a "dynamic hyperuniformity," characterized by a set of fluctuation relations, in which current, mass, and tagged-particle displacement fluctuations are shown to have a precise quantitative relationship with the density-dependent activity (or its derivative). In particular, the relation, D_{s}(ρ[over ¯])=a(ρ[over ¯])/ρ[over ¯], between the self-diffusion coefficient D_{s}(ρ[over ¯]), activity a(ρ[over ¯]) and density ρ[over ¯] explains a previous simulation observation [Eur. Phys. J. B 72, 441 (2009)10.1140/epjb/e2009-00367-0] that, near criticality, the self-diffusion coefficient in the Manna sandpile has the same scaling behavior as the activity.

Self-organized quantization and oscillations on continuous fixed-energy sandpiles

Atmospheric self-organization and activator-inhibitor dynamics in biology provide examples of checkerboardlike spatiotemporal organization. We study a simple model for local activation-inhibition processes. Our model, first introduced in the context of atmospheric moisture dynamics, is a continuous-energy and non-Abelian version of the fixed-energy sandpile model. Each lattice site is populated by a nonnegative real number, its energy. Upon each timestep all sites with energy exceeding a unit threshold redistribute their energy at equal parts to their nearest neighbors. The limit cycle dynamics gives rise to a complex phase diagram in dependence on the mean energy μ: For low μ, all dynamics ceases after few redistribution events. For large μ, the dynamics is well-described as a diffusion process, where the order parameter, spatial variance σ, is removed. States at intermediate μ are dominated by checkerboardlike period-two phases which are however interspersed by much more complex phases of far longer periods. Phases are separated by discontinuous jumps in σ or ∂_{μ}σ-akin to first- and higher-order phase transitions. Overall, the energy landscape is dominated by few energy levels which occur as sharp spikes in the single-site density of states and are robust to noise.

Emergence of oscillations in fixed-energy sandpile models on complex networks

Fixed-energy sandpile (FES) models, introduced to understand the self-organized criticality, show a continuous phase transition between absorbing and active phases. In this work, we study the dynamics of the deterministic FES models on random networks. We observe that close to absorbing transition the density of active nodes oscillates and nodes topple in synchrony. The deterministic toppling rule and the small-world property of random networks lead to the emergence of sustained oscillations. The amplitude of oscillations becomes larger with increasing the value of network randomness. The bifurcation diagram for the density of active nodes is obtained. We use the activity-dependent rewiring rule and show that the interplay between the network structure and the FES dynamics leads to the emergence of a bistable region with a first-order transition between the absorbing and active states. Furthermore during the rewiring, the ordered activation pattern of the nodes is broken, which causes the oscillations to disappear.

On the tail of the branching random walk local time

Consider a critical branching random walk on \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathbb Z^d$$\end{document}Zd, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$d\ge 1$$\end{document}d≥1, started with a single particle at the origin, and let L(x) be the total number of particles that ever visit a vertex x. We study the tail of L(x) under suitable conditions on the offspring distribution. In particular, our results hold if the offspring distribution has an exponential moment.

A New Look at Vaccine Strategies Against PPRV Focused on Adenoviral Candidates

Peste des petits ruminants virus (PPRV) is a virus that mainly infects goats and sheep causing significant economic loss in Africa and Asia, but also posing a serious threat to Europe, as recent outbreaks in Georgia (2016) and Bulgaria (2018) have been reported. In order to carry out the eradication of PPRV, an objective set for 2030 by the Office International des Epizooties (OIE) and the Food and Agriculture Organization of the United Nations (FAO), close collaboration between governments, pharmaceutical companies, farmers and researchers, among others, is needed. Today, more than ever, as seen in the response to the SARS-CoV2 pandemic that we are currently experiencing, these goals are feasible. We summarize in this review the current vaccination approaches against PPRV in the field, discussing their advantages and shortfalls, as well as the development and generation of new vaccination strategies, focusing on the potential use of adenovirus as vaccine platform against PPRV and more broadly against other ruminant pathogens.

Density relaxation in conserved Manna sandpiles

We study relaxation of long-wavelength density perturbations in a one-dimensional conserved Manna sandpile. Far from criticality where correlation length ξ is finite, relaxation of density profiles having wave numbers k→0 is diffusive, with relaxation time τ_{R}∼k^{-2}/D with D being the density-dependent bulk-diffusion coefficient. Near criticality with kξ≳1, the bulk diffusivity diverges and the transport becomes anomalous; accordingly, the relaxation time varies as τ_{R}∼k^{-z}, with the dynamical exponent z=2-(1-β)/ν_{⊥}<2, where β is the critical order-parameter exponent and ν_{⊥} is the critical correlation-length exponent. Relaxation of initially localized density profiles on an infinite critical background exhibits a self-similar structure. In this case, the asymptotic scaling form of the time-dependent density profile is analytically calculated: we find that, at long times t, the width σ of the density perturbation grows anomalously, σ∼t^{w}, with the growth exponent ω=1/(1+β)>1/2. In all cases, theoretical predictions are in reasonably good agreement with simulations.

Asymptotic Height Distribution in High-Dimensional Sandpiles

We give an asymptotic formula for the single-site height distribution of Abelian sandpiles on $$\mathbb {Z}^d$$ Z d as $$d \rightarrow \infty $$ d → ∞ , in terms of $$\mathsf {Poisson}(1)$$ Poisson ( 1 ) probabilities. We provide error estimates.

Mean-field avalanche size exponent for sandpiles on Galton–Watson trees

We show that in Abelian sandpiles on infinite Galton–Watson trees, the probability that the total avalanche has more than t topplings decays as $$t^{-1/2}$$t-1/2. We prove both quenched and annealed bounds, under suitable moment conditions. Our proofs are based on an analysis of the conductance martingale of Morris (Probab Theory Relat Fields 125:259–265, 2003), that was previously used by Lyons et al. (Electron J Probab 13(58):1702–1725, 2008) to study uniform spanning forests on $${\mathbb {Z}^d}$$Zd, $$d\ge 3$$d≥3, and other transient graphs.

Harmonic dynamics of the abelian sandpile

The abelian sandpile is a cellular automaton which serves as the archetypical model to study self-organized criticality, a phenomenon occurring in various biological, physical, and social processes. Its recurrent configurations form an abelian group, whose identity is a fractal composed of self-similar patches. Here, we analyze the evolution of the sandpile identity under harmonic fields of different orders. We show that this evolution corresponds to periodic cycles through the abelian group characterized by the smooth transformation and apparent conservation of the patches constituting the identity. The dynamics induced by second- and third-order harmonics resemble smooth stretchings and translations, respectively, while the ones induced by fourth-order harmonics resemble magnifications and rotations. Based on an extensive analysis of these sandpile dynamics on domains of different size, we conjecture the existence of several scaling limits for infinite domains. Furthermore, we show that the space of harmonic functions provides a set of universal coordinates identifying configurations between different domains, which directly implies that the sandpile group admits a natural renormalization. Finally, we show that the harmonic fields can be induced by simple Markov processes and that the corresponding stochastic dynamics show remarkable robustness. Our results suggest that harmonic fields might split the sandpile group into subsets showing different critical coefficients and that it might be possible to extend the fractal structure of the identity beyond the boundaries of its domain.

Woven Nematic Defects, Skyrmions, and the Abelian Sandpile Model

We show that a fixed set of woven defect lines in a nematic liquid crystal supports a set of nonsingular topological states which can be mapped on to recurrent stable configurations in the Abelian sandpile model or chip-firing game. The physical correspondence between local skyrmion flux and sandpile height is made between the two models. Using a toy model of the elastic energy, we examine the structure of energy minima as a function of topological class and show that the system admits domain wall skyrmion solitons.

Self-organized criticality and pattern emergence through the lens of tropical geometry

Tropical geometry, an established field in pure mathematics, is a place where string theory, mirror symmetry, computational algebra, auction theory, and so forth meet and influence one another. In this paper, we report on our discovery of a tropical model with self-organized criticality (SOC) behavior. Our model is continuous, in contrast to all known models of SOC, and is a certain scaling limit of the sandpile model, the first and archetypical model of SOC. We describe how our model is related to pattern formation and proportional growth phenomena and discuss the dichotomy between continuous and discrete models in several contexts. Our aim in this context is to present an idealized tropical toy model (cf. Turing reaction-diffusion model), requiring further investigation.

Impact of network topology on self-organized criticality

The general mechanisms behind self-organized criticality (SOC) are still unknown. Several microscopic and mean-field theory approaches have been suggested, but they do not explain the dependence of the exponents on the underlying network topology of the SOC system. Here, we first report the phenomena that in the Bak-Tang-Wiesenfeld (BTW) model, sites inside an avalanche area largely return to their original state after the passing of an avalanche, forming, effectively, critically arranged clusters of sites. Then, we hypothesize that SOC relies on the formation process of these clusters, and present a model of such formation. For low-dimensional networks, we show theoretically and in simulation that the exponent of the cluster-size distribution is proportional to the ratio of the fractal dimension of the cluster boundary and the dimensionality of the network. For the BTW model, in our simulations, the exponent of the avalanche-area distribution matched approximately our prediction based on this ratio for two-dimensional networks, but deviated for higher dimensions. We hypothesize a transition from cluster formation to the mean-field theory process with increasing dimensionality. This work sheds light onto the mechanisms behind SOC, particularly, the impact of the network topology.

Statistical investigation of avalanches of three-dimensional small-world networks and their boundary and bulk cross-sections

In many situations we are interested in the propagation of energy in some portions of a three-dimensional system with dilute long-range links. In this paper, a sandpile model is defined on the three-dimensional small-world network with real dissipative boundaries and the energy propagation is studied in three dimensions as well as the two-dimensional cross-sections. Two types of cross-sections are defined in the system, one in the bulk and another in the system boundary. The motivation of this is to make clear how the statistics of the avalanches in the bulk cross-section tend to the statistics of the dissipative avalanches, defined in the boundaries as the concentration of long-range links (α) increases. This trend is numerically shown to be a power law in a manner described in the paper. Two regimes of α are considered in this work. For sufficiently small αs the dominant behavior of the system is just like that of the regular BTW, whereas for the intermediate values the behavior is nontrivial with some exponents that are reported in the paper. It is shown that the spatial extent up to which the statistics is similar to the regular BTW model scales with α just like the dissipative BTW model with the dissipation factor (mass in the corresponding ghost model) m^{2}∼α for the three-dimensional system as well as its two-dimensional cross-sections.

Hydrodynamics, density fluctuations, and universality in conserved stochastic sandpiles

We study conserved stochastic sandpiles (CSSs), which exhibit an active-absorbing phase transition upon tuning density ρ. We demonstrate that a broad class of CSSs possesses a remarkable hydrodynamic structure: There is an Einstein relation σ^{2}(ρ)=χ(ρ)/D(ρ), which connects bulk-diffusion coefficient D(ρ), conductivity χ(ρ), and mass fluctuation, or scaled variance of subsystem mass, σ^{2}(ρ). Consequently, density large-deviations are governed by an equilibrium-like chemical potential μ(ρ)∼lna(ρ), where a(ρ) is the activity in the system. By using the above hydrodynamics, we derive two scaling relations: As Δ=(ρ-ρ_{c})→0^{+}, ρ_{c} being the critical density, (i) the mass fluctuation σ^{2}(ρ)∼Δ^{1-δ} with δ=0 and (ii) the dynamical exponent z=2+(β-1)/ν_{⊥}, expressed in terms of two static exponents β and ν_{⊥} for activity a(ρ)∼Δ^{β} and correlation length ξ∼Δ^{-ν_{⊥}}, respectively. Our results imply that conserved Manna sandpile, a well studied variant of the CSS, belongs to a distinct universality-not that of directed percolation (DP), which, without any conservation law as such, does not obey scaling relation (ii).

Absorbing phase transitions in deterministic fixed-energy sandpile models

We investigate the origin of the difference, which was noticed by Fey et al. [Phys. Rev. Lett. 104, 145703 (2010)PRLTAO0031-900710.1103/PhysRevLett.104.145703], between the steady state density of an Abelian sandpile model (ASM) and the transition point of its corresponding deterministic fixed-energy sandpile model (DFES). Being deterministic, the configuration space of a DFES can be divided into two disjoint classes such that every configuration in one class should evolve into one of absorbing states, whereas no configurations in the other class can reach an absorbing state. Since the two classes are separated in terms of toppling dynamics, the system can be made to exhibit an absorbing phase transition (APT) at various points that depend on the initial probability distribution of the configurations. Furthermore, we show that in general the transition point also depends on whether an infinite-size limit is taken before or after the infinite-time limit. To demonstrate, we numerically study the two-dimensional DFES with Bak-Tang-Wiesenfeld toppling rule (BTW-FES). We confirm that there are indeed many thresholds. Nonetheless, the critical phenomena at various transition points are found to be universal. We furthermore discuss a microscopic absorbing phase transition, or a so-called spreading dynamics, of the BTW-FES, to find that the phase transition in this setting is related to the dynamical isotropic percolation process rather than self-organized criticality. In particular, we argue that choosing recurrent configurations of the corresponding ASM as an initial configuration does not allow for a nontrivial APT in the DFES.

Optimization by Self-Organized Criticality

Self-organized criticality (SOC) is a phenomenon observed in certain complex systems of multiple interacting components, e.g., neural networks, forest fires, and power grids, that produce power-law distributed avalanche sizes. Here, we report the surprising result that the avalanches from an SOC process can be used to solve non-convex optimization problems. To generate avalanches, we use the Abelian sandpile model on a graph that mirrors the graph of the optimization problem. For optimization, we map the avalanche areas onto search patterns for optimization, while the SOC process receives no feedback from the optimization itself. The resulting method can be applied without parameter tuning to a wide range of optimization problems, as demonstrated on three problems: finding the ground-state of an Ising spin glass, graph coloring, and image segmentation. We find that SOC search is more efficient compared to other random search methods, including simulated annealing, and unlike annealing, it is parameter free, thereby eliminating the time-consuming requirement to tune an annealing temperature schedule.

Scaling limit of the odometer in divisible sandpiles

In a recent work Levine et al. (Ann Henri Poincaré 17:1677–1711, 2016. 10.1007/s00023-015-0433-x) prove that the odometer function of a divisible sandpile model on a finite graph can be expressed as a shifted discrete bilaplacian Gaussian field. For the discrete torus, they suggest the possibility that the scaling limit of the odometer may be related to the continuum bilaplacian field. In this work we show that in any dimension the rescaled odometer converges to the continuum bilaplacian field on the unit torus.

Mechanisms of self-organized criticality in social processes of knowledge creation

In online social dynamics, a robust scale invariance appears as a key feature of collaborative efforts that lead to new social value. The underlying empirical data thus offers a unique opportunity to study the origin of self-organized criticality (SOC) in social systems. In contrast to physical systems in the laboratory, various human attributes of the actors play an essential role in the process along with the contents (cognitive, emotional) of the communicated artifacts. As a prototypical example, we consider the social endeavor of knowledge creation via Questions and Answers (Q&A). Using a large empirical data set from one of such Q&A sites and theoretical modeling, we reveal fundamental characteristics of SOC by investigating the temporal correlations at all scales and the role of cognitive contents to the avalanches of the knowledge-creation process. Our analysis shows that the universal social dynamics with power-law inhomogeneities of the actions and delay times provides the primary mechanism for self-tuning towards the critical state; it leads to the long-range correlations and the event clustering in response to the external driving by the arrival of new users. In addition, the involved cognitive contents (systematically annotated in the data and observed in the model) exert important constraints that identify unique classes of the knowledge-creation avalanches. Specifically, besides determining a fine structure of the developing knowledge networks, they affect the values of scaling exponents and the geometry of large avalanches and shape the multifractal spectrum. Furthermore, we find that the level of the activity of the communities that share the knowledge correlates with the fluctuations of the innovation rate, implying that the increase of innovation may serve as the active principle of self-organization. To identify relevant parameters and unravel the role of the network evolution underlying the process in the social system under consideration, we compare the social avalanches to the avalanche sequences occurring in the field-driven physical model of disordered solids, where the factors contributing to the collective dynamics are better understood.

Fractal frontiers of bursts and cracks in a fiber bundle model of creep rupture

We investigate the geometrical structure of breaking bursts generated during the creep rupture of heterogeneous materials. Using a fiber bundle model with localized load sharing we show that bursts are compact geometrical objects; however, their external frontiers have a fractal structure which reflects their growth dynamics. The perimeter fractal dimension of bursts proved to have the universal value 1.25 independent of the external load and of the amount of disorder in the system. We conjecture that according to their geometrical features, breaking bursts fall in the universality class of loop-erased self-avoiding random walks with perimeter fractal dimension 5/4 similar to the avalanches of Abelian sand pile models. The fractal dimension of the growing crack front along which bursts occur proved to increase from 1 to 1.25 as bursts gradually cover the entire front.

Annealed and quenched disorder in sand-pile models with local violation of conservation

In this paper we consider the Bak, Tang, and Wiesenfeld (BTW) sand-pile model with local violation of conservation through annealed and quenched disorder. We have observed that the probability distribution functions of avalanches have two distinct exponents, one of which is associated with the usual BTW model and another one which we propose to belong to a new fixed point; that is, a crossover from the original BTW fixed point to a new fixed point is observed. Through field theoretic calculations, we show that such a perturbation is relevant and takes the system to a new fixed point.

Percolation transitions in the survival of interdependent agents on multiplex networks, catastrophic cascades, and solid-on-solid surface growth

We present an efficient algorithm for simulating percolation transitions of mutually supporting viable clusters on multiplex networks (also known as "catastrophic cascades on interdependent networks"). This algorithm maps the problem onto a solid-on-solid-type model. We use this algorithm to study interdependent agents on duplex Erdös-Rényi (ER) networks and on lattices with dimensions 2, 3, 4, and 5. We obtain surprising results in all these cases, and we correct statements in the literature for ER networks and for two-dimensional lattices. In particular, we find that d=4 is the upper critical dimension and that the percolation transition is continuous for d≤4 but-at least for d≠3-not in the universality class of ordinary percolation. For ER networks we verify that the cluster statistics is exactly described by mean-field theory but find evidence that the cascade process is not. For d=5 we find a first-order transition as for ER networks, but we find also that small clusters have a nontrivial mass distribution that scales at the transition point. Finally, for d=2 with intermediate-range dependency links we propose a scenario that differs from that proposed in W. Li et al. [Phys. Rev. Lett. 108, 228702 (2012)].

Statistical investigation of the cross sections of wave clusters in the three-dimensional Bak-Tang-Wiesenfeld model

We consider the three-dimensional (3D) Bak-Tang-Wiesenfeld model in a cubic lattice. Along with analyzing the 3D problem, the geometrical structure of the two-dimensional (2D) cross section of waves is investigated. By analyzing the statistical observables defined in the cross sections, it is shown that the model in that plane (named as 2D-induced model) is in the critical state and fulfills the finite-size scaling hypothesis. The analysis of the critical loops that are interfaces of the 2D-induced model is of special importance in this paper. Most importantly, we see that their fractal dimension is D(f)=1.387±0.005, which is compatible with the fractal dimension of the external perimeter of geometrical spin clusters of 2D critical Ising model. Some hyperscaling relations between the exponents of the model are proposed and numerically confirmed. We then address the problem of conformal invariance of the mentioned domain walls using Schramm-Lowener evolution (SLE). We found that they are described by SLE with the diffusivity parameter κ=2.8±0.2, nearly consistent with observed fractal dimension.

Unambiguous reconstruction of network structure using avalanche dynamics

A robust method for inferring the structure of networks is presented based on the one-to-one correspondence between the expected composition of cascades of bursts of activity, called crackling noise or avalanches, and the weight matrix. Using a model of neuronal avalanches as a paradigmatic example, we derive this correspondence exactly by calculating the closed-form expression of the joint probability distribution of avalanche sizes obtained by counting separately the number of elements active in each subnetwork during avalanches.

Vaccination with recombinant adenoviruses expressing the peste des petits ruminants virus F or H proteins overcomes viral immunosuppression and induces protective immunity against PPRV challenge in sheep

Peste des petits ruminants (PPR) is a highly contagious disease of small ruminants caused by the Morbillivirus peste des petits ruminants virus (PPRV). Two recombinant replication-defective human adenoviruses serotype 5 (Ad5) expressing either the highly immunogenic fusion protein (F) or hemagglutinin protein (H) from PPRV were used to vaccinate sheep by intramuscular inoculation. Both recombinant adenovirus vaccines elicited PPRV-specific B- and T-cell responses. Thus, neutralizing antibodies were detected in sera from immunized sheep. In addition, we detected a significant antigen specific T-cell response in vaccinated sheep against two different PPRV strains, indicating that the vaccine induced heterologous T cell responses. Importantly, no clinical signs and undetectable virus shedding were observed after virulent PPRV challenge in vaccinated sheep. These vaccines also overcame the T cell immunosuppression induced by PPRV in control animals. The results indicate that these adenovirus constructs could be a promising alternative to current vaccine strategies for the development of PPRV DIVA vaccines.

Stochastic dynamics of lexicon learning in an uncertain and nonuniform world

We study the time taken by a language learner to correctly identify the meaning of all words in a lexicon under conditions where many plausible meanings can be inferred whenever a word is uttered. We show that the most basic form of cross-situational learning--whereby information from multiple episodes is combined to eliminate incorrect meanings--can perform badly when words are learned independently and meanings are drawn from a nonuniform distribution. If learners further assume that no two words share a common meaning, we find a phase transition between a maximally efficient learning regime, where the learning time is reduced to the shortest it can possibly be, and a partially efficient regime where incorrect candidate meanings for words persist at late times. We obtain exact results for the word-learning process through an equivalence to a statistical mechanical problem of enumerating loops in the space of word-meaning mappings.

Power spectrum of mass and activity fluctuations in a sandpile

We consider a directed Abelian sandpile on a strip of size 2×n, driven by adding a grain randomly at the left boundary after every T timesteps. We establish the exact equivalence of the problem of mass fluctuations in the steady state and the number of zeros in the ternary-base representation of the position of a random walker on a ring of size 3^{n}. We find that while the fluctuations of mass have a power spectrum that varies as 1/f for frequencies in the range 3^{-2n}≪f≪1/T, the activity fluctuations in the same frequency range have a power spectrum that is linear in f.

Avalanche frontiers in the dissipative Abelian sandpile model and off-critical Schramm-Loewner evolution

Avalanche frontiers in Abelian sandpile model (ASM) are random simple curves whose continuum limit is known to be a Schramm-Loewner evolution with diffusivity parameter κ=2. In this paper we consider the dissipative ASM and study the statistics of the avalanche and wave frontiers for various rates of dissipation. We examine the scaling behavior of a number of functions, such as the correlation length, the exponent of distribution function of loop lengths, and the gyration radius defined for waves and avalanches. We find that they do scale with the rate of dissipation. Two significant length scales are observed. For length scales much smaller than the correlation length, these curves show properties close to the critical curves, and the corresponding diffusivity parameter is nearly the same as the critical limit. We interpret this as the ultraviolet limit where κ=2 corresponding to c=-2. For length scales much larger than the correlation length, we find that the avalanche frontiers tend to self-avoiding walk, and the corresponding driving function is proportional to the Brownian motion with the diffusivity parameter κ=8/3 corresponding to a field theory with c=0. We interpret this to be the infrared limit of the theory or at least a crossover.

Numerical study of the correspondence between the dissipative and fixed-energy Abelian sandpile models

We consider the Abelian sandpile model (ASM) on the square lattice with a single dissipative site (sink). Particles are added one by one per unit time at random sites and the resulting density of particles is calculated as a function of time. We observe different scenarios of evolution depending on the value of initial uniform density (height) h(0). During the first stage of the evolution, the density of particles increases linearly. Reaching a critical density ρ(c)(h(0)), the system changes its behavior and relaxes exponentially to the stationary state of the ASM with density ρ(s). Considering initial heights -1 ≤ h(0) ≤ 4, we observe a dramatic decrease of the difference ρ(c)(h(0)) - ρ(s) when h(0) is zero or negative. In parallel with the ASM, we consider the conservative fixed energy sandpile (FES). The extensive Monte Carlo simulations show that the threshold density ρ(th)(h(0)) of the FES converges rapidly to ρ(s) for h(0) < 1.

Abelian deterministic self-organized-criticality model: complex dynamics of avalanche waves

The aim of this study is to investigate a wave dynamics and a size scaling of avalanches which were created by the mathematical model [J. Černák, Phys. Rev. E 65, 046141 (2002)]. Numerical simulations were carried out on a two-dimensional lattice L×L in which two constant thresholds E(c)(I) = 4 and E(c)(II) > E(c)(I) were randomly distributed. The density of sites c of the thresholds E(c)(II) and threshold E(c)(II) are parameters of the model. Autocorrelations of avalanche size waves, Hurst exponents, avalanche structures, and avalanche size moments were determined for several densities c and thresholds E(c)(II). The results show correlated avalanche size waves and multifractal scaling of avalanche sizes not only for specific conditions, densities c = 0.0,1.0 and thresholds 8 ≤ E(c)(II) ≤ 32, in which relaxation rules were precisely balanced, but also for more general conditions, densities 0.0 < c < 1.0 and thresholds 8 ≤ E(c)(II) ≤ 32, in which relaxation rules were unbalanced. The results suggest that the hypothesis of a precise relaxation balance could be a specific case of a more general rule.

Approach to criticality in sandpiles

A popular theory of self-organized criticality predicts that the stationary density of the Abelian sandpile model equals the threshold density of the corresponding fixed-energy sandpile. We recently announced that this "density conjecture" is false when the underlying graph is any of Z2, the complete graph K(n), the Cayley tree, the ladder graph, the bracelet graph, or the flower graph. In this paper, we substantiate this claim by rigorous proof and extensive simulations. We show that driven-dissipative sandpiles continue to evolve even after a constant fraction of the sand has been lost at the sink. Nevertheless, we do find (and prove) a relationship between the two models: the threshold density of the fixed-energy sandpile is the point at which the driven-dissipative sandpile begins to lose a macroscopic amount of sand to the sink.

Statistical properties of directed avalanches

A two-dimensional directed stochastic sandpile model is studied both numerically and analytically. One of the known analytical approaches is extended by considering general stochastic toppling rules. The probability density distribution for the first-passage time of stochastic process described by a nonlinear Langevin equation with power-law dependence of the diffusion coefficient is obtained. Large-scale Monte Carlo simulations are performed with the aim to analyze statistical properties of the avalanches, such as the asymmetry between the initial and final stages, scaling of voids and the width of the thickest branch. Comparison with random walks description is drawn and different plausible scenarios for the avalanche evolution and the scaling exponents are suggested.

Driving sandpiles to criticality and beyond

A popular theory of self-organized criticality relates driven dissipative systems to systems with conservation. This theory predicts that the stationary density of the Abelian sandpile model equals the threshold density of the fixed-energy sandpile. We refute this prediction for a wide variety of underlying graphs, including the square grid. Driven dissipative sandpiles continue to evolve even after reaching criticality. This result casts doubt on the validity of using fixed-energy sandpiles to explore the critical behavior of the Abelian sandpile model at stationarity.

Renormalization-group approach to the Manna sandpile

In this paper a renormalization group (RG) scheme for the q -state Manna model is proposed based on the stochastic characteristic and the similarity of topplings at different scales. A full enumeration of the RG evolution events inside a 2 x 2 RG cell for q=2 , 3, and 4 was carried out. A fixed point analysis shows that the resulting height probabilities are very close to the results obtained by the numerical simulations. The calculations of the toppling number exponent tau and the dynamical exponent z are also provided. It was found that the RG values of tau and z for q=4 are very close to the simulation values.

Patterned and disordered continuous Abelian sandpile model

We study critical properties of the continuous Abelian sandpile model with anisotropies in toppling rules that produce ordered patterns on it. Also, we consider the continuous directed sandpile model perturbed by a weak quenched randomness, study critical behavior of the model using perturbative conformal field theory, and show that the model has a random fixed point.

Branching process in a stochastic extremal model

We considered a stochastic version of the Bak-Sneppen model (SBSM) of ecological evolution where the number M of sites mutated in a mutation event is restricted to only two. Here the mutation zone consists of only one site and this site is randomly selected from the neighboring sites at every mutation event in an annealed fashion. The critical behavior of the SBSM is found to be the same as the BS model in dimensions d=1 and 2. However on the scale-free graphs the critical fitness value is nonzero even in the thermodynamic limit but the critical behavior is mean-field like. Finally M has been made even smaller than two by probabilistically updating the mutation zone, which also shows the original BS model behavior. We conjecture that a SBSM on any arbitrary graph with any small branching factor greater than unity will lead to a self-organized critical state.