arxiv: v1 [math.pr] 15 Oct 2014

Transcription

1 A Fourier approach to pathwise stochastic integration arxiv:141.46v1 [math.pr] 15 Oct 214 Massimiliano Gubinelli CEREMADE & CNRS UMR 7534 Université Paris-Dauphine and Institut Universitaire de France Nicolas Perkowski CEREMADE & CNRS UMR 7534 Université Paris-Dauphine October 16, 214 Abstract Peter Imkeller Institut für Mathematik Humboldt-Universität zu Berlin We develop a Fourier approach to rough path integration, based on the series decomposition of continuous functions in terms of Schauder functions. Our approach is rather elementary, the main ingredient being a simple commutator estimate, and it leads to recursive algorithms for the calculation of pathwise stochastic integrals, both of Itô and of Stratonovich type. We apply it to solve stochastic differential equations in a pathwise manner. Contents 1 Introduction 2 2 Preliminaries Ciesielski s isomorphism Young integration and rough paths Paradifferential calculus and Young integration Paradifferential calculus with Schauder functions Young s integral and its different components M.G. is supported by a Junior fellowship of the Institut Universitaire de France IUF) and by the ANR Project ECRU ANR-9-BLAN-114-1) N.P. is supported by the Fondation Sciences Mathématiques de Paris FSMP) and by a public grant overseen by the French National Research Agency ANR) as part of the Investissements d Avenir program reference: ANR-1-LABX-98), and acknowledges generous support from Humboldt-Universität zu Berlin, where a major part of this work was completed. 1

2 4 Paracontrolled paths and pathwise integration beyond Young Paracontrolled paths A basic commutator estimate Pathwise integration for paracontrolled paths Pathwise Itô integration 22 6 Construction of the Lévy area Hypercontractive processes Continuous martingales Pathwise stochastic differential equations 34 A Regularity for Schauder expansions with affine coefficients 36 1 Introduction The theory of rough paths [Lyo98] has recently been extended to a multiparameter setting [Hai14, GIP12]. While [Hai14] has a much wider range of applicability, both approaches allow to solve many interesting SPDEs that were well out of reach with previously existing methods; for example the continuous parabolic Anderson in dimension two[hai14, GIP12], the three-dimensional stochastic quantization equation [Hai14, CC13], the KPZ equation [Hai13, GP14], or the three-dimensional stochastic Navier Stokes equation [ZZ14b, ZZ14a]. Our methods developed in [GIP12] are based on harmonic analysis, on Littlewood-Paley decompositions of tempered distributions, and on a simple commutator lemma. This requires a non-negligible knowledge of Littlewood-Paley theory and Besov spaces, while at the same time the application to classical rough path SDEs is not quite straightforward. That is why here we develop the approach of [GIP12] in the slightly different language of Haar / Schauder functions, which allows us to communicate our basic ideas while requiring only very basic knowledge in analysis. Moreover, in the Haar Schauder formulation the application to SDEs poses no additional technical challenges. It is a classical result of Ciesielski [Cie6] that C α := C α [,1],R d ), the space of α Hölder continuous functions on [,1] with values in R d, is isomorphic to l R d ), the space of bounded sequences with values in R d. The isomorphism gives a Fourier decomposition of a Hölder-continuous function f as f = p,m H pm,df G pm, where H pm ) are the Haar functions and G pm ) are the Schauder functions. Ciesielski proved that a continuous function f is in C α [,1],R d ) if and only if the coefficients H pm,df ) p,m decay rapidly enough. Following Ciesielski s work, similar isomorphisms have been developed for many Fourier and wavelet bases, showing that the regularity of a function is encoded in the decay of its coefficients in these bases; see for example Triebel [Tri6]. But until this day, the isomorphism based on Schauder functions plays a special role in stochastic analysis, because the coefficients in the Schauder basis have the pleasant property that they are just rescaled second order increments of f. So if f is a stochastic process with known distribution, then also the distribution of its coefficients in the Schauder basis is known 2

3 explicitly. A simple application is the Lévy-Ciesielski construction of Brownian motion. An incomplete list with further applications will be given below. Another convenient property of Schauder functions is that they are piecewise linear, and thereforetheir iterated integrals G pms)dg qn s), can beeasily calculated. Thismakes them an ideal tool for our purpose of studying integrals. Indeed, given two continuous functions f and g on [,1] with values in LR d,r n ), the space of linear maps from R d to R n, and R d respectively, we can formally define t fs)dgs) := p,m H pm,df H qn,dg q,n t G pm s)dg qn s). In this paper we study, under which conditions this formal definition can be made rigorous. We start by observing that the integral introduces a bounded operator from C α C β to C β if and only if α+β > 1. Obviously, here we simply recover Young s integral [You36]. In our study of this integral, we identify different components: t fs)dgs) = Sf,g)t)+π < f,g)t)+lf,g)t), where S is the symmetric part, π < the paraproduct, and Lf,g) the Lévy area. The operators S and π < are defined for f C α and g C β for arbitrary α,β >, and it is only the Lévy area which requires α + β > 1. Considering the regularity of the three operators, we have Sf,g) C α+β, π < f,g) C β, and Lf,g) C α+β whenever the latter is defined. Therefore, in the Young regime fs)dgs) π <f,g) C α+β. We will also see that for sufficiently smooth functions F we have Ff) C α but Ff) π < DFf),f) C 2α. So both fs)dgs) and Ff) are given by a paraproduct plus a smoother remainder. This leads us to call a function f C α paracontrolled by g if there exists a function f g C β such that f π < f g,g) C α+β. Our aim is then to construct the Lévy area Lf,g) for α < 1/2 and f paracontrolled by g. If β > 1/3, then the term Lf π < f g,g),g) is well defined, and it suffices to make sense of the term Lπ < f g,g),g). This is achieved with the following commutator estimate: Lπ <f g,g),g) f g s)dlg,g)s) f g β g β g β. 3β Therefore, the integral fs)dgs) can be constructed for all f that are paracontrolled by g, provided that Lg,g) can be constructed. In other words, we have found an alternative formulation of Lyons [Lyo98] rough path integral, at least for Hölder continuous functions of Hölder exponent larger than 1/3. Since we approximate f and g by functions of bounded variation, our integral is of Stratonovich type, that is it satisfies the usual integration by parts rule. We also consider a non-anticipating Itô type integral, that can essentially be reduced to the Stratonovich case with the help of the quadratic variation. The last remaining problem is then to construct the Lévy area Lg, g) for suitable stochastic processes g. We construct it for certain hypercontractive processes. For continuous martingales that possess sufficiently many moments we give a construction of the Itô iterated integrals that allows us to use them as integrators for our pathwise Itô integral. Below we give some pointers to the literature, and we introduce some basic notations which we will use throughout. In Section 2 we recall some details on Ciesielski s isomorphism, and we give a short overview on rough paths and Young integration. In Section 3 3

4 we develop a paradifferential calculus in terms of Schauder functions, and we examine the different components of Young s integral. In Section 4 we construct the rough path integral based on Schauder functions. Section 5 develops the pathwise Itô integral. In Section 6 we construct the Lévy area for suitable stochastic processes. And in Section 7 we apply our integral to solve both Itô type and Stratonovich type SDEs in a pathwise way. Relevant literature Starting with the Lévy-Ciesielski construction of Brownian motion, Schauder functions have been a very popular tool in stochastic analysis. They can be used to prove in a comparatively easy way that stochastic processes belong to Besov spaces; see for example Ciesielski, Kerkyacharian, and Roynette [CKR93], Roynette [Roy93], and Rosenbaum [Ros9]. Baldi and Roynette [BR92] have used Schauder functions to extend the large deviation principle for Brownian motion from the uniform to the Hölder topology; see also Ben Arous and Ledoux[BL94] for the extension to diffusions, Eddahbi, N zi, and Ouknine[ENO99] for the large deviation principle for diffusions in Besov spaces, and Andresen, Imkeller, and Perkowski [AIP13] for the large deviation principle for a Hilbert space valued Wiener process in Hölder topology. Ben Arous, Grădinaru, and Ledoux [BGL94] use Schauder functions to extend the Stroock-Varadhan support theorem for diffusions from the uniform to the Hölder topology. Lyons and Zeitouni [LZ99] use Schauder functions to prove exponential moment bounds for Stratonovich iterated integrals of a Brownian motion conditioned to stay in a small ball. Gantert [Gan94] uses Schauder functions to associate to every sample path of the Brownian bridge a sequence of probability measures on path space, and continues to show that for almost all sample paths these measures converge to the distribution of the Brownian bridge. This shows that the law of the Brownian bridge can be reconstructed from a single typical sample path. Concerning integrals based on Schauder functions, there are three important references: Roynette [Roy93] constructs a version of Young s integral on Besov spaces and shows that in the one dimensional case the Stratonovich integral FW s)dw s, where W is a Brownian motion, and F C 2, can be defined in a deterministic manner with the help of Schauder functions. Roynette also constructs more general Stratonovich integrals with the help of Schauder functions, but in that case only almost sure convergence is established, where the null set depends on the integrand, and the integral is not a deterministic operator. Ciesielski, Kerkyacharian, and Roynette [CKR93] slightly extend the Young integral of [Roy93], and simplify the proof by developing the integrand in the Haar basis and not in the Schauder basis. They also construct pathwise solutions to SDEs driven by fractional Brownian motions with Hurst index H > 1/2. Kamont [Kam94] extends the approach of [CKR93] to define a multiparameter Young integral for functions in anisotropic Besov spaces. Ogawa [Oga84, Oga85] investigates an integral for anticipating integrands he calls noncausal starting from a Parseval type relation in which integrand and Brownian motion as integrator are both developed by a given complete orthonormal system in the space of square integrable functions on the underlying time interval. This concept is shown to be strongly related to Stratonovich type integrals see Ogawa [Oga85], Nualart, Zakai [NZ89]), and used to develop a stochastic calculus on a Brownian basis with noncausal SDE Ogawa [Oga7]). Rough paths have been introduced by Lyons [Lyo98], see also [Lyo95, LQ96, LQ97] for previous results. Lyons observed that solution flows to SDEs or more generally ordinary differential equations ODEs) driven by rough signals) can be defined in a pathwise, continuous way if paths are equipped with sufficiently many iterated integrals. More precisely, if a path 4

5 has finite p variation for some p 1, then one needs to associate p iterated integrals to it to obtain an object which can be taken as the driving signal in an ODE, such that the solution to the ODE depends continuously on that signal. Gubinelli [Gub4, Gub1] simplified the theory of rough paths by introducing the concept of controlled paths, on which we will strongly rely in what follows. Roughly speaking, a path f is controlled by the reference path g if the small scale fluctuations of f look like those of g. Good monographs on rough paths are [LQ2, LCL7, FV1b, FH14]. Notation and conventions. Throughout the paper, we use the notation a b if there exists a constant c >, independent of the variables under consideration, such that a c b, and we write a b if a b and b a. If we want to emphasize the dependence of c on the variable x, then we write ax) x bx). For a multi-index µ = µ 1,...,µ d ) N d we write µ = µ µ d and µ = µ / µ 1 x 1 µ d x d. DF or F denote the total derivative of F. For k N we denote by D k F the k-th order derivative of F. We also write x for the partial derivative in direction x. 2 Preliminaries 2.1 Ciesielski s isomorphism Let us briefly recall Ciesielski s isomorphism between C α [,1],R d ) and l R d ). The Haar functions H pm,p N,1 m 2 p ) are defined as 2 p, t [ m 1 2, 2m 1 ) p 2, p+1 H pm t) := 2 p, t [ 2m 1, m ) 2 p+1 2, p, otherwise. When completed by H 1, the Haar functions are an orthonormal basis of L 2 [,1],dt). For convencience of notation, we also define H p for p 1. The primitives of the Haar functions are called Schauder functions and they are given by G pm t) := t H pms)ds for t [,1], p N, m 2 p. More explicitly, G t) = t and for p N, 1 m 2 p 2 p/2 t m 1 ) [ 2, t m 1 p 2, 2m 1 ) p 2, p+1 G pm t) = 2 p/2 t m ) [ 2, t 2m 1, m ) p 2 p+1 2, p, otherwise. Since every G pm satisfies G pm ) =, we are only able to expand functions f with f) = in termsofthisfamily G pm ). Therefore, wecomplete G pm ) oncemore, bydefiningg 1 t) := 1 for all t [,1]. To abbreviate notation, we define the times t i pm, i =,1,2, as t pm := m 1 2 p, t 1 pm := 2m 1 2 p+1, t2 pm := m 2 p, for p N and 1 m 2 p. Further, we set t 1 :=, t1 1 :=, t2 1 := 1, and t :=, t 1 := 1, t2 := 1, as well as ti p := for p 1 and i =,1,2. The definition of ti 1 and ti for i 1 is rather arbitrary, but the definition for i = 1 simplifies for example the statement of Lemma 2.1 below. 5

6 For f C[,1],R d ), p N, and 1 m 2 p, we write H pm,df :=2 p [ ) )) ) ))] 2 f t 1 pm f t pm f t 2 pm f t 1 pm [ ) ) )] 2f t 1 pm f t pm f t 2 pm =2 p 2 and H,df := f1) f) as well as H 1,df := f). Note that we only defined G 1 and not H 1. Lemma 2.1. For f: [,1] R d, the function f k := H 1,df G 1 + H,df G + k 2 p p= m=1 H pm,df G pm = k 2 p p= 1m= H pm,df G pm is the linear interpolation of f between the points t 1 1,t1,t1 pm, p k,1 m 2p. If f is continuous, then f k ) converges uniformly to f as k. Ciesielski [Cie6] observed that if f is Hölder-continuous, then the series f k ) converges absolutely and the speed of convergence can be estimated in terms of the Hölder norm of f. The norm C α is defined as where we introduced the notation f s,t f C α := f + sup s<t 1 t s α, f s,t := ft) fs). Lemma 2.2 [Cie6]). Let α,1). A continuous function f : [,1] R d is in C α if and only if sup p,m 2 pα 1/2) H pm,df <. In this case sup2 pα 1/2) H pm,df f α and 1) p,m 2 p f f N 1 = H pm,df G pm f α 2 αn. p=n m= Before we continue, let us slightly change notation. We want to get rid of the factor 2 p/2 in 1), and therefore we define for p N and m 2 p the rescaled functions χ pm := 2 p 2H pm and ϕ pm := 2 p 2G pm, as well as ϕ 1 := G 1 1. Then we have for p N and 1 m 2 p ϕ pm t) = ϕ pm t 1 pm) = 2 p 2 t 1 pm t pm 2 p 2 ds = 2 p 2m 1 2 p+1 2m 2 ) 2 p+1 = 1 2, so that ϕ pm 1 for all p,m. The expansion of f in terms of ϕ pm ) is given by f k = k p= 2 p m= f pmϕ pm, where f 1 := f1), and f := f1) f) and for p N and m 1 f pm := 2 p χ pm,df = 2f t 1 pm) f t pm ) f t 2 pm ) = ft pm,t 1 pm f t 1 pm,t 2 pm. We write χ pm,df := 2 p f pm for all values of p,m), despite not having defined χ 1. 6

7 Definition 2.3. For α > and f: [,1] R d the norm α is defined as f α := sup2 pα f pm, pm and we write C α := C α R d ) := { } f C[,1],R d ) : f α <. The space C α is isomorphic to l R d ), in particular it is a Banach space. For α,1), Ciesielski s isomorphism Lemma 2.2) states that C α = C α [,1],R d ). Moreover, it can be shown that C 1 is the Zygmund space of continuous functions f satisfying 2fx) fx+h) fx h) h. But for α > 1, there is no reasonable identification of C α with a classical function space. For example if α 1,2), the space C α [,1],R d ) consists of all continuously differentiable functions f with α 1) Hölder continuous derivative Df. Since the tent shaped functions ϕ pm are not continuously differentiable, even an f with a finite Schauder expansion is generally not in C α. The a priori requirement of f being continuous can be relaxed, but not much. Since the coefficients f pm ) evaluate the function f only in countably many points, a general f will not be uniquely determined by its expansion. But for example it would suffice to assume that f is càdlàg. Littlewood-Paley notation. We will employ notation inspired from Littlewood-Paley theory. For p 1 and f C[,1]) we define 2 p p f := pm ϕ pm and S p f := m=f q f. q p We will occasionally refer to p f) as the Schauder blocks of f. Note that C α = {f C[,1],R d ) : 2 pα p f ) p l < }. 2.2 Young integration and rough paths Here we present the main concepts of Young integration and of rough path theory. The results presented in this section will not be applied in the remainder of this chapter, but we feel that it could be useful for the reader to be familiar with the basic concepts of rough paths, since it is the main inspiration for the constructions developed below. Young s integral [You36] allows to define fdg for f C α, g C β, and α+β > 1. More precisely, let f C α and g C β be given, let t [,1], and let π = {t,...,t N } be a partition of [,t], i.e. = t < t 1 < < t N = t. Then it can be shown that the Riemann sums ft k )gt k+1 ) gt k )) := t k π N 1 k= ft k )gt k+1 ) gt k )) converge as the mesh size max k=,...,n 1 t k+1 t k tends to zero, and that the limit does not depend on the approximating sequence of partitions. We denote the limit by t fs)dgs), 7

8 and we define t s fr)dgr) := t fr)dgr) s fr)dgr). The function t t fs)dgs) is uniquely characterized by the fact that t fr)dgr) fs)gt) gs)) t s α+β f α g β s for all s,t [,1]. The condition α+β > 1 is sharp, in the sense that there exist f,g C 1/2, and a sequence of partitions π n ) n N with mesh size going to zero, for which the Riemann sums t k π n ft k )gt k+1 ) gt k )) do not converge as n tends to. The condition α + β > 1 excludes one of the most important examples: we would like to take g as a sample path of Brownian motion, and f = Fg). Lyons theory of rough paths [Lyo98] overcomes this restriction by stipulating the existence of basic integrals and by defining a large class of related integrals as their functionals. Here we present the approach of Gubinelli [Gub4]. Let α 1/3,1) and assume that we are given two functions v,w C α, as well as an associated Riemann integral I v,w s,t = t vr)dwr) that satisfies the estimate s Φ v,w s,t := Iv,w s,t vs)w s,t t s 2α. 2) The remainder Φ v,w is often incorrectly) called the area of v and w. This name has its origin in the fact that its antisymmetric part 1/2Φ v,w s,t Φ w,v s,t ) corresponds to the algebraic area spanned by the curve vr),wr)) : r [s,t]) in the plane R 2. If α 1/2, then the integral I v,w cannot be constructed using Young s theory of integration, and also I v,w is not uniquely characterized by 2). But let us assume nonetheless that we are given such an integral I v,w satisfying 2). A function f C α is controlled by v C α if there exists f v C α, such that for all s,t [,1] f s,t f v s v s,t t s 2α. 3) Propostion 2.4 [Gub4], Theorem 1). Let α > 1/3, let v,w C α, and let I v,w satisfy 2). Let f and g be controlled by v and w respectively, with derivatives f v and g w. Then there exists a unique function If,g) = fs)dgs) that satisfies for all s,t [,1] If,g) s,t fs)g s,t f v s)g w s)φ v,w s,t t s 3α. If π n ) is a sequence of partitions of [,t], with mesh size going to zero, then ) If,g)t) = lim ft k )g tk,t k+1 +ft v k gt w k Φ v,w t k,t k+1. n t k π n The integral If, g) coincides with the Riemann-Stieltjes integral and with the Young integral, whenever these are defined. Moreover, the integral map is self-consistent, in the sense that if we consider v and w as paracontrolled by themselves, with derivatives v v = w w 1, then Iv,w) = I v,w. The only remaining problem is the construction of the integral I v,w. This is usually achieved with probabilistic arguments. If v and w are Brownian motions, then we can for example use Itô or Stratonovich integration to define I v,w. Already in this simple example we see that the integral I v,w is not unique if v and w are outside of the Young regime. It is possible to go beyond α > 1/3 by stipulating the existence of higher order iterated integrals. For details see [Gub1] or any book on rough paths, such as [LQ2, LCL7, FV1b, FH14]. 8

9 3 Paradifferential calculus and Young integration In this section we develop the basic tools that will be required for our rough path integral in terms of Schauder functions, and we study Young s integral and its different components. 3.1 Paradifferential calculus with Schauder functions Here we introduce a paradifferential calculus in terms of Schauder functions. Paradifferential calculus is usually formulated in terms of Littlewood-Paley blocks and was initiated by Bony [Bon81]. For a gentle introduction see [BCD11]. We will need to study the regularity of p,m u pmϕ pm, where u pm are functions and not constant coefficients. For this purpose we define the following space of sequences of functions. Definition 3.1. If u pm : p 1, m 2 p ) is a family of affine functions of the form u pm : [t pm,t 2 pm] R d, we set for α > u pm ) A α := sup2 pα u pm, p,m where it is understood that u pm := max t [t pm,t 2 pm ] u pm t). The space A α := A α R d ) is then defined as { } A α := u pm ) p 1, m 2 p : u pm C[t pm,t 2 pm],r d ) is affine and u pm ) A α <. In Appendix A we prove the following regularity estimate: Lemma 3.2. Let α,2) and let u pm ) A α. Then p,m u pmϕ pm C α, and α u pm ϕ pm u pm ) A α. p,m Let us introduce a paraproduct in terms of Schauder functions. Lemma 3.3. Let β,2), let v C[,1],LR d,r n )), and w C β R d ). Then π < v,w) := S p 1 v p w C β R n ) and π < v,w) β v w β. 4) p= Proof. We have π < v,w) = p,m u pmϕ pm with u pm = S p 1 v) [t pm,t 2 pm ] w pm. For every p,m), the function S p 1 v) [t pm,t 2 pm ] is the linear interpolation of v between t pm and t 2 pm. As S p 1 v) [t pm,t 2 pm]w pm 2 pβ v w β, the statement follows from Lemma 3.2. Remark 3.4. If v C α R) and w C β R), we can decompose the product vw into three components, vw = π < v,w)+π > v,w)+π v,w), where π > v,w) := π > w,v) and π v,w) := p pv p w, and we have the estimates π > v,w) α v α w, and π v,w) α+β v α w β whenever α+β,2). However, we will not use this. 9

10 The paraproduct allows us to paralinearize nonlinear functions. We allow for a smoother perturbation, which will come in handy when constructing global in time solutions to SDEs. Propostion 3.5. Let α,1/2), β,α], let v C α R d ), w C α+β, and F C 1+β/α b R d,r). Then Fv +w) π < DFv +w),v) α+β F 1+β/α C 1+ v α ) 1+β/α 1+ w α+β ). 5) b If F C 2+β/α b, then Fv) π < DFv),v) depends on v in a locally Lipschitz continuous way: Fv) π < DFv),v) Fu) π < DFu),u)) α+β F 2+β/α C 1+ v α + u α ) 1+β/α v u α. 6) b Proof. First note that Fv+w) F, which implies therequired estimate for p,m) = 1,) and p,m) =,). For all other values of p,m) we apply a Taylor expansion: Fv +w)) pm = DFvt 1 pm)+wt 1 pm))v pm +R pm, where R pm 2 pα+β) F 1+β/α C v 1+β/α α + w α+β ). Subtracting π < DFv),v) gives b Fv +w) π < DFv +w),v) = [DFvt 1 pm)+wt 1 pm)) S p 1 DFv +w)) [t pm,t 2 pm ] ]v pm ϕ pm +R. pm Now S p 1 DFv+w)) [t pm,t 2 pm] is the linear interpolation of DFv+w) between t pm and t 2 pm, so according to Lemma 3.2 it suffices to note that [DFvt 1 pm )+wt1 pm )) S p 1DFv +w)) [t pm,t 2 pm ] ]v pm 2 pβ DFv +w) C β2 pα v α 2 pα+β) F 1+β/α C 1+ v α + w α ) β/α v α. b The local Lipschitz continuity is shown in the same way. Remark 3.6. Since v has compact support, it actually suffices to have F C 1+β/α without assuming boundedness. Of course, then the estimates in Proposition 3.5 have to be adapted. Remark 3.7. The same proof shows that if f is controlled by v in the sense of Section 2.1, i.e. f s,t = f v s)v s,t +R s,t with f v C α and R s,t R 2α t s 2α, then f π < f v,v) C 2α. 3.2 Young s integral and its different components In this section we construct Young s integral using the Schauder expansion. If v C α and w C β, then we formally define vs)dws) := v pm w qn ϕ pm s)dϕ qn s) = p,m q,n p,q p vs)d q ws). We show that this definition makes sense provided that α + β > 1, and we identify three components of the integral that behave quite differently. This will be our starting point towards an extension of the integral beyond the Young regime. In a first step, let us calculate the iterated integrals of Schauder functions. 1

11 Lemma 3.8. Let p > q. Then 1 ϕ pm s)dϕ qn s) = 2 p 2 χ qn t pm ) 7) for all m,n. If p = q, then 1 ϕ pms)dϕ pn s) =, except if p = q =, in which case the integral is bounded by 1. If p < q, then for all m,n) we have 1 If p = 1, then the integral is bounded by 1. ϕ pm s)dϕ qn s) = 2 q 2 χ pm t qn ). 8) Proof. The cases p = q and p = 1 are easy, so let p > q. Since χ qn χ qn t pm ) on the support of ϕ pm, we have 1 ϕ pm s)dϕ qn s) = χ qn t pm) 1 If p < q, then integration by parts and 7) imply 8). ϕ pm s)ds = χ qn t pm)2 p 2. Next we estimate the coefficients of iterated integrals in the Schauder basis. Lemma 3.9. Let i,p 1, q, j 2 i, m 2 p, n 2 q. Then 2 i χ ij,d ϕ pm χ qn ds) 2 2i p q)+p+q, 9) except if p < q = i. In this case we only have the worse estimate 2 i χ ij,d ϕ pm χ qn ds) 1. 1) Proof. We have χ 1,d ϕ pmχ qn ds) = for all p,m) and q,n). So let i. If i < p q, then χ ij is constant on the support of ϕ pm χ qn, and therefore Lemma 3.8 gives 2 i χ ij,ϕ pm χ qn ϕ pm,χ qn 2 p+q 2p q) = 2 2i p q)+p+q. Now let i > q. Then χ qn is constant on the support of χ ij, and therefore another application of Lemma 3.8 implies that 2 i χ ij,ϕ pm χ qn 2 i 2 q 2 p+i 2p i) = 2 2i p q)+p+q. The only remaining case is i = q p, in which t 2 2 i χ ij,ϕ pm χ qn 2 i ij ϕ pm s)ds ϕ pm 1. t ij 11

12 Corollary 3.1. Let i,p 1 and q. Let v C[,1],LR d,r n )) and w C[,1],R d ). Then i p vs)d q ws)) 2 i p q) i+p+q p v q w, 11) except if i = q > p. In this case we only have the worse estimate i Proof. The case i = 1 is easy, so let i. We have i p vs)d q ws)) p v q w. 12) ) p vs)d q ws) = j,m,n v pm w qn 2 i χ ij,ϕ pm χ qn ϕ ij. For fixed j, there are at most 2 i p q) i non-vanishing terms in the double sum. Hence, we obtain from Lemma 3.9 that v pm w qn 2 i χ ij,ϕ pm χ qn ϕ ij 2 i p q) i p v q w 2 2i p q)+p+q +1 i=q>p ) m,n = 2 i p q) i+p+q +1 i=q>p ) p v q w. Corollary Let i,p,q 1. Let v C[,1],LR d,r n )) and w C[,1],R d ). Then for p q i we have i p v q w) 2 i p q) i+p+q p v q w, 13) except if i = q > p or i = p > q, in which case we only have the worse estimate If p > i or q > i, then i p v q w). i p v q w) p v q w. 14) Proof. The case p = 1 or q = 1 is easy. Otherwise we apply integration by parts and note that the estimates 11) and 12) are symmetric in p and q. If for example p > i, then p v)t k ij ) = for all k,j, which implies that i p v q w) =. The estimates 11) and 12) allow us to identify different components of the integral vs)dws). More precisely, 12) indicates that the series p<q pvs)d q ws) is rougher than the remainder p q pvs)d q ws). Integration by parts gives p<q p vs)d q ws) = π < v,w) p<q v pm w qn ϕ qn s)dϕ pm s). This motivates us to decompose the integral into three components, namely p,q m,n p vs)d q ws) = Lv,w)+Sv,w)+π < v,w). 12

13 Here L is defined as the antisymmetric Lévy area we will justify the name below by showing that L is closely related to the Lévy area of certain dyadic martingales): The symmetric part S is defined as Lv,w) := v pm w qn v qn w pm ) ϕ pm dϕ qn p>q m,n = ) p vds p 1 w ds p 1 v) p w. p Sv,w) := m,n 1 = m,n 1 v m w n ϕ m dϕ n + v pm w pm ϕ pm dϕ pm p 1 m v m w n ϕ m dϕ n + 1 p v p w, 2 and π < is the paraproduct defined in 4). As we observed in Lemma 3.3, π < v,w) is always well defined, and it inherits the regularity of w. Let us study S and L. Lemma Let α,β,1) be such that α+β > 1. Then L is a bounded bilinear operator from C α C β to C α+β. Proof. We only argue for p pvds p 1 w, the term ds p 1v) p w can be treated with the same arguments. Corollary 3.1 more precisely 11)) implies that i p ) p vds p 1 w ) i p vd q w + ) i p vd q w p i q<p p>i q<p p i q<p2 2i+p+q 2 pα v α 2 qβ w β + ) 2 i+q 2 pα v α 2 qβ w β p>i q<p α+β 2 iα+β) v α w β, whereweused1 α < and1 β < andforthesecondserieswealsousedthatα+β > 1. Unlike the Lévy area L, the symmetric part S is always well defined. It is also smooth. Lemma Let α,β,1). Then S is a bounded bilinear operator from C α C β to C α+β. Proof. This is shown using the same arguments as in the proof of Lemma In conclusion, the integral consists of three components. The Lévy area Lv, w) is only defined if α+β > 1, but then it is smooth. The symmetric part Sv,w) is always defined and smooth. And the paraproduct π < v,w) is always defined, but it is rougher than the other components. To summarize: p 1 13

14 Theorem 3.14 Young s integral). Let α,β,1) be such that α+β > 1, and let v C α and w C β. Then the integral Iv,dw) := p,q satisfies Iv,dw) β v α w β and Lévy area and dyadic martingales p vd q w = Lv,w)+Sv,w)+π < v,w) C β Iv,dw) π < v,w) α+β v α w β. 15) Here we show that the Lévy area Lv,w)1) can be expressed in terms of the Lévy area of suitable dyadic martingales. To simplify notation, we assume that v) = w) =, so that we do not have to bother with the components v 1 and w 1. We define a filtration F n ) n on [,1] by setting F n = σχ pm : p n, m 2 p ), we set F = n F n, and we consider the Lebesgue measure on [,1],F). On this space, the process M n = n 2 p p= m= χ pm, n N, is a martingale. For any continuous function v : [,1] R with v) =, the process M v n = n 2 p p= m= 2 p χ pm,dv χ pm = n 2 p p= m= v pm χ pm = t S n v, n N, is a martingale transform of M, and therefore a martingale as well. Since it will be convenient later, we also define F 1 = {,[,1]} and M 1 v = for every v. Assume now that v and w are continuous real-valued functions with v) = w) =, and that the Lévy area Lv,w)1) exists. Then it is given by Lv,w)1) = = = = p 1 v pm w qn v qn w pm ) p= q= m,n 1 ϕ pm s)χ qn s)ds p 1 1 v pm w qn v qn w pm )2 p χ qn s)1 [t pm,t 2 pm ) s)ds ϕ pm,1 p= q= m,n p 1 1 v pm w qn v qn w pm )2 p χ qn s)χ 2 pm s)ds2 p 2 p= q= m,n p p 2 p= q= m,n m v pm w qn v qn w pm )χ qn s)χ pm s)χ pm s)ds, where in the last step we used that χ pm and χ pm have disjoint support for m m. The p th Rademacher function or square wave) is defined for p 1 as r p t) := 2 p m =1 2 p χ pm t). 14

15 The martingale associated to the Rademacher functions is given by R := and R p := p k=1 r k for p 1. Let us write Mp v = Mv p Mv p 1 and similarly for Mw and R and all other discrete time processes that arise. This notation somewhat clashes with the expression p v for the dyadic blocks of v, but we will only use it in the following lines, where we do not directly work with dyadic blocks. The quadratic covariation of two dyadic martingales is defined as [M,N] n := n k= M k N k, and the discrete time stochastic integral is defined as M N) n := n k= M k 1 N k. Writing E ) for the integral 1 ds, we obtain Lv,w)1) = = = p 1 2 p 2 E Mp M v q w R p Mq M v p w ) R p p= q= 2 p 2 E Mp 1 w ) ) Mv p Mv p 1 Mw p Rp p= p= ) 2 p 2 E [M w M v M v M w,r] p. Hence, Lv,w)1) is closely related to the Lévy area 1/2M w M v M v M w ) of the dyadic martingale M v,m w ). 4 Paracontrolled paths and pathwise integration beyond Young In this section we construct a rough path integral in terms of Schauder functions. 4.1 Paracontrolled paths We observed in Section 3 that for w C α and F C 1+β/α b we have Fw) π < DFw),w) C α+β. In Section 3.2 we observed that if v C α, w C β and α + β > 1, then the Young integral Iv,dw) satisfies Iv,dw) π < v,w) C α+β. Hence, in both cases the function under consideration can be written as π < f w,w) for a suitable f w, plus a smooth remainder. We make this our definition of paracontrolled paths: Definition 4.1. Let α > and v C α R d ). For β,α] we define { } Dv β := DvR β n ) := f,f v ) C α R n ) C β LR d,r n )) : f = f π < f v,v) C α+β R n ). If f,f v ) D β v, then f is called paracontrolled by v. The function f v is called the derivative of f with respect to v. Abusing notation, we write f D β v when it is clear from the context what the derivative f v is supposed to be. We equip D β v with the norm f v,β := f v β + f α+β. If v C α and f, fṽ) D β ṽ, then we also write d D βf, f) := f v fṽ β + f f α+β. Example 4.2. Let α,1) and v C α. Then Proposition 3.5 shows that Fv) Dv β for every F C 1+β/α b, with derivative DFv). 15

16 Example 4.3. Let α+β > 1 and v C α, w C β. Then by 15), the Young integral Iv,dw) is in Dw α, with derivative v. Example 4.4. If α+β < 1and v C α, then f,f v ) Dv β if and only if f s,t fsv v s,t t s α+β and in that case f v +sup s t fs,t v t s β +sup f s,t fs vv s,t s t t s α+β f v,β 1+ v α ). Indeed we have f v sv s,t π < f v,v) s,t t s α+β f v β v α, which can be shown using similar arguments as for Lemma B.2 in [GIP12]. In other words, for α,1/2) the space D α v coincides with the space of controlled paths defined in Section 2.2. The following commutator estimate, the analog of Theorem 2.3 of [Bon81] in our setting, will be useful for establishing some stability properties of D β v. Lemma 4.5. Let α,β,1), and let u C[,1],LR n ;R m )), v C α LR d ;R n )), and w C β R d ). Then Proof. We have π < u,π < v,w)) π < uv,w) α+β u v α w β. π < u,π < v,w)) π < uv,w) = p,ms p 1 uπ < v,w)) pm S p 1 uv)w pm )ϕ pm and [S p 1 uπ < v,w)) pm S p 1 uv)w pm ] [t pm,t 2 pm ] is affine. By Lemma 3.2 it suffices to control [S p 1 uπ < v,w)) pm S p 1 uv)w pm ] [t pm,t 2 pm]. The cases p,m) = 1,) and p,m) =,) are easy, so let p and m 1. For r < q < p we denote by m q and m r the unique index in generation q and r respectively for which χ pm ϕ qmq and similarly for r. We apply Lemma 3.9 to obtain for q < p S q 1 v q w) pm = v rmr w qmq 2 p χ pm,dϕ rmr ϕ qmq ) r<q = v rmr w qmq 2 p χ pm,χ rmr ϕ qmq +ϕ rmr χ qmq r<q v α w β 2 rα 2 qβ 2 p 2 2p+r+p+q 2 2p+q2 α β) v α w β. r<q Hence S p 1 u q<ps q 1 v q w) pm ) [t pm,t 2 pm ] u v α w β 2 pα+β). If p < q, then q wt k pm) = for all k and m, and therefore S q 1 v q w) pm =, so that it only remains to bound [S p 1 us p 1 v p w) pm S p 1 uv)w pm ] t pm,t 2 pm ]. We have p wt pm) = p wt 2 pm) = and p wt 1 pm) = w pm /2. On[t pm,t 2 pm], thefunctions p 1 v isgivenbythelinear interpolation of vt pm ) and vt2 pm ), and therefore S p 1v p w) pm = 1 2 vt pm ) + vt2 pm ))w pm, 16

17 leading to [S p 1 us p 1 v p w) pm S p 1 uv)w pm ] [t pm,t 2 pm ] [ w pm ut pm )+ t ) pm vt pm )+vt 2 pm t 2 pm t t pm,t pmu ) 2 pm 2 uv)t pm)+ t pm [t t 2 t pm pmuv) t pm,tpm)] 2 u v α w β 2 pα+β), where the last step follows by rebracketing. pm,t2 pm ] As a consequence, we can show that paracontrolled paths are stable under the application of smooth functions. Corollary 4.6. Let α,1), β,α], v C α, and f Dv β with derivative f v. Let F C 1+β/α b. Then Ff) Dv β with derivative DFf)f v, and Ff) v,β F 1+β/α C 1+ v α ) 1+β/α 1+ f v,β )1+ f v ) 1+β/α. b Moreover, there exists a polynomial P which satisfies for all F C 2+β/α b, ṽ C α, f D β ṽ, and M = max{ v α, ṽ α, f v,β, f ṽ,β } the bound d D βff),f f)) PM) F 2+β/α C d D βf, f)+ u ũ α ). b Proof. The estimate for DFf)f v β is straightforward. For the remainder we apply Proposition 3.5 and Lemma 4.5 to obtain Ff) α+β Ff) π < DFf),f) α+β + π < DFf),f ) α+β + π < DFf),π < f v,v)) π < DFf)f v,v) α+β F 1+β/α C 1+ π < f v,v) α ) 1+β/α 1+ f α+β ) b + F C 1 f b v,β + F C 1 f v b β v α F 1+β/α C 1+ f v ) 1+β/α 1+ v α ) 1+β/α 1+ f v,β ). b The difference Ff) F f) is treated in the same way. When solving differential equations it will be crucial to have a bound which is linear in f v,β. The superlinear dependence on f v will not pose any problem as we will always have f v = F f) for some suitable f, so that for bounded F we get Ff) v,β F,v 1+ f v,β. 4.2 A basic commutator estimate Here we prove the commutator estimate which will be the main ingredient in the construction of the integral If,dg), wheref is paracontrolled by v andg is paracontrolled by w, andwhere we assume that the integral Iv, dw) exists. 17

18 Propostion 4.7. Let α,β,γ,1), and assume that α + β + γ > 1 and β + γ < 1. Let f C α, v C β, and w C γ. Then the commutator Cf,v,w) := Lπ < f,v),w) If,dLv,w)) 16) := lim [LS Nπ < f,v)),s N w) If,dLS N v,s N w))] N [ = lim p π < f,v))s)d q ws) d q π < f,v)))s) p ws) N p N q<p fs) p vs)d q ws) converges in C α+β+γ ε for all ε >. Moreover, Cf,v,w) α+β+γ f α v β w γ. Proof. We only argue for the first difference in 16), i.e. for X N := [ p π < f,v))s)d q ws) p N q<p ) ] fs)d q v)s) p ws) ] fs) p vs)d q ws). 17) The second difference can be handled using the same arguments. First we prove that X N ) converges uniformly, then weshowthat X N α+β+γ stays uniformlybounded. Thiswill imply the desired result, since bounded sets in C α+β+γ are relatively compact in C α+β+γ ε. To prove uniform convergence, note that X N X N 1 = [ ] N π < f,v))s)d q ws) fs) N vs)d q ws) q<n = [ N i f j v)s)d q ws) q<n j N i<j ] j i f N v)s)d q ws), 18) j N i j where for the second term it is possible to take the infinite sum over j outside of the integral because j jg converges uniformly to g and because q w is a finite variation path. We also used that N i f j v) = whenever i > N or j > N. Only very few terms in 18) cancel. Nonetheless these cancellations are crucial, since they eliminate most terms for which we only have the worse estimate 14) in Corollary We obtain X N X N 1 = q<n j<n i<j q<n j>n q<n j>n i<j N i f j v)s)d q ws) q<n j i f N v)s)d q ws) N N f N v)s)d q ws) j j f N v)s)d q ws). 19) 18

19 Note that t q w 2 q q w. Hence, an application of Corollary 3.11, where we use 13) for the first three terms and 14) for the fourth term, yields X N X N 1 f α v β w γ [ q<n 2 2N+i+j 2 iα 2 jβ 2 q1 γ) j<n i<j + 2 Nα+β) 2 q1 γ) + q<n q<n + ] 2 jα 2 Nβ 2 q1 γ) q<n j>n 2 2j+i+N 2 iα 2 Nβ 2 q1 γ) j>n i<j f α v β w γ 2 Nα+β+γ 1), 2) where in the last step we used α,β,γ < 1. Since α + β + γ > 1, this gives us the uniform convergence of X N ). Next let us show that X N α+β+γ f α v β w γ for all N. Similarly to 19) we obtain for n N n X N = p N n[ q<p j<p i<j j>p and therefore by Corollary 3.1 i j p i f j v)s)d q ws) j i f p v)s)d q ws) ], p p f p v)s)d q ws) n X N p [ 2 n p) n+p+q p i f j v) q w q<p j<p i<j +2 n p) n+p+q p p f p v) q w + 2 n j) n+j+q j i f p v) q w ]. j>p i j Now we apply Corollary 3.11, where for the last term we distinguish the cases i < j and i = j. Using that 1 γ >, we get [ n X N f α v β w γ 2 p1 γ) 2 n p) n+p 2 2p 2 i1 α) 2 j1 β) p j<p i<j f α v β w γ 2 nα+β+γ), where we used both that α+β +γ > 1 and that β +γ < n p) n+p 2 pα 2 pβ + 2 n j) n+j 2 2j+i1 α)+p1 β) j>p i<j + ] 2 n j) n+j 2 jα pβ j>p 19

20 Remark 4.8. If β+γ = 1, we can apply Proposition 4.7 with β ε to obtain that Cf,v,w) C α+β+γ ε for every sufficiently small ε >. If β +γ > 1, then we are in the Young setting and there is no need to introduce the commutator. For later reference, we collect the following result from the proof of Proposition 4.7: Lemma 4.9. Let α,β,γ,f,v,w be as in Proposition 4.7. Then Cf,v,w) LS N π < f,v)),s N w) If,dLS N v,s N w)) 2 Nα+β+γ 1) f α v β w γ. Proof. Simply sum up 2) over N. 4.3 Pathwise integration for paracontrolled paths In this section we apply the commutator estimate to construct the rough path integral under the assumption that the Lévy area exists for a given reference path. Theorem 4.1. Let α 1/3,1), β,α] and assume that 2α+β > as well as α+β 1. Let v C α R d ) and assume that the Lévy area Lv,v) := lim N LSN v k,s N v l ) ) 1 k d,1 l d converges uniformly and that sup N LS N v,s N v) 2α <. Let f D α vlr d,r m )). Then IS N f,ds N v) converges in C α ε for all ε >. Denoting the limit by If,dv), we have If,dv) α f v,β v α + v 2 α + Lv,v) 2α ). Moreover, If,dv) D α v with derivative f and If,dv) v,α f v,β 1+ v 2 α + Lv,v) 2α ). Proof. If β +γ > 1, everything follows from the Young case, Theorem 3.14, so let β +γ < 1. We decompose IS N f,ds N v) = SS N f,s N v)+π < S N f,s N v)+ls N f,s N v) +[LS N π < f v,v),s N v) If v,dls N v,s N v))]+if v,dls N v,s N v)). Convergence then follows from Proposition 4.7 and Theorem The limit is given by If,dv) = Sf,v)+π < f,v)+lf,v)+cf v,v,v)+if v,dlv,v)), from where we easily deduce the claimed bounds. Remark Since If,dv) = lim N S NfdS N v, the integral is a local operator in the sense that If,dv) is constant on every interval [s,t] for which f [s,t] =. In particular we can estimate If,dv) [,t] using only f [,t] and f v [,t]. For fixed v and Lv,v), the map f If,dv) is linear and bounded from D β v to D α v, and this is what we will need to solve differential equations driven by v. But we can also estimate the speed of convergence of IS N f,ds N v) to If,dv), measured in uniform distance rather than in C α : 2

21 Corollary Let α 1/3,1/2] and let β,v,f be as in Theorem 4.1. Then we have for all ε,2α+β 1) IS N f,ds N v) If,dv) ε 2 N2α+β 1 ε) f v,β v α + v 2 ) α + f v β LS N v,s N v) Lv,w) 2α ε. Proof. We decompose IS N f,ds N v) as described in the proof of Theorem 4.1. This gives us for example the term π < S N f f,s N v)+π < f,s N v v) ε S N f f v α + f f α S N v v ε for all ε >. From here it is easy to see that π < S N f f,s N g)+π < f,s N g g) 2 Nα ε) f α v α 2 Nα ε) f v,β v α + v 2 α). But now β α 1/2 and therefore α 2α+β 1. Let us treat one of the critical terms, say LS N f,s N v) Lf,v). Since 2α+β ε > 1, we can apply Lemma 3.12 to obtain LS N f,s N v) Lf,v) LS N f f,s N v) 1+ε + Lf,S N v v) 1+ε Lemma 4.9 gives ε S N f f 1+ε α v α + f α+β S N v v 1+ε α β 2 Nα+β 1+ε α)) f α+β v α +2 Nα 1+ε α β)) f α+β v α 2 N2α+β 1 ε) f α+β v α. LS N π < f v,v),s N v) Lπ < f v,v),v) 2 N2α+β 1) f v β v 2 α + If v,dls N v,s N v)) If v,dlv,v)). The second term on the right hand side can be estimated using the continuity of the Young integral, and the proof is complete. Remark In Lemma 4.9 we saw that the rate of convergence of LS N π < f v,v),s N v) If v,dls N v,s N v)) Lπ < f v,v),v) If v,dlv,v))) is in fact 2 N2α+β 1) when measured in uniform distance, and not just 2 N2α+β 1 ε). It is possible to show that this optimal rate is attained by the other terms as well, so that IS N f,ds N v) If,dv) 2 N2α+β 1) f v,β v α + v 2 ) α + f v β LS N v,s N w) Lv,w) 2α ε. Since this requires a rather lengthy calculation, we decided not to include the arguments here. Since we approximate f and g by the piecewise smooth functions S N f and S N g when defining the integral If, dg), it is not surprising that we obtain a Stratonovich type integral: Propostion Let α 1/3,1) and v C α R d ). Let ε > be such that 2+ε)α > 1 and let F C 2+ε R d,r). Then for all t [,1]. Fvt)) Fv)) = IDFv),dv)t) := lim N IS NDFv),dS N v)t) 21

22 Proof. The function S N v is Lipschitz continuous, so that integration by parts gives FS N vt)) FS N v)) = IDFS N v),ds N v)t). The left hand side converges to Fvt)) Fv)). It thus suffices to show that IS N DFv) DFS N v),ds N v) converges to zero. By continuity of the Young integral, Theorem 3.14, it suffices to show that lim N S N DFv) DFS N v) α1+ε ) = for all ε < ε. Recall that S N v is the linear interpolation of v between the points t 1 pm) for p N and m 2 p, and therefore p DFS N v) = p DFv) = p S N DFv) for all p N. For p > N and 1 m 2 p we apply a first order Taylor expansion to both terms and use the ε Hölder continuity of D 2 F to obtain [S N DFv) DFS N v)] pm C F 2 pα1+ε) S N v α for a constant C F >. Therefore, we get for all ε ε which completes the proof. S N DFv) DFS N v) α1+ε ) F 2 Nαε ε ) v α, Remark Note that here we did not need any assumption on the area Lv,v). The reason are cancellations that arise due to the symmetric structure of the derivative of DF, the Hessian of F. Proposition 4.14 was previously obtained by Roynette [Roy93], except that there v is assumed to be one dimensional and in the Besov space B 1/2 1,. 5 Pathwise Itô integration In the previous section we saw that our pathwise integral If,dv) is of Stratonovich type, i.e. it satisfies the usual integration by parts rule. But in applications it may be interesting to have an Itô integral. Here we show that a slight modification of If,dv) allows us to treat non-anticipating Itô-type integrals. A natural approximation of a non-anticipating integral is given for k N by I Itô k 2k f,dv)t) := = m=1 2 k ft km )vt2 km t) vt km t)) f pm v qn ϕ pm t km )ϕ qnt 2 km t) ϕ qnt km t)). m=1 p,q m,n Let us assume for the moment that t = m2 k for some m 2 k. In that case we obtain for p k or q k that ϕ pm t km )ϕ qnt 2 km t) ϕ qnt km t)) =. For p,q < k, both ϕ pm and ϕ qn are affine functions on [t km t,t2 km t], and for affine u and w and s < t we have us)wt) ws)) = t s ur)dwr) 1 2 [ut) us)][wt) ws)]. 22

23 Hence, we conclude that for t = m2 k I Itô k f,dv)t) = IS k 1f,dS k 1 v)t) 1 2 [f,v] kt), 21) where [f,v] k is the k th dyadic approximation of the quadratic covariation [f,v], i.e. [f,v] k t) := 2 k m=1 [ft 2 km t) ft km t)][vt2 km t) vt km t)]. From now on we study the right hand side of 21) rather than Ik Itô f,dv), which is justified by the following remark. Remark 5.1. Let α,1). If f C[,1]) and v C α, then Ik IS Itô f,dv) k 1 f,ds k 1 v) 1 ) 2 [S k 1f,S k 1 v] k 2 kα f v α. This holds because both functions agree in all dyadic points of the form m2 k, and because between those points the integrals can pick up mass of at most f 2 kα v α. We write [v,v] := [v i,v j ]) 1 i,j d and Lv,v) := Lv i,v j )) 1 i,j d, and similarly for all expressions of the same type. Theorem 5.2. Let α,1/2) and let β α be such that 2α+β > 1. Let v C α R d ) and f DvLR β d ;R n )). Assume that LS k v,s k v)) converges uniformly, with uniformly bounded C 2α norm. Also assume that [v,v] k ) converges uniformly. Then Ik Itô f,dv)) converges uniformly to a limit I Itô f,dv) = If,dv) 1/2[f,v] which satisfies I Itô f,dv) f v,β v α + v 2 α + Lv,v) 2α + [v,v] ), and where the quadrativ variation [f,v] is given by [f,v] = d f v s)d[v,v]s) := j,l=1 ) f ij ) v,l s)d[v j,v l ]s), 22) 1 i n where f ij ) v,l is the l th component of the v derivative of f ij. For ε,3α 1) the speed of convergence can be estimated by I Itô f,dv) Ik Itô f,dv) ε 2 k2α+β 1 ε) f v,β v α + v 2 ) α + f v β LS k 1 v,s k 1 v) Lv,v) 2α + f v [v,v] k [v,v]. Proof. By Remark 5.1, it suffices to show our claims for IS k 1 f,ds k 1 v) 1/2[f,v] k. The statements for the integral IS k 1 f,ds k 1 g) follow from Theorem 4.1 and Corollary So let us us concentrate on the quadratic variation [f,v] k. Recall from Example 4.4 that f D β v if and only if R f s,t = f s,t f v s)w s,t satisfies R f s,t t s α+β. Hence [f,v] i k t) = m = m ft km t,t2 km tv t km t,t2 km t ) i R f t km t,t2 km tv t km t,t2 km t ) i + d j,l=1 f ij ) v,l t km t)vl t km t,t2 km tvj t km t,t2 t. km m 23

24 It is easy to see that the first term on the right hand side is bounded by R f ) i t km t,t2 tv t 2 km km t,t2 km t k2α+β 1) f v,β v α + v 2 α ). m For the second term, let us fix l and j. Then the sum over m is just the integral of f ij ) v,l with respect to the signed measure µ k t = m δ t vj km t km t,t2 km tvl t km t,t2 t. Decomposing µk t km into a positive and negative part as µ k t = 1 [ δ 4 t [v j +v l ) km t km t,t 2 km t]2 δ t [v j v l ) km t km t,t 2 t]2] km m m and similarly for dµ t = d[v j,v l ] t we can estimate 1 f ij ) v,l s)µ k tds) 1 f ij ) v,l s)µ t ds) f v [v i +v j ] k [v i +v j ] + [v i v j ] k [v i v j ] ) f v [v,v] k [v,v], wherewewrite[u] := [u,u]andsimilarlyfor[u] k. Byassumptiontherighthandsideconverges to zero, from where we get the uniform convergence of [f,g] k to [f,g]. Remark 5.3. We calculate the pathwise Itô integral I Itô f,dv) as limit of nonanticipating Riemann sums involving only f and v. This is interesting for applications of mathematical finance, because the integral process has a natural interpretation as capital obtained from investing. The classical rough path integral, see Proposition 2.4, is obtained via compensated Riemann sums that explicitly depend on f v and I Itô v,dv). Remark 5.4. We calculate the pathwise Itô integral I Itô f,dv) as limit of nonanticipating Riemann sums involving only f and v. The classical rough path integral, see Proposition 2.4, is obtained via compensated Riemann sums that depend explicitly on the derivative f v and the iterated integrals of v. For applications in mathematical finance, it is more convenient to have an integral that is the limit of nonanticipating Riemann sums, because this can be interpreted as capital process obtained from investing. Note that [v,v] is always a continuous function of bounded variation, but a priori it is not clear whether it is in C 2α. Under this additional assumption we have the following stronger result. Corollary 5.5. In addition to the conditions of Theorem 5.2, assume that also [v,v] C 2α. Then I Itô f,dv) D α v with derivative f, and I Itô f,dv) v,α f v,β 1+ v 2 α + Lv,v) 2α + [v,v] 2α ). Proof. This is a combination of Theorem 4.1 and the explicit representation 22) together with the continuity of the Young integral, Theorem The term IS k 1 f,ds k 1 v) has the pleasant property that if we want to refine our calculation by passing from k to k + 1, then we only have to add the additional term IS k 1 f,d k v) + I k f,ds k v). For the quadratic variation [f,v] k this is not exactly true. But [f,v] k m2 k ) = [S k 1 f,s k 1 v] k m2 k ) for m =,...,2 k, and there is a recursive way of calculating [S k 1 f,s k 1 v] k : 24