Dispersion Relations for Thermally Excited Waves in Plasma Crystals

Thermally excited waves in a Plasma crystal were numerically simulated using a Box_Tree code. The code is a Barnes_Hut tree code proven effective in modeling systems composed of large numbers of particles. Interaction between individual particles was assumed to conform to a Yukawa potential. Particle charge, mass, density, Debye length and output data intervals are all adjustable parameters in the code. Employing a Fourier transform on the output data, dispersion relations for both longitudinal and transverse wave modes were determined. These were compared with the dispersion relations obtained from experiment as well as a theory based on a harmonic approximation to the potential. They were found to agree over a range of 0.9<k<5, where k is the shielding parameter, defined by the ratio between interparticle distance a and dust Debye length lD. This is an improvement over experimental data as current experiments can only verify the theory up to k = 1.5.