Changed the definition of the external potential: U = -Z*alpha/r