I'd like to prove the LSZ formula, but there is a specific step that is bugging me a lot. I know there are many subtleties in its derivation, but I'm not worrying about this right now: I'm trying to understand the naive proof, so to speak.
You can find an example of the usual proof here: http://isites.harvard.edu/fs/docs/icb.topic473482.files/06-LSZ.pdf
My question is: how to get eqs. (24-25) from eq. (22-23). In (22-23) the time-ordering symbol is to the left of $(\partial^2+m^2)\phi$ and in (24-25) it is to the right of the differential operator.
I'm asking how to go from $T(\partial^2+m^2)\phi\cdots$ to $(\partial^2+m^2) T\phi\cdots$. I feel that this cannot be done in general, because the symbol $T$ will introduce some $\Theta(x_0-y_0)$ functions, which when differentiated will give rise to some deltas.
This step makes no sense to me... Is my question legitimate? Is there anything from the proof I'm missing?
If I think of specific examples, I find different results depending on whether the $T$ symbol is to the right or to the left of the Klein-Gordon differential operator. This means that the order is important... so, which is the right order? Should the $T$ symbol be placed to the right or to the left of the KG operator? In the begining of the proof the $T$ symbol is always to the left, and in the end, it is always to the right.
This same problem appears on many proofs online, such as
http://isites.harvard.edu/fs/docs/icb.topic473482.files/06-LSZ.pdf (eqs. 22-25)
http://www2.ph.ed.ac.uk/~egardi/MQFT_2013/MQFT_2013_lecture_2.pdf (eqs. 32-33)
Srednicki's book, page 51 (eqs. 5.14-5.15) online: http://web.physics.ucsb.edu/~mark/qft.html
etc.
Edit: I was asked to post a self-contained question, so I'll write the delails here:
Take the in state to be $|i\rangle\propto a^\dagger_1(-\infty)a^\dagger_2(-\infty)|0\rangle$ and the out state to be $|f\rangle\propto a^\dagger_3(+\infty)a^\dagger_4(+\infty)|0\rangle$. Then the transition amplitude is $$ \langle i|f\rangle \propto \langle 0|a_2(-\infty)a_1(-\infty)a^\dagger_3(+\infty)a^\dagger_4(+\infty)|0\rangle= $$ $$ =\langle 0|Ta_2(-\infty)a_1(-\infty)a^\dagger_3(+\infty)a^\dagger_4(+\infty)|0\rangle $$
Now write $a(+\infty)-a(-\infty)\propto \int \mathrm d x\ \mathrm e^{\cdots}(\partial^2+m^2)\phi$; all terms with $a$'s or $a^\dagger$'s annihilate the vacuum, so that the only remaining term is $$ \langle i|f\rangle \propto \langle 0|T \int \mathrm dx_1\mathrm dx_2\mathrm dx_3 \mathrm dx_4\ \mathrm e^{\cdots}(\partial_1^2+m^2)(\partial_2^2+m^2)(\partial_3^2+m^2)(\partial_4^2+m^2) \phi_1\phi_2\phi_3\phi_4 |0\rangle $$
This is usually writen as $$ \langle i|f\rangle \propto \int \mathrm dx_1\mathrm dx_2\mathrm dx_3 \mathrm dx_4\ \mathrm e^{\cdots}(\partial_1^2+m^2)(\partial_2^2+m^2)(\partial_3^2+m^2)(\partial_4^2+m^2) \langle 0|T\phi_1\phi_2\phi_3\phi_4 |0\rangle $$ but this step is problematic, because the $T$ symbol and the KG operator don't commute, right?