We present an accurate method to calculate the Gouy Phase Shift of arbitrary beams, combining the Fresnel-Kirchhoff diffraction formula, or, alternatively, the Angular Spectrum Representation with the transverse momenta spreading formulation. Our approach was assessed by calculating the Gouy phase as a function of the distance to the beam waist of various well-known beams, and the results are compared with the expected analytical values. Good agreement is observed for optical beams that satisfy the paraxial approximation. Afterwards, the Gouy phase shifts of optical beams diffracted from different waveguides end faces are also calculated, and the results are discussed in this manuscript.
Daniel et al. (Thu,) studied this question.