The surface integral formulation is a flexible, multiscale and accurate tool to simulate light scattering on nanostructures. Its generalization to periodic arrays is introduced in this paper. The general electromagnetic scattering problem is reduced to a discretizated model using the Method of Moments on the surface of the scatterers in the unit cell. The study of the resonances of an array of bowtie antennas illustrates the main features of the method. When placed into an array, the bowtie antennas show additional resonances compared to those of an individual antenna. Using the surface integral formulation, we are able to investigate both near-field and far-field properties of these resonances, with a high level of accuracy. (C) 2010 Elsevier B.V. All rights reserved.