“…In our case, in addition, we deal with higher order terms, such as the third order ones (r 2 z and z 3 ones). The theory is more complex and generally no more linked to the classical world with for instance negative values for the Wigner function [82,86]. Fortunately, if the extra terms are considered as a perturbation L 1 on the Lagrangian L = L 0 + L 1 , the phase shift δφ introduced into the final wavefunction, by the perturbation L 1 , is given simply (to the first order) by the integral of the perturbation along the classical unperturbed path Γ 0 [87][88][89]:…”