sqrt(1 - (proton mass * light speed^2)^2 / (2000 J)^2)