Return to Article Details On the enhanced structures and rough paths integrals for time-changed paths

On the Enhanced Structures and Rough Path Integrals for Time-Changed PathsThanks:  Department of Mathematics, Tunis El Manar University, Faculty of Sciences of Tunis, Laboratory of Mathematical Analysis and Applications, LR11ES11, El Manar I, Tunisia, e-mail: bedhiafi.mounir@yahoo.fr

Mounir Bedhiafi
Date: October 28, 2025; accepted: April 08, 2026; published online: August 31, 2026.
Abstract.

In this paper, we study the effects of a discontinuous time change Φ on a continuous path x. While rough integration in this setting poses significant challenges, we construct an enhancement for the time-changed path xΦ, providing a structured framework for discontinuous-time rough paths. We then apply our approach to Lévy processes, demonstrating its effectiveness in specific stochastic models, including the variance gamma process and the normal inverse Gaussian process.

A key focus of our work is the study of rough path integrals under time changes. We establish conditions under which the Lyons integral of the continuously time-changed rough path XΦ coincides with the Φ-time-changed rough path integral driven by the original rough path X. As a central result, we extend Terry Lyons’ Itô formula to a time-changed Itô formula, offering a novel perspective on rough path calculus under time deformations.

These findings enhance the theoretical understanding of rough integration under time changes, with potential applications in stochastic analysis and mathematical finance.

Key words and phrases: 
Rough paths, time change, Itô formula, Lévy process, variance gamma process, normal inverse Gaussian process.
2005 Mathematics Subject Classification
28C05, 60G07, 60G51

1. Introduction

The theory of rough paths, pioneered by Terry Lyons, provides a powerful framework for analyzing differential equations driven by highly irregular signals, including paths with low regularity and stochastic processes with jumps. A fundamental aspect of rough path theory is the construction of pathwise integrals, enabling the robust study of differential equations beyond classical Itô and Stratonovich calculus.

A natural question arises when considering time-changed rough paths, where the driving signal undergoes a time transformation. Time changes are fundamental in stochastic analysis, appearing in models where volatility varies dynamically or where time is distorted due to external influences. This motivates the study of rough integration in the context of time-changed paths.

In the case of discontinuous time changes, additional challenges emerge. Unlike the continuous-time setting, rough paths enhanced over a discontinuous time change require additional structure. We address this issue by constructing a suitable enhancement of the time-changed path xΦ, ensuring a meaningful interpretation within the rough path framework.

We then apply our results to Lévy processes, which serve as natural candidates for studying rough integration under discontinuous time changes. Specifically, we provide explicit enhancements for key examples such as the variance gamma process and the normal inverse Gaussian process, demonstrating the applicability of our theoretical framework to real-world stochastic models.

A central focus of this work is the relationship between different formulations of rough path integration under time changes. We investigate whether the Lyons integral of a time-changed rough path XΦ coincides with the Φ-time-changed integral of X. Under suitable conditions, assuming the time change Φ is continuous, we establish this equivalence and derive a novel time-changed Itô formula, extending Lyons’ results to time-deformed rough paths.

This paper is structured as follows: In Section 2, we review the necessary background on rough paths and time changes.

In Section 3, we construct enhancements for discontinuous time-changed paths, apply our results to Lévy processes, and provide concrete examples.
In Section 4, we study rough path integrals under time changes and establish conditions under which different integration formulations coincide.

2. Preliminaries

2.1. Rough Paths

In this section, with reference to [5, 7, 8, 9, 10], we gather notations, definitions, and recapitulate some fundamental concepts.

Let (E,d) be a metric space. A path x:[0,T]E is said to be:

  • (i)

    α-Hölder continuous if

    |x|α-Höl,[0,T]:=sup0s<tTd(xs,xt)|ts|α<.
  • (ii)

    of finite p-variation for some p>0 if

    |x|p-var,[0,T]:=suptiD([0,T])(id(xti,xti+1)p)1p<,

    where D([0,T]) denotes the set of all finite dissections of [0,T].

Let V be a Euclidean space. For each n, we denote by Tn(V) the truncated tensor algebra:

Tn(V)=VV2Vn.

A multiplicative functional in Tn(V) is a map 𝕏=(1,𝕏1,,𝕏n) defined on the simplex

Δ={(s,t)[0,T]20stT,T>0}

with values in Tn(V), satisfying the Chen identity:

𝕏s,u=𝕏s,t𝕏t,u,for all stu.

A multiplicative functional 𝕏=(1,𝕏1,,𝕏n) in Tn(V) has finite total p-variation if

supDl|𝕏tl1tli|pi<+,i=1,,n,

where supD runs over the set D of all finite dissections of [0,T].

Multiplicative functionals with finite total p-variation in T[p](V) are called rough paths (of roughness p), and their set is denoted by Ωp(V). The space Ωp(V) is endowed with the p-variation metric:

dp(𝕏,𝕐)=maxsupD0i[p](l|𝕏tl1tli𝕐tl1tli|pi)ip.

An enhanced path 𝕏 of degree n associated with a path x in some Banach space is a multiplicative functional 𝕏s,t=(1,𝕏s,t1,,𝕏s,tn) from Δ into Tn(V) such that 𝕏s,t1=xs,t. The path x is called the trace of X, and 𝕏 is called the enhancement of x.

Definition .

A control function, is a continuous non-negative function w on ΔT which is super-additive in the sense that

w(s,t)+w(t,u)w(s,u)stu,(s,t)[0,T]2,

and for which w(t,t)=0 for all t[0,T].

Lemma ([9, Lemma 1.10, p. 7]).

Let ω be a control. Let x:JE be a continuous path. Assume that, for some p1 and for all (s,t)ΔT, one has

|xsxt|pω(s,t).

Then, for all (s,t)ΔT,

xp-var;[s,t]ω(s,t)1/p.

When the conclusion of this lemma holds, we say that the p-variation of x is controlled by ω.

A rough path 𝕏Ωp(V) is called a smooth rough path if txt𝕏0,t1 is a continuous path with finite variation and 𝕏s,ti is the i-th iterated integral of the path xt over the interval [s,t], for i=1,,[p].

That is,

𝕏s,ti=s<t1<<ti<tdxt1dxti(s,t)Δ.

A multiplicative functional 𝕏Ωp(V) is called geometric rough path if there is a sequence 𝕏(n) of smooth rough path in Ωp(V) such that

dp(𝕏(n),𝕏)0asn+.

We denote by GΩp(V) the set of all geometric rough paths of roughness p in T[p](V).

A function 𝕏:ΔT[p](V) is called an almost rough path (of roughness p) if it is of finite p-variation, 𝕏s,t0=1, and, for some control ω and some constant θ>1,

|(𝕏s,t𝕏t,u)i𝕏s,ui|ω(s,u)θ,

for all (s,t),(t,u)Δ and i=1,,[p].

In [10, p. 41], it was shown that if 𝕏 a θ-almost p-rough path controlled by ω, then there exists a unique p-rough path 𝕏^:ΔT[p](V) and a constant K which depends only on p,θ and ω(0,T) such that

sup0s<tT𝕏^s,ti𝕏s,tiw(s,t)θ<K

and the p-variation of 𝕏^ is controlled by Kω.

2.2. Rough Path Integration

In the following, we assume that:

  • α:VL(V,W) is a function that maps elements of V linearly to W-valued one-forms on V.

  • α possesses all continuous and bounded derivatives dkα up to order [p].

  • We denote αi=di1α for i=1,2,,[p].

2.2.1. Integration of Rough Paths of Degree Two

Definition .

Let 𝕏Ωp(V). The integral of the one-form α against the rough path 𝕏, denoted by

α(𝕏)𝑑𝕏,

is the unique p-rough path in T2(W) associated with the almost rough path 𝕐, where:

𝕐s,t1=α1(𝕏s)(𝕏s,t1)+α2(𝕏s)(𝕏s,t2),

and

𝕐s,t2=α1(𝕏s)α1(𝕏s)(𝕏s,t2).

With 𝕏s𝕏0,s.

2.2.2. Itô’s Formula

Lemma ([10, Lemma 5.4.1, p. 135]).

Let 𝕏Ωp(V), and define

𝕐s,t=(1,𝕏s,t1,𝕏s,t2+As,t).

Then 𝕐Ωp(V) if and only if A is additive and possesses finite p2-variation. We denote 𝕐 by 𝕏A.

Theorem 1.

Let 𝕏Ωp(V), and let A be a continuous path in V2 which possesses finite p2 variation. Then

stα(𝕏A)d𝕏A1 = stα(𝕏)d𝕏1+stα2(𝕏r)(dAr)stα(𝕏A)d𝕏A2
= stα(𝕏)d𝕏2+stα1(𝕏r)α1(𝕏r)(dAr)
+st[srα2(𝕏r)(dAr)]dr+sts,rα2(𝕏r)(dAr)
+st(suα2(𝕏r)dAr)α2(𝕏u)dAu.

Here the integrals involving A are Young’s integrals, and s,t=stα(𝕏)d𝕏1.

2.2.3. Stochastic Integration

Theorem 2.

Let (Wt)t0 be a d-dimensional Brownian motion, and let

𝕏s,t=(1,𝕏s,t1,𝕏s,t2)

be the geometric rough path associated with W. Then:

(α(𝕏)𝑑𝕏)s,tstα(𝕏)𝑑𝕏=(1,stα(𝕏)d𝕏1,stα(𝕏)d𝕏2),

where:

Zs,t=stα(𝕏)d𝕏1=stα(Wr)dWr,

and

stα(𝕏)d𝕏2=s<t1<t2<tdZt1dZt2.

Here, the integrals are Stratonovich integrals.

2.2.4. Integration Against Geometric Rough Paths

Definition .

Let 𝕏GΩp(V). The integral α(𝕏)𝑑𝕏 is the unique rough path in Ωp(W) associated with the quasi α-differential almost rough path 𝕐, where:

𝕐s,ti=l=(l1,,li),1lj[p]αl1αli(𝕏s)(πΠlπ1𝕏s,t|l|).

2.3. Time-Changed Young Integration and the Rough Path Approach

It is well known (see, for example, [2]) that if x and y are, respectively, α-Hölder and β-Hölder continuous functions on [0,T] such that α+β>1, and if Φ is an increasing γ-Hölder continuous function from [0,T] into itself satisfying (α+β)γ>1, then for any 0abT, we have the following identity for the Young integral:

ab(xΦ)(u)d(yΦ)(u)=Φ(a)Φ(b)x(u)𝑑y(u).

This result provides a fundamental relationship between integration and time changes in the Young integration framework.

Terry Lyons extended these ideas by developing a general theory of rough paths, which gives meaning to integrals of the form

Zt=Z0+0tα(xs)dxs,

where x is a path of finite p-variation in a Banach space V. To define such an integral rigorously, one needs to lift x to a path X of finite p-variation in the free nilpotent group of V. This lifting procedure enables the application of rough path integration techniques, which generalize Young’s integration by incorporating iterated integrals.

In this paper, we consider a path x with finite p-variation and a time change Φ, and we investigate the integral:

ab(xΦ)(u)d(xΦ)(u),for 0abT,

where the integral is interpreted in the rough path sense using Lyons’ theory.

3. Enhanced time-changed paths

In [8], T. Lyons developed a theory that provides a rigorous framework for solving the integral equation

(1) Yt=Y0+0tf(xs)dxs,

where x is a continuous path of finite p-variation in a Banach space V with p>1. This requires lifting x to an enhanced path 𝕏 of finite p-variation in the free nilpotent group of V, allowing the integral to be defined using rough path techniques.

Building on this foundation, Peter Friz [4] extended Lyons’ approach to handle cases where the driving process x is a discontinuous path of finite p-variation in a Banach space V. In this work, we further explore this setting by considering the integral equation (1) under a time change Φ.

Given a time change Φ, we define the transformed process xΦ(t) and denote it by xtΦ. Our goal is to provide a rigorous interpretation of the equation:

(2) Yt=Y0+0tf(xsΦ)dxsΦ.

To achieve this, we employ Lyons’ rough path integration framework, which necessitates constructing an enhancement of the time-changed process xΦ. The development that follows is dedicated to establishing this enhancement and its properties.

Definition .

A time change Φ is a family Φs,s0, of stopping times such that the map sΦs is almost surely increasing and right-continuous.

3.1. Finite discontinuities

Let V be a Banach space and let x:[0,T]V be a continuous path. Assume that x has a geometric enhancement 𝕏, that is,

𝕏s,t=(1,𝕏s,t1,𝕏s,t2),

where:

𝕏s,t1=xtxs,

and

𝕏s,t2,i,j=limn+s<ln<t(xlnixsi)(xl+1njxlnj).
Theorem 3.

Let s,t[0,T] and let Φ be a time change with a finite number of discontinuities a1,a2,,ar in [s,t]. Assume that akln for all k=1,,r. Then,

𝕏s,tΦ,2=𝕏Φ(s),Φ(t)2k=1r[𝕏Φ(ak),Φ(ak+)2+(xΦ(ak)xΦ(s))(xΦ(ak+)xΦ(ak))].
Proof.

By definition,

𝕏s,tΦ,2,i,j=st(xΦ,i(u)xΦ,i(s))dxΦ,j(u).

Taking the limit as n+ over a partition,

𝕏s,tΦ,2,i,j =limn+[l/n<a1(xΦ(l/n)ixΦ(s)i)(xΦ((l+1)/n)jxΦ(l/n)j)
+k=1r1ak<l/n<ak+1(xΦ(l/n)ixΦ(s)i)(xΦ((l+1)/n)jxΦ(l/n)j)
+l/n>ar(xΦ(l/n)ixΦ(s)i)(xΦ((l+1)/n)jxΦ(l/n)j)].

Using integration properties,

𝕏s,tΦ,2 =Φ(s)Φ(a1)(xuxΦ(s))dxu+k=1r1Φ(ak+)Φ(ak+1)(xuxΦ(s))dxu
+Φ(ar+)Φ(t)(xuxΦ(s))dxu.

Rewriting the sums using enhancement properties,

𝕏s,tΦ,2 =𝕏Φ(s),Φ(a1)2+k=1r1𝕏Φ(ak+),Φ(ak+1)2+𝕏Φ(ar+),Φ(t)2
+k=1r1(xΦ(ak+)xΦ(s))(xΦ(ak+1)xΦ(ak+))
+(xΦ(ar+)xΦ(s))(xΦ(t)xΦ(ar+)).

Using the repeated Chen identity,

𝕏Φ(s),Φ(t)2 =𝕏Φ(s),Φ(a1)2+k=1r𝕏Φ(ak),Φ(ak+)2+k=1r1𝕏Φ(ak+),Φ(ak+1)2+𝕏Φ(ar+),Φ(t)2
+k=1r(xΦ(ak)xΦ(s))(xΦ(ak+)xΦ(ak)).

Rearranging the terms completes the proof. ∎

3.2. Countable discontinuities

Theorem 4.

Let x:[0,T]V be a continuous path admitting a geometric rough path lift 𝕏=(1,𝕏1,𝕏2)GΩp(V) with 2p<3. Let Φ be an increasing càdlàg time change with countably many discontinuities (ak)k in (s,t). Assume that there exists a control ω such that

|𝕏u,v1|ω(u,v)1p,|𝕏u,v2|ω(u,v)2p,

and

s<ak<tω(Φ(ak),Φ(ak+))2p<.

Then the second level of the enhanced time-changed path xΦ is given by

𝕏s,tΦ,2=𝕏Φ(s),Φ(t)2s<ak<t[𝕏Φ(ak),Φ(ak+)2+(xΦ(ak)xΦ(s))(xΦ(ak+)xΦ(ak))],

and 𝕏Φ defines a geometric rough path of roughness p.

Proof.

The proof proceeds by approximation. We construct a sequence (Φn)n1 of increasing càdlàg time changes, each having only finitely many discontinuities, such that ΦnΦ pointwise. For each Φn, we apply Theorem 3 (finite jump case) to obtain an explicit expression for the second level 𝕏Φn,2. We then show that the sequence of rough paths (𝕏Φn)n1 is Cauchy with respect to the p-variation metric and hence converges to a limiting rough path, denoted by 𝕏Φ. Finally, we identify this limit and verify that it satisfies the announced formula, and that it retains the geometric rough path structure.

Step 1: Approximation of Φ. Let (ak) be the discontinuity points of Φ in (s,t), ordered increasingly. Define

ΔΦ(ak):=Φ(ak+)Φ(ak).

Since Φ is increasing,

k=1ΔΦ(ak)Φ(t)Φ(s)<.

For each n, define

Φn(u):=Φ(u)k>n:akuΔΦ(ak).

Then Φn is increasing, right-continuous, has finitely many jumps (a1,,an), and

Φn(u)Φ(u),u.

Step 2: Finite jump formula. By the finite discontinuity case, for each n,

𝕏s,tΦn,2=𝕏Φn(s),Φn(t)2k=1n[𝕏Φ(ak),Φ(ak+)2+(xΦ(ak)xΦn(s))(xΦ(ak+)xΦ(ak))].

Step 3: Convergence of the second level. We show that (𝕏Φn,2) is Cauchy uniformly. Let m>n. Then

𝕏s,tΦm,2𝕏s,tΦn,2=Am,n+Bm,n+Cm,n,

where:

Am,n:= 𝕏Φm(s),Φm(t)2𝕏Φn(s),Φn(t)2,
Bm,n:= k=n+1m[𝕏Φ(ak),Φ(ak+)2+(xΦ(ak)xΦm(s))(xΦ(ak+)xΦ(ak))],
Cm,n:= (xΦm(s)xΦn(s))k=1n(xΦ(ak+)xΦ(ak)).

Estimate of Am,n. By continuity of (u,v)𝕏u,v2,

Am,n0.

Estimate of Bm,n. Using the control,

|𝕏Φ(ak),Φ(ak+)2| ω(Φ(ak),Φ(ak+))2/p,
|xΦ(ak+)xΦ(ak)| ω(Φ(ak),Φ(ak+))1/p.

Hence

|Bm,n|Ck=n+1mω(Φ(ak),Φ(ak+))2/p.

By assumption, the series converges, so Bm,n0.

Estimate of Cm,n. We have

|Cm,n||xΦm(s)xΦn(s)|k=1n|xΦ(ak+)xΦ(ak)|.

Using the control,

k=1|xΦ(ak+)xΦ(ak)|k=1ω(Φ(ak),Φ(ak+))1/p.

Since ω1/p1+ω2/p for small ω, the series converges. Thus the sum is bounded, and since Φn(s)Φ(s),

Cm,n0.

Step 4: Convergence in dp. The same estimates applied to increments over partitions show that

dp(𝕏Φn,𝕏Φm)0.

Since (Ωp(V),dp) is complete, there exists a limit 𝕏Φ.

Step 5: Identification of the limit. Passing to the limit,

𝕏s,tΦ,2=𝕏Φ(s),Φ(t)2s<ak<t[𝕏Φ(ak),Φ(ak+)2+(xΦ(ak)xΦ(s))(xΦ(ak+)xΦ(ak))].

Step 6: Chen identity. Each 𝕏Φn satisfies Chen’s identity. Passing to the limit preserves it, hence 𝕏Φ is a geometric rough path. ∎

The summability condition ensures that the correction terms define a path of finite p/2-variation.

Theorem 5.

Let x:[0,T]V be a continuous path admitting a geometric rough path lift X=(1,X1,X2)GΩp(V) with 2p<3. Assume that there exists a control ω such that

|Xs,t1|ω(s,t)1p,|Xs,t2|ω(s,t)2p.

Let Φ be an increasing càdlàg time change with jump times (ak) in (s,t). Assume that

s<ak<tω(Φ(ak),Φ(ak+))2p<.

Define

As,t:= s<ak<tXΦ(ak),Φ(ak+)2,
Γs,t:= s<ak<t(xΦ(ak)xΦ(s))(xΦ(ak+)xΦ(ak)),

and

X~s,t2:=XΦ(s),Φ(t)2As,tΓs,t.

Then the path

Xs,tΦ:=(1,XΦ(s),Φ(t)1,X~s,t2)

defines a geometric rough path of finite p-variation.

Proof.

Step 1: Additivity. For s<u<t, the sums defining As,t and Γs,t split over disjoint sets of indices, hence

As,t=As,u+Au,t,Γs,t=Γs,u+Γu,t.

Thus both A and Γ are additive.

Step 2: p/2-variation of A. Let D={ti} be a partition of [0,T]. Then

i|Ati,ti+1|p2=i|ti<ak<ti+1XΦ(ak),Φ(ak+)2|p2.

Using the triangle inequality and the control on X2, we obtain

i|Ati,ti+1|p2i(ti<ak<ti+1ω(Φ(ak),Φ(ak+))2p)p2.

Since p21, there exists a constant Cp>0 such that

(xk)p2Cpxkp2.

Hence

i|Ati,ti+1|p2Cpkω(Φ(ak),Φ(ak+)).

By assumption, this sum is finite, so A has finite p2-variation.

Step 3: p/2-variation of Γ. We estimate

|Γs,t|s<ak<t|xΦ(ak)xΦ(s)||xΦ(ak+)xΦ(ak)|.

Using the control,

|xuxv|ω(u,v)1p,

we obtain

|Γs,t|s<ak<tω(Φ(s),Φ(ak))1pω(Φ(ak),Φ(ak+))1p.

Since ω is a control,

ω(Φ(s),Φ(ak))ω(Φ(s),Φ(t)),

hence

|Γs,t|ω(Φ(s),Φ(t))1ps<ak<tω(Φ(ak),Φ(ak+))1p.

Using Hölder’s inequality and the summability assumption, we deduce that Γ has finite p2-variation.

Step 4: Control of X~2. We have

X~s,t2=XΦ(s),Φ(t)2As,tΓs,t.

Each term has finite p2-variation, hence so does X~2.

Step 5: Chen identity. Since X satisfies Chen’s identity and A, Γ are additive, we obtain

X~s,t2=X~s,u2+X~u,t2+XΦ(s),Φ(u)1XΦ(u),Φ(t)1.

Thus XΦ satisfies Chen’s identity.

Conclusion. The path XΦ has finite p-variation and satisfies Chen’s identity, hence it is a geometric rough path. ∎

3.3. Algebraic structure in the case p3

Theorem 6.

Let x:[0,T]V be a continuous path admitting a geometric rough path enhancement X. Let s<t and let Φ be a time change with finitely many discontinuities a1,,ar in (s,t). Then, for any k1, the enhanced time-changed rough path satisfies

(3) Xs,tΦ,k=XΦ(s),Φ(t)km=1r[XΦ(am),Φ(am+)k+j=1k1XΦ(s),Φ(am)jXΦ(am),Φ(am+)kj].
Proof.

By definition, the iterated integrals of the time-changed path satisfy

Xs,tΦ,k=stXs,uΦ,k1dxuΦ.

Approximating the integral by Riemann sums over partitions (tln) of [s,t], we write

Xs,tΦ,k=limn[ tln<a1XΦ(s),Φ(tln)k1(xΦ(tl+1n)xΦ(tln))
+m=1r1am<tln<am+1XΦ(s),Φ(tln)k1(xΦ(tl+1n)xΦ(tln))
+tln>arXΦ(s),Φ(tln)k1(xΦ(tl+1n)xΦ(tln))].

Passing to the limit yields

Xs,tΦ,k =Φ(s)Φ(a1)XΦ(s),uk1dxu+m=1r1Φ(am+)Φ(am+1)XΦ(s),uk1dxu+Φ(ar+)Φ(t)XΦ(s),uk1dxu.

Using the multiplicative (Chen) property of the rough path, we write for u[Φ(am+),Φ(am+1)]:

(4) XΦ(s),uk1=j=0k1XΦ(s),Φ(am)jXΦ(am),uk1j.

Substituting this decomposition into the integrals, we obtain

Xs,tΦ,k =XΦ(s),Φ(a1)k+m=1r1XΦ(am+),Φ(am+1)k+XΦ(ar+),Φ(t)k
+m=1r1j=1k1XjΦ(s),Φ(am)XkjΦ(am),Φ(am+1)
+j=1k1XjΦ(s),Φ(ar)XkjΦ(ar),Φ(t).

Now, applying Chen’s identity to the partition

s<a1<a1+<<ar<ar+<t,

we have

XΦ(s),Φ(t)k =XΦ(s),Φ(a1)k+m=1rXΦ(am),Φ(am+)k+m=1r1XΦ(am+),Φ(am+1)k+XΦ(ar+),Φ(t)k
+m=1rj=1k1XjΦ(s),Φ(am)XkjΦ(am),Φ(am+).

Rearranging terms yields exactly

Xs,tΦ,k=XΦ(s),Φ(t)km=1r[XΦ(am),Φ(am+)k+j=1k1XΦ(s),Φ(am)jXΦ(am),Φ(am+)kj].

Finally, since this representation is obtained by subtracting local jump contributions from a Chen decomposition, it follows that XΦ inherits Chen’s identity from X. ∎

The previous results provide an explicit construction of the enhanced time-changed path. We now show that this construction is in fact canonical.

3.4. Uniqueness of the Enhanced Time-Changed Path

The constructions in  Theorem 3Theorem 5 provide an explicit enhancement of the time-changed path xΦ:=xΦ. A natural question is whether this enhancement is uniquely determined by the original rough path 𝕏 and the time change Φ.

Continuous time changes. If Φ is continuous, then the enhancement is given by

𝕏s,tΦ=𝕏Φ(s),Φ(t).

In this case, 𝕏Φ is a geometric rough path and coincides with the canonical lift of xΦ. Uniqueness follows from the universal limit property of geometric rough paths: the lift is obtained as the limit of smooth approximations, and any other candidate with the same trace and finite p-variation must coincide with it (see [5, Proposition 7.10, p. 132]).

Discontinuous time changes. When Φ has jumps, the path xΦ is no longer continuous, and the resulting object is not a geometric rough path in the classical sense. Nevertheless, one can still define an enhancement satisfying Chen’s identity and suitable p-variation bounds. The formulas in Theorem 3 and Theorem 4 are dictated by the following structural requirements:

  1. 1)

    Chen’s identity: the enhanced path must satisfy the multiplicative (Chen) relation on all intervals.

  2. 2)

    Consistency away from jumps: on any interval containing no discontinuity of Φ, the enhancement must reduce to the canonical composition 𝕏Φ(s),Φ(t).

  3. 3)

    Finite p-variation: the second level must satisfy the p/2-variation bounds ensured by  Theorem 5.

Uniqueness argument. These conditions determine the enhancement uniquely. Indeed, consider two candidates 𝕏Φ and 𝕏~Φ satisfying the above properties, and define their difference at level two:

As,t:=𝕏s,tΦ,2𝕏~s,tΦ,2.

Then A satisfies:

  • As,t=0 whenever (s,t) contains no discontinuity of Φ,

  • an additive (Chen-type) relation: As,t=As,u+Au,t,

  • finite p/2-variation.

Hence A is supported on the jump intervals and must be a purely atomic additive functional. The p/2-variation constraint forces the total contribution of these atoms to be summable, and the Chen relation then uniquely determines their values. This yields A0.

Relation with Lévy rough paths. This uniqueness property is consistent with the construction of rough paths over Lévy processes developed by P. K. Friz and Atul Shekhar, where a canonical lift is obtained from the Lévy–Khintchine triplet under a Blumenthal–Getoor condition. In that framework, the summability of jumps ensures both existence and uniqueness of the enhancement. In our setting,  Theorem 8 shows that the condition β<1/p guarantees the required summability of the jump corrections, and hence the uniqueness of the time-changed enhancement follows as a consequence.

Conclusion. The enhancement constructed in  Theorem 3Theorem 5 is therefore the unique rough path over xΦ satisfying Chen’s identity, finite p-variation bounds, and consistency with the original rough path away from the discontinuities of Φ.

3.5. Application on Lévy Processes

Proposition ([3, Proposition  3.5, p.  75]).

If X is a Lévy process, then the quantity

J([0,t]×A)=#{s[0,t]:ΔXsA}

defines a Poisson random measure on +×({0}) with intensity measure dtdν, where

ν(A)=𝔼[J([0,1]×A)].

3.5.1. Poisson Integral

Let J be a Poisson random measure with intensity measure dtdν on +×(d{0}).

Let f:dn be a measurable function, and let A(d{0}) be such that

ν(A)<+.

Then, for all t0 and ωΩ, the Poisson integral of f is defined by:

Af(Z)J(t,𝑑Z)=ZAf(Z)J(t,{Z}).
Remark .

1) If J is associated with a Lévy process, then

Af(Z)J(t,𝑑Z)=0stf(ΔXs)1A(ΔXs).

2) If ν is the Lévy measure and 0A¯, then ν(A)<+, ensuring that Af(Z)J(t,𝑑Z) is well-defined.

Theorem 7.

Let A be a Borel set such that 0A¯. Then:

1) For all t0, the integral Af(Z)J(t,𝑑Z) follows a compensated Poisson law characterized by:

(5) 𝔼[exp(iu,Af(Z)J(t,𝑑Z))]=exp(tA(eiu,f(Z)1)ν(𝑑Z)).

2) If fL1(A,ν), then:

𝔼[Af(Z)J(t,𝑑Z)]=tAf(Z)ν(𝑑Z).
Proof.

1) Suppose f(z) is a step function of the form:

f(z)=j=1ncj1Aj(z),

where cjn and Aj are disjoint Borel sets. Then,

𝔼[exp(iu,Af(Z)J(t,𝑑Z))] =𝔼[exp(ij=1nJ(t,Aj)u,cj)]
=j=1n𝔼[exp(iJ(t,Aj)u,cj)]
=j=1nexp(t(eiu,cj1)ν(Aj))
=exp(tA(eiu,f(Z)1)ν(𝑑Z)).

The general case follows by approximating f with step functions.

2) The expectation result follows by differentiating equation (5). ∎

Proposition .

Let W be a d-dimensional Brownian motion and let Φ be an independent subordinator with Lévy measure νΦ. Define the subordinate Brownian motion

Xt=WΦt.

Let βX and βΦ denote the Blumenthal–Getoor indices of X and Φ, respectively. Then

βX=2βΦ.
Proof.

Step 1: Lévy measure of the subordinated process. It is well known that X is a Lévy process whose Lévy measure νX is given by

νX(A)=0(WsA)νΦ(𝑑s),A(d{0}).

Step 2: Characterization of the Blumenthal–Getoor index. By definition (see [11]),

βX=inf{α>0:|x|1|x|ανX(𝑑x)<}.

We compute

|x|1|x|ανX(𝑑x)=0(|x|1|x|α(Ws𝑑x))νΦ(𝑑s).

Step 3: Scaling of Brownian motion. By Brownian scaling,

Ws=dsW1,

hence

|x|1|x|α(Wsdx)=𝔼[|Ws|α𝟏{|Ws|1}]=sα/2𝔼[|W1|α𝟏{|W1|s1/2}].

As s0, the truncation becomes irrelevant and we obtain the estimate

𝔼[|Ws|α𝟏{|Ws|1}]Cαsα/2,

for some constant Cα>0.

Step 4: Reduction to the subordinator. Thus,

|x|1|x|ανX(dx)<01sα/2νΦ(ds)<.

Step 5: Identification of the index. By the definition of the Blumenthal–Getoor index of Φ,

βΦ=inf{γ>0:01sγνΦ(𝑑s)<}.

Comparing with the previous condition, we see that

|x|1|x|ανX(dx)<α2>βΦ.

Hence

βX=2βΦ,

which completes the proof. ∎

Theorem 8.

Let W be a d-dimensional Brownian motion and let Φ be an independent subordinator. Define the Lévy process Xt=WΦt. Let β be the Blumenthal–Getoor index of X, and let 2p<3. Assume that

β<2p.

Then the enhanced time-changed path 𝕏Φ, defined via the formula of Theorem 3 and Theorem 4, is a geometric rough path of finite p-variation.

Proof.

Step 1: Blumenthal–Getoor index and jump summability. For the subordinate Brownian motion Xt=WΦt, the Blumenthal–Getoor index satisfies

β=2βΦ,

where βΦ is the index of the subordinator The condition β<2/p implies

βΦ<1p.

Let α such that

β<α<2p.

By the definition of the Blumenthal–Getoor index,

sT|ΔXs|α<a.s.

Step 2: Control of the jump contributions. We now relate this to the control required in Theorem 4Theorem 5. For the Brownian rough path 𝕏 (Stratonovich lift), there exists a control ω such that

|𝕏u,v1|ω(u,v)1/p,|𝕏u,v2|ω(u,v)2/p.

Moreover, for the time change Φ,

ω(Φ(ak),Φ(ak+))|ΔΦ(ak)|.

Since Xt=WΦt, the jumps of X satisfy

|ΔXak||ΔΦ(ak)|1/2.

Thus

ω(Φ(ak),Φ(ak+))2/p|ΔXak|4/p.

Step 3: Summability condition. Since α<2/p, we have

|ΔXak|4/pC|ΔXak|αfor small jumps.

Hence

kω(Φ(ak),Φ(ak+))2/p<a.s.

Step 4: Application of Theorem 4Theorem 5. The summability condition required in Theorem 4-5 is therefore satisfied. It follows that the correction terms define a path of finite p/2-variation. Consequently, the enhanced path 𝕏Φ is a geometric rough path of finite p-variation. ∎

3.5.2. Examples of Enhanced Time-Changed Paths

In this subsection, we present some examples of enhanced time-changed paths, focusing on two specific cases of subordinators. A Lévy process Xt on d is characterized by its Lévy-Khintchine triplet (a,σ,Π), where

𝔼[eizXt]=etψ(z),zd,

with the characteristic exponent

ψ(z)=12zσz+iaz+d(eizx1izx1{|x|1})Π(dx).

3.5.3. Enhanced Variance Gamma Process

The (λ,c)-gamma process (Γt)t0 is a Lévy process with probability density function

Pt(x)=λctΓ(ct)xct1eλx1{x>0}.

It has the Lévy-Khintchine triplet (a,σ,Π) given by

a=01xΠ(dx),σ=0,Π(dx)=cx1eλxdx.

The Blumenthal-Getoor index of the gamma process is 0, since:

|x|1|x|αΠ(𝑑x) =01xαceλxx𝑑x=c01eλxx1α𝑑x.

Since eλx behaves like 1 near x=0, the integral is finite for all α>0, implying that the index of the gamma process is 0, and thus it has finite variation.

Let (Wt)t0 be a d-dimensional Brownian motion. The variance gamma process is defined as:

(Vt)t0:=(WΓt)t0.

By Theorem 4, the enhancement of the variance gamma process is:

𝕍s,t= (1,WΓtWΓs,Γs<t1<t2<ΓtdWt1dWt2
s<ak<t[Γ(ak)<t1<t2<Γ(ak+)dWt1dWt2
+(WΓ(ak)WΓ(s))(WΓ(ak+)WΓ(ak))]).

3.5.4. Enhanced Inverse Gaussian Process

Let (Bt)t0 be a standard Brownian motion. Define the first passage time:

Is=inf{t>0:Bt+bt>s}.

The process (It)t0 is called the inverse Gaussian process. It has the Lévy-Khintchine triplet (a,σ,Π) given by:

a =2sb0b12πey2/2dy,
σ =0,
Π(dx) =s12πx3eb2x/21{x>0}dx.

The Blumenthal-Getoor index of the inverse Gaussian process is 12, since:

|x|1|x|αΠ(𝑑x) =12π01xα1x3/2eb2x/2dx.

Near x=0, the term xαx3/2eb2x/2 behaves like xα3/2, and the integral is finite if and only if α>12. Thus, the inverse Gaussian process has finite variation. Let (Wt)t0 be a d-dimensional Brownian motion. The normal inverse Gaussian (NIG) process is defined as:

(Nt)t0:=(WIt)t0.

By Theorem 4, the enhancement of the normal inverse Gaussian process is:

s,t= (1,WItWIs,Is<t1<t2<ItdWt1dWt2
s<ak<t[I(ak)<t1<t2<I(ak+)dWt1dWt2
+(WI(ak)WI(s))(WI(ak+)WI(ak))]).

3.6. Worked Examples

To illustrate the construction of enhanced time-changed paths and to verify the conditions of  Theorem 3, Theorem 4, Theorem 5, and Theorem 8, we provide several explicit examples.

Example (Deterministic time change with a single jump).

Let x(t)=t on [0,2] and define the time change

Φ(t)=t+𝟏[1,2](t),

which has a single discontinuity at a1=1 with Φ(1)=1 and Φ(1+)=2. Then xΦ(t)=x(Φ(t))=Φ(t). For s=0, t=2,  Theorem 3 gives

𝕏0,2Φ,2=𝕏1,22[𝕏1,22+(x(1)x(1))(x(2)x(1))]=0.

Direct computation using the definition of the iterated integral confirms the result:

02(xΦ(u)xΦ(0))dxΦ(u)=01u𝑑u+1+2(u1)𝑑u+jump term=0.

The correction term exactly cancels the contribution from the jump, verifying the formula.

Example (Countable discontinuities – convergence verification).

Let Φ(t)=t+k=12k𝟏[12k,1](t) on [0,1]. This time change has countably many discontinuities accumulating at t=1; the jump sizes are ΔΦ(12k)=2k. Let x(t)=t. For any s<t, the sum

s<ak<t[𝕏Φ(ak),Φ(ak+)2+(xΦ(ak)xΦ(s))(xΦ(ak+)xΦ(ak))]

converges absolutely because k(ΔΦ(ak))2<. By taking partitions that refine at each discontinuity, the finite-discontinuity approximations converge uniformly, and the limit satisfies the Chen identity. Hence  Theorem 4 applies.

Example (Variance gamma process – numerical check).

Let Γt be a gamma process with shape c=1 and rate λ=1, and let Wt be a one-dimensional Brownian motion. Consider the time interval [0,1] and suppose Γ has a single jump at t=0.5 of size ΔΓ=0.2. Using the enhancement formula from  Theorem 8,

𝕍0,1= (1,WΓ1W0,Γ0Γ1dWsdWs
[Γ(0.5)Γ(0.5+)dWsdWs+(WΓ(0.5)W0)(WΓ(0.5+)WΓ(0.5))]).

The Blumenthal–Getoor index of Γ is 0. For p=2.5, we have β=0<1/p, so  Theorem 8 guarantees that 𝕍 is a finite p-variation rough path. A direct simulation confirms that the Chen identity holds for the three intervals [0,0.5], [0.5,1], and [0,1].

Example (Normal inverse Gaussian process – parameter limitation).

Let It be an inverse Gaussian process with parameters b=1, s=1. Its Blumenthal–Getoor index is β=12. Choose p=2.5; then β=0.51/2.5=0.4, so the hypothesis of  Theorem 8 is not satisfied. Indeed, the sum |ΔIt|1/p may diverge, and the enhancement would require additional regularization (e.g., using a different lift or restricting to a subclass of time changes with β<1/p). For p=4, we have β=0.50.25, again failing the condition.

Remark (Connection with Friz–Shekhar).

The above result is consistent with the general construction of rough paths associated with Lévy processes developed by Friz and Shekhar [4]. In their work, a canonical geometric rough path lift is constructed directly from the Lévy–Khintchine triplet, without relying on a time-change representation. More precisely, for a Lévy process X with Blumenthal–Getoor index β<1/p, they show that X admits a natural rough path lift with finite p-variation. In this framework, the summability of the jumps follows from the classical property:

0sT|ΔXs|α<a.s. for all α>β.

Our approach is complementary: instead of constructing the lift directly, we start from a given rough path and study its stability under discontinuous time changes. The summability condition appearing in Theorem 4 and Theorem 5 is recovered here as a consequence of the Blumenthal–Getoor condition β<1/p. In particular, Section 3.5.1 shows that the time-change procedure is compatible with the Friz–Shekhar construction, and yields the same class of p-variation rough paths.

4. Integrals Driven by Enhanced Continuous Time-Changed Paths

In this section, we develop a calculus for integrals driven by continuous time-changed rough paths. We establish conditions under which the enhancement of a time-changed rough path preserves its rough structure and show how this framework applies to rough integration.

Lemma .

Let x be a continuous finite p-variation path and let Φ be a continuous time change. Then xΦ is also a finite p-variation path.

Proof.

We have:

|xΦ|p-var,[0,T] =suptiD([0,T])(id(xΦti,xΦti+1)p)1p
suptiD([0,T])(id(xti,xti+1)p)1p
=|x|p-var,[0,T]<.

Lemma .

Let Φ:[a,b][0,T] be a continuous time change. Assume that Dn={t1n,,tknn} is a sequence of finite partitions of [a,b] such that the mesh size

|Dn|=sup1ikn|ti+1ntin|

satisfies |Dn|0 as n. Then, for each n, the set

D~n={Φ(t1n),,Φ(tknn)}

forms a partition of [Φ(a),Φ(b)], and its mesh size satisfies

|D~n|0asn.
Proof.

Since Φ is continuous and increasing, it maps a partition of [a,b] into a partition of [Φ(a),Φ(b)]. We now show that the mesh size of D~n converges to zero.

By the uniform continuity of Φ on the compact interval [a,b] (which follows from Heine’s theorem), for every ε>0, there exists δ>0 such that:

|x1x2|<δ|Φ(x1)Φ(x2)|<ε.

Since |Dn|0, there exists N such that for all nN,

sup1ikn|ti+1ntin|<δ.

Applying the uniform continuity condition, we obtain:

sup1ikn|Φ(ti+1n)Φ(tin)|<ε.

Since ε is arbitrary, we conclude that |D~n|0 as n. ∎

Theorem 9.

Let x:[0,T]V be a continuous path of finite p-variation (p>1) and let 𝕏Ωp(V) be a geometric rough path lifting x. Let ϕ:[s,t][0,T] be a continuous, increasing time change. Then the rough path integral satisfies

stα(xϕ)(u)d(xϕ)(u)=ϕ(s)ϕ(t)α(x(u))𝑑x(u),

where the left‑hand side is the Lyons integral of the time‑changed rough path 𝕏ϕ (with 𝕏s,tϕ=𝕏ϕ(s),ϕ(t)) and the right‑hand side is the Lyons integral of the original rough path 𝕏 over the transformed time interval.

Proof.

We assume 2p<3. Let 𝕏=(1,𝕏s,t1,𝕏s,t2) be the geometric rough path associated with x. In particular case of Theorem 3, since x is continuous, we have

𝕏Φ=(1,𝕏Φ(s),Φ(t)1,𝕏Φ(s),Φ(t)2).

Define =(1,s,t1,s,t2) by:

s,t1 =α1(𝕏(Φ(s)))(𝕏Φ(s),Φ(t)1)+α2(𝕏(Φ(s)))(𝕏Φ(s),Φ(t)2),
s,t2 =α1(𝕏(Φ(s)))α1(𝕏(Φ(s)))(𝕏Φ(s),Φ(t)2).

Let Dn={t0n=s,,tknn=t} be a sequence of partitions of [s,t] with mesh |Dn|0.

By Section 4, the image partitions D~n={ϕ(tin)} have mesh tending to zero.

By Lyons’ extension theorem [10, Theorem 5.2.1, p. 118], Z defines an almost rough path. The corresponding rough path ~=(1,~s,t1,~s,t2) satisfies:

~s,t1 =lim|D~n|0l=1rα1(𝕏Φ(tl1)Big)(𝕏Φ(tl1),Φ(tl)1)+α2(𝕏Φ(tl1))(𝕏Φ(tl1),Φ(tl)2)
=lim|Dn|0l=1rα1(𝕏tl1)(𝕏tl1,tl1)+α2(𝕏tl1)(𝕏tl1,tl2)
=Φ(s)Φ(t)α(𝕏)d𝕏1.

Similarly, for the second iterated integral:

~s,t2 =lim|D~n|0l=1rα1(𝕏Φ(tl1))α1(𝕏Φ(tl1))(𝕏Φ(tl1),Φ(tl)2)
=lim|Dn|0l=1rα1(𝕏tl1)α1(𝕏tl1)(𝕏tl1,tl2)
=Φ(s)Φ(t)α(𝕏)d𝕏2.

For p3, the proof extends similarly using the higher-order iterated integrals. Thus, we obtain the desired result. ∎

Corollary .

Let Wt be a d-dimensional Brownian motion, and let

𝕏s,t=(1,𝕏s,t1,𝕏s,t2)

be the geometric rough path associated with W. Suppose Φ is a continuous time change. Then:

Zs,t:=stα(𝕏Φ)d𝕏Φ1=Φ(s)Φ(t)α(Wr)dWr,

and:

stα(𝕏Φ)d𝕏Φ2=Φ(s)<t1<t2<Φ(t)dZt1dZt2.

4.1. Time-Changed Itô Formula

Lemma .

Let 𝕏Ωp(V), and let

𝕐s,t=(1,𝕏s,t1,𝕏s,t2+As,t)

where A is an additive process with finite p2-variation. Let Φ be a continuous time change. Then the time-changed process:

YΦ=(1,XΦ(s),Φ(t)1,XΦ(s),Φ(t)2+AΦ(s),Φ(t))

belongs to Ωp(V).

Proof.

Define As,tΦ=AΦ(s),Φ(t). Then, using the additivity property,

As,uΦ+Au,tΦ=AΦ(s),Φ(u)+AΦ(u),Φ(t),

which shows that AΦ remains additive. Since

𝕏s,tΦ:=(1,𝕏Φ(s),Φ(t)1,𝕏Φ(s),Φ(t)2)Ωp(V),

applying Section 2.2.2, we conclude that YΦΩp(V). ∎

Theorem 10.

Let 𝕏Ωp(V), and let A be a continuous path in V2 with finite p2-variation. Then the following formulas hold:

(6) stα(𝕐Φ)d𝕐Φ1 =Φ(s)Φ(t)α(𝕏)d𝕏1+Φ(s)Φ(t)α2(𝕏r)dAr.
stα(𝕐Φ)d𝕐Φ2 =Φ(s)Φ(t)α(𝕏)d𝕏2+Φ(s)Φ(t)α1(𝕏r)α1(𝕏r)dAr
+Φ(s)Φ(t)[Φ(s)rα2(𝕏r)dAr]dr
+Φ(s)Φ(t)Φ(s),rα2(𝕏r)dAr
(7) +Φ(s)Φ(t)[Φ(s)uα2(𝕏r)dAr]α2(𝕏u)dAu.

Here, the integrals involving A are Young integrals, and

s,t=stα(𝕏)d𝕏1.
Proof.

By definition of the integral, we have:

stα(𝕐Φ)d𝕐Φ1 =limm(D)0l[α1(𝕏Φ(tl1))𝕏Φ(tl1),Φ(tl)1
+α2(𝕏Φ(tl1))(𝕏Φ(tl1),Φ(tl)2+AΦ(tl1),Φ(tl))].

Changing variables and considering a finer partition D~, we obtain:

stα(𝕐Φ)d𝕐Φ1=
=limm(D~)0l[α1(𝕏tl1)𝕏tl1,tl1+α2(𝕏tl1Big)(𝕏tl1,tl2+Atl1,tl)]
=Φ(s)Φ(t)α(𝕏)d𝕏1+Φ(s)Φ(t)α2(𝕏r)dAr.

Similarly, for the second integral:

stα(𝕐Φ)d𝕐Φ2
= limm(D)0lα1(𝕏Φ(tl1))α1(𝕏Φ(tl1))(𝕏Φ(tl1),Φ(tl)2+AΦ(tl1),Φ(tl))
+(Φ(s),Φ(tl1)+Φ(s)Φ(tl1)α2(𝕏r)dAr)(Φ(tl1),Φ(tl)
+Φ(tl1)Φ(tl)α2(𝕏r)dAr).

Applying the partition refinement and limits, we get:

stα(𝕐Φ)d𝕐Φ2 =Φ(s)Φ(t)α(𝕏)d𝕏2+Φ(s)Φ(t)α1(𝕏r)α1(𝕏r)dAr
+Φ(s)Φ(t)Φ(s),rα2(𝕏r)dAr
+Φ(s)Φ(t)[Φ(s)uα2(𝕏r)dAr]du
+Φ(s)Φ(t)[Φ(s)uα2(𝕏r)dAr]α2(𝕏u)dAu.

Acknowledgements.

The author would like to thank the anonymous referee for a thorough and insightful report. The referee’s detailed comments and suggestions have significantly improved the presentation and the mathematical precision of the manuscript. In particular, they contributed to clarifying the arguments in Theorem 45, strengthening the treatment of countable discontinuities, and making explicit the connection with the Blumenthal–Getoor index in the finite p-variation framework. Any remaining inaccuracies are entirely the responsibility of the author.

References

  • [2] L. Coutin, Rough paths via sewing lemma, ESAIM Probab. Stat., 16 (2012), pp. 479–526.
  • [3] R. Cont and P. Tankov, Financial Modelling with Jump Processes, Chapman & Hall/CRC Financial Mathematics Series, Chapman & Hall/CRC, Boca Raton, FL, 2004.
  • [4] P.K. Friz and A. Shekhar, General rough integration, Lévy rough paths and a Lévy–Khintchine-type formula, Ann. Probab., 45 (2017) no. 4, pp. 2707–2765.
  • [5] P.K. Friz and N.B. Victoir, Multidimensional Stochastic Processes as Rough Paths, Cambridge Studies in Advanced Mathematics, vol. 120, Cambridge University Press, Cambridge, 2010.
  • [6] K. Kobayashi, Stochastic calculus for a time-changed semimartingale and the associated stochastic differential equations, J. Theoret. Probab., 24 (2011), pp. 789–820.
  • [7] A. Lejay, An introduction to rough paths, Séminaire de Probabilités XXXVII, Lecture Notes in Math., vol. 1832, Springer, Berlin, 2003, pp. 1–59.
  • [8] T.J. Lyons, Differential equations driven by rough signals, Rev. Mat. Iberoamericana, 14 (1998) no. 2, pp. 215–310.
  • [9] T.J. Lyons, M. Caruana and T. Lévy, Differential Equations Driven by Rough Paths, Lecture Notes in Mathematics, vol. 1908, Springer, Berlin, 2007.
  • [10] T.J. Lyons and Z. Qian, System Control and Rough Paths, Oxford Mathematical Monographs, Oxford University Press, Oxford, 2002.
  • [11] K. Sato, Lévy Processes and Infinitely Divisible Distributions, Cambridge Studies in Advanced Mathematics, vol. 68, Cambridge University Press, Cambridge, 1999.