QFT and its notations I know hardly anything about quantum field theory (QFT) but I'm giving a try to understand some ideas of it. As far as I understand, in QFT one is interested in studying measures such as:
\begin{eqnarray}
d\mu(\varphi) \propto e^{-S(\varphi)}d\varphi \tag{1}\label{1}
\end{eqnarray}
where $S$ is a given action and $d\varphi$ an 'a priori' measure. In particular, one of the most interesting actions is given by:
\begin{eqnarray}
S(\varphi) = \int_{\mathbb{R}^{d}}\bigg{[}\frac{1}{2}\langle \varphi, (-\Delta+m^{2})\varphi \rangle + g\varphi(x)^{4} dx\bigg{]}. \tag{2}\label{2}
\end{eqnarray}
In view of (\ref{2}), expressions such as:
\begin{eqnarray}
d\mu(\varphi) \propto e^{-S_{0}(\varphi)}e^{-\int_{\mathbb{R}^{d}}g\varphi(x)^{4}}d\varphi \equiv e^{-g\int_{\mathbb{R}^{d}}g\varphi(x)^{4}dx}d\mu_{G}(\varphi) \tag{3}\label{3}
\end{eqnarray}
are very common in the theory, where the term $e^{-g\int_{\mathbb{R}^{d}}\varphi(x)^{4}dx}$ is usually viewed as a perturbation of a Gaussian measure. This being said, I'd like to clarify something. As far as I understand, the measure $d\mu_{G}$ in (\ref{3}) is a Gaussian measure on $\mathcal{S}'(\mathbb{R}^{d})$ induced by a positive-definite function $e^{-\frac{1}{2}B(f,Gf)}$ where $B$ is a quadratic form induced by a bilinear and continuous map $B$ defined on $\mathcal{S}(\mathbb{R}^{d})\times \mathcal{S}(\mathbb{R}^{d})$  in terms of the Green's function $G$ of the massive Laplacian $-\Delta + m^{2}$. Thus, the representation $-g\int_{\mathbb{R}^{d}}\varphi(x)^{4}dx$ is only formal (what does it mean?), since $\varphi$ is supposed to be a function on $\mathcal{S}'(\mathbb{R}^{d})$. It turns out that expressions such as (\ref{3}) end up being very confusing to non-experienced students like me, once we're often misguided to wrong conclusions or get wrong understandings of the subject. I think, however, that there is a deeper reason for such notations than just being a simpler way of writting things down. I think this has to be with the fact that one is interested in seeing $\varphi(x)$ as 'tempered distributional'-valued random variables, but I'm not sure. So I'd like to clarify these concepts. Why do quantum field theorists use pointwise notation to express distributions and what can one benefit from it?
 A: The quick answer is that (3) written by physicists is not to be taken too seriously by mathematicians. However, it is a statement of a goal or research problem which is to find a rigorous definition/construction of what (3) is trying to say.
One has an injective continuous linear map $\iota:\mathscr{S}(\mathbb{R}^d)\rightarrow \mathscr{S}'(\mathbb{R}^d)$ given by
$$
\iota(\varphi)=\left(f\longmapsto \int_{\mathbb{R}^d}\varphi(x)f(x)\ d^dx\right)
$$
where $f$ is a generic element of $\mathscr{S}(\mathbb{R}^d)$. Via this map it is natural to identify $\mathscr{S}(\mathbb{R}^d)$ with a subset of $\mathscr{S}'(\mathbb{R}^d)$. Moreover, this subset is dense (in fact sequentially dense) in $\mathscr{S}'(\mathbb{R}^d)$. The action $S(\varphi)$ in (2) is perfectly well defined for $\varphi$ in the subset $\mathscr{S}(\mathbb{R}^d)$ but not for $\varphi\in \mathscr{S}'(\mathbb{R}^d)\backslash\mathscr{S}(\mathbb{R}^d)$.
Unfortunately, the Gaussian measure $d\mu_G$ is not supported on $\mathscr{S}(\mathbb{R}^d)$ but on the much bigger space $\mathscr{S}'(\mathbb{R}^d)$.
One can be a bit more precise and work in a weighted Sobolev or Besov space of exponent $\alpha$ but the latter would be negative except for $d=1$.
So one has to introduce a regularization in order to remove the regularization. This is explained in my answer to
A roadmap to Hairer's theory for taming infinities
and to
https://physics.stackexchange.com/questions/372306/what-is-the-wilsonian-definition-of-renormalizability
A: To make an interacting (i.e. not purely quadratic) QFT at all meaningful, you have to impose a regulator. The most transparent regulator in many ways, and the only known regulator that allows to address the full non-perturbative content of the QFT, is a lattice regulator. Replacing $\mathbb{R}^d$ by $\mathbb{Z}^d$ (or a finite subset such as $\Lambda=\mathbb{Z}^d\cap[0;L]^d$ with suitable boundary conditions) makes the path integral measure $\mathrm{d}\varphi=\prod_{x\in\Lambda}\mathrm{d}\varphi_x$ unambiguously defined and removes any problems with pointwise operations on the fields appearing inside $S(\varphi)$. Whether a continuum limit then exists is of course a non-trivial (and generally unsolved) problem.
To answer your explicit question: Quantum field theorists (if such a generalization is permitted), at least when writing expressions such as the action of $\varphi^4$ theory, tend to think not so much in terms of distributions as in terms of fields in the regularized theory, for which pointwise products are unproblematic.
It should also be noted that a QFT is the quantization of a classical field theory, and that is governed by an action that involves pointwise products; I'd consider that sufficient motivation to keep the pointwise notation (which moreover makes the Poincaré invariance of the continuum action apparent, something which a more distribution-centric approach would likely obscure).
