
Preprint
in the form of a circuit model can become a tool for designing and predicting these behaviors. This work focuses
on microwave and terahertz frequencies, where the permittivity of the metals becomes high enough to justify perfect
electric conductor (PEC) approximation. The transmission line model is considered in this paper due to its appropriate
accuracy and high flexibility, which can be easily generalized to multilayer structures.
Modeling a waveguide structure in the form of a transmission line has been researched for a long time and has numerous
applications in obtaining scattering parameters [21, 22]. The transmission line model is used in a wide frequency range
from microwave to infrared [23
–
27]. With the help of the transmission line, the array of apertures in a PEC film has
already been investigated [28
–
31]. In the first developed circuit model, the array of holes was analyzed by considering
the thickness of the PEC film. In this case, the circuit model was created by considering only the dominant mode inside
the apertures [32, 33].
In previous works, a circuit model has been developed for rectangular apertures by considering only
T E01
as the
dominant mode [33, 34]. In addition to two-dimensional arrays, one-dimensional arrays of slits have also been
investigated by considering the
T EM
mode as the dominant mode using the circuit model [32, 35]. If the thickness of
the aperture is large enough that the excited high-order modes inside the hole are damped from the first opening to
the second opening of the hole, the model has high accuracy. In these models, if the structure’s thickness decreases,
the error increases due to high-order modes; additionally, the error increases more rapidly for circular-shape apertures
compared to rectangular apertures. Due to the widespread applications of metasurfaces and frequency selective surfaces,
we need to create models for thin arrays.
Frequency selective surfaces are thin periodic surfaces that are important both in terms of industrial applications as well
as in terms of fundamental science. The array of apertures perforated in a PEC thin film is a kind of FSS, and many
powerful circuit models have been proposed for these structures that are suitable for better design and understanding
of the behavior of these structures [36, 37]. In these structures, the circuit model is created with the help of the
spatial profile of the transverse electric field at the opening of the aperture. Previous works have provided profiles
for rectangular [37], cross [28], and annular holes [37, 38]. Initially, the circuit model was developed for a structure
consisting of a single layer of the intended FSS [25]; then, the circuit model was improved for the case where the FSS is
embedded in a layered medium [39, 40]. Structures with multiple dielectric layers and multiple FSS layers were then
examined using circuit models [28, 37], but all of these circuit models were developed regardless of the FSS’s thickness.
It should be noted that effect of thickness can cause a considerable error; however, the thickness is about one-fiftieth of
the structure’s period.
In this paper, an extended circuit model is developed, which is much faster than full-wave simulations and is much
more accurate than its predecessors. The circuit model is formulated in such a way that it takes into account high-order
modes inside the holes; thus, it has high accuracy in low thickness structures. A closed-form expression is presented for
the spatial profile of the transverse electrical field on the circular aperture which means that, unlike previous models, we
do not need to extract the field profile from the full-wave simulators. It is shown that the previous electric field profile
for the square hole has a large error when the width of the square hole is large (the width of the square is greater than
half the period) [28], and this problem is well solved by presenting the new profile. Unlike square and circular holes, an
array of PEC pillars always supports propagating mode that makes interesting properties [16], and these structures are
also analyzed analytically with the help of the proposed circuit model in this work.
The paper is organized as follows: Section 2 is devoted to describing how the proposed circuit model is developed,
and three different cases: single-mode, multi-mode, and modified single-mode circuit models are described. In section
3, the low-thickness structures are investigated. In section 4, the proposed circuit model is extended to multilayered
structures; finally, a conclusion is described in section 5.
2 The proposed circuit model
Periodic structures can be analyzed similar to other waveguide structures by considering a unit cell as a waveguide;
thus, the formulation presented can be extended to any waveguide structure, but we focus only on the square array of
sub-wavelength holes. In addition to the previous circuit models, environmental models for the desired structure have
also been investigated [16, 41]. Similarly, the isotropic/anisotropic environmental model is created only by considering
the dominant mode inside the aperture, which means that these methods also have a large error for a thin array of holes.
The following formulation is presented by considering all the modes inside the hole.
Figure 1 shows the schematic of the structures which are illuminated by a normal incident plane wave. The upper
environment of the structure is region I (
n1
), inside the holes is region II (
n2
), and the lower environment is region III
(
n3
). For simplicity’s sake, we first assume the structures depicted in Fig. 1 is semi-infinite, i.e.,
L→ ∞
. Impinging of
a plane-wave on a semi-infinite structure excites aperture modes in region II, and the Floquet modes in region I. The
2