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 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 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 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

 ex Proof of Theorem 2.7. First, we show (2.12). Because ∆− Ξi = Ξi 1 + 1 Ξi and − ⊗ ⊗ g (Π)Ξi = E[ξ˙i(0)] = 0, we have

Πˆ z(Ξi) = Πz(Ξi).

ex n n We calculate ∆− ( (Ξi) ). Because there exists no tree τ (Ξi) such that − I ex ⊂ I τ < 0 and g (Π) ˜− τ = 0, we have | | A 6 n − ex ex n Πˆ ( (Ξ ) ) = (g (Π) ˜− Π )∆− ( (Ξ ) ) z I i A ⊗ z I i − ex n = (g (Π) ˜− Π )(1 (Ξ ) ) A ⊗ z ⊗ I i = Π ( (Ξ )n). z I i Next, we show (2.13). By using (2.16) in Lemma 2.9, we have

n n n − ex l n−l n Πˆ (Ξ (Ξ ) )= (g (Π) ˜− Π )(Ξ (Ξ ) (Ξ ) ) + Π (Ξ (Ξ ) ) z iI j l A ⊗ z iI j ⊗ I j z iI j Xl=1   − ex n−1 n = n(g (Π) ˜− Π )(Ξ (Ξ ) (Ξ ) ) + Π (Ξ (Ξ ) ) A ⊗ z iI j ⊗ I j z I j = Π (Ξ (Ξ )n) nE[ξ˙ (0)ξ (0)]Π (Ξ (Ξ )n−1), z iI j − i j z iI j and thus we deduce the conclusion.  BPHZ RENORMALIZATION IN GAUSSIAN ROUGH PATHS 9

3. Application to a Rough Volatility Model In this section, we provide an example of applying BHPZ renormalization to rough paths. As discussed in [1] and [2], we set a regularity structure compatible with a rough volatility model in a financial market. In [1], the following integral is considered: T H (3.1) f(Ws )dWs, Z0 H where Wt is a Brownian motion and Wt is a fractional Brownian motion with Hurst parameter H (0, 1/2]. Here, we define the fractional Brownian motion in ∈ t Riemann–Liouville form, that is, W H := √2H KH (t s)dW where KH (t) = t 0 − s t H−1/2. | | R Remark 3.1. We do not directly define the model using fractional Brownian motion W H as a driving noise because W H does not possess stationarity. Instead, as in [2], we use a stationary noise Wˆ to approximate W H . Here, the process Wˆ is defined as Wˆ = Kˆ H ξ, where ξ is a white noise on R, and the kernel Kˆ H is smooth on R 0 R ∗and satisfies the following properties: \{ } → (1) Kˆ H = KH on [0,T ] and supp Kˆ H [0, 2T ]; (2) there exists a constant C > 0 such⊂ that for k =0, 1, 2, we have ∂kKˆ H (u) C ∂kKH (u) . | |≤ | | We construct a regularity structure (A, T, G) as follows. We define a vector space T as • S := Ξ, Ξ (Ξ)ˆ ,..., Ξ (Ξ)ˆ M , 1, (Ξ)ˆ ,..., (Ξ)ˆ M , { I I I I } T := S , h i where M := min m N; (m + 1)(H κ) 1/2 κ> 0 . We define a degree{ ∈ in T as − − − } • | · | Ξ (Ξ)ˆ m := 1/2 κ + m(H κ), m 0; | I | − − − ≥ (Ξ)ˆ m := m(H κ), m> 0; |I | − 1 := 0, | | and we define a set A := τ ; τ T . We define a structure group{| |G :=∈ Γ} ; h R as • { h ∈ } Γ 1 = 1, Γ Ξ = Ξ, Γ (Ξ)ˆ = (Ξ)ˆ + h1, h h hI I Γ (τ τ ′)=Γ (τ)Γ (τ ′), τ,τ ′ S. h · h h ∈ The next step is to construct the model on (A, T, G). We prepare the notation Wn : R2 R defined as → t Wn(s,t) := (Wˆ (r) Wˆ (s))ndξ(r), s t; − ≤ Zs n n Wn(s,t) := (Wˆ (t) Wˆ (s))iWn−i(t,s), s t. − i − ≥ i=0 X   10 H. ITO

Definition 3.2. We set the naive model (Π, Γ) in (A, T, G) as Π 1 =1, Π Ξ= ξ, Π (Ξ)ˆ n = (Wˆ ( ) Wˆ (s))n, z z zI · − d Π Ξ (Ξ)ˆ n = t Wn(s,t) , z I 7→ dt   and Γ 1 = 1, Γ Ξ = Ξ, Γ Ξ (Ξ)ˆ = (Ξ)+(ˆ Wˆ (t) Wˆ (s))1, ts ts ts I I − Γ ττ ′ =Γ τΓ τ ′, for τ, τ ′ S with ττ ′ S, ts ts ts ∈ ∈ R d Wn Wn where ξ is a white noise on , and dt (s,t) is the derivative of (s,t) in a distributional sense. Following the procedure given in Section 2, we generate a model with smoothing noises. We fix a mollifier ̺ : R R and set ̺ ( ) = ̺( /ε)/ε, ξε = ̺ ξ, and → ε · · ε ∗ ξˆε = ̺ Wˆ for ε> 0. ε ∗ Definition 3.3. We define the smooth model (Πε, Γε) in (A, T, G) as Πε 1 =1, Πε Ξ= ξ˙ε, Π (Ξ)ˆ n = (ξˆε( ) ξˆε(z))n, z z zI · − Πε Ξ (Ξ)ˆ n = ξ˙ε(ξˆε( ) ξˆε(z))n, z I · − and Γε 1 = 1, Γε Ξ = Ξ, Γε Ξ (Ξ)ˆ = (Ξ)+(ˆ ξˆε(t) ξˆε(s))1, ts ts ts I I − Γε ττ ′ =Γε τΓε τ ′, for τ, τ ′ S with ττ ′ S. ts ts ts ∈ ∈ Equations (2.10) and (2.11) are well-defined in each symbol τ T so that we can ∈ construct the BPHZ-renormalized model (Πˆ ε, Γˆε) in (A, T, G) as in Definition 2.5. By Theorem 2.7, we have the following theorem. Theorem 3.4. We have (3.2) Πˆ ε1 =1, Πˆ εΞ= ξ˙ε, Πˆ ε (Ξ)ˆ n = (ξˆε( ) ξˆε(z))n, z z zI · − (3.3) Πˆ εΞ (Ξ)ˆ n = ΠεΞ (Ξ)ˆ n nE[ξ˙ε(0)ξˆε(0)]ΠεΞ (Ξ)ˆ n−1, z I z I − z I and we also have Γˆε =Γε. ε ε Proof. By applying Theorem 2.7 to the case of d = 2, ξ1 = ξ , and ξ2 = ξˆ , we deduce the conclusion immediately.  This result corresponds to the renormalized model in [1] and [2]. In addition, BPHZ renormalization enables the model to converge to the naive model. Theorem 3.5. We have the following convergence in the norm of model : (3.4) Πˆ ε Π 0 as ε 0. ||| − ||| → → Proof. It is sufficient to prove the following claim: there exists κ> 0 such that (Π (Ξ) Πˆ ε(Ξ))(ϕλ) . λ−1/2−κεκ, s − s s

(Π ( (Ξ)ˆ Πˆ ε( (Ξ))))(ˆ ϕλ) . λH−κεκ, s I − s I s

(Π (Ξ (Ξ)ˆ n) Πˆ ε(Ξ (Ξ)ˆ n))( ϕλ) . λ−1/2+nH−κεκ, s I − s I s

BPHZ RENORMALIZATION IN GAUSSIAN ROUGH PATHS 11

λ where ϕs (t) := ϕ((t s)/λ)/λ. The proof of these inequalities is similar to that in [1]. −  Remark 3.6. In [1], the renormalized model (Π¯ ε, Γ¯ε) satisfies the following. For t [0,T ], ∈ Π¯ ε (Ξ (Ξ)n)(t) = Π Ξ (Ξ)n nE[W˙ ε(t)W H,ε(t)]Π Ξ (Ξ)n−1(t), z I z I − z I ˙ ε d ε H,ε ε H where W (t) := dt (̺ W )(t) and W (t) = (̺ W )(t). The reason for this slight difference is that∗ the model use raw noises∗ such as Brownian motion and fractional Brownian motion but not the white noises ξ and Wˆ = Kˆ H ξ. ∗ Appendix A. Trees and Forests In this article, trees and forests from graph theory are repeatedly used as symbols of regularity structures. Definition A.1. Let N be an abstract finite set and E be the subset of N N. We define maps s,t : E N as e = (v , v ) s(e)= v ,t(e)= v . × → 1 2 7→ 1 2 (1) We say that the pair T = (N, E) is a graph. The elements of N are called vertices or nodes, and the elements of E are called edges. If we want to emphasize T , then denote the sets of vertices and edges as NT and ET , respectively. (2) The pair (N, E) is called a simple directed graph if for any (v1, v2) E we have (v , v ) / E. For a simple directed graph (N, E) and (v , v ) ∈ E, v 2 1 ∈ 1 2 ∈ 1 is called the parent of v2, and v2 is called the child of v1. (3) We say that (N, E) is connected if for any v, v′ N there exists a sequence ∈′ of edges e1,...,eN such that s(e1)= v, t(eN )= v , and t(ei)= s(ei+1) (1 i N 1). ≤ (4) We≤ say− that (N, E) have a cycle if there exists a set of edges e ,...,e E 1 n ∈ such that one has t(ei)= s(ei+1) (1 i n 1) and t(en)= s(e1). (5) Let a simple directed graph T = (N,≤ E) be≤ connected− and have no cycle. If there exists a unique node v N such that t(e) = v for any e E, then the graph (N, E) is called a rooted∈ tree and the unique6 node v ∈ N is called ∈ the root of (N, E). We denote by ̺T the root of T . (6) Let L be a finite set, (N, E) be a rooted tree, and Type : E L be a map. The triplet (N,E, Type) is called a typed tree. → (7) For any rooted trees T ,T , we define T = T T as 1 2 1 · 2 N := N ′ N ′ ̺ , T T1 ⊔ T2 ⊔{ T } E := E′ E′ (̺ , v); v N N T T1 ⊔ T2 ⊔{ T ∈ T1 ⊔ T2 such that (̺ , v) E or (̺ , v) E T1 ∈ T1 T2 ∈ T2 } ′ ′ ′ where NTi := NTi ̺Ti , ETi := e = (v, v ) ETi ; v = ̺Ti ,i = 1, 2. This operation T ,T\{ T} T is called{ a tree product.∈ 6 } 1 2 7→ 1 · 2 Throughout this article, we assume that a rooted tree has no node such that it has two or more parents, that is, for any rooted tree τ and any x Nτ , we have #(t−1( x )) 1. ∈ { } ≤ Example A.2. A rooted tree can be described as placing the root on the bottom as follows: T1 = 12 H. ITO

T2 =

The tree product of T1 and T2 is expressed as

T T = 1 · 2

However, the graphs T3 and T4 defined below are not rooted trees because T3 has two candidates for the root and T4 has a loop.

T3 =

T4 =

Definition A.3. (1) A graph in which each connected component is a tree is called a forest. In addition, the forest is called a rooted forest if each connected component of it is a rooted tree. (2) Let L be a finite set. The triplet (N,E, Type) such that (N, E) is a rooted forest and Type : E L is called a typed rooted forest. → Regarding the empty set as a forest, a vector space generated by the totality of rooted forests is a unital algebra∅ in which the product is the disjoint union and the unit is the empty set. This product is called a forest product. We denote by 1 the unit of a forest product unless confusion occurs. Definition A.4. A colored forest is a pair (F, Fˆ) satisfying the following conditions.

(1) F = (EF ,NF , Type) is a typed rooted forest. (2) The map Fˆ : E N N 0 satisfies the following: for any e = (x, y) F ⊔ F → ∪{ } ∈ E such that Fˆ(e) = 0, we have Fˆ(x)= Fˆ(y)= Fˆ(e). F 6 Fˆ is called the coloring of F . For i> 0, we set Fˆ = (Eˆ , Nˆ ), Eˆ = Fˆ−1(i) E , Nˆ = Fˆ−1(i) N . i i i i ∩ F i ∩ F Finally, we denote by CF the set of all colored forests.

Appendix B. BPHZ Renormalization of Γ First, we introduce a positive renormalization, which is the foundationofΓ. A map ∆ : T T T is defined as + → ⊗ ∆ 1 = 1 1, ∆ Ξ = Ξ 1 + 1 Ξ , + ⊗ + i i ⊗ ⊗ i ∆ = 1 + 1 , +I I⊗ ⊗ I ∆ (Ξ )= (Ξ ) 1 + Ξ + 1 (Ξ ), +I i I i ⊗ I⊗ i ⊗ I i ∆ (τ τ ′) = ∆ τ ∆ τ ′ + · + · + for τ, τ ′ S such that τ τ ′ S, and the domain is extended by imposing linearity. Note that∈ a product on T· T∈can be defined as (τ τ ) (τ ′ τ ′ ) := (τ τ ′ ) (τ τ ′ ). ⊗ 1 ⊗ 2 · 1 ⊗ 2 1 · 1 ⊗ 2 · 2 BPHZ RENORMALIZATION IN GAUSSIAN ROUGH PATHS 13

Definition B.1. (1) A subspace J− is defined as

J− := τ T ; τ = σ σ,¯ σ, σ¯ T, σ =1 , σ 0 , h{ ∈ · ∈ 6 2 | |≤ }i where is a tree product. (2) A quotient· space generated by the above subspace is defined by

T+ := T/J−, ex and pG : T T+ denotes the canonical projection. (3) A map ∆ex →: T T T is defined as ∆ex := (Id pex)∆ . + → ⊗ + + ⊗ + + ex As defined, the vector space T+ and the map ∆+ can be given as T = 1, n, (Ξ )n;1 i d, n N + h{ I I i ≤ ≤ ∈ } and ∆ex1 = 1 1, ∆exΞ = Ξ 1, + ⊗ + i i ⊗ ∆ex = 1 + 1 , + I I⊗ ⊗ I ∆ex (Ξ )= (Ξ ) 1 + 1 (Ξ ). + I i I i ⊗ ⊗ I i Definition B.2. For s,t R, a map γ : T R is defined as ∈ ts + → γ 1 =1, γ = Π( )(t) Π( )(s), ts tsI I − I γ (Ξ )= Π(Ξ )(t) Π(Ξ )(s), tsI i i − i γ (τ τ ′)= γ τ γ τ ′ ts · ts · ts for τ, τ ′ S such that τ τ ′ T . ∈ · ∈ + Note that Γ = (Id γ )∆ex. The BPHZ renormalization of Γ is as follows. ts ⊗ ts + Definition B.3. Let a pair (Π, Γ) be a model on the regularity structure (A, T, G). The BPHZ-renormalized model (Πˆ , Γ)ˆ is defined by − ex ex Πˆ τ = (g (Π) ˜− Π )∆− τ, s A ⊗ z − ex ex Γˆ τ = (Id γ (g (Π) Id)∆− )∆ τ. ts ⊗ ts ⊗ + The next proposition is the main subject of this subsection. Proposition B.4. We have (B.1) Γˆ τ =Γ τ, τ T. ts ts ∈ Proof. The case where τ = 1 or = Ξi is obvious. We consider the case where τ = n. Then, we have I ∆ex n = ( 1 + 1 )n + I I⊗ ⊗ I n n = ( i n−l). l I ⊗ I Xl=0   ex l l The edge has positive degree = 1, so we have ∆− = 1 . Thus, we have I |I| I ⊗ I n n − ex n i n−l Γˆ = (Id γ (g (Π) Id)∆− ) ( ) tsI ⊗ ts ⊗ l I ⊗ I l=0   n X n = (Id γ ) ( i n−l) ⊗ ts l I ⊗ I Xl=0   =Γ n. tsI 14 H. ITO

Next, we show the claim when τ = Ξ (Ξ )n. Then, we have iI j ∆exΞ (Ξ )n = ∆exΞ ∆ex (Ξ )n + iI j + i · + I j = (Ξ 1) ( (Ξ ) 1 + 1 (Ξ ))n i ⊗ · I j ⊗ ⊗ I j n n = Ξ (Ξ )l (Ξ )n−l. l iI j ⊗ I j Xl=0   − ex ex l l For l 1, we have (g (Π) ˜− Id)∆− ( (Ξj ) ) = (Ξj ) because there exists no ≥′ m A ′ ⊗ I− ex I′ tree τ (Ξ ) such that τ < 0 and g (Π) ˜− τ = 0. Thus, we have ⊂ I j | | A 6 n n − ex ex n l n−l Γˆ (Ξ (Ξ ) ) = (Id γ (g (Π) ˜− Id)∆− ) (Ξ (Ξ ) (Ξ ) ) zz¯ iI j ⊗ zz¯ A ⊗ l iI j ⊗ I j Xl=0   n n = (Id γ ) (Ξ (Ξ )l (Ξ )n−l) ⊗ zz¯ l iI j ⊗ I j Xl=0   =Γ (Ξ (Ξ )n). zz¯ iI j 

References

[1] C. Bayer, P. K. Friz, P. Gassiat, J. Martin, and B. Stemper. A regularity structure for rough volatility. Mathematical Finance, 30(3):782–832, 2020. [2] C. Bellingeri, P. K. Friz, and M. Gerencs´er. Singular paths spaces and applications. arXiv preprint arXiv:2003.03352, 2020. [3] Y. Bruned. Renormalisation from non-geometric to geometric rough paths. arXiv preprint arXiv:2007.14385, 2020. [4] Y. Bruned, I. Chevyrev, and P. K. Friz. Examples of renormalized SDEs. In International Conference on Stochastic Partial Differential Equations and Related Fields, pages 303–317. Springer, 2016. [5] Y. Bruned, I. Chevyrev, P. K. Friz, and R. Preiss. A rough path perspective on renormalization. Journal of Functional Analysis, 277(11):108–283, 2019. [6] Y. Bruned, M. Hairer, and L. Zambotti. Algebraic renormalisation of regularity structures. Inventiones mathematicae, 215(3):1039–1156, 2019. [7] A. Chandra and M. Hairer. An analytic BPHZ theorem for regularity structures. arXiv preprint arXiv:1612.08138, 2016. [8] M. Hairer. A theory of regularity structures. Inventiones mathematicae, 198(2):269–504, 2014. [9] S. Janson. Gaussian Hilbert spaces. Number 129 in Cambridge tracts in mathematics. Cambridge University Press, Cambridge, United Kingdom, 1997.

Graduate School of Engineering Science, Osaka University, Osaka 560-8531, Japan Email address: [email protected]