To approximate this Gaussian potential in both the radial and axial directions of the beam, the intensity profile must be expanded to second order in z and r for r = 0 and z = 0 respectively and equated to the harmonic potential 1 2 m ( ω ω z 2 z 2 + ω ω r 2 r 2 ) . These expansions are evaluated assuming fixed power.