We consider the extrinsic geometry of surfaces in simply isotropic space, a three-dimensional space equipped with a rank 2 metric of index zero. Since the metric is degenerate, a surface normal cannot be unequivocally defined based on metric properties only. To understand the contrast between distinct choices of an isotropic Gauss map, here we study surfaces with a Gauss map whose coordinates are eigenfunctions of the surface Laplace-Beltrami operator. We take into account two choices, the so-called minimal and parabolic normals, and show that when applied to simply isotropic invariant surfaces the condition that the coordinates of the corresponding Gauss map are eigenfunctions leads to planes, certain cylinders, or surfaces with constant isotropic mean curvature. Finally, we also investigate (non-necessarily invariant) surfaces with harmonic Gauss map and show this characterizes constant mean curvature surfaces.