I am following the derivation steps of the LSZ reduction formula as done in the lecture notes by Timo Weigand the section on S matrix in the full interacting theory, here:
http://www.thphys.uni-heidelberg.de/~weigand/QFT1-13-14/SkriptQFT1.pdf
On page 50 the equation 2.56 below does not make sense to me.
$$ \langle p_1,\dots,p_n, out| a^{\dagger}_{out}(q_1)|q_2,\dots, q_r\rangle \sqrt(2E_{q_1})=$$ $$=\sum_{k=1}^{n}2E_{p_k}(2\pi)^3\delta^{3}(p_k-q_1)\langle p_1,\dots,\hat{p_k},\dots, p_n, out|q_2,\dots,q_r, in\rangle ,\tag{2.56}$$
where $\hat{p_k}$ has to be taken out.
For instance when I expand for n=2 the right side by replacing the delta function with commutation involving annihilation and creation of "out" operators I do not obtain the left side of the equation.
I mean using $$(2\pi)^3\delta^{3}(p_k-q_1)=[a(p_k), a^{\dagger}(q_1)], \qquad \langle p_1,p_2|/\sqrt(2E_{p2})=\langle p_1|a_{out}(p_2),$$
but this does not show that the two sides are equal.
Also do not see why this equation describes a process where one of the in- and outgoing states are identical and do not participate in scattering.