“…Few years later, exploiting the QFT advances [22,23], these shifts were calculated up to the fifth order in the fine structure constant α [24,25]. However, in those calculations and in many of the others that followed, the final results were obtained in the semiclassical approximation: a crucial role was played by the value at the origin of the non relativistic Coulomb wave function, raising a long lasting question on the opportunity of adopting some kind of smearing of the corresponding δ-function [26,27]. In parallel with the investigations of the relativistic few body systems, the interest increased also for the analysis of their non-quantum counterpart, mainly based on the group symmetries and on the constraint theory [28][29][30]: the hope was to be able to catch the essence of the relationship between canonicity and relativistic covariance, in order to have sound guidelines to the quantization procedure.…”