Altered growth and death in dilution-based viral predation assays

Viral lysis of phytoplankton is one of the most common forms of death on Earth. Building on an assay used extensively to assess rates of phytoplankton loss to predation by grazers, lysis rates are increasingly quantified through dilution-based techniques. In this approach, dilution of viruses and hosts are expected to reduce infection rates and thus increase host net growth rates (i.e., accumulation rates). The difference between diluted and undiluted host growth rates is interpreted as a measurable proxy for the rate of viral lytic death. These assays are usually conducted in volumes ≥ 1 L. To increase throughput, we implemented a miniaturized, high-throughput, high-replication, flow cytometric microplate dilution assay to measure viral lysis in environmental samples sourced from a suburban pond and the North Atlantic Ocean. The most notable outcome we observed was a decline in phytoplankton densities that was exacerbated by dilution, instead of the increased growth rates expected from lowered virus-phytoplankton encounters. We sought to explain this counterintuitive outcome using theoretical, environmental, and experimental analyses. Our study shows that, while die-offs could be partly explained by a 'plate effect' due to small incubation volumes and cells adhering to walls, the declines in phytoplankton densities are not volume-dependent. Rather, they are driven by many density- and physiology-dependent effects of dilution on predation pressure, nutrient limitation, and growth, all of which violate the original assumptions of dilution assays. As these effects are volume-independent, these processes likely occur in all dilution assays that our analyses show to be remarkably sensitive to dilution-altered phytoplankton growth and insensitive to actual predation pressure. Incorporating altered growth as well as predation, we present a logical framework that categorizes locations by the relative dominance of these mechanisms, with general applicability to dilution-based assays.

Depinning in the quenched Kardar-Parisi-Zhang class. I. Mappings, simulations, and algorithm

Depinning of elastic systems advancing on disordered media can usually be described by the quenched Edwards-Wilkinson equation (qEW). However, additional ingredients such as anharmonicity and forces that cannot be derived from a potential energy may generate a different scaling behavior at depinning. The most experimentally relevant is the Kardar-Parisi-Zhang (KPZ) term, proportional to the square of the slope at each site, which drives the critical behavior into the so-called quenched KPZ (qKPZ) universality class. We study this universality class both numerically and analytically: by using exact mappings we show that at least for d=1,2 this class encompasses not only the qKPZ equation itself, but also anharmonic depinning and a well-known class of cellular automata introduced by Tang and Leschhorn. We develop scaling arguments for all critical exponents, including size and duration of avalanches. The scale is set by the confining potential strength m^{2}. This allows us to estimate numerically these exponents as well as the m-dependent effective force correlator Δ(w), and its correlation length ρ:=Δ(0)/|Δ^{'}(0)|. Finally, we present an algorithm to numerically estimate the effective (m-dependent) elasticity c, and the effective KPZ nonlinearity λ. This allows us to define a dimensionless universal KPZ amplitude A:=ρλ/c, which takes the value A=1.10(2) in all systems considered in d=1. This proves that qKPZ is the effective field theory for all these models. Our work paves the way for a deeper understanding of depinning in the qKPZ class, and in particular, for the construction of a field theory that we describe in a companion paper.

Alternative stable ecological states observed after a biological invasion

Although biological invasions play an important role in ecosystem change worldwide, little is known about how invasions are influenced by local abiotic stressors. Broadly, abiotic stressors can cause large-scale community changes in an ecosystem that influence its resilience. The possibility for these stressors to increase as global changes intensify highlights the pressing need to understand and characterize the effects that abiotic drivers may have on the dynamics and composition of a community. Here, we analyzed 26 years of weekly abundance data using the theory of regime shifts to understand how the structure of a resident community of dung beetles (composed of dweller and tunneler functional groups) responds to climatic changes in the presence of the invasive tunneler Digitonthophagus gazella. Although the community showed an initial dominance by the invader that decreased over time, the theory of regime shifts reveals the possibility of an ecological transition driven by climate factors (summarized here in a climatic index that combines minimum temperature and relative humidity). Mid and low values of the driver led to the existence of two alternative stable states for the community structure (i.e. dominance of either dwellers or tunnelers for similar values of the climatic driver), whereas large values of the driver led to the single dominance by tunnelers. We also quantified the stability of these states against climatic changes (resilience), which provides insight on the conditions under which the success of an invasion and/or the recovery of the previous status quo for the ecosystem are expected. Our approach can help understand the role of climatic changes in community responses, and improve our capacity to deal with regime shifts caused by the introduction of exotic species in new ecosystems.

Unconstrained coevolution of bacterial size and the latent period of plastic phage

Viruses play critical roles in the dynamics of microbial communities. Lytic viruses, for example, kill significant fractions of autotrophic and heterotrophic microbes daily. The dynamic interplay between viruses and microbes results from an overlap of physiological, ecological, and evolutionary responses: environmental changes trigger host physiological changes, affecting the ecological interactions of host and virus and, ultimately, the evolutionary pressures influencing the two populations. Recent theoretical work studied how the dependence of viral traits on host physiology (viral plasticity) affects the evolutionarily stable host cell size and viral infection time emerging from coevolution. Here, we broaden the scope of the framework to consider any coevolutionary outcome, including potential evolutionary collapses of the system. We used the case study of Escherichia coli and T-like viruses under chemostat conditions, but the framework can be adapted to any microbe-virus system. Oligotrophic conditions led to smaller, lower-quality but more abundant hosts, and infections that were longer but produced a reduced viral offspring. Conversely, eutrophic conditions resulted in fewer but larger higher-quality hosts, and shorter but more productive infections. The virus influenced host evolution decreasing host size more noticeably for low than for high dilution rates, and for high than for low nutrient input concentration. For low dilution rates, the emergent infection time minimized host need/use, but higher dilution led to an opportunistic strategy that shortened the duration of infections. System collapses driven by evolution resulted from host failure to adapt quickly enough to the evolving virus. Our results contribute to understanding the eco-evolutionary dynamics of microbes and virus, and to improving the predictability of current models for host-virus interactions. The large quantitative and qualitative differences observed with respect to a classic description (in which viral traits are assumed to be constant) highlights the importance of including viral plasticity in theories describing short- and long-term host-virus dynamics.

Resource availability and heterogeneity shape the self-organisation of regular spatial patterning

Explaining large-scale ordered patterns and their effects on ecosystem functioning is a fundamental and controversial challenge in ecology. Here, we coupled empirical and theoretical approaches to explore how competition and spatial heterogeneity govern the regularity of colony dispersion in fungus-farming termites. Individuals from different colonies fought fiercely, and inter-nest distances were greater when nests were large and resources scarce-as expected if competition is strong, large colonies require more resources and foraging area scales with resource availability. Building these principles into a model of inter-colony competition showed that highly ordered patterns emerged under high resource availability and low resource heterogeneity. Analysis of this dynamical model provided novel insights into the mechanisms that modulate pattern regularity and the emergent effects of these patterns on system-wide productivity. Our results show how environmental context shapes pattern formation by social-insect ecosystem engineers, which offers one explanation for the marked variability observed across ecosystems.

Evolutionarily Stable Coevolution Between a Plastic Lytic Virus and Its Microbial Host

Hosts influence and are influenced by viral replication. Cell size, for example, is a fundamental trait for microbial hosts that can not only alter the probability of viral adsorption, but also constrain the host physiological processes that the virus relies on to replicate. This intrinsic connection can affect the fitness of both host and virus, and therefore their mutual evolution. Here, we study the coevolution of bacterial hosts and their viruses by considering the dependence of viral performance on the host physiological state (viral plasticity). To this end, we modified a standard host-lytic phage model to include viral plasticity, and compared the coevolutionary strategies emerging under different scenarios, including cases in which only the virus or the host evolve. For all cases, we also obtained the evolutionary prediction of the traditional version of the model, which assumes a non-plastic virus. Our results reveal that the presence of the virus leads to an increase in host size and growth rate in the long term, which benefits both interacting populations. Our results also show that viral plasticity and evolution influence the classic host quality-quantity trade-off. Poor nutrient environments lead to abundant low-quality hosts, which tends to increase viral infection time. Conversely, richer nutrient environments lead to fewer but high-quality hosts, which decrease viral infection time. Our results can contribute to advancing our understanding of the microbial response to changing environments. For instance, both cell size and viral-induced mortality are essential factors that determine the structure and dynamics of the marine microbial community, and therefore our study can improve predictions of how marine ecosystems respond to environmental change. Our study can also help devise more reliable strategies to use phage to, for example, fight bacterial infections.

Shape of species climate response curves affects community response to climate change

Understanding how community composition is reshaped by changing climate is important for interpreting and predicting patterns of community assembly through time or across space. Community composition often does not perfectly correspond to expectations from current environmental conditions, leading to community-climate mismatches. Here, we combine data analysis and theory development to explore how species climate response curves affect the community response to climate change. We show that strong mismatches between community and climate can appear in the absence of demographic delays or limited species pools. Communities simulated using species response curves showed temporal changes of similar magnitude to those observed in natural communities of fishes and plankton, suggesting no overall delays in community change despite substantial unexplained variation from community assembly and other processes. Our approach can be considered as a null model that will be important to use when interpreting observed community responses to climate change and variability.

Impact of Lytic Phages on Phosphorus- vs. Nitrogen-Limited Marine Microbes

Lytic viruses kill almost 20% of marine bacteria every day, re-routing nutrients away from the higher trophic levels of the marine food web and back in the microbial loop. Importantly, the effect of this inflow of key elements on the ecosystem depends on the nutrient requirements of bacteria as well as on the elemental composition of the viruses that infect them. Therefore, the influence of viruses on the ecosystem could vary depending on which nutrient is limiting. In this paper, we considered an existing multitrophic model (nutrient, bacteria, zooplankton, and viruses) that accounts for nitrogen limitation, and developed a phosphorus-limited version to assess whether the limiting nutrient alters the role of viruses in the ecosystem. For both versions, we evaluated the stationary state of the system with and without viruses. In agreement with existing results, nutrient release increased with viruses for nitrogen–limited systems, while zooplankton abundance and export to higher trophic levels decreased. We found this to be true also for phosphorus-limited systems, although nutrient release increased less than in nitrogen-limited systems. The latter supports a nutrient-specific response of the ecosystem to viruses. Bacterial concentration decreased in the phosphorus-limited system but increased in most nitrogen-limited cases due to a switch from mostly bottom-up to entirely top-down control by viruses. Our results also show that viral concentration is best predicted by a power-law of bacterial concentration with exponent different from 1. Finally, we found a positive correlation between carbon export and viruses regardless of the limiting nutrient, which led us to suggest viral abundance as a predictor of carbon sink.

Evolution in the Debian GNU/Linux software network: analogies and differences with gene regulatory networks

Biological networks exhibit intricate architectures deemed to be crucial for their functionality. In particular, gene regulatory networks, which play a key role in information processing in the cell, display non-trivial architectural features such as scale-free degree distributions, high modularity and low average distance between connected genes. Such networks result from complex evolutionary and adaptive processes difficult to track down empirically. On the other hand, there exists detailed information on the developmental (or evolutionary) stages of open-software networks that result from self-organized growth across versions. Here, we study the evolution of the Debian GNU/Linux software network, focusing on the changes of key structural and statistical features over time. Our results show that evolution has led to a network structure in which the out-degree distribution is scale-free and the in-degree distribution is a stretched exponential. In addition, while modularity, directionality of information flow, and average distance between elements grew, vulnerability decreased over time. These features resemble closely those currently shown by gene regulatory networks, suggesting the existence of common adaptive pathways for the architectural design of information-processing networks. Differences in other hierarchical aspects point to system-specific solutions to similar evolutionary challenges.

Temporal phosphate gradients reveal diverse acclimation responses in phytoplankton phosphate uptake

Phytoplankton face environmental nutrient variations that occur in the dynamic upper layers of the ocean. Phytoplankton cells are able to rapidly acclimate to nutrient fluctuations by adjusting their nutrient-uptake system and metabolism. Disentangling these acclimation responses is a critical step in bridging the gap between phytoplankton cellular physiology and community ecology. Here, we analyzed the dynamics of phosphate (P) uptake acclimation responses along different P temporal gradients by using batch cultures of the diatom Phaeodactylum tricornutum. We employed a multidisciplinary approach that combined nutrient-uptake bioassays, transcriptomic analysis, and mathematical models. Our results indicated that cells increase their maximum nutrient-uptake rate (V_max) both in response to P pulses and strong phosphorus limitation. The upregulation of three genes coding for different P transporters in cells experiencing low intracellular phosphorus levels supported some of the observed V_max variations. In addition, our mathematical model reproduced the empirical V_max patterns by including two types of P transporters upregulated at medium-high environmental and low intracellular phosphorus levels, respectively. Our results highlight the existence of a sequence of acclimation stages along the phosphate continuum that can be understood as a succession of acclimation responses. We provide a novel conceptual framework that can contribute to integrating and understanding the dynamics and wide diversity of acclimation responses developed by phytoplankton.

Ecological and Evolutionary Consequences of Viral Plasticity

Viruses use the host machinery to replicate, and their performance thus depends on the host's physiological state. For bacteriophages, this link between host and viral performance has been characterized empirically and with intracellular theories. Such theories are too detailed to be included in models that study host-phage interactions in the long term, which hinders our understanding of systems that range from pathogens infecting gut bacteria to marine phage shaping the oceans. Here, we combined data and models to study the short- and long-term consequences that host physiology has on bacteriophage performance. We compiled data showing the dependence of lytic-phage traits on host growth rate (referred to as viral phenotypic plasticity) to deduce simple expressions that represent such plasticity. Including these expressions in a standard host-phage model allowed us to understand mechanistically how viral plasticity affects emergent evolutionary strategies and the population dynamics associated with different environmental scenarios including, for example, nutrient pulses or host starvation. Moreover, we show that plasticity on the offspring number drives the phage ecological and evolutionary dynamics by reinforcing feedbacks between host, virus, and environment. Standard models neglect viral plasticity, which therefore handicaps their predictive ability in realistic scenarios. Our results highlight the importance of viral plasticity to unravel host-phage interactions and the need of laboratory and field experiments to characterize viral plastic responses across systems.

Cell adhesion and fluid flow jointly initiate genotype spatial distribution in biofilms

Biofilms are microbial collectives that occupy a diverse array of surfaces. It is well known that the function and evolution of biofilms are strongly influenced by the spatial arrangement of different strains and species within them, but how spatiotemporal distributions of different genotypes in biofilm populations originate is still underexplored. Here, we study the origins of biofilm genetic structure by combining model development, numerical simulations, and microfluidic experiments using the human pathogen Vibrio cholerae. Using spatial correlation functions to quantify the differences between emergent cell lineage segregation patterns, we find that strong adhesion often, but not always, maximizes the size of clonal cell clusters on flat surfaces. Counterintuitively, our model predicts that, under some conditions, investing in adhesion can reduce rather than increase clonal group size. Our results emphasize that a complex interaction between fluid flow and cell adhesiveness can underlie emergent patterns of biofilm genetic structure. This structure, in turn, has an outsize influence on how biofilm-dwelling populations function and evolve.

Eco-evolutionary Red Queen dynamics regulate biodiversity in a metabolite-driven microbial system

The Red Queen Hypothesis proposes that perpetual co-evolution among organisms can result from purely biotic drivers. After more than four decades, there is no satisfactory understanding as to which mechanisms trigger Red Queen dynamics or their implications for ecosystem features such as biodiversity. One reason for such a knowledge gap is that typical models are complicated theories where limit cycles represent an idealized Red Queen, and therefore cannot be used to devise experimental setups. Here, we bridge this gap by introducing a simple model for microbial systems able to show Red Queen dynamics. We explore diverse biotic sources that can drive the emergence of the Red Queen and that have the potential to be found in nature or to be replicated in the laboratory. Our model enables an analytical understanding of how Red Queen dynamics emerge in our setup, and the translation of model terms and phenomenology into general underlying mechanisms. We observe, for example, that in our system the Red Queen offers opportunities for the increase of biodiversity by facilitating challenging conditions for intraspecific dominance, whereas stasis tends to homogenize the system. Our results can be used to design and engineer experimental microbial systems showing Red Queen dynamics.

A theoretical foundation for multi-scale regular vegetation patterns

Self-organized regular vegetation patterns are widespread and thought to mediate ecosystem functions such as productivity and robustness, but the mechanisms underlying their origin and maintenance remain disputed. Particularly controversial are landscapes of overdispersed (evenly spaced) elements, such as North American Mima mounds, Brazilian murundus, South African heuweltjies, and, famously, Namibian fairy circles. Two competing hypotheses are currently debated. On the one hand, models of scale-dependent feedbacks, whereby plants facilitate neighbours while competing with distant individuals, can reproduce various regular patterns identified in satellite imagery. Owing to deep theoretical roots and apparent generality, scale-dependent feedbacks are widely viewed as a unifying and near-universal principle of regular-pattern formation despite scant empirical evidence. On the other hand, many overdispersed vegetation patterns worldwide have been attributed to subterranean ecosystem engineers such as termites, ants, and rodents. Although potentially consistent with territorial competition, this interpretation has been challenged theoretically and empirically and (unlike scale-dependent feedbacks) lacks a unifying dynamical theory, fuelling scepticism about its plausibility and generality. Here we provide a general theoretical foundation for self-organization of social-insect colonies, validated using data from four continents, which demonstrates that intraspecific competition between territorial animals can generate the large-scale hexagonal regularity of these patterns. However, this mechanism is not mutually exclusive with scale-dependent feedbacks. Using Namib Desert fairy circles as a case study, we present field data showing that these landscapes exhibit multi-scale patterning-previously undocumented in this system-that cannot be explained by either mechanism in isolation. These multi-scale patterns and other emergent properties, such as enhanced resistance to and recovery from drought, instead arise from dynamic interactions in our theoretical framework, which couples both mechanisms. The potentially global extent of animal-induced regularity in vegetation-which can modulate other patterning processes in functionally important ways-emphasizes the need to integrate multiple mechanisms of ecological self-organization.

Eluding catastrophic shifts

Transitions between regimes with radically different properties are ubiquitous in nature. Such transitions can occur either smoothly or in an abrupt and catastrophic fashion. Important examples of the latter can be found in ecology, climate sciences, and economics, to name a few, where regime shifts have catastrophic consequences that are mostly irreversible (e.g., desertification, coral reef collapses, and market crashes). Predicting and preventing these abrupt transitions remains a challenging and important task. Usually, simple deterministic equations are used to model and rationalize these complex situations. However, stochastic effects might have a profound effect. Here we use 1D and 2D spatially explicit models to show that intrinsic (demographic) stochasticity can alter deterministic predictions dramatically, especially in the presence of other realistic features such as limited mobility or spatial heterogeneity. In particular, these ingredients can alter the possibility of catastrophic shifts by giving rise to much smoother and easily reversible continuous ones. The ideas presented here can help further understand catastrophic shifts and contribute to the discussion about the possibility of preventing such shifts to minimize their disruptive ecological, economic, and societal consequences.

Ecological feedbacks. Termite mounds can increase the robustness of dryland ecosystems to climatic change

Self-organized spatial vegetation patterning is widespread and has been described using models of scale-dependent feedback between plants and water on homogeneous substrates. As rainfall decreases, these models yield a characteristic sequence of patterns with increasingly sparse vegetation, followed by sudden collapse to desert. Thus, the final, spot-like pattern may provide early warning for such catastrophic shifts. In many arid ecosystems, however, termite nests impart substrate heterogeneity by altering soil properties, thereby enhancing plant growth. We show that termite-induced heterogeneity interacts with scale-dependent feedbacks to produce vegetation patterns at different spatial grains. Although the coarse-grained patterning resembles that created by scale-dependent feedback alone, it does not indicate imminent desertification. Rather, mound-field landscapes are more robust to aridity, suggesting that termites may help stabilize ecosystems under global change.

Quenched disorder forbids discontinuous transitions in nonequilibrium low-dimensional systems

Quenched disorder affects significantly the behavior of phase transitions. The Imry-Ma-Aizenman-Wehr-Berker argument prohibits first-order or discontinuous transitions and their concomitant phase coexistence in low-dimensional equilibrium systems in the presence of random fields. Instead, discontinuous transitions become rounded or even continuous once disorder is introduced. Here we show that phase coexistence and first-order phase transitions are also precluded in nonequilibrium low-dimensional systems with quenched disorder: discontinuous transitions in two-dimensional systems with absorbing states become continuous in the presence of quenched disorder. We also study the universal features of this disorder-induced criticality and find them to be compatible with the universality class of the directed percolation with quenched disorder. Thus, we conclude that first-order transitions do not exist in low-dimensional disordered systems, not even in genuinely nonequilibrium systems with absorbing states.

Evolutionary comparison between viral lysis rate and latent period

Marine viruses shape the structure of the microbial community. They are, thus, a key determinant of the most important biogeochemical cycles in the planet. Therefore, a correct description of the ecological and evolutionary behavior of these viruses is essential to make reliable predictions about their role in marine ecosystems. The infection cycle, for example, is indistinctly modeled in two very different ways. In one representation, the process is described including explicitly a fixed delay between infection and offspring release. In the other, the offspring are released at exponentially distributed times according to a fixed release rate. By considering obvious quantitative differences pointed out in the past, the latter description is widely used as a simplification of the former. However, it is still unclear how the dichotomy "delay versus rate description" affects long-term predictions of host-virus interaction models. Here, we study the ecological and evolutionary implications of using one or the other approaches, applied to marine microbes. To this end, we use mathematical and eco-evolutionary computational analysis. We show that the rate model exhibits improved competitive abilities from both ecological and evolutionary perspectives in steady environments. However, rate-based descriptions can fail to describe properly long-term microbe-virus interactions. Moreover, additional information about trade-offs between life-history traits is needed in order to choose the most reliable representation for oceanic bacteriophage dynamics. This result affects deeply most of the marine ecosystem models that include viruses, especially when used to answer evolutionary questions.

Dynamic model of flexible phytoplankton nutrient uptake

The metabolic machinery of marine microbes can be remarkably plastic, allowing organisms to persist under extreme nutrient limitation. With some exceptions, most theoretical approaches to nutrient uptake in phytoplankton are largely dominated by the classic Michaelis-Menten (MM) uptake functional form, whose constant parameters cannot account for the observed plasticity in the uptake apparatus. Following seminal ideas by earlier researchers, we propose a simple cell-level model based on a dynamic view of the uptake process whereby the cell can regulate the synthesis of uptake proteins in response to changes in both internal and external nutrient concentrations. In our flexible approach, the maximum uptake rate and nutrient affinity increase monotonically as the external nutrient concentration decreases. For low to medium nutrient availability, our model predicts uptake and growth rates larger than the classic MM counterparts, while matching the classic MM results for large nutrient concentrations. These results have important consequences for global coupled models of ocean circulation and biogeochemistry, which lack this regulatory mechanism and are thus likely to underestimate phytoplankton abundances and growth rates in oligotrophic regions of the ocean.

Temporal Griffiths phases

Disorder is an unavoidable ingredient of real systems. Spatial disorder generates Griffiths phases (GPs) which, in analogy to critical points, are characterized by a slow relaxation of the order parameter and divergences of quantities such as the susceptibility. However, these singularities appear in an extended region of the parameter space and not just at a (critical) point, i.e., there is generic scale invariance. Here, we study the effects of temporal disorder, focusing on systems with absorbing states. We show that for dimensions d≥2 there are Temporal Griffiths phases (TGPs) characterized by generic power-law scaling of some magnitudes and generic divergences of the susceptibility. TGPs turn out to be a counterpart of GPs, but with space and time playing reversed roles. TGPs constitute a unifying concept, shedding light on the nontrivial effects of temporal disorder.

Cusps, self-organization, and absorbing states

Elastic interfaces embedded in (quenched) random media exhibit metastability and stick-slip dynamics. These nontrivial dynamical features have been shown to be associated with cusp singularities of the coarse-grained disorder correlator. Here we show that annealed systems with many absorbing states and a conservation law but no quenched disorder exhibit identical cusps. On the other hand, similar nonconserved systems in the directed percolation class are also shown to exhibit cusps but of a different type. These results are obtained both by a recent method to explicitly measure disorder correlators and by defining an alternative new protocol inspired by self-organized criticality, which opens the door to easily accessible experimental realizations.

Simplest nonequilibrium phase transition into an absorbing state

We study in further detail particle models displaying a boundary-induced absorbing state phase transition [Deloubrière and van Wijland Phys. Rev. E 65, 046104 (2002) and Barato and Hinrichsen, Phys. Rev. Lett. 100, 165701 (2008)]. These are one-dimensional systems consisting of a single site (the boundary) where creation and annihilation of particles occur, and a bulk where particles move diffusively. We study different versions of these models and confirm that, except for one exactly solvable bosonic variant exhibiting a discontinuous transition and trivial exponents, all the others display nontrivial behavior, with critical exponents differing from their mean-field values, representing a universality class. Finally, the relation of these systems with a (0+1)-dimensional non-Markovian process is discussed.

Confirming and extending the hypothesis of universality in sandpiles

Stochastic sandpiles self-organize to an absorbing-state critical point with scaling behavior different from directed percolation (DP) and characterized by the presence of an additional conservation law. This is usually called the C-DP or Manna universality class. There remains, however, an exception to this universality principle: a sandpile automaton introduced by Maslov and Zhang, which was claimed to be in the DP class despite the existence of a conservation law. We show, by means of careful numerical simulations as well as by constructing and analyzing a field theory, that (contrarily to what was previously thought) this sandpile is also in the C-DP or Manna class. This confirms the hypothesis of universality for stochastic sandpiles and gives rise to a fully coherent picture of self-organized criticality in systems with conservation. In passing, we obtain a number of results for the C-DP class and introduce a strategy to easily discriminate between DP and C-DP scaling.

Absorbing states and elastic interfaces in random media: two equivalent descriptions of self-organized criticality

We elucidate a long-standing puzzle about the nonequilibrium universality classes describing self-organized criticality in sandpile models. We show that depinning transitions of linear interfaces in random media and absorbing phase transitions (with a conserved nondiffusive field) are two equivalent languages to describe sandpile criticality. This is so despite the fact that local roughening properties can be radically different in the two pictures, as explained here. Experimental implications of our work as well as promising paths for future theoretical investigations are also discussed.

Sticky grains do not change the universality class of isotropic sandpiles

We revisit the sandpile model with "sticky" grains introduced by Mohanty and Dhar [Phys. Rev. Lett. 89, 104303 (2002)] whose scaling properties were claimed to be generically in the universality class of directed percolation for both isotropic and directed models. While for directed models this conclusion is unquestionable, for isotropic models we present strong evidence that the asymptotic scaling in the self-organized regime (in which a stationary critical state exists in the limit of slow driving and vanishing dissipation rate) is, like other stochastic sandpiles, generically in the Manna universality class. This conclusion is drawn from extensive Monte Carlo simulations, and is strengthened by the analysis of the Langevin equations (proposed by the same authors to account for this problem), argued to converge upon coarse-graining to the well-established set of Langevin equations for the Manna class.