Samsung Advanced Institute of Technology

Fast and accurate electromagnetic simulation of large-area metasurfaces remains a major ob-
stacle towards automating their design. In this paper, we propose a metasurface simulation dis-
tribution strategy which achieves a linear reduc-tion in the simulation time with the number
of compute node. Combining this distribution strategy with a GPU-based implementation of
the Transition-matrix method, we perform accu-rate simulations and adjoint sensitivity analysis of
large-area metasurfaces. We demonstrate ability to perform a distributed simulation of large-area
(over 600lambda x 600lambda), while accurately accounting for scatterer-scatterer interactions significantly be-
yond the locally periodic approximation.