Next Article in Journal
Switching-Jumps-Dependent Quasi-Synchronization Criteria for Fractional-Order Memrisive Neural Networks
Next Article in Special Issue
Approximating the First Passage Time Density of Diffusion Processes with State-Dependent Jumps
Previous Article in Journal
Fuzzy Analysis of Artificial Drug Transmission Model with Optimal Control
Previous Article in Special Issue
Ergodic Stationary Distribution and Threshold Dynamics of a Stochastic Nonautonomous SIAM Epidemic Model with Media Coverage and Markov Chain
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

On the Absorbing Problems for Wiener, Ornstein–Uhlenbeck, and Feller Diffusion Processes: Similarities and Differences

by
Virginia Giorno
and
Amelia G. Nobile
*,†
Dipartimento di Informatica, Università degli Studi di Salerno, Via Giovanni Paolo II n. 132, Fisciano, 84084 Salerno, Italy
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Fractal Fract. 2023, 7(1), 11; https://doi.org/10.3390/fractalfract7010011
Submission received: 2 November 2022 / Revised: 19 December 2022 / Accepted: 20 December 2022 / Published: 24 December 2022
(This article belongs to the Special Issue Stochastic Modeling in Biological System)

Abstract

:
For the Wiener, Ornstein–Uhlenbeck, and Feller processes, we study the transition probability density functions with an absorbing boundary in the zero state. Particular attention is paid to the proportional cases and to the time-homogeneous cases, by obtaining the first-passage time densities through the zero state. A detailed study of the asymptotic average of local time in the presence of an absorbing boundary is carried out for the time-homogeneous cases. Some relationships between the transition probability density functions in the presence of an absorbing boundary in the zero state and between the first-passage time densities through zero for Wiener, Ornstein–Uhlenbeck, and Feller processes are proven. Moreover, some asymptotic results between the first-passage time densities through zero state are derived. Various numerical computations are performed to illustrate the role played by parameters.

1. Introduction and Background

Diffusion models are widely used to describe dynamical systems in economics, finance, biology, genetics, physics, engineering, neuroscience, queueing, and other fields (cf. Bailey [1], Ricciardi [2], Gardiner [3], Stirzaker [4], Janssen et al. [5], Pavliotis [6]). In various applications, it is useful to consider diffusion processes with linear infinitesimal drift and linear infinitesimal variance. This class incorporates Wiener, Ornstein–Uhlenbeck, and Feller diffusion processes. In population dynamics, these processes can be used to describe the growth of a population and the zero state represents the absorbing extinction threshold. With this aim, we study the absorbing problem for linear diffusion processes.
In the remaining part of this section, we shall briefly review some background results on the absorbing problems that will be used in the next sections for Wiener, Ornstein–Uhlenbeck and Feller diffusion processes.
Let { Z ( t ) , t t 0 } be a time-inhomogeneous diffusion (TNH-D) process with state-space D = ( r 1 , r 2 ) , which satisfies the stochastic differential equation
d Z ( t ) = ζ 1 [ Z ( t ) , t ] d t + ζ 2 [ Z ( t ) , t ] d W ( t ) , Z ( t 0 ) = x 0 ,
with ζ 1 ( x , t ) and ζ 2 ( x , t ) denoting, respectively, the infinitesimal drift and the infinitesimal variance of Z ( t ) and where W ( t ) is a standard Brownian motion. Often, D = ( , + ) , with ± unattainable endpoints, but in some cases Z ( t ) is confined to the state space D = ( 0 , + ) and in the zero state is imposed an absorbing condition.
When the endpoints ± of D are unattainable boundaries, the transition probability density function (PDF) f Z ( x , t | x 0 , t 0 ) is the solution of the backward Kolmogorov equation (cf. Dynkin [7])
f Z ( x , t | x 0 , t 0 ) t 0 + ζ 1 ( x 0 , t 0 ) f Z ( x , t | x 0 , t 0 ) x 0 + 1 2 ζ 2 ( x 0 , t 0 ) 2 f Z ( x , t | x 0 , t 0 ) x 0 2 = 0 ,
with the initial delta condition lim t 0 t f Z ( x , t | x 0 , t 0 ) = δ ( x x 0 ) . In the backward Kolmogovov Equation (1), the forward variables x and t are constant and enter only through the initial and boundary conditions.
We remark that the PDF f Z ( x , t | x 0 , t 0 ) is also solution of a forward Kolmogorov equation, also known as the Fokker–Planck equation (cf. Dynkin [7]), in which the backward variables x 0 and t 0 are essentially constant. In this paper, we choose to use the Kolmogorov backward equation because we will address absorption problems. Indeed, if one is interested to the first-passage time distribution through a fixed state S as a function of the initial position x 0 , then the backward Kolmogorov equation provides the most appropriate method (cf. Cox and Miller [8]).
For a diffusion process Z ( t ) , the first-passage time (FPT) problem can be reduced to estimate the density of the random variable
T Z ( S | x 0 , t 0 ) = inf t t 0 { t : Z ( t ) S } , Z ( t 0 ) = x 0 < S , inf t t 0 { t : Z ( t ) S } , Z ( t 0 ) = x 0 > S ,
which describes the FPT of Z ( t ) through the state S starting from Z ( t 0 ) = x 0 S .
The FPT problem plays an important role in various biological applications. For instance, in the context of population dynamics the FPT problem is suitable to model population’s extinction or persistence (see Bailey [1], Ricciardi [2], Allen [9,10]).
Let g Z ( S , t | x 0 , t 0 ) = d P { T Z ( S | x 0 , t 0 ) t } / d t be the FPT density, being P { T Z ( S | x 0 , t 0 ) t } the distribution function of the random variable T Z ( S | x 0 , t 0 ) . If the endpoints of D are unattainable boundaries, the densities f Z ( x , t | x 0 , t 0 ) and g Z ( S , t | x 0 , t 0 ) are related by the following renewal equation (cf. Blake and Lindsey [11]):
f Z ( x , t | x 0 , t 0 ) = t 0 t g Z ( S , τ | x 0 , t 0 ) f Z ( x , t | S , τ ) d τ , ( x 0 < S x ) or ( x S < x 0 ) .
Equation (2) indicates that any sample path that reaches x S [ x S ] , after starting from x 0 < S [ x 0 > S ] at time t 0 , must necessarily cross S for the first time at some intermediate instant τ ( t 0 , t ) .
For diffusion processes, closed form expressions for FPT densities through constant boundaries are not available, except in some special cases (see Ricciardi et al. [12], Ding and Rangarajan [13], Molini et al. [14], Giorno and Nobile [15], Masoliver [16]). In particular, closed form expressions are available in the following cases: (i) the Wiener process through an arbitrary constant boundary; (ii) the Ornstein–Uhlenbeck process through the boundary in which the drift vanishes; and (iii) the Feller process through the zero state. In the literature many efforts have been devoted to determining the asymptotic behavior of FPT density and its moments for large boundaries or large times and to search efficient numerical and simulation methods to estimate the FPT densities (cf. Ricciardi et al. [12], Linetsky [17]). Furthermore, the FPT problems play a relevant role also in the context of fractional processes (see, for instance, Guo et al. [18], Wiese [19], Abundo [20], Leonenko and Pirozzi [21]).
For a TNH-D process Z ( t ) confined to interval ( 0 , + ) , with 0 absorbing boundary and + unattainable boundary, we denote with
a Z ( x , t | x 0 , t 0 ) = x P { Z ( t ) x ; Z ( θ ) > 0 , θ < t | Z ( t 0 ) = x 0 } , x > 0 , x 0 > 0
the PDF of Z ( t ) with an absorbing condition in the zero state. The PDF a Z ( x , t | x 0 , t 0 ) satisfies the Kolmogorov Equation (1) with the initial condition lim t 0 t a Z ( x , t | x 0 , t 0 ) = δ ( x x 0 ) and the absorbing condition lim x 0 0 a Z ( x , t | x 0 , t 0 ) = 0 .
The densities f Z ( x , t | x 0 , t 0 ) , g Z ( 0 , t | x 0 , t 0 ) , and a Z ( x , t | x 0 , t 0 ) are related by the following integral equations (cf. Siegert [22]):
a Z ( x , t | x 0 , t 0 ) = f Z ( x , t | x 0 , t 0 ) t 0 t g Z ( 0 , θ | x 0 , t 0 ) f Z ( x , t | 0 , θ ) d θ , x 0 > 0 , x > 0 ,
0 + a Z ( x , t | x 0 , t 0 ) d x + t 0 t g Z ( 0 , θ | x 0 , t 0 ) d θ = 1 , x 0 > 0 .
In the context of population dynamics, the first integral in (4) gives the survival probability, i.e., the probability that the trajectories of the process Z ( t ) are not absorbed in the zero state in ( t 0 , t ) . Moreover, from (4) one obtains the FPT density
g Z ( 0 , t | x 0 , t 0 ) = t 0 + a Z ( x , t | x 0 , t 0 ) d x , x 0 > 0 ,
and the ultimate FPT probability of Z ( t ) through the zero-state
P Z ( 0 | x 0 , t 0 ) = t 0 + g Z ( 0 , τ | x 0 , t 0 ) d τ = 1 lim t + 0 + a Z ( x , t | x 0 , t 0 ) d x , x 0 > 0 .
In population dynamics, g Z ( 0 , t | x 0 , t 0 ) in Equation (5) represents the density of the time required to reach the zero state for the first time (extinction density); instead, P Z ( 0 | x 0 , t 0 ) in Equation (6) provides the probability that the population will become extinct sooner or later.
For a TNH-D process Z ( t ) , the local time L ( t , x | t 0 ) at an interior state x D is a random variable defined as (cf. Karlin and Taylor [23], Aït-Sahalia and Park [24]):
L ( t , x | t 0 ) = lim ε 0 1 2 ε t 0 t 1 { | Z ( θ ) x | ε } d θ , t > t 0 ,
where, for ε > 0 , we have set
1 { | Z ( θ ) x | ε } = 1 , | Z ( θ ) x | ε , 0 , otherwise .
The asymptotic average of the local time in the presence of an absorbing boundary in the zero state, for x > 0 and x 0 > 0 is:
L Z ( x | x 0 , t 0 ) = lim t + E [ L ( t , x | t 0 ) | Z ( t 0 ) = x 0 ] = t 0 + a Z ( x , θ | x 0 , t 0 ) d θ .
For a time-homogeneous diffusion (TH-D) process Z ( t ) one has ζ 1 ( x , t ) = ζ 1 ( x ) and ζ 2 ( x , t ) = ζ 2 ( x ) for all t. In this case, the classification of the endpoints of the state space D , due to Feller [25,26], is based on integrability properties of the functions
h Z ( x ) = exp 2 x ζ 1 ( u ) ζ 2 ( u ) d u , s Z ( x ) = 2 ζ 2 ( u ) h Z ( u ) , x D ,
called scale function and speed density, respectively. Such functions allow us to determine the FPT moments for TH-D processes thanks to the Siegert formula (cf. Siegert [22]). Specifically, if Z ( t ) is a TH-D process with state space D = ( r 1 , r 2 ) , for n = 1 , 2 , it results in
  • for x 0 < S , if P Z ( S | x 0 ) = 0 + g Z ( S , t | x 0 ) d t = 1 and if r 1 z s Z ( u ) d u converges one has:
    t n ( Z ) ( S | x 0 ) = 0 + t n g Z ( S , t | x 0 ) d t = n x 0 S d z h Z ( z ) r 1 z s Z ( u ) t n 1 ( S | u ) d u , x 0 < S ,
  • for x 0 > S , if P Z ( S | x 0 ) = 1 , and if z r 2 s Z ( u ) d u converges one has
    t n ( Z ) ( S | x 0 ) = 0 + t n g Z ( S , t | x 0 ) d t = n S x 0 d z h Z ( z ) z r 2 s Z ( u ) t n 1 ( S | u ) d u , x 0 > S ,
with t 0 ( Z ) ( S | x 0 ) = P Z ( S | x 0 ) .
In the sequel, for the FPT of TH-D process Z ( t ) we denote by
Var ( Z ) ( S | x 0 ) = t 2 ( Z ) ( S | x 0 ) [ t 1 ( Z ) ( S | x 0 ) ] 2 , Cv ( Z ) ( S | x 0 ) = Var ( Z ) ( S | x 0 ) t 1 ( Z ) ( S | x 0 ) , Σ ( Z ) ( S | x 0 ) = t 3 ( Z ) ( S | x 0 ) 3 t 1 ( Z ) ( S | x 0 ) t 2 ( Z ) ( S | x 0 ) + 2 [ t 1 ( Z ) ( S | x 0 ) ] 3 [ Var ( Z ) ( S | x 0 ) ] 3 / 2 ·
the variance, the coefficient of variation, and the skewness, respectively.
For a TH-D process in ( 0 , + ) , with 0 absorbing boundary, if x 0 > 0 and x > 0 the asymptotic average of the local time is (cf. Giorno and Nobile [27]):
L Z ( x | x 0 ) = s Z ( x ) 0 x 0 x h Z ( z ) d z , + unattainable , nonattracting , s Z ( x ) P Z ( 0 | x x 0 ) 0 x 0 x h Z ( z ) d z , + unattainable , attracting ,
where x 0 x = min ( x 0 , x ) and x 0 x = max ( x 0 , x ) .
For a TH-D process Z ( t ) , in the sequel we denote by
φ λ ( Z ) ( x | x 0 ) = 0 + e λ t φ Z ( x , t | x 0 ) d t
the Laplace transform (LT) of the function φ Z ( x , t | x 0 ) .

Plan of the Paper

In Section 2, we consider the time-inhomogeneous Wiener (TNH-W) process X ( t ) , with infinitesimal drift and infinitesimal variance A 1 ( t ) = β ( t ) and A 2 ( t ) = σ 2 ( t ) , respectively. For β ( t ) = γ σ 2 ( t ) , with γ R , we determine the PDF a X ( x , t | x 0 , t 0 ) and the FPT density g X ( 0 , t | x 0 , t 0 ) . Furthermore, for the time-homogeneous Wiener (TH-W) process, the FPT moments through a boundary S R and the asymptotic average of the local time are studied.
In Section 3, we take into account the time-inhomogeneous Ornstein–Uhlenbeck (TNH-OU) process Y ( t ) , with infinitesimal drift and infinitesimal variance B 1 ( x , t ) = α ( t ) x + β ( t ) and B 2 ( t ) = σ 2 ( t ) , respectively. For β ( t ) = γ σ 2 ( t ) e A ( t | 0 ) , with γ R and A ( t | 0 ) = 0 t α ( u ) d u , we determine a Y ( x , t | x 0 , t 0 ) and g Y ( 0 , t | x 0 , t 0 ) . Moreover, for the TH-OU process, the FPT mean through a constant boundary and the asymptotic average of the local time are evaluated.
In Section 4, we consider the time-inhomogeneous Feller (TNH-F) process Z ( t ) with infinitesimal drift and infinitesimal variance C 1 ( x , t ) = α ( t ) x + β ( t ) and C 2 ( x , t ) = 2 r ( t ) x , respectively, with an absorbing boundary in the zero-state. For β ( t ) = ξ r ( t ) , with 0 ξ < 1 , we obtain a Z ( x , t | x 0 , t 0 ) and g Z ( 0 , t | x 0 , t 0 ) . Furthermore, for the TH-F process, the FPT mean through a constant boundary and the asymptotic average of the local time are examined.
We remark that time-inhomogeneous Wiener, Ornstein–Uhlenbeck and Feller diffusion processes are used in biological systems to model the growth of a population. In such a context, α ( t ) represents the growth intensity function and β ( t ) denotes the immigration/emigration intensity function. The functions σ 2 ( t ) (in Wiener and Ornstein–Uhlenbeck processes) and r ( t ) (in the Feller process) are the noise intensity functions and take into account the environmental fluctuations.
In Section 2, Section 3 and Section 4, by using Siegert Formulas (10) and (11), extensive computation are performed with MATHEMATICA to obtain the mean, the variance, the coefficient of variation, and the skewness of FPT for the TH-W, TH-OU, and TH-F processes for various choices of parameters. For these processes, some considerations on the asymptotic average of the local time in the presence of an absorbing boundary in the zero state are also made.
In Section 5, for β ( t ) = r ( t ) / 2 , some relationships between the PDF in the presence of an absorbing boundary in the zero state and between the FPT densities through zero for Wiener, Ornstein–Uhlenbeck and Feller processes are proved. Moreover, for β ( t ) = ξ r ( t ) ( 0 < ξ < 1 ) some asymptotic results for large times between the FPT densities are provided.

2. Wiener-Type Diffusion Process

Let { X ( t ) , t t 0 } , t 0 0 , be a TNH-W process, having infinitesimal drift and infinitesimal variance
A 1 ( t ) = β ( t ) , A 2 ( t ) = σ 2 ( t ) ,
with the state space R , where β ( t ) R and σ ( t ) > 0 are continuous functions.
The Wiener process arises as the mathematical limit of other stochastic processes, such as random walks (see Knight [28]). This process has been originally used in physics to model the motion of particles suspended in a fluid and it is still used as a mathematical model for various random phenomena in applied mathematics, economics, quantitative finance, evolutionary biology, and physics.
For t t 0 , the PDF of X ( t ) is normal,
f X ( x , t | x 0 , t 0 ) = 1 2 π V X ( t | t 0 ) exp x M X ( t | x 0 , t 0 ) 2 2 V X ( t | t 0 ) , x , x 0 R ,
with
M X ( t | x 0 , t 0 ) = x 0 + t 0 t β ( u ) d u , V X ( t | t 0 ) = t 0 t σ 2 ( u ) d u .
We now consider the TNH-W process X ( t ) , having infinitesimal moments given in (13), restricted to the state space ( 0 , + ) with 0 absorbing boundary; we denote by a X ( x , t | x 0 , t 0 ) its PDF. For the Wiener process X ( t ) in the presence of an absorbing boundary in the zero state, we analyze two cases: the proportional case with β ( t ) = γ σ 2 ( t ) , being γ R and σ ( t ) > 0 , and the time-homogeneous case.

2.1. Proportional Case for the Wiener Process

Proposition 1.
Let β ( t ) = γ σ 2 ( t ) , with γ R and σ ( t ) > 0 in (13). For the TNH-W process X ( t ) one has
a X ( x , t | x 0 , t 0 ) = f X ( x , t | x 0 , t 0 ) e 2 γ x f X ( x , t | x 0 , t 0 ) , x > 0 , x 0 > 0 ,
with f X ( x , t | x 0 , t 0 ) given in (14).
Proof. 
If β ( t ) = γ σ 2 ( t ) , from (14) the following symmetry relation holds,
f X ( x , t | 0 , t 0 ) = e 2 γ x f X ( x , t | 0 , t 0 ) , x R ,
so that from (3) one has
a X ( x , t | x 0 , t 0 ) = f X ( x , t | x 0 , t 0 ) e 2 γ x t 0 t g X ( 0 , τ | x 0 , t 0 ) f X ( x , t | 0 , τ ) d τ , x 0 > 0 , x > 0 .
Hence, by virtue of the renewal Equation (2), Equation (15) follows from (16). □
From (15), for β ( t ) = γ σ 2 ( t ) , with γ R and σ ( t ) > 0 , one explicitly obtains
a X ( x , t | x 0 , t 0 ) = 1 2 π V X ( t | t 0 ) [ exp x x 0 γ V X ( t | t 0 ) 2 2 V X ( t | t 0 ) e 2 γ x exp x + x 0 + γ V X ( t | t 0 ) 2 2 V X ( t | t 0 ) ] , x 0 > 0 , x > 0 .
We note that Equation (17) for t 0 = 0 is in agreement with Equation (25) in Molini et al. [14].
Proposition 2.
Under the assumptions of Proposition 1, for the TNH-W process X ( t ) one has
g X ( 0 , t | x 0 , t 0 ) = x 0 σ 2 ( t ) 2 π V X ( t | t 0 ) 3 exp x 0 + γ V X ( t | t 0 ) 2 2 V X ( t | t 0 ) , x 0 > 0 .
Moreover, if lim t + V X ( t | t 0 ) = + , the ultimate FPT probability of X ( t ) through zero is
P X ( 0 | x 0 , t 0 ) = t 0 + g X ( 0 , t | x 0 , t 0 ) d t = 1 , γ 0 , e . 2 γ x 0 , γ > 0 , x 0 > 0 .
Proof. 
For x 0 > 0 , from (17) one obtains
0 + a X ( x , t | x 0 , t 0 ) d x = 1 2 1 + Erf x 0 + γ V X ( t | t 0 ) 2 V X ( t | t 0 ) e 2 γ x 0 Erfc x 0 γ V X ( t | t 0 ) 2 V X ( t | t 0 ) ,
where Erf ( x ) = ( 2 / π ) 0 x e z 2 d z denotes the error function and Erfc ( x ) = 1 Erf ( x ) is the complementary error function. Hence, due to (5) and recalling (20), Equation (18) follows. Finally, if lim t + V X ( t | t 0 ) = + , Equation (19) follows, making use of (20) in (6) and by noting that
lim t + Erf x 0 + γ V X ( t | t 0 ) 2 V X ( t | t 0 ) = 1 , γ < 0 , 0 , γ = 0 , 1 , γ > 0 , lim t + Erfc x 0 γ V X ( t | t 0 ) 2 V X ( t | t 0 ) = 0 , γ < 0 , 1 , γ = 0 , 2 , γ > 0 ,
for any x 0 . □
Equation (19) shows that if β ( t ) = γ σ 2 ( t ) , with γ R and σ ( t ) > 0 in (13), the first-passage for the Wiener process through zero is a sure event for γ > 0 and x 0 > 0 .

2.2. Time-Homogeneous Case for the Wiener Process

We consider the TH-W process, obtained from (13) by setting β ( t ) = β and σ 2 ( t ) = σ 2 , with β R and σ > 0 . When β > 0 ( β < 0 ) the end point is a nonattracting (attracting) natural boundary and the end point + is an attracting (nonattracting) natural boundary. Instead, for β = 0 the end points and + are nonattracting natural boundaries. The scale function and the speed density, defined in (9) for the TH-W process X ( t ) are
h X ( x ) = exp 2 β σ 2 x , s X ( x ) = 2 σ 2 exp 2 β σ 2 x ,
respectively.
The FPT density of the TH-W process X ( t ) through the constant boundary S starting from x 0 is
g X ( S , t | x 0 ) = | S x 0 | σ 2 π t 3 exp ( S x 0 β t ) 2 2 σ 2 t , S x 0
and the ultimate FPT probability is
P X ( S | x 0 ) = 0 + g X ( S , t | x 0 ) d t = 1 , β = 0 or β ( S x 0 ) > 0 , exp 2 β ( S x 0 ) σ 2 , β ( S x 0 ) < 0 .
For β ( S x 0 ) > 0 , the FPT moments of the TH-W process X ( t ) are finite and from (22) one has
t n ( X ) ( S | x 0 ) = 2 | S x 0 | σ 2 π S x 0 β n 1 / 2 exp β ( S x 0 ) σ 2 K n 1 / 2 [ β ( S x 0 ) σ 2 ] , n = 1 , 2 ,
where K ν ( z ) denotes the modified Bessel function of the third kind, which can be expressed in terms of the modified Bessel function of first kind I ν ( z ) (see Abramowitz and Stegun [29], p. 375, n. 9.6.2),
K ν ( z ) = π 2 I ν ( z ) I ν ( z ) sin ( ν π ) , I ν ( z ) = k = 0 + 1 k ! Γ ( ν + k + 1 ) z 2 2 k + ν ,
where Γ ( ν ) = 0 + y ν 1 e y d y , with Re ν > 0 , is the Euler gamma function.
In particular, for β ( S x 0 ) > 0 the first three FPT moments of the TH-W process X ( t ) are
t 1 ( X ) ( S | x 0 ) = S x 0 β , t 2 ( X ) ( S | x 0 ) = S x 0 β 2 1 + σ 2 β ( S x 0 ) , t 3 ( X ) ( S | x 0 ) = S x 0 β 3 1 + 3 σ 2 β ( S x 0 ) + 3 σ 4 β 2 ( S x 0 ) 2 .
In Table 1 and Table 2, the mean t 1 ( X ) ( S | x 0 ) , the variance Var ( X ) ( S | x 0 ) , the coefficient of variation Cv ( X ) ( S | x 0 ) , and the skewness Σ ( X ) ( S | x 0 ) of the FPT are listed for x 0 = 4 , σ = 1 and some choices of β and S.
As shown in Table 1 and Table 2, for the TH-W process X ( t ) the coefficient of variation and the skewness of the FPT decrease when S moves away from x 0 .
Moreover, by setting β ( t ) = β , σ 2 ( t ) = σ 2 and γ = β / σ 2 in (17), for the TH-W process X ( t ) one has
a X ( x , t | x 0 ) = 1 2 π σ 2 t exp x x 0 β t 2 2 σ 2 t exp 2 β x σ 2 exp x + x 0 + β t 2 2 σ 2 t
with x 0 > 0 and x > 0 .
Proposition 3.
For the TH-W process X ( t ) , the asymptotic average of the local time is
L X ( x | x 0 ) = 0 + a X ( x , t | x 0 ) d t = 1 | β | exp β ( x x 0 ) σ 2 [ exp | β ( x x 0 ) | σ 2 exp | β ( x + x 0 ) | σ 2 ] , β 0 , | x + x 0 | σ 2 | x x 0 | σ 2 = 2 ( x 0 x ) σ 2 , β = 0 ,
with x 0 x = min ( x 0 , x ) and x 0 x = max ( x 0 , x ) .
Proof. 
Because + is a nonattracting boundary for β 0 and attracting for β > 0 , Equation (26) follows from (12) making use of (21) and (23). □
From (26), for β R and σ > 0 one has lim x 0 L X ( x | x 0 ) = 0 and
lim x + L X ( x | x 0 ) = 0 , β < 0 , 2 x 0 σ 2 , β = 0 , 1 e 2 β x 0 / σ 2 β , β > 0 .
We note that L X ( x | x 0 ) tends to zero as x increases if β < 0 , and it approaches a positive value when β 0 .
In Figure 1, the asymptotic average of the local time for the TH-W process X ( t ) is plotted for x 0 = 4 , σ = 1 and some choices of β . We note that L X ( x | x 0 ) tends to zero only if β < 0 , otherwise approaches to a positive value.

3. Ornstein–Uhlenbeck-Type Diffusion Process

Let { Y ( t ) , t t 0 } , t 0 0 , be a TNH-OU process, having infinitesimal drift and infinitesimal variance
B 1 ( x , t ) = α ( t ) x + β ( t ) , B 2 ( t ) = σ 2 ( t ) , x R ,
with state space R , where α ( t ) R , β ( t ) R , σ ( t ) > 0 are continuous functions. Note that when α ( t ) = 0 for all t, the process Y ( t ) identifies with the TNH-W process X ( t ) having infinitesimal moments (13).
Although the Ornstein–Uhlenbeck process has been originally used in physics to model the velocity of a Brownian particle (see Uhlenbeck and Ornstein [30]), it finds many applications in several scientific fields. In particular, the Ornstein–Uhlenbeck process is frequently proposed as a stochastic model for the single neuronal activity (see Ricciardi and Sacerdote [31], Lánský and Ditlevsen [32]). A wide field of applications of the Ornstein–Uhlenbeck process lies also in mathematical finance to model the evolution of the interest rate of financial markets (cf. Vasicek [33], Hull and White [34]).
The PDF of Y ( t ) is normal,
f Y ( x , t | x 0 , t 0 ) = 1 2 π V Y ( t | t 0 ) exp [ x M Y ( t | x 0 , t 0 ) ] 2 2 V Y ( t | t 0 ) , x , x 0 R ,
with
M Y ( t | x 0 , t 0 ) = x 0 e A ( t | t 0 ) + t 0 t β ( θ ) e A ( t | θ ) d θ , V Y ( t | t 0 ) = t 0 t σ 2 ( θ ) e 2 A ( t | θ ) d θ ,
being
A ( t | t 0 ) = t 0 t α ( θ ) d θ .
We now consider the TNH-OU process Y ( t ) , having infinitesimal moments given in (27), restricted to the state space ( 0 , + ) , with 0 absorbing boundary, and denote by a Y ( x , t | x 0 , t 0 ) its PDF. For the TNH-OU process Y ( t ) with 0 absorbing boundary, we take into account two cases: the proportional case in which β ( t ) = γ σ 2 ( t ) e A ( t | 0 ) , with γ R , α ( t ) R and σ ( t ) > 0 , and the time-homogeneous case.

3.1. Proportional Case for the Ornstein-Uhlenbeck Process

Proposition 4.
Let β ( t ) = γ σ 2 ( t ) e A ( t | 0 ) , with γ R , α ( t ) R , σ ( t ) > 0 in (27) and A ( t | 0 ) defined in (30). For the TNH-OU process Y ( t ) one has
a Y ( x , t | x 0 , t 0 ) = f Y ( x , t | x 0 , t 0 ) exp 2 γ x e A ( t | 0 ) f Y ( x , t | x 0 , t 0 ) , x > 0 , x 0 > 0 ,
with f Y ( x , t | x 0 , t 0 ) given in (28).
Proof. 
By choosing β ( t ) = γ σ 2 ( t ) e A ( t | 0 ) , from (28) the following symmetry relation holds,
f Y ( x , t | 0 , t 0 ) = exp 2 γ x e A ( t | 0 ) f Y ( x , t | 0 , t 0 ) , x R ,
so that from (3) one obtains
a Y ( x , t | x 0 , t 0 ) = f Y ( x , t | x 0 , t 0 ) exp 2 γ x e A ( t | 0 ) t 0 t g Y ( 0 , τ | x 0 , t 0 ) f Y ( x , t | 0 , τ ) d τ
for x 0 > 0 and x > 0 . Hence, by virtue of the renewal Equation (2), Equation (31) follows from (32). □
From (31), if β ( t ) = γ σ 2 ( t ) e A ( t | 0 ) , for x 0 > 0 and x > 0 one obtains
a Y ( x , t | x 0 , t 0 ) = 1 2 π V Y ( t | t 0 ) [ exp x x 0 e A ( t | t 0 ) γ e A ( t | 0 ) V Y ( t | t 0 ) 2 2 V Y ( t | t 0 ) exp 2 γ x e A ( t | 0 ) exp x + x 0 e A ( t | t 0 ) + γ e A ( t | 0 ) V Y ( t | t 0 ) 2 2 V Y ( t | t 0 ) ] .
Proposition 5.
Under the assumptions of Proposition 4, for the TNH-OU process Y ( t ) one has
g Y ( 0 , t | x 0 , t 0 ) = x 0 σ 2 ( t ) e A ( t | t 0 ) 2 π [ V Y ( t | t 0 ) ] 3 exp x 0 e A ( t | t 0 ) + γ e A ( t | 0 ) V Y ( t | t 0 ) 2 2 V Y ( t | t 0 ) , x 0 > 0 .
Furthermore, if lim t + [ e 2 A ( t | t 0 ) V Y ( t | t 0 ) ] = + , the ultimate FPT probability for x 0 > 0 is
P Y ( 0 | x 0 , t 0 ) = t 0 + g Y ( 0 , t | x 0 , t 0 ) d t = 1 , γ 0 , exp 2 γ x 0 e A ( t 0 | 0 ) , γ > 0 .
Proof. 
Recalling (33), one obtains
0 + a Y ( x , t | x 0 , t 0 ) d x = 1 2 [ 1 + Erf x 0 e A ( t | t 0 ) + γ e A ( t | 0 ) V Y ( t | t 0 ) 2 V Y ( t | t 0 ) exp 2 γ x 0 e A ( t 0 | 0 ) Erfc x 0 e A ( t | t 0 ) γ e A ( t | 0 ) V Y ( t | t 0 ) 2 V Y ( t | t 0 ) ] , x 0 > 0 .
By virtue of (5) and recalling (36), Equation (34) follows. Moreover, under the assumption lim t + [ e 2 A ( t | t 0 ) V Y ( t | t 0 ) ] = + , Equation (35) follows, making use of (36) in (6) by noting that
lim t + Erf x 0 e A ( t | t 0 ) + γ e A ( t | 0 ) V Y ( t | t 0 ) 2 V Y ( t | t 0 ) = 1 , γ < 0 , 0 , γ = 0 , 1 , γ > 0 , lim t + Erfc x 0 e A ( t | t 0 ) γ e A ( t | 0 ) V Y ( t | t 0 ) 2 V Y ( t | t 0 ) ] = 0 , γ < 0 , 1 , γ = 0 , 2 , γ > 0 ,
for any x 0 . □

3.2. Time-Homogeneous Case for the Ornstein–Uhlenbeck Process

We consider the TH-OU process Y ( t ) , by setting in (27) α ( t ) = α , β ( t ) = β , σ 2 ( t ) = σ 2 , with α 0 , β R and σ > 0 . The end points and + are nonattracting natural boundaries for α < 0 and attracting natural boundaries for α > 0 . The scale function and the speed density, defined in (9), for the TH-OU process Y ( t ) are
h Y ( x ) = exp α σ 2 x 2 + 2 β α x , s Y ( x ) = 2 σ 2 exp α σ 2 x 2 + 2 β α x ,
respectively. The LT of f Y ( x , t | x 0 ) is
f λ ( Y ) ( x | x 0 ) = 2 λ | α | 1 σ π | α | Γ λ 2 | α | Γ 1 2 + λ 2 | α | exp | α | 2 σ 2 x + β α 2 x 0 + β α 2 × D λ | α | 2 | α | σ x 0 x + β α D λ | α | 2 | α | σ x 0 x + β α , α < 0 , 2 λ α σ π α Γ 1 + λ 2 α Γ 1 2 + λ 2 α exp α 2 σ 2 x 0 + β α 2 x + β α 2 × D λ α 1 2 α σ x 0 x + β α D λ α 1 2 α σ x 0 x + β α , α > 0 ,
where D ν ( z ) is the parabolic cylinder function defined as (cf. Gradshteyn and Ryzhik [35], p. 1028, no. 9.240). We have
D ν ( z ) = 2 ν / 2 e z 2 / 4 π Γ 1 ν 2 Φ ν 2 , 1 2 ; z 2 2 z 2 π Γ ν 2 Φ 1 ν 2 , 3 2 ; z 2 2
in terms of Kummer’s confluent hypergeometric function
Φ ( a , c ; x ) = 1 + n = 1 + ( a ) n ( c ) n x n n ! ,
with ( a ) 0 = 1 and ( a ) n = a ( a + 1 ) ( a + n 1 ) for n = 1 , 2 , In the following, we will make use of the relations (cf. Gradshteyn and Ryzhik [35], p. 1030, no. 9.251 and no. 9.254).
D 0 ( x ) = e x 2 / 4 , D 1 ( x ) = x e x 2 / 4 , D 1 ( x ) = π 2 e x 2 / 4 Erfc x 2 .
For the TH-OU process, taking the Laplace transform in (2) and recalling (38), for x 0 S one has
g λ ( Y ) ( S | x 0 ) = exp { | α | 2 σ 2 x 0 + β α 2 S + β α 2 } D λ | α | sign ( x 0 S ) 2 | α | σ x 0 + β α D λ | α | sign ( x 0 S ) 2 | α | σ S + β α , α < 0 , exp { α 2 σ 2 x 0 + β α 2 S + β α 2 } D λ α 1 sign ( x 0 S ) 2 α σ x 0 + β α D λ α 1 sign ( x 0 S ) 2 α σ S + β α , α > 0 ,
where sign ( z ) denotes the sign function that returns 1 if z < 0 , + 1 if z > 0 and 0 otherwise. Moreover, by setting λ = 0 in (41) and recalling (40), for x 0 S one has
P Y ( S | x 0 ) = 0 + g Y ( S , t | x 0 ) d t = 1 , α < 0 , Erfc α σ x 0 + β α Erfc α σ S + β α , α > 0 ,
so that the first passage through the state S is a sure event for α < 0 .
The inverse LT of g λ ( Y ) ( S | x 0 ) can be obtained in closed form only if S = β / α .
Proposition 6.
For the TH-OU process, the FPT density through the boundary S = β / α is
g Y β α , t | x 0 = 2 e α t | x 0 + β / α | σ π α e 2 α t 1 3 / 2 exp α e 2 α t ( x 0 + β / α ) 2 σ 2 ( e 2 α t 1 ) , x 0 β / α ,
and the ultimate FPT probability is
P Y β α | x 0 = 0 + g Y β α , t | x 0 d t = 1 , α < 0 , Erfc α σ x 0 + β α , α > 0 .
Proof. 
Because
D ν ( 0 ) = 2 ν / 2 π Γ 1 ν 2 ,
from (41) for α 0 and x 0 β / α one has
g λ ( Y ) β α | x 0 = exp { | α | 2 σ 2 x 0 + β α 2 } 2 λ 2 | α | π Γ ( 1 2 + λ 2 | α | ) D λ | α | 2 | α | σ | x 0 + β α | , α < 0 , exp { α 2 σ 2 x 0 + β α 2 2 λ 2 α + 1 2 π Γ ( 1 + λ 2 α D λ α 1 2 α σ | x 0 + β α | , α > 0 .
Equation (43) follows by taking the inverse LT of (45) and making use of the following result (cf. Erdèlyi et al. [36], p. 290, no. 9):
0 + e p t e t / 2 ( e t 1 ) ν + 1 / 2 exp γ 2 4 ( e t 1 ) D 2 ν γ 1 e t d t = 2 p + ν Γ ( p + ν ) D 2 p ( γ ) , Re p > 0 .
Moreover, by setting λ = 0 in (45) and recalling (40), one obtains (44). □
When α < 0 , the FPT moments through S starting from x 0 can be evaluated by making use of Siegert Formulas (10) and (11) with r 1 = and r 2 = + . In particular, for n = 1 and α < 0 one has
t 1 ( Y ) ( S | x 0 ) = 1 | α | { π 2 Erfi | α | σ x 0 S + β α Erfi | α | σ x 0 S + β α + ψ 1 | α | σ S + β α ψ 1 | α | σ x 0 + β α } , x 0 S , ,
where
Erfi ( z ) = 2 π 0 z e u 2 d u = 2 π k = 0 + z 2 k + 1 ( 2 k + 1 ) k ! , ψ 1 ( z ) = k = 0 + 2 k z 2 k + 2 ( k + 1 ) ( 2 k + 1 ) ! ! ·
Furthermore, for α < 0 from (10) and (11) one obtains (cf. Ricciardi et al. [12])
lim S + t n ( Y ) ( S | x 0 ) [ t 1 ( Y ) ( S | x 0 ) ] n = n ! ( x 0 < S )
for n = 1 , 2 , so that for α < 0 the FPT density of the Ornstein–Uhlenbeck process exhibits an exponential asymptotic behavior as the boundary moves away from the starting point.
In Table 3 and Table 4, the mean t 1 ( Y ) ( S | x 0 ) , the variance Var ( Y ) ( S | x 0 ) , the coefficient of variation Cv ( Y ) ( S | x 0 ) , and the skewness Σ ( Y ) ( S | x 0 ) of the FPT, obtained by using (10) and (11), are listed for x 0 = 4 , α = 0.02 , σ = 1 and some choices of β and S.
From Table 3, we note that for the TH-OU process Y ( t ) the coefficient of variation approaches the value 1 and the skewness approaches the value 2 for large boundaries. Hence, when α < 0 the FPT density of the TH-OU process exhibits an exponential behavior for large boundaries S, such that S > x 0 .
From Table 4, we note that for the TH-OU process Y ( t ) the coefficient of variation and the skewness decreases as S decreases.
Moreover, taking the Laplace transform in (3) one has
a λ ( Y ) ( x | x 0 ) = f λ ( Y ) ( x | x 0 ) g λ ( Y ) ( 0 | x 0 ) f λ ( Y ) ( x | 0 ) , x 0 > 0 , x > 0 ,
so that, recalling (38) and (41), one can obtain the LT of a Y ( x , t | x 0 ) for the TH-OU process in ( 0 , + ) with 0 absorbing boundary.
Proposition 7.
Let Y ( t ) be a TH-OU process.
  • For α < 0 , one has
    L Y ( x | x 0 ) = 0 + a Y ( x , t | x 0 ) d t = 1 σ π | α | exp | α | σ 2 x + β α 2 × Erfi | α | σ x 0 x + β α Erfi | α | σ β α , x 0 > 0 , x > 0 .
  • For α > 0 , it results in
    L Y ( x | x 0 ) = 0 + a Y ( x , t | x 0 ) d t = 1 σ π α exp α σ 2 x + β α 2 Erfc α σ x 0 x + β α Erfc α σ β α × Erf α σ x 0 x + β α Erf α σ β α , x 0 > 0 , x > 0 .
Proof. 
Because + is a nonattracting boundary for α < 0 and attracting for α > 0 , Equations (47) and (48) follow from (12) making use of (37) and (42). □
From (47) and (48), for α 0 , β R and σ > 0 one obtains lim x 0 L Y ( x | x 0 ) = 0 and lim x + L Y ( x | x 0 ) = 0 .
In Figure 2, the asymptotic average of the local time for the TH-OU process Y ( t ) is plotted for x 0 = 4 , σ = 1 and some choices of α and β .

4. Feller-Type Diffusion Process

Let { Z ( t ) , t t 0 } , t 0 0 be a TNH-F process having infinitesimal drift and infinitesimal variance
C 1 ( x , t ) = α ( t ) x + β ( t ) , C 2 ( x , t ) = 2 r ( t ) x ,
with state space ( 0 , + ) , where α ( t ) R , β ( t ) R , r ( t ) > 0 continuous functions.
We point out that the processes (27) and (49) have identical infinitesimal drifts; instead, the infinitesimal variances are different in terms of the involved noise intensity functions.
The TNH-F process is used to describe the growth of a population (cf. Feller [37], Giorno and Nobile [38]) and the number of customers in queueing models (cf. Di Crescenzo and Nobile [39]). This process is also applied in mathematical finance to study stochastic volatility and interest rates (see Tian and Zhang [40], Cox et al. [41], Di Nardo and D’Onofrio [42]) and in neurobiology to model the input–output behavior of single neurons (see Ditlevsen and Lánský [43], D’Onofrio et al. [44]).
We consider the TNH-F process Z ( t ) , having infinitesimal moments (49), with an absorbing condition placed in the zero state and we denote with a Z ( x , t | x 0 , t 0 ) its PDF. We assume that α ( t ) R , β ( t ) R , r ( t ) > 0 , β ( t ) ξ r ( t ) , with 0 ξ < 1 . For the TNH-F process Z ( t ) with an absorbing boundary in zero, we consider two cases: the proportional case in which β ( t ) = ξ r ( t ) , with 0 ξ < 1 and r ( t ) > 0 , and the time-homogeneous case.

4.1. Proportional Case for the Feller Process

We assume that α ( t ) R , r ( t ) > 0 and β ( t ) = ξ r ( t ) , with 0 ξ < 1 , in (49). As proven in Giorno and Nobile [45] one has
a Z ( x , t | x 0 , t 0 ) = e A ( t | t 0 ) Γ ( 2 ξ ) 1 R ( t | t 0 ) 2 ξ x 0 1 ξ exp x 0 R ( t | t 0 ) , x = 0 , e A ( t | t 0 ) R ( t | t 0 ) x 0 x ( 1 ξ ) / 2 exp x 0 + x e A ( t | t 0 ) R ( t | t 0 ) × exp 1 ξ 2 A ( t | t 0 ) I 1 ξ 2 x x 0 e A ( t | t 0 ) R ( t | t 0 ) , x > 0 ,
with A ( t | t 0 ) given in (30), I ν ( z ) defined in (24) and
R ( t | t 0 ) = t 0 t r ( θ ) e A ( θ | t 0 ) d θ .
Proposition 8.
Let α ( t ) R , r ( t ) > 0 and β ( t ) = ξ r ( t ) , with 0 ξ < 1 , in (49). For the TNH-F process Z ( t ) one has
g Z ( 0 , t | x 0 , t 0 ) = 1 Γ ( 1 ξ ) r ( t ) e A ( t | t 0 ) R ( t | t 0 ) x 0 R ( t | t 0 ) 1 ξ exp x 0 R ( t | t 0 ) , x 0 > 0 ,
with R ( t | t 0 ) given in (51). Moreover, it results in
P Z ( 0 | x 0 , t 0 ) = t 0 + g Z ( 0 , t | x 0 , t 0 ) d t = 1 , lim t + R ( t | t 0 ) = + , Γ 1 ξ , x 0 c Γ 1 ξ , lim t + R ( t | t 0 ) = c .
Proof. 
From (50), one has (cf. Erdèlyi et al. [36], p. 197, no. 19)
0 + a Z ( x , t | x 0 , t 0 ) d x = 1 Γ ( 1 ξ ) γ 1 ξ , x 0 R ( t | t 0 ) , 0 ξ < 1 ,
where Γ ( ν ) is the Euler gamma function and γ ( ν , z ) = 0 z y ν 1 e y d y , with ν > 0 , is the incomplete gamma function. Hence, due to (5) and recalling (54), Equation (52) follows. Finally, Equation (53) is obtained, making use of (54) in (6). □
We point out that the general TNH-F process with an absorbing boundary in zero is considered in Giorno and Nobile [45], Masoliver and Perelló [46], Masoliver [47] and Lavigne and Roques [48].

4.2. Time-Homogeneous Case for the Feller Process

Let Z ( t ) be the TH-F process, obtained by setting α ( t ) = α , β ( t ) = β and r ( t ) = r in (49). From (9), the scale function and the speed density of the TH-F process Z ( t ) are
h Z ( x ) = x β / r exp α x r , s Z ( x ) = x β / r 1 r exp α x r ,
respectively. As proven by Feller, the state 0 is an exit boundary for β 0 , regular for 0 < β < r and entrance for β 0 . Furthermore, the end point + is a nonattracting natural boundary for α 0 and an attracting natural boundary for α > 0 . In the sequel, we assume that α R , β R , r > 0 , with β < r , and an absorbing condition is set in the zero-state.
As proven in Giorno and Nobile [45], for a TH-F process Z ( t ) having β R , r > 0 , with β < r , one has
  • If α = 0 one has
    a Z ( x , t | x 0 ) = 1 Γ ( 2 β / r ) 1 r t 2 β / r x 0 1 β / r exp x 0 r t , x = 0 , 1 r t x 0 x ( 1 β / r ) / 2 exp x + x 0 r t I 1 β / r 2 x x 0 r t , x > 0 .
  • If α 0 one obtains
    a Z ( x , t | x 0 ) = e α t Γ ( 2 β / r ) α e α t r ( e α t 1 ) 2 β / r x 0 1 β / r exp α x 0 e α t r ( e α t 1 ) , x = 0 , α e α ( 1 β / r ) t / 2 r ( e α t 1 ) x 0 x ( 1 β / r ) / 2 exp α ( x + x 0 e α t ) r ( e α t 1 ) I 1 β / r 2 α x x 0 e α t r ( e α t 1 ) , x > 0 .
Proposition 9.
Let α R , β R , r > 0 , with β < r . For the TH-F process Z ( t ) , with x 0 > 0 , one has
g Z ( 0 , t | x 0 ) = 1 t Γ ( 1 β / r ) x 0 r t 1 β / r exp x 0 r t , α = 0 , 1 Γ ( 1 β / r ) α e α t 1 α x 0 e α t r ( e α t 1 ) 1 β / r exp α x 0 e α t r ( e α t 1 ) , α 0
and
P Z ( 0 | x 0 ) = 0 + g Z ( 0 , t | x 0 ) d t = 1 , α 0 , Γ 1 β r , α x 0 r Γ 1 β / r , α > 0 .
Proof. 
From (56) and (57), one obtains (cf. Erdèlyi et al. [36], p. 197, no. 19)
0 + a Z ( x , t | x 0 ) d x = 1 Γ ( 1 β / r ) γ 1 β r , x 0 r t , α = 0 , 1 Γ ( 1 β / r ) γ 1 β r , α x 0 e α t r ( e α t 1 ) , α 0 .
Making use of (60) in (5), Equation (58) follows. Finally, by virtue of (6) and (60), we obtain the FPT probability (59). □
By applying the Siegert Formula (11) with r 2 = + and recalling (55), for α = 0 and β < r one has that the FPT mean t 1 ( Z ) ( 0 | x 0 ) diverges, whereas for α < 0 and β < r one obtains
t 1 ( Z ) ( S | x 0 ) = 1 | α | | α | S / r | α | x 0 / r x β / r e x Γ β r , x d x , x 0 > S 0 .
In Table 5, the mean t 1 ( Z ) ( S | x 0 ) , the variance Var ( Z ) ( S | x 0 ) , the coefficient of variation Cv ( Z ) ( S | x 0 ) , and the skewness Σ ( Z ) ( S | x 0 ) of the FPT, obtained by using the Siegert Formula (11), are listed for x 0 = 4 and some choices of S, with α = 0.02 , β = 0.1 , 0 , 0.1 and r = 0.5 .
As shown in Table 5, for the TH-F process Z ( t ) the mean and the variance of the FPT increases as S decreases; instead, the coefficient of variation and the skewness decrease as S decreases.
Proposition 10.
Let Z ( t ) be a TH-F process having β R , r > 0 , with β < r .
  • If α 0 , for x 0 > 0 and x > 0 one has
    L Z ( x | x 0 ) = 0 + a Z ( x , t | x 0 ) d t = 1 r r | α | x 1 β / r e | α | x / r 0 | α | ( x 0 x ) / r y β / r e y d y , α < 0 1 r 1 1 β / r x 0 x x 1 β / r , α = 0 .
  • If α > 0 , for x 0 > 0 and x > 0 one obtains
    L Z ( x | x 0 ) = 0 + a Z ( x , t | x 0 ) d t = 1 r r α x 1 β / r e α x / r γ 1 β r , α ( x 0 x ) r Γ 1 β r , α ( x 0 x ) r Γ 1 β r ·
Proof. 
Because + is a nonattracting boundary for α 0 and attracting for α > 0 , Equations (61) and (62) follow from (12), making use of (55) and (59). □
From (61) and (62), for α R , β R , r > 0 , with β < r , one has lim x + L Z ( x | x 0 ) = 0 and
lim x 0 L Z ( x | x 0 ) = 1 r 1 1 β / r α 0 , 1 r 1 1 β / r Γ 1 β r , α x 0 r Γ 1 β r , α > 0 .
Therefore, for the TH-F process the asymptotic average of local time tend to zero as x increases, whereas it is positive for x 0 .
In Figure 3, the asymptotic average of the local time for the TH-F process Z ( t ) is plotted for x 0 = 4 , r = 0.5 and some choices of α and β .

5. Relationships and Asymptotic Results

In this section, for β ( t ) = r ( t ) / 2 some relationships between the PDF in the presence of an absorbing boundary in the zero state and between the FPT densities through zero for Wiener, Ornstein–Uhlenbeck and Feller processes are proven; moreover, for β ( t ) = ξ r ( t ) ( 0 < ξ < 1 ) some asymptotic results for large times between the FPT densities are provided.

5.1. Relations between the Transition Densities with an Absorbing Boundary in the Zero State

We consider the TNH-F process (49) with β ( t ) = r ( t ) / 2 in the presence of an absorbing boundary in the zero state, and we show that its PDF can be related to the PDF of the Wiener and of the Ornstein–Uhlenbeck processes with an absorbing boundary in the zero state.
Proposition 11.
Let Z ( t ) be a TNH-F process with C 1 ( t ) = r ( t ) / 2 and C 2 ( x , t ) = 2 r ( t ) x , where r ( t ) > 0 , and let X ( t ) be a TNH-W process with A 1 = 0 and A 2 ( t ) = r ( t ) / 2 . One has
a Z ( x , t | x 0 , t 0 ) = 1 2 x a X ( x , t | x 0 , t 0 ) , x 0 > 0 , x > 0 ,
g Z ( 0 , t | x 0 , t 0 ) = g X ( 0 , t | x 0 , t 0 ) , x 0 > 0 .
Proof. 
For the TNH-F process Z ( t ) , by setting α ( t ) = 0 and β ( t ) = r ( t ) / 2 in (50) and in Proposition 9, recalling that
I 1 / 2 ( x ) = 2 π sinh ( x ) x , γ 1 2 , x = π Erf ( x ) ,
one has
a Z ( x , t | x 0 , t 0 ) = 2 x 0 π [ R ˜ ( t | t 0 ) ] 3 exp x 0 R ˜ ( t | t 0 ) , x = 0 , 1 2 π R ˜ ( t | t 0 ) x exp x x 0 2 R ˜ ( t | t 0 ) exp x + x 0 2 R ˜ ( t | t 0 ) , x > 0 ,
and
g Z ( 0 , t | x 0 , t 0 ) = r ( t ) x 0 π [ R ˜ ( t | t 0 ) ] 3 exp x 0 R ˜ ( t | t 0 ) , x 0 > 0 ,
where R ˜ ( t | t 0 ) = t 0 t r ( θ ) d θ . Furthermore, for the TNH-W process X ( t ) with β ( t ) = 0 and σ 2 ( t ) = r ( t ) / 2 , one has V X ( t | t 0 ) = R ˜ ( t | t 0 ) / 2 . Then, (63) and (64) follow by comparing (66) and (67) with (17) and (18), respectively. □
Under the assumptions of the Proposition 13, one has L Z ( x | x 0 , t 0 ) = L X ( x | x 0 , t 0 ) / ( 2 x ) for x > 0 , x 0 > 0 and, if lim t + R ˜ ( t | t 0 ) = + , one obtains t n ( Z ) ( 0 | x 0 , t 0 ) = t n ( X ) ( 0 | x 0 , t 0 ) for n = 1 , 2 , with x 0 > 0 .
Proposition 12.
Let Z ( t ) be a TNH-F process with C 1 ( x , t ) = α ( t ) x + r ( t ) / 2 and C 2 ( x , t ) = 2 r ( t ) x , where α ( t ) is not always zero and r ( t ) > 0 , and let Y ( t ) be a TNH-OU process with B 1 ( x , t ) = α ( t ) x / 2 and B 2 ( t ) = r ( t ) / 2 . One has
a Z ( x , t | x 0 , t 0 ) = 1 2 x a Y ( x , t | x 0 , t 0 ) , x 0 > 0 , x > 0 ,
g Z ( 0 , t | x 0 , t 0 ) = g Y ( 0 , t | x 0 , t 0 ) , x 0 > 0 .
Proof. 
For the TNH-F process Z ( t ) , by setting β ( t ) = r ( t ) / 2 in (50) and in Proposition 9, recalling (65), one obtains
a Z ( x , t | x 0 , t 0 ) = 2 e A ( t | t 0 ) x 0 π [ R ( t | t 0 ) ] 3 exp x 0 R ( t | t 0 ) , x = 0 , e A ( t | t 0 ) / 2 2 π R ( t | t 0 ) x [ exp x e A ( t | t 0 ) x 0 2 R ( t | t 0 ) exp x e A ( t | t 0 ) + x 0 2 R ( t | t 0 ) ] , x > 0 ,
and
g Z ( 0 , t | x 0 , t 0 ) = r ( t ) e A ( t | t 0 ) x 0 π [ R ( t | t 0 ) ] 3 exp x 0 R ( t | t 0 ) , x 0 > 0 ,
with A ( t | t 0 ) and R ( t | t 0 ) given in (30) and (51), respectively. Moreover, in the TNH-OU process Y ( t ) we set β ( t ) = 0 , σ 2 ( t ) = r ( t ) / 2 and we change α ( t ) into α ( t ) / 2 , so that, by virtue of (29) and (51), one has V Y ( t | t 0 ) = R ( t | t 0 ) e A ( t | t 0 ) / 2 . Then, (68) and (69) follow by comparing (70) and (71) with (33) and (34), respectively. □
Under the assumptions of Proposition 14, one has L Z ( x | x 0 , t 0 ) = L Y ( x | x 0 , t 0 ) /
( 2 x ) for x > 0 , x 0 > 0 and, if lim t + R ( t | t 0 ) = + , one obtains t n ( Z ) ( 0 | x 0 , t 0 ) = t n ( Y ) ( 0 | x 0 , t 0 ) for n = 1 , 2 , with x 0 > 0 .

5.2. Asymptotic Behaviors between the FPT Densities

In this section, for β ( t ) = ξ r ( t ) ( 0 < ξ < 1 ) some asymptotic results for large times between the FPT densities of TNH-W, TNH-OU and TNH-F processes are shown.
Proposition 13.
Let Z ( t ) be a TNH-F process with C 1 ( t ) = ξ r ( t ) and C 2 ( x , t ) = 2 r ( t ) x , where r ( t ) > 0 , 0 < ξ < 1 , and let X ( t ) be a TNH-W process with A 1 = 0 and A 2 ( t ) = ξ r ( t ) . If lim t + R ˜ ( t | t 0 ) = + , and one has
lim t + g Z ( 0 , t | x 0 , t 0 ) [ R ˜ ( t | t 0 ) ] 1 / 2 ξ g X ( 0 , t | x 0 1 ξ , t 0 ) = 2 π ξ Γ ( 1 ξ ) , x 0 > 0 .
Proof. 
Recalling (18) and (67) and noting that V X ( t | t 0 ) = ξ R ˜ ( t | t 0 ) , one has
g Z ( 0 , t | x 0 , t 0 ) [ R ˜ ( t | t 0 ) ] 1 / 2 ξ g X ( 0 , t | x 0 1 ξ , t 0 ) = 2 π ξ Γ ( 1 ξ ) exp x 0 R ˜ ( t | t 0 ) + x 0 2 ( 1 ξ ) 2 ξ R ˜ ( t | t 0 ) ,
from which, under the assumption lim t + R ˜ ( t | t 0 ) = + , Equation (72) follows. □
Proposition 14.
Let Z ( t ) be a TNH-F process having C 1 ( x , t ) = α ( t ) x + ξ r ( t ) and C 2 ( x , t ) = 2 r ( t ) x , with α ( t ) not always zero, r ( t ) > 0 , 0 < ξ < 1 , and let Y ( t ) be a TNH-OU process with B 1 ( x , t ) = α ( t ) x / 2 and B 2 ( t ) = ξ r ( t ) . If lim t + R ( t | t 0 ) = + , and one has
lim t + g Z ( 0 , t | x 0 , t 0 ) [ R ( t | t 0 ) ] 1 / 2 ξ g Y ( 0 , t | x 0 1 ξ , t 0 ) = 2 π ξ Γ ( 1 ξ ) , x 0 > 0 .
Proof. 
Making use of (34) and (52) and noting that V Y ( t | t 0 ) = ξ R ( t | t 0 ) e A ( t | t 0 ) , one obtains
g Z ( 0 , t | x 0 , t 0 ) [ R ( t | t 0 ) ] 1 / 2 ξ g Y ( 0 , t | x 0 1 ξ , t 0 ) = 2 π ξ Γ ( 1 ξ ) exp x 0 R ( t | t 0 ) + x 0 2 ( 1 ξ ) 2 ξ R ( t | t 0 ) ,
from which, recalling that lim t + R ( t | t 0 ) = + , Equation (73) follows. □

6. Conclusions

For the Wiener, Ornstein–Uhlenbeck, and Feller processes, we analyze the transition densities in the presence of an absorbing boundary in the zero state and the FPT problem to the zero state. Particular attention is dedicated to the proportional cases and to the time-homogeneous cases, by achieving the FPT densities through the zero state. Extensive computation are performed with MATHEMATICA to obtain the mean, the variance, the coefficient of variation and the skewness of FPT for TH-W, TH-OU and TH-F processes. Moreover, for these processes, a detailed study of the asymptotic average of local time with an absorbing boundary in the zero-state is carried out.
In Table 6, a summary containing the conditions and the most important equations numbering in Section 2, Section 3 and Section 4 concerning the absorbing problem for Wiener, Ornstein–Uhlenbeck and Feller diffusion processes is given.
As shown in Table 6, by setting β ( t ) = 0 in TNH-W, TNH-OU and TNH-F processes, the PDF in the presence of an absorbing boundary in the zero state and the FPT density through zero are given in closed form. Moreover, in TH-W, TH-F processes, the previous densities are obtainable, whereas for the TH-OU process only the LT is available.
The knowledge of the PDF in the presence of an absorbing boundary in the zero state is of interest in the context of biological systems because it allows us to evaluate the survival probabilities (20), (36) and (54) for Wiener, Ornstein–Uhlenbeck and Feller processes, respectively. Moreover, such PDF allows one to get information on the FPT density through zero (extinction density) (18), (34), and (52) and on the probability of extinction (19), (35) and (53) of the considered processes. Furthermore, the asymptotic average of the local time for TH-W, TH-OU, and TH-F processes provides information on the average of the sojourn time in the various states before the absorption occurs in the zero state.
The results of Section 5 show that the same FPT density through the zero-state (extinction density) may correspond to different diffusion processes with modified initial states.

Author Contributions

Conceptualization, V.G. and A.G.N.; methodology, V.G. and A.G.N.; validation, V.G. and A.G.N.; formal analysis, V.G. and A.G.N.; investigation, V.G. and A.G.N.; supervision, V.G. and A.G.N. All authors have read and agreed to the published version of the manuscript.

Funding

This research is partially supported by MIUR - PRIN 2017, Project “Stochastic Models for Complex Systems”.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The authors are members of the research group GNCS of INdAM.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
PDFTransition Probability Density Function
FPTFirst Passage Time
TNH-DTime Inhomogeneous Diffusion
TNH-WTime Inhomogeneous Wiener
TNH-OUTime Inhomogeneous Ornstein-Uhlenbeck
TNH-FTime Inhomogeneous Feller
TH-DTime Homogeneous Diffusion
TH-WTime Homogeneous Wiener
TH-OUTime Homogeneous Ornstein-Uhlenbeck
TH-FTime Homogeneous Feller

References

  1. Bailey, N.T.J. The Elements of Stochastic Processes with Applications to the Natural Sciences; John Wiley & Sons, Inc.: New York, NY, USA, 1964. [Google Scholar]
  2. Ricciardi, L.M. Diffusion Processes and Related Topics in Biology; Springer: Berlin/Heidelberg, Germany, 1977. [Google Scholar]
  3. Gardiner, C.W. Handbook of Stochastic Methods: For Physics, Chemistry and the Natural Sciences; Springer: Berlin/Heidelberg, Germany, 1985. [Google Scholar]
  4. Stirzaker, D. Stochastic Processes and Models; Oxford University Press: Oxford, UK, 2005. [Google Scholar]
  5. Janssen, J.; Manca, O.; Manca, R. Applied Diffusion Processes from Engineering to Finance; Wiley-ISTE: London, UK, 2013. [Google Scholar]
  6. Pavliotis, G.A. Stochastic Processes and Applications. Diffusion Processes, the Fokker-Planck and Langevin Equations; Springer: New York, NY, USA, 2014. [Google Scholar]
  7. Dynkin, E.B. Kolmogorov and the theory of Markov processes. Ann. Probab. 1989, 17, 822–832. [Google Scholar] [CrossRef]
  8. Cox, D.R.; Miller, H.D. The Theory of Stochastic Processes; Chapman & Hall/CRC: Boca Raton, FL, USA, 1996. [Google Scholar]
  9. Allen, L.J.S. An Introduction to Stochastic Processes with Applications to Biology, 2nd ed.; Chapman and Hall/CRC: Lubbock, TX, USA, 2010. [Google Scholar]
  10. Allen, L.J.S. Stochastic Population and Epidemic Models Persistence and Extinction; Springer: Cham, Switzerland, 2015. [Google Scholar]
  11. Blake, I.F.; Lindsey, W.C. Level-Crossing Problems for Random Processes. IEEE Trans. Inf. Theory 1973, 19, 295–315. [Google Scholar] [CrossRef]
  12. Ricciardi, L.M.; Di Crescenzo, A.; Giorno, V.; Nobile, A.G. An outline of theoretical and algorithmic approaches to first passage time problems with applications to biological modeling. Math. Jpn. 1999, 50, 247–322. [Google Scholar]
  13. Ding, M.; Rangarajan, G. First Passage Time Problem: A Fokker-Planck Approach. In New Directions in Statistical Physics; Wille, L.T., Ed.; Springer: Berlin/Heidelberg, Germany, 2004; pp. 31–46. [Google Scholar]
  14. Molini, A.; Talkner, P.; Katul, G.G.; Porporato, A. First passage time statistics of Brownian motion with purely time dependent drift and diffusion. Physica A 2011, 390, 1841–1852. [Google Scholar] [CrossRef] [Green Version]
  15. Giorno, V.; Nobile, A.G. On the first-passage time problem for a Feller-type diffusion process. Mathematics 2021, 9, 2470. [Google Scholar] [CrossRef]
  16. Masoliver, J. Random Processes: First-Passage and Escape; World Scientific: Singapore, 2018. [Google Scholar]
  17. Linetsky, V. Computing hitting time densities for CIR and OU diffusions. Applications to mean-reverting models. J. Comput. Finance 2004, 7, 1–22. [Google Scholar] [CrossRef] [Green Version]
  18. Guo, G.; Chen, B.; Zhao, X.; Zhao, F.; Wang, Q. First passage time distribution of a modified fractional diffusion equation in the semi-infinite interval. Physica A 2015, 433, 279–290. [Google Scholar] [CrossRef]
  19. Wiese, K.J. First passage in an interval for fractional Brownian motion. Phys. Rev. E 2019, 99, 032106. [Google Scholar] [CrossRef] [Green Version]
  20. Abundo, M. An inverse first-passage problem revisited: The case of fractional Brownian motion, and time-changed Brownian motion. Stoch. Anal. Appl. 2019, 37, 708–716. [Google Scholar] [CrossRef]
  21. Leonenko, N.; Pirozzi, E. First passage times for some classes of fractional time-changed diffusions. Stoch. Anal. Appl. 2022, 40, 735–763. [Google Scholar] [CrossRef]
  22. Siegert, A.J.F. On the first passage time probability problem. Phys. Rev. 1951, 81, 617–623. [Google Scholar] [CrossRef]
  23. Karlin, S.; Taylor, H.W. A Second Course in Stochastic Processes; Academic Press: New York, NY, USA, 1981. [Google Scholar]
  24. Aït-Sahalia, Y.; Park, J.Y. Stationarity-based specification tests for diffusions when the process is nonstationary. J. Econom. 2012, 169, 279–292. [Google Scholar] [CrossRef]
  25. Feller, W. Two singular diffusion problems. Ann. Math. 1951, 5, 173–182. [Google Scholar] [CrossRef]
  26. Feller, W. The parabolic differential equations and the associated semi-groups of transformations. Ann. Math. 1952, 55, 468–518. [Google Scholar] [CrossRef]
  27. Giorno, V.; Nobile, A.G. Asymptotic average of the local time for one-dimensional diffusion processes. Lect. Notes Semin. Interdiscip. Mat. 2015, 12, 139–160. [Google Scholar]
  28. Knight, F.B. Essentials of Brownian Motion and Diffusions; American Mathematical Society: Providence, RI, USA, 1981. [Google Scholar]
  29. Abramowitz, I.A.; Stegun, M. Handbook of Mathematical Functions; Dover Publications: Mineola, NY, USA, 1972. [Google Scholar]
  30. Uhlenbeck, G.E.; Ornstein, L.S. On the theory of the Brownian motion. Phys. Rev. 1930, 36, 823–841. [Google Scholar] [CrossRef]
  31. Ricciardi, L.M.; Sacerdote, L. The Ornstein-Uhlenbeck process as a model for neuronal activity. I. Mean and variance of the firing time. Biol. Cybern. 1979, 35, 1–9. [Google Scholar] [CrossRef]
  32. Lánský, P.; Ditlevsen, S. A review of the methods for signal estimation in stochastic diffusion leaky integrate-and-fire neuronal models. Biol. Cybern. 2008, 99, 253–262. [Google Scholar] [CrossRef]
  33. Vasicek, O. An equilibrium characterization of the term structure. J. Financ. Econ. 1977, 5, 177–188. [Google Scholar] [CrossRef]
  34. Hull, J.; White, A. Pricing interest-rate-derivative securities. Rev. Financ. Stud. 1990, 3, 573–592. [Google Scholar] [CrossRef] [Green Version]
  35. Gradshteyn, I.S.; Ryzhik, I.M. Table of Integrals, Series and Products; Academic Press Inc.: New York, NY, USA, 2014. [Google Scholar]
  36. Erdèlyi, A.; Magnus, W.; Oberthettinger, F.; Tricomi, F.G. Tables of Integral Transforms; Mc Graw-Hill: New York, NY, USA, 1954; Volume 1. [Google Scholar]
  37. Feller, W. Diffusion processes in genetics. In Proceedings of the Second Berkeley Symposium on Mathematical Statistics and Probability, Berkeley, CA, USA, 31 July–12 August 1950; pp. 227–246. [Google Scholar]
  38. Giorno, V.; Nobile, A.G. Time-inhomogeneous Feller-type diffusion process in population dynamics. Mathematics 2021, 9, 1879. [Google Scholar] [CrossRef]
  39. Di Crescenzo, A.; Nobile, A.G. Diffusion approximation to a queueing system with time-dependent arrival and service rates. Queueing Syst. 1995, 19, 41–62. [Google Scholar] [CrossRef]
  40. Tian, Y.; Zhang, H. Skew CIR process, conditional characteristic function, moments and bond pricing. Appl. Math. Comput. 2018, 329, 230–238. [Google Scholar] [CrossRef]
  41. Cox, J.C.; Ingersoll, J.E., Jr.; Ross, S.A. A theory of the term structure of interest rates. Econometrica 1985, 53, 385–407. [Google Scholar] [CrossRef]
  42. Di Nardo, E.; D’Onofrio, G. A cumulant approach for the first-passage-time problem of the Feller square-root process. Appl. Math. Comput. 2021, 391, 125707. [Google Scholar] [CrossRef]
  43. Ditlevsen, S.; Lánský, P. Estimation of the input parameters in the Feller neuronal model. Phys. Rev. E 2006, 73, 061910. [Google Scholar] [CrossRef]
  44. D’Onofrio, G.; Lánský, P.; Pirozzi, E. On two diffusion neuronal models with multiplicative noise: The mean first-passage time properties. Chaos 2018, 28, 043103. [Google Scholar] [CrossRef]
  45. Giorno, V.; Nobile, A.G. Time-inhomogeneous Feller-type diffusion process with absorbing boundary condition. J. Stat. Phys. 2021, 183, 1–27. [Google Scholar] [CrossRef]
  46. Masoliver, J.; Perelló, J. First-passage and escape problems in the Feller process. Phys. Rev. E 2012, 86, 041116. [Google Scholar] [CrossRef] [Green Version]
  47. Masoliver, J. Extreme values and the level-crossing problem: An application to the Feller process. Phys. Rev. E 2014, 89, 042106. [Google Scholar] [CrossRef] [Green Version]
  48. Lavigne, F.; Roques, L. Extinction times of an inhomogeneous Feller diffusion process: A PDF approach. Expo. Math. 2021, 39, 137–142. [Google Scholar] [CrossRef]
Figure 1. L X ( x | x 0 ) , given in (26), with x 0 = 4 , σ = 1 and some choices of β . In (a) β = 0.1 , 0 , 0.1 and in (b) β = 0.2 , 0 , 0.2 .
Figure 1. L X ( x | x 0 ) , given in (26), with x 0 = 4 , σ = 1 and some choices of β . In (a) β = 0.1 , 0 , 0.1 and in (b) β = 0.2 , 0 , 0.2 .
Fractalfract 07 00011 g001
Figure 2. L Y ( x | x 0 ) , given in Proposition 7, with x 0 = 4 , σ = 1 , and some choices of β . In (a) α = 0.02 and in (b) α = 0.02 .
Figure 2. L Y ( x | x 0 ) , given in Proposition 7, with x 0 = 4 , σ = 1 , and some choices of β . In (a) α = 0.02 and in (b) α = 0.02 .
Fractalfract 07 00011 g002
Figure 3. L Z ( x | x 0 ) , given in Proposition 12, with x 0 = 4 , r = 0.5 , and some choices of β . In (a) α = 0.02 , in (b) α = 0.02 and in (c) α = 0 .
Figure 3. L Z ( x | x 0 ) , given in Proposition 12, with x 0 = 4 , r = 0.5 , and some choices of β . In (a) α = 0.02 , in (b) α = 0.02 and in (c) α = 0 .
Fractalfract 07 00011 g003aFractalfract 07 00011 g003b
Table 1. For the Wiener process, with A 1 ( x ) = β and A 2 ( x ) = 1 , the mean, the variance, the coefficient of variation, and the skewness of FPT are listed for x 0 = 4 , β = 0.1 , 0.2 and for increasing values the boundary S > x 0 .
Table 1. For the Wiener process, with A 1 ( x ) = β and A 2 ( x ) = 1 , the mean, the variance, the coefficient of variation, and the skewness of FPT are listed for x 0 = 4 , β = 0.1 , 0.2 and for increasing values the boundary S > x 0 .
S t 1 ( X ) ( S | x 0 ) Var ( X ) ( S | x 0 ) Cv ( X ) ( S | x 0 ) Σ ( X ) ( S | x 0 )
β = 0.1 10096096,0000.3227490.968246
5004960496,0000.1419900.425971
10009960996,0000.1002010.300602
150014,9601,496,0000.08175870.245276
200019,9601,996,0000.07078150.212344
250024,9602,496,0000.06329620.189889
300029,9602,996,0000.05777360.173321
β = 0.2 10048012,0000.2282180.684653
500248062,0000.1004020.301207
10004980124,5000.07085250.212558
15007480187,0000.05781220.173436
20009980249,5000.05005010.150150
250012,480312,0000.04475720.134272
300014,980374,5000.04085210.122556
Table 2. As in Table 1, with x 0 = 4 , σ = 1 , β = 0.1 , 0.2 and for decreasing values the boundary S [ 0 , x 0 ) .
Table 2. As in Table 1, with x 0 = 4 , σ = 1 , β = 0.1 , 0.2 and for decreasing values the boundary S [ 0 , x 0 ) .
S t 1 ( X ) ( S | x 0 ) Var ( X ) ( S | x 0 ) Cv ( X ) ( S | x 0 ) Σ ( X ) ( S | x 0 )
β = 0.1 3.555004.4721413.4164
3.01010003.162289.48683
2.51515002.581997.74597
2.02020002.236076.7082
1.52525002.06.0
1.03030001.825745.47723
0.53535001.690315.07093
0.04040001.581144.74342
β = 0.2 3.52.562.53.162289.48683
3.051252.236076.7082
2.57.5187.51.825745.47723
2.0102501.581144.74342
1.512.5312.51.414214.24264
1.0153751.290993.87298
0.517.5437.51.195233.58569
0.0205001.118033.3541
Table 3. For the TH-OU process, with B 1 ( x ) = 0.02 x + β and B 2 ( x ) = 1 , the mean, the variance, the coefficient of variation, and the skewness of FPT are listed for x 0 = 4 , β = 0.1 , 0 , 0.1 and for increasing values the boundary S > x 0 .
Table 3. For the TH-OU process, with B 1 ( x ) = 0.02 x + β and B 2 ( x ) = 1 , the mean, the variance, the coefficient of variation, and the skewness of FPT are listed for x 0 = 4 , β = 0.1 , 0 , 0.1 and for increasing values the boundary S > x 0 .
S t 1 ( Y ) ( S | x 0 ) Var ( Y ) ( S | x 0 ) Cv ( Y ) ( S | x 0 ) Σ ( Y ) ( S | x 0 )
β = 0.1 5 1.491996 × 10 2 1.265053 × 10 5 2.3838933.759894
10 3.97436 × 10 3 1.842427 × 10 7 1.0800102.020738
15 1.005474 × 10 5 1.017251 × 10 10 1.0030971.999075
20 7.036678 × 10 6 4.951899 × 10 13 1.0000421.986165
25 1.413375 × 10 9 1.997625 × 10 18 0.99999921.970391
β = 0 5 3.077237 × 10 1 5.966963 × 10 3 2.5102434.446831
10 4.475225 × 10 2 2.545499 × 10 5 1.1273832.098001
15 4.272683 × 10 3 1.855231 × 10 7 1.0080882.000674
20 1.008457 × 10 5 1.017263 × 10 10 1.0001361.993045
25 7.036975 × 10 6 4.951901 × 10 13 1.0000001.980451
β = 0.1 5 1.161050 × 10 1 7.316673 × 10 2 2.3297324.818854
10 1.162808 × 10 2 1.533674 × 10 4 1.0650212.181767
15 5.330310 × 10 2 2.639197 × 10 5 0.96379232.002620
20 4.358192 × 10 3 1.856168 × 10 7 0.98855831.997674
25 1.009312 × 10 5 1.017264 × 10 10 0.99928941.998622
Table 4. As in Table 3 with x 0 = 4 , α = 0.02 , σ = 1 , β = 0.1 , 0 , 0.1 and for decreasing values the boundary S [ 0 , x 0 ) .
Table 4. As in Table 3 with x 0 = 4 , α = 0.02 , σ = 1 , β = 0.1 , 0 , 0.1 and for decreasing values the boundary S [ 0 , x 0 ) .
S t 1 ( Y ) ( S | x 0 ) Var ( Y ) ( S | x 0 ) Cv ( Y ) ( S | x 0 ) Σ ( Y ) ( S | x 0 )
β = 0.1 3.52.32181938.350142.6671986.633648
3.04.74092281.061321.8990834.708105
2.57.264863128.77041.5619983.861617
2.09.902021182.22691.3632723.362382
1.512.66172242.31661.2294163.026467
1.015.55435310.09041.1321192.782932
0.518.59156386.80001.0578582.597858
0.021.7864473.94320.99925822.452724
β = 0 3.53.763743130.82653.038986.851246
3.07.756523281.91042.1646544.855690
2.512.00269457.19131.7814373.978597
2.016.53010661.52331.5559553.462116
1.521.37079900.92131.4045033.115650
1.026.561691182.8821.2948362.865643
0.532.145581516.8041.2115562.676884
0.038.172191914.5491.1462682.530113
β = 0.1 3.57.745813697.73323.4101836.676430
3.016.236151550.6032.4253104.719137
2.525.584272600.6851.9932913.859455
2.035.925673903.4241.7390743.354977
1.547.423275532.5041.5684453.018737
1.060.274067586.6971.4450932.778458
0.574.7175110199.601.3516662.599416
0.091.0464313553.661.2786902.462534
Table 5. For the TH-F process, with C 1 ( x ) = 0.02 x + β and C 2 ( x ) = x , the mean, the variance, the coefficient of variation, and the skewness of FPT are listed for x 0 = 4 , β = 0.1 , 0 , 0.1 and for decreasing values the boundary S [ 0 , x 0 ) .
Table 5. For the TH-F process, with C 1 ( x ) = 0.02 x + β and C 2 ( x ) = x , the mean, the variance, the coefficient of variation, and the skewness of FPT are listed for x 0 = 4 , β = 0.1 , 0 , 0.1 and for decreasing values the boundary S [ 0 , x 0 ) .
S t 1 ( Z ) ( S | x 0 ) Var ( Z ) ( S | x 0 ) Cv ( Z ) ( S | x 0 ) Σ ( Z ) ( S | x 0 )
β = 0.1 3.51.39262043.551364.73879913.55745
3.02.85900889.019223.3000969.501619
2.54.412074136.55212.6485397.684201
2.06.069704186.29772.2487246.585614
1.57.858374238.37331.9646995.823477
1.09.821571292.76901.7421345.249299
0.512.04587348.94651.5507474.792382
0.014.86611401.34131.3475964.429951
β = 0 3.51.70212169.834334.90957813.0398
3.03.512494144.38783.4209739.119134
2.55.452850224.42062.7473117.356073
2.07.554206310.94092.3342656.285290
1.59.864220405.35532.0410575.536382
1.012.46511509.7821.8113244.964136
0.515.53300627.87811.6131784.494585
0.019.91651768.91711.3922804.068225
β = 0.1 3.52.126455115.10395.04532112.41484
3.04.415487240.99263.5157938.659792
2.56.904328380.06982.8236466.964637
2.09.647790535.80212.3992445.929435
1.512.73337713.57522.0978615.199456
1.016.31983922.79971.8613944.631900
0.520.774051184.0081.6563664.155107
0.028.393021607.0701.4119063.602509
Table 6. Summary containing conditions and the most important equations numbering in Section 2, Section 3 and Section 4 for Wiener, Ornstein–Uhlenbeck and Feller diffusion processes.
Table 6. Summary containing conditions and the most important equations numbering in Section 2, Section 3 and Section 4 for Wiener, Ornstein–Uhlenbeck and Feller diffusion processes.
ConditionsResults—Equations Numbering
W i e n e r   p r o c e s s A 1 ( t ) = β ( t ) A 2 ( t ) = σ 2 ( t ) ( β ( t ) R , σ ( t ) > 0 ) β ( t ) = γ σ 2 ( t ) ( γ R , σ ( t ) > 0 ) a X ( x , t | x 0 , t 0 ) ( 17 ) g X ( 0 , t | x 0 , t 0 ) ( 18 ) P X ( 0 | x 0 , t 0 ) ( 19 )
β ( t ) = β , σ 2 ( t ) = σ 2 ( β R , σ > 0 ) g X ( S , t | x 0 ) ( 22 ) P X ( S | x 0 ) ( 23 ) a X ( x , t | x 0 ) ( 25 ) L X ( x | x 0 ) ( 26 )
O r n s t e i n U h l e n b e c k   p r o c e s s B 1 ( x , t ) = α ( t ) x + β ( t ) B 2 ( t ) = σ 2 ( t ) ( α ( t ) R , β ( t ) R , σ ( t ) > 0 ) β ( t ) = γ σ 2 ( t ) e A ( t | 0 ) ( γ R , α ( t ) R , σ ( t ) > 0 ) a Y ( x , t | x 0 , t 0 ) ( 33 ) g Y ( 0 , t | x 0 , t 0 ) ( 34 ) P Y ( 0 | x 0 , t 0 ) ( 35 )
α ( t ) = α , β ( t ) = β , σ 2 ( t ) = σ 2 ( α 0 , β R , σ > 0 ) g λ ( Y ) ( S | x 0 ) ( 41 ) P Y ( S | x 0 ) ( 42 ) a λ ( Y ) ( x | x 0 ) ( 46 ) L Y ( x | x 0 ) ( 47 ) , ( 48 )
F e l l e r   p r o c e s s C 1 ( x , t ) = α ( t ) x + β ( t ) C 2 ( x , t ) = 2 r ( t ) x ( α ( t ) R , β ( t ) R , r ( t ) > 0 , ( β ( t ) ξ r ( t ) , 0 ξ < 1 ) β ( t ) = ξ r ( t ) ( 0 ξ < 1 , r ( t ) > 0 ) a Z ( x , t | x 0 , t 0 ) ( 50 ) g Z ( 0 , t | x 0 , t 0 ) ( 52 ) P Z ( 0 | x 0 , t 0 ) ( 53 )
α ( t ) = α , β ( t ) = β , r ( t ) = r ( α R , β R , r > 0 , β < r ) a Z ( x , t | x 0 ) ( 56 ) , ( 57 ) g Z ( 0 , t | x 0 ) ( 58 ) P Z ( 0 | x 0 ) ( 59 ) L Z ( x | x 0 ) ( 61 ) , ( 62 )
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

Giorno, V.; Nobile, A.G. On the Absorbing Problems for Wiener, Ornstein–Uhlenbeck, and Feller Diffusion Processes: Similarities and Differences. Fractal Fract. 2023, 7, 11. https://doi.org/10.3390/fractalfract7010011

AMA Style

Giorno V, Nobile AG. On the Absorbing Problems for Wiener, Ornstein–Uhlenbeck, and Feller Diffusion Processes: Similarities and Differences. Fractal and Fractional. 2023; 7(1):11. https://doi.org/10.3390/fractalfract7010011

Chicago/Turabian Style

Giorno, Virginia, and Amelia G. Nobile. 2023. "On the Absorbing Problems for Wiener, Ornstein–Uhlenbeck, and Feller Diffusion Processes: Similarities and Differences" Fractal and Fractional 7, no. 1: 11. https://doi.org/10.3390/fractalfract7010011

APA Style

Giorno, V., & Nobile, A. G. (2023). On the Absorbing Problems for Wiener, Ornstein–Uhlenbeck, and Feller Diffusion Processes: Similarities and Differences. Fractal and Fractional, 7(1), 11. https://doi.org/10.3390/fractalfract7010011

Article Metrics

Back to TopTop