In a laser-driven atomic quantum system, a continuous state couples to a discrete state resulting in quantum interference that provides a transmission peak within a broad absorption profile the so-called electromagnetically induced transparency (EIT). In the field of plasmonic metamaterials, the sub-wavelength metallic structures play a role similar to atoms in nature. The interference of their near-field coupling at plasmonic resonance leads to a plasmon induced transparency (PIT) that is analogous to the EIT of atomic systems. A sensitive control of the PIT is crucial to a range of potential applications such as slowing light and biosensor. So far, the PIT phenomena often arise from the electric resonance, such as an electric dipole state coupled to an electric quadrupole state. Here we report the first three-dimensional photonic metamaterial consisting of an array of erected U-shape plasmonic gold nanostructures that exhibits PIT phenomenon with magnetic dipolar interaction between magnetic metamolecules. We further demonstrate using a numerical simulation that the coupling between the different excited pathways at an intermediate resonant wavelength allows for a π phase shift resulting in a destructive interference. A classical RLC circuit was also proposed to explain the coupling effects between the bright and dark modes of EIT-like electromagnetic spectra. This work paves a promising approach to achieve magnetic plasmon devices.
Plasmonic metamaterials are composed of artificial metamolecules exhibiting unusual optical properties such as negative refraction index [1, 2], and toroidal dipolar response [3, 4] that can lead to applications that are otherwise unattainable in nature, such as sub-diffraction imaging [5, 6], and optical and spectrum manipulation [7–10]. These unique characteristics also enable us to explore tunability and nonlinearity that either do not exist in nature or are too weak for practical applications [11, 12]. The magnetic and electric responses of artificial metamolecules play a similar role to those of atoms in electromagnetic media. Split ring resonators (SRRs), commonly constructed building blocks of plasmonic metamaterials, have been proposed to produce electric as well as magnetic dipolar response [13–16]. The dipolar response of such SRR structures can be excited by an incident wave with either electric field perpendicular to two prongs or the magnetic field passing through the gap of SRR . While the electric dipolar response can be readily explored, the magnetic counterpart has proven to be much harder to excite because of the orientation in induced magnetic dipoles of planar SRRs, which has largely prevented a number of important applications that rely upon magnetic interactions [17, 18]. This limitation can be overcome with SRR structures arranged in a particular orientation that enables magnetic dipolar interaction between neighboring SRRs for a specific incident wave [14, 16, 19–21].
Electromagnetically induced transparency (EIT) is a quantum interference effect in a material between two pathways induced by two different incident lasers [22, 23]. By introducing a pumping laser to the material, the absorption intensity of the probe laser can be dramatically reduced, producing a controllable transparency window in electromagnetic spectra. Many electromagnetic theorems and classical models such as mass-spring systems and RLC circuits were used to analyze the mechanism of EIT [24–27]. Such a phenomenon has also been observed in plasmonic metamaterials where plasmon induced transparency (PIT) works very similarly to that of EIT, except that in the case of PIT, the mechanism is considered as the plasmon coupling of the different pathways of metamolecules between a bright element (radiative mode, which can couple to incident light) and a dark element (subradiative mode, which can hardly couple to incident light) which leads to sharp resonance and suppresses the radiative loss . Therefore, PIT provides a promising pathway towards realizing of ultra-sensitive biosensors [29, 30], slow light propagation [31, 32], and nanoruler of sub-microstructures or bio-molecular [33, 34].
Up to now, due to fabrication limitation [35, 36], most plasmonic metamaterials exhibiting PIT have been demonstrated in planar or multi-layered structures by electric coupling rather than by magnetic coupling [28, 30, 37–40]. Recently, a number of plasmonic structures with magnetic resonances have also been proposed because of their potential advantages [41, 42]. For example, Liu et al. demonstrated the propagation of electromagnetic energy via magnetic plasmons is capable of supporting low-loss propagation of plasmon modes [17, 43]. The utilization of magnetic coupling will provide an additional degree of design freedom for plasmonic devices based on magnetic plasmon material.
In this paper, we numerically and experimentally demonstrate the first PIT phenomenon in three-dimensional U-shaped SRR metamolecules in the optical region. The PIT is created by the magnetic resonance where the dark mode is excited through coupling with the bright mode via magnetic dipolar response using the broken geometrical symmetry. The transmittance intensity at the PIT peak can be controlled precisely with the manipulation of structural symmetry. We find that the coupling strength between the bright and dark modes has a direct relationship with the degree and broadness of the transparency window in transmittance spectra. Using numerical simulation, we demonstrate the π phase shift between the different excitation paths, providing a clear evidence of the destructive interference at the PIT peak. Finally, we develop a coupled circuit to explain the coupling strength between the different elements in a unit cell.
Materials and methods
Using the electron beam double exposure process with alignment technology , we fabricated three erected U-shaped gold resonance rings in a unit cell covering a total area of 75×75 μm2 on fused silica. For precision alignment of the e-beam double exposure process, two golden cross alignment marks (size 9×150 μm) with 100 nm thickness were first fabricated. A 200 nm-thick PMMA layer was spin-coated at 4000 rpm on fused silica wafer and then baked on a hot plate for 3 min at 180°C. The Espacer (Kokusai Eisei Co., Showa Denso Group, Japan) is spin-coated at 1500 rpm over the PMMA layer. Espacer is an organic polymer which can eliminate the static charging problem during the e-beam exposure process. An ELS-7000 electron beam lithography system (Elionix Inc., Tokyo, Japan) at the acceleration voltage of 100 keV with 30 picoamp (pA) of current is used for exposure. The bottom rods of erected SRR structures were defined in positive resist (495k PMMA) after the first e-beam exposure and lift-off process. Subsequently, the two prongs of erected SRR structure are made by second e-beam exposure and lift-off process.
All of the transmittance spectra were measured using the VERTEX 70 Fourier-transform infared spectrometer (Bruker Co., Ettlingen, Germany) with a Bruker HYPERION 1000 infrared microscope (15× Cassegrain objective, numerical aperture NA=0.4, near-infrared polarizer, and an InGaAs detector). An iris was used to collect the incident light to a square area of about 75×75 μm2.
Results and discussions
The geometries of a three-dimensional plasmonic metamolecule are illustrated in Figure 1A along with the polarization state of an incident light. The metamolecule (unit cell) consists of three erected U-shaped gold split-ring resonators, including one in the middle along the vertical direction (VSRR) and a pair on the side along the horizontal direction (HSRR). The gap G between the VSRR and HSRR is fixed at 30 nm and the lateral displacement between the central axes of VSRR and HSRR pair is defined as D (see Figure 1B). The feature sizes of each single SRR structure is presented in Figure 1C. The VSRR, supports a bright mode whose magnetic resonance can be excited by a normal incident light with electric field polarization along the y direction. Owing to the strong coupling between the incident light and VSRR, the bright mode suffers a significant far field scattering loss. On the contrary, the HSRR pair supporting only dark modes can be viewed as non-radiative elements because of its much weaker coupling with the incident light. The coupling effect between the incident light, VSRR and HSRR can be studied by introducing asymmetry between the SRRs with the tuning of the displacement parameter D.
In order to investigate the coupling effects between VSRR and HSRRs, we fabricated three different kinds of metamolecules corresponding to displacement parameters D=0 nm, 30 nm, and 60 nm. Figure 1D presents oblique views in the case of D=30 nm, showing the accurate alignment for the fabrication process (the misalignment is controlled to within 10 nm). To experimentally investigate the optical properties of these PIT metamolecules, near-infrared transmittance spectra are measured by Fourier-transform infrared spectrometer in the case of normal illumination with the magnetic field passing through the gap of middle SRR. The measured transmittance spectra are presented by colored curves in Figure 2. In the case of D=0 nm (dark yellow line in Figure 2), a typical transmittance profile with a maximum power absorption can be observed. The origin of this plasmonic resonance is contributed by the magnetic resonance of VSRR (we will present a more detailed analysis on this in the next section). By introducing structural asymmetry with D=30 nm (blue line in Figure 2), a transmittance peak appears within the resonance profile of D=0 nm case as a result of the near-field coupling between VSRR and HSRRs. Furthermore, the transmittance peak becomes more pronounced for D=60 nm (orange line in Figure 2).
Now, we present numerical simulation results for the various coupling effects in the plasmonic metamolecule unit cells. All of the simulation spectra and field distribution are carried out by solving the three-dimensional Maxwell equation with the commercial COMSOL Multiphysics software based on the finite element method (FEM). A single metamolecule pattern was simulated (with periodic boundary conditions) under normal-incidence illumination by linearly-polarized light in the wavelength range from 950 nm to 1750 nm. The refractive index of fused silica wafer is fixed at 1.4584. The permittivity of bulk gold in the near infrared regime is described by the Drude-Lorentz model with plasmon frequency ωp=8.997 eV and damping constant Γp=0.14 eV. The simulation results are calculated using a damping constant that is two times larger than the bulk value because of the surface defect scattering and grain effects. Figure 3 shows the simulation transmittance (Figure 3A), absorbance (Figure 3B) as a function of incident wavelength, and the magnetic field intensity (Figure 3C) for the three cases: D=0 nm, 30 nm and 60 nm. The absorbance spectra were performed using A=1-T-R, where R is reflectance. All the simulation results are in qualitative agreement with our experiment, i.e. the structural asymmetry allows for the excitation of the dark modes which are otherwise not allowed to be excited, resulting in PIT spectra. The small difference between our experimental and simulation results is likely to arise from the deviation of fabricated patterns from that of the simulation design as well as the fact that the sample surface is not as smooth as the simulation condition. Therefore, there is a slight difference compared with the simulation (Figure 3A) and measured (Figure 2) transmittance spectra.
Similarly, in the case of D=0 nm, due to the plasmonic resonance of the VSRR at the middle of metamolecule unit cell, a single transmittance dip can be observed. The plasmonic resonance of VSRR shows a strong magnetic dipolar response resulting from a ring-like current oscillation on the surface of VSRR. The schematic of the magnetic field lines for the induced magnetic dipole of VSRR is shown in the inset of Figure 3A. In this case, the HSRRs cannot be excited either by the y-polarized incident light or indirectly by the induced magnetic dipole of VSRR, since the VSRR induced magnetic field is parallel to the gap of HSRRs. By introducing structural asymmetry with D=30 nm, a transmittance peak appears (see the second row of Figure 3) within the resonance profile for D=0 nm. It reveals that increasing structural asymmetry will enhance the coupling effect between the VSRR and HSRR pair. The structural asymmetry, allows for a component of the magnetic field induced by the plasmonic resonance of VSRR to pass through the gap of HSRRs. Subsequently, according to Lenz’s law, an additional magnetic field will be induced in the HSRR pair to oppose the change of magnetic flux going through their gaps, resulting in the magnetic resonance of HSRRs. It is worth mentioning that the VSRR induced current oscillations in the two HSRRs are opposite (magnetic quadrupole) of each other (see the purple arrows of the direction of magnetic field in Figure 3C). As a result, the HSRRs supporting only dark modes exhibit a higher quality factor than that of VSRR (bright mode) because the scattering power of the induced magnetic dipoles in the HSRRs will cancel each other out. The PIT phenomenon is very difficult to observe without the high quality factor of the dark mode even when the coupling between bright and dark modes is strong. The high quality factor of the HSRRs contributes to the fact that we were able to observe and simulate a sharp transmittance peak for the case of D = 60 nm.
There are two different coupling pathways between incident light and our metamolecules. In order to investigate the phase relation between two excited pathways, the oscillating magnitude of induced magnetic field for both VSRR and HSRRs have been calculated and carried by FEM simulation. Figure 4A shows the schematic diagram (top view) with the direction of magnetic field Hx and Hy for D=60 nm. The magnitudes of induced magnetic field in x-component Hx (green arrow in Figure 4A) and y-component Hy (red arrow in Figure 4A) are extracted from the middle point of VSRR and HSRRs, respectively. Figures 4B–D show the phase revolution corresponding to Hx and Hy at wavelength λ=1285 nm, 1195 nm and 1395 nm, respectively. We find a 90° phase shift between the induced magnetic field of VSRR with respect to that of HSRR at wavelength λ=1285 nm (Figure 4B). We define the oscillating phase at the first excitation pathway (∣Incident light〉→∣VSRR〉) to be zero due to the excitation of VSRR under normal illumination. Because of the near-field coupling between VSRR and HSRRs, the second excitation pathway (∣Incident light〉→∣VSRR〉→∣HSRRs〉→∣VSRR〉) experiences a 180° phase delay as the excited HSRRs couples energy back to the VSRR resulting from π/2 phase shift between VSRR and HSRRs. As a result, a 180° phase shift between the first and second pathway leads to the destructive interference that results in a transmission peak at wavelength λ=1285 nm, similar to the EIT phenomenon in atomic physics . On the other hand, the phase shift between VSRR and HSRRs is 140° and 30° at wavelength λ=1195 nm and 1395 nm, respectively, as shown in Figures 4C and 4D. It reveals that the transmittance spectra exhibit visible dips rather than peaks without a 180° phase shift between two excited pathways.
For investigating the contribution of magnetic and magneto- electric coupling effect between VSRR and HSRRs, the enhancement of scattered power is calculated, as shown in Figure 5. The magnitude of the scattered power of electric and magnetic multipoles can be defined as [3, 4]:
electric quadrupole moment:
magnetic quadrupole moment:
where c is the speed of light, j is the induced volume current density and α, β = x, y, z. Subsequently, the scattered power of the multipole moments is calculated from the induced currents [3, 4]:
The scattered power enhancement is defined as
To provide a quantitative analysis of the PIT phenomena, we describe the coupling mechanism using a model of three coupled RLC circuits (see Figure 6A). The VSRR is described by the first circuit driven by Vin(t) = Ve-iωt which describes the in-coupling effect of the incident light. The HSRRs, on the other hand, can be excited via the induced magnetic field of VSRR, thus the coupling mechanism between the first circuit and other two circuits is accounted for as mutual inductance M21 and M31. With the identical geometric dimensions of the three SRRs inside a metamolecule unit cell, their inductance and capacitance are set as equal L1=L2=L3=L and C1=C2=C3=C. The Lagrangian of the coupled circuit system can be written as
where, Ii(t)=dQi(t)/dt, where i=1, 2, 3 is the current flowing in the circuit i. Substituting Eq. (6) to the Euler-Lagrangian equations
Here, γi=Ri/L (i=1, 2, 3) is the energy dissipation of each circuit,
Subsequently, the energy dissipation as a function of frequency can be obtained by solving the Eqs. (7), (8) and (9):
The results calculated using Eq. (10) are shown in Figure 3B as red-dotted curves with δ,γ1 and γ2 taken as fitting parameters to reproduce the simulation results where the detuning parameter δ is approximately 1.6 THz and γ1 and γ2 are roughly 35.8 THz and 20.5 THz, respectively. These values provide evidence to support our expectation of loss reduction in the HSRR pair because of the anti-parallel magnetic dipoles induced in the pair.
For enhancing the near-field coupling between the VSRR and HSRRs, the structural asymmetry must be introduced. However, not only the near-field coupling but the lateral distance between VSRR and HSRRs will be increased by introducing the structural asymmetry. Due to the near-field depending strength between plasmonic metamaterials strongly depend on the distance between each other, there will be a threshold of distance parameter D to generate the largest transmittance intensity at PIT peak. It will provide more practical applications based on the PIT characterizations so the evolution of PIT spectra can be better understood. Therefore, in order to gain a deeper realization of the evolution of the PIT spectra, the PIT transmission intensity peak as a function of the distance parameter D from 20 nm to 140 nm is shown in Figure 6B (black curve with spherical dots) for investigation. As expected, the intensity of transmittance peak increases initially with D, because the asymmetry allowed for coupling between VSRR and HSRRs. Further increase of D beyond a critical point near 80 nm actually leads to a decrease of the PIT transmittance because of the reduced coupling strength as a result of the separation increase, which is reflected by the coupling parameter Ω also shown in Figure 6B (blue curve with spherical dots).
A plasmon induced transparency in three-dimensional photonic metamolecules based on magnetic interaction has been demonstrated and investigated in the optical region. The excitation of such a magnetic quadrupole is a result of the fundamental electromagnetic theorem (Lenz’s law) which indicates that it will induce a magnetic field to oppose the change in magnetic flux. Because of the particular orientation of erected U-shaped SRRs, the strong magnetic dipolar interaction between them provides multiple pathways for near-field coupling in plasmonic metamolecules by introducing structural asymmetry. Using numerical simulation, we have proved that destructive interference takes place between different excited pathways of bright and dark elements which leads to the observed PIT phenomenon. Furthermore, to generate the maximum transmittance intensity in PIT spectra, the lateral distance between bright and dark modes is also a key factor that needs to be considered carefully. Without the proper distance between bright and dark modes, the PIT phenomenon is very difficult to observe even when the introducing structural asymmetry is very strong. This study of optical properties associated with magnetic resonance will provide important insight towards practical applications using magnetic plasmon material [30, 31, 34, 44, 45].
The authors thank financial aids from National Taiwan University, National Science Council under grant numbers 100-2923-M-002-007-MY3, 101-3113-P-002-021 and 101-2112-M-002-023. Authors are grateful to EPSRC, UK and the Royal Society, London, the National Center for Theoretical Sciences, Taipei Office, the Molecular Imaging Center of National Taiwan University and National Center for High-Performance Computing, Taiwan for their support.
 Huang Y-W, Chen WT, Wu PC, Fedotov V, Savinov V, Ho YZ, Chau Y-F, Zheludev NI, Tsai DP. Design of plasmonic toroidal metamaterials at optical frequencies. Opt Express 2012;20:1760–68.10.1364/OE.20.001760Search in Google Scholar PubMed
 Kadic M, Guenneau S, Enoch S, Huidobro PA, Martín-Moreno L, García-Vidal FJ, Renger J, Quidant R. Transformation plasmonics. Nanophotonics 2012;1:51–64.10.1515/nanoph-2012-0011Search in Google Scholar
 Chen WT, Wu PC, Chen CJ, Weng CJ, Lee H-C, Yen T-J, Kuan C-H, Mansuripur M, Tsai DP. Manipulation of multidimensional plasmonic spectra for information storage. Appl Phys Lett 2011;98:171106–3.10.1063/1.3584020Search in Google Scholar
 Ly-Gagnon DS, Balram KC, White JS, Wahl P, Brongersma ML, Miller DAB. Routing and photodetection in subwavelength plasmonic slot waveguides. Nanophotonics 2012;1:9–16.10.1515/nanoph-2012-0002Search in Google Scholar
 Castelletto S, Li X, Gu M. Frontiers in diffraction unlimited optical methods for spin manipulation, magnetic field sensing and imaging using diamond nitrogen vacancy defects. Nanophotonics 2012;1:139–153.10.1515/nanoph-2012-0001Search in Google Scholar
 Yen TJ, Padilla WJ, Fang N, Vier DC, Smith DR, Pendry JB, Basov DN, Zhang X. Terahertz magnetic response from artificial materials. Science 2004;303:1494–96.10.1126/science.1094025Search in Google Scholar PubMed
 Linden S, Enkrich C, Wegener M, Zhou JF, Koschny T, Soukoulis CM. Magnetic response of metamaterials at 100 terahertz. Science 2004;306:1351–53.10.1126/science.1105371Search in Google Scholar PubMed
 Liu N, Mukherjee S, Bao K, Brown LV, Dorfmüller J, Nordlander P, Halas NJ. Magnetic plasmon formation and propagation in artificial aromatic molecules. Nano Lett 2011;12:364–69.10.1021/nl203641zSearch in Google Scholar PubMed
 Chen WT, Chen C-J, Wu PC, Sun S, Zhou L, Guo G-Y, Hsiao CT, Yang K-Y, Zheludev NI, Tsai DP. Optical magnetic response in three-dimensional metamaterial of upright plasmonic meta-molecules. Opt Express 2011;19:12837–42.10.1364/OE.19.012837Search in Google Scholar PubMed
 Chen CC, Hsiao CT, Sun S, Yang K-Y, Wu PC, Chen WT, Tang YH, Chau Y-F, Plum E, Guo G-Y, Zheludev NI, Tsai DP. Fabrication of three dimensional split ring resonators by stress-driven assembly method. Opt Express 2012;20:9415–20.10.1364/OE.20.009415Search in Google Scholar PubMed
 Field JE, Hahn KH, Harris SE. Observation of electromagnetically induced transparency in collisionally broadened lead vapor. Phys Rev Lett 1991;67:3062–65.10.1103/PhysRevLett.67.3062Search in Google Scholar PubMed
 Davis TJ, Gó’mez DE, Vernon KC. Simple model for the hybridization of surface plasmon resonances in metallic nanoparticles. Nano Lett 2010;10:2618–25.10.1021/nl101335zSearch in Google Scholar PubMed
 Dong Z-G, Liu H, Cao J-X, Li T, Wang S-M, Zhu S-N, Zhang X. Enhanced sensing performance by the plasmonic analog of electromagnetically induced transparency in active metamaterials. Appl Phys Lett 2010;97:114101.10.1063/1.3488020Search in Google Scholar
 Liu N, Weiss T, Mesch M, Langguth L, Eigenthaler U, Hirscher M, Sönnichsen C, Giessen H. Planar metamaterial analogue of electromagnetically induced transparency for plasmonic sensing. Nano Lett 2010;10:1103–07.10.1021/nl902621dSearch in Google Scholar PubMed
 Yannopapas V, Paspalakis E, Vitanov NV. Electromagnetically induced transparency and slow light in an array of metallic nanoparticles. Physical Review B 2009;80:035104.10.1103/PhysRevB.80.035104Search in Google Scholar
 Soukoulis CM, Wegener M. Past achievements and future challenges in the development of three-dimensional photonic metamaterials. Nat Photon 2011;5:523–30.10.1038/nphoton.2011.154Search in Google Scholar
 Liu N, Langguth L, Weiss T, Kästel J, Fleischhauer M, Pfau T, Giessen H. Plasmonic analogue of electromagnetically induced transparency at the drude damping limit. Nat Mater 2009;8:758–62.10.1038/nmat2495Search in Google Scholar PubMed
 Papasimakis N, Fedotov VA, Zheludev NI, Prosvirnin SL. Metamaterial analog of electromagnetically induced transparency. Phys Rev Lett 2008;101:253903.10.1103/PhysRevLett.101.253903Search in Google Scholar PubMed
 Lassiter JB, Sobhani H, Knight MW, Mielczarek WS, Nordlander P, Halas NJ. Designing and deconstructing the fano lineshape in plasmonic nanoclusters. Nano Lett 2011;12:1058–62.10.1021/nl204303dSearch in Google Scholar PubMed
 Maier SA, Brongersma ML, Kik PG, Meltzer S, Requicha AAG, Atwater HA. Plasmonics – a route to nanoscale optical devices. Adv Mater 2001;13:1501–05.10.1002/1521-4095(200110)13:19<1501::AID-ADMA1501>3.0.CO;2-ZSearch in Google Scholar
 Liu N, Mukherjee S, Bao K, Li Y, Brown LV, Nordlander P, Halas NJ. Manipulating magnetic plasmon propagation in metallic nanocluster networks. ACS Nano 2012;6:5482–88.10.1021/nn301393xSearch in Google Scholar
©2012 by Science Wise Publishing & De Gruyter Berlin Boston