Why should field operators satisfy the classical equations of motion?

To quantize a scalar field theory with the action:

$$S=\int \mathrm d^Dx\mathscr{L}(\phi,\partial_\mu\phi)=\int \mathrm dx^0L(\phi,\partial_0\phi)$$

we promote $\phi(\vec{x})$ and $\pi=\frac{\delta L}{\delta(\partial_0\phi)}=\pi(\vec{x})$ to be field operator at fixed time $x^0$ (with the canonical commutation relation $[\phi(\vec{x}),\pi(\vec{y})]=\delta^{D-1}(\vec{x}-\vec{y})$), then we use the classical equation of motion to find $\phi$ and $\pi$ at all time.

What principles lie behind this? Why can we assume that the field operators satisfy the classical equation of motion?

---

By example if we have a classical theory with $S=\int \,\mathrm dtL(\phi,\partial_t\phi)$, momentum $\pi=\frac{\delta L}{\delta(\partial_t\phi)}$ hamiltonian $H=\pi\phi-L$, the classical equation of motion is:

$$\frac{\mathrm d\phi}{\mathrm dt}=\frac{\delta L}{\delta(\partial_t\phi)}$$

When we do the canonical quantization, we promote $[\phi,\pi]=\delta(\vec{x}-\vec{y})$ and the equation of motion for $\phi$ is

$$\frac{\mathrm d\phi}{\mathrm dt}=i[H,\phi]$$

Why should $i[H,\phi]$ has the same form as $\frac{\delta L}{\delta(\partial_t\phi)}$? This fact comes from where?

Answer

How to see it in canonical quantization: All operators $\mathcal{O}$ in a quantum theory fulfill the Heisenberg equations of motion

$$\frac{\mathrm{d}}{\mathrm{d}t}\mathcal{O}(t) = \mathrm{i}[H,\mathcal{O}(t)]$$ where $H$ is the Hamiltonian density and which is exactly the quantum version of the classical Hamiltonian equations of motion. So the fields, as they are operators, indeed must obey the classical equations of motion.

How to see it in path integral quantization: Write $\phi'(x) = \phi(x) + \epsilon\delta(y-x)$ and observe the path integral measure is invariant under this. Expand the integrand to first order in $\epsilon$, and deduce

$$\int\left(\frac{\delta S}{\delta\phi(x)} + J(x)\right)\mathrm{e}^{\mathrm{i}S[\phi]+J\phi}\mathcal{D}\phi = 0$$ which is known as the Schwinger-Dyson equation. Setting $J=0$ gives $\delta S/\delta\phi = 0$ (inside the path integral, which is the PI version of "as an operator equation"), which is exactly the classical equation of motion.