arXiv:2103.07044v1 [math.PR] 12 Mar 2021 se[],s h ogZkiapoiaino 11 svlal o fin for valuable is (1.1) of approximation mollifier rou Wong–Zakai any with for the market However, financial so mathematics. a in [1]), pricing option (see in emerges (1.1) Integral osnteiti eea,where general, ( in exist not does euaiyta rwinmotion Brownian than regularity (1.1) function where rwinmto ihHrtparameter Hurst with motion Brownian eie ogZkiapoiainb digatr opeetthe prevent 1.1. to Theorem term a adding by approximation Wong–Zakai revised (1.2) olwn kthilsrtsterltosbtenteconvergence the between relations the illustrates sketch following ̺ ε ecnie h olwn totp tcatcitga:for integral: Itˆo-type stochastic following the consider We 2020 h euti rvdb h hoyo euaiysrcue yHair by structures regularity of theory the by proved is result The PZRNRAIAINI ASINRUHPATHS ROUGH GAUSSIAN IN RENORMALIZATION BPHZ ∗ W W ahmtc ujc Classification. Subject Mathematics ) og ahstigb osdrn t plcto orough to application its is considering finance. renormalization by mathematical the setting path that rough demonstrate a we Here mo p the singul also of rough structure. of convergence We the theory a assures the and in in equations differential role model structure. central BPHZ-renormalized a regularity plays the renormalization and of of theory expression view plain path in path rough rough a between of renormalization (BPHZ) Zimmermann Abstract. t t f Z hsfiuei u ofatoa rwinmotion Brownian fractional to due is failure This . sasadr n-iesoa rwinmto,and motion, Brownian one-dimensional standard a is 0 : T R f → ( W Let R t H,ε ecntutapoeuefrBogoliubov–Parasiuk–Hepp– for procedure a construct We , C ) dW ε ( t W = ) t ε t − H E Z = ε 0 lim [ → W T √ 1. ˙ ̺ 0 Z C 2 AAIEITO HAYAHIDE ε ε W Z 0 ( Introduction H ( ε T t 0 · ( t ) = ) T ae ta.[]dmntae h following the demonstrated [1] al. et Bayer . t f Z W ) f rmr 03,Scnay6L0 91-10. 60L20, Secondary 60L30, Primary 0 ( f W t H,ε ( ′ H 1 ̺ W | ( ̺ t t W ( H : · − ( t H,ε ∈ /ε ) t t R H,ε )] dW s (0 ) hn ehv h olwn limit: following the have we Then, . | ) → /ε H dW ) , t dt − 1 , , R / 1 W t / ε → h limit the , )gvnby given 2) 2 t dW H,ε Z 0 T s ( = . rsohsi partial stochastic ar e nteregularity the in del f ( ̺ W ahter and theory path loeetv in effective also ε > T t H ∗ W t.BPHZ ath. W h relation the ) W dW rvd a provide t H t n smooth a and 0 H H fteintegral the of t ) safractional a is divergence. ε , aiglower having t hvolatility gh and , r[] The [8]. er → W ancial 0 . t ε = 2 H. ITO and the regularity structure. → Πε R / ΠB,ε ε 0 / Π O O O
T T T f(W H,ε)dW ε ✤ / f(W H,ε)dW ε Cε / f(W H )dW 0 t t 0 t t − 0 t t We constructR a regularity structureR and a model Πε compatibleR with the integral T H,ε ε 0 f(Wt )dWt. Because model Π does not converge to any model simply , we revise it by adding some terms to obtain model ΠB,ε. This procedure Πε ΠB,ε R is called renormalization, and ΠB,ε is called a renormalized model. We7→ can show that the renormalized model ΠB,ε converges to a model Π in the norm of model . We can also show that the renormalized model ΠB,ε and the model Π T H,ε ε T C ε ′ H,ε correspond to the revised integral 0 f(Wt )dWt 0 (t)f (Wt )dt and the T H − B,ε Itˆointegral 0 f(Wt )dWt, respectively.R The renormalizedR model Π can be written explicitly as R ΠB,ε(1) = Πε(1), ΠB,ε(Ξ) = Πε(Ξ), ΠB,ε( (Ξ)ˆ n) = Πε( (Ξ)ˆ n), s s s s s I s I (1.3) ΠB,ε(Ξ (Ξ)ˆ n) = Πε(Ξ (Ξ)n) nC ε(t)Πε(Ξ (Ξ)n−1), n 1. s I s I − s I ≥ This renormalization stems from probabilistic and analytic results such as a decomposition by the Wick product and a transformation into the Skorokhod integral; see, for example, [1, Lemma 3.9] and [1, Lemma 3.11]. Meanwhile, Bruned et al. [6] demonstrated an algebraic procedure for renormalization in regularity structures, known as Bogoliubov–Parasiuk–Hepp– Zimmermann (BPHZ) renormalization, and Chandra and Hairer [7] showed that a model subjected to BPHZ renormalization converges to some model under desired conditions. Thus, BPHZ renormalization is key in the theory of regularity structures. Also, study of the relation between rough paths and regularity structures gave rise to discussion on the renormalization of rough paths, given that they are a counterpart of regularity structures ([3], [4], [5]). To supplement the current trend in the theory of regularity structures, herein we focus on BPHZ renormalization of a regularity structure corresponding to a rough path generated by Gaussian processes. We also consider an application of BPHZ renormalization to the Wong–Zakai approximation of integral (1.1). The central theorem in this article is as follows, which corresponds to Theorem 2.7.
Theorem 1.2. Let symbols Ξ1,..., Ξd correspond to Gaussian processes ξ1,...,ξd, respectively. Let Π be a model expressing the rough path generated by ξ1,...,ξd, and let Πˆ be the BPHZ-renormalized model of Π. Assume that for any a ,...,a ,b ,...,b R, d a ξ˙ (0)+ d b (0)ξ is a normal distribution with 1 d 1 d ∈ i=1 i i j=1 j j zero mean. For 1 i d and n 1, we have ≤ ≤ P ≥ P (1.4) Πˆ (1) = Π (1), Πˆ (Ξ) = Π (Ξ), Πˆ ( (Ξ )n) = Π ( (Ξ )n). z z z z z I i z I i In addition, if Ξ (Ξ )n < 0 for 1 i, j d, and n N, then we have | iI j | ≤ ≤ ∈ (1.5) Πˆ (Ξ (Ξ )n) = Π (Ξ (Ξ )n) nE[ξ˙ (0)ξ (0)]Π ( (Ξ )n−1). z iI j z iI j − i j z I j The result of Theorem 1.2 demonstrates a plain and handy expression of BPHZ renormalization under the assumption of a rough path generated by Gaussian processes, while Bruned et al. [6] provided BPHZ renormalization in a more general BPHZ RENORMALIZATION IN GAUSSIAN ROUGH PATHS 3 setting. We apply Theorem 1.2 to the Wong–Zakai approximation, leading to the following theorem that corresponds to Theorem 3.4.
Theorem 1.3. The processes ξε and ξˆε correspond to Brownian motion and fractional Brownian motion convoluted by a mollifier, respectively. Then, we have
Πˆ ε 1 = Πε (1), Πˆ εΞ = Πε (Ξ), Πˆ ε ( (Ξ)ˆ n) = Πε( (Ξ)ˆ n), z z z z z I z I (1.6) Πˆ ε Ξ (Ξ)ˆ n = Πε Ξ (Ξ)ˆ n nE[ξ˙ε(0)ξˆε(0)]ΠεΞ (Ξ)ˆ n−1, n 1. z I z I − z I ≥ Comparing (1.3) with (1.6), we realize that BPHZ renormalization is mostly similar to the counterpart given by Bayer et al. [1]. In addition, Theorem 3.5 shows that the BPHZ-renormalized model Πˆ ε converges to a naive model Π, leading to BPHZ renormalization also being effective in the Wong–Zakai approximation in integral (1.1). This article is organized as follows. In Section 2, we discuss the construction of a model on the regularity structure and define the BPHZ renormalization of that model. In particular, we provide a simple expression of the BPHZ-renormalized model in symbols corresponding to a geometric rough path, provided that the noises are Gaussian. In Section 3, we consider applying these results to a rough volatility model in a financial market, and we demonstrate the correspondence between the renormalization in [1] and BPHZ renormalization.
2. Regularity Structure of Gaussian Rough Paths In this section, we construct a regularity structure of Gaussian rough paths, which reflects a regularity structure expressing integral (1.1). We also consider BPHZ renormalization in the regularity structure and provide a plain expression of the renormalization. Let ξ¯i (1 i d) be a stochastic process on R such that it has αi-continuity ≤ ≤ ε ¯ almost surely, where αi (0, 1). Also, let ξi := ̺ε ξi, where ̺ε = ̺( /ε)/ε, ε > ∈ ε ∗ · 0 and ̺ is a smoothing mollifier. We denote ξi simply by ξi unless ε must be emphasized. Now, we construct a regularity structure that can express the integrals
T T T (2.1) ξ dξ , ξ2ξ ,..., ξndξ ,..., 1 i, j d. i j i j i j ≤ ≤ Z0 Z0 Z0 2.1. Symbols and Tree Expression. First, we introduce symbols corresponding to the integrands and integrations in (2.1), which are fundamental elements in the regularity structure. A finite set L and a map : L R are defined respectively as | · | → (2.2) L := , Ξ ; i =1, 2,...,d , {I i } (2.3) := 1, Ξ := α 1. |I| | i| i − Remark 2.1. The set L corresponds to the symbols that express kernels and noises, so we set L = , Ξ ,..., Ξ . Also, corresponds to the convolution of the {I 1 d} I kernel K with Ξ1,..., Ξd to ξ¯1,..., ξ¯d, respectively. The degree corresponds to the regularity of the noise or kernel, so we set = 1 and Ξ = α 1 where α (0, 1). |I| | i| i − i ∈ We construct a triplet (A, T, G) as follows. 4 H. ITO
We define a set S and a vector space T as • (2.4) S := 1, Ξ , n, Ξ n, (Ξ ), Ξ (Ξ )n;1 i, j d, n N , { i I iI I i iI j ≤ ≤ ∈ } T := S , h i where S is a vector space spanned by elements in S. A set Ah i R is given by • ⊂ A := τ ; τ S , {| | ∈ } where a map : S R is defined as (2.3) and | · | → τ = τ ′ , | | ′ | | Xτ where the sum is over all τ ′ L emerging in τ. For example, ∈ Ξ (Ξ ) = Ξ + + Ξ = α + α 1, | iI j | | i| |I| | j | i j − Ξ (Ξ )n = Ξ + n( + Ξ )= α + nα 1. | iI j | | i| |I| | j | i j − A group G := Γ ; h = (h ,h ,...,h ) Rd+1 is given by • { h 0 1 d ∈ } (2.5) Γ 1 = 1, Γ Ξ = Ξ , Γ (Ξ )= (Ξ )+ h 1, Γ = + h 1, h h i i hI i I i i hI I 0 (2.6) Γ (τ τ ′)=Γ τ Γ τ ′ h · h · h for τ, τ ′ S such that τ τ ′ S. ∈ · ∈ We can easily show that the triplet (A, T, G) is a regularity structure. An element in S can be viewed as a rooted tree by the relations shown below (see Appendix A for details of rooted trees).
Ξi Ξi I ←→ I ←→
Ξi Ξj (Ξ ) Ξi (Ξj ) I i ←→ I ←→ Ξi I I n−times
Ξzj }| Ξ{ j Ξ (Ξ )n ··· iI j ←→ Ξi I I
2.2. Model on the Regularity Structure. We construct a model on the regularity structure (A, T, G) to translate an abstract symbol in S into a concrete function. We define Π : T C∞(R) as → Π(1)=1, Π(Ξ )(t)= ξ˙ (t), Π( (Ξ ))(t)= ξ (t), Π( )(t)= t, i i I i i I Π(τ τ ′)= Π(τ) Π(τ ′). · · Now, we define a model on (A, T, G) by using the map Π. The pair (Π, Γ) is given by the following. BPHZ RENORMALIZATION IN GAUSSIAN ROUGH PATHS 5
For s R, a map Π : T R is defined as • ∈ s → Πs1 =1, ΠsΞi(t)= Π(Ξi)(t), Π ( )(t)= Π( )(t) Π( )(s), s I I − I Π ( (Ξ ))(t)= Π( (Ξ ))(t) Π( (Ξ ))(s), s I i I i − I i Π (τ τ ′)(t) = Π (τ)(t) Π (τ ′)(t) s · s · s for τ, τ ′ S such that τ τ ′ S, and the domain is extended by imposing linearity.∈ · ∈ For s,t R, a map Γ : T T is defined as • ∈ ts → Γts1 = 1, Γts(Ξi) = Ξi, Γ ( )= + (Π( )(t) Π( (s)))1, ts I I I − I Γ ( (Ξ )) = (Ξ ) + (Π( (Ξ )(t) Π( (Ξ ))(s))1, ts I i I i I i − I i Γ (τ τ ′)=Γ (τ) Γ (τ ′) ts · ts · ts for τ, τ ′ S such that τ τ ′ S, and the domain is extended by imposing linearity.∈ · ∈
2.3. BPHZ Renormalization. We revise the model (Π, Γ) defined in Section 2.2 along the lines of the discussion in [6]. First, we introduce negative renormalization in a regularity structure. As follows, we define a set T− that plays a central role in negative renormalization:
(2.7) S− := τ τ ; τ S 1 ,n N 1 , { 1 ⊔···⊔ n i ∈ \{ } ∈ }∪{ } (2.8) Tˆ− := S− , h i where is the disjoint union in set theory. By convention, we replace with , which is⊔ called a forest product in [6]. • ⊔
Definition 2.2. (1) Let τ Tˆ−. Let A be a subforest of τ, that is, A = 1 or ∈ A consists of edges included in τ. Then, we define R τ Tˆ− as a forest A ∈ obtained by trimming edges of A from τ. RAτ is called a contraction of τ by A. (2) We define a map ∆− : Tˆ− Tˆ− Tˆ− as → ⊗ (2.9) ∆−τ := A RAτ, τ S−, ⊂ ∈ ⊗ ∈ A τX;A S− where A τ is a subforest of τ and R τ is the contraction of τ by A, and ⊂ A n n ∆− aiτi := ai∆−τi, i=1 ! i=1 X X where a ,...,a R and τ ,...,τ S−. The restriction ∆− : T 1 n ∈ 1 n ∈ |T → Tˆ− T denotes simply ∆− unless confusion occurs. ⊗ Example 2.3. (1) Let τ = Ξ (Ξ ), A = Ξ , and B = Ξ Ξ . Then, we have 1I 2 1 1 • 2 R τ = (Ξ ), R τ = A I 2 B I because RAτ and RBτ are obtained by trimming the edge Ξ1 and Ξ1 Ξ2, respectively. • 6 H. ITO
(2) We have
∆−Ξ =1 Ξ + Ξ 1, 1 ⊗ 1 1 ⊗ ∆−( (Ξ ))=1 (Ξ )+ Ξ + Ξ + (Ξ ) 1, I 1 ⊗ I 1 I⊗ 1 1 ⊗ I I 1 ⊗ ∆−(Ξ (Ξ ))=1 Ξ (Ξ ) + Ξ (Ξ ) + Ξ Ξ + (Ξ ) Ξ 1I 1 ⊗ 1I 1 1 ⊗ I 1 1 ⊗ 1I I 1 ⊗ 1 + (Ξ ) Ξ + Ξ Ξ + Ξ Ξ + Ξ (Ξ ). I 1 ⊗ 1 1I⊗ 1 1 • 1 ⊗ I 1I 1 Definition 2.4. (1) A subspace J T− is defined as + ⊂ J := τ S−; τ = σ σ,¯ σ, σ¯ S−, σ 0 . + h{ ∈ • ∈ | |≥ }i ex (2) A quotient space T− is defined as T− := Tˆ−/J+, and p− : Tˆ− T− denotes the canonical projection. → ex ex ex (3) We define a map ∆− : T T− T as ∆− := (p− Id)∆−. → ⊗ ⊗ ex The map ∆− takes forests including a tree with negative degree, which is ill-posed in taking a limit, from a whole tree. BPHZ renormalization involves finding ill-posed symbols and adding extra terms there to prevent divergence, so ex ∆− becomes a part of BPHZ renormalization and plays a role in finding these symbols. Definition 2.5. For the random linear map Π : T C∞(R) defined above, we define the following. → − (1) A character g (Π) on Tˆ− is defined inductively as g−(Π)(τ)= E[Πτ(0)], τ T, ∈ − ′ − − ′ ′ g (Π)(τ τ )= g (Π)(τ) g (Π)(τ ), τ,τ Tˆ−. • · ∈ ex (2) A twisted negative antipode ˜− : T− Tˆ− is defined recursively as A → ex ˜− 1=1, A ex ex ex ex ex ˜− τ = ˆ − ( ˜− Id)(∆− i− τ τ 1), τ T− 1 , A −M A ⊗ − ⊗ ∈ \{ } ex ′ ′ ex where ˆ − maps (τ, τ ) Tˆ− Tˆ− to τ τ Tˆ−, and i− : T− Tˆ− is a canonicalM injective. ∈ × • ∈ → (3) A pair (Πˆ , Γ)ˆ is defined as
− ex ex (2.10) Πˆ τ = (g (Π) ˜− Π )∆− τ, z A ⊗ z (2.11) Γˆzz¯τ =Γzz¯τ. This pair is actually a model on T by [6, Theorem 6.16], and we call this the BPHZ-renormalized model.
Remark 2.6. In [6], Γˆzz¯ is given by
− ex ex Γˆ τ = (Id γ (g (Π) Id)∆− )∆ τ, zz¯ ⊗ zz¯ ⊗ + so it is not obvious that Γˆ is equal to Γ. For the proof of that equality, see the Appendix B. Now, we give a plain expression of the BPHZ-renormalized model in a branched rough path, which is the central result herein. BPHZ RENORMALIZATION IN GAUSSIAN ROUGH PATHS 7
Theorem 2.7. Assume that for any a ,...,a ,b ,...,b R, d a ξ˙ (0) + 1 d 1 d ∈ i=1 i i d b ξ (0) is a normal distribution with zero mean. Then, for any 1 i d, we j=1 j j P have ≤ ≤ P (2.12) Πˆ (Ξ ) = Π (Ξ ), Πˆ ( (Ξ )n) = Π ( (Ξ )n). z i z i z I i z I i In addition, if Ξ (Ξ )n < 0 for 1 i, j d, and n N, then we have | iI j | ≤ ≤ ∈ (2.13) Πˆ (Ξ (Ξ )n) = Π (Ξ (Ξ )n) nE[ξ˙ (0)ξ (0)]Π ( (Ξ )n−1). z iI j z iI j − i j z I j Before proving this theorem, we prepare the following lemmas.
Lemma 2.8. Let random variables X1,...,Xn be centered jointly normal variables. Then we have n
(2.14) E Xi = [Xik Xjk ] , "i=1 # Y X Yk where the sum is taken over the set of disjoint partitions ik, jk of 1,...,n . In particular, the left-hand side of the above expression is identical{ } to zero{ if n is} odd. Proof. See [9, Theorem 1.28].
Lemma 2.9. If Ξ (Ξ )n < 0 for 1 i, j d, and n 2, then we have | iI j | ≤ ≤ ≥ (2.15) n − ex n n − ex ex l n−l g (Π) ˜− (Ξ (Ξ ) )= g (Π) ˆ − ( ˜− Id)(Ξ (Ξ ) (Ξ ) ) A iI j l M A ⊗ iI j ⊗ I j Xl=0 (2.16) =0.
− ex Proof. First, we show (2.15). Because g (Π) and ˜− are multiplicative in terms of the forest product , for 1 i d we have A • ≤ ≤ − ex − ex − ex g (Π) ˜− (Ξ τ)= g (Π) ˜− (Ξ ) g (Π) ˜− (τ) A i • A i · A − − ex = g (Π)(Ξ ) g (Π) ˜− (τ) − i · A (2.17) =0.
ex ex n ex ex n (p) (p) (p) We denote ∆− i− (Ξ (Ξ ) ) by ∆− i− (Ξ (Ξ ) )= τ τ , where τ T− iI j iI j p 1 ⊗ 2 1 ∈ (p) ˆ n (p) (p) and τ2 T−. Because Ξi (Ξj ) = τ1 + τ2 P, we have the decomposition (p) ′ ∈ ′ | I | | | | | τ1 = τ1 τm such that it satisfies either of the following conditions (a) and (b): •···• ′ ′ (a) for any 1 s m, τs = Ξi or τs = Ξj ; ≤ ≤ ′ (b) we have m = 1 and τ = Ξ (Ξ )l (1 l n). 1 iI j ≤ ≤ By this result and (2.17), we have
− ex n ex ex ex ex n n g (Π) ˜− (Ξ (Ξ) )= ˆ − ( ˜− Id)(∆− i− (Ξ (Ξ ) ) Ξ (Ξ ) 1) A iI −M A ⊗ iI j − iI j ⊗ n n − ex ex l n−l = g (Π) ˆ − ( ˜− Id)(Ξ (Ξ ) (Ξ ) ), l M A ⊗ iI j ⊗ I j Xl=0 which is identical to (2.15). 8 H. ITO
We show (2.16) by induction. We begin by considering the case of n = 2. By (2.15), we have
− ex 2 − ex ex − 2 g (Π) ˜− (Ξ (Ξ ) )= 2g (Π) ˆ − ( ˜− Id)(Ξ (Ξ ) (Ξ )) g (Π)(Ξ (Ξ ) ) A iI j − M A ⊗ iI j ⊗ I j − iI j =2g−(Π)(Ξ (Ξ ) (Ξ )) g−(Π)(Ξ (Ξ )2) iI j • I j − iI j =2E[ξ˙ (0)ξ (0)]E[ξ (0)] E[ξ˙ (0)ξ (0)2]. i j j − i j 2 By Theorem 2.8, we have E[ξj (0)] = E[ξ˙i(0)ξj (0) ] = 0 so that − ex g (Π) ˜− (Ξi (Ξj )) = 0. Next, we consider (2.16) with n n0. Then, by (2.15), we haveA I ≤
n0+1 − ex n0+1 n0 +1 − ex ex l n0+1−l g (Π) ˜− Ξ (Ξ ) = g (Π) ˆ − ( ˜− Id)(Ξ (Ξ ) (Ξ ) ) A iI j − l M A ⊗ iI j ⊗ I j l=0 X g−(Π)(Ξ (Ξ )n0+1). − iI j Using the assumption by induction, we have
− ex n0+1 − − ex n0 g (Π) ˜− Ξ (Ξ ) = (n + 1)(g (Π) g (Π))( ˜− Id)(Ξ (Ξ ) (Ξ ) ) A iI j − 0 ⊗ A ⊗ iI j ⊗ I j g−(Π)(Ξ (Ξ )n0+1) − iI j = (n + 1)E[ξ˙ (0)ξ (0)]E[ξ (0)n0 ] E[ξ˙ (0)ξ (0)n0+1]. 0 i j j − i j
n0+1 n0+1 If n0 +1 is even, then we have E[ξj (0) ]= E[ξ˙i(0)ξj (0) ]=0. If n0 +1 is odd, n0+1 n0+1 then by Lemma 2.8 we have (n0 + 1)E[ξ˙i(0)ξj (0)]E[ξj (0) ]= E[ξ˙i(0)ξj (0) ]. − ex n0+1 Thus, we have g (Π) ˜− Ξ (Ξ ) = 0. A iI j