Search

The importance of a strict quarantine has been widely debated during the COVID-19 epidemic even from the purely epidemiological point of view. One argument against strict lockdown measures is that once the strict quarantine is lifted, the epidemic comes back, and so the cumulative number of infected individuals during the entire epidemic will stay the same. We consider an SIR model on a network and follow the disease dynamics, modeling the phases of quarantine by changing the node degree distribution. We show that the system reaches different steady states based on the history: the outcome of the epidemic is path-dependent despite the same final node degree distribution. The results indicate that two-phase route to the final node degree distribution (a strict phase followed by a soft phase) are always better than one phase (the same soft one) unless all the individuals have the same number of connections at the end (the same degree); in the latter case, the overall number of infected is indeed history-independent. The modeling also suggests that the optimal procedure of lifting the quarantine consists of releasing nodes in the order of their degree - highest first., Comment: 6 pages, 4 figures, accepted to EPL (Europhysics Letters)

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_f\sim 1/N$ of the front can be calculated, in the leading order, via a perturbation theory in $1/N \ll 1$, where $N\gg 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., Comment: 10 pages, 7 figures

Continuum models of epidemics do not take into account the underlying microscopic network structure of social connections. This drawback becomes extreme during quarantine when most people dramatically decrease their number of social interactions, while others (like cashiers in grocery stores) continue maintaining hundreds of contacts per day. We formulate a two-level model of quarantine. On a microscopic level, we model a single neighborhood assuming a star-network structure. On a mesoscopic level, the neighborhoods are placed on a two-dimensional lattice with nearest neighbors interactions. The modeling results are compared with the COVID-19 data for several counties in Michigan (USA) and the phase diagram of parameters is identified., Comment: 5 pages, 6 figures (figure 3 has two panels)

The granular Leidenfrost effect (B. Meerson et al, Phys. Rev. Lett. {\bf 91}, 024301 (2003), P. Eshuis et al, Phys. Rev. Lett. {\bf 95}, 258001 (2005)) is the levitation of a mass of granular matter when a wall below the grains is vibrated giving rise to a hot granular gas below the cluster. We find by simulation that for a range of parameters the system is bistable: the levitated cluster can occasionally break and give rise to two clusters and a hot granular gas above and below. We use techniques from the theory of rare events to compute the mean transition time for breaking to occur. This requires the introduction of a two-component reaction coordinate.

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., Comment: 5 pages, 5 figures

We consider population dynamics on a network of patches, each of which has a the same local dynamics, with different population scales (carrying capacities). It is reasonable to assume that if the patches are coupled by very fast migration the whole system will look like an individual patch with a large effective carrying capacity. This is called a "well-mixed" system. We show that, in general, it is not true that the well-mixed system has the same dynamics as each local patch. Different global dynamics can emerge from coupling, and usually must be figured out for each individual case. We give a general condition which must be satisfied for well-mixed systems to have the same dynamics as the constituent patches., Comment: 4 pages

Many populations in nature are fragmented: they consist of local populations occupying separate patches. A local population is prone to extinction due to the shot noise of birth and death processes. A migrating population from another patch can dramatically delay the extinction. What is the optimal migration rate that minimizes the extinction risk of the whole population? Here we answer this question for a connected network of model habitat patches with different carrying capacities., Comment: 7 pages, 3 figures, accepted for publication in PRL, appendix contains supplementary material

Thermal wall is a convenient idealization of a rapidly vibrating plate used for vibrofluidization of granular materials. The objective of this work is to incorporate the Knudsen temperature jump at thermal wall in the Navier-Stokes hydrodynamic modeling of dilute granular gases of monodisperse particles that collide nearly elastically. The Knudsen temperature jump manifests itself as an additional term, proportional to the temperature gradient, in the boundary condition for the temperature. Up to a numerical pre-factor of order unity, this term is known from kinetic theory of elastic gases. We determine the previously unknown numerical pre-factor by measuring, in a series of molecular dynamics (MD) simulations, steady-state temperature profiles of a gas of elastically colliding hard disks, confined between two thermal walls kept at different temperatures, and comparing the results with the predictions of a hydrodynamic calculation employing the modified boundary condition. The modified boundary condition is then applied, without any adjustable parameters, to a hydrodynamic calculation of the temperature profile of a gas of inelastic hard disks driven by a thermal wall. We find the hydrodynamic prediction to be in very good agreement with MD simulations of the same system. The results of this work pave the way to a more accurate hydrodynamic modeling of driven granular gases., Comment: 7 pages, 3 figures

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., Comment: 6 pages, 6 figures

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., Comment: 11 figures, submitted to PRE

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.

We present a stochastic model which describes fronts of cells invading a wound. In the model cells can move, proliferate, and experience cell-cell adhesion. We find several qualitatively different regimes of front motion and analyze the transitions between them. Above a critical value of adhesion and for small proliferation large isolated clusters are formed ahead of the front. This is mapped onto the well-known ferromagnetic phase transition in the Ising model. For large adhesion, and larger proliferation the clusters become connected (at some fixed time). For adhesion below the critical value the results are similar to our previous work which neglected adhesion. The results are compared with experiments, and possible directions of future work are proposed., Comment: to appear in Journal of Statistical Physics

We investigate shear-induced crystallization in a very dense flow of mono-disperse inelastic hard spheres. We consider a steady plane Couette flow under constant pressure and neglect gravity. We assume that the granular density is greater than the melting point of the equilibrium phase diagram of elastic hard spheres. We employ a Navier-Stokes hydrodynamics with constitutive relations all of which (except the shear viscosity) diverge at the crystal packing density, while the shear viscosity diverges at a smaller density. The phase diagram of the steady flow is described by three parameters: an effective Mach number, a scaled energy loss parameter, and an integer number m: the number of half-oscillations in a mechanical analogy that appears in this problem. In a steady shear flow the viscous heating is balanced by energy dissipation via inelastic collisions. This balance can have different forms, producing either a uniform shear flow or a variety of more complicated, nonlinear density, velocity and temperature profiles. In particular, the model predicts a variety of multi-layer two-phase steady shear flows with sharp interphase boundaries. Such a flow may include a few zero-shear (solid-like) layers, each of which moving as a whole, separated by fluid-like regions. As we are dealing with a hard sphere model, the granulate is fluidized within the "solid" layers: the granular temperature is non-zero there, and there is energy flow through the boundaries of the "solid" layers. A linear stability analysis of the uniform steady shear flow is performed, and a plausible bifurcation diagram of the system, for a fixed m, is suggested. The problem of selection of m remains open., Comment: 11 pages, 7 eps figures, to appear in PRE

We present a discrete stochastic model which represents many of the salient features of the biological process of wound healing. The model describes fronts of cells invading a wound. We have numerical results in one and two dimensions. In one dimension we can give analytic results for the front speed as a power series expansion in a parameter, p, that gives the relative size of proliferation and diffusion processes for the invading cells. In two dimensions the model becomes the Eden model for p near 1. In both one and two dimensions for small p, front propagation for this model should approach that of the Fisher-Kolmogorov equation. However, as in other cases, this discrete model approaches Fisher-Kolmogorov behavior slowly., Comment: 16 pages, 7 figures

In this work, we study the in-vitro dynamics of the most malignant form of the primary brain tumor: Glioblastoma Multiforme. Typically, the growing tumor consists of the inner dense proliferating zone and the outer less dense invasive region. Experiments with different types of cells show qualitatively different behavior. Wild-type cells invade a spherically symmetric manner, but mutant cells are organized in tenuous branches. We formulate a model for this sort of growth using two coupled reaction-diffusion equations for the cell and nutrient concentrations. When the ratio of the nutrient and cell diffusion coefficients exceeds some critical value, the plane propagating front becomes unstable with respect to transversal perturbations. The instability threshold and the full phase-plane diagram in the parameter space are determined. The results are in a good agreement with experimental findings for the two types of cells., Comment: 4 pages, 4 figures

Equations of granular hydrostatics are used to compute the phase diagram of the recently discovered van der Waals-like phase separation in a driven granular gas. The model two-dimensional system consists of smooth hard disks in a rectangular box, colliding inelastically with each other and driven by a "thermal" wall at zero gravity. The spinodal line and the critical point of the phase separation are determined. Close to the critical point the spinodal and binodal (coexistence) lines are determined analytically. Effects of finite size of the confining box in the direction parallel to the thermal wall are investigated. These include suppression of the phase separation by heat conduction in the lateral direction and a change from supercritical to subcritical bifurcation., Comment: 14 pages, 8 figures, to appear in PRE

We discovered an oscillatory instability in a system of inelastically colliding hard spheres, driven by two opposite "thermal" walls at zero gravity. The instability, predicted by a linear stability analysis of the equations of granular hydrodynamics, occurs when the inelasticity of particle collisions exceeds a critical value. Molecular dynamic simulations support the theory and show a stripe-shaped cluster moving back and forth in the middle of the box away from the driving walls. The oscillations are irregular but have a single dominating frequency that is close to the frequency at the instability onset, predicted from hydrodynamics., Comment: 7 pages, 4 figures, to appear in Europhysics Letters

The Navier-Stokes granular hydrodynamics is employed for determining the threshold of thermal convection in an infinite horizontal layer of granular gas. The dependence of the convection threshold, in terms of the inelasticity of particle collisions, on the Froude and Knudsen numbers is found. A simple necessary condition for convection is formulated in terms of the Schwarzschild's criterion, well-known in thermal convection of (compressible) classical fluids. The morphology of convection cells at the onset is determined. At large Froude numbers, the Froude number drops out of the problem. As the Froude number goes to zero, the convection instability turns into a recently discovered phase separation instability., Comment: 6 pages, 6 figures. An extended version. A simple and universal necessary criterion for convection presented

Symmetry-breaking instability of a laterally uniform granular cluster (strip state) in a prototypical driven granular gas is investigated. The system consists of smooth hard disks in a two-dimensional box, colliding inelastically with each other and driven, at zero gravity, by a "thermal" wall. The limit of nearly elastic particle collisions is considered, and granular hydrodynamics with the Jenkins-Richman constitutive relations is employed. The hydrodynamic problem is completely described by two scaled parameters and the aspect ratio of the box. Marginal stability analysis predicts a spontaneous symmetry breaking instability of the strip state, similar to that predicted recently for a different set of constitutive relations. If the system is big enough, the marginal stability curve becomes independent of the details of the boundary condition at the driving wall. In this regime, the density perturbation is exponentially localized at the elastic wall opposite to the thermal wall. The short- and long-wavelength asymptotics of the marginal stability curves are obtained analytically in the dilute limit. The physics of the symmetry-breaking instability is discussed., Comment: 11 pages, 14 figures

We investigate parametric autoresonance: a persisting phase locking which occurs when the driving frequency of a parametrically excited nonlinear oscillator slowly varies with time. In this regime, the resonant excitation is continuous and unarrested by the oscillator nonlinearity. The system has three characteristic time scales, the fastest one corresponding to the natural frequency of the oscillator. We perform averaging over the fastest time scale and analyze the reduced set of equations analytically and numerically. Analytical results are obtained by exploiting the scale separation between the two remaining time scales which enables one to use the adiabatic invariant of the perturbed nonlinear motion., Comment: 8 pages, 10 figures