Equations of motion for higher spin quantum fields In the canonical quantization of a massive Lagrangian, we obtain equations of motion via the Euler-Lagrange formalism. As I have read, these can become rather long but decouple into the Klein-Gordon (or Dirac) plus a set of lower-order equations. E.g., for the spin-1 Proca field:
$$ (\Box + m^2) A^\mu - \partial^\mu (\partial^\nu A_\nu) = 0 \Rightarrow (\Box + m^2)A^\mu = 0,~ \partial^\mu A_\mu = 0. $$
As I see it, there is no obvious way of going back, i.e., deriving the coupled equations from the decoupled ones.
When constructing quantum fields of spin $s$ via irreducible representations of the Poincaré group (let's call it the Weinberg approach), we end up with completely symmetric tensors (or spinors, in the half-integer spin case) with $s$ Lorentz indices. As was shown in this answer, it follows that these fields fulfill a number of differential equations, namely:
$$ (\Box + m^2) \Phi_{\mu_1\cdots\mu_s} = 0 \text{ or } (i\not\partial + m) \Phi_{\mu_1\cdots\mu_s} = 0  ~,~~ \partial^{\mu_1} \Phi_{\mu_1\cdots\mu_s} = 0 \text{ or } \gamma^{\mu_1} \Phi_{\mu_1\cdots\mu_s} = 0~,~~ \eta^{\mu_i\mu_j} \Phi_{\mu_1\cdots\mu_s} = 0, $$
depending on whether it is a tensor or tensorial spinor. So here, we directly obtain the decoupled equations of motion.
My question now is:
Are the "coupled" equation of motion an artifact of the Lagrange formalism or is there a way to derive them from the Weinberg approach?
I am curious about this because time-ordered products of the fields do have a special relation to the coupled equations of motion. They are the respective greens function. E.g., for the Proca field:
$$ \Big(\eta^{\mu\nu} (\Box + m^2) - \partial^\mu \partial^\nu \Big) \langle T A_\nu(x) A_\rho(x') \rangle = -i \eta^\mu_\rho \delta(x-x').  $$
This does not necessarily follow from the decoupled equations, even though we have $D(\partial) \langle T A_\nu(x) A_\rho(x') \rangle \propto \delta(x-x')$ for $D = (\Box + m^2), \partial^\mu, ... $.
I am happy about any insight on how these properties are related!

[Edit:]
With some help, I was able to narrow down what I am looking for. Let me try to (hopefully) clarify.
It is always possible to combine the "decoupled" equations such that the resulting equations are equivalent to the ones we started with. The equations of motion resulting from the Lagrangian are special because they allow to describe a deformation with an interaction Lagrangian $\mathcal{L}$ via a simple Dyson-series(?) (due to renormalization, this holds only when coupling linealy to an external potential, more general interactions will deviate), e.g. for spin-0 or spin-1:
$$ (\Box + m^2)\phi_{int} = - \Big( \frac{\partial \mathcal{L}}{\partial \phi}\Big)_{int},~ \Big(\eta^\mu_\nu (\Box + m^2) - \partial^\mu \partial_\nu \Big) A^\nu = - \Big( \frac{\partial \mathcal{L}}{\partial A_\mu}\Big)_{int}. $$
The $_{int}$ here denotes a perturbative series in the propagators, here the time-ordered products (this corresponds to the equation for the propagator of $A^\mu$ I wrote above).
So my question would rather be:
Is there a generalization of this formalism to higher-spin fields? Are there still equations of motion and time-ordered products that act as their Green's functions?
Thanks again for any help!
 A: Several comments.

*

*This is not really a well-defined question. Say I have field equations $\mathcal A\phi=\mathcal B\phi=0$ for some differential operators $\mathcal A,\mathcal B$. This system is obviously equivalent to $(\mathcal A+\mathcal B)\phi=(\mathcal A-\mathcal B)\phi=0$. It is also equivalent to $(\mathcal A\phi)^2+(\mathcal B\phi)^2=0$. Which of these presentations is "coupled" and which is "decoupled"? There is no answer: they are all equivalent. If I write down an arbitrary differential equation, how do I determine whether it is "coupled" or "decoupled"? Until you give a formal definition of what it means to be "(de)coupled", I don't think your question really makes sense. Right now it is merely based on aesthetics. One presentation "looks better" than the other, but this is not a precise notion.


*Most PDEs do not admit a Lagrangian description, and those that do, we in general don't know its Lagrangian to begin with. Figuring out a suitable Lagrangian given a set of differential equations is a rather non-trivial task. And this Lagrangian is never unique. So in general there is no "canonical approach". The equations of motion are well-defined; the Lagrangian is most often not, and when it is it is not unique. In this sense, the "Weinberg" approach is the most general one.


*In the case of arbitrary spin fields, the equations of motion are known as the Fronsdal equations. I think these were first written down in 1978. A valid Lagrangian was obtained a year later, 1979. This Lagrangian was not even guaranteed to exist to begin with. In fact, in the interacting case we still don't know of a Lagrangian description, and I believe most people suspect it doesn't exist (at least in the traditional sense; there are some proposals in terms of "master fields" similar in spirit to String field theory).


*In an arbitrary theory, the propagator always satisfies the equations of motion up to contact terms. As the "coupled" and "decoupled" forms are equivalent, the propagator will satisfy both, so both approaches are equally as natural. Again, I don't think there is a well-defined distinction between what you call "coupled" and "decoupled"; they are literally the same equations, written in different form.


*In the case of Proca you are missing a term, $\eta^{\mu\nu}\partial^2-\partial^\mu\partial^\nu+\eta^{\mu\nu}m^2$. The mass term is essential: by taking a derivative you get $m^2(\partial\cdot \Delta)\sim 0$, so for $m\neq0$ you get the standard transversality condition. In the massless case this condition disappears. But anyway, I am not sure what point you are trying to make regarding the Proca field. The two presentations of the equations of motion are equivalent so they lead to the same propagator, even if one approach is perhaps less transparent in terms of what degrees of freedom are dynamical. The on-shell (Weinberg) approach was designed precisely to make these d.o.f. manifest. (By contrast, the off-shell (Lagrangian) approach was designed to make locality manifest, and also because it is straightforward to quantise).
A: Let's discuss for simplicity only free massive fields. Usually fields can interpolate for different representations of the Poincaré group. Equations of motion are equivalent to projectors on states you want to describe (including the on-shell condition). For instance

*

*Scalar field. It can interpolate only for spin $0$ states, therefore, equation of motion is
\begin{equation}
(\Box+m^2)\phi=0.
\end{equation}


*Vector field. It can interpolate for spin $0$ and spin $1$ states. The corresponding projectors (you can easily convince yourself) are
\begin{equation}
\Pi^L_{\mu \nu} = \frac{\partial_\mu \partial_\nu}{\Box}, ~~\text{and} ~~
\Pi^T_{\mu \nu} = \eta _{\mu \nu}-\Pi^L_{\mu \nu}=\eta _{\mu \nu}-\frac{\partial_\mu \partial_\nu}{\Box}
\end{equation}

*

*To get spin $1$ we thus have to impose
\begin{equation}
\Pi^L_{\mu \nu} A^\nu = 0 ~~ \text{and} ~~ (\Box+m^2)\Pi^T_{\mu \nu} A^\nu = 0.
\end{equation}
Obviously, due to properties of projectors (orthogonality), a linear combination of the two constraints will be equivalent to the two equations. In other words, we can write
\begin{equation}
\left [ \alpha \Pi^L_{\mu \nu} + (\Box+m^2)\Pi^T_{\mu \nu} \right ] A^\nu = 0.
\end{equation}
Is it possible now to obtain this equation from a local Lagrangian? Yes, using the definition of projectors we get for $\alpha = m^2$
\begin{equation}
m^2\frac{\partial_\mu \partial_\nu}{\Box} + (\Box+m^2) \left ( \eta _{\mu \nu} - \frac{\partial_\mu \partial_\nu}{\Box}\right )= \eta_{\mu \nu}(\Box+m^2)-\partial_\mu \partial_\nu,
\end{equation}
leading to
\begin{equation}
0 = (\Box+m^2)A_\mu-\partial_\mu (\partial A) =m^2A_\mu - \partial^\nu F_{\mu \nu},
\end{equation}
which is clearly derivable from the Lagrangian
\begin{equation}
\mathcal L_1 = -\frac{1}{4}F^2+\frac{1}{2}m^2 A^2.
\end{equation}

*To get spin $0$ we need to impose
\begin{equation}
\Pi^T_{\mu \nu} A^\nu = 0 ~~ \text{and} ~~ (\Box+m^2)\Pi^L_{\mu \nu} A^\nu = 0.
\end{equation}
Similarly, considering a linear combination
\begin{equation}
m^2\Pi^T_{\mu \nu} + (\Box+m^2)\Pi^L_{\mu \nu} = m^2\eta_{\mu \nu}+\partial_\mu \partial_\nu,
\end{equation}
resulting in
\begin{equation}
0 = m^2 A_\mu + \partial _\mu (\partial A),
\end{equation}
which can be obtained from
\begin{equation}
\mathcal L_0 = \frac{1}{2}(\partial A)^2-\frac{1}{2}m^2 A^2.
\end{equation}
In all these cases getting the Lagrangian is simple, you just multiply eom by the field itself (it does not spoil the orthogonality).



*You can do the same for say Dirac field which interpolates only for spin $1/2$ states. In this case the Dirac equation is just a constraint (and there is no need to impose $p^2=m^2$).
Similar procedures can be done for any spin.
