Understanding localized surface plasmon resonance with propagative surface plasmon polaritons in optical nanogap antennas

Hongwei Jia; Fan Yang; Ying Zhong; Haitao Liu

doi:doi:10.1364/PRJ.4.000293

Photonics Research, 2016, 4 (6): 06000293, Published Online: Dec. 23, 2016

Understanding localized surface plasmon resonance with propagative surface plasmon polaritons in optical nanogap antennas Download： 1189次

论文大纲

Hongwei Jia ^1,2Fan Yang ³Ying Zhong ⁴Haitao Liu ^1,*

Author Affiliations

¹ Key Laboratory of Optical Information Science and Technology, Ministry of Education, Institute of Modern Optics, Nankai University, Tianjin 300071, China

² SZU-NUS Collaborative Innovation Center for Optoelectronic Science & Technology, Key Laboratory of Optoelectronic Devices and Systems of Ministry of Education and Guangdong Province, College of Optoelectronic Engineering, Shenzhen University, Shenzhen 518060, China

³ State Key Laboratory of Low Dimensional Quantum Physics, Department of Physics, Tsinghua University, Beijing 100084, China

⁴ State Key Laboratory of Precision Measuring Technology and Instruments, Tianjin University, Tianjin 300072, China

Surface plasmons Surface plasmons Resonance Resonance Subwavelength structures Subwavelength structures

Abstract

The plasmonic nanogap antenna is an efficient radiating or receiving optical device. The resonance behavior of optical antennas is commonly attributed to the excitation of a localized surface plasmon resonance (LSPR), which can be theoretically defined as the quasi-normal mode (QNM). To clarify the physical origin of the LSPR, we build up an analytical model of the LSPR by considering a multiple scattering process of propagative surface plasmon polaritons (SPPs) on the antenna arms. The model can comprehensively reproduce the complex eigenfrequency and the field distribution of QNMs of the antenna, unveiling that the LSPR arises from a Fabry–Perot resonance of SPPs. By further applying the complex pole expansion theorem of meromorphic functions, the field of the antenna under illumination by a nearby dipole emitter can be analytically expanded with QNMs, which well predicts the frequency response of the enhancement factor of radiation. The present model establishes explicit relations between the concepts of the LSPR and the propagative SPP and integrates the advantages of the Fabry–Perot and QNM formalisms of nanogap antennas.

1. INTRODUCTION

Resonant optical antennas are important devices that can efficiently accelerate the radiation of emitters such as molecules or quantum dots placed nearby ^{[1–6" target="_self" style="display: inline;">–6]} or, reciprocally, allow giant enhancement of the near field under far-field illumination ^{[7–12" target="_self" style="display: inline;">–12]}. By controlling the geometrical parameters to fulfill resonance conditions, optical antennas can be widely used in nonlinear optics ^[13,14], white light supercontinuum generation ^[7], single-emitter fluorescence enhancement ^[1,3], and surface enhanced Raman scattering ^{[1517" target="_self" style="display: inline;">–17]}. The enhancement of field or emission is especially stringent in nanogaps formed at the mouth of metal nanoparticles that are nearly in contact, so that nanogap antennas have many established applications ranging from electron tunneling microscopy, nanocatalysis, and Raman spectroscopy to disruptive electronics and light emission ^[10]. The resonance enhancement is commonly attributed to an excitation of localized surface plasmon resonance (LSPR) as the illumination frequency matches the eigenfrequency of the LSPR ^{[12,18–21" target="_self" style="display: inline;">–21]}.

For the nanogap antenna composed of two strongly coupled individual nanowires, the LSPR is classified into bonding and antibonding modes with opposite field symmetries ^[18,19]. Theoretically, the LSPR can be defined as the quasi-normal mode (QNM) by treating the metallic nanoantenna as an open cavity ^[21]. The QNM is a solution of source-free Maxwell’s equations with an eigenfrequency of complex value due to the cavity loss ^[22]. The formalism of QNMs is a powerful tool able to provide an analytical description of the frequency response of open or lossy resonators ^{[18–23" target="_self" style="display: inline;">–23]} and shows great advantages over standard full-wave solvers that repeat the entire calculation for different frequencies ^[12]. However, presently the QNMs of antennas are commonly obtained via full-wave numerical calculations ^{[20,21,24–29" target="_self" style="display: inline;">–29]}. Too few reports provide a comprehensive analytical description of QNMs (for instance, for some special cavities such as one-dimensional cavities and microspheres ^[30,31], the eigenfrequency and field of the QNM are analytical), which blocks a further understanding of the physical origin of QNMs and thus an efficient design of antenna devices.

To achieve a physical understanding of the resonance behaviors of nanoantennas, different theoretical approaches have been reported. In the lumped-elements circuit model, the antenna is treated as a circuit composed of resistors, inductors, and capacitors, so that the resonance is intuitively understood as the resonance of the AC current in the circuit ^[32,33]. The temporal coupled-mode theory is established by considering the coupling between the eigenmodes and the input and output ports of the resonator, and it can be readily used for analyzing the Fano or Lorentzian line shape of the frequency response of antenna arrays ^[34,35]. However, in these theories, surface plasmon polaritons (SPPs) do not appear in an explicit form. The SPP can be defined as a waveguide mode that propagates along a metal–dielectric interface ^{[36–39" target="_self" style="display: inline;">–39]} and is conceptually different from the LSPR. A counterexample can be found for the analysis of metallic-sphere dimmers by showing how SPPs propagate at the metal surface and are slowed down at the dimmer mouth to build a hot spot ^[40]. SPP Fabry–Perot models have been developed for antennas composed of a single wire, in which SPPs are bouncing back and forth between the arm facets that act as mirrors ^{[41–44" target="_self" style="display: inline;">–44]}. Recently, the Fabry–Perot model has been extended to nanogap antennas with two arms by including the SPP hopping from one arm to the other and the SPP reflection at the gap ^[5,11]. The extended model, which relies on the ab initio calculation of a few elementary SPP scattering coefficients, has been shown to comprehensively and accurately reproduce many antenna characteristics, for example, scattering cross section, near-field enhancement factor, Purcell factor, and far-field radiation pattern. However, since the excitation and scattering processes of SPPs depend on the illumination frequency in an implicit way, the Fabry–Perot model cannot provide an analytical description of the frequency response of the antenna.

In this paper, we combine the Fabry–Perot and QNM formalisms, capitalizing on the intuitive force of the former and the frequency-analyticity of the latter, to propose an insightful and effective theoretical treatment of nanogap antennas. For that purpose, we build up an SPP Fabry–Perot model for the QNM of nanogap antennas by considering a multiple scattering process of SPPs on the antenna arms. The existence of slightly damped QNMs that cause the resonance nature of nanogap antennas is demonstrated with the model by seeking the solutions of two transcendental equations. The model sets a solid electromagnetic foundation for the intuitive picture that the LSPR (i.e., QNM) of the antenna actually arises from the Fabry–Perot resonance of SPPs at the complex eigenfrequency. With a few assumptions on the frequency dependence of the field and with the use of complex pole expansion of meromorphic functions (Mittag–Leffler theorem), the field of the nanogap antenna excited by a nearby point emitter can be expanded upon the basis of QNMs, thus providing an analytical description of the frequency response of the field. For calculating the Purcell factor, the model provides a new analytical expression of the mode volume of QNM in terms of the SPP scattering coefficients, which avoids the difficulty in calculating the mode volume caused by the divergence of the QNM field at infinity that was overcome recently ^[20]. The present model establishes explicit relations between the concepts of the LSPR and the propagative SPP, and it integrates the advantages of the Fabry–Perot and QNM formalisms in understanding the resonance of nanogap antennas, which may inspire new design strategies of antenna devices with different geometries ^{[1,45–50" target="_self" style="display: inline;">–50]}.

2. SPP MODEL OF LSPR FOR NANOGAP ANTENNAS

As shown in Fig. 1(a), the considered nanogap antenna is formed by two arms of gold nanowires (with length $L$ ) separated by a nanogap (with width $w = 0.03 μm$ ). The nanowires have a square cross section with a side length $D = 0.04 μm$ . The antenna is surrounded by air (refractive index $n_{a} = 1$ ) without a substrate for simplicity. For solving QNMs at complex frequencies, the gold permittivity at complex frequencies is obtained through an analytical continuation by using an analytical expression obtained with a polynomial fitting of the experimental data of the gold permittivity at real frequencies ^[51]. Next we will try to derive analytical expressions for the field and the eigenfrequency of the QNMs based on a multiple scattering formalism of SPPs ^[5]. Since the transversal size of the antenna arms is much smaller than the wavelength, only one fundamental SPP mode is propagative and bounded on the antenna arms that are treated as $x$ -invariant waveguides ^[36,37]. All other waveguide modes on antenna arms are either evanescent or unbounded and thus are neglected in the model. In view of the multiple scattering processes of SPPs as sketched in Fig. 1(a), the field of QNM of the antenna can be expressed as where $Ψ = [E, H]$ denotes both the electric- and the magnetic-field vectors. $Ψ_{SPP, +}^{G, s}$ [Fig. 1(b)] and $Ψ_{SPP, +}^{R, s}$ [Fig. 1(c)] denote the fields scattered at the nanogap and at the right termination of the antenna for an incident right-going SPP, respectively. $Ψ_{SPP, -}^{G, s}$ and $Ψ_{SPP, -}^{L, s}$ are defined similarly for an incident left-going SPP. $u = \exp (i k_{0, c} n_{eff} L)$ is the phase shift of the SPP accumulated over one antenna arm ( $k_{0, c} = ω_{c} / c$ , $ω_{c}$ being the complex angular frequency of the QNM, $c$ being the light speed in vacuum, and $n_{eff}$ being the complex effective index of the SPP). Equation (1) is written with the view that the field of the QNM contains four parts: the two fields scattered at the two terminations for incident SPPs propagating away from the gap (with coefficients $a_{1}$ and $a_{2}$ and a damping $u$ ) and the two fields scattered at the nanogap for incident SPPs propagating toward the gap (with coefficients $b_{1}$ and $b_{2}$ and a damping $u$ ). The scattered fields ( $Ψ_{SPP, +}^{G, s}$ , $Ψ_{SPP, +}^{R, s}$ , $Ψ_{SPP, -}^{G, s}$ , and $Ψ_{SPP, -}^{L, s}$ ) and $n_{eff}$ can be calculated independently with the full-wave aperiodic Fourier modal method (a-FMM) ^[5,52,53]. The efficiency and accuracy of the a-FMM compared with other full-wave solvers of Maxwell’s equations such as the finite element method or finite difference time domain method has been confirmed ^[54]. In Eq. (1), the unknown amplitude coefficients $a_{1}$ , $a_{2}$ , $b_{1}$ , and $b_{2}$ of SPPs [Fig. 1(a)] can be determined by solving a set of coupled SPP equations: where $ρ$ and $τ$ denote the SPP reflection and transmission coefficients at the nanogap ^[38] [Fig. 1(b)], and $r$ is the SPP reflection coefficient at the antenna termination ^[39] [Fig. 1(c)]. Here $ρ$ , $τ$ , and $r$ are calculated in a rigorous manner as the scattering matrix elements ^[5,55] with the full-wave a-FMM ^[52,53] without relying on any assumption or any fitting of rigorous data, which ensures a solid electromagnetic foundation of the present model. Equations (2) can be understood intuitively. For the first equation, for instance, the coefficient $a_{1}$ results from two contributions: the first contribution from the reflection ( $ρ$ ) of an incident right-going SPP (with coefficient $b_{1}$ and damping $u$ ) on the left arm, and the second contribution from the transmission ( $τ$ ) of an incident left-going SPP on the right arm (with coefficient $b_{2}$ and damping $u$ ). Note that no excitation terms appear in Eqs. (2), since the QNM is the eigensolution of Maxwell’s equations without source. Equations (2) form a set of homogeneous linear equations of $a_{1}$ , $a_{2}$ , $b_{1}$ , and $b_{2}$ . To ensure the existence of nontrivial solutions of Eqs. (2), the determinant of the coefficient matrix should be zero, which yields Equations (3) form two transcendental equations for determining the complex eigenfrequencies of QNMs. Since $ρ$ , $τ$ , and $r$ are independent of the antenna arm length $L$ and $u = \exp (i k_{0, c} n_{eff} L)$ depends on $L$ analytically, Eqs. (3) could be more computationally efficient in solving the QNM eigenfrequencies than full-wave calculations, especially for large values of $L$ . With Eqs. (3a) or (3b) inserted into Eqs. (2), we can obtain the SPP coefficients $a_{1} = a_{2}$ , $b_{1} = b_{2}$ , or $a_{1} = - a_{2}$ , $b_{1} = - b_{2}$ , respectively. The two solutions correspond to QNMs with different symmetries of field, termed as the bonding mode and antibonding mode, respectively ^[18,19]. With the two solutions inserted into Eq. (1), we finally obtain the analytical expressions for the field distribution of the bonding and antibonding QNMs (with normalization $a_{1} = u^{- 1}$ ): Equations (4) can be further simplified at some specific positions such as in the nanogap or near the antenna terminations, and details can be found in Appendix D. To seek the solution, Eqs. (2) can be rewritten as where $M$ and $N$ are integers that correspond to different orders of QNMs, and $\arg ()$ denotes the argument. In view of the propagative nature of the SPP mode, $n_{eff}$ is approximately a real number, so that $Im (k_{0, c}) \approx \ln (| r | | ρ \pm τ |) / (2 n_{eff} L)$ . It has been shown that $| r | \approx 1$ for the SPP reflection at the antenna termination ^[39], and $| ρ \pm τ | \approx 1$ for the SPP transmission and reflection at the nanogap due to energy conservation ^[38]. Thus we have $Im (k_{0, c}) \approx 0$ , implying the existence of slightly damped QNMs that cause the resonance nature of the nanogap antenna. In fact, a crude evaluation of $k_{0, c}$ can be obtained by calculating the right side of Eqs. (5) [denoted by $g (k_{0, c})$ ] at a certain frequency, in view that $ρ$ , $τ$ , $r$ , $n_{eff}$ , and thus $g (k_{0, c})$ are slowly varying functions of $k_{0, c}$ . The accuracy of this evaluation can be further improved with the iteration formula $k_{0, c}^{m} = g (k_{0, c}^{m - 1})$ ^[11]. To seek the solution of Eqs. (5) in a rigorous manner, we can solve the nonlinear equation $f (k_{0, c}) = k_{0, c} - g (k_{0, c}) = 0$ with classical numerical techniques such as the linear interpolation method ^[56]. By setting two initial values $k_{0, c}^{0}$ and $k_{0, c}^{1}$ and approximating $f (k_{0, c})$ as a linear form $f (k_{0, c}) = p + q k_{0, c}$ , a sequence of $k_{0, c}^{m}$ can be obtained by solving $p + q k_{0, c}^{m} = 0$ , which approaches the exact solution rapidly as $m$ increases. Equations (5) are superior to Eqs. (3) for seeking solutions since for designated integers $M$ or $N$ , the solution of Eqs. (5) is commonly unique and thus not sensitive to the initial values $k_{0, c}^{0}$ and $k_{0, c}^{1}$ . While Eqs. (3) have an infinite number of solutions, so that the solved eigenfrequency is the one closest to the initial values. Taking the real part of both sides of Eqs. (5) after multiplying by $2 n_{eff} L$ , we obtain Equations (6) show that at the complex eigenfrequency $ω_{c}$ of QNMs, the phase shift of SPPs that propagate back and forth over one round on the antenna is multiples of $2 π$ , indicating that the QNM (i.e., LSPR) originates from a Fabry–Perot resonance of SPPs at the eigenfrequency.

(a) Sketch of the nanogap antenna. The antenna is composed of two gold nanowire arms of length L (with a square cross section of side length D=40 nm) separated by a nanogap (gap width w=30 nm). a1, a2, b1, and b2 are the complex amplitude coefficients of SPPs at the complex eigenfrequency of QNMs. (b) and (c) definitions of SPP scattering coefficients ρ, τ, and r and scattered fields ΨSPP,+G,s and ΨSPP,+R,s used in the model.

Fig. 1. (a) Sketch of the nanogap antenna. The antenna is composed of two gold nanowire arms of length $L$ (with a square cross section of side length $D = 40 nm$ ) separated by a nanogap (gap width $w = 30 nm$ ). $a_{1}$ , $a_{2}$ , $b_{1}$ , and $b_{2}$ are the complex amplitude coefficients of SPPs at the complex eigenfrequency of QNMs. (b) and (c) definitions of SPP scattering coefficients $ρ$ , $τ$ , and $r$ and scattered fields $Ψ_{SPP, +}^{G, s}$ and $Ψ_{SPP, +}^{R, s}$ used in the model.

下载图片查看所有图片

Now we check the validity of the model through a comparison to the full-wave numerical data. We first reproduce the eigenfrequency $ω_{c}$ of the QNM of the antenna by solving Eqs. (5), as shown in Table 1 for antenna lengths $L = 0.2$ and 0.6 μm. For solving the QNM with the full-wave a-FMM, a point emitter with complex frequency is placed in the vicinity of the antenna. In view of the divergence of the excited field $F (ω)$ as the emitter frequency $ω$ approaches the eigenfrequency of QNMs, the QNMs can be obtained by solving a nonlinear equation $1 / F (ω) = 0$ (with the linear interpolation method) ^[11,56]. The solution of $1 / F (ω) = 0$ is commonly not unique and depends sensitively on the initial value, which causes difficulty in finding all QNMs within the wavelength range of interest. However, as mentioned earlier, this difficulty can be largely overcome by solving Eqs. (5) of the model due to their uniqueness of solution. It is seen that the eigenfrequencies obtained with the model agree well with the full-wave a-FMM results. Only for the $M = 1$ order QNM does the predicted eigenfrequency slightly deviate from the a-FMM result (especially the imaginary part), which is due to the impact of surface waves other than SPPs that are neglected in the model (see Appendix A for details).

Table 1. Complex Eigenfrequencies $ω_{c}$ ( $λ_{c} = 2 π c / ω_{c}$ ) of QNMs

$L$ (μm)	$M$ or $N$	$λ_{c}$ (SPP model)	$λ_{c}$ (a−FMM)
0.2	$M = 1$	$0.9776 + 0.0574 i$	$0.9635 + 0.0925 i$
0.2	$N = 1$	$0.8068 + 0.0226 i$	$0.8089 + 0.0233 i$
0.6	$M = 1$	$2.4262 + 0.1789 i$	$2.3928 + 0.2823 i$
	$N = 1$	$1.9468 + 0.0921 i$	$1.9525 + 0.0976 i$
	$M = 2$	$1.0830 + 0.0427 i$	$1.0920 + 0.0455 i$
	$N = 2$	$0.9782 + 0.0260 i$	$0.9783 + 0.0251 i$
	$M = 3$	$0.7707 + 0.0166 i$	$0.7716 + 0.0159 i$
	$N = 3$	$0.7410 + 0.0141 i$	$0.7410 + 0.0141 i$
	$M = 4$	$0.6638 + 0.0115 i$	$0.6636 + 0.0117 i$
	$N = 4$	$0.6522 + 0.0103 i$	$0.6522 + 0.0103 i$

查看所有表

Then we will reproduce the field of QNMs with Eqs. (4) of the model. The results for QNMs of different orders with antenna length $L = 0.6 μm$ are shown in Fig. 2 (obtained in plane $y = 0$ ). Note that the fields of QNMs calculated with the a-FMM are normalized by setting the left-going SPP coefficient on the left arm to be $u^{- 1}$ , which is consistent with the QNM normalization in Eqs. (4) of the model. The QNM field distributions predicted by the model (right column) agree well with the full-wave results (left column). This is true even for the $M = 1$ order QNM, for which the model cannot predict the eigenfrequency as accurate (Table 1), and the reason lies in the consistent normalization of the QNM field for the model and for the a-FMM. In Figs. 2(a)–2(f), the first two parts show the distributions of $Re (E_{x})$ for $M = 1$ , 2, 3 and $N = 1$ , 2, 3 orders of QNMs, exhibiting different field symmetries of the bonding and antibonding modes as indicated earlier with Eqs. (4). Other parts of Fig. 2 show the distributions of the electric field amplitude ( $| E | = \sqrt{{| E_{x} |}^{2} + {| E_{y} |}^{2} + {| E_{z} |}^{2}}$ ). It is seen that for the bonding modes ( $M = 1$ , 2, …), the field in the nanogap is very strong, but the field in the gap is nearly zero for the antibonding modes ( $N = 1$ , 2, …). The number of field nodes on the antenna arms increases with the increase of the resonance order $M$ or $N$ , which is due to an interference of the two counterpropagating SPPs on each arm of the antenna. The field of QNMs for antenna length $L = 0.2 μm$ is provided in Appendix B.

Field distributions of QNMs for antenna arm length L=0.6 μm. The left and right columns show the results obtained with the full-wave a-FMM and with the SPP model, respectively. (a)–(c) show the bonding QNMs for M=1, 2, 3, respectively, where the electric-field amplitude is defined as |E|=|Ex|2+|Ey|2+|Ez|2. (d)–(f) show the antibonding QNMs for N=1, 2, 3, respectively.

Fig. 2. Field distributions of QNMs for antenna arm length $L = 0.6 μm$ . The left and right columns show the results obtained with the full-wave a-FMM and with the SPP model, respectively. (a)–(c) show the bonding QNMs for $M = 1$ , 2, 3, respectively, where the electric-field amplitude is defined as $| E | = \sqrt{{| E_{x} |}^{2} + {| E_{y} |}^{2} + {| E_{z} |}^{2}}$ . (d)–(f) show the antibonding QNMs for $N = 1$ , 2, 3, respectively.

下载图片查看所有图片

(a) Definitions of the amplitude coefficients (c1, c2, d1, and d2) of SPPs excited by an x-polarized electric point source located at the center of the nanogap. (b) Definitions of the SPP excitation coefficient β and scattered field used in the model. (c) Enhancement factor F of radiation plotted as a function of illumination frequency ω/(2πc). Results are obtained with the full-wave a-FMM (blue circles) and the SPP model (red solid curve) for antenna length L=0.6 μm.

Fig. 3. (a) Definitions of the amplitude coefficients ( $c_{1}$ , $c_{2}$ , $d_{1}$ , and $d_{2}$ ) of SPPs excited by an $x$ -polarized electric point source located at the center of the nanogap. (b) Definitions of the SPP excitation coefficient $β$ and scattered field used in the model. (c) Enhancement factor $F$ of radiation plotted as a function of illumination frequency $ω / (2 π c)$ . Results are obtained with the full-wave a-FMM (blue circles) and the SPP model (red solid curve) for antenna length $L = 0.6 μm$ .

下载图片查看所有图片

(a) Definitions of the amplitude coefficients (e1, e2, f1, and f2) of SPPs excited by an x-polarized electric point source located near the antenna termination (with a distance of 15 nm). (b) Definitions of the SPP excitation coefficient α and scattered field ΨSourceT used in the model. (c) Enhancement factor F of radiation plotted as a function of illumination frequency ω/(2πc). Results are obtained with the full-wave a-FMM (blue circles) and the SPP model (red-solid curve), with antenna length L=0.6 μm.

Fig. 4. (a) Definitions of the amplitude coefficients ( $e_{1}$ , $e_{2}$ , $f_{1}$ , and $f_{2}$ ) of SPPs excited by an $x$ -polarized electric point source located near the antenna termination (with a distance of 15 nm). (b) Definitions of the SPP excitation coefficient $α$ and scattered field $Ψ_{Source}^{T}$ used in the model. (c) Enhancement factor $F$ of radiation plotted as a function of illumination frequency $ω / (2 π c)$ . Results are obtained with the full-wave a-FMM (blue circles) and the SPP model (red-solid curve), with antenna length $L = 0.6 μm$ .

下载图片查看所有图片

SPP field and the residual field on the surface of the antenna arms for different orders of QNMs. (a1) and (a2) correspond to M=1 and M=2 orders of bonding QNMs, respectively. (b1) and (b2) correspond to N=1 and N=2 orders of antibonding QNMs, respectively. The results are obtained for antenna length L=0.6 μm.

Fig. 5. SPP field and the residual field on the surface of the antenna arms for different orders of QNMs. (a1) and (a2) correspond to $M = 1$ and $M = 2$ orders of bonding QNMs, respectively. (b1) and (b2) correspond to $N = 1$ and $N = 2$ orders of antibonding QNMs, respectively. The results are obtained for antenna length $L = 0.6 μm$ .

下载图片查看所有图片

SPP field and the residual field on the surface of the antenna arms at different resonance peaks of the enhancement factor F of radiation plotted in Figs. 3(c) and 4(c). (a1) and (a2) show the results at resonances corresponding to the M=1 and M=2 orders of QNMs for the case that the source is located at the center of the nanogap [see Fig. 3(c)]. (b1), (b2), (c1), and (c2) show the results at the resonances corresponding to M=1, 2 and N=1, 2 orders of QNMs for the case that the source is located near the antenna termination [see Fig. 4(c)]. The results are obtained for antenna length L=0.6 μm.

Fig. 6. SPP field and the residual field on the surface of the antenna arms at different resonance peaks of the enhancement factor $F$ of radiation plotted in Figs. 3(c) and 4(c). (a1) and (a2) show the results at resonances corresponding to the $M = 1$ and $M = 2$ orders of QNMs for the case that the source is located at the center of the nanogap [see Fig. 3(c)]. (b1), (b2), (c1), and (c2) show the results at the resonances corresponding to $M = 1$ , 2 and $N = 1$ , 2 orders of QNMs for the case that the source is located near the antenna termination [see Fig. 4(c)]. The results are obtained for antenna length $L = 0.6 μm$ .

下载图片查看所有图片

Field distributions of QNMs for antenna length L=0.2 μm. The left and right columns show the results obtained with the full-wave a-FMM and with the SPP model, respectively. (a) and (b) show the real part of Ex component and the amplitude of the electric field (|E|=|Ex|2+|Ey|2+|Ez|2) for the M=1 order QNM. (c) and (d) show the results for the N=1 order QNM.

Fig. 7. Field distributions of QNMs for antenna length $L = 0.2 μm$ . The left and right columns show the results obtained with the full-wave a-FMM and with the SPP model, respectively. (a) and (b) show the real part of $E_{x}$ component and the amplitude of the electric field ( $| E | = \sqrt{{| E_{x} |}^{2} + {| E_{y} |}^{2} + {| E_{z} |}^{2}}$ ) for the $M = 1$ order QNM. (c) and (d) show the results for the $N = 1$ order QNM.

下载图片查看所有图片

3. REPRODUCING ENHANCEMENT FACTOR OF RADIATION WITH THE MODEL

In previous theories that treat the antenna arms as microcavities of SPPs ^[5,41,43], the enhancement of the radiation from a dipole emitter in the vicinity of the antenna is attributed to a Fabry–Perot resonance of SPPs. Another interpretation of the radiation enhancement is based on a resonant excitation of the LSPR (i.e., QNM) ^{[18–21" target="_self" style="display: inline;">–21]}. In this section, we will apply the SPP model of the LSPR established in Section 2 to reproduce the enhancement factor of the radiation of the nanogap antenna, thus establishing explicit relations between the two interpretations.

3.2 A. Radiation Enhancement for an Emitter Located in the Nanogap

We first consider the case that an $x$ -polarized electric dipole source (with a current density expressed with Dirac function $J = δ (x, y, z) x$ , $x$ being a unit vector along the $x$ direction) is located at the center of the gap of the antenna, as sketched in Fig. 3(a). The dipole source is used to represent a molecule ^{[13" target="_self" style="display: inline;">–3]} or a quantum dot ^[4,6] in practical situations. The enhancement factor $F = Γ_{tot} / Γ_{air}$ is used to characterize the acceleration of the spontaneous emission/decay rate of a dipole emitter in the vicinity of nanostructures ^[57], where $Γ_{tot}$ and $Γ_{air}$ represent the spontaneous emission rates with and without the presence of nanostructures, respectively. For the present case, $Γ_{tot}$ can be calculated with the field at the source position, $Γ_{tot} = - Re [E_{x} (r_{0})] / 2$ ( $r_{0}$ being the coordinate of the source), and $Γ_{air} = η_{vac} k_{0}^{2} n_{a} / (12 π)$ ( $k_{0} = ω / c$ , $n_{a} = 1$ being the refractive index of air, and $η_{vac}$ being the wave impedance in vacuum). The Purcell factor, which denotes the highest acceleration of the spontaneous emission/decay rate when changing the frequency, the position, and the polarization of the emitter, is required to achieve high values for applications such as high-speed quantum information processing ^[58] and fluorescence sensing of molecules ^[59]. To obtain the field under the excitation of the dipole emitter, we use $c_{1}$ , $c_{2}$ , $d_{1}$ , and $d_{2}$ to denote the amplitude coefficients of SPPs that propagate back and forth on the antenna arms, as sketched in Fig. 3(a). To solve the unknown SPP coefficients, a set of coupled-mode equations can be written ^[5]: $c_{1} = β + d_{1} u ρ + d_{2} u τ,$ (7a) $c_{2} = β + d_{2} u ρ + d_{1} u τ,$ (7b) $d_{1} = c_{1} u r,$ (7c) $d_{2} = c_{2} u r,$ (7d)

where $β$ denotes the excitation coefficient of SPPs by the dipole source [Fig. 3(b)], which can be calculated as the scattering-matrix element ^[5,53,55] or with the Lorentzian reciprocity relation by sending an incident SPP [see Eq. (C5) in Appendix C]. Other scattering coefficients ( $ρ$ , $τ$ , and $r$ ) have been defined earlier in Fig. 1. Equations (7) can be understood in a way similar to Eqs. (2). They can be analytically solved as $c_{1} = c_{2} = \frac{β}{1 - u^{2} r (ρ + τ)},$ (8a) $d_{1} = d_{2} = \frac{β u r}{1 - u^{2} r (ρ + τ)} .$ (8b)Then the field excited by the emitter can be expressed as ^[5] $Ψ = Ψ_{Source}^{G} + c_{1} u Ψ_{SPP, -}^{L, s} + c_{2} u Ψ_{SPP, +}^{R, s} + d_{1} u Ψ_{SPP, +}^{G, s} + d_{2} u Ψ_{SPP, -}^{G, s},$ (9)where $Ψ_{Source}^{G}$ denotes the field excited by the point source for an infinitely long nanowire with a nanogap [see Fig. 3(b)], and other terms ( $Ψ_{SPP +}^{G, s}$ , $Ψ_{SPP -}^{G, s}$ , $Ψ_{SPP +}^{R, s}$ and $Ψ_{SPP -}^{L, s}$ ) are defined earlier in Fig. 1. Equation (9) is written similar to Eq. (1), except for the first term representing the source excitation. Note that all the quantities on the right side of Eq. (9) depend on the illumination frequency (taking real values) in an implicit form. To seek the frequency response in an analytical form, next we consider to expand the excited field expressed with Eq. (9) upon the basis of QNMs ^[20,60,61]. Equations (3) and (8) show that as the illumination frequency $ω$ approaches the complex eigenfrequency $ω_{c}$ of QNMs, the SPP coefficients $c_{1}$ , $c_{2}$ , $d_{1}$ , and $d_{2}$ approach infinity, and resultantly, the excited field $Ψ$ in Eq. (9) diverges to infinity. Then we assume that $Ψ$ is a meromorphic function of $ω$ , with the $ω_{c}$ being the first-order pole of $Ψ$ , and that $Ψ$ is bounded as $| ω | \to \infty$ in the complex plane. This assumption has been shown to be valid for the field inside or near nanogrooves ^[62] and will be shown to be valid for the present case of nanogap antennas. According to the complex pole expansion of meromorphic functions (Mittag–Leffler theorem) ^{[6264" target="_self" style="display: inline;">–64]}, the field $Ψ$ can be expressed analytically with respect to $ω$ : $Ψ (ω) = Ψ (0) + \sum_{n} \frac{ω / ω_{c, n}}{ω - ω_{c, n}} p_{n},$ (10)where $ω_{c, n}$ denotes the eigenfrequency of the $n$ th order QNM [ $n$ can be $M$ or $N$ in Eqs. (5)]. The term $Ψ (0) \approx 0$ in view that at $ω = 0$ the scattering effect of the finite-size antenna vanishes for an infinitely-large wavelength. $p_{n} = \lim_{ω \to ω_{c, n}} (ω - ω_{c, n}) Ψ (ω)$ is the residue of $Ψ$ at the pole $ω_{c, n}$ . Note that the way of deriving the expansion of Eq. (10) is essentially different from previous methods (such as the Lorentzian reciprocity theorem ^[20] and the two component expansion method ^[60,61]) that depend on the assumption of the completeness of QNMs. With Eqs. (8) and (9), we obtain $p_{n} = \lim_{ω \to ω_{c, n}} \frac{β u (ω - ω_{c, n})}{1 - u^{2} r (ρ + τ)} \times (Ψ_{SPP, -}^{L, s} + Ψ_{SPP, +}^{R, s} + u r Ψ_{SPP, +}^{G, s} + u r Ψ_{SPP, -}^{G, s}),$ (11)where $\lim_{ω \to ω_{c, n}} (ω - ω_{c, n}) Ψ_{Source}^{G} (ω)$ is zero since $Ψ_{Source}^{G} (ω)$ is a nonresonant term. The summation in the last bracket in Eq. (11) is nothing else than the field of the bonding QNM expressed with Eq. (4a), showing that for the emitter at the center of the nanogap, only the bonding QNMs are excited (i.e., $n = M$ ). The limit in Eq. (11) (with $n = M$ ) is the type of 0/0 and can be calculated with L’Hospital’s rule, $\lim_{ω \to ω_{c, M}} \frac{β u (ω - ω_{c, M})}{1 - u^{2} r (ρ + τ)} = \frac{{(β u)}_{ω = ω_{c, M}}}{{[- \frac{4 π i ω n_{eff} L}{c} (\frac{1}{ω} + \frac{1}{n_{eff}} \frac{\partial n_{eff}}{\partial ω}) - \frac{1}{r} \frac{\partial r}{\partial ω} - \frac{1}{ρ + τ} \frac{\partial (ρ + τ)}{\partial ω}]}_{ω = ω_{c, M}}} .$ (12)With Eqs. (10)–(12), we finally obtain $Ψ (ω) = \sum_{M} η_{M} (ω) Ψ_{QNM, M}^{b},$ (13)where $Ψ_{QNM, M}^{b}$ is the field of the $M$ th-order bonding QNM at eigenfrequency $ω_{c, M}$ [expressed with Eq. (4a)], and $η_{M}$ is a complex expansion coefficient expressed as $η_{M} (ω) = \frac{ω / ω_{c, M}}{ω - ω_{c, M}} \times \frac{{(β u)}_{ω = ω_{c, M}}}{{[- \frac{4 π i ω n_{eff} L}{c} (\frac{1}{ω} + \frac{1}{n_{eff}} \frac{\partial n_{eff}}{\partial ω}) - \frac{1}{r} \frac{\partial r}{\partial ω} - \frac{1}{ρ + τ} \frac{\partial (ρ + τ)}{\partial ω}]}_{ω = ω_{c, M}}} .$ (14)Equation (14) shows that the expansion coefficients of QNMs can be analytically expressed with the SPP scattering coefficients at the eigenfrequency $ω_{c, M}$ . Note that the scattering coefficients ( $β$ , $r$ , $ρ$ , $τ$ ) and the complex effective index ( $n_{eff}$ ) of the SPP are calculated at the eigenfrequency (with the full-wave a-FMM), so that they are independent of the illumination frequency $ω$ . Also note that when changing the position and polarization of the point emitter, in Eqs. (13) and (14) only the SPP excitation coefficient $β$ changes and all other parameters do not change. And with the use of reciprocity theorem, $β$ for different emitter positions and polarizations can be obtained with only one full-wave calculation [see Eq. (C5) in Appendix C]. This virtue with the use of reciprocity is analogous to that with the use of reciprocity between a point emitter and QNMs ^[20]. Equations (13) and (14) provide an analytical description of the frequency response of nanogap antennas along with an intuitive multiple scattering picture of SPPs. The selection of the QNMs in the expansion of Eq. (13) is determined by the frequency range of interests.

For the case that the antenna resonance is dominated by a single QNM, the Purcell factor can be expressed as ^[20] $F_{P} = 3 {(4 π^{2})}^{- 1} {(λ_{0} / n_{a})}^{3} Re (Q / V)$ , where $λ_{0}$ is the vacuum wavelength at resonance, $n_{a}$ is the refractive index of environment ( $n_{a} = 1$ in air), and $Q = - Re (ω_{c}) / [2 Im (ω_{c})]$ and $V$ are the quality factor and the complex-valued mode volume of the dominant QNM, respectively. With Eq. (14), the mode volume $V$ can be expressed analytically with the SPP scattering coefficients (a detailed derivation can be found in Appendix C). This avoids the difficulty caused by the divergence of the QNM field at infinity when calculating the mode volume as an integral of the QNM field over the whole space. This difficulty has been recently overcome by introducing perfectly matched layers to treat the divergence of the QNM field at infinity ^[20].

Now we check the validity of the model in predicting the enhancement factor $F$ of radiation. The results obtained with the full-wave a-FMM and with the model [Eqs. (13) and (14)] are shown in Fig. 3(c) with blue circles and red solid curves, respectively, with antenna length $L = 0.6 μm$ . It is seen that the model predictions agree well with the a-FMM data except for an observable deviation at the resonance corresponding to the $M = 1$ order QNM, which is due to the impact of surface waves other than SPPs (see Appendix A for details). As shown explicitly by Eqs. (13) and (14), the enhancement factor $F$ of emission peaks at $ω = Re (ω_{c, M})$ , which is confirmed by the results in Fig. 3(c). In view of $ω_{c, M} \approx Re (ω_{c, M})$ , since $Im (ω_{c, M})$ is a small number (as shown in Table 1), Eq. (6a) at $ω = ω_{c, M}$ then approximately becomes $2 k_{0} Re (n_{eff}) L + \arg (r) + \arg (ρ + τ) = 2 M π$ at $ω = Re (ω_{c, M})$ . The latter is the already derived phase-matching condition for predicting resonance ^[5] and is related to a resonant excitation (i.e., constructive interference) of SPPs at real frequencies. In this sense, the resonant excitation of QNMs and that of SPPs for explaining the enhanced radiation are logically unified. For design tasks, this phase-matching condition provides an analytical expression to determine the antenna arm length $L$ if the real resonance frequency $Re (ω_{c, M})$ is prescribed.

3.3 B. Radiation Enhancement for an Emitter Located Near the Antenna Termination

As has been shown in Section 3.A, the antibonding modes do not contribute to the radiation of the emitter located in the nanogap. This can be also understood with the reciprocity between the source and the QNM ^[20] in view that the electric field of antibonding QNMs in the nanogap is nearly zero. To excite both the bonding and the antibonding modes, we place an $x$ -polarized point emitter near one termination of the antenna, as shown in Fig. 4(a), where the distance between the emitter and the antenna termination is 15 nm. The unknown complex amplitude coefficients of SPPs are denoted by $e_{1}$ , $e_{2}$ , $f_{1}$ , and $f_{2}$ . They satisfy a set of coupled-mode equations, $f_{1} = e_{1} u r,$ (15a) $f_{2} = α + e_{2} u r,$ (15b) $e_{1} = f_{1} u ρ + f_{2} u τ,$ (15c) $e_{2} = f_{2} u ρ + f_{1} u τ,$ (15d)where $α$ denotes the excitation coefficient of the SPP by the emitter ^[43] [see Fig. 4(b)] and can be calculated as the scattering-matrix element ^[5,55] or with the Lorentzian reciprocity theorem [see Eqs. (C9) in Appendix C]. To solve Eqs. (15), we reformulate the unknown coefficients as $e_{1} + e_{2} = 2 e_{+}$ , $e_{1} - e_{2} = 2 e_{-}$ , $f_{1} + f_{2} = 2 f_{+}$ , and $f_{1} - f_{2} = 2 f_{-}$ , and the analytical solution can be obtained, $e_{+} = \frac{α u (ρ + τ) / 2}{1 - u^{2} r (ρ + τ)},$ (16a) $e_{-} = - \frac{α u (ρ - τ) / 2}{1 - u^{2} r (ρ - τ)},$ (16b) $f_{+} = \frac{α / 2}{1 - u^{2} r (ρ + τ)},$ (16c) $f_{-} = - \frac{α / 2}{1 - u^{2} r (ρ - τ)} .$ (16d)Then the excited field of the antenna can be expressed as $Ψ = Ψ_{Source}^{T} + e_{1} u Ψ_{SPP, -}^{L, s} + e_{2} u Ψ_{SPP, +}^{R, s} + f_{1} u Ψ_{SPP, +}^{G, s} + f_{2} u Ψ_{SPP, -}^{G, s} .$ (17)Equation (17) is written in a way similar to Eq. (1), except for the first term representing the source excitation [as sketched in Fig. 4(b)]. Inserting Eqs. (16) into Eq. (17), we obtain $Ψ = Ψ_{Source}^{T} + \frac{α / (2 r)}{1 - u^{2} r (ρ + τ)} [u^{2} r (ρ + τ) Ψ_{SPP, -}^{L, s} + u^{2} r (ρ + τ) Ψ_{SPP, +}^{R, s} + u r Ψ_{SPP, +}^{G, s} + u r Ψ_{SPP, -}^{G, s}] + \frac{- α / (2 r)}{1 - u^{2} r (ρ - τ)} [u^{2} r (ρ - τ) Ψ_{SPP, -}^{L, s} - u^{2} r (ρ - τ) Ψ_{SPP, +}^{R, s} + u r Ψ_{SPP, +}^{G, s} - u r Ψ_{SPP, -}^{G, s}] .$ (18)Then applying the complex pole expansion of meromorphic functions to Eq. (18) [with the use of Eqs. (3) and (4) for QNMs, as has been done for deriving Eqs. (13) and (14)], we finally obtain $Ψ (ω) = \sum_{M} ξ_{M} (ω) Ψ_{QNM, M}^{b} + \sum_{N} ζ_{N} (ω) Ψ_{QNM, N}^{a},$ (19)where $ξ_{M}$ and $ζ_{N}$ are complex expansion coefficients for the bonding ( $Ψ_{QNM, M}^{b}$ ) and antibonding QNMs ( $Ψ_{QNM, N}^{a}$ ), respectively, and are expressed as $ξ_{M} (ω) = \frac{ω / ω_{c, M}}{ω - ω_{c, M}} \times \frac{{(α / 2 r)}_{ω = ω_{c, M}}}{{[- \frac{4 π i ω n_{eff} L}{c} (\frac{1}{ω} + \frac{1}{n_{eff}} \frac{\partial n_{eff}}{\partial ω}) - \frac{1}{r} \frac{\partial r}{\partial ω} - \frac{1}{ρ + τ} \frac{\partial (ρ + τ)}{\partial ω}]}_{ω = ω_{c, M}}},$ (20a) $ζ_{N} (ω) = \frac{ω / ω_{c, N}}{ω - ω_{c, N}} \times \frac{{(α / 2 r)}_{ω = ω_{c, N}}}{{[\frac{4 π i ω n_{eff} L}{c} (\frac{1}{ω} + \frac{1}{n_{eff}} \frac{\partial n_{eff}}{\partial ω}) + \frac{1}{r} \frac{\partial r}{\partial ω} + \frac{1}{ρ - τ} \frac{\partial (ρ - τ)}{\partial ω}]}_{ω = ω_{c, N}}} .$ (20b)Equation (19) shows that for an emitter located near the termination of the antenna, both the bonding and the antibonding QNMs contribute to the radiation ^[18,19]. Similar to the calculation of $β$ in Eq. (14), with the use of reciprocity, $α$ in Eqs. (20) for different emitter positions and polarizations can be obtained with only one full-wave calculation [see Eq. (C9) in Appendix C].

To check the validity of the SPP model, we plot the enhancement factor $F$ as a function of frequency in Fig. 4(c). The results are obtained with the full-wave a-FMM (blue circles) and with Eqs. (19) and (20) of the model (red solid curves) for antenna length $L = 0.6 μm$ . Good agreement can be seen in the figure, except for a deviation at the resonance corresponding to the $M = 1$ order QNM, which is due to the impact of surface waves other than SPPs (see Appendix A). Similar to the analysis in Section 3.A, Eqs. (19) and (20) show that the enhancement factor $F$ peaks at $ω = Re (ω_{c, M})$ and $ω = Re (ω_{c, N})$ , corresponding to the excitation of bonding and antibonding QNMs, respectively.

4. CONCLUSION

We build up an analytical model for the LSPR (i.e., QNM) of nanogap antennas by considering the multiple-scattering processes of SPPs that propagate back and forth on the antenna arms. The model integrates the advantages of the QNM approach and the Fabry–Perot model of SPPs and logically unifies the interpretations of the resonance behaviors of nanogap antennas with the two approaches. The existence of slightly damped QNMs for nanogap antennas is demonstrated with the model by seeking the solutions of two transcendental equations with the complex frequency as the unknown. The latter can be sometimes more computationally efficient than full-wave solvers of QNMs due to their uniqueness of solutions and analyticity with respect to the antenna arm length. The model demonstrates the fact that the QNM originates from a Fabry–Perot resonance of SPPs at the complex eigenfrequency. The bonding and the antibonding QNMs that possess different symmetries of field are explicitly discriminated with different analytical expressions. The predicted eigenfrequency slightly deviates from the full-wave a-FMM results for the lowest-order bonding QNM, unveiling the impact of surface waves other than SPPs to the antenna resonance.

With the assumption that the field of the antenna is a meromorphic function of frequency and by using the complex pole expansion theorem (Mittag–Leffler theorem), the field of the antenna under excitation by a nearby point emitter can be expanded with QNMs. With the model, the expansion coefficients and the mode volume of QNMs for calculating the Purcell factor are analytically expressed with the scattering coefficients of the SPP mode at the complex eigenfrequency of QNMs. The enhancement factor of radiation is shown to peak when the illumination frequency matches the real part of the complex eigenfrequency of QNMs, which is shown to be equivalent to a phase-matching condition that represents the Fabry–Perot resonance of SPPs. The present analysis of the LSPR with the SPP model establishes explicit relations between the concepts of LSPR and propagative SPP. The method can be extended to the analysis of other antenna structures, such as antennas with circular cross section ^[37,39], $L$ -shaped antennas ^[48], cross antennas ^[45], split ring antennas ^[46,47], or antenna arrays ^[1]. For antennas with arms of varying cross section such as bowtie antennas ^[49] or diabolo antennas ^[50], it could be possible to extend the present model by defining a SPP mode that propagates along the antenna arms adiabatically if the variation of the cross section is slow enough.

5 Acknowledgment

Acknowledgment. Dr. Philippe Lalanne is acknowledged for providing helpful comments.

References

[1] R. M. Bakker, V. P. Drachev, Z. Liu, H. K. Yuan, R. H. Pedersen, A. Boltasseva, J. J. Chen, J. Irudayaraj, A. V. Kildishev, V. M. Shalaev. Nanoantenna array-induced fluorescence enhancement and reduced lifetimes. New J. Phys., 2008, 10: 125022.

[2] O. L. Muskens, V. Giannini, J. A. Sanchez-Gil, J. Gómez Rivas. Strong enhancement of the radiative decay rate of emitters by single plasmonic nanoantennas. Nano Lett., 2007, 7: 2871-2875.

[3] A. Kinkhabwala, Z. Yu, S. Fan, Y. Avlasevich, K. Müllen, W. E. Moerner. Large single-molecule fluorescence enhancements produced by a bowtie nanoantenna. Nat. Photonics, 2009, 3: 654-657.

[4] A. V. Akimov, A. Mukherjee, C. L. Yu, D. E. Chang, A. S. Zibrov, P. R. Hemmer, H. Park, M. D. Lukin. Generation of single optical plasmons in metallic nanowires coupled to quantum dots. Nature, 2007, 450: 402-406.

[5] H. Jia, H. Liu, Y. Zhong. Role of surface plasmon polaritons and other waves in the radiation of resonant optical dipole antennas. Sci. Rep., 2015, 5: 8456.

[6] X. W. Chen, V. Sandoghdar, M. Agio. Coherent interaction of light with a metallic structure coupled to a single quantum emitter: from superabsorption to cloaking. Phys. Rev. Lett., 2013, 110: 153605.

[7] P. Mühlschlegel, H. J. Eisler, O. J. F. Martin, B. Hech, D. W. Pohl. Resonant optical antennas. Science, 2005, 308: 1607-1609.

[8] Z. Liu, E. Li, V. M. Shalaev, A. V. Kildishev. Near field enhancement in silver nanoantenna–superlens systems. Appl. Phys. Lett., 2012, 101: 021109.

[9] H. Fischer, O. J. Martin. Polarization sensitivity of optical resonant dipole antennas. J. Eur. Opt. Soc., 2008, 3: 08018.

[10] L. Novotny, N. Van Hulst. Antennas for light. Nat. Photonics, 2011, 5: 83-90.

[11] H. Jia, P. Lalanne, H. Liu. Comprehensive surface-wave description for the nano-scale energy concentration with resonant dipole antennas. Plasmonics, 2016, 11: 1025-1033.

[12] R. Esteban, G. Aguirregabiria, A. G. Borisov, Y. M. Wang, P. Nordlander, G. W. Bryant, J. Aizpurua. The morphology of narrow gaps modifies the plasmonic response. ACS Photon., 2015, 2: 295-305.

[13] T. Hanke, G. Krauss, D. Träutlein, B. Wild, R. Bratschitsch, A. Leitenstorfer. Efficient nonlinear light emission of single gold optical antennas driven by few-cycle near-infrared pulses. Phys. Rev. Lett., 2009, 103: 257404.

[14] H. Harutyunyan, G. Volpe, R. Quidant, L. Novotny. Enhancing the nonlinear optical response using multifrequency gold-nanowire antennas. Phys. Rev. Lett., 2012, 108: 217403.

[15] W. Zhang, H. Fischer, T. Schmid, R. Zenobi, O. J. Martin. Mode-selective surface-enhanced Raman spectroscopy using nanofabricated plasmonic dipole antennas. J. Phys. Chem. C, 2009, 113: 14672-14675.

[16] F. Jäckel, A. A. Kinkhabwala, W. E. Moerner. Gold bowtie nanoantennas for surface-enhanced Raman scattering under controlled electrochemical potential. Chem. Phys. Lett., 2007, 446: 339-343.

[17] K. Höflich, M. Becker, G. Leuchs, S. Christiansen. Plasmonic dimer antennas for surface enhanced Raman scattering. Nanotechnology, 2012, 23: 185303.

[18] M. Liu, T. W. Lee, S. K. Gray, P. Guyot-Sionnest, M. Pelton. Excitation of dark plasmons in metal nanoparticles by a localized emitter. Phys. Rev. Lett., 2009, 102: 107401.

[19] J. S. Huang, J. Kern, P. Geisler, P. Weinmann, M. Kamp, A. Forchel, P. Biagioni, B. Hecht. Mode imaging and selection in strongly coupled nanoantennas. Nano Lett., 2010, 10: 2105-2110.

[20] C. Sauvan, J. P. Hugonin, I. S. Maksymov, P. Lalanne. Theory of the spontaneous optical emission of nanosize photonic and plasmon resonators. Phys. Rev. Lett., 2013, 110: 237401.

[21] J. Yang, H. Giessen, P. Lalanne. Simple analytical expression for the peak-frequency shifts of plasmonic resonances for sensing. Nano Lett., 2015, 15: 3439-3444.

[22] E. S. C. Ching, P. T. Leung, A. M. van den Brink, W. M. Suen, S. S. Tong, K. Young. Quasinormal-mode expansion for waves in open systems. Rev. Mod. Phys., 1998, 70: 1545-1554.

[23] Q. Bai, M. Perrin, C. Sauvan, J. P. Hugonin, P. Lalanne. Efficient and intuitive method for the analysis of light scattering by a resonant nanostructure. Opt. Express, 2013, 21: 27371-27382.

[24] R. C. Ge, P. T. Kristensen, J. F. Young, S. Hughes. Quasinormal mode approach to modelling light-emission and propagation in nanoplasmonics. New J. Phys., 2014, 16: 113048.

[25] C. Sauvan, J. P. Hugonin, R. Carminati, P. Lalanne. Modal representation of spatial coherence in dissipative and resonant photonic systems. Phys. Rev. A, 2014, 89: 043825.

[26] R. C. Ge, J. F. Young, S. Hughes. Quasi-normal mode approach to the local-field problem in quantum optics. Optica, 2015, 2: 246-249.

[27] P. T. Leung, Y. T. Liu, C. Y. Tam, K. Young. Numerical solution for quasinormal modes for potentials with exponential tails. Phys. Lett. A, 1998, 247: 253-260.

[28] J. Wiersig. Boundary element method for resonances in dielectric microcavities. J. Opt. A, 2003, 5: 53-60.

[29] M. B. Doost, W. Langbein, E. A. Muljarov. Resonant-state expansion applied to three-dimensional open optical systems. Phys. Rev. A, 2014, 90: 013834.

[30] P. T. Leung, S. Y. Liu, K. Young. Completeness and orthogonality of quasinormal modes in leaky optical cavities. Phys. Rev. A, 1994, 49: 3057-3067.

[31] E. A. Muljarov, W. Langbein, R. Zimmermann. Brillouin–Wigner perturbation theory in open electromagnetic systems. Europhys. Lett., 2010, 92: 50010.

[32] M. S. Eggleston, K. Messer, L. Zhang, E. Yablonovitch, M. C. Wu. Optical antenna enhanced spontaneous emission. Proc. Natl. Acad. Sci. USA, 2015, 112: 1704-1709.

[33] A. Alù, N. Engheta. Input impedance, nanocircuit loading, and radiation tuning of optical nanoantennas. Phys. Rev. Lett., 2008, 101: 043901.

[34] A. Artar, A. A. Yanik, H. Altug. Directional double Fano resonances in plasmonic hetero-oligomers. Nano Lett., 2011, 11: 3694-3700.

[35] R. Adato, H. Altug. In-situ ultra-sensitive infrared absorption spectroscopy of biomolecule interactions in real time with plasmonic nanoantennas. Nat. Commun., 2013, 4: 2154.

[36] VassalloC., Optical Waveguide Concepts (Elsevier, 1991).

[37] D. E. Chang, A. S. Sørensen, P. R. Hemmer, M. D. Lukin. Strong coupling of single emitters to surface plasmons. Phys. Rev. B, 2007, 76: 035420.

[38] H. Liu. Coherent-form energy conservation relation for the elastic scattering of a guided mode in a symmetric scattering system. Opt. Express, 2013, 21: 24093-24098.

[39] R. Gordon. Reflection of cylindrical surface waves. Opt. Express, 2009, 17: 18621-18629.

[40] A. Aubry, D. Y. Lei, S. A. Maier, J. B. Pendry. Interaction between plasmonic nanoparticles revisited with transformation optics. Phys. Rev. Lett., 2010, 105: 233901.

[41] T. H. Taminiau, F. D. Stefani, N. F. van Hulst. Optical nanorod antennas modeled as cavities for dipolar emitters: evolution of sub- and super-radiant modes. Nano Lett., 2011, 11: 1020-1024.

[42] L. Douillard, F. Charra, Z. Korczak, R. Bachelot, S. Kostcheev, G. Lerondel, P. M. Adam, P. Royer. Short range plasmon resonators probed by photoemission electron microscopy. Nano Lett., 2008, 8: 935-940.

[43] C. Liu, H. Liu, Y. Zhong. Impact of surface plasmon polaritons and other waves on the radiation of a dipole emitter close to a metallic nanowire antenna. Opt. Express, 2014, 22: 25539-25549.

[44] J. Dorfmuller, R. Vogelgesang, W. Khunsin, C. Rockstuhl, C. Etrichm, K. Kern. Plasmonic nanowire antennas: experiment, simulation, and theory. Nano Lett., 2010, 10: 3596-3603.

[45] P. Biagioni, J. S. Huang, L. Duò, M. Finazzi, B. Hecht. Cross resonant optical antenna. Phys. Rev. Lett., 2009, 102: 256801.

[46] J. Scheuer. Ultra-high enhancement of the field concentration in split ring resonators by azimuthally polarized excitation. Opt. Express, 2011, 19: 25454-25464.

[47] E. Cubukcu, S. Zhang, Y. S. Park, G. Bartal, X. Zhang. Split ring resonator sensors for infrared detection of single molecular monolayers. Appl. Phys. Lett., 2009, 95: 043113.

[48] S. C. Jiang, X. Xiong, Y. S. Hu, Y. H. Hu, G. B. Ma, R. W. Peng, C. Sun, M. Wang. Controlling the polarization state of light with a dispersion-free metastructure. Phys. Rev. X, 2014, 4: 021026.

[49] X. Zhang, C. J. Chung, S. Wang, H. Subbaraman, Z. Pan, Q. Zhan, R. T. Chen. Integrated broadband bowtie antenna on transparent silica substrate. IEEE Antennas Wireless Propag. Lett., 2016, 15: 1377-1381.

[50] Z. Pan, J. Guo. Enhanced optical absorption and electric field resonance in diabolo metal bar optical antennas. Opt. Express, 2013, 21: 32491-32500.

[51] PalikE. D., Handbook of Optical Constants of Solids Part II (Academic, 1985).

[52] J. P. Hugonin, P. Lalanne. Perfectly matched layers as nonlinear coordinate transforms: a generalized formalization. J. Opt. Soc. Am. A, 2005, 22: 1844-1849.

[53] The calculation is performed with in-house software, LiuH., DIF CODE for Modeling Light Diffraction in Nanostructures (Nankai University, 2010).

[54] M. Besbes, J. P. Hugonin, P. Lalanne, S. van Haver, O. T. A. Janssen, A. M. Nugrowati, M. Xu, S. F. Pereira, H. P. Urbach, A. S. van de Nes, P. Bienstman, G. Granet, A. Moreau, S. Helfert, M. Sukharev, T. Seideman, F. I. Baida, B. Guizal, D. Van Labeke. Numerical analysis of a slit-groove diffraction problem. J. Eur. Opt. Soc., 2007, 2: 07022.

[55] L. Li. Formulation and comparison of two recursive matrix algorithms for modeling layered diffraction gratings. J. Opt. Soc. Am. A, 1996, 13: 1024-1035.

[56] Y. Li, H. Liu, H. Jia, F. Bo, G. Zhang, J. Xu. Fully-vectorial modeling of cylindrical microresonators with aperiodic Fourier modal method. J. Opt. Soc. Am. A, 2014, 31: 2459-2466.

[57] E. M. Purcell. Spontaneous emission probabilities at radio frequencies. Phys. Rev., 1946, 69: 37-38.

[58] E. Knill, R. Laflamme, G. J. Milburn. A scheme for efficient quantum computation with linear optics. Nature, 2001, 409: 46-52.

[59] M. Agio. Optical antennas as nanoscale resonators. Nanoscale, 2012, 4: 692-706.

[60] P. T. Leung, S. S. Tong, K. Young. Two-component eigenfunction expansion for open systems described by the wave equation I: completeness of expansion. J. Phys. A, 1997, 30: 2139-2151.

[61] P. T. Leung, S. S. Tong, K. Young. Two-component eigenfunction expansion for open systems described by the wave equation II: linear space structure. J. Phys. A, 1997, 30: 2153-2162.

[62] F. Yang, H. Liu, H. Jia, Y. Zhong. Analytical description of quasi normal mode in resonant plasmonic nano-cavities. J. Opt., 2016, 18: 035003.

[63] ArfkenG. B., Mathematical Methods for Physicists (Elsevier, 2005).

[64] M. B. Doost, W. Langbein, E. A. Muljarov. Resonant state expansion applied to two-dimensional open optical systems. Phys. Rev. A, 2013, 87: 043827.

[65] F. van Beijnum, C. Rétif, C. B. Smiet, H. Liu, P. Lalanne, M. P. van Exter. Quasi-cylindrical wave contribution in experiments on extraordinary optical transmission. Nature, 2012, 492: 411-414.

1. INTRODUCTION

2. SPP MODEL OF LSPR FOR NANOGAP ANTENNAS

3. REPRODUCING ENHANCEMENT FACTOR OF RADIATION WITH THE MODEL

3.2 A. Radiation Enhancement for an Emitter Located in the Nanogap

3.3 B. Radiation Enhancement for an Emitter Located Near the Antenna Termination

4. CONCLUSION

5 Acknowledgment

Hongwei Jia, Fan Yang, Ying Zhong, Haitao Liu. Understanding localized surface plasmon resonance with propagative surface plasmon polaritons in optical nanogap antennas[J]. Photonics Research, 2016, 4(6): 06000293.

Understanding localized surface plasmon resonance with propagative surface plasmon polaritons in optical nanogap antennas Download： 1189次

1. INTRODUCTION

2. SPP MODEL OF LSPR FOR NANOGAP ANTENNAS

Table 1. Complex Eigenfrequencies $ω_{c}$ ( $λ_{c} = 2 π c / ω_{c}$ ) of QNMs

3. REPRODUCING ENHANCEMENT FACTOR OF RADIATION WITH THE MODEL

3.2 A. Radiation Enhancement for an Emitter Located in the Nanogap

3.3 B. Radiation Enhancement for an Emitter Located Near the Antenna Termination

4. CONCLUSION

5 Acknowledgment

Article Outline

关于本站 Cookie 的使用提示

全站搜索

Understanding localized surface plasmon resonance with propagative surface plasmon polaritons in optical nanogap antennas Download： 1189次

1. INTRODUCTION

2. SPP MODEL OF LSPR FOR NANOGAP ANTENNAS

Table 1. Complex Eigenfrequencies ωc (λc=2πc/ωc) of QNMs

3. REPRODUCING ENHANCEMENT FACTOR OF RADIATION WITH THE MODEL

3.2 A. Radiation Enhancement for an Emitter Located in the Nanogap

3.3 B. Radiation Enhancement for an Emitter Located Near the Antenna Termination

4. CONCLUSION

5 Acknowledgment

Article Outline

相关论文

相关资讯

关于本站 Cookie 的使用提示

全站搜索

Table 1. Complex Eigenfrequencies $ω_{c}$ ( $λ_{c} = 2 π c / ω_{c}$ ) of QNMs