## Simulation of the continuous-time random walk using subordination schemes

Danhua Jiang, Yuanze Hong, and Wanli Wang ∗

School of Mathematical Sciences, Zhejiang University of Technology, Hangzhou 310023, China

The continuous time random walk model has been widely applied in various fields, including physics, biology, chemistry, finance, social phenomena, etc. In this work, we present an algorithm that utilizes a subordinate formula to generate data of the continuous time random walk in the long time limit. The algorithm has been validated using commonly employed observables, such as typical fluctuations of the positional distribution, rare fluctuations, the mean and the variance of the position, and breakthrough curves with time-dependent bias, demonstrating a perfect match.

## I. INTRODUCTION

The continuous-time random walk (CTRW) [1-7] is a stochastic process for a random walk that jumps from one position to another one. As an extension of the discrete random walk process, it was originally discussed by Montroll and Weiss [2], where a random walk is subordinated to a renewal process. Waiting times and displacements of the random particles are independent and identically distributed (IID) [4]. For the waiting times, it may be drawn from a heavy or narrow-tailed distribution. Different combinations of waiting times and displacements distributions yield different types of diffusion statistics, ranging from sub-diffusion, normal diffusion, and super diffusion. The CTRW model has attracted great attention from researchers due to its wide applications. For example, it was used in many physical and chemical phenomena, including the transport of amorphous materials [8], contamination in disorder systems [9], the U.S. dollar in Deutsche mark future exchange [10], diffusion of polymers in attractive nano-particle polymer mixtures [11] and so on.

The key issue is to discuss and simulate the positional distribution of the CTRW model. Two alternative approaches are the fractional diffusion equation [4, 12] and Langevin equations [13-15]. For example, Fogedby used decouple Langevin equations concerning power law step and waiting time distributions [16] to illustrate the CTRW model. One advantage is that the drift force field of the system is natural to be added. This method was further extended in Ref. [17, 18]. See also discussions of coupled cases [19]. Here we use the statistics of the number of renewals in the long time limit, eliminating the need to generate waiting times statistics.

This manuscript focuses on the case when waiting time has a finite mean and an infinite variance, and the displacement follows a Gaussian distribution with a nonzero mean. However, the methods we provided are much more general. Due to the asymmetric displacement distribution, sometimes it was called the biased CTRW model. The biased CTRW has garnered significant attention and a vast number of specialized literatures [20-25].

∗ contact author: wanliiwang@163.com

When the bias of the system or the mean of the displacement is zero, the process shows normal diffusion. While, when the bias is added, CTRW shows the enhanced or supper diffusion [1], i.e.,

$$\langle x ^ { 2 } \rangle - \langle x \rangle ^ { 2$$

with 1 &lt; α &lt; 2. In the long time limit, the positional distribution follows the asymmetric L´vy stable law [26] e according to the Montroll-Weiss equation [3, 4]. The fat tail is a key characteristic of the asymmetric L´vy dise tribution. Generally, when seeking data on CTRW, we obtain a random count representing the number of renewals occurring within the time interval from zero to t , then use subordinated processes to get the particle's position. A drawback is that it takes a long time to find the number of renewals when the observation time is long. Thus, it would be interesting to find a new way to generate the number of renewals and the position of particles.

Recall that in the long time limit, the number of renewals follows the asymmetric L´vy stable law [27]. e By instinct, one possible way is to generate the number of renewals based on the mentioned L´ evy stable law. The L´ evy stable law is perfectly correct, but it does have some drawbacks. For example, the far tail of L´vy distribue tion tends to infinity, which indicates that the number of the renewals tends toward negative infinity. In turn, the L´ evy stable law gives an infinite mean square displacement [27, 28] (MSD) of the number of renewals. It is certainly not true, as the far tail should have a cutoff [28] for a finite observation time t . In addition, as suggested in [28], the MSD of the number is governed by an infinite density. On the other hand, the positional distribution in the context of simulations converges to the L´ evy stable law rather slowly, especially when the system's bias is weak. In some sense, the convergence problem was solved in [29] using the fractional advectiondiffusion-asymmetry equation. To be more exact, the solution of the equation is the convolution of the L´vy stae ble law and a Gaussian distribution. Unfortunately, the solution of the fractional advection-diffusion-asymmetry equation yields a divergent MSD [29, 30].

Our goal in this manuscript is to construct a PDF to describe the number of renewals, which is not only valid for the typical fluctuations of the number of renewals but also effective for rare events. The main idea is that we use

the L´vy stable law and the non-normalized density to ile lustrate the central part and the far tail of the number of renewals, respectively. Below, we will demonstrate that these two laws exhibit distinct scaling limits on the number of renewals and then utilize commonly employed observables to illustrate the effectiveness of the constructed PDF.

The organization of the manuscript is as follows. In Sec. II, we give the definition of the CTRW model and the algorithm to generate the number of renewals in the long time limit. Specifically, the modified L´ evy stable law is introduced. We further check the theoretical prediction using the widely used observables in Sec. III, including typical fluctuations of the position, large deviations of the position describing the far tail, the mean and the variance of the position, and breakthrough curves with the time-dependent bias. Finally, we conclude the manuscript with a discussion in Sec. IV.

## II. MODEL AND ALGORITHM

Below, we present the definition of the CTRW model and the algorithm for simulating data on the number of renewals and positions.

## A. Model

In this subsection, we introduce the CTRW model, which is defined as follows. A random walk starts from its initial position. The particle waits at its initial position for time τ 1 drawn from φ τ ( ), and then makes a jump to x 1 with x 1 following f ( χ ). The particle traps at x 1 for time τ 2 and then jumps to the position x 1 + x 2 with x 2 generated from f ( χ ). Then, the process is renewed. All along the manuscript, the waiting times τ i are IID random variables, and the same holds for χ i . The number of renewals, denoted as N t , obeys the relation ∑ N t i =1 τ i + B t = , where t B t is the backward recurrence [27, 28] and t is the total observation time.

All along the manuscript, waiting times follow a power law distribution [28, 31, 32]

$$\begin{array} {$$

with 1 &lt; α &lt; 2. As mentioned, waiting times τ have a finite mean and an infinite variance. Performing the Laplace transform with respect to φ τ ( ), we have ̂ φ s ( ) = ∫ ∞ 0 exp( -st φ τ ) ( ) dτ ∼ 1 -〈 〉 τ s +( τ 0 ) α | Γ(1 -α ) | s α with 〈 τ 〉 = ατ / α 0 ( -1) being the mean of waiting times.

For displacements, we focus on Gaussian distribution with a finite mean a and a finite variance σ 2

$$f ( x ) = \frac { 1 } { \sqrt { 2 \sigma ^ { 2$$

Let Q N t ( ) be the PDF of the number of renewals and f ( x N | ) be the conditional probability of finding the particle on the position x conditioned it made N jumps. The PDF of the position follows

$$P ( x, t ) & = \sum _ { N = 0 } ^ { \infty }$$

Eq. (3) describes the particle at position x after N steps on conditional that exactly N steps were made up to time t . Eq. (3) can be interpreted as the so-called integral formula of subordination [3, 29, 33, 34], i.e., a random walk is subordinated to a renewal process. The subordination scheme has been employed in mathematics and physics within various contexts, including fractional kinetic equations [31, 35, 36], Jacobi stochastic volatility models [37], and financial models [38]. Below, we will use Eq. (3) to generate the data of positions.

Recall that displacements are IID random variables, thus f ( x N | ) obeys

$$f ( x | N ) = \frac { 1 } { \sqrt { 2 \sigma ^$$

As described in [26-28], in the long time limit, the number of renewals follows the asymmetric L´vy distribution e

$$Q _ { t } ( \xi ) \sim \mathcal { L } _ { \alpha } ( \xi )$$

with ξ = ( N -t/ τ 〈 〉 ) / t/t ( ¯ ) 1 /α and ¯ = t 〈 τ 〉 1+ α / τ ( α 0 | Γ(1 -α ) ). | Eq. (5) is plotted by the red solid line in Fig. 1. Here the symbol L α ( ξ ) is the L´ evy stable distribution, whose Fourier transform is ∫ ∞ -∞ exp( -ikξ ) L α ( ξ dx ) = exp(( -ik ) α ). Typical fluctuations Eq. (5) describes the case when N -t/ τ 〈 〉 is of the order of t 1 /α , i.e., the central part of the distribution. Using the asymptotic behavior of the L´ evy distribution, the far left tail of N follows Q ξ t ( ) ∼ -( ξ ) - -1 α / Γ( -α ). It indicates that Eq. (5) does not give the information of the MSD [28].

## B. Algorithm

According to Eqs. (3, 4, 5), in the long time limit, the positional distribution follows

$$P ( x, t ) \sim & \frac { 1 } { ( t / \overline { t } ) ^ { 1 / \alpha } } \int _ { 0 } ^ { \infty } \mathcal { L } _ { \alpha } \left ( \frac { N - t / \langle \tau \rangle } { ( t / \overline { t } ) ^ { 1 / \alpha } } \right ) \\ & \frac { \exp \left ( - \frac { ( x - a N ) ^ { 2 } } { 2 \sigma ^ { 2 } N } \right ) } { \sqrt { 2 \pi \sigma ^ { 2 } N } } \mathrm d N,$$

FIG. 1: Distribution of the number of renewals. The red solid line describes typical fluctuations of the number of renewals, shown by the L´ evy law Eq. (5) and rare events Eq. (9) are plotted by the black dashed line. Symbols plotted by circles are simulations of the renewal process for 2 × 10 7 realizations. We choose t = 1000 , τ 0 = 0 1, and . α = 3 2. /

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00299v2-figure-9b248851-ab36-4ea4-8ed4-5fd91324e20a.jpg)

Changing the random variables, the above equation leads to

$$\stackrel { \text{----} } { \text{----} } \circ \, \stackrel { \text{----} } { \text{----} } \, \stackrel { \text{----} } { \text{----} } \, \stackrel { \text{----} } { \text{----} } \, \stackrel { \text{----} } { \text{----} } \, \stackrel { \text{----} } { \text{----} } \, \stackrel { \text{----} } { \text{----} } \, \stackrel { \text{----} } { \text{----} } \, \stackrel { \text{----} } { \text{----} } \, \stackrel { \text{----} } { \text{----} } \, \stackrel { \text{----} } { \text{----} } \, \stackrel { \text{----} } { \text{----} } \, \stackrel { \text{--} } { \text{----} } \, \stackrel { \text{----} } { \text{----} } \, \stackrel { \text{----} } { \text{----} } \, \stackrel { \text{----} } { \text{----} } \, \stackrel { \text{----} } { \text{----} } \, \stackrel { \text{----} } { \text{----} } \, \stackrel { \text{----} } { \text{----} } \, \stackrel { \text{----} } { \text{----} } \, \stackrel { \text{----} } { \text{----} } \, \stackrel { \text{----} } { \text{----} } \, \stackrel { \text{----} } { \text{----} } { \text{----} } \, \stackrel { \text{----} } { \text{----} } \, \stackrel { \text{----} } { \text{----} } \, \stackrel { \text{----} } { \text{----} } \, \stackrel { \text{----} } { \text{----} } \, \stackrel { \text{----} } { \text{----} } \, \stackrel { \text{----} } { \text{----} } \, \stackrel { \text{----} } { \text{----} } \, \stackrel { \text{----} } { \text{----} } \, \stackrel { \text{----} } { \text{----} } \, \stackrel { \text{----} } { \end{ \text{----} } { \text{----} } \, \textrel{ \text{----} } { \text{----} } \, \textrel{ \text{----} } { \text{----} } \, \textrel{ \text{----} } { \text{----} } \, \textrel{ \text{----} } { \text{----} } \, \textrel{ \text{----} } { \text{----} } \, \textrel{ \text{----} } { \text{----} } \, \textrel{ \text{----} } { \text{----} } \, \textrel{ \text{----} } { \text{----} } \, \textrel{ \text{----} } { \text{----} } \, \textrel{ \text{----} } { \text{----} } \, \textrel{ \text{----} } { \end{ \text{----} } { \text{----} } \, \textrel{ \text{----} } { \text{----} } \, \textrel{ \text{----} } { \text{----} } \, \textrel{ \text{----} } { \text{----} } \, \textrel{ \text{----} } { \text{----} } \, \textrel{ \text{----} } { \text{----} } \, \textrel{ \text{----} } { \text{----} } \, \textrel{ \text{----} } { \text{----} } \, \textrel{ \text{----} } { \text{----} } \, \textrel{ \text{----} } { \text{----} } \, \textrel { \text{----} } { \text{----} } \, \textrel{ \text{----} } { \text{----} } \, \textrel{ \text{----} } { \text{----} } \, \textrel{ \text{----} } { \text{----} } \, \textrel{ \text{----} } { \text{----} } \, \textrel{ \text{----} } { \text{----} } \, \textrel{ \text{----} } { \text{----} } \, \textrel{ \text{----} } { \text{----} } \, \textrel{ \text{----} } { \text{----} } \, \textrel{ \text{----} } { \text{----} } \, \textrel{ \text{----} } { \textrel{ \text{----} } \, \textrel{ \text{----} } { \text{----} } \, \textrel{ \text{----} } { \text{----} } \, \textrel{ \text{----} } { \text{----} } \, \textrel{ \text{----} } { \text{----} } \, \textrel{ \text{----} } { \text{----} } \, \textrel{ \text{----} } { \text{----} } \, \textrel{ \text{----} } { \text{----} } \, \textrel{ \text{----} } { \text{----} } \, \textrel{ \text{----} } { \text{----} } \, \textrel{ \text{----} } { \text{----} } \, \ \textrel{ \text{----} } { \text{----} } \, \textrel{ \text{----} } { \text{----} } \, \textrel{ \text{----} } { \text{----} } \, \textrel{ \text{----} } { \text{----} } \, \textrel{ \text{----} } { \text{----} } \, \textrel{ \text{----} } { \text{----} } \, \textrel{ \text{----} } { \text{----} } \, \textrel{ \text{----} } { \text{----} } \, \textrel{ \text{----} } { \text{----} } \, \textrel{ \text{----} } { \text{----} } \, \textrel{ \text{----} \, \textrel{ \text{----} } { \text{----} } \, \textrel{ \text{----} } { \text{----} } \, \textrel{ \text{----} } { \text{----} } \, \textrel{ \text{----} } { \text{----} } \, \textrel{ \text{----} } { \text{----} } \, \textrel{ \text{----} } { \text{----} } \, \textrel{ \text{----} } { \text{----} } \, \textrel{ \text{----} } { \text{----} } \, \textrel{ \text{----} } { \text{----} } \, \textrel{ \text{----} } { \text{----} } \, \textrel{ \text{--} } { \text{----} } \, \textrel{ \text{----} } { \text{----} } \, \textrel{ \text{----} } { \text{----} } \, \textrel{ \text{----} } { \text{----} } \, \textrel{ \text{----} } { \text{----} } \, \textrel{ \text{----} } { \text{----} } \, \textrel{ \text{----} } { \text{----} } \, \textrel{ \text{----} } { \text{----} } \, \textrel{ \text{----} } \, \textrel{ \text{----} } { \text{----} } \, \textrel{ \text{----} } { \text{----} } \, \textrel{ \text{----} } { \text{----} } \, \textrel{ \text{----} } \, \textrel{ \text{----} } { \text{----} } \, \textrel{ \text{----} } { \text{----} } \, \textrel{ \text{----} } { \text{----} } \, \textrel{ \text{----} } { \text{----} } \, \textrel{ \text{----} } { \text{----} } \, \textrel{ \text{----} } { \text{----} } \, \textrel{ \text{----} } { \text{----} } \, \textrel{ \text{----} } \, \textrel{ \text{----} } { \text{----} } \, \textrel{ \text{----} } { \text{----} } \, \textrel{ \text{---- } } { \text{----} } \, \textrel{ \text{----} } { \text{----} } \, \textrel{ \text{----} } { \text{----} } \, \textrel{ \text{----} } { \text{----} } \, \textrel{ \text{----} } { \text{----} } \, \textrel{ \text{----} } { \text{----} } \, \textrel{ \text{----} } { \text{----} } \, \textrel{ \text{----} } { \text{----} } \, \textrel{ \text{----} } { \text{----} } \, \textrel{ \text{----} } { \text{----} } \, \textrel{ \text{----} } \, \textrel{ \text{----} } { \text{----} } \, \textrel{ \text{----} } \, \textrel{ \text{----} } { \text{----} } \, \textrel{ \text{----} } { \text{----} } \, \textrel{ \text{----} } { \text{----} } \, \textrel{ \text{----} } \, \textrel{ \text{----} } { \text{----} } \, \textrel{ \text{----} } { \text{----} } \, \textrel{ \text{----} } { \text{----} } \, \textrel{ \text{----} } { \text{----} } \, \textrel{ \text{----} } \, \textrel{ \text{----} } { \text{----} } \, \textrel{ \text{----} } { \text{----} } \, \textrel{ \text{----} } { \text{----} } \, \textrel{ \text{----} } { \text{----} } \, \textrel{ \text{----} } { \text{----} } \, \textrel{ \text{----} } { \text{----} } \, \textrel{ \text{---- } } { \textrel{ \text{----} } \, \textrel{ \text{----} } { \text{----} } \, \textrel{ \text{----} } { \text{----} } \, \textrel{ \text{----} } { \text{----} } \, \textrel{ \text{----} } { \text{----} } \, \textrel{ \text{----} } { \text{----} } \, \textrel{ \text{----} } { \text{----} } \, \textrel{ \text{----} } { \text{----} } \, \textrel{ \text{----} } { \text{----} } \, \textrel{ \text{----} } { \text{----} } \, \textrel{ \text{----} } { \text{--} } \, \textrel{ \text{----} } { \text{----} } \, \textrel{ \text{----} } { \text{----} } \, \textrel{ \text{----} } { \text{----} } \, \textrel{ \text{----} } { \text{----} } \, \textrel{ \text{----} } \, \textrel{ \text{----} } { \textrel{ \text{----} } \, \textrel{ \text{----} } \, \textrel{ \text{----} } { \textrel{ \text{----} } \, \textrel{ \text{----} } \, \textrel{ \text{----} } { \textrel{ \text{----} } \, \textrel{ \text{----} } \, \textrel{ \text{----} } { \textrel{ \text{----} } \, \textrel{ \text{----} } { \textrel{ \text{----} } \, \textrel{ \text{----} } { \textrel{ \text{----} } \, \textrel{ \text{----} } { \textrel{ \text{----} } \, \textrel{ \text{----} } \, \textrel{ \text{----} } { \textrel{ \text{----} } \, \textrel{ \text{----} } \, \textrel{ \text{----} } { \text{----} } \, \textrel{ \text{----} } { \text{----} } \, \textrel{ \text{----} } { \textrel{ \text{----} } \, \textrel{ \text{----} } { \textrel{ \text{----} } \, \textrel{ \text{----} } \, \textrel{ \text{----} } { \text{----} } \, \textrel{ \text{----} } { \textrel{ \text{----} } \, \textrel{ \text{----} } { \textrel{ \text{----} } \, \textrel{ \text{----} } \, \textrel{ \text{----} } \, \textrel{ \text{----} } { \textrel{ \text{----} } \, \textrel{ \text{----} } { \textrel{ \text{----} } \, \textrel{ \text{----} } { \text{----} } \, \textrel{ \text{----} } { \text{----} } \, \textrel{ \text{----} } { \text{----} } \, \textrel{ \text{----} } { \text{----} } \, \textrel{ \text{----} } { \textrel{ \text{----} } \, \textrel{ \text{----} } { \text{----} } \, \textrel{ \text{----} } { \textrel{ \text{----} } \, \textrel{ \text{----} } { \textrel{ \text{----} } \, \textrel{ \text{----} } { \textrel{ \text{----} } \, \textrel{ \text{----} } { \textrel{ \text{----} } \, \textrel{ \text{----} } { \textrel{ \text{----} } \, \textrel{ \text{----} } { \textrel{ \text{----} } \, \textrel{ \text{----} }$$

Rare fluctuations of N were investigated, describing the scaling when N -t/ τ 〈 〉 is of the order of t . According to the result given in [28], large deviations of N obey

Eq. (7) is the main result of this section used to generate the random variables. Note that L α ( ξ ) in Eq. (7) should be modified, denoted as L ∗ α ( ξ ). Otherwise, it will not accurately account for the MSD and the far tail of the position distribution. For that, we generate the statistics of L ∗ α ( ξ ), stemming from typical fluctuations and the infinite density [39] and then use it to get the data of the position. Afterwards, positions of particles are drawn from a Gaussian distribution with the mean at/ τ 〈 〉 + aξ t/t ( ¯ ) 1 /α and the variance δ 2 ( t/ τ 〈 〉 + ( ξ t/t ¯ ) 1 /α ) 2 .

$$Q _ { t } ( \epsilon ) & \sim \alpha ( \tau _ { 0 } ) ^ { \alpha } \left [ t ( - \epsilon \langle \tau \rangle ) ^ { - \alpha - 1 } + \frac { 1 - \alpha } { \alpha } ( - \epsilon \langle \tau \rangle ) ^ { - \alpha } \right ] ^ { 2 } ( 8 ) ^ { 2 } \quad \text{As} \quad \text{parts}, \\ \dots \quad \text{--} \quad \dots \quad \text{--} \quad \text{--} \quad \text{--} \quad \text{--} \quad \text{--} \quad \text{--} \quad \text{--} \quad \text{--} \quad \text{--} \quad \text{--} \quad \text{--} \quad \text{--} \quad \text{--} \quad \text{--} \quad \text{--} \quad \text{--} \quad \text{--} \quad \text{--} \quad \text{--} \quad \text{--} \quad \text{--} \quad \text{--} \quad \text{--} \quad \text{--} \quad \text{--} \quad \text{--} \quad \text{--} \quad \text{--} \quad \text{--} \quad \end{cases}$$

with /epsilon1 = N -t/ τ 〈 〉 . Consider the random variable, ξ = /epsilon1/ ( t/t ¯ ) 1 /α , we have

$$Q _ { t } ( \xi ) & \sim \alpha ( \tau _ { 0 } ) ^ { \alpha } \left ( \frac { t } { t } \right ) ^ { 1 / \alpha } & \quad \text{for $p$} \inf i n \text{if $are dr} \\ & \quad \times \left [ \frac { t } { ( - \xi ( \frac { t } { t } ) ^ { 1 / \alpha } \langle \tau \rangle ) ^ { \alpha + 1 } } + \frac { ( 1 - \alpha ) / \alpha } { ( - \xi ( \frac { t } { t } ) ^ { 1 / \alpha } \langle \tau \rangle ) ^ { \alpha } } \right ]. & \quad \text{begin{array} { c } \text{$P$} \\ \text{$Bell$} \\ \text{ables} } \end{array}$$

See the black dashed line in Fig. 1. Recall that Eq. (9) is valid for small ξ corresponding to small N . As shown in Fig. 1, Eq. (5) is valid for the central part of the number of renewals together with the right tail, and the non-normalized density Eq. (9) works for the left far tail. By instinct, we can utilize the effective part of Eqs. (5) and (9) to generate data of the number of renewals.

In order to generate the random variable based on Eq. (9), we construct a PDF by adding a condition z on Eq. (9),

$$\Big ] \quad & Q ^ { I } _ { t } ( \xi ) = \alpha ( \tau _ { 0 } ) ^ { \alpha } \left [ \frac { \bar { t } } { ( - \xi \langle \tau \rangle ) ^ { \alpha + 1 } } + \frac { ( 1 - \alpha ) / \alpha } { ( - \xi ( \tau \rangle ) ^ { \alpha } ( \frac { t } { t } ) ^ { 1 - \frac { 1 } { \alpha } } } \right ] \ ( 1 0 ) \\ & _ { - }$$

with ¯ b &lt; ξ ≤ z &lt; 0. Here, z is an undetermined constant. As N ≥ 0, the bottom limit of ξ should be great than -t 1 -1 /α ¯ t 1 /α / τ 〈 〉 , i.e.,

$$\bar { b } = - \left ( \frac { \alpha } { \Gamma ( 2 - \alpha ) } \right ) ^ { 1 / \alpha } \left ( \frac { t } { \tau _ { 0 } } \right ) ^ { 1 - \frac { 1 } { \alpha } }. \quad \ \ ( 1 1 )$$

When 1 &lt; α &lt; 2, in the long time limit, ¯ b →-∞ . While for a finite observation time t , ξ can not tend to infinity. Thus, the modified L´vy stable law e L ∗ α ( ξ ) should have a cutoff at the tail.

We further introduce the cumulative distribution function, describing the probability of ξ when ξ is less than y

$$h ( y ) = \int _ { \bar { b } } ^ { y } Q _ { t } ^ { I } ( \xi ) \mathrm d \xi. \quad \text{(1)} \\ -$$

Performing the integral from ¯ to b z , Eq. (12) leads to

$$\stackrel { \_ \_ } { \text{ot} } \quad \frac { \tau _ { 0 } z ( ( \alpha - 1 ) t ) ^ { 1 / \alpha } | \Gamma ( 1 - \alpha ) | ^ { 1 / \alpha } + t ( \alpha \tau _ { 0 } ) ^ { 1 / \alpha } } { t ( \alpha \tau _ { 0 } ) ^ { 1 / \alpha } ( - z ) ^ { \alpha } | \Gamma ( 1 - \alpha ) | } = 1. \quad ( 1 3 )$$

Eq. (13) can be used to determine the upper limit z defined in Eq. (10). When α = 3 2, using a cubic equa-/ tion, the formula for the roots of Eq. (13) is easy to find. While, if α = 5 4, Eq. (13) reduces to a 5-th degree poly-/ nomial. Thus, according to the Abel-Ruffini theorem, the equation concerning z shown by Eq. (13) is generally unsolvable. For that, to have a universal program valid for all kinds of α , we have to abandon some methods that are z -dependent.

As shown in Fig. 2, the modified L´vy stable has three e parts, labeled as A B , , and C , respectively. Part A represents the probability for the normalization, describing the fade probability of random variables equaling to ¯ . b For part B , we utilize the information of the modified infinite density Eq. (10). In part C , random variables ξ are drawn from the L´ evy stable law Eq. (5). The dot D is the intersection of the part B and the part C .

Below, we calculate the probability of random variables located in regions B and C, respectively. For the

FIG. 2: Plot of the modified asymmetric L´vy stable distribue tion. Here, the black dot marked by A describes the modified survival probability Eq. (16), the part B stems from the modified infinite density Eq. (10), the part C is related to typical fluctuations Eq. (5), and the black dot D denoted as ξ ∗ is the crossover point between typical fluctuations and rare fluctuations.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00299v2-figure-8ef4249a-d0ae-448f-941e-a1fedea97213.jpg)

probability of rare fluctuations in region B, we have

$$P _ { \text{rare} } ( \xi ^ { * } ) & = h ( \xi ^ { * } ) \\ & = \frac { \tau _ { 0 } z ( ( \alpha - 1 ) t | \Gamma ( 1 - \alpha ) | ) ^ { 1 / \alpha } + t ( \alpha \tau _ { 0 } ) ^ { 1 / \alpha } } { t ( \alpha \tau _ { 0 } ) ^ { 1 / \alpha } ( - \xi ^ { * } ) ^ { \alpha } | \Gamma ( 1 - \alpha ) | }, \quad \text{In} \ \ \text{ obtain} \ \.$$

as mentioned before, ξ ∗ is the transition point of the region B and the region C. Recall that part C is illustrated by the L´ evy stable law, the probability in the region C follows

$$P _ { \text{typical} } ( \xi ^ { * } ) = \int _ { \xi ^ { * } } ^ { \infty } \mathcal { L } ( \xi ) d \xi. \quad \quad ( 1 5 ) \\ \quad \. \. \. \.$$

When t → ∞ , P rare ( ξ ∗ ) + P typical ( ξ ∗ ) tends to unity. While for a finite t , P rare ( ξ ∗ ) + P typical ( ξ ∗ ) &lt; 1. Therefore, we further introduce the modified or fake survival probability describing the non-traveling particles

$$P _ { \text{survival} } ( \xi ^ { * } ) = 1 - P _ { \text{rare} } ( \xi ^ { * } ) - P _ { \text{typical} } ( \xi ^ { * } ). \quad ( 1 6 ) \quad \text{tiona}$$

In other words, it is related to scenarios where the number of renewals is zero, which leads to the non-moving particles in the context of the CTRW model.

Now we discuss the details of generating random variables ξ drawn from L ∗ α ( ξ ). Firstly, we generate a random variable drawn from a uniform distribution, denoted as η . Based on the value of η , we determine which part it belongs to. If η &gt; P survival ( ξ ∗ ) + P rare ( ξ ∗ ), we generate a random variable ξ from L ( ξ ) under the condition that ξ ≥ ξ ∗ . Instructions for generating a L´vy stable distrie bution can be found on [40, 41] or the MATLAB Central File Exchange, as detailed in the work by Veillette [42]. When P survival ( ξ ∗ ) &lt; η &lt; P survival ( ξ ∗ ) + P rare ( ξ ∗ ), ξ are drawn from the modified infinity density Eq. (10) using the accept/reject algorithm. If η &lt; P survival ( ξ ), the random variable is fixed, i.e., ξ = ¯ . b Note that the MSD is sensitive to lagging particles far behind the mean [43]. Therefore, when simulating the position for the MSD, the random variables η are larger than P survival ( ξ ∗ ) using the theory given in Eq. (10). Throughout the manuscript, the theories marked in legends of figures refer to the algorithm discussed in this subsection. Specifically, we generate ξ drawn from the modified L´ evy law L ∗ α ( ξ ) and then use it to determine the position statistics given by the expression

$$J ( x ) = \frac { \exp \left ( - \frac { ( x - a \frac { t } { ( \tau ) } - a \xi ( \frac { t } { t } ) ^ { 1 / \alpha } ) ^ { 2 } } { 2 \sigma ^ { 2 } ( t / \langle \tau \rangle + \xi ( t / t ) ^ { 1 / \alpha } ) } \right ) } { \sqrt { 2 \sigma ^ { 2 } \pi ( \frac { t } { \langle \tau \rangle } + \xi ( \frac { t } { t } ) ^ { 1 / \alpha } ) } }. \\ \text{e. } J ( x ) \text{ represents the distribution of the positio}$$

Here, J x ( ) represents the distribution of the position x at the time t within the context of the CTRW model, where ξ are drawn from L ∗ α ( ξ ).

## III. THEORETICAL PREDICTIONS COMPARED WITH ASYMPTOTIC BEHAVIORS AND SIMULATIONS

In this section, we generate the data of the position obtained from the method outlined in Sec. II B to obtain typical fluctuations of the position, rare fluctuations, the mean of position, the variance of the position, and breakthrough curves. Additionally, we generate the positional distribution data using simulations and outline the wellknown asymptotic behaviors for comparative analysis.

## A. Typical fluctuations

Several authors have studied typical fluctuations of the position for the biased CTRW; see Refs. [21, 26, 29, 43, 44]. In the long-time approximation, the scaled positional distribution follows the asymmetric L´ evy stable distribution [26, 45]

$$P ( \zeta ) \sim \mathcal { L } _ { \alpha } ( \zeta ) \text{ \quad \ \ } \text{(1)}$$

with ζ = ( x -at/ τ 〈 〉 ) /l ( ) t 〉 and l ( ) t = a t/t ( ) 1 /α . In Eq. (17), it is crucial to emphasize that a must be nonzero; otherwise, the positional distribution would yield a Gaussian distribution. In the general scenario, P x, t ( ) is described by the convolution of a L´vy stable e

distribution and a Gaussian distribution [29]

$$& \text{distribution and a Gaussian distribution [29]} \\ & P ( x, t ) \sim \int _ { - \infty } ^ { \infty } \mathcal { L } _ { \alpha } ( \xi ) \frac { \exp \left ( - \frac { ( x - \frac { a t } { ( \tau ) } - a \xi ( \frac { t } { t } ) ^ { 1 / \alpha } ) ^ { 2 } } { 2 \sigma ^ { 2 } \frac { t } { ( \tau ) } } \right ) } { \sqrt { 2 \pi \sigma ^ { 2 } \frac { t } { ( \tau ) } } } d \xi, \\ & \text{describing typical fluctuations when $x-a t/\langle\tau\rangle \sim t^{1/\alpha}$,} \quad \widehat { \underline { s } } }$$

describing typical fluctuations when x -at/ τ 〈 〉 ∼ t 1 /α , i.e., the central part of the positional distribution. As shown in Fig. 3, Eq. (18) fails for the left tail of the positional distribution. While the theoretical prediction given in the manuscript is valid for both tails. It indicates that if we use the modified L´ evy stable law in Eq. (7), Eq. (7) does work perfectly. In addition, we plot the positional distribution at a short time t , see Figs. 5 and 6. For this case, the positional distribution follows a nearly Gaussian distribution, affirming the validity of our findings. All along the manuscript, the theoretical predictions are based on the algorithm described in Sec. II B, which uses the modified L´vy stable law to e generate statistics of the number of renewals.

FIG. 3: Distribution of the position for different biases. Here, the red symbols describe simulations and the green symbols indicate theoretical predictions. For the theories, we use the algorithm provided in Sec. II B. The solid lines stem from Eq. (18) describing typical fluctuations of x . We choose α = 3 2 / , σ = 1 , t = 1000, and τ 0 = 0 1. .

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00299v2-figure-0e1d4527-bfe6-4a53-af04-1c3d66e6a21f.jpg)

## B. Rare fluctuations

To proceed, we discuss the limit of the Montroll-Weiss equation (20) but when | s / k | | | is fixed. Specifically, s and k are small and in the same order. As a complementary to the L´vy law Eq. (17), the rare fluctuations illustrate e the far tail of the positional distribution. Following the result given in [43], rare fluctuations are

$$P ( x, t ) \sim \frac { ( \tau _ { 0 } ) ^ { \alpha } } { a t ^ { \alpha } } ( \alpha ( 1 - \frac { x / a } { t / \langle \tau \rangle } ) ^ { - 1 - \alpha } - ( \alpha - 1 ) ( 1 - \frac { x / a } { t / \langle \tau \rangle } ) ^ { - \alpha } ), \quad \text{We} \\ ( 1 9 ) \quad \langle x ( t ) \rangle.$$

FIG. 4: Same as in Fig. 3 in linear scale.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00299v2-figure-bb603943-3a5b-443d-aa97-97786402c7fa.jpg)

FIG. 5: Positional distribution with a short observation time t , i.e., t = 2 10, and 50. , The symbols are theoretical predictions, and the corresponding simulations are plotted using the solid lines. Here we choose α = 3 2 / , τ 0 = 0 1 . , δ = 1 , a = 1, and 2 × 10 7 realizations.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00299v2-figure-5ae0a394-7c8b-4015-8aa2-dd5bf3738aca.jpg)

describing the scaling when x is of the order of the time t . Eq. (19) gives a prediction for particles that are near the initial position. In other words, the mentioned slowingmoving particles stem from these longest single waiting times drawn from Eq. (1). As shown in Fig. 7, our theory is effective for both the positional distribution's central part and far tails, whereas the L´vy stable law Eq. (17) e tends to infinity.

## C. Moments of the position

We investigate the mean of the position denoted by 〈 x t ( ) 〉 , which refers to the sum of all positions divided by

FIG. 6: Same as in Fig. 5 in linear scale.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00299v2-figure-b867686c-9fdd-490a-a8fe-ec95b4e45c9a.jpg)

FIG. 7: Simulations (' /square ') of the positional distribution are compared to rare fluctuations Eq. (19) (solid line), the L´vy e law Eq. (17) in terms of x (dashed line), and theory (' ◦ '). The parameters are a = 1, σ = 1, t = 1000, α = 3 2, / τ 0 = 0 1, . and 2 × 10 7 realizations.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00299v2-figure-061e2b82-67a1-4097-bd17-4aa6e1dd47df.jpg)

the number of particles. It is a measure that reflects the tendency of the position. Let us start from the MontrollWeiss equation, which in Fourier-Laplace space is as follows [3, 4]

$$\widetilde { \widehat { P } } ( k, s ) = \frac { 1 - \widehat { \phi } ( s ) } { s } \frac { 1 } { 1 - \widetilde { f } ( k ) \widehat { \phi } ( s ) }. \quad \ ( 2 0 )$$

FIG. 8: Mean of the position under various biases are depicted. Lines correspond to the asymptotic prediction given by Eq. (21). The parameters are σ = 1, α = 3 2, / τ 0 = 0 1, . and 2 × 10 7 realizations.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00299v2-figure-b23149d1-7d6d-4632-a284-40c42c603f95.jpg)

FIG. 9: Same as Fig. 8 in log-log scale.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00299v2-figure-0051b4d4-6f4b-4321-8d2b-60ca0590fcff.jpg)

Using the relation

$$\langle \hat { x } ^ { q } ( s ) \rangle = i ^ { q } \frac { \partial ^ { q } \widetilde { \widehat { P } } ( k, s ) } { \partial k ^ { q } } | _ { k = 0 }, \\ \intertext { r time limit F a. } \text{$2$} \, \text{$7$} \, \text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text$$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\font$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text	$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text
$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text#$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text}$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text'$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text`$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text $\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text<$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\texts$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text1$$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\test$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text51$\text$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\textc$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text.$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$$$$

in the long time limit Eq. (20) yields

$$\langle x ( t ) \rangle \sim a \frac { t } { \langle \tau \rangle }.$$

Recall that the waiting time has a finite mean 〈 τ 〉 , then for a fixed t , the mean of the number of renewals is 〈 N t ( ) 〉 ∼ t/ τ 〈 〉 [27, 28]. Thus 〈 x t ( ) 〉 ∼ 〈/triangle x 〉〈 N t ( ) 〉 ∼ at/ τ 〈 〉 is obtained again, where 〈/triangle x 〉 is the mean of the displacements. As shown in Figs. 8 and 9, the theoretical prediction grows linearly with the observation time t , showing a perfect match with simulations and Eq. (21).

The MSD is the most commonly employed quantifier for diffusion processes. Due to the effect of the bias, the biased CTRW shows superdiffusion [22]. Fascinatingly, the MSD is controlled by these sluggish particles that significantly lag far behind the average position, as opposed to the swift particles near the average. In that sense, the rare events, governed by the infinite density, are responsible for the enhanced diffusion [43]. Mathematically, the variance of the position grows as

$$\text{Var} ( x ) & = \langle x ^ { 2 } \rangle - \langle x \rangle ^ { 2 } \\ & \sim \frac { 2 a ^ { 2 } \tau _ { 0 } ^ { \alpha } } { ( 2 - \alpha ) ( 3 - \alpha ) \langle \tau \rangle ^ { 3 } } t ^ { 3 - \alpha } + \sigma ^ { 2 } \frac { t } { \langle \tau \rangle }, \quad ( 2 2 ) \quad \text{and} \, t \\ \text{movement to the solution of ite potential extension in}$$

measuring the time evolution of its spatial extension. In the long time limit, the leading term of Eq. (22) is t 3 -α , and when the observation time is short, the correction term is of importance. As shown in Fig. 10, theoretical results are consistent with simulations. In particular, when the bias is weak, for example, a = 0 5, our the-. ory gives a perfect prediction if compared with Eq. (22). In linear-linear scale, small deviations exist between theoretical predictions and simulation results; see Fig. 11. The reason is as follows. When the bias is strong, particles are pushed by the bias and move quickly. In this scenario, the asymptotic L´vy stable law given by Eq. (17) e applies. Consequently, the modified infinite density becomes insignificant, as its probability is very low compared to the weak bias case (see Fig. (3)). Therefore, the algorithm provided in Sec. II B essentially utilizes Eq. (17) with a cutoff at the tail to compute the MSD. Besides, fluctuations of the far tail of the positional distribution are large when the bias are strong. It leads to deviations, but it does not like the results given by Eq. (18), which predicts an infinite MSD.

## D. Breakthrough curves

Now we further discuss the positional distribution, but from aspects of breakthrough curves, which are essential tools in measuring the contamination in porous materials [9, 46-49]. We assume that all the particles are injected immediately when we start to observe the process, wherein the positional concentration is measured at a fixed position x b . To be more exact, the breakthrough curves, the probability of particles being at x b at a specific observation time t , are measured in terms of the distribution.

Based on the discussions in [9, 29], we delve into the time-dependent bias to validate the program outlined

FIG. 10: Plot of the Var( x ). The solid lines describes the asymptotic behavior Eq. (22), together with simulations (' ◦ ') and theories ('+'). The parameters are the same as in Fig. 8.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00299v2-figure-fd19b793-9ffd-45ad-b47d-0bf2d16e6d9f.jpg)

FIG. 11: Same as in Fig. 10 in linear-linear scale.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00299v2-figure-ec06e6d6-9d78-4a5c-969e-6853d16f179e.jpg)

in Sec. II. The process entails four states. In state (i), particles are promptly injected; thereafter, displacements follow a Gaussian distribution with a mean of a and a variance of σ 2 from time zero until t a . All along the process, the variance of displacements remains the same. Subsequently, in state (ii), the bias of the system is enhanced from time t a to t b , i.e., the displacements are drawn from Gaussian distribution with the mean 4 aη/ (1 + η ), where η &gt; 2 / 3. State (iii) witnesses a bias of the system represented by a/ (2 / 3 + η ) over the time interval ( t b , t c ), constituting the weakest bias among the four states. Ultimately, the bias reverts to a after t c until the end of the observation time. The corresponding theoretical prediction was considered using the fractional advection-diffusion-asymmetry equation in Ref. [29], where the main strategy is that the final position of each state is regarded as the initial position of the next process. In the long time limit, the positional

with and

$$P ( x, t ) \sim \int _ { - \infty } ^ { \infty } \frac { 1 } { \sqrt { 4 \pi c _ { m 1 } } } \exp \left [ - \frac { ( x - y - c _ { m 2 } ) ^ { 2 } } { 4 c _ { m 1 } } \right ] \frac { 1 } { ( c _ { m 3 } ) ^ { 1 / \alpha } } \mathcal { L } _ { \alpha } \left [ \frac { y } { ( c _ { m 3 } ) ^ { 1 / \alpha } } \right ] \mathrm d y$$

$$c _ { m 1 } = \begin{cases} \begin{array} { c } t \frac { \sigma ^ { 2 } } { 2 ( \tau ) }, & m = 1 ; \\ t \frac { \sigma ^ { 2 } } { 2 ( \tau ) }, & m = 2 ; \\ t \frac { \sigma ^ { 2 } } { 2 ( \tau ) }, & m = 3 ; \\ t \frac { \sigma ^ { 2 } } { 2 ( \tau ) }, & m = 4, \end{array} \\ a _ { 1 } \frac { t } { ( \tau ) }, & \quad. & m = 1 ; \end{cases}$$

$$c _ { m 2 } = \begin{cases} \begin{array} { c } a _ { 1 } \frac { t } { ( \tau ) }, \\ [ a _ { 1 } t _ { 1 } + a _ { 2 } ( t - t _ { 1 } ) ] \frac { 1 } { ( \tau ) }, \\ [ a _ { 1 } t _ { 1 } + a _ { 2 } ( t _ { 2 } - t _ { 1 } ) + a _ { 3 } ( t - t _ { 2 } ) ] \frac { 1 } { ( \tau ) }, \\ [ a _ { 1 } t _ { 1 } + a _ { 2 } ( t _ { 2 } - t _ { 1 } ) + a _ { 3 } ( t _ { 3 } - t _ { 2 } ) + a _ { 4 } ( t - t _ { 3 } ) ] \frac { 1 } { ( \tau ) }, \ m = 4, \end{array} \\ \begin{array} { c } a ^ { \alpha } t \frac { 1 } { m } & m = 1. \end{array} \end{cases}$$

$$c _ { m 3 } = \begin{cases} \begin{array} { c } a _ { 1 } ^ { \alpha } t \frac { 1 } { 4 }, \\ [ a _ { 1 } ^ { \alpha } t _ { 1 } + a _ { 2 } ^ { \alpha } ( t - t _ { 1 } ) ] \, \frac { 1 } { 4 }, \\ [ a _ { 1 } ^ { \alpha } t + a _ { 2 } ^ { \alpha } ( t _ { 2 } - t _ { 1 } ) + a _ { 3 } ^ { \alpha } ( t - t _ { 2 } ) ] \, \frac { 1 } { 4 }, \\ [ a _ { 1 } ^ { \alpha } t _ { 1 } + a _ { 2 } ^ { \alpha } ( t _ { 2 } - t _ { 1 } ) + a _ { 3 } ^ { \alpha } ( t _ { 3 } - t _ { 2 } ) + a _ { 4 } ^ { \alpha } ( t - t _ { 3 } ) ] \, \frac { 1 } { 4 }, \ m = 4. \end{array} \\ \end{cases}$$

with m = 1 2 3 4 being the number of the state. This is , , , plotted in Fig. 12 for different times t . Note that when 0 &lt; t &lt; t 1 , the solution Eq. (23) reduces to Eq. (7) in the main text.

## IV. CONCLUSION

In this manuscript, we introduced a way to generate statistics of the positional distribution that is consistent with the CTRW model for a long observation time t . The main idea is that we use the modified L´ evy stable law L ∗ α ( ξ ) to simulate the number of renewals, which yields a perfect prediction for the positional distribution. In some sense, the method avoids the generating of IID random waiting times, which saves a significant amount of time. Mathematically, we can use the following formula to predict the positional distribution

$$\text{act one posterior asurrituation} & & \text{prov} \\ P ( x, t ) & \sim \int _ { - \frac { t } { ( \tau ) } ( \tilde { t } / t ) ^ { 1 / \alpha } } ^ { \infty } \mathcal { L } _ { \alpha } ^ { * } ( \xi ) & & \text{the n} \\ & \quad \times \frac { \exp \left ( - \frac { ( x - a \frac { t } { ( \tau ) } - a \xi ( \frac { t } { t } ) ^ { 1 / \alpha } ) ^ { 2 } } { 2 \sigma ^ { 2 } ( t / ( \tau ) + \xi ( t / t ) ^ { 1 / \alpha } ) } \right ) } { \sqrt { 2 \sigma ^ { 2 } \pi ( \frac { t } { ( \tau ) } + \xi ( \frac { t } { t } ) ^ { 1 / \alpha } ) } } d \xi, & & \text{case}$$

with

$$\mathcal { L } _ { \alpha } ^ { * } ( \xi ) = \begin{cases} \mathcal { L } _ { \alpha } ^ { * } ( \xi ), \ \xi \geq \xi ^ { * } ; \\ Q _ { t } ^ { I } ( \xi ), \ \bar { b } < \xi < \xi ^ { * } \quad & ( 2 8 ) \\ \bar { b }, \quad & \text{else.} \end{cases} \\ \text{method is not limited to examples discussed in the} \,. \,. \, \sim \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,. \,.$$

The method is not limited to examples discussed in the manuscript. On the one hand, it is valid for the case of waiting times having an infinite mean, where the number of renewals follows the one-sided L´vy stable distrie bution [27]. On the other hand, the displacements can be extended to situations when steps have a finite mean and a finite variance. When discussing the MSD, deviations exist between simulations and theoretical predictions. However, we did not find such deviations for the positional distribution and the mean of the position. This indicates that, for the biased CTRW model, the MSD is sensitive to the details of the positional distribution. Improving the presented work is a task left for the future.

Recently, researchers have been increasingly drawn to the nearly exponential decay observed in stochastic processes [11, 50-70]. Unlike the approach investigated here, the exponential decay is easy to detect at a short time scale, such as, the total observation time is half. Consequently, Eq. (27) fails. How to extend the present approach is left for future work.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00299v2-figure-6036f62d-86e6-4df4-83fe-c8139537547b.jpg)

t

FIG. 12: Breakthrough curves with time-dependent bias. The theory is based on the algorithm given in Sec. II, asymptotic behavior is given in Eq. (23), and simulations are obtained from 2 × 10 7 trajectories. Here we choose t a = t b / 2 = t c / 3 = 100, t = 1000, η = 4, and σ = 5.

- [1] M. F. Shlesinger, J. Stat. Phys. 10 , 421 (1974).
- [2] E. W. Montroll and G. H. Weiss, J. Math. Phys. 6 , 167 (1965).
- [3] J.-P. Bouchaud and A. Georges, Phys. Rep. 195 , 127 (1990).
- [4] R. Metzler and J. Klafter, Phys. Rep. 339 , 1 (2000).
- [5] G. Gradenigo, E. Bertin, and G. Biroli, Phys. Rev. E 93 , 060105(R) (2016).
- [6] R. Kutner and J. Masoliver, Eur. Phys. J. B 90 , 50 (2017).
- [7] J. Liu, P. Zhu, J.-D. Bao, and X. Chen, Phys. Rev. E 105 , 014122 (2022).
- [8] H. Scher and E. W. Montroll, Phys. Rev. B 12 , 2455 (1975).
- [9] A. Nissan, I. Dror, and B. Berkowitz, Water Resour. Res. 53 , 3760 (2017).
- [10] J. Masoliver, M. Montero, and G. H. Weiss, Phys. Rev. E 67 , 021112 (2003).
- [11] M. Hu, H. Chen, H. Wang, S. Burov, E. Barkai, and D. Wang, ACS Nano 17 , 21708 (2023).
- [12] E. Abdel-Rehim and R. Gorenflo, J. Comput. Appl. Math. 222 , 274 (2008).
- [13] M. Magdziarz, J. Stat. Phys. 135 , 763 (2009).
- [14] T. Akimoto and E. Yamamoto, Phys. Rev. E 93 , 062109 (2016).
- [15] S. N. Majumdar and H. Orland, J Stat Mech: Theory Exp 2015 , P06039 (2015).
- [16] H. C. Fogedby, Phys. Rev. E 50 , 1657 (1994).
- [17] D. Kleinhans and R. Friedrich, Phys. Rev. E 76 , 061102 (2007).
- [18] R. Co and C. Villagonzalo, In Proceedings of the 37th Samahang Pisika ng Pilipinas Physics Conference , SPP-

## Acknowledgments

WW expresses gratitude to Eli Barkai for the discussion that inspired this work. The work is supported by the National Natural Science Foundation of China under Grant No. 12105243 and 12101555, and the Zhejiang Province Natural Science Foundation LQ22A050002.

2019-PB-48 (Tagbilaran City, Philippines, 2019).

- [19] M. Magdziarz and M. Teuerle, Commun. Nonlinear Sci. Numer. Simul. 20 , 489 (2015).
- [20] G. Bel and E. Barkai, J. Phys.: Condens. Matter 17 , S4287 (2005).
- [21] R. Burioni, G. Gradenigo, A. Sarracino, A. Vezzani, and A. Vulpiani, Commun. Theor. Phys. 62 , 514 (2014).
- [22] R. Hou, A. G. Cherstvy, R. Metzler, and T. Akimoto, Phys. Chem. Chem. Phys. 20 , 20827 (2018).
- [23] T. Akimoto, A. G. Cherstvy, and R. Metzler, Phys. Rev. E 98 , 022105 (2018).
- [24] T. M. Michelitsch, F. Polito, and A. P. Riascos, Fractal Fract. 4 , 51 (2020).
- [25] S. Burov, W. Wang, and E. Barkai, Exponential tails and asymmetry relations for the spread of biased random walks. ArXiv:2209.03410.
- [26] M. Kotulski, J. Stat. Phys. 81 , 777 (1995).
- [27] C. Godr` eche and J. M. Luck, J. Stat. Phys. 104 , 489 (2001).
- [28] W. Wang, J. H. P. Schulz, W. H. Deng, and E. Barkai, Phys. Rev. E 98 , 042139 (2018).
- [29] W. Wang and E. Barkai, Phys. Rev. Lett. 125 , 240606 (2020).
- [30] W. Wang and E. Barkai, Journal of Physics A: Mathematical and Theoretical 57 , 035203 (2024).
- [31] J. Klafter and I. M. Sokolov, First Steps in Random Walks: From Tools to Applications (Oxford University Press, Oxford, 2011).
- [32] W. Deng, R. Hou, W. Wang, and P. Xu, Modeling Anomalous Diffusion: From Statistics to Mathematics (World Scientific, Singapore, 2020).
- [33] J. Klafter and G. Zumofen, J. Phys. Chem. 98 , 7366

(1994).

- [34] E. Barkai, Phys. Rev. E 63 , 046118 (2001).
- [35] A. I. Saichev and G. M. Zaslavsky, Chaos: An Interdisciplinary Journal of Nonlinear Science 7 , 753 (1997).
- [36] A. Chechkin and I. M. Sokolov, Phys. Rev. E 103 , 032133 (2021).
- [37] Z. Tong and A. Liu, Physica A 593 , 126941 (2022).
- [38] P. Semeraro, Mathematics and Financial Economics 16 , 685 (2022).
- [39] E. Aghion, D. A. Kessler, and E. Barkai, Phys. Rev. Lett. 118 , 260601 (2017).
- [40] J. M. Chambers, C. L. Mallows, and B. W. Stuck, J. Amer. Statist. Assoc. 71 , 340 (1976).
- [41] A. Janicki and A. Weron, Mathematics and Computers in Simulation 39 , 9 (1995).
- [42] M. Veillette, Matlab Central File Exchange, retreived October 10 , 2012 (2012).
- [43] W. Wang, A. Vezzani, R. Burioni, and E. Barkai, Phys. Rev. Res. 1 , 033172 (2019).
- [44] A. Vezzani, E. Barkai, and R. Burioni, Phys. Rev. E 100 , 012108 (2019).
- [45] R. Burioni, G. Gradenigo, A. Sarracino, A. Vezzani, and A. Vulpiani, J. Stat. Mech. Theory Exp. 2013 , P09022 (2013).
- [46] B. Berkowitz, A. Cortis, M. Dentz, and H. Scher, Rev. Geophys. 44 , RG2003 (2006).
- [47] B. Berkowitz and H. Scher, Adv. Water Resour. 32 , 750 (2009).
- [48] R. de Moraes, L. Ozelim, and A. Cavalcante, Sustainability 14 (2022).
- [49] T. J. Doerries, A. V. Chechkin, R. Schumer, and R. Metzler, Phys. Rev. E 105 , 014105 (2022).
- [50] W. K. Kegel and A. van Blaaderen, Science 287 , 290 (2000).
- [51] J. Masoliver, M. Montero, and J. M. Porr`, Physica A a 283 , 559 (2000).
- [52] E. R. Weeks, J. C. Crocker, A. C. Levitt, A. Schofield, and D. A. Weitz, Science 287 , 627 (2000).
- [53] P. Chaudhuri, L. Berthier, and W. Kob, Phys. Rev. Lett.
- 99 , 060604 (2007).
- [54] B. Wang, S. M. Anthony, S. C. Bae, and S. Granick, Proc. Natl. Acad. Sci. U.S.A. 106 , 15160 (2009).
- [55] S. Hapca, J. W. Crawford, and I. M. Young, J. R. Soc. Interface 6 , 111 (2009).
- [56] K. C. Leptos, J. S. Guasto, J. P. Gollub, A. I. Pesci, and R. E. Goldstein, Phys. Rev. Lett. 103 , 198103 (2009).
- [57] C. Eisenmann, C. Kim, J. Mattsson, and D. A. Weitz, Phys. Rev. Lett. 104 , 035502 (2010).
- [58] T. Toyota, D. A. Head, C. F. Schmidt, and D. Mizuno, Soft Matter 7 , 3234 (2011).
- [59] M. J. Skaug, J. Mabry, and D. K. Schwartz, Phys. Rev. Lett. 110 , 256101 (2013).
- [60] C. Xue, X. Zheng, K. Chen, Y. Tian, and G. Hu, J. Phys. Chem. 7 , 514 (2016).
- [61] D. Wang, H. Wu, and D. K. Schwartz, Phys. Rev. Lett. 119 , 268001 (2017).
- [62] R. Jeanneret, D. O. Pushkin, V. Kantsler, and M. Polin, Nat. Commun. 7 , 12518 (2016).
- [63] A. V. Chechkin, F. Seno, R. Metzler, and I. M. Sokolov, Phys. Rev. X 7 , 021002 (2017).
- [64] A. G. Cherstvy, S. Thapa, C. E. Wagner, and R. Metzler, Soft Matter 15 , 2526 (2019).
- [65] P. Witzel, M. G¨ otz, Y. Lanoisel´e, e T. Franosch, D. S. Grebenkov, and D. Heinrich, Biophys. J. 117 , 203 (2019).
- [66] K. Shin, S. Song, Y. H. Song, S. Hahn, J.-H. Kim, G. Lee, I.-C. Jeong, J. Sung, and K. T. Lee, J. Phys. Chem. 10 , 3071 (2019).
- [67] R. K. Singh, J. Mahato, A. Chowdhury, A. Sain, and A. Nandi, J. Chem. Phys. 152 , 024903 (2020).
- [68] C. Mejia-Monasterio, S. Nechaev, G. Oshanin, and O. Vasilyev, New J. Phys. 22 , 033024 (2020).
- [69] C. Xue, X. Shi, Y. Tian, X. Zheng, and G. Hu, Nano Letters 20 , 3895 (2020).
- [70] L. Defaveri, E. Barkai, and D. A. Kessler, Stretchedexponential relaxation in weakly-confined Brownian systems through large deviation theory (2023).
Danhua Jiang Yuanze Hong Wanli Wang
Simulation of the continuous-time random walk using subordination schemes
The continuous time random walk model has been widely applied in various
fields, including physics, biology, chemistry, finance, social phenomena, etc.
In this work, we present an algorithm that utilizes a subordinate formula to
generate data of the continuous time random walk in the long time limit. The
algorithm has been validated using commonly employed observables, such as
typical fluctuations of the positional distribution, rare fluctuations, the
mean and the variance of the position, and breakthrough curves with
time-dependent bias, demonstrating a perfect match.