We present a study of the performance enhancement of a quantum dot infrared photodetector (QDIP), by means of complementary split-ring resonator (CSRR) nano-antennae. The QDIP is based on an asymmetric heterostructure containing a single layer of self-assembled InAs/GaAs quantum dots (QDs). The proximity of the QD plane to the top contact layer is exploited for the coupling with the near-field of the CSRR modes. The co-existence of the CSRR LC mode, at λLC = 7.4 μm, and of non-localized Bragg-like modes, is observed for the two-dimensional array of nano-antennae implemented on the QDIP. At λLC and a temperature T = 10 K, the antenna coupled device is characterized by a responsivity of 44 μA/W and a specific detectivity D* = 1.5 × 108Jones. For the highly localized LC mode, enhancements of a factor 1.7 in responsivity and 2.1 in specific detectivity are observed. Within the sub-wavelength LC mode effective surface, normalizing the overall response to the active surface of the detector, a responsivity enhancement of ∼19 is estimated, showing the potentiality of this approach for the realization of high-performance QDIPs working at normal incidence.

The inter-sublevel transitions in semiconductor quantum dots (QDs)1 can be exploited to realize photodetectors for the mid-infrared (mid-IR) spectral range.2–4 Several advantages have been envisaged for the quantum dot infrared photodetector (QDIP) with respect to the quantum well infrared photodetector (QWIP),5 such as:6 possibility of normal-incidence excitation; reduced thermal generation and relaxation rates of electrons; low dark current. Several comprehensive reviews4,7–10 have been written in order to discuss the issues and potentialities of detectors based on QDs. Among several limitations, the low absorption strength, for light polarized along in-plane directions, associated to transitions from the ground QD shell to the outermost shells or to the continuum, strongly limits the QDIP responsivity at normal incidence. Many efforts have been focused on the implementation on QDIPs of optically resonant structures, such as photonic11,12 or plasmonic crystals,13–17 aimed at increasing the coupling efficiency with normally incident light. Several systems are mentioned in Ref. 7. Another approach typically used for the enhancement of the detector performance is based on nanoantennae capable of concentrating the field in low-sized volumes.18 In Ref. 19 performance improvements of a QWIP, based on an array of patch antenna micro-cavities, were presented. In Ref. 20, the resonant coupling between a complementary split-ring-resonator (CSRR) metamaterial and the active region of a quantum cascade detector (QCD) was shown. A similar CSRR antenna concept, as we show in this letter, can be used for enhancing the responsivity of a QDIP.

For our study, as a proof of principle, we used a simple single-layer QDIP asymmetric structure analogous to the one presented in Refs. 21–23. The latter is based on a single layer of self-assembled InAs/GaAs QDs (2.4ML grown at 0.088ML/s), coupled via a GaAs(5 nm)/Al0.33Ga0.67As(3.5 nm)/GaAs(2 nm) structure to an In0.21Ga0.79As QW. By means of scanning electron microscopy (SEM), a QD density nQD = 2.1 × 1010 cm−2 was estimated.

The electronic structure of the lens shaped QDs was computed by means of Nextnano++,24 using the experimentally determined QD geometry and the inter-band spectrum obtained from the low-temperature photoluminescence (not shown).23 The conduction band (CB) profile of the structure is shown in Fig. 1, and was computed using the method described in Refs. 22 and 23. In the inset of Fig. 1, it is possible to observe how the structure was engineered in order to position the Fermi level (red dashed line) above the s shell of the QDs. The electrons sitting on the QD s shell can be excited by means of mid-IR optical radiation towards the continuum above the GaAs band edge, generating a photocurrent. This type of detector works for zero (photovoltaic operation)25 and negative values of the applied voltage Vbias.

FIG. 1.

The CB profile of the detector heterostructure, obtained by solving self-consistently the 1D Schrödinger and Poisson equations within the structure. The inset shows the active region of the QDIP consisting of a QD layer, coupled via a GaAs/AlGaAs/GaAs multi-barrier structure to the InGaAs QW. The Fermi level, at 0 V applied bias, lies above the s shell of the QDs.

FIG. 1.

The CB profile of the detector heterostructure, obtained by solving self-consistently the 1D Schrödinger and Poisson equations within the structure. The inset shows the active region of the QDIP consisting of a QD layer, coupled via a GaAs/AlGaAs/GaAs multi-barrier structure to the InGaAs QW. The Fermi level, at 0 V applied bias, lies above the s shell of the QDs.

Close modal

The layer was processed in square mesas of side lengths between 20 μm and 750 μm. The 300 nm thick n+GaAs top contact was thinned, by means of wet chemical etching, down to 50 nm. A layer of Ge/Au/Ni/Au was deposited on the bottom contact by electron beam deposition and alloyed by rapid thermal annealing. On the top contact, a layer of Ti(10 nm)/Au(100 nm) was deposited, on which the nanoantenna array was defined, using a fabrication technique based on an electron beam lithography step followed by Ar+ ion-beam etching for the metal patterning.

The nano-antenna design is based on the split-ring resonator (SRR) geometry.26 In particular, the electric SRR,27,28 consisting of two inductive loops (L) interrupted by a central capacitive gap (C) (Fig. 2(a)), is used in its complementary (CSRR) design.29,30 Such a system can be effectively treated as an LC circuit, showing a resonant behaviour at λ L C L C 31–33 (LC mode). The unit cell of the CSRR 2D square array and a SEM micrograph of a resonator obtained after the metal etching procedure are shown in Fig. 2(a).

FIG. 2.

In (a), on the left side, the unit cell of the 2D CSRR array, with side length d = 1.5 μm, and, on the right side, the SEM micrograph of a resonator obtained from the metal etching process. In (b) a schematic representation of the reference (left) and the antenna-coupled (right) devices, excited at a normal incidence. In (c) the dark current of the reference (black dashed line) and the CSRR coupled (red solid line) devices, measured at T = 10 K, is shown as a function of the applied bias voltage. In (d) the spectral response of the reference device (black dotted and dot-dashed lines) and the CSRR coupled device (blue solid and red dashed lines) are shown. The measurements were performed at T = 10 K and Vbias = 0 V. The responsivity spectra are obtained exciting the devices with light polarized along the directions ⊥ (blue solid and black dash-dotted lines) and (red dashed and black dotted lines) to the central gap of the CSRR.

FIG. 2.

In (a), on the left side, the unit cell of the 2D CSRR array, with side length d = 1.5 μm, and, on the right side, the SEM micrograph of a resonator obtained from the metal etching process. In (b) a schematic representation of the reference (left) and the antenna-coupled (right) devices, excited at a normal incidence. In (c) the dark current of the reference (black dashed line) and the CSRR coupled (red solid line) devices, measured at T = 10 K, is shown as a function of the applied bias voltage. In (d) the spectral response of the reference device (black dotted and dot-dashed lines) and the CSRR coupled device (blue solid and red dashed lines) are shown. The measurements were performed at T = 10 K and Vbias = 0 V. The responsivity spectra are obtained exciting the devices with light polarized along the directions ⊥ (blue solid and black dash-dotted lines) and (red dashed and black dotted lines) to the central gap of the CSRR.

Close modal

Besides the antenna-based devices, reference devices were processed, for which the top contact was left uncovered. Metallic pads of Ti(10 nm)/Pt(50 nm)/Au(200 nm) were finally deposited on the top contacts: they define the ohmic contacts and facilitate the wire bonding procedure. A schematic representation of the devices is shown in Fig. 2(b). For the processed devices, the distance separating the QD plane from the CSRR/GaAs interface is zQD ≃ 90 nm. The proximity between the QD plane and the CSRRs is exploited for the near-field coupling between the antenna modes and the QDIP active region. The performance of an antenna coupled device is compared with the one of a reference device without resonators (mesa surface Smesa = 750 μm × 750 μm). As it can be observed in Fig. 2(c), the dark current has the same behavior, as a function of the applied bias, on both devices. Due to the presence of the AlGaAs stopping layer, between the QW and the QDs, and the thick GaAs spacer, between the QD layer and the bottom contact, low values of dark current were observed on both devices: at Vbias= −30 mV, current densities of 2.12 × 10−10A/cm2 and 7.1 × 10−11A/cm2 were measured on the reference and the antenna coupled devices, respectively. These values are lower than the ones reported for tunnel QDIPs34,35 and QD QCDs.36 

While the spectral response of the detectors was characterized by means of FTIR spectroscopy, using a normal incidence configuration for the device excitation, as shown in Fig. 2(b), the absolute responsivity was measured exciting the device by means of a quantum cascade laser (QCL) emitting at λ = 6.8 μm, light polarized along the direction ⊥ to the central capacitive gap of the CSRR. Using a procedure similar to the one described in Ref. 25, knowing the value of the responsivity at λ = 6.8 μm, the photocurrent spectra could be calibrated.

The results are shown in Fig. 2(d), both for the reference and the antenna coupled devices, at T = 10 K and Vbias = 0 V. The presence of the resonators strongly affects the spectral response of the photodetector and, as shown in Fig. 2(d), a significant polarization anisotropy is observed. The QD photoresponse is peaked at λQD = 6.6 μm (see black dotted and dash-dotted curves in Fig. 2(d)). The spectral signature of the CSRRs is observed at λLC = 7.4 μm for the LC mode (polarization ⊥ to the capacitive gap). Several further resonances are observable, for both polarization directions, between 3.5 μm and 5.5 μm, which are associated to the excitation of Bragg-like modes, at wavelengths λ i j = ( ε G a A s ε A u ) / ( ε G a A s + ε A u ) × d / i 2 + j 2 (with i, j integers), related to the unit cell side length d.37–39 These modes have the same nature as the ones observed in metallic periodic structures already adopted for QDIPs.13–17 The splitting of the (±1, 0) and (0, ±1) resonances is attributed to slight unintentional deviations, associated to the experimental configuration, from a perfectly normal incidence condition, which introduce in-plane wavevector components.37 Other effects on the resonance positions are associated to the shape of the sub-wavelength holes.40,41

The values of the responsivity R and the quantum efficiency ξ = hc/() × R (where q is the electron charge), measured on the two devices at λLC = 7.4 μm and λ01 = 4.9 μm, are listed in Table I. The quantum efficiency ξ of the two detectors is expressed as

(1)
(2)

where NQD is the electron sheet density in the QD layer; the quantities Gref and GCSRR are coefficients quantifying the optical coupling between the incident radiation and the QDIP active region; σip and σz are the values of the cross-sections associated to the photogeneration process,22 for light polarized along the in-plane (ip) and growth (z) directions. For this type of QDIP, we estimated experimentally a ratio σz/σip ≃ 4.5. The CSRRs, for normally incident light (ip polarized), generate field components along both the ip and z directions, producing a more efficient coupling with the QDs. The experimentally estimated enhancement factors ηCSRR(λ) = RCSRR(λ)/Rref(λ) for incident light polarized along the directions and ⊥ to the central capacitive gap are shown in Figs. 3(a) and 3(b), respectively. At λLC = 7.4 μm, a responsivity enhancement ηLC = 1.7 is estimated. A significantly higher value η01 = 9.8 is estimated at λ01 = 4.9 μm. In order to study the spectral properties of ηCSRR, the 2D CSRR array was simulated by means of CST Microwave Studio®42 using the actual dimensions shown in Fig. 2(a). The simulated transmission spectra are shown in Figs. 3(a) and 3(b) (green dashed lines) for the two polarization directions of the incident field. A good agreement is found between the positions of the maxima in transmission and the ones for ηCSRR(λ). The field enhancement (FE) distributions computed for the (0, 1) and the LC modes are shown in Figs. 3(a) and 3(b), respectively. The CSRR optical coupling coefficient, corresponding to the field component El, is given by

(3)

where FE l = | E l | / | E i n | and Ein is the incident field. In Fig. 3(d), the FEip and FEz distributions of the LC mode are shown. The significantly higher value of ηCSRR observed at λ01 is explained by the larger effective surface associated to the (0, 1) mode, with respect to the localized LC mode (see Figs. 3(a) and 3(b)), therefore involving a higher number of QDs in the photogeneration process. For the reference sample, G i p r e f 0.7 is the coupling coefficient for normal incidence excitation, having considered the effects of the reflection at the vacuum/GaAs interface. Using the computed FE distributions in Eqs. (1) and (2), η L C c o m p = 3.9 and η 01 c o m p = 19.3 were calculated for the modes at λLC and λ01, respectively. The measured enhancement factors were found to be ∼2 times smaller as the computed ones. This discrepancy can be explained taking into account, for the reference sample, the effects of stray reflections within the substrate, especially at the bottom surface; this contribution is significantly lower for the antenna coupled device, due to the low transmission of the 2D CSRR interface. Moreover, the wet etching technique adopted for the top contact thickness reduction can easily produce inhomogeneity on the top surface of the mesas: regions with a top contact thickness larger than the expected value can be present within the device mesa. The LC mode, as it can be observed from the FE distributions shown in Fig. 3(b), is localized within the “plates” of the central capacitive gap and an effective surface SLC can be defined, within which the electric field E a b s ( x , y , z Q D ) | E x | 2 + | E y | 2 + | E z | 2 E a b s m a x / e ; the ratio between the areas of the unit cell Sunit = d2 and SLC is ≃ 11.3. A normalized responsivity enhancement can be calculated for the LC mode as: η L C n o r m = 11.3 × η L C 19 . As shown in Fig. 3(c), despite their proximity to the CSRR antennae, the QDs are coupled to the “tail” of the LC FE distribution ( FE a b s m a x = 4.5 at z = 90 nm). At z = 0 ( FE a b s m a x = 44 ), a responsivity enhancement factor ∼ 12 times higher is expected.

TABLE I.

Summary of the values of the responsivity (R) and of the quantum efficiency (ξ) at λLC = 7.4 μm (pol ⊥) and λ01 = 4.9 μm (pol ) for the reference (ref) and the antenna-coupled (CSRR) devices.

Mode RCSRR (μA/W) ξCSRR Rref (μA/W) ξref
λLC (pol ⊥)  44.1  7.4 × 10−6  25.6  4.3 × 10−6 
λ01 (pol 51.7  8.8 × 10−6  5.3  9 × 10−7 
Mode RCSRR (μA/W) ξCSRR Rref (μA/W) ξref
λLC (pol ⊥)  44.1  7.4 × 10−6  25.6  4.3 × 10−6 
λ01 (pol 51.7  8.8 × 10−6  5.3  9 × 10−7 
FIG. 3.

The estimated enhancement factors ηCSRR(λ) and the simulated 2D CSRR transmission spectra (green dashed line) are shown, for incident light polarized along the directions (in (a)) and ⊥ (in (b)) to the CSRR capacitive gap. The computed distributions of the absolute field enhancement (FEabs) at z = 90 nm, for the modes at λ01 and at λLC, are shown at the bottom of (a) and (b), respectively. In (c) the computed FEabs distribution of the LC mode, within the xz plane, is shown. In (d) the in-plane (FEip) and along-growth (FEz) field enhancement distributions are shown, at zQD = 90 nm, for the LC mode. From the distribution of FEabs it is possible to obtain the area of the LC mode surface SLCd2/11.3 (indicated in (b)).

FIG. 3.

The estimated enhancement factors ηCSRR(λ) and the simulated 2D CSRR transmission spectra (green dashed line) are shown, for incident light polarized along the directions (in (a)) and ⊥ (in (b)) to the CSRR capacitive gap. The computed distributions of the absolute field enhancement (FEabs) at z = 90 nm, for the modes at λ01 and at λLC, are shown at the bottom of (a) and (b), respectively. In (c) the computed FEabs distribution of the LC mode, within the xz plane, is shown. In (d) the in-plane (FEip) and along-growth (FEz) field enhancement distributions are shown, at zQD = 90 nm, for the LC mode. From the distribution of FEabs it is possible to obtain the area of the LC mode surface SLCd2/11.3 (indicated in (b)).

Close modal

In order to estimate the specific detectivity D*, we measured the background current IBG generated by the devices without a cryo-shield (NA = 1) and in absence of an incident optical signal: values of 1.9 nA and 2.9 nA were measured for the antenna-coupled and reference devices, respectively. From the study of the temperature dependence of Idark and IBG, TBLIP ≃ 55 K was estimated for both devices. Analogously to what was shown in Ref. 25, taking into account the thermal noise and the dark and BG contributions to the shot noise,43,44 we estimated D * = R S m e s a / 2 q ( I d a r k + I B G ) + 4 k B T / R Q D I P (where R Q D I P is the QDIP differential resistance). For λ = λLC, at Vbias= −30 mV and T = 10 K , D C S R R * = 1.5 × 10 8 J o n e s , was estimated for the CSRR coupled device and D r e f * = 0.7 × 10 8 J o n e s , for the reference device. The enhancement of the specific detectivity (∼2.1) confirms the performance improvements observed for the responsivity.

The relatively low values of R and D* are inherently related to the detector structure based on a single layer of QDs, which has a significantly lower efficiency with respect to multi-layer QDIPs, for which responsivity values more than five orders of magnitude higher are reached.45–47 A significant increase of the quantum efficiency can be achieved by including an extraction layer, enhancing the escape of the photoexcited electrons from the QDs. In Ref. 25, for a similar structure based on InGaAs capped QDs, a quantum efficiency about sixty times higher was reported, for photovoltaic operation.

In conclusion, we have demonstrated the responsivity enhancement on a single-layer QDIP, driven by the coupling of the detector active region with the localized LC mode of CSRR antennae. Within the LC active surface, a normalized responsivity enhancement of ∼19 was estimated. The great potentiality of such approach can be exploited for the realization of a next generation of devices with higher values of responsivity and detectivity. Etching techniques can be used to pattern the top mesa surface allowing both to remove inactive QDs and to position the detector active region within the capacitor plates” of the resonators, below the top surface level. This would lead to a significant increase of the responsivity and the reduction of the dark current.19 An improved SRR near-field antenna coupling approach could be suitable also for the implementation on multiple-layer QDIPs. Within the mid-IR spectral range considered in this letter, only a limited number of QD layers could be included within the LC mode volume. For far-infrared and terahertz QDIPs, due to the larger size of the resonators, a more efficient coupling could be attained for structures based on a larger number of QD layers.

This work was financed by the Swiss National Science Foundation (SNF). We acknowledge the support from the technical staff of the FIRST-lab and the Binnig and Rohrer Nanotechnology Center (BRNC).

1.
D.
Bimberg
,
M.
Grundmann
, and
N. N.
Ledentsov
,
Quantum Dot Heterostructures
(
John Wiley & Sons
,
1999
), p.
328
.
2.
J.
Phillips
,
K.
Kamath
, and
P.
Bhattacharya
,
Appl. Phys. Lett.
72
,
2020
(
1998
).
3.
D.
Pan
,
E.
Towe
, and
S.
Kennerly
,
Appl. Phys. Lett.
73
,
1937
(
1998
).
4.
H. C.
Liu
, “
Quantum dot infrared photodetector
,”
Opto-Electron. Rev.
11
,
1
(
2003
).
5.
H.
Schneider
and
H. C.
Liu
,
Quantum Well Infrared Photodetectors: Physics and Applications
(
Springer
,
2006
), p.
248
.
6.
J.
Phillips
,
J. Appl. Phys.
91
,
4590
(
2002
).
7.
J.
Wu
,
S.
Chen
,
A.
Seeds
, and
H.
Liu
,
J. Phys. D: Appl. Phys.
48
,
363001
(
2015
).
8.
9.
P.
Martyniuk
and
A.
Rogalski
,
Prog. Quantum Electron.
32
,
89
(
2008
).
10.
A.
Barve
,
S.
Lee
,
S.
Noh
, and
S.
Krishna
,
Laser Photonics Rev.
4
,
738
(
2010
).
11.
K. T.
Posani
,
V.
Tripathi
,
S.
Annamalai
,
N. R.
Weisse-Bernstein
,
S.
Krishna
,
R.
Perahia
,
O.
Crisafulli
, and
O. J.
Painter
,
Appl. Phys. Lett.
88
,
151104
(
2006
).
12.
R. V.
Shenoi
,
D. A.
Ramirez
,
Y.
Sharma
,
R. S.
Attaluri
,
J.
Rosenberg
,
O. J.
Painter
, and
S.
Krishna
, in
Optical Engineering + Applications
, edited by
E. W.
Taylor
and
D. A.
Cardimona
(
International Society for Optics and Photonics
,
2007
), p.
67130P
.
13.
C.-C.
Chang
,
Y. D.
Sharma
,
Y.-S.
Kim
,
J. A.
Bur
,
R. V.
Shenoi
,
S.
Krishna
,
D.
Huang
, and
S.-Y.
Lin
,
Nano Lett.
10
,
1704
(
2010
).
14.
C.-Y.
Chang
,
H.-Y.
Chang
,
C.-Y.
Chen
,
M.-W.
Tsai
,
Y.-T.
Chang
,
S.-C.
Lee
, and
S.-F.
Tang
,
Appl. Phys. Lett.
91
,
163107
(
2007
).
15.
S.
Hellström
,
Z.-H.
Chen
,
Y.
Fu
,
M.
Qiu
,
R.
Soltanmoradi
,
Q.
Wang
, and
J. Y.
Andersson
,
Appl. Phys. Lett.
96
,
231110
(
2010
).
16.
S. C.
Lee
,
S.
Krishna
, and
S. R. J.
Brueck
,
Opt. Express
17
,
23160
(
2009
).
17.
S. J.
Lee
,
Z.
Ku
,
A.
Barve
,
J.
Montoya
,
W.-Y.
Jang
,
S. R. J.
Brueck
,
M.
Sundaram
,
A.
Reisinger
,
S.
Krishna
, and
S. K.
Noh
,
Nat. Commun.
2
,
286
(
2011
).
18.
L.
Novotny
and
N.
van Hulst
,
Nature Photonics
5
,
83
90
(
2011
).
19.
Y. N.
Chen
,
Y.
Todorov
,
B.
Askenazi
,
A.
Vasanelli
,
G.
Biasiol
,
R.
Colombelli
, and
C.
Sirtori
,
Appl. Phys. Lett.
104
,
031113
(
2014
).
20.
A.
Benz
,
M.
Krall
,
S.
Schwarz
,
D.
Dietze
,
H.
Detz
,
A. M.
Andrews
,
W.
Schrenk
,
G.
Strasser
, and
K.
Unterrainer
,
Sci. Rep.
4
,
4269
(
2014
).
21.
L.
Nevou
,
V.
Liverini
,
P.
Friedli
,
F.
Castellano
,
A.
Bismuto
,
H.
Sigg
,
F.
Gramm
,
E.
Müller
, and
J.
Faist
,
Nat. Phys.
7
,
423
(
2011
).
22.
G.
Cerulo
,
L.
Nevou
,
V.
Liverini
,
F.
Castellano
, and
J.
Faist
,
J. Appl. Phys.
112
,
043702
(
2012
).
23.
G.
Cerulo
,
V.
Liverini
,
M.
Beck
,
F.
Gramm
,
H.
Sigg
, and
J.
Faist
“A study of the performance of optically driven electron pumps based on quantum dots with different spectral properties” (unpublished).
24.
Nextnano++ is a software for the simulation of electronic and optoelectronic semiconductor nanodevices and materials, produced by nextnano GmbH, http://www.nextnano.com.
25.
L.
Nevou
,
V.
Liverini
,
F.
Castellano
,
a.
Bismuto
, and
J.
Faist
,
Appl. Phys. Lett.
97
,
023505
(
2010
).
26.
J.
Pendry
,
A.
Holden
,
D.
Robbins
, and
W.
Stewart
,
IEEE Trans. Microwave Theory Tech.
47
,
2075
(
1999
).
27.
W. J.
Padilla
,
M. T.
Aronsson
,
C.
Highstrete
,
M.
Lee
,
A. J.
Taylor
, and
R. D.
Averitt
,
Phys. Rev. B
75
,
041102
(
2007
).
28.
D.
Schurig
,
J. J.
Mock
, and
D. R.
Smith
,
Appl. Phys. Lett.
88
,
041109
(
2006
).
29.
H.-T.
Chen
,
J. F.
O'Hara
,
A. J.
Taylor
,
R. D.
Averitt
,
C.
Highstrete
,
M.
Lee
, and
W. J.
Padilla
,
Opt. Express
15
,
1084
(
2007
).
30.
F.
Falcone
,
T.
Lopetegi
,
M. A. G.
Laso
,
J. D.
Baena
,
J.
Bonache
,
M.
Beruete
,
R.
Marqués
,
F.
Martín
, and
M.
Sorolla
,
Phys. Rev. Lett.
93
,
197401
(
2004
).
31.
J.
Baena
,
J.
Bonache
,
F.
Martin
,
R.
Sillero
,
F.
Falcone
,
T.
Lopetegi
,
M.
Laso
,
J.
Garcia-Garcia
,
I.
Gil
,
M.
Portillo
, and
M.
Sorolla
,
IEEE Trans. Microwave Theory Tech.
53
,
1451
(
2005
).
32.
J. F.
O'Hara
,
E.
Smirnova
,
A. K.
Azad
,
H. T.
Chen
, and
A. J.
Taylor
,
Act. Passive Electron. Compon.
2007
,
49691
.
33.
D.
Zhu
,
M.
Bosman
, and
J. K. W.
Yang
,
Opt. Express
22
,
9809
(
2014
).
34.
A. V.
Barve
,
S. Y.
Shah
,
J.
Shao
,
T. E.
Vandervelde
,
R. V.
Shenoi
,
W.-Y.
Jang
, and
S.
Krishna
,
Appl. Phys. Lett.
93
,
131115
(
2008
).
35.
X. H.
Su
,
J.
Yang
,
P.
Bhattacharya
,
G.
Ariyawansa
, and
A. G. U.
Perera
,
Appl. Phys. Lett.
89
,
031117
(
2006
).
36.
A. V.
Barve
and
S.
Krishna
,
Appl. Phys. Lett.
100
,
021105
(
2012
).
37.
T. W.
Ebbesen
,
H. J.
Lezec
,
H. F.
Ghaemi
,
T.
Thio
, and
P. A.
Wolff
,
Nature
391
,
667
(
1998
).
38.
W. L.
Barnes
,
A.
Dereux
, and
T. W.
Ebbesen
,
Nature
424
,
824
(
2003
).
39.
C.
Genet
and
T. W.
Ebbesen
,
Nature
445
,
39
(
2007
).
40.
K. J. K.
Koerkamp
,
S.
Enoch
,
F. B.
Segerink
,
N. F.
van Hulst
, and
L.
Kuipers
,
Phys. Rev. Lett.
92
,
183901
(
2004
).
41.
K. L.
van der Molen
,
K. J.
Klein Koerkamp
,
S.
Enoch
,
F. B.
Segerink
,
N. F.
van Hulst
, and
L.
Kuipers
,
Phys. Rev. B
72
,
045421
(
2005
).
42.
CST is a tool for the 3D simulation for electromagnetic design and analysis https://www.cst.com/.
43.
E.
Rosencher
and
B.
Vinter
,
Optoelectronics
(
Cambridge University Press
,
2002
).
44.
J.-M.
Liu
,
Photonic Devices
(
Cambridge University Press
,
Cambridge
,
2005
).
45.
S.
Chakrabarti
,
A.
Stiff-Roberts
,
P.
Bhattacharya
,
S.
Gunapala
,
S.
Bandara
,
S.
Rafol
, and
S.
Kennerly
,
IEEE Photonics Technol. Lett.
16
,
1361
(
2004
).
46.
H. S.
Ling
,
S. Y.
Wang
,
C. P.
Lee
, and
M. C.
Lo
,
Appl. Phys. Lett.
92
,
193506
(
2008
).
47.
J.
Szafraniec
,
S.
Tsao
,
W.
Zhang
,
H.
Lim
,
M.
Taguchi
,
A. A.
Quivy
,
B.
Movaghar
, and
M.
Razeghi
,
Appl. Phys. Lett.
88
,
121102
(
2006
).