The formation of three-dimensional (3D) dust clusters within a complex plasma modeled as a spatially confined Yukawa system is simulated using the box_tree code. Similar to unscreened Coulomb clusters, the occurrence of concentric shells with characteristic occupation numbers was observed. Both the occupation numbers and radii were found to depend on the Debye length. Ground and low energy meta-stable states of the shielded 3D Coulomb clusters were determined for 4<N<20. The structure and energy of the clusters in different states was analyzed for various Debye lengths. Structural phase transitions, including inter-shell structural phase transitions and intra-shell structural phase transitions, were observed for varying Debye length and the critical value for transitions calculated.

The vertical as well as horizontal oscillation modes in a thermally excited two-dimensional (2D) dust coulomb cluster were investigated for particle numbers between and using both a box_tree simulation and an analytical method. The horizontal mode spectra is shown to agree with published results while the vertical mode spectra obtained from the box_tree simulation and the analytical method are shown to agree with one another. The maximum frequency of the vertical modes is shown to be the vertical oscillation frequency of the whole system acting as a solid plane and the frequency is shown to decrease as the mode number l decreases. The oscillation patterns of the vertical modes are also investigated. For clusters with large numbers of particles, the high frequency modes are shown to have oscillation patterns similar in shape to Bessel-Fourier functions with various index m and n. For specified values of n, modes with higher m (except for) have lower frequencies and for specified m, modes with higher n have lower frequencies. For low frequency modes, the largest amplitude particle motion is concentrated in a few inner rings with the outer rings remaining almost motionless. This is in contrast to the horizontal modes where the strongest motion of the particles is concentrated in the inner rings at high frequency. For clusters with small numbers of particles (), these modes where the strongest vertical motion is concentrated in the inner rings does not exist resulting in oscillation magnitudes for all the modes appearing in the shape of Bessel-Fourier functions.

A self-consistent three-dimensional model for a complex (dusty) plasma is used to study the effects of multiple-sized dust grains in a dust crystal. In addition to the interparticle forces, which interact through a Yukawa potential, the model includes the effects of gravity, the variation of the sheath potential above the powered electrode, and a radial confining potential. Simulations studied various ratios of a mix of 6.5- and 8.9-micron monodisperse particles and compared their correlation functions, electric potential energy of the crystal formations, and the dispersion relations for in-plane and out-of-plane dust lattice wave (DLW) modes for two different sheath thicknesses. In the 7 mm sheath, the particles formed two layers in the vertical direction by size, and acted as a two-layer crystal with weak correlation between the layers. In the 3 mm sheath, the particles formed an essentially monolayer crystal; however the crystal dynamics showed some characteristics of a bilayer crystal.