Showing 1–50 of 70 results for author: Pruessner, G

Memoryless Chemotaxis with Discrete Cues

Authors: Jacob Knight, Paula García-Galindo, Johannes Pausch, Gunnar Pruessner

Abstract: A wide array of biological systems can navigate in shallow gradients of chemoattractant with remarkable precision. Whilst previous approaches model such systems using coarse-grained chemical density profiles, we construct a dynamical model consisting of a chemotactic cell responding to discrete cue particles. For a cell without internal memory, we derive an effective velocity with which the cell a… ▽ More A wide array of biological systems can navigate in shallow gradients of chemoattractant with remarkable precision. Whilst previous approaches model such systems using coarse-grained chemical density profiles, we construct a dynamical model consisting of a chemotactic cell responding to discrete cue particles. For a cell without internal memory, we derive an effective velocity with which the cell approaches a point source of cue particles. We find that the effective velocity becomes negative beyond some homing radius, which represents an upper bound on the distance within which chemotaxis can be reliably performed. This work lays the foundation for the analytical characterisation of more detailed models of chemotaxis. △ Less

Submitted 18 December, 2023; originally announced December 2023.

Comments: 8 pages, 5 figures

Correlation, crossover and broken scaling in the Abelian Manna Model

Authors: Letian Chen, Hoai Nguyen Huynh, Gunnar Pruessner

Abstract: The role of correlations in self-organised critical (SOC) phenomena is investigated by studying the Abelian Manna Model (AMM) in two dimensions. Local correlations of the debris left behind after avalanches are destroyed by re-arranging particles on the lattice between avalanches, without changing the one-point particle density. It is found that the spatial correlations are not relevant to small a… ▽ More The role of correlations in self-organised critical (SOC) phenomena is investigated by studying the Abelian Manna Model (AMM) in two dimensions. Local correlations of the debris left behind after avalanches are destroyed by re-arranging particles on the lattice between avalanches, without changing the one-point particle density. It is found that the spatial correlations are not relevant to small avalanches, while changing the scaling of the large (system-wide) ones, yielding a crossover in the model's scaling behaviour. This crossover breaks the simple scaling observed in normal SOC. △ Less

Submitted 6 November, 2023; originally announced November 2023.

Field Theory of Active Brownian Particles in Potentials

Authors: Ziluo Zhang, Lili Fehértói-Nagy, Maria Polackova, Gunnar Pruessner

Abstract: The Active Brownian Particle (ABP) model exemplifies a wide class of active matter particles. In this work, we demonstrate how this model can be cast into a field theory in both two and three dimensions. Our aim is manifold: we wish both to extract useful features of the system, as well as to build a framework which can be used to study more complex systems involving ABPs, such as those involving… ▽ More The Active Brownian Particle (ABP) model exemplifies a wide class of active matter particles. In this work, we demonstrate how this model can be cast into a field theory in both two and three dimensions. Our aim is manifold: we wish both to extract useful features of the system, as well as to build a framework which can be used to study more complex systems involving ABPs, such as those involving interaction. Using the two-dimensional model as a template, we calculate the mean squared displacement exactly, and the one-point density in an external potential perturbatively. We show how the effective diffusion constant appears in the barometric density formula to leading order, and determine the corrections to it. We repeat the calculation in three dimensions, clearly a more challenging setup. Comparing different ways to capture the self-propulsion, we find that its perturbative treatment results in more tractable derivations without loss of exactness, where this is accessible. △ Less

Submitted 23 December, 2022; originally announced December 2022.

Field theories of active particle systems and their entropy production

Authors: Gunnar Pruessner, Rosalba Garcia-Millan

Abstract: Active particles that translate chemical energy into self-propulsion can maintain a far-from-equilibrium steady state and perform work. The entropy production measures how far from equilibrium such a particle system operates and serves as a proxy for the work performed. Field theory offers a promising route to calculating entropy production, as it allows for many interacting particles to be consid… ▽ More Active particles that translate chemical energy into self-propulsion can maintain a far-from-equilibrium steady state and perform work. The entropy production measures how far from equilibrium such a particle system operates and serves as a proxy for the work performed. Field theory offers a promising route to calculating entropy production, as it allows for many interacting particles to be considered simultaneously. Approximate field theories obtained by coarse-graining or smoothing that draw on additive noise can capture densities and correlations well, but they generally ignore the microscopic particle nature of the constituents, thereby producing spurious results for the entropy production. As an alternative we demonstrate how to use Doi-Peliti field theories, which capture the microscopic dynamics, including reactions and interactions with external and pair potentials. Such field theories are in principle exact, while offering a systematic approximation scheme, in the form of diagrammatics. We demonstrate how to construct them from a Fokker-Planck equation (FPE) of the single-particle dynamics and show how to calculate entropy production of active matter from first principles. This framework is easily extended to include interaction. We use it to derive exact, compact and efficient general expressions for the entropy production for a vast range of interacting particle systems. These expressions are independent of the underlying field theory and can be interpreted as the spatial average of the local entropy production. They are readily applicable to numerical and experimental data. In general, any pair interaction draws at most on the three point, equal time density and an n-point interaction on the (2n-1)-point density. We illustrate the technique in a number of exact, tractable examples, including some with pair-interaction. △ Less

Submitted 21 November, 2022; originally announced November 2022.

Comments: 10 page main text, no figures; 49 pages supplement, two figures

Exact solution of a boundary tumbling particle system in one dimension

Authors: Connor Roberts, Gunnar Pruessner

Abstract: We derive the fully time-dependent solution to a run-and-tumble model for a particle which has tumbling restricted to the boundaries of a one-dimensional interval. This is achieved through a field-theoretic perturbative framework by exploiting an elegant underlying structure of the perturbation theory. We calculate the particle densities, currents and variance as well as characteristics of the bou… ▽ More We derive the fully time-dependent solution to a run-and-tumble model for a particle which has tumbling restricted to the boundaries of a one-dimensional interval. This is achieved through a field-theoretic perturbative framework by exploiting an elegant underlying structure of the perturbation theory. We calculate the particle densities, currents and variance as well as characteristics of the boundary tumbling. The analytical findings, in agreement with Monte-Carlo simulations, show how the particle densities are linked to the scale of diffusive fluctuations at the boundaries. The generality of our approach suggests it could be readily applied to similar problems described by Fokker-Planck equations containing localised reaction terms. △ Less

Submitted 13 June, 2022; originally announced June 2022.

Comments: 17 pages, 4 figures

Particle Entity in the Doi-Peliti and Response Field Formalisms

Authors: Marius Bothe, Luca Cocconi, Zigan Zhen, Gunnar Pruessner

Abstract: We introduce a procedure to test a theory for point particle entity, that is, whether said theory takes into account the discrete nature of the constituents of the system. We then identify the mechanism whereby particle entity is enforced in the context of two field-theoretic frameworks designed to incorporate the particle nature of the degrees of freedom, namely the Doi-Peliti field theory and th… ▽ More We introduce a procedure to test a theory for point particle entity, that is, whether said theory takes into account the discrete nature of the constituents of the system. We then identify the mechanism whereby particle entity is enforced in the context of two field-theoretic frameworks designed to incorporate the particle nature of the degrees of freedom, namely the Doi-Peliti field theory and the response field theory that derives from Dean's equation. While the Doi-Peliti field theory encodes the particle nature at a very fundamental level that is easily revealed, demonstrating the same for Dean's equation is more involved and results in a number of surprising diagrammatic identities. We derive those and discuss their implications. These results are particularly pertinent in the context of active matter, whose surprising and often counterintuitive phenomenology rests wholly on the particle nature of the agents and their degrees of freedom as particles. △ Less

Submitted 30 May, 2022; v1 submitted 19 May, 2022; originally announced May 2022.

Optimal Ratchet Potentials for Run-and-Tumble particles

Authors: Zigan Zhen, Gunnar Pruessner

Abstract: Run-and-Tumble particles, mimicking the behaviour of microorganisms like E. coli, are a paradigmatic model of active matter. Due to self-propulsion, their random and undirected motion can be rectified in a ratchet potential. Using perturbative field theory, we determine the shape of the potential that produces the maximum particle current as a function of the particles' parameters. Run-and-Tumble particles, mimicking the behaviour of microorganisms like E. coli, are a paradigmatic model of active matter. Due to self-propulsion, their random and undirected motion can be rectified in a ratchet potential. Using perturbative field theory, we determine the shape of the potential that produces the maximum particle current as a function of the particles' parameters. △ Less

Submitted 8 April, 2022; originally announced April 2022.

Scaling of Entropy Production under Coarse Graining in Active Disordered Media

Authors: Luca Cocconi, Guillaume Salbreux, Gunnar Pruessner

Abstract: Entropy production plays a fundamental role in the study of non-equilibrium systems by offering a quantitative handle on the degree of time-reversal symmetry breaking. It depends crucially on the degree of freedom considered as well as on the scale of description. It was hitherto unknown how the entropy production at one resolution of the degrees of freedom is related to the entropy production at… ▽ More Entropy production plays a fundamental role in the study of non-equilibrium systems by offering a quantitative handle on the degree of time-reversal symmetry breaking. It depends crucially on the degree of freedom considered as well as on the scale of description. It was hitherto unknown how the entropy production at one resolution of the degrees of freedom is related to the entropy production at another resolution. This relationship is of particular relevance to coarse grained and continuum descriptions of a given phenomenon. In this work, we derive the scaling of the entropy production under iterative coarse graining on the basis of the correlations of the underlying microscopic transition rates. Our approach unveils a natural criterion to distinguish equilibrium-like and genuinely non-equilibrium macroscopic phenomena based on the sign of the scaling exponent of the entropy production per mesostate. △ Less

Submitted 18 March, 2022; v1 submitted 11 September, 2021; originally announced September 2021.

Comments: 21 pages, 6 figures

Field theory of survival probabilities, extreme values, first passage times, and mean span of non-Markovian stochastic processes

Authors: Benjamin Walter, Gunnar Pruessner, Guillaume Salbreux

Abstract: We provide a perturbative framework to calculate extreme events of non-Markovian processes, by mapping the stochastic process to a two-species reaction diffusion process in a Doi-Peliti field theory combined with the Martin-Siggia-Rose formalism. This field theory treats interactions and the effect of external, possibly self-correlated noise in a perturbation about a Markovian process, thereby pro… ▽ More We provide a perturbative framework to calculate extreme events of non-Markovian processes, by mapping the stochastic process to a two-species reaction diffusion process in a Doi-Peliti field theory combined with the Martin-Siggia-Rose formalism. This field theory treats interactions and the effect of external, possibly self-correlated noise in a perturbation about a Markovian process, thereby providing a systematic, diagrammatic approach to extreme events. We apply the formalism to Brownian Motion and calculate its survival probability distribution subject to self-correlated noise. △ Less

Submitted 8 April, 2023; v1 submitted 8 September, 2021; originally announced September 2021.

Comments: 24 pages, 4 figures; added figures, comments and references

Journal ref: Phys. Rev. Res. 4 (2022) 043197

Noise can lead to exponential epidemic spreading despite $R_0$ below one

Authors: Johannes Pausch, Rosalba Garcia-Millan, Gunnar Pruessner

Abstract: Branching processes are widely used to model evolutionary and population dynamics as well as the spread of infectious diseases. To characterize the dynamics of their growth or spread, the basic reproduction number $R_0$ has received considerable attention. In the context of infectious diseases, it is usually defined as the expected number of secondary cases produced by an infectious case in a comp… ▽ More Branching processes are widely used to model evolutionary and population dynamics as well as the spread of infectious diseases. To characterize the dynamics of their growth or spread, the basic reproduction number $R_0$ has received considerable attention. In the context of infectious diseases, it is usually defined as the expected number of secondary cases produced by an infectious case in a completely susceptible population. Typically $R_0>1$ indicates that an outbreak is expected to continue and to grow exponentially, while $R_0<1$ usually indicates that an outbreak is expected to terminate after some time. In this work, we show that fluctuations of the dynamics in time can lead to a continuation of outbreaks even when the expected number of secondary cases from a single case is below $1$. Such fluctuations are usually neglected in modelling of infectious diseases by a set of ordinary differential equations, such as the classic SIR model. We showcase three examples: 1) extinction following an Ornstein-Uhlenbeck process, 2) extinction switching randomly between two values and 3) mixing of two populations with different $R_0$ values. We corroborate our analytical findings with computer simulations. △ Less

Submitted 1 September, 2021; originally announced September 2021.

Comments: 17 pages, 10 figures

Field Theory of Free Run and Tumble Particles in d Dimensions

Authors: Ziluo Zhang, Gunnar Pruessner

Abstract: In this paper, Doi-Peliti field theory is used to describe the motion of free Run and Tumble particles in arbitrary dimensions. After deriving action and propagators, the mean square displacement and the corresponding entropy production at stationarity are calculated in this framework. We further derive the field theory of free Active Brownian Particles in two dimensions for comparison. In this paper, Doi-Peliti field theory is used to describe the motion of free Run and Tumble particles in arbitrary dimensions. After deriving action and propagators, the mean square displacement and the corresponding entropy production at stationarity are calculated in this framework. We further derive the field theory of free Active Brownian Particles in two dimensions for comparison. △ Less

Submitted 11 June, 2021; originally announced June 2021.

Field Theory of free Active Ornstein-Uhlenbeck Particles

Authors: Marius Bothe, Gunnar Pruessner

Abstract: We derive a Doi-Peliti field theory for free active Ornstein-Uhlenbeck particles, or, equivalently, free inertial Brownian particles, and present a way to diagonalise the Gaussian part of the action and calculate the propagator. Unlike previous course-grained approaches this formulation correctly tracks particle identity and yet can easily be expanded to include potentials and arbitrary reactions. We derive a Doi-Peliti field theory for free active Ornstein-Uhlenbeck particles, or, equivalently, free inertial Brownian particles, and present a way to diagonalise the Gaussian part of the action and calculate the propagator. Unlike previous course-grained approaches this formulation correctly tracks particle identity and yet can easily be expanded to include potentials and arbitrary reactions. △ Less

Submitted 22 January, 2021; v1 submitted 18 January, 2021; originally announced January 2021.

Journal ref: Phys. Rev. E 103, 062105 (2021)

Run-and-tumble motion: field theory and entropy production

Authors: Rosalba Garcia-Millan, Gunnar Pruessner

Abstract: Run-and-tumble motion is an example of active motility where particles move at constant speed and change direction at random times. In this work we study run-and-tumble motion with diffusion in a harmonic potential in one dimension via a path integral approach. We derive a Doi-Peliti field theory and use it to calculate the entropy production and other observables in closed form. All our results a… ▽ More Run-and-tumble motion is an example of active motility where particles move at constant speed and change direction at random times. In this work we study run-and-tumble motion with diffusion in a harmonic potential in one dimension via a path integral approach. We derive a Doi-Peliti field theory and use it to calculate the entropy production and other observables in closed form. All our results are exact. △ Less

Submitted 4 December, 2020; originally announced December 2020.

Comments: 22 pages, 6 figures

Entropy production in exactly solvable systems

Authors: Luca Cocconi, Rosalba Garcia-Millan, Zigan Zhen, Bianca Buturca, Gunnar Pruessner

Abstract: The rate of entropy production by a stochastic process quantifies how far it is from thermodynamic equilibrium. Equivalently, entropy production captures the degree to which detailed balance and time-reversal symmetry are broken. Despite abundant references to entropy production in the literature and its many applications in the study of non-equilibrium stochastic particle systems, a comprehensive… ▽ More The rate of entropy production by a stochastic process quantifies how far it is from thermodynamic equilibrium. Equivalently, entropy production captures the degree to which detailed balance and time-reversal symmetry are broken. Despite abundant references to entropy production in the literature and its many applications in the study of non-equilibrium stochastic particle systems, a comprehensive list of typical examples illustrating the fundamentals of entropy production is lacking. Here, we present a brief, self-contained review of entropy production and calculate it from first principles in a catalogue of exactly solvable setups, encompassing both discrete- and continuous-state Markov processes, as well as single- and multiple-particle systems. The examples covered in this work provide a stepping stone for further studies on entropy production of more complex systems, such as many-particle active matter, as well as a benchmark for the development of alternative mathematical formalisms. △ Less

Submitted 8 October, 2020; originally announced October 2020.

Comments: 33 pages, 13 figures

First passage time distribution of active thermal particles in potentials

Authors: Benjamin Walter, Gunnar Pruessner, Guillaume Salbreux

Abstract: We introduce a perturbative method to calculate all moments of the first-passage time distribution in stochastic one-dimensional processes which are subject to both white and coloured noise. This class of non-Markovian processes is at the centre of the study of thermal active matter, that is self-propelled particles subject to diffusion. The perturbation theory about the Markov process considers t… ▽ More We introduce a perturbative method to calculate all moments of the first-passage time distribution in stochastic one-dimensional processes which are subject to both white and coloured noise. This class of non-Markovian processes is at the centre of the study of thermal active matter, that is self-propelled particles subject to diffusion. The perturbation theory about the Markov process considers the effect of self-propulsion to be small compared to that of thermal fluctuations. To illustrate our method, we apply it to the case of active thermal particles (i) in a harmonic trap (ii) on a ring. For both we calculate the first-order correction of the moment-generating function of first-passage times, and thus to all its moments. Our analytical results are compared to numerics. △ Less

Submitted 11 February, 2021; v1 submitted 29 May, 2020; originally announced June 2020.

Comments: 25 pages, 10 figures

Journal ref: Phys. Rev. Research 3, 013075 (2021)

Dynamically accelerated cover times

Authors: Gcina Maziya, Luca Cocconi, Gunnar Pruessner, Nicholas Moloney

Abstract: Among observables characterising the random exploration of a graph or lattice, the cover time, namely the time to visit every site, continues to attract widespread interest. Much insight about cover times is gained by mapping to the (spaceless) coupon-collector problem, which amounts to ignoring spatio-temporal correlations, and an early conjecture that the limiting cover time distribution of regu… ▽ More Among observables characterising the random exploration of a graph or lattice, the cover time, namely the time to visit every site, continues to attract widespread interest. Much insight about cover times is gained by mapping to the (spaceless) coupon-collector problem, which amounts to ignoring spatio-temporal correlations, and an early conjecture that the limiting cover time distribution of regular random walks on large lattices converges to the Gumbel distribution in $d \geq 3$ was recently proved rigorously. Furthermore, a number of mathematical and numerical studies point to the robustness of the Gumbel universality to modifications of the \textit{spatial} features of the random search processes (e.g.\ introducing persistence and/or intermittence, or changing the graph topology). Here we investigate the robustness of the Gumbel universality to dynamical modification of the \textit{temporal} features of the search, specifically by allowing the random walker to "accelerate" or "decelerate" upon visiting a previously unexplored site. We generalise the mapping mentioned above by relating the statistics of cover times to the roughness of $1/f^αあるふぁ$ Gaussian signals, leading to the conjecture that the Gumbel distribution is but one of a family of cover time distributions, ranging from Gaussian for highly accelerated cover, to exponential for highly decelerated cover. While our conjecture is confirmed by systematic Monte Carlo simulations in dimensions $d > 3$, our results for acceleration in $d=3$ challenge the current understanding of the role of correlations in the cover time problem. △ Less

Submitted 6 July, 2020; v1 submitted 27 April, 2020; originally announced April 2020.

Comments: 10 pages, 15 figures

Journal ref: Phys. Rev. Research 2, 023421 (2020)

Time-dependent branching processes: a model of oscillating neuronal avalanches

Authors: Johannes Pausch, Rosalba Garcia-Millan, Gunnar Pruessner

Abstract: Recently, neuronal avalanches have been observed to display oscillations, a phenomenon regarded as the co-existence of a scale-free behaviour (the avalanches close to criticality) and scale-dependent dynamics (the oscillations). Ordinary continuous-time branching processes with constant extinction and branching rates are commonly used as models of neuronal activity, yet they lack any such time-dep… ▽ More Recently, neuronal avalanches have been observed to display oscillations, a phenomenon regarded as the co-existence of a scale-free behaviour (the avalanches close to criticality) and scale-dependent dynamics (the oscillations). Ordinary continuous-time branching processes with constant extinction and branching rates are commonly used as models of neuronal activity, yet they lack any such time-dependence. In the present work, we extend a basic branching process by allowing the extinction rate to oscillate in time as a new model to describe cortical dynamics. By means of a perturbative field theory, we derive relevant observables in closed form. We support our findings by quantative comparison to numerics and qualitative comparison to available experimental results. △ Less

Submitted 6 March, 2020; originally announced March 2020.

Volume explored by a branching random walk on general graphs

Authors: Ignacio Bordeu, Saoirse Amarteifio, Rosalba Garcia-Millan, Benjamin Walter, Nanxin Wei, Gunnar Pruessner

Abstract: Branching processes are used to model diverse social and physical scenarios, from extinction of family names to nuclear fission. However, for a better description of natural phenomena, such as viral epidemics in cellular tissues, animal populations and social networks, a spatial embedding---the branching random walk (BRW)---is required. Despite its wide range of applications, the properties of the… ▽ More Branching processes are used to model diverse social and physical scenarios, from extinction of family names to nuclear fission. However, for a better description of natural phenomena, such as viral epidemics in cellular tissues, animal populations and social networks, a spatial embedding---the branching random walk (BRW)---is required. Despite its wide range of applications, the properties of the volume explored by the BRW so far remained elusive, with exact results limited to one dimension. Here we present analytical results, supported by numerical simulations, on the scaling of the volume explored by a BRW in the critical regime, the onset of epidemics, in general environments. Our results characterise the spreading dynamics on regular lattices and general graphs, such as fractals, random trees and scale-free networks, revealing the direct relation between the graphs' dimensionality and the rate of propagation of the viral process. Furthermore, we use the BRW to determine the spectral properties of real social and metabolic networks, where we observe that a lack of information of the network structure can lead to differences in the observed behaviour of the spreading process. Our results provide observables of broad interest for the characterisation of real world lattices, tissues, and networks. △ Less

Submitted 1 September, 2019; originally announced September 2019.

Comments: 11 pages, 6 figures

Critical Density of the Abelian Manna Model via a Multi-type Branching Process

Authors: Nanxin Wei, Gunnar Pruessner

Abstract: A multi-type branching process is introduced to mimic the evolution of the avalanche activity and determine the critical density of the Abelian Manna model. This branching process incorporates partially the spatio-temporal correlations of the activity, which are essential for the dynamics, in particular in low dimensions. An analytical expression for the critical density in arbitrary dimensions is… ▽ More A multi-type branching process is introduced to mimic the evolution of the avalanche activity and determine the critical density of the Abelian Manna model. This branching process incorporates partially the spatio-temporal correlations of the activity, which are essential for the dynamics, in particular in low dimensions. An analytical expression for the critical density in arbitrary dimensions is derived, which significantly improves the results over mean-field theories, as confirmed by comparison to the literature on numerical estimates from simulations. The method can easily be extended to lattices and dynamics other than those studied in the present work. △ Less

Submitted 25 September, 2019; v1 submitted 20 February, 2019; originally announced February 2019.

Comments: 5 pages, 2 figures, 2 tables

Journal ref: Phys. Rev. E 100, 032116 (2019)

Uncovering novel phase transitions in dense dry polar active fluids using a lattice Boltzmann method

Authors: David Nesbitt, Gunnar Pruessner, Chiu Fan Lee

Abstract: The dynamics of dry active matter have implications for a diverse collection of biological phenomena spanning a range of length and time scales, such as animal flocking, cell tissue dynamics, and swarming of inserts and bacteria. Uniting these systems are a common set of symmetries and conservation laws, defining dry active fluids as a class of physical system. Many interesting behaviours have bee… ▽ More The dynamics of dry active matter have implications for a diverse collection of biological phenomena spanning a range of length and time scales, such as animal flocking, cell tissue dynamics, and swarming of inserts and bacteria. Uniting these systems are a common set of symmetries and conservation laws, defining dry active fluids as a class of physical system. Many interesting behaviours have been observed at high densities, which remain difficult to simulate due to the computational demand. Here, we show how two-dimensional dry active fluids in a dense regime can be studied using a simple modification of the lattice Boltzmann method. We apply our method on a model that exhibits motility-induced phase separation, and an active model with contact inhibition of locomotion, which has relevance to collective cell migration. For the latter, we uncover multiple novel phase transitions: two first-order and one potentially critical. We further support our simulation results with an analytical treatment of the hydrodynamic equations obtained via a Chapman-Enskog coarse-graining procedure. △ Less

Submitted 4 January, 2021; v1 submitted 1 February, 2019; originally announced February 2019.

Comments: 25 pages, 4 figures

Journal ref: New Journal of Physics (2021)

Is Actin Filament and Microtubule Growth Reaction- or Diffusion-Limited?

Authors: Johannes Pausch, Gunnar Pruessner

Abstract: Inside cells of living organisms, actin filaments and microtubules self-assemble and dissemble dynamically by incorporating actin or tubulin from the cell plasma or releasing it into their tips' surroundings. Such reaction-diffusion systems can show diffusion- or reaction-limited behaviour. However, neither limit explains the experimental data: while the offset of the linear relation between growt… ▽ More Inside cells of living organisms, actin filaments and microtubules self-assemble and dissemble dynamically by incorporating actin or tubulin from the cell plasma or releasing it into their tips' surroundings. Such reaction-diffusion systems can show diffusion- or reaction-limited behaviour. However, neither limit explains the experimental data: while the offset of the linear relation between growth speed and bulk tubulin density contradicts the diffusion limit, the surprisingly large variance of the growth speed rejects a pure reaction limit. In this Letter, we accommodate both limits and use a Doi-Peliti field-theory model to estimate how diffusive transport is perturbing the chemical reactions at the filament tip. Furthermore, a crossover bulk density is predicted at which the limiting process changes from chemical reactions to diffusive transport. In addition, we explain and estimate larger variances of the growth speed. △ Less

Submitted 29 August, 2018; originally announced August 2018.

MSC Class: 92Bxx

Field-theoretic approach to the universality of branching processes

Authors: Rosalba Garcia-Millan, Johannes Pausch, Benjamin Walter, Gunnar Pruessner

Abstract: Branching processes are widely used to model phenomena from networks to neuronal avalanching. In a large class of continuous-time branching processes, we study the temporal scaling of the moments of the instant population size, the survival probability, expected avalanche duration, the so-called avalanche shape, the $n$-point correlation function and the probability density function of the total a… ▽ More Branching processes are widely used to model phenomena from networks to neuronal avalanching. In a large class of continuous-time branching processes, we study the temporal scaling of the moments of the instant population size, the survival probability, expected avalanche duration, the so-called avalanche shape, the $n$-point correlation function and the probability density function of the total avalanche size. Previous studies have shown universality in certain observables of branching processes using probabilistic arguments, however, a comprehensive description is lacking. We derive the field theory that describes the process and demonstrate how to use it to calculate the relevant observables and their scaling to leading order in time, revealing the universality of the moments of the population size. Our results explain why the first and second moment of the offspring distribution are sufficient to fully characterise the process in the vicinity of criticality, regardless of the underlying offspring distribution. This finding implies that branching processes are universal. We illustrate our analytical results with computer simulations. △ Less

Submitted 15 December, 2018; v1 submitted 25 August, 2018; originally announced August 2018.

Comments: 12 pages, 8 figures

Journal ref: Phys. Rev. E 98, 062107 (2018)

Linear stability analysis of morphodynamics during tissue regeneration in plants

Authors: Anne-Mieke Reijne, Gunnar Pruessner, Giovanni Sena

Abstract: One of the key characteristics of multicellular organisms is the ability to establish and maintain shapes, or morphologies, under a variety of physical and chemical perturbations. A quantitative description of the underlying morphological dynamics is a critical step to fully understand the self-organising properties of multicellular systems. Although many powerful mathematical tools have been deve… ▽ More One of the key characteristics of multicellular organisms is the ability to establish and maintain shapes, or morphologies, under a variety of physical and chemical perturbations. A quantitative description of the underlying morphological dynamics is a critical step to fully understand the self-organising properties of multicellular systems. Although many powerful mathematical tools have been developed to analyse stochastic dynamics, rarely these are applied to experimental developmental biology. Here, we take root tip regeneration in the plant model system Arabidopsis thaliana as an example of robust morphogenesis in living tissue and present a novel approach to quantify and model the relaxation of the system to its unperturbed morphology. By generating and analysing time-lapse series of regenerating root tips captured with confocal microscopy, we are able to extract and model the dynamics of key morphological traits at cellular resolution. We present a linear stability analysis of its Markovian dynamics, with the stationary state representing the intact root in the space of morphological traits. We find that the resulting eigenvalues can be classified into two groups, suggesting the co-existence of two distinct temporal scales during the process of regeneration. We discuss the possible biological implications of our specific results and suggest future experiments to further probe the self-organising properties of living tissue. △ Less

Submitted 21 June, 2018; v1 submitted 17 June, 2018; originally announced June 2018.

Comments: 18 pages, 11 figures, 1 table. Typos fixed, added one ref, improved some figure labels

Correlations and hyperuniformity in the avalanche size of the Oslo Model

Authors: Rosalba Garcia-Millan, Gunnar Pruessner, Luke Pickering, Kim Christensen

Abstract: Certain random processes display anticorrelations resulting in local Poisson-like disorder and global order, where correlations suppress fluctuations. Such processes are called hyperuniform. Using a map to an interface picture we show via analytic calculations that a sequence of avalanche sizes of the Oslo Model is hyperuniform in the temporal domain with the minimal exponent $λらむだ=0$. We identify th… ▽ More Certain random processes display anticorrelations resulting in local Poisson-like disorder and global order, where correlations suppress fluctuations. Such processes are called hyperuniform. Using a map to an interface picture we show via analytic calculations that a sequence of avalanche sizes of the Oslo Model is hyperuniform in the temporal domain with the minimal exponent $λらむだ=0$. We identify the conserved quantity in the interface picture that gives rise to the hyperuniformity in the avalanche size. We further discuss the fluctuations of the avalanche size in two variants of the Oslo Model. We support our findings with numerical results. △ Less

Submitted 5 August, 2018; v1 submitted 30 September, 2017; originally announced October 2017.

Comments: 11 pages, 7 figures

MSC Class: 60G99

Journal ref: EPL (Europhysics Letters), 122(5), 50003 (2018)

Edge instability in incompressible planar active fluids

Authors: David Nesbitt, Gunnar Pruessner, Chiu Fan Lee

Abstract: Interfacial instability is highly relevant to many important biological processes. A key example arises in wound healing experiments, which observe that an epithelial layer with an initially straight edge does not heal uniformly. We consider the phenomenon in the context of active fluids. Improving upon the approximation used in J. Zimmermann, M. Basan and H. Levine, Euro. Phys. J.: Special Topics… ▽ More Interfacial instability is highly relevant to many important biological processes. A key example arises in wound healing experiments, which observe that an epithelial layer with an initially straight edge does not heal uniformly. We consider the phenomenon in the context of active fluids. Improving upon the approximation used in J. Zimmermann, M. Basan and H. Levine, Euro. Phys. J.: Special Topics 223, 1259 (2014), we perform a linear stability analysis on a two dimensional incompressible hydrodynamic model of an active fluid with an open interface. We categorise the stability of the model and find that for experimentally relevant parameters, fingering stability is always absent in this minimal model. Our results point to the crucial importance of density variation in the fingering instability in tissue regeneration. △ Less

Submitted 12 June, 2017; originally announced June 2017.

Comments: 9 pages, 4 figures

Journal ref: Phys. Rev. E 96, 062615 (2017)

Statistical mechanics of exploding phase spaces: Ontic open systems

Authors: Henrik Jeldtoft Jensen, Roozbeh H. Pazuki, Gunnar Pruessner, Piergiulio Tempesta

Abstract: The volume of phase space may grow super-exponentially ("explosively") with the number of degrees of freedom for certain types of complex systems such as those encountered in biology and neuroscience, where components interact and create new emergent states. Standard ensemble theory can break down as we demonstrate in a simple model reminiscent of complex systems where new collective states emerge… ▽ More The volume of phase space may grow super-exponentially ("explosively") with the number of degrees of freedom for certain types of complex systems such as those encountered in biology and neuroscience, where components interact and create new emergent states. Standard ensemble theory can break down as we demonstrate in a simple model reminiscent of complex systems where new collective states emerge. We present an axiomatically defined entropy and argue that it is extensive in the micro-canonical, equal probability, and canonical (max-entropy) ensemble for super-exponentially growing phase spaces. This entropy may be useful in determining probability measures in analogy with how statistical mechanics establishes statistical ensembles by maximising entropy. △ Less

Submitted 25 July, 2018; v1 submitted 7 September, 2016; originally announced September 2016.

Comments: 18 pages, 1 figure. This version is significantly expanded and improved as a result of referee comments. This is the manuscript accepted for publication in J Phys A

Spatio-temporal Correlations in the Manna Model in one, three and five dimensions

Authors: Gary Willis, Gunnar Pruessner

Abstract: Although the paradigm of criticality is centred around spatial correlations and their anomalous scaling, not many studies of Self-Organised Criticality (SOC) focus on spatial correlations. Often, integrated observables, such as avalanche size and duration, are used, not least as to avoid complications due to the unavoidable lack of translational invariance. The present work is a survey of spatio-t… ▽ More Although the paradigm of criticality is centred around spatial correlations and their anomalous scaling, not many studies of Self-Organised Criticality (SOC) focus on spatial correlations. Often, integrated observables, such as avalanche size and duration, are used, not least as to avoid complications due to the unavoidable lack of translational invariance. The present work is a survey of spatio-temporal correlation functions in the Manna Model of SOC, measured numerically in detail in d=1,3 and 5 dimensions and compared to theoretical results, in particular relating them to "integrated" observables such as avalanche size and duration scaling, that measure them indirectly. Contrary to the notion held by some of SOC models organising into a critical state by re-arranging their spatial structure avalanche by avalanche, which may be expected to result in large, non-trivial, system-spanning spatial correlations in the quiescent state (between avalanches), correlations of inactive particles in the quiescent state have a small amplitude that does not increase with the system size, although they display (noisy) power law scaling over a range linear in the system size. Self-organisation, however, does take place as the (one-point) density of inactive particles organises into a particular profile that is asymptotically independent of the driving location, also demonstrated analytically in one dimension. Activity and its correlations, on the other hand, display non-trivial long-ranged spatio-temporal scaling with exponents that can be related to established results, in particular avalanche size and duration exponents. The correlation length and amplitude are set by the system size (confirmed analytically for some observables), as expected in systems displaying finite size scaling. In one dimension, we find some surprising inconsistencies of the dynamical exponent. A (spatially extended) mean ... △ Less

Submitted 30 July, 2017; v1 submitted 2 August, 2016; originally announced August 2016.

Comments: 42 pages, 28 figures (added references and some more data and discussion)

Comment on "Anomalous Discontinuity at the Percolation Critical Point of Active Gels"

Authors: Gunnar Pruessner, Chiu Fan Lee

Abstract: In their recent work Sheinman et al. [Phys. Rev. Lett. 114, 098104 (2015)] introduce a variation of percolation which they call no-enclaves percolation (NEP). The main claims are 1) the salient physics captured in NEP is closer to what happens experimentally; 2) The Fisher exponent of NEP, is $τたう=1.82(1)$; 3) Due to the different Fisher exponent, NEP constitutes a universality class distinct from r… ▽ More In their recent work Sheinman et al. [Phys. Rev. Lett. 114, 098104 (2015)] introduce a variation of percolation which they call no-enclaves percolation (NEP). The main claims are 1) the salient physics captured in NEP is closer to what happens experimentally; 2) The Fisher exponent of NEP, is $τたう=1.82(1)$; 3) Due to the different Fisher exponent, NEP constitutes a universality class distinct from random percolation (RP). While we fully agree with 1) and found NEP to be a very interesting variation of random percolation, we disagree with 2) and 3). We will demonstrate that $τたう$ is exactly $2$, directly derivable from RP, and thus there is no foundation of a new universality class. △ Less

Submitted 10 May, 2016; originally announced May 2016.

Comments: Comment on Sheinman et al. [Phys. Rev. Lett. 114, 098104 (2015)], 2 pages, 1 figure, reply to appear here as well

Journal ref: Physical Review Letters 116, 189801 (2016)

Directed Abelian sandpile with multiple downward neighbors

Authors: Deepak Dhar, Gunnar Pruessner, Paul Expert, Kim Christensen, Nicky Zachariou

Abstract: We study the directed Abelian sandpile model on a square lattice, with $K$ downward neighbors per site, $K > 2$. The $K=3$ case is solved exactly, which extends the earlier known solution for the $K=2$ case. For $K>2$, the avalanche clusters can have holes and side-branches and are thus qualitatively different from the $K=2$ case where avalance clusters are compact. However, we find that the criti… ▽ More We study the directed Abelian sandpile model on a square lattice, with $K$ downward neighbors per site, $K > 2$. The $K=3$ case is solved exactly, which extends the earlier known solution for the $K=2$ case. For $K>2$, the avalanche clusters can have holes and side-branches and are thus qualitatively different from the $K=2$ case where avalance clusters are compact. However, we find that the critical exponents for $K>2$ are identical with those for the $K=2$ case, and the large scale structure of the avalanches for $K>2$ tends to the $K=2$ case. △ Less

Submitted 11 February, 2016; v1 submitted 19 November, 2015; originally announced November 2015.

Comments: Accepted for publication in PRE

Journal ref: Phys. Rev. E 93, 042107 (2016)

Percolation with trapping mechanism drives active gels to the critically connected state

Authors: Chiu Fan Lee, Gunnar Pruessner

Abstract: Cell motility and tissue morphogenesis depend crucially on the dynamic remodelling of actomyosin networks. An actomyosin network consists of an actin polymer network connected by crosslinker proteins and motor protein myosins that generate internal stresses on the network. A recent discovery shows that for a range of experimental parameters, actomyosin networks contract to clusters with a power-la… ▽ More Cell motility and tissue morphogenesis depend crucially on the dynamic remodelling of actomyosin networks. An actomyosin network consists of an actin polymer network connected by crosslinker proteins and motor protein myosins that generate internal stresses on the network. A recent discovery shows that for a range of experimental parameters, actomyosin networks contract to clusters with a power-law size distribution [Alvarado J. et al. (2013) Nature Physics 9 591]. Here, we argue that actomyosin networks can exhibit robust critical signature without fine-tuning because the dynamics of the system can be mapped onto a modified version of percolation with trapping (PT), which is known to show critical behaviour belonging to the static percolation universality class without the need of fine-tuning of a control parameter. We further employ our PT model to generate experimentally testable predictions. △ Less

Submitted 27 April, 2016; v1 submitted 27 August, 2015; originally announced August 2015.

Comments: 7 pages, 6 figures. To appear in Physical Review E

Journal ref: Phys. Rev. E 93, 052414 (2016)

25 Years of Self-Organized Criticality: Numerical Detection Methods

Authors: R. T. James McAteer, Markus J. Aschwanden, Michaila Dimitropoulou, Manolis K. Georgoulis, Gunnar Pruessner, Laura Morales, Jack Ireland, Valentyna Abramenko

Abstract: The detection and characterization of self-organized criticality (SOC), in both real and simulated data, has undergone many significant revisions over the past 25 years. The explosive advances in the many numerical methods available for detecting, discriminating, and ultimately testing, SOC have played a critical role in developing our understanding of how systems experience and exhibit SOC. In th… ▽ More The detection and characterization of self-organized criticality (SOC), in both real and simulated data, has undergone many significant revisions over the past 25 years. The explosive advances in the many numerical methods available for detecting, discriminating, and ultimately testing, SOC have played a critical role in developing our understanding of how systems experience and exhibit SOC. In this article, methods of detecting SOC are reviewed; from correlations to complexity to critical quantities. A description of the basic autocorrelation method leads into a detailed analysis of application-oriented methods developed in the last 25 years. In the second half of this manuscript space-based, time-based and spatial-temporal methods are reviewed and the prevalence of power laws in nature is described, with an emphasis on event detection and characterization. The search for numerical methods to clearly and unambiguously detect SOC in data often leads us outside the comfort zone of our own disciplines - the answers to these questions are often obtained by studying the advances made in other fields of study. In addition, numerical detection methods often provide the optimum link between simulations and experiments in scientific research. We seek to explore this boundary where the rubber meets the road, to review this expanding field of research of numerical detection of SOC systems over the past 25 years, and to iterate forwards so as to provide some foresight and guidance into developing breakthroughs in this subject over the next quarter of a century. △ Less

Submitted 26 June, 2015; originally announced June 2015.

Comments: Space Science Review series on SOC

25 Years of Self-Organized Criticality: Concepts and Controversies

Authors: Nicholas Watkins, Gunnar Pruessner, Sandra Chapman, Norma Bock Crosby, Henrik Jensen

Abstract: Introduced by the late Per Bak and his colleagues, self-organized criticality (SOC) has been one of the most stimulating concepts to come out of statistical mechanics and condensed matter theory in the last few decades, and has played a significant role in the development of complexity science. SOC, and more generally fractals and power laws, have attacted much comment, ranging from the very posit… ▽ More Introduced by the late Per Bak and his colleagues, self-organized criticality (SOC) has been one of the most stimulating concepts to come out of statistical mechanics and condensed matter theory in the last few decades, and has played a significant role in the development of complexity science. SOC, and more generally fractals and power laws, have attacted much comment, ranging from the very positive to the polemical. The other papers in this special issue (Aschwanden et al, 2014; McAteer et al, 2014; Sharma et al, 2015) showcase the considerable body of observations in solar, magnetospheric and fusion plasma inspired by the SOC idea, and expose the fertile role the new paradigm has played in approaches to modeling and understanding multiscale plasma instabilities. This very broad impact, and the necessary process of adapting a scientific hypothesis to the conditions of a given physical system, has meant that SOC as studied in these fields has sometimes differed significantly from the definition originally given by its creators. In Bak's own field of theoretical physics there are significant observational and theoretical open questions, even 25 years on (Pruessner, 2012). One aim of the present review is to address the dichotomy between the great reception SOC has received in some areas, and its shortcomings, as they became manifest in the controversies it triggered. Our article tries to clear up what we think are misunderstandings of SOC in fields more remote from its origins in statistical mechanics, condensed matter and dynamical systems by revisiting Bak, Tang and Wiesenfeld's original papers. △ Less

Submitted 20 April, 2015; originally announced April 2015.

Comments: 52 pages, 3 figures

Journal ref: Space.Sci.Rev. 198 (2016) 3-44

The perils of thresholding

Authors: Francesc Font-Clos, Gunnar Pruessner, Anna Deluca, Nicholas R. Moloney

Abstract: The thresholding of time series of activity or intensity is frequently used to define and differentiate events. This is either implicit, for example due to resolution limits, or explicit, in order to filter certain small scale physics from the supposed true asymptotic events. Thresholding the birth-death process, however, introduces a scaling region into the event size distribution, which is chara… ▽ More The thresholding of time series of activity or intensity is frequently used to define and differentiate events. This is either implicit, for example due to resolution limits, or explicit, in order to filter certain small scale physics from the supposed true asymptotic events. Thresholding the birth-death process, however, introduces a scaling region into the event size distribution, which is characterised by an exponent that is unrelated to the actual asymptote and is rather an artefact of thresholding. As a result, numerical fits of simulation data produce a range of exponents, with the true asymptote visible only in the tail of the distribution. This tail is increasingly difficult to sample as the threshold is increased. In the present case, the exponents and the spurious nature of the scaling region can be determined analytically, thus demonstrating the way in which thresholding conceals the true asymptote. The analysis also suggests a procedure for detecting the influence of the threshold by means of a data collapse involving the threshold-imposed scale. △ Less

Submitted 22 October, 2014; originally announced October 2014.

Comments: 16 pages, 10 figures

A field-theoretic approach to the Wiener Sausage

Authors: Stefan Nekovar, Gunnar Pruessner

Abstract: The Wiener Sausage, the volume traced out by a sphere attached to a Brownian particle, is a classical problem in statistics and mathematical physics. Initially motivated by a range of field-theoretic, technical questions, we present a single loop renormalised perturbation theory of a stochastic process closely related to the Wiener Sausage, which, however, proves to be exact for the exponents and… ▽ More The Wiener Sausage, the volume traced out by a sphere attached to a Brownian particle, is a classical problem in statistics and mathematical physics. Initially motivated by a range of field-theoretic, technical questions, we present a single loop renormalised perturbation theory of a stochastic process closely related to the Wiener Sausage, which, however, proves to be exact for the exponents and some amplitudes. The field-theoretic approach is particularly elegant and very enjoyable to see at work on such a classic problem. While we recover a number of known, classical results, the field-theoretic techniques deployed provide a particularly versatile framework, which allows easy calculation with different boundary conditions even of higher momenta and more complicated correlation functions. At the same time, we provide a highly instructive, non-trivial example for some of the technical particularities of the field-theoretic description of stochastic processes, such as excluded volume, lack of translational invariance and immobile particles. The aim of the present work is not to improve upon the well-established results for the Wiener Sausage, but to provide a field-theoretic approach to it, in order to gain a better understanding of the field-theoretic obstacles to overcome. △ Less

Submitted 29 February, 2016; v1 submitted 28 July, 2014; originally announced July 2014.

Comments: 45 pages, 3 Figures, Springer style

Journal ref: J Stat Phys (2016) 163: 604

Opinion formation models on a gradient

Authors: Michael T. Gastner, Nikolitsa Markou, Gunnar Pruessner, Moez Draief

Abstract: Statistical physicists have become interested in models of collective social behavior such as opinion formation, where individuals change their inherently preferred opinion if their friends disagree. Real preferences often depend on regional cultural differences, which we model here as a spatial gradient $g$ in the initial opinion. The gradient does not only add reality to the model. It can also r… ▽ More Statistical physicists have become interested in models of collective social behavior such as opinion formation, where individuals change their inherently preferred opinion if their friends disagree. Real preferences often depend on regional cultural differences, which we model here as a spatial gradient $g$ in the initial opinion. The gradient does not only add reality to the model. It can also reveal that opinion clusters in two dimensions are typically in the standard (i.e.\ independent) percolation universality class, thus settling a recent controversy about a non-consensus model. However, using analytical and numerical tools, we also present a model where the width of the transition between opinions scales $\propto g^{-1/4}$, not $\propto g^{-4/7}$ as in independent percolation, and the cluster size distribution is consistent with first-order percolation. △ Less

Submitted 26 January, 2015; v1 submitted 6 June, 2014; originally announced June 2014.

Comments: 12 pages, 8 figures, version accepted by PLoS ONE, online supplement added as appendix

Journal ref: PLoS ONE 9(12): e114088 (2014)

25 Years of Self-Organized Criticality: Solar and Astrophysics

Authors: Markus J. Aschwanden, Norma Crosby, Michaila Dimitropoulou, Manolis Georgoulis, Stefan Hergarten, James MdAteer, Alexander V. Milovanov, Shin Mineshige, Laura Morales, Naoto Nishizuka, Gunnar Pruessner, Raul Sanchez, Surja Sharma, Antoine Strugarek, Vadim Uritsky

Abstract: Shortly after the seminal paper {\sl "Self-Organized Criticality: An explanation of 1/f noise"} by Bak, Tang, and Wiesenfeld (1987), the idea has been applied to solar physics, in {\sl "Avalanches and the Distribution of Solar Flares"} by Lu and Hamilton (1991). In the following years, an inspiring cross-fertilization from complexity theory to solar and astrophysics took place, where the SOC conce… ▽ More Shortly after the seminal paper {\sl "Self-Organized Criticality: An explanation of 1/f noise"} by Bak, Tang, and Wiesenfeld (1987), the idea has been applied to solar physics, in {\sl "Avalanches and the Distribution of Solar Flares"} by Lu and Hamilton (1991). In the following years, an inspiring cross-fertilization from complexity theory to solar and astrophysics took place, where the SOC concept was initially applied to solar flares, stellar flares, and magnetospheric substorms, and later extended to the radiation belt, the heliosphere, lunar craters, the asteroid belt, the Saturn ring, pulsar glitches, soft X-ray repeaters, blazars, black-hole objects, cosmic rays, and boson clouds. The application of SOC concepts has been performed by numerical cellular automaton simulations, by analytical calculations of statistical (powerlaw-like) distributions based on physical scaling laws, and by observational tests of theoretically predicted size distributions and waiting time distributions. Attempts have been undertaken to import physical models into the numerical SOC toy models, such as the discretization of magneto-hydrodynamics (MHD) processes. The novel applications stimulated also vigorous debates about the discrimination between SOC models, SOC-like, and non-SOC processes, such as phase transitions, turbulence, random-walk diffusion, percolation, branching processes, network theory, chaos theory, fractality, multi-scale, and other complexity phenomena. We review SOC studies from the last 25 years and highlight new trends, open questions, and future challenges, as discussed during two recent ISSI workshops on this theme. △ Less

Submitted 25 March, 2014; originally announced March 2014.

Comments: 139 pages, 28 figures, Review based on ISSI workshops "Self-Organized Criticality and Turbulence" (2012, 2013, Bern, Switzerland)

Mapping multiplicative to additive noise

Authors: Katy J. Rubin, Gunnar Pruessner, Grigorios A. Pavliotis

Abstract: The Langevin formulation of a number of well-known stochastic processes involves multiplicative noise. In this work we present a systematic mapping of a process with multiplicative noise to a related process with additive noise, which may often be easier to analyse. The mapping is easily understood in the example of the branching process. In a second example we study the random neighbour (or infin… ▽ More The Langevin formulation of a number of well-known stochastic processes involves multiplicative noise. In this work we present a systematic mapping of a process with multiplicative noise to a related process with additive noise, which may often be easier to analyse. The mapping is easily understood in the example of the branching process. In a second example we study the random neighbour (or infinite range) contact process which is mapped to an Ornstein-Uhlenbeck process with absorbing wall. The present work might shed some light on absorbing state phase transitions in general, such as the role of conditional expectation values and finite size scaling, and elucidate the meaning of the noise amplitude. While we focus on the physical interpretation of the mapping, we also provide a mathematical derivation. △ Less

Submitted 3 January, 2014; originally announced January 2014.

Comments: 22 pages, 4 figures, IOP style

Journal ref: J. Phys. A: Math. Theor. 47(19) 195001 (2014)

SOC computer simulations

Authors: Gunnar Pruessner

Abstract: The following chapter provides an overview of the techniques used to understand Self-Organised Criticality (SOC) by performing computer simulations. Those are of particular significance in SOC, given its very paradigm, the BTW (Bak-Tang-Wiesenfeld) sandpile, was introduced on the basis of a process that is conveniently implemented as a computer program. The chapter is divided into three sections:… ▽ More The following chapter provides an overview of the techniques used to understand Self-Organised Criticality (SOC) by performing computer simulations. Those are of particular significance in SOC, given its very paradigm, the BTW (Bak-Tang-Wiesenfeld) sandpile, was introduced on the basis of a process that is conveniently implemented as a computer program. The chapter is divided into three sections: In the first section a number of key concepts are introduced, followed by four brief presentations of SOC models which are most commonly investigated or which have played an important part in the development of the field as a whole. The second section is concerned with the basics of scaling with particular emphasis of its role in numerical models of SOC, introducing a number of basic tools for data analysis such as binning, moment analysis and error estimation. The third section is devoted to numerical methods and algorithms as applied to SOC models, addressing typical computational questions with the particular application of SOC in mind. The present chapter is rather technical, but hands-on at the same time, providing practical advice and even code snippets (in C) wherever possible. △ Less

Submitted 14 January, 2013; originally announced January 2013.

Comments: 57 pages, 5 figures, chapter 7 of e-book Self-Organized Criticality Systems, edited by M Aschwanden, OpenAcademicPress

Doppler synchronization of pulsating phases by time delay

Authors: Gunnar Pruessner, Seng Cheang, Henrik Jeldtoft Jensen

Abstract: Synchronization by exchange of pulses is a widespread phenomenon, observed in flashing fireflies, applauding audiences and the neuronal network of the brain. Hitherto the focus has been on integrate-and-fire oscillators. Here we consider entirely analytic time evolution. Oscillators exchange narrow but finite pulses. For any non-zero time lag between the oscillators complete synchronization occurs… ▽ More Synchronization by exchange of pulses is a widespread phenomenon, observed in flashing fireflies, applauding audiences and the neuronal network of the brain. Hitherto the focus has been on integrate-and-fire oscillators. Here we consider entirely analytic time evolution. Oscillators exchange narrow but finite pulses. For any non-zero time lag between the oscillators complete synchronization occurs for any number of oscillators arranged in interaction networks whose adjacency matrix fulfils some simple conditions. The time to synchronization decreases with increasing time lag. △ Less

Submitted 12 December, 2012; originally announced December 2012.

Comments: 4 pages and 4 figures

The average avalanche size in the Manna Model and other models of self-organised criticality

Authors: Gunnar Pruessner

Abstract: The average avalanche size can be calculated exactly in a number of models of self-organised criticality (SOC). While the calculation is straight-forward in one dimension, it is more involved in higher dimensions and further complicated by the presence of different boundary conditions and different forms of external driving. Amplitudes of the leading order are determined analytically and evaluated… ▽ More The average avalanche size can be calculated exactly in a number of models of self-organised criticality (SOC). While the calculation is straight-forward in one dimension, it is more involved in higher dimensions and further complicated by the presence of different boundary conditions and different forms of external driving. Amplitudes of the leading order are determined analytically and evaluated to obtain analytical references for numerical work. A subtle link exists between the procedure to calculate the average avalanche size and the field theory of SOC. △ Less

Submitted 12 September, 2012; v1 submitted 9 August, 2012; originally announced August 2012.

Comments: 7 pages SVJour style, one figure, one table of amplitudes; Revision: added references and clarifications

Mounding in Epitaxial Surface Growth

Authors: Edward Sherman, Gunnar Pruessner

Abstract: For the past two decades the Villain-Lai-Das Sarma equation has served as the theoretical framework for conserved surface growth processes, such as molecular-beam epitaxy. However some phenomena, such as mounding, are yet to be fully understood. In the following, we present a systematic analysis of the full, original Villain-Lai-Das Sarma equation showing that mound forming terms, which should hav… ▽ More For the past two decades the Villain-Lai-Das Sarma equation has served as the theoretical framework for conserved surface growth processes, such as molecular-beam epitaxy. However some phenomena, such as mounding, are yet to be fully understood. In the following, we present a systematic analysis of the full, original Villain-Lai-Das Sarma equation showing that mound forming terms, which should have been included initially on symmetry grounds, are generated under renormalisation. A number of widely studied Langevin equations are recovered as limits or trivial fixed points of the full theory. △ Less

Submitted 13 April, 2012; originally announced April 2012.

Comments: 4 pages, 1 figure

Universality and a numerical εいぷしろん-expansion of the Abelian Manna model below upper critical dimension

Authors: Hoai Nguyen Huynh, Gunnar Pruessner

Abstract: The Abelian Manna model of self-organized criticality is studied on various three-dimensional and fractal lattices. The exponents for avalanche size, duration and area distribution of the model are obtained by using a high-accuracy moment analysis. Together with earlier results on lower-dimensional lattices, the present results reinforce the notion of universality below the upper critical dimensio… ▽ More The Abelian Manna model of self-organized criticality is studied on various three-dimensional and fractal lattices. The exponents for avalanche size, duration and area distribution of the model are obtained by using a high-accuracy moment analysis. Together with earlier results on lower-dimensional lattices, the present results reinforce the notion of universality below the upper critical dimension and allow us to determine the the coefficients of an εいぷしろん-expansion. Rescaling the critical exponents by the lattice dimension and incorporating the random walker dimension, a remarkable relation is observed, satisfied by both regular and fractal lattices. △ Less

Submitted 22 May, 2012; v1 submitted 16 January, 2012; originally announced January 2012.

Comments: 6 pages, 3 figures, 6 tables, submitted to PRE

Journal ref: Phys. Rev. E 85, 061133 (2012)

Abelian Manna model on various lattices in one and two dimensions

Authors: Hoai Nguyen Huynh, Gunnar Pruessner, Lock Yue Chew

Abstract: We perform a high-accuracy moment analysis of the avalanche size, duration and area distribution of the Abelian Manna model on eight two-dimensional and four one-dimensional lattices. The results provide strong support to establish universality of exponents and moment ratios across different lattices and a good survey for the strength of corrections to scaling which are notorious in the Manna univ… ▽ More We perform a high-accuracy moment analysis of the avalanche size, duration and area distribution of the Abelian Manna model on eight two-dimensional and four one-dimensional lattices. The results provide strong support to establish universality of exponents and moment ratios across different lattices and a good survey for the strength of corrections to scaling which are notorious in the Manna universality class. The results are compared against previous work done on Manna model, Oslo model and directed percolation. We also confirm hypothesis of various scaling relations. △ Less

Submitted 2 June, 2011; originally announced June 2011.

Comments: 16 pages, 3 figures, 7 tables, submitted to Journal of Statistical Mechanics

Journal ref: J. Stat. Mech. (2011) P09024

The Edwards-Wilkinson equation with drift and Neumann boundary conditions

Authors: Seng Cheang, Gunnar Pruessner

Abstract: The well known scaling of the Edwards-Wilkinson equation is essentially determined by dimensional analysis. Once a drift term is added, more sophisticated reasoning is required, which initially suggests that the drift term dominates over the diffusion. However, the diffusion term is dangerously irrelevant and the resulting scaling in fact non-trivial. In the present article we compare the resultin… ▽ More The well known scaling of the Edwards-Wilkinson equation is essentially determined by dimensional analysis. Once a drift term is added, more sophisticated reasoning is required, which initially suggests that the drift term dominates over the diffusion. However, the diffusion term is dangerously irrelevant and the resulting scaling in fact non-trivial. In the present article we compare the resulting scaling of the Edwards-Wilkinson equation with drift and Neumann boundary conditions to the published case with Dirichlet boundary conditions. △ Less

Submitted 5 August, 2010; v1 submitted 31 July, 2010; originally announced August 2010.

Comments: 20 pages, IOP style article, in version 2 minor mistakes and typos corrected

Journal ref: J. Phys. A: Math. Theor. 44, 065003 (2011)

The Abelian Manna model on two fractal lattices

Authors: Hoai Nguyen Huynh, Lock Yue Chew, Gunnar Pruessner

Abstract: We analyze the avalanche size distribution of the Abelian Manna model on two different fractal lattices with the same dimension d_g=ln(3)/ln(2), with the aim to probe for scaling behavior and to study the systematic dependence of the critical exponents on the dimension and structure of the lattices. We show that the scaling law D(2-tau)=d_w generalizes the corresponding scaling law on regular latt… ▽ More We analyze the avalanche size distribution of the Abelian Manna model on two different fractal lattices with the same dimension d_g=ln(3)/ln(2), with the aim to probe for scaling behavior and to study the systematic dependence of the critical exponents on the dimension and structure of the lattices. We show that the scaling law D(2-tau)=d_w generalizes the corresponding scaling law on regular lattices, in particular hypercubes, where d_w=2. Furthermore, we observe that the lattice dimension d_g, the fractal dimension of the random walk on the lattice d_w, and the critical exponent D, form a plane in 3D parameter space, i.e. they obey the linear relationship D=0.632(3) d_g + 0.98(1) d_w - 0.49(3). △ Less

Submitted 30 June, 2010; originally announced June 2010.

Comments: 4 pages, 3 figures, 3 tables, submitted to PRE as a Brief Report

Journal ref: Phys. Rev. E 82, 042103 (2010)

Comment on: "Avalanches and Non-Gaussian Fluctuations of the Global Velocity of Imbibition Fronts"

Authors: Gunnar Pruessner

Abstract: Comment on R. Planet, S. Santucci, J. Ortin, Phys. Rev. Lett. 102, 094502 (2009) about the rescaling of the data and the data collapse. Reply to be found here. Comment on R. Planet, S. Santucci, J. Ortin, Phys. Rev. Lett. 102, 094502 (2009) about the rescaling of the data and the data collapse. Reply to be found here. △ Less

Submitted 17 June, 2010; originally announced June 2010.

Tuning- and order parameter in the SOC ensemble

Authors: Ole Peters, Gunnar Pruessner

Abstract: The one-dimensional Oslo model is studied under self-organized criticality (SOC) conditions and under absorbing state (AS) conditions. While the activity signals the phase transition under AS conditions by a sudden increase, this is not the case under SOC conditions. The scaling parameters of the activity are found to be identical under SOC and AS conditions, but in SOC the activity lacks a pick… ▽ More The one-dimensional Oslo model is studied under self-organized criticality (SOC) conditions and under absorbing state (AS) conditions. While the activity signals the phase transition under AS conditions by a sudden increase, this is not the case under SOC conditions. The scaling parameters of the activity are found to be identical under SOC and AS conditions, but in SOC the activity lacks a pickup. △ Less

Submitted 11 December, 2009; originally announced December 2009.

Comments: 4 pages, 2 figures

Reply to "Comment on `Self-organized Criticality and Absorbing States: Lessons from the Ising Model'"

Authors: Gunnar Pruessner, Ole Peters

Abstract: In [Braz. J. Phys. 30, 27 (2000)] Dickman et al. suggested that self-organized criticality can be produced by coupling the activity of an absorbing state model to a dissipation mechanism and adding an external drive. We analyzed the proposed mechanism in [Phys. Rev. E 73, 025106R (2006)] and found that if this mechanism is at work, the finite-size scaling found in self-organized criticality will… ▽ More In [Braz. J. Phys. 30, 27 (2000)] Dickman et al. suggested that self-organized criticality can be produced by coupling the activity of an absorbing state model to a dissipation mechanism and adding an external drive. We analyzed the proposed mechanism in [Phys. Rev. E 73, 025106R (2006)] and found that if this mechanism is at work, the finite-size scaling found in self-organized criticality will depend on the details of the implementation of dissipation and driving. In the preceding comment [Phys. Rev. E XX, XXXX (2008)] Alava et al. show that one avalanche exponent in the AS approach becomes independent of dissipation and driving. In our reply we clarify their findings and put them in the context of the original article. △ Less

Submitted 4 April, 2008; originally announced April 2008.

Comments: 4 pages, REVTeX (draft)

On the scaling of probability density functions with apparent power-law exponents less than unity

Authors: Kim Christensen, Nadia Farid, Gunnar Pruessner, Matthew Stapleton

Abstract: We derive general properties of the finite-size scaling of probability density functions and show that when the apparent exponent \tautilde of a probability density is less than 1, the associated finite-size scaling ansatz has a scaling exponent τたうequal to 1, provided that the fraction of events in the universal scaling part of the probability density function is non-vanishing in the thermodynami… ▽ More We derive general properties of the finite-size scaling of probability density functions and show that when the apparent exponent \tautilde of a probability density is less than 1, the associated finite-size scaling ansatz has a scaling exponent τたうequal to 1, provided that the fraction of events in the universal scaling part of the probability density function is non-vanishing in the thermodynamic limit. We find the general result that τたう>=1 and τたう>=\tautilde. Moreover, we show that if the scaling function G(x) approaches a non-zero constant for small arguments, \lim_{x-> 0} G(x) > 0, then τたう=\tautilde. However, if the scaling function vanishes for small arguments, \lim_{x-> 0} G(x) = 0, then τたう=1, again assuming a non-vanishing fraction of universal events. Finally, we apply the formalism developed to examples from the literature, including some where misunderstandings of the theory of scaling have led to erroneous conclusions. △ Less

Submitted 30 January, 2008; originally announced January 2008.

Comments: 6 pages, 2 figures, EPJB style

Equivalence of conditional and external field ensembles in absorbing state phase transitions

Authors: Gunnar Pruessner

Abstract: I comment on the relation between two sampling methods for absorbing state models. It is shown that a certain ensemble without external field conditional to activity coincides with the unconditional ensemble for sufficiently small external field. The actual physical processes involved are identical and the derivation of the identity of the scaling behaviour relies on a single (established) scali… ▽ More I comment on the relation between two sampling methods for absorbing state models. It is shown that a certain ensemble without external field conditional to activity coincides with the unconditional ensemble for sufficiently small external field. The actual physical processes involved are identical and the derivation of the identity of the scaling behaviour relies on a single (established) scaling law. While the conditional ensemble by construction does not contain information about the system with large external field, it contains all information about the limit of vanishing external field and about the vicinity of the critical point: Finite size scaling as well as critical scaling in the temperature-like variable or in (small) external field. △ Less

Submitted 6 December, 2007; originally announced December 2007.

Comments: 6 pages (single column), 2 figures, RevTeX

Journal ref: Phys. Rev. E 76, 061103 (2007)