The detailed electric field distribution is necessary for calculating the trapping potential above the waveguide. An approximation for the y-dependence of the fields is given by the analytically-known solution for the slab waveguide; this gives exponential decay, as already discussed. To get knowledge of the actual potential shape in the x-y plane, we used a full-vector finite element calculation.