Next Article in Journal
Assessment of Fractal Synchronization during an Epileptic Seizure
Next Article in Special Issue
Generalized Kinetic Equations with Fractional Time-Derivative and Nonlinear Diffusion: H-Theorem and Entropy
Previous Article in Journal
A Hierarchical Multi-Task Learning Framework for Semantic Annotation in Tabular Data
Previous Article in Special Issue
Fisher and Shannon Functionals for Hyperbolic Diffusion
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

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

by
Katarzyna Górska
1,
Francisco J. Sevilla
2,
Guillermo Chacón-Acosta
3 and
Trifce Sandev
4,5,6,*
1
Institute of Nuclear Physics, Polish Academy of Science, ul. Radzikowskiego 152, PL-31342 Kraków, Poland
2
Instituto de Física, Universidad Nacional Autónoma de México, Apdo. Postal 20-364, Ciudad de México 01000, Mexico
3
Departamento de Matemáticas Aplicadas y Sistemas, Universidad Autónoma Metropolitana-Cuajimalpa, Vasco de Quiroga 4871, Santa Fe, Cuajimalpa, Ciudad de México 05348, Mexico
4
Research Center for Computer Science and Information Technologies, Macedonian Academy of Sciences and Arts, Bul. Krste Misirkov 2, 1000 Skopje, Macedonia
5
Institute of Physics, Faculty of Natural Sciences and Mathematics, Ss. Cyril and Methodius University, Arhimedova 3, 1000 Skopje, Macedonia
6
Department of Physics, Korea University, Seoul 02841, Republic of Korea
*
Author to whom correspondence should be addressed.
Entropy 2024, 26(8), 665; https://doi.org/10.3390/e26080665
Submission received: 29 June 2024 / Revised: 29 July 2024 / Accepted: 30 July 2024 / Published: 5 August 2024
(This article belongs to the Special Issue Theory and Applications of Hyperbolic Diffusion and Shannon Entropy)

Abstract

:
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.

1. Introduction

Due to economic and military reasons in the XIX-th century, the rapid communication between the countries belonging to the British Empire was an important matter. To send a message from one part of the Empire to another was a true challenge, for its solution contributed to many famous physicists of that time, namely J. Maxwell, Lord Kelvin, and O. Heaviside. Nevertheless, only Heaviside solved this problem. He proposed the so-called telegrapher’s (telegraph) equation (TE) to describe the current propagation inside the telegrapher’s cable [1]. His equation is a hyperbolic equation which, in ( 1 + 1 ) -dimensional space, has the following form
t 2 p 0 ( x , t ) + τ 1 t p 0 ( x , t ) = v 2 x 2 p 0 ( x , t ) , x R and t R + ,
where τ is a time parameter measured in seconds, and v is the propagation velocity measured in seconds/meter. The boundary conditions are set to zero at infinity, while the initial conditions are
p 0 ( x , t = 0 ) = δ ( x x 0 ) , t p 0 ( x , t ) | t = 0 = 0 .
These boundary and initial conditions, allowed C. R. Cattaneo [2,3] and P. Vernotte [4] to employ Equation (1) to describe the heat transport. Thus, the physical interpretation of τ and v is changing. Here, τ modifies the Fourier law which connects the heat current j ( x , t ) and the temperature T ( x , t ) such that we have
j ( x , t ) + τ 1 t j ( x , t ) x T ( x , t ) .
The parameter v is associated with heat conductivity. Recently, Equation (1) has been used to characterize the diffusion process with finite propagation velocity v [5], where v = K / τ , K is the diffusion coefficient, and the Fourier law is replaced by Fick’s law. In this case, the solution of the telegrapher’s equation is a probability density distribution (PDF) and, for arbitrary initial conditions, is presented in Ref. [6] (Equation (102) on p. 303 and/or [7] Equation (7.4.28)). This solution for the initial conditions given by Equation (2) reduces to the form
p 0 ( x , t ) = 1 2 e t / ( 2 τ ) δ v t | x x 0 | + 1 4 v τ e t / ( 2 τ ) Θ v t | x x 0 | × I 0 v 2 t 2 | x x 0 | 2 2 v τ + v t v 2 t 2 | x x 0 | 2 I 1 v 2 t 2 | x x 0 | 2 2 v τ ,
derived in Ref. [8]. We use the standard notation for which Θ ( · ) denotes the Heaviside step function, and I ν ( z ) is the modified Bessel function of the first kind of order ν = 0 , 1 . The corresponding mean squared displacement (MSD) reads
x 2 ( t ) 0 = R x 2 p 0 ( x , t ) d x = 2 K τ t τ + e t / τ 1 ,
which leads to the ballistic time dependence x 2 ( t ) 0 v 2 t 2 in the short-time regime and to the linear time dependence x 2 ( t ) 0 2 K t in the long-time one.
Equation (1) characterizes the telegraph process whose corresponding Langevin equation reads
x ˙ ( t ) = v ζ ( t ) ,
where v denotes the constant particle speed, and ζ ( t ) is a stationary dichotomic Markov process that jumps between two states; ± 1 , with a mean rate ν [9,10,11,12] ( τ = 1 2 ν , is a time scale that corresponds to the inverse mean sojourn time for each state). As shown in [13,14] subordinating this process by the Lévy process, we can find the fractional telegrapher’s equation of type I (abbreviated as FTE-I), namely
τ μ D t 2 μ C p 1 ( x , t ) + D t μ C p 1 ( x , t ) = K x 2 p 1 ( x , t ) ,
where μ ( 0 , 1 ] and D t μ C h ( t ) is the fractional derivative in the Caputo sense (see Appendix A). Equation (5) can be interpreted as an anomalous diffusion equation, which can be derived by using the continuity equation with an appropriate modified constitutive relation [15].
Despite the compelling generalization given by Equation (5), whose higher fractional derivative is always twice the lower one, it is possible to postulate different kinds of fractional telegrapher’s equations. For instance, we focus on one of these possibilities called type II (abbreviated as FTE-II), which has the form (this form corresponds to Equation (17) in Ref. [15])
τ D t μ + 1 C p 2 ( x , t ) + D t μ C p 2 ( x , t ) = K x 2 p 2 ( x , t ) ,
where μ ( 0 , 1 ] . In this case, the higher fractional derivative is always larger than 1. Equations (5) and (6) are completed with the same initial and boundary conditions (2) as for the telegrapher’s Equation (1). We can also find different fractional generalisations of the TE [5,10,15,16,17,18], including non-Markovian discrete time versions of the telegraph process [19], TE in random media [20,21,22], telegraph processes with random velocities [23], etc.
In this paper, we consider the stochastic resetting of the PDFs p 1 ( x , t ) and p 2 ( x , t ) that solve Equations (5) and (6), respectively. The stochastic dynamics of the particle, which are initially located at x = x 0 at t = 0 , under the effects of stochastic resetting, are described by the renewal equation [24]
p r ( x , t ) = ϱ ( t ) p j ( x , t ) + 0 t ρ ( t ) p r ( x , t t ) d t ,
where j = 0 , 1 , 2 , indicates the probability density corresponding to the different telegrapher’s equations considered here; ϱ ( t ) = t ρ ( t ) d t is the probability that no renewal has taken place up to time t. We consider here the stochastic Poissonian resetting with probability ρ ( t ) = r e r t and from this, we have ϱ ( t ) = e r t , where r is the resetting rate. The renewal Equation (7) for Poissonian resetting has the form [24,25,26]:
p r , j ( x , t ) = e r t p j ( x , t ) + 0 t r e r t p j ( x , t ) d t ,
which, in the Laplace space, reads
p ^ r , j ( x , s ) = s + r s p ^ j ( x , s + r ) , j = 0 , 1 , 2 .
Furthermore, we assume that the particle is instantaneously reset to the initial position. These assumptions offer the possibility of analytical calculations, as well a description of the experimental results for the mean-first passage time under stochastic resetting, by using holographic optical tweezers [27] or laser traps [28].
The paper is organized as follows. In Section 2, we consider the stochastic resetting of FTE-I. For that purpose, we find its solution in two ways based on integral decomposition techniques. This decomposition is rooted in recognizing the Brownian motion or the telegraph process. Next, we will reset the process described by FTE-I and calculate the first-passage time problem and survival probability. In Section 3, we repeat all procedures with the stochastic resetting for FTE-II. The paper is summarized in Section 4. The paper contains four Appendices.

2. FTE-I under Resetting

The PDF p 1 ( x , t ) can be found in two ways, which come from the Efros theorem (Appendix C) applied to p 1 ( x , t ) . In the Laplace–Fourier (LF) space, it is written as
p ^ ˜ 1 ( κ , s ) = s 1 ( τ μ s 2 μ + s μ ) τ μ s 2 μ + s μ + K κ 2 e i κ x 0 = [ M ^ 1 ( s ) ] 1 s [ M ^ 1 ( s ) ] 1 + K κ 2 e i κ x 0 , μ ( 0 , 1 ) .
Here,
M ^ 1 ( s ) = τ μ s 1 μ s μ + τ μ ,
and s ÷ t and κ ÷ x are the Laplace and Fourier coordinates, respectively. The Efros theorem allows us to make the integral decomposition of Equation (10), for which we represent p 1 ( x , t ) as the integral 0 N ( x , ξ ) f 1 ( ξ , t ) d ξ with N ( x , ξ ) the Gaussian
N ( x , ξ ) = 1 2 π K ξ exp ( x x 0 ) 2 4 K ξ
and
f 1 ( ξ , t ) = L 1 [ M ^ 1 ( s ) ] 1 e ξ s / M ^ 1 ( s ) ; t .
The function f 1 ( ξ , t ) for μ ( 0 , 1 / 2 ] is non-negative, such that it can be named the PDF of the leading process. Another possibility is to present p 1 ( x , t ) as the solution of the telegrapher’s equation p 0 ( x , ξ ) given by Equation (3) and h ( ξ , t ) = L 1 [ ( τ s ) μ 1 exp ( ξ τ μ 1 s μ ) ; t ] being for μ ( 0 , 1 ] the PDF of the leading process (for details, see Section 2.2). Hence, the solution obtained in this way is called the subordination approach. Both techniques in the presence of stochastic resetting will be described below, and we show that they lead to the same results.

2.1. Resetting of FTE-I—The First Possibility

Let us begin by examining the outcomes that can be derived from the second equality in Equation (10). In the short-time regime, s 1 / τ , we have that Equation (10) is approximated by
p ^ ˜ 1 f W E ( κ , s ) τ μ s 2 μ 1 τ μ s 2 μ + K κ 2 e i κ x 0 , μ ( 0 , 1 ) ,
which, after inverting the Laplace and Fourier transforms, we obtain the fractional wave equation (fWE)
τ μ D t 2 μ C p 1 f W E ( x , t ) = K x 2 p 1 f W E ( x , t ) .
Transport in this regime is described by the fWE. In the opposite regime, i.e., in the long-time regime, s 1 / τ , Equation (10) is asymptotically approximated by
p ^ ˜ 1 f D E ( κ , s ) s μ 1 s μ + K κ 2 e i κ x 0 , μ ( 0 , 1 ) ,
which corresponds to the standard fractional diffusion equation (fDE) after inverting the Laplace and Fourier transforms, i.e.,
D t μ C p 1 f D E ( x , t ) = K x 2 p 1 f D E ( x , t ) .
In the general case, upon performing the Fourier transform inversion of Equation (10), we can express it in Laplace space as:
s p ^ 1 ( x , s ) δ ( x x 0 ) = K M ^ 1 ( s ) x 2 p ^ 1 ( x , s ) ,
which, after taking the inverse Laplace transform, can be written in t-space as the generalized diffusion equation
t p 1 ( x , t ) = K 0 t M 1 ( t t ) x 2 p 1 ( x , t ) d t , M 1 ( t ) = τ μ t 2 μ 2 E μ , 2 μ 1 ( t / τ ) μ .
The series representation of the two-parameter Mittag–Leffler function E μ , ν ( z ) is given by Equation (A1) for δ = 1 . After multiplying the last equation by x 2 and integrating over space we are led to first-order differential equations for the MSD, which can be integrated straightforwardly to give x 2 ( t ) 1 = 2 K 0 t 0 t M 1 ( s ) d s d t . We emphasize that the procedure of expressing the hyperbolic-like Equation (5) as the parabolic-like Equation (19) is only formal and it is possible only for the diffusion-like initial conditions (2). Formally, we can present the solution of Equation (19) as
p 1 ( x , t ) = 0 N ( x , u ) f 1 ( u , t ) d u ,
where N ( x , u ) and f 1 ( u , t ) are given by Equations (12) and (13), respectively. The same results can be obtained using Equation (A3) resulting from the Efros theorem in which G ^ ( s ) = [ M ^ 1 ( s ) ] 1 , q ^ ( s ) = s / M ^ 1 ( s ) , and g ^ ( x , s ) = N ^ ( x , s ) . In the LF space, it figures as
p ^ 1 ( x , s ) = [ M ^ 1 ( s ) ] 1 N ^ x , s / M ^ 1 ( s ) .
Note that f 1 ( u , t ) is a non-negative function for μ ( 0 , 1 / 2 ] . That comes from the Bernstein theorem [29] and the fact that f ^ 1 ( u , s ) , for this range of μ , is a completely monotonic function (CMF), i.e., the non-negative function whose derivatives exist and alternate, see Appendix B. In consequence, we can say that p 1 ( x , t ) , given by Equation (20) for μ ( 0 , 1 / 2 ] is a PDF expressed by the subordination approach, in which the leading process, described by f 1 ( u , t ) , subordinates the parent process characterized by the normal distribution N ( x , u ) [13,14]. The function f ^ 1 ( u , s ) for μ ( 1 / 2 , 1 ] is not a CMF. Then, according to the Bernstein theorem, f 1 ( u , t ) is negative or contains negative parts. However, p ^ 1 ( x , t ) calculated from Equation (21) for μ ( 0 , 1 ] is a CMF, yielding that its inverse is non-negative and can be called a PDF.
The corresponding MSD x 2 ( t ) 1 can be expressed in terms of the MSD of the normal distribution x 2 ( t ) N = 2 K t , we have
x ^ 2 ( s ) 1 = [ M ^ 1 ( s ) ] 1 x ^ 2 ( s / M ^ 1 ( s ) ) N = 2 K τ μ s μ 1 s μ + τ μ ,
from where, by inverse Laplace transform, one finds
x 2 ( t ) 1 = 2 K τ μ ( t / τ ) 2 μ E μ , 2 μ + 1 ( t / τ ) μ ,
which is a known result [5]. From (22), we have that x 2 ( t ) 1 ( 2 K / τ μ ) t 2 μ / Γ ( 2 μ + 1 ) in the short-time regime, and x 2 ( t ) 1 2 K t μ / Γ ( μ + 1 ) in the long-time regime; see Figure 1, curves marked with circle symbols for μ = 1 / 4 , 1/2, and 3/4.
If we reset p 1 ( x , t ) to the initial position according to the stochastic Poissonian resetting ρ ( λ ) = r e r λ , where r is the resetting rate, then the corresponding p r , 1 ( x , t | x 0 ) p r , 1 ( x , t ) follows the renewal Equation (8) for j = 1 . In Laplace space, it reads Equation (9) for j = 1 , which, after employing Equation (21), gives
p ^ r , 1 ( x , s ) = s + r s 1 M ^ 1 ( s + r ) N ^ x , s + r M ^ 1 ( s + r ) .
From here, we conclude that in the long-time limit, the system approaches the non-equilibrium stationary state (NESS) p r , 1 st ( x ) = lim t p r , 1 ( x , t ) = lim s 0 s p ^ r , 1 ( x , s ) = r p ^ 1 ( x , r ) , which is explicitly given by
p r , 1 st ( x ) = r M ^ 1 ( r ) N ^ x , r M ^ 1 ( r ) = 1 2 r μ ( τ μ r μ + 1 ) K exp r μ ( τ μ r μ + 1 ) K | x x 0 | .
The NESS is given by the Laplace distribution with a stationary variance σ 1 st = K / [ r μ ( τ μ r μ + 1 ) ] around the average value x 0 . For μ = 1 , the NESS is analogous to Equation (27) of Ref. [8]. For τ = 0 , we recover the result for the NESS for an anomalous diffusion process under resetting [24,30,31,32,33], i.e., p r , 1 st ( x ) = r μ / 2 / ( 2 K ) exp r μ / K | x x 0 | . We want to point out the exponential decay of the NESS tails, which is a characteristic induced by the stochastic Poissonian resetting process.
For the MSD in the case of resetting, we have
x ^ 2 ( s ) r , 1 = s + r s x ^ 2 ( s + r ) 1 = s + r s 1 M ^ 1 ( s + r ) x ^ 2 s + r M ^ 1 ( s + r ) N = 2 K τ μ s 1 ( s + r ) μ ( s + r ) μ + τ μ ,
from which we obtain
x 2 ( t ) r , 1 = 2 K τ μ 1 0 t e r t ( t / τ ) 2 μ 1 E μ , 2 μ ( t / τ ) μ d t .
In the short-time limit it behaves as x 2 ( t ) r , 1 t 2 μ , and in the long-time limit it approaches the constant value,
x 2 ( t ) r , 1 = 2 K r μ ( r τ ) μ + 1 .
The graphical representation of the MSD (25) for different values of μ and r is given in Figure 1, where we observe the asymptotic behavior of the MSD for short- and long-time limits.

2.2. Resetting of FTE-I—The Second Possibility

As mentioned in the Introduction, the FTE-I for μ ( 0 , 1 ] can be obtained from the telegraph process subordinated by Lévy noise. Here, we demonstrate this approach, starting from Langevin’s equations:
x ˙ ( u ) = v ζ ( u ) and t ˙ ( u ) = ξ ( u ) ,
where ζ ( u ) represents the same dichotomic noise as in the standard TE, and ξ ( u ) is a Lévy stable noise with Lévy index in Laplace space given by Ψ ^ ( s ) = τ μ 1 s μ , μ ( 0 , 1 ) . Therefore, the process t ( u ) = 0 u ξ ( u ) d u is a stable Lévy motion with the characteristic function given by the stretched exponential function Φ ^ μ ( u , s ) = exp ( u τ μ 1 s μ ) . Its inverse Laplace transform is denoted as Φ μ ( u , t ) and known as a one-sided Lévy stable distribution whose series form was found by H. Pollard in [34]. Its representation through the Fox H function was later found by R. Hilfer [35]. In [36], it was presented in the Meijer G form and the finite series of generalized hypergeometric functions. The corresponding PDF p 1 ( x , t ) of this subordinated telegraph process can be found from the subordination integral [13,14,37,38,39]
p 1 ( x , t ) = 0 p 0 ( x , u ) h ( u , t ) d u ,
where h ( u , t ) = u Θ ( t t ( u ) ) is the subordination function which, in Laplace space, reads as
h ^ ( u , s ) = 1 s u 0 δ ( t t ( u ) ) e s t d t = s 1 u e s t ( u ) = s 1 u Φ ^ μ ( u , s ) = ( τ s ) μ 1 e u τ μ 1 s μ .
Inverting (28), we get the explicit time dependence of the subordination function
h ( u , t ) = L 1 [ h ^ ( u , s ) ; t ] = t μ u Φ μ ( u , t ) .
This subordination is also obtained from the first equality of Equation (10) by employing the Efros theorem in which G ^ ( s ) = ( τ s ) μ 1 , q ^ ( s ) = τ μ 1 s μ , and g ^ ( x , s ) = p ^ 0 ( x , s ) [13,14]. Then, from Equation (A3), we have
p ^ 1 ( x , s ) = ( τ s ) μ 1 p ^ 0 ( x , τ μ 1 s μ ) ,
whose inverse Laplace transform gives Equation (27). From Equation (29), it follows that the MSD reads
x ^ 2 ( s ) 1 = ( τ s ) μ 1 x ^ 2 ( τ μ 1 s μ ) 0 = 2 K τ μ s μ 1 s μ + τ μ ,
where x 2 ( t ) is the MSD (4) for the standard TE. Therefore, we obtain the same MSD as in Equation (22), as expected. Thus, making the Poissonian resetting as in the previous subsection, we obtain Equation (25); in the long-time limit, this gives Equation (26). Moreover, NESS, given by Equation (23) for this kind of subordination, is also equal to Equation (24).
For μ = 1 , we recover the known results for the standard telegrapher’s equation with stochastic resetting [8].

2.3. First-Passage Time Problem

From Equations (18) and (9), we obtain
s p ^ r , 1 ( x , s ) δ ( x x 0 ) = s s + r M ^ 1 ( s + r ) x 2 p ^ r , 1 ( x , s ) ,
where, by inverse Laplace transform, the generalized diffusion equation
t p r , 1 ( x , t ) = K 0 t M 1 ( t t ) r M ¯ 1 ( t t ) x 2 p r , 1 ( x , t ) d t ,
is obtained, with memory function M 1 ( t ) r M ¯ 1 ( t ) , where
M ¯ 1 ( t ) = e r t 0 t M 1 ( t ) d t .
Let us write the corresponding backward equation for the survival probability Q 1 ( x 0 , t ) , which will give the probability of the particle starting at x 0 > 0 to reach the target at the origin. Thus, we have
D t 2 μ C Q 1 ( x 0 , t ) + τ μ D t μ C Q 1 ( x 0 , t ) = v 2 x 0 2 Q 1 ( x 0 , t ) ,
with initial conditions
Q 1 ( x 0 , 0 ) = 1 , t Q 1 ( x 0 , t ) | t = 0 = 0 ,
and boundary conditions Q 1 ( 0 , t ) = 0 and Q 1 ( , t ) = 0 . By Laplace transform, we find that
Q ^ 1 ( x 0 , s ) = 1 s 1 exp τ μ K s μ + τ μ s μ x 0 .
From here, we can calculate the first-passage time density P 1 ( t ) = d d t Q 1 ( x 0 , t ) , i.e., P ^ 1 ( s ) = 1 s Q ^ 1 ( x 0 , s ) , see [40,41,42], which, in the Laplace space, reads
P ^ 1 ( s ) = exp τ μ K s μ + τ μ s μ x 0 .
If we consider the process as a random search, then we can calculate the efficiency of the search, defined as the number of visited targets divided by the average number of steps needed. If there is a single target, it can be calculated as the inverse of the MFPT [43]
E 1 = 1 t = 0 P 1 ( t ) t d t = 0 P ^ 1 ( s ) d s = 1 2 1 / μ μ Γ 1 + 1 2 μ 2 K x 0 1 / μ exp x 0 2 K τ μ ,
whose derivation is presented in Appendix D.
If we further consider exponential resetting to the telegraph process, then for the survival probability, one finds [41,42]
Q ^ 1 , r ( x 0 , s ) = Q ^ 1 ( x 0 , s + r ) 1 r Q ^ 1 ( x 0 , s + r ) ,
from where we derive the MFPT
T 1 , r ( x 0 ) = 0 t [ t Q 1 , r ( x 0 , t ) ] d t = Q 1 , r ( x 0 , s = 0 ) = 1 r exp r μ ( τ μ r μ + 1 ) K x 0 1 .
We see that the MFPT in the limits r 0 and r diverges, so there is an optimal resetting rate r * for which MFPT is minimal, i.e.,
r T 1 , r ( x 0 ) | r = r * = 0 ,
from where we have
μ ( 1 + 2 r * μ τ μ ) 2 ( r * μ τ μ + 1 ) r * μ ( τ μ r * μ + 1 ) K x 0 = 1 exp r * μ ( τ μ r * μ + 1 ) K x 0 ,
i.e.,
1 e ξ 1 = r * d ξ 1 d r * , ξ 1 = r * μ ( τ μ r * μ + 1 ) K x 0 .
For τ = 0 it reduces to [44]
μ 2 r * μ K x 0 = 1 exp r * μ K x 0 .
For μ = 1 , we recover the MFPT for the telegraph process under resetting [8], see also [45],
T r ( x 0 ) = 1 r exp r ( τ r + 1 ) K x 0 1 .
while for τ = 0 , the result for the subdiffusive search [44,46]
T r ( x 0 ) = 1 r exp r μ K x 0 1 ,
and for τ = 0 and μ = 1 , the Brownian search with [47]
T r ( x 0 ) = 1 r exp r K x 0 1 .
The dependence of the MFPT on the resetting rate r is shown in Figure 2, while the optimal resetting rate versus parameter μ is shown in Figure 3. From Equation (32), for a given exponent μ , we can numerically find the optimal value r for which the MFPT is minimum, see Figure 3. From Figure 2, we see that by increasing μ , the resetting rate r * , for which the MFPT is minimum, decreases. Therefore, we need a higher resetting rate in order to reset the particle, which is stacked due to the long-tailed waiting time, given by the fractional exponent μ (the lower fractional exponent μ means a longer waiting time, and therefore, we need a higher resetting rate for the particle to reach the target in a shorter time).

3. FTE-II under Resetting

By the same methods used in the previous section, we have that the Laplace transform of Equation (6) looks analogous to Equation (12), but with the difference that instead of M ^ 1 ( s ) , we have to use M ^ 2 ( s ) defined as
M ^ 2 ( s ) = 1 τ s 1 μ s + τ 1 and M 2 ( t ) = L 1 [ M ^ 2 ( s ) ; t ] = t μ 1 τ E 1 , μ t τ .
Proceeding similarly as in the previous case of FTE-I to obtain Equation (18), we have that
s p ^ 2 ( x , s ) δ ( x x 0 ) = K M ^ 2 ( s ) x 2 p ^ 2 ( x , s ) .
In the short-time regime, we have that the memory function M ^ 2 ( s ) 1 / ( s μ τ ) ; therefore, Laplace inversion of Equation (33) is possible since μ ( 0 , 1 ] . Thus, we have the fWE
τ D t μ + 1 C p 2 f W E ( x , t ) = K x 2 p 2 f W E ( x , t ) .
In the long-time regime M ^ 2 ( s ) s 1 μ we recover the fDE
D t μ C p 2 f D E ( x , t ) = K x 2 p 2 f D E ( x , t ) .
We point out that M ^ 2 ( s ) is not a completely Bernstein function (CBF) since its algebraic inverse is not a Stieltjes function (SF). Formally, we can take the inverse Laplace transform of Equation (33) with initial condition (2), but, due to the results given in Ref. [48], it will not be a well-posed Cauchy problem. Hence, we limit our consideration only to the Laplace space in which the solution of Equation (33) can be written as
p ^ 2 ( x , s ) = [ M ^ 2 ( s ) ] 1 N ^ x , s / M ^ 2 ( s ) .
The MSD corresponding to p 2 ( x , t ) becomes
x ^ 2 ( s ) 2 = 1 M ^ 2 ( s ) x ^ 2 s / M ^ 2 ( s ) N = 2 K τ s μ 1 s + τ 1 ,
from which, by inverse Laplace transform, one finds [5] (Equation (24)), namely
x 2 ( t ) 2 = 2 K τ μ t / τ 1 + μ E 1 , 2 + μ t / τ .
In the case of stochastic resetting, we express the renewal equation given by Equation (8) for j = 2 in Laplace space as
p ^ r , 2 ( x , s | x 0 ) = s + r s p ^ 2 ( x , s + r ) = s + r s 1 M ^ 2 ( s + r ) N ^ x , s + r M ^ 2 ( s + r ) ,
which, in the long-time limit, approaches the NESS equal to
p r , 2 st ( x ) = r M ^ 2 ( r ) N ^ x , r M ^ 2 ( r ) = 1 2 r μ ( r τ + 1 ) K exp r μ ( r τ + 1 ) K | x x 0 | .
Therefore, the NESS is given by the Laplace distribution, as well, with a stationary variance σ 2 st = K / [ r μ ( r τ + 1 ) ] around the average value x 0 . For μ = 1 , NESS reduces to the one for the standard telegrapher’s process under resetting [8]. For τ = 0 , it turns to the NESS for an anomalous diffusion process under resetting [24,30,31,32,33].
The corresponding MSD in the case of resetting becomes
x ^ 2 ( s ) r , 2 = 2 K τ s 1 ( s + r ) μ ( s + r ) + τ 1 ,
that is
x 2 ( t ) r , 2 = 2 K τ μ 1 0 t e r t t τ μ E 1 , 1 + μ t τ d t .
From here, for the short-time limit, we have the behavior x 2 ( t ) r , 2 t 1 + μ , while for the long-time limit we observe saturation
x 2 ( t ) r , 2 = 2 K r μ r τ + 1 .
This crossover dynamics of the MSD (38) for different values of μ and r is shown in Figure 4.

First-Passage Time Problem

We can first solve the backward equation for the survival probability without resetting and then directly find the survival probability with resetting. It reads
τ D t 1 + μ C Q 2 ( x 0 , t ) + D t μ C Q 2 ( x 0 , t ) = K x 2 Q 2 ( x 0 , t ) ,
with the same initial and boundary conditions as before, we find the survival probability
Q ^ 2 ( x 0 , s ) = 1 s 1 exp r μ ( τ r + 1 ) K x 0 ,
from where for the first-passage time density we have
P ^ 2 ( s ) = 1 s Q ^ 2 ( x 0 , s ) = exp s μ ( τ s + 1 ) K x 0 .
The efficiency then becomes
E 2 = 0 exp s μ ( τ s + 1 ) K x 0 d s = 1 μ τ π 4 K x 0 2 2 H 2 , 1 1 , 2 4 K τ μ x 0 2 1 μ 3 2 , 1 μ , 2 , 1 + 1 μ 1 , 1 ,
where H p , q m , n is the Fox H function in which the argument is given by [ ( 4 K τ μ ) / x 0 2 ] 1 / μ , and the upper and lower list of parameters reads [ 3 / 2 , 1 / μ ] , [ 2 , 1 + 1 / μ ] , and [ 1 , 1 ] , respectively. The derivation of the lower formula in Equation (39) is presented in Appendix D.
In the presence of resetting, we can directly find the survival probability from the corresponding one without resetting, from where for the MFPT we obtain
T 2 , r ( x 0 ) = 1 r exp r μ ( τ r + 1 ) K x 0 1 .
For μ = 1 , we obtain the known result for the MFPT in the case of a standard telegraph process. We can find the optimal resetting rate r for which MFPT is minimal, i.e.,
r T 2 , r ( x 0 ) | r = r * = 0 ,
from where we have
μ + r * τ + μ r * τ 2 ( r * τ + 1 ) r * μ ( τ r * + 1 ) K x 0 = 1 exp r * μ ( τ r * + 1 ) K x 0 ,
i.e.,
1 e ξ 2 = r d ξ 2 d r , ξ 2 = r * μ ( τ r * + 1 ) K x 0 .
The dependence of the MFPT on the resetting rate is shown in Figure 5, while the changing of the optimal resetting rate by changing parameter μ is demonstrated in Figure 6. We see that by increasing parameter μ , the optimal resetting rate decreases, the behavior which was also observed for the FTE-I.

4. Summary

In this paper, we considered two distinct kinds of fractional telegrapher equations. In the absence of resetting, we found that both generalizations of the telegrapher equations describe fractional-ballistic transport in the short-time regime, and transit to fractional diffusion transport in the long-time one. FTE-I can be obtained in two ways based on the integral decomposition method, both ways lead to the same results. The methods consist of presenting the solution for the FTE-I from either the normal distribution or from the PDF of the telegraph process. The solution of FTE-II is presented only by use of the normal distribution in the integral decomposition.
Further, we analyzed these two fractional telegraphic processes in the presence of Poissonian resetting, which means that after a random time drawn from an exponential distribution, the particle is reset to the initial position. We found that in the long-time limit, due to the resetting in both cases, the particle reaches non-equilibrium stationary states while the MSDs saturate. In both cases, if the resetting rate is large enough, the nonequilibrium stationary distribution is determined by the short-time regime of the solution of the corresponding FTE, while if the resetting rate is small, the stationary distribution is determined by the long-time regime of the corresponding solution of the FTE.
We also considered the first-passage time problem, and for both cases, we calculated the survival probability, efficiency, and MFPT. It is shown that there is an optimal resetting rate for which the MFPT is minimal. This optimal resetting rate depends on the anomalous diffusion parameter μ . Additionally, the efficiency in the first case decreases faster with μ than in the second case.
Analysis of the generalized telegraph processes under non-instantaneous resetting [49,50,51], partial resetting [52,53], and resetting in an interval [54,55,56] are left for future investigation. The first-passage time problem could also be of interest in the case of resetting to multiple [57] and random positions [58]. Finding the connection between the Shannon and Fisher functionals [59] could also be of interest to future research.

Author Contributions

Conceptualization, K.G., F.J.S., G.C.-A. and T.S.; methodology, K.G., F.J.S., G.C.-A. and T.S.; software, K.G., F.J.S., G.C.-A. and T.S.; validation, K.G., F.J.S., G.C.-A. and T.S.; formal analysis, K.G., F.J.S., G.C.-A. and T.S.; investigation, K.G., F.J.S., G.C.-A. and T.S.; writing—original draft preparation, K.G., F.J.S., G.C.-A. and T.S.; writing—review and editing, K.G., F.J.S., G.C.-A. and T.S.; visualization, K.G., F.J.S., G.C.-A. and T.S.; supervision, K.G., F.J.S., G.C.-A. and T.S. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by UNAM-PAPIIT IN112623. TS is supported by the German Science Foundation (DFG, Grant number ME 1535/12-1) and by the Alliance of International Science Organizations (Project No. ANSO-CR-PP-2022-05). TS is also supported by the Alexander von Humboldt Foundation. K.G. acknowledges the NCN Research Grant Preludium Bis 2 No. UMO-2020/39/O/ST2/01563.

Data Availability Statement

No new data were created or analyzed in this study. Data sharing is not applicable to this article.

Acknowledgments

K.G. and T.S. acknowledge the hospitality and support from the Lorentz Center in Leiden.

Conflicts of Interest

The authors declare no conflicts of interest.

Abbreviations

The following abbreviations are used in this manuscript:
TETelegrapher’s Equation
FTEFractional Telegrapher’s Equation
MSDMean Squared Displacement
MFPTMean First-Passage Time

Appendix A. Fractional Calculus and Mittag–Leffler Functions

The fractional derivative in the Caputo sense reads
D t ν C f ( t ) = 0 t f ( n ) ( ξ ) ( t ξ ) ν + 1 n d ξ Γ ( n ν ) , ν ( n 1 , n ) , n N .
The three-parameter Mittag–Leffler function is defined by [60]
E ρ , β δ ( z ) = n = 0 ( δ ) n Γ ( ρ n + β ) z n n ! ,
Its Laplace transform reads
L t β 1 E ρ , β δ λ t ρ = s ρ δ β ( s ρ + λ ) δ ,
where | λ / s ρ | < 1 . Note that for δ = 1 , it becomes the two parameter Mittag–Leffler function, E ρ , β 1 ( z ) = E ρ , β ( z ) , and for β = δ = 1 , it becomes the one parameter Mittag–Leffler function, E ρ , 1 1 ( z ) = E ρ , 1 ( z ) = E ρ ( z ) .
The Fox H function is defined through the Mellin transform as follows
H p , q m , n z | [ a p , A p ] [ b q , B q ] : = γ L i = 1 m Γ ( b i + B i s ) i = 1 n Γ ( 1 a i A i s ) i = n + 1 p Γ ( a i + A i s ) i = m + 1 q Γ ( 1 b i B i s ) x s d s 2 π i ,
where the parameters are subject to conditions
z 0 , 0 m q , 0 n p ; a i C , A i > 0 , i = 1 , , p ; b i C , B i > 0 , i = 1 , , q ; [ a p , A p ] = ( a 1 , A 1 ) , , ( a p , A p ) ; [ b q , B q ] = ( b 1 , B 1 ) , , ( b q , B q ) ; ( a p ) = a 1 , a 2 , , a p ; ( b q ) = b 1 , b 2 , , b q .

Appendix B. Completely Monotone and Bernstein Functions

Among the non-negative functions, we can distinguish the complete monotone and Bernstein functions. Based on [29], their definitions and subclasses are given below.
A function c ^ : ( 0 , ) R is a completely monotone function (CMF) if its derivatives exist and ( 1 ) n c ^ ( n ) ( s ) 0 , for all n N 0 .
The important feature of CMF is that the Bernstein theorem uniquely connects it with a non-negative function c ( t ) through the Laplace integral
c ^ ( s ) = 0 e s t c ( t ) d t .
A (non-negative) Stieltjes function (SF) is a function g ^ : ( 0 , ) [ 0 , ) which can be written as
g ^ ( s ) = A s + B + 0 σ ( d t ) s + t ,
where A , B 0 and σ is a measure on ( 0 , ) such that 0 σ ( d t ) / ( 1 + t ) < . Note that the class of SFs belongs to the class of CMFs, and thus, it is a decreasing function. Usually, it is assumed that A = 0 and B = 0 , which simplifies the calculations and ensures its integrability at zero and vanishing at infinity.
A function b ^ : ( 0 , ) R is a Bernstein function (BF) if its all derivatives exist, b ^ ( s ) is non-negative, b ^ ( s ) is a CMF. If a BF b ^ : ( 0 , ) ( 0 , ) , and, additionally, b ^ ( s ) / s is an SF, then it is called a completely Bernstein function (CBF).
The characteristics of these functions and their properties can be referenced in [29]. The key features highlighted in this paper include (i) the product or sum of CMFs results in another CMF; (ii) the composition of CMF and BF(CBF) yields a CMF; (iii) for the CBF, there exists a corresponding partner b ^ such that b ^ ( s ) b ^ ( s ) = s ; and (iv) the algebraic inverse of SF (CBF) is CBF (SF).

Appendix C. Efros Theorem

The Efros theorem [61,62,63,64,65] generalizes the convolution theorem for the Laplace transform. It states the following:
Theorem A1.
If G ^ ( s ) and q ^ ( s ) are analytic functions, and
L [ g ( x , ξ ) ; s ] = g ^ ( x , s )
as well as
L [ f ( ξ , t ) ; s ] = 0 f ( ξ , t ) e s t d t = G ^ ( s ) e ξ q ^ ( s ) ,
exist, then
G ^ ( s ) g ^ ( x , q ^ ( s ) ) = 0 0 g ( x , ξ ) f ( ξ , t ) d ξ e s t d t .
From Efros theorem appears that
L 1 [ G ^ ( s ) g ^ ( x , q ^ ( s ) ) ; t ] = 0 g ( x , ξ ) f ( ξ , t ) d ξ .

Appendix D. Derivation of the Efficiencies E 1 and E 2

The derivation of E 1 and E 2 base on using the Laplace transform of the Lévy–Smirnov distribution Φ 1 / 2 ( σ ) = exp [ 1 / ( 4 σ ) ] / ( 2 π σ 3 / 2 ) , this is
e a z 1 / 2 = 0 d t e z t a 2 Φ 1 / 2 ( t a 2 ) = 0 d t e z t a 2 π t 3 / 2 e a 2 / ( 4 t ) .
Inserting it into Equation (30) with z = τ μ s 2 μ + s μ and into Equation (39) with z = τ s μ + 1 + s μ , and next changing the order of integrals we get
E j = 0 d t x 0 2 π K t 3 / 2 e x 0 2 / ( 4 K t ) Int j , j = 1 , 2 .
The auxiliary integral Int j for the efficiencies E j respectively read
Int 1 = 0 d s e t τ μ s 2 μ t s μ and Int 2 = 0 d s e t τ s μ + 1 t s μ .
At first, we calculate Int 1 by setting s 2 μ = y , which can be expressed as
Int 1 = 1 2 μ 0 y 1 2 μ 1 e t τ μ y t y 1 / 2 = 1 μ Γ 1 + 1 2 μ e t / ( 8 τ μ ) D 1 / μ t 2 τ μ ,
where we used [66] (Equation (2.2.1.16)). The parabolic cylinder function is standard denoted as D 1 / μ ( · ) . Then, we substitute Equation (A5) into Equation (A4) and change t 1 / 2 onto ξ . That leads to
E 1 = x 0 μ π K 2 1 / μ Γ 1 + 1 2 μ 0 d ξ ξ 1 μ 2 e x 0 2 4 K ξ 2 + ξ 2 8 τ μ D 1 / μ ξ 2 τ μ ,
which, after using [67] (Equation (2.11.4.4)), gives the lower formula in Equation (30).
Proceeding analogically, we derive the explicit form of E 2 . Now, setting s 1 / ( 1 + μ ) = y , we can present Int 2 as
Int 2 = 1 1 + μ 0 d y y μ 1 + μ e t τ y t y μ / ( 1 + μ ) .
This integral can be calculated for fractional μ = l / k , such that l < k by employing [66] (Equation (2.2.1.22)) and then transforming the results into the more transparent language of the Fox H function according to formula [68] (Equation (8.3.2.22)). That gives
Int 2 = 1 μ ( t τ ) 1 1 + μ H 1 , 1 1 , 1 τ t 1 / μ [ 1 , 1 + 1 μ ] [ 1 1 + μ , 1 ] ,
which substitutes into Equation (A4) for j = 2 and after changing t 1 onto ξ yields
E 2 = x 0 2 μ K π τ 1 1 + μ 0 ξ 3 2 + 1 1 + μ e x 0 2 4 K ξ H 1 , 1 1 , 1 τ ξ 1 / μ [ 1 , 1 + 1 μ ] [ 1 1 + μ , 1 ] d ξ = 1 μ τ π 4 K x 0 2 2 4 K τ μ τ x 0 2 1 1 + μ H 2 , 1 1 , 2 4 K τ μ x 0 2 1 / μ [ 3 2 1 1 + μ , 1 μ ] , [ 1 , 1 + 1 μ ] [ 1 1 + μ , 1 ] ,
where we used [68] (Equation (2.25.2.3)). Next, applying [68] (Equation (8.3.1.8)) allows us to get the lower formula in Equation (39).
A comparison of the expressions for the efficiencies in each case as a function of μ is shown in Figure A1. In the FTE-I case, the efficiency decreases with μ , ranging from high values for small μ to low values for μ close to one. In contrast, in the FTE-II case, although efficiency decreases with μ , it is maintained close to 1 for the whole range of values.
Figure A1. Comparison of the efficiencies for FTE-I and FTE-II from Equations (30) and (39), for K = 1 , x 0 = 2 , τ = 1 .
Figure A1. Comparison of the efficiencies for FTE-I and FTE-II from Equations (30) and (39), for K = 1 , x 0 = 2 , τ = 1 .
Entropy 26 00665 g0a1

References

  1. Heaviside, O. On induction between parallel wires. J. Soc. Telegr. Eng. 1880, 9, 427–458. [Google Scholar] [CrossRef]
  2. Cattaneo, C.R. Sulla conduzione del calore. Atti. Sem. Mat. Fis. Univ. Modena 1948, 3, 83. [Google Scholar]
  3. Cattaneo, C.R. Sur une forme de l’equation de la chaleur eliminant le paradoxe d’une propagation instantanee’. C. R. Acad. Sci. Paris 1958, 247, 431. [Google Scholar]
  4. Vernotte, P. Les paradoxes de la théories continue de l’equation de la chaleur. C. R. Acad. Sci. Paris 1958, 246, 3154–3155. [Google Scholar]
  5. Górska, K.; Horzela, A.; Lenzi, E.K.; Pagnini, G.; Sandev, T. Generalized Cattaneo (telegrapher’s) equations in modeling anomalous diffusion phenomena. Phys. Rev. E 2020, 102, 022128. [Google Scholar] [CrossRef] [PubMed]
  6. Stratton, J.A. Electromagnetic Theory; McGraw-Hill Book Co.: New York, NY, USA, 1941. [Google Scholar]
  7. Morse, P.M.; Feshbach, H. Methods of Theoretical Physics; McGraw-Hill Book Co.: New York, NY, USA, 1953. [Google Scholar]
  8. Masoliver, J. Telegraphic processes with stochastic resetting. Phys. Rev. E 2019, 99, 012121. [Google Scholar] [CrossRef] [PubMed]
  9. Masoliver, J.; Weiss, G.H. Finite-velocity diffusion. Eur. J. Phys. 1996, 17, 190. [Google Scholar] [CrossRef]
  10. Masoliver, J.; Lindenberg, K. Continuous time persistent random walk: A review and some generalizations. Eur. J. Phys. B 2017, 90, 107. [Google Scholar] [CrossRef]
  11. Kac, M. A stochastic model related to the telegrapher’s equation. Rocky Mt. J. Math. 1974, 4, 497. [Google Scholar] [CrossRef]
  12. Weiss, G.H. Some applications of persistent random walks and the telegrapher’s equation. Phys. A 2002, 311, 381. [Google Scholar] [CrossRef]
  13. Górska, K.; Horzela, A. Subordination and memory dependent kinetics in diffusion and relaxation phenomena. Fract. Calc. Appl. Anal. 2023, 26, 480. [Google Scholar] [CrossRef]
  14. Górska, K. Integral decomposition for the solutions of the generalized Cattaneo equation. Phys. Rev. E 2021, 104, 024113. [Google Scholar] [CrossRef] [PubMed]
  15. Compte, A.; Metzler, R. The generalized Cattaneo equation for the description of anomalous transport processes. J. Phys. A Math. Gen. 1997, 30, 7277–7289. [Google Scholar] [CrossRef]
  16. Kosztołowicz, T. Cattaneo-type subdiffusion-reaction equation. Phys. Rev. E 2014, 90, 042151. [Google Scholar] [CrossRef] [PubMed]
  17. D’Ovidio, M.; Polito, F. Fractional Diffusion—Telegraph Equations and Their Associated Stochastic Solutions. Theory Probab. Appl. 2018, 62, 552–574. [Google Scholar] [CrossRef]
  18. Awad, E.; Metzler, R. Crossover dynamics from superdiffusion to subdiffusion: Models and solutions. Fract. Calc. Appl. Anal. 2020, 23, 55–102. [Google Scholar] [CrossRef]
  19. Michelitsch, T.M.; Polito, F.; Riascos, A.P. Squirrels can remember little: A random walk with jump reversals induced by a discrete-time renewal process. Commun. Nonlinear Sci. Numer. Simul. 2023, 118, 107031. [Google Scholar] [CrossRef]
  20. Cáceres, M.O. Finite-velocity diffusion in random media. J. Stat. Phys. 2020, 179, 729–747. [Google Scholar] [CrossRef]
  21. Cáceres, M.O. Localization of plane waves in the stochastic telegrapher’s equation. Phys. Rev. E 2022, 105, 014110. [Google Scholar] [CrossRef]
  22. Cáceres, M.O.; Nizama, M. Stochastic telegrapher’s approach for solving the random Boltzmann-Lorentz gas. Phys. Rev. E 2022, 105, 044131. [Google Scholar] [CrossRef]
  23. Stadje, W.; Zacks, S. Telegraph processes with random velocities. J. Appl. Probab. 2004, 41, 665–678. [Google Scholar] [CrossRef]
  24. Masó-Puigdellosas, A.; Campos, D.; Méndez, V. Transport properties and first-arrival statistics of random motion with stochastic reset times. Phys. Rev. E 2019, 99, 012141. [Google Scholar] [CrossRef] [PubMed]
  25. Evans, M.R.; Majumdar, S.N. Diffusion with resetting in arbitrary spatial dimension. J. Phys. A Math. Theor. 2014, 47, 285001. [Google Scholar] [CrossRef]
  26. Bodrova, A.S.; Chechkin, A.V.; Sokolov, I.M. Scaled Brownian motion with renewal resetting. Phys. Rev. E 2019, 100, 012120. [Google Scholar] [CrossRef] [PubMed]
  27. Tal-Friedman, O.; Pal, A.; Sekhon, A.; Reuveni, S.; Roichman, Y. Experimental realization of diffusion with stochastic resetting. J. Phys. Chem. Lett. 2020, 11, 7350–7355. [Google Scholar] [CrossRef] [PubMed]
  28. Besga, B.; Bovon, A.; Petrosyan, A.; Majumdar, S.N.; Ciliberto, S. Optimal mean first-passage time for a Brownian searcher subjected to resetting: Experimental and theoretical results. Phys. Rev. Res. 2020, 2, 032029. [Google Scholar] [CrossRef]
  29. Schilling, R.L.; Song, R.; Vondraček, Z. Bernstein Functions: Theory and Applications; Walter de Gruyter: Berlin, Germany, 2009. [Google Scholar]
  30. Masoliver, J.; Montero, M. Anomalous diffusion under stochastic resettings: A general approach. Phys. Rev. E 2019, 100, 042103. [Google Scholar] [CrossRef] [PubMed]
  31. Kuśmierz, Ł.; Gudowska-Nowak, E. Subdiffusive continuous-time random walks with stochastic resetting. Phys. Rev. E 2019, 99, 052116. [Google Scholar] [CrossRef] [PubMed]
  32. Stanislavsky, A.; Weron, A. Optimal non-Gaussian search with stochastic resetting. Phys. Rev. E 2021, 104, 014125. [Google Scholar] [CrossRef] [PubMed]
  33. Singh, R.K.; Gorska, K.; Sandev, T. General approach to stochastic resetting. Phys. Rev. E 2022, 105, 064133. [Google Scholar] [CrossRef] [PubMed]
  34. Pollard, H. The representation of e−xλ as a Laplace integral. Bull. Am. Math. Soc. 1946, 52, 908. [Google Scholar] [CrossRef]
  35. Hilfer, R. H-function representations for stretched exponential relaxation and non-Debye susceptibilities in glassy systems. Phys. Rev. E 2002, 65, 061510. [Google Scholar] [CrossRef] [PubMed]
  36. Penson, K.A.; Górska, K. Exact and Explicit Probability Densities for One-Sided Lévy Stable Distributions. Phys. Rev. Lett. 2010, 105, 210604. [Google Scholar] [CrossRef] [PubMed]
  37. Fogedby, H.C. Langevin equations for continuous time Lévy flights. Phys. Rev. E 1994, 50, 1657. [Google Scholar] [CrossRef] [PubMed]
  38. Baule, A.; Friedrich, R. Joint probability distributions for a class of non-Markovian processes. Phys. Rev. E 2005, 71, 026101. [Google Scholar] [CrossRef] [PubMed]
  39. Chechkin, A.V.; Hofmann, M.; Sokolov, I.M. Continuous-time random walk with correlated waiting times. Phys. Rev. E 2009, 80, 031112. [Google Scholar] [CrossRef] [PubMed]
  40. Redner, S. A Guide to First-Passage Processes; Cambridge University Press: Cambridge, UK, 2001. [Google Scholar]
  41. Bray, A.J.; Majumdar, S.N.; Schehr, G. Persistence and first-passage properties in nonequilibrium systems. Adv. Phys. 2013, 62, 225–361. [Google Scholar] [CrossRef]
  42. Pal, A.; Kundu, A.; Evans, M.R. Diffusion under time-dependent resetting. J. Phys. A Math. Theor. 2016, 49, 225001. [Google Scholar] [CrossRef]
  43. Palyulin, V.V.; Chechkin, A.V.; Metzler, R. Space-fractional Fokker–Planck equation and optimization of random search processes in the presence of an external bias. J. Stat. Mech. 2014, 2014, P11031. [Google Scholar] [CrossRef]
  44. Pal, A.; Stojkoski, V.; Sandev, T. Random resetting in search problems. arXiv 2023, arXiv:2310.12057. [Google Scholar]
  45. Sandev, T.; Kocarev, L.; Metzler, R.; Chechkin, A. Stochastic dynamics with multiplicative dichotomic noise: Heterogeneous telegrapher’s equation, anomalous crossovers and resetting. Chaos Solitons Fractals 2022, 165, 112878. [Google Scholar] [CrossRef]
  46. Shkilev, V. Continuous-time random walk under time-dependent resetting. Phys. Rev. E 2017, 96, 012126. [Google Scholar] [CrossRef] [PubMed]
  47. Evans, M.E.; Majumdar, S.N. Diffusion with stochastic resetting. Phys. Rev. Lett. 2011, 106, 160601. [Google Scholar] [CrossRef] [PubMed]
  48. Kochubei, A.N. General Fractional Calculus, Evolution Equations, and Renewal Processes. Integral Equ. Oper. Theory 2011, 71, 583–600. [Google Scholar] [CrossRef]
  49. Bodrova, A.S.; Sokolov, I.M. Resetting processes with noninstantaneous return. Phys. Rev. E 2020, 101, 052130. [Google Scholar] [CrossRef]
  50. Tucci, G.; Gambassi, A.; Majumdar, S.N.; Schehr, G. First-passage time of run-and-tumble particles with noninstantaneous resetting. Phys. Rev. E 2022, 106, 044127. [Google Scholar] [CrossRef] [PubMed]
  51. Radice, M. One-dimensional telegraphic process with noninstantaneous stochastic resetting. Phys. Rev. E 2021, 104, 044126. [Google Scholar] [CrossRef] [PubMed]
  52. Tal-Friedman, O.; Roichman, Y.; Reuveni, S. Diffusion with partial resetting. Phys. Rev. E 2022, 106, 054116. [Google Scholar] [CrossRef]
  53. Di Bello, C.; Chechkin, A.V.; Hartmann, A.K.; Palmowski, Z.; Metzler, R. Time-dependent probability density function for partial resetting dynamics. New J. Phys. 2023, 25, 082002. [Google Scholar] [CrossRef]
  54. Christou, C.; Schadschneider, A. Diffusion with resetting in bounded domains. J. Phys. A Math. Theor. 2015, 48, 285003. [Google Scholar] [CrossRef]
  55. Pal, A.; Prasad, V. First passage under stochastic resetting in an interval. Phys. Rev. E 2019, 99, 032123. [Google Scholar] [CrossRef] [PubMed]
  56. Tucci, G.; Gambassi, A.; Gupta, S.; Roldán, É. Controlling particle currents with evaporation and resetting from an interval. Phys. Rev. Res. 2020, 2, 043138. [Google Scholar] [CrossRef]
  57. Julián-Salgado, P.; Dagdug, L.; Boyer, D. Diffusion with two resetting points. Phys. Rev. E 2024, 109, 024134. [Google Scholar] [CrossRef] [PubMed]
  58. Mendez, V.; Flaquer-Galmés, R.; Campos, D. First-passage time of a Brownian searcher with stochastic resetting to random positions. Phys. Rev. E 2024, 109, 044134. [Google Scholar] [CrossRef] [PubMed]
  59. Cáceres, M.O.; Nizama, M.; Pennini, F. Fisher and Shannon Functionals for Hyperbolic Diffusion. Entropy 2023, 25, 1627. [Google Scholar] [CrossRef] [PubMed]
  60. Prabhakar, T.R. A singular integral equation with a generalized Mittag Leffler function in the kernel. Yokohama Math. J. 1971, 19, 7. [Google Scholar]
  61. Efros, A.M. The application of the operational calculus to the analysis. in Russian. Mat. Sborni 1935, 42, 699. [Google Scholar]
  62. Włodarski, L. Sur une formule de Efros. Stud. Math. 1952, 13, 183. [Google Scholar] [CrossRef]
  63. Graf, U. Applied Laplace Transforms and z-Transforms for Sciences and Engineers; Birkhäuser: Basel, Switzerland, 2004. [Google Scholar]
  64. Górska, K.; Penson, K.A. Lévy stable distributions via associated integral transform. J. Math. Phys. 2012, 53, 053302. [Google Scholar] [CrossRef]
  65. Apelblat, A.; Mainardi, F. Application of the Efros theorem to the function represented by the inverse Laplace transform of s−μexp(-sν). Symmetry 2021, 13, 354. [Google Scholar] [CrossRef]
  66. Prudnikov, A.P.; Brychkov, Y.A.; Marichev, O.I. Integrals and Series. Direct Laplace Transforms; Gordon and Breach: London, UK, 1992; Volume 4. [Google Scholar]
  67. Prudnikov, A.P.; Brychkov, Y.A.; Marichev, O.I. Integrals and Series. Special Functions; Taylor & Francis: Abingdon, UK, 1998; Volume 2. [Google Scholar]
  68. Prudnikov, A.P.; Brychkov, Y.A.; Marichev, O.I. Integrals and Series. More Special Functions; Gordon and Breach: London, UK, 1990; Volume 3. [Google Scholar]
Figure 1. MSD (25) under the effects of stochastic resetting for different values of the resetting rate, r = { 0 , 0.1 , 1.0 , 10.0 } . The case r = 0 corresponds to the MSD (22).
Figure 1. MSD (25) under the effects of stochastic resetting for different values of the resetting rate, r = { 0 , 0.1 , 1.0 , 10.0 } . The case r = 0 corresponds to the MSD (22).
Entropy 26 00665 g001
Figure 2. MFPT (31) for x 0 = 2 , τ = 1 , K = 1 , μ = 1 / 2 (blue solid line), μ = 3 / 4 (red dashed line), μ = 1 (black dotted line).
Figure 2. MFPT (31) for x 0 = 2 , τ = 1 , K = 1 , μ = 1 / 2 (blue solid line), μ = 3 / 4 (red dashed line), μ = 1 (black dotted line).
Entropy 26 00665 g002
Figure 3. Optimal resetting rate r * versus μ , obtained by numerically solving Equation (32).
Figure 3. Optimal resetting rate r * versus μ , obtained by numerically solving Equation (32).
Entropy 26 00665 g003
Figure 4. MSD (38) under the effects of stochastic resetting for different values of μ and the resetting rate, r = { 0 ,   0.1 ,   1.0 ,   10.0 } . The case r = 0 corresponds to the MSD (36).
Figure 4. MSD (38) under the effects of stochastic resetting for different values of μ and the resetting rate, r = { 0 ,   0.1 ,   1.0 ,   10.0 } . The case r = 0 corresponds to the MSD (36).
Entropy 26 00665 g004
Figure 5. MFPT (40) for x 0 = 2 , τ = 1 , K = 1 , μ = 1 / 2 (blue solid line), μ = 3 / 4 (red dashed line), μ = 1 (black dotted line).
Figure 5. MFPT (40) for x 0 = 2 , τ = 1 , K = 1 , μ = 1 / 2 (blue solid line), μ = 3 / 4 (red dashed line), μ = 1 (black dotted line).
Entropy 26 00665 g005
Figure 6. Optimal resetting rate r * versus μ , obtained by numerically solving Equation (41).
Figure 6. Optimal resetting rate r * versus μ , obtained by numerically solving Equation (41).
Entropy 26 00665 g006
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Górska, K.; Sevilla, F.J.; Chacón-Acosta, G.; Sandev, T. Fractional Telegrapher’s Equation under Resetting: Non-Equilibrium Stationary States and First-Passage Times. Entropy 2024, 26, 665. https://doi.org/10.3390/e26080665

AMA Style

Górska K, Sevilla FJ, Chacón-Acosta G, Sandev T. Fractional Telegrapher’s Equation under Resetting: Non-Equilibrium Stationary States and First-Passage Times. Entropy. 2024; 26(8):665. https://doi.org/10.3390/e26080665

Chicago/Turabian Style

Górska, Katarzyna, Francisco J. Sevilla, Guillermo Chacón-Acosta, and Trifce Sandev. 2024. "Fractional Telegrapher’s Equation under Resetting: Non-Equilibrium Stationary States and First-Passage Times" Entropy 26, no. 8: 665. https://doi.org/10.3390/e26080665

APA Style

Górska, K., Sevilla, F. J., Chacón-Acosta, G., & Sandev, T. (2024). Fractional Telegrapher’s Equation under Resetting: Non-Equilibrium Stationary States and First-Passage Times. Entropy, 26(8), 665. https://doi.org/10.3390/e26080665

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop