Arxiv:2004.02267V1 [Math.PR] 5 Apr 2020 E Words: Key Norm Literature

arXiv:2004.02267v1 [math.PR] 5 Apr 2020 e words: Key norm literature. for in results process corresponding studied the well provide assoc also Lévy density We the and process. function Green Lévy the dens and and stable density ordinator tempered potential of a temp behavior with normal asymptotic of motion the theory Brownian potential the the time-changing study by we article, this In aeyGaiainladEetoantcter.Teterm i The origin theory. its Electromagnetic has from and theory derived Gravitational This be namely 2009). to Helms, were e.g. nature (see of equation 1 forces of physics fundamental the the in origin that its has theory” “potential term The etc. 2020) al. et (Gupta 2 Mainard Polito, process and 2003; (Orsingher Laskin, process Poisson e.g. space-fractional (see proc process Student Poisson 1997), fractional al. et Lap alt (Mandelbrot fractional models to multifractal 1998), required is (Barndorff-Nielsen, it proce process time gamma Gaussian variance same include the processes at subordinated i and known where process models parent stochastic develop the to of way convenient a calle or is process subordinator tion stochastic called some Lévy process of non-decreasing S time other the biology. changing and by physics tained statistical gett finance, are in processes applications stochastic subordinated years recent In Introduction 1 etr8,SSNgr ujb100,India 140306, Punjab Nagar, SAS 81, Sector b a eateto ahmtclSine,Ida nttt fS of Institute Indian Sciences, Mathematical of Department eateto ahmtc,Ida nttt fTechnology of Institute Indian Mathematics, of Department oeta hoyo omlTmee tbeProcess Stable Tempered Normal of Theory Potential eprdsal uodntr ´v est,ptnilde potential Lévy density, subordinator, stable Tempered rnKumar Arun Abstract a n as Verma Harsh and 1 1)adtmee pc-rcinlPoisson space-fractional tempered and 012) n nraigatnindet hi wide their to due attention increasing ing ineEuainadRsac Mohali, Research and Education cience briae tcatcpoessaeob- are processes stochastic ubordinated h ne rcs.Nt htsubordina- that Note process. inner the oa,Rpaa,Pna 401 India 140001, Punjab Rupnagar, Ropar, se HyeadLoek,20) time- 2005), Leonenko, and (Heyde esses t etr,weeadmnn eifwas belief dominant a where century, 9th 2011), al. et Meerschaert 2004; al. et i t soitdwt eprdsal sub- stable tempered with associated ity h w elkontere fphysics of theories known well two the n aemto Kzbwk ta.2006), al. et (Kozubowski motion lace s(aa ta.19) omlinverse normal 1998), al. et (Madan ss ae ihtenra eprdstable tempered normal the with iated livreGusa rcs hc sa is which process Gaussian inverse al srqie orti oeproperties some retain to required is t sprn rotrpoesb some by process outer or parent as d oeta ucinwsfis associated first was function potential rdsal rcs hc sobtained is which process stable ered oetaswihstse Laplace’s satisfied which potentials rsm hrceitc.Sm well Some characteristics. some er uodntr rcsl,w study we Precisely, subordinator. b st,Genfunction. Green nsity, with the work done in moving a point charge from one point of space to other in the presence of 1 another external charge particle. The basic potential function varies as d , where d is the distance between the particles and the dimension of the space is greater than or equal to 3. This function has a property that it satisfies Laplace equation and such functions are called harmonic functions. From a mathematical point of view, potential theory is basically the study of harmonic functions (Armitage and Gardiner, 2001). Potential theory has intimate connection with probability theory (see e.g. Doob, 2001). One important connection is that the transition function of a Markov process can be used to define the Green function related to the potential theory. In this paper, we study the potential theory of normal tempered stable process and nornal inverse Gaussian process by realizing them as Brownian motion time-changed tempered stable subordiantor and inverse Gaussian process respectively. Precisely, we study the asymptotic behavior of potential density and Lévy density associated with different subordinators and also Green function and Lévy density associated with the subordinate processes. In easy language, potential measure represents the average time stay of a subordinator in a Borel subset of real numbers and Lévy measure quantifies the density of the num- ber of jumps per unit time of the subordinator. The Green function denoted by G(x,y)= G(x y) − for the Markov process is the expected amount of time spent at y by the Markov process started at x (see e.g. Doob, 2001; Liggett, 2010). The potential measure may be of interest to an investor, who is concern for the average time the stock prices stay in a particular price range.

In this paper, we will find the asymptotic behavior of potential density and Lévy density associated with the tempered stable subordinator and the inverse Gaussian subordinator. Further, we will find the Green function and Lévy density associated with the Brownian motion directed with these subordinators. More precisely, we find the asymptoric behaviour of Green function and Lévy density associated with normal tempered stable process and nornal inverse Gaussian process.

2 Preliminaries

In this section, we state some important theorems and definitions which will be used later to prove other results. We start with Tauberian theorems. For more details and proof of these theorems please refer to Theorems 1.7.1, 1.7.2, 1.7.1′ of Bingham et al. (1987). Tauberian theorems are used to find the asymptotic behaviour at infinity and zero of a real valued non-decreasing function by using the asymptotic behaviours of its Laplace-Stieltjes transform. If f and g are two function then f g means f converges asymptotically to 1. ∼ g Definition 2.1 (Slowly varying function). Let l be a positive measurable function, defined on some neighbourhood [x, ) of infinity, and satisfying l(λx) 1 as (x ) λ> 0; then l is said to be ∞ l(x) → → ∞ ∀ slowly varying at infinity (for details see Bingham et al. 1987, p.6).

Theorem 2.1. (Karamata’s Tauberian theorem) Let U be a non-decreasing right continuous function on R with U(x) = 0 x< 0. If l slowly varies and c 0, ρ 0, the following are equivalent: ∀ ≥ ≥

2 ρ a. As x (respectively x 0+), U(x) cx l(x) , →∞ → ∼ Γ(1+ρ) b. As s 0+ (respectively x ), U¯(s) cs ρl 1 . → →∞ ∼ − s Theorem 2.2. (Karamata’s monotone density theorem) Let U(x) = x u(y)dy. If U(x) cxρl(x) 0 ∼ as x , where c R, ρ R, l be a slowly varying function, and u is ultimately monotone, then → ∞ ∈ ∈ R u(x) cρxρ 1l(x). ∼ − Deﬁnition 2.2 (Subordinator). A subordinator is a 1-dimensional non-decreasing L´evy process. A subordinator S = (S(t) : t 0) is characterized by its Laplace transform ≥ sS(t) E e− = exp ( tφ(s)) . − h i The function φ is called the Laplace exponent of the subordinator (see Applebaum, 2009).

Deﬁnition 2.3. (Potential measure and density) The potential measure of a subordinator S(t) is deﬁned by E ∞ U(A)= 1(S(t) A)dt , ∈ Z0 where A is a Borel subset of (0, ). The potential measure has a density which is called the potential ∞ density of the subordinator (see e.g. Sikic et al. 2006).

Observe U(A) is the expected time the subordinator S(t) spends in the set A. Measure U has the following Laplace-Stieltjes transform

∞ st ∞ sS(t) ∞ tφ(s) 1 U¯(s)= e− dU(t)= E e− dt = e− dt = . (2.1) φ(s) Z0 Z0 Z0 If the Laplace exponent φ is a complete Bernstein function, then potential measure has a density u and u is completely monotone. If the subordinator has a potential density u(t) then, its Laplace transfrom is given by ∞ st 1 u¯(t)= e− u(t)dt = . (2.2) φ(s) Z0 If the L´evy measure µ of φ has a completely monotone density µ(t), i.e., ( 1)nDnµ 0 for every − ≥ non-negative integer n 1, which is equivalent to saying that φ is a complete Bernstein function. ≥ Let B(t)bea d-dimensional Brownian motion with transitional density given by

2 d x y d p(t,x,y) = (4πt)− 2 exp | − | , x,y R , t> 0. (2.3) − 4t ! ∈ Deﬁnition 2.4. (Green function for Markov process) The Green function for the Markov process is deﬁned by G(x,y)= ∞ p(t,x,y)dt, Z0 where p(t,x,y) is the transition function of the Markov process. It is the expected amount of time spent at y by the process started at x (Liggett, 2010).

3 Let B(t) be d-dimensional Brownain motion and S(t) be a subordinator then the process deﬁned by X(t) = B(S(t)) is called a subordinated process to Brownian motion . The Green function of X is given by (Sikic et al. 2006) G(x)= ∞ p(t, 0,x)u(t)dt, (2.4) Z0 where u(t) is the potential density of the subordinator S and p(t, 0,x) is the transition density of Brownian motion. Further, the L´evy density of X is given by (Sikic et al. 2006)

J(x)= ∞ p(t, 0,x)µ(t)dt, (2.5) Z0 where µ(t) is the L´evy density of S.

3 Tempered stable subordinator

Let S (t) be the tempered stable subordinator (TSS) with index β (0, 1) and tempering param- λ,β ∈ eter λ > 0. Note that, TSS are obtained by exponential tempering in the distribution of stable subordinator (see e.g. Rosinski, 2007). The marginal pdf for TSS Sλ,β(t) is given by

λx+λβt g (x,t)= e− f (x,t), s> 0, β (0, 1), (3.6) λ,β β ∈ where fβ(x,t) is the pdf of a stable subordinator. The advantage of tempered stable process over a stable process is that its all moments exist and its density is also inﬁnitely divisible. However, the process is not self-similar. The Laplace transform of the TSS is given by

β β sSλ,β(t) t((s+λ) λ ) E(e− )= e− − .

The L´evy density associated with the tempered stable subordinator is given by (see e.g. Rosi´nski, 2007; Kumar and Vellaisamy, 2015)

ce βx µ(x)= − , x> 0, (3.7) xλ+1 where c = β/Γ(1 β). Using direct computation, we can show that ( 1)nDnµ 0, n 1 for TSS. − − ≥ ≥ Hence the corresponding potential measure has a potential density. We next discuss the asymptotic behavior of the potential density of the TSS.

Theorem 3.1. For β (0, 1) and λ> 0, the potential density u(x) for TSS have following asymptotic ∈ forms

λ 1 1. u(x) λx − , as x 0+, ∼ Γ(1+λ) → 1 2. u(x) λ 1 , as x . ∼ λβ − →∞ Proof.

4 1. The Laplace exponent of the Tempered stable subordinator is given by

φ(s) = (s + β)λ βλ, − where λ (0, 1) and β > 0. Using (2.1), the Laplace transform of the potential measure U of ∈ tempered stable subordinator will be

¯ 1 1 1 + U(s)= = λ λ 1 , s 0 . φ(s) (s + β) βλ ∼ λβ − s → − Therefore, by Theorem 2.1, the potential measure is 1 U(x) λ 1 x, x . ∼ λβ − Γ(2) →∞

Hence, by Theorem 2.2, the potential density 1 1 u(x) λ 1 = λ 1 , x . ∼ λβ − Γ(2) λβ − →∞

2. As s , we have φ(s) sλ, therefore, U¯(s) 1 . Thus, by Theorem 2.1, we have →∞ ∼ ∼ sλ xλ U(x) , x 0+. ∼ Γ(1 + λ) →

Using Theorem 2.2, λxλ 1 u(x) − , x 0+. ∼ Γ(1 + λ) →

In next result, we ﬁnd the exact form of the potential density.

Theorem 3.2. For λ (0, 1) and β > 0, we have ∈ λ sin(πλ) ∞ exp( xu)u u(x) = exp( βx) − du. − π u2λ + β2λ 2uλβλ cos(πλ) Z0 − Proof. Let φ(s) = f(s + β) = (s + β)λ βλ be the Laplace exponent of the tempered stable − subordinator, where f(s)= sλ βλ. Then using the properties of inverse Laplace transform we have −

1 1 1 1 − u¯(s)= (t) = exp( βx) − (x), L φ(s) − L f(s) where 1[¯u(s)] is the inverse Laplace transform of the potential density related to the tempered L− stable subordinator. Now to find the inverse Laplace transform of the potential density we will first find 1 1 (x) using complex inversion formula (Schiff, 1999). Since s = 0 is a branch point of L− f(s) h i

5 Figure 1: Contour C anti-clockwise

1 f(s) , so we take a branch cut along non-positive real line to make the function single valued as shown in the Figure (1). Inside and on the contour the function is analytic so by Cauchy’s theorem

1 exp(sx) ds = 0. 2πi f(s) ZC Now 1 exp(sx) 1 exp(sx) 1 exp(sx) 1 exp(sx) ds = ds + ds + ds 2πi f(s) 2πi f(s) 2πi f(s) 2πi f(s) ZC ZAB ZBC ZCD 1 exp(sx) 1 exp(sx) 1 exp(sx) + ds + ds + ds = 0. 2πi f(s) 2πi f(s) 2πi f(s) ZDE ZEF ZFA also exp(sx) exp(sx) exp(sx) ds = ds = ds = 0, f(s) f(s) f(s) ZAB ZCD ZEF (see Schiﬀ, 1999 for details). We know from Schiﬀ (1999) that as r 0 and R → →∞

1 exp(sx) 1 1 ds = L− (t), (3.8) 2πi f(s) f(s) ZFA

6 therefore, 1 1 1 exp(sx) 1 exp(sx) − (x)= ds + ds . (3.9) L f(s) − 2πi f(s) 2πi f(s) ZBC ZDE Consider, r exp(sx) exp(sx) − exp(sx) ds = λ λ ds = λ λ ds. BC f(s) BC s β R − s β Z Z − Z− − Let s = ueiπ then ds = du, thus, − exp(sx) r exp( ux) ds = − du. sλ βλ −uλeiπλ βλ ZBC − ZR − Taking the limits r 0 and R on both side of the above equation we get → →∞ exp(sx) ∞ exp( ux) ds = − du. (3.10) sλ βλ uλ exp(iπλ) βλ ZBC − Z0 − Now R exp(sx) exp(sx) − exp(sx) ds = λ λ ds = λ λ ds, DE f(s) DE s β r − s β Z Z − Z− − For DE, take s = ue iπ, which gives ds = du, − − exp(sx) R exp( xu) ds = − du. sλ βλ −uλe iπλ βλ ZDE − Zr − − Again we take the limits r 0 and R on both side of the above equation we get → →∞ exp(sx) ∞ exp( xu) ds = − du. (3.11) sλ βλ −uλ exp( iπλ) βλ ZDE − Z0 − − From equations (3.9), (3.10) and (3.11), we get

1 1 1 ∞ exp( xu) ∞ exp( xu) − (x)= − du − du L f(s) −2πi uλ exp(iπλ) βλ − uλ exp( iπλ) βλ Z0 − Z0 − − 1 1 1 = ∞ exp( xu) du −2πi − uλ exp(iπλ) βλ − uλ exp( iπλ) βλ Z0 − − − λ 1 ∞ 2i sin(πλ)u = exp( xu) − du −2πi − u2λ + β2λ 2 cos(πλ)uλ Z0 −λ sin(πλ) ∞ exp( xu)u = − du π u2λ + β2λ 2uλβλ cos(πλ) Z0 − Thus, the potential density is

λ 1 sin(πλ) ∞ exp( xu)u u(x)= − [¯u(s)](x) = exp( βx) − du. L − π u2λ + β2λ 2uλβλ cos(πλ) Z0 −

7 4 Normal tempered stable process

In this section, we compute the asymptotic behavior of the Green function and L´evy density associated with the Brownian motion time-changed by tempered stable subordinator which is called the normal tempered stable process see Barndorﬀ-Nielsen and Shephard (2001). Let B(t) be the

Brownian motion with transition probability given in (2.3) and Sλ,β(t) be the TSS, then the time change process X(t)= B(Sλ,β(t)) is called normal tempered stable process.

Proposition 4.1. For λ (0, 1) and β > 0, we have ∈ d 2λ Γ( −2 ) 2λ d + 1. G(x) d x − , x 0 , ∼ π 2 4λΓ(λ) | | | |→

d 2 Γ( −2 ) 2 d 2. G(x) d x − , x . 2 λ 1 ∼ 4π λβ − | | | |→∞ Proof.

1. As s , φ(s) sλ, then it follows directly from Theorem 3.1 of Rao et al. (2006) that →∞ ∼ d 2λ Γ −2 2λ d + G(x) d x − , x 0 . ∼ π 24λΓ(λ)| | | |→

2. As s 0+, φ(s) λβλ 1s, then it follows directly from Theorem 3.3 of Rao et al. (2006) that → ∼ − d 2 1 Γ −2 2 d G(x) d 2 x − , x ∼ 2 2 λ 1 Γ | | | |→∞ π 2 λβ − 2 or d 2 Γ −2 2 d G(x) d x − , x . ∼ 2 λ 1 | | | |→∞ 4π λβ −

We now find the asymptotic behaviour of the Lévy density J(x) of the subordinated process. We can not use the same method as in Sikic et al. (2006) because the asymptotic behavior of Lévy density µ(x) of the Lévy measure µ associated with the tempered stable subordinator does not follow assumption 2 of Lemma 3.1 so we can’t use the results given in Sikic et al. (2006).

Theorem 4.1. For λ (0, 1) and β > 0, we have ∈ 4λ+1 d (2λ+d) 1. J(x) d λ + 2 c x − , x 0, ∼ π 2 | | | |→

2λ d+1 2λ+d 1 −2 4 − 2λ+d 1 2 β 2 − 2. J(x) d 1 c x − exp √β x , x . ∼ π −2 | | − | | | |→∞

8 Proof. Using (2.5) and (3.7), the L´evy density associated with the subordinated process 2 ∞ ∞ d ( x ) c exp( βt) J(x)= p(t, 0,x)µ(t)dt = (4πt)− 2 exp | | − dt − 4t tλ+1 Z0 Z0 2 c ∞ d ( x ) (λ+1+ 2 ) = d t− exp | | βt dt d 2 − 4t − 2 π Z0 2 c ∞ d 1 ( x ) (λ+1+ 2 ) = d t− exp | | + 2βt dt. d 2 −2 2t 2 π Z0 x y x dy | | | | Let t = 2√β , then dt = 2√β , thus, (λ+1+ d ) c ∞ x y − 2 1 √β x y x dy J(x)= d | | exp β x y + | | | | d 2 2√β −2 | | y 2√β 2 π Z0 (λ+ d ) p c x − 2 ∞ d 1 1 (λ+1+ 2 ) = d | | y− exp β x y + dy d 2 2√β −2 | | y 2 π Z0 (λ d ) 1 (λ+ d ) p 2 − 2 β 2 2 d ∞ d 1 1 (λ+ 2 ) (λ+1+ 2 ) = d c x − y− exp β x y + dy 2 | | −2 | | y π Z0 (λ d ) 1 (λ+ d ) p 2 − 2 β 2 2 d (λ+ 2 ) = d c x − 2K ν (ω), (4.12) π 2 | | − where 1 ∞ (λ+1+ d ) 1 1 K ν(ω)= y− 2 exp β x y + dy, − 2 −2 | | y Z0 pd is a modiﬁed Bessel function of third type with ν = λ + and ω = √β x . Further, K ν(ω) = 2 | | − Kν (ω) (Jorgensen, 1982). 1. As x 0, ω 0 thus, using A.7 of Jorgensen (1982), we get | |→ → d λ+ d 1 (λ+ d ) K (ω) λ λ + 2 2 − ( β x )− 2 . ν ∼ 2 | | We get p (λ d ) 1 (λ+ d ) 2 − 2 β 2 2 d d d d (λ+ 2 ) λ+ 2 +1 (λ+ 2 ) J(x) d c x − 2λ λ + 2 ( β x )− ∼ 2 | | 2 | | π λ+1 p 4 d (2λ+d) d λ + c x − , x 0. ∼ 2 2 | | | |→ π 2. As x , ω , thus from Jorgensen (1982), we get | |→∞ →∞ 1 π 2 K (ω) exp β x β x − . ν ∼ 2 − | | | | r Thus, p p (λ d ) 1 (λ+ d ) 1 2 − 2 β 2 2 d π 2 (λ+ 2 ) − J(x) d c x − 2 exp β x β x ∼ 2 | | 2 − | | | | π r 2λ d+1 2λ+d 1 p p 2 −2 β 4 − 2λ+d 1 2 − d 1 c x − exp β x , x . ∼ π −2 | | − | | | |→∞ p 9 5 Inverse Gaussian Subordinator

The PDF of a random variable following an inverse Gaussian (IG) distribution with parameters δ and λ is given by Applebaum (2009)

δλ 2 δe 1 ( δ +λ2x) g(x)= e− 2 x , x> 0,δ,λ > 0. √2πx3 IG distributions are inﬁnitely divisible and the corresponding L´evy process is called IG subordinator. Let S(t) be the IG subordinator which has alternative representation

S(t) = inf s> 0 : B(s)+ λs = δt , (5.13) { } where B(t) is the standard one-dimensional Brownian motion. The Laplace transform of the inverse Gaussian subordinator is given by

sS(t) 2 E[e− ] = exp tδ 2s + λ λ . − − h p i The L´evy density for IG subrodinator is

2 δ λ x µ(x)= e− 2 , x> 0. (5.14) √2πx3 Next, we discuss about the asymptotic behaviour of potential density of the IG subordinator.

Theorem 5.1. For δ,λ > 0, we have

1. u(x) 1 , x 0+, ∼ δ√2πx → 2. u(x) λ , x . ∼ δ →∞ Proof.

1. The Laplace exponent of the Inverse Gaussian subordinator is given by

φ(s)= δ 2s + λ2 λ , − p where λ (0, 1) and β > 0. Therefore, the Laplace transform of the potential measure U of ∈ tempered stable subordinator will be 1 U¯(s)= . δ √2s + λ2 λ − 1 1 √ 2 2s 2 + 2s 2 s Since δ 2s + λ λ = δλ 1+ λ2 1 and as s 0 1+ λ2 1+ λ2 , therefore as − 1 − → ∼ + s 2 h δ i s 0 , φ(s) δλ 1+ 2 1 = s, thus, → ∼ λ − λ h i λ U¯(s) , s 0+. ∼ δs →

10 Therefore, by Theorem 2.1, we have λx λx U(x) = , x , ∼ δΓ(2) δ →∞ and hence, by Theorem 2.2, the potential density is given by λ u(x) , x . ∼ δ →∞

2. As s , φ(s) δ√2s, therefore, U¯(s) 1 . Thus, by Theorem 2.1, we have →∞ ∼ ∼ δ√2s 1 √2x U(x) √x = , x 0+. ∼ √ 3 δ√π → δ 2Γ 2 Therefore, by Theorem 2.2, we have 1 u(x) , x 0+. ∼ δ√2πx →

Next we ﬁnd the exact form of the potential density associated with the inverse Gaussian subordinator.

Theorem 5.2. For δ,λ > 0, the potential density of IG subordinator is given by

λ2x 1 exp 2 λ λ√x u(x)= − + erfc . (5.15) √  √πx √ − √  2δ 2 2   2 2 Proof. Let φ(s)= f(s+ λ )= √2δ s + λ λ be the Laplace exponent of the inverse Gaussian 2 2 − √2 q subordinator, where f(s)= √2δ √s λ . Then using the properties of inverse Laplace transform, − √2 we get 2 1 1 λ x 1 1 u(s)= − u¯(s)= (x) = exp − (x). (5.16) L φ(s) − 2 L f(s) Using formula 128 on page 16 of Poularikas (1998)

2 1 1 1 1 1 1 λ λ x λ√x − (x)= − (x)= + exp erfc L f(s) L √2δ √s λ  √2δ √πx √2 2 − √2 − √2   (5.17)

Using (5.16) and (5.17), we get

2 exp λ x 1 1 2 λ λ√x u(x)= − [¯u(s)] = − + erfc . L √  √πx √ − √  2δ 2 2  

11 Corollary 5.1. The asymptotic behaviour in Theorem 5.1 can be proved using (5.15).

Proof. Note that

λ2x 1 exp 2 λ λ√x u(x)= − + erfc √2δ  √πx √2 − √2  1  λ2x λ√πx λ√x 1 = exp + erfc , x 0+, √ − 2 √ − √ ∼ √ → δ 2πx 2 2 δ 2πx since erfc(0) = 1. Further for large x, we have (see e.g. Abramowitz and Stegun, 1992)

x2 e− ∞ n 1 3 5 (2n 1) erfc(x) = 1+ ( 1) · · · · ·2 n − . (5.18) x√π " − (2x ) # nX=1 Using, erfc( x) = 2 erfc(x) and (5.18), we have − − λ u(x) , x . ∼ δ →∞

Next we compute the asymptotic behavior of Green function and L´evy density associated with the normal inverse Gaussian process which is obtained by subordinating Brownian motion with inverse Gaussian subordinator.

6 Normal Inverse Gaussian Process

Let B(t) be the d-dimensional Brownian motion with transition probability given in (2.3) and S(t) be the IG subordinator deﬁned in (5.13). The process deﬁned by

Y (t)= B(S(t)), is called a d-dimensional normal inverse Gaussian process without drift (see e.g. Barndorﬀ-Nielsen, 1998). The asymptotic form of associated L´evy density of this d-dimensional Normal inverse Gaussian process Y (t) can be calculated similar to Theorem (4.1) and is given by

Proposition 6.1. For δ,λ > 0, we have

d 1 1 Γ( −2 ) 1 d + 1. G(x) 3 d 1 x − , x 0 , ∼ 2 2 π 2 δ Γ( 2 ) | | | |→

12 λ d 2 2 d 2. G(x) d Γ −2 x − , x . ∼ 4π 2 δ | | | |→∞ Proof.

1. As s , φ(s) δ√2s, then it follows directly from Theorem 3.1 of Rao et al. (2006) that →∞ ∼ d 1 1 Γ −2 1 d + G(x) 3 d 1 x − , x 0 . ∼ 2 2 | | | |→ 2 π δ Γ 2 2. As s 0+, φ(s) δ s, then using Theorem 3.3 of Rao et al. (2006), it follows → ∼ λ

λ d 2 2 d G(x) d Γ − x − , x . ∼ 2 2 | | | |→∞ 4π δ

7 Conclusions

In this paper, we mainly computed the asymptotic behavior of potential density and L´evy density associated with the tempered stable subordinator and inverse Gaussian subordinator and also the asymptotic behavior of Green function and L´evy density associated with the Brownian motion timed changed by tempered stable subordinator and inverse Gaussian process.

Acknowledgments: A. Kumar is supported by Science and Engineering Research Board (SERB), India project MTR/2019/000286. This work is part of H. Verma MS thesis work.

References

[1] Abramowitz, M., Stegun, I.A. (Eds.), 1992. Handbook of Mathematical Functions with Formu- las, Graphs and Mathematical Tables. Dover, New York.

[2] Applebaum, D., 2009. L´evy Processes and Stochastic Calculus (second ed.). Cambridge Univer- sity Press, Cambridge.

[3] Armitage, D.H., Gardiner, S.J., 2001. Classical Potential Theory. Springer-Verlag, London.

[4] Barndorﬀ-Nielsen, O.E., 1998. Processes of normal inverse Gaussian type. Finance Stochast. 2, 41-68.

[5] Barndorﬀ-Nielsen O.E., Shephard N., 2001. Normal modiﬁed stable processes. Theory Probab. Math. Statist. 65, 1-19.

[6] Bingham, N.H., Goldie, C.M., Teugels, J.L., 1987. Regular Variation. Cambridge University Press, Cambridge.

13 [7] Doob, J.L., 2001. Classical Potential Theory and its Probabilistic Counterpart. Springer-Verlag, Berlin Heidelberg.

[8] Gupta, N., Kumar, A. and Leonenko, N.N., 2020. Tempered Fractional Poisson Processes and Fractional Equations with Z-Transform. Stoch. Anal. Appl. To Appear.

[9] Helms, L.L., 2009. Potential Theory. Springer-Verlag, London.

[10] Heyde, C.C., Leonenko, N.N., 2005. Student processes. Adv. Appl. Probab. 37, 342-365.

[11] Jorgensen, B., 1982. Statistical Properties of the Generalized Inverse Gaussian Distribution. Springer-Verlag, New York.

[12] Kumar, A., Vellaisamy, P., 2015. Inverse tempered stable subordinators. Statist. Probab. Lett. 103, 134-141.

[13] Kozubowski, T.J., Meerschaert, M.M., Podgorski, K., 2006. Fractional Laplace motion. Adv. Appl. Prob. 38, 451-464.

[14] Laskin, N., 2003. Fractional Poisson process. Commun. Nonlinear Sci. Numer. Simul. 8, 201-213.

[15] Liggett, T.M., 2010. Continuous Time Markov Processes: An Introduction. American Mathe- matical Society.

[16] Madan, D.B., Carr, P., Chang, E.C., 1998. The variance gamma process and option pricing. Eur. Fin. Rev. 2, 74-105.

[17] Mandelbrot, B.B., Fisher, A., Calvet, L., 1997. A multifractal model of asset returns. Cowles Foundation discussion paper no. 1164.

[18] Mainardi, F., Gorenﬂo, R., Scalas, E., 2004. A fractional generalization of the Poisson processes. Vietnam J. Math. 32, 53-64.

[19] Meerschaert, M.M., Nane, E., Vellaisamy, P., 2011. The fractional Poisson process and the inverse stable subordinator. Electron. J. Probab., 16, 1600-1620.

[20] Orsingher, E., Polito, F., 2012. The space-fractional Poisson process. Statist. Probab. Lett. 82, 852-858.

[21] Poularikas, A.D., 1998. Handbook of Formulas and Tables for Signal Processing. Springer Science & Business Media.

[22] Rao, M., Song, R., Vondracek, Z., 2006. Green function estimates and Harnack inequality for subordinate Brownian motions. Potential Anal. 25, 1-27.

[23] Rosi´nski, J., 2007. Tempering stable processes. Stochastic Process Appl. 117, 677-707.

14 [24] Schiﬀ, J.L., 1999. The Laplace Transform: Theory and Applications. Springer-Verlag New York.

[25] Sikic, H., Song, R., Vondracek, Z., 2006. Potential Theory of Geometric Stable Processes. Probab. Theory Relat. Fields. 135, 547-575.