## 3 Results
### Analysis for Scenario without a Dielectric Layer  
In Figure 6 the bias markings of the center column in different directions can be seen. In order to investigate the effect of electromagnetic coupling for eSRR in detail, the electromagnetic characteristics of the eSRR device with different values of U for zero bias Q are discussed. The effect of electromagnetic coupling was examined in the TE and TM modes. At Q = 0 \(\upmu\)m and U = 0 \(\upmu\)m, in TM mode the resonant frequency is 1.934 THz, as shown in Figure 7a. When changing the value of U to 2 \(\upmu\)m there is increase of this frequency to 1.935 THz. At U = 6 \(\upmu\)m there is the lowest resonant frequency equal to 1.931.5 THz. Then, by increasing the shift to U = 7 \(\upmu\)m, the resonant frequency has value 1.935.5 THz. After that there is a small increase at U = 8 \(\upmu\)m to 1.933 THz. Increasing the U value to 11 \(\upmu\)m shows increasing resonance frequency value to 1.935 THz and then there are only two bigger increases, at U = 13 \(\upmu\)m with frequency 1.936 THz and U = 15 \(\upmu\)m with frequency equal to 1.937 THz. In this mode, the stability of the resonant frequency to the change in column offset can be seen, compared to other figures. At Q = 0 \(\upmu\)m and U = 8 \(\upmu\)m, the resonant frequency is 1.933 THz. By shifting the I-shaped metal bar to the left and thus manipulating the Q value by 2 \(\upmu\)m, it can be observed, as in Figure 7b, that there is a small change of the resonant frequency with value equal to 1.935 THz. At the value of Q = 5  
Figure 5: Figure showing the planes of: (**a**) perfect electric conductor, perfect magnetic conductor; (**b**) port, impedance boundary condition, transition boundary condition. in TM mode two or even three resonances appeared with most geometry alignments and were sensitive to geometry change, in TE mode with some geometries there were multiple resonances that began to overlap. The readings of the main resonant frequency that appeared with each eSRR geometry was more complicated, which is not to say that this resonance maintained its stability longer than with TM mode. In the results shown in Figure 8 the resonances dependent on the change of position were included and those that showed their presence at each Q and U value. In Figure 8, the change in resonant frequency with tunable eSRR geometry in TE mode is shown. As can be seen in Figure 8a, the resonant frequency for Q = 0 \(\upmu\)m and U = 0 \(\upmu\)m was 1.949 THz. The resonant frequency, when shifted to a height of 2 \(\upmu\)m, has increased to 1.96 THz. In this mode, the sensitivity of the resonant frequency to a change in the geometry of  
Figure 8: Resonance frequency values while changing the I-shaped position in TE mode: (**a**) change in resonance frequency value when moving the center I-column metal bar up in 2 \(\upmu\)m increments to a maximum excursion of 16 \(\upmu\)m; (**b**) change in resonance frequency value when moving the center column right and left in 2 \(\upmu\)m up to a maximum excursion of 6 \(\upmu\)m for a rightward shift and \(-6\)\(\upmu\)m for a leftward shift. Shifting in this direction is done for a value of U equal to 8 \(\upmu\)m, i.e., at the center position of the moving I-shaped column. the unit cell can be seen. By changing the geometry, the resonant frequency increased to 1.993 THz at Q = 0 \(\upmu\)m and U = 9 \(\upmu\)m, U= 10 \(\upmu\)m. By increasing the position to Q = 0 \(\upmu\)m and U = 12 \(\upmu\)m, resonant frequency began to fall to the value of 1.99 THz. The decrease in resonant frequency has continued and in the position of Q = 0 \(\upmu\)m and U = 16 \(\upmu\)m reached 1.97 THz. Figure 8 only contains the values of one resonance, which appeared before shifting, at the zero position of the geometry. There are more than one, two, or even three resonances in the TE mode. When moving the I-shaped metallic bar, over a height of 10, four resonances started to build up, which even had a similar amplitude.