# Double Fano resonances in a composite metamaterial possessing tripod plasmonic resonances

###### Abstract

By embedding four-rod resonators inside double-split ring resonators superlattice, a planar composite metamaterial possessing tripod plasmonic resonances is fabricated. Double Fano resonances are observed where a common subradiant driven oscillator is coupled with two superradiant oscillators. As a classical analogue of four-level tripod atomic system, the extinction spectrum of the composite metamaterial exhibits a coherent effect based on double Fano resonances. Transfer of the absorbed power between two orthogonal superradiant oscillators is shown to be mediated by the common subradiant oscillator.

###### pacs:

81.05.Xj, 78.67.Pt## I Introduction

Control of Fano resonance has been a focus of research in application of artificially structured nano-sized materials including quantum dots, nano-wires, photonic-crystal slabs, plasmonic nanostructures, as well as metamaterials. verellen2009fano (); 8.RevModPhys.82.2257 (); luk2010fano (); rahmani2013fano () Metamaterials can be designed to possess a broad bright mode and a non-resonant dark mode, which couples each other to exhibit a destructive interference effect such as classical analogue of electromagnetically induced transparency (EIT) phenomenon, corresponding to Fano resonance with the asymmetry parameter equal to zero 15.PhysRevLett.101.047401 (); 1.nature-material-Stutgart (); 14.nano-letters-Stutgart1 ().

As a further development to engineer metamaterials possessing a higher sensitivity as well as a broad spectral response for sensor applications, multiple Fano resonances are investigated in plasmonic metamaterials. Multiple Fano resonances can be achieved in two different schemes. The first scheme is to employ one bright mode and two or more dark modes. Examples are abundant where a bright dipole couples with two dark dipoles or quadrupoles: one bright dipole resonator and a pair of split ring resonators (SRRs) as the dark elementsliu2012electromagnetically (), double symmetrical U-shaped SRRs with a nanorod between the two SRRswang2013double (), plasmonic clusters of heptamerliu2012multiple () and pentamerliu2013tuning (), and 3D metamaterialsliu2011three (); davis2012analytical (); artar2011multispectral ().

The second scheme is to employ one dark mode and two or more bright modes, which already has been utilized in a four-level tripod atomic system to achieve quantum coherent effects such as polarization photon gaterebic2004polarization (), enhanced cross-phase modulationli2008enhanced (), and two slow-pulse matchingmacrae2008matched (). In the four-level tripod atomic system, two different low-lying energy levels and , bright states, are coupled to a common excited state by weak probe and trigger lasers, while the third low-lying energy level , dark state, is coupled to the excited state by a strong pumping laser to provide pathways to both probe and trigger lasers necessary for double EITs. While a quantum coherence of and is maintained, double EIT-based coherent effects are observed. As a classical analogue, a theoretical analysis of Fano correlation effect is reported for a plasmon-exciton-plasmon interaction in an artificial hybrid molecular system, one quantum dot placed in the gap of two identical metal nanoparticleshe2013fano (). A classical model has already been introduced, where two bright oscillators with different resonance frequencies are coupled to one dark oscillatorbai2012classical ().

In this work, we experimentally demonstrate a classical analogue of double Fano resonances in a planar composite metamaterial possessing tripod plasmonic resonances, where a common subradiant dipole corresponding to the dark mode is coupled with two superradient dipoles corresponding to two bright modes. The composite metamaterial is structured such that four-rod resonators (FRR) are embedded inside double-split ring resonators (DSRR) superlattice. See Fig. 1. Two dipole resonances of DSRR are superradiant modes corresponding to the excitations of and , and one dipole resonance of FRR is subradiant mode corresponding to the excitation of . Proximity of the inner diameter of DSRR and the rod-length of FRR permits near-field coherent couplings, and the composite metamaterial is a tripod metamaterial system possessing two superradiant dipole resonances of DSRR coupled coherently with one common subradiant dipole resonance of FRR in near-field, similar to a tripod atomic system with four atomic levels coupled coherently by coherent optical fields.

Important finding is that double Fano resonances in the composite metamaterial are correlated, showing up as a transfer of the absorbed power from one superradiant dipole to the other superradient dipole in DSRR.

In Sec.2, we introduce a planar composite metamaterial possessing tripod plasmonic resonances, and THz extinction spectrum of the fabricated sample is presented along with the finite-difference-time-domain (FDTD) simulated spectrum. In Sec.3, double Fano resonances are analyzed analytically by a set of coupled equations of motion of two superradiant oscillators and one subradiant oscillators to identify coherent processes of the energy transfer between two superradiant oscillators. In particular, the phase evolutions of superradiant and subradiant oscillators are closely examined. In Sec.4, the multi-polar nature of near-field coherent interaction is examined by electric field distribution obtained from FDTD theoretical simulation. Sec.5 summarizes the major findings with conclusion.

In Appendix the general feature of plasmonic Fano resonance is examined where both superradiant and subradiant oscillators are externally driven. Analysis based on two coupled oscillators model leads to the fact that the characteristic asymmetric Fano resonance formula of plasmonic structure is kept the same with a modification in the asymmetry parameter .

## Ii Sample design, fabrication, and its terahertz spectra

### ii.1 Sample design and fabrication

Double split ring resonator (DSRR) is composed of two split rings, and we can construct two different kinds of superlattice arrays of DSRR by arraying them in square lattice either symmetrically or asymmetrically. We compare extinction spectra of symmetric and asymmetric superlattice arrays in Fig. 2, as discussed in kang2010optical (). When the gap orientation of DSRR is along -axis corresponding to the symmetric superlattice array, two distinct resonances of a low energy and a high-energy take place in Fig. 2(a). While the high-energy resonance is excited for -polarization, the low-energy resonance is excited for -polarization. From the extinction spectra of the symmetric superlattice array we find that DSRR exhibits two orthogonal dipolar resonances, a low-energy and a high-energy , which is equivalent to the notion that one DSRR accommodates two independent dipole oscillators with resonances at and kang2010optical ().

In the asymmetric superlattice array of Fig. 2(b), DSRRs are arrayed in a square lattice such that the gap orientation of DSRR is rotated by angles alternately with respect to -axis, which allows simultaneous excitations of both low-energy and high-energy resonances of DSRR upon incidence of a linearly polarized electromagnetic wave. The amount of relative excitations of and dipole oscillators depends on the angle between the polarization direction and the gap orientation of DSRR17.JHWoo (). We pay a particular attention to the extinction spectrum for -polarization, i.e., the black curve in Fig. 2(b). We observe that the high-energy is excited as well as the low-energy in the asymmetric superlattice array. Refer to Fig. 7(a) for a discussion of the high-energy excitation. In the study of double Fano resonances associated with two bright modes and one dark mode, it is necessary to have simultaneous excitations of two bright modes. In this respect, the asymmetric superlattice array provides an appropriate meta-structure where both high-energy and low-energy are excited as bright modes for -polarization. On the other hand, FRR possesses a polarization-independent broad resonance, , which will be examined more carefully later in connection with double Fano resonances.

Metamaterials are fabricated by a standard photo-lithography and subsequent lift-off process of Au film with micrometer resolution, which is prepared from e-beam evaporation of 10nm thick titanium adhesion layer and 200nm thick gold layers on 500m thick -doped Si wafer. We call DSRR metamaterial as MM, FRR metamaterial as MM, and the composite metamaterial as MM, respectively, and optical microscopy images are shown in Fig. 3 (a), (b), and (c), with structural dimensions. MM is an asymmetric superlattice array of DSRRs with the inner radius =14 m, the outer radius =18 m, the width of circle line 4m, which are arranged with the gap orientation and the lattice constant 50 m. Previous studies showed that MM exhibits two polarization-dependent dipole resonances with different Q-valueskang2010optical (); 17.JHWoo (); lee2013reflection (). MM is an array of FRRs with , =26 m, =4 m, and the lattice constant 50 m, possessing a 4-fold rotational symmetry. The relative orientation of rods, i.e., , is set such that the electromagnetic response of MM is the same regardless of the polarization angle of a normally incident electromagnetic wave in an effective medium approximation. That is, MM is isotropic with respect to the polarization direction. In MM, the relative orientation of split-rings in DSRR and rods in FRR results in non-symmetric dipolar interactions between rod of FRR and split-ring of DSRR for the upper and lower arc split-rings.

### ii.2 Terahertz spectrum measurement and FDTD simulation spectrum

It is important to distinguish the electromagnetic characteristics of Au metal in the spectral range of THz and visible regimes. As well-known, in THz frequency regime the real part of dielectric constant of Au has a very large negative value with a small positive imaginary part, which leads that Au behaves as a perfect electric conductor and Au cannot support surface plasmons, electromagnetic surface excitations localized near the surface. However, there exist THz meta-resonances in metamaterials prepared in Au film, originating from an electromagnetic resonance in subwavelength structures, which is called plasmonic resonancependry1998low (); pendry1999magnetism (); pendry2004mimicking (); ma2011plasmon (); wang2013double (); wallauer2013fano ().

Time-domain polarized terahertz transmission measurements are carried out by a Teraview TPS Spectra 3000 Spectrometer with a resolution of 1.2 cm at room temperature in vacuumteraview (). Sample size is and THz beam diameter is about 5.

FDTD simulations of electromagnetic responses in MM, MM, and MM are carried out by commercial software FDTD solution Lumericallumerical (). Thickness of Au film in metamaterial samples is 200 . However, in the FDTD simulation we adopted the thickness of Au film as = 4 , which satisfies the relation , with the setting of an override maximum mesh step size of dz = 1 and minimum mesh step size 0.00025 . As far as the relation holds, the thickness of Au is not relevant, since Au film itself does not support any resonances and THz resonances originate from the subwavelength structure of metamateriallumerical (); chen2006active (); fdtd ().

A unit cell of metamaterial has been set as the unit cell of FDTD simulation and the periodic boundary conditions are adopted in - and -directions and the perfectly matched layer boundary condition is adopted in -direction. In order to identify whether the broad high-energy resonance of MM stems from a coupling between neighboring DSRRs, we examined FDTD simulated extinction spectrum by varying the horizontal distance between DSRRs. And we found that there always exists the broad high-energy resonance with the identical spectral shape except for a decrease in the magnitude of extinction simply owing to a reduced number of meta-particles per unit area.

For the dielectric permittivity of -doped silicon substrate we employed the values obtained from THz spectrum measurement, and the dielectric permittivity of Au is adopted from CRC handbook weast1988crc (), which gives a typical value of at 2.0 THz.

Figure 4 shows the polarization-dependent optical density plots of (a) measured and (b) simulated THz extinction spectra of MM for a normally incident wave with -polarization, , (red curve) and -polarization, , (blue curve). Also shown are the extinction spectra of MM with -polarization, , (solid black curve) and -polarization, , (dashed black curve) and the extinction spectrum of MM (gray curve). MM exhibits polarization-independent broad resonance (gray curve), and the composite metamaterial MM shows polarization-dependent resonances (red and blue curve), which is associated with polarization-dependent meta-resonances of MM. Comparison of Fig. 4(a) and (b) shows that there is an excellent agreement between measured and simulated extinction spectra.

In the study of double Fano resonances taking place in MM, we closely examine the extinction spectra with -polarization, , which are re-plotted in Fig. 5. As a digression we deconvolute the extinction spectra of MM to 2 Lorentzian resonances. Figure 6 shows how the extinction spectrum of MM is composed of Lorentzian resonances. We find that the broad extinction spectrum of MM is a linear sum of one resonance with a width corresponding to (green curve) and the other resonance with a width corresponding to (cyan curve). Comparison of the sum of 2 Lorentzian resonances (solid gray curve) and the measured extinction spectrum (dashed gray curve) shows a good agreement. As will be discussed later, the Lorentzian resonance takes part in giving rise to double Fano resonances, while the Lorentzian resonance stays as a background absorption since the spectral location of resonance is far away from bright mode resonances of MM. Refer to Fig. 7 and its discussion as well as Tab. 1.

Let’s examine the characteristics of MM and MM spectra in Fig. 5(a). MM (black) possesses one low-energy resonance with a large oscillator strength at = 1.4THz and the other high-energy resonance with a small oscillator strength stretching with spanning from 1.8THz to 3.0THz. On the other hand, MM (gray) possesses a broad resonance which is a linear sum of = 1.8THz and = 2.3THz.

Comparison of the linear absorption spectra of MM and MM shows two notable features. First, the total absorption cross section of MM is much larger than that of MM, indicating that the dipolar interaction with an incident electric field is much stronger in MM than in MM. Therefore, MM is composed of a strongly driven MM oscillator and a weakly driven MM oscillator with the two interacting by a near-field coupling. Second, the entire spectrum of MM is contained under the envelope of MM spectrum, with resonance in MM positioned between resonance and resonance in MM.

Next, we examine the difference in extinction spectrum between MM and MM qualitatively by paying attention to the low-frequency and high-frequency peaks in Fig. 5(a). At the low-frequency peak, MM (red) spectrum shows a dip at 1.6THz, the spectral position being on the blue side of MM (black) low-energy resonance = 1.4THz. On the other hand, at the high-frequency peak, the extinction is enhanced almost two-fold at 2.4THz in MM (red) spectrum when compared with MM (black) spectrum. These features at low- and high-frequency peaks show that there exists a correlation between changes in the low-frequency and high-frequency peaks. We note that the spectral positions of dip at 1.6THz and enhancement in the extinction at 2.4THz are located on the red and blue sides of MM (gray) resonance =1.8THz, respectively.

## Iii Double Fano resonances

### iii.1 Correlation of double Fano resonances

Now we examine the correlation between changes in the low-frequency and high-frequency peaks in terms of double Fano resonances analytically. Let’s call the superradiant dipole with resonance and the superradiant dipole with resonance in DSRR as and , respectively, and the subradiant dipole with resonance in FRR as . In the composite metamaterial MM, DSRRs are arrayed with the orientation angle , and upon normal incidence of -polarization wave, both and are excited.

Noting that a pair of superradiant dipoles at the upper and lower arcs of DSRR interact with subradiant dipole of FRR symmetrically, dipolar interactions between the superradiant dipoles & and the subradiant dipole are expressed as

(1) |

where and are the coupling strengths of and with , respectively. is the frequency associated with the characteristic coherent coupling between superradiant dipoles and subradiant dipole.

As shown in Appendix, an asymmetric Fano lineshape takes place even when the subradiant oscillator is externally driven. In the composite metamaterial MM, FRR resonator corresponds to the externally driven subradiant oscillator. In double Fano resonances we have two superradiant oscillators, and from Eq. (5) and (6) in Appendix the coupled equations of motion are straightforward with superradiant and superradiant and subradiant externally driven by a normally incident -polarization wave.

(2) |

(3) |

(4) |

Note that the subradiant dipole is externally driven with a driving amplitude of .

In the absence of near-field couplings with the subradiant oscillator , i.e., , two superradiant oscillators and are excited independently with the driving amplitudes of and , respectively. When near-field couplings of and with are turned on, i.e., and , there occurs an interplay between and through coupling . That is, oscillator and oscillator reside in one and the same DSRR, and Fano resonances in and are correlated. This is similar to the control of EITs in a tripod atomic system where the absorption spectrum of the hyperfine field shows a dependence on the powers of the Zeeman EIT signalmacrae2008matched (). A classical example is to control the absorbed power spectrum of one metal nanoparticle by manipulating the distance between quantum dot and the other metal nanoparticle in a system composed of one quantum dot placed in the gap of two identical metal nanoparticleshe2013fano ().

Analytical expression of double Fano resonances by the coupled equations of motion Eq. (2)-(4) allows to identify how the two independent orthogonal superradiant oscillators interact indirectly via near-field couplings with the common subradiant oscillator.

In Fig. 7 we compare the spectra of absorbed power obtained from Eq. (2)-(4) with the measured extinction spectra. The spectra of absorbed power of both isolated and coupled dipoles are calculated from the relation , . When near-field couplings are absent, i.e., , the red, blue, and green curves in Fig. 7(a) correspond to the absorbed powers in the isolated Lorentzian oscillators with resonances at and of DSRR in MM and of FRR in MM, respectively, We set the parameters as shown in Table 1 for resonance frequency in sec, damping in sec, and quality factor Q of isolated Lorentzian oscillators , , and of Fig. 7(a). Note that we adopt the Lorentzian resonance and with in Fig. 6 as the subradiant mode . The quality factor Q of oscillator of FRR, the isolated dipole, is the highest, asserting that oscillator is subradiant.

u(t), | red | v(t), | blue | , | green | |||
---|---|---|---|---|---|---|---|---|

47 | 15 | 3.1 | 70 | 25 | 2.8 | 60 | 14 | 4.3 |

f | g | |

1.0 | 1.0 | 0.1 |

In the presence of near-field couplings, i.e., and , two independent superradiant oscillators and residing in DSRR couple with a common subradiant oscillator of FRR via near-field couplings. With the parameters listed in Table 2, we obtain the absorbed powers of each coupled oscillator as shown in Fig. 7(b). Gray, pink, and dark brown curves correspond to the absorbed powers of and oscillators of DSRR and oscillator of FRR in MM, respectively, in the presence of near-field couplings. Black curve is the sum of gray, pink, and dark brown curves. Near-field couplings of oscillator (red) and oscillator (blue) with the subradiant oscillator (green) give rise to the characteristic asymmetric Fano resonances of the coupled dipoles in Fig. 7(b). Gray curve of Fig. 7(b) shows that there occurs a Fano dip in the high energy side of resonance of isolated (red) of Fig. 7(a). Pink curve of Fig. 7(b) shows that there occurs a Fano dip in the low energy side of resonance of isolated (blue) of Fig. 7(a). While the dipolar excitations in oscillator and oscillator are independent, the resulting double Fano resonances in the black curve show that the decrease and increase in the absorbed powers of DSRR are correlated by the common oscillator.

For comparison, in Fig. 7(c) and (d), we plotted the measured spectra of MM and MM in purple and green curves, respectively, along with MM spectrum in black curve. We find that the relation of MM MM spectra and MM spectrum is well described by the coupled equations of motion Eq. (2)-(4). In other words, in the composite MM possessing tripod plasmonic resonances, the correlation between Fano dip and extinction enhancement stems from an interplay of and oscillators in double Fano resonances, mediated by coupling coherently coupled to the common narrow oscillator.

### iii.2 Phases in double Fano resonances

In Fano resonance, a destructive interference between superradiant and subradiant modes can be understood by examining how phases of superradiant and subradiant modes evolve upon increasing the driving frequencyjoe2006classical (). In the case that the resonance frequency of subradiant mode is higher than that of superradiant mode, the superradiant mode experiences a phase of when going through its own superradiant resonance and the phase keeps on increasing, and when passing through the subradiant resonance the phase experiences a rapid decrease to the minimum phase of a value less than and recovers back to its original value, eventually accumulating to a total value of . The phase evolution of subradiant mode, however, is very different. The subradiant mode does not experience a decrease in the phase but keeps on accumulating the phase to a total value of by going through the superradiant resonance and its own subradiant resonance. Hence near the subradiant resonance there occurs a spectral region where superradiant and subradiant modes oscillate out-of-phase, which gives rise to a destructive interference for Fano resonance.

When double Fano resonances take place in tripod plasmonic resonances, on the other hand, two superradiant modes are coupled to the common subradiant mode, and the phase evolutions of two superradiant modes are not independent of each other but are correlated in a way determined by the relative spectral positions of the superradiant resonances with respect to the subradiant resonance.

In Fig. 8 we closely examine phase and amplitude of dipole oscillators comparing with the absorbed powers , where =, =, and =, where is an external force driving dipole oscillator, . Figure 8(a-1), (a-2), and (a-3) show phase , and Fig. 8(b-1), (b-2), and (b-3) show amplitude , and Fig. 8(c-1), (c-2), and (c-3) show absorbed power of dipole oscillator, . In all the figures, as a reference, we show phase , amplitude , and absorbed power of isolated dipoles in the absence of near field couplings, , in red, blue, and green curves. Red and blue curves correspond to the superradiant isolated dipoles and , respectively, while green curve corresponds to the subradiant isolated dipole with the parameters of Lorentzian oscillators shown in Table 1.

We consider 3 different cases of near-field couplings, i.e., (left panel), (center panel), and (right panel). In the presence of coupling between superradiant modes and subradiant mode red, blue, and green curves in the plots correspond to the isolated dipoles From a comparison of red and gray curves in Fig. 8 (a-1), we find that the phase evolution of (gray) exhibits a dip near the resonance of subradiant (dark brown), a feature characteristic of a single Fano resonance. Comparison of green and dark brown curves shows that the phase of subradiant (dark brown) increases to when passing the resonance of superradiant (red) which takes place before its own resonance of (green). In Fig. 8 (a-2) we observe a similar behavior in the phase evolution of (pink), but the phase of subradiant (dark brown) increases to by passing its own resonance of (green) before it comes across the resonance of superradiant (blue). We note that the spectral positions of phase of the subradiant dipole are located below and above its resonance when coupled to and , respectively. It is evident that the phase dips in superradiant and superradiant of Fig. 8 (a-1) and (a-2) are responsible for single Fano resonance as seen in the extinction spectra of Fig. 8 (c-1) (gray) and (c-2) (pink). Black curves are the extinction spectra of MM when a single Fano resonance takes place. That is, black curve is the sum of gray, blue, and dark brown curves in Fig. 8 (c-1), and is the sum of red, pink, and dark brown curves in Fig. 8 (c-2).

A numerical integration of absorbed powers shows that the sum of areas under red, blue, and green curves is identical to the sum of areas under gray, pink, and dark brown curves. Negative values of the absorbed powers by the coupled dipoles (pink) and (dark brown) which has the additional phase shifts stand for a transfer of powers to the other dipole, which is observed in Fano resonance in plasmonic structuresPhysRevLett.108.077404 (). However, we note that the amplitudes of the coupled dipoles and are always positive-definite.

In double Fano resonances, however, the phases of and are not independent of each other since a common subradiant dipole is coupled to both and . In particular, the phase of (pink) is strongly affected by the presence of , while the phase of (gray) is not much different from that in single Fano resonance, as can be seen in Fig. 8 (a-1), (a-2), and (a-3). Most importantly, a phase increase to of subradiant (dark brown) takes place when passing the resonance of superradiant (red), even though is also coupled to superradiant at the same time. As a result, the phase of (pink) is pulled down sharply near the resonance of (green), which leads to a delay in recovering back to its original value, when compared with what happens in a single Fano resonance. See blue and pink curves at the phase in Fig. 8(a-2) and (a-3). We also examine the amplitude evolution of the superradiant dipole (pink) in Fig. 8 (b-2) and (b-3). We find that the spectral position of amplitude maximum at the low frequency is red-shifted to a frequency below the resonance of superradiant dipole (red).

These two changes in phase and amplitude of , namely, a delay in phase recovery and a red-shift of amplitude maximum, are mainly responsible for the difference in extinction spectra (black curves) between single Fano resonances of Fig. 8 (c-1) & (c-2) and double Fano resonances of Fig. 8 (c-3). In Fig. 8 (c-3) we can identify how double Fano resonances take place by examining the extinction spectra of isolated (red), isolated (blue), and MM (black) where and are coupled to (green). In region 1⃝, both enhancement and red-shift of extinction in as seen in the change from red curve to black curve originate from the red-shift of amplitude maximum of dipole . In region 3⃝, on the other hand, enhancement and blue-shift of extinction in as seen in the change from blue curve to black curve is due to the delay in phase recovery of dipole . These two changes give rise to a Fano dip wider than that in a single Fano resonance in region 2⃝. These features are examples of coherent effects in double Fano resonances resulting from near-field couplings of two superradiant oscillators and to a common subradiant oscillator in a tripod metamaterial system. In the extinction spectra of Fig. 8 (c-3) the coherent effects show up as an energy transfer from low-energy oscillator to high-energy oscillator in DSRR, which was mentioned in the discussion of Fig. 7.

## Iv Electric field distributions in double Fano resonances

One advantage of the coupled equations of motion Eq. (2)-(4) in describing double Fano resonances is that two orthogonal dipolar oscillations and residing in DSRR of MM can be identified. On the other hand, an FDTD simulation provides the information on the electric field distribution of MM, which allows us to identify the electric multipolar nature of coherent coupling between DSRR and FRR in MM.

All the FDTD simulation is conducted with a commercial software FDTD solution Lumerical lumerical (). We plot distributions of -component of electric fields, which are shown in Fig. 9 for a normally incident wave with horizontal polarization (-axis). In a planar metamaterial, it is known that -components of the electric field, normal to the planar metamaterial surface ( plane), provide the information on planar charge distributions of the plamonic excitationsgu2012active (); shen2012photoexcited (); He2013371 (); Bitzer:11 (). First of all, we note that excitations in DSRR and FRR of the composite MM are dipolar, similar to dipolar double Fano resonances reported in Ref. artar2011directional (), where a single dipolar resonance is coupled two dipoles placed side-by-side, and double Fano resonances are observed in opposite ends of the spectrum with respect to a single dipolar resonance. In region 1⃝, dipolar excitation of oscillator is dominant in DSRR, and the dipoles in DSRR and FRR are anti-parallel forming a bonding hybridization. In region 3⃝, on the other hand, dipolar excitation of oscillator is dominant in DSRR, and the dipoles in DSRR and FRR are parallel forming an anti-bonding hybridization.

In region 2⃝, a Fano dip takes place at the frequency 1.6THz. In order to show explicitly the excitation of subradiant mode, distributions of z-component of electric fields in the isolated DSRR, the isolated FRR, and the composite MM are shown in Fig. 10. Comparison of Figs. 10(a) and (c) shows that the induced dipole of DSRR is significantly suppressed by a coherent coupling with FRR, giving rise to the Fano dip in the spectrum of DSRR in MM. Most importantly, the induced dipole of FRR in Fig. 10(c) is opposite to that of FRR in Fig. 10(b). While the FDTD simulation provides the electric field distribution, thereby, clarifies the multipolar nature of excitations of double Fano resonances, it cannot distinguish two orthogonal superradiant dipolar excitations in DSRR of MM or DSRR of MM.

By comparing Fano resonances dipole-quadrupole and DSRR-FRR structures, we clarify the Fano mode of dipole excitations in FRR. Figure 11 shows schematics of Fano mode excitations. The left panel represents the uncoupled (a) dipole-quadrupole structure and (c) DSRR-FRR structure. While a quadrupole is excited in two parallel rods in dipole-quadrupole structure, a weak dipole is excited in FRR in DSRR-FRR structure. A Fano dip in dipole-quadrupole structure comes from a suppression of dipole excitation owing to a coherent coupling with the quadrupole excitation in two parallel rods as shown in Fig. 11(b). In DSRR-FRR structure, on the other hand, a dipole excitation in subradiant FRR induced by a near-field coupling with DSRR suppresses a superradiant mode dipolar excitation in DSRR of MM as shown in Fig. 11(d).

We briefly examine how double Fano resonances behave as the relative orientation of DSRR and FRR changes by FDTD simulation. As shown in Fig. 12, MM with exhibits a deeper Fano dip at the spectral location of a lower energy than MM with , while MM with possesses double Fano resonances similar to MM with . We suppose that the highest electric field confinement takes place near edges of DSRR and FRR for MM with .

## V Summary

By embedding four-rod resonators inside double-split ring resonators superlattice, a composite metamaterial system is introduced, where a common subradiant oscillator in four-rod resonator is coupled with two superradiant oscillators in double-split ring resonator to give rise to double Fano resonances. As a classical analogue of four-level tripod atomic system, the extinction spectrum of the composite metamaterial exhibits a double Fano-based coherent effect. Double Fano resonances are correlated owing to the conjugation of two coupling strengths and , which shows up as a transfer of absorbed power from one superradiant dipole to the other superradiant dipole in double-split ring resonators. Double Fano resonances in a tripod metamaterial system can be utilized in designing plasmonic resonances of metamaterial to engineer coherent effects in metamaterials.

## Acknowledgments

This work is supported by Quantum Metamaterial Research Center program and Global Frontier Program of Center for Advanced Meta-Materials (Ministry of Science, ICT and Future Planning, National Research Foundation, Republic of Korea).

## References

- [1] Niels Verellen, Yannick Sonnefraud, Heidar Sobhani, Feng Hao, Victor V Moshchalkov, Pol Van Dorpe, Peter Nordlander, and Stefan A Maier. Fano resonances in individual coherent plasmonic nanocavities. Nano letters, 9(4):1663–1667, 2009.
- [2] Andrey E. Miroshnichenko, Sergej Flach, and Yuri S. Kivshar. Fano resonances in nanoscale structures. Rev. Mod. Phys., 82(3):2257–2298, Aug 2010.
- [3] Boris Luk’yanchuk, Nikolay I Zheludev, Stefan A Maier, Naomi J Halas, Peter Nordlander, Harald Giessen, and Chong Tow Chong. The fano resonance in plasmonic nanostructures and metamaterials. Nature materials, 9(9):707–715, 2010.
- [4] Mohsen Rahmani, Boris Luk’yanchuk, and Minghui Hong. Fano resonance in novel plasmonic nanostructures. Laser & Photonics Reviews, 7(3):329–349, 2013.
- [5] Shuang Zhang, Dentcho A. Genov, Yuan Wang, Ming Liu, and Xiang Zhang. Plasmon-induced transparency in metamaterials. Phys. Rev. Lett., 101(4):047401, Jul 2008.
- [6] Na Liu, Lutz Langguth, Thomas Weiss, Jurgen Kastel, Michael Fleischhauer, Tilman Pfau, and Harald Giessen. Plasmonic analogue of electromagnetically induced transparency at the drude damping limit. Nat. Mater., 8(9):758–762, 2009.
- [7] N. Liu, T. Weiss, M. Mesch, L. Langguth, U. Eigenthaler, M. Hirscher, C. Sonnichsen, and H. Giessen. Planar metamaterial analogue of electromagnetically induced transparency for plasmonic sensing. Nano Letters, 10(4):1103–1107, 2010.
- [8] Xiaojun Liu, Jianqiang Gu, Ranjan Singh, Yingfang Ma, Jun Zhu, Zhen Tian, Mingxia He, Jiaguang Han, and Weili Zhang. Electromagnetically induced transparency in terahertz plasmonic metamaterials via dual excitation pathways of the dark mode. Applied Physics Letters, 100(13):131101–131101, 2012.
- [9] Junqiao Wang, Chunzhen Fan, Jinna He, Pei Ding, Erjun Liang, and Qianzhong Xue. Double fano resonances due to interplay of electric and magnetic plasmon modes in planar plasmonic structure with high sensing sensitivity. Optics express, 21(2):2236–2244, 2013.
- [10] Shao-Ding Liu, Zhi Yang, Rui-Ping Liu, and Xiu-Yan Li. Multiple fano resonances in plasmonic heptamer clusters composed of split nanorings. Acs Nano, 6(7):6260–6271, 2012.
- [11] Shao-Ding Liu, Ming-Jiang Zhang, Wen-Jie Wang, and Yun-Cai Wang. Tuning multiple fano resonances in plasmonic pentamer clusters. Applied Physics Letters, 102(13):133105–133105, 2013.
- [12] Na Liu, Mario Hentschel, Thomas Weiss, A Paul Alivisatos, and Harald Giessen. Three-dimensional plasmon rulers. Science, 332(6036):1407–1410, 2011.
- [13] Timothy J Davis, Mario Hentschel, Na Liu, and Harald Giessen. Analytical model of the three-dimensional plasmonic ruler. ACS nano, 6(2):1291–1298, 2012.
- [14] Alp Artar, Ahmet A Yanik, and Hatice Altug. Multispectral plasmon induced transparency in coupled meta-atoms. Nano letters, 11(4):1685–1689, 2011.
- [15] S Rebić, D Vitali, C Ottaviani, P Tombesi, M Artoni, F Cataliotti, and R Corbalan. Polarization phase gate with a tripod atomic system. Physical Review A, 70(3):032317, 2004.
- [16] Shujing Li, Xudong Yang, Xuemin Cao, Chunhong Zhang, Changde Xie, and Hai Wang. Enhanced cross-phase modulation based on a double electromagnetically induced transparency in a four-level tripod atomic system. Physical review letters, 101(7):073602, 2008.
- [17] Andrew MacRae, Geoff Campbell, and AI Lvovsky. Matched slow pulses using double electromagnetically induced transparency. Optics letters, 33(22):2659–2661, 2008.
- [18] Yong He and Ka-Di Zhu. Fano correlation effect of optical response due to plasmon–exciton–plasmon interaction in an artificial hybrid molecule system. JOSA B, 30(4):868–873, 2013.
- [19] Zhengyang Bai, Chao Hang, and Guoxiang Huang. Classical analogs of double electromagnetically induced transparency. Optics Communications, 291:253–258, 2013.
- [20] Boyoung Kang, JH Woo, E Choi, Hyun-Hee Lee, ES Kim, J Kim, Tae-Jong Hwang, Young-Soon Park, DH Kim, and JW Wu. Optical switching of near infrared light transmission in metamaterial-liquid crystal cell structure. Optics Express, 18(16):16492–16498, 2010.
- [21] J.H. Woo, E.S. Kim, E. Choi, Boyoung Kang, Hyun-Hee Lee, J. Kim, Y.U. Lee, Tae Y. Hong, Jae H. Kim, and J.W. Wu. Cryogenic temperature measurement of thz meta-resonance in symmetric metamaterial superlattice. Opt. Express, 19:4384–4392, 2011.
- [22] YU Lee, EY Choi, JH Woo, ES Kim, and JW Wu. Reflection resonance switching in metamaterial twisted nematics cell. Optics express, 21(15):17492–17497, 2013.
- [23] JB Pendry, AJ Holden, DJ Robbins, and WJ Stewart. Low frequency plasmons in thin-wire structures. Journal of Physics: Condensed Matter, 10(22):4785, 1998.
- [24] John B Pendry, A J_ Holden, DJ Robbins, and WJ Stewart. Magnetism from conductors and enhanced nonlinear phenomena. Microwave Theory and Techniques, IEEE Transactions on, 47(11):2075–2084, 1999.
- [25] JB Pendry, L Martin-Moreno, and FJ Garcia-Vidal. Mimicking surface plasmons with structured surfaces. Science, 305(5685):847–848, 2004.
- [26] Yingfang Ma, Zhongyang Li, Yuanmu Yang, Ran Huang, Ranjan Singh, Shuang Zhang, Jianqiang Gu, Zhen Tian, Jiaguang Han, and Weili Zhang. Plasmon-induced transparency in twisted fano terahertz metamaterials. Optical Materials Express, 1(3):391–399, 2011.
- [27] J Wallauer and M Walther. Fano line shape and phase reversal in a split-ring resonator based metamaterial. Physical Review B, 88(19):195118, 2013.
- [28] See http://www.teraview.com.
- [29] See http://www.lumerical.com.
- [30] Hou-Tong Chen, Willie J Padilla, Joshua MO Zide, Arthur C Gossard, Antoinette J Taylor, and Richard D Averitt. Active terahertz metamaterial devices. Nature, 444(7119):597–600, 2006.
- [31] See http://docs.lumerical.com/en/metamaterialsthzchen.html.
- [32] Robert C Weast, Melvin J Astle, and William H Beyer. CRC handbook of chemistry and physics, volume 69. CRC press Boca Raton, FL, 1988.
- [33] Yong S Joe, Arkady M Satanin, and Chang Sub Kim. Classical analogy of fano resonances. Phys. Scr, 74:259–266, 2006.
- [34] Martin Frimmer, Toon Coenen, and A. Femius Koenderink. Signature of a fano resonance in a plasmonic metamolecule’s local density of optical states. Phys. Rev. Lett., 108:077404, Feb 2012.
- [35] Jianqiang Gu, Ranjan Singh, Xiaojun Liu, Xueqian Zhang, Yingfang Ma, Shuang Zhang, Stefan A Maier, Zhen Tian, Abul K Azad, Hou-Tong Chen, et al. Active control of electromagnetically induced transparency analogue in terahertz metamaterials. Nature communications, 3:1151, 2012.
- [36] Xiaopeng Shen and Tie Jun Cui. Photoexcited broadband redshift switch and strength modulation of terahertz metamaterial absorber. Journal of Optics, 14(11):114012, 2012.
- [37] Xun jun He, Jian min Wang, Xiao hua Tian, Jiu xing Jiang, and Zhao xin Geng. Dual-spectral plasmon electromagnetically induced transparency in planar metamaterials based on bright-dark coupling. Optics Communications, 291(0):371 – 375, 2013.
- [38] Andreas Bitzer, Alex Ortner, Hannes Merbold, Thomas Feurer, and Markus Walther. Terahertz near-field microscopy of complementary planar metamaterials: Babinet’s principle. Opt. Express, 19(3):2537–2545, Jan 2011.
- [39] Alp Artar, Ahmet Ali Yanik, and Hatice Altug. Directional double fano resonances in plasmonic hetero-oligomers. Nano letters, 11(9):3694–3700, 2011.
- [40] U. Fano. Effects of configuration interaction on intensities and phase shifts. Phys. Rev., 124:1866–1878, Dec 1961.
- [41] Benjamin Gallinet and Olivier J. F. Martin. Ab initio theory of fano resonances in plasmonic nanostructures and metamaterials. Phys. Rev. B, 83:235427, 2011.
- [42] Benjamin Gallinet and Olivier J. F. Martin. Relation between near–field and far–field properties of plasmonic fano resonances. Opt. Express, 19(22):22167–22175, Oct 2011.
- [43] CL Garrido Alzar, MAG Martinez, and P Nussenzveig. Classical analog of electromagnetically induced transparency. American Journal of Physics, 70:37, 2002.
- [44] B Gallinet, T Siegfried, H Sigg, P Nordlander, OJ Martin, et al. Plasmonic radiance: probing structure at the angstrom scale with visible light. Nano letters, 13(2):497–503, 2013.
- [45] F López-Tejeira, R Paniagua-Domínguez, R Rodríguez-Oliveros, and JA Sánchez-Gil. Fano-like interference of plasmon resonances at a single rod-shaped nanoantenna. New Journal of Physics, 14(2):023035, 2012.
- [46] Feng Hao, Yannick Sonnefraud, Pol Van Dorpe, Stefan A Maier, Naomi J Halas, and Peter Nordlander. Symmetry breaking in plasmonic nanocavities: subradiant lspr sensing and a tunable fano resonance. Nano letters, 8(11):3983–3988, 2008.
- [47] Feng Hao, Elin M Larsson, Tamer A Ali, Duncan S Sutherland, and Peter Nordlander. Shedding light on dark plasmons in gold nanorings. Chemical Physics Letters, 458(4):262–266, 2008.
- [48] Benjamin Gallinet and Olivier JF Martin. Influence of electromagnetic interactions on the line shape of plasmonic fano resonances. ACS nano, 5(11):8999–9008, 2011.
- [49] Andrea Lovera, Benjamin Gallinet, and Peter Nordlander. Mechanisms of fano resonances in coupled plasmonic systems. ACS Nano, 10, 2013.

## Appendix: Fano resonance in two coupled oscillators driven simultaneously

We show that the characteristic asymmetry feature of Fano resonance is retained in plasmonic structure when both superradiant and subradiant oscillators are externally driven.

Originally Fano resonance is introduced to explain an asymmetric peak in excitation spectra observed near autoionization discrete energy level overlapping with continuum energy level of helium.[40] While the discrete state (dark mode) is not directly excited, a rapid phase change takes place in the phase of a slowly varying-phase background continuum (bright mode) near the discrete resonance. Hence in-phase and out-of-phase interferences take place between continuum and discrete resonances, giving rise to an asymmetric line shape in the continuum spectrum. The sequence of constructive and destructive interferences depends on whether the discrete resonance falls on the red or blue side of the background continuum resonance. The asymmetric Fano line shape is characterized by 3 parameters, the resonance shift , the resonance width , and the asymmetry parameter .

Differently from a dark mode in atomic spectroscopy, there exits an intrinsic non-radiative Joule loss in a plamonic dark mode when electromagnetically driven, and Gallinet et al. [41] obtained a general formula of Fano-like asymmetric line shape in metallic plasmonic structures, where a modulation damping parameter is introduced, a measure of the decrease in Fano resonance contrast [42]. It is also shown that a classical two coupled oscillators model [43] with a lossy dark oscillator provides the same general formula when the complex amplitude of the bright oscillator is expanded near the dark resonance [44].

Furthermore, even without resorting to a plasmonic dark mode, a far-field Fano-like asymmetric line shape is observed in plasmonic structures. The occurrence of a Fano-like interference in a single nano-rod plasmonic structures is attributed to a spatial overlap of two non-orthogonal adjacent plasmonic modes [45]. And in a nonconcentric ring/disk cavity plasmonic structure, Fano line shape still persists when dark multipolar resonances are optically excited by a retardation effect at a non-normal incidence angle [46, 47]. In a heptamer plasmonic nanostructure composed of one center disc with surrounding hexamer ring, an examination of local density of optical states by a selective e-beam excitation of center and ring showed that a far-field interference between the dressed superradiant and subradiant eigenmodes fully describes Fano resonance [34]. In fact, asymmetric Fano line shapes observed in a variety of plasmonic structures are fitted to the general Fano formula with an excellent agreement [48]. These studies lead to the notion that in plasmonic structures an interference between superradiant and subradient eigenmodes overlapping both spectrally and spatially gives rise to an asymmetric Fano line shape resonance [45, 49].

Now we raise the question whether asymmetric Fano line shape is retained even when both superradiant and subradiant oscillators are driven simultaneously in the presence of a coherent coupling. We examine how the broad Lorentzian spectrum of a superradiant oscillator is modified when coupled with a driven subradiant narrow oscillator by adopting a simple mechanical two coupled oscillators model. Two coupled oscillators are described by a set of equations of motion.

(5) | |||||

(6) |

Narrow spectrum of the subradiant oscillator is located within broad spectrum of the superradiant oscillator . Putting , we obtain the amplitude of the superradiant oscillator near resonance of narrow subradiant oscillator .

where and . An external driving of the subradiant oscillator contributes to the superradiant amplitude via a coupling . Following procedures in Ref. [44], the superradiant oscillator spectrum has the following relation.

where

Formula of asymmetric Fano resonance line shape reported in Ref. [44] is kept the same as a closed form with a modification in the asymmetry parameter from to .

We observe that a dark mode is not necessarily prerequisite for an asymmetric Fano resonance to take place in plasmonic structures. The dark mode is simply an example of non-driven narrow resonance. The characteristic asymmetry feature of Fano resonance in plasmonic structures is retained even in the presence of an external driving of narrow oscillator as far as a dispersive coherent coupling between two oscillators ensures a -phase shift in the broad superradiant oscillator near the narrow subradiant oscillator resonance.