In Landau Classical theory of fields, in the section of elastic collision(chapter 2, section 13), he mentions that for $m_1>m_2$, the scattering angle $\theta_i$ cannot exceed a maximum value, and by elementary computations, the value is
$$\sin\theta=\frac{m_2}{m_1}$$
The $4$ momentums are:
$$p_1=\left(E_1,p_1,0,0\right)$$ $$p_2=\left(m_2,0,0,0\right)$$ $$p’_1=\left(E’_1,p_1-p_x,-p_y,0\right)$$ $$p’_2=\left(E’_2,p_x,p_y,0\right)$$
From conservation of 4-momentum,
$$E_1+m_2=E’_1+E’_2$$
And from the magnitude of 4-momentum,
$$m_1^2+p_1^2=E_1^2$$
$$m_1^2+p_1^2+p_x^2+p_y^2={E’}_1^2+2p_1p_x$$
$$m_2^2+p_x^2+p_y^2={E’}_2^2$$
By considering$\vec{p}_1\cdot\vec{p}’_1$, i get $$p_1^2-p_1p_x=p_1\sqrt{p_1^2+p_x^2-2p_1p_x+p_y^2}\cos\theta$$
After this I’m not exactly sure how to maximise $\theta$, maximising $\theta$ is equivalent to maximising $p_y$, but not exactly sure what’s bounding it from above. How should I go about finding the maxima of $\theta$ with these relations?
I’ve considered if $p_x=0$, but that isn’t possible(conservation of energy)
Update 1:
By boosting into the center of mass frame(C-frame) and back to the lab frame(L-frame), i got something similar:
CM velocity: $v=\frac{-p_1}{E_1+m_2}$
Initial L-frame 4-momentum: $$p_{1L}=\left(E_1,p_1,0,0\right)$$ $$p_{2L}=\left(m_2,0,0,0\right)$$
Initial C-frame 4-momentum: $$p_{1C}=\gamma\left(E_1+vp_1,p_1+vE_1,0,0\right)$$ $$p_{2C}=\gamma\left(m_2,vm_2,0,0\right)$$
Final C-frame 4-momentum(for a angle of $\frac{\pi}{2}$, largest possible angle): $$p'_{1C}=\gamma\left(E_1+vp_1,0,-vm_2,0\right)$$ $$p'_{2C}=\gamma\left(m_2,0,vm_2,0\right)$$
Final L-frame 4-momentum $$p'_{1L}=\gamma\left(\gamma\left(E_1+vp_1\right),-\gamma \left(vE_1+v^2p_1\right),-vm_2,0\right)$$ $$p'_{2L}\gamma\left(\gamma m_2,-\gamma vm_2,vm_2,0\right)$$
From the magnitude of $p'_1$ in the C-frame, we get $$\left(\gamma E_1+\gamma vp_1\right)^2=m_1^2+\gamma^2v^2m_2^2$$ Now using $p'_1$ in the L-frame, we find $\sin^2\theta$ $$\sin^2\theta=\frac{p^{'2}_{1y}}{p^{'2}_{1x}+p^{'2}_{1y}}=\frac{\gamma^2v^2m_2^2}{\gamma^4v^2\left(E_1+vp_1\right)^2+\gamma^2v^2m_2^2}=\frac{m_2^2}{\gamma^2\left(E_1+vp_1\right)^2+m_2^2}$$ $$=\frac{m_2^2}{m_1^2+\gamma^2v^2m_2^2+m_2^2}=\frac{m_2^2}{m_1^2+\gamma^2m_2^2}$$
Not sure if I messed up somewhere and got an extra $\gamma^2m_2^2$ term