Potentials of stable processes

Transcription

1 Potentials of stable processes A. E. Kyprianou A. R. Watson 5th December 23 arxiv:32.222v [math.pr] 4 Dec 23 Abstract. For a stable process, we give an explicit formula for the potential measure of the process killed outside a bounded interval and the joint law of the overshoot, undershoot and undershoot from the maximum at exit from a bounded interval. We obtain the equivalent quantities for a stable process reflected in its infimum. The results are obtained by exploiting a simple connection with the Lamperti representation and exit problems of stable processes. AMS 2 subject classifications: 6G52, 6G8, 6G5. Keywords and phrases: Lévy processes, stable processes, reflected stable processes, hitting times, positive self-similar Markov processes, Lamperti representation, potential measures, resolvent measures Introduction and results For a Lévy process X, the measure U A (x, dy) = E x ½[X t dy]½[ s t : X s A] dt, called the potential (or resolvent) measure of X killed outside A, is a quantity of great interest, and is related to exit problems. The main cases where the potential measure can be computed explicitly are as follows. If X is a Lévy process with known Wiener Hopf factors, it can be obtained when A is half-line or R; see [2, Theorem VI.2]. When X is a totally asymmetric Lévy process with known scale functions, it can be obtained for A a bounded interval, a half-line or R; see [9, Section 8.4]. Finally, [] details a technique to obtain a potential measure for a reflected Lévy process killed outside a bounded interval from the same quantity for the unreflected process. University of Bath, UK. a.kyprianou@bath.ac.uk CIMAT, Mexico. alexander.watson@cimat.mx. This work was begun while the second author was a Ph.D. student at the University of Bath, and completed during a position at CIMAT.

2 In this note, we consider the case where X is a stable process and A is a bounded interval. We compute the measure U [,], from which U A may be obtained for any bounded interval A via spatial homogeneity and scaling; and from this we compute the joint law at first exit of [, ] of the overshoot, undershoot and undershoot from the maximum. Furthermore, we give the potential measure and triple law also for the process reflected in its infimum. The potential measure has been already been computed when X is symmetric; see Blumenthal et al. [4, Corollary 4] and references therein, as well as Baurdoux []. We extend these results to asymmetric stable processes with jumps on both sides. The essential observation is that a potential for X with killing outside a bounded interval may be converted into a potential for the Lamperti transform of X, say ξ, with killing outside a half-line. To compute this potential in a half-line, it is enough to know the killing rate of ξ and the solution of certain exit problems for X. The results for the reflected process are obtained via the work of Baurdoux []. We now give our results. Some facts we will rely on are summarised in section 2, and proofs are given in section 3. We work with the (strictly) stable process X with scaling parameter α and positivity parameter ρ, which is defined as follows. For (α, ρ) in the set A = {(α, ρ) : α (, ), ρ (, )} {(α, ρ) = (, /2)} {(α, ρ) : α (, 2), ρ ( /α, /α)}, let X, with probability laws (P x ) x R, be the Lévy process with characteristic exponent { c θ α ( iβ tan πα Ψ(θ) = 2 sgn θ) α (, 2) \ {}, c θ α =, θ R, where c = cos(πα(ρ /2)) and β = tan(πα(ρ /2))/ tan(πα/2). This Lévy process has absolutely continuous Lévy measure with density where c x (α) ½[x > ] c x (α) ½[x < ], x R, c = Γ(α ) Γ(αρ)Γ( αρ), c = Γ(α ) Γ(αˆρ)Γ( αˆρ) and ˆρ = ρ. The parameter set A and the characteristic exponent Ψ represent, up a multiplicative constant in Ψ, all (strictly) stable processes which jump in both directions, except for Brownian motion and the symmetric Cauchy processes with non-zero drift. The normalisation is the same as that in [8], and when X is symmetric, that is when ρ = /2, the normalisation agrees with that of [4]. We remark that the quantities we are interested in can also be derived in cases of one-sided jumps: either X is a subordinator, in which case the results are trivial, or X is a spectrally one-sided Lévy process, in which case the potentials in question may be assembled using the theory of scale functions; see [9, Theorem 8.7 and Exercise 8.2]. 2

3 The choice α and ρ as parameters is explained as follows. X satisfies the α-scaling property, that under P x, the law of (cx tc α) t is P cx, () for all x R, c >. The second parameter satisfies ρ = P (X t > ). Having defined the stable process, we proceed to our results. Let σ [,] = inf{t : X t / [, ]}, and define the killed potential measure and potential density σ [,] U (x, dy) := U [,] (x, dy) = E x ½ [X t dy] dt = u (x, y) dy, provided the density u exists. Theorem. For < x, y <, (x y)α u (x, y) = (y x)α y(x) xy x(y) yx s αρ (s ) αˆρ ds, y < x, s αˆρ (s ) αρ ds, x < y. When X is symmetric, this reduces, by spatial homogeneity and scaling of X, and substituting in the integral, to [4, Corollary 4]. With very little extra work, Theorem yields an apparently stronger result. Let and write τ = inf{t : X t < }; X t = sup X s, t, s t τ E x ½[X t dy, X t dz] dt = u(x, y, z) dy dz, if the right-hand side exists. Then we have the following. Corollary 2. For x >, y [, z), z > x, u(x, y, z) = Proof. Rescaling, we obtain τ E x xαˆρ y αρ (z x)αρ (z y) αˆρ z α dy dz. (2) ½[X t dy, X t z] dt = z α u (x/z, y/z), and the density is found by differentiating the right-hand side in z. 3

4 From this density, one may recover the following hitting distribution, which originally appeared in Kyprianou et al. [, Corollary 5]. Let τ = inf{t : X t > }. Corollary 3. For u [, x), v (u, ], y, P x ( X τ du, X τ dv, X τ = Γ(α ) Γ(αˆρ)Γ( αρ) dy, τ < τ ) x αˆρ ( v) αρ ( u x) αρ (v u) αˆρ ( u) α (v y) α du dv dy. (3) Proof. Following the proof of [2, Proposition III.2], one may show that the left-hand side of (3) is equal to u(x, v, u)π(v y), where π is the Lévy density of X. Remark 4. The proof of Corollary 3 suggests an alternative derivation of Theorem. Since the identity (3) is already known, one may deduce u(x, y, z) from it by following the proof backwards. The potential u (x, y) without X may then be obtained via integration. However, in section 3 we offer instead a self-contained proof based on well-known hitting distributions for the stable process. Now let Y denote the stable process X reflected in its infimum, that is, Y t = X t X t, t, where X t = inf{x s, s t} for t. Y is a self-similar Markov process. Let T = inf{t > : Y t > } denote the first passage time of Y above the level, and define T R (x, dy) = E x ½ [Y t dy] dt = r (x, y) dy, provided that the density on the right-hand side exists. We may then use the results of Baurdoux [] to find r. Note that, as Y is self-similar, R suffices to deduce the potential of Y killed at first passage above any level. Theorem 5. For < y <, Hence, for < x, y <, r (x, y) = r (, y) = Γ(α) yαρ ( y) αˆρ. [(x y) α y(x) xy y αρ ( y) αˆρ x s αρ (s ) αˆρ ds ] t αρ ( t) αˆρ dt, y < x, x(y) [(y x) α yx s αˆρ (s ) αρ ds x y αρ ( y) αˆρ ] t αρ ( t) αˆρ dt, x < y. 4

5 Writing E x ½[Y t dy, Y t dz] dt = r(x, y, z) dy dz, where Y t is the supremum of Y up to time t, we obtain the following corollary, much as we had for X. Corollary 6. For y (, z), z, r(, y, z) = αˆρ Γ(α) yαρ (z y) αˆρ, and for x >, y (, z), z x, [ r(x, y, z) = yαρ (z y) αˆρ x αˆρ (z x) αρ z α x ] z αˆρ t αρ ( t) αˆρ dt We also have the following corollary, which is the analogue of Corollary 3. Corollary 7. For u (, ], v (u, ), y, P ( Y T du, Y T dv, Y T = dy) α αˆρ Γ(αρ)Γ( αρ) and for x, u [, x), v (u, ), y, P x ( Y T Γ(α ) = Γ(αˆρ)Γ( αρ) du, Y T dv, Y T dy) ( v) αρ (v u) αˆρ (v y) α [ x αˆρ ( u x) αρ ( u) α αˆρ x u ( v) αρ (v u) αˆρ (v y) α du dv dy, ] t αρ ( t) αˆρ dt du dv dy. The marginal in dv dy appears in Baurdoux [, Corollary 3.5] for the case where X is symmetric and x =. The marginal in dy is given in Kyprianou [] for the process reflected in the supremum; this corresponds to swapping ρ and ˆρ. However, unless x =, it appears to be difficult to integrate in Corollary 7 and obtain the expression found in []. Finally, one may integrate in Theorem 5 and obtain the expected first passage time for the reflected process. Corollary 8. For x, [ E x [T ] = x ] x αˆρ ( x) αρ αˆρ t αρ ( t) αˆρ dt. Γ(α ) In particular, E [T ] = Γ(α) Γ(αρ)Γ(αˆρ ). Γ(α ) 5

6 2 The Lamperti representation We will calculate potentials related to X by appealing to the Lamperti transform [2, 4]. Recall that a process Y with probability measures (P x ) x> is a positive self-similar Markov process (pssmp) if it is a standard Markov process (in the sense of [3]) with state space [, ) which has zero as an absorbing state and satisfies the scaling property: under P x, the law of (cy tc α) t is P cx, for all x, c >. The Lamperti transform gives a correspondence between pssmps and killed Lévy processes, as follows. Let S(t) = t (Y u) α du; this process is continuous and strictly increasing until Y reaches zero. Let T be its inverse. Then, the process ξ s = log Y T (s), s is a Lévy process, possibly killed at an independent exponential time, and termed the Lamperti transform of Y. Note that ξ = log x when Y = x, and one may easily see from the definition of S that e αξ S(t) ds(t) = dt. A simple example of the Lamperti transform in action is given by considering the process X. Define τ = inf{t : X t < }, and let P x (X t ) = P x (X t, t < τ ), t, x >. The process X with laws (P x) x> is a pssmp. Caballero and Chaumont [5] gives explicitly the generator of its Lamperti transform, whose laws we denote (P y) y R, finding in particular that it is killed at rate 3 Proofs q := c /α = Γ(α) Γ(αˆρ)Γ( αˆρ). (4) To avoid the proliferation of symbols, we generally distinguish processes only by the measures associated with them; the exception is that self-similar processes will be distinguished from processes obtained by Lamperti transform. Thus, the time τ = inf{t : X t > } always refers to the canonical process of the measure it appears under, and will be used for self-similar processes; and S = inf{s : ξ s > }, and S = inf{s : ξ s < } will likewise be used for processes obtained by Lamperti transform. 6

7 Proof of Theorem. Our proof makes use of the pssmp (X, P ) and its Lamperti transform (ξ,p ), both defined in section 2. Let < x, y <. Then σ [,] U (x, dy) = E x ½ [X t dy] dt = E x τ ½ [X t dy] dt, using nothing more than the definition of (X, P ). We now use the Lamperti representation to relate this to (ξ,p ). This process is killed at the rate q given in (4), and so it may be represented as an unkilled Lévy process (ξ,p) which is sent to some cemetery state at the independent exponental time e q. We now make the following calculation, in which S and T are the time-changes used in the Lamperti transform as described in section 2, U (x, dy) = E log(x) T (S ) S = y α E log(x) S = y α Ê log(/x) [ ] ½ e ξ S(t) dy e αξ S(t) ds(t) [ ] ½ e ξs dy ½ [e q > s] ds ½ [ξ s log(/dy)] e qs ds, where Ê refers to the dual Lévy process. Examining the proof of Theorem VI.2 in Bertoin [2] reveals that, for any a >, S Ê a ½ [ξ s ] e qs ds = q [, ) ˆP (ξ eq dw) ˆP (ξ eq dz)½ [a w z ], [,a] where for each t, ξ t = sup{ξ s : s t} and ξ t = inf{ξ s : s t}. Then, provided that the measures ˆP (ξ eq ) and ˆP (ξ eq ) possess respective densities g S and g I (as we will shortly see they do), it follows that for a >, S Ê a ½ [ξ s dv] e qs ds = dv q a (av) dz g I (z)g S (v a z). We may apply this result to our potential measure U in order to find its density, giving u (x, y) = q yα x t g y I (log t )g S (log(tx/y)) dt. (5) x It remains to determine the densities g S and g I of the measures ˆP (ξ eq ) and ˆP (ξ eq ). These can be related to functionals of X by the Lamperti transform: ˆP (ξ eq ) = P (ξ eq ) = P ( log X τ ) ˆP (ξ eq ) = P (ξ eq ) = P ( log X τ ), (6) 7

8 where X and X are defined in the obvious manner. The laws of the rightmost random variables in (6) are available explicitly, as we now show. For the law of X τ, we transform it into an overshoot problem and make use of Example 7 in Doney and Kyprianou [6], as follows. We omit the calculation of the integral, which uses [7, 8.38.]. P (X τ dy) = ˆP ( X τ dy) = K y dv du (v y) αρ (v u) (α) ( y) αˆρ dy = sin(παˆρ) y αˆρ ( y) αˆρ dy. (7) π For the law of X τ, consider the following calculation. P (X τ y) = P (τ y < τ ) = P /y(τ < τ ). This final quantity depends on the solution of the two-sided exit problem for the stable process; it is computed in Rogozin [3], where it is denoted f (/y, ). Note that [3] contains a typographical error: in Lemma 3 of that work and the discussion after it, the roles of q (which is ρ in our notation) and q should be swapped. In the corrected form, we have P (X τ Γ(α) y) = Γ(α) = /y y u αˆρ ( u) αρ du t α (t ) αρ dt, (8) which gives us the density. Now we substitute (7) and (8) into (5): u (x, y) = xαˆρ y αρ x t α (t ) αρ y x ( t y x) αˆρ dt. The expression in the statement follows by a short manipulation of this integral. Proof of Theorem 5. According to Baurdoux [, Theorem 4.], since X is regular upwards, we have the following formula for r (, y): r (, y) = lim z u (z, y) P z (τ < τ ). We have found u above, and as we already mentioned, we have from Rogozin [3] that P x (τ < τ ) = Γ(α) x t αˆρ ( t) αρ dt. 8

9 We may then make the following calculation, using l Hôpital s rule on the second line since the integrals converge, r (, y) = Γ(α) yα lim z z(y) yz z s αˆρ (s ) αρ ds t αˆρ ( t) αρ dt = z αˆρ ( y) αˆρ (y z) αˆρ Γ(α) yα lim [ ] z z = Γ(α) yαρ ( y) αˆρ. z αˆρ z z [ z(y) ] yz Finally, the full potential density r (x, y) follows simply by substituting in the following formula, from the same theorem in []: r (x, y) = u (x, y) P x (τ < τ )r (, y). References [] E. J. Baurdoux. Some excursion calculations for reflected Lévy processes. ALEA Lat. Am. J. Probab. Math. Stat., 6:49 62, 29. ISSN [2] J. Bertoin. Lévy processes, volume 2 of Cambridge Tracts in Mathematics. Cambridge University Press, Cambridge, 996. ISBN [3] R. M. Blumenthal and R. K. Getoor. Markov processes and potential theory. Pure and Applied Mathematics, Vol. 29. Academic Press, New York, 968. [4] R. M. Blumenthal, R. K. Getoor, and D. B. Ray. On the distribution of first hits for the symmetric stable processes. Trans. Amer. Math. Soc., 99:54 554, 96. ISSN [5] M. E. Caballero and L. Chaumont. Conditioned stable Lévy processes and the Lamperti representation. J. Appl. Probab., 43(4): , 26. ISSN doi:.239/jap/ [6] R. A. Doney and A. E. Kyprianou. Overshoots and undershoots of Lévy processes. The Annals of Applied Probability, 6():9 6, 26. doi:.24/ [7] I. S. Gradshteyn and I. M. Ryzhik. Table of integrals, series, and products. Elsevier/Academic Press, Amsterdam, seventh edition, 27. ISBN Translated from the Russian. Translation edited and with a preface by Alan Jeffrey and Daniel Zwillinger. [8] A. Kuznetsov and J. C. Pardo. Fluctuations of stable processes and exponential functionals of hypergeometric Lévy processes. arxiv:2.87v [math.pr], 2. 9

10 [9] A. E. Kyprianou. Introductory lectures on fluctuations of Lévy processes with applications. Universitext. Springer-Verlag, Berlin, 26. ISBN [] A. E. Kyprianou. First passage of reflected strictly stable processes. ALEA Lat. Am. J. Probab. Math. Stat., 2:9 23, 26. ISSN [] A. E. Kyprianou, J. C. Pardo, and V. Rivero. Exact and asymptotic n-tuple laws at first and last passage. Ann. Appl. Probab., 2(2): , 2. ISSN doi:.24/9-aap626. [2] J. Lamperti. Semi-stable Markov processes. I. Z. Wahrscheinlichkeitstheorie und Verw. Gebiete, 22:25 225, 972. [3] B. A. Rogozin. The distribution of the first hit for stable and asymptotically stable walks on an interval. Theory of Probability and its Applications, 7(2): , 972. ISSN 4585X. doi:.37/735. [4] J. Vuolle-Apiala. Itô excursion theory for self-similar Markov processes. Ann. Probab., 22(2): , 994. ISSN