 \varphi_{j} - \left %(\omega_{0}^2/\omega\right ) \sin \varphi_{j}, 
 \end{aligned}$$ and the solution of Eq. (\[maxwell-laplace\]) reads $$\begin{aligned}
\label{electric}
{\bf E}({\bf r},\omega ) = \frac{4\pi d_{0}q \omega^{2}}{c^2} \sum_{j}{\bf G}({\bf r},{\bf r}_{j},\omega) \cdot {\bf e}_{j} v_{j}, \end{aligned}$$ where ${\bf G}({\bf r},{\bf r}',\omega)$ is the electric field Green diadic *in the presence of* NP. From Eq. (\[displ\]), for photon frequency close to dipoles frequency, $\omega \approx \omega_{0}$, we obtain a coupled system of equations for normalized displacements, $$\begin{aligned}
\label{system}
\sum_{k} \Bigl[(\omega_{0}-\omega ) \delta_{jk} +  \Sigma _{jk}\Bigr] v_{k} = 
%\frac{v_{0j}}{2\omega_{0}}=
\frac{-i}{2} e^{-i\varphi_{j}},\end{aligned}$$ where $\Sigma _{jk}=\Delta_{jk}-\frac{i}{2}\Gamma_{jk}$ is the complex *self-energy matrix*, given by $$\begin{aligned}
\label{self-general}
\Sigma _{jk}(\omega)= - \frac{2\pi q^2 \omega_{0}}{m c^2} \, {\bf e}_{j}\cdot  {\bf G}({\bf r}_{j},{\bf r}_{k};\omega) \cdot {\bf e}_{k}.\end{aligned}$$ The system Eq. (\[system\]) determines eigenstates of the ensemble of $N$ dipoles in common radiation and plasmon field, while complex eigenvalues of the self-energy matrix give eigenstates frequency shifts with respect to $\omega_{0}$ and their decay rates. In the absence of NP, real and imaginary parts of the self-energy matrix are dipole-dipole interaction and radiation coupling between dipoles $j$ and $k$, given by (in lowest order in $kr_{jk}$) $$\begin{aligned}
\label{dipole-dipole}
&\Delta_{jk}^{0}=\frac{3\Gamma_{0
\varphi_{j } - \left% (\omega_{0}^2/\omega\right) \sin \varphi_{j }, 
  \end{aligned}$$ and the solution of Eq.   (\[maxwell - laplace\ ]) reads $ $ \begin{aligned } 
 \label{electric } 
 { \bf E}({\bf r},\omega) = \frac{4\pi d_{0}q \omega^{2}}{c^2 } \sum_{j}{\bf G}({\bf r},{\bf r}_{j},\omega) \cdot { \bf e}_{j } v_{j }, \end{aligned}$$ where $ { \bf G}({\bf r},{\bf r}',\omega)$ is the electric field Green diadic * in the presence of * NP. From Eq.   (\[displ\ ]), for photon frequency close to dipole frequency, $ \omega \approx \omega_{0}$, we receive a coupled system of equations for normalize displacements, $ $ \begin{aligned } 
 \label{system } 
 \sum_{k } \Bigl[(\omega_{0}-\omega) \delta_{jk } +   \Sigma _ { jk}\Bigr ] v_{k } = 
% \frac{v_{0j}}{2\omega_{0}}= 
 \frac{-i}{2 } e^{-i\varphi_{j}},\end{aligned}$$ where $ \Sigma _ { jk}=\Delta_{jk}-\frac{i}{2}\Gamma_{jk}$ is the complex * self - energy matrix *, give by $ $ \begin{aligned } 
 \label{self - general } 
 \Sigma _ { jk}(\omega)= - \frac{2\pi q^2 \omega_{0}}{m c^2 } \, { \bf e}_{j}\cdot   { \bf G}({\bf r}_{j},{\bf r}_{k};\omega) \cdot { \bf e}_{k}.\end{aligned}$$ The system Eq.   (\[system\ ]) determine eigenstates of the ensemble of $ N$ dipoles in coarse radiation and plasmon field, while complex eigenvalue of the self - energy matrix give eigenstates frequency shifts with respect to $ \omega_{0}$ and their decay rates. In the absence of NP, real and imaginary part of the self - energy matrix are dipole - dipole interaction and radiation coupling between dipole $ j$ and $ k$, given by (in lowest order in $ kr_{jk}$) $ $ \begin{aligned } 
 \label{dipole - dipole } 
 & \Delta_{jk}^{0}=\frac{3\Gamma_{0
 \vagphi_{j} - \left %(\omega_{0}^2/\omega\rinht ) \sin \varphi_{j}, 
 \end{almgned}$$ ahd the sulution of Eq. (\[maxwell-laplace\]) ceadw $$\begun{aligned}
\label{electric}
{\cf E}({\bf r},\olega ) = \feac{4\pm d_{0}q \omega^{2}}{c^2} \sum_{j}{\ur G}({\bf r},{\ny r}_{j},\ojcga) \cbov {\bf e}_{j} v_{j}, \end{akigned}$$ whese ${\bf G}({\bf r},{\bf s}',\ooeya)$ is the electric field Green diadis *in thr oresence of* NP. Frok Eq. (\[djspl\]), for photon frequency close to dipolev frequency, $\okega \approx \omega_{0}$, we obtaij a foupled system of fquations fir njemalized disolacements, $$\begin{alignes}
\label{system}
\sum_{k} \Bigl[(\omega_{0}-\omega ) \delca_{jk} +  \Sigmq _{hk}\Blcr] v_{k} = 
%\frac{t_{0j}}{2\omegw_{0}}=
\frac{-i}{2} e^{-i\varini_{j}},\end{dligned}$$ where $\Sigma _{jl}=\Denta_{hk}-\frac{i}{2}\Gamma_{jk}$ is the complex *self-energi matrix*, ghvzn by $$\begin{aligned}
\labwl{welf-ggneran}
\Sigoq _{jy}(\omtga)= - \rrac{2\pi q^2 \kmega_{0}}{m c^2} \, {\gf e}_{j}\cdot  {\vf G}({\bf r}_{j},{\bf r}_{k};\omega) \cqit {\bf e}_{k}.\end{aljgned}$$ Ehq system Eq. (\[system\]) determines eigenstatev or the ensemble of $N$ dipiles in common radiatlon and pjasmon field, while complex eigenvalues of the selx-enerjy maurlw gixw figenstates frequency shifts with respect to $\jjeba_{0}$ and their deccy rates. In the sbdemse of NP, real and imagjnary parts of the self-engrgy mqtrix are dipple-dipole interaction and rqdiation couilint between dipoles $l$ and $k$, giveu by (im lowrst order in $kr_{jk}$) $$\begin{clignes}
\label{dipolf-dipole}
&\Demga_{jk}^{0}=\frac{3\Gamma_{0
\varphi_{j} - \left %(\omega_{0}^2/\omega\right ) \sin \varphi_{j}, the of Eq. reads $$\begin{aligned} \label{electric} \frac{4\pi \omega^{2}}{c^2} \sum_{j}{\bf G}({\bf r}_{j},\omega) \cdot {\bf v_{j}, \end{aligned}$$ where ${\bf G}({\bf r},{\bf is the electric field Green diadic *in the presence of* NP. From Eq. for photon frequency close to dipoles frequency, $\omega \approx \omega_{0}$, we obtain a system equations normalized $$\begin{aligned} \label{system} \sum_{k} \Bigl[(\omega_{0}-\omega ) \delta_{jk} + \Sigma _{jk}\Bigr] v_{k} = %\frac{v_{0j}}{2\omega_{0}}= \frac{-i}{2} e^{-i\varphi_{j}},\end{aligned}$$ where $\Sigma is the complex *self-energy matrix*, given by $$\begin{aligned} \Sigma _{jk}(\omega)= - \frac{2\pi \omega_{0}}{m c^2} \, {\bf e}_{j}\cdot G}({\bf r}_{k};\omega) \cdot e}_{k}.\end{aligned}$$ system (\[system\]) determines eigenstates the ensemble of $N$ dipoles in common radiation and plasmon field, while complex eigenvalues of the self-energy give eigenstates with respect $\omega_{0}$ their rates. In the NP, real and imaginary parts of are dipole-dipole interaction and radiation coupling between dipoles and $k$, by (in lowest order in $kr_{jk}$) \label{dipole-dipole} &\Delta_{jk}^{0}=\frac{3\Gamma_{0
 \varphi_{j} - \left %(\omega_{0}^2/\omega\righT ) \sin \varphi_{J}, 
 \end{aLigNed}$$ AnD the SoluTion of Eq. (\[maxwelL-LaplAce\]) reads $$\begin{aligned}
\laBel{elEcTRic}
{\bF e}({\bF r},\omeGa ) = \frac{4\pI D_{0}q \OMEga^{2}}{C^2} \sUm_{J}{\bf g}({\bF R},{\bF r}_{j},\omEga) \Cdot {\bf e}_{J} v_{j}, \end{aligNed}$$ WhEre ${\bf G}({\bf r},{\bf r}',\OMeGa)$ is the eleCtrIc field Green DiaDic *in tHe PreSEnce oF* NP. from EQ. (\[displ\]), FOr photOn frequenCy CLose to DIpoles fREQuEncy, $\Omega \approx \omega_{0}$, wE ObTAin a coupled sysTem of eQuATiONS foR noRmalized diSpLacemENts, $$\begiN{AlIGNEd}
\lABel{system}
\sum_{k} \bigl[(\omega_{0}-\omEGa ) \dElta_{jk} +  \siGma _{JK}\Bigr] v_{K} = 
%\frac{V_{0j}}{2\OMegA_{0}}=
\frac{-i}{2} e^{-i\varPhi_{j}},\End{aligneD}$$ where $\sIgma _{jk}=\DELta_{jk}-\frAc{i}{2}\GamMa_{jK}$ is The cOMpLeX *seLf-ENerGY mAtrIX*, giVen by $$\begIn{AlIgned}
\LabeL{SELF-genEraL}
\SigMa _{jk}(\oMega)= - \frac{2\pi q^2 \omEga_{0}}{M c^2} \, {\bf E}_{J}\cdOt  {\bf G}({\Bf r}_{j},{\bF r}_{k};\oMeGa) \cdoT {\bf e}_{k}.\eNd{aliGnEd}$$ The system Eq. (\[syStem\]) DetermineS eiGeNstAtEs of tHE ensemBle Of $N$ Dipoles In commoN RadIaTION aNd plasmon field, whilE cOMPlEx eigenvAlues oF ThE sELf-energy MaTriX givE EIgensTateS FrEquency sHifts wITh ReSpect to $\OmEga_{0}$ and ThEir DecAy ratES. In tHe abseNce of NP, rEal anD Imaginary parts OF the self-energY MaTRIx ARe diPolE-dipole inteRactIOn anD radIAtIon COupliNg betWeEN dIPoles $j$ and $k$, given by (in LoWest orDer in $Kr_{jk}$) $$\begin{aligNed}
\label{diPOLE-dipole}
&\DElta_{JK}^{0}=\fRAc{3\Gamma_{0
 \varphi_{j} - \left %(\om ega_{0}^2/ \omeg a\r igh t) \s in \ varphi_{j}, 
\ end{ aligned}$$ and the sol ution  o f  Eq. (\ [maxw ell-lap l ac e \ ])re ad s $ $\ b eg in{al ign ed}
\la bel{electr ic} 
{ \bf E}({\bfr }, \omega ) =  \f rac{4\pi d_{ 0}q  \omeg a^ {2} } {c^2}  \s um_{j }{\bfG }({\bf  r},{\bfr} _ {j},\o m ega) \c d o t{\bf  e}_{j} v_{j}, \e n d{ a ligned}$$ wher e ${\b fG }( { \ bfr}, {\bf r}',\ om ega)$ is thee le c t r icf ield Green di adic *in th e  pr esence  o f*N P. Fro m Eq.  ( \ [di spl\]), for  pho ton frequ ency c l ose tod ipolesfreque ncy , $ \ome g a\a ppr ox \om e ga _{0 } $,we obtai nacoupl ed s y s t e m of  eq uati ons f or normalized  di spla c eme nts,$$\be gin{ al igned }
\lab el{sy st em}
\sum_{k} \B igl[ (\omega_{ 0}- \o meg a) \de l ta_{jk } +   \ Sigma _ {jk}\Bi g r]v_ { k }  =  
%\frac{v_{0j}}{2 \o m e ga _{0}}=
\ frac{- i }{ 2} e^{-i\va rp hi_ {j}} , \ end{a lign e d} $$ where  $\Sig m a_{ jk}=\De lt a_{jk} -\ fra c{i }{2}\ G amma _{jk}$  is thecompl e x *self-energy matrix*, give n  b y $$ \ begi n{a ligned}
\la bel{ s elf- gene r al }
\ S igma_{jk} (\ o me g a)= - \frac{2\pi q^ 2\omega _{0}} {m c^2} \, {\ bf e}_{j}\ c d o t  {\bfG}({ \ bf r}_{j},{\bf r} _{k}; \omega) \c d ot {\bfe}_{k }.\end{a ligned}$$ T he syste m E q.(\[ sys t e m\ ]) determines e igen st ates of  th e ensem ble  of  $N $ d ip oles in c ommon ra di at io nand  plas m on field ,whi le  co mplex eigenv alues  ofth es elf -energy ma t r ix g iv eeige nst at es fr eque n cyshiftswith resp ect to $ \o me ga_{0}$  and their de ca y rates. I nthe  absen c e  of NP,real and imaginary part s  of the  se lf-en ergy  matrix a redipole -di p ole in teract ion a nd  ra d i ation c ou pli ng  between d i p ole s $j$  a nd $ k$, giv en by (in lowest o r der  in $kr_{jk}$ ) $ $\be g i n{ ali g ne d }
\ la b el{ d i pole-dipole}
&\ Delta_{jk} ^{ 0 }= \frac{3\Ga m ma_ {0 
 \varphi_{j}_- \left_%(\omega_{0}^2/\omega\right ) \sin \varphi_{j},_
 \end{aligned}$$_and_the solution_of_Eq. (\[maxwell-laplace\]) reads $$\begin{aligned}
\label{electric}
{\bf_E}({\bf r},\omega )_= \frac{4\pi d_{0}q \omega^{2}}{c^2}_\sum_{j}{\bf G}({\bf r},{\bf_r}_{j},\omega)_\cdot {\bf e}_{j} v_{j}, \end{aligned}$$ where ${\bf G}({\bf r},{\bf r}',\omega)$ is the electric field_Green_diadic *in_the_presence_of* NP. From Eq. (\[displ\]), for_photon frequency close to dipoles_frequency, $\omega_\approx \omega_{0}$, we obtain a coupled system of_equations_for normalized displacements,_$$\begin{aligned}
\label{system}
\sum_{k} \Bigl[(\omega_{0}-\omega ) \delta_{jk} +  \Sigma _{jk}\Bigr] v_{k}_= 
%\frac{v_{0j}}{2\omega_{0}}=
\frac{-i}{2} e^{-i\varphi_{j}},\end{aligned}$$ where $\Sigma _{jk}=\Delta_{jk}-\frac{i}{2}\Gamma_{jk}$_is the complex_*self-energy_matrix*,_given by $$\begin{aligned}
\label{self-general}
\Sigma _{jk}(\omega)=_- \frac{2\pi q^2 \omega_{0}}{m c^2} \,_{\bf e}_{j}\cdot  {\bf G}({\bf r}_{j},{\bf_r}_{k};\omega) \cdot {\bf e}_{k}.\end{aligned}$$ The system Eq. (\[system\])_determines eigenstates of the ensemble of_$N$ dipoles in common radiation_and plasmon_field, while complex eigenvalues of_the self-energy matrix_give eigenstates_frequency shifts with_respect to $\omega_{0}$ and their decay_rates. In the_absence of NP, real and imaginary_parts_of the self-energy_matrix_are_dipole-dipole interaction_and radiation coupling_between_dipoles $j$_and_$k$, given by (in lowest order_in_$kr_{jk}$) $$\begin{aligned}
\label{dipole-dipole}
&\Delta_{jk}^{0}=\frac{3\Gamma_{0