Actually,the exact formulas for Ex,Ey,Ez are the same as separable Helmholtz eqations in Cartesian coordinates. The Gaussian packet(beam) form is just a kind of B.C. You can simply superpose plane waves(in coordinate) to get a Gaussian shape. In this sense, kx,ky are generally non-zero,so Ez...