Next Article in Journal
Compact Photonic Crystal Polarization Beam Splitter Based on the Self-Collimation Effect
Next Article in Special Issue
Electromagnetic Wave Scattering from a Moving Medium with Stationary Interface across the Interluminal Regime
Previous Article in Journal
Perspectives on UVC LED: Its Progress and Application
Previous Article in Special Issue
Free-Space Nonreciprocal Transmission Based on Nonlinear Coupled Fano Metasurfaces
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Fundamentals of Lossless, Reciprocal Bianisotropic Metasurface Design

Department of Electrical Engineering and Computer Science, University of Michigan, Ann Arbor, MI 48109, USA
*
Author to whom correspondence should be addressed.
Current address: 1301 Beal Avenue, EECS Building, Room 3238, Ann Arbor, MI 48109, USA.
Photonics 2021, 8(6), 197; https://doi.org/10.3390/photonics8060197
Submission received: 30 April 2021 / Revised: 29 May 2021 / Accepted: 29 May 2021 / Published: 4 June 2021
(This article belongs to the Special Issue Advances in Complex Media Electromagnetics)

Abstract

:
Lossless, reciprocal bianisotropic metasurfaces have the ability to control the amplitude, phase, and polarization of electromagnetic wavefronts. However, producing the responses that are necessary for achieving this control with physically realizable surfaces is a challenging task. Here, several design approaches for bianisotropic metasurfaces are reviewed that produce physically realizable metasurfaces using cascaded impedance sheets. In practice, three or four impedance sheets are often used to realize bianisotropic responses, which can result in narrowband designs that require the unit cells to be optimized in order to improve the performance of the metasurface. The notion of a metasurface quality factor is introduced for three-sheet metasurfaces to address these issues in a systematic manner. It is shown that the quality factor can be used to predict the bandwidth of a homogeneous metasurface, and it can also be used to locate problematic unit cells when designing inhomogeneous metasurfaces. Several design examples are provided to demonstrate the utility of the quality factor, including an impedance matching layer with maximal bandwidth and a gradient metasurface for plane wave refraction. In addition to these examples, several metasurfaces for polarization control are also reported, including an isotropic polarization rotator and an asymmetric circular polarizer.

1. Introduction

Metasurfaces are the two-dimensional analogue of metamaterials, which interact with electromagnetic fields at a surface rather than throughout a volume [1]. They are often realized using electrically thin layers that consist of 2D arrays of subwavelength-spaced meta-atoms that can be homogenized. This allows for metasurfaces to be modeled using surface boundary conditions, called the generalized sheet transition conditions (GSTCs), which determine the interaction with an incident field through quasi-static surface polarizabilities [2,3,4]. If the incident fields are time-varying, then the surface polarizabilities produce equivalent electric and magnetic polarization currents. These equivalent currents can be related to surface admittances and impedances [5,6,7,8,9,10,11]. The mechanisms by which electric and magnetic surface currents are induced on a metasurface determine its classification as an electric, magnetic, electric and magnetic, or bianisotropic metasurface. Purely electric metasurfaces only contain electric polarizabilites, which interact with only the electric field to produce electric currents. Purely magnetic metasurfaces only contain magnetic polarizabilities, which only interact with the magnetic field to produce magnetic currents. Electric and magnetic surfaces contain both electric and magnetic polarizabilities in a single sheet. Finally, a metasurface that is bianisotropic can also contain electric and magnetic polarizabilities, as well as electro-magnetic and magneto-electric polarizabilities, i.e., magnetic polarization due to an electric field and electric polarization due to an magnetic field.
Bianisotropic metasurfaces provide the metasurface designer with the most degrees of freedom, which makes them useful for the extreme manipulation of electromagnetic fields. Bianisotropic metasurfaces include a wide range of both reciprocal and non-reciprocal responses. However, the focus here will be on design methods for reciprocal bianisotropic metasurfaces. Reciprocal bianisotropic metasurfaces can be split into two main classes: chiral and omega. Chiral metasurfaces contain meta-atoms that have broken mirror symmetry. This results in electric fields inducing magnetic currents along the impinging electric field and magnetic fields inducing electric currents along the impinging magnetic field. These chiral responses alter the polarization state of the incident wave. On the other hand, omega metasurfaces contain meta-atoms with broken directional symmetry. This results in electric fields inducing magnetic currents that are orthogonal to the impinging electric field, and magnetic fields inducing electric currents orthogonal to the impinging magnetic field. This leads to an asymmetric scattering response from omega metasurfaces, which maintains the polarization state.
The applications for chiral and omega bianisotropic metasurfaces fall into two main categories: those that guide and radiate electromagnetic waves and those that control reflection and transmission from a surface. Guided-wave bianisotropic metasurfaces shape fields along the surface through guided or leaky waves, and they can be used to produce desired radiation patterns [12,13]. Whereas, metasurfaces that control reflection and transmission interact with incident wavefronts to manipulate the amplitude, phase, and polarization of the scattered fields. There are many design synthesis methods and realizations of planar and cylindrical bianisotropic metasurfaces that control reflection and transmission at frequencies that range from microwave to optical using both composite (metal/dielectric) and all-dielectric metasurfaces [5,6,7,11,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30]. This is by no means a complete representation of all the work in bianisotropic metasurfaces. For a more complete review of the literature, see [31].
In this paper, the design and synthesis methods that are presented in [5,23] are reviewed, and several design examples are provided. Additionally, a definition for the quality factor of a three-sheet metasurface is provided, which can be used to estimate the bandwidth of a homogeneous metasurface. This is demonstrated through the design of an impedance matching metasurface with maximal bandwidth. In addition to improving bandwidth, the quality factor can also aid designers in improving the performance of inhomogeneous metasurfaces. This is demonstrated by using the quality factor to guide the selection of appropriate unit cells in the design of a gradient metasurface for plane wave refraction.

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 η 1 = μ 1 / ϵ 1 and of region 2 is η 2 = μ 2 / ϵ 2 . The metasurface is at the z = 0 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 2 × 2 matrix.
S nm = S n m x x S n m x y S n m y x S n m y y
When viewed from region 1, S 11 is the reflection coefficient and S 21 is the transmission coefficient. Similarly, when viewed from region 2, the reflection coefficient is S 22 and the transmission coefficient is S 12 .
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 2 × 2 tensors: the electric sheet admittance tensor Y , the magnetic sheet impedance tensor Z , 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.
J M = Y χ γ Z E avg H avg = z ^ × H 2 ¯ H 1 ¯ z ^ × E 2 ¯ E 1 ¯
The variables Y , χ , γ , and Z 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 E = [ E x E y ] T and the magnetic field vector is H = [ H x H y ] T (the surface current quantities J and M are similarly defined), where the averaged fields are E avg = ( E 1 + E 2 ) / 2 and H avg = ( H 1 + H 2 ) / 2 . The electric admittance tensor is defined as
Y = Y x x Y x y Y y x Y y y
with the other parameters being similarly defined.
For a reciprocal metasurface, Y = Y T , γ = χ T , and Z = Z T [32]. Imposing isotropy on the surface parameters results in
Y = Y I , χ = χ n , γ = γ n , Z = Z I
where I is the 2 × 2 identity matrix and n = 0 1 1 0 .
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
E = E 1 E 2 = S 11 S 12 S 21 S 22 E 1 + E 2 + = S E +
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
J x = Y 2 ( E 1 x + E 2 x ) + χ 2 ( H 1 y + H 2 y ) = H 2 y + H 1 y
M y = γ 2 ( E 1 x + E 2 x ) + Z 2 ( H 1 y + H 2 y ) = E 2 x + E 1 x
From Equations (6) and (7), we obtain four equations by considering the illumination from region 1 ( E 2 + = 0 ) and region 2 ( E 1 + = 0 ) separately. These four equations relate the S-parameters to the surface parameters of the metasurface. In each case, E 1 and E 2 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.
1 2 Y χ γ Z = 1 η 1 S 11 η 1 S 21 η 2 1 η 2 S 22 η 2 S 12 η 1 1 + S 11 S 21 1 S 22 + S 12 1 + S 11 + S 21 1 + S 22 + S 12 1 η 1 S 11 η 1 + S 21 η 2 1 η 2 + S 22 η 2 S 12 η 1 1
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
S 11 = 1 σ Y + Z 1 η 1 η 2 + 1 4 η 1 ( 2 γ ) ( 2 χ ) + Y Z 1 4 η 2 ( 2 + γ ) ( 2 + χ ) + Y Z
S 12 = 1 σ 1 2 η 2 ( 2 γ ) ( 2 + χ ) Y Z
S 21 = 1 σ 1 2 η 1 ( 2 + γ ) ( 2 χ ) Y Z
S 22 = 1 σ Y + Z 1 η 1 η 2 1 4 η 1 ( 2 γ ) ( 2 χ ) + Y Z 1 4 η 2 ( 2 + γ ) ( 2 + χ ) + Y Z
σ = Y + Z 1 η 1 η 2 + 1 4 η 1 ( 2 γ ) ( 2 χ ) + Y Z + 1 4 η 2 ( 2 + γ ) ( 2 + χ ) + Y Z
In the case of a lossless metasurface, the surface parameters Y = j B and Z = j X are purely imaginary, while γ = χ = R are real quantities [32]. In this case, Equations (9)–(12) can be further simplified, as shown in Equations (14)–(17). Note that S 21 = S 12 only when η 1 = η 2 , and S 11 = S 22 only when R = 0 and η 1 = η 2 .
S 11 = 1 σ j X η 1 η 2 B + 1 4 η 1 ( 2 R ) 2 B X 1 4 η 2 ( 2 + R ) 2 B X
S 12 = 1 σ 1 2 η 2 4 R 2 + B X
S 21 = 1 σ 1 2 η 1 4 R 2 + B X
S 22 = 1 σ j X η 1 η 2 B 1 4 η 1 ( 2 R ) 2 B X 1 4 η 2 ( 2 + R ) 2 B X
σ = j X η 1 η 2 + B + 1 4 η 1 ( 2 R ) 2 B X + 1 4 η 2 ( 2 + R ) 2 B X
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, S 21 = S 12 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 | S 21 | | S 12 | .
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
η 1 η 2 | S 21 | 2 = η 2 η 1 | S 12 | 2
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
η 1 η 2 S 21 = η 2 η 1 S 12 .
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
P a v g = 1 2 E 1 + E 2 + + E 1 E 2 T H 1 + H 2 + H 1 H 2 * = 0
By applying the plane wave relation between the electric and magnetic fields, and expressing E 1 and E 2 in terms of the S-parameters from Equation (5), Equation (21) becomes
P a v g = 1 2 R e E + + S E + T 1 / η E + 1 / η S E + * = 0
where
1 / η = 1 / η 1 0 0 1 / η 2 , E + = E 1 + E 2 + .
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.
1 = | S 11 | 2 + η 1 η 2 | S 21 | 2
1 = | S 22 | 2 + η 1 η 2 | S 21 | 2
0 = | S 11 | cos ϕ S 11 ϕ S 21 + ϕ E 1 + ϕ E 2 + + | S 22 | cos ϕ S 21 ϕ S 22 + ϕ E 1 + ϕ E 2 +
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 | S 21 | is chosen, as is commonly the case, then | S 11 | = | S 22 | , and a phase constraint is obtained from Equation (26)
ϕ S 11 ϕ S 21 = ϕ S 21 ϕ S 22 + π
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 Y = j B , Z = j X , and γ = χ = R . 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 ( ε r = 9.4 ), 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, Z i n . Because the two media have different wave impedances the metasurface must transform the wave impedance of the transmitted wave, Z L , to that of the incident wave, Z i n . 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 ( ε r = 9.4 ) 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 Z L = | Z L | e j ϕ L to a source impedance Z i = | Z i n | e j ϕ i n with an arbitrary transmission phase ϕ 21 [36]. Consider a general lossless two-port Z-matrix in order to determine the Z-matrix that provides the desired functionality,
V 1 V 2 = j X 11 X 12 X 21 X 22 I 1 I 2
Imposing the impedance boundary conditions and enforcing power conservation on (28) produces the following system of equations,
1 r v e j ϕ 21 = j X 11 X 12 X 21 X 22 1 r v e j ϕ 21
where r v 2 = Z L Z i n cos ϕ i n cos ϕ L and ϕ 21 = V 2 V 1 . Splitting (29) into its real and imaginary components allows for the elements of the Z-matrix to be solved in terms of Z L , Z i n and ϕ 21 ,
X 11 X 12 X 21 X 22 = | Z i n | cos ( ϕ 21 ϕ L ) | Z i n | r v cos ( ϕ L ) | Z i n | r v cos ( ϕ L ) | Z L | cos ( ϕ 21 + ϕ i n ) csc ( ϕ 21 + ϕ i n ϕ L )
From (30), it is clear that the required two-port network is reciprocal, since X 12 = X 21 , 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),
Z s 1 = j Z 0 sin ( β d ) cos ( β d ) + ( X 12 + X 22 det Z ) Z 0 sin ( β d )
Z s 2 = j Z 0 2 sin 2 ( β d ) X 12 det Z + X 12 Z 0 sin ( 2 β d )
Z s 3 = j Z 0 sin ( β d ) cos ( β d ) + ( X 12 + X 11 det Z ) Z 0 sin ( β d ) ,
where det Z is the determinant of the Z-matrix and β and Z 0 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,
Q = ω 0 2 W e P d ,
where ω 0 is the angular resonant frequency, W e is the average electric energy stored in the network at ω 0 , and P d 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
Q = ω 0 2 Z i n ( C s 1 + β d 2 ω 0 Z 0 ) + R i n t ( C s 2 + β d ω 0 Z 0 ) + Z L ( C s 3 + β d 2 ω 0 Z 0 ) ,
where, R i n t = Z i n + Z L + Z i n Z L cos ϕ 21 sin 2 ϕ 21 ( Z 0 sin β d ) 2 Z i n Z L , and C s i is the capacitance of the ith impedance sheet (if the sheet is inductive, then C s i = 0 ). If the load impedances are not purely real, then the imaginary part of the load can be absorbed into either Z s 1 or Z s 3 , and (35) can still be used. The quality factor, Q, of the metasurface will be used to approximate the fractional bandwidth, F B W = B W / f 0 , where B W 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 Z i n = 377 Ω . Let us assume that Z L = 123 Ω , and the spacers are free-space with a thickness d = λ 0 / 20 . 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 ϕ 21 = 68 . 5 . 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: Z s 1 = 1 / ( j ω C s 1 ) = j 468.9 Ω , Z s 2 = 1 / ( j ω C s 2 ) = j 641.9 Ω , and Z s 3 = j ω L s 3 = j38.5 k Ω . The metasurface performance is simulated in Ansys HFSS while using dispersive impedance sheets that correspond to the following lumped elements: C s 1 = 33.9 f F , C s 2 = 24.8 f F , and L s 3 = 612.7 n H . 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 ε r = 9.4 , 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 ε r , its quality factor is given by,
Q = ω 0 2 Z i n ( C s 1 + ε r ε 0 d 2 ) + R i n t ( C s 2 + ε r ε 0 d ) + Z L ( C s 3 + ε r ε 0 d 2 ) .
If the dielectric spacer is electrically thin, then the sheet capacitances can be approximated as,
C s 1 = 1 ω 1 ω μ 0 d + X 22 + X 12 det ( Z ) ε r ε 0 d 2
C s 2 = 1 ω 2 ω μ 0 d + det ( Z ) X 12 1 ( ω μ 0 d ) 2 ε r ε 0 d
C s 3 = 1 ω 1 ω μ 0 d + X 11 + X 12 det ( Z ) ε r ε 0 d 2 ,
when the sheet impedance Z s i 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, R i n t 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 C s i , ε r ε 0 d 2 , and ε r ε 0 d . When considering the terms C s 1 + ε r ε 0 d 2 , C s 2 + ε r ε 0 d , and C s 3 + ε r ε 0 d 2 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 ϕ 21 = 0 , 180 , and 360 , 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 70 at a frequency of f 0 = 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 ( k i = k sin ( θ i ) ) to produce the desired refracted wavenumber ( k t = k sin ( θ t ) ), where k is the wavenumber in the surrounding medium. Therefore, the metasurface must impart transverse momentum that is equal to Δ k = k t k i . Practically, this is realized by discretizing the metasurface into N sub-wavelength unit cells that are of size D = 2 π N max ( k i , k t ) , each possessing a transmission phase ϕ j , such that Δ ϕ = ϕ j + 1 ϕ j = Δ k D . 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 d = λ 0 / 25 . 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 ϕ 21 = 180 .
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 Q ( ϕ 21 ) 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 d < λ / 40 . 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 2 × 2 matrix of Equation (1)). Note that z ^ × E x E y T = n E x E y T . Expressing the surface parameters in terms of the S-parameters gives [5]
1 2 Y χ γ Z = I η 1 S 11 η 1 S 21 η 2 I η 2 S 22 η 2 S 12 η 1 n I + S 11 S 21 n I + S 22 S 12 I + S 11 + S 21 I + S 22 + S 12 n I η 1 S 11 η 1 + S 21 η 2 n I η 2 S 22 η 2 + S 12 η 1 1
Equation (40) can also be re-arranged so that the S-parameters are expressed in terms of surface parameters.
S 11 S 12 S 21 S 22 = I η 1 + Y 2 χ n 2 η 1 I η 2 + Y 2 + χ n 2 η 2 n + γ 2 Z n 2 η 1 n + γ 2 + Z n 2 η 2 1 I η 1 Y 2 χ n 2 η 1 I η 2 Y 2 + χ n 2 η 2 n γ 2 Z n 2 η 1 n γ 2 + Z n 2 η 2
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 4 × 4 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)
S 11 S 12 S 21 S 22 = S 11 S 12 S 21 S 22 T
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:
S 11 S 12 S 21 S 22 T S 11 S 12 S 21 S 22 * = 1 0 0 0 0 1 0 0 0 0 1 0 0 0 0 1 .
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 Y = Y T , γ = χ T , and Z = Z T for a reciprocal metasurface [32]. There are three free entries in both of the Y and Z matrices, and four free entries in the γ or χ matrix. Moreover, for the metasurface to be lossless, Y and Z 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 Y , 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 Y x x , Y y y and Y x y = Y y x . 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 M is defined as:
E L + E L = M E R + E R .
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, S spec , is stipulated based on the required application. It is then converted to a wave matrix, as follows [23]:
M spec = I 0 S 11 , spec S 12 , spec S 21 , spec S 22 , spec 0 I 1 ,
where 0 represents a 2 × 2 null matrix. This wave matrix M spec is known and it is set as the synthesis goal. It is worth noting that, if S 21 , spec has a zero determinant, taking the inverse matrix in (45) becomes invalid. In this case, a perturbation can be introduced into S 21 , spec 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 E 1 ± to E 4 ± is
M casc = M sheet ( 1 ) M dielectric ( 2 ) M sheet ( 2 ) M dielectric ( 3 ) M sheet ( 3 )
in which M sheet ( 1 ) , M sheet ( 2 ) and M sheet ( 3 ) are the admittance tensors that need to be solved. By setting M spec = M casc , and with some algebraic manipulation, one can find the admittance tensor of the second sheet Y 2 [23]:
e Y 2 = 1 A 2 ( ( e I ) M spec ( e I ) ( e t 1 Φ 2 t 2 Φ 3 t 3 e ) I ) .
The symbol ⊗ in (47) represents the Kronecker product of matrices and A 2 is some constant scalar. The matrix t 1 contains information regarding the dielectric interface where the first sheet is located, Φ 2 carries the phase information of the second dielectric spacer ( t 2 , t 3 , and Φ 3 are similarly defined), and e is a constant matrix [23]:
t 1 = 1 2 η 2 η 2 + η 1 η 2 η 1 η 2 η 1 η 2 + η 1 , Φ 2 = e j φ 2 0 0 e j φ 2 , e = 1 1 1 1 .
Similarly, the admittance tensors of the first and third sheets Y 1 and Y 3 can be expressed in terms of Y 2 ,
e Y 1 = 1 A 1 [ M spec ( e I ) ( t 1 Φ 2 t 2 Φ 3 t 3 e ) I η 2 2 ( t 1 Φ 2 e Φ 3 t 3 e ) Y 2 ] · ( I ( I + B 1 A 1 Y 2 ) 1 ) e Y 3 = 1 A 3 ( I ( I + B 3 A 3 Y 2 ) 1 ) · [ ( e I ) M spec ( e t 1 Φ 2 t 2 Φ 3 t 3 ) I η 2 2 ( e t 1 Φ 2 e Φ 3 t 3 ) Y 2 ]
where A 1 , B 1 , A 3 , and B 3 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.

5. Design Examples

In this section, several polarization-converting design examples are shown and discussed in order to illustrate the broad applicability of bianisotropic metasurfaces. For the details of the analysis and synthesis procedure, please refer to [23]. All of the devices considered here are realized by cascading several electric sheets and, thus, can be realized in practice using subwavelength patterned surfaces. The admittances of the electric sheets can be characterized analytically [37,38,44], or through full-wave extraction methods. However, evanescent coupling resulting from the fine features of the patterning may shift the unit cell’s response when the impedance sheets are cascaded to produce the bianisotropic unit cells. In order to account for this, the impedance sheets will need to be designed, such that the desired response of the unit cell is maintained. Several examples of metasurfaces demonstrating the feasibility of this approach at microwave and millimeter wave frequencies can be found in [5,7,18,45].

5.1. Asymmetric Circular Polarizer

An asymmetric circular polarizer is the first device presented here. The device converts a right-hand circularly polarized incident wave to a left-hand circularly polarized transmitted wave. On the contrary, when the incident wave is left-hand circularly polarized, it is totally reflected. Figure 12a provides an illustration of the operation [23]. Accordingly, the device has S-parameters:
S 12 , spec = S 21 , spec T = e j ϕ 2 1 j j 1 , S 11 , spec = S 22 , spec = e j ϕ 2 1 j j 1 .
for normal incidence.
Let us synthesize this asymmetric circular polarizer while using the wave matrix approach [23]. This specified S-matrix leads to a S 21 , spec with a vanishing determinant. Hence, a perturbation is required in this case, as noted in the previous section:
S 12 , spec = S 21 , spec T = e j ϕ 2 1 j j e j 1 .
This device can be realized by cascading three sheets. Figure 12b shows a unit cell of the device. In this example, we stipulate both dielectric spacers to have a dielectric constant ε 2 = ε 3 = 5 . It is assumed that the phase delay ϕ = 0 and the electrical lengths of both dielectric spacers are φ 2 = φ 3 = 2 π / 5 . By substituting these parameters into the design equations, (47) and (49), the following sheet admittance tensors are obtained,
Y 1 = j η 0 0.73 1.00 1.00 0.72 , Y 2 = j η 0 1268.31 5.52 5.52 1.43 , Y 3 = Y 1 .
Synthesized admittance values (52) are modeled in Ansys HFSS as anisotropic boundary conditions for full-wave verification of the design. At frequencies other than the design frequency, ω 0 , the sheet admittances are assumed to obey Foster’s reactance theorem. The eigenvalues of each sheet are first found by diagonalizing the tensors (52). We assume a capacitive frequency dependence if the resulting susceptances B 0 are positive,
B c ( ω ) = ω ω 0 B 0 , if B 0 > 0 .
On the other hand, negative susceptances are assumed to possess an inductive response
B l ( ω ) = ω 0 ω B 0 , if B 0 < 0 .
The assumptions, (53) and (54), are, in fact, reasonable since the admittances are realized by simple metallic patterns. Figure 12c displays the simulated frequency response of the asymmetric circular polarizer, where the transmission characteristics match the specified performance at the operating frequency (51).

5.2. Polarization Rotator

The second device considered in this section is a reflectionless polarization rotator that rotates the polarization of any linearly polarized incident wave by 90 , as shown in Figure 13a.
The stipulated S-parameters for this example are,
S 12 , spec = S 21 , spec T = e j ϕ 0 1 1 0 , S 11 , spec = S 22 , spec = e j ϕ 0 0 0 0 .
As discussed earlier, four sheets are required to reconstruct the stipulated S-matrix. With the unit cell that is shown in Figure 13b, it is assumed that ϕ = π / 4.5 , ε 2 = ε 3 = ε 4 = 3.5 , and φ 2 = φ 3 = φ 4 = π / 5 . Following the synthesis procedure that is outlined in [23], the admittance tensors of each sheet are calculated to be:
Y 1 = j η 0 5.01 0.77 0.77 0.13 , Y 2 = j η 0 9.30 0 0 1.00 Y 3 = j η 0 7.59 7.77 7.77 2.71 , Y 4 = j η 0 2.57 1.30 1.30 2.57 .
The metasurface unit cell was simulated in Ansys HFSS, and Figure 13c plots the resulting frequency response. Again, the transmission characteristics match the specified performance at the design frequency.

6. Conclusions

In this paper, two design procedures for realizing reciprocal bi-isotropic and bianisotropic metasurfaces using cascaded impedance sheets were reviewed. The design procedures use generalized sheet transition conditions (GSTCs) to relate bianisotropic surface parameters to the scattering parameters of the cascaded sheet impedances. Such approaches allow for metasurfaces with arbitrary lossless, reciprocal, and bianisotropic responses to be realized in practice. These design methods were then used to realize several examples of practical devices with phase and polarization control. Specifically, they were used to realize an asymmetric circular polarizer and a reflectionless, polarization rotator.
In addition to these design procedures, the quality factor for metasurfaces that are composed of three impedance sheets was introduced. The quality factor was shown to predict the bandwidth of a homogeneous metasurface. This was demonstrated through the design of an impedance matching metasurface with maximal bandwidth. It was also shown that the quality factor could be used to improve the performance of inhomogeneous metasurfaces. This was demonstrated through the design of a gradient metasurface for plane wave refraction. The unit cell quality factor was used to identify cells that degraded the overall metasurface performance, and it was used to select alternative unit cells that improved the overall performance. Such an approach can be used to balance manufacturability and performance trade-offs for metasurface devices.

Author Contributions

Conceptualization, A.G.; methodology, L.S., C.-W.L. and B.O.R.; software, L.S., C.-W.L. and B.O.R.; validation, L.S., C.-W.L. and B.O.R.; investigation, L.S., C.-W.L. and B.O.R.; resources, A.G.; writing—original draft preparation, L.S., C.-W.L. and B.O.R.; writing—review and editing, A.G., L.S., C.-W.L. and B.O.R.; visualization, L.S., C.-W.L. and B.O.R.; supervision, A.G.; project administration, A.G.; funding acquisition, A.G. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Office of Naval Research under Grant No. N00014-18-1-2536, the National Science Foundation under the Grant Opportunities for Academic Liaison with Industry (GOALI) program under Grant 1807940, the National Science Foundation Graduate Research Fellowship Program under Grant No. DGE 1256260, and the Chia-Lun Lo Fellowship at the University of Michigan.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Here, the approximations that are used to calculate the quality factor of a metasurface with electrically thin spacers are derived. First, we will show that the internal matching resistance, R i n t , is independent of the spacer’s dielectric constant, ε r . For an electrically thin spacer β d < < 1 . Using second-order small angle approximations, sin β d β d and cos β d 1 ( β d ) 2 2 , in the expression for R i n t yields,
R i n t = Z i n + Z L + Z i n Z L cos ϕ 21 sin 2 ϕ 21 ( ω μ 0 d ) 2 Z i n Z L .
Which, does not depend on the dielectric constant of the spacer so R i n t is independent of the spacer’s dielectric constant.
Next, we will derive the approximate expressions for the impedance sheets (31)–(33). Converting the sin ( 2 β d ) term in (32) to 2 sin ( β d ) cos ( β d ) and making the same small argument approximations in (31)–(33),
Z s 1 = j ω μ 0 d 1 ω 2 μ 0 ε r ε 0 d 2 2 + ( X 12 + X 22 det Z ) ω μ 0 d
Z s 2 = j ( ω μ 0 d ) 2 X 12 det Z + 2 X 12 ω μ 0 d ( 1 ω 2 μ 0 ε r ε 0 d 2 2 )
Z s 3 = j ω μ 0 d 1 ω 2 μ 0 ε r ε 0 d 2 2 + ( X 12 + X 11 det Z ) ω μ 0 d .
Because we only need to consider capactive sheets in the calculation for the quality factor, we can set (A2)–(A4) equal to the impedance of a capacitive element Z s i = 1 / j ω C s i . The following expressions are produced by setting (A2)–(A4) equal to 1 / j ω C s i and solving for the sheet capacitances,
C s 1 = 1 ω 1 ω μ 0 d + X 22 + X 12 det ( Z ) ε r ε 0 d 2
C s 2 = 1 ω 2 ω μ 0 d + det ( Z ) X 12 1 ( ω μ 0 d ) 2 ε r ε 0 d
C s 3 = 1 ω 1 ω μ 0 d + X 11 + X 12 det ( Z ) ε r ε 0 d 2 .

References

  1. Holloway, C.L.; Kuester, E.F.; Gordon, J.A.; O’Hara, J.; Booth, J.; Smith, D.R. An Overview of the Theory and Applications of Metasurfaces: Two-dimensional Equivalents of Metamaterials. IEEE Antennas Propag. Mag. 2012, 54, 10–35. [Google Scholar] [CrossRef]
  2. Idemen, M. Straightforward derivation of Boundary Conditions on Sheet Simulating Anisotropic Thin Layer. Electron. Lett. 1988, 24, 663–665. [Google Scholar] [CrossRef]
  3. Idemen, M. Universal Boundary Relations of the Electromagnetic Field. J. Phys. Soc. Jpn. 1990, 59, 71–80. [Google Scholar] [CrossRef]
  4. Kuester, E.F.; Mohamed, M.A.; Piket-May, M.; Holloway, C.L. Averaged Transition Conditions for Electromagnetic Fields at a Metafilm. IEEE Trans. Antennas Propag. 2003, 51, 2641–2651. [Google Scholar] [CrossRef]
  5. Pfeiffer, C.; Grbic, A. Bianisotropic Metasurfaces for Optimal Polarization Control: Analysis and Synthesis. Phys. Rev. Appl. 2014, 2, 044011. [Google Scholar] [CrossRef]
  6. Epstein, A.; Eleftheriades, G.V. Arbitrary Power-Conserving Field Transformations with Passive Lossless Omega-Type Bianisotropic Metasurfaces. IEEE Trans. Antennas Propag. 2016, 64, 3880–3895. [Google Scholar] [CrossRef] [Green Version]
  7. Chen, M.; Abdo-Sánchez, E.; Epstein, A.; Eleftheriades, G.V. Theory, design, and experimental verification of a reflectionless bianisotropic Huygens’ metasurface for wide-angle refraction. Phys. Rev. B 2018, 97, 125433. [Google Scholar] [CrossRef] [Green Version]
  8. Pfeiffer, C.; Grbic, A. Metamaterial Huygens’ Surfaces: Tailoring Wave Fronts with Reflectionless Sheets. Phys. Rev. Lett. 2013, 110, 197401. [Google Scholar] [CrossRef]
  9. Selvanayagam, M.; Eleftheriades, G.V. Discontinuous electromagnetic fields using orthogonal electric and magnetic currents for wavefront Manipulation. Opt. Express 2013, 21, 14409–14429. [Google Scholar] [CrossRef]
  10. Epstein, A.; Eleftheriades, G.V. Huygens’ metasurfaces via the equivalence principle: Design and applications. J. Opt. Soc. Am. B 2016, 33, A310–A50. [Google Scholar] [CrossRef]
  11. Achouri, K.; Salem, M.A.; Caloz, C. General Metasurface Synthesis Based on Susceptibility Tensors. IEEE Trans. Antennas Propag. 2015, 63, 2977–2991. [Google Scholar] [CrossRef] [Green Version]
  12. Tierney, B.B.; Grbic, A. Controlling Leaky waves with 1-D cascaded metasurfaces. IEEE Trans. Antennas Propag. 2018, 66, 2143–2146. [Google Scholar] [CrossRef]
  13. Abdo-Sánchez, E.; Chen, M.; Epstein, A.; Eleftheriades, G.V. A Leaky-Wave Antenna With Controlled Radiation Using a Bianisotropic Huygens’ Metasurface. IEEE Trans. Antennas Propag. 2019, 67, 108–120. [Google Scholar] [CrossRef] [Green Version]
  14. Kwon, D.H.; Werner, P.L.; Werner, D.H. Optical planar chiral metamaterial designs for strong circular dichroism and polarization rotation. Opt. Express 2008, 16, 11802–11807. [Google Scholar] [CrossRef] [PubMed]
  15. Kwon, D.H.; Werner, D.H.; Kildishev, A.V.; Shalaev, V.M. Material parameter retrieval procedure for general bi-isotropic metamaterials and its application to optical chiral negative-index metamaterial design. Opt. Express 2008, 16, 11822–11829. [Google Scholar] [CrossRef] [Green Version]
  16. Shaltout, A.; Shalaev, V.; Kildishev, A. Homogenization of bi-anisotropic metasurfaces. Opt. Express 2013, 21, 21941–21950. [Google Scholar] [CrossRef]
  17. Niemi, T.; Karilainen, A.O.; Tretyakov, S.A. Synthesis of Polarization Transformers. IEEE Trans. Antennas Propag. 2013, 614, 3102–3111. [Google Scholar] [CrossRef]
  18. Pfeiffer, C.; Zhang, C.; Ray, V.; Guo, L.J.; Grbic, A. High Performance Bianisotropic Metasurfaces: Asymmetric Transmission of Light. Phys. Rev. Lett. 2014, 113, 023902. [Google Scholar] [CrossRef]
  19. Wu, C.; Arju, N.; Kelp, G.; Fan, J.A.; Dominguez, J.; Gonzales, E.; Tutuc, E.; Brener, I.; Shvets, G. Spectrally selective chiral silicon metasurfaces based on infrared Fano resonances. Nat. Commun. 2014, 5, 3892. [Google Scholar] [CrossRef]
  20. Tretyakov, S.A. Metasurfaces for general transformations of electromagnetic fields. Philos. Trans. R. Soc. A 2015, 373, 20140362. [Google Scholar] [CrossRef]
  21. Alaee, R.; Albooyeh, M.; Rahimzadegan, A.; Mirmoosa, M.S.; Kivshar, Y.S.; Rockstuhl, C. All-dielectric reciprocal bianisotropic nanoparticles. Phys. Rev. B 2015, 92, 245130. [Google Scholar] [CrossRef] [Green Version]
  22. Odit, M.; Kapitanova, P.; Belov, P.; Alaee, R.; Rockstuhl, C.; Kivshar, Y.S. Experimental realisation of all-dielectric bianisotropic metasurfaces. Appl. Phys. Lett. 2016, 108, 221903. [Google Scholar] [CrossRef]
  23. Ranjbar, A.; Grbic, A. Analysis and synthesis of cascaded metasurfaces using wave matrices. Phys. Rev. B 2017, 95, 205114. [Google Scholar] [CrossRef] [Green Version]
  24. Kwon, D.H.; Ptitcyn, G.; Díaz-Rubio, A.; Tretyakov, S.A. Transmission Magnitude and Phase Control for Polarization-Preserving Reflectionless Metasurfaces. Phys. Rev. Appl. 2018, 9, 034005. [Google Scholar] [CrossRef] [Green Version]
  25. Ranjbar, A.; Grbic, A. Broadband, Multiband, and Multifunctional All-Dielectric Metasurfaces. Phys. Rev. Appl. 2019, 11, 054066. [Google Scholar] [CrossRef]
  26. Sipus, Z.; Bosiljevac, M.; Grbic, A. Modelling cascaded cylindrical metasurfaces using sheet impedances and a transmission matrix formulation. IET Microw. Antennas Propag. 2018, 12, 1041–1047. [Google Scholar] [CrossRef]
  27. Sipus, Z.; Eres, Z.; Barbaric, D. Modeling Cascaded Cylindrical Metasurfaces with Spatially-Varying Impedance Distribution. Radioengineering 2019, 28, 505–511. [Google Scholar] [CrossRef]
  28. Li, J.; Diaz-Rubio, A.; Jia, Z.; Shen, C.; Tretyakov, S.; Cummer, S. Highly Efficient Generation of Angular Momentum with Cylindrical Bianisotropic Metasurfaces. Phys. Rev. Appl. 2019, 11, 024016. [Google Scholar] [CrossRef] [Green Version]
  29. Xu, G.; Eleftheriades, G.V.; Hum, S.V. Discrete-Fourier-Transform-Based Framework for Analysis and Synthesis of Cylindrical Omega-Bianisotropic Metasurfaces. Phys. Rev. Appl. 2020, 14, 064055. [Google Scholar] [CrossRef]
  30. Lin, C.W.; Grbic, A. Analysis and Synthesis of Cascaded Cylindrical Metasurfaces using a Wave Matrix Approach. IEEE Trans. Antennas Propag. 2021. [Google Scholar] [CrossRef]
  31. Asadchy, V.S.; Díaz-Rubio, A.; Tretyakov, S.A. Bianisotropic metasurfaces: Physics and applications. Nanophotonics 2018, 7, 1069–1094. [Google Scholar] [CrossRef]
  32. Kong, J.A. Theorems of bianisotropic media. Proc. IEEE 1972, 60, 1036. [Google Scholar] [CrossRef]
  33. Pozar, D.M. Microwave Engineering, 4th ed.; Wiley: New York, NY, USA, 2011; pp. 191–194. [Google Scholar]
  34. Vehmas, J.; Hrabar, S.; Tretyakov, S. Omega transmission lines with applications to effective medium models of metamaterials. Phys. Rev. Appl. 2014, 115, 134905. [Google Scholar] [CrossRef] [Green Version]
  35. Dorrah, A.H.; Chen, M.; Eleftheriades, G.V. Bianisotropic Huygens’ Metasurface for Wideband Impedance Matching Between Two Dielectric Media. IEEE Trans. Antennas Propag. 2018, 66, 4729–4742. [Google Scholar] [CrossRef]
  36. Wong, J.P.; Epstein, A.; Eleftheriades, G.V. Reflectionless Wide-Angle Refracting Metasurfaces. IEEE Antennas Wirel. Propag. Lett. 2016, 15, 1293–1296. [Google Scholar] [CrossRef]
  37. Luukkonen, O.; Simovski, C.; Granet, G.; Goussetis, G.; Lioubtchenko, D.; Räisänen, A.V.; Tretyakov, S. Simple and Accurate Analytical Model of Planar Grids and High-Impedance Surfaces Comprising Metal Strips or Patches. IEEE Trans. Antennas Propag. 2008, 56, 1624–1632. [Google Scholar] [CrossRef] [Green Version]
  38. Langley, R.J.; Parker, E.A. Equivalent circuit model for arrays of square loops. Electron. Lett. 1982, 18, 294–296. [Google Scholar] [CrossRef]
  39. Yu, N.; Genevet, P.; Kats, M.A.; Aieta, F.; Tetienne, J.P.; Capasso, F.; Gaburro, Z. Light Propagation with Phase Discontinuities: Generalized Laws of Reflection and Refraction. Science 2011, 334, 333–337. [Google Scholar] [CrossRef] [Green Version]
  40. Raeker, B.O.; Grbic, A. Compound Metaoptics for Amplitude and Phase Control of Wave Fronts. Phys. Rev. Lett. 2019, 122, 113901. [Google Scholar] [CrossRef] [Green Version]
  41. Dorrah, A.H.; Eleftheriades, G.V. Bianisotropic Huygens’ Metasurface Pairs for Nonlocal Power-Conserving Wave Transformations. IEEE Antennas Wirel. Propag. Lett. 2018, 17, 1788–1792. [Google Scholar] [CrossRef]
  42. Popov, V.; Boust, F.; Tretyakov, S.; Burokur, S.N. Omega-bianisotropic metasurface for converting a propagating wave into a surface wave. Phys. Rev. B 2019, 100, 125103. [Google Scholar] [CrossRef] [Green Version]
  43. Budhu, J.; Grbic, A. Perfectly Reflecting Metasurface Reflectarrays: Mutual Coupling Modeling Between Unique Elements Through Homogenization. IEEE Trans. Antennas Propag. 2021, 69, 122–134. [Google Scholar] [CrossRef]
  44. Chiotellis, N.; Grbic, A. Analytical modeling of tensor metasurfaces. J. Opt. Soc. Am. B 2021, 33, A51–A60. [Google Scholar] [CrossRef]
  45. Yang, F.; Raeker, B.O.; Nguyen, D.T.; Miller, J.D.; Xiong, Z.; Grbic, A.; Ho, J.S. Antireflection and Wavefront Manipulation with Cascaded Metasurfaces. Phys. Rev. Appl. 2020, 14, 064044. [Google Scholar] [CrossRef]
Figure 1. The geometry of a metasurface between two regions with different material properties. The equivalent surface current densities J (electric) and M (magnetic) describe the interaction of the metasurface with the tangential fields. Under illumination by a normally incident plane wave, each region can contain two plane waves that are denoted by + for a wave propagating toward the surface or − for a wave propagating away from the surface.
Figure 1. The geometry of a metasurface between two regions with different material properties. The equivalent surface current densities J (electric) and M (magnetic) describe the interaction of the metasurface with the tangential fields. Under illumination by a normally incident plane wave, each region can contain two plane waves that are denoted by + for a wave propagating toward the surface or − for a wave propagating away from the surface.
Photonics 08 00197 g001
Figure 2. A metasurface at the interface between air and alumina half-spaces. The metasurface is used to impedance match a normally incident plane wave traveling from the region of air into the alumina.
Figure 2. A metasurface at the interface between air and alumina half-spaces. The metasurface is used to impedance match a normally incident plane wave traveling from the region of air into the alumina.
Photonics 08 00197 g002
Figure 3. L and T circuit network topologies used for impedance matching in circuit theory.
Figure 3. L and T circuit network topologies used for impedance matching in circuit theory.
Photonics 08 00197 g003
Figure 4. Bi-isotropic metasurface realized using three impedance sheets that are separated by dielectric spacers with thickness d.
Figure 4. Bi-isotropic metasurface realized using three impedance sheets that are separated by dielectric spacers with thickness d.
Photonics 08 00197 g004
Figure 5. The quality factor, fractional bandwidth, and the magnitude of the frequency response for metasurfaces that provide impedance matching with six different transmission phases. (a) The quality factor is minimized when the transmission phase is 68 . 5 (red circle). (b) The fractional bandwidth is maximized at 68 . 5 (red circle). (c) Plots of the transmission amplitude over frequency for several transmission phases and the maximum bandwidth is observed when the transmission phase is 68 . 5 , as predicted by the quality factor.
Figure 5. The quality factor, fractional bandwidth, and the magnitude of the frequency response for metasurfaces that provide impedance matching with six different transmission phases. (a) The quality factor is minimized when the transmission phase is 68 . 5 (red circle). (b) The fractional bandwidth is maximized at 68 . 5 (red circle). (c) Plots of the transmission amplitude over frequency for several transmission phases and the maximum bandwidth is observed when the transmission phase is 68 . 5 , as predicted by the quality factor.
Photonics 08 00197 g005
Figure 6. Plots of the (a) transmission and (b) reflection magnitudes for the interface with the metasurface ( ϕ 21 = 68 . 5 ), a quarter-wave transformer, and with no impedance matching (bare interface). The simulations of the metasurface were performed in Ansys HFSS. The metasurface has a bandwidth that is comparable to a quarter-wave transformer.
Figure 6. Plots of the (a) transmission and (b) reflection magnitudes for the interface with the metasurface ( ϕ 21 = 68 . 5 ), a quarter-wave transformer, and with no impedance matching (bare interface). The simulations of the metasurface were performed in Ansys HFSS. The metasurface has a bandwidth that is comparable to a quarter-wave transformer.
Photonics 08 00197 g006
Figure 7. (a) The depiction of an inhomogeneous, bi-isotropic metasurface implemented as a three-sheet cascade in free-space. (b) Plots of the sheet reactances for different transmission phases. The solid circles indicate the values used for the linear phase gradient and the empty squares indicate the sheet values used for the perturbed phase gradient. (c) Full-wave simulation results for the real part of the electric field using the metasurface with a linear phase gradient. (d) Full-wave simulation results for the real part of the electric field using the metasurface with a perturbed phase gradient.
Figure 7. (a) The depiction of an inhomogeneous, bi-isotropic metasurface implemented as a three-sheet cascade in free-space. (b) Plots of the sheet reactances for different transmission phases. The solid circles indicate the values used for the linear phase gradient and the empty squares indicate the sheet values used for the perturbed phase gradient. (c) Full-wave simulation results for the real part of the electric field using the metasurface with a linear phase gradient. (d) Full-wave simulation results for the real part of the electric field using the metasurface with a perturbed phase gradient.
Photonics 08 00197 g007
Figure 8. (a) The quality factor of the metasurface unit cells versus transmission phase. (b) The first derivative of the quality factor with respect to the transmission phase. The solid black circles indicate the values corresponding to the linear phase gradient and the hollow red squares indicate the adjusted values used for the perturbed phase gradient.
Figure 8. (a) The quality factor of the metasurface unit cells versus transmission phase. (b) The first derivative of the quality factor with respect to the transmission phase. The solid black circles indicate the values corresponding to the linear phase gradient and the hollow red squares indicate the adjusted values used for the perturbed phase gradient.
Photonics 08 00197 g008
Figure 9. A comparison of the transmission phases used for the original and perturbed phase gradients. The solid black circles indicate the transmission phases used for the linear phase gradient. The empty red squares indicate the transmission phases that are used for the perturbed phase gradient.
Figure 9. A comparison of the transmission phases used for the original and perturbed phase gradients. The solid black circles indicate the transmission phases used for the linear phase gradient. The empty red squares indicate the transmission phases that are used for the perturbed phase gradient.
Photonics 08 00197 g009
Figure 10. An illustration of the wave matrix and the constitutive blocks of the cascaded structure. (a) The definition of a wave matrix. (b) A metasurface interface between two dielectric media. (c) A dielectric spacer.
Figure 10. An illustration of the wave matrix and the constitutive blocks of the cascaded structure. (a) The definition of a wave matrix. (b) A metasurface interface between two dielectric media. (c) A dielectric spacer.
Photonics 08 00197 g010
Figure 11. A cascaded metasurface structure consists of three sheets with only electric responses ( γ = χ = Z = 0 ).
Figure 11. A cascaded metasurface structure consists of three sheets with only electric responses ( γ = χ = Z = 0 ).
Photonics 08 00197 g011
Figure 12. An asymmetric circular polarizer. (a) Polarization-converting operation of the metasurface [23]. (b) Unit cell of the asymmetric circular polarizer. (c) Simulated transmission coefficients.
Figure 12. An asymmetric circular polarizer. (a) Polarization-converting operation of the metasurface [23]. (b) Unit cell of the asymmetric circular polarizer. (c) Simulated transmission coefficients.
Photonics 08 00197 g012
Figure 13. A linear polarization rotator. (a) Polarization-converting operation of the metasurface [23]. (b) Unit cell of the polarization rotator. (c) Simulated transmission coefficients.
Figure 13. A linear polarization rotator. (a) Polarization-converting operation of the metasurface [23]. (b) Unit cell of the polarization rotator. (c) Simulated transmission coefficients.
Photonics 08 00197 g013
Table 1. The unit cell transmission phases ( ϕ 21 ) used in the design of the gradient metasurface for plane wave refraction. The original phase gradient corresponds to the linear phase gradient. The perturbed phase gradient corresponds to the adjusted phases used to improve the performance of the metasurface.
Table 1. The unit cell transmission phases ( ϕ 21 ) used in the design of the gradient metasurface for plane wave refraction. The original phase gradient corresponds to the linear phase gradient. The perturbed phase gradient corresponds to the adjusted phases used to improve the performance of the metasurface.
Unit Cell ϕ 21 (Original) ϕ 21 (Perturbed)
1 18 31
2 54 54
3 90 90
4 126 126
5 162 147
6 198 216
7 234 234
8 270 270
9 306 306
10 342 330
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Szymanski, L.; Raeker, B.O.; Lin, C.-W.; Grbic, A. Fundamentals of Lossless, Reciprocal Bianisotropic Metasurface Design. Photonics 2021, 8, 197. https://doi.org/10.3390/photonics8060197

AMA Style

Szymanski L, Raeker BO, Lin C-W, Grbic A. Fundamentals of Lossless, Reciprocal Bianisotropic Metasurface Design. Photonics. 2021; 8(6):197. https://doi.org/10.3390/photonics8060197

Chicago/Turabian Style

Szymanski, Luke, Brian O. Raeker, Chun-Wen Lin, and Anthony Grbic. 2021. "Fundamentals of Lossless, Reciprocal Bianisotropic Metasurface Design" Photonics 8, no. 6: 197. https://doi.org/10.3390/photonics8060197

APA Style

Szymanski, L., Raeker, B. O., Lin, C. -W., & Grbic, A. (2021). Fundamentals of Lossless, Reciprocal Bianisotropic Metasurface Design. Photonics, 8(6), 197. https://doi.org/10.3390/photonics8060197

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop