Velocity fluctuations of noisy reaction fronts propagating into a metastable state: testing theory in stochastic simulations

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.

Pattern Formation of Glioma Cells: Effects of Adhesion

We investigate clustering of malignant glioma cells. emph{In vitro} experiments in collagen gels identified a cell line that formed clusters in a region of low cell density, whereas a very similar cell line (which lacks an important mutation) did not cluster significantly. We hypothesize that the mutation affects the strength of cell-cell adhesion. We investigate this effect in a new experiment, which follows the clustering dynamics of glioma cells on a surface. We interpret our results in terms of a stochastic model and identify two mechanisms of clustering. First, there is a critical value of the strength of adhesion; above the threshold, large clusters grow from a homogeneous suspension of cells; below it, the system remains homogeneous, similarly to the ordinary phase separation. Second, when cells form a cluster, we have evidence that they increase their proliferation rate. We have successfully reproduced the experimental findings and found that both mechanisms are crucial for cluster formation and growth.

A generalized Cahn-Hilliard equation for biological applications

Recently we considered a stochastic discrete model which describes fronts of cells invading a wound cite{KSS}. In the model cells can move, proliferate, and experience cell-cell adhesion. In this work we focus on a continuum description of this phenomenon by means of a generalized Cahn-Hilliard equation (GCH) with a proliferation term. As in the discrete model, there are two interesting regimes. For subcritical adhesion, there are propagating pulled fronts, similarly to those of Fisher-Kolmogorov equation. The problem of front velocity selection is examined, and our theoretical predictions are in a good agreement with a numerical solution of the GCH equation. For supercritical adhesion, there is a nontrivial transient behavior, where density profile exhibits a secondary peak. To analyze this regime, we investigated relaxation dynamics for the Cahn-Hilliard equation without proliferation. We found that the relaxation process exhibits self-similar behavior. The results of continuum and discrete models are in a good agreement with each other for the different regimes we analyzed.

Hydrodynamics of fluid-solid coexistence in dense shear granular flow

We consider dense rapid shear flow of inelastically colliding hard disks. Navier-Stokes granular hydrodynamics is applied accounting for the recent finding cite{Luding,Khain} that shear viscosity diverges at a lower density than the rest of constitutive relations. New interpolation formulas for constitutive relations between dilute and dense cases are proposed and justified in molecular dynamics (MD) simulations. A linear stability analysis of the uniform shear flow is performed and the full phase diagram is presented. It is shown that when the inelasticity of particle collision becomes large enough, the uniform sheared flow gives way to a two-phase flow, where a dense solid-like striped cluster is surrounded by two fluid layers. The results of the analysis are verified in event-driven MD simulations, and a good agreement is observed.