Source author record

Erwin Frey

Erwin Frey appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

68works
17topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

68 published item(s)

preprint2022arXiv

Bridging scales in a multiscale pattern-forming system

Self-organized pattern formation is vital for many biological processes. Reaction-diffusion models have advanced our understanding of how biological systems develop spatial structures, starting from homogeneity. However, biological processes inherently involve multiple spatial and temporal scales and transition from one pattern to another over time, rather than progressing from homogeneity to a pattern. To deal with such multiscale systems, coarse-graining methods are needed that allow the dynamics to be reduced to the relevant degrees of freedom at large scales, but without losing information about the patterns at the small scales. Here, we present a semi-phenomenological approach which exploits mass-conservation in pattern formation, and enables to reconstruct information about patterns from the large-scale dynamics. The basic idea is to partition the domain into distinct regions (coarse-grain) and determine instantaneous dispersion relations in each region, which ultimately inform about local pattern-forming instabilities. We illustrate our approach by studying the Min system, a paradigmatic model for protein pattern formation. By performing simulations, we first show that the Min system produces multiscale patterns in a spatially heterogeneous geometry. This prediction is confirmed experimentally by in vitro reconstitution of the Min system. Using a recently developed theoretical framework for mass-conserving reaction-diffusion systems, we show that the spatiotemporal evolution of the total protein densities on large-scales reliably predicts the pattern-forming dynamics. Our approach provides an alternative and versatile theoretical framework for complex systems where analytical coarse-graining methods are not applicable, and can in principle be applied to a wide range of systems with an underlying conservation law.

preprint2021arXiv

Subdiffusive Activity Spreading in the Diffusive Epidemic Process

The diffusive epidemic process is a paradigmatic example of an absorbing state phase transition in which healthy and infected individuals spread with different diffusion constants. Using stochastic activity spreading simulations in combination with finite-size scaling analyses we reveal two qualitatively different processes that characterize the critical dynamics: subdiffusive propagation of infection clusters and diffusive fluctuations in the healthy population. This suggests the presence of a strong-coupling regime and sheds new light on a longstanding debate about the theoretical classification of the system.

preprint2020arXiv

Cytosolic flow induced symmetry breaking in a conceptual polarity model

Important cellular processes, such as cell motility and cell division, are coordinated by cell polarity, which is determined by the non-uniform distribution of certain proteins. Such protein patterns form via an interplay of protein reactions and protein transport. Since Turing's seminal work, the formation of protein patterns resulting from the interplay between reactions and diffusive transport has been widely studied. Over the last few years, increasing evidence shows that also advective transport, resulting from cytosolic and cortical flows, is present in many cells. However, it remains unclear how and whether these flows contribute to protein-pattern formation. To address this question, we use a minimal model that conserves the total protein mass to characterize the effects of cytosolic flow on pattern formation. Combining a linear stability analysis with numerical simulations, we find that membrane-bound protein patterns propagate against the direction of cytoplasmic flow with a speed that is maximal for intermediate flow speed. We show that the mechanism underlying this pattern propagation relies on a higher protein influx on the upstream side of the pattern compared to the downstream side. Furthermore, we find that cytosolic flow can change the membrane pattern qualitatively from a peak pattern to a mesa pattern. Finally, our study shows that a non-uniform flow profile can induce pattern formation by triggering a regional lateral instability.

preprint2020arXiv

Stochastic Yield Catastrophes and Robustness in Self-Assembly

A guiding principle in self-assembly is that, for high production yield, nucleation of structures must be significantly slower than their growth. However, details of the mechanism that impedes nucleation are broadly considered irrelevant. Here, we analyze self-assembly into finite-sized target structures employing mathematical modeling. We investigate two key scenarios to delay nucleation: (i) by introducing a slow activation step for the assembling constituents and, (ii) by decreasing the dimerization rate. These scenarios have widely different characteristics. While the dimerization scenario exhibits robust behavior, the activation scenario is highly sensitive to demographic fluctuations. These demographic fluctuations ultimately disfavor growth compared to nucleation and can suppress yield completely. The occurrence of this stochastic yield catastrophe does not depend on model details but is generic as soon as number fluctuations between constituents are taken into account. On a broader perspective, our results reveal that stochasticity is an important limiting factor for self-assembly and that the specific implementation of the nucleation process plays a significant role in determining the yield.

preprint2020arXiv

Topological phase transition in coupled rock-paper-scissor cycles

A hallmark of topological phases is the occurrence of topologically protected modes at the system`s boundary. Here we find topological phases in the antisymmetric Lotka-Volterra equation (ALVE). The ALVE is a nonlinear dynamical system and describes, e.g., the evolutionary dynamics of a rock-paper-scissors cycle. On a one-dimensional chain of rock-paper-scissor cycles, topological phases become manifest as robust polarization states. At the transition point between left and right polarization, solitonic waves are observed. This topological phase transition lies in symmetry class $D$ within the "ten-fold way" classification as also realized by 1D topological superconductors.

preprint2019arXiv

Pattern localization to a domain edge

The formation of protein patterns inside cells is generically described by reaction-diffusion models. The study of such systems goes back to Turing, who showed how patterns can emerge from a homogenous steady state when two reactive components have different diffusivities (e.g. membrane-bound and cytosolic states). However, in nature, systems typically develop in a heterogeneous environment, where upstream protein patterns affect the formation of protein patterns downstream. Examples for this are the polarization of Cdc42 adjacent to the previous bud-site in budding yeast, and the formation of an actin-recruiter ring that forms around a PIP3 domain in macropinocytosis. This suggests that previously established protein patterns can serve as a template for downstream proteins and that these downstream proteins can 'sense' the edge of the template. A mechanism for how this edge sensing may work remains elusive. Here we demonstrate and analyze a generic and robust edge-sensing mechanism, based on a two-component mass-conserving reaction-diffusion (McRD) model. Our analysis is rooted in a recently developed theoretical framework for McRD systems, termed local equilibria theory. We extend this framework to capture the spatially heterogeneous reaction kinetics due to the template. This enables us to graphically construct the stationary patterns in the phase space of the reaction kinetics. Furthermore, we show that the protein template can trigger a regional mass-redistribution instability near the template edge, leading to the accumulation of protein mass, which eventually results in a stationary peak at the template edge. We show that simple geometric criteria on the reactive nullcline's shape predict when this edge-sensing mechanism is operational. Thus, our results provide guidance for future studies of biological systems, and for the design of synthetic pattern forming systems.

preprint2016arXiv

A nonequilibrium diffusion and capture mechanism ensures tip-localization of regulating proteins on dynamic filaments

Diffusive motion of regulatory enzymes on biopolymers with eventual capture at a reaction site is a common feature in cell biology. Using a lattice gas model we study the impact of diffusion and capture for a microtubule polymerase and a depolymerase. Our results show that the capture mechanism localizes the proteins and creates large-scale spatial correlations. We develop an analytic approximation that globally accounts for relevant correlations and yields results that are in excellent agreement with experimental data. Our results show that diffusion and capture operates most efficiently at cellular enzyme concentrations which points to in vivo relevance.

preprint2016arXiv

Active Curved Polymers form Vortex Patterns on Membranes

Recent in vitro experiments with FtsZ polymers show self-organization into different dynamic patterns, including structures reminiscent of the bacterial Z-ring. We model FtsZ polymers as active particles moving along chiral, circular paths by Brownian dynamics simulations and a Boltzmann approach. Our two conceptually different methods point to a generic phase behavior. At intermediate particle densities, we find self-organization into vortex structures including closed rings. Moreover, we show that the dynamics at the onset of pattern formation is described by a generalized complex Ginzburg-Landau equation.

preprint2016arXiv

Binary Mixtures of Particles with Different Diffusivities Demix

The influence of size differences, shape, mass and persistent motion on phase separation in binary mixtures has been intensively studied. Here we focus on the exclusive role of diffusivity differences in binary mixtures of equal-sized particles. We find an effective attraction between the less diffusive particles, which are essentially caged in the surrounding species with the higher diffusion constant. This effect leads to phase separation for systems above a critical size: A single close-packed cluster made up of the less diffusive species emerges. Experiments for testing of our predictions are outlined.

preprint2015arXiv

Emergence and Persistence of Collective Cell Migration on Small Circular Micropatterns

The spontaneous formation of vortices is a hallmark of collective cellular activity. Here, we study the onset and persistence of coherent angular motion (CAMo) as a function of the number of cells $N$ confined in circular micropatterns. We find that the persistence of CAMo increases with $N$ but exhibits a pronounced discontinuity accompanied by a geometric rearrangement of cells to a configuration containing a central cell. Computer simulations based on a generalized Potts model reproduce the emergence of vortex states and show in agreement with experiment that their stability depends on the interplay of spatial arrangement and internal polarization of neighboring cells. Hence, the distinct migrational states in finite size ensembles reveal significant insight into the local interaction rules guiding collective migration.

preprint2015arXiv

Evolutionary games of condensates in coupled birth-death processes

Condensation phenomena arise through a collective behaviour of particles. They are observed in both classical and quantum systems, ranging from the formation of traffic jams in mass transport models to the macroscopic occupation of the energetic ground state in ultra-cold bosonic gases (Bose-Einstein condensation). Recently, it has been shown that a driven and dissipative system of bosons may form multiple condensates. Which states become the condensates has, however, remained elusive thus far. The dynamics of this condensation are described by coupled birth-death processes, which also occur in evolutionary game theory. Here, we apply concepts from evolutionary game theory to explain the formation of multiple condensates in such driven-dissipative bosonic systems. We show that vanishing of relative entropy production determines their selection. The condensation proceeds exponentially fast, but the system never comes to rest. Instead, the occupation numbers of condensates may oscillate, as we demonstrate for a rock-paper-scissors game of condensates.

preprint2015arXiv

How turbulence regulates biodiversity in systems with cyclic competition

Cyclic, nonhierarchical interactions among biological species represent a general mechanism by which ecosystems are able to maintain high levels of biodiversity. However, species coexistence is often possible only in spatially extended systems with a limited range of dispersal, whereas in well-mixed environments models for cyclic competition often lead to a loss of biodiversity. Here we consider the dispersal of biological species in a fluid environment, where mixing is achieved by a combination of advection and diffusion. In particular, we perform a detailed numerical analysis of a model composed of turbulent advection, diffusive transport, and cyclic interactions among biological species in two spatial dimensions and discuss the circumstances under which biodiversity is maintained when external environmental conditions, such as resource supply, are uniform in space. Cyclic interactions are represented by a model with three competitors, resembling the children's game of rock-paper-scissors, whereas the flow field is obtained from a direct numerical simulation of two-dimensional turbulence with hyperviscosity. It is shown that the space-averaged dynamics undergoes bifurcations as the relative strengths of advection and diffusion compared to biological interactions are varied.

preprint2015arXiv

Quantifying protein diffusion and capture on filaments

The functional relevance of regulating proteins is often limited to specific binding sites such as the ends of microtubules or actin-filaments. A localization of proteins on these functional sites is of great importance. We present a quantitative theory for a diffusion and capture process, where proteins diffuse on a filament and stop diffusing when reaching the filament's end. It is found that end-association after one-dimensional diffusion is the main source for tip-localization of such proteins. As a consequence, diffusion and capture is highly efficient in enhancing the reaction velocity of enzymatic reactions, where proteins and filament ends are to each other as enzyme and substrate. We show that the reaction velocity can effectively be described within a Michaelis-Menten framework. Together one-dimensional diffusion and capture beats the (three-dimensional) Smoluchowski diffusion limit for the rate of protein association to filament ends.

preprint2015arXiv

The Emergence of Cooperation from a Single Mutant during Microbial Life-Cycles

Cooperative behavior is widespread in nature, even though cooperating individuals always run the risk to be exploited by free-riders. Population structure effectively promotes cooperation given that a threshold in the level of cooperation was already reached. However, the question how cooperation can emerge from a single mutant, which cannot rely on a benefit provided by other cooperators, is still puzzling. Here, we investigate this question for a well-defined but generic situation based on typical life-cycles of microbial populations where individuals regularly form new colonies followed by growth phases. We analyze two evolutionary mechanisms favoring cooperative behavior and study their strength depending on the inoculation size and the length of a life-cycle. In particular, we find that population bottlenecks followed by exponential growth phases strongly increase the survival and fixation probabilities of a single cooperator in a free-riding population.

preprint2015arXiv

Traffic Dynamic Instability

Here we study a driven lattice gas model for microtubule depolymerizing molecular motors, where traffic jams of motors induce stochastic switching between microtubule growth and shrinkage. We term this phenomenon \enquote{traffic dynamic instability} because it is reminiscent of microtubule dynamic instability [T. Mitchison and M. Kirschner, Nature 312, 237 (1984)]. The intermittent dynamics of growth and shrinking emerges from the interplay between the arrival of motors at the microtubule tip, motor induced depolymerization, and motor detachment from the tip. The switching dynamics correlates with low and high motor density on the lattice. This leads to an effectively bistable particle density in the system. A refined domain wall theory predicts this transient appearance of different phases in the system. The theoretical results are supported by stochastic simulations.

preprint2014arXiv

Chemical warfare and survival strategies in bacterial range expansions

Dispersal of species is a fundamental ecological process in the evolution and maintenance of biodiversity. Limited control over ecological parameters has hindered progress in understanding of what enables species to colonise new area, as well as the importance of inter-species interactions. Such control is necessary to construct reliable mathematical models of ecosystems. In our work, we studied dispersal in the context of bacterial range expansions and identified the major determinants of species coexistence for a bacterial model system of three Escherichia coli strains (toxin producing, sensitive, and resistant). Genetic engineering allowed us to tune strain growth rates and to design different ecological scenarios (cyclic and hierarchical). We found that coexistence of all strains depended on three strongly interdependent factors: composition of inoculum, relative strain growth rates, and effective toxin range. Robust agreement between our experiments and a thoroughly calibrated computational model enabled us to extrapolate these intricate interdependencies in terms of phenomenological biodiversity laws. Our mathematical analysis also suggested that cyclic dominance between strains is not a prerequisite for coexistence in competitive range expansions. Instead, robust three-strain coexistence required a balance between growth rates and either a reduced initial ratio of the toxin-producing strain, or a sufficiently short toxin range.

preprint2014arXiv

Defect-Mediated Phase Transitions in Active Soft Matter

How do topological defects affect the degree of order in active matter? To answer this question we investigate an agent-based model of self-propelled particles, which accounts for polar alignment and short-ranged repulsive interactions. For strong alignment forces we find collectively moving polycrystalline states with fluctuating networks of grain boundaries. In the regime where repulsive forces dominate, the fluctuations generated by the active system give rise to quasi-long-range transitional order, but---unlike thermal system---without creating topological defects.

preprint2014arXiv

Dynamics of a Semiflexible Polymer or Polymer Ring in Shear Flow

Polymers exposed to shear flow exhibit a rich tumbling dynamics. While rigid rods rotate on Jeffery orbits, flexible polymers stretch and coil up during tumbling. Theoretical results show that in both of these asymptotic regimes the tumbling frequency f_c in a linear shear flow of strength γscales as a power law Wi^(2/3) in the Weissenberg number Wi=γτ, where τis a characteristic time of the polymer's relaxational dynamics. For flexible polymers these theoretical results are well confirmed by experimental single molecule studies. However, for the intermediate semiflexible regime the situation is less clear. Here we perform extensive Brownian dynamics simulations to explore the tumbling dynamics of semiflexible polymers over a broad range of shear strength and the polymer's persistence length l_p. We find that the Weissenberg number alone does not suffice to fully characterize the tumbling dynamics, and the classical scaling law breaks down. Instead, both the polymer's stiffness and the shear rate are relevant control parameters. Based on our Brownian dynamics simulations we postulate that in the parameter range most relevant for cytoskeletal filaments there is a distinct scaling behavior with f_c τ*=Wi^(3/4) f_c (x) with Wi=γτ* and the scaling variable x=(l_p/L)(Wi)^(-1/3); here τ* is the time the polymer's center of mass requires to diffuse its own contour length L. Comparing these results with experimental data on F-actin we find that the Wi^(3/4) scaling law agrees quantitatively significantly better with the data than the classical Wi^(2/3) law. Finally, we extend our results to single ring polymers in shear flow, and find similar results as for linear polymers with slightly different power laws.

preprint2014arXiv

Effective two-dimensional model does not account for geometry sensing by self-organized proteins patterns - Supplementary document

Here we provide a thorough discussion of the model for Min protein dynamics proposed by Schweizer et al. [11]. The manuscript serves as supplementary document for our letter to the editor to appear in PNAS. Our analysis is based on the original COMSOL simulation files that were used for the publication. We show that all computational data in Schweizer et al. rely on exploitation of simulation artifacts and various unmentioned modifications of model parameters that strikingly contradict the experimental setup and experimental data. We find that the model neither accounts for MinE membrane interactions nor for any observed MinDE protein patterns. All conclusions drawn from the computational model are void. There is no evidence at all that persistent MinE membrane binding has any role in geometry sensing.

preprint2014arXiv

Mobility-Dependent Selection of Competing Strategy Associations

Standard models of population dynamics focus on the the interaction, survival, and extinction of the competing species individually. Real ecological systems, however, are characterized by an abundance of species (or strategies, in the terminology of evolutionary-game theory) that form intricate, complex interaction networks. The description of the ensuing dynamics may be aided by studying associations of certain strategies rather than individual ones. Here we show how such a higher-level description can bear fruitful insight. Motivated from different strains of colicinogenic Escherichia coli bacteria, we investigate a four-strategy system which contains a three-strategy cycle and a neutral alliance of two strategies. We find that the stochastic, spatial model exhibits a mobility-dependent selection of either the three-strategy cycle or of the neutral pair. We analyze this intriguing phenomenon numerically and analytically.

preprint2014arXiv

Molecular Mechanisms for Microtubule Length Regulation by Kinesin-8 and XMAP215 Proteins

The cytoskeleton is regulated by a plethora of enzymes that influence the stability and dynamics of cytoskeletal filaments. Molecular motors of the kinesin-8 protein family depolymerise microtubules in a length-dependent manner, and experimental and theoretical evidence suggest a role for kinesin-8 in the dynamic regulation of microtubules. However, so far the detailed molecular mechanisms how these molecular motors interact with the growing microtubule tip remain elusive. Here we investigate two interaction scenarios for kinesin-8 and the microtubule tip. We give a comprehensive analysis of regimes where length-regulation is possible and characterise how the stationary length depends on the biochemical rates and the bulk concentrations of the various proteins. For a neutral scenario, where microtubules grow irrespective of whether the microtubule tip is occupied by a molecular motor, length regulation is possible only for a narrow range of biochemical rates and limited to small polymerisation rates. In contrast, for an inhibition scenario, where the presence of a motor at the microtubule tip inhibits microtubule growth, the regime of length regulation is extremely broad and includes high growth rates. These results also apply to situations where polymerising enzymes, like XMAP215, and kinesin-8 mutually exclude each other from the microtubule tip. We also investigate the stochastic dynamics of the two scenarios. While for the neutral scenario length is tightly controlled, length dynamics is intermittent for the inhibition scenario and exhibits extended periods of microtubule growth and shrinkage, reminiscent of microtubule dynamic instability. On a broader perspective, the set of models established in this work quite generally suggests that mutual exclusion of molecules at the ends of cytoskeletal filaments is an important factor for filament dynamics and regulation.

preprint2014arXiv

Numerical treatment of the Boltzmann equation for self-propelled particle systems

Kinetic theories constitute one of the most promising tools to decipher the characteristic spatio-temporal dynamics in systems of actively propelled particles. In this context, the Boltzmann equation plays a pivotal role, since it provides a natural translation between a particle-level description of the system's dynamics and the corresponding hydrodynamic fields. Yet, the intricate mathematical structure of the Boltzmann equation substantially limits the progress toward a full understanding of this equation by solely analytical means. Here, we propose a general framework to numerically solve the Boltzmann equation for self-propelled particle systems in two spatial dimensions and with arbitrary boundary conditions. We discuss potential applications of this numerical framework to active matter systems, and use the algorithm to give a detailed analysis to a model system of self-propelled particles with polar interactions. In accordance with previous studies, we find that spatially homogeneous isotropic and broken symmetry states populate two distinct regions in parameter space, which are separated by a narrow region of spatially inhomogeneous, density-segregated moving patterns. We find clear evidence that these three regions in parameter space are connected by first order phase transitions, and that the transition between the spatially homogeneous isotropic and polar ordered phases bears striking similarities to liquid-gas phase transitions in equilibrium systems. Within the density segregated parameter regime, we find a novel stable limit-cycle solution of the Boltzmann equation, which consists of parallel lanes of polar clusters moving in opposite directions, so as to render the overall symmetry of the system's ordered state nematic, despite purely polar interactions on the level of single particles.

preprint2014arXiv

Periodic vs. intermittent adaptive cycles in quasispecies co-evolution

We study an abstract model for the co-evolution between mutating viruses and the adaptive immune system. In sequence space, these two populations are localized around transiently dominant strains. Delocalization or error thresholds exhibit a novel interdependence because immune response is conditional on the viral attack. An evolutionary chase is induced by stochastic fluctuations and can occur via periodic or intermittent cycles. Using simulations and stochastic analysis, we show how the transition between these two dynamic regimes depends on mutation rate, immune response, and population size.

preprint2014arXiv

Specialization and Bet Hedging in Heterogeneous Populations

Phenotypic heterogeneity is a strategy commonly used by bacteria to rapidly adapt to changing environmental conditions. Here, we study the interplay between phenotypic heterogeneity and genetic diversity in spatially extended populations. By analyzing the spatio-temporal dynamics, we show that the level of mobility and the type of competition qualitatively influence the persistence of phenotypic heterogeneity. While direct competition generally promotes persistence of phenotypic heterogeneity, specialization dominates in models with indirect competition irrespective of the degree of mobility.

preprint2013arXiv

A Critical Assessment of the Boltzmann Approach for Active Systems

Generic models of propelled particle systems posit that the emergence of polar order is driven by the competition between local alignment and noise. Although this notion has been confirmed employing the Boltzmann equation, the range of applicability of this equation remains elusive. We introduce a broad class of mesoscopic collision rules and analyze the prerequisites for the emergence of polar order in the framework of kinetic theory. Our findings suggest that a Boltzmann approach is appropriate for weakly aligning systems but is incompatible with experiments on cluster forming systems.

preprint2013arXiv

Coexistence and Survival in Conservative Lotka-Volterra Networks

Analyzing coexistence and survival scenarios of Lotka-Volterra (LV) networks in which the total biomass is conserved is of vital importance for the characterization of long-term dynamics of ecological communities. Here, we introduce a classification scheme for coexistence scenarios in these conservative LV models and quantify the extinction process by employing the Pfaffian of the network's interaction matrix. We illustrate our findings on global stability properties for general systems of four and five species and find a generalized scaling law for the extinction time.

preprint2013arXiv

Cooperative effects enhance the transport properties of molecular spider teams

Molecular spiders are synthetic molecular motors based on DNA nanotechnology. While natural molecular motors have evolved towards very high efficiency, it remains a major challenge to develop efficient designs for man-made molecular motors. Inspired by biological motor proteins like kinesin and myosin, molecular spiders comprise a body and several legs. The legs walk on a lattice that is coated with substrate which can be cleaved catalytically. We propose a novel molecular spider design in which n spiders form a team. Our theoretical considerations show that coupling several spiders together alters the dynamics of the resulting team significantly. Although spiders operate at a scale where diffusion is dominant, spider teams can be tuned to behave nearly ballistic, which results in fast and predictable motion. Based on the separation of time scales of substrate and product dwell times, we develop a theory which utilises equivalence classes to coarse-grain the micro-state space. In addition, we calculate diffusion coefficients of the spider teams, employing a mapping of an n-spider team to an n-dimensional random walker on a confined lattice. We validate these results with Monte Carlo simulations and predict optimal parameters of the molecular spider team architecture which makes their motion most directed and maximally predictable.

preprint2013arXiv

Fluctuation Effects in the Pair-Annihilation Process with Lévy Dynamics

We investigate the density decay in the pair-annihilation process A+A->0 in the case when the particles perform anomalous diffusion on a cubic lattice. The anomalous diffusion is realized via Lévy flights, which are characterized by long-range jumps and lead to superdiffusive behavior. As a consequence, the critical dimension depends continuously on the control parameter of the Lévy flight distribution. This instance is used to study the system close to the critical dimension by means of the nonperturbative renormalization group theory. Close to the critical dimension, the assumption of well-stirred reactants is violated by anticorrelations between the particles, and the law of mass action breaks down. The breakdown of the law of mass action is known to be caused by long-range fluctuations. We identify three interrelated consequences of these fluctuations. First, despite being a nonuniversal quantity and thus depending on the microscopic details, the renormalized reaction rate can be approximated by a universal law close to the critical dimension. The emergence of universality relies on the fact that long-range fluctuations suppress the influence of the underlying microscopic details. Second, as criticality is approached, the macroscopic reaction rate decreases such that the law of mass action loses its significance. And third, additional nonanalytic power law corrections complement the analytic law of mass action term. An increasing number of those corrections accumulate and give an essential contribution as the critical dimension is approached. We test our findings for two implementations of Lévy flights that differ in the way they cross over to the normal diffusion.

preprint2013arXiv

Global attractors and extinction dynamics of cyclically competing species

Transitions to absorbing states are of fundamental importance in non-equilibrium physics as well as ecology. In ecology, absorbing states correspond to the extinction of species. We here study the spatial population dynamics of three cyclically interacting species. The interaction scheme comprises both direct competition between species as in the cyclic Lotka-Volterra model, and separated selection and reproduction processes as in the May-Leonard model. We show that the dynamic processes leading to the transient maintenance of biodiversity are closely linked to attractors of the nonlinear dynamics for the overall species' concentrations. The characteristics of these global attractors change qualitatively at certain threshold values of the mobility, and depend on the relative strength of the different types of competition between species. They give information about the scaling of extinction times with the system size and thereby the stability of biodiversity. We define an effective free energy as the negative logarithm of the probability to find the system in a specific global state before reaching one of the absorbing states. The global attractors then correspond to minima of this effective energy landscape and determine the most probable values for the species' global concentrations. As in equilibrium thermodynamics, qualitative changes in the effective free energy landscape indicate and characterize the underlying non-equilibrium phase transitions. We provide the complete phase diagrams for the population dynamics, and give a comprehensive analysis of the spatio-temporal dynamics and routes to extinction in the respective phases.

preprint2013arXiv

Language change in a multiple group society

The processes leading to change in languages are manifold. In order to reduce ambiguity in the transmission of information, agreement on a set of conventions for recurring problems is favored. In addition to that, speakers tend to use particular linguistic variants associated with the social groups they identify with. The influence of other groups propagating across the speech community as new variant forms sustains the competition between linguistic variants. With the utterance selection model, an evolutionary description of language change, Baxter et al. [Phys. Rev. E 73, 046118 (2006)] have provided a mathematical formulation of the interactions inside a group of speakers, exploring the mechanisms that lead to or inhibit the fixation of linguistic variants. In this paper, we take the utterance selection model one step further by describing a speech community consisting of multiple interacting groups. Tuning the interaction strength between groups allows us to gain deeper understanding about the way in which linguistic variants propagate and how their distribution depends on the group partitioning. Both for the group size and the number of groups we find scaling behaviors with two asymptotic regimes. If groups are strongly connected, the dynamics is that of the standard utterance selection model, whereas if their coupling is weak, the magnitude of the latter along with the system size governs the way consensus is reached. Furthermore, we find that a high influence of the interlocutor on a speaker's utterances can act as a counterweight to group segregation.

preprint2013arXiv

Long-Range and Many-Body Effects in Coagulation Processes

We study the problem of diffusing particles which coalesce upon contact. With the aid of a non-perturbative renormalization group, we first analyze the dynamics emerging below the critical dimension two, where strong fluctuations imply anomalously slow decay. Above two dimensions, the long-time, low-density behavior is known to conform with the law of mass action. For this case, we establish an exact mapping between the physics at the microscopic scale (lattice structure, particle shape and size) and the macroscopic decay rate in the law of mass action. In addition, we identify a term violating this classical law. It originates in long-range and many-particle fluctuations and is a simple, universal function of the macroscopic decay rate.

preprint2013arXiv

Mobility, fitness collection, and the breakdown of cooperation

The spatial arrangement of individuals is thought to overcome the dilemma of cooperation: When cooperators engage in clusters they might share the benefit of cooperation while being more protected against non-cooperating individuals, which benefit from cooperation but save the cost of cooperation. This is paradigmatically shown by the spatial prisoner's dilemma model. Here, we study this model in one and two spatial dimensions, but explicitly take into account that in biological setups fitness collection and selection are separated processes occurring mostly on vastly different time scales. This separation is particularly important to understand the impact of mobility on the evolution of cooperation. We find that even small diffusive mobility strongly restricts cooperation since it enables non-cooperative individuals to invade cooperative clusters. Thus, in most biological scenarios, where the mobility of competing individuals is an irrefutable fact, the spatial prisoner's dilemma alone cannot explain stable cooperation but additional mechanisms are necessary for spatial structure to promote the evolution of cooperation. The breakdown of cooperation is analyzed in detail. We confirm the existence of a phase transition, here controlled by mobility and costs, which distinguishes between purely cooperative and non-cooperative absorbing states. While in one dimension the model is in the class of the Voter Model, it belongs to the Directed Percolation (DP) universality class in two dimensions.

preprint2013arXiv

Role of particle conservation in self-propelled particle systems

Actively propelled particles undergoing dissipative collisions are known to develop a state of spatially distributed coherently moving clusters. For densities larger than a characteristic value clusters grow in time and form a stationary well-ordered state of coherent macroscopic motion. In this work we address two questions: (i) What is the role of the particles' aspect ratio in the context of cluster formation, and does the particle shape affect the system's behavior on hydrodynamic scales? (ii) To what extent does particle conservation influence pattern formation? To answer these questions we suggest a simple kinetic model permitting to depict some of the interaction properties between freely moving particles and particles integrated in clusters. To this end, we introduce two particle species: single and cluster particles. Specifically, we account for coalescence of clusters from single particles, assembly of single particles on existing clusters, collisions between clusters, and cluster disassembly. Coarse-graining our kinetic model, (i) we demonstrate that particle shape (i.e. aspect ratio) shifts the scale of the transition density, but does not impact the instabilities at the ordering threshold. (ii) We show that the validity of particle conservation determines the existence of a longitudinal instability, which tends to amplify density heterogeneities locally, and in turn triggers a wave pattern with wave vectors parallel to the axis of macroscopic order. If the system is in contact with a particle reservoir this instability vanishes due to a compensation of density heterogeneities.

preprint2013arXiv

Stability of localized wave fronts in bistable systems

Localized wave fronts are a fundamental feature of biological systems from cell biology to ecology. Here, we study a broad class of bistable models subject to self-activation, degradation and spatially inhomogeneous activating agents. We determine the conditions under which wave-front localization is possible and analyze the stability thereof with respect to extrinsic perturbations and internal noise. It is found that stability is enhanced upon regulating a positional signal and, surprisingly, also for a low degree of binding cooperativity. We further show a contrasting impact of self-activation to the stability of these two sources of destabilization.

preprint2013arXiv

Understanding Collective Dynamics of Soft Active Colloids by Binary Scattering

Collective motion in actively propelled particle systems is triggered on the very local scale by nucleation of coherently moving units consisting of just a handful of particles. These units grow and merge over time, ending up in a long-range ordered, coherently-moving state. So far, there exists no bottom-up understanding of how the microscopic dynamics and interactions between the constituents are related to the system's ordering instability. In this paper, we study a class of models for propelled colloids allowing an explicit treatment of the microscopic details of the collision process. Specifically, the model equations are Newtonian equations of motion with separate force terms for particles' driving, dissipation and interaction forces. Focusing on dilute particle systems, we analyze the binary scattering behavior for these models, and determine-based on the microscopic dynamics-the corresponding collision-rule, i.e., the mapping of pre-collisional velocities and impact parameter on post-collisional velocities. By studying binary scattering we also find that the considered models for active colloids share the same principle for parallel alignment: the first incoming particle (with respect to the center of collision) is aligned to the second particle as a result of the encounter. This behavior is distinctively different to alignment in non-driven dissipative gases. Moreover, the obtained collision rule lends itself as a starting point to apply kinetic theory for propelled particle systems in order to determine the phase boundary to a long-range ordered, coherently-moving state. The microscopic origin of the collision rule offers the opportunity to quantitatively scrutinize the predictions of kinetic theory for propelled particle systems through direct comparison with multi-particle simulations.

preprint2012arXiv

Evolutionary and Population Dynamics: A Coupled Approach

We study the interplay of population growth and evolutionary dynamics using a stochastic model based on birth and death events. In contrast to the common assumption of an independent population size, evolution can be strongly affected by population dynamics in general. Especially for fast reproducing microbes which are subject to selection, both types of dynamics are often closely intertwined. We illustrate this by considering different growth scenarios. Depending on whether microbes die or stop to reproduce (dormancy), qualitatively different behaviors emerge. For cooperating bacteria, a permanent increase of costly cooperation can occur. Even if not permanent, cooperation can still increase transiently due to demographic fluctuations. We validate our analysis via stochastic simulations and analytic calculations. In particular, we derive a condition for an increase in the level of cooperation.

preprint2012arXiv

Growth dynamics and the evolution of cooperation in microbial populations

Microbes providing public goods are widespread in nature despite running the risk of being exploited by free-riders. However, the precise ecological factors supporting cooperation are still puzzling. Following recent experiments, we consider the role of population growth and the repetitive fragmentation of populations into new colonies mimicking simple microbial life-cycles. Individual-based modeling reveals that demographic fluctuations, which lead to a large variance in the composition of colonies, promote cooperation. Biased by population dynamics these fluctuations result in two qualitatively distinct regimes of robust cooperation under repetitive fragmentation into groups. First, if the level of cooperation exceeds a threshold, cooperators will take over the whole population. Second, cooperators can also emerge from a single mutant leading to a robust coexistence between cooperators and free-riders. We find frequency and size of population bottlenecks, and growth dynamics to be the major ecological factors determining the regimes and thereby the evolutionary pathway towards cooperation.

preprint2012arXiv

Microtubule Length-Regulation by Molecular Motors

Length-regulation of microtubules (MTs) is essential for many cellular processes. Molecular motors like kinesin 8, which move along MTs and also act as depolymerases, are known as key players in MT dynamics. However, the regulatory mechanisms of length control remain elusive. Here, we investigate a stochastic model accounting for the interplay between polymerization kinetics and motor-induced depolymerization. We determine the dependence of MT length and variance on rate constants and motor concentration. Moreover, our analyses reveal how collective phenomena lead to a well-defined MT length.

preprint2012arXiv

Nucleation-induced transition to collective motion in active systems

While the existence of polar ordered states in active systems is well established, the dynamics of the self-assembly processes are still elusive. We study a lattice gas model of self-propelled elongated particles interacting through excluded volume and alignment interactions, which shows a phase transition from an isotropic to a polar ordered state. By analyzing the ordering process we find that the transition is driven by the formation of a critical nucleation cluster and a subsequent coarsening process. Moreover, the time to establish a polar ordered state shows a power-law divergence.

preprint2012arXiv

On the Validity of the Law of Mass Action in Three-Dimensional Coagulation Processes

Diffusion limited reactions are studied in detail on the classical coalescing process. We demonstrate how, with the aid of a recent renormalization group approach, fluctuations can be integrated systematically. We thereby obtain an exact relation between the microscopic physics (lattice structure, particle shape and size) and the macroscopic decay rate in the law of mass action. Moreover, we find a strong violation of the law of mass action. The corresponding term in the kinetic equations originates in long wavelength fluctuations and is a universal function of the macroscopic decay rate.

preprint2011arXiv

Crowding of molecular motors determines microtubule depolymerization

Assembly and disassembly dynamics of microtubules (MTs) is tightly controlled by MT associated proteins. Here, we investigate how plus-end-directed depolymerases of the kinesin-8 family regulate MT depolymerization dynamics. Employing an individual-based model, we reproduce experimental findings. Moreover, crowding is identified as the key regulatory mechanism of depolymerization dynamics. Our analysis gives two qualitatively distinct regimes. For motor densities above a particular threshold, a macroscopic traffic jam emerges at the plus-end and the MT dynamics become independent of the motor concentration. Below this threshold, microscopic traffic jams at the tip arise which cancel out the effect of the depolymerization kinetics such that the depolymerization speed is solely determined by the motor density. Because this density changes over the MT length, length-dependent regulation is possible. Remarkably, motor cooperativity does not affect the depolymerization speed but only the end-residence time of depolymerases.

preprint2011arXiv

Current reversal and exclusion processes with history-dependent random walks

A class of exclusion processes in which particles perform history-dependent random walks is introduced, stimulated by dynamic phenomena in some biological and artificial systems. The particles locally interact with the underlying substrate by breaking and reforming lattice bonds. We determine the steady-state current on a ring, and find current-reversal as a function of particle density. This phenomenon is attributed to the non-local interaction between the walkers through their trails, which originates from strong correlations between the dynamics of the particles and the lattice. We rationalize our findings within an effective description in terms of quasi-particles which we call front barriers. Our analytical results are complemented by stochastic simulations.

preprint2011arXiv

Driven Transport on Parallel Lanes with Particle Exclusion and Obstruction

We investigate a driven two-channel system where particles on different lanes mutually obstruct each others motion extending an earlier model by Popkov and Peschel [1]. This obstruction may occur in biological contexts due to steric hinderance where motor proteins carry cargos by `walking' on microtubules. Similarly, the model serves as a description for classical spin transport where charged particles with internal states move unidirectionally on a lattice. Three regimes of qualitatively different behavior are identified depending on the strength of coupling between the lanes. For small and large coupling strengths the model can be mapped to a one-channel problem, whereas a new and rich phase behavior emerges for intermediate ones. We derive an approximate but quantitatively accurate theoretical description in terms of a one-site cluster approximation, and obtain insight into the phase behavior through the current-density relations combined with an extremal-current principle. Our results are confirmed by stochastic simulations. [1] V. Popkov and I. Peschel, Phys. Rev. E 64, 026126 (2001).

preprint2011arXiv

Frozen steady states in active systems

Even simple active systems can show a plethora of intriguing phenomena and often we find complexity were we would have expected simplicity. One striking example is the occurrence of a quiescent or absorbing state with frozen fluctuations that at first sight seems to be impossible for active matter driven by the incessant input of energy. While such states were reported for externally driven systems through macroscopic shear or agitation, the investigation of frozen active states in inherently active systems like cytoskeletal suspensions or active gels is still at large. Using high density motility assay experiments, we demonstrate that frozen steady states can arise in active systems if active transport is coupled to growth processes. The experiments are complemented by agent-based simulations which identify the coupling between self-organization, growth and mechanical properties to be responsible for the pattern formation process.

preprint2011arXiv

Range Expansion with Mutation and Selection: Dynamical Phase Transition in a Two-Species Eden Model

The colonization of unoccupied territory by invading species, known as range expansion, is a spatially heterogeneous non-equilibrium growth process. We introduce a two-species Eden growth model to analyze the interplay between uni-directional (irreversible) mutations and selection at the expanding front. While the evolutionary dynamics leads to coalescence of both wild-type and mutant clusters, the non-homogeneous advance of the colony results in a rough front. We show that roughening and domain dynamics are strongly coupled, resulting in qualitatively altered bulk and front properties. For beneficial mutations the front is quickly taken over by mutants and growth proceeds Eden-like. In contrast, if mutants grow slower than wild-types, there is an antagonism between selection pressure against mutants and growth by merging of mutant domains with an ensuing absorbing state phase transition to an all-mutant front. We find that surface roughening has a marked effect on the critical properties of the absorbing state phase transition. While reference models, which keep the expanding front flat, exhibit directed percolation critical behavior, the exponents of the two-species Eden model strongly deviate from it. In turn, the mutation-selection process induces an increased surface roughness with exponents distinct from that of the classical Eden model.

preprint2011arXiv

Social dynamics with peer support on heterogeneous networks: The "mafia model"

Human behavior often exhibit a scheme in which individuals adopt indifferent, neutral, or radical positions on a given topic. The mechanisms leading to community formation are strongly related with social pressure and the topology of the contact network. Here, we discuss an approach to model social behavior which accounts for the protection by alike peers proportional to their relative abundance in the closest neighborhood. We explore the ensuing non-linear dynamics emphasizing the role of the specific structure of the social network, modeled by scale-free graphs. We find that both coexistence of opinions and consensus on the default position are possible stationary states of the model. In particular, we show how these states critically depend on the heterogeneity of the social network and the specific distribution of external control elements.

preprint2010arXiv

Buckling of stiff polymer rings in weak spherical confinement

Confinement is a versatile and well-established tool to study the properties of polymers either to understand biological processes or to develop new nano-biomaterials. We investigate the conformations of a semiflexible polymer ring in weak spherical confinement imposed by an impenetrable shell. We develop an analytic argument for the dominating polymer trajectory depending on polymer flexibility considering elastic and entropic contributions. Monte Carlo simulations are performed to assess polymer ring conformations in probability densities and by the shape measures asphericity and nature of asphericity. Comparison of the analytic argument with the mean asphericity and the mean nature of asphericity confirm our reasoning to explain polymer ring conformations in the stiff regime, where elastic response prevails.

preprint2010arXiv

Coexistence in a One-Dimensional Cyclic Dominance Process

Cyclic (rock-paper-scissors-type) population models serve to mimic complex species interactions. Focusing on a paradigmatic three-species model with mutations in one dimension, we observe an interplay between equilibrium and non-equilibrium processes in the stationary state. We exploit these insights to obtain asymptotically exact descriptions of the emerging reactive steady state in the regimes of high and low mutation rates. The results are compared to stochastic lattice simulations. Our methods and findings are potentially relevant for the spatio-temporal evolution of other non-equilibrium stochastic processes.

preprint2010arXiv

Entropy production of cyclic population dynamics

Entropy serves as a central observable in equilibrium thermodynamics. However, many biological and ecological systems operate far from thermal equilibrium. Here we show that entropy production can characterize the behavior of such nonequilibrium systems. To this end we calculate the entropy production for a population model that displays nonequilibrium behavior resulting from cyclic competition. At a critical point the dynamics exhibits a transition from large, limit-cycle like oscillations to small, erratic oscillations. We show that the entropy production peaks very close to the critical point and tends to zero upon deviating from it. We further provide analytical methods for computing the entropy production which agree excellently with numerical simulations.

preprint2010arXiv

Error thresholds for self- and cross-specific enzymatic replication

The information content of a non-enzymatic self-replicator is limited by Eigen's error threshold. Presumably, enzymatic replication can maintain higher complexity, but in a competitive environment such a replicator is faced with two problems related to its twofold role as enzyme and substrate: as enzyme, it should replicate itself rather than wastefully copy non-functional substrates, and as substrate it should preferably be replicated by superior enzymes instead of less-efficient mutants. Because specific recognition can enforce these propensities, we thoroughly analyze an idealized quasispecies model for enzymatic replication, with replication rates that are either a decreasing (self-specific) or increasing (cross-specific) function of the Hamming distance between the recognition or "tag" sequences of enzyme and substrate. We find that very weak self-specificity suffices to localize a population about a master sequence and thus to preserve its information, while simultaneous localization about complementary sequences in the cross-specific case is more challenging. A surprising result is that stronger specificity constraints allow longer recognition sequences, because the populations are better localized. Extrapolating from experimental data, we obtain rough quantitative estimates for the maximal length of the recognition or tag sequence that can be used to reliably discriminate appropriate and infeasible enzymes and substrates, respectively.

preprint2010arXiv

Evolutionary game theory in growing populations

Existing theoretical models of evolution focus on the relative fitness advantages of different mutants in a population while the dynamic behavior of the population size is mostly left unconsidered. We here present a generic stochastic model which combines the growth dynamics of the population and its internal evolution. Our model thereby accounts for the fact that both evolutionary and growth dynamics are based on individual reproduction events and hence are highly coupled and stochastic in nature. We exemplify our approach by studying the dilemma of cooperation in growing populations and show that genuinely stochastic events can ease the dilemma by leading to a transient but robust increase in cooperation

preprint2010arXiv

Excluded volume effects on semiflexible ring polymers

Two-dimensional semiflexible polymer rings are studied both by imaging circular DNA adsorbed on a mica surface and by Monte Carlo simulations of phantom polymers as well as of polymers with finite thickness. Comparison of size and shape of the different models over the full range of flexibilities shows that excluded volume caused by finite thickness induces an anisotropic increase of the main axes of the conformations, a change of shape, accomplished by an enhanced correlation along the contour.

preprint2010arXiv

Longitudinal Response of Confined Semiflexible Polymers

The longitudinal response of single semiflexible polymers to sudden changes in externally applied forces is known to be controlled by the propagation and relaxation of backbone tension. Under many experimental circumstances, realized, e.g., in nano-fluidic devices or in polymeric networks or solutions, these polymers are effectively confined in a channel- or tube-like geometry. By means of heuristic scaling laws and rigorous analytical theory, we analyze the tension dynamics of confined semiflexible polymers for various generic experimental setups. It turns out that in contrast to the well-known linear response, the influence of confinement on the non-linear dynamics can largely be described as that of an effective prestress. We also study the free relaxation of an initially confined chain, finding a surprising superlinear t^(9/8) growth law for the change in end-to-end distance at short times.

preprint2010arXiv

Persistent memory for a Brownian walker in a random array of obstacles

We show that for particles performing Brownian motion in a frozen array of scatterers long-time correlations emerge in the mean-square displacement. Defining the velocity autocorrelation function (VACF) via the second time-derivative of the mean-square displacement, power-law tails govern the long-time dynamics similar to the case of ballistic motion. The physical origin of the persistent memory is due to repeated encounters with the same obstacle which occurs naturally in Brownian dynamics without involving other scattering centers. This observation suggests that in this case the VACF exhibits these anomalies already at first order in the scattering density. Here we provide an analytic solution for the dynamics of a tracer for a dilute planar Lorentz gas and compare our results to computer simulations. Our result support the idea that quenched disorder provides a generic mechanism for persistent correlations irrespective of the microdynamics of the tracer particle.

preprint2010arXiv

Predictive Modeling of Non-Viral Gene Transfer

In non-viral gene delivery, the variance of transgenic expression stems from the low number of plasmids successfully transferred. Here, we experimentally determine Lipofectamine- and PEI-mediated exogenous gene expression distributions from single cell time-lapse analysis. Broad Poisson-like distributions of steady state expression are observed for both transfection agents, when used with synchronized cell lines. At the same time, co-transfection analysis with YFP- and CFP-coding plasmids shows that multiple plasmids are simultaneously expressed, suggesting that plasmids are delivered in correlated units (complexes). We present a mathematical model of transfection, where a stochastic, two-step process is assumed, with the first being the low-probability entry step of complexes into the nucleus, followed by the subsequent release and activation of a small number of plasmids from a delivered complex. This conceptually simple model consistently predicts the observed fraction of transfected cells, the cotransfection ratio and the expression level distribution. It yields the number of efficient plasmids per complex and elucidates the origin of the associated noise, consequently providing a platform for evaluating and improving non-viral vectors.

preprint2010arXiv

Quantitative predictions on auxin-induced polar distribution of PIN proteins during vein formation in leaves

The dynamic patterning of the plant hormone auxin and its efflux facilitator the PIN protein are the key regulator for the spatial and temporal organization of plant development. In particular auxin induces the polar localization of its own efflux facilitator. Due to this positive feedback auxin flow is directed and patterns of auxin and PIN arise. During the earliest stage of vein initiation in leaves auxin accumulates in a single cell in a rim of epidermal cells from which it flows into the ground meristem tissue of the leaf blade. There the localized auxin supply yields the successive polarization of PIN distribution along a strand of cells. We model the auxin and PIN dynamics within cells with a minimal canalization model. Solving the model analytically we uncover an excitable polarization front that triggers a polar distribution of PIN proteins in cells. As polarization fronts may extend to opposing directions from their initiation site we suggest a possible resolution to the puzzling occurrence of bipolar cells, such we offer an explanation for the development of closed, looped veins. Employing non-linear analysis we identify the role of the contributing microscopic processes during polarization. Furthermore, we deduce quantitative predictions on polarization fronts establishing a route to determine the up to now largely unknown kinetic rates of auxin and PIN dynamics.

preprint2010arXiv

The Localization Transition of the Two-Dimensional Lorentz Model

We investigate the dynamics of a single tracer particle performing Brownian motion in a two-dimensional course of randomly distributed hard obstacles. At a certain critical obstacle density, the motion of the tracer becomes anomalous over many decades in time, which is rationalized in terms of an underlying percolation transition of the void space. In the vicinity of this critical density the dynamics follows the anomalous one up to a crossover time scale where the motion becomes either diffusive or localized. We analyze the scaling behavior of the time-dependent diffusion coefficient D(t) including corrections to scaling. Away from the critical density, D(t) exhibits universal hydrodynamic long-time tails both in the diffusive as well as in the localized phase.

preprint2009arXiv

Scaling and universality in coupled driven diffusive models

Inspired by the physics of magnetohydrodynamics (MHD) a simplified coupled Burgers-like model in one dimension (1d), a generalization of the Burgers model to coupled degrees of freedom, is proposed to describe 1dMHD. In addition to MHD, this model serves as a 1d reduced model for driven binary fluid mixtures. Here we have performed a comprehensive study of the universal properties of the generalized d-dimensional version of the reduced model. We employ both analytical and numerical approaches. In particular, we determine the scaling exponents and the amplitude-ratios of the relevant two-point time-dependent correlation functions in the model. We demonstrate that these quantities vary continuously with the amplitude of the noise cross-correlation. Further our numerical studies corroborate the continuous dependence of long wavelength and long time-scale physics of the model on the amplitude of the noise cross-correlations, as found in our analytical studies. We construct and simulate lattice-gas models of coupled degrees of freedom in 1d, belonging to the universality class of our coupled Burgers-like model, which display similar behavior. We use a variety of numerical (Monte-Carlo and Pseudospectral methods) and analytical (Dynamic Renormalization Group, Self-Consistent Mode-Coupling Theory and Functional Renormalization Group) approaches for our work. The results from our different approaches complement one another. Possible realizations of our results in various nonequilibrium models are discussed.

preprint2009arXiv

Statics and Dynamics of the Wormlike Bundle Model

Bundles of filamentous polymers are primary structural components of a broad range of cytoskeletal structures, and their mechanical properties play key roles in cellular functions ranging from locomotion to mechanotransduction and fertilization. We give a detailed derivation of a wormlike bundle model as a generic description for the statics and dynamics of polymer bundles consisting of semiflexible polymers interconnected by crosslinking agents. The elastic degrees of freedom include bending as well as twist deformations of the filaments and shear deformation of the crosslinks. We show that a competition between the elastic properties of the filaments and those of the crosslinks leads to renormalized effective bend and twist rigidities that become mode-number dependent. The strength and character of this dependence is found to vary with bundle architecture, such as the arrangement of filaments in the cross section and pretwist. We discuss two paradigmatic cases of bundle architecture, a uniform arrangement of filaments as found in F-actin bundles and a shell-like architecture as characteristic for microtubules. Each architecture is found to have its own universal ratio of maximal to minimal bending rigidity, independent of the specific type of crosslink induced filament coupling; our predictions are in reasonable agreement with available experimental data for microtubules. Moreover, we analyze the predictions of the wormlike bundle model for experimental observables such as the tangent-tangent correlation function and dynamic response and correlation functions. Finally, we analyze the effect of pretwist (helicity) on the mechanical properties of bundles. We predict that microtubules with different number of protofilaments should have distinct variations in their effective bending rigidity.

preprint2007arXiv

Velocity oscillations in actin-based motility

We present a simple and generic theoretical description of actin-based motility, where polymerization of filaments maintains propulsion. The dynamics is driven by polymerization kinetics at the filaments' free ends, crosslinking of the actin network, attachment and detachment of filaments to the obstacle interfaces and entropic forces. We show that spontaneous oscillations in the velocity emerge in a broad range of parameter values, and compare our findings with experiments.

preprint2006arXiv

Floppy modes and non-affine deformations in random fiber networks

We study the elasticity of random fiber networks. Starting from a microscopic picture of the non-affine deformation fields we calculate the macroscopic elastic moduli both in a scaling theory and a self-consistent effective medium theory. By relating non-affinity to the low-energy excitations of the network (``floppy-modes'') we achieve a detailed characterization of the non-affine deformations present in fibrous networks.

preprint2006arXiv

Generic principles of active transport

Nonequilibrium collective motion is ubiquitous in nature and often results in a rich collection of intringuing phenomena, such as the formation of shocks or patterns, subdiffusive kinetics, traffic jams, and nonequilibrium phase transitions. These stochastic many-body features characterize transport processes in biology, soft condensed matter and, possibly, also in nanoscience. Inspired by these applications, a wide class of lattice-gas models has recently been considered. Building on the celebrated {\it totally asymmetric simple exclusion process} (TASEP) and a generalization accounting for the exchanges with a reservoir, we discuss the qualitative and quantitative nonequilibrium properties of these model systems. We specifically analyze the case of a dimeric lattice gas, the transport in the presence of pointwise disorder and along coupled tracks.

preprint2006arXiv

Stochastic effects on biodiversity in cyclic coevolutionary dynamics

Finite-size fluctuations arising in the dynamics of competing populations may have dramatic influence on their fate. As an example, in this article, we investigate a model of three species which dominate each other in a cyclic manner. Although the deterministic approach predicts (neutrally) stable coexistence of all species, for any finite population size, the intrinsic stochasticity unavoidably causes the eventual extinction of two of them.

preprint2005arXiv

Brownian motion: a paradigm of soft matter and biological physics

This is a pedagogical introduction to Brownian motion on the occasion of the 100th anniversary of Einstein's 1905 paper on the subject. After briefly reviewing Einstein's work in its contemporary context, we pursue some lines of further developments and applications in soft condensed matter and biology. Over the last century Brownian motion became promoted from an odd curiosity of marginal scientific interest to a guiding theme pervading all of the modern (live) sciences.

preprint1995arXiv

Critical Dynamics of Magnets

We review our current understanding of the critical dynamics of magnets above and below the transition temperature with focus on the effects due to the dipole--dipole interaction present in all real magnets. Significant progress in our understanding of real ferromagnets in the vicinity of the critical point has been made in the last decade through improved experimental techniques and theoretical advances in taking into account realistic spin-spin interactions. We start our review with a discussion of the theoretical results for the critical dynamics based on recent renormalization group, mode coupling and spin wave theories. A detailed comparison is made of the theory with experimental results obtained by different measuring techniques, such as neutron scattering, hyperfine interaction, muon--spin--resonance, electron--spin--resonance, and magnetic relaxation, in various materials. Furthermore we discuss the effects of dipolar interaction on the critical dynamics of three--dimensional isotropic antiferromagnets and uniaxial ferromagnets. Special attention is also paid to a discussion of the consequences of dipolar anisotropies on the existence of magnetic order and the spin--wave spectrum in two--dimensional ferromagnets and antiferromagnets. We close our review with a formulation of critical dynamics in terms of nonlinear Langevin equations.

preprint1995arXiv

Crossover scaling functions and an extended minimal subtraction scheme

A field theoretic renormalization group method is presented which is capable of dealing with crossover problems associated with a change in the upper critical dimension. The method leads to flow functions for the parameters and coupling constants of the model which depend on the set of parameters which characterize the fixed point landscape of the underlying problem. Similar to Nelson's trajectory integral method any vertex function can be expressed as a line integral along a renormalization group trajectory, which in the field theoretic formulation are given by the characteristics of the corresponding Callan-Symanzik equation. The field theoretic renormalization automatically leads to a separation of the regular and singular parts of all crossover scaling functions. The method is exemplified for the crossover problem in magnetic phase transitions, percolation problems and quantum phase transitions. The broad applicability of the method is emphasized.

preprint1994arXiv

Two-Loop Renormalization Group Analysis of the Burgers-Kardar-Parisi-Zhang Equation

A systematic analysis of the Burgers--Kardar--Parisi--Zhang equation in $d+1$ dimensions by dynamic renormalization group theory is described. The fixed points and exponents are calculated to two--loop order. We use the dimensional regularization scheme, carefully keeping the full $d$ dependence originating from the angular parts of the loop integrals. For dimensions less than $d_c=2$ we find a strong--coupling fixed point, which diverges at $d=2$, indicating that there is non--perturbative strong--coupling behavior for all $d \geq 2$. At $d=1$ our method yields the identical fixed point as in the one--loop approximation, and the two--loop contributions to the scaling functions are non--singular. For $d>2$ dimensions, there is no finite strong--coupling fixed point. In the framework of a $2+ε$ expansion, we find the dynamic exponent corresponding to the unstable fixed point, which describes the non--equilibrium roughening transition, to be $z = 2 + {\cal O} (ε^3)$, in agreement with a recent scaling argument by Doty and Kosterlitz. Similarly, our result for the correlation length exponent at the transition is $1/ν= ε+ {\cal O} (ε^3)$. For the smooth phase, some aspects of the crossover from Gaussian to critical behavior are discussed.