GTPhFields

GTPhFields[file name,kpoint,band number,Gvectors,geometry]
calculates the electromagnetic field in a photonic crystal.

DetailsDetails

  • The matrix representation of the photonic master equation is constructed with GTPhMasterEquation. Afterwards the eigenvalues and eigenvectors are calculated at certain k points by means of GTBands and stored.
  • The stored eigenvalues in file name are read for k point and band number. To construct the field also the reciprocal lattice vectors (G vectors) have to be provided.
  • The list geometry = {xmin,xmax,nx,ymin,ymax,ny} defines a rectangular grid. The field is calculated at nx mesh points in [xmin,xmax] in x-direction and ny mesh points in [ymin,ymax] in y-direction.
  • is calculated for TM polarization and for TE polarization.
  • Click for copyable input
    is restricted to two - dimensional photonic crystals.
  • See: W. Hergert, M. Geilhufe, Group Theory in Solid State Physics and Photonics. Problem Solving with Mathematica, chapter 10.4

ExamplesExamplesopen allclose all

Basic Examples  (1)Basic Examples  (1)

First load the package:

In[1]:=
Click for copyable input
In[2]:=
Click for copyable input

This example is discussed in Joannopoulos, Meade, Winn, Photonic Crystals, Princeton University Press 1995, p. 68.

Cylindrical pillars of permittivity estruc in a background eback with a radius ra are considered. Filling factor and radius define the list geom.

In[3]:=
Click for copyable input

The symmetry points Γ,X, and M are considered:

In[4]:=
Click for copyable input

The reciprocal lattice will be defined and reordered in shells.

In[5]:=
Click for copyable input

The matrix representation of the master equation will be set up. The calculation will done for TM polarization, if no option is given.

TM polarization:

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

The band 2 at Γ investigated

In[8]:=
Click for copyable input
Out[8]=

TE polarization:

In[9]:=
Click for copyable input
In[10]:=
Click for copyable input

The band 2 at M is investigated

In[11]:=
Click for copyable input
In[12]:=
Click for copyable input
Out[12]=

The band 3 at M is investigated

In[13]:=
Click for copyable input
Out[13]=