The entrance law of the excursion measure of the reflected process for some classes of L\'evy processes

We provide integral formulae for the Laplace transform of the entrance law of the reflected excursions for symmetric L\'evy processes in terms of their characteristic exponent. For subordinate Brownian motions and stable processes we express the density of the entrance law in terms of the generalized eigenfunctions for the semigroup of the process killed when exiting the positive half-line. We use the formulae to study in-depth properties of the density of the entrance law such as asymptotic behavior of its derivatives in time variable.


Introduction
It follows from excursion theory that the trajectories of a Lévy process can be decomposed using the excursions of the process reflected in its past infimum. This result justifies the importance of knowing the excursion measure of the reflected process and more particularly, the entrance law of this measure. There are also several interesting applications of this entrance law. First it is directly related to the potential measure of the time space ladder height process, see Lemma 1 in [3]. Moreover it provides a useful expression of the distribution density of the supremum of the Lévy process before fixed times, [3], [4]. More recently it has been involved in the study of the probability of creeping through curves of Lévy processes, [6].
In this article we obtain integral representations of the densities and the Laplace transforms of the entrance laws of the reflected excursions for two classes of real valued Lévy processes. The first class consists of symmetric Lévy processes, with a particular emphasis on subordinate Brownian motions, when the Lévy measure of the underlying subordinator has a completely monotone density. The other class is that of stable processes. The presented formulae for symmetric processes are given in terms of the corresponding Lévy-Kchintchin exponent Ψ(ξ) and the related generalized eigenfunctions introduced by M. Kwaśnicki in [14]. In the stable case, we based the calculations on the generalized eigenfunctions studied recently by A. Kuznetsov and M. Kwaśnicki in [13]. Then we used the formulae obtained for the entrance law densities to derive corresponding integral representations for supremum densities. Although the theory of Lévy processes is very rich and abounds in numerous general relationships, as those coming from the Wiener-Hopf factorizations, there are few examples where the explicit representations of the related densities are available. Apart from Brownian motion and Cauchy process, some series representations were recently found in [8], [12], [7] in the case of stable processes. A different approach was presented in [16], where the theory of Kwaśnicki's generalized eigenfunctions were used to described the first passage time density through a barrier for subordinate Brownian motions with regular Lévy measures. This concept was generalized to non-symmetric stable processes in [13]. In the present paper we stay in this framework and show that a similar approach leads to integral representations of the entrance law density, the supremum density or the density of joint distribution of the process itself and its supremum. Then we apply the obtained formulae to study the asymptotic behavior of the derivatives in time variable of the entrance law densities of the reflected excursions. Let us finally mention that these formulae can be used to perform numerical simulations and study in-depth properties of the process coupled with its past supremum.

Preliminaries
Let (X, P) be the real valued Lévy process whose characteristic exponent Ψ(ξ) is characterized in terms of the Lévy triplet (a, σ 2 , Π) by the Lévy-Kchintchin formula We write P x for the law of the process starting from x ∈ R. We denote by X * = −X the dual process and P * x stands for its law with respect to P x . The past supremum and past infimum of X before a deterministic time t ≥ 0 are For given t > 0 we write f t (dx) = P(X t ∈ dx) for the corresponding distribution and f t (x) stands for its density with respect to the Lebesgue measure on (0, ∞) whenever it exists. Recalling that the reflected processes X − X and X − X are Markovian, we write L t and L * t for their local times at 0 respectively, where these are normalized in the following way We write n (and n * ) for the Itô measure of the excursions away from 0 of the reflected process X −X (resp. X −X). Our main objects of studies are the corresponding entrance laws defined by q t (dx) = n(X t ∈ dx, t < ζ), q * t (dx) = n * (X t ∈ dx, t < ζ), t > 0, where ζ is the life time of the generic excursion and q t (x), q * t (x) denotes the densities on (0, ∞) of q t (dx) and q * t (dx), whenever they exist. In this paper, it will always be assumed that 0 is regular for both half-lines (−∞, 0) and (0, ∞). In this case, the double Laplace transform of q t (dx) is given by where κ(z, ξ) is the Laplace exponent of the ladder process (L −1 t , H t ), t < L(∞). Here L −1 t denotes the right continuous inverse of L t (ladder time process) and the ladder-height process is defined by H t = X L −1 t . Analogous relations hold for q * t (dx) and the Laplace exponent κ * (z, ξ) for the ladder process ((L * t ) −1 , H * t ). Formula (2.1) actually shows that q s (dx)ds is the potential measure of (L −1 , H). We denote by h the renewal function of the ladder height process H, that is In the light of Theorem 6 in [3], the entrance laws q t (dx) and q * t (dx) seem to be basic objects in the study of the supremum distributions. More precisely, under our assumption that 0 is regular for both negative and positive half-lines, the representation (4.4) from [3] reads as which, in particular, implies Finally, for x > 0, we denote by Q * x the law of the processes killed when exiting the positive half-line, i.e.
The law Q x is defined in the same way, but with respect to the dual process. The corresponding semigroups are defined as , for non-negative Borel functions f . We also write q * t (x, dy), q t (x, dy) and q * t (x, y), q t (x, y) for the corresponding transition probability measures and their densities whenever they exist. Recall that whenever q * t (x, ·) and q t (x, ·) are absolutely continuous, the duality relation holds q * t (x, y) = q t (y, x).

Symmetric Lévy processes and subordinated Brownian motions
This section is devoted to symmetric Lévy processes with some addition regularity assumptions on the Lévy-Kchintchin exponent Ψ(ξ) presented in details below. We also exclude compound Poisson processes from our considerations. Note that the symmetry assumptions simplify the general exposure presented in Preliminaries, where, roughly speaking, we can remove the notation with * . Moreover, the ladder time process is the 1/2-stable subordinator for every symmetric Lévy process, which implies that Finally, we recall the integral representation of the Laplace exponent of the ladder process where, in the symmetric case, Ψ(ξ) is a real-valued function. Our first result gives the expression for the Laplace transform of q t (dx) (for fixed t > 0) in the case of symmetric Lévy processes with increasing Lévy-Khintchin exponent. This is an analogue of Theorem 4.1 in [15], where the corresponding formula for X t was derived. Note that even though the formulae for the Laplace transforms of q t (dx) and P(X t ∈ dx) seem to be similar, passing from one to the other by using (2.3) and (3.1) is not straightforward. Theorem 1. Let (X, P) be a symmetric Lévy process that is not a compound Poisson process. Assume that the Lévy-Khintchin exponent Ψ(ξ) of (X, P) is increasing in ξ > 0. Then Proof. The proof is based on the same idea as in the proof of Theorem 4.1 in [15] with a slight modification of the arguments. For the completeness of the exposure and the convenience of the reader we present it below. We put ψ(ξ) = Ψ( √ ξ) for ξ > 0 and define Obviously, for fixed ξ > 0, the function ϕ(ξ, z) is a holomorphic function of z, which is positive for z > 0. Note also that lim z→0+ ϕ(ξ, z) = 0 (by monotone convergence) and lim z→∞ ϕ(ξ, z) = 1 (by dominated convergence). Moreover, as it was shown in [15], that for Im z > 0 we have Thus Arg( √ zϕ(ξ, z)) ∈ (0, π) for Im z > 0. This is equivalent to h ξ (z) = ϕ(ξ, z)/ √ z being a Stieltjes function (for fixed ξ). In general, a function g(z) is said to be a Stielties function if The constants and a measure appearing in the definitions of Stielties functions are given by Note that the last limit is understood in the sense of weak limit of measures. Since the constants appearing in the representation (3.4) for Stielties function h ξ (z) are zero. Moreover, for z = ψ(λ 2 ) we get Therefore, by (3.5), for every z > 0 we have Thus, the Laplace transform of the right-hand side of (3.3) is equal to 1/κ(z, ξ) and the theorem follows from uniqueness of the Laplace transform.
From now on, for the rest of the section, we will follow the approach presented in [16] and restrict our consideration to the case where (X, P) is a subordinate Brownian motion whose underlying subordinator has a complete monotone density. The process (X, P) has the latter form if and only if its characteristic exponent Ψ(ξ) can be written as Ψ(ξ) = ψ(ξ 2 ) for a complete Bernstein function ψ (see Proposition 2.3 in [14] where a 1 ≥ 0, a 2 ≥ 0 and µ(dζ) is a Radon measure on positive half-line such that min(ζ −1 , ζ −2 )µ(dζ) is finite. As in the Stielties function representation, the above-given constants and the measure µ are determined by suitable limits as follows The spectral theory of subordinate Brownian motion on a half-line was developed by M. Kwaśnicki in [14], where the generalized eigenfunctions F λ (x) of the transition semigroup Q t of the process (X, P) killed upon leaving the half-line [0, ∞) were constructed. Some additional properties of F λ (x) were also studied in [16]. For a fixed CBF ψ and λ > 0 the generalized eigenfunctions of Q t with eigenvalue e −tψ(λ 2 ) are given by where the phase shift ϑ λ belongs to [0, π/2) and is given by Recall the following upper-bounds (Proposition 4.3 and Proposition 4.5 in [16] respectively) The function G λ is the Laplace transform of the finite measure Here ψ + denotes the holomorphic extension of ψ in the complex upper half-plane. The Laplace transform of F λ (x) is given by Recall also the following estimates Proposition 5.4 in [16] states that for unbounded ψ such that lim sup λ→0 + ϑ λ < π/2 we have the following limiting behavior and the convergence is locally uniform in x ≥ 0.
The functions F λ (x) were used to find the integral representations for the density function of τ − 0 and its derivatives (see Theorem 1.5 in [16]). In the next Theorem we show that an analogous representation can be obtained for the density of the entrance law.
Proof. The proofs of both formulae are direct consequences of the integral representation (3.15), the relations (2.2), (2.3) and (3.1) together with the Fubini's theorem, which can be applied due to the integral condition (3.14).
The representation (3.15) enables to compute the derivatives of q t (x) and examine its behavior in two asymptotic regimes: as t goes to infinity and x goes to 0. It is described in the following theorem.

17)
where ψ −1 denotes the inverse of ψ, and the convergence is locally uniform in x. This also holds for α 0 = 1 with the additional assumption

19)
where p t denotes the density of the transition semigroup of (X, P).

using (3.20) and dominated convergence we obtain
Finally, the expression can by bounded for every t > t 1 by with some c 3 = c 3 (n, t 0 , t 1 ) > 0, which together with the regularity of ψ −1 at zero and estimates (3.20) implies that 1 vanishes uniformly in x, as t → ∞. Collecting all together we arrive at Because the justification of the fact that under assumption from point (b) we have follows in the same way as in the proof of Theorem 1.7 in [16], we omit the proof. Note that using (3.14) we can rewrite the last integral as where the last equality follows simply by integration by parts. Finally, the regular behavior of ψ at infinity implies that e −tψ(λ 2 ) is in L 1 (R, dλ), which in particular means that the transition probability density is given by the inverse Fourier transform Combining all together we get (3.19), which ends the proof.
In addition to numerical applications of our results, they can be used to obtain more transparent representations as in the following example related to the Cauchy process. Proposition 1. For the symmetric Cauchy process, i.e. ψ(ξ) = √ ξ, we have Then the density f t (x) of the past supremum at time t of (X, P) can be derived from the above expression together with (2.3) and (3.1).
Proof. Since ψ(ξ) = ξ 1/2 , ψ ′ (ξ) = 1/(2 √ ξ) the formula (3.15) reads as where we used the scaling property F λ (x) = F 1 (λx) = F x (λ). By the Plancherel's theorem we get, for fixed b ∈ (0, t), that where f (x) = √ 2x/π. The Laplace transform of f can easily be computed as follows Formula (3.11) gives Substituting u = z/s in the last integral we get Finally, the function studied in details in [11], is holomorphic in the region. We recall (see (3.13) in [11]) that and (see (4.1) in [11]) that where σ(z) = 1 for Im(z) > 0 and σ(z) = −1 for Im(z) < 0. Consequently, defining (for fixed x) the function of complex variable z it is easy to see that G x (z) is a meromorphic function on {z ∈ C : Re(z) < t} \ (−∞, 0] with single poles at ix and −ix. To evaluate the integral (3.21) we integrate G x over the (positively oriented) curve consisting of (see  Figure 1. The contour of integration. First we compute the residua of G x at ix and −ix. By (3.22), we have Since (t ± ix) 3/2 = (t 2 + x 2 ) 3/4 e ±3i/2 arctan(x/t) , we arrive at Using the relation (3.23) we obtain where the last integral, after substituting y = −xu, is equal to Using the bounds (3.12) we can write It implies that the integrals of G x (z) over γ 1 , γ 2 , γ 5 and γ 6 vanish as n goes to infinity.
Since G x (z) is bounded in the neighborhood of 0 (Re(z) > 0), the same holds for the integral over the semi-circle γ 7 . Now we can finish the computations by applying the residue theorem in order to get Taking into account (3.25) and (3.24) and dividing both sides by 2 √ π lead to the result.
Remark 2. It is also possible to find similar formula for the entrance law density of the symmetric α-stable process with index α ∈ (0, 1). Using the scaling property F λ (x) = F 1 (λx) and writing t (u) is the density of the α-stable subordinator we obtain The inner integral can be evaluated similarly as in Proposition 1.

Stable processes
For the rest of the paper we focus on stable processes and use the theory of the corresponding generalized eigenfunctions developed in [13]. We assume that X is a stable process with characteristic exponent Ψ(ξ) = |ξ| α e πiα(1/2−ρ)sign(ξ) , ξ ∈ R.
We exclude spectrally one-sided processes from our considerations, i.e. we assume that α ∈ (0, 1] and ρ ∈ (0, 1) or α ∈ (1, 2], but then we assume that ρ ∈ (1 − 1/α, 1/α). We write ρ * = 1 − ρ and define non-symmetric analogous of F 1 (x) defined in Section 3 for stable processes as follows The function S 2 (z) = S 2 (z; α) is the double sine function uniquely determined by the following functional equations together with the normalizing condition S 2 ((1 + α)/2) = 1 (see [9], [10] and Appendix A in [13] for equivalent definitions and further properties). We define F * (x) and G * (x) by the same formulae as in (4.1) and (4.2) but with ρ replaced by ρ * (and consequently ρ * replaced by (ρ * ) * = ρ). Note that whenever ρ > 1/2 the oscillations of F coming from the sine function are multiplied by the exponentially decreasing factor, but then F * oscillates exponentially, when x → ∞ and the situation is reversed for ρ < 1/2. The behaviour of F at zero is described by (see the proof of Lemma 2.8 in [13]) Although the constant √ α 2 S 2 (αρ) Γ(1+αρ * ) was not specified in [13], using (1.10) and (1.19) from [13], we obtain that Consequently, using the Karamata's Tauberian theorem and the Monotone Density Theorem we obtain (4.3). Moreover, if ρ > 1/2 then Even though the functions F and F * do not simultaneously belong to L 2 (0, ∞) (for ρ = 1/2), they can be understood as the generalized eigenfunctions of the semigroups Q * t and Q t , respectively (see Theorem 1.3 in [13]). Moreover, using Theorem 1.1 in [13] the transition probability density of the process X killed when exiting the positive half-line is given by whenever α > 1. Note that the restriction on α ensures that the exponential oscillations of F (or F * ) are suppressed by the factor e −tλ α , which makes the integral convergent. The formula (4.5) is the analogue of the integral representation for subordinate Brownian motions presented in [14]. Note also that assuming ρ = 1/2 we have F (x) = F * (x) = F 1 (x), where F λ (x) is the generalized eigenfunction defined in Section 3 for symmetric α-stable process.
By duality, we have the corresponding integral representation for q t (x) with F * (x) and ρ * replaced by F (x) and ρ.
The analogue of Theorem 3 can now be proved.