The calculation of the permittivity of a composite via the FEM requires knowledge of the permittivity ε(r) at all points in the RVE. The bulk matrix response can be calculated via quantum dynamics or obtained by experiment. However, the response of the inclusions and that of the interface layers needs to be determined: the former may differ significantly from either bulk metal or nanoparticle-in-vacuum values, whereas the latter is completely unknown.