# Fluctuations Around a Homogenised Semilinear Random PDE

Martin Hairer, Etienne Pardoux

► To cite this version:

Martin Hairer, Etienne Pardoux. Fluctuations Around a Homogenised Semilinear Random PDE. Archive for Rational Mechanics and Analysis, 2021, 239 (1), pp.151-217. 10.1007/s00205-020-01574-8. hal-03119140

HAL Id: hal-03119140

https://hal.science/hal-03119140

Submitted on 25 Jan 2021

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

# Fluctuations Around a Homogenised Semilinear Random PDE

MARTIN HAIRER & ÉTIENNE PARDOUX

Communicated by F. OTTO

## Abstract

We consider a semilinear parabolic partial differential equation in $R_+ \times [0, 1]^d$, where $d = 1, 2$ or $3$, with a highly oscillating random potential and either homogeneous Dirichlet or Neumann boundary condition. If the amplitude of the oscillations has the right size compared to its typical spatiotemporal scale, then the solution of our equation converges to the solution of a deterministic homogenised parabolic PDE, which is a form of law of large numbers. Our main interest is in the associated central limit theorem. Namely, we study the limit of a properly rescaled difference between the initial random solution and its LLN limit. In dimension $d = 1$, that rescaled difference converges as one might expect to a centred Ornstein–Uhlenbeck process. However, in dimension $d = 2$, the limit is a non-centred Gaussian process, while in dimension $d = 3$, before taking the CLT limit, we need to subtract at an intermediate scale the solution of a deterministic parabolic PDE, subject (in the case of Neumann boundary condition) to a *non-homogeneous* Neumann boundary condition. Our proofs make use of the theory of regularity structures, in particular of the very recently developed methodology allowing to treat parabolic PDEs with boundary conditions within that theory.

## Contents

<table><tbody><tr><td>1. Introduction</td><td>152</td></tr><tr><td>2. Assumptions on the Noise</td><td>157</td></tr><tr><td>2.1. Coalescence Trees</td><td>157</td></tr><tr><td>2.2. Justification</td><td>159</td></tr><tr><td>3. Law of Large Numbers</td><td>162</td></tr><tr><td>4. Central Limit Theorem</td><td>166</td></tr><tr><td>4.1. Decomposition of the Solution</td><td>167</td></tr><tr><td>4.2. Definition of the Ambient Regularity Structure</td><td>169</td></tr><tr><td>4.3. Description of the Models</td><td>171</td></tr><tr><td>4.4. Description of v<sub>e</sub><sup>(0)</sup></td><td>178</td></tr></tbody></table>
<table><tbody><tr><td>4.5. Description of v<sub>ε</sub><sup>(1)</sup></td><td>180</td></tr><tr><td>4.6. Formulation of the Fixed Point Problem</td><td>183</td></tr><tr><td>5. Convergence of Models</td><td>188</td></tr><tr><td>5.1. Cumulant Homogeneity Assignments</td><td>188</td></tr><tr><td>5.2. Power-Counting Conditions</td><td>191</td></tr><tr><td>5.3. Special Bounds</td><td>192</td></tr><tr><td>5.4. Tightness and Convergence for the Noise</td><td>193</td></tr><tr><td>5.5. Boundary Term</td><td>196</td></tr><tr><td>5.6. Bounds on the Function Q</td><td>204</td></tr><tr><td>6. Auxiliary Results</td><td>205</td></tr><tr><td>References</td><td>215</td></tr></tbody></table>

# 1. Introduction

Fix $D = [0, 1]^d$ with $d \le 3$, and consider the family of functions $u_\epsilon : [0, T] \times D \to \mathbf{R}$ solving the PDE

$$ \partial_t u_\epsilon(t, x) = \Delta u_\epsilon(t, x) + H(u_\epsilon(t, x)) + G(u_\epsilon(t, x))\eta_\epsilon(t, x), \quad u_\epsilon(0, x) = u_0(x), \qquad (1.1) $$

endowed with either Dirichlet boundary conditions $u_\epsilon(t, x) = 0$ for $x \in \partial D$ or Neumann boundary conditions $\langle n(x), \nabla u_\epsilon(t, x) \rangle = 0$, where $n$ denotes the outward facing unit vector normal to the boundary of $D$. The driving noise $\eta_\epsilon$ appearing in this equation is given by

$$ \eta_{\epsilon}(t, x) = \epsilon^{-1} \eta(\epsilon^{-2} t, \epsilon^{-1} x), \qquad (1.2) $$

where $\eta(t, x)$ is a stationary centred random field, which we do not assume Gaussian, but with relatively good mixing properties (see Assumption 2.1 for details) and moments of all orders after testing against a test function. Note that $\eta_\epsilon$ is scaled by $\epsilon^{-1}$ rather than $\epsilon^{-(d+2)/2}$, so the noise from [10, 14] (which were restricted to $d=1$) has been multiplied by $\epsilon^{d/2}$.

In the case when $G$ is linear and $H = 0$, this problem has been well studied. See for example [1] for the case where furthermore $\eta$ is Gaussian and constant in time [note that in our case $m = 2$ so the exponent $\alpha$ appearing there would equal 1 in our case, as used in (1.2)], [4] for a similar result in the non-Gaussian case.

Although we will allow $\eta$ to be a generalised random field in $d=1$, we assume throughout that there exist locally integrable functions $\kappa_p : (\mathbf{R}^{d+1})^p \to \mathbf{R}$ that are continuous outside of the big diagonal $\Delta_p = \{(z_1, \dots, z_p) \in (\mathbf{R}^{d+1})^p : \exists i \ne j \text{ with } z_i = z_j\}$ and such that, for any $C_0^\infty$ test functions $\phi_1, \dots, \phi_p$, the joint cumulant $\kappa_p(\phi_1, \dots, \phi_p)$ of $\eta(\phi_1), \dots, \eta(\phi_p)$ satisfies

$$ \kappa_p(\phi_1, \dots, \phi_p) = \int \kappa_p(z_1, \dots, z_p) \phi_1(z_1) \cdots \phi_p(z_p) \mathrm{d}z_1 \cdots \mathrm{d}z_p. $$

(By stationarity, the functions $\kappa_p$ only depend on the differences of their arguments.) Here and below, we always use the convention that $z$ (resp. $z_i$, $\bar{z}$, etc) denotes a space-time coordinate given by $z = (t, x)$ (resp. $z_i = (t_i, x_i)$, $\bar{z} = (\bar{t}, \bar{x})$, etc). We
furthermore normalise our problem by assuming that the covariance of $\eta$ integrates to 1 in the sense that

$$ \int_{\mathbb{R}^{d+1}} \kappa_2(0, z) \, dz = 1. \qquad (1.3) $$

In particular, we assume that $\kappa_2(0, \cdot)$ is absolutely integrable, but this will in any case follow from Assumption 2.1. Furthermore, define the constants

$$ \begin{aligned} \circled{0} = \int P(z) \kappa_2(0, z) \, dz, & \qquad \circled{0} = \int P(z) x\kappa_2(0, z) \, dz, \\ \circled{00} = \int P(z) P(z') \kappa_3(0, z, z') \, dz \, dz', & \qquad \circled{00} = \int P(z) P(z' - z) \kappa_3(0, z, z') \, dz \, dz', \end{aligned} $$

where $P$ denotes the heat kernel, that is the fundamental solution to the heat equation on the whole space. Here and below, symbols drawn in red denote fixed constants, while symbols drawn in blue will later denote basis vectors of a suitable regularity structure associated to our problem. In dimension $d=1$, our assumptions on $\kappa_2$ will guarantee that the integral $\circled{0}$ converges absolutely, while in dimensions 2 and 3 our assumptions on $\kappa_2$ and $\kappa_3$ will guarantee that all of these integrals converge absolutely.

The scaling $\varepsilon^{-1}$ chosen in $\eta_\varepsilon$ is such that $u_\varepsilon$ converges as $\varepsilon \to 0$ to a limit $u^{(0)}$, which is our first result. Indeed, writing

$$ H_{\eta}(u) = H(u) + \circled{0} G'(u)G(u), \qquad (1.4) $$

we have the following “law of large numbers”:

**Theorem 1.1.** Let $u_\varepsilon$ be as above and let $u^{(0)}$ be the (local) solution to the deterministic PDE

$$ \partial_t u^{(0)} = \Delta u^{(0)} + H_{\eta}(u^{(0)}), \qquad (1.5) $$

with the same initial condition $u_0 \in C^\alpha$ as (1.1) [for some arbitrary $\alpha \in (0, 1]$] and with homogeneous Dirichlet (resp. Neumann) boundary condition. In the case of Dirichlet boundary conditions, we impose that $u_0$ vanishes on the boundary.

Assume that the functions $G, H: \mathbb{R} \to \mathbb{R}$ are of class $C^5$ and $C^4$ respectively, that the driving field $\eta$ satisfies Assumption 2.1, and let furthermore $T > 0$ be such that the (possible) explosion time for $u^{(0)}$ is greater than $T$. Then, in probability and uniformly over $[0, T] \times D$, $u_\varepsilon$ converges to $u^{(0)}$ as $\varepsilon \to 0$.

The proof of this result will be given in Section 3. Our main quantity of interest however are the fluctuations of $u_\varepsilon$ around its limit $u^{(0)}$. One interesting feature of this problem is that in order to see these fluctuations, it is not sufficient to recenter around $u^{(0)}$. Instead, as soon as $d \ge 2$, a suitable first-order correction $u^{(1)}$ living at scale $\varepsilon$ has to be subtracted first. Our precise “central limit theorem” then takes the following form:
**Theorem 1.2.** Let $u_0$ be such that its extension to all of $\mathbf{R}^d$ by reflections<fnref n="1" /> is of class $C^3$, let $u_\epsilon$ and $u^{(0)}$ be as above and assume $G, H, T$ and $\eta$ are as in Theorem 1.1. Let furthermore $u^{(1)} = 0$ for $d = 1$ and, for $d \in \{2, 3\}$, let $u^{(1)}$ be the solution to

$$ 
\begin{aligned}
\partial_t u^{(1)} &= \Delta u^{(1)} + H'_{\eta}(u^{(0)})u^{(1)} + \Psi(u^{(0)}, \nabla u^{(0)}), \\
\Psi(u, p) &= \frac{1}{2} \mathring{\mathring{\mathring{G}}}^2(G^2 G'')(u) + \mathring{\mathring{\mathring{G}}}^2(G(G')^2)(u) + (G')^2(u) \langle \mathring{\mathring{\mathring{G}}}, p \rangle, \quad (1.6)
\end{aligned}
$$

with zero initial condition. (Note that $\mathring{\mathring{\mathring{G}}}$ is an $\mathbf{R}^d$-valued constant.) In the case of Neumann boundary condition, we furthermore impose that $\langle \nabla u^{(1)}(t, x), n(x) \rangle = c(x)GG'(u^{(0)}(t, x))$ for $x \in \partial D$, where $c$ is an explicit function on the boundary of $D$ which is constant on each of its faces [the precise values of $c$ on each face will be given in (5.16)], while we impose homogeneous Dirichlet boundary conditions otherwise.

Then, in law and in $C^\alpha([0, T] \times D)$ for any $\alpha < 1 - \frac{d}{2}$, one has

$$ \lim_{\epsilon \to 0} \frac{u_{\epsilon} - u^{(0)} - \epsilon u^{(1)}}{\epsilon^{d/2}} = v, \quad (1.7) $$

where $v$ is the Gaussian process solving

$$ \partial_t v = \Delta v + H'_{\eta}(u^{(0)})v + G(u^{(0)})\xi, \quad (1.8) $$

endowed with homogeneous Dirichlet (resp. Neumann) boundary condition and 0 initial condition, and $\xi$ denotes a standard space-time white noise.

*Proof.* Combining Proposition 4.18 with Proposition 4.16 and (4.1) shows that if we set $v_\epsilon = \epsilon^{-d/2}(u_\epsilon - u^{(0)} - \epsilon u^{(1)})$ then we do indeed have $\lim_{\epsilon \to 0} v_\epsilon = v$ (weakly in $C^{(d-2)/2-\kappa}$ on $[0, T]$). The limit $v$ is identified as the solution to (1.8) by combining the second part of Proposition 4.16 with Lemma 4.10. $\square$

**Remark 1.3.** If all we were interested in is the law of large numbers, then the conditions of Assumption 2.1 on $\eta$ could easily be weakened.

**Remark 1.4.** In the case of Neumann boundary conditions, it may appear paradoxal that, even though $u_\epsilon$, $u^{(0)}$ and $v$ all satisfy homogeneous boundary conditions, $u^{(1)}$ does not! This phenomenon is very similar to the presence of the “boundary renormalisation” that can appear in the context of singular SPDEs [9]. There is no contradiction since the convergence $v_\epsilon \to v$ takes place in a very weak topology in which the notion of “normal derivative at the boundary” is meaningless in a pointwise sense. (A very simple example displaying a similar phenomenon is $n^{-1/2} \sin(nx)$, whose derivative at the origin diverges like $\sqrt{n}$ while that of its limit vanishes.)

<footnote n="1">¹ Here we perform the reflections consistent with the reflection principle for our choice of boundary conditions.</footnote>
**Remark 1.5.** Regarding the precise meaning of the equation fulfilled by $u^{(1)}$ in the case of Neumann boundary condition, denote by $\delta_{\partial D}$ the distribution on $\mathbf{R} \times \partial D$ given by

$$ \delta_{\partial D}(\phi) = \int_{\mathbf{R}} \int_{\partial D} \phi(t, x) \, dx \, dt, $$

where the integration over the faces of $\partial D$ are performed against the two-dimensional Lebesgue measure. With this notation, the solution to any equation of the form

$$ \partial_t u = \Delta u + f \text{ in } D, \quad \langle \nabla u(t, x), n(x) \rangle = g(x) \text{ on } \partial D, \quad u(0, x) = u_0(x), \tag{1.9} $$

(and in particular the equation determining $u^{(1)}$) is defined as the solution to the integral equation

$$ u(z) = \int_D P_{\text{Neu}}(z, (0, x')) u_0(x') \, dx' + \int_{\mathbf{R}^4} P_{\text{Neu}}(z, z') (f \mathbf{1}_+^D + g \delta_{\partial D}) (dz'), \tag{1.10} $$

where $P_{\text{Neu}}$ denotes the homogeneous Neumann heat kernel, with the convention that $g(t, x) = 0$ for $t \le 0$, and where $\mathbf{1}_+(t, x) = \mathbf{1}_{\{t \ge 0\}}\mathbf{1}_{\{x \in D\}}$. Here, we used the notation $\int \phi(z) \eta(dz)$ for the usual pairing between a distribution $\eta$ and a suitable test function $\phi$. To see that solutions to (1.10) and (1.9) do indeed coincide if $f$ and $g$ are sufficiently regular for the solution to be differentiable up to the boundary, it suffices to note that the mild formulation is equivalent to the weak formulation, see for example [20], with the term $g\delta_{\partial D}$ appearing as the boundary term when integrating by parts.

**Remark 1.6.** In dimension $d=1$, the term $u^{(1)}$ in (1.7) is of course redundant. In dimension $d=2$, it is still the case that $\varepsilon^{-d/2}(u_\varepsilon - u^{(0)})$ converges to a limit, but this limit is not centred anymore. In higher dimensions, additional corrections appear. We expect to have a result of the form

$$ \lim_{\varepsilon \to 0} \varepsilon^{-d/2} \left( u_{\varepsilon} - \sum_{k=0}^{\lfloor d/2 \rfloor} \varepsilon^k u^{(k)} \right) = v, \tag{1.11} $$

where $u^{(0)}$ is as above and the $\bar{u}^{(k)}$ satisfy an equation of the type

$$ \partial_t u^{(k)} = \Delta u^{(k)} + H'_\eta(u^{(0)})u^{(k)} + \Psi_k, $$

for some inhomogeneity $\Psi_k$ depending on the $u^{(\ell)}$ for $\ell < k$ and some of their derivatives. Since $v$ has vanishing expectation, we expect to also have

$$ \lim_{\varepsilon \to 0} \varepsilon^{-d/2} (u_\varepsilon - \mathbf{E}u_\varepsilon) = v, $$

so that the $u^{(k)}$ provide an expansion of $\mathbf{E}u_\varepsilon$ in powers of $\varepsilon$. Note however that the techniques used in this paper do not provide moment bounds on the solution, so that even in $d \le 3$ this would require some additional work.
**Remark 1.7.** The form of the terms appearing in the successive correctors as well as the constants multiplying them can in principle be derived from [2, Eq. 2.12] which describes the form of the counterterms $\Upsilon$ associated to a given tree $\tau$. The recursion given there suggests a correspondence $\circ \sim G(u)$ with incoming edges corresponding to functional derivatives with respect to $u$ (so $\circ \sim G'(u)$ and $\circ \sim G''(u)$) and powers of $X$ corresponding to formal directional derivatives, so $\otimes \sim G'(u)\nabla u$. This shows a priori that the counterterm multiplying $\circled$ for example must be of the form $G'(u)G(u)$, etc. The numerical constants multiplying these terms do however differ from those appearing in [2] since their meaning is slightly different: in [2] we use counterterms to “recenter” the original equation in order to obtain a finite limit while here we leave the original equation untouched and compute its “centering”. If we leave aside the behaviour at the boundary, this in principle allows to guess the general form of the equations for the $u^{(k)}$ appearing in (1.11) for any dimension.

The most surprising part of Theorem 1.2 is surely the fact that in the case when $u_\epsilon$ has homogeneous Neumann boundary conditions, even though $v$ and $u^{(0)}$ both also have homogeneous boundary conditions, $u^{(1)}$ does not, which seems to contradict (1.7). This is of course not a contradiction but merely suggests that if we write $v_\epsilon$ for the expression appearing under the limit in (1.7), then $v_\epsilon$ exhibits a kind of boundary layer. Note also that the statement that “$v$ satisfies homogeneous boundary conditions” only makes sense in terms of the integral equation that it solves since $v$ itself is not differentiable at the boundary. (It is not even a function!)

Before we proceed, let us give a heuristic explanation for the appearance of this boundary layer. Consider the simplest case $H = 0$, $G(u) = u$ and $u_0 > 0$, in which case we can consider the Hopf-Cole transform $h_\epsilon = \log u_\epsilon$, yielding

$$ \partial_t h_\epsilon = \Delta h_\epsilon + |\nabla h_\epsilon|^2 + \eta_\epsilon. $$

To leading order, one would expect the right hand side to behave like $|\nabla h_\epsilon|^2 \simeq E|\nabla Z_\epsilon|^2$, where $Z_\epsilon$ solves $\partial_t Z_\epsilon = \Delta Z_\epsilon + \eta_\epsilon$ endowed with homogeneous Neumann boundary conditions. It turns out that, in the interior of the domain, one has

$$ \lim_{\epsilon \to 0} E|\nabla Z_\epsilon|^2 = \circled, $$

which allows one to “guess” the correct limit $u^{(0)}$. On the boundary however $\nabla Z_\epsilon = 0$, so that one expects $E|\nabla Z_\epsilon|^2 - \circled$ to be of order $O(1)$ in a layer of width $O(\epsilon)$ around $\partial D$. When going to the next scale, this results in a boundary correction of order $O(\epsilon^{-1})$ in this boundary layer, which precisely scales like a surface measure on the boundary. Remark 1.5 shows that the net effect of this correction is to modify the boundary condition.

The remainder of the article is structured as follows. First, in Section 2, we formulate our main assumption on the driving noise $\eta$ and we show that this assumption is “reasonable” by exhibiting an explicit class of examples for which it is satisfied. In Section 3, we then show that the law of large numbers holds. Although this could probably be shown by “classical” means without too much effort, we will use the theory of regularity structures because it shortens the argument and allows
us to introduce some results and notation that will be of use later on. In Section 4, we then show that the central limit theorem holds. The main tool in this proof is the convergence of a certain “model” for an appropriate regularity structure as well as refinements of the type of boundary estimates first considered in [9]. The convergence of the model is given in Section 5. “Appendix B” is devoted to the proof of a result showing that the operations of “convolution by a singular kernel” and “multiplication by a smooth function” almost commute, modulo a much smoother remainder, a fact that will undoubtedly sound familiar to anyone acquainted with microlocal analysis. “Appendix C” contains a version of the reconstruction theorem that is purpose-built to allow us to deal with modelled distribution that have very singular boundary behaviour and goes beyond the version obtained in [9]. This appendix was written in collaboration with Máté Gerencsér.

## 2. Assumptions on the Noise

In this section, we formulate our precise assumptions on the driving noise and we show that they are satisfied for example by a mollified Poisson process. In a nutshell, we want to assume that correlations are bounded by $||z - \bar{z}||^{-2\epsilon}$ at small scales and $||z - \bar{z}||^{-2\bar{c}}$ at large scales with $c = \frac{1}{2} - \delta$ and $\bar{c} = \frac{d+2}{2} + \delta$ for some $\delta \in (0, \frac{1}{2})$. However, we also want to encode the fact that higher order cumulants behave “better” than what is obtained from simply using the Cauchy–Schwartz inequality. Note that our assumptions are trivially satisfied by any continuous Gaussian process with correlations that decay at least like $||z - \bar{z}||^{-2\bar{c}}$. Here and below, $||\cdot||$ will always denote the parabolic distance between space-time points. It will be convenient (in particular in “Appendix A”) to make sure that $||\cdot||$ is smooth away from the origin, so we set for example $||z||^4 = ||(t, x)||^4 = |x|^4 + |t|^2$.

### 2.1. Coalescence Trees

In order to formulate this precisely, we need a simplified version of the construction of [15, Appendix A]. Given any configuration $(z_1, \dots, z_p)$ of $p$ points in $\mathbb{R}^{d+1}$ with all distances distinct, we associate to it a binary tree $T$ in the following way. Consider Kruskal’s algorithm [17] for constructing the minimal spanning tree of the complete graph with vertices $\{z_1, \dots, z_p\}$ and edge-weights given by their (parabolic) distances. One way of formalising this is the following. Consider the set $\mathcal{P}_p$ of partitions of $\Omega_p = \{1, \dots, p\}$. We define a distance $d_z$ between subsets of $\Omega_p$ as the Hausdorff distance induced by $\{z_1, \dots, z_p\}$, namely

$$d_z(A, B) = \max \left\{ \sup_{i \in A} \inf_{j \in B} ||z_i - z_j||, \sup_{j \in B} \inf_{i \in A} ||z_i - z_j|| \right\}.$$

We then define a map $K: \mathcal{P}_p \to \mathcal{P}_p$ in the following way. If $\pi = \{\Omega_p\}$, then $K(\pi) = \pi$. Otherwise, let $A \neq B \in \pi$ be such that $d_z(A, B) \le d_z(C, D)$ for all $C, D \in \pi$. Thanks to our assumption on the $z_i$, this pair is necessarily unique. We then set

$$K(\pi) = (\pi \setminus \{A, B\}) \cup \{A \cup B\},$$
that is $K(\pi)$ is obtained by coalescing the two sets $A$ and $B$ in the partition $\pi$. The vertices of $T$ are then given by $V_T = \bigcup_{n \ge 0} K^n(\{\{1\}, \dots, \{p\}\})$, that is $V_T$ consists of all the blocks of those partitions. The set $V_T$ comes with a natural partial order given by inclusion: $A \le B$ if and only if $A \supset B$. The (directed) edge set $E_T \in V_T \times V_T$ of $T$ is then given by the Hasse diagram of $(V_T, \le): (A, B) \in E_T$ if and only if $A < B$ and there is no $C \in V_T$ such that $A < C < B$. It is easy to verify that $T$ is a binary tree and that its leaves are precisely given by the singletons. It will be convenient to also add to $V_T$ a “point at infinity” $\otimes$ which is connected to $\Omega_p$ by an edge $(\Omega_p, \otimes)$ and to view $\otimes$ as the minimal element of $V_T$.

We write $\mathring{V}_T = V_T \setminus \{\{1\}, \dots, \{p\}, \otimes\}$ for the interior nodes. Each interior node $A \in \mathring{V}_T$ has exactly two children $A_1$ and $A_2$ such that $(A, A_i) \in E_T$ for $i = 1, 2$. We then define an integer labelling $\mathbf{n}: \mathring{V}_T \to \mathbf{Z}$ by $\mathbf{n}(A) = -\lceil \log_2 d_z(A_1, A_2) \rceil$. We will always view $\mathbf{n}$ as a function on all of $V_T$ with values in $\mathbf{Z} \cup \{\pm\infty\}$ by setting $\mathbf{n}(\otimes) = -\infty$ and $\mathbf{n}(\{i\}) = +\infty$ for $i = 1, \dots, p$. Note now that if $A, B$ and $C$ are three disjoint sets, then

$$d_z(A, B) \le \min\{d_z(A, C), d_z(B, C)\} \Rightarrow d_z(A, B) \le d_z(A \cup B, C).$$

As a consequence, the map $\mathbf{n}$ is monotone increasing on $V_T$. Furthermore, as in [15, Eq. A.15], there exist constants $c, C$ depending only on $p$ such that

$$c2^{-\mathbf{n}(\{i\} \wedge \{j\})} \le \|z_i - z_j\| \le C2^{-\mathbf{n}(\{i\} \wedge \{j\})},$$

for all $i, j$.

Given a configuration of points $z = (z_1, \dots, z_p) \in (\mathbf{R}^{d+1})^p$, we now write $\mathbf{t}_z = (T, \mathbf{n})$ for the corresponding data constructed as above. We furthermore define a function $\rho: \mathbf{R}_+ \to \mathbf{R}_+$ by

$$\rho(r) = r^{-\bar{c}} \wedge r^{-c}.$$

(Beware that $\rho$ is an upper bound for the square root of the covariance between two points.) We then assume that the following bound holds:

**Assumption 2.1.** With the notations as above, for any $p \ge 2$ and any $\{k_i\}_{i=1}^p \subset \mathbf{Z}_+^{d+1}$, the $p$th joint cumulant for $\eta$ satisfies the bound

$$\left| \left( \prod_{i=1}^{p} D_i^{k_i} \right) \kappa_p(z_1, \dots, z_p) \right| \lesssim \rho(2^{-\mathbf{n}(\Omega_p)}) \prod_{A \in \mathring{V}_T} \rho(2^{-\mathbf{n}(A)}) \prod_{i=1}^{p} 2^{|k_i|} \mathbf{n}(i^\top), \quad (2.1)$$

uniformly over all $z \in (\mathbf{R}^{d+1})^p$. (Recall that $\Omega_p$ is the root of the tree $T$.) Here and below, the length of the multiindex $k$ should be interpreted in the parabolic sense, namely $|k| = 2k_0 + \sum_{i=1}^d k_i$.

In dimensions $d \in \{2, 3\}$, we furthermore assume that $\eta: \Omega \times \mathbf{R}^{d+1} \to \mathbf{R}$ is a measurable function with $\mathbf{E}|\eta(0)|^p < \infty$ for $p = (d+2)/c$.

**Remark 2.2.** The additional condition that $\eta$ takes values in $L^p$ for sufficiently high $p$ is mainly technical and could probably be dropped with some additional effort. It will be used to bound $\hat{R}_\epsilon^{(d)}$ in the proof of Proposition 4.16. The exponent $(d+2)/c$ is consistent with the condition on the correlation function in the sense that this is the lowest value of $p$ for which $L_{\text{loc}}^p \subset C^{-c}$.
Note also that the cumulants $\kappa_p^{(\epsilon)}$ of the rescaled process $\eta_\epsilon$ satisfy

$$ \kappa_p^{(\epsilon)}(z_1, \dots, z_p) = \epsilon^{-p} \kappa_p(S_\epsilon z_1, \dots, S_\epsilon z_p), \quad (2.2) $$

where $S_\epsilon(t, x) = (t/\epsilon^2, x/\epsilon)$.

## 2.2. Justification

We claim that the assumption on the noise is rather weak on the ground that many natural constructions yield stationary random processes that satisfy it. We provide details for the following example:

**Proposition 2.3.** Let $\theta: \mathbf{R}^{d+1} \to \mathbf{R}$ be smooth away from 0 and such that for all $k \in \mathbf{Z}_+^d$, $|D^k\theta(z)| \lesssim \|z\|^{-2\bar{c}-|k|}$ for $\|z\| > 1$ and $|D^k\theta(z)| \lesssim \|z\|^{-\epsilon-|k|}$ for $\|z\| \le 1$. Let $\mu$ be a Poisson point measure over $\mathbf{R}^{d+1}$ with intensity 1 and set $\eta = \mu \star \theta$, then $\eta$ satisfies the above assumption.

Before proving this proposition, let us first establish a property of joint cumulants of integrals of deterministic functions with respect to a Poisson point measure.

**Lemma 2.4.** Let $p \ge 1$ and let $f_1, \dots, f_p$ be elements of $L^1(\mathbf{R}^{d+1}) \cap L^p(\mathbf{R}^{d+1})$, and again $\mu$ be a Poisson point measure over $\mathbf{R}^{d+1}$ with intensity 1. Then the joint cumulant $\kappa_p(\mu(f_1), \dots, \mu(f_p))$ of the random variables $\mu(f_1), \dots, \mu(f_p)$ satisfies

$$ \kappa_p(\mu(f_1), \dots, \mu(f_p)) = \int_{\mathbf{R}^{d+1}} f_1(z) \times \dots \times f_p(z) \, dz. $$

*Proof.* It is sufficient to prove the result in case there exist disjoint Borel subsets $A_1, \dots, A_k$ of $\mathbf{R}^{d+1}$ with finite Lebesgue measure such that for $1 \le i \le p$,

$$ f_i(z) = \sum_{j=1}^{k} a_{i,j} \mathbf{1}_{A_j}(z). $$

In this case, however, the result follows readily from the fact that the joint cumulant is $p$-linear, and that the joint cumulant of a collection of random variables which can be split into two mutually independent subcollections vanishes, see for example property (iii) in [19, p. 32]. □

*Proof of Proposition 2.3.* It follows from Lemma 2.4 that

$$ \kappa_p(z_1, \dots, z_p) = \int_{\mathbf{R}^{d+1}} \theta(z_1 - z) \cdots \theta(z_p - z) \, dz, \quad (2.3) $$

so it remains to obtain a bound on this integral. We now consider $z_1, \dots, z_p$ to be fixed and we shall make use of the labelled tree $(T, \mathbf{n})$ built from these points as above. We are first going to treat the simpler case with all $k_i$ vanishing and then show how the argument can be modified to deal with the general case.
**Case 1.** The case with all $k_i = 0$. For every edge $e = (\varepsilon, \bar{e}) \in E_T$ and every $n \in \mathbb{Z}$ with $\mathbf{n}(\varepsilon) < n < \mathbf{n}(\bar{e})$, we define the domain

$$D_{(e,n)} = \{z \in \mathbb{R}^{d+1} : c^{-1}2^{-n_i} \le \|z - z_i\| \le c2^{-n_i}, \forall i \in \{1, \dots, p\}\},$$

where

$$n_i = \begin{cases} n & \text{if } \{i\} \ge \bar{e}, \\ \mathbf{n}(\{i\} \wedge \bar{e}) & \text{otherwise.} \end{cases}$$

It is possible to convince oneself that, provided that the constant $c$ appearing in the definition of $D_{(e,n)}$ is sufficiently large, one has

$$\bigcup_{e \in E_T} \bigcup_{\mathbf{n}(\varepsilon) < n < \mathbf{n}(\bar{e})} D_{(e,n)} = \mathbb{R}^{d+1} \setminus \{z_1, \dots, z_p\}, \quad |D_{(e,n)}| \le (2c)^{d+2} 2^{-(d+2)n}.$$

As a consequence, the integral appearing in (2.3) is bounded by some constant times

$$\sum_{e \in E_T} \sum_{\mathbf{n}(\varepsilon) < n < \mathbf{n}(\bar{e})} 2^{-(d+2)n} \prod_{i=1}^{p} \bar{\rho}(2^{-n_i}), \quad \bar{\rho}(r) = r^{-c} \wedge r^{-2\bar{c}}. \quad (2.4)$$

We first use the fact that $\bar{\rho}$ is decreasing to conclude that, for $\mathbf{n}(\varepsilon) < n < \mathbf{n}(\bar{e})$, one has the bound

$$2^{-(d+2)n} \prod_{i=1}^{p} \bar{\rho}(2^{-n_i}) \le 2^{-(d+2)n} \bar{\rho}(2^{-n}) \prod_{v \in \hat{V}_T} \bar{\rho}(2^{-\mathbf{n}(v)}). \quad (2.5)$$

This can be seen as follows. Write $\{v_1, \dots, v_k\}$ for the (possibly empty) set of nodes in $\hat{V}_T$ lying on the shortest path joining $\bar{e}$ to $\otimes$ (not including $\bar{e}$ and $\otimes$ themselves). We then have, for every $j = 1, \dots, k$,

$$\prod_{i: \{i\} \wedge \bar{e} = v_j} \bar{\rho}(2^{-n_i}) = \prod_{i: \{i\} \wedge \bar{e} = v_j} \bar{\rho}(2^{-\mathbf{n}(v_j)}) \le \prod_{v \in \hat{V}_T : v \wedge \bar{e} = v_j} \bar{\rho}(2^{-\mathbf{n}(v)}),$$

since the number of factors appearing in each term is the same. Similarly, we have

$$\prod_{i: \{i\} \ge \bar{e}} \bar{\rho}(2^{-n}) \le \bar{\rho}(2^{-n}) \prod_{v \in \hat{V}_T : v \ge \bar{e}} \bar{\rho}(2^{-\mathbf{n}(v)}),$$

hence (2.5). Since $\bar{\rho} \le \rho \wedge \rho^2$, it follows from (2.5) that

$$2^{-(d+2)n} \prod_{i=1}^{p} \bar{\rho}(2^{-n_i}) \le 2^{-(d+2)n} \bar{\rho}(2^{-n}) \rho(2^{-\mathbf{n}(\Omega_p)}) \prod_{v \in \hat{V}_T} \rho(2^{-\mathbf{n}(v)}).$$

It remains to observe that $\sum_{n \in \mathbb{Z}} 2^{-(d+2)n} \bar{\rho}(2^{-n}) = \sum_{n \in \mathbb{Z}} 2^{(c-d-2)n} \wedge 2^{(2\bar{c}-d-2)n} < \infty$, so that (2.4) is indeed bounded by the required expression.
**Case 2.** Note that we actually showed that the expression (2.3) with $\theta$ replaced by $\bar{\rho}$ is bounded by the right hand side of (2.1) with $k_i = 0$. To obtain the general case, it therefore suffices to show that

$$ \int_{\mathbf{R}^{d+1}} \prod_{i=1}^{d} D^{k_i} \theta(z_i - z) dz \lesssim \prod_{i=1}^{p} 2^{|k_i|} \mathbf{n}(i^\uparrow) \int_{\mathbf{R}^{d+1}} \prod_{i=1}^{d} \bar{\rho}(z_i - z) dz. \quad (2.6) $$

Let $\chi \in C^\infty(\mathbf{R}^{d+2})$ be such that

$$ \chi(z) = \begin{cases} 1, & \text{on } B(0, 1/4); \\ 0, & \text{on } B(0, 1/2)^c. \end{cases} $$

For $1 \le j \le p$, we define $\chi_j(z) = \chi(2^{\mathbf{n}(j^\uparrow)}(z-z_j))$, and $\chi_0(z) = 1 - \sum_{j=1}^p \chi_j(z)$. It is clear that

$$ \begin{aligned} & \int_{\mathbf{R}^{d+1}} \prod_{i=1}^{d} D^{k_i} \theta(z_i - z) dz \\ &= \sum_{j=0}^{p} \int_{\mathbf{R}^{d+1}} \chi_j(z) \prod_{i=1}^{d} D^{k_i} \theta(z_i - z) dz \\ &= \int_{\mathbf{R}^{d+1}} \chi_0(z) \prod_{i=1}^{d} D^{k_i} \theta(z_i - z) dz \\ &\quad + \sum_{j=1}^{p} \int_{\mathbf{R}^{d+1}} \theta(z_j - z) D^{k_j} \left( \chi_j(\cdot) \prod_{i \ne j} D^{k_i} \theta(z_i - \cdot) \right) (z) dz. \end{aligned} $$

We note that for $z$ in the support of $\chi_0$, for $1 \le i \le p$, $2\|z_i - z\| \ge 2^{-\mathbf{n}(i^\uparrow)}$,

$$ \begin{aligned} |D^{k_i} \theta(z_i - z)| &\lesssim \bar{\rho}(\|z_i - z\|) \cdot \|z_i - z\|^{-|k_i|} \\ &\lesssim \bar{\rho}(\|z_i - z\|) \cdot 2^{|k_i|} \mathbf{n}(i^\uparrow), \end{aligned} $$

thus yielding (2.6) as required.

To bound the final term, we note that its integrand can be written as a finite sum of terms of the form

$$ M(z) = \theta(z_j - z) D^k \chi_j(z) \prod_{i \ne j} D^{k_i + k_{j,i}} \theta(z_i - z), $$

where $k, k_{j,i} \in \mathbf{Z}_+^{d+1}$ and $k + \sum_{i \ne j} k_{j,i} = k_j$. Each of these terms is bounded above by the indicator function of the support of $\chi_j$ times

$$ \bar{\rho}(\|z - z_j\|) 2^{|k|} \mathbf{n}(j^\uparrow) \prod_{i \ne j} \bar{\rho}(\|z - z_i\|) \|z - z_i\|^{-|k_i| - |k_{j,i}|}. $$

Since for $z$ in the support of $\chi_j$ and $i \ne j$, $2\|z - z_j\| \le \|z_i - z_j\|$, so that

$$ \|z_i - z_j\| \le \|z_i - z\| + \|z - z_j\| \le \|z_i - z\| + \frac{1}{2} \|z_j - z_i\|, $$
one has $2\|z - z_i\| \ge \|z_j - z_i\| \ge 2^{-n(i^\uparrow)} \wedge 2^{-n(j^\uparrow)}$. Combining all of these bounds, we finally obtain

$$ 
\begin{aligned}
|M(z)| &\le \bar{\rho}(\|z - z_j\|) 2^{|k|\|n(j^\uparrow)\|} \prod_{i \ne j} \bar{\rho}(\|z - z_i\|) 2^{-|k_i|\|n(i^\uparrow)\|} 2^{-|k_{j,i}|\|n(j^\uparrow)\|} \\
&\le \prod_i \bar{\rho}(\|z - z_i\|) 2^{|k_i|\|n(i^\uparrow)\|},
\end{aligned}
$$

at which point we apply again (2.6) to obtain the required bound. □

### 3. Law of Large Numbers

The aim of this section is to use a simplified variant of the arguments in [14] to show that Theorem 1.1 holds.<fnref n="2" /> Although it would probably not be much more involved to obtain this proof by usual techniques, we give a proof using regularity structures. The main reason is that this allows us to introduce in a simpler setting a number of notions and notations that will be useful in the proof of our main result later on.

Before we turn to the proof proper, let us comment on the way in which we deal with the Neumann boundary conditions. Writing $P_{\text{Neu}}$ for the Neumann heat kernel and using the notation $z = (t, x)$ (and similarly for $z'$), we rewrite (1.1) as an integral equation

$$ u_{\varepsilon}(z) = \int_{D_t} P_{\text{Neu}}(z, z') (H(u_{\varepsilon}(z')) + G(u_{\varepsilon}(z')) \eta_{\varepsilon}(z')) dz' \\ + \int_D P_t^{\text{Neu}}(x, x') u_0(x') dx', $$

where $D_t = [0, t] \times D$ and $P_t^{\text{Neu}}(x, x') = P_{\text{Neu}}((t, x), (0, x'))$. We also fix an arbitrary time horizon $T \le 1$ which is not a restriction since the argument can be iterated.

Following [9], we then construct two functions $K$ on $\mathbf{R}^{d+1}$ and $K_\partial$ on $\mathbf{R}^{d+1} \times \mathbf{R}^{d+1}$ such that $K$ is compactly supported, $K_\partial$ is supported on a strip of finite width around the diagonal, and the identity

$$ P_{\text{Neu}}(z, z') = K_{\text{Neu}}(z, z'), $$

holds for $z, z' \in [0, 1] \times D$, where we set

$$ K_{\text{Neu}}(z, z') = K(z - z') + K_{\partial}(z, z'). \quad (3.1) $$

See “Appendix A” for more details on the construction of these kernels and a proof that this can be done in a way that is compatible with the results of [9, 11] that we will use in our argument.

<footnote n="2">² Simplified except for the treatment of the boundary conditions, which leads to non-trivial complications.</footnote>
**Remark 3.1.** We make no claim on the values of $K$ and $K_\partial$ for arguments outside of $[0, 1] \times D$. This is because these will always be integrated against functions that are supported on $[0, 1] \times D$ and only the values of the result inside the domain will matter.

We choose $K$ in such a way that it coincides with the whole space heat kernel $P$ on the (parabolic) ball of radius 1 and is compactly supported in the ball of radius 2. We furthermore choose $K$ in such a way that it annihilates polynomials of degree up to 3, is invariant under the transformation $(t, x) \mapsto (t, -x)$, and is such that the sum of its reflections agrees with the Neumann heat kernel on $[0, 1] \times D$. (See “Appendix A” for more details.) For example, we can choose $K$ as in [14]. The kernel $K_\partial$ is a correction term that encodes the effect of the boundary condition. Regarding our regularity structure, we then proceed as if there was no boundary condition whatsoever: we construct models defined on the whole space that are translation invariant and we use convolution with $K$ as our integration operator. We then define an operator $P_{\text{Neu}}$ on modelled distributions by setting

$$ P_{\text{Neu}} = \mathcal{K} + \tilde{\mathcal{K}}_{\partial}, \quad \text{where} \quad \tilde{\mathcal{K}}_{\partial} = \mathcal{L}_2 K_{\partial} \mathcal{R}, \qquad (3.2) $$

and $\mathcal{K}$ is built in exactly the same way as in [11, Sec. 4]. Note that $\tilde{\mathcal{K}}_\partial$ encodes the effect of the boundary condition. (There is a completely analogous definition in the case of Dirichlet boundary conditions.)

Here, $\mathcal{L}_\gamma: C^\gamma \to D^\gamma$ denotes the “Taylor lift” given by

$$ (\mathcal{L}_\gamma f)(z) = \sum_{|k| \le \gamma} \frac{f^{(k)}(z)}{k!} X^k, \qquad (3.3) $$

where $z = (t, x)$ and $k$ denotes a multiindex in $\mathbb{N}^{1+d}$.

We now have the preliminaries in place to turn to the proof of Theorem 1.1.

**Proof of Theorem 1.1.** We use a strategy similar to that in [14, 15], combining this with results from [9] to deal with the boundary conditions. We refer to [7, 8, 12] for introductions to the theory of regularity structures, as well as to [11] for details. In our present context, we use the regularity structure obtained by extending the usual polynomial structure with parabolic scaling with a symbol $\Xi_1$ of degree $-1 - \kappa$ representing the driving noise $\eta_\epsilon$, as well as an abstract integration operator $\mathcal{I}$ of order 2 representing convolution with the (singular part of) heat kernel. As usual, we will often use graphical representations for the basis vectors in our regularity structure(s), and we decree that $\circ$ is our symbolic representation for $\Xi_1$. (The reason for introducing the “accent” representing the index “1” will become clear later on where more general notations of this type are needed.) Although our goal is to consider (1.1) on the bounded domain $D \subset \mathbb{R}^d$, we construct the models for our regularity structure on the whole of $\mathbb{R} \times \mathbb{R}^d$.

With notations almost identical to those in [14] and the formula (3.19) there, it would then be natural to consider a fixed point problem of the type

$$ U = P_{\text{Neu}} \mathbf{1}_{+}^{D} (\hat{H}_{\eta}(U) + \hat{G}(U) \circ) + P_{\text{Neu}} u_0, \qquad (3.4) $$
where $\mathbf{1}_+^D$ denotes the indicator function of the space-time domain $\mathbf{R}_+ \times D$. Leaving considerations regarding the precise spaces of modelled distributions in which this equation makes sense aside for the moment, it is straightforward to see as in [11] that if we solve (3.4) for the renormalised lift of $\eta_\epsilon$, that is the admissible model such that

$$ \hat{\Pi}_{\epsilon} \circ, = \eta_{\epsilon}, \quad \hat{\Pi}_{\epsilon} \circled_{\epsilon} = \eta_{\epsilon}(K \star \eta_{\epsilon}) - \circled_{\epsilon}, \qquad (3.5) $$

then the function $u_\epsilon = \mathcal{R}_\epsilon U$ actually solves (1.1).

Indeed, iterating (3.4), we see that any solution $U$ to such a fixed point problem is necessarily of the form

$$ U = u\mathbf{1} + G(u) \circled_{\epsilon} + \nabla u X, $$

for some continuous functions $u$ and $\nabla u$. (This is purely notational, $\nabla u$ is not the gradient of $u$, but can be interpreted as a kind of “renormalised gradient”.) In particular, the factor multiplying $\mathcal{P}_{\text{Neu}}\mathbf{1}_+^D$ in the right hand side of (3.4) is given by

$$ L \stackrel{\text{def}}{=} H_{\eta}(u)\mathbf{1} + G(u) \circ, + G'(u)G(u) \circled_{\epsilon} + G'(u)u' \otimes, \qquad (3.6) $$

where we projected onto terms of negative (or vanishing) degree. At this point, it then suffices to note that the application of the reconstruction operator to $L$ yields

$$ \begin{aligned} (\mathcal{R}_{\epsilon} L)(z) &= (\hat{\Pi}_{\epsilon}^{\circ} L(z))(z) = H_{\eta}(u(z)) + G(u(z))\eta_{\epsilon}(z) - (G'G)(u(z))\circled_{\epsilon} \\ &= H(u(z)) + G(u(z))\eta_{\epsilon}(z), \end{aligned} $$

as required.

The problem with the argument outlined above is that since $\deg \circ, < -1$, the behaviour of the modelled distribution $\mathbf{1}_+^D L$ near $\partial D$ is such that the reconstruction operator is not a priori well-defined on it, see [9, Secs 4.1 & 4.2]. This is for precisely the same reason why the restriction of a generic distribution $\zeta \in C^\alpha$ to a “nice” domain $D$ is only well-defined if $\alpha > -1$. (For $\alpha \le -1$ there are non-zero distributions with support contained in $\partial D$.)

Before we tackle this problem, recall the definition of the spaces $D^{\gamma, \eta}$ as in [11, Sec. 6] (the hyperplane $P$ being given by the time slice $t=0$) as well as the spaces $D^{\gamma, w}$ as in [9, Sec. 4] (in which case $P_0$ is again the time 0 slice while $P_1 = \mathbf{R} \times \partial D$). These two spaces are distinguished by the fact that $\eta$ is a real exponent while $w$ denotes a triple of exponents describing the singular behaviour near $t=0$, $\partial D$ and the intersection of both regions respectively. It will be convenient to use the notation $(\alpha)_3$ with $\alpha \in \mathbf{R}$ as a shorthand for the triple $(\alpha, \alpha, \alpha)$.

The idea then is the following. First, we introduce a new symbol $\diamond$, also of degree $-1-\kappa$, but representing the function $\eta_\epsilon \mathbf{1}_{\mathbf{R} \times D}$ instead of representing $\eta_\epsilon$ and we add to our regularity structure the symbols $X \diamond$, and $\circled_{\epsilon}$. Write then $V$ for the sector spanned by $\circ$, $\circled_{\epsilon}$, and $X \circ$, $\hat{V}$ for the sector spanned by $\diamond$, $\circled_{\epsilon}$, and $X \diamond$, and $\iota: V \to \hat{V}$ for the linear map with $\iota \circ, = \diamond$, and similarly for the remaining basis vectors. We will furthermore only ever consider models $\hat{\Pi}$ with the property that

$$ (\hat{\Pi}_{\iota \tau})(\phi) = \begin{cases} 0 & \text{if } \text{supp}\phi \subset \mathbf{R} \times D^c, \\ (\hat{\Pi}_{\tau})(\phi) & \text{if } \text{supp}\phi \subset \mathbf{R} \times D, \end{cases} \qquad (3.7) $$
for all $\tau \in V$. Since $\iota$ commutes with the structure group and preserves degrees, it follows that $F \mapsto \iota F$ is continuous from $D^{\gamma,w}(V)$ to $D^{\gamma,w}(\hat{V})$ for all choices of exponents $\gamma$ and $w$ and, for $\gamma > 0$, the local reconstruction operator $\tilde{R}$ (which yields a distribution on $\mathbf{R} \times (\mathbf{R}^d \setminus \partial D)$ and is always well-defined) satisfies

$$ (\tilde{R}_{\iota F})(\phi) = \begin{cases} 0 & \text{if } \text{supp}\phi \subset \mathbf{R} \times D^c, \\ (\tilde{R}F)(\phi) & \text{if } \text{supp}\phi \subset \mathbf{R} \times D. \end{cases} $$

The reason for the introduction of these extra symbols is that we would like to interpret (3.4) as a fixed point problem in the space $D^{1+2\kappa,(2\kappa)_3}$ with values in $V_0$ (for small enough $\kappa$), where $V_0$ is spanned by $\mathbf{1}, \mathbf{0}',$ and $X$. The problem now is that, for $F \in D^{1+2\kappa,(2\kappa)_3}$, we have $F \circ_\iota \in D^{\kappa,(\kappa-1)_3}$, but we lose some regularity at the boundary when multiplying by the indicator function of our domain (see [9]), so that we only have $\mathbf{1}_+^D F \circ_\iota \in D^{\kappa,(\kappa-1)_3}$. Since the boundary index is now below $-1$, it follows that the reconstruction operator of [9] is not well-defined on $\mathbf{1}_+^D F \circ_\iota$. By Theorem C.5, it is however perfectly well-defined on $\mathbf{1}_+ F \diamond_\iota = \mathbf{1}_{+\iota}(F \circ_\iota)$ since only the temporal singularity index is below $-1$ (but above $-2$). Furthermore, one has the identity $\mathcal{R}(\mathbf{1}_+ F \diamond_\iota) = \tilde{R}(\mathbf{1}_+^D F \circ_\iota)$ for test functions whose support does not intersect the boundary of $D$.

Recall now that, given a modelled distribution $F$ and a distribution $\zeta$ agreeing with $\tilde{R}F$ outside the boundary of $D$, [9, Lem 4.12] defines a modelled distribution $\mathcal{K}^\zeta F$ with improved regularity and such that $\mathcal{R}\mathcal{K}^\zeta F = K \star \zeta$. Furthermore, the map $(\zeta, F) \mapsto \mathcal{K}^\zeta F$ is Lipschitz continuous in the natural topologies. This allows us to define an operator $\mathcal{P}_D: D^{\kappa,(\kappa-1)_3}(V) \to D^{\kappa+2,(1-\kappa)_3}(V)$ by setting

$$ \mathcal{P}_D: F \mapsto \mathcal{K}^{\tilde{R}\iota F} \mathbf{1}_+^D F + \tilde{K}_{\partial\iota} \mathbf{1}_{+\iota} F. $$

Our discussion suggests that, instead of (3.4), we should consider the fixed point problem

$$ U = \mathcal{P}_{\text{Neu}} \mathbf{1}_+^D \hat{H}_\eta(U) + \mathcal{P}_D \hat{G}(U) \circ_\iota + \mathcal{P}_{\text{Neu}} u_0, \quad (3.8) $$

which admits unique local solutions in $D^{1+2\kappa,(2\kappa)_3}(V)$ by [9], which are continuous with respect to admissible models on the full regularity structure satisfying furthermore the consistency condition (3.7). Here, we need to choose $\kappa$ small enough to guarantee that $P_{\text{Neu}} u_0$ does indeed belong to the space $C^{\kappa+2,(1-\kappa)_3}$, which is possible thanks to our assumption that $u_0$ itself is Hölder continuous for some positive exponent.

Retracing the discussion given at the beginning of the proof, but now with the renormalised model $\hat{\mathbf{\Pi}}_\epsilon$ such that, in addition to (3.5), one has $\hat{\mathbf{\Pi}}_\epsilon \iota \tau = \mathbf{1}_{\mathbf{R} \times D} \hat{\mathbf{\Pi}}_\epsilon \tau$ for $\tau \in V$, we conclude that for $\epsilon > 0$, solutions to (3.8) coincide with those of (1.1). We now refer to Theorem 4.8 which shows that the sequence of models $\hat{\mathbf{\Pi}}_\epsilon$ converges, as $\epsilon \to 0$, to a limiting model $\hat{\mathbf{\Pi}}$ such that $\hat{\mathbf{\Pi}} \circ_\iota = \hat{\mathbf{\Pi}} \diamond_\iota = 0$, extended canonically to the whole regularity structure. It follows immediately that the solution $\bar{U}$ to (3.4) with the model $\hat{\mathbf{\Pi}}$ is such that $u^{(0)} = \mathcal{R}\bar{U}$ does indeed solve (1.5) as claimed. $\square$
Remark 3.2. Note that (3.7) does not force us to set $\hat{\Pi}_{\varepsilon} \diamondsuit_\iota = \mathbf{1}_{\mathbf{R}\times D}\hat{\Pi}_{\varepsilon} \circ_\iota = \mathbf{1}_{\mathbf{R}\times D}\eta_{\varepsilon}$, but we could have added a sufficiently regular distribution supported on $\partial D$. This however would break the identity

$$ \mathcal{R}_{\varepsilon} \mathcal{P}_{D} G = P_{\text{Neu}} \mathbf{1}_{+}^{D} \mathcal{R}_{\varepsilon} G, $$

on $[0, 1] \times D$ and would therefore modify the boundary condition of the resulting solution.

## 4. Central Limit Theorem

We now turn to the proof of the main result of this article, Theorem 1.2. We will mainly focus on the case of Neumann boundary conditions in dimension $d = 3$, which is the most interesting (and technically most difficult) case. We set

$$ v_{\varepsilon} = \frac{u_{\varepsilon} - u^{(0)} - \varepsilon u^{(1)}}{\varepsilon^{d/2}}, \quad \xi_{\varepsilon} = \varepsilon^{-d/2} \eta_{\varepsilon}, \quad \sigma_{\varepsilon} = \varepsilon^{d/2} \eta_{\varepsilon}. \quad (4.1) $$

With this notation, we then have in the case $d \in \{2, 3\}$

$$
\begin{align*}
\partial_t v_{\varepsilon} &= \Delta v_{\varepsilon} + H'_{\eta}(u^{(0)})v_{\varepsilon} + \varepsilon^{-\frac{d}{2}}(H(u_{\varepsilon}) - H_{\eta}(u^{(0)}) - H'_{\eta}(u^{(0)})(u_{\varepsilon} - u^{(0)})) \\
&\quad + G(u_{\varepsilon})\xi_{\varepsilon} - \varepsilon^{1-\frac{d}{2}}\Psi(u^{(0)}, \nabla u^{(0)}) \\
&= \Delta v_{\varepsilon} + H'_{\eta}(u^{(0)})v_{\varepsilon} + G(u^{(0)})\xi_{\varepsilon} \\
&\quad + \varepsilon^{-\frac{d}{2}}(H_{\eta}(u_{\varepsilon}) - H_{\eta}(u^{(0)}) - H'_{\eta}(u^{(0)})(u_{\varepsilon} - u^{(0)})) - \varepsilon^{-\frac{d}{2}}\mathring{\circ}(GG')(u^{(0)}) \\
&\quad - \varepsilon^{-\frac{d}{2}}\mathring{\circ}((G')^2 + GG'')(u^{(0)})(u_{\varepsilon} - u^{(0)}) + G'(u^{(0)})(u_{\varepsilon} - u^{(0)})\xi_{\varepsilon} \\
&\quad + \frac{1}{2}G''(u^{(0)})(u_{\varepsilon} - u^{(0)})^2\xi_{\varepsilon} - \varepsilon^{1-\frac{d}{2}}\Psi(u^{(0)}, \nabla u^{(0)}) + R_{1,\varepsilon},
\end{align*}
$$

where, setting

$$ w_{\varepsilon} = \varepsilon^{-\frac{d}{2}}(u_{\varepsilon} - u^{(0)}) = v_{\varepsilon} + \varepsilon^{1-\frac{d}{2}}u^{(1)}, $$

we have the explicit expression for the remainder term

$$
\begin{align*}
R_{1,\varepsilon} = -\varepsilon^{d/2} w_{\varepsilon}^2 &\int_0^1 \mathring{\circ}(GG')''(u^{(0)} + s\varepsilon^{d/2}w_{\varepsilon})(1-s) \, ds \\
&\quad + \frac{\varepsilon^{\frac{3d}{2}}}{2} w_{\varepsilon}^3 \int_0^1 G^{(3)}(u^{(0)} + s\varepsilon^{d/2}w_{\varepsilon})(1-s)^2 \, ds \, \xi_{\varepsilon}.
\end{align*}
$$

Furthermore, due to the non-vanishing boundary condition of $u^{(1)}$ in the Neumann case, $v_{\varepsilon}$ is then endowed with the inhomogeneous boundary condition

$$ \langle \nabla v_{\varepsilon}, n(x) \rangle = -\varepsilon^{1-d/2} c(x) GG'(u^{(0)}(t, x)). $$

We now also incorporate the first part of the second line into the remainder, so that we can write

$$ \partial_t v_{\varepsilon} = \Delta v_{\varepsilon} + H'_{\eta}(u^{(0)})v_{\varepsilon} + G(u^{(0)})\xi_{\varepsilon} - \varepsilon^{-\frac{d}{2}}\mathring{\circ}(GG')(u^{(0)}) $$
$$ - \stackrel{\circ}{\circ} (GG')'(u^{(0)})(v_{\varepsilon} + \varepsilon^{1-\frac{d}{2}}u^{(1)}) + G'(u^{(0)})(v_{\varepsilon} + \varepsilon^{1-\frac{d}{2}}u^{(1)})\eta_{\varepsilon} \\ + \frac{1}{2}G''(u^{(0)})(v_{\varepsilon} + \varepsilon^{1-\frac{d}{2}}u^{(1)})^2\sigma_{\varepsilon} - \varepsilon^{1-\frac{d}{2}}\Psi(u^{(0)}, \nabla u^{(0)}) + \hat{R}_{\varepsilon}^{(d)}(v_{\varepsilon}, \varepsilon^{\alpha}\xi_{\varepsilon}), \quad (4.2) $$

where the remainder term $\hat{R}_{\varepsilon}^{(d)}$ is given by

$$ \begin{aligned} \hat{R}_{\varepsilon}^{(d)}(v, \zeta) &= R_{\varepsilon}^{(d)}(v + \varepsilon^{1-d/2}u^{(1)}, \zeta), \\ R_{\varepsilon}^{(d)}(w, \zeta) &= \varepsilon^{d/2}w^2 \int_0^1 H''(u^{(0)} + s\varepsilon^{d/2}w)(1-s) \, ds \\ &\quad + \frac{\varepsilon^{\frac{3d}{2}-\alpha}}{2} w^3 \int_0^1 G^{(3)}(u^{(0)} + s\varepsilon^{d/2}w)(1-s)^2 \, ds \zeta. \end{aligned} \quad (4.3) $$

Note that here we have $H$ appearing in (4.3) rather than $H_{\eta}$ and that the two are related by (1.4). In dimension $d=1$, we set $u^{(1)} = 0$ so that $\hat{R}_{\varepsilon}^{(d)} = R_{\varepsilon}^{(d)}$ and $v_{\varepsilon} = (u_{\varepsilon} - u^{(0)})/\sqrt{\varepsilon}$, and we obtain in the same way the slightly simpler expression

$$ \begin{aligned} \partial_t v_{\varepsilon} ={}& \Delta v_{\varepsilon} + H'_{\eta}(u^{(0)})v_{\varepsilon} + G(u^{(0)})\xi_{\varepsilon} - \varepsilon^{-\frac{1}{2}} \stackrel{\circ}{\circ} (GG')(u^{(0)}) \\ & - \stackrel{\circ}{\circ} (GG')'(u^{(0)})v_{\varepsilon} + G'(u^{(0)})v_{\varepsilon}\eta_{\varepsilon} + \frac{1}{2}G''(u^{(0)})v_{\varepsilon}^2\sigma_{\varepsilon} + \hat{R}_{\varepsilon}^{(1)}(v_{\varepsilon}, \varepsilon^{\alpha}\xi_{\varepsilon}) \end{aligned} \quad (4.4) $$

(The reason why the term containing $\Psi$ does not appear in this expression is because this was generated by $\partial_t u^{(1)}$ which vanishes by definition in dimension one.) The exponent $\alpha$ appearing in this expression is of course arbitrary, but allowing to tune it will be convenient when expressing this as a fixed point problem.

### 4.1. Decomposition of the Solution

In order to show that $v_{\varepsilon}$ converges to a limit, it will be convenient to break it into a sum of three terms. The first term will be a straightforward approximation to the stochastic heat equation with noise strength $G(u^{(0)})$ and homogeneous boundary condition. The second term will converge to 0, but incorporates the diverging boundary condition, which is used to compensate a resonance appearing in its right hand side. The final term will be a remainder that is sufficiently regular to be dealt with by the techniques of [9]. For this, we write

$$ v_{\varepsilon} = v_{\varepsilon}^{(0)} + v_{\varepsilon}^{(1)} + \bar{v}_{\varepsilon}, $$

and, with the convention that $G$ and its derivatives are always evaluated at $u^{(0)}$, we set

$$ \partial_t v_{\varepsilon}^{(0)} = \Delta v_{\varepsilon}^{(0)} + G\xi_{\varepsilon}, \quad (4.5) $$

$$ \partial_t v_{\varepsilon}^{(1)} = \Delta v_{\varepsilon}^{(1)} + G'\left(v_{\varepsilon}^{(0)}\eta_{\varepsilon} - \varepsilon^{-\frac{d}{2}}\stackrel{\circ}{\circ} G - \varepsilon^{1-\frac{d}{2}}G'\left(\stackrel{\circ}{\circ}, \nabla u^{(0)}\right)\right) $$
$$ + \frac{1}{2} G''((v_{\varepsilon}^{(0)})^2 \sigma_{\varepsilon} - 2 \stackrel{\circ}{\varepsilon} G v_{\varepsilon}^{(0)} - \varepsilon^{1-\frac{d}{2}} \stackrel{\circ}{\varepsilon} \stackrel{\circ}{\varepsilon} G^2), \quad (4.6) $$

endowed with the boundary conditions on $\partial D$

$$ \langle \nabla v_{\varepsilon}^{(0)}, n \rangle = 0, \quad \langle \nabla v_{\varepsilon}^{(1)}, n \rangle = -\varepsilon^{1-\frac{d}{2}} c GG', $$

as well as vanishing initial conditions. The reason why $v_{\varepsilon}^{(1)}$ will actually converge to 0 despite the diverging boundary condition when $d \ge 2$ is the following. Consider the function $\Pi_{\varepsilon} \upharpoonright$ defined by

$$ (\Pi_{\varepsilon} \upharpoonright)(z) = \int_{\mathbb{R} \times D} K_{\text{Neu}}(z, z') \xi_{\varepsilon}(z') dz' - \int_{\mathbb{R}^{d+1}} K(z - z') \xi_{\varepsilon}(z') dz', \quad (4.7) $$

Then, we will see in (4.24) that the behaviour of $v_{\varepsilon}^{(0)}$ is locally very well described by that of

$$ G(u^{(0)})(\Pi_{\varepsilon} \upharpoonright + \Pi_{\varepsilon} \downharpoonright), \quad (4.8) $$

where $\Pi_{\varepsilon} \upharpoonright = K \star \xi_{\varepsilon}$ as usual. This implies that the behaviour of the term $G'(u^{(0)})v_{\varepsilon}^{(0)}\eta_{\varepsilon}$ appearing in the right hand side of the equation for $v_{\varepsilon}^{(1)}$ is well described locally by that of

$$ (GG')(u^{(0)})(\Pi_{\varepsilon} \upharpoonright + \Pi_{\varepsilon} \downharpoonright), \quad (4.9) $$

Where $\circ$ denotes multiplication by $\eta_{\varepsilon}$ as previously, while $\diamond$, denotes multiplication by $\eta_{\varepsilon} \mathbb{1}_{\mathbb{R} \times D}$. (This will be formalised later on.) We will see that, up to vanishingly small errors, $\Pi_{\varepsilon} \upharpoonright \approx \varepsilon^{-\frac{d}{2}} \stackrel{\circ}{\varepsilon}$, while $\Pi_{\varepsilon} \downharpoonright \approx \varepsilon^{1-\frac{d}{2}} c \delta_{\partial D}$ for a suitable constant $c$ (in fact a different constant for each face of $\partial D$ in general), so that the first term in (4.9) is cancelled up to small errors by the term $-\varepsilon^{-\frac{d}{2}} \stackrel{\circ}{\varepsilon} (GG')(u^{(0)})$ appearing in the equation for $v_{\varepsilon}^{(1)}$, while the second term in (4.9) is cancelled by the term $-\varepsilon^{1-\frac{d}{2}} c GG' \delta_{\partial D}$ created by the boundary condition.

Note that this argument does not see much of a difference between the Neumann and Dirichlet cases. Indeed, if we want the right hand side of the equation for $v_{\varepsilon}^{(1)}$ to converge to a limiting distribution for the latter, we also need to add a diverging (when $d = 3$) term proportional to $\delta_{\partial D}$. The difference is that $K_{\text{Dir}}(z, z')$ vanishes for $z' \in \partial D$, so that this term has no influence on $v_{\varepsilon}^{(1)}$ in the Dirichlet case.

The idea now is to proceed as follows:

* In a first step, we describe in Section 4.2 a regularity structure that is sufficiently rich to allow us to give precise control on the behaviour of $v_{\varepsilon}^{(0)}$, $v_{\varepsilon}^{(1)}$ and $v_{\varepsilon}$. As already alluded to, this will in particular include symbols representing non-stationary space-time stochastic processes, but we will try to keep these to an absolute minimum. We also describe there the renormalisation procedure which allows us to construct a suitable (random) model.
* We then make precise the description (4.8) for $v_{\varepsilon}^{(0)}$ by expressing it as a modelled distribution in this regularity structure, which in particular contains the two symbols ⌢ and ⌣. The main conclusion in (4.24) will be that the presence of ⌣ allows us to express $v_{\varepsilon}^{(0)}$ as a modelled distribution with a good behaviour near $\partial D$. (At this stage, we could of course also use one symbol only and define our model using the Neumann heat kernel only, but the decompositions (4.8) and (4.9) are convenient for the remainder of the argument and to be able to reuse existing results.)

* In a second step, we show that if we set

$$ \hat{\Pi}_{\varepsilon} \circled{=} \Pi_{\varepsilon} \circled{=} -\varepsilon^{-\frac{d}{2}} \circled{,} \quad \hat{\Pi}_{\varepsilon} \text{bold} \circled{=} \Pi_{\varepsilon} \text{bold} \circled{=} -\varepsilon^{1-\frac{d}{2}} c \delta_{\partial D}, $$

then $\hat{\Pi}_{\varepsilon}$ converges to a limiting model as $\varepsilon \to 0$, which furthermore has good restriction properties to $D$, uniformly in $\varepsilon$. Furthermore, since $\hat{\Pi}_{\varepsilon}$ is only singular near the boundary of $D$, we can describe $v_{\varepsilon}^{(0)}$ by another modelled distribution with worse behaviour near the boundary, but which only uses “translation invariant symbols” in its description, see Lemma 4.10, which allows us to give a description of $v_{\varepsilon}^{(1)}$ in terms of such symbols in Lemma 4.12.

* We set up a sufficiently large regularity structure so that we can formulate a fixed point problem for the remainder $\tilde{v}_{\varepsilon}$ and control its behaviour as $\varepsilon \to 0$, see Propositions 4.16 and 4.18. Combining this with the convergence of the corresponding renormalised model which is performed in Theorem 4.8 but relies crucially on the next section, we are finally able to conclude.

## 4.2. Definition of the Ambient Regularity Structure

We start by defining a regularity structure that is sufficiently large to allow us to perform the steps mentioned above and in particular to formulate (4.2) as a fixed point problem for a modelled distribution $V$. This fixed point problem will be chosen precisely in such a way that if we take as our model the renormalised lift of the noise $\eta_{\varepsilon}$, then the corresponding counterterms are precisely such that if $V$ solves the fixed point problem, then $RV$ solves (4.2).

Let $\Xi_i^j$ be new symbols representing $\xi_{\varepsilon}\varepsilon^{\alpha(i,j)}$ with

$$ \alpha(i, j) \stackrel{\text{def}}{=} \left(1 - \frac{d}{2}\right) j + \frac{d}{2} i, \quad 0 \le j \le i \le 2. \qquad (4.10) $$

(Beware that $j$ is simply a superscript in $\Xi_i^j$, not a power.) In the graphical notation analogous to [14], we will use “accents” to denote the upper and lower indices on $\Xi = \Xi_0^0 = \circ$, so for example $\Xi_2 = \circ_{\mu}$, $\Xi_1^1 = \circ_{\nu}$, etc. We will also sometimes write $\circ_i^j$ instead of $\Xi_i^j$. The degree of these symbols is given by

$$ \deg \Xi_i^j = \frac{\delta - 1}{2} \wedge \left( \alpha(i, j) - \frac{d+2}{2} - \kappa \right). \qquad (4.11) $$

The main reason why we never consider these “noises” as having degree larger than $\frac{\delta-1}{2}$ is that we want to view them as “noise types” and so the structure group
should act trivially on them. A byproduct of this choice is that it allows us to deal with driving noises $\eta$ that are themselves unbounded.

We now build a regularity structure extending the one built in Section 3, in which (4.35) can be formalised, by applying the framework of [5, Sec. 5]. We take as our basic “noise types” the noises $\circled_i^j$ with $0 \le j \le i \le 2$, as well as an additional noise type $\diamond$ which will be used to represent the noise $\xi_\epsilon$, restricted to $\mathbf{R} \times \mathcal{D}$. We also introduce two “edge types” $\mathcal{I}$ representing convolution by a suitable cutoff $\mathcal{K}$ of the standard heat kernel in the whole space and $\mathcal{I}$ representing the integral operator with kernel $K_{\text{Neu}}$ as in (3.1) (see “Appendix A” for a precise definition), both having degree 2. It will be convenient to also introduce edge types $\mathcal{I}_i$ and $\mathcal{I}_i$ with $i = 1, \dots, d$ of degree 3 representing the integral operators with kernel $K_i(x, y) = (y_i - x_i)K(y-x)$ and $K_{\text{Neu},i}(x, y) = (y_i - x_i)K_{\text{Neu}}(y-x)$ respectively. Finally, we introduce a “virtual” edge types $\hat{\mathcal{I}}$ of degree 2 which will allow us to produce a rule (in the technical sense of [5, Sec. 5]) generating the relevant trees containing non-translation invariant symbols $\mathcal{I}$ and $\diamond$, but without cluttering our regularity structure with unneeded symbols.

The rule used to generate the regularity structure is then given by $R(\circled_i^j) = \{1\}$ and

$$ 
\begin{align} 
R(\mathcal{I}) &= \{1, \mathcal{I}, \mathcal{I}_i, \mathcal{I}_i \circled_i\} \cup \{\mathcal{I}^k \circled_i^j : k \le i - j\}, \\ 
R(\mathcal{I}_i) &= \{1, \circ\}, \quad R(\mathcal{I}) = R(\mathcal{I}_i) = \{1, \diamond\}, \\ 
R(\hat{\mathcal{I}}) &= \{1, \diamond, \mathcal{I} \diamond, \mathcal{I}, \mathcal{I}_i, \diamond, \diamond, \mathcal{I}, \mathcal{I} \diamond, \mathcal{I}_i \diamond\}. \tag{4.12}
\end{align} 
$$

Recall that, given a collection $\mathcal{L}$ of “edge types” (in our case, these are the “integration” types $\mathcal{I}$, $\mathcal{I}$, etc as well as the “noise” types $\circled_i^j$ and $\diamond_i^j$ which are also interpreted as edges for the purpose of this discussion), a “rule” is a map $R: \mathcal{L} \to \mathcal{P}(\hat{\mathcal{P}}(\mathcal{L}))\setminus\{\emptyset\}$, where $\mathcal{P}(\mathcal{A})$ denotes the powerset of a set $\mathcal{A}$, $\hat{\mathcal{P}}(\mathcal{L})$ denotes the set of non-empty multisets with elements in $\mathcal{L}$.$^3$ In other words, an element $R(t)$ is a collection of “node types”, with each node type being a collection of edge types, with repetitions allowed. In (4.12), we use the identification between a multiset and a formal monomial, that is $\circled_2^2$ denotes the multiset with one copy of $\circled_0$ and two copies of $\mathcal{I}$.

The basis vectors for our regularity structure are rooted trees $(V, E, \rho)$ with edges $e \in E$ labelled by $\mathcal{L}$ and nodes $v \in V$ labelled by $\mathbb{N}^d$, denoting polynomial factors, with the convention that $0 \in \mathbb{N}$. By convention, edges are oriented away from the root. Any node $v \in V$ then has a “type” $\mathcal{N}(v) \in \hat{\mathcal{P}}(\mathcal{L})$ given by the collection of the types of the outgoing edges adjacent to $v$. If $v \in V\setminus\{\rho\}$, there is a unique edge coming into $v$ and we write $t(v) \in \mathcal{L}$ for the type of the incoming edge. A tree is said to “conform” to a rule $R$ if, whenever $v \in V\setminus\{\rho\}$, one has $\mathcal{N}(v) \in R(t(v))$, while $\mathcal{N}(\rho) \in \bigcup_{t \in \mathcal{L}} R(t)$.

**Remark 4.1.** Note that $\hat{\mathcal{I}}$ never appears inside any $R(t)$, so that a tree conforming to the rule $R$ is not allowed to have any edge of type $\hat{\mathcal{I}}$. The only reason for its

<footnote>$^3$ We ignore the possibility of having “derivatives” on our edges as in [5], as these do not occur in this article.</footnote>
presence is to allow the root of a conforming tree to be of type $\mathcal{N}(\rho) \in R(\hat{\mathcal{I}})$, since otherwise the tree $\mathcal{G}$, which will be used later on would not be conforming to our rule. Actually, the symbol $\mathcal{G}$, which was introduced in a completely ad hoc manner so far, will then be interpreted as $\mathcal{G} = \mathcal{G} - \mathcal{G}$,

We assign degrees to the components appearing here by (4.11) as well as

$$ \deg \diamond = \deg \circ, $$

$$ \deg \diamond = \deg \circ, $$

$$ \deg \mathcal{I} = \deg \hat{\mathcal{I}} = \deg \mathcal{I} = 2, $$

$$ \deg \mathcal{I}_i = \deg \mathcal{I}_i = 3. $$

It is straightforward to verify that this rule is subcritical and complete. We henceforth denote by $(\mathcal{T}, \mathcal{G})$ the (reduced in the terminology of [5, Sec. 6.4]) regularity structure generated by the rule $\mathcal{R}$.

### 4.3. Description of the Models

Throughout this article, we will use the notation $\Pi$ (possibly with additional decorations) for a continuous linear map $\Pi: \mathcal{T} \to \mathcal{D}'(\mathbb{R}^{d+1})$ such that there exists a (necessarily unique) admissible model $(\Pi, \Gamma)$ related to $\Pi$ by [11, Sec. 8.3]. Here, we associate the kernel $K$ to $\mathcal{I}$ and $\hat{\mathcal{I}}$, $K_i$ to $\mathcal{I}_i$, and $K_{\text{Neu}}$, $K_{\text{Neu},i}$ to $\mathcal{I}$, $\mathcal{I}_i$ respectively.

We henceforth make a slight abuse of terminology and call $\Pi$ itself a model. The following notion of an “admissible” model is a slight strengthening of the usual one to our context which essentially states that our “square” symbols are the restrictions of the “round” symbols to $\mathbb{R} \times \mathcal{D}$.

**Definition 4.2.** A model $\Pi$ for $(\mathcal{T}, \mathcal{G})$ is *admissible* if it is admissible in the sense of [11, Def. 8.29] for the kernels listed above and furthermore, for any $\tau \in \mathcal{T}$ and $i, j$ such that $\circled_i^j \tau, \circled_i^j \tau \in \mathcal{T}$, one has

$$ (\Pi \diamond \circled_i^j \tau)(\phi) = (\Pi \circ \circled_i^j \tau)(\phi) $$

for test functions $\phi$ with $\text{supp}\phi \subset \mathbb{R} \times \mathcal{D}$ and $(\Pi \diamond \circled_i^j \tau)(\phi) = 0$ for test functions with $\text{supp}\phi \subset \mathbb{R} \times \mathcal{D}^c$.

**Remark 4.3.** Note furthermore that all basis vectors $\hat{\tau} \in \mathcal{T}$ with $\deg \hat{\tau} < -1$ are of the form $\circled_i^j \tau$ (or similar with $\circ$ replaced by $\diamond$) for some $i, j$ and $\tau$ as in Definition 4.2. As a consequence, admissible models in our sense admit a canonical decomposition $\Pi = \Pi^+ + \Pi^-$ as in Assumption C.1 by setting $\Pi^+ \circled_i^j \tau = \Pi \diamond \circled_i^j \tau$. This is a crucial fact which allows us to use Theorem C.5.

We also define the following notion:

**Definition 4.4.** Consider the rule $R_\circ$ given by $R_\circ(\circled_i^j) = R_\circ(\mathcal{I}) = R(\mathcal{I}_i) = \{\mathbf{1}\}$ and

$$ R_\circ(\mathcal{I}) = \{\mathbf{1}, \mathcal{I}, \mathcal{I}_i, \mathcal{I}_i \circ\} \cup \{\circled_i^j \mathcal{I}^k : k \le i - j\}, \\ R_\circ(\mathcal{I}_i) = \{\mathbf{1}, \circ\}, \quad R_\circ(\hat{\mathcal{I}}) = \{\mathbf{1}, \mathcal{I}\}. $$

We write $\mathcal{T}_\circ \subset \mathcal{T}$ for the sector spanned by trees conforming to $R_\circ$ and we call $\mathcal{T}_\circ$ the *translation invariant sector*.
Given any admissible model for $(\mathcal{T}, \mathcal{G})$, we then write $\mathcal{K}$, $\mathcal{K}_i$, $\mathcal{K}_{\text{Neu}}$, $\mathcal{K}_{\text{Neu},i}$ for the corresponding integration operators as defined in [11, Sec. 5]. We now provide a complete description of the renormalised models $\hat{\Pi}_\epsilon$ we will use for this regularity structure. From now on, whenever we do not explicitly mention a model on $(\mathcal{T}, \mathcal{G})$, we assume that we talk about the specific (random) model $\hat{\Pi}_\epsilon$, and not a general admissible model. In particular, all the variants of the spaces $\mathcal{D}^\nu$ used below will be the spaces associated to this model. In order to specify the $\epsilon$-dependence, we will write $\mathcal{R}_\epsilon$ for the reconstruction operator associated to the model $\hat{\Pi}_\epsilon$. Theorem 4.8 shows that $\hat{\Pi}_\epsilon$ converges to some limiting model $\Pi$ as $\epsilon \to 0$, giving rise to a reconstruction operator $\mathcal{R}_0$. Whenever we simply write $\mathcal{R}$, it denotes the reconstruction operator for a generic admissible model.

Writing $1^D$ for the indicator function of $\mathbf{R} \times \mathcal{D}$, we first define a family of non-renormalised (random) models $\Pi_\epsilon$ as the canonical lift for the “noise” given by

$$ \Pi_{\epsilon} \Xi_{i}^{j} = \xi_{\epsilon} \epsilon^{\alpha(i, j)}, \quad \Pi_{\epsilon} \diamond = 1^{D} \xi_{\epsilon}, \quad \Pi_{\epsilon} \diamond_{i} = 1^{D} \eta_{\epsilon}, $$

with $\alpha(i, j)$ as in (4.10). We then define an intermediate model $\tilde{\Pi}_\epsilon$ related to $\Pi_\epsilon$ by

$$ \begin{align} \tilde{\Pi}_{\epsilon} \circled_{i} &= \Pi_{\epsilon} \circled_{i} - \epsilon^{-\frac{d}{2}} \circled_{i} 1^{D} - \sum_{i=1}^{d} \epsilon^{1-\frac{d}{2}} (c_{i,0} \delta_{\partial_{i,0}} \mathcal{D} + c_{i,1} \delta_{\partial_{i,1}} \mathcal{D}), \nonumber \\ \tilde{\Pi}_{\epsilon} \circled_{i} &= \Pi_{\epsilon} \circled_{i} - \epsilon^{-\frac{d}{2}} \circled_{i} 1^{D}, \nonumber \\ \tilde{\Pi}_{\epsilon} \circled_{i} &= \Pi_{\epsilon} \circled_{i} - \epsilon^{-\frac{d}{2}} \circled_{i} 1^{D} \Pi_{\epsilon} X - \epsilon^{1-\frac{d}{2}} \circled_{i} 1^{D}, \tag{4.13} \end{align} $$

as well as $\tilde{\Pi}_{\epsilon} \diamond_{i} \mathcal{T}_{i}(\circled) = \tilde{\Pi}_{\epsilon}(\circled_{i} \circled_{i} \circled)$. We furthermore impose admissibility which forces us to set $\tilde{\Pi}_{\epsilon} X \tau = \Pi_{\epsilon} X \cdot \tilde{\Pi}_{\epsilon} \tau$. One can verify that the only remaining basis vector of $\mathcal{T}$ of negative degree and not belonging to the translation invariant sector is $\circled$. However, this element will never be needed for our considerations, so we do not need to specify the action of $\tilde{\Pi}_{\epsilon}$ on it.

We now construct a renormalised model $\hat{\Pi}_\epsilon$ from $\tilde{\Pi}_\epsilon$ by applying a slight modification of the BPHZ renormalisation procedure [5,6] to the translation invariant sector (which can be viewed as a regularity structure in its own right, generated by the rule $\mathcal{R}_\circ$). Writing $\mathcal{T}_-$ for the subspace of $\mathcal{T}_\circ$ consisting of symbols of strictly negative degree, we will define $\hat{\Pi}_\epsilon$ by an expression of the form

$$ \hat{\Pi}_{\epsilon} \tau = (g_{\epsilon} \otimes \tilde{\Pi}_{\epsilon}) \Delta^{-} \tau, $$

where $\Delta^{-}$ is a certain linear operator from $\mathcal{T}$ into $\text{Alg}\mathcal{T}_{-} \otimes \mathcal{T}$, with $\text{Alg}\mathcal{T}_{-}$ the free unital algebra generated by $\mathcal{T}_{-}$, and $g_\epsilon$ is a character of $\text{Alg}\mathcal{T}_{-}$, which is canonically identified with an element of the dual space $\mathcal{T}^*$. In the BPHZ renormalisation procedure, one should choose $g_\epsilon$ of the form

$$ g_{\epsilon}^{\text{BPHZ}} = (\mathbf{E}\tilde{\Pi}_{\epsilon}) \circ \mathcal{A}, \tag{4.14} $$
[In general, we should also set $C(\tau) = 0$ if $\tau$ is of the form $\tau = \mathcal{I}(\tau')$ for some $\tau'$, but the only symbol of negative degree of this type appearing in this work is $\mathring{\gamma}$ which is already covered by (4.16).] Other constants that will be relevant for our analysis (in dimension $d = 3$) are given by $\mathring{\gamma} = 0$ as well as

$$ \begin{aligned} \mathring{\gamma} &= \int P(z) P(z') \kappa_2(z, z') \mathrm{d}z \mathrm{d}z', & \mathring{\gamma} &= \int x P(z) P(z') \kappa_2(z, z') \mathrm{d}z \mathrm{d}z', \\ \mathring{\gamma} &= \int P(-z) P(-z') P(z-z'') \kappa_3(z, z', z'') \mathrm{d}z \mathrm{d}z' \mathrm{d}z'', & & \\ \mathring{\gamma} &= \int P(-z) P(-z') P(z'') \kappa_3(z, z', z'') \mathrm{d}z \mathrm{d}z' \mathrm{d}z''. & & \end{aligned} $$

The convergence of integrals corresponding to $\mathring{\gamma}$, $\mathring{\gamma}''$ and $\mathring{\gamma}'''$ in dimension 3 can easily be verified by using our assumption on the cumulants and the self-similarity of the heat kernel. The convergence of the integral for $\mathring{\gamma}$ is more subtle since $\deg \mathring{\gamma} = 0$. As a consequence, although $\kappa_2(z, z')$ decays fast enough when $\|z - z'\|$ is large, the function $z \mapsto xP^2(z)$ is homogeneous of (parabolic) degree $-5$ and is therefore not absolutely integrable at large scales. However, since it is odd under $(t, x) \mapsto (t, -x)$, additional cancellations occur and the integral should be interpreted as

$$ \mathring{\gamma} = \frac{1}{2} \int x (P(z') - P(z)) P(z) \kappa_2(z, z') \mathrm{d}z \mathrm{d}z', \quad (4.17) $$

which does converge absolutely, so we set $g_{\epsilon}(\mathring{\gamma}) = -\mathring{\gamma}$ in dimension 3.

In dimension 2, $\deg \mathring{\gamma} > 0$, but $\deg \mathring{\gamma} = 0$ and the expectation of $\tilde{\Pi}_{\epsilon} \mathring{\gamma}$ diverges logarithmically and the expression $\mathring{\gamma}$ given above fails to converge. We then have no choice but to set

$$ g_{\epsilon}(\mathring{\gamma}) = -\mathring{\gamma}_{\epsilon}, \quad \mathring{\gamma}_{\epsilon} \stackrel{\text{def}}{=} \int K_{\epsilon}(z) K_{\epsilon}(z') \kappa_2(z, z') \mathrm{d}z \mathrm{d}z', $$

where $K_{\epsilon}(t, x) = \epsilon^2 K(\epsilon^2 t, \epsilon x)$. We then have the following preliminary result:

**Proposition 4.6.** The model $\hat{\Pi}_{\epsilon}$ restricted to the translation invariant sector converges as $\epsilon \to 0$ to the BPHZ model $\hat{\Pi}$ such that $\hat{\Pi} \circ = \xi$ and $\hat{\Pi} \Xi_i^j = 0$ for $i + j > 0$.

*Proof.* Consider the “BPHZ model” $\Pi_{\epsilon}^{\text{BPHZ}}$ on $\mathcal{T}_{\circ}$ given by

$$ \Pi_{\epsilon}^{\text{BPHZ}} \tau = (g_{\epsilon}^{\text{BPHZ}} \otimes \Pi_{\epsilon}) \Delta^{-} \tau, $$

with $g_{\epsilon}^{\text{BPHZ}}$ defined as in (4.14). Thanks to Proposition 5.5, our noises are such that the norm of [6, Def. A.18] is finite, uniformly in $\epsilon$, for the cumulant homogeneity described in (5.3).

It therefore follows from [6, Thm 2.33] that $\Pi_{\epsilon}^{\text{BPHZ}}$ converges to $\hat{\Pi}$. Since furthermore the action of the character group of $\mathcal{T}_{\circ}$ on the space of admissible models is continuous [5], it suffices to show that one can write

$$ \hat{\Pi}_{\epsilon} = (\delta g_{\epsilon} \otimes \Pi_{\epsilon}^{\text{BPHZ}}) \Delta^{-}, \quad (4.18) $$

for some $\delta g_{\epsilon} \in \mathcal{T}_{\circ}^*$ with $\lim_{\epsilon \to 0} \delta g_{\epsilon} = 0$. For this, the following result is useful:
**Lemma 4.7.** Let $g, \bar{g} \in \mathcal{T}_*^*$ such that furthermore $g(\tau) = \bar{g}(\tau) = 0$ for every $\tau$ of the type (4.16). Then, one has the identity

$$ (g \otimes (\bar{g} \otimes \mathrm{id})\Delta^-)\Delta^- \tau = ((g + \bar{g}) \otimes \mathrm{id})\Delta^- \tau, $$

for all $\tau \in \mathcal{T}$.

*Proof.* It was shown in [5, 11] that

$$ (g \otimes (\bar{g} \otimes \mathrm{id})\Delta^-)\Delta^- \tau = ((g \otimes \bar{g})\Delta^- \otimes \mathrm{id})\Delta^- \tau, $$

where $\Delta^- : \mathrm{Alg}\mathcal{T}_- \to \mathrm{Alg}\mathcal{T}_- \otimes \mathrm{Alg}\mathcal{T}_-$ is an extraction/contraction operator defined just like above, but extended multiplicatively to $\mathrm{Alg}\mathcal{T}_-$ and such that only those terms are kept that actually belong to $\mathrm{Alg}\mathcal{T}_- \otimes \mathrm{Alg}\mathcal{T}_-$ (that is every factor needs to be of negative degree on both sides of the tensor product). Inspection of the list (4.15) reveals that in our case, the only situation in which we have a “subsymbol” of negative degree appearing in any of our symbol in such a way that the contracted symbol is still of negative degree is when the subsymbol contains only one noise. We conclude that

$$ (g \otimes \bar{g})\Delta^- \tau = (g \otimes \bar{g})(\tau \otimes \mathbf{1} + \mathbf{1} \otimes \tau) = g(\tau) + \bar{g}(\tau), $$

and the claim follows. $\square$

We conclude from Lemma 4.7 that (4.18) holds with $\delta g_\epsilon = g_\epsilon - g_\epsilon^{\mathrm{BPHZ}}$, so that it remains to show that $\lim_{\epsilon \to 0} \delta g_\epsilon(\tau) = 0$ for every $\tau \in \mathcal{T}_-$. For elements $\tau$ of the form $\tau = X^k \Xi_i^j$ we have $g_\epsilon(\tau) = g_\epsilon^{\mathrm{BPHZ}}(\tau) = 0$. For all other elements $\tau$ with $\mathrm{deg}\tau \le 0$, a simple scaling argument shows that $g_\epsilon^{\mathrm{BPHZ}}(\tau)$ is given by the exact same formula as $g_\epsilon(\tau)$, except that all instances of the heat kernel $P$ are replaced by $K_\epsilon$, where

$$ K_\epsilon(t, x) \stackrel{\mathrm{def}}{=} \epsilon^d K(\epsilon^2 t, \epsilon x). $$

Note that $K_\epsilon$ coincides with $P$ in a parabolic ball of radius $\mathcal{O}(1/\epsilon)$ around the origin and vanishes outside of another ball of radius $\mathcal{O}(1/\epsilon)$.

This in particular shows that

$$ \delta g_{\epsilon}(\mathring{\omega}) = \epsilon^{-d/2} \int (P(z) - K_{\epsilon}(z)) \kappa_{2}(0, z) \mathrm{d}z. \quad (4.19) $$

Since $\kappa_2(0, z)$ decreases like $\|z\|^{-2\bar{c}} = \|z\|^{-(d+2+2\delta)}$ for large $z$ and $P$ decreases like $\|z\|^{-d}$, it follows that

$$ |\delta g_{\epsilon}(\mathring{\omega})| \lesssim \epsilon^{d/2+2\delta}, $$

which of courses converges to 0. The symbols $\tau$ differing from $\mathring{\omega}$ only by the placement of their accents then also converge since $\delta g_\epsilon(\tau)$ is given by the same expression as (4.19), except for being multiplied by a higher power of $\epsilon$.
Turning now to $\mathring{\mathcal{O}}^\circ$ (which only appears when $d \in \{2, 3\}$), it follows from [15, Lem. 6.8] that $\delta g_\epsilon(\mathring{\mathcal{O}}^\circ)$ is a sum of terms of the form

$$ \epsilon^{2-d} \int (P(z) - K_\epsilon(z)) G(z) \, dz, $$

where $|G(z)| \lesssim (1 + \|z\|)^{2-2\bar{c}} = (1 + \|z\|)^{-(d+2\delta)}$. It follows that $|\delta g_\epsilon(\mathring{\mathcal{O}}^\circ)| \lesssim \epsilon^{2\delta}$ as desired, and $\delta g_\epsilon(\mathring{\mathcal{O}}^\circ)$ is controlled in the same way by a higher power. In dimension 3, $\deg \mathring{\mathcal{O}}^\circ = 0$ and it was shown in (4.17) that $\mathring{\mathcal{O}}^\circ$ converges absolutely, which immediately implies that $\delta g_\epsilon(\mathring{\mathcal{O}}^\circ) \to 0$.

To deal with the symbol $\mathring{\mathcal{O}}^\circ$ (again with $d \in \{2, 3\}$), we first note that Assumption 2.1 implies the bound

$$ |\kappa_3(z_1, z_2, z_3)| \lesssim \rho(\|z_1 - z_2\|) \rho(\|z_2 - z_3\|) \rho(\|z_1 - z_3\|). $$

We also note that for any $\kappa \in [0, d]$ one has the bound

$$ |P(z) - K_{\epsilon}(z)| \lesssim \epsilon^{\kappa} (1 \wedge \|z\|^{\kappa-d}). \quad (4.20) $$

This allows us to make use of [13, Thm 4.3]. Since $\deg \mathring{\mathcal{O}}^\circ = \frac{2-d}{2}$ we apply the bound (4.20) with $\kappa = \delta + \frac{d-2}{2}$ which, in the notation of [13], yields a bound of the type

$$ |\delta g_{\epsilon}(\mathring{\mathcal{O}}^\circ)| \le \epsilon^{\delta} |\Pi^{\bar{K}, \bar{R}} \Gamma|, $$

for some $\epsilon$-dependent kernel assignment $(\bar{K}, \bar{R}) \in \mathcal{K}_0^- \times \mathcal{K}_0^+$ with bounds that are independent of $\epsilon$ and the Feynman diagram

<figure><img src="image_3.png" /></figure>

Here, the first coordinate of the label for each edge denotes its small-scale degree assignment while the second coordinate denotes its large-scale degree assignment. It is straightforward to verify that the small-scale degree assignment for this diagram satisfies the assumption of [13, Prop. 2.3], so that it does not require renormalisation. Furthermore, the large-scale degree assignment is seen to satisfy the assumption of [13, Thm 4.3], which guarantees that the integral converges absolutely and is bounded independently of $\epsilon$, so that $|\delta g_\epsilon(\mathring{\mathcal{O}}^\circ)| \le \epsilon^\kappa$.

The symbols $\mathring{\mathcal{O}}^\circ$ and $\mathring{\mathcal{O}}_\circ$ (in dimension 3) which have vanishing degree $\deg$, can be dealt with using the same technique, leading to the bound $|\delta g_\epsilon(\mathring{\mathcal{O}}^\circ)| + |\delta g_\epsilon(\mathring{\mathcal{O}}_\circ)| \lesssim \epsilon^{1/2}$ by using the Feynman diagrams

<figure><img src="image_4.png" /></figure>
for bounding ⌢ and

<figure><img src="image_5.png" /></figure>

for bounding ⌢. All three are easily seen to satisfy both the small-scale and large-scale integrability conditions.

The remaining three symbols $\tau \in \{\stackrel{\circ}{\tau}_n, \stackrel{\circ}{\tau}_u, \stackrel{\circ}{\tau}_o\}$ in the list (4.15) are all such that $\deg \tau > 0$, so we need to show that $g_{\varepsilon}^{\text{BPHZ}}(\tau) \to 0$. This can in principle be shown again by using the bounds from [13]. A “cheaper” way of showing that $g_{\varepsilon}^{\text{BPHZ}}(\tau) \to 0$ is to note that in all three cases we can make use of a combination of Proposition 5.5 (used in the same way as in the proof of Proposition 5.9) and (4.14) to conclude that one can build a regularity structure $\hat{T}$ extending $T_0$ (by adding additional “noises” representing $\eta_{\varepsilon}^{(\alpha)}$ for suitable choices of $\alpha$) such that, for every $\tau$, one can find $\kappa > 0$ and a symbol $\hat{\tau} \in \hat{T}$ such that $\deg \hat{\tau} > 0$ and such that

$$g_{\varepsilon}^{\text{BPHZ}}(\tau) = \varepsilon^{\kappa} \mathbf{E}(\Pi_{\varepsilon}^{\text{BPHZ}} \hat{\tau})(\phi),$$

for some suitable fixed test function $\phi$. Since we know from [6, Thm 2.33] that the BPHZ renormalised model $\Pi_{\varepsilon}^{\text{BPHZ}}$ converges (so in particular remains uniformly bounded), we conclude that $g_{\varepsilon}^{\text{BPHZ}}(\tau) \to 0$ as required. $\square$

**Theorem 4.8.** The random models $\hat{\Pi}_{\varepsilon}$ converge weakly to a limiting admissible model $\hat{\Pi}$ which, on the translation invariant sector, is given by the BPHZ lift of

$$\hat{\Pi} \circ = \xi, \quad \hat{\Pi} \Xi_i^j = 0, \quad i + j > 0.$$

For the remaining symbols, it is given by the unique admissible model such that

$$\hat{\Pi} \diamond = \xi \mathbf{1}_{\mathbb{R} \times D}, \quad \hat{\Pi} \diamond, = 0,$$

as well as $\hat{\Pi} \tau = 0$ for any symbol $\tau$ containing the noise $\diamond$.

**Remark 4.9.** Note that only symbols $\Xi_i^j$ with $i \ge j$ appear in our regularity structure.

*Proof.* Convergence on the translation invariant sector was already shown in Proposition 4.6, so it only remains to consider the non-translation invariant symbols of negative degree. In dimension 3, these are $\diamond$, $\diamond$, $\diamond$, $\uparrow$, $\downarrow$, $\circled$, $\circled$, $\circled$, and $\circled$. (There is also the symbol $\diamond$, $\mathcal{I}_i(\circ)$, but applying the model to it yields the exact same distribution as when applying it to $\circled, \downarrow$.)

The convergence on the remainder of the regularity structure is shown in the next section, but we collect the various parts of the proof here. Convergence of
$\hat{\Pi}_{\epsilon} \diamond$, $\hat{\Pi}_{\epsilon} \diamond$, and $\hat{\Pi}_{\epsilon} \diamond$ to $\xi \mathbf{1}^D$ and $0$ in $C^{-\frac{5}{2}-\kappa}$, $C^{-1-\kappa}$ and $C^{-\frac{3}{2}-\kappa}$ respectively follows from Corollary 5.13 and Corollary 5.6.

Convergence of $\hat{\Pi}_{\epsilon} \diamond$, essentially follows from [6, Thm 2.31], noting that the bound for $\tau = \diamond$ does not require any derivative of the test function in this case, so that we immediately obtain the required bound by noting that $(\hat{\Pi}_{\epsilon} \diamond)(\phi) = (\hat{\Pi}_{\epsilon} \diamond)(\mathbf{1}_{\mathbb{R} \times D}\phi)$. The reason why this is so is that the only point in the proof where derivatives of the test function could potentially appear is in the bound [6, Eq. A.29] in the proof of [6, Thm A.32]. By the definition of $R(S)$, these derivatives can only hit a test function in a situation where $\tau$ contains a connected subtree containing its root and of degree less than $-\frac{d+2}{2}$. This is not the case for $\diamond$.

Convergence of $\hat{\Pi}_{\epsilon} \diamond$, $\hat{\Pi}_{\epsilon} \diamond$, and $\hat{\Pi}_{\epsilon} \diamond$, also follows in the same way. Regarding $\diamond$, we note that $\diamond = \bigcirc + \diamond$, by (4.23) and we already obtained convergence of $\hat{\Pi}_{\epsilon} \diamond$. Convergence of $\hat{\Pi}_{\epsilon} \bigcirc$, is the content of Theorem 5.20. The convergence of $\hat{\Pi}_{\epsilon} \uparrow$ follows from Corollary A.3, combined with the usual Schauder estimates for integration against $K$. $\square$

Most of Section 5 is devoted to filling in the missing parts in the proof of Theorem 4.8, namely the proofs of Theorem 5.11 and Theorem 5.20.

### 4.4. Description of $v_{\epsilon}^{(0)}$

We now have the notation in place to formalise the discussion given above regarding the local behaviour of $v_{\epsilon}^{(0)}$ and $v_{\epsilon}^{(1)}$. Recall the definition (4.5) of $v_{\epsilon}^{(0)}$ which we rewrite in integral form as

$$v_{\epsilon}^{(0)} = K_{\text{Neu}}(G(u^{(0)}) \mathbf{1}_{+}^{D}\xi_{\epsilon}). \quad (4.21)$$

Depending on context, we will model $v_{\epsilon}^{(0)}$ by three different modelled distributions $V_{\epsilon}^{(0)}$, $\tilde{V}_{\epsilon}^{(0)}$ and $\hat{V}_{\epsilon}^{(0)}$.

Regarding $V_{\epsilon}^{(0)}$, we use Proposition 6.1, which guarantees that one can find $\Phi_{\epsilon} \in C^{\frac{3}{2}-\kappa, -\frac{1}{2}-\kappa}$ such that, setting$^4$

$$V_{\epsilon}^{(0)} = \Phi_{\epsilon} + \mathcal{L}G(u^{(0)}) K_{\text{Neu}}(\mathbf{1}_{+} \diamond) + G'(u^{(0)}) \partial_{i} u^{(0)} K_{\text{Neu},i}(\mathbf{1}_{+} \diamond), \quad (4.22)$$

where the integration operators $K_{\text{Neu}}$ and $K_{\text{Neu},i}$ should be interpreted as the natural extensions of the corresponding integration operators described in Proposition 6.1, one has

$$V_{\epsilon}^{(0)} \in D^{\frac{3}{2}-\kappa, -\frac{1}{2}-\kappa}, \quad \mathcal{R}_{\epsilon} V_{\epsilon}^{(0)} \upharpoonright \mathbf{R}_{+} \times D = v_{\epsilon}^{(0)},$$

provided that we consider admissible models $\hat{\Pi}_{\epsilon}$ with $\hat{\Pi}_{\epsilon} \diamond = \xi_{\epsilon} \mathbf{1}_{\mathbb{R} \times D}$. Note that $V_{\epsilon}^{(0)}$ depends on $\epsilon$ via the choice of model $\hat{\Pi}_{\epsilon}$. Furthermore, Theorem 4.8 and Lemma 5.14 guarantee that the bounds on $V_{\epsilon}^{(0)}$ are uniform over $\epsilon$ and that one has

<footnote>$^4$ with the convention $\mathcal{L} = \mathcal{L}_1$ for $\mathcal{L}_1$ the Taylor lift (3.3)</footnote>
$V_{\varepsilon}^{(0)} \to V^{(0)}$ in $D^{\frac{3}{2}-\kappa, -\frac{1}{2}-\kappa}$ with respect to the renormalised models $\hat{\Pi}_{\varepsilon}$ and the limiting model $\hat{\Pi}$.

It is natural at this stage, as already mentioned earlier, to define the symbol ♭ as the element of $\mathcal{T}$ of degree $-\frac{1}{2} - \kappa$ given by

$$ \mathfrak{A} = \mathfrak{A} - \mathfrak{A}, \qquad (4.23) $$

which is indeed consistent with (4.7). With this notation, Propositions A.2 and A.4 combined with (4.22) guarantee the existence of a function $\tilde{\Phi}_{\varepsilon} \in C^{\frac{3}{2}-\kappa, w}$ with $w = (-\frac{1}{2}-\kappa, \frac{1}{2}-\kappa, -\frac{1}{2}-\kappa)$ in the sense of [9, Def. 3.2] such that if we set $\tilde{V}_{\varepsilon}^{(0)} = \tilde{V}_{\varepsilon}^{(0,1)} + \tilde{V}_{\varepsilon}^{(0,2)}$ with

$$ \begin{aligned} \tilde{V}_{\varepsilon}^{(0,1)} &= \mathcal{L}G(u^{(0)})\mathfrak{A}, \\ \tilde{V}_{\varepsilon}^{(0,2)} &= \mathcal{L}G(u^{(0)})\mathfrak{A} + G'(u^{(0)})\partial_i u^{(0)} K_i(\mathfrak{A}) + \mathcal{L}\tilde{\Phi}_{\varepsilon}, \end{aligned} \qquad (4.24) $$

then we have $\tilde{V}_{\varepsilon}^{(0,1)} \in D^{\frac{3}{2}-\kappa}$, $\tilde{V}_{\varepsilon}^{(0,2)} \in D^{\frac{3}{2}-\kappa, w}$ and $(\mathcal{R}_{\varepsilon}\tilde{V}_{\varepsilon}^{(0)})(\phi) = (\mathcal{R}_{\varepsilon}V_{\varepsilon}^{(0)})(\phi)$ for all test functions $\phi$ supported in $\mathbf{R}_{+} \times D$. As before, all these objects converge in the limit $\varepsilon \to 0$ provided that the underlying models converge. To see this, we note that we can choose

$$ \begin{aligned} \tilde{\Phi}_{\varepsilon} ={}& \Phi_{\varepsilon} + \mathbf{1}_{+}G(u^{(0)})K_{\text{Neu}}((\mathbf{1}_{+} - 1)\mathbf{1}_{D}\xi_{\varepsilon}) \\ & + \mathbf{1}_{+}^{D}G'(u^{(0)})\partial_{i}u^{(0)} K_{i} * ((\mathbf{1}_{+}^{D} - 1)\xi_{\varepsilon}) \\ & + \mathbf{1}_{+}^{D}G'(u^{(0)})\partial_{i}u^{(0)} K_{\partial,i}(\mathbf{1}_{+}^{D}\xi_{\varepsilon}), \end{aligned} $$

and then apply Proposition A.4 to bound the first term and Proposition A.2 to bound the remaining two terms.

On the other hand, we define $\hat{V}_{\varepsilon}^{(0)}$ by setting

$$ \hat{V}_{\varepsilon}^{(0)} = \mathcal{K}(\mathcal{L}_2 G(u^{(0)})\mathbf{1}_{+\circ}) + \mathcal{L}\mathbf{1}_{+}^{D}(K_{\partial}(G(u^{(0)})\xi_{\varepsilon}\mathbf{1}_{+}^{D}) - K(G(u^{(0)})\xi_{\varepsilon}\mathbf{1}_{+}^{D^{\varepsilon}})), $$

where we use the convention that $u^{(0)}$ is extended outside of $\mathbf{R}_{+} \times D$ in any way that makes it globally $C^3$. (For positive times, this is possible since the extension of $u_0$ to the whole space by suitable reflections is of class $C^3$ by our assumptions. For negative times, this is possible by Whitney's extension theorem [21].) Here we made a slight abuse of notation: the operator $\mathcal{K}$ should be interpreted in the sense of [9, Sec. 4.5] with $\hat{\mathcal{R}}(\mathcal{L}_2 G(u^{(0)})\mathbf{1}_{+\circ}) = G(u^{(0)})\mathbf{1}_{+}\xi_{\varepsilon}$, which converges as $\varepsilon \to 0$ by an argument very similar to that of Proposition 6.1, combined with the fact that $G(u^{(0)})$ is $C^3$.

Note that we have

$$ \mathcal{R}_{\varepsilon}\hat{V}_{\varepsilon}^{(0)} = \mathcal{R}_{\varepsilon}\tilde{V}_{\varepsilon}^{(0)} = \mathcal{R}_{\varepsilon}V_{\varepsilon}^{(0)} = v_{\varepsilon}^{(0)} \quad \text{on } \mathbf{R}_{+} \times D. $$

The modelled distribution $\hat{V}_{\varepsilon}^{(0)}$ exhibits rather singular behaviour near the boundary of the domain, but by Proposition A.2 it converges as $\varepsilon \to 0$ in $D^{1,w}$ with $w = (-\frac{1}{2}-\kappa)_3$ [we use again the notation $(\eta)_3 = (\eta, \eta, \eta)]$. It has the advantage however of not involving the integration map $\mathcal{I}$ and the restricted noise $\diamond$, so it is purely described in terms of the “translation invariant” part of the regularity structure, that is the sector generated by $\mathcal{I}$ and $\circ$. We summarise the above discussion with the following statement.
**Lemma 4.10.** We have $\hat{V}_{\varepsilon}^{(0)} \in D_{-\frac{1}{2}-\kappa}^{1,(-\frac{1}{2}-\kappa)3}$, $\mathcal{R}_{\varepsilon}\hat{V}_{\varepsilon}^{(0)} = v_{\varepsilon}^{(0)}$ on $\mathbf{R}_{+} \times D$, and $\hat{V}_{\varepsilon}^{(0)}$ is of the form

$$ \hat{V}_{\varepsilon}^{(0)} = \mathbf{1}_{+}(G(u^{(0)}) \diamond + G'(u^{(0)}) \nabla u^{(0)} \diamond) + \Phi_{\varepsilon}^{(0)} \mathbf{1}, $$

for some continuous function $\Phi_{\varepsilon}^{(0)}$. Furthermore, $\lim_{\varepsilon \to 0} \hat{V}_{\varepsilon}^{(0)} = \hat{V}^{(0)} \in D^{1,(-\frac{1}{2}-\kappa)3}$ with $\mathcal{R}_0 \hat{V}^{(0)} = P_{\text{Neu}}(\mathbf{1}_{+}^D G(u^{(0)}))\xi$.

### 4.5. Description of $v_{\varepsilon}^{(1)}$

Recall that $v_{\varepsilon}^{(1)}$ was defined in (4.6) as the solution to an inhomogeneous linear equation with inhomogeneous boundary conditions (et least in the Neumann case). Regarding $v_{\varepsilon}^{(1)}$, we would like to describe it by a modelled distribution $V_{\varepsilon}^{(1)}$ given by

$$ V_{\varepsilon}^{(1)} = \mathbf{1}_{+}^{D} \mathcal{P}_{\text{Neu}} \mathbf{1}_{+} (\mathcal{L}(G'(u^{(0)})) \tilde{V}_{\varepsilon}^{(0)} \diamond_{\varepsilon}) + v_{\varepsilon}^{(1,*)} \mathbf{1}, \quad (4.25) $$

with $\mathcal{P}_{\text{Neu}}$ as in (3.2) and where $v_{\varepsilon}^{(1,*)}$ would be given by the solution to

$$ \begin{aligned} \partial_t v_{\varepsilon}^{(1,*)} &= \Delta v_{\varepsilon}^{(1,*)} + \frac{1}{2} G''(u^{(0)}) W_{1,\varepsilon}, \\ W_{1,\varepsilon} &= (v_{\varepsilon}^{(0)})^2 \sigma_{\varepsilon} - 2 \mathring{\circ} G(u^{(0)}) v_{\varepsilon}^{(0)} - \varepsilon^{1-\frac{d}{2}} \mathring{\circ} \mathring{\circ} G^2(u^{(0)}), \end{aligned} \quad (4.26) $$

endowed with homogeneous Neumann boundary conditions. The reason why the identity $\mathcal{R}_{\varepsilon}V_{\varepsilon}^{(1)} = v_{\varepsilon}^{(1)}$ holds (on $\mathbf{R}_{+} \times D$ as usual) is as follows. By (4.24),

$$ \tilde{V}_{\varepsilon}^{(0)} \diamond_{\varepsilon} = \mathcal{L}(G(u^{(0)}))(\mathring{\circ}_{\varepsilon} + \mathring{\circ}_{\varepsilon}) + G'(u^{(0)})\partial_{\varepsilon}u^{(0)} \diamond_{\varepsilon} \mathcal{K}_{\varepsilon}(\mathring{\circ}) + \mathcal{L}\tilde{\Phi}_{\varepsilon} \diamond_{\varepsilon}. \quad (4.27) $$

It then follows from the definition (4.13) of the renormalised model and the fact that $\mathcal{R}_{\varepsilon}(\tilde{V}_{\varepsilon}^{(0)}) = v_{\varepsilon}^{(0)}$ that, applying the reconstruction operator to this expression yields

$$ \begin{aligned} \mathcal{R}_{\varepsilon}(\tilde{V}_{\varepsilon}^{(0)} \diamond_{\varepsilon}) ={}& \mathbf{1}^{D} v_{\varepsilon}^{(0)} \eta_{\varepsilon} - \sum_{i=1}^{d} G(u^{(0)}) \varepsilon^{1-\frac{d}{2}} (c_{i,0} \delta_{\partial_{i,0}D} + c_{i,1} \delta_{\partial_{i,1}D}) \\ & - \varepsilon^{-\frac{d}{2}} \mathring{\circ} \mathbf{1}^{D} G(u^{(0)}) - \varepsilon^{1-\frac{d}{2}} \mathbf{1}^{D} G'(u^{(0)}) (\mathring{\circ}, \nabla u^{(0)}). \end{aligned} $$

We then note that the two terms appearing on the second line, when multiplied by another factor of $G(u^{(0)})$, are exactly the additional two terms appearing on the first line of (4.6), while the singular term involving Dirac masses on the boundary of $D$, when hit by $K_{\text{Neu}}$, is responsible for the non-homogeneous boundary conditions.

The problem with such a definition is that it yields a description of $V_{\varepsilon}^{(1)}$ in terms of symbols involving $\diamond$, while the general convergence results of [6] require translation invariance of the noise objects, which is not the case here. If we were to try to improve the situation by replacing $\mathbf{1}_{+} \diamond$, by $\mathbf{1}_{+}^{D} \mathring{\circ}$ in (4.25), then we immediately run into the problem that the behaviour of this modelled distribution on the boundary of $D$ is too singular for the general results of [9] to apply. This is
not just a technicality: this singular behaviour is precisely what is responsible for
the additional boundary renormalisation!

Instead, we define $V_{\varepsilon}^{(1)}$ as a sum of terms that is “equivalent” to the definition (4.25) in the sense that they reconstruct the same function, but such that each of the terms can be controlled in a slightly different, situation-specific, way.

We first deal with the boundary correction by setting

$$
V_{\varepsilon}^{(1,0)} = \mathbf{1}_{+}^{D} \mathcal{L} K_{\partial} (\mathbf{1}_{+} G'(u^{(0)}) \mathcal{R}_{\varepsilon} (\tilde{V}_{\varepsilon}^{(0)} \diamond)) . \quad (4.28)
$$

Since $\tilde{V}_{\varepsilon}^{(0)} \diamond, \in D^{\frac{1}{2}-2\kappa, w}$ with $w = (-\frac{3}{2}-2\kappa, -\frac{1}{2}-2\kappa, -\frac{3}{2}-2\kappa)$ and since it belongs to a sector of regularity $-\frac{3}{2}-2\kappa$, its reconstruction belongs to $C^{-\frac{3}{2}-2\kappa}$. It then follows from Proposition A.2 that $V_{\varepsilon}^{(1,0)} \in D^{2,(0)3}$.

We now break up $\tilde{V}_{\varepsilon}^{(0)}$ in (4.25) as in (4.24) and deal with the first term. By Proposition 6.3, we can find $V_{\varepsilon}^{(1,1)} \in D^{2-2\kappa, \tilde{w}}$ with $\tilde{w} = (\frac{1}{2}-2\kappa, \frac{1}{2}-2\kappa, 0)$ of the form

$$
V_{\varepsilon}^{(1,1)} = \mathbf{1}_{+}^{D} G'(G(u^{(0)})(\hat{\Pi}_{\varepsilon} \perp)) \circ + \Phi_{\varepsilon}^{(1,1)}, \quad (4.29)
$$

with $\Phi_{\varepsilon}^{(1,1)}$ taking values in the classical Taylor polynomials, and such that

$$
\mathcal{R}_{\varepsilon} V_{\varepsilon}^{(1,1)} = K \left( (G'G)(u^{(0)}) \mathbf{1}_{+} \hat{\Pi}_{\varepsilon} \perp \right) = \mathcal{R}_{\varepsilon} K \mathbf{1}_{+} \left( \mathcal{L}(G'(u^{(0)})) \tilde{V}_{\varepsilon}^{(0,1)} \diamond \right).
$$

The second term is dealt with similarly. As a consequence of Proposition 6.4 with $g_1 = G'G(u^{(0)})$, $g_{2,i} = G'(u^{(0)})^2 \partial_i u^{(1)}$ and $g_3 = G'(u^{(0)}) \tilde{\Phi}_\varepsilon$ (with $\tilde{\Phi}_\varepsilon$ as in (4.25)), we can find $V_\varepsilon^{(1,2)} \in D^{2-2\kappa, \tilde{w}}$ such that

$$
\mathcal{R}_{\varepsilon} V_{\varepsilon}^{(1,2)} = K \mathbf{1}_{+} (G'G(u^{(0)})) \hat{\Pi}_{\varepsilon} \circ, \\
\qquad + G'(u^{(0)})^2 \partial_i u^{(1)} \mathcal{R}_{\varepsilon}(\diamond, K_i(\circ)) + G'(u^{(0)}) \tilde{\Phi}_{\varepsilon} \hat{\Pi}_{\varepsilon} \diamond,
$$

and such that furthermore $V_{\varepsilon}^{(1,2)}$ takes values in the translation invariant sector and
is of the form

$$
V_{\varepsilon}^{(1,2)} = \mathbf{1}_{+}^{D} \left( (G'G)(u^{(0)}) \circ_{\diamond} + \tilde{\Phi}_{\varepsilon} \circ_{\diamond} \right) + \Phi_{\varepsilon}^{(1,2)}, \quad (4.30)
$$

for some $\Phi_{\varepsilon}^{(1,2)}$ taking values in the Taylor polynomials. In order to define $V_{\varepsilon}^{(1,3)}$,
we make use of the following lemma:

**Lemma 4.11.** Let $\phi_{\varepsilon}$ be such that on $\mathbf{R}_{+} \times D$ one has the identity

$$
v_{\varepsilon}^{(0)} = G(u^{(0)}) \hat{\Pi}_{\varepsilon} \circ + \phi_{\varepsilon} \mathbf{1},
$$

and one has $\phi_{\varepsilon}(t, x) = 0$ for $t < 0$ or $x \notin D$. Then, for any $\alpha \in [0, 1)$, one has the bound $\mathbf{E}\|\phi_{\varepsilon}\|_{\alpha+\frac{1}{2}-\kappa, w} \lesssim \varepsilon^{-\alpha}$ with $w = (\alpha - \frac{1}{2} - \kappa)_3$.
*Proof.* We decompose $\phi_\epsilon$ as

$$ \phi_\epsilon = K_\partial (\mathbf{1}_+^D G(u^{(0)}) \xi_\epsilon) + G(u^{(0)}) K((1 - \mathbf{1}_+^D) \xi_\epsilon) + (K(\mathbf{1}_+^D G(u^{(0)}) \xi_\epsilon) - G(u^{(0)}) K(\mathbf{1}_+^D \xi_\epsilon)), $$

and we treat the three terms separately. The first two terms are estimated by combining Proposition A.2 with Corollary 5.10.

The bound on the last term follows from combining Proposition 5.9 with Corollary B.5. To apply the latter, we set $\theta = \kappa$ (small enough) and $\chi = \alpha - \frac{5}{2} - \kappa$, which yields a bound in $C^{\kappa+2}$ on

$$ K(\mathbf{1}_+^D G(u^{(0)}) \xi_\epsilon) - G(u^{(0)}) K(\mathbf{1}_+^D \xi_\epsilon) - \sum_i G'(u^{(0)}) \partial_i u^{(0)} K_i(\mathbf{1}_+^D \xi_\epsilon). $$

Since $K_i$ gains three derivatives, the term $K_i(\mathbf{1}_+^D \xi_\epsilon)$ itself satisfies the required bound and we are done. $\square$

Recalling $W_{1,\epsilon}$ as defined in (4.26), it follows from Lemma 4.11 combined with the definition of the renormalised model that we can rewrite it as

$$ W_{1,\epsilon} = G^2(u^{(0)}) \hat{\Pi}_{\epsilon} \circ_{\omega} + 2G(u^{(0)}) \phi_{\epsilon} \hat{\Pi}_{\epsilon} \circ_{\omega} + \phi_{\epsilon}^2 \hat{\Pi}_{\epsilon} \circ_{\omega}. \quad (4.31) $$

As a consequence of Proposition 5.9 (with $\alpha = \kappa$), combined with Lemma 4.11 (with $\alpha = \frac{1}{2}$), we conclude that one has

$$ \lim_{\epsilon \to 0} \| W_{1,\epsilon} \|_{-\frac{1}{2}-3\kappa} = 0. \quad (4.32) $$

Indeed, the reconstruction theorem [9, Thm 4.9] and the multiplication rules [9, Lem. 4.3] imply that if $g \in C^{\gamma,(\eta)_3}$ for $\gamma \le 0$ and $\gamma > 0$, and $\zeta \in C^\beta$ with $\beta \le 0$ and $\gamma + \beta > 0$ then, provided that $\eta + \beta > -1$, one has $g\zeta \in C^{\eta+\beta}$. (View $\zeta$ as the constant function in a regularity structure containing only one symbol of degree $\beta$ and apply the reconstruction theorem to $g\zeta$.)

Since $-\frac{1}{2} - 3\kappa > -1$, we can multiply such a distribution by the indicator function of $\mathbf{R}_+ \times D$. It follows that, setting

$$ V_{\epsilon}^{(1,3)} \stackrel{\text{def}}{=} \frac{1}{2} \mathcal{L} K_{\text{Neu}} (\mathbf{1}_{D}^{+} G''(u^{(0)}) W_{1,\epsilon}), \quad (4.33) $$

we have $\lim_{\epsilon \to 0} V_{\epsilon}^{(1,3)} = 0$ in $D^{\frac{3}{2}-3\kappa}$ and furthermore $R_{\epsilon} V_{\epsilon}^{(1,3)} = v_{\epsilon}^{(1,\star)}$. Combining these definitions, we set

$$ \hat{V}_{\epsilon}^{(1)} = V_{\epsilon}^{(1,0)} + V_{\epsilon}^{(1,1)} + V_{\epsilon}^{(1,2)} + V_{\epsilon}^{(1,3)}. $$

Summarising this discussion, one has the following result:

**Lemma 4.12.** We have $\hat{V}_{\epsilon}^{(1)} \in D^{\frac{3}{2}-3\kappa,(0)_3}$, $R_{\epsilon} \hat{V}_{\epsilon}^{(1)} = v_{\epsilon}^{(1)}$ and $\hat{V}_{\epsilon}^{(1)}$ is of the form

$$ \hat{V}_{\epsilon}^{(1)} = \mathbf{1}_{+}^{D} G' G(u^{(0)}) \circ_{\omega} + \mathbf{1}_{+}^{D} G'(u^{(0)}) \Phi_{\epsilon}^{(0)} \circ_{\omega} + \Phi_{\epsilon}^{(1)}, \quad (4.34) $$

for some $\Phi_{\epsilon}^{(0)}$ taking values in the Taylor polynomials, and where $\Phi_{\epsilon}^{(0)}$ is as in Lemma 4.10. Furthermore, $\lim_{\epsilon \to 0} \hat{V}_{\epsilon}^{(1)} = \hat{V}^{(1)} \in D^{\frac{3}{2}-3\kappa,(0)_3}$ with $R_0 \hat{V}^{(1)} = 0$.
*Proof.* Collecting (4.28), (4.29), (4.30) and (4.33), we see that (4.34) holds, but with $\Phi_{\varepsilon}^{(0)}$ replaced by $G(u^{(0)})(\hat{\mathbf{\Pi}}_{\varepsilon} \mathbf{1}) + \tilde{\Phi}_{\varepsilon}$. These two expression are seen to coincide on $\mathbf{R}_{+} \times D$ by comparing Lemma 4.10 with (4.24).

The only statement we haven't shown yet is that $\mathcal{R}_0 \hat{V}^{(1)} = 0$. Since we already know by (4.32) that $W_{1,\varepsilon}$ converges to 0 and since $\mathcal{R}_{\varepsilon} \hat{V}_{\varepsilon}^{(1)} = \mathcal{R}_{\varepsilon} V_{\varepsilon}^{(1)}$, it remains by (4.25) to show that $\lim_{\varepsilon \to 0} \mathcal{R}_{\varepsilon} (\tilde{V}_{\varepsilon}^{(0)} \diamond_{\varepsilon}) = 0$. This in turn is immediate from (4.27) when combined with Theorem 4.8 which guarantees that the limiting model vanishes on all accented symbols. $\square$

## 4.6. Formulation of the Fixed Point Problem

Introduce now a modelled distribution $\overline{V}_{\varepsilon}$ and, using the shorthand $\hat{V}_{\varepsilon} = \hat{V}_{\varepsilon}^{(0)} + \hat{V}_{\varepsilon}^{(1)} + \overline{V}_{\varepsilon}$, consider the fixed point problem

$$ \overline{V}_{\varepsilon} = \mathcal{P}_{\text{Neu}} \mathbf{1}_{+}^{D} (H'_{\eta}(u^{(0)}) \hat{V}_{\varepsilon} + \mathcal{L}(G'(u^{(0)}))(\hat{V}_{\varepsilon} - \hat{V}_{\varepsilon}^{(0)}) \circ_{\varepsilon} + \mathcal{L}(G'(u^{(0)})u^{(1)}) \circ_{\varepsilon} \\ + \frac{1}{2} \mathcal{L}(G''(u^{(0)}))(\hat{V}_{\varepsilon}^2 - (\hat{V}_{\varepsilon}^{(0)})^2) \circ_{u} + \mathcal{L}(G''(u^{(0)})u^{(1)}) \hat{V}_{\varepsilon} \circ_{u} \\ + \frac{1}{2} \mathcal{L}(G''(u^{(0)})(u^{(1)})^2) \circ_{\varepsilon}^{\mathfrak{a}}) + \mathcal{L} \mathcal{P}_{\text{Neu}} \mathbf{1}_{+}^{D} (\hat{R}_{\varepsilon}^{(d)}(\mathcal{R}_{\varepsilon} \hat{V}_{\varepsilon}, \zeta) + \tilde{R}_{\varepsilon}^{(d)}(\tilde{\zeta})), \quad (4.35) $$

where $\tilde{R}_{\varepsilon}^{(d)} = 0$ for $d = 1$ and

$$ \begin{aligned} \tilde{R}_{\varepsilon}^{(2)}(\tilde{\zeta}) &= \varepsilon^{\kappa} \frac{1}{2} G^2 G''(\zeta) \circ_{\varepsilon}^{\mathfrak{a}} \tilde{\zeta}, \\ \tilde{R}_{\varepsilon}^{(3)}(\tilde{\zeta}) &= \varepsilon^{\kappa} (G''G) \left( \frac{1}{2} \zeta \circ_{\varepsilon}^{\mathfrak{a}} G + \varepsilon G' (\zeta \circ_{\varepsilon}^{\mathfrak{a}} G + \langle \zeta \circ_{\varepsilon}^{\mathfrak{a}}, \nabla u^{(0)} \rangle + \zeta \circ_{\varepsilon}^{\mathfrak{a}} u^{(1)}) \right) \tilde{\zeta}. \end{aligned} \quad (4.36) $$

**Remark 4.13.** We will set this up as a fixed point problem in the space $D^{\frac{3}{2}-3\kappa,(0)3}$. Since $\deg \circ_{\varepsilon} < -1$ and $\deg \circ_{\varepsilon} < -1$ (in $d \in \{2, 3\}$ for the latter), this forces us to rely on Theorem C.5 for the reconstruction of the right hand side of (4.35) and to combine this with [9, Lem. 4.12] to provide an interpretation for the integration operator $\mathcal{K}$ appearing in the definition (3.2) of $\mathcal{P}_{\text{Neu}}$.

**Remark 4.14.** Recall that the definition (4.3) of the remainder $\hat{R}_{\varepsilon}^{(d)}$ involves an arbitrary exponent $\alpha$. We henceforth fix a choice $\alpha = \alpha(d)$ depending on the dimension, namely

$$ \alpha(1) = \frac{5}{4}, \quad \alpha(2) = \frac{9}{4}, \quad \alpha(3) = \frac{11}{4}. \quad (4.37) $$

All further statements about $\hat{R}_{\varepsilon}^{(d)}$ hold for this particular choice.
We claim that with this definition and provided that we consider the renormalised model constructed in Section 4.3, (4.35) admits a unique solution in $D^{\frac{3}{2}-3\kappa,(0)3}$ and, provided that we set $\zeta = \epsilon^\alpha \xi_\epsilon$ and $\bar{\zeta} = \epsilon^{-1-\kappa}\sigma_\epsilon$, one has $v_\epsilon = \mathcal{R}_\epsilon \hat{V}_\epsilon$. The reason for the appearance of $\tilde{R}_\epsilon^{(d)}$ is to cancel out some additional unwanted terms arising from the renormalisation procedure. Before this, we formulate a technical lemma, where we write $\|\theta\|_\alpha$ for the $C^\alpha$ norm of the function / distribution $\theta$ on $D_T = [0, T] \times D$ with $T$ as in Theorem 1.2.

**Lemma 4.15.** Let $w, \bar{w}$ with $\|w\|_{L^\infty} + \|\bar{w}\|_{L^\infty} \le \epsilon^{-d/2}$ on the domain $D_T$ and let $\kappa \in (0, \frac{1}{4})$. Writing $X = \|\zeta\|_{-\frac{1}{2}+2\kappa} + 1$, one has the bounds

$$ \|R_{\epsilon}^{(1)}(w, \zeta)\|_{-\frac{1}{2}+2\kappa} \lesssim \epsilon^{1/4}(1 + \|w\|_{\frac{1}{2}-\kappa})^3 X, \\ \|R_{\epsilon}^{(1)}(w, \zeta) - R_{\epsilon}^{(1)}(\bar{w}, \zeta)\|_{-\frac{1}{2}+2\kappa} \lesssim \epsilon^{1/4}\|w - \bar{w}\|_{\frac{1}{2}-\kappa}(1 + \|w\|_{\frac{1}{2}-\kappa} + \|\bar{w}\|_{\frac{1}{2}-\kappa})^3 X, $$

for some proportionality constants depending only on $u^{(0)}$, $G$ and $H$. In dimensions 2 and 3, we set $X = \|\zeta\|_{L^p} + 1$ (for any fixed $p \in [1, \infty)$) and we have the bounds

$$ \|R_{\epsilon}^{(d)}(w, \zeta)\|_{L^p} \lesssim \epsilon^{\kappa}(1 + \epsilon^{\beta}\|w\|_{L^{\infty}})^3 X, \\ \|R_{\epsilon}^{(d)}(w, \zeta) - R_{\epsilon}^{(d)}(\bar{w}, \zeta)\|_{L^p} \lesssim \epsilon^{\kappa+\beta}\|w - \bar{w}\|_{L^{\infty}}(1 + \epsilon^{\beta}\|w\|_{L^{\infty}} + \epsilon^{\beta}\|\bar{w}\|_{L^{\infty}})^3 X, \\ \text{with } \beta(2) = \frac{1}{4} - \frac{\kappa}{3} \text{ and } \beta(3) = \frac{7}{12} - \frac{\kappa}{3} \text{ for } \kappa \text{ sufficiently small.} $$

*Proof.* The case of dimensions 2 and 3 is straightforward to verify since all bounds are uniform. In dimension 1, the first term of (4.3) is easy to bound. To bound the second term, we use the fact that composition with a smooth function is a (locally) Lipschitz continuous operation in $C^{\frac{1}{2}-\kappa}$, combined with the fact that the product is continuous as a bilinear map from $C^{\frac{1}{2}-\kappa} \times C^{-\frac{1}{2}+2\kappa}$ into $C^{-\frac{1}{2}+2\kappa}$, see [3] or [8, Thm 13.16]. $\square$

**Proposition 4.16.** Fix an initial condition $u_0$, a final time $T < 1$ and nonlinearities $G$ and $H$, all as in Theorem 1.2, as well as the random model $\hat{\Pi}_\epsilon$ as defined in Section 4.3. Choose $\zeta \in L^p$ with $p = (d+2)/c$ (for $d \in \{2, 3\}$) or $\zeta \in C^{2\kappa-\frac{1}{2}}$ (for $d=1$), as well as $\bar{\zeta} \in C^{-\frac{1}{2}-2\kappa}$ (for $d \in \{2, 3\}$). Then, the fixed point problem (4.35) admits a unique local solution $\bar{V}_\epsilon$ in $D^{\frac{3}{2}-3\kappa,(0)3}$. Furthermore, bounds on the solution are uniform over $\epsilon \in [0, 1]$ and over $\zeta, \bar{\zeta}$ in bounded balls in their respective spaces.

Furthermore, for $\epsilon = 0$ and $\hat{\Pi}$ as in Theorem 4.8, the solution $\bar{V}$ is such that $\bar{v} = \mathcal{R}_0 \bar{V}$ solves

$$ \partial_t \bar{v} = \Delta \bar{v} + H'_{\eta}(u^{(0)})(\bar{v} + \mathcal{R}_0 \hat{V}^{(0)}), \quad (4.38) $$

with homogeneous boundary conditions, where $\hat{V}^{(0)}$ is as in Lemma 4.10. In particular, $\lim_{\epsilon \to 0} \hat{V}_\epsilon = \hat{V}$ is such that $\mathcal{R}_0 \hat{V}$ coincides with the process $v$ defined in (1.8).

**Remark 4.17.** In the case $d=1$, there is no condition on $\bar{\zeta}$ since the fixed point problem does not depend on it.
**Proof.** We first consider the case $d = 3$. Note first that, for any $C^4$ function $\tilde{G}$, we have $\mathcal{L}(\tilde{G}(u^{(0)}, u^{(1)})) \in D^{2,(0)3}$. Since $\hat{V}_\varepsilon^{(0)} \in D^{1,(-\frac{1}{2}-\kappa)(3)}$ by Lemma 4.10 and $\hat{V}_\varepsilon^{(1)} \in D^{\frac{3}{2}-3\kappa,(0)3}$ by Lemma 4.12, it then follows from [9, Lem. 4.3] that, for $\overline{V}_\varepsilon \in D^{\frac{3}{2}-3\kappa,(0)3}$, all the terms appearing after $P_{\text{Neu}}1_+^D$ in the right hand side of (4.35) belong to $D^{\frac{1}{2}-4\kappa,(-\frac{3}{2}-2\kappa)(3)}$, provided that $\kappa$ is sufficiently small. In particular, the operator $\mathcal{K}$ (defined as described in Remark 4.13) maps this continuously into $D^{2,(0)3}$, with arbitrarily small norm for small time intervals.

Furthermore, the reconstruction operator of Theorem C.5 continuously maps the space $D^{\frac{1}{2}-4\kappa,(-\frac{3}{2}-2\kappa)(3)}$ into $C^{-\frac{3}{2}-2\kappa}$, which is then mapped continuously into $C^{2,(0)3}$ by $K_\partial$ by Proposition A.2, and therefore into $D^{2,(0)3}$ by the Taylor lift $\mathcal{L}$, again with arbitrarily small norm for small time intervals as a consequence of the bound (A.3) which also holds for $K_\partial$.

Note now that by Corollary 5.6, we have $\mathbf{E}\|\xi_\varepsilon\|_{\kappa-2} \lesssim \varepsilon^{-\frac{1}{2}-\kappa}$. As a consequence of (4.21), we conclude from this that $\mathbf{E}|v_\varepsilon^{(0)}|_{L^\infty} \lesssim \varepsilon^{-\frac{1}{2}-\kappa}$. Since $\beta > \frac{1}{2}$, it follows from Lemma 4.15 that, for $p = (d+2)/c$,

$$ \| \hat{R}_\varepsilon^{(d)}(\mathcal{R}_\varepsilon \hat{V}_\varepsilon, \zeta) \|_{-c} \lesssim \| R_\varepsilon^{(d)}(\mathcal{R}_\varepsilon \hat{V}_\varepsilon + \varepsilon^{-\frac{1}{2}} u^{(1)}, \zeta) \|_{L^p} \lesssim \| \zeta \|_{L^p}, $$

uniformly over bounded sets for the model $\hat{\Pi}_\varepsilon$ and over bounded sets for $\overline{V}_\varepsilon + \hat{V}_\varepsilon^{(1)}$ in $D^{\frac{3}{2}-3\kappa,(0)3}$. Since $u^{(0)}$ and $u^{(1)}$ are bounded in $C^1$, it is immediate from (4.36) that one has a bound of the type

$$ \| \tilde{R}_\varepsilon^{(3)}(\bar{\zeta}) \|_{-\frac{1}{2}-2\kappa} \lesssim \| \bar{\zeta} \|_{-\frac{1}{2}-2\kappa}. $$

In particular, the argument of $P_{\text{Neu}}$ appearing in the last term on the right hand side of (4.35) is mapped continuously by $P_{\text{Neu}}$ into $C^{\frac{3}{2}-3\kappa,(0)3}$, again with arbitrarily small norm when considering a short enough time interval. Furthermore, all of these expressions are locally Lipschitz continuous (with similar bounds) as a function of $\overline{V}_\varepsilon$ in $D^{\frac{3}{2}-3\kappa,(0)3}$ and of the model $\hat{\Pi}_\varepsilon$, uniformly over $\varepsilon \in [0, 1]$ which yields the first claim over a short enough time (but bounded from below independently of $\varepsilon$) interval. This can be maximally extended as usual, and the claim follows from the fact that we know a priori that solutions to (4.38) do not explode.

The second claim is straightforward by simply setting $\varepsilon = 0$ and applying the reconstruction operator to both sides of (4.35). The case of $d = 2$ is virtually identical, noting in particular that even though $\mathcal{Q}_\varepsilon$ diverges in this case, it only does so logarithmically and is therefore compensated by the factor $\varepsilon^\kappa$ in (4.36). We leave the verification of the case $d = 1$ to the reader. $\square$

**Proposition 4.18.** Let $\zeta_\varepsilon = \varepsilon^\alpha \xi_\varepsilon$, $\bar{\zeta}_\varepsilon = \varepsilon^{2-d-\kappa} \sigma_\varepsilon$, and define $\hat{\Pi}_\varepsilon$ as in Section 4.3. Then, the assumptions of Proposition 4.16 are satisfied and we have $\zeta_\varepsilon, \bar{\zeta}_\varepsilon \to 0$ in their respective spaces. Furthermore, for any $\varepsilon > 0$, the modelled distribution $\hat{V}_\varepsilon$ constructed in Proposition 4.16 is such that $\mathcal{R}_\varepsilon \hat{V}_\varepsilon$ coincides with the process $v_\varepsilon$ defined in (4.1).

**Proof.** We first show that the assumptions of Proposition 4.16 are satisfied. The fact that the random models $\hat{\Pi}_\varepsilon$ are uniformly bounded (in probability) as $\varepsilon \to 0$
and converge in probability to $\hat{\Pi}$ is the content of Theorem 4.8. By the second part of Assumption 2.1 combined with stationarity, we furthermore see that

$$ \mathbf{E}\|\zeta_{\varepsilon}\|_{L_p}^p = \mathbf{E}\|\varepsilon^{\alpha}\xi_{\varepsilon}\|_{L_p}^p = \varepsilon^{p\alpha - \frac{(d+2)p}{2}} T\mathbf{E}|\eta(0)|^p \lesssim \varepsilon^{p/4}, $$

when $d \in \{2, 3\}$. For $d = 1$, we have

$$ \mathbf{E}\|\zeta_{\varepsilon}\|_{2\kappa - \frac{1}{2}} = \|\varepsilon^{(1/4)}\|_{2\kappa - \frac{1}{2}} \le \|\varepsilon^{(1/2-3\kappa)}\|_{2\kappa - \frac{1}{2}}, $$

which converges to 0 in probability by Corollary 5.6. We also conclude from Corollary 5.6 and our definitions that, for $d \in \{2, 3\}$, $\|\bar{\zeta}_{\varepsilon}\|_{-\frac{1}{2}-2\kappa} = \|\eta_{\varepsilon}^{(\frac{d}{2}-1+\kappa)}\|_{-\frac{1}{2}-2\kappa} \to 0$ in probability.

It remains to show that solutions coincide with $v_{\varepsilon}$. This is a special case of the general result obtained in [2] and could in principle also be obtained in a way similar to [14]. We present a short derivation here in order to remain reasonably self-contained.

The powercounting of the various symbols appearing in our structure depends on the dimension, so we first restrict ourselves to the case $d = 3$, which is the one with the largest number of terms of negative degree appearing. Combining (4.35) with Lemmas 4.10 and 4.12, we conclude that if we take for $\bar{V}_{\varepsilon}$ any solution to (4.35), there exist functions $v$ and $\nabla v$ such that, for $\Phi_{\varepsilon}^{(0)}$ and $\Phi_{\varepsilon}^{(1)}$ as in Lemma 4.12, the following identities hold on $\mathbf{R}_+ \times D$:

$$ \hat{V}_{\varepsilon}^{(0)} = G^{\circled{0}} + \Phi_{\varepsilon}^{(0)} \mathbf{1} + G^{\circled{0}} \nabla u^{(0)}^{\circled{0}}, $$

$$ \hat{V}_{\varepsilon}^{(1)} = \Phi_{\varepsilon}^{(1)} \mathbf{1} + GG^{\circled{0}} \circled{0} + G^{\circled{0}} \Phi_{\varepsilon}^{(0)} \circled{0}, $$

$$ \hat{V}_{\varepsilon} = G^{\circled{0}} + v \mathbf{1} + GG^{\circled{0}} \circled{0} + G^{\circled{0}} \nabla u^{(0)}^{\circled{0}} + G^{\circled{0}} u^{(1)}^{\circled{0}} \\ + \frac{G''G^2}{2} \circled{0}^{\circled{0}} + G^{\circled{0}} v^{\circled{0}} + \nabla v \mathbf{X}. $$

Developing the argument of $\mathcal{P}_{\text{Neu}}\mathbf{1}_{+}^D$ in (4.35) up to order 0, we conclude that it is given by

$$ \begin{aligned} H'_{\eta} G^{\circled{0}} + H'_{\eta} v \mathbf{1} + \frac{1}{2} G'' \nabla^2 u^{(0)} \Xi X^2 + G'(v - \Phi_{\varepsilon}^{(0)}) \circled{0} + G(G')^2 \circled{0}^{\circ} \\ + (G')^2 \nabla u^{(0)}^{\circled{0}} + (G')^2 u^{(1)}^{\circled{0}} + (G')^2 v^{\circled{0}} + G' \nabla v \circled{0} + \frac{1}{2} G^2 G' G''^{\circled{0}} \\ + GG'' \nabla u^{(0)}^{\circled{0}} + G'' v \nabla u^{(0)}^{\circled{0}} + G' u^{(1)}^{\circled{0}} + G' \nabla u^{(1)}^{\circled{0}} \\ + G'' \nabla u^{(0)} u^{(1)}^{\circled{0}} + \frac{1}{2} G^2 G''^{\circled{0}} + \frac{1}{2} G'' v^2 \circled{0} + \frac{1}{2} G''(u^{(1)})^2 \circled{0}^{\circ} + GG'' u^{(1)}^{\circled{0}} \\ + GG'' \nabla u^{(1)}^{\circled{0}} + GG'' v^{\circled{0}} + G' G''(u^{(1)})^2 \circled{0}^{\circ} + G^2 G' G''^{\circled{0}} + GG' G'' \nabla u^{(0)}^{\circled{0}} \\ + GG' G'' u^{(1)}^{\circled{0}} + G'' u^{(1)} v \circled{0} + GG' G'' u^{(1)}^{\circled{0}} + G' G'' u^{(1)} \nabla u^{(0)}^{\circled{0}}. \end{aligned} \quad (4.39) $$
At this point, we apply the results of [2]. Comparing [2, Eq. 2.20] with [2, Def. 3.20] and [2, Thm. 3.25], we see that each term appearing on the right hand side generates a counterterm for the renormalised equation. Each of these terms is of the form $\hat{F}(v, \nabla v, u^{(0)}, \nabla u^{(0)}, u^{(1)}, \nabla u^{(1)})\tau$ for some function $\hat{F}$ and some symbol $\tau$. The counterterm generated by any such term is then obtained precisely by simply replacing $\tau$ by the corresponding renormalisation constant and by interpreting the first two arguments of $\hat{F}$ as the value and gradient of the actual solution (after reconstruction).

**Remark 4.19.** One may worry that we are not quite in the framework of [2] because of the special treatment of $\hat{V}_\epsilon^{(0)}$ and $\hat{V}_\epsilon^{(1)}$. This however is due to purely analytical reasons that only affect the boundary behaviour. The computation of the renormalisation terms on the other hand is a purely algebraic affair which is not affected by this. The boundary conditions of $v_\epsilon$ however *are* affected by our decomposition and need to be determined separately.

It follows that, in dimension 3, the solution $v_\epsilon = \mathcal{R}_\epsilon \hat{V}_\epsilon$ to the fixed point problem with the renormalised model satisfies in $\mathbf{R}_+ \times D$ the PDE

$$ 
\begin{aligned}
\partial_t v_\epsilon ={}& \Delta v_\epsilon + H'_\eta v_\epsilon + G \xi_\epsilon + (v_\epsilon + \epsilon^{-1/2} u^{(1)}) G' \eta_\epsilon + \frac{1}{2} (v_\epsilon + \epsilon^{-1/2} u^{(1)})^2 G'' \sigma_\epsilon \\
& - \epsilon^{-3/2} G' G'' - \epsilon^{-1/2} G(G')^2 G''' - \epsilon^{-1/2} \langle \nabla u^{(0)}, (G')^2 \rangle + GG'' \sigma_\epsilon \\
& - (v_\epsilon + \epsilon^{-1/2} u^{(1)})(GG')' \frac{\epsilon^{-1/2}}{2} G^2 G''' - \frac{\epsilon^{-1}}{2} G^2 G'' \sigma_\epsilon \\
& - G''(G^2 G''' + GG' \langle \nabla u^{(0)} \rangle + GG' u^{(1)} \langle \nabla u^{(0)} \rangle) \sigma_\epsilon \\
& + \hat{R}_\epsilon^{(3)}(v_\epsilon, \zeta) + \tilde{R}_\epsilon^{(3)}(\bar{\zeta}).
\end{aligned}
$$

Furthermore, both $\mathcal{R}_\epsilon \hat{V}_\epsilon$ and $v_\epsilon^{(0)} = \mathcal{R}_\epsilon \hat{V}_\epsilon^{(0)}$ have homogeneous boundary conditions, so that the boundary conditions of $v_\epsilon$ coincide with those of $v_\epsilon^{(1)} = \mathcal{R}_\epsilon \hat{V}_\epsilon^{(1)}$.

By (4.36) and since we chose $\bar{\zeta} = \epsilon^{-1-\kappa}\sigma_\epsilon$, there is a cancellation between $\tilde{R}_\epsilon(\bar{\zeta})$ and some of the other terms appearing in this equation. Since furthermore $\circled{0}$, we obtain

$$ 
\begin{aligned}
\partial_t v_\epsilon ={}& \Delta v_\epsilon + H'_\eta v_\epsilon + G \xi_\epsilon + (v_\epsilon + \epsilon^{-1/2} u^{(1)}) G' \eta_\epsilon + \frac{1}{2} (v_\epsilon + \epsilon^{-1/2} u^{(1)})^2 G'' \sigma_\epsilon \\
& - \epsilon^{-3/2} G' G'' - \epsilon^{-1/2} G(G')^2 G''' - \epsilon^{-1/2} \langle \nabla u^{(0)}, (G')^2 \rangle \\
& - (v_\epsilon + \epsilon^{-1/2} u^{(1)})(GG')' \frac{\epsilon^{-1/2}}{2} G^2 G'' \sigma_\epsilon + \hat{R}_\epsilon^{(3)}(v_\epsilon, \zeta),
\end{aligned}
$$

which, when combining with the definition of $\Psi$ given in (1.6) and the fact that $\zeta = \epsilon^\alpha \xi_\epsilon$, precisely coincides with (4.2). Since its initial condition and boundary condition coincide as well, this completes the proof of the claim.

In dimension 2, a similar argument (but taking less terms into account) yields

$$ 
\begin{aligned}
\partial_t v_\epsilon ={}& \Delta v_\epsilon + H'_\eta v_\epsilon + G \xi_\epsilon + (v_\epsilon + u^{(1)}) G' \eta_\epsilon + \frac{1}{2} (v_\epsilon + u^{(1)})^2 G'' \sigma_\epsilon \\
& - \epsilon^{-1} G' G'' - G(G')^2 G''' - \langle \nabla u^{(0)}, (G')^2 \rangle + GG'' \sigma_\epsilon
\end{aligned}
$$
$$ -(v_{\varepsilon} + u^{(1)})(GG')^{\bigcirc} - \frac{1}{2}G^2G''^{\bigcirc} - \frac{1}{2}G^2G''^{\bigcirc\bigcirc} \sigma_{\varepsilon} + \hat{R}_{\varepsilon}^{(2)}(v_{\varepsilon}, \zeta) + \tilde{R}_{\varepsilon}^{(2)}(\bar{\zeta}). $$

Again, the term $\tilde{R}_{\varepsilon}^{(2)}$ precisely cancels the term proportional to $\bigcirc\bigcirc\sigma_{\varepsilon}$, so that this again coincides with (4.2). In dimension 1, an even simpler argument shows that

$$ \partial_t v_{\varepsilon} = \Delta v_{\varepsilon} + H'_{\eta} v_{\varepsilon} + G \xi_{\varepsilon} + v_{\varepsilon} G' \eta_{\varepsilon} + \frac{1}{2} v_{\varepsilon}^2 G'' \sigma_{\varepsilon} - \varepsilon^{-1/2} G' G^{\bigcirc} - v_{\varepsilon} (GG')^{\bigcirc} + \hat{R}_{\varepsilon}^{(1)}(v_{\varepsilon}, \zeta), $$

which again coincides with (4.2), noting that in this case one has $u^{(1)} = 0$. $\square$

## 5. Convergence of Models

In order to show convergence of the models, we apply the general result of [6, Thm 2.31]. This result shows that if one considers the “BPHZ lifts” of a sequence of smooth and stationary stochastic processes $\xi_n$ as given in [5, Thm 6.17] then, provided that one has uniform bounds of a suitable “norm” of $\xi_n$ and under a few relatively weak additional algebraic assumptions, the resulting sequence of models converges to a limit, provided that $\xi_n \to \xi$ weakly in probability.

### 5.1. Cumulant Homogeneity Assignments

In this section, we define

$$ \eta_{\varepsilon}^{(\alpha)}(t, x) = \varepsilon^{-\alpha} \eta(\varepsilon^{-2}t, \varepsilon^{-1}x), \qquad (5.1) $$

and we often use $z = (t, x)$ for space-time coordinates. The exponents $\alpha$ will always be chosen in $(c, \frac{d+2}{2}]$. Our aim is to obtain a suitable bound independent of $\varepsilon$ for joint cumulants of the form

$$ \kappa_p(\eta_{\varepsilon}^{(\alpha_1)}(z_1), \dots, \eta_{\varepsilon}^{(\alpha_p)}(z_p)). $$

Given a finite collection of at least two space-time points $z = \{z_a\}_{a \in A}$, we again consider the corresponding labelled binary tree $t_z = (T, n)$ as in Section 2.1, with the leaves of $T$ identified with the index set $A$. Recall that the nodes $V_T$ of $T$ are given by subsets of $A$, with inner nodes $\hat{V}_T$ given by subsets with at least two elements and the root of $T$ given by $A$ itself.

Recall from [6, Def. A.14] the following definition:

**Definition 5.1.** A “consistent cumulant homogeneity” consists, for each finite index set $A$, each binary tree $T$ over $A$ as above, and each choice of indices $\alpha: A \to [c, \frac{d+2}{2}]$, a function $c_T^{(\alpha)}: \hat{V}_T \to \mathbf{R}_+$ satisfying furthermore

* For every $B \subset A$, $\sum_{v \in \hat{V}_T: v \cap B \neq \emptyset} c_T^{(\alpha)}(v) \geq \sum_{a \in B} \alpha_a$.

* For every $u \in \hat{V}_T$, $\sum_{v \in \hat{V}_T: v \subset u} c_T^{(\alpha)}(v) \leq \sum_{a \in u} \alpha_a$.
• If $|A| \ge 3$, then for every $u \in \hat{V}_T$ with $|u| \le 3$, one has $\sum_{v \in \hat{V}_T: v \subset u} c_T^{(\alpha)}(v) < (d+2)(|u| - 1)$.

**Remark 5.2.** Applying the first two conditions with $B = A$ and $u = A$ respectively implies in particular that $\sum_{v \in \hat{V}_T} c_T^{(\alpha)}(v) = \sum_{a \in A} \alpha_a$.

We will now display a consistent cumulant homogeneity such that, for every finite set $A$, space-time points $\{z_a\}_{a \in A}$, and choices $\alpha: A \to [c, \frac{d+2}{2}]$, We have the bound

$$|\kappa_A(\{\eta_{\varepsilon}^{(\alpha_a)}(z_a)\}_{a \in A})| \lesssim \prod_{u \in \hat{V}_T} 2^{c_T^{(\alpha)}(u)\mathbf{n}(u)}, \quad (T, \mathbf{n}) = \mathbf{t}_z. \quad (5.2)$$

We claim that one possible choice is obtained by setting

$$c_T^{(\alpha)}(u) = \sum_{a \in A} \alpha_a 2^{-\vartheta(u,a)}, \quad (5.3)$$

where

$$\vartheta(u, a) = \begin{cases} |\{v \in \hat{V}_T : a \in v \subset u\}| & \text{if } a \in u \text{ and } u \neq A, \\ |\{v \in \hat{V}_T : a \in v\}| - 1 & \text{if } u = A, \\ +\infty & \text{if } a \notin u. \end{cases} \quad (5.4)$$

**Proposition 5.3.** *The choice (5.3) is a consistent cumulant homogeneity.*

*Proof.* The first two conditions of Definition 5.1 follow immediately from the structure of the formula (5.3), in particular the facts that $2^{-\vartheta(u,a)}$ is positive, vanishes for $a \notin u$, and is such that $\sum_{u \in \hat{V}_T} 2^{-\vartheta(u,a)} = 1$.

Regarding the last condition, the case $|u| = 3$ follows from the fact that the second condition holds and $\alpha_a \le \frac{d+2}{2}$. The case $|u| = 2$ follows from the condition $|A| \ge 3$ which guarantees that the corresponding sum is bounded by $\frac{d+2}{2}$ since $\vartheta(u, a) = 1$ in this case. $\square$

**Lemma 5.4.** *Setting $\underline{\alpha} = \inf_{a \in A} \alpha_a$ and $\bar{\alpha} = \sup_{a \in A} \alpha_a$, one has $c_T^{(\alpha)}(u) \in [\underline{\alpha}, \bar{\alpha}]$ for $u \in \hat{V}_T \setminus \{A\}$ and $c_T^{(\alpha)}(A) \in [2\underline{\alpha}, 2\bar{\alpha}]$.*

*Proof.* By convexity of (5.3), it suffices to consider the case where $\alpha_a = 1$ for all $a$. We proceed by induction on the size of $A$. When $|A| = 2$, one has $\hat{V}_T = \{A\}$, so that $\vartheta(A, a) = 0$ and therefore $c_T^{(\alpha)}(A) = 2$ as claimed.

Assume now that $|A| > 2$ and fix a binary tree $T$ over $A$. Write $A_1$ and $A_2$ for the children of $A$ in $\hat{V}_T$, so that $A = A_1 \sqcup A_2$. We distinguish two cases. In the first case, $|A_1| \wedge |A_2| \ge 2$, so that the tree $T$ can naturally be thought of as two trees $T_1, T_2$ over $A_1, A_2$ joined by their roots. By (5.4) (in particular the fact that there is an additional $-1$ at the root), we then have

$$c_T^{(\alpha)}(A_1) = \frac{1}{2} c_{T_1}^{(\alpha)}(A_1), \quad c_T^{(\alpha)}(A_2) = \frac{1}{2} c_{T_2}^{(\alpha)}(A_2),$$
$$c_T^{(\alpha)}(A) = \frac{1}{2}(c_{T_1}^{(\alpha)}(A_1) + c_{T_2}^{(\alpha)}(A_2)),$$

while $c_T^{(\alpha)}(u) = c_{T_1}^{(\alpha)}(u)$ for all other $u \in \hat{V}_T$, with $i \in \{1, 2\}$ depending on whether $u \subset A_1$ or $u \subset A_2$. We conclude by using the induction hypothesis, which implies that $c_{T_1}^{(\alpha)}(A_1) = c_{T_2}^{(\alpha)}(A_2) = 2$.

In the second case, we have $|A_1| = 1$ and $|A_2| \ge 2$ (or vice-versa), the case $|A| = 2$ having already been dealt with. In this case, the tree $T$ consists of a subtree $T_2$ over $A_2$ as before, with an additional root vertex $A$ and single extra leaf. In this case, we have

$$c_T^{(\alpha)}(A_2) = \frac{1}{2}c_{T_2}^{(\alpha)}(A_2), \quad c_T^{(\alpha)}(A) = 1 + \frac{1}{2}c_{T_2}^{(\alpha)}(A_2),$$

whence we conclude as before.

□

**Proposition 5.5.** *Under Assumption 2.1, the bound (5.2) holds for the choice (5.3).*

*Proof.* It follows from Assumption 2.1 and (5.1) that

$$|\kappa_A(\{\eta_{\varepsilon}^{(\alpha_a)}(z_a)\}_{a \in A})| \lesssim \varepsilon^{c_{\varepsilon}(\mathbf{n}(A)) - \sum_{a \in A} \alpha_a} 2^{c_{\varepsilon}(\mathbf{n}(A))\mathbf{n}(A)} \prod_{u \in \hat{V}_T} \varepsilon^{c_{\varepsilon}(\mathbf{n}(u))} 2^{c_{\varepsilon}(\mathbf{n}(u))\mathbf{n}(u)}, \quad (5.5)$$

where

$$c_{\varepsilon}(n) = \begin{cases} \frac{c}{\bar{c}} & \text{if } n \ge \log_2 \frac{1}{\varepsilon}, \\ \bar{c} & \text{otherwise.} \end{cases} \quad (5.6)$$

Let now $\hat{c}: \hat{V}_T \to [\underline{c}, \bar{c}]$ be any map such that $\hat{c}(A) + \sum_{u \in \hat{V}_T} \hat{c}(u) = \sum_{a \in A} \alpha_a$, and rewrite (5.5) as

$$|\kappa_A(\{\eta_{\varepsilon}^{(\alpha_a)}(z_a)\}_{a \in A})| \lesssim \varepsilon^{c_{\varepsilon}(\mathbf{n}(A)) - \hat{c}(A)} 2^{c_{\varepsilon}(\mathbf{n}(A))\mathbf{n}(A)} \prod_{u \in \hat{V}_T} \varepsilon^{c_{\varepsilon}(\mathbf{n}(u)) - \hat{c}(u)} 2^{c_{\varepsilon}(\mathbf{n}(u))\mathbf{n}(u)}. \quad (5.7)$$

We now note that (5.6) implies that

$$\varepsilon^{c_{\varepsilon}(n)} - \hat{c} 2^{c_{\varepsilon}(n)n} = (\varepsilon 2^n)^{c_{\varepsilon}(n)} - \hat{c} 2^{\hat{c}n} \le 2^{\hat{c}n},$$

for every $\hat{c} \in [\underline{c}, \bar{c}]$, uniformly over $n \in \mathbb{Z}$ and $\varepsilon \in (0, 1]$. Inserting this into (5.7) immediately yields that, uniformly in $\varepsilon$, one has

$$|\kappa_A(\{\eta_{\varepsilon}^{(\alpha_a)}(z_a)\}_{a \in A})| \lesssim 2^{\hat{c}(A)\mathbf{n}(A)} \prod_{u \in \hat{V}_T} 2^{\hat{c}(u)\mathbf{n}(u)}. \quad (5.8)$$

Since the map $c_T^{(\alpha)}$ is of the desired type by Lemma 5.4 [modulo the additional factor 2 at the root which is taken care of explicitly in (5.8)], the claim follows. □
**Corollary 5.6.** For any $\alpha \in (c, \frac{d+2}{2})$, one has $\eta_{\varepsilon}^{(\alpha)} \to 0$ in probability in $C^{\beta}$ for every $\beta < -\alpha$. In particular, $\eta_{\varepsilon} \to 0$ in probability in $C^{\beta}$ for every $\beta < -1$ and $\zeta_{\varepsilon} \to 0$ in probability in $C^{\beta}$ for every $\beta < \frac{d-2}{2}$. The same holds for $\eta_{\varepsilon}^{(\alpha)} \mathbf{1}_A$ for any fixed Borel set $A$.

*Proof.* The first statement is an immediate consequence of [6, Thm 2.31]. The fact that we can multiply $\eta_{\varepsilon}^{(\alpha)}$ by an arbitrary indicator function follows from the fact that these bounds do not involve the derivative of the test function in this case. $\square$

## 5.2. Power-Counting Conditions

By (5.3)–(5.4), the quantity $\mathcal{H}_{c,D}(A)$ defined in [6, Def. A.24] can be estimated by

$$ \mathcal{H}_{c,D}(A) = \sum_{a \in D} 2^{-\partial(D,a)} \alpha_a \geq 2^{-(|D|-1)} \sum_{a \in D} \alpha_a. \quad (5.9) $$

This is because for $D \in \mathring{V}_T$ with $D \neq A$, one always has $\partial(D, a) \leq |D| - 1$. Furthermore, by [6, Rem. 2.28], one has $j_D(A) \geq \frac{d+2}{2} - \kappa$ for all typed sets $A$ and $D$. These observations allow us to conclude the following:

**Lemma 5.7.** All the decorated trees generated by the rule $R_{\circ}$ of Definition 4.4 are c-super regular in the terminology of [6, Def. A.27].

**Remark 5.8.** Since our scaling and degree assignment are fixed throughout and since we consider all cumulants, that is we choose $\mathcal{L}_{\text{Cum}}$ to contain all possible cumulants, we omit these from our notation.

*Proof.* We reduce ourselves to considering symbols of negative degree since the claim then follows automatically for the remaining ones. These symbols are listed in (4.15). Note also that the definition of c-super regularity is non-trivial only for trees that contain at least three noises, so that it suffices to consider the symbols

$$ \mathring{\circ}, \mathring{\circ}_{\omega}, \mathring{\circ}_{\omega'}, \mathring{\circ}_{\omega''}, \mathring{\circ}_{\omega'''}, \mathring{\circ}_{\omega''''}, \mathring{\circ}_{\omega'''''}, \mathring{\circ}_{\omega''''''}, \mathring{\circ}_{\omega'''''''} \quad (5.10) $$

By (5.9), it is sufficient to verify that for every subtree $\tau$ containing $k \geq 2$ instances of a noise one has the bound

$$ \deg \tau + \left( \frac{d+2}{2} \wedge 2^{-(k-1)} \sum_{i=1}^{k} |\deg \Xi_i| \right) - \kappa \geq 0, $$

where $\Xi_i$ denotes the $i$th noise appearing in $\tau$. This can be seen simply by inspection of the list (5.10). $\square$
### 5.3. Special Bounds

For some of the symbols in our regularity structure, we will bounds that are stronger than what is suggested by the degrees of the symbols in question. In this statement, $\phi$ denotes an arbitrary measurable function with

$$ \sup_z |\phi(z)| \le 1, \quad \supp\phi \subset B(0, 1). \qquad (5.11) $$

**Proposition 5.9.** Let $\delta > 0$ be as in Section 2. For $\alpha \in [-\delta, 1]$ and $\beta = -\frac{1}{2} - \alpha$, we have the bounds

$$ (\hat{\Pi}_{\epsilon} \circ_{\mathcal{O}_u}^{\circ}) (\phi_z^{\lambda}) \lesssim \epsilon^{\alpha} \lambda^{\beta}, \quad (\hat{\Pi}_{\epsilon} \circ_{\mathcal{O}_u}^{\circ}) (\phi_z^{\lambda}) \lesssim \epsilon^{\frac{1}{2}+\alpha} \lambda^{\beta}, \quad (\hat{\Pi}_{\epsilon} \circ_{\mathcal{O}_u}) (\phi_z^{\lambda}) \lesssim \epsilon^{1+\alpha} \lambda^{\beta}, $$

where we write $X \lesssim Y$ as a shorthand for the existence, for every $p \ge 1$, of a constant $C$ such that $E|X|^p \le C Y^p$, uniformly over all $\lambda \le 1$, $\epsilon \le 1$ and $\phi$ as in (5.11).

*Proof.* Recall that, by Proposition 5.5, we can view any $\eta_{\epsilon}^{(\theta)}$ as in (5.1) for $\theta \in [\frac{1}{2}-\delta, \frac{5}{2}]$ as a “noise” of regularity $-\theta$ whose “norm”, as measured by [6, Def. A.18] with respect to the cumulant homogeneity just described remains uniformly bounded as $\epsilon \to 0$.

In particular, we can write

$$ \Pi_{\epsilon} \circ_{\mathcal{O}_u}^{\circ} = \epsilon^{\alpha} (K * \eta_{\epsilon}^{(2)})^2 \eta_{\epsilon}^{(\frac{1}{2}+\alpha)}, $$

and we can apply [6, Thm 2.31], showing that the BPHZ renormalisation of this term satisfies the required bounds. Recall that $\hat{\Pi}_{\epsilon} \circ_{\mathcal{O}_u}^{\circ}$ doesn't quite agree with the BPHZ renormalisation, but the error between the two is given by $2\epsilon^{3/2}\delta g_{\epsilon}(\circ_{\mathcal{O}_u}) \Pi_{\epsilon} \circ$ with $\delta g_{\epsilon}$ as in (4.19), which is easily seen to satisfy the required bounds.

The other two terms can be dealt with similarly by writing

$$ \Pi_{\epsilon} \circ_{\mathcal{O}_u}^{\circ} = \epsilon^{\frac{1}{2+\alpha}} (K * \eta_{\epsilon}^{(2)}) \eta_{\epsilon}^{(\frac{1}{2}+\alpha)}, \quad \Pi_{\epsilon} \circ_{\mathcal{O}_u} = \epsilon^{1+\alpha} \eta_{\epsilon}^{(\frac{1}{2}+\alpha)}, $$

thus concluding the proof. $\square$

**Corollary 5.10.** For every $\kappa > 0$, every $\alpha \in [0, 1)$, and every Borel set $A \subset \mathbb{R}^{d+1}$, one has the bound $(E\|1_A \xi_{\epsilon}\|_{C^w}^p)^{1/p} \lesssim \epsilon^{-\alpha}$ with $w = \alpha - \frac{5}{2} - \kappa$.

*Proof.* This follows from the third bound of Proposition 5.9 by Kolmogorov's criterion, using the fact that $\hat{\Pi}_{\epsilon} \circ_{\mathcal{O}_u} = \epsilon^3 \xi_{\epsilon}$ and that we can move the multiplication by $1_A$ onto the test function. $\square$
### 5.4. Tightness and Convergence for the Noise

In this section, we show that the convergence announced in Theorem 4.8 holds. As usual, convergence is obtained by first showing tightness and then identification of the limiting distribution. More precisely, we prove the following:

**Theorem 5.11.** One has $\xi_{\varepsilon} \to \xi$ weakly in $C^{\alpha}$ for every $\alpha < -\frac{d+2}{2}$.

It turns out that this statement is a relatively straightforward consequence of the following proposition. Writing $\kappa_p(X)$ for the $p$th cumulant of a real-valued random variable $X$, one has the following:

**Proposition 5.12.** For $p \ge 2$, we have the bound

$$|\kappa_p(\langle \eta, \phi_0^\lambda \rangle)| \lesssim (\lambda^{-\bar{c}p} + \lambda^{-(d+2)(p-1)}) \wedge \lambda^{-\varepsilon p},$$

uniformly over all $\lambda \in \mathbf{R}_+$ and all $\phi$ as in (5.11).

Before we turn to the proof of Proposition 5.12, let us show how to deduce Theorem 5.11 from it. First, we have the following corollary:

**Corollary 5.13.** For any $\kappa \le \frac{1}{2}$ and $p \ge 2$, one has the bounds

$$\mathbf{E}|\langle \xi_{\varepsilon}, \phi_0^{\lambda} \rangle|^p \lesssim \lambda^{-\frac{d+2}{2}p},$$

uniformly over $\lambda \le 1$, $\varepsilon \le 1$ and $\phi$ as in (5.11).

*Proof.* By simple rescaling, it is straightforward to see that the required bound is equivalent to the bound

$$\mathbf{E}|\langle \eta, \phi_0^{\bar{\lambda}} \rangle|^p \lesssim \bar{\lambda}^{-\frac{d+2}{2}p}, \quad (5.12)$$

uniformly over $\bar{\lambda} \le 1/\varepsilon$ and $\varepsilon \le 1$.

For even integers $p \ge 2$, Proposition 5.12 implies that $|\kappa_p(\langle \eta, \phi_0^\lambda \rangle)| \lesssim \lambda^{-\frac{d+2}{2}p} \wedge \lambda^{-\varepsilon p}$, so that.

$$\mathbf{E}|\langle \eta, \phi_0^{\bar{\lambda}} \rangle|^p \lesssim \lambda^{-\frac{d+2}{2}p} \wedge \lambda^{-\varepsilon p}, \quad (5.13)$$

It remains to observe that (5.13) implies that $\mathbf{E}|\langle \eta, \phi_0^{\bar{\lambda}} \rangle|^p \lesssim \bar{\lambda}^{-\bar{c}p}$, uniformly over all $\bar{\lambda} > 0$, for any $c \in [c, \frac{d+2}{2}]$. Since (5.12) is of this form and our assumptions guarantee that the values of $c$ appearing there fall into the correct interval, this concludes the proof. □

*Proof of Theorem 5.11.* Using [11, Eq. 10.4] (which is nothing but an analogue of Kolmogorov's continuity criterion) and the compactness of the embeddings $C^{\alpha} \subset C^{\beta}$ for $\alpha > \beta$ (over any bounded domain), it follows from Corollary 5.13 that the laws $\xi_{\varepsilon}$ are tight in $C^{\alpha}$ for every $\alpha < -\frac{d+2}{2}$.

It remains to show that every limit $\xi$ of a convergent subsequence of $\xi_{\varepsilon}$ is space-time white noise. It thus suffices to show that, for every $\phi \in C_0^{\infty}$ with $\int \phi^2(z) dz =$
1, $\xi(\phi)$ is centred normal with variance 1. Since all moments of $\langle \xi_{\epsilon}, \phi \rangle$ remain bounded as $\epsilon \to 0$ by Corollary 5.13, one has

$$ \begin{aligned} \kappa_p(\xi(\phi)) &= \lim_{\epsilon \to 0} \kappa_p(\langle \xi_{\epsilon}, \phi \rangle) = \lim_{\epsilon \to 0} \epsilon^{-\frac{(d+2)p}{2}} \kappa_p(\langle \eta, \phi_0^{1/\epsilon} \rangle) \\ &\lesssim \lim_{\epsilon \to 0} \epsilon^{-\frac{(d+2)p}{2}} (\epsilon^{\bar{c}p} + \epsilon^{(d+2)(p-1)}). \end{aligned} $$

For $p \ge 3$, this vanishes, thus showing that $\xi(\phi)$ is Gaussian. It clearly has zero mean since this is already the case for $\langle \xi_{\epsilon}, \phi \rangle$. Its variance is given by

$$ \mathbf{E}|\langle \xi_{\epsilon}, \phi \rangle|^2 = \epsilon^{-(d+2)} \int \kappa_2(z, \bar{z}) \phi_0^{1/\epsilon}(z) \phi_0^{1/\epsilon}(\bar{z}) \, dz \, d\bar{z}. $$

Note now that, for every $\epsilon > 0$, every $\delta \in [0, 1]$ and any $z, \bar{z} \in \mathbf{R}^{d+1}$, one has the bound

$$ |\phi_0^{1/\epsilon}(z) - \phi_0^{1/\epsilon}(\bar{z})| \lesssim \epsilon^{d+2} (\epsilon \|z - \bar{z}\|)^{\delta}. $$

Recalling that $\int \phi^2 = 1$ and that $\int \kappa_2(0, z) \, dz = 1$ by (1.3), we thus obtain

$$ \begin{aligned} |\mathbf{E}|\langle \xi_{\epsilon}, \phi \rangle|^2 - 1 &= |\epsilon^{-(d+2)} \int \kappa_2(z, \bar{z}) \phi_0^{1/\epsilon}(z) (\phi_0^{1/\epsilon}(\bar{z}) - \phi_0^{1/\epsilon}(z)) \, dz \, d\bar{z}| \\ &\lesssim \epsilon^{\delta} \int |\kappa_2(z, \bar{z}) \phi_0^{1/\epsilon}(z)| \|z - \bar{z}\|^{\delta} \, dz \, d\bar{z} \\ &\lesssim \epsilon^{\delta} \int (\|z'\|^{-2\bar{c}} \wedge \|z'\|^{-2\epsilon}) \|z'\|^{\delta} \, dz' \lesssim \epsilon^{\delta}. \end{aligned} $$

Here, we used the fact that the integral of $\phi_0^{1/\epsilon}$ is independent of $\epsilon$, as well as our assumptions (1.3) and (2.1) on the covariance function of $\eta$. In particular, we choose $\delta$ as in the definition of $\bar{c}$, see Section 2, whence $2\bar{c} - \delta > d + 2$ which guarantees integrability at infinity. (Integrability at 0 is guaranteed by $2\underline{c} - \delta = d - \delta < d + 2$.) $\square$

**Proof of Proposition 5.12.** We expand the expression for the cumulant as

$$ \kappa_p(\langle \eta, \phi_0^{\lambda} \rangle) = \int \kappa_p(z_1, \dots, z_p) \phi_0^{\lambda}(z_1) \dots \phi_0^{\lambda}(z_p) \, dz_1 \dots dz_p. \quad (5.14) $$

In order to bound this integral, we use a simplified version of the type of multiscale analysis used in [15, 16]. Let us recall how this works.

We now write $\mathbf{t} = (T, \mathbf{n})$ for a generic binary tree, together with a scale assignment as above (that is we enforce the fact that $\mathbf{n}$ is monotone) and $\mathbf{T}$ for the set of all such $\mathbf{t}$. Given $\mathbf{t} \in \mathbf{T}$, we write $V_{\mathbf{t}}/E_{\mathbf{t}}$ for the vertex/edge set of the corresponding tree and $\mathbf{n}_{\mathbf{t}}$ for the scale assignment. In particular, $\mathbf{n}_{\mathbf{t}}(\Omega_{\mathbf{p}})$ denotes the scale assignment for the root, which controls the diameter of the set $\{z_1, \dots, z_{\mathbf{p}}\}$ in $\mathbf{R}^{d+1}$. The number $p$ will always be considered fixed, so we do not include it explicitly in our notation. Denoting by $\mathcal{T}$ the map $\mathcal{T}: (z_1, \dots, z_{\mathbf{p}}) \mapsto \mathbf{t} \in \mathbf{T}$ defined at the start of Section 2, we set $D_{\mathbf{t}} = \mathcal{T}^{-1}(\mathbf{t})$ for the set of all configurations of points
$z \in (\mathbf{R}^{d+1})^p$ giving rise to a given combinatorial data. Then, it was shown in [15, Lem. A.13] that, for every bounded Borel set $U \in \mathbf{R}^{d+1}$ one has the bound

$$|D_t \cap \{z : z_1 \in U\}| \lesssim |U| \prod_{A \in \hat{V}_t} 2^{-(d+2)\mathbf{n}_t(A)},$$

where $|\cdot|$ denotes Lebesgue measure. Furthermore, by construction, $\bigcup_{t \in T} D_t$ is of full Lebesgue measure.

Without loss of generality, we can restrict ourselves to the case where the support of $\phi$ has diameter bounded by 1 in the parabolic distance. With all of this notation at hand, we then bound (5.14) by

$$|\kappa_p((\eta, \phi_0^\lambda))| \lesssim \lambda^{-(d+2)p} \sum_{t \in T} |D_t \cap \{z : z_i \in \text{supp} \phi_0^\lambda \ \forall i\}| \sup_{z \in D_t} |\kappa_p(z_1, \dots, z_p)|.$$

We simplify this expression as follows. First, we note that

$$D_t \cap \{z : z_i \in \text{supp} \phi_0^\lambda \ \forall i\} = \emptyset$$

as soon as $2^{-n(\Omega_p)} \ge \lambda$, since the support of $\phi_0^\lambda$ is bounded by $\lambda$, so that we can restrict the sum above to those $t$ satisfying $2^{-n(\Omega_p)} \le \lambda$. Furthermore, one has

$$|D_t \cap \{z : z_i \in \text{supp} \phi_0^\lambda \ \forall i\}| \le |D_t \cap \{z : z_1 \in \text{supp} \phi_0^\lambda \}| \lesssim \lambda^{d+2} \prod_{A \in \hat{V}_t} 2^{-(d+2)\mathbf{n}_t(A)}.$$

Combining this with Assumption 2.1, we conclude that

$$|\kappa_p((\eta, \phi_0^\lambda))| \lesssim \lambda^{(d+2)(1-p)} \sum_{t \in T} \mathbf{1}_{2^{-n(\Omega_p)} \le \lambda} 2^{c(\Omega_p) \mathbf{n}_t(\Omega_p)} \prod_{A \in \hat{V}_t} 2^{(c(A)-(d+2))\mathbf{n}_t(A)}.$$

We treat separately the cases $\lambda \le 1$ and $\lambda \ge 1$. In the former case, we can apply [15, Lem. A.10] with the distinguished vertex $v_*$ appearing there equal to the root $\Omega_p$. The first condition appearing there is then satisfied by the fact that $2_c - (d+2) < 0$ by assumption, while the second condition is empty and therefore trivially satisfied. Note that the shape $T$ of the tree is fixed in [15, Lem. A.10], while we also sum over all possible shapes, but since there are finitely many of them for any fixed $p$, this just yields an additional prefactor. We thus obtain the bound

$$\begin{aligned} |\kappa_p((\eta, \phi_0^\lambda))| & \lesssim \lambda^{(d+2)(1-p)} \sum_{t \in T} \mathbf{1}_{2^{-n(\Omega_p)} \le \lambda} 2^{c \mathbf{n}_t(\Omega_p)} \prod_{A \in \hat{V}_t} 2^{(c-(d+2))\mathbf{n}_t(A)} \\ & \lesssim \lambda^{(d+2)(1-p)-c} \prod_{A \in \hat{V}_t} \lambda^{d+2-c}. \end{aligned}$$

Since $\hat{V}_t$ contains exactly $p-1$ elements (the tree $T$ is binary and has $p$ leaves), we finally obtain

$$|\kappa_p((\eta, \phi_0^\lambda))| \lesssim \lambda^{-pc},$$

as required.
The case $\lambda > 1$ is split into two subcases. If $p \ge 3$, we use the fact that $|\kappa_p(z)| \lesssim \bar{\kappa}_{p,\bar{c}}(z)$ uniformly over all $z$, so that, using [15, Lem. A.10] as above, we have

$$|\kappa_p(\langle \eta, \phi_0^\lambda \rangle)| \lesssim \lambda^{(d+2)(1-p)} \sum_{\mathbf{t} \in \mathcal{T}} \mathbf{1}_{2^{-n}(\Omega_p) \le \lambda} 2^{\bar{c} \cdot \mathbf{n}_t(\Omega_p)} \prod_{A \in \hat{V}_\mathbf{t}} 2^{(\bar{c} - (d+2)) \cdot \mathbf{n}_t(A)} \lesssim \lambda^{-p\bar{c}}.$$

Note that, in order to be able to apply that result, we need to verify that, for every subtree $\hat{\mathbf{t}}$ of $\mathbf{t}$ spanned by some subset of its leaves satisfies $\sum_{A \in \hat{V}_\mathbf{t}} \alpha(A) < 0$, where $\alpha(A) = 2\bar{c} - (d+2)$ if $A = \Omega_p$ and $\alpha(A) = \bar{c} - (d+2)$ otherwise. This is of course trivially satisfied as soon as $\hat{\mathbf{t}} \ne \mathbf{t}$ since $\bar{c} < d+2$ by assumption. The exponent at the root however is given by $2\bar{c} - (d+2)$, which is positive, but since $p \ge 3$, the sum of all exponents is given by $2\bar{c} - (d+2) + (p-2)(\bar{c} - (d+2)) = p\bar{c} - (d+2)(p-1)$, which is indeed negative for $p \ge 3$ and $\delta < \frac{d+2}{6}$, which we assume w.l.o.g..

It remains to consider the case $p = 2$ and $\lambda > 1$. In this case, the above computation reduces to

$$|\kappa_2(\langle \eta, \phi_0^\lambda \rangle)| \lesssim \lambda^{-(d+2)} \left( \sum_{1 \le 2^{-n} \le \lambda} 2^{(2\bar{c} - (d+2))n} + \sum_{2^{-n} < 1} 2^{(2\bar{c} - (d+2))n} \right) \lesssim \lambda^{-(d+2)},$$

as required, thus concluding the proof. $\square$

<figure><img src="image_6.png" /></figure>

We also use the following convergence results which do not strictly speaking follow from Theorem 5.11 since multiplication by an indicator function (even that of a hypercube) is not a continuous operation on $C^\beta$ for $\beta < 0$:

**Lemma 5.14.** Let $\xi_{\varepsilon,D} = \xi_{\varepsilon} \mathbf{1}_{\mathbf{R} \times D}$ and $\xi_{\varepsilon,D}^{+} = \xi_{\varepsilon} \mathbf{1}_{\mathbf{R}_{+} \times D}$. Then, $\xi_{\varepsilon}, \xi_{\varepsilon,D}$, and $\xi_{\varepsilon,D}^{+}$ jointly weakly converge to limits $\xi, \xi_D$ and $\xi_D^+$.

*Proof.* The proof is identical to that of Theorem 5.11, using the fact that all bounds we used are uniform over test functions as in (5.11), so that multiplying them by the indicator function of some domain changes nothing. $\square$

<figure><img src="image_7.png" /></figure>

### 5.5. Boundary Term

Recall that we have set

$$\hat{\Pi}_{\varepsilon} \diamond, = \Pi_{\varepsilon} \diamond, - \sum_{i=1}^{d} \varepsilon^{1-\frac{d}{2}} (c_{i,0} \delta_{\partial_{i,0} D} + c_{i,1} \delta_{\partial_{i,1} D}), \quad (5.15)$$

where $\partial_{i,0} D = \{x \in D : x_i = 0\}$ and $\partial_{i,1} D = \{x \in D : x_i = 1\}$ and the constants are given by

$$c_{i,j} = \int_{\mathbf{R}_+} Q_{i,j}(s) \, ds, \quad (5.16)$$
where the function $Q_{i,j}$ is defined as follows. For $i = 1, \dots, d$, write $\iota_i : \mathbf{R} \times \mathbf{R}^d \to \mathbf{R}^{d+1}$ for the map given by

$$ \iota_i^{-1}(t, x) = (x_i, v), \quad v = (t, x^{(i)}), $$

where $x^{(i)} \in \mathbf{R}^{d-1}$ denotes the vector obtained from $x$ by deleting the $i$th coordinate. With this notation, we then set

$$ Q_{i,0}(s) = \int_{\mathbf{R}_+ \times \mathbf{R}^d} (P \circ \iota_i)(s + \beta, v) [(\kappa_2 \circ \iota_i)(s - \beta, v) - (\kappa_2 \circ \iota_i)(s + \beta, v)] \, d\beta \, dv, $$

$$ Q_{i,1}(s) = \int_{\mathbf{R}_+ \times \mathbf{R}^d} (P \circ \iota_i)(s + \beta, v) [(\kappa_2 \circ \iota_i)(\beta - s, v) - (\kappa_2 \circ \iota_i)(-s - \beta, v)] \, d\beta \, dv. $$

**Remark 5.15.** We will show in Lemma 5.21 that both the integrands in the definition of $Q_{i,j}$ and the functions $Q_{i,j}$ themselves are integrable, so that these expressions are all finite.

**Remark 5.16.** The formula given above is valid for the case of Neumann boundary conditions. In the case of Dirichlet boundary conditions, a similar formula holds, but the precise values of the constants do not matter since they do not affect the solutions.

Note that although $P \circ \iota_i \neq P$, it does not depend on $i$, while $\kappa_2 \circ \iota_i$ does depend on $i$ in general, since we do not assume that the driving noise is isotropic. We henceforth write $\hat{P} = P \circ \iota_i$. One of the main results of this section is that the renormalised model on $\diamond$ vanishes in a suitable sense as $\varepsilon \to 0$. We first provide a bound on its expectation, which requires the bulk of the work. For the formulation of this result, we write $B_0^1$ for the set of all test functions $\phi \in C_0^1$ with support contained in the parabolic ball of radius 1 and such that $\max\{\|\phi\|_\infty, \|D\phi\|_\infty\} \le 1$.

**Proposition 5.17.** With $\hat{\Pi}_\varepsilon \diamond$ defined as in (5.15), one has for $d \in \{2, 3\}$ and any $\kappa > 0$ small enough,

$$ |\mathbf{E}(\hat{\Pi}_\varepsilon \diamond)(\phi_z^\lambda)| \lesssim \varepsilon^\kappa \lambda^{-\frac{d}{2}-\kappa}, \qquad (5.17) $$

uniformly over $\varepsilon, \lambda \le 1$, $z \in \mathbf{R}^{d+1}$, and $\phi \in B_0^1$.

Before we turn to the proof, we introduce a number of notations and preliminary bounds. Write $\mathcal{G}$ for the reflection group generated by “elementary” reflections across the $2d$ planes containing the faces of $D$. The group $\mathcal{G}$ is naturally identified with $\mathbf{Z}^d$ (as a set, not as a group!) since for each $k \in \mathbf{Z}^d$ there exists exactly one element $R_k \in \mathcal{G}$ mapping $k+D$ into $D$. We also write $\mathcal{G} \ni R \mapsto (-1)^R \in \{-1, 1\}$ for the group homomorphism mapping the elementary reflections to $-1$. We will write $\Lambda : \mathbf{R}^{d+1} \to \mathbf{R} \times D$ for the map such that $\Lambda \upharpoonright \mathbf{R} \times (k+D) = \text{id} \times \mathcal{R}_k$ and $S : \mathbf{R}^{d+1} \to \{-1, 1\}$ by $S \upharpoonright \mathbf{R} \times (k+D) = (-1)^{R_k}$.
With these notations, it follows from Proposition A.1 that the truncated Neumann and Dirichlet heat kernels are such that for $z \in \mathbf{R} \times D$ one has the identities

$$ 
\begin{aligned}
\int_{\mathbf{R} \times D} K_{\text{Neu}}(z, z') f(z') \, dz' &= \int_{\mathbf{R}^{d+1}} K(z - z') f(\Lambda(z')) \, dz', \\
\int_{\mathbf{R} \times D} K_{\text{Dir}}(z, z') f(z') \, dz' &= \int_{\mathbf{R}^{d+1}} S(z') K(z - z') f(\Lambda(z')) \, dz'.
\end{aligned} 
\quad (5.18) $$

(Note that $S$ is ill-defined on the measure zero set consisting of the reflection planes, but since it always appears in an integral this does not matter.)

**Lemma 5.18.** Let $U \subset \mathbf{R}^{d+1}$ and let $\Phi: U \to \mathbf{R}^{d+1}$ be a diffeomorphism between $U$ and its image such that $\|D\Phi^{-1}\|$ is bounded uniformly over $\Phi(U)$. Assume furthermore that $z$ is such that $\|z - z'\| \ge \lambda$ for all $z' \in U$. Then,

$$ \int_U K(z - z') |\kappa_2^{(\varepsilon)}(z, \Phi(z'))| \, dz' \lesssim \varepsilon^{\frac{d}{2}+\kappa} \lambda^{-\frac{d}{2}-\kappa}. $$

for all $\kappa \in [0, \frac{d}{2}]$.

*Proof.* We write $U = U_1 \cup U_2$ where

$$ U_1 = \{z' \in U : \|z - z'\| \ge \|z - \Phi(z')\|\}. $$

Since $\|z - z'\| \ge \lambda$ for $z' \in U$ by assumption, we have the bound $K(z - z') \lesssim \lambda^{-d}$. On $U_1$, it follows from the definition that $K(z - z') \lesssim \|z - \Phi(z')\|^{-d}$. As a consequence, we obtain the bound

$$ 
\begin{aligned}
\int_{U_1} K(z - z') |\kappa_2^{(\varepsilon)}(z, \Phi(z'))| \, dz' &\lesssim \int_{\Phi(U_1)} (\|z - z'\|^{-d} \wedge \lambda^{-d}) |\kappa_2^{(\varepsilon)}(z, z')| \, dz' \\
&\lesssim \int_{\mathbf{R}^{d+1}} (\|z'\|^{-d} \wedge (\lambda/\varepsilon)^{-d}) \rho^2(\|z'\|) \, dz' \lesssim 1 \wedge \varepsilon^d \lambda^{-d} \\
&\le \varepsilon^{\frac{d}{2}+\kappa} \lambda^{-\frac{d}{2}-\kappa},
\end{aligned}
 $$

as claimed. On $U_2$ on the other hand, we use the fact that $\rho$ is a decreasing function, so that

$$ 
\begin{aligned}
\int_{U_2} K(z - z') |\kappa_2^{(\varepsilon)}(z, \Phi(z'))| \, dz' &\lesssim \varepsilon^{-2} \int_{U_2} (\|z - z'\|^{-d} \wedge \lambda^{-d}) \rho^2(\|S_{\varepsilon}(z - \Phi(z'))\|) \, dz' \\
&\lesssim \varepsilon^{-2} \int_{\mathbf{R}^{d+1}} (\|z - z'\|^{-d} \wedge \lambda^{-d}) \rho^2(\|S_{\varepsilon}(z - z')\|) \, dz',
\end{aligned}
 $$

which is then bounded exactly as above. $\square$

<figure><img src="image_8.png" /></figure>
**Lemma 5.19.** For any fixed $c > 0$, one has

$$ \int_{\|z\| \ge c} |\kappa_2^{(\varepsilon)}(0, z)| \, dz \lesssim \varepsilon^{d+2\delta}, $$

uniformly over $\varepsilon \in (0, 1]$.

*Proof.* We can assume that $\varepsilon < c$, so that

$$ |\kappa_2^{(\varepsilon)}(0, z)| \lesssim \varepsilon^{d+2+2\delta} \|z\|^{-d-2-2\delta} $$

by Assumption 2.1, and the bound follows at once.

<figure><img src="image_9.png" /></figure>

*Proof of Proposition 5.17.* We now consider the Neumann case, the Dirichlet case follows from a virtually identical calculation. We start by bounding the expectation of $(\hat{\Pi}_\varepsilon)_{\varepsilon, \varepsilon}(\phi_\varepsilon^\lambda)$. By the reflection principle (5.18), the correction due to the Neumann boundary conditions is given by

$$ \int_{\mathbf{R} \times D} K_{\text{Neu}}(z, z') \xi_{\varepsilon}(z') \, dz' - \int_{\mathbf{R}^{d+1}} K(z - z') \xi_{\varepsilon}(z') \, dz' \\ = \int_{\mathbf{R}^{d+1}} K(z - z') (\xi_{\varepsilon}(\Lambda z') - \xi_{\varepsilon}(z')) \, dz'. $$

It then follows from the definitions that

$$ \mathbf{E}(\hat{\Pi}_{\varepsilon})_{\varepsilon}(\psi) = \varepsilon^{-\frac{d}{2}} \int_{\mathbf{R} \times D} \psi(z) \int_{\mathbf{R}^{d+1}} K(z - z') (\kappa_2^{(\varepsilon)}(z, \Lambda z') - \kappa_2^{(\varepsilon)}(z, z')) \, dz' \, dz. \quad (5.19) $$

Writing $\lambda = \lambda_\psi$ for the diameter of $\text{supp}\psi$, we aim to give a bound of the form

$$ |\mathbf{E}(\hat{\Pi}_{\varepsilon})_{\varepsilon}(\psi)| \lesssim \varepsilon^{\kappa} \lambda^{-\frac{d}{2}-\kappa} (\lambda^{d} \|\psi\|_{\infty} + \lambda^{d+1} \|D_{x}\psi\|_{\infty}), \quad (5.20) $$

for any $\kappa > 0$ small enough, uniformly over $\varepsilon, \lambda \le 1$. We restrict ourselves to the case $\lambda_\psi \le \frac{1}{8}$ and we set

$$ W_{\psi} = \{z : \exists z' \in \text{supp}\psi \quad \|z - z'\| \le \lambda_{\psi}\}. $$

We furthermore assume for the moment that

$$ W_{\psi} \subset \mathbf{R} \times [-3/4, 3/4]^d \quad (5.21) $$

and that, whenever $W_\psi \cap \partial D = \emptyset$, one has $W_\psi \subset D$. We will see later that it is always possible to reduce oneself to this case. Finally, for $x \in \mathbf{R}^d$, we write $|x|_0$ for the number of vanishing coordinates of $x$ and we set $d_\psi = \sup\{|x|_0 : (t, x) \in W_\psi\}$. We treat the different cases separately.

**The case** $d_\psi = 0$. In this case $W_\psi \cap \partial D = \emptyset$, so that $\text{supp}\psi \cap \partial D = \emptyset$ and in particular $(\hat{\Pi}_\varepsilon)_{\varepsilon, \varepsilon}(\psi) = (\Pi_\varepsilon)_{\varepsilon, \varepsilon}(\psi)$. We then have

$$ |\mathbf{E}(\hat{\Pi}_{\varepsilon})_{\varepsilon}(\psi)| \le \varepsilon^{-\frac{d}{2}} \int_{\mathbf{R} \times D} |\psi(z)| \int_{\mathbf{R} \times D^c} K(z - z') (|\kappa_2^{(\varepsilon)}(z, \Lambda z')| $$
$$+ |\kappa_2^{(\varepsilon)}(z, z')|) \, dz' \, dz. \qquad (5.22)$$

We break the inner integral into a finite sum of integrals over $\mathbf{R} \times (k + D)$, since $K$ has compact support and $z \in \mathbf{R} \times D$. Since we can restrict $z$ to the support of $\psi$, we have $|x - x'| \ge \lambda$ for all $x' \in k + D$ with $k \ne 0$ by the definition of $W_\psi$. We can then apply Lemma 5.18 with $U = \mathbf{R} \times (k + D)$ and $\Phi = \Lambda$ for the first term, while $\Phi = \text{id}$ for the second term. It follows that $|\mathbf{E}(\hat{\Pi}_\varepsilon \mathfrak{g})(\psi)|$ is bounded by

$$\varepsilon^\kappa \lambda^{-\frac{d}{2}-\kappa} \int |\psi(z)| \, dz, \qquad (5.23)$$

which is indeed bounded by the right hand side of (5.20).

**The case** $d_\psi = 1$. In this case, the support of $\psi$ is located near one of the faces of $D$ (say $\partial_{i,0}D$), but its distance to the other faces is at least $\lambda$. Write $\mathfrak{R}$ for the element of $G$ which corresponds to reflection around the plane containing $\partial_{i,0}D$ and $\pi_i : \mathbf{R}^{d+1} \to \mathbf{R}^{d+1}$ for the orthogonal projection onto that plane. We also write $E_i = \{(t, x) : x_i < 0\}$.

We first note that

$$\mathbf{E}(\Pi_{\varepsilon} \mathfrak{g})(\psi) = \varepsilon^{-\frac{d}{2}} \int_{\mathbf{R} \times D} \psi(z) \int_{E_i} K(z-z') (\kappa_2^{(\varepsilon)}(z, \mathfrak{R}z') - \kappa_2^{(\varepsilon)}(z, z')) \, dz' \, dz + \mathcal{R}, \qquad (5.24)$$

where $|\mathcal{R}|$ is bounded by (5.23). Indeed, the integrands in (5.19) and (5.24) vanish on $D$ and coincide on $\mathfrak{R}D$. Their integral over the complement of these two regions is then bounded exactly as before by applying Lemma 5.18 to finitely many translates of $D$. By Lemma 5.19, we can furthermore replace $K$ by $P$ so that

$$\mathbf{E}(\Pi_{\varepsilon} \mathfrak{g})(\psi) = \varepsilon^{-\frac{d}{2}} \int_{\mathbf{R} \times D} \psi(z) \int_{E_i} P(z-z') (\kappa_2^{(\varepsilon)}(z, \mathfrak{R}z') - \kappa_2^{(\varepsilon)}(z, z')) \, dz' \, dz + \tilde{\mathcal{R}}, \qquad (5.25)$$

where $|\tilde{\mathcal{R}}|$ is bounded by (5.23). The reason for this is that $P - K$ is supported outside of an annulus of radius 1 and $\psi$ is supported at a distance of at most 1/4 from the reflection plane of $\mathfrak{R}$, and one has $\|z - \mathfrak{R}z'\| \ge 1/2$ for all $z, z'$ with $z \in \text{supp}\psi$ and $\|z - z'\| \ge 1$.

On the other hand, we claim that

$$\varepsilon^{1-\frac{d}{2}} c_{i,0} \delta_{\partial_{i,0}D}(\psi) = \varepsilon^{-\frac{d}{2}} \int_{\mathbf{R} \times D} \psi(\pi_i z) \int_{E_i} P(z-z') (\kappa_2^{(\varepsilon)}(z, \mathfrak{R}z') - \kappa_2^{(\varepsilon)}(z, z')) \, dz' \, dz.$$

To see that this is the case, we first perform the change of variables $z = \iota_i(s, v)$ and similarly for $z'$, and we note that $\pi_i \iota_i(s, v) = \iota_i(0, v)$, so that the right hand side is given by

$$\varepsilon^{-\frac{d}{2}} \int_{\mathbf{R}_+ \times \mathbf{R}^d} \psi(\iota_i(0, v)) \int_{\mathbf{R}_- \times \mathbf{R}^d} P(\iota_i(s - s', v - v'))$$
$$
\begin{align*}
& (\kappa_2^{(\varepsilon)}(\iota_i(s+s', v-v')) - \kappa_2^{(\varepsilon)}(\iota_i(s-s', v-v'))) \, d(s', v') \, d(s, v) \\
&= \varepsilon^{-\frac{d}{2}} \int_{\mathbf{R}^d} \psi(\iota_i(0, v)) \, dv \int_{\mathbf{R}_+} \int_{\mathbf{R}_- \times \mathbf{R}^d} P(\iota_i(s-s', v')) \\
& \qquad (\kappa_2^{(\varepsilon)}(\iota_i(s+s', v')) - \kappa_2^{(\varepsilon)}(\iota_i(s-s', v'))) \, d(s', v') \, ds \\
&= \varepsilon^{1-\frac{d}{2}} \int_{\mathbf{R}^d} \psi(\iota_i(0, v)) \, dv \int_{\mathbf{R}_+} \int_{\mathbf{R}_- \times \mathbf{R}^d} P(\iota_i(s-s', v')) \\
& \qquad (\kappa_2(\iota_i(s+s', v')) - \kappa_2(\iota_i(s-s', v'))) \, d(s', v') \, ds
\end{align*}
$$

Performing the substitution $s' \mapsto -\beta$ and comparing to the definition of $Q_{i,0}$, we
conclude that

$$
\mathbf{E}(\hat{\mathbf{\Pi}}_{\varepsilon} \diamond) (\psi) = \varepsilon^{-\frac{d}{2}} \int_{\mathbf{R}_+} Q_{i,0} \left(\frac{s}{\varepsilon}\right) \int_{\mathbf{R}^d} \left( (\psi \circ \iota_i)(s, u) - \psi \circ \iota_i)(0, u) \right) du \, ds + \tilde{R}. \quad (5.26)
$$

By Lemma 5.21, the function $Q_{i,0}$ satisfies the bounds

$$
\sup_{s \in \mathbf{R}_+} |Q_{i,0}(s)| < \infty, \quad \int_{\mathbf{R}_+} |Q_{i,0}(s)| \, ds < \infty, \quad \int_{\mathbf{R}_+} |Q_{i,0}(s)| \, s \, ds < \infty.
$$

Note now that the integral over $u$ is restricted to the projection of the support of $\psi$
which is of volume at most $\lambda^{d+1}$ (since $u$ consists of $d-1$ spatial variables and
one time variable). Bounding $\psi$ by its supremum, it follows that

$$
\mathbf{E}(\hat{\mathbf{\Pi}}_{\varepsilon} \diamond) (\psi) \lesssim \varepsilon^{1-\frac{d}{2}} \lambda^{d+1} \|\psi\|_{\infty} \int_{\mathbf{R}_+} |Q_{i,0}(s)| \, ds + |R|.
$$

On the other hand, we can bound $|(\psi \circ \iota_i)(s, u) - \psi \circ \iota_i)(0, u)|$ by $|s| \|D_x \psi\|_{\infty}$,
which similarly yields

$$
\mathbf{E}(\hat{\mathbf{\Pi}}_{\varepsilon} \diamond) (\psi) \lesssim \varepsilon^{2-\frac{d}{2}} \lambda^{d+1} \|D_x \psi\|_{\infty} \int_{\mathbf{R}_+} |Q_{i,0}(s)| s \, ds + |R|.
$$

Combining these and choosing $\psi = \phi_z^\lambda$ so that $\|\psi\|_\infty \lesssim \lambda^{-d-2}$ and $\|D_x\psi\|_\infty \lesssim \lambda^{-d-3}$ yields a bound of the order $\varepsilon^{-\frac{d}{2}}(\varepsilon\lambda^{-1} \wedge \varepsilon^2\lambda^{-2})$, which does imply (5.17) since $-\frac{d}{2} - \kappa \in (-2, -1)$.

The case $d_{\psi} = 2$. We claim that this case can be obtained as a consequence of the bounds for the cases $d_{\psi} = 0$ and $d_{\psi} = 1$. We consider the case of dimension $d = 2$ (but the computation below is done in a way that keeps track of dimension and applies to $d = 3$ as well) so that by (5.21)

$$
W_{\psi} \cap \{(t, x) : x_1 = x_2 = 0\} \neq \emptyset.
$$

We then fix a smooth function $\chi : \mathbf{R} \to \mathbf{R}_+$ such that $\text{supp}\chi \subset [0, 2]$ and such that $\sum_{k \in \mathbb{Z}} \chi_k(x) = 1$, where $\chi_k(x) = \chi(x-k)$. For some fixed constant $c > 1$ and integers $n, k_i$ and $\ell$, we then set

$$
\chi_{n,k,l}(z) = \chi_n(-\log|x|/\log 2) \chi_{k_1}(c2^n x_1) \chi_{k_2}(c2^n x_2) \chi_l(c^2 2^{2n} t).
$$
By choosing $c$ sufficiently large, we can guarantee that, for every $n, k$ and $\ell$, the function $\chi_{n,k,\ell}$ is such that $d_{\chi_{n,k,\ell}} \in \{0, 1\}$. Furthermore, there are only finitely many values of $k$ (independently of $n$ and $\ell$) for which $\chi_{n,k,\ell} \neq 0$. This is because $\chi_{n,k,\ell}(S_{2^n}z)$ is independent of $n$. Fix now a test function of the form $\psi = \phi_{z_0}^\lambda$ and write

$$ \psi_{n,k,\ell}(z) = \phi_{z_0}^{\lambda}(z) \chi_{n,k,\ell}(z). $$

By construction, one has $\psi_{n,k,\ell} = 0$ for $n$ such that $2^{-n} \ge 2\lambda$, so that

$$ \lambda \psi_{n,k,\ell} \le 2^{-n}, \quad \|D^m \psi_{n,k,\ell}\|_{\infty} \lesssim \lambda^{-d-2} 2^{n|m|}. $$

Applying the bounds we already obtained for $d\psi \in \{0, 1\}$, we conclude that

$$ |\mathbf{E}(\hat{\mathbf{\Pi}}_{\varepsilon})_{\mathcal{S}_l}(\psi_{n,k,\ell})| \lesssim \varepsilon^{\kappa} 2^{n(\frac{d}{2}+\kappa)} \lambda^{-d-2} 2^{-n(d+2)}. $$

For any given $n$, the number of values for $k$ and $\ell$ leading to non-vanishing $\lambda\psi_{n,k,\ell}$ is of the order of $(\lambda 2^n)^d$, so that we eventually obtain the bound

$$ \begin{aligned} |\mathbf{E}(\hat{\mathbf{\Pi}}_{\varepsilon})_{\mathcal{S}_l}(\psi)| & \lesssim \varepsilon^{\kappa} \sum_{2^{-n} \ge 2\lambda} (\lambda 2^n)^d 2^{n(\frac{d}{2}+\kappa)} \lambda^{-d-2} 2^{-n(d+2)} \\ &= \varepsilon^{\kappa} \lambda^{-2} \sum_{2^{-n} \ge 2\lambda} 2^{n(\frac{d}{2}+\kappa-2)} \lesssim \varepsilon^{\kappa} \lambda^{-\frac{d}{2}-\kappa}, \end{aligned} $$

as claimed. The case of dimension $d=3$ is identical, the only difference being that $\ell$ now has two components. Note that this calculation breaks in $d=4$ where the sum over $n$ diverges. This suggests that in this case one would have to add to $\hat{\mathbf{\Pi}}_{\varepsilon}$ an additional correction term that charges faces of codimension 2.

**The case** $d\psi = 3$. This is relevant only for $d=3$, we shall however keep track of $d$ in our calculation to illustrate how this would behave in higher dimensions. We then proceed in the same way as for the case $d\psi = 2$, making use this time of the fact that we already have the required bound for all test functions with $d\psi < 3$. This time, we have

$$ W_{\psi} \cap \{(t, x) : x_1 = x_2 = x_3 = 0\} \neq \emptyset, $$

and we set in a fashion similar to above,

$$ \chi_{n,k,\ell}(z) = \chi_n(-\log|x|/\log 2)\chi_{k_1}(c2^n x_1)\chi_{k_2}(c2^n x_2)\chi_{k_3}(c2^n x_3)\chi_{\ell}(c^2 2^{2n}t). $$

This time, for any given $n$, the number of values for $k$ and $\ell$ leading to non-vanishing $\lambda\psi_{n,k,\ell}$ is of the order of $(\lambda 2^n)^{d-1}$, which then yields similarly to above

$$ \begin{aligned} |\mathbf{E}(\hat{\mathbf{\Pi}}_{\varepsilon})_{\mathcal{S}_l}(\psi)| & \lesssim \varepsilon^{\kappa} \sum_{2^{-n} \ge 2\lambda} (\lambda 2^n)^{d-1} 2^{n(\frac{d}{2}+\kappa)} \lambda^{-d-2} 2^{-n(d+2)} \\ &= \varepsilon^{\kappa} \lambda^{-3} \sum_{2^{-n} \ge 2\lambda} 2^{n(\frac{d}{2}+\kappa-3)} \lesssim \varepsilon^{\kappa} \lambda^{-\frac{d}{2}-\kappa}. \end{aligned} $$

Note that this time the series actually converges for all $d < 6$.
To conclude, we justify the assumption (5.21) and the fact that, for $W_\psi \cap \partial D = \emptyset$, one has $W_\psi \subset D$. Indeed, by our assumption on $\lambda_\psi$, it is always possible to enforce this by applying a finite number of reflections around the planes $\{x : x_i = 1/2\}$. If $\text{supp}\psi$ intersects $\partial_{i,1}D$ for example, then we reflect around $x_i = 1/2$ to have $\text{supp}\psi$ intersect $\partial_{i,0}D$ instead. The only effect of this reflection is that, in order to obtain the same answer, we only need to reflect the noise $\eta$ around that plane. The effect of this operation on its covariance function is to change the sign of the $i$th spatial coordinate of its argument, which is why we obtain $Q_{i,1}$ rather than $Q_{i,0}$ in (5.26). $\square$

We now have the main ingredients in place to prove the main result of this section.

**Theorem 5.20.** With $\hat{\Pi}_\varepsilon \diamond$, defined as in (5.15), one has for $d \in \{2, 3\}$ and any $\varepsilon > 0$ small enough,

$$ \| (\hat{\Pi}_\varepsilon \diamond)(\phi_\varepsilon^\lambda) \|_{L^p} \lesssim \varepsilon^\kappa \lambda^{-\frac{d}{2}-\kappa}, \quad (5.27) $$

uniformly over $\varepsilon, \lambda \le 1$, $z \in \mathbf{R}^{d+1}$, and $\phi \in \mathcal{B}_0^1$.

*Proof.* Writing $\psi = \phi_\varepsilon^\lambda$, the triangle inequality yields

$$ \| (\hat{\Pi}_\varepsilon \diamond)(\psi) \|_{L^p} \le \| (\hat{\Pi}_\varepsilon \diamond)(\psi) - \mathbf{E}(\hat{\Pi}_\varepsilon \diamond)(\psi) \|_{L^p} + |\mathbf{E}(\hat{\Pi}_\varepsilon \diamond)(\psi)|, $$

and we already obtained the required bound on the second term in Proposition 5.17, so we focus on the first one. Furthermore, $\hat{\Pi}_\varepsilon \diamond$, differs from $\Pi_\varepsilon \diamond$, by a deterministic quantity, so that we only need to bound

$$ \| (\Pi_\varepsilon \diamond)(\psi) - \mathbf{E}(\Pi_\varepsilon \diamond)(\psi) \|_{L^p}. $$

By (4.7) combined with (5.18), this random variable equals

$$ \varepsilon^{-\frac{d}{2}} \int_{\mathbf{R} \times D} \psi(z) \int_{\mathbf{R}^{d+1}} K(z-z') (\varepsilon_\varepsilon(z)\varepsilon_\varepsilon(\Lambda z') - \varepsilon_\varepsilon(z)\varepsilon_\varepsilon(z')) dz' dz (5.28) $$

This time, we will not need to exploit the cancellation between these two terms on $\mathbf{R} \times D$, so we simply bound both terms separately. The second term equals $(\hat{\Pi}_\varepsilon \diamond)(\psi)$, which is bounded by the right hand side of (5.27) by Theorem 4.8.

For the first term, we use the fact that $K$ is compactly supported, so that it can be bounded by a finite sum of terms of the type

$$ \varepsilon^{-\frac{d}{2}} \int_{\mathbf{R} \times D} \psi(z) \int_{\mathbf{R} \times D} K(z - Rz') : \varepsilon_\varepsilon(z) \varepsilon_\varepsilon(z') : dz' dz, $$

with $R \in \mathcal{G}$.

The expectation of the $p$th power of this expression is given by a multiple integral with the integrand given by a sum of terms, each of which is a product of heat kernels and of cumulants. At this stage, we note that the bound in [6] does not exploit any further cancellations, so we can put absolute values everywhere, bound $K(z - Rz')$ by $\|z - z'\|^{-d}$, and use the bounds from that paper. $\square$
5.6. Bounds on the Function Q

It remains to prove

**Lemma 5.21.** *Under Assumption 2.1,*

$$
\sup_{s \in \mathbf{R}_+} |Q_{i,0}(s)| + \int_{\mathbf{R}_+} |Q_{i,0}(s)| (1+s) \, ds < \infty,
$$

and similarly for $Q_{i,1}$.

*Proof.* We fix $i$ and we simply write $Q$ instead of $Q_{i,0}$. In all the estimates below we will use repeatedly the inequality $|(\kappa_2 \circ t_i)(s-\beta, v)| \ge |(\kappa_2 \circ t_i)(s+\beta, v)| \lesssim t^{-\beta} \wedge t^{-\bar{c}}$.
We first estimate $\sup_{s \in \mathbf{R}_+} |Q(s)|$. Here $\beta \in \mathbf{R}_+$, $v = (t, x)$, with $x \in \mathbf{R}^{d-1}$.

$$
\begin{align*}
|Q(s)| & \lesssim \int_0^\infty \int_0^\infty \int_{\mathbf{R}^{d-1}} \hat{P}(s+\beta, v) t^{-\bar{c}} \wedge t^{-\bar{c}} \, d\beta \, dv \\
& \lesssim \int_0^\infty \int_0^\infty \frac{e^{-(s+\beta)^2/(4t)}}{\sqrt{t}} t^{-\bar{c}} \wedge t^{-\bar{c}} \, d\beta \, dt \\
& \lesssim \int_s^\infty \int_0^\infty \frac{e^{-1/t}}{\sqrt{t}} |a| \left( (a^2 t)^{-\bar{c}} \wedge (a^2 t)^{-\bar{c}} \right) da \, dt \\
& \lesssim \int_0^\infty \frac{e^{-1/t}}{t^{1-\delta}} \, dt \int_0^{t^{-1/2}} a^{2\delta} \, da + \int_0^\infty \frac{e^{-1/t}}{t^{\bar{c}+1/2}} \, dt \int_{t^{-1/2}}^\infty \frac{da}{a^{2\bar{c}-1}} \\
& \lesssim \int_0^\infty \frac{e^{-1/t}}{t^{3/2}} \, dt + \int_0^\infty \frac{e^{-1/t}}{t^{3/2}} \, dt < \infty,
\end{align*}
$$

and the bound does not depend upon $s$. Here, to go from the second to the third line, we set $a = s + \beta$ and we performed the substitution $t \mapsto a^2 t/4$.

It remains to estimate $\int_{\mathbf{R}_+} |Q(s)| (1+s) \, ds$. Again $v = (t, x)$, with $x \in \mathbf{R}^{d-1}$.

$$
\begin{align*}
\int_{\mathbf{R}_+} |Q(s)| (1+s) \, ds & \lesssim \int_0^\infty \int_0^\infty \int_{\mathbf{R}^d} (1+s) \hat{P}(s+\beta, v) t^{-\bar{c}} \wedge t^{-\bar{c}} \, dv \, d\beta \, ds \\
& \lesssim \int_0^\infty \int_{\mathbf{R}^d} (a+a^2) \hat{P}(a, v) t^{-\bar{c}} \wedge t^{-\bar{c}} \, dv \, da \\
& \lesssim \int_0^\infty \int_0^\infty (a^2+a^3) \frac{e^{-1/t}}{\sqrt{t}} (a^2 t)^{-\bar{c}} \wedge (a^2 t)^{-\bar{c}} \, dt \, da \\
& \lesssim \int_0^\infty \frac{e^{-1/t}}{t^{1-\delta}} \, dt \int_0^{t^{-1/2}} (a^{1+2\delta} + a^{2+2\delta}) \, da \\
& \phantom{\lesssim \int_0^\infty \frac{e^{-1/t}}{t^{1-\delta}} \, dt \int_0^{t^{-1/2}} (a^{1+2\delta} + a^{2+2\delta}) \, da} + \int_0^\infty \frac{e^{-1/t}}{t^{\bar{c}+1/2}} \, dt \int_{t^{-1/2}}^\infty (a^{2-2\bar{c}} + a^{3-2\bar{c}}) \, da \\
& \lesssim \int_0^\infty \left( \frac{1}{t^2} + \frac{1}{t^{5/2}} \right) e^{-1/t} \, dt \\
& \phantom{\lesssim \int_0^\infty \left( \frac{1}{t^2} + \frac{1}{t^{5/2}} \right) e^{-1/t} \, dt} + \int_0^\infty \left( \frac{1}{t^2} + \frac{1}{t^{5/2}} \right) e^{-1/t} \, dt < \infty,
\end{align*}
$$

thus concluding the proof.

□
## 6. Auxiliary Results

In this section, we collect a number of results that are more or less straightforward consequences of known results, specialised to our setting. Throughout this section, we assume that we are working with the regularity structure defined in Section 4.2 and that $\xi_{\epsilon}$ is defined as in (4.1) and satisfies Assumption 2.1.

We will write $C^{\gamma}$ as a shorthand for $D^{\gamma}(\bar{T})$, where $\bar{T}$ is the sector spanned by the Taylor polynomials, and similarly for $C^{\gamma, \eta}$, etc. Note that for $\gamma \notin \mathbb{N}$ this is consistent with the usual definition of $C^{\gamma}$.

**Proposition 6.1.** Let $\zeta^{+}, \zeta \in C^{-\frac{5}{2}-\kappa}$ be such that $\zeta^{+}(\phi) = 0$ for $\phi$ supported in $\{t < 0\}$ and $\zeta^{+}(\phi) = \zeta(\phi)$ for $\phi$ supported in $\{t > 0\}$ and let $\Pi$ be an admissible model with $\Pi \diamond = \zeta$. Write $K1_{+} \diamond = K\hat{\zeta}1_{+} \diamond$ for $K\hat{\zeta}$ defined as in [9, Sec. 4.5]. Define $K_{\text{Neu}}1_{+} \diamond$, $K_i1_{+} \diamond$ and $K_{\text{Neu},i}1_{+} \diamond$ analogously and, given $\Phi$, let $V^{(0)}$ be given by (4.22) for some $u^{(0)} \in C^{3-\kappa}$.

Then, setting $\bar{\gamma} = \frac{3}{2} - \kappa$ and $\bar{\eta} = -\frac{1}{2} - \kappa$, there exists a choice $\Phi \in C^{\bar{\gamma}, \bar{\eta}}$ such that $RV^{(0)} = K_{\text{Neu}}(G(u^{(0)})\zeta^{+})$ and $V^{(0)}$ belongs to $D^{\bar{\gamma}, \bar{\eta}}$. In particular, if $\zeta^{+} \in C^{\alpha}$ for some $\alpha > -1$ is supported in $R_{+} \times D$ then, for $t \in [0, 1]$, $RV^{(0)}$ coincides with the solution of $\partial_t v = \Delta v + G(u^{(0)})\zeta^{+}$ with vanishing initial condition, endowed with Neumann (respect. Dirichlet) boundary conditions.

If furthermore $\zeta_n \to \zeta$ in $C^{-\frac{5}{2}-\kappa}$ and $\Pi_n \to \Pi$ as admissible models, then one has $\|V_n^{(0)}\|_V V^{(0)} \|\bar{\gamma}, \bar{\eta} - 0$.

**Remark 6.2.** In principle, the model $\Pi$ does contain non-trivial information through its action on $\hat{\gamma}$. This is because the kernel $K_{\theta}$ is not 2-regularising (condition 5.5 in [11, Ass. 5.1] fails to be satisfied), so that the extension theorem [11, Thm 5.14] cannot be applied here.

*Proof.* We aim to apply Corollary B.6. Let $\gamma = \frac{1}{2} - 2\kappa$, $\eta = -\frac{5}{2} - \kappa$ (so that in particular $\gamma - \eta = 3 - \kappa$), and let $B = C^{\gamma-\eta}$ (on which $C^{\gamma-\eta}$ then acts canonically by multiplication) with the injection $\iota g = (L_{\gamma-\eta}g)1_{+} \diamond \in D^{\gamma, \eta}$. (This is actually independent of the model $\Pi$.) Given $g \in B$, we set

$$ \hat{R}_g = g \zeta^+, $$

which is consistent with the reconstruction operator by our assumption on $\zeta^+$. We are therefore in the setting of Corollary B.6 provided that we set $K_0 = K_{\text{Neu}}$.

This guarantees that we can find $\Phi \in C^{\gamma+2, \eta+2}$ with the desired properties. The continuity as a function of $\zeta^+$ and the model $\Pi$ follows from the corresponding continuity statement in Corollary B.6. $\square$

**Proposition 6.3.** For every $g \in C^{2-\kappa}$ one can find $\Phi$ taking values in the Taylor polynomials such that, setting

$$ V = 1_{+}^{D} g(\hat{\Pi}_{\epsilon})^{\flat} + \Phi, $$

one has $V \in D^{2-2\kappa, \bar{w}}$ with $\bar{w} = (\frac{1}{2} - 2\kappa, \frac{1}{2} - 2\kappa, 0)$, and $R_{\epsilon}V = K(1_{+}g\hat{\Pi}_{\epsilon})^{\flat}$.
**Proof.** We make use of [9, Lem. 4.12] and Corollary B.6. For this, similarly to above, we set $\mathcal{B} = \mathcal{C}^{2-\kappa}$ and, for $g \in \mathcal{B}$, we set

$$ \iota g \stackrel{\text{def}}{=} \mathbf{1}_{+}^{\mathcal{D}} \mathcal{L}_{2-\kappa} (g \hat{\mathbf{\Pi}}_{\varepsilon} \mathbf{1}) \circ, \quad \hat{\mathcal{R}} g \stackrel{\text{def}}{=} \mathbf{1}_{+} g \hat{\mathbf{\Pi}}_{\varepsilon} \mathbf{1}. $$

Note that as a consequence of Proposition A.2 one has $\hat{\mathbf{\Pi}}_{\varepsilon} \mathbf{1} \in \mathcal{C}^{\tilde{\gamma}, \tilde{w}}$ for any $\tilde{\gamma} > 0$ and for $\tilde{w} = (-\frac{1}{2}-\kappa)_3$, where $(\eta)_3 \stackrel{\text{def}}{=} (\eta, \eta, \eta)$. Since $\deg \circ = -1-\kappa$, it follows from [9, Lem. 4.3] that $\iota g \in \mathcal{D}^{\tilde{\gamma}, w}$ for $\gamma = 1-2\kappa$ and for $w = (-\frac{3}{2}-2\kappa)_3$. It then follows from [9, Lem. 4.12] that one can find a modelled distribution $V \in \mathcal{D}^{2-2\kappa, \tilde{w}}$ with $\tilde{w} = (\frac{1}{2}-2\kappa, \frac{1}{2}-2\kappa, 0)$ of the form

$$ V = \mathbf{1}_{+}^{\mathcal{D}} (g \hat{\mathbf{\Pi}}_{\varepsilon} \mathbf{1}) \circ + \hat{\Phi}, $$

with $\hat{\Phi}$ taking values in the Taylor polynomials and such that

$$ \mathcal{R}_{\varepsilon} V = K \hat{\mathcal{R}} g = K (\mathbf{1}_{+} g \hat{\mathbf{\Pi}}_{\varepsilon} \mathbf{1}), $$

thus concluding the proof.

□

**Proposition 6.4.** Equip the regularity structure $(\mathcal{T}, \mathcal{G})$ with an admissible model $\Pi$ in the sense of Definition 4.2 and let $G$ be a modelled distribution of the form

$$ G = (\mathcal{L}(g_1) \circ + \mathcal{L}(g_{2,i}) K_i(\circ) + \mathcal{L}(g_3)), $$

for some functions $g_1, g_2 \in \mathcal{C}^{2-\kappa}$ and $g_3 \in \mathcal{C}^{\frac{3}{2}-\kappa, w}$ with $w = (-\frac{1}{2}-\kappa, \frac{1}{2}-\kappa, \frac{1}{2}-\kappa)$. Then, there exists a unique modelled distribution which we call $\mathcal{R}(\diamond, G)$ such that $\mathcal{R}(\diamond, G) = \hat{\mathcal{R}}(\diamond, G)$ on test functions whose support does not intersect $\mathbf{R} \times \partial D$ and such that

$$ (\mathcal{R}(\diamond, G) - \Pi_x(\diamond, G(x))) (\psi_x^{\lambda}) \lesssim \lambda^{-\frac{1}{2}-2\kappa} $$

locally uniformly over $x \in \mathbf{R} \times D$ and uniformly over $\lambda \in (0, 1]$.

Furthermore, there exists $V \in \mathcal{D}^{2-2\kappa, \tilde{w}}$ with $\tilde{w} = (\frac{1}{2}-2\kappa, \frac{1}{2}-2\kappa, 0)$ taking values in the translation invariant sector, and such that $\mathcal{R}V = K\mathcal{R}(\diamond, G)$. The function $V$ is of the form

$$ V = \mathbf{1}_{+} \mathcal{I}(\circ, G) + \Phi, $$

with $\Phi$ taking values in the Taylor polynomials, and the map $(g_1, g_2, g_3, \Pi) \mapsto V$ is uniformly Lipschitz continuous on bounded sets.

**Proof.** We use again the same strategy of proof as in Proposition 6.3, but this time we take as our space $\mathcal{B}$ the space of triples $g = (g_1, g_{2,i}, g_3)$ as in the statement of the proposition and we set

$$ \iota g = \mathbf{1}_{+} \circ G, \quad \hat{\mathcal{R}} g = \mathcal{R}(\diamond, G), $$

where $\mathcal{R}$ is the reconstruction operator given by Theorem C.5. This time, we have $\deg \mathbf{1} = -\frac{3}{2}-2\kappa$ and $\deg \circ = -1-\kappa$, so that it follows from [9, Lem. 4.3] that $\iota g, \diamond, G \in \mathcal{D}^{\tilde{\gamma}, w}$ for $\gamma = \frac{1}{2}-2\kappa$ and $w = (-\frac{3}{2}-2\kappa, -\frac{1}{2}-2\kappa, -\frac{3}{2}-2\kappa)$.
This shows that Theorem C.5 can indeed be applied to this situation since furthermore our admissible models are such that $\hat{\Pi}_{\varepsilon}(\diamond, \tau)(\phi) = 0$ as soon as $\phi$ is supported outside of $\mathbf{R} \times \mathcal{D}$. The remainder of the proof then follows from an application of [9, Lem. 4.12] in the same way as in the proof of Proposition 6.3. $\square$

**Remark 6.5.** Note that the results of [9] do not apply here since $\deg \mathcal{O}$, and $\deg \mathcal{O}$, are both strictly below $-1$. Our saving grace is that the coefficients are sufficiently well-behaved near the boundary of the domain.

*Acknowledgements.* The authors gratefully acknowledge financial support from the Leverhulme Trust through a leadership award, from the European Research Council through a consolidator grant, Project 615897, and from the Royal Society through a research professorship. We are also grateful to Máté Gerencsér for allowing us to use his ideas for the results of “Appendix C” and to Ajay Chandra for discussions on the application of the results of [6].

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

**Publisher's Note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

## Appendix A: Extension of the Kernel

We fix a function $K: \mathbf{R}^{d+1} \to \mathbf{R}$ which is smooth, non-anticipative (that is which is supported on positive times), even in the spatial variable, agrees with the heat kernel on $[0, 1] \times [-1, 1]^d$, and is supported on $[0, 2] \times [-2, 2]^d$. As usual, we write $K = \sum_{n \ge 0} \tilde{K}_n$ with $\tilde{K}_n$ satisfying the conditions of Definition B.1 with $\beta = 2$. We furthermore write

$$K_i(x, y) = (y - x)_i K(y - x).$$

In order to implement integration against the heat kernel with Neumann boundary conditions, we set

$$K_{\text{Neu}}(z, z') = \sum_{n \ge 0} \sum_{R \in \mathcal{G}} \phi(2^n \|z - z'\|) \tilde{K}_n(z - R(z')), \quad (A.1)$$

(recall that the reflection group $\mathcal{G}$ was defined in Section 4.1) where $\psi: \mathbf{R}_+ \to \mathbf{R}_+$ is a smooth function such that $\phi(r) = 1$ for $r \le 2$ and $\phi(r) = 0$ for $r \ge 3$. The kernel $K_{\text{Dir}}$ is defined similarly.

Furthermore, it is obvious that Assumption B.1 is satisfied with $\beta = 2$, since this is the case for $K$ itself. Note that this would not be true if it weren't for the presence of the cutoff $\phi$ in (A.1) since the kernel without cutoff has singularities at $z = R(z')$ for all reflections $R \in \mathcal{G}$.
**Proposition A.1.** One has $K_{\text{Neu}}(z, z') = \sum_{R \in \mathcal{G}} K(z - R(z'))$ for $z, z' \in ([0, 1] \times D)^2$.

*Proof.* Recall that $\tilde{K}_n$ is supported in the ball of radius $2^{-n}$ and note that $\|z - R(z')\| \ge \|z - z'\|$ whenever both $z$ and $z'$ belong to $\mathbf{R} \times D$. As a consequence, for every $z, z' \in ([0, 1] \times D)^2$, every $n \ge 0$ and every $R \in \mathcal{G}$, one has either $\tilde{K}_n(z - R(z')) = 0$ or $\phi(2^n \|z - z'\|) = 1$, so that one can replace $\phi$ by 1 in (A.1). $\square$

In order to state the main result of this appendix, we write $\tilde{D} = \mathbf{R}_+ \times D$ as a shorthand.

**Proposition A.2.** Let $K, K_{\partial}, K_i$ and $K_{\partial,i}$ be as above and let $\zeta, \zeta^c \in C^\alpha$ with $\alpha \le -2$ be such that $\zeta$ is supported on $\tilde{D}$ and $\zeta^c$ is supported on its complement. Then, restricted to $\tilde{D}$, both $K\zeta^c$ and $K_{\partial}\zeta$ belong to $C^{\gamma,w}$ for $w = (\alpha + 2)_3$ and any $\gamma > 0$. Furthermore, when restricted to $\tilde{D}^c$, $K\zeta$ belongs to $C^{\gamma,w}$. The same statements hold when $\alpha \le -3$ with $K$ and $K_{\partial}$ replaced by $K_i$ and $K_{\partial,i}$ respectively and $w = (\alpha + 3)_3$.

*Proof.* Consider the case $z \notin \tilde{D}$, so that

$$(D^k K \zeta)(z) = \sum_{n \ge 0} \langle \zeta, D^k \tilde{K}_n(z - \cdot) \rangle.$$

We now note that $\langle \zeta, D^k \tilde{K}_n(z - \cdot) \rangle = 0$ whenever $2^{-n} \le d(z, \tilde{D})$, while in general

$$|\langle \zeta, D^k \tilde{K}_n(z - \cdot) \rangle| \lesssim 2^{-n(\alpha+2-|k|)}, \quad (A.2)$$

by definition of $C^\alpha$. It immediately follows that one has

$$|(D^k K \zeta)(z)| \lesssim d(z, \tilde{D})^{\alpha+2-|k|}, \quad (A.3)$$

so that one has indeed $K\zeta \in C^{\gamma,(\alpha+2)_3}$ for every $\gamma > 0$. The case of $K\zeta^c$, $K_i\zeta$ and $K_i\zeta^c$ is dealt with in exactly the same way.

Consider now the case $z \in \tilde{D}$ and define

$$\tilde{K}_n^R(z, z') = \phi(2^n|z - z'|) \tilde{K}_n(z - R(z')).$$

We then have the identity

$$(D^k K_{\partial \zeta})(z) = \sum_{n \ge 0} \sum_{R \in \mathcal{G} \setminus \{\text{id}\}} \langle \zeta, D^k \tilde{K}_n^R(z, \cdot) \rangle. \quad (A.4)$$

It follows again from the definition of $C^\alpha$ that the bound (A.2) holds with $\tilde{K}_n(z - \cdot)$ replaced by $\tilde{K}_n^R(z, \cdot)$.

We also note that $\tilde{K}_n^R(z, \cdot) = 0$ unless there exists a point $z'$ such that one has on one hand $\|z - z'\| \lesssim 2^{-n}$ and on the other hand $\|z - R(z')\| \lesssim 2^{-n}$. In particular, there exists a constant $C$ such that, if $d(z, \partial \tilde{D}) \ge C2^{-n}$, one has $\tilde{K}_n^R(z, \cdot) = 0$ unless $R = \text{id}$ (which is excluded from the sum (A.4)). If on the other hand $d(z, \partial \tilde{D}) \le C2^{-n}$, then $\tilde{K}_n^R(z, \cdot)$ is only non-zero for at most eight different reflections $R$. Combining these observations shows as before that the bound (A.3) holds with $K\zeta$ replaced by $K_{\partial}\zeta$. The proof for $K_{\partial}$ replaced by $K_{\partial,i}$ is virtually identical. $\square$

**Corollary A.3.** For $\alpha \in (-3, -2]$ and $\zeta \in C^\alpha$ supported on $\mathbf{R}_+ \times D$, one has $K_{\partial}\zeta \in C^{\alpha+2}$.

*Proof.* This follows immediately from Proposition A.2, combined with the fact that $C^{\gamma,(\eta)_3} \subset C^\eta$ whenever $\gamma \ge \eta > -1$. $\square$
**Proposition A.4.** Let $\zeta \in C^\alpha$ with $\alpha \le -2$ be supported on $\mathbf{R}_- \times \mathbf{R}^d$. Then, restricted to $\mathbf{R}_+ \times \mathbf{R}^d$, $K_{\text{Neu}\zeta}$ belongs to $C^{\gamma,w}$ for $w = (\alpha + 2, \gamma, \alpha + 2)$ and any $\gamma > 0$.

*Proof.* With the same notations as above, we have for $z = (t, x)$,

$$(D^k K_{\text{Neu}\zeta})(z) = \sum_{n \ge 0} \sum_{R \in \mathcal{G}} \langle \zeta, D^k \tilde{K}_n^R(z, \cdot) \rangle.$$

As before, the summand vanishes as soon as $2^{-n} \lesssim \sqrt{|t|}$ and, for all $z \in \mathbf{R}_+ \times \mathbf{R}^d$ and $n \ge 0$, only at most a fixed number of test functions $\tilde{K}_n^R(z, \cdot)$ are non-vanishing, so that we obtain for $t > 0$ the bound

$$|(D^k K_{\text{Neu}\zeta})(t, x)| \lesssim t^{\frac{\alpha+2}{2}-k_0}.$$

This then implies the required bound by the definition of the spaces $C^{\gamma,w}$.

<figure><img src="image_10.png" /></figure>

## Appendix B: Integration and Multiplication by Smooth Functions Almost Commute

We will use a form of the multilevel Schauder theorem of [11, Sec. 5] with slightly weaker assumptions on the kernel $K$. In this section, we fix a space-time scaling $s$ as in [11]. (In our case this would be the parabolic scaling $s = (2, 1, \dots, 1)$.)

**Definition B.1.** A function $K: \mathbf{R}^d \times \mathbf{R}^d \to \mathbf{R}$ is said to be $\beta$-regularising if it can be decomposed as

$$K(x, y) = \sum_{n \ge 0} \tilde{K}_n(x, y), \qquad (B.1)$$

where the functions $\tilde{K}_n$ have the following properties:

* For all $n \ge 0$, the map $\tilde{K}_n$ is supported in the set $\{(x, y) : \|x - y\|_s \le 2^{-n}\}$.

* For any two multiindices $k$ and $l$, there exists a constant $C$ such that the bound

$$|D_x^k D_y^l \tilde{K}_n(x, y)| \le C 2^{|s|-\beta+|l|_s+|k|_s} n, \qquad (B.2)$$

holds uniformly over all $n \ge 0$ and all $x, y \in \mathbf{R}^d$.

**Remark B.2.** This is identical to [11, Ass. 5.1], but with the last condition absent. It turns out that the only place where the third condition of [11, Ass. 5.1] is ever used in [11, Sec. 5] is in the proof of [11, Lemma 5.19], which in turn is only used for the proof of the extension theorem, [11, Thm 5.14].

We assume that we are given a regularity structure of the type studied in [5, 11], endowed with a family of integration maps $\mathcal{I}_k$ for multiindices $k \in \mathbf{N}^d$. We also assume that we are given a model $(\Pi, \Gamma)$ which is admissible for the collection of kernels $K_k$ such that

$$K_k(x, y) = (y - x)^k K_0(x, y). \qquad (B.3)$$

We will furthermore assume that our regularity structure contains the polynomial structure on $\mathbf{R}^d$ (for some fixed scaling) and that admissible models satisfy the usual identity $\Pi_x X^k \tau = (\leftarrow x)^k \Pi_x \tau$ for every $\tau \in \mathcal{T}$. Assuming that $\mathcal{I}_0$ is of order $\beta$, we assume that $\mathcal{I}_k$ is of order $\beta + |k|_s$, which is compatible with (B.3) in the sense that if $K_0$ is $\beta$-regularising in the sense of Definition B.1, then $K_k$ is $(\beta + |k|_s)$-regularising.
**Lemma B.3.** Let $K_0$ be a $\beta$-regularising kernel for some $\beta > 0$ and let $(\Pi, \Gamma)$ be an admissible model for the collection of kernels $K_k$ given in (B.3). Then, one has the identities

$$ \begin{aligned} \Pi_x \mathcal{I}_k(X^\ell \tau) &= \sum_{m \le \ell} \binom{\ell}{m} \Pi_x X^{\ell-m} \mathcal{I}_{k+m}(\tau), \\ \Gamma_{xy} \mathcal{I}_k(X^\ell \tau) - \mathcal{I}_k(\Gamma_{xy}(X^\ell \tau)) &= \sum_{m \le \ell} \binom{\ell}{m} \Gamma_{xy} X^{\ell-m} (\Gamma_{xy} \mathcal{I}_{k+m}(\tau) - \mathcal{I}_{k+m}(\Gamma_{xy} \tau)). \end{aligned} \quad (B.4) $$

*Proof.* We first consider the identity for $\Pi_x$. We can assume that $k=0$ since the general case then follows at once by simply setting $\tilde{\mathcal{I}}_\ell = \mathcal{I}_{k+\ell}$. We also assume without loss of generality that $\tau$ is homogeneous of degree $\alpha$, so that

$$ (\Pi_x \mathcal{I}_0(X^\ell \tau))(y) = \int (z-x)^\ell K(y,z) (\Pi_x \tau)(dz) \\ - \sum_k \frac{1}{k!} \int (z-x)^\ell (y-x)^k D_x^k K(x,z) (\Pi_x \tau)(dz), \quad (B.5) $$

where the sum is constrained by $|k|_S < \alpha + \beta + |\ell|_S$. On the other hand, one has

$$ \begin{aligned} & \sum_m \binom{\ell}{m} (\Pi_x X^m \mathcal{I}_n(\tau))(y) \\ &= \sum_m \binom{\ell}{m} (y-x)^m \int (z-y)^{\ell-m} K(y,z) (\Pi_x \tau)(dz) \\ & \quad - \sum_{m,k} \binom{\ell}{m} (y-x)^m \int \frac{(y-x)^k}{k!} D_x^k ((z-x)^{\ell-m} K(x,z)) (\Pi_x \tau)(dz), \end{aligned} $$

where the sum this time is constrained by $|k|_S < \alpha + \beta + |\ell - m|_S$. The first term in this expression clearly equals the first term of (B.5) by the binomial theorem, so we only need to consider the second term. The integrand can be written as

$$ \begin{aligned} & \sum_{m,k} \binom{\ell}{m} (y-x)^m \frac{(y-x)^k}{k!} D_x^k ((z-x)^{\ell-m} K(x,z)) \\ &= \sum_{m,k,p} \frac{(-1)^p \ell!}{m!(\ell-m-p)!p!(k-p)!} (y-x)^{k+m} (z-x)^{\ell-m-p} D_x^{k-p} K(x,z), \end{aligned} $$

where $p$ is constrained by $p \le k \wedge (\ell - m)$. Setting $q = k - p$ and $r = k + m$, this can be written as

$$ \sum_{q,r,p} \frac{(-1)^p \ell!}{(r-q-p)!(\ell-r+q)!p!q!} (y-x)^r (z-x)^{\ell-r+q} D_x^q K(x,z), $$

where the sum is constrained by the fact that all three variables are positive multiindices and furthermore

$$ |r|_S < \alpha + \beta + |\ell|_S, \quad p+q \le r, \quad r \le \ell+q. $$
It follows that for any fixed values of $q$ and $r$ the above sum vanishes, except when $r-q=0$, so that it equals

$$ \sum_q \frac{1}{q!} (y-x)^q (z-x)^\ell D_x^q K(x, z), $$

constrained by $|q|_S < \alpha + \beta + |\ell|_S$, which is precisely the integrand appearing in the second term of (B.5).

The second identity immediately follows from the first one. Indeed, both sides take values in the Taylor polynomials by the definition of an admissible model. Furthermore, the first identity implies that

$$ \begin{aligned} \Pi_x \mathcal{I}(\Gamma_{xy}(X^\ell \tau)) &= \sum_{p \le \ell} \binom{\ell}{p} (x-y)^{\ell-p} \Pi_x \mathcal{I}(X^p \Gamma_{xy} \tau) \\ &= \sum_{p \le \ell} \sum_{m \le p} \binom{\ell}{p} \binom{p}{m} (x-y)^{\ell-p} \Pi_x X^{p-m} \mathcal{I}_m(\Gamma_{xy} \tau) \\ &= \sum_{m \le \ell} \sum_{n \le \ell-m} \binom{\ell}{m} \binom{\ell-m}{n} (x-y)^{\ell-m-n} \Pi_x X^n \mathcal{I}_m(\Gamma_{xy} \tau) \\ &= \sum_{m \le \ell} \binom{\ell}{m} \Pi_x ((\Gamma_{xy} X^{\ell-m}) \mathcal{I}_m(\Gamma_{xy} \tau)). \end{aligned} \quad (B.6) $$

As a consequence of combining this with the first identity of (B.4), both sides of the second identity of (B.4) are equal after applying $\Pi_x$ to them. Since furthermore both sides belong to the space of Taylor polynomials on which $\Pi_x$ is injective, the claim follows. $\square$

**Corollary B.4.** In the context of Lemma B.3, write $K_k$ for the integration operator associated to $\mathcal{I}_k$ as in [11] and write $\mathcal{L}_\gamma$ for the Taylor lift $\mathcal{C}^\gamma \to \mathcal{D}^\gamma$. Then, for every $F \in \mathcal{D}_\chi^\gamma$ with $\chi \le 0 < \gamma$ and $g \in \mathcal{C}^{\theta-\chi}$ with $\theta \in (0, \gamma]$ (and such that $\theta + \beta \notin \mathbb{Z}$), one can find a function $\phi \in \mathcal{C}^{\theta+\beta}$ such that, setting

$$ G = \mathcal{L}_{\theta+\beta}\phi + \sum_{|\ell|_S < \theta-\chi} \frac{1}{\ell!} \mathcal{L}_{\theta-\chi-|\ell|_S} (D^\ell g) K_\ell(F), \quad (B.7) $$

one has $G \in \mathcal{D}^{\theta+\beta}$ and $RG = K(gRF)$.

*Proof.* A straightforward calculation virtually identical to (B.6) shows that, as a consequence of the first identity of Lemma B.3, one has

$$ \Pi_x \mathcal{I}(\mathcal{L}_{\theta-\chi}(g)(x)F(x)) = \sum_{|\ell|_S < \theta-\chi} \frac{1}{\ell!} D^\ell g(x) \Pi_x \mathcal{I}_\ell(F(x)). \quad (B.8) $$

For $F \in \mathcal{D}_\chi^\gamma$, write $N_k(F) = K_k(F) - \mathcal{I}_k(F)$, so that $N_k(F)$ takes values in the Taylor polynomials of degree at most $\gamma + \beta$, and set

$$ \Phi = N_0(\mathcal{L}_{\theta-\chi}(g)F) - \sum_{|\ell|_S < \theta-\chi} \frac{1}{\ell!} \mathcal{L}_{\theta-\chi-|\ell|_S} (D^\ell g) N_\ell(F). $$

We then set

$$ G = \Phi + \sum_{|\ell|_S < \theta-\chi} \frac{1}{\ell!} \mathcal{L}_{\theta-\chi-|\ell|_S} (D^\ell g) K_\ell(F) $$
$$= \mathcal{N}_0(F) + \sum_{|\ell|_s < \theta - \chi} \frac{1}{\ell!} \mathcal{L}_{\theta-\chi-|\ell|_s} (D^\ell g) \mathcal{I}_\ell(F).$$

Combining this with (B.8), it follows that, setting $\tilde{G} = \mathcal{K}_0(\mathcal{L}_{\theta-\chi}(g)F)$, one has

$$\Pi_x G(x) = \Pi_x \mathcal{N}_0(\mathcal{L}_{\theta-\chi}(g)F)(x) + \Pi_x \mathcal{I}(\mathcal{L}_{\theta-\chi}(g)(x)F(x)) = \Pi_x \tilde{G}(x). \quad (B.9)$$

Combining this with the second identity of Lemma B.3 and writing $\tilde{Q}$ for the projection onto the Taylor polynomials, we conclude that

$$\tilde{Q}(G(x) - \Gamma_{xy}G(y)) = \tilde{Q}(\tilde{G}(x) - \Gamma_{xy}\tilde{G}(y)).$$

Since furthermore $\tilde{G}$ and each of the terms $\mathcal{L}_{\theta-\chi-|\ell|_s}(D^\ell g) \mathcal{K}_\ell(F)$ belong to $\mathcal{D}^{\theta+\beta}$, we conclude that one necessarily has $G \in \mathcal{D}^{\theta+\beta}$. This in turn implies that $\Phi \in \mathcal{D}^{\theta+\beta}$, so that it is the lift of a function $\phi \in \mathcal{C}^{\theta+\beta}$. By (B.9), we furthermore have $\mathcal{R}G = \mathcal{R}\tilde{G} = K(g\mathcal{R}F)$, thus concluding the proof. $\square$

One simple but very useful corollary of this result can be formulated as follows:

**Corollary B.5.** Let $\zeta \in \mathcal{C}^\chi$ with $\chi \le 0$, let $K$ and $K_k$ be as above, and let $g \in \mathcal{C}^{\theta-\chi}$ with $\theta > 0$ and $\theta + \beta \notin \mathcal{Z}$. Then,

$$K(g\zeta) - \sum_{|\ell|_s < \theta - \chi} \frac{D^\ell g}{\ell!} K_\ell \zeta \in \mathcal{C}^{\theta+\beta}.$$

*Proof.* It suffices to consider the case of a regularity structure with symbol $\Xi$ (plus Taylor polynomials and their products with $\Xi$) and model mapping $\Xi$ to $\zeta$. We then apply the reconstruction operator to both sides of (B.7) with $F = \Xi$. $\square$

In our context, we will need an analogous result, but for the spaces $\mathcal{D}^{\gamma, \eta}$ of [11, Sec. 6]. Furthermore, we will need to be able to cover situations in which [11, Prop. 6.16] does not apply because we consider elements taking values in a sector of regularity below $-2$, so that the reconstruction theorem [11, Prop. 6.9] fails. We therefore make use instead of [9, Lem. 4.12] which, given $F \in \mathcal{D}^{\gamma, \eta}(V)$ with $V$ of regularity $\alpha \ge \eta$, allows us to specify a distribution $\zeta \in \mathcal{C}^\eta$ such that $(\mathcal{R}F)(\phi) = \zeta(\phi)$ for every test function $\phi$ whose support does not intersect the plane $P = \{(t, x) : t = 0\}$. We then have the following result.

**Corollary B.6.** Let $\gamma > 0$, let $V$ be a sector of regularity $\alpha \le 0$, and let $w = (\eta, \sigma, \mu)$ with $\eta, \sigma, \mu \le \alpha$ and $\eta \le \sigma \wedge \mu$, $\eta + \beta > -2$, $\sigma + \beta > -1$. For each admissible model $\Pi$, let $B$ be a Banach space equipped with a bounded map $\iota: B \to \mathcal{D}^{\gamma, w}(V)$. (Since the latter depends on $\Pi$ in general, this can also be the case for $B$ and/or $\iota$. In this case, we assume that $\iota$ is bounded independently of the underlying model.) Let furthermore $\hat{\mathcal{R}}: B \to \mathcal{C}^\eta$ be a continuous linear operator such that $(\hat{\mathcal{R}}F)(\phi) = (\mathcal{R}\iota F)(\phi)$ for every $F \in B$ and every test function $\phi$ whose support does not intersect the two boundaries. We furthermore assume that we have a continuous bilinear map

$$C^{\gamma-\eta} \times B \to B \quad (g, F) \mapsto gF$$

such that $\iota(gF) = \mathcal{L}_{\gamma-\eta}(g)\iota F$ and such that

$$\hat{\mathcal{R}}(gF) = g\hat{\mathcal{R}}(F),$$

where the right hand side is meaningful thanks to the fact that $\gamma > 0$. Then, with $\tilde{\gamma}$ and $\tilde{w}$ as in [9, Lem. 4.12], one can find $\phi \in \mathcal{C}^{\tilde{\gamma}, \tilde{w}}$ such that the modelled distribution $G$ given by (B.7) (with suitably defined $K_\ell$, see the proof) belongs to $\mathcal{D}^{\tilde{\gamma}, \tilde{w}}$ and satisfies $\mathcal{R}G = K_0(g\hat{\mathcal{R}}F)$.