Complex modes and effective refractive index in 3D periodic arrays of plasmonic nanospheres

We characterize the modes with complex wavenumber for both longitudinal and transverse polarization states (with respect to the mode traveling direction) in three dimensional (3D) periodic arrays of plasmonic nanospheres, including metal losses. The Ewald representation of the required dyadic periodic Green's function to represent the field in 3D periodic arrays is derived from the scalar case, which can be analytically continued into the complex wavenumber space. We observe the presence of one longitudinal mode and two transverse modes, one forward and one backward. Despite the presence of two modes for transverse polarization, we notice that the forward one is "dominant" (i.e., it contributes most to the field in the array). Therefore, in case of transverse polarization, we describe the composite material in terms of a homogenized effective refractive index, comparing results from (i) modal analysis, (ii) Maxwell Garnett theory, (iii) Nicolson-Ross-Weir retrieval method from scattering parameters for finite thickness structures (considering different thicknesses, showing consistency of results), and (iv) the fitting of the fields obtained through HFSS simulations. The agreement among the different methods justifies the performed homogenization procedure in case of transverse polarization.