Abstract-for the j th particle includes electrostatic forces due to particleparticle interaction φ jl , a confinement potential , gas friction −ν g m pṙ j , and random Brownian force ζ g j due to gas molecules. Here, m p is the particle mass and ν g is the gas friction constant. We model the interparticle potential as a pairwise Yukawa potential, φ jl (r jl ) = (Q 2 /4πε 0 r jl )e −r jl /λ D for identical dust particles of charge Q with a screening length λ D due to the electrons and ions. The collection of dust particles can be characterized by two dimensionless parameters: the Coulomb coupling parameter = Q 2 /4πε 0 ak B T d and the screening parameter κ = a/λ D . Here a = (3/4πn d ) 1/3 is the Wigner-Seitz radius, n d is the number density of dust particles, and T d is the particle kinetic temperature. For the coupling parameter > 1, the dust component is said to be strongly coupled, and dust particles can self-organize like atoms in a solid or liquid and sustain waves [2], [3].Our simulation parameters are for the PK-4 instrument [4]. We use microsphere dust particles of radius 3.43 μm and m p = 2.55 × 10 −13 kg, with neon gas at 50 Pa pressure and 0.03 eV temperature so that ν g = 51 s −1 . We assume Q = −8520e, n d = 3 × 10 4 cm −3 , and λ D = 8.3 × 10 −3 cm. The characteristic interparticle distance is a = 0.020 cm, so that κ = 2.4. The characteristic time for particle motion is ω p = 157 rad/s, whereWe chose T d = 8.3 eV, corresponding to ≈ 63. For these values of and κ, the collection of dust particles is predicted Results shown in Fig. 1(a) reveal the structural arrangement of the dust particles at a time during the simulation. This image was prepared by plotting the simulated particles in a 3-D coordinate system, with a sphere representing each particle. In this structure, each particle is in a cage defined by its nearest neighbors, but the structure is irregular, not crystalline. A video can be seen at [6] showing similar particles as shown in Fig. 1(a) from a rotating viewpoint.To quantify the order of the 3-D structure, we calculate the pair correlation function g(r ) [7], [8]. For this liquid, g(r ) has only one distinctive peak in Fig. 1(b), indicating short-range translational order.We characterize the dynamics using a wave spectrum. We start by using the particle position r j (t) and velocitẏ r j (t) to calculate the time series of the so-called longitudinal current, for a specified wave vector kThe spectral power | J L (k, ω)| 2 is then computed as the square modulus of the Fourier transformation in time of J L (k, t).Results in Fig. 1(c) show that, as expected, the spectral power is concentrated along a curved band. The band has a great width in the ω-k space due to damping, arising from gas friction and the viscous motion of dust particles.The band of spectral power in Fig. 1(c) corresponds to a real dispersion relation curve, which we plot in Fig. 1(c) as a dotted line. We determined this dispersion curve as the peak of the spectral power |J L (k, ω)| 2 ; to reduce the uncertainty, we compu...