2. Scattering from Bi-Isotropic Metasurfaces
In this section, we describe the scattering performance of an omega-type bi-isotropic metasurface that is illuminated by a normally incident plane wave. The scattering analysis of bi-isotropic metasurfaces provided in this section follows that introduced in [
5]. We consider a metasurface at a planar boundary between two regions of space, as shown in
Figure 1, where the intrinsic wave impedance of region 1 is
and of region 2 is
. The metasurface is at the
plane separating the two regions, and it is illuminated by normally incident plane waves.
The interaction between the metasurface and an illuminating plane wave can be described via scattering parameters (S-parameters), which comprise the ratio between the scattered plane wave electric field and the incident plane wave electric field. The ratio of scattered electric field in region
n to the incident electric field in region
m for different polarizations is given as a
matrix.
When viewed from region 1,
is the reflection coefficient and
is the transmission coefficient. Similarly, when viewed from region 2, the reflection coefficient is
and the transmission coefficient is
.
In general, a bianisotropic metasurface can be modeled as a two-dimensional array of polarizable particles [
4]. For time-varying illuminating fields, the polarizabilities can be effectively characterized using equivalent surface impedances [
5,
6,
7,
8,
9,
10,
11]. The equivalent surface currents can then be related to the averaged, tangential electric, and magnetic fields using surface parameters that are represented as
tensors: the electric sheet admittance tensor
, the magnetic sheet impedance tensor
, and the magneto-electric coupling tensors
and
. With these parameters, the electric and magnetic surface currents that are induced on the metasurface can be related to the average tangential fields and compared to the boundary conditions across the metasurface.
The variables
,
,
, and
relate the
x- and
y-polarized averaged field components to the
x- and
y-polarized current density components that are induced on the metasurface. The various electric field vectors are
and the magnetic field vector is
(the surface current quantities
and
are similarly defined), where the averaged fields are
and
. The electric admittance tensor is defined as
with the other parameters being similarly defined.
For a reciprocal metasurface,
,
, and
[
32]. Imposing isotropy on the surface parameters results in
where
is the
identity matrix and
.
Restricting the metasurface to omega-type bi-isotropy precludes polarization conversion by the metasurface. Therefore, the response for each polarization is identical. This allows us to analyze the metasurface as a two-port network for a single polarization, rather than as a four-port network when all of the polarizations were considered. The two-port S-parameters relate the electric field of the incident and reflected plane waves as
In order to calculate the S-parameters of the metasurface, consider an
x-polarized plane wave, as shown in
Figure 1. Assuming that the surface is isotropic, the boundary conditions of Equation (
2) simplify to
From Equations (
6) and (7), we obtain four equations by considering the illumination from region 1
and region 2
separately. These four equations relate the S-parameters to the surface parameters of the metasurface. In each case,
and
are expressed in terms of the S-parameters and the illuminating electric field. Additionally, the assumption of plane wave illumination allows us to express the magnetic field quantities in terms of the electric field and wave impedance of each region. These four equations are simplified and assembled into a matrix equation to express the surface parameters in terms of S-parameters.
The form of Equation (
8) is convenient for calculating the surface parameters that will implement the desired S-parameters. However, re-arranging Equation (
8) and simplifying to solve the S-parameter quantities provides
In the case of a lossless metasurface, the surface parameters
and
are purely imaginary, while
are real quantities [
32]. In this case, Equations (9)–(12) can be further simplified, as shown in Equations (14)–(17). Note that
only when
, and
only when
and
.
We can also determine the limitations that are placed on the S-parameters when passive, lossless, and reciprocal restrictions are enforced. For a bi-isotropic metasurface, the S-parameters represent a two-port network, as described in Equation (
5). Each element is a complex number, so there are eight total variables (four real and four imaginary quantities). For a reciprocal network,
when the port impedances are the same. This relationship shows that both of the transmission coefficients are the same in amplitude and phase. However, a different relationship is needed for the case of the bi-isotropic metasurface, since the port impedances are different. Reciprocity is satisfied when the transmission phase shift and transmitted power are the same for each direction of illumination. When the port impedances are not equal, the electric field amplitude will change, depending on the wave impedance of the medium in order to satisfy the reciprocity conditions, so
.
In order to determine the reciprocity relationship for a bi-isotropic metasurface, consider two cases: (i) where the metasurface is only illuminated from region 1 and transmitted power is determined in region 2, and (ii) the metasurface is only illuminated from region 2 and transmission measured in region 1. By equating the transmitted power in both cases, we arrive at
While Equation (
19) provides a relationship between the transmission coefficient magnitudes, reciprocity also requires that the transmission phase be the same. Applying this and assuming the wave impedance of each region is real, we arrive at
Note that Equations (15) and (16) satisfy this relationship, since the metasurface parameters were restricted to be reciprocal.
In order to enforce the lossless condition, the time-average power that is absorbed by the metasurface must be zero. This is calculated as
By applying the plane wave relation between the electric and magnetic fields, and expressing
and
in terms of the S-parameters from Equation (
5), Equation (
21) becomes
where
Simplifying Equation (
22) and utilizing the reciprocity relationship in (
20) results in three equations that must be satisfied in order to implement a lossless and reciprocal metasurface.
These three equations under-define the six independent scattering matrix variables. Consequently, three variables can be chosen freely without violating the lossless and reciprocal conditions. Specifically, Equations (24) and (25) provide the ability to choose one amplitude of the scattering matrix. If
is chosen, as is commonly the case, then
, and a phase constraint is obtained from Equation (26)
where two phase shifts of the S-parameters can be freely chosen.
Therefore, for a bi-isotropic metasurface to be both lossless and reciprocal, three degrees of freedom exist in its S-parameters: one S-parameter amplitude and two S-parameter phases. These three degrees of freedom are set through the design choices of the metasurface. It is worth recalling that enforcing lossless and reciprocal behavior in the surface parameters for the bi-isotropic metasurface results in , , and . Thus, three distinct surface parameters can be chosen to achieve three desired scattering properties.
3. Bi-Isotropic Metasurfaces: Bandwidth and Quality Factor
In practice, bi-isotropic metasurfaces typically rely on resonant structures to produce the strong field interactions that are required to perform the desired field transformations. However, the use of resonances places inherent limitations on the bandwidth. In this section, the relationship between matching networks and bi-isotropic metasurfaces is considered, and the quality factor of a metasurface realized using three impedance sheets is defined. We demonstrate that the quality factor can be used as a metric to predict the metasurface’s bandwidth and identify unit cells that degrade the performance of inhomogeneous metasurfaces.
We consider the following example to understand the relationship between impedance matching networks and bi-isotropic metasurfaces. Suppose that there is a planar interface between air and alumina (
), as in
Figure 2, and the goal is to maximize the power that is transferred across the interface. Because the intrinsic wave impedances of the media are real, this amounts to minimizing the amplitude of the reflected wave. To do this, the input impedance of the metasurface must be equal to the wave impedance of the incident wave,
. Because the two media have different wave impedances the metasurface must transform the wave impedance of the transmitted wave,
, to that of the incident wave,
. In this scenario, the metasurface acts as an impedance matching layer. Here, the impedance matching layer is analogous to an impedance matching network from circuit theory, like an L or T-network, as shown in
Figure 3. From circuit theory, it is known that a complex load impedance can be matched to a complex source impedance using either an L-, T-, or
-network [
33]. The L-network contains two degrees of freedom allowing for the real and imaginary components of the input impedance to be matched. For an L-network, the solution is unique (all the degrees of freedom are used) and no other characteristics of the impedance match, such as its bandwidth or the transmission phase, can be controlled. Adding a third degree of freedom to the L-network produces a T- or
-network. This additional degree of freedom can be used to control the bandwidth or the transmission phase. Bi-isotropic metasurfaces are like T-matching networks for fields [
34]. They have three degrees of freedom that allow for impedance matching with phase or bandwidth control [
35]. To illustrate this idea, we consider a metasurface that impedance matches a normally incident plane wave on an air-alumina (
) interface over a maximum bandwidth, as shown in
Figure 2.
In order to design the impedance matching metasurface, recall that a bi-isotropic metasurface can be viewed as a two-port network that controls one scattering amplitude and two scattering phases. Therefore, designing a lossless, reflectionless, and bi-isotropic metasurface is equivalent to designing a lossless two-port impedance matrix (Z-matrix) that impedance matches a load impedance
to a source impedance
with an arbitrary transmission phase
[
36]. Consider a general lossless two-port Z-matrix in order to determine the Z-matrix that provides the desired functionality,
Imposing the impedance boundary conditions and enforcing power conservation on (
28) produces the following system of equations,
where
and
. Splitting (
29) into its real and imaginary components allows for the elements of the Z-matrix to be solved in terms of
,
and
,
From (
30), it is clear that the required two-port network is reciprocal, since
, and it has three degrees of freedom.
Three cascaded sheet impedances, as shown in
Figure 4, can be used to realize a metasurface with a Z-matrix given by (
30), as in [
5]. Expressing
Figure 4 in terms of its Z-matrix, and solving for the necessary impedance sheets, results in the following expressions for the sheets in terms of the elements of (
30),
where
is the determinant of the Z-matrix and
and
are the wavenumber and wave impedance of the dielectric spacers, respectively. Once the input and load impedances, spacer thickness, and the transmission phase are specified, (
30)–(
33) can be used to determine the necessary impedance sheets to implement the metasurface.
In order to maximize the bandwidth of the impedance match, a method is needed for comparing the metasurface’s bandwidth for different transmission phases. Here, an expression for the metasurface’s quality factor as a function of the transmission phase is derived for this purpose. The quality factor of a three-sheet metasurface is defined as,
where
is the angular resonant frequency,
is the average electric energy stored in the network at
, and
is the power dissipated in the network. In order to calculate the quality factor using (
34), the impedance sheets (
31)–(
33) are expressed in terms of lumped capacitances and inductances. The dielectric spacers in the metasurface are assumed to be electrically thin, so they can be modeled as lumped
-networks. Therefore, if the dielectric spacers are electrically thin and the source and load impedances are purely real, then the quality factor of the metasurface can be expressed as
where,
, and
is the capacitance of the
ith impedance sheet (if the sheet is inductive, then
). If the load impedances are not purely real, then the imaginary part of the load can be absorbed into either
or
, and (
35) can still be used. The quality factor,
Q, of the metasurface will be used to approximate the fractional bandwidth,
, where
is the 3 dB bandwidth of each unit cell. However, due to the presence of multiple resonances this approximation is only valid when the resonances are well separated in frequency.
The quality factor expression (
35) can be used to maximize the bandwidth of an impedance matching bi-isotropic metasurface. For a normally incident plane wave, the relevant impedance is
= 377
. Let us assume that
= 123
, and the spacers are free-space with a thickness
. Using (
35) to calculate the quality factor and the fractional bandwidth versus transmission phase produces
Figure 5.
Figure 5a plots the quality factor which is minimized at a transmission phase of
.
Figure 5b shows that this transmission phase is predicted to maximize the bandwidth, and
Figure 5c shows that it produces the maximum 3dB bandwidth. The metasurface with this transmission phase is composed of the following impedance sheets:
=
468.9
,
=
=
641.9
, and
=
=
j38.5 k
. The metasurface performance is simulated in Ansys HFSS while using dispersive impedance sheets that correspond to the following lumped elements:
= 33.9
,
= 24.8
, and
= 612.7
.
Figure 6 shows the transmission magnitudes from this simulation, where they are compared to a quarter-wave transformer and the bare interface without any impedance matching. The metasurface has a size and bandwidth comparable to a quarter-wave transformer. However, it does not require the realization of a medium with the dielectric constant
, which can be heavy and challenging to manufacture.
In practice, a metasurface’s impedance sheets are typically realized using subwavelength metal or dielectric patterning on support structures. Here, we will consider the impact this has on the bandwidth of the metasurface. First, we will consider the effect of using subwavelength patterned sheets. Subwavelength unit cells that are non-resonant exhibit a response of either a capacitive or inductive sheet impedance [
37,
38]. This indicates that modeling the patterned sheets as impedance sheets should not result in a significant bandwidth reduction when the metasurface is realized in practice. However, it may be necessary to modify the design to include additional impedance sheets to avoid extreme impedance values that are difficult to realize in practice. If impedance sheets cannot be realized at the design frequency due to manufacturing difficulties, an alternative design approach may be required, such as using detuned resonant elements. Their responses will be more narrowband.
We will also consider the effect of using a dielectric spacer as the support structure. For a metasurface that is designed using a non-magnetic dielectric spacer with a relative permittivity
, its quality factor is given by,
If the dielectric spacer is electrically thin, then the sheet capacitances can be approximated as,
when the sheet impedance
is capacitive (see
Appendix A). Otherwise, the impedance sheet is inductive and it can be ignored in the calculation of the quality factor. Additionally, when the spacer is electrically thin,
does not depend on the dielectric constant (see
Appendix A). Accordingly, the only terms in (
36) that depend on the permittivity are the capacitance terms
,
, and
. When considering the terms
,
, and
individually with (
37)–(
39) in them, it becomes apparent that the quality factor is unaffected by the dielectric spacer for capacitive impedance sheets. However, if any of the impedance sheets are inductive, then the dielectric spacer will increase the quality factor, thereby reducing the metasurface’s bandwidth.
In addition to bandwidth information, the quality factor also provides information that can guide the design of inhomogeneous metasurfaces where local periodicity is assumed. Obtaining good performance from a metasurface that is designed assuming local periodicity requires that neighboring unit cells produce fields that are approximately the same, i.e., the fields vary smoothly along the surface without large discontinuities in the amplitude or phase. In this work, it has been found that the quality factor and its first derivative with respect to transmission phase can help the designer to select unit cells that satisfy the assumption of local periodicity.
The quality factor, as given by (
35), is divergent at transmission phases near
, and
, indicating that the unit cells that are required to achieve these transmission phases possess large quality factors. Large quality factors are associated with strong resonances that are sensitive to perturbations in the surrounding environment and are lossy when realized in practice. Therefore, these unit cells should be avoided. Additionally, areas where (
35) is not smooth (i.e., points where the first derivative is discontinuous or undefined) indicate transmission phases where the reactance of at least one of the impedance sheets changes sign. These points should also be avoided because they identify transmission phases where the required reactance values display asymptotic behavior. This introduces rapid variations in the values of the impedance sheets and fields in the metasurface that invalidate the assumption of local periodicity.
In order to see how this information can be used, consider a metasurface embedded in free-space that refracts a normally incident plane wave to
at a frequency of
= 10 GHz. This requires a gradient metasurface: an inhomogeneous metasurface that imposes a linear phase gradient on an impinging wave-front to produce reflection or refraction in a desired direction [
39]. Refraction requires the metasurface to alter the transverse wavenumber of an incident plane wave (
) to produce the desired refracted wavenumber (
), where
k is the wavenumber in the surrounding medium. Therefore, the metasurface must impart transverse momentum that is equal to
. Practically, this is realized by discretizing the metasurface into
N sub-wavelength unit cells that are of size
, each possessing a transmission phase
, such that
. Each unit cell must be reflectionless in order to maximize the transmitted power into the refracted wave. This means that impedance matching and phase control are required, so (
31)–(
33) can be used to design the unit cells of the metasurface.
For this example, the metasurface will have 10 unit cells per transverse wavelength (in free-space) and the spacers will be assumed to be free-space with a thickness
. As a first pass at the design, the metasurface is designed to impose a linear phase gradient with the unit cell transmission phases that are shown in
Table 1. The required sheet impedances, as shown in
Figure 7b, are solved using (
31)–(
33) and one period (10 unit cells) of the metasurface is simulated in COMSOL using periodic boundary conditions.
Figure 7c shows the results.
The metasurface that is designed using this phase gradient exhibits significant reflections and the transmitted wave is not purely refracted. A slight perturbation of the linear phase gradient can be used to improve the performance. The appropriate perturbed phase gradient is found using the quality factor and its first derivative with respect to transmission phase. To find the problematic transmission phases in the original design, plots of the quality factor and its first derivative are shown in
Figure 8. By inspecting the plots, four unit cells with problematic transmission phases are identified: 1, 5, 6, and 10. Unit cells 1, 5, and 6 are problematic, because they are near points where (
35) is not smooth, and unit cell 10 is problematic due to its large quality factor. The problematic transmission phases are adjusted to improve the performance of the metasurface, as shown in
Table 1 and
Figure 9. These phase shifts reduce the maximum unit cell quality factor by approximately 10 and force the reactance of each impedance sheet to change sign only once at
.
The metasurface is redesigned with the modified transmission phases and
Figure 7b shows the required sheet impedances. Ten unit cells of the metasurface are again simulated in COMSOL using periodic boundary conditions, and the results are shown in
Figure 7d. We see that the redesigned metasurface performs significantly better than the analytical design. This indicates that avoiding transmission phases that require a large quality factor or exist near non-smooth or asymptotic regions of
can improve the performance of gradient metasurfaces that are designed using the local periodicity assumption.
Violations of local periodicity (like those discussed above) can present challenges when realizing inhomogeneous metasurfaces where local periodicity has been assumed. Issues that arise from these violations have been implicitly handled in the literature in a variety of ways. In [
7], the phase gradient was altered to improve the metasurface’s performance by reducing the transmission losses. On the other hand [
36,
40,
41], made the spacers between the sheets extremely thin
. This generally increases the quality factor of the unit cells, but it has the benefit of shifting the transmission phases where all three impedance sheets transition from capactive to inductive to occur at the same point. This means that shrinking the spacings makes it easier to select transmission phases that avoid regions where (
35) is not smooth. Consequently, extremely thin spacings can improve the design performance at the expense of increasing manufacturing difficulties and producing higher quality factors: lower bandwidths. Alternatively, PEC [
42] or PMC [
29] baffles have been used to eliminate inter-cell coupling to validate the assumption of local periodicity. However, in practice, the use of PEC baffles presents a manufacturing challenge and PMC baffles cannot be realized. These examples indicate a design trade-off between manufacturability and performance when realizing inhomogeneous metasurfaces. Using the quality factor, as shown in this section, provides an alternative way to improve design performance. It can be used to systematically identify the problematic unit cells and adjust them where possible to allow for the trade-off between performance and manufacturability to be balanced. An alternative to this approach is to avoid the assumption of local periodicity and model interactions between unique unit cells through homogenization and integral equations, as reported in [
43].
4. Scattering from Bianisotropic Metasurfaces
While the scattering of plane waves was analyzed for the simplified case of a bi-isotropic metasurface shown in
Section 2, it is worthwhile to consider the general case where isotropy is not assumed. Following the general process of
Section 2, the boundary conditions of Equation (
2) can be expressed in terms of the S-parameters (where each S-parameter term is the
matrix of Equation (
1)). Note that
. Expressing the surface parameters in terms of the S-parameters gives [
5]
Equation (
40) can also be re-arranged so that the S-parameters are expressed in terms of surface parameters.
Analyzing the degrees of freedom helps to determine the number of surface parameters required to realize a specificed S-matrix, as in the bi-isotropic case that is discussed in
Section 2. In the bianisotropic case, both of the polarizations are taken into account, which leads to a
scattering matrix and 16 complex numbers as its entries. In most cases, reciprocity is desired for metasurfaces, which results in a symmetric S-matrix (assuming the port impedances are identical)
which indicates that only 10 out of the 16 entries are actually independent. Because each complex number contains its real part and imaginary part, there are 20 free variables in total under the reciprocal condition. The S-matrix also has to be unitary if we further require the metasurface to be lossless:
By incorporating the reciprocal condition (
42) into the lossless condition (
43), one can expand (
43) into 10 different equations, which impose 10 additional restrictions on the 20 free variables. Consequently, for a reciprocal and lossless bianisotropic metasurface, there are 10 degrees of freedom in total.
A similar conclusion can be drawn by considering the surface parameters. Recall that
,
, and
for a reciprocal metasurface [
32]. There are three free entries in both of the
and
matrices, and four free entries in the
or
matrix. Moreover, for the metasurface to be lossless,
and
must have purely imaginary entries, while
and
must have purely real ones [
32]. Again, it can be seen that the total degrees of freedom of the system is 10.
In practice, several sheets are usually cascaded and separated by dielectric spacers to form bianisotropic metasurfaces. Typically, these sheets only possess electric responses that are characterized by admittance tensors
, since they can be readily realized using metallic patterns. For bi-isotropic metasurfaces, or in the case where only a single polarization is of concern, we know that three sheets are enough to realize a specified response. However, the situation becomes more complicated for bianisotropic metasurfaces. When both of the polarizations are involved, a single lossless, reciprocal electric sheet provides three degrees of freedom under lossless and reciprocal conditions, i.e., the imaginary numbers
,
and
. Therefore, at most, four sheets are required to realize an arbitrary reciprocal and lossless bianisotropic metasurface with 10 degrees of freedom. Although many bianisotropic metasurfaces can be realized with only three electric sheets, there are some cases in which introducing a fourth sheet is necessary. Examples include the polarization rotators in [
5,
23]. The fourth sheet not only provides the required degree of freedom, but also enhances the operational bandwidth.
A network analysis technique, known as the wave matrix approach, was adopted in [
23] in order to synthesize a cascaded sheet design. Wave matrices relate the forward and backward propagating electric fields on one side of the scatterer to those on the other side. For an arbitrary scatterer that is shown in
Figure 10a, the wave matrix
is defined as:
Similar to the S-matrices, wave matrices contain information regarding the incident and reflected waves. The advantage of using wave matrices is that they significantly simplify the analysis of cascaded structures, such as ABCD matrices. The wave matrix of a cascaded structure can be obtained by simply multiplying the wave matrices of its constitutive blocks. In our multi-layer metasurfaces, these blocks include metasurface interfaces across two dielectric media and dielectric spacers, as illustrated in
Figure 10b,c, respectively. Their corresponding wave matrices can be derived from the boundary conditions, and they are explicitly shown in [
23].
The procedure for synthesizing a reciprocal and lossless S-matrix is briefly outlined here. First, the desired S-matrix,
, is stipulated based on the required application. It is then converted to a wave matrix, as follows [
23]:
where
represents a
null matrix. This wave matrix
is known and it is set as the synthesis goal. It is worth noting that, if
has a zero determinant, taking the inverse matrix in (
45) becomes invalid. In this case, a perturbation can be introduced into
to alleviate this problem. For simplicity, it is assumed that this S-matrix can be realized by cascading three sheets. Accordingly,
Figure 11 displays the targeted structure and the cascaded wave matrix that relates
to
is
in which
,
and
are the admittance tensors that need to be solved. By setting
, and with some algebraic manipulation, one can find the admittance tensor of the second sheet
[
23]:
The symbol ⊗ in (
47) represents the Kronecker product of matrices and
is some constant scalar. The matrix
contains information regarding the dielectric interface where the first sheet is located,
carries the phase information of the second dielectric spacer (
,
, and
are similarly defined), and
is a constant matrix [
23]:
Similarly, the admittance tensors of the first and third sheets
and
can be expressed in terms of
,
where
,
,
, and
are the constants explicitly calculated in [
23]. A more complicated synthesis procedure involving four cascaded sheets is also discussed in [
23], but the main idea follows the three-sheet case shown here.