Fractional Telegrapher's Equation under Resetting: Non-Equilibrium Stationary States and First-Passage Times.

We consider two different time fractional telegrapher's equations under stochastic resetting. Using the integral decomposition method, we found the probability density functions and the mean squared displacements. In the long-time limit, the system approaches non-equilibrium stationary states, while the mean squared displacement saturates due to the resetting mechanism. We also obtain the fractional telegraph process as a subordinated telegraph process by introducing operational time such that the physical time is considered as a Lévy stable process whose characteristic function is the Lévy stable distribution. We also analyzed the survival probability for the first-passage time problem and found the optimal resetting rate for which the corresponding mean first-passage time is minimal.

Random Walks on Comb-like Structures under Stochastic Resetting.

We study the long-time dynamics of the mean squared displacement of a random walker moving on a comb structure under the effect of stochastic resetting. We consider that the walker's motion along the backbone is diffusive and it performs short jumps separated by random resting periods along fingers. We take into account two different types of resetting acting separately: global resetting from any point in the comb to the initial position and resetting from a finger to the corresponding backbone. We analyze the interplay between the waiting process and Markovian and non-Markovian resetting processes on the overall mean squared displacement. The Markovian resetting from the fingers is found to induce normal diffusion, thereby minimizing the trapping effect of fingers. In contrast, for non-Markovian local resetting, an interesting crossover with three different regimes emerges, with two of them subdiffusive and one of them diffusive. Thus, an interesting interplay between the exponents characterizing the waiting time distributions of the subdiffusive random walk and resetting takes place. As for global resetting, its effect is even more drastic as it precludes normal diffusion. Specifically, such a resetting can induce a constant asymptotic mean squared displacement in the Markovian case or two distinct regimes of subdiffusive motion in the non-Markovian case.

Random Walks on Networks with Centrality-Based Stochastic Resetting.

We introduce a refined way to diffusely explore complex networks with stochastic resetting where the resetting site is derived from node centrality measures. This approach differs from previous ones, since it not only allows the random walker with a certain probability to jump from the current node to a deliberately chosen resetting node, rather it enables the walker to jump to the node that can reach all other nodes faster. Following this strategy, we consider the resetting site to be the geometric center, the node that minimizes the average travel time to all the other nodes. Using the established Markov chain theory, we calculate the Global Mean First Passage Time (GMFPT) to determine the search performance of the random walk with resetting for different resetting node candidates individually. Furthermore, we compare which nodes are better resetting node sites by comparing the GMFPT for each node. We study this approach for different topologies of generic and real-life networks. We show that, for directed networks extracted for real-life relationships, this centrality focused resetting can improve the search to a greater extent than for the generated undirected networks. This resetting to the center advocated here can minimize the average travel time to all other nodes in real networks as well. We also present a relationship between the longest shortest path (the diameter), the average node degree and the GMFPT when the starting node is the center. We show that, for undirected scale-free networks, stochastic resetting is effective only for networks that are extremely sparse with tree-like structures as they have larger diameters and smaller average node degrees. For directed networks, the resetting is beneficial even for networks that have loops. The numerical results are confirmed by analytic solutions. Our study demonstrates that the proposed random walk approach with resetting based on centrality measures reduces the memoryless search time for targets in the examined network topologies.

Income inequality and mobility in geometric Brownian motion with stochastic resetting: theoretical results and empirical evidence of non-ergodicity.

We explore the role of non-ergodicity in the relationship between income inequality, the extent of concentration in the income distribution, and income mobility, the feasibility of an individual to change their position in the income rankings. For this purpose, we use the properties of an established model for income growth that includes 'resetting' as a stabilizing force to ensure stationary dynamics. We find that the dynamics of inequality is regime-dependent: it may range from a strictly non-ergodic state where this phenomenon has an increasing trend, up to a stable regime where inequality is steady and the system efficiently mimics ergodicity. Mobility measures, conversely, are always stable over time, but suggest that economies become less mobile in non-ergodic regimes. By fitting the model to empirical data for the income share of the top earners in the USA, we provide evidence that the income dynamics in this country is consistently in a regime in which non-ergodicity characterizes inequality and immobility. Our results can serve as a simple rationale for the observed real-world income dynamics and as such aid in addressing non-ergodicity in various empirical settings across the globe. This article is part of the theme issue 'Kinetic exchange models of societies and economies'.

Non-Markovian SIR epidemic spreading model of COVID-19.

We introduce non-Markovian SIR epidemic spreading model inspired by the characteristics of the COVID-19, by considering discrete- and continuous-time versions. The distributions of infection intensity and recovery period may take an arbitrary form. By taking corresponding choice of these functions, it is shown that the model reduces to the classical Markovian case. The epidemic threshold is analytically determined for arbitrary functions of infectivity and recovery and verified numerically. The relevance of the model is shown by modeling the first wave of the epidemic in Italy, Spain and the UK, in the spring, 2020.

Generalised Geometric Brownian Motion: Theory and Applications to Option Pricing.

Classical option pricing schemes assume that the value of a financial asset follows a geometric Brownian motion (GBM). However, a growing body of studies suggest that a simple GBM trajectory is not an adequate representation for asset dynamics, due to irregularities found when comparing its properties with empirical distributions. As a solution, we investigate a generalisation of GBM where the introduction of a memory kernel critically determines the behaviour of the stochastic process. We find the general expressions for the moments, log-moments, and the expectation of the periodic log returns, and then obtain the corresponding probability density functions using the subordination approach. Particularly, we consider subdiffusive GBM (sGBM), tempered sGBM, a mix of GBM and sGBM, and a mix of sGBMs. We utilise the resulting generalised GBM (gGBM) in order to examine the empirical performance of a selected group of kernels in the pricing of European call options. Our results indicate that the performance of a kernel ultimately depends on the maturity of the option and its moneyness.

Fractional heterogeneous telegraph processes: Interplay between heterogeneity, memory, and stochastic resetting.

Fractional heterogeneous telegraph processes are considered in the framework of telegrapher's equations accompanied by memory effects. The integral decomposition method is developed for the rigorous treating of the problem. Exact solutions for the probability density functions and the mean squared displacements are obtained. A relation between the fractional heterogeneous telegrapher's equation and the corresponding Langevin equation has been established in the framework of the developed subordination approach. The telegraph process in the presence of stochastic resetting has been studied, as well. An exact expression for both the nonequilibrium stationary distributions/states and the mean squared displacements are obtained.

Ornstein-Uhlenbeck process and generalizations: Particle dynamics under comb constraints and stochastic resetting.

The Ornstein-Uhlenbeck process is interpreted as Brownian motion in a harmonic potential. This Gaussian Markov process has a bounded variance and admits a stationary probability distribution, in contrast to the standard Brownian motion. It also tends to a drift towards its mean function, and such a process is called mean reverting. Two examples of the generalized Ornstein-Uhlenbeck process are considered. In the first one, we study the Ornstein-Uhlenbeck process on a comb model, as an example of the harmonically bounded random motion in the topologically constrained geometry. The main dynamical characteristics (as the first and the second moments) and the probability density function are studied in the framework of both the Langevin stochastic equation and the Fokker-Planck equation. The second example is devoted to the study of the effects of stochastic resetting on the Ornstein-Uhlenbeck process, including stochastic resetting in the comb geometry. Here the nonequilibrium stationary state is the main question in task, where the two divergent forces, namely, the resetting and the drift towards the mean, lead to compelling results in the cases of both the Ornstein-Uhlenbeck process with resetting and its generalization on the two-dimensional comb structure.

Geometric Brownian motion under stochastic resetting: A stationary yet nonergodic process.

We study the effects of stochastic resetting on geometric Brownian motion with drift (GBM), a canonical stochastic multiplicative process for nonstationary and nonergodic dynamics. Resetting is a sudden interruption of a process, which consecutively renews its dynamics. We show that, although resetting renders GBM stationary, the resulting process remains nonergodic. Quite surprisingly, the effect of resetting is pivotal in manifesting the nonergodic behavior. In particular, we observe three different long-time regimes: a quenched state, an unstable state, and a stable annealed state depending on the resetting strength. Notably, in the last regime, the system is self-averaging and thus the sample average will always mimic ergodic behavior establishing a stand-alone feature for GBM under resetting. Crucially, the above-mentioned regimes are well separated by a self-averaging time period which can be minimized by an optimal resetting rate. Our results can be useful to interpret data emanating from stock market collapse or reconstitution of investment portfolios.

Continuous time random walks under Markovian resetting.

We investigate the effects of Markovian resetting events on continuous time random walks where the waiting times and the jump lengths are random variables distributed according to power-law probability density functions. We prove the existence of a nonequilibrium stationary state and finite mean first arrival time. However, the existence of an optimum reset rate is conditioned to a specific relationship between the exponents of both power-law tails. We also investigate the search efficiency by finding the optimal random walk which minimizes the mean first arrival time in terms of the reset rate, the distance of the initial position to the target, and the characteristic transport exponents.

Hitting times in turbulent diffusion due to multiplicative noise.

We study a distribution of times of the first arrivals to absorbing targets in turbulent diffusion, which is due to a multiplicative noise. Two examples of dynamical systems with a multiplicative noise are studied. The first one is a random process according to inhomogeneous diffusion, which is also known as a geometric Brownian motion in the Black-Scholes model. The second model is due to a random processes on a two-dimensional comb, where inhomogeneous advection is possible only along the backbone, while Brownian diffusion takes place inside the branches. It is shown that in both cases turbulent diffusion takes place as the one-dimensional random process with the log-normal distribution in the presence of absorbing targets, which are characterized by the Lévy-Smirnov distribution for the first hitting times.

Reaction and ultraslow diffusion on comb structures.

A two-dimensional (2D) comb model is proposed to characterize reaction-ultraslow diffusion of tracers both in backbones (x direction) and side branches (y direction) of the comblike structure with two memory kernels. The memory kernels include Dirac delta, power-law, and logarithmic and inverse Mittag-Leffler (ML) functions, which can also be considered as the structural functions in the time structural derivative. Based on the comb model, ultraslow diffusion on a fractal comb structure is also investigated by considering spatial fractal geometry of the backbone volume. The mean squared displacement (MSD) and the corresponding concentration of the tracers, i.e., the solution of the comb model, are derived for reactive and conservative tracers. For a fractal structure of backbones, the derived MSDs and corresponding solutions depend on the backbone's fractal dimension. The proposed 2D comb model with different kernel functions is feasible to describe ultraslow diffusion in both the backbone and side branches of the comblike structure.

Anomalous diffusion on a fractal mesh.

We present exact analytical results for properties of anomalous diffusion on a fractal mesh. The fractal mesh structure is a direct product of two fractal sets, one belonging to a main branch of backbones, the other to the side branches of fingers. Both fractal sets are constructed on the entire (infinite) y and x axes. We suggest a special algorithm in order to construct such sets out of standard Cantor sets embedded in the unit interval. The transport properties of the fractal mesh are studied, in particular, subdiffusion along the backbones is obtained with the dispersion relation 〈x^{2}(t)〉≃t^{ß}, where the transport exponent ß<1 is determined by the fractal dimensions of both backbone and fingers. Superdiffusion with ß>1 has been observed as well when the environment is controlled by means of a memory kernel.

Fractional diffusion on a fractal grid comb.

A grid comb model is a generalization of the well known comb model, and it consists of N backbones. For N=1 the system reduces to the comb model where subdiffusion takes place with the transport exponent 1/2. We present an exact analytical evaluation of the transport exponent of anomalous diffusion for finite and infinite number of backbones. We show that for an arbitrarily large but finite number of backbones the transport exponent does not change. Contrary to that, for an infinite number of backbones, the transport exponent depends on the fractal dimension of the backbone structure.

Distributed-order diffusion equations and multifractality: Models and solutions.

We study distributed-order time fractional diffusion equations characterized by multifractal memory kernels, in contrast to the simple power-law kernel of common time fractional diffusion equations. Based on the physical approach to anomalous diffusion provided by the seminal Scher-Montroll-Weiss continuous time random walk, we analyze both natural and modified-form distributed-order time fractional diffusion equations and compare the two approaches. The mean squared displacement is obtained and its limiting behavior analyzed. We derive the connection between the Wiener process, described by the conventional Langevin equation and the dynamics encoded by the distributed-order time fractional diffusion equation in terms of a generalized subordination of time. A detailed analysis of the multifractal properties of distributed-order diffusion equations is provided.