Bloch wave – MoM method

See also: Photonic crystal and Metamaterial

Electromagnetism

Electricity Magnetism
Electrostatics Electric charge Static electricity Electric field Conductor Insulator Triboelectricity Electrostatic discharge Induction Coulomb's law Gauss's law Electric flux / potential energy Electric dipole moment Polarization density
Magnetostatics Ampère's law Magnetic field Magnetization Magnetic flux Biot–Savart law Magnetic dipole moment Gauss's law for magnetism
Electrodynamics Lorentz force law Electromagnetic induction Faraday's law Lenz's law Displacement current Magnetic potential Maxwell's equations Electromagnetic field Electromagnetic radiation Maxwell tensor Poynting vector Liénard–Wiechert potential Jefimenko's equations Eddy current London equations Mathematical descriptions of the electromagnetic field
Electrical network Electric current Electric potential Voltage Resistance Ohm's law Series circuit Parallel circuit Direct current Alternating current Charged current Neutral current Electromotive force Capacitance Inductance Impedance Resonant cavities Waveguides
Covariant formulation Electromagnetic tensor (stress–energy tensor) Four-current Electromagnetic four-potential
Scientists Ampère Coulomb Faraday Gauss Heaviside Henry Hertz Lorentz Maxwell Tesla Volta Weber Ørsted

Bloch wave – MoM is a first principles technique for determining the photonic band structure of triply periodic electromagnetic media such as photonic crystals. This technique uses the method of moments (MoM) in combination with a Bloch wave expansion of the electromagnetic field in the structure. It is based on the 3-dimensional spectral domain method (Kastner [1987]), specialized to periodic media. This approach is very efficient in terms of the number of plane waves needed for good convergence and is analogous to the spectral domain MoM method commonly used for analyzing 2D periodic structures such as frequency selective surfaces (FSS). In both cases, the field is expanded as a set of eigenfunction modes (either a Bloch wave in 3D or a discrete plane wave - aka Floquet mode - spectrum in 2D), and an integral equation is enforced on the surface of the scatterers in each unit cell. In the FSS case, the unit cell is 2-dimensional and in the photonic crystal case, the unit cell is 3-dimensional.

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 infinite 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 an orthogonal 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 unit cell dimensions 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 (actually, the propagation constant of whatever medium the metallic scatterers are embedded in), 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.

The Bloch wave expansions in equations (2.1) are nothing more than exponential Fourier series multiplied by the cell-to-cell propagation factor: $e^{j(\alpha _{0}x+\beta _{0}y+\gamma _{0}z)}.$ The Bloch wave expansions are chosen since any field solution within an infinite periodic volume must have the same periodicity as the medium itself, or stated another way, the fields in neighboring cells must be identical up to a (real or complex) propagation factor.

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)

is the tensor Green's function in the spectral domain. Note that the spatial domain convolution has been transformed into a simple multiplication in the spectral domain, consistent with the convolution theorem for Fourier transforms.

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 components of electric field are actually 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)~~~~~~~~~~~~~~~~~(4.1)

Substituting (4.1) 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}}~~~(4.2)

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 (see Scott [1998], available on researchgate.net). 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 (4.2) is identical to the basic EFIE for a freestanding PEC FSS (Scott [1989] or Frequency selective surface eq. (4.2)), the only difference being the stronger singularity in 3D which accelerates convergence of the triple sums (and the fact that the vectors are now 3-dimensional).

We see from (4.2) that the EFIE can become singular whenever the free space wavenumber is exactly equal to one of the wave numbers in any of the 3 periodic coordinate directions. This can happen for example when the free space wavelength exactly equals the lattice spacing. This is a statistically rare occurrence in computational practice and corresponds to a propagation anomaly similar to a Wood's anomaly for gratings.

The advantages of this method over both the plane wave expansion method and the FEM method include both mathematical and computational simplicity. Electric current unknowns are only needed on the surface of the scatterers in the unit cell, so the matrix eigenvalue problem can be as small as 1x1 for simple scatterers at low frequencies. As a result, this equation can be easily solved for relatively complex structures even on a low-end PC.

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. Basically we have to try different values of frequency with pre-selected values of propagation constant and propagation direction - through trial and error - until we find a combination that produces the desired zero field condition on the surface of the array of scatterers. Equation (4.2) has been used to compute bands in various types of doped and undoped photonic crystals (Scott[1998], Scott [2002], both available on researchgate.net).

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, available on researchgate.net
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, available on researchgate.net

External links

This article is issued from Wikipedia - version of the 11/28/2016. The text is available under the Creative Commons Attribution/Share Alike but additional terms may apply for the media files.