There are some problems in trying to define a quasi-probability distribution in phase space for quantum fields, à la Wigner.
The foremost complication is given by the lack of a universal reference measure in an infinite dimensional vector (e.g. Hilbert) space. For finite dimensional vector spaces, there is the Lebesgue measure, and thus it makes sense to define a probability (or quasi-probability) distribution with respect to the Lebesgue measure. The Wigner function provides such a way, consistently with taking the classical limit (in the limit in fact the Wigner function becomes a true probability distribution).
The second method cited by the OP is trying to do something of this kind, however by using a mathematically ill-defined or at least ambiguous object in the form of the reference path measure "$\mathcal{D}\phi(x)$. What is the measure chosen? The Gaussian (free) one? Some equilibrium measure for an interacting field theory? These measures behave very differently, and are very tricky to define even in the simplest cases. Even if the measure is unambiguously chosen, it would not be good for all quantum states: roughly speaking, the gaussian free measure is not suited to describe interacting states (e.g. the interacting vacua).
There is however a way to define an object similar, at least in spirit, to the Wigner function also for quantum fields. A quantum state (be it of a theory with finitely or infinitely many degrees of freedom) is a noncommutative probability. The Wigner way of characterizing this noncommutativity is by defining a phase space distribution that is not necessarily positive. Another way, that can be generalized more easily to bosonic theories, is to look at Fourier transforms of measures. Let us remain finite dimensional for the moment. The well-known Bochner theorem states that there is a bijection between finite measures and the continuous functions of positive type whose value in zero coincides with the total mass of the measure. The bijection map is given by the Fourier transform (here written directly in phase space terms for convenience):
$$ \hat{\mu}(\xi,\eta)= \int_{\mathbb{R}^{2d}} e^{i(\xi\cdot q + \eta\cdot p)} \mathrm{d}\mu(q,p)\;. $$
If one now quantizes the character $e^{i(\xi\cdot q + \eta\cdot p)}$ to an operator, and replaces the integral with a trace w.r.t. a quantum state $\rho_{\hbar}$ (density matrix) obtains a noncommutative Fourier transform: the quantization of the character is the so-called Weyl operator $W_{\hbar}(\xi,\eta)$, thus
$$\hat{\rho}_{\hbar}(\xi,\eta)=\mathrm{tr}[\rho_{\hbar} W_{\hbar}(\xi,\eta)]\; .$$
It turns out that the noncommutative Fourier transform is very nice, in particular behaves very similarly to the Fourier transform of a measure: it is continuous, its value in zero is the trace of the state, and it is of almost positive type. With respect to positive type functions, there is a quantum correction in the form of relative phases of the type $e^{i\hbar \sigma(\xi,\eta)}$, where $\sigma(\xi,\eta)$ is the canonical symplectic form associated to the (dual of the) phase space. This is analogous to the non-positivity of the Wigner function, and the appearence of the symplectic form is a clear signature of the phase space structure. Finally, taking suitably the limit $\hbar\to 0$, there exists a classical measure $\mu$ such that the noncommutative Fourier transform converges to the classical one:
$$\lim_{\hbar\to 0} \hat{\rho}_{\hbar}(\xi,\eta) = \hat{\mu}(\xi,\eta)\; .$$
Therefore, the quantum description correctly reduces to the classical one in the limit $\hbar\to 0$.
The advantage of this formulation with respect to Wigner's one is that the extension to infinitely many degrees of freedom is possible, without (too many) difficulties. Bosonic quantum fields are indeed defined starting from the Weyl operators, and their canonical commutation relations. In QFTs it is customary to take the complexification of the phase space as the reference space (roughly speaking, taking $z=q+ip$ and $\bar{z}=q-ip$ as canonical variables). In addition, as we have seen the character is defined with respect to the dual space with variables $(\xi,\eta)$ (let us call $\zeta = \xi +i\eta$ its complexification). However, while in finite dimensions any vector space is isomorphic to its dual this is not true for infinite dimensional spaces (in addition, there are several possible dual spaces). To distinguish between the phase space of fields, and its dual space, let us call $z$ the (classical) field and $\zeta$ a test function. So the algebra of Weyl operators, encoding the canonical commutation relations, is given by identifying the space of test functions. Such space shall be infinite dimensional, and endowed with a symplectic form that enters in the CCR (as it entered in the noncommutative Fourier transform). This symplectic form is usually defined starting from the canonical poisson bracket between a field and its momentum. In the complexified version, the standard choice is to define an inner product on the infinite dimensional complex space of test functions, from which the symplectic form is obtained by taking the imaginary part of the inner product (and considering the complex vector space as a real space with "twice many" variables). A standard example is to take the test functions in the Schwartz space of rapidly decreasing functions $\mathscr{S}$ (with inner product given by the standard $L^2$ inner product), and the tempered distributions $\mathscr{S}'$ as the (complexified) phase space of classical fields. This choice amounts to the canonical poisson bracket
$$\{\bar{z}(x),z(y)\}= i \delta(x-y)\; .$$
Now, given a bosonic quantum field state $\rho_{\hbar}$, the noncommutative Fourier transform is again defined as
$$\hat{\rho}_{\hbar}(\zeta)= \mathrm{tr}[\rho_{\hbar} W_{\hbar}(\zeta)]\;,$$
the same as for finite dimensions. Again, a noncommutative form of the Bochner theorem holds (with some technical restrictions): there is a bijection, given by the noncommutative Fourier transform, between regular$^{1}$ states and complex-valued functionals of test functions that are continuous when restricted to finite dimensional subspaces, whose value in zero coincides with the trace of the state, and that are of almost positive type. An analogous theorem, in the commutative case, holds for so-called cylindrical measures (an object that is more general than measures, and exists only in infinite dimensional spaces; this object is the correct classical counterpart of a quantum state in QFT). So again, the parallel between classical and noncommutative probability is done effectively. Furthermore, as in the finite dimensional case, taking suitably the limit $\hbar\to 0$ there exists a cylindrical measure $m$ such that
$$\lim_{\hbar\to 0} \hat{\rho}_{\hbar}(\zeta) = \hat{m}(\zeta)\; .$$
Let me conclude with a small remark on fermionic theories. Since for fermions the Weyl operators (exponentials) are not available, one should use the noncommutative version of the moments of a probability distribution (essentially the trace of arbitrary products of smeared fields). In addition, for fermions there are currently no results available concerning the limit $\hbar\to 0$, and recovering a classical field theory.
References.
Noncommutative Fourier Transform and Bochner Theorem:
- I. Segal. Foundations of the theory of dyamical systems of infinitely many
degrees of freedom. II. Canad. J. Math., 13:1–18, 1961.
- M. Merkli. The ideal quantum gas, in Open Quantum Systems I. Lecture Notes in Mathematics. 2006
Classical limit of states in bosonic QFTs: limit of Fourier transforms:
- Z. Ammari, F. Nier. Mean field limit for bosons and infinite
dimensional phase-space analysis.
Ann. Henri Poincaré, 9(8):1503–
1574, 2008.
- M. Falconi. Cylindrical Wigner Measures. Documenta Math. 23:1677–1756, 2018.
1: Regular states are defined rather abstractly, but are essentially the only type of physically relevant quantum states, I won't make a lengthy and probably unnecessary discussion about it here.