Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userVelocity fluctuations of noisy reaction fronts propagating into a metastable state: testing theory in stochastic simulations
128
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
The position of a reaction front, propagating into a metastable state, fluctuates because of the shot noise of reactions and diffusion. A recent theory [B. Meerson, P.V. Sasorov, and Y. Kaplan, Phys. Rev. E 84, 011147 (2011)] gave a closed analytic expression for the front diffusion coefficient in the weak noise limit. Here we test this theory in stochastic simulations involving reacting and diffusing particles on a one-dimensional lattice. We also investigate a small noise-induced systematic shift of the front velocity compared to the prediction from the spatially continuous deterministic reaction-diffusion equation.
rate research
Read More
The empirical velocity of a reaction-diffusion front, propagating into an unstable state, fluctuates because of the shot noises of the reactions and diffusion. Under certain conditions these fluctuations can be described as a diffusion process in the reference frame moving with the average velocity of the front. Here we address pushed fronts, where the front velocity in the deterministic limit is affected by higher-order reactions and is therefore larger than the linear spread velocity. For a subclass of these fronts -- strongly pushed fronts -- the effective diffusion constant $D_fsim 1/N$ of the front can be calculated, in the leading order, via a perturbation theory in $1/N ll 1$, where $Ngg 1$ is the typical number of particles in the transition region. This perturbation theory, however, overestimates the contribution of a few fast particles in the leading edge of the front. We suggest a more consistent calculation by introducing a spatial integration cutoff at a distance beyond which the average number of particles is of order 1. This leads to a non-perturbative correction to $D_f$ which even becomes dominant close to the transition point between the strongly and weakly pushed fronts. At the transition point we obtain a logarithmic correction to the $1/N$ scaling of $D_f$. We also uncover another, and quite surprising, effect of the fast particles in the leading edge of the front. Because of these particles, the position fluctuations of the front can be described as a diffusion process only on very long time intervals with a duration $Delta t gg tau_N$, where $tau_N$ scales as $N$. At intermediate times the position fluctuations of the front are anomalously large and non-diffusive. Our extensive Monte-Carlo simulations of a particular reacting lattice gas model support these conclusions.
In genetic circuits, when the mRNA lifetime is short compared to the cell cycle, proteins are produced in geometrically-distributed bursts, which greatly affects the cellular switching dynamics between different metastable phenotypic states. Motivated by this scenario, we study a general problem of switching or escape in stochastic populations, where influx of particles occurs in groups or bursts, sampled from an arbitrary distribution. The fact that the step size of the influx reaction is a-priori unknown, and in general, may fluctuate in time with a given correlation time and statistics, introduces an additional non-demographic step-size noise into the system. Employing the probability generating function technique in conjunction with Hamiltonian formulation, we are able to map the problem in the leading order onto solving a stationary Hamilton-Jacobi equation. We show that bursty influx exponentially decreases the mean escape time compared to the usual case of single-step influx. In particular, close to bifurcation we find a simple analytical expression for the mean escape time, which solely depends on the mean and variance of the burst-size distribution. Our results are demonstrated on several realistic distributions and compare well with numerical Monte-Carlo simulations.
We discuss similarities and differences between systems of interacting players maximizing their individual payoffs and particles minimizing their interaction energy. Long-run behavior of stochastic dynamics of spatial games with multiple Nash equilibria is analyzed. In particular, we construct an example of a spatial game with three strategies, where stochastic stability of Nash equilibria depends on the number of players and the kind of dynamics.
Surfaces serve as highly efficient catalysts for a vast variety of chemical reactions. Typically, such surface reactions involve billions of molecules which diffuse and react over macroscopic areas. Therefore, stochastic fluctuations are negligible and the reaction rates can be evaluated using rate equations, which are based on the mean-field approximation. However, in case that the surface is partitioned into a large number of disconnected microscopic domains, the number of reactants in each domain becomes small and it strongly fluctuates. This is, in fact, the situation in the interstellar medium, where some crucial reactions take place on the surfaces of microscopic dust grains. In this case rate equations fail and the simulation of surface reactions requires stochastic methods such as the master equation. However, in the case of complex reaction networks, the master equation becomes infeasible because the number of equations proliferates exponentially. To solve this problem, we introduce a stochastic method based on moment equations. In this method the number of equations is dramatically reduced to just one equation for each reactive species and one equation for each reaction. Moreover, the equations can be easily constructed using a diagrammatic approach. We demonstrate the method for a set of astrophysically relevant networks of increasing complexity. It is expected to be applicable in many other contexts in which problems that exhibit analogous structure appear, such as surface catalysis in nanoscale systems, aerosol chemistry in stratospheric clouds and genetic networks in cells.
Field theory tools are applied to analytically study fluctuation and correlation effects in spatially extended stochastic predator-prey systems. In the mean-field rate equation approximation, the classic Lotka-Volterra model is characterized by neutral cycles in phase space, describing undamped oscillations for both predator and prey populations. In contrast, Monte Carlo simulations for stochastic two-species predator-prey reaction systems on regular lattices display complex spatio-temporal structures associated with persistent erratic population oscillations. The Doi-Peliti path integral representation of the master equation for stochastic particle interaction models is utilized to arrive at a field theory action for spatial Lotka-Volterra models in the continuum limit. In the species coexistence phase, a perturbation expansion with respect to the nonlinear predation rate is employed to demonstrate that spatial degrees of freedom and stochastic noise induce instabilities toward structure formation, and to compute the fluctuation corrections for the oscillation frequency and diffusion coefficient. The drastic downward renormalization of the frequency and the enhanced diffusivity are in excellent qualitative agreement with Monte Carlo simulation data.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
