Propagation of Gaussian beams in free space

the computational domain
n = 1.0, w = 15.1 µm, vacuum wavelength 1.0 µm,
TE polarization, Mx = Mz = 150.

QUEP simulations

Ey(x,z,t) Ey(x,z,t) Ey(x,z,t) Ey(x,z,t) Ey(x,z,t) Ey(x,z,t) Ey(x,z,t) Ey(x,z,t) Ey(x,z,t) Ey(x,z,t)