Extraordinary Optical Transmission through Single Sub-Wavelength Slot Nano Antennas

This paper investigates the electromagnetic transmission through the sub-wavelength slot model on a metal film for TM- and TE-polarized light. The influence of several parameters such as the slot width, the metal film thickness and the polarization of the incident field is investigated using FDTD method. 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 EOT has been observed as being symmetrically punctured 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. 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, are such examples.

tonics do not have, i.e. the size of electronics and the speed of photonics [1]. 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 [2] [3], 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 [4]; a simple thin film of metal can serve as an optical lens [5]; and sliver films with nano-scale holes show extraordinarily optical transmission (EOT) [6] 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 [7] [8] [9] [10] [11]. 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 [12] [13] [14] and information and communications technologies [15] and sensing with high spectral and spatial precision [16] [17], coupling devices [18], enhanced solar cells [19] [20] [21], efficient optical sources [22] [23] [24] and ultrafast acousto-magneto-plasmonics [25] detectors [26], disease treatment [27] [28], nano laser [29], 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) [30], Finite-Difference Time-Domain (FDTD) [31], Beam propagation method (BPM) [32], Method of Moment (MOM) [33] etc. Due to the known limitations in accuracy achievable using BPM and MOM [32] [33], 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. Optics and Photonics Journal 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.

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 ε . (Au) is the second most significant metal after silver for plasmonic applications [34]. 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 [35].
where ω is the angular frequency of the incident light, є ∞ the dielectric constant at an infinite angular frequency and equals 3.7, p ω the bulk plasma frequency and equals 16 1.3709 10 rad s * , and γ the electron collision frequency equals 13 4.059 10 rad s * [36]. 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), [37] [38] [39] 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.

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 finite-difference-time-domain (FDTD) method [40] 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 sufficient for convergence of the numerical results. Calculations in FDTD are performed according to the steps listed in Figure 2 [41] [42].

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.

TE-Polarization
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 (H x , H y and E z ) 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 |E z | field distribution in Figure 3), while H x at this region is a maximum. This is due to the fact that because Maxwell's equations require E z 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 J z surface current. In contrast, continuity at the interface is not required for H x because of the existence of the surface current J z [43].
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 E z creates a surface current J z along the slot walls, which, in turn produce an H y on the walls' exterior surfaces as seen in |H y | 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 p y (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

TM-Polarization
The   being supported by the appropriates distribution of surface charges on these walls. Figure 6 shows simulated field allocation of |E x |, |E y | and |H z | 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 J x 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 |E y | 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 [45]. 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 E y distribution field displayed in Figure 6.  For the sub-wavelength slot with thickness of t = 500 nm, only one strong dipole is observed at the bottom of the slot (see |E y | 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 [43].  nm strong electric dipoles once again at the sharp edges of the slot (see |E y | and |E x | 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 |H z | 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 |E x | and |E y | 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 |H z | field distribution in Figure 8.
In Figure 9 we show the simulated field distribution of the simulated pointing vector p y (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   In Figure 11 we summarize the study of the energy flux (i.   waves at wavelength 800 nm will not transmit, while it will transmit for width 500 nm.

Conclusion
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