Abstract
Given an arbitrarily complicated object, it is often difficult to say immediately how it interacts with a specific illumination. Optically small objects, e.g., spheres, can often be modeled as electric dipoles, but which multipole moments are excited for larger particles possessing a much more complicated shape? The Tmatrix answers this question, as it contains the entire information about how an object interacts with any electromagnetic illumination. Moreover, a multitude of interesting properties can be derived from the Tmatrix such as the scattering cross section for a specific illumination and information about symmetries of the object. Here, we present a method to calculate the Tmatrix of an arbitrary object numerically, solely by illuminating it with multiple plane waves and analyzing the scattered fields. Calculating these fields is readily done by widely available tools. The finite element method is particularly advantageous, because it is fast and efficient. We demonstrate the Tmatrix calculation at four examples of relevant optical nanostructures currently at the focus of research interest. We show the advantages of the method to obtain useful information, which is hard to access when relying solely on full wave solvers.
Introduction
Recent advances in nanofabrication technology made the creation of large volumes of particles with complicated geometries possible [16]. The latter constitute the base for nanomaterials with advanced properties [7]. This also triggered the need for efficient computational tools to back up experimental findings with simulations that allow to understand the underlying principles that cause the properties of the respective nanomaterials [810]. With current possibilities it is fairly straightforward to calculate the scattered fields of an arbitrary object illuminated with a plane wave or more complicated illumination scenarios. They all rely on the numerical solution of Maxwell’s equations while considering a given distribution of material in space. Examples for such numerical routines are finite element solvers [1113], finitedifference timedomain methods [1416], discrete dipole approximation [17,18], or similar tools. This yields valuable information about the single structure. On the other hand, there are efficient methods to calculate the scattering of large clusters of spherical particles, for example the extended Mie theory [19,20]. This consideration of many particles is the prerequisite to study the emergence of properties if many of these particles are brought together to form artificial nanomaterials.
The principles of these two distinct approaches can be combined to the Tmatrix method [2123]. Using this method, the Tmatrix, which contains the information about how an isolated object scatters the electromagnetic field, is calculated and subsequently used in a multiplescattering algorithm, similar to the multiMie Method [20,24], to calculate the scattered fields of a large cluster of arbitrary particles.
There are several analytical and semianalytical ways to calculate the Tmatrix of simple particles with high symmetries [21,2527]. Here, we present a method to calculate the Tmatrix of an arbitrarily complex particle. Our approach only requires the object to be illuminated in a sequence of simulations with multiple plane waves. From these simulations the scattering coefficients and, ultimately, the Tmatrix, are deduced from the scattered fields. These fields can be calculated with established methods, such as the finite element method but, in general, any other of the aforementioned methods can be considered as well. The resulting Tmatrix is then independent of the form and direction of the external illumination. It only depends on the material and geometry of the object. The scattering coefficients that can then be calculated for any arbitrary illumination are vital for the calculation of all further quantities, such as the scattered fields and the scattering cross section for a specific illumination. The availability of the Tmatrix simplifies particularly the analysis of quantities that require the evaluation of the same response of the object to a larger number of different illumination fields. A referential example therefore is the calculation of the force map for a complicated illumination in space exerted on the particle. There, the response of the same particle to the field at each spatial location of the illumination needs to be calculated. This can be a tedious task if every time a new full wave solution to Maxwell’s equations is required. The situation is much more simplified if the Tmatrix of the object is known. The coefficients can also be used to determine the strength of the individual multipole contributions to the scattering for example to check if the dipole approximation is applicable. This is particularly useful in the context of the homogenization, i.e., the assignment of effective material parameters to a medium. We emphasize the fact that when only the optical response of a given particle to a given illumination is needed, it does not make much sense to calculate the Tmatrix first, as it requires the solution of multiple fullwave problems. However, whenever more complicated objects are of interest, or if the response of a given structure to many different illuminations is needed, as it is often required in the calculation of force and torque, it is expected that the computation of the Tmatrix pays off. Furthermore, as we will point out when we investigate the example cases, the Tmatrix itself is of scientific value, as it contains all the information available concerning how an object interacts with light.
In the following section we outline the multipole decomposition of the fields [28]. This is the theoretical basis for the Tmatrix calculation. The numerical calculation of the Tmatrix is demonstrated in the next section. Finally, we provide four numerical examples in the form of interesting but simple objects the Tmatrix of which we calculate to demonstrate the usefulness of the method.
Multipole expansion of electromagnetic fields
We want do decompose the fields produced by scattering off of an illuminated object into contributions from different multipoles. We start by decomposing the total external fields surrounding the object into incident and scattered fields
Note that this is always possible for theoretical considerations, but experimentally a strict distinction is not possible. In the following we largely skip the space and frequency dependency of the fields for simplicity but these dependencies are always assumed.
Then, the fields are expanded into vector spherical harmonic functions
with the normalization factor
Here, the index l represents the underlying function: l = 1,2 corresponds to spherical Bessel functions of first z^{(1)} = j_{n} and second kind z^{(2)} = y_{n}, respectively, while l = 3,4 corresponds to Hankel functions of first and second kind , respectively [29,30]. The incident field has to be finite in the center of the coordinate system. Thus, we have to take Bessel functions and the expansion reads as
The scattered field needs to satisfy the outgoing Sommerfeld radiation condition and is expanded with Hankel functions
Here, k is the wave number in the surrounding medium, which is characterized by the permittivity ε(ω) and the permeability μ(ω). The complex expansion coefficients a_{nm}(ω) and b_{nm}(ω) are called scattering coefficients and the coefficients p_{nm}(ω) and q_{nm}(ω) are called incident coefficients. Together they contain all relevant information about the interaction of the particle with the given illumination. The physical meaning of the individual terms is also straightforward. For example in the expansion of the scattered field, the first term on the right hand side expresses the field scattered due to the induced electric multipole moments. The second term expresses the field scattered due to the induced magnetic multipole moments.
In order to be able to perform numerical calculations, we need to truncate the infinite sums in all expansions to a finite number N. This number is the multipole order that we take into account, N = 1 corresponding to dipoles, N = 2 to quadrupoles and so on. Thus, by choosing the maximum order we want to consider, we can strongly influence the computation time and final accuracy.
The scattering coefficients can be used to calculate the total scattering cross section of the system
Here, η is a scaling factor, depending on the illumination. For a plane wave we have η = 1/(k^{2}E_{0}^{2}), where E_{0} is the amplitude of the incident plane wave. In order to obtain more information, e.g., about a resonance, we can modify the equation to calculate the contributions of each multipolar order
where n^{′} is the multipolar index and c_{nm} stands for a_{nm} if we want to consider the electric part or b_{nm} if we consider the magnetic part, respectively. For example the electric and magnetic dipole contributions, which are usually the strongest contributions for electromagnetically small particles, read as
Thus, if we know the scattering coefficients of an object, we know exactly how the different multipoles contribute to the total scattering. That is particularly valuable if we want to design a specific multipole resonance of a nanoparticle or if we want to investigate a given resonance, e.g., from experimental measurements.
Tmatrix calculation
In the previous section we showed the decomposition of the fields into scattering and incident coefficients. Now we introduce a link between them that represents the interaction of the object with the illumination: the Tmatrix
Here a and b are concatenated vectors that contain the scattering coefficients of the outgoing wave and p and q are concatenated vectors containing the coefficients of the incident field.
The Tmatrix contains the entire near and far field information on how a specific object interacts with any illumination. It depends on the geometry and material composition of the structure in question. The Tmatrix of a sphere is known analytically as the Mie coefficients. They form the Tmatrix if they are ordered at the diagonal of the matrix. For arbitrary objects, however, the Tmatrix is generally dense and is typically ordered in the following way
where T_{ee} describes coupling of electric incidence modes to electric scattered modes, T_{mm} describes the magnetic coupling and the other parts correspond to mixing of electric and magnetic modes. These parts of the Tmatrix can be ordered to show the coupling between the multipolar modes
Then Tnn' is the part that represents the coupling of order n′ to n. For a sphere, T_{11} would contain the first Mie coefficient a_{1} on the diagonal, T_{22} the second and so on.
Please note that the vector spherical harmonic functions and in their presented form are normalized to the unit sphere. Therefore, also the Tmatrix is normalized. In many cases it is crucial to use a normalized form of the Tmatrix, e.g., when applying rotations [23,31] or for comparing the absolute value of different entries.
To perform simulations and to analyze a given object we need to calculate the Tmatrix numerically. There are several established methods to calculate it, the first being the extended boundary condition method, originally introduced by Waterman [21,32,33]. Here, only homogeneous and isotropic particles can be considered and the method is quite time and resource consuming. A similar, more advanced possibility is the point matching method [34,35]. There, the coefficients of the incident, scattered and internal fields are related, but without the expensive surface integrations needed for the extended boundary condition method. Another established strategy to compute the Tmatrix of an arbitrary object is to excite it with pure vector spherical harmonic functions to extract exactly one line of the Tmatrix per illumination [22]. There, the fields are calculated with an integral equation solver. This method has the advantage that the incident field can be expressed in an exact way with a single incident coefficient. But the actual implementation of such an illumination with available tools is cumbersome. Similar methods have been proposed to retrieve the polarizability tensor of an object [3639]. These techniques provide valuable information with comparatively little effort, but are restricted to dipolar objects.
Here, we propose a method that is easy to implement practically and, depending on the applied Maxwell solver, very efficient. We illuminate the particle with several different plane waves and extract the expansion coefficients by projecting the scattered fields onto vector spherical harmonic functions. To get the scattered field of the object for every plane wave, however, we need an additional tool that solves Maxwell’s equations for an arbitrary particle under plane wave illumination. Such tools are widely available. We chose here the finite element solver JCMsuite [12] because it is especially fast for the simulation of several different illuminations of an object at the same frequency. But in general, any fullwave solver that provides the scattered fields of an object could be used.
The scattering coefficients can be obtained from the scattered fields by Equation 13 and similarly the incident coefficients by Equation 14.
The fields are evaluated at a virtual sphere of fixed radius R that completely encloses the object. Analytically, the chosen radius has no influence on the result. Yet, in numerical realizations the choice may effect the result due to the finite sampling. It is advisable to use a small R, to keep the computational domain small, without touching the actual object. Note that we need to take the full details of the scattered field into account.
Now we get a set of and corresponding for each illumination k. Then we can construct a system of equations to calculate the numerical Tmatrix T_{K} from Equation 10
which we invert to find the Tmatrix. In order to get an invertible system, we need to calculate the response from K illuminations, where K is equal to or greater than the rank d of the Tmatrix. This rank depends on the multipole order N we want to take into account
Using more illuminations is not strictly necessary, but will improve the accuracy of the result, as we will see in the following examples. The inversion of the rectangular matrix can then be done for example with QR decomposition [40]. Note that the resulting Tmatrix is always a square matrix in our formalism.
Results and Discussion
Single sphere
As a first check to prove that the method delivers correct results, we investigate the most simple object: a single isolated sphere. We have the correct analytical solution for a single sphere available. The Tmatrix has only the known Miecoefficients on the diagonal and all other entries are zero.
Let us consider a single dielectric sphere with a radius of 100 nm and a relative permittivity of 16 in vacuum. The Miecoefficients of this sphere for the first two orders are nonnegligible at 600 THz. Such a highpermittivity sphere is nowadays at the focus of interest since it sustains a notable electric and magnetic dipole moment in the visiblelight range.
To calculate the Tmatrix numerically we illuminate the sphere with the appropriate number of plane waves K. The wave vectors are chosen, such that they are evenly distributed in all directions. They can be considered as normal vectors of the surface of an imaginary sphere. They are distributed with the same nearest neighbor distance, as is shown in Figure 1. The polarization of the electric field is then chosen randomly for each wave.
To check the convergence of our method we increase the number of illuminations used to calculate the Tmatrix. We introduce
as a measure for the deviation of the Tmatrix. Here, (i) represents an index that runs over all matrix components and T_{Mie} the analytically known Tmatrix. This quantity is visualized in Figure 2 in a double logarithmic plot.
We clearly see that the coefficients can be reproduced with sufficient precision. The deviation drops at first if we increase the number of illuminations, then a plateau is reached. The remaining deviation can be explained by numerical errors introduced by the finite meshing of the sphere in the Maxwell solver as can be seen by considering a second calculation with a finer mesh.
Dimer
We now consider a slightly more complex object to further investigate the convergence of our method. The dimer, i.e., two coupled spheres, is a suitable example, because we can calculate the Tmatrix semianalytically as a reference. We get the Tmatrix of the single spheres directly from the Mie coefficients and construct the total Tmatrix by translating the coefficients with addition theorems [41]. Please note that this solution is not completely accurate, because we truncate the infinite sums at N. Additionally, we calculate the Tmatrix from the illumination with plane waves with a finite element solver, according to Equaiton Equation 15 and compare the results.
The object we investigate consists of two strongly coupled silver spheres with radius r = 30 nm and a centertocenter distance of d = 63 nm embedded in glass with ε = 2.25 [42]. We take established experimental data for the dispersive permittivity of silver [43]. Such objects can be fabricated in large quantities by self assembly methods, e.g., by connecting commercially available metal nanospheres with a linker molecule [44].
We set N = 2, because the higher orders do not contribute notably. For general objects we can construct a condition to find the optimal multipole order. We perform the Tmatrix calculation with increasing multipolar order and check the contribution of each order by summing the corresponding Tmatrix entries
If we reach an order that has only a negligible contribution compared to the dominant ones, we stop the procedure.
As depicted in Figure 3, the scattering cross sections obtained from the different Tmatrix methods agree very well. We can clearly identify three different resonances: a broad electric dipole resonance at 680 THz, and two resonances with magnetic dipole and electric quadrupole contributions at 620 and 710 THz, respectively. While the broad resonance is connected to the eigenmode of the single sphere, the two sharper peaks can be explained by a hybridization caused by the coupling of the two spheres and coupling of higher order multipole modes, respectively [42]. Here the multipole contributions are calculated according to Equation 8.
We get additional insights by looking at the entries of the Tmatrix in Figure 4 at a discrete frequency. The way the Tmatrix is presented in that figure will be the same in all the following figures showing the Tmatrix. Note that the entries are complex numbers and we always display the absolute values. At first we notice that the dominant entries are connected to the electric dipole and, to a lesser extent, to the magnetic dipole and electric quadrupole components . Additionally, we see the aforementioned cross coupling between magnetic dipole and electric quadrupole in the submatrix and all permutations. This cross coupling also became evident in the scattering cross section.
Now, we investigate the behavior when we increase the number of illuminations K in the Tmatrix calculation. We consider the Tmatrix with N = 4 at the frequency 600 THz, where several multipole orders contribute to the total scattering. As a reference, we take the calculation based on Mie coefficients. In Figure 5 we can observe the expected behavior of the error, the deviation goes down exponentially as we increase the number of plane wave illuminations. So the more complicated structure actually benefits from multiple additional illuminations. In Figure 6 we see the Tmatrix entries of two different numbers of illuminations. We clearly see that the numerical noise, i.e., the parts of the matrix that should be zero, is going down significantly.
In conclusion, we can state that the presented method delivers correct results and that it is highly beneficial to use a few illuminations more than absolutely necessary.
Sandwich particle
Another interesting and in a sense similar object is the sandwich particle. It consists of two metallic disks, usually separated by a dielectric spacer layer. This particle is preferably fabricated by topdown procedures, such as electron beam lithography [5,45].
The investigated object has a disk radius of 60 nm, a disk height of 30 nm and a gap thickness of 10 nm. The metal disks consist of gold [43] and the spacer is made of SiO_{2} where we assumed a constant permittivity of ε = 2.13 and the surrounding is vacuum. Similarly to the dimer case, the two separated disks support a magnetic dipole resonance at 336 THz, as can be seen in Figure 7. Additionally, there is a strong and broad electric dipole resonance at 520 THz that is, again, attributed to the resonance of a single disk.
We can observe a magnetic resonance that is weaker than in the case of the dimer, but the big advantage is that the peak is well isolated from other resonances. Furthermore, we observe no higher multipoleorder contributions at all in the investigated frequency region. This is particularly beneficial for homogenization, because we can apply the Clausius–Mosotti relation, which is only valid in dipole approximation. The different geometry makes it now possible to shift the electric and magnetic dipole resonances almost independently, because we have more parameters. The electric mode depends strongly on the radius of the disks, while the magnetic mode can be shifted by changing the thickness of the gold and spacer layers, respectively.
Again, we can get a deeper understanding of the response, by looking at the Tmatrix. Figure 8 shows the entries of the matrix of the presented structure at the position of the magnetic resonance, 336 THz. We notice a strong similarity to the previous case of the dimer. This is because the object has almost exactly the same geometrical symmetries. However, we also see some notable differences. Firstly, the particle can be described in very good approximation as a dipole, because the parts corresponding to higher multipoles are negligible. Furthermore, the cross coupling between the magnetic dipole and the electric quadrupole is significantly lower than in the previous case. This is the cause of the well isolated magnetic resonance we saw in the scattering cross section in Figure 7 in contrast to the mixed magnetic dipole and electric quadrupole resonances of the dimer.
Additionally, we can deduce the preferential orientation of the induced multipoles from the pattern of the Tmatrix. Let us investigate the electric dipolar part and call the entries
Now connects the incident coefficient p_{1j} to the scattering coefficient a_{1i}. The same is true for all the other parts of the Tmatrix, with adjustments for the higherorder parts. The entry is special because it is connected to the zaxis. The zaxis is distinguished in the definition of the vector spherical harmonic functions. This means that if this entry has a large absolute value, the corresponding object will exhibit a strong electric dipole in zdirection, . Conversely the dipole moments in the x–yplane are proportional to linear combinations of entries with ±1 indices. This also translates directly to the magnetic part of the Tmatrix. Considering this, we investigate again Figure 8. We see that the electric dipole is almost isotropic. This becomes clear when we compare the Tmatrix to that of a sphere, as stated above, we would get = = in this completely isotropic case. The magnetic dipole, however, is much stronger in the x–y plane because of the geometry of the investigated object with the gap. The magnetic dipole moment is induced by antiparallel currents in the metal disks. Please note that these considerations are completely independent of the illumination, because we investigate the Tmatrix directly.
Helix
Now we can utilize the biggest advantage of the proposed method, namely that we can calculate the Tmatrix of an arbitrarily complex object. To showcase this we compute the Tmatrix of a metal helix. In this case we can not make use of inherent symmetries of the object a priori and have no analytic solution available.
Calculating the Tmatrix is particularly interesting for the present object, because numerous electromagnetic properties can be deduced from its entries, such as the duality conservation and electromagnetic chirality [46]. The helix was optimized to express a strong electromagnetic chirality at a specific frequency and, thus, shows a strong contrast in the scattering cross section for two opposite circular polarized incident fields. The major radius of the helix is 6.48 μm, the pitch is 8.52 μm, the minor radius is 0.8 μm, and the helix is surrounded by vacuum. The operating frequency of the system, i.e., the position of the resonance is at 1.5 THz. As we see in Figure 9, there is a strong difference in the scattering cross section if we illuminate with righthanded or lefthanded circularly polarized light.
For the following considerations it is beneficial to introduce shortly the concept of duality symmetry. In free space, Maxwell’s equations are invariant under the transformation [47]
where E is the electric field, H is the magnetic field, Z is the impedance and θ is an arbitrary constant angle. The symmetry that is corresponding to this transformation invariance is duality. Helicity is defined as the product of the total angular momentum J and the direction of the linear momentum of the wave P/P[48]
It has the two eigenvalues 1 and −1, and the corresponding eigenstates [49]
For example, a circularly polarized plane wave is a state of pure helicity and it is preserved if duality symmetry applies.
We can see directly that the object is approximately dual at the resonance frequency of 1.5 THz, because the magnetic and electric dipole contributions are balanced and no higher multipole contributions arise [50]. This means that the interaction of the particle with the incident light does not change the helicity of the incident waves.
Looking at the Tmatrix entries in Figure 10 we see, like in the previous cases, that the object can be described mainly by the dipole contributions. The interesting aspect here is that the electric and magnetic dipoles have almost the same strength and are both aligned along the zaxis at the same frequency. Furthermore the mixing of electric and magnetic dipole moments is also very strong. This means that an electric dipole is induced by a magnetic dipole excitation and vice versa. All these properties lead to the approximate duality symmetry of the object.
As stated above, we can calculate measures for the duality breaking, electromagnetic chirality and the scattering cross section contrast for general plane wave illuminations from the Tmatrix entries. This is done by transforming the Tmatrix from its current form to the helicity basis
The submatrices represent now the coupling of the helicity + and − states. In Figure 11 the entries of the transformed matrix are depicted for a lefthanded and a righthanded helix. We can clearly distinguish the two cases, because each helix interacts mainly only with one helicity. The left handed helix interacts with helicity − () states and the right handed interacts with helicity + ().
To further analyze the Tmatrix of the object we can perform a singular value decomposition (svd) of the submatrices [46]. The values are ordered to get the vectors
and
from which we calculate the quantities that are related to helicity
Here, is the duality breaking, C can be interpreted as the averaged total scattering cross section independent of a specific illumination, Δ is the cross section contrast for different circular polarizations, and χ is a measure for the electromagnetic chirality [46]. The values for the chosen helix are shown in Figure 12.
We clearly see the anticipated minimum in the duality breaking. This means that the helix is approximately dual at the resonance frequency. Simultaneously the chirality reaches a maximum of 0.9, which means that the object responds differently to the two opposite polarizations of illuminating radiation. This is also evident from the strong cross section contrast.
The presented quantities, and possibly others, are only accessible because we calculated the Tmatrix of the object. This information makes it also possible to optimize such objects numerically to achieve the desired behavior.
Conclusion
We demonstrated an algorithm to calculate the Tmatrix of an arbitrary object. We showcased numerical examples by calculating the Tmatrix of four objects. Especially the calculation of the silver helix demonstrates the strength of the approach. We investigated properties that are usually not available with established full wave solvers.
The Tmatrix is calculated with the help of widely available and well established tools to calculate the scattered fields upon plane wave illumination. The Tmatrix contains the electromagnetic scattering information independently of the illumination. By investigating the matrix, we get valuable information about the multipolar composition of the scattering cross section and electromagnetic chirality. This is extremely useful for understanding the nature of a given response.
Furthermore, the Tmatrix can be used in an extended multiscattering formalism to calculate the response of a whole cluster of objects efficiently. This will speed up calculations by a considerable amount.
Acknowledgements
This work was supported by the German Science Foundation within project RO 3640/31. We also acknowledge partial financial support by the Deutsche Forschungsgemeinschaft through CRC 1173. We would like to thank JCMwave for providing us access to the FEM solver JCMsuite. M.F. also acknowledges the support by the Karlsruhe School of Optics & Photonics (KSOP).
References

Loudet, J.C.; Barois, P.; Poulin, P. Nature 2000, 407, 611–613. doi:10.1038/35036539
Return to citation in text: [1] 
Gantzounis, G.; Stefanou, N.; Yannopapas, V. J. Phys.: Condens. Matter 2005, 17, 1791. doi:10.1088/09538984/17/12/003
Return to citation in text: [1] 
Frein, S.; Boudon, J.; Vonlanthen, M.; Scharf, T.; Barberá, J.; SüssFink, G.; Bürgi, T.; Deschenaux, R. Helv. Chim. Acta 2008, 91, 2321–2337. doi:10.1002/hlca.200890253
Return to citation in text: [1] 
Wojcik, M.; Lewandowski, W.; Matraszek, J.; Mieczkowski, J.; Borysiuk, J.; Pociecha, D.; Gorecka, E. Angew. Chem., Int. Ed. 2009, 48, 5167–5169. doi:10.1002/anie.200901206
Return to citation in text: [1] 
Lehr, D.; Dietrich, K.; Helgert, C.; Käsebier, T.; Fuchs, H.J.; Tünnermann, A.; Kley, E.B. Opt. Lett. 2012, 37, 157–159. doi:10.1364/OL.37.000157
Return to citation in text: [1] [2] 
Asadchy, V.; Albooyeh, M.; Tretyakov, S. J. Opt. Soc. Am. B 2016, 33, A16–A20. doi:10.1364/JOSAB.33.000A16
Return to citation in text: [1] 
Alaee, R.; Albooyeh, M.; Rahimzadegan, A.; Mirmoosa, M. S.; Kivshar, Y. S.; Rockstuhl, C. Phys. Rev. B 2015, 92, 245130. doi:10.1103/PhysRevB.92.245130
Return to citation in text: [1] 
Li, L. J. Opt. Soc. Am. A 1997, 14, 2758–2767. doi:10.1364/JOSAA.14.002758
Return to citation in text: [1] 
Modinos, A.; Stefanou, N.; Yannopapas, V. Opt. Express 2001, 8, 197–202. doi:10.1364/OE.8.000197
Return to citation in text: [1] 
Oskooi, A. F.; Roundy, D.; Ibanescu, M.; Bermel, P.; Joannopoulos, J. D.; Johnson, S. G. Comput. Phys. Commun. 2010, 181, 687–702. doi:10.1016/j.cpc.2009.11.008
Return to citation in text: [1] 
Ippolito, S. J.; KalantarZadeh, K.; Powell, D. A.; Wlodarski, W. A. 3dimensional finite element approach for simulating acoustic wave propagation in layered SAW devices. In Ultrasonics, 2003 IEEE Symposium on, IEEE Publishing: Piscataway, NJ, U.S.A., 2003; pp 303–306.
Return to citation in text: [1] 
Pomplun, J.; Burger, S.; Zschiedrich, L.; Schmidt, F. Phys. Status Solidi B 2007, 244, 3419–3434. doi:10.1002/pssb.200743192
Return to citation in text: [1] [2] 
Moreau, A.; Ciracì, C.; Mock, J. J.; Hill, R. T.; Wang, Q.; Wiley, B. J.; Chilkoti, A.; Smith, D. R. Nature 2012, 492, 86–89. doi:10.1038/nature11615
Return to citation in text: [1] 
Yang, P.; Liou, K. N.; Mishchenko, M. I.; Gao, B.C. Appl. Opt. 2000, 39, 3727–3737. doi:10.1364/AO.39.003727
Return to citation in text: [1] 
Kuwata, H.; Tamaru, H.; Esumi, K.; Miyano, K. Appl. Phys. Lett. 2003, 83, 4625–4627. doi:10.1063/1.1630351
Return to citation in text: [1] 
Nordlander, P.; Oubre, C.; Prodan, E.; Li, K.; Stockman, M. I. Nano Lett. 2004, 4, 899–903. doi:10.1021/nl049681c
Return to citation in text: [1] 
Mackowski, D. W. J. Opt. Soc. Am. A 2002, 19, 881–893. doi:10.1364/JOSAA.19.000881
Return to citation in text: [1] 
Draine, B. T. The discrete dipole approximation for light scattering by irregular targets. In Light Scattering by Nonspherical Particles: Theory, Measurements, and Applications; Mishchenko, M. I.; Hovenier, J. W.; Travis, L. D., Eds.; Academic Press: Cambridge, MA, U.S.A., 2000; p 131.
Return to citation in text: [1] 
Xu, Y.l. Appl. Opt. 1995, 34, 4573–4588. doi:10.1364/AO.34.004573
Return to citation in text: [1] 
Mühlig, S.; Rockstuhl, C.; Pniewski, J.; Simovski, C. R.; Tretyakov, S. A.; Lederer, F. Phys. Rev. B 2010, 81, 075317. doi:10.1103/PhysRevB.81.075317
Return to citation in text: [1] [2] 
Waterman, P. C. Proc. IEEE 1965, 53, 805–812. doi:10.1109/PROC.1965.4058
Return to citation in text: [1] [2] [3] 
Gimbutas, Z.; Greengard, L. J. Comput. Phys. 2013, 232, 22–32. doi:10.1016/j.jcp.2012.01.041
Return to citation in text: [1] [2] 
Fruhnert, M.; Monti, A.; FernandezCorbaton, I.; Alù, A.; Toscano, A.; Bilotti, F.; Rockstuhl, C. Phys. Rev. B 2016, 93, 245127. doi:10.1103/PhysRevB.93.245127
Return to citation in text: [1] [2] 
Auger, J.C.; Stout, B. J. Quant. Spectrosc. Radiat. Transfer 2003, 79–80, 533–547. doi:10.1016/S00224073(02)003060
Return to citation in text: [1] 
Mackowski, D. W.; Mishchenko, M. I. J. Opt. Soc. Am. A 1996, 13, 2266–2278. doi:10.1364/JOSAA.13.002266
Return to citation in text: [1] 
Stefanou, N.; Yannopapas, V.; Modinos, A. Comput. Phys. Commun. 2000, 132, 189–196. doi:10.1016/S00104655(00)001314
Return to citation in text: [1] 
Yannopapas, V.; Vanakaras, A. G. Phys. Rev. B 2011, 84, 085119. doi:10.1103/PhysRevB.84.085119
Return to citation in text: [1] 
Bohren, C. F.; Huffman, D. R. Absorption and scattering of light by small particles; John Wiley & Sons: New York, NY, U.S.A., 2008.
Return to citation in text: [1] 
Mie, G. Ann. Phys. 1908, 330, 377–445. doi:10.1002/andp.19083300302
Return to citation in text: [1] 
Mischenko, M. I.; Travis, L. D.; Lacis, A. A. Scattering, absorption, and emission of light by small particles; Cambridge University Press: Cambridge, United Kingdom, 2002.
Return to citation in text: [1] 
Mishchenko, M. I.; Travis, L. D.; Mackowski, D. W. J. Quant. Spectrosc. Radiat. Transfer 1996, 55, 535–575. doi:10.1016/00224073(96)000027
Return to citation in text: [1] 
Tsang, L.; Kong, J.; Shin, R. Theory of microwave remote sensing; WileyInterscience: New York, NY, U.S.A., 1985.
Return to citation in text: [1] 
Mishchenko, M. I.; Dlugach, J. M. Planet. Space Sci. 1993, 41, 173–181. doi:10.1016/00320633(93)900568
Return to citation in text: [1] 
Nieminen, T. A.; RubinszteinDunlop, H.; Heckenberg, N. R. J. Quant. Spectrosc. Radiat. Transfer 2003, 79, 1019–1029. doi:10.1016/S00224073(02)003369
Return to citation in text: [1] 
Yannopapas, V. J. Opt. Soc. Am. B 2014, 31, 631–636. doi:10.1364/JOSAB.31.000631
Return to citation in text: [1] 
Sersic, I.; Tuambilangana, C.; Kampfrath, T.; Koenderink, A. F. Phys. Rev. B 2011, 83, 245102. doi:10.1103/PhysRevB.83.245102
Return to citation in text: [1] 
Asadchy, V. S.; Faniayeu, I. A.; Ra’di, Y.; Tretyakov, S. A. Photonics Nanostruct.  Fundam. Appl. 2014, 12, 298–304. doi:10.1016/j.photonics.2014.04.004
Return to citation in text: [1] 
Liu, X.X.; Zhao, Y.; Alù, A. IEEE Trans. Antennas Propag. 2016, 64, 2301–2310. doi:10.1109/TAP.2016.2546958
Return to citation in text: [1] 
Yazdi, M.; Komjani, N. J. Opt. Soc. Am. B 2016, 33, 491–500. doi:10.1364/JOSAB.33.000491
Return to citation in text: [1] 
Trefethen, L. N.; Bau, D., III. Numerical linear algebra; Siam, 1997.
Return to citation in text: [1] 
Stein, S. Q. Appl. Math. 1961, 15–24. doi:10.1090/qam/120407
Return to citation in text: [1] 
Mühlig, S.; Menzel, C.; Rockstuhl, C.; Lederer, F. Metamaterials 2011, 5, 64–73. doi:10.1016/j.metmat.2011.03.003
Return to citation in text: [1] [2] 
Johnson, P. B.; Christy, R. W. Phys. Rev. B 1972, 6, 4370–4379. doi:10.1103/PhysRevB.6.4370
Return to citation in text: [1] [2] 
Fruhnert, M.; Kretschmer, F.; Geiss, R.; Perevyazko, I.; CiallaMay, D.; Steinert, M.; Janunts, N.; Sivun, D.; Hoeppener, S.; Hager, M. D.; Pertsch, T.; Schubert, U. S.; Rockstuhl, C. J. Phys. Chem. C 2015, 119, 17809–17817. doi:10.1021/acs.jpcc.5b04346
Return to citation in text: [1] 
Hicks, E. M.; Zou, S.; Schatz, G. C.; Spears, K. G.; Van Duyne, R. P.; Gunnarsson, L.; Rindzevicius, T.; Kasemo, B.; Käll, M. Nano Lett. 2005, 5, 1065–1070. doi:10.1021/nl0505492
Return to citation in text: [1] 
FernandezCorbaton, I.; Fruhnert, M.; Rockstuhl, C. Phys. Rev. X 2016, 6, 031013. doi:10.1103/PhysRevX.6.031013
Return to citation in text: [1] [2] [3] 
Calkin, M. G. Am. J. Phys. 1965, 33, 958–960. doi:10.1119/1.1971089
Return to citation in text: [1] 
Tung, W.K. Group theory in physics; World Scientific Publishing Co Pte Ltd: Singapore, 1985.
Return to citation in text: [1] 
Lipkin, D. M. J. Math. Phys. 1964, 5, 696–700. doi:10.1063/1.1704165
Return to citation in text: [1] 
FernandezCorbaton, I.; Fruhnert, M.; Rockstuhl, C. ACS Photonics 2015, 2, 376–384. doi:10.1021/ph500419a
Return to citation in text: [1]
48.  Tung, W.K. Group theory in physics; World Scientific Publishing Co Pte Ltd: Singapore, 1985. 
1.  Loudet, J.C.; Barois, P.; Poulin, P. Nature 2000, 407, 611–613. doi:10.1038/35036539 
2.  Gantzounis, G.; Stefanou, N.; Yannopapas, V. J. Phys.: Condens. Matter 2005, 17, 1791. doi:10.1088/09538984/17/12/003 
3.  Frein, S.; Boudon, J.; Vonlanthen, M.; Scharf, T.; Barberá, J.; SüssFink, G.; Bürgi, T.; Deschenaux, R. Helv. Chim. Acta 2008, 91, 2321–2337. doi:10.1002/hlca.200890253 
4.  Wojcik, M.; Lewandowski, W.; Matraszek, J.; Mieczkowski, J.; Borysiuk, J.; Pociecha, D.; Gorecka, E. Angew. Chem., Int. Ed. 2009, 48, 5167–5169. doi:10.1002/anie.200901206 
5.  Lehr, D.; Dietrich, K.; Helgert, C.; Käsebier, T.; Fuchs, H.J.; Tünnermann, A.; Kley, E.B. Opt. Lett. 2012, 37, 157–159. doi:10.1364/OL.37.000157 
6.  Asadchy, V.; Albooyeh, M.; Tretyakov, S. J. Opt. Soc. Am. B 2016, 33, A16–A20. doi:10.1364/JOSAB.33.000A16 
14.  Yang, P.; Liou, K. N.; Mishchenko, M. I.; Gao, B.C. Appl. Opt. 2000, 39, 3727–3737. doi:10.1364/AO.39.003727 
15.  Kuwata, H.; Tamaru, H.; Esumi, K.; Miyano, K. Appl. Phys. Lett. 2003, 83, 4625–4627. doi:10.1063/1.1630351 
16.  Nordlander, P.; Oubre, C.; Prodan, E.; Li, K.; Stockman, M. I. Nano Lett. 2004, 4, 899–903. doi:10.1021/nl049681c 
34.  Nieminen, T. A.; RubinszteinDunlop, H.; Heckenberg, N. R. J. Quant. Spectrosc. Radiat. Transfer 2003, 79, 1019–1029. doi:10.1016/S00224073(02)003369 
35.  Yannopapas, V. J. Opt. Soc. Am. B 2014, 31, 631–636. doi:10.1364/JOSAB.31.000631 
11.  Ippolito, S. J.; KalantarZadeh, K.; Powell, D. A.; Wlodarski, W. A. 3dimensional finite element approach for simulating acoustic wave propagation in layered SAW devices. In Ultrasonics, 2003 IEEE Symposium on, IEEE Publishing: Piscataway, NJ, U.S.A., 2003; pp 303–306. 
12.  Pomplun, J.; Burger, S.; Zschiedrich, L.; Schmidt, F. Phys. Status Solidi B 2007, 244, 3419–3434. doi:10.1002/pssb.200743192 
13.  Moreau, A.; Ciracì, C.; Mock, J. J.; Hill, R. T.; Wang, Q.; Wiley, B. J.; Chilkoti, A.; Smith, D. R. Nature 2012, 492, 86–89. doi:10.1038/nature11615 
22.  Gimbutas, Z.; Greengard, L. J. Comput. Phys. 2013, 232, 22–32. doi:10.1016/j.jcp.2012.01.041 
8.  Li, L. J. Opt. Soc. Am. A 1997, 14, 2758–2767. doi:10.1364/JOSAA.14.002758 
9.  Modinos, A.; Stefanou, N.; Yannopapas, V. Opt. Express 2001, 8, 197–202. doi:10.1364/OE.8.000197 
10.  Oskooi, A. F.; Roundy, D.; Ibanescu, M.; Bermel, P.; Joannopoulos, J. D.; Johnson, S. G. Comput. Phys. Commun. 2010, 181, 687–702. doi:10.1016/j.cpc.2009.11.008 
23.  Fruhnert, M.; Monti, A.; FernandezCorbaton, I.; Alù, A.; Toscano, A.; Bilotti, F.; Rockstuhl, C. Phys. Rev. B 2016, 93, 245127. doi:10.1103/PhysRevB.93.245127 
31.  Mishchenko, M. I.; Travis, L. D.; Mackowski, D. W. J. Quant. Spectrosc. Radiat. Transfer 1996, 55, 535–575. doi:10.1016/00224073(96)000027 
7.  Alaee, R.; Albooyeh, M.; Rahimzadegan, A.; Mirmoosa, M. S.; Kivshar, Y. S.; Rockstuhl, C. Phys. Rev. B 2015, 92, 245130. doi:10.1103/PhysRevB.92.245130 
21.  Waterman, P. C. Proc. IEEE 1965, 53, 805–812. doi:10.1109/PROC.1965.4058 
32.  Tsang, L.; Kong, J.; Shin, R. Theory of microwave remote sensing; WileyInterscience: New York, NY, U.S.A., 1985. 
33.  Mishchenko, M. I.; Dlugach, J. M. Planet. Space Sci. 1993, 41, 173–181. doi:10.1016/00320633(93)900568 
20.  Mühlig, S.; Rockstuhl, C.; Pniewski, J.; Simovski, C. R.; Tretyakov, S. A.; Lederer, F. Phys. Rev. B 2010, 81, 075317. doi:10.1103/PhysRevB.81.075317 
24.  Auger, J.C.; Stout, B. J. Quant. Spectrosc. Radiat. Transfer 2003, 79–80, 533–547. doi:10.1016/S00224073(02)003060 
28.  Bohren, C. F.; Huffman, D. R. Absorption and scattering of light by small particles; John Wiley & Sons: New York, NY, U.S.A., 2008. 
46.  FernandezCorbaton, I.; Fruhnert, M.; Rockstuhl, C. Phys. Rev. X 2016, 6, 031013. doi:10.1103/PhysRevX.6.031013 
21.  Waterman, P. C. Proc. IEEE 1965, 53, 805–812. doi:10.1109/PROC.1965.4058 
22.  Gimbutas, Z.; Greengard, L. J. Comput. Phys. 2013, 232, 22–32. doi:10.1016/j.jcp.2012.01.041 
23.  Fruhnert, M.; Monti, A.; FernandezCorbaton, I.; Alù, A.; Toscano, A.; Bilotti, F.; Rockstuhl, C. Phys. Rev. B 2016, 93, 245127. doi:10.1103/PhysRevB.93.245127 
29.  Mie, G. Ann. Phys. 1908, 330, 377–445. doi:10.1002/andp.19083300302 
30.  Mischenko, M. I.; Travis, L. D.; Lacis, A. A. Scattering, absorption, and emission of light by small particles; Cambridge University Press: Cambridge, United Kingdom, 2002. 
19.  Xu, Y.l. Appl. Opt. 1995, 34, 4573–4588. doi:10.1364/AO.34.004573 
20.  Mühlig, S.; Rockstuhl, C.; Pniewski, J.; Simovski, C. R.; Tretyakov, S. A.; Lederer, F. Phys. Rev. B 2010, 81, 075317. doi:10.1103/PhysRevB.81.075317 
50.  FernandezCorbaton, I.; Fruhnert, M.; Rockstuhl, C. ACS Photonics 2015, 2, 376–384. doi:10.1021/ph500419a 
17.  Mackowski, D. W. J. Opt. Soc. Am. A 2002, 19, 881–893. doi:10.1364/JOSAA.19.000881 
18.  Draine, B. T. The discrete dipole approximation for light scattering by irregular targets. In Light Scattering by Nonspherical Particles: Theory, Measurements, and Applications; Mishchenko, M. I.; Hovenier, J. W.; Travis, L. D., Eds.; Academic Press: Cambridge, MA, U.S.A., 2000; p 131. 
21.  Waterman, P. C. Proc. IEEE 1965, 53, 805–812. doi:10.1109/PROC.1965.4058 
25.  Mackowski, D. W.; Mishchenko, M. I. J. Opt. Soc. Am. A 1996, 13, 2266–2278. doi:10.1364/JOSAA.13.002266 
26.  Stefanou, N.; Yannopapas, V.; Modinos, A. Comput. Phys. Commun. 2000, 132, 189–196. doi:10.1016/S00104655(00)001314 
27.  Yannopapas, V.; Vanakaras, A. G. Phys. Rev. B 2011, 84, 085119. doi:10.1103/PhysRevB.84.085119 
46.  FernandezCorbaton, I.; Fruhnert, M.; Rockstuhl, C. Phys. Rev. X 2016, 6, 031013. doi:10.1103/PhysRevX.6.031013 
36.  Sersic, I.; Tuambilangana, C.; Kampfrath, T.; Koenderink, A. F. Phys. Rev. B 2011, 83, 245102. doi:10.1103/PhysRevB.83.245102 
37.  Asadchy, V. S.; Faniayeu, I. A.; Ra’di, Y.; Tretyakov, S. A. Photonics Nanostruct.  Fundam. Appl. 2014, 12, 298–304. doi:10.1016/j.photonics.2014.04.004 
38.  Liu, X.X.; Zhao, Y.; Alù, A. IEEE Trans. Antennas Propag. 2016, 64, 2301–2310. doi:10.1109/TAP.2016.2546958 
39.  Yazdi, M.; Komjani, N. J. Opt. Soc. Am. B 2016, 33, 491–500. doi:10.1364/JOSAB.33.000491 
12.  Pomplun, J.; Burger, S.; Zschiedrich, L.; Schmidt, F. Phys. Status Solidi B 2007, 244, 3419–3434. doi:10.1002/pssb.200743192 
43.  Johnson, P. B.; Christy, R. W. Phys. Rev. B 1972, 6, 4370–4379. doi:10.1103/PhysRevB.6.4370 
46.  FernandezCorbaton, I.; Fruhnert, M.; Rockstuhl, C. Phys. Rev. X 2016, 6, 031013. doi:10.1103/PhysRevX.6.031013 
42.  Mühlig, S.; Menzel, C.; Rockstuhl, C.; Lederer, F. Metamaterials 2011, 5, 64–73. doi:10.1016/j.metmat.2011.03.003 
5.  Lehr, D.; Dietrich, K.; Helgert, C.; Käsebier, T.; Fuchs, H.J.; Tünnermann, A.; Kley, E.B. Opt. Lett. 2012, 37, 157–159. doi:10.1364/OL.37.000157 
45.  Hicks, E. M.; Zou, S.; Schatz, G. C.; Spears, K. G.; Van Duyne, R. P.; Gunnarsson, L.; Rindzevicius, T.; Kasemo, B.; Käll, M. Nano Lett. 2005, 5, 1065–1070. doi:10.1021/nl0505492 
43.  Johnson, P. B.; Christy, R. W. Phys. Rev. B 1972, 6, 4370–4379. doi:10.1103/PhysRevB.6.4370 
44.  Fruhnert, M.; Kretschmer, F.; Geiss, R.; Perevyazko, I.; CiallaMay, D.; Steinert, M.; Janunts, N.; Sivun, D.; Hoeppener, S.; Hager, M. D.; Pertsch, T.; Schubert, U. S.; Rockstuhl, C. J. Phys. Chem. C 2015, 119, 17809–17817. doi:10.1021/acs.jpcc.5b04346 
42.  Mühlig, S.; Menzel, C.; Rockstuhl, C.; Lederer, F. Metamaterials 2011, 5, 64–73. doi:10.1016/j.metmat.2011.03.003 
© 2017 Fruhnert et al.; licensee BeilsteinInstitut.
This is an Open Access article under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
The license is subject to the Beilstein Journal of Nanotechnology terms and conditions: (http://www.beilsteinjournals.org/bjnano)