Elsevier Science Home
Computer Physics Communications Program Library
Full text online from Science Direct
Programs in Physics & Physical Chemistry
CPC Home

[Licence| Download | New Version Template] adij_v2_0.tar.gz(189 Kbytes)
Manuscript Title: A program for calculating photonic band structures, Green's functions and transmission/relfection coefficients using a non-orthogonal FDTD method.
Authors: A.J. Ward, J.B. Pendry
Program title: ONYX Version 2.1
Catalogue identifier: ADIJ_v2_0
Distribution format: tar.gz
Journal reference: Comput. Phys. Commun. 128(2000)590
Programming language: Fortran.
Computer: Digital Alpha 250.
Operating system: UNIX V3.2D-1 (Rev. 41).
RAM: 23M words
Word size: 32
Keywords: Electrostatics, Discretised Maxwell's Equations, Order (N) method, Finite difference time domain (FDTD) method, Photonic band structure, Green's function, Density of states, Transmission and reflection coefficients.
Classification: 10.

Nature of problem:
Efficient calculation of either photonic dispersion relationships, Green's functions or transmission and reflection coefficients for photons in complex dielectric structures.

Solution method:
A discretisation of Maxwell's equations in both the space and time domains which leads to finite difference equations connecting the electric and magnetic fields at one time step to those at the next. After using these equations to find the response in the time domain to a particular initial set of fields, we perform a Fourier transform to obtain the response in the frequency domain. From this we can easily extract dispersion relationship information, or alternatively, by setting the initial fields to be a delta function, we can obtain the Green's function for the system under consideration. In addition, by projecting onto a complete basis set of plane waves we can find the transmission and reflection matrices for the scattering system.

The complexity of the dielectric structure that the method can be applied to is limited only by the computer time and memory available. Both time and memory requirements scale linearly with the system size. One restriction on the method is that the dielectric permittivity and magnetic permeability must both be independent of frequency. This means it cannot treat some problems, typically those involving metals.

Unusual features:
Option to work with non-orthogonal co-ordinate systems.

Running time:
Highly dependent on the system under consideration. For the test transmission calculation given, 120 seconds on a Digital Alpha 250 workstation.