  An analytical model for homogenization of Metamaterials (MM) with interacting
Metaatoms (MA) is developed based on multipole approach for bulk media. The
interaction is assumed to be near field type. i.e. no retardation in lateral
direction between adjacent MAs is assumed. The interaction takes place between
the adjacent MAs only in lateral (perpendicular to the propagation) direction;
the considered MM is supposed to consist of non-interacting identical layers
with periodically spaced MAs. It is shown, that the interaction between MAs
leads to significant changes of dispersion characteristics, in particular,
appearance of the spatial dispersion is emphasized. The results indicate
increase of the available real part of the k vector values and increase of
imaginary part by the approximately same percentage, which leads to the final
conclusion that the effect of coupling does not provide extra opportunities for
the resolution enhancement.
