<?xml version="1.0" encoding="UTF-8"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD Journal Publishing DTD v2.3 20070202//EN" "journalpublishing.dtd">
<article article-type="review-article" dtd-version="2.3" xml:lang="EN" xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">
<front>
<journal-meta>
<journal-id journal-id-type="publisher-id">Front. Nanotechnol.</journal-id>
<journal-title>Frontiers in Nanotechnology</journal-title>
<abbrev-journal-title abbrev-type="pubmed">Front. Nanotechnol.</abbrev-journal-title>
<issn pub-type="epub">2673-3013</issn>
<publisher>
<publisher-name>Frontiers Media S.A.</publisher-name>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="publisher-id">886636</article-id>
<article-id pub-id-type="doi">10.3389/fnano.2022.886636</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Nanotechnology</subject>
<subj-group>
<subject>Review</subject>
</subj-group>
</subj-group>
</article-categories>
<title-group>
<article-title>Plasmonic Nanostructures for Optically Induced Movement</article-title>
<alt-title alt-title-type="left-running-head">Balestrieri et al.</alt-title>
<alt-title alt-title-type="right-running-head">Plasmonic Nanostructures for Induced Thrust</alt-title>
</title-group>
<contrib-group>
<contrib contrib-type="author" corresp="yes">
<name>
<surname>Balestrieri</surname>
<given-names>Sergio</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
<xref ref-type="corresp" rid="c001">&#x2a;</xref>
<uri xlink:href="https://loop.frontiersin.org/people/1699804/overview"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Zito</surname>
<given-names>Gianluigi</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Coppola</surname>
<given-names>Giuseppe</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<xref ref-type="fn" rid="fn1">
<sup>&#x2020;</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/1549484/overview"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Iodice</surname>
<given-names>Mario</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<xref ref-type="fn" rid="fn1">
<sup>&#x2020;</sup>
</xref>
</contrib>
</contrib-group>
<aff id="aff1">
<sup>1</sup>
<institution>Institute of Science Applied and System Intelligence-National Research Council</institution>, <addr-line>Naples</addr-line>, <country>Italy</country>
</aff>
<aff id="aff2">
<sup>2</sup>
<institution>Dipartimento di Fisica</institution>, <institution>Universit&#xe0; di Napoli Federico II</institution>, <addr-line>Napoli</addr-line>, <country>Italy</country>
</aff>
<author-notes>
<fn fn-type="edited-by">
<p>
<bold>Edited by:</bold> <ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/284664/overview">Carlo Rizza</ext-link>, University of L&#x27;Aquila, Italy</p>
</fn>
<fn fn-type="edited-by">
<p>
<bold>Reviewed by:</bold> <ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/1424711/overview">Arindam Dasgupta</ext-link>, Missouri University of Science and Technology, United States</p>
<p>
<ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/1159623/overview">Linhan Lin</ext-link>, Tsinghua University, China</p>
</fn>
<corresp id="c001">&#x2a;Correspondence: Sergio Balestrieri, <email>sergio.balestrieri@na.isasi.cnr.it</email>
</corresp>
<fn fn-type="equal" id="fn1">
<label>
<sup>&#x2020;</sup>
</label>
<p>These authors have contributed equally to this work and share senior authorship</p>
</fn>
<fn fn-type="other">
<p>This article was submitted to Nanophotonics, a section of the journal Frontiers in Nanotechnology</p>
</fn>
</author-notes>
<pub-date pub-type="epub">
<day>05</day>
<month>05</month>
<year>2022</year>
</pub-date>
<pub-date pub-type="collection">
<year>2022</year>
</pub-date>
<volume>4</volume>
<elocation-id>886636</elocation-id>
<history>
<date date-type="received">
<day>28</day>
<month>02</month>
<year>2022</year>
</date>
<date date-type="accepted">
<day>04</day>
<month>04</month>
<year>2022</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#xa9; 2022 Balestrieri, Zito, Coppola and Iodice.</copyright-statement>
<copyright-year>2022</copyright-year>
<copyright-holder>Balestrieri, Zito, Coppola and Iodice</copyright-holder>
<license xlink:href="http://creativecommons.org/licenses/by/4.0/">
<p>This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.</p>
</license>
</permissions>
<abstract>
<p>Optical forces generated at the nanoscale using electric field gradients have proven to be a powerful tool for trapping and moving nano-objects in a variety of application fields ranging from aerospace engineering to biology and medicine. Typically, to achieve this optical effect plasmonic resonant cavities that combine localized surface plasmon resonances and propagative surface plasmon polaritons are used. Indeed, these structures allow to engineer the distribution of the excited field hotspots, so inducing a precise movement of the nanoparticles interacting with the plasmonic field. In this paper, starting from the theoretical analysis of the surface plasmons, the potentialities of plasmonic nanostructures are reviewed, analysing the geometric conformation designed according to the application. The configurations with the most interesting performance, among those mentioned in the literature, are described in detail, examining their main characteristics and limitations. Finally, the future development and prospects of these plasmonic nanostructures are discussed.</p>
</abstract>
<kwd-group>
<kwd>plasmonic structures</kwd>
<kwd>optical force</kwd>
<kwd>photokinetic</kwd>
<kwd>optical propulsion</kwd>
<kwd>nanostructures</kwd>
</kwd-group>
</article-meta>
</front>
<body>
<sec id="s1">
<title>Introduction</title>
<p>The interaction of a light beam with a microscopic object can generate forces and torques, allowing manipulation even if the source of forces is at a considerable distance from the point of application. Each feature of the light (such as amplitude, phase, polarization, spatial distribution, etc) can induce surprising effect on the object manipulation. Absorption and scattering give rise to radiation pressure through the transfer of the linear momentum density of light. The momentum density, in turn, can be influenced by the angular momentum density of light spin (<xref ref-type="bibr" rid="B22">Ruffner and Grier, 2012</xref>) and/or gradients in the phase of the wavefronts (<xref ref-type="bibr" rid="B19">Roichman et al., 2008</xref>). Gradients in light intensity generate dipole forces that can be used for optically trapping a particle in a well-defined position. (<xref ref-type="bibr" rid="B14">Mansuripur, 2014</xref>). Properly structured light beams can transport small objects against the direction of propagation (<xref ref-type="bibr" rid="B2">Chen et al., 2011</xref>; <xref ref-type="bibr" rid="B17">Novitsky et al., 2011</xref>; <xref ref-type="bibr" rid="B32">Wang et al., 2013</xref>).</p>
<p>In this article, we analyze the advantages of generating a force on small objects using surface plasmonic nanostructures that enhance the optical radiation features. Indeed, the design of the surface plasmons, that is the excitation resulting from the resonant interaction in a structure composed of a metal on a dielectric substrate, allows to generate an electromagnetic field with appropriately structured characteristics, such as: concentration of high energies in small areas (beyond the diffraction limit), uniform or non-uniform spatial amplitude distribution, polarisation-dependent behaviour, and so on. Plasmonic structures are capable of inducing forces for the manipulation of small objects; in this review we will focus on the dynamics generated by these structures. In particular, we will observe that nanoparticles placed in proximity to a metal structure capable to excite a plasmonic phenomenon can induce small displacements around an equilibrium point. Moreover, it is also possible to design particular metal structures that generate an appropriate electric field distribution capable of inducing a propulsion. The latter capability could be very useful in new emerging applicative sectors such as medical-biological or space applications. Indeed, with the recent development of micro- and nanosatellites, interest in propulsion driven by a remote optical source has increased. The nanoplasmonic approach could allow generating a micro/nanosatellite propulsion without the need to add power on the satellite, unlike the classical approach (such as chemical or electric propulsion (<xref ref-type="bibr" rid="B6">Gassend et al., 2009</xref>; <xref ref-type="bibr" rid="B24">Sathiyanathan et al., 2011</xref>; <xref ref-type="bibr" rid="B1">Cervone et al., 2016</xref>; <xref ref-type="bibr" rid="B10">Lemmer, 2017</xref>; <xref ref-type="bibr" rid="B30">Tummala and Dutta, 2017</xref>). To provide a self-consistent understanding of the potential of plasmonic nanostructures for optically induced motion, the rest of the paper is organized as follows: in the first section, the theoretical principles of plasmon-enhanced light&#x2013;matter interactions are discussed, highlighting the phenomena to be exploited for force generation. In the second chapter, the optical trapping is reviewed, starting from Askin&#x2019;s method and analyzing the advantages of using plasmonic nanostructures in generate static or dynamic equilibrium. In the third chapter, nanostructures able of generating force fields for long-range particle motion are reported. In addition, some applications and the relative induced thrust parameters, are illustrated. Finally, conclusions and some future perspectives are presented.</p>
</sec>
<sec id="s2">
<title>Theoretical Principles</title>
<p>In this section, we will analyse how, from a theoretical point of view, a force can be generated by an electromagnetic field of an optical radiation and which optically induced phenomena allow the amplification of this force. Given an electromagnetic wave, whose associated fields can be described by means of Maxwell&#x2019;s equations, and an object located inside the light beam, it is possible to define the Maxwell stress tensor as a tensor that encloses all the momentum per unit time resulting from the interaction of the object with the electromagnetic fields, i.e., (<xref ref-type="bibr" rid="B18">Novothny and Hecht, 2006</xref>):<disp-formula id="e1">
<mml:math id="m1">
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="italic">T</mml:mi>
<mml:mrow>
<mml:mi mathvariant="italic">ij</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mi mathvariant="normal">&#x3f5;</mml:mi>
<mml:mn>0</mml:mn>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="italic">E</mml:mi>
<mml:mi mathvariant="italic">i</mml:mi>
</mml:msub>
<mml:msub>
<mml:mi mathvariant="italic">E</mml:mi>
<mml:mi mathvariant="italic">j</mml:mi>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi mathvariant="italic">&#x3b4;</mml:mi>
<mml:mrow>
<mml:mi mathvariant="italic">ij</mml:mi>
</mml:mrow>
</mml:msub>
<mml:msup>
<mml:mrow>
<mml:mrow>
<mml:mo>&#x7c;</mml:mo>
<mml:mrow>
<mml:mover accent="true">
<mml:mi mathvariant="italic">E</mml:mi>
<mml:mo>&#x2192;</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mo>&#x7c;</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x2b;</mml:mo>
<mml:mfrac>
<mml:mn>1</mml:mn>
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="italic">&#x3bc;</mml:mi>
<mml:mn>0</mml:mn>
</mml:msub>
</mml:mrow>
</mml:mfrac>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="italic">B</mml:mi>
<mml:mi mathvariant="italic">i</mml:mi>
</mml:msub>
<mml:msub>
<mml:mi mathvariant="italic">B</mml:mi>
<mml:mi mathvariant="italic">j</mml:mi>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi mathvariant="italic">&#x3b4;</mml:mi>
<mml:mrow>
<mml:mi mathvariant="italic">ij</mml:mi>
</mml:mrow>
</mml:msub>
<mml:msup>
<mml:mrow>
<mml:mrow>
<mml:mo>&#x7c;</mml:mo>
<mml:mrow>
<mml:mover accent="true">
<mml:mi mathvariant="italic">B</mml:mi>
<mml:mo>&#x2192;</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mo>&#x7c;</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(1)</label>
</disp-formula>Where <inline-formula id="inf1">
<mml:math id="m2">
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="italic">E</mml:mi>
<mml:mi mathvariant="italic">i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula id="inf2">
<mml:math id="m3">
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="italic">E</mml:mi>
<mml:mi mathvariant="italic">j</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> are the electric field components; <inline-formula id="inf3">
<mml:math id="m4">
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="italic">B</mml:mi>
<mml:mi mathvariant="italic">i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula id="inf4">
<mml:math id="m5">
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="italic">B</mml:mi>
<mml:mi mathvariant="italic">j</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> are the magnetic induction field components; <inline-formula id="inf5">
<mml:math id="m6">
<mml:mrow>
<mml:mrow>
<mml:mo>&#x7c;</mml:mo>
<mml:mrow>
<mml:mover accent="true">
<mml:mi mathvariant="italic">E</mml:mi>
<mml:mo>&#x2192;</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mo>&#x7c;</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula id="inf6">
<mml:math id="m7">
<mml:mrow>
<mml:mo>&#x7c;</mml:mo>
<mml:mrow>
<mml:mover accent="true">
<mml:mi mathvariant="italic">B</mml:mi>
<mml:mo>&#x2192;</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mo>&#x7c;</mml:mo>
</mml:mrow>
</mml:math>
</inline-formula> are the amplitudes of the electric field and the magnetic induction field, respectively; <inline-formula id="inf7">
<mml:math id="m8">
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="italic">&#x3f5;</mml:mi>
<mml:mn>0</mml:mn>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula id="inf8">
<mml:math id="m9">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3bc;</mml:mi>
<mml:mn>0</mml:mn>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> are the dielectric permittivity in vacuum and the magnetic permeability in vacuum, respectively.</p>
<p>It is possible from Maxwell&#x2019;s stress tensor to derive the force applied to an object resulting from the interaction with electromagnetic fields as described (<xref ref-type="bibr" rid="B18">Novothny and Hecht, 2006</xref>):<disp-formula id="e2">
<mml:math id="m10">
<mml:mrow>
<mml:mrow>
<mml:mo>&#x2329;</mml:mo>
<mml:mover accent="true">
<mml:mi mathvariant="italic">F</mml:mi>
<mml:mo>&#x2192;</mml:mo>
</mml:mover>
<mml:mo>&#x232A;</mml:mo>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:munder>
<mml:mstyle displaystyle="true">
<mml:mo>&#x222b;</mml:mo>
</mml:mstyle>
<mml:mrow>
<mml:mi mathvariant="italic">&#x3b4;V</mml:mi>
</mml:mrow>
</mml:munder>
<mml:mrow>
<mml:mo>&#x2329;</mml:mo>
<mml:mover accent="true">
<mml:mi mathvariant="italic">T</mml:mi>
<mml:mo>&#x21C9;</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mrow>
<mml:mrow>
<mml:mo>(</mml:mo>
</mml:mrow>
<mml:mrow>
<mml:mrow>
<mml:mover accent="true">
<mml:mi mathvariant="italic">r</mml:mi>
<mml:mo>&#x2192;</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mo>,</mml:mo>
<mml:mi mathvariant="italic">t</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x232A;</mml:mo>
</mml:mrow>
<mml:mo>&#x22c5;</mml:mo>
<mml:mrow>
<mml:mover accent="true">
<mml:mi mathvariant="italic">n</mml:mi>
<mml:mo>&#x2192;</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mover accent="true">
<mml:mi mathvariant="italic">r</mml:mi>
<mml:mo>&#x2192;</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mi mathvariant="italic">da</mml:mi>
</mml:mrow>
</mml:math>
<label>(2)</label>
</disp-formula>Where <inline-formula id="inf9">
<mml:math id="m11">
<mml:mo>&#x2329;</mml:mo>
<mml:mover accent="true">
<mml:mi mathvariant="italic">F</mml:mi>
<mml:mo>&#x2192;</mml:mo>
</mml:mover>
<mml:mo>&#x232A;</mml:mo>
</mml:math>
</inline-formula> is the time-averaged force on the object dependent on the electromagnetic field, <inline-formula id="inf10">
<mml:math id="m12">
<mml:mrow>
<mml:mo>&#x2329;</mml:mo>
<mml:mrow>
<mml:mover accent="true">
<mml:mi mathvariant="italic">T</mml:mi>
<mml:mo>&#x21C9;</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mrow>
<mml:mover accent="true">
<mml:mi mathvariant="italic">r</mml:mi>
<mml:mo>&#x2192;</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mo>,</mml:mo>
<mml:mi mathvariant="italic">t</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x232A;</mml:mo>
</mml:mrow>
</mml:math>
</inline-formula> is the time-averaged Maxwell tensor, <inline-formula id="inf11">
<mml:math id="m13">
<mml:mrow>
<mml:mrow>
<mml:mover accent="true">
<mml:mi mathvariant="italic">n</mml:mi>
<mml:mo>&#x2192;</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mover accent="true">
<mml:mrow>
<mml:mi mathvariant="italic">r</mml:mi>
</mml:mrow>
<mml:mo>&#x2192;</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula> is the surface normal and <inline-formula id="inf12">
<mml:math id="m14">
<mml:mrow>
<mml:mi mathvariant="italic">&#x3b4;V</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula> is the surface of the object under consideration.</p>
<p>If the dimensions of the object are much smaller than the wavelength of the incident beam, some approximations (described in (<xref ref-type="bibr" rid="B18">Novothny and Hecht, 2006</xref>)) in the resolution of the Maxwell&#x2019;s tensor can be applied. So the force term reported in equation (<xref ref-type="disp-formula" rid="e2">Eq. 2</xref>) can be redefined in the following way:<disp-formula id="e3">
<mml:math id="m15">
<mml:mrow>
<mml:mrow>
<mml:mo>&#x2329;</mml:mo>
<mml:msup>
<mml:mo>F</mml:mo>
<mml:mo>&#x2192;</mml:mo>
</mml:msup>
<mml:mo>&#x232A;</mml:mo>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:msup>
<mml:mi mathvariant="italic">&#x3b1;</mml:mi>
<mml:mi mathvariant="bold-italic">&#x2032;</mml:mi>
</mml:msup>
</mml:mrow>
<mml:mn>4</mml:mn>
</mml:mfrac>
<mml:mrow>
<mml:msup>
<mml:mo>&#x2207;</mml:mo>
<mml:mo>&#x2192;</mml:mo>
</mml:msup>
</mml:mrow>
<mml:mo>&#xa0;</mml:mo>
<mml:msup>
<mml:mrow>
<mml:mrow>
<mml:mo>&#x7c;</mml:mo>
<mml:mrow>
<mml:mover accent="true">
<mml:mi mathvariant="italic">E</mml:mi>
<mml:mo>&#x2192;</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mo>&#x7c;</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mo>&#x2b;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:msup>
<mml:mi mathvariant="italic">&#x3b1;</mml:mi>
<mml:mrow>
<mml:mi mathvariant="bold-italic">&#x2032;&#x2032;</mml:mi>
</mml:mrow>
</mml:msup>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:mfrac>
<mml:msup>
<mml:mrow>
<mml:mrow>
<mml:mo>&#x7c;</mml:mo>
<mml:mrow>
<mml:mover accent="true">
<mml:mi mathvariant="italic">E</mml:mi>
<mml:mo>&#x2192;</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mo>&#x7c;</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mi mathvariant="bold">&#xa0;</mml:mi>
<mml:mrow>
<mml:msup>
<mml:mo>&#x2207;</mml:mo>
<mml:mo>&#x2192;</mml:mo>
</mml:msup>
</mml:mrow>
<mml:mi mathvariant="italic">&#x3d5;</mml:mi>
</mml:mrow>
</mml:math>
<label>(3)</label>
</disp-formula>Where <inline-formula id="inf13">
<mml:math id="m16">
<mml:mrow>
<mml:msup>
<mml:mi mathvariant="italic">&#x3b1;</mml:mi>
<mml:mi mathvariant="bold-italic">&#x2032;</mml:mi>
</mml:msup>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula id="inf14">
<mml:math id="m17">
<mml:mrow>
<mml:msup>
<mml:mi mathvariant="italic">&#x3b1;</mml:mi>
<mml:mrow>
<mml:mi mathvariant="bold-italic">&#x2032;&#x2032;</mml:mi>
</mml:mrow>
</mml:msup>
</mml:mrow>
</mml:math>
</inline-formula> are the real and imaginary parts of the polarisability, respectively; <inline-formula id="inf15">
<mml:math id="m18">
<mml:mi>&#x3d5;</mml:mi>
</mml:math>
</inline-formula> is the phase function relative to the electric field, and <inline-formula id="inf16">
<mml:math id="m19">
<mml:mrow>
<mml:mo>&#x7c;</mml:mo>
<mml:mrow>
<mml:mover accent="true">
<mml:mi mathvariant="italic">E</mml:mi>
<mml:mo>&#x2192;</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mo>&#x7c;</mml:mo>
</mml:mrow>
</mml:math>
</inline-formula> is the electric field amplitude. From this formula the force induced by an electromagnetic field on microscopic objects is composed of two terms. The former, called gradient force, depends on the variation of the electric field amplitude and the dispersive part of complex polarizability. The second term, called radiation pressure, is related to the phase variations of the electric field and the dissipative part of the complex polarizability. In general, the gradient force tends to accelerate the particles orthogonally to the direction of propagation of the incident electromagnetic field, whereas the radiation pressure tends to accelerate the particles along the propagation direction of the incident field (<xref ref-type="bibr" rid="B18">Novothny and Hecht, 2006</xref>). This pressure is formed only in particular configurations of the particle (e.g., metal particle) as it is also related to dissipation.</p>
<p>Because of these behaviours, gradient forces are commonly used in the literature to develop structures capable of inducing motion in particles, and this means designing structures capable of generating fields of non-uniform amplitude. Among the various approaches described in the literature (<xref ref-type="bibr" rid="B28">Soukoulis, 2002</xref>; <xref ref-type="bibr" rid="B5">Fujimaki et al., 2008</xref>) to generate electric fields with a spatial distribution of non-uniform amplitude, one of the most effective is surface plasmonics. This phenomenon occurs when an electromagnetic wave irradiates an interface between a metallic material and a dielectric material (<xref ref-type="bibr" rid="B13">Maier, 2007</xref>). Surface plasmonics manifests itself in two different ways called Surface Plasmons Polaritons (SPP) and Localized Surface Plasmon (LSP).</p>
<p>This phenomenon related to coherent oscillations of delocalised electrons occurs when an electromagnetic wave radiates an interface between a metallic material and a dielectric material. In fact, when light and freely moving electrons (as in a metal) interact, the oscillating electric-magnetic field of the light forces the free electrons to move together with the electromagnetic oscillations. This movement of electrons produces an additional electromagnetic field. Two types of surface plasmonics can be considered: Surface Plasmons Polaritons (SPP) and Localized Surface Plasmons (LSP).</p>
<p>SPPs are waves that propagate along the metal-dielectric interface and are highly localised, indeed the electric field in the direction orthogonal to the interface is evanescent. To obtain an SPP analytically, the Maxwell&#x2019;s equations on an interface between a metal (<inline-formula id="inf17">
<mml:math id="m20">
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="normal">&#x3f5;</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>) and a dielectric material (<inline-formula id="inf18">
<mml:math id="m21">
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="normal">&#x3f5;</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>) defined along an xz plane has to be considered. Maxwell&#x2019;s equations with continuity conditions for electromagnetic fields in TM polarization return a propagating wave with a <inline-formula id="inf19">
<mml:math id="m22">
<mml:mrow>
<mml:mi mathvariant="italic">&#x3b2;</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula> constant defined as:<disp-formula id="e4">
<mml:math id="m23">
<mml:mrow>
<mml:mrow>
<mml:mi mathvariant="italic">&#x3b2;</mml:mi>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="italic">k</mml:mi>
<mml:mn>0</mml:mn>
</mml:msub>
</mml:mrow>
<mml:msqrt>
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="italic">&#x3f5;</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
<mml:msub>
<mml:mi mathvariant="italic">&#x3f5;</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="italic">&#x3f5;</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi mathvariant="italic">&#x3f5;</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:msqrt>
</mml:mrow>
</mml:math>
<label>(4)</label>
</disp-formula>where <inline-formula id="inf20">
<mml:math id="m24">
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="italic">k</mml:mi>
<mml:mn>0</mml:mn>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the propagation constant in vacuum. On the other hand, the Maxwell&#x2019;s equations with continuity conditions for the fields defined in TE polarization doesn&#x2019;t return any solutions. Thus, SPPs are only obtainable in TM polarization.</p>
<p>In order to generate an SPP, the x-long component of the momentum of the free photon (<inline-formula id="inf21">
<mml:math id="m25">
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="italic">k</mml:mi>
<mml:mrow>
<mml:msub>
<mml:mn>0</mml:mn>
<mml:mi>x</mml:mi>
</mml:msub>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>) has to coincide with the propagation constant of the SPP. However, in formula (<xref ref-type="disp-formula" rid="e4">Eq. 4</xref>) we see that <inline-formula id="inf22">
<mml:math id="m26">
<mml:mrow>
<mml:mi mathvariant="italic">&#x3b2;</mml:mi>
<mml:mo>&#x3e;</mml:mo>
<mml:msub>
<mml:mi mathvariant="italic">k</mml:mi>
<mml:mrow>
<mml:msub>
<mml:mn>0</mml:mn>
<mml:mi>x</mml:mi>
</mml:msub>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>. This implies that an SPP can only exist if there is a momentum mismatch between its momentum and that of free space, otherwise the wave cannot be confined within the surface. Thus, SPP can never be generated by direct illumination, but specific excitation techniques must be used, such as (<xref ref-type="bibr" rid="B13">Maier, 2007</xref>): prism, grating or near field light source.</p>
<p>On the other hand, Localized Surface Plasmon are non-propagating excitations of conduction electrons relative to metallic nanostructures coupled to the electromagnetic field. LSPs are derived from Mie theory, assuming that the nanostructure has spatial dimensions much smaller than the wavelength of the incident beam (quasi-static field approximation) (<xref ref-type="bibr" rid="B13">Maier, 2007</xref>). This assumption implies that the equations of electrostatics, i.e., Laplace&#x2019;s equations, can be used in the evaluation of the electric field. Once solved, a polarizability in the following form is obtained:<disp-formula id="e6">
<mml:math id="m27">
<mml:mrow>
<mml:mi mathvariant="italic">&#x3b1;</mml:mi>
<mml:mo>&#x2243;</mml:mo>
<mml:mi mathvariant="italic">V</mml:mi>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="italic">&#x3f5;</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi mathvariant="italic">&#x3f5;</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="italic">&#x3f5;</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:mn>2</mml:mn>
<mml:msub>
<mml:mi mathvariant="italic">&#x3f5;</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:math>
<label>(6)</label>
</disp-formula>Where V is the volume of the nanostructure. Thus, for polarizability can be determined a resonance (<inline-formula id="inf23">
<mml:math id="m28">
<mml:mrow>
<mml:mi mathvariant="italic">&#x3b1;</mml:mi>
<mml:mo>&#x2192;</mml:mo>
<mml:mi>&#x221e;</mml:mi>
<mml:mo>.</mml:mo>
</mml:mrow>
</mml:math>
</inline-formula>) at a given frequency of the incident field by the Froehlich condition:<disp-formula id="e7">
<mml:math id="m29">
<mml:mrow>
<mml:mo>&#x211c;</mml:mo>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="italic">&#x3f5;</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi mathvariant="italic">&#x3c9;</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>2</mml:mn>
<mml:msub>
<mml:mi mathvariant="italic">&#x3f5;</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
</mml:mrow>
</mml:math>
<label>(7)</label>
</disp-formula>
</p>
<p>Thus, at resonance there is an amplification of electromagnetic fields.</p>
<p>In addition, unlike SPPs, there are no conditions limiting the generation of a LSP, so direct illumination can also be used.</p>
</sec>
<sec id="s3">
<title>Plasmonic Force Confinement</title>
<p>In the previous chapter, it was explained from a theoretical point of view how a force can be generated from an electromagnetic field, and that in order to obtain intense forces, the mechanisms of surface plasmonics can be used, which allow electric fields with high amplitudes and non-uniform spatial distribution to be achieved. This chapter describes how surface plasmonics phenomena can be induced by nanostructures suitably designed to achieve the desired electric field distributions.</p>
<p>A typical example are optical tweezers, structures capable of trapping particles and keeping them in stable equilibrium at a given location in space. Historically, optical tweezers were proposed by Askin (<xref ref-type="bibr" rid="B14">Mansuripur, 2014</xref>) in 1970. Their operation in a macroscopic case is defined as follows: A spherical particle of dielectric material (index of refraction <inline-formula id="inf24">
<mml:math id="m30">
<mml:mrow>
<mml:msub>
<mml:mi>n</mml:mi>
<mml:mi>b</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>) is immersed in a liquid (refractive index <inline-formula id="inf25">
<mml:math id="m31">
<mml:mrow>
<mml:msub>
<mml:mi>n</mml:mi>
<mml:mi>l</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>). A laser is focused on the dielectric particle at a point far from the centre of the sphere. Due to the transmission effect, the momentum of the incoming beam (<inline-formula id="inf26">
<mml:math id="m32">
<mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mover>
<mml:mi mathvariant="bold-italic">p</mml:mi>
<mml:mo stretchy="true">&#x2192;</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mi>n</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula>) will have a different direction than the momentum of the outgoing beam (<inline-formula id="inf27">
<mml:math id="m33">
<mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mover>
<mml:mi mathvariant="bold-italic">p</mml:mi>
<mml:mo stretchy="true">&#x2192;</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mrow>
<mml:mi>o</mml:mi>
<mml:mi>u</mml:mi>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula>); this implies the formation of (<inline-formula id="inf28">
<mml:math id="m34">
<mml:mrow>
<mml:mi mathvariant="italic">&#x394;</mml:mi>
<mml:mover accent="true">
<mml:mi mathvariant="bold-italic">p</mml:mi>
<mml:mo>&#x2192;</mml:mo>
</mml:mover>
<mml:mo>&#x3d;</mml:mo>
<mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mover>
<mml:mi mathvariant="bold-italic">p</mml:mi>
<mml:mo stretchy="true">&#x2192;</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mrow>
<mml:mi>o</mml:mi>
<mml:mi>u</mml:mi>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mover>
<mml:mi mathvariant="bold-italic">p</mml:mi>
<mml:mo stretchy="true">&#x2192;</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mi>n</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula>) such that the force obtained from this <inline-formula id="inf29">
<mml:math id="m35">
<mml:mrow>
<mml:mi mathvariant="italic">&#x394;</mml:mi>
<mml:mover accent="true">
<mml:mi mathvariant="bold-italic">p</mml:mi>
<mml:mo>&#x2192;</mml:mo>
</mml:mover>
</mml:mrow>
</mml:math>
</inline-formula> allows the particle to rise balancing with the force of gravity (<xref ref-type="fig" rid="F1">Figure 1</xref>).</p>
<fig id="F1" position="float">
<label>FIGURE 1</label>
<caption>
<p>Optical trapping method used by Askin.</p>
</caption>
<graphic xlink:href="fnano-04-886636-g001.tif"/>
</fig>
<p>Moreover, if the particle is moved sideways relative to the equilibrium point, a pullback force is generated that returns the particle to its original position, making the equilibrium stable.</p>
<p>If the particle size is reduced to the order of microns or less, Askin&#x2019;s method is unable to maintain the stable equilibrium of the particle due to effects that disturb its effectiveness (<xref ref-type="bibr" rid="B9">Juan et al., 2011</xref>); in particular:<list list-type="simple">
<list-item>
<p>&#x2022; The reduction of the pullback force related to its dependence on the volume of the particle.</p>
</list-item>
<list-item>
<p>&#x2022; The reduction in viscous friction due to the dependence on the surface area of the particle.</p>
</list-item>
<list-item>
<p>&#x2022; The effect of temperature, which rises when the particle is microscopic, and thus generating chaotic motions around the equilibrium point, delocalising the particle.</p>
</list-item>
</list>
</p>
<p>An increase in the intensity of the incident laser beam can be used to mitigate the delocalization effect and other microscopic effects.</p>
<p>The main limitation of this technique is related to the size of the particles. When this latter reaches a value close to half the wavelength of the incident beam, diffractive effects appear which tend to disperse the beam of light and prevent it from focusing on the particle.</p>
<p>There are some techniques that can overcome the limitations associated with diffractive effects. For example, optoelectronic tweezers (<xref ref-type="bibr" rid="B33">Wu, 2007</xref>), in which light creates virtual electrodes and therefore, by coupling with a non-uniform electric field, generates a region in which the particles are confined. A similar effect is achieved with the opto-thermophoretic tweezers technique (<xref ref-type="bibr" rid="B11">Zheng et al., 2018</xref>), which creates a temperature gradient by optical heating, allowing for manipulation of displacement and trapping of particles.</p>
<p>These techniques allow to have a good control on micro and nanoparticles, even if the trapping will be extremely weak; on the other hand, Surface Plasmons Polaritons (<xref ref-type="bibr" rid="B26">Shuller et al., 2010</xref>; <xref ref-type="bibr" rid="B9">Juan et al., 2011</xref>; <xref ref-type="bibr" rid="B3">Daly et al., 2015</xref>) generate extremely concentrated fields beyond the diffractive limit, succeeding to trap a nanoparticle with extreme efficiency, however, making some dynamic operations more complex than in previous cases, in which they were easily performed.</p>
<p>In general, when a plasmonic optical tweezer is irradiated by an electromagnetic wave at resonance frequency, a symmetrical field distribution (hotspot) is generated in a specific area of the structure such that the resulting forces tend to trap the nanoparticles in the centre of the hotspot. A typical example of a plasmonic optical tweezer is shown in (<xref ref-type="bibr" rid="B34">Yang et al., 2011</xref>) (<xref ref-type="fig" rid="F2">Figure 2</xref>), where a cylindrical silicon waveguide is placed at a distance (<italic>g</italic>) from a metallic substrate.</p>
<fig id="F2" position="float">
<label>FIGURE 2</label>
<caption>
<p>Image showing the silicon waveguide separated by a sub-micrometer gap to a metal substrate, showing the optical force applied to the nanoparticles (Complete with permission from <xref ref-type="bibr" rid="B34">Yang et al. (2011)</xref>. Copyright 2011 American Chemical Society).</p>
</caption>
<graphic xlink:href="fnano-04-886636-g002.tif"/>
</fig>
<p>When the system is illuminated from the <italic>z</italic> axis, an electric field hotspot is generated in the space between the guide and the metal substrate, and the associated force is such that the nanoparticles are trapped in this gap.</p>
<p>In addition to generating symmetrical field distributions that tend to lock the nanoparticles at a defined point, plasmonic optical tweezers can also be designed to achieve equilibrium dynamically. That is, generating a stationary movement around a given point. An example of this type of structure is described in (<xref ref-type="bibr" rid="B7">Huft et al., 2017</xref>) and illustrated in <xref ref-type="fig" rid="F3">Figure 3A</xref>, where an optical tweezer is constructed using a plasmonic bull&#x2019;s eye configuration to transform the incident wave into converging plasmonic waves.</p>
<fig id="F3" position="float">
<label>FIGURE 3</label>
<caption>
<p>
<bold>(A)</bold> the plasmonic bull&#x2019;s eye with the direction of the plasmonic waves used to trap the particle in the centre and some trajectories that can be generated with time-dependent beams; <bold>(B)</bold> centred hot spot simulated with equation (<xref ref-type="disp-formula" rid="e9">Eq. 9</xref>); <bold>(C)</bold> shifted phase profile hot spot simulated with equation (<xref ref-type="disp-formula" rid="e9">Eq. 9</xref>); <bold>(D)</bold> Normalized Force plot for a Bull&#x2019;s Eye device and a conventional optical tweezers for x-component Force.; <bold>(E)</bold> Normalized Force plot for a Bull&#x2019;s Eye device and a conventional optical tweezers for z-component force. (Complete with permission from <xref ref-type="bibr" rid="B7">Huft et al. (2017)</xref>. Copyright 2017 American Chemical Society).</p>
</caption>
<graphic xlink:href="fnano-04-886636-g003.tif"/>
</fig>
<p>In order for the particles to be trapped in the centre of the structure, the plasmonic waves have to interfere constructively in the centre and to achieve this effect, a phase shift of &#x3c0; is required in the generation of the various SPPs. This implies that it is necessary to modulate the beam in such a way that the interference of the various SPPs occurs at different positions from the centre of the structure. To achieve this, the phase function (<inline-formula id="inf30">
<mml:math id="m36">
<mml:mi>&#x3d5;</mml:mi>
</mml:math>
</inline-formula>) of the incident beam must depend on both the polar angle of the bull&#x2019;s eye (<inline-formula id="inf31">
<mml:math id="m37">
<mml:mrow>
<mml:mi>&#x3b8;</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:math>
</inline-formula> and the spatial point where the particle is to be trapped (<italic>x</italic>, <italic>y</italic>), i.e.<disp-formula id="e8">
<mml:math id="m38">
<mml:mrow>
<mml:mi mathvariant="italic">&#x3d5;</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi mathvariant="italic">x</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi mathvariant="italic">y</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi mathvariant="italic">&#x3b8;</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mrow>
<mml:mo>{</mml:mo>
<mml:mrow>
<mml:mi mathvariant="italic">&#xa0;</mml:mi>
<mml:mtable>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:mn>2</mml:mn>
<mml:mi mathvariant="italic">&#x3c0;</mml:mi>
<mml:msqrt>
<mml:mrow>
<mml:msup>
<mml:mrow>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi mathvariant="normal">cos</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi mathvariant="italic">&#x3b8;</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x2b;</mml:mo>
<mml:mfrac>
<mml:mi mathvariant="italic">x</mml:mi>
<mml:mi mathvariant="italic">a</mml:mi>
</mml:mfrac>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mo>&#x2b;</mml:mo>
<mml:msup>
<mml:mrow>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi mathvariant="normal">sin</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi mathvariant="italic">&#x3b8;</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x2b;</mml:mo>
<mml:mfrac>
<mml:mi mathvariant="italic">y</mml:mi>
<mml:mi mathvariant="italic">a</mml:mi>
</mml:mfrac>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mi mathvariant="italic">&#xa0;</mml:mi>
</mml:mrow>
</mml:msqrt>
</mml:mrow>
</mml:mtd>
<mml:mtd>
<mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:mfrac>
<mml:mi mathvariant="italic">&#x3c0;</mml:mi>
<mml:mn>2</mml:mn>
</mml:mfrac>
<mml:mo>&#x3c;</mml:mo>
<mml:mi mathvariant="italic">&#x3b8;</mml:mi>
<mml:mo>&#x3c;</mml:mo>
<mml:mfrac>
<mml:mi mathvariant="italic">&#x3c0;</mml:mi>
<mml:mn>2</mml:mn>
</mml:mfrac>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:mi mathvariant="italic">&#x3c0;</mml:mi>
<mml:mo>&#x2b;</mml:mo>
<mml:mn>2</mml:mn>
<mml:mi mathvariant="italic">&#x3c0;</mml:mi>
<mml:msqrt>
<mml:mrow>
<mml:msup>
<mml:mrow>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi mathvariant="normal">cos</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi mathvariant="italic">&#x3b8;</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x2b;</mml:mo>
<mml:mfrac>
<mml:mi mathvariant="italic">x</mml:mi>
<mml:mi mathvariant="italic">a</mml:mi>
</mml:mfrac>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mo>&#x2b;</mml:mo>
<mml:msup>
<mml:mrow>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi mathvariant="normal">sin</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi mathvariant="italic">&#x3b8;</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x2b;</mml:mo>
<mml:mfrac>
<mml:mi mathvariant="italic">y</mml:mi>
<mml:mi mathvariant="italic">a</mml:mi>
</mml:mfrac>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mi mathvariant="bold-italic">&#xa0;</mml:mi>
</mml:mrow>
</mml:msqrt>
</mml:mrow>
</mml:mtd>
<mml:mtd>
<mml:mrow>
<mml:mfrac>
<mml:mi mathvariant="italic">&#x3c0;</mml:mi>
<mml:mn>2</mml:mn>
</mml:mfrac>
<mml:mo>&#x3c;</mml:mo>
<mml:mi mathvariant="italic">&#x3b8;</mml:mi>
<mml:mo>&#x3c;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mn>3</mml:mn>
<mml:mi mathvariant="italic">&#x3c0;</mml:mi>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:mfrac>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(8)</label>
</disp-formula>Where <italic>a</italic> is the period of the Bull&#x2019;s eye pattern which is chosen equal to the wavelength of the SPP at the interface between silver and water. To verify this behaviour, the authors carried out some simulations using COMSOL Multiphysics (<xref ref-type="bibr" rid="B7">Huft et al., 2017</xref>) in which the plasmonic ring was modelled assuming a set of sources along the ring with a field defined as:<disp-formula id="e9">
<mml:math id="m39">
<mml:mrow>
<mml:msubsup>
<mml:mi mathvariant="italic">E</mml:mi>
<mml:mi mathvariant="italic">z</mml:mi>
<mml:mrow>
<mml:mi mathvariant="italic">sp</mml:mi>
</mml:mrow>
</mml:msubsup>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi mathvariant="italic">&#x3c1;</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi mathvariant="italic">&#x3b8;</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:msup>
<mml:mi mathvariant="italic">e</mml:mi>
<mml:mrow>
<mml:mi mathvariant="italic">i</mml:mi>
<mml:msub>
<mml:mi mathvariant="italic">k</mml:mi>
<mml:mrow>
<mml:mi mathvariant="italic">sp</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mi mathvariant="italic">&#x3c1;&#xa0;</mml:mi>
</mml:mrow>
</mml:msup>
</mml:mrow>
<mml:mrow>
<mml:msqrt>
<mml:mi mathvariant="italic">&#x3c1;</mml:mi>
</mml:msqrt>
</mml:mrow>
</mml:mfrac>
<mml:msup>
<mml:mi mathvariant="italic">e</mml:mi>
<mml:mrow>
<mml:mi mathvariant="italic">i</mml:mi>
<mml:mi>&#x3d5;</mml:mi>
</mml:mrow>
</mml:msup>
<mml:mi mathvariant="normal">cos</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi mathvariant="italic">&#x3b8;</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(9)</label>
</disp-formula>Where <inline-formula id="inf32">
<mml:math id="m40">
<mml:mi>&#x3c1;</mml:mi>
</mml:math>
</inline-formula> is the distance from the source point, <inline-formula id="inf33">
<mml:math id="m41">
<mml:mrow>
<mml:msub>
<mml:mi>k</mml:mi>
<mml:mrow>
<mml:mi>s</mml:mi>
<mml:mi>p</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the wave number of the SPP and <inline-formula id="inf34">
<mml:math id="m42">
<mml:mrow>
<mml:mi>&#x3b8;</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>&#x3d5;</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula> are the same as in <xref ref-type="disp-formula" rid="e8">Eq. 8</xref>.</p>
<p>The simulations show that hot-spots (<xref ref-type="fig" rid="F3">Figures 3B,C</xref>) are formed whose position depends on the phase of the incedent beam, confirming the assumptions defined above. <xref ref-type="fig" rid="F3">Figure 3D</xref> shows the difference between the x-long components of the resulting forces in a Bull&#x2019;eye device case (orange) compared to a conventional optical tweezer (blue). It can be seen that the force generated by the Bull&#x2019;s Eye is much more intense than in the conventional case; the same result is also obtained if we compare (see <xref ref-type="fig" rid="F3">Figure 3E</xref>) the components along z of the forces of the Bull&#x27;Eye (orange) with respect to the conventional case (blue). For this reason, the plasmonic Bull&#x2019;s Eye has more effective trapping than the conventional case.</p>
<p>Moreover, the results shown that using a beam with a time-varying phase, the trapping site moves, generating a motion of the nanoparticle. The motions obtained are such as to generate either linear trajectories along a specific direction, or circular trajectories (some examples are shown in <xref ref-type="fig" rid="F3">Figure 3A</xref>).</p>
<p>
<xref ref-type="bibr" rid="B12">Liaw et al. (2018)</xref> proposed another particular case of a plasmonic optical tweezer with dynamic equilibrium obtained using an array of gold dimers separated by a gap (see <xref ref-type="fig" rid="F4">Figure 4</xref>). These are irradiated by a polarised field along x, generating forces and relative moments on polystyrene nanoparticles in the vicinity of the system.</p>
<fig id="F4" position="float">
<label>FIGURE 4</label>
<caption>
<p>Configuration of a 1D gold dimer array trapping a nanoparticle using a normally incident Gaussian beam. The dimer gap, lattice constant and displacement along y are denoted respectively by <inline-formula id="inf35">
<mml:math id="m43">
<mml:mrow>
<mml:mi mathvariant="italic">d</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi mathvariant="italic">c</mml:mi>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi mathvariant="italic">D</mml:mi>
<mml:mi mathvariant="italic">y</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>. The Gaussian beam is polarised along x and has parameters <inline-formula id="inf36">
<mml:math id="m44">
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="italic">h</mml:mi>
<mml:mi mathvariant="italic">f</mml:mi>
</mml:msub>
<mml:mi mathvariant="italic">&#xa0;e&#xa0;</mml:mi>
<mml:msub>
<mml:mi mathvariant="italic">w</mml:mi>
<mml:mn>0</mml:mn>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>. [Reproduced from <xref ref-type="bibr" rid="B12">Liaw et al. (2018)</xref>].</p>
</caption>
<graphic xlink:href="fnano-04-886636-g004.tif"/>
</fig>
<p>In <xref ref-type="fig" rid="F5">Figure 5A</xref> can be seen the optical force generated by the system; to achieve the equilibrium this optical force is compensated by the friction force generated by the liquid in which the nanoparticles are immersed. Two types of equilibrium can be generated: a contact equilibrium in which the nanoparticles tend to adhere to a gold dimer, and a non-contact equilibrium in which the dimers generate a swirling motion around a point in space known as the stagnation point. The resulting moments are the combination of the angular momentum generated by the Maxwell tensor and the momentum generated by friction with the liquid. In <xref ref-type="fig" rid="F5">Figure 5B</xref> illustrates the generation of two nanoparticle spins, the direction of which depends on the vortex&#x2019;s torsion. In particular, if the vortex rotates clockwise (<xref ref-type="fig" rid="F5">Figure 5A</xref>), there will be an anti-clockwise nanoparticle spin (<xref ref-type="fig" rid="F5">Figure 5B</xref>), whereas if the vortex rotates anti-clockwise (<xref ref-type="fig" rid="F5">Figure 5A</xref>), there will be a clockwise nanoparticle spin (<xref ref-type="fig" rid="F5">Figure 5B</xref>). Moreover, the results show that if the dimer array is moved away from the focal point of the beam, the nanoparticles initially tend to follow the array (<xref ref-type="bibr" rid="B12">Liaw et al., 2018</xref>). However, there is a threshold length which, if exceeded, generates an additional vortex leading to trapping of the nanoparticles away from the dimeric matrix. Finally, since these structures are plasmonic in nature, it is also possible to obtain an increase in vortex intensity by reducing the distance or gap between the dimers.</p>
<fig id="F5" position="float">
<label>FIGURE 5</label>
<caption>
<p>
<bold>(A)</bold> Optical vortices manifested in a streamline map of the optical force field. In particular, this streamline map depicts the field of forces (<inline-formula id="inf37">
<mml:math id="m45">
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="italic">F</mml:mi>
<mml:mi mathvariant="italic">y</mml:mi>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi mathvariant="italic">F</mml:mi>
<mml:mi mathvariant="italic">z</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>) applied to a polystyrene nanoparticle centred at different places in space (<inline-formula id="inf38">
<mml:math id="m46">
<mml:mrow>
<mml:mi mathvariant="italic">y</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi mathvariant="italic">z</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula>) in the yz plane at x &#x3d; 0. The arrows and colours represent the direction and magnitude of the optical force respectively. The dotted lines represent the physical limit of the polystyrene nanoparticle considered here with a radius of 100&#xa0;nm. From this image we therefore have evidence of the presence of both contact and non-contact modes for optical trapping. The stagnation points of the non-contact modes are shown at <inline-formula id="inf39">
<mml:math id="m47">
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="italic">x</mml:mi>
<mml:mn>0</mml:mn>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mn>0</mml:mn>
<mml:mo>,</mml:mo>
<mml:mo>&#xb1;</mml:mo>
<mml:mn>361</mml:mn>
<mml:mi mathvariant="italic">nm</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi mathvariant="bold-italic">&#xa0;</mml:mi>
<mml:mn>207</mml:mn>
<mml:mi mathvariant="italic">nm</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula>. <bold>(B)</bold> Map of nanoparticle spin moments. The colouring represents the intensity of the moment. [Adapted from <xref ref-type="bibr" rid="B12">Liaw et al. (2018)</xref>].</p>
</caption>
<graphic xlink:href="fnano-04-886636-g005.tif"/>
</fig>
</sec>
<sec id="s4">
<title>Plasmonic Optical Thruster</title>
<p>In the previous paragraph it was reviewed how optical tweezers can produce symmetrical field distributions such that the resulting forces are able to keep a nanoparticle in equilibrium, while in the following a number of proposals for achieving a long-range movement of nanoparticle will be presented. There is the possibility of obtaining thrust forces using radiation-matter interaction effects; in particular, we mention solar sail (<xref ref-type="bibr" rid="B29">Swartzlander, 2017</xref>), (<xref ref-type="bibr" rid="B36">Davis and Mead, 2011</xref>), laser ablation (<xref ref-type="bibr" rid="B23">Sasoh, et al., 2017</xref>) and the thermo-optical effect observed in (<xref ref-type="bibr" rid="B16">Ndukaife et al., 2015</xref>). In this paragraph, however, we will observe that it is possible to induce a long-range force in microscopic objects by means of a gradient force induced by specific field distributions. It will also be shown that this effect is achieved very effectively thanks to surface plasmonic phenomenon.</p>
<p>A simple plasmonic device has been presented in (<xref ref-type="bibr" rid="B31">Wang et al., 2009</xref>), consisting of a metal layer placed on the base of a prism with the tip pointing downwards (this type of solution is also known as the Kreshmann configuration (<xref ref-type="bibr" rid="B13">Maier, 2007</xref>) see <xref ref-type="fig" rid="F6">Figure 6</xref>). Gold nanoparticles are placed on this structure in aqueous solution. If the particles are free from any constraint they tend to move only for thermal effects. These effects generate a chaotic motion (Brownian motion) in all directions of space. When an incident beam illuminates the prism at a certain angle, nanoparticles tend to acquire a drift velocity in the direction parallel to the air-metal interface (<italic>x</italic>-axis) and in the direction orthogonal to it (<italic>z</italic>-axis); this is sketched in <xref ref-type="fig" rid="F6">Figure 6</xref>.</p>
<fig id="F6" position="float">
<label>FIGURE 6</label>
<caption>
<p>Interface of plasmonic coupling with gold nanoparticle (xz plane) (Complete with permission from <xref ref-type="bibr" rid="B31">Wang et al. (2009)</xref>. Copyright 2009 American Chemical Society).</p>
</caption>
<graphic xlink:href="fnano-04-886636-g006.tif"/>
</fig>
<p>However, the drift only occurs if the incident beam has a TM polarization. This suggests that the cause of the drift phenomenon is attributable to the formation of an SPP. The hypothesis is confirmed by the contextual observation of a chaotic motion along the <italic>y</italic>-axis. In fact, the SPPs are such that they generate field distributions only along the interface and orthogonal to the interface (see <xref ref-type="fig" rid="F6">Figure 6</xref>) so this implies that forces are obtained only along the two aforesaid directions, leaving the <italic>y</italic>-axis free from any constraint. <xref ref-type="fig" rid="F7">Figure 7</xref> presents the results of the analysis of the optical forces (and the induced velocity) acting on the nanoparticles, as function either of the distance between the metal surface and the aqueous solution and the intensity of the optical beam.</p>
<fig id="F7" position="float">
<label>FIGURE 7</label>
<caption>
<p>
<bold>(A)</bold> The components of the optical force applied to the nanoparticle as a function of the distance between the particle and the metal substrate. <bold>(B)</bold> The velocity of the particle measured as the intensity of the incident beam changes. Complete with permission from <xref ref-type="bibr" rid="B31">Wang et al. (2009)</xref>. Copyright 2009 American Chemical Society.</p>
</caption>
<graphic xlink:href="fnano-04-886636-g007.tif"/>
</fig>
<p>In particular, it has been observed that as the intensity of the incident beam increases, the drift speed tends to increase, implying that the optical force also increases as the intensity increases. On the other hand, as the distance between the metal and the nanostructures increases, the force tends to decrease, with different behaviour depending on the force component observed. In fact, the force along z tends to decay exponentially as the distance varies, while that along x decays in an almost linear manner. Finally, the authors performed a comparison between the propulsion generation obtained through the formation of an SPP, and that obtained through total internal reflection, i.e., using the same structure but without the metal layer. The results show that the fields obtained by SPP are much more intense, although they can only be used if the excitation beam illuminates the structure at the right resonance angle.</p>
<p>A more complex configuration capable of generating propulsion with an output speed that can be used for medical-biological applications is well explained in (<xref ref-type="bibr" rid="B27">S&#xf8;ndergaard et al., 2010</xref>; <xref ref-type="bibr" rid="B25">Shalin and Sukhov, 2013</xref>). The structure, as can be seen in <xref ref-type="fig" rid="F8">Figure 8</xref>, consists of a V-shape cavity made of a gold film deposited on a glass substrate.</p>
<fig id="F8" position="float">
<label>FIGURE 8</label>
<caption>
<p>Sketch of the V-shape plasmonic structure, illuminated by a TM-polarized light. Complete with permission of <xref ref-type="bibr" rid="B27">S&#xf8;ndergaard et al. (2010)</xref>. Copyright 2010 American Chemical Society.</p>
</caption>
<graphic xlink:href="fnano-04-886636-g008.tif"/>
</fig>
<p>When the structure is irradiated with TM-polarized light, an intense electric field is generated at the bottom of the trench; this is due to the generation of a resonant LSP. Once the resonance is excited, two SPPs are formed along the edges of the structure by near field excitation. The coupling of the two Surface Plasmon implies a significant amplification of the field. The amplification factor is a function of the geometrical parameters and in particular an increase can be obtained if the angle (<inline-formula id="inf40">
<mml:math id="m48">
<mml:mi>&#x3b8;</mml:mi>
</mml:math>
</inline-formula>) of the V-shape is reduced. In fact, when the two edges of the V-shaped trench approach each other, the interference effects between the two SPPs are amplified. Even the resonance is highly dependent on the depth (h) and angle (&#x3b8;) of the V-shape. In particular, the wavelength of the resonance increases (red shift) if the trench depth h increases or its angle &#x3b8; decreases. In particular, the wavelength of the resonance increases (red shift) if the trench depth h increases or its angle <inline-formula id="inf41">
<mml:math id="m49">
<mml:mi>&#x3b8;</mml:mi>
</mml:math>
</inline-formula> decreases. On the other hand, the amplification factor of the electric field tends to decrease when the resonances are at long wavelengths. This behaviour can be explained by considering that gold at long wavelengths tends to behave like an ideal metal, absorbing most of the SPPs that form along the edges. The plasmonic structure generates an electromagnetic field with an asymmetric distribution, with more amplitude at the bottom of the V-shape and less at the top of the trench, at the surface of the metal layer. This distribution can be explained by taking into account that SPPs tend to reduce their effectiveness as they propagate along the edge of the trench. Starting from the spatial distribution of the electromagnetic field and integrating the Maxwell tensor on the surface of a nanoparticle (<xref ref-type="disp-formula" rid="e2">Eq. 2</xref>), the forces acting on the particle and its outward ejection can be determined. The force obtained from the Maxwell tensor can be written as in equation (<xref ref-type="disp-formula" rid="e3">Eq. 3</xref>) with a dominant gradient term with respect to the radiation pressure. In particular, the gradient force is such that the nanoparticle will move along the <italic>z</italic>-axis (i.e., along the height of the V-shape) with a velocity derived from formula (<xref ref-type="disp-formula" rid="e2">Eq. 2</xref>), that takes on the following value (<xref ref-type="bibr" rid="B25">Shalin and Sukhov, 2013</xref>):<disp-formula id="e10">
<mml:math id="m50">
<mml:mrow>
<mml:mrow>
<mml:mover accent="true">
<mml:mi mathvariant="italic">v</mml:mi>
<mml:mo>&#x2192;</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:msqrt>
<mml:mrow>
<mml:mfrac>
<mml:mn>3</mml:mn>
<mml:mrow>
<mml:mn>2</mml:mn>
<mml:msub>
<mml:mi mathvariant="italic">&#x3c1;</mml:mi>
<mml:mrow>
<mml:mi mathvariant="italic">eff</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mfrac>
<mml:mrow>
<mml:mover accent="true">
<mml:mi mathvariant="italic">&#x3b1;</mml:mi>
<mml:mo>&#x2dc;</mml:mo>
</mml:mover>
</mml:mrow>
<mml:msub>
<mml:mi mathvariant="italic">&#x3f5;</mml:mi>
<mml:mn>0</mml:mn>
</mml:msub>
</mml:mrow>
</mml:msqrt>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mrow>
<mml:mover accent="true">
<mml:mi mathvariant="italic">E</mml:mi>
<mml:mo>&#x2192;</mml:mo>
</mml:mover>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(10)</label>
</disp-formula>where <inline-formula id="inf42">
<mml:math id="m51">
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="italic">&#x3c1;</mml:mi>
<mml:mrow>
<mml:mi mathvariant="italic">eff</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the effective density of the nanoparticle whereas <inline-formula id="inf43">
<mml:math id="m52">
<mml:mrow>
<mml:mover accent="true">
<mml:mi mathvariant="italic">&#x3b1;</mml:mi>
<mml:mo>&#x2dc;</mml:mo>
</mml:mover>
</mml:mrow>
</mml:math>
</inline-formula> represents the polarizability per unit volume that is a function of the dielectric permittivity relative to the substrate and metal. The equation (<xref ref-type="disp-formula" rid="e10">Eq. 10</xref>) is derived by considering that the force is entirely a gradient force and that the potential generated by it is entirely transformed into kinetic energy.</p>
<p>Numerical simulations (<xref ref-type="bibr" rid="B25">Shalin and Sukhov, 2013</xref>) were carried out on the structure and the velocity of a silver nanoparticle with a radius of 8&#xa0;nm was estimated. The estimated velocity is 0.1&#xa0;m/s when illuminating the plasmonic structure by a TM-polarized wave having a power density of 1&#x2013;2 <inline-formula id="inf44">
<mml:math id="m53">
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:mi mathvariant="bold-italic">mW</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi mathvariant="bold-italic">&#x3bc;</mml:mi>
<mml:msup>
<mml:mi mathvariant="bold-italic">m</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:math>
</inline-formula> and a wavelength in the visible range (required for the SPPs formation). The speeds obtained are such as to allow the transport of drugs and genetic material in the cells and are therefore interesting for possible applications in the medical-biological sector. However, a single structure is not capable of inducing high velocities in nanoparticles, so the propulsion described above can be used for applications requiring the movement of small objects in extremely small areas.</p>
<p>Rovey&#x2019;s group has designed a plasmonic structure that allows higher values of thrust and velocity: (<xref ref-type="bibr" rid="B20">Rovey et al., 2015</xref>; <xref ref-type="bibr" rid="B8">Jiao et al., 2016</xref>; <xref ref-type="bibr" rid="B15">Maser et al., 2016</xref>; <xref ref-type="bibr" rid="B35">Zhang et al., 2019</xref>). This is not due to an increase in the force generated by the single structure, but to the possibility of being able to couple a large number of elementary structures in such a way as to form an array and thus obtain a force resulting from the sum of the elementary contributions related to the single structures. Thus, the total force is greatly amplified and can be used in numerous applications, including space propulsion. The system (see <xref ref-type="fig" rid="F9">Figure 9</xref>) consists of two trapezoidal gold structures separated by a gap and deposited on a glass substrate.</p>
<fig id="F9" position="float">
<label>FIGURE 9</label>
<caption>
<p>Sketch of the asymmetrical trapezoidal plasmonic structure. Adapted from <xref ref-type="bibr" rid="B8">Jiao et al. (2016)</xref>
</p>
</caption>
<graphic xlink:href="fnano-04-886636-g009.tif"/>
</fig>
<p>When the structures are uniformly illuminated, a resonance, due to the generation of an LSP, is formed in the gap between the two trapezoids. The resonance wavelength depends on the geometrical parameters of the gold structures, i.e. the width of the two bases of the trapezium, its length, its thickness and the gap. Once the resonance has been induced, the near-field excitation forms two SPPs along the edges of the two trapezoids that couple within the gap and, as in the structure discussed in the previous paragraph, generate a strongly amplified and asymmetric field distribution. By calculating the forces through the Maxwell tensor (<xref ref-type="disp-formula" rid="e1">Eq. 1</xref>), we obtain a distribution that depends on the position along y of the system (see <xref ref-type="fig" rid="F10">Figure 10A</xref>), such that a nanoparticle undergoes a strong initial acceleration followed by a deceleration zone of lesser intensity (see <xref ref-type="fig" rid="F10">Figure 10B</xref>).</p>
<fig id="F10" position="float">
<label>FIGURE 10</label>
<caption>
<p>
<bold>(A)</bold> A cut in the xy plane of the trapezoidal structure, showing the direction of the ejection force; <bold>(B)</bold> the graph of the force as the y coordinate of the nanoparticle changes. Adapted from <xref ref-type="bibr" rid="B20">Rovey et al. (2015)</xref>.</p>
</caption>
<graphic xlink:href="fnano-04-886636-g010.tif"/>
</fig>
<p>Preliminary experimental results for this configuration are also reported in refs. (<xref ref-type="bibr" rid="B20">Rovey et al., 2015</xref>; <xref ref-type="bibr" rid="B8">Jiao et al., 2016</xref>; <xref ref-type="bibr" rid="B15">Maser et al., 2016</xref>; <xref ref-type="bibr" rid="B35">Zhang et al., 2019</xref>), confirming the functionality expected from theory and simulations. In general, it is observed that structures with larger horizontal asymmetries (i.e., along the <italic>y</italic>-axis) have more amplified electric fields. Therefore, in order to induce further asymmetry, the oblique side of the trapezoid can be modified in a stepwise manner.</p>
<p>
<xref ref-type="fig" rid="F11">Figure 11</xref> shows a comparison of the field distribution between a trapezoidal structure with linear and stepped oblique sides (<xref ref-type="bibr" rid="B35">Zhang et al., 2019</xref>).</p>
<fig id="F11" position="float">
<label>FIGURE 11</label>
<caption>
<p>
<bold>(A)</bold> Electric field distribution in a trapezoidal structure the wave is polarised along x at wavelength 532&#xa0;nm; <bold>(B)</bold> electric field distribution of a modified oblique-sided trapezoidal structure with the wave polarised along x and at wavelength 580&#xa0;nm. &#x00A9; 2019 IEEE. Reprinted, with permission, from <xref ref-type="bibr" rid="B35">Zhang et al. (2019)</xref>.</p>
</caption>
<graphic xlink:href="fnano-04-886636-g011.tif"/>
</fig>
<p>As mentioned above, in general, the output speed induced by a single plasmonic structure is relatively low. In order to obtain more efficient performance useful for practical applications, an array of the above elementary structures can be constructed by organising an array through a so-called end-to-end configuration (<xref ref-type="fig" rid="F12">Figure 12A</xref>). In this configuration, the minor base of a trapezoid is adjacent to the major base of the next structure. Using this spatial configuration, it is possible to obtain a continuous sequence of pushes on the nanoparticle, which implies a considerable increase in the force applied and, consequently, in the output speed.</p>
<fig id="F12" position="float">
<label>FIGURE 12</label>
<caption>
<p>
<bold>(A)</bold> A sketch of the array configuration in end-to-end mode; <bold>(B)</bold> the thruster hypothesised in <xref ref-type="bibr" rid="B20">Rovey et al. (2015)</xref> to make full use of sunlight. Adapted from <xref ref-type="bibr" rid="B20">Rovey et al. (2015)</xref>.</p>
</caption>
<graphic xlink:href="fnano-04-886636-g012.tif"/>
</fig>
<p>The plasmonic structure was arranged in a multilayer array configuration, where each layer has a different resonance wavelength, to function properly even with broadband illumination such as solar radiation (<xref ref-type="bibr" rid="B20">Rovey et al., 2015</xref>) (see <xref ref-type="fig" rid="F12">Figure 12B</xref>). In fact, with this approach, each layer of the stack is able to utilise a specific part of the solar spectrum, increasing the overall efficiency. The nanoparticles are contained in nano-channels between the active layers, so that the resulting force applied to the nanoparticles will be a combination of the forces obtained from the arrays above and below each channel.</p>
<p>Thanks to its estimated performance, this plasmonic configuration has been proposed as a thruster for space applications. The values of thrust, specific impulse and ejection velocity of a nanoparticle for the thruster based on the multilayer array configuration are estimated using the following relationships:<disp-formula id="equ1">
<mml:math id="m54">
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:msubsup>
<mml:mi mathvariant="italic">v</mml:mi>
<mml:mi mathvariant="italic">i</mml:mi>
<mml:mn>2</mml:mn>
</mml:msubsup>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:mfrac>
<mml:mo>&#x3d;</mml:mo>
<mml:munderover>
<mml:mstyle displaystyle="true">
<mml:mo>&#x222b;</mml:mo>
</mml:mstyle>
<mml:mn>0</mml:mn>
<mml:mi mathvariant="italic">L</mml:mi>
</mml:munderover>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mi mathvariant="italic">F</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi mathvariant="italic">above</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi mathvariant="italic">y</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mrow>
<mml:mi mathvariant="italic">F</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi mathvariant="italic">below</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi mathvariant="italic">y</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mi mathvariant="italic">m</mml:mi>
</mml:mfrac>
<mml:mi mathvariant="italic">dy</mml:mi>
</mml:mrow>
</mml:math>
</disp-formula>
<disp-formula id="e11">
<mml:math id="m55">
<mml:mrow>
<mml:mtable>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:mrow>
<mml:mover accent="true">
<mml:mi mathvariant="italic">T</mml:mi>
<mml:mo>&#x2192;</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mo>&#x2211;</mml:mo>
<mml:mi mathvariant="italic">mf</mml:mi>
<mml:msub>
<mml:mrow>
<mml:mover accent="true">
<mml:mi mathvariant="italic">v</mml:mi>
<mml:mo>&#x2192;</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mi mathvariant="italic">i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mtd>
<mml:mtd>
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="italic">I</mml:mi>
<mml:mrow>
<mml:mi mathvariant="italic">sp</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mo>&#x007C;</mml:mo>
<mml:mover accent="true">
<mml:mi mathvariant="italic">T</mml:mi>
<mml:mo>&#x2192;</mml:mo>
</mml:mover>
<mml:mo>&#x007C;</mml:mo>
</mml:mrow>
<mml:mrow>
<mml:mi mathvariant="italic">Nmf</mml:mi>
<mml:msub>
<mml:mi mathvariant="italic">g</mml:mi>
<mml:mn>0</mml:mn>
</mml:msub>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:math>
<label>(11)</label>
</disp-formula>where <italic>m</italic> is the mass of the nanoparticle, <italic>N</italic> is the number of arrays in a single layer (the number of rows in <xref ref-type="fig" rid="F12">Figure 12B</xref>), <italic>f</italic> is the ejection rate of the nanoparticle, <italic>L</italic> is the acceleration length, and <inline-formula id="inf45">
<mml:math id="m56">
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mi mathvariant="italic">F</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi mathvariant="italic">above</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mrow>
<mml:mi mathvariant="italic">F</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi mathvariant="italic">below</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> are the force distributions generated by the arrays located above and below the nanoparticle under consideration, respectively. The results obtained from the simulations show propulsion values of around 400&#x2013;600&#xa0;nN and specific impulse around 20 s, with an exit velocity of 3&#x2013;4&#xa0;m/s (<xref ref-type="bibr" rid="B20">Rovey et al., 2015</xref>). Obviously, these parameters are highly dependent on the area of the plasmonic lens, which determines the actual value of the optical power density illuminating the structure.</p>
<p>However, as can be seen in <xref ref-type="fig" rid="F13">Figure 13</xref>, the array presents a potential limitation (<xref ref-type="bibr" rid="B21">Rovey et al., 2019</xref>), namely when the elementary structures are organised in array with an end-to-end configuration, the distribution of the field supported by the plasmonic structure changes compared to the case of the isolated elementary structure, and with it, obviously, the spatial distribution of the force acting on the nanoparticles also changes. In fact, in the array configuration, the distribution of forces is such that the particle must overcome a potential barrier to pass from one elementary structure to the next. The intensity of this potential barrier is such that the kinetic energy of the particle is not sufficient to allow its overcoming in the short space in which the particle acquires velocity. For this reason, the particle needs an additional initial thrust that is sufficiently high to allow it to overcome the potential barrier and continue its motion.</p>
<fig id="F13" position="float">
<label>FIGURE 13</label>
<caption>
<p>Y-component of the force applied to a glass particle of radius 50&#xa0;nm as a function of the y-coordinate of a set of trapezoidal structures as in <xref ref-type="bibr" rid="B20">Rovey et al. (2015)</xref>, placed in an end-to-end configuration with the following parameters: height 350&#xa0;nm, major base 120&#xa0;nm, minor base 60&#xa0;nm, gap 40&#xa0;nm, thickness 50&#xa0;nm. The light is polarized along x and at wavelength 730&#xa0;nm and the amplitude of the incident electric field is such that it is normalized.</p>
</caption>
<graphic xlink:href="fnano-04-886636-g013.tif"/>
</fig>
</sec>
<sec sec-type="conclusion" id="s5">
<title>Conclusion</title>
<p>The growing demand for the manipulation of microscopic objects using long-range techniques has increasingly required the development of techniques capable of spatially intensifying and concentrating electromagnetic fields with the aim of improving the electromagnetic forces acting on objects. In this context, one of the most interesting techniques is based on the use of the surface plasmonics phenomena that have the capability to amplify the electric field in areas much smaller than the diffractive limit. Research in recent years has been trying to exploit this ability of surface plasmonics to realize structures that can control the motion and trajectory of a nanoparticle. Properly designed plasmonic structures have been conceived and characterised to generate hotspots, i.e., symmetrical field distributions, such that the resulting forces are able to trap particles in a specific spatial region. Thanks to the potential of plasmonic structures to engineer the spatial distribution of generated electromagnetic fields, plasmonic phenomena are recently being studied to induce also studied to induce long-range movement of nanoparticles. Some interesting applications have already been imagined for these structures, such as the transport of molecules and genetic material within cells or the propulsion of micro/nanosatellites. The latter appears to be particularly fascinating and in fact there are already a good number of publications on the optimisation of usable configurations. However, there are still limitations to the practical and efficient use of these structures; in particular, the presence of a potential barrier which inhibits the achievement of a large exit velocity. To address this limitation, a new design of the structure or, as already demonstrated for dynamic optical tweezers, a suitable modification of the incident beam could be useful. It is therefore clear that resonant plasmonic structures, due to their properties and design versatility, are an interesting field of research for the induction of microscopic particle motion by light radiation and in particular for the development of future optical propulsion techniques useful for the next generation of micro/nanosatellite constellations.</p>
</sec>
</body>
<back>
<sec id="s6">
<title>Author Contributions</title>
<p>GC and MI conceived the idea of the paper. SB and GZ formalized the theoretical section. SB detailed the optical tweezer section. GC, MI and SB described the plasmonic propulsion section. All the authors analyzed the data and made additional effort in preparing the paper.</p>
</sec>
<sec id="s7">
<title>Funding</title>
<p>&#x201c;PON ARS01_01181 PM3&#x2014;Multi Mission Modular Platform&#x201d; Industrial Research and Experimental Development Projects in the 12 areas of specialisation identified by the National Research Programme 2015-2020&#x2014;CUP B66G18000740005.</p>
</sec>
<sec sec-type="COI-statement" id="s8">
<title>Conflict of Interest</title>
<p>The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.</p>
</sec>
<sec sec-type="disclaimer" id="s9">
<title>Publisher&#x2019;s Note</title>
<p>All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.</p>
</sec>
<ref-list>
<title>References</title>
<ref id="B1">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Cervone</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Zandbergen</surname>
<given-names>B.</given-names>
</name>
<name>
<surname>Guerrieri</surname>
<given-names>D. C.</given-names>
</name>
<name>
<surname>De Athayde Costa e Silva</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Krusharev</surname>
<given-names>I.</given-names>
</name>
<name>
<surname>van Zeijl</surname>
<given-names>H.</given-names>
</name>
</person-group> (<year>2016</year>). <article-title>Green Micro-Resistojet Research at Delft University of Technology: The New Frontiers</article-title>. <source>CEAS Space J.</source> <volume>9</volume> (<issue>1</issue>), <fpage>111</fpage>&#x2013;<lpage>125</lpage>. <pub-id pub-id-type="doi">10.1007/s12567-016-0135-3</pub-id> </citation>
</ref>
<ref id="B2">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Chen</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Ng</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Lin</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Chan</surname>
<given-names>C.</given-names>
</name>
</person-group> (<year>2011</year>). <article-title>Optical Pulling Force</article-title>. <source>Nat. Photon.</source> <volume>5</volume>, <fpage>531</fpage>&#x2013;<lpage>534</lpage>. <pub-id pub-id-type="doi">10.1038/nphoton.2011.153</pub-id> </citation>
</ref>
<ref id="B3">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Daly</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Sergides</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Chormaic</surname>
<given-names>S. N.</given-names>
</name>
</person-group> (<year>2015</year>). <article-title>Optical Trapping and Manipulation of Micrometer and Submicrometer Particles</article-title>. <source>Laser Photon. Rev</source> <volume>9</volume> (<issue>3</issue>), <fpage>309</fpage>&#x2013;<lpage>329</lpage>. <pub-id pub-id-type="doi">10.1002/lpor.201500006</pub-id> </citation>
</ref>
<ref id="B5">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Fujimaki</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Rockstuhl</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Awazu</surname>
<given-names>K.</given-names>
</name>
<name>
<surname>Tominaga</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Fukuda</surname>
<given-names>N.</given-names>
</name>
<etal/>
</person-group> (<year>2008</year>). <article-title>The Design of Evanescent-Field-Coupled Waveguide-Mode Sensors</article-title>. <source>Nanotechnology</source> <volume>19</volume>, <fpage>095503</fpage>. <pub-id pub-id-type="doi">10.1088/0957-4484/19/9/095503</pub-id> </citation>
</ref>
<ref id="B6">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Gassend</surname>
<given-names>B.</given-names>
</name>
<name>
<surname>Vel&#xe1;squez-Garc&#xed;a</surname>
<given-names>L. F.</given-names>
</name>
<name>
<surname>Akinwande</surname>
<given-names>A. I.</given-names>
</name>
<name>
<surname>Mart&#xed;nez-S&#xe1;nchez</surname>
<given-names>M.</given-names>
</name>
</person-group> (<year>2009</year>). <article-title>A Microfabricated Planar Electrospray Array Ionic Liquid Ion Source with Integrated Extractor</article-title>. <source>J. Microelectromechanical Syst.</source> <volume>18</volume> (<issue>3</issue>), <fpage>679</fpage>&#x2013;<lpage>694</lpage>. <comment>Available at: <ext-link ext-link-type="uri" xlink:href="http://hdl.handle.net/1721.1/59347">http://hdl.handle.net/1721.1/59347</ext-link>
</comment>. <pub-id pub-id-type="doi">10.1109/jmems.2009.2015475</pub-id> </citation>
</ref>
<ref id="B7">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Huft</surname>
<given-names>P. R.</given-names>
</name>
<name>
<surname>Kolbow</surname>
<given-names>J. D.</given-names>
</name>
<name>
<surname>Thweatt</surname>
<given-names>J. T.</given-names>
</name>
<name>
<surname>Lindquist</surname>
<given-names>N. C.</given-names>
</name>
</person-group> (<year>2017</year>). <article-title>Holographic Plasmonic Nanotweezers for Dynamic Trapping and Manipulation</article-title>. <source>Nano Lett.</source> <volume>17</volume>, <fpage>7920</fpage>&#x2013;<lpage>7925</lpage>. <pub-id pub-id-type="doi">10.1021/acs.nanolett.7b04289</pub-id> </citation>
</ref>
<ref id="B8">
<citation citation-type="confproc">
<person-group person-group-type="author">
<name>
<surname>Jiao</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Lin</surname>
<given-names>E.</given-names>
</name>
<name>
<surname>Gaofeng</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Zhao</surname>
<given-names>Q.</given-names>
</name>
</person-group> (<year>2016</year>). &#x201c;<article-title>Nano Optical Propeller Based on Localized Field Intensity Enhancement of Surface Plasmons</article-title>,&#x201d; in <conf-name>SPIE</conf-name>. <comment>Chengdu, China, Yongkun Ding</comment>. <pub-id pub-id-type="doi">10.1117/12.2268294</pub-id> </citation>
</ref>
<ref id="B9">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Juan</surname>
<given-names>M. L.</given-names>
</name>
<name>
<surname>Rrighini</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Quidant</surname>
<given-names>R.</given-names>
</name>
</person-group> (<year>2011</year>). <article-title>Plasmon Nano-Optical Tweezers</article-title>. <source>Nat. Photon.</source> <volume>5</volume>, <fpage>349</fpage>&#x2013;<lpage>356</lpage>. <pub-id pub-id-type="doi">10.1038/nphoton.2011.56</pub-id> </citation>
</ref>
<ref id="B10">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Lemmer</surname>
<given-names>K.</given-names>
</name>
</person-group> (<year>2017</year>). <article-title>Propulsion for CubeSat</article-title>. <source>Acta Astronautica</source> <volume>134</volume>, <fpage>231</fpage>&#x2013;<lpage>243</lpage>. <pub-id pub-id-type="doi">10.1016/j.actaastro.2017.01.048</pub-id> </citation>
</ref>
<ref id="B11">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Li</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Lin</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Inoue</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Zheng</surname>
<given-names>Y.</given-names>
</name>
</person-group> (<year>2018</year>). <article-title>Opto-Thermophoretic Tweezers and Assembly</article-title>. <source>J. Micro- Nano-Manufactoring</source> <volume>6</volume>, <fpage>040801</fpage>. <pub-id pub-id-type="doi">10.1115/1.4041615</pub-id> </citation>
</ref>
<ref id="B12">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Liaw</surname>
<given-names>J.-W.</given-names>
</name>
<name>
<surname>Chien</surname>
<given-names>C-W.</given-names>
</name>
<name>
<surname>Liu</surname>
<given-names>K-C.</given-names>
</name>
<name>
<surname>Ku</surname>
<given-names>Y-C.</given-names>
</name>
<name>
<surname>Kuo</surname>
<given-names>M-K.</given-names>
</name>
</person-group> (<year>2018</year>). <article-title>3D Optical Vortex Trapping of Plasmonic Nanostructure</article-title>. <source>Scientific Rep.</source> <volume>8</volume>, <fpage>12673</fpage>. <pub-id pub-id-type="doi">10.1038/s41598-018-30948-y</pub-id> </citation>
</ref>
<ref id="B13">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Maier</surname>
<given-names>S. A.</given-names>
</name>
</person-group> (<year>2007</year>). <source>Plasmonics: Fundamental and Applications</source>. <edition>I ed.</edition> <publisher-loc>New York, NY</publisher-loc>: <publisher-name>Springer</publisher-name>. </citation>
</ref>
<ref id="B14">
<citation citation-type="confproc">
<person-group person-group-type="author">
<name>
<surname>Mansuripur</surname>
<given-names>M.</given-names>
</name>
</person-group> (<year>2014</year>). &#x201c;<article-title>Mechanical Effects of Light on Material Media: Radiation Pressure and the Linear and Angular Momenta of Photons</article-title>,&#x201d; in <conf-name>SPIE at The University of Arizona, Tucson, US, Harrison H. Barret</conf-name>. <pub-id pub-id-type="doi">10.1117/12.2063432</pub-id> </citation>
</ref>
<ref id="B15">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Maser</surname>
<given-names>J. N.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Deng</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Yang</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Rovey</surname>
<given-names>J. L.</given-names>
</name>
</person-group> (<year>2016</year>). <article-title>Transmission Spectrum of Asymmetric Nanostructures for Plasmonic Space Propulsion</article-title>. <source>J. Spacecraft Rockets</source> <volume>53</volume> (<issue>5</issue>), <fpage>998</fpage>&#x2013;<lpage>1000</lpage>. <pub-id pub-id-type="doi">10.2514/1.A33576</pub-id> </citation>
</ref>
<ref id="B16">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ndukaife</surname>
<given-names>J. C.</given-names>
</name>
<name>
<surname>Kildishev</surname>
<given-names>A. V.</given-names>
</name>
<name>
<surname>Nnanna</surname>
<given-names>A. G. A.</given-names>
</name>
<name>
<surname>Shalaev</surname>
<given-names>V. M.</given-names>
</name>
<name>
<surname>Wereley</surname>
<given-names>S. T.</given-names>
</name>
<name>
<surname>Boltasseva</surname>
<given-names>A.</given-names>
</name>
</person-group> (<year>2015</year>). <article-title>Long-Range and Rapid Transport of Individual Nano-Objects by a Hybrid Electrothermoplasmonic Nanotweezer</article-title>. <source>Nat. Nanotechnologies</source> <volume>11</volume> (<issue>1</issue>), <fpage>53</fpage>&#x2013;<lpage>59</lpage>. <pub-id pub-id-type="doi">10.1038/NNANO.2015.248</pub-id> </citation>
</ref>
<ref id="B17">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Novitsky</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Qiu</surname>
<given-names>C.-W.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>H.</given-names>
</name>
</person-group> (<year>2011</year>). <article-title>Single Gradientless Light Beam Drags Particles as Tractor Beams</article-title>. <source>Phys. Rev. Lett.</source> <volume>107</volume> (<issue>20</issue>), <fpage>203601</fpage>. <pub-id pub-id-type="doi">10.1103/PhysRevLett.107.203601</pub-id> </citation>
</ref>
<ref id="B18">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Novothny</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Hecht</surname>
<given-names>B.</given-names>
</name>
</person-group> (<year>2006</year>). <source>Principles of Nano-Optics</source>. <edition>I ed.</edition> <publisher-loc>New York</publisher-loc>: <publisher-name>Cambridge University Press</publisher-name>. </citation>
</ref>
<ref id="B19">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Roichman</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Sun</surname>
<given-names>B.</given-names>
</name>
<name>
<surname>Roichman</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Amato-Grill</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Grier</surname>
<given-names>D. G.</given-names>
</name>
</person-group> (<year>2008</year>). <article-title>Optical Forces Arising from Phase Gradients</article-title>. <source>Phys. Rev. Lett.</source> <volume>100</volume>, <fpage>013602</fpage>. <pub-id pub-id-type="doi">10.1103/PhysRevLett.100.013602</pub-id> </citation>
</ref>
<ref id="B20">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Rovey</surname>
<given-names>J. L.</given-names>
</name>
<name>
<surname>Friz</surname>
<given-names>P. D.</given-names>
</name>
<name>
<surname>Hu</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Glascock</surname>
<given-names>M. S.</given-names>
</name>
<name>
<surname>Yang</surname>
<given-names>X.</given-names>
</name>
</person-group> (<year>2015</year>). <article-title>Plasmonic Force Space Propulsion</article-title>. <source>J. Spacecraft Rockets</source> <volume>52</volume> (<issue>4</issue>), <fpage>1163</fpage>&#x2013;<lpage>1168</lpage>. <pub-id pub-id-type="doi">10.2514/1.A33155</pub-id> </citation>
</ref>
<ref id="B21">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Rovey</surname>
<given-names>J. L.</given-names>
</name>
<name>
<surname>Joshua</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Jie</surname>
<given-names>G.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>Experimental Demonstration and System Analysis for Plasmonic Force Propulsion</article-title>. <comment>No. HQ-E-DAA-TN75812.: Nasa Report</comment>. </citation>
</ref>
<ref id="B22">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ruffner</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Grier</surname>
<given-names>D.</given-names>
</name>
</person-group> (<year>2012</year>). <article-title>Optical Forces and Torques in Non Uniform Beams of Light</article-title>. <source>Phys. Rev. Lett.</source> <volume>108</volume>, <fpage>173602</fpage>. <pub-id pub-id-type="doi">10.1103/PhysRevLett.108.173602</pub-id> </citation>
</ref>
<ref id="B23">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Sasoh</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Tsuruta</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Katagiri</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Dondelewski</surname>
<given-names>O.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>B.</given-names>
</name>
</person-group> (<year>2017</year>). <source>Characteristics of Ablation Impulse Induced by Repetitive Laser Pulse Irradiations</source>. <publisher-loc>Darmstadt, Germany</publisher-loc>: <publisher-name>T. Flohrer; F. Schmitz</publisher-name>. </citation>
</ref>
<ref id="B24">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Sathiyanathan</surname>
<given-names>K.</given-names>
</name>
<name>
<surname>Lee</surname>
<given-names>R.</given-names>
</name>
<name>
<surname>Chesser</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Dubois</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Stowe</surname>
<given-names>R.</given-names>
</name>
<name>
<surname>Farinaccio</surname>
<given-names>R.</given-names>
</name>
<etal/>
</person-group> (<year>2011</year>). <article-title>Solid Propellant Microthruster Design for Nanosatellite Applications</article-title>. <source>J. Propulsion Power</source> <volume>27</volume> (<issue>6</issue>), <fpage>1288</fpage>&#x2013;<lpage>1294</lpage>. <pub-id pub-id-type="doi">10.251/1.B34109</pub-id> </citation>
</ref>
<ref id="B25">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Shalin</surname>
<given-names>A. S.</given-names>
</name>
<name>
<surname>Sukhov</surname>
<given-names>S. V.</given-names>
</name>
</person-group> (<year>2013</year>). <article-title>Plasmonic Nanostructures as Accelerators for Nanoparticles: Optical Nanocannon</article-title>. <source>Plasmonics</source> <volume>8</volume>, <fpage>625</fpage>&#x2013;<lpage>629</lpage>. <pub-id pub-id-type="doi">10.1007/s11468-012-9447-0</pub-id> </citation>
</ref>
<ref id="B26">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Shuller</surname>
<given-names>J. A.</given-names>
</name>
<name>
<surname>Barnard</surname>
<given-names>E. S.</given-names>
</name>
<name>
<surname>Cai</surname>
<given-names>W.</given-names>
</name>
<name>
<surname>Jun</surname>
<given-names>Y. C.</given-names>
</name>
<name>
<surname>White</surname>
<given-names>J. S.</given-names>
</name>
<name>
<surname>Brongersma</surname>
<given-names>M. L.</given-names>
</name>
</person-group> (<year>2010</year>). <article-title>Plasmonic for Extreme Light Concentration and Manipulation</article-title>. <source>Natur. Mater.</source> <volume>9</volume>, <fpage>193</fpage>&#x2013;<lpage>204</lpage>. <pub-id pub-id-type="doi">10.1038/nmat2630</pub-id> </citation>
</ref>
<ref id="B27">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>S&#xf8;ndergaard</surname>
<given-names>T.</given-names>
</name>
<name>
<surname>Bozhevolnyi</surname>
<given-names>S. I.</given-names>
</name>
<name>
<surname>Beermann</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Novikov</surname>
<given-names>S. M.</given-names>
</name>
</person-group> (<year>2010</year>). <article-title>Resonant Plasmon Nanofocusing by Closed Tapered Gaps</article-title>. <source>Nano Lett.</source> <volume>10</volume>, <fpage>291</fpage>&#x2013;<lpage>295</lpage>. <pub-id pub-id-type="doi">10.1021/nl903563e</pub-id> </citation>
</ref>
<ref id="B28">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Soukoulis</surname>
<given-names>C. M.</given-names>
</name>
</person-group> (<year>2002</year>). <article-title>The History and a Review of the Modeling and Fabrication of Photonic Crystal</article-title>. <source>Nanotechnology</source> <volume>13</volume> (<issue>3</issue>), <fpage>420</fpage>. <pub-id pub-id-type="doi">10.1088/0957-4484/13/3/335</pub-id> </citation>
</ref>
<ref id="B29">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Swartzlander</surname>
<given-names>J. G. A.</given-names>
</name>
</person-group> (<year>2017</year>). <article-title>Radiation Pressure on a Diffractive Sailcraft</article-title>. <source>J. Opt. Soc. America B</source> <volume>34</volume> (<issue>6</issue>), <fpage>C25</fpage>&#x2013;<lpage>C30</lpage>. <pub-id pub-id-type="doi">10.1364/josab.34.000c25</pub-id> </citation>
</ref>
<ref id="B30">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Tummala</surname>
<given-names>A. R.</given-names>
</name>
<name>
<surname>Dutta</surname>
<given-names>A.</given-names>
</name>
</person-group> (<year>2017</year>). <article-title>An Overview of Cube-Satellite Propulsion Technologies and Trends</article-title>. <source>Aerospace</source> <volume>4</volume>, <fpage>58</fpage>. <pub-id pub-id-type="doi">10.3390/aerospace4040058</pub-id> </citation>
</ref>
<ref id="B36">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>W. Davis</surname>
<given-names>E.</given-names>
</name>
<name>
<surname>Franklin</surname>
<given-names>B.</given-names>
</name>
<name>
<surname>Mead</surname>
<given-names>J.</given-names>
</name>
</person-group> (<year>2011</year>). <article-title>Review of Laser Lightcraft Propulsion System</article-title>. <comment>American Istitute of Physics, Mariland, US, A.V. Pakhomov.</comment> </citation>
</ref>
<ref id="B31">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wang</surname>
<given-names>K.</given-names>
</name>
<name>
<surname>Schonbrun</surname>
<given-names>E.</given-names>
</name>
<name>
<surname>Crozier</surname>
<given-names>K. B.</given-names>
</name>
</person-group> (<year>2009</year>). <article-title>Propulsion of Gold Nanoparticles with Surface Plasmon Polaritons: Evidence of Enhanced Optical Force from Near-Field Coupling Between Gold Particle and Gold Film</article-title>. <source>Nano Lett.</source> <volume>9</volume> (<issue>7</issue>), <fpage>2623</fpage>&#x2013;<lpage>2629</lpage>. <pub-id pub-id-type="doi">10.1021/nl900944y</pub-id> </citation>
</ref>
<ref id="B32">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wang</surname>
<given-names>N.</given-names>
</name>
<name>
<surname>Chen</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Liu</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Lin</surname>
<given-names>Z.</given-names>
</name>
</person-group> (<year>2013</year>). <article-title>Dynamical and Phase-Diagram Study on Stable Optical Pulling Force in Bessel Beams</article-title>. <source>Phys. Rev. A.</source> <volume>87</volume>, <fpage>063812</fpage>. <pub-id pub-id-type="doi">10.1103/physreva.87.063812</pub-id> </citation>
</ref>
<ref id="B33">
<citation citation-type="confproc">
<person-group person-group-type="author">
<name>
<surname>Wu</surname>
<given-names>M. C.</given-names>
</name>
</person-group> (<year>2007</year>). &#x201c;<article-title>Optoelectronic Tweezers for Nanomanipulation</article-title>,&#x201d; in <conf-name>2007 International Nano-Optoelectronics Workshop</conf-name>, <conf-loc>Bejing, China</conf-loc>, <conf-date>29 July-11 Aug. 2007</conf-date> (<publisher-name>IEEE</publisher-name>). <pub-id pub-id-type="doi">10.1109/INOW.2007.4302858</pub-id> </citation>
</ref>
<ref id="B34">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Yang</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Liu</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Oulton</surname>
<given-names>R. F.</given-names>
</name>
<name>
<surname>Yin</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>X.</given-names>
</name>
</person-group> (<year>2011</year>). <article-title>Optical Forces in Hybrid Plasmonics Waveguides</article-title>. <source>Nano Lett.</source> <volume>11</volume>, <fpage>321</fpage>&#x2013;<lpage>328</lpage>. <pub-id pub-id-type="doi">10.1021/nl103070n</pub-id> </citation>
</ref>
<ref id="B35">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhang</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Jiao</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Tong</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Ma</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Zhao</surname>
<given-names>Q.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>Study on Optical Field Distribution of Asymmetric Nanostructures Based on Surface Plasmon</article-title>. <source>IEEE Photon. J.</source> <volume>11</volume> (<issue>2</issue>), <fpage>5700108</fpage>. <pub-id pub-id-type="doi">10.1109/JPHOT.2019.2904335</pub-id> </citation>
</ref>
</ref-list>
</back>
</article>