BPHZ RENORMALIZATION IN GAUSSIAN ROUGH PATHS HAYAHIDE ITO Abstract. We construct a procedure for Bogoliubov–Parasiuk–Hepp– Zimmermann (BPHZ) renormalization of a rough path in view of the relation between rough path theory and regularity structure. We also provide a plain expression of the BPHZ-renormalized model in a rough path. BPHZ renormalization plays a central role in the theory of singular stochastic partial differential equations and assures the convergence of the model in the regularity structure. Here we demonstrate that the renormalization is also effective in a rough path setting by considering its application to rough path theory and mathematical finance. 1. Introduction We consider the following Itˆo-type stochastic integral: for T > 0 and a smooth function f : R R, → T H (1.1) f(Wt )dWt, Z0 H where Wt is a standard one-dimensional Brownian motion, and Wt is a fractional Brownian motion with Hurst parameter H (0, 1/2) given by ∈ t W H = √2H t s H−1/2dW . t | − | s Z0 Integral (1.1) emerges in option pricing in a financial market with rough volatility (see [1]), so the Wong–Zakai approximation of (1.1) is valuable for financial mathematics. However, for any mollifier ̺ : R R, the limit → T H,ε ε lim f(Wt )dWt ε→0 Z0 arXiv:2103.07044v1 [math.PR] 12 Mar 2021 ε H,ε ε H ε does not exist in general, where ̺ ( ) = ̺( /ε)/ε, Wt = (̺ W )t, and Wt = (̺ε W ) . This failure is due to fractional· · Brownian motion∗ W H having lower ∗ t t regularity than Brownian motion Wt. Bayer et al. [1] demonstrated the following revised Wong–Zakai approximation by adding a term to prevent the divergence. Theorem 1.1. Let C ε(t)= E[W˙ ε(t)W H,ε(t)]. Then, we have the following limit: T T T (1.2) f(W H,ε)dW ε C ε(t)f ′(W H,ε)dt f(W H )dW , ε 0. t t − t → t t → Z0 Z0 Z0 The result is proved by the theory of regularity structures by Hairer [8]. The following sketch illustrates the relations between the convergence of the integral 2020 Mathematics Subject Classification. Primary 60L30, Secondary 60L20, 91-10. 1 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 τ.