Abstract
We model the response of nanoscale Ag prolate spheroids to an external uniform static electric field using simulations based on the discrete dipole approximation, in which the spheroid is represented as a collection of polarizable subunits. We compare the results of simulations that employ subunit polarizabilities derived from the Clausius–Mossotti relation with those of simulations that employ polarizabilities that include a local environmental correction for subunits near the spheroid’s surface [Rahmani et al. Opt Lett 27: 2118 (2002)]. The simulations that employ corrected polarizabilities give predictions in very good agreement with exact results obtained by solving Laplace’s equation. In contrast, simulations that employ uncorrected Clausius–Mossotti polarizabilities substantially underestimate the extent of the electric field “hot spot” near the spheroid’s sharp tip, and give predictions for the field enhancement factor near the tip that are 30 to 50% too small.
Keywords:
Metallic nanoparticles; Optical properties; SimulationNano Express
The electrical and optical properties of noble metal nanoparticles have attracted considerable scientific interest for many decades. Over a century ago, for example, Mie [1]—building on even earlier work by Lorenz [2] and possibly others—attributed the colors of colloidal suspensions of Au nanoparticles [3] to the nanoparticles’ visiblewavelength optical scattering properties. Interest in the optical properties of noble metal nanoparticles has risen dramatically in recent years with the recognition that these properties, if understood in sufficient detail, can be harnessed to create nanoscale photonic devices and sensors.
The discrete dipole approximation [4,5] (also called the coupled dipole approximation) is one of several numerical methods that have been developed to simulate the response of a small particle to an incident electromagnetic (EM) field. In simulations based on the discrete dipole approximation (DDA), a nanoparticle is modeled as a regular (typically cubic) lattice of polarizable subunits. The incident EM field induces dipole moments in each subunit; these dipole moments in turn generate local fields that further polarize nearby subunits. Once the subunits’ induced dipole moments are mutually selfconsistent, the electromagnetic and optical properties of the dipole lattice are taken to mimic those of the real nanoparticle. The assumption that only dipolar interactions among subunits and between the subunits and the external field need be considered, an assumption that is implicit in DDAbased simulations, is generally thought to be a reasonable one provided that the subunits are small enough so that the electric field is nearly constant across an individual subunit; this assumption is frequently tested by comparing the results obtained from simulations at two or more levels of discretization.
The connection between the lattice of dipoles and the real nanoparticle is made through the choice of the polarizability tensor α of the polarizable subunits. In the original formulation [4] of the DDA approach, α was assumed to be an isotropic, diagonal tensor defined by the Clausius–Mossotti (CM) relation
where ρ is the number density of the polarizable subunits and ε is the nanoparticle’s dielectric constant; this relation is exact for an infinite cubic lattice of subunits in a zerofrequency external electric field [6]. For finite (nonzero) frequency external EM fields, a radiative reaction [6] correction to the zerofrequency polarizability tensor defined by Eq. 1 ensures that the optical theorem holds for the dipole lattice [5]. Other finite frequency corrections to the static polarizability given by Eq. 1 can be derived from an analysis of the dispersion relation for electromagnetic waves propagating along a lattice of polarizable points [7].
Real nanoparticles, of course, have surfaces, and hence cannot be represented as infinite lattices; consequently, the use of polarizabilities defined by Eq. 1 in DDAbased simulations of nanoparticles represents an additional approximation [8], one which persists even when the subunits are very small, which is not remediated by radiative reaction corrections or other finite frequency corrections, and which seems to be especially severe for materials whose dielectric constant has a large imaginary component [9]. Recent work [9,10] suggests that the use of subunit polarizabilities that properly account for the anisotropic local environment of DDA subunits near surfaces can increase substantially the accuracy with which highly averaged farfield quantities, such as absorption and scattering crosssections, can be computed using DDAbased methods. In this letter, we employ these corrected polarizabilities in DDAbased simulations of nanoscale Ag prolate spheroids in homogeneous static electric fields; we find that the new polarizabilities, which include a local environmental correction (LEC) to the CM polarizabilities, also substantially improve the description of spatially resolved nearfield quantities, such as localized electric field enhancement factors, computed in these simulations.
We begin by summarizing some exact results obtained by solving Laplace’s equation for a homogeneous prolate spheroid in a uniform static external field [11,12]; these are the benchmarks against which we assess the DDAbased simulations. We consider a prolate spheroid with major semiaxis c (henceforth assumed to coincide with the spacefixed z axis) and minor semiaxis a. The surface of the spheroid is one member of a family of confocal surfaces defined by the parameter ξ. These surfaces satisfy the equation
the surface of the spheroid corresponds to ξ = 0. If such a spheroid, with dielectric constant ε, is immersed in a medium with dielectric constant ε_{m} and exposed to a uniform static electric field parallel to the spacefixed z axis, the electrical potential at any point outside the spheroid is given by
where s = (ε − ε_{m})/ε_{m} and L_{z} (ξ) is the dimensionless integral
The integral L_{z}(0) can be computed analytically:
where the spheroid’s eccentricityThe potential inside the spheroid is given by
It is clear from this equation that the field inside the spheroid is uniform and parallel to the external field. In addition, the polarization P (dipole moment per unit volume) inside the spheroid is uniform and is given by
where ε_{0} is the absolute permittivity of free space.
The electric field outside the spheroid isfor points on the z axis,
The quantity in square brackets in this equation is the onaxis electric field enhancement factor, which we henceforth denote as F. It has the value F_{tip} = (1 + s)/[1 + sL_{z}(0)] at the spheroid’s tip (x, y, z) = (0, 0, c), and approaches F = 1 as z → ∞. Large F_{tip} values can be achieved when the quantity 1 + sL_{z} (0), which is controlled by the spheroid’s aspect ratio c/a and dielectric constant, is small in magnitude.
We now turn to our DDAbased simulations. The technical aspects of these simulations have been extensively reviewed [5,13]; we therefore report only those computational details that are specific to the simulations presented here. We model a spheroid as a collection of N contiguous cubic subunits, with edges of length d, centered at the positions (xyz) = (n_{x}dn_{y}dn_{z}d); here, (n_{x}n_{y}n_{z}) is an integer triple that satisfies where n_{max} is an integer that determines the discretization level of the spheroid. The edge length d is chosen so that the volume enclosed by the collection of cubic subunits is equal to the spheroid volume. The linear algebraic equations that determine the dipole moments m_{j} of the individual subunits (here j is an index that distinguishes individual subunits) are solved using the complexarithmetic implementation of the GMRES algorithm described by Frayssé et al. [14]; we terminate the algorithm and record the dipole moments m_{j} once the normwise backward error drops below 10^{−6}. We obtain the wavelengthdependent dielectric function of Ag via linear interpolation of the data points compiled by Lynch and Hunter [15]; as our main goal in the present work is not to provide results for comparison with experiment, but to compare the accuracy of the results obtained in simulations with and without the local environmental correction to the polarizabilities, we neglect finitesize corrections to the dielectric constant that arise from electronic scattering from the spheroid surface [16]. Henceforth, we set ε_{m} = 1 (corresponding to vacuum as the medium surrounding the spheroid) and E_{0} = 1 au; all of the results we report are scaled by 1/E_{0}, so the numerical value of E_{0} is ultimately irrelevant.
First, we examine the polarization P induced in a metallic nanoparticle by a uniform static external electric field. We consider a prolate spheroid with a = 10 nm, c = 40 nm, and dielectric constant ε = 12.26 + 0.84i (corresponding to an excitation wavelength of λ ≈ 570 nm). For this aspect ratio and dielectric constant, the quantity 1 + sL_{z} (0) is purely imaginary and small in magnitude: 1 + sL_{z} (0) ≈ 0.0633i. We use the DDA to simulate this spheroid at several levels of discretization, ranging from N = 6041 subunits (d = 1.405 nm) to N = 24679 subunits (d = 0.879 nm). We divide the dipole moment m_{j} of each subunit by the subunit volume d^{3} to obtain the polarization P_{j} for each subunit; we then divide the magnitude of this vector by the magnitude of the exact polarization vector defined in Eq. 7 to obtain a dimensionless relative polarizationfor each subunit. This quantity has the valuewhen the magnitude of a subunit’s dipole moment m_{j} is consistent with the exact uniform polarization given by Eq. 7.
For DDAbased simulations employing CM polarizabilities, Fig. 1 shows how the mean and standard deviation ofevaluated over the N subunits in a given spheroid, depend on the subunit edge length d. We see that for all of the spheroids considered here, the meanvalue differs considerably from the value. For each spheroid, the standard deviation of thevalues is about 0.2, indicating that the polarization within the spheroid is rather nonuniform—in contrast to the exact result given by Eq. 7—and does not become more uniform as the subunits become smaller; for two of the spheroids, Fig. 2 depicts graphically the large subunittosubunit variations inthat are observed using CM polarizabilities. By comparison, when we use the subunit polarizabilities of [8] that include the LEC, our DDAbased simulations produce subunit dipole moments that give for each subunit in the spheroid, indicating that the magnitudes of the dipole moments are in exact agreement with Eq. 7. This is not much of a surprise, because the corrected polarizabilities given by Rahmani et al. [8] are defined so that, when used in DDAbased simulations, they reproduce exactly the positiondependent polarization inside an object immersed in a static external electric field [8,9]. What Figs. 1 and 2 show is that DDAbased simulations that use CM polarizabilities may fail in this regard, and that this failure is not simply a result of the discretization that necessarily accompanies the DDA.
Figure 1. Mean (boxes) and SD (circles) of the relative subunit polarizationsas a function of subunit edge length d, for DDAbased simulations of a prolate Ag spheroid with semiaxes a = 10 nm and c = 40 nm; the simulations employ uncorrected CM polarizabilities
Figure 2. Relative polarizationsderived from DDAbased simulations of a prolate Ag spheroid with semiaxes a = 10 nm and c = 40 nm; the simulations employ uncorrected CM polarizabilities. Each subunit is represented by a square colored according to the relative polarization scale shown at the bottom of the figure. The upper panel gives the results for a spheroid modeled using N = 6041 subunits of edge length d = 1.405 nm; the lower panel gives the results for N = 24679 subunits of edge length d = 0.879 nm. Only subunits with x = 0, y ≥ 0, and z ≥ 0 are shown
We now examine the nearfield properties of the Ag nanoparticle, focusing on the localized enhancement of the applied electric field near the nanoparticle’s surface. Figure 3a compares the magnitude F of the exact onaxis electric field enhancement factor F [which is a complex quantity because ε is complex; see Eq. 8] near the spheroid’s sharp tip with the DDAbased results obtained for N = 24679 subunits using both CM and LEC polarizabilities. The enhancement factor computed using LEC polarizabilities is in good agreement with the exact result, even at points within 0.4 nm (which is less than onehalf of the subunit separation d) of the spheroid’s surface (the electric field varies discontinuously across the spheroid’s surface, and no DDAbased simulation will be able to model this discontinuous change; it is therefore unreasonable to expect these simulations to give accurate F values just outside the spheroid’s surface). By contrast, the simulation that employs CM polarizabilities substantially underestimates F. In Fig. 3b, we show how the values of F computed at z = 41 nm (1 nm away from the sharp tip) vary with d over the range of discretizations considered here; although the F values computed using CM polarizabilities vary slightly as d decreases, it appears that very small subunits will be needed before the CM result approaches the exact one. On the other hand, the F values computed using LEC polarizabilities are within a few percent of the exact result at all levels of discretization.
Figure 3. Magnitude of F, the onaxis electric field enhancement factor, for a prolate Ag spheroid with semiaxes a = 10 nm and c = 40 nm. a Dependence of F on position z; the point z = 40 nm is at the spheroid’s sharp tip. Solid line gives the exact result of Eq. 8; boxes and circles give the results of DDAbased simulations with N = 24679 subunits employing LEC and CM polarizabilities, respectively. b Dependence of F at z = 41 nm on the edge length d of the DDA subunits. The dotted line at F = 102.7 gives the exact result of Eq. 8; boxes and circles give the results of DDAbased simulations employing LEC and CM polarizabilities, respectively
To gain more insight into the relative performance of DDAbased simulations employing CM and LEC polarizabilities, we use the simulations to compute the electric field enhancement factor in the vicinity of the spheroid’s sharp tip, and compare these enhancement factors to reference results obtained by numerically differentiating the exact electrical potential Φ_{out} defined in Eq. 3. To partially mitigate the discretization effects that are inherent in DDAbased simulations, we rotationally average the field enhancement factor obtained from these simulations by computing it on ten evenly spaced dihedral planes containing the spacefixed z axis and then averaging the enhancement factors obtained for each dihedral plane.
Fig. 4 shows, for N = 24679 subunits, how the magnitudes of the enhancement factors computed using DDAbased simulations compare with the reference results derived from Eq. 3. Although neither of the DDAbased simulations can predict accurately the field enhancement factors at the surface of the spheroid (because the electric field varies discontinuously across the spheroid’s surface, as previously noted), the shape, size, and internal structure of the spheroid’s nearfield hot spot are modeled fairly well by the simulations that employ LEC polarizabilities. The DDAbased simulations that employ CM polarizabilities, by contrast, yield a hot spot that is too small and whose peak intensity is too low.
Figure 4. Magnitude of the electric field enhancement factor near the sharp tip of a prolate Ag spheroid with semiaxes a = 10 nm and c = 40 nm. Panels marked “REF” are reference results obtained by differentiating Eq. 3; “CM” and “LEC” indicate the results of DDAbased simulations employing CM and LEC polarizabilities, respectively. The DDAbased simulations use N = 24679 subunits of edge length d = 0.879 nm to model the spheroid; the colored squares in the legend at the bottom of the figure are scaled to have the same edge length
In summary, we have modeled the response of a nanoscale Ag prolate spheroid to an external electric field using DDAbased simulations that employ subunit polarizabilities that either include or omit a local environmental correction. We invoke the electrostatic approximation, in which the incident field is assumed to be spatially uniform and static, but the spheroid’s dielectric constants is taken from the wavelengthdependent dielectric function of bulk Ag; this allows us to compare the predictions of the DDAbased simulations to exact results obtained by solving Laplace’s equation for prolate spheroids in a uniform static external field. We have chosen a dielectric constant for the spheroid that maximizes the electric field enhancement factor at the spheroid’s sharp tip. The predictions of DDAbased simulations that employ LEC polarizabilities are much closer to the exact results than are those of DDAbased simulations that employ CM polarizabilities; simulations using CM polarizabilities yield a nearfield hot spot that is too small and field enhancement factors that are too low. We therefore conclude that DDAbased simulations of metallic nanoparticles that employ uncorrected CM polarizabilities may give inaccurate predictions of the particle’s spatially resolved nearfield properties, even at locations some distance away from the particle’s surface.
Acknowledgments
This work was supported by grants from the University of Tennessee Honors Program (A.E.D.) and the US Department of Energy (R.J.H.). R.J.H. thanks L. Blocker, L. Dixon, and E. Read (University of Tennessee Libraries) for bibliographic assistance.
Open Access
This article is distributed under the terms of the Creative Commons Attribution Noncommercial License which permits any noncommercial use, distribution, and reproduction in any medium, provided the original author(s) and source are credited.
References

Ann. Phys. Leipz.. 1908, 25:377.
COI number [1:CAS:528:DyaD1cXhtFWqtQ%3D%3D]; Bibcode number [1908AnP...330..377M]
Publisher Full Text 
Phil. Trans. R. Soc. Lond.. 1857, 147:145. Publisher Full Text

Astrophys. J.. 1973, 186:705.
Bibcode number [1973ApJ...186..705P]
Publisher Full Text 
Astrophys. J.. 1988, 333:848.
COI number [1:CAS:528:DyaL1cXmtFyrtbs%3D]; Bibcode number [1988ApJ...333..848D]
Publisher Full Text 
Jackson JD: Classical electrodynamics. Wiley, New York; 1999.

Astrophys. J.. 1993, 405:685.
Bibcode number [1993ApJ...405..685D]
Publisher Full Text 
Rahmani A, Chaumet PC, Bryant GW:
Opt. Lett.. 2002, 27:2118.
Bibcode number [2002OptL...27.2118R]
PubMed Abstract  Publisher Full Text 
Rahmani A, Chaumet PC, Bryant GW:
Astrophys. J.. 2004, 607:873.
Bibcode number [2004ApJ...607..873R]
Publisher Full Text 
J. Opt. Soc. Am. A. 2004, 10:2023.
Bibcode number [2004JOSAA..21.2023C]
Publisher Full Text 
Stratton JA: Electromagnetic theory. McGraw Hill, New York; 1941.

Bohren CF, Huffman DR: Absorption and scattering of light by small particles. Wiley, New York; 1983.

J. Opt. Soc. Am. A. 1994, 11:1491.
Bibcode number [1994OSAJ...11.1491D]
Publisher Full Text 
Frayssé V, Girard L, Gratton S, Langou J:
ACM Trans. Math. Soft.. 2005, 21:228. Publisher Full Text

Lynch DW, Hunter WR: Handbook of optical constants of solids. Edited by Pali ED. Academic, San Diego; 1985.
Part II, Subpart 1

J. Chem. Phys.. 1983, 79:6120.
Bibcode number [1983JChPh..79.6130K]
Publisher Full Text