Open Access

This thesis is devoted to experiments on three-dimensional dust clouds which are confined in low temperature plasmas. Such ensembles of highly electrically charged micrometer-sized particles reveal fascinating physics, such as self-excited density waves and vortices. At the same time, these systems are challenging for experimental approaches due to their three-dimensional character. In this thesis, new optical diagnostics for dusty plasmas have been developed and, in combination with existing techniques, have been used to study these 3D dusty plasmas on different size and time scales.

During the past decade, various physical properties of the Yukawa ball, like structure and energy states, were unraveled using experiments. However, the dynamical features served further attention. Therefore, the main aim of my thesis was to investigate and understand how a finite system-represented by Yukawa clusters-evolves from a solid, crystalline structure to a liquid-like system, how it behaves in this phase and in what manner the reordering back into the solid state can be described. As a method of choice to reach this goal, laser heating has been proven successful. Moreover, the special importance of wakefields for dust clusters confined at low neutral gas pressure was addressed. Melting of finite dust clouds can be induced in two ways, either by altering the properties of the ambient plasma or by laser heating. The latter was shown to be a generic melting scenario, allowing to estimate a critical coupling parameter at the melting point. Moreover, the melting transition of finite 3D dust systems was found to be a two-step process where angular order is lost before the radial order starts to diminish at higher energies. Next, the mode dynamics of finite 3D dust ensembles in the solid and the liquid phase was studied. Crystal and fluid modes revealed the main spectral properties of the system. The normal modes are mainly suited to describe crystalline states. Fluid modes were excited naturally and via laser heating, with excitation frequencies almost independent of the coupling parameter in the solid and the liquid-like regime. Tuning the plasma parameters can be used to vary the particle-particle interaction via the ion focus. Both methods, even though assuming equilibrium situations, allowed to hint at these wakefields. The corresponding peaks in the fluid and normal mode spectra were no eigenmodes, confirming the nonequilibrium character of the ion focusing effect. First steps to extend the normal mode theory to achieve the dynamics of wake-affected nonequilibrium dust clusters were presented. Statistical quantities were obtained evaluating long-run experiments and transport coeffcients for finite dust systems were calculated via the instantaneous normal mode technique. Diffusion was found considerably higher for 3D than for 2D dust clusters. Using the configurational entropy, we have shown that in 2D and 3D disorder increases with increasing size of the system, in agreement with simulations. The temperature dependence of the configurational entropy differs for 2D and 3D dust clouds, with a threshold behavior found for finite 2D ensembles only. Finally, using instantaneous normal modes to reveal the total fraction of unstable modes, the predictive connection of Keyes (Phys Rev E 62, p7905, 2000), between transport and disorder was tested and verified for 2D, but not for 3D clusters. The reason for this has to be left open. Finally, laser-mediated recrystallization processes of finite 3D dust clouds were investigated. First, the temporal evolution of the Coulomb coupling parameter was traced during heating and recrystallization. A cooling rate has been determined from the initial phase of recrystallization. This cooling rate is lower than damping by the neutral gas, in agreement with simulations. We have observed a large fraction of metastable states for the final cluster configurations. Further, we have revealed that the time scale for the correlation buildup in the finite 3D ensemble was on even slower scales than cooling. Thus, different time scales can be attributed to the fast emergence of the shells and to the slower individual ordering within the shells.

This thesis constitutes a computational study of charge and ion drag force on micron-sized dust particles immersed in rf discharges. Knowledge of dust parameters like dust charge, floating potential, shielding and ion drag force is very crucial for explaining complex laboratory dusty plasma phenomena, such as void formation in microgravity experiments and wakefield formation in the sheaths. Existing theoretical models assume standard distribution functions for plasma species and are applicable over a limited range of flow velocities and collisionality. Kinetic simulations are suitable tools for studying dust charging and drag force computation. The main aim of this thesis is to perform three dimensional simulations using a Particle-Particle-Particle-Mesh ($P^3M$) model to understand how the dust parameters vary for different positions of dust in rf discharges and how these parameters on a dust evolve in the presence of neighboring dust particles. At first, rf discharges in argon have been modelled using a three-dimensional PIC-MCC code for the discharge conditions relevant to the dusty plasma experiments. All necessary elastic and inelastic collisions have been considered. The plasma background is found collisional, charge-exchange collisions between ions and neutrals being dominant. Electron and ion distributions are non-Maxwellian. The dominant heating mechanism is Ohmic. Then, simulations have been done to compute the dust parameters for various sizes of dust located at different positions in the rf discharges. Dust charge and floating potential in the presheath are slightly larger than the values in the bulk due to the higher electron flux to the dust particle in the presheath. From presheath to the sheath the charge and floating potential values decrease due to the decrease of the electron current to the dust. A linear dependence of dust potential on dust size has been found, which results in a nonlinear dependence of the dust charge with the dust size when the particle is assumed to be a spherical capacitor. This has been verified by independently counting the charges collected by the dust. %where indeed it has been noted that the dust charge %scales nonlinearly with the dust size. The computed dust parameters are also compared with theoretical models. Simulated dust floating potentials are comparable to values obtained from Allen-Boyd-Reynolds (ABR) and Khrapak models, but much smaller than the values obtained from Orbit Motion Limited (OML) model. The dust potential distribution behaves Debye-H\"{u}ckel-like. The shielding lengths are in between ion and electron Debye lengths. % indicating shielding by both ions and electrons. Further, the orbital drag force is typically larger than the collection drag force. The total drag force for the collisional case is larger than for the collisionless case and it scales nonlinearly with the dust size. The collection drag values and size-scaling agrees with Zobnin's model. The charging and drag force computation is then extended to two and multiple static dust particles in the rf discharge to study the influence of neighboring dust particles on the dust parameters. Initially, the dust parameters on two dust particles are computed for various interparticle separation distances and for dust particles placed at different locations in the rf discharge. It is observed that for dust separations larger than the shielding length the dust parameters for the two dust particles match with the single dust particle values. As the dust separation is equal to or less than the shielding length the ion drag force increases due to the buildup of a parallel drag force component. However, the main dust properties like charge, potential, vertical component of ion drag are not affected considerably. This is attributed to the smaller collection impact parameter values compared to the dust separation. %This is because the %collection impact parameter values in the sheath and the presheath are smaller %than the smallest dust separation and in case of the dust in the bulk, the %collection impact parameter is comparable with the dust separation. Then the dust charges on multiple dust particles located at different positions in the discharge and arranged along the discharge axis are also computed. It is found that the charges of the multiple dust particles in the bulk or presheath do not differ much from the single particle values at that location. But the dust charges of multiple dust particles located in the sheath drastically differ from the single dust parameter values. Due to ion focusing from dust particles in the upper layers, the ion current increases to dust particles in the lower layers resulting in smaller charge values. This is as well the case where dust particles are vertically aligned as in the standard experiments of dusty plasmas. In conclusion, this work used a fully kinetic (PIC and MD or $P^3M$) model to study the physics of dust charging in rf plasmas. Our simulations revealed that the dust parameters vary considerably from the bulk to the sheath. The CX collisions increase flux to the dust thereby affecting the dust parameters and their scaling with dust size. Also, a dust particle affects the charging dynamics of its neighbor only when their separation is within the shielding length. In the plasma sheath, ion focussing can cause great reduction in dust charges.