Unluckily, the fundamental laws of diffraction limit the size of photonics, while the interconnect RC delay-time issues limit the speed of semiconductor electronics. Over the last two decades, it has gradually become clear that the plasmonic devices can play in the future technologies to complement conventional photonics and electronics. Plasmonics offers exactly what electronics and photonics do not have, i.e. the size of electronics and the speed of photonics . In addition, plasmonic devices might logically interface with comparative-speed photonic devices and with comparative-size electronic components, so, expanding the collaboration between these technologies. Therefore it can serve as a bridge between similar-size nano-electronics and similar-speed dielectric photonics.
Plasmonics technologies have grown rapidly over the last two decades due to the exciting physical properties of Surface Plasmon Polaritons (SPPs). SPPs are quantized charge density oscillations occurring at the interface between the free electron gas of a metal and a dielectric material. The nature of these quasi-particles is that of an electromagnetic wave trapped at the surface between such two media. Motions of the charge density drive the electromagnetic wave, which can propagate for relatively large distances across the surface of the interface (in the order of micrometers, or even millimeters relying upon the materials and the frequency of operation). Even though surface-enhanced Raman spectroscopy (SERS), being the first applications of surface Plasmon metallic nanostructures, was discovered in the 1970s   , the field of plasmonics started to rapidly spread in the early 2000s. Around then it was exhibited that metallic nanowires can guide light well below the diffraction limit  ; a simple thin film of metal can serve as an optical lens  ; and sliver films with nano-scale holes show extraordinarily optical transmission (EOT)  The most essential applications in the plasmonics field appear to depend on one key property of engineered metallic structures that is the ability to guide and concentrate light at the nano-scale     . It is envisaged that plasmonics, combined with nano photonics, has the potential to provide ultra-small optoelectronic components having higher speeds and greater bandwidths. Investigations involve building, manipulating, and characterizing optically active nanostructures to innovate new capabilities for the nano-scale high-spatial-resolution imaging    and information and communications technologies  and sensing with high spectral and spatial precision   , coupling devices  , enhanced solar cells    , efficient optical sources    and ultrafast acousto-magneto-plasmonics  detectors  , disease treatment   , nano laser  , and many other applications. In spite of the fact that the basics for SPPs are well established, SPPs based structures are complex. Therefore, it is frequently not possible to achieve analytical solutions that characterize the operation of a device; in this way the numerical modeling becomes an essential tool to analyze and investigate the operation of such structures. Numerous numerical techniques have been developed, such as Finite Element Method (FEM)  , Finite-Difference Time-Domain (FDTD)  , Beam propagation method (BPM)  , Method of Moment (MOM)  etc. Due to the known limitations in accuracy achievable using BPM and MOM   , when modeling discrete structures such as slots with sharp structural transitions, it is not used heavily. In this paper, the numerical models have been primarily based on FDTD. In this paper, the light transmission through a single sub-wavelength slot is discussed. The influence of several parameters such as the slot width, the metal film thickness and the polarization of the incident field are investigated using FDTD method. The simulation results show that a sub-wavelength slot in the gold metal films has extraordinary optical transmission properties (EOT). The EOT is observed as being symmetrically perforated and characterized by the appearance of a series of transmission peaks and dips in the transmission spectrum This paper is organized as follows Section 2 introduces the proposed model and its dimension and materials parameters Section 3 describes the FDTD method where the parameters needed to simulate this model is presented. The results are presented in Section 4. Finally, Section 5 summarizes our concluded remark.
2. The Single Sub-Wavelength Slot Model
In the next sections we will study the conditions of propagation inside an infinitely long metallic slot. A schematic of the system is represented in Figure 1(a) thin dielectric layer (air) of width w and thickness t and depth L permittivity is sandwiched between two semi-infinite metal areas (goldAu) of permittivity . (Au) is the second most significant metal after silver for plasmonic applications . It is widely chosen for many applications over silver because of its stability. Moreover, Au is an easily accessible material that can be evaporated to a thin smooth layer on a surface. Because of these two reasons, throughout this paper, we will concentrate our focus on gold. Au dispersive permittivity can be denoted by the Drude model .
where is the angular frequency of the incident light, the dielectric constant at an inﬁnite angular frequency and equals 3.7, the bulk plasma frequency and equals , and the electron collision frequency equals . A simplified picture of the 2D slot in a metal layer Figure 1 is that of a waveguide, where the nano-slot forms a vertical metal insulator-metal (MIM),    waveguide of finite length with thickness t. The 2D Finite Difference Time Domain (FDTD) method allows us to directly use experimental data for the frequency-dependent dielectric constant of metals such as gold including both the real and imaginary parts, with no further approximation.
3. FDTD Model
In this section, we investigate the nano slot properties of the nano slot waveguide in a two dimensional (2D) plasmonic waveguide with Au-air-Au interfaces using a ﬁnite-diﬀerence-time-domain (FDTD) method  and perfectly matched-layer (PML) as an absorbing boundary condition (ABC). In our numerical calculations with FDTD, the nano slot waveguide are described by a spatial discretization mesh. A spatial mesh size of Δx = Δy = 9 nm in the FDTD algorithm is found suﬃcient for convergence of the numerical results. Calculations in FDTD are performed according to the steps listed in Figure 2  .
Figure 1. Schematic of the single sub-wavelength slot structure and parameters’ definitions used for the structure dimensions and incident light field directions.
Figure 2. Layout of FDTD method .
This system is made up of a 2D box (along the x & y-axes) with propagation along the y-axis. UPML boundary conditions are applied at both boundaries in y and x of the box in order to eliminate the reﬂection of outgoing waves . 50 nm width for the PML is found to be thick enough to fully absorb outgoing waves. The structure is supposed to be inﬁnite along the z-direction. The equations of motion are solved with a time integration step seconds to avoid the problems of numerical dispersion and stability. The number of time steps is equal to 2, which is the necessary tested time for good convergence of the numerical calculations.
The wavelength of the incident light, λo = 830 nm, which lies in the near-infrared range.
4. Numerical Results
After the presentation of our analysis for the propagation characteristics of the electromagnetic waves throw the nano slots. The structure presented in Figure 1 is studied in two different cases; the TE and TM. Where, the direction of propagation is the y direction and the wavelength of the incident light equals λo = 830 nm.
The structure display in Figure 1 is simulated using FDTD technique under the TE polarized incident wave. Figure 2 shows the simulated EM field distributions (Hx, Hy and Ez) for a sub-wavelength slot with W = 350 nm (i.e. <1/2λo) and 300 nm film thickness. The weak EM fields down to the hole demonstrate that no light transmits through the sub-wavelength slot.
The electric field at the top metal surface is approximately zero (see the |Ez| field distribution in Figure 3), while Hx at this region is a maximum. This is due to the fact that because Maxwell’s equations require Ez to be continuous at the interface, and negligible inside the metal. As seen in Figure 2, very little electric field is needed to oscillate the surface conduction electrons, and thus to establish the Jz surface current. In contrast, continuity at the interface is not required for Hx because of the existence of the surface current Jz .
As seen in Figure 3, the incident electric and magnetic fields carry the EM energy in through the sub-wavelength aperture, while the metallic film sends a strong reflected beam back. For a metal film with no aperture, these counter-propagating beams form a standing wave above the upper surface. The penetrating Ez creates a surface current Jz along the slot walls, which, in turn produce an Hy on the walls’ exterior surfaces as seen in |Hy| field distribution.
When W < 0.5 λo, very little light transmits through the sub-wavelength slot and the incident EM wave is reflected back. These simulations for TE-polarized incident light revealed the existence of a cut-off at around W ≈ 0.5λo for TE-polarized waves through the sub-wavelength slot. In contrast, as shown in Figure 4, when W = 500 nm (i.e., > ½λo), there is little attenuation and the EM waves transmit through the slot.
In Figure 5, we show the simulated field distribution of the pointing vector py (i.e., energy flowing along the y direction) under TE-polarized incident light for two values of the sub wavelength slot width; W = 350 nm (i.e., below cut-off) and W = 500 nm. Decreasing the slot-width below the slot cut-off reduces the transmission as shown in Figure 5(a)  .
The simulations described in this section show the demonstration of a cut-off for TM-polarized incident light with extraordinary light transmission, even for an aperture as narrow as W = 0.01λo. Transmission is seen to be the result of strong electric and magnetic fields that propagate along the sub-wavelength slot walls,
Figure 3. Simulated field distribution of |Ez|, |Hx| and |Hy| with TE-polarized incident light for film thickness t = 300 nm and slot width W = 350 nm.
Figure 4. Simulated field distribution of |Hx|, |Hy| and |Ez| with TE-polarization incident light for film thickness t = 300 nm and slot width W = 500 nm.
being supported by the appropriates distribution of surface charges on these walls. Figure 6 shows simulated field allocation of |Ex|, |Ey| and |Hz| for a 100 nm slot width in a 300 nm Au film thickness.
As for TE illumination, very little Ex is needed on the top metal surface to sustain the Jx surface current which supports the magnetic field Hz immediately above the surface. The reflected Ex and Hz interfere with the incident fields to produce standing waves above the top surface. The surface current stops at the edges of the slot, giving rise to accumulated charges at the slot corners (see |Ey| field distribution in Figure 6). These oscillating charges, on opposite edges of the slot, act as an electric dipole. These surface charges play a role in enhancing transmission through the slot . Inside the slot, the surface charges and currents carry the travelling beam along the slot; it creates a second electric dipole as seen in the Ey distribution field displayed in Figure 6.
Figure 5. Simulated field distribution of |py| with TE-polarized incident light for a film thickness t = 300 nm and slot width (a) W = 350 nm and (b) W = 500 nm.
Figure 6. Simulation field distribution of |Hz|, |Ex| and |Ey| for a film thickness of t = 300 nm and aperture width of W = 100 nm with TM-polarized incident light.
For the sub-wavelength slot with thickness of t = 500 nm, only one strong dipole is observed at the bottom of the slot (see |Ey| distribution field in Figure 7).
The charges that produce the dipole at the top of the slot have diminished, and the transmission efficiency is substantially reduced compared with the 300 nm film thickness. It appears that the distractive interference between two counter-propagating beams within the slot is responsible for the reduced transmission efficiency in this case . The metal film thickness of t = 300 nm and 500 nm are not the only interesting cases with large transmission. At t = 600
Figure 7. Simulated field distribution of |Hz|, |Ex| and |Ey| for a film thickness of t = 500 nm and aperture width of W = 100 nm with TM-polarized incident light.
nm strong electric dipoles once again at the sharp edges of the slot (see |Ey| and |Ex| fields allocation in Figure 8), and the great transmission is subsequently observed. The difference with the case of t = 600 nm is that, in the middle of the metal film the charges disappear while the wall surface current in that region reaches a maximum (see the |Hz| field distribution in Figure 8). It is evident that a strong current along the slit walls creates the top and bottom dipoles.
As shown in the |Ex| and |Ey| fields distribution in (Figure 8). The charges are accumulated in three places on each side of the slot walls: in the middle of the film and at the top and bottom corners of the slot. Also, the current on the slot walls is apparently stronger than that on the top surface as seen in the |Hz| field distribution in Figure 8.
In Figure 9 we show the simulated field distribution of the simulated pointing vector py (i.e., energy flowing along the −y direction) under TM-polarized incident light for the three values of the metal film thickness; t = 300 nm, 500 nm and 600 nm for 100 nm sub-wavelength slot width.
In Figure 10 we summarize the study of the energy flux (i.e., power) transmitted for an 830 nm TM-polarized incident light for different sub-wavelength slot thicknesses. The results are shown for a 100 nm wide slot while the Au film thickness varies from 50 nm to 600 nm with 20 nm increments.
The normalized transmitted power throughput from the sub-wavelength slot oscillates periodically as the Au thickness increases shows series of Fabry-Pérot like resonance. The transmission is maximal due to constructive interference when thickness t gives a length phase of even integer of π/2 and is minimal when thickness t corresponds to a length phase of odd integer of π/2. As seen from Figure 10, the maximum normalized transmission power is always greater than the unity (i.e. extraordinary transmission). These results indicate that some of the incident photons outside of the sub-wavelength slot area are either captured and “funneled” through the slot to the transmission region or the transmission is
Figure 8. Simulated field distribution of |Hz|, |Ex| and |Ey| for a film thickness of t = 600 nm and aperture width of W = 100 nm with TM-polarized incident light.
Figure 9. Simulated field distribution of |py| with TM-polarized incident light for a slot width W = 100 nm and film thickness of (a) t = 300 nm, (b) t = 500 nm, (c) t = 600 nm.
enhanced through the SPP waves which contribute to the final transmission power.
In Figure 11 we summarize the study of the energy flux (i.e. power) transmitted for wavelengths from a 500 nm to 1500 nm TM-polarized incident light. The results are shown for different slot widths; 300, 400, and 500 nm. It is well known from the electromagnetics theory that the cut-off frequency of the waveguide depends mainly on the width of the waveguide. It is clear from this figure that for each structure there is a resonance wavelength. This resonance wavelength increases as the width of the slot increases. In other words, as the width of the slot is reduced the cut off frequency increases while the cut off wavelength is decreased. In these interpretations, we considered the slot as a waveguide. Each waveguide has a dominant cutoff wavelength. For example at width equals 300 nm
Figure 10. Normalized transmission power versus sub-wavelength slit thickness with 100 nm width.
Figure 11. Normalized transmission power versus sub-wavelength for different widthes.
waves at wavelength 800 nm will not transmit, while it will transmit for width 500 nm.
The transmission through the sub-wavelength slot model on a metal film, namely gold, for TM- and TE-polarized light has been investigated in this paper. The FDTD simulation’s results have shown that the sub-wavelength slot in the metal films has extraordinary optical transmission (EOT) properties for TM-polarized light. The effect of different parameters has been studied such as width, thickness, and polarization of the incident wave. In these interpretations, we considered the slot as a waveguide. Each waveguide has a dominant cutoff wavelength. The EOT has been observed as being symmetrically perforated and characterized by the appearance of a series of transmission peaks and dips in the transmission field. By varying the slot thickness we have investigated the effect of the Fabry-Pérot like resonance in the sub wavelength slot. This component can be a key element in many applications such as high-spatial-resolution imaging and information and communications technologies and sensing with high spectral and spatial precision, enhanced solar cells, efficient optical sources and detectors, disease treatment which are such examples.
 Jeanmaire, D.L. and Van Duyne, R.P. (1977) Surface Raman Spectroelectrochemistry: Part I. Heterocyclic, Aromatic, and Aliphatic Amines Adsorbed on the Anodized Silver Electrode. Journal of Electroanalytical Chemistry and Interfacial Electrochemistry, 84, 1-20.
 Fleischmann, M., Hendra, P.J. and McQuillan, A.J. (1974) Raman Spectra of Pyridine Adsorbed at a Silver Electrode. Chemical Physics Letters, 26, 163-166.
 Takahara, J., Yamagishi, S., Taki, H., Morimoto, A. and Kobayashi, T. (1997) Guiding of a One-Dimensional Optical Beam with Nanometer Diameter. Optics Letters, 22, 475-477.
 Schuller, J.A., Barnard, E.S., Cai, W., Jun, Y.C., White, J.S. and Brongersma, M.L. (2010) Plasmonics for Extreme Light Concentration and Manipulation. Nature Materials, 9, 193.
 Krajnik, B., Czechowski, N., Piatkowski, D., Mackowski, S., Hofmann, E., Pichler, S. and Heiss, W. (2013) Influence of Plasmon Excitations in Au Nanoparticles upon Fluorescence and Photostability of Photosynthetic Complexes. Optics and Photonics Journal, 3, 1.
 Yao, J., Stewart, M.E., Maria, J., Lee, T.W., Gray, S.K., Rogers, J.A. and Nuzzo, R.G. (2008) Seeing Molecules by Eye: Surface Plasmon Resonance Imaging at Visible Wavelengths with High Spatial Resolution and Submonolayer Sensitivity. Angewandte Chemie International Edition, 47, 5013-5017.
 Shan, X., Diez-Pérez, I., Wang, L., Wiktor, P., Gu, Y., Zhang, L. and Li, J. (2012) Imaging the Electrocatalytic Activity of Single Nanoparticles. Nature Nanotechnology, 7, 668.
 Jung, K.Y., Teixeira, F.L. and Reano, R.M. (2007) Au/SiO2 Nanoring Plasmon Waveguides at Optical Communication Band. Journal of Light Wave Technology, 25, 2757-2765.
 Pala, R.A., White, J., Barnard, E., Liu, J. and Brongersma, M.L. (2009) Design of Plasmonic Thin-Film Solar Cells with Broadband Absorption Enhancements. Advanced Materials, 21, 3504-3509.
 Omelyanovich, M., Makarov, S., Milichko, V. and Simovski, C. (2016) Enhancement of Perovskite Solar Cells by Plasmonic Nanoparticles. Materials Sciences and Applications, 7, 836-847.
 Yu, N., Fan, J., Wang, Q.J., Pflügl, C., Diehl, L., Edamura, T., Capasso, F., et al. (2008) Small-Divergence Semiconductor Lasers by Plasmonic Collimation. Nature Photonics, 2, 564-570.
 Oulton, R.F., Sorger, V.J., Zentgraf, T., Ma, R.M., Gladden, C., Dai, L. and Zhang, X. (2009) Plasmon Lasers at Deep Subwavelength Scale. Nature, 461, 629-632.
 Gao, N., Huang, K., Li, J., Li, S., Yang, X. and Kang, J. (2012) Surface-Plasmon-Enhanced Deep-UV Light Emitting Diodes Based on AlGaN Multi-Quantum Wells. Scientific Reports, 2, Article No. 816.
 Ren, F.F., Ang, K.W., Ye, J., Yu, M., Lo, G.Q. and Kwong, D.L. (2011) Split Bull’s Eye Shaped Aluminum Antenna for Plasmon-Enhanced Nanometer Scale Germanium Photodetector. Nano Letters, 11, 1289-1293.
 Hirsch, L.R., Stafford, R.J., Bankson, J.A., Sershen, S.R., Rivera, B., Price, R.E., West, J.L., et al. (2003) Nanoshell-Mediated Near-Infrared Thermal Therapy of Tumors under Magnetic Resonance Guidance. Proceedings of the National Academy of Sciences, 100, 13549-13554.
 Li, Z., Wang, Y., He, J., Feng, D., Gu, E. and Li, W. (2016) An Ultraviolet Hybrid Plasmonic Waveguide for Nanolaser Applications. Optics and Photonics Journal, 6, 19-23.
 Yee, K. (1966) Numerical Solution of Initial Boundary Value Problems Involving Maxwell’s Equations in Isotropic Media. IEEE Transactions on Antennas and Propagation, 14, 302-307.
 Van Roey, J., Van der Donk, J. and Lagasse, P.E. (1981) Beam-Propagation Method: Analysis and Assessment. Journal of the Optical Society of America, 71, 803-810.
 Noual, A., Akjouj, A., Pennec, Y., Gillet, J.N. and Djafari-Rouhani, B. (2009) Modeling of Two-Dimensional Nanoscale Y-Bent Plasmonic Waveguides with Cavities for Demultiplexing of the Telecommunication Wavelengths. New Journal of Physics, 11, Article ID: 103020.
 Ordal, M.A., Bell, R.J., Alexander, R.W., Long, L.L. and Querry, M.R. (1985) Optical Properties of Fourteen Metals in the Infrared and Far Infrared: Al, Co, Cu, Au, Fe, Pb, Mo, Ni, Pd, Pt, Ag, Ti, V, and W. Applied Optics, 24, 4493-4499.
 Keleshtery, M.H., Kaatuzian, H. and Mir, A. (2016) Analysis and Investigation of Slow Light Based on Plasmonic Induced Transparency in Metal-Dielectric-Metal Ring Resonator in a Waveguide System with Different Geometrical Designs. Optics and Photonics Journal, 6, 177-184.
 Ghosh, S., Brahmachari, K. and Ray, M. (2012) Experimental Investigation of Surface Plasmon Resonance Using Tapered Cylindrical Light Guides with Metal-Dielectric Interface. Journal of Sensor Technology, 2, 48-54.
 Taheri, A.N. and Kaatuzian, H. (2013) Simulation and Design of a Submicron Ultrafast Plasmonic Switch Based on Nonlinear Doped Silicon MIM Waveguide. Journal of Computer and Communications, 1, 23-26.
 Mashade, E., Bakry, M. and Abdel Aleem, M.N. (2009) Analysis of Ultra-Short Pulse Propagation in Nonlinear Optical Fiber. Progress in Electromagnetics Research, 12, 219-241.
 Nady Abdul Aleem, M. (2017) Mode Splitting Based on the Coupling between Modes of Two Nanodisks Cavities and a Plasmonic Waveguide. Progress in Electromagnetics Research, 55, 179-188.
 Xie, Y., Zakharian, A.R., Moloney, J.V. and Mansuripur, M. (2004) Transmission of Light through Slit Apertures in Metallic Films. Optics Express, 12, 6106-6121.
 Schouten, H.F., Visser, T.D., Lenstra, D. and Blok, H. (2003) Light Transmission through a Subwavelength Slit: Waveguiding and Optical Vortices. Physical Review E, 67, Article ID: 036608.
 Mansuripur, M., Xie, Y., Zakharian, A.R. and Moloney, J.V. (2005) Transmission of Light through Slit Apertures in Metallic Films. IEEE Transactions on Magnetics, 41, 1012-1015.