Source author record

Igor M. Sokolov

Igor M. Sokolov appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

44works
17topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

44 published item(s)

preprint2022arXiv

Universal fluctuations and ergodicity of generalized diffusivity on critical percolation clusters

Despite a long history and a clear overall understanding of properties of random walks on an incipient infinite cluster in percolation, some important information on it seems to be missing in the literature. In the present work, we revisit the problem by performing massive numerical simulations for (sub)diffusion of particles on such clusters. Thus, we discuss the shape of the probability density function (PDF) of particles' displacements, and the way it converges to its long-time limiting scaling form. Moreover, we discuss the properties of the mean squared displacement (MSD) of a particle diffusing on the infinite cluster at criticality. This one is known not to be self-averaging. We show that the fluctuations of the MSD in different realizations of the cluster are universal, and discuss the properties of the distribution of these fluctuations. These strong fluctuations coexist with the ergodicity of subdiffusive behavior in the time domain. The dependence of the relative strength of fluctuations in time-averaged MSD on the total trajectory length (total simulation time) is characteristic for diffusion in a percolation system and can be used as an additional test to distinguish this process with disorder-induced memory from processes with otherwise similar behavior, like fractional Brownian motion with the same value of the Hurst exponent.

preprint2021arXiv

Can a patchy model describe the potential spread of West Nile virus in Germany?

In 2018, West Nile Virus (WNV) was detected for the first time in Germany. Since the first detection, 36 human cases and 175 cases in horses and birds are detected. The transmission cycle of West Nile Virus includes birds and mosquitoes and -- as dead-end hosts -- people and horses. Spatial dissemination of the disease is caused by the movements of birds and mosquitoes. While the activity and movement of mosquitoes are depending mainly on temperature, in the birds there is a complex movement pattern caused by local birds and long range dispersal birds. To this end, we have developed a metapopulation network model framework to delineate the potential spatial distribution and spread of WNV across Germany as well as to evaluate the risk throughout our proposed network model. Our model facilitates the interconnection amongst the vector, local birds and long range dispersal birds contact networks. We have assumed different distance dispersal kernels models for the vector and avian populations with the intention to include short and long range dispersal. The model includes spatial variation of mosquito abundance and the movements to resemble the reality.

preprint2021arXiv

Emergence of polarized ideological opinions in multidimensional topic spaces

Opinion polarization is on the rise, causing concerns for the openness of public debates. Additionally, extreme opinions on different topics often show significant correlations. The dynamics leading to these polarized ideological opinions pose a challenge: How can such correlations emerge, without assuming them a priori in the individual preferences or in a preexisting social structure? Here we propose a simple model that qualitatively reproduces ideological opinion states found in survey data, even between rather unrelated, but sufficiently controversial, topics. Inspired by skew coordinate systems recently proposed in natural language processing models, we solidify these intuitions in a formalism of opinions unfolding in a multidimensional space where topics form a non-orthogonal basis. Opinions evolve according to the social interactions among the agents, which are ruled by homophily: two agents sharing similar opinions are more likely to interact. The model features phase transitions between a global consensus, opinion polarization, and ideological states. Interestingly, the ideological phase emerges by relaxing the assumption of an orthogonal basis of the topic space, i.e. if topics thematically overlap. Furthermore, we analytically and numerically show that these transitions are driven by the controversialness of the topics discussed, the more controversial the topics, the more likely are opinion to be correlated. Our findings shed light upon the mechanisms driving the emergence of ideology in the formation of opinions.

preprint2021arXiv

Ticks on the run: A mathematical model of Crimean-Congo Haemorrhagic Fever (CCHF)-key factors for transmission

Crimean-Congo haemorrhagic fever (CCHF) is a tick-borne zoonotic disease caused by the Crimean-Congo hemorrhagic fever virus (CCHFV). Ticks belonging to the genus \textit{Hyalomma} are the main vectors and reservoir for the virus. It is maintained in nature in an endemic vertebrate-tick-vertebrate cycle. CCHFV is prevalent in wide geographical areas including Asia, Africa, South-Eastern Europe and the Middle East. Over the last decade, several outbreaks of CCHFV have been observed in Europe, mainly in Mediterranean countries. Due to the high case/fatality ratio of CCHFV in human sometimes, it is of great importance for public health. Climate change and the invasion of CCHFV vectors in Central Europe suggest that the establishment of the transmission in Central Europe may be possible in future. We developed a compartment-based nonlinear Ordinary Differential Equation (ODE) system to model the disease transmission cycle including blood sucking ticks, livestock and human. Sensitivity analysis of the basic reproduction number $R_0$ shows that decreasing in the tick survival time is an efficient method to eradicate the disease. The model supports us in understanding the influence of different model parameters on the spread of CCHFV. Tick to tick transmission through co-feeding and the CCHFV circulation through trasstadial and transovarial stages are important factors to sustain the disease cycle. The proposed model dynamics are calibrated through an empirical multi-country analysis and multidimensional scaling reveals the disease-parameter sets of different countries burdened with CCHF are different. This necessary information may help us to select most efficient control strategies.

preprint2020arXiv

Brownian motion under non-instantaneous resetting in higher dimensions

We consider Brownian motion under resetting in higher dimensions for the case when the return of the particle to the origin occurs at a constant speed. We investigate the behavior of the probability density function (PDF) and of the mean-squared displacement (MSD) in this process. We study two different resetting protocols: exponentially distributed time intervals between the resetting events (Poissonian resetting) and resetting at fixed time intervals (deterministic resetting). We moreover discuss a general problem of the invariance of the PDF with respect to the return speed, as observed in the one-dimensional system for Poissonian resetting, and show, that this one dimensional situation is the only one in which such an invariance can be found. However, the invariance of the MSD can still be observed in higher dimensions.

preprint2020arXiv

Continuous time random walks under power-law resetting

We study continuous time random walks (CTRW) with power law distribution of waiting times under resetting which brings the walker back to the origin, with a power-law distribution of times between the resetting events. Two situations are considered. Under complete resetting, the CTRW after the resetting event starts anew, with a new waiting time, independent of the prehistory. Under incomplete resetting, the resetting of the coordinate does not influence the waiting time until the next jump. We focus on the behavior of the mean squared displacement (MSD) of the walker from its initial position, on the conditions under which the probability density functions of the walker's displacement show universal behavior, and on this universal behavior itself. We show, that the behavior of the MSD is the same as in the scaled Brownian motion (SBM), being the mean field model of the CTRW. The intermediate asymptotics of the probability density functions (PDF) for CTRW under complete resetting (provided they exist) are also the same as in the corresponding case for SBM. For incomplete resetting, however, the behavior of the PDF for CTRW and SBM is vastly different.

preprint2020arXiv

Inhomogeneous parametric scaling and variable-order fractional diffusion equations

We discuss the derivation and the solutions of integro-differential equations (variable-order time-fractional diffusion equations) following as continuous limits for lattice continuous time random walk schemes with power-law waiting-time probability density functions, whose parameters are position-dependent. We concentrate on subdiffusive cases and discuss two situations as examples: A system consisting of two parts with different exponents of subdiffusion, and a system in which the subdiffusion exponent changes linearly from one end of the interval to another one. In both cases we compare the numerical solutions of generalized master equations describing the process on the lattice with the corresponding solutions of the continuous equations, which follow by exact solution of the corresponding equations in the Laplace domain with subsequent numerical inversion using the Gaver-Stehfest algorithm.

preprint2020arXiv

Modeling echo chambers and polarization dynamics in social networks

Echo chambers and opinion polarization recently quantified in several sociopolitical contexts and across different social media, raise concerns on their potential impact on the spread of misinformation and on openness of debates. Despite increasing efforts, the dynamics leading to the emergence of these phenomena stay unclear. We propose a model that introduces the dynamics of radicalization, as a reinforcing mechanism driving the evolution to extreme opinions from moderate initial conditions. Inspired by empirical findings on social interaction dynamics, we consider agents characterized by heterogeneous activities and homophily. We show that the transition between a global consensus and emerging radicalized states is mostly governed by social influence and by the controversialness of the topic discussed. Compared with empirical data of polarized debates on Twitter, the model qualitatively reproduces the observed relation between users' engagement and opinions, as well as opinion segregation in the interaction network. Our findings shed light on the mechanisms that may lie at the core of the emergence of echo chambers and polarization in social media.

preprint2020arXiv

Resetting processes with non-instanteneous return

We consider a random two-phase process which we call a reset-return one. The particle starts its motion at the origin. The first, displacement, phase corresponds to a stochastic motion of a particle and is finished at a resetting event. The second, return, phase corresponds to the particle's motion towards the origin from the position it attained at the end of the displacement phase. This motion towards the origin takes place according to a given equation of motion. The whole process is a renewal one. We provide general expressions for the stationary probability density function of the particle's position and for the mean hitting time in one dimension. We perform explicit analysis for the Brownian motion during the displacement phase and three different types of the return motion: return at a constant speed, return at a constant acceleration with zero initial speed and return under the action of a harmonic force. We assume that the waiting times for resetting events follow an exponential distribution, or that resetting takes place at a constant pace. For the first two types of return motion and the exponential resetting the stationary probability density function of the particle's position is invariant under return speed (acceleration), while no such invariance is found for deterministic resetting, and for exponential resetting with return under the action of the harmonic force. We discuss necessary conditions for such invariance of the stationary PDF of the positions with respect to the properties of the return process, and demonstrate some additional examples when this invariance does or does not take place.

preprint2020arXiv

Unexpected crossovers in correlated random-diffusivity processes

The passive and active motion of micron-sized tracer particles in crowded liquids and inside living biological cells is ubiquitously characterised by "viscoelastic" anomalous diffusion, in which the increments of the motion feature long-ranged negative and positive correlations. While viscoelastic anomalous diffusion is typically modelled by a Gaussian process with correlated increments, so-called fractional Gaussian noise, an increasing number of systems are reported, in which viscoelastic anomalous diffusion is paired with non-Gaussian displacement distributions. Following recent advances in Brownian yet non-Gaussian diffusion we here introduce and discuss several possible versions of random-diffusivity models with long-ranged correlations. While all these models show a crossover from non-Gaussian to Gaussian distributions beyond some correlation time, their mean squared displacements exhibit strikingly different behaviours: depending on the model crossovers from anomalous to normal diffusion are observed, as well as unexpected dependencies of the effective diffusion coefficient on the correlation exponent. Our observations of the strong non-universality of random-diffusivity viscoelastic anomalous diffusion are important for the analysis of experiments and a better understanding of the physical origins of "viscoelastic yet non-Gaussian" diffusion.

preprint2020arXiv

Voltage distribution in a non-locally but globally electroneutral confined electrolyte medium: applications for nanophysiology

The distribution of voltage in sub-micron cellular domains remains poorly understood. In neurons, the voltage results from the difference in ionic concentrations which are continuously maintained by pumps and exchangers. However, it not clear how electro-neutrality could be maintained by an excess of fast moving positive ions that should be counter balanced by slow diffusing negatively charged proteins. Using the theory of electro-diffusion, we study here the voltage distribution in a generic domain, which consists of two concentric disks (resp. ball) in two (resp. three) dimensions, where a negative charge is fixed in the inner domain. When global but not local electro-neutrality is maintained, we solve the Poisson-Nernst-Planck equation both analytically and numerically in dimension 1 (flat) and 2 (cylindrical) and found that the voltage changes considerably on a spatial scale which is much larger than the Debye screening length, which assumes electro-neutrality. The present result suggests that long-range voltage drop changes are expected in neuronal microcompartments, probably relevant to explain the activation of far away voltage-gated channels located on the surface.

preprint2016arXiv

Anomalous transport in cellular flows: The role of initial conditions and aging

We consider the diffusion-advection problem in two simple cellular flow models (often invoked as examples for subdiffusive tracer's motion) and concentrate on the intermediate time range, in which the tracer's motion indeed may show subdiffusion. We have performed extensive numerical simulations of the systems under different initial conditions, and show that the pure intermediate-time subdiffusion regime is only evident when the particles start at the border between different cells, i.e. at the separatrix, and is less pronounced or absent for other initial conditions. The motion moreover shows quite peculiar aging properties which are also mirrored in the behavior of the time-averaged mean squared displacement for single trajectories. This kind of behavior is due to the complex motion of tracers trapped inside the cell, and is absent in classical models based on continuous time random walks (CTRW) with no dynamics in the trapped state.

preprint2016arXiv

Diffusion of Active Particles With Stochastic Torques Modeled as $α$-Stable Noise

We investigate the stochastic dynamics of an active particle moving at a constant speed under the influence of a fluctuating torque. In our model the angular velocity is generated by a constant torque and random fluctuations described as a Lévy-stable noise. Two situations are investigated. First, we study white Lévy noise where the constant speed and the angular noise generate a persistent motion characterized by the persistence time $τ_D$. At this time scale the crossover from ballistic to normal diffusive behavior is observed. The corresponding diffusion coefficient can be obtained analytically for the whole class of symmetric $α$-stable noises. As typical for models with noise-driven angular dynamics, the diffusion coefficient depends non-monotonously on the angular noise intensity. As second example, we study angular noise as described by an Ornstein-Uhlenbeck process with correlation time $τ_c$ driven by the Cauchy white noise. We discuss the asymptotic diffusive properties of this model and obtain the same analytical expression for the diffusion coefficient as in the first case which is thus independent on $τ_c$. Remarkably, for $τ_c>τ_D$ the crossover from a non-Gaussian to a Gaussian distribution of displacements takes place at a time $τ_G$ which can be considerably larger than the persistence time $τ_D$.

preprint2016arXiv

Effective medium approximation for lattice random walks with long-range jumps

We consider the random walk on a lattice with random transition rates and arbitrarily long-range jumps. We employ Bruggeman's effective medium approximation (EMA) to find the disorder averaged (coarse-grained) dynamics. The EMA procedure replaces the disordered system with a cleverly guessed reference system in a self-consistent manner. We give necessary conditions on the reference system and discuss possible physical mechanisms of anomalous diffusion. In case of a power-law scaling between transition rates and distance, lattice variants of Levy-ights emerge as the effective medium, and the problem is solved analytically, bearing the effective anomalous diffusivity. Finally, we discuss several example distributions, and demonstrate very good agreement with numerical simulations.

preprint2016arXiv

Reconstructing Interaction Potentials in Thin Films from Real-Space Images

We demonstrate that an inverse Monte Carlo approach allows to reconstruct effective interaction potentials from real-space images. The method is exemplified on monomolecular ethanol-water films imaged with scanning force microscopy (SFM), which provides the spatial distribution of the molecules. Direct Monte Carlo simulations with the reconstructed potential allow for obtaining characteristics of the system which are unavailable in the experiment, such as the heat capacity of the monomolecularly thin film, and for a prediction of the critical temperature of the demixing transition.

preprint2015arXiv

Quantifying the non-ergodicity of scaled Brownian motion

We examine the non-ergodic properties of scaled Brownian motion, a non-stationary stochastic process with a time dependent diffusivity of the form $D(t)\simeq t^{α-1}$. We compute the ergodicity breaking parameter EB in the entire range of scaling exponents $α$, both analytically and via extensive computer simulations of the stochastic Langevin equation. We demonstrate that in the limit of long trajectory lengths $T$ and short lag times $Δ$ the EB parameter as function of the scaling exponent $α$ has no divergence at $α=1/2$ and present the asymptotes for EB in different limits. We generalise the analytical and simulations results for the time averaged and ergodic properties of scaled Brownian motion in the presence of ageing, that is, when the observation of the system starts only a finite time span after its initiation. The approach developed here for the calculation of the higher time averaged moments of the particle displacement can be applied to derive the ergodic properties of other stochastic processes such as fractional Brownian motion.

preprint2015arXiv

Robust linear regression with broad distributions of errors

We consider the problem of linear fitting of noisy data in the case of broad (say $α$-stable) distributions of random impacts ("noise"), which can lack even the first moment. This situation, common in statistical physics of small systems, in Earth sciences, in network science or in econophysics, does not allow for application of conventional Gaussian maximum-likelihood estimators resulting in usual least-squares fits. Such fits lead to large deviations of fitted parameters from their true values due to the presence of outliers. The approaches discussed here aim onto the minimization of the width of the distribution of residua. The corresponding width of the distribution can either be defined via the interquantile distance of the corresponding distributions or via the scale parameter in its characteristic function. The methods provide the robust regression even in the case of short samples with large outliers, and are equivalent to the normal least squares fit for the Gaussian noises. Our discussion is illustrated by numerical examples.

preprint2014arXiv

Absence of Anderson localization of light in a random ensemble of point scatterers

As discovered by Philip Anderson in 1958, strong disorder can block propagation of waves and lead to the localization of wave-like excitations in space. Anderson localization of light is particularly exciting in view of its possible applications for random lasing or quantum information processing. We show that, surprisingly, Anderson localization of light cannot be achieved in a random three-dimensional ensemble of point scattering centers that is the simplest and widespread model to study the multiple scattering of waves. Localization is recovered if the vector character of light is neglected. This shows that, at least for point scatterers, the polarization of light plays an important role in the Anderson localization problem.

preprint2014arXiv

Anomalous Diffusion of Self-Propelled Particles in Directed Random Environments

We theoretically study the transport properties of self-propelled particles on complex structures, such as motor proteins on filament networks. A general master equation formalism is developed to investigate the persistent motion of individual random walkers, which enables us to identify the contributions of key parameters: the motor processivity, and the anisotropy and heterogeneity of the underlying network. We prove the existence of different dynamical regimes of anomalous motion, and that the crossover times between these regimes as well as the asymptotic diffusion coefficient can be increased by several orders of magnitude within biologically relevant control parameter ranges. In terms of motion in continuous space, the interplay between stepping strategy and persistency of the walker is established as a source of anomalous diffusion at short and intermediate time scales.

preprint2014arXiv

Codifference as a practical tool to measure interdependence

Correlation and spectral analysis represent the standard tools to study interdependence in statistical data. However, for the stochastic processes with heavy-tailed distributions such that the variance diverges, these tools are inadequate. The heavy-tailed processes are ubiquitous in nature and finance. We here discuss codifference as a convenient measure to study statistical interdependence, and we aim to give a short introductory review of its properties. By taking different known stochastic processes as generic examples, we present explicit formulas for their codifferences. We show that for the Gaussian processes codifference is equivalent to covariance. For processes with finite variance these two measures behave similarly with time. For the processes with infinite variance the covariance does not exist, however, the codifference is relevant. We demonstrate the practical importance of the codifference by extracting this function from simulated as well as from real experimental data. We conclude that the codifference serves as a convenient practical tool to study interdependence for stochastic processes with both infinite and finite variances as well.

preprint2014arXiv

Dynamical localization and eigenstate localization in trap models

The one-dimensional random trap model with a power-law distribution of mean sojourn times exhibits a phenomenon of dynamical localization in the case where diffusion is anomalous: The probability to find two independent walkers at the same site, as given by the participation ratio, stays constant and high in a broad domain of intermediate times. This phenomenon is absent in dimensions two and higher. In finite lattices of all dimensions the participation ratio finally equilibrates to a different final value. We numerically investigate two-particle properties in a random trap model in one and in three dimensions, using a method based on spectral decomposition of the transition rate matrix. The method delivers a very effective computational scheme producing numerically exact results for the averages over thermal histories and initial conditions in a given landscape realization. Only a single averaging procedure over disorder realizations is necessary. The behavior of the participation ratio is compared to other measures of localization, as for example to the states' gyration radius, according to which the dynamically localized states are extended. This means that although the particles are found at the same site with a high probability, the typical distance between them grows. Moreover the final equilibrium state is extended both with respect to its gyration radius and to its Lyapunov exponent. In addition, we show that the phenomenon of dynamical localization is only marginally connected with the spectrum of the transition rate matrix, and is dominated by the properties of its eigenfunctions which differ significantly in dimensions one and three.

preprint2014arXiv

Model of lateral diffusion in ultrathin layered films

We consider the diffusion of markers in a layered medium, with the lateral diffusion coefficient being the function of hight. We show that the probability density of the lateral displacements follows one-dimensional Batchelor's equation with time-dependent diffusion coefficient governed by the particles' redistribution in height. For the film of a finite thickness the resulting mean squared displacement exhibits superdiffusion at short times and crosses over to normal diffusion at long times. The approach is used for description of experimental results on inhomogeneous molecular diffusion in thin liquid films deposited on solid surfaces.

preprint2014arXiv

Spectral properties of fractional Fokker-Plank operator for the Lévy flight in a harmonic potential

We present a detailed analysis of the eigenfunctions of the Fokker-Planck operator for the Lévy-Ornstein-Uhlenbeck process, their asymptotic behavior and recurrence relations, explicit expressions in coordinate space for the special cases of the Ornstein-Uhlenbeck process with Gaussian and with Cauchy white noise and for the transformation kernel, which maps the fractional Fokker-Planck operator of the Cauchy-Ornstein-Uhlenbeck process to the non-fractional Fokker-Planck operator of the usual Gaussian Ornstein-Uhlenbeck process. We also describe how non-spectral relaxation can be observed in bounded random variables of the Lévy-Ornstein-Uhlenbeck process and their correlation functions.

preprint2014arXiv

Tracer diffusion inside fibrinogen layers

We investigate the obstructed motion of tracer (test) particles in crowded environments by carrying simulations of two-dimensional Gaussian random walk in model fibrinogen monolayers of different orientational ordering. The fibrinogen molecules are significantly anisotropic and therefore they can form structures where orientational ordering, similar to the one observed in nematic liquid crystals, appears. The work focuses on the dependence between level of the orientational order (degree of environmental crowding) of fibrinogen molecules inside a layer and non-Fickian character of the diffusion process of spherical tracer particles moving within the domain. It is shown that in general particles motion is subdiffusive and strongly anisotropic, and its characteristic features significantly change with the orientational order parameter, concentration of fibrinogens and radius of a diffusing probe.

preprint2013arXiv

Disentangling sources of anomalous diffusion

We show that some important properties of subbdiffusion of unknown origin (including those of mixed origin) can be easily assessed when findeng the "fundamental moment" of the corresponding process, i.e., the one which is additive in time. In subordinated processes, the index of the fundamental moment is inherited from the parent process and its time-dependence from the leading one. In models of particle's motion in disordered potentials, the index is governed by the structural part of the disorder while the time dependence is given by its energetic part.

preprint2013arXiv

Scaled Brownian motion as a mean field model for continuous time random walks

We consider scaled Brownian motion (sBm), a random process described by a diffusion equation with explicitly time-dependent diffusion coefficient $D(t) = D_0 t^{α- 1}$ (Batchelor's equation) which, for $α< 1$, is often used for fitting experimental data for subdiffusion of unclear genesis. We show that this process is a close relative of subdiffusive continuous-time random walks and describes the motion of the center of mass of a cloud of independent walkers. It shares with subdiffusive CTRW its non-stationary and non-ergodic properties. The non-ergodicity of sBm does not however go hand in hand with strong difference between its different realizations: its heterogeneity ("ergodicity breaking") parameter tends to zero for long trajectories.

preprint2012arXiv

A Test for Determining a Subdiffusive Model in Ergodic Systems from Single Trajectories

Experiments on particles' motion in living cells show that it is often subdiffusive. This subdiffusion may be due to trapping, percolation-like structures, or viscoelatic behavior of the medium. While the models based on trapping (leading to continuous-time random walks) can easily be distinguished from the rest by testing their non-ergodicity, the latter two cases are harder to distinguish. We propose a statistical test for distinguishing between these two based on the space-filling properties of trajectories, and prove its feasibility and specificity using synthetic data. We moreover present a flow-chart for making a decision on a type of subdiffusion for a broader class of models.

preprint2012arXiv

Active Particles Forced by an Asymmetric Dichotomous Angle Drive

We analyze the dynamics of particles in two dimensions with constant speed and a stochastic switching angle dynamics defined by a correlated dichotomous Markov process (telegraph noise) plus Gaussian white noise. We study various cases of the asymptotic diffusional motion of the particle which is characterized by the effective diffusion coefficient. Expressions for this coefficient are derived and discussed in dependence on the correlation time and the intensity of the noise. The situation with a given mean curvature is of special interest since a non-monotonic behavior of the effective diffusion coefficient as function of the noise intensity and correlation time is found. A timescale matching condition for maximal diffusion is formulated.

preprint2012arXiv

Bulk-mediated diffusion on a planar surface: full solution

We consider the effective surface motion of a particle that intermittently unbinds from a planar surface and performs bulk excursions. Based on a random walk approach we derive the diffusion equations for surface and bulk diffusion including the surface-bulk coupling. From these exact dynamic equations we analytically obtain the propagator of the effective surface motion. This approach allows us to deduce a superdiffusive, Cauchy-type behavior on the surface, together with exact cutoffs limiting the Cauchy form. Moreover we study the long-time dynamics for the surface motion.

preprint2012arXiv

Diffusion of small particles in a solid polymeric medium

We analyze diffusion of small particles in a solid polymeric medium taking into account a short range particle-polymer interaction. The system is modeled by a particle diffusion on a ternary lattice where the sites occupied by polymer segments are blocked, the ones forming the hull of the chains correspond to the places at which the interaction takes place, and the rest are voids, in which the diffusion is free. In the absence of interaction the diffusion coefficient shows only a weak dependence on the polymer chain length and its behavior strongly resembles usual site percolation. In presence of interactions the diffusion coefficient (and especially its temperature dependence) shows a non-trivial behavior depending on the sign of interaction and on whether the voids and the hulls of the chains percolate or not. The temperature dependence may be Arrhenius-like or strongly non-Arrhenius, depending on parameters. The analytical results obtained within the effective medium approximation are in qualitative agreement with those of Monte Carlo simulations.

preprint2012arXiv

Impact of inertia on biased Brownian transport in confined geometries

We consider the impact of inertia on biased Brownian motion of point particles in a two-dimensional channel with sinusoidally varying width. If the time scales of the problem separate, the adiabatic elimination of the transverse degrees of freedom leads to an effective description for the motion along the channel given by the potential of mean force. The possibility of such description is intimately connected with equipartition. Numerical simulations show that in the presence of external bias the equipartition may break down leading to non-monotonic dependence of mobility on external force and several other interesting effects.

preprint2012arXiv

Normal and anomalous diffusion in random potential landscapes

A relation between the effective diffusion coefficient in a lattice with random site energies and random trasition rates and the macroscopic conductivity in a random resistor network allows for elucidating possible sources of anomalous diffusion in random potential models. We show that subdiffusion is only possible either if the mean Boltzmann factor in the corresponding potential diverges or if the percolation concentration in the system is equal to unity (or both), and that superdiffusion is impossible in our system under any condition. We show also other useful applications of this relation.

preprint2011arXiv

Effective surface motion on a reactive cylinder of particles that perform intermittent bulk diffusion

In many biological and small scale technological applications particles may transiently bind to a cylindrical surface. In between two binding events the particles diffuse in the bulk, thus producing an effective translation on the cylinder surface. We here derive the effective motion on the surface, allowing for additional diffusion on the cylinder surface itself. We find explicit solutions for the number of adsorbed particles at one given instant, the effective surface displacement, as well as the surface propagator. In particular sub- and superdiffusive regimes are found, as well as an effective stalling of diffusion visible as a plateau in the mean squared displacement. We also investigate the corresponding first passage and first return problems.

preprint2011arXiv

Harmonic oscillator under Levy noise: Unexpected properties in the phase space

A harmonic oscillator under influence of the noise is a basic model of various physical phenomena. Under Gaussian white noise the position and velocity of the oscillator are independent random variables which are distributed according to the bivariate Gaussian distribution with elliptic level lines. The distribution of phase is homogeneous. None of these properties hold in the general Lévy case. Thus, the level lines of the joint probability density are not elliptic. The coordinate and the velocity of the oscillator are strongly dependent, and this dependence is quantified by introducing the corresponding parameter ("width deficit"). The distribution of the phase is inhomogeneous and highly nontrivial.

preprint2011arXiv

The Unequal Twins - Probability Distributions Aren't Everything

It is the common lore to assume that knowing the equation for the probability distribution function (PDF) of a stochastic model as a function of time tells the whole picture defining all other characteristics of the model. We show that this is not the case by comparing two exactly solvable models of anomalous diffusion due to geometric constraints: The comb model and the random walk on a random walk (RWRW). We show that though the two models have exactly the same PDFs, they differ in other respects, like their first passage time (FPT) distributions, their autocorrelation functions and their aging properties.

preprint2010arXiv

Kramers escape driven by fractional Brownian motion

We investigate the Kramers escape from a potential well of a test particle driven by fractional Gaussian noise with Hurst exponent 0<H<1. From a numerical analysis we demonstrate the exponential distribution of escape times from the well and analyze in detail the dependence of the mean escape time as function of H and the particle diffusivity D. We observe different behavior for the subdiffusive (antipersistent) and superdiffusive (persistent) domains. In particular we find that the escape becomes increasingly faster for decreasing values of H, consistent with previous findings on the first passage behavior. Approximate analytical calculations are shown to support the numerically observed dependencies.

preprint2010arXiv

Lévy Ratchet in a Weak Noise Limit: Theory and Simulation

We study the motion of a particle embedded in a time independent periodic potential with broken mirror symmetry and subjected to a Lévy noise possessing Lévy stable probability law (Lévy ratchet). We develop analytical approach to the problem based on the asymptotic probabilistic method of decomposition proposed by P. Imkeller and I. Pavlyukevich [J. Phys. A {\bf39}, L237 (2006); Stoch. Proc. Appl. {\bf116}, 611 (2006)]. We derive analytical expressions for the quantities characterizing the particle motion, namely the splitting probabilities of first escape from a single well, the transition probabilities and the particle current. A particular attention is devoted to the interplay between the asymmetry of the ratchet potential and the asymmetry (skewness) of the Lévy noise. Intensive numerical simulations demonstrate a good agreement with the analytical predictions for sufficiently small intensities of the Lévy noise driving the particle.

preprint2010arXiv

On the relation between event-based and time-based current statistics

Current statistics can be calculated in various ways. Event-based approaches use the statistics of the number of events occuring during a given time. Time-based approaches use the statistics of the time needed to reach a given number of events. By analyzing normal as well as anomalous statistics of nonequilibrium currents through a two level system in contact with two different reservoirs, we investigate the conditions under which these different statistics do or do not yield identical predictions. We rely on the continuous time random walk formulation introduced in our earlier work [Phys. Rev. E 77, 051119 (2008)].

preprint2010arXiv

Quantum memory for light via stimulated off-resonant Raman process: beyond the three-level Lambda-scheme approximation

We consider a quantum memory scheme based on the conversion of a signal pulse into a long-lived spin coherence via stimulated off-resonant Raman process. For a storing medium consisting of alkali atoms, we have calculated the Autler-Townes resonance structure created by a strong control field. By taking into account the upper hyperfine states of the D1 optical transition, we show important deviations from the predictions of the usual three-level Lambda-scheme approximation and we demonstrate an enhancement of the process for particular detunings of the control. We estimate the memory efficiency one can obtain using this configuration.

preprint2010arXiv

Stationary states in single-well potentials under symmetric Levy noises

We discuss the existence of stationary states for subharmonic potentials $V(x) \propto |x|^c$, $c<2$, under action of symmetric $α$-stable noises. We show analytically that the necessary condition for the existence of the steady state is $c>2-α$. These states are characterized by heavy-tailed probability density functions which decay as $P(x) \propto x^{-(c+α-1)}$ for $|x| \to \infty$, i.e. stationary states posses a heavier tail than the corresponding $α$-stable law. Monte Carlo simulations confirm the existence of such stationary states and the form of the tails of corresponding probability densities.

preprint2008arXiv

Bulk-mediated surface diffusion on a cylinder: propagators and crossovers

We consider the effective surface motion of a particle that freely diffuses in the bulk and intermittently binds to that surface. From an exact approach we derive various regimes of the effective surface motion characterized by physical rates for binding/unbinding and the bulk diffusivity. We obtain a transient regime of superdiffusion and, in particular, a saturation regime characteristic for the cylindrical geometry. This saturation, however, in a finite system is not terminal but eventually turns over to normal surface diffusion. The first passage behavior of particles to the cylinder surface is derived. Consequences for actual systems are discussed.

preprint2008arXiv

Front propagation in A+B -> 2A reaction under subdiffusion

We consider an irreversible autocatalytic conversion reaction A+B -> 2A under subdiffusion described by continuous time random walks. The reactants' transformations take place independently on their motion and are described by constant rates. The analog of this reaction in the case of normal diffusion is described by the Fisher-Kolmogorov-Petrovskii-Piskunov (FKPP) equation leading to the existence of a nonzero minimal front propagation velocity which is really attained by the front in its stable motion. We show that for subdiffusion this minimal propagation velocity is zero, which suggests propagation failure.

preprint2001arXiv

Fractional Klein-Kramers equation for superdiffusive transport: normal versus anomalous time evolution in a differential L{é}vy walk model

We introduce a fractional Klein-Kramers equation which describes sub-ballistic superdiffusion in phase space in the presence of a space-dependent external force field. This equation defines the differential L{é}vy walk model whose solution is shown to be non-negative. In the velocity coordinate, the probability density relaxes in Mittag-Leffler fashion towards the Maxwell distribution whereas in the space coordinate, no stationary solution exists and the temporal evolution of moments exhibits a competition between Brownian and anomalous contributions.