A layer-multiple-scattering method for phononic crystals and heterostructures of such

Abstract

We present a computer program to calculate the frequency band structure of an infinite phononic crystal, and the transmission, reflection and absorption of elastic waves by a slab of this crystal. The crystal consists of a stack of identical slices parallel to a given surface; the slice may consist of multilayers of non-overlapping spheres of given periodicity parallel to the surface and homogeneous plates. The elastic coefficients of the various components of the crystal may be complex functions of the frequency. Program summary Title of program: MULTEL Catalogue identifier:ADUT Program summary URL: http://cpc.cs.qub.ac.uk/summaries/ADUT Program obtainable from: CPC Program Library, Queen's University of Belfast, N. Ireland Computers: Pentium 4 PC, HP J-5600 Installation: University of Athens, Section of Solid State Physics Operating systems: MS Windows, LINUX, HP-UX Programming language used: FORTRAN 77 Memory required to execute with typical data: 0.9M words Number of bits in a word: 32 Number of processors used: 1 Has the code been vectorized or parallelized?: No Distribution format:tar.gz. Number of lines in distributed program, including test data, etc.: 8676 Number of bytes in distributed program, including test data, etc.: 61 459 Nature of physical problem: Calculation of the complex band structure associated with a given surface of a phononic crystal, and of the transmission, reflection and absorption coefficients of elastic waves 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. Method of solution: Solution of the equations of elasticity using multiple-scattering techniques. Restrictions on the complexity of the program: The structures that can be considered consist of parallel planes of non-overlapping spheres of given two-dimensional periodicity and uniform plates. Typical running time: For the given test run, about 0.1 s per frequency for the band structure calculation and 0.1 s per frequency for the transmission/reflection/absorption calculation, on a Pentium 4 (3.00 GHz) PC.