  Surface plasmons have been used recently to generate heat nanosources, the
intensity of which can be tuned, for example, with the wavelength of the
excitation radiation. In this paper, we present versatile analytical and
numerical investigations for the three-dimensional computation of the
temperature rise in complex planar arrays of metallic particles. In the
particular case of elongated particles sustaining transverse and longitudinal
plasmon modes, we show a simple temperature rise control of the surrounding
medium when turning the incident polarization. This formalism is then used for
designing novel thermoplasmonic metasurfaces for the nanoscale remote control
of heat flux and temperature gradients.
