 has just three nonzero eigenvalues. Indeed, let us introduce new decay matrices as $$\gamma_{\mu\nu}^{r}=\dfrac{N\Gamma^{r}}{3} B_{\mu\nu},$$ where $$B_{\mu\nu}=\dfrac{3}{N}\sum_{j} e_{j}^{\mu}e_{j}^{\nu}$$ is $3\times 3$ matrix in coordinate space with ${\rm Tr}\hat{B}=3$. It is easy to see that ${\rm Tr}\bigl [ (\hat{\Gamma}^{r} )^{n}\bigr]={\rm Tr}\bigl [ (\hat{\gamma}^{r} )^{n}\bigr]$ for any integer $n$, i.e., the $N\times N$ matrix $\Gamma_{jk}^{r}$ has only three non-zero eigenvalues coinciding with those of matrix $\gamma_{\mu\nu}^{r}$ $$\label{gamma-coop}
\Gamma_{\mu}^{r}=\dfrac{N\Gamma^{r}}{3}\,\lambda_{\mu},$$ where $\lambda_{\mu}\sim 1$ are eigenvalues of $B_{\mu\nu}$. Note that the decay rates of the remaining $N-3$ subradiant states vanish in the long wave approximation used here; they acquire finite values in the next order in $(kr_{0})^{2}$.

Let us turn to non-radiative coupling, described by matrix $\Gamma_{jk}^{nr}$. Its diagonal elements, $\Gamma_{jj}^{nr}$, describe nonradiative energy exchange between excited dipole and NP plasmon modes with all angular momenta, as indicated by polarizabilities $\alpha''_{l}$ in Eq. (\[gamma-perp\]). The non-diagonal elements of $\Gamma_{jk}^{nr}$ describe a process by which a plasmon nonradiatively excited in the NP by dipole $k$ transfers its energy to another dipole $j$. In general, due to high-$l$ plasmons involved in nonradiative coupling, the eigenstates of $\Gamma_{jk}^{nr}$ are *different* from those of plasmon-enhanced radiative coupling $\Gamma_{jk}^{r}$ which contains only dipole ($l=1$) plasmon mode. Therefore, the eigenstates of full decay matrix $\Gamma_{jk}=\Gamma_{jk}^{nr}+\Gamma_{jk}^{nr}$ are *not* pure superradiant and subradiant states but their admixtures. However, the high-$l$ plasmon contribution to $\Gamma_{jk}^{nr}$ is significant only at
has just three nonzero eigenvalues. Indeed, let us insert newfangled decay matrices as $ $ \gamma_{\mu\nu}^{r}=\dfrac{N\Gamma^{r}}{3 } B_{\mu\nu},$$ where $ $ B_{\mu\nu}=\dfrac{3}{N}\sum_{j } e_{j}^{\mu}e_{j}^{\nu}$$ is $ 3\times 3 $ matrix in coordinate space with $ { \rm Tr}\hat{B}=3$. It is easy to witness that $ { \rm Tr}\bigl [ (\hat{\Gamma}^{r }) ^{n}\bigr]={\rm Tr}\bigl [ (\hat{\gamma}^{r }) ^{n}\bigr]$ for any integer $ n$, i.e., the $ N\times N$ matrix $ \Gamma_{jk}^{r}$ has only three non - zero eigenvalue coinciding with those of matrix $ \gamma_{\mu\nu}^{r}$ $ $ \label{gamma - coop } 
 \Gamma_{\mu}^{r}=\dfrac{N\Gamma^{r}}{3}\,\lambda_{\mu},$$ where $ \lambda_{\mu}\sim 1 $ are eigenvalues of $ B_{\mu\nu}$. bill that the decay rates of the remaining $ N-3 $ subradiant state vanish in the long wave estimate used here; they acquire finite values in the next order in $ (kr_{0})^{2}$. 

 Let us turn to non - radiative yoke, described by matrix $ \Gamma_{jk}^{nr}$. Its diagonal elements, $ \Gamma_{jj}^{nr}$, trace nonradiative energy exchange between aroused dipole and NP plasmon modes with all angular momentum, as indicated by polarizabilities $ \alpha''_{l}$ in Eq.   (\[gamma - perp\ ]). The non - diagonal elements of $ \Gamma_{jk}^{nr}$ describe a process by which a plasmon nonradiatively activated in the NP by dipole $ k$ transfers its energy to another dipole $ j$. In general, due to high-$l$ plasmons involved in nonradiative coupling, the eigenstates of $ \Gamma_{jk}^{nr}$ are * different * from those of plasmon - enhanced radiative coupling $ \Gamma_{jk}^{r}$ which contains only dipole ($ l=1 $) plasmon mode. Therefore, the eigenstates of full decay matrix $ \Gamma_{jk}=\Gamma_{jk}^{nr}+\Gamma_{jk}^{nr}$ are * not * pure superradiant and subradiant states but their admixtures. However, the high-$l$ plasmon contribution to $ \Gamma_{jk}^{nr}$ is significant merely at
 had just three nonzero eigtnvalues. Indeed, lgt us invroduce new decxy matrices as $$\gamma_{\mu\nu}^{r}=\dfrec{N\Gqmma^{r}}{3} B_{\mu\nu},$$ where $$B_{\mu\nu}=\dfraz{3}{N}\sum_{j} e_{j}^{\lu}e_{j}^{\nu}$$ iw $3\tinws 3$ matrix in coordinate dpacz xith ${\rm Tr}\hat{B}=3$. Lt is easy do see that ${\rm Tf}\bngl [ (\hat{\Gamma}^{r} )^{n}\bigr]={\rm Tr}\bigl [ (\hat{\gaima}^{r} )^{n}\bogg]$ for any inteder $m$, i.e., fhe $N\times N$ matrix $\Gamma_{jk}^{r}$ has ohly thrte non-zero eigenvakues coinciding with those of latrix $\gamma_{\mu\nu}^{r}$ $$\pabel{gamma-ciop}
\Gwnma_{\mu}^{r}=\dfrac{N\Eamma^{r}}{3}\,\lambda_{\mu},$$ where $\lzmbda_{\mu}\sim 1$ are eigenvalues of $C_{\mu\nu}$. Note that tye fgcay rates oh the gemaining $N-3$ smnradiatt statrs vanish in tme loig wqve approximation usev here; they acquire sinite vanuzs in the next order un $(kr_{0})^{2}$.

Lej us durn ro von-daviafive clupming, descrjbed by mateix $\Gamma_{jk}^{nr}$. Its doadinal elements, $\Gamma_{tj}^{gr}$, describe nonradiative energy exchangt betseen excited dipole and NP plasmon modes witj all angtlar momenta, as indicated by polarizabilities $\alpva''_{l}$ ii Dq. (\[gcnma-pefp\]). Hhe non-diagonal elements of $\Gamma_{jk}^{nr}$ describq a pgocess by which a plasmon nomrwdowtively excitgd in tks HP by dipole $k$ trajsfers yts ebergy to wnotner dipole $j$. In general, due to high-$l$ plcsmins involved in nouradiative cuuplong, tne eigenstates of $\Gamma_{lk}^{nr}$ ade *differenh* from thkre of plasmon-enhxncvd rddiative coupling $\Gamma_{jk}^{r}$ which coitainx only aipoke ($l=1$) pjasmon modf. Thevafore, the eigenstahes oy fuln decay mahrix $\Gamma_{jk}=\Gamma_{jk}^{nr}+\Gamma_{jk}^{nr}$ ece *not* pure sopesragiant anb subrsdiant states but their admnxtures. Kowevef, the high-$m$ plasmmn contribueion to $\Gamma_{bn}^{nr}$ is signihicant onjy ar
has just three nonzero eigenvalues. Indeed, let new matrices as B_{\mu\nu},$$ where $$B_{\mu\nu}=\dfrac{3}{N}\sum_{j} in space with ${\rm It is easy see that ${\rm Tr}\bigl [ (\hat{\Gamma}^{r} Tr}\bigl [ (\hat{\gamma}^{r} )^{n}\bigr]$ for any integer $n$, i.e., the $N\times N$ matrix has only three non-zero eigenvalues coinciding with those of matrix $\gamma_{\mu\nu}^{r}$ $$\label{gamma-coop} \Gamma_{\mu}^{r}=\dfrac{N\Gamma^{r}}{3}\,\lambda_{\mu},$$ $\lambda_{\mu}\sim are of Note that the decay rates of the remaining $N-3$ subradiant states vanish in the long wave used here; they acquire finite values in the order in $(kr_{0})^{2}$. Let turn to non-radiative coupling, described matrix Its diagonal $\Gamma_{jj}^{nr}$, nonradiative exchange between excited and NP plasmon modes with all angular momenta, as indicated by polarizabilities $\alpha''_{l}$ in Eq. (\[gamma-perp\]). The elements of a process which plasmon excited in the dipole $k$ transfers its energy to In general, due to high-$l$ plasmons involved in coupling, the of $\Gamma_{jk}^{nr}$ are *different* from those plasmon-enhanced radiative coupling $\Gamma_{jk}^{r}$ which contains only dipole plasmon mode. Therefore, the eigenstates of full decay matrix $\Gamma_{jk}=\Gamma_{jk}^{nr}+\Gamma_{jk}^{nr}$ are *not* pure superradiant and but their admixtures. However, high-$l$ plasmon contribution $\Gamma_{jk}^{nr}$ significant at
 has just three nonzero eigenvAlues. IndeeD, let uS inTroDuCe neW decAy matrices as $$\gaMMa_{\mu\Nu}^{r}=\dfrac{N\Gamma^{r}}{3} B_{\mu\nu},$$ whEre $$B_{\mU\nU}=\DfraC{3}{n}\sUm_{j} e_{j}^{\Mu}e_{j}^{\nu}$$ iS $3\TiMES 3$ maTrIx In cOoRDiNate sPacE with ${\rm tr}\hat{B}=3$. It is EasY tO see that ${\rm Tr}\BIgL [ (\hat{\Gamma}^{r} )^{N}\biGr]={\rm Tr}\bigl [ (\haT{\gaMma}^{r} )^{n}\bIgR]$ foR Any inTegEr $n$, i.e., The $N\tiMEs N$ matRix $\Gamma_{jK}^{r}$ HAs only THree non-ZERo EigeNvalues coinciding WItH Those of matrix $\gAmma_{\mu\Nu}^{R}$ $$\LaBEL{gaMma-Coop}
\Gamma_{\mU}^{r}=\Dfrac{n\gamma^{r}}{3}\,\lAMbDA_{\MU},$$ whERe $\lambda_{\mu}\sim 1$ Are eigenvalUEs oF $B_{\mu\nu}$. noTe tHAt the dEcay rAtES of The remaininG $N-3$ suBradiant sTates vANish in tHE long waVe apprOxiMatIon uSEd HeRe; tHeY AcqUIrE fiNIte Values in ThE nExt orDer iN $(KR_{0})^{2}$.

lEt us TurN to nOn-radIative couplinG, deScriBEd bY matrIx $\GamMa_{jk}^{Nr}$. its diAgonal ElemeNtS, $\Gamma_{jj}^{nr}$, descrIbe nOnradiatiVe eNeRgy ExChangE BetweeN exCitEd dipolE and NP pLAsmOn MODEs With all angular momeNtA, AS iNdicated By polaRIzAbILities $\alPhA''_{l}$ iN Eq. (\[gAMMa-perP\]). The NOn-Diagonal ElemenTS oF $\GAmma_{jk}^{nR}$ dEscribE a ProCesS by whICh a pLasmon NonradiaTivelY Excited in the NP BY dipole $k$ transFErS ITs ENergY to Another dipoLe $j$. IN GeneRal, dUE tO hiGH-$l$ plaSmons InVOlVEd in nonradiative couPlIng, the EigenStates of $\Gamma_{Jk}^{nr}$ are *difFEREnt* from tHose OF pLAsmon-enhanced rAdiatIve couplinG $\gamma_{jk}^{r}$ Which Contains Only dipolE ($L=1$) Plasmon mOde. theRefOre, THE eIgenstates of fULL decAy Matrix $\GAmmA_{jk}=\GammA_{jk}^{Nr}+\GAmmA_{jk}^{Nr}$ Are *not* purE superraDiAnT aNd SubRadiaNT states bUt TheIr AdmIxturES. HowevEr, the High-$L$ pLaSMon ContribUTiON To $\GaMmA_{jK}^{nr}$ iS siGnIficaNt onLY at
 has just three nonzero ei genvalues.  Inde ed,  le tus i ntro duce new decay matr ices as $$\gamma_{\mu\ nu}^{ r} = \dfr a c{ N\Gam ma^{r}} { 3} B _{\ mu \n u}, $$ wh ere $ $B_ {\mu\nu }=\dfrac{3 }{N }\ sum_{j} e_{j } ^{ \mu}e_{j}^ {\n u}$$ is $3\t ime s 3$ m at rix in co ord inate  space with $ {\rm Tr}\ ha t {B}=3$ .  It ise a sy  tosee that ${\rm Tr } \b i gl [ (\hat{\Ga mma}^{ r} )^ { n }\b igr ]={\rm Tr} \b igl [ (\hat{\ g am m a } ^{r }  )^{n}\bigr]$  for any in t ege r $n$,  i .e. ,  the $ N\tim es N$matrix $\Ga mma_ {jk}^{r}$  has o n ly thre e  non-ze ro eig env alu es c o in ci din gw ith th ose ofmatrix $ \g am ma_{\ mu\n u } ^ { r}$$$\ labe l{gam ma-coop}
\Gam ma_ {\mu } ^{r }=\df rac{N \Gam ma ^{r}} {3}\,\ lambd a_ {\mu},$$ where$\la mbda_{\mu }\s im  1$  a re ei g envalu esof$B_{\mu \nu}$.N ote  t h a t  t he decay rates ofth e re maining$N-3$s ub ra d iant sta te s v anis h in th e lo n gwave app roxima t io nused he re ; they  a cqu ire  fini t e va lues i n the ne xt or d er in $(kr_{0} ) ^{2}$.

Let u s  t u r nt o no n-r adiative co upli n g, d escr i be d b y  matr ix $\ Ga m ma _ {jk}^{nr}$. Its dia go nal el ement s, $\Gamma_{j j}^{nr}$,d e s cribe no nrad i at i ve energy exch angebetween ex c ited dip ole a nd NP pl asmon mod e s  with al l a ngu lar  mo m e nt a, as indicat e d  bypo larizab ili ties $\ alp ha' '_{ l}$  i n Eq. (\[ gamma-pe rp \] ).  T henon-d i agonal e le men ts  of  $\Ga m ma_{jk }^{nr }$ d es cr i bea proce s sb y  whi ch  a  pla smo nnonra diat i vel y excit ed in the  NP by d ip ol e $k$ t ransfers itsen ergy to an ot her  dipol e $j$. Ingeneral, due to high-$l $  plasmo nsinvol vedin nonrad iat ive co upl i ng, th e eige nstat es  of $ \Gamm a _ {j k}^ {n r}$ are *d i f fer ent*fr om t hose of  plasmon-enhancedr adi ative couplin g $ \Gam m a _{ jk} ^ {r } $ w hi c h c o n tains only dipo le ($l=1$)  p l as mon mode.T her ef ore, th e eigen state s  of ful l decay m atrix $\G am ma_{ j k }=\ Gamma_{jk} ^{nr}+\G amma_{jk} ^ {nr}$ ar e *no t*pure s up err adian t ands ubr adian t stat es  but t heirad mixtures . However, the high-$l$  plasm on co ntr ibution t o $ \ Gam ma_{jk}^{ nr}$  is signif ica ntonlyat
 has_just three_nonzero eigenvalues. Indeed, let_us introduce_new_decay matrices_as_$$\gamma_{\mu\nu}^{r}=\dfrac{N\Gamma^{r}}{3} B_{\mu\nu},$$ where_$$B_{\mu\nu}=\dfrac{3}{N}\sum_{j} e_{j}^{\mu}e_{j}^{\nu}$$ is_$3\times 3$ matrix in_coordinate space with_${\rm_Tr}\hat{B}=3$. It is easy to see that ${\rm Tr}\bigl [ (\hat{\Gamma}^{r} )^{n}\bigr]={\rm Tr}\bigl [_(\hat{\gamma}^{r}_)^{n}\bigr]$ for_any_integer_$n$, i.e., the $N\times N$_matrix $\Gamma_{jk}^{r}$ has only three_non-zero eigenvalues_coinciding with those of matrix $\gamma_{\mu\nu}^{r}$ $$\label{gamma-coop}
\Gamma_{\mu}^{r}=\dfrac{N\Gamma^{r}}{3}\,\lambda_{\mu},$$ where_$\lambda_{\mu}\sim_1$ are eigenvalues_of $B_{\mu\nu}$. Note that the decay rates of the_remaining $N-3$ subradiant states vanish in_the long wave_approximation_used_here; they acquire finite_values in the next order in_$(kr_{0})^{2}$.

Let us turn to non-radiative coupling,_described by matrix $\Gamma_{jk}^{nr}$. Its diagonal elements,_$\Gamma_{jj}^{nr}$, describe nonradiative energy exchange between_excited dipole and NP plasmon_modes with_all angular momenta, as indicated_by polarizabilities $\alpha''_{l}$_in Eq. (\[gamma-perp\])._The non-diagonal elements_of $\Gamma_{jk}^{nr}$ describe a process by_which a plasmon_nonradiatively excited in the NP by_dipole_$k$ transfers its_energy_to_another dipole_$j$. In general,_due_to high-$l$_plasmons_involved in nonradiative coupling, the eigenstates_of_$\Gamma_{jk}^{nr}$ are *different* from those of plasmon-enhanced_radiative coupling $\Gamma_{jk}^{r}$ which_contains_only dipole ($l=1$) plasmon_mode. Therefore, the eigenstates of_full decay matrix $\Gamma_{jk}=\Gamma_{jk}^{nr}+\Gamma_{jk}^{nr}$ are *not*_pure superradiant_and subradiant_states but their admixtures. However, the high-$l$ plasmon contribution to $\Gamma_{jk}^{nr}$_is significant only at