Abstract
A blueshift tunable metamaterial (MM) exhibiting a doublenegative refractive index based on a topological insulator (bismuth selenide, Bi_{2}Se_{3}) has been demonstrated in the nearinfrared (NIR) spectral region. The potential of Bi_{2}Se_{3} as a dielectric interlayer of the multilayer MM is explored. The optical response of elliptical nanohole arrays penetrating through Au/Bi_{2}Se_{3}/Au films is numerically investigated using the finite difference time domain (FDTD) method. The blueshift tuning range of the MM is as high as 370 nm (from 2,140 to 1,770 nm) after switching the Bi_{2}Se_{3} between its trigonal and orthorhombic states.
Background
Metamaterials (MMs) are artificially engineered composites that attract considerable interests due to their exceptional electromagnetic properties, which are not typically found in nature, such as negative refractive index and cloaking [14]. These MMs with various subwavelength resonant elements have offered magnetic and/or electric resonant responses to incident electromagnetic radiation, scalable from the microwave frequencies up to the terahertz and optical ones [57]. Particularly, nanohole resonators embedded in metaldielectricmetal (MDM) multilayers are frequently used as building blocks of negativerefractiveindex MMs [811], owing to the coupling between surface plasmons counterpropagating on the two closely spaced interfaces which results in a closed loop of the electric currents. This gives rise to magnetic dipolar resonances between the two coupled metal layers, while the continuous metallic strip parts provide the electric resonance moments [12,13]. All these features make the nanohole array perforating through MDM films become a strong candidate for developing threedimensional negativeindex MMs [14,15].
One of the obstacles in this progress is the resonance responses of MMs to the impinge light which are usually fixed once the dimension of the structure is determined, thus making the MMs possess a limited bandwidth. However, for many applications (switching, modulation, filtering, etc.), it would be highly desirable to tune the MM resonances over a wide bandwidth. To this end, tunable photonic MMs, the spectral range of which can be controlled by changing the dielectric environment of the resonator with liquid crystals (LCs) [1618]; phase transition materials [19,20]; and optical pumping [21,22] have been discussed recently. However, the challenge is to develop tunable MDMMMs in the nearinfrared (NIR) regime. It is due to the fact that frequency tunability of the MDMMM mainly requires for the interlayer dielectric material to possess a tunable effective dielectric constant in the NIR region, hence limiting the choice of the active materials. Here, we take a different approach to actively tune the resonant frequency of the MDMMMs in the NIR regions by using bismuth selenide (Bi_{2}Se_{3}) as the dielectric layer.
Recently, a rising Dirac material  topological insulators (TIs)  had been intensively researched in condensed matter physics [23,24]. In analogy to the optoelectronic applications of graphene, a thin layer of TIs has been theoretically predicted to be a promising material for broadband and highperformance optoelectronic devices such as photodetectors, terahertz lasers, waveguides, and transparent electrodes [25]. Among these TIs, Bi_{2}Se_{3} is a particularly interesting compound due to its relatively large bulk band gap and a simple surface state consisting of a single Dirac conelike structure [26,27]. Study of the dielectric function reveals that the optical dielectric constant of Bi_{2}Se_{3} can be very different for the trigonal and orthorhombic phases in the NIR regime [28]. Bi_{2}Se_{3} exhibits a number of means through which their dielectric properties can be altered [2833]. Herein, structural phase transition between trigonal and orthorhombic states, which is achieved by a high pressure and temperature, is proposed and studied as a means to change the intrinsic effective dielectric properties of the MDMMMs [28].
Here, we numerically demonstrate a blueshift tunable nanometerscale MM consisting of an elliptical nanohole array (ENA) embedded in the MDM multilayers where the dielectric core layer is a Bi_{2}Se_{3} composite. Under a high pressure of 2 to 4.3 Pa at 500°C, Bi_{2}Se_{3} occurring in trigonal phase undergoes a transition to orthorhombic phase and features a large change in the values of the effective dielectric constant [28]. Accordingly, a massive blueshift of the resonant response (from 2,140 to 1,770 nm) of a Bi_{2}Se_{3}based MDMENA is achieved in the NIR region. Our proposed blueshift tunable negativeindex MM provides greater flexibility in the practical application and has a potential of enabling efficient switches and modulators in the NIR region.
Methods
The proposed MDMENA suspended in a vacuum is shown in Figure 1, with the coordinate axes and the polarization configuration of the normally incident light. The structure consists of trilayers of Au/Bi_{2}Se_{3}/Au. The thickness of each Au layer is 30 nm, and the thickness of the Bi_{2}Se_{3} layer is 60 nm. The metamaterial parameters are optimized for the maximum sensitivity of the resonance to a change in the refractive index of the Bi_{2}Se_{3} core dielectric layer in the NIR spectral range. The element resonator is shown in Figure 1b, where the pitch of the elliptical holes is L = 400 nm, the diameters of the elliptical holes are d_{1} = 240 nm and d_{2} = 120 nm, and β is a crosssectional plane of the structure. The zaxis is normal to the structure surface, and the xy plane is parallel to the structure surface. This simulated structure is periodically extended along the x and y axes. The tunable optical properties of the structure are calculated using 3D EM Explorer Studio [34], a commercial finite difference time domain (FDTD) code. In the simulation, a simple Drudetype model for Au permittivity was used, which is a good approximation to experimental values in the NIR region. The dielectric properties of Au as given by Johnson and Christy are used [35]. A plane wave source is simulated at normal incidence to the structure. The computational domain (400 nm × 400 nm × 1,000 nm) has a perfectly matched layer (PML), absorbing boundaries in the z direction and periodic boundaries in the xy plane [36]. A uniform FDTD mesh size is adopted. The mesh size is the same along all Cartesian axes: ∆x = ∆y = ∆z = 2 nm, which is sufficient to minimize the numerical errors arising from the FDTD method.
Figure 1. Schematic of the proposed structure. (a) Schematic of the MDM structure consisting of a 60nmthick Bi_{2}Se_{3} dielectric layer between two 30nmthick Au films perforated with a square array of elliptical holes suspended in air. The lattice constant is L = 400 nm, and hole diameters are d_{1} = 240 nm and d_{2} = 120 nm. (b) Illustration of the square lattice of ENA.
The topological insulator material Bi_{2}Se_{3} was selected due to its significantly different optical properties between the trigonal and orthorhombic phases. The real (ϵ_{1}) and imaginary (ϵ_{2}) parts of the dielectric function for the different structural phases of Bi_{2}Se_{3} were obtained from the published data in [28]; the NIR spectral region is shown in Figure 2. A large change in the dielectric function across the NIR is obtained after switching Bi_{2}Se_{3} from trigonal to its orthorhombic phase.
Figure 2. Dielectric constant of the Bi_{2}Se_{3}. (a) Real part of dielectric function ϵ_{1}(ω) for trigonal and orthorhombic phases of Bi_{2}Se_{3}. (b) Imaginary part of dielectric function ϵ_{2}(ω) for trigonal and orthorhombic phases of Bi_{2}Se_{3}.
After the complex coefficients of transmission and reflection are obtained by the 3D EM Explorer Studio, in which T_{a} is the amplitude and φ_{a} is the phase of the transmission coefficient, and R_{a} is the amplitude and φ_{ra} is the phase of the reflection coefficient, the effective optical parameters can be extracted using the Fresnel formula [37].
For an equivalent isotropic homogenous slab of thickness h surrounded by semiinfinite media with refractive index n_{1} and n_{3} under normal incidence, we have
The socalled material parameters ϵ_{eff} and μ_{eff} of a single layer of a doublefishnet negativeindex metamaterial are extracted using the wellknown NicholsonRossWeir (NRW) method [3840]. Therefore, once n_{eff} and η are evaluated, the effective permittivity and permeability are calculated using
where n_{eff} is the effective refractive index, η is the impedance, h is the thickness of the structure, k = ω/c, c is the speed of light, m is an arbitrary integer, and n_{1} = n_{3} = 1 since the structure is suspended in a vacuum. The signs of n_{eff} and η and the value of m are resolved by the passivity of the metamaterial that requires the signs of the real part of impedance η and imaginary part of effective index n_{eff} to be positive, i.e., Real(η) > 0, Imag(n_{eff}) > 0 which is consistent with the study described in [39,40]. We then apply this extraction approach to determine the change in the optical response of the structure when the phase of Bi_{2}Se_{3} is switched between its trigonal and orthorhombic states.
Results and discussion
The ENA has a lower transmittance for spolarized light due to the electric field's orientation with respect to the metallic stripe width [12]; hence, the polarization of the incident wave was set to be ppolarized. As shown in Figure 1a, s polarization means that the incident electric field vector is parallel to the long axis of the ENA, and the incident electric field vector perpendicular to the long axis of the ENA is then denoted by p polarization. We first investigate the transmittance T = t^{2} and reflectance R = r^{2} of the structure for p polarization in Figure 3. Structures with a different dielectric constant of Bi_{2}Se_{3} (shown in Figure 2) were modeled to investigate the effect of the phase change of Bi_{2}Se_{3} on the position and amplitude of the spectrums. It can be seen that the resonance wavelength blueshifts from 2,140 to 1,770 nm when the structural phase of Bi_{2}Se_{3} switches from trigonal to orthorhombic. The structure is impedancematched, hence possessing a low reflectance corresponding to the dips in reflectance of Figure 3b for different forms of Bi_{2}Se_{3}.
Figure 3. Transmittance and reflectance. 3D FDTD simulation of (a) spectrum of transmittance and (b) spectrum of reflectance, for the different phases of the Bi_{2}Se_{3} dielectric layer, where the light source is p polarization at normal incidence angle.
In Figure 4, the transmission (t) and reflection(r) phases are demonstrated. The transmission phase exhibits a dip around the resonance, indicating that the light is advanced in phase at the resonance, characteristic of a lefthanded material [41]. Importantly, changing the structural phase of the Bi_{2}Se_{3} offers transmission and reflection phase tunability which implies tunable effective constitutive parameters in the structure.
Figure 4. Transmission and reflection phase. 3D FDTD simulation of (a) phase of transmission and (b) phase of reflection, for the different phases of the Bi_{2}Se_{3} dielectric layer, where the light source is p polarization at normal incidence angle.
Taking into account the subwavelength thickness of the structure, the extracted n_{eff} can be retrieved from the transmission and reflection coefficients shown in Figure 5. For the MM with the trigonal Bi_{2}Se_{3} dielectric layer, the negativeindex band extends from 1,880 to 2,420 nm with a minimum value of the real part of the refractive index Real(n_{eff}) = 7. Regarding losses, the figure of merit (FOM) defined as is taken to show the overall performance of the MM, where Imag(n_{eff}) is the imaginary part of the refractive index. As shown in Figure 5c, the FOM for the trigonal phase is 2.7 at the operating wavelength of 2,080 nm. The negativeindex band of the orthorhombic Bi_{2}Se_{3}based MM extends from 1,600 to 2,214 nm having a minimum value of Real(n_{eff}) = 3.2. The FOM is 1.2 at the resonant wavelength of 1,756 nm. Furthermore, the bandwidth of Real(n_{eff}) becomes wider for the orthorhombic Bi_{2}Se_{3} film in Figure 5a due to increased damping of the plasmon resonance but at the cost of an accompanied lower value of FOM in Figure 5c.
Figure 5. Effective index and figure of merit. 3D FDTD simulation of (a) real part of n_{eff}, (b) imaginary part of n_{eff}, and (c) figure of merit for the different phases of the Bi_{2}Se_{3} dielectric layer, where the light source is p polarization at normal incidence angle.
The refractive index is expressed in terms of the real and imaginary parts of the permeability μ_{eff} and permittivity ϵ_{eff}. However, the sign of the real part of the permeability μ_{eff}: Real(μ_{eff}) determines the relative magnitudes of the imaginary and real parts of the refractive index [41]. To achieve a negative index with a small loss, a negative Real(μ_{eff}) is required. Therefore, we have simulated μ_{eff} and ϵ_{eff} for the structure as shown in Figure 6. For the trigonal and orthorhombic phases of Bi_{2}Se_{3}, Real(μ_{eff}) has a Fanotype line shape and Im(μ_{eff}) has a Lorentzian line shape in the region of the negative index. Moreover, a doublenegative MM can be achieved when Real(μ_{eff}) and Real(ϵ_{eff}) simultaneously reach negative values over a wide frequency range and thus a reduced loss. The maximum negative Real(μ_{eff}) decreases with the phase transition from the trigonal to orthorhombic, hence resulting in the smaller value of the maximum negative Real(n_{eff}) at the orthorhombic phase.
Figure 6. Permittivity and permeability. 3D FDTD simulation of (a) the real part of μ_{eff}, (b) the imaginary part of μ_{eff}, (c) the real part of ϵ_{eff}, and (d) the imaginary part of ϵ_{eff} for the different phases of the Bi_{2}Se_{3} dielectric layer, where the light source is p polarization at normal incidence angle.
This magnetic negative response can be explained looking at the current and field distribution at the resonance wavelengths. Figure 7 shows the current and total magnetic field intensity for the magnetic resonant wavelengths of 2,140 and 1,770 nm at the β plane shown in Figure 1. In the field maps of Figure 7, the arrows show the currents, whereas the color shows the intensity of the magnetic field. It clearly shows that the antiparallel currents are excited at opposite internal metallic interfaces, closed by an electric displacement current J_{D}. Therefore, these virtual current loops between two Au layers on the β plane give rise to magnetic resonant responses of negative Re(μ_{eff}) that interact strongly with the incident magnetic field at which the total magnetic field intensity H is strongly localized in the Bi_{2}Se_{3} dielectric layer between the top and bottom Au layers [42].
Figure 7. Magnetic field intensity and displacement current. A map of the normalized total magnetic field intensity distribution (color bar) and displacement current (J_{D}) (arrows) along the plane β: (a) at a 2,140nm resonance wavelength for trigonal Bi_{2}Se_{3} and (b) at a 1,770nm resonance wavelength for orthorhombic Bi_{2}Se_{3}, where the light source is p polarization at normal incidence angle.
Specifically, H for the orthorhombic phase shown in Figure 7b is weaker than the trigonal phase shown in Figure 7a. It depicts that the MM based on orthorhombic phase has a smaller magnetic dipolar moment than the trigonal phase and thus smaller FOM.
To further understand the negativeindex resonance in the metamaterials, it is useful to study the dispersion of the surface plasmon polariton (SPP) modes within the multilayer structure. Both the internal and external SPP modes in the multilayer metamaterials are similar to those of the same structure without resonant elements, i.e., MDM films [42], where the internal SPP mode resonates in the inner surfaces of the metal layers and the external SPP mode resonates in the outer surfaces of the metal layers. Therefore, the SPP dispersion relation of the multilayer metamaterial can be approximately approached by that of the MDM structure. In Figure 8, we have calculated the SPP mode dispersion relation of the AuBi_{2}Se_{3}Au sheets with the top Au film thickness t_{1} = 30 nm, middle Bi_{2}Se_{3} film thickness t_{2} = 60 nm, and bottom Au film thickness t_{3} = 30 nm. The transmittance spectrum of the multilayer metamaterials is also depicted together with the dispersion relation of the AuBi_{2}Se_{3}Au films.
Figure 8. Dispersion relation of the structure. Representation of the dispersion relation of the AuBi_{2}Se_{3}Au trilayer (left) and the transmittance of the multilayer metamaterials (right) for both (a) trigonal Bi_{2}Se_{3} and (b) orthorhombic Bi_{2}Se_{3}.
Recalling the coupling condition from light to SPP modes [42], it can be seen that the (1,1) internal resonance of the AuBi_{2}Se_{3}Au trilayer is excited at 2,350 nm associated with the trigonal Bi_{2}Se_{3} in Figure 8a. This internal SPP resonance blueshifts to 2,010 nm when the trigonal state changes to the orthorhombic state as shown in Figure 8b. We also observe that the two internal (1,1) modes which appear at 2,350 and 2,010 nm in the simple MDM structure do not perfectly match the two absorbance peaks at the resonance wavelengths of 2,140 and 1,770 nm in the multilayer metamaterials for both the trigonal and orthorhombic phases, respectively. This difference is because the dispersion relation of the SPP modes used as matching condition does not include the resonant squares, which cause a resonance shift [42].
Conclusions
In conclusion, this work numerically demonstrates the tunable optical properties of an ENA perforated through Au/Bi_{2}Se_{3}/Au trilayers. We present that the MDMENA can be improved to exhibit a substantial frequency tunability of the intrinsic resonance in the NIR spectral region by selecting Bi_{2}Se_{3} as the active dielectric material. Particularly, the resonant transmission, reflection, and the effective constitutive parameters of the Bi_{2}Se_{3}coupled multilayer MM can be massively blueshifted by transiting the phase of the Bi_{2}Se_{3} film from the trigonal to orthorhombic. This may offer an innovative and practical paradigm for the development of tunable photonic devices. We expect that our results will facilitate further experimental studies of the tunable MMs and make this technique suitable for tuning of plasmon resonance in the optical regime.
Competing interests
The authors declare that they have no competing interests.
Authors' contributions
TC conceived the idea of using topological insulator for tuning the resonance in the metamaterials, designed the metamaterial, and wrote the manuscript. SW carried out the simulations and prepared the figures. Both authors read and approved the final manuscript.
Acknowledgements
We acknowledge the financial support from National Natural Science Foundation of China (grant nos. 61172059, 51302026), PhD Programs Foundation of the Ministry of Education of China (grant no. 20110041120015), Postdoctoral Gathering Project of Liaoning Province (grant no. 2011921008), and The Fundamental Research for the Central University (grant no. DUT12JB01).
References

Qiu CW, Gao L: Resonant light scattering by small coated nonmagnetic spheres: magnetic resonances, negative refraction and prediction.

Shalaev VM: Optical negativeindex metamaterials.
Nat Photonics 2007, 1:4148. Publisher Full Text

Soukoulis CM, Wegener M: Past achievements and future challenges in the development of threedimensional photonic metamaterials.

Zhou S, Huang X, Li Q, Xie YM: A study of shape optimization on the metallic nanoparticles for thinfilm solar cells.
Nanoscale Res Lett 2013, 8:447. PubMed Abstract  BioMed Central Full Text

Liaw JW, Chen HC, Kuo MK: Plasmonic Fano resonance and dip of AuSiO_{2}Au nanomatryoshka.
Nanoscale Res Lett 2013, 8:468. PubMed Abstract  BioMed Central Full Text  PubMed Central Full Text

Zhang S, Fan W, Panoiu NC, Malloy KJ, Osgood RM, Brueck SRJ: Experimental demonstration of nearinfrared negativeindex metamaterials.
Phys Rev Lett 2005, 95:137404. PubMed Abstract  Publisher Full Text

Li T, Li JQ, Wang FM, Wang QJ, Liu H, Zhu SN, Zhu YY: Exploring magnetic plasmon polaritons in optical transmission through hole arrays perforated in trilayer structures.
Appl Phys Lett 2007, 90:251112. Publisher Full Text

Valentine J, Zhang S, Zentgraf T, UlinAvila E, Genov DA, Bartal G, Zhang X: Threedimensional optical metamaterial with a negative refractive index.
Nature 2008, 455:376379. PubMed Abstract  Publisher Full Text

Minovich A, Neshev DN, Powell DA, Shadrivov IV, Lapine M, Hattori HT, Tan HH, Jagadish C, Kivshar YS: Tilted response of fishnet metamaterials at nearinfrared optical wavelengths.

Zhang S, Fan W, Panoiu NC, Malloy KJ, Osgood RM, Brueck SRJ: Demonstration of metaldielectric negativeindex metamaterials with improved performance at optical frequencies.
J Opt Soc Am B 2006, 23:434438. Publisher Full Text

Cao T, Cryan MJ: Study of incident angle dependence for dualband double negativeindex material using elliptical nanohole arrays.

Pendry JB, Holden A, Robbins D, Stewart W: Magnetism from conductors and enhanced nonlinear phenomena.
IEEE Trans Microw Theory Tech 1999, 47(11):20752084. Publisher Full Text

Smith DR, Padilla WJ, Vier DC, NematNasser SC, Schultz S: Composite medium with simultaneously negative permeability and permittivity.
Phys Rev Lett 2000, 84:41844187. PubMed Abstract  Publisher Full Text

Zhao Q, Kang L, Du B, Li B, Zhou J, Tang H, Liang X, Zhang B: Electrically tunable negative permeability metamaterials based on nematic liquid crystals.
Appl Phys Lett 2007, 90:011112. Publisher Full Text

Wang X, Kwon DH, Werner DH, Khoo IC, Kildishev AV, Shalaev VM: Tunable optical negativeindex metamaterials employing anisotropic liquid crystals.
Appl Phys Lett 2007, 91:143122. Publisher Full Text

Minovich A, Neshev DN, Powell DA, Shadrivov IV, Kivshar YS: Tunable fishnet metamaterials infiltrated by liquid crystals.
Appl Phys Lett 2010, 96:193103. Publisher Full Text

Dicken MJ, Aydin K, Pryce IM, Sweatlock LA, Boyd EM, Walavalkar S, Ma J, Atwater HA: Frequency tunable nearinfrared metamaterials based on VO_{2} phase transition.
Opt Express 2009, 17:1833018339. PubMed Abstract  Publisher Full Text

Driscoll T, Kim HT, Chae BG, Kim BJ, Lee YW, Jokerst NM, Smith DR, Ventra MD, Basov DN: Memory metamaterials.
Science 2009, 325:15181521. PubMed Abstract  Publisher Full Text

Chen HT, O'Hara JF, Azad AK, Taylor AJ, Averitt RD, Shrekenhamer DB, Padilla WJ: Experimental demonstration of frequencyagile terahertz metamaterials.
Nat Photon 2008, 2:295298. Publisher Full Text

Hu XY, Zhang YB, Fu YL, Yang H, Gong QH: Lowpower and ultrafast alloptical tunable nanometerscale photonic metamaterials.
Adv Mater 2011, 23:42954300. PubMed Abstract  Publisher Full Text

Hasan MZ, Kane CL: Topological insulators.
Rev Mod Phys 2010, 82:3045. Publisher Full Text

Qi XY, Zhang SC: Topological insulators and superconductors.
Rev Mod Phys 2011, 83:1057. Publisher Full Text

Zhang X, Wang J, Zhang SC: Topological insulators for highperformance terahertz to infrared applications.

Hsieh D, Xia Y, Qian D, Wray L, Dil JH, Meier F, Osterwalder J, Patthey L, Checkelsky JG, Ong NP, Fedorov AV, Lin H, Bansil A, Grauer D, Hor YS, Cava RJ, Hasan MZ: A tunable topological insulator in the spin helical Dirac transport regime.
Nature 2009, 460:1101. PubMed Abstract  Publisher Full Text

Pan ZH, Vescovo E, Fedorov AV, Gardner D, Lee YS, Chu S, Gu GD, Valla T: Electronic structure of the topological insulator Bi_{2}Se_{3} using angleresolved photoemission spectroscopy: evidence for a nearly full surface spin polarization.
Phys Rev Lett 2011, 106:257004. PubMed Abstract  Publisher Full Text

Sharma Y, Srivastava P: Firstprinciples study of electronic and optical properties of Bi_{2}Se_{3} in its trigonal and orthorhombic phases.

Shao LH, Ruther M, Linden S, Essig S, Busch K, Weissmüller J, Wegener M: Electrochemical modulation of photonic metamaterials.
Adv Mater 2010, 22:51735177. PubMed Abstract  Publisher Full Text

Peng H, Dang W, Cao J, Chen Y, Wu D, Zheng W, Li H, Shen ZX, Liu Z: Topological insulator nanostructures for nearinfrared transparent flexible electrodes.
Nat Chem 2012, 4:281286. PubMed Abstract  Publisher Full Text

Dordevic SV, Wolf MS, Stojilovic N, Lei H, Petrovic C: Signatures of charge inhomogeneities in the infrared spectra of topological insulators Bi_{2}Se_{3}, Bi_{2}Te_{3} and Sb_{2}Te_{3}.
J Phys Condens Matter 2013, 25:075501. PubMed Abstract  Publisher Full Text

Hafiz MM, ElShazly O, Kinawy N: Reversible phase change in Bi_{x}Se_{100x} chalcogenide thin films for using as optical recording medium.
Appl Surf Sci 2001, 171:231241. Publisher Full Text

Zhao J, Liu H, Ehm L, Dong D, Chen Z, Gu G: Highpressure phase transitions, amorphization, and crystallization behaviors in Bi_{2}Se_{3}.
J Phys Condens Matter 2013, 25:125602. PubMed Abstract  Publisher Full Text

EM Explorer http://www.emexplorer.net/ webcite

Johnson PB, Christy RW: Optical constants of the noble metals.
Phys Rev B 1972, 6:43704379. Publisher Full Text

Berenger JP: Threedimensional perfectly matched layer for the absorption of electromagnetic waves.
J Comput Phys 1996, 127:363379. Publisher Full Text

Born M, Wolf E, Bhatia AB: Principles of Optics. Cambridge: Cambridge University Press; 1997:6170.

Nicolson AM, Ross GF: Measurement of the intrinsic properties of materials by timedomain techniques.

Smith DR, Schultz S, Markos P, Soukoulis CM: Determination of effective permittivity and permeability of metamaterials from reflection and transmission coefficients.

Chen XD, Grzegorczyk TM, Wu B, Pacheco JJ, Kong JA: Robust method to retrieve the constitutive effective parameters of metamaterials.

Zhang S, Fan W, Malloy KJ, Brueck SRJ: Nearinfrared double negative metamaterials.
Opt Express 2005, 13:49224930. PubMed Abstract  Publisher Full Text

Ortuño R, GarcíaMeca C, RodríguezFortuño FJ, Martí J, Martínez A: Role of surface plasmon polaritons on optical transmission through double layer metallic hole arrays.