Heating and phase transitions of dust-plasma crystals in a flowing plasma

A three dimensional particle simulation is used to study hexagonal dust crystals in Bowing plasmas. The flowing velocity upsilon(0) is mesothermal (upsilon(Ti)<upsilon(0)<upsilon(Te)), which is a typical situation for dust particles confined in industrial plasmas and plasma experiments. Unlike for the Yukawa system (Debye-Huckel shielded dust), we find parameter regimes where a hexagonal dust crystal is stable. This is due to the wakefield generated around the individual particles which causes attractive interparticle forces along the axis of the plasma flow. These attractive forces cause a rather strong binding of the crystal since the simulation shows that solid-fluid transitions occur at a much lower Gamma value (higher dust temperature) than for the Yukawa system. Stability of a hexagonal dust crystal is found to depend strongly upon the gas neutral pressure since stable crystal structure is obtained for dust-neutral collision frequencies gamma(d) above some threshold value gamma(d0). For gamma(d)<gamma(d0) the ion flow will excite crystal waves (phonons) where the amplitude of these waves grows as a function of time. This instability will either saturate due to nonlinear lattice waves or continue until the dust crystal structure eventually melts.