## Abstract

Magneto-optical spectra (Faraday effect) for nanometer-scale particles and collections of particles are calculated using a modification of the discrete dipole approximation (DDA). The approximation is used as a finite-element method for non-spherical particles whose dimensions are on the order of or smaller than the incident light wavelength, λ. Also, we use the approximation to calculate scattering from arrangements of spherical nanoparticles with diameters ≪λ. We propose that for scattering from subwavelength magnetic particles, the specific Faraday rotation should be defined as the difference in optical extinction for left- and right-circularly polarized light. We apply the model to calculations of Faraday rotation from magnetite nanowires as well as a binary (two-component) nanoparticle arrangement. Enhancements in Faraday rotation are predicted for composites containing both noble metal and ferrite nanoparticles.

©2006 Optical Society of America

## 1. Introduction

Electromagnetic scattering from nanometer-scale objects is currently being investigated both experimentally and theoretically for the understanding and prediction of novel near-field effects. This is particularly relevant due to the current research on nanomaterials and the use of nanometer-scale optical probes like near-field scanning optical microscopy (NSOM). Typically, near-field optical studies are performed on noble-metal nanostructures where resonances from surface plasmons play an important role in near-field coupling and propagation effects. Previous studies have focused on surface-enhanced Raman scattering (SERS)[1, 2], anomalous transmission through nanowire gratings[3] and plasmonics[4–7]. In contrast, little attention has been focused on other potentially interesting nanometer-scale optical effects. In this paper, we examine the magneto-optical scattering from nanometer-scale structures using a discrete dipole approximation (DDA). In the simplest form, the approximation replaces individual nanoparticles with a single, electromagnetic dipole radiator. The optical response of a collection of particles is calculated self consistently by calculating the response of each dipole to the incident field plus the scattered field from all the other dipoles.[8, 9] Used in this way, the approximation is often called the coupled dipole approximation (CDA). In addition, the approximation is used as a finite element computation technique where a complicated, non spherical object (target) is discretized into an array of interacting dipole radiators[10, 11]. The electric field, or alternatively, the scattering matrix is then calculated by summing the contribution from each dipole element.

Earlier, we reported on the magneto-optical response (Faraday rotation) from collections of Fe_{3}O_{4} (magnetite) nanoparticles[12]. We found that the spectral Faraday rotation was dependent on the nanoparticle spacing. We were able to qualitatively model this behavior using the DDA formulated in the left-right circular polarization basis corresponding to the normal propagating modes for a magneto-optical medium. Here, we extend those calculations with a more general formulation of the DDA using a non-diagonal dielectric tensor appropriate to gyrotropic media and applicable to arbitrary magnetization direction. In addition, we state a general expression for effective Faraday rotation in terms of scattering matrix elements derived from the calculated polarizabilities. We demonstrate the use of this model by calculating the magneto-optical response of extended nanowire structures and nanoparticle assemblies containing two types of nanoparticles, one a magnetic oxide and one a noble metal.

## 2. Theory

#### 2.1 Basic equations

In the DDA approximation, either the spherical nanoparticle, with diameter 2*a*≪*λ*, or a similarly small element of an extended structure with dimensions on the order*λ* is represented by an oscillating dipole with polarizability tensor **α**
*j* which is determined by the material properties. The dipole moment of the element located at **r**
_{j}
is **P**
_{j}
=α*j*
**E**
_{loc,j} where the local field **E**
_{loc,j} is the incident field **E**
_{inc,j} plus the field radiated by all the other dipoles,

The incident field is given by **E**
_{inc,j}=**E**
_{0} exp[*i*
**κ·r**
_{j}
-*ωt*] where *ω* is the frequency of radiation and κ is the wavevector with magnitude *κ=ω/c*.

Substituting the expression for the dipole moment and rearranging yields,

where the radiation from a dipole has been written **E**
_{dp,k}=**A**
_{jk}
**P**
_{k}
. Using the well-known expression for radiation of a dipole located at position
_{k}
r evaluated at point
_{j}
r [13, 14],

where **n**
_{jk}
is a unit vector pointing from *r*_{k}
to *r*_{j}
, and *r*_{jk}
=|**r**
_{j}
-**r**
_{k}
|. Equations (2) and (3) are simplified by defining **A**
_{jj}
=**α**1-
_{j}
so that the equations can be written as a single matrix equation[14]

For *N* dipole elements, Ã is a 3*N*×3*N* matrix which describes the interactions between dipoles, _{inc}
**Ẽ** is a 3*N*×1 vector, and **P̃** is a 3*N*×1 vector of unknown dipole moments. Notice that **A**
_{jk}
, defined in Eq. (3), is itself a 3×3 matrix which becomes the (*j, k*) element of Ã.

Key to the implementation of the DDA is the connection between the dielectric properties of the material given by the dielectric tensor ε and the polarizability of the individual elements, **α**
_{j}
. The optical interactions with a material which exhibits magneto-optical effects (a gyrotropic material) are described by its dielectric tensor,

where *Q*≡*ε*_{xy}
is the magneto-optical Voigt parameter and *m*_{x}*, m*_{y}
and *m*_{z}
are the direction cosines of the magnetization vector. The induced (uniform) polarization of a dielectric sphere valid for arbitrary anisotropy is given by[15, 16],

where 1 is the 3×3 identity matrix and *a* is the radius of the sphere This allows the identification of the polarizability as[15]

where (*ε*+21)^{-1} is understood as matrix inversion. Draine et al. have established corrections to the basic formula which include a radiation reaction correction for finite wavelength[17]. From Eq. (7) a Clausius-Mossotti formula can be derived by summing the contributions from individual spherical polarizable elements[16]. Corrections which insure that an infinite lattice of polarizable points has the same dispersion relation as the macroscopic material are also given by Draine[17]. A note of caution: Eq. (7) is exact for infinite wavelength only for an isolated sphere or a dipole on an infinite lattice. For scatters with high aspect ratios (or otherwise arbitrary shapes) Eq. (7) should be modified to account for local field corrections[18]. Here, we wish to highlight the differences between the implementation of the DDA for isotropic and gyrotropic media, so we use the most basic form of the polarizability, Eq. (7).

In our magneto-optical experiments to date we have used the polar geometry in which the magnetic field is parallel (or antiparallel) to the direction of light propagation, in other words, *m*_{z}
=1 and *m*_{x}
=*m*_{y}
=0 in Eq. (5). The equations are valid for arbitrary magnetization; the present implementation of the DDA will accommodate non-zero *m*_{x}
or *m*_{y}
without difficulty. We point out that the connection between magnetic field and optical field is through the magnetization-dependent dielectric tensor. In fact, at optical frequencies, the magnetic permeability is equal to unity and the effects of the magnetic field are characterized by the dielectric tensor through its dependence on the magnetization (see, for example, Ref. [19]).

#### 2.2 Scattering matrix elements

The macroscopic magneto-optical effects, Kerr rotation for reflection and Faraday rotation for transmission, are defined in terms of Fresnel coefficients[20]. These familiar Fresnel equations are approximations which apply to optical interactions with surfaces larger than many square wavelengths[21]. More generally, and perhaps more fundamentally, optical properties are manifestations of absorption and reradiation of light from individual elements of the material. For this reason, a slightly different definition is needed for magneto-optical effects arising from materials whose dimensions are on the order of or smaller than an optical wavelength. We define the complex Faraday rotation angle as the difference in the scattering extinction coefficients for left and right circularly polarized light. The Faraday rotation, *θ* and ellipticity, *η* (per unit length), are given by

The extinction coefficients are calculated for a nanostructure (or collection of nanostructures) and the observable Faraday rotation is interpreted as being the result of a slab of a given thickness containing a given volume density of nanostructures.

In general, the calculation of the extinction coefficients requires the calculation of dipole moments from Eq. (4) for two orthogonal polarizations of incident electric field. The calculation of the scattering matrix elements using the DDA is given explicitly by Draine[14]. In our case, we take the both the incident and scattered fields to be the ẑ direction (forward scattering), which matches our experimental geometry. Using Draine’s notation, the scattering matrix elements are

where ê is a unit vector for the polarization and *m* and *m*′ denote the incident and scattered fields, respectively. ${\mathbf{P}}_{j}^{\mathit{\left(}m\mathit{\right)}}$
is the calculated dipole moment of the *j*
^{th} particle for incident polarization ê
_{m}
. Since both the incident and scattered fields are in the ẑ direction, the polarization unit vectors are simply x̂ and ŷ. Basically, we calculate the dipole moments for incident *x* polarization and evaluate the scattering into the *y* polarization and show that this is equivalent to the difference in extinction coefficients for left and right circularly polarized light. The Faraday rotation is calculated from the *f*_{xy}
matrix element,

Substituting for *f*_{xy}
from Eq. (9),

We now transform into the left-right circular polarization basis to show that this is indeed equivalent to Eq. (8). The usual transformation from (*x,y,z*) to (*r,l,z*), where r and l stand for right and left respectively, is given by ê
_{r}
=1/√2(ê*x*-iê
_{y}
) and êl=1/√2(ê
_{x}
+iê
_{y}
). In this basis the *x* polarization is written

and the *y* unit vector is

Substituting this into Eq. (11) yields,

which is equivalent to difference in extinction coefficients for left and right circularly polarized light, *C*_{ext,l}*-C*_{ext,r}
.

The calculation proceeds as follows. For a given geometry, a lattice is established with individual polarizabilities given by Eq. (7); values of the dielectric tensor are taken from the literature. Equation (4) is solved with incident *x*-polarized electric field which yields **P**
^{(x)}
*j*. The Faraday rotation is calculated using Eqs. (11). Calculations involving few dipoles (<30) were performed using Mathematica[22]. Calculations involving more than 30 dipoles (up to several hundred) were performed using MATLAB using the available complex-conjugate gradient (CCG) matrix solver[23].

## 3. Results

We first apply the DDA calculations to magnetite nanowires. The model for the nanowires is similar to that used in Ref. [10] and is shown in Fig. 1. The diagonal and off diagonal components of the dielectric tensor are taken from the literature[24, 25]. The Faraday rotation is calculated for nanowires with 2 nm×2 nm cross section and lengths from 2 nm to 100 nm.

The results of these calculations are shown in Fig. 2 for incident *x* polarization and nanowires aligned on the *x* axis and *z* axis. The magneto-optics calculations assume that the externally-applied magnetic field is large enough to saturate the magnetization of the material. In order to highlight differences in nanowire orientation, we assume that the magnetization is in the *z* direction in both Fig. 2 (a) and (b).

The absolute Faraday rotation is dependent on the volume of each scatterer as well as the number density of scatterers and thickness of the sample. We have normalized the calculations in Fig. 2 so that the Faraday rotation at the positive peak (~500 nm) is unity. The spectrum of the 2×2×2 element particle is the same as a single nanoparticle. Little change in the spectral Faraday rotation is noted for nanowires when the nanowire is aligned with polarization of the incident field. However, the Faraday rotation spectrum changes dramatically when the nanowire is aligned in the direction of light propagation (z direction). The overall effect is a blue shift of the spectrum with increasing nanowire length.

We also performed calculations on a noble-metal/magnetic nanoparticle system to illustrate possible near-field effects. There are some experimental reports of enhancements in Faraday and Kerr effect for magnetic materials physically located near noble-metal nanoparticles [26–28] or thin films[29–31]. The electromagnetic field is concentrated near the noble metal surface when the optical frequency is resonant with surface plasmon. It is theorized that this is the origin of the observed enhancements and may be considered to be the magneto-optical counterpart of surface-enhanced Raman scattering[29, 32].

Specifically, we calculate the response of assemblies of silver and cobalt ferrite (CoFe_{2}O_{4}) nanoparticles separated by a few nanometers. This system was chosen because of the recent experimental report of magneto-optic enhancements in coupled silver-cobalt ferrite nanoparticles[27]. However, complete *spectral* magneto-optical data is not available, only Faraday rotation at selected laser wavelengths.

For these calculations, each particle is represented by a single dipole. Dielectric tensor data is taken from the literature for CoFe_{2}O_{4}.[33, 34] and silver[35] We used the experimental values for the particle radii, 7 nm for CoFe_{2}O_{4} and 3 nm for silver[27]. Our results are shown in Fig. 3 for two-particle systems. A background dielectric constant of 1.7 was chosen so that the plasmon peak of silver was shifted to ~421 nm as the data in Ref. [27] indicate. In our calculations, the presence of the near-by noble metal particle results in enhancements in the Faraday rotation near the plasmon peak and is dependent on the particle geometry.

Figure 4 shows the optical extinction and magneto-optical response of a three-dimensional binary array: CoFe_{2}O_{4} and silver arranged in a rock-salt structure. Binary assemblies of noble metal and magnetic particles have recently been demonstrated. These assemblies can have a variety of mesoscopic arrangements isostructural with rock salt, Cu_{3}Au and AlB_{2}, for example[36]. The optical extinction spectrum is dependent, of course, on interparticle spacing, even splitting into two peaks at the smallest separation. Previous calculations have shown line broadening, splitting and even line narrowing in regular arrays of noble metal nanoparticles[8]. In our calculations, the primary effect in the magneto-optical spectrum is at or near the surface plasmon wavelength of the noble-metal particles and results in an increase in the magnitude of the Faraday rotation of up to a factor of five. The spectrum is also dependent on the particle arrangement. For example, in the calculations presented here, the ferrite particles are located on the unit cube corners; slightly different spectra are obtained if the positions of the ferrite and silver particles are interchanged.

## 4. Discussion

These calculations provide a tool for modeling the magneto-optical scattering from non-spherical structures as well as provide a framework for the basic, qualitative understanding of the near-field coupling effects on magneto-optical scattering. Incidentally, none of these effects are predicted by any effective-medium type model or any mean-field theory which does not include near-field interactions. The calculations for the nanowires predict shape and orientation effects. The effects are most evident for the wires aligned with the *z* axis. This suggests that the effects are due to an anisotropy induced by the wire. The incident *x*-polarized optical wave consists of equal parts left- and right-circular polarizations. For bulk materials or nanowires aligned in the *x-y* plane, these are the normal modes of propagation. However, for the wires aligned with the *z* axis, these modes no longer propagate without change of polarization state; there is a depolarization of the light. Note that magneto-optical effects arise due differences in the real or imaginary parts of the index of refraction for circularly polarized light Anisotropy in the *x-y* plane could effect *linear* but not *circular* birefringence or dichroism.

For the binary nanoparticle arrangements, the presence of the noble-metal and subsequent concentration of the electromagnetic field perturbs the polarization of the magnetic nanoparticle and, hence, the reradiation of the optical field. As would be expected, the most profound changes are observed near resonances of the silver particle (or particle arrangement). We predict that such effects would be observed for particle spacings on the order of 20 nm or less. In Ref. [27], the silver and cobalt ferrite particles were attached forming dimers (two-particle assemblies). The optical absorption shows a peak at ~421 nm, but the primary differences between the measured Faraday rotation for cobalt ferrite alone and cobalt ferrite with silver particles attached occurs at much longer wavelengths. Other reports of magneto-optical enhancements resulting from composites with noble-metal nanoparticles contain Faraday rotation data at even fewer wavelengths. Without more complete spectral Faraday rotation, we can only speculate on the absence of an observable effect near the plasmon resonance.

## 5. Conclusion

We have developed a model for magneto-optical scattering from nanometer-scale structures based on a modification of the discrete dipole approximation. The model can be used as finite element method for the calculation of the magneto-optical response of non-spherical nanostructures with dimensions approaching the wavelength of light. The model can also be used to calculate magneto-optical scattering from collections of nanoparticles, where each nanoparticle is represented by a radiating dipole. The particle collections need not be homogeneous; in fact, we apply the model to binary nanoparticle systems which consist of a noble-metal and magnetic oxide nanoparticles. We calculate enhancements of the Faraday rotation near the plasmon resonances of the metal nanoparticles. These enhancements are evident for interparticle spacings on the order of 20 nm or less.

## Acknowledgments

We gratefully acknowledge the support of DARPA through grant HR0011-04-1-0029 and the Louisiana State Board of Regents.

## References and links

**1. **N. Felidj, J. Aubard, G. Levi, J. R. Krenn, M. Salerno, G. Schider, B. Lamprecht, A. Leitner, and F. R. Aussenegg, “Controlling the optical response of regular arrays of gold particles for surface enhanced Raman scattering,” Phys. Rev. B. **65**, 075419 (2002). [CrossRef]

**2. **W.-H. Yang, G. C. Schatz, and R. P. Van Duyne, “Discrete dipole approximation for calculating extinction and Raman intensities for small particles with arbitrary shapes,” J. Chem. Phys. **103**, 869 (1995). [CrossRef]

**3. **F. J. Garcia-Vidal and L. Martin-Moreno, “Transmission and focusing of light in one-dimensional periodically nanostructured materials,” Phys. Rev. B. **66**, 155412 (2002). [CrossRef]

**4. **M. L. Brongersma, J. W. Hartman, and H. A. Atwater, “Electromagnetic energy transfer and switching in nanoparticle chain arrays below the diffraction limit,” Phys. Rev. B. **62**, R16356 (2000). [CrossRef]

**5. **S. A. Maier, M. L. Brongersma, P. G. Kik, S. Meltzer, A. A. G. Requicha, and H. A. Atwater, “Plasmonics - A route to nanoscale optical devices,” Adv. Mater. **13**, 1501 (2001). [CrossRef]

**6. **S. A. Maier, P. G. Kik, M. L. Brongersma, H. A. Atwater, S. Meltzer, A. A. G. Requicha, and B. K. Koel
, “Observation of coupled plasmon-polariton modes of plasmon waveguides for electromagnetic energy transport below the diffraction limit,” in *Mat. Res. Soc. Symp. Proc*., vol. 722,
R. B. Wehrspohn, S. Noda, C. Soukoulis, and R. Marz, Eds. (Warrendale, PA, MRS, 2002), p. L6.2.

**7. **V. M. Shalaev, V. A. Podolskiy, and A. K. Sarychev
, “Plasmonic nanophotonics: Manipulating light and sensing molecules,” in *Proceedings of SPIE vol. 4806, Complex Mediums III: Beyond Linear Isotropic Dielectrics*,
A. Lakhtakia, G. Dewar, and M. W. McCall, Eds. (Bellingham, WA, SPIE, 2002), p. 32.

**8. **L. L. Zhao, K. L. Kelly, and G. C. Schatz, “The extinction spectra of silver nanoparticle arrays: Influence of array structure on plasmon resonance wavelength and width,” J. Phys. Chem. B **107**, 7343 (2003). [CrossRef]

**9. **A. A. Lazarides and G. C. Schatz, “DNA-linked metal nanosphere materials: Structural basis for the optical properties,” J. Phys. Chem. B **104**, 460 (2000). [CrossRef]

**10. **V. A. Podolskiy, A. K. Sarychev, and V. M. Shalaev, “Plasmon modes and negative refraction in metal nanowire composites,” Opt. Express **11**, 735 (2003). [CrossRef] [PubMed]

**11. **V. A. Podolskiy, A. K. Sarychev, E. E. Narimanov, and V. M. Shalaev, “Resonant light interaction with plasmonic nanowire systems,” J. Opt. A: Pure Appl. Opt. **7**, S32 (2005). [CrossRef]

**12. **D. A. Smith, Y. A. Barnakov, B. L. Scott, S. A. White, and K. L. Stokes, “Magneto-optical spectra of closely-spaced magnetite nanoparticles,” J. Appl. Phys. **97**, 10M504 (2005). [CrossRef]

**13. **J. D. Jackson, *Classical Electrodynamics*, 3rd ed. (New York, John Wiley and Sons, Inc., 1999).

**14. **B. T. Draine, “The discrete dipole approximation and its application to interstellar graphite grains,” Astrophys. J. **333**, 848 (1988). [CrossRef]

**15. **A. H. Sihvola, “Dielectric polarizability of a sphere with arbitrary anisotropy,” Opt. Lett. **19**, 430 (1994). [CrossRef] [PubMed]

**16. **K. Hinsen, A. Bratz, and B. U. Felderhof, “Anisotropic dielectric tensor and the Hall effect in a suspension of spheres,” Phys. Rev. B. **23**, 14995 (1992). [CrossRef]

**17. **B. T. Draine and J. Goodman, “Beyond Clausius-Mossotti: Wave propagation on a polarizable point lattice and the discrete dipole approximation,” Astrophys. J. **405**, 685 (1993). [CrossRef]

**18. **A. Rahmani, P. C. Chaumet, and G. W. Bryant, “On the importance of local-field corrections for polarizable particles on a finite lattice: Application to the discrete dipole approximation,” Astrophys. J. **607**, 873 (2004). [CrossRef]

**19. **P. S. Pershan, “Magneto-optical effects,” J. Appl. Phys. **38**, 1482 (1967). [CrossRef]

**20. **A. K. Zvezdin and V. A. Kotov, *Modern Magnetooptics and Magnetooptical Materials* (London, Institute of Physics Publishing, 1997). [CrossRef]

**21. **C. F. Bohren and D. R. Huffman, *Absorption and Scattering of Light by Small Particles* (New York, John Wiley and Sons, 1983).

**22. **Wolfram Research, Inc., *Mathematica, Version 5.1*, Champaign, IL (2004).

**23. **The MathWorks Inc, *MATLAB version 6.5.1*, Natick, Massachusetts (2003).

**24. **A. Schlegel, S. F. Alvarado, and P. Wachter, “Optical properties of magnetite (Fe_{3}O_{4}),” J. Phys. C: Solid State Phys. **12**, 1157 (1979). [CrossRef]

**25. **X. Zhang, J. Schoenes, and P. Wachter, “Kerr-effect and dielectric tensor elements of magnetite (Fe_{3}O_{4}) between 0.5 and 4.3 eV,” Solid State Commun. **39**, 189 (1981). [CrossRef]

**26. **M. R. Pufall, A. Berger, and S. Schultz, “Measurement of the scattered light magneto-optical Kerr effect from plasmon-resonant Ag particles near a magnetic film,” J. Appl. Phys. **81**, 5689 (1997). [CrossRef]

**27. **Y. Li, Q. Zhang, A. V. Nurmikko, and S. Sun, “Enhanced magneto-optical response in dumbbell-like Ag-CoFe_{2}O_{4} nanoparticle pairs,” Nano Lett. **5**, 1689 (2005). [CrossRef] [PubMed]

**28. **G. Shemer and G. Markovich, “Enhancement of magneto-optical effects in magnetite near gold surfaces,” J. Phys. Chem. B **106**, 9195 (2002). [CrossRef]

**29. **V. I. Safarov, V. A. Kosobukin, C. Hermann, G. Lampel, and J. Peretti, “Magneto-optical effects enhanced by surface plasmons in metallic multilayer films,” Phys. Rev. Lett. **73**, 3584 (1994). [CrossRef] [PubMed]

**30. **T. Katayama, Y. Suzuki, H. Awano, Y. Nishihara, and N. Koshizuka, “Enhancement of the magneto-optical Kerr rotation in Fe/Cu bilayered films,” Phys. Rev. Lett. **60**, 1426 (1988). [CrossRef] [PubMed]

**31. **C. Hermann, V. A. Kosobukin, G. Lampel, J. Peretti, V. I. Safarov, and P. Bertrand, “Surface-enhanced magneto-optics in metallic multilayer films,” Phys. Rev. B. **64**, 235422 (2001). [CrossRef]

**32. **V. A. Kosobukin, “Magneto-optics via the near field,” Surf. Science **406**, 32 (1998). [CrossRef]

**33. **W. F. J. Fontijn, P. J. van der Zaag, and L. F. Feiner, “A consistent interpretation of the magneto-optical spectra of spinel-type ferrites,” J. Appl. Phys. **85**, 5100 (1999). [CrossRef]

**34. **W. D. Martens, W. L. Peeters, H. M. van Noort, and M. Erman, “Optical, magneto-optical and Mossbauer spectroscopy on Co^{3+} substituted cobalt ferrite,” J. Phys. Chem. Solids **46**, 411 (1985). [CrossRef]

**35. **P. B. Johnson and R. W. Christy, “Optical constants of noble metals,” Phys. Rev. B **6**, 4370 (1972). [CrossRef]

**36. **E. V. Shevchenko, D. V. Talapin, C. B. Murray, and S. O’Brien, “Structural characterization of self-assembled multifunctional binary nanoparticle superlattices,” J. Am. Chem. Soc. **128**, 3620 (2006). [CrossRef] [PubMed]