This is part (b) of Schwartz's Problem 14.3 in his Quantum Field Theory and the Standard Model textbook.
Suppose that we have a real scalar field operator $\hat{\Phi}(x^0,\mathbf{x})$ with conjugate momentum field $\hat{\Pi}(x^0,\mathbf{x}) := \partial_0 \hat{\Phi}(x^0,\mathbf{x})$. These operators satisfy the equal-time commutation relations $$ [ \hat{\Phi}(x^0,\mathbf{x}), \hat{\Phi}(x^0,\mathbf{y}) ] = 0 \\ [ \hat{\Pi}(x^0,\mathbf{x}), \hat{\Pi}(x^0,\mathbf{y}) ] = 0 \\ [ \hat{\Phi}(x^0,\mathbf{x}), \hat{\Pi}(x^0,\mathbf{y}) ] = i \delta^{(3)}(\mathbf{x} - \mathbf{y}) $$ At some initial time $t=0$ we can define simultaneous (orthonormal) eigenstates of the Schrodinger-picture operator $\hat{\Phi}(0,\mathbf{x})$ and $\hat{\Pi}(0,\mathbf{x})$ which satisfy $$ \hat{\Phi}(0,\mathbf{x}) | \phi_a \rangle = \phi_a(\mathbf{x}) | \phi_a \rangle \ \ \ \ \mathrm{and} \ \ \ \ \hat{\Pi}(0,\mathbf{x}) | \pi_a \rangle = \pi_a(\mathbf{x}) | \pi_a \rangle $$
Question: The task of Schwartz's problem 14.3(b) here is to use the (third) commutation relation to show that $\hat{\Pi}(0,\mathbf{x})$ acts on eigenstates of $\hat{\Phi}(0,\mathbf{x})$ as the variational derivative $- i \delta/ \delta\phi_{a}(\mathbf{x})$.
My Attempt: I think that this means to show that $\langle \phi_a | \hat{\Pi}(0,\mathbf{x}) | \zeta \rangle = - i \dfrac{\delta}{\delta \phi_a(\mathbf{x})} \langle \phi_a | \zeta \rangle$ for any state $|\zeta \rangle$ in the Fock space.
So far what I have shown that the third commutation relation implies that $[ \hat{\Phi}(x^0,\mathbf{x}), \hat{\Pi}(x^0,\mathbf{y})^n ] = i n \delta^{(3)}(\mathbf{x} - \mathbf{y})\hat{\Pi}(x^0,\mathbf{y})^{n-1}$ for any $n \geq 1$. From this it follows that for any number $\epsilon$ we have $[ \hat{\Phi}(x^0,\mathbf{x}), e^{ - i \epsilon \hat{\Pi}(x^0,\mathbf{y}) } ] = \epsilon e^{ - i \epsilon \hat{\Pi}(x^0,\mathbf{y}) } \delta^{(3)}(\mathbf{x} - \mathbf{y})$.
From here, applying this commutator to a field eigenstate $|\phi_a \rangle$ yields $$ \hat{\Pi}(x^0, \mathbf{x}) \big( e^{ - i \epsilon \hat{\Phi}(x^0,\mathbf{y}) } | \phi_a \rangle \big) = \big( \phi_a(\mathbf{x}) + \epsilon \delta^{(3)}(\mathbf{x} - \mathbf{y} ) \big) \big( e^{ - i \epsilon \hat{\Phi}(x^0,\mathbf{y}) } | \phi_a \rangle \big) $$ This is where I get stuck though. I was hoping to use the above to show that $\big( e^{ - i \epsilon \hat{\Phi}(x^0,\mathbf{y}) } | \phi_a \rangle \big) \propto | \phi_a + \epsilon \rangle$, but the extra $\delta$-function doesn't seem to make this work (even without the $\delta$-function problem, this would still just be a proportionality, up to some phase). From there my idea was to consider the inner product $\langle \phi_a | e^{- i \epsilon \hat{\Pi}(0,\mathbf{x}) } | \zeta \rangle$ and take the limit $\epsilon$ to prove the result.