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] adim_v2_0.gz(35 Kbytes)
Manuscript Title: MULTEM 2: a new version of the program for transmission and band- structure calculations of photonic crystals.
Authors: N. Stefanou, V. Yannopapas, A. Modinos
Program title: MULTEM, version 2
Catalogue identifier: ADIM_v2_0
Distribution format: gz
Journal reference: Comput. Phys. Commun. 132(2000)189
Programming language: Fortran.
Computer: HP 9000/778/B180L, SGI Indy IP22/150MHz, Pentium PC.
Operating system: HP-UX 10.20, IRIX6.2, Red Hat LINUX 5.2 (GNU compiler), MS Windows.
RAM: 1M words
Word size: 32
Keywords: Photonic crystals, Complex photonic band structure, Transmission and reflection coefficients, Multiple scattering of electromagnetic waves.
Classification: 10.

Nature of problem:
Calculation of the complex band structure associated with a given surface of a photonic crystal, and of the transmission, reflection and absorption coefficients of light by a slab of the crystal parallel to the given surface. We note that the ordinary frequency band structure of the infinite crystal is contained within the complex band structure of any surface of the crystal.

Solution method:
Solution of Maxwell's equations using multiple-scattering techniques.

Summary of revisions:
  1. In the expansion of the electromagnetic field into plane waves we now use a basis of transverse plane waves, thus excluding from the beginning unphysical longitudinal solutions. In this way the dimensions of the reflection/transmission matrices are reduced by a factor of 3/2. For this purpose, subroutines DLMKG, HOSLAB, and PLW have been modified, with consequent amendments in other subroutines and in the main program.
  2. In relation to a transmission calculation, the new version of the program makes it possible to have different semi-infinite homogeneous media on the left and right sides of the slab of the crystal. For this purpose appropriate changes have been made in the main program and in subroutine SCAT. For KEMB=1, in the last two lines of the input data file for a transmission calculation, we give the values of the dielectric constants and magnetic permeabilities of the media on the left and right sides of the slab.
  3. The matrix-multiplication operations in subroutine PAIR are written in a more efficient way thus leading to a considerable speed-up in the execution of the program.>
  4. A few mistakes of minor importance and inconsistencies, leading to warning or error messages in some Fortran compilers, have been corrected.

The structures that can be considered consist of parallel planes of non-overlapping spheres of given two-dimensional periodicity and uniform plates.

Running time:
For the same problem (test run in Ref. [1]) the time required on an HP 9000/778/B180L workstation by the present version is about 1.5s (2.4s) per frequency for the band-structure (transmission) calculation, instead of 20.3s (49.6s) required by the previous version. The input and results are the same and are not reproduced here.

[1] N. Stefanou, V. Yannopapas, and A. Modinos, Comput. Phys. Commun. 113 (1998) 49.