GTPhPermittivityMatrix

GTPhPermittivityMatrix[reciprocal basis,cutoff,objects,background]
gives the complete matrix i.e. the Fourier transform of . The reciprocal lattice vectors are defined by the reciprocal basis and vectors with < cutoff for a structure in the unit cell defined by objects. The structure is embedded in a dielectric background.

DetailsDetails

  • Two methods of the construction of the master equation, i.e. the eigenvalue problem to calculate photonic band structures are possible. First, the Fourier transform of ϵ(r) is calculated on the fly during the construction of the eigenvalue problem, second the is calculated beforehand. GTPhPermittivityMatrix allows to calculate the matrix in two ways. First, it can be calculated directly by Fourier transform of epsilon^(-1)(r), the so called direct method. In the second method, so called Hos method (method of matrix inversion), the Fouriert ransform of ϵ(r) is calculated. The corresponding matrix is inverted.
  • The following options can be given:
  • GODCMethod"Direct"Method to calculate the matrix
    GOPixelFalseDecides the type of calculation: pixelmap or list of objects
    GOVerboseFalseControls the output of additional information
  • See: W. Hergert, M. Geilhufe, Group Theory in Solid State Physics and Photonics. Problem Solving with Mathematica, chapter 10.4
  • K. M. Ho, C. T. Chan, C. M. Soukulis, Existance of a Photonic Gap in Periodic Dielectric structures. Phys. Rev. Lett. 65, 3152 (1990).

ExamplesExamplesopen allclose all

Basic Examples  (1)Basic Examples  (1)

First load the package:

In[1]:=
Click for copyable input

As an example dielectric rods in a square lattice in air are considered. The structure is defined:

In[2]:=
Click for copyable input

The Fourier transform of is calculated.

In[3]:=
Click for copyable input

The eigenvalue problem is constructed for the TM polarization.

In[4]:=
Click for copyable input

The path in the Brillouin zone is fixed.

In[5]:=
Click for copyable input

The band structure is calculated and plotted.

In[6]:=
Click for copyable input
In[7]:=
Click for copyable input
Out[7]=