We follow the procedure proposed by Ref. to derive the correct form of optical force density for two-dimensional metamaterial systems. We use a mesoscopic lattice model as a numerically accurate reference, and compare its results with those obtained from an macroscopic effective media model using two different stress tensors. We show that only the Helmholtz stress tensor is able to give the just answer.
The paper is available at here.