## 1.

## Introduction

Electromagnetic scattering from nanometer-scale objects has long been a topic of large interest and relevance to fields from astrophysics or meteorology to biophysics, medicine, and material science.^{1, 2, 3, 4} In the last few years, small particles with resonant magnetic properties are being explored as constitutive elements of new metamaterials and devices. Magnetic effects, however, cannot be easily exploited in the visible or infrared regions due to intrinsic natural limitations of optical materials and the quest for magnetic plasmons and magnetic resonant structures at optical frequencies^{5, 6, 7} has then been mainly focused on metallic structures. The unavoidable problems of losses and saturation effects inherent to these metamaterials in the optical and near-infrared regimes have stimulated the study of high-permittivity particles as their constitutive elements:^{8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18} For very large permittivities, small spherical particles present well defined sharp resonances;^{1} either electric or magnetic resonant responses can then be tuned by choosing the appropriate sphere radius. However, the applicability of conventional homogenization procedures to determine the constitutive parameters of composites of resonant spheres cannot be implemented without due regard to their limitations on volume fraction and size parameter.^{19}

In the presence of both electric and magnetic properties, the scattering characteristics of a small object present marked differences with respect to pure electric or magnetic responses. Even in the simplest case of small or dipolar scatterers, remarkable scattering effects of magnetodielectric particles were theoretically established by Kerker^{20} concerning suppression or minimization of either forward or backward scattering. Intriguing applications in scattering cancellation and cloaking^{21, 22, 23} and magneto-optical systems,^{24, 25, 26, 27} together with the unusual properties of the optical forces on magnetodielectric particles,^{28, 29} have renewed interest in the field.

The striking characteristics of the scattering diagram of small (Rayleigh) magnetodielectric particles^{20, 30, 31} were obtained assuming arbitrary values of relative permittivity and permeability.

Nevertheless, no concrete example of such particles that might present those interesting properties in the visible or infrared regions had been proposed. Very recently, it has been shown^{32, 33} that submicrometer silicon spheres present dipolar magnetic and electric responses, characterized by their respective first-order Mie coefficient, in the near-infrared, in such a way that either of them can be selected by choosing the illumination wavelength. In a later work, it has also been shown^{29, 34} that Si spheres constitute such a previously quested real example of a dipolar particle with either electric and/or magnetic response, of consequences both for their emitted intensity and behavior under electromagnetic forces.

These properties should not be restricted to Si particles, but should also apply to other dielectric materials with relatively moderate refraction index. In the present work, we discuss the effects associated to the interference between electric and magnetic dipoles in germanium spheres. The paper is structured as follows. Based on the exact Mie theory, in Sec. 2 we show that both the extinction cross section and the scattering diagrams of submicrometer Ge spheres in the infrared region can be well described by dipolar electric and magnetic fields, being quadrupolar and higher order contributions negligible in this frequency range. Specifically, the scattering diagrams calculated at the generalized Kerker's (GK) conditions are shown to be equivalent to those previously reported^{20, 31} for hypothetical magnetodielectric particles. Following previous work regarding the peculiar properties of optical forces at GK conditions,^{29} in Sec. 3 we analyze the consequences of the strong scattering anisotropy on the radiation pressure on Ge particles.

## 2.

## Extinction Cross Sections and Scattering Anisotropy of Submicrometer Germanium Spheres

Consider a nonabsorbing dielectric sphere of radius *a* and permittivity [TeX:] \documentclass[12pt]{minimal}\begin{document}$\epsilon _p=m_p^2$\end{document}
${\epsilon}_{p}={m}_{p}^{2}$
in an otherwise uniform medium with real relative permittivity ε and refractive index [TeX:] \documentclass[12pt]{minimal}\begin{document}$m=\sqrt{\epsilon \mu }$\end{document}
$m=\sqrt{\epsilon \mu}$
. The relative permeability μ* _{p}* of the sphere and of the surrounding medium, μ, are both assumed to be μ

*= μ = 1 (even though in some of the equations to be written in this paper, μ will explicitly be included for the sake of comprehensiveness of the theory) .*

_{p}Under plane wave illumination, and assuming linearly polarized light, the incident wave is described by

## 1

[TeX:] \documentclass[12pt]{minimal}\begin{document} \begin{equation} E = E_0 {\bf u}_x e^{ikz} e^{-i\omega t}, \quad B = B_0 {\bf u}_y e^{ikx} e^{-i\omega t} \end{equation} \end{document} $E={E}_{0}{\mathbf{u}}_{x}{e}^{ikz}{e}^{-i\omega t},\phantom{\rule{1em}{0ex}}B={B}_{0}{\mathbf{u}}_{y}{e}^{ikx}{e}^{-i\omega t}$*k*=

*m*ω/

*c*=

*m*2π/λ, λ being the wavelength in vacuum and

*B*

_{0}=

*mE*

_{0}. The field scattered by the sphere can be decomposed into a multipole series (that which is called Mie's expansion) characterized by the {

*a*

*} electric and {*

_{n}*b*

*} magnetic Mie coefficients (being*

_{n}*a*

_{1}and

*b*

_{1}proportional to the electric and magnetic dipoles,

*a*

_{2}and

*b*

_{2}to the quadrupoles, and so on). Mie coefficients are related to the scattering phase-shifts ξ

*and β*

_{n}*through*

_{n}^{1}

## 2

[TeX:] \documentclass[12pt]{minimal}\begin{document} \begin{eqnarray} a_n &=& \frac{1}{2}(1-e^{-2i\xi _n }) = i \sin \xi _n e^{-i\xi _n}, \end{eqnarray} \end{document} $\begin{array}{ccc}\hfill {a}_{n}& =& \frac{1}{2}(1-{e}^{-2i{\xi}_{n}})=i\mathrm{sin}{\xi}_{n}{e}^{-i{\xi}_{n}},\hfill \end{array}$## 3

[TeX:] \documentclass[12pt]{minimal}\begin{document} \begin{eqnarray} b_n &=& \frac{1}{2}(1-e^{-2i\beta _n }) = i \sin \beta _n e^{-i\beta _n}. \end{eqnarray} \end{document} $\begin{array}{ccc}\hfill {b}_{n}& =& \frac{1}{2}(1-{e}^{-2i{\beta}_{n}})=i\mathrm{sin}{\beta}_{n}{e}^{-i{\beta}_{n}}.\hfill \end{array}$*Q*

_{ext}, and scattering,

*Q*

*, cross sections of a dielectric sphere, expressed in the Mie coefficients, read*

_{S}## 4

[TeX:] \documentclass[12pt]{minimal}\begin{document} \begin{eqnarray} Q_{\rm{ext}} &=& \frac{2\pi }{k^2}\sum _{n=1}^\infty \left(2n+1\right) \mbox{Re}\lbrace a_n + b_n \rbrace , \end{eqnarray} \end{document} $\begin{array}{ccc}\hfill {Q}_{\mathrm{ext}}& =& \frac{2\pi}{{k}^{2}}\sum _{n=1}^{\infty}\left(2n+1\right)\text{Re}\{{a}_{n}+{b}_{n}\},\hfill \end{array}$## 5

[TeX:] \documentclass[12pt]{minimal}\begin{document} \begin{eqnarray} Q_S &=& \frac{2\pi }{k^2} \sum _{n=1}^\infty (2n+1) \lbrace |a_n|^2+|b_n|^2 \rbrace . \end{eqnarray} \end{document} $\begin{array}{ccc}\hfill {Q}_{S}& =& \frac{2\pi}{{k}^{2}}\sum _{n=1}^{\infty}(2n+1)\left\{\right|{a}_{n}{|}^{2}+|{b}_{n}{|}^{2}\}.\hfill \end{array}$*m*,

*m*

*,*

_{p}*Q*

_{ext}=

*Q*

*.*

_{S}In the small particle limit (*x* ≡ *ka* ≪ 1) and large particle permittivities (*m** _{p}*/

*m*≫ 1) the extinction cross section presents characteristic sharp resonance peaks. At each resonance, the extinction cross section is of the order of λ

^{2}and it is independent of the particle size or refractive index.

^{1}Interesting applications of well defined Mie resonances, associated to low loss and large dielectric constants, are accessible for different materials at microwave and terahertz frequencies. However, as

*m*

*decreases there is an increasing overlap between the cross section peaks, and the resonant character weakens until it disappears. Since usually nonabsorbing materials present low refractive index in the near-infrared and visible frequency ranges, Mie resonances of small particles in these regimes have not been considered in detail. However, a recent analysis of the cross section of submicrometer dielectric particles*

_{p}^{33}show that well defined resonances can be found for materials with relative refractive index as low as

*m*

*∼ 3*

_{p}*m*.

Germanium submicrometer particles are a good candidate to explore the effects associated to the interference between electric and magnetic responses. In the micrometer wavelength regime, within the transparent region of germanium (λ ≳ 1.4 μm), the refraction index can be well approximated by a real constant [TeX:] \documentclass[12pt]{minimal}\begin{document}$m_p \approx \sqrt{16} = 4$\end{document}
${m}_{p}\approx \sqrt{16}=4$
(see, for example, Ref. 35). The calculated exact extinction (or scattering) cross section of a 240 nm radius Ge sphere in vacuum (*m* = 1) is plotted in Fig. 1. Although there is an overlap between the first dipolar peaks, the first dipolar magnetic resonance (at λ ≈ 2000 nm) and the electric dipolar resonance (at λ ≈ 1500 nm) are still very well defined. As we will see, this overlap plays a key role in determining the peculiar scattering diagrams of Ge particles. Interestingly, for wavelengths larger than λ ≈ 1400 nm, the cross section is completely determined by the first *b*_{1} and *a*_{1} coefficients. Therefore, Eq. 4 becomes

## 6

[TeX:] \documentclass[12pt]{minimal}\begin{document} \begin{eqnarray} Q_{\rm{ext}}=Q_S \approx \frac{6\pi }{k^2} \mbox{Re}\lbrace a_1 + b_1 \rbrace \end{eqnarray} \end{document} $\begin{array}{c}\hfill {Q}_{\mathrm{ext}}={Q}_{S}\approx \frac{6\pi}{{k}^{2}}\text{Re}\{{a}_{1}+{b}_{1}\}\end{array}$*a*

_{1}and

*b*

_{1}are the exact (nonapproximated) dipolar electric and magnetic Mie coefficients, respectively.

^{1}In other words, in this regime submicrometer Ge particles can be treated as dipolar particles in the sense that only the dipolar term of Eq. 4 contributes to the scattered field.

Dipolar particles are usually characterized by their electric and/or magnetic complex polarizabilities, α* _{e}* and α

*, which may be written in the form*

_{m}^{28, 33}

## 7

[TeX:] \documentclass[12pt]{minimal}\begin{document} \begin{eqnarray} \displaystyle\alpha _e &=& 3i \epsilon a_1 /(2 k ^3) = \frac{\alpha _e^{(0)}}{ 1- i \displaystyle\frac{2}{3 \epsilon } k ^3 \alpha _e^{(0)}}, \end{eqnarray} \end{document} $\begin{array}{ccc}\hfill {\displaystyle {\alpha}_{e}}& =& 3i\epsilon {a}_{1}/\left(2{k}^{3}\right)=\frac{{\alpha}_{e}^{\left(0\right)}}{{\displaystyle 1-i\frac{2}{3\epsilon}{k}^{3}{\alpha}_{e}^{\left(0\right)}}},\hfill \end{array}$## 8

[TeX:] \documentclass[12pt]{minimal}\begin{document} \begin{eqnarray} \displaystyle\alpha _m &=& 3 i b_1 /(2\mu k ^3)=\frac{\alpha _m^{(0)}}{1- i \displaystyle\frac{2}{3} \mu k ^3 \alpha _m^{(0)}}. \end{eqnarray} \end{document} $\begin{array}{ccc}\hfill {\displaystyle {\alpha}_{m}}& =& 3i{b}_{1}/\left(2\mu {k}^{3}\right)=\frac{{\alpha}_{m}^{\left(0\right)}}{{\displaystyle 1-i\frac{2}{3}\mu {k}^{3}{\alpha}_{m}^{\left(0\right)}}}.\hfill \end{array}$## 9

[TeX:] \documentclass[12pt]{minimal}\begin{document} \begin{eqnarray} \alpha _e^{(0)} = - \frac{3 \epsilon }{2 k^3} \tan \xi _1, \quad \alpha _m^{(0)} = - \frac{3}{2 \mu k^3} \tan \beta _1. \end{eqnarray} \end{document} $\begin{array}{c}\hfill {\alpha}_{e}^{\left(0\right)}=-\frac{3\epsilon}{2{k}^{3}}\mathrm{tan}{\xi}_{1},\phantom{\rule{1em}{0ex}}{\alpha}_{m}^{\left(0\right)}=-\frac{3}{2\mu {k}^{3}}\mathrm{tan}{\beta}_{1}.\end{array}$For a pure electric or a pure magnetic dipole, in the absence of interferences, the far field radiation pattern is symmetrically distributed between forward and backward scattering. However, when we consider the coherent contribution of both electric and magnetic dipoles, the radiation pattern is mainly distributed in the forward or backward region according to whether [TeX:] \documentclass[12pt]{minimal}\begin{document}$\Re (\alpha _{e} \alpha _{m}^{*})$\end{document}
$\Re \left({\alpha}_{e}{\alpha}_{m}^{*}\right)$
is positive or negative, respectively.^{28, 29} Interestingly, at the GK conditions |ε^{−1}α* _{e}*|

^{2}= |μα

*|*

_{m}^{2}, the scattered intensity should be independent of the incident polarization angle:

## 10

[TeX:] \documentclass[12pt]{minimal}\begin{document} \begin{eqnarray} \frac{dQ_S}{d\Omega } (\theta )= k ^4 | \epsilon ^{-1} \alpha _e\|^2 ( 1 + \cos ^2 \theta ) + 2 k^4 \frac{\mu }{\epsilon } \Re (\alpha _e \alpha _m^* ) \cos \theta . \end{eqnarray} \end{document} $\begin{array}{c}\hfill \frac{d{Q}_{S}}{d\Omega}\left(\theta \right)={k}^{4}{|{\epsilon}^{-1}{\alpha}_{e}\Vert}^{2}(1+{\mathrm{cos}}^{2}\theta )+2{k}^{4}\frac{\mu}{\epsilon}\Re \left({\alpha}_{e}{\alpha}_{m}^{*}\right)\mathrm{cos}\theta .\end{array}$i) The intensity in the backscattering direction can be exactly zero whenever *a*_{1} = *b*_{1} or

## 11

[TeX:] \documentclass[12pt]{minimal}\begin{document} \begin{eqnarray} \epsilon ^{-1} \alpha _e = \mu \alpha _m; \quad \frac{dQ_S}{d\Omega }(180^\circ ) = 0 . \end{eqnarray} \end{document} $\begin{array}{c}\hfill {\epsilon}^{-1}{\alpha}_{e}=\mu {\alpha}_{m};\phantom{\rule{1em}{0ex}}\frac{d{Q}_{S}}{d\Omega}\left({180}^{\circ}\right)=0.\end{array}$^{20}ε

*/ε = μ*

_{p}*/μ. Equation 11 can be considered as a generalized first Kerker condition for arbitrary dipolar particles.*

_{p}ii) Although the intensity cannot be exactly zero in the forward direction (causality imposes ℑ{α* _{e}*}, ℑ{α

*} > 0), in absence of particle absorption, the forward intensity presents a minimum at*

_{m}^{29}

## 12

[TeX:] \documentclass[12pt]{minimal}\begin{document} \begin{eqnarray} &\displaystyle\Re \lbrace \epsilon ^{-1} \alpha _e \rbrace = -\Re \lbrace \mu \alpha _m\rbrace \quad \mbox{ and } \quad \Im \lbrace \epsilon ^{-1} \alpha _e \rbrace = \Im \lbrace \mu \alpha _m\rbrace & \nonumber \\ &\displaystyle\frac{dQ_S}{d\Omega }(0^\circ ) = \frac{16}{9} k^{10} |\epsilon ^{-1} \alpha _e|^4.& \end{eqnarray} \end{document} $\begin{array}{ccc}& {\displaystyle \Re \left\{{\epsilon}^{-1}{\alpha}_{e}\right\}=-\Re \left\{\mu {\alpha}_{m}\right\}\phantom{\rule{1em}{0ex}}\phantom{\rule{4.pt}{0ex}}\text{and}\phantom{\rule{4.pt}{0ex}}\phantom{\rule{1em}{0ex}}\Im \left\{{\epsilon}^{-1}{\alpha}_{e}\right\}=\Im \left\{\mu {\alpha}_{m}\right\}}& \hfill \\ & {\displaystyle \frac{d{Q}_{S}}{d\Omega}\left({0}^{\circ}\right)=\frac{16}{9}{k}^{10}{\left|{\epsilon}^{-1}{\alpha}_{e}\right|}^{4}.}\end{array}$*a*

_{1}and

*b*

_{1}can be approximated by the first terms in the expansions of these coefficients in powers of the size parameter,

*x*) this happens at the second Kerker condition:

^{20, 31}ε

*/ε = −(μ*

_{p}*/μ − 4)/(2μ*

_{p}*/μ + 1). Although the original derivation*

_{p}^{20}was obtained in the quasi-static approximation (ℑ{α} ≈ 0, thus leading to

*dQ*(0°)/

*d*Ω = 0), the actual intensity for a very small particle goes as ∼(

*ka*)

^{10}which, for Rayleigh particles, would be negligible.

^{23}However, the derivation of the GK conditions 11, 12 with unusual scattering properties was obtained

^{29}with the unique assumption that the radiation fields are well described by dipolar electric and magnetic fields (i.e., using the exact expressions for

*a*

_{1}and

*b*

_{1}which adequately describe the scattered field). This goes well beyond the Rayleigh limit and should apply to any dipolar particle. Specifically, the GK conditions also apply to purely dielectric spheres (μ

*= 1) providing that their scattering properties may be described by the two first terms in the Mie expansion.*

_{p}Figure 2 shows the real and imaginary parts of the polarizabilities [Eq. 8] for a Ge sphere with *a* = 240 nm. One sees the values of λ at which ℑ{α* _{e}*} = ℑ{α

*}, which are where the first and second GK conditions hold for these polarizabilities.*

_{m}In order to check the validity of the dipolar approach and confirm the predictions at the GK conditions, we have computed the exact scattering diagram from the full Mie expansion. The numerical results are shown in Fig. 3. While the backward intensity drops to zero at the first GK condition wavelength, at the second condition, although the most of the intensity goes backward, the scattering diagram presents a very small peak in the forward direction. As expected from the extinction cross section, the far field pattern is fully consistent with the dipolar approximation. In particular, at the GK conditions it is independent of the incoming polarization. Dielectric spheres and, in particular, lossless Ge particles in the near-infrared, then constitute a realizable laboratory to observe such interesting scattering effects.

## 3.

## Effects on Optical Forces

Let us now discuss some of the consequences of the GK condition on optical forces. The theory of the force on a dipolar magnetic particle has recently been developed.^{28, 36} This includes pure dielectric particles which can be well described by its first two electric and magnetic Mie coefficients.^{28} The time averaged force on a dipolar particle can be written as the sum of three terms:^{28}

## 13

[TeX:] \documentclass[12pt]{minimal}\begin{document} \begin{eqnarray} \langle {\bf F}\rangle &=& \langle {\bf F}_{e}\rangle + \langle {\bf F}_{m}\rangle + \langle {\bf F}_{e-m}\rangle \nonumber \\[4pt] &=& {\bf u}_z F_0 \left[ \frac{1}{a^3}\Im \lbrace \epsilon ^{-1} \alpha _e + \mu \alpha _m \rbrace - \frac{2k^3}{3a^3}\frac{\mu }{\epsilon }\Re (\alpha _e \alpha _m^{\ast })\right], \end{eqnarray} \end{document} $\begin{array}{ccc}\hfill \u27e8\mathbf{F}\u27e9& =& \u27e8{\mathbf{F}}_{e}\u27e9+\u27e8{\mathbf{F}}_{m}\u27e9+\u27e8{\mathbf{F}}_{e-m}\u27e9\hfill \\ & =& {\mathbf{u}}_{z}{F}_{0}\left[\frac{1}{{a}^{3}}\Im \{{\epsilon}^{-1}{\alpha}_{e}+\mu {\alpha}_{m}\}-\frac{2{k}^{3}}{3{a}^{3}}\frac{\mu}{\epsilon}\Re \left({\alpha}_{e}{\alpha}_{m}^{*}\right)\right],\hfill \end{array}$*F*

_{0}≡ ε

*ka*

^{3}|

*E*

_{0}|

^{2}/2. The first two terms, 〈

**F**

*〉 and 〈*

_{e}**F**

*〉, correspond to the forces on the induced pure electric and magnetic dipoles, respectively. 〈*

_{m}**F**

*〉, due to the interaction between both dipoles,*

_{ − m}^{28, 36}is related to the asymmetry in the scattered intensity distribution [cf. the last term in Eq. 10].

^{28}

At the first GK condition, Eq. 11, the interference term of Eq. 13 cancels out the magnetic contribution and we obtain 〈**F**〉 = 〈**F*** _{e}*〉. At the second GK condition, Eq. 12, where the backscattering is enhanced, 〈

**F**〉 = 3〈

**F**

*〉. Notice that at both GK conditions the scattering cross section is exactly the same; however, the radiation pressures differ by a factor of 3. These properties are illustrated in Fig. 4, where we show the different contributions to the total time averaged force on a submicrometer Ge particle.*

_{e}The strong peak in the radiation pressure force is mainly dominated by the first “magnetic” Mie resonance. This is striking and in contrast with all previous beliefs about optical forces on dipolar dielectric particles, that assumed that these forces would solely be described by the electric polarizability. It is also common to assume that for dielectric particles the real part of the polarizability is much larger than its imaginary part. As a matter of fact, this is behind the development of optical tweezers, in which gradient forces (that are proportional to ℜ{α* _{e}*}), dominate over the radiation pressure or scattering force contribution (which is proportional to ℑ{α

*}).*

_{e}^{37}However, as the size of the particle increases, and for any dielectric particle, there is a crossover from electric to magnetic response as we approach the first Mie resonance, the point at which the response is absolutely dominated by the magnetic dipole. Moreover, just at the resonance, and in absence of absorption, ℜ{α

*} = 0 and ℑ{α*

_{m}*} = 3/(μ2*

_{m}*k*

^{3}). Then, the radiation pressure contribution of the magnetic term dominates the total force 〈

**F**〉 ≃ 〈

**F**

*〉 ≈ (ε|*

_{m}*E*

_{0}|

^{2}/2)[3/(2

*k*

^{2})]. Namely, in resonance the radiation pressure force presents a strong peak, the maximum force being independent of both material parameters and particle radius.

## 4.

## Conclusions

In summary, we have predicted that real small dielectric particles made of nonmagnetic materials present interesting scattering properties similar to those previously reported for somewhat hypothetical magnetodielectric particles,^{20} resulting from an interplay between real and imaginary parts of both electric and magnetic polarizabilities. The exact scattering diagram, computed from the full Mie expansion, of submicrometer Ge particles in the infrared, was shown to be consistent with the expected result for dipolar electric and magnetic scattering. Then we showed that these unusual scattering effects also affect the radiation pressure on these small particles; specifically, the generalized Kerker's conditions have been tested on Ge spheres. Submicro-meter Ge particles constitute an excellent laboratory to observe such new scattering phenomena and force effects. Being Ge permittivitys is higher than Si, the present work extends the range of some possible applications (previously suggested for silicon particles), to systems where the host medium presents refractive index larger than vacuum.

The extraordinary scattering properties discussed here will strongly affect the dynamics of particle confinement in optical traps and vortex lattices^{38} governed by both gradient and curl forces.^{28, 39} The interference between electric and magnetic dipoles suggest also intriguing possibilities regarding resonant coupling between the scattered dipolar field and guided modes in confined geometries.^{40} We do believe that our results will stimulate further experimental and theoretical work in different directions, from optical trapping and particle manipulation to cloaking and the design of optical metamaterials based on lossless dielectric particles.

## Acknowledgments

We appreciate interesting discussions with P. Albella, L. S. Froufe-Pérez, A. García-Etxarri, and J. M. Saiz. This work has been supported by the EU NMP3-SL-2008-214107-Nanomagma, the Spanish MICINN Consolider NanoLight (CSD2007-00046), FIS2010-21984, FIS2009-13430-C01-C02, and FIS2007-60158, as well as by the Comunidad de Madrid Microseres-CM (S2009/TIC-1476). Work by R.G.-M. was supported by the MICINN “Juan de la Cierva” Fellowship.