Elastic Proton-Proton Scattering at tree level

Hello im trying compute the Cross section of an elastic proton-proton scattering and my results are beyond terrrible and i can't figure out where my mistakes lie. I'm going to describe what i tried to do and hopefully someone can help me figure out where i mess up, im not asking you to do these calculations on your own but just see if you can see an obvious mistake in my reasoning or results. Im also aware that the problem could lie entirely on my computer code but i checked it over and over again and can't spot a mistake, so im thinking my error lies in the Physics. I also did not take into account the electromagnetic interaction because it should be negligible compared to the nuclear force. Thank you for reading this in advance

I use the following Lagrangians to Describe the Interactions of the protons with the exchange mesons:

$$$$\begin{array}{c} \mathcal{L}_{\mathrm{NNps}}=-g_{p s} \bar{\psi} i \gamma_{5} \psi \varphi_{\mathrm{ps}} \\ \mathcal{L}_{\mathrm{NNs}}=-g_{s} \bar{\psi} \psi \varphi_{\mathrm{s}} \\ \mathcal{L}_{\mathrm{NNv}}=-g_{v} \bar{\psi} \gamma_{\mu} \psi \varphi_{\mathrm{v}}^{\mu}-\frac{f_{\mathrm{v}}}{4 M_{p}} \bar{\psi} \sigma_{\mu \nu} \psi \cdot\left(\partial^{\mu} \varphi_{\mathrm{v}}^{\nu}-\partial^{\nu} \varphi_{\mathrm{v}}^{\mu}\right) \end{array}$$$$

Each proton described by the Dirac Spinor $$$$u\left(\boldsymbol{q}, \lambda\right)=\sqrt{E+m}\left(\begin{array}{c} 1 \\ \frac{2 \lambda q}{E+m} \end{array}\right)\left|\lambda\right\rangle$$$$

What i did is to calculate the Matrix Elements for each of these interaction in the C.O.M frame like this (i will only include the calculations on one vertex of the feynman diagrams)

$$N=E+m$$ for all the calculations that follow

For the Pseudoscalar mesons: $$$$ig_{ps}N\left(\frac{2\lambda q}{N}-\frac{2\lambda' q}{N}\right) \left\langle\lambda'\mid \lambda\right\rangle$$$$

For the Scalar mesons: $$$$g_{s}N\left(1-\frac{4\lambda\lambda'q^2}{N^2}\right) \left\langle\lambda'\mid \lambda\right\rangle$$$$

For the Vector Mesons: $$$$j^{0}=(f_\mathrm{v}+g_\mathrm{v})N\left(1+\frac{4\lambda\lambda'q^2}{N^2}\right)\left\langle\lambda'\mid \lambda\right\rangle-\frac{f_\mathrm{v}}{2m}(q'^0+q^0)\left(1-\frac{4\lambda\lambda'q^2}{N^2}\right)\left\langle\lambda'\mid \lambda\right\rangle$$$$ $$$$j^{k}=(f_\mathrm{v}+g_\mathrm{v})N\left(\frac{2\lambda q}{N}+\frac{2\lambda' q}{N}\right)\left\langle\lambda'\mid\sigma^k\mid\lambda\right\rangle-\frac{f_\mathrm{v}}{2m}(q'^k+q^k)\left(1-\frac{4\lambda\lambda'q^2}{N^2}\right)\left\langle\lambda'\mid \lambda\right\rangle$$$$

What i did next is to find general expressions for $$$$\left\langle\lambda'\mid \lambda\right\rangle$$$$ and $$$$\left\langle\lambda'\mid \sigma^k \mid \lambda\right\rangle$$$$ So my PC could then sum over all posible helicity combinations. I added all 16 different matrix Elements like this: \begin{aligned} \left\langle\left|\mathcal{M}_{f i}\right|^{2}\right\rangle &=\frac{1}{4}\left(\left|\mathcal{M}_{R R}\right|^{2}+\left|\mathcal{M}_{R L}\right|^{2}+\left|\mathcal{M}_{L R}\right|^{2}+\left|\mathcal{M}_{L L}\right|^{2}\right) \\ &=\frac{1}{4}\left(\left|\mathcal{M}_{R R \rightarrow R R}\right|^{2}+\left|\mathcal{M}_{R R \rightarrow R L}\right|^{2}+\cdots+\left|\mathcal{M}_{R L \rightarrow R R}\right|^{2}+\cdots\right) \end{aligned} Then i took into account both Feynman Diagrams (Picture from David Tong Lectures on Quantum Field Theory)

So the Matrix Element is $$$$\mathcal{M}=\mathcal{M}_1-\mathcal{M}_2$$$$ Where the Matrix Elements Correspond to each Feynman diagrams and the minus sign is used because of Fermi-Dirac Statistics. For each of the Feynman diagrams i summed for the six exchange Bosons listed below with the following properties and i used their respective propagator in computing the Matrix Elements: $$$$\begin{array}{lllllllc} \hline \hline \text { Meson } & J^{P} & I^{G} & g_{\alpha}^{2} / 4 \pi & f_{\alpha} / g_{\alpha} & m_{\alpha}(\mathrm{MeV}) & \Lambda_{\alpha}(\mathrm{GeV}) & 2 n_{\alpha} \\ \hline \pi & 0^{-} & 1^{-} & 14.4 & & 138.03 & 1.7 & 2 \\ \eta & 0^{-} & 0^{+} & 3 & & 548.8 & 1.5 & 2 \\ \rho & 1^{-} & 1^{+} & 0.9 & 6.1 & 769 & 1.85 & 4 \\ \omega & 1^{-} & 0^{-} & 24.5 & 0 & 782.6 & 1.85 & 4 \\ \delta & 0^{+} & 1^{-} & 2.488 & & 983 & 2.0 & 2 \\ \sigma^{a} & 0^{+} & 0^{+} & 8.9437 & & 550 & 1.9 & 2 \\ & & & (18.3773) & & (720) & (2.0) & (2) \\ \hline \end{array}$$$$ At each vertex is used a form factor like this: $$$$\mathcal{F}_{\alpha}\left[\left(\boldsymbol{q}^{\prime}-\boldsymbol{q}\right)^{2}\right]=\left[\frac{\Lambda_{\alpha}^{2}-m_{\alpha}^{2}}{\Lambda_{\alpha}^{2}+\left(\boldsymbol{q}^{\prime}-\boldsymbol{q}\right)^{2}}\right]^{n_{\alpha}}$$$$

Lastly i use the following equation to compute the Cross section $$$$\sigma=\frac{1}{64 \pi^{2} s} \frac{\mathrm{p}_{f}}{\mathrm{p}_{i}} \int\left|\mathcal{M}_{f i}\right|^{2} \mathrm{~d} \Omega$$$$

Where the cross section must include a factor of $$1/2$$ because each final state may only be counted once in integrating over the angles

To my mind all of this should yield reasonable results. Not exact since this is a tree level treatment of the scattering but my results as i mentioned before do not even come close to the experimental results. I know that for Moller Scattering (Elastic electron scattering) tree level treatment is not enough since the integral is diverging. Is this one case where tree level treatment also doesn't work? Again i thank you for reading all this. The Graph below is my attempt where you can see my results where the Green points are the experimental findings and the blue curve is the cross section i calculated.

• What do you mean by "does not come close"? 10%, factor 2, factor 10? What energy range? Tree level is not exact, of course, even at threshold. However, it should not be way off (i.e. you may miss by factors of 2, but not factors of 10). – Thomas May 21 at 15:32
• Im gonna add my curve so you can see what i mean, you can see thatmy curve is way off. Now i know that i shouldnt except exact results but shouldn't their at least be some resemblance? – Μπαμπης Ποζουκιδης May 21 at 17:25
• You are looking at fairly high energies, of order 10 GeV. In this regime it is not the case that mesons are the right degrees of freedom. The couplings and form-factors you use were presumably extracted from elastic scattering below the inelatstic threshold, $\sqrt{s}<(2m_N+m_\pi)$. – Thomas May 21 at 18:51
• Make that your answer so I can award you the bounty – Μπαμπης Ποζουκιδης May 22 at 7:51

At low energy, using $$\sqrt{s}$$ is a little misleading, because $$s$$ is dominated by the rest masses. More relevant is the relative three-momentum in the CM frame. If that exceed a few hundred MeV, then relativistic effects and in-elastic channels become important. Also, for momenta bigger than a GeV or so you start to resolve the quark-gluon substructure, and a meson-exchange description is not appropriate.
You may think that at high energies, $$\sqrt{s}>10$$ GeV or so, perturbative QCD is reliable. Unfortunately, this is not the case. Both total and elastic nucleon-nucleon cross sections are dominated by small angle scattering, and pQCD is not applicable. There are semi-phenomenological models based on pomerons and Regge trajectories, but as theories of total cross sections these are not grounded in QCD.