Design of plasmonic directional antennas via evolutionary optimization
Peter R. Wiecha Clément Majorel Christian Girard Aurélien Cuche Vincent Paillard Otto L. Muskens Arnaud Arbouet













































^1Physics and Astronomy, Faculty of Engineering and Physical Sciences, University of Southampton, Southampton, UK

^2CEMES-CNRS, Université de Toulouse, CNRS, UPS, Toulouse, France


^*p.wiecha@soton.ac.uk











We demonstrate inverse design of plasmonic nanoantennas for directional light scattering. 
Our method is based on a combination of full-field electrodynamical simulations via the Green dyadic method and evolutionary optimization (EO).
Without any initial bias, we find that the geometries reproducibly found by EO, work on the same principles as radio-frequency antennas. 
We demonstrate the versatility of our approach by designing various directional optical antennas for different scattering problems.
EO based nanoantenna design has tremendous potential for a multitude of applications like nano-scale information routing and processing or single-molecule spectroscopy.
Furthermore, EO can help to derive general design rules and to identify inherent physical limitations for photonic nanoparticles and metasurfaces.










The vast opportunities of light-based nanotechnology lead to tremendous research efforts on light-matter interaction at sub-wavelength scales, and in particular gave rise to the broad field of plasmonics, as localized surface plasmon-polariton (LSP) resonances in metallic nanostructures can strongly confine far-field radiation <cit.>.
By variations of the particle geometry or its material it is possible to tailor manifold optical properties like resonance positions <cit.>, polarization conversion <cit.>, optical chirality <cit.>, localized heat generation <cit.> or nonlinear optical effects <cit.> with many applications like field enhanced spectroscopy <cit.> or refractive index sensing <cit.>.

Another optical functionality that is in the focus of the plasmonics community is the directional routing of light.
Not only individual nanostructures were demonstrated for directional scattering of far-field light <cit.>, color-routing <cit.>, quantum emitter radiation steering <cit.>, electro-luminescence <cit.> or directional non-linear emission <cit.>.
Also metasurfaces for directional scattering and color-routing have been proposed <cit.>.
Such photonic devices imply applications like on-chip light routing and de-multiplexing <cit.> or sensing <cit.>.







To obtain functionalities such as directional scattering of light, a careful design of the nano-antenna is crucial. 
In the conventional approach a reference geometry is chosen from qualitative considerations whose properties are subsequently studied systematically. However, imposing a geometry from the start risks excluding interesting solutions to the problem.
A promising alternative can be evolutionary optimization (EO) strategies. 
By mimicking natural selection through a cycle of reproduction, evaluation and selection (see Fig. <ref>(a)), EO is able to find global optima for complex non-analytical problems <cit.>.
EO has recently been applied very successfully to various problems in nano-optics like maximization of local field enhancement <cit.> or of optical resonances and structural color <cit.>.
Other applications include the design of carbon nanotube field-emission sources <cit.>, dielectric anti-reflection coatings <cit.>, optical cloaking <cit.>, multi focal-spot flat lenses <cit.> or photonic power-splitters <cit.>.

In this paper we use EO to design directional optical antennas. 
We couple an evolutionary algorithm to full-field electro-dynamical simulations using the Green dyadic method (GDM).
With a randomly initialized and almost unconstrained geometric model, our approach reproducibly yields antennas resembling the classical radio-frequency (RF) Yagi-Uda layout <cit.>.
We find, that the working principle of these plasmonic antennas is indeed equivalent to their RF counter-parts.
We demonstrate the versatility of our method on different directional scattering problems such as varying angles of scattering, double-focusing, or quantum emitter steering. 
We finally show that evolutionary optimization can yield general design rules and furthermore can be used to identify and analyze physical limitations in nanoantenna design.







§ OPTIMIZATION METHOD, MODEL AND PROBLEM


Evolutionary optimization is a heuristic approach to search for an optimum solution of a complex problem, which cannot be optimized by classical means such as steepest descent methods. 
EO considers a certain number of randomly generated parametersets describing a problem. Parameters may be for instance the positions or sizes of sub-constituents of a nanostructure.
The ensemble of parametersets is called the population.
In a first step, each parameterset (also called an individual) is evaluated, in our case through a numerical simulation of the optical scattering (more details below).
The evaluation criterion is in general referred to as the fitness function.
In a selection step, the best solutions of the population (i.e. with the highest fitness values) are chosen for the following reproduction step.
Therein, randomly chosen parameters are mutated (randomly modified) and interchanged between individuals, mimicking evolutionary processes in nature and yielding a new population of parametersets.
This cyclic process is repeated multiple times (see Fig. <ref>(a)) until some stop criterion like a timelimit is met, leading to a close-to-optimum solution to the problem.

As EO algorithm we use the implementation of self-adaptive differential evolution (“jDE”) in the “pyGMO” toolkit <cit.>.
It takes as input the directionality of scattering, calculated using “pyGDM” <cit.>, our own python implementation of the Green Dyadic Method (GDM) <cit.>. 
By including a surface propagator we can take into account a glass substrate <cit.>.
Details on our EO-GDM technique can be found elsewhere <cit.>.
We want to note, that any numerical approach for solving Maxwell's equations can be used together with EO (possible other methods are for example finite difference time domain <cit.> or mode expansion techniques <cit.>).
With the GDM, we use a volume integral method which is particularly convenient when a relatively small nanostructure is placed in a large domain, where the GDM provides fast convergence.


A scheme of the geometric model to be optimized is shown in Fig. <ref>(b).
40 gold blocks B_i, each 40 × 40 × 40nm^3 in size, are placed in air (n_env=1) on a glass substrate (n_s = 1.5). For gold, we use the tabulated dispersion from Ref. johnson_optical_1972.
Each of the blocks are modeled by 2 × 2 × 2 dipoles. 
A discussion on the discretization step and the elementary block size and shape can be found in the appendix (Fig. <ref>).
The optimization target parameters are the positions (x_i, y_i) of the blocks, which are bound to an area of 1 × 1µm^2.
Additionally, the positions of the blocks are constrained to lie on a grid with steps of the block-size (40nm), to avoid the problem of partially overlapping elements. 
This results in a total of 25× 25 possible positions.
We note that blocks with identical positions are treated as a single block, hence the amount of material is not strictly constant.
With 40 gold blocks, we have 80 free parameters and (25× 25)! / (25× 25 - 40)! ≈ 10^111 possible arrangements.


A plane wave with λ_0 = 800nm, linearly polarized along OX, is normally incident onto the gold blocks. 
Using the GDM, we calculate at first the average intensity per solid angle I_direct, scattered through a solid angle Ω_direct of polar and azimuthal dimensions Δφ = 30 ^∘ and Δθ = 45 ^∘, centered at φ = 45 ^∘ and θ = 0 ^∘ (diagonal between X and Z axis), as illustrated by the green window in Fig. <ref>(c).
Second we calculate I_rest, the average intensity per solid angle, scattered towards the remaining hemisphere (Ω_rest, red in Fig. <ref>(c)).
The goal of the EO is to maximize the ratio between scattering per solid angle to the target “window” and to the rest of the hemisphere:


    R_direct = ( _Ω_direct I(Ω) dΩ) / Ω_direct( _Ω_rest I(Ω) dΩ) / Ω_rest
    					= I_directI_rest  .









§ DIRECTIONAL SCATTERING OF A PLANE WAVE




We found that for small populations (N_pop≪ 100), the optimization does not always converge to the best observed directionality ratio of R_direct≈ 12 (10.8dB).
Small populations seem to offer insufficient genetic diversity and tend to get stuck in local extrema.
The improvement during 50 generations of an evolution with a population size of N_pop = 500 is shown in figure <ref>(a).
Blue dots correspond to the fitness of the individuals, the green and red lines indicate the best and average fitness, respectively.
The best solutions during successive iterations and their far-field backscattering patterns are shown in Fig. <ref>(b)-<ref>(g) and Fig. <ref>(h) shows the best solution after the 50^th generation.
To verify the convergence and reproducibility, we repeated the same optimization with random initial populations, yielding always similar antenna morphologies with directionality ratios R_direct≈ 10.8 dB (see Appendix, Fig. <ref>).






At first sight we observe that the optimized nano-antenna visually resembles a Yagi-Uda RF antenna <cit.>.
An analysis of the phase of the electric field E_x component reveals a difference of π between the central and the outer parts of the nanostructure (Fig. <ref>(i)).
We identify three main functional elements, as shown in Fig. <ref>(j):
A driving feed in the antenna center (blue), a reflector on the left (green) and a director on the right (red).
As we explain in more detail in the appendix (see also Fig. <ref>), we want to note that this is not a categorical analogy to RF antenna design, but rather a simplified picture giving an idea about the roles of the different parts of the structure.

The intrinsic phase differences from resonance de-tuning of the elements in combination with the adequate distances between reflector, feed and director, result in cancellation of the scattered fields in backward direction and in constructive interference in forward direction.
In fact, this working principle is equivalent to the RF Yagi-Uda antenna, except that in our case the scattering is occurring at an angle with respect to the antenna axis.
Actually it has been shown earlier, that RF antenna concepts can work also at optical frequencies, despite large plasmonic losses <cit.>.



Nano Yagi-Uda antennas have been demonstrated at several occasions <cit.>.
Through an analysis of the impact of a removal or insertion of a gold block at every possible position (Fig. <ref>(k)) we find that the possible further increase in directionality is small. Furthermore we observe that fusing single blocks in the reflector and director, or cutting through the feed element kills the directionality. 
This is because changing the length of an element would result in a spectral shift of its LSP resonance, which would break the correct relative phases and scattered field amplitudes required for directional scattering.











 §.§ Different directional scattering objectives



In a second step, we run several optimizations for directional scattering problems, shown in figure <ref>.
The top row polar plots indicate the target solid angles Ω_direct by green color, which are of varying size and position.
The bottom plots show top views of the optimized gold block arrangements.
The center plots show the scattering radiation patterns of the optimization outcome.
In order to demonstrate that directionality is not obtained mainly through the index contrast at the substrate, we also show optimizations with different target solid angles for a vacuum environment (n_env=1, Figs. <ref>(a)-<ref>(d)).
Figs. <ref>, <ref>(e) and <ref>(f) show optimizations for plasmonic structures on a dielectric substrate (n_s=1.5), where in <ref>(e) and <ref>(f) furthermore scattering into the substrate was maximized. In all other cases scattering goes towards the air hemisphere.

While EO yields good directionality ratios in all cases, we observe that the scattering intensity and scattering efficiency Q_scat (ratio of scattering and geometric cross sections) are generally very low. 
The highest scattering intensities are obtained for large target solid angles and small deflection angles (Fig. <ref>(a)) or for scattering into a substrate at the critical angle (Fig. <ref>(e)). But even those cases have very small scattering efficiencies well below 1.
As a first conclusion we record that it is a physically demanding task to steer an incident plane wave unidirectionally at large angles via a planar gold nanostructure.

For comparison we ran EOs maximizing the absolute backscattered intensity I_direct towards polar angles of 45^∘ and 90^∘, shown in Figs. <ref>(a) and <ref>(b), respectively.
In both cases, the largest absolute scattered intensity towards the target solid angle is found for a structure built of several resonant gold elements, giving 2-3 orders of magnitude higher overall scattering efficiencies. 
But most of the light is backscattered with no deflection, directionality is completely lost.
Still, the scattered intensity I_direct towards the target solid angles in Fig. <ref> is several times higher compared to Fig. <ref> and Fig. <ref>(b), where the ratio R_direct was maximized. The latter antennas are built from detuned elements, resulting in a very weak absolute scattering.
This supports our conclusion that using planar plasmonic nanostructures, it is physically almost impossible to redirect a normally incident plane wave towards the direction of its linear polarization. 
Large scattering intensities seem to imply that almost no directionality is obtained.

We note that Fig. <ref> shows structures on a substrate. We also performed the same optimization in a vacuum environment, yielding similar results.







We provide a visualization online (Visualization 1), showing the progress of EO of a large-area (5 × 5µm^2), intensity optimized scatterer (analog to Fig. <ref>(a)). It results in an assembly of several resonant nanoantennas, arranged in a grating-like way. Despite the large allowed area and consequently much more degrees of freedom, also this optimization fails to combine a good directionality ratio with high scattering intensity.



 §.§ Discussion: Weak scattering due to phase difference requirement


The weak scattering at high deflection angles is a direct result of the physical properties of the plane wave illumination. The phase of the illumination field is constant across the entire planar structure, which means that it is not possible to obtain phase differences only through positioning of plasmonic elements. 
On the other hand, the directional antenna requires a proper phase difference between antenna feed and director elements.
To achieve this, the EO algorithm choses very small blocks as director elements, which are driven by the plane wave below resonance (ω_res≫ω_0), resulting in a phase of -π / 2 between illumination and response.
As driving feed the optimizer choses a large rod having its LSP resonance in the infrared (ω_res≪ω_0), which gives a phase of +π / 2 with respect to the driving plane wave. 
In consequence, the feed and director structures have a relative phase of π with respect to each other but are also entirely off-resonance and thus scatter only very weakly. This results in an overall small scattering intensity, so most of the incident light passes just through the plasmonic structure without interaction. 
The weak interaction with the incident plane wave is confirmed by the very low absorption efficiencies of the optimized antennas in Figs 2 and 3.

As we see in Figs. <ref>(a) and <ref>(b), if the intensity is target of the optimization, the resulting structures look fundamentally different: The dispersed small blocks, that give the directionality-optimized antennas the fragmented look, disappear. The EO algorithm now uses all the available material to construct as many resonant nano-rods as possible, yielding a strong scattering cross section but very weak directionality.










We finally want to note, that a drawback of the GDM with regards to fabrication robustness is, that diagonally neighboring metallic blocks are electrically isolated in the simulations. 
In an actual sample however, there are high chances that diagonal neighbor blocks are electrically connected at their touching point due to unavoidable imperfections in the fabrication process. Such electric connections can have a significant impact on localized surface plasmon resonances and in consquence have the potential to change the optical response of the structure.
To assess the robustness of our results against fabricational imperfections, we ran simulations, connecting diagonally touching elements in the optimized antennas using additional dipoles, which didn't have a severe impact in these cases.
To generally overcome this shortcoming, dielectric / plasmonic hybrid structures may be designed which can yield inherent phase differences without requiring precise tailoring of the relative LSP resonances <cit.>.












§ QUANTUM EMITTER DIRECTIONAL ANTENNA





As found above, plane wave directional scattering at high scattering efficiency is physically difficult to obtain. 
On the other hand, it has been demonstrated that plasmonic nanostructures can be used as efficient directional antennas for quantum emitters, being also able to boost the decay of a quantum transition thanks to the Purcell effect <cit.>.
In fact, in the case of a local dipole transition as light source, there is not the problem that the illumination excites every part of the antenna simultaneously and with constant phase. This makes it significantly easier to find a solution.






Following the configuration reported in curto_unidirectional_2010, we replace the plane wave illumination by a dipole transition along OY, oscillating at λ_0 = 800nm, positioned at the origin of the coordinate system. 
We prohibit the EO algorithm to put a gold block at this location, which results in a minimum distance of 20nm between gold and emitter. The emission is to be maximized towards OX along the critical angle into the n_s=1.5 substrate.
The EO converges towards an antenna with a directionality ratio as high as R_direct≈ 13.9dB, shown in Fig. <ref>(a).
Ref. curto_unidirectional_2010 uses a slightly different definition of directionality, therefore we simulate the antenna “YU145” of Ref. curto_unidirectional_2010 with the GDM, using the same emitter-feed distance of 20nm as for the EO. We find a directionality of R_direct, ref.≈ 9.2dB (Fig. <ref>(c)).
While this is lower compared to the optimized antenna, we see in Fig. <ref>(d) that the emission intensity of the optimized antenna is around two orders of magnitude weaker.
We re-run the EO using as optimization target the absolute intensity of directed emission (Fig. <ref>(b)).
While directionality is moderately reduced to 11.9dB, the signal from the intensity-optimized directional nanoantenna (orange line in Fig. <ref>(d)) is now almost one order of magnitude higher than that of the reference. 
Having a closer look at the geometry in Fig. <ref>(b), we observe that the EO algorithm designed the antenna such that the emitter is surrounded by gold to maximize the Purcell effect, which strongly boosts the emitted intensity.
Finally we find that an increased directionality comes at the cost of a reduced spectral full width at half maximum (FWHM) of R_direct, hence leads to a narrower operation window (see inset Fig. <ref>(d)).











 §.§ Dipole emitter orientation


To study the role of the orientation of the quantum emitter's dipole transition, we run antenna optimizations for X, Y and Z oriented dipoles, shown in Figs. <ref>(a)-<ref>(c).
To keep the symmetry of the system as high as possible, we assume a vacuum environment and maximize scattering towards OX (solid angle as used in Fig. <ref>(b)).
In case of X- and Y-oriented dipoles (Figs. <ref>(a)-<ref>(b)), the optimization organizes the gold blocks in a compact arrangement, forming a strongly directional antenna.
On the other hand, for a dipole pointing out of the structure plane (Z direction), the limited height of the planar gold structure (H=40nm) impedes obtaining a significant optical response of the gold blocks at the emission wavelength of λ_0=800nm, hence directionality remains weak (Fig. <ref>(c)).




 §.§ Antenna symmetry


In contrast to the intuitive expectation that symmetric antennas should perform best, EO finds mostly rather fragmented antennas of low symmetry.
To assess the role of symmetry, we re-run the quantum emitter antenna optimizations for different dipole-orientation (Figs. <ref>(a)-<ref>(c)), imposing a mirror symmetry around the X axis.
The results shown in Fig. <ref>(d)-<ref>(f) suggest, that symmetric antennas give a similar performance in the cases of dipole transitions aligned along Y and Z, as depicted in Fig. <ref>(e), respectively (f), even though the symmetric designs do not yield noticeable improvement over an unconstrained geometric model.

In the case of an X-oriented dipole however, the symmetric EO completely fails to design a directional antenna (Fig. <ref>(d)).
By inspecting the antenna in Fig. <ref>(a), we notice that it is designed such that the polarization of the X-oriented dipole emitter is converted by a gold element to a plasmonic current along Y, effectively generating an accordingly oriented plasmonic dipole emitter.
The scattered light of this LSP dipolar resonance is then guided towards the target solid angle via a Yagi-Uda-type assembly of structures around the feed.
By imposing an axial symmetry at the dipole position, this local polarization conversion is physically impossible. 
Instead, in Fig. <ref>(d) the EO algorithm designs a quadrupolar feed element, directing the light close to the desired solid angle. However, due to symmetry, emission is zero towards X.











§ CONCLUSIONS







In conclusion, we demonstrated how evolutionary optimization can be used to design directional plasmonic nanoantennas.
In our approach, within the limits of possible geometries built with 40 equal gold blocks, the design of the structures is completely free.

We showed, that EO automatically finds a nanostructure layout which corresponds to well known design-principles from radio-frequency antennas.
By comparing the performance of directionality- and intensity-optimized planar gold nanostructures, we found that high directionality ratios for scattering of normal incidence plane waves to large angles comes necessarily at the cost of weak scattering intensities.
We demonstrated that in contrast to weak directional plane wave scattering, EO can design very efficient directional nanoantennas for quantum emitters. 
Via an appropriate fitness function, EO finds an antenna which not only yields very good directionality but also exploits at best the Purcell effect for quantum emitter signal enhancement. 
Through an analysis of the EO results for different fitness functions we furthermore found, that high directionality ratios come at the cost of spectrally less broadband working windows.
Finally, through an EO-based analysis of symmetry conditions we found that, dependent on the orientation of the quantum emitter dipole transition, it can be necessary to design and asymmetric nano-antenna for good performance.

We conclude, that from the results of EO general design rules can be derived for the tailoring of photonic nanostructures. 
EO can also be used as a benchmark tool in nano-photonic design and can reveal intrinsic physical limitations of nano-optical configurations.
We foresee that EO has manifold further promising applications beyond directional scattering, for examples in field enhanced spectroscopy or non-linear nanophotonics.






§ SUPPLEMENTARY MATERIALS


Visualization 1: Movie showing the convergence process of EO for an intensity-optimized, plane wave illuminated, directional optical antenna.





§ APPENDIX




 §.§ Reproducibility

To test the convergence and the reproducibility of the EO, we re-run the optimization shown in main text Fig. <ref> several times with different, randomly initialized starting populations.
As shown in Fig. <ref>, the final antennas converge always to the same maximum directionality ratio of R_direct≈ 12.
The geometries of the solutions are in all cases similar, leading to the same phase distribution.








 §.§ Functional elements of plane wave directional antenna


It is clear that the parasitic elements around the feed lead to the directional scattering effect in the plane-wave antennas (as demonstrated in Fig. <ref>(a) and <ref>(b) below). 
On the other hand it is difficult to ascribe to specific clusters of parasitic elements the unambiguous role of a “reflector” or a “director”.
To assess whether it is justified to call the agglomeration of blocks on the left of the feed in Fig. <ref> an equivalent of a “reflector” and the agglomeration of blocks on the right a “director”, we simulated the feed alone, the feed + “reflector” and the feed + “director”, shown in Fig. <ref>(b)-<ref>(d).





While we observe no directional scattering for the isolated feed (Fig. <ref>(b)), the reflector as well as director alone suffice to impose a directional effect on the feed element (Figs. <ref>(c)-<ref>(d)), in analogy to what would be expected in an RF Yagi-Uda antenna.
On the other hand, if only a part of for example the reflector cluster is used, the direction of scattering becomes tilted with respect to the target angle (Fig. <ref>(e)). The full ensemble of “reflector” blocks is necessary to obtain the desired scattering effect.
We conclude that it is adequate to consider the left and right agglomerations of parasitic elements equivalent to a “reflector” and a “director”. 

As a final note, we want to stress that this is not a categorical analogy. 
A direct adaptation of the structure to RF antennas is not easily possible, we rather consider our analogy a simplified description of the role of the complex arrangements of parasitic elements.






 §.§ Discretization, size and shape of elementary blocks


In principle any kind of constituents can be used with the EO-GDM method to compose the structure.
However, larger building blocks or elementary blocks of some complex shape constrain the generality of possible geometries. 
On the other hand, small blocks allow to perform a more complete search of the solution space, but at the cost of slower convergence.
The choice of size (and shape) of the elementary blocks is therefore always a compromise between how flexible the geometric model should be and the speed of convergence of the optimization.
If constraints to the attainable diversity of geometries are acceptable, only few but larger building blocks and/or blocks of more specific shape like cuboids, rods or polygonal structures might be chosen as elementary unit.
In view of a possible experimental realization of the structures, it makes also sense to use building blocks which are not smaller then the structural resolution that can be obtained in fabrication.






Concerning the discretization, we used steps of 20nm which gives good optimization speed because the full structure consists of only 40× (2× 2× 2) = 320 dipoles.
In order to verify that this rather coarse discretization is a good approximation, we re-calculated the structure shown in figure 2 with finer step-sizes.
As shown in Fig. <ref>, simulations with steps of 20nm (2× 2× 2 dipoles per block), 10nm  (4× 4× 4 dipoles per block) and 6.67nm (6× 6× 6 dipoles per block) yield very similar results, so our approximation using steps of 20nm is good.

















§ FUNDING

This work was supported by Programme Investissements d'Avenir under the program ANR-11-IDEX-0002-02, reference ANR-10-LABX-0037-NEXT and by the computing facility center CALMIP of the University of Toulouse under grant P12167.
PW received funding by the German Research Foundation (DFG) through a Research Fellowship (WI 5261/1-1). OM acknowledges support through EPSRC grant EP/M009122/1.




§ ACKNOWLEDGMENTS

We gratefully thank N. Bonod for fruitful discussions.
All data supporting this study are openly available from the University of Southampton repository (DOI: 10.5258/SOTON/D0992).
























10
	
	
	muhlschlegel_resonant_2005
	P. Mühlschlegel, H.-J. Eisler, O. J. F. Martin, B. Hecht, and D. W. Pohl,
	Resonant Optical Antennas, Science
	308, 1607–1609 (2005).
	
	tan_plasmonic_2014
	S. J. Tan, L. Zhang, D. Zhu, X. M. Goh, Y. M. Wang, K. Kumar, C.-W. Qiu, and
	J. K. W. Yang, Plasmonic Color Palettes for Photorealistic
				Printing with Aluminum Nanostructures, Nano
			Letters 14, 4023–4029 (2014).
	
	wiecha_polarization_2017
	P. R. Wiecha, L.-J. Black, Y. Wang, V. Paillard, C. Girard, O. L. Muskens, and
	A. Arbouet, Polarization conversion in plasmonic nanoantennas for
		metasurfaces using structural asymmetry and mode hybridization,
	Scientific Reports 7, 40906 (2017).
	
	valev_chirality_2013
	V. K. Valev, J. J. Baumberg, C. Sibilia, and T. Verbiest, Chirality
		and Chiroptical Effects in Plasmonic Nanostructures:
		Fundamentals, Recent Progress, and Outlook,
	Advanced Materials 25, 2517–2534 (2013).
	
	baffou_heat_2009
	G. Baffou, R. Quidant, and C. Girard, Heat generation in plasmonic
		nanostructures: Influence of morphology, Applied
			Physics Letters 94, 153109 (2009).
	
	girard_designing_2018
	C. Girard, P. R. Wiecha, A. Cuche, and E. Dujardin, Designing
		thermoplasmonic properties of metallic metasurfaces,
	Journal of Optics 20, 075004 (2018).
	
	kauranen_nonlinear_2012
	M. Kauranen and A. V. Zayats, Nonlinear plasmonics,
	Nature Photonics 6, 737–748 (2012).
	
	mesch_nonlinear_2016
	M. Mesch, B. Metzger, M. Hentschel, and H. Giessen, Nonlinear
		Plasmonic Sensing, Nano Letters 16,
	3155–3159 (2016).
	
	vercruysse_directional_2014
	D. Vercruysse, X. Zheng, Y. Sonnefraud, N. Verellen, G. Di Martino, L. Lagae,
	G. A. E. Vandenbosch, V. V. Moshchalkov, S. A. Maier, and P. Van Dorpe,
	Directional Fluorescence Emission by Individual V-Antennas
				Explained by Mode Expansion, ACS Nano
	8, 8232–8241 (2014).
	
	wersall_directional_2014
	M. Wersäll, R. Verre, M. Svedendahl, P. Johansson, M. Käll, and
	T. Shegai, Directional Nanoplasmonic Antennas for
		Self-Referenced Refractometric Molecular Analysis,
	The Journal of Physical Chemistry C 118,
	21075–21080 (2014).
	
	vercruysse_unidirectional_2013
	D. Vercruysse, Y. Sonnefraud, N. Verellen, F. B. Fuchs, G. Di Martino,
	L. Lagae, V. V. Moshchalkov, S. A. Maier, and P. Van Dorpe,
	Unidirectional Side Scattering of Light by a
		Single-Element Nanoantenna, Nano Letters
	13, 3843–3849 (2013).
	
	abass_insights_2016
	A. Abass, P. Gutsche, B. Maes, C. Rockstuhl, and E. R. Martins,
	Insights into directional scattering: From coupled dipoles to
		asymmetric dimer nanoantennas, Optics Express
	24, 19638 (2016).
	
	yao_controlling_2016
	K. Yao and Y. Liu, Controlling Electric and Magnetic
				Resonances for Ultracompact Nanoantennas with Tunable
				Directionality, ACS Photonics  (2016).
	
	shegai_bimetallic_2011
	T. Shegai, S. Chen, V. D. Miljković, G. Zengin, P. Johansson, and
	M. Käll, A bimetallic nanoantenna for directional colour
		routing, Nature Communications 2, 481
	(2011).
	
	curto_unidirectional_2010
	A. G. Curto, G. Volpe, T. H. Taminiau, M. P. Kreuzer, R. Quidant, and N. F. van
	Hulst, Unidirectional Emission of a Quantum Dot Coupled to a
		Nanoantenna, Science 329, 930–933
	(2010).
	
	maksymov_optical_2012
	I. S. Maksymov, I. Staude, A. E. Miroshnichenko, and Y. S. Kivshar,
	Optical Yagi-Uda nanoantennas,
	Nanophotonics 1, 65–81 (2012).
	
	hancu_multipolar_2014
	I. M. Hancu, A. G. Curto, M. Castro-López, M. Kuttge, and N. F. van
		Hulst, Multipolar Interference for Directed Light
				Emission, Nano Letters 14, 166–171
	(2014).
	
	dregely_imaging_2014
	D. Dregely, K. Lindfors, M. Lippitz, N. Engheta, M. Totzeck, and H. Giessen,
	Imaging and steering an optical wireless nanoantenna link,
	Nature Communications 5, 4354 (2014).
	
	ramezani_hybrid_2015
	M. Ramezani, A. Casadei, G. Grzela, F. Matteini, G. Tütüncüoglu,
	D. Rüffer, A. Fontcuberta i Morral, and J. Gómez Rivas,
	Hybrid Semiconductor Nanowire–Metallic Yagi-Uda
				Antennas, Nano Letters 15, 4889–4895
	(2015).
	
	gurunarayanan_electrically_2017
	S. P. Gurunarayanan, N. Verellen, V. S. Zharinov, F. James Shirley, V. V.
	Moshchalkov, M. Heyns, J. Van de Vondel, I. P. Radu, and P. Van Dorpe,
	Electrically Driven Unidirectional Optical Nanoantennas,
	Nano Letters 17, 7433–7439 (2017).
	
	r.__kullock_directed_2018
	R. Kullock, P. Grimm, M. Ochs, and B. Hecht, Directed emission
		by electrically driven optical antennas, in Quantum Sensing and
		Nano Electronics and Photonics XV,  vol. 10540 (2018), pp.
	1054012–10540–9.
	
	xiong_compact_2016
	X. Y. Z. Xiong, L. J. Jiang, W. E. I. Sha, Y. H. Lo, and W. C. Chew,
	Compact Nonlinear Yagi-Uda Nanoantennas,
	Scientific Reports 6, 18872 (2016).
	
	lindfors_imaging_2016
	K. Lindfors, D. Dregely, M. Lippitz, N. Engheta, M. Totzeck, and H. Giessen,
	Imaging and Steering Unidirectional Emission from Nanoantenna
				Array Metasurfaces, ACS Photonics 3,
	286–292 (2016).
	
	yan_fano-resonance-assisted_2017
	C. Yan, K.-Y. Yang, and O. J. F. Martin, Fano-resonance-assisted
		metasurface for color routing, Light: Science &
			Applications 6, e17017 (2017).
	
	guo_highbit_2017
	R. Guo, M. Decker, F. Setzpfandt, X. Gai, D.-Y. Choi, R. Kiselev,
	A. Chipouline, I. Staude, T. Pertsch, D. N. Neshev, and Y. S. Kivshar,
	High–bit rate ultra-compact light routing with
		mode-selective on-chip nanoantennas, Science
			Advances 3, e1700007 (2017).
	
	sivanandam_introduction_2008
	S. Sivanandam and S. Deepa, Introduction to Genetic Algorithms
	(Springer, Heidelberg, 2008).
	
	forestiere_particle-swarm_2010
	C. Forestiere, M. Donelli, G. F. Walsh, E. Zeni, G. Miano, and L. Dal Negro,
	Particle-swarm optimization of broadband nanoplasmonic arrays,
	Optics Letters 35, 133 (2010).
	
	feichtner_evolutionary_2012
	T. Feichtner, O. Selig, M. Kiunke, and B. Hecht, Evolutionary
		Optimization of Optical Antennas, Physical
			Review Letters 109, 127701 (2012).
	
	forestiere_inverse_2016
	C. Forestiere, Y. He, R. Wang, R. M. Kirby, and L. Dal Negro, Inverse
		Design of Metal Nanoparticles' Morphology,
	ACS Photonics 3, 68–78 (2016).
	
	ginzburg_resonances_2011
	P. Ginzburg, N. Berkovitch, A. Nevet, I. Shor, and M. Orenstein,
	Resonances On-Demand for Plasmonic Nano-Particles,
	Nano Letters 11, 2329–2333 (2011).
	
	bigourdan_design_2014
	F. Bigourdan, F. Marquier, J.-P. Hugonin, and J.-J. Greffet, Design of
		highly efficient metallo-dielectric patch antennas for single-photon
		emission, Optics Express 22, 2337 (2014).
	
	mirzaei_superscattering_2014
	A. Mirzaei, A. E. Miroshnichenko, I. V. Shadrivov, and Y. S. Kivshar,
	Superscattering of light optimized by a genetic algorithm,
	Applied Physics Letters 105, 011109 (2014).
	
	wiecha_evolutionary_2017
	P. R. Wiecha, A. Arbouet, C. Girard, A. Lecestre, G. Larrieu, and V. Paillard,
	Evolutionary multi-objective optimization of colour pixels based on
		dielectric nanoantennas, Nature Nanotechnology
	12, 163–169 (2017).
	
	gonzalez-alcalde_optimization_2018
	A. K. González-Alcalde, R. Salas-Montiel, H. Mohamad, A. Morand,
	S. Blaize, and D. Macías, Optimization of all-dielectric
		structures for color generation, Applied Optics
	57, 3959–3967 (2018).
	
	chen_optimal_2007
	P. Y. Chen, C. H. Chen, J. S. Wu, H. C. Wen, and W. P. Wang, Optimal
		design of integrally gated CNT field-emission devices using a genetic
		algorithm, Nanotechnology 18, 395203
	(2007).
	
	good_general_2016
	B. L. Good, S. Simmons, and M. Mirotznik, General optimization of
		tapered anti-reflective coatings, Optics Express
	24, 16618 (2016).
	
	mirzaei_all-dielectric_2015
	A. Mirzaei, A. E. Miroshnichenko, I. V. Shadrivov, and Y. S. Kivshar,
	All-Dielectric Multilayer Cylindrical Structures for
		Invisibility Cloaking, Scientific Reports
	5, 9574 (2015).
	
	hu_evolutionary_2017
	J. Hu, X. Ren, A. N. Reed, T. Reese, D. Rhee, B. Howe, L. J. Lauhon, A. M.
	Urbas, and T. W. Odom, Evolutionary Design and Prototyping
		of Single Crystalline Titanium Nitride Lattice Optics,
	ACS Photonics 4, 606–612 (2017).
	
	molesky_inverse_2018
	S. Molesky, Z. Lin, A. Y. Piggott, W. Jin, J. Vucković, and A. W.
	Rodriguez, Inverse design in nanophotonics,
	Nature Photonics 12, 659 (2018).
	
	yagi_beam_1928
	H. Yagi, Beam Transmission of Ultra Short Waves,
	Proceedings of the Institute of Radio Engineers
	16, 715–740 (1928).
	
	biscani_global_2010
	F. Biscani, D. Izzo, and C. H. Yam, A Global Optimisation Toolbox
		for Massively Parallel Engineering Optimisation,
	arXiv:1004.3824 [cs, math]  (2010).
	
	wiecha_pygdmpython_2018
	P. R. Wiecha, pyGDM—A python toolkit for full-field
		electro-dynamical simulations and evolutionary optimization of
		nanostructures, Computer Physics Communications
	233, 167–192 (2018).
	
	girard_shaping_2008
	C. Girard, E. Dujardin, G. Baffou, and R. Quidant, Shaping and
		manipulation of light fields with bottom-up plasmonic structures,
	New Journal of Physics 10, 105016 (2008).
	
	girard_generation_1995
	C. Girard, A. Dereux, O. J. F. Martin, and M. Devel, Generation of
		optical standing waves around mesoscopic surface structures: Scattering
		and light confinement, Physical Review B
	52, 2889–2898 (1995).
	
	johnson_optical_1972
	P. B. Johnson and R. W. Christy, Optical Constants of the Noble
				Metals, Physical Review B 6, 4370–4379
	(1972).
	
	hofmann_design_2007
	H. F. Hofmann, T. Kosako, and Y. Kadoya, Design parameters for a
		nano-optical Yagi–Uda antenna, New
			Journal of Physics 9, 217 (2007).
	
	novotny_effective_2007
	L. Novotny, Effective Wavelength Scaling for Optical
				Antennas, Physical Review Letters 98,
	266802 (2007).
	
	zhang_unidirectional_2019
	T. Zhang, J. Xu, Z.-L. Deng, D. Hu, F. Qin, and X. Li, Unidirectional
		Enhanced Dipolar Emission with an Individual Dielectric Nanoantenna,
	Nanomaterials 9, 629 (2019).
	
	francs_plasmonic_2016
	G. C. des Francs, J. Barthes, A. Bouhelier, J. C. Weeber, A. Dereux, A. Cuche,
	and C. Girard, Plasmonic Purcell factor and coupling efficiency
		to surface plasmons. Implications for addressing and controlling optical
		nanosources, Journal of Optics 18, 094005
	(2016).
	









