In this paper, we investigate numerically the characteristics of surfaceplasmon polaritons (SPPs) sustained by two-dimensional arrays of metallic pillars protruding out of planar metal surfaces at terahertz (THz) frequencies. Various shapes of the pillars
