Topics in Stochastic Differential Equations and Rough Path Theory
vorgelegt von Dipl.-Math. Joscha Diehl Von der Fakult¨atII – Mathematik und Naturwissenschaften der Technischen Universit¨atBerlin zur Erlangung des akademischen Grades
Doktor der Naturwissenschaften – Dr. rer. nat. – genehmigte Dissertation
Promotionsausschuss Vorsitzender: Prof. Dr. Dietmar H¨omberg Berichter: Prof. Dr. Peter Friz Prof. Dr. Terry Lyons
Tag der wissenschaftlichen Aussprache: 22. M¨arz2012
Berlin 2012 D 83 Acknowledgment I want to express my deep gratitude towards my advisor Prof. Peter Friz, who set me on an interesting collection of problems. He was always ready to discuss any piece of mathematics and his unyielding enthusiasm was a constant source of joy. I want to thank Prof. Terry Lyons for accepting to be co-examiner of this thesis. Financial support by the DFG through the IRTG “Stochastic Models of Complex Processes” as well as the SPP1324 is gratefully acknowledged. Contents
1 Introduction 2
2 Stochastic differential equations with rough drift 5 2.1 Notation ...... 5 2.2 Main results ...... 5 2.3 The Young case ...... 21 2.4 Applications ...... 32 2.4.1 Robust filtering ...... 33 2.4.2 Pathwise stochastic control ...... 40 2.4.3 Mixed stochastic differential equations ...... 45
3 Backward stochastic differential equations with rough driver 46 3.1 Notation ...... 47 3.2 Main results ...... 47 3.3 The Markovian Setting - Connection to rough PDEs ...... 58 3.4 Connection to BDSDEs ...... 64 3.5 Technical results ...... 66 3.5.1 Comparison for BSDEs ...... 66 3.5.2 Properties of rough flows ...... 70 3.5.3 Comparison for PDEs ...... 71
1 1. Introduction
1 Introduction
Rough path theory, introduced by Lyons [52], is a recent development to treat solutions of differential equations driven by paths that are very irregular in time. More precisely, the theory gives meaning to differential equation of the form
Z t Yt = Y0 + V (Yr)dXr, (1) 0 under appropriate assumptions on the deterministic path X and the vector field V , which go 1 d way beyond the case X ∈ C ([0,T ], R ). Heuristically, the theory can deal with α-H¨oldercontinuous paths of some order α > 0. The regime α > 1/2 is the well understood Young theory ([62, 49]). For α ∈ (1/3, 1/2] it turns d out that the information of the R -valued path alone is not sufficient to build a satisfying theory. An α-H¨olderrough path then takes values in a larger metric space which, in addition to the path itself, contains “second order information” which can be though of as the iterated integrals of the path against itself. 1 This information has to be given exogenously, since it is not canonically determined, and has to satisfy a certain algebraic compatibility condition. On the space of all such paths a rough path “norm” 2 (similar to the usual α-H¨oldernorm d of paths in R ) is introduced, and all elements with finite norm are denoted α-H¨olderrough paths. Distance between paths is measured using this “norm”, giving rise to a rough path metric. 3 There are then at least three different approaches to solutions of rough differential equations (RDEs) of the form (1), which lead essentially to the same results. Lyons’ original idea was to define some sort of rough integration first and consider the integral equation as a fixed point problem. One of the key results is stability: with respect to rough path metric the solution mapping of RDEs is continuous in the driving signal (Lyons’ limit theorem). In particular, if X is a geometric rough path, which by definition means that there exists a sequence of smooth paths Xn converging to X in rough path metric, then in particular the solution to the RDE driven by X is the limit of classical ODE solutions driven by the approximating smooth paths Xn . The approach by Davie [21] defines a solution locally, taking inspiration from (higher order) Euler approximations. In this thesis we shall adopt the approach put forward by Friz and Victoir [33]. It takes Lyons’ limit theorem as starting point, and defines a solution to an RDE driven by a geometric rough path as the limit (if it exists uniquely) of ODE solutions driven by approximating smooth paths. Let us note that the theory can treat many stochastic processes as driving signals, continuous semimartingales being the prime example (see e.g. [51, 33]). Owing to stability properties of rough path theory, many known results from stochastic analysis can then be easily recovered (existence and regularity of stochastic flows, support theorems, large deviation principles, see e.g. [46, 33]). A plethora of new results have also been made possible. Let us mention new numerical algorithms for SDEs [35], SDEs driven by non-semimartingales [17, 33, 40], H¨ormandertype theorems for SDEs driven by Gaussian (rough) paths [11, 3] and ergodicity for SDEs driven by fractional Brownian motion [39]. As we pointed out, an important (or even defining) feature of rough path theory is, that it establishes continuity in the driving signal for time inhomogeneous differential equations with
1For α < 1/3, more and more additional information (more “iterated integrals”) has to be given a priori, related to the size of 1/α. 2The space is not a vector space. 3We note that one can actually choose from a family of metrics; this shall be of no concern here.
2 1. Introduction respect to a certain metric. This thesis investigates whether the principle of stability can be extended to other objects that are classically well-defined when taking a smooth path η as input. In particular we consider the solution mapping for SDEs with an additional drift Z t Z t Z t St = S0 + a(Sr)dr + b(Sr)dBr + c(Sr)dηr, (2) 0 0 0 BSDEs with an additional driver Z T Z T Z T Yt = ξ + f(r, Yr,Zr)dr + H(Yr)dηr − ZrdWr, (3) t t t and parabolic PDES of the form
2 ∂tu + F (t, u, Du, D u) + H(u, Du)η ˙t = 0. (4)
Note that if η is smooth, these equations possess a unique solution (under appropriate as- sumptions on the coefficients). Section 2 covers our results on equation (2). We establish stability in rough path metric with respect to η, which makes it possible to define a solution to an SDE with rough drift for arbitrary geometric rough paths η. Deterministic rough path estimates do not easily mix with stochastic estimates and we tackle the problem using three different approaches. The first one consist in the application of a flow transform, which converts the equation into one that depends on the driving signal in a harmless way; in particular there is no explicit dependence on the derivative anymore. A similar technique is used in the section on BSDEs with rough driver, but the flow transform leads to more technical problems there. The second approach interprets equation (2) directly as an RDE. For this to work one needs to understand η and the Brownian motion B together as a joint rough path. There is a canonical way of doing this though using stochastic integration and formal integration by parts, at least for the case α > 1/3. The main problem in this approach is showing exponential integrability of solutions, a property that we need for applications. Lastly, in the Young regime α > 1/2 we are able to solve the integral equation directly. The main technical difficulty here stems from the fact that the deterministic estimates available for the Young integral are not easily combined with the stochastic (not pathwise) estimates available for the stochastic integral. The motivation and main application for equations of type (2) is a problem in the theory of nonlinear filtering. As conjectured by Lyons [52], we show that in order to do robust filtering for multidimensional observation, correlated with the signal, one needs to treat the observation as an element of a rough path space. This is shown in Section 2.4.1 using SDEs with rough drift. We refer to the introduction of that section for a thorough explanation and a historical account concerning the problem of robust filtering. In Section 2.4.2 we consider an application to (pathwise) stochastic control. The resulting Hamilton-Jacobi-Bellman equation is of main interest here, since it gives a stochastic repre- sentation for a class of “rough partial differential equations” (4). We note that rough path theory has been applied to the field of partial differential equations (PDEs) in several papers [36, 27, 61]. We work with the approach put forward by Friz et al [9, 10, 32], which treats rough PDEs as the limit of PDEs in which the rough path is replaced by an approximating sequence of smooth paths. The aforementioned stability of rough path theory is also a key feature in the theory of viscosity solutions to PDEs ([18]). It is then natural to interpret solutions to the approximating PDEs in the viscosity sense, and this is excatly what Friz et al did. This also turns out to be convenient for this particular application, since viscosity solutions are the right solution concept for PDEs originating from optimal control. Let us
3 1. Introduction mention that we can, in a Brownian setting (that is, when formally replacing η by a Wiener process), recover results by Buckdahn and Ma [8], under somewhat weaker assumptions and significantly shorter proofs. We note that (formally), when plugging in a fractional Brownian motion with Hurst param- eter H > 1/2 for the path η in (2), one gets the solution to a mixed Brownian-fractional Brownian SDE. Such equations are treated for example in [63] (with only one-dimensional Brownian motion and stronger assumptions on the vector fields). This statement is made rigorous in Theorem 33. Backward stochastic differential equations (BSDEs) were introduced by Bismut in 1973. In [4] he applied linear BSDEs to stochastic optimal control. In 1990 Pardoux and Peng [57] then considered non-linear equations. A solution to a BSDE with driver f and random variable 2 ξ ∈ L (FT ) is an adapted pair of processes (Y,Z) in suitable spaces, satisfying
Z T Z T Yt = ξ + f(r, Yr,Zr)ds − ZrdWr, t ≤ T. t t Under appropriate conditions on f and ξ they showed the existence of a unique solution to such an equation. Again we are interested in the (formal) equation (3), when η is a rough path. In the case of η beeing a smooth path, one is back in the classical setting. In Section 3 we will show that under appropriate conditions, if the smooth path converge to a rough path, then the corresponding solutions also converge. We use a flow transformation to prove the limit theorem. As opposed to the case of afore- mentioned SDEs, this leads to technical difficulties. The transformed BSDE falls outside the classical Lipschitz framework, since it is quadratic in the control variable, and we have to use stability results by Kobylanski [43]. We note that formally, when plugging in (the rough path lift of) a Brownian motion for the rough path in (3), one gets the solution to a backward doubly stochastic differential equation (BDSDE), a notion introduced in [55]. This statement is made rigorous in Section 3.4. When the randomness of the parameters (f, ξ) of a BSDE comes from the state of a standard SDE, the system it is called a forward-backward stochastic differential equation (FBSDE). The solutions of FBSDE’s are connected to (viscosity) solutions of partial differential equa- tions (PDEs) in what is sometimes called a nonlinear Feynman-Kac formula, see e.g. [54]. We extend this result to BSDEs with rough drift, thereby arriving at a nonlinear Feynman- Kac formula for a class of rough RDEs in Section 3.3. At the same time we get existence, uniqueness and stability for a class of equations, that has not been treated before in the literature.
4 2. Stochastic differential equations with rough drift
2 Stochastic differential equations with rough drift
The aim of this part is to give meaning to an equation of the form Z t Z t Z t ¯ St = S0 + a(Sr)dr + b(Sr)dBr + c(Sr)dηr, 0 0 0 in the case where η is a rough path. As mentioned in the introduction, we proceed by continuously extending the solution mapping, which is of course classically well defined if η is (the lift of) a smooth path. After introducing some notation, we present and prove the main results in Section 2.2. In Section 2.3 we improve, in the case of a Young driving path, on the regularity assumptions. A specific feature of this case is, that the SDE with rough drift actually satisfies the integral equation. Finally, in Section 2.4 we present some applications.
2.1 Notation
γ 4 m n Lip is the set of γ-Lipschitz functions a : R → R where m and n are chosen according to the context.
[p] dY ∼ d d G (R ) = R ⊕ so(dY ) is the free nilpotent group of step [p] over R . Here [p] denotes the largest integer not larger then p. 5
0,α 0,α−H¨ol [1/α] dY C := C0 ([0, t],G (R )) is the set of geometric α-H¨olderrough paths η : [0, t] → [1/α] d G (R Y ) starting at 0. We shall use the non-homogeneous metric ρα−H¨ol on this space.
0,q−var 0,q−var [q] dY C := C0 ([0, t],G (R )) is the set of geometric q-variation rough paths η : [0, t] → [q] d G (R Y ) starting at 0. ¯ ¯ ¯ ¯ 0 0 ¯ (Ω, F, (Ft)t≥0, P) will be a filtered probability space. Let S = S (Ω) denote the space of d adapted, continuous processes in R S , with the topology of uniform convergence in probability. For q ≥ 1 we denote by Sq = Sq(Ω)¯ the space of processes X ∈ S0 such that
1/q ¯ q ||X||Sq := E[sup |Xs| ] < ∞. s≤t
2.2 Main results
¯ ¯ ¯ ¯ Let (Ω, F, (Ft)t≥0, P) be a filtered probability space carrying a dB-dimensional Brownian motion B¯ and a bounded dS-dimensional random vector S0 independent of B¯. n d n In the following α ∈ (0, 1]. Let η : [0, t] → R Y be smooth paths, such that η → η in 0,α n α-H¨older,for some η ∈ C and let S be a dS-dimensional process which is the unique solution to the classical SDE Z t Z t Z t n n n ¯ n n St = S0 + a(Sr )dr + b(Sr )dBr + c(Sr )dηr , 0 0 0 4In the sense of E. Stein, i.e. bounded k-th derivative for k = 0,..., bγc and γ − bγc-H¨oldercontinuous bγc-th derivative, where bγc is the largest integer strictly smaller then γ. 5This is the ”correct” state space for a geometric 1/p-H¨olderrough path; the space of such paths subject to 1/p-H¨olderregularity (in rough path sense) yields a complete metric space under 1/p-H¨olderrough path metric. Technical details of geometric rough path spaces can be found e.g. in Section 9 of [33].
5 2. Stochastic differential equations with rough drift where we assume that
1 dS 1 dS γ+2 dS (a1) α ∈ (0, 1], γ > 1/α, a ∈ Lip (R ), b1, . . . , bdB ∈ Lip (R ) and c1, . . . , cdY ∈ Lip (R )
1 dS 1 dS (a1’) α ∈ (0, 1], γ > 1/α, a ∈ Lip (R ), b1, . . . , bdB ∈ Lip (R ) and c1, . . . , cdY ∈ γ+3 d Lip (R S )).
1 dS γ dS (a2) α ∈ (1/3, 1/2), γ > 1/α, a ∈ Lip (R ), b1, . . . , bdB ∈ Lip (R ) and c1, . . . , cdY ∈ γ d Lip (R S )).
Theorem 1. Under assumption (a1), (a1’) or (a2), there exists a dS-dimensional process S∞ ∈ S0 such that
Sn → S∞, in S0.
In addition, the limit Ξ(η) := S∞ only depends on η and not on the approximating sequence. Moreover, for all q ≥ 1, η ∈ C0,α it holds that Ξ(η) ∈ Sq and the corresponding mapping Ξ: C0,α → Sq is locally uniformly continuous (and locally Lipschitz under assumption (a1’) or (a2)).
Following Theorem 1, we say that Ξ(η) is a solution of the SDE with rough drift
Z t Z t Z t ¯ Ξ(η)t = S0 + a(Ξ(η)r)dr + b(Ξ(η)r)dBr + c(Ξ(η)r)dηr. (5) 0 0 0 The following result establishes some of the salient properties of solutions of SDEs with rough ˆ ¯ drift. Let (Ω, F, P0) carry, a dY -dimensional Brownian motion Y and let Ω = Ω × Ω be the ˆ ¯ product space, with product measure P := P0 ⊗ P. Let S be the unique solution on this probability space to the SDE
Z t Z t Z t St = S0 + a(Sr)dr + b(Sr)dB¯r + c(Sr) ◦ dYr. (6) 0 0 0 Denote by Y the rough path lift of Y (i.e. the enhanced Brownian Motion over Y ).
Theorem 2. Under assumption (a1), (a1’) or (a2) we have that
• For every R > 0, q ≥ 1 ¯ sup E[exp(q|Ξ(η)|∞;[0,t])] < ∞. (7) ||η||α−H¨ol • For P0 − a.e. ω ¯ P[Ss(ω, ·) = Ξ(Y (ω))s(·), s ≤ t] = 1. (8) We split the proof in two cases. The first case uses a flow transformation, similar to what is done in Section 3 for BSDEs with rough drift, and works for assumption (a1) or (a1’). In the second case, under assumption (a2), we create a joint rough path lift of the Brownian motion B¯ and the deterministic rough path η, and then solve a classical RDE with this driving signal. The major problem here, is to show some integrability properties. 6 2. Stochastic differential equations with rough drift Proof of Case 1 Proof of Theorem 1. Let η ∈ C0,α be the lift of a smooth path η. Let Sη be the unique solution of the SDE Z t Z t Z t η η η ¯ η St = S0 + a(Sr )dr + b(Sr )dBr + c(Sr )dηr. 0 0 0 ˜η η −1 η η Define S := (φ ) (t, St ), where φ is the ODE flow Z t η η φ (t, x) = x + c(φ (r, x))dηr. (9) 0 By Lemma 3, we have that S˜η satisfies the SDE Z t Z t ˜η η ˜η ˜η ˜η ¯ St = S0 + a˜ (r, Sr )dr, + b (r, Sr )dBr (10) 0 0 witha ˜η, ˜bη defined as in Lemma 3. This equation makes sense, even if η is a generic rough path in C0,α, (in which case (9) is now really an RDE). Indeed, since the first two derivatives of φη and its inverse are bounded (Proposition 11.11 in [33]) we have thata ˜η(t, ·), ˜bη(t, ·) are also in Lip1. Hence by Theorem V.7 in [59], there exists a unique strong solution to (10). We define the mapping introduced in Theorem 1 as η ˜η Ξ(η)t := φ (t, St ). To show continuity of the mapping we restrict ourselves to the case q = 2. Moreover we shall γ+3 dS assume c1, . . . , cdY ∈ Lip (R ), and we will hence prove the local Lipschitz property of the respective maps. 1 2 0,α 1 2 Let η , η ∈ C with |η |α−H¨ol, |η |α−H¨ol < R. By Lemma 6 we have ¯ ˜1 ˜2 2 1/2 1 2 E[sup |Ss − Ss | ] ≤ CLem6(R)ρα−H¨ol(η , η ). s≤t Hence ¯ 1 2 2 1/2 ¯ 1 ˜1 2 ˜2 2 1/2 E[sup |Ξ(η )s − Ξ(η )s| ] = E[sup |φ (s, Ss ) − φ (s, Ss )| ] s≤t s≤t ¯ 1 ˜1 1 ˜2 2 1/2 ¯ 1 ˜2 2 ˜2 2 1/2 ≤ E[sup |φ (s, Ss ) − φ (s, Ss )| ] + E[sup |φ (s, Ss ) − φ (s, Ss )| ] s≤t s≤t ¯ 1 ˜1 1 ˜2 2 1/2 1 1 ≤ E[sup |φ (s, Ss ) − φ (s, Ss )| ] + sup |φ (s, x) − φ (s, x)| d s≤t s≤t,x∈R S ¯ ˜1 ˜2 2 1/2 1 2 ≤ K(R)E[sup |Ss − Ss | ] + CLem7(R)ρα−H¨ol(η , η ) s≤t 1 2 ≤ C1ρα−H¨ol(η , η ), as desired, where C1 = KLem7(R)CLem6(R) + CLem7(R), where KLem7(R) and CLem7(R) are the constants from Lemma 7. 7 2. Stochastic differential equations with rough drift Proof of Theorem 2. In order to show (7), pick k ∈ 1, . . . , dS. We first note that by simply scaling (the coefficients of) Sη it is sufficient to argue for q = 1. And consider the k-th component of Ξ. Then (k) η η ˜(k);η E[exp(|Ξ (η)|∞;[0,t])] ≤ E[exp(|Dψ |∞ |φ (0,S0)| + |S |∞;[0,t] )] η η η ˜(k);η ≤ exp(|Dψ |∞ sup |φ (0, x))E[exp(|Dψ |∞ sup |St |)] |x|≤|S0|L∞ s≤t η η ≤ exp(|Dψ |∞ sup |φ (0, x)) |x|≤|S0|L∞ η ˜(k);η η ˜(k);η × E[exp(|Dψ |∞ sup Ss )] + E[− exp(|Dψ |∞ sup Ss )] s≤t s≤t η η = exp(|Dψ |∞ sup |φ (0, x)) |x|≤|S0|L∞ η ˜(k);η η ˜(k);η × E[sup exp(|Dψ |∞Ss )] + E[− sup exp(|Dψ |∞Ss )] . s≤t s≤t Now, only the boundedness of the last two terms remains to be shown, for η bounded. By applying Itˆo’s formula we get that Z t Z t ˜(k);η ˜(k);η ˜(k);η ˜(k);η ˜(k);η exp(St ) = 1 + exp(Sr )dSr + exp(Sr )dhS ir 0 0 Z t dB Z t ˜(k);η η ˜(k);η X ˜(k);η ˜η ˜(k);η ¯i = 1 + exp(Sr )˜ak (Sr )dr + exp(Sr )bki(Sr )dBr 0 i=1 0 dB Z t X ˜(k);η ˜η ˜(k);η 2 + exp(Sr )|bki(Sr )| dr. i=1 0 Hence the process exp(S˜(k);η) satisfies a SDE with Lipschitz coefficients and by an application of Gronwall’s Lemma and Burkholder-Davis-Gundy inequality (see also Lemma V.2 in [59]) one arrives at ¯ η ˜(k);η sup sup E[exp(|Dψ |∞St )] ≤ C exp(C2), (11) |η|α−H¨ol η ˜η 2 C2 := sup |a˜ |∞ + |b |∞, η:||η||α−H¨ol ¯ η ˜(k);η which is finite because of Lemma 4. One argues analogously for sups≤t E[− exp(|Dψ |∞St )], which then gives (7). Now, for the correspondence to an SDE solution let Ω be the additional probability space as given in the statement. Let S be the solution to the SDE (6). In Section 3 in [8] it was shown (see also Theorem 2 in [1]), that if we let Θ be the stochastic (Stratonovich) flow Z t Θ(ω; t, x) = x + c(Θ(ω; r, x)) ◦ dYr(ω), 0 8 2. Stochastic differential equations with rough drift ˆ −1 ˆ then with St := Θ (t, St) we have P-a.s. Z s Z s ˆ B¯ ˆ ˆ ˆ ¯ ˆ B¯ Ss(ω, ω ) = S0 + aˆ(r, Sr)dr + b(r, Sr)dBr, s ∈ [0, t], P a.e. (ω, ω ). (12) 0 0 Here, componentwise, X 1 X X aˆ(t, x) := ∂ Θ−1(t, Θ(t, x))a (Θ(t, x)) + ∂ Θ−1(t, Θ(t, x)) b (Θ(t, x))b (Θ(t, x)), i xk i k 2 xj xk i jl kl k j,k l ˆ X −1 b(t, x)ij := ∂xk Θi (t, x)bkj(Θ(t, x)). k Especially, by a Fubini type theorem (e.g. Theorem 3.4.1 in [5]), there exists Ω0 with P0(Ω0) = ¯ 1 such that for ω ∈ Ω0 equation (12) holds true P a.s.. Let Y ∈ C0,α be the enhanced Brownian motion over Y . We can then construct ω-wise the Y(ω) ˜Y(ω) rough flow φ as given in (9). By the very definition of Ξ we know that St (ω) := Y(ω) −1 (φ ) (ω; t, Ξ(ω)t) satisfies the SDE Z t Z t ˜Y(ω) ˆY(ω) ˜Y(ω) ˆY(ω) ˜Y(ω) ¯ ¯ B¯ St = S0 + b (r, Sr )dr + b (r, Sr )dBr, Pa.e. ω , (13) 0 0 where Y(ω) X Y(ω) −1 Y(ω) Y(ω) a˜ (t, x)i := ∂xk (φ )i (t, φ (t, x))ak(t, φ (t, x)) k 1 X X + ∂ (φY(ω))−1(t, φY(ω)(t, x)) b (t, φY(ω)(t, x))b (t, φY(ω)(t, x)), 2 xj xk i jl kl j,k l ˜Y(ω) X Y(ω) −1 Y(ω) b (t, x)ij := ∂xk (φ )i (t, x)bkj(t, φ (t, x)). k It is a classical rough path result (see e.g. section 17.5 in [33]), that there exists Ω1 with Y P (Ω1) = 1 such that for ω ∈ Ω1, we have φY(ω)(·, ·) = Θ(ω; ·, ·). Y(ω) Y(ω) Hence for ω ∈ Ω1 we have thata ˆ =a ˜ , ˆb = ˜b . Hence for for ω ∈ Ω0 ∩ Ω1 the processes ˆ ˜Y(ω) ¯ 6 St(ω, ·), St (·) satisfy the same Lipschitz SDE (with respect to P). By strong uniqueness ¯ we hence have for ω ∈ Ω0 ∩ Ω1 that P a.s. ˆ ˜Y(ω) Ss(ω, ·) = Ss (·), s ≤ t. Hence for ω ∈ Ω0 ∩ Ω1 ¯ Ss(ω, ·) = Ξ(Y(ω))(·)s, s ≤ t, P a.s. 6 Here one has to argue, that fixing ω in equation (13) gives (P0-a.s.) the solution to the respective SDE on ΩB¯ . This is similar to the arguments given in Case 2 below, so we omit them. 9 2. Stochastic differential equations with rough drift Lemma 3. Let η be a smooth dY -dim path η and S be the solution of the the following classical SDE Z t Z t Z t St = S0 + a(Sr)dr + b(Sr)dB¯r + c(Sr)dηr 0 0 0 ¯ R t PdS R t i where B is a dB-dimensional Brownian motion, 0 c(Sr)dηr := i=1 0 ci(Sr)η ˙rdr, a ∈ 1 dS 1 dS 3 dS ∞ ¯ dS Lip (R ), b1, . . . , bdB ∈ Lip (R ), c1, . . . , cdY ∈ Lip (R ), and S0 ∈ L (Ω; R ) inde- pendent of B¯. Consider the flow Z t φ(t, x) = x + c(φ(r, x))dηr. (14) 0 −1 Then S˜t := φ (t, St) satisfies the following SDE Z t Z t S˜t = S0 + a˜(r, S˜r)dr + ˜b(r, S˜r)dB¯r, 0 0 where we define componentwise X 1 X X a˜(t, x) := ∂ φ−1(t, φ(t, x))a (φ(t, x)) + ∂ φ−1(t, φ(t, x)) b (φ(t, x))b (φ(t, x)), i xk i k 2 xj xk i jl kl k j,k l ˜ X −1 b(t, x)ij := ∂xk φi (t, φ(t, x))bkj(φ(t, x)). k Proof. Denote ψ(t, x) := φ−1(t, x). Then Z t ψ(r, x) = x − ∂xψ(r, x)c(x)dηr. 0 By Itˆo’sformula Z t X Z t ψi(t, St) − ψi(0,S0) = ∂tψi(r, Sr)dr + ∂xj ψi(r, Sr)dSj(r) 0 j 0 X 1 Z t + ∂ ψ (r, S )dhS ,S i 2 xj xk i r k j r j,k 0 Z t Z t X X X ¯ = ∂xj ψi(r, Sr)aj(Sr)dr + ∂xj ψi(r, Sr) bjk(Sr)dBk(r) j 0 j 0 k X 1 Z t X + ∂ ψ (S ) b (S )b (S )dr. 2 xj xk i r kl r jl r j,k 0 l Lemma 4. Consider for a rough path η ∈ C0,α the analogously to Lemma 3 transformed coefficients, a˜η, ˜bη, i.e. consider the rough flow Z t η φ(t, x) = φ (t, x) = x + c(φ(r, x))dηr. (15) 0 10 2. Stochastic differential equations with rough drift and define X 1 X X a˜η(t, x) := ∂ φ−1(t, φ(t, x))a (φ(t, x)) + ∂ φ−1(t, φ(t, x)) b (φ(t, x))b (φ(t, x)), i xk i k 2 xj xk i jl kl k j,k l ˜η X −1 b (t, x)ij := ∂xk φi (t, φ(t, x))bkj(φ(t, x)). k Then for every R > 0 there exists KLem4 = KLem4(R) < ∞ such that η sup |a˜ |∞ ≤ KLem4, η:|η|α−H¨ol 1 2 1 2 sup |a˜ (t, x) − a˜ (t, x)| ≤ KLem4(R)ρα−H¨ol(η , η ), t,x ˜1 ˜2 1 2 sup |b (t, x) − b (t, x)| ≤ KLem4(R)ρα−H¨ol(η , η ). t,x Proof. This is a straightforward calculation using Lemma 7 and the properties of a, b. The following is a standard result for continuous dependence of SDEs on parameters. Lemma 5. Let a˜i(t, x), ˜bi(t, x), i = 1, 2 be bounded and uniformly Lipschitz in x. Let S˜i be the corresponding unique solutions to the SDEs Z t Z t ˜i i ˜i ˜i ˜i ¯ St = S0 + a˜ (r, Sr)dr + b (r, Sr)dBr, i = 1, 2. 0 0 Assume 1 1 sup |Da˜ (s, ·)|∞, sup |D˜b (s, ·)|∞ < K < ∞, s≤t s≤t sup a˜1(r, x) − a˜2(r, x) , sup ˜b1(r, x) − ˜b2(r, x) < ε < ∞. r,x r,x Then there exists CLem5 = CLem5(K) such that ˜1 ˜2 2 1/2 E[sup |Ss − Ss | ] ≤ CLem5 ε, s≤t Proof. This is a straightforward application of Ito’s formula and Burkholder-Davis-Gundy inequality. We now apply the previous Lemma to our concrete setting. 11 2. Stochastic differential equations with rough drift Lemma 6. Let η1, η2 ∈ C0,α and let S˜1, S˜2 be the corresponding unique solutions to the SDEs Z t Z t ˜i i ˜i ˜i ˜i ¯ St = S0 + a˜ (r, Sr)dr + b (r, Sr)dBr, i = 1, 2, 0 0 where a˜i, ˜bi are given as in Lemma 4. 1 2 Assume R > max{|η |α−H¨ol, |η |α−H¨ol}. Then there exists CLem6 = CLem6(R). ˜1 ˜2 2 1/2 1 2 E[sup |Ss − Ss | ] ≤ CLem6ρα−H¨ol(η , η ). s≤t 1 2 Proof. Fix R > 0. Let ||η ||α−H¨ol, ||η ||α−H¨ol < R. From Lemma 4 we know that ˜1 ˜2 1 2 sup |b (t, x) − b (t, x)| ≤ KLem4(R)ρα−H¨ol(η , η ), t,x Analogously, we get 1 2 1 2 sup |a˜ (t, x) − a˜ (t, x)| ≤ L2ρα−H¨ol(η , η ), t,x for a L2 = L2(R). 1 Lemma 7. Let α ∈ (0, 1). Let γ > α ≥ 1, k ∈ {1, 2,... } and assume that V = (V1,...,Vd) γ+k e e is a collection of Lip -vector fields on R . Write n = (n1, . . . , ne) ∈ N and assume |n| := n1 + ··· + ne ≤ k. Then, for all R > 0 there exist C = C(R, |V |Lipγ+k ),K = K(R, |V |Lipγ+k )) such that if 1 2 α−H¨ol [p] d i x , x ∈ C ([0, t],G (R )) with maxi ||x ||α−H¨ol;[0,t] ≤ R then 1 2 1 2 sup ∂nπ(V )(0, y0; x ) − ∂nπ(V )(0, y0; x ) ≤ Cρα−H¨ol(x , x ), e α−H¨ol;[0,t] y0∈R 1 −1 2 −1 1 2 sup ∂nπ(V )(0, y0; x ) − ∂nπ(V )(0, y0; x ) ≤ Cρα−H¨ol(x , x ), e α−H¨ol;[0,t] y0∈R 1 sup ∂nπ(V )(0, y0; x ) ≤ K, e α−H¨ol;[0,t] y0∈R 1 −1 sup ∂nπ(V )(0, y0; x ) ≤ K. e α−H¨ol;[0,t] y0∈R Proof. The fact that V ∈ Lipγ+k (instead of just Lipγ+k−1) entails that the derivatives up to γ order k are unique non-explosive solutions to RDEs with Liploc vector fields (see Section 11 in [33]). Localization (uniform for driving paths bounded in α-H¨oldernorm) then yields the desired results. 12 2. Stochastic differential equations with rough drift Proof of Case 2 Proof of Theorem 1. We will first define the solution mapping and only in the end we will see that it is the limit of classical SDEs. B¯ Write η = π1 (η) for the projection of η to the first level. We define a path λ = λ η, ω 2 d taking values in G (R ), d := dY + dB. The first level is given as ¯ λ(1) = η, B¯ ωB (2);i,j (2);i,j (2);i,j B¯ The elements λ of the second level are defined as η for i, j ∈ {1, . . . , d1}, B¯ ω for i, j ∈ {d1 + 1, . . . , d1 + d2} and the remaining second-level entries as Z t (2);i,j i ¯j λt = ηudBu , for i ∈ {1, . . . , d1} , j ∈ {1 + d1, . . . , d1 + d2} , 0 Z t (2);i,j j ¯i j ¯i λ = ηt Bt − ηudBu , otherwise. 0 −1 The Campbell-Baker-Hausdorff formula (Theorem 7.24 in [33]) then shows that λs,t := λs ⊗ λt = exp ηs,t, B¯s,t + As,t where 1 R t i j R t j i η dηu − ηs,udη , for i, j ∈ (1, ··· , dY ) 2 s s,u s u t j−d t j−d 1 R ¯i−dY ¯ Y R ¯ Y ¯i−dY i,j 2 s Bs,u dBu − s Bs,u dBu , for i, j ∈ (dY + 1, ··· , dY + dB) so (d) 3 As,t := R t i ¯j−dY 1 i ¯j−dY η dBu − η B , for i ∈ (1, ··· , dY ), j ∈ (dY + 1, ··· , dY + dB) s s,u 2 s,t s,t t j j i−d R ¯i−dY 1 ¯ Y − 0 ηudBu + 2 ηs,tBs,t , for i ∈ (dY + 1, ··· , dY + dB), j ∈ (1, ··· , dY ) 2 d B¯ 0,α0−H¨ol 2 d That is, λt takes values in G (R ). We demonstrate that λ η, ω ∈ C G (R ) 0 (P-a.s.) for α < α by showing that |λ|α0 even has Gaussian tails. It is clear that the first level of λ is α-H¨older; to deal with the second level note that P ij |Ast| ∼ ij Ast and that 2α−H¨olderregularity already holds for the first two cases, that is i, j ∈ {1, . . . , d1} and i, j ∈ {d1 + 1, . . . , d1 + d2} (in fact even a Gauss tail via the Fernique estimate for rough path norms). Now for the remaining case we have by definition that Z t 1 ij i j i j As,t ≤ ηs,udBu + ηs,tBs,t s 2 j j and since ηi B ≤ |η| B |t − s|2α it just remains to treat R t ηi dB . But s,t s,t α−H¨ol α−H¨ol s s,u u q R t i j R t i 2 since for each s < t, s ηs,udBu has the same distribution as s ηs,u du.Z for some fixed Z ∼ N (0, 1), we also have r 2 2 t 2 R t i j R i η dB ηs,u du s s,u u s exp η =Law exp ηZ2 2α 2α (t − s) (t − s) and by the elementary estimate R t ηi 2 du ≤ |η| (t − s)2α+1 we can conclude by s s,u α−H¨ol taking sups 13 2. Stochastic differential equations with rough drift for some κ > 0. By Theorem A.19 in [33] this already yields the desired 2α0-H¨olderregularity R t i j 0 of η dB for any α < α. Putting everything together, we have shown that kλk 0 < s s,u u α −H¨ol ∞. We can actually deduce for every R > 0 the existence of b = b(α0, α, T, R) > 0 such that 2 sup E exp(b|λ|α0−H¨ol) < ∞, (16) |η|α≤R α0−H¨ol 2 d 0 So indeed λ ∈ C G (R ) , P − a.s. for all α < α. Hence (see e.g. Corollary 8.24 in 0,α0−H¨ol 2 d 0 [33]) λ ∈ C G (R ) , P − a.s. for all α < α. Using α0 < α large enough, such that γ > 1/α0, we can then use standard existence and uniqueness for the solution of the RDE Z t Z t λ λ λ St = S0 + ˚a(Sr )dr + V (Sr )dλr, (17) 0 0 k k PdB k 7 where V = (c, b) and ˚a (x) := a (x) + i=1hDbi (x), bi(x)i. We then define the mapping λ as Ξ(η) := St . We proceed to show continuity of the mapping of a rough path η to the lift λ(η). Fix α0 < α 0 0 such that γ > 1/α . We shall use Theorem A.13 in [33]. Let then q ≥ q0(α, α ), as given 0,α in this Theorem. We show that for η, η¯ ∈ C with |η|α−H¨ol, |η¯|α−H¨ol ≤ R we have for the corresponding lifts λ, λ¯ that ¯ q 1/q E[ρα0 (λ, λ) ] ≤ CLipρα(η, η¯), 0 where CLip = CLip(q, R, α, α ). ¯ Denote ε := ρα0 (λ, λ). By (16) exists a constant C1 = C1(q, R) such that q ¯ ¯ q αq E[|d(λs, λt)| ], E[|d(λs, λt)| ] ≤ C1|t − s| . Moreover ¯ q q E[|π1(λs,t − λs,t)| ] = |ηs,t − η¯s,t| ≤ εq|t − s|αq. and (again the constants C may only depend on R and q) 1 [|π (λ − λ¯ )|q/2] = [| π (λ − λ¯ ) ⊗ π (λ − λ¯ ) + A − A¯ |q/2] E 2 s,t s,t E 2 1 s,t s,t 1 s,t s,t s,t s,t q/2 αq ¯ q/2 ≤ Cε |t − s| + CE[|As,t − As,t| ] Also, since n ≥ N Z t Z t Z t Z t i j j i i j j i q/2 q/2 αq | ηs,udηu − ηs,udηu − η¯s,udη¯u − η¯s,udη¯u | ≤ ε |t − s| . s s s s 7This is a plus the Stratonovich corrector for the dB¯ integral. Note that we only have ˚a ∈ Lip1 so we have to use results on RDEs with drift (e.g. Theorem 12.6 and Theorem 12.10 in [33]) to get existence of a unique solution. 14 2. Stochastic differential equations with rough drift And finally Z t 1 Z t 1 i ¯j−dY i ¯j−dY i ¯j−dY i ¯j−dY q/2 E[| ηs,udBu − ηs,tBs,t − η¯s,udBu − η¯s,tBs,t | ] s 2 s 2 Z t i i ¯j−dY q/2 i ¯j−dY i ¯j−dY q/2 ≤ CE[| ηs,u − η¯s,udBu | ] + CE[|ηs,tBs,t − η¯s,tBs,t | ] s Z t i i 2 q/4 i i q/2 ¯j−dY q/2 ≤ C| |ηs,u − η¯s,u| du| + C|ηs,t − ηs,t| E[|Bs,t | ] s ≤ Cεq/2(t − s)αq/2+q/4 + Cεq/2(t − s)αq/2(t − s)q/4 ≤ Cεq/2(t − s)αq. Hence ¯ q/2 q/2 αq E[|π2(λs,t − λs,t)| ] ≤ C2ε |t − s| 1/q 1/(2q) So with M = M(R, q) := max{1,C1 ,C2 } we have for all q ≥ 1 ¯ q 1/q α E[|d(λs, λt)| ] ≤ M|t − s| , ¯ q 1/q α E[|π1(λs,t − λs,t)| ] ≤ εM|t − s| , ¯ q/2 2/q 2 2α E[|π2(λs,t − λs,t)| ] ≤ εM |t − s| . Hence by Theorem A.13 (i) in [33] there exists q large enough and a K = K(α, α0, T, q) such that ¯ q |π1(λs,t − λs,t| 1/q E[ sup α0 ] ≤ εKM, s ¯ !q ¯ !q q 1/q π1(λs,t − λs,t) 1/q π2(λs,t − λs,t) 1/q 0 ¯ E[ρα −H¨ol(λ, λ) ] ≤ E[ sup α0 ] + E[ sup 2α0 ] ≤ CLipε, s6=t |t − s| s6=t |t − s| as desired. This consideration was for q large enough, but since Lr is Lipschitz continuously embedded in Lq for r > q, the result follows for all q. We can now proceed to show local Lipschitzness of the solution mapping Ξ. Let still η, η¯ ∈ C0,α, Let α0 < α with γ > 1/α0 and let p0 := 1/α0. By Theorem 8 in [13] we have the deterministic estimate 8, that for some C∗, any β ∈ (0, 1) there exists C = C(C∗, β) such that ¯ ∗ ρp0−var(Ξ(η), Ξ(η¯)) ≤ Cρp0−var(λ, λ) exp(−N(β, C , υ) log(1 − β)), where υ, N are given as in Lemma 8 below. 8This estimate modifies the estimates in [33] in a way that is essential to get integrability in our context. 15 2. Stochastic differential equations with rough drift Hence, for every q ≥ 1, q 1/q q 1/q E[sup |Ξ(η)t − Ξ(η¯)t| ] ≤ E[ρp0−var(Ξ(η), Ξ(η¯)) ] t≤T ¯ 2q 1/(2q) ∗ 2q 1/(2q) ≤ CE[ ρp0−var(λ, λ) ] E[[exp(−N(β, C , υ) log(1 − β))] ] ¯ 2q 1/(2q) ∗ 2q 1/(2q) ≤ CE[ ρα0−H¨ol(λ, λ) ] E[[exp(−N(β, C , υ) log(1 − β))] ] ∗ 2q 1/(2q) ≤ Cρα0−H¨ol(η, η¯)E[[exp(−N(β, C , υ) log(1 − β))] ] . By Lemma 8 we know that for every q ≥ 1, R > 0 there exists β0(q, R) > 0 such that for |η|α−H¨ol, |η¯|α−H¨ol < R, we have for β ≤ β0 that ∗ q 1/q E[[exp(−N(β, C , υ) log(1 − β))] ] is bounded. This yields the desired local Lipschitzness of Ξ. Proof of Theorem 2. Fix R > 0. We show ¯ (k) sup E[exp(qΞ (η)t)] < ∞. (18) ||η||α−H¨ol By Theorem 10.14 in [33] we have p0 |S|p0−var;[s,t] ≤ C ||λ||p0−var;[s,t] ∨ ||λ||p0−var;[s,t] . Define the control p0 p0 υ(s, t) := ||λ||p0−var;[s,t] ∨ ||λ||p0−var;[s,t] . Then for any partition D of [0, t] X ||S||∞;[0,t] ≤ |S0| + ||S − Sti ||∞;[ti,ti+1] ti∈D X ≤ |S0| + ||S||p−var;[ti,ti+1] ti∈D X 1/p0 ≤ |S0| + C υ(ti, ti+1) . ti∈D Then, using for a > 0 the partition t0(a) := 0, p0 ti+1(a) := inf{s ≥ ti(a): υ(σi(a), s) ≥ a } ∨ t, we get p0 ||S||∞;[0,t] ≤ |S0| + Ca (sup{n : tn < t} + 1) . The exponential integrability (18) then follows as in the proof of Lemma 8. We now show correspondence to the SDE solution. On Ωˆ let EBM(Y, B¯) be the enhanced Brownian motion lift of (Y, B¯). We want to compare it to the mapping λEBM(Y ), where we 16 2. Stochastic differential equations with rough drift plug the enhanced Brownian motion lift of Y into the above constructed mapping. Note, that this composition is in general not a random variable on Ω,ˆ since joint measurability is not clear. We know by definition that π1(λ(EBM(Y ))) is measurable and coincides with π1(EBM(Y, B¯)) Pˆ-a.s. We now consider only one component of the second level, the others follow similarly. So without loss of generality let i ∈ {1, . . . , dY }, j ∈ {dY + 1, . . . , dY + dB}. Then Z t ¯ (2);i,j i ¯j−dY ˆ EBM(Y, B)t = Yr dBr P − a.s. 0 Now by results on stochastic integrals we have n t 2 Z ¯ h ¯ ¯ i i ¯j−dY B X k j−dY B j−dY B Y dB (ω, ω ) = lim Y n (ω) B n (ω ) − B n (ω ) r r n→∞ (k−1)/2 t k/2 t (k−1)/2 t 0 k=1 on Aˆ ⊂ Ωˆ with Pˆ[Aˆ] = 1. By a Fubini type theorem (e.g. Theorem 3.4.1 in [5]), there exists ˚ ˚ ˆ B¯ B¯ ˆ B¯ ˆ a Ω ⊂ Ω such that for ω ∈ Ω we have that Aω := {ω :(ω, ω ) ∈ A} satisfy P [Aω] = 1. On the other hand for fixed ω ∈ Ω Z t (2);i,j i ¯j−d1 λt (EBM(Y )(ω)) = Yr (ω)dBr . 0 Now n t 2 m Z ¯ h ¯ ¯ i i ¯j−dY B X k j−dY B j−dY B Y (ω)dB (ω ) = lim Y nm (ω) B n (ω ) − B n (ω ) r r m→∞ (k−1)/2 t k/2 m t (k−1)/2 m t 0 k=1 B Bˆ Bˆ for ω ∈ Dω ⊂ Ω with P [Dω] = 1. Here, Dω as well as the subsequence (nm)m will depend on ω. B So for ω ∈ ˚Ω, ω ∈ Aˆω ∩ Dω we have that n t 2 Z ¯ h ¯ ¯ i i ¯j−dY B X k j−dY B j−dY B Y dB (ω, ω ) = lim Y n (ω) B n (ω ) − B n (ω ) r r n→∞ (k−1)/2 t k/2 t (k−1)/2 t 0 k=1 2nm X k h j−dY B¯ j−dY B¯ i = lim Y nm (ω) B n (ω ) − B n (ω ) m→∞ (k−1)/2 t k/2 m t (k−1)/2 m t k=1 t Z ¯ i ¯j−dY B = Yr (ω)dBr (ω ). 0 Here, the second equality follows from the fact that the sum already converges along n, so it will converge along any subsequence (nm)m. ˚ B¯ ˆ Y Noting that for ω ∈ Ω we have P [Aω ∩ Dω] = 1 we get that for P -a.e. ω ∈ Ω ¯ B¯ EBM(Y, B)· = λ·(EBM(Y )(ω)), P − a.s. (19) Now by classical rough paths results (see e.g. Section 17.5 in [33]), the RDE solution to (17) driven by EBM(Y, B¯) corresponds Pˆ-a.s. to the SDE solution of (6). Combining with (19) we arrive at (8). 17 2. Stochastic differential equations with rough drift The results in Section 17.5 in [33] also yield, that if η is (the lift of) a smooth path η, then Ξ(η) actually solves (5) as a classical SDE. Together with continuity of the mapping Ξ this yields the first claim of Theorem 1. Lemma 8. Let η, η¯ ∈ C0,α, with corresponding lifts (as in the proof of Case 2 above) λ, λ¯. 0 0 0 0 Let |η|α, |η¯|α < R and Let C > 0. Let α < α with γ > 1/α and let p := 1/α . Define p0 ¯ p0 υ(s, t) := ||λ||p0−var;[s,t] + ||λ||p0−var;[s,t], & ' X N(x, y, υ) := f(x, y) sup υ(ti, ti+1) , −1 υ(ti,ti+1)≤f(x,y) i where yp0 f(x, y) := , ψ−1(x)p ψ(x) := x exp(x). For every q ≥ 1 there exists β0 = β0(R, C, q) > 0,K = K(R, C, q) > 0 such that for β < β0 || exp(−N(β, C, υ) log(1 − β))||Lq ≤ K. Proof. Define −1 −1 M(β, C, υ) := sup{n : τ1(f(β, C) ) + ··· + τn(f(β, C) ) ≤ T }, where τ0(a) := 0, τi+1(a) := inf{t ≥ 0 : υ(σi(a), σi(a) + t) ≥ a}, with the convention inf ∅ = ∞ and ∞ − ∞ = ∞, where i X σi(a) := τj(a). j=1 By the argument in Proposition 4.6 in [12] one sees that N(β, C, υ) ≤ 2M(β, C, υ) + 1, so we can focus on M. First −1 −1 P[M(β, C, υ) > n] = P[τ1(f(β, C) ) + ··· + τn(f(β, C) ) ≤ T ] −1 −1 ≤ P[τ1(f(β, C) ) ∧ φ(β) + ··· + τn(f(β, C) ) ∧ φ(β) ≤ T ], where we take φ to be the function given in Lemma 9. Define −1 mi := mi(β) := E[τi(f(β, C) ) ∧ φ(β)], n 1 X S¯ := τ (f(β, C)−1) ∧ φ(β) − m . n n i i i=1 18 2. Stochastic differential equations with rough drift then n T 1 X [M(β, C, υ) > n] ≤ [S¯ ≤ − m ]. P P n n n i i=1 1 Pn Now Lemma 9 shows, that 0 < m(β) ≤ n i=1 mi for all n, where m(β) := κφ(β). Then, if n ≥ n0(T, β) := 2T/m(β), 1 · · · ≤ [S¯ ≤ − m(β)]. P n 2 1 If ε0 = 2 m(β), then ¯ ... ≤ P[Sn ≤ −ε0] 2ε2n2 ≤ 2 exp(− 0 ). nφ(β)2 by Hoeffding’s inequality (see e.g. Theorem 2 in [41]). Now ∞ X log n [exp(−q log(1 − β)M(β, C, υ))] ≤ [M(β, C, υ) ≥ b c] E P q(− log(1 − β)) n=1 Now for n ≥ n1 := exp((n0 + 1)(− log(1 − β))q), one has by the above, log n 2ε2 log n [M(β, C, υ) ≥ b c] ≤ 2 exp(− 0 b c) P q(− log(1 − β)) φ(β)2 q(− log(1 − β)) 2ε2 2ε2 log n ≤ 2 exp( 0 ) exp(− 0 ) φ(β)2 φ(β)2 q(− log(1 − β)) 2 2ε2 2ε − 0 = 2 exp( 0 )n φ(β)2q(− log(1−β)) φ(β)2 1 κ 1 − 2 = 2 exp( κ)n q(− log(1−β)) 2 Finally, pick β0 = β0(q) such that 1 κ 2 ≥ 2, q(− log(1 − β)) for β ≤ β0. Then for β ≤ β0 1 [exp(−q log(1 − β)M(β, C, υ))] ≤ n (β) + 2 exp( κ)ζ(2), E 1 2 and n1(β) is bounded for |η|α−H¨ol, |η¯|α−H¨ol bounded, by Lemma 9. This finishes the proof. 19 2. Stochastic differential equations with rough drift Lemma 9. In the setting of the previous Theorem, there exists φ = φη,η¯ : R+ → R+, such that φ(x) > 0 for x 6= 0, and κ > 0 such that for all n ≥ 1 1 X [τ (f(β, C)−1) ∧ φ(β)] ≥ κφ(β), n E i i=1 and such that for x > 0, φ(x) is bounded away from 0 for |η|α−H¨ol, |η¯|α−H¨ol bounded. Proof. First for every i −1 −1 E[τi(f(β, C) ) ∧ φ(β)] ≥ P[τi(f(β, C) ) ≥ φ(β)] · φ(β) −1 = 1 − P[τi(f(β, C) ) < φ(β)] · φ(β) −1 ≥ 1 − P[τi(f(β, C) ) ≤ φ(β)] · φ(β) −1 −1 −1 = 1 − P[υ(σi(f(β, C) ), σi(f(β, C) ) + φ(β)) ≥ f(β, C) ] · φ(β). Now, by Markov inequality, −1 −1 −1 P[υ(σi(f(β, C) ), σi(f(β, C) ) + φ(β)) ≥ f(β, C) ] −1 −1 ≤ f(β, C)E[υ(σi(f(β, C) ), σi(f(β, C) ) + φ(β))]. Hence −1 −1 −1 E[τi(f(β, C) ) ∧ φ(β)] ≥ 1 − f(β, C)E[υ(σi(f(β, C) ), σi(f(β, C) ) + φ(β))] · φ(β). Summing up we get 1 X [τ (f(β, C)−1) ∧ φ(β)] n E i i=1 n 1 X ≥ 1 − f(β, C) [υ(σ (f(β, C)−1), σ (f(β, C)−1) + φ(β))] · φ(β) n E i i (20) i=1 n 1 X = 1 − f(β, C) [υ(σ (f(β, C)−1), σ (f(β, C)−1) + φ(β))] · φ(β). n E i i i=1 Now n 1 X 1 − f(β, C) [υ(σ (f(β, C)−1), σ (f(β, C)−1) + φ(β))] n E i i i=1 ≥ 1 − f(β, C) sup E[υ(τ, τ + φ(β))] , τ where the supremum is over all stopping times τ. We define 1 1 φ(β) := sup{a : sup E[υ(τ, τ + a)] ≤ }. 2 τ 2f(β, C) We claim that φ > 0 for β > 0. 20 2. Stochastic differential equations with rough drift Indeed, note that p0 ¯ p0 sup E[υ(τ, τ + a)] = sup E[||λ||p0−var;[τ,τ+a] + ||λ||p0−var;[τ,τ+a]] τ τ p0 ¯ p0 ≤ sup E[a||λ||α0−H¨ol;[0,T ] + a||λ||α0−H¨ol;[0,T ]] τ p0 ¯ p0 = E[||λ||α0−H¨ol;[0,T ] + ||λ||α0−H¨ol;[0,T ]]a, which shows that φ(β) > 0 for β > 0. By (16) we have in particular that φ(β) is bounded away from 0 for |η|α−H¨ol, |η¯|α−H¨ol bounded. Then by (20) 1 X −1 E[τi(f(β, C) ) ∧ φ(β)] ≥ 1 − f(β, C) sup E[υ(τ, τ + φ(β))] · φ(β) n τ i=1 1 ≥ φ(β), 2 so κ := 1/2 does the job. 2.3 The Young case In Section 2.2 we showed existence, uniqueness and stability of an SDE with rough drift under the assumption that the vector fields in front of the Brownian motion are Lip1 and the vector fields in front of the rough path are Lipγ+2, if the rough path is α-H¨olderand γ > 1/α. The aim of this section is to show that in the case of α > 1/2 actually Lipγ regularity suffices in front of the rough path, which is the regularity one would need when solving a standard RDE. Moreover the vector fields in front of the Brownian motion do not have to be bounded anymore, as one would hope when thinking of the standard assumption for classical SDEs. The results are most conveniently derived when working in variation spaces. In this section [q] d d d q ∈ (1, 2), so of course G (R Y ) = R Y and we are dealing with classical paths in R Y . Let γq > q. Assume γq dS (a3) a Lipschitz, b1, . . . , bdB Lipschitz and c1, . . . , cdY ∈ Lip (R ) ¯ ¯ ¯ ¯ Let as before (Ω, F, (Ft)t≥0, P) be a filtered probability space carrying a dB-dimensional Brownian motion B¯ and a bounded dS-dimensional random vector S0 independent of B¯. Let n d n 0,q−var η : [0, t] → R Y be smooth paths, such that η → η in q-variation for some η ∈ C and n let S be a dS-dimensional process which is the unique solution to the classical SDE Z t Z t Z t n n n ¯ n n St = S0 + a(Sr )dr + b(Sr )dBr + c(Sr )dηr , 0 0 0 ∞ 0 Theorem 10. Under assumption (a3) there exists a dS-dimensional process S ∈ S such that Sn → S∞, in S0. In addition, the limit Ξ(η) := S∞ only depends on η and not on the approximating sequence. 21 2. Stochastic differential equations with rough drift Moreover, for all r ≥ 1, η ∈ C0,q−var it holds that Ξ(η) ∈ Sr and the corresponding mapping Ξ: C0,q−var → Sr is locally uniformly continuous. 9 Moreover, S∞ ∈ C0,p−var, for all p > 2 and satisfies the integral equation Z t Z t Z t ¯ St = S0 + a(Sr)dr + b(Sr)dBr + c(Sr)dηr, (21) 0 0 0 where the last integral is a (pathwise) Young integral. It is the only process that satisfies this equation. We prove this in four steps. First we show uniqueness of a solution to the integral equation, if it exists (Theorem 12). Then we show weak existence of a solution under the assumption c ∈ Lip1. (Theorem 13). Using standard methods from strong existence for standard SDEs, we can then deduce from these two results the existence of a unique strong solution (Theorem 14). We finally show stability of the solution with respect to the Young path η, which yields also the identification as a limit, as in the statement of the theorem. For notational convenience we assume a ≡ 0. Remark 11. The main reason why we were not able to employ a contraction mapping argu- ment directly, is the fact that the deterministic Young inequality only gives the locally Lipschitz estimate Z Z ¯ || H(Y )dη − H(Y )dη||q−var;[0,T ] ¯ ¯ ≤ C ||Y ||p−var;[0,T ] + ||Y ||p−var;[0,T ] ||Y − Y ||p−var;[0,T ]||η||q−var;[0,T ]. If we apply the expectation operator to both sides of this inequality it becomes evident that a naive Gronwall argument will fail in a stochastic setting. Theorem 12 (Pathwise Uniqueness). Let T > 0. Let q ∈ (1, 2) and η ∈ C0,q−var. γq Let γq > q and c = (c1, . . . , cdY ) be a collection of Lip vector fields, let b = (b1, . . . , bdB ) be a collection of Lipschitz (not necessarily bounded) vector fields. Let ξ be a random variable independent of B. Then there exists at most one solution S on [0,T ] of (21). Proof. Let p ∈ (2, 3) such that 1/p + 1/q > 1. Let S, S¯ be two solutions of (21) on [0,T ] (if the solutions are only valid up to a stopping time, the following argument shows that they coincide up to this stopping time). For T > 0, n ≥ 0 define the stopping times T τ0 := 0, T T T τ := inf{t ≥ τ : ||S|| T ≥ 1 or ||S¯|| T ≥ 1} ∧ T ∧ (τ + T). n+1 n p−var;[τn ,t] p−var;[τn ,t] n Since S, S¯ solve (21), their p-variation is a.s. finite (this follows e.g. from Theorem 14.9 in [33]). Hence for every T > 0, the stopping times exhaust the interval [0,T ]. 9 Actually we show that the mapping is continuous into the space of all processes X with norm ||X||p,r = 1/r E[||X||p−var;[0,T ]r ] < ∞ for all p > 2, r ≥ 1. 22 2. Stochastic differential equations with rough drift γq Assume that S T = S¯ T . First (here it is essential that c ∈ Lip ) τn τn Z Z || c(S)dη − c(S¯)dη|| T T p−var;[τn ,τn+1] Z Z ≤ || c(S)dη − c(S¯)dη|| T T q−var;[τn ,τn+1] ≤ C 1 + max{||S|| T T , ||S¯|| T T } ||η|| T T ||S − S¯|| T T p−var;[τn ,τn+1] p−var;[τn ,τn+1] q−var;[τn ,τn+1] p−var;[τn ,τn+1] ¯ ≤ C2 sup ||η||q−var;[t,t+T] ||S − S||p−var;[τ T,τ T ]. t≤T n n+1 Choosing T small enough we get Z · Z · ¯ r 1/r ¯ r 1/r E[|| c(S)dη − c(S)dη||p−var;[τ T,τ T ]] ≤ 1/3E[||S − S||p−var;[τ T,τ T ]] . 0 0 n n+1 n n+1 On the other hand, define the martingale B˜ := 1 T (t) B − B . Using the BDG t [τn ,∞) t∧τn+1 τn inequality for p-variation (Theorem 14.12 in [33]), we get Z Z ¯ r 1/r ¯ ˜ r 1/r E[|| b(S) − b(S)dB|| T T ] = E[|| b(S) − b(S)dB||p−var;[0,T ]] p−var;[τn ,τn+1] Z T r/2 ¯ 2 1/r ≤ CE[ |b(Sr) − b(Sr)| dr ] 0 n !r/2 Z τi+1 ¯ 2 1/r = CE[ |b(Sr) − b(Sr)| dr ] n τi 1/2 ¯ r 1/r ≤ CT E[ sup |Sv − Sv| ] n n τi ≤v≤τi+1 1/2 ¯ r 1/r ≤ T CE[||S − S|| T T ] , p−var;[τn ,τn+1] (where we need T ≤ 1 to bound sup with p-var; also we use the fact that S T = S¯ T ). τn τn Choosing T small enough we get Z Z ¯ r 1/r ¯ r 1/r E[|| b(S)dB − b(S)dB|| T T ] ≤ 1/3E[||S − S|| T T ] . p−var;[τn ,τn+1] p−var;[τn ,τn+1] Hence ¯ r 1/r E[||S − S|| T T ] p−var;[τn ,τn+1] Z · Z · Z · Z · ¯ ¯ r 1/r = E[|| c(Sr)dηr + b(Sr)dBr − c(Sr)dηr − b(Sr)dBr||p−var;[τ T,τ T ]] 0 0 0 0 n n+1 ¯ r 1/r ≤ 2/3E[||S − S|| T T ] . p−var;[τn ,τn+1] Hence (note that the right hand side is finite, because of the stopping time) S = S¯ on T T [τn , τn+1]. This is true for all n, hence S = S¯ on [0,T ]. Theorem 13 (Weak Existence). Let T > 0. Let q ∈ (1, 2) and η ∈ C0,q−var. γq Let γq > q and c = (c1, . . . , cdY ) be a collection of Lip vector fields, let b = (b1, . . . , bdB ) be a collection of Lipschitz Lip1 (i.e. Lipschitz and bounded) vector fields. 23 2. Stochastic differential equations with rough drift n Let µ be a probability measure on R . Then there exists a probability space carrying a Brownian motion B˜ and an adapted process S˜ such that on [0,T ] Z t Z t ˜ ˜ ˜ ˜ ˜ St = S0 + b(Sr)dBr + c(Sr)dηr, (22) 0 0 and S˜0 ∼ µ. n 4 n Proof. Let c be a sequence of collections of Lip vector fields such that |c − c|Lipγq → 0. Let a Brownian motion B be given. Let ξ be a random variable, independent of B, with ξ ∼ µ. For every n there exists a solution Sn to Z t Z t n n n n St = ξ + b(Sr )dBr + c (Sr )dηr. 0 0 Indeed, by Theorem 1 there exists a solution to corresponding SDE with rough drift, which so far only formally satisfies the integral equation. But looking at the proof, the solution is the back-transformation of an SDE with transformed coefficients. Applying Theorem 19 then gives, that Sn actually satisfies the integral equation. n 0,p Let p− < p ∈ (2, 3) such that 1/p + 1/q > 1. We will show that S is tight in C . Let T > 0. First (using BDG, exploiting boundedness of b, when applying Kolmogorov for the stochastic integral)