No Arabic abstract
We analyse four approaches to elimination of a fast variable, which are applicable to systems like passive Brownian particles: (i) moment formalism, (ii) corresponding cumulant formalism, (iii) Hermite function basis, (iv) formal `cumulants for the Hermit function basis. The accuracy and its strong order are assessed. The applicability and performance of two first approaches are also demonstrated for active Brownian particles.
We analyze two approaches to elimination of a fast variable (velocity) in stochastic systems: moment and cumulant formalisms. With these approaches, we obtain the corresponding Smoluchovski-type equations, which contain only the coordinate/phase variable. The adiabatic elimination of velocity in terms of cumulants and moments requires the first three elements. However, for the case of small inertia, the corrected Smoluchowski equation in terms of moments requires five elements, while in terms of cumulants the same first three elements are sufficient. Compared to the method based on the expansion of the velocity distribution in Hermite functions, the considered approaches have comparable efficiency, but do not require individual mathematical preparation for the case of active Brownian particles, where one has to construct a new basis of eigenfunctions instead of the Hermite ones.
Freidlin-Wentzell theory of large deviations can be used to compute the likelihood of extreme or rare events in stochastic dynamical systems via the solution of an optimization problem. The approach gives exponential estimates that often need to be refined via calculation of a prefactor. Here it is shown how to perform these computations in practice. Specifically, sharp asymptotic estimates are derived for expectations, probabilities, and mean first passage times in a form that is geared towards numerical purposes: they require solving well-posed matrix Riccati equations involving the minimizer of the Freidlin-Wentzell action as input, either forward or backward in time with appropriate initial or final conditions tailored to the estimate at hand. The usefulness of our approach is illustrated on several examples. In particular, invariant measure probabilities and mean first passage times are calculated in models involving stochastic partial differential equations of reaction-advection-diffusion type.
Avalanching systems are treated analytically using the renormalization group (in the self-organized-criticality regime) or mean-field approximation, respectively. The latter describes the state in terms of the mean number of active and passive sites, without addressing the inhomogeneity in their distribution. This paper goes one step further by proposing a kinetic description of avalanching systems making use of the distribution function for clusters of active sites. We illustrate application of the kinetic formalism to a model proposed for the description of the avalanching processes in the reconnecting current sheet of the Earth magnetosphere.
We study the dynamics of predator-prey systems where prey are confined to a single region of space and where predators move randomly according to a power-law (Levy) dispersal kernel. Site fidelity, an important feature of animal behaviour, is incorporated in the model through a stochastic resetting dynamics of the predators to the prey patch. We solve in the long time limit the rate equations of Lotka-Volterra type that describe the evolution of the two species densities. Fixing the demographic parameters and the Levy exponent, the total population of predators can be maximized for a certain value of the resetting rate. This optimal value achieves a compromise between over-exploitation and under-utilization of the habitat. Similarly, at fixed resetting rate, there exists a Levy exponent which is optimal regarding predator abundance. These findings are supported by 2D stochastic simulations and show that the combined effects of diffusion and resetting can broadly extend the region of species coexistence in ecosystems facing resources scarcity.
It is shown that a small system in thermodynamic equilibrium with a finite thermostat can have a q-exponential probability distribution which closely depends on the energy nonextensivity and the particle number of the thermostat. The distribution function will reduce to the exponential one at the thermodynamic limit. However, the nonextensivity of the system should not be neglected.