1. Introduction
We consider the stochastic nonlinear beam equation (SNLB) on 
$\mathbb{T}^4 = (\mathbb{R}/\mathbb{Z})^4$ with additive space-time white noise:
\begin{equation}
\begin{cases}
\partial_t^2 u + \Delta^2 u \pm u^k = \xi , \\
(u, \partial_t u)\vert_{t=0} = (u_0, u_1) \in \mathcal{H}^s(\mathbb{T}^4),
\end{cases}
\end{equation} where 
$u: \mathbb{R}_+\times\mathbb{T}^4 \to \mathbb{R}$, Δ2 denotes the bi-harmonic operator, 
$k\ge2$ is a natural number, ξ is a space-time white-noise on 
$\mathbb{R}_+\times \mathbb{T}^4$, and 
$\mathcal{H}^s(\mathbb{T}^4) = H^s(\mathbb{T}^4) \times H^{s-2}(\mathbb{T}^4)$. We refer to the Eq. (1.1) with ‘+’ as defocusing and with ‘−’ as focusing.
The deterministic beam equation appears in the literature under various names, such as the fourth-order wave equation, the extensible beam/plate equation, and the Bretherton equation. In the one-dimensional setting, it was first derived by Bretherton in [Reference Bretherton8] to describe the weak interaction between dispersive waves and it has a variety of applications in physics and mechanics; see [Reference Peletier and Troy50] and references therein. We also refer to the non-local model derived by Woinowsky-Krieger [Reference Woinowsky-Krieger57] to describe the vibration of a clamped extensible beam.
Our main goal is to establish low regularity well-posedness of (1.1) on 
$\mathbb{T}^4$ with space-time white-noise, which is of analytical interest due to the roughness of the noise. This study on 
$\mathbb{T}^d$ for 
$d=1,2,3$ was pursued in [Reference Mosincat, Pocovnicu, Tolomeo and Wang38, Reference Tolomeo53, Reference Tolomeo54]. We also mention the results in [Reference Brzeźniak, Maslowski and Seidler9, Reference Brzeźniak, Ondreját and Seidler10, Reference Chow13] for a non-local version of (1.1) with multiplicative noises. For the study of the deterministic nonlinear beam equation, we refer the interested readers to [Reference Lebeau31, Reference Pausader48, Reference Pausader49] and references therein.
The main difficulty in studying (1.1) on 
$\mathbb{T}^4$ comes from the roughness of the noise ξ. To illustrate this, we first consider the mild formulation of (1.1):
\begin{equation}
u(t) = S(t)(u_0,u_1) \mp \int_0^t \frac{\sin((t-t')\Delta)}{\Delta} u^k(t')\,dt' + \Psi(t),
\end{equation}where S(t) denotes the linear propagator
\begin{equation}
S(t)(u_0,u_1) = \cos(t\Delta) u_0 + \frac{\sin(t\Delta)}{\Delta} u_1,
\end{equation} with the understanding that 
$\frac{\sin(t \cdot 0)}{0} = t $, and Ψ is the stochastic convolution which solves the linear stochastic beam equation on 
$\mathbb{T}^4$:
\begin{equation}
\partial_t^2 \Psi + \Delta^2 \Psi = \xi.
\end{equation}More precisely, Ψ is given by
\begin{align}
\Psi (t) = \int_0^t \frac{\sin ((t - t') \Delta)}{\Delta} dW (t'),
\end{align} where W denotes a cylindrical Wiener process on 
$L^2 (\mathbb{T}^4)$:
\begin{align}
W(t,x) \stackrel{\textrm{def}}{=} \sum_{n\in \mathbb{Z}^{4} } \beta_n(t) e_n(x),
\end{align} withFootnote 1 
$
e_n(x) = e^{2\pi in\cdot x},
$ and 
$\{\beta_n\}_{n \in \mathbb{Z}^4 }$ a family of mutually independent complex-valued Brownian motions conditioned to 
$\beta_{-n} = \overline{\beta_n}$, 
$n\in\mathbb{Z}^4$, with variance 
$\operatorname{Var}(\beta_n(t))=t$. One can show that W lies almost surely inFootnote 2 
$C^{\alpha}(\mathbb{R}_+;H^{-2-\varepsilon}(\mathbb{T}^4))$ for any 
$\alpha \lt \frac{1}{2}$ and ɛ > 0. Therefore, due to the two degrees of spatial smoothing of the linear beam equation, it follows that 
$\Psi(t) \in H^{-\varepsilon}(\mathbb{T}^4) \setminus L^2(\mathbb{T}^4)$ almost surely, for any ɛ > 0, thus it is merely a distribution. Consequently, we expect the solution u to (1.2) to also only be a distribution and thus the product uk is classically ill-defined. To overcome this difficulty, we closely follow the work of Gubinelli-Koch-Oh [Reference Gubinelli, Koch and Oh22] for wave equations (see also [Reference Oh and Thomann47]), and construct solutions 
$u=\Psi+v$ which solve a suitably renormalized version of (1.1).
We now detail this renormalization procedure. We first smooth the noise ξ in (1.1) via Fourier truncation and consider the truncated stochastic convolution 
$\Psi_N$ given by
\begin{align*}
\Psi_N(t, x) = \pi_N \Psi(t, x)
= \sum_{\substack{n \in \mathbb{Z}^4 \\ |n| \leq N}} e_n(x)
\int_0^t \frac{\sin ((t - t') |n|^2)}{|n|^2 } d \beta_n (t'),
\end{align*} where πN denotes the frequency truncation onto 
$\{|n| \le N\}$. Then, for each fixed 
$x \in \mathbb{T}^4$ and 
$t \geq 0$, it follows from the Ito isometry that the random variable 
$\Psi_N(t, x)$ is a mean-zero real-valued Gaussian random variable with variance
\begin{align}
\sigma_N(t) &= \mathbb{E} \big[\Psi_N^2(t, x)\big]
\sim t \log N,
\end{align} which is independent of 
$x \in \mathbb{T}^4$.
Let uN be the solution to SNLB (1.1) with the regularized noise 
$\pi_N \xi$, which satisfies the mild formulation (1.2) with the truncated stochastic convolution 
$\Psi_N$. Motivated by (1.2), we introduce the first order expansion [Reference Bourgain6, Reference Da Prato and Debussche17, Reference McKean and McKean36]:
\begin{equation}
u_N = \Psi_N + v_N,
\end{equation}where the remainder vN solves the following nonlinear beam equation:
\begin{equation}
\partial_t^2 v_N + \Delta^2 v_N \pm \sum_{\ell=0}^k {k\choose \ell} \Psi_N^\ell v_N^{k-\ell} = 0.
\end{equation} Unfortunately, due to (1.7), the monomials 
$\Psi_N^\ell$ in (1.9) do not have good limiting behavior as 
$N\to\infty$. Instead, we define the Wick-ordered power 
$\mathcal{W}_\sigma ( \Psi^\ell_N ) $ as
\begin{align}
\mathcal{W}_\sigma ( \Psi^\ell_N(t, x) ) \, \stackrel{\text{def}}{=} H_\ell(\Psi_N(t, x);\sigma_N(t)),
\end{align} where 
$H_\ell(x; \sigma )$ is the Hermite polynomial of degree 
$\ell$, which can be shown to converge to a limit 
$\mathcal{W}_\sigma (\Psi^\ell) $ in 
$L^p(\Omega; C([0,T];W^{-\varepsilon,\infty}(\mathbb{T}^4)))$, for any 
$1\le p \lt \infty$ and ɛ > 0 as 
$N \to \infty$; see 
$\S$ 2.2. We then consider the Wick renormalized version of (1.9)
\begin{align}
\partial_t^2 v_N + \Delta^2 v_N \pm \sum_{\ell=0}^k {k\choose \ell} \mathcal{W}_\sigma (\Psi_N^\ell) v_N^{k-\ell} = 0,
\end{align} which converges, as 
$N\to\infty$, to the following equation:
\begin{align}
\partial_t^2 v + \Delta^2 v \pm \sum_{\ell=0}^k {k\choose \ell} \mathcal{W}_\sigma (\Psi^\ell) v^{k - \ell} = 0.
\end{align} Lastly, from (1.8) and (2.3) below, we can define the Wick-ordered nonlinearity 
$\mathcal{W}_\sigma (u_N^k) $ as
\begin{equation*}
\mathcal{W}_\sigma (u_N^k (t, x)) \stackrel{\textrm{def}}{=} H_k(u_N(t, x);\sigma_N(t)) = \sum_{\ell=0}^k {k\choose \ell} \mathcal{W}_\sigma ( \Psi_N^\ell (t,x) )
v_N^{k-\ell} (t, x).
\end{equation*} Consequently, if vN solves (1.11), then 
$u_N = \Psi_N + v_N$ satisfies the following truncated Wick renormalized SNLB:
\begin{align}
\partial_t^2 u_N + \Delta^2 u_N \pm \mathcal{W}_\sigma (u_N^k) = \pi_N \xi.
\end{align} Similarly, with 
$u=\Psi + v$ for some suitable v, we define the Wick-ordered nonlinearity as
\begin{align}
\mathcal{W}_\sigma ( u^k) \stackrel{\textrm{def}}{=} \sum_{\ell=0}^k {k\choose \ell} \mathcal{W}_\sigma ( \Psi^\ell ) v^{k - \ell},
\end{align} and so if v solves (1.12), then 
$u = \Psi + v$ solves the following Wick renormalized SNLB:
\begin{align}
\partial_t^2 u + \Delta^2 u \pm \mathcal{W}_\sigma (u^k) = \xi.
\end{align}Before stating our first main result on local well-posedness of (1.15), let us discuss the scaling critical regularity associated to the deterministic nonlinear beam equation (NLB):
\begin{align}
\partial_t^2 u + \Delta^2 u \pm u^k = 0.
\end{align} On 
$\mathbb{R}^4$, (1.16) enjoys the following scaling symmetry: if u is a solution to (1.16) then 
$u_\lambda(t,x) \stackrel{\textrm{def}}{=} \lambda^{\frac{4}{k-1}} u(\lambda^2 t, \lambda x)$ is also a solution to (1.16). This induces the scaling critical Sobolev index 
$s_\text{scaling} = 2 - \frac{4}{k-1}$, i.e., the homogeneous Sobolev 
$\dot{H}^{s}(\mathbb{R}^4)$-norm with 
$s=s_{\text{scaling}}$ is invariant under the scaling. Moreover, for a given integer 
$k \ge 2$, we define s crit by
\begin{align}
s_\text{crit} \stackrel{\textrm{def}}{=} \max(s_\text{scaling}, 0) = \max\bigg(2 - \frac{4}{k-1}, 0\bigg),
\end{align} where the restriction 
$s_{\text{crit}} \ge0$ appears in making sense of the powers of u. Although the scaling symmetry does not extend to 
$\mathbb{T}^4$, the numerology still plays an important role in predicting local well-posedness issues. In particular, our aim is to show that the SNLB (1.1) is locally well-posed in the scaling (sub)critical Sobolev spaces 
$\mathcal{H}^s(\mathbb{T}^4)$ with 
$s \geq s_{\text{crit}}$. In fact, we show pathwise local well-posedness of (1.15) in the subcritical regime for 
$s \gt s_{\text{crit}}$ and all order nonlinearities 
$k\ge2$, and also in the critical case (
$s=s_{\text{crit}}$) for 
$k\ge4$.
Theorem 1.1 Given an integer 
$k \geq 2$, let 
$s_\textrm{crit}$ be as in (1.17). Then, the Wick renormalized SNLB (1.15) is pathwise locally well-posed in 
$\mathcal{H}^s(\mathbb{T}^4)$ for
\begin{equation*} \textrm{(i) } k \geq 4:
\ s \geq s_\textrm{crit}
\qquad \text{or} \qquad
\textrm{(ii) } k = 2, 3:
\ s \gt s_\textrm{crit}.
\end{equation*} More precisely, given any 
$(u_0, u_1) \in \mathcal{H}^s(\mathbb{T}^4)$, there exists an almost surely positive stopping time 
$T = T( \omega, u_0, u_1) \gt 0$ such that there exists a unique solution u to (1.15) on 
$[0, T]$ with 
$(u, \partial_t u)\vert_{t=0} = (u_0,u_1)$ and
\begin{align*}
u \in \Psi + C([0, T]; H^{s'}(\mathbb{T}^4))
\subset C([0, T]; H^{-\varepsilon}(\mathbb{T}^4))
\end{align*} for any ɛ > 0, where 
$s' = \min(s, 2-\varepsilon)$.
The solution u in Theorem 1.1 is understood as 
$u=\Psi +v$ where we construct 
$(v,\partial_t v) \in C([0,T]; \mathcal{H}^{s'}(\mathbb{T}^4))$ with v solving the following Duhamel formulation:
\begin{align}
v(t)=
S(t) (u_0, u_1)
\mp \int_{0}^t \frac{\sin ((t-t') \Delta )}{\Delta} \mathcal{W}_\sigma (u^k(t')) dt',
\end{align} for 
$\mathcal{W}_\sigma (u^k) $ and S(t) as in (1.14) and (1.3), respectively. The main ingredient in proving Theorem 1.1 in the (almost) critical regime are the Strichartz estimates for the beam equation. In the Euclidean setting, by exploiting the formal decomposition
\begin{equation*}\partial_t^2 +\Delta^2 = (i\partial_t + \Delta)(-i\partial_t +\Delta), \end{equation*} which sheds light on the relation between the beam equation and the Schrödinger equation, and the analysis of oscillatory integrals, Pausader [Reference Pausader48, Reference Pausader49] established Strichartz estimates for the beam equation. However, in contrast to the wave equation, the lack of finite speed of propagation poses difficulties in transferring these estimates from the Euclidean to the periodic setting. Instead, we exploit the connection between the operator S(t) in (1.3) appearing in (1.18) and the free Schrödinger operators 
$e^{\pm i t\Delta}$ via the periodic Schrödinger Strichartz estimates in [Reference Bourgain and Demeter7, Reference Killip and Vişan27] from the 
$\ell^2$-decoupling theory. See 
$\S$ 3 for details.
(i) In Theorem 1.1, we cannot reach the critical regularity
$s=s_{\text{crit}}=0$ for the quadratic and cubic SNLB (1.15), $k=2,3$. This restriction comes from the sharp Strichartz estimates for Schrödinger (see Lemma 3.1), where the endpoint p = 3 is not included, which is needed for our argument in the critical setting for $k=2,3$. Strichartz estimates for p = 3 are known to only hold with a derivative loss [Reference Bourgain4, Reference Bourgain and Demeter7], which prevents us from taking s = 0. Thus our result is sharp with respect to the method. It may be possible to reach the critical regularity in these cases by using the Up-Vp spaces introduced in [Reference Koch and Tataru28].(ii) The proof of Theorem 1.1 can be easily adapted to show local well-posedness of the truncated Wick-ordered SNLB (1.13), uniformly in N. In fact, it follows that for
$(u_0,u_1)\in \mathcal{H}^s(\mathbb{T}^4)$, there exists an almost surely positive stopping time $T=T(\omega, u_0,u_1) \gt 0$ independent of N and a unique solution $u_N \in \Psi_N + C([0,T]; H^{s'}(\mathbb{T}^4))$ to (1.13). Moreover, we can show that uN converges to the corresponding solution u to (1.15). We note that although this seems to depend on regularizing by πN, one can consider a different regularization procedure, such as mollification. Indeed, one can show that the Wick-ordered monomials $\mathcal{W}_\sigma (\Psi^k) $ are independent of the choice of mollifier, and thus so is the renormalized nonlinearity (1.14). See [Reference Gubinelli, Koch and Oh22, Remark 1.2] for further discussion.
Our next goal is to extend the solutions constructed in Theorem 1.1 globally-in-time. We restrict our attention to the defocusing case (‘+’ sign in (1.1)) and odd-ordered nonlinearities, as the energies corresponding to the deterministic NLB equation are sign definite in this setting. First, we construct pathwise global-in-time solutions for the cubic defocusing Wick-ordered (1.15) by adapting the hybrid method of Gubinelli-Koch-Oh-Tolomeo [Reference Gubinelli, Koch, Oh and Tolomeo24] to the beam equation. Then, we use Bourgain’s invariant measure argument to show almost sure global well-posedness and invariance of the Gibbs measure for the defocusing damped Wick renormalized SNLB with odd-power nonlinearities.
We first consider the cubic Wick renormalized SNLB (1.15) in the defocusing case, with k = 3 and ‘+’ sign. In Theorem 1.1, for s > 0 we constructed a solution 
$u=\Psi +v$ where the remainder v solves
\begin{equation}
\partial_t^2 v + \Delta^2 v + \mathcal{W}_\sigma (u^3) =0 ,
\end{equation} and 
$\mathcal{W}_\sigma (u^3) $ is given in (1.14). A consequence of the (deterministic) contraction argument used to show Theorem 1.1 is the following (almost sure) blow-up alternative: either the solution v exists globally in time or there exists some finite time 
$T_* = T_*(\omega) \gt 0$ such that
\begin{align}
\lim_{t \nearrow T_*}\| \vec v(t) \|_{\mathcal{H}^{s'}} = \infty,
\end{align} where 
$\vec v = (v, \partial_t v)$ and 
$s' = \min(s, 2- \varepsilon)$ for any small ɛ > 0.
To globalize solutions, we must control the growth of the norm in (1.20). In the parabolic setting, there are various results where deterministic arguments have been adapted to the stochastic setting to directly control the growth of norms of solutions; see [Reference Gubinelli and Hofmanová21, Reference Moinat and Weber37, Reference Mourrat and Weber39, Reference Mourrat and Weber40]. Unfortunately, for (1.15), due to the lack of a strong smoothing effect, such arguments do not apply. Instead, even in the deterministic setting, we must consider conservation laws. For the deterministic cubic nonlinear beam equation (NLB):
\begin{equation*}
\partial_t^2 v + \Delta^2 v +
v^3 = 0,
\end{equation*}the associated energy
\begin{align}
E(\vec v ) = \frac{1}{2}\int_{\mathbb{T}^4} (\Delta v)^2 dx
+
\frac{1}{2}\int_{\mathbb{T}^4} (\partial_t v)^2dx
+ \frac1{4} \int_{\mathbb{T}^4} v^{4} dx,
\end{align} gives control over the 
$\mathcal{H}^2(\mathbb{T}^4)$-norm of 
$\vec v$, as this quantity is conserved for sufficiently regular solutions. Unfortunately, when adding noise to the equation and considering a solution v to (1.19), two problems arise: (i) the energy 
$E(\vec v)$ is not conserved under the dynamics of (1.19), and (ii) since 
$\vec v\in \mathcal{H}^{s'}(\mathbb{T}^4) \setminus \mathcal{H}^{2}(\mathbb{T}^4)$ for 
$s'=\min(s,2-\varepsilon)$ for any ɛ > 0, the energy 
$E(\vec{v})$ is actually infinite.
In the context of the two-dimensional cubic stochastic nonlinear wave equation, Gubinelli-Koch-Oh-Tolomeo [Reference Gubinelli, Koch, Oh and Tolomeo24] introduced a new hybrid method to overcome these difficulties, by combining the I-method of Colliander-Keel-Staffilani-Takaoka-Tao [Reference Colliander, Keel, Staffilani, Takaoka and Tao15, Reference Colliander, Keel, Staffilani, Takaoka and Tao16] and the Gronwall-type globalization argument by Burq-Tzvetkov [Reference Burq and Tzvetkov12]. See also [Reference Forlano19, Reference Tolomeo55] for other instances of this method. To establish our next main result, we adapt this argument to show pathwise global well-posedness of (1.19).
Theorem 1.3 Let 
$s \gt \frac 74$. Then, the defocusing cubic Wick renormalized SNLB (1.19) is globally well-posed in 
$\mathcal{H}^s(\mathbb{T}^4)$. More precisely, given any 
$(u_0, u_1) \in \mathcal{H}^s(\mathbb{T}^4)$, the solution v to the equation (1.19) exists globally in time and 
$(v, \partial_t v) \in C(\mathbb{R}_+; \mathcal{H}^{s'}(\mathbb{T}^4))$, almost surely, for 
$s' = \min(s, 2- \varepsilon)$ for any small ɛ > 0.
We briefly detail the ideas of the proof of Theorem 1.3. For simplicity, let 
$\frac74 \lt s \lt 2$ so that 
$s'=s$. In view of the blow-up alternative (1.20), our main goal is to control the 
$H^s(\mathbb{T}^4)$-norm of the solution v to (1.19), where the conservation of 
$E(\vec v)$ is not useful. Instead, the I-method is based on studying the growth of a modified energy obtained from 
$E(\vec v)$ which controls the 
$H^s(\mathbb{T}^4)$-norm of v. In particular, for 
$N\in\mathbb{N}$, we consider 
$E (I \vec{v})$ where 
$I=I_N$ denotes the I-operator, a Fourier operator with a smooth, radially symmetric, non-increasing multiplier mN given by
\begin{equation}
m_N(\xi)=
\begin{cases}
1, & \text{if } |\xi|\le N, \\
\Big(\frac{N}{|\xi|}\Big)^{2-s}, & \text{if } |\xi|\ge2N.
\end{cases}
\end{equation} Note that 
$If \in H^2(\mathbb{T}^4)$ if and only if 
$f \in H^s(\mathbb{T}^4)$; see (4.1)-(4.2).
To study the growth of 
$E(I \vec v)$, we consider the following I-SNLB:
\begin{equation*}
\partial_t^2 I v + \Delta^2 I v +
I \mathcal{W}_\sigma (u^3) = 0,
\end{equation*} where 
$\mathcal{W}_\sigma (u^3) $ is as in (1.14). Unfortunately, the modified energy 
$E(I\vec v)$ is not conserved under the flow of I-SNLB, and by direct computation we obtain
\begin{equation}
\begin{aligned}
E(I\vec v)(t_2) - E(I\vec v)(t_1)
& = \int_{t_1}^{t_2} \int_{\mathbb{T}^4} (\partial_t I v) \big\{- I (v^3) + (Iv)^3\big\} dx dt' \\
&
-3
\int_{t_1}^{t_2}\int_{\mathbb{T}^4} (\partial_t I v) \big\{I(v^2 \Psi) + I(v \mathcal{W}_\sigma ( \Psi^2 ) ) \big\} dx dt' \\
&
-
\int_{t_1}^{t_2} \int_{\mathbb{T}^4} (\partial_t I v) I(\mathcal{W}_\sigma ( \Psi^3 ) ) dx dt' ,
\end{aligned}
\end{equation} for 
$0\le t_1 \lt t_2$. The first term, due to the I-operator, requires a certain (deterministic) commutator estimate; see Lemma 4.1. The difficulty in the remaining contributions comes from the roughness of Ψ, which is handled by exploiting a finer regularity property of 
$I\Psi$ combined with commutator estimates and a Gronwall-type argument. Finally, due to the growth of the modified energy 
$E(I_N \vec v)$, we iterate the argument above over time-intervals of fixed size, but with an increasing sequence Nk of parameters for the I-operator, extending the solution to (1.19) globally-in-time. See 
$\S$ 4 for details.
(i) There is a gap between the global well-posedness result for the Wick-ordered cubic SNLB (1.19) in Theorem 1.3 for
$s \gt \frac74$ and the local well-posedness threshold s > 0 from Theorem 1.1. The technical assumption of $s \gt \frac74$ comes from controlling the growth of the energy (Proposition 4.5) and that of the chosen sequence of parameters Nk in a way that allows for an iterative argument (see (4.38) and (4.40)). We do not believe this restriction to be sharp, and it may be possible to improve it by refining the I-method part of the argument. However, we do not pursue this issue in this paper.(ii) At this point, we do not know how to extend pathwise global well-posedness of the Wick ordered defocusing SNLB (1.15) to a (super-)quintic nonlinearity. As mentioned earlier, the method of proof for Theorem 1.3 is partially based on the Gronwall-type globalization argument by Burq-Tzvetkov [Reference Burq and Tzvetkov12], which only applies to the cubic case. Indeed, the main restriction comes from the term
\begin{equation*} \int_{\mathbb{T}^4} (\partial_t v) (v^2 \Psi) dx, \end{equation*}appearing on the second contribution on the right-hand side of (1.23), where we dropped the I-operator for simplicity. In order to estimate this contribution by a power of the energy
$E(\vec{v})$ in (1.21), by Cauchy-Schwarz inequality, we must place $\partial_t v$ in $L^2(\mathbb{T}^4)$, which implies that $v^2 \Psi$ is also in $L^2(\mathbb{T}^4)$. Consequently, we obtain the L 4-norm of v, which is also controlled by the energy $E(\vec{v})$. However, one can see that the same argument fails for the analogous term for higher order nonlinearities; see, for example, the case-by-case analysis in [Reference Liu33, Section 5]. To deal with a (super-)quintic nonlinearity, one needs to exploit some other ideas such as those in [Reference Latocca30, Reference Oh and Pocovnicu43], but we choose not to pursue this issue in this paper.(iii) A standard application of the I-method results in a polynomial growth bound (in time) on the Sobolev norm of a solution. See, for example, [Reference Colliander, Keel, Staffilani, Takaoka and Tao16, Section 6]. The hybrid argument used for Theorem 1.3 yields a double exponential growth bound on the
$\mathcal{H}^s$-norm of the solution; see Remark 4.6 below. It may be possible to improve this double exponential bound, but we expect that one can obtain at best a polynomial growth bound for SNLB (1.19) due to the polynomial growth (in time) of the stochastic convolution Ψ. One can compare this situation with the damped case in the next subsection, where the invariant measure argument yields a logarithmic growth bound; see remark 1.7(i) below.
Lastly, we restrict our attention to the following (defocusing) stochastic damped nonlinear beam equation (SdNLB):
\begin{align}
\partial_t^2 u + \partial_t u + (1 - \Delta)^2 u + u^k = \sqrt 2\xi,
\end{align} for 
$k \in 2 \mathbb{N} + 1$. By modifying the proof of Theorem 1.3, we can show global well-posedness for the damped dynamics (1.24) when k = 3, after renormalization, but we do not know how to extend this deterministic argument to higher nonlinearities. Instead, we consider a probabilistic approach and establish almost sure global well-posedness of (1.24) and invariance of the Gibbs measure 
$\vec \rho$ via Bourgain’s invariant measure argument [Reference Bourgain5, Reference Bourgain6], where 
$\vec \rho$ is formally given by
\begin{align}
\text{"}d\vec{\rho} (u, \partial_t u ) = Z^{-1}e^{-E(u,\partial_t u )}du\,d(\partial_t u)\text{"}.
\end{align} Here, 
$E(u, \partial_t u)$ denotes the energy (or Hamiltonian) of the deterministic undamped defocusing NLB (1.16):
\begin{align}
E(u, \partial_t u ) = \frac{1}{2}\int_{\mathbb{T}^4} [( 1- \Delta) u]^2 dx
+
\frac{1}{2}\int_{\mathbb{T}^4} (\partial_t u)^2dx
+ \frac1{k+1} \int_{\mathbb{T}^4} u^{k+1} dx.
\end{align} We can understand SdNLB (1.24) as a superposition of the defocusing NLB dynamics (1.16) and the Ornstein-Uhlenbeck dynamics (for the component 
$\partial_t u$):
\begin{align*}
\partial_t (\partial_t u) = - \partial_t u +
\sqrt 2 dW.
\end{align*} The latter leaves the Gibbs measure 
$\vec{\rho}$ invariant, which is also expected to hold under the dynamics of NLB (1.16) due to its Hamiltonian structure; see [Reference Oh, Seong and Tolomeo45] and [Reference Xia58, Chapter 3]. Therefore, we expect 
$\vec{\rho}$ to also be invariant under SdNLB (1.24).
Moreover, this invariance is also inferred from the stochastic quantization viewpoint. In fact, (1.24) is the so-called canonical stochastic quantization equation of the 
$\Phi^{k + 1}_4$-model; see [Reference Ryang, Saito and Shigemoto51]. We thus refer to (1.24) as the hyperbolic 
$\Phi_4^{k + 1}$-model, which is of importance in constructive quantum field theory. The invariance of the Gibbs measure is also related to other applications in physics such as the study of equilibrium states, couplings of fields, and scattering of particles; see [Reference Aizenman1–Reference Aizenman and Graham3, Reference Feldman, Magnen, Rivasseau and Sénéor18, Reference Gawȩdzki and Kupiainen20, Reference Hara25, Reference Sokal52] and references therein. See also [Reference Gubinelli, Koch and Oh22, Reference Gubinelli, Koch, Oh and Tolomeo24, Reference Liu, Tzvetkov and Wang35, Reference Oh, Robert and Tzvetkov44, Reference Tolomeo54] for further results on wave-like 
$\Phi_d^{k+1}$-models.
Our first step is to rigorously construct the measure 
$\vec{\rho}$, since the expression in (1.25) is only formal. We want to define 
$\vec{\rho}$ as a weighted Gaussian measure of the form
\begin{align}
\text{"}d\vec{\rho} (u, \partial_t u) = Z^{-1} e^{
- \frac1{k+1} \int_{\mathbb{T}^4} u^{k+1} dx}
d\vec{\mu}_2 (u, \partial_t u)\text{''},
\end{align} where 
$\vec{\mu}_2 = \mu_2 \otimes \mu_0$ and µs denotes a Gaussian measure on periodic distributions given by
\begin{align}
d \mu_s
= Z_s^{-1} e^{-\frac{1}{2} \| u\|_{{H}^{s}}^2} du
& = Z_s^{-1} \prod_{n \in \mathbb{Z}^4}
e^{-\frac{1}{2} \langle n \rangle^{2s} |\widehat u(n)|^2}
d\widehat u(n),
\end{align} for 
$s\in\mathbb{R}$. Note that µ 0 corresponds to the white noise measure. More precisely, 
$\vec{\mu}_2$ is defined as the induced probability measure under the map 
$
\omega \in \Omega \longmapsto
(X^1(\omega), X^2(\omega)),$ where 
$X^1(\omega)$ and 
$X^2(\omega)$ are given by
\begin{equation}
X^1(\omega) = \sum_{n \in \mathbb{Z}^4} \frac{g_n(\omega)}{\langle n \rangle^2}e_n
\qquad\text{and}\qquad
X^2(\omega) = \sum_{n \in \mathbb{Z}^4} h_n(\omega)e_n.
\end{equation} Here, 
$\{g_n,h_n\}_{n\in\mathbb{Z}^4}$ denotes a family of independent standard complex-valued Gaussian random variables conditioned so that 
$\overline{g_n}=g_{-n}$ and 
$\overline{h_n}=h_{-n}$, 
$n \in \mathbb{Z}^4$. The main difficulty in making sense of (1.27) comes from the rough support of the base Gaussian measure 
$\vec{\mu}_2$, namely 
$\mathcal{H}^{-\varepsilon}(\mathbb{T}^4) \setminus \mathcal{H}^0 (\mathbb{T}^4)$ for any ɛ > 0; see [Reference Burq and Tzvetkov11, Lemma B.1]. Since the typical element in the support of 
$\vec{\mu}_2$ is merely a distribution, the term 
$\int_{\mathbb{T}^4} u^{k+1} \,dx $ in (1.27) is ill-defined and a renormalization is needed in rigorously constructing 
$\vec{\rho}$.
Similarly to the local theory for SNLB (1.1), where we introduced a renormalization based on the logarithmically diverging variance of Ψ in (1.7), here the same difficulty appears due to the roughness of the support of 
$\vec{\mu}_2$. In fact, for 
$N\in\mathbb{N}$, the typical element X 1 in the support of µ 2 satisfies
\begin{align}
\alpha_N \stackrel{\textrm{def}}{=} \mathbb{E}\big[\big(\pi_N X^1(x)\big)^2\big]
=\sum_{\substack{n\in\mathbb{Z}^4\\|n|\leq N}}\frac{1}{\langle n \rangle^4}
\sim \log N ,
\end{align} which is independent of both 
$t\in\mathbb{R}_+$ and 
$x\in\mathbb{T}^4$. We then define the Wick renormalized truncated potential energy
\begin{align}
R_N(u) =
\frac{1}{k+1}
\int_{\mathbb{T}^4} \mathcal{W}_\alpha ( (\pi_N u)^{k+1} ) dx,
\end{align} where the Wick-ordered power 
$\mathcal{W}_\alpha ((\pi_N u)^{k+1}) $ is defined by
\begin{align}
\mathcal{W}_\alpha ( (\pi_N u)^{k+1}(t,x) )
\, \stackrel{\textrm{def}}{=} H_{k+1} (\pi_N u(t,x); \alpha_N).
\end{align} One can show that 
$\{R_N \}_{N \in \mathbb{N}}$ forms a Cauchy sequence in 
$L^p(\mu_2)$ for any finite 
$p \geq 1$, from which we conclude that there exists a limiting random variable R(u) given by
\begin{align}
\lim_{N\to\infty} R_N(u) \stackrel{\textrm{def}}{=} R(u) =
\frac{1}{k+1}
\int_{\mathbb{T}^4}
\mathcal{W}_\alpha (u^{k+1}(x) ) dx.
\end{align} See [Reference Oh and Thomann46, Proposition 1.1] and [Reference Li, Oh and Zheng32, Proposition 3.4] for details. We then construct the Gibbs measure 
$\vec{\rho}$ as the limit of the following truncated Gibbs measures
\begin{align}
d\vec{\rho}_N(u,\partial_t u )= Z_N^{-1}e^{-R_N(u)}d\vec{\mu}_2(u,\partial_t u).
\end{align}Proposition 1.5. Given any 
$1 \leq p \lt \infty$, we have
\begin{equation}
\lim_{N\rightarrow\infty} e^{-R_N(u)} = e^{-R(u)}
\qquad \text{in } L^p( \mu_2).
\end{equation} Consequently, the truncated Gibbs measure 
$\vec{\rho}_N$ in (1.34) converges, in the sense of (1.35), to a limiting Gibbs measure 
$\vec{\rho}$ given by
\begin{align}
d\vec{\rho}(u,\partial_t u) = Z^{-1} e^{-R(u)}d\vec{\mu}_2(u,\partial_t u).
\end{align} We now sketch the proof of Proposition 1.5. From an application of Nelson’s estimate, we obtain uniform in N integrability of the truncated density; for any 
$1\le p \lt \infty$,
\begin{equation}
\sup_{N\in \mathbb{N}} \big\| e^{-R_N(u)} \big\|_{L^p( \mu_2)} \lt \infty.
\end{equation} See, for example, [Reference Oh and Thomann46, Proposition 1.2] and [Reference Li, Oh and Zheng32, Proposition 3.6]. Combining the uniform bound (1.37) with a convergence in measure deduced from (1.33), we obtain (1.35); see, for example, [Reference Tzvetkov56, Remark 3.8] and [Reference Li, Oh and Zheng32, (3.32)]. This allows us to construct the Gibbs measure 
$\vec{\rho}$ in (1.36), which is mutually absolutely continuous with respect to the base Gaussian measure 
$\vec{\mu}_2$.
We can now consider the dynamical problem for the 
$\Phi^{k+1}_4$-model (1.24). In particular, we consider the following truncated Wick renormalized SdNLB
\begin{align}
\partial_t^2 u_N + \partial_t u_N +(1-\Delta)^2 u_N
+
\pi_N\big( \mathcal{W}_\alpha ( (\pi_N u)^{k} ) \big)
= \sqrt{2} \xi,
\end{align} and show almost sure global well-posedness and invariance of the Gibbs measure 
$\vec{\rho}$ for the limiting equation:
\begin{align}
\partial_t^2 u + \partial_t u +(1-\Delta)^2 u
+
\mathcal{W}_\alpha (u^k)
= \sqrt{2} \xi.
\end{align}Theorem 1.6 Let 
$k \in 2 \mathbb{N} + 1$.
The Wick renormalized SdNLB (1.39) is almost surely globally well-posed with respect to the Gibbs measure 
$\vec{\rho}$ in (1.36) and the Gibbs measure 
$\vec{\rho}$ is invariant under the dynamics. More precisely, there exists a non-trivial stochastic process 
$(u,\partial_t u)\in C(\mathbb{R}_+;\mathcal{H}^{-\varepsilon}(\mathbb{T}^4))$ for any ɛ > 0 such that, given any T > 0, the solution 
$(u_N,\partial_t u_N)$ to the renormalized truncated SdNLB (1.38) with random initial data 
$(u_N, \partial_t u _N)|_{t = 0}$ distributed according to the truncated Gibbs measure 
$\vec{\rho}_N$ in (1.34), converges in probability to some stochastic process 
$(u,\partial_t u)$ in 
$C([0,T];\mathcal{H}^{-\varepsilon}(\mathbb{T}^4))$. Moreover, the law of 
$(u(t),\partial_t u(t))$ is given by the renormalized Gibbs measure 
$\vec{\rho}$ in (1.36) for any 
$t\ge 0$.
By using Bourgain’s invariant measure argument, due to the convergence of 
$\vec{\rho}_N$ to 
$\vec{\rho}$, Theorem 1.6 follows once we construct the limiting process 
$(u,\partial_t u)$ locally-in-time with a good approximation property for the solution uN to (1.38) and establish invariance of the truncated measures 
$\vec{\rho}_N$ under (1.38). The former follows from adapting the proof of Theorem 1.1 to the damped models (1.38)-(1.39), while the latter exploits the Hamiltonian structure of the truncated system (1.38). See 
$\S$ 5 for details.
(i) Let
$(u, \partial_t u)$ be the limiting process constructed in Theorem 1.3. Then, as a consequence of Bourgain’s invariant measure argument, one can obtain the following logarithmic growth bound (in time):
\begin{align*}
\| (u(t), \partial_t u (t)) \|_{\mathcal{H}^{- \varepsilon}} \leq C(\omega) \big( \log (1 + t) \big)^{\frac{k}{2}},
\end{align*}for any
$t \geq 0$. For details, see [Reference Oh, Robert and Tzvetkov44].(ii) The local well-posedness in Theorem 1.1 can be easily adapted to the Wick renormalized SNLB with damped massive linear part
$\partial_t^2 u + \partial_t u + (1-\Delta)^2u$, which we detail in $\S$ 5. We choose to consider the massive linear part $(1-\Delta)^2$ instead of Δ2 to avoid a problem at the zero-th frequency when constructing the Gibbs measure $\vec{\rho}$, as in [Reference Gubinelli, Koch, Oh and Tolomeo24, Reference Oh and Thomann46].(iii) The Gaussian measure µ 2 is the log-correlated Gaussian free field on
$\mathbb{T}^4$ studied in [Reference Oh, Seong and Tolomeo45], and thus the SdNLB dynamics (1.24) are associated with this log-correlated Gibbs measure. Our construction of $\vec{\rho}$ in (1.36) is valid for $k\in 2\mathbb{N}+1$ and with a plus sign in front of the potential energy in (1.26). However, in the case of a focusing quartic interaction (i.e., with a minus sign in front of the potential energy and k = 3 in (1.26)), the authors in [Reference Oh, Seong and Tolomeo45] obtained a non-normalizability result for the corresponding measure and established its exact divergence rate; see [Reference Oh, Seong and Tolomeo45, Theorem 1.4].
2. Preliminaries
In this section, we introduce notations and recall basic lemmas. For 
$a, b \gt 0$, we use 
$a \lesssim b$ to denote that there exists a constant C > 0 such that 
$a \leq Cb$. We write 
$a \sim b$ if 
$a \lesssim b$ and 
$b \lesssim a$. When writing the norm of a space-time function, we usually use short-hand notation, such as 
$L_I^q L_x^r = L^q (I; L^r (\mathbb{T}^4))$ for a given time interval 
$I \subset \mathbb{R}_+$. We will also use the notation 
$L^q_T L^r_x = L^q([0,T]; L^r(\mathbb{T}^4))$ for T > 0.
2.1. Deterministic tools
We first introduce some function spaces. For 
$s \in \mathbb{R}$, we define the L 2-based Sobolev space 
$H^s(\mathbb{T}^4)$ via the norm:
\begin{align*}
\| f \|_{H^s}= \big\| \langle n \rangle^s \widehat f (n) \big\|_{\ell_n^2},
\end{align*} where 
$\langle \cdot \rangle=(1+|\cdot|^2)^\frac12$ and 
$\widehat{f}$ denotes the spatial Fourier transform of f. For 
$1 \leq p \leq \infty$, we define the Lp-based Sobolev space 
$W^{s, p} (\mathbb{T}^4)$ via the norm:
\begin{align}
\| f \|_{W^{s, p}} = \big\| \mathcal{F}^{-1} \big( \langle n \rangle^s \widehat f (n) \big) \big\|_{L^p},
\end{align} where 
$\mathcal{F}^{-1}$ denotes the inverse Fourier transform. Note that 
$H^s (\mathbb{T}^4) = W^{s, 2} (\mathbb{T}^4)$.
We now introduce notation for Littlewood-Paley projections. Let 
$\phi : \mathbb{R} \to [0, 1]$ be a smooth bump function with 
$\text{supp}\phi\subset[-\frac 85, \frac 85]$ and 
$\phi \equiv 1$ on 
$[-\frac 54, \frac 54]$. For 
$\xi \in \mathbb{R}^4$, we define
\begin{align*}
\varphi_1 (\xi) = \phi(|\xi|) \qquad \text{and} \qquad
\varphi_{N} (\xi) = \phi \Big( \frac{|\xi|}{N} \Big) - \phi \Big( \frac{2|\xi|}{N} \Big),
\end{align*} for 
$N \geq 2$ a dyadic number. For a dyadic number 
$N \geq 1$, we define the Littlewood-Paley projector PN as the Fourier multiplier operator with the symbol φN. Then,
\begin{equation*}f=\sum_{N\geq1\;\text{ dyadic}}{\mathbf P}_Nf.\end{equation*}We also write
\begin{equation*}{\mathbf P}_{\leq N}f=\sum_{1\leq M\leq N\;\text{ dyadic}}{\mathbf P}_Mf.\end{equation*}Next, we recall the following Christ–Kiselev lemma. For a proof, see [Reference Christ and Kiselev14, Reference Hassell, Tao and Wunsch26].
Lemma 2.1. Let 
$X,Y$ be Banach spaces and 
$K(s,t): X \to Y$ be an operator-valued kernel from X to Y. Suppose that we have the estimate
\begin{equation*}\bigg\| \int_{-\infty}^{t_0} K(s,t) f(s) \, ds \bigg\|_{L^q([t_0, \infty); Y)} \lesssim \| f\|_{L^p(\mathbb{R};X)},\end{equation*} for some 
$ 1 \leq p \lt q \leq \infty$, all 
$t_0\in\mathbb{R}$, and all 
$f\in L^p\big((-\infty, t_0); X\big)$. Then, we have
\begin{equation*}\bigg\| \int_{-\infty}^t K(s,t) f(s) \, ds \big\|_{L^q(\mathbb{R}; Y)} \lesssim \|f\|_{L^p(\mathbb{R};X)}. \end{equation*}Note that the assumption in the above lemma is satisfied in particular if we have
\begin{equation*}\bigg\| \int_{\mathbb{R}} K(s,t) f(s) \, ds \bigg\|_{L^q(\mathbb{R};Y)} \lesssim \|f\|_{L^p(\mathbb{R};X)}.\end{equation*}Lastly, we recall the following product estimates. See for example [Reference Gubinelli, Koch and Oh22, Lemma 3.4].
Lemma 2.2. Let 
$0\leq s \leq 1$.
(i) Suppose that 
$1 \lt p_j,q_j,r \lt \infty$, 
$\frac{1}{p_j} + \frac{1}{q_j} = \frac{1}{r}$, 
$j=1,2$. Then, we have
\begin{align*}
\| \langle \nabla \rangle^{-s} (fg) \|_{L^r(\mathbb{T}^4)} & \lesssim \| f \|_{L^{p1}(\mathbb{T}^4)}\| \langle \nabla \rangle^{s} (g) \|_{L^{q1}(\mathbb{T}^4)} \\
& \quad + \|\langle \nabla \rangle^{-s} (f) \|_{L^{p2}(\mathbb{T}^4)}\|g\| \| \langle \nabla \rangle^{-s} (fg) \|_{L^{q2}(\mathbb{T}^4)}.
\end{align*} (ii) Suppose that 
$1 \lt p,q,r \lt \infty$ satisfy the scaling condition 
$\frac1p + \frac1q = \frac1r + \frac{s}{4}$. Then, we have
\begin{equation*}
\| \langle \nabla \rangle^{-s} (fg) \|_{L^r(\mathbb{T}^4)} \lesssim \| \langle \nabla \rangle^{-s} f \|_{L^p(\mathbb{T}^4)} \| \langle \nabla \rangle^s g \|_{L^q(\mathbb{T}^4)}.
\end{equation*}2.2. Tools from stochastic analysis
In the following, we first review some basic facts on Hermite polynomials. See, for example, [Reference Kuo29, Reference Nualart41]. We define the kth Hermite polynomial 
$H_k(x; \sigma)$ with variance σ > 0 via the following generating function:
\begin{equation}
e^{tx - \frac{1}{2}\sigma t^2} = \sum_{k = 0}^\infty \frac{t^k}{k!} H_k(x;\sigma),
\end{equation} for 
$t, x \in \mathbb{R}$. We list the first few Hermite polynomials for readers’ convenience:
\begin{align*}
& H_0(x; \sigma) = 1,
\qquad
H_1(x; \sigma) = x,
\qquad
H_2(x; \sigma) = x^2 - \sigma, \\
& H_3(x; \sigma) = x^3 - 3\sigma x,
\qquad
H_4(x; \sigma) = x^4 - 6\sigma x^2 +3\sigma^2.
\end{align*} From (2.2), we obtain the following identities for any 
$k \in \mathbb{N}$ and 
$x,y\in\mathbb{R}$:
\begin{align}
H_k(x+y; \sigma ) &=
\sum_{\ell = 0}^k
\binom{k}{\ell}
x^{k - \ell} H_\ell(y; \sigma).
\end{align} We now recall the regularities of the stochastic convolutions and their Wick-powers introduced in 
$\S$ 1. Let Ψ be the stochastic convolution defined in (1.5) and 
${\Psi}^{\mathrm d}$ be the stochastic convolution associated with SdNLB (1.24), namely the solution to the linear stochastic damped beam equation:
\begin{align}
\begin{cases}
\partial_t^2 {\Psi}^{\mathrm d} + \partial_t{\Psi}^{\mathrm d} +(1-\Delta)^2 {\Psi}^{\mathrm d} = \sqrt{2}\xi\\
({\Psi}^{\mathrm d},\partial_t{\Psi}^{\mathrm d})|_{t=0}=(u_0^\omega, u_1^\omega),
\end{cases}
\end{align} with initial data with law 
$\mathcal{L}(u_0^\omega, u_1^\omega) = \vec{\mu}_2$. See 
$\S$ 5 for further details on 
${\Psi}^{\mathrm d}$. Then, using standard stochastic analysis with the Wiener chaos estimate, we have the following regularity and convergence result. For an analogous proof, we refer the readers to [Reference Liu and Oh34, Lemma 2.1]. See also [Reference Gubinelli, Koch and Oh22, Reference Gubinelli, Koch, Oh and Tolomeo24].
Lemma 2.3. Let 
$Z = \Psi$ or 
${\Psi}^{\mathrm d}$, 
$\ell \in \mathbb{N}$, T > 0, and 
$1\le p \lt \infty$. For 
$ \mathcal{W} ( Z_N^\ell )
\, = \, \mathcal{W} ((\pi_N Z)^\ell ) $ denoting the truncated Wick power defined in (1.10) or (1.32), respectively. Then, 
$\{ \mathcal{W} ( Z^\ell_N) \}_{N\in\mathbb{N}}$ is a Cauchy sequence in 
$L^p\big(\Omega; C([0,T];W^{-\varepsilon, \infty}(\mathbb{T}^4)) \big)$. Moreover, denoting the limit by 
$\mathcal{W} (Z^\ell) $, we have 
$\mathcal{W} (Z^\ell ) \in C\big([0,T]; W^{-\varepsilon,\infty}(\mathbb{T}^4)\big)$ almost surely, with the following tail estimate for any 
$1 \le q \lt \infty$, 
$T\ge1$, and λ > 0:
\begin{align*}
P\Big( \| \mathcal{W} (Z^\ell ) \|_{L^q_T W^{-\varepsilon, \infty}_x} \gt \lambda\Big)
\leq C\exp\bigg(-c \frac{\lambda^{\frac{2}{\ell}}}{T^{1 + \frac{2}{q \ell}}}\bigg).
\end{align*} When 
$q = \infty$, we also have the following tail estimate:
\begin{align}
P\Big( \|\mathcal{W} ( Z^\ell ) \|_{L^\infty ([j, j+1]; W^{-\varepsilon, \infty}_x)} \gt \lambda\Big)
\leq C\exp\bigg(-c \frac{\lambda^{\frac{2}{\ell}}}{j+1}\bigg)
\end{align} for any 
$j \in \mathbb{Z}_{\ge 0}$ and λ > 0.
In order to prove Theorem 1.3, we need the following finer regularity property of Ψ. For a proof, see [Reference Gubinelli, Koch, Oh and Tolomeo24, Lemma 2.4].
Lemma 2.4. Let Ψ be as in (1.5) and fix 
$ 0 \lt s \lt 2$. Then, given any 
$x \in \mathbb{T}^4$ and 
$t \in \mathbb{R}_+$, 
$I\Psi(t, x)$ is a mean-zero Gaussian random variable with variance bounded by 
$C_0 t \log N$, where the constant C 0 is independent of 
$x \in \mathbb{T}^4$ and 
$t \in \mathbb{R}_+$.
3. Local well-posedness of SNLB
In this section, we show pathwise local well-posedness of the Wick renormalized SNLB (1.15) in Theorem 1.1. In 
$\S$ 3.1, we show homogeneous and inhomogeneous Strichartz estimates for the linear beam operators appearing in the mild formulation (1.18). We then apply these in 
$\S$ 3.2 to show Theorem 1.1 via a contraction mapping argument.
3.1. Strichartz estimates
To obtain Strichartz estimates for the beam equation, we need the following sharp Strichartz estimates for the linear Schrödinger propagator 
$e^{\pm it\Delta}$ due to Bourgain-Demeter [Reference Bourgain and Demeter7] and Killip-Vişan [Reference Killip and Vişan27].
Lemma 3.1. Let 
$3 \lt p\leq\infty$ and 
$ N\ge1$ be a dyadic integer. Then, we have
\begin{align*}
\big\| e^{\pm it\Delta} \mathbf{P}_{\leq N} f \|_{L^p_{t,x}([0,1]\times \mathbb{T}^4)} \lesssim N^{2-\frac6p} \|f\|_{L^2_x(\mathbb{T}^4)}.
\end{align*} For 
$1\le q,r \lt \infty$, we define the index 
$s_{q,r}$ as follows
\begin{equation}
s_{q,r} \stackrel{\textrm{def}}{=} 2 - \tfrac2q - \tfrac4r.
\end{equation}We then obtain the following estimate.
Lemma 3.2. For 
$3 \lt q\leq r \lt \infty$ and 
$s_{q,r}$ in (3.1), we have
\begin{equation*} \|e^{\pm it\Delta}f\|_{L_t^q([0,1]; L^r_x(\mathbb{T}^4))}\lesssim \|f\|_{H^{s_{q,r}}(\mathbb{T}^4)}.
\end{equation*}Proof. Let 
$I=[0,1]$. We start by writing 
$f = \sum_{N} \mathbf{P}_N f$ from Littlewood-Paley decomposition. From Bernstein’s inequality and Lemma 3.1, we have
\begin{align*}
\|e^{\pm it\Delta} \mathbf{P}_N f \|_{L^q_I L^r_x} & \lesssim N^{\frac4q - \frac4r} \|e^{\pm it\Delta} \mathbf{P}_{\leq 2N} ( \mathbf{P}_N f) \|_{L^q_I L^q_x} \lesssim N^{2 - \frac2q - \frac4r} \|\mathbf{P}_N f \|_{L^2_x}.
\end{align*}Using the Littlewood-Paley theorem, Minkowski’s inequality, and the above estimate, we obtain
\begin{align*}
\| e^{\pm it\Delta} f \|_{L^q_I L^r_x} & \lesssim \bigg( \sum_{N \geq 1 \text{ dyadic}} \|e^{\pm it\Delta} \mathbf{P}_N f \|_{L^q_I L^r_x}^2 \bigg)^{1/2} \\
&\lesssim \bigg(\sum_{N \geq 1 \text{ dyadic}} N^{2s_{q,r}} \|\mathbf{P}_N f \|^2_{L^2_x} \bigg)^{1/2} \\
& \sim \|f\|_{H^{s_{q,r}}}.
\end{align*}as desired.
From the definition of the linear beam operator S(t) in (1.3), the fact that 
$e^{\pm i t \Delta}$ are isometries in 
$H^s(\mathbb{T}^4)$ for any 
$s\in\mathbb{R}$, and Lemma 3.2, we obtain the following homogeneous Strichartz estimates for the linear beam operator.
Lemma 3.3. Let S(t) be the linear operator in (1.3), 
$0 \lt T\le 1$, 
$3 \lt q \leq r \lt \infty$, and 
$s\geq s_{q,r}$. Then, we have
\begin{align*}
\| S(t)(u_0,u_1)\|_{L_T^\infty H_x^s} + \| S(t)(u_0,u_1)\|_{L_T^q L_x^r}&\lesssim \|(u_0,u_1)\|_{\mathcal{H}^s}.
\end{align*} We now establish the following inhomogeneous Strichartz estimate, using a 
$TT^*$-argument.
Lemma 3.4. For 
$3 \lt q\leq r \lt \infty$ and 
$s_{q,r}$ as in (3.1), we have
\begin{align*}
\bigg\|\int_0^t\frac{\sin((t-t') \Delta)}{\Delta}F(t')dt' \bigg\|_{L_t^q([0,1];L^r_x(\mathbb{T}^4))} &\lesssim \|F\|_{L_t^1H_x^{s_{q,r}-2}([0,1]\times\mathbb{T}^4)}.
\end{align*}Proof. Note that the zero-th frequency of F can be estimated easily, and so we can assume that F has mean zero below. Let 
$I=[0,1]$. First note that
\begin{align*}
\bigg\| \int_0^t &\frac{\sin((t-t') \Delta )}{\Delta} F(t') \, dt' \bigg\|_{L^q_I L^r_x} \\
&\lesssim \bigg\| \int_0^t \frac{e^{i(t-t')\Delta}}{\Delta} F(t') \, dt' \bigg\|_{L^q_I L^r_x} + \bigg\| \int_0^t \frac{e^{-i(t-t')\Delta}}{\Delta} F(t') \, dt' \bigg\|_{L^q_I L^r_x}.
\end{align*} Thus, we focus on estimating the first term, as the estimate for the second term follows from an analogous strategy. The operator T defined by 
$Tu_0 = e^{it\Delta}u_0$ is a bounded operator 
$T: H^{s_{q,r}} \to L^q_I L^r_x$ from Lemma 3.2. Note that we have
\begin{align*}
\langle Tu_0,G \rangle_{t,x}& = \int_0^1 \int_{\mathbb{T}^4} e^{it\Delta} u_0(x) \overline{G(t, x)} \, dx \, dt\\
& = \sum_{n \in \mathbb{Z}^4} \widehat{u}_0(n) \overline{\int_0^1 e^{-it|n|^2} \widehat{G}(t,n) \, dt} = \langle u_0, T^* G \rangle_{x},
\end{align*} where the dual operator 
$T^*$ is given by
\begin{equation*}T^* G = \int_0^1 e^{-it \Delta} G(t,\cdot) \, dt,\end{equation*} which in turn is bounded from 
$L^{q'}_I L^{r'}_x$ to 
$H^{-s_{q,r}}_x$. From the trivial boundness of 
$T: H^s_x \to L^\infty_I H^s_x$ for any 
$s\in\mathbb{R}$, we conclude that 
$T^*: L^1_I H^{s_{q,r}}_x \to H^{s_{q,r}}_x$ is also bounded. Consequently, we have that 
$TT^*: L^1_I H^{s_{q,r}}_x \to L^q_t L^r_x$ and
\begin{equation*}TT^* G = \int_0^1 e^{i(t-t')\Delta} G(t', \cdot) \, dt'.\end{equation*}From the Christ–Kiselev lemma (Lemma 2.1), we get that
\begin{equation*}\bigg\| \int_0^t e^{i(t-t')\Delta} G(t', \cdot) \, dt' \bigg\|_{L^q_I L^r_x} \lesssim \|G\|_{L^1_I H^{s_{q,r}}_x},\end{equation*} and by choosing 
$G = \frac{1}{\Delta} F$, we obtain the intended result.
3.2. Proof of Theorem 1.1
In this subsection, we prove Theorem 1.1 by constructing a solution 
$u=\Psi+v$ where Ψ denotes the stochastic convolution solving (1.4) and the remainder v solves (1.12). In particular, we consider the following mild formulation for v:
\begin{align}
v(t) = S(t) (u_0, u_1) \mp \int_0^t \frac{\sin ((t - t') \Delta)}{\Delta} \sum_{\ell = 0}^k \binom{k}{\ell} \Xi_\ell \, v^{k - \ell} (t') dt'
\end{align} for given initial data 
$(u_0, u_1)$ and a source 
$(\Xi_0, \Xi_1, \dots, \Xi_k)$ with the understanding that 
$\Xi_0 \equiv 1$, where S(t) is the linear propagator as defined in (1.3). Given 
$s, \varepsilon \in \mathbb{R}$, we define the space 
$\mathcal{X}^{s,\varepsilon} (\mathbb{T}^4)= \mathcal{H}^s (\mathbb{T}^4) \times \big( C([0, 1]; W^{-\varepsilon, \infty} (\mathbb{T}^4)) \big)^{\otimes k}$ with the following norm for 
$\pmb{\Xi} = (u_0, u_1, \Xi_1, \dots, \Xi_k) \in \mathcal{X}^{s, \varepsilon} (\mathbb{T}^4)$:
\begin{align*}
\| \pmb{\Xi} \|_{\mathcal{X}^{s, \varepsilon}} = \| (u_0, u_1) \|_{\mathcal{H}^s} + \sum_{j = 1}^k \| \Xi_j \|_{C([0, 1]; W^{-\varepsilon, \infty})}.
\end{align*} Moreover, we introduce our solution space 
$X^{s,q,r}(T)$ for 
$s\in\mathbb{R}$ and 
$1 \le q,r\le \infty$:
\begin{align*}
X^{s, q, r} (T) \stackrel{\textrm{def}}{=} C([0, T]; H^s (\mathbb{T}^4)) \cap L^q ([0, T]; L^r (\mathbb{T}^4)).
\end{align*} The local well-posedness in Theorem 1.1 follows from local well-posedness of (3.2) and Lemma 2.3, which states that the random enhanced data set 
$(u_0, u_1, \Psi, \mathcal{W}_\sigma (\Psi^2) , \dots, \mathcal{W}_\sigma (\Psi^k) )$ almost surely belongs to 
$\mathcal{X}^{s, \varepsilon}(\mathbb{T}^4)$ for any ɛ > 0. We then show the following deterministic result for (3.2).
Proposition 3.5. Given an integer 
$k \geq 2$, let s crit be as defined in (1.17). Then, the mild formulation (3.2) is locally well-posed in 
$\mathcal{X}^{s, \varepsilon}(\mathbb{T}^4)$ for
\begin{equation*} \textrm{(i) } k \geq 4:
\ s \geq s_\textrm{crit}
\qquad \text{or} \qquad
\textrm{(ii) } k = 2, 3:
\ s \gt s_\textrm{crit},
\end{equation*}and ɛ > 0 sufficiently small. More precisely, given an enhanced data set
\begin{align*}
\pmb{\Xi} = (u_0, u_1, \Xi_1, \dots, \Xi_k) \in \mathcal{X}^{s, \varepsilon}(\mathbb{T}^4),
\end{align*} there exist 
$T = T (\pmb{\Xi}) \in (0, 1]$ and a unique solution v to the mild formulation (3.2) in the class 
$X^{s', q, r} (T)$ for 
$s' = \min (s, 2 - \varepsilon)$ and for some appropriate 
$1 \le q,r\le \infty$.
Proof. We define the map Γ by
\begin{align}
\begin{split}
\Gamma [v] (t) &\stackrel{\textrm{def}}{=} S(t) (u_0, u_1) \mp \int_0^t \frac{\sin ((t - t') \Delta)}{\Delta} \sum_{\ell = 0}^k \binom{k}{\ell} (\Xi_\ell v^{k - \ell} )(t') dt' \\
&\stackrel{\textrm{def}}{=} S(t) (u_0, u_1) \mp \sum_{\ell = 0}^k \binom{k}{\ell} \mathcal{I} \big( \Xi_\ell v^{k - \ell} \big) (t),
\end{split}
\end{align}and consider the following three cases.
Case 1: 
$k \geq 4$ and 
$s \gt s_{\text{crit}}$.
Let ɛ > 0 sufficiently small and 
$(q, r) = (k - 1 + \theta, 2k - 2)$ for θ > 0 such that 
$s'\ge s_{q,r} \gt s_{\text{crit}}$ for 
$s_{q,r}$ in (3.1). For 
$\ell = 0$, by Lemma 3.4, Sobolev’s inequality and Hölder’s inequality, we obtain
\begin{align}
\big\| \mathcal I \big( v^k \big) \big\|_{X^{s', q,r} (T)}
&\lesssim \|v^k\|_{L^1_T H^{s'-2}_x} \lesssim \|v^k\|_{L^1_T L^{\frac{4}{4-s'}}_x} \lesssim \|v\|_{L^\infty_T L^{\frac{4}{2-s'}}_x} \|v^{k-1}\|_{L^1_T L^{2}_x} \nonumber\\
& \lesssim \|v\|_{L^\infty_T H^{s'}_x} \|v\|^{k-1}_{L^{k-1}_T L^{2k-2}_x} \lesssim T^{\eta} \|v\|^k_{X^{s', q,r}(T)}
\end{align} for some η > 0. For 
$1 \leq \ell \leq k - 1$, proceeding as before, with Lemma 2.2 (ii) and Lemma 2.2 (i) repetitively, we obtain
\begin{align}
\big\| \mathcal I \big( \Xi_\ell v^{k - \ell} \big) \big\|_{X^{s', q,r} (T)}
&
= \big\| \langle \nabla \rangle^{-\varepsilon} (\Xi_\ell v^{k-\ell}) \big\|_{L^1_T H^{s'-2+\varepsilon}_x} \nonumber\\
&\lesssim
\big\| \langle \nabla \rangle^{-\varepsilon} (\Xi_\ell v^{k-\ell}) \big\|_{L^1_T L^{\frac{4}{4-s'-\varepsilon }}_x}\nonumber\\
&\lesssim \| \langle \nabla \rangle^{-\varepsilon} \Xi_\ell \|_{L^\infty_T L^{\frac{4}{\varepsilon}}_x } \big\| \langle \nabla \rangle^\varepsilon v^{k-\ell} \big\|_{L^1_T L^{\frac{4}{4-s'-\varepsilon }}_x}\nonumber\\
&\lesssim \| \Xi_\ell \|_{L^\infty_T W^{-\varepsilon, \infty}_x } \| \langle \nabla \rangle^\varepsilon v \|_{L^\infty_T L^{\frac{4}{2-s'+\varepsilon }}_x}
\| v \|^{k-\ell-1}_{L^{k-\ell-1}_T L^{\frac{2(k-\ell-1)}{1-\varepsilon }}_x} \nonumber\\
&\lesssim \| \Xi_\ell \|_{L^\infty_T W^{-\varepsilon, \infty}_x } \| v \|_{L^\infty_T H^{s'}_x }
T^\eta \| v \|^{k-\ell-1}_{L^{q}_T L^{r}_x} \nonumber\\
&\lesssim T^\eta \| \Xi_\ell \|_{L^\infty_T W^{-\varepsilon, \infty}_x } \| v \|^{k-\ell}_{X^{s', q,r}(T)}
\end{align} for some η > 0 and ɛ > 0 sufficiently small. Lastly, for 
$\ell = k$, by Lemma 3.4, since 
$s' \lt 2$, we have
\begin{align}
\begin{split}
\big\| \mathcal I ( \Xi_k ) \big\|_{X^{s', q,r} (T)}
& \lesssim \| \Xi_k \|_{L^1_T H^{s'-2}_x} \lesssim T \| \Xi_k \|_{L_T^\infty W_x^{-\varepsilon, \infty}}.
\end{split}
\end{align}By Lemma 3.3, (3.3), (3.4), (3.5), and (3.6), we have
\begin{align*}
\| \Gamma[v] \|_{X^{s', q,r} (T)}
&\lesssim \| \pmb{\Xi} \|_{\mathcal{X}^{s , \varepsilon}} + T^\eta \big[ \| \pmb{\Xi} \|_{\mathcal{X}^{s , \varepsilon}}^k + \| v \|_{X^{s', q,r} (T)}^k \big].
\end{align*}A straightforward modification of the above steps yields the following difference estimate:
\begin{align*}
\| \Gamma[v_1] - \Gamma[v_2] \|_{X^{s', q,r} (T)} \lesssim T^\eta & \big[ \| \pmb{\Xi} \|_{\mathcal{X}^{s , \varepsilon}}^k + \| v_1 - v_2 \|_{X^{s', q,r} (T)}\nonumber\\ &\quad \big( \| v_1 \|_{X^{s', q,r} (T)}^{k - 1} + \| v_2 \|_{X^{s', q,r} (T)}^{k - 1} \big) \big].
\end{align*} Then, by 
$T = T (\| \pmb{\Xi} \|_{\mathcal{X}^{s, \varepsilon}}) \gt 0$ sufficiently small, the local well-posedness of (3.2) on 
$[0, T]$ follows from a contraction mapping argument.
Case 2: 
$k = 2, 3$ and 
$s \gt s_{\text{crit}} = 0$.
In this case, we take 
$(q,r) = (3+\theta, 3+\theta)$ for 
$0 \lt \theta\le\frac{3s'}{2-s'}$ which guarantees that 
$s' \ge s_{q,r}$. For 
$\ell = 0$, proceeding as in (3.4), we have
\begin{align}
\big\| \mathcal I \big( v^k \big) \big\|_{X^{s', q,r} (T)} & \lesssim \| v^k\|_{L^1_T H^{s'-2}_x} \lesssim \| v^k\|_{L^1_T L^{\frac{4}{4-s'}}_x} \nonumber\\
& \lesssim \|v\|_{L^k_T L_x^{\frac{4k}{4-s'}}} \lesssim T^\eta \|v\|_{X^{s',q,r}(T)},
\end{align} for some η > 0, since 
$\frac{4k}{4-s'} \le \frac{6}{2-s'}$ for 
$k=2,3$. For 
$1 \leq \ell \leq k - 1$, noticing that 
$\frac{2 (k - \ell - 1)}{1 - \varepsilon} \lt 3 + \theta$ for 
$k \leq 3$ and ɛ > 0 sufficiently small, we proceed as in (3.5) to obtain
\begin{align}
\begin{split}
\big\| \mathcal I \big( \Xi_\ell v^{k - \ell} \big) \big\|_{X^{s', q,r} (T)} = T^\eta \| \Xi_\ell \|_{L_T^\infty W_x^{-\varepsilon, \infty}} \| v \|_{X^{s', q,r} (T)}^{k - \ell}
\end{split}
\end{align}for some η > 0. By Lemma 3.3, (3.3), (3.7), (3.8), and (3.6), we have
\begin{align*}
\| \Gamma[v] \|_{X^{s', q,r} (T)} \lesssim \| \pmb{\Xi} \|_{\mathcal{X}^{s , \varepsilon}} + T^\eta \| \pmb{\Xi} \|_{\mathcal{X}^{s , \varepsilon}}^k + T^\eta \| v \|_{X^{s', q,r} (T)}^k.
\end{align*}Similar steps yield a difference estimate and we conclude the argument as in Case 1.
Case 3: 
$k\ge4$ and 
$s = s_{\text{crit}}$.
In this case, we take 
$
(q, r) = (k , \frac{2k(k-1)}{k+1})$ so that 
$s'=s=s_{\text{crit}} = s_{q,r}$ and 
$q,r \gt 3$. By proceeding as in Case 1, the estimates (3.5) and (3.6) hold, but we can only show (3.4) without the gain of T η on the right-hand side. Thus, we have
\begin{align}
\| \Gamma [v] \|_{L^\infty_T H^{s'}_x} &\leq C \| (u_0, u_1) \|_{\mathcal{H}^{s}} + C T^\eta \| \pmb{\Xi} \|_{\mathcal{X}^{s, \varepsilon}} \| v \|_{L^\infty_T H_x^{s'}} \sum_{\ell = 0}^{k - 2} \| v \|_{L_T^{q} L_x^{r}}^\ell \nonumber\\
&\quad + C T \|\pmb\Xi \|_{\mathcal{X}^{s,\varepsilon}} + C \| v \|_{L^q_T L^r_x}^{k} ,
\end{align}
\begin{align}
\| \Gamma [v] \|_{L_T^{q} L_x^{r}} &\leq C \| S(t) (u_0, u_1) \|_{L_T^q L_x^{r}} + C T^\eta \| \pmb{\Xi} \|_{\mathcal{X}^{s, \varepsilon}} \| v \|_{L^\infty_T H_x^{s'}} \sum_{\ell = 0}^{k - 2} \| v \|_{L_T^{q} L_x^{r}}^\ell \nonumber\\
&\quad + C T \|\pmb\Xi \|_{\mathcal{X}^{s,\varepsilon}} + C \| v \|_{L^q_T L^r_x}^{k} ,
\end{align} for some C > 0 and η > 0. We now define the set 
$B_{a, b, T}$ as
\begin{align*}
B_{a, b, T} \stackrel{\textrm{def}}{=} \big\{v \in X^{s', q,r} (T): \| v \|_{L^\infty_T H^{s'}_x} \leq a \text{ and } \| v \|_{L_T^{q} L_x^{r}} \leq b \big\}.
\end{align*} Suppose that 
$\| (u_0, u_1) \|_{\mathcal{H}^s} \leq A$ for some A > 0. We let 
$a = 4CA$ and 
$0 \lt b \leq 1$ small enough such that
\begin{align}
C b^{k} \le \min(\tfrac a 4, \tfrac b 4).
\end{align} By dominated convergence theorem, we can let 
$T = T(u_0, u_1) \gt 0$ be small enough so that
\begin{align}
\| S(t) (u_0, u_1) \|_{L_T^{q} L_x^{r}} \leq \tfrac{b}{4C}.
\end{align}Choosing T smaller, if necessary, we also assume that
\begin{align}
CT^\eta (k - 1) \| \pmb{\Xi} \|_{\mathcal{X}^{s, \varepsilon}} \leq \min (\tfrac{1}{4a}, \tfrac{1}{4b}) \qquad \text{and} \qquad CT \| \pmb{\Xi} \|_{\mathcal{X}^{s, \varepsilon}} \leq \min (\tfrac a 4, \tfrac b4)
\end{align} Combining (3.9), (3.10), (3.11), (3.12), and (3.13), we know that for 
$v \in B_{a, b, T}$, we have
\begin{equation*}\Arrowvert\Gamma\lbrack v\rbrack\Arrowvert_{L_T^\infty H_x^{s'}}\leq a\qquad\text{and}\qquad\Arrowvert\Gamma\lbrack v\rbrack\Arrowvert_{L_T^qL_x^r}\leq b,\end{equation*} so that Γ maps 
$B_{a,b,T}$ to 
$B_{a,b,T}$. By further shrinking b and T if necessary, we can use similar steps to obtain
\begin{align*}
\| \Gamma[v_1] - \Gamma[v_2] \|_{X^{s', q,r} (T)} \leq \tfrac 12 \| v_1 - v_2 \|_{X^{s', q, r} (T)},
\end{align*} so that Γ is a contraction map on 
$B_{a,b,T}$. We can then conclude the proof of local well-posedness of (3.2).
Remark 3.6. Note that in Case 3 above, to extend the argument to cover the critical regularity s = 0 for 
$k=2,3$ (even without the noise terms), we would need to find suitable 
$q,r$ such that 
$s_{q,r}=0$ with 
$s_{q,r}$ in (3.1). However, we can easily see that this requires that q > 3 which implies that r < 3 and vice-versa, thus the Strichartz estimates in Lemmas 3.3–3.4 do not apply. Moreover, since these are derived from the sharp Strichartz estimates in Lemma 3.1 which are known to fail at the endpoint p = 3 [Reference Bourgain4], the argument above is insufficient to reach critical regularity for quadratic and cubic nonlinearities.
4. Pathwise global well-posedness of the cubic SNLB
In this section, we show pathwise global well-posedness of the Wick-ordered cubic SNLB (1.19) via the hybrid argument in [Reference Gubinelli, Koch, Oh and Tolomeo24]. We restrict our attention to 
$0 \lt s \lt 2$, since the result for 
$s\ge 2$ follows from the same argument. In 
$\S$ 4.1, we first show some preliminary estimates involving the I-operator, and establish commutator estimates to control (1.23). We then prove Theorem 1.3 in 
$\S$ 4.2.
4.1. Commutator estimates and other preliminaries
We recall the definition of the I-operator with Fourier multiplier mN in (1.22). In the following, we fix 
$N \in \mathbb{N}$. From the definition of the I-operator and the Littlewood-Paley theorem, we have that
\begin{align}
\|f\|_{H^s}\lesssim \|I f\|_{H^2} &\lesssim N^{2-s} \|f\|_{H^s},
\end{align}
\begin{align}
\|I f\|_{W^{s_0 + s_1, p}} &\lesssim N^{s_1}\|f\|_{W^{s_0, p}},
\end{align} for any 
$s_0 \in \mathbb{R}$, 
$0 \leq s_1 \leq 2- s$, and 
$1 \lt p \lt \infty$. For simplicity, we will use the notations
\begin{align}
f_{\lesssim N} \stackrel{\textrm{def}}{=} \pi_{\frac{N}{3}} f
\qquad \text{and}\qquad f_{ \gt rsim N} \stackrel{\textrm{def}}{=} \pi_{\frac{N}{3}}^\perp f \stackrel{\textrm{def}}{=} f - f_{\lesssim N},
\end{align} where πN denotes the projection onto frequencies 
$\{|n| \le N\}$.
We first go over some basic commutator estimates in the following lemmas.
Lemma 4.1. Let 
$\frac 43 \le s \lt 2$. Then, for 
$k=1,2,3$, we have
\begin{align*}
\|(If)^k - I(f^k)\|_{L^2} \lesssim N^{-2+k(2-s)} \|If\|^k_{H^2}.
\end{align*}Proof. By the definition of the I-operator and (4.3), we have 
$I(f_{\lesssim N}^k) = f_{\lesssim N}^k$ for 
$ k = 1, 2, 3$. Thus, we obtain
\begin{align}
(If)^k - I(f^k)
&= \big(I(f_{\lesssim N} + f_{ \gt rsim N})\big)^k - I\big((f_{\lesssim N} + f_{ \gt rsim N})^k\big) \nonumber\\
&= \big(f_{\lesssim N} + I(f_{ \gt rsim N})\big)^k - I\big((f_{\lesssim N} + f_{ \gt rsim N})^k\big)\nonumber\\
&= f_{\lesssim N}^k - I\big(f_{\lesssim N}^k\big) +
\sum_{j=0}^{k-1} {k \choose j}
\Big(f_{\lesssim N}^j(If_{ \gt rsim N})^{k-j} - I\big(f_{\lesssim N}^j f_{ \gt rsim N}^{k-j}\big)
\Big) \nonumber\\
&= \sum_{j=0}^{k-1} {k \choose j}
\Big(f_{\lesssim N}^j(If_{ \gt rsim N})^{k-j} - I\big(f_{\lesssim N}^j f_{ \gt rsim N}^{k-j}\big)
\Big).
\end{align} We first consider the case when 
$1 \leq j \leq k - 1$. We let 
$1 \lt q \lt \infty$ sufficiently large and δ > 0 small such that 
$\frac{1}{2} = \frac{j}{q} + \frac{1}{2 + \delta}$. Then, by Hölder’s and Sobolev’s inequalities, we have
\begin{align}
\begin{split}
\|f_{\lesssim N}^j (If_{ \gt rsim N})^{k-j}\|_{L^2}
& \le
\|f_{\lesssim N}\|_{L^{q}}^{j}
\|If_{ \gt rsim N}\|_{L^{(2+\delta) (k - j)}}^{k-j} \\
& \lesssim
\|f_{\lesssim N} \|_{H^{2}}^{j}
\|If_{ \gt rsim N}\|_{H^{2 - \frac{4}{(2+\delta) (k - j)}}}^{k-j}\\
& \lesssim N^{-\frac{4}{2+\delta}} \|If\|_{H^2}^k.
\end{split}
\end{align}Similarly, using the boundedness of the multiplier mN and (4.1), we have
\begin{align}
\begin{split}
\big\|I\big(f_{\lesssim N}^j f_{ \gt rsim N}^{k-j}\big)\big\|_{L^2}
& \lesssim \|f_{\lesssim N}^j f_{ \gt rsim N}^{k-j}\|_{L^2} \\
& \le \|f_{\lesssim N}\|_{L^q}^j \|f_{ \gt rsim N}\|_{L^{(2+\delta)(k-j)}}^{k - j} \\
& \lesssim
\|I f_{\lesssim N}\|_{H^2}^j
\| f_{ \gt rsim N} \|_{H^{2 - \frac{4}{(2+\delta) (k - j)}}}^{k-j}\\
& \lesssim
N^{(k-j)(2-s)-\frac{4}{2+\delta}}
\|I f_{\lesssim N}\|_{H^2}^j
\| f_{ \gt rsim N} \|_{H^s}^{k-j}\\
& \lesssim N^{-2 + k(2-s)} \|If \|_{H^2}^k.
\end{split}
\end{align} Here, we have used the fact that 
$2 - \frac{4}{(2 + \delta) (k - j)} \leq s$, which is guaranteed by 
$ \frac 43 \le s \lt 2$. When j = 0, similar estimates to (4.5) and (4.6) hold with 
$q = \infty$ and δ = 0. Therefore, the desired estimate follows from (4.4), (4.5), and (4.6).
Lemma 4.2. Let 
$0 \lt s \lt 2$ and 
$0 \lt \gamma \lt 1$. Given 
$\delta = \delta(s) \gt 0$ sufficiently small, there exist small 
$\gamma_0 = \gamma_0(\delta) \gt 0$ and large 
$p = p(\delta) \gg1$ such that
\begin{align*}
\|(If)(Ig) - I(fg)\|_{L^2}
\lesssim N^{- \frac{1-\gamma}{2} + \delta } \|f\|_{H^{2-\gamma}}\|g\|_{W^{-\gamma_0, p}}
\end{align*} for any sufficiently large 
$N \gg1$.
Proof. By writing 
$f = f_{\lesssim N^{\frac{1}{2}}} + f_{ \gt rsim N^{\frac12}}$ and 
$g= g_{\lesssim N} + g_{ \gt rsim N}$, we have
\begin{align*}
\begin{split}
(If)(Ig) - I(fg)
&= \Big\{(If_{\lesssim N^{\frac12}}) (Ig_{\lesssim N})
- I(f_{\lesssim N^{\frac12}}g_{\lesssim N})\Big\}\\
&
+ \Big\{(If_{\lesssim N^{\frac12}})(Ig_{ \gt rsim N}) -
I(f_{\lesssim N^{\frac12}}g_{ \gt rsim N})\Big\}\\
&
+ (If_{ \gt rsim N^{\frac12}})(Ig) - I(f_{ \gt rsim N^{\frac12}}g)\\
& =: B_1 + B_2 + B_3 + B_4.
\end{split}
\end{align*} Since the Fourier support of 
$f_{\lesssim N^\frac12 } g_{\lesssim N}$ is contained in 
$\{|n| \le \frac23 N\}$, then 
$B_1 \equiv 0$.
For B 2, note that for 
$(n_1,n_2) \in \Lambda_n \stackrel{\textrm{def}}{=} \{(x,y)\in \mathbb{Z}^4 \times \mathbb{Z}^4: \ n= x +y , \, |x| \le \frac{N^\frac12}{3}, \, |y| \gt \frac N 3 \}$, by considering the sub-regions 
$|n_2|\ge 3N$ and 
$|n_2| \lt 3N$, from the mean value theorem and the definition in (1.22), we get
\begin{align*}
|m_N(n_1+n_2)-m_N(n_2)| \lesssim N^{2-s}|n_2|^{-3+s}
|n_1|.
\end{align*} From the above, the fact that 
$m_N(n_1) \equiv 1$ on 
$\Lambda_n$, and Cauchy-Schwarz inequality, we have
\begin{align*}
\begin{split}
\| B_2\|_{L^2}
& = \bigg\|\sum_{(n_1, n_2) \in \Lambda_n }
\big(m(n_2)-m(n_1+n_2)\big)\widehat f (n_1)\widehat g(n_2)\bigg\|_{\ell^2_n}\\
& \lesssim N^{2-s}
\Bigg\| \sum_{(n_1, n_2) \in \Lambda_n} \frac{1}{\langle n_1 \rangle^{1-\gamma} \langle n_2 \rangle^{3-s-\delta}} \langle n_1 \rangle^{2-\gamma} |\widehat{f}(n_1)| \frac{|\widehat{g}(n_2)|}{\langle n_2 \rangle^\delta}
\Bigg\|_{\ell^2_n}\\
& \lesssim N^{-\frac12 + \delta}
\|f\|_{H^{2-\gamma}}\| g\|_{H^{-\delta} }.
\end{split}
\end{align*}As for B 3, by Hölder’s inequality, Sobolev’s embedding, and (4.2), we have
\begin{align*}
\begin{split}
\|B_3\|_{L^2}
& \le \|If_{ \gt rsim N^{\frac12}}\|_{L^2} \|Ig\|_{L^\infty} \lesssim N^{-1+\frac{\gamma}{2}} \|f \|_{H^{2-\gamma}}\|Ig\|_{W^{5\delta,\delta^{-1}}} \\
& \lesssim N^{-1+\frac{\gamma}{2}+ 6 \delta} \|f\|_{H^{2-\gamma}}\|g\|_{W^{-\delta,\delta^{-1}}}
\end{split}
\end{align*} for 
$\delta = \delta(s) \gt 0$ sufficiently small.
Lastly, by (4.2) and Lemma 2.2 (ii), we have
\begin{align*}
\begin{split}
\|B_4\|_{L^2}
& \lesssim N^{4\delta} \| f_{ \gt rsim N^\frac12} g \|_{H^{-4\delta}} \lesssim N^{4\delta} \| f_{ \gt rsim N^\frac{1}{2}}\|_{H^{4\delta}} \|g\|_{W^{-4\delta , \delta^{-1}}}\\
& \lesssim N^{-1+\frac {\gamma}2 + 6\delta} \| f\|_{H^{2-\gamma}} \|g\|_{W^{-4\delta , \delta^{-1}}}
\end{split}
\end{align*} for 
$\delta = \delta(s) \gt 0$ sufficiently small.
We now show the following commutator estimate using Lemma 4.1 and Lemma 4.2.
Lemma 4.3. Let 
$\frac 32 \lt s \lt 2$ and 
$k = 1, 2$. Given 
$\delta = \delta(s) \gt 0$ sufficiently small, there exist small 
$\gamma_0 = \gamma_0(\delta) \gt 0$ and 
$p = p(\delta) \gg1$ such that
\begin{align*}
\| I(f^k g) - (If)^k I g \|_{L^2}
\lesssim
N^{-\frac{1-k(2-s)}{2} + \delta } \|If\|_{H^2}^k \|g\|_{W^{-\gamma_0 ,p}}
\end{align*} for sufficiently large 
$N \gg1$.
Proof. Using triangle inequality, we have
\begin{align*}
\| I ( f^kg ) - (If)^k Ig\|_{L^2}
&\le \| I(f^kg) -I(f^k)Ig\|_{L^2}
+ \big\|\big(I(f^k) - (If)^k\big)Ig\big\|_{L^2}\\ & \quad =: D_1 + D_2.
\end{align*}For D 1, by Sobolev’s inequality and Lemma 2.2(i), we have
\begin{align}
\|f^{k}\|_{H^{2-k(2-s)}}
\lesssim \|f^k\|_{W^{s,\frac{4}{2+(k-1)(2-s)}}}
\lesssim \|f\|_{H^{s}} \|f\|_{L^{\frac4{2-s}}}^{k-1} \lesssim\|f\|_{H^{s}}^k.
\end{align} Thus, by Lemma 4.2 with 
$\gamma = k(2-s)$, (4.7), and (4.1), there exists δ > 0 sufficiently small such that
\begin{align*}
\begin{split}
\|D_1\|_{L^2}
& \lesssim N^{-\frac{1-k(2-s)}2 + \delta}
\|f\|^k_{H^{s}} \|g\|_{W^{-\gamma_0, p }} \lesssim N^{-\frac{1-k(2-s)}2 + \delta}
\|If\|^k_{H^{2}} \|g\|_{W^{-\gamma_0, p }}
\end{split}
\end{align*} for some small 
$\gamma_0 = \gamma_0 (\delta) \gt 0$ and large 
$p = p (\delta) \gg1 $.
Lastly, by Hölder’s inequality, emma 4.1, Sobolev embedding, and (4.2), we have
\begin{align*}
\begin{split}
\|D_2\|_{L^2} &
\le \|I(f^k) - (If)^k \|_{L^2}
\|Ig\|_{L^\infty} \\
& \lesssim N^{-2+k(2-s)} \|If\|_{H^2}^k \|Ig\|_{W^{5\delta,\delta^{-1}}} \\
& \lesssim N^{-2+k(2-s) + 6\delta} \|If\|_{H^2}^k
\|g\|_{W^{-\delta,\delta^{-1}}}
\end{split}
\end{align*} given that 
$\delta = \delta(s) \gt 0$ is sufficiently small.
We conclude this subsection by showing the following estimates, which will be useful in estimating the second and third terms in (1.23).
Lemma 4.4. (i) Let 
$0 \lt s \lt 2$ and 
$k = 0,1$. Then, for any 
$0 \le \lambda \le 2-s$, we have
\begin{align*}
\bigg|\int_{\mathbb{T}^4} (\partial_t I v(t) )(I v(t))^k I w (t) \, dx \bigg|
\lesssim N^\lambda \big(1 + \big[E(I \vec v )(t)\big]^{\frac34} \big) \| w(t) \|_{W_x^{-\lambda, 4} }
\end{align*} for any 
$t \geq 0$, where E is the energy defined in (1.21).
(ii) We have
\begin{align*}
& \quad \bigg|\int_{t_1}^{t_2} \int_{\mathbb{T}^4} (\partial_t Iv) (Iv)^2 I w \, dx dt \bigg| \\ & \lesssim
\|I w\|_{L^{\eta^{-1}}_{[t_1, t_2]} L^{\eta^{-1}}_x } \int_{t_1}^{t_2} \bigg(E^{\frac{1 + \eta}{1 - 2 \eta}}(I\vec v)(t) + \frac{\eta}{(t-t_1)^\frac12}\bigg) dt
,
\end{align*} uniformly in 
$0 \lt \eta \lt \frac 18$ and 
$t_2 \ge t_1 \ge 0$.
Proof. (i) By Hölder’s inequality, (1.21), and (4.2), we have
\begin{align*}
\bigg|\int_{\mathbb{T}^4} (\partial_t Iv(t)) (Iv(t))^k Iw(t) \, dx\bigg|
& \lesssim\|\partial_t Iv(t)\|_{L^2} \|Iv (t)\|_{L^4}^k \|I w(t)\|_{L^4} \\
& \lesssim N^\lambda \big[ E(I\vec v) (t) \big]^{\frac12 + \frac k4}
\|w(t)\|_{W^{-\lambda , 4}}.
\end{align*}(ii) From Hölder’s inequality, Sobolev inequality, and (1.23), we have
\begin{align}
\bigg|\int_{t_1}^{t_2} \int_{\mathbb{T}^4} (\partial_t Iv) (Iv)^2 I w \, dx dt\bigg|
&\le \int_{t_1}^{t_2}
\|\partial_t Iv\|_{L^2_x}
\|Iv\|_{L^4_x} \|Iv\|_{L^{\frac 4 {1-4\eta}}_x} \|Iw \|_{L_x^{\eta^{-1}}} dt \nonumber\\
& \leq \int_{t_1}^{t_2} \big[ E(I \vec v )(t)\big]^{\frac34} \|I v\|_{L_x^{\frac{4}{1-4\eta}}} \| I w\|_{L_x^{\eta^{-1}}} \, dt
\end{align}for η > 0. By Sobolev inequality, interpolation, and (1.21), we get
\begin{equation}
\begin{aligned}
\|I v\|_{L^{\frac 4 {1-4\eta}}_x}
&\lesssim \| I v \|_{W^{8\eta, \frac{4}{1 + 4\eta}}_x} \lesssim \|Iv \|_{H^2_x}^{4\eta} \|I v\|_{L^4-x}^{1-4\eta } \le \big[ E(I\vec v) \big]^{\frac{1+4\eta}{4} },
\end{aligned}
\end{equation} uniformly in 
$0 \lt \eta \lt \frac18$. From (4.8), (4.9), and Hölder’s inequality, we get
\begin{align}
\bigg|\int_{t_1}^{t_2} \int_{\mathbb{T}^4} (\partial_t Iv) (Iv)^2 I w \, dx dt\bigg| & \lesssim \bigg( \int_{t_1}^{t_2} \big[ E(I \vec v) \big]^{\frac{1+\eta}{1-\eta}} \, dt \bigg)^{1-\eta} \nonumber\\
& \quad \| I w\|_{L_{[t_1,t_2]}^{\eta^{-1}} L^{\eta^{-1}}_x }.
\end{align}To estimate the first factor in (4.10), let
\begin{align*}
p = \frac{1-\eta}{1-2 \eta}, &&
q = \frac{1}{1-2\eta}, && p' = \frac {1-\eta}{\eta}, &&
q' = \frac{1}{2\eta},
\end{align*} where 
$p',q'$ are the Hölder conjugates of 
$p,q$, respectively. By Hölder’s and Young’s inequalities, we have
\begin{align*}
\bigg(\int_{t_1}^{t_2} f (t) dt \bigg)^{1-\eta}
&\le \bigg(\int_{t_1}^{t_2} |f(t) |^p dt \bigg)^{\frac{1-\eta}{p}} (t_2 - t_1)^{\frac{1-\eta}{p'}}\\
&\le \frac 1{q}\bigg(\int_{t_1}^{t_2} |f(t)|^pdt \bigg)^{\frac{q(1-\eta)}{p}} + \frac{1}{q'} (t_2 - t_1)^{\frac{q'(1-\eta)}{p'}} \\
&= (1-2\eta) \int_{t_1}^{t_2} |f(t) |^{\frac{1-\eta}{1-2 \eta}} \, dt + 2\eta (t_2 - t_1)^\frac12.
\end{align*}Thus, we obtain
\begin{align*}
\bigg(\int_{t_1}^{t_2} E^{\frac{1+\eta}{1-\eta}}(I\vec v)(t)dt \bigg)^{1-\eta}
\lesssim
\int_{t_1}^{t_2} \bigg(E^{\frac{1+\eta}{1-2\eta}}(I\vec v)(t) + \frac{\eta}{(t-t_1)^\frac12}\bigg) dt.
\end{align*}Combining the above estimates gives the intended estimate.
4.2. Proof of Theorem 1.3
In this subsection, we construct a solution to the Wick renormalized cubic SNLB (1.19) on the time interval 
$[0, T]$ for any given 
$T \gg 1$. The argument is based on that in [Reference Gubinelli, Koch, Oh and Tolomeo24].
We first fix 
$\frac32 \lt s \lt 2$, 
$N\gg1$ sufficiently large, and T > 0, and establish an estimate for the growth of the modified energy 
$E(t)=E(I_N\vec v) (t)$ on the time interval 
$[0, T]$. Note that by (1.21) and Hölder’s inequality, we have
\begin{equation}\| Iv\|_{H_x^2}^2=\|Iv\|_{L_x^2}^2+\|\Delta(Iv)\|_{L^2_x}^2\leq2E^\frac12(t)+2E(t).\end{equation}Then, by (1.23), Cauchy-Schwarz inequality, Lemmas 4.1, 4.3–4.4, and (4.11), we have
\begin{align}
E(t_2) - E(t_1)
&\lesssim
\int_{t_1}^{t_2} N^{-2+3(2-s)} \big( 1 + E^2(t) \big) dt \nonumber\\
&\quad + \sum_{k=1}^2 \int_{t_1}^{t_2} N^{-\frac{1-k(2-s)}2 + \delta}
\big( 1 + E^\frac {k+1} 2(t) \big) \| \mathcal{W}_\sigma ( \Psi^{3-k}(t) ) \|_{W^{-\gamma_0, p }_x} dt \nonumber\\
&\quad + \sum_{k=2}^3 \int_{t_1}^{t_2} N^\lambda \big( 1 + E^{\frac34}(t) \big)
\| \mathcal{W}_\sigma ( \Psi^{k}(t)) \!\|_{W_x^{-\lambda, 4}} dt \nonumber\\
&\quad +
\Bigg\{\int_{t_1}^{t_2} \bigg(E^{\frac{1 + \eta}{1 - 2 \eta}}(t) + \frac{\eta}{(t-t_1)^\frac12}\bigg) dt\Bigg\}
\|I \Psi\|_{L^{\eta^{-1}}_{[t_1, t_2]} L^{\eta^{-1}}_x }
\end{align} for any 
$t_2 \geq t_1 \geq 0$, where 
$\gamma_0 = \gamma_0 (\delta) \gt 0$ is sufficiently small, 
$p=p(\delta)\gg1$ sufficiently large, 
$0 \leq \lambda \leq 2 - s$, and 
$0 \lt \eta \lt \frac 18$.
Before proceeding to the iterative argument, we introduce some notations. Given 
$j \in \mathbb{Z}_{\ge 0}$, we define 
$V_j = V_j(\omega)$ by
\begin{align*}
V_j = \max_{k = 1, 2} \| \mathcal{W}_\sigma ( \Psi^{3-k}) \|_{L^\infty_{[j, j+1]}W^{-\gamma_0, p}_x}
+ \max_{k = 0, 1}\|\mathcal{W}_\sigma ( \Psi^{3-k} ) \|_{L^\infty_{[j, j+1]}W^{-\lambda, 4}_x}
\end{align*} and define 
$V = V(\omega)$ by
\begin{align}
e^{V^\frac{1}{3}} = \sum_{j = 0}^\infty e^{-j K} e^{V_j^\frac{1}{3}}
\end{align}for some K > 0 large enough. Note that by applying (2.5) in Lemma 2.3 and letting K > 0 be sufficiently large, we have
\begin{align*}
\mathbb{E}\big[
e^{V^\frac{1}{3}} \big] = \sum_{j = 0}^\infty e^{-j K} \mathbb{E}\big[e^{V_j^\frac{1}{3}}\big]
\leq \sum_{j = 0}^\infty e^{-j K} e^{c (j+1)} \lt \infty,
\end{align*} so that V is almost surely finite. Also, for T > 0, we define 
$M_T = M_T(\omega)$ as follows
\begin{align}
M_T = \max_{k = 1, 2} \| \mathcal{W}_\sigma ( \Psi^{3-k} ) \|_{L^\infty_TW^{-\gamma_0, p}_x}
+ \max_{k = 0, 1}\| \mathcal{W} ( \Psi^{3-k} ) \|_{L^\infty_TW^{-\lambda, 4}_x} .
\end{align} From (4.13) we have that 
$ V_j^\frac{1}{3} \leq V^\frac{1}{3} + j K$, and therefore
\begin{align}
M_T = \max_{j \leq T} V_j \lesssim V + K^3 T^3.
\end{align} Furthermore, we define 
$R = R(\omega)$ by
\begin{align}
R = 1+ \sum_{N = 1}^\infty \sum_{j = 1}^\infty
e^{- j K \log N}
\int_0^j \int_{\mathbb{T}^4}e^{|I_N \Psi (t,x) |} dx dt.
\end{align}Then, by using Lemma 2.4 and taking K > 0 possibly larger, we have
\begin{align*}
\mathbb{E}[ R ]
& = 1+\sum_{N = 1}^\infty \sum_{j = 1}^\infty
e^{- j K \log N}
\int_0^j \int_{\mathbb{T}^4}\mathbb{E} \Big[e^{|I_N \Psi (t,x) |}\Big] dx dt\\
& \lesssim \sum_{N = 1}^\infty \sum_{j = 1}^\infty
e^{- j K \log N}
j e^{c j\log N} \lt \infty.
\end{align*} Therefore, 
$1\le R(\omega) \lt \infty$ almost surely.
In the following, we fix 
$\omega\in\Omega$, where Ω is the full probability set where for all 
$\omega\in\Omega$ we have 
$V(\omega), R(\omega) \lt \infty$, and prove pathwise well-posedness of (1.19) on Ω. We first need the following crucial result.
Proposition 4.5. Let 
$\frac 32 \lt s \lt 2$, 
$T \ge T_0 \gg 1$, and 
$N \in \mathbb{N}$ with N > 10. Let 
$V = V(\omega) \lt \infty$ and 
$R = R(\omega) \lt \infty$ be as in (4.13) and (4.15). Then, there exist 
$0 \lt \alpha\le 2s-3$ and 
$0 \lt \beta \lt \alpha$ such that if
\begin{align}
E(t_0) \le N^\beta
\end{align} for some 
$0 \le t_0 \lt T$, then there exists small 
$\tau = \tau (s, T, K, \omega) \gt 0 $ such that
\begin{align*}
E(t) \le N^\alpha
\end{align*} for any t satisfying 
$t_0 \le t \le \min(T, t_0 + \tau)$.
Proof. By replacing E(t) by 
$E(t) + 1$, we can assume that 
$E(t) \geq 1$. Then, from (4.12) with (4.14), we have
\begin{align}
\begin{split}
E(t) \, - \, E(t_0)
&\lesssim
\, \int_{t_0}^{t} N^{-2+3(2-s)} E^2(t') dt' \\
&\quad + M_T \sum_{k=1}^2 \int_{t_0}^{t} N^{-\frac{1-k(2-s)}2 + \delta}
E^\frac {k+1} 2(t') dt' \\
&\quad + M_T \int_{t_0}^{t} N^\lambda E^{\frac34}(t')
dt' \\
&\quad +
\Bigg\{\int_{t_0}^{t} \bigg(E^{1+c\eta}(t') + \frac{\eta}{(t'-t_0)^\frac12}\bigg) dt'\Bigg\}
\|I \Psi\|_{L^{\eta^{-1}}_{[t_0, t]} L^{\eta^{-1}}_x }
\end{split}
\end{align} for any 
$t \geq t_0$ and for 
$c= \frac{3}{1-2\eta} \gt 0$.
We assume that (4.17) holds for some 
$0 \leq t_0 \lt T$. By the continuity in time of E(t) and (4.17) with 
$\alpha \gt \beta$, there exists 
$t_1 \gt t_0$ sufficiently close to t 0 such that
\begin{align}
\max_{t_0 \leq \tau \leq t} E(\tau) \leq 100 N^{\alpha}
\end{align} for any 
$t_0 \leq t \leq t_1$, where 
$\alpha \gt \beta$ is to be determined later. Note that at this point, t 1 depends on t 0. This issue will be dealt with later.
Let 
$\eta = \frac{1}{n}$ for some 
$n \in \mathbb{N}$. We note from (4.16) and 
$n! \le n^n$ that
\begin{align*}
\| I_N \Psi\|_{L^n_{[t_0, t], x}}^n
& = \int_{t_0}^{t}\int_{\mathbb{T}^4}
| I_N \Psi(x, t)|^n
dx dt
\leq n! \int_0^T
\int_{\mathbb{T}^4}
e^{| I_N \Psi(x, t)|}
dx dt \\
& \le
n! e^{K T \log N} R \leq n^n e^{K T \log N} R.
\end{align*}We now choose
\begin{align*}
n \sim K T \log N + c \log (100 N^\alpha ) \sim K T \log N \gg 1,
\end{align*} where we may have to take 
$K\gg1$ larger.
Then, due to (4.19) and 
$\eta = n^{-1}$, we can estimate the last term on the right-hand side of (4.18) as
\begin{align}
\begin{split}
& \Bigg\{\int_{t_0}^{t} \bigg(E^{1+c\eta}(t') + \frac{\eta}{(t'-t_0)^\frac12}\bigg) dt'\Bigg\}
\|I \Psi\|_{L^{\eta^{-1}}_{[t_0, t]} L^{\eta^{-1}}_x}\\
&
\leq
\int_{t_0}^{t} \bigg(E(t') n e^{\frac 1 n (K T\log N + c \log (100N^\alpha )) } R^{\frac{1}{n}} + \frac{e^{\frac 1n K T \log N} R^\frac{1}{n}}{(t'-t_0)^\frac12}\bigg) dt'\\
&
\lesssim
\int_{t_0}^{t} \bigg( \big[K R T \log N \big]E(t') + \frac{R}{(t'-t_0)^\frac12}\bigg) dt',
\end{split}
\end{align} where we used that 
$R = R(\omega) \geq 1$.
Next, we define F by
\begin{align}
F(t) \stackrel{\textrm{def}}{=} \max_{t_0 \le \tau \le t} E(\tau) - E(t_0) + N^\beta \geq E(t).
\end{align}Then, by (4.19), we have
\begin{align}
N^\beta \leq F(t) \leq 200 N^\alpha
\end{align} for 
$t_0 \le t \le t_1$. In particular, we have 
$
\log F(t) \sim \log N.
$ Moreover, from (4.22), we have
\begin{align}
\begin{cases}
N^{-2+3(2-s)} F^2(t) \lesssim N^{-\alpha}F^2(t) \leq 200 F(t), \\
N^{-\frac{1-2(2-s)}2 + \delta}F^\frac{3}{2}(t)
\lesssim N^{-\frac \alpha2 }F^\frac 32(t) \leq \sqrt{200} F(t),\\
N^{-\frac{1-(2-s)}2 + \delta} F(t) \leq F(t), \\
N^\lambda F^{\frac34}(t) \lesssim
N^\lambda F^{-\frac{1}{4}}(t) F(t)
\le F(t),
\end{cases}
\end{align}provided that
\begin{align}
\alpha \le \min(3s-4,2s-3), \quad
\delta \le \min\big(\tfrac{2s - 3 - \alpha}{2}, \tfrac {s-1}2\big),
\quad \text{and}\quad \lambda \le \tfrac\beta 4,
\end{align}which requires that
\begin{equation*}
s \gt \max\Big(\frac43, \frac32\Big)=\frac32.
\end{equation*}Hence, by (4.21), (4.18), (4.23), (4.20), (4.22), and (4.15), we obtain
\begin{align}
F(t) & - F(t_0) \nonumber\\
& = \max_{t_0 \leq \tau \leq t} E (\tau) - E (t_0) \nonumber\\
& \lesssim
(1 + M_T) \int_{t_0}^{t} F(t') dt'
+
\int_{t_0}^{t} \bigg(K R T F(t') \log F(t') + \frac{R}{(t'-t_0)^\frac12}\bigg) dt'\nonumber\\
& \lesssim
(1 + V + K^3 T^3) \int_{t_0}^{t} F(t') dt'
+
\int_{t_0}^{t} \bigg(K R T F(t') \log F(t') + \frac{R}{(t'-t_0)^\frac12}\bigg) dt'\nonumber\\
& \lesssim
(1 + V + K R T)
\int_{t_0}^{t} F(t') (\log F(t') + K^2 T^2) dt' + 2 R (t - t_0)^{\frac{1}{2}}
\end{align} for any 
$t_0 \le t \le t_1$ such that (4.22) holds. Denoting by 
$C_0 = C_0 (s)$ the implicit constant in (4.25), we define G by
\begin{align}
G(t) = F(t) - 2 C_0 R (t - t_0)^\frac{1}{2}.
\end{align} Let us pick 
$t_*(s,R) \gt 0$ such that
\begin{align}
2 C_0 R (t - t_0)^\frac{1}{2} \ll 1,
\end{align}sufficiently small so that
\begin{align}
F(t) \leq 5^{\frac{\alpha - \beta}{2}} G(t).
\end{align} Then, from (4.26) and (4.28), we get that 
$F(t) \sim G(t)$, which combined with (4.25) gives
\begin{align}
\begin{split}
G(t) - G(t_0)
\le C (1 + V + K R T)
\int_{t_0}^{t} G(t') (\log G(t') + K^2 T^2) dt'
\end{split}
\end{align} for any 
$t_0 \le t \le \min( t_1, t_0 + t_*(s, R))$ and some C > 0.
Now, note that the equation
\begin{equation*} \partial_t H(t) =\widetilde{C} H(t) (\log H(t) + K^2 T^2)\end{equation*} has a solution 
$ H(t) = \exp\big( \exp(\widetilde{C} t) (\log H(0) + K^2T^2) - K^2T^2\big). $ Then, by comparison, we deduce from (4.29) that
\begin{align}
G(t) \leq \exp \Big( e^{C(1 + V + K R T) (t - t_0) } (\log G(t_0) + K^2T^2) - K^2T^2\Big),
\end{align}for some constant C > 0.
Recall from (4.26) and (4.21) that 
$G(t_0) = N^\beta$. Then, if the condition
\begin{align}
e^{C(1 + V + K R T) (t - t_0) } (\beta \log N + K^2 T^2) \le \alpha \log N + K^2 T^2 - \frac{\alpha - \beta}{2} \log 5,
\end{align} holds for 
$t_0 \leq t \leq \min(t_1, t_0 + t_* (s, V, R, T, K))$ (where 
$t_* (s, V, R, T , K) \gt 0$ will be specified later), the bound (4.30) implies
\begin{align}
G(t) \le 5^{\frac{\beta - \alpha}{2}} N^\alpha,
\end{align} for any 
$t_0 \le t \le \min( t_1, t_0 + t_*(s, V, R, T))$. Then, we conclude from (4.21), (4.26), (4.27), and (4.28) that
\begin{align}
E(t) \le F(t) \le N^\alpha,
\end{align} for any 
$t_0 \le t \le \min( t_1, t_0 + t_*(s, V, R, T, K))$. This in turn guarantees that the conditions (4.19) and (4.22) are met. Therefore, by a standard continuity argument, we conclude that the bounds (4.32) and (4.33) hold for any t with 
$t_0 \le t \le t_0 + t_*(s, V, R, T, K)$ such that the condition (4.31) holds.
Finally, let us consider the condition (4.31). Let 
$\alpha = \alpha(s) \gt \beta = \beta(s) $ be such that the conditions in (4.24) hold. Since 
$\alpha \gt \beta$, there exists 
$t_{**}(s, V, R, T, K)$ such that, for 
$0\le \tau\le t_{**}$, we have
\begin{align}
\alpha - e^{C(1 + V + R + KRT) \tau }\beta \geq \frac{\alpha - \beta}{2} \gt 0.
\end{align} Then, since N > 10, by choosing 
$0\le \tau \le t_{**}$ sufficiently small such that
\begin{align}
e^{C(1 + V + R + KRT) \tau } - 1
\leq \frac{\frac{\alpha - \beta}{2} \log N - \frac{\alpha - \beta}{2}\log 5}{K^2T^2},
\end{align} we can guarantee that the condition (4.31) is satisfied for 
$t_0 \le t \le t_0+\tau$, and hence so is (4.33). This concludes the proof of Proposition 4.5.
We now present the proof of Theorem 1.3. Fix 
$\frac 74 \lt s \lt 2$, 
$T\gg1 $, 
$\omega \in \Omega$ such that 
$V = V(\omega) \lt \infty$ and 
$R = R(\omega) \lt \infty$, and let the parameters 
$\alpha, \beta, \tau$ be as in Proposition 4.5.
Fix 
$N_0 \gg 1$ which is to be determined later. Then, for 
$k \in \mathbb{Z}_{\geq 0}$, define an increasing sequence 
$\{N_k\}_{k \in \mathbb{Z}_{\ge 0}}$ by setting
\begin{align}
N_{k } = N_0^{\sigma^{k}},
\end{align} for some 
$\sigma = \sigma(s) \gt 1$ sufficiently large satisfying
\begin{equation}
\alpha \ll \sigma (\beta - 2 (2-s)) \quad \text{and} \quad 2\alpha \ll \sigma \beta,
\end{equation}which implies
\begin{align}
N_{k+1}^{2(2-s)} N_k^\alpha + N_k^{2\alpha} \ll N_{k+1}^\beta.
\end{align} From (4.37) and the assumptions on 
$\alpha,\beta$ in Proposition 4.5, we have
\begin{align*}
2(2-s) \lt \beta \lt \alpha \le 2s-3,
\end{align*}which imposes the further restriction on s:
\begin{align*}
s \gt \max\Big(\frac74, \frac32\Big) = \frac74.
\end{align*} Suppose that for some 
$k \in \mathbb{Z}_{\geq 0}$ and 
$t \geq 0 $, it holds that
\begin{align}
E(I_{N_k} \vec v)(t) \leq N^\alpha_k.
\end{align}Then, by (1.21), (4.1), Sobolev inequality, (4.39), and (4.38), we have
\begin{align}
\begin{split}
E(I_{N_{k+1}}\vec v)(t)
& \lesssim N_{k+1}^{2(2-s)} \|\vec v\|_{\mathcal{H}^{s}}^2 + \|v\|_{H^{1}_x}^4 \\
& \lesssim N_{k+1}^{2(2-s)} \|I_{N_k} \vec v\|_{\mathcal{H}^2_x}^2 + \|I_{N_k} v\|_{H^2_x}^4\\
& \lesssim N_{k+1}^{2(2-s)} E(I_{N_k}\vec v)(t) + E(I_{N_k}\vec v)^2(t) \\
& \lesssim N_{k+1}^{2(2-s)} N_k^\alpha+ N_k^{2\alpha}
\ll N_{k+1}^\beta.
\end{split}
\end{align} We are now ready to implement an iterative argument. Given 
$(u_0, u_1) \in \mathcal{H}^s(\mathbb{T}^4)$, choose 
$N_0 = N_0(u_0, u_1, s) \gg 1$ such that
\begin{align}
E(I_{N_{0}}\vec v) (0) \leq N_{0}^\beta.
\end{align} By applying Proposition 4.5, there exists 
$\tau = \tau(s, T, K, \omega) \gt 0$ such that
\begin{align*}
E(I_{N_{0}}\vec v) (t) \leq N_{0}^\alpha,
\end{align*} for any 
$0 \le t \le \tau$. By (4.39) and (4.40), this then implies
\begin{align*}
E(I_{N_{1}}\vec v) (\tau)\leq N_{1}^\beta.
\end{align*}Applying Proposition 4.5 once again, we in turn obtain
\begin{align*}
E(I_{N_{1}}\vec v) (t) \leq N_{1}^\alpha,
\end{align*} for 
$0 \le t \le 2 \tau$. By (4.39) and (4.40), this then implies
\begin{align*}
E(I_{N_{2}}\vec v) (2\tau) \leq N_{2}^\beta.
\end{align*} By iterating this argument 
$\big[\frac{T}{\tau}\big] + 1$ times, we obtain a solution v to the renormalized cubic SNLB (1.19) on the time interval 
$[0, T]$. Since the choice of 
$T \gg1 $ was arbitrary, this proves global well-posedness of (1.19).
Remark 4.6. From the argument above, we can also establish a growth bound on the Sobolev norm of the solution v to SNLB (1.19). Namely, for 
$T\gg1 $ and with the same choice of parameters, we have
\begin{align*}
\| \vec v(t)\|_{\mathcal{H}^s} \lesssim
\big( 1 + E(I_{N_{k}}\vec v) (t)\big)^\frac{1}{2} \leq N_{k}^\frac{\alpha}{2}
\end{align*} for any 
$ 0\le t \le T $ such that 
$k \tau \leq t \leq (k+1) \tau$ for some 
$k \in \mathbb{Z}_{\ge 0}$. Then, by (4.36), we have
\begin{align}
\| \vec v(t)\|_{\mathcal{H}^s} \lesssim
\exp\big( \tfrac{\alpha}{2} \sigma^k \log N_{0}\big)
\le
\exp\big( \tfrac{\alpha}{2} \log N_{0}
\cdot \exp \big(\tfrac{(\log \sigma) t }{\tau}\big)\big)
\end{align} for 
$0 \leq t \leq T$. Moreover, in view of (4.41), we choose 
$N_0 \in \mathbb{N}$ such that 
$1 + E(I_{N_{0}}\vec v) (0) \sim N_{0}^\beta$, so that by (4.1) and the fact that 
$\beta \gt 2 (2 - s)$, we have
\begin{align}
\log N_0 \sim \log
\big(2 + \| \vec v(0)\|_{\mathcal{H}^s}\big).
\end{align} In order to iteratively apply Proposition 4.5 
$ \frac{T}{\tau}$-many times to reach the target time T, we need to guarantee the condition (4.35). By taking
\begin{align}
\tau \sim_{s, V, R, K} T^{-1},
\end{align} the condition (4.34) holds. Thus, in view of (4.36) with 
$k \sim \frac T\tau$, the condition (4.35) becomes
\begin{align*}
0 \lt C_0
\leq \frac{\frac{\alpha - \beta}{2} \sigma^{T^2} \log N_0 - \frac{\alpha - \beta}{2} \log 5}{K^2 T^2},
\end{align*} which holds true for any sufficiently large 
$T\gg1$. Finally, from (4.42), (4.43), and (4.44), we conclude the following double exponential bound for any 
$t\ge0$
\begin{align*}
\| \vec v(t)\|_{\mathcal{H}^s}
\le
C \exp\Big( c \log
\big(2 + \| \vec v(0)\|_{\mathcal{H}^s}\big)
\cdot e^{C(\omega) t^2}\Big).
\end{align*}5. Almost sure global well-posedness of the hyperbolic $\Phi_4^{k + 1}$-model
In this section, we prove Theorem 1.6, i.e. almost sure global well-posedness of the renormalized SdNLB (1.39) and invariance of the corresponding Gibbs measure (1.36). Due to the convergence of 
$\vec{\rho}_N$ to 
$\vec{\rho}$, the invariance of 
$\vec{\rho}_N$ under the truncated SdNLB dynamics (1.38), and Bourgain’s invariant measure argument [Reference Bourgain5, Reference Bourgain6], Theorem 1.6 follows once we construct the limiting process 
$(u,\partial_t u)$ locally-in-time with a good approximation property for the solution uN to (1.38). Furthermore, since 
$\vec{\rho}$ is mutually absolutely continuous with respect to 
$\vec{\mu}_2$, it suffices to study the renormalized SdNLB (1.38) and (1.39) with the Gaussian random initial data 
$(u_0^\omega, u_1^\omega)$ with 
$\mathcal{L}(u_0^\omega, u_1^\omega) = \vec{\mu}_2$.
We first detail how to adapt the proof of Theorem 1.1 to show local well-posedness of (1.38) and (1.39), uniformly in the truncation N, and then show invariance of the truncated Gibbs measure 
$\vec{\rho}_N$ in (1.34) under the dynamics of the truncated SdNLB (1.38).
As in 
$\S$ 3, to construct solutions for SdNLB (1.38)-(1.39), we proceed with a first order expansion centered around the stochastic convolution 
${\Psi}^{\mathrm d}$ which solves (2.4). By defining the operator 
$\mathcal{D}(t)$ as
\begin{equation*}
\mathcal{D}(t) \stackrel{\textrm{def}}{=} e^{-\frac{t}2}\frac{\sin\big(t [[ \nabla ]]^2\big)}{[[ \nabla ]]^2} \qquad \text{with} \qquad [[ \nabla ]] \stackrel{\textrm{def}}{=} \Big( \langle \nabla \rangle^4 - \frac 14 \Big)^{1/4},
\end{equation*} the stochastic convolution 
${\Psi}^{\mathrm d}$ which solves the stochastic damped linear beam equation in (2.4) can be expressed as
\begin{align*}
{\Psi}^{\mathrm d} (t)
=\partial_t\mathcal{D}(t)u_0^\omega + \mathcal{D}(t)(u_0^\omega+u_1^\omega)+ \sqrt{2}\int_0^t\mathcal{D}(t - t')dW(t'),
\end{align*} where W is a cylindrical Wiener process on 
$L^2 (\mathbb{T}^4)$ as in (1.6). A direct but tedious computation shows that 
${\Psi}^{\mathrm d}_N=\pi_N{\Psi}^{\mathrm d}$ is a mean-zero real-valued Gaussian random variable with variance
\begin{align*}
\mathbb{E} \big[{\Psi}^{\mathrm d}_N(t, x)^2\big] = \mathbb{E}\big[\big(\pi_N{\Psi}^{\mathrm d}(t, x)\big)^2\big]
= \alpha_N
\end{align*} for any 
$t\in \mathbb{R}_+$, 
$x\in\mathbb{T}^4$, and 
$N \in\mathbb{N}$, where αN is as in (1.30). Unlike 
$\sigma_N(t)$ in (1.7), the variance αN is independent of time t. This is due to the fact that the massive Gaussian free field µ 2 is invariant under the dynamics of (2.4).
Let uN be the solution to (1.38) with 
$\mathcal{L}\big((u_N, \partial_t u_N)|_{t=0}\big) = \vec{\mu}_2$. Then, we write uN as
\begin{align}
u_N = v_N + {\Psi}^{\mathrm d}
= (v_N + {\Psi}^{\mathrm d}_N) + \pi_N^\perp {\Psi}^{\mathrm d},
\end{align} where 
$\pi_N^\perp=\text{Id}-\pi_N$. Note that the dynamics of the truncated Wick-ordered SdNLB (1.38) decouple into the linear dynamics for the high frequency part given by 
$\pi_N^\perp {\Psi}^{\mathrm d}$ and the nonlinear dynamics for the low frequency part 
$\pi_N u_N$:
\begin{align}
\partial_t^2 \pi_N u_N + \partial_t \pi_N u_N +(1-\Delta)^2 \pi_N u_N
+
\pi_N\big( \mathcal{W}_\alpha ( (\pi_N u)^{k} ) \big)
= \sqrt{2} \pi_N \xi.
\end{align} Then, by (2.3), the remainder term 
$v_N = \pi_N u_N - {\Psi}^{\mathrm d}_N $ satisfies the following equation:
\begin{align}
\begin{cases}
\partial_t^2 v_N + \partial_t v_N +(1-\Delta)^2 v_N +
\sum\limits_{\ell=0}^k {k\choose \ell} \pi_N \big( \mathcal{W}_\alpha (({\Psi}^{\mathrm d}_N)^\ell ) v_N^{k-\ell}\big)
=0,\\
(v_N,\partial_t v_N)|_{t = 0}=(0,0),
\end{cases}
\end{align} where the Wick power 
$
\mathcal{W}_\alpha (({\Psi}^{\mathrm d}_N)^\ell )
\, \stackrel{\textrm{def}}{=} H_{\ell} ({\Psi}^{\mathrm d}_N; \alpha_N)
$
converges to a limit, denoted by 
$\mathcal{W}_\alpha ( ({\Psi}^{\mathrm d})^\ell ) \,$, in 
$C([0,T];W^{-\varepsilon,\infty}(\mathbb{T}^4))$ for any ɛ > 0 and T > 0, almost surely (and also in 
$L^p(\Omega)$ for any 
$p \lt \infty$); see emma 2.3. Thus, we formally obtain the limiting equation:
\begin{align}
\begin{cases}
\partial_t^2 v + \partial_t v +(1-\Delta)^2v +
\sum\limits_{\ell=0}^k {k\choose \ell} \mathcal{W}_\alpha ( ({\Psi}^{\mathrm d})^\ell ) v^{k-\ell}
=0,\\
(v,\partial_t v)|_{t = 0}=(0,0).
\end{cases}
\end{align} We now detail how to modify the proof of Theorem 1.1 to show local well-posedness of (5.3)–(5.4), uniformly in 
$N\in\mathbb{N}$. Note that v is a solution to (5.4) if and only if 
$w = e^{\frac{t}{2}} v$ satisfies the following equation:
\begin{align*}
\partial_t^2 w + (1 - \Delta)^2 w - \frac 14 w + e^{\frac{t}{2}} \sum_{\ell = 0}^k \binom{k}{\ell} \mathcal{W}_\alpha (({\Psi}^{\mathrm d})^\ell) (e^{-\frac{t}{2}} w)^{k - \ell} = 0.
\end{align*} The terms in the mild formulation corresponding to the w-equation can be treated as in Proposition 3.5, except for the one coming from 
$\frac34 w - 2 \Delta w$ term. However, this term can be viewed as a perturbation thanks to the two degrees of smoothing in the integral Duhamel operator, and the analogue of Proposition 3.5 follows. The same argument allows us to show local well-posedness of (5.3) where the time of existence depends only on the stochastic convolution 
${\Psi}^{\mathrm d}$ and its Wick-powers, but not on 
$N\in\mathbb{N}$.
Now, it remains to show the invariance of the truncated Gibbs measure 
$\vec{\rho}_N$ under the truncated SdNLB dynamics (1.38) in the following proposition. In fact, the rest of the proof of Theorem 1.6 follows from a standard application of Bourgain’s invariant measure argument, whose details we omit. See, for example, [Reference Oh, Robert and Tzvetkov44] for further details.
Proposition 5.1. Let 
$N \in \mathbb{N}$. Then, the truncated SdNLB Eq (1.38) is almost surely globally well-posed with respect to the random initial data distributed by the truncated Gibbs measure 
$\vec{\rho}_N$ in (1.34). Moreover, the truncated Gibbs measure 
$\vec{\rho}_{N}$ (1.34) is invariant under the dynamics of (1.38). More precisely, denoting by uN the global solution to truncated SdNLB Eq (1.38), we have 
$\mathcal{L} (u_N (t), \partial_t u_N (t)) = \vec{\rho}_N$ for any 
$t \in \mathbb{R}_+$.
Proof. The idea of the proof has already appeared in [Reference Gubinelli, Koch, Oh and Tolomeo24, Reference Liu and Oh34, Reference Oh, Okamoto and Tolomeo42, Reference Oh, Robert and Tzvetkov44] and so we only sketch the key steps. Given 
$N \in \mathbb{N}$, we define 
$\vec{\mu}_{2, N}$ and 
$\vec{\mu}^\perp_{2,N}$ to be the marginal probability measures on 
$\pi_N \mathcal{H}^{-\varepsilon} (\mathbb{T}^4)$ and 
$\pi_N^\perp \mathcal{H}^{-\varepsilon} (\mathbb{T}^4)$, respectively. In other words, recalling X 1 and X 2 in (1.29), 
$\vec{\mu}_{2, N}$ and 
$\vec{\mu}^\perp_{2,N}$ are the induced probability measures under the maps 
$\omega\in\Omega\mapsto(\pi_N X^1 (\omega), \pi_N X^2 (\omega))$ and 
$\omega\in\Omega\mapsto (\pi_N^\perp X^1 (\omega), \pi_N^\perp X^2 (\omega))$, respectively. Then, with 
$\vec{\mu}_2 = \vec{\mu}_{2, N} \otimes \vec{\mu}_{2, N}^\perp$ and (1.34), we can write
\begin{align}
\vec{\rho}_N = \vec{\nu}_N \otimes \vec{\mu}_{2, N}^\perp,
\end{align} where the measure 
$\vec{\nu}_N$ is given by
\begin{align*}
d \vec{\nu}_N = Z_N^{-1} R_N (u) d \vec{\mu}_{2, N},
\end{align*}with the density RN as in (1.31).
We recall the decomposition (5.1). Since the high frequency part 
$\pi_N^\perp u_N = \pi_N^\perp {\Psi}^{\mathrm d}$ satisfies
\begin{align}
\partial_t^2 \pi_N^\perp {\Psi}^{\mathrm d} + \partial_t \pi_N^\perp {\Psi}^{\mathrm d} + (1 - \Delta)^2 \pi_N^\perp {\Psi}^{\mathrm d} = \sqrt{2} \pi_N^\perp \xi,
\end{align} the dynamics of 
$\pi_N^\perp{\Psi}^{\mathrm d}$ are linear and thus we can study the evolution of each frequency on the Fourier side to conclude that the Gaussian measure 
$\vec{\mu}_{2, N}^\perp$ is invariant under the dynamics of (5.6). In fact, a tedious but direct computation shows that
\begin{equation*}\mathbb{E} \big[ |\widehat {\Psi}^{\mathrm d} (t, n) |^2\big] = \frac{1}{\langle n \rangle^4} \quad \text{and} \quad \mathbb{E} \big[ |\widehat{\partial_t {\Psi}^{\mathrm d}} (t, n)|^2 \big] = 1,\end{equation*} for any 
$t \in \mathbb{R}_+$ and 
$n \in \mathbb{Z}^4$, so that 
$\mathcal{L} ({\Psi}^{\mathrm d} (t), \partial_t {\Psi}^{\mathrm d} (t)) = \vec{\mu}_{2}$ for any 
$t \in \mathbb{R}_+$.
We now consider the low frequency part 
$\pi_N u_N$, which solves (5.2). Denoting 
$(u_{1, N}, u_{2, N}) = (\pi_N u_N, \partial_t \pi_N u_N)$, we can write (5.2) in the following Ito formulation:
\begin{align}
& d \begin{pmatrix} u_{1, N} \\ u_{2, N} \end{pmatrix}
+ \Bigg\{\begin{pmatrix} 0 & -1 \\ (1 - \Delta)^2 & 0 \end{pmatrix}
\begin{pmatrix} u_{1, N} \\ u_{2, N} \end{pmatrix}
+ \begin{pmatrix} 0 \\ \pi_N \mathcal{W}_\alpha (u_{1, N}^k) \end{pmatrix} \Bigg\} dt \nonumber\\& \quad
= \begin{pmatrix} 0 \\ -u_{2, N} dt + \sqrt{2} \pi_N dW \end{pmatrix}.
\end{align} This shows that the generator 
$\mathcal{L}^N$ of the Markov semigroup for (5.7) can be written as 
$\mathcal{L}^N = \mathcal{L}^N_1 + \mathcal{L}^N_2$, where 
$\mathcal{L}_1^N$ corresponds to the (deterministic) NLB with truncated nonlinearity
\begin{align}
\begin{split}
d \begin{pmatrix}
u_{1,N} \\ u_{2,N}
\end{pmatrix}
+ \Bigg\{
\begin{pmatrix}
0 & -1\\
(1-\Delta)^2 & 0
\end{pmatrix}
\begin{pmatrix}
u_{1,N} \\ u_{2,N}
\end{pmatrix}
+
\begin{pmatrix}
0 \\ \pi_N\mathcal{W}_\alpha ( u_{1,N}^k )
\end{pmatrix}
\Bigg\} dt = 0,
\end{split}
\end{align} while 
$\mathcal{L}_2^N$ corresponds to the Ornstein-Uhlenbeck process:
\begin{align}
d u_{2,N}
=
- u_{2,N} dt + \sqrt 2\pi_N dW.
\end{align} The invariance of 
$ \vec{\nu}_N $ under the dynamics of (5.8) follows from Liouville’s theorem and the conservation of the Hamiltonian 
$E_N(u_{1,N}, u_{2,N})$ under the dynamics of (5.8), where
\begin{align*}
E_N (u_{1,N}, u_{2,N} ) & = \frac{1}{2}\int_{\mathbb{T}^4} | (1-\Delta) u_{1,N} |^2 dx +
\frac{1}{2}\int_{\mathbb{T}^4} (u_{2,N})^2dx \\
& \quad + \frac{1}{k+1} \int \mathcal{W}_\alpha (u_{1,N}^{k+1}) \,dx.
\end{align*} Hence, we have 
$(\mathcal{L}^N_1)^* \vec{\nu}_N = 0$, where 
$(\mathcal{L}^N_1)^*$ denotes the adjoint of 
$\mathcal{L}^N_1$. Regarding (5.9), we recall that the Ornstein-Uhlenbeck process preserves the standard Gaussian measure. Thus, 
$ \vec{\nu}_N $ is also invariant under the dynamics of (5.9), since the measure 
$\vec{\nu}_N$ on the second component is the white noise µ 0 (see (1.28) with s = 0 and projected onto the low frequencies 
${|n|\leq N}$). Hence, we have 
$(\mathcal{L}^N_2)^* \vec{\nu}_N = 0$, and so
\begin{equation*}(\mathcal{L}^N)^* \vec{\nu}_N =
(\mathcal{L}^N_1)^* \vec{\nu}_N + (\mathcal{L}^N_2)^* \vec{\nu}_N = 0.\end{equation*} This shows invariance of 
$\vec{\nu}_N$ under (5.7) and hence under (5.2).
Therefore, invariance of the truncated Gibbs measure 
$\vec{\rho}_N$ in (1.34) under the truncated SdNLB dynamics (1.38) follows from (5.5) and the invariance of 
$\vec{\nu}_N$ and 
$\vec{\mu}_{2, N}^\perp$ under (5.7) and (5.6) respectively.
Acknowledgements
The authors would like to thank Tadahiro Oh for proposing the problem and for his support during the project. The authors also thank Chenmin Sun for helpful discussions. A.C., G.L., and R.L. were supported by the European Research Council (Grant No. 864138 “SingStochDispDyn”). G.L. was also supported by the EPSRC New Investigator Award (Grant No. EP/S033157/1). The authors are also grateful to the anonymous referees for their comments.
Open access
