TY - JOUR A1 - Petreska, Irina A1 - Pejov, Ljupco A1 - Sandev, Trifce A1 - Kocarev, Ljupčo A1 - Metzler, Ralf T1 - Tuning of the dielectric relaxation and complex susceptibility in a system of polar molecules: a generalised model based on rotational diffusion with resetting JF - Fractal and fractional N2 - The application of the fractional calculus in the mathematical modelling of relaxation processes in complex heterogeneous media has attracted a considerable amount of interest lately. The reason for this is the successful implementation of fractional stochastic and kinetic equations in the studies of non-Debye relaxation. In this work, we consider the rotational diffusion equation with a generalised memory kernel in the context of dielectric relaxation processes in a medium composed of polar molecules. We give an overview of existing models on non-exponential relaxation and introduce an exponential resetting dynamic in the corresponding process. The autocorrelation function and complex susceptibility are analysed in detail. We show that stochastic resetting leads to a saturation of the autocorrelation function to a constant value, in contrast to the case without resetting, for which it decays to zero. The behaviour of the autocorrelation function, as well as the complex susceptibility in the presence of resetting, confirms that the dielectric relaxation dynamics can be tuned by an appropriate choice of the resetting rate. The presented results are general and flexible, and they will be of interest for the theoretical description of non-trivial relaxation dynamics in heterogeneous systems composed of polar molecules. KW - rotational diffusion KW - memory kernel KW - Fokker-Planck equation KW - non-exponential relaxation KW - autocorrelation function KW - complex KW - susceptibility Y1 - 2022 U6 - https://doi.org/10.3390/fractalfract6020088 SN - 2504-3110 VL - 6 IS - 2 PB - MDPI AG, Fractal Fract Editorial Office CY - Basel ER - TY - JOUR A1 - Stojkoski, Viktor A1 - Sandev, Trifce A1 - Basnarkov, Lasko A1 - Kocarev, Ljupco A1 - Metzler, Ralf T1 - Generalised geometric Brownian motion BT - theory and applications to option pricing JF - Entropy N2 - 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. KW - geometric Brownian motion KW - Fokker– Planck equation KW - Black– Scholes model KW - option pricing Y1 - 2020 U6 - https://doi.org/10.3390/e22121432 SN - 1099-4300 VL - 22 IS - 12 PB - MDPI CY - Basel ER - TY - JOUR A1 - Sandev, Trifce A1 - Iomin, Alexander A1 - Kocarev, Ljupco T1 - Hitting times in turbulent diffusion due to multiplicative noise JF - Physical review : E, Statistical, nonlinear and soft matter physics N2 - 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 Levy-Smirnov distribution for the first hitting times. Y1 - 2020 U6 - https://doi.org/10.1103/PhysRevE.102.042109 SN - 2470-0045 SN - 2470-0053 VL - 102 IS - 4 PB - American Institute of Physics CY - Woodbury, NY ER - TY - JOUR A1 - Basnarkov, Lasko A1 - Tomovski, Igor A1 - Sandev, Trifce A1 - Kocarev, Ljupčo T1 - Non-Markovian SIR epidemic spreading model of COVID-19 JF - Chaos, solitons & fractals : applications in science and engineering ; an interdisciplinary journal of nonlinear science N2 - 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. KW - Epidemic spreading models KW - Non-Markovian processes KW - COVID-19 KW - SIR model Y1 - 2022 U6 - https://doi.org/10.1016/j.chaos.2022.112286 SN - 0960-0779 SN - 1873-2887 VL - 160 PB - Elsevier CY - Oxford [u.a.] ER - TY - JOUR A1 - Sandev, Trifce A1 - Domazetoski, Viktor A1 - Kocarev, Ljupco A1 - Metzler, Ralf A1 - Chechkin, Aleksei T1 - Heterogeneous diffusion with stochastic resetting JF - Journal of physics : A, Mathematical and theoretical N2 - We study a heterogeneous diffusion process (HDP) with position-dependent diffusion coefficient and Poissonian stochastic resetting. We find exact results for the mean squared displacement and the probability density function. The nonequilibrium steady state reached in the long time limit is studied. We also analyse the transition to the non-equilibrium steady state by finding the large deviation function. We found that similarly to the case of the normal diffusion process where the diffusion length grows like t (1/2) while the length scale xi(t) of the inner core region of the nonequilibrium steady state grows linearly with time t, in the HDP with diffusion length increasing like t ( p/2) the length scale xi(t) grows like t ( p ). The obtained results are verified by numerical solutions of the corresponding Langevin equation. KW - heterogeneous diffusion KW - Fokker-Planck equation KW - Langevin equation KW - stochastic resetting KW - nonequilibrium stationary state KW - large deviation function Y1 - 2022 U6 - https://doi.org/10.1088/1751-8121/ac491c SN - 1751-8113 SN - 1751-8121 VL - 55 IS - 7 PB - IOP Publ. Ltd. CY - Bristol ER - TY - JOUR A1 - Peng, Junhao A1 - Sandev, Trifce A1 - Kocarev, Ljupco T1 - First encounters on Bethe lattices and Cayley trees JF - Communications in nonlinear science & numerical simulation N2 - In this work we consider the first encounter problems between a fixed and/or mobile target A and a moving trap B on Bethe lattices and Cayley trees. The survival probabilities (SPs) of the target A on the both kinds of structures are considered analytically and compared. On Bethe lattices, the results show that the fixed target will still prolong its survival time, whereas, on Cayley trees, there are some initial positions where the target should move to prolong its survival time. The mean first encounter time (MFET) for mobile target A is evaluated numerically and compared with the mean first passage time (MFPT) for the fixed target A. Different initial settings are addressed and clear boundaries are obtained. These findings are helpful for optimizing the strategy to prolong the survival time of the target or to speed up the search process on Cayley trees, in relation to the target's movement and the initial position configuration of the two walkers. We also present a new method, which uses a small amount of memory, for simulating random walks on Cayley trees. (C) 2020 Elsevier B.V. All rights reserved. KW - Random walks KW - Survival probability KW - Mean first encounter time KW - Bethe KW - lattices KW - Cayley trees Y1 - 2021 U6 - https://doi.org/10.1016/j.cnsns.2020.105594 SN - 1007-5704 SN - 1878-7274 VL - 95 PB - Elsevier CY - Amsterdam ER - TY - JOUR A1 - Stojkoski, Viktor A1 - Jolakoski, Petar A1 - Pal, Arnab A1 - Sandev, Trifce A1 - Kocarev, Ljupco A1 - Metzler, Ralf T1 - Income inequality and mobility in geometric Brownian motion with stochastic resetting: theoretical results and empirical evidence of non-ergodicity JF - Philosophical transactions of the Royal Society A: Mathematical, physical and engineering sciences N2 - 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'. KW - income inequality KW - income mobility KW - geometric Brownian motion KW - non-ergodicity KW - stochastic resetting Y1 - 2022 U6 - https://doi.org/10.1098/rsta.2021.0157 SN - 1364-503X SN - 1471-2962 VL - 380 IS - 2224 PB - Royal Society CY - London ER - TY - JOUR A1 - Stojkoski, Viktor A1 - Sandev, Trifce A1 - Kocarev, Ljupco A1 - Pal, Arnab T1 - Autocorrelation functions and ergodicity in diffusion with stochastic resetting JF - Journal of physics : A, Mathematical and theoretical N2 - Diffusion with stochastic resetting is a paradigm of resetting processes. Standard renewal or master equation approach are typically used to study steady state and other transport properties such as average, mean squared displacement etc. What remains less explored is the two time point correlation functions whose evaluation is often daunting since it requires the implementation of the exact time dependent probability density functions of the resetting processes which are unknown for most of the problems. We adopt a different approach that allows us to write a stochastic solution for a single trajectory undergoing resetting. Moments and the autocorrelation functions between any two times along the trajectory can then be computed directly using the laws of total expectation. Estimation of autocorrelation functions turns out to be pivotal for investigating the ergodic properties of various observables for this canonical model. In particular, we investigate two observables (i) sample mean which is widely used in economics and (ii) time-averaged-mean-squared-displacement (TAMSD) which is of acute interest in physics. We find that both diffusion and drift-diffusion processes with resetting are ergodic at the mean level unlike their reset-free counterparts. In contrast, resetting renders ergodicity breaking in the TAMSD while both the stochastic processes are ergodic when resetting is absent. We quantify these behaviors with detailed analytical study and corroborate with extensive numerical simulations. Our results can be verified in experimental set-ups that can track single particle trajectories and thus have strong implications in understanding the physics of resetting. KW - autocorrelations KW - ergodicity KW - diffusion KW - stochastic resetting Y1 - 2022 U6 - https://doi.org/10.1088/1751-8121/ac4ce9 SN - 1751-8113 SN - 1751-8121 VL - 55 IS - 10 PB - IOP Publ. Ltd. CY - Bristol ER -