Abstract

The paper presents experimental and numerical investigations on the terahertz (THz) surface plasmon propagation in a periodically patterned doped silicon substrate. Silicon substrates are periodically patterned with 2D array of vertical structures forming a plasmonic waveguide. The waveguide configurations are found to support resonant surface modes at certain frequencies which can occur anywhere depending on the structural parameters. The 2D pattern of vertical structures is observed to affect the THz surface plasmon propagation along the waveguide configuration. The periodicities are changed in both directions to examine the change in amplitude and cut-off frequencies of the resonant surface modes. The results are confirmed independently through finite element method based numerical simulations and compared with theory. The present study should find applications in facilitating the development of futuristic ultrahigh speed networks and devices operating at the terahertz frequencies.

1. Introduction

Recent years have witnessed significant developments in the field of plasmonic metamaterials. Metamaterials are the artificially designed subwavelength scale structures that have the ability to alter electromagnetic radiation propagation in a controlled fashion [17]. The artificially designed structures have been shown to guide and manipulate surface electromagnetic waves whose properties can be controlled with structural parameters [8, 9]. The guided wave devices utilizing plasmonic metamaterials have led to several applications which include nanoimaging [10], nanoscale photonic components [11, 12], slow light systems [13], and biosensing [1416]. There has been immense interest in developing guided wave components which can operate at terahertz frequencies. Terahertz frequency regime of the electromagnetic spectrum has the potential to significantly improve the efficiency and speed of the devices [17, 18]. In this direction, plasmonic metamaterials have been investigated to play an important role, and guided wave components and devices operating at terahertz frequencies have been developed [1921]. The metals have been favorable choice for this purpose as they exhibit negligible dielectric and ohmic losses [22]. Williams et al. have demonstrated that a two-dimensional array of groove pattern with metal coating can support highly confined terahertz surface modes [23]. The metal coating is considerably thick (more than skin depth) to reduce the propagation losses. Zhu et al. also later on showed that when a thin metal sheet of stainless steel is periodically patterned with the one-dimensional array of apertures, then it supports high confined terahertz surface modes at certain frequencies which can be defined by the geometrical parameters of the apertures [24]. Subsequently, a significant amount of work has been noticed in this direction with different shapes of the metallic corrugations [2529]. It is important to mention that these experimental investigations on guiding surface electromagnetic waves are inspired by Pendry’s seminal work on plasmonic metamaterials [30, 31].

In spite of an extensive research work in the field of plasmonic metamaterials using metals, we have not been able to achieve active waveguide components such as modulators, ultrafast switches, filters, and active resonators. This is because of the difficulty in tuning the properties of the metal based waveguides. In order to overcome this issue, recently Kumar et al. demonstrate the planar plasmonic terahertz waveguides based on heavily doped silicon [32]. These waveguides can be more readily fabricated using photolithographic techniques in a clean room and properties of the substrates can be altered by changing the carrier concentration when externally stimulated, that is, by means of a laser or photo illumination. It has been found that a doped silicon can exhibit metallic properties when one reaches the dopant concentration of the order of 2.2 × 1019 cm−3 [32, 33]. In assuming doped silicon with high conductivity, the classical Drude model approach is used under the free carrier approximation. It estimates dielectric constant values reasonably well, justifying our approximation. For instance, for the heavily doped silicon considered (dopant concentration ~1020 cm−3;  mΩcm), the model predicts a relative permittivity of at 1 THz, which is large enough to be accurately modeled as a perfect electric conductor. For comparison, the relative permittivity of Al at this frequency is calculated to be . However, recent investigations have questioned its validity. It is believed that a modified model taking into account the electron-electron and electron-ion interactions may explain discrepancy [34]. In heavily doped silicon, one can possibly design and fabricate even complex geometries and accomplish plasmonic response [35]. In this direction, waveguides with geometries such as long V-shape structures [36] and concave pyramidal troughs have been demonstrated more recently [37]. Further, recently Kumar et al. demonstrated that a one-dimensional array of pillars forms a plasmonic waveguide and supports highly confined terahertz modes [32]. In certain applications such as modulators, sensors, and resonators, one may be required to have more interaction of the coupled THz wave with the resonant structures in order to achieve more sensitivity and responsivity; therefore, two- or more dimensional configurations could be effective. Further, a more pronounced resonant plasmonic behavior is anticipated in multidimensionally structured configurations. Therefore, a thorough investigation in this direction is prerequisite.

In this paper, the terahertz radiation propagation on a boron doped silicon surface structured with subwavelength scale square and rectangular shape vertical structures is examined both numerically and experimentally. The vertical structures are patterned in two dimensions and consequences of these structures on 1D terahertz surface plasmon propagation are investigated. A number of devices comprising vertical structures are fabricated with various different geometrical parameters. The THz time domain spectroscopy, a powerful tool for getting both amplitude and phase information, is used to characterize the THz transmission response. The patterned surfaces act as an effective media for the propagation of terahertz surface plasmons at a particular frequency. The antiresonant frequencies of the terahertz radiation modes are found to be dependent on the geometry of the structures and hence a desired frequency mode can be made to propagate by properly choosing the geometrical parameters of the structures. The effect of change of periodicity in both dimensions is examined numerically. The results are validated with a theoretical model. The paper is divided as follows. In Section 2, the fabrication of the samples and experimental setup, basically a modified THz time domain spectrometer used to characterize the samples, are discussed. In Section 3, experimental observations in conjunction with the numerical simulations are examined. This section also presents a detailed study on the behavior of antiresonant frequencies of the fundamental modes with geometrical parameters of the structures. The results are summarized in Section 4.

2. Experimental Details

The experimental part mainly consists of fabrication and characterization. The 2D vertical structures are fabricated on a heavily doped silicon substrate using advanced fabrication facilities. The silicon when heavily doped exhibits metallic properties making it a suitable candidate for the plasmonic guided wave applications at terahertz frequencies. The waveguide structures with subwavelength shape vertical rectangular structures are fabricated on a highly -doped silicon surface (carrier density ~5 × 1020 cm−3) on an area of 5 cm × 4 mm. The photo lithographically patterned silicon surfaces are etched through deep reactive ion etching (DRIE) with an etch rate of 1.5 μm/s. The schematic diagram of a patterned waveguide is shown in Figure 1(a). The inset shows a bigger view of the pattern with all the parameters being clearly defined. In the figure, “” is the length, “” is the width, “” is the height, and and stand for the periodicity of the vertical structures in - and -directions, respectively. A scanning electron micrograph (SEM) of a portion of a fabricated waveguide pattern is shown in Figure 1(b). At one end of the patterned surface, a rectangular groove is fabricated in a second step of photolithography. The rectangular groove is patterned with a width to length () ratio of . The groove is again etched using deep reactive ion etching to a depth of 100 μm. The structured waveguide patterns for their plasmonic response are characterized using modified THz time domain spectroscopy setup. The schematic of the experimental setup is shown in Figure 2. An optical beam at 800 nm wavelength from the amplified laser system is split into 80 : 20 in power for use as pump and probe beam, respectively. The femtosecond pulsed beam from the amplified laser system with average power of ~1.2 W, repetition rate of 1 kHz, and temporal pulse duration of ~43 fs is used for the generation of terahertz through a 1 mm thick InAs via Photo-Dember effect. The s-polarized pump beam incident on a nonlinear crystal leads to broadband terahertz radiation generation which is collected and collimated using off-axis parabolic mirrors. The patterned substrate in the setup is placed in such a way that the terahertz beam from the second parabolic mirror focuses on the coupling groove fabricated at one end of the structured waveguide. The spot size of the THz beam focusing on coupling groove is approximately 3 mm. By hitting free space THz at an angle, one may anticipate a free space signal arriving at the detection crystal which can interfere with the surface plasmons. In order to minimize this, we used a metal film in the middle of waveguide held at a gap of ~2 mm above the surface as shown in Figure 2. One may always use more promising approaches to avoid unsought free space signal. This includes a curved sheet for separating free space THz from the terahertz surface plasmons as proposed by Jeon and Grischkowsky [38]. The incidence of input terahertz at 90° with respect to coupling groove can also help in minimizing the unsought THz signal arriving at the detection crystal [37]. The terahertz wave once coupled to the waveguide has the electric field component pointing normal to the surface and magnetic field along the surface in the -direction. The terahertz surface plasmon polaritons (SPPs) after propagating along the patterned structures are detected by a (110) ZnTe crystal of thickness 1 mm with crystal face oriented orthogonal to the -axis. We used the electrooptic sampling method to measure the specific vector component of the SPPs. The detection system utilizes differential detection for better sensitivity.

3. Results and Discussions

In Figure 3, both the measured and simulated waveguide transmission spectra are shown for the vertical subwavelength scale structures designed across the 5 cm × 4 mm silicon substrate. The experimental results are obtained using THz-TDS spectroscopy discussed above. The corresponding structures are numerically simulated using finite element method based electromagnetic simulation package, that is, CST Microwave Studio. Figure 3(a) shows the frequency domain spectra corresponding to the patterned vertical rectangular structures having parameters: length “” = 300 μm, width “” = 200 μm, lattice constant “” = “” = 400 μm, and height “” = 100 μm. The blue traces correspond to the experiment, while green traces correspond to the simulation. A discrete THz source was used at one end of the waveguide for simulation and probe detector was placed at the other end. The results from both the simulation and experiment are in fairly good agreement and are found to support resonant modes. Both spectra do show a resonant mode at 0.36 THz, which we call as the fundamental plasmonic mode supported by the structured substrate. The experiment does show a higher order mode as well at 1 THz frequency which is missing in the simulated spectra. Further, it is noted that experimental spectral features are broader than the simulated resonant features. It is important to mention that, in simulations, highly doped silicon is considered as perfect electrical conductor (PEC) which is a good approximation at THz frequencies; however, in practice, the assumption may slightly differ, resulting in broader spectral features in the experiment. Further, geometrical defects and greater losses, namely, scattering, diffraction, and so forth, associated with the fabricated structures, add to the discrepancies. One may also see that, on the higher frequency side, experimental observations do not comply with the simulation results. This is believed to be due to the inaccuracy in the simulation at the shorter wavelengths resulting from the smaller mesh size which is of the order of . Further, greater quantization errors are anticipated at the higher frequency side that is beyond the first Brillouin boundary which is given by . In our case, this is calculated to be 0.38 THz. The experimental and simulation observations were further verified theoretically. For the fundamental TM-like mode, neglecting the diffraction effect under the assumption of , the dispersion relation is given by [31] where is the wave number in the direction of propagation. is the free space wave number, is the height of the pillars, and is the ratio of width to periodicity for the vertical structures. We calculated the antiresonant frequency at for the “” = 300 μm, “” = 200 μm, “” = “” = 400 μm, and “” = 100 μm and it turned out to be 0.36 THz, which is same as the observed antiresonant frequency of the fundamental mode. Later in this paper, we will see that a variation of periodicity “” does not affect the fundamental resonant frequency and hence the analytical modal given by (1) accurately predicts the frequency of the fundamental terahertz mode.

In Figure 3(b), another patterned configuration of vertical rectangular structures with a height of “” = 300 μm is examined. Other parameters of the structured substrate are same as the previous pattern. Again, it is apparent from the figure that both the experimental and simulation spectra exhibit resonant behavior and antiresonant frequencies of the fundamental mode are found to be in good agreement. The blue traces correspond to the experiment; however, red traces represent numerically simulated frequency domain spectra. The antiresonant frequency of the fundamental mode is observed to be 0.19 THz. One may see the higher order modes as well in the figure which are more pronounced in case of numerical simulation. The discrepancy in the experiment and simulation on the higher frequency side can be understood again from the explanation given for earlier plot.

In Figure 4(a), the effect of variation of height “” of the vertical square shape structures on the terahertz surface plasmon propagation is examined numerically. The figure shows the numerically calculated frequency domain spectra for the periodically patterned structured surface for the different heights of the vertical structures. The periodicities and are the same, that is, 400 μm in both directions. The other parameters are “” = 200 μm and “” = 200 μm. The red traces in the figure correspond to the = 100 μm, while blue, green, and brown color traces represent frequency domain spectra corresponding to the vertical structures of μm, 300 μm, and 400 μm, respectively. One may note that the amplitude of the waveguide transmission decreases as the height of the structures increases. Hence, one is required to fabricate structures with shorter height for the efficient transmission; however, frequency should be kept in mind. Further, it may be noted that as the height increases the antiresonance frequency goes down. It can be explained from (1) which clearly shows a dependence of the dispersion relation on the height of the structures. The frequency of saturation in the dispersion relation can be approximated with antiresonant frequency of the fundamental mode which changes with height. One may also note that waveguide geometries with taller structures would result in lower group velocity of the fundamental modes and the lower cut-off frequencies. Therefore, one can design 2D plasmonic waveguides keeping in mind the amplitude as well as the frequency of the terahertz waves. In Figure 4(b), an electric field profile across the structured surface in the plane is shown at 0.36 THz. One may note that there is a strong field confinement along the vertical square shape structures at the resonance. The field profile at the nonresonant frequency was also examined (not shown here) and no confinement along the surface was seen. In Figures 5(a) and 5(b), terahertz amplitude transmissions for a patterned substrate with different periodicity in the - and -directions are shown. Unlike the groove patterns as examined in the literature in the past [25], the periodicity strongly influences the resonant mode propagation in the -direction as shown in Figure 5(a). Figure 5(a) shows red, green, and blue trances which correspond to the = 300 μm, 400 μm, and 500 μm, respectively. The other parameters of the structures are “” = 200 μm, “” = 200 μm, and “” = 400 μm. With smaller periodicity, the coupling between the structures is stronger and one see a sharp resonance as can be seen in the red color. When the periodicity increases, coupling between successive structures becomes weaker and resonant effect is not pronounced. One may observe a broadband response of the terahertz radiation propagation with periodicity on the order of wavelength scale. However, the amplitude transmission is very small in that case and there is poor confinement for the terahertz radiations. In Figure 5(b), effect of periodicity variation in the -direction is examined on the resonant behavior of the patterned surface. It is observed that the fundamental resonant frequency is not influenced by the periodicity in this case. One may certainly observe a change in the amplitude of the fundamental resonance with a change in the periodicity. As periodicity is increased, the amplitude goes down. This is because of the weak coupling between the structures and hence lower contribution to the surface plasmon propagation in the -direction from the structures in the -direction. The antiresonant frequencies are caused by the Fano interference of discrete and continuum spectrums. The TM mode propagation in the -direction and associated cut-off frequency are dependent on the periodicity in the -direction as can be noticed from (1). For THz surface plasmons, field in the -direction is negligible and periodicity is not expected to change the antiresonant frequency. Further, the effect of the lateral length (in -direction), that is, “,” and the width “” of the structures on the propagation properties of the terahertz radiation propagation is investigated. The numerically examined results are shown in Figure 5(c). The red triangular dots represent the numerically calculated antiresonant frequencies for different values of the width “” of the vertical square shape structures. The other parameters of the structures are “” = 200 μm, height “” = 100 μm, and “” = “” = 400 μm. The calculated data points are fitted with the straight line equation which is represented in the red color. One may note that straight line fits fairly well with the data implying that a variation in the width does not leave any effect on the antiresonant frequency of the fundamental mode. Similarly, the effect of length is examined in the same figure. The blue circular dots represent the antiresonant frequencies of the fundamental mode for various different lengths of the rectangular structures. The width and other parameters are kept constant in this case and are taken as “” = 200 μm, height “” = 85 μm, and “” = “” = 400 μm. The data points are again fitted with the straight line equation which is represented by the blue line. Again it fits the data very well and hence one can infer that change in the length of the vertical structures does not leave any effect on the antiresonant frequency of the fundamental mode. In Figure 5(d), a comparison between the transverse confinement of THz surface plasmon modes in 1D and 2D structured waveguides is shown for parameters: “” = 200 μm, “” = 200 μm, height “” = 100 μm, and “” = “” = 400 μm. Using 3D numerical simulation software, the field was measured in the transverse direction (-axis) in the middle of the waveguides using discrete detection ports. The data points are exponentially fitted. In case of 1D sample, the decay length was observed ~420 μm, while in 2D sample it came out ~550 μm. Therefore, transverse confinement of the terahertz surface plasmon decreases when waveguide is patterned in both spatial directions. The information can be used to design the plasmonic metamaterials based terahertz active and passive components and can prove to be the building blocks in developing the terahertz networks.

4. Conclusions

In this paper, 2D periodically patterned vertical structures on a silicon substrate are examined for terahertz surface plasmon polariton propagation. Heavily doped silicon substrates are periodically patterned with 2D array of structures resulting in plasmonic waveguide configurations. The experimental measurements through a THz time domain spectrometer show resonant terahertz mode propagation at certain frequencies defined by the geometrical parameters of the structures. Numerically simulated structures with similar shape and sizes yield results which are consistent with the experimental observations. An increase in the height of the structures results in a decrease in the amplitude as well as antiresonant frequency. The alterations in width and length of the structures do not affect antiresonant behavior. In order to understand the effect of structures in 2D on the THz surface plasmon propagation properties, periodicity is changed in both spatial directions and resonant behavior is examined. An alteration in periodicity in -direction significantly affects the frequency and amplitude of the resonant mode; however, cut-off frequency remains unchanged in other spatial directions. Further, a comparison in the transverse confinement of the 1D and 2D structured samples is made. The decay length of the THz surface plasmons increases under 2D patterning signifying a decrease in the transverse confinement. The study is anticipated to greatly felicitate the designing of futuristic devices and networks on silicon surfaces capable of operating at terahertz frequencies.

Conflict of Interests

The author declares that there is no conflict of interests regarding the publication of this paper.

Acknowledgment

The author would like to thank the University of Maryland for the characterization and fabrication facilities.