As for the Heisenberg picture... it's still there! You can switch to a Heisenberg picture through a unitary transformation. In this picture the operators evolve through time as:
[tex]A_H(t) = e^{iHt} A_H(0) e^{-iHt}[/tex]
The operator [tex]e^{-iHt)[/tex] is the time-evolution operator. In the Schroedinger picture the states transform in the usual way:
[tex]|\psi_S(t)\rangle = e^{-iHt}|\psi_S(t)\rangle[/tex]
And here is something you do not hear every day (but is certainly still true!). The equation of motion of the time evolution operator is nothing but Schroedingers equation!
[tex]\frac{d}{dt} e^{-iHt} = -i H e^{-iHt}[/tex]
Yes, the Schroedinger equation still plays a role in relativistic quantum field theory!
In fact, for fields operators we can apply the time evolution for operators. Take for instance the conjugate field [tex]\pi(x,t)[/tex]
[tex]\pi(x,t) = e^{iHt} \pi(x,0) e^{-iHt}[/tex]
which -- when taking the time derivative -- reduces to Heisenbergs equations of motion:
[tex]\partial_t \phi(x,t) = i[H,\pi(x,t)][/tex]
But for a free scalar field theory, i.e. KG theory, this commutator can be computed explicitly. Namely, the Hamiltonian is given by:
[tex]H = \frac{1}{2} \int d^3 x' \left(\pi(x,0)^2 + (\nabla\phi(x,t))^2 +m^2\phi(x,t)^2\right)[/tex]
Now, the equal-time commutation relations should be familiar
[tex][\phi(x,t),\pi(x',t)] = i \delta(x-x')[/tex]
with on the right hand side the three dimensional spatial delta function. All other commutators are zero. The commutator we are after is then...
[tex][H,\pi(x,t)] = i\int d^3x' \left(\nabla^2 \phi(x',t) + m^2\phi(x',t)\right)\delta(x-x') = \nabla^2 \phi(x',t) + m^2\phi(x',t)[/tex]
You need to apply some partial integration to obtain this, but it's quite straightforward. Anyways, now, this commutator determines the time-evolution of the conjugate field. So we have:
[tex]\partial_t \pi(x,t) = \nabla^2 \phi(x,t) - m^2\phi(x,t)[/tex]
This is simply the Heisenberg equation of motion for the conjugate field. But remember, the conjugate field [tex]\pi(x,t)[/tex] is nothing but the time derivative of the original KG field [tex]\phi(x,t)[/tex]. And so we end up with...
[tex]\partial_t^2 \phi(x,t) - \nabla^2 \phi(x,t) + m^2\phi(x,t) = 0[/tex]
Ha, that should look familiar!