## Abstract

We explore the passive optical sorting of plasmon nanoparticles and investigate the optimal wavelength and optimal beam shape of incident field. The condition for optimal wavelength is found by maximising the nanoparticle separation whilst minimising the temperature increase in the system. We then use the force optical eigenmode (FOEi) method to find the beam shape of incident electromagnetic field, maximising the force difference between plasmon nanoparticles. The maximum force difference is found with respect to the whole sorting region. The combination of wavelength and beam shape study is demonstrated for a specific case of gold nanoparticles of radius 40nm and 50nm respectively. The optimum wavelength for this particular situation is found to be above 700nm. The optimum beam shape depends upon the size of sorting region and ranges from plane-wave illumination for infinite sorting region to a field maximising gradient force difference in a single point.

© 2011 OSA

## 1. Introduction

The preparation of highly monodispersed colloidal solutions of plasmon nanoparticles is crucial for any application where a narrow size distribution is required particularly to exploit their plasmonic properties. This includes biomolecular sensors [1], localised heaters [2] and photo-thermal imaging [3]. As most of these applications require colloids prepared in sterile way, the application of optical sorting offers an ideal solution not only because of the non-contact nature of sorting but also because of the exceptional sensitivity on particle size, shape and refractive index. Two main optical sorting approaches exist – namely *active* and *passive*. Active sorting techniques use fluorescent signals [4], optical switches [5, 6] and real-time computation [7, 8] for particle recognition that then subsequently acts as a trigger for another part of the system where optical (or other) forces are used to separate particles of different properties into separate streams. However, the need to use a trigger may add complexity and be impractical in many instances. For this reason, the last decade has seen the emergence of passive sorting methods which rely entirely on the different physical response of various particles to an extended optical field, commonly referred to as an optical potential energy landscape. Such passive sorting offers exceptional size and refractive index sensitivity [9, 10] and has been demonstrated in a number of geometries both with and without the presence of microfluidic flow. Sorting of particles in closed chambers with static fluid has been realised both by means of moving interference pattern [11, 12] or Bessel modes [13]. However, in the majority of sorting applications, laminar fluid flow perpendicular to the optical forces is employed and the separation of particles is either realised solely by scattering force differences [14] in aperiodic optical patterns or by means of structured light fields [10, 15–19] creating periodical optical potential energy landscapes. In all these methods, the determination of the precise form of the applied optimal field, such that the maximum sorting sensitivity on size, shape and refractive index is achieved, remains an open question. To date it is also to be noted that passive optical sorting has been mainly considered solely with respect to micron-sized dielectric particles and cellular media. It is intriguing to consider how passive optical sorting may be extended to the domain of plasmon nanoparticles.

Here, we study passive optical sorting of plasmon nanoparticles and present a general two step approach that can be used to design the optimal illumination for sorting plasmon nanoparticles. In the first step, we find the optimal wavelength of illumination such that the separation is achieved with minimum temperature increase in the system. This is an important consideration in plasmonic systems as the excessive heat increases diffusion and convective effects, which is counter-productive in any sorting application. In the second step, the optimisation of illumination shape is realised using our force optical eigenmode (FOEi) [20] method, which can be readily applied to determine the optimal laser illumination field for sorting. The approach extends previous studies of forces on plasmon nanoparticles [21, 22] by considering the shape of the field, its wavelength and heating in the system at the same time. This paper is divided into three sections. In the first section, we present the FOEi method and derive formulae leading to optimised beam shape for sorting. The second part focuses on exploiting the plasmon resonances in the system and finding the optimal wavelength for sorting. The third part uses the optimal wavelength as input and extends the FOEi method to optimise the force difference over the whole sorting region.

## 2. Description of FOEi method

Our aim is to optimise the incident electromagnetic field **E**
* _{inc}* in a way that maximises the exerted force

*F*

**=**

^{u}**F**·

**u**in a specified direction

**u**(Fig. 1). Since our incident field of angular frequency

*ω*can be decomposed into a sum of

*μ*monochromatic plane waves (

*e*), we can write (using summation over repeating indices)

^{iωt}*a*are the complex expansion coefficients and ${\mathbf{\text{E}}}_{inc}^{\mu}$ are the incident plane waves. The

^{μ}*a*coefficients modulate both phase and amplitude of each and single plane wave independently. The scattered field generated upon interaction with the particle has the same expansion coefficients, so that we can write ${\mathbf{\text{E}}}_{sca}\hspace{0.17em}=\hspace{0.17em}{a}^{\mu}\hspace{0.17em}{\mathbf{\text{E}}}_{sca}^{\mu}$. The final field

^{μ}**E**, which is a sum of incident and scattered field, can then be written as

**E**

*is a solution of the scattering problem for an incident field given by plane wave ${\mathbf{\text{E}}}_{inc}^{\mu}$.*

^{μ}The optical-cycle averaged electromagnetic force in the direction **u** is given by

*n*is outward unit normal to an element d

_{j}*s*of the curve

*C*enclosing the particle for which we optimise the force and 〈·〉 is optical-cycle average. The Maxwell stress tensor,

*σ*, can be written for the final field

_{ij}**E**as [23]

*a*we obtain

^{μ}*M*are given by the line integrals in Eq. (7) and

^{μν}**a**is the vector form of

*a*. We remark that the matrix is Hermitian (

^{μ}**M**=

**M**

^{†}) and thus its eigenvalues are real. This means that the force

*F*

**is in a symmetric sesquilinear form, which is just an extension of quadratic form to complex numbers. Any symmetric sesquilinear form can be visualised as an ellipsoid with the length of principal axes equal to the eigenvalues**

^{u}*λ*of the matrix

_{n}**M**. This has far reaching implications for our optimisation process since the surface of the ellipsoid generated by the symmetric sesquilinear form is extremized at the end points of principal axes. This means that finding the eigenvalues

*λ*of the matrix

_{n}**M**and selecting the largest one from the set extremizes our problem. The eigenvector (force optical eigenmode) ${\mathbf{\text{a}}}_{n}^{max}$ corresponding to maximum eigenvalue ${\lambda}_{n}^{max}$ (given by ${\mathbf{\text{Ma}}}_{n}^{max}\hspace{0.17em}=\hspace{0.17em}{\lambda}_{n}^{max}\hspace{0.17em}{\mathbf{\text{a}}}_{n}^{max}$), then provides the necessary information about amplitude and phase of incident plane waves ${\mathbf{\text{E}}}_{inc}^{\mu}$ in Eq. (1) so that the force is maximised. Experimentally, the optimised ${\mathbf{\text{a}}}_{n}^{max}$ can be created using spatial light modulator (offering control of both phase and amplitude) in the system placed in conjugate plane [20] with respect to the back-focal plane of a microscope objective.

We remark that the above described method optimises the force on one type of particle at a single point only. However, the method can be easily extended to provide the optimised illumination for the force difference over the whole sorting region for two types of nanoparticles. This is discussed later in the paper.

*Numerical considerations:* We use COMSOL Multiphysics v4.1 RF module in scattering formulation to calculate the total field solutions **E**
* ^{μ}* for corresponding incident plane waves
${\mathbf{\text{E}}}_{inc}^{\mu}$. We first find the solutions

**E**

*for particle*

^{μ}*p*

_{1}. We subsequently use the solutions

**E**

*to find the elements ${M}_{1}^{\mu \nu}$ of matrix*

^{μ}**M**

_{1}and determine the eigenvalues and corresponding eigenvectors. The principal eigenvector gives the optimum force for particle

*p*

_{1}in a single position. The same procedure is repeated for a second type of particle,

*p*

_{2}, delivering matrix

**M**

_{2}. The matrices

**M**

_{1}and

**M**

_{2}then encode all the information about interactions of the incident fields with the particles. Finding the matrix elements

*M*of matrix

^{μν}**M**is computationally very intensive as combinations of

*N*(

*N*+ 1)/2 solutions need to be constructed and integrated over a sphere boundary. Here

*N*denotes number of plane waves in the angular spectrum representation and thus the number of pixels on spatial light modulator. As the azimuthal discretization of angular spectrum representation increases the number of combinations in 3D significantly, we have restricted the simulations to 2D to illustrate the method.

We can obtain educated estimates of 3D values from 2D values by extruding the 2D circle of radius *r* by *d* such that it creates a cylinder with a volume equal to the volume of the sphere with the same radius. The extrusion factor *d* is given by

## 3. Plasmon resonances in the system

We choose as our testing system gold nanoparticles [24] of radius *r*
_{1} = 50nm and *r*
_{2} = 40nm. We consider a substrate of glass with refractive index *n _{g}* = 1.5. The particles are assumed to be dispersed in water with

*n*= 1.33 (Fig. 1).

_{w}Plasmonic resonances offer exceptional sensitivity on size. However, for nanoparticles of very similar sizes, the force difference generated solely by the plasmon resonance in the system is still rather small. In case of dielectrics, increasing the intensity offers simple solution in such a situation, however plasmonic resonances are associated with non-negligible heat generation and as such increasing the intensity may produce increased diffusion rates and convective currents, which will interfere with sorting efforts. It is thus beneficial to find a wavelength for which the sorting effects are maximised and heating is minimised. Further, multiple laser wavelengths in a counter-propagating geometry with carefully adjusted powers can be employed to exploit the small plasmon resonance differences of nanoparticles. Our method can solve for optimum illumination for each and single wavelength separately, however, the use of multiple laser wavelengths is expensive and increases complexity of the system. It is therefore of interest to optimize the force differences using a single laser wavelength.

We choose the p-polarisation for the incident plane waves as the kind of plasmon resonance supported by the sphere appears in 2D for p-polarisation. The s-polarisation would only induce movement of electrons along the infinite cylinder. Figure 2(a) shows the scattering *Q _{sca}* and absorption

*Q*efficiencies for

_{abs}*r*

_{1}= 50nm gold nanoparticle. Note that the 3D efficiencies calculated from Mie theory and the corresponding 2D efficiencies (transformed to 3D) for p-polarisation follow very similar pattern, which differs only in amplitude and a slight blue shift of 2D resonance peaks with respect to 3D resonances. Figure 2(b) shows the 3D forces along the substrate acting on the gold nanoparticles calculated for plane-wave incident at near critical angle of

*θ*= 64°. The slight shift in resonances due to the different sizes of nanoparticles creates a force difference improvement with a peak around 550nm (black curve in Fig. 2(b)).

If we do not take into account the proximity of the spheres to the surface (Faxen’s correction [25]) and assume we are in a low Reynolds number regime, then the drag force is given by the Stokes equation. Neglecting inertial effects we can equal optical and drag force and obtain the expression for the particle speed in 3D as

where the force*F*

_{2D}was calculated in 2D using Eq. (8) and the dynamic viscosity of water at

*T*= 20°

*C*is

*η*= 1.002 × 10

^{−3}Pa · s. Figure 2(c) shows the speed difference generated by force differences in Fig. 2(b) and the average temperature increase as the particle enters the sorting field. The estimate of average temperature increase was calculated using [26] (neglecting particle movement and proximity of glass surface)

*Q*

_{1}and

*Q*

_{2}are 3D powers of heat generation in nanoparticles and

*κ*

_{0}= 0.6W · m

^{−1}· K

^{−1}is the thermal conductivity of water. This temperature increase is reached very quickly as the particle enters the sorting field. Results in Fig. 2(c) and Fig. 2(d) suggest, that heating at wavelengths above 700nm is quite low and the sorting speed remains high. For this reason we set the vacuum wavelength for our method to

*λ*

_{0}= 700nm. Please note that a slight (30nm) blue shift in the resonances of 2D case does not have significant impact on the choice of this wavelength for the 3D scenario.

## 4. Optimising the force difference in region of interest (ROI)

Using the optimum wavelength from previous section we can proceed with optimisation of the beam shape. We introduce a 10nm separation between the lowest point of particle and interface, which closely mimics a typical experimental situation. The set of *N* incident plane waves defined by **k**
* _{θ}* vectors (

*θ*= 〈−70°,...,70°〉 with a step of 2°) is used for discretization (see Fig. 1). The limits of

*θ*correspond to experimental limitation for

*NA*= 1.4 oil immersion objective. The goal is to optimise the force difference along

**u**= (1, 0).

Using matrices **M**
_{1} and **M**
_{2}, the equation for force difference in a single point for our choice of particles is

*x*direction and the second for −

*x*direction. In subsequent discussions we always choose the solution optimising the force in +

*x*direction. The field locally optimising the force difference for gold nanoparticles of our choice is on Fig. 3. Notice that the field creates a very strong field gradient in the +

*x*direction around point

*x*= 0, where we want to maximise the force difference for our testing particles. Also notice that the back focal plane pattern corresponding to this field has significant contributions from plane waves propagating in the opposite −

*x*direction. Although this might seem surprising, we need to realise that the final goal of our method is to interfere the plane waves in such a way to create the strongest gradient in +

*x*direction. Apparently the counter-propagating waves increase the number of degrees of freedom for efficient interference leading to strong gradient and are thus utilised automatically by the FOEi method. It also make sense that the increased intensity at the back focal plane appears for near critical angle plane waves as those plane waves contribute the most to the intensity near the glass/water interface. Note that the phase at the back focal plane is also significantly altered to maximise the force difference.

So far our method optimises the force difference locally. To expand this approach to a larger region we need to optimise Δ*F* over a certain range, in our case line segment defined by *x* = 〈*−l,l*〉. Displacing the particle in *x*-direction causes the particle to experience different relative phases between the fields **E**
* ^{μ}*. Since we use combination of solutions to calculate matrix

**M**this relative phase can be taken into account using

*k*= (2

_{w}*π*/

*λ*

_{0})

*n*and ${\theta}_{t}^{\mu}$ is the angle of the plane wave after the glass/water interface for each incident plane wave ${\mathbf{\text{E}}}_{inc}^{\mu}$.

_{w}*M*is the matrix calculated for particle at position

^{μν}*x*= 0. Using the translation relations for matrix

*M*we obtain

^{μν}*l*. Finding the optimum illumination for a ROI of any size is thus very efficient.

Optimal illumination for ROI sizes of *l* = 500nm (Fig. 4) and *l* = 5*μ*m (Fig. 5) differs significantly from the single point optimised problem (Fig. 3). The optimised field corresponds in its bulk to the focusing of light into ROI. The solution is quite close to the Gaussian beam send to the edge of the back focal plane of the objective. However, the phase for plane waves above critical angle is significantly modulated and the intensity profile is not entirely Gaussian. The width of the beam at the back-focal plane optimising the *l* = 5*μ*m situation is also noticeble smaller than for the case of *l* = 500nm. This is a direct consequence of Eq. (14). As we increase *l*, the off-diagonal terms in matrix
${\mathbf{\text{R}}}_{12}^{\mathbf{\text{u}}}\hspace{0.17em}(l)$ become less important due to the behaviour of the sinc function as *l* increases. In the limit *l* → ∞ only the diagonal terms remain. This means that the eigenmodes (eigenvectors of
${\mathbf{\text{R}}}_{12}^{\mathbf{\text{u}}}\hspace{0.17em}(l)$) in this case correspond to single plane waves as defined in our initial set. The eigenmode (plane wave) with maximum eigenvalue optimises our problem for an infinitely large ROI. The solution found by the FOEi method for *l* = 100mm (Fig. 6) is the plane wave near the critical angle. As the phase of *a ^{μ}* for zero amplitude |

*a*| is not well defined, it is not displayed in the graph. The result validates that the FOEi method is working correctly, as the near critical angle plane wave provides the highest intensity and force difference at the interface for infinite system.

^{μ}The force difference is significantly increased in the cases of small ROI sizes compared with plane wave illuminated system that optimises the sorting for infinite sorting space (Fig. 7(a)). The bulk of this improvement is due to increased intensity of light in ROI, but the periodic pattern of forces indicates a more complex response of the system. This periodic pattern is not related to the discretization of **k**-space, where we expect periodicity to appear around 14*μ*m.

Naturally, we do not wish the sorting ROI to become too small as the experimental realisation would become increasingly complicated. It is interesting to look at the dependence of Δ*F* on the size of ROI. To show the improvement compared to infinite system, we normalise Δ*F* by Δ*F _{pw}*, where Δ

*F*is the force difference for optimised infinite system (pw stands for plane wave). The result (Fig. 7(b)) indicates that the gain is significant for a wide range of experimentally interesting ROI sizes. The dip around

_{pw}*l*= 14

*μ*m is present due to discretization of

**k**-space described above. The increase in ratio Δ

*F*/Δ

*F*for

_{pw}*l*> 14

*μ*m is then equivalent to the formation of second beam focus in ROI due to onset of periodicity.

*Discussion of results:* The FOEi method is capable of finding the optimal beam shape for illumination such that the force difference is maximised over the whole sorting region. The computationally intensive calculation of matrices **M**
_{1} and **M**
_{2} is compensated by the fact that the same matrices can be used for finding optimal illumination for any size of sorting region. We note that even though the solution does not optimise the vertical force pointing towards the substrate, we found that this is the case for all our solutions. However, the sign of this force is wavelength and particle size dependent and as such the attractive vertical force is not a general feature of the method. Further, the vertical force is not constant in the sorting region, which might introduce some modulation of force difference due to the Faxen correction. To resolve this one may minimise the vertical force and use an auxiliary beam with constant vertical force over the whole sorting region. This would restrict the diffusion of particles in vertical direction in more controlled way. It is possible to modify the method to simultaneously optimise for several parameters of the system. In our case, the full optimised solution for sorting applications of plasmon nanoparticles would involve simultaneous maximalisation of force difference in ROI, minimisation of vertical force, and minimisation of heating. Such a problem reduces to finding matrices (operators) for all parameters of interest and choosing the eigenmodes optimising for such a set of parameters. It is very interesting, for plasmonic sorting in general, to find the beam shape of the field maximising the force difference and minimising the heating. This would also clearly identify the contribution of focusing to the overall improvement of force difference. This is a focus of our ongoing research.

## 5. Conclusion

We successfully optimised the illumination for sorting gold nanoparticles using our two step approach. Firstly, we found the optimal wavelength maximising the nanoparticle separation and minimising the temperature increase in the system. This is an important consideration in plasmonic systems as the excessive heat increases diffusion and convective effects. Secondly, we found the optimum beam shape of the illumination field for sorting using the method of force optical eigenmodes (FOEi). The applicability of the method was numerically demonstrated for the special case of sorting gold nanoparticles of different size. We plan to extend our approach and perform simultaneous optimisation of several parameters of interest for sorting applications, e.g., minimised heating, maximised force difference along substrate and minimised vertical force. This involves finding operators for all of parameters of interest and choosing the eigenmodes optimising them. This will be subject of further work along with the efficient extension of the FOEi method to 3D case.

## Acknowledgments

We thank the UK Engineering and Physical Sciences Research Council for funding, KD is a Royal Society-Wolfson Merit Award Holder.

## References and links

**1. **G. Raschke, S. Kowarik, T. Franzl, C. Sonnichsen, T. A. Klar, J. Feldmann, A. Nichtl, and K. Kurzinger, “Biomolecular recognition based on single gold nanoparticle light scattering,” Nano Lett. **3**, 935–938 (2003). [CrossRef]

**2. **J. Stehr, C. Hrelescu, R. A. Sperling, G. Raschke, M. Wunderlich, A. Nichtl, D. Heindl, K. Kurzinger, W. J. Parak, T. A. Klar, and J. Feldmann, “Gold nanostoves for microsecond dna melting analysis,” Nano Lett. **8**, 619–623 (2008). [CrossRef] [PubMed]

**3. **D. Boyer, P. Tamarat, A. Maali, B. Lounis, and M. Orrit, “Photothermal imaging of nanometer-sized metal particles among scatterers,” Science **297**, 1160–1163 (2002). [CrossRef] [PubMed]

**4. **D. W. Galbraith, M. T. Anderson, and L. A. Herzenberg, “Flow cytometric analysis and facs sorting of cells based on gfp accumulation,” Methods Cell Biol. **58**, 315–341 (1999). [CrossRef] [PubMed]

**5. **T. N. Buican, M. J. Smyth, H. A. Crissman, G. C. Salzman, C. C. Stewart, and J. C. Martin, “Automated single-cell manipulation and sorting by light trapping,” Appl. Opt. **26**, 5311–5316 (1987). [CrossRef] [PubMed]

**6. **K. Grujic, O. G. Helleso, J. P. Hole, and J. S. Wilkinson, “Sorting of polystyrene microspheres using a y-branched optical waveguide,” Opt. Express **13**, 1–7 (2005). [CrossRef] [PubMed]

**7. **S. C. Chapin, V. Germain, and E. R. Dufresne, “Automated trapping, assembly, and sorting with holographic optical tweezers,” Opt. Express **14**, 13095–13100 (2006). [CrossRef] [PubMed]

**8. **P. J. Rodrigo, R. L. Eriksen, V. R. Daria, and J. Gluckstad, “Interactive light-driven and parallel manipulation of inhomogeneous particles,” Opt. Express **10**, 1550–1556 (2002). [PubMed]

**9. **M. Pelton, K. Ladavac, and D. G. Grier, “Transport and fractionation in periodic potential-energy landscapes,” Phys. Rev. E **70**, 031108 (2004). [CrossRef]

**10. **K. Xiao and D. G. Grier, “Multidimensional optical fractionation of colloidal particles with holographic verification,” Phys. Rev. Lett. **104**, 028302 (2010). [CrossRef] [PubMed]

**11. **T. Cizmar, M. Siler, M. Sery, P. Zemanek, V. Garces-Chavez, and K. Dholakia, “Optical sorting and detection of submicrometer objects in a motional standing wave,” Phys. Rev. B **74**, 035105 (2006). [CrossRef]

**12. **I. Ricardez-Vargas, P. Rodriguez-Montero, R. Ramos-Garcia, and K. Volke-Sepulveda, “Modulated optical sieve for sorting of polydisperse microparticles,” Appl. Phys. Lett. **88**, 121116 (2006). [CrossRef]

**13. **L. Paterson, E. Papagiakoumou, G. Milne, V. Garces-Chavez, S. A. Tatarkova, W. Sibbett, F. J. Gunn-Moore, P. E. Bryant, A. C. Riches, and K. Dholakia, “Light-induced cell separation in a tailored optical landscape,” Appl. Phys. Lett. **87**, 123901 (2005). [CrossRef]

**14. **R. F. Marchington, M. Mazilu, S. Kuriakose, V. Garces-Chavez, P. J. Reece, T. F. Krauss, M. Gu, and K. Dholakia, “Optical deflection and sorting of microparticles in a near-field optical geometry,” Opt. Express **16**, 3712–3726 (2008). [CrossRef] [PubMed]

**15. **M. P. MacDonald, G. C. Spalding, and K. Dholakia, “Microfluidic sorting in an optical lattice,” Nature **426**, 421–424 (2003). [CrossRef] [PubMed]

**16. **P. T. Korda, M. B. Taylor, and D. G. Grier, “Kinetically locked-in colloidal transport in an array of optical tweezers,” Phys. Rev. Lett. **89**, 128301 (2002). [CrossRef] [PubMed]

**17. **K. Ladavac, K. Kasza, and D. G. Grier, “Sorting mesoscopic objects with periodic potential landscapes: optical fractionation,” Phys. Rev. E **70**, 010901 (2004). [CrossRef]

**18. **G. Milne, D. Rhodes, M. MacDonald, and K. Dholakia, “Fractionation of polydisperse colloid with acousto-optically generated potential energy landscapes,” Opt. Lett. **32**, 1144–1146 (2007). [CrossRef] [PubMed]

**19. **Y. Y. Sun, L. S. Ong, and X. C. Yuan, “Composite-microlens-array-enabled microfluidic sorting,” Appl. Phys. Lett. **89**, 141108 (2006). [CrossRef]

**20. **M. Mazilu, J. Baumgartl, S. Kosmeier, and K. Dholakia, “Optical eigenmodes; exploiting the quadratic nature of the energy flux and of scattering interactions,” Opt. Express **19**, 933–945 (2011). [CrossRef] [PubMed]

**21. **A. S. Zelenina, R. Quidant, G. Badenes, and M. Nieto-Vesperinas, “Tunable optical sorting and manipulation of nanoparticles via plasmon excitation,” Opt. Lett. **31**, 2054–2056 (2006). [CrossRef] [PubMed]

**22. **R. Quidant, S. Zelenina, and M. Nieto-Vesperinas, “Optical manipulation of plasmonic nanoparticles,” Appl. Phys. A: Mater. Sci. Process. **89**, 233–239 (2007). [CrossRef]

**23. **J. D. Jackson, *Classical Electrodynamics* (John Wiley & Sons, Inc., 1999).

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

**25. **K. Svoboda and S. M. Block, “Biological applications of optical forces,” Annu. Rev. Biophys. Biomol. Struct. **23**, 247–285 (1994). [CrossRef] [PubMed]

**26. **G. Baffou, R. Quidant, and C. Girard, “Heat generation in plasmonic nanostructures: influence of morphology,” Appl. Phys. Lett. **94**, 153109 (2009). [CrossRef]