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
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 , localised heaters  and photo-thermal imaging . 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 , 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 . 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  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)  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 (eiωt), we can write (using summation over repeating indices)
The optical-cycle averaged electromagnetic force in the direction u is given by23] Eq. (2), we can rewrite Eq. (4) as Eq. (6) into Eq. (3) gives 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 u 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 λn of the matrix 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 λn of the matrix M and selecting the largest one from the set extremizes our problem. The eigenvector (force optical eigenmode) corresponding to maximum eigenvalue (given by ), then provides the necessary information about amplitude and phase of incident plane waves in Eq. (1) so that the force is maximised. Experimentally, the optimised can be created using spatial light modulator (offering control of both phase and amplitude) in the system placed in conjugate plane  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 . We first find the solutions E μ for particle p 1. We subsequently use the solutions E μ to find the elements 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  of radius r 1 = 50nm and r 2 = 40nm. We consider a substrate of glass with refractive index ng = 1.5. The particles are assumed to be dispersed in water with nw = 1.33 (Fig. 1).
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 Qsca and absorption Qabs efficiencies for 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 ) 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 asEq. (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  (neglecting particle movement and proximity of glass surface) 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 isFig. 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
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 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 ) 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 ΔFpw, where ΔFpw 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 l = 14μm is present due to discretization of k-space described above. The increase in ratio ΔF/ΔFpw for 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.
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.
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]
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]
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]
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]
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]
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]
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]
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]
26. G. Baffou, R. Quidant, and C. Girard, “Heat generation in plasmonic nanostructures: influence of morphology,” Appl. Phys. Lett. 94, 153109 (2009). [CrossRef]