A Meyer-Itô formula for stable processes via fractional calculus

The infinitesimal generator of a one-dimensional strictly α\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha $$\end{document}-stable process can be represented as a weighted sum of (right and left) Riemann-Liouville fractional derivatives of order α\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha $$\end{document} and one obtains the fractional Laplacian in the case of symmetric stable processes. Using this relationship, we compute the inverse of the infinitesimal generator on Lizorkin space, from which we can recover the potential if α∈(0,1)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha \in (0,1)$$\end{document} and the recurrent potential if α∈(1,2)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha \in (1,2)$$\end{document}. The inverse of the infinitesimal generator is expressed in terms of a linear combination of (right and left) Riemann-Liouville fractional integrals of order α\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha $$\end{document}. One can then state a class of functions that give semimartingales when applied to strictly stable processes and state a Meyer-Itô theorem with a non-zero (occupational) local time term, providing a generalization of the Tanaka formula given by Tsukada [1]. This result is used to find a Doob-Meyer (or semimartingale) decomposition for |Xt-x|γ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$|X_t - x|^{\gamma }$$\end{document} with X a recurrent strictly stable process of index α\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha $$\end{document} and γ∈(α-1,α)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\gamma \in (\alpha -1,\alpha )$$\end{document}, generalizing the work of Engelbert and Kurenok [2] to the asymmetric case.


Introduction and statement of the results
One might argue that the connection between fractional calculus and stable processes can be further strengthened even though links between the fractional Laplacian and symmetric stable processes are often alluded to (cf.[MS19,LPGea20] with Remark 1.5).The oldest references that relate fractional calculus and stable random variables are the seminal work of Feller [Fel52], which uses fractional calculus to compute a series for stable densities and the articles of Gorenflo and Mainardi (cf.[GM98,MPG07]), which identify a correspondence between stable characteristic functions and the Fourier transform of fractional derivatives.More recent references are the book of Meerschaert and Sikorskii [MS19] and the article of Kolokoltsov [Kol15], where the infinitesimal generator of a stable process and various transformations are written in terms of different types of fractional derivatives.
One objective of this work is to present a natural application of fractional calculus to (one-dimensional and asymmetric) stable processes by inverting their infinitesimal generator.The inversion is valid in the so called Lizorkin space.Multiple consequences include a generalization of the celebrated Tanaka formula for Brownian motion into the stable setting as first obtained by Tsukada [Tsu19].This will follow from constructing a function which the generator transforms into the δ distribution.More generally, one can define a class of functions whose image under the generator is a signed measure.One obtains semimartingales when applying functions of this class to stable processes; the semimartingale decomposition gives us a version of Meyer-Itô formula for discontinuous semimartingales (cf.Protter [Pro04, IV.7]) which features a non-zero local time term.This allows for a concrete semimartingale decompositions for power functions applied to stable processes which were recently obtained for symmetric stable processes in Engelbert and Kurenok [EK19].
This work is based on several results, from both probability theory and fractional calculus, so let us first state the basic elements we will need.
Definition 1.1 (Strictly stable process).A Lévy process pX t q tě0 is called a strictly stable process with index of stability α P p0, 2qz t1u if, for any c ą 0: X ct d " c 1{α X t We will only consider strictly stable processes in this paper, excluding the cases when the index is 1 or 2, corresponding to the symmetric Cauchy process and Brownian motion, which have been studied with different techniques.Stable processes belong to the class of Lévy processes (c.f.[Ber96]) and their properties are deferred to the next section.
According to [Sat99](Chapter 3) there exist some constants c ´, c `ě 0, not both zero, such that the Lévy measure ν of X, which describes the jumps of X and is given by νpAq " Ep#tt P r0, 1s : X t ´Xt´P Auq, satisfies: νpdhq " `c´1 l thă0u `c`1 l thą0u ˘dh |h| α`1 .Stable processes can then be constructed using a Poisson random measure N with intensity ds νpdhq if α P p0, 1q or the compensated Poisson random measure Ñ when α P p1, 2q by means of the Lévy-Itô decomposition: .
In fact, note that in the recurrent case when α P p1, 2q, X t is integrable for any t and X is a martingale (whenever X 0 is deterministic).In both cases, we will write that X " S α pc ´, c `q when we refer to a strictly stable process with such parameters.The infinitesimal generator L of X can be defined as the derivative at zero of the semigroup on an adequate class of functions.Indeed, recall that if φ : R Ñ R belongs to the Schwartz space SpRq of rapidly decreasing functions, we have The behavior and further properties of the process X differ substantially whether α P p0, 1q or α P p1, 2q, so they will be studied separately (as above, there are many differences in these cases, such as the transient/recurrent dichotomy, the polar/non-polar character of zero, or the bounded vs unbounded variation of the sample paths).Nevertheless, in both cases we get a representation of their infinitesimal generator and its inverse in terms of fractional operators.The fractional operators we will use are the Riemann-Liouville's, these definitions and further properties can be consulted in [SKM93, Ch. 2].
(1) In fact, both the fractional integral and derivative are defined under L p assumptions depending on α, but we restrict to Schwartz space so that their Fourier transforms are well defined.
(2) If α ą 0, we will use the following notation for the fractional integrals and derivatives: (3) If α P N, then the left fractional operators I α ´and D α ´, correspond to the iterated integral and classical differential operators of order α.Fractional operators can be regarded as "nice" interpolations between their corresponding integer neighbors.(4) On an adequate domain (the so called Lizorkin space, to be introduced), they satisfy the group property with respect to composition: for α, β P R, Hence, both the left and right fractional operators of order α can be inverted by fractional operators of order ´α.(5) Again on Lizorkin space, we have that W β φ Ñ W α φ as β Ñ α, as follows from the expressions of the Fourier transforms in Proposition 2.6.
With some algebraic manipulations, the infinitesimal generator of a strictly α-stable process can be written as a linear combination of left and right Riemann-Liouville fractional derivatives of order α.For a detailed proof see for example the article of Kolokoltsov [Kol15] (Section 2), or the book of Meerschaert and Sikorskii [MS19] (Section 2.2).
Remark 1.5.This representation is consistent with the case α " 2 and c ´" c `, which corresponds to the Brownian motion, and its infinitesimal generator is the Laplacian ∆.In the case α P p0, 2qz t1u and c ´" c `, corresponding to a symmetric strictly α-stable process, the infinitesimal generator is given by the fractional Laplacian ´p´∆q α{2 .
The semigroup property of fractional operators is no longer enough to invert the infinitesimal generator, since we are lacking an expresion for the composition of left and right fractional operators.The result of this computation is stated in the forthcoming Proposition 1.7.
The main problem working in the fractional calculus framework is the domain of definition of these operators; Schwartz space is not invariant under fractional operators (cf.[SKM93], section 8.2).Since we are seeking for the inverse of the infinitesimal generator, it is useful to have a space which remains invariant under the action of the Riemann-Liouville fractional operators.This kind of space has been thoroughly studied by Lizorkin [Liz69,Liz71], Samko, Kilbas and Marichev [SKM93] and Rubin [Rub96,Rub15].
Then, the space of functions whose Fourier transforms are in Ψ is called the Lizorkin space and is defined by Φ " tφ P SpRq |F rφs P Ψ u .
In the Lizorkin space, compositions of fractional operators are well defined and therefore fractional integrals are the inverses of fractional derivatives.In general, to invert the generator, we need to see how crossed compositions are computed.The following result is stated, without proof, for fractional integrals in the article of Feller [Fel52].
Working in the Lizorkin space and using the last result we can compute the inverse of the infinitesimal generator of a stable process: Theorem 1 (Inverse of the Infinitesimal Generator).Let α P p0, 2qz t1u, c ´, c `ě 0, not both zero.Consider X " S α pc ´, c `q with infinitesimal generator L.Then, L is invertible in Φ and for every φ P Φ (2) and the constants M i as defined in Proposition 1.4.
Note that Lizorkin space is known to be dense in the space of continuous functions vanishing at infinity and in L p (cf. [Liz71] and [Sam95]), so that the above inversion formula is quite general.As an application of equation (2), the following known results can be recovered: (1) For the case α P p0, 1q, the Lévy process X is transient.Therefore, its potential corresponds to the inverse of the negative of the infinitesimal generator, p´Lq ´1.The above theorem recovers the expression of Sato [Sat72,Example 5.4].(2) For the case α P p1, 2q, the Lévy process X is recurrent and its classical potential is infinite.Nevertheless, Port [Por67] defined the recurrent potential for stable processes (by an appropriate compensated kernel) and computed it explicitly.As Sato [Sat72] notes, for a wide class of Lévy processes, the limit lim λÑ0 pλ ´Lq ´1 corresponds to a potential (classical or recurrent).On Lizorkin space, where we can explicitly compute an inversion thanks to the above theorem, Port's computation can be recovered.(3) A heuristic explanation of the function involved in the Tanaka formula for strictly stable processes given by Tsukada in [Tsu19] can be given as follows.Note that the Itô formula for Lévy processes (see Proposition 2.4) tells us that for any Schwartz function f , writing g " Lf , we have where M f is a martingale whose explicit expression is only needed later.Formally, if g equals the Dirac δ distribution, the last summand equals the time that X spends at x on r0, ts, which is one guiding principle behind the construction of the local time of X at x. Hence, if LF " δ (which will be given a sense in Section 2 and a proof in Lemma 3.2), then the local time should equal F px `Xq ´M F .Our formula for L ´1 allows us to guess a solution to LF " δ as a linear combination κ ´px ´qα´1 `κ`p x `qα´1 , which is exactly the formula of Tsukada.That κ ´‰ κ `in general is a manifestation of the asymmetry in the jumps of X.To state the Tanaka and Meyer-Itô formulae for stable processes, we need more preliminaries concerning definition of local time and an important class of admissible functions.We now consider α P p1, 2q for states to be recurrent and local time to be non trivial.
Definition 1.8 (Occupational local time).Consider a family of random variables with two indices tL a t pXq : a P R, t ě 0u.We will call it an occupational local time of a process X if the occupation time formula is satisfied for any positive Borel measurable function f : R Ñ r0, 8q: The fact that this local time exists for recurrent stable processes, as well as being jointly continuous in time and space, was established by Boylan [Boy64] and Barlow [Bar88].See the textbook account in [Ber96, Ch.V].
The following definition corresponds to the function that appears in the Tanaka formula given by Tsukada in [Tsu19], but we will write it in our notation.Definition 1.9.For every fixed α P p1, 2q, c ´, c `ě 0, not both zero, we define the function F " F α,c´,c`b y: ´`2c `c´c os απs It is intentional that κ ´accompanies x `because of Lemma 3.2.As we have remarked and will prove after Proposition 3.2, F is a weak solution to the Poisson equation LF " δ.If we consider adequate measures µ for which the convolution F ˚µ is well defined, we could regard f pxq " pF ˚µqpxq as a solution to Lf " µ.
Definition 1.10 (The class C α,c´,c`) .For every fixed α P p1, 2q and c ´, c `ě 0, not both zero, the class C α is defined as The integrability condition on µ implies that the convolution is well defined and pointwise finite.The Meyer-Itô formula will feature functions f " F ˚µ P C α,c´,c`w here µ is finite and of compact support.Later, in Theorem 3, we will consider convolutions where µ is a non compactly supported measure.This class of functions is quite large.Indeed, it contains the absolute value function and functions of the type |x| γ for γ P pα ´1, αq (cf.Lemma 3.2).Therefore, differences of convex functions are contained in C α,c´,c`.The case γ " α ´1 is special in that we can only prove its membership to C α,c´,c`i n the symmetric case.
Recall that the Meyer-Itô theorem for semimartingales, for example from [Pro04](Theorem 70), gives a semimartingale decomposition for |X| which contains a semimartingale local time term.However, the latter is zero for a strictly stable process.For functions in the class C α,c´,c`w e prove the following occupational Meyer-Itô theorem, with a non-zero local time term.
Theorem 2 (Occupational Meyer-Itô formula).Let α P p1, 2q, c ´, c `ě 0, not both zero, and consider a strictly stable process X " S α pc ´, c `q.Let f " F ˚µ P C α,c´,c`a nd furthermore assume that µ is finite and compactly supported.Then, (4) is a martingale and L a t pXq is the occupational local time at a up to time t of X.
The novel part in this result is the representation of the semimartingale in terms of an occupational local time.
‚ In the limiting case α " 2, we have F ˘pxq " x ˘and the corresponding class C can be identified with that of differences of convex functions (cf.[KS91, Thm.6.22]).‚ For recurrent symmetric stable process, that is α P p1, 2q and c ´" c `" c ą 0, we have F pxq " κ α,c |x| α´1 for some constant κ α,c (cf.[SY07, Corollary 1]).‚ The Tanaka formula of Tsukada [Tsu19], corresponds to the case where f " F " F ˚δ. ‚ The compact support hypothesis of µ is sufficient to ensure the integrability of all the terms in (4).Since strictly stable processes have finite κ-moments for κ P p´1, αq, for non compactly supported measures µ, we would at least need to verify (or assume) the integrability of f pX t q in L 1 pPq.
In general, we cannot handle the case when µ is not compactly supported, due to the integrability restrictions of strictly stable processes.Nevertheless, in the following particular case, we obtain a generalization of the works of Salminen and Yor in [SY07] and Engelbert and Kurenok [EK19] from the symmetric to the general case.Formally, the result would follow from applying Theorem 2 to the infinite measure µpdyq " |y| γ´α rk ´1 l yą0 `k`1 l yă0 s dy.Recall the definition of the constants M ˘in Proposition 1.4.
Theorem 3 (Power decomposition).Let α P p1, 2q and c ´, c `ě 0 not both zero, and consider a strictly stable process X " S α pc ´, c `q. Then for all x P R and γ P pα ´1, αq we have the decomposition where k ˘:" k ˘pα, γ, c ´, c `q are given by k ´"

Γpγ `1q
Γpγ ´α `1q Note that the last integral in (5) could be written in terms of the local time as The main result of Engelbert and Kurenok [EK19] is that this decomposition corresponds to a submartingale, thus providing the Doob-Meyer decomposition for |X t ´x| γ , when X is a symmetric stable process.However, if asymmetry in the jumps of the stable process is allowed, this decomposition will not be in general a submartingale.By direct inspection, the last term of the decomposition will correspond to an increasing process if and only if k ˘ě 0.
The constants k ˘have been found and used by Fournier [Fou13] by other means and in a different context.Fournier proved pathwise uniqueness for SDEs driven by an asymmetric strictly stable process and, in order to use the Gronwall inequality, he defined a constant βpa, cq P pα ´1, 1q, where a " cospπαq and c " c ´{c `, assuming 0 ă c ´ă c `.Then, he proved that k `" 0 for γ " βpa, cq.We will prove in Lemma 3.7 that, in fact, both k ˘are non negative for all γ ě βpa, cq and otherwise one of them is negative.
Corollary 1.12.Let a " cospπαq and c " pc ´^c `q{pc ´_ c `q. Then the power decomposition in Theorem 3 for the process |X t ´x| γ is a submartingale if γ P rβpa, cq, αq; whereas, for γ P pα ´1, βpa, cqq it is a semimartingale, whose finite variation part is not monotone.
The organization of the paper is as follows.In section 2 we state known and preliminary results regarding strictly stable processes and fractional calculus that we need for the main results.Section 3 contains proofs of the main results, examining the crossed composition Proposition 1.7, the Inversion Theorem 1, the Meyer-Itô Theorem 2 and finally the Power Decomposition Theorem 3.

Preliminaries: stable processes and fractional operators
Fractional calculus has been studied almost since the invention of calculus.One of the most famous applications is the solution to the tautochrone problem by Abel (cf.[PMT17]).Even though many mathematicians have contributed to the formalization of the field; it was Marcel Riesz who systematized several results in terms of non-local operator theory (cf.[Rie40]).The book of Samko, Kilbas and Marichev [SKM93] will be our main reference for the theory of fractional calculus in what follows.As been pointed out in the introduction, the connection between fractional calculus and stable processes will appear naturally by means of their infinitesimal generator.
In this section we state the preliminaries, regarding stable processes and fractional operators, we will need in order to prove the results outlined in the previous section.
Following Applebaum [App09] (Theorem 1.2.14 and 2.4.16),we state the Lévy-Khintchine formula and the Lévy-Itô decomposition for the special case of strictly stable processes.
Corollary 2.1 (Lévy-Khintchine formula).Let X " S α pc ´, c `q with c ´, c `ě 0, not both zero.Then its characteristic exponent, φpuq :" 1 t log E `eiuXt ˘with u P R, can be written as Moreover, it can be proved (cf.Applebaum [App09] Theorem 1.2.21) that in the case α P p0, 2qz t1u the characteristic exponent of a stable process X of index α is equal to: ¯¯ı .
Here we have another parametrization of a stable process in terms of the skewness and scale parameters, denoted by X " S α pβ, σq.We can recover the pc ´, c `q parametrization solving: Remark 2.2.The characteristic exponent of a strictly stable process and the Fourier transform of the fractional operators are intrinsically related as we will see in Remark 2.7.
The following result concerns the finiteness of moments for stable processes.For the first part, the proof can be consulted in [Tsu19] and the second one in [Ber96].
Proposition 2.3.Let α P p0, 2qz t1u, c ´, c `ě 0, not both zero, and consider a strictly stable process X " S α pc ´, c `q.Then, the following bounds are satisfied: (1) For all t ą 0, x P R and 0 ă γ ă 1, where Spα, γq is a constant which depends on α and γ and is independent of x.
(2) For all t ą 0 and 0 ď γ ă α, The following proposition is a corollary of Theorem 4.4.7 in Applebaum [App09], it is a version of Itô's formula (termed predictable in [SY07]) for stable processes.
Proposition 2.4 (Itô's formula).Let X " S α pc ´, c `q with c ´, c `ě 0, not both zero, and α P p0, 1q and f P C 2 1`,b .Then for any t ě 0, with probability 1 we have rf pX s `hq ´f pX s qs νpdhqds, and in the case α P p1, 2q we have As pointed out by Engelbert and Kurenok [EK19] in their Remark 1.1, it is a common mistake to state the Itô formula in terms of the infinitesimal generator L, since functions in C 2 1`,b are not in its domain.So, when we have a function f P C 2 1`,b , we define: f pX s `hq ´f pX s q ´hf 1 pX s q ‰ νpdhqds.
, then it coincides with the infinitesimal generator, so that L can be considered as an extension of the infinitesimal generator in the class C 2 1`,b .In the next proposition we rewrite the definition of fractional derivative depending on the index α, this representation is called the generator form.The fact that they are equivalent can be found on the book of Meerschaert and Sikorskii [MS19] and the article of Kolokoltsov [Kol15].
Proposition 2.5 (Generator form).Let f P SpRq and α P p0, 2qzt1u.Then the generator form of the left and right fractional derivatives are as follow: From this generator form follows that the infinitesimal generator of strictly stable processes can be seen as a weighted sum of fractional derivatives given in Proposition 1.4.Now we focus on the properties of the fractional operators that will lead us to the proof of the Inversion Theorem 1.
The main reason to use the Lizorkin space Φ, defined in the introduction 1.6, is that the Fourier transform of fractional operators applied to functions in Φ behaves well.First, recall that for f P SpRq the Fourier transform of f is defined by: Proposition 2.6 (Fourier transform of fractional operators).Let f P Φ and α ě 0, then the Fourier transforms of the Riemann-Liouville fractional operators of index α satisfy the following identities: The proof of this proposition can be found in the book of Samko, Kilbas and Marichev [SKM93] Lemma 8.1.Note that one of the main features of the Lizorkin space Φ, is that the Fourier transforms are well behaved near zero in such a way that the product p˘iuq ´αF rf spuq is well defined.
Remark 2.7.If we take the principal branch of logarithm, we have " |u| α ´cos ´απ 2 ¯˘i sgnpuq sin ´απ 2 ¯¯, for all u, α P R.These are precisely the characteristic functions of the one sided stable processes, see equation ( 6) with σ " 1 and β " ˘1.
Now we are ready to prove Proposition 1.7, regarding the crossed composition of Riemann-Liouville operators.First, note that from the definition of the Riemann-Liouville operators and their Fourier transforms it is easy to verify that composition of operators of the same side, left or right, commute and satisfy the semigroup property.However, the composition of crossed operators, left with right or vice versa, is not as direct as in the previous case.
Proof of Proposition 1.7.Since the Fourier transform characterizes a function φ P Φ, we will prove that the Fourier transform of both sides of the statement coincide.First, using the Fourier transform of fractional operators in Proposition 2.6 we have: Then, for the LHS of equation (1) we have: These trigonometric relations are proved in lemma A.1, finishing the proof.
Finally, we will be interested in some distributions acting on the Lizorkin space of test functions.For the definition of the action of Riemann-Liouville operators on distributions we refer to [SKM93](Section 8.1) and Rubin [Rub96](section 3).
Definition 2.8.Let f P Φ 1 and α P R. The distributions W α ´f and W α `f are defined by duality: for any φ P Φ and pg, φq denotes the evaluation of the distribution g on the function φ.
Note that δ belongs to Φ 1 (and, indeed, any Schwartz distribution) since Φ is contained in S. The infinitesimal generator L associated to X " S α pc ´, c `q corresponds to a linear combination of fractional derivatives as in Proposition 1.4.Then, if we denote its dual operator by L, this corresponds to the infinitesimal generator associated to X " S α pc `, c ´q.So that, if we take f P Φ 1 , then for any φ P Φ we have (8) pLf, φq " ´f, Lφ ¯.
In the next section we are going to use these results to prove the those outlined in the Introduction.

Proof of the main results
The objective of this section is to prove the Inversion Theorem 1, the occupational Meyer-Itô formula stated as Theorem 2 and the Doob-Meyer/semimartingale decomposition of Theorem 3 together with Corollary 1.12.
The representation in Proposition 1.4 is crucial to work out the Inversion Theorem 1.While it is well understood that the left (right) fractional derivatives is the inverse of the left (right) fractional integral, at least in Lizorkin space, the action of the crossed compositions, as far as we know, has not been reported yet.For instance, in the symmetric case, where c ´" c `, the infinitesimal generator corresponds to the fractional Laplacian ´p´∆q α{2 and its inverse operator is known in the literature as the Riesz potential (cf.[SKM93]); however, since in this case the left and right fractional derivatives merge into the fractional Laplacian, a crossed composition does not appear.
Before the proof of the Inversion Theorem 1, we will prove one more lemma regarding the composition of fractional derivatives and integrals.Since we are taking functions in the Lizorkin space, these compositions are well defined.
Lemma 3.1 (Fractional compositions).Let φ P Φ and α ą 0 with α R N, then the compositions of fractional derivatives and integrals of order α satisfy: Proof.The first two equations as well as the fact that all the compositions of fractional operators commute follow from Proposition 2.6.For the last equation, we use Proposition 1.7 with λ " α and µ Ñ ´α, to get the result.This last limit can be taken since the composition groups µ Ñ W µ ˘are continuous on Lizorkin space.Using equation (1) twice to obtain both cross compositions, we have: Where we used that W 0 is the identity operator as in Definition 1.2.
We are ready to prove the invertibility of the infinitesimal generator L in Φ and its expression as a weighted sum of fractional integrals of order α P p0, 2qz t1u.

Proof. (Inversion Theorem 1) Define the operator G as
Gφ pxq " K ´I α ´φ pxq `K`I α `φ pxq , we will prove that G pLφq " L pGφq " φ, so that L is invertible and L ´1 " G.By our definition of G, we have: Substituting the values of K ´and K `, and defining M " M 2 ´`M 2 ``2M ´M`c ospαπq to temporarily ease notation, and using Lemma 3.1, we get ´`M 2 ``2M ´M`c ospαπq M φ pxq " φ pxq .
We conclude that pG ˝Lqφ " φ and analogous computations prove that pL ˝Gqφ " φ.
Considering the following generalized functions in the dual space Φ 1 , using Definition 2.8 we will prove an important relationship between the Dirac δ distribution and the power functions, which are strongly related with the strictly stable processes.Moreover, the following lemma could be regarded as the key result to obtain Tanaka type formulae.Lemma 3.2.If λ ą 0, then, the (generalized) functions f λ `pxq :" x λ 1l txą0u and f λ ´pxq :" |x| λ 1l txă0u belong to Φ 1 and f λ `pxq " Γpλ `1qI λ`1 ´δ pxq , f λ ´pxq " Γpλ `1qI λ`1 `δ pxq .Therefore, Proof.The computation of I α ˘δ is found in [SKM93, Ch. 2 §8,p.153].It follows from Definition 2.8, the Inversion Theorem 1 and the previous Lemma 3.2 that: ´pxq.
If we substitute the values of K ´and K `in terms of α, c ´and c `we will get that L ´1pδq " F α,c´,c`i n the sense of Φ 1 distributions.
Thus, Theorem 1 provides an insight to the function that satisfies the Tanaka formula.The class of convolutions f " F α,c´,c`˚µ in Definition 1.10, is defined in such a way that the distribution induced by the measure µ coincides with Lf , in the sense of Φ 1 distributions.As a consequence, µ can be considered as the extension of Lf from the Lizorkin space to the class C c of continuous functions with compact support.A precise version of this is contained in the following lemma.It is here that the completely balanced averages of Lizorkin play a fundamental rôle: they constitute a way to approximate δ and other distributions from within Lizorkin space.
Finally, if µ is a finite measure with compact support, then φ Þ Ñ pLf, φq extends by continuity to φ Þ Ñ pµ, φq from Φ to C c with the topology of uniform convergence.
Proof.Let φ P Φ.Since α ´1 P p0, 1q, then x Þ Ñ x α´1 is subadditive on r0, 8q.Hence, Lizorkin, in [Liz71] (cf.after Definition 1.6), gives an approximation of δ in Φ 1 by means of a collection of functions κ β P Φ with the following property.If φ P C c , then φ β :" κ β ˚φ Ñ φ uniformly on compact sets; note that φ β P Φ.Indeed, Lizorkin writes κ β " κ 1 β ´κ2 β where κ 1 β is a centered Gaussian density of variance 2β 2 .Hence κ 1 β ˚φ Ñ φ uniformly if φ P C c .On the other hand, the proof of Theorem 1 [Liz71, Ch.II §4] tells us that κ 2 β ˚φ Ñ 0 uniformly on compact sets since φ is integrable.Hence, Φ is dense in C c If µ is finite and of compact support then it also has a finite moment of order α´1 and so, by the previous paragraph, LpF ˚µq " µ in Φ 1 .The bounded linear functional φ Þ Ñ pµ, φq on C c coincides with φ Þ Ñ pLf, φq on L, so that, by denseness, the latter extends uniquely by continuity to C c .The result in Lemma 3.3 with the Brownian motion case, where Lf pxq " 1 2 ∆f pxq, C 2,c,c corresponds to the class of differences of convex functions, whose second derivative are signed measures.
The following results are inspired by the work of Tsukada [Tsu19], which will be generalized by a well-known procedure to construct approximations of a function, smoothing it with mollifiers (cf.[KS91], Theorem 6.22), allowing us to use Itô formula (2.4).
A positive real function ρ P C 8 c , with support in r´1, 1s and integral equal to one, is said to be a mollifier.Then, if we consider a sequence of functions given by ρ n pxq " nρpnxq for all n P N, this sequence converges weakly to the Dirac δ distribution in the sense of Schwartz distributions, that is ˇˇˇż 8 ´8 ρ n pxqφpxqdx ´φp0q ˇˇˇÝ Ñ 0, as n Ñ 8, for all φ P S.
Let C 8 1`,b be the family of continuous functions with bounded derivatives of any order greater than or equal to one.We are going to use some bounds for the function F α,c´,c`a s well as of its increments, for a proof of the following results we refer to [Tsu19](cf.equation (3.9) of the proof of Theorem 3.1 and the proof of Lemma 3.1 in that reference).For fixed α, c ´and c `, to ease the notation, we are going to write F instead of F α,c´,c`w hen there is no confusion with the parameters.Also, recall the constants κ ˘in the definition of F and write κ " κ ´_ κ `, so that 0 ď F pxq ď κ|x| α´1 .Lemma 3.4.Let α P p1, 2q, c ´, c `ě 0, not both zero, and consider a strictly stable process X " S α pc ´, c `q.Consider the function F α,c´,c`i n equation (3), then the following results are satisfied: (1) Let pρ n q ně1 as above, then F n :" F α,c´,c`˚ρ n P C 8 1`,b for all n P N and F n Ñ F , uniformly on compact sets as n Ñ 8.
(2) Let |h| ď 1, a P R, s ą 0 and ǫ 0 ď pα ´1q ^p2 ´αq, then we have: where c 1 " 20κ 2 and the constant Sp¨, ¨q as in Proposition 2.3, and the same bound holds if we replace F by F n .Moreover, this bound satisfies: (3) Let |h| ą 1, a P R and s ą 0, then we have: where c 2 " 4κ and the same bound holds if we replace F by F n .Moreover, this bound satisfies: The following result is a corollary of Lemma 3.4 and it will be useful in several steps of the Meyer-Itô theorem's proof.
Corollary 3.5.Under the assumptions of Lemma 3.4, let f P C α,c´,c`, such that f " F ˚µ P C α,c´,c`w ith µ a finite Radon measure and consider f n " f ˚ρn for n P N. Then we have: and the same bounds are satisfied if we replace f with f n .These bounds are an elements of L 1 pp0, tq ˆA, Bpp0, tq ˆAq, Leb bνqq, with A " r´1, 1szt0u and A " r´1, 1s c respectively.
These results follow from the Lemma 3.4 and an application of a Jensen-like inequality for finite measures.
Proof of the Occupational Meyer-Itô Formula (Theorem 2).Without loss of generality, we asume that µ is actually a positive measure, which was assumed to be finite with compact support and, therefore, with moments of order α and 2pα ´1q.Then, we have the representation: ´8 F α,c´,c`p x ´aq µ pdaq .Consider the sequences F n " F ˚ρn and f n " f ˚ρn " F ˚ρn ˚µ as the infinitely differentiable approximations of F and f by the sequence tρ n u ně0 , with n P N, and we have that f n Ñ f uniformly on compact sets ([EG15] Theorem 4.1: Properties of mollifiers).
Since f n P C 8 1`,b Ă C 2 , using Itô's formula (Proposition 2.4) we have: (10) where the last two terms are rf n pX s´`h q ´fn pX s´q s Ñ pds, dhq and Moreover, since the behavior of M n t is different depending on the size of the jumps, we will consider M n t " M 1,n t `M 2,n t , where rf n pX s´`h q ´fn pX s´q s Ñ pds, dhq.
In a similar fashion, we define M t " M 1 t `M 2 t , by replacing f n with f .The proof consists in establishing the following steps: Step 1: f pX t q and f n pX t q are in L 1 pPq and f n pX t q Ñ f pX t q in L 1 .
Step 2: M 1 and M 1,n are square integrable martingales and M 1,n Step 3: M 2 and M 2,n are integrable martingales and M 2,n t Ñ M 2 t in L 1 .
Step 4: V n t Ñ ş L a t µpdaq in L 1 .Let's begin with Step 1. First, we provide a bound for f pxq and f n pxq in terms of x and which does not depend on n.Using that α ´1 P p0, 1q, we have that x Þ Ñ x α´1 is subadditive on r0, 8q, so that 0 ď f n pxq " which is finite for any x P R by the assumptions on µ and does not depend on n.
By similar arguments we have that For the squared difference, using a Jensen-like inequality for finite measures, we have, Then, similar arguments give and these bounds are independent of n and belong to L 1 pPq since 0 ă 2α ´2 ă α and µ is a finite measure with a moment of order 2α ´2.We can conclude that f n pX t q and f pX t q are elements of L 2 pPq.Moreover, by dominated convergence, we get |f n pX t q ´f pX t q| 2 ı " 0, so that f n pX t q Ñ f pX t q in L 2 pPq, which implies Step 1's assertions.Let's move to Step 2. In this case we are considering the jumps smaller than one, i.e. h ď 1.To prove that M 1,n is a square integrable martingale, according to Ikeda and Watanabe ([IW89] section II.3), we need to show that: Since the integrand is positive and pX , BpX qq-measurable with X " pΩˆr´1, 1szt0uˆr0, tsq, by the Fubini theorem (cf.[Kal02] Theorem 1.27), it suffices to prove the finiteness in any order of integration.Then, using the bound in Corollary 3.5 for |h| ď 1 and Lemma 3.4, we have: The result for m 1 t follows from Corollary 3.5 in a similar fashion.Hence, M 1 is also a square integrable martingale.
In order to prove the convergence of M 1,n t Ñ M 1 t in L 2 pPq, first note that according to Corollary 3.5 we have Thus, pM 1 n q ně1 is dominated in L 1 pp0, tq ˆr´1, 1szt0u, Bpp0, tq ˆr´1, 1szt0uq, Leb bνqq.We know that pM 1,n t ´M 1 t q is a square integrable martingale for any n P N, then using Itô's isometry ([App09] p. 223) and dominated convergence theorem for the sequence pM 1 n q ně1 we have: |f n pX s´`h q ´fn pX s´q ´pf pX s´`h q ´f pX s´q q| 2 ı νpdhqds |f n pX s´`h q ´fn pX s´q ´pf pX s´`h q ´f pX s´q q| 2 ı νpdhqds " 0.
The convergence to zero of the last equation is a consequence of equation (13 Step 2. For Step 3, we are considering the jumps greater than one, i.e. h ą 1.To prove that M 2,n t is a martingale, following Ikeda and Watanabe ([IW89] section II.3) we must show: Since the integrand is positive and pX , BpX qq-measurable with X " pΩ ˆr´1, 1s c ˆr0, tsq, by the Fubini theorem it suffices to prove the finiteness in any order of integration.Then, using the bound in Corollary 3.5 for |h| ą 1 and Lemma 3.4, we have: The result for m 2 t follows by the same bounds in Corollary 3.5, so that M 2 t is also a martingale.As in the previous step, to prove the convergence of M 2,n Ñ M 2 in L 1 pPq, first note that according to Corollary 3.5 we have M 2 n :" E r|f n pX s´`h q ´fn pX s´q ´pf pX s´`h q `f pX s´q q|s ď E r|f n pX s´`h q ´fn pX s´q |s `E r|f pX s´`h q ´f pX s´q |s ď 2µpRqc 2 |h| α´1 .
Thus, pM 2 n q ně1 is dominated in L 1 pp0, tq ˆr´1, 1s c , Bpp0, tq ˆr´1, 1s c q, Leb bνqq.We know that pM 2,n t ´M 2 t q is a stochastic integral with respect to a Poisson random measure for any n P N, then using Campbell's theorem ([Kin93] section 3.2) and dominated convergence theorem for the sequence pM 2 n q ně1 we have: E r|f n pX s´`h q ´fn pX s´q ´pf pX s´`h q ´f pX s´q q|s νpdhqds.
The convergence to zero of the last equation is a consequence of equation (13) in Step 1.So that M 2,n t Ñ M 2 t in L 1 pPq, ending with Step 3. By Step 2 and Step 3 we conclude that M n and M in equation ( 10) are martingales and Finally, for Step 4, we have from equation (10) that: f n pX t q ´fn pX 0 q ´M n t L 1 pPq Ñ f pX t q ´f pX 0 q ´Mt , as n Ñ 8, so that the limit lim nÑ8 V n t pX t q P L 1 pPq.We just need to verify that this limit coincides with the one stated in the theorem.
We know that f n " F ˚pρ n ˚µq P C 8 1`,b XC α,c´,c`i s positive and measurable and that ρ n ˚µ is a finite measure with compact support.Then Lf n is well defined, positive and measurable as well.So, by the occupation formula we have: Since L a t pωq P C c for almost all ω P Ω, Lemma 3.3 tells us that ´8 L a t pµ ˚ρn qpdaq, and since ρ n Ñ δ weakly as n Ñ 8, then pµ ˚ρn q Ñ µ weakly as n Ñ 8 as well.Hence, ˇˇˇż 8 ´8 L a t pµ ˚ρn qpdaq ´ż 8 ´8 L a t µpdaq ˇˇˇÑ 0, as n Ñ 8.
Steps 1-4 finish the proof of Theorem 2.
For a first application, we have the Tanaka formula for asymmetric strictly stable processes.
Corollary 3.6 (Tanaka formula).Let α P p1, 2q, c ´, c `ě 0, not both zero, and consider a strictly stable process X " S α pc ´, c `q.Then, the Tanaka formula is satisfied: (14) F α,c´,c`p X t ´aq " F α,c´,c`p X 0 ´aq `M a t pXq `La t pXq, where L a t pXq is the occupational local time at a up to time t of X and M a t pXq is a square integrable martingale given by M a t pXq " rF α,c´,c`p X s´´a `hq ´F α,c´,c`p X s´´a qs Ñ pds, dhq.
Proof.Consider the unitary measure concentrated in a, that is δ a pEq " 1 if a P E and zero otherwise, with f pxq " pF α,c´,c`˚δ a q pxq " F α,c´,c`p x´aq, using the occupational Meyer-Itô theorem we have: F α,c´,c`p X t ´aq " F α,c´,c`p X 0 ´aq `ż t 0 ż R0 rF α,c´,c`p X s´´a `hq ´F α,c´,c`p X s´´a qs Ñ pds, dhq `ż 8 ´8 L x t pXq δ a pdxq , " F α,c´,c`p X 0 ´aq `M a t `La t pXq .
We turn our attention to the power decomposition of Theorem 3. Our first step will be to explicitly compute the infinitesimal generator of the power functions in Lemma 3.2.=Let α P p1, 2q, c ´, c `ě 0 not both zero and α ´1 ă γ ă α.From Lemma 3.2 we know that f γ belong to Φ 1 and can be identified with the following fractional integrals: Before we prove Theorem 3, we need to undestand the constants k ˘pα, γ, c ´, c `q that are used there.They play an important role in the bounded variation part of the power decomposition (5), because in order to be an increasing process, both need to be positive.The following lemma states the critical exponent γ from which both k ˘pα, γ, c ´, c `q are positive.Recall the definition of c in Corollary 1.12.
We will prove in a similar way the change of signs of k `pα, γ, c ´, c `q.Note that, as in the previous case, we just need to analyze the change of signs of the function: h `pγq :" c sin p´απq ´sin pγπq ´c sin ppγ ´αqπq .
Finally, when c `ă c ´, just note that since k `pα, γ, c ´, c `q " k ´pα, γ, c `, c ´q we can use the same proof.
We are ready to prove the power decomposition theorem.These results are a generalization of the works of Salminen and Yor [SY07] and of Engelbert and Kurenok [EK19].The proof of the decomposition uses the Tanaka formula for asymmetric stable processes (14) and relies on the representation of the infinitesimal generator of a power function given in Lemma 3.2.Note that in [SY07] it was easy to find the measure which could recover the power decomposition in the symmetric case and for the generalization we made direct use of fractional calculus to find the relevant measure needed for the asymmetric case.
Taking the Tanaka formula ( 14) at the level a and integrating both sides by µ x pdaq (the measure µ translated by x) we have: Note that the representation of f as a member of the Class C α,c´,c`i s precisely F ˚µ.We will now use a version of Fubini's theorem for compensated Poisson random measures and apply it to the small jumps of M a pXq above.See [MR15, Lemma A.1.2].We need to verify some integrability assumptions to apply it, which are (15) and ( 16) below.Applying the Fubini theorem, we get Using the occupational formula for the local time, the last integral is equivalent to This finishes the proof modulo showing that the first integral is a martingale and the applicability of Fubini's theorem.The proof of the martingale character will follow the ideas of [EK19, Section 3].Incidentally, the same argument will justify the application of Fubini's theorem above.We can identify two cases depending on the size of the jump: In order to prove that M 1,γ is a square integrable martingale, according to Ikeda and Watanabe ([IW89] section II.Take c " c ´_ c `, then the intensity measure ν c pdhq " c|h| ´α´1 dh is greater than the intensity measure νpdhq, corresponding to X t , and if we consider the change of variable h " pX s´´x qu we have: is finite for all t ě 0.
It remains to check that the second integral is finite.Consider the auxiliary function gpuq " |1 `u| γ , and note that for any u P p´1, 1q we have that gpuq " p1 `uq γ , which is differentiable.By the mean value theorem we can choose u ˚P p´1, 0q and u ˚P p0, 1q such that: This corresponds to: p1 `uq γ ´1 " # γp1 `u˚q γ´1 u ´1 ă u ă 0, γp1 `u˚qγ´1 u 0 ă u ă 1.
We get the following bound for any u P p´1, 1q: |p1 `uq γ ´1| ď γc 1 pγq|u|, where c 1 pγq " maxpp1 `u˚q γ´1 , p1 `u˚qγ´1 q.Then, we have that So that m 1,γ t for any t ě 0 and M 1,γ is a square integrable martingale.Now, to prove that M 2,γ is a martingale, according to Ikeda  ||1 `u| γ ´1| c |u| α`1 du ă 8, since γ ´α P p´1, 0q and this moment of X t is finite for any t ě 0, the expectation is finite.To see the last integral is finite, just note that ||1 `u| γ ´1| behaves like |u| γ as |u| Ñ 8.Then, we have that m 2,γ t is finite for any t ě 0 and we can conclude that M 2,γ is a martingale.This allow us to conclude that M γ " M 1,γ `M 2,γ is a martingale.
Finally, we state when this power decomposition is a submartingale or just a semimartingale.
Proof.(of Corollary 1.12) From the Lemma 3.7 and Theorem 3 the last integral is a non decreasing process if and only if γ P rβpa, cq, αq, by Lemma 3.7, so that we get a Doob-Meyer decomposition for |X t ´x| γ .In the other case, γ P pα ´1, βpa, cqq, this results in a semimartingale instead of a submartingale.

Appendix A. Trigonometric results
The following trigonometric result is used in the proof of the composition of crossed fractional operators.Proof.Using the trigonometric identities for the sum of angles we start from the LHS: The following lemma is used to analyze the constant βpα, cq in Theorem 1.12.
Lemma A.2.The functions h ˘of Lemma 3.7 have minimum period 2.
Proof.Let f ˘pxq " h ˘px{2πq, so that we now wish to prove that the minimum period of f ȋs 2π.First, note that f ˘is a solution to f 2 `f " 0. Second, all solutions to the above ODE are given by a cos `b sin.Finally, we assert that the minimum period of the above linear combination is 2π as long as a and b are not both zero.Let us assume that a ‰ 0. If f˘p xq " f ˘px `pq for some p, by equating initial conditions at zero, we obtain a " a cos p `b sin p and b " ´a sin p `b cos p.
By substituting the value for b obtained in the second equation in the first and cancelling a, since it is non-zero, we get 1 ´cos 2 p " sin 2 ppq " p1 ´cos pq 2 .
Expanding the square, we get cos p " cos 2 p from which p " 2kπ.The case when b ‰ 0 is handled similarly.