|
Abstract
We first show that the Gaussian free field is the stationary one-time (resp. one-point) marginal of the parabolic (resp. elliptic) stochastic quantisation equation. We then decompose it along a scale parameter and take as primary objects the stochastic observables: martingales in the filtration of finite-resolution measurements. As examples we study Wick powers as observables in every dimension, and their products (operator product expansion). Interactions are introduced in Lecture 4. Exercises are collected at the end of each section. These notes follow the slides of Lecture 3.
Introduction 1
Structure 2
Conventions 2
4Stochastic quantisation 2
4.1The GFF as stationary marginal of parabolic stochastic quantisation 2
Langevin dynamics 2
History 2
4.2The GFF as one-point marginal of elliptic stochastic quantisation 4
Exercises 5
5Setting: covariance profiles and observables 6
5.1Covariance profiles 6
5.2Stochastic observables 7
Observables as equivalence classes 7
Random variables as observables 7
Exercises 8
6Observables of the Gaussian free field and their products 9
6.1Wick powers and exponentials 9
Point evaluations 9
Wick powers 9
Exponentials 9
6.2Products of observables 9
Brownian deltas 10
Products of point evaluations 10
Operator product expansion 10
Higher orders 11
Locality and factorisation 11
Outlook 12
Exercises 12
References 12
We start from stochastic quantisation (Section 4): the GFF is the stationary marginal of a Langevin dynamics in an auxiliary time, and the one-point marginal of an elliptic equation in two additional dimensions. This is a first way to describe a Euclidean field by an equation for the field itself.
We then change the point of view (Section 5): instead of a single random distribution, we look at the field observed at every finite resolution, and take as primary objects the stochastic observables of [BCG26]: martingales in the filtration of finite-resolution measurements.
Sections, equations and exercises keep the numbering of the whole course: references to Sections 1–3 are to the notes of Lectures 1–2.
Stochastic quantisation describes a Euclidean field as the marginal of a stochastic dynamics. For the GFF the dynamics are linear and everything can be computed explicitly; we give two versions, parabolic and elliptic.
| \(\displaystyle \partial_t \Phi = - (m^2 - \Delta) \Phi - V' (\Phi) + \sqrt{2} \hspace{0.17em} \xi \qquad \text{on } \mathbb{R}_+ \times \mathcal{M},\) | (4.1) |
with \(\xi\) space-time white noise. This is (parabolic) stochastic quantisation. For \(V = 0\) it is the stochastic heat equation, an infinite-dimensional Ornstein–Uhlenbeck process.
Space-time white noise on \(\mathbb{R} \times \mathcal{M}\) is a centred
Gaussian family \((\xi (u))_{u \in L^2 (\mathbb{R} \times
\mathcal{M})}\) with \(\mathbb{E} [\xi (u) \xi (v)] = \langle u, v
\rangle\), formally \(\mathbb{E} [\xi (t, x) \xi (s, y)] = \delta (t -
s) \delta (x - y)\). Let \(W_t (h) = \xi (\mathbb{1}_{[0, t]} \otimes
h)\), \(t \geqslant 0\), a cylindrical Brownian motion. For a (weak)
solution \(\Phi\) of
| \(\displaystyle \mathcal{L}F (\varphi) = \sum_{i, j} \langle f_i, f_j \rangle \hspace{0.17em} \partial_{ij} g (\ldots) - \sum_i \varphi (Af_i) \hspace{0.17em} \partial_i g (\ldots) .\) | (4.2) |
Proposition
Proof. We may assume \(F = g (\varphi (f_1), \ldots)\) and \(G = k (\varphi (f_1), \ldots)\) with the same \(f_i\). Apply (3.10) of Lecture 2 to the cylinder function \(k \hspace{0.17em} \partial_i g\) with \(h = f_i\) and sum over \(i\): \(\sum_i \mathbb{E} [k \hspace{0.17em} \partial_i g \hspace{0.17em} \varphi (Af_i)] = \sum_{i, j} \langle f_i, f_j \rangle \mathbb{E} [\partial_j k \hspace{0.17em} \partial_i g + k \hspace{0.17em} \partial_{ij} g]\). Insert this into \(\mathbb{E} [G\mathcal{L}F]\).\(\Box\)
So the static characterisation (DS) and the dynamic one (invariance under the Langevin flow) are two faces of the same integration by parts formula.
Theorem
| \(\displaystyle \Phi_t (f) := \sqrt{2} \hspace{0.17em} \xi (K_{t, f}), \qquad K_{t, f} (r, x) := \mathbb{1}_{\{r < t\}} (\mathrm{e}^{- (t - r) A} f) (x) .\) | (4.3) |
\(\Phi\) is a centred Gaussian process with \(\mathbb{E} [\Phi_t (f) \Phi_s (g)] = \langle f, \mathrm{e}^{- |t - s|A} Cg \rangle\). In particular \((\Phi_t)_{t \in \mathbb{R}}\) is stationary and, for every \(t\), \(\Phi_t\) has a version with values in \(\mathcal{S}' (\mathcal{M})\) whose law is \(\mu_m\).
\(\Phi\) is a mild and weak solution of
Let \(\Psi_0\) be an \(\mathcal{S}' (\mathcal{M})\)-valued random variable independent of \(\xi\), and let \(\Psi\) be the mild solution started from \(\Psi_0\),
Then \(\mathbb{E} [\mathrm{e}^{\mathrm{i} \Psi_t (f)}] \to \mathrm{e}^{- \frac{1}{2} \langle f, Cf \rangle}\) as \(t \to \infty\) for every \(f \in \mathcal{S} (\mathcal{M})\). Moreover \(\mu_m\) is invariant (\(\Psi_0 \sim \mu_m\) implies \(\Psi_t \sim \mu_m\) for all \(t\)), and it is the unique invariant probability measure of the dynamics.
Proof. (a) Since \(A\) is self-adjoint and commutes with \(\mathrm{e}^{- uA}\),
In particular \(K_{t, f} \in L^2 (\mathbb{R} \times \mathcal{M})\), and multiplying by \(2\) gives the covariance. Stationarity follows. For fixed \(t\), \(f \mapsto \Phi_t (f)\) is linear and \(\mathbb{E} [\Phi_t (f)^2] = \langle f, Cf \rangle \leqslant m^{- 2} \|f\|_{L^2}^2\), so it is continuous from \(\mathcal{S} (\mathcal{M})\) to \(L^2 (\mathbb{P})\). By the regularisation theorem for random linear functionals on nuclear spaces (Itô–Nawata, Mitoma) it has an \(\mathcal{S}' (\mathcal{M})\)-valued version, which is centred Gaussian with covariance \(C\), i.e. has law \(\mu_m\). On the torus one can see this directly: the coordinates \(\Phi_t (e_n)\) are independent stationary OU processes \(\mathrm{d} X = - \lambda_n X \mathrm{d} t + \sqrt{2} \mathrm{d} B\) with variance \(\lambda_n^{- 1}\), and (3.3) of Lecture 2 applies.
(b) The mild formula is the identity \(K_{t, h} = K_{s, \mathrm{e}^{- (t - s) A} h} + \mathbb{1}_{(s, t)} (r) \mathrm{e}^{- (t - r) A} h\) together with the linearity of \(\xi\). The weak formulation follows from the mild one by a stochastic Fubini argument (Exercise 4.3).
(c) Since \(\Psi_0\) is independent of \(\xi\),
For \(m > 0\), \(\mathrm{e}^{- tA} f \to 0\) in \(\mathcal{S} (\mathcal{M})\) as \(t \to \infty\), so \(\Psi_0 (\mathrm{e}^{- tA} f) \to 0\) a.s. because \(\Psi_0 \in \mathcal{S}' (\mathcal{M})\). By dominated convergence the first factor tends to \(1\), while the second tends to \(\mathrm{e}^{- \frac{1}{2} \langle f, Cf \rangle}\). If \(\Psi_0 \sim \mu_m\), the first factor equals \(\mathrm{e}^{- \frac{1}{2} \langle f, \mathrm{e}^{- 2 tA} Cf \rangle}\) and the product is exactly \(\mathrm{e}^{- \frac{1}{2} \langle f, Cf \rangle}\), which proves invariance. If \(\nu\) is invariant, start from \(\Psi_0 \sim \nu\): then \(\text{Law} (\Psi_t) = \nu\) for all \(t\), so the characteristic functional of \(\nu\) equals that of \(\mu_m\), and \(\nu = \mu_m\).\(\Box\)
|
Figure 4.1. Parabolic stochastic quantisation of the GFF on a circle of length \(8\), \(m = 1\), started from \(\Phi_0 = 0\). Left: a sample of the solution. Right: each Fourier mode is an Ornstein–Uhlenbeck process, whose variance relaxes to \(\lambda_k^{- 1}\) at rate \(2 \lambda_k\) (Theorem 4.2(c)). Dots: Monte Carlo estimates with \(4000\) samples. |
Remark
Parisi and Sourlas [PS79] discovered that the elliptic SPDE in \(d + 2\) dimensions
| \(\displaystyle (m^2 - \Delta_{\mathbb{R}^2 \times \mathcal{M}}) \Psi + V' (\Psi) = \zeta, \qquad \zeta \text{white noise on } \mathbb{R}^2 \times \mathcal{M},\) | (4.4) |
is related to the \(d\)-dimensional Gibbs measure: formally, the law of \(\Psi (z, \cdot)\), \(z \in \mathbb{R}^2\), is \(\propto \exp \left( - 4 \pi \left[ \tfrac{1}{2} \langle \varphi, A \varphi \rangle + \int V (\varphi) \right] \right) \mathcal{D} \varphi\) (dimensional reduction). The underlying mechanism is a hidden supersymmetry [PS79, KLP84]. In [KLP84] the dimensional reduction of the supersymmetric representation is proved rigorously, but not the supersymmetric representation of the solution of the SPDE; the full nonlinear statement was proved for \(d = 0\) in [ADG20] and for some two-dimensional models (polynomial interactions on \(\mathbb{T}^2\), exponential interactions on \(\mathbb{R}^2\)) in [ADG21]; elliptic \(\Phi^4\) equations in \(4\) and \(5\) dimensions are solved in [GH19]. For \(V = 0\) the statement is elementary, and it explains the factor \(4 \pi\). We normalise the noise so that the slice is exactly \(\mu_m\).
Theorem
The equation \(\textsf{L} \Psi = \sqrt{4 \pi} \hspace{0.17em} \zeta\) has a unique \(\mathcal{S}' (\mathbb{R}^2 \times \mathcal{M})\)-valued solution, \(\Psi (u) = \sqrt{4 \pi} \hspace{0.17em} \zeta (\textsf{L}^{- 1} u)\).
For \(z \in \mathbb{R}^2\) and \(f \in \mathcal{S} (\mathcal{M})\) the function \(K_{z, f} := \textsf{L}^{- 1} (\delta_z \otimes f)\) lies in \(L^2 (\mathbb{R}^2 \times \mathcal{M})\), and \(z \mapsto K_{z, f}\) is continuous. The Wiener integrals \(\Psi (z, f) := \sqrt{4 \pi} \hspace{0.17em} \zeta (K_{z, f})\) satisfy \(\Psi (v \otimes f) = \int_{\mathbb{R}^2} v (z) \Psi (z, f) \hspace{0.17em} \mathrm{d} z\) a.s. for \(v \in \mathcal{S} (\mathbb{R}^2)\).
For all \(z \in \mathbb{R}^2\) and \(f, g \in \mathcal{S} (\mathcal{M})\), \(\mathbb{E} [\Psi (z, f) \Psi (z, g)] = \langle f, Cg \rangle\). Hence \(\Psi (z, \cdot)\) has an \(\mathcal{S}' (\mathcal{M})\)-valued version with law \(\mu_m\), and the field \((\Psi (z, \cdot))_{z \in \mathbb{R}^2}\) is stationary and isotropic in \(z\).
Proof. (a) \(\textsf{L}\) is a continuous bijection of \(\mathcal{S} (\mathbb{R}^2 \times \mathcal{M})\), so \(\Psi (u) := \sqrt{4 \pi} \zeta (\textsf{L}^{- 1} u)\) defines an \(\mathcal{S}'\)-valued variable with \(\Psi (\textsf{L} u) = \sqrt{4 \pi} \zeta (u)\). If \(D \in \mathcal{S}'\) satisfies \(D (\textsf{L} u) = 0\) for all \(u \in \mathcal{S}\), then \(D = 0\) since \(\textsf{L}\) is onto: the solution is unique.
(b), (c) Take the Fourier transform in \(z\) only, with measure \(\mathrm{d} q / (2 \pi)^2\). The transform of \(K_{z, f}\) is the \(L^2 (\mathcal{M})\)-valued function \(q \mapsto \mathrm{e}^{- \mathrm{i} q \cdot z} (|q|^2 + A)^{- 1} f\). By Plancherel, the spectral theorem for \(A\) and Tonelli,
using the scalar identity
In particular \(K_{z, f} \in L^2\), and its continuity in \(z\) follows by dominated convergence. Since the transform of \(\int v (z) K_{z, f} \mathrm{d} z\) is \(\hat{v} (q) (|q|^2 + A)^{- 1} f\), the Bochner integral \(\int v (z) K_{z, f} \mathrm{d} z\) equals \(\textsf{L}^{- 1} (v \otimes f)\). The Wiener integral \(\zeta\) is an isometry \(L^2 \to L^2 (\mathbb{P})\) and commutes with Bochner integrals, which gives (b). Multiplying by \(4 \pi\) gives the covariance in (c). The \(\mathcal{S}' (\mathcal{M})\)-valued version is obtained as in Theorem 4.2. Stationarity and isotropy hold because \(\mathbb{E} [\Psi (z, f) \Psi (z', g)] = 4 \pi \int \mathrm{e}^{\mathrm{i} q \cdot (z - z')} \langle f, (|q|^2 + A)^{- 2} g \rangle \frac{\mathrm{d} q}{(2 \pi)^2}\) depends only on \(|z - z' |\).\(\Box\)
Remark
Exercise
Exercise
Exercise
Exercise
In Lectures 1–2 a Euclidean field was a probability measure on \(\mathcal{S}' (\mathcal{M})\), i.e. a random distribution \(\varphi\), and local operators were functions of \(\varphi\). This viewpoint hides the scale-by-scale organisation of the Wilson–Kadanoff picture [WK74], and forces us to work with singular objects from the start. Already for the GFF, the Wick power \(: \varphi^n\): is not a random distribution when \(n (d - 2) \geqslant d\), since its regularisations are not even tight (Proposition 3.11 of Lecture 2).
From now on we index everything by a resolution \(t \in [0, \infty)\). The field observed at resolution \(t\) is a smooth function \(\varphi_t\), and the information available at resolution \(t\) is the \(\sigma\)-algebra \(\mathcal{F}_t\) generated by \((\varphi_s)_{s \leqslant t}\). A physical quantity is then not a single random variable but a consistent family of finite-resolution estimates of it, i.e. a martingale. This is the notion of stochastic observable of [BCG26], on which this lecture is based; see also [Pol84, BK87, Sal99, BBD24] for the renormalisation group side and [BG20, BG23] for the variational side.
In Section 4 the variable \(t\) was the auxiliary time of a stochastic dynamics; from now on \(t\) is a resolution.
Definition
Our standard example is the heat profile
| \(\displaystyle Q_t := C \hspace{0.17em} \mathrm{e}^{- A / t^2}, \qquad \dot{Q}_t = \frac{2}{t^3} \hspace{0.17em} \mathrm{e}^{- A / t^2},\) | (5.1) |
which smooths the Green function with the heat kernel at time \(t^{- 2}\). On \(\mathbb{R}^d\) the kernels are explicit,
| \(\displaystyle \dot{Q}_t (x - y) = \frac{2 t^{d - 3} \mathrm{e}^{- m^2 / t^2}}{(4 \pi)^{d / 2}} \hspace{0.17em} \mathrm{e}^{- t^2 |x - y|^2 / 4},\) | (5.2) |
so \(\dot{Q}_t\) is concentrated on distances \(\lesssim 1 / t\), and in Fourier space \(Q_t\) retains the modes \(|k| \lesssim t\): \(t\) is a momentum scale, or resolution. Other profiles are useful, even when they violate some smoothness requirement of Definition 5.1. A sharp Fourier cutoff is \(\hat{Q}_t (k) = \mathbb{1}_{|k| \leqslant t} (|k|^2 + m^2)^{- 1}\), which is not smooth in \(t\). In a finite range decomposition, \(\dot{Q}_t (x, y) = 0\) for \(|x - y| \geqslant 1 / t\) (Remark 6.4). We write
For \(d \geqslant 2\), \(c_t \to \infty\) as \(t \to \infty\) (Exercise 5.2).
Let \((W_t)_{t \geqslant 0}\) be a cylindrical Brownian motion on \(L^2 (\mathcal{M})\), \(\mathbb{E} [W_t (x) W_s (y)] = (t \wedge s) \delta (x - y)\), and
| \(\displaystyle X_t := \int_0^t \dot{Q}_s^{1 / 2} \hspace{0.17em} \mathrm{d} W_s .\) | (5.3) |
Proposition
So the GFF is the terminal value of a Gaussian martingale of smooth fields: the field at resolution \(t\) (Figure 5.1). The proof is a variant of Theorem 3.2 of Lecture 2 and is left as Exercise 5.1.
Fix a filtration \((\mathcal{F}_t)_{t \geqslant 0}\), for the moment that of \(X\), and write \(\mathbb{E}_t := \mathbb{E} [\hspace{0.17em} \cdot \hspace{0.17em} |\mathcal{F}_t]\).
Definition
Every integrable \(\mathcal{F}_{\infty}\)-measurable random variable \(Y\) gives a closed observable \(O_t =\mathbb{E}_t [Y]\), with \(O_t \to Y\). The point of the definition is that interesting observables need not be closed: a martingale may diverge as \(t \to \infty\) (for instance \(X_t (x)\) for \(d \geqslant 2\)), and still carry a perfectly meaningful family of finite-resolution estimates. The next lemma is elementary but basic.
Lemma
Proof. (a) \(O^1_s - O^2_s =\mathbb{E}_s [O^1_t - O^2_t]\) for every \(t \geqslant s\), and conditional expectation is an \(L^1\) contraction. (b) Likewise, \(\|O_s -\mathbb{E}_s N_t \|_{L^1} \leqslant \|O_t - N_t \|_{L^1}\).\(\Box\)
Example
are observables without a terminal value, and they are not even bounded in \(L^1\): \(\mathbb{E} |B_t^2 - t| = t\mathbb{E} |Z^2 - 1| \to \infty\), \(Z \sim N (0, 1)\). They are the prototypes of the Wick powers of Section 6: at a fixed point, \(X_t (x)\) is a continuous Gaussian martingale with quadratic variation \(c_t\), i.e. a Brownian motion run with the clock \(c_t\). The exponential martingale
makes the loss of terminal information explicit: \(Z_t \to 0\) almost surely, while \(\mathbb{E}Z_t = 1\) for every \(t\). The almost sure limit is zero, the observable is not: the mass escapes to rarer and rarer paths, and the limit forgets all the finite-resolution expectations.
Proposition
is an observable, and it is the only observable with \(\|N_t - O_t \|_{L^1} \to 0\).
Proof. \(O_t = N_0 + M_t +\mathbb{E}_t \left[ \int_0^{\infty} \beta_s \mathrm{d} s \right]\) is a sum of martingales, and \(\|O_t - N_t \|_{L^1} \leqslant \mathbb{E} \int_t^{\infty} | \beta_s | \mathrm{d} s \to 0\). Uniqueness is Lemma 5.4(a).\(\Box\)
The correction is the expected drift still to come. For instance, for \(N_t = B_t + \mathrm{e}^{- t} \sin B_t\) Itô's formula gives \(\beta_t = - \frac{3}{2} \mathrm{e}^{- t} \sin B_t\), and since \(\mathbb{E}_t \sin B_s = \mathrm{e}^{- (s - t) / 2} \sin B_t\),
so the observable in the class of \(N\) is \(O_t = B_t\). This correction is the germ of the forward–backward equations of Lecture 4: there the forecast is obtained from the field in the same way, by adding the expected drift still to come.
Example
Here the resolution runs over \([0, T)\) (reparametrise by \((T - t)^{- 1}\) to get \([T^{- 1}, \infty)\)). The heat semigroup gives \(\mathbb{E}_s \Pi^y_t = \Pi^y_s\), and \(\int g (y) \Pi^y_t \mathrm{d} y =\mathbb{E}_t g (B_T)\) for bounded \(g\). At fixed \(y\), \(\Pi^y_t \to 0\) almost surely as \(t \uparrow T\), but \(\mathbb{E} \Pi^y_t = p_T (y) > 0\). So \(\delta (B_T - y)\) is not a random variable, but it is an observable: the family of its conditional densities. The same construction works for any Markov process with a transition density.
So an observable is determined by its asymptotic behaviour. A germ for \(O\) is a family of functionals \(\mathring{O}_t\) on smooth fields such that \(\|O_t - \mathring{O}_t (X_t)\|_{L^1} \to 0\). The observable field \(O (x)\) is local if it admits germs of the form
Germs are not unique, and this freedom reflects the freedom of choosing counterterms. By Lemma 5.4(b), a germ determines the observable.
Exercise
Exercise
Exercise
Exercise
Proposition
is a continuous martingale, with \(\mathrm{d} \hspace{0.17em} : X^n_t : (x) = n \hspace{0.17em} : X^{n - 1}_t : (x) \hspace{0.17em} \mathrm{d} X_t (x)\) and
| \(\displaystyle \mathbb{E} \left[ : X^n_t : (x) \hspace{0.17em} : X^k_t : (y) \right] = \delta_{nk} \hspace{0.17em} n! \hspace{0.17em} Q_t (x, y)^n .\) | (6.1) |
In particular, for \(x \neq y\), \(\lim_t \mathbb{E} [: X^n_t : (x) : X^n_t : (y)] = n! \hspace{0.17em} G (x - y)^n\).
Proof. By Itô's formula, using \(\mathrm{d} \langle X (x), X (x) \rangle_t = \dot{c}_t \hspace{0.17em} \mathrm{d} t\),
The covariance is Lemma 3.8 of Lecture 2, applied to the Gaussian vector \((X_t (x), X_t (y))\).\(\Box\)
Compare with Theorem 3.9 of Lecture 2: there \(: \varphi^n\): existed as a random distribution only for \(n (d - 2) < d\), and by Proposition 3.11 not at all otherwise. As observables, Wick powers exist in every dimension and need no smearing. When \(n (d - 2) \geqslant d\) and \(0 \leqslant f \mathrel{\not\equiv} 0\), the smeared martingale \(t \mapsto : X_t^n : (f)\) is unbounded in \(L^2\) and, by the Paley–Zygmund argument of Lecture 2, not even tight as \(t \to \infty\): it has no terminal value, and the observable is the whole family, not a limit. The singularity has not disappeared. It has moved from the observables themselves to their products at coincident points, as we now see.
The pointwise product of two martingales is not a martingale. The natural definition of the product keeps the asymptotics.
Definition
Proposition
| \(\displaystyle \mathbb{E} \int_0^{\infty} | \mathrm{d} \langle O^1, O^2 \rangle_s | < \infty .\) | (6.2) |
Then \(O^1 \star O^2\) exists and
| \(\displaystyle (O^1 \star O^2)_t = O^1_t O^2_t +\mathbb{E}_t \left[ \int_t^{\infty} \mathrm{d} \langle O^1, O^2 \rangle_s \right] .\) | (6.3) |
In particular \(\langle O^1 \star O^2 \rangle =\mathbb{E} [O^1_0 O^2_0] +\mathbb{E} \langle O^1, O^2 \rangle_{\infty}\).
Proof. \(P_t := O^1_t O^2_t - \langle O^1, O^2
\rangle_t\) is a martingale on bounded intervals, and \(N_t :=
\mathbb{E}_t \langle O^1, O^2 \rangle_{\infty}\) is a closed martingale
by
Since \(B_t \sim N (0, t)\), \(\mathbb{E}p_{(T - t) / 2} (B_t - \bar{y}) = p_{(T + t) / 2} (\bar{y})\), and, the product being nonnegative, its expectation is its \(L^1\) norm: for \(y \neq z\)
so \(\Pi^y_t \Pi^z_t\) is asymptotically equivalent to the zero observable, and by Definition 6.2 the product exists and \(\Pi^y \star \Pi^z = 0\). At \(y = z\) the exponential suppression is lost, and \(\mathbb{E} (\Pi^y_t)^2 = \kappa_{T - t} p_{(T + t) / 2} (y) \to \infty\): the square does not exist, since a martingale \(O\) with \(\|O_t - (\Pi^y_t)^2 \|_{L^1} \to 0\) would force \(\mathbb{E} (\Pi^y_t)^2 \to \mathbb{E}O_0 < \infty\). The square can, however, be resolved at a fixed coarser resolution. Let \(u = T - \varepsilon\) be a fine resolution and \(s < u\) a coarser one. Since \(p_{\varepsilon} (w)^2 = \kappa_{\varepsilon} p_{\varepsilon / 2} (w)\),
by Taylor expansion in the time variable and \(\partial_{\tau} p_{\tau} = \frac{1}{2} \partial_y^2 p_{\tau}\). This is an operator product expansion: singular coefficients (\(\kappa_{\varepsilon} \to \infty\), while \(\varepsilon \kappa_{\varepsilon}, \varepsilon^2 \kappa_{\varepsilon} \to 0\)) times the observables \(\Pi^y\), \(\partial_y^2 \Pi^y\), \(\ldots\). The order of the limits matters: the remainder is small only after conditioning at a fixed resolution \(s\). Without conditioning, \((\Pi^y_{T - \varepsilon})^2 - \kappa_{\varepsilon} \Pi^y_{T - \varepsilon} = \kappa_{\varepsilon} [p_{\varepsilon / 2} - p_{\varepsilon}] (B_{T - \varepsilon} - y)\) has \(L^1\) norm of order \(\varepsilon^{- 1 / 2}\).
| \(\displaystyle (\textbf{X} (x) \star \textbf{X} (y))_t = X_t (x) X_t (y) + G (x - y) - Q_t (x - y), \qquad \langle \textbf{X} (x) \star \textbf{X} (y) \rangle = G (x - y),\) | (6.4) |
which is the correct two-point function. We never had to discuss the regularity of \(X_{\infty}\). For \(x = y\) and \(d \geqslant 2\) the product does not exist: its expectation would be \(\lim_t \mathbb{E}X_t (x)^2 = \lim_t c_t = \infty\).
| \(\displaystyle (\textbf{X} (x) \star \textbf{X} (y))_t - G (x - y) \hspace{0.27em} \longrightarrow \hspace{0.27em} X_t (x)^2 - c_t = : X^2_t : (x) .\) | (6.5) |
Since both sides are martingales, this reads \(\textbf{X} (x) \star \textbf{X} (y) = G (x - y) \mathbb{1} + : \textbf{X}^2 : (x) + o (1)\), the remainder tending to zero at each fixed resolution \(t\): the Wick square is the regular part of the product, and \(G (x - y)\) is its singular coefficient. More generally (Exercise 6.2), for \(x \neq y\),
where \(: \textbf{X}^p (x) \textbf{X}^q (y) :_t := H_p (X_t (x) ; c_t) H_q (X_t (y) ; c_t)\) minus its contractions at scales \(\leqslant t\). This is Wick's theorem written at the level of observables. The coefficients \(G (x - y)^j\) carry all the singularities, and the observables on the right have regular limits as \(y \to x\).
Subtracting \(G\) leaves the bilocal Wick product \(X_t (x_-) X_t (x_+) - Q_t (x_-, x_+)\), which is smooth in \(\rho\) at fixed \(t\), and the symmetry \(\rho \mapsto - \rho\) kills the odd orders. Each coefficient is a finite sum of Wick bilinears, hence an observable. The constant in the remainder depends on \(t\); the limit \(t \to \infty\) is a separate question. A systematic study of the operator product expansion in this framework is due to P. Paulovics (in preparation).
Remark
Exercise
Exercise
Exercise
Exercise
Exercise
S. Albeverio, F.C. De Vecchi, M. Gubinelli. Elliptic stochastic quantization. Ann. Probab. 48 (2020), 1693–1741.
S. Albeverio, F.C. De Vecchi, M. Gubinelli. The elliptic stochastic quantization of some two dimensional Euclidean QFTs. Ann. Inst. Henri Poincaré Probab. Stat. 57 (2021), 2372–2414.
S. Albeverio, M. Röckner. Stochastic differential equations in infinite dimensions: solutions via Dirichlet forms. Probab. Theory Related Fields 89 (1991), 347–386.
R. Bauerschmidt, T. Bodineau, B. Dagallier. Stochastic dynamics and the Polchinski equation: an introduction. Probab. Surveys 21 (2024), 200–290.
R. Bauerschmidt, D.C. Brydges, G. Slade. Introduction to a Renormalisation Group Method. Lecture Notes in Math. 2242, Springer, 2019.
I. Bailleul, I. Chevyrev, M. Gubinelli. Stochastic observables and Wilson–Itô diffusions. In preparation.
N. Barashkov, M. Gubinelli. A variational method for \(\Phi^4_3\). Duke Math. J. 169 (2020), 3339–3415.
N. Barashkov, M. Gubinelli. On the variational method for Euclidean quantum fields in infinite volume. Probab. Math. Phys. 4 (2023), 761–801.
D.C. Brydges, T. Kennedy. Mayer expansions and the Hamilton–Jacobi equation. J. Stat. Phys. 48 (1987), 19–49.
R. Catellier, K. Chouk. Paracontrolled distributions and the 3-dimensional stochastic quantization equation. Ann. Probab. 46 (2018), 2621–2679.
P.H. Damgaard, H. Hüffel. Stochastic quantization. Phys. Rep. 152 (1987), 227–398.
P.H. Damgaard, H. Hüffel (eds.). Stochastic Quantization. World Scientific, 1988.
G. Da Prato, A. Debussche. Strong solutions to the stochastic quantization equations. Ann. Probab. 31 (2003), 1900–1916.
M. Gubinelli, M. Hofmanová. Global solutions to elliptic and parabolic \(\Phi^4\) models in Euclidean space. Comm. Math. Phys. 368 (2019), 1201–1266.
M. Gubinelli, M. Hofmanová. A PDE construction of the Euclidean \(\Phi^4_3\) quantum field theory. Comm. Math. Phys. 384 (2021), 1–75.
M. Gubinelli, P. Imkeller, N. Perkowski. Paracontrolled distributions and singular PDEs. Forum Math. Pi 3 (2015), e6.
M. Gubinelli. Stochastic quantisation. In Encyclopedia of Mathematical Physics, 2nd ed., Academic Press, 2025, 648–667.
M. Hairer. A theory of regularity structures. Invent. Math. 198 (2014), 269–504.
G. Jona-Lasinio, P.K. Mitter. On the stochastic quantization of field theory. Comm. Math. Phys. 101 (1985), 409–436.
A. Klein, L.J. Landau, J.F. Perez. Supersymmetry and the Parisi–Sourlas dimensional reduction: a rigorous proof. Comm. Math. Phys. 94 (1984), 459–482.
J.-C. Mourrat, H. Weber. The dynamic \(\Phi^4_3\) model comes down from infinity. Comm. Math. Phys. 356 (2017), 673–753.
J. Polchinski. Renormalization and effective Lagrangians. Nuclear Phys. B 231 (1984), 269–295.
G. Parisi, N. Sourlas. Random magnetic fields, supersymmetry, and negative dimensions. Phys. Rev. Lett. 43 (1979), 744–745.
G. Parisi, Y.-S. Wu. Perturbation theory without gauge fixing. Sci. Sinica 24 (1981), 483–496.
M. Salmhofer. Renormalization: An Introduction. Springer, 1999.
K.G. Wilson, J. Kogut. The renormalization group and the \(\epsilon\) expansion. Phys. Rep. 12 (1974), 75–199.