optical theorem and the concept of radiation damping to the full 6x6 tensorial polarizability of arbitrary magneto-electric scatterers. In this case, the work done per unit cycle by the incident field $\bm{E}_\mathrm{in}$ and $\bm{H}_\mathrm{in}$ to drive $\bm{p}$ and $\bm{m}$ is equal to $${W}=\quad \ll\mathrm{Re}\bm{E}_{in}\cdot\mathrm{Re}\frac{d\bm{p}}{dt}+\mathrm{Re}\bm{H}_{in}\cdot\mathrm{Re}\frac{d\bm{m}}{dt}\gg
$$ which evaluates to $${W}= \frac{2\pi}{Z} k \mbox{Im} \left[ \left(
\begin{array} {c c}
\bm{E}_\mathrm{in} & \bm{H}_\mathrm{in}
\end{array}\right)^{*}
\boldsymbol{\alpha} \left(
\begin{array} {c}
\bm{E}_\mathrm{in}\\
\bm{H}_\mathrm{in}
\end{array} \right)\right],
$$ where $(\cdot)^{*}$ indicates complex conjugate. The power per solid angle radiated by the induced dipoles in a direction specified by a unit vector $\bm{\hat{r}}$ is easily found by calculating the far-field Poynting vector from Eq. (\[Eq:dipoleG\]). The result is composed of three terms: $$\frac{dP}{d\Omega}=\frac{dP_p}{d\Omega}+\frac{dP_m}{d\Omega}+\frac{k^4}{2Z}\mathrm{Re}(\bm{p}\times
\bm{m})\cdot {\hat{\bm{r}}},\ \label{Eq:radiationpattern}$$ The first term in Eq. (\[Eq:radiationpattern\]) represents the scattered radiation of just the electric dipole $\bm{p}$, which integrates to a total scattered power given by Larmor’s formula. The second term in Eq. (\[Eq:radiationpattern\]) represents the radiation pattern of just the magnetic dipole $\bm{m}$, again given by Larmor’s formula. Note that both terms simply represent the well known $\sin^2\theta$ donut shaped radiation pattern for $\bm{p}$ and $\bm{m}$. The third term, however, can completely change the radiation pattern, as it contains the interference between the 