Fundamental Role of Periodicity and Geometric Shape to Resonant Terahertz Transmission

Optical properties of the plasmonic structures based on the surface of variously structured metal plates became a subject of intense research ever since Ebbesen et al. reported the enhanced optical transmission through subwavelength hole arrays [1]. Their experimental results provided evidence that the enhanced optical properties originate from the coupling between the incident light and surface plasmon modes. Studies of such plasmonic structures and its mechanisms have a long history in the development of electromagnetism, though it has been an intense subject of research in recent years.


Introduction
Optical properties of the plasmonic structures based on the surface of variously structured metal plates became a subject of intense research ever since Ebbesen et al. reported the enhanced optical transmission through subwavelength hole arrays [1]. Their experimental results provided evidence that the enhanced optical properties originate from the coupling between the incident light and surface plasmon modes. Studies of such plasmonic structures and its mechanisms have a long history in the development of electromagnetism, though it has been an intense subject of research in recent years.
In the 1900s, Wood and Rayleigh reported the anomalous optical properties of reflection gratings [2,3], which are intimately connected with the excitation of surface bound waves on the metal-dielectric interfaces [4]. Hessel and Oliner presented two types of the anomalies occurring in the reflection grating surface: a Rayleigh wavelength type attributed to the usual propagating surface plasmon modes and a resonance type which is related to the guided waves supported by the grating itself [5]. On the other hand, since the 1970s, the optical properties of transmission gratings or grids have been occasionally studied in the microwave, terahertz (THz), and infrared regions [6][7][8][9][10]. In particular, Lochbihler and Depine in 1993 presented a theoretical approach for calculating the fields diffracted by the gratings made of highly conducting wires [11]. Here they calculated Maxwell's equations based on modal expansions inside the metal and Rayleigh expansions outside the metal, applying the theoretical treatments by means of surface impedance boundary condition method [11,12].
This research field reinvigorated after the Ebbesen's pioneering work which was focused on the understanding of mechanisms of the enhanced transmission through periodically perforated holes or slits on metallic plates. The one-dimensional plasmonic structures with periodic arrays of slits have been studied by some groups with theoretical and experimental approaches in the optical region [13][14][15][16][17][18]. These results strongly support the existence of two mechanisms of generating enhanced optical transmission: the geometric shape controlled mainly by the slit thickness and the periodicity determining the excitation of coupled surface plasmon polaritons (SPPs) on both surfaces of the perforated metallic structures. In particular the former called as a Fabry-Perot-like behavior was studied theoretically in a single narrow slit and clearly verified by strong resonant transmission of microwave radiation [19,20].
On the other hand in two-dimensional plasmonic structures with periodic arrays of holes, the geometric shape effect has been underestimated in spite of many researches about the strong transmission enhancement in the optical region [21][22][23][24][25][26]. Recently, the geometric shape resonance however was emerged as one of the dominant mechanisms of the strong transmission enhancement [27][28][29]. In particular, Koerkamp et al. already noted that the strong optical transmission is strongly influenced by a hole shape and a shape change from circular to rectangular [30]. The importance of the geometric shape was also confirmed by studies of strong enhancement of the light transmission through random arrays of holes [30][31][32] and single apertures [33][34][35][36], in particular a single rectangular hole theoretically studied by Garcia-Vidal et al. [37].
Unlike most typical systems in optical and infrared regions, THz time-domain spectroscopy system based on coherent THz radiation is ideally suited for phase-sensitive, broadband transmission measurements [38][39][40][41]. Furthermore, long wavelengths in THz region, 3000 μm to 30 μm (0.1 THz to 10 THz) in general, permit relatively simple fabrication of free-standing samples which can exclude complex surface modes by the substrate. For these reasons, the coherent THz waves have been generally used to study SPP behaviors [42][43][44][45][46][47][48][49][50][51][52][53][54][55]. The SPPs in the strictest sense however do not exist on perfectly conducting metal-dielectric interfaces since surface electromagnetic fields on the dielectric side cannot be strongly confined to the interface. However, recently the concept of the SPP has been extended even to include perfectly good metals with corrugations, which is referred to simply as designer (or spoof) SPP [56][57][58]. Pendry et al. reported the electromagnetic surface excitations localized near the surface which are governed by an effective permittivity of the same plasma form [56].
In the THz region, the most important thing is that the perfect transmission, of up to nearunity, at specific frequencies can be achieved since most metals become perfect conductor and therefore have extremely small ohmic loss fraction. Until now, a few groups have theoretically predicted the perfect transmission [59][60][61][62][63][64] and experimentally observed the near-unity transmission [65][66][67] in the plasmonic structures with periodic arrays of slits and holes. However, despite great progresses in realizing the perfect transmission, we do not have the clear understanding of what are the relative contributions of the geometric shape and the periodicity effects toward the perfect transmission in specific conditions.
Here, we first report on the importance of the combined effects of the geometric shape and periodicity in enhancing THz transmission through one-dimensional periodic arrays of slits. Theoretical predictions based on perfect conductor model show that the measured transmission peaks lie within the broad geometric shape band. We also report that the perfect transmission can be realized by the geometric shape resonance appearing in spectral region below first Rayleigh minimum. Calculated near-field distributions confirm that the perfect transmission is caused by the concentration of the electric and magnetic fields within the silts. Furthermore, we find that, despite the existence of only the geometric shape effect in random arrays of slits of relatively thick thicknesses, the enhanced transmission comes to the perfect transmission in the long wavelength region. Angle dependent transmission spectra and relative phase measurements clarify the relative roles of the geometric shape and periodicity toward the perfect transmission.

Experimental setup
In our experiments, we used a standard THz time-domain spectroscopy system with a spectral range from 0.1 to 2.5 THz, based on a femtosecond Ti:sapphire laser which generates optical pulses with average power of 600 mW and temporal width of about 120 fs set at a center wavelength of 760 nm. Experimental setup is shown in Fig. 1. A photoconducting antenna method is used to generate THz pulses. A crystal for the emitter is undoped, <110> orientated semi-insulating GaAs (INGCRYS Laser System Ltd.) with dark resistivity in range 5.3×10 7 to 5.6×10 7 Ohm cm. Two metallic electrodes are painted by silver paint on one side of the crystal, having a separation of about 0.4 mm which is among a large aperture emitter. Unlike a typical biased system of the photoconductive antenna, an ac bias voltage of 300 V and 50 kHz square wave was applied to the electrodes of the crystal [68]. To make this system, we used a high voltage pulse generator (DEI Model PVX-4150) which generates the squared high-voltage sources by combining DC high voltage source from a power supply (F. u. G. Model MCL 140-650) with square function generated using a pulse/function generator. Optical pulses generated from the femtosecond Ti:sapphire laser are divided into two parts, pump and probe beams, by a beamsplitter. The pump beam is focused on the surface of the emitter by using a 25 mm focal lens and immediately generates the THz wave which is essentially single cycle electromagnetic pulses. The generated THz wave is collimated by using an off-axis parabolic mirror at a focal distance. The degree of collimation is of the order of ∆θ ~ λ/D ~1/100 radian because our beam size D is about 5 cm and the THz wave is a point-like source. The collimated THz wave with the beam diameter of about 5 cm impinges on a reference box or samples. The angle-dependent transmission spectra are measured by tilting the sample stage in range -10 to 50 degrees. The THz pulses transmitted through the samples meet the optical probe beam again on an electro-optic crystal (ZnTe). The polarization change of the optical probe beam due to the Pockels effect is detected by a home-built differential photo-detector and provides us information of both amplitude and phase of the THz wave [69,70]. The signals are automatically acquired by a synchronized system (as shown in Fig. 1) of three components, a motion controller which tilts the sample stage, a data acquisition system, and a delay stage which controls time delay between pump and probe optical pulses. The whole setup is also enclosed in a box which is purged with dry nitrogen gas to reduce the absorption effects of water vapor.
We performed our experiments as follows. First, the collimated THz wave impinges on a reference box made of a 2 cm by 2 cm square hole punctured on a piece of aluminum. We measure the transmission spectrum through the reference box. We then mount our samples with the plasmonic structures of periodic or random arrays of slits perforated on aluminum plates, right on top of the reference box, and measure the transmission spectrum. Our system provides the time-resolved trace of the transmitted THz electric field with subpicosecond temporal resolution as an original signal. The fast Fourier transforms of the measured time-domain THz pulses give us information on both spectral amplitude and phase. The normalization is carried out over the entire frequency simply by dividing the spectral amplitude of the transmitted THz pulse with one of the reference beam transmitting through the reference box.

Sample fabrication
The samples were fabricated by a micro-drilling and a femtosecond laser machining methods. The latter is based on laser ablation which is performed by amplified femtosecond pulses. The pulses have energy of up to 1 mJ centered at 800 nm and a repetition rate of 1 kHz. The positions of the focused laser pulses are accurately controlled by using a galvanometer scanner (Scanlab AG, Germany) and a mechanical shutter system. The line edge roughness of the perforated slits is controlled to make successful slit structures, having the LER 3σ value of less than 3 μm which is two orders of magnitude less than the THz wavelengths of interest. The achieved high quality originates from the superiority of femtosecond laser machining system. The samples manufactured by this method have a fixed period d of 500 μm and different thicknesses h of 17, 50, 75 and 153 μm respectively. The silt widths a of these samples are 78, 80, 85, and 83 μm respectively, which are properly designed to be subwavelength [as shown in Fig. 2(a)]. On the other hand, the samples with thicknesses of 200 and 400 μm were fabricated by the micro-drilling method by using a screw with a diameter of 100 μm which is a minimum size in diameter of available and commercial screw. These thick plates were chosen for explicitly observing Fabry-Perot-like modes inside the slits. In particular the samples with thickness of 400 μm were considered suitable for randomization of the slit structures, based on the different frozen-phonon fashion of neighboring silts [as shown in Fig. 2(b)]. The position of the ith slit by the frozen-phonon fashion is given by where 0 x is the position of the basis of coordinate and i x Δ is the random amount of the ith slit which is getting by generating random values [72]. We displaced the position of each slit with random amounts of ± d/4 and ± d/2 which, in this chapter, are called as Random Type I and Random Type II.
Comparison between transmission properties through random and periodic arrays of slits reveal that the transmission spectra strongly affected by the geometric shape show similar spectral waveforms in both random and periodic arrays, since the geometric shape resonance depends on intrinsic properties of the structure factors not the structural arrangement of the slits. The geometric shape in the one-dimensional plasmonic structures is determined by the thickness h of the slits, providing the Fabry-Perot-like modes 2 where c is the speed of the light in vacuum and n is the integer mode index. In contrast, the periodicity effect might be found in transmission spectra through the periodic arrays not random, having different constructive interference patterns and different far-field transmissions. Therefore the relative contribution of the geometric shape and periodicity is directly compared in this chapter.

Theoretical framework
Herein, we are first interested in theoretical approach of THz electromagnetic wave transmission through the periodic arrays of slits, which may allow us to improve our intuitive understanding of the underlying physics and to find the mechanisms enhancing the transmission. To calculate the transmittance, the metallic structures are considered as a perfect conductor so that finite conductivity effects are neglected. This is a good approximation in the THz region, since the real and imaginary parts of dielectric constant of most metals are of the order of -30000 and 100000 respectively. We also consider a simple model system with rectangular shaped slits as shown in Fig. 2(a). In this system, the slits are along the y axis and a direction of propagation of the incident THz wave is the z axis which is normal to the sample surface. The polarization of the incident plane wave is placed along the x axis (TM polarization). Here, there are three different regions: the incident and reflection region (region I), the metal and perforated slit region (region II), and the transmission region (region III).
Let ( , ) f x z be the spatial function of the component along the z axis of the magnetic field. The fields can be expressed in terms of Rayleigh expansions in region I and III and modal expansion in region II [11,46]. The boundary matching is only carried out at the metaldielectric interfaces of I-II and II-III since which are semi-infinite. Here we use the single mode approximation inside the slits since the slit width is much smaller than the wavelength. The magnetic fields are described as ( 2 for the continuity of the normal derivative of the field at 2 for the continuity of the field at 2 z h = , and for the continuity of the normal derivative of the field at 2 z h = .
With the orthogonality, we apply some of integral equation to obtain simple solution with proper projection and integration. Multiplying a term where the quantities n P and n Q appearing in Eqs.
0 1 t a n ( 2) and 0 2 2 where W is defined as For normal incidence, zero-order transmittance for the periodic arrays of slits is expressed as follows

Enhanced transmission in thin metal plates
Theoretical calculations are first compared with experimental results of samples with the thicknesses of 17, 75 and 153 μm, which are carefully chosen so that the peak positions due to the geometric shape resonance are placed on the higher spectral region than the first Rayleigh minimum. Measured THz time traces and corresponding spectral amplitudes obtained by the fast Fourier transform method are shown in Fig. 3(a) and 3(b) respectively. The top curves in Fig. 3(a) and 3(b) show the reference signals transmitted through the reference box. The single cycle THz pulses with a temporal width of about 2 ps in time domain are transformed to the spectral pulses with a well-formed shape centered at 0.5 THz and a broad spectral width. In our experimental setup, we effectively removed multiple reflections of the THz pulses occurring from the emitter, the electro-optic crystal and a beamsplitter, which can help assure good quality of the measured time-resolved THz signals.
The thinnest sample with the thickness of 17 μm does not have the peak position of the enhanced transmission locating inside the presented spectral region. However, the two curves for the thicknesses of 75 and 153 μm in Fig. 3(a) have significant long-time oscillations with the periods of 0.87 and 1.80 picoseconds respectively, which also appears in the Fourier transformed spectra as peaks at 1.14 and 0.56 THz respectively (Fig. 3(b), gray arrows in bottom two curves). Here the peak positions of the enhanced transmission appear at the frequency regions slightly below the Rayleigh minima R f c d = , which indicate that the transmission properties are related to the periodicity effect. On the other hand, the peak positions shift toward longer wavelength region as the sample thickness increases, keeping the typical characteristics of Fabry-Perot-like resonance. These observations indicate that the enhanced transmissions are attributed to the combined effects of the geometric shape resonance and the periodicity.
Theoretical calculations shown in Fig. 3(c) are in well agreement with the experimental results except for extremely sharp transmission peaks immediately below the first Rayleigh wavelength. In particular in the cases of the two thinnest samples, the peak linewidths have almost infinitesimal values less than several gigahertz, which cannot be experimentally observed because of the finite temporal range of time domain signals and finite sample size. Theoretically predicted peak suddenly appears at 0.56 THz for h=153 μm as shown in bottom curve of Fig. 3(b). This appearance results from the approach of the geometric shape resonance toward the first Rayleigh minimum at which the first diffracted orders become evanescent. The independent contribution of Fabry-Perot-like resonance can draw out from single-slit approximation which can be achieved by carrying the distance between slits to be infinite.  At this point, we can expect that perfect transmission phenomenon will be appeared under the condition of overlapping between the Fabry-Perot-like resonance band and totally evanescent spectral regions below the first Rayleigh wavelength. The first diffracted waves touch down on the metal surface at the wavelength of the first Rayleigh minimum, at which wavelength eventually all higher diffracted orders become evanescent which is the reason why the total incident energy can be converged to the zero-order transmission amplitude. We therefore design a thicker sample to realize the perfect transmission.

Perfect transmission in thick metallic plates
The sample with the thickness of 200 μm eventually has broad Fabry-Perot-like resonance band centered at 0.45 THz which is located in the spectral region below the first Rayleigh minimum appearing at 0.6 THz as shown in Fig. 5(a). Angle dependent transmission amplitudes [ Fig. 5(b)] show the perfect transmission (in this sample, over 99% at normal incidence) centered at 0.5 THz. Despite the two transmission minima lines caused by Rayleigh frequencies for the ±1 diffraction orders, the perfect transmission appears at a broad incident angle range up to 10° because of the strong contribution of the Fabry-Perotlike resonance. Theoretical calculation with Eq. (20) is in good agreement with experimental result at normal incidence.
The perfect transmission phenomenon necessarily appears when the peak position of the enhanced transmission is located anywhere in the spectral region below the first Rayleigh frequencies, since all the diffracted orders of the transmittance Tn except for zeroth order become zero. A very useful theoretical condition for the perfect transmission can be therefore derived by setting the zeroth order reflection coefficient at zero as following Below the first Rayleigh frequencies, W can be considered as since all the reflected diffracted orders except for zeroth order become evanescent. Where γ is an imaginary part of kW . We can then write the condition of the perfect transmission as  Near-field distributions of THz waves above the metal surfaces can also give a key for understanding the perfect transmission phenomenon. The enhancement factor of the near electric fields onto the slits can be intuitively predicted by using the language of the surface impedance as where slit E , metal E and 0 H are the electric fields along the x direction at the slits and the metal and the incident magnetic field respectively. Since the effective surface impedance becomes one under the condition of the perfect transmission and there is no tangential component of the electric field on the metal surface, we can write Eq. (25) as The degree of the electric field enhancement is not only proportional to the incident magnetic field, but also inversely proportional to the areal sample coverage a d β = For the case of the areal sample coverage of 0.2 β = ( 100 a = and 500 d = μm), the theoretical results show the enhancement factor of 5 on the entrance of the silts as shown in Fig. 6(a). At a specific frequency at which the perfect transmission is realized, the incident waves are focused on the entrance of the slits due to the strong evanescent waves built up by a touchdown of all orders of diffracted modes on metal surface, followed by coupled with the Fabry-Perot-like resonant modes inside the slits, and then the focused waves reached at the exit of the slits are completely converted to the transmission.
These results suggest a possibility of realizing even stronger enhancement of the electric field strength onto the slits. We therefore fabricated a sample with the areal sample coverage of 0.1 β = by increasing the period. (Reducing the slit width is actually more effective, but the technical limit of the micro-drilling method to fabricate thick metallic plates is 100 μm.) Calculated near-field distribution of the sample is given by Fig. 6(b), showing the enhancement factor of 10 corresponding to the inverse of the areal sample coverage. The near-field enhancement necessarily induces the transmission enhancement into the farfield. Measured angle dependent transmission spectra shown in Fig. 7 show not only angle-independent, broad spectral bands near 0.3 and 0.6 THz caused by the geometric shape resonance, but also a strong transmission peak with the value of about 0.85 at 0.26 THz at normal incidence. Though the peak value does not reach the THz transparency because of a narrow spectral bandwidth, the enhancement factor compared with the areal sample coverage of 0.1 is over eight. Controlling the areal sample coverage by using the period or width is therefore extremely important in realizing a high power near-field THz source.

THz transparency by geometric shape
At this point, it is important to understand that what the relative contributions of the periodicity and the geometric shape resonance toward the perfect transmission are in a given situation. From the experimental viewpoint, the best method to compare the difference between the relative contributions of two mechanisms is to fabricate the samples with random arrays and to investigate its transmission properties. The calculated transmission spectra as shown in Fig. 8 play an important role in determining the appropriate sample parameters. For the convenience of comparison, the transmission spectra are plotted versus the first Fabry-Perot frequency, fc=c/2h instead of the slit thickness.
In the cases of the relatively thinner metal plates, there are two peak maxima of the transmission spectra: one corresponds to the first Fabry-Perot-like mode, fc=c/2h, and the other corresponds to the so-called zeroth Fabry-Perot-like mode appearing at the long wavelength region near 0 THz, which is originated from the geometric shape effect of the semi-infinite slits in a thin metal film. In the long wavelength region near 0 THz, an aluminum plate with the thickness of several hundreds of microns can be actually regarded as an optically thin metal film. With increasing the thickness of the aluminum plates, the peak maxima originating from higher Fabry-Perot-like modes such as fc=2c/2h, 3c/2h, and etc are added consecutively, maintaining the perfect transmission.
All curves of the Fabry-Perot-like modes show that the curves approach to straight lines of the classical Fabry-Perot positions (Fig. 8, dotted lines) as getting away from the wavelength of the first Rayleigh minimum of 0.6 THz. The calculated results may therefore present that, when the peak maxima approach to the long wavelength limit, the contribution by the coupled evanescent surface modes presented in Eq. (19) becomes weak. In our experiments, we therefore fabricated the samples of periodic and random arrays of slits with the thickness of 400 μm since in which there are three different peak maxima originating from the zeroth ( Fig. 8, C), first (Fig. 8, B) and second (Fig. 8, A) Fabry-Perot-like modes and each peak maximum may be stimulated by different contributions of the periodicity and the geometric shape. In fact, as the peak position is getting away from the wavelength of the first Rayleigh minimum toward longer wavelength region, the geometric shape effect becomes dominant.   Figure 9 shows the normalized transmission amplitudes of periodic and random arrays of silts at normal incidence. Here the samples have same slit width a of about 100 μm and same areal sample coverage of β=0.2. As predicted in theoretical calculations, the calculated transmission spectrum (Fig. 9, a solid line) shows three peak maxima at 0.56, 0.30 and near zero THz, corresponding to the second, first and zeroth Fabry-Perot-like modes respectively. The measured transmission spectrum (Fig. 9, a dark gray line) of the periodic arrays of slits agrees well with the calculated one, except slight disagreement at the peak maximum of the second Fabry-Perot-like mode. Here the experimental linewidth of the second Fabry-Perotlike mode does not completely reproduce the theoretical one since which is too sharp to experimentally observe in real sample.
We then consider two random structures with the different frozen-phonon fashion of neighboring slits. The Random Type I and II samples were fabricated by displacing the position of each slit of periodic arrays with random amounts of ± d/4 and ± d/2 respectively. The transmission spectrum of the Random Type I (Fig. 9, a bright gray line) shows a slightly lower peak value at 0.30 THz and a considerably decreased peak value at 0.56 THz, but overall the spectral shape of the transmission spectra are well matched, particularly reproducing the peak value due to the zeroth Fabry-Perot-like mode. The measured transmission spectrum of the Random Type II (Fig. 9, a dotted line) clearly exhibits the spectral properties of the perfectly random arrays of slits, which does not show any more the transmission minimum at 0.6 THz originating from the first Rayleigh minimum due to the periodicity effect only. On the other hand, although the peak value appearing near 0.3 THz is slightly decreased with further increasing the fraction of random amounts, the spectral shape and peak position can be not only considered unchangeable, the resonant peak also exhibits an enhancement factor larger than 3.5. Furthermore, the resonant peaks of the zeroth Fabry-Perot-like mode are more dramatic where all the transmission spectra in theory and experiments are perfectly matched, showing the transmittance more than 80%. The relative phase of transmitted THz waveforms not only improves the understanding for the relative contributions of the geometric shape and the periodicity but also makes more clearly what the dominant mechanism in enhancing the THz transmission is. When the geometric shape effect becomes dominant, there will be a little difference between the random and periodic arrays of slits, as shown in Fig. 10. In addition, the geometric shape does not depend on the incident angle of the THz waves. Shown in Fig. 11 are the angle dependent transmission amplitudes for the samples. For the sample of periodic arrays of slits, there are two dominant transmission minima lines crossing at 0.6 THz at zero angle which is the Rayleigh wavelengths corresponding to the ±1 diffracted orders as The Rayleigh wavelengths determined by the grating period distinctly appear as transmission minima lines over all the incident angle of the THz waves, generating an angle-dependent and strongly enhanced transmission peak near 5.4 THz. With increasing the randomness, the transmission minima lines of the Rayleigh wavelengths become indistinct as shown in Fig. 11(b) and eventually disappear as shown in Fig. 11(c). Similarly, the angle dependent transmission peak appearing at the sample of periodic arrays of slits also disappears since its strong enhancement is caused by the strong field accumulation of the evanescent surface waves by the periodic structures. On the contrary, the enhanced peaks appearing near 0.3 THz are essentially angle independent, making the strong enhancement. Figure 11. Normalized angle dependent transmission amplitudes for three types of the samples: periodic arrays (a), random type I (b), and random type II (c).

Conclusion
In conclusion, we have demonstrated two mechanisms for perfect transmission appearing at specific frequencies through the periodic and random arrays of slits. The theoretical results show that the perfect transmission can be realized in spectral region below the first Rayleigh minimum determined by the periodicity. Under the condition for the perfect transmission, the symmetric electric-and magnetic-field profiles inside the slits are excited and the incident energy is funneled onto the slits with the enhancement factor proportional to the inverse of the areal sample coverage. The measured transmission spectra through the slits of thick metal plates show the enhanced transmission of perfect transmission, which is attributed by the increase of the spectral linewidth resulting from the approach of the geometric shape resonance to the spectral region below the first Rayleigh minimum. The resonant bands of the geometric shape presented by the single silt approximation and the measured angle dependent transmission spectra not only are angle independent but also determine the spectral peak positions of the perfect transmission.
The random arrays of slits designed to remove the effect of the periodicity clearly show the dominant contribution of the geometric shape resonance toward the enhanced transmission.
In the long-wavelength region far from the first Rayleigh minimum, the geometric shape resonance becomes dominant and the transmission enhancement factors are almost equal to the cases of the periodic arrays of slits. The strong transmission and local field enhancements open a possibility of potential applications in perfect transmission of solid materials such as semiconductors and biological materials requiring the high power THz sources and in designing photonic devices such as transmission filters, low-pass frequency filters, and transmissive waveguides requiring high transmission at desired frequencies.