Print version ISSN 0103-9733
Braz. J. Phys. vol.31 no.1 São Paulo Mar. 2001
Multiple scattering of electromagnetic waves in disordered magnetic media: localization parameter, energy transport velocity and diffusion constant
F. A. Pinheiro1, A. S. Martinez2, and L. C. Sampaio1
1Centro Brasileiro de Pesquisas Físicas/CNPq,
Rua Dr. Xavier Sigaud, 150, CEP 22290-180, Rio de Janeiro, RJ, Brazil
2Faculdade de Filosofia, Ciências e Letras de Ribeirão Preto, Universidade de São Paulo,
Av. Bandeirantes 3900, CEP 14040-901, Ribeirão Preto, SP, Brazil
Received on 16 August, 2000
We review some of our recent results concerning the single and multiple eletromagnetic scattering by magnetic spherical particles. For a single electromagnetic scattering we show that the magnetic contribution alters, when compared to nonmagnetic scattering, the behavior of the cross sections and mean cosine of the scattering angle (ácos qñ). For ferromagnetic particles, resonances may occur even in the small-particle limit when the particle radius is much smaller than the wavelength. The resonances increase the cross-sections while ácos qñ is diminished, and even may become negative. Several quantities such the Ioffe-Regel parameter for localization are calculated for the multiple scattering regime. We show that magnetic scattering favors the observation of localization of electromagnetic waves in three dimensions. Further, this is also verified for dynamical experiments, where we show that the diffusion constant can be very small. Since the magnetic permeability of the scatterers can vary significatively around the Curie-Weiss critical temperature, experiments should be done on the microwave region and the scatterers should be soft ferrites. Some aspects of such experiment are presented.
The interest in studying the propagation of electromagnetic (EM) waves in random media has experienced a considerable increase over the last years. After the observation of the coherent backscattering effect,  the optical counterpart of electronic weak localization, it has been soon realized that many electronic effects in condensed matter physics have their analog in EM wave systems. Among these, are the photonic Hall effect,  the anisotropic light diffusion  and the localization of light. Although these phenomena characterize a close analogy between EM waves and electrons, dynamical situations concerning EM and electron propagation reveal fundamental differences between them.  In fact, van Albada et al.[6, 7] have shown that, the energy transport velocity vE, which characterizes the dynamics of light diffusion through the disordered media, is not equal to the phase velocity vp as in the electronic case. This is because a correction factor in vE (associated with the stored EM energy inside each scatterer due to Mie resonances) should be taken into account in the optical case, in contrast to electronic systems since for electron-impurity scattering this correction is absent.
However, in all multiple scattering phenomena reported so far, including the destruction of the coherent backscattering effect by magneto-optical active materials,[8, 9] it is assumed that the scatterers are nonmagnetic. The present paper is thus devoted to investigate some aspects of single and multiple scattering and diffusion of EM waves through magnetic random media.[10, 11, 12]
This paper is organized as follows. In Sec. II we treat both analytically and numerically the EM scattering by a single magnetic scatterer, assumed to be a sphere. It is important to point that several of the main properties of multiple scattering in magnetic media are essentially due to unusual features which characterizes magnetic single scattering, such as preferential backscattering and resonance effects even in the small-particle limit. In Sec. III these results concerning single scattering will be inserted in the multiple scattering context, where we will investigate the consequences of the introduction of magnetic scatterers in the localization parameter and coherent backscattering effect, in the energy transport velocity and in the diffusion constant. We have shown that the presence of magnetic scatterers induces a global decrease in the localization parameter as well as oscillations in the diffusion constant. Sec. IV is devoted to the discussion of some experimental aspects and possible applications of our results. Exploring the fact that the magnetic permeability of magnetic scatterers is strongly dependent of an external magnetic field Hext and typically follows the Curie-Weiss law, we suggest that both the localization parameter and the diffusion constant can be tuned by varying either Hext or the temperature. Finally, some comments and conclusions are addressed in Sec.V.
II Single Scattering
In order to investigate the influence of magnetic scatterers on the localization parameter k* (where k is the wavenumber and * is the transport mean free path), it is essential to consider the single scattering of EM waves by magnetic particles. In single Mie scattering problem, a plane wave with wavelength l is scattered by a homogeneous magnetic sphere of radius a, complex magnetic permeability ms and complex electric permittivity s, leading to a complex refraction index ns = (s ms)1/2. The medium has a real electric permittivity m (it does not absorb nor emit light) and vacuum magnetic permeability m0, which leads to a real refractive index nm = (m mm)1/2. The relative refractive index is m = ns/nm = (m)1/2, where = ¢ + i¢¢ and m = m¢ + im¢¢ are the complex relative electric permittivity and the complex relative magnetic permeability, respectively. Both the scatterers and the medium are considered to be nonmagneto-optical active.
Although in standard texts on EM scattering [13, 14, 15,] the relative magnetic permeability m is usually assumed to be unitary, this is no longer valid for magnetic particles. In ferromagnetic materials m can assume a rather large spectrum of values. For instance, for a nickel bulk at room temperature, m is of the order of 106 for incident radiation in the microwave range.
The quantities of interest in single scattering are the total cross-section s and anisotropy factor ácosqñ, which can be expressed in terms of the Mie coefficients an and bn. [13, 14, 15] In the small-particle limit (ka << 1), we have, in the lowest order in ka,
where = m/m = (/m)1/2 and yn(z) and yn¢(z) are the Ricatti-Bessel function and its derivative with respect to the argument, respectively. It is important to stress the fact that small-particle limit does not mean pointlike scatterers.
To obtain the Rayleigh limit, one has also to ensure that the size of the scatterer is small compared to the wavelength inside the particle (|m|ka << 1), which leads to a1 = g(1-)/(2+) and b1 = g(1-m)/(2+m), where g = -2i(ka)3/3. Notice that all the dependence on m is in the b1 term, which vanishes in the nonmagnetic case. In this case, the leading term in the b1 expansion is of the order of (ka)5. The anisotropy factor in terms of lowest order in ka is ácosqñ = Â(a1)/(|a1|2 + |b1|2). Notice that in magnetic scattering ácosqñ does not vanish, in contrast to the nonmagnetic case.
It is important to point out that, when |m| 1, which is the case of ferromagnetic materials, the small-particle limit must be treated in a different form. In this case, the limit |m|ka 1 must be considered using yn(z) ~ cos[z - (n + 1)/2p] in Eqs. (1), leading to:
where p(m, ka) = ka tan(mka). Differently from the Rayleigh limit, the a1 term also has a magnetic contribution for large values of |m|. The periodic function tan(mka), which depends on m, is responsible for resonances in ácosqñ. Furthermore, these oscillations are due to the real part of m, since tan(ix) = itanh(x) which is nonperiodic.
To analyze in more detail the problem of single magnetic scattering, we numerically calculate the value of ácosqñ as a function of the size parameter ka for dielectric ¢¢ = 0 scatterers for three different values of the real part of the relative magnetic permeability m¢, as is exhibited in Fig. 1 (for details, see Ref. ). In the small-particle limit (ka << 1), the scattering is isotropic (ácosqñ = 0) for nonmagnetic scatterers (m = 1). The appearance of usual Mie resonances in ácosqñ is observed as ka increases. This situation changes drastically when magnetic scatterers are considered. Notice that EM scattering by magnetic particles is characterized by the nonvanishing value of ácosqñ even in the small-particle limit, as it can be seen by m¢ = 2 and m¢ = 100 curves. Furthermore, we observe the presence of resonance phenomena for large values of m in ácosqñ even in the small-particle limit, which is absent in the nonmagnetic case. This is in agreement with our analytical result [see Eqs. (2) and comments that follow].
To further investigate EM scattering by magnetic particles in the small-particle limit, in Fig. 2 ácosqñ is exhibited as function of both the real (m¢) and imaginary (m¢¢) part of m for ka = 0.63. ácosqñ has an oscillatory dependence on m¢, confirming the presence of resonance phenomena. As m¢ increases, the real part of m also increases. A large value of m leads to the building up of standing waves inside the scatterer, which are responsible for the observed resonance phenomena. In addition, notice that ácosqñ can be negative for some
|Figure 1. The scattering anisotropy factor ácosqñ, plotted as a function of the size parameter ka for three values of the real part of the relative magnetic permeability: nonmagnetic case m¢ = 1 (dashed curve), m¢ = 2 (dotted curve) and m¢ = 100 (solid curve). The other parameters used are: ¢ = 1.4161, ¢¢ = 0 and m¢¢ = 0. [reprinted from Pinheiro, Martinez and Sampaio, Phys. Rev. Lett. 84, 1435 (2000), with permission from the American Physical Society].|
|Figure 2. ácosqñ plotted as a function of both the real m¢ (solid curve) and the imaginary m¢¢ (dashed curve) parts of m in the small-particle limit ka = 0.63 and in the large-particle limit ka = 63 (inset). The relative dielectric constant is = ¢ = 1.4161. [reprinted from Pinheiro, Martinez and Sampaio, Phys. Rev. Lett. 84, 1435 (2000), with permission from the American Physical Society].|
values of m, meaning a predominant backward scattering. This preferential backscattering in the small-particle limit, which is quite unusual in EM scattering by nonmagnetic particles, is due to the significant contribution given by magnetic dipole radiation in magnetic scattering. On the other hand, notice the absence of resonances when m has a pure imaginary value: ácosqñ is a very slowly varying function of m¢¢. We numerically obtain the limiting value ácosqñ ~ -0.3 when m¢¢ ® ¥. This means that a high level of magnetic dissipation of energy contributes only with a constant negative value to ácosqñ in the small-particle limit. In the inset of Fig. 2 ácosqñ is plotted as a function of m¢ and m¢¢ for a size parameter ka = 63, in the so-called large-particle limit. We observe the typical Mie resonances for the real magnetic permeability case, as expected, but the scattering is always preferentially in the forward direction, i.e., ácosqñ > 0.
III Multiple Scattering
An important quantity in multiple scattering is the mean length that the wavevector loose memory of the incidence direction. This length is given by the transport mean free path * = /(1 - ácosqñ), where = 1/(Fst) is the mean free path, F is the number of particle per volume and st is the total scattering cross section. The localization (Ioffe-Regel) parameter is defined as k* and strong localization is expect when k* ~ 1.
III.1 Localization parameter
In the following, using the above results concerning single magnetic scattering, we numerically calculate 1/(k*), which is proportional to the width of the backscattering cone, in a magnetic medium for both the small and large particle limits. We only consider the scalar aspect of light propagation. In Fig. 3 we show the dependence of 1/(k*) on both m¢ and m¢¢ in the small-particle limit for dielectric scatterers. Notice the global decrease of the localization parameter induced by the presence of magnetic (m ¹ 1) scatterers. Furthermore, while m¢¢ gives an almost constant contribution to the cone width, m¢ contributes to a fast oscillation. It is also important to point out the wide spectrum of values involved in this oscillation: the cone width can vary over three orders of magnitude when |m| varies from 1 to 100. This interesting pattern of the angular profile of the backscattering cone is essentially due to the fact that the scattering anisotropy can oscillate between positive and negative values as m varies.
Figure 3. The inverse of the localization parameter 1/(k*), proportional to the backscattering cone width, plotted as a function of both m¢ (solid curve) and m¢¢ (dashed curve) in the small-particle limit ka = 0.63 for = ¢ = 1.4161. The wavelength is l = 633 nm and the concentration of scatterers F = 0.01. [reprinted from Pinheiro, Martinez and Sampaio, Phys. Rev. Lett. 84, 1435 (2000), with permission from the American Physical Society].
The oscillatory behavior for real values of m is also present in the large-particle limit (see Fig. 4). Here, not only we observe a global decrease in the localization parameter, when compared with the nonmagnetic case, but also the fact that k * decreases as both m¢ and m¢¢ increases until it achieves a saturation minimum value.
|Figure 4. The inverse of the localization parameter 1/(k*) plotted as a function of both m¢ (solid curve) and m¢¢ (dashed curve) in the large-particle limit ka = 63 for = ¢ = 1.4161. The wavelength is l = 633 nm and the concentration of scatterers F = 0.01. [reprinted from Pinheiro, Martinez and Sampaio, Phys. Rev. Lett. 84, 1435 (2000), with permission from the American Physical Society].|
III.2 Energy transport velocity and diffusion constant
The unusual resonant behavior exhibited by EM scattering by magnetic particles in the small-particle limit discussed in Section 2 has also important implications in the dynamics of the diffusion process. In fact, the diffusion constant D = vE*/3 of EM propagation is strongly affected in the presence of magnetic scatterers, not only through the transport mean free path *, as discussed in Section 3.1, but also through the "transport velocity" vE. The fact that single scattering by magnetic particles exhibits a characteristic resonant behavior, even in the small-particle limit, will introduce an extra time delay in EM propagation, causing a decrease of vE if compared to the nonmagnetic case. To solve a discrepancy between stationary and dynamic experiments, van Albada et al.,[6, 7] have developed a scalar theory using the low-density approximation, which takes into account a renormalization mechanism of D for classical waves. According to this theory, the transport velocity vE for spherical dielectric Mie scatterers is given by:
where c0 is the velocity of propagation in the vacuum, f is the packing fraction and vp = c0/(1+fC)1/2 is the phase velocity. The phase angles an(ka) and bn(ka) are defined in terms of the Mie scattering coefficients an(ka) and bn(ka). The parameter C(ka) is defined according to:
|Figure 5. The energy transport velocity vE (normalized by c0), plotted as a function of ka in the small-particle limit for three values of the (real) magnetic permeability: m = 1 (dashed curve), m = 103 (bold solid curve) and m = 105 (solid curve). The packing fraction is f = 0.36 and = ¢ = 7.29. [reprinted from Pinheiro, Martinez and Sampaio, Phys. Rev. Lett. ?? , with permission from the American Physical Society].|
From Eqs. (2) and (3), we have analytically calculated vE for ferromagnetic Mie scatterers in the small-particle limit, where resonance phenomena are present in a singular and unusual way. In Fig. 5, vE is exhibited as a function of the size parameter ka in the small-particle limit for three different values of m for a packing fraction f = 0.36, the same used in experiments with nonmagnetic scatterers TiO2. For comparison, the curve m = 1 (nonmagnetic case), calculated using the well-known expressions for a1 and b1 in the Rayleigh limit,[13, 14, 15] is also shown. Notice that the absence of resonances in the nonmagnetic case (m = 1) produces an almost constant value of vE. This situation changes drastically when magnetic scatterers are taken into account. In this case, vE exhibits sharp drops even in the small-particle limit, in clear contrast to the nonmagnetic case where these drops only occur in the intermediate regime (ka ~ 1) and in the large-particle limit (ka 1). In addition, notice that vE can even vanish for some values of ka. As pointed out by van Albada et al.,[6, 7] these sharp drops in vE are essentially due to resonances in the scatterers cross-section, which cause the building up of standing waves inside the scatterers, introducing an extra time delay in EM propagation.
The influence of the presence of magnetic scatterers on the diffusion constant D in the small-particle (ka = 0.01) limit can be seen in Fig. 6, where D, vE and * are exhibited as a function of the relative magnetic permeability m. Here one can notice a striking consequence of the introduction of magnetic scatterers in the dynamic properties of EM propagation through disordered media: the presence of oscillations in the diffusion constant D. For some values of m, in particular, one can observe the situation of vanishing diffusion (D = 0) of electromagnetic propagation. In addition, it is interesting to point out the damping in the diffusion constant oscillations induced by the increase of the values of m.
|Figure 6. The transport velocity vE (solid curve), transport mean free path * (dashed curve) and diffusion constant D (bold solid curve), plotted as a function of the (real) ( for dielectric ferromagnetic scatterers in the small-particle limit (ka = 0.01). In order to exhibit dimensionless quantities vE, * and D are multiplied by 1/c0, 3f/4a and 3f/4ac0, respectively. For comparison, vE is multiplied by the factor 107. The other parameters used are: ¢ = 7.29 and f = 0.36. [reprinted from Pinheiro, Martinez and Sampaio, Phys. Rev. Lett. ?? , with permission from the American Physical Society].|
IV Experimental Aspects
Let us now briefly discuss possible situations where our results could be experimentally verified. Most of the soft ferromagnets and ferrimagnets display high magnetic permeability tensor with m being strongly dependent on an external applied magnetic field Hext and on the temperature T, i.e., m = m(Hext, T). These facts make all the mentioned properties concerning both single and multiple scattering of EM waves by magnetic particles tunable mechanisms: they can be controlled by the variation of an external parameter, such as Hext or T. For instance, soft ferrites present large values of m at microwave frequencies generally below 100 of MHz and without loss, which in turn can be easily controlled by the application of an external dc magnetic bias field. For the majority of ferrites, in high frequencies, m becomes very small, which turns the utilization of this tunable scattering mechanism fairly difficult at optical frequencies. This kind of restriction to the magnetic permeability have also been employed previously in other contexts involving the interaction between electromagnetic radiation and magnetic materials, such as reported by Sigalas et al., who have recently discussed the possibility of designing photonic crystals with soft magnetic materials.
The magnetic permeability of the vast majority of ferrites typically follows the Curie-Weiss law above the critical temperature, i.e., m µ (T - Tc)-1. Thus, a precise tuning of m can be performed in this temperature range. As a result, the utilization of magnetic scatterers following the Curie-Weiss susceptibility law offers the possibility to tune the value of both the localization parameter k* and the diffusion constant D by varying the sample temperature. In the large-particle limit, in particular, this fact could be explored in alternative experimental investigations of Anderson localization of EM waves since, as the temperature approaches the Curie-Weiss critical temperature, one can set k* to a minimum value (see Sec. III).
Finally, we suggest that the possibility of controlling the oscillations in the diffusion constant (Sec. III) offers the theoretical basis for the implementation of a magneto-optical device, able to control the transmission of EM waves through disordered magnetic media. Such a device would operate as an optical switch controlled by either the temperature or by an external dc magnetic field (for device applications exploring the transport properties of EM waves through random media see, for instance, Ref. ).
Summarizing, we have studied the influence of magnetic scatterers on single and multiple electromagnetic scattering, and their implications in the coherent backscattering effect, localization parameter, energy transport velocity and diffusion constant. We have shown that, in the small-particle limit, not only single scattering exhibits a forward-backward spatial asymmetry, but also an unusual resonance effect. In multiple scattering context, these unusual aspects concerning single magnetic scattering lead to several interesting consequences, such as the global decrease of the localization parameter, sharp drops in the energy transport velocity even in the small-particle limit and oscillations in the diffusion constant. Furthermore, exploring the fact that the magnetic permeability of magnetic scatterers is strongly dependent on an external magnetic field Hext and typically follows the Curie-Weiss law, we have suggested that both the localization parameter and the diffusion constant can be tuned by varying either Hext or the temperature.
We thank I. S. Oliveira and L. G. Guimarães for fruitful discussions. This work was supported by the Brazilian agencies CAPES and CNPq.
 G. L .J. A. Rikken and B. A. van Tiggelen, Nature 381, 54 (1996). [ Links ]
 B. A. van Tiggelen, R. Maynard and A. Heiderich Phys. Rev. Lett. 77, 639 (1996); [ Links ]H. Stark and T. C. Lubensky, ibid. 77, 2229 (1996).
 D. S. Wiersma, P. Bartolini, A. Lagendijk and R. Righini, Nature 390, 671 (1997). [ Links ]
 B. A. van Tiggelen and A. Lagendjik, Phys. Rep. 270, 143 (1996). [ Links ]
 M. P. van Albada, B. A. van Tiggelen, A. Lagendijk and A. Tip, Phys. Rev. Lett. 66, 3132 (1991); [ Links ]
 B. A. van Tiggelen, A. Lagendijk, M. P. van Albada, and A. Tip, Phys. Rev. B 45, 12 233 (1992). [ Links ]
 R. Lenke and G. Maret, Phys. Scripta 49, 605 (1993). [ Links ]
 A. S. Martinez and R. Maynard, Phys. Rev. B 50, 3714 (1994). [ Links ]
 F. A. Pinheiro, A. S. Martinez and L. C. Sampaio, Phys. Rev. Lett. 84, 1435 (2000). [ Links ]
 F. A. Pinheiro, A. S. Martinez and L. C. Sampaio, Electromagnetic Scattering by Small Magnetic Particles (to appear in: Journal of Magnetism and Magnetic Materials). [ Links ]
 F. A. Pinheiro, A. S. Martinez and L. C. Sampaio, Vanishing of Energy Transport Velocity and Diffusion Constant Electromagnetic Waves in Disordered Magnetic Media (to appear in: Phys. Rev. Lett.). [ Links ]
 H. C. van de Hulst, Light Scattering by Small Particles (Dover, New York, 1980). [ Links ]
 M. Kerker, The Scattering of Light and Other Electromagnetic Radiation (Academic Press, New York, 1969). [ Links ]
 C. F. Bohren and D. R. Huffman, Absorption and Scattering of Light by Small Particles (John Wiley, New York, 1983). [ Links ]
 C. W. Chen, Magnetism and Metallurgy of Soft Magnetic Materials (Dover, New York, 1986). [ Links ]
 M. Kerker, D. S. Wang and C. L. Giles, J. Opt. Soc. Am. 73, 765 (1983). [ Links ]
 A. Z. Genack and N. Garcia, J. Opt. Soc. Am. B 10, 408 (1993). [ Links ]
 M. M. Sigalas, C. M. Soukoulis, R. Biswas, and K. M. Ho, Phys. Rev. B 56, 959 (1997). [ Links ]