3-dimensional eigenmodal analysis of plasmonic nanostructures

doi:10.1364/OE.20.005481

« Show journal navigation

3-dimensional eigenmodal analysis of plasmonic nanostructures

Hua Guo, Benedikt Oswald, and Peter Arbenz »View Author Affiliations

Optics Express, Vol. 20, Issue 5, pp. 5481-5500 (2012)
http://dx.doi.org/10.1364/OE.20.005481

View Full Text Article

Acrobat PDF (1179 KB)

Journal Search
Article Lookup

Article Tools

Citations

Alert me when this article is cited
Export Citation/Save

Article Navigation

▸ Article Outline
– Abstract
1. Introduction
2. Formulation of the problem
3. Numerical solution
4. Validation and application of the algorithm
5. Discussion and conclusions
– References and links

Article Tools
Full Text
Article Info
References
Cited By
Figures
Metrics
Download PDF

Viewing Equations in the
HTML Article

Optimal Viewing Recommendations

Full Text
Article Info
References (50)
Cited By
Figures (9)
Metrics

Abstract

We introduce a 3-dimensional electromagnetic eigenmodal algorithm for the theoretical analysis of resonating nano-optical structures. The method, a variant of the Jacobi–Davidson algorithm, solves the electric field vector wave, or curl-curl, equation for the electromagnetic eigenmodes of resonant optical structures with a finite element method. In particular, the method includes transparent boundary conditions that enable the analysis of resonating structures in unbounded space. We demonstrate the performance of the method. First, we calculate the modes of several dielectric resonator antennas and compare them to theoretically determined results. Second, we calculate the modes of a nano-cuboid and compare them to theoretically determined results. Third, we numerically analyze spherical nanoparticles and compare the result to the theoretical Mie solution. Fourth, we analyze optical dipole antenna configurations in order to assess the method’s capability for solving technologically relevant problems.

1. Introduction

The interaction of nano-meter structured, metallic, aka. plasmonic, devices with light, especially from lasers, has attracted considerable interest during the last few years. Consequentially, a significant body of literature exists, see e.g. [1

L. Novotny and N. van Hulst, “Antennas for light,” Nat. Photon. 5, 83–90 (2011). [CrossRef]

–4

W. L. Barnes, A. Dereux, and T. W. Ebbesen, “Surface plasmon subwavelength optics,” Nature 424, 824–830 (2003). [CrossRef] [PubMed]

]. The concept of the antenna, in particular, transferred from the microwave to the optical region, has attracted enormous interest since the conversion of propagating light into localized, electromagnetic energy, and vice-versa, holds the promise of novel technological applications in many fields, such as in sensing, detection manipulation and communication [5

H. Aouani, O. Mahboub, N. Bonod, E. Devaux, E. Popov, H. Rigneault, T. W. Ebbesen, and J. Wenger, “Bright unidirectional fluorescence emission of molecules in a nanoaperture with plasmonic corrugations,” Nano Lett. 11, 637–644 (2011). [CrossRef] [PubMed]

–8

M. F. Garcia-Parajo, “Optical antennas focus in on biology,” Nat. Photon. 2, 201–203, (2008). [CrossRef]

]. Certainly, it is well beyond the scope of this study to give a comprehensive overview. Simultaneously, fabrication technology has advanced impressively and is now capable of producing ultra-smooth structures with justifiable effort [9

P. Nagpal, N. C. Lindquist, S.-H. Oh, and D. J. Norris, “Ultrasmooth patterned metals for plasmonics and metamaterials,” Science 325, 594–597 (2009). [CrossRef] [PubMed]

Since the behavior of nano-optical devices is sometimes counter-intuitive and the fabrication processes are often time-consuming, the theoretical analysis of these devices has become indispensable in order to identify promising candidates and to assess their performance before resources are committed in the clean room.

Generally speaking, theoretical methods are classified into analytical [10

G. W. Hanson and A. B. Yakovlev, Operator Theory for Electromagnetics – An Introduction (Springer, 2002). [PubMed]

], semi-analytical [11

Ch. Hafner, The Generalized Multipole Technique for Computational Electromagnetics (Artech House Books, Boston, 1990).

] and numerical approaches [12

J. Jin, The Finite Element Method in Electromagnetics (John Wiley, New York, 2002).

–14

J. L. Volakis, A. Chatterjee, and L. C. Kempel, Finite Element Method for Electromagnetics – Antennas, Microwave Circuits and Scattering Applications (IEEE Press, New York, 1998).

]. While analytical methods can provide detailed insight into the behavior of devices and allow for relatively simple parameter sweeps, quite often their applicability is limited to simplified geometries and, thus, their predictions are of limited validity when technologically relevant problems are studied. Numerical methods on the other hand allow for the detailed modeling of delicate geometrical details, almost fully preserving the reality of the device. Thus, they are the methods of choice for detailed performance analysis. We note however that such detailed wealth of prediction comes at a price in the form of significantly increased computational effort and substantial pre- and post-processing work. In between these methods is the class of semi-analytical methods that are usually more flexible with respect to geometry, but on the other hand still provide analytical insight. Among the many semi-analytical methods we particularly mention the multiple multipole (MMP) method, originated by Hafner [11

Ch. Hafner, The Generalized Multipole Technique for Computational Electromagnetics (Artech House Books, Boston, 1990).

] in the 1980s. The MMP method subdivides geometry into subdomains for which the analytical solutions are available. It then matches the superposition of these scaled analytical solutions on the interfaces between the subdomains by formulating the mismatch as an error energy. The resulting, overdetermined linear system is then solved and the mismatch on the subdomain interfaces is a direct measure for accuracy.

Here, we use the finite element method (FEM) [12

J. Jin, The Finite Element Method in Electromagnetics (John Wiley, New York, 2002).

] in the time-harmonic regime to calculate electromagnetic modes associated with resonant nano-optical devices. We derive and solve the corresponding electromagnetic eigenmodal problem. The resonance wavelength, the decay rate and the quality factor (Q factor) are derived from the (complex) eigenvalue, and the electro-magnetic fields are derived from the calculated eigenvector. The radiated power, dissipated energy and stored energy of the system can be computed as well. We demonstrate the performance of the method by solving four different problems, namely:

The dielectric resonator antenna (DRA) [15
A. Okaya and L. F. Barash, “The dielectric microwave resonator,” Proc. IRE 50, 2081–2092 (1962). [CrossRef]
,16
R. K. Mongia and A. Ittipiboon, “Theoretical and experimental investigations on rectangular dielectric resonator antennas,” IEEE Trans. Antennas Propag. 45, 1348–1356 (1997). [CrossRef]
]: the problem is solved in the microwave region. A theoretical model, applicable under specific conditions, is available to benchmark the method.
The single nano-cuboid: the problem is solved in the optical region. We compare the numerical solution with the theoretical model proposed in (i).
The single spherical nanoparticle: we compare the numerical solution with the Mie solution [17
C. Bohren and D. Huffmann, Absorption and Scattering of Light by Small Particles (John Wiley, New York, 1983).
].
The optical dipole antenna of finite lateral dimensions: we realistically model the dipole antenna with two smoothly rounded arms. This geometry has traditionally been studied via the thin wire approximation [18
S. Ramo, J. R. Whinnery, and T. V. Duzer, Fields and Waves in Communication Electronics (John Wiley, New York, 1993).
]. The FDTD method [19
A. Dhawan, S. J. Norton, M. D. Gerhold, and T. Vo-Dinh, “Comparison of FDTD numerical computations and analytical multipole expansion method for plasmonics-active nanosphere dimers,” Opt. Express 17, 9688–9703 (2009). [CrossRef] [PubMed]
] experiences intrinsic difficulties when modeling rounded objects. Other work [20
X. Cui and D. Erni, “The influence of particle shapes on the optical response of nearly touching plasmonic nanoparticle dimers,” J. Comput. Theor. Nanosci. 47, 1610–1615 (2010). [CrossRef]
–22
J. Smajic, C. Hafner, L. Raguin, K. Tavzarashvili, and M. Mishrikey, “Comparison of numerical methods for the analysis of plasmonic structures,” J. Comput. Theor. Nanosci. 6, 763–774 (2009). [CrossRef]
] is restricted to 2-dimensional geometries. In [23
J. Smajic, C. Hafner, K. Tavzarashvili, and R. Vahldieck, “Numerical analysis of channel plasmon polaritons enhanced optical antennas,” J. Comput. Theor. Nanosci. 5, 725–734 (2008). [CrossRef]
, 24
C. G. Khoury, S. J. Norton, and T. Vo-Dinh, “Plasmonics of 3-D nanoshell dimers using multipole expansion and finite element method,” ACS Nano 3, 2776–2788 (2009). [CrossRef] [PubMed]
], the 3-D FEM is restricted to relatively small-scale problems due to the huge memory consumption and excessive computation time, further suffering from relatively low accuracy. In [25
A. M. Kern and O. J. F. Martin, “Modeling near-field properties of plasmonic nanoparticles: a surface integral approach,” in Plasmonic: Nanoimaging, Nanofabrication, and their Applications V, V. M. Shalaev and D. P. Tsai, eds., Proc. SPIE 7395, 739518 (2009).
, 26
A. M. Kern and O. J. F. Martin, “Excitation and reemission of molecules near realistic plasmonic nanostructures,” Nano Lett. 11, 482–487 (2011). [CrossRef] [PubMed]
], the surface integral equation (SIE) method experiences difficulties when modeling the substrate. Here, we study a technologically relevant geometrical discretization of the optical dipole antenna, free of the aforementioned limitations. Both bright and dark modes [27
M. W. Chu, V. Myroshnychenko, C. H. Chen, J. P. Deng, C. Y. Mou, and F. J. García de Abajo, “Probing bright and dark surface-plasmon modes in individual and coupled noble metal nanoparticles using an electron beam,” Nano Lett. 9, 399–404 (2009). [CrossRef]
] are investigated.

2. Formulation of the problem

The electromagnetic background of our simulations originates in the Maxwell equations combined with proper boundary conditions. In the time-harmonic regime, after eliminating the magnetic field H(x), the electric field E(x) satisfies

\begin{array}{l} \nabla \times (μ_{r}^{- 1} \nabla \times E (x)) + {i k}_{0} σ Z_{0} E (x) - k_{0}^{2} ε_{r} E (x) = 0, & \nabla \cdot (ε_{r} E (x)) = 0, & x \in Ω . \end{array}

(1)

Here, Ω ⊂ 𝕉³ is a bounded domain,

k_{0} (= \tilde{ω} \sqrt{μ_{0} ε_{0}})

is the complex wavenumber in free space; ω̃ (= ω + iα) is the complex angular frequency, ω is the angular frequency and α is the exponential decay rate; μ_r, ε_r are relative magnetic permeability and relative electric permittivity, respectively; σ is the ohmic loss of the material;

Z_{0} (= \sqrt{μ_{0} / ε_{0}})

is the characteristic impedance of free space.

We model the material loss of the nano metallic structure through the complex permittivity ε_r, normally with negative real part, instead of ohmic loss σ. The dielectric function is not generally constant, but ε_r depends on the frequency ω, or, equivalently, the real part of k₀ [28

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

]. In this paper, we use experimental data from Johnson and Christy [28

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

] and employ a piecewise linear interpolation. We obtain the dielectric function within the spectral range from 0.64 eV to 6.60 eV [28

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

]. The substrate’s, e.g., glass, or the environment’s, e.g. a covering layer of water, dielectric functions are not constant either [29

G. M. Hale and M. R. Querry, “Optical constants of water in the 200 nm to 200 μm wavelength region,” Appl. Opt. 12, 555–563 (1973). [CrossRef] [PubMed]

]. However, since the variation is usually small within a certain optical spectral range, we assume that the dielectric properties of the substrate and the environment are constant.

The formulation of the nano-optical device simulation uses a bounded domain Ω = Ω₁ ∪ Ω₂ ∪ Ω₃. It encloses the nano-optical device Ω₁, its substrate Ω₂, and the environment Ω₃, cf. Fig. 1. We solve the eigenproblem Eq. (1), where

ε_{r} = {\begin{array}{l} ε_{r} (Re (k_{0})), & x \in Ω_{1}, \\ ε_{sub}, & x \in Ω_{2}, \\ ε_{env}, & x \in Ω_{3} . \end{array}

(2)

ε_sub and ε_env are the constant relative permittivity for the substrate and the environment, respectively. ε_r depends on Re(k₀) in the subdomain Ω₁.

Fig. 1 A 3-D view of a typical computational domain Ω = Ω₁ ∪ Ω₂ ∪ Ω₃. Ω₁ is the optical device, Ω₂ is the substrate, Ω₃ is the environment. Γ₁₃ is the interface between Ω₁ and Ω₂ ∪ Ω₃, and Γ₂₃ is the interface between Ω₂ and Ω₁ ∪ Ω₃. Γ is the boundary of Ω.

On the interface Γ₁₃ the tangential component of E(x) must be continuous. On the artificial surface Γ we use the 1st order absorbing, aka. transparent, boundary condition (ABC) [12

J. Jin, The Finite Element Method in Electromagnetics (John Wiley, New York, 2002).

]

n \times \nabla \times E (x) = - i k n \times (n \times E (x)),

(3)

where

k = k_{0} \sqrt{μ_{r} ε_{r}}

is the wavenumber of the surrounding medium. We let Ω be a sphere of radius R_b. R_b should be as large as possible to minimize the artificial reflections on Γ introduced by the ABC, which arise in the case of non-orthogonal incidence of out-going wave. On the other hand, R_b should be small to limit the size of the discretized system. By numerical experiment we fix R_b such that a change of R_b does not result in a significant deviation of the computed resonant wavelength, cf. the discussion in Sections 4.3. Following this rule, we observe that the minimum distance between Γ and the device Ω₁ is at least (and usually larger than) 1/3 of the resonance wavelength.

We solve Eqs. (1)–(3) to determine the electric field E(x) and the wavenumber k₀. The corresponding magnetic field is obtained by

H (x) = \frac{1}{- i \tilde{ω} μ_{0} μ_{r}} \nabla \times E (x) = \frac{i}{k_{0} Z_{0} μ_{r}} \nabla \times E (x) .

(4)

The resonant frequency is f = ω/2π (ω = Re(ω̃)). The decay rate α (= Im(ω̃)) is directly connected to the loss. In analogy to an oscillating system with damping, the quality factor is defined as

Q_{1} = | ω / 2 α | .

(5)

We compute the system’s total energy U, stored energy U_s, dissipated energy U_d, and the radiated power P_r of the dispersive metallic nano-structure, respectively. Under the condition of no magnetic dispersion (μ_r = 1), the total energy is [30

R. Loudon, “The propagation of electromagnetic energy through an absorbing dielectric,” J. Phys. A: Gen. Phys. 3, 233245 (1970). [CrossRef]

, 31

R. Ruppin, “Electromagnetic energy density in a dispersive and absorptive material,” Phys. Lett. A 299, 309–312 (2002). [CrossRef]

]

U = \frac{ε_{0}}{4} \int_{Ω} (Re (ε_{r}) + \frac{2 ω Im (ε_{r})}{γ_{e}}) {| E (x) |}^{2} d x + \frac{μ_{0}}{4} \int_{Ω} {| H (x) |}^{2} d x .

(6)

where γ_e is the damping frequency of the dispersive metal with γ_e = τ⁻¹ where τ is the relaxation time. τ ≈ 31fs for silver and τ ≈ 9.3fs for gold [28

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

If (i) the metal dielectric function can be described by the Drude model [17

C. Bohren and D. Huffmann, Absorption and Scattering of Light by Small Particles (John Wiley, New York, 1983).

], and (ii) ωτ ≫ 1, we obtain the stored energy [32

T. G. Philbin, “Electromagnetic energy momentum in dispersive media,” Phys. Rev. A 83, 013823 (2011). [CrossRef]

, 33

F. D. Nunes, T. C. Vasconcelos, M. Bezerra, and J. Weiner, “Electromagnetic energy density in dispersive and dissipative media,” J. Opt. Soc. Am. B 28, 1544–1552 (2011). [CrossRef]

] in Eq. (7), equivalent to the Brillouin formula [34

L. Brillouin, Wave Propagation and Group Velocity (Academic Press, New York, 1960).

U_{s} = \frac{ε_{0}}{4} \int_{Ω} (Re (ε_{r}) + ω \frac{d (Re (ε_{r}))}{d ω} {| E (x) |}^{2} d x + \frac{μ_{0}}{4} \int_{Ω} {| H (x) |}^{2} d x .

(7)

We note that computing the stored energy for nano-metallic structures is still a matter of intense, scientific debate [32

T. G. Philbin, “Electromagnetic energy momentum in dispersive media,” Phys. Rev. A 83, 013823 (2011). [CrossRef]

,33

F. D. Nunes, T. C. Vasconcelos, M. Bezerra, and J. Weiner, “Electromagnetic energy density in dispersive and dissipative media,” J. Opt. Soc. Am. B 28, 1544–1552 (2011). [CrossRef]

]. We emphasize that Eq. (7) holds only if both assumptions (i) and (ii) are satisfied. The second assumption, i.e., ωτ ≫ 1, holds for all eigenmodes computed in this paper, e.g., ωτ > 100 for silver spheres, Section 4.3, and ωτ > 20 for gold dipole antennas, Section 4.4. While the Drude model is a good fit for dispersive silver, still, there is a difference between the analytical model and the experimental data from [28

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

], particularly noticeable in the range of 2.5eV to 3.5eV [33

F. D. Nunes, T. C. Vasconcelos, M. Bezerra, and J. Weiner, “Electromagnetic energy density in dispersive and dissipative media,” J. Opt. Soc. Am. B 28, 1544–1552 (2011). [CrossRef]

], cf. Section 4.3. For the optical dispersion of gold, the Drude model is a rough fit. In particular, the difference between the analytical model’s imaginary part of ε_r and the experimental data is considerable for the photon energy above 2 eV [35

A. Vial, A. S. Grimault, D. Macías, D. Barchiesi, and M. L. de la Chapelle, “Improved analytical fit of gold dispersion: application to the modeling of extinction spectra with a finite-difference time-domain method,” Phys. Rev. B 71, 085416 (2005). [CrossRef]

]. Therefore, in the spectral range above 2eV the calculation of the stored energy via Eq. (7) is invalid for gold, cf. discussion in Section 4.4.2. However, a study on how to accurately evaluate the stored energy for metals is definitely beyond the scope of this paper. Nevertheless, we employ Eq. (7) in order to assess its accuracy by comparing it to results derived via an alternative route. Dissipated energy [33

F. D. Nunes, T. C. Vasconcelos, M. Bezerra, and J. Weiner, “Electromagnetic energy density in dispersive and dissipative media,” J. Opt. Soc. Am. B 28, 1544–1552 (2011). [CrossRef]

] is then computed as

U_{d} = \frac{ε_{0}}{2} \int_{Ω_{1}} Im (ε_{r}) {| E (x) |}^{2} d x .

(8)

Radiated power is calculated via the Poynting theorem [18

S. Ramo, J. R. Whinnery, and T. V. Duzer, Fields and Waves in Communication Electronics (John Wiley, New York, 1993).

P_{r} = \frac{1}{2} \int_{Γ} Re (E (x) \times H {(x)}^{*}) \cdot d x .

(9)

Therefore the system’s quality factor can also be defined by

Q_{2} = \frac{ω U_{s}}{ω U_{d} + P_{r}} .

(10)

We comment that Q₂ will appear to be a less reliable definition than Q₁, since U_s in (7) may not be accurate enough for nano-metallic structures. Therefore, In our numerical examples, we list both Q₁ and Q₂ for comparison purposes. The radiative quantum yield η is evaluated as the ratio of radiation loss and total loss, i.e.,

η = \frac{P_{r}}{P_{r} + ω U_{d}} .

(11)

In the special case of a non-dispersive medium that has no dissipative loss, e.g., the dielectric resonator antenna (DRA) with constant Re(ε_r) and Im(ε_r) = 0, (7) and (8) lead to

\begin{matrix} U_{s} = \frac{ε_{0}}{2} \int_{Ω} Re (ε_{r}) {| E (x) |}^{2} d x, & U_{d} = 0, \end{matrix}

(12)

such that Q₂ = ωU_s/P_r and Q₁ = Q₂.

3. Numerical solution

3.1. The finite element method

The finite element method (FEM) is a suitable method for the nano-optical devices simulation, since such devices exhibit delicate geometries with a wide span of geometrical scales. In order to apply the finite element method we need a weak form of (1). We first note that for sufficiently smooth vector functions E(x) and f(x) that satisfy the boundary condition (3) we have [12

J. Jin, The Finite Element Method in Electromagnetics (John Wiley, New York, 2002).

, 18

S. Ramo, J. R. Whinnery, and T. V. Duzer, Fields and Waves in Communication Electronics (John Wiley, New York, 1993).

]

\begin{array}{l} \int_{Ω} f (x) \cdot (\nabla \times μ_{r}^{- 1} \nabla \times E (x)) d x & = & \int_{Ω} μ_{r}^{- 1} \nabla \times f \cdot \nabla \times E d x - \int_{Γ} μ_{r}^{- 1} n \cdot (f \times \nabla \times E) d x \\ = & \int_{Ω} μ_{r}^{- 1} \nabla \times f \cdot \nabla \times E d x + i \int_{Γ} k μ_{r}^{- 1} (n \times f) \cdot (n \times E) d x . \end{array}

Similarly, for any sufficiently smooth scalar function q(x) vanishing on Γ we get

\int_{Ω} ε_{r} (\nabla \cdot E (x)) q (x) d x = - \int_{Ω} ε_{r} E (x) \cdot (\nabla q (x)) d x = 0.

Therefore, a natural weak form of (1)–(3) is:

Find k₀ ∈ 𝔺 and E ∈ V, E ≠ 0, such that for all f ∈ V and all q ∈ W

\begin{array}{l} \int_{Ω} μ_{r}^{- 1} \nabla \times f \cdot \nabla \times E d x + i k_{0} [\int_{Γ} \sqrt{\frac{ε_{r}}{μ_{r}}} (n \times f) \cdot (n \times E) d x + \int_{Ω} σ Z_{0} f \cdot E d x] + k_{0}^{2} \int_{Ω} ε_{r} f \cdot E d x = 0, \\ \int_{Ω} ε_{r} E \cdot \nabla q d x = 0. \end{array}

(13)

Here, V denotes the functions in H(curl;Ω) that satisfy the boundary condition (3) and

W = H_{0}^{1} (Ω)

. For details on these function spaces see, e.g., [13

P. Monk, Finite Element Methods for Maxwell’s Equations (Oxford University Press, Oxford, 2003). [CrossRef] [PubMed]

] or [12

J. Jin, The Finite Element Method in Electromagnetics (John Wiley, New York, 2002).

We discretize (13) by the Ritz-Galerkin method [12

J. Jin, The Finite Element Method in Electromagnetics (John Wiley, New York, 2002).

, 13

P. Monk, Finite Element Methods for Maxwell’s Equations (Oxford University Press, Oxford, 2003). [CrossRef] [PubMed]

] employing appropriate finite element subspaces of V and W. To that end we triangulate Ω by tetrahedra with Gmsh [36

C. Geuzaine and J.-F. Remacle, “Gmsh: A 3-D finite element mesh generator with built-in pre- and postprocessing facilities,” Int. J. Numer. Methods Eng. 79, 1309–1331 (2009). [CrossRef]

]. The triangulation has to respect the interfaces Γ₁₃ and Γ₂₃. The mesh must be fine in the vicinity of the surface of the optical device Ω₁, but it can be coarse in the far-field zone. This strategy maintains efficiency and accuracy while, on the other hand, it reduces the computational cost.

The electric/magnetic vector functions in V are then approximated by Nédélec (edge) elements, while the scalar functions in W are approximated by Lagrange (nodal) finite elements [13

P. Monk, Finite Element Methods for Maxwell’s Equations (Oxford University Press, Oxford, 2003). [CrossRef] [PubMed]

]. This approach avoids the generation of spurious eigensolutions, and imposing the boundary conditions is straightforward [12

J. Jin, The Finite Element Method in Electromagnetics (John Wiley, New York, 2002).

Let the vector functions N_i, 1 ≤ i ≤ n, be the Nédélec basis functions, while the scalar functions N_ℓ, 1 ≤ ℓ ≤ m, denote the Lagrange basis functions. Then, Eq. (13) yield a constrained complex quadratic eigenvalue problem

\begin{array}{l} T (λ) x = A x + λ R x - λ^{2} M x = 0, \\ C^{T} x = 0, \end{array}

(14)

where λ (= k₀) is the eigenvalue and x is the eigenvector. The matrices A, R, M, and C in (14) have entries

\begin{array}{l} a_{i j} = \int_{Ω} μ_{r}^{- 1} (\nabla \times N_{i}) \cdot (\nabla \times N_{j}) d x, & 1 \leq i, j \leq n, \\ m_{i j} = \int_{Ω} ε_{r} N_{i} \cdot N_{j} d x, & 1 \leq i, j \leq n, \\ r_{i j} = i \int_{Γ} \sqrt{ε_{r} / μ_{r}} (n \times N_{i}) \cdot (n \times N_{J}) d x + i \int_{Ω} σ Z_{0} N_{i} \cdot N_{j} d x, & 1 \leq i, j \leq n, \\ c_{i ℓ} = \int_{Ω} ε_{r} N_{i} (x) \cdot \nabla N_{ℓ} (x) d x, & 1 \leq i \leq n, 1 \leq ℓ \leq m, \end{array}

(15)

respectively. Since ε_r and μ_r are element-wise constant, A, M, and R are composed of real symmetric element matrices multiplied by a complex factor. Thus, they are complex symmetric.

3.2. The eigensolver

Instead of solving the eigenvalue problem (14) we use an appropriate linearization [37

F. Tisseur and K. Meerbergen, “The quadratic eigenvalue problem,” SIAM Rev . 43, 235–286 (2001). [CrossRef]

]. Defining x₂ = λx₁, (14) is transformed into a constrained linear eigenproblem of the form

\begin{matrix} 𝒜 x \equiv [\begin{matrix} A & O \\ O & I \end{matrix}] (\begin{matrix} x_{1} \\ x_{2} \end{matrix}) = λ [\begin{matrix} R & M \\ I & O \end{matrix}] (\begin{matrix} x_{1} \\ x_{2} \end{matrix}) \equiv λ ℬ x, & 𝒞^{T} x \equiv {[\begin{matrix} C & O \\ O & C \end{matrix}]}^{T} x = 0 . \end{matrix}

(16)

The identity matrix I in (16) could be replaced by any nonsingular matrix X, i.e., we could write Xx₂ = λXx₁. Choosing X = M would make both 𝒜 and 𝒨 complex symmetric. Since complex symmetry does not necessarily increase numerical simplicity [38

P. Arbenz and M. E. Hochstenbach, “A Jacobi–Davidson method for solving complex symmetric eigenvalue problems,” SIAM J. Sci. Comput. 25, 1655–1673 (2004). [CrossRef]

] we set X = I.

The linearization of the quadratic eigenvalue problem allows us to use well-known eigen-solvers. For solving large sparse complex eigenvalue problems the Jacobi–Davidson QZ algorithm (JDQZ) is appropriate [39

Z. Bai, J. Demmel, J. Dongarra, A. Ruhe, and H. van der Vorst, Templates for the Solution of Algebraic Eigenvalue Problems: A Practical Guide (SIAM, Philadelphia PA, 2000). [CrossRef]

, 40

D. R. Fokkema, G. L. G. Sleijpen, and H. A. Van der Vorst, “Jacobi–Davidson style QR and QZ algorithms for the partial reduction of matrix pencils,” SIAM J. Sci. Comput. 20, 94–125 (1996). [CrossRef]

]. JDQZ computes a partial QZ decomposition

\begin{matrix} 𝒜 Q = {Z T}_{A}, & ℬ Q = {Z T}_{B}, & 𝒞^{T} Q = O, \end{matrix}

(17)

where Q and Z are 2n × r matrices, and T_A and T_B are upper-triangular r × r matrices. The quotients of corresponding diagonal elements of T_A and T_B provide the eigenvalues. We are looking for a few eigenvalues closest to a prescribed target value τ. The initial target is chosen to be an estimate of the expected eigenvalue. More details on selecting the initial target are given in Section 4.1. The columns of Q are the Schur vectors from which the desired eigenvectors are extracted. Z is a so-called shadow space. For details on eigenvalue/eigenvector extraction and on restarting we refer to [39

Z. Bai, J. Demmel, J. Dongarra, A. Ruhe, and H. van der Vorst, Templates for the Solution of Algebraic Eigenvalue Problems: A Practical Guide (SIAM, Philadelphia PA, 2000). [CrossRef]

]. Here, we only discuss the most time consuming step of the Jacobi–Davidson algorithm, the expansion of the search space. To that end, the so-called correction equation has to be solved,

\begin{matrix} (ℐ - \tilde{Z} {\tilde{Z}}^{*}) (𝒜 - \tilde{λ} 𝒝) (ℐ - \tilde{Q} {\tilde{Q}}^{*}) t = r & {\tilde{Q}}^{*} t = 0, & 𝒞^{T} t = 0, \end{matrix}

(18)

where the solution t is used to span the search space. λ̃ is the actual best eigenvalue approximation, Q̃ is the matrix Q expanded by the corresponding Schur vector u. Z̃ is the matrix Z expanded by the new shadow vector p. So, Q̃ and Z̃ are 2n × (r + 1) matrices.

The solution t of (18) is not needed to high accuracy. Therefore, we can approximately solve the correction equation by executing a few steps of a preconditioned Krylov space method. The Jacobi–Davidson method has the advantages of short execution time and low memory consumption, which is crucial for solving large scale eigenproblems. To impose the divergence-free condition in (14), we construct an appropriate projector to assert that each vector in the search space is in the null space of 𝒞^T [41

R. Geus, “The Jacobi–Davidson algorithm for solving large sparse symmetric eigenvalue problems.” PhD Thesis No. 14734, ETH Zurich 2002.

–43

P. Arbenz and R. Geus, “Multilevel preconditioned iterative eigensolvers for Maxwell eigenvalue problems,” Appl. Numer. Math. 54, 107–121 (2005). [CrossRef]

The global matrix of the correction Eq. (18) has a 2 × 2 block structure, cf. (16). We use the block Gauss–Seidel preconditioner on the global matrix. To solve with the diagonal (1,1) block, we use the Jacobi preconditioner or a hierarchical basis preconditioner [42

P. Arbenz, M. Bečka, R. Geus, U. Hetmaniuk, and T. Mengotti, “On a parallel multilevel preconditioned Maxwell eigensolver,” Parallel Comput. 32, 157–165 (2006). [CrossRef]

, 43

P. Arbenz and R. Geus, “Multilevel preconditioned iterative eigensolvers for Maxwell eigenvalue problems,” Appl. Numer. Math. 54, 107–121 (2005). [CrossRef]

]. The diagonal (2,2) block is related to the identity matrix, so it is trivial to solve.

3.3. Implementation

We have implemented the algorithm in our software package femaxx [41

R. Geus, “The Jacobi–Davidson algorithm for solving large sparse symmetric eigenvalue problems.” PhD Thesis No. 14734, ETH Zurich 2002.

] which is continuously being developed at the Swiss Federal Institute of Technology (ETH) Zurich and the Paul Scherrer Institute (PSI) in order to solve large scale electromagnetic eigenmodal problems. The code has been efficiently parallelized to profit from today’s widespread, low-cost distributed memory parallel computers. It targets solving various kinds of electromagnetic eigenmodal problems based on the finite element method in 3-dimensional space. Both linear and quadratic Nédélec elements are available. The femaxx eigensolver is based on the Trilinos software framework [44

Trilinos Project Home Page, http://trilinos.sandia.gov/.

] which defines basic parallel objects such as vectors and sparse matrices. The real-symmetric Jacobi–Davidson eigensolver has been implemented with efficient multilevel preconditioning [42

P. Arbenz, M. Bečka, R. Geus, U. Hetmaniuk, and T. Mengotti, “On a parallel multilevel preconditioned Maxwell eigensolver,” Parallel Comput. 32, 157–165 (2006). [CrossRef]

,43

P. Arbenz and R. Geus, “Multilevel preconditioned iterative eigensolvers for Maxwell eigenvalue problems,” Appl. Numer. Math. 54, 107–121 (2005). [CrossRef]

]. We have now added a complex Jacobi–Davidson eigensolver. At present, Trilinos does not support a robust complex package [44

Trilinos Project Home Page, http://trilinos.sandia.gov/.

]. Thus, our complex eigensolver treats real and imaginary parts of all quantities separately. The femaxx postprocessor evaluates the electromagnetic field at any spatial location and visualizes the result with Paraview [45

Paraview Home Page, http://www.paraview.org/.

4. Validation and application of the algorithm

To validate the eigenmodal solver algorithm and, in particular, its combination with an absorbing boundary condition (ABC), we analyze three problems, namely (i) the dielectric resonator antenna (DRA), (ii) the nano-cuboid, and (iii) a spherical nanoparticle with variable radius. Eventually, in addition to the benchmark cases, we apply the algorithm to the analysis of the optical dipole antenna.

Unless otherwise specified, quadratic Nédélec elements are used with each element having 20 local degrees of freedom. We use spatially highly resolved tetrahedral meshes. In particular, for the nano-optical application examples, the size of an element in the region of the resonating structure is in the order of a few nanometers, or smaller. In combination with second order basis functions we thus employ a highly accurate field approximation. Therefore, the error caused by the spatial discretization is generally negligible.

All simulations were carried out on the Cray XT5 at the Swiss National Supercomputing Centre (CSCS) [46

Home Page of the Swiss National Supercomputing Centre (CSCS), http://www.cscs.ch/.

]. We stress that femaxx generally operates on any reasonable distributed memory parallel computer. Naturally, given the size of the numerical problem at hand, femaxx indeed profits from increased compute power but it does by no means enforce the use of a supercomputer.

4.1. Dielectric resonator antenna

We analyze rectangular dielectric resonator antennas (DRAs) in the microwave region. DRAs have been analyzed theoretically [15

A. Okaya and L. F. Barash, “The dielectric microwave resonator,” Proc. IRE 50, 2081–2092 (1962). [CrossRef]

,16

R. K. Mongia and A. Ittipiboon, “Theoretical and experimental investigations on rectangular dielectric resonator antennas,” IEEE Trans. Antennas Propag. 45, 1348–1356 (1997). [CrossRef]

]. We study DRAs with various dimensions a, b, d, and dielectric function ε_r. The meshes usually contain ≈ 80′000 tetrahedra for each DRA, so that the discretization counts ≈ 600′000 degrees of freedom (dof).

The radius of the computational domain Ω is R_b = 30 mm, see Fig. 2(a). Theoretical solutions from [16

R. K. Mongia and A. Ittipiboon, “Theoretical and experimental investigations on rectangular dielectric resonator antennas,” IEEE Trans. Antennas Propag. 45, 1348–1356 (1997). [CrossRef]

] and numerical results for the

T E_{111}^{z}

mode are given in Table 1. For all simulations, we set τ to a real value, the equivalent of 4.8 GHz. Each simulation completes in ≈ 3 minutes using 64 cores on the Cray XT5. The normalized residual (||𝒜x − λ𝒝x||₂/||x||₂) of the converged eigenpair is in the order of 10⁻⁵.

Fig. 2 (a) a 3-D sketch of the DRA (or nano-cuboid). (b): a DRA with a = b = 10 mm, d = 4 mm and ε_r = 20.0, resonance at 6.545GHz. The electric field distribution |E| of the

T E_{111}^{z}

mode in the antenna and its surrounding region is shown. (c): a gold cuboid with a = b = 100 nm, d = 40 nm, resonance at 935.1THz (where ε_r = −1.2308 + i5.8458). It shows the electric field distribution |E| of the

T E_{111}^{z}

mode in the cuboid and the surrounding region. (b) and (c) are visualized on the xy-plane through the center of the antenna. (d) 3 components of the electric field of the DRA (in (b)) plotted along the y-axis. (e) 3 components of the electric field of the gold cuboid (in (c)) plotted along the y-axis. Note that, E_x, E_y, and E_z are complex. For E_y and E_z (approximately zero), only real parts are shown. E_x is scaled such that ||Re(E_x)||_∞ = 1.

Table 1

T E_{111}^{z}

mode of DRAs (f : resonant frequency computed by femaxx; f_MI: theoretically determined resonant frequency [16

R. K. Mongia and A. Ittipiboon, “Theoretical and experimental investigations on rectangular dielectric resonator antennas,” IEEE Trans. Antennas Propag. 45, 1348–1356 (1997). [CrossRef]

]; Diff.: relative difference of f and f_MI; Q_MI: theoretically determined the Q factor [16

R. K. Mongia and A. Ittipiboon, “Theoretical and experimental investigations on rectangular dielectric resonator antennas,” IEEE Trans. Antennas Propag. 45, 1348–1356 (1997). [CrossRef]

].)

ε_r	a (mm)	b (mm)	d (mm)	f (GHz)	f_MI (GHz)	Diff. (%)	Q₁	Q₂	Q_MI

79.46	7.45	7.45	2.98	4.644	4.346	6.9	176.8	174.4	128.7
37.84	8.60	8.60	2.58	6.221	5.934	4.8	31.1	31.3	37.9
37.84	7.45	7.45	3.51	5.614	5.337	5.2	50.8	51.0	44.0
20.0	10	10	4	6.545	6.409	2.1	13.3	13.4	18.1

We note that the theory in [16

R. K. Mongia and A. Ittipiboon, “Theoretical and experimental investigations on rectangular dielectric resonator antennas,” IEEE Trans. Antennas Propag. 45, 1348–1356 (1997). [CrossRef]

] is based on a simplified model restricted to Ω₁, i.e. the DRA, assuming perfect magnetic boundary conditions (PMC) E·n = 0 on the surface Γ₁₃ of the DRA. In femaxx Γ₁₃ is treated as an interface where just the continuity of the tangential component of E is enforced. E · n need not be zero.

As expected, Q₁ matches very well with Q₂. However, the different treatment of the DRA walls causes a noticeable difference between f and f_MI and also between Q₁ and Q_MI. We note that, interestingly, our frequencies f are closer to the experimental data cited in [16

R. K. Mongia and A. Ittipiboon, “Theoretical and experimental investigations on rectangular dielectric resonator antennas,” IEEE Trans. Antennas Propag. 45, 1348–1356 (1997). [CrossRef]

] than f_MI. Thus, the model implemented in femaxx appears to be more realistic.

We select the last DRA as an example for the convergence behavior of the solution w.r.t. the number N of tetrahedral mesh elements. We let N be 8′137, 28′150, 87′380 and 645′868, respectively. The computed resonant frequencies are, respectively, 6.548 GHz, 6.545 GHz, 6.545 GHz and 6.545 GHz, while the corresponding quality factors Q₁ are, respectively, 13.9, 13.4, 13.3 and 13.3. Thus, a tetrahedral mesh containing ≈ 80′000 elements for each DRA in Table 1 is good enough.

The electric field distribution |E| of the last DRA in Table 1 (a = b = 10 mm, d = 4 mm) is shown in Fig. 2(b). Here, the number of tetrahedra of the mesh is 645′868. According to the theoretical model in [16

R. K. Mongia and A. Ittipiboon, “Theoretical and experimental investigations on rectangular dielectric resonator antennas,” IEEE Trans. Antennas Propag. 45, 1348–1356 (1997). [CrossRef]

], the electric field inside the DRA along the y-axis satisfies E_y = E_z = 0 and E_x = Ak_y sin(k_yy) where A is a complex constant and k_y = π/b is the wavenumber. We plot the three components of the electric field along the y-axis in Fig. 2(d). They are in good agreement with the theoretical model, except that E_x varies a bit more than half a cycle inside the DRA. This, again, can be ascribed to the differing treatments of the DRA surface.

For DRA’s, in the microwave region of the electromagnetic spectrum, the ε_r of the employed materials is generally constant. However, in the optical region of the electromagnetic spectrum the dielectric permittivity is in general not constant but depends on the real part of eigenvalue Re(k₀), cf. (2). The dispersive material property implies that the matrices M and C in (15) also depend on Re(k₀). So, in the nano-optical region, the eigenproblem (14) becomes a ‘truly’ non-linear, i.e. not linearizable, eigenproblem. At present, femaxx is capable of solving linear and quadratic eigenproblems. Nevertheless, we can address this non-linear eigenproblem via solving a sequence of quadratic eigenproblems. We start with an estimate λ⁽⁰⁾ for the eigenvalue of the expected eigenmode. By means of Re(λ⁽⁰⁾), we determine the dielectric permittivity ε_r = ε_r(Re(λ⁽⁰⁾)), based on [28

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

], so that the matrices M(Re(λ⁽⁰⁾)), and C(Re(λ⁽⁰⁾)) can be constructed. Then, we solve the quadratic eigenvalue problem (14) and consequently obtain an improved eigenvalue estimate λ⁽¹⁾. This procedure is repeated to yield further estimates λ⁽²⁾, λ⁽³⁾, etc. We terminate the iteration as soon as two consecutive estimates are close enough. In what follows, we use this iterative procedure. We first validate our method with the nano-cuboid and nano-sphere geometries. Then we use femaxx for the analysis of specific optical devices. In the individual quadratic eigenvalue problems we use τ = Re (λ⁽^j⁾) as the target.

4.2. Cuboid

We investigate the

T E_{111}^{z}

mode of a nano-cuboid and reuse the sketch in Fig. 2(a). We let a = b = 100 nm, d = 40 nm, R_b = 300 nm and assume gold and silver, respectively, for the material of the nano-cuboid. We also reuse the theoretical DRA model [15

A. Okaya and L. F. Barash, “The dielectric microwave resonator,” Proc. IRE 50, 2081–2092 (1962). [CrossRef]

,16

R. K. Mongia and A. Ittipiboon, “Theoretical and experimental investigations on rectangular dielectric resonator antennas,” IEEE Trans. Antennas Propag. 45, 1348–1356 (1997). [CrossRef]

] to compute the resonance frequency f_MI and the electromagnetic field. The theoretically determined wavenumber k₀ is complex. The quality factor Q_MI is also theoretically determined via Eq. (5). We comment that Mongia’s model was developed in the microregion of the electromagnetic spectrum. Thus, strictly speaking, it may not be used as a reference or benchmark in the optical region where the dielectric permittivity of metals is highly dispersive and exhibits considerable dielectric loss. Nevertheless, we believe that, at present, there is no better analytical model available than Mongia’s for cuboid geometries. Therefore, we use Mongia’s model as a signpost into the optical region. Fully aware that we stretch Mongia’s model beyond its original limits of applicability, we consider it useful to help acquire at least qualitative insight into the mode structure of the nano-cuboid. The results are shown in Table 2. The mesh contains 645′868 tetrahedra leading to 4′120′676 degrees of freedom. For all simulations, the initial value of τ (= Re(λ⁽⁰⁾)) is equivalent of 1000 THz. The eigensolver calculation time for a single run is less than 10 minutes with 256 cores, while total computation time for one eigenmode is then below 40 minutes for 3–4 runs. The normalized residual of the converged eigenpair is around 10⁻⁵.

Table 2 Numerical figures of merit for the

T E_{111}^{z}

mode of the nano-cuboid where λ : resonant wavelength computed by femaxx; f : resonant frequency computed by femaxx; f_MI: theoretically determined resonant frequency [16

R. K. Mongia and A. Ittipiboon, “Theoretical and experimental investigations on rectangular dielectric resonator antennas,” IEEE Trans. Antennas Propag. 45, 1348–1356 (1997). [CrossRef]

]; Diff.: relative difference of f and f_MI; Q_MI: theoretically determined Q factor [16

R. K. Mongia and A. Ittipiboon, “Theoretical and experimental investigations on rectangular dielectric resonator antennas,” IEEE Trans. Antennas Propag. 45, 1348–1356 (1997). [CrossRef]

material	λ (nm)	f (THz)	f_MI (THz)	Diff. (%)	Q₁	Q₂	Q_MI	η

gold	320.6	935.1	739.8	26.4	0.47	0.82	0.39	5.9 × 10⁻⁶
silver	200.7	1494	1107	35.0	0.56	0.79	0.55	2.5 × 10⁻⁷

The electric field distribution |E| of the gold cuboid is shown in Fig. 2(c). The three components of E, plotted along the y-axis, are displayed in Fig. 2(e). Again they are in good agreement with the fields of the theoretical model. In Fig. 2(e), E_y ≈ 0 and E_z ≈ 0. E_x varies by a bit more than half a cycle inside the cuboid which indicates that the PMC is not strictly satisfied on the surfaces of the cuboid. We also notice that the field is strong inside the cuboid and drops to almost zero in the far-field zone, see Fig. 2(e). Therefore, the mode experiences rather high dissipative loss, but radiates only scant electromagnetic power. Thus, η is very small.

We comment that the dielectric function ε_r is not at all constant for gold and silver in the optical region of the spectrum, The dispersive material property, in combination with the different treatment of the DRA wall, largely explains the noticeable difference between f and f_MI. Compared to dielectric resonator antennas, the quality factors (Q₁, Q₂, or Q_MI) of the nano-cuboid are very small, significantly below 1. This is attributed to high dissipative loss. The difference between Q₁ and Q₂ is noticeable, cf. the discussion in Section 2.

4.3. Sphere

We compute the resonant, plasmonic mode of a silver sphere, hovering in vacuum, with varying radius, i.e. R = 25,30,35,40,60 and 70 nm. Thus, ε_sub = ε_env = 1.0. The radius of the computational domain is R_b = 300 nm. For each sphere, the mesh contains about 250′000 tetrahedra. For all simulations, the initial value of τ is equivalent of 3.00eV. In this example we had problems with the quadratic elements. The computation of E_res1 for R = 25 nm and of E_res2 (see below) for R = 60 nm failed to converge such that we had to resort to linear elements. In general, quadratic elements deliver more accurate solutions for fewer elements but they are more demanding w.r.t. iterative solver schemes. The eigensolver calculation time for a single run is less than 10 minutes with 256 cores (or 32 cores) if quadratic (or linear) elements are used. The normalized residual of the converged eigenpair is around 5 × 10⁻⁴. The total computation time for one eigenmode is thus below 1 hour for 5–6 runs. We study the energy range from 2 eV to 4 eV, in wavelength units, 310 nm to 620 nm. We use λ for the numerically determined resonance wavelength and the resonance energy E = hc/λ, with h the Plank constant and c the speed of light in vacuum, respectively.

For comparison, we also compute the theoretical Mie solution [17

C. Bohren and D. Huffmann, Absorption and Scattering of Light by Small Particles (John Wiley, New York, 1983).

], for each radius and for the same energy range. When R = 25, 30, 35, or 40 nm, the peaks of scattering (Q_sca) and absorption efficiency (Q_abs) occur at almost the same position,

E_{res 1}^{mie}

. In Fig. 3(a) we show the Mie solution for R = 30 nm. The full-linewidth-at-half-maximum (FWHM) energy width ΔE_ext of the extinction efficiency, Q_ext = Q_sca + Q_abs, can be evaluated. The quality factor is computed as

Q_{res1}^{mie} = E_{res 1}^{mie} / Δ E_{ext}

[47

C. Sönnichsen, T. Franzl, T. Wilk, G. von Plessen, and J. Feldmann, “Drastic reduction of plasmon damping in gold nanorods,” Phys. Rev. Lett. 88, 077402 (2002) [CrossRef] [PubMed]

Fig. 3 Mie solutions of silver spheres with radii of (a) 30 nm and (b) 60 nm.

When R = 60 or 70 nm, the peaks of Q_sca and Q_abs are separated from each other, such that there are two distinct resonances close together. Figure 3(b) shows the Mie solution for R = 60 nm. We first study the resonance

E_{res 1}^{mie}

where Q_sca reaches maximum. Let ΔE_sca be the FWHM of Q_sca, and the quality factor of this resonance

Q_{res 1}^{mie} = E_{res 1}^{mie} / Δ E_{sca}

. However, in such cases where two peaks overlap, Q_sca deviates from a Lorentzian shape [48

R. Fuchs and K. L. Kliewer, “Optical modes of vibration in an ionic crystal spheres,” J. Opt. Soc. Am. 58, 319–330 (1968) [CrossRef]

], so that evaluating

Q_{res 1}^{mie}

by FWHM ΔE_sca is not accurate.

The radiative quantum yield η^mie, computed with the Mie solution, is defined as the ratio between scattering and extinction efficiency, i.e., η^mie = Q_sca/Q_ext.

In Table 3 and Fig. 4, we show the numerical results, compared with Mie solutions. In general, the resonance (E_res1 and

E_{res 1}^{mie}

) and radiative quantum yield (η and η^mie) show good agreements. The differences are mainly attributed to the 1st order ABC employed by femaxx, while the Mie model is unbounded. Quality factors (Q₁ and

Q_{res 1}^{mie}

) also match when R = 25,30,35 or 40 nm. When R = 60 or 70 nm, two peaks of resonances (

E_{res 1}^{mie}

and

E_{res 2}^{mie}

) exist and are close together,

Q_{res 1}^{mie}

is not accurate, see discussion above, and Q₁ is more reliable. As noted in Section 2, in the range of 2.5eV to 3.5eV, the difference of silver ε_r between the Drude model and the experimental data is noticeable. Therefore Q₂ is not accurate either. Note that E_res1 approximates a 3-fold degenerate resonance. The corresponding computed mode splits in three simple modes since the mesh does not respect the spherical symmetry of the problem. For R = 60 nm, the relative differences among the three computed resonances are below 1%.

Fig. 4 Numerical and Mie solutions of Q factors and radiative quantum yield (subfigure) of silver spheres with varying radius R.

Table 3 Results of the numerical analysis and the semi-analytical Mie solutions of the silver sphere hovering in vacuum, where: E_res1: resonant energy computed by femaxx; E_diff: the relative difference between E_res1 and

E_{res 1}^{mie}

. Q_diff: the relative difference between Q₁ and

Q_{res 1}^{mie}

R (nm)	λ (nm)	E_res1 (eV)	$E_{res 1}^{mie}$ (eV)	E_diff (%)	η	η^mie	Q₁	Q₂	$Q_{res 1}^{mie}$	Q_diff (%)
25	368.7	3.36	3.41	1.5	0.57	0.48	21.5	23.2	22.7	5.3
30	373.4	3.32	3.37	1.5	0.64	0.63	17.0	20.7	16.9	0.6
35	381.0	3.25	3.33	2.4	0.78	0.75	10.1	19.3	12.3	17.9
40	388.4	3.19	3.25	1.9	0.81	0.86	9.6	17.1	8.8	9.1
60	408.9	3.03	3.02	0.3	0.88	0.95	6.0	12.3	3.0	100
70	440.6	2.81	2.81	0	0.95	0.98	5.3	8.6	2.1	152

We use R = 60nm as an example (see Fig. 3(b)), and study the resonance

E_{res 2}^{mie}

where Q_abs reaches the maximum. (Our computations indicate that E_res2 is 5-fold degenerate.) The numerical value of this mode is E_res2 = 3.48eV and Q₁ = 8.1. Let the angular frequency, dissipated energy and radiated powers of resonance E_res2 be ω_res2,

U_{d}^{res 2}

and

P_{r}^{res 2}

, respectively. Further, let the corresponding quantities for resonance E_res1 be ω_res1,

U_{d}^{res 1}

and

P_{r}^{res 1}

. The electric fields are scaled so that ∫_Ω Re (ε_r)|E(x)|² dx are equal. We define

r_{d} = (ω_{res 2} U_{d}^{res 2}) / (ω_{res 1} U_{d}^{res 1})

and

r_{r} = P_{r}^{res 2} / P_{r}^{res 1}

. Then r_d = 8.9 and r_r = 0.21. Resonance E_res2 has higher dissipative loss since Q_abs reaches maximum, while resonance E_res1 has higher radiation loss since Q_sca reaches maximum. The Mie solution of this mode is

E_{res 2}^{mie} = 3.45 eV

. Again, evaluating

Q_{res 2}^{mie} (= 17.2)

by FWHM ΔE_abs is not accurate. The electric field plots of the two modes, when R = 60nm, are shown in Figures 5(a) and 5(b).

Fig. 5 (a) and (b) show the electric field distributions (|E|) in the vicinity of a silver sphere on the xy-plane. The radius of the sphere is 60 nm. (a): the mode when Q_sca reaches maximum, E_res = 3.03 eV; (b): the mode when Q_abs reaches maximum, E_res = 3.48 eV.

We use the resonance E_res1 of the silver sphere (R = 60nm) to study the influence of the 1st order ABC on the accuracy of the result. If we let the boundary radius R_b vary from 250 to 300nm (keeping the number of tetrahedra around 250′000), the deviations of the computed resonant wavelength are as small as 3nm. For R_b = 300nm, the computed resonant wavelength (408.9nm) is closest to the resonant wavelength of the Mie solution (≈ 410.5nm). For smaller R_b, the error gets larger. For R_b = 200nm, e.g., the computed resonant wavelength is 385.0nm.

4.4. An optical dipole antenna fabricated with gold dielectric permittivity

We study the electromagnetic near-field of the gold nano-optical dipole antenna investigated in [25

A. M. Kern and O. J. F. Martin, “Modeling near-field properties of plasmonic nanoparticles: a surface integral approach,” in Plasmonic: Nanoimaging, Nanofabrication, and their Applications V, V. M. Shalaev and D. P. Tsai, eds., Proc. SPIE 7395, 739518 (2009).

]. The geometry is shown in Fig. 6(a) and a 3-D sketch is in 6(b). The antenna has two arms, each of which has dimensions a = b = 40 mm, l = 100 nm. The corners of each arm are rounded, with a radius of curvature = 5 nm. The gap g between the two arms is 20 nm wide but will be varied in order to study the electromagnetic field enhancement as a function of the gap width. The radius R_b of the spherical computational domain is 400 nm.

Fig. 6 The nano-optical dipole antenna. (a) geometry; (b) 3-D sketch.

We study three different dielectric material arrangements: (1) the antenna hovers in vacuum, thus ε_sub = ε_env = 1.0; (2) the antenna resides on a silica substrate, thus ε_sub = 2.25,ε_env = 1.0; (3) the antenna resides on a silica substrate and is covered with water, thus ε_sub = 2.25, ε_env = 1.77 [29

G. M. Hale and M. R. Querry, “Optical constants of water in the 200 nm to 200 μm wavelength region,” Appl. Opt. 12, 555–563 (1973). [CrossRef] [PubMed]

]. We understand this arrangement to be a relatively detailed model for an optical antenna used in biochemical sensing where many processes are studied in aqueous solutions.

We start by considering a single arm in vacuum whose resonant wavelength is 623.5 nm. When bringing two such arms into close proximity, the interaction of the two fundamental eigenmodes lead to mode splitting. When these two arms come close enough to each other, the so-called bright, aka. optically active, mode and the dark mode split. The modes can then be well distinguished in the scattering spectra [27

M. W. Chu, V. Myroshnychenko, C. H. Chen, J. P. Deng, C. Y. Mou, and F. J. García de Abajo, “Probing bright and dark surface-plasmon modes in individual and coupled noble metal nanoparticles using an electron beam,” Nano Lett. 9, 399–404 (2009). [CrossRef]

, 49

S. C. Yang, H. Kobori, C. L. He, M. H. Lin, H. Y. Chen, C. Li, M. Kanehara, T. Teranishi, and S. Gwo, “Plasmon hybridization in individual gold nanocrystal dimers: direct observation of bright and dark modes,” Nano Lett. 10, 632–637 (2010). [CrossRef] [PubMed]

4.4.1. Bright mode

We compute the bright mode within the spectral range of 400 nm to 1000 nm. For all simulations, the initial value of τ is equivalent of 660nm. The results are presented in Table 4. The fields of model (1), an antenna in vacuum, are evaluated in the sample domain Ω_sample = 400nm×200nm×200 nm. Both magnetic and electric fields are shown in Fig. 7(a) and 7(b), respectively. A close-up view of the electric field in the gap region is shown in Fig. 7(c). The results in Table 4 confirm the influence of the substrate onto the antenna’s resonant wavelength. With increasing permittivity ε_sub the antenna resonance is red-shifted; with increasing ε_env the resonance wavelength is also red-shifted [50

H. Fischer and O. J. F. Martin, “Engineering the optical response of plasmonic nanoantennas,” Opt. Express 16, 9144–9154 (2008). [CrossRef] [PubMed]

]. Each single run of the quadratic eigensolver requires 10 to 15 minutes. The total computation time for one bright mode is then below 90 minutes for 5–6 runs. The normalized residual of the converged eigenpair is around 5 × 10⁻⁴. The memory per core required by the eigensolver is 33.38 MB for the model (1), and 28.65 MB for models (2) and (3). Thus, the method is particularly attractive when compared to other 3-dimensional FEM approaches [22

J. Smajic, C. Hafner, L. Raguin, K. Tavzarashvili, and M. Mishrikey, “Comparison of numerical methods for the analysis of plasmonic structures,” J. Comput. Theor. Nanosci. 6, 763–774 (2009). [CrossRef]

, 24

C. G. Khoury, S. J. Norton, and T. Vo-Dinh, “Plasmonics of 3-D nanoshell dimers using multipole expansion and finite element method,” ACS Nano 3, 2776–2788 (2009). [CrossRef] [PubMed]

]. From Figs. 7(b) and 7(c) we note that a high-intensity electric field is generated in the vicinity of the antenna surface. It is stronger in the gap (Fig. 7(d)), and it is strongest at the inner corners of the antenna.

Fig. 7 (a)–(c): Numerical analysis of a gold, nano-optical dipole antenna (model (1)). The gap width is 20 nm and the resonance is at 665.7 nm; (a) magnetic and (b) electric field distribution (|H| and |E|) in the vicinity of the antenna surface; (c) electric field distribution |E| around the gap; (a), (b), and (c) are evaluated on the xy-plane through the center of the antenna arms. (d) plots the electric field amplitudes |E| along the x-axis, associated with 3 material arrangements.

Table 4 Numerical analysis of a gold nano-optical antenna (#p: number of cores; N: number of tetrahedra; dofs: degrees of freedom; λ : resonant wavelength computed by femaxx; Q₁ and Q₂: quality factors by Eq. (5) and Eq. (10), respectively;

model	ε_sub	ε_env	#p	N	dofs	λ (nm)	Q₁	Q₂

(1) in vacuum	1.0	1.0	300	340’269	2’162’346	665.7	20.0	19.8
(2) on substrate	2.25	1.0	600	597’408	3’791’532	717.5	12.4	13.6
(3) on substrate, water-covered	2.25	1.77	600	597’408	3’791’532	801.7	10.6	14.3

We now investigate how the gap width g affects the resonance wavelength λ and the field intensity |E|² in the gap. We use the dielectric arrangement of model (2). Let the gap width g be 20 nm, 10 nm, 5 nm, and 1 nm, respectively. The corresponding resonance wavelengths are then 717.5 nm, 751.9 nm, 796.3 nm, and 913.5 nm, respectively. A red shift, i.e., increasing wavelength, of the antenna resonance occurs with decreasing gap widths. Since the gap widths are significantly smaller than the resonance wavelengths, we can approximate the gap region with an electrostatic regime. Thus, we understand the gap region, including those portions of the dipole antenna, that immediately border it, to form a plate capacitor whose capacitance is given as C = ε₀ε_rA/d. Here, A, d, and ε_r correspond to the plate area, the distance between the plates and the dielectric property of the medium between the plates, respectively; in particular, ε_r = 1 since the gap is in vacuum. Then, it becomes clear that, with decreasing gap width, capacitance increases. A capacitance increase on the other hand is roughly proportional to a higher value of the dielectric permittivity and thus, for constant length of the dipole arms, the resonance wavelength increases since the resonance frequency is reduced. We comment that adding a capacitor parallel to an antenna’s terminal contacts is a well-known technique in microwave electronics in order to geometrically shorten the antenna, thus mimicking higher dielectric permittivity in the vicinity of the antenna.

On the other hand, when decreasing the gap width g = 10 nm, 5 nm, and 1 nm, the corresponding field intensities |E|² increase and are 3.4, 9.6, and 56 times larger than when a gap width g = 20 nm is employed, see Fig. 8(b). Decreasing the gap width thus dramatically increases the field intensity in the gap. An example field distribution, for a 5 nm gap width, is shown in Fig. 8(a). Due to the substrate effect the field is not symmetric with respect to the x-axis.

Fig. 8 Numerical analysis of a gold nano-optical antenna (model (2)): (a) gap width = 5 nm, resonance at 796.3 nm. It shows the electric field distribution |E| (visualized on the xy-plane) in the vicinity of the antenna surface; (b) plots of the electric field amplitude |E| along the x-axis with varying gap width g.

4.4.2. Dark mode

We study the dark mode of the dipole antenna in vacuum, model (1). The respective charge profiles [27

] of the bright and the dark modes are shown in Fig. 9(a). Let the angular frequency, dissipative energy and radiated powers of the dark mode be ω_dark,

U_{d}^{dark}

and

P_{r}^{dark}

, respectively; further, let the corresponding quantities for the bright mode be ω_bright,

U_{d}^{bright}

and

P_{r}^{bright}

. Again, the electric fields are scaled so that the integrals of the type ∫_Ω Re(ε_r)|E(x)|² dx are equal. We define

r_{d} = (ω_{dark} U_{d}^{dark}) / (ω_{bright} U_{d}^{bright})

and

r_{r} = P_{r}^{dark} / P_{r}^{bright}

. The radiative quantum yield η is computed via Eq. (11). Associated results are listed in Table 5. The mesh contains about 350′000 tetrahedra which results in more than 2′000′000. For all simulations of the dark modes, the initial value of τ is equivalent of 550nm. The eigensolver calculation time for a single run is 10 to 15 minutes with 256 cores, and the normalized residual of the converged eigenpair is around 5 × 10⁻⁴ for the bright mode and 10⁻⁴ for the dark mode. The total computation time for one eigenmode is then below 90 minutes for 5–6 runs.

Fig. 9 (a) Charge profiles; (b) the dark mode of the gold nano-optical antenna (gap width g = 10 nm) in vacuum: resonance at 557.4 nm. The electric field distribution |E| (visualized on the xy-plane) in the vicinity of the antenna surface is shown; (c) the gold nano-optical antenna (gap width g = 5 nm) in vacuum. The electric field amplitude |E| along the x-axis is plotted, associated with the bright mode (resonance at 713.9 nm) and the dark mode (resonance at 532.1 nm), respectively.

Table 5 Dark modes vs. bright modes

g (nm)	dark mode				bright mode				r_d	r_r
g (nm)	λ (nm)	η	Q₁	Q₂	λ (nm)	η	Q₁	Q₂	r_d	r_r
20	582.0	2.4 × 10⁻³	0.89	13.3	665.7	0.27	20.0	19.8	10.1	0.067
10	557.4	7.0 × 10⁻⁴	0.78	9.7	689.8	0.29	23.3	20.3	18.2	0.032
5	532.1	2.3 × 10⁻⁴	0.65	7.4	713.9	0.26	21.5	20.7	44.9	0.029

When the width g of the gap between the two arms is reduced, the single arm’s degenerate fundamental mode with its resonance at 623.5 nm splits into a bright and a dark mode. With decreasing gap width, the bright mode red-shifts, due to the effect of increased capacity for smaller width, and the dark mode shifts to blue. The dark mode radiates less, i.e. (r_r < 0.1), which implies that it cannot efficiently radiate energy into the far-field, thus earning its dark reputation. On the other hand, the dark mode experiences higher dissipative loss inside the two gold arms, i.e. r_d > 10. Its dissipative loss approximates almost 100% of the total loss, namely η < 0.005, and thus the quality factor Q₁ of the dark mode is smaller than the bright mode.

The resonances of all dark modes are above 2eV, i.e. below 620nm; on the other hand, the resonances of all bright modes are all below 2eV, i.e. above 620nm. Therefore, computing the stored energy with Eq. (7) is a good approximation for bright modes, but significantly less so for dark modes, cf. the discussion in section 2. As a consequence, Q₂ is in good agreement with Q₁ for bright modes, but is significantly off for dark modes.

An exemplary field distribution for the dark mode, associated with a 10 nm gap width, is shown in Fig. 9(b). In Fig. 9(c) the field plots of bright and dark modes for a 5 nm gap width are displayed. The field associated with a dark mode is not only generated in the vicinity of the antenna surface, but also inside the gold arms, corresponding to a large

U_{d}^{dark}

. The field intensity is quite low in the gap.

5. Discussion and conclusions

We have implemented a 3-dimensional finite element method for the full-wave numerical analysis of electromagnetic resonators, surrounded by free-space. Specifically, there are no cavity walls, and the resonator structures need not be enclosed within perfect electric conductor (PEC) boundary conditions. Instead, we need to apply an absorbing boundary condition (ABC) scheme. The associated quadratic eigenvalue problem is linearized so that we can employ our current implementation of the Jacobi–Davidson QZ algorithm.

We have analyzed four different electromagnetic problems. The dielectric resonator antenna (DRA), analyzed in the microwave region of the spectrum, clearly demonstrates that the proposed Jacobi–Davidson algorithm is capable of calculating the dominant mode with excellent accuracy. The differences between the theoretically determined resonance frequencies, electric field and Q factor and their numerical counterparts are attributed to the fact that, in the analytical model, the perfect magnetic conductor (PMC) condition, is imposed on all surfaces of DRA, while femaxx does not need this simplifying assumption. The theoretical model of the DRA is reused as a signpost into the optical region for simulations of a nano-cuboid. The

T E_{111}^{z}

mode of the nano-cuboid has a field pattern similar to the DRA, but exhibits a rather small quality factor due to high dissipative loss in metals. There is a noticeable difference between the numerically determined resonance wavelength and quality factor and the analytically determined resonance wavelength and quality factor, respectively. On the other hand, the electromagnetic field solutions demonstrate amazing agreement, in particular when we consider that Mongia’s model has been stretched considerably beyond its limits of applicability. The analysis of electromagnetic scattering off spheres represents a widely used benchmark in nano-optics which a computational electrodynamics code must be able to solve. In general, we observe good agreement between the theoretically determined resonance frequency, radiative quantum yield and Q factor and their numerical counterparts.

Thus, now that femaxx has passed a series of benchmarks, we apply it to more complex geometries and configurations. One of them, the optical dipole antenna, represents an important, elementary building block in many detection and sensing devices, especially for the life and materials sciences. We observe excellent agreement with previously published results by Kern and Martin [25

], who simulate a dipole antenna with 20 nm gap width. We emphasize that our algorithm allows for the simple inclusion of background media and also media covering the antenna, e.g. an aqueous solution covering the antenna. This corresponds to a setup that is representative for many experimental arrangements. We note that the surface integral equation (SIE) method is also capable of modeling background media but at the cost of considerably increased complexity since it must include additional Green’s functions to accommodate, e.g., a substrate. Meanwhile, the finite element method ‘absorbs’ all this additional geometrical and material complexity into the corresponding mesh. On the other hand, the finite element method needs a transparent boundary condition of high-quality which should be placed as close as possible to the resonating structure, thus ensuring that precious mesh elements are not wasted but are put to good use when modeling fine geometrical details. We also note that the dark mode, with little radiated power, is another ideal test for our approach. A precise understanding of the dark modes is important for plasmonic applications [49

]. Indeed, we have found the dark mode, starting from the original single arm configuration with a degenerate fundamental mode, which also strengthens the confidence into our solver. In conclusion, our approach is flexible enough to analyze nano-optical devices with almost arbitrary geometry or material arrangements. The femaxx code is therefore a useful instrument for the analysis of technologically relevant nano-optical device concepts. We continue to develop the code.

Acknowledgments

It is our pleasure to thank Olivier Martin, Andreas Kern (EPFL, NAM) and Arya Fallahi (PSI) for sharing thoughts on the theoretical analysis of the nanosphere and optical dipole antenna. We also thank an anonymous reviewer for the thoughtful and supportive review. The work of the first author (H. Guo) was supported in part by grant no. 200021-117978 of the Swiss National Science Foundation.

References and links

1.	L. Novotny and N. van Hulst, “Antennas for light,” Nat. Photon. 5, 83–90 (2011). [CrossRef]
2.	L. Novotny, “Effective wavelength scaling for optical antennas,” Phys. Rev. Lett. 98, 266802 (2007). [CrossRef] [PubMed]
3.	P. Bharadwaj, B. Deutsch, and L. Novotny, “Optical antennas,” Adv. Opt. Photon. 1, 438–483 (2009). [CrossRef]
4.	W. L. Barnes, A. Dereux, and T. W. Ebbesen, “Surface plasmon subwavelength optics,” Nature 424, 824–830 (2003). [CrossRef] [PubMed]
5.	H. Aouani, O. Mahboub, N. Bonod, E. Devaux, E. Popov, H. Rigneault, T. W. Ebbesen, and J. Wenger, “Bright unidirectional fluorescence emission of molecules in a nanoaperture with plasmonic corrugations,” Nano Lett. 11, 637–644 (2011). [CrossRef] [PubMed]
6.	D. W. Pohl, S. G. Rodrigo, and L. Novotny, “Stacked optical antennas,” Appl. Phys. Lett. 98, 023111 (2011). [CrossRef]
7.	A. Kinkhabwala, Z. Yu, S. Fan, Y. Avlasevich, K. Müllen, and W. E. Moerner, “Large single-molecule fluorescence enhancements produced by a bowtie nanoantenna,” Nat. Photon. 3, 654–657 (2009). [CrossRef]
8.	M. F. Garcia-Parajo, “Optical antennas focus in on biology,” Nat. Photon. 2, 201–203, (2008). [CrossRef]
9.	P. Nagpal, N. C. Lindquist, S.-H. Oh, and D. J. Norris, “Ultrasmooth patterned metals for plasmonics and metamaterials,” Science 325, 594–597 (2009). [CrossRef] [PubMed]
10.	G. W. Hanson and A. B. Yakovlev, Operator Theory for Electromagnetics – An Introduction (Springer, 2002). [PubMed]
11.	Ch. Hafner, The Generalized Multipole Technique for Computational Electromagnetics (Artech House Books, Boston, 1990).
12.	J. Jin, The Finite Element Method in Electromagnetics (John Wiley, New York, 2002).
13.	P. Monk, Finite Element Methods for Maxwell’s Equations (Oxford University Press, Oxford, 2003). [CrossRef] [PubMed]
14.	J. L. Volakis, A. Chatterjee, and L. C. Kempel, Finite Element Method for Electromagnetics – Antennas, Microwave Circuits and Scattering Applications (IEEE Press, New York, 1998).
15.	A. Okaya and L. F. Barash, “The dielectric microwave resonator,” Proc. IRE 50, 2081–2092 (1962). [CrossRef]
16.	R. K. Mongia and A. Ittipiboon, “Theoretical and experimental investigations on rectangular dielectric resonator antennas,” IEEE Trans. Antennas Propag. 45, 1348–1356 (1997). [CrossRef]
17.	C. Bohren and D. Huffmann, Absorption and Scattering of Light by Small Particles (John Wiley, New York, 1983).
18.	S. Ramo, J. R. Whinnery, and T. V. Duzer, Fields and Waves in Communication Electronics (John Wiley, New York, 1993).
19.	A. Dhawan, S. J. Norton, M. D. Gerhold, and T. Vo-Dinh, “Comparison of FDTD numerical computations and analytical multipole expansion method for plasmonics-active nanosphere dimers,” Opt. Express 17, 9688–9703 (2009). [CrossRef] [PubMed]
20.	X. Cui and D. Erni, “The influence of particle shapes on the optical response of nearly touching plasmonic nanoparticle dimers,” J. Comput. Theor. Nanosci. 47, 1610–1615 (2010). [CrossRef]
21.	J. M. McMahon, A. I. Henry, K. L. Wustholz, M. J. Natan, R. G. Freeman, R. P. Van Duyne, and G. C. Schatz, “Gold nanoparticle dimer plasmonics: finite element method calculations of the electromagnetic enhancement to surface-enhanced Raman spectroscopy,” Anal. Bioanal. Chem. 394, 1819–1825 (2009). [CrossRef] [PubMed]
22.	J. Smajic, C. Hafner, L. Raguin, K. Tavzarashvili, and M. Mishrikey, “Comparison of numerical methods for the analysis of plasmonic structures,” J. Comput. Theor. Nanosci. 6, 763–774 (2009). [CrossRef]
23.	J. Smajic, C. Hafner, K. Tavzarashvili, and R. Vahldieck, “Numerical analysis of channel plasmon polaritons enhanced optical antennas,” J. Comput. Theor. Nanosci. 5, 725–734 (2008). [CrossRef]
24.	C. G. Khoury, S. J. Norton, and T. Vo-Dinh, “Plasmonics of 3-D nanoshell dimers using multipole expansion and finite element method,” ACS Nano 3, 2776–2788 (2009). [CrossRef] [PubMed]
25.	A. M. Kern and O. J. F. Martin, “Modeling near-field properties of plasmonic nanoparticles: a surface integral approach,” in Plasmonic: Nanoimaging, Nanofabrication, and their Applications V, V. M. Shalaev and D. P. Tsai, eds., Proc. SPIE 7395, 739518 (2009).
26.	A. M. Kern and O. J. F. Martin, “Excitation and reemission of molecules near realistic plasmonic nanostructures,” Nano Lett. 11, 482–487 (2011). [CrossRef] [PubMed]
27.	M. W. Chu, V. Myroshnychenko, C. H. Chen, J. P. Deng, C. Y. Mou, and F. J. García de Abajo, “Probing bright and dark surface-plasmon modes in individual and coupled noble metal nanoparticles using an electron beam,” Nano Lett. 9, 399–404 (2009). [CrossRef]
28.	P. B. Johnson and R. W. Christy, “Optical constants of the noble metals,” Phys. Rev. B 6, 4370–4379 (1972). [CrossRef]
29.	G. M. Hale and M. R. Querry, “Optical constants of water in the 200 nm to 200 μm wavelength region,” Appl. Opt. 12, 555–563 (1973). [CrossRef] [PubMed]
30.	R. Loudon, “The propagation of electromagnetic energy through an absorbing dielectric,” J. Phys. A: Gen. Phys. 3, 233245 (1970). [CrossRef]
31.	R. Ruppin, “Electromagnetic energy density in a dispersive and absorptive material,” Phys. Lett. A 299, 309–312 (2002). [CrossRef]
32.	T. G. Philbin, “Electromagnetic energy momentum in dispersive media,” Phys. Rev. A 83, 013823 (2011). [CrossRef]
33.	F. D. Nunes, T. C. Vasconcelos, M. Bezerra, and J. Weiner, “Electromagnetic energy density in dispersive and dissipative media,” J. Opt. Soc. Am. B 28, 1544–1552 (2011). [CrossRef]
34.	L. Brillouin, Wave Propagation and Group Velocity (Academic Press, New York, 1960).
35.	A. Vial, A. S. Grimault, D. Macías, D. Barchiesi, and M. L. de la Chapelle, “Improved analytical fit of gold dispersion: application to the modeling of extinction spectra with a finite-difference time-domain method,” Phys. Rev. B 71, 085416 (2005). [CrossRef]
36.	C. Geuzaine and J.-F. Remacle, “Gmsh: A 3-D finite element mesh generator with built-in pre- and postprocessing facilities,” Int. J. Numer. Methods Eng. 79, 1309–1331 (2009). [CrossRef]
37.	F. Tisseur and K. Meerbergen, “The quadratic eigenvalue problem,” SIAM Rev . 43, 235–286 (2001). [CrossRef]
38.	P. Arbenz and M. E. Hochstenbach, “A Jacobi–Davidson method for solving complex symmetric eigenvalue problems,” SIAM J. Sci. Comput. 25, 1655–1673 (2004). [CrossRef]
39.	Z. Bai, J. Demmel, J. Dongarra, A. Ruhe, and H. van der Vorst, Templates for the Solution of Algebraic Eigenvalue Problems: A Practical Guide (SIAM, Philadelphia PA, 2000). [CrossRef]
40.	D. R. Fokkema, G. L. G. Sleijpen, and H. A. Van der Vorst, “Jacobi–Davidson style QR and QZ algorithms for the partial reduction of matrix pencils,” SIAM J. Sci. Comput. 20, 94–125 (1996). [CrossRef]
41.	R. Geus, “The Jacobi–Davidson algorithm for solving large sparse symmetric eigenvalue problems.” PhD Thesis No. 14734, ETH Zurich 2002.
42.	P. Arbenz, M. Bečka, R. Geus, U. Hetmaniuk, and T. Mengotti, “On a parallel multilevel preconditioned Maxwell eigensolver,” Parallel Comput. 32, 157–165 (2006). [CrossRef]
43.	P. Arbenz and R. Geus, “Multilevel preconditioned iterative eigensolvers for Maxwell eigenvalue problems,” Appl. Numer. Math. 54, 107–121 (2005). [CrossRef]
44.	Trilinos Project Home Page, http://trilinos.sandia.gov/.
45.	Paraview Home Page, http://www.paraview.org/.
46.	Home Page of the Swiss National Supercomputing Centre (CSCS), http://www.cscs.ch/.
47.	C. Sönnichsen, T. Franzl, T. Wilk, G. von Plessen, and J. Feldmann, “Drastic reduction of plasmon damping in gold nanorods,” Phys. Rev. Lett. 88, 077402 (2002) [CrossRef] [PubMed]
48.	R. Fuchs and K. L. Kliewer, “Optical modes of vibration in an ionic crystal spheres,” J. Opt. Soc. Am. 58, 319–330 (1968) [CrossRef]
49.	S. C. Yang, H. Kobori, C. L. He, M. H. Lin, H. Y. Chen, C. Li, M. Kanehara, T. Teranishi, and S. Gwo, “Plasmon hybridization in individual gold nanocrystal dimers: direct observation of bright and dark modes,” Nano Lett. 10, 632–637 (2010). [CrossRef] [PubMed]
50.	H. Fischer and O. J. F. Martin, “Engineering the optical response of plasmonic nanoantennas,” Opt. Express 16, 9144–9154 (2008). [CrossRef] [PubMed]

OCIS Codes
(140.4780) Lasers and laser optics : Optical resonators
(240.6680) Optics at surfaces : Surface plasmons
(260.3910) Physical optics : Metal optics
(260.5740) Physical optics : Resonance
(220.4241) Optical design and fabrication : Nanostructure fabrication

ToC Category:
Optics at Surfaces

History
Original Manuscript: January 18, 2012
Revised Manuscript: February 10, 2012
Manuscript Accepted: February 10, 2012
Published: February 21, 2012

Virtual Issues
Vol. 7, Iss. 4 Virtual Journal for Biomedical Optics

Citation
Hua Guo, Benedikt Oswald, and Peter Arbenz, "3-dimensional eigenmodal analysis of plasmonic nanostructures," Opt. Express 20, 5481-5500 (2012)
http://www.opticsinfobase.org/oe/abstract.cfm?URI=oe-20-5-5481

Sort: Author | Year | Journal | Reset

References

L. Novotny and N. van Hulst, “Antennas for light,” Nat. Photon.5, 83–90 (2011). [CrossRef]
L. Novotny, “Effective wavelength scaling for optical antennas,” Phys. Rev. Lett.98, 266802 (2007). [CrossRef] [PubMed]
P. Bharadwaj, B. Deutsch, and L. Novotny, “Optical antennas,” Adv. Opt. Photon.1, 438–483 (2009). [CrossRef]
W. L. Barnes, A. Dereux, and T. W. Ebbesen, “Surface plasmon subwavelength optics,” Nature424, 824–830 (2003). [CrossRef] [PubMed]
H. Aouani, O. Mahboub, N. Bonod, E. Devaux, E. Popov, H. Rigneault, T. W. Ebbesen, and J. Wenger, “Bright unidirectional fluorescence emission of molecules in a nanoaperture with plasmonic corrugations,” Nano Lett.11, 637–644 (2011). [CrossRef] [PubMed]
D. W. Pohl, S. G. Rodrigo, and L. Novotny, “Stacked optical antennas,” Appl. Phys. Lett.98, 023111 (2011). [CrossRef]
A. Kinkhabwala, Z. Yu, S. Fan, Y. Avlasevich, K. Müllen, and W. E. Moerner, “Large single-molecule fluorescence enhancements produced by a bowtie nanoantenna,” Nat. Photon.3, 654–657 (2009). [CrossRef]
M. F. Garcia-Parajo, “Optical antennas focus in on biology,” Nat. Photon.2, 201–203, (2008). [CrossRef]
P. Nagpal, N. C. Lindquist, S.-H. Oh, and D. J. Norris, “Ultrasmooth patterned metals for plasmonics and metamaterials,” Science325, 594–597 (2009). [CrossRef] [PubMed]
G. W. Hanson and A. B. Yakovlev, Operator Theory for Electromagnetics – An Introduction (Springer, 2002). [PubMed]
Ch. Hafner, The Generalized Multipole Technique for Computational Electromagnetics (Artech House Books, Boston, 1990).
J. Jin, The Finite Element Method in Electromagnetics (John Wiley, New York, 2002).
P. Monk, Finite Element Methods for Maxwell’s Equations (Oxford University Press, Oxford, 2003). [CrossRef] [PubMed]
J. L. Volakis, A. Chatterjee, and L. C. Kempel, Finite Element Method for Electromagnetics – Antennas, Microwave Circuits and Scattering Applications (IEEE Press, New York, 1998).
A. Okaya and L. F. Barash, “The dielectric microwave resonator,” Proc. IRE50, 2081–2092 (1962). [CrossRef]
R. K. Mongia and A. Ittipiboon, “Theoretical and experimental investigations on rectangular dielectric resonator antennas,” IEEE Trans. Antennas Propag.45, 1348–1356 (1997). [CrossRef]
C. Bohren and D. Huffmann, Absorption and Scattering of Light by Small Particles (John Wiley, New York, 1983).
S. Ramo, J. R. Whinnery, and T. V. Duzer, Fields and Waves in Communication Electronics (John Wiley, New York, 1993).
A. Dhawan, S. J. Norton, M. D. Gerhold, and T. Vo-Dinh, “Comparison of FDTD numerical computations and analytical multipole expansion method for plasmonics-active nanosphere dimers,” Opt. Express17, 9688–9703 (2009). [CrossRef] [PubMed]
X. Cui and D. Erni, “The influence of particle shapes on the optical response of nearly touching plasmonic nanoparticle dimers,” J. Comput. Theor. Nanosci.47, 1610–1615 (2010). [CrossRef]
J. M. McMahon, A. I. Henry, K. L. Wustholz, M. J. Natan, R. G. Freeman, R. P. Van Duyne, and G. C. Schatz, “Gold nanoparticle dimer plasmonics: finite element method calculations of the electromagnetic enhancement to surface-enhanced Raman spectroscopy,” Anal. Bioanal. Chem.394, 1819–1825 (2009). [CrossRef] [PubMed]
J. Smajic, C. Hafner, L. Raguin, K. Tavzarashvili, and M. Mishrikey, “Comparison of numerical methods for the analysis of plasmonic structures,” J. Comput. Theor. Nanosci.6, 763–774 (2009). [CrossRef]
J. Smajic, C. Hafner, K. Tavzarashvili, and R. Vahldieck, “Numerical analysis of channel plasmon polaritons enhanced optical antennas,” J. Comput. Theor. Nanosci.5, 725–734 (2008). [CrossRef]
C. G. Khoury, S. J. Norton, and T. Vo-Dinh, “Plasmonics of 3-D nanoshell dimers using multipole expansion and finite element method,” ACS Nano3, 2776–2788 (2009). [CrossRef] [PubMed]
A. M. Kern and O. J. F. Martin, “Modeling near-field properties of plasmonic nanoparticles: a surface integral approach,” in Plasmonic: Nanoimaging, Nanofabrication, and their Applications V, V. M. Shalaev and D. P. Tsai, eds., Proc. SPIE7395, 739518 (2009).
A. M. Kern and O. J. F. Martin, “Excitation and reemission of molecules near realistic plasmonic nanostructures,” Nano Lett.11, 482–487 (2011). [CrossRef] [PubMed]
M. W. Chu, V. Myroshnychenko, C. H. Chen, J. P. Deng, C. Y. Mou, and F. J. García de Abajo, “Probing bright and dark surface-plasmon modes in individual and coupled noble metal nanoparticles using an electron beam,” Nano Lett.9, 399–404 (2009). [CrossRef]
P. B. Johnson and R. W. Christy, “Optical constants of the noble metals,” Phys. Rev. B6, 4370–4379 (1972). [CrossRef]
G. M. Hale and M. R. Querry, “Optical constants of water in the 200 nm to 200 μm wavelength region,” Appl. Opt.12, 555–563 (1973). [CrossRef] [PubMed]
R. Loudon, “The propagation of electromagnetic energy through an absorbing dielectric,” J. Phys. A: Gen. Phys.3, 233245 (1970). [CrossRef]
R. Ruppin, “Electromagnetic energy density in a dispersive and absorptive material,” Phys. Lett. A299, 309–312 (2002). [CrossRef]
T. G. Philbin, “Electromagnetic energy momentum in dispersive media,” Phys. Rev. A83, 013823 (2011). [CrossRef]
F. D. Nunes, T. C. Vasconcelos, M. Bezerra, and J. Weiner, “Electromagnetic energy density in dispersive and dissipative media,” J. Opt. Soc. Am. B28, 1544–1552 (2011). [CrossRef]
L. Brillouin, Wave Propagation and Group Velocity (Academic Press, New York, 1960).
A. Vial, A. S. Grimault, D. Macías, D. Barchiesi, and M. L. de la Chapelle, “Improved analytical fit of gold dispersion: application to the modeling of extinction spectra with a finite-difference time-domain method,” Phys. Rev. B71, 085416 (2005). [CrossRef]
C. Geuzaine and J.-F. Remacle, “Gmsh: A 3-D finite element mesh generator with built-in pre- and postprocessing facilities,” Int. J. Numer. Methods Eng.79, 1309–1331 (2009). [CrossRef]
F. Tisseur and K. Meerbergen, “The quadratic eigenvalue problem,” SIAM Rev. 43, 235–286 (2001). [CrossRef]
P. Arbenz and M. E. Hochstenbach, “A Jacobi–Davidson method for solving complex symmetric eigenvalue problems,” SIAM J. Sci. Comput.25, 1655–1673 (2004). [CrossRef]
Z. Bai, J. Demmel, J. Dongarra, A. Ruhe, and H. van der Vorst, Templates for the Solution of Algebraic Eigenvalue Problems: A Practical Guide (SIAM, Philadelphia PA, 2000). [CrossRef]
D. R. Fokkema, G. L. G. Sleijpen, and H. A. Van der Vorst, “Jacobi–Davidson style QR and QZ algorithms for the partial reduction of matrix pencils,” SIAM J. Sci. Comput.20, 94–125 (1996). [CrossRef]
R. Geus, “The Jacobi–Davidson algorithm for solving large sparse symmetric eigenvalue problems.” PhD Thesis No. 14734, ETH Zurich 2002.
P. Arbenz, M. Bečka, R. Geus, U. Hetmaniuk, and T. Mengotti, “On a parallel multilevel preconditioned Maxwell eigensolver,” Parallel Comput.32, 157–165 (2006). [CrossRef]
P. Arbenz and R. Geus, “Multilevel preconditioned iterative eigensolvers for Maxwell eigenvalue problems,” Appl. Numer. Math.54, 107–121 (2005). [CrossRef]
Trilinos Project Home Page, http://trilinos.sandia.gov/ .
Paraview Home Page, http://www.paraview.org/ .
Home Page of the Swiss National Supercomputing Centre (CSCS), http://www.cscs.ch/ .
C. Sönnichsen, T. Franzl, T. Wilk, G. von Plessen, and J. Feldmann, “Drastic reduction of plasmon damping in gold nanorods,” Phys. Rev. Lett.88, 077402 (2002) [CrossRef] [PubMed]
R. Fuchs and K. L. Kliewer, “Optical modes of vibration in an ionic crystal spheres,” J. Opt. Soc. Am.58, 319–330 (1968) [CrossRef]
S. C. Yang, H. Kobori, C. L. He, M. H. Lin, H. Y. Chen, C. Li, M. Kanehara, T. Teranishi, and S. Gwo, “Plasmon hybridization in individual gold nanocrystal dimers: direct observation of bright and dark modes,” Nano Lett.10, 632–637 (2010). [CrossRef] [PubMed]
H. Fischer and O. J. F. Martin, “Engineering the optical response of plasmonic nanoantennas,” Opt. Express16, 9144–9154 (2008). [CrossRef] [PubMed]

Aouani, H.
Arbenz, P.
Avlasevich, Y.
Bai, Z.
Barash, L. F.
Barchiesi, D.
Barnes, W. L.
Becka, M.
Bezerra, M.
Bharadwaj, P.
Bohren, C.
Bonod, N.
Brillouin, L.
Chatterjee, A.
Chen, C. H.
Chen, H. Y.
Christy, R. W.
Chu, M. W.
Cui, X.
de la Chapelle, M. L.
Demmel, J.
Deng, J. P.
Dereux, A.
Deutsch, B.
Devaux, E.
Dhawan, A.
Dongarra, J.
Duzer, T. V.
Ebbesen, T. W.
Erni, D.
Fan, S.
Feldmann, J.
Fischer, H.
Fokkema, D. R.
Franzl, T.
Freeman, R. G.
Fuchs, R.
García de Abajo, F. J.
Garcia-Parajo, M. F.
Gerhold, M. D.
Geus, R.
Geuzaine, C.
Grimault, A. S.
Gwo, S.
Hafner, C.
Hafner, Ch.
Hale, G. M.
Hanson, G. W.
He, C. L.
Henry, A. I.
Hetmaniuk, U.
Hochstenbach, M. E.
Huffmann, D.
Ittipiboon, A.
Jin, J.
Johnson, P. B.
Kanehara, M.
Kempel, L. C.
Kern, A. M.
Khoury, C. G.
Kinkhabwala, A.
Kliewer, K. L.
Kobori, H.
Li, C.
Lin, M. H.
Lindquist, N. C.
Loudon, R.
Macías, D.
Mahboub, O.
Martin, O. J. F.
McMahon, J. M.
Meerbergen, K.
Mengotti, T.
Mishrikey, M.
Moerner, W. E.
Mongia, R. K.
Monk, P.
Mou, C. Y.
Müllen, K.
Myroshnychenko, V.
Nagpal, P.
Natan, M. J.
Norris, D. J.
Norton, S. J.
Novotny, L.
Nunes, F. D.
Oh, S.-H.
Okaya, A.
Philbin, T. G.
Pohl, D. W.
Popov, E.
Querry, M. R.
Raguin, L.
Ramo, S.
Remacle, J.-F.
Rigneault, H.
Rodrigo, S. G.
Ruhe, A.
Ruppin, R.
Schatz, G. C.
Sleijpen, G. L. G.
Smajic, J.
Sönnichsen, C.
Tavzarashvili, K.
Teranishi, T.
Tisseur, F.
Vahldieck, R.
van der Vorst, H.
Van der Vorst, H. A.
Van Duyne, R. P.
van Hulst, N.
Vasconcelos, T. C.
Vial, A.
Vo-Dinh, T.
Volakis, J. L.
von Plessen, G.
Weiner, J.
Wenger, J.
Whinnery, J. R.
Wilk, T.
Wustholz, K. L.
Yakovlev, A. B.
Yang, S. C.
Yu, Z.

ACS Nano

C. G. Khoury, S. J. Norton, and T. Vo-Dinh, “Plasmonics of 3-D nanoshell dimers using multipole expansion and finite element method,” ACS Nano3, 2776–2788 (2009). [CrossRef] [PubMed]

Adv. Opt. Photon.

P. Bharadwaj, B. Deutsch, and L. Novotny, “Optical antennas,” Adv. Opt. Photon.1, 438–483 (2009). [CrossRef]

Anal. Bioanal. Chem.

J. M. McMahon, A. I. Henry, K. L. Wustholz, M. J. Natan, R. G. Freeman, R. P. Van Duyne, and G. C. Schatz, “Gold nanoparticle dimer plasmonics: finite element method calculations of the electromagnetic enhancement to surface-enhanced Raman spectroscopy,” Anal. Bioanal. Chem.394, 1819–1825 (2009). [CrossRef] [PubMed]

Appl. Numer. Math.

P. Arbenz and R. Geus, “Multilevel preconditioned iterative eigensolvers for Maxwell eigenvalue problems,” Appl. Numer. Math.54, 107–121 (2005). [CrossRef]

Appl. Opt.

G. M. Hale and M. R. Querry, “Optical constants of water in the 200 nm to 200 μm wavelength region,” Appl. Opt.12, 555–563 (1973). [CrossRef] [PubMed]

Appl. Phys. Lett.

D. W. Pohl, S. G. Rodrigo, and L. Novotny, “Stacked optical antennas,” Appl. Phys. Lett.98, 023111 (2011). [CrossRef]

IEEE Trans. Antennas Propag.

R. K. Mongia and A. Ittipiboon, “Theoretical and experimental investigations on rectangular dielectric resonator antennas,” IEEE Trans. Antennas Propag.45, 1348–1356 (1997). [CrossRef]

Int. J. Numer. Methods Eng.

C. Geuzaine and J.-F. Remacle, “Gmsh: A 3-D finite element mesh generator with built-in pre- and postprocessing facilities,” Int. J. Numer. Methods Eng.79, 1309–1331 (2009). [CrossRef]

J. Comput. Theor. Nanosci.

J. Smajic, C. Hafner, L. Raguin, K. Tavzarashvili, and M. Mishrikey, “Comparison of numerical methods for the analysis of plasmonic structures,” J. Comput. Theor. Nanosci.6, 763–774 (2009). [CrossRef]
J. Smajic, C. Hafner, K. Tavzarashvili, and R. Vahldieck, “Numerical analysis of channel plasmon polaritons enhanced optical antennas,” J. Comput. Theor. Nanosci.5, 725–734 (2008). [CrossRef]
X. Cui and D. Erni, “The influence of particle shapes on the optical response of nearly touching plasmonic nanoparticle dimers,” J. Comput. Theor. Nanosci.47, 1610–1615 (2010). [CrossRef]

J. Opt. Soc. Am.

R. Fuchs and K. L. Kliewer, “Optical modes of vibration in an ionic crystal spheres,” J. Opt. Soc. Am.58, 319–330 (1968) [CrossRef]

J. Opt. Soc. Am. B

F. D. Nunes, T. C. Vasconcelos, M. Bezerra, and J. Weiner, “Electromagnetic energy density in dispersive and dissipative media,” J. Opt. Soc. Am. B28, 1544–1552 (2011). [CrossRef]

J. Phys. A: Gen. Phys.

R. Loudon, “The propagation of electromagnetic energy through an absorbing dielectric,” J. Phys. A: Gen. Phys.3, 233245 (1970). [CrossRef]

Nano Lett.

A. M. Kern and O. J. F. Martin, “Excitation and reemission of molecules near realistic plasmonic nanostructures,” Nano Lett.11, 482–487 (2011). [CrossRef] [PubMed]
M. W. Chu, V. Myroshnychenko, C. H. Chen, J. P. Deng, C. Y. Mou, and F. J. García de Abajo, “Probing bright and dark surface-plasmon modes in individual and coupled noble metal nanoparticles using an electron beam,” Nano Lett.9, 399–404 (2009). [CrossRef]
H. Aouani, O. Mahboub, N. Bonod, E. Devaux, E. Popov, H. Rigneault, T. W. Ebbesen, and J. Wenger, “Bright unidirectional fluorescence emission of molecules in a nanoaperture with plasmonic corrugations,” Nano Lett.11, 637–644 (2011). [CrossRef] [PubMed]
S. C. Yang, H. Kobori, C. L. He, M. H. Lin, H. Y. Chen, C. Li, M. Kanehara, T. Teranishi, and S. Gwo, “Plasmon hybridization in individual gold nanocrystal dimers: direct observation of bright and dark modes,” Nano Lett.10, 632–637 (2010). [CrossRef] [PubMed]

Nat. Photon.

L. Novotny and N. van Hulst, “Antennas for light,” Nat. Photon.5, 83–90 (2011). [CrossRef]
A. Kinkhabwala, Z. Yu, S. Fan, Y. Avlasevich, K. Müllen, and W. E. Moerner, “Large single-molecule fluorescence enhancements produced by a bowtie nanoantenna,” Nat. Photon.3, 654–657 (2009). [CrossRef]
M. F. Garcia-Parajo, “Optical antennas focus in on biology,” Nat. Photon.2, 201–203, (2008). [CrossRef]

Nature

W. L. Barnes, A. Dereux, and T. W. Ebbesen, “Surface plasmon subwavelength optics,” Nature424, 824–830 (2003). [CrossRef] [PubMed]

Opt. Express

Parallel Comput.

P. Arbenz, M. Bečka, R. Geus, U. Hetmaniuk, and T. Mengotti, “On a parallel multilevel preconditioned Maxwell eigensolver,” Parallel Comput.32, 157–165 (2006). [CrossRef]

Phys. Lett. A

R. Ruppin, “Electromagnetic energy density in a dispersive and absorptive material,” Phys. Lett. A299, 309–312 (2002). [CrossRef]

Phys. Rev. A

T. G. Philbin, “Electromagnetic energy momentum in dispersive media,” Phys. Rev. A83, 013823 (2011). [CrossRef]

Phys. Rev. B

P. B. Johnson and R. W. Christy, “Optical constants of the noble metals,” Phys. Rev. B6, 4370–4379 (1972). [CrossRef]
A. Vial, A. S. Grimault, D. Macías, D. Barchiesi, and M. L. de la Chapelle, “Improved analytical fit of gold dispersion: application to the modeling of extinction spectra with a finite-difference time-domain method,” Phys. Rev. B71, 085416 (2005). [CrossRef]

Phys. Rev. Lett.

C. Sönnichsen, T. Franzl, T. Wilk, G. von Plessen, and J. Feldmann, “Drastic reduction of plasmon damping in gold nanorods,” Phys. Rev. Lett.88, 077402 (2002) [CrossRef] [PubMed]
L. Novotny, “Effective wavelength scaling for optical antennas,” Phys. Rev. Lett.98, 266802 (2007). [CrossRef] [PubMed]

Proc. IRE

A. Okaya and L. F. Barash, “The dielectric microwave resonator,” Proc. IRE50, 2081–2092 (1962). [CrossRef]

Science

P. Nagpal, N. C. Lindquist, S.-H. Oh, and D. J. Norris, “Ultrasmooth patterned metals for plasmonics and metamaterials,” Science325, 594–597 (2009). [CrossRef] [PubMed]

SIAM J. Sci. Comput.

P. Arbenz and M. E. Hochstenbach, “A Jacobi–Davidson method for solving complex symmetric eigenvalue problems,” SIAM J. Sci. Comput.25, 1655–1673 (2004). [CrossRef]
D. R. Fokkema, G. L. G. Sleijpen, and H. A. Van der Vorst, “Jacobi–Davidson style QR and QZ algorithms for the partial reduction of matrix pencils,” SIAM J. Sci. Comput.20, 94–125 (1996). [CrossRef]

SIAM Rev

F. Tisseur and K. Meerbergen, “The quadratic eigenvalue problem,” SIAM Rev. 43, 235–286 (2001). [CrossRef]

Other

L. Brillouin, Wave Propagation and Group Velocity (Academic Press, New York, 1960).
A. M. Kern and O. J. F. Martin, “Modeling near-field properties of plasmonic nanoparticles: a surface integral approach,” in Plasmonic: Nanoimaging, Nanofabrication, and their Applications V, V. M. Shalaev and D. P. Tsai, eds., Proc. SPIE7395, 739518 (2009).
G. W. Hanson and A. B. Yakovlev, Operator Theory for Electromagnetics – An Introduction (Springer, 2002). [PubMed]
Ch. Hafner, The Generalized Multipole Technique for Computational Electromagnetics (Artech House Books, Boston, 1990).
J. Jin, The Finite Element Method in Electromagnetics (John Wiley, New York, 2002).
P. Monk, Finite Element Methods for Maxwell’s Equations (Oxford University Press, Oxford, 2003). [CrossRef] [PubMed]
J. L. Volakis, A. Chatterjee, and L. C. Kempel, Finite Element Method for Electromagnetics – Antennas, Microwave Circuits and Scattering Applications (IEEE Press, New York, 1998).
C. Bohren and D. Huffmann, Absorption and Scattering of Light by Small Particles (John Wiley, New York, 1983).
S. Ramo, J. R. Whinnery, and T. V. Duzer, Fields and Waves in Communication Electronics (John Wiley, New York, 1993).
R. Geus, “The Jacobi–Davidson algorithm for solving large sparse symmetric eigenvalue problems.” PhD Thesis No. 14734, ETH Zurich 2002.
Z. Bai, J. Demmel, J. Dongarra, A. Ruhe, and H. van der Vorst, Templates for the Solution of Algebraic Eigenvalue Problems: A Practical Guide (SIAM, Philadelphia PA, 2000). [CrossRef]
Trilinos Project Home Page, http://trilinos.sandia.gov/ .
Paraview Home Page, http://www.paraview.org/ .
Home Page of the Swiss National Supercomputing Centre (CSCS), http://www.cscs.ch/ .

Cited By

Alert me when this paper is cited

OSA is able to provide readers links to articles that cite this paper by participating in CrossRef's Cited-By Linking service. CrossRef includes content from more than 3000 publishers and societies. In addition to listing OSA journal articles that cite this paper, citing articles from other participating publishers will also be listed.

Click here to see a list of articles that cite this paper