Frequency selective surface

See also: Photonic crystal and Metamaterial

A Frequency Selective Surface (FSS) is a doubly-periodic, planar structure designed to pass or reflect electromagnetic energy based on its frequency. In this sense, an FSS is a type of Optical filter in which the filtering is accomplished by virtue of the regular, periodic arrangement of planar scatterers on the surface of the FSS. Frequency selective surfaces have been most commonly used in the radio frequency region of the electromagnetic spectrum and find use in applications as diverse as microwave ovens, Antenna radomes, Metamaterials, and Stealth technology. Sometimes frequency selective surfaces are referred to periodic surfaces and are a 2-dimensional analog of the new periodic volumes known as Photonic crystals.

Field equations for 3D PEC photonic crystal structures

For perfectly electrically conducting (PEC) structures admitting only electric current sources J, the electric field E is related to the vector magnetic potential A via the well-known relation:

${\mathbf {E}}~=~-jk\eta \left[{\mathbf {A}}+{\frac {1}{k^{2}}}\nabla (\nabla \bullet {\mathbf {A}})\right]~~~~~~~~~~~~~~(1.1)$

and the vector magnetic potential is in turn related to the source currents via:

$\nabla ^{2}{\mathbf {A}}+k^{2}{\mathbf {A}}~=~-{\mathbf {J}}~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~(1.2)$

where

$k^{2}~=~\omega ^{2}(\mu \epsilon )~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~(1.3)$

Bloch wave expansion of the fields

To solve equations (1.1) and (1.2) within the periodic volume, we may assume a Bloch wave expansion for all currents, fields and potentials:

${\mathbf {J}}(x,y,z)~=~\sum _{mnp}~{\mathbf {J}}(\alpha _{m},\beta _{n},\gamma _{p})~e^{j(\alpha _{m}x+\beta _{n}y+\gamma _{p}z)}~~~~~(2.1a)$

${\mathbf {E}}(x,y,z)~=~\sum _{mnp}~{\mathbf {E}}(\alpha _{m},\beta _{n},\gamma _{p})~e^{j(\alpha _{m}x+\beta _{n}y+\gamma _{p}z)}~~~~(2.1b)$

${\mathbf {A}}(x,y,z)~=~\sum _{mnp}~{\mathbf {A}}(\alpha _{m},\beta _{n},\gamma _{p})~e^{j(\alpha _{m}x+\beta _{n}y+\gamma _{p}z)}~~~(2.1c)$

where for simplicity, we assume a cubic lattice in which α only depends on m, β only depends on n and γ only depends on p. In the equations above,

$\alpha _{m}~=~k_{0}~\sin \theta _{0}~\cos \phi _{0}~+~{\frac {2m\pi }{l_{x}}}~~~~~~~~~~~(2.2a)$

$\beta _{n}~=~k_{0}~\sin \theta _{0}~\sin \phi _{0}~+~{\frac {2n\pi }{l_{y}}}~~~~~~~~~~~~~(2.2b)$

$\gamma _{p}~=~k_{0}~\cos \theta _{0}~+~{\frac {2p\pi }{l_{z}}}~~~~~~~~~~~~~~~~~~~~~~(2.2c)$

and,

$k_{0}~=~2\pi /\lambda ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~(2.3)$

where l_x, l_y, l_z are the dimensions of the unit cell in the x,y,z directions respectively, λ is the effective wavelength in the crystal and θ₀, φ₀ are the directions of propagation in spherical coordinates. Note that k in equations (1.1) and (1.2) comes from the time derivative in Maxwell's equations and is the free space propagation constant, proportional to frequency as we see in equation (1.3). On the other hand, k₀ in the equations above comes from our assumed Bloch wave solution given by equations (2.1) & (2.2). As a result, it represents the propagation constant in the periodic medium. These two k's, i.e. the free space propagation constant and the propagation constant of the Bloch wave, are in general different thereby allowing for dispersion in our solution.

Integral equation for PEC media

Substituting equations (2.1) into (1.1) and (1.2) yields the spectral domain Greens function relating the radiated electric field to its source currents:

${\mathbf {E}}(\alpha _{m},\beta _{n},\gamma _{p})~=~{\frac {jk\eta }{k^{2}-\alpha _{m}^{2}-\beta _{n}^{2}-\gamma _{p}^{2}}}~{\mathbf {G}}_{mnp}~{\mathbf {J}}(\alpha _{m},\beta _{n},\gamma _{p})~~~~~~~~~~~~~~~~~~~~~~~~~(3.1)$

where,

${\mathbf {G}}_{mnp}~=~\left({\begin{matrix}1-{\frac {\alpha _{m}^{2}}{k^{2}}}&-{\frac {\alpha _{m}\beta _{n}}{k^{2}}}&-{\frac {\alpha _{m}\gamma _{p}}{k^{2}}}\\-{\frac {\alpha _{m}\beta _{n}}{k^{2}}}&1-{\frac {\beta _{n}^{2}}{k^{2}}}&-{\frac {\beta _{n}\gamma _{p}}{k^{2}}}\\-{\frac {\alpha _{m}\gamma _{p}}{k^{2}}}&-{\frac {\beta _{n}\gamma _{p}}{k^{2}}}&1-{\frac {\gamma _{p}^{2}}{k^{2}}}\end{matrix}}\right)~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~(3.2)$

With this, the electric field boundary condition on the surface of PEC material within a unit cell becomes:

$~\sum _{mnp}~{\frac {1}{k^{2}-\alpha _{m}^{2}-\beta _{n}^{2}-\gamma _{p}^{2}}}~{\mathbf {G}}_{mnp}~{\mathbf {J}}(\alpha _{m},\beta _{n},\gamma _{p})~e^{j(\alpha _{m}x+\beta _{n}y+\gamma _{p}z)}~=~{\mathbf {0}}~~~~~~~~~~~~~~~~(3.3)$

Since we are seeking characteristic modes (eigenmodes) of the structure, there is no impressed E-field on the RHS of this electric field integral equation (EFIE). Equation (3.3) is not strictly correct, since only the tangential electric field is zero on the surface of the PEC scatterer. This inexactness will be resolved presently when we test with the current basis functions, defined as residing on the surface of the scatterer.

Method of Moments solution

As is usual in the method of moments, we assume an expansion for the source currents over some known set of basis functions with unknown weighting coefficients J_j:

$~{\mathbf {J}}(x,y,z)~=~\sum _{j}~J_{j}~{\mathbf {J}}_{j}(x,y,z)~~~~~~~~~~~~~~~~~(3.4)$

Substituting (3.4) into (3.3) and then testing the resulting equation with the i-th current basis function (i.e., dotting from the left and integrating over the domain of the i-th current basis function, thereby completing the quadratic form) produces the i-th row of the matrix eigenvalue problem as:

$~\sum _{j}~J_{j}~\left[~\sum _{mnp}~{\frac {{\mathbf {J}}_{i}(-\alpha _{m},-\beta _{n},-\gamma _{p})~{\mathbf {G}}_{mnp}~{\mathbf {J}}_{j}(\alpha _{m},\beta _{n},\gamma _{p})}{k^{2}-\alpha _{m}^{2}-\beta _{n}^{2}-\gamma _{p}^{2}}}\right]~=~{\mathbf {0}}~~~(3.5)$

This matrix equation is very simple to implement and requires only that the 3D FT of the basis functions be computed, preferably in closed form. With this method, computing bands of a 3D photonic crystal is as easy as computing reflection and transmission from a 2D periodic surface. In fact, equation (3.5) is identical to the basic EFIE for PEC FSS (Scott [1989]), the only difference being the stronger singularity in 3D which accelerates convergence of the triple sums.

Computing bands

To compute bands of the crystal (i.e. k-k₀ diagrams), we may assume values for (k₀, θ₀, φ₀) and then search for those values of k which drive the determinant of the impedance matrix to zero. Equation (3.5) has been used to compute bands in various types of doped and undoped photonic crystals (Scott[1998], Scott [2002]).

See also

• Bloch wave
• Photonic crystal
• Metamaterial
• Fourier optics
• Bragg diffraction

References

• Kastner, Raphael (1987), On the Singularity of the Full Spectral Green's Dyad, IEEE Trans. on Antennas and Propagation, vol. AP-35, No. 11, pp. 1303–1305
• Scott, Craig (1989), The Spectral Domain Method in Electromagnetics, Artech House, ISBN 0-89006-349-4
• Scott, Craig (1998), Analysis, Design and Testing of Integrated Structural Radomes Built Using Photonic Bandgap Structures, 1998 IEEE Aerospace Conf. Aspen CO, pp. 463 - 479
• Scott, Craig (2002), Spectral Domain Analysis of Doped Electromagnetic Crystal Radomes Using the Method of Moments, 2002 IEEE Aerospace Conf. Big Sky MT, paper #504, pp. 2-957 - 2-963