Simulation and Optimization of Nonperiodic Plasmonic NanoParticles
 Author: Akhlaghi Majid, Emami Farzin, Sadeghi Mokhtar Sha, Yazdanypoor Mohammad
 Organization: Akhlaghi Majid; Emami Farzin; Sadeghi Mokhtar Sha; Yazdanypoor Mohammad
 Publish: Journal of the Optical Society of Korea Volume 18, Issue1, p82~88, 25 Feb 2014

ABSTRACT
A binarycoupled dipole approximation (BCDA) is described for designing metal nanoparticles with nonperiodic structures in one, two, and three dimensions. This method can be used to simulate the variation of near and farfield properties through the interactions of metal nanoparticles. An advantage of this method is in its combination with the binary particle swarm optimization (BPSO) algorithm to find the best array of nanoparticles from all possible arrays. The BPSO algorithm has been used to design an array of plasmonic nanospheres to achieve maximum absorption, scattering, and extinction coefficient spectra. In BPSO, a swarm consists of a matrix with binary entries controlling the presence (‘1’) or the absence (‘0’) of nanospheres in the array. This approach is useful in optical applications such as solar cells, biosensors, and plasmonic nanoantennae, and optical cloaking.

KEYWORD
Plasmonic nano particles , Optimization algorithm , Binary PSO

I. INTRODUCTION
Metal nanoparticles have broad applications in electronics, photonics, chemical sensing, and imaging [13]. Recently there has been a growing interest in plasmonic nanoparticles. They have strong effects on light in the visible and invisible regions of the photonic spectrum for applications such as Raman scattering [4], radiative rate enhancement [5], solar cells [6], and optical biosensors [7]. Since the optical properties of these particles depend on their sizes and shapes, one of the desired goals is to control the shapes of metal and semiconductor nanoparticles. Shape control has been successfully demonstrated for gold nanoparticles using nonionic surfactants, silver underpotential deposition, and nanoreplica molding [811]. Nanoreplica molding has been demonstrated as a lowcost method for manufacturing a variety of devices comprised of nano structured surfaces. Recently a plasmonic nanodome array was fabricated by the nanoreplica molding process [10]. In addition to the recent interest in shape control of nanoparticles, the optical properties of noble metal particles, with their intense colors, have fascinated scientists since the turn of this century. For example, spherical gold nanoparticles show a strong absorption band in the visible region of the electromagnetic spectrum at about 520 nm [12]. In addition to the shape, size, and material of nanoparticles, the properties of transmitted light strongly depend on the localized positions and gaps between nanoparticles. Plasmonic nanoparticles with periodic structure have been reported [13, 14]. One of the most promising plasmonics nanoparticle platforms is studying the effect of deterministic aperiodic structure of nanoparticles on the properties of transmitted light [15]. Such structure, which is intermediate between a disordered system and a periodic one, enables unique control and manipulation of spatially localized plasmonic states over broadband frequency and angular spectra [16].
In structures with periodic plasmonic nanoparticles the characteristics of transmitted light can be calculated efficiently using the coupleddipole approximation (CDA), originally developed by Purcell and PennyPacker [17] and improved by Draine [18]. In this paper, an alternative design method called the binary coupled dipole approximation (BCDA) is described. Utilizing this method, the effect of a nonperiodic structure of nanoparticles on light properties could be studied. A big advantage of this method is its application to optimization problems by utilizing binary evolutionary computational techniques, such as binary genetic algorithms and binary particle swarm optimization. Moreover, higher absorption, scattering, and extinction coefficients are obtained using BPSO and BCDA methods. These improvements can be helpful in optical applications such as solar cells, optical bio sensors, and plasmonic nanoantennae.
II. METHODS
2.1. Binary Couple Dipole Approximation
A schematic diagram of a threedimensional array of arbitrarily shaped objects periodically arranged in the x y z plane is shown in Fig. 1(a). In the first time step the x, y, and z positions of each particle are calculated. To assign the position of each particle, it is assumed that there is a cube filled with nanoparticles equally distant from each other. The Matlab code for determining the positions of N particles is presented in Fig. 1(b). In this code, r is a particle’s position vector, N1, N2, and N3 represent the odd number of particles in the i, j, and k directions, and d represents the distance between the particles. B(i, j, k) is a threedimensional array of binary numbers. The values 1 and 0 in the B matrix refer respectively to the presence and absence of particles in the desired locations. Next, the incident electric field is simulated. The scattered light may be illuminated by any form of incident electric field, the simplest of which is a plane wave. To calculate the Efield of each dipole at rj, the time harmonic component iωt of the Efield is left out. The Local field arising from incident light with polar (θ) and azimuthal (φ) angles at each particle is:
Where:
For an incident field with ppolarization, the following can be written:
and for an incident field with spolarization:
When the applied field is parallel to one of the principle axes, the polarizability α is [19]:
where V is particle volume, ε_{m} is the dielectric function of the surrounding medium, L is the shape factor, and . One advantage of the BCDA method is studying the effect of an ellipsoidal particle on light. For oblate and prolate spheroids (b=c), the following analytical expression can be given for L1 as a function of eccentricity e [19]:
where a, b, and c are seminal excess of an ellipsoid (a>b> c) and, for a sphere, L_{1}=1/3.
The dipole moment induced in a single particle by a local electric field is given by:
Here, is the induced dipole moment, α_{i} is the polarizability of the particle centered at ,
E _{loc} is the local electric field, and ε_{0} is the permittivity of free space. The local field arises from two sources: the first term is incident light, , and the second term is superposition of the retarded fields from each of the other N1 radiating dipoles in the array. Combining these terms leads to the local field as follows:Here, , k is wave vector, and Aij are 3x3 matrices representing the interaction of two particles i and j. Now, 3N unknown dipole moments Pj should be solved for in the following exactly determined system of 3N linear equations:
where A is a 3N×3N matrix containing N_N of A_{ij} 3x3 tensors, where N is the number of dipoles and both p and E_{inc} are 3N vectors (i.e. each of N particles is represented by a 3vector). These matrices are shown in Fig. 3. When this set of 3N complex linear equations is solved (
A .P =E _{inc}), the P array of selfconsistent dipole moments is obtained. The optical properties may be then calculated from this dipole array.2.2. FarField Properties of Light
As explained in the previous section, optical absorption, scattering, and extinction can be directly calculated from the dipole array. The extinction and scattering coefficients are [12]:
and the absorption coefficient is Q_{abs}=Q_{ext}Q_{sca}.
2.3. Binary Particle Swarm Optimization Algorithm
The PSO algorithm is an optimization procedure inspired by a colony such as birds which can improve its behaviors [20]. Any element of this colony is called a particle and moves in an ndimensional space, correcting its trajectory based on the previous actions of itself and its neighboring particles. For each particle, velocity and displacement are updated based on the following relations [21]:
where
k is the number of the current iteration,n is the number of particles, w is the inertia weight,c _{1 }andc _{2} are acceleration parameters, and finallyr _{1} andr _{2} are random parameters between 0 and 1. The best position for thei th particle which has been stored so far is represented as [22]:All the (
p_best ) are evaluated by a fitness function. The best particle among allp_best is represented as (Gbest_best_value ), which in a minimization problem is the smallest member of the vector while for a maximization problem it is the largest member of . PSO was designed for continuous problems, but cannot deal with discrete problems. A new version of PSO, called Binary PSO (BPSO), was introduced by Kennedy and Eberhart in 1997 and applied to discrete binary variables. After that, many optimization problems in various areas were solved by this method. The position in BPSO is represented by a binary vector and the velocity is still a floatingpoint vector; however, velocity is used to determine whether the probability changes from 0 to 1 or from 1 to 0 when the positions of particles are being updated. The equation for updating the positions is then replaced with:In order to achieve maximum absorption, extinction, and scattering coefficient spectra, BPSO has been used to design a nonperiodic array of plasmonic nanospheres. The presence or absence of nanospheres could be controlled by a matrix with binary entries. As seen, BPSO is an algorithm that minimizes a cost function. To use the BPSO algorithm for maximizing the extinction coefficient spectrum, the algorithm should minimize the following function:
and for maximizing the absorption and scattering coefficient spectra:
where i is the wavelength number and Qext, Qabs, and Qsca are the extinction, absorption, and scattering coefficients respectively.
III. NUMERICAL RESULTS
The BCDA method is applied and the absorption spectrum is plotted for periodic and nonperiodic arrays of gold nanospheres with 1 nm diameter and 2 nm edgeto edge separations in 1, 2, and 3 dimensions. Figure 4(b) shows the absorption spectrum as a function of wavelength in the range of 300600 nm for a 1D nonperiodic array. The present simulations show that the absorption spectrum strongly depends on localized positions of nanospheres. On the other hand, further enhancement in absorption coefficient is possible for nonperiodic structure compare to the periodic array (Fig. 4(a)). These results are quite consistent with the results obtained using Mie theory [12]. Figure 5(a) shows this possibility for a 2D nonperiodic array. One benefit of increased material absorption coefficients is in solar cell. Then the effect of nonperiodic arrays of gold nanospheres on the absorption spectrum with a 3D nonperiodic array is investigated. Figure 5(b) shows the absorption spectrum as a function of wavelength for nonperiodic array in 3D coordinates. Figure 6(b) shows the extinction spectrum of a 2D nonperiodic array with 3x5 particles compared to that of periodic array (Figure 6(a)). The large extinction coefficient observed here was due to the interaction of metal nanoparticles via their strong nearfield coupling. In this case the plasmon resonance of a nanoparticle acts as a light resonator and couples light to other resonators in the vicinity [23]. The BCDA method is applied and extinction coefficient spectra are plotted for 1, 2, 3, and 4 layers of gold nanospheres (5×5 array) with 1 nm diameter and 2 nm edgetoedge separations. Figure 7 shows these extinction spectra in the wavelength range of 300600 nm. Our simulations show that the extinction coefficient spectrum strongly depends on the number of layers of nanospheres, and further enhancement of the extinction coefficient is possible with multi periodic layers.
IV. OPTIMIZATION RESULTS
Binary PSO and BCDA are used to optimize the farfield properties of transmitted light for optical applications. Our goal is to maximize Qext by optimizing 15 binary particles in 2D arrays (3×5). Figure 8(a) shows the extinction spectrum of the optimized array (3×5) compared to that of the periodic array (Fig. 6(a)). The maximum extinction in the entire range of 300600 nm is 9.8, which occurs at 560 nm. The large extinction coefficient observed here is due to the interaction of metal nanoparticles via their strong nearfield coupling. The convergence characteristic of the BPSO method is shown in Fig. 8(b). As seen, this algorithm had good performance in terms of convergence speed.
Next, the BPSO algorithm is applied to maximize absorption and scattering coefficients by optimizing the selection of binary particles. Figures 9(a) and Figures 9(b) show the absorption and scattering spectra of the optimized array compared to those of the periodic array. One possible application of this method can be reduction of the scattering factor to have an optical invisibility cloak based on scattering cancellation. In addition to sensor, solar cell, and cloaking applications, BPSO can enhance the Raman scattering factor.
V. CONCLUSION
A numerical solution called the binarycoupled dipole approximation (BCDA) was described for simulating the effects of nonperiodic plasmonic nanoparticles on light properties. This method is useful for optimization problems to find the best nanoparticle array from all possible arrays. Finally, it was shown that using a nonperiodic structure could lead to higher extinction, absorption, and scattering coefficients. This approach can be useful in optical applications such as solar cells, bio sensor, and plasmonic nanoantennae.

[FIG. 1.] (a) 3dimensional array of nanoparticles. (b) Matlab code for determining the position of N number of nanoparticles.

[]

[]

[]

[]

[]

[FIG. 2.] Prolate (left, a>b=c) and oblate (right, a<b=c) ellipsoids.

[FIG. 3.] A, P, Einc matrix from equation 10.

[]

[]

[]

[]

[]

[]

[]

[]

[]

[]

[]

[]

[]

[]

[]

[FIG. 4.] (a) QAbs for a 1D periodic array with B=[1 1 1 1 1] . (b) QAbs for a 1D non periodic array with B=[1 0 1 0 0].

[FIG. 5.] (a) Absorption coefficient for a 2D nonperiodic structure with B=[0111 0;01001;10110;01100;01110], (b) Absorption coefficient for 3D nonperiodic structure with B = [10001; 00010;00100;01011;10100], B(:,:,2) = [10101;0101 0;0 0 0 0 0; 1 0 0 1 0; 0 0 1 0 0].

[FIG. 6] (a) Extinction spectrum of a periodic 2D array with 3*5 particles. (b) Extinction spectrum of nonperiodic 2D array with 3*5 particles with B = [0 1 0;0 1 0;1 0 0;1 0 0;1 0 1].

[FIG. 7.] Extinction coefficient spectra for (a) a 1layer periodic array, (b) a 2layer periodic array. (c) a 3layer periodic array, and (d) a 4layer periodic array.

[FIG. 8.] (a) Extinction spectrum of the optimized array (3*5) with B = [ 0 1 0;0 1 1;0 0 1;1 0 0;1 0 0], (b) Convergence characteristic of BPSO.

[FIG. 9.] (a) Absorption spectrum of the optimized array (3*5) with B=[0 0 1;101;1 10;010;000]. (b) Scattering spectrum of the optimized array (3*5) with B = [0 01; 0 0 1;100;100;000].