arxiv: v2 [math.pr] 4 Jul 2008

Transcription

1 Applied Probability Trust (7 January 4) OLD AND NEW EXAMPLES OF SCALE FUNCTIONS FOR SPECTRALLY NEGATIVE LÉVY PROCESSES F. HUBALEK, Vienna University of Technology A. E. KYPRIANOU, The University of Bath arxiv:8.393v [math.pr] 4 Jul 8 Abstract We give a review of the state of the art with regard to the theory of scale functions for spectrally negative Lévy processes. From this we introduce a general method for generating new families of scale functions. Using this method we introduce a new family of scale functions belonging to the Gaussian Tempered Stable Convolution (GTSC) class. We give particular emphasis to special cases as well as cross-referencing their analytical behaviour against known general considerations. Keywords: Scale functions, Spectrally negative Lévy processes, Mittag-Leffler functions, Wiener-Hopf factorization Mathematics Subject Classification: Primary 6G5, 6J75 Secondary 6G99. Spectrally negative Lévy processes and scale functions Let X = {X t : t } be a Lévy process defined on a filtered probability space (Ω, F, F, P), where {F t : t } is the filtration generated by X satisfying the usual conditions. For R denote by P the law of X when it is started at and write simply P = P. Accordingly we shall write E and E for the associated epectation operators. In this paper we shall assume throughout that X is spectrally negative meaning here that it has no positive jumps and that it is not the negative of a subordinator. It is well known that the latter allows us to talk about the Laplace eponent ψ(θ) := log E[e θx ] for R(θ) where in particular we have the Lévy-Khintchine representation ψ(θ) = aθ + σ θ + (e θ θ {> } )Π(d) () (,) where a R, σ is the Gaussian coefficient and Π is a measure concentrated on (, ) satisfying (,) ( )Π(d) <. The, so-called, Lévy triple (a, σ, Π) completely characterises the process X. For later reference we also introduce the function Φ : [, ) [, ) as the right inverse of ψ on (, ) so that for all q Φ(q) = sup{θ : ψ(θ) = q}. () Postal address: Vienna University of Technology, Financial and Actuarial Mathematics, Wiedner Hauptstraße 8/5, A-4 Vienna, Austria. fhubalek@fam.tuwien.ac.at Postal address: Department of Mathematical Sciences, The University of Bath, Claverton Down, Bath BA 7AY. a.kyprianou@bath.ac.uk

2 F. HUBALEK AND A. E. KYPRIANOU Note that it is straightforward to show that ψ is a strictly conve function which is zero at the origin and tends to infinity at infinity and hence there are at most two solutions of the equation ψ(θ) = q. Suppose now we define the stopping times for each R τ + = inf{t > : X t > } and τ = inf{t > : X t < }. (3) A fluctuation identity with a long history concerns the probability that X eits an interval [, a] (where a > ) into (a, ) before eiting into (, ) when issued at [, a]. In particular it is known that E (e qτ+ a {τ + a <τ }) = W (q) () (4) W (q) (a) where (, a], q and the function W (q) : R [, ) is defined up to a multiplicative constant as follows. On (, ) we have W (q) () = and otherwise W (q) is a continuous function (right continuous at the origin) with Laplace transform e θ W (q) ()d = ψ(θ) q for R(θ) > Φ(q). (5) The functions {W (q) : q } are known as scale functions and for convenience and consistency with eisting literature we write W in place of W (). Identity (4) eemplifies the relation between scale functions for q = and the classical ruin problem. Indeed setting q = we have that P (τ < τ+ a ) = W()/W(a) and so assuming that ψ (+) > and taking limits as a it is possible to deduce that W( ) = ψ (+) and hence P (τ < ) = ψ (+)W(). (6) It is in this contet of ruin theory that scale functions make their earliest appearance in the works of Zolotarev (964), Takács (966) and then later either eplicitly or implicitly in the work of Emery (973) and Korolyuk (974, 975), Suprun (976) and Rogers (99). The real value of scale functions as a class with which one may epress a whole range of fluctuation identities for spectrally negative Lévy processes became apparent in the work of Chaumont (994, 996) of Bertoin (996, 997) and an ensemble of subsequent articles; see for eample Lambert (), Avram et al. (4), Pistorius (3, 4, 5, 7), Kyprianou and Palmowski (6), Doney and Kyprianou (6) and Doney (99, 5a, 7). Moreover with the advent of these new fluctuation identities and a better understanding of the analytical properties of the function W (q) came the possibility of revisiting and solving a number of classical and modern problems from applied probability, but now with the underlying source of randomness being a general spectrally negative Lévy processes. For eample, in the theory of optimal stopping Avram et al. (, 4); Alili and Kyprianou (5) and Kyprianou (6), in the theory of optimal control Avram et al. (7) and Loeffen (7), in the theory of queuing and storage models Dube et al. (4) and Bekker et al. (7), in the theory of branching processes Bingham (976) and Lambert (7), in the theory of insurance risk and ruin Chiu and Yin (5), Klüppelberg et al. (4), Klüppelberg and Kyprianou (6), Doney and Kyprianou (6), Kyprianou and Palmowski (7), Renaud and Zhou (7), Loeffen (7) and Kyprianou et al. (8), in the theory of credit risk Hilberink and Rogers () and Kyprianou and Surya (7) and in the theory of fragmentation Krell (7).

3 3 Although scale functions are now firmly embedded within the theory of spectrally negative Lévy processes and their applications, and although there is a reasonable understanding of how they behave analytically (see for eample the summary in Kyprianou (6)), one of their main failings from a practical point of view is that there are limited number of eamples which are available off-the-shelf for modeling purposes. The original purpose of this paper was to introduce a previously unknown family of scale functions. However, in doing so, we came across a number of additional eamples which were seemingly less well known in the literature or are themselves new and concurrent with this article. We shall therefore initially spend some time documenting known eamples of scale functions. It seemed then sensible to also document known general analytical properties of scale functions and this we have also done in Section 3. Following that we shall epose a new family of eplicit eamples of scale functions as well specifying in detail the associated Lévy processes. By eplicit we mean that, as in the below, known eamples, we are able to give eact analytical epressions for the scale functions. By specifying in detail we mean that we shall describe the Lévy triple (a, σ, Π) of the associated spectrally negative Lévy process and in particular it will turn out that a density for the Lévy measure can be obtained and we are able to identify it as a known function. Moreover, it turns out that the associated Lévy processes discussed may also be described as the independent sums of other, well known Lévy processes.. Known eamples of scale functions At the time of writing, and to the best of the authors knowledge, there are essentially si main eamples of eplicit classes; three of which deal with the case of a compound Poisson process with negative jumps and positive drift.. Brownian motion with drift. The associated Laplace eponent is written ψ(θ) = σ θ /+µθ, where σ > and µ R then, for, W (q) () = qσ + µ e µ/σ sinh( σ qσ + µ) (7) for q, where in the case that q = µ = the above epression is to be taken in the limiting sense.. Spectrally negative stable processes with stability parameter β (, ). When ψ(θ) = θ β with β (, ), we have for W (q) () = β β E β,(q β ) (8) for q. where E β, (z) = k zk /Γ( + βk) is the Mittag-Leffler function. There also eists an epression for the scale function of the aforementioned β-stable process but now with a strictly positive drift (but only for the case q = ) implicitly in the paper of Furrer (998). Indeed, by considering the epression there for the ruin probability one may etract the following, W() = c ( Eβ ( c β ) ) (9) for. Note that by taking c one obtains an epression which agrees with (8), namely β /Γ(β).

4 4 F. HUBALEK AND A. E. KYPRIANOU 3. Generally speaking it is possible to deal with a strictly positive linear drift minus a compound Poisson process process with positive jumps with a rational Laplace transform. The calculations are however rather etensive. See for eample Mordecki and Lewis (5) where the scale function can be etracted from their epression for the ruin probability. In the special case that the drift is denoted c and jumps arrive at rate λ and are eponentially distributed with parameter µ such that c λ/µ >, that is to say ψ(θ) = cθ λ( µ(µ+ θ) ), we have for the case q = and, W() = ( + λ ) c cµ λ ( λ) e(µ c ). () The epression for the scale function in the case that q > with eponentially distributed jumps, and indeed for the slightly more general case where an independent Gaussian component is added in, can in principle be etracted from Kou and Wang (3). 4. This eample is taken from Asmussen () and again comes from the case of a known ruin probability giving a scale function for the case q = and positive drift. Consider a spectrally negative compound Poisson process whose jumps are eactly of size α (, ), with arrival rate λ > and with positive drift c > such that c λα >. In that case, ψ(θ) = cθ λ( e αθ ) and for we have W() = c /α n= where /α is the integer part of /α. e λ(αn )/c n! ( ) n λ (αn ) n () c 5. The net eample comes from Abate and Whitt (999) who themselves refer to older work of Boma and Cohen (998). Consider a unit-rate linear drift minus a compound Poisson process of rate λ > with jumps whose distribution F has Laplace transform given by e θ θ F(d) = (µ + θ)( + θ) for R(θ) The latter Laplace transform corresponds to a random variable whose mean is equal to µ > and whose tail distribution takes the form F(, ) = ( + )η() π, (3) when µ = and otherwise when µ ( ) F(, ) = (η() µη(µ )), (4) µ where η() = e erfc( ). It is assumed that the underlying Lévy process drifts to. Since E(X ) = λ/µ the latter assumption is tantamount to λ/µ <. Said another way, we are interested in a spectrally negative Lévy process with Laplace eponent given by λθ ψ(θ) = θ (µ + θ)( + θ). (5) ()

5 5 For the scale function with q = it is known that for W() = ( λ/µ ) (ν η(ν λ/µ ν ν ) ν η(ν )). (6) and ν, = + µ ± ( + µ ) ( λ ) µ. (7) µ 6. The final two eamples, both scale functions only for the case q =, appeared very recently in the theory of positive self-similar Markov processes, see Chaumont et al. (7). See also Chaumont and Caballero (6) for the origin of the underling Lévy processes. The first eample is the scale function which belongs to a spectrally negative Lévy process with no Gaussian component, whose Lévy measure takes the form Π(dy) = where β (, ) and whose Laplace eponent takes the form e(β )y (e y dy, y < (8) ) β+ ψ(θ) = Γ(θ + β) Γ(θ )Γ(β) (9) for R(θ). Note that ψ (+) < and hence the process drifts to. In that case it was found that for W() = ( e ) β e. () The second eample is the scale function associated with the aforementioned Lévy process when conditioned to drift to. It follows that there is still no Gaussian component and the Lévy measure takes the form Π(dy) = and the associated Laplace eponent is given by e βy (e y dy, y < () ) β+ ψ(θ) = Γ(θ + β) Γ(θ)Γ(β) for R(θ). The scale function is then given for by W() = ( e ) β. () 3. Known analytic properties of scale functions. Although at the time of writing further eamples other than those above are lacking, there are a collection of general properties known for scale functions, mostly concerning their behaviour at and. For later reference in this tet and to give credibility to some we review them briefly here. As usual, (a, σ, Π) denotes the Lévy triple of a general spectrally negative Lévy process.

6 6 F. HUBALEK AND A. E. KYPRIANOU Smoothness The following facts are taken from Lambert (), Chan and Kyprianou (7), Kyprianou et al. (8) and Doney (5b). It is known that if X has paths of bounded variation then, for all q, W (q) (, ) C (, ) if and only if Π has no atoms. In the case that X has paths of unbounded variation, it is known that, for all q, W (q) (, ) C (, ). Moreover if σ > then C (, ) may be replaced by C (, ). It was also noted by Renming Song (see the remarks in Chan and Kyprianou (7)) that if Π(, ) is completely monotone then W (q) (, ) C (, ). Concavity and conveity It was shown in Loeffen (7) that the latter assumption that Π(, ) is completely monotone also implies that W (q) () is conve for q >. Note in particular, the latter implies that there eists an a such that W (q) is concave on (, a ) and conve on (a, ). In the case that ψ (+) and q = the argument in Chan and Kyprianou (7) shows that a = and W is necessarily concave. In Kyprianou et al. (8) it is shown that if Π(, ) has a density on (, ) which is non-increasing and log-conve then for each q, the scale function W (q) () and its first derivative are conve beyond some finite value of. Continuity at the origin For all q, { W (q) (+) = c if σ > or (,) ( )Π(d) = if σ = and (,) ( )Π(d) <, (3) where c = a (,) Π(d). Derivative at the origin For all q, /σ if σ > W (q) (+) = if σ = and Π(, ) = (q + Π(, ))/c if σ = and Π(, ) <. (4) Behaviour at for q = As we have W() { /ψ (+) if ψ (+) > e Φ() /ψ (Φ()) if ψ (+) <. (5) When E(X ) = a number of different asymptotic behaviours may occur. For eample, if writing ψ(θ) = θφ(θ) (which is possible by the Wiener-Hopf factorization), we have φ (+) < then W() /φ (+) as. Behaviour at for q > As we have and thus there is asymptotic eponential growth. W (q) () e Φ(q) /ψ (Φ(q)) (6) 4. Methodology for new eamples For any given spectrally negative Lévy process, the scale functions are intimately connected to descending ladder height process and this forms the key to constructing new eamples. Therefore we shall briefly review the connection between scale functions and the descending ladder height process before describing the common methodology that leads to the new eamples of scale functions. For a more detailed account of this connection, the reader is referred to the books of Bertoin (996), Kyprianou (6) or Doney (7).

7 7 It is straightforward to show that the process X X := {X t X t : t }, where X t := inf s t X s, is a strong Markov process with state space [, ). Following standard theory of Markov local times (cf. Chapter IV of Bertoin (996)), it is possible to construct a local time at zero for X X which we henceforth refer to as L = {L t : t }. Its inverse process, L := {L t : t } where L t = inf{s > : L s > t}, is a (possibly killed) subordinator. Sampling X at L we recover the points of minima of X. If we define H t = X L when t < with H t = otherwise, then it is known that the process H = {H t : t } is L t a (possibly killed) subordinator. The latter is known as the descending ladder height process. Moreover, if Υ is the Lévy measure of H then Υ(, ) = e Φ() e Φ()u Π(, u)du for >, (7) see for eample Vigon (a). Further, the subordinator has a drift component if and only if σ > in which case the drift is necessarily equal to σ /. The killing rate of H is given by the constant E(X ). The starting point for the relationship between the descending ladder height process and scale functions is given by the Wiener-Hopf factorization. In Laplace form for spectrally negative Lévy processes this takes the appearance where ψ(θ) = (θ Φ())φ(θ) (8) φ(θ) = log E(e θh ), ψ(θ) = log E(e θx ), (9) and R(θ). With this form of the Wiener-Hopf factorization in mind, we appeal principally to two techniques.. We choose the process H, or equivalently φ(θ), so that the Laplace inversion of (5) may be performed. In some cases, when the process does not drift to (or equivalently ψ (+) ), it can be worked to ones advantage that, after an integration by parts, one also has that e θ W(d) = (3) φ(θ) for R(θ) >.. We choose the process H to be such that its semigroup P(H t d) is known in eplicit form and then make use of the following identity (see for eample Bertoin (996) or Kyprianou (6)), whenever X does not drift to. dt P(H t d) = W(d),, (3) Naturally, forcing a choice of descending Ladder height process, or equivalently φ, requires one to know that a spectrally negative Lévy process, X, eists whose Laplace eponent respects the factorization (8). The net Theorem provides the necessary justification. Indeed, to some etent it shows how to construct a Lévy process with a given descending ladder height process as well as a prescribed ascending ladder height process. For spectrally negative Lévy processes, the ascending ladder height process is a (killed) linear unit drift and has only a single parameter, namely Φ(), the killing rate. Hence the construction in the below theorem offers a parameter ϕ

8 8 F. HUBALEK AND A. E. KYPRIANOU which plays the role of Φ(). Since both ascending and descending ladder height process cannot both be killed one also sees the parameter restriction where κ is the killing rate of the descending ladder height process. ϕκ = (3) Theorem. Suppose that H is a subordinator, killed at rate κ, with Lévy measure which is absolutely continuous with non-increasing density and drift ζ. Suppose further that ϕ is given such that ϕκ =. Then there eists a spectrally negative Lévy process X, henceforth referred to as the parent process, such that for all, P(τ + < ) = e ϕ and whose descending ladder height process is precisely the process H. The Lévy triple (a, σ, Π) of the parent process is uniquely identified as follows. The Gaussian coefficient is given by σ = ζ. The Lévy measure is given by Π(, ) = ϕυ(, ) + dυ() d. (33) Finally a = (, ) Π(d) κ (34) if ϕ = and otherwise when ϕ > we can establish the value of a from the equation aϕ = σ ϕ + (e ϕ ϕ {> } )Π(d). (35) (,) In all cases, the Laplace eponent of the parent process is also given by for θ where φ(θ) = log E(e θh ). ψ(θ) = (θ ϕ)φ(θ) (36) Proof. The proof is reasonably self evident given the preceding account of the ladder height process. One needs only the additional information that if X is any spectrally negative Lévy process with Lévy measure Π then E(X ) = a + Π(d). (37) (, ) Moreover when Φ() > then necessarily the descending ladder height process has no killing and when the descending ladder height process is killed then Φ() =. Further, for all P(τ + < ) = e Φ(). The idea of working backwards through the Wiener-Hopf factorization as we have done above can also be found in Bertoin et al. (4) and Vigon (b). Note that it is more practical to describe the parent process in terms of the triple (σ, Π, ψ) than the triple (a, σ, Π) and we shall frequently do this in the sequel. It is also worth making an observation for later reference concerning the path variation of the process X for a given a descending ladder height process H. Corollary. Given a killed subordinator H satisfying the conditions of the previous Theorem, (i) the parent process has paths of unbounded variation if and only if Υ(, ) = or ζ >,

9 9 (ii) if Υ(, ) = λ < then the parent process necessarily decomposes in the form X t = (κ + λ ζϕ)t + ζb t S t (38) where B = {B t : t } is a Brownian motion, S = {S t : t } is an independent driftless subordinator with Lévy measure ν satisfying ν(, ) = Π(, ). Proof. The path variation of X follows directly from (33) and the fact that σ = ζ. Also using (33), the Laplace eponent of the decomposition (38) can be computed as follows with the help of an integration by parts; (κ + λ ζϕ)θ + ζθ ϕθ = (κ + Υ(, ) ζϕ)θ + ζθ ϕ ( = (θ ϕ) κ + θζ + e θ Υ(, )d θ ( e θ ) dυ d ()d e θdυ d ()d ( e θ ) dυ d ()d θ ). e θdυ d ()d This agrees with the Laplace eponent ψ(θ) = (θ ϕ)φ(θ) of the parent process constructed in Theorem. 5. The Gaussian tempered stable convolution class In this paper we introduce a new family of spectrally negative Lévy processes from which our new eamples of scale functions will emerge. We call them Gaussian tempered stable convolution, GTSC for short. When there is no Gaussian part we call the distribution a tempered stable convolution and write TSC for short. 5.. The tempered stable ladder process The starting point for the construction is a tempered stable subordinator plus a linear drift, possibly killed, that will play the role of the descending ladder height process for the GTSC parent process. Some references for tempered stable distributions and tempered stable Lévy processes are Steutel and van Harn (4), Cont and Tankov (3) and Schoutens (3). Note, that several names and origins, and many different parameterisations are used for the tempered stable distributions and processes. In our parametrization, the tempered stable subordinator involves three parameters, the stability parameter α <, the tempering parameter γ, and the scaling parameter c >. When α it is required that γ >. Furthermore we might add a linear drift with rate ζ to the process, and possibly kill the process at rate κ. The Laplace eponent of the tempered stable subordinator, and henceforth the descending ladder height process is thus taken to be φ(θ) = κ + ζθ + cγ( α)(γ α (γ + θ) α ), R(θ) > γ, (39) and the associated Lévy measure is given by Υ(d) = c α e γ d ( > ). (4) For α < the process has infinite activity. For α = the epression (39) is to be understood in a limiting sense, i.e., φ(θ) = κ + ζθ c log(γ/(γ + θ)), and the process is simply a (killed) Note that we use Γ(z) as a meromorphic function with simple poles at the non-positive integers. By the functional equation of the gamma function we have Γ( α) = α Γ( α), and Γ( α) < for α (, ).

10 F. HUBALEK AND A. E. KYPRIANOU gamma subordinator (with drift). If α < the process is a compound Poisson process with intensity parameter cγ( α)γ α and gamma distributed jumps with α degrees of freedom and eponential parameter γ >. 5.. The associated parent process We may now invoke Theorem to construct the associated GTSC process. This introduces another parameter ϕ. To be able to apply the theorem, we need a decreasing Lévy density, and thus we have to restrict the scaling parameter to α <. The theorem tells us that the Laplace eponent of the parent process is ψ(θ) = (κ ϕζ)θ + ζθ + c(θ ϕ)γ( α)(γ α (γ + θ) α ). (4) for R(θ) > γ, and obviously for α = we understand the above epression in the limiting sense so that ψ(θ) = (κ+ζϕ)θ+ζθ c(θ+ϕ)log (γ/γ + θ). It is important to recall here and throughout the remainder of the paper that κϕ =. The corresponding Lévy measure given by (ϕ + γ) (α + ) Π(d) = c ( ) α+ eγ d + c ( ) α+ eγ d (4) for <. This indicates that the jump part is the result of the independent sum of two spectrally negative tempered stable processes with stability parameters α and α +. We also note from Theorem that σ = ζ, indicating the presence of a Gaussian component. This motivates the choice of terminology Gaussian Tempered Stable Convolution. If < α < the jump component is the sum of an infinite activity negative tempered stable subordinator and an independent spectrally negative tempered stable process with infinite variation. If α = the jump component is the sum of a spectrally negative infinite variation tempered stable process with stability parameter and eponential parameter γ and the negative of a gamma subordinator with eponential parameter γ. If α < the jump part of the parent process is the independent sum tempered stable subordinator with stability parameter +α and eponential parameter γ, and an independent negative compound Poisson subordinator with jumps from a gamma distribution with α degrees of freedom and eponential parameter γ. In the etreme case α =, the parent process has negative jumps which are compound Poisson and eponentially distributed with parameter γ. 6. Evaluating GTSC scale and q-scale functions 6.. The case with rational nonzero stability parameter Theorem. Suppose α = m/n with m Z, n Z, < m < n, and gcd( m, n) =. Let us consider the GTSC process with Laplace eponent ψ(θ) = (θ ϕ)[κ + ζθ + cγ( α)(γ α (γ + θ) α )]. (43) Let m + = ma(m, ) and m = ma( m, ). Then the polynomial f q (z) = (z n γ ϕ) [ (κ + ζγ + cγ( α)γ α )z m + ζz n+m cγ( α)z m+ ) ] qz m (44) has at least one real root. Let l denote the number of distinct roots, r,..., r l the roots, arranged such that r is the largest real root, and µ,..., µ l their multiplicities. Let A kj for k =,...,l, j =,...,µ k denote the coefficients in the partial fraction decomposition z l m f q (z) = k= µ k j= A kj. (45) (z r k ) j+

11 Then we have and Proof. We have W (q) () = e γ l k= Φ(q) = r n γ (46) µ k j= for z Z with arg(z) < π/n, and conversely, A kj j! (j+)/n E (j) (r n, k n ). (47) n f q (z) = z m (ψ(z n γ) q) (48) ψ(θ) q m (θ + γ) n = f q ((θ + γ) n) (49) for θ C\(, γ]) such that the denominator on the left hand side does not vanish. We know from the theory of scale functions that θ = Φ(q) is a real root of ψ(θ) q =, and there is no larger real root of that equation. Thus z = (Φ(q) + γ) n is a real root of f q (z) = and there is no larger real root of that equation. This shows (46). The left hand side of (49) is the Laplace transform of W (q) (), and we know it is an analytic function for R(θ) > Φ(q). Using the partial fraction decomposition we get e θ W (q) ()d = l k= µ k j= A kj, (5) ((γ + θ) n r k ) j+ and this equation is valid for any q and R(θ) > Φ(q). We recognize the Laplace transform of eponentially tilted derivatives of Mittag-Leffler functions on the right hand side, see for eample (Jacob, 5, Prop.7..9, p.359), and obtain (47). Generally speaking this is rather an awkward formula to work with analytically. The main strength of the epression lies with it being a simple matter to program into a package such as MATLAB or Mathematica. It is instructive to revisit the (known) analytical properties of W() and W (q) () listed in section 3 and check them directly for our eplicit epression as this will yield in most cases more detail. Smoothness: Since the Mittag-Leffler functions are entire functions, we see that W (q) (), which is a priori defined for (, ), admits an analytic continuation to C \ (, ], and thus W (q) is C on (, ). Behaviour at zero: We obtain the behaviour at zero from the power series epansion of the derivatives of the Mittag-Leffler function. Eploiting algebraic relations of the r k and A kj several terms cancel and we obtain for ζ > W (q) () ζ, W (q) () ζ ( ). (5) A similar argument applies for ζ =. If < α < we get α W (q) () cγ( α)γ( + α), α W (q) (), ( ). (5) cγ( α)γ(α)

12 F. HUBALEK AND A. E. KYPRIANOU In the case that < α < and ζ =. we get again by cancellation of terms and limiting behaviour at zero of Mittag-Leffler functions W (q) () κ + cγ α Γ( α), W (q) () c α (κ + cγ α, ( ). (53) Γ( α)) Note, another way to establish these last two results is by Karamata s Tauberian Theorem and the Monotone Density Theorem respectively from the asymptotics of ψ(θ) as θ. This approach shows, that the asymptotics above also hold for irrational α (, ). Behaviour at infinity: The behaviour at infinity can be obtained from the asymptotics of the Mittag-Leffler functions at infinity. The dominating contribution comes from the term with r. It is useful to recall that f q (z) = ψ(z n γ) q and thus f q (z) = nzn ψ (z n γ). Let us consider first q =. Note, that we have ψ (+) = κ ϕ(ζ + cγ α Γ( α)). Suppose ψ (+) >, which implies that ϕ =. Then r = γ n is a single root of f (z). After some elementary simplifications we obtain W() ( ) (54) κ in agreement with the theory. Suppose now ψ (+) <. Then r = (ϕ + γ) n is a simple root of f (z). As above, using the asymptotics of the Mittag-Leffler functions at infinity, we see that the term with r dominates all other terms. After elementary simplifications we obtain W() e ϕ κ + ζϕ + cγ( α)(γ α (γ + ϕ) α ) ( ). (55) Consider the third case, ψ (+) =. This can only happen when κ = and ϕ =. Then r = γ n is a double root of f (z). The asymptotics show, that the contribution from r dominates and we obtain W() ζ + cγ α ( ). (56) Γ( α) Finally in the case q >, little more can be said than (6). 6.. The case with general nonzero stability parameter α (, ) When we take q = and ζ = and α (, )\{} (without the restriction of begin a rational number) one may obtain much cleaner epressions for the scale function than the formulation in Theorem. Theorem 3. Suppose α (, ) \ {} and consider a TSC process (without Gaussian component) and Laplace eponent Then If < α < then if < α <, then W() = W() = cγ( α) ψ(θ) = (θ ϕ)[κ + cγ( α)(γ α (γ + θ) α )]. (57) eϕ e ϕ κ + cγ( α)γ α + cγ( α)eϕ (κ + cγ( α)γ α ) ( κ + cγ( α)γ e (γ+ϕ)y y α α y α ) E α,α dy, (58) cγ( α) ( ) cγ( α)y e (γ+ϕ)y y α α E α, α κ + cγ( α)γ α dy. (59)

13 3 Proof. This follows from the known Laplace transforms e θ α E α,α (λ α )d = θ α λ (6) and e θ λ α E α, α (λ α )d = λ, (6) λ θα valid for α > resp. α < together with the the well-known rules for Laplace transforms concerning, primitives and tilting. With these closed-form epressions, the eplicit asymptotics discussed after Theorem are immediate by inspection. Note in the particular case that κ = cγ( α)γ α the epression for the scale function in the case < α < reduces simply to W() = κ γ α Γ(α) yα e γy dy. (6) Note that although the Laplace inversions involved in the proof of Theorem 3 are straightforward, it is not clear that the resulting epressions for W are scale functions without the presence of Theorem The case with α = : inverse Gaussian descending ladder height process We will assume now that ζ =, ϕ =, κ =. Substituting α /, c δ/ π, γ γ / reveals that the tempered stable subordinator we have been working with is in fact the familiar inverse Gaussian IG(δ, γ)-subordinator when α = /. It turns out for this case that we also get some cleaner epressions than the one given in Theorem. When the descending ladder height process is the IG(δ, γ) subordinator, the jump part of the parent process is a superposition of a negative IG subordinator and a spectrally negative tempered stable process with stability parameter 3/. Following the general approach in Theorem we relate ψ(θ) to a polynomial, which, in the IG parametrization of this section, is f q (z) = ψ(z γ /) q. The first simplification is, that f q (z) is a polynomial of degree 4 if ζ >, and degree 3, if ζ =. Thus completely eplicit, elementary (and somewhat lengthy) epressions of its roots can be given in terms of radicals by the formulæ of Ferrari resp. Cardano. The second simplification concerns the Mittag-Leffler functions. As α = / the Mittag-Leffler function simplifies to an epression involving the more familiar (complementary) error function. We find that f (z) has a positive double a negative single root, and after some elementary simplifications, W() = [ ( ( + γ )erfc γ ) / + γ ] /πe γ. (63) δγ Actually, in this special case the result can be verified faster by calculating the Laplace transform of (63). As a side remark, we note, that W () = 3/ e γ /( πδ) and thus we see, that W is indeed concave, as we already know from the theory. Let us now consider q > and put q = 6 7 δγ3. (64)

14 4 F. HUBALEK AND A. E. KYPRIANOU If < q < q then we have three simple real roots r, r, r 3, that can be obtained by Cardano s formula, and [ ] W (q) () = e γ r e r f q(r ) erfc( r r e r ) + f q(r ) erfc( r r 3 e r 3 ) + f q(r 3 ) erfc( r 3 ). (65) If q > q we have a simple real root r and two comple conjugate simple roots r and r 3. The same formulas hold as in the previous case. Alternatively, we could in the comple case write [ }] W (q) () = e γ r e r f q(r ) erfc( r {r e r ) + R f q(r ) erfc( r ). (66) If q = q then we have a simple positive and a double negative root, and after some easy simplifications, [ ( W (q) () = 6γ 36δγ π e γ + 5e 8 9 γ erfc 5γ ) 3 e 4 9 γ (5 + γ )erfc ( γ 3 )]. (67) 6.4. The stable case Because of their definition via their Laplace transform, it is immediate that scale functions are continuous in the parameters of the parent process. Taking γ in the Laplace eponent of the tempered stable ladder height process shows us that the parent process is an α + stable process when κ = ϕ = ζ = and c = /Γ( α). Taking limits as γ in the epression for the scale function in the first part of Theorem 3 gives us the known result W() = α /Γ(α + ) The case α = : gamma descending ladder height process When α = we can only treat the special case with q =, κ =, ζ = and ϕ =. Theorem 4. W() = where P(a, ) is the regularized lower incomplete gamma function P(ct, γ)dt, (68) Proof. We may appeal to the second method in (3) using resolvents to note that the scale function density satisfies W () = c e γ ϕ( log(γ)) (69) Eplicitly r = γ r = γ r 3 = γ 3 + δγ 3 3q 8δ 3 γ 3 +7δ q+3 3 7δ 4 q 6δ 5 γ 3 q 3 (+i 3)δγ 3 3q 8δ 3 γ 3 +7δ q+3 3 7δ 4 q 6δ 5 γ 3 q 3 ( i 3)δγ 3 3q 8δ 3 γ 3 +7δ q+3 3 7δ 4 q 6δ 5 γ 3 q If q = δγ 3 / those formulae simplify, as one root is then zero. + q 3 8δ 3 γ 3 +7δ q+3 3 7δ 4 q 6δ 5 γ 3 q 3 δ q ( i 3) 3 8δ 3 γ 3 +7δ q+3 3 7δ 4 q 6δ 5 γ 3 q 6 δ q (+i 3) 3 8δ 3 γ 3 +7δ q+3 3 7δ 4 q 6δ 5 γ 3 q 6 δ.

15 5 where ϕ(θ) = e θ d, (7) Γ() which is the Laplace transform of the reciprocal gamma function. 3 As X has infinite variation, we have W(+) =. Integrating (69) yields the result Remarks on the cases α = / and α = The case α = / is again covered in Theorem above, but several simplifications occur. The jump part of the ladder process is in this case a compound Poisson process with jumps from the χ -distribution. Again, we can avoid Mittag-Leffler functions and use the complimentary error function instead. With regard to the roots of the polynomial equation we may say the following. If ζ > the polynomial f q (z) has degree 5 and its roots typically cannot be epressed in terms of radicals. If ζ = the polynomial is of degree 3 and Cardano s formula can be used. If in addition κ = and ϕ = then we have r = γ a double and r = γ a single root. We can proceed as is Section 6.3 to obtain some simple closed form epressions for W() and W (q) (). When α = the jump part of the parent process is a compound Poisson process with negative eponential jumps, and the ladder process is a compound Poisson process with eponential jumps. The Laplace transform of the scale functions are rational, and the scale functions can be obtained directly from the partial fraction decomposition. The results are implicitly contained in Mordecki and Lewis (5) and Kou and Wang (3) and we eclude the calculations here. 7. Numerical illustrations Figure contains graphs of the GTSC scale function W() for stability parameters α = /4, /3, /, /3, 3/4 in si cases. Figure contains graphs of the GTSC q-scale functions W (q) () with q =, but otherwise for the same parameters. Below we introduce si cases considered in terms of the classification of the parent process. Case A, κ =, ϕ =, ζ =, c =, γ = : The parent process is oscillating, has no diffusion part and infinite variation jumps. The ladder process is an infinite activity pure jump subordinator, has no linear drift, and is not killed. Case B, κ =, ϕ =, ζ =, c =, γ = : The parent process drifts to +, has no diffusion part and infinite variation jumps. The ladder process is an infinite activity pure jump subordinator killed at rate κ =. Case C, κ =, ϕ =, ζ =, c =, γ = : The parent process drifts to, has no diffusion part and infinite variation jumps. The ladder process is an infinite activity pure jump subordinator. Case D, κ =, ϕ =, ζ =, c =, γ = : The parent process is oscillating, has no linear drift, but a diffusion part and infinite variation jumps. The ladder process is an infinite activity pure jump subordinator plus a linear drift. Case E, κ =, ϕ =, ζ =, c =, γ = : The parent process drifts to +, has a Gaussian component and infinite variation jumps. The ladder process is an infinite activity subordinator plus linear drift, killed at unit rate. 3 The value at θ = is called the Fransén-Robertson constant and denoted by ϕ() = F. A lot of analytical and numerical material on ϕ(θ) can be found in Fransén and Wrigge (984). That paper contains also several references to related work going back to Paley, Wiener, Hardy and Ramanujan. Note, that the reciprocal gamma distribution in Fransén and Wrigge (984) is not to be confused with the inverse gamma distribution, which is often called reciprocal gamma distribution as well.

16 6 F. HUBALEK AND A. E. KYPRIANOU Case F, κ =, ϕ =, ζ =, c =, γ = : The parent process drifts to, has a diffusion part and infinite variation jumps. The ladder process is an infinite activity pure jump subordinator plus a linear drift. Now let us discuss the graphs in view of the theoretical properties listed in Section 3. Smoothness: All graphs look smooth, as the theory predicts. Concavity and conveity: For q = we know in cases A, B, D and E the graph is predicted to be concave. In cases C and F the graph is predicted to be conve-concave. Indeed this is the appearance in Figure. For q = we know, that all graphs are concave-conve. We observe this behaviour in Figure. In some cases, it is necessary to inspect the shape of the graph more closely in the neighbourhood of the origin in order to see concavity. Behaviour at zero: For all cases we observe W() = in Figure and W (q) () = in Figure, in agreement with formula (3). For the cases A C (no diffusion part) we observe W () = +, for the cases D E (nonzero diffusion part) we observe W () = in Figure. We observe the same behaviour, W (q) () = resp. W (q) () =, in Figure. This is in agreement with formula (4). Behaviour at infinity: Firstly let us consider the case q =. In cases A and D we have ψ (+) = and thus, according to (56), asymptotic linear growth as. In cases B and E we have κ > and hence ψ (+) > so that W() /κ as. In cases C and F we have ϕ > and hence ψ () < and so there is eponential growth of W() as according to the second case in formula (5). This corresponds to what is observed in the graphs in Figure. (The behaviour for becomes more prominent when plotting < <.) For the case q >, all graphs in Figure appear ehibit eponential growth as, in agreement with formula (6). 8.. Etensions and ramifications 8. Additional Remarks Our methods apply also to the integrated q-scale functions, see (Kyprianou, 6, Theorem 8., p.4) Z (q) () = + q W (q) (y)dy, (7) for the Gaussian tempered stable convolution class with rational stability parameter. Also we could treat more complicated Gaussian tempered stable convolutions with φ(θ) = κ + ζθ + n i= c i Γ( α i )(γ αi i (γ i + θ) αi ) (7) with α,..., α n rational numbers in [, ). The key observation is simply, that any Laplace transform, that is of the form f(θ α ) with f a rational function and α some number, can be inverted in terms of the partial fraction decomposition of f using Mittag-Leffler functions and its derivatives. In other very recent work, Kyprianou and Rivero (7) have looked at families of scale functions which are the result of choosing the descending ladder height process to have a Laplace eponent which is a special or complete Bernstein function.

17 q=, Case A κ=, ϕ=, ζ=, c=, γ= α=/4 α=/3 α=/ α=/3 α=3/ q=, Case D κ=, ϕ=, ζ=, c=, γ= α=/4 α=/3 α=/ α=/3 α=3/4..6 W().8 W() W() q=, Case B κ=, ϕ=, ζ=, c=, γ= α=/4 α=/3 α=/ α=/3 α=3/4 W() q=, Case E κ=, ϕ=, ζ=, c=, γ= α=/4 α=/3 α=/ α=/3 α=3/ W() q=, Case C κ=, ϕ=, ζ=, c=, γ= α=/4 α=/3 α=/ α=/3 α=3/4 W() q=, Case F κ=, ϕ=, ζ=, c=, γ= α=/4 α=/3 α=/ α=/3 α=3/ Figure : Scale functions W() for the GTSC class with < α <

18 8 F. HUBALEK AND A. E. KYPRIANOU W (q) () q=, Case A κ=, ϕ=, ζ=, c=, γ= α=/4 α=/3 α=/ α=/3 α=3/4 W (q) () q=, Case D κ=, ϕ=, ζ=, c=, γ= α=/4 α=/3 α=/ α=/3 α=3/ W (q) () q=, Case B κ=, ϕ=, ζ=, c=, γ= α=/4 α=/3 α=/ α=/3 α=3/4 W (q) () q=, Case E κ=, ϕ=, ζ=, c=, γ= α=/4 α=/3 α=/ α=/3 α=3/ W (q) () q=, Case C κ=, ϕ=, ζ=, c=, γ= α=/4 α=/3 α=/ α=/3 α=3/4 W (q) () q=, Case F κ=, ϕ=, ζ=, c=, γ= α=/4 α=/3 α=/ α=/3 α=3/ Figure : Scale functions W (q) () for the GTSC class with q = and < α <

19 9 8.. Numerical inversion of the Laplace transform While our formulas provide closed form epressions for rational α in terms of partial fraction decompositions and Mittag-Leffler functions, those epression can be rather lengthy, involve the zeros of polynomials of high degree, and, eperience shows them to be numerically problematic when α = m/n unless n is quite small. Alternatively one may always consider numerical inversion of the Laplace transform of W (q), which is etremely simple for the GTSC class. Moreover this approach works also for irrational α. We use the Bromwich inversion integral for the Laplace transform, W (q) () = πi r+i r i e θ dθ ψ(θ) q (73) where r > Φ(q). Let us consider the integrand for θ with R(θ) r. The term e θ is bounded. If ζ > we have ψ(θ) ζθ, if ζ = and < α < we have ψ(θ) cγ( α)θ α+ and (73) is a proper Lebesgue integral. Otherwise it must be interpreted as a principal value integral. Some remarks how to handle this situation and how to improve the numerics and a piece of Mathematica code, that evaluates the scale function by (straightforward, and numerically unsophisticated) integration, are in the etended report. For further material on the numerical evaluation of scale function see also (Surya, 7, Chapter7, p.ff) Eamples of applications Let us briefly conclude by re-enforcing the genuine importance of establishing eplicit eamples of scale functions from the point of view of modelling by giving three classical eamples from the theory of applied probability. It will suffice to consider the role of -scale functions. We shall do so in the light of some of the eplicit eamples above. Eample. Recent literature suggests that a more modern approach to the theory of ruin should involve considering, in place of the classical Cramér-Lundberg process, a spectrally negative Lévy process which drifts to for the risk process. As noted earlier, for any such given Lévy process with Laplace eponent ψ, the quantity ψ (+)W() is the probability of ruin when the Lévy process is issued from at time. Consider the following parent process where {S (i) t by X t = (κ + λ)t S () t S () t, t (74) : t }, i =,, are subordinators whose Lévy measures, ν (i), i =,, are identified ν () (d) = λ γν+ Γ(ν) ν e γ d and ν () (d) = λ γν Γ(ν) ( ν)ν e γ d (75) for >. Following the interpretation of Lévy risk processes in Klüppelberg and Kyprianou (6) and Huzak et al. (4) we may think of the parent process as the result of premiums collected at rate λ+κ to offset the result of two competing claim processes. The first has claims arriving at a compound Poisson rate λ which are gamma distributed and the second has an infinite intensity of small claims whose aggregate behaviour is similar to that of a stable subordinator with inde ν. Moreover, these claims are interlaced with larger claims of finite intensity, the tail of whose distribution are similar to those of gamma distributed jumps.

20 F. HUBALEK AND A. E. KYPRIANOU Recalling the discussion in Section 5. we have from the second part of Theorem 3 that the probability of ruin from an initial reserve > is conveniently given by λ + κ ργν λ + κ y ν e γy E ν,ν (ργ ν y ν )dy (76) Eample. If instead we consider the parent process in the previous eample reflected in its supremum, then we are looking at the workload of a so called M/Π/ queue (cf Dube et al. (4)). That is to say, a queue in which work is processed with a constant rate, say µ, and which arrives in packets according to a Poisson point process such that a job of size > arrives in the interval (t, t + dt) with probability Π(d)dt + o(dt) where ( )Π(d) <. For the case at hand µ = λ + κ and Π(d) = ν () (d) + ν () (d) and there is again the possibility of interpreting the incoming work as the result of two competing mechanisms. In this case, the stationary distribution of the workload is given by one minus the ruin probability given above. Eample 3. In the very recent work of Loeffen (7), the classical dividend control problem of De Finetti s dividend problem has been considered in the contet of the driving risk process being a general spectrally negative Lévy process denoted as usual X. Here the objective is to find the optimal strategy and value to the following control problem v() = sup π ( σ π E e qt dl π t where q >, σ π = inf{t > : X t L π t < } and Lπ = {L π t : t } is the process of dividends paid out associated with the strategy π and the supremum is taken over all strategies π such that L π is a non-decreasing, left-continuous adapted process which starts at zero and for which ruin cannot occur by a dividend pay-out. In Loeffen (7), the author proves the remarkable fact that, under the assumption that the dual of the underlying Lévy process has a completely monotone density, for each q >, the associated q-scale function, W (q), has a first derivative which is strictly conve on (, ). Moreover, the point a := inf{a : W (q) (a) W (q) () for all } is the level of the barrier which characterizes the optimal strategy. The latter being that dividends are paid out in such a way that the aggregate process has the dynamics of the underlying Lévy process reflected at the barrier a. The value function of this strategy is given by W (q) () for a W v() = (q) (a) (78) a + W(q) (a) for > a W (q) (a) If we are to take any of the parent processes mentioned in this paper as the underlying spectrally negative Lévy process, then it is immediately obvious from (4) that they all have the property that their dual has a Lévy density which is completely monotone. In such a case one may observe the value a graphically (see for eample Figure ) as well as compute it numerically to a reasonable degree of eactness with the help of software such as Mathematica. Acknowledgments This work was initiated when both authors were attending the mini-workshop Lévy processes and Related Topics in Modelling in Oberwolfach. We are grateful to the organizers of this workshop and MFO for the opportunity it has provided us. We would also like to thank Ronnie Loeffen for hinting that previously unnoticed scale functions are to be found in the references of Furrer (998), Asmussen () and Abate and Whitt (999). Thanks also go to Victor Rivero who commented on earlier versions of this paper. ) (77)

21 References [] J. Abate and W. Whitt. Eplicit M/G/ waiting-time distributions for a class of long-taile service-time distributions. Operations Research Letters, 5:5 3, 999. [] L Alili and A.E. Kyprianou. Some remarks on first passage of Lévy processes, the American put and pasting principles. The Annals of Applied Probability, 5(3):6 8, 5. [3] S. Asmussen. Ruin probabilities. World Scientific Publishing,. [4] F. Avram, T. Chan, and M. Usabel. On the valuation of constant barrier options under spectrally one-sided eponential Lévy models and Carr s approimation for American puts. Stochastic Processes and their Applications, :75 7,. [5] F. Avram, A.E. Kyprianou, and M. R. Pistorius. Eit problems for spectrally negative Lévy processes and applications to (Canadized) Russian options. The Annals of Applied Probability, 4():5 38, 4. [6] F. Avram, Z. Palmowski, and M. R. Pistorius. On the optimal dividend problem for a spectrally negative Lévy process. Annals of Applied Probability, 7:56 8, 7. [7] R. Bekker, O. Boma, and O. Kella. Queues with delays in two-state strategies and Lévy input. EURANDOM report 7 8, 7. Submitted for publication. [8] J. Bertoin. Lévy processes, volume of Cambridge Tracts in Mathematics. Cambridge University Press, Cambridge, 996. [9] J. Bertoin. Eponential decay and ergodicity of completely asymmetric Lévy processes in a finite interval. The Annals of Applied Probability, 7():56 69, 997. [] J. Bertoin, B. Roynette, and M. Yor. Some connections between (sub)critical branching mechanisms and bernstein functions. 4. arxiv:math/43v [math.pr]. [] N. H. Bingham. Continuous branching processes and spectral positivity. Stochastic Processes and their Applications, 4(3):7 4, 976. [] O. J. Boma and J.W. Cohen. The M/G/ queue with heavy-tailed service-time distribution. IEEE Journal on Selected Areas in Communications, 6: , 998. [3] T. Chan and A.E. Kyprianou. Smoothness of scale functions for spectrally negative Lévy processes. Preprint, 7. [4] L. Chaumont. Sur certains processus de Lévy conditionnés à rester positifs. Stochastics and stochastics reports, 47:, 994. [5] L. Chaumont. Conditionings and path decompositions for Lévy processes. Stochastics and stochastics reports, 64:39 54, 996. [6] L. Chaumont and M.E. Caballero. Conditioned stable Lévy processes and the Lamperti representation. Journal of Applied Probability, 43: , 6. [7] L. Chaumont, A. E. Kyprianou, and J.C. Pardo. Wiener-Hopf factorization and some eplicit identities associated with positive self-similar Markov processes. Preprint, 7. [8] S.N. Chiu and C. Yin. Passage times for a spectrally negative Lévy process with applications to risk theory. Bernoulli, (3):5 5, 5.

22 F. HUBALEK AND A. E. KYPRIANOU [9] R. Cont and P. Tankov. Financial modelling with Jump Processes. Chapman & Hall, CRC Press., 3. [] A.D. Doney and A.E. Kyprianou. Overshoots and undershoots of Lévy processes. The Annals of Applied Probability, 6():9 6, 6. [] R. A. Doney. Some ecursion calculations for spectrally one-sided Lévy processes. In Séminaire de Probabilités XXXVIII, volume 857 of Lecture Notes in Math., pages 5 5. Springer, Berlin, 5a. [] R.A. Doney. Hitting probabilities for spectrally positive Lévy processes. Journal of the London Mathematical Society. Second Series, 44(3): , 99. [3] R.A. Doney. Fluctuation theory for Lévy processes. In École d été de probabilités de Saint-Flour, XXXV 5, volume 897 of Lecture Notes in Math. Springer, Berlin, 7. [4] R.A. Doney. Some ecursion calculations for spectrally one-sided Lévy processes. Séminaire de Probabilités, XXXVIII:5 5, 5b. [5] P. Dube, F. Guillemin, and R. R. Mazumdar. Scale functions of Lévy processes and busy periods of finite-capacity M/GI/ queues. Journal of Applied Probability, 4(4):45 56, 4. [6] D. J. Emery. Eit problem for a spectrally positive process. Advances in Applied Probability, 5: 498 5, 973. ISSN [7] A. Fransén and S. Wrigge. Calculation of the moments and the moment generating function for the reciprocal gamma distribution. Mathematics of Computation, 4(66):6 66, 984. [8] H. Furrer. Risk processes perturbed by α-stable Lévy motion. Scandinavian. Actuarial Journal, :59 74, 998. [9] B. Hilberink and L. C. G. Rogers. Optimal capital structure and endogenous default. Finance and Stochastics, 6():37 63,. [3] M. Huzak, M. Perman, H. Šikić, and Z. Vondraček. Ruin probabilities for competing claim processes. Journal of Applied Probability, 4:679 69, 4. [3] N. Jacob. Pseudo differential operators and Markov processes. Vol.III. Imperial College Press, London, 5. Markov processes and applications. [3] C. Klüppelberg and A.E. Kyprianou. On etreme ruinous behaviour of Lévy insurance risk processes, 6. [33] C. Klüppelberg, A. E. Kyprianou, and R. A. Maller. Ruin probabilities and overshoots for general Lévy insurance risk processes. The Annals of Applied Probability, 4(4):766 8, 4. [34] V. S. Korolyuk. Boundary problems for a compound Poisson process. Theory Probab. Appl., 9: 4, 974. [35] V. S. Korolyuk. On ruin problems for a compound Poisson process. Theory Probab. Appl., : , 975. [36] S. Kou and H. Wang. First passage times of a jump diffusion process. Advances of Applied Probability, 35:54 53, 3.