Plane Wave Expansion Code for: Acoustic Logic Gates Implemented Using a Phase-Controlling Phononic Crystal

The plane wave expansion (PWE) method  code used in the "Acoustic Logic Gates Implemented Using a Phase-Controlling Phonic Crystal" Masters Thesis  granted by the University of Arizona. The code is written Matlab (worked up to R2011b) and is based on a FORTRAN code originally maintained by N. Swinteck. The program is designed to calculate the band structure of a Phononic/Sonic Crystal using the PWE method commonly used in electromagnetics. The results are displayed in terms of reduced wave number coordinates. 

The user should define the parameters at the beginning of the script PWE_Zmodes.m which includes:

%% Lattice Properties
lat_par: lattice parameter spacing [L:meters]
radius: radius of cylinder inclusions [L:meters]
fillfract: filling fraction
latticekey: Type of lattice, 1-> Square, 2-> Hexagonal
size: Number of repeating units

%% Material Properties 
rho_matr: Density of matrix material [M:kilograms]
speed_matr: speed of sound in matrix material [L/T : meters/seconds]
rho_inc: density of lattice inclusion material [M: kilograms]
speed_inc: speed of sound in inclusion material [L/T: meters/seconds]

%% Number of bands to calculate
number_passbands: number of bandstructure eigenvalues to calculate

The sampling of in reciprocal space is fixed to the irreducible symmetry directions for both a square and hexagonal lattice.The reciprocal lattice is calculated using the subroutine in reciplat.m. The code can be executed with the GNU Octave program however plotting options have to be changed to make it work with Octave's graphing interface.