Environmental vs demographic variability in stochastic predator-prey models

In contrast to the neutral population cycles of the deterministic mean-field Lotka--Volterra rate equations, including spatial structure and stochastic noise in models for predator-prey interactions yields complex spatio-temporal structures associated with long-lived erratic population oscillations. Environmental variability in the form of quenched spatial randomness in the predation rates results in more localized activity patches. Population fluctuations in rare favorable regions in turn cause a remarkable increase in the asymptotic densities of both predators and prey. Very intriguing features are found when variable interaction rates are affixed to individual particles rather than lattice sites. Stochastic dynamics with demographic variability in conjunction with inheritable predation efficiencies generate non-trivial time evolution for the predation rate distributions, yet with overall essentially neutral optimization.Comment: 28 pages, 10 figures, Proceedings paper of the STATPHYS 25 conferenc

Pinning time statistics for vortex lines in disordered environments

We study the pinning dynamics of magnetic flux (vortex) lines in a disordered type-II superconductor. Using numerical simulations of a directed elastic line model, we extract the pinning time distributions of vortex line segments. We compare different model implementations for the disorder in the surrounding medium: discrete, localized pinning potential wells that are either attractive and repulsive or purely attractive, and whose strengths are drawn from a Gaussian distribution; as well as continuous Gaussian random potential landscapes. We find that both schemes yield power law distributions in the pinned phase as predicted by extreme-event statistics, yet they differ significantly in their effective scaling exponents and their short-time behavior.Comment: 7 pages, 5 figures, to appear in Phys. Rev.

Reconstructing the gradient source position from steady-state fluxes to small receptors

Recovering the position of a source from the fluxes of diffusing particles through small receptors allows a biological cell to determine its relative position, spatial localization and guide it to a final target. However, how a source can be recovered from point fluxes remains unclear. Using the Narrow Escape Time approach for an open domain, we compute the diffusion fluxes of Brownian particles generated by a steady-state gradient from a single source through small holes distributed on a surface in two dimensions. {We find that the location of a source can be recovered when there are at least 3 receptors and the source is positioned no further than 10 cell radii away}, but this condition is not necessary in a narrow strip. The present approach provides a computational basis for the first step of direction sensing of a gradient at a single cell level.This work was supported by EPSRC grant no EP/K032208/1. U.D. was supported by a Junior Interdisciplinary Fellowship via Wellcome Trust grant number 105602/Z/14/Z and a Herchel Smith Postdoctoral Fellowship. D.H. team is supported by a FRM grant

Disordered vortex matter out of equilibrium: a Langevin molecular dynamics study

We discuss the use of Langevin molecular dynamics in the investigation of the non-equilibrium properties of disordered vortex matter. Our special focus is set on values of system parameters that are realistic for disordered high-$T_c$ superconductors such as YBCO. Using a discretized elastic line model, we study different aspects of vortices far from thermal equilibrium. On the one hand we investigate steady-state properties of driven magnetic flux lines in a disordered environment, namely the current-voltage characteristics, the gyration radius, and the pinning time statistics. On the other hand we study the complex relaxation processes and glassy-like dynamics that emerge in type-II superconductors due to the intricate competition between the long-range vortex-vortex repulsion and flux pinning due to randomly placed point defects. To this end we consider different types of sudden perturbations: temperature, magnetic field, and external current quenches.Comment: 15 pages, 7 figures, to appear in Molecular Simulation for a special issue on "Nonequilibrium Systems

Relaxation dynamics of vortex lines in disordered type-II superconductors following magnetic field and temperature quenches

We study the effects of rapid temperature and magnetic field changes on the non-equilibrium relaxation dynamics of magnetic vortex lines in disordered type-II superconductors by employing an elastic line model and performing Langevin molecular dynamics simulations. In a previously equilibrated system, either the temperature is suddenly changed, or the magnetic field is instantaneously altered which is reflected in adding or removing flux lines to or from the system. The subsequent aging properties are investigated in samples with either randomly distributed point-like or extended columnar defects, which allows to distinguish the complex relaxation features that result from either type of pinning centers. One-time observables such as the radius of gyration and the fraction of pinned line elements are employed to characterize steady-state properties, and two-time correlation functions such as the vortex line height autocorrelations and their mean-square displacement are analyzed to study the non-linear stochastic relaxation dynamics in the aging regime.Comment: 17 pages, 14 figures, to appear in Phys. Rev. E (2015

Mixed analytical-stochastic simulation method for the recovery of a Brownian gradient source from probability fluxes to small windows.

Is it possible to recover the position of a source from the steady-state fluxes of Brownian particles to small absorbing windows located on the boundary of a domain? To address this question, we develop a numerical procedure to avoid tracking Brownian trajectories in the entire infinite space. Instead, we generate particles near the absorbing windows, computed from the analytical expression of the exit probability. When the Brownian particles are generated by a steady-state gradient at a single point, we compute asymptotically the fluxes to small absorbing holes distributed on the boundary of half-space and on a disk in two dimensions, which agree with stochastic simulations. We also derive an expression for the splitting probability between small windows using the matched asymptotic method. Finally, when there are more than two small absorbing windows, we show how to reconstruct the position of the source from the diffusion fluxes. The present approach provides a computational first principle for the mechanism of sensing a gradient of diffusing particles, a ubiquitous problem in cell biology

Critical scaling and aging near the flux-line-depinning transition

We utilize Langevin molecular dynamics simulations to study dynamical critical behavior of magnetic flux lines near the depinning transition in type-II superconductors subject to randomly distributed attractive point defects. We employ a coarse-grained elastic line Hamiltonian for the mutually repulsive vortices and purely relaxational kinetics. In order to infer the stationary-state critical exponents for the continuous non-equilibrium depinning transition at zero temperature T = 0 and at the critical driving current density j_c, we explore two-parameter scaling laws for the flux lines' gyration radius and mean velocity as functions of the two relevant scaling fields T and j - j_c. We also investigate critical aging scaling for the two-time height auto-correlation function in the early-time non-equilibrium relaxation regime to independently measure critical exponents. We provide numerical exponent values for the distinct universality classes of non-interacting and repulsive vortices