BZPEERPapers, people, evidence and research pathsDiscover
Research tools
Network
Opportunities
Account
Source author record
Gunnar Pruessner appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.
Catalog footprint
Research graph
Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.
BZPEER is loading the nearby papers, people, topics and institutions for this page.
Published work
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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).
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.